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,621 records · Page 90Linked to original sources
PubMed2005

Application of molecular dynamics simulations to spin-labeled oligonucleotides.

The EPR study of spin labeled macromolecules has provided insight into structural and dynamical properties of DNA, proteins, and related systems. While spin labeling has been useful, it is experimentally difficult to determine if the spin label significantly alters the conformation of the macromolecule to which it is attached. Molecular modeling has proven to be a powerful tool for studying structure and dynamics of biologically important molecules. Here, we have conducted molecular dynamics (MD) studies of spin labeled oligonucleotides (ONs) bearing a five (5sp) or six (6sp) membered ring nitroxide, and the corresponding unmodified ON using the suite of programs contained in Amber 5.0 with the Cornell et al. 94 force field (Cornell, W. D., Cieplak, P., Bayly, C. I., Gould, I. R., Merz, Jr., K. M. Ferguson, D. M., Spellmeyer, D. C., Fox, T., Caldwell, J. W., and Kollman, P. A. A Second Generation Force Field for the Simulation of Proteins and Nucleic Acids. J. Am. Chem. Soc. 117, 5179-5197 (1995)). Quantum mechanical calculations employing the B3LYP method with the standard 6-31G* basis set using Gaussian98 were performed and, together with available crystallographic data for analogous nitroxides, new parameters for the nitrogen, oxygen, nitroxide alpha-carbon, and sp-hybridized carbon atoms have been developed suitable for the Cornell et al. 94 force field. MD simulations on the double-stranded (ds) spin labeled ONs, along with the corresponding unmodified analogues, have been studied over the course of 4 ns and conformational properties of all ONs are described based on the analysis of the trajectories. The spin labels were found to alter the global conformation of the ONs to which they were attached to accommodate the spin labels. The major changes include widening the major groove, decreasing helical twist, and more negative X-displacement of the base pairs. The magnitude of the effect was dependent on the specific structure of the spin label. Average and 'most representative' structures derived from the molecular dynamics simulations correlate with the experimental data on the spin labeled ONs.

Computer Simulation↗
PubMed1989

Molecular dynamics simulation as a complement to diffraction in the study of disorder in crystals.

It is not uncommon that a crystal structure contains a minority component which does not share the overall space-group of the remainder of the structure, but has a lower (or no) effective space-group symmetry. A conventional diffraction treatment will inevitably be made within the space-group of the majority of the system, and will therefore result in an average structure also with this symmetry. We here indicate the value of molecular dynamics (MD) simulation within a supercell of the structure as a complementary tool in providing information concerning the local structure for the minority component. Such information, though illustrative, is often the most interesting and the physically most significant aspect of the structure. It can also be crucial information for the testing of new theory, and in the design of new materials. The approach is illustrated for the case of the Na+ ion distribution in Na+ beta"-alumina, and for the experimentally more inaccessible mixed-ion system Na+/Ba2+ beta"-alumina.

Aluminum↗
PubMed2004

Probing flexibility and "induced-fit" phenomena in aldose reductase by comparative crystal structure analysis and molecular dynamics simulations.

Aldose reductase is a promising target for the treatment of diabetic complications, and as such, has become the focus of various drug design projects. As revealed by a survey of available crystal structures, the protein shows pronounced induced-fit effects upon ligand binding. Although helping to explain the enzyme's substrate promiscuity, phenomena of this kind are still responsible for significant complications in structure-based design efforts directed to aldose reductase. Accordingly, a deeper understanding of the principles governing conformational alterations in this enzyme would be of utmost practical importance. As a first step in addressing this issue, molecular dynamics (MD) simulations have been carried out. The ultrahigh resolution crystal structure of aldose reductase complexed with inhibitor IDD594 served as ideal starting point for a set of different simulations of nanosecond time scale: the native complexed state with bound inhibitor, the uncomplexed state (after removal of the inhibitor) at standard temperature, and the uncomplexed state at elevated temperature. The reference simulation of the complex exhibits extraordinary stability of the overall fold, whereas two distinct conformational substates are found for the binding-site region. In contrast, already at standard temperature pronounced changes are observed in the binding region during the simulation of the uncomplexed state. Leu300, for example, closes the access to the pocket opened by IDD594. On the other hand, conformations around the catalytic site are highly conserved, with the His110-Tyr48-NADP+ orientation being stabilized by a water molecule. Detailed analysis of the trajectories allows to reveal a set of distinct conformational substates that may prove useful as alternative structural templates in virtual screening for new aldose reductase inhibitors.

