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 127 records · Page 7Linked to original sources
PubMed1995

Structural fluctuation of methyl N,N'-diacetyl-beta-D-chitobioside in vacuo and in aqueous solution: molecular dynamics simulations and proton NMR spectroscopy.

Molecular dynamics (MD) simulations of methyl N,N'-diacetyl-beta-D-chitobioside (GlcNAc beta-(1-->4)GlcNAc beta-OMe) have been performed both in vacuo and in aqueous solution with the explicit inclusion of the solvent water molecules. The beta-(1-->4) glycosidic linkage fluctuates considerably, over a range of +/- 10 degrees, in each of the MD simulations in vacuo and in aqueous solution. The intra- and inter-residue hydrogen bonds in vacuo are replaced by intermolecular hydrogen bonds with the solvent water molecules in aqueous solution. Multiple conformations (gg and gt) exist for the exocyclic hydroxymethyl groups. The results of the MD simulations are compared with those of 1H-1H nuclear Overhauser effect measurements.

Carbohydrate Conformation↗
PubMed2006

Coarse-grained molecular dynamics simulations of membrane proteins and peptides.

Molecular dynamics (MD) simulations provide a valuable approach to the dynamics, structure, and stability of membrane-protein systems. Coarse-grained (CG) models, in which small groups of atoms are treated as single particles, enable extended (>100 ns) timescales to be addressed. In this study, we explore how CG-MD methods that have been developed for detergents and lipids may be extended to membrane proteins. In particular, CG-MD simulations of a number of membrane peptides and proteins are used to characterize their interactions with lipid bilayers. CG-MD is used to simulate the insertion of synthetic model membrane peptides (WALPs and LS3) into a lipid (PC) bilayer. WALP peptides insert in a transmembrane orientation, whilst the LS3 peptide adopts an interfacial location, both in agreement with experimental biophysical data. This approach is extended to a transmembrane fragment of the Vpu protein from HIV-1, and to the coat protein from fd phage. Again, simulated protein/membrane interactions are in good agreement with solid state NMR data for these proteins. CG-MD has also been applied to an M3-M4 fragment from the CFTR protein. Simulations of CFTR M3-M4 in a detergent micelle reveal formation of an alpha-helical hairpin, consistent with a variety of biophysical data. In an I231D mutant, the M3-M4 hairpin is additionally stabilized via an inter-helix Q207/D231 interaction. Finally, CG-MD simulations are extended to a more complex membrane protein, the bacterial sugar transporter LacY. Comparison of a 200 ns CG-MD simulation of LacY in a DPPC bilayer with a 50 ns atomistic simulation of the same protein in a DMPC bilayer shows that the two methods yield comparable predictions of lipid-protein interactions. Taken together, these results demonstrate the utility of CG-MD simulations for studies of membrane/protein interactions.

Amino Acid Sequence↗
PubMed2006

Pathways of H2 toward the active site of [NiFe]-hydrogenase.

Hydrogenases catalyze the reversible oxidation of molecular hydrogen (H(2)), but little is known about the diffusion of H(2) toward the active site. Here we analyze pathways for H(2) permeation using molecular dynamics (MD) simulations in explicit solvent. Various MD simulation replicates were done, to improve the sampling of the system states. H(2) easily permeates hydrogenase in every simulation and it moves preferentially in channels. All H(2) molecules that reach the active site made their approach from the side of the Ni ion. H(2) is able to reach distances of <4 A from the active site, although after 6 A permeation is difficult. In this region we mutated Val-67 into alanine and perform new MD simulations. These simulations show an increase of H(2) inside the protein and at lower distances from the active site. This valine can be a control point in the H(2) access to the active center.

Binding Sites↗
PubMed1998

Molecular dynamic simulations of environment and sequence dependent DNA conformations: the development of the BMS nucleic acid force field and comparison with experimental results.

