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.

SourceResearch connectionFrom yearThrough year

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

At least 775 records · Page 43Linked to original sources

PubMed1999

Molecular dynamics simulations of the complex between human U1A protein and hairpin II of U1 small nuclear RNA and of free RNA in solution.

RNA-protein interactions are essential to a wide range of biological processes. In this paper, a 0.6-ns molecular dynamics simulation of the sequence-specific interaction of human U1A protein with hairpin II of U1 snRNA in solution, together with a 1.2-ns simulation of the free RNA hairpin, is reported. Compared to the findings in the x-ray structure of the complex, most of the interactions remained stable. The nucleotide U8, one of the seven conserved nucleotides AUUGCAC in the loop region, was unusually flexible during the simulation, leading to a loss of direct contacts with the protein, in contrast to the situation in the x-ray structure. Instead the sugar-phosphate backbone of nucleotide C15 was found to form several interactions with the protein. Compared to the NMR structure of U1A protein complexed with the 3'-untranslated region of its own pre-mRNA, the protein core kept the same conformation, and in the two RNA molecules the conserved AUUGCAC of the loop and the closest CG base pair were located in very similar positions and orientations, and underwent very similar interactions with the protein. Therefore, a common sequence-specific interaction mechanism was suggested for the two RNA substrates to bind to the U1A protein. Conformational analysis of the RNA hairpin showed that the conformational changes of the RNA primarily occurred in the loop region, which is just involved in the sites of binding to the protein and in agreement with experimental observation. Both the loop and stem of the RNA became more ordered upon binding to the protein. It was also demonstrated that the molecular dynamics method could be successfully used to simulate the dynamical behavior of a large RNA-protein complex in aqueous solution, thus opening a path for the exploration of the complex biological processes involving RNA at a molecular level.

Amino Acid Sequence↗

PubMed1996

The new program OPAL for molecular dynamics simulations and energy refinements of biological macromolecules.

A new program for molecular dynamics (MD) simulation and energy refinement of biological macromolecules, OPAL, is introduced. Combined with the supporting program TRAJEC for the analysis of MD trajectories, OPAL affords high efficiency and flexibility for work with different force fields, and offers a user-friendly interface and extensive trajectory analysis capabilities. Salient features are computational speeds of up to 1.5 GFlops on vector supercomputers such as the NEC SX-3, ellipsoidal boundaries to reduce the system size for studies in explicit solvents, and natural treatment of the hydrostatic pressure. Practical applications of OPAL are illustrated with MD simulations of pure water, energy minimization of the NMR structure of the mixed disulfide of a mutant E. coli glutaredoxin with glutathione in different solvent models, and MD simulations of a small protein, pheromone Er-2, using either instantaneous or time-averaged NMR restraints, or no restraints.

Computer Simulation↗

PubMed1997

Molecular dynamics simulations of the two disaccharides of hyaluronan in aqueous solution.

Hyaluronan is an unusually stiff polymer when in aqueous solution, which has important consequences for its biological function. Molecular dynamics simulations of hyaluronan disaccharides have been performed, with explicit inclusion of water, to determine the molecular basis of this stiffness, and to investigate the dynamics of the glycosidic linkages. Our simulations reveal that stable sets of hydrogen bonds frequently connect the neighboring residues of hyaluronan. Water caging around the glycosidic linkage was observed to increase the connectivity between sugars, and further constrain them. This, we propose, explains the unusual stiffness of polymeric hyaluronan. It would allow the polysaccharide to maintain local secondary structure, and occupy large solution domains consistent with the visco-elastic nature of hyaluronan. Simulations in water showed no significant changes on inclusion of the exoanomeric effect. This, we deduced, was due to hyaluronan disaccharides ordering first shell water molecules. In some cases these waters were observed to transiently induce conformational change, by breaking intramolecular hydrogen bonds.

Acetylglucosamine↗

PubMed2002

Molecular dynamics simulations of the first steps of the reaction catalyzed by HIV-1 protease.

The mechanism of the first steps of the reaction catalyzed by HIV-1 protease was studied through molecular dynamics simulations. The potential energy surface in the active site was generated using the approximate valence bond method. The approximate valence bond (AVB) method was parameterized based on density functional calculations. The surrounding protein and explicit water environment was modeled with conventional, classical force field. The calculations were performed based on HIV-1 protease complexed with the MVT-101 inhibitor that was modified to a model substrate. The protonation state of the catalytic aspartates was determined theoretically. Possible reaction mechanisms involving the lytic water molecule are accounted for in this study. The modeled steps include the dissociation of the lytic water molecule and proton transfer onto Asp-125, the nucleophilic attack followed by a proton transfer onto peptide nitrogen. The simulations show that in the active site most preferable energetically are structures consisting of ionized or polarized molecular fragments that are not accounted for in conventional molecular dynamics. The mobility of the lytic water molecule, the dynamics of the hydrogen bond network, and the conformation of the aspartates in the active center were analyzed.