Acetates↗
PubMed2006

Mixed atomistic and coarse-grained molecular dynamics: simulation of a membrane-bound ion channel.

The recently developed multiscale coarse-graining (MS-CG) method (Izvekov, S.; Voth, G. A. J. Phys. Chem. B 2005, 109, 2469; J. Chem. Phys. 2005, 123, 134105) is used to build a mixed all-atom and coarse-grained (AA-CG) model of the gramicidin A (gA) ion channel embedded in a dimyristoylphosphatidylcholine (DMPC) lipid bilayer and water environment. In this model, the gA peptide was described in full atomistic detail, while the lipid and water molecules were described using coarse-grained representations. The atom-CG and CG-CG interactions in the mixed AA-CG model were determined using the MS-CG method. Molecular dynamics (MD) simulations were performed using the resulting AA-CG model. The results from simulations of the AA-CG model compare very favorably to those from all-atom MD simulations of the entire system. Since the MS-CG method employs a general and systematic approach to obtain effective interactions from the underlying all-atom models, the present approach to rigorously develop mixed AA-CG models has the potential to be extended to many other systems.

Computer Simulation↗
PubMed1996

Molecular dynamics simulation of hydrated phospholipid bilayers.

Understanding of microscopic behaviour of biological membrane is crucial for designing of molecules to control transport properties of the membranes. Phospholipid-water forms a good model system to study ligand induced structural and dynamical changes in membrane. The review has its main focus on molecular dynamics (MD) simulation of phospholipid bilayers. A brief summary of the current status of structure of phospholipid membranes based on different physico-chemical measurements is given. We discuss here mainly results of MD simulations in the recent years on hydrated phospholipid bilayers and their interaction with ligands. Simulation parameters as: choice of initial system, force fields, protocols for simulation are compared. Main results on: order parameter, head group and chain conformation, water penetration profile, chain tilts, pair-correlation function between atoms of lipid and water, diffusion of ions and ligands are discussed. The review gives application and limitation of MD method for studying lipid water system.

Computer Simulation↗
PubMed2006

Synthetic polymers and biomembranes. How do they interact? Atomistic molecular dynamics simulation study of PEO in contact with a DMPC lipid bilayer.

The understanding of interactions of poly(ethylene glycol) (PEG) or poly(ethylene oxide) (PEO) with biological interfaces has important technological application in industry and in medicine. In this paper, structural and dynamical properties of PEO at the dimyristoylphospatidylcholine (DMPC) bilayer/water interface have been investigated by molecular dynamics (MD) and steered molecular dynamics (SMD) simulations. The structural properties of a PEO chain in bulk water, at the water/vacuum interface, and in the presence of the membrane were compared with available experimental data. The presence of a barrier for the PEO penetration into the DMPC bilayer has been found. A qualitative estimation of the barrier provided a value equal to approximately 19 kJ/mol, that is, 7 times the value of kT at 310 K.

Dimyristoylphosphatidylcholine↗
PubMed2004

Monte Carlo versus molecular dynamics simulations in heterogeneous systems: an application to the n-pentane liquid-vapor interface.

The Monte Carlo (MC) and molecular dynamics (MD) methodologies are now well established for computing equilibrium properties in homogeneous fluids. This is not yet the case for the direct simulation of two-phase systems, which exhibit nonuniformity of the density distribution across the interface. We have performed direct MC and MD simulations of the liquid-gas interface of n-pentane using a standard force-field model. We obtained density and pressure components profiles along the direction normal to the interface that can be very different, depending on the truncation and long range correction strategies. We discuss the influence on predicted properties of different potential truncation schemes implemented in both MC and MD simulations. We show that the MD and MC profiles can be made in agreement by using a Lennard-Jones potential truncated via a polynomial function that makes the first and second derivatives of the potential continuous at the cutoff distance. In this case however, the predicted thermodynamic properties (phase envelope, surface tension) deviate from experiments, because of the changes made in the potential. A further readjustment of the potential parameters is needed if one wants to use this method. We conclude that a straightforward use of bulk phase force fields in MD simulations may lead to some physical inconsistencies when computing interfacial properties.

