SEARCH · PubMed Health

Results for “Molecular Dynamic simulation”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8Linked to original sources

PubMed1993

Molecular dynamics simulations of a lipid bilayer and of hexadecane: an investigation of membrane fluidity.

Molecular dynamics simulations of a fluid-phase dipalmitoyl phosphatidylcholine lipid bilayer in water and of neat hexadecane are reported and compared with nuclear magnetic resonance spin-lattice relaxation and quasi-elastic neutron scattering data. On the 100-picosecond time scale of the present simulations, there is effectively no difference in the reorientational dynamics of the carbons in the membrane interior and in pure hexadecane. Given that the calculated fast reorientational correlation times and the "microscopic" lateral diffusion of the lipids show excellent agreement with the experimental results, it is concluded that the apparently high viscosity of the membrane is more closely related to molecular interactions on the surface rather than in the interior.

1,2-Dipalmitoylphosphatidylcholine↗

PubMed1990

Molecular dynamics simulations of phospholipases A2.

An extensive molecular dynamics study of phospholipases A2 from pancreatic bovine and Crotalus atrox venom has shown that the well-conserved homologous core of the phospholipases A2, including the so called catalytic network, is very stable during the course of the calculations. The fluctuations which occur are located in segments which have significantly different three-dimensional conformations in the two phospholipases A2 studied, suggesting that a particularly stable core conformation gives rise to a large homologous family of similar three-dimensional structure. The calcium ion, which exhibits a crucial structural role in the monomeric phospholipases A2, appears not to be required to stabilize the C.atrox dimer. Moreover, the behaviour of the dimeric structure during the dynamics raises the question of a possible dissociation of the two subunits into functional monomers.

Animals↗

PubMed2023

A new horizon in the phosphorylated sites of AGA: the structural impact of C163S mutation in aspartylglucosaminuria through molecular dynamics simulation.

Aspartylglucosaminuria (AGU) is a lysosomal storage disorder caused by insufficient aspartylglucosaminidase (AGA) activity leading to chronic neurodegeneration. We utilized the PhosphoSitePlus tool to identify the AGA protein's phosphorylation sites. The phosphorylation was induced on the specific residue of the three-dimensional AGA protein, and the structural changes upon phosphorylation were studied via molecular dynamics simulation. Furthermore, the structural behaviour of C163S mutation and C163S mutation with adjacent phosphorylation was investigated. We have examined the structural impact of phosphorylated forms and C163S mutation in AGA. Molecular dynamics simulations (200 ns) exposed patterns of deviation, fluctuation, and change in compactness of Y178 phosphorylated AGA protein (Y178-p), T215 phosphorylated AGA protein (T215-p), T324 phosphorylated AGA protein (T324-p), C163S mutant AGA protein (C163S), and C163S mutation with Y178 phosphorylated AGA protein (C163S-Y178-p). Y178-p, T215-p, and C163S mutation demonstrated an increase in intramolecular hydrogen bonds, leading to greater compactness of the AGA forms. Principle component analysis (PCA) and Gibbs free energy of the phosphorylated/C163S mutation structures exhibit transition in motion/orientation than Wild type (WT). T215-p may be more dominant among these than the other studied phosphorylated forms. It might contribute to hydrolyzing L-asparagine functioning as an asparaginase, thereby regulating neurotransmitter activity. This study revealed structural insights into the phosphorylation of Y178, T215, and T324 in AGA protein. Additionally, it exposed the structural changes of the C163S mutation and C163S-Y178-p of AGA protein. This research will shed light on a better understanding of AGA's phosphorylated mechanism.Communicated by Ramaswamy H. Sarma.

Molecular Dynamics Simulation↗

PubMed1997

Molecular dynamics simulations of the unfolding of barnase in water and 8 M aqueous urea.

