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 181 records · Page 10Linked to original sources
PubMed2005

Effect of hexafluoroisopropanol alcohol on the structure of melittin: a molecular dynamics simulation study.

The molecular mechanism by which HFIP stabilizes the alpha-helical structure of peptides is not well understood. In the present study, we use melittin as a model to gain insight into the details of the atomistic interactions of HFIP with the peptide. We have performed extensive comparative molecular dynamics simulations (up to 100 nsec) in the absence and in the presence of HFIP. In agreement with recent NMR experiments, the simulations show rapid loss of tertiary structure in water at pH 2 but much higher helicity in 35% HFIP. The MD simulations also indicate that melittin adopts a highly dynamic global structure in 35% HFIP solution with two alpha-helical segments sampling a wide range of angular orientations. The analysis of the HFIP distribution shows the tendency of HFIP to aggregate around the peptide, increasing the local cosolvent concentration to more than two times that in the bulk concentration. The correlation of local peptide structure with HFIP coating suggests that displacement of water at the peptide surface is the main contribution of HFIP in stabilizing the secondary structure of melittin. Finally, a stabilizing effect promoted by the presence of counter-ions was also observed in the simulations.

Computer Simulation↗
PubMed1997

Insight into the stabilization of A-DNA by specific ion association: spontaneous B-DNA to A-DNA transitions observed in molecular dynamics simulations of d[ACCCGCGGGT]2 in the presence of hexaamminecobalt(III).

BACKGROUND: Duplex DNA is more than a simple information carrier. The sequence-dependent structure and its inherent deformability, in concert with the subtle modulating effects of the environment, play a crucial role in the regulation and packaging of DNA. Recent advances in force field and simulation methodologies allow molecular dynamics simulations to now represent the specific effects of the environment. An understanding of the environmental dependence of DNA structure gives insight into how histones are able to package DNA, how various proteins are able to bind and modulate nucleic acid structure and will ultimately aid the design of molecules to package DNA for more effective gene therapy. RESULTS: Molecular dynamics simulations of d[ACCCGCGGGT]2 in solution in the presence of hexaamminecobalt(III) [Co(NH3)6(3+)] show stabilization of A-DNA and spontaneous B-DNA to A-DNA transitions, which is consistent with experimental results from NMR and Raman spectroscopic and X-ray crystallographic studies. In the absence of Co(NH3)6(3+), A-DNA to B-DNA transitions are observed instead. In addition to their interaction with the guanines in the major groove, Co(NH3)6(3+) ions bridge opposing strands in the bend across the major groove, probably stabilizing A-DNA. CONCLUSIONS: The simulation methods and force fields have advanced to a sufficient level that some representation of the environment can be seen in nanosecond length molecular dynamics simulations. These simulations suggest that, in addition to the general explanation of A-DNA stabilization by dehydration, hydration and ion association in the major groove stabilize A-DNA.

Cobalt↗
PubMed2005

Force probe molecular dynamics simulations.

Many proteins are molecular nano-machines, which perform their biological function via well-coordinated structural transitions. Often, these motions occur on much slower time scales than those accessible to conventional molecular dynamics techniques, which are limited to submicrosecond time scales by current computer technology. This is also true for ligand binding and unbinding reactions. Force probe simulations (or steered molecular dynamics) provide a powerful means to overcome this limitation, and thus to get insight into the atomistic mechanisms that underlie biological functions such as ligand binding. This chapter provides a basic introduction into this method. It further sketches a simple nonequilibrium statistical mechanics treatment that shows how to relate the results of force probe simulations to atomic force microscopy (AFM) or optical tweezer experiments. As an example, enforced unbinding simulations of streptavidin/biotin complexes are detailed.

Biomechanical Phenomena↗
PubMed1993

Molecular dynamics simulation of Pf1 coat protein.

