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 811 records · Page 45Linked to original sources

PubMed2001

Understanding beta-hairpin formation by molecular dynamics simulations of unfolding.

We have studied the mechanism of formation of a 16-residue beta-hairpin from the protein GB1 using molecular dynamics simulations in an aqueous environment. The analysis of unfolding trajectories at high temperatures suggests a refolding pathway consisting of several transient intermediates. The changes in the interaction energies of residues are related with the structural changes during the unfolding of the hairpin. The electrostatic energies of the residues in the turn region are found to be responsible for the transition between the folded state and the hydrophobic core state. The van der Waals interaction energies of the residues in the hydrophobic core reflect the behavior of the radius of gyration of the core region. We have examined the opposing influences of the protein-protein (PP) energy, which favors the native state, and the protein-solvent (PS) energy, which favors unfolding, in the formation of the beta-hairpin structure. It is found that the behavior of the electrostatic components of PP and PS energies reflects the structural changes associated with the loss of backbone hydrogen bonding. Relative changes in the PP and PS van der Waals interactions are related with the disruption of the hydrophobic core of a protein. The results of the simulations support the hydrophobic collapse mechanism of beta-hairpin folding.

Animals↗

PubMedDate not supplied

Molecular dynamics simulations to determine the effect of supercritical carbon dioxide on the structural integrity of hen egg white lysozyme.

In this study, various molecular dynamics simulations were conducted to investigate the effect of supercritical carbon dioxide on the structural integrity of hen egg white lysozyme. The analyses of backbone root-mean-square deviation, radius of gyration, and secondary structure stability all show that supercritical CO(2) exhibits the ability to increase the stability of this protein, probably as a result of the solvent with less polarity, where hydrophobic interactions stabilizing the native structure are weakened and simultaneously the local hydrogen bonds are strengthened, resulting in stabilization of the secondary structures. The hydrophobic cores in the alpha- and beta-domains also play an important role in preventing this protein from thermal unfolding. As supercritical CO(2) has been attractive for biomedical applications because of the advantages of mild critical condition, nonflammability, nontoxity, and the purity of the resulting products, the structural stabilizing effect found in this study strongly suggests that it is possible to increase the thermostability of hen egg white lysozyme by pretreatment with supercritical CO(2), leading to better industrial applications of this protein.

Carbon Dioxide↗

PubMed2005

Some insights into the stereochemistry of inhibition of macrophage migration inhibitory factor with 2-fluoro-p-hydroxycinnamate and its analogues from molecular dynamics simulations.

Macrophage migration inhibitory factor (MIF) exhibits tautomerase activity on phenylpyruvate and has E-stereochemistry preference. To investigate the binding modes of its competitive inhibitors and evaluate their binding affinities, molecular dynamics simulations together with MM-PBSA (molecular mechanics Poisson-Boltzmann surface area) analysis were performed on MIF complexed with (E)-2-fluoro-p-hydroxycinnamate and five analogues. Pro-1 was discovered to form a bifurcated hydrogen bond between its protonated nitrogen and carboxylate oxygens of E-ligands and Tyr-36. No hydrogen bonds were found between Pro-1 and Z-ligands. This distinct binding characteristic of E- and Z-ligands with Pro-1 may be the main factor for the large difference in their binding affinities, which is consistent with the previous report that Pro-1 is essential for the catalytic activity of MIF. MM-PBSA analysis revealed that energy components including van der Waals, electrostatic, and hydrophobic interactions are in favor of binding, among which electrostatic interactions are predominant to the binding affinity difference.

Cinnamates↗

PubMed2005

Thermal stability of benzorods: molecular-dynamics simulations.

Thermal stability of benzorods 2C6-20C6, which are obtained by stacking n (n=2-20) dehydrogenated benzene, have been investigated by molecular-dynamics simulations. It has been found that these structures assume a geometrical form depending on the number of dehydrogenated benzene layers, and they are stable under heat treatment up to elevated temperatures with a dependence on length.

Benzene Derivatives↗

PubMed2006

Molecular dynamics simulations of nanoparticles in dense isotropic nematogens: the role of matrix-induced long-range repulsive interactions.

