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,351 records · Page 75Linked to original sources
PubMedDate not supplied

[Study of the compact denatured state of a protein by molecular dynamics simulation].

The native state can be considered as a unique conformation of the protein molecule with the lowest free energy of residue contacts. In this case, all other conformations correspond to the denatured state. The degree of their compactness varies significantly. Under folding conditions, the compact denatured state rather than the random coil is in equilibrium with native protein. The balance between the main forces of protein folding, the solvophobic interactions and conformational entropy, suggests that some properties of the compact denatured state are close to those of native protein, whereas other properties are close to those of the random coil. To investigate the molecular structure of the compact denatured state, the method of molecular dynamics simulation seems to be very useful.

Protein Conformation↗
PubMed2005

Dynamics of poly(ethylene oxide) in a blend with poly(methyl methacrylate): a quasielastic neutron scattering and molecular dynamics simulations study.

In this paper, we have addressed the question of the dynamic miscibility in a blend characterized by very different glass-transition temperatures, Tg, for the components: poly(ethylene oxide) and poly(methyl methacrylate) (PEO/PMMA). The combination of quasielastic neutron scattering with isotopic labeling and fully atomistic molecular dynamics simulations has allowed us to selectively investigate the dynamics of the two components in the picosecond-10 nanoseconds scale at temperatures close and above the Tg of the blend. The main focus was on the PEO component, i.e., that of the lowest Tg, but first we have characterized the dynamics of the other component in the blend and of the pure PEO homopolymer as reference. In the region investigated, the dynamics of PMMA in the blend is strongly affected by the alpha-methyl rotation; an additional process detected in the experimental window 65 K above the blend-Tg can be identified as the merged alphabeta process of this component that shows strong deviations from Gaussian behavior. On the other hand, pure PEO displays entropy driven dynamics up to very large momentum transfers. Such kind of motion seems to freeze when the PEO chains are in the blend. There, we have directly observed a very heterogeneous and moreover confined dynamics for the PEO component. The presence of the hardly moving PMMA matrix leads to the creation of little pockets of mobility where PEO can move. The characteristic size of such confined islands of mobility might be estimated to be of approximately 1 nm. These findings are corroborated by the simulation study, which has been an essential support and guide in our data analysis procedure.

Journal Article↗
PubMed2004

Equilibrium structure and folding of a helix-forming peptide: circular dichroism measurements and replica-exchange molecular dynamics simulations.

We have performed experimental measurements and computer simulations of the equilibrium structure and folding of a 21-residue alpha-helical heteropeptide. Far ultraviolet circular dichroism spectroscopy is used to identify the presence of helical structure and to measure the thermal unfolding curve. The observed melting temperature is 296 K, with a folding enthalpy of -11.6 kcal/mol and entropy of -39.6 cal/(mol K). Our simulations involve 45 ns of replica-exchange molecular dynamics of the peptide, using eight replicas at temperatures between 280 and 450 K, and the program CHARMM with a continuum solvent model. In a 30-ns simulation started from a helical structure, conformational equilibrium at all temperatures was reached after 15 ns. This simulation was used to calculate the peptide melting curve, predicting a folding transition with a melting temperature in the 330-350 K range, enthalpy change of -10 kcal/mol, and entropy change of -30 cal/(mol K). The simulation results were also used to analyze the peptide structural fluctuations and the free-energy surface of helix unfolding. In a separate 15-ns replica-exchange molecular dynamics simulation started from the extended structure, the helical conformation was first attained after approximately 2.8 ns, and equilibrium was reached after 10 ns of simulation. These results showed a sequential folding process with a systematic increase in the number of hydrogen bonds until the helical state is reached, and confirmed that the alpha-helical state is the global free-energy minimum for the peptide at low temperatures.

Circular Dichroism↗
PubMed2001

Molecular dynamics simulations of the conformational changes of the glutamate receptor ligand-binding core in the presence of glutamate and kainate.

Excitatory synaptic transmission is mediated by ionotropic glutamate receptors (iGluRs) through the induced transient opening of transmembrane ion channels. The three-dimensional structure of the extracellular ligand-binding core of iGluRs shares the overall features of bacterial periplasmic binding proteins (PBPs). In both families of proteins, the ligand-binding site is arranged in two domains separated by a cleft and connected by two peptide stretches. PBPs undergo a typical hinge motion of the two domains associated with ligand binding that leads to a conformational change from an open to a closed form. The common architecture suggests a similar closing mechanism in the ligand-binding core of iGluRs induced by the binding of specific agonists. Starting from the experimentally determined kainate-bound closed form of the S1S2 GluR2 construct, we have studied by means of molecular dynamics simulations the opening motion of the ligand-binding core in the presence and in the absence of both glutamate and kainate. Our results suggest that the opening/closing interdomain hinge motions are coupled to conformational changes in the insertion region of the transmembrane segments. These changes are triggered by the interaction of the agonists with the essential Glu 209 residue. A plausible mechanism for the coupling of agonist binding to channel gating is discussed.