The results of molecular dynamics simulations of Pf1 coat protein are described and compared to experimental NMR data on both the membrane bound and structural forms of this viral coat protein. Molecular dynamics simulations of the 46 residue coat protein and related model sequences were performed according to a simple protocol. The simulations were initiated with the polypeptides in a completely uniform alpha helical conformation in a dielectric continuum (epsilon = 2) and the motions of individual residues were followed as a function of time by monitoring the angular fluctuations of amide NH bond vectors. The simulations of Pf1 coat protein were able to identify the same mobile and structured segments found in experimental NMR studies of the membrane bound form of the protein (Shon, K.-J., Y. Kim, L. A. Colnago, and S. J. Opella. 1991. Science (Wash. DC). 252:1303-1305). Significantly, in addition to mobile amino and carboxyl terminal regions, a mobile internal loop was found that connects the rigid hydrophobic and amphipathic helices in the protein. NMR experiments show that this mobile loop is present in both the viral and membrane bound forms of the protein and that it plays a role in viral assembly (Nambudripad, R., W. Stark, S. J. Opella, and L. Makowski. 1991. Science (Wash. DC) 252:1305-1308). The results of simulations of several alanine based 46 residue polypeptides with some of the charged residues present in the Pf1 coat protein sequence suggest that interactions between the Asp 14 and Asp 18 sidechains and the peptide backbone are responsible for the formation of the mobile loop. The agreement between the results of the calculations presented here and the previously reported NMR experiments suggest that molecular dynamics simulations might be useful in the prediction of the secondary structure and dynamics of individual residues in membrane and structural proteins with predominantly alpha helical secondary structure.

Bacteriophages↗
PubMed1994

Molecular dynamics simulations of HIV-1 protease with peptide substrate.

Molecular dynamics simulations of human immunodeficiency virus (HIV)-1 protease with a model substrate were used to test if there is a stable energy minimum for a proton that is equidistant from the four delta oxygen atoms of the two catalytic aspartic acids. The crystal structure of HIV-1 protease with a peptidic inhibitor was modified to model the peptide substrate Ser-Gln-Asn-Tyr-Pro-Ile-Val-Gln for the starting geometry. A proton was positioned between the two closet oxygen atoms of the two catalytic aspartic acids, and close to the carbonyl oxygen of the scissile bond in the substrate. All crystallographic water molecules were included. Two molecular dynamics simulations were run: 30 ps with united-atom potentials and 40 ps using the more accurate all-atom potentials. The molecular dynamics used a new algorithm that increased the speed and allowed the elimination of a cut-off for non-bonded interactions and the inclusion of an 8 A shell of water molecules in the calculations. The overall structure of the protease dimer, including the catalytic aspartic acids, was stable during the course of the molecular dynamics simulations. The substrate and a water molecule, that is an important component of the binding site, were stable during the simulation using all-atom potentials, but more mobile when united-atom potentials were used. A Poincare map representation showed that the positions of the proton and its coordinating oxygen atoms were stable for 93% of both simulations, although many of the buried and poorly accessible water molecules exchanged with solvent. The proton has a stable minimum energy position and maintains coordination with all four delta oxygen atoms of the two catalytic aspartic acids and the carbonyl oxygen of the scissile bond of the substrate. Therefore, a loosely bound hydrogen ion at this position will not be rapidly exchanged with solvent, and will rebond to either a catalytic aspartic acid or possibly the substrate. The implications for the reaction mechanism are discussed.

Amino Acid Sequence↗
PubMed2005

Coupling between lysozyme and glycerol dynamics: microscopic insights from molecular-dynamics simulations.

We explore possible molecular mechanisms behind the coupling of protein and solvent dynamics using atomistic molecular-dynamics simulations. For this purpose, we analyze the model protein lysozyme in glycerol, a well-known protein-preserving agent. We find that the dynamics of the hydrogen bond network between the solvent molecules in the first shell and the surface residues of the protein controls the structural relaxation (dynamics) of the whole protein. Specifically, we find a power-law relationship between the relaxation time of the aforementioned hydrogen bond network and the structural relaxation time of the protein obtained from the incoherent intermediate scattering function. We demonstrate that the relationship between the dynamics of the hydrogen bonds and the dynamics of the protein appears also in the dynamic transition temperature of the protein. A study of the dynamics of glycerol as a function of the distance from the surface of the protein indicates that the viscosity seen by the protein is not the one of the bulk solvent. The presence of the protein suppresses the dynamics of the surrounding solvent. This implies that the protein sees an effective viscosity higher than the one of the bulk solvent. We also found significant differences in the dynamics of surface and core residues of the protein. The former is found to follow the dynamics of the solvent more closely than the latter. These results allowed us to propose a molecular mechanism for the coupling of the solvent-protein dynamics.

