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
• denaturace nukleových kyselin MeSH
• draslík chemie   MeSH
• G-kvadruplexy * MeSH
• kationty MeSH
• párování bází MeSH
• promotorové oblasti (genetika) * MeSH
• protoonkogenní proteiny c-kit genetika   MeSH
• simulace molekulární dynamiky MeSH
• sodík chemie   MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Názvy látek
• draslík MeSH
• kationty MeSH
• protoonkogenní proteiny c-kit MeSH
• sodík 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 chemie   MeSH
• G-kvadruplexy * MeSH
• jednovláknová DNA chemie   MeSH
• lidé MeSH
• simulace molekulární dynamiky * MeSH
• telomery chemie   MeSH
• Check Tag
• lidé MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Názvy látek
• DNA MeSH
• jednovláknová DNA 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 chemie   MeSH
• G-kvadruplexy * MeSH
• guanin chemie   MeSH
• konformace nukleové kyseliny MeSH
• ligandy MeSH
• RNA chemie   MeSH
• rozpouštědla chemie   MeSH
• simulace molekulární dynamiky * MeSH
• telomery chemie   MeSH
• vodíková vazba MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• přehledy MeSH
• Research Support, N.I.H., Extramural MeSH
• Research Support, U.S. Gov't, Non-P.H.S. MeSH
• Názvy látek
• DNA MeSH
• guanin MeSH
• ligandy MeSH
• RNA MeSH
• rozpouštědla 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
• glukosa-6-fosfát analogy a deriváty   metabolismus   MeSH
• glukosamin analogy a deriváty   metabolismus   MeSH
• glutaminfruktosa-6-fosfáttransaminasa (izomerizující) genetika   MeSH
• katalytická doména * MeSH
• koenzymy metabolismus   MeSH
• konformace nukleové kyseliny MeSH
• molekulární sekvence - údaje MeSH
• protony * MeSH
• RNA katalytická chemie   genetika   metabolismus   MeSH
• sekvence nukleotidů MeSH
• simulace molekulární dynamiky * MeSH
• Thermoanaerobacter enzymologie   MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Research Support, N.I.H., Extramural MeSH
• Názvy látek
• glucosamine 6-phosphate MeSH Prohlížeč
• glukosa-6-fosfát MeSH
• glukosamin MeSH
• glutaminfruktosa-6-fosfáttransaminasa (izomerizující) MeSH
• koenzymy MeSH
• protony * MeSH
• RNA katalytická 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 chemie   MeSH
• G-kvadruplexy * MeSH
• guanin chemie   MeSH
• ligandy MeSH
• molekulární modely MeSH
• počítačová simulace MeSH
• termodynamika MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Názvy látek
• DNA MeSH
• guanin MeSH
• ligandy 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
• guanin chemie   MeSH
• konformace nukleové kyseliny * MeSH
• molekulární modely * MeSH
• počítačová simulace * MeSH
• thymidin chemie   MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Názvy látek
• guanin MeSH
• thymidin 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
• časové faktory MeSH
• DNA chemie   MeSH
• G-kvadruplexy MeSH
• guanin chemie   MeSH
• ionty MeSH
• kationty MeSH
• kinetika MeSH
• konformace nukleové kyseliny MeSH
• molekulární konformace MeSH
• molekulární modely MeSH
• oligonukleotidy chemie   MeSH
• sodík chemie   MeSH
• software MeSH
• teplota MeSH
• termodynamika MeSH
• vodíková vazba MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Research Support, U.S. Gov't, Non-P.H.S. MeSH
• Názvy látek
• DNA MeSH
• guanin MeSH
• ionty MeSH
• kationty MeSH
• oligonukleotidy MeSH
• sodík 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
• bakteriální RNA chemie   MeSH
• chybné párování bází MeSH
• denaturace nukleových kyselin MeSH
• druhová specificita MeSH
• Escherichia coli chemie   MeSH
• hořčík chemie   MeSH
• konformace nukleové kyseliny MeSH
• makromolekulární látky MeSH
• molekulární modely * MeSH
• párování bází * MeSH
• počítačová simulace MeSH
• pohyb těles MeSH
• RNA ribozomální 5S chemie   MeSH
• RNA rostlin chemie   MeSH
• rozpouštědla chemie   MeSH
• sodík chemie   MeSH
• Spinacia oleracea chemie   MeSH
• vazebná místa MeSH
• voda chemie   MeSH
• Publikační typ
• časopisecké články MeSH
• hodnotící studie MeSH
• práce podpořená grantem MeSH
• Research Support, U.S. Gov't, Non-P.H.S. MeSH
• Research Support, U.S. Gov't, P.H.S. MeSH
• srovnávací studie MeSH
• validační studie MeSH
• Názvy látek
• bakteriální RNA MeSH
• hořčík MeSH
• makromolekulární látky MeSH
• RNA ribozomální 5S MeSH
• RNA rostlin MeSH
• rozpouštědla MeSH
• sodík MeSH
• voda 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 chemie   MeSH
• matematika MeSH
• molekulární modely MeSH
• nukleární magnetická rezonance biomolekulární * MeSH
• oligodeoxyribonukleotidy chemie   MeSH
• radioizotopy uhlíku MeSH
• rozpouštědla MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Research Support, U.S. Gov't, Non-P.H.S. MeSH
• Research Support, U.S. Gov't, P.H.S. MeSH
• Názvy látek
• DNA MeSH
• oligodeoxyribonukleotidy MeSH
• radioizotopy uhlíku MeSH
• rozpouštědla 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 chemie   MeSH
• konformace nukleové kyseliny MeSH
• molekulární modely MeSH
• párování bází MeSH
• polydeoxyribonukleotidy chemie   MeSH
• pružnost MeSH
• puriny chemie   MeSH
• vodíková vazba MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Názvy látek
• DNA MeSH
• polydeoxyribonukleotidy MeSH
• puriny MeSH