We have carried out molecular dynamics simulation studies of binary mixtures of spherical nanoparticles (NPs) in a matrix of dense isotropic rod-shaped nematogens, with the size of the nematogen length being similar to that of the NP diameter. NPs at even low concentrations were found to shift the isotropic-nematic (I-N) transition significantly to higher pressure at a given temperature, indicative of long-range perturbation of the nematogenic matrix by the NPs. The NPs were found to be dispersed in the dense isotropic nematogenic matrix over a wide range of NP concentrations due to long-range (compared with the molecular size of the nematogens) repulsion caused by NP-induced local order fluctuations and reduced local orientational correlation in the isotropic nematogenic matrix, in contrast to the phase separation predicted and observed in other studies where the particles were much larger or smaller than the nematogens. Furthermore, since the repulsion observed in the NP-nematogen mixtures is only microscopically long range (on the order of about ten molecular lengths of the nematogens), globally ordered clustering observed in mixtures of colloidal particles in nematic matrices resulting from macroscopically long-range interaction is not observed in our simulations.

Journal Article↗

PubMed1998

Interaction between dihydropyridines and phospholipid bilayers: a molecular dynamics simulation.

Interaction of the calcium-channel antagonist dihydropyridines (DHPs), lacidipine and nifedipine, with a phospholipid bilayer was studied using 600 ps molecular dynamic simulations. We have constructed a double layer membrane model composed of 42 dimirystoyl-phosphatidylcholine molecules. The DHP molecules locate at about 7 A from the centre of the membrane, inducing an asymmetry in the bilayer. While lacidipine did not induce significant local perturbations as judged by the gauche-trans isomerisation rate, nifedipine significantly decreased this rate, probably by producing a local rigidity of the membrane in the vicinity of the DHP.

Calcium Channel Blockers↗

PubMed1998

Molecular dynamics simulations on parallel and antiparallel C.G*G triplexes.

Molecular dynamics (MD) studies have been carried out on the Hoogsteen hydrogen bonded parallel and the reverse Hoogsteen hydrogen bonded antiparallel C.G*G triplexes. Earlier, the molecular mechanics studies had shown that the parallel structure was energetically more favourable than the antiparallel structure. To characterize the structural stability of the two triplexes and to investigate whether the antiparallel structure can transit to an energetically more favourable structure, due to the local fluctuations in the structure during the MD simulation, the two structures were subjected to 200ps of constant temperature vacuum MD simulations at 300K. Initially no constraints were applied to the structures and it was observed that for the antiparallel triplex, the structure showed a large root mean square deviation from the starting structure within the first 12ps and the N4-H41--O6 hydrogen bond in the WC duplex got distorted due to a high propeller twist and a moderate increase in the opening angle in the basepairs. Starting from an initial value of 30 degrees , helical twist of the average structure from this simulation had a value of 36 degrees , while the parallel structure stabilized at a twist of 33 degrees. In spite of the hydrogen bond distortions in the antiparallel triplex, it was energetically comparable to the parallel triplex. To examine the structural characteristics of an undistorted structure, another MD simulation was performed on the antiparallel triplex by constraining all the hydrogen bonds. This structure stabilized at an average twist of 33 degrees. In the course of the dynamics though the energy of the molecule - compared to the initial structure - improved, it did not become comparable to the parallel structure. Energy minimization studies performed in the presence of explicit water and counterions also showed the two structures to be equally favourable energetically. Together these results indicate that the parallel C.G*G triplex with Hoogsteen hydrogen bonds also represents a stereochemically and energetically favourable structure for this class of triplexes.

Carbohydrates↗

PubMed2000

Glycosaminoglycan conformation: do aqueous molecular dynamics simulations agree with x-ray fiber diffraction?

Glycosaminoglycan-protein interactions are biologically important and require an appreciation of glycan molecular shape in solution, which is presently unavailable. In previous studies we found strong similarity between aqueous molecular dynamics (MD) simulations and published x-ray diffraction refinements of hyaluronan. We have applied a similar approach here to chondroitin and dermatan, attempting to clarify some of the issues raised by the x-ray diffraction literature relating to chondroitin and dermatan sulfate. We predict that chondroitin has the same beta(1-->4) linkage conformation as hyaluronan, and that their average beta(1-->3) conformations differ. This is explained by changes in hydrogen-bonding across this linkage, resulting from its axial hydroxyl, causing a different sampling of left-handed helices in chondroitin (2.5- to 3.5-fold) as compared with hyaluronan (3.0- to 4.0-fold). Few right-handed helices, which lack intramolecular hydrogen-bonds, were sampled during our MD simulations. Thus, we propose that the 8-fold helix observed in chondroitin-6-sulfate, represented in the literature as an 8(3) helix (right-handed), though it has never been refined, is more likely to be 8(5) (left-handed) helix. Molecular dynamics simulations implied that (4)C(1) and (2)S(O), but not (1)C(4), forms of iduronate could be used in refinements of dermatan x-ray fiber diffraction patterns. Current models of 8-fold dermatan sulfate chains containing (4)C(1) iduronate refine to right-handed helices, which possess no intramolecular hydrogen-bonds. However, MD simulations predict that models containing (2)S(O) iduronate could provide better (8(5) helix) starting structures for refinement. Thus, the 8-fold dermatan sulfate refinement (8(3) helix) could be in error.

