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 1,261 records · Page 70Linked to original sources

PubMed2001

Molecular dynamics simulation of the human U2B" protein complex with U2 snRNA hairpin IV in aqueous solution.

A 2200-ps molecular dynamics (MD) simulation of the U2 snRNA hairpin IV/U2B" complex was performed in aqueous solution using the particle mesh Ewald method to consider long-range electrostatic interactions. To investigate the interaction and recognition process between the RNA and protein, the free energy contributions resulting from individual amino acids of the protein component of the RNA/protein complex were calculated using the recently developed glycine-scanning method. The results revealed that the loop region of the U2 snRNA hairpin IV interacted mainly with three regions of the U2B" protein: 1) beta 1-helix A, 2) beta 2-beta 3, and 3) beta 4-helix C. U2 snRNA hairpin IV bound U2B" in a similar orientation as that previously described for U1 snRNA with the U1A' protein; however, the details of the interaction differed in several aspects. In particular, beta 1-helix A and beta 4-helix C in U2B" were not observed to interact with RNA in the U1A' protein complex. Most of the polar and charged residues in the interacting regions had larger mutant free energies than the nonpolar residues, indicating that electrostatic interactions were important for stabilizing the RNA/protein complex. The interaction was further stabilized by a network of hydrogen bonds and salt bridges formed between RNA and protein that was maintained throughout the MD trajectory. In addition to the direct interactions between RNA and the protein, solvent-mediated interactions also contributed significantly to complex stability. A detailed analysis of the ordered water molecules in the hydration of the RNA/protein complex revealed that bridged water molecules reside at the interface of RNA and protein as long as 2100 ps in the 2200-ps trajectory. At least 20 bridged water molecules, on average, contributed to the instantaneous stability of the RNA/protein complex. The stabilizing interaction energy due to bridging water molecules was obtained from ab initio Hartree-Fock and density functional theory calculations.

Autoantigens↗

PubMed2006

Structural investigation of syringomycin-E using molecular dynamics simulation and NMR.

Syringomycin-E (SR-E) is a cyclic lipodepsinonapeptide produced by certain strains of the bacterium Pseudomonas syringae pv. syringae. It shows inhibitory effects against many fungal species, including human pathogens. Its primary biological target is the plasma membrane, where it forms channels comprised of at least six SR-E molecules. The high-resolution structure of SR-E and the structure of the channels are currently not known. In this paper, we investigate in atomic detail the molecular features of SR-E in water by NMR and in water and octane by molecular dynamics simulation (MD). We built a model of the peptide and examined its structure in water and octane in 200 ns MD simulations both with and without distance restraints derived from NMR NOE data. The resulting trajectories show good agreement with the measured NOEs and circular dichroism data from the literature and provide atomistic models of SR-E that are an important step toward a better understanding of the antifungal and antibacterial activity of this peptide.

Biophysical Phenomena↗

PubMed1998

Unfolding of titin immunoglobulin domains by steered molecular dynamics simulation.

Titin, a 1-microm-long protein found in striated muscle myofibrils, possesses unique elastic and extensibility properties in its I-band region, which is largely composed of a PEVK region (70% proline, glutamic acid, valine, and lysine residue) and seven-strand beta-sandwich immunoglobulin-like (Ig) domains. The behavior of titin as a multistage entropic spring has been shown in atomic force microscope and optical tweezer experiments to partially depend on the reversible unfolding of individual Ig domains. We performed steered molecular dynamics simulations to stretch single titin Ig domains in solution with pulling speeds of 0.5 and 1.0 A/ps. Resulting force-extension profiles exhibit a single dominant peak for each Ig domain unfolding, consistent with the experimentally observed sequential, as opposed to concerted, unfolding of Ig domains under external stretching forces. This force peak can be attributed to an initial burst of backbone hydrogen bonds, which takes place between antiparallel beta-strands A and B and between parallel beta-strands A' and G. Additional features of the simulations, including the position of the force peak and relative unfolding resistance of different Ig domains, can be related to experimental observations.

Amino Acid Sequence↗

PubMed2004

Molecular dynamics simulations of transducin: interdomain and front to back communication in activation and nucleotide exchange.