Journal Article↗
PubMed1994

X-ray crystal structure and molecular dynamics simulation of bovine pancreas phospholipase A2-n-dodecylphosphorylcholine complex.

The crystal structure of n-dodecylphosphorylcholine (n-C12PC)-bovine pancreas phospholipase A2 (PLA2) complex provided the following structural characteristics: (1) the dodecyl chain of n-C12PC was located at the PLA2 N-terminal helical region by hydrophobic interactions, which corresponds to the binding pocket of 2-acyl fatty acid chain (beta-chain) of the substrate phospholipid, (2) the region from Lys-53 to Lys-56 creates a choline-receiving pocket of n-C12PC and (3) the N-terminal group of Ala-1 shifts significantly toward the Tyr-52 OH group by the binding of the n-C12PC inhibitor. Since the accuracy of the X-ray analysis (R = 0.275 at 2.3 A resolution) was insufficient to establish these important X-ray insights, the complex structure was further investigated through the molecular dynamics (MD) simulation, assuming a system in aqueous solution at 310K. The MD simulation covering 176 ps showed that the structural characteristics observed by X-ray analysis are intrinsic and also stable in the dynamic state. Furthermore, the MD simulation made clear that the PLA2 binding pocket is large enough to permit the conformational fluctuation of the n-C12PC hydrocarbon chain.

Animals↗
PubMed2005

Multiple loop conformations of peptides predicted by molecular dynamics simulations are compatible with nuclear magnetic resonance.

The affinity and selectivity of protein-protein interactions can be fine-tuned by varying the size, flexibility, and amino acid composition of involved surface loops. As a model for such surface loops, we study the conformational landscape of an octapeptide, whose flexibility is chemically steered by a covalent ring closure integrating an azobenzene dye into and by a disulfide bridge additionally constraining the peptide backbone. Because the covalently integrated azobenzene dyes can be switched by light between a bent cis state and an elongated trans state, six cyclic peptide models of strongly different flexibilities are obtained. The conformational states of these peptide models are sampled by NMR and by unconstrained molecular dynamics (MD) simulations. Prototypical conformations and the free-energy landscapes in the high-dimensional space spanned by the phi/psi angles at the peptide backbone are obtained by clustering techniques from the MD trajectories. Multiple open-loop conformations are shown to be predicted by MD particularly in the very flexible cases and are shown to comply with the NMR data despite the fact that such open-loop conformations are missing in the refined NMR structures.

Aminobenzoates↗
PubMed2006

Atomistic molecular dynamics simulation of the temperature and pressure dependences of local and terminal relaxations in cis-1,4-polybutadiene.

The dynamics of cis-1,4-polybutadiene (cis-1,4-PB) over a wide range of temperature and pressure conditions is explored by conducting atomistic molecular dynamics (MD) simulations with a united atom model on a 32-chain C128 cis-1,4-PB system. The local or segmental dynamics is analyzed in terms of the dipole moment time autocorrelation function (DACF) of the simulated polymer and its temperature and pressure variations, for temperatures as low as 195 K and pressures as high as 3 kbars. By Fourier transforming the DACF, the dielectric spectrum, epsilon* = epsilon' + i epsilon" = epsilon*omega, is computed and the normalized epsilon"/epsilon(max)" vs omega/omega(max) plot is analyzed on the basis of the time-temperature and time-pressure superposition principles. The relative contribution of thermal energy and volume to the segmental and chain relaxation processes are also calculated and evaluated in terms of the ratio of the activation energy at constant volume to the activation energy at constant pressure, Q(V)/Q(P). Additional results for the temperature and pressure dependences of the Rouse times describing terminal relaxation in the two polymers show that, in the regime of the temperature and pressure conditions covered here, segmental and chain relaxations are influenced similarly by the pressure and temperature variations. This is in contrast to what is measured experimentally [see, e.g., G. Floudas and T. Reisinger, J. Chem. Phys. 111, 5201 (1999); C. M. Roland et al.,J. Polym. Sci. Part B, 41, 3047 (2003)] for other, chemically more complex polymers that pressure has a stronger influence on the dynamics of segmental mode than on the dynamics of the longest normal mode, at least for the regime of temperature and pressure conditions covered in the present MD simulations.

