Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “Molecular Dynamic simulation”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,657 records · Page 92Linked to original sources
PubMed2000

X-Ray crystal structure and molecular dynamics simulations of silver hake parvalbumin (Isoform B).

Parvalbumins constitute a class of calcium-binding proteins characterized by the presence of several helix-loop-helix (EF-hand) motifs. In a previous study (Revett SP, King G, Shabanowitz J, Hunt DF, Hartman KL, Laue TM, Nelson DJ, 1997, Protein Sci 7:2397-2408), we presented the sequence of the major parvalbumin isoform from the silver hake (Merluccius bilinearis) and presented spectroscopic and structural information on the excised "EF-hand" portion of the protein. In this study, the X-ray crystal structure of the silver hake major parvalbumin has been determined to high resolution, in the frozen state, using the molecular replacement method with the carp parvalbumin structure as a starting model. The crystals are orthorhombic, space group C2221, with a = 75.7 A, b = 80.7 A, and c = 42.1 A. Data were collected from a single crystal grown in 15% glycerol, which served as a cryoprotectant for flash freezing at -188 degrees C. The structure refined to a conventional R-value of 21% (free R 25%) for observed reflections in the range 8 to 1.65 A [1 > 2sigma(I)]. The refined model includes an acetylated amino terminus, 108 residues (characteristic of a beta parvalbumin lineage), 2 calcium ions, and 114 water molecules per protein molecule. The resulting structure was used in molecular dynamics (MD) simulations focused primarily on the dynamics of the ligands coordinating the Ca2+ ions in the CD and EF sites. MD simulations were performed on both the fully Ca2+ loaded protein and on a Ca2+ deficient variant, with Ca2+ only in the CD site. There was substantial agreement between the MD and X-ray results in addressing the issue of mobility of key residues in the calcium-binding sites, especially with regard to the side chain of Ser55 in the CD site and Asp92 in the EF site.

Animals↗
PubMed1998

Molecular dynamics simulation of interaction of histone-like protein of mycobacterium tuberculosis (Hlpmt) and histone of clostridium pasteurianum (DBHclopa) with 35 based paired GC rich U-bend DNA.

Three dimensional structure of the first ninety two residues of histone like protein from Mycobacterium tuberculosis (Hlpmt) and histone of Clostridium pasteurianum (DBHclopa) are obtained here, on the basis of amino acid sequences of the two proteins, making use of secondary structure prediction programs, sequence search and HOMOLOGY based modeling tools available on Internet. The proteins were docked to a 35 base paired GC rich U bend DNA (U35DNA). Structures of proteins Hlpmt and DBHclopa; U35DNA; and complexes: Hlpmt-U35DNA and DBHclopa-U35DNA were optimized by molecular mechanics (MM) and simulated for 260 pico seconds (ps) in vacuum by molecular dynamics (MD) technique using AMBER 4.0 package with Cornell et al force field. The proteins, when simulated alone, showed compaction. DBHclopa showed larger compaction compared with Hlpmt. U35DNA when simulated alone straightened out and assumed a B-form. In the complexes, Hlpmt showed same order of compaction as in absence of DNA, while DBHclopa showed reduced compaction. In the presence of Hlpmt two ends of helicoidal axis of U35DNA came closer, but slightly out of plane, indicative of its role in overwinding and packaging double stranded DNA. DBHclopa did not give rise to DNA overwinding. The results show architectural role of Hlpmt and DBHclopa in DNA packaging and its sequence dependence.

Amino Acid Sequence↗
PubMed2001

Molecular dynamics simulations of the denaturation and refolding of an RNA tetraloop.