Computer Simulation↗
PubMed2003

Molecular dynamics simulations of the NGF-TrkA domain 5 complex and comparison with biological data.

The nerve growth factor (NGF) is an important pharmacological target for Alzheimer's and other neurodegenerative diseases. Its action derives partly from its binding to the tyrosine kinase A receptor (TrkA). Here we study energetics and dynamics of the NGF-TrkA complex by carrying out multinanosecond molecular dynamics simulations, accompanied by electrostatic calculations based on the Poisson-Boltzmann equation. Our calculations, which are based on the x-ray structure of the complex, suggest that some of the mutations affecting dramatically the affinity of the complex involve residues that form highly favorable, direct or water-mediated hydrogen bond interactions at the ligand-receptor interface and, in some cases, that also critically participate to the large-scale motions of the complex. Furthermore, our calculations offer a rationale for the small effect on binding affinity observed upon specific mutations involving large changes in electrostatics (i.e., the charged-to-neutral mutations). Finally, these calculations, used along with the mutagenesis data, provide a basis for designing new peptides that mimic NGF in TrkA binding function.

Binding Sites↗
PubMed2001

Carbohydrate-protein recognition: molecular dynamics simulations and free energy analysis of oligosaccharide binding to concanavalin A.

Carbohydrate ligands are important mediators of biomolecular recognition. Microcalorimetry has found the complex-type N-linked glycan core pentasaccharide beta-GlcNAc-(1-->2)-alpha-Man-(1-->3)-[beta-GlcNAc-(1-->2)-alpha-Man-(1-->6)]-Man to bind to the lectin, Concanavalin A, with almost the same affinity as the trimannoside, Man-alpha-(1-->6)-[Man-alpha-(1-->3)]-Man. Recent determination of the structure of the pentasaccharide complex found a glycosidic linkage psi torsion angle to be distorted by 50 degrees from the NMR solution value and perturbation of some key mannose-protein interactions observed in the structures of the mono- and trimannoside complexes. To unravel the free energy contributions to binding and to determine the structural basis for this degeneracy, we present the results of a series of nanosecond molecular dynamics simulations, coupled to analysis via the recently developed MM-GB/SA approach (Srinivasan et al., J. Am. Chem. Soc. 1998, 120:9401-9409). These calculations indicate that the strength of key mannose-protein interactions at the monosaccharide site is preserved in both the oligosaccharides. Although distortion of the pentasaccharide is significant, the principal factor in reduced binding is incomplete offset of ligand and protein desolvation due to poorly matched polar interactions. This analysis implies that, although Concanavalin A tolerates the additional 6 arm GlcNAc present in the pentasaccharide, it does not serve as a key recognition determinant.

Acetylation↗
PubMed2003

Molecular dynamics simulations of lignin peroxidase in solution.

The dynamical and structural properties of lignin peroxidase and its Trp171Ala mutant have been investigated in aqueous solution using molecular dynamics (MD) simulations. In both cases, the enzyme retained its overall backbone structure and all its noncovalent interactions in the course of the MD simulations. Very interestingly, the analysis of the MD trajectories showed the presence of large fluctuations in correspondence of the residues forming the heme access channel; these movements enlarge the opening and facilitate the access of substrates to the enzyme active site. Moreover, steered molecular dynamics docking simulations have shown that lignin peroxidase natural substrate (veratryl alcohol) can easily approach the heme edge through the access channel.

Benzyl Alcohols↗
PubMed2003

A molecular dynamics simulation study of oriented DNA with polyamine and sodium counterions: diffusion and averaged binding of water and cations.