Molecular dynamics simulations of barnase have been conducted both in water and in 8 M urea solution for 500 ps at 25 degrees C and for 2000 ps at 85 degrees C. The final structure of the aqueous simulation at room temperature matches closely the structure obtained by NMR and the experimentally observed protections from isotopic exchange. The comparison of the structures generated by the aqueous simulation at 85 degrees C reveals a trajectory composed of groups of geometrically related structures separated by narrow regions of rapid change in structure. The first of these regions displays changes in backbone rmsd to the crystal structure and solvent-accessible area suggestive of a transition state, while the properties observed during the final 300 ps of the simulation are consistent with a stable intermediate. These assignments were confirmed by calculation of the "progress along the reaction coordinate" phi-values using an empirical equation based on a linear response method. The pathway of unfolding defined in this fashion agrees well with the experimental results of site-directed mutagenesis in terms of secondary structure content of the transition state and the intermediate and reproduces the relative stability of the different elements of secondary structure. The results of the simulations in urea suggest a mechanism at the molecular level for its well-known enhancement of the denaturation of proteins. The analysis of radial distribution functions shows that the first solvation shell of the protein is enriched in urea relative to the bulk solvent. The displacement of water molecules allows greater exposure of hydrophobic side chains, as witnessed particularly in the analysis of solvent-accessible surface areas at the higher temperature. Almost all urea molecules in the first shell form at least one hydrogen bond with the protein. They provide a more favorable environment for accommodation of the remaining water molecules, and they facilitate the separation of secondary structure elements by acting as a bridge between groups previously forming intraprotein hydrogen bonds.

Bacterial Proteins↗

PubMed1997

Molecular dynamics simulations of the N-linked oligosaccharide of the lectin from Erythrina corallodendron.

Molecular dynamics simulations have been used to model the flexibility of the seven-sugar oligosaccharide of the lectin from Erythrina corallodendron in three separate simulations: one of the isolated oligosaccharide in vacuo, one of the oligosaccharide in solution and one of the oligosaccharide linked to the protein in solution. Adiabatic conformational energy maps were prepared for each of the disaccharide linkages as a means of interpreting the observed dynamics and conformational averages in terms of intramolecular energy. The inclusion of aqueous solvent molecules appears to be necessary to reproduce the experimental conformational behavior, which also cannot be predicted well from conformational energy maps for the disaccharide linkages alone. The crystallographically determined conformation does not appear to be induced by the crystal dimerization, but is rather stable in solution. The build-up of fluctuations along the successive linkages of the oligosaccharide is significant and would be sufficient to prevent branch residues from being located in most crystal structure determinations. Good general agreement between the calculated solution structure and the average structure determined by NMR was found for most of the oligosaccharide linkages.

Carbohydrate Conformation↗

PubMed1992

FOCUS: a program for analyzing molecular dynamics simulations, featuring digital signal-processing techniques.

FOCUS is a program for analyzing molecular dynamics simulations. It enables the researcher to monitor structural and energetic properties during the trajectory, and to calculate the corresponding statistical averages, correlation functions and Fourier transforms. In addition to these conventional methods, the program also utilizes novel methods based on digital signal-processing techniques to characterize the various motions. The characteristic frequencies in the system are revealed by the frequency distribution function g(v), which is calculated from the Fourier transform of the atomic coordinates. A filtering technique is employed to remove uninteresting motion (e.g., high-frequency bond stretching) while retaining and focusing on important motion (e.g., low-frequency conformational motion). The filtering technique enables fast display of slow events without getting a blurry or jittery picture due to the high-frequency motions. Another new way for analyzing the motion is by extracting "characteristic modes" and associated frequencies. This yields a pictorial description of the oscillatory motions in a manner analogous to normal mode analysis.

Acetamides↗

PubMed1994

Molecular dynamics simulations of trp apo- and holorepressors: domain structure and ligand-protein interaction.

