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 865 records · Page 48Linked to original sources
PubMed2004

Mechanism of glucose oxidation by quinoprotein soluble glucose dehydrogenase: insights from molecular dynamics studies.

We have generated 3 ns molecular dynamic (MD) simulations, in aqueous solution, of the bacterial soluble glucose dehydrogenase enzyme.PQQ.glucose complex and intermediates formed in PQQ reduction. In the MD structure of enzyme.PQQ.glucose complex the imidazole of His144 is hydrogen bonded to the hydroxyl hydrogen of H[bond]OC1(H) of glucose. The tightly hydrogen-bonded triad Asp163-His144-glucose (2.70 and 2.91 A) is involved in proton abstraction from glucose concerted with the hydride transfer from the C1[bond]H of glucose to the >C5[double bond]O quinone carbon of PQQ. The reaction is assisted by Arg228 hydrogen bonding to the carbonyl oxygen of >C5[double bond]O. The rearrangement of [bond](H)C5(O-)[bond]C4([double bond]O)[bond] of II to [bond]C5(OH)[double bond]C4(OH)[bond] of PQQH(2) hydroquinone is assisted by general acid protonatation of the >C4[double bond]O oxygen by protonated His144 and hydrogen bonds of Arg228 to the oxyanion O5. The continuous hydrogen bonding of the amide side chain of Asn229 to >C4[double bond]O4 oxygen and that of the O5 oxygen of the cofactor to Wat89 is observed throughout the entire reaction.

Binding Sites↗
PubMed1993

Molecular dynamics simulations of poly(dA).poly(dT): comparisons between implicit and explicit solvent representations.

The program AMBER 3.0 has been used to generate molecular dynamics trajectories of a poly(dA).poly(dT) decamer. The simulations were performed using different methods to treat solvent effects. Results of a simulation including 18 counterions NH4+ and 4109 water molecules under (N, P, T) conditions were compared to simulation runs with implicit solvent representation in which solvent screening effects were represented by the use of a sigmoidal distance-dependent dielectric function. In the latter case, the system was simulated under microcanonical (N, V, E) and canonical (N, V, T) conditions. For the fully hydrated system simulation, a preequilibration protocol was developed since it was observed that long and progressive periods of heating and equilibration on the overall system were necessary in order to avoid energetic collisions between the solute and the solvent molecules, leading to severe irreversible deformation of the solute. A detailed analysis of DNA conformations, sugar puckers, and stability of the hydrogen bonds, Watson-Crick and three-center H bonds, is reported. The results show that DNA remains essentially in the B conformer with a tendency in the hydrated model to adopt a slightly distorted, unwound, and stretched conformation in comparison to standard B-DNA. Concerning sugar puckers, the mean pseudorotation phases of the adenine residues are systematically higher than those of the thymine residues, except in the case of the hydrated model for which a articular behavior is observed for the adenine strand. In this case, the terminal bases oscillate between C2'-endo and O4'-endo and the central ones stay in the C3'-endo domain. The mean lifetimes of the internal Watson-Crick H-bond (A) HN6...O4(T) are also dependent on the base pairs included in the calculation, excepted for the implicit solvent simulation at constant temperature. The three-center H bonds have very small mean lifetimes in all three cases of MD simulation. In the minor groove of the hydrated model, a spine of hydration is found as observed by x-ray crystallography and other theoretical simulations. On the basis of the rms deviations, it appears that the fully hydrated simulation has not reached a plateau at the end of the run, while the implicit simulation at constant energy seems to have converged. At constant temperature, very large oscillations in rms deviations are observed.

Base Sequence↗
PubMed2006

Balancing solvation and intramolecular interactions: toward a consistent generalized Born force field.

