Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “MD 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 613 records · Page 34Linked to original sources
PubMed1996

The aqueous solution structure of a lipoteichoic acid from Streptococcus pneumoniae strain R6 containing 2,4-diamino-2,4,6-trideoxy-galactose: evidence for conformational mobility of the galactopyranose ring.

The 2D-NOESY spectra for the per-N-acetylated and the native lipoteichoic acid (LTA) oligomer from Streptococcus pneumoniae strain R6 clearly indicate a difference in conformation of the 2,4,6-trideoxy-galactopyranose ring. Whereas the 2,4-N-acetylated Gal24N adopts the usual 4C1 chair conformation, the native 2-N-acetyl-4-amino Gal24N exhibits conformational mobility with comparable populations in the 4C1 chair and 5S1 skew conformations, as determined using MD simulation for the partial trisaccharide Me-beta-D-Glc6P-(1-->3)-alpha-D-Gal24N-[6-PC]-(1-->4)-alpha- D-galNAc and from the intra-ring NOE effects. 31P-NMR spectra point to a strong electrostatic or hydrogen-bonding interaction between the free 4-NH2 group on the Gal24N and the negatively charged diester phosphate group between adjacent pentasaccharide repeating-units [Ribitol-(5-->6)-beta-D-Glc6P]. Molecular modelling and MD simulation experiments confirmed that such an interaction was feasible with the Gal24N galactopyranose ring in the inverted B1.4 or skew 5S1 conformation.

Acetylgalactosamine↗
PubMed2005

Allostery of actin filaments: molecular dynamics simulations and coarse-grained analysis.

The structural and mechanical properties of monomeric actin (G-actin), the trimer nucleus, and actin filaments (F-actins) are determined as a function of the conformation of the DNase I-binding loop (DB loop) by using all-atom molecular dynamics simulations and coarse-grained (CG) analysis. Recent x-ray structures of ADP-bound G-actin (G-ADP) by Otterbein et al. [Otterbein, L. R., Graceffa, P. & Dominguez, R. (2001) Science 293, 708-711] and ATP-bound G-actin (G-ATP) by Graceffa and Dominguez [Graceffa, P. & Dominguez, R. (2003) J. Biol. Chem. 278, 34172-34180] indicate that the DB loop of actin does not have a well defined secondary structure in the ATP state but folds into an alpha-helix in the ADP state. MD simulations and CG analysis indicate that such a helical DB loop significantly weakens the intermonomer interactions of actin assemblies and thus leads to a wider, shorter, and more disordered filament. The computed persistence lengths of F-actin composed of G-ATP (16 microm) and of G-ADP (8.5 microm) agree well with the experimental values for the two states. Therefore, the loop-to-helix transition of the DB loop may be one of the factors that lead to the changes in structural and mechanical properties of F-actin after ATP hydrolysis. This result may provide a direct connection between the conformational changes of an actin monomer and the structural and mechanical properties of the cytoskeleton. The information provided by MD simulations also helps to understand the possible origin of the special features of actin dynamics.

Actin Cytoskeleton↗
PubMed2006

Molecular dynamics of unstable motions and capillary instability in liquid nanojets.

We present an investigation of the capillary instability of nanometer-sized surface-tension-driven flow using molecular dynamics (MD) simulations with Lennard-Jones fluid. Unstable motions of a liquid nanojet are successfully simulated and it is found that the thermal fluctuation, which is significant in a nanoscale system, is the most important factor for various breakup scenarios of a nanojet. The nanojet diameter at the nozzle outlet is varied to show the effect of size on the rupture phenomena and the formation of small droplets. Numerical results for the rupture time and the growth rate of spherical droplets are compared with those of various classical linear instability theories. Even though the MD simulation results for the growth rate of the droplets are close to those predicted by the classical instability theories, the former are shown to be independent of the wave number unlike the latter. Therefore, the classical continuum-based theories may not be applicable to studying the instability of nanoscale systems.

Journal Article↗
PubMed1995

