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

Molecular simulations for the conformational assessment of a porphyrin-fullerene dyad in different environments.

Conformational space of a porphyrin-fullerene dyad with the donor and acceptor connected by a relatively flexible linker is studied by molecular dynamics simulations in both non-polar and polar solvents, as well as in vacuum. The most probable conformations obtained from the vacuum MD simulations were optimized with semi-empirical (SE) and density functional theory (DFT) methods and the extent of the structural changes is assessed. The computational results indicate the co-existence of different conformers in both polar and nonpolar solvents showing agreement with experimental results. The most probable vacuum conformations at 300 K are similar to the ones at 0 K, while the structures most often observed in the solvents show less compact conformations. Optimization with SE and DFT calculations leads to structures, which represent relatively well the folded conformations in solvent, which validates the electronic structure calculations relevant to describing photoinduced electron-transfer in H2P-O34-C60.

Computer Simulation↗
PubMed2006

Simulations of a protein translocation pore: SecY.

SecY is the central channel protein of the SecYEbeta translocon, the structure of which has been determined by X-ray diffraction. Extended (15 ns) MD simulations of the isolated SecY protein in a phospholipid bilayer have been performed to explore the relationship between protein flexibility and the mechanisms of channel gating. In particular, principal components analysis of the simulation trajectory has been used to probe the intrinsic flexibility of the isolated SecY protein in the absence of the gamma-subunit (SecE) clamp. Analysis and visualization of the principal eigenvectors support a "plug and clamshell" model of SecY channel gating. The simulation results also indicate that hydrophobic gating at the central pore ring prevents leakage of water and ions through the channel in the absence of a translocating peptide.

Archaeal Proteins↗
PubMed2005

Absolute and relative entropies from computer simulation with applications to ligand binding.

A comparison between two related methods, Schlitter's formula and quasiharmonic analysis, for calculating absolute entropies from the covariance matrix of atomic fluctuations using molecular dynamics (MD) simulations is presented. Calculations for a set of organic compounds in the gas phase are compared to the corresponding statistical thermodynamics results for translational and rotational entropies and to experimental data for vibrational entropies. Encouraging agreement is obtained for translational entropies, but for the rotational contribution, both methods fail to reproduce the theoretically calculated values. Absolute and relative vibrational entropies are found to be better reproduced using quasiharmonic analysis compared to Schlitter's formula. For rotational entropies, we propose a method based on the variances in Euler angles, which gives good agreement with theory. Alternative methods for estimating translational entropies based on principal root mean-square (rms) fluctuations of the center of mass are also presented, and these reproduce theoretically calculated values well. These methodologies are applied to the binding of benzene to T4-lysozyme, where close agreement with the literature is obtained for translational and rotational entropies.

Bacteriophage T4↗
PubMed2002

A procedure for refining a coiled coil protein structure using x-ray fiber diffraction and modeling.

We describe a combined use of experimental and simulation techniques to configure side chains in a coiled coil structure. As already demonstrated in a previous work, x-ray diffraction patterns from hard alpha-keratin fibers in the 5.15 A meridian zone reflect the global configuration of the chi(1) dihedral angle of the coiled coil side chains. Molecular simulations, such as energy minimization and molecular dynamics, and rotameric representation in the PDB, are used here on a heterodimeric coiled coil to investigate the dihedral angle distribution along the sequence. Different procedures have been used to build the structure, the quality assessment was based on the agreement between the simulated diffraction patterns and the experimental ones in the fingerprint region of coiled coils (5.15 A). The best one for building a realistic coiled coil structure consists of placing the side chains using molecular dynamics (MD) simulations, followed by side chain positioning using SMD or SCWRL procedures. The side chains and the backbone are equilibrated during the MD until they reach an equilibrium state for the t/g(+) ratio. Positioning the side chains on the resulting backbone, using the above procedures, gives rise to a well-defined 5.15 A meridian reflection.

Algorithms↗
PubMed2001

Finite-size effects in molecular dynamics simulations: intermediate scattering function and velocity of sound. III. Theory and application to a model krypton fluid.