Molecular dynamics simulations of the apo- and holo-forms of the trp-repressor protein were performed under extensively solvated conditions in order to elucidate their dynamic structures and ligand-protein interactions. The root mean square fluctuations calculated from the trajectories agreed with those calculated from X-ray temperature factors. Distance, distance fluctuation, and dynamic cross-correlation maps were drawn to provide information on the dynamic structures and communications among the domains. A three-domain format has been proposed for the crystal structure (Zhang et al., Nature 327:591-597, 1987); namely, helices A-C and F of both subunits make up a central core, and D and E of each subunit forms a DNA binding head. The results of the simulations were mostly consistent with the three-domain format. However, helix F was more flexible and freer than other parts of the central core. The turn DE, the helix-turn-helix DNA binding motif, was free from interactions and correlations with other domains in both forms of the repressor. A comparison of the simulations of the aporepressor and holorepressor showed that tryptophan binding made the DNA-binding helix D more flexible but helix F less flexible. Several amino acid residues in contact with the bound tryptophan were identified as making concerted motions with it. Interaction energies between the corepressor and the amino acid residues of the protein were analyzed; the results were mostly consistent with the mutational experiments.

Amino Acid Sequence↗

PubMed1997

Molecular dynamics simulation study of DNA dodecamer d(CGCGAATTCGCG) in solution: conformation and hydration.

A molecular dynamics simulation of the dodecamer duplex d(CGCGAATTCGCG) using the particle mesh Ewald sum assumed a B-conformation remarkably close to the observed X-ray structure. The Ewald summation method effectively eliminates the usual "cut-off" of long-range interactions and allowed us to evaluate the full effect of the electrostatic forces. This simulation showed remarkable agreement with the Dickerson X-ray structure in both average structure and B-factors; within the EcoRI site itself, the rms deviation between the average theoretical and observed structures was 1.1 A. The width of the minor groove fluctuated between a wide and narrow configuration with the latter corresponding closely to the X-ray structure. The simulation also suggested a strong sequence-dependent signature on the minor groove width in both wide and narrow conformers. Hydration shells in both the major and minor grooves were observed. The "spine of hydration" in the minor groove was clear. In the major groove the first hydration shell appears to be a ribbon-like structure that reproduces the principal features of observed X-ray structures; subtle variations of this hydration pattern suggest sequence dependencies. Sequence-dependent features were also examined for helical and other geometric parameters. The successful reproduction of many experimentally observed fine structural features shows that the Ewald summation significantly improves the fidelity of the calculations.

Crystallography, X-Ray↗

PubMed1993

Molecular dynamics simulations and rigid body (TLS) analysis of aspartate carbamoyltransferase: evidence for an uncoupled R state.

In the R form of ATCase complexed with the bisubstrate analogue, N-(phosphonacetyl)-L-aspartate, large temperature factors are reported for the allosteric domains of the regulatory chains. We studied the conformational flexibility of the holoenzyme with molecular dynamics simulations and rigid body (TLS) analysis. The results of the molecular dynamics simulations suggest that, although local atomic fluctuations account for the temperature factors of the catalytic and zinc domains, they do not account for the large temperature factors of the allosteric regions. However, the temperature factors of the allosteric domains can be satisfactorily analyzed using a rigid body model. The simulations and rigid body analysis support the idea that the allosteric regions are mechanically uncoupled from the rest of the enzyme in the PALA structure. Implications of this uncoupling for allosteric regulation are discussed.

Allosteric Site↗

PubMed1997

Molecular dynamics simulation in solvent of the estrogen receptor protein DNA binding domain in complex with a non-consensus estrogen response element DNA sequence.