Spectral densities of nitrogen nuclei in Escherichia coli ribonuclease HI obtained by 15N NMR relaxation and molecular dynamics.

Spectral densities of the 15N amide in Escherichia coli ribonuclease HI, obtained from NMR relaxation experiments, were compared with those calculated using a molecular dynamics (MD) simulation. All calculations and comparisons assumed that the auto-correlation function describing the internal motions of the molecule was independent of the auto-correlation function associated with overall rotational diffusion. Comparisons were limited to those residues for which the auto-correlation function of internal motions rapidly relaxed and reached a steady state within 205 ps. The results show the importance of frequency components as well as amplitudes of internal motions in order to obtain a meaningful comparison of MD simulations with NMR data.

Computer Simulation↗
PubMed2006

G-quadruplexes can maintain their structure in the gas phase.

Several very extended (0.5-1 micros) molecular dynamics (MD) simulations of parallel and antiparallel G-quadruplex DNA strongly suggest that in the presence of suitable cations the quadruplex not only remains stable in the gas phase, but also displays a structure that closely resembles that found in extended (25-ns long) trajectories in aqueous solution. In the absence of the crucial cations, the trajectories become unstable and in general the quadruplex structure is lost. To our knowledge, this is the first physiologically relevant structure of DNA for which very large MD simulations suggest that the structure in water and in the gas phase are indistinguishable.

Computer Simulation↗
PubMed2005

Generalized Langevin theory on the dynamics of simple fluids under external fields.

A theory on the time development of the density and current fields of simple fluids under an external field is formulated through the generalized Langevin formalism. The theory is applied to the linear solvation dynamics of a fixed solute regarding the solute as the external field on the solvent. The solute-solvent-solvent three-body correlation function is taken into account through the hypernetted-chain integral equation theory, and the time correlation function of the random force is approximated by that in the absence of the solute. The theoretical results are compared with those of molecular-dynamics (MD) simulation and the surrogate theory. As for the transient response of the density field, our theory is shown to be free from the artifact of the surrogate theory that the solvent can penetrate into the repulsive core of the solute during the relaxation. We have also found a large quantitative improvement of the solvation correlation function compared with the surrogate theory. In particular, the short-time part of the solvation correlation function is in almost perfect agreement with that from the MD simulation, reflecting that the short-time expansion of the theoretical solvation correlation function is exact up to t(2) with the exact three-body correlation function. A quantitative improvement is found in the long-time region, too. Our theory is also applied to the force-force time correlation function of a fixed solute, and similar improvement is obtained, which suggests that our present theory can be a basis to improve the mode-coupling theory on the solute diffusion.

Journal Article↗
PubMed2005

Common semiopen conformations of Mg2+-free Ras, Rho, Rab, Arf, and Ran proteins combined with GDP and their similarity with GEF-bound forms.

A computational study was performed on the Mg(2+)-free conformations of the small guanine nucleotide-binding proteins (GNBPs): Ras, Rho, Rab, Arf, and Ran, which were complexed with GDP. Molecular dynamics (MD) simulation was executed for each complex for the duration of 3.0 ns to investigate the effects of Mg(2+) ions on the GNBPs' structure. The results indicated that all Mg(2+)-free GNBPs formed a groove between the switch region and the nucleotide-binding site. In some GNBP families, the release of Mg(2+) was reported to play an important role in binding the guanine nucleotide-exchanging factor (GEF) promoting the GDP/GTP exchange reaction. Interestingly, the grooves, which appeared in the MD simulations, were similar to the grooves experimentally observed in the GNBP-GEF complex. We also calculated the Mg(2+)-bound GNBPs to compare with the Mg(2+)-free forms. No groove was observed in the Mg(2+)-bound GNBPs. These results demonstrated a regulatory role of Mg(2+) ion to prepare a template for the GEF binding. Moreover, the results suggested that the release of Mg(2+) ion lead to the GEF-GNBP binding.

Amino Acid Sequence↗
PubMed2000

Structure/thermodynamics relationships of lectin-saccharide complexes: the Erythrina corallodendron case.