We describe a method for obtaining the intermediate scattering function I(Q,t) from a computer simulation: it is an extension of our earlier calculation [Salacuse, Denton, and Egelstaff, Phys. Rev. E 53, 2382 (1996)] for the t-->0 limit. We use this approach to obtain I(Q,t) for low Q and t from molecular dynamics (MD) simulations of a model krypton fluid whose atoms interact via a truncated Aziz pair potential, and the results are compared over their range of validity to I(Q,t) determined by the standard MD method and also by a time expansion approach. In its range of validity our approach is much more efficient than the standard MD method; however, it covers a restricted range of t due to the movement of density fluctuations (sound waves) through the simulated fluid which produces an anomaly in the time behavior of I(Q,t). By analyzing I(Q=0,t) the velocity of sound in the simulation is determined, and the results compare favorably with published experimental results for the sound velocity of liquid krypton.

Journal Article↗
PubMed1997

Secondary structure of an isolated P-region from the voltage-gated sodium channel: a molecular modelling/dynamics study.

Conformational studies of synthetic peptides corresponding to the pore-forming regions of voltage-gated sodium channels show a high tendency for beta-sheet conformation when interacting with lipid vesicles, as revealed by circular dichroism and infrared spectroscopy. These observations have guided our choice of possible molecular models for the P-region peptide of domain II of voltage-gated sodium channels: three alternative beta-hairpins, with differing turn assignments, or an alpha-helical hairpin. After generation of models by distance geometry-based methods, molecular dynamics (MD) simulations were run. in the absence of explicit solvent molecules but employing three different dielectric constants, to explore possible conformational preferences. The simulations in the different dielectric environments suggest that a 4-residue turn with the sequence LCGE yields more stable beta-hairpins. The MD results suggest that the SS1 part of the peptide may be more stable as an alpha-helix, whereas the SS2 part tends to adopt a beta-conformation.

Journal Article↗
PubMed2002

Molecular dynamics simulation of proton transport through the influenza A virus M2 channel.

The structural and dynamical properties of a solvated proton in the influenza A virus M2 channel are studied using a molecular dynamics (MD) simulation technique. The second-generation multi-state empirical valence bond (MS-EVB2) model was used to describe the interaction between the excess proton and the channel environment. Solvation structures of the excess proton and its mobility characteristics along the channel were determined. It was found that the excess proton is capable of crossing the channel gate formed by the ring of four histidine residues even though the gate was only partially open. Although the hydronium ion itself did not cross the channel gate by traditional diffusion, the excess proton was able to transport through the ring of histidine residues by hopping between two water molecules located at the opposite sides of the gate. Our data also indicate that the proton diffusion through the channel may be correlated with the changes in channel conformations. To validate this observation, a separate simulation of the proton in a "frozen" channel has been conducted, which showed that the proton mobility becomes inhibited.

Biological Transport↗
PubMed1994

Kinked structures of isolated nicotinic receptor M2 helices: a molecular dynamics study.

The pore-lining M2 helix of the nicotinic acetylcholine receptor exhibits a pronounced kink when the corresponding ion channel is in a closed conformation [N. Unwin (1993) Journal of Molecular Biology, Vol. 229, pp. 1101-1124]. We have performed molecular dynamics simulations of isolated 22-residue M2 helices in order to identify a possible molecular origin of this kink. In order to sample a wide range of conformational space, a simulated annealing protocol was used to generate five initial M2 helix structures, each of which was subsequently used as the basis of 300 ps MD simulations. Two helix sequences (M2 alpha and M2 delta) were studied in this manner, resulting in a total of ten 300 ps trajectories. Kinked helices present in the trajectories were identified and energy minimized to yield a total of five different stable kinked structures. For comparison, a similar molecular dynamics simulation of a Leu23 helix yielded no stable kinked structures. In four of the five kinked helices, the kink was stabilized by H bonds between the helix backbone and polar side-chain atoms. Comparison with data from the literature on site-directed mutagenesis of M2 residues suggests that such polar side-chain to main-chain H bonds may also contribute to kinking of M2 helices in the intact channel protein.

Amino Acid Sequence↗
PubMed2003

Comparison of various implicit solvent models in molecular dynamics simulations of immunoglobulin G light chain dimer.

