The 22-mer c-kit promoter sequence folds into a parallel-stranded quadruplex with a unique structure, which has been elucidated by crystallographic and NMR methods and shows a high degree of structural conservation. We have carried out a series of extended (up to 10 μs long, ∼50 μs in total) molecular dynamics simulations to explore conformational stability and loop dynamics of this quadruplex. Unfolding no-salt simulations are consistent with a multi-pathway model of quadruplex folding and identify the single-nucleotide propeller loops as the most fragile part of the quadruplex. Thus, formation of propeller loops represents a peculiar atomistic aspect of quadruplex folding. Unbiased simulations reveal μs-scale transitions in the loops, which emphasizes the need for extended simulations in studies of quadruplex loops. We identify ion binding in the loops which may contribute to quadruplex stability. The long lateral-propeller loop is internally very stable but extensively fluctuates as a rigid entity. It creates a size-adaptable cleft between the loop and the stem, which can facilitate ligand binding. The stability gain by forming the internal network of GA base pairs and stacks of this loop may be dictating which of the many possible quadruplex topologies is observed in the ground state by this promoter quadruplex.

• MeSH
• Nucleic Acid Denaturation MeSH
• Potassium chemistry   MeSH
• G-Quadruplexes * MeSH
• Cations MeSH
• Base Pairing MeSH
• Promoter Regions, Genetic * MeSH
• Proto-Oncogene Proteins c-kit genetics   MeSH
• Molecular Dynamics Simulation MeSH
• Sodium chemistry   MeSH
• Publication type
• Journal Article MeSH
• Research Support, Non-U.S. Gov't MeSH
• Names of Substances
• Potassium MeSH
• Cations MeSH
• Proto-Oncogene Proteins c-kit MeSH
• Sodium MeSH

Explicit solvent molecular dynamics simulations have been used to complement preceding experimental and computational studies of folding of guanine quadruplexes (G-DNA). We initiate early stages of unfolding of several G-DNAs by simulating them under no-salt conditions and then try to fold them back using standard excess salt simulations. There is a significant difference between G-DNAs with all-anti parallel stranded stems and those with stems containing mixtures of syn and anti guanosines. The most natural rearrangement for all-anti stems is a vertical mutual slippage of the strands. This leads to stems with reduced numbers of tetrads during unfolding and a reduction of strand slippage during refolding. The presence of syn nucleotides prevents mutual strand slippage; therefore, the antiparallel and hybrid quadruplexes initiate unfolding via separation of the individual strands. The simulations confirm the capability of G-DNA molecules to adopt numerous stable locally and globally misfolded structures. The key point for a proper individual folding attempt appears to be correct prior distribution of syn and anti nucleotides in all four G-strands. The results suggest that at the level of individual molecules, G-DNA folding is an extremely multi-pathway process that is slowed by numerous misfolding arrangements stabilized on highly variable timescales.

• MeSH
• DNA chemistry   MeSH
• G-Quadruplexes * MeSH
• DNA, Single-Stranded chemistry   MeSH
• Humans MeSH
• Molecular Dynamics Simulation * MeSH
• Telomere chemistry   MeSH
• Check Tag
• Humans MeSH
• Publication type
• Journal Article MeSH
• Research Support, Non-U.S. Gov't MeSH
• Names of Substances
• DNA MeSH
• DNA, Single-Stranded MeSH

The article reviews the application of biomolecular simulation methods to understand the structure, dynamics and interactions of nucleic acids with a focus on explicit solvent molecular dynamics simulations of guanine quadruplex (G-DNA and G-RNA) molecules. While primarily dealing with these exciting and highly relevant four-stranded systems, where recent and past simulations have provided several interesting results and novel insight into G-DNA structure, the review provides some general perspectives on the applicability of the simulation techniques to nucleic acids.

• MeSH
• DNA chemistry   MeSH
• G-Quadruplexes * MeSH
• Guanine chemistry   MeSH
• Nucleic Acid Conformation MeSH
• Ligands MeSH
• RNA chemistry   MeSH
• Solvents chemistry   MeSH
• Molecular Dynamics Simulation * MeSH
• Telomere chemistry   MeSH
• Hydrogen Bonding MeSH
• Publication type
• Journal Article MeSH
• Research Support, Non-U.S. Gov't MeSH
• Review MeSH
• Research Support, N.I.H., Extramural MeSH
• Research Support, U.S. Gov't, Non-P.H.S. MeSH
• Names of Substances
• DNA MeSH
• Guanine MeSH
• Ligands MeSH
• RNA MeSH
• Solvents MeSH