Journal Article↗
PubMed2006

Speciation of the curium(III) ion in aqueous solution: a combined study by quantum chemistry and molecular dynamics simulation.

The structures of aquo complexes of the curium(III) ion have been systematically studied using quantum chemical and molecular dynamics (MD) methods. The first hydration shell of the Cm3+ ion has been calculated using density functional theory (DFT), with and without inclusion of the conductor-like polarizable continuum medium (CPCM) model of solvation. The calculated results indicate that the primary hydration number of Cm3+ is nine, with a Cm-O bond distance of 2.47-2.48 A. The calculated bond distances and the hydration number are in excellent agreement with available experimental data. The inclusion of a complete second hydration shell of Cm3+ has been investigated using both DFT and MD methods. The presence of the second hydration shell has significant effects on the primary coordination sphere, suggesting that the explicit inclusion of second-shell effects is important for understanding the nature of the first shell. The calculated results indicate that 21 water molecules can be coordinated in the second hydration shell of the Cm3+ ion. MD simulations within the hydrated-ion model suggest that the second-shell water molecules exchange with the bulk solvent with a lifetime of 161 ps.

Journal Article↗
PubMed2006

Catalytic mechanism of cyclophilin as observed in molecular dynamics simulations: pathway prediction and reconciliation of X-ray crystallographic and NMR solution data.

Cyclophilins are proteins that catalyze X-proline cis-trans interconversion, where X represents any amino acid. Its mechanism of action has been investigated over the past years but still generates discussion, especially because until recently structures of the ligand in the cis and trans conformations for the same system were lacking. X-ray crystallographic structures for the complex cyclophilin A and HIV-1 capsid mutants with ligands in the cis and trans conformations suggest a mechanism where the N-terminal portion of the ligand rotates during the cis-trans isomerization. However, a few years before, a C-terminal rotating ligand was proposed to explain NMR solution data. In the present study we use molecular dynamics (MD) simulations to generate a trans structure starting from the cis structure. From simulations starting from the cis and trans structures obtained through the rotational pathways, the seeming contradiction between the two sets of experimental data could be resolved. The simulated N-terminal rotated trans structure shows good agreement with the equivalent crystal structure and, moreover, is consistent with the NMR data. These results illustrate the use of MD simulation at atomic resolution to model structural transitions and to interpret experimental data.

Catalysis↗
PubMed2006

Molecular dynamics simulations of arachidonic acid complexes with COX-1 and COX-2: insights into equilibrium behavior.

The cyclooxygenase (COX) enzymes are responsible for the committed step in prostaglandin biosynthesis, the generation of prostaglandin H(2). As a result, these enzymes are pharmacologically important targets for nonsteroidal antiinflammatory drugs, such as aspirin and newer COX-2 selective inhibitors. The cyclooxygenases are functional homodimers, and each subunit contains both a cyclooxygenase and a peroxidase active site. These enzymes are quite interesting mechanistically, as the conversion of arachidonic acid to prostaglandin H(2) requires two oxygenation and two cyclization reactions, resulting in the formation of five new chiral centers with nearly absolute regio- and stereochemical fidelity. We have used molecular dynamics (MD) simulations to investigate the equilibrium behavior of both COX-1 and COX-2 enzyme isoforms with bound arachidonate. These simulations were compared with reference simulations of arachidonate in solution to explore the effect of enzyme on substrate conformation and positioning in the active site. The simulations suggest that the substrate has greater conformational freedom in the COX-2 active site, consistent with the larger COX-2 active site volume observed in X-ray crystal structures. The simulations reveal different conformational behavior for arachidonate in each subunit over the course of extended equilibrium MD simulations. The simulations also provide detailed information for several protein channels that might be important for oxygen and water transport to or from active sites or for intermediate trafficking between the cyclooxygenase and peroxidase active sites. The detailed comparisons for COX-1 versus COX-2 active site structural fluctuations may also provide useful information for design of new isozyme-selective inhibitors.