The efficient and accurate characterization of solvent effects is a key element in the theoretical and computational study of biological problems. Implicit solvent models, particularly generalized Born (GB) continuum electrostatics, have emerged as an attractive tool to study the structure and dynamics of biomolecules in various environments. Despite recent advances in this methodology, there remain limitations in the parametrization of many of these models. In the present work, we demonstrate that it is possible to achieve a balanced implicit solvent force field by further optimizing the input atomic radii in combination with adjusting the protein backbone torsional energetics. This parameter optimization is guided by the potentials of mean force (PMFs) between amino acid polar groups, calculated from explicit solvent free energy simulations, and by conformational equilibria of short peptides, obtained from extensive folding and unfolding replica exchange molecular dynamics (REX-MD) simulations. Through the application of this protocol, the delicate balance between the competing solvation forces and intramolecular forces appears to be better captured, and correct conformational equilibria for a range of both helical and beta-hairpin peptides are obtained. The same optimized force field also successfully folds both beta-hairpin trpzip2 and mini-protein Trp-Cage, indicating that it is quite robust. Such a balanced, physics-based force field will be highly applicable to a range of biological problems including protein folding and protein structural dynamics.

Amino Acid Sequence↗
PubMed2005

Structure of a cyclodextrin functionalized anionic clay: XRD analysis, spectroscopy, and computer simulations.

Carboxy-methyl beta-cyclodextrin (CMCD) cavities have been intercalated within the galleries of anionic clay, Mg-Al layered double hydroxide (LDH). The cyclodextrin functionalized LDH has been reported to adsorb neutral and nonpolar guest molecules. X-ray diffraction, IR, and Raman vibrational spectroscopy and (13)C CPMAS NMR have been used to characterize the confined CMCD molecules, whereas molecular dynamics simulations have been used to probe the interlayer arrangement and orientation of the intercalated species. Spectroscopic measurements as well as MD simulations show that there is no significant change in the geometry of the CMCD cavity on intercalation. Within the galleries of the anionic clay, the CMCD anions are arranged as bilayers with the carboxy methyl substituents, located at the narrower opening of the bucket-like cyclodextrin toroid, anchored to the LDH sheet. This arrangement leaves the wider opening of the CMCD anion facing away from the layers allowing the interior of the cyclodextrin cavity to be accessible to guest molecules. Finally, the hydrophobicity of the anchored cyclodextrin cavity has been characterized using fluorescence from pyrene included within it.

Anions↗
PubMed2005

Effect of hexafluoroisopropanol alcohol on the structure of melittin: a molecular dynamics simulation study.

The molecular mechanism by which HFIP stabilizes the alpha-helical structure of peptides is not well understood. In the present study, we use melittin as a model to gain insight into the details of the atomistic interactions of HFIP with the peptide. We have performed extensive comparative molecular dynamics simulations (up to 100 nsec) in the absence and in the presence of HFIP. In agreement with recent NMR experiments, the simulations show rapid loss of tertiary structure in water at pH 2 but much higher helicity in 35% HFIP. The MD simulations also indicate that melittin adopts a highly dynamic global structure in 35% HFIP solution with two alpha-helical segments sampling a wide range of angular orientations. The analysis of the HFIP distribution shows the tendency of HFIP to aggregate around the peptide, increasing the local cosolvent concentration to more than two times that in the bulk concentration. The correlation of local peptide structure with HFIP coating suggests that displacement of water at the peptide surface is the main contribution of HFIP in stabilizing the secondary structure of melittin. Finally, a stabilizing effect promoted by the presence of counter-ions was also observed in the simulations.

Computer Simulation↗
PubMed2006

Understanding macroscopic diffusion of adsorbed molecules in crystalline nanoporous materials via atomistic simulations.

The diffusion rates of molecules inside nanoporous materials lie at the heart of many large-scale industrial applications of these materials. Quantitatively describing this diffusion, particularly diffusion of chemical mixtures in situations leading to net mass transport, remains challenging. Molecular dynamics (MD) simulations can play an important complementary role to experiments in this area. This Account describes applications of MD to diffusion in nanoporous materials with a particular focus on macroscopic diffusion, that is, diffusion involving mass transport. These methods have made useful contributions to developing mixing theories for predicting multicomponent diffusion from single-component data and to screening new classes of materials for practical applications.

Journal Article↗
PubMed2005

Cyan fluorescent protein: molecular dynamics, simulations, and electronic absorption spectrum.

