Computer simulations of coupled problems in geological and geochemical systems

The finite element method is used to simulate coupled problems, which describe the related physical and chemical processes of ore body formation and mineralization, in geological and geochemical systems. The main purpose of this paper is to illustrate some simulation results for different types of modelling problems in pore-fluid saturated rock masses. The aims of the simulation results presented in this paper are: (1) getting a better understanding of the processes and mechanisms of ore body formation and mineralization in the upper crust of the Earth; (2) demonstrating the usefulness and applicability of the finite element method in dealing with a wide range of coupled problems in geological and geochemical systems; (3) qualitatively establishing a set of showcase problems, against which any numerical method and computer package can be reasonably validated. (C) 2002 Published by Elsevier Science B.V.

Stochastic gauges in quantum dynamics for many-body simulations

Quantum dynamics simulations can be improved using novel quasiprobability distributions based on non-orthogonal Hermitian kernel operators. This introduces arbitrary functions (gauges) into the stochastic equations. which can be used to tailor them for improved calculations. A possible application to full quantum dynamic simulations of BEC's is presented. (C) 2001 Elsevier Science B.V. All rights reserved.

The surface accessibility of the glycine receptor M2-M3 loop is increased in the channel open state

Mutations in the extracellular M2-M3 loop of the glycine receptor (GlyR) alpha1 subunit have been shown previously to affect channel gating. In this study, the substituted cysteine accessibility method was used to investigate whether a structural rearrangement of the M2-M3 loop accompanies GlyR activation. All residues from R271C to V277C were covalently modified by both positively charged methanethiosulfonate ethyltrimethylammonium (MTSET) and negatively charged methanethiosulfonate ethylsulfonate (MTSES), implying that these residues form an irregular surface loop. The MTSET modification rate of all residues from R271C to K276C was faster in the glycine-bound state than in the unliganded state. MTSES modification of A272C, L274C, and V277C was also faster in the glycine-bound state. These results demonstrate that the surface accessibility of the M2-M3 loop is increased as the channel transitions from the closed to the open state, implying that either the loop itself or an overlying domain moves during channel activation.

Simulations of thermal Bose fields in the classical limit

We demonstrate that the time-dependent projected Gross-Pitaevskii equation (GPE) derived earlier [M. J. Davis, R. J. Ballagh, and K. Burnett, J. Phys. B 34, 4487 (2001)] can represent the highly occupied modes of a homogeneous, partially-condensed Bose gas. Contrary to the often held belief that the GPE is valid only at zero temperature, we find that this equation will evolve randomized initial wave functions to a state describing thermal equilibrium. In the case of small interaction strengths or low temperatures, our numerical results can be compared to the predictions of Bogoliubov theory and its perturbative extensions. This demonstrates the validity of the GPE in these limits and allows us to assign a temperature to the simulations unambiguously. However, the GPE method is nonperturbative, and we believe it can be used to describe the thermal properties of a Bose gas even when Bogoliubov theory fails. We suggest a different technique to measure the temperature of our simulations in these circumstances. Using this approach we determine the dependence of the condensate fraction and specific heat on temperature for several interaction strengths, and observe the appearance of vortex networks. Interesting behavior near the critical point is observed and discussed.

Database technologies for l-system simulations in virtual plant applications on bioinformatics

One of the most important advantages of database systems is that the underlying mathematics is rich enough to specify very complex operations with a small number of statements in the database language. This research covers an aspect of biological informatics that is the marriage of information technology and biology, involving the study of real-world phenomena using virtual plants derived from L-systems simulation. L-systems were introduced by Aristid Lindenmayer as a mathematical model of multicellular organisms. Not much consideration has been given to the problem of persistent storage for these simulations. Current procedures for querying data generated by L-systems for scientific experiments, simulations and measurements are also inadequate. To address these problems the research in this paper presents a generic process for data-modeling tools (L-DBM) between L-systems and database systems. This paper shows how L-system productions can be generically and automatically represented in database schemas and how a database can be populated from the L-system strings. This paper further describes the idea of pre-computing recursive structures in the data into derived attributes using compiler generation. A method to allow a correspondence between biologists' terms and compiler-generated terms in a biologist computing environment is supplied. Once the L-DBM gets any specific L-systems productions and its declarations, it can generate the specific schema for both simple correspondence terminology and also complex recursive structure data attributes and relationships.

Nerve growth factor overexpression and autocrine loop in breast cancer cells

We show here that nerve growth factor (NGF), the canonical neurotrophic factor, is synthesized and released by breast cancer cells. High levels of NGF transcript and protein were detected in breast cancer cells by reverse transcription-PCR, Western blotting, ELISA assay and immunohistochemistry. Conversely, NGF production could not be detected in normal breast epithelial cells at either the transcriptional or protein level. Confocal analysis indicated the presence of NGF within classical secretion vesicles. Breast cancer cell-produced NGF was biologically active, as demonstrated by its ability to induce the neuronal differentiation of embryonic neural precursor cells. Importantly, the constitutive growth of breast cancer cells was strongly inhibited by either NGF-neutralizing antibodies or K-252a, a pharmacological inhibitor of NGF receptor TrkA, indicating the existence of an NGF autocrine loop. Together, our data demonstrate the physiological relevance of NGF in breast cancer and its potential interest as a marker and therapeutic target.

