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

Phi-value analysis by molecular dynamics simulations of reversible folding.

In Phi-value analysis, the effects of mutations on the folding kinetics are compared with the corresponding effects on thermodynamic stability to investigate the structure of the protein-folding transition state (TS). Here, molecular dynamics (MD) simulations (totaling 0.65 ms) have been performed for a large set of single-point mutants of a 20-residue three-stranded antiparallel beta-sheet peptide. Between 57 and 120 folding events were sampled at near equilibrium for each mutant, allowing for accurate estimates of folding/unfolding rates and stability changes. The Phi values calculated from folding and unfolding rates extracted from the MD trajectories are reliable if the stability loss upon mutation is larger than approximately 0.6 kcal/mol, which is observed for 8 of the 32 single-point mutants. The same heterogeneity of the TS of the wild type was found in the mutated peptides, showing two possible pathways for folding. Single-point mutations can induce significant TS shifts not always detected by Phi-value analysis. Specific nonnative interactions at the TS were observed in most of the peptides studied here. The interpretation of Phi values based on the ratio of atomic contacts at the TS over the native state, which has been used in the past in MD and Monte Carlo simulations, is in agreement with the TS structures of wild-type peptide. However, Phi values tend to overestimate the nativeness of the TS ensemble, when interpreted neglecting the nonnative interactions.

Computer Simulation↗
PubMed2005

A molecular dynamics study on universal properties of polymer chains in different solvent qualities. Part I. A review of linear chain properties.

This paper investigates the conformational and scaling properties of long linear polymer chains. These investigations are done with the aid of Monte Carlo (MC) and molecular dynamics (MD) simulations. Chain lengths that comprise several orders of magnitude to reduce errors of finite size scaling, including the effect of solvent quality, ranging from the athermal limit over the theta-transition to the collapsed state of chains are investigated. Also the effect of polydispersity on linear chains is included which is an important issue in the real fabrication of polymers. A detailed account of the hybrid MD and MC simulation model and the exploited numerical methods is given. Many results of chain properties in the extrapolated limit of infinite chain lengths are documented and universal properties of the chains within their universality class are given. An example of the difference between scaling exponents observed in actual solvents and those observed in the extremes of "good solvents" and "theta-solvents" in simulations is provided by comparing simulation results with experimental data on low density polyethylene. This paper is concluded with an outlook on the extension of this study to branched chain systems of many different branching types.

Journal Article↗
PubMed2026

Stable simulations do not guarantee functional engagement: a case study of off-target prediction for Seladelpar and Zanamivir.

Identifying off-target interactions of approved drugs is important to anticipate side effects and uncover repurposing opportunities. Computational pipelines combining structural homology, structure prediction, and molecular dynamics (MD) simulations offer a promising strategy, but it remains unclear whether stable, control-like MD trajectories reliably indicate functional engagement. We examined this in a case study of two approved drugs. Using the Evolutionary Classification of Protein Domains (ECOD) framework to select candidate off-targets, we modeled each drug-protein complex as two independent AlphaFold3 models and simulated both by MD, for Seladelpar (a PPARδ agonist) and Zanamivir, an influenza neuraminidase inhibitor that also inhibits human Sialidase-2 (NEU2). Candidates were ranked by the similarity of global MD descriptors to the on-target control. For Seladelpar, the three top-ranked candidates (FXR, RARγ, ERRγ) were tested experimentally; the Zanamivir set was analyzed computationally only. None showed measurable activity in reporter or thermal shift assays, despite stable simulations and descriptor values comparable to the control. Including PPARα and PPARγ as weak-positive comparators, these descriptors did not rank genuine interactions closer to the control than inactive candidates. Residue-level comparison with experimental structures showed the predicted poses reproduced only part of the canonical contacts. Where experimental drug-bound structures existed, AlphaFold3 reproduced the pose for PPARα but not PPARγ, and its per-model confidence did not track pose accuracy. Within this case study, the specific global descriptors examined reflect complex stability rather than functional engagement, which does not mean MD-based approaches cannot make this distinction.

Zanamivir↗
PubMed2001

Finite-size dependence of the bridge function extracted from molecular dynamics simulations.

The bridge function for liquid sodium at T=373 K is obtained by using the mean spherical approximation to extrapolate the pair distribution function (PDF), calculated in molecular dynamics (MD) simulations, beyond the half simulation box length for two sizes of the MD system. The bridge function is found to strongly depend on the total number of particles used in the simulation cell. This dependency leads to a spurious maximum of the static structure factor at long wavelengths, obtained from the reference hypernetted-chain approximation (RHNC) with the MD system used as a reference system (RHNC-MD). A simple self-consistent procedure, proposed to account for the finite-size effects in the bridge function, allows one to efficiently correct the RHNC-MD static structure factor for all unphysical manifestations.

