SEARCH · PubMed Health

Results for “Molecular Dynamic simulation”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

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

At least 253 records · Page 14Linked to original sources

PubMed1993

Molecular dynamics simulations of heme reorientational motions in myoglobin.

Molecular dynamics simulations of 2-ns duration were performed on carbonmonoxymyoglobin and deoxymyoglobin in vacuo to study the reorientational dynamics of the heme group. The heme in both simulations undergoes reorientations of approximately 5 degrees amplitude on a subpicosecond time scale, which produce a rapid initial decay in the reorientational correlation function to about 0.99. The heme also experiences infrequent changes in average orientation of approximately 10 degrees amplitude, which lead to a larger slow decay of the reorientational correlation function over a period of hundreds of picoseconds. The simulations have not converged with respect to these infrequent transitions. However, an estimate of the order parameter for rapid internal motions of the heme from those orientations which are sampled by the simulations suggests that the subnanosecond orientational dynamics of the heme accounts for at least 30% of the unresolved initial anisotropy decay observed in the nanosecond time-resolved optical absorption experiments on myoglobin reported by Ansari et al. in a companion paper (Ansari, A., C.M. Jones, E.R. Henry, J. Hofrichter, and W.A. Eaton. 1992. Biophys. J. 64:852-868.). A more complete sampling of the accessible heme orientations would most likely increase this fraction further. The simulation of the liganded molecule also suggests that the conformational dynamics of the CO ligand may contribute significantly to discrepancies between the ligand conformation as probed by x-ray diffraction and by infrared-optical photoselection experiments. The protein back-bone explores multiple conformations during the simulations, with the largest structural changes appearing in the E and F helices, which are in contact with the heme. The variations in the heme orientation correlate with the conformational dynamics of the protein on a time scale of hundreds of picoseconds, suggesting that the heme orientation may provide a useful probe of dynamical processes in the protein.

Animals↗

PubMed2005

Molecular dynamics simulation of swollen membrane of perfluorinated ionomer.

Molecular dynamics simulations of the swollen membrane of perfluorinated ionomer, which is composed of poly(tetrafluoroethylene) backbones and perfluosulfonic pendant side chains, have been undertaken to analyze the static and dynamic properties of the water and the side chain in the membrane. The calculations were carried out for four different water contents, 5, 10, 20 and 40 wt %, at 358.15 K and 0.1 MPa. The results are summarized as follows: (1) The sulfonic acid is the unique site to which water molecules can bind, and the other sites in the pendant side chain have no bound water even at high water concentration. (2) Sulfonic acids aggregate in the short range within 4.6-7.7 A despite the electrostatic repulsion between them. In such aggregates, a water molecule bridges two sulfonic acids. (3) Pendant side chains prefer to orient perpendicular to the hydrophilic/hydrophobic interface, and long-range correlation of side chain orientations is observed at 20 and 40 wt % water uptake membranes. (4) In a low water uptake membrane, the dynamics of water is substantially restricted due to strong attractive interactions with acidic sites. In contrast, at high water content, even the water locating near the sulfonic acid is relatively mobile. The short residence time of the bound water reveals that such water can frequently exchange position with relatively free water, which locates in the center of water cluster, in highly swollen membranes.

Journal Article↗

PubMed2006

Molecular dynamics simulations of the melting of aluminum nanoparticles.

Molecular dynamics simulations are performed to determine the melting points of aluminum nanoparticles of 55-1000 atoms with the Streitz-Mintmire [Phys. Rev. B 1994, 50, 11996] variable-charge electrostatic plus potential. The melting of the nanoparticles is characterized by studying the temperature dependence of the potential energy and Lindemann index. Nanoparticles with less than 850 atoms show bistability between the solid and liquid phases over temperature ranges below the point of complete melting. The potential energy of a nanoparticle in the bistable region alternates between values corresponding to the solid and liquid phases. This bistability is characteristic of dynamic coexistence melting. At higher temperatures, only the liquid state is stable. Nanoparticles with more than 850 atoms undergo a sharp solid-liquid-phase transition characteristic of the bulk solid phase. The variation of the melting point with the effective nanoparticle radius is also determined.

Journal Article↗

PubMed2004

Molecular dynamic simulation methods for anisotropic liquids.