Tetraloops are very abundant structural elements of RNA that are formed by four nucleotides in a hairpin loop which is closed by a double stranded helical stem with some Watson-Crick base pairs. A tetraloop r(GCGAAGGC) was identified from the crystal structure of the central domain of 16S rRNA (727-730) in the Thermus thermophilus 30S ribosomal complex. The crystal structure of the 30S complex includes a total of 104 nucleotides from the central domain of the 16S rRNA and three ribosomal proteins S6, S15 and S18. Independent biochemical experiments have demonstrated that protein S15 plays the role in initiating the formation of the central domain of this complex. In the crystal, the tetraloop interacts with the protein S15 at two sites: one of them is associated with hydrogen bond interactions between residue His50 and nucleotide G730, and the other is associated with the occurrence of residue Arg53 beside A728. This paper uses molecular dynamics (MD) simulations to investigate the protein-dependent structural stability of the tetraloop and demonstrates the folding pathway of this tetraloop via melting denaturation and its subsequent refolding. Three important results are derived from these simulations: (i) The stability of nucleotide A728 appears to be protein dependent. Without the interaction with S15, A728 flips away from stacking with A729. (ii) The melting temperature demonstrated by the simulations is analogous to the results of thermodynamic experiments. In addition, the simulated folding of the tetraloop is stepwise: the native shape of the backbone is formed first; this is then followed by the formation of the Watson- Crick base pairs in the stem; and finally the hydrogen bonds and base stacking in the loop are formed. (iii) The tetraloop structure is similar to the crystal structure at salt concentrations of 0.1 M and 1.0 M used for the simulations, but the refolded structure at 0.1 M salt is more comparable to the crystal structure than at 1.0 M. The results from the simulations using both the Generalized Born continuum model and explicit solvent model (Particle Mesh Ewald) generate a similar pathway for unfolding/refolding of the tetraloop.

Arginine↗
PubMed2005

Solvent decompositions and physical properties of decomposition compounds in Li-ion battery electrolytes studied by DFT calculations and molecular dynamics simulations.

The density functional theory (DFT) calculations have been performed for the reduction decompositions of solvents widely used in Li-ion secondary battery electrolytes, ethylene carbonate (EC), propylene carbonate (PC), dimethyl carbonates (DMC), ethyl methyl carbonate (EMC), and diethyl carbonate (DEC), including a typical electrolyte additive, vinylene carbonate (VC), at the level of B3LYP/6-311+G(2d,p), both in the gas phase and solution using the polarizable conductor calculation model. In the gas phase, the first electron reduction for the cyclic carbonates and for the linear carbonates is found to be exothermic and endothermic, respectively, while the second electron reduction is endothermic for all the compounds examined. On the contrary, in solution both first and second electron reductions are exothermic for all the compounds. Among the solvents and the additive examined, the likelihood of undergoing the first electron reduction in solution was found in the order of EC > PC > VC > DMC > EMC > DEC with EC being the most likely reduced. VC, on the other hand, is most likely to undergo the second electron reduction among the compounds, in the order of VC > EC > PC. Based on the results, the experimentally demonstrated effectiveness of VC as an excellent electrolyte additive was discussed. The bulk thermodynamic properties of two dilithium alkylene glycol dicarbonates, dilithium ethylene glycol dicarbonate (Li-EDC) and dilithium 1,2-propylene glycol dicarbonate (Li-PDC), as the major component of solid-electrolyte interface (SEI) films were also examined through molecular dynamics (MD) simulations in order to understand the stability of the SEI film. It was found that film produced from a decomposition of EC, modeled by Li-EDC, has a higher density, more cohesive energy, and less solubility to the solvent than the film produced from decomposition of PC, Li-PDC. Further, MD simulations of the interface between the decomposition compound and graphite suggested that Li-EDC has more favorable interactions with the graphite surface than Li-PDC. The difference in the SEI film stability and the behavior of Li-ion battery cycling among the solvents were discussed in terms of the molecular structures.

Journal Article↗
PubMed2005

Ligand-binding regulation of LXR/RXR and LXR/PPAR heterodimerizations: SPR technology-based kinetic analysis correlated with molecular dynamics simulation.