The dynamic events that underlie the nucleotide exchange process for the Galpha subunit of transducin (Galpha(t)) were studied with nanosecond time-scale molecular dynamics simulations. The modeled systems include the active and inactive forms of the wild-type Galpha(t) and three of its mutants (GDP-bound form only): F332A, A322S, and Q326A that are known to exhibit various degrees of enhancement of their basal and receptor-catalyzed rates of nucleotide exchange (150-fold, 70-fold and WT-like, respectively). The results of these computational experiments reveal a number of nucleotide-dependent structural and dynamic changes (involving the alpha(B)-alpha(C) loop, the inter-domain orientation of the helical and GTPase domains and the alpha(5) helix) that were not observed in the various crystal structures of Galpha(t). Notably, the results show the existence of a front to back communication device (involving the beta(2)-beta(3) hairpin, the alpha(1) helix and the alpha(5) helix), strategically located near all elements susceptible to be involved in receptor-mediated activation/nucleotide exchange. The wild-type simulations suggest that the dynamic interplay between the elements of this device would be critical for the activation of the Galpha(t) subunit. This inference is confirmed by the results of the computational experiments on the mutants that show that even in their GDP-bound forms, the A322S and F332A mutants acquire an "active-like" structure and dynamics phenotype. The same is not true for the Q326A mutant whose structural and dynamic properties remain similar to those of the GDP-bound WT. Taken together the results suggest a nucleotide exchange mechanism, analogous to that found in the Arf family GTPases, in which a partially activated state, achievable from a receptor-mediated action of the front to back communication device either by displacement of the C-terminal alpha(5) helix, of the N-terminal alpha(N) helix, or of the Gbetagamma subunit, could precede the dissociation of GDP from the native Galpha subunit.

Enzyme Activation↗

PubMed2005

Molecular dynamics simulations of evolved collective motions of atoms in the myosin motor domain upon perturbation of the ATPase pocket.

A crucial point for mechanical force generation in actomyosin systems is how the energy released by ATP hydrolysis in the myosin motor domain gives rise to the movement of the myosin head along the actin filament. We assumed the signal of the ATP hydrolysis to be transmitted as modulated atomic vibrations from the nucleotide-binding site throughout the myosin head, and carried out 1-ns all-atom molecular dynamics simulations for that signal transmission. We distributed the released energy to atoms located around the ATPase pocket as kinetic energies and examined how the effect of disturbance extended throughout the motor domain. The result showed that the disturbance signal extended over the motor domain in 150 ps and induced slowly varying collective motions of atoms at the actin-binding site and the junction with the neck, both of which are relevant to the movement of the myosin head along the actin filament. We also performed a principal component analysis of thermal atomic motions for the motor domain, and the first principal component was consistent with the response to the disturbance given to the ATPase pocket.

Adenosine Triphosphatases↗

PubMed2001

At nonzero temperatures, stacked structures of methylated nucleic acid base pairs and microhydrated nonmethylated nucleic acid base pairs are favored over planar hydrogen-bonded structures: a molecular dynamics simulations study.

The dynamic structure of all ten possible nucleic acid (NA) base pairs and methylated NA base pairs hydrated by a small number of water molecules (from 1 to 16) was determined by using molecular dynamics simulations in the NVE microcanonical and NVT canonical ensembles with the Cornell force field (W. D. Cornell, P. Cieplak, C. I. Bayly, I. R. Gould, K. M. Merz, D. M. Ferguson, D. C. Spellmeyer, T. Fox, J. E. Caldwell, P. Kollman, J. Am. Chem. Soc. 1995, 117, 5179). The presence of one water molecule does not affect the structure of any hydrogen-bonded (H-bonded) nonmethylated base pair. An equal population of H-bonded and stacked structures of adenine...adenine, adenine...guanine and adenine... thymine pairs is reached if as few as two water molecules are present, while obtaining equal populations of these structures in the case of adenine...cytosine, cytosine...thymine, guanine... guanine and guanine...thymine required the presence of four water molecules, and in the case of guanine...cytosine, six. A comparable population of planar, H-bonded and stacked structures for cytosine...cytosine and thymine... thymine base pairs was only obtained if at least eight water molecules hydrated a pair. Methylation of bases changed the situation dramatically and stacked structures were favoured over H-bonded ones even in the absence of water molecules in most cases. Only in the case of methyl cytosine...methyl cytosine, methyl guanine...methyl guanine and methyl guanine...methyl cytosine pairs were two, two or six water molecules, respectively, needed in order to obtain a comparable population of planar, H-bonded and stacked structures. We believe that these results give clear evidence that the preferred stacked structure of NA base pairs in the microhydrated environment, and also apparently in a regular solvent, is due to the hydrophilic interaction of a small number of water molecules. In the case of methylated bases, it is also due to the fact that the hydrogen atoms most suitable for the formation of H-bonds have been replaced by a methyl group. A preferred stacked structure is, thus, not due to a hydrophobic interaction between a large bulk of water molecules and the base pair, as believed.