The dynamics and electronic absorption spectrum of enhanced cyan fluorescent protein (ECFP), a mutant of green fluorescent protein (GFP), have been studied by means of a 1 ns molecular dynamics (MD) simulation. The two X-ray conformations A' and B' of ECFP were considered. The chromophore was assumed to be neutral, and all titratable residues were taken in their standard protonation state at neutral pH. The protein was embedded in a box of water molecules (and counterions). The first result is that the two conformations A' and B' are found to be stable all along the simulation. Then, an analysis of the hydrogen-bond networks shows strong differences between the two conformations in the surroundings of the nitrogen atom of the indolic part of the chromophore. This is partly due to the imperfection in the beta barrel near the His148 residue, which allows the access of one solvent molecule inside the protein in conformation A'. Finally, quantum mechanical calculations of the electronic transition energies of the chromophore in the charge cloud of the protein and solvent water molecules were performed using the TDDFT method on 160 snapshots extracted every 5 ps of the MD trajectories. It is found that conformations A' and B' exhibit very similar spectra despite different H-bond networks involving the chromophore. This similarity is related to the weak charge transfer involved in the electronic transition and the weak electrostatic field created by ECFP near the chromophore, within the hypotheses made in the present simulation.

Computer Simulation↗
PubMed2004

Comparison of properties of Aib-rich peptides in crystal and solution: a molecular dynamics study.

In order to study the differences of the structural properties of Aib-rich peptides in solution and in the crystalline state, molecular dynamics (MD) simulations of the Aib-containing peptide II (pBrBz-(Aib)5-Leu-(Aib)2-OMe) were performed in the crystalline state, starting from two different conformers obtained experimentally by X-ray diffraction. The structural properties as derived from X-ray crystallography (e.g., torsional angles and hydrogen bonds) are well-reproduced in both constant-volume and constant-pressure simulations, although the force-field parameters used result in a too-high density of the crystals. Through comparison with the results from previous MD and nuclear magnetic resonance (NMR) studies of the very similar peptide I (Z-(Aib)s-Leu-(Aib)2-OMe) in dimethylsulfoxide (DMSO) solution, it is found that, in the crystal simulation, the conformational distribution of peptide II is much narrower than that in the solution simulation of peptide. I. This leads to a significant difference in 3 [symbol: see text] (HN, HC alpha) coupling constant values, in agreement with experimental data, whereas the NOE intensities or proton-proton distance bounds appear insensitive to the difference in conformational distribution. For small peptides the differences between their conformational distribution in the crystalline form and in solution may be much larger than for proteins, a fact which should be kept in mind when interpreting molecular properties in the solution state by using X-ray crystallographic data.

Bromides↗
PubMed1990

Conformational and helicoidal analysis of the molecular dynamics of proteins: "curves," dials and windows for a 50 psec dynamic trajectory of BPTI.

A new procedure for the graphic analysis of molecular dynamics (MD) simulations on proteins is introduced, in which comprehensive visualization of results and pattern recognition is greatly facilitated. The method involves determining the conformational and helicoidal parameters for each structure entering the analysis via the method "Curves," developed for proteins by Sklenar, Etchebest, and Lavery (Proteins: Structure, Function Genet. 6:46-60, 1989) followed by a novel computer graphic display of the results. The graphic display is organized systematically using conformation wheels ("dials") for each torsional parameter and "windows" on the range values assumed by the linear and angular helicoidal parameters, and is present in a form isomorphous with the primary structure per se. The complete time evolution of dynamic structure can then be depicted in a set of four composite figures. Dynamic aspects of secondary and tertiary structure are also provided. The procedure is illustrated with an analysis of a 50 psec in vacuo simulation on the 58 residue protein, bovine pancreatic trypsin inhibitor (BPTI), in the vicinity of the local minimum on the energy surface corresponding to a high resolution crystal structure. The time evolution of 272 conformational and 788 helicoidal parameters for BPTI is analyzed. A number of interesting features can be discerned in the analysis, including the dynamic range of conformational and helicoidal motions, the dynamic extent of 2 degrees structure motifs, and the calculated fluctuations in the helix axis. This approach is expected to be useful for a critical analysis of the effects of various assumptions about force field parameters, truncation of potentials, solvation, and electrostatic effects, and can thus contribute to the development of more reliable simulation protocols for proteins. Extensions of the analysis to present differential changes in conformational and helicoidal parameters is expected to be valuable in MD studies of protein complexes with substrates, inhibitors, and effectors and in determining the nature of structural changes in protein-protein interactions.

