Guanine quadruplex (GQ) is a noncanonical nucleic acid structure formed by guanine-rich DNA and RNA sequences. Folding of GQs is a complex process, where several aspects remain elusive, despite being important for understanding structure formation and biological functions of GQs. Pulling experiments are a common tool for acquiring insights into the folding landscape of GQs. Herein, we applied a computational pulling strategy─steered molecular dynamics (SMD) simulations─in combination with standard molecular dynamics (MD) simulations to explore the unfolding landscapes of tetrameric parallel GQs. We identified anisotropic properties of elastic conformational changes, unfolding transitions, and GQ mechanical stabilities. Using a special set of structural parameters, we found that the vertical component of pulling force (perpendicular to the average G-quartet plane) plays a significant role in disrupting GQ structures and weakening their mechanical stabilities. We demonstrated that the magnitude of the vertical force component depends on the pulling anchor positions and the number of G-quartets. Typical unfolding transitions for tetrameric parallel GQs involve base unzipping, opening of the G-stem, strand slippage, and rotation to cross-like structures. The unzipping was detected as the first and dominant unfolding event, and it usually started at the 3'-end. Furthermore, results from both SMD and standard MD simulations indicate that partial spiral conformations serve as a transient ensemble during the (un)folding of GQs.

• MeSH
• biomechanika MeSH
• DNA chemie   MeSH
• G-kvadruplexy * MeSH
• mechanické jevy MeSH
• simulace molekulární dynamiky * MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Názvy látek
• DNA 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

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

We show using polyacrylamide gel electrophoresis that guanine+adenine repeat strands of DNA associate into homoduplexes at neutral pH and moderate ionic strength. The homoduplexes melt in a cooperative way like the Watson-Crick duplex, although they contain no Watson-Crick base pair. Guanine is absolutely needed for the homoduplex formation and the homoduplex stability increases with the guanine content of the repeat. The present results have implications for the nature of the first replicators, as well as regarding forces stabilizing the duplexes of DNA.

• MeSH
• adenin chemie   MeSH
• cirkulární dichroismus * MeSH
• DNA chemie   MeSH
• elektroforéza v polyakrylamidovém gelu metody   MeSH
• guanin chemie   MeSH
• konformace nukleové kyseliny MeSH
• nukleotidy chemie   MeSH
• oligonukleotidy chemie   MeSH
• párování bází * MeSH
• pyrimidiny chemie   MeSH
• tandemové repetitivní sekvence * MeSH
• teplota MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Názvy látek
• adenin MeSH
• DNA MeSH
• guanin MeSH
• nukleotidy MeSH
• oligonukleotidy MeSH
• pyrimidine MeSH Prohlížeč
• pyrimidiny MeSH