Aspartic Acid↗

PubMed1992

A molecular dynamics simulation of a polyamine-induced conformational change of DNA. A possible mechanism for the B to Z transition.

A 75ps molecular dynamics simulation has been performed on a fully solvated complex of spermine with the B DNA decamer (dGdC)5.(dGdC)5. The simulation indicates a possible mechanism by which polyamines might induce the formation of a left-handed helix, the B to Z transition. Spermine was initially located in the major groove, hydrogen bonded to the helix. During the simulation the ligand migrates deeper into the DNA, maintaining strong hydrogen bonding to the central guanine bases and destroying the Watson-Crick base pairing with their respective cytosines. Significant rotation of these and other cytosine bases was observed, in part due to interactions of the helix with the aminopropyl chains of spermine. An intermediate BII conformation might be of importance in this process.

Base Sequence↗

PubMed2003

Molecular dynamics simulations of charged and neutral lipid bilayers: treatment of electrostatic interactions.

Molecular dynamics (MD) simulations complement experimental methods in studies of the structure and dynamics of lipid bilayers. The choice of algorithms employed in this computational method represents a trade-off between the accuracy and real calculation time. The largest portion of the simulation time is devoted to calculation of long-range electrostatic interactions. To speed-up evaluation of these interactions, various approximations have been used. The most common ones are the truncation of long-range interactions with the use of cut-offs, and the particle-mesh Ewald (PME) method. In this study, several multi-nanosecond cut-off and PME simulations were performed to establish the influence of the simulation protocol on the bilayer properties. Two bilayers were used. One consisted of neutral phosphatidylcholine molecules. The other was a mixed lipid bilayer consisting of neutral phosphatidylethanolamine and negatively charged phosphatidylglycerol molecules. The study shows that the cut-off simulation of a bilayer containing charge molecules generates artefacts; in particular the mobility and order of the charged molecules are vastly different from those determined experimentally. In the PME simulation, the bilayer properties are in general agreement with experimental data. The cut-off simulation of bilayers containing only uncharged molecules does not generate artefacts, nevertheless, the PME simulation gives generally better agreement with experimental data.

Algorithms↗

PubMed2004

Diffusion and electrophoretic mobility of single-stranded RNA from molecular dynamics simulations.

Hydrodynamic properties of small single-stranded RNA homopolymers with three and six nucleotides in free solution are determined from molecular dynamics simulations in explicit solvent. We find that the electrophoretic mobility increases with increasing RNA length, consistent with experiment. Diffusion coefficients of RNA, corrected for finite-size effects and solvent viscosity, agree well with those estimated from experiments and hydrodynamic calculations. The diffusion coefficients and electrophoretic mobilities satisfy a Nernst-Einstein relation in which the effective charge of RNA is reduced by the charge of transiently bound counterions. Fluctuations in the counterion atmosphere are shown to enhance the diffusive spread of RNA molecules drifting along the direction of the external electric field. As a consequence, apparent diffusion coefficients measured by capillary zone electrophoresis can be significantly larger than the actual values at certain experimental conditions.

Computer Simulation↗

PubMed2000

Molecular dynamics simulation shows large volume fluctuations of proteins.

In this paper we present a new approach to study the volume fluctuations of proteins. From a 1 ns molecular dynamics simulation, the volume fluctuation of human lysozyme has been calculated. We used two different ways for the calculation. In the first one, the volume fluctuation is extracted directly from the trajectory. For the second one, a newly developed formalism based on principal component analysis is used. The r.m.s. volume fluctuations obtained from the two analyses agree well with each other. The isothermal intrinsic compressibility was found to be larger than the one reported by experiment. The difference is discussed and suggested to exist in the assumed uncertainty of the compressibility of hydrated water to deduce the isothermal intrinsic compressibility from the experimental value. Spectral analysis shows that low-frequency dynamics dominate the total volume fluctuation. The same aspect is found in the study using principal component analysis. This low-frequency region is related to large and slow motions of proteins. Therefore a long time dynamics simulation is necessary to describe the volume fluctuations of proteins.

