Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “molecular dynamics”

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 127 records · Page 7Linked to original sources
PubMed2003

Mechanisms of selectivity in channels and enzymes studied with interactive molecular dynamics.

Interactive molecular dynamics, a new modeling tool for rapid investigation of the physical mechanisms of biological processes at the atomic level, is applied to study selectivity and regulation of the membrane channel protein GlpF and the enzyme glycerol kinase. These proteins facilitate the first two steps of Escherichia coli glycerol metabolism. Despite their different function and architecture the proteins are found to employ common mechanisms for substrate selectivity: an induced geometrical fit by structurally homologous binding sites and an induced rapid dipole moment reversal. Competition for hydrogen bonding sites with water in both proteins is critical for substrate motion. In glycerol kinase, it is shown that the proposed domain motion prevents competition with water, in turn regulating the binding of glycerol.

Aquaporins↗
PubMed1997

Role of water on unfolding kinetics of helical peptides studied by molecular dynamics simulations.

Molecular dynamics simulations have been carried out with four polypeptides, Ala13, Val(13), Ser13, and Ala4Gly5Ala4, in vacuo and with explicit hydration. The unfolding of the polypeptides, which are initially fully alpha-helix in conformation, has been monitored during trajectories of 0.3 ns at 350 K. A rank of Ala < Val < Ser < Gly is found in the order of increasing rate of unwinding. The unfolding of Ala13 and Val(13) is completed in hundreds of picoseconds, while that of Ser13 is about one order of magnitude faster. The helix content of the peptide containing glycine residues falls to zero within a few picoseconds. Ramachandran plots indicate quite distinct equilibrium distributions and time evolution of dihedral angles in water and in vacuum for each residue type. The unfolding of polyalanine and polyvaline helices is accelerated due to solvation. In contrast, polyserine is more stable in water compared to vacuum, because its side chains can form intramolecular hydrogen bonds with the backbone more readily in vacuum, which disrupts the helix. Distribution functions of the spatial and angular position of water molecules in the proximity of the polypeptide backbone polar groups reveal the stabilization of the coiled structures by hydration. The transition from helix to coil is characterized by the appearance of a new peak in the probability distribution at a specific location characteristic of hydrogen bond formation between water and backbone polar groups. No significant insertion of water molecules is observed at the precise onset of unwinding, while (i, i+3) hydrogen bond formation is frequently detected at the initiation of alpha-helix unwinding.

Biophysical Phenomena↗
PubMed2003

Conformations of terminal sialyloligosaccharide fragments--a molecular dynamics study.

Molecular dynamics simulations have been performed to understand the conformational features of the terminal sialyloligosaccharide fragments NeuNAc alpha(2-3)Gal, NeuNAc alpha(2-6)Gal, NeuNAc alpha(2-8)NeuNAc and NeuNAc alpha(2-9)NeuNAc. The conformational regions A(i), B(i) and C(i) were identified in the Ramachandran plot. Analysis of the 1000 ps trajectories collected through simulation (2000 ps in the case of NeuNAc alpha (2-9)NeuNAc) revealed that these molecules have conformational propensity in region B(i). The occurrence of these molecules in the common conformational space leads to a structural similarity between them. This structural similarity may be an essential requirement for the neuraminidase activity towards sialyloligosaccharides. The local change in the conformation of the active site residues of neuraminidases may contribute for the specificity differences between different linkages of sialyloligosaccharides. A highly conserved water-mediated hydrogen bond observed in these structures between the sugar residues, acts as an additional stabilizing force.

Carbohydrate Conformation↗
PubMed1996

Temperature dependence of the redox potential of rubredoxin from Pyrococcus furiosus: a molecular dynamics study.