Algorithms↗
PubMed1995

Conformational modeling of elastin tetrapeptide Boc-Gly-Leu-Gly-Gly-NMe by molecular dynamics simulations with improvements to the thermalization procedure.

Molecular Dynamics Simulations (MD) at Constant-Temperature or Constant-total Energy for the conformational Global-Minimum (GM) of elastin tetrapeptide Boc-Gly-Leu-Gly-Gly-NMe have been performed. The thermalization problem concerning the initial state of Constant-Temperature MD has been solved developing two effective strategies. In the first one, the run starts from the room-temperature state reached by Molecular Dynamics Simulated Annealing (SA). In the second one, one starts from the annealed-state at low-temperature and performs a long constant-low-temperature run until the initial conformer is perfectly equilibrated. Then, the low-temperature equilibrated-state is used as initial state for MD at room-temperature. heuristic criteria on order to define the onset of steady-state have been established monitoring the hystories of collective parameters (e.g., the total energy, temperature, end-to-end distance, etc.) and their amplitude fluctuations. Moreover, the equilibrium between the system and the heat bath is verified analyzing the total linear momentum conservation by the time evolution of center mass velocity. The slow drift of total energy during Constant-total Energy MD has been corrected using a loose coupling between the system and the heat bath. Moreover, we have verified that the roto-translational motions do not affect significantly the properties of molecular vibrations. The librations of peptide unit inside the type II beta-turn [Gly1]C = 0 ... HN[Gly4], previously detected, were confirmed. Large -Gly-Gly- chain motions were identified and modeled as fluctuations occurring between the tetratepeptide GM and the saddle-point corresponding to the transition state of the conversion toward the extended-chain conformation. All these peptide motions could contribute to the elasticity mechanism of elastin.

Amino Acid Sequence↗
PubMed2006

Model of an asymmetric DPPC/DPPS membrane: effect of asymmetry on the lipid properties. A molecular dynamics simulation study.

The study of asymmetric lipid bilayers is of a crucial importance due to the great number of biological process in which they are involved such as exocytosis, intracellular fusion processes, phospholipid-protein interactions, and signal transduction pathway. In addition, the loss of this asymmetry is a hallmark of the early stages of apoptosis. In this regard, a model of an asymmetric lipid bilayer composed of DPPC and DPPS was simulated by molecular dynamics simulation. Thus, the asymmetric membrane was modeled by 264 lipids, of which 48 corresponded to DPPS- randomly distributed in the same leaflet with 96 DPPC. In the other leaflet, 120 DPPC were placed without DPPS-. Due to the presence of a net charge of -1 for the DPPS- in physiological conditions, 48 Na+ were introduced into the system to balance the charge. To ascertain whether the presence of the DPPS- in only one of the two leaflets perturbs the properties of the DPPC in the other leaflet composed only of DPPC, different properties were studied, such as the atomic density of the different components across the membrane, the electrostatic potential across the membrane, the translational diffusion of DPPC and DPPS, the deuterium order parameters, lipid hydration, and lipid-lipid charge bridges. Thus, we obtained that certain properties such as the surface area lipid molecule, lipid head orientation, order parameter, translational diffusion coefficient, or lipid hydration of DPPC in the leaflet without DPPS remain unperturbed by the presence of DPPS in the other leaflet, compared with a DPPC bilayer. On the other hand, in the leaflet containing DPPS, some of the DPPC properties were strongly affected by the presence of DPPS such as the order parameter or electrostatic potential.

1,2-Dipalmitoylphosphatidylcholine↗
PubMed2002

Ultrafast spectroscopy reveals subnanosecond peptide conformational dynamics and validates molecular dynamics simulation.