Journal Article↗
PubMed2003

Insights from molecular dynamics simulations into pH-dependent enantioselective hydrolysis of ibuprofen esters by Candida rugosa lipase.

An interesting observation was found during our continued studies on the hydrolysis of ibuprofen esters by Candida rugosa lipase (CRL). An important role is played by pH in the stereospecific hydrolysis of these esters. The flap region of CRL plays a significant role in the access of the substrate to the active site of the enzyme. At pH 5.6, 48% of the methyl ester and 5% of the butyl ester of ibuprofen were hydrolysed in 5.5 h, whereas at pH 7.2, 9% of methyl ester and 45% of the butyl ester of ibuprofen was hydrolysed in a identical reaction time using CRL. This lead us to assume that CRL prefers the methyl ester of ibuprofen as a substrate at an acidic pH and the butyl ester of ibuprofen at a neutral pH. Therefore, in order to understand the role of pH in the substrate selection by CRL for the esters of ibuprofen we used the crystallographic coordinates of the open form of the CRL (1CRL) for molecular dynamics (MD) simulations under acidic and neutral conditions for 2 ns using GROMACS. The final structures obtained after simulation in acidic and neutral conditions were compared with the energy-minimized structure, and the root-mean-square deviations (r.m.s.ds) were calculated. The r.m.s.d. of the CRL flap at neutral pH was found to be greater than that of the CRL flap at acidic pH. The extent to which the flap opens at neutral pH allowed the bulkier substrate, the butyl ester of ibuprofen, to diffuse into the active site and provides the best enzyme-substrate fit for this specific substrate. At acidic pH there is a decreased opening of the flap thereby accommodating a more compact substrate, namely the methyl ester of ibuprofen. Thus, simulation experiments using MD provide reasonable insight for the pH-dependent substrate selectivity of CRL in aqueous environments.

Amino Acids↗
PubMed1987

Molecular dynamics simulation of crystalline beta-cyclodextrin dodecahydrate at 293 K and 120 K.

Molecular dynamics (MD) simulations for crystalline beta-cyclodextrin dodecahydrate (beta-CD) at two different temperatures, 293 K and 120 K, have been performed using the GROMOS program package. The calculated structural properties are compared to those obtained from neutron diffraction studies of this system at the quoted temperatures. The simulation was carried out over a period of 20 ps on four unit cells containing 8 beta-CD molecules and 96 water molecules, whereby all atoms were allowed to move. At room temperature, the experimental positions of the (non-hydrogen) glucose atoms are reproduced within 0.034 nm, a value which is smaller than the experimental (0.041 nm) or simulated (0.049 nm) overall root mean square (rms) positional fluctuation. The corresponding numbers for the low temperature study are 0.046 nm, 0.019 nm and 0.022 nm. At both temperatures the experimentally observed degree of anisotropy of the atomic motions is also found in the simulations. The comparison of a variety of structural properties leads to the conclusion that the molecular model and force field used are able to simulate the cyclodextrin system very well. Experimentally observed differences in properties as a function of number of glucose units in the CD molecule (alpha-CD, 6 versus beta-CD, 7) and as a function of temperature are qualitatively reproduced by the simulations.

Carbohydrate Conformation↗
PubMed2006

Molecular dynamics simulations on the Escherichia coli ammonia channel protein AmtB: mechanism of ammonia/ammonium transport.

Molecular dynamics (MD) simulations have been performed at the atomic level to study the ammonium/ammonia transport across the Escherichia coli AmtB membrane protein. Although ammonia primarily exists in the form of NH(4)(+) in aqueous solution, the recent X-ray structure determination of AmtB reveals that the ammonium/ammonia transporter proteins are ammonia-conducting channels rather than ammonium ion transporters [Khademi, S.; et al. Science 2004, 305, 1587; Zheng, L.; et al. Proc. Natl. Acad. Sci. U.S.A. 2004, 101, 17090]. Our simulations showed that the entrance of NH(4)(+) into the periplasmic recruitment vestibule requires only 3.1 kcal/mol of energy. This is consistent with the X-ray crystal structure, where one NH(4)(+) is captured in the binding vestibule. In this vestibule, NH(4)(+) loses one water of hydration, but the loss is compensated by a hydrogen bond, first with the backbone carbonyl oxygen of Phe161 then with the hydroxyl group of Ser219, as well as the stabilizing pi-cation interactions with the aromatic rings of Trp148 and Phe107 in the AmtB protein. In the end of this recruitment vestibule, the phenyl ring of Phe107 dynamically switches to an open state. This is correlated with a slight rotation and shifting of the indole ring of Trp148, which eventually creates a slot for the initially buried carboxylate group of Asp160 to become exposed to the bulk solvent. A hydrogen bond wire between NH(4)(+) and the carboxylate group of Asp160 via two water molecules was observed. Thus, Asp160 is most likely the proton acceptor from NH(4)(+). This explains the high conservation of Asp160 in Amt proteins and why the D160A mutant would completely quench the activity of AmtB [Javelle, A.; et al. J. Biol. Chem. 2004, 279, 8530; Marini, A. M.; et al. Curr. Genet. 2006, 49, 364]. Once NH(4)(+) deprotonates, the phenyl ring of Phe215 rotates to open, and the subsequent passage of NH(3) through the channel is straightforward.