Four different molecular dynamics (MD) simulations have been performed for ordered DNA decamers, d(5'-ATGCAGTCAG)*d(5'-TGACTGCATC). The counterions were the two natural polyamines spermidine3+ (Spd3+) and putrescine2+ (Put2+), the synthetic polyamine diaminopropane2+ (DAP2+) and Na+. The simulation set-up corresponds to an infinite array of parallel DNA mimicking the state in oriented DNA fibers or crystals. This work describes general properties of polyamine and Na+ binding to DNA. Simulated diffusion coefficients show satisfactory agreement with experimental NMR diffusion data of comparable systems. The interaction of the polyamines with DNA is dynamic in character and the cations mostly form short-lived contacts with the electronegative binding sites of DNA. Polyamines, Na+ and water interact most frequently with the charged phosphate atoms with preference for association from the minor groove side with O1P over O2P. There is a strong anti-correlation in the cation binding to the electronegative groups of DNA, i.e. the presence of a cation near one of the DNA sites repels other cations from binding to this and to the other sites separated by <7.5 A from each other. In contrast to the other polyamines, DAP2+ is able to form 'bridges' connecting neighboring phosphate groups along the DNA strand. A small fraction of DAP2+ and Put2+ can be found in the major grooves, while Spd3+ is absent there. The results of the MD simulations reveal principal differences in the polyamine-DNA interactions between the natural (Spd3+, Put2+ and spermine4+) and synthetic (DAP2+) polyamines.

Binding Sites↗
PubMed2005

Characterization of the stereochemical selectivity of beta-hairpin formation by molecular dynamics simulations.

The stability of secondary structure motifs found in proteins is influenced by the choice of the configuration of the chiral centers present in the amino acid residues (i.e., D vs L). Experimental studies showed that the structural properties of the tetrapeptide (L)V(L)P(L)A(L)L (all-L) are drastically altered upon mutating the L-proline and the L-alanine by their d-enantiomers [J. Am. Chem. Soc. 1996, 118, 6975]. The all-L diastereomer is unstructured, experiencing little or no beta-hairpin formation, while the (L)V(D)P(D)A(L)L peptide exhibits a substantial population of beta-hairpin conformation. In this study, we perform molecular dynamics simulations to investigate the folding propensity of these two model peptides. The results confirm the experimental findings, namely, that the presence of d-amino acids in the loop region strongly induces beta-hairpin formation (a population increase from about 1.5% to 50% is observed). The major factor determining the different behavior is found to be the large difference in energy between the two diastereomers, approximately 22 kJ/mol, when they adopt a beta-hairpin structure. The higher energy observed for the all-L peptide is a consequence of none-ideal hydrogen bond formation and of steric repulsions. The results suggest that selective incorporation of D-amino acids in proteins can be used to enhance certain secondary structure elements. The kinetic behavior of the folding process observed in the simulations is also investigated. We find that the decay rate of the folded structure fits to a biexponential function, suggesting that the folding/unfolding process of a beta-hairpin is governed by two different mechanisms.

Amino Acid Substitution↗
PubMed2006

Exceptional thermodynamic stability of DNA duplexes modified by nonpolar base analogues is due to increased stacking interactions and favorable solvation: Correlated ab initio calculations and molecular dynamics simulations.

The geometries of DNA hexamer (5'-GGAACC-3') and DNA 13-mer (5'-GCGTACACATGCG-3') have been determined by molecular dynamics (MD) simulations using an empirical force field. The central canonical base pair was replaced by a pair of nonpolar base analogues, 2,2'-bipyridyl and 3-methylisocarbostyril. The stabilization energy of the model system (model A) consisting of a central base pair (base-analogue pair) and two neighboring base pairs was determined by the RI-MP2 method using an extended aug-cc-pVDZ basis set. The geometry of the model was averaged from structures determined by MD simulations. The role of the solvent was covered by the COSMO continuum solvent model and calculations were performed for a larger model system (model B) which also contained a sugar-phosphate backbone. The total stabilization energies of the unperturbed system and the system perturbed by a base-analogue pair (model A) were comparable to the stability of both duplexes experimentally determined. This is due to large stacking interaction energy of the base-analogue self-pair which compensates for the missing hydrogen-bonding energy of the replaced adenine...thymine base pair. The selectivity of the base-analogue pair was reproduced (model B) when their desolvation energy was included with the interaction energy of both strands determined by the approximate SCC-DFTB-D method.

Base Pairing↗
PubMed2002

Constant normal pressure, constant surface tension, and constant temperature molecular dynamics simulation of hydrated 1,2-dilignoceroylphosphatidylcholine monolayer.