Molecular dynamic (MD) simulations using the BMS nucleic acid force field produce environment and sequence dependent DNA conformations that closely mimic experimentally derived structures. The parameters were initially developed to reproduce the potential energy surface, as defined by quantum mechanics, for a set of small molecules that can be used as the building blocks for nucleic acid macromolecules (dimethyl phosphate, cyclopentane, tetrahydrofuran, etc.). Then the dihedral parameters were fine tuned using a series of condensed phase MD simulations of DNA and RNA (in zero added salt, 4M NaCl, and 75% ethanol solutions). In the tuning process the free energy surface for each dihedral was derived from the MD ensemble and fitted to the conformational distributions and populations observed in 87 A- and B-DNA x-ray and 17 B-DNA NMR structures. Over 41 nanoseconds of MD simulations are presented which demonstrate that the force field is capable of producing stable trajectories, in the correct environments, of A-DNA, double stranded A-form RNA, B-DNA, Z-DNA, and a netropsin-DNA complex that closely reproduce the experimentally determined and/or canonical DNA conformations. Frequently the MD averaged structure is closer to the experimentally determined structure than to the canonical DNA conformation. MD simulations of A- to B- and B- to A-DNA transitions are also shown. A-DNA simulations in a low salt environment cleanly convert into the B-DNA conformation and converge into the RMS space sampled by a low salt simulation of the same sequence starting from B-DNA. In MD simulations using the BMS force field the B-form of d(GGGCCC)2 in a 75% ethanol solution converts into the A-form. Using the same methodology, parameters, and conditions the A-form of d(AAATTT)2 correctly converts into the B-DNA conformation. These studies demonstrate that the force field is capable of reproducing both environment and sequence dependent DNA structures. The 41 nanoseconds (nsec) of MD simulations presented in this paper paint a global picture which suggests that the DNA structures observed in low salt solutions are largely due to the favorable internal energy brought about by the nearly uniform screening of the DNA electrostatics. While the conformations sampled in high salt or mixed solvent environments occur from selective and asymmetric screening of the phosphate groups and DNA grooves, respectively, brought about by sequence induced ion and solvent packing.

Adenine↗
PubMed2005

Comparison of the glass transition temperature and fragility parameter of isomalto-olygomer predicted by molecular dynamics simulations with those measured by differential scanning calorimetry.

The purpose of this study is to examine whether molecular dynamics (MD) simulations using a commercially available software for personal computers can estimate the glass transition temperature (Tg) of amorphous systems containing pharmaceutically-relevant excipients. MD simulations were carried out with an amorphous matrix model constructed from isomaltoheptaose, and the Tg estimated from the calculated density versus temperature profile was compared with the Tg measured by differential scanning calorimetry (DSC) for freeze-dried isomalto-oligomer having an average molecular weight close to that of isomaltoheptaose. The Tg values determined by DSC were lower by 10 to 20 K than those extrapolated from the Tg values estimated by MD simulation. Fragility parameter was estimated to be 56 and 51 from MD simulation and from DSC measurement, respectively. Thus, the results suggest that MD simulation can provide approximate estimates for the Tg and fragility parameter of amorphous formulations. However, a reduction of the cooling rate, achievable by sufficiently elongating the simulation duration, is necessary for more accurate estimation.

Algorithms↗
PubMed2001

Folding study of an Aib-rich peptide in DMSO by molecular dynamics simulations.

To evaluate the ability of molecular dynamics (MD) simulations using atomic force-fields to correctly predict stable folded conformations of a peptide in solution, we show results from MD simulations of the reversible folding of an octapeptide rich in alpha-aminoisobutyric acid (2-amino-2-methyl-propanoic acid, Aib) solvated in di-methyl-sulfoxide (DMSO). This solvent generally prevents the formation of secondary structure, whereas Aib-rich peptides show a high propensity to form secondary structural elements, in particular 3(10)- and alpha-helical structures. Aib is, moreover, achiral, so that Aib-rich peptides can form left- or right-handed helices depending on the overall composition of the peptide, the temperature, and the solvation conditions. This makes the system an interesting case to study the ensembles of peptide conformations as a function of temperature by MD simulation. Simulations involving the folding and unfolding of the peptide were performed starting from two initial structures, a right-handed alpha-helical structure and an extended structure, at three temperatures, 298 K, 340 K, and 380 K, and the results are compared with experimental nuclear magnetic resonance (NMR) data measured at 298 K and 340 K. The simulations generally reproduce the available experimental nuclear Overhauser effect (NOE) data, even when a wide range of conformations is sampled at each temperature. The importance of adequate statistical sampling in order to reliably interpret the experimental data is discussed.

Dimethyl Sulfoxide↗
PubMed2006

Phenol-benzene complexation dynamics: quantum chemistry calculation, molecular dynamics simulations, and two dimensional IR spectroscopy.