Biological Evolution↗
PubMed2005

Molecular dynamics studies of the full-length integrase-DNA complex.

We have carried out a molecular dynamics (MD) simulation of full-length HIV-1 integrase (IN) dimer complexed with viral DNA with the aim of gaining information about the enzyme motion and investigating the movement of the catalytic flexible loop (residues 140-149) thought to be essential in the catalytic mechanism of IN. During the simulation, we observed quite a different behavior of this region in the presence or absence of the viral DNA. In particular, the MD results underline the crucial role of the residue Tyr143 in the mechanism of integration of viral DNA into the host chromosome. The present findings confirm the experimental data (e.g., site-directed mutagenesis experiments) showing that the loop is involved in the integration reactions and its mobility is correlated with the catalytic activity of HIV-1 integrase.

Catalytic Domain↗
PubMed2004

Chloroethylnitrosourea-derived ethano cytosine and adenine adducts are substrates for Escherichia coli glycosylases excising analogous etheno adducts.

Exocyclic ethano DNA adducts are saturated etheno ring derivatives formed mainly by therapeutic chloroethylnitrosoureas (CNUs), which are also mutagenic and carcinogenic. In this work, we report that two of the ethano adducts, 3,N4-ethanocytosine (EC) and 1,N6-ethanoadenine (EA), are novel substrates for the Escherichia coli mismatch-specific uracil-DNA glycosylase (Mug) and 3-methyladenine DNA glycosylase II (AlkA), respectively. It has been shown previously that Mug excises 3,N4-ethenocytosine (epsilonC) and AlkA releases 1,N6-ethenoadenine (epsilonA). Using synthetic oligonucleotides containing a single ethano or etheno adduct, we found that both glycosylases had a approximately 20-fold lower excision activity toward EC or EA than that toward their structurally analogous epsilonC or epsilonA adduct. Both enzymes were capable of excising the ethano base paired with any of the four natural bases, but with varying efficiencies. The Mug activity toward EC could be stimulated by E. coli endonuclease IV and, more efficiently, by exonuclease III. Molecular dynamics (MD) simulations showed similar structural features of the etheno and ethano derivatives when present in DNA duplexes. However, also as shown by MD, the stacking interaction between the EC base and Phe 30 in the Mug active site is reduced as compared to the epsilonC base, which could account for the lower EC activity observed in this study.

Adenine↗
PubMed1998

Studies on the conformational behaviour of GlcNAc-Man3-GlcNAc2 oligosaccharides using molecular dynamics simulations.

Three-dimensional structures of the natural substrate unit for the enzyme N-acetylglucosamine-transferase II, GIcNAc-Man3-GlcNAc2, were investigated by molecular modelling methods. Molecular dynamics (MD) and molecular mechanics calculations on two hexasaccharides, namely GlcNAc-Man3-GlcNAc2-Asn and GlcNAc-Man3-GlcNAc2-OMe were performed by the Biosym/MSI software using the CVFF and CFF95 force fields in vacuum. The MD simulations were calculated for 3 ns at different simulation temperatures and for two values of dielectric constant, epsilon=1 and epsilon=4. From each 3 ns trajectory, 3050 structures have been optimized. The local minima obtained have been clustered into families exhibiting similar values of glycosidic torsional angles phi, psi, and omega. The influence of the simulation conditions and force fields used on the conformational behaviour and structure of the title oligosaccharides is discussed.

Acetylglucosamine↗
PubMed2003

Computational and experimental studies of the catalytic mechanism of Thermobifida fusca cellulase Cel6A (E2).