A constant normal pressure, constant surface tension, and constant temperature (NP(N)gammaT) molecular dynamics (MD) simulation of the liquid condensed phase of a 1,2-dilignoceroylphosphatidylcholine (DLGPC) monolayer has been performed at 293.15 K. A DLGPC molecule has two saturated 24-carbon acyl chains, giving the hydrocarbon core thickness of the monolayer approximately 28 A, which is close to the hydrocarbon core thickness of a membrane of a living system. NP(N)gammaT ensemble was used to reproduce the experimental observations, such as area/lipid, because surface tension is an essential factor in determining the monolayer structure. Data analysis on DLGPC/water monolayer shows that various liquid condensed-phase properties of the monolayer have been well reproduced from the simulation, indicating that surface tension 22.9 mN/M used in the simulation is an appropriate condition for the condensed-phase NP(N)gammaT simulation. The simulation results suggest that this long-chain phospholipid monolayer shares many structural characteristics with typical short-chain 1,2-diacylphosphatidylcholine systems, such as DPPC/water monolayer in the condensed phase and DPPC/water bilayer in the gel phase. Furthermore, it was found that DLGPC/water monolayer has almost completely rotationally disordered acyl chains, which have not been observed so far in short-chain 1,2-diacylphosphatidylcholine/water bilayers. This study indicates the good biological relevance of the DLGPC/water monolayer which might be useful in protein/lipid studies to reveal protein structure and protein/lipid interactions in a membrane environment.

Acylation↗
PubMed2002

Experiments and molecular-dynamics simulation of elastic waves in a plasma crystal radiated from a small dipole source.

The radiation of elastic waves from a localized source is observed experimentally in a two-dimensional plasma crystal. An initial shear stress applied by a laser forms a small dipole source. The emerging complex wave pattern is shown to consist of outgoing compressional and shear wave pulses. Subsequent structures are identified as inward-going waves due to the finite size of the source region, which reappear on the opposite side. The compressional wave forms a trailing wave train due to strong dispersion, while the nondispersive shear wave evolves into a vortex-antivortex pair on either side. The experiments are compared with a molecular-dynamics simulation.

Journal Article↗
PubMed2001

Dimer interface of transmembrane domains for neu/erbB-2 receptor dimerization and transforming activation: a model revealed by molecular dynamics simulations.

The specific point mutation Val-->Glu664 within the transmembrane domain of the neu/erbB-2 receptor is associated with increased receptor dimerization and increased receptor tyrosine kinase activity resulting in malignant transformation of cells. It is well established that Glu and residues in proximity are necessary for receptor dimerization but many studies suggest that other intramembrane constraints, not yet elucidated, are determinant for transformation. In this work, we investigated dimer models both to understand the structural role of the Glu mutation in the transmembrane domain association and to determine helix-helix contacts required for oncogenic transformation. Different types of helix-helix association based on data resulting from Cys mutational studies of the full wild receptor and spectroscopic data of transmembrane neu peptides have been explored by molecular dynamics simulations. The study leads to propose a model for the dimeric association of the transmembrane domains of the oncogenic neu receptor showing left-handed interactions of the two helices stabilized by symmetrical hydrogen bonding interactions involving the Glu side chain on one helix and the facing carbonyl of Ala661 on the second helix. Contacting residues observed in the symmetric interface explain the transforming activity or the non transforming activity of many neu mutants. Moreover the left-handed coiled coil structure is fully consistent with recent results proving the role of rotational linkage of the transmembrane domain with the kinase domain. Comparison between the predicted dimer model and those presumed from experiments strongly suggests helix flexibility in the extracellular juxtamembrane region.

Amino Acid Sequence↗
PubMed2006

Flexibility in a drug transport accessory protein: molecular dynamics simulations of MexA.

Drug resistance in gram-negative bacteria may be conferred via efflux through a tripartite complex of an inner membrane pump, an outer membrane pore, and a periplasmic adaptor protein. These are AcrB, TolC, and AcrA, respectively, in Escherichia coli. In Pseudomonas aerugonisa, their homologs are MexB, OprM, and MexA. Defining the interdomain dynamics of the adaptor protein is essential to understanding the mechanism of complex formation. Extended (25 ns) molecular dynamics simulations of MexA have been performed to determine such interdomain dynamics. Analysis of conformational drift demonstrates substantial motions of the three domains of MexA relative to one another. Principal components analysis reveals a hinge-bending motion and rotation of the alpha-helical hairpin relative to the other domains to be the two dominant motions. These two motions provide an element of considerable flexibility which is likely to be exploited in the adaptor function of MexA.

Bacterial Outer Membrane Proteins↗
PubMed2001

Molecular dynamics simulation of the calmodulin-trifluoperazine complex in aqueous solution.