Molecular dynamics (MD) simulations and quantum mechanical electronic structure calculations are used to investigate the nature and dynamics of the phenol-benzene complex in the mixed solvent, benzene/CCl4. Under thermal equilibrium conditions, the complexes are continuously dissociating and forming. The MD simulations are used to calculate the experimental observables related to the phenol hydroxyl stretching mode, i.e., the two dimensional infrared vibrational echo spectrum as a function of time, which directly displays the formation and dissociation of the complex through the growth of off-diagonal peaks, and the linear absorption spectrum, which displays two hydroxyl stretch peaks, one for the complex and one for the free phenol. The results of the simulations are compared to previously reported experimental data and are found to be in quite reasonable agreement. The electronic structure calculations show that the complex is T shaped. The classical potential used for the phenol-benzene interaction in the MD simulations is in good accord with the highest level of the electronic structure calculations. A variety of other features is extracted from the simulations including the relationship between the structure and the projection of the electric field on the hydroxyl group. The fluctuating electric field is used to determine the hydroxyl stretch frequency-frequency correlation function (FFCF). The simulations are also used to examine the number distribution of benzene and CCl4 molecules in the first solvent shell around the phenol. It is found that the distribution is not that of the solvent mole fraction of benzene. There are substantial probabilities of finding a phenol in either a pure benzene environment or a pure CCl4 environment. A conjecture is made that relates the FFCF to the local number of benzene molecules in phenol's first solvent shell.

Benzene↗
PubMed2006

Local order in aqueous NaCl solutions and pure water: X-ray scattering and molecular dynamics simulations study.

The microstructures of pure water and aqueous NaCl solutions over a wide range of salt concentrations (0-4 m) under ambient conditions are characterized by X-ray scattering and molecular dynamics (MD) simulations. MD simulations are performed with the rigid SPC water model as a solvent, while the ions are treated as charged Lennard-Jones particles. Simulated data show that the first peaks in the O...O and O...H pair correlation functions clearly decrease in height with increasing salt concentration. Simultaneously, the location of the second O...O peak, the signature of the so-called tetrahedral structure of water, gradually disappears. Consequently, the degree of hydrogen bonding in liquid water decreases when compared to pure fluid. MD results also show that the hydration number around the cation decreases as the salt concentration increases, which is most likely because some water molecules in the first hydration shell are occasionally substituted by chlorine. In addition, the fraction of contact ion pairs increases and that of solvent-separated ion pairs decreases. Experimental data are analyzed to deduce the structure factors and the pair correlation functions of each system. X-ray results clearly show a perturbation of the association structure of the solvent and highlight the appearance of new interactions between ions and water. A model of intermolecular arrangement via MD results is then proposed to describe the local order in each system, as deduced from X-ray scattering data.

Computer Simulation↗
PubMed1996

Conformational sampling by NMR solution structures calculated with the program DIANA evaluated by comparison with long-time molecular dynamics calculations in explicit water.

The NMR solution structure of bovine pancreatic trypsin inhibitor (BPTI) obtained by distance geometry calculations with the program DIANA is compared with groups of conformers generated by molecular dynamics (MD) simulations in explicit water at ambient temperature and pressure. The MD simulations started from a single conformer and were free or restrained either by the experimental NOE distance restraints or by time-averaged restraints; the groups of conformers were collected either in 10 ps intervals during 200 ps periods of simulation, or in 50 ps intervals during a 1 ns period of simulation. Overall, these comparisons show that the standard protein structure determination protocol with the program DIANA provides a picture of the protein structure that is in agreement with MD simulations using "realistic" potential functions over a nanosecond timescale. For well-constrained molecular regions there is a trend in the free MD simulation of duration 1 ns that the sampling of the conformation space is slightly increased relative to the DIANA calculations. In contrast, for surface-exposed side-chains that are less extensively constrained by the NMR data, the DIANA conformers tend to sample larger regions of conformational space than conformers selected from any of the MD trajectories. Additional insights into the behavior of surface side-chains come from comparison of the MD runs of 200 ps or 1 ns duration. In this time range the sampling of conformation space by the protein surface depends strongly on the length of the simulation, which indicates that significant side-chain transitions occur on the nanosecond timescale and that much longer simulations will be needed to obtain statistically significant data on side-chain dynamics.

Amino Acid Sequence↗
PubMed2004

Theoretical investigations on Azotobacter vinelandii ferredoxin I: effects of electron transfer on protein dynamics.