Mutational analysis of the large periplasmic loop 7-8 of the putrescine transporter PotE in Escherichia coli

The PotE protein is a putrescine-ornithine antiporter found in many gram-negative bacteria. It is a member of the APA family of transporters and has 12 predicted alpha-helical transmembrane spanning segments (TMS). While the substrate binding site has previously been mapped to a region near the surface of the cytoplasmic lipid layer, no structural feature within the periplasmic domains of PotE have been shown to be important for function. We examined the role of the only large outer loop, situated between transmembrane spanning segment 7 and 8, in putrescine uptake. Deletion of the highly conserved amino acids in the region closest to transmembrane spanning segment 7 produced a protein with little activity. Glycine-scanning mutagenesis of this region showed that Val(249) and Leu(254) were required for optimal transporter function. The V249G mutant transported putrescine at a lower maximal rate compared to wild-type (WT) but with the same substrate binding affinity. In contrast, the L254G mutant had a higher substrate affinity. A series of Val(249) mutants indicated that the hydrophobicity of this residue, which is located at or near the membrane surface, is important for PotE function. Secondary structure predictions of the large outer loop indicated the presence of a hydrophobic alpha-helix in the centre with a hydrophobic region at each end suggesting that the loop was not entirely exposed to the aqueous periplasmic space. The study shows that loop 7-8 is important for PotE function, possibly by forming a re-entrant loop in the channel of the transporter. (C) 2003 Elsevier Ltd. All rights reserved.

Ion trap simulations of quantum fields in an expanding universe

We propose an experiment in which the phonon excitation of ion(s) in a trap, with a trap frequency exponentially modulated at rate kappa, exhibits a thermal spectrum with an Unruh temperature given by k(B)T=h kappa. We discuss the similarities of this experiment to the response of detectors in a de Sitter universe and the usual Unruh effect for uniformly accelerated detectors. We demonstrate a new Unruh effect for detectors that respond to antinormally ordered moments using the ion's first blue sideband transition.

Roles of transmembrane domain 2 and the first intracellular loop in human noradrenaline transporter function: pharmacological and SCAM analysis

The aim was to investigate the roles of transmembrane domain 2 and the adjacent region of the first intracellular loop in determining human noradrenaline transporter (hNET) function by pharmacological and substituted-cysteine accessibility method (SCAM) analyses. It was first necessary to establish a suitable background NET for SCAM. Alanine mutants of endogenous hNET cysteines, hC86A, hC131A and hC339A, were examined and showed no marked effects on expression or function. hNET and the mutants were also resistant to methanethiosulfonate (MTS), ethylammonium (MTSEA) and MTStrimethylammonium (MTSET). Hence, wild-type hNET is an appropriate background for production of cysteine mutants for SCAM. Pharmacological investigation showed that all mutants except hT99C and hL109C showed reduced cell-surface expression, while all except hM107C showed a reduction in functional activity. The mutations did not markedly affect the apparent affinities of substrates, but apparent affinities of cocaine were decreased 7-fold for hP97C and 10-fold for hF101C and increased 12-fold for hY98C. [H-3]Nisoxetine binding affinities were decreased 13-fold for hP97C and 5-fold for hF101C. SCAM analysis revealed that only hL102C was sensitive to 1.25 mM MTSEA, and this sensitivity was protected by noradrenaline, nisoxetine and cocaine. The results suggest that this region of hNET is important for interactions with antidepressants and cocaine, but it is probably not involved in substrate translocation mechanisms.

A picrotoxin-specific conformational change in the glycine receptor M2-M3 loop.