Liver X receptor (LXR) and peroxisome proliferator-activated receptor (PPAR) are two members of nuclear receptors involved in the nutrient metabolisms of dietary fatty acid and cholesterol. They are found to be of cross-talk function in that LXR regulates fatty acid synthesis and PPAR controls fatty acid degradation. LXRs (LXRalpha and LXRbeta) function by forming obligate heterodimers with the retinoid X receptor (RXR), and subsequently binding to specific DNA response elements within the regulatory regions of their target genes. In this work, the kinetic features concerning LXR/RXR and LXR/PPAR interactions have been fully investigated using surface plasmon resonance (SPR) technology. It is found that LXRs could bind to all the three PPAR subtypes, PPARalpha, PPARgamma and PPARdelta with different binding affinities, and such receptor/receptor interactions could be regulated by ligand binding. Moreover, molecular dynamics (MD) simulations were performed on six typical complex models. The results revealed that ligands may increase the interaction energies between the receptor interfaces of the simulated receptor/receptor complexes. The MD results are in agreement with the SPR data. Further analyses on the MD results indicated that the ligand binding might increase the hydrogen bonds between the interfaces of the receptor/receptor complex.

DNA-Binding Proteins↗
PubMed1988

Microfolding: conformational probability map for the alanine dipeptide in water from molecular dynamics simulations.

A direct attack on the protein-folding problem has been initiated with the free energy perturbation methods of molecular dynamics. The complete conformational probability map for the alanine dipeptide is presented. This work uses the SPC model for the explicit hydration of the dipeptide. Free energy differences for the four observed minima (beta, alpha R, alpha L, C7ax) are given, and the free energy barriers between minima are outlined.

Dipeptides↗
PubMed2001

Molecular dynamics simulations on the mycotoxin fumonisin B1.

The solution conformational properties of the mycotoxin fumonisin B(1) have been studied using molecular dynamics methodology. Fumonisins have been shown to inhibit sphinganine (sphingosine) N-acyltransferase (ceramide synthase) and show a wide range of toxic effects in many animals. This study of the solution properties of fumonisin B(1) attempts to add to the structural models necessary for the understanding of the binding and activity properties. The computational method uses a box with periodic boundaries, filled with explicit TIP3P water molecules, the substrate fumonisin B(1), and selected counterions for charge neutrality. The starting structure of fumonisin B(1) is added to the box by excluding water molecules. The explicit image method using 12-A cutoffs is applied to the system and molecular dynamics are carried out on different starting conformations at 300 K in 100-picosecond (ps) steps. Examination of the resulting equilibrated conformations suggests that the structure is relatively extended and that previous computational studies in vacuo, showing a compact folded structure, may not be consistent with the solution structure.

Acyltransferases↗
PubMed2003

Insights from molecular dynamics simulations into pH-dependent enantioselective hydrolysis of ibuprofen esters by Candida rugosa lipase.

An interesting observation was found during our continued studies on the hydrolysis of ibuprofen esters by Candida rugosa lipase (CRL). An important role is played by pH in the stereospecific hydrolysis of these esters. The flap region of CRL plays a significant role in the access of the substrate to the active site of the enzyme. At pH 5.6, 48% of the methyl ester and 5% of the butyl ester of ibuprofen were hydrolysed in 5.5 h, whereas at pH 7.2, 9% of methyl ester and 45% of the butyl ester of ibuprofen was hydrolysed in a identical reaction time using CRL. This lead us to assume that CRL prefers the methyl ester of ibuprofen as a substrate at an acidic pH and the butyl ester of ibuprofen at a neutral pH. Therefore, in order to understand the role of pH in the substrate selection by CRL for the esters of ibuprofen we used the crystallographic coordinates of the open form of the CRL (1CRL) for molecular dynamics (MD) simulations under acidic and neutral conditions for 2 ns using GROMACS. The final structures obtained after simulation in acidic and neutral conditions were compared with the energy-minimized structure, and the root-mean-square deviations (r.m.s.ds) were calculated. The r.m.s.d. of the CRL flap at neutral pH was found to be greater than that of the CRL flap at acidic pH. The extent to which the flap opens at neutral pH allowed the bulkier substrate, the butyl ester of ibuprofen, to diffuse into the active site and provides the best enzyme-substrate fit for this specific substrate. At acidic pH there is a decreased opening of the flap thereby accommodating a more compact substrate, namely the methyl ester of ibuprofen. Thus, simulation experiments using MD provide reasonable insight for the pH-dependent substrate selectivity of CRL in aqueous environments.