Structural, energetic, and dynamical studies of Azotobacter vinelandii ferredoxin I are presented for native and mutant forms. The protein contains two iron-sulfur clusters, one of which ([3Fe-4S]) is believed to play a central role in the electron-coupled proton transfer. Different charge sets for the [3Fe-4S] cluster in its reduced and oxidized state are calculated with broken symmetry ab initio density functional theory methods and used in molecular dynamics (MD) simulations. The validity of the ab initio calculations is assessed by comparing partially optimized structures of the [3Fe-4S] clusters with x-ray structures. Possible proton transfer pathways between the protein and the iron-sulfur cluster are examined by both MD simulations and ab initio calculations. The MD simulations identify three main-chain hydrogen atoms--HN(13), HN(14), and HN(16)--that are within H-bonding distance of the [3Fe-4S] cluster throughout the MD simulations. They could thus play a role in the proton transfer from the protein to the iron-sulfur cluster. By contrast, the HD2(15) atom of the Asp-15 is seldom close enough to the [3Fe-4S] cluster to transfer a proton. Poisson-Boltzmann calculations indicate that there is a low, but nonzero probability, that Asp-15 is protonated at pH 7; this is a requirement for it to serve as a proton donor. Ab initio calculations with a fragment model for the protein find similar behavior for the transfer of a proton from the OH of the protonated side chain and the main-chain NH of Asp-15. The existence of a stable salt bridge between Asp-15 and Lys-84 in the D15E mutant, versus its absence in the wild-type, has been suggested as the cause of the difference in the rate of proton transfer. Extensive MD simulations were done to test this idea; the results do not support the proposal. The present findings, together with the available data, serve as the basis for an alternative proposal for the mechanism of the coupled electron-proton transfer reaction in ferredoxin I.

Azotobacter vinelandii↗
PubMedDate not supplied

[Generalized ordering parameter S2 for N-H peptide bonding as measure of the conformational flexibility of proteins: comparison of algorithms of S2 calculation from molecular dynamics simulation data].

NMR spectroscopy in combination with molecular dynamics (MD) simulation technique is a powerful tool for investigation of biomolecule conformations. In the present study we compare the two algorithms of calculation of the generalized order parameter S2 for the peptide N-H bond from MD simulation data. The program based on this algorithm has been written and the effect of an adjustible parameter on results of calculation was studied. HIV-1 protease MD simulation was performed during 7.36 ns using GROMACS 3.1.4. software package. The calculated order parameter values are in good agreement with those obtained by other authors using NMR.

Algorithms↗
PubMed2001

Successful molecular dynamics simulation of two zinc complexes bridged by a hydroxide in phosphotriesterase using the cationic dummy atom method.

I report herein two 2.0 ns (1.0 fs time step) MD simulations of two zinc complexes bridged by a hydroxide in phosphotriesterase (PTE) employing the nonbonded method and the cationic dummy atom method that uses virtual atoms to impose orientational requirement for zinc ligands. The cationic dummy atom method was able to simulate the four-ligand coordination of the two zinc complexes in PTE. The distance (3.39 +/- 0.07A) between two nearby zinc ions in the time-average structure of PTE derived from the MD simulation using the cationic dummy atoms matched that in the X-ray structure (3.31 +/- 0.001A). Unequivocally, the time-average structure of PTE was able to fit into the experimentally determined difference electron density map of the corresponding X-ray structure. The results demonstrate the practicality of the cationic dummy atom method for MD simulations of zinc proteins bound with multiple zinc ions. In contrast, a 2.0 ns (1.0 fs time step) MD simulation using the nonbonded method revealed a striking difference in the active site between the X-ray structure and the time-average structure that was unable to fit into the density map of PTE. The results suggest that caution should be used in the MD simulations using the nonbonded method.

Algorithms↗
PubMed2006

Blocking of the nicotinic acetylcholine receptor ion channel by chlorpromazine, a noncompetitive inhibitor: A molecular dynamics simulation study.

A large series of pharmacological agents, distinct from the typical competitive antagonists, block in a noncompetitive manner the permeability response of the nicotinic acetylcholine receptor (nAChR) to the neurotransmitter acetylcholine. Taking the neuroleptic chlorpromazine (CPZ) as an example of such agents, the blocking mechanism of noncompetitive inhibitors to the ion channel pore of the nAChR has been explored at the atomic level using both conventional and steered molecular dynamics (MD) simulations. Repeated steered MD simulations have permitted calculation of the free energy (approximately 36 kJ/mol) of CPZ binding and identification of the optimal site in the region of the serine and leucine rings, at approximately 4 A from the pore entrance. Coulomb and the Lennard-Jones interactions between CPZ and the ion channel as well as the conformational fluctuations of CPZ were examined to assess the contribution of each to the binding of CPZ to the nAChR. The MD simulations disclose a dynamic interaction of CPZ binding to the nAChR ionic channel. The cationic ammonium head of CPZ forms strong hydrogen bonds with Glu262 (alpha), Asp268 (beta), Glu272 (beta), Ser276 (beta), Glu280 (delta), Gln271 (gamma), Glu275 (gamma), and Asn279 (gamma) nAChR residues. Finally, the conventional MD simulation of CPZ at its identified binding site demonstrates that the binding of CPZ not only blocks ion transport through the channel but also markedly inhibits the conformational transitions of the channel, necessary for nAChR to carry out its biological function.