Mutagenesis experiments suggest that Asp79 in cellulase Cel6A (E2) from Thermobifida fusca has a catalytic role, in spite of the fact that this residue is more than 13 A from the scissile bond in models of the enzyme-substrate complex built upon the crystal structure of the protein. This suggests that there is a substantial conformational shift in the protein upon substrate binding. Molecular mechanics simulations were used to investigate possible alternate conformations of the protein bound to a tetrasaccharide substrate, primarily involving shifts of the loop containing Asp79, and to model the role of water in the active site complex for both the native conformation and alternative low-energy conformations. Several alternative conformations of reasonable energy have been identified, including one in which the overall energy of the enzyme-substrate complex in solution is lower than that of the conformation in the crystal structure. This conformation was found to be stable in molecular dynamics simulations with a cellotetraose substrate and water. In simulations of the substrate complexed with the native protein conformation, the sugar ring in the -1 binding site was observed to make a spontaneous transition from the (4)C(1) conformation to a twist-boat conformer, consistent with generally accepted glycosidase mechanisms. Also, from these simulations Tyr73 and Arg78 were found to have important roles in the active site. Based on the results of these various MD simulations, a new catalytic mechanism is proposed. Using this mechanism, predictions about the effects of changes in Arg78 were made which were confirmed by site-directed mutagenesis.

Actinomycetales↗
PubMed1995

Conformational modeling of elastin tetrapeptide Boc-Gly-Leu-Gly-Gly-NMe by molecular dynamics simulations with improvements to the thermalization procedure.

Molecular Dynamics Simulations (MD) at Constant-Temperature or Constant-total Energy for the conformational Global-Minimum (GM) of elastin tetrapeptide Boc-Gly-Leu-Gly-Gly-NMe have been performed. The thermalization problem concerning the initial state of Constant-Temperature MD has been solved developing two effective strategies. In the first one, the run starts from the room-temperature state reached by Molecular Dynamics Simulated Annealing (SA). In the second one, one starts from the annealed-state at low-temperature and performs a long constant-low-temperature run until the initial conformer is perfectly equilibrated. Then, the low-temperature equilibrated-state is used as initial state for MD at room-temperature. heuristic criteria on order to define the onset of steady-state have been established monitoring the hystories of collective parameters (e.g., the total energy, temperature, end-to-end distance, etc.) and their amplitude fluctuations. Moreover, the equilibrium between the system and the heat bath is verified analyzing the total linear momentum conservation by the time evolution of center mass velocity. The slow drift of total energy during Constant-total Energy MD has been corrected using a loose coupling between the system and the heat bath. Moreover, we have verified that the roto-translational motions do not affect significantly the properties of molecular vibrations. The librations of peptide unit inside the type II beta-turn [Gly1]C = 0 ... HN[Gly4], previously detected, were confirmed. Large -Gly-Gly- chain motions were identified and modeled as fluctuations occurring between the tetratepeptide GM and the saddle-point corresponding to the transition state of the conversion toward the extended-chain conformation. All these peptide motions could contribute to the elasticity mechanism of elastin.

Amino Acid Sequence↗
PubMed1996

Observation of the A-DNA to B-DNA transition during unrestrained molecular dynamics in aqueous solution.

A large challenge in molecular dynamics (MD) simulations of proteins and nucleic acids is to find the correct "experimental" geometry when a simulation is started a significant distance away from it. In this study, we have carried out four unrestrained approximately 1 ns length MD trajectories in aqueous solution on the DNA duplex d(CCAACGTTGG)2, two beginning in a canonical A-DNA structure and two beginning in a canonical B-DNA structure. As judged by root-mean-squared coordinate deviations, average structures computed from all four of the trajectories converge to within approximately 0.8 to 1.6 angstroms (all atoms) of each other, which is 1.3 to 1.7 angstroms (all atoms of the central six residues from each strand) and 3.1 to 3.6 angstroms (all atoms) away from the B-DNA-like X-ray structure reported for this sequence. To our knowledge, this is the first example of multiple nanosecond molecular dynamics trajectories with full representation of DNA charges, solvent and long range electrostatics that demonstrate both internal consistency (two different starting structures and four different trajectories lead to a consistent average structure) and considerable agreement with the X-ray crystal structure of this sequence and NMR data on duplex DNA in aqueous solution. This internal consistency of structure for a given sequence suggests that one can now begin to realistically examine sequence-dependent structural effects in DNA duplexes using molecular dynamics.