Femtosecond time-resolved spectroscopy on model peptides with built-in light switches combined with computer simulation of light-triggered motions offers an attractive integrated approach toward the understanding of peptide conformational dynamics. It was applied to monitor the light-induced relaxation dynamics occurring on subnanosecond time scales in a peptide that was backbone-cyclized with an azobenzene derivative as optical switch and spectroscopic probe. The femtosecond spectra permit the clear distinguishing and characterization of the subpicosecond photoisomerization of the chromophore, the subsequent dissipation of vibrational energy, and the subnanosecond conformational relaxation of the peptide. The photochemical cis/trans-isomerization of the chromophore and the resulting peptide relaxations have been simulated with molecular dynamics calculations. The calculated reaction kinetics, as monitored by the energy content of the peptide, were found to match the spectroscopic data. Thus we verify that all-atom molecular dynamics simulations can quantitatively describe the subnanosecond conformational dynamics of peptides, strengthening confidence in corresponding predictions for longer time scales.

Computer Simulation↗
PubMed2026

Validating the potential mechanism and therapeutic effect of Qinlian Jiangxia decoction in the treatment of type 2 diabetes mellitus complicated with hyperlipidemia through network pharmacology, molecular docking, molecular dynamics simulation, andexperiments.

OBJECTIVE: To investigate the mechanism of action of Qinlian Jiangxia decoction (, QLJXD) in the treatment of type 2 diabetes mellitus (T2DM) complicated by hyperlipidemia using network pharmacology, molecular docking, molecular dynamics simulation and in vivo experiments. METHODS: Drug components, targets and disease targets were identified using databases such as TCM systems pharmacology database and analysis platform and GeneCards. The intersecting targets were subjected to protein-protein interaction analysis using the search tool for the retrieval of interacting genes/proteins database. Subsequently, Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analysis of the intersecting targets were conducted using the Metascape platform to identify core components and targets. The results were validated using molecular docking, molecular dynamics simulations and in vivo experiments. RESULTS: QLJXD contains 76 active ingredients and 136 disease targets. The core ingredients are quercetin, β-sitosterol, wogonin and baicalein, while the core targets are fatty acid binding protein 4 (FABP4) and peroxisome proliferative activated receptor gamma (PPARG). Molecular docking and molecular dynamics simulations revealed that the core ingredients bound well to the core targets. Animal experiments demonstrated that QLJXD effectively inhibited the expression of FABP4 and increased the expression of PPARG, thereby enhancing disorders of glycolipid metabolism. CONCLUSION: The putative therapeutic efficacy of QLJXD in the management of T2DM complicated with hyperlipidemia may be ascribed to the synergistic actions of multiple components, such as quercetin, β-sitosterol, wogonin, and baicalein, which collectively modulate FABP4 and PPARG molecular targets.

Molecular Docking Simulation↗
PubMed1992

Molecular dynamics simulations of the whey protein beta-lactoglobulin.

Molecular dynamics simulations have been used to model the motions and conformational behavior of the whey protein bovine beta-lactoglobulin. Simulations were performed for the protein by itself and complexed to a single retinol ligand located in a putative interior binding pocket. In the absence of the retinol ligand, the backbone loops around the opening of this interior pocket shifted inward to partially close off this cavity, similar to the shifts observed in previously reported molecular dynamics simulations of the uncomplexed form of the homologous retinol binding protein. The protein complexed with retinol does not exhibit the same conformational shifts. Conformational changes of this type could serve as a recognition signal allowing in vivo discrimination between the free and retinol complexed forms of the beta-lactoglobulin molecule. The unusual bending of the single alpha-helix observed in the simulations of retinol binding protein were not observed in the present calculations.

Animals↗
PubMed2002

The Trp cage: folding kinetics and unfolded state topology via molecular dynamics simulations.

Using over 75 mus of molecular dynamics simulation, we have generated several thousand folding simulations of the 20-residue Trp cage at experimental temperature and solvent viscosity. A total of 116 independent folding simulations reach RMSDcalpha values below 3 A RMSDcalpha, some as close as 1.4 A RMSDcalpha. We estimate a folding time of 5.5+/-3.5 mus, a rate that is in reasonable agreement with experimental kinetics. Finally, we characterize both the folded and unfolded ensemble under native conditions and note that the average topology of the unfolded ensemble is very similar to the topology of the native state.

