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,585 records · Page 88Linked to original sources
PubMed2006

Blocking of the nicotinic acetylcholine receptor ion channel by chlorpromazine, a noncompetitive inhibitor: A molecular dynamics simulation study.

A large series of pharmacological agents, distinct from the typical competitive antagonists, block in a noncompetitive manner the permeability response of the nicotinic acetylcholine receptor (nAChR) to the neurotransmitter acetylcholine. Taking the neuroleptic chlorpromazine (CPZ) as an example of such agents, the blocking mechanism of noncompetitive inhibitors to the ion channel pore of the nAChR has been explored at the atomic level using both conventional and steered molecular dynamics (MD) simulations. Repeated steered MD simulations have permitted calculation of the free energy (approximately 36 kJ/mol) of CPZ binding and identification of the optimal site in the region of the serine and leucine rings, at approximately 4 A from the pore entrance. Coulomb and the Lennard-Jones interactions between CPZ and the ion channel as well as the conformational fluctuations of CPZ were examined to assess the contribution of each to the binding of CPZ to the nAChR. The MD simulations disclose a dynamic interaction of CPZ binding to the nAChR ionic channel. The cationic ammonium head of CPZ forms strong hydrogen bonds with Glu262 (alpha), Asp268 (beta), Glu272 (beta), Ser276 (beta), Glu280 (delta), Gln271 (gamma), Glu275 (gamma), and Asn279 (gamma) nAChR residues. Finally, the conventional MD simulation of CPZ at its identified binding site demonstrates that the binding of CPZ not only blocks ion transport through the channel but also markedly inhibits the conformational transitions of the channel, necessary for nAChR to carry out its biological function.

Acetylcholine↗
PubMed2005

Anionic clays containing anti-inflammatory drug molecules: comparison of molecular dynamics simulation and measurements.

Three representative nonsteroidal anti-inflammatory drug molecules, Ibuprofen, Diclofenac, and Indomethacin, have been intercalated within the galleries of an anionic clay, Mg-Al layered double hydroxide (LDH). X-ray diffraction, IR and Raman vibrational spectroscopy and (13)C cross-polarization magic-angle spinning NMR have been used to characterize the confined drug molecules, while molecular dynamics (MD) simulations were used to probe the interlayer structure, arrangement, orientation, and geometry of the intercalated species. All three drug molecules are arranged as bilayers in the interlamellar space of the anionic clay. But while the structure of the intercalated Ibuprofen is identical to that of the molecule outside the layers, spectroscopy as well as MD simulation shows that there is a change in the geometry of Diclofenac and Indomethacin upon confinement within the galleries of the LDH. The change in geometry of Diclofenac and Indomethacin upon intercalation is shown to originate from the electrostatic interaction between the electronegative chlorine atoms on the drug molecule and the positively charged metal hydroxide sheets of the anionic clay. It is shown that these changes in the geometry of the intercalated drug molecules allow for the observed interlayer spacing to be realized without the bilayers having to interdigitate, which would otherwise have been necessary if the structure of the drug molecules had remained identical to that outside the layers. Comparisons of experimental measurements with simulation have provided a more detailed understanding of the geometry and organization of flexible drug molecules confined in the anionic clay.

Journal Article↗
PubMed2003

Myoglobin-CO substate structures and dynamics: multidimensional vibrational echoes and molecular dynamics simulations.

Spectrally resolved infrared stimulated vibrational echo data were obtained for sperm whale carbonmonoxymyoglobin (MbCO) at 300 K. The measured dephasing dynamics of the CO ligand are in agreement with dephasing dynamics calculated with molecular dynamics (MD) simulations for MbCO with the residue histidine-64 (His64) having its imidazole epsilon nitrogen protonated (N(epsilon)-H). The two conformational substate structures B(epsilon) and R(epsilon) observed in the MD simulations are assigned to the spectroscopic A(1) and A(3) conformational substates of MbCO, respectively, based on the agreement between the experimentally measured and calculated dephasing dynamics for these substates. In the A(1) substate, the N(epsilon)-H proton and N(delta) of His64 are approximately equidistant from the CO ligand, while in the A(3) substate, the N(epsilon)-H of His64 is oriented toward the CO, and the N(delta) is on the surface of the protein. The MD simulations show that dynamics of His64 represent the major source of vibrational dephasing of the CO ligand in the A(3) state on both femtosecond and picosecond time scales. Dephasing in the A(1) state is controlled by His64 on femtosecond time scales, and by the rest of the protein and the water solvent on longer time scales.