Amino Acids↗
PubMed2006

Dynamic models of G-protein coupled receptor dimers: indications of asymmetry in the rhodopsin dimer from molecular dynamics simulations in a POPC bilayer.

Based on the growing evidence that G-protein coupled receptors (GPCRs) form homo- and hetero-oligomers, models of GPCR signaling are now considering macromolecular assemblies rather than monomers, with the homo-dimer regarded as the minimal oligomeric arrangement required for functional coupling to the G-protein. The dynamic mechanisms of such signaling assemblies are unknown. To gain some insight into properties of GPCR dimers that may be relevant to functional mechanisms, we study their current structural prototype, rhodopsin. We have carried out nanosecond time-scale molecular dynamics (MD) simulations of a rhodopsin dimer and compared the results to the monomer simulated in the same type of bilayer membrane model composed of an equilibrated unit cell of hydrated palmitoyl-oleoyl-phosphatidyl choline (POPC). The dynamic representation of the homo-dimer reveals the location of structural changes in several regions of the monomeric subunits. These changes appear to be more pronounced at the dimerization interface that had been shown to be involved in the activation process [Proc Natl Acad Sci USA 102:17495, 2005]. The results are consistent with a model of GPCR activation that involves allosteric modulation through a single GPCR subunit per dimer.

Animals↗
PubMed2003

Dissipative particle dynamics: a useful thermostat for equilibrium and nonequilibrium molecular dynamics simulations.

We discuss dissipative particle dynamics as a thermostat to molecular dynamics, and highlight some of its virtues: (i) universal applicability irrespective of the interatomic potential; (ii) correct and unscreened reproduction of hydrodynamic correlations; (iii) stabilization of the numerical integration of the equations of motion; and (iv) the avoidance of a profile bias in boundary-driven nonequilibrium simulations of shear flow. Numerical results on a repulsive Lennard-Jones fluid illustrate our arguments.

Journal Article↗
PubMed2004

Influence of a fluorobenzene nucleobase analogue on the conformational flexibility of RNA studied by molecular dynamics simulations.

Chemically modified bases are frequently used to stabilize nucleic acids, to study the driving forces for nucleic acid structure formation and to tune DNA and RNA hybridization conditions. In particular, fluorobenzene and fluorobenzimidazole base analogues can act as universal bases able to pair with any natural base and to stabilize RNA duplex formation. Although these base analogues are compatible with an A-form RNA geometry, little is known about the influence on the fine structure and conformational dynamics of RNA. In the present study, nano-second molecular dynamics (MD) simulations have been performed to characterize the dynamics of RNA duplexes containing a central 1'-deoxy-1'-(2,4-difluorophenyl)-beta-D-ribofuranose base pair or opposite to an adenine base. For comparison, RNA with a central uridine:adenine pair and a 1'-deoxy-1'-(phenyl)-beta-D-ribofuranose opposite to an adenine was also investigated. The MD simulations indicate a stable overall A-form geometry for the RNAs with base analogues. However, the presence of the base analogues caused a locally enhanced mobility of the central bases inducing mainly base pair shear and opening motions. No stable 'base-paired' geometry was found for the base analogue pair or the base analogue:adenine pairs, which explains in part the universal base character of these analogues. Instead, the conformational fluctuations of the base analogues lead to an enhanced accessibility of the bases in the major and minor grooves of the helix compared with a regular base pair.

Adenine↗
PubMed2002

Structure of docosahexaenoic acid-containing phospholipid bilayers as studied by (2)H NMR and molecular dynamics simulations.