Molecular dynamics (MD) simulations of Erythrina corallodendron lectin binding to a monosaccharide, alpha-galactose, and a disaccharide, N-acetyl lactosamine, have been performed in order to investigate the relationship between structure and thermodynamics. A simulated annealing protocol has been used to generate ensembles of structures for the two complexes, from which both qualitative and quantitative information on binding dynamics have been extracted. The ensembled averaged lectin-saccharide interaction enthalpy is equivalent for both sugars, whereas the calculation based on the X-ray structures does show a difference. Within large statistical errors, the calculated 'binding enthalpy' is also the same for the two systems. These errors arise largely from terms involving solvent and are a typical limitation of current MD simulations. Significant qualitative differences in binding between the two complexes are, however, observed over the ensembles. These could be important for unraveling the structure/thermodynamic relationship. Stated simply, there are a greater number of binding options available to the disaccharide compared to the monosaccharide. The implications of alternative binding states on thermodynamic parameters and the 'breaking of enthalpy-entropy compensation' are discussed. The role of solvent in lectin-saccharide complex formation is suggested to be significant.

Amino Sugars↗
PubMed2003

Glass transition of hard sphere systems: molecular dynamics and density functional theory.

The glass transition of a hard sphere system is investigated within the framework of the density functional theory (DFT). Molecular dynamics (MD) simulations are performed to study the dynamical behavior of the system on the one hand and to provide the data to produce the density field for the DFT on the other hand. Energy landscape analysis based on the DFT shows that there appears a metastable (local) free energy minimum representing an amorphous state as the density is increased. This state turns out to become stable, compared with the uniform liquid, at some density around which we also observe a sharp slowing down of the alpha relaxation in the MD simulations.

Journal Article↗
PubMed1994

A crystallographic re-investigation into the structure of Streptomyces griseus proteinase A reveals an acyl-enzyme intermediate.

A molecular dynamics (MD) simulation on the contents of two asymmetric units of the crystal structure of the bacterial serine proteinase, Streptomyces griseus proteinase A (SGPA), resulted in four dihydrogen phosphate anions migrating to form a cluster in a solvent region far removed from the protein surface. In an effort to provide experimental verification for this unexpected phenomenon, native SGPA crystals were soaked in 2.0 M KH2AsO4; intensity data were collected and difference electron density maps examined for evidence of H2AsO4- ion clusters. These maps did not corroborate the cluster observed in the MD simulation. They did, however, show positive electron density features located in the active site cavity that could be interpreted as a tetrapeptide, Gly-Ala-Ser-(beta-OH)Asp, covalently bonded to O gamma of Ser195 as an acyl-enzyme intermediate. There was also a peak of electron density corresponding to an ideal position for the hydrolytic water in the deacylation reaction. This density is centred 3.1 A from His57 N epsilon 2 and 2.7 A from the carbonyl-carbon atom of the planar acyl-ester bond of the complex. The carbonyl-oxygen atom of the ester is located in the oxyanion pocket and participates in hydrogen bonds with Gly193 NH (2.6 A) and Ser195 NH (3.0 A).

Acylation↗
PubMed2026

Exploring the Potential Molecular Targets of Cyanidin-3-O-glucoside for Type 2 Diabetes Mellitus Treatment.