The glmS catalytic riboswitch is part of the 5'-untranslated region of mRNAs encoding glucosamine-6-phosphate (GlcN6P) synthetase (glmS) in numerous gram-positive bacteria. Binding of the cofactor GlcN6P induces site-specific self-cleavage of the RNA. However, the detailed reaction mechanism as well as the protonation state of the glmS reactive form still remains elusive. To probe the dominant protonation states of key active site residues, we carried out explicit solvent molecular dynamic simulations involving various protonation states of three crucial active site moieties observed in the available crystal structures: (i) guanine G40 (following the Thermoanaerobacter tengcongensis numbering), (ii) the GlcN6P amino/ammonium group, and (iii) the GlcN6P phosphate moiety. We found that a deprotonated G40(-) seems incompatible with the observed glmS active site architecture. Our data suggest that the canonical form of G40 plays a structural role by stabilizing an in-line attack conformation of the cleavage site A-1(2'-OH) nucleophile, rather than a more direct chemical role. In addition, we observe weakened cofactor binding upon protonation of the GlcN6P phosphate moiety, which explains the experimentally observed increase in K(m) with decreasing pH. Finally, we discuss a possible role of cofactor binding and its interaction with the G65 and G1 purines in structural stabilization of the A-1(2'-OH) in-line attack conformation. On the basis of the identified dominant protonation state of the reaction precursor, we propose a hypothesis of the self-cleavage mechanism in which A-1(2'-OH) is activated as a nucleophile by the G1(pro-R(p)) nonbridging oxygen of the scissile phosphate, whereas the ammonium group of GlcN6P acts as the general acid protonating the G1(O5') leaving group.

• MeSH
• Glucose-6-Phosphate analogs & derivatives   metabolism   MeSH
• Glucosamine analogs & derivatives   metabolism   MeSH
• Glutamine-Fructose-6-Phosphate Transaminase (Isomerizing) genetics   MeSH
• Catalytic Domain * MeSH
• Coenzymes metabolism   MeSH
• Nucleic Acid Conformation MeSH
• Molecular Sequence Data MeSH
• Protons * MeSH
• RNA, Catalytic chemistry   genetics   metabolism   MeSH
• Base Sequence MeSH
• Molecular Dynamics Simulation * MeSH
• Thermoanaerobacter enzymology   MeSH
• Publication type
• Journal Article MeSH
• Research Support, Non-U.S. Gov't MeSH
• Research Support, N.I.H., Extramural MeSH
• Names of Substances
• glucosamine 6-phosphate MeSH Browser
• Glucose-6-Phosphate MeSH
• Glucosamine MeSH
• Glutamine-Fructose-6-Phosphate Transaminase (Isomerizing) MeSH
• Coenzymes MeSH
• Protons * MeSH
• RNA, Catalytic MeSH

This review provides a critical assessment of the advantages and limitations of modeling methods available for guanine quadruplex (G-DNA) molecules. We characterize the relations of simulations to the experimental techniques and explain the actual meaning and significance of the results. The following aspects are discussed: pair-additive approximation of the empirical force fields, sampling limitations stemming from the simulation time and accuracy of description of base stacking, H-bonding, sugar-phosphate backbone and ions by force fields. Several methodological approaches complementing the classical explicit solvent molecular dynamics simulations are commented on, including enhanced sampling methods, continuum solvent methods, free energy calculations and gas phase simulations. The successes and pitfalls of recent simulation studies of G-DNA are demonstrated on selected results, including studies of cation interactions and dynamics of G-DNA stems, studies of base substitutions (inosine, thioguanine and mixed tetrads), analysis of possible kinetic intermediates in folding pathway of a G-DNA stem and analysis of loop regions of G-DNA molecules.

• MeSH
• DNA chemistry   MeSH
• G-Quadruplexes * MeSH
• Guanine chemistry   MeSH
• Ligands MeSH
• Models, Molecular MeSH
• Computer Simulation MeSH
• Thermodynamics MeSH
• Publication type
• Journal Article MeSH
• Research Support, Non-U.S. Gov't MeSH
• Names of Substances
• DNA MeSH
• Guanine MeSH
• Ligands MeSH

A computational analysis of d(GGGGTTTTGGGG)(2) guanine quadruplexes containing either lateral or diagonal four-thymidine loops was carried out using molecular dynamics (MD) simulations in explicit solvent, locally enhanced sampling (LES) simulations, systematic conformational search, and free energy molecular-mechanics, Poisson Boltzmann, surface area (MM-PBSA) calculations with explicit inclusion of structural monovalent cations. The study provides, within the approximations of the applied all-atom additive force field, a qualitatively complete analysis of the available loop conformational space. The results are independent of the starting structures. Major conformational transitions not seen in conventional MD simulations are observed when LES is applied. The favored LES structures consistently provide lower free energies (as estimated by molecular-mechanics, Poisson Boltzmann, surface area) than other structures. Unfortunately, the predicted optimal structure for the diagonal loop arrangement differs substantially from the atomic resolution experiments. This result is attributed to force field deficiencies, such as the potential misbalance between solute-cation and solvent-cation terms. The MD simulations are unable to maintain the stable coordination of the monovalent cations inside the diagonal loops as reported in a recent x-ray study. The optimal diagonal and lateral loop arrangements appear to be close in energy although a proper inclusion of the loop monovalent cations could stabilize the diagonal architecture.

• MeSH
• Guanine chemistry   MeSH
• Nucleic Acid Conformation * MeSH
• Models, Molecular * MeSH
• Computer Simulation * MeSH
• Thymidine chemistry   MeSH
• Publication type
• Journal Article MeSH
• Research Support, Non-U.S. Gov't MeSH
• Names of Substances
• Guanine MeSH
• Thymidine MeSH

The formation of a cation-stabilized guanine quadruplex (G-DNA) stem is an exceptionally slow process involving complex kinetics that has not yet been characterized at atomic resolution. Here, we investigate the formation of a parallel stranded G-DNA stem consisting of four strands of d(GGGG) using molecular dynamics simulations with explicit inclusion of counterions and solvent. Due to the limitations imposed by the nanosecond timescale of the simulations, rather than watching for the spontaneous formation of G-DNA, our approach probes the stability of possible supramolecular intermediates (including two-, three-, and four-stranded assemblies with out-of-register base pairing between guanines) on the formation pathway. The simulations suggest that "cross-like" two-stranded assemblies may serve as nucleation centers in the initial formation of parallel stranded G-DNA quadruplexes, proceeding through a series of rearrangements involving trapping of cations, association of additional strands, and progressive slippage of strands toward the full stem. To supplement the analysis, approximate free energies of the models are obtained with explicit consideration of the integral cations. The approach applied here serves as a prototype for qualitatively investigating other G-DNA molecules using molecular dynamics simulation and free-energy analysis.

• MeSH
• Time Factors MeSH
• DNA chemistry   MeSH
• G-Quadruplexes MeSH
• Guanine chemistry   MeSH
• Ions MeSH
• Cations MeSH
• Kinetics MeSH
• Nucleic Acid Conformation MeSH
• Molecular Conformation MeSH
• Models, Molecular MeSH
• Oligonucleotides chemistry   MeSH
• Sodium chemistry   MeSH
• Software MeSH
• Temperature MeSH
• Thermodynamics MeSH
• Hydrogen Bonding MeSH
• Publication type
• Journal Article MeSH
• Research Support, Non-U.S. Gov't MeSH
• Research Support, U.S. Gov't, Non-P.H.S. MeSH
• Names of Substances
• DNA MeSH
• Guanine MeSH
• Ions MeSH
• Cations MeSH
• Oligonucleotides MeSH
• Sodium MeSH

Explicit solvent and counterion molecular dynamics simulations have been carried out for a total of >80 ns on the bacterial and spinach chloroplast 5S rRNA Loop E motifs. The Loop E sequences form unique duplex architectures composed of seven consecutive non-Watson-Crick basepairs. The starting structure of spinach chloroplast Loop E was modeled using isostericity principles, and the simulations refined the geometries of the three non-Watson-Crick basepairs that differ from the consensus bacterial sequence. The deep groove of Loop E motifs provides unique sites for cation binding. Binding of Mg(2+) rigidifies Loop E and stabilizes its major groove at an intermediate width. In the absence of Mg(2+), the Loop E motifs show an unprecedented degree of inner-shell binding of monovalent cations that, in contrast to Mg(2+), penetrate into the most negative regions inside the deep groove. The spinach chloroplast Loop E shows a marked tendency to compress its deep groove compared with the bacterial consensus. Structures with a narrow deep groove essentially collapse around a string of Na(+) cations with long coordination times. The Loop E non-Watson-Crick basepairing is complemented by highly specific hydration sites ranging from water bridges to hydration pockets hosting 2 to 3 long-residing waters. The ordered hydration is intimately connected with RNA local conformational variations.

• MeSH
• RNA, Bacterial chemistry   MeSH
• Base Pair Mismatch MeSH
• Nucleic Acid Denaturation MeSH
• Species Specificity MeSH
• Escherichia coli chemistry   MeSH
• Magnesium chemistry   MeSH
• Nucleic Acid Conformation MeSH
• Macromolecular Substances MeSH
• Models, Molecular * MeSH
• Base Pairing * MeSH
• Computer Simulation MeSH
• Motion MeSH
• RNA, Ribosomal, 5S chemistry   MeSH
• RNA, Plant chemistry   MeSH
• Solvents chemistry   MeSH
• Sodium chemistry   MeSH
• Spinacia oleracea chemistry   MeSH
• Binding Sites MeSH
• Water chemistry   MeSH
• Publication type
• Journal Article MeSH
• Evaluation Study MeSH
• Research Support, Non-U.S. Gov't MeSH
• Research Support, U.S. Gov't, Non-P.H.S. MeSH
• Research Support, U.S. Gov't, P.H.S. MeSH
• Comparative Study MeSH
• Validation Study MeSH
• Names of Substances
• RNA, Bacterial MeSH
• Magnesium MeSH
• Macromolecular Substances MeSH
• RNA, Ribosomal, 5S MeSH
• RNA, Plant MeSH
• Solvents MeSH
• Sodium MeSH
• Water MeSH

Unrestrained 5-20-ns explicit-solvent molecular dynamics simulations using the Cornell et al. force field have been carried out for d[GCG(N)11GCG]2 (N, purine base) considering guanine*cytosine (G*C), adenine*thymine (A*T), inosine*5-methyl-cytosine (I*mC), and 2-amino-adenine*thymine (D*T) basepairs. The simulations unambiguously show that the structure and elasticity of N-tracts is primarily determined by the presence of the amino group in the minor groove. Simulated A-, I-, and AI-tracts show almost identical structures, with high propeller twist and minor groove narrowing. G- and D-tracts have small propeller twisting and are partly shifted toward the A-form. The elastic properties also differ between the two groups. The sequence-dependent electrostatic component of base stacking seems to play a minor role. Our conclusions are entirely consistent with available experimental data. Nevertheless, the propeller twist and helical twist in the simulated A-tract appear to be underestimated compared to crystallographic studies. To obtain further insight into the possible force field deficiencies, additional multiple simulations have been made for d(A)10, systematically comparing four major force fields currently used in DNA simulations and utilizing B and A-DNA forms as the starting structure. This comparison shows that the conclusions of the present work are not influenced by the force field choice.

• MeSH
• DNA chemistry   MeSH
• Nucleic Acid Conformation MeSH
• Models, Molecular MeSH
• Base Pairing MeSH
• Polydeoxyribonucleotides chemistry   MeSH
• Elasticity MeSH
• Purines chemistry   MeSH
• Hydrogen Bonding MeSH
• Publication type
• Journal Article MeSH
• Research Support, Non-U.S. Gov't MeSH
• Names of Substances
• DNA MeSH
• Polydeoxyribonucleotides MeSH
• Purines MeSH

A two-dimensional, quantitative J-correlation NMR experiment for precise measurements of the proton-carbon vicinal coupling constants 3J(C2)/4-H1' and 3J(C6)/8-H1' in uniformly 13C-labeled nucleic acids is presented. To reduce loss of signal due to 1H-13C dipole-dipole relaxation, a multiple-quantum constant time experiment with appropriately incorporated band selective 1H and 13C pulses was applied. The experiment is used to obtain the 3J(C2)/4-H1' and 3J(C6)/8-H1' coupling constants in a uniformly 13C, 15N-labeled [d(G4T4G4)]2 quadruplex. The measured values and glycosidic torsion angles in the G-quadruplex, obtained by restrained molecular dynamics with explicit solvent using the previously published restraints, along with selected data from the literature are used to check and modify existing parameters of the Karplus equations. The parameterizations obtained using glycosidic torsion angles derived from the original solution and recently determined X-ray structures are also compared.

• MeSH
• DNA chemistry   MeSH
• Mathematics MeSH
• Models, Molecular MeSH
• Nuclear Magnetic Resonance, Biomolecular * MeSH
• Oligodeoxyribonucleotides chemistry   MeSH
• Carbon Radioisotopes MeSH
• Solvents MeSH
• Publication type
• Journal Article MeSH
• Research Support, Non-U.S. Gov't MeSH
• Research Support, U.S. Gov't, Non-P.H.S. MeSH
• Research Support, U.S. Gov't, P.H.S. MeSH
• Names of Substances
• DNA MeSH
• Oligodeoxyribonucleotides MeSH
• Carbon Radioisotopes MeSH
• Solvents MeSH