Animals↗
PubMed1995

Interaction of an amphiphilic peptide with a phospholipid bilayer surface by molecular dynamics simulation study.

Corticotropin-releasing factor (CRF) is the principal neuroregulator of adrenocorticotropic hormone (ACTH) secretion. Previous experiments have demonstrated that CRF binds avidly to the surface of single egg phosphatidylcholine vesicles and its amphiphilic secondary structure might play an important role in the function. In this study, the interaction of the residues 13-41 in human CRF with the surface of a DOPC bilayer was investigated by molecular dynamics (MD) simulation in order to understand the role of the membrane surface in the formation of the amphiphilic alpha helix as well as to determine the effects of the peptide on the lipid bilayer. The model used included 60 DOPC molecules, 1 helical peptide (CRF13-41) on the bilayer surface, and explicit waters of solvation in the lipid polar head group regions, together with constant-volume periodic boundary conditions in three dimensions. The MD simulation was carried out for 510 ps. In addition, CRF13-41, initially in a helical form, was simulated in vacuo as a control. The results indicate that while it was completely unstable in vacuo, the peptide helical form was generally maintained on the bilayer surface, but with distortions near the terminal ends. The peptide was confined to the bilayer headgroup/water region, similar to that reported from neutron diffraction measurement of tripeptides bound to the phosphatidylcholine bilayer surface (Ref 1). The amphiphilicity of the peptide matched that of the bilayer headgroup environment, with the hydrophilic side oriented toward water and the hydrophobic side making contact with the bilayer hydrocarbon core. These results support the hypothesis that the amphiphilic environment of a membrane surface is important in the induction of peptide amphiphilic alpha-helical secondary structure. Two major effects of the peptide on the lipids were found: the first CH2 segment in the lipid chains was significantly disordered and the lipid headgroup distribution was broadened towards the water region.

Amino Acid Sequence↗
PubMed2001

Finite-size effects in molecular dynamics simulations: intermediate scattering function and velocity of sound. III. Theory and application to a model krypton fluid.

We describe a method for obtaining the intermediate scattering function I(Q,t) from a computer simulation: it is an extension of our earlier calculation [Salacuse, Denton, and Egelstaff, Phys. Rev. E 53, 2382 (1996)] for the t-->0 limit. We use this approach to obtain I(Q,t) for low Q and t from molecular dynamics (MD) simulations of a model krypton fluid whose atoms interact via a truncated Aziz pair potential, and the results are compared over their range of validity to I(Q,t) determined by the standard MD method and also by a time expansion approach. In its range of validity our approach is much more efficient than the standard MD method; however, it covers a restricted range of t due to the movement of density fluctuations (sound waves) through the simulated fluid which produces an anomaly in the time behavior of I(Q,t). By analyzing I(Q=0,t) the velocity of sound in the simulation is determined, and the results compare favorably with published experimental results for the sound velocity of liquid krypton.

Journal Article↗
PubMed2002

Molecular dynamics simulation of proton transport through the influenza A virus M2 channel.