Amino Acid Sequence↗
PubMed2005

Partition of protein solvation into group contributions from molecular dynamics simulations.

Linear response theory coupled to molecular dynamics simulations with an explicit solvent representation is used to derive fractional contributions of amino acid residues to the solvation of proteins. The new fractional methods developed here are compared with standard approaches based on empirical 1D and 3D statistical potentials, as well as with estimates obtained from the analysis of classical molecular interaction potentials. The new fractional methods, which have a clear physical basis and explicitly account for the effects due to protein structure and flexibility, provide an accurate picture of the contribution to solvation of different regions of the protein.

Computer Simulation↗
PubMed2004

Dynamic water networks in cytochrome C oxidase from Paracoccus denitrificans investigated by molecular dynamics simulations.

We present a molecular dynamics study of cytochrome c oxidase from Paracoccus denitrificans in the fully oxidized state, embedded in a fully hydrated dimyristoylphosphatidylcholine lipid bilayer membrane. Parallel simulations with different levels of protein hydration, 1.125 ns each in length, were carried out under conditions of constant temperature and pressure using three-dimensional periodic boundary conditions and full electrostatics to investigate the distribution and dynamics of water molecules and their corresponding hydrogen-bonded networks inside cytochrome c oxidase. The majority of the water molecules had residence times shorter than 100 ps, but a few water molecules are fixed inside the protein for up to 1.125 ns. The hydrogen-bonded network in cytochrome c oxidase is not uniformly distributed, and the degree of water arrangement is variable. The average number of solvent sites in the proton-conducting K- and D-pathways was determined. In contrast to single water files in narrow geometries we observe significant diffusion of individual water molecules along these pathways. The highly fluctuating hydrogen-bonded networks, combined with the significant diffusion of individual water molecules, provide a basis for the transfer of protons in cytochrome c oxidase, therefore leading to a better understanding of the mechanism of proton pumping.

Computer Simulation↗
PubMed1998

Structural and dynamic properties of the homodimeric hemoglobin from Scapharca inaequivalvis Thr-72-->Ile mutant: molecular dynamics simulation, low temperature visible absorption spectroscopy, and resonance Raman spectroscopy studies.

Molecular dynamics simulations, low temperature visible absorption spectroscopy, and resonance Raman spectroscopy have been performed on a mutant of the Scapharca inaequivalvis homodimeric hemoglobin, where residue threonine 72, at the subunit interface, has been substituted by isoleucine. Molecular dynamics simulation indicates that in the Thr-72-->Ile mutant several residues that have been shown to play a role in ligand binding fluctuate around orientations and distances similar to those observed in the x-ray structure of the CO derivative of the native hemoglobin, although the overall structure remains in the T state. Visible absorption spectroscopy data indicate that in the deoxy form the Soret band is less asymmetric in the mutant than in the native protein, suggesting a more planar heme structure; moreover, these data suggest a similar heme-solvent interaction in both the liganded and unliganded states of the mutant protein, at variance with that observed in the native protein. The "conformation sensitive" band III of the deoxy mutant protein is shifted to lower energy by >100 cm-1 with respect to the native one, about one-half of that observed in the low temperature photoproducts of both proteins, indicating a less polar or more hydrophobic heme environment. Resonance Raman spectroscopy data show a slight shift of the iron-proximal histidine stretching mode of the deoxy mutant toward lower frequency with respect to the native protein, which can be interpreted in terms of either a change in packing of the phenyl ring of Phe-97, as also observed from the simulation, or a loss of water in the heme pocket. In line with this latter interpretation, the number of water molecules that dynamically enters the intersubunit interface, as calculated by the molecular dynamics simulation, is lower in the mutant than in the native protein. The 10-ns photoproduct for the carbonmonoxy mutant derivative has a higher iron-proximal histidine stretching frequency than does the native protein. This suggests a subnanosecond relaxation that is slowed in the mutant, consistent with a stabilization of the R structure. Taken together, the molecular dynamics and the spectroscopic data indicate that the higher oxygen affinity displayed by the Thr-72-->Ile mutant is mainly due to a local perturbation in the dimer interface that propagates to the heme region, perturbing the polarity of the heme environment and propionate interactions. These changes are consistent with a destabilization of the T state and a stabilization of the R state in the mutant relative to the native protein.