Polyunsaturated phospholipids are known to be important with regard to the biological functions of essential fatty acids, for example, involving neural tissues such as the brain and retina. Here we have employed two complementary structural methods for the study of polyunsaturated bilayer lipids, viz. deuterium ((2)H) NMR spectroscopy and molecular dynamics (MD) computer simulations. Our research constitutes one of the first applications of all-atom MD simulations to polyunsaturated lipids containing docosahexaenoic acid (DHA; 22:6 cis-Delta(4,7,10,13,16,19)). Structural features of the highly unsaturated, mixed-chain phospholipid, 1-palmitoyl-2-docosahexaenoyl-sn-glycero-3-phosphocholine (PDPC), have been studied in the liquid-crystalline (L(alpha)) state and compared to the less unsaturated homolog, 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC). The (2)H NMR spectra of polyunsaturated bilayers are dramatically different from those of less unsaturated phospholipid bilayers. We show how use of MD simulations can aid in interpreting the complex (2)H NMR spectra of polyunsaturated bilayers, in conjunction with electron density profiles determined from small-angle X-ray diffraction studies. This work clearly demonstrates preferred helical and angle-iron conformations of the polyunsaturated chains in liquid-crystalline bilayers, which favor chain extension while maintaining bilayer flexibility. The presence of relatively long, extended fatty acyl chains may be important for solvating the hydrophobic surfaces of integral membrane proteins, such as rhodopsin. In addition, the polyallylic DHA chains have a tendency to adopt back-bended (hairpin-like) structures, which increase the interfacial area per lipid. Finally, the material properties have been analyzed in terms of the response of the bilayer to mechanical stress. Simulated bilayers of phospholipids containing docosahexaenoic acid were less sensitive to the applied surface tension than were saturated phospholipids, possibly implying a decrease in membrane elasticity (area elastic modulus, bending rigidity). The above features distinguish DHA-containing lipids from saturated or monounsaturated lipids and may be important for their biological modes of action.

Computer Simulation↗
PubMed2004

Molecular dynamics simulations of papilloma virus E2 DNA sequences: dynamical models for oligonucleotide structures in solution.

The specificity of papilloma virus E2 protein-DNA binding depends critically upon the sequence of a region of the DNA not in direct contact with the protein, and represents one of the simplest known examples of indirect readout. A detailed characterization of this system in solution is important to the further investigation hypothesis of a structural code for DNA recognition by regulatory proteins. In the crystalline state, the E2 DNA oligonucleotide sequence, d(ACCGAATTCGGT), exhibits three different structural forms. We report herein studies of the structure of E2 DNA in solution based on a series of molecular dynamics (MD) simulations including counterions and water, utilizing both the canonical and various crystallographic structures as initial points of departure. All MDs converged on a single dynamical structure of d(ACCGAATTCGGT) in solution. The predicted structure is in close accord with two of the three crystal structures, and indicates that a significant kink in the double helix at the central ApT step in the other crystal molecule may be a packing effect. The dynamical fine structure was analyzed on the basis of helicoidal parameters. The calculated curvature in the sequence was found to originate primarily from YPR steps in the regions flanking the central AATT tract. In order to study the role of structural adaptation of the DNA in the binding process, a subsequent simulation on the 16-mer cognate sequence d(CAACCGAATTCGGTTG) was initiated from the crystallographic coordinates of the bound DNA in the crystal structure of the protein DNA complex. MD simulations starting with the protein-bound form relaxed rapidly back to the dynamical structure predicted from the previous simulations on the uncomplexed DNA. The MD results show that the bound form E2 DNA is a dynamically unstable structure in the absence of protein, and arises as a consequence of both structural changes intrinsic to the sequence and induced by the interaction with protein.

Animals↗
PubMed2006

Toward understanding the inactivation mechanism of monooxygenase P450 BM-3 by organic cosolvents: a molecular dynamics simulation study.

