The programming language python in earth system simulations

-scale vary from a planetary scale and million years for convection problems to 100km and 10 years for fault systems simulations. Various techniques are in use to deal with the time dependency (e.g. Crank-Nicholson), with the non-linearity (e.g. Newton-Raphson) and weakly coupled equations (e.g. non-linear Gauss-Seidel). Besides these high-level solution algorithms discretization methods (e.g. finite element method (FEM), boundary element method (BEM)) are used to deal with spatial derivatives. Typically, large-scale, three dimensional meshes are required to resolve geometrical complexity (e.g. in the case of fault systems) or features in the solution (e.g. in mantel convection simulations). The modelling environment escript allows the rapid implementation of new physics as required for the development of simulation codes in earth sciences. Its main object is to provide a programming language, where the user can define new models and rapidly develop high-level solution algorithms. The current implementation is linked with the finite element package finley as a PDE solver. However, the design is open and other discretization technologies such as finite differences and boundary element methods could be included. escript is implemented as an extension of the interactive programming environment python (see www.python.org). Key concepts introduced are Data objects, which are holding values on nodes or elements of the finite element mesh, and linearPDE objects, which are defining linear partial differential equations to be solved by the underlying discretization technology. In this paper we will show the basic concepts of escript and will show how escript is used to implement a simulation code for interacting fault systems. We will show some results of large-scale, parallel simulations on an SGI Altix system. Acknowledgements: Project work is supported by Australian Commonwealth Government through the Australian Computational Earth Systems Simulator Major National Research Facility, Queensland State Government Smart State Research Facility Fund, The University of Queensland and SGI.

Pilot-scale extraction of PHB from recombinant E-coli by homogenization and centrifugation

A new method of poly-beta-hydroxybutyrate (PHB) extraction from recombinant E. coli is proposed, using homogenization and centrifugation coupled with sodium hypochlorite treatment. The size of PHB granules and cell debris in homogenates was characterised as a function of the number of homogenization passes. Simulation was used to develop the PHB and cell debris fractionation system, enabling numerical examination of the effects of repeated homogenization and centrifuge-feedrate variation. The simulation provided a good prediction of experimental performance. Sodium hypochlorite treatment was necessary to optimise PHB fractionation. A PHB recovery of 80% at a purity of 96.5% was obtained with the final optimised process. Protein and DNA contained in the resultant product were negligible. The developed process holds promise for significantly reducing the recovery cost associated with PHB manufacture.

Prediction of MHC class II-binding peptides using an evolutionary algorithm and artificial neural network

Motivation: Prediction methods for identifying binding peptides could minimize the number of peptides required to be synthesized and assayed, and thereby facilitate the identification of potential T-cell epitopes. We developed a bioinformatic method for the prediction of peptide binding to MHC class II molecules. Results: Experimental binding data and expert knowledge of anchor positions and binding motifs were combined with an evolutionary algorithm (EA) and an artificial neural network (ANN): binding data extraction --> peptide alignment --> ANN training and classification. This method, termed PERUN, was implemented for the prediction of peptides that bind to HLA-DR4(B1*0401). The respective positive predictive values of PERUN predictions of high-, moderate-, low- and zero-affinity binder-a were assessed as 0.8, 0.7, 0.5 and 0.8 by cross-validation, and 1.0, 0.8, 0.3 and 0.7 by experimental binding. This illustrates the synergy between experimentation and computer modeling, and its application to the identification of potential immunotheraaeutic peptides.

Determination of pindolol enantiomers in human plasma and urine by simple liquid-liquid extraction and high-performance liquid chromatography

A simple method for the measurement of pindolol enantiomers by HPLC is presented. Alkalinized serum or urine is extracted with ethyl acetate and the residue remaining after evaporation of the organic layer is then derivatised with (S)-(-)-alpha-methylbenzyl isocyanate. The diastereoisomers of derivatised pindolol and metoprolol (internal standard) are separated by high-performance liquid chromatography (HPLC) using a C-18 silica column and detected using fluorescence (excitation lambda: 215 nm, emission lambda: 320 nm). The assay displays reproducible linearity for pindolol enantiomers with a correlation coefficient of r(2) greater than or equal to 0.998 over the concentration range 8-100 ng ml(-1) for plasma and 0.1-2.5 mu g ml(-1) for urine. The coefficient of variation for accuracy and precision of the quality control samples for both plasma and urine are consistently