Base Pairing↗

PubMed2002

Parallel-in-time molecular-dynamics simulations.

While there have been many progress in the field of multiscale simulations in the space domain, in particular, due to efficient parallelization techniques, much less is known in the way to perform similar approaches in the time domain. In this paper we show on two examples that, provided we can describe in a rough but still accurate way the system under consideration, it is indeed possible to parallelize molecular dynamics simulations in time by using the recently introduced pararealalgorithm. The technique is most useful for ab initio simulations.

Journal Article↗

PubMed1999

Side-chain conformational disorder in a molten globule: molecular dynamics simulations of the A-state of human alpha-lactalbumin.

Molten globules are compact, partially folded forms of proteins consisting of an ensemble of interconverting conformers with disorder in the side-chain packing across the ensemble. Using insights from experimental data a strategy has been devised to describe this side-chain disorder in a molten globule. Structures of human alpha-lactalbumin have been generated with significantly different side-chain packings to those observed in the native protein by changing all the chi1 torsion angles. Molecular dynamics simulations have been carried out starting from these structures using conditions under which the molten globule is seen experimentally, i.e. at low pH and in the absence of calcium. In each simulation, low energy conformers were generated which are compact and retain substantial secondary structure and an overall native fold, but have very different side-chain contacts; there are all-atom root-mean-square differences of up to 0.5 nm between the different structures. The structures from the different simulations taken together provide an initial description of possible contributors to the ensemble of conformers in the molten globule state which is consistent, at least in general terms, with experimental data.

Calcium↗

PubMed2005

Release of ADP from the catalytic subunit of protein kinase A: a molecular dynamics simulation study.

Substrate phosphorylation by cAMP-dependent-protein kinase A (protein kinase A, PKA) has been studied extensively. Phosphoryl transfer was found to be fast, whereas ADP release was found to be the slow, rate-limiting step. There is also evidence that ADP release may be preceded by a partially rate-limiting conformational change. However, the atomic details of the conformational change and the mode of ADP release are difficult to obtain experimentally. In this work, we studied ADP release from PKA by carrying out molecular dynamics simulations with different pulling forces applied to the ligand. The detailed ADP release pathway and the associated conformational changes were analyzed. The ADP release process was found to involve a swinging motion with the phosphate of ADP anchored to the Gly-rich loop, so that the more buried adenine base and ribose ring came out before the phosphate. In contrast to the common belief that a hinge-bending motion was responsible for the opening of the ligand-binding cleft, our simulations showed that the small lobe exhibited a large amplitude "rocking" motion when the ligand came out. The largest conformational change of the protein was observed at about the first quarter time point along the release pathway. Two prominent intermediate states were observed in the release process.

Adenosine Diphosphate↗

PubMed2003

A box-counting-based algorithm for computing Shannon entropy in molecular dynamics simulations.

A box-counting-based algorithm (SEBC) has been developed for the numerical computation of the Shannon entropy from samples of continuous functions. Its performance was tested by applying it to several samples of known continuous distribution functions. The results obtained with SEBC reproduced those obtained by analytical or numerical integration. SEBC was also employed for computing the Shannon entropies of the steric energy, Sh(E(S)), of several amino acids from their in vacuo NVE molecular dynamics simulations using the AMBER-4 force field. The results obtained correlate linearly with the experimental standard thermodynamic entropies of these compounds. This work points to the possibility of introducing straightforward and reliable calculations of thermodynamic entropies from empirical linear relationships with Sh(E(S)) obtained from MD simulations.

Journal Article↗

PubMed2006

Selectivity principle of the ligand escape process from a two-gate tunnel in myoglobin: molecular dynamics simulation.

We proposed a selectivity principle for the ligand escape process from two fluctuating bottlenecks in a cavity with a multigate inside a myoglobin pocket. Our previous analytical theory proposed a fluctuating bottleneck model for a Brownian particle passing through two gates on a cavity surface of an enzyme protein and has determined the escape rate in terms of the time-dependent gate function and the competition effect. It illustrated that with two (or more than two) gates on a cavity surface the gate modulation, which is controlled by protein fluctuation, dominates the ligand escape pathway. We have performed a molecular dynamics simulation to investigate the selectivity principle of the ligand escape process from two-gate tunnel in myoglobin. The simulation results confirm our theoretical conjecture. It indicates that the escape process is actually entropy driven, and the ligand escape pathway is chosen via the gate modulation. This suggests an interesting intrinsic property, that is, the oxymyoglobin tertiary structure is favorable to the departure of the ligand from one direction rather than through a biased random walk.