Base Sequence↗
PubMed2001

Molecular dynamics simulations of the insertion of two ideally amphipathic lytic peptides LK(15) and LK(9) in a 1,2-dimyristoylphosphatidylcholine monolayer.

We present here the results of 1-ns molecular dynamics (MD) simulations of two ideally amphipathic lytic peptides, namely LK(15) and LK(9), in a 1,2-dimyristoylphosphatidylcholine monolayer with two different cross-sectional areas per lipid of 80 A(2) (loose film) and 63 A(2) (tight standard film). These peptides are lytic, ideally amphipathic with a minimalist composition L(i)K(j) and the following sequences: H(2)N-KLLKLLLKLLLKLLK-CO-Ph for LK(15) and H(2)N-KLKLKLKLK-CO-Ph for LK(9). From experimental data, LK(15) exhibits an alpha-helical secondary structure, whereas LK(9) was found to form antiparallel beta-sheets at the interface of a DMPC monolayer. Whatever the specific lipid surface is, the two peptides exhibit very different behavior: the alpha-helix inserts deeply into the monolayer whereas the beta-sheeted peptide stays at the surface within the upper polar part of the monolayer. In all cases, a loose monolayer (80 A(2)) results in noticeable artifacts whereas a monolayer with standard specific surface leads to very reliable behavior well in accordance with experimental data. Despite their different insertion depth, the two peptides exhibit identical lytic efficiency. This is very likely a direct consequence of the very strong Van der Waals interactions between the fatty alkyl chains of the lipids and the highly lipophilic lower part of the peptide, resulting in an identical thinning of the two monolayers.

Amino Acid Sequence↗
PubMed2002

Theoretical calculation of the coiled-coil stability in water in the context of its possible use as a molecular rack.

The coiled-coil stability and rigidity may be of importance for molecular electronics (electronically bistable molecules). The coiled-coil binding free energy has been calculated using molecular dynamics (MD). The energy has been computed as a difference of the appropriate free energies; derived for the coiled-coil and isolated alpha-helices separately. All MD simulations have been performed using an explicit model of the solvent, whereas the continuum solvent approach has been applied to analyze the MD trajectories. The computed stability of the coiled-coil is of the order of -87 kcal/mol, i.e., about -1.2 kcal/mol per amino acid residue, and arises mainly from the electrostatic interactions and hydrophobic effect. The entropy term has been roughly estimated to be of the order of -22 kcal/mol. This assures that coiled-coil polypeptides may be used as a stable molecular scaffolding.

Amino Acid Sequence↗
PubMed2005

Molecular dynamics simulations of the periplasmic ferric-hydroxamate binding protein FhuD.

FhuD is a periplasmic binding protein (PBP) that, under iron-limiting conditions, transports various hydroxamate-type siderophores from the outer membrane receptor (FhuA) to the inner membrane ATP-binding cassette transporter (FhuBC). Unlike many other PBPs, FhuD possesses two independently folded domains that are connected by an alpha-helix rather than two or three central beta-strands. Crystal structures of FhuD with and without bound gallichrome have provided some insight into the mechanism of siderophore binding as well as suggested a potential mechanism for FhuD binding to FhuB. Since the alpha-helix connecting the two domains imposes greater rigidity on the structure relative to the beta-strands in other 'classical' PBPs, these structures reveal no large conformational change upon binding a hydroxamate-type siderophore. Therefore, it is difficult to explain how the inner membrane transporter FhuB can distinguish between ferrichrome-bound and ferrichrome-free FhuD. In the current study, we have employed a 30 ns molecular dynamics simulation of FhuD with its bound siderophore removed to explore the dynamic behavior of FhuD in the substrate-free state. The MD simulation suggests that FhuD is somewhat dynamic with a C-terminal domain closure of 6 degrees upon release of its siderophore. This relatively large motion suggests differences that would allow FhuB to distinguish between ferrichrome-bound and ferrichrome-free FhuD.

Amino Acid Sequence↗