Trifluoperazine (TFP) has been widely studied in relation to its mode of binding and its inactivation of calmodulin (CaM). Most studies in solution have indicated that CaM has two high-affinity binding sites for TFP. The crystal structure of the 1:4 CaM-TFP complex (CaM-4TFP) shows that three TFP molecules bind to the C-domain of CaM, and that one TFP molecule binds to the N-domain. In contrast, the crystal structure of the 1:1 CaM-TFP complex (CaM-1TFP) shows that one TFP molecule binds to the C-domain. It has been thought that the binding of one TFP molecule to the C-domain is followed by binding to the N-domain. The crystal structure of the 1:2 CaM-TFP complex (CaM-2TFP), moreover, has recently been determined, showing that two TFP molecules bind to the C-domain. In order to determine the structure of the CaM-TFP complex and to clarify the interaction between CaM and TFP in solution, we performed a molecular dynamics simulation of the CaM-TFP complex in aqueous solution starting from the CaM-4TFP crystal structure. The obtained solution structure is very similar to the CaM-2TFP crystal structure. The computer simulation showed that the binding ability of the secondary binding site of the C-domain is higher than that of the primary binding site of the N-domain.

Binding Sites↗
PubMed2006

Possible pathway for ubiquinone shuttling in Rhodospirillum rubrum revealed by molecular dynamics simulation.

In the last decade, the structures of many components of the photosynthetic apparatus of purple bacteria, as well as the mutual organization of these components within the purple membrane, were resolved. One key question that emerged concerned the assembly of the core complex consisting of the reaction center (RC) and the light-harvesting 1 (LH1) complex. In some species, like Rhodobacter sphaeroides, the ring-shaped LH1 complex was found to be open, whereas other species, like Rhodospirillum rubrum, have a closed ring surrounding the reaction center. This poses the question of how the ubiquinone molecule that transports electrons and protons from the RC to the cytochrome bc(1) complex overcomes the apparent barrier of the LH1 ring. In this study, we investigated how, in the case of a closed LH1 ring, the ubiquinone molecule diffuses through the LH1 ring. For this purpose, the LH1 structure of R. rubrum was modeled and the potential of mean force along the diffusion pathway through the LH1 was determined by steered molecular-dynamics simulations. The potential was reconstructed using the fluctuation theorem in combination with the stiff spring approximation. An upper limit for the mean first-passage time for diffusion of ubiquinone through the LH1 ring, based on a worst-case scenario potential, was calculated as approximately 8 x 10(-3) s, which is still in agreement with known turnover rates of RC and RC-LH1 complexes in the range of approximately 1000 Hz.

Biological Transport↗
PubMed2006

Molecular dynamics simulations of surface-initiated melting of nitromethane.

The melting of nitromethane initiated at solid-vacuum interfaces has been investigated using molecular dynamics nvt simulations with a realistic force field [D. C. Sorescu et al., J. Phys. Chem. B 104, 8406 (2000)]. The calculated melting point (251+/-5 K) is in good agreement with experiment (244.73 K) and values obtained previously (approximately 255.5 and 266.5+/-8 K) using other simulation methods [P. M. Agrawal et al., J. Chem. Phys. 119, 9617 (2003)]. Analyses of the molecular orientations and diffusion during the simulations as functions of the distance from the exposed surfaces show that the melting is a direct crystal-to-liquid transition, in which the molecules first gain rotational freedom, then mobility. There is a slight dependence of the melting temperature on the exposed crystallographic face.

Journal Article↗
PubMed2004

Molecular Dynamics Simulations of Micelle Formation around Dimeric Glycophorin A Transmembrane Helices.

Insertion and formation of membrane proteins involves the interaction of protein helices with one another in lipid environments. Researchers have studied glycophorin A (GpA) transmembrane helices embedded in sodium dodecyl sulfate (SDS) micelles to identify contacts significant for helix dimerization. However, a detailed picture of the conformation and dynamics of the GpA-SDS system cannot be obtained solely through experiment. Molecular dynamics simulations of SDS and a GpA dimer can provide an atomic-level picture of SDS aggregation and helix association. We report 2.5-ns simulations of GpA wild-type and mutants in a preformed micelle as well as a 32-ns simulation showing the formation of a complete micelle around wild-type GpA from an initially random placement of SDS molecules in an aqueous environment. In the latter case, an initial instability of GpA helices in water is reversed after the helices become surrounded by SDS. The properties of the spontaneously formed micelle surrounding the GpA are indistinguishable from those of the preformed micelle surrounding the GpA dimer.

Amino Acid Substitution↗