The structural and dynamical properties of a solvated proton in the influenza A virus M2 channel are studied using a molecular dynamics (MD) simulation technique. The second-generation multi-state empirical valence bond (MS-EVB2) model was used to describe the interaction between the excess proton and the channel environment. Solvation structures of the excess proton and its mobility characteristics along the channel were determined. It was found that the excess proton is capable of crossing the channel gate formed by the ring of four histidine residues even though the gate was only partially open. Although the hydronium ion itself did not cross the channel gate by traditional diffusion, the excess proton was able to transport through the ring of histidine residues by hopping between two water molecules located at the opposite sides of the gate. Our data also indicate that the proton diffusion through the channel may be correlated with the changes in channel conformations. To validate this observation, a separate simulation of the proton in a "frozen" channel has been conducted, which showed that the proton mobility becomes inhibited.

Biological Transport↗
PubMed2001

Calculation of NMR-relaxation parameters for flexible molecules from molecular dynamics simulations.

Comparatively small molecules such as peptides can show a high internal mobility with transitions between several conformational minima and sometimes coupling between rotational and internal degrees of freedom. In those cases the interpretation of NMR relaxation data is difficult and the use of standard methods for structure determination is questionable. On the other hand, in the case of those system sizes, the timescale of both rotational and internal motions is accessible by molecular dynamics (MD) simulations using explicit solvent. Thus a comparison of distance averages ([r(-6)](-1/6) or [r(-3)](1/3)) over the MD trajectory with NOE (or ROE) derived distances is no longer necessary, the (back)calculation of the complete spectra becomes possible. In the present study we use two 200 ns trajectories of a heptapeptide of beta-amino acids in methanol at two different temperatures to obtain theoretical ROESY spectra by calculating the exact spectral densities for the interproton vectors and the full relaxation matrix. Those data are then compared with the experimental ones. This analysis permits to test some of the assumptions and approximations that generally have to be made to interpret NMR spectra, and to make a more reliable prediction of the conformational equilibrium that leads to the experimental spectrum.

Computer Simulation↗
PubMed2005

Dynamic mechanism for the autophosphorylation of CheA histidine kinase: molecular dynamics simulations.

The two-component system (TCS) is an important signal transduction component for most bacteria. This signaling pathway is mediated by histidine kinases via autophosphorylation between P1 and P4 domains. Taking chemotaxis protein CheA as a model of TCS, the autophosphorylation mechanism of the TCS histidine kinases has been investigated in this study by using a computational approach integrated homology modeling, ligand-protein docking, protein-protein docking, and molecular dynamics (MD) simulations. Four nanosecond-scale MD simulations were performed on the free P4 domain, P4-ATP, P4-TNPATP, and P1-P4-ATP complexes, respectively. Upon its binding to the binding pocket of P4 with a folded conformation, ATP gradually extends to an open state with help from a water molecule. Meanwhile, ATP forms two hydrogen bonds with His413 and Lys494 at this state. Because of the lower energy of the folded conformations, ATP shrinks back to its folded conformations, leading to the rupture of the hydrogen bond between ATP and Lys494. Consequently, Lys494 moves away from the pocket entrance, resulting in an open of the ATP lid of P4. It is the open state of P4 that can bind tightly to P1, where the His45 of P1 occupies a favorable position for its autophosphorylation from ATP. This indicates that ATP is not only a phosphoryl group donor but also an activator for CheA phosphorylation. Accordingly, a mechanism of the autophosphorylation of CheA is proposed as that the ATP conformational switch triggers the opening of the ATP lid of P4, leading to P1 binding tightly, and subsequently autophosphorylation from ATP to P1.

Adenosine Triphosphate↗
PubMed2000

Molecular dynamics simulations highlight mobile regions in proteins: A novel suggestion for converting a murine V(H) domain into a more tractable species.