Computer Simulation↗

PubMed1996

Molecular dynamics simulations of N-terminal peptides from a nucleotide binding protein.

Molecular dynamics (MD) simulations of N-terminal peptides from lactate dehydrogenase (LDH) with increasing length and individual secondary structure elements were used to study their stability in relation to folding. Ten simulations of 1-2 ns of different peptides in water starting from the coordinates of the crystal structure were performed. The stability of the peptides was compared qualitatively by analyzing the root mean square deviation (RMSD) from the crystal structure, radius of gyration, secondary and tertiary structure, and solvent accessible surface area. In agreement with earlier MD studies, relatively short (< 15 amino acids) peptides containing individual secondary structure elements were generally found to be unstable; the hydrophobic alpha 1-helix of the nucleotide binding fold displayed a significantly higher stability, however. Our simulations further showed that the first beta alpha beta supersecondary unit of the characteristic dinucleotide binding fold (Rossmann fold) of LDH is somewhat more stable than other units of similar length and that the alpha 2-helix, which unfolds by itself, is stabilized by binding to this unit. This finding suggests that the first beta alpha beta unit could function as an N-terminal folding nucleus, upon which the remainder of the polypeptide chain can be assembled. Indeed, simulations with longer units (beta-alpha-beta-alpha and beta-alpha-beta-alpha beta-beta) showed that all structural elements of these units are rather stable. The outcome of our studies is in line with suggestions that folding of the N-terminal portion of LDH in vivo can be a cotranslational process that takes place during the ribosomal peptide synthesis.

Amino Acid Sequence↗

PubMed1999

Solution conformation of phakellistatin 8 investigated by molecular dynamics simulations.

Phakellistatin 8 is a cyclic decapeptide that inhibits cancer cell growth and has sequence and structure similar to antamanide. In molecular dynamics simulations of phakellistatin 8 in water, the decapeptide ring undergoes a conformational change from the saddle-like crystal structure to a more elongated conformation by a transition of the Tyr9 main chain from the alpha L to an extended structure. This is coupled to the loss of the NH9-O6 beta-turn hydrogen bond and the transient dissociation of the Pro7-Tyr9 side-chain packing. Furthermore, the water molecule acting as a transannular bridge forms an additional hydrogen bond with phakellistatin 8, namely with the NH group of Val5 besides those already present in the crystal structure, i.e., with the NH of Ile10 and the CO of Leu6. The alpha-turn hydrogen bond between the Phe4 amide hydrogen and the Ile10 carbonyl oxygen is always present. The solution conformations of the two cyclic decapeptides are similar, in particular in the region involving the NH4-O10 alpha turn of phakellistatin 8 and the NH5-O1 alpha turn of antamanide. The simulation results suggest that in aqueous solution the conformation of phakellistatin 8 is more extended than in the crystalline state, and on a nanosecond time scale phakellistatin 8 is more flexible than antamanide.

Antineoplastic Agents↗

PubMed2004

Molecular dynamics simulation of transmembrane polypeptide orientational fluctuations.

The orientation and motion of a model lysine-terminated transmembrane polypeptide were investigated by molecular dynamics simulation. Recent 2H NMR studies of synthetic polypeptides with deuterated alanine side chains suggest that such transmembrane polypeptides undergo fast, axially symmetric reorientation about the bilayer normal but have a preferred average azimuthal orientation about the helix axis. In this work, interactions that might contribute to this behavior were investigated in a simulated system consisting of 64 molecules of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) and one alpha-helical polypeptide with the sequence acetyl-KK-(LA)11-KK-amide. In one simulation, initiated with the peptide oriented along the bilayer normal, the system was allowed to evolve for 8.5 ns at 1 atm of pressure and a temperature of 55 degrees C. A second simulation was initiated with the peptide orientation chosen to match a set of experimentally observed alanine methyl deuteron quadrupole splittings and allowed to proceed for 10 ns. Simulated alanine methyl group orientations were found to be inequivalent, a result that is consistent with 2H NMR observations of specifically labeled polypeptides in POPC bilayers. Helix tilt varied substantially over the durations of both simulations. In the first simulation, the peptide tended toward an orientation about the helix axis similar to that suggested by experiment. In the second simulation, orientation about the helix axis tended to return to this value after an excursion. These results provide some insight into how interactions at the bilayer surface can constrain reorientation about the helix axis while accommodating large changes in helix tilt.

Alanine↗

PubMed2002