Ammonia↗
PubMed2002

Molecular dynamics simulations of a calmodulin-peptide complex in solution.

We have performed an 4-ns MD simulation of calmodulin complexed with a target peptide in explicit water, under realistic conditions of constant temperature and pressure, in the presence of a physiological concentration of counterions and using Ewald summation to avoid truncation of long-range electrostatic forces. During the simulation the system tended to perform small fluctuations around a structure similar to, but somewhat looser than the starting crystal structure. The calmodulin-peptide complex was quite rigid and did not exhibit any large amplitude domain motions such as previously seen in apo- and calcium-bound calmodulin. We analyzed the calmodulin-peptide interactions by calculating buried surface areas, CHARMM interaction energies and continuum model interaction free energies. In the trajectory, the protein surface area buried by contact with the peptide is 1373 A(2) approximately evenly divided between the calmodulin N-terminal, C-terminal and central linker regions. A majority of this buried surface, 803 A(2), comes from nonpolar residues, in contrast to the protein as a whole, for which the surface is made up of mostly polar and charged groups. Our continuum calculations indicate that the largest favorable contribution to peptide binding comes from burial of molecular surface upon complex formation. Electrostatic contributions are favorable but smaller in the trajectory structures, and actually unfavorable for binding in the crystal structure. Since nonpolar groups make up most of buried surface of the protein, our calculations suggest that the hydrophobic effect is the main driving force for binding the helical peptide to calmodulin, consistent with thermodynamic analysis of experimental data. Besides the burial of nonpolar surface area, secondary contributions to peptide binding come from burial of polar surface and electrostatic interactions. In the nonpolar interactions a crucial role is played by the nine methionines of calmodulin. In the electrostatic interactions the negatively charged protein residues and positively charged peptide residues play a dominant role.

Algorithms↗
PubMed2004

Molecular dynamics simulations of the lipid bilayer edge.

Phospholipid bilayers have been intensively studied by molecular dynamics (MD) simulation in recent years. The properties of bilayer edges are important in determining the structure and stability of pores formed in vesicles and biomembranes. In this work, we use molecular dynamics simulation to investigate the structure, dynamics, and line tension of the edges of bilayer ribbons composed of pure dimyristoylphosphatidylcholine (DMPC) or palmitoyl-oleoylphosphatidylethanolamine (POPE). As expected, we observe a significant reorganization of lipids at and near the edges. The treatment of electrostatic effects is shown to have a qualitative impact on the structure and stability of the edge, and significant differences are observed in the dynamics and structure of edges formed by DMPC and palmitoyl-oleoylphosphatidylethanolamine. From the pressure anisotropy in the simulation box, we calculate a line tension of approximately 10-30 pN for the DMPC edge, in qualitative agreement with experimental estimates for similar lipids.

Anisotropy↗
PubMed2005

Reproducible polypeptide folding and structure prediction using molecular dynamics simulations.

The folding of a polypeptide from an extended state to a well-defined conformation is studied using microsecond classical molecular dynamics (MD) simulations and replica exchange molecular dynamics (REMD) simulations in explicit solvent and in vacuo. It is shown that the solvated peptide folds many times in the REMD simulations but only a few times in the conventional simulations. From the folding events in the classical simulations we estimate an approximate folding time of 1-2 micros. The REMD simulations allow enough sampling to deduce a detailed Gibbs free energy landscape in three dimensions. The global minimum of the energy landscape corresponds to the native state of the peptide as determined previously by nuclear magnetic resonance (NMR) experiments. Starting from an extended state it takes about 50 ns before the native structure appears in the REMD simulations, about an order of magnitude faster than conventional MD. The calculated melting curve is in good qualitative agreement with experiment. In vacuo, the peptide collapses rapidly to a conformation that is substantially different from the native state in solvent.