Computer Simulation↗

PubMed1999

Forced unfolding of fibronectin type 3 modules: an analysis by biased molecular dynamics simulations.

Titin, an important constituent of vertebrate muscles, is a protein of the order of a micrometer in length in the folded state. Atomic force microscopy and laser tweezer experiments have been used to stretch titin molecules to more than ten times their folded lengths. To explain the observed relation between force and extension, it has been suggested that the immunoglobulin and fibronectin domains unfold one at a time in an all-or-none fashion. We use molecular dynamics simulations to study the forced unfolding of two different fibronectin type 3 domains (the ninth, 9Fn3, and the tenth, 10Fn3, from human fibronectin) and of their heterodimer of known structure. An external biasing potential on the N to C distance is employed and the protein is treated in the polar hydrogen representation with an implicit solvation model. The latter provides an adiabatic solvent response, which is important for the nanosecond unfolding simulation method used here. A series of simulations is performed for each system to obtain meaningful results. The two different fibronectin domains are shown to unfold in the same way along two possible pathways. These involve the partial separation of the "beta-sandwich", an essential structural element, and the unfolding of the individual sheets in a stepwise fashion. The biasing potential results are confirmed by constant force unfolding simulations. For the two connected domains, there is complete unfolding of one domain (9Fn3) before major unfolding of the second domain (10Fn3). Comparison of different models for the potential energy function demonstrates that the dominant cohesive element in both proteins is due to the attractive van der Waals interactions; electrostatic interactions play a structural role but appear to make only a small contribution to the stabilization of the domains, in agreement with other studies of beta-sheet stability. The unfolding forces found in the simulations are of the order of those observed experimentally, even though the speed of the former is more than six orders of magnitude greater than that used in the latter.

Dimerization↗

PubMed2005

Out-of-plane motions in open sliding clamps: molecular dynamics simulations of eukaryotic and archaeal proliferating cell nuclear antigen.

Sliding clamps are ring-like multimeric proteins that encircle duplex DNA and serve as mobile DNA-bound platforms that are essential for efficient DNA replication and repair. Sliding clamps are placed on DNA by clamp loader complexes, in which the clamp-interacting elements are organized in a right-handed spiral assembly. To understand how the flat, ring-like clamps might interact with the spiral interaction surface of the clamp loader complex, we have performed molecular dynamics simulations of sliding clamps (proliferating cell nuclear antigen from the budding yeast, humans, and an archaeal species) in which we have removed one of the three subunits so as to release the constraint of ring closure. The simulations reveal significant structural fluctuations corresponding to lateral opening and out-of-plane distortions of the clamp, which result principally from bending and twisting of the beta-sheets that span the intermolecular interfaces, with smaller but similar contributions from beta-sheets that span the intramolecular interfaces within each subunit. With the integrity of these beta-sheets intact, the predominant fluctuations seen in the simulations are oscillations between lateral openings and right-handed spirals. The tendency for clamps to adopt a right-handed spiral conformation implies that once opened, the conformation of the clamp can easily match the spiraling of clamp loader subunits, a feature that is intrinsic to the recognition of DNA and subsequent hydrolysis of ATP by the clamp-bound clamp loader complex.

Archaeal Proteins↗

PubMedDate not supplied

A molecular dynamics simulation of the (dG)6 . (dC)6 minihelix including counterions and water.

The results of a 60 ps molecular dynamics (MD) simulation of (dG)6.(dC)6 including 10 Na+ counterions and 292 water molecules are presented. All backbone angles and helix parameters for the hexamer are reported in this paper along with trajectory plots of selected angles. Hydrogen bonding between the bases along the helical axis was observed to fluctuate with time, showing the dynamic nature of the base-pairing interaction. These fluctuations gave rise to unusual hydrogen-bonding patterns. Good intrastrand base stacking and no interstrand base stacking were also observed. The hexamer minihelix retains an essentially B-DNA conformation throughout the entire simulation even though some helix parameters and backbone angles do not have strict B-DNA values. The most striking feature obtained from the simulation was a high propeller twist, which resulted in a narrow minor groove for the minihelix. It is proposed that (dG)n.(dC)n sequences are resistant to DNAase I because of this narrow minor groove in dilute aqueous solution.

Base Sequence↗

PubMed2005

Pressure denaturation of apomyoglobin: a molecular dynamics simulation study.