Molecular dynamics simulations and KcsA channel gating.

The gating mechanism of a bacterial potassium channel, KcsA, has been investigated via multi-nanosecond molecular dynamic simulations of the channel molecules embedded in a fully solvated palmitoyloleoylphosphatidylcholine bilayer. Four events are seen in which a cation (K(+) or, in one case, Na(+)) initially present in the central cavity exits through the intracellular mouth (the presumed gate) of the channel. Whilst in the cavity a cation interacts with the sidechain T107 O gamma atom of one of the subunits prior to its exit from the channel. Secondary structure analysis as a function of time reveals a break in the helicity of one of the M2 helices. This break is expected to lend flexibility to the helices, enabling them to "open" (minimum pore radius >0.13 nm) and "close" (minimum pore radius <0.13 nm) the channel. Fluctuations in the pore radius at the intracellular gate region are of the order of 0.05 nm, with an average radius in the region of the gate of ca. 0.1 nm. However, around the time of exit of a cation, the pore widens to about 0.15 nm. The distances between the C alpha atoms of the inner helices M2 reveal a coupled increase and decrease between the opposite pair of helices at about the time of exit of the ion. This suggests a breathing motion of the M2 helices that may form the basis for a gating mechanism.

Bacterial Proteins↗

PubMed2003

Molecular dynamics simulations of A. T-rich oligomers: sequence-specific binding of Na+ in the minor groove of B-DNA.

Molecular dynamics simulations have been employed to probe the sequence-specific binding of sodium ions to the minor groove of B-DNA of three A. T-rich oligomers having identical compositions but different orders of the base pairs: C(AT)(4)G, CA(4)T(4)G, and CT(4)A(4)G. Recent experimental investigations, either in crystals or in solution, have shown that monovalent cations bind to DNA in a sequence-specific mode, preferentially in the narrow minor groove regions of uninterrupted sequences of four or more adenines (A-tracts), replacing a water molecule of the ordered hydration structure, the hydration spine. Following this evidence, it has been hypothesized that in A-tracts these events may be responsible for structural peculiarities such as a narrow minor groove and a curvature of the helix axis. The present simulations confirm a sequence specificity of the binding of sodium ions: Na(+) intrusions in the first layer of hydration of the minor groove, with long residence times, up to approximately 3 ns, are observed only in the minor groove of A-tracts but not in the alternating sequence. The effects of these intrusions on the structure of DNA depend on the ion coordination: when the ion replaces a water molecule of the spine, the minor groove becomes narrower. Ion intrusions may also disrupt the hydration spine modifying the oligomer structure to a large extent. However, in no case intrusions were observed to locally bend the axis toward the minor groove. The simulations also show that ions may reside for long time periods in the second layer of hydration, particularly in the wider regions of the groove, often leading to an opening of the groove.

Adenine↗

PubMed2006

Modeling the 3D structure of wheat subtilisin/chymotrypsin inhibitor (WSCI). Probing the reactive site with two susceptible proteinases by time-course analysis and molecular dynamics simulations.

Comparative modeling and time-course hydrolysis experiments have been applied to investigate two enzyme-inhibitor complexes formed between the wheat subtilisin-chymotrypsin inhibitor (WSCI) and two susceptible proteinases. WSCI represents the first case of a wheat protein inhibitor active against animal chymotrypsins and bacterial subtilisins. The model was created using as template structure that of the CI-2A inhibitor from barley (PDB code: 2CI2), which shares 87% sequence identity with WSCI. Under these conditions of high similarity, the comparative modeling approach can be successfully applied. We predicted the WSCI 3D model and used it to investigate enzyme-inhibitor complex systems. Experimental observations indicated that chymotrypsin, but not subtilisin, in addition to cleavage at the primary reactive site Met48-Glu49, is able to hydrolyze a second peptide bond between Phe58 and Val59. Here, we report on cleavage of the peptide bond at the inhibitor's reactive site (Met48-Glu49) determined using time-course hydrolysis experiments; the same event was investigated for both subtilisin/WSCI and chymotrypsin/WSCI complexes using molecular dynamics simulations. The molecular details of the initial inhibitor-enzyme interactions, as well as of the changes observed during the simulations, allow us to speculate on the different fates of the two WSCI-proteinase complexes.

Amino Acid Sequence↗

PubMed2003

Nonspecific interaction forces at water-membrane interface by forced molecular dynamics simulations.