The V(H) region of the murine antibody 1F7 has been identified as a single-domain chorismate mutase, but a tendency to denature and aggregate has hampered its biochemical characterization. Standard mutagenesis approaches targeting antibody chain dimerization areas have been exhausted. We describe a new approach to the problem, where we use molecular dynamics (MD) simulations to find the differences between the untractable protein and the known soluble V(H) domain from a llama antibody. MD simulations of proteins yield information on the relative stability and fluctuations of parts of the proteins. By comparing simulation results of two related proteins their differences in stability and fluctuations can be analyzed and may suggest mutations aimed at (de)stabilization of one of the two proteins. For the mouse versus llama simulations, this approach highlights an untried area in the protein which shows increased fluctuations. The replacement of this eight-residue segment with the corresponding llama sequence gave a chimeric mutant which shows significantly decreased fluctuations. We see this as a general scheme to generate suggestions for mutagenesis experiments, not only obviously generalizable to other immunoglobulin domains, but to other protein systems as well.

Amino Acid Sequence↗
PubMed2005

Slip behavior in liquid films on surfaces of patterned wettability: comparison between continuum and molecular dynamics simulations.

We investigate the behavior of the slip length in Newtonian liquids subject to planar shear bounded by substrates with mixed boundary conditions. The upper wall, consisting of a homogenous surface of finite or vanishing slip, moves at a constant speed parallel to a lower stationary wall, whose surface is patterned with an array of stripes representing alternating regions of no shear and finite or no slip. Velocity fields and effective slip lengths are computed both from molecular dynamics (MD) simulations and solution of the Stokes equation for flow configurations either parallel or perpendicular to the stripes. Excellent agreement between the hydrodynamic and MD results is obtained when the normalized width of the slip regions, a/sigma greater than or approximately equal O (10) , where sigma is the (fluid) molecular diameter characterizing the Lennard-Jones interaction. In this regime, the effective slip length increases monotonically with a/sigma to a saturation value. For a/sigma less than or approximately O (10) and transverse flow configurations, the nonuniform interaction potential at the lower wall constitutes a rough surface whose molecular scale corrugations strongly reduce the effective slip length below the hydrodynamic results. The translational symmetry for longitudinal flow eliminates the influence of molecular scale roughness; however, the reduced molecular ordering above the wetting regions of finite slip for small values of a/sigma increases the value of the effective slip length far above the hydrodynamic predictions. The strong correlation between the effective slip length and the liquid structure factor representative of the first fluid layer near the patterned wall illustrates the influence of molecular ordering effects on slip in noninertial flows.

Journal Article↗
PubMed2005

Cyan fluorescent protein: molecular dynamics, simulations, and electronic absorption spectrum.

The dynamics and electronic absorption spectrum of enhanced cyan fluorescent protein (ECFP), a mutant of green fluorescent protein (GFP), have been studied by means of a 1 ns molecular dynamics (MD) simulation. The two X-ray conformations A' and B' of ECFP were considered. The chromophore was assumed to be neutral, and all titratable residues were taken in their standard protonation state at neutral pH. The protein was embedded in a box of water molecules (and counterions). The first result is that the two conformations A' and B' are found to be stable all along the simulation. Then, an analysis of the hydrogen-bond networks shows strong differences between the two conformations in the surroundings of the nitrogen atom of the indolic part of the chromophore. This is partly due to the imperfection in the beta barrel near the His148 residue, which allows the access of one solvent molecule inside the protein in conformation A'. Finally, quantum mechanical calculations of the electronic transition energies of the chromophore in the charge cloud of the protein and solvent water molecules were performed using the TDDFT method on 160 snapshots extracted every 5 ps of the MD trajectories. It is found that conformations A' and B' exhibit very similar spectra despite different H-bond networks involving the chromophore. This similarity is related to the weak charge transfer involved in the electronic transition and the weak electrostatic field created by ECFP near the chromophore, within the hypotheses made in the present simulation.

Computer Simulation↗
PubMed2006

Molecular dynamics simulations of peptide carboxylate hydration.