The present study tests performance of different solvation models applied to molecular dynamics simulation of a large, dimeric protein molecule. Analytical Continuum Electrostatics (ACE) with two different parameter sets, older V98 and new V01, and Effective Energy Function (EEF) are employed in molecular dynamics simulation of immunoglobulin G (IgG) light chain dimer and variable domain of IgG light chain. Results are compared with explicit solvent and distance dependent dielectric constant (DDE) calculations. The overall analysis shows that the EEF method yields results comparable to explicit solvent simulations; however, the stability of simulations is lower. On the other hand, the ACE_V98 model does not seem to achieve the accuracy or stability expected in nanosecond timescale MD simulation for the studied systems. The ACE_V01 model greatly improves stability of the calculation; nonetheless, changes in radius of gyration and solvent accessible surface of the studied systems may indicate that the parameter set still needs to be improved if the method is supposed to be used for simulations of large, polymeric proteins. Additionally, electrostatic contribution to the solvation free energy calculated in the ACE model is compared with a numerical treatment of the dielectric continuum model. Wall clock time of all simulations is compared. It shows that EEF calculation is six times faster than corresponding ACE and 50 times faster than explicit solvent simulations.

Computer Simulation↗
PubMed2004

Beta-hairpin conformation of fibrillogenic peptides: structure and alpha-beta transition mechanism revealed by molecular dynamics simulations.

Understanding the conformational transitions that trigger the aggregation and amyloidogenesis of otherwise soluble peptides at atomic resolution is of fundamental relevance for the design of effective therapeutic agents against amyloid-related disorders. In the present study the transition from ideal alpha-helical to beta-hairpin conformations is revealed by long timescale molecular dynamics simulations in explicit water solvent, for two well-known amyloidogenic peptides: the H1 peptide from prion protein and the Abeta(12-28) fragment from the Abeta(1-42) peptide responsible for Alzheimer's disease. The simulations highlight the unfolding of alpha-helices, followed by the formation of bent conformations and a final convergence to ordered in register beta-hairpin conformations. The beta-hairpins observed, despite different sequences, exhibit a common dynamic behavior and the presence of a peculiar pattern of the hydrophobic side-chains, in particular in the region of the turns. These observations hint at a possible common aggregation mechanism for the onset of different amyloid diseases and a common mechanism in the transition to the beta-hairpin structures. Furthermore the simulations presented herein evidence the stabilization of the alpha-helical conformations induced by the presence of an organic fluorinated cosolvent. The results of MD simulation in 2,2,2-trifluoroethanol (TFE)/water mixture provide further evidence that the peptide coating effect of TFE molecules is responsible for the stabilization of the soluble helical conformation.

Animals↗
PubMed2004

OmpT: molecular dynamics simulations of an outer membrane enzyme.

Five molecular dynamics simulations (total duration >25 ns) have been performed on the Escherichia coli outer membrane protease OmpT embedded in a dimyristoylphosphatidylcholine lipid bilayer. Globally the protein is conformationally stable. Some degree of tilt of the beta-barrel is observed relative to the bilayer plane. The greatest degree of conformational flexibility is seen in the extracellular loops. A complex network of fluctuating H-bonds is formed between the active site residues, such that the Asp210-His212 interaction is maintained throughout, whereas His212 and Asp83 are often bridged by a water molecule. This supports a catalytic mechanism whereby Asp83 and His212 bind a water molecule that attacks the peptide carbonyl. A configuration yielded by docking calculations of OmpT simulation snapshots and a model substrate peptide Ala-Arg-Arg-Ala was used as the starting point for an extended Huckel calculation on the docked peptide. These placed the lowest unoccupied molecular orbital mainly on the carbon atom of the central C=O in the scissile peptide bond, thus favoring attack on the central peptide by the water held by residues Asp83 and His212. The trajectories of water molecules reveal exchange of waters between the intracellular face of the membrane and the interior of the barrel but no exchange at the extracellular mouth. This suggests that the pore-like region in the center of OmpT may enable access of water to the active site from below. The simulations appear to reveal the presence of specific lipid interaction sites on the surface of the OmpT barrel. This reveals the ability of extended MD simulations to provide meaningful information on protein-lipid interactions.

Binding Sites↗
PubMed2005

Reorganization and conformational changes in the reduction of tetraheme cytochromes.

Molecular dynamics simulation (MD) constitutes an alternative to time-consuming experiments for studying conformational changes. We apply MD on a redox system where experimental information exists for the fully oxidized and fully reduced states: tetraheme cytochrome c3. Instead of doing one simulation for each state, we apply 10 4-ns replicas for both states, which provides robust statistics to characterize the redox changes. Besides these long simulations, we perform 120 short ones (50 ps), where an equilibrated oxidized state is perturbed to a reduced state. This allows the application of a nonequilibrium method, the subtraction technique, which makes it possible to characterize the different timescales of conformational changes. Reduction induces conformational changes in the N-terminus and on the loops spanning residues 36-42 and 88-93, which correlate very well with experiments, demonstrating the applicability of this methodology. We also analyze the effect of reduction on hydrogen bonds, solvent accessible surface and bound water, the changes being found to involve the hemes and propionate groups. Redox-induced protonation is also investigated, by protonating the propionates D from hemes I and IV. Although this change in the former does not have major conformational consequences, it induces in the latter conformational changes beyond the ones obtained with reduction.