We investigated protein/DNA interactions, using molecular dynamics simulations computed between a 10 Angstom water layer model of the estrogen receptor (ER) protein DNA binding domain (DBD) amino acids and DNA of a non-consensus estrogen response element (ERE) consisting of 29 nucleotide base pairs. This ERE nucleotide sequence occurs naturally upstream of the Xenopus laevis Vitelligenin A1 gene. The ER DBD is encoded by three exons. Namely, exons 2 and 3 which encode the two zinc binding motifs and a sequence of exon 4 which encodes a predicted alpha helix. We generated a computer model of the ER DBD using atomic coordinates derived from the average of 30 nuclear magnetic resonance (NMR) spectroscopy coordinate sets. Amino acids on the carboxyl end of the ER DBD were disordered in both X-ray crystallography and NMR determinations and no coordinates were reported. This disordered region includes 10 amino acids of a predicted alpha helix encoded in exon 4 at the exon 3/4 splice junction. These amino acids are known to be important in DNA binding and are also believed to function as a nuclear translocation signal sequence for the ER protein. We generated a computer model of the predicted alpha helix consisting of the 10 amino acids encoded in exon 4 and attached this helix to the carboxyl end of the ER DBD at the exon 3/4 splice junction site. We docked the ER DBD model within the DNA major groove halfsites of the 29 base pair non-consensus ERE and flanking nucleotides. We constructed a solvated model with the ER DBD/ERE complex surrounded by a ten Angstrom water layer and conducted molecular dynamics simulations. Hydrogen bonding interactions were monitored. In addition, van der Waals and electrostatic interaction energies were calculated. Amino acids of the ER DBD DNA recognition helix formed both direct and water mediated hydrogen bonds at cognate codon-anticodon nucleotide base and backbone sites within the ERE DNA right major groove halfsite. Amino acids of the ER DBD exon 4 encoded predicted alpha helix formed direct and water mediated H-bonds with base and backbone sites of their cognate codon-anticodon nucleotides within the minor grooves flanking the ERE DNA major groove halfsites. These interactions together induced bending of the DNA into the protein.

Animals↗

PubMed1993

Molecular dynamics simulation of a phospholipase A2-substrate complex.

We have used knowledge of the three-dimensional structure of phospholipids and phospholipases A2 together with biochemical data, computer graphics modelling and a 48 ps molecular dynamics simulation to predict the structure of a phospholipase A2-substrate complex. There is remarkable similarity between this predicted structure of enzyme-substrate complex and the structure that can be deduced from the observed enzyme-inhibitor complex. Molecular dynamics simulation highlights the importance of the calcium-ion in substrate binding and the persistence of the His-48 to water-hydrogen bond is compatible with the proposed role of this water molecule as the nucleophile in catalysis.

Animals↗

PubMed2026

Molecular docking, molecular dynamics simulation, and enzyme inhibitory studies of vitamin K family members on aldose reductase.

Aldose reductase (AR) is a key enzyme in the polyol pathway and plays a major role in the progression of secondary complications of diabetes. Despite extensive efforts to develop natural and synthetic aldose reductase inhibitors (ARIs), most candidates have shown limited clinical efficacy, highlighting the need for more potent and selective inhibitors. In this study, we have systematically evaluated the inhibitory potential of vitamin K family members (vitamin K1, vitamin K2, and vitamin K3) using molecular docking, protein-ligand interaction analysis, molecular dynamics simulations, and enzyme kinetics. Docking analysis predicted that vitamin K2 has the highest binding affinity for AR. Subsequent molecular dynamics simulations revealed that both vitamin K1 and vitamin K2 formed stable complexes with the protein, exhibiting comparable RMSD (∼0.5 Å difference), similar RMSF profiles, and reduced radius of gyration, indicating compact and stable binding. Interaction analysis demonstrated that ligand binding is predominantly driven by hydrophobic interactions, with vitamin K2 forming a higher number of hydrophobic contacts, while vitamin K1 exhibited slightly more hydrogen bonding. Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) results further supports stronger binding of vitamin K2 (-56 kcal/mol) compared to vitaminK1 (-51 kcal/mol). Consistent with these findings, enzyme kinetics showed a slightly lower Ki value for vitamin K2 than vitamin K1. In contrast, vitamin K3 failed to maintain stable binding and moved out of the active site during simulation. Overall, the study highlights that hydrophobic interaction-driven stabilization plays a key role in ligand binding, and identifies vitamin K1 and vitamin K2 as promising inhibitors against AR, with vitamin K2 exhibiting more favourable hydrophobic interactions and binding stability.