INTRODUCTION: This study aims to elucidate the multi-target molecular mechanism of cyanidin-3-O-glucoside (C3G) in treating Type 2 Diabetes (T2DM) through network pharmacology methods. METHODS: The study was designed to predict the targets of C3G through public databases and to screen for T2DM-related targets. Protein-protein interaction (PPI) network analysis, Gene Ontology (GO), and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis were performed on the common targets. Core targets were further validated through molecular docking and molecular dynamics (MD) simulations. RESULTS: This research identified a total of 57 potential targets of C3G in the treatment of T2DM. Subsequent PPI analysis identified ALB (Degree=43), AKT1 (Degree=41), and TNF (Degree=41) as the top three hub proteins. Pathway analysis indicated significant involvement in the insulin signaling pathway (P = 4.205×10-9), AMPK signaling pathway (P = 9.582×10-7), and FoxO signaling pathway (P = 1.315×10-6). Molecular docking revealed strong binding affinities between C3G and NOS3 (-9.5 kcal/mol), PPARG (-9.0 kcal/mol), TNF (-8.5 kcal/mol), and INSR (-8.4 kcal/mol). MD simulations further confirmed that the C3G-target complex has excellent binding stability. DISCUSSION: C3G may intervene in the pathological progression of T2DM by regulating key pathways such as insulin sensitivity, inflammatory responses, and oxidative stress. Further studies suggest that INSR and NOS3 may be new targets through which C3G exerts its effects, but their specific mechanisms and in vivo biological functions still need to be elucidated by subsequent experiments. CONCLUSION: C3G may intervene in the progression of T2DM in a multi-pathway synergistic manner by targeting key molecules such as INSR and NOS3.

Anthocyanins↗
PubMed2004

Promoting vibrations in human purine nucleoside phosphorylase. A molecular dynamics and hybrid quantum mechanical/molecular mechanical study.

Crystallographic studies of human purine nucleoside phosphorylase (hPNP) with several transition-state (TS) analogues in the immucillin family showed an unusual geometric arrangement of the atoms O-5', O-4', and O(P), the nucleophilic phosphate oxygen, lying in a close three-oxygen stack. These observations were corroborated by extensive experimental kinetic isotope effect analysis. We propose that protein-facilitated dynamic modes in hPNP cause this stack, centered on the ribosyl O-4' oxygen, to squeeze together and push electrons toward the purine ring, stabilizing the oxacarbenium character of the TS. As the N-ribosidic bond is cleaved during the reaction, the pK(a) values of N-7 and O-6 increase by the electron density expelled by the oxygen-stack compression toward the purine ring. Increased electron density in the purine ring improves electrostatic interactions with nearby residues and facilitates the abstraction of a proton from a solvent proton or an unidentified general acid, making the purine a better leaving group, and accelerating catalysis. Classical and mixed quantum/classical molecular dynamics (MD) simulations of the Michaelis complex of hPNP with the substrates guanosine and phosphate were performed to assess the existence of protein-promoting vibrations (PPVs). Analogous simulations were performed for the substrates in aqueous solution. In the catalytic site, the O-5', O-4', and O(P) oxygens vibrate at frequencies of ca. 125 and 465 cm(-1), as opposed to 285 cm(-1) in the absence of hPNP. The hybrid quantum mechanical/molecular mechanical method was used to assess whether this enzymatic vibration pushing the oxygens together is coupled to the reaction coordinate, and thus has a direct positive impact on catalysis. The potential energy surface for the phosphorolysis reaction for several snapshots taken from the classical MD simulation showed substantial differences in oxygen compression. Our calculations showed the existence of PPVs coupled to the reaction coordinate, which effect electronic alterations in the active site by pushing the three oxygen centers together in proximity, and accelerate substrate turnover in the phosphorolysis reaction catalyzed by hPNP.

Animals↗
PubMed2004

Computational analysis of the transient movement of helices in sensory rhodopsin II.

MD simulation of sensory rhodopsin II was executed for three intermediates (ground-state, K-state, M-state) appearing in its photocycle. We observed a large displacement of the cytoplasmic side of helixF only in M-state among the three intermediates. This displacement was transmitted to TM2, and the cytoplasmic side of TM2 rotated clockwise. These transient movements are in agreement with the results of an EPR experiment. That is, the early stage of signal transduction in a sRII-HtrII complex was successfully reproduced by the in silico MD simulation. By analyzing the structure of the sRII-HtrII complex, the following findings about the photocycle of sRII were obtained: (1) The hydrogen bonds between helixF and other helices determine the direction of the movement of helixF; (2) three amino acids (Arg162, Thr189, Tyr199) are essential for sRII-HtrII binding and contribute to the motion transfer from sRII to HtrII; (3) after the isomerization of retinal, a major conformational change of retinal was caused by proton transfer from Schiff base to Asp75, which, in turn, triggers the steric collision of retinal with Trp171. This is the main reason for the movement of the cytoplasmic side of helixF.