Carbohydrate Conformation↗

PubMed2025

Exploring potential targets and molecular mechanisms of traumatic brain injury exacerbated by Benzo(a)pyrene via network toxicology and molecular dynamics simulation.

Benzo(a)pyrene (BaP) is a common environmental pollutant from combustion sources that promotes oxidative stress, neuroinflammation and disruption of blood-brain barrier (BBB). However, its contribution to worsening traumatic brain injury (TBI) remains unclear. In this study, we aimed to assess the contribution of BaP to secondary injury in TBI. By integrating data from e.g., the Comparative Toxicogenomics Database, GeneCards, and Online Mendelian Inheritance in Man, 121 overlapping core targets were identified between BaP and TBI. Enrichment analyses via Gene Ontology and Kyoto Encyclopedia of Genes and Genomes, combined with protein-protein interaction networks and topological algorithms (degree, closeness centrality, betweenness centrality, average shortest path length, topological coefficient and partner of multi-edged node pairs), highlighted five hub genes (TP53, EGFR, AKT1, ACTB, and TNF) implicated in mitogen-activated protein kinase signaling, oxidative stress, and neuroinflammation. Molecular docking showed strong binding affinities of BaP to these hub proteins, with energies from -9.3 to -12.1 kcal/mol, tighter than co-crystal ligands and existing protein-binding drugs. Molecular dynamics simulations confirmed interaction stability through low root-mean-square deviation (< 0.5 nm), fluctuation, and radius of gyration values. Calculation of binding free energies using MM-PBSA validated the strong binding affinity between BaP and binding pockets of each hub genes. Toxicity prediction analysis revealed an oral LD50 of 316 mg/kg for BaP, with high probabilities for neurotoxicity, BBB permeability, carcinogenicity, and mutagenicity, associated with aryl hydrocarbon receptor activation. These findings reveal a "neurovascular homeostasis disruption" network underlying BaP-exacerbated TBI pathology and highlight potential targets to reduce pollution-related risks in TBI management.

Benzo(a)pyrene↗

PubMed2006

Molecular-dynamics simulations of binary structure II hydrogen and tetrahydrofurane clathrates.

The binary structure II hydrogen and tetrahydrofurane (THF) clathrates are studied with molecular-dynamics simulations. Simulations are done at pressures of 120 and 1.013 bars for temperatures ranging from 100 to 273 K. For the small cages of the structure II unit cell, H2 guest molecule occupancies of 0, 16 (single occupancy), and 32 (double occupancy) are considered. THF occupancies of 0-8 in the large cages are studied. For cases in which THF does not occupy all large cages in a unit cell, the remaining large cages can be occupied with sets of four H2 guest molecules. The unit-cell volumes and configurational energies are compared in the different occupancy cases. Increasing the small cage occupancy leads to an increase in the unit-cell volume and thermal-expansion coefficient. Among simulations with the same small cage occupancy, those with the large cages containing 4H2 guests have the largest volumes. The THF guest molecules have a stabilizing effect on the clathrate and the configurational energy of the unit cell decreases linearly as the THF content increases. For binary THF + H2 clathrates, the substitution of the THF molecules in the large cages with sets of 4H2 molecules increases the configurational energy. For the binary clathrates, various combinations of THF and H2 occupancies have similar configurational energies.

Journal Article↗

PubMed2004

Exploring the counterion atmosphere around DNA: what can be learned from molecular dynamics simulations?

The counterion distribution around a DNA dodecamer (5'-CGCGAATTCGCG-3') is analyzed using both standard and novel techniques based on state of the art molecular dynamics simulations. Specifically, we have explored the population of Na(+) in the minor groove of DNA duplex, and whether or not a string of Na(+) can replace the spine of hydration in the narrow AATT minor groove. The results suggest that the insertion of Na(+) in the minor groove is a very rare event, but that when once the ion finds specific sites deep inside the groove it can reside there for very long periods of time. According to our simulation the presence of Na(+) inside the groove does not have a dramatic influence in the structure or dynamics of the duplex DNA. The ability of current MD simulations to obtain equilibrated pictures of the counterion atmosphere around DNA is critically discussed.

Binding Sites↗

PubMed1999

Molecular dynamics simulations of human glutathione transferase P1-1: analysis of the induced-fit mechanism by GSH binding.