Aqueous solvation of carboxylate groups, as present in the glycine zwitterion and the dipeptide aspartylalanine, is studied employing a force-field that includes distributed multipole electrostatics and induction contributions (Amoebapro: P. Ren and J. W. Ponder, J. Comput. Chem., 2002, 23, 1497; P. Ren and J. W. Ponder, J. Phys. Chem. B, 2003, 107, 5933; J. W. Ponder and D. A. Case, Adv. Protein Chem., 2003, 66, 27). Radial and orientation distribution functions, as well as hydration numbers, are calculated and compared with existing simulation data derived from Car-Parrinello molecular dynamics (CPMD), and also distributed-charge force-fields. Connections are also made with experimental data for solvation of carboxylates in water. Our findings show that Amoebapro yields carboxylate solvation properties in very good agreement with CPMD results, significantly closer agreement than can be obtained from traditional force-fields. We also demonstrate that the influence of solvation on the conformation of the dipeptide is markedly different using Amoebapro compared with the other force-fields.

Carboxylic Acids↗
PubMed2005

Molecular dynamics simulations and free energy calculations of base flipping in dsRNA.

The family of adenosine deaminases acting on RNA (ADARs) targets adenosines in RNA that is mainly double stranded. Some substrates are promiscuously deaminated whereas others, such as the mammalian glutamate receptor B (gluR-B) pre-mRNA, are more selectively deaminated. Many DNA/RNA-base modification enzymes use a base flipping mechanism to be able to reach their target base and it is believed that ADARs function in a similar way. In this study we used molecular dynamics (MD) simulations to describe two sites on the gluR-B pre-mRNA, the selectively targeted R/G site and the nontargeted 46 site, in an attempt to explain the substrate specificity. We used regular MD and also a forced base flipping method with umbrella sampling to calculate the free energy of base opening. Spontaneous opening of the mismatched adenosine was observed for the R/G site but not for the 46 site.

Adenosine↗
PubMed2006

Structural dynamics of the mitochondrial ADP/ATP carrier revealed by molecular dynamics simulation studies.

The mitochondrial adenosine diphosphate/adenosine triphosphate (ADP/ATP) carrier has been recently crystallized in complex with its specific inhibitor carboxyatractyloside (CATR). In the crystal structure, the six-transmembrane helix bundle that defines the nucleotide translocation pathway is closed on the matrix side due to sharp kinks in the odd-numbered helices. The closed conformation is further sealed by the loops protruding into the matrix that interact through an intricate network of charge-pairs. To gain insight into its structural dynamics we performed molecular dynamics (MD) simulation studies of the ADP/ATP carrier with and without its cocrystallized inhibitor. The two trajectories sampled a conformational space around two different configurations characterized by distinct salt-bridge networks with a significant shift from inter- to intrarepeat bonding on the matrix side in the absence of CATR. Analysis of the geometrical parameters defining the transmembrane helices showed that even-numbered helices can undergo a face rotation, whereas odd-numbered helices can undergo a change in the wobble angle with a conserved proline acting as molecular hinge. Our results provide new information on the dynamical properties of the ADP/ATP carrier and for the first time yield a detailed picture of a stable carrier conformation in absence of the inhibitor.

Animals↗
PubMed1999

Incoherent neutron scattering of copper azurin: a comparison with molecular dynamics simulation results.

The low-frequency dynamics of copper azurin has been studied at different temperatures for a dry and deuterium hydrated sample by incoherent neutron scattering and the experimental results have been compared with molecular dynamics (MD) simulations carried out in the same temperature range. Experimental Debye-Waller factors are consistent with a dynamical transition at approximately 200 K which appears partially suppressed in the dry sample. Inelastic and quasielastic scattering indicate that hydration water modulates both vibrational and diffusive motions. The low-temperature experimental dynamical structure factor of the hydrated protein shows an excess of inelastic scattering peaking at about 3 meV and whose position is slightly shifted downwards in the dry sample. Such an excess is reminiscent of the "boson peak" observed in glass-like materials. This vibrational peak is quite well reproduced by MD simulations, although at a lower energy. The experimental quasielastic scattering of the two samples at 300 K shows a two-step relaxation behaviour with similar characteristic times, while the corresponding intensities differ only by a scale factor. Also, MD simulations confirm the two-step diffusive trend, but the slow process seems to be characterized by a decay faster than the experimental one. Comparison with incoherent neutron scattering studies carried out on proteins having different structure indicates that globular proteins display common elastic, quasielastic and inelastic features, with an almost similar hydration dependence, irrespective of their secondary and tertiary structure.