The effect of pressure on the structure and mobility of Sperm Wale Apomyoglobin was studied by Molecular Dynamics computer simulation at 1 bar and 3 kbar (1 atm=1.01325 bar=101.325 kPa). The results are in good agreement with the available experimental data, allowing further analysis of other features of the effect of pressure on the protein solution. From the analysis of Secondary Structures (SS) along the trajectories it is observed that alpha-helixes are favoured under pressure at the expense of bends, turns and 3-helixes. The studies of mobility show that although the general mobility is restricted under pressure this is not true for some particular residues. The studies of tertiary structure show important conformational changes. The evolution of the Solvent Accessed Surface (SAS) with pressure shows a notorious increase due almost completely to a biased raise in the hydrophobic area exposed, which consequently shows that the hydrophobic interaction is considerably weaker under high hydrostatic pressure conditions.

Animals↗

PubMed2002

Interaction of bombolitin II with a membrane-mimetic environment: an NMR and molecular dynamics simulation approach.

Bombolitins are five natural heptadecapeptides originally isolated from the venom of a bumblebee. They induce lysis of erythrocytes and liposomes and increase the activity of phospholipase A(2) (PLA(2)), that plays an important role in the early steps of the inflammatory process. It has been proposed that PLA(2) activation depends on the alteration of the physical state of the membrane. Bombolitin II folds into an alpha-helix in a membrane mimicking environment constituted by sodium dodecyl sulfate micelles (Macromol. Chem. Phys., 196 (1995) 2827). In the present work, the topological orientation of the peptide relative to the micelle was determined, using three spin probes localized in different positions of the water/micelle system. The reduction in intensity of the 1H NMR signals clearly demonstrated that the peptide is located on the surface of the micelle, with its helical axis parallel to the interface. Only a small portion of the helix is exposed to the aqueous environment. Results from NMR experiments were confirmed by molecular dynamics simulations, performed using a two-phase water/decane simulation cell. The timescale for the reorientation of the peptide was between 120 and 450 ps, depending on the starting position of the peptide.

Molecular Mimicry↗

PubMed2005

Early events in protein aggregation: molecular flexibility and hydrophobicity/charge interaction in amyloid peptides as studied by molecular dynamics simulations.

In a previous article (Zbilut et al., Biophys J 2003;85:3544-3557), we demonstrated how an aggregation versus folding choice could be approached considering hydrophobicity distribution and charge. In this work, our aim is highlighting the mutual interaction of charge and hydrophobicity distribution in the aggregation process. Use was made of two different peptides, both derived from a transmembrane protein (amyloid precursor protein; APP), namely, Abeta(1-28) and Abeta(1-40). Abeta(1-28) has a much lower aggregation propensity than Abeta(1-40). The results obtained by means of molecular dynamics simulations show that, when submitted to the most "aggregation-prone" environment, corresponding to the isoelectric point and consequently to zero net charge, both peptides acquire their maximum flexibility, but Abeta(1-40) has a definitely higher conformational mobility than Abeta(1-28). The absence of a hydrophobic "tail," which is the most mobile part of the molecule in Abeta(1-40), is the element lacking in Abeta(1-28) for obtaining a "fully aggregating" phenotype. Our results suggest that conformational flexibility, determined by both hydrophobicity and charge effect, is the main mechanistic determinant of aggregation propensity.

Amyloid↗

PubMed2004

A molecular dynamics simulation study of the alpha-relaxation in a 1,4-polybutadiene melt as probed by the coherent dynamic structure factor.

The dynamic coherent structure factor Scoh(q,t) for a 1,4-polybutadiene (PBD) melt has been investigated using atomistic molecular dynamics simulations. The relaxation of Scoh(q,t) at q = 1.44 angstroms(-1) and q = 2.72 angstroms(-1), corresponding to the first and second peaks in the static structure factor for PBD, was studied in detail over a wide range of temperature. It was found that time-temperature superposition holds for the alpha-relaxation for both q values over a wide temperature range and that the alpha-relaxation can be well described by a stretched (Kohlrauch-William-Watts) exponential with temperature independent but q dependent amplitude and stretching exponent. The alpha-relaxation times for both q values were found to exhibit the same non-Arrhenius temperature dependence, indicating that the same physical processes are responsible for relaxation on both length scales. The alpha-relaxation time was found to depend strongly upon the dynamical range of data utilized in determining the relaxation time, accounting for qualitative discrepancies between alpha-relaxation times reported here and those extracted for PBD from experimentally measured Scoh(q,t).

Journal Article↗