Computer Simulation↗
PubMed2004

Validation of the GROMOS force-field parameter set 45Alpha3 against nuclear magnetic resonance data of hen egg lysozyme.

The quality of molecular dynamics (MD) simulations of proteins depends critically on the biomolecular force field that is used. Such force fields are defined by force-field parameter sets, which are generally determined and improved through calibration of properties of small molecules against experimental or theoretical data. By application to large molecules such as proteins, a new force-field parameter set can be validated. We report two 3.5 ns molecular dynamics simulations of hen egg white lysozyme in water applying the widely used GROMOS force-field parameter set 43Alpha1 and a new set 45Alpha3. The two MD ensembles are evaluated against NMR spectroscopic data NOE atom-atom distance bounds, (3)J(NHalpha) and (3)J(alphabeta) coupling constants, and (15)N relaxation data. It is shown that the two sets reproduce structural properties about equally well. The 45Alpha3 ensemble fulfills the atom-atom distance bounds derived from NMR spectroscopy slightly less well than the 43Alpha1 ensemble, with most of the NOE distance violations in both ensembles involving residues located in loops or flexible regions of the protein. Convergence patterns are very similar in both simulations atom-positional root-mean-square differences (RMSD) with respect to the X-ray and NMR model structures and NOE inter-proton distances converge within 1.0-1.5 ns while backbone (3)J(HNalpha)-coupling constants and (1)H-(15)N order parameters take slightly longer, 1.0-2.0 ns. As expected, side-chain (3)J(alphabeta)-coupling constants and (1)H-(15)N order parameters do not reach full convergence for all residues in the time period simulated. This is particularly noticeable for side chains which display rare structural transitions. When comparing each simulation trajectory with an older and a newer set of experimental NOE data on lysozyme, it is found that the newer, larger, set of experimental data agrees as well with each of the simulations. In other words, the experimental data converged towards the theoretical result.

Animals↗
PubMed1994

A new method for predicting binding affinity in computer-aided drug design.

A new semi-empirical method for calculating free energies of binding from molecular dynamics (MD) simulations is presented. It is based on standard thermodynamic cycles and on a linear approximation of polar and non-polar free energy contributions from the corresponding MD averages. The method is tested on a set of endothiapepsin inhibitors and found to give accurate results both for absolute as well as relative free energies.

Aspartic Acid Endopeptidases↗
PubMed2003

Structure and dynamics of metal ions in solution: QM/MM molecular dynamics simulations of Mn(2+) and V(2+).

Structural and dynamical properties of the transition metal ions V(2+) and Mn(2+) in aqueous solution, resulting from combined quantum mechanical (QM)/molecular mechanical (MM) molecular dynamics (MD) simulations have been compared. The necessity of polarization functions on the ligand's oxygen for a satisfactory description of such ions in aqueous solution is shown using V(2+) as test case. Radial distribution functions, coordination number distributions, and several angle distributions were pursued for a detailed structural comparison of the first hydration shells. Dynamical properties, such as the librational and vibrational motions of water molecules were evaluated by means of velocity autocorrelation functions. Approximative normal coordinate analyses were employed to calculate the rotational frequencies and vibrational motions around the three principal axes. The very low exchange rates for the first shell water exchanges only allow an investigation of the water exchange processes in the second shell, which take place within the picosecond range.

Cations, Divalent↗
PubMed2006

Assessment of two theoretical methods to estimate potentiometric titration curves of peptides: comparison with experiment.

