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

Genomic sequences susceptible to form G-quadruplexes (G4s) are always flanked by other nucleotides, but G4 formation in vitro is generally studied with short synthetic DNA or RNA oligonucleotides, for which bases adjacent to the G4 core are often omitted. Herein, we systematically studied the effects of flanking nucleotides on structural polymorphism of 371 different oligodeoxynucleotides that adopt intramolecular G4 structures. We found out that the addition of nucleotides favors the formation of a parallel fold, defined as the 'flanking effect' in this work. This 'flanking effect' was more pronounced when nucleotides were added at the 5'-end, and depended on loop arrangement. NMR experiments and molecular dynamics simulations revealed that flanking sequences at the 5'-end abolish a strong syn-specific hydrogen bond commonly found in non-parallel conformations, thus favoring a parallel topology. These analyses pave a new way for more accurate prediction of DNA G4 folding in a physiological context.

• MeSH
• cirkulární dichroismus MeSH
• DNA genetika   ultrastruktura   MeSH
• G-kvadruplexy * MeSH
• konformace nukleové kyseliny MeSH
• nukleotidy chemie   genetika   MeSH
• oligonukleotidy chemie   genetika   MeSH
• polymorfismus genetický genetika   MeSH
• RNA genetika   ultrastruktura   MeSH
• simulace molekulární dynamiky MeSH
• vodíková vazba MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Názvy látek
• DNA MeSH
• nukleotidy MeSH
• oligonukleotidy MeSH
• RNA MeSH

With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.

• MeSH
• DNA chemie   MeSH
• katalýza MeSH
• konformace nukleové kyseliny * MeSH
• počítačová simulace MeSH
• RNA chemie   MeSH
• simulace molekulární dynamiky * MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• přehledy MeSH
• Research Support, N.I.H., Extramural MeSH
• Názvy látek
• DNA MeSH
• RNA MeSH

We have carried out a series of extended unbiased molecular dynamics (MD) simulations (up to 10 μs long, ∼162 μs in total) complemented by replica-exchange with the collective variable tempering (RECT) approach for several human telomeric DNA G-quadruplex (GQ) topologies with TTA propeller loops. We used different AMBER DNA force-field variants and also processed simulations by Markov State Model (MSM) analysis. The slow conformational transitions in the propeller loops took place on a scale of a few μs, emphasizing the need for long simulations in studies of GQ dynamics. The propeller loops sampled similar ensembles for all GQ topologies and for all force-field dihedral-potential variants. The outcomes of standard and RECT simulations were consistent and captured similar spectrum of loop conformations. However, the most common crystallographic loop conformation was very unstable with all force-field versions. Although the loss of canonical γ-trans state of the first propeller loop nucleotide could be related to the indispensable bsc0 α/γ dihedral potential, even supporting this particular dihedral by a bias was insufficient to populate the experimentally dominant loop conformation. In conclusion, while our simulations were capable of providing a reasonable albeit not converged sampling of the TTA propeller loop conformational space, the force-field description still remained far from satisfactory.

• MeSH
• DNA chemie   genetika   metabolismus   MeSH
• G-kvadruplexy * MeSH
• lidé MeSH
• sekvence nukleotidů MeSH
• shluková analýza MeSH
• simulace molekulární dynamiky * MeSH
• telomery genetika   MeSH
• voda chemie   MeSH
• Check Tag
• lidé MeSH
• Publikační typ
• časopisecké články MeSH
• Názvy látek
• DNA MeSH
• voda MeSH

DNA G-hairpins are potential key structures participating in folding of human telomeric guanine quadruplexes (GQ). We examined their properties by standard MD simulations starting from the folded state and long T-REMD starting from the unfolded state, accumulating ∼130 μs of atomistic simulations. Antiparallel G-hairpins should spontaneously form in all stages of the folding to support lateral and diagonal loops, with sub-μs scale rearrangements between them. We found no clear predisposition for direct folding into specific GQ topologies with specific syn/anti patterns. Our key prediction stemming from the T-REMD is that an ideal unfolded ensemble of the full GQ sequence populates all 4096 syn/anti combinations of its four G-stretches. The simulations can propose idealized folding pathways but we explain that such few-state pathways may be misleading. In the context of the available experimental data, the simulations strongly suggest that the GQ folding could be best understood by the kinetic partitioning mechanism with a set of deep competing minima on the folding landscape, with only a small fraction of molecules directly folding to the native fold. The landscape should further include non-specific collapse processes where the molecules move via diffusion and consecutive random rare transitions, which could, e.g. structure the propeller loops.

• MeSH
• DNA chemie   MeSH
• G-kvadruplexy * MeSH
• kationty chemie   MeSH
• lidé MeSH
• Oxytricha genetika   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
• kationty MeSH

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