Cytochrome P450 BM-3 from Bacillus megaterium is an extensively studied enzyme for industrial applications. A major focus of current protein engineering research is directed to improving the catalytic performance of P450 BM-3 toward nonnatural substrates of industrial importance in the presence of organic solvents or cosolvents. For the latter reason, it is important to study the effect of organic cosolvent molecules on the structure and dynamics of the enzyme, in particular, the effect of cosolvent molecules on the active site's structure and dynamics. In this paper, we have studied, using molecular dynamics (MD) simulations, the F87A mutant of P450 BM-3 in the presence of DMSO as cosolvent, to understand the role of the F87A substitution for its catalytic activity. This mutant exhibits an altered regioselectivity and substrate specificity compared with wild-type; however, it has lower tolerance toward DMSO. The simulation results offer an explanation for the DMSO sensitivity of the F87A mutant. Our simulation results show that the F87 side chain prevents the disturbance of the water molecule bound to the heme iron by DMSO molecules. The absence of the phenyl ring in F87A mutant promotes interactions of the DMSO molecule with the heme iron resulting in water displacement by DMSO at the catalytic heme center.

Bacillus megaterium↗
PubMed1986

Molecular dynamics simulations of the holo and apo forms of retinol binding protein. Structural and dynamical changes induced by retinol removal.

The effects of removing retinol from the X-ray structure of holo-retinol binding protein are studied using the molecular dynamics technique. Structural and dynamical properties emerging from an 80 ps simulation of the apo form, for which no crystallographic structure is available, are compared with the results of a 70 ps trajectory of the holo-protein. Dynamical stationarity is attained after roughly 30 ps, and the resulting average structure is proposed as a reasonable model of the apo-protein. Conformational changes are observed for the loops at the beta-barrel entrance during the non-equilibrium part of the apo-trajectory. Tryptophan labelling experiments and retinoid reconstitution experiments point towards this part of the molecule as being involved in prealbumin binding. Structural changes in this region may therefore explain the differences in prealbumin affinity between the apo and holo forms. Furthermore, a change in the position of the alpha-helix, corresponding to a pivot around its C terminus, is observed for the apo-protein. The resulting conformation of the alpha-helix is found to be similar to that in apo-beta-lactoglobulin, which also can bind retinol and for which a crystal structure exists. The results from the holo simulation are compared to the crystallographic data and show good agreement. The dynamics of the secondary and tertiary structural elements are analysed and compared for the two forms. The beta-barrel is found to be extremely cooperative in its atomic motions in both simulations, and the top and bottom beta-sheets perform collective fluctuations with respect to each other in the low-frequency limit of the simulations. The dynamics of the alpha-helical region presents clear differences between the two forms; while the holo-protein has a well-defined spectrum for the longitudinal stretching mode, the apo form displays a fairly large bending of the alpha-helix at several points of the trajectory.

Computer Simulation↗
PubMedDate not supplied

Molecular dynamics simulations and experimental studies of binding and mobility of 7-tert-butyldimethylsilyl-10-hydroxycamptothecin and its 20(S)-4-aminobutyrate ester in DMPC membranes.

The enhanced permeability and retention of liposomes in solid tumors makes liposomal formulations attractive for the targeting of various antitumor agents. This study explores the binding, orientation, and dynamic properties of a potent topoisomerase I inhibitor, 7-tert-butyldimethylsilyl-10-hydroxycamptothecin (DB-67), and its 20(S)-4-aminobutyrate ester prodrug (DB-67-AB) in DMPC liposomes by molecular dynamics (MD) simulations and experimental studies. MD simulations of an all-atom and fully hydrated liquid-crystalline bilayer (2 x 36 DMPC lipids) containing single molecules of DB-67 and DB-67-AB were conducted for up to 50 ns. Membrane/water partition coefficients for DB-67 and DB-67-AB vs pH were determined by ultracentrifugation. Fluorescence spectra and/or steady-state anisotropies were measured in various solvents and in DMPC liposomes. Kinetics for the reversible DB-67 lactone ring-opening in the presence and absence of DMPC liposomes were determined by HPLC with fluorescence detection. During the entire simulation time both DB-67 and DB-67-AB were located on the bilayer membrane near the polar ester groups of DMPC. The average depth of penetration for DB-67 and DB-67-AB was similar (12.4-13.2 A) with the prodrug's protonated amino group strongly solvated by surface water and lipid phosphate groups. Binding and fluorescence experiments revealed only a modest reduction in the binding affinity upon attachment of the ionized 4-aminobutyrate group onto DB-67. The binding microenvironment polarity resembles that of a polar solvent such as EtOH and DMSO. Kinetics experiments confirmed that DB-67 lactone hydrolysis is inhibited in the presence of DMPC liposomes, consistent with the reduced exposure of its lactone ring to water, as observed in the simulations. Both bound DB-67 and bound DB-67-AB have nonrandom orientations and reduced mobility in the membrane, especially for diffusion normal to the bilayer surface, and rotational relaxation, both of which are > or =2 orders of magnitude slower than in bulk water. MD simulations correctly predicted the high binding affinities for DB-67-AB to DMPC bilayers, protection of bound DB-67 toward lactone hydrolysis, and the lack of a substantial reduction in binding for the 20(S)-4-aminobutyrate prodrug of DB-67.