Methods of molecular dynamics simulations for anisotropic molecules are presented. The new methods, with an anisotropic factor in the cell dynamics, dramatically reduce the artifacts related to cell shapes and overcome the difficulties of simulating anisotropic molecules under constant hydrostatic pressure or constant volume. The methods are especially effective for anisotropic liquids, such as smectic liquid crystals and membranes, of which the stacks of layers are compressible (elastic in direction perpendicular to the layers) while the layer itself is liquid and only elastic under uniform compressive force. The methods can also be used for crystals and isotropic liquids as well.

Journal Article↗

PubMed2005

Molecular-dynamics simulations of methane hydrate dissociation.

Nonequilibrium molecular-dynamics simulations have been carried out at 276.65 K and 68 bar for the dissolution of spherical methane hydrate crystallites surrounded by a liquid phase. The liquid was composed of pure water or a water-methane mixture ranging in methane composition from 50% to 100% of the corresponding theoretical maximum for the hydrate and ranged in size from about 1600 to 2200 water molecules. Four different crystallites ranging in size from 115 to 230 water molecules were used in the two-phase systems; the nanocrystals were either empty or had a methane occupation from 80% to 100% of the theoretical maximum. The crystal-liquid systems were prepared in two distinct ways, involving constrained melting of a bulk hydrate system or implantation of the crystallite into a separate liquid phase. The breakup rates were very similar for the four different crystal sizes investigated. The method of system preparation was not found to affect the eventual dissociation rates, despite a lag time of approximately 70 ps associated with relaxation of the liquid interfacial layer in the constrained melting approach. The dissolution rates were not affected substantially by methane occupation of the hydrate phase in the 80%-100% range. In contrast, empty hydrate clusters were found to break up significantly more quickly. Our simulations indicate that the diffusion of methane molecules to the surrounding liquid layer from the crystal surface appears to be the rate-controlling step in hydrate breakup. Increasing the size of the liquid phase was found to reduce the initial delay in breakup. We have compared breakup rates computed using different long-range electrostatic methods. Use of the Ewald, minimum image, and spherical cut-off techniques led to more rapid dissociation relative to the Lekner method.

Journal Article↗

PubMed2006

Structure change of liquid GaSb under pressure: an ab initio molecular-dynamics simulation.

We have performed ab initio molecular-dynamics simulation of liquid GaSb (l-GaSb) up to 20.0 GPa. The calculated structure factors are consistent with the recent experimental results, and the partial structure parameters show that the structure of l-GaSb under pressure contracts nonuniformly. In the whole calculated pressure region, the contraction of l-GaSb can be divided into three substages: 1.8-5.4, 5.4-10.0, and 10.0-20.0 GPa. It is further confirmed by analyzing the bond-angle distributions of Ga-Ga-Ga and Sb-Sb-Sb that the rearrangement of Sb atoms under pressure plays a crucial role in the structure change of l-GaSb.

Journal Article↗

PubMed2006

Molecular dynamics simulations of water droplets on polymer surfaces.

Molecular dynamics simulations were used to study the wetting of polymer surfaces with water. Contact angles of water droplets on crystalline and two amorphous polyethylene (PE) and poly(vinyl chloride) (PVC) surfaces were extracted from atomistic simulations. Crystalline surfaces were produced by duplicating the unit cell of an experimental crystal structure, and amorphous surfaces by pressing the bulk polymer step by step at elevated temperature between two repulsive grid surfaces to a target density. Different-sized water droplets on the crystalline PE surface revealed a slightly positive line tension on the order of 10(-12)-10(-11) N, whereas droplets on crystalline PVC did not yield a definite line tension. Microscopic contact angles produced by the simple point charge (SPC) water model were mostly a few degrees smaller than those produced by the extended SPC model, which, as the model with lowest bulk energy, presents an upper boundary for contact angles. The macroscopic contact angle for the SPC model was 94 degrees on crystalline PVC and 113 degrees on crystalline PE. Amorphicity of the surface increased the water contact angle on PE but decreased it on PVC, for both water models. If the simulated contact angles on crystalline and amorphous surfaces are combined in proportion to the crystallinity of the polymer in question, simulated values in relatively good agreement with measured values are obtained.

Journal Article↗

PubMed2002

Molecular dynamics simulations of vibrated granular gases.