PCR bias toward the wild-type k-ras and p53 sequences: Implications for PCR detection of mutations and cancer diagnosis

PCR-based cancer diagnosis requires detection of rare mutations in k-ras, p53 or other genes. The assumption has been that mutant and wild-type sequences amplify with near equal efficiency, so that they are eventually present in proportions representative of the starting material. Work factor IX suggests that this assumption is invalid for one case of near-sequence identity To test the generality of this phenomenon and its relevance to cancer diagnosis, primers distant from point mutations in p53 and k-ras were used to amplify, wild-type and mutant sequences from these genes. A substantial bias against PCR amplification of mutants was observed for two regions of the p53 gene and one region of k-ras. For kras and p53, bias was observed when the wild-type and mutant sequences were amplified separately or when mixed in equal proportions before PCR. Bias was present with proofreading and non-proofreading polymerases. Mutant and wild-type segments of the factor V cystic fibrosis transmembrane conductance regulator and prothrombin genes were amplified and did not exhibit PCR bias. Therefore, the assumption of equal PCR efficiency for point mutant and wild-type sequences is invalid in several systems. Quantitative or diagnostic PCR will require validation for each locus, and enrichment strategies may be needed to optimize detection of mutants.

A protocol for the fluorometric quantification of mGFP5-ER and sGFP(S65T) in transgenic plants

The Green Fluorescent Protein (GFP) from Aequorea victor-in has begun to be used as a reporter protein in plants. It is particularly useful as GFP fluorescence can be detected in a non-destructive manner, whereas detection of enzyme-based reporters often requires destruction of the plant tissue. The use of GFP as a reporter enables transgenic plant tissues to be screened in vivo at any growth stage. Quantification of GFP in transgenic plant extracts will increase the utility of GFP as a reporter protein. We report herein the quantification of a mGFP5-ER Variant in tobacco leaf extracts by UV excitation and a sGFP(S65T) variant in sugarcane leaf and callus extracts by blue light excitation using the BioRad VersaFluor(TM) Fluorometer System or the Labsystems Fluoroskan Ascent FL equipped with a narrow band emission filter (510 +/- 5 nm). The GFP concentration in transgenic plant extracts was determined from a GFP-standard series prepared in untransformed plant extract with concentrations ranging from 0.1 to 4 mu g/ml of purified rGFP. Levels of sgfp(S65T) expression, driven by the maize ubiquitin promoter, in sugarcane calli and leaves ranged up to 0.525 mu g and 2.11 mu g sGFP(S65T) per mg of extractable protein respectively. In tobacco leaves the expression of mgfPS-ER, driven by the cauliflower mosaic virus (CaMV) 35S promoter, ranged up to 7.05 mu g mGFP5-ER per mg extractable protein.

Quantification of free mycophenolic acid by high-performance liquid chromatography-atmospheric pressure chemical ionisation tandem mass spectrometry

To facilitate the investigation of free mycophenolic acid concentrations we developed a high-performance liquid chromatography tandem mass spectrometry method using indomethacin as an internal standard. Free drug was isolated from plasma samples (500 mul) using ultrafiltration, The analytes were extracted from the ultrafiltrate (200 mul) using C-18 solid-phase extraction. Detection was by selected reactant monitoring of mycophenolic acid (m/z 318.9-->190.9) and the internal standard (m/z 356.0-->297.1) with an atmospheric pressure chemical ionisation interface. The total chromatographic analysis time was 12 min. The method was found to be linear over the range investigated, 2.5-200 mug/l (r>0.990, n=6). The relative recovery of the method for the control samples studied (7.5, 40.0 and 150 mug/l) ranged from 95 to 104%. The imprecision of the method, expressed in terms of intra- and inter-day coefficients of variation, was