Molecular dynamics simulations are used to evaluate the temperature dependent differences in structure, solvation, and energies for the iron-sulfur protein rubredoxin from the hyperthermophilic archebacterium Pyrococcus furiosus to understand the unusual temperature dependence of its redox potential [Adams, M. W. W. (1992) Adv. Inorg. Chem. 38, 341-396]. Simulations of both redox states performed at 295 and 363 K reveal that almost no backbone structure alteration occurs at the higher temperature and that the radius of gyration of the protein is temperature and redox state independent. The most striking change is that the penetration of the redox site by solvent molecules in the reduced from at 295 K, which was also seen in simulations of the reduced form of the mesophilic Clostridium pasteurianum rubredoxin at 295 K (Yelle, R. B., et al. (1995) Proteins 22, 154-167], is no longer seen to a significant extent in either redox state at 363 K. Comparing 295 to 363 K, the calculated change in the electrostatic potential of about -300 mV and in the negative of the potential energy of about -550 meV is consistent with the observed change in redox potential of -160 mV. Moreover, the calculated change is in the wrong direction if the penetrating water is excluded. These results show that changing solvent accessibility may be responsible for the temperature dependence of the redox potential of P. furiosus rubredoxin.

Archaea↗
PubMed1999

Molecular Shape and Solvation of the Lacunar, Saddle-Shaped, and Planar Metal Cyclidene Complexes: Molecular Dynamics Studies.

Molecular dynamics simulations have been used to study the three-dimensional distribution of methanol solvent molecules around three cobalt(II) cyclidene complexes differing in details of their ligand structures. The ligands are a planar unbridged 14-membered macrocycle in Co([14]Cyc), a saddle-shaped unbridged 16-membered macrocycle in Co([16]Cyc), and a lacunar bridged 16-membered macrocycle in Co(C6[16]Cyc). All three complexes contain five-coordinate cobalt(II) with the metal ion bound to four nitrogen donor atoms from the macrocycle and one nitrogen donor from an axial methylimidazole. Distinctly different solvation patterns are exemplified for the three complexes by the positions of the maximum in the Co-O pair distribution function (at r(Co-O) = 2.5, 3.7, and 4.5 Å for Co([14]Cyc), Co([16]Cyc), and Co(C6[16]Cyc), respectively) and by the number of methanol molecules within the macrocyclic cleft (1.5, 0.7, and 0.4 molecules at a Co-O distance of 5.25 Å in the "cavity", respectively). Analysis of the anisotropic solvent structure reveals the presence of a methanol molecule directly above the cobalt(II) center, at a distance of ca. 2.5 Å, for planar Co([14]Cyc), and the absence of solvent from such close proximity to the metal ion for the remaining complexes. The bridge further protects the sides of the cavity from the solvent. The width of an empty cavity of Co(C6[16]Cyc) shrinks by 0.3 Å in methanol solution, as compared to vacuum simulations. These results confirm the experimentally based (decrease in absolute value of enthalpies and entropies of dioxygen binding) suggestion that extensive solvation of the cobalt(II) center reduces its accessibility to incoming small molecules.

Journal Article↗
PubMed2002

Microwave-driven zeolite-guest systems show athermal effects from nonequilibrium molecular dynamics.

Nonequilibrium molecular dynamics simulations show that steady-state systems obtained by microwave heating are qualitatively different from those at thermal equilibrium. This difference arises because energy transfer from hotter to colder species is not efficient enough to equilibrate the distribution of energy. Under nonequilibrium conditions, we found that microwave radiation can selectively heat methanol in a binary mixture of methanol-benzene adsorbed in faujasite zeolite. The difference in steady-state temperatures follows the trend Tmethanol > Tbenzene > Tzeolite, which is qualitatively consistent with recent experimental results.

Journal Article↗
PubMed2005

Liquid-vapor interface of methanol-water mixtures: a molecular dynamics study.