Antineoplastic Agents↗
PubMed2004

Molecular dynamics simulations of gas separations using faujasite-type zeolite membranes.

Gas separations with faujasite zeolite membranes have been examined using the method of molecular dynamics. Two binary mixtures are investigated, oxygen/nitrogen and nitrogen/carbon dioxide. These mixtures have been found experimentally to exhibit contrasting behavior. In O(2)/N(2) mixtures the ideal selectivity (pure systems) is higher than the mixture selectivity, while in N(2)/CO(2) the mixture selectivity is higher than the ideal selectivity. One of the key goals of this work was to seek a fundamental molecular level understanding of such divergent behavior. Our simulation results (using previously developed intermolecular models for both the gases and zeolites investigated) were found to replicate this experimental behavior. By examining the loading of the membranes and the diffusion rates inside the zeolites, we have been able to explain such contrasting behavior of O(2)/N(2) and N(2)/CO(2) mixtures. In the case of O(2)/N(2) mixtures, the adsorption and loading of both O(2) and N(2) in the membrane are quite competitive, and thus the drop in the selectivity in the mixture is primarily the result of oxygen slowing the diffusion of nitrogen and nitrogen somewhat increasing the diffusion of oxygen when they pass through the zeolite pores. In N(2)/CO(2) systems, CO(2) is rather selectively adsorbed and loaded in the zeolite, leaving very little room for N(2) adsorption. Thus although N(2) continues to have a higher diffusion rate than CO(2) even in the mixture, there are so few N(2) molecules in the zeolite in mixtures that the selectivity of the mixture increases significantly compared to the ideal (pure system) values. We have also compared simulation results with hydrodynamic theories that classify the permeance of membranes to be either due to surface diffusion, viscous flow, or Knudsen diffusion. Our results show surface diffusion to be the dominant mode, except in the case of N(2)/CO(2) binary mixtures where Knudsen diffusion also makes a contribution to N(2) transport.

Journal Article↗
PubMed2006

Free energy surfaces of miniproteins with a betabetaalpha motif: replica exchange molecular dynamics simulation with an implicit solvation model.

Designed miniproteins with a betabetaalpha motif, such as BBA5, 1FSD, and 1PSV can serve as a benchmark set to test the validity of all-atom force fields with computer simulation, because they contain all the basic structural elements in protein folding. Unfortunately, it was found that the standard all-atom force fields with the generalized Born (GB) implicit solvation model tend to produce distorted free energy surfaces for the betabetaalpha proteins, not only because energetically those proteins need to be described by more balanced weights of the alpha- and beta-strands, but also because the GB implicit solvation model suffers from overestimated salt bridge effects. In an attempt to resolve these problems, we have modified one of the standard all-atom force fields in conjunction with the GB model, such that each native state of the betabetaalpha proteins is in its free energy minimum state with reasonable energy barriers separating local minima. With this modified energy model, the free energy contour map in each protein was constructed from the replica exchange molecular dynamics REMD simulation. The resulting free energy surfaces are significantly improved in comparison with previous simulation results and consistent with general views on small protein folding behaviors with realistic topology and energetics of all three proteins.

Calorimetry↗