We report here a 1-ns molecular dynamics simulation on the ligand-free monomer of human glutathione transferase P1-1 in bulk water. The average conformation obtained from the last 500 ps of simulation is taken as a model for the apo-structure of this protein and compared to the available crystallographic data. Remarkable changes in the tertiary structure take place during the simulation and are ascribed to the removal of the ligand. They support an induced fit mechanism occurring upon glutathione binding, whose major features can be described in detail. A portion of helix 2 (residues 42-50), which participates in the formation of the active site, undergoes the most prominent conformational changes. Other protein segments, such as the C-terminal loop and helix 4, also show relevant structural rearrangements. All these transitions cause a significant shielding from the solvent of the hydrophobic binding site of the co-substrate, whose exposed surface goes from 4.6 nm(2) in the holo-structure to about 3.1 nm(2) in the apo-conformation. The results of this simulation are consistent with numerous experimental observations previously obtained on GST P1-1 and provide new insights for their explanation at the molecular level. Proteins 1999;37:1-9.

Binding Sites↗

PubMed2006

Ab initio molecular dynamics simulation of the structure and proton transport dynamics of methanol-water solutions.

Ab initio molecular dynamics simulations are employed to study the structural and proton transport properties of methanol-water mixtures. Structural characteristics analyzed at two different methanol mole fractions (X(M) = 0.25 and X(M) = 0.5) reveal enhanced structuring of water as the methanol mole fraction increases in agreement with recent neutron diffraction experiments. The simulations reveal the existence of separate hydrogen-bonded water and methanol networks, also in agreement with the neutron diffraction data. The addition of a single proton to the X(M) = 0.5 mixture leads to an anomalous structural or Grotthuss-type diffusion mechanism of the charge defect in which water-to-water, methanol-to-water, and water-to-methanol proton transfer reactions play the dominant role with methanol-to-methanol transfers being much less significant. Unlike in bulk water, where coordination number fluctuations drive the proton transport process, suppression of the coordination number of waters in the first solvation shell of the defect diminish the importance of coordination number fluctuations as a driving force in the structural diffusion process. The charge defect is found to reside preferentially at the interface between water and methanol networks. The length of the ab initio molecular dynamics run (approximately 120 ps), allowed the diffusion constant of the charge defect to be computed, yielding a value of D = 4.2 x 10(-5) cm2/s when deuterium masses are assigned to all protons in the system. The relation of this value to excess proton diffusion in bulk water is discussed. Finally, a kinetic theory is introduced to identify the relevant time scales in the proton transfer/transport process.

Journal Article↗

PubMed2005

Transmembrane peptides from tyrosine kinase receptor. Mutation-related behavior in a lipid bilayer investigated by molecular dynamics simulations.

Polar mutations in transmembrane alpha helices may alter the structural details of the hydrophobic sequences and control intermolecular contacts. We have performed molecular dynamics simulations on the transmembrane domain of the proto-oncogenic and the oncogenic forms of the Neu receptor in a fluid DMPC bilayer to test whether the Glu mutation which replaces the Val residue at position 664 may alter the helical structure and its insertion in the membrane. The simulations show that the wild and the mutant forms of the transmembrane domain have a different behavior in the bilayer. The native transmembrane sequence is found to be more flexible than in the presence of the Glu mutation, characterized by a tendency to pi deformation to accommodate the helix length to the membrane thickness. The mutant form of this domain does not evidence helical deformation in the present simulation. Hydrophobic matching is achieved both by a larger helix tilt and a vertical shift of the helix towards the membrane interface, favoring the accessibility of the Glu side chain to the membrane environment. A rapid exchange of hydrogen bond interactions with the surrounding water molecules and the lipid headgroups is observed. The difference in the behavior between the two peptides in a membrane environment was also observed experimentally. Both simulation and experimental results agree with the hypothesis that water may act as an intermediate for the formation of cross links between the facing Glu side chains stabilizing the dimer.

Amino Acid Sequence↗

PubMed1997

Protein hinge bending as seen in molecular dynamics simulations of native and M61 mutant T4 lysozymes.