Molecular dynamics simulations were carried out to investigate the structural and thermodynamic properties and variations in the dipole moments of the liquid-vapor interfaces of methanol-water mixtures. Various methanol-water compositions were simulated at room temperature. We found that methanol tends to concentrate at the interface, and the computed surface tension shows a composition dependence that is consistent with experimental measurements. The methanol molecule shows preferred orientation near the interface with the methyl group pointing into the vapor phase. The methanol in the mixture is found to have larger dipole moments than that of pure liquid methanol. The strong local field induced by the surrounding water molecules is partly the reason for this difference. The dependence of hydrogen-bonding patterns between methanol and water on the interface and the composition of the mixture is also discussed in the paper.

Computer Simulation↗
PubMed2005

Cation mobility and the sorption of chloroform in zeolite NaY: molecular dynamics study.

Molecular dynamics simulations at temperatures of 270, 330, and 390 K have been carried out to address the question of cation migration upon chloroform sorption in sodium zeolite Y. The results show that sodium cations located in different sites exhibit different types of mobility. These may be summarized as follows: (1) SII cations migrate toward the center of the supercage upon sorption, due to interactions with the polar sorbate molecules. (2) SI' cations hop from the sodalite cage into the supercage to fill vacant SII sites. (3) SI' cations migrate to other SI' sites within the same sodalite cage. (4) SI cations hop out of the double six-rings into SI' sites. In some instances, concerted motion of cations is observed. Furthermore, former SI' and SI cations, having crossed to SII sites, may then further migrate within the supercage, as in (1). The cation motion is dependent on the level of sorbate loading, with 10 molecules per unit cell not being enough to induce significant cation displacements, whereas the sorption of 40 molecules per unit cell results in a number of cations being displaced from their original positions. Further rearrangement of the cation positions is observed upon evacuation of the simulation cell, with some cations reverting back to sites normally occupied in bare NaY.

Journal Article↗
PubMed2005

Characterization of titanium dioxide nanoparticles using molecular dynamics simulations.

Molecular dynamics simulations of titanium dioxide nanoparticles in the three commonly occurring phases (anatase, brookite, and rutile) are reported. The structural properties inferred by simulated X-ray diffraction patterns of the nanoparticles were investigated. The titanium-oxygen bond length as a function of size, phase, and temperature was determined and was found to be dependent on the coordination environment of the titanium and independent of phase and size. The equilibrium Ti-O bond length is 1.86 A for a four-coordinated titanium ion, 1.92 A for a five-coordinated titanium ion, and 1.94 A for an octahedral titanium ion. Smaller nanoparticles are characterized by a higher fraction of titanium ions that are four and five coordinated, due to the larger surface area-to-volume ratios. The surface energies for anatase, rutile, and brookite particles were reported. The surface energy of the nanoparticle increases and approaches a constant value as the particle gets bigger. The surface energies of small rutile particles are higher than that for anatase particles of a similar size, consistent with anatase being the more stable phase of nanocrystalline titanium dioxide.

Journal Article↗
PubMed2005

Structure of an associating polymer melt in a narrow slit by molecular dynamics simulation.

Molecular dynamics simulation has been used to study the equilibrium properties of a generic coarse-grained polymer melt with associating terminal groups, confined in a narrow slit by two atomically smooth walls. Simulations were carried out as a function of wall separation and attracting strength as well as polymer end-end interaction strength. We find that confinement has an important effect on the melt properties. In particular, strongly attracting walls can produce radical changes in chain conformation, the nature of the transient network, and the structure of the aggregates formed by the associating terminals.

Journal Article↗
PubMed2006

Structure of the liquid-vacuum interface of room-temperature ionic liquids: a molecular dynamics study.