Extensive chromosomal variation associated with taxon divergence and host specificity in the gall-inducing scale insect Apiomorpha munita (Schrader) (Hemiptera : Sternorrhyncha : Coccoidea : Eriococcidae)

Apiomorpha Rubsaamen (Hemiptera: Coccoidea: Eriococcidae) is one of the most chromosomally diverse of all animal genera. There is extensive karyotypic variation within many of the morphologically defined species, including A. munita (Schrader) which is here reported to have diploid chromosome counts ranging from 6 to more than 100. Each of the three morphologically defined subspecies of A. munita also displays considerable chromosomal variation: A. m. tereticornuta Gullan (2n =6, 8, 20, 22 or 24), A. m. malleensis Gullan (2n =6, 20, 22, 24 or 26), and A. m. munita (Schrader) (2n=54 or >100). Apiomorpha munita appears to occur only on eucalypts of the informal subgenus Symphyomyrtus, with each of the subspecies of A. munita restricted to discrete symphyomyrt sections. Several different karyotypic forms within each subspecies of A. munita appear to be restricted to only one or a few eucalypt species or series. The association between apparent host specificity and chromosomal rearrangements in A. munita suggests that both may be playing an active role in taxon divergence in Apiomorpha. (C) 2001 The Linnean Society of London.

Expression, purification, crystallization, and NMR studies of the helicase interaction domain of Escherichia coli DnaG primase

in Escherichia coli, the DnaG primase is the RNA polymerase that synthesizes RNA primers at replication forks. It is composed of three domains, a small N-terminal zinc-binding domain, a larger central domain responsible for RNA synthesis, and a C-terminal domain comprising residues 434-581 [DnaG(434-581)] that interact with the hexameric DnaB helicase. Presumably because of this interaction, it had not been possible previously to express the C-terminal domain in a stably transformed E coli strain. This problem was overcome by expression of DnaG(434-581) under control of tandem bacteriophage gimel-promoters, and the protein was purified in yields of 4-6 mg/L of culture and studied by NMR. A TOCSY spectrum of a 2 mM solution of the protein at pH 7.0, indicated that its structured core comprises residues 444-579. This was consistent with sequence conservation among most-closely related primases. Linewidths in a NOESY spectrum of a 0.5 mM sample in 10 mM phosphate, pH 6.05, 0.1 M NaCl, recorded at 36 degreesC, indicated the protein to be monomeric. Crystals of selenomethionine-substituted DnaG(434-581) obtained by the hanging-drop vapor-diffusion method were body-centered tetragonal, space group I4(1)22, with unit cell parameters a = b 142.2 Angstrom, c = 192.1 Angstrom, and diffracted beyond 2.7 Angstrom resolution with synchrotron radiation. (C) 2003 Elsevier Inc. All rights reserved.

Computational methods for prediction of T-cell epitopes - a framework for modelling, testing, and applications

Computational models complement laboratory experimentation for efficient identification of MHC-binding peptides and T-cell epitopes. Methods for prediction of MHC-binding peptides include binding motifs, quantitative matrices, artificial neural networks, hidden Markov models, and molecular modelling. Models derived by these methods have been successfully used for prediction of T-cell epitopes in cancer, autoimmunity, infectious disease, and allergy. For maximum benefit, the use of computer models must be treated as experiments analogous to standard laboratory procedures and performed according to strict standards. This requires careful selection of data for model building, and adequate testing and validation. A range of web-based databases and MHC-binding prediction programs are available. Although some available prediction programs for particular MHC alleles have reasonable accuracy, there is no guarantee that all models produce good quality predictions. In this article, we present and discuss a framework for modelling, testing, and applications of computational methods used in predictions of T-cell epitopes. (C) 2004 Elsevier Inc. All rights reserved.

Purification by reflux electrophoresis of whey proteins and of a recombinant protein expressed in Dictyostelium discoideum