Arachidonic Acid↗
PubMed2005

Molecular-dynamics simulations of alkaline-earth metal cations in water by atom-bond electronegativity equalization method fused into molecular mechanics.

Intermolecular potential for alkaline-earth metal (Be(2+), Mg(2+), and Ca(2+)) cations in water has been derived using the atom-bond electronegativity equalization method fused into molecular mechanics (ABEEM/MM), and it is consistent with what was previously applied to the hydration study of the monovalent cations. Parameters for the effective interaction between a cation and a water molecule were determined, reproducing the ab initio results. The static, dynamic, and thermodynamic properties of Be(2+)(aq), Mg(2+)(aq), and Ca(2+)(aq) were studied using these potential parameters. Be(2+) requires a more complicated form of the potential function than Mg(2+) and Ca(2+) in order to obtain better fits. Strong influences of the twofold charged cations on the structures of the hydration shells and some other properties of aqueous ionic solutions are discussed and compared with the results of a previous study of monovalent cations in water. At the same time, comparative study of the hydration properties of each cation is also discussed. This work demonstrates that ABEEM/MM provides a useful tool in the exploration of the hydration of double-charged cations in water.

Algorithms↗
PubMed2000

Molecular dynamics simulations of bovine cathepsin B and its complex with CA074.

To promote our better understanding of the dynamic stability of the bovine cathepsin B structure, which is characterized by an extra disulfide bond at Cys148-Cys252 from the other species, and of the binding stability of CA074 (a cathepsin B-specific inhibitor), molecular dynamics (MD) simulations were performed for the enzyme and its CA074 complex, assuming a system in aqueous solution at 300 K. The MD simulation covering 400 ps indicated that the existence of a Cys148-Cys252 disulfide bond increases the conformational flexibility of the occluding loop, although the conformational stability of the overall structure is little affected. The structural characteristics of the complex elucidated by X-ray analysis were suggested to be also intrinsic and stable in the dynamic state; the hydrogen bonding/electrostatic interactions between the main and side chains of CA074 and the Sn and Sn' subsites of the enzyme were maintained throughout the MD simulation. Furthermore, the simulation made clear that the binding of CA074 significantly restricted the conformational flexibility of the substrate binding region, especially the occluding loop, of cathepsin B. Statistical analyses during the simulation suggest that the selectivity of CA074 for cathepsin B stems from the tight P1'-S1' and P2'-S2' interactions, assisted in particular by double hydrogen bonds between the carboxyl two oxygens of the CA074 C-terminus and the imidazole NH groups of His110 and His111 residues.

Animals↗
PubMed1994

Molecular dynamics simulations of oligosaccharides and their conformation in the crystal structure of lectin-carbohydrate complex: importance of the torsion angle psi for the orientation of alpha 1,6-arm.