Animals↗
PubMed2001

Comparison of different schemes to treat long-range electrostatic interactions in molecular dynamics simulations of a protein crystal.

Eight molecular dynamics simulations of a ubiquitin crystal unit cell were performed to investigate the effect of different schemes to treat the long-range electrostatic interactions as well as the need to include counter ions. A crystal system was chosen as the test system, because the higher charge density compared with a protein in solution makes it more sensitive to the way of treating the electrostatic interactions. Three different schemes of treating the long-range interactions were compared: straight cutoff, reaction-field approximation, and a lattice-sum method (P3M). For each of these schemes, two simulations were performed, one with and one without the counter ions. Two additional simulations with a reaction-field force and different initial placements of the counter ions were performed to examine the effect of the initial positions of the ions. The inclusion of long-range electrostatic interactions using either a reaction-field or a lattice-sum method proved to be necessary for the simulation of crystals. These two schemes did not differ much in their ability to reproduce the crystallographic structure. The inclusion of counter ions, on the other hand, seems not necessary for obtaining a stable simulation. The initial positions of the ions have a visible but small effect on the simulation.

Amino Acid Sequence↗
PubMed2003

Self-consistent-field modeling of hydrated unsaturated lipid bilayers in the liquid-crystal phase and comparison to molecular dynamics simulations.

A molecular-level self-consistent-field (SCF) theory is applied to model the lipid bilayer structures composed of 1-stearoyl-2-oleoyl-sn-glycero-3-phosphatidylcholine (18:0/18:1 omega 9cis PC) and 1-stearoyl-2-docosahexaenoyl-sn-glycero-3-phosphatidylcholine (18:0/22:6 omega 3cis PC). As compared to earlier attempts to model (saturated) PC membranes several additional features are implemented: (i) A water model is used which correctly leads to low water concentration in the bilayers. (ii) Free volume is allowed for, which is important to obtain bilayers in the fluid state. (iii) A polarization term is included in the segment potentials; this new feature corrects for a minor thermodynamic inconsistency present in (all) earlier results for charged bilayers. (iv) The CH3 groups in the lipid molecules are assumed to have twice the volume of a CH2 group; this leads to stable noninterdigitated bilayers. (v) A cis double bond is simulated by forcing gauche conformations along the sn-2 acyl chain. Results of an all-atom molecular dynamics (MD) simulation, using the collision dynamics method, on the same system are presented. Both SCF and MD prove, in accordance with experimental facts, that acyl unsaturation effectively reduces the length of the chain which counteracts interdigitation. It is also found that the phosphatidylcholine head group is lying almost flat on the membrane surface and the water penetrates into the bilayer upto the glycerol backbone units. From the SCF results it further followed that the free volume is not exactly evenly distributed over the bilayer. There is a small increase in free volume in the center of the bilayer as well as in the glycerol backbone region.

Journal Article↗
PubMed1991

Molecular dynamics simulation on a network of workstations using a machine-independent parallel programming language.

Molecular dynamics simulations investigate local and global motion in molecules. Several parallel computing approaches have been taken to attack the most computationally expensive phase of molecular simulations, the evaluation of long range interactions. This paper develops a straightforward but effective algorithm for molecular dynamics simulations using the machine-independent parallel programming language, Linda. The algorithm was run both on a shared memory parallel computer and on a network of high performance Unix workstations. Performance benchmarks were performed on both systems using two proteins. This algorithm offers a portable cost-effective alternative for molecular dynamics simulations. In view of the increasing numbers of networked workstations, this approach could help make molecular dynamics simulations more easily accessible to the research community.

Algorithms↗