Computer Simulation↗
PubMed2005

Unfolding and extraction of a transmembrane alpha-helical peptide: dynamic force spectroscopy and molecular dynamics simulations.

An atomic force microscope (AFM) was used to visualize CWALP(19)23 peptides ((+)H(3)N-ACAGAWWLALALALALALALWWA-COO(-)) inserted in gel-phase DPPC and DSPC bilayers. The peptides assemble in stable linear structures and domains. A model for the organization of the peptides is given from AFM images and a 20 ns molecular dynamics (MD) simulation. Gold-coated AFM cantilevers were used to extract single peptides from the bilayer through covalent bonding to the cystein residue. Experimental and simulated force curves show two distinct force maxima. In the simulations these two maxima correspond to the extraction of the two pairs of tryptophan residues from the membrane. Unfolding of the peptide precedes extraction of the second distal set of tryptophans. To probe the energies involved, AFM force curves were obtained from 10 to 10(4) nm/s and MD force curves were simulated with 10(8)-10(11) nm/s pulling velocities (V). The velocity relationship with the force, F, was fitted to two fluctuation adhesive potential models. The first assumes the pulling produces a constant bias in the potential and predicts an F approximately ln (V) relationship. The second takes into account the ramped bias that the linker feels as it is being driven out of the adhesion complex and scales as F approximately (ln V)2/3.

1,2-Dipalmitoylphosphatidylcholine↗
PubMed2006

Folding free-energy landscape of a 10-residue mini-protein, chignolin.

Chignolin is an artificial mini-protein composed of 10 residues (GYDPETGTWG) that has been shown to cooperatively fold into a beta-hairpin structure in water. We extensively explored the conformational space of chignolin using a 180-ns multicanonical molecular dynamics (MD) simulation and analyzed its folding free-energy landscape. In the MD trajectory, we found structures that satisfy 99% of the experimental restraints and are quite close to the experimentally determined structures with C(alpha) root-mean-square-deviations of less than 0.5 Angstroms. These structures formed a large cluster in the conformational space with the largest probability of existence, agreeing well with the experiment.

Cluster Analysis↗
PubMed2004

All-atom Monte Carlo simulation of GCAA RNA folding.

We report a detailed all-atom simulation of the folding of the GCAA RNA tetraloop. The GCAA tetraloop motif is a very common and thermodynamically stable secondary structure in natural RNAs. We use our simulation methods to study the folding behavior of a 12-base GCAA tetraloop structure with a four-base helix adjacent to the tetraloop proper. We implement an all-atom Monte Carlo (MC) simulation of RNA structural dynamics using a Go potential. Molecular dynamics (MD) simulation of RNA and protein has realistic energetics and sterics, but is extremely expensive in terms of computational time. By coarsely treating non-covalent energetics, but retaining all-atom sterics and entropic effects, all-atom MC techniques are a useful method for the study of protein and now RNA. We observe a sharp folding transition for this structure, and in simulations at room temperature the state histogram shows three distinct minima: an unfolded state (U), a more narrow intermediated state (I), and a narrow folded state (F). The intermediate consists primarily of structures with the GCAA loop and some helix hydrogen bonds formed. Repeated kinetic folding simulations reveal that the number of helix base-pairs forms a simple 1D reaction coordinate for the I-->N transition.

Base Sequence↗
PubMed2005

Anesthetic interaction with ketosteroid isomerase: insights from molecular dynamics simulations.

The nature and the sites of interactions between anesthetic halothane and homodimeric Delta5-3-ketosteroid isomerase (KSI) are characterized by flexible ligand docking and confirmed by 1H-15N NMR. The dynamics consequence of halothane interaction and the implication of the dynamic changes to KSI function are studied by multiple 5-ns molecular dynamics simulations in the presence and absence of halothane. Both docking and MD simulations show that halothane prefer the amphiphilic dimeric interface to the hydrophobic active site of KSI. Halothane occupancy at the dimer interface disrupted the intersubunit hydrogen bonding formed either directly through side chains of polar residues or indirectly through the mediation of the interfacial water molecules. Moreover, in the presence of halothane, the exchange rate of the bound waters with bulk water was increased. Halothane perturbation to the dimer interface affected the overall flexibility of the active site. This action is likely to contribute to the halothane-induced reduction of the KSI activity. The allosteric halothane modulation of the dynamics-function relationship of KSI without direct competition at the enzymatic active sites may be generalized to offer a unifying explanation of anesthetic action on a diverse range of multidomain neuronal proteins that are potentially relevant to clinical general anesthesia.

Anesthetics, Inhalation↗
PubMed1998