The conformation of the heptasaccharide Man-alpha 1,6-(Man-alpha 1,3)(Xyl-beta 1,2)-Man-beta 1,4-GlcNAc2-beta 1,4-(L- Fuc-alpha 1,3)-GlcNAc1, the carbohydrate moiety of Erythrina corallodendron lectin (EcorL), the hexasaccharide Man-alpha 1,6-(Man-alpha 1,3) (GlcNAc-beta 1,4)-Man-beta 1,4-GlcNAc-beta 1,4-GlcNAc and their disaccharide fragments have been studied by molecular dynamics (MD) simulations for 1000 ps with different initial conformations. In the isolated heptasaccharide, the most frequently accessed conformation during MD has a psi value of 180 degrees around Man-alpha 1,6-Man linkage. This conformation is stabilized by the formation of a hydrogen bond between the carbonyl oxygen of GlcNAc2 with the O3/O4 hydroxyls of the alpha 1,6-linked mannose residue. The conformation of the heptasaccharide found in the crystal structure of the EcorL-lactose complex (Shaanan et al., Science, 254, 862, 1991), that has a psi value of approximately 76 degrees around Man-alpha 1,6-Man linkage, is accessed, although less frequently, during MD of the isolated oligosaccharide. The phi, psi, chi = 58 degrees, -134 degrees, -60 degrees conformation around Man-alpha 1,6-Man fragment observed in the crystal structure of the Lathyrus ochrus lectin complexed with a biantennary octasaccharide (Table I in Homans, S.W., Glycobiology, 3, 551, 1993) has also been accessed in the present MD simulations. These psi values for the alpha 1,6-linkage, which are observed in the protein-carbohydrate crystal structures and are accessed in the MD simulations, though occasionally, have not been predicted from NMR studies. Furthermore, these different values of psi lead to significantly different orientations of the alpha 1,6-arm for the same value of chi. This contrasts with the earlier predictions that only different values of chi can bring about significant changes in the orientation of the alpha 1,6-arm. The MD simulations also show that the effects of bisecting GlcNAc or beta 1,2-xylose are very similar on the alpha 1,3-arm and slightly different on the alpha 1,6-arm.

Carbohydrate Conformation↗
PubMed2003

Molecular dynamics simulation of dimeric and monomeric forms of human prion protein: insight into dynamics and properties.

A central theme in prion protein research is the detection of the process that underlies the conformational transition from the normal cellular prion form (PrP(C)) to its pathogenic isoform (PrP(Sc)). Although the three-dimensional structures of monomeric and dimeric human prion protein (HuPrP) have been revealed by NMR spectroscopy and x-ray crystallography, the process underlying the conformational change from PrP(C) to PrP(Sc) and the dynamics and functions of PrP(C) remain unknown. The dimeric form is thought to play an important role in the conformational transition. In this study, we performed molecular dynamics (MD) simulations on monomeric and dimeric HuPrP at 300 K and 500 K for 10 ns to investigate the differences in the properties of the monomer and the dimer from the perspective of dynamic and structural behaviors. Simulations were also undertaken with Asp178Asn and acidic pH, which is known as a disease-associated factor. Our results indicate that the dynamics of the dimer and monomer were similar (e.g., denaturation of helices and elongation of the beta-sheet). However, additional secondary structure elements formed in the dimer might result in showing the differences in dynamics and properties between the monomer and dimer (e.g., the greater retention of dimeric than monomeric tertiary structure).

Binding Sites↗
PubMed2006

Molecular dynamics simulation studies of a protein-RNA complex with a selectively modified binding interface.

The RNA recognition motif (RRM) is one of the most common RNA binding domains. We have investigated the contribution of three highly conserved aromatic amino acids to RNA binding by the N-terminal RRM of the U1A protein. Recently, we synthesized a modified base (A-4CPh) in which a phenyl group is tethered to adenine using a linker of 4 methylene groups. The substitution of this base for adenine in the target RNA selectively stabilizes the complex formed with a U1A protein in which one of the conserved aromatic amino acids is replaced with Ala (Phe56Ala). In this article, we report molecular dynamics (MD) simulations that probe the structural consequences of the substitution of A-4CPh for adenine in the wild type and Phe56Ala U1A-RNA complexes and in the free RNA. The simulations suggest that A-4CPh stabilizes the complex formed with Phe56Ala by adopting a folded conformation in which the tethered phenyl group fills the site occupied by the phenyl group of Phe56 in the wild-type complex. In contrast, an extended conformation of A-4CPh is predicted to be most stable in the complex formed with the wild-type protein. The calculations indicate A-4CPh is in an extended conformation in the free RNA. Therefore, preorganizing the structure of the phenyl-tethered base for binding may improve both the affinity and specificity of the RNA containing A-4CPh for the Phe56Ala U1A protein. Taken together, the previous experimental work and the calculations reported here suggest a general design strategy for altering RRM-RNA complex stability.

Computer Simulation↗