Acetylcholine↗
PubMed1998

Molecular dynamics simulations on parallel and antiparallel C.G*G triplexes.

Molecular dynamics (MD) studies have been carried out on the Hoogsteen hydrogen bonded parallel and the reverse Hoogsteen hydrogen bonded antiparallel C.G*G triplexes. Earlier, the molecular mechanics studies had shown that the parallel structure was energetically more favourable than the antiparallel structure. To characterize the structural stability of the two triplexes and to investigate whether the antiparallel structure can transit to an energetically more favourable structure, due to the local fluctuations in the structure during the MD simulation, the two structures were subjected to 200ps of constant temperature vacuum MD simulations at 300K. Initially no constraints were applied to the structures and it was observed that for the antiparallel triplex, the structure showed a large root mean square deviation from the starting structure within the first 12ps and the N4-H41--O6 hydrogen bond in the WC duplex got distorted due to a high propeller twist and a moderate increase in the opening angle in the basepairs. Starting from an initial value of 30 degrees , helical twist of the average structure from this simulation had a value of 36 degrees , while the parallel structure stabilized at a twist of 33 degrees. In spite of the hydrogen bond distortions in the antiparallel triplex, it was energetically comparable to the parallel triplex. To examine the structural characteristics of an undistorted structure, another MD simulation was performed on the antiparallel triplex by constraining all the hydrogen bonds. This structure stabilized at an average twist of 33 degrees. In the course of the dynamics though the energy of the molecule - compared to the initial structure - improved, it did not become comparable to the parallel structure. Energy minimization studies performed in the presence of explicit water and counterions also showed the two structures to be equally favourable energetically. Together these results indicate that the parallel C.G*G triplex with Hoogsteen hydrogen bonds also represents a stereochemically and energetically favourable structure for this class of triplexes.

Carbohydrates↗
PubMed1994

Molecular dynamics simulation of equine infectious anemia virus Tat protein in water and in 40% trifluoroethanol.

Two molecular dynamics (MD) simulations were performed in order to increase the understanding of the dependence of protein conformation on solvent environment. The protein used for these simulations is the transcriptional activator of the equine infectious anemia virus (EIAV-Tat). The structure of this protein has been determined by nuclear magnetic resonance (NMR) in aqueous solution (Willbold et al., Science 264, 1584 (1994)) and in 40% (v/v) trifluoroethanol (TFE) (Sticht et al., Eur. J. Biochem., submitted) showing considerable differences in the stability of the secondary structure elements. In order to investigate the influence of the solvent MD simulations (300 K: 200 ps) were carried out in water and in a solvent containing 40% (v/v) TFE. In both simulations the structure as determined in 40% TFE by NMR showing three-helices and a tight type II turn, was used as the initial structure. The MD simulations clearly indicate a decreased stability of the secondary structure elements in aqueous environment as made obvious by larger atomic motions and stronger fluctuations in the length of the hydrogen bonds. Complete unfolding of the helices was not observed on a 200 ps timescale. The root mean square deviation (RMSD) values of the backbone atoms after 200 ps simulation compared to the starting structure underline the strong influence of the solvent on the protein stability. This RMSD value is 1.95 A for the simulation in water and 1.29 A for the simulation in TFE/water. This result supports the notion that TFE acts as a secondary structure inducing and stabilizing solvent. The differences apparent from the MD simulations are in good agreement with the data derived from NMR measurements, showing the relevance of MD as a method for estimating conformational and dynamical properties of proteins.

Amino Acid Sequence↗
PubMed1996

Molecular dynamics simulations of adipocyte lipid-binding protein: effect of electrostatics and acyl chain unsaturation.