Brownian and essential dynamics studies of the HIV-1 integrase catalytic domain.

The three-dimensional structure of the active site region of the enzyme HIV-1 integrase is not unambiguously known. This region includes a flexible peptide loop that cannot be well resolved in crystallographic determinations. Here we present two different computational approaches with different levels of resolution and on different time-scales to understand this flexibility and to analyze the dynamics of this part of the protein. We have used molecular dynamics simulations with an atomic model to simulate the region in a realistic and reliable way for 1 ns. It is found that parts of the loop wind up after 300 ps to extend an existing helix. This indicates that the helix is longer than in the earlier crystal structures that were used as basis for this study. Very recent crystal data confirms this finding, underlining the predictive value of accurate MD simulations. Essential dynamics analysis of the MD trajectory yields an anharmonic motion of this loop. We have supplemented the MD data with a much lower resolution Brownian dynamics simulation of 600 ns length. It provides ideas about the slow-motion dynamics of the loop. It is found that the loop explores a conformational space much larger than in the MD trajectory, leading to a "gating"-like motion with respect to the active site.

Catalytic Domain↗
PubMed2001

alpha-, beta-, and gamma-Cyclodextrin dimers. Molecular modeling studies by molecular mechanics and molecular dynamics simulations.

The alpha-, beta-, and gamma-cyclodextrin (CyD) dimers have been studied by molecular mechanics (MM) and molecular dynamics (MD) calculations, and the relative stability of dimers and the involved molecular interactions have been determined. Three possible orientations were considered for the alpha-, beta-, and gamma-CyD dimers: the head-to-head, the head-to-tail, and the tail-to-tail. In vacuo MM calculations were used to obtain the most stable arrangements, and MD simulations were performed over all energy minima obtained for each dimer. Results from MD always show head-to-head orientation as the most stable as a result of the larger number of intermolecular hydrogen bonds present.

Journal Article↗
PubMed2005

Cutoff size need not strongly influence molecular dynamics results for solvated polypeptides.

The correct treatment of van der Waals and electrostatic nonbonded interactions in molecular force fields is essential for performing realistic molecular dynamics (MD) simulations of solvated polypeptides. The most computationally tractable treatment of nonbonded interactions in MD utilizes a spherical distance cutoff (typically, 8-12 A) to reduce the number of pairwise interactions. In this work, we assess three spherical atom-based cutoff approaches for use with all-atom explicit solvent MD: abrupt truncation, a CHARMM-style electrostatic shift truncation, and our own force-shifted truncation. The chosen system for this study is an end-capped 17-residue alanine-based alpha-helical peptide, selected because of its use in previous computational and experimental studies. We compare the time-averaged helical content calculated from these MD trajectories with experiment. We also examine the effect of varying the cutoff treatment and distance on energy conservation. We find that the abrupt truncation approach is pathological in its inability to conserve energy. The CHARMM-style shift truncation performs quite well but suffers from energetic instability. On the other hand, the force-shifted spherical cutoff method conserves energy, correctly predicts the experimental helical content, and shows convergence in simulation statistics as the cutoff is increased. This work demonstrates that by using proper and rigorous techniques, it is possible to correctly model polypeptide dynamics in solution with a spherical cutoff. The inherent computational advantage of spherical cutoffs over Ewald summation (and related) techniques is essential in accessing longer MD time scales.

Circular Dichroism↗
PubMed2003

Insights into the acylation mechanism of class A beta-lactamases from molecular dynamics simulations of the TEM-1 enzyme complexed with benzylpenicillin.

Herein, we present results from molecular dynamics MD simulations ( approximately 1 ns) of the TEM-1 beta-lactamase in aqueous solution. Both the free form of the enzyme and its complex with benzylpenicillin were studied. During the simulation of the free enzyme, the conformation of the Omega loop and the interresidue contacts defining the complex H-bond network in the active site were quite stable. Most interestingly, the water molecule connecting Glu166 and Ser70 does not exchange with bulk solvent, emphasizing its structural and catalytic relevance. In the presence of the substrate, Ser130, Ser235, and Arg244 directly interact with the beta-lactam carboxylate via H-bonds, whereas the Lys234 ammonium group has only an electrostatic influence. These interactions together with other specific contacts result in a very short distance ( approximately 3 A) between the attacking hydroxyl group of Ser70 and the beta-lactam ring carbonyl group, which is a favorable orientation for nucleophilic attack. Our simulations also gave insight into the possible pathways for proton abstraction from the Ser70 hydroxyl group. We propose that either the Glu166 carboxylate-Wat1 or the substrate carboxylate-Ser130 moieties could abstract a proton from the nucleophilic Ser70.

Acylation↗