Archaeal Proteins↗
PubMed1997

Cyclic HIV-1 protease inhibitors derived from mannitol: synthesis, inhibitory potencies, and computational predictions of binding affinities.

Ten C2-symmetric cyclic urea and sulfamide derivatives have been synthesized from L-mannonic gamma-lactone and D-mannitol. The results of experimental measurement of their inhibitory potencies against HIV-1 protease were compared to calculated free energies of binding derived from molecular dynamics (MD) simulations. The compounds were selected, firstly, to enable elucidation of the role of stereochemistry for binding affinity (1a-d) and, secondly, to allow evaluation of the effects of variation in the link to the P1 and P1' phenyl groups on affinity (1a and 2-5). Thirdly, compounds with hydrogen bond-accepting or-donating groups attached to the phenyl groups in the P2 and P2' side chains (6 and 7) were selected. Binding free energies were estimated by a linear response method, whose predictive power for estimating binding affinities from MD simulations was demonstrated.

Anti-HIV Agents↗
PubMed2005

Solvent decompositions and physical properties of decomposition compounds in Li-ion battery electrolytes studied by DFT calculations and molecular dynamics simulations.

The density functional theory (DFT) calculations have been performed for the reduction decompositions of solvents widely used in Li-ion secondary battery electrolytes, ethylene carbonate (EC), propylene carbonate (PC), dimethyl carbonates (DMC), ethyl methyl carbonate (EMC), and diethyl carbonate (DEC), including a typical electrolyte additive, vinylene carbonate (VC), at the level of B3LYP/6-311+G(2d,p), both in the gas phase and solution using the polarizable conductor calculation model. In the gas phase, the first electron reduction for the cyclic carbonates and for the linear carbonates is found to be exothermic and endothermic, respectively, while the second electron reduction is endothermic for all the compounds examined. On the contrary, in solution both first and second electron reductions are exothermic for all the compounds. Among the solvents and the additive examined, the likelihood of undergoing the first electron reduction in solution was found in the order of EC > PC > VC > DMC > EMC > DEC with EC being the most likely reduced. VC, on the other hand, is most likely to undergo the second electron reduction among the compounds, in the order of VC > EC > PC. Based on the results, the experimentally demonstrated effectiveness of VC as an excellent electrolyte additive was discussed. The bulk thermodynamic properties of two dilithium alkylene glycol dicarbonates, dilithium ethylene glycol dicarbonate (Li-EDC) and dilithium 1,2-propylene glycol dicarbonate (Li-PDC), as the major component of solid-electrolyte interface (SEI) films were also examined through molecular dynamics (MD) simulations in order to understand the stability of the SEI film. It was found that film produced from a decomposition of EC, modeled by Li-EDC, has a higher density, more cohesive energy, and less solubility to the solvent than the film produced from decomposition of PC, Li-PDC. Further, MD simulations of the interface between the decomposition compound and graphite suggested that Li-EDC has more favorable interactions with the graphite surface than Li-PDC. The difference in the SEI film stability and the behavior of Li-ion battery cycling among the solvents were discussed in terms of the molecular structures.

Journal Article↗
PubMed1998

Inherent DNA curvature and flexibility correlate with TATA box functionality.

Four 1.5 ns molecular dynamics (MD) simulations were performed on the d(GCTATAAAAGGG).d(CCCTTTTATAGC) double helix dodecamer bearing the Adenovirus major late promoter TATA element and three iso-composition mutants for which physical and biochemical data are available from the same laboratory. Three of these DNA sequences experimentally induce tight binding with the TATA box binding protein (TBP) and induce high transcription rates; the other DNA sequence induces much lower TBP binding and transcription. The x-ray crystal structures have previously shown that the duplex DNA in DNA-TBP complexes are highly bent. We performed and analyzed MD simulations for these four DNAs, whose experimental structures are not available, in order to address the issue of whether inherent DNA structure and flexibility play a role in establishing these observed preferences. A comparison of the experimental and simulated results demonstrated that DNA duplex sequence-dependent curvature and flexibility play a significant role in TBP recognition, binding, and transcriptional activation.