Molecular dynamics (MD) simulations have been performed on adipocyte lipid-binding protein, using the apo and holo forms, bound with stearic and oleic acid. The contribution of electrostatics to protein dynamics and ligand stabilization was assayed by perturbing the electrostatic charge of Arg106 and Arg126 (positive-->neutral) and the fatty acid (132H) headgroup (negative-->neutral). MD simulations for charged holo forms demonstrated significantly greater electrostatic binding energy and a more stabilized hydrogen bond network than simulations performed using neutral forms. Electrostatics, however, appeared to have little effect on fatty acid behavior, e.g., fluctuation of the dihedral head group; number of dihedral transitions within the acyl chain; and change in the end-to-end distance for fatty acid. Instead, fatty acid behavior appeared to be dictated by the presence or absence of an unsaturated bond within the acyl chain. A significantly greater number of transitions were observed during MD simulations in oleic than stearic acid. In addition, significantly greater fluctuation was observed for oleic acid, within the C2 headgroup and C9 and C11 dihedrals (which lie adjacent to the olefin bond of oleic acid). The dynamic behavior of the acyl chain may thereby be more a property of van der Waals contact, and the degree of acyl chain unsaturation, than a function of electrostatics. In the absence of fatty acid, an increase in distance between guanidino carbon centered atoms of Arg126 and Arg106 was observed during MD simulations of the charged apo form. This effect not observed with the neutral apo form or in any of the holo complexes and, presumbably, was a result of repulsion between the negatively charged arginine sidechains. Conserved waters reflected substantially lower mean-square displacement (msd) in all simulations, except the neutral apo form. This suggests that the presence of either charged amino acids or lipid provides increased order for water within the binding pocket. These results provide a dynamic perspective of the interactive nature within the FABP binding pocket regulated in a complex manner by the electrostatics within the binding cavity, acyl chain structure and behavior, and water energetics.

Carrier Proteins↗
PubMed2005

Electronic relaxation dynamics of Ni2+-ion aqueous solution: molecular-dynamics simulation.

Electronic relaxation dynamics of Ni2+-ion aqueous solution is investigated using molecular-dynamics (MD) simulations with the model-effective Hamiltonian developed previously. The nonadiabatic transition rates from the first three excited states to the ground state are evaluated by the golden rule formula with the adiabatic MD simulations. The MD simulations with the fewest-switch surface-hopping method are also carried out to obtain a more detailed description of the electronic relaxation dynamics among the excited states. We found out that the transitions among the three excited states are very fast, in the order of 10 fs, while the transition between the excited and ground states is slow, about 800 ps. These findings are consistent with the time scales of energy dissipation detected by the transient lens experiment. In both simulations, we explore the effects of the quantum decoherence, where the decoherence functions are derived by the energy-gap dynamics with the displaced harmonic-oscillator model.

Journal Article↗
PubMed2004

Molecular dynamics simulations of DNA curvature and flexibility: helix phasing and premelting.

Recent studies of DNA axis curvature and flexibility based on molecular dynamics (MD) simulations on DNA are reviewed. The MD simulations are on DNA sequences up to 25 base pairs in length, including explicit consideration of counterions and waters in the computational model. MD studies are described for ApA steps, A-tracts, for sequences of A-tracts with helix phasing. In MD modeling, ApA steps and A-tracts in aqueous solution are essentially straight, relatively rigid, and exhibit the characteristic features associated with the B'-form of DNA. The results of MD modeling of A-tract oligonucleotides are validated by close accord with corresponding crystal structure results and nuclear magnetic resonance (NMR) nuclear Overhauser effect (NOE) and residual dipolar coupling (RDC) structures of d(CGCGAATTCGCG) and d(GGCAAAAAACGG). MD simulation successfully accounts for enhanced axis curvature in a set of three sequences with phased A-tracts studied to date. The primary origin of the axis curvature in the MD model is found at those pyrimidine/purine YpR "flexible hinge points" in a high roll, open hinge conformational substate. In the MD model of axis curvature in a DNA sequence with both phased A-tracts and YpR steps, the A-tracts appear to act as positioning elements that make the helix phasing more precise, and key YpR steps in the open hinge state serve as curvature elements. Our simulations on a phased A-tract sequence as a function of temperature show that the MD simulations exhibit a premelting transition in close accord with experiment, and predict that the mechanism involves a B'-to-B transition within A-tracts coupled with the prediction of a transition in key YpR steps from the high roll, open hinge, to a low roll, closed hinge substate. Diverse experimental observations on DNA curvature phenomena are examined in light of the MD model with no serious discrepancies. The collected MD results provide independent support for the "non-A-tract model" of DNA curvature. The "junction model" is indicated to be a special case of the non-A-tract model when there is a Y base at the 5' end of an A-tract. In accord with crystallography, the "ApA wedge model" is not supported by MD.

Base Sequence↗