Azurin↗
PubMed2005

Molecular dynamics simulation of clustered DNA damage sites containing 8-oxoguanine and abasic site.

Clustered DNA damage sites induced by ionizing radiation have been suggested to have serious consequences to organisms, such as cancer, due to their reduced probability to be repaired by the enzymatic repair machinery of the cell. Although experimental results have revealed that clustered DNA damage sites effectively retard the efficient function of repair enzymes, it remains unclear as to what particular factors influence this retardation. In this study, approaches based on molecular dynamics (MD) simulation have been applied to examine conformational changes and energetic properties of DNA molecules containing clustered damage sites consisting of two lesioned sites, namely 7,8-dihydro-8-oxoguanine (8-oxoG) and apurinic/apyrimidinic (AP) site, located within a few base pairs of each other. After 1 ns of MD simulation, one of the six DNA molecules containing a clustered damage site develops specific characteristic features: sharp bending at the lesioned site and weakening or complete loss of electrostatic interaction energy between 8-oxoG and bases located on the complementary strand. From these results it is suggested that these changes would make it difficult for the repair enzyme to bind to the lesions within the clustered damage site and thereby result in a reduction of its repair capacity.

Algorithms↗
PubMed2007

Coarse-grained protein molecular dynamics simulations.

A limiting factor in biological science is the time-scale gap between experimental and computational trajectories. At this point, all-atom explicit solvent molecular dynamics (MD) are clearly too expensive to explore long-range protein motions and extract accurate thermodynamics of proteins in isolated or multimeric forms. To reach the appropriate time scale, we must then resort to coarse graining. Here we couple the coarse-grained OPEP model, which has already been used with activated methods, to MD simulations. Two test cases are studied: the stability of three proteins around their experimental structures and the aggregation mechanisms of the Alzheimer's Abeta16-22 peptides. We find that coarse-grained isolated proteins are stable at room temperature within 50 ns time scale. Based on two 220 ns trajectories starting from disordered chains, we find that four Abeta16-22 peptides can form a three-stranded beta sheet. We also demonstrate that the reptation move of one chain over the others, first observed using the activation-relaxation technique, is a kinetically important mechanism during aggregation. These results show that MD-OPEP is a particularly appropriate tool to study qualitatively the dynamics of long biological processes and the thermodynamics of molecular assemblies.

Amyloid beta-Peptides↗
PubMed1999

Evaluation of Interaction Forces between Macroparticles in Simple Fluids by Molecular Dynamics Simulation.

The present article provides the description of the solvation forces between large spheres in a fluid. The molecular dynamics (MD) method was applied to the relatively simple systems in which a pair of structureless macroparticles, either solvophobic or solvophilic, is immersed in a simple fluid of two types, either a soft-sphere or a Lennard-Jones fluid. When a pair of solvophobic macroparticles was in the attractive Lennard-Jones fluid, no dense layer of the solvent particles formed near the surface of the macroparticles and the strong attractive forces were induced between them. In the other combinations of macroparticles and fluids, the dense layers formed and the solvation forces oscillated, exhibiting the attraction and repulsion, whose periodic distance was about the diameter of solvent particles. Our results agreed well with those of the other simulation and theoretical studies with respect to the solvent density profile near a macroparticle and the force-distance profile between macroparticles. The benefit of our approach would be the simplicity in specifying or finding the bulk condition that is in equilibrium with the thin film of molecules between large surfaces. The present method can be applied straightforward to macroparticles immersed in mixtures and complex fluids described by the bead-spring model, to which the conventional grand canonical ensemble Monte Carlo (GCEMC) method is hardly accessible. Copyright 1999 Academic Press.

Journal Article↗