Adenoviridae↗
PubMed2003

Homology modeling of the central catalytic domain of insertion sequence ISLC3 isolated from Lactobacillus casei ATCC 393.

The tertiary structure of the central catalytic domain of insertion sequence ISLC3 isolated from Lactobacillus casei ATCC 393 was predicted using the homology modeling approach. The novel insertion sequence was isolated by us from the template bacteriophage phiA3 of L.casei ATCC 393. The number of amino acid residues of the ISLC3 central catalytic domain was 116 and was treated as the query sequence. There were five Web-available threading methods used to find some primary structure templates for the query sequence. These primary templates were further screened using the SWISS-MODEL Protein Modeling Server and the default parameter settings therein to give six final structure templates. All of these final structure templates were the integrase (IN) protein of retroviruses. Multiple sequence alignment using these IN sequences against the query one revealed the signature DDE motif. Based on the structures of these final templates, the structure of the query sequence was constructed using the InsightII/Discover/Homology programs. A metal ion, Mg(2+), was inserted into the center of the putative catalytic pocket formed by the DDE residues of the predicted structure in the final rounds of refinement by molecular dynamics (MD) simulations. The structure with a metal ion included was designated with Mg and that without a metal ion was designated free Mg. The average exposed surface area of some hydrophobic residues of both the predicted free Mg and with Mg structures were computed and compared with those computed for the six structure templates. Whereas the predicted with Mg structure was slightly more exposed than the predicted free Mg structure, the former appeared to be more stable than the latter, as revealed by the lower conformation energy recorded for the former during the structure refinement by MD simulations. To verify further the predicted structures, the coordinates of both predicted structures were fed into the ERRAT Protein Verification Server. It was found that the quality of the predicted with Mg structure was much better than that of the free Mg structure. The validation results also indicated that regions of the predicted with Mg structure that can be rejected at the 95% confidence level were approximately 20% whereas those which can be rejected at the same level for the six structure templates were approximately 10%. The predicted with Mg structure was also docked into a short oligonucleotide representing the substrate of the ISLC3 transposase using the DOCK_4.0.2 program. It was found that both Glu140 and Asp68 residues of the DDE motif of the predicted with Mg structure were able to form hydrogen bonds with the DNA substrate, which was similar to what was observed in a docking study using the retrovirus IN 1asu and its DNA substrate.

Amino Acid Sequence↗
PubMed2004

Constant-pH molecular dynamics study of protonation-structure relationship in a heptapeptide derived from ovomucoid third domain.

Molecular dynamics (MD) simulations with implicit solvent and variable protonation states for titratable residues at constant pH are performed for a short peptide derived from ovomucoid third domain (OMTKY3), acetyl-Ser-Asp-Asn-Lys-Thr-Tyr-Gly-amide (residues 26-32 of OMTKY3). Nuclear magnetic resonance (NMR) measurements indicate that the pK(a) for Asp is 3.6. However, if the charge on Lys is neutralized by acetylation, then the pK(a) for Asp is 4.0. These pK(a)'s, and therefore the Asp-Lys interaction, are insensitive to changes in ionic strength. The constant-pH MD simulations for both variants of the heptapeptide yield Asp pK(a) values that are 0.6-0.9 pH units greater than experimental values, but the difference between the variants that is observed in the NMR experiments is reproduced much better. Moreover, the simulations suggest that Asp-Lys interactions do not dominate the behavior of this heptapeptide, even for normal Lys residue where there is a possibility of forming a salt bridge between negatively charged Asp and positively charged Lys. This is consistent with the experimentally observed independence of Asp pK(a) values with respect to ionic strength. Another important result of the simulations with variable protonation states is that they lead to ensembles of the heptapeptide structures that are different from those derived from simulations with fixed protonation states. It should be stressed that these results are for structures generated entirely by computer simulations without any restrictions imposed by experimental data.

Algorithms↗