Aldose Reductase↗

PubMed1994

Molecular dynamics simulations of an enzyme surrounded by vacuum, water, or a hydrophobic solvent.

We report on molecular dynamics simulations of a medium-sized protein, a lipase from Rhizomucor miehei, in vacuum, in water, and in a nonpolar solvent, methyl hexanoate. Depending on force field and solvent, the molecular dynamics structures obtained as averages over 150 ps had root-mean-square deviations in the range of 1.9 to 3.6 A from the crystal structure. The largest differences between the structures were in hydrogen bonding and exposed surface areas of the protein. The surface area increased in both solvents and became smaller in vacuum. The change of surface exposure varied greatly between different residues and occurred in accordance with the hydrophobicity of the residue and the nature of the solvent. The fluctuations of the atoms were largest in the external loops and agreed well with crystallographic temperature factors. Root-mean-square fluctuations were significantly smaller in the nonpolar solvents than they were in water, which is in accordance with the notion that proteins become more rigid in nonpolar solvents. In methyl hexanoate a partial opening of the lid covering the active site occurred, letting a methyl hexanoate molecule approach the active site.

Computer Simulation↗

PubMed1997

Transmembrane helix structure, dynamics, and interactions: multi-nanosecond molecular dynamics simulations.

To probe the fundamentals of membrane/protein interactions, all-atom multi-nanosecond molecular dynamics simulations were conducted on a single transmembrane poly(32)alanine helix in a fully solvated dimyristoyphosphatidylcholine (DMPC) bilayer. The central 12 residues, which interact only with the lipid hydrocarbon chains, maintained a very stable helical structure. Helical regions extended beyond these central 12 residues, but interactions with the lipid fatty-acyl ester linkages, the lipid headgroups, and water molecules made the helix less stable in this region. The C and N termini, exposed largely to water, existed as random coils. As a whole, the helix tilted substantially, from perpendicular to the bilayer plane (0 degree) to a 30 degrees tilt. The helix experienced a bend at its middle, and the two halves of the helix at times assumed substantially different tilts. Frequent hydrogen bonding, of up to 0.7 ns in duration, occurred between peptide and lipid molecules. This resulted in correlated translational diffusion between the helix and a few lipid molecules. Because of the large variation in lipid conformation, the lipid environment of the peptide was not well defined in terms of "annular" lipids and on average consisted of 18 lipid molecules. When compared with a "neat" bilayer without peptide, no significant difference was seen in the bilayer thickness, lipid conformations or diffusion, or headgroup orientation. However, the lipid hydrocarbon chain order parameters showed a significant decrease in order, especially in those methylene groups closest to the headgroup.

Computer Simulation↗

PubMed1998

Deducing polymeric structure from aqueous molecular dynamics simulations of oligosaccharides: predictions from simulations of hyaluronan tetrasaccharides compared with hydrodynamic and X-ray fibre diffraction data.

Molecular dynamics simulations of the two hyaluronan tetrasaccharides in water predict that over a period of 500 ps, their central linkages populate a single primary minima. Over the same period the peripheral linkages explore this minima, but also a secondary minima. Structures constructed using the primary minima were found to be extended left-handed helices of axial rise per disaccharide (h) 0.8 to 1.0 nm and 2.8 to 4.5 disaccharides per turn (n), in good agreement with n=3 and n=4 helices found by X-ray fibre diffraction studies. We have used the predicted average conformation from molecular dynamics to calculate the translational diffusion coefficients of the oligosaccharide series up to decasaccharide, and compared these with experimental measurements obtained using the method of capillary dispersion. Our calculated values are found to be in good agreement with experiment beyond the size of a tetrasaccharide. A partial digest of hyaluronan in the molecular mass range 10 to 100 kDa was fractionated by gel chromatography. Molecular weights were determined by in-line laser light-scattering measurements, and the translational diffusion coefficients of selected fractions were determined by dynamic laser light-scattering. A similar experiment was performed on hyaluronan with a molecular mass greater than 1MDa. The data suggest a change from rod-like to stiff coil behaviour beyond a molecular weight of 10 kDa. We have also examined the conformations available using the secondary minima, found at the peripheral linkages. In contrast to the extended structures previously described we have found left and right-handed helices with high values of n (5-10) and low values of h. Although there is no experimental evidence for these structures, they are of interest as, over short stretches, they would introduce folds, loops, and turns into the hyaluronan molecule. Such shapes may play an important role in the hydrodynamics of hyaluronan and its interaction with lipids and proteins.