Protein purification that combines the use of molecular mass exclusion membranes with electrophoresis is particularly powerful as it uses properties inherent to both techniques. The use of membranes allows efficient processing and is easily scaled up, while electrophoresis permits high resolution separation under mild conditions. The Gradiflow apparatus combines these two technologies as it uses polyacrylamide membranes to influence electrokinetic separations. The reflux electrophoresis process consists of a series of cycles incorporating a forward phase and a reverse phase. The forward phase involves collection of a target protein that passes through a separation membrane before trailing proteins in the same solution. The forward phase is repeated following clearance of the membrane in the reverse phase by reversing the current. We have devised a strategy to establish optimal reflux separation parameters, where membranes are chosen for a particular operating range and protein transfer is monitored at different pH values. In addition, forward and reverse phase times are determined during this process. Two examples of the reflux method are described. In the first case, we describe the purification strategy for proteins from a complex mixture which contains proteins of higher electrophoretic mobility than the target protein. This is a two-step procedure, where first proteins of higher mobility than the target protein are removed from the solution by a series of reflux cycles, so that the target protein remains as the leading fraction. In the second step the target protein is collected, as it has become the leading fraction of the remaining proteins. In the second example we report the development of a reflux strategy which allowed a rapid one-step preparative purification of a recombinant protein, expressed in Dictyostelium discoideum. These strategies demonstrate that the Gradiflow is amenable to a wide range of applications, as the protein of interest is not necessarily required to be the leading fraction in solution. (C) 1997 Elsevier Science B.V.

Ethical considerations for drug abuse epidemiological research

The mortality and morbidity caused by alcohol, tobacco and illicit drug misuse represents a significant public health burden (Ezzati et al., 2002). A key part of the public health response is the collection of epidemiological and social science data to define at-risk populations to identify opportunities for intervention and to evaluate the effectiveness of policies in preventing or treating drug misuse and drug-related harm. The systematic use of epidemiological and social science research methods to study illicit drug use is barely 40 years old in the United States and United Kingdom, which have pioneered this approach. Because of the sensitive nature of epidemiological research on illicit drug use a unique set of ethical challenges need to be explicitly addressed by the field. Although ethics guidelines have been proposed (Council for International Organizations of Medical Sciences, 1991), scholarship on the ethics of epidemiology is scant, and consensus on core values not yet achieved (Coughlin, 2000).

Implicit Taylor methods for stiff stochastic differential equations

In this paper we discuss implicit Taylor methods for stiff Ito stochastic differential equations. Based on the relationship between Ito stochastic integrals and backward stochastic integrals, we introduce three implicit Taylor methods: the implicit Euler-Taylor method with strong order 0.5, the implicit Milstein-Taylor method with strong order 1.0 and the implicit Taylor method with strong order 1.5. The mean-square stability properties of the implicit Euler-Taylor and Milstein-Taylor methods are much better than those of the corresponding semi-implicit Euler and Milstein methods and these two implicit methods can be used to solve stochastic differential equations which are stiff in both the deterministic and the stochastic components. Numerical results are reported to show the convergence properties and the stability properties of these three implicit Taylor methods. The stability analysis and numerical results show that the implicit Euler-Taylor and Milstein-Taylor methods are very promising methods for stiff stochastic differential equations.

Development of ligase-assisted spacer addition for the measurement of microsatellites

Conventional methods for detecting differences in microsatellite repeat lengths rely on electrophoretic fractionation on long denaturing polyacrylamide gels, a time-consuming and labor-intensive method. Therefore, there is a need for the development of new and rapid approaches to routinely detect such length polymorphisms. The advent of techniques allowing the coupling of DNA molecules to solid surfaces has provided new prospects in the area of mutation. We describe here the development and optimization of the ligase-assisted spacer addition (LASA) method, a novel and rapid procedure based on an ELISA format to measure microsatellite repeat lengths. The LASA assay was successfully applied to a set of 11 bird samples to assess its capability as a genotyping method.