Nonspecific interactions are the main driving forces for the behavior of molecules with great affinity for biologic membranes. To investigate not only the molecular details of these interactions but to estimate their magnitude as well, the theoretical method of Forced Molecular Dynamics Simulations, based on the Atomic Force Spectroscopy experimental technique, was applied. In this approach, an additional one-dimensional elastic force, representing the cantilever probe, was incorporated to the force field of a Molecular Dynamics computational program. This force represents a spring fixed on one end to a selected atom of the molecule; the other end of the spring is displaced at constant velocity to pull the molecule out of the membrane. The force experimented by the molecule due to the spring, is proportional to the spring elongation relative to its equilibrium position. This value is registered during the entire simulation, and its maximum value will determine the molecule-membrane interaction force. Nonexplicit medium simulations were carried out. Polar and apolar media were considered according to their polarizability degree and a specific dielectric constant value was assigned. In this approach, the membrane was considered as the apolar region limited by two flat surfaces with a polar aqueous medium. The potential energy discontinuity at the interfaces was smoothed by considering the polarization-induced effects using the image method. The results of this methodology are presented using a small system, a single Alanine amino acid model, which enables extended simulations in a microsecond time scale. The confinement of this amino acid at the interface reduces its degrees of freedom and forces it to adopt one of the six defined conformations. A correlation between these stable structures at the water-membrane interface and the interaction force value was determined.

Alanine↗

PubMed1997

Effect of periodic box size on aqueous molecular dynamics simulation of a DNA dodecamer with particle-mesh Ewald method.

The particle-mesh Ewald (PME) method is considered to be both efficient and accurate for the evaluation of long-range electrostatic interactions in large macromolecular systems being studied by molecular dynamics simulations. This method assumes "infinite" periodic boundary conditions resembling the symmetry of a crystal environment. Can such a "solid-state" method accurately portray a macromolecular solute such as DNA in solution? To address this issue, we have performed three 1500-ps PME molecular dynamics (MD) simulations, each with a different box size, on the d(CGCGA6CG)-(CGT6CGCG) DNA dodecamer. The smallest box had the DNA solvated by a layer of water molecules of at least 5 A along each orthogonal direction. The intermediate size box and the largest box had the DNA solvated by a layer of water molecules of at least 10 A and 15 A, respectively, along each orthogonal direction. The intermediate size box in the present study is similar to the box size currently chosen by most workers in the field. Based on a comparison of RMSDs and curvature for this single DNA dodecamer sequence, the larger two box sizes do not appear to afford any extra benefit over the smallest box. The implications of this finding are briefly discussed.

Base Sequence↗

PubMed2005

n-Heptane under pressure: structure and dynamics from molecular simulations.

Atomistic molecular dynamics simulations have been performed in the isothermal-isobaric ensemble to explore the phase behavior of n-heptane. Motivated by recent high-pressure spectroscopic experiments on n-heptane, the present work aims at understanding the liquid-solid and the alluded to solid-solid transitions upon increasing pressure. Starting from the stabilized solid phase at 300 K and 10 kbar, we have investigated the range of these two transitions by a gradual decrease and increase of pressure, respectively. Although the solid-liquid transition has clear signatures such as the formation of gauche defects along the molecular backbone, the present model does not show any sign of a first-order solid-solid transition at high pressures. However, interesting changes in the environment around methyl groups and in their dynamics are observed. These have been substantiated by calculations of the vibrational density of states obtained from a normal-mode analysis and from the simulation trajectory.

Journal Article↗

PubMed1999

Steered molecular dynamics simulations of force-induced protein domain unfolding.

Steered molecular dynamics (SMD), a computer simulation method for studying force-induced reactions in biopolymers, has been applied to investigate the response of protein domains to stretching apart of their terminal ends. The simulations mimic atomic force microscopy and optical tweezer experiments, but proceed on much shorter time scales. The simulations on different domains for 0.6 nanosecond each reveal two types of protein responses: the first type, arising in certain beta-sandwich domains, exhibits nanosecond unfolding only after a force above 1,500 pN is applied; the second type, arising in a wider class of protein domain structures, requires significantly weaker forces for nanosecond unfolding. In the first case, strong forces are needed to concertedly break a set of interstrand hydrogen bonds which protect the domains against unfolding through stretching; in the second case, stretching breaks backbone hydrogen bonds one by one, and does not require strong forces for this purpose. Stretching of beta-sandwich (immunoglobulin) domains has been investigated further revealing a specific relationship between response to mechanical strain and the architecture of beta-sandwich domains.

Computer Simulation↗