Carbohydrate Conformation↗

PubMed1996

The application of different solvation and electrostatic models in molecular dynamics simulations of ubiquitin: how well is the X-ray structure "maintained"?

We present molecular dynamics simulations on ubiquitin with explicit solvent molecules and investigate the influence of different force fields [Weiner et al. (J. Am. Chem. Soc. 106:765-784, 1984; J. Comput. Chem. 7:230-252, 1986) vs. Cornell et al. (J. Am. Chem. Soc. 117:5179-5197, 1995)], different treatments of the long-range electrostatic interaction (8 A cutoff vs. particle mesh Ewald), and different solvation models (periodic box vs. small shell of water molecules) on the structure and the dynamics of the protein. Structural data are monitored by atomic root mean square deviations (RMSDs) from the crystal structure, the radius of gyration, the solvent-accessible surface area, and the pattern of the backbone hydrogen bonds. The dynamic behavior is assessed by the atomic fluctuations and the order parameters of the N-H backbone vectors. With the Cornell et al. force field and a periodic box model, the simulated structures stay much closer to the experimental X-ray structure than with the older Weiner et al. force field. A further improvement of the simulation is found when the electrostatic interaction is evaluated with the particle mesh Ewald method; after 1.2 ns of simulation the backbone RMSD amounts to only 1.13 A. The analysis of the dynamic parameters shows that this good structural agreement is not due to a damping of internal motion in the protein. For a given length of simulation time, the shell models achieve an agreement between simulated and experimental structures that is comparable to the best models that employ a periodic box of solvent models. However, compared with the box models, the fluctuations of the protein atoms in the shell models are smaller, and only with simulation times as long as 2 ns do they become of comparable size to the experimental ones.

Computer Simulation↗

PubMed1987

Anisotropy and anharmonicity of atomic fluctuations in proteins: analysis of a molecular dynamics simulation.

Positional probability density functions (pdf) for the atomic fluctuations are determined from a molecular dynamics simulation for hen egg-white lysozyme. Most atoms are found to have motions that are highly anisotropic but only slightly anharmonic. The largest deviations from harmonic motion are in the direction of the largest rms fluctuations in the local principal axis frame. Backbone atoms tend to be more nearly harmonic than sidechain atoms. The atoms with the largest anharmonicities tend to have pdfs with multiple peaks, each of which is close to harmonic. Several model pdfs are evaluated on the basis of how well they fit probability densities from the dynamics simulations when parameterized in terms of the moments of the distribution. Gram-Charlier and Edgeworth perturbation expansions, which have been successful in describing the motions of small molecules in crystals, are shown to be inadequate for the distributions found in the dynamics of proteins. Multipeaked distribution functions are found to be more appropriate.

Animals↗

PubMed1998

Molecular Dynamic Simulations of Single, Interacting, and Sheared Double Layers.

This article presents the results of molecular dynamic simulations of double layers developing on a flat mineral surface in aqueous sodium chloride solution. This part depicts the simulation results compared to classical double layer theory. Ion distribution, solvent orientation, and potential curve were derived for this propose. Dynamical values were calculated in the steady state and used to assess the idea of two parts of a double layer, stagnant and mobile, separated by a shear plane. Copyright 1998 Academic Press.

Journal Article↗