We present molecular dynamics simulations of monodisperse or bidisperse inelastic granular gases driven by vibrating walls, in two dimensions (without gravity). Because of the energy injection at the boundaries, a situation often met experimentally, density and temperature fields display heterogeneous profiles in the direction perpendicular to the walls. A general equation of state for an arbitrary mixture of fluidized inelastic hard spheres is derived and successfully tested against numerical data. Single-particle velocity distribution functions with non-Gaussian features are also obtained, and the influence of various parameters (inelasticity coefficients, density, etc.) are analyzed. The validity of a recently proposed random restitution coefficient model is assessed through the study of projected collisions onto the direction perpendicular to that of energy injection. For the binary mixture, the nonequipartition of translational kinetic energy is studied and compared both to experimental data and to the case of homogeneous energy injection ("stochastic thermostat"). The rescaled velocity distribution functions are found to be very similar for both species.

Journal Article↗

PubMed2006

Molecular dynamics simulation of electro-osmotic flows in rough wall nanochannels.

We performed equilibrium and nonequilibrium molecular dynamics simulation to study electro-osmotic flows inside charged nanochannels with different types of surface roughness. We modeled surface roughness as a sequence of two-dimensional subnanoscale grooves and ridges (step function-type roughness) along the flow direction. The amplitude, spatial period, and symmetry of surface roughness were varied. The amplitude of surface roughness was on the order of the Debye length. The walls have uniform negative charges at the interface with fluids. We included only positive ions (counterions) for simplicity of computation. For the smooth wall, we compared our molecular dynamics simulation results to the well-known Poisson-Boltzmann theory. The density profiles of water molecules showed "layering" near the wall. For the rough walls, the density profiles measured from the wall are similar to those for the smooth wall except near where the steps are located. Because of the layering of water molecules and the finite size effect of ions and the walls, the ionic distribution departs from the Boltzmann distribution. To further understand the structure of water molecules and ions, we computed the polarization density. Near the wall, its z component dominates the other components, indicating the preferred orientation ("ordering") of water molecules. Especially, inside the groove for the rough walls, its maximum is 10% higher (stronger ordering) than for the smooth wall. The dielectric constant, computed with a Clausius-Mosotti-type equation, confirmed the ordering near the wall and the enhanced ordering inside the groove. The residence time and the diffusion coefficient, computed using the velocity autocorrelation function, showed that the diffusion of water and ions along the direction normal to the wall is significantly reduced near the wall and further decreases inside the groove. Along the flow direction, the diffusion of water and ions inside the groove is significantly lowered while it is similar to the bulk value elsewhere. We performed nonequilibrium molecular dynamics simulation to compute electro-osmotic velocities and flow rates. The velocity profiles correspond to those for overlapped electric double layers. For the rough walls, velocity inside the groove is close to zero, meaning that the channel height is effectively reduced. The flow rate was found to decrease as the period of surface roughness decreases or the amplitude of surface roughness increases. We defined the zeta potential as the electrostatic potential at the location of a slip plane. We computed the electrostatic potential with the ionic distribution and the dielectric constant both from our molecular dynamics simulation. We estimated the slip plane from the velocity profile. The zeta potential showed the same trend as the flow rate: it decreases with an increasing amplitude and a decreasing period of surface roughness.

Journal Article↗

PubMed1997

Model-free methods of analyzing domain motions in proteins from simulation: a comparison of normal mode analysis and molecular dynamics simulation of lysozyme.

Model-free methods are introduced to determine quantities pertaining to protein domain motions from normal mode analyses and molecular dynamics simulations. For the normal mode analysis, the methods are based on the assumption that in low frequency modes, domain motions can be well approximated by modes of motion external to the domains. To analyze the molecular dynamics trajectory, a principal component analysis tailored specifically to analyze interdomain motions is applied. A method based on the curl of the atomic displacements is described, which yields a sharp discrimination of domains, and which defines a unique interdomain screw-axis. Hinge axes are defined and classified as twist or closure axes depending on their direction. The methods have been tested on lysozyme. A remarkable correspondence was found between the first normal mode axis and the first principal mode axis, with both axes passing within 3 A of the alpha-carbon atoms of residues 2, 39, and 56 of human lysozyme, and near the interdomain helix. The axes of the first modes are overwhelmingly closure axes. A lesser degree of correspondence is found for the second modes, but in both cases they are more twist axes than closure axes. Both analyses reveal that the interdomain connections allow only these two degrees of freedom, one more than provided by a pure mechanical hinge.

Animals↗

PubMed2006

Balancing an accurate representation of the molecular surface in generalized born formalisms with integrator stability in molecular dynamics simulations.