We compared the ability of two theoretical methods of pH-dependent conformational calculations to reproduce experimental potentiometric titration curves of two models of peptides: Ac-K5-NHMe in 95% methanol (MeOH)/5% water mixture and Ac-XX(A)7OO-NH2 (XAO) (where X is diaminobutyric acid, A is alanine, and O is ornithine) in water, methanol (MeOH), and dimethyl sulfoxide (DMSO), respectively. The titration curve of the former was taken from the literature, and the curve of the latter was determined in this work. The first theoretical method involves a conformational search using the electrostatically driven Monte Carlo (EDMC) method with a low-cost energy function (ECEPP/3 plus the SRFOPT surface-solvation model, assumming that all titratable groups are uncharged) and subsequent reevaluation of the free energy at a given pH with the Poisson-Boltzmann equation, considering variable protonation states. In the second procedure, molecular dynamics (MD) simulations are run with the AMBER force field and the generalized Born model of electrostatic solvation, and the protonation states are sampled during constant-pH MD runs. In all three solvents, the first pKa of XAO is strongly downshifted compared to the value for the reference compounds (ethylamine and propylamine, respectively); the water and methanol curves have one, and the DMSO curve has two jumps characteristic of remarkable differences in the dissociation constants of acidic groups. The predicted titration curves of Ac-K5-NHMe are in good agreement with the experimental ones; better agreement is achieved with the MD-based method. The titration curves of XAO in methanol and DMSO, calculated using the MD-based approach, trace the shape of the experimental curves, reproducing the pH jump, while those calculated with the EDMC-based approach and the titration curve in water calculated using the MD-based approach have smooth shapes characteristic of the titration of weak multifunctional acids with small differences between the dissociation constants. Nevertheless, quantitative agreement between theoretically predicted and experimental titration curves is not achieved in all three solvents even with the MD-based approach, which is manifested by a smaller pH range of the calculated titration curves with respect to the experimental curves. The poorer agreement obtained for water than for the nonaqueous solvents suggests a significant role of specific solvation in water, which cannot be accounted for by the mean-field solvation models.

Algorithms↗
PubMed2004

Pivotal role of Gly 121 in dihydrofolate reductase from Escherichia coli: the altered structure of a mutant enzyme may form the basis of its diminished catalytic performance.

The structure and folding of dihydrofolate reductase (DHFR) from Escherichia coli and the mutant G121V-DHFR, in which glycine 121 in the exterior FG loop was replaced with valine, were studied by molecular dynamics simulations and CD and fluorescence spectroscopy. The importance of residue 121 for the chemical step during DHFR catalysis had been demonstrated previously. High-temperature MD simulations indicated that while DHFR and G121V-DHFR followed similar unfolding pathways, the strong contacts between the M20 loop and the FG loop in DHFR were less stable in the mutant. These contacts have been proposed to be involved in a coupled network of interactions that influence the protein dynamics and promote catalysis [Benkovic, S. J., and Hammes-Schiffer, S. (2003) Science 301, 1196-1202]. CD spectroscopy of DHFR and G121V-DHFR indicated that the two proteins existed in different conformations at room temperature. While the thermally induced unfolding of DHFR was highly cooperative with a midpoint at 51.6 +/- 0.7 degrees C, G121V-DHFR exhibited a gradual decrease in its level of secondary structure without a clear melting temperature. Temperature-induced unfolding and renaturation from the urea-denatured state revealed that both proteins folded via highly fluorescent intermediates. The formation of these intermediates occurred with relaxation times of 149 +/- 4.5 and 256 +/- 13 ms for DHFR and G121V-DHFR, respectively. The fluorescence intensity for the intermediates formed during refolding of G121V-DHFR was approximately twice that of the wild-type. While the fluorescence intensity then slowly decayed for DHFR toward a state representing the native protein, G121V-DHFR appeared to be trapped in a highly fluorescent state. These results suggest that the reduced catalytic activity of G121V-DHFR is the consequence of nonlocal structural effects that may result in a perturbation of the network of promoting motions.

Catalysis↗
PubMed1991

Conformational study of cyclolinopeptide A. A distance geometry and molecular dynamics approach.

The conformation of cyclolinopeptide A, c(Pro-Pro-Phe-Phe-Leu-Ile-Ile-Leu-Val), a naturally occurring peptide with remarkable cytoprotective activity, has been investigated by means of distance geometry calculations and molecular dynamics simulations. The starting points for all the calculations were an X-ray structure and other structures obtained from distance geometry calculations based on NMR data. Restrained and unrestrained molecular dynamics simulations are reported in vacuo and in CCl4. Structural and dynamic properties are investigated and compared with those experimentally determined. The conformation obtained from the MD simulations which best reproduces the NMR parameters is at the same time one of the most stable ones and is also fairly similar to the crystal structure. An explanation for the occurrence of multiple conformations in solution at room temperature is given.

Amino Acid Sequence↗