A dynamical model of interdomain "hinge bending" of T4 lysozyme in aqueous solution has been developed on the basis of molecular dynamics (MD) simulation. The MD model study provides a description of the conformational reorganization expected to occur for the protein in aqueous solution as compared to the crystalline environment. Three different 500 ps molecular dynamics simulations were calculated, each using a distinctly different crystal conformation of T4 lysozyme as the starting points of the MD simulations. Crystal structures of wild-type lysozyme and "open" and "closed" forms of M61 variant structures were analyzed in this study. Large-scale, molecular-conformational rearrangements were observed in all three simulations, and the largest structural change was found for the open form of the M61 allomorph. All three simulated proteins had closed relative to the wild-type crystal structure, and the closure of the "jaws" of the active site cleft occurred gradually over the time course of the trajectories. The time average MD structures, calculated over the final 50 ps of each trajectory, had all adapted to conformations more similar to each other than to their incipient crystal forms. Using a similar MD protocol on cytochrome P450BM-3 [M. D. Paulsen and R. L. Ornstein (1995) Proteins: Structure Function and Genetics, Vol. 27, pp. 237-243] we have found that the opposite type of motion relative to the starting crystal structure, that is, the open form of the crystal structure, had opened to a greater degree relative to the incipient crystal structure form. Therefore we do not believe that either result is merely a simulation artifact, but rather the protein dynamics are due to protein relaxation in the absence of crystal packing forces in the simulated solution environments.

Bacteriophage T4↗

PubMed2005

On the influence of charged side chains on the folding-unfolding equilibrium of beta-peptides: a molecular dynamics simulation study.

The influence of charged side chains on the folding-unfolding equilibrium of beta-peptides was investigated by means of molecular dynamics simulations. Four different peptides containing only negatively charged side chains, positively charged side chains, both types of charged side chains (with the ability to form stabilizing salt bridges) or no charged side chains were studied under various conditions (different simulation temperatures, starting structures and solvent environment). The NMR solution structure in methanol of one of the peptides (A) has already been published; the synthesis and NMR analysis of another peptide (B) is described here. The other peptides (C and D) studied herein have hitherto not been synthesized. All four peptides A-D are expected to adopt a left-handed 3(14)-helix in solution as well as in the simulations. The resulting ensembles of structures were analyzed in terms of conformational space sampled by the peptides, folding behavior, structural properties such as hydrogen bonding, side chain-side chain and side chain-backbone interactions and in terms of the level of agreement with the NMR data available for two of the peptides. It was found that the presence of charged side chains significantly slows down the folding process in methanol solution due to the stabilization of intermediate conformers with side chain-backbone interactions. In water, where the solvent competes with the solute-solute polar interactions, the folding process to the 3(14)-helix is faster in the simulations.

Hydrogen Bonding↗

PubMed2006

Structural change in lipid bilayers and water penetration induced by shock waves: molecular dynamics simulations.

The structural change of a phospholipid bilayer in water under the action of a shock wave is numerically studied with unsteady nonequilibrium molecular dynamics simulations. The action of shock waves is modeled by the momentum change of water molecules, and thereby we demonstrate that the resulting collapse and rebound of the bilayer are followed by the penetration of water molecules into the hydrophobic region of the bilayer. The high-speed phenomenon that occurs during the collapse and rebound of the bilayer is analyzed in detail, particularly focusing on the change of bilayer thickness, the acyl chain bend angles, the lateral fluidity of lipid molecules, and the penetration rate of water molecules. The result shows that the high-speed phenomenon can be divided into two stages: in the first stage the thickness of bilayer and the order parameter are rapidly reduced, and then in the second stage they are recovered relatively slowly. It is in the second stage that water molecules are steadily introduced into the hydrophobic region. The penetration of water molecules is enhanced by the shock wave impulse and this qualitatively agrees with a recent experimental result.

1,2-Dipalmitoylphosphatidylcholine↗

PubMed2006

Membrane assembly of simple helix homo-oligomers studied via molecular dynamics simulations.

The assembly of simple transmembrane helix homo-oligomers is studied by combining a generalized Born implicit membrane model with replica exchange molecular dynamics simulations to sample the conformational space of various oligomerization states and the native oligomeric conformation. Our approach is applied to predict the structures of transmembrane helices of three proteins--glycophorin A, the M2 proton channel, and phospholamban--using only peptide sequence and the native oligomerization state information. In every case, the methodology reproduces native conformations that are in good agreement with available experimental structural data. Thus, our method should be useful in the prediction of native structures of transmembrane domains of other peptides. When we ignore the experimental constraint on the native oligomerization state and attempt de novo prediction of the structure and oligomerization state based only on sequence and simple energetic considerations, we identify the pentamer as the most stable oligomer for phospholamban. However, for the glycophorin A and the M2 proton channels, we tend to predict higher oligomers as more stable. Our studies demonstrate that reliable predictions of the structure of transmembrane helical oligomers can be achieved when the observed oligomerization state is imposed as a constraint, but that further efforts are needed for the de novo prediction of both structure and oligomeric state.

Cell Membrane↗