Different integrator time steps in NVT and NVE simulations of protein and nucleic acid systems are tested with the GBMV (Generalized Born using Molecular Volume) and GBSW (Generalized Born with simple SWitching) methods. The simulation stability and energy conservation is investigated in relation to the agreement with the Poisson theory. It is found that very close agreement between generalized Born methods and the Poisson theory based on the commonly used sharp molecular surface definition results in energy drift and simulation artifacts in molecular dynamics simulation protocols with standard 2-fs time steps. New parameters are proposed for the GBMV method, which maintains very good agreement with the Poisson theory while providing energy conservation and stable simulations at time steps of 1 to 1.5 fs.

Computer Simulation↗

PubMed2005

Supercooled water in PVA matrixes. II. A molecular dynamics simulation study and comparison with QENS results.

Molecular dynamics (MD) simulations were carried out to elucidate the dynamic behavior of water confined in poly(vinyl alcohol), PVA, hydrogels. Model topology is supported by experimental network parameters, and simulation results are compared to an incoherent quasielastic neutron scattering (QENS) investigation carried out on PVA hydrogels. From the QENS dynamic scattering law (part I), a random jump model was adopted for the description of water diffusion to extract a microscopic diffusion coefficient and a residence time between two "jumps". In the present work, consistently with this framework, water diffusion parameters as diffusion coefficients and residence times have been evaluated using the mean square displacement of water in a time window of 10 ps and the time autocorrelation function of water hydrogen bonds. The calculated parameters are in good agreement with the experimental ones, giving confidence to this approach. Further developments are in progress to take into account a more realistic description of hydrogel structure in the molecular dynamics simulations.

Journal Article↗

PubMed1994

Molecular dynamics simulation of the gramicidin channel in a phospholipid bilayer.

A molecular dynamics simulation of the gramicidin A channel in an explicit dimyristoyl phosphatidylcholine bilayer was generated to study the details of lipid-protein interactions at the microscopic level. Solid-state NMR properties of the channel averaged over the 500-psec trajectory are in excellent agreement with available experimental data. In contrast with the assumptions of macroscopic models, the membrane/solution interface region is found to be at least 12 A thick. The tryptophan side chains, located within the interface, are found to form hydrogen bonds with the ester carbonyl groups of the lipids and with water, suggesting their important contribution to the stability of membrane proteins. Individual lipid-protein interactions are seen to vary from near 0 to -50 kcal/mol. The most strongly interacting conformations are short-lived and have a nearly equal contribution from both van der Waals and electrostatic energies. This approach for performing molecular dynamics simulations of membrane proteins in explicit phospholipid bilayers should help in studying the structure, dynamics, and energetics of lipid-protein interactions.

Gramicidin↗

PubMed2004

Methods for molecular dynamics simulations of protein folding/unfolding in solution.

All atom molecular dynamics simulations have become a standard method for mapping equilibrium protein dynamics and non-equilibrium events like folding and unfolding. Here, we present detailed methods for performing such simulations. Generic protocols for minimization, solvation, simulation, and analysis derived from previous studies are also presented. As a measure of validation, our water model is compared with experiment. An example of current applications of these methods, simulations of the ultrafast folding protein Engrailed Homeodomain are presented including the experimental evidence used to verify their results. Ultrafast folders are an invaluable tool for studying protein behavior as folding and unfolding events measured by experiment occur on timescales accessible with the high-resolution molecular dynamics methods we describe. Finally, to demonstrate the prospect of these methods for folding proteins, a temperature quench simulation of a thermal unfolding intermediate of the Engrailed Homeodomain is described.

Computer Simulation↗

PubMed2000

Changes at the floor of the peptide-binding groove induce a strong preference for proline at position 3 of the bound peptide: molecular dynamics simulations of HLA-A*0217.