The external loop linking the M2 and M3 transmembrane domains is crucial for coupling agonist binding to channel gating in the glycine receptor chloride channel (GlyR). A substituted cysteine accessibility scan previously showed that glycine activation increased the surface accessibility of 6 contiguous residues (Arg(271) Lys(276)) toward the N-terminal end of the homomeric alpha 1 GlyR M2 - M3 loop. In the present study we used a similar approach to determine whether the allosteric antagonist, picrotoxin, could impose conformational changes to this domain that cannot be induced by varying agonist concentrations alone. Picrotoxin slowed the reaction rate of a sulfhydryl-containing compound ( MTSET) with A272C, S273C, and L274C. Before interpreting this as a picrotoxin-specific conformational change, it was necessary to eliminate the possibility of steric competition between picrotoxin and MTSET. Accordingly, we showed that picrotoxin and the structurally unrelated blocker, bilobalide, were both trapped in the R271C GlyR in the closed state and that a point mutation to the pore-lining Thr(6') residue abolished inhibition by both compounds. We also demonstrated that the picrotoxin dissociation rate was linearly related to the channel open probability. These observations constitute a strong case for picrotoxin binding in the pore. We thus conclude that the picrotoxin-specific effects on the M2 - M3 loop are mediated allosterically. This suggests that the M2 - M3 loop responds differently to the occupation of different binding sites.

The adsorption of water in finite carbon pores

Grand canonical Monte Carlo simulations were applied to the adsorption of SPCE model water in finite graphitic pores with different configurations of carbonyl functional groups on only one surface and several pore sizes. It was found that almost all finite pores studied exhibit capillary condensation behaviour preceded by adsorption around the functional groups. Desorption showed the reverse transitions from a filled to a near empty pore resulting in a clear hysteresis loop in all pores except for some of the configurations of the 1.0nm pore. Carbonyl configurations had a strong effect on the filling pressure of all pores except, in some cases, in 1.0nm pores. A decrease in carbonyl neighbour density would result in a higher filling pressure. The emptying pressure was negligibly affected by the configuration of functional groups. Both the filling and emptying pressures increased with increasing pore size but the effect on the emptying pressure was much less. At pressures lower than the pore filling pressure, the adsorption of water was shown to have an extremely strong dependence on the neighbour density with adsorption changing from Type IV to Type III to linear as the neighbour density decreased. The isosteric heat was also calculated for these configurations to reveal its strong dependence on the neighbour density. These results were compared with literature experimental results for water and carbon black and found to qualitatively agree.

First-principles quantum simulations of dissociation of molecular condensates: Atom correlations in momentum space

We investigate the quantum many-body dynamics of dissociation of a Bose-Einstein condensate of molecular dimers into pairs of constituent bosonic atoms and analyze the resulting atom-atom correlations. The quantum fields of both the molecules and atoms are simulated from first principles in three dimensions using the positive-P representation method. This allows us to provide an exact treatment of the molecular field depletion and s-wave scattering interactions between the particles, as well as to extend the analysis to nonuniform systems. In the simplest uniform case, we find that the major source of atom-atom decorrelation is atom-atom recombination which produces molecules outside the initially occupied condensate mode. The unwanted molecules are formed from dissociated atom pairs with nonopposite momenta. The net effect of this process-which becomes increasingly significant for dissociation durations corresponding to more than about 40% conversion-is to reduce the atom-atom correlations. In addition, for nonuniform systems we find that mode mixing due to inhomogeneity can result in further degradation of the correlation signal. We characterize the correlation strength via the degree of squeezing of particle number-difference fluctuations in a certain momentum-space volume and show that the correlation strength can be increased if the signals are binned into larger counting volumes.

Mimicking the action of GroEL in molecular dynamics simulations: Application to the refinement of protein structures

Bacterial chaperonin, GroEL, together with its co-chaperonin, GroES, facilitates the folding of a variety of polypeptides. Experiments suggest that GroEL stimulates protein folding by multiple cycles of binding and release. Misfolded proteins first bind to an exposed hydrophobic surface on GroEL. GroES then encapsulates the substrate and triggers its release into the central cavity of the GroEL/ES complex for folding. In this work, we investigate the possibility to facilitate protein folding in molecular dynamics simulations by mimicking the effects of GroEL/ES namely, repeated binding and release, together with spatial confinement. During the binding stage, the (metastable) partially folded proteins are allowed to attach spontaneously to a hydrophobic surface within the simulation box. This destabilizes the structures, which are then transferred into a spatially confined cavity for folding. The approach has been tested by attempting to refine protein structural models generated using the ROSETTA procedure for ab initio structure prediction. Dramatic improvements in regard to the deviation of protein models from the corresponding experimental structures were observed. The results suggest that the primary effects of the GroEL/ES system can be mimicked in a simple coarse-grained manner and be used to facilitate protein folding in molecular dynamics simulations. Furthermore, the results Sur port the assumption that the spatial confinement in GroEL/ES assists the folding of encapsulated proteins.

Molecular dynamics simulations of the hydrophobin SC3 at a hydrophobic/hydrophilic interface

Hydrophobins are small (similar to 100 aa) proteins that have an important role in the growth and development of mycelial fungi. They are surface active and, after secretion by the fungi, self-assemble into amphipathic membranes at hydrophobic/hydrophilic interfaces, reversing the hydrophobicity of the surface. In this study, molecular dynamics simulation techniques have been used to model the process by which a specific class I hydrophobin, SC3, binds to a range of hydrophobic/ hydrophilic interfaces. The structure of SC3 used in this investigation was modeled based on the crystal structure of the class II hydrophobin HFBII using the assumption that the disulfide pairings of the eight conserved cysteine residues are maintained. The proposed model for SC3 in aqueous solution is compact and globular containing primarily P-strand and coil structures. The behavior of this model of SC3 was investigated at an air/water, an oil/water, and a hydrophobic solid/water interface. It was found that SC3 preferentially binds to the interfaces via the loop region between the third and fourth cysteine residues and that binding is associated with an increase in a-helix formation in qualitative agreement with experiment. Based on a combination of the available experiment data and the current simulation studies, we propose a possible model for SC3 self-assembly on a hydrophobic solid/water interface.