Molecular dynamics simulations for the liquid-vacuum interface of the ionic liquid 1-ethyl-3-methylimidazolium nitrate (EMIM+/NO3-) were performed for both electronically polarizable and nonpolarizable potential energy surfaces. The interfacial structural properties, such as the oscillation in the number density profile, the orientational ordering, and the local clustering in the interfacial region, were calculated. The simulations with both the polarizable and nonpolarizable model demonstrate the existence of an inhomogeneous interfacial structure normal to the surface layer. It was found for both models that the ethyl tail group on EMIM+ is likely to protrude outward from the surface. In the outmost surface layer, the cation is likely to lie on the surface with the imidazolium ring parallel to the interface, while there is a second region with enhanced density from that in the bulk where the cation preferably slants with the imidazolium ring tending to be perpendicular to the surface. The results also reveal that the electronic polarization effect is important for the ionic liquid interface. It is found that the cation is likely to be segregated at the ionic liquid surface for the polarizable model, while for the nonpolarizable model, the anion is found to be more likely to exhibit such behavior. The surface tension of the polarizable model (58.5 +/- 0.5 mN/m) is much smaller than that of the nonpolarizable model (82.7 +/- 0.6 mN/m), in better agreement with extrapolated experimental measurements on similar ionic liquid systems.

Journal Article↗
PubMed2006

LiTFSI structure and transport in ethylene carbonate from molecular dynamics simulations.

Molecular dynamics (MD) simulations using a many-body polarizable force field were performed on ethylene carbonate (EC) doped with lithium bistrifluoromethanesulfonamide (LiTFSI) salt as a function of temperature and salt concentration. At 313 K Li+ was coordinated by 2.7-3.2 EC carbonyl oxygen atoms and 0.67-1.05 TFSI- oxygen atoms at EC:Li = 10 and 20 salt concentrations. In completely dissociated electrolytes, however, Li+ was solvated by approximately 3.8 carbonyl oxygen atoms from EC on average. The probability of ions to participate ion aggregates decreased exponentially with an increase in the size of the aggregate. Ion and solvent self-diffusion coefficients and conductivity predicted by MD simulations were in good agreement with experiments. Approximately half of the charge was transported by charged ion aggregates with the other half carried by free (uncomplexed by counterion) ions. Investigation of the Li+ transport mechanism revealed that contribution from the Li+ diffusion together with its coordination shell to the total Li+ transport is similar to the contribution arising from Li+ exchanging solvent molecules in its first coordination shell with solvents from the outer shells.

Journal Article↗
PubMed2006

Preferential solvation of phenol in binary solvent mixtures. A molecular dynamics study.

Molecular dynamics computer simulations were carried out to study the preferential solvation of phenol in equimolar acetonitrile-water and ethanol-water binary mixtures. Two water models were used to investigate the model dependence of preferential solvation. The results are compared to recent intermolecular 1H NOESY experiments reported on the same systems. In the case of acetonitrile-water the local mole fraction obtained from simulations agrees quite well with experiments. In the case of ethanol-water there was a qualitative difference, which was observed for both water models. However, when comparing the degree of preferential solvation of the two cosolvents ethanol and acetonitrile with each of the two water models, the trend obtained from the simulations agrees with experimental data.

Journal Article↗
PubMed2006

Interaction of the sugars trehalose, maltose and glucose with a phospholipid bilayer: a comparative molecular dynamics study.

Molecular dynamics simulations are used to investigate the interaction of the sugars trehalose, maltose, and glucose with a phospholipid bilayer at atomic resolution. Simulations of the bilayer in the absence or in the presence of sugar (2 molal concentration for the disaccharides, 4 molal for the monosaccharide) are carried out at 325 and 475 K. At 325 K, the three sugars are found to interact directly with the lipid headgroups through hydrogen bonds, replacing water at about one-fifth to one-quarter of the hydrogen-bonding sites provided by the membrane. Because of its small size and of the reduced topological constraints imposed on the hydroxyl group locations and orientations, glucose interacts more tightly (at identical effective hydroxyl group concentration) with the lipid headgroups when compared to the disaccharides. At high temperature, the three sugars are able to prevent the thermal disruption of the bilayer. This protective effect is correlated with a significant increase in the number of sugar-headgroups hydrogen bonds. For the disaccharides, this change is predominantly due to an increase in the number of sugar molecules bridging three or more lipid molecules. For glucose, it is primarily due to an increase in the number of sugar molecules bound to one or bridging two lipid molecules.

Computer Simulation↗
PubMed2006

Cation migration upon adsorption of methanol in NaY and NaX faujasite systems: a molecular dynamics approach.

Molecular dynamics simulations have been carried out to address the question of cation migration upon adsorption of methanol in NaY and NaX faujasite systems as a function of the loading. For NaY, it has been shown that, at low and intermediate loadings, SII cations can migrate toward the center of the supercage due to strong interactions with the adsorbates, followed by a hopping of SI' from the sodalite cage into the supercage to fill the vacant SII site. A SI' cation can also migrate across the double six ring and takes a SI' vacant position. SI cations mainly remain trapped in their initial sites whatever the loading. At high loading, only limited motions are observed for SII cations due to steric effects induced by the presence of adsorbates within the supercage. For NaX, the SIII' cations which occupy the most accessible adsorption sites are significantly moving upon coordination to the methanol molecules; the extent of this mobility exhibits a maximum for 48 methanol molecules per unit cell before decreasing at higher loadings due to steric hindrance. In addition, the SI' and SII cations remain almost trapped in their initial sites whatever the loading. Indeed, the most probable migration mechanism involves SIII' cation displacements into nearby SIII' sites.

Journal Article↗
PubMed2004

The structure of phosphine-functionalised silsesquioxane-based dendrimers: a molecular dynamics study.

Molecular dynamics modelling has been used to simulate the structures of phopshine-functionalised, silsesquioxane-based dendrimers whose excellent catalytic properties have been previously demonstrated. The effect of changing the chemical composition of the dendrimer branches is simulated. The results indicate that adding a methylene unit to a branch increases the overall size of the dendrimer but replacing a methylene unit with an oxygen atom decreases the size of the dendrimer. The size and shape of the dendrimers have also been simulated on changing the temperature and polarity of the solvent. The distribution of phosphine groups on the exterior has also been modelled and this suggests that all are available for bonding to catalytic metals in all the compounds.

Journal Article↗
PubMed2005

A molecular dynamics computer simulation study of room-temperature ionic liquids. II. Equilibrium and nonequilibrium solvation dynamics.

The molecular dynamics (MD) simulation study of solvation structure and free energetics in 1-ethyl-3-methylimidazolium chloride and 1-ethyl-3-methylimidazolium hexafluorophosphate using a probe solute in the preceding article [Y. Shim, M. Y. Choi and H. J. Kim, J. Chem. Phys. 122, 044510 (2005)] is extended to investigate dynamic properties of these liquids. Solvent fluctuation dynamics near equilibrium are studied via MD and associated time-dependent friction is analyzed via the generalized Langevin equation. Nonequilibrium solvent relaxation following an instantaneous change in the solute charge distribution and accompanying solvent structure reorganization are also investigated. Both equilibrium and nonequilibrium solvation dynamics are characterized by at least two vastly different time scales--a subpicosecond inertial regime followed by a slow diffusive regime. Solvent regions contributing to the subpicosecond nonequilibrium relaxation are found to vary significantly with initial solvation configurations, especially near the solute. If the solvent density near the solute is sufficiently high at the outset of the relaxation, subpicosecond dynamics are mainly governed by the motions of a few ions close to the solute. By contrast, in the case of a low local density, solvent ions located not only close to but also relatively far from the solute participate in the subpicosecond relaxation. Despite this difference, linear response holds reasonably well in both ionic liquids.

Journal Article↗
PubMed2005

Melting of icosahedral gold nanoclusters from molecular dynamics simulations.

Molecular dynamics simulations show that gold clusters with about 600-3000 atoms crystallize into a Mackay icosahedron upon cooling from the liquid. A detailed surface analysis shows that the facets on the surface of the Mackay icosahedral gold clusters soften but do not premelt below the bulk melting temperature. This softening is found to be due to the increasing mobility of vertex and edge atoms with temperature, which leads to inter-layer and intra-layer diffusion, and a shrinkage of the average facet size, so that the average shape of the cluster is nearly spherical at melting.

Journal Article↗