We report on molecular dynamics simulations of major histocompatibility complex (MHC)-peptide complexes. Class I MHC molecules play an important role in cellular immunity by presenting antigenic peptides to cytotoxic T cells. Pockets in the peptide-binding groove of MHC molecules accommodate anchor side chains of the bound peptide. Amino acid substitutions in MHC affect differences in the peptide-anchor motifs. HLA-A*0217, human MHC class I molecule, differs from HLA-A*0201 only by three amino acid residues substitutions (positions 95, 97, and 99) at the floor of the peptide-binding groove. A*0217 showed a strong preference for Pro at position 3 (p3) and accepted Phe at p9 of its peptide ligands, but these preferences have not been found in other HLA-A2 ligands. To reveal the structural mechanism of these observations, the A*0217-peptide complexes were simulated by 1000 ps molecular dynamics at 300 K with explicit solvent molecules and compared with those of the A*0201-peptide complexes. We examined the distances between the anchor side chain of the bound peptide and the pocket, and the rms fluctuations of the bound peptides and the HLA molecules. On the basis of the results from our simulations, we propose that Pro at p3 serves as an optimum residue to lock the dominant anchor residue (p9) tightly into pocket F and to hold the peptide in the binding groove, rather than a secondary anchor residue fitting optimally the complementary pocket. We also found that Phe at p9 is used to occupy the space created by replacements of three amino acid residues at the floor within the groove. These findings would provide a novel understanding in the peptide-binding motifs of class I MHC molecules.

Amino Acid Motifs↗

PubMed1999

Defining the transmembrane helix of M2 protein from influenza A by molecular dynamics simulations in a lipid bilayer.

Integral membrane proteins containing at least one transmembrane (TM) alpha-helix are believed to account for between 20% and 30% of most genomes. There are several algorithms that accurately predict the number and position of TM helices within a membrane protein sequence. However, these methods tend to disagree over the beginning and end residues of TM helices, posing problems for subsequent modeling and simulation studies. Molecular dynamics (MD) simulations in an explicit lipid and water environment are used to help define the TM helix of the M2 protein from influenza A virus. Based on a comparison of the results of five different secondary structure prediction algorithms, three different helix lengths (an 18mer, a 26mer, and a 34mer) were simulated. Each simulation system contained 127 POPC molecules plus approximately 3500-4700 waters, giving a total of approximately 18,000-21,000 atoms. Two simulations, each of 2 ns duration, were run for the 18mer and 26mer, and five separate simulations were run for the 34mer, using different starting models generated by restrained in vacuo MD simulations. The total simulation time amounted to 11 ns. Analysis of the time-dependent secondary structure of the TM segments was used to define the regions that adopted a stable alpha-helical conformation throughout the simulation. This analysis indicates a core TM region of approximately 20 residues (from residue 22 to residue 43) that remained in an alpha-helical conformation. Analysis of atomic density profiles suggested that the 18mer helix revealed a local perturbation of the lipid bilayer. Polar side chains on either side of this region form relatively long-lived H-bonds to lipid headgroups and water molecules.

Amino Acid Sequence↗

PubMed1991

Mass-weighted molecular dynamics simulation of the protein-ligand complex of rhizopuspepsin and inhibitor.

The mass-weighted molecular dynamics simulation method was developed previously for sampling the multidimensional conformational space of linear and cyclic polypeptides and studying their conformational flexibility. Herein results from molecular dynamics simulations of the protein-ligand complex of the aspartyl protease rhizopuspepsin and a polypeptide inhibitor are reported. The dihedral conformational space sampling for the linear peptide inhibitor in situ was found to be increased in the mass-weighted simulation as in other molecular systems previously studied. More significantly, the physical space of the enzyme binding pocket was also sampled efficiently in the simulations and multiple binding sites were identified for the inhibitor. These results suggest that it may be possible now to study, by computer simulations, the putative initial enzyme-inhibitor complex suggested experimentally from the time-dependent kinetics of enzyme inhibition by slow-binding inhibitors (Morrison, J. F., and C. T. Walsh. 1988. Adv. Enzymol. 61:201), and/or conformational substates in protein-ligand complexes suggested in the study of reassociation dynamics of myoglobin and carbon monoxide following photolysis (Austin, R. H., K. W. Beeson, L. Eisenstein, H. Frauenfelder, and I. C. Gunsalus. 1975. Biochemistry. 14:5355). Moreover, the intermediate binding steps and the molecular flexibility of the inhibitor shown in the MWMD simulation may have crucial roles in the ligand binding process.

Amino Acid Sequence↗

PubMed2002

Molecular dynamics simulation of brittle fracture in silicon.

Brittle fracture in silicon is simulated with molecular dynamics utilizing a modified embedded atom method potential. The simulations produce propagating crack speeds that are in agreement with previous experimental results over a large range of fracture energy. The dynamic fracture toughness is found to be equal to the energy consumed by creating surfaces and lattice defects in agreement with theoretical predictions. The dynamic fracture toughness is approximately 1/3 of the static strain energy release rate, which results in a limiting crack speed of 2/3 of the Rayleigh wave speed.

Journal Article↗