Coarse-Grained Simulations Complemented by Atomistic Molecular Dynamics Provide New Insights into Folding and Unfolding of Human Telomeric G-Quadruplexes

Molecular dynamics simulations reveal the parallel stranded d(GGGA)3GGG DNA quadruplex folds via multiple paths from a coil-like ensemble

G-quadruplexes (G4s) are non-canonical nucleic acid structures that fold through complex processes. Characterization of the G4 folding landscape may help to elucidate biological roles of G4s but is challenging both experimentally and computationally. Here, we achieved complete folding of a three-quartet parallel DNA G4 with (GGGA)3GGG sequence using all-atom explicit-solvent enhanced-sampling molecular dynamics (MD) simulations. The simulations suggested early formation of guanine stacks in the G-tracts, which behave as semi-rigid blocks in the folding process. The folding continues via the formation of a collapsed compact coil-like ensemble. Structuring of the G4 from the coil then proceeds via various cross-like, hairpin, slip-stranded and two-quartet ensembles and can bypass the G-triplex structure. Folding of the parallel G4 does not appear to involve any salient intermediates and is a multi-pathway process. We also carried out an extended set of simulations of parallel G-hairpins. While parallel G-hairpins are extremely unstable when isolated, they are more stable inside the coil structure. On the methodology side, we show that the AMBER DNA force field predicts the folded G4 to be less stable than the unfolded ensemble, uncovering substantial force-field issues. Overall, we provide unique atomistic insights into the folding landscape of parallel-stranded G4 but also reveal limitations of current state-of-the-art MD techniques.

RNA G-quadruplex folding is a multi-pathway process driven by conformational entropy

The kinetics of folding is crucial for the function of many regulatory RNAs including RNA G-quadruplexes (rG4s). Here, we characterize the folding pathways of a G-quadruplex from the telomeric repeat-containing RNA by combining all-atom molecular dynamics and coarse-grained simulations with circular dichroism experiments. The quadruplex fold is stabilized by cations and thus, the ion atmosphere forming a double layer surrounding the highly charged quadruplex guides the folding process. To capture the ionic double layer in implicit solvent coarse-grained simulations correctly, we develop a matching procedure based on all-atom simulations in explicit water. The procedure yields quantitative agreement between simulations and experiments as judged by the populations of folded and unfolded states at different salt concentrations and temperatures. Subsequently, we show that coarse-grained simulations with a resolution of three interaction sites per nucleotide are well suited to resolve the folding pathways and their intermediate states. The results reveal that the folding progresses from unpaired chain via hairpin, triplex and double-hairpin constellations to the final folded structure. The two- and three-strand intermediates are stabilized by transient Hoogsteen interactions. Each pathway passes through two on-pathway intermediates. We hypothesize that conformational entropy is a hallmark of rG4 folding. Conformational entropy leads to the observed branched multi-pathway folding process for TERRA25. We corroborate this hypothesis by presenting the free energy landscapes and folding pathways of four rG4 systems with varying loop length.

Studying the Dynamics of a Complex G-Quadruplex System: Insights into the Comparison of MD and NMR Data

Molecular dynamics (MD) simulations are coming of age in the study of nucleic acids, including specific tertiary structures such as G-quadruplexes. While being precious for providing structural and dynamic information inaccessible to experiments at the atomistic level of resolution, MD simulations in this field may still be limited by several factors. These include the force fields used, different models for ion parameters, ionic strengths, and water models. We address various aspects of this problem by analyzing and comparing microsecond-long atomistic simulations of the G-quadruplex structure formed by the human immunodeficiency virus long terminal repeat (HIV LTR)-III sequence for which nuclear magnetic resonance (NMR) structures are available. The system is studied in different conditions, systematically varying the ionic strengths, ion numbers, and water models. We comparatively analyze the dynamic behavior of the G-quadruplex motif in various conditions and assess the ability of each simulation to satisfy the nuclear magnetic resonance (NMR)-derived experimental constraints and structural parameters. The conditions taking into account K⁺-ions to neutralize the system charge, mimicking the intracellular ionic strength, and using the four-atom water model are found to be the best in reproducing the experimental NMR constraints and data. Our analysis also reveals that in all of the simulated environments residues belonging to the duplex moiety of HIV LTR-III exhibit the highest flexibility.

Computer-aided comprehensive explorations of RNA structural polymorphism through complementary simulation methods

While RNA folding was originally seen as a simple problem to solve, it has been shown that the promiscuous interactions of the nucleobases result in structural polymorphism, with several competing structures generally observed for non-coding RNA. This inherent complexity limits our understanding of these molecules from experiments alone, and computational methods are commonly used to study RNA. Here, we discuss three advanced sampling schemes, namely Hamiltonian-replica exchange molecular dynamics (MD), ratchet-and-pawl MD and discrete path sampling, as well as the HiRE-RNA coarse-graining scheme, and highlight how these approaches are complementary with reference to recent case studies. While all computational methods have their shortcomings, the plurality of simulation methods leads to a better understanding of experimental findings and can inform and guide experimental work on RNA polymorphism.

Biasing RNA Coarse-Grained Folding Simulations with Small-Angle X-ray Scattering Data

RNA molecules can easily adopt alternative structures in response to different environmental conditions. As a result, a molecule's energy landscape is rough and can exhibit a multitude of deep basins. In the absence of a high-resolution structure, small-angle X-ray scattering data (SAXS) can narrow down the conformational space available to the molecule and be used in conjunction with physical modeling to obtain high-resolution putative structures to be further tested by experiments. Because of the low resolution of these data, it is natural to implement the integration of SAXS data into simulations using a coarse-grained representation of the molecule, allowing for much wider searches and faster evaluation of SAXS theoretical intensity curves than with atomistic models. We present here the theoretical framework and the implementation of a simulation approach based on our coarse-grained model HiRE-RNA combined with SAXS evaluations "on-the-fly" leading the simulation toward conformations agreeing with the scattering data, starting from partially folded structures as the ones that can easily be obtained from secondary structure prediction-based tools. We show on three benchmark systems how our approach can successfully achieve high-resolution structures with remarkable similarity with the native structure recovering not only the overall shape, as imposed by SAXS data, but also the details of initially missing base pairs.

RNA Modeling with the Computational Energy Landscape Framework

The recent advances in computational abilities, such as the enormous speed-ups provided by GPU computing, allow for large scale computational studies of RNA molecules at an atomic level of detail. As RNA molecules are known to adopt multiple conformations with comparable energies, but different two-dimensional structures, all-atom models are necessary to better describe the structural ensembles for RNA molecules. This point is important because different conformations can exhibit different functions, and their regulation or mis-regulation is linked to a number of diseases. Problematically, the energy barriers between different conformational ensembles are high, resulting in long time scales for interensemble transitions. The computational potential energy landscape framework was designed to overcome this problem of broken ergodicity by use of geometry optimization. Here, we describe the algorithms used in the energy landscape explorations with the OPTIM and PATHSAMPLE programs, and how they are used in biomolecular simulations. We present a recent case study of the 5'-hairpin of RNA 7SK to illustrate how the method can be applied to interpret experimental results, and to obtain a detailed description of molecular properties.

Theoretical investigation of conformational deviation of the human parallel telomeric G-quadruplex DNA in the presence of different salt concentrations and temperatures under confinement

Various experimental reports address the stability of G-quadruplex DNA inside a close confinement such as α-hemolysin, nanocavity water pool and different metal-organic-frameworks (MOFs). To understand the conformational change of G-quadruplex DNA at the atomistic level, we have carried out a total of 40 μs simulation run under both non-polar and polar confinement conditions. To investigate the dynamics, we have considered two different KCl salt concentrations, i.e., 0.47 M (minimal salt concentration) and higher than 2 M (higher salt concentration), at two distinct temperatures, 300 K and 350 K. Here, we have observed that the human telomeric G-quadruplex DNA deviates more from its crystal structure at minimal salt concentration under both non-polar and polar confinement conditions. Besides, the loop regions deviate and fluctuate more compared to the other regions, i.e., sugar-phosphate backbone and tetrad regions. The presence of K+ ions is found to be primarily responsible for this phenomenon. From the spatial density function (SDF) plots, a higher density of K+ ions is observed in the backbone region. Furthermore, from the residue-wise first solvation shell estimation, we have noticed that the K+ ions mainly accumulate in the tetrad region under both non-polar and polar confinement conditions due to which the tetrad regions are more rigid than the loop regions. Higher salt concentration results in increased rigidity of the G-quadruplex DNA. Our study provides valuable insight into the conformational deviation of the G-quadruplex DNA under nanoconfinement conditions.

Folding Dynamics of DNA G-Quadruplexes Probed by Millisecond Temperature Jump Circular Dichroism

G-quadruplexes play important roles in cellular regulatory functions, but despite significant experimental and theoretical efforts, their folding mechanisms remain poorly understood. In this context, we developed a T-jump experiment to access the thermal denaturation and renaturation dynamics of short intramolecular G-quadruplexes in vitro, on the time scale of a few hundred milliseconds. With this new setup, we compared the thermal denaturation and renaturation kinetics of three antiparallel topologies made of the human telomeric sequences d[(5'-GGG(TTAGGG)₃-3']/Na⁺ and d[5'-AGGG(TTAGGG)₃-3']/Na⁺ and the thrombin-binding aptamer sequence d[5'-GGTTGGTGTGGTTGG-3']/K⁺, with those of the parallel topology made of the human CEB25 minisatellite d[5'-AAGGGTGGGTGTAAGTGTGGGTGGGT-3']/Na⁺. In all cases, exponential kinetics of the order of several hundred milliseconds were observed. Measurements performed for different initial temperatures revealed distinct denaturation and renaturation dynamics, ruling out a simple two-state mechanism. The parallel topology, in which all guanines adopt an anti conformation, displays much slower dynamics than antiparallel topologies associated with very low activation barriers. This behavior can be explained by the constrained conformational space due to the presence of the single-base propeller loops that likely hinders the movement of the coiled DNA strand and reduces the contribution of the entropy during the renaturation process at high temperatures.

Learning to Model G-Quadruplexes: Current Methods and Perspectives

G-quadruplexes have raised considerable interest during the past years for the development of therapies against cancer. These noncanonical structures of DNA may be found in telomeres and/or oncogene promoters, and it has been observed that the stabilization of such G-quadruplexes may disturb tumor cell growth. Nevertheless, the mechanisms leading to folding and stabilization of these G-quadruplexes are still not well established, and they are the focus of much current work in this field. In seminal works, stabilization was observed to be produced by cations. However, subsequent studies showed that different kinds of small molecules, from planar and nonplanar organic molecules to square-planar and octahedral metal complexes, may also lead to the stabilization of G-quadruplexes. Thus, the comprehension and rationalization of the interaction of these small molecules with G-quadruplexes are also important topics of current interest in medical applications. To shed light on the questions arising from the literature on the formation of G-quadruplexes, their stabilization, and their interaction with small molecules, synergies between experimental studies and computational works are needed. In this review, we mainly focus on in silico approaches and provide a broad compilation of different leading studies carried out to date by different computational methods. We divide these methods into twomain categories: (a) classical methods, which allow for long-timescale molecular dynamics simulations and the corresponding analysis of dynamical information, and (b) quantum methods (semiempirical, quantum mechanics/molecular mechanics, and density functional theory methods), which allow for the explicit simulation of the electronic structure of the system but, in general, are not capable of being used in long-timescale molecular dynamics simulations and, therefore, give a more static picture of the relevant processes.

An in silico investigation of the binding modes and pathway of APTO-253 on c-KIT G-quadruplex DNA

The stability of c-KIT G-quadruplex DNA via ligands has been a significant concern in the growing field of cancer therapy. Thus, it is very important to understand the mechanism behind the high binding affinity of the small drug molecules on the c-KIT G-quadruplex DNA. In this study, we have investigated the binding mode and pathway of the APTO-253 ligand on the c-KIT G-quadruplex DNA employing a total of 10 μs all atom molecular dynamics simulations and further 8.82 μs simulations via the umbrella sampling method using both OL15 and BSC1 latest force fields for DNA structures. From the cluster structure analysis, mainly three binding pathways i.e., top, bottom and side loop stacking modes are identified. Moreover, RMSD, RMSF and 2D-RMSD values indicate that the c-KIT G-quadruplex DNA and APTO-253 molecules are stable throughout the simulation run. Furthermore, the number of hydrogen bonds in each tetrad and the distance between the two central K⁺ cations confirm that the c-KIT G-quadruplex DNA maintains its conformation in the process of complex formation with the APTO-253 ligand. The binding free energies and the minimum values in the potential of mean forces suggest that the binding processes are energetically favorable. Furthermore, we have found that the bottom stacking mode is the most favorable binding mode among all the three modes for the OL15 force field. However, for the BSC1 force field, both the top and bottom binding modes of the APTO-253 ligand in c-KIT G-quadruplex DNA are comparable to each other. To investigate the driving force for the complex formation, we have noticed that the van der Waals (vdW) and π-π stacking interactions are mainly responsible. Our detailed studies provide useful information for the discovery of novel drugs in the field of stabilization of G-quadruplex DNAs.

RNAxplorer: Harnessing the Power of Guiding Potentials to Sample RNA Landscapes

Motivation

Predicting the folding dynamics of RNAs is a computationally difficult problem, first and foremost due to the combinatorial explosion of alternative structures in the folding space. Abstractions are therefore needed to simplify downstream analyses, and thus make them computationally tractable. This can be achieved by various structure sampling algorithms. However, current sampling methods are still time consuming and frequently fail to represent key elements of the folding space.

Method

We introduce RNAxplorer, a novel adaptive sampling method to efficiently explore the structure space of RNAs. RNAxplorer uses dynamic programming to perform an efficient Boltzmann sampling in the presence of guiding potentials, which are accumulated into pseudo-energy terms and reflect similarity to already well-sampled structures. This way, we effectively steer sampling toward underrepresented or unexplored regions of the structure space.

Results

We developed and applied different measures to benchmark our sampling methods against its competitors. Most of the measures show that RNAxplorer produces more diverse structure samples, yields rare conformations that may be inaccessible to other sampling methods and is better at finding the most relevant kinetic traps in the landscape. Thus, it produces a more representative coarse graining of the landscape, which is well suited to subsequently compute better approximations of RNA folding kinetics.

Availabilityand implementation

https://github.com/ViennaRNA/RNAxplorer/.

Supplementary information

Supplementary data are available at Bioinformatics online.

Insights into the Kinetic Partitioning Folding Dynamics of the Human Telomeric G-Quadruplex from Molecular Simulations and Machine Learning

The human telomeric DNA G-quadruplex follows a kinetic partitioning folding mechanism. The underlying folding landscape potentially has many minima separated by high free-energy barriers. However, using current theoretical models to characterize this complex folding landscape has remained a challenging problem. In this study, by developing a hybrid atomistic structure-based model that merges structural information on the hybrid-1, hybrid-2, and chair-type G-quadruplex topologies, we investigated a kinetic partitioning folding process of human telomeric DNA involving three native folds. The model was validated as it reproduced the experimental observation that the hybrid-1 conformation is the major fold and the hybrid-2 conformation is kinetically more accessible. A three-step mechanism was revealed for the formation of the hybrid-1 conformation, while a two-step mechanism was demonstrated for the formation of hybrid-2 and chair-type conformations. Likewise, a class of state in which structures adopted inappropriate combinations of syn/anti guanine nucleotides was found to greatly slow down the folding process. In addition, by employing the XGBoost machine learning algorithm, three interatom distances and six dihedral angles were identified as essential internal coordinates to represent the low-dimensional folding landscape. The strategy of coupling the multibasin model and the machine learning algorithm may be useful to investigate the conformational dynamics of other multistate biomolecules.

Delicate Balance of Non-Covalent Forces Govern the Biocompatibility of Graphitic Carbon Nitride towards Genetic Materials

Despite a plethora of suggested technological and biomedical applications, the nanotoxicity of two-dimensional (2D) graphitic carbon nitride (g-C₃ N₄ ) towards biomolecules remains elusive. To address this issue, we employ all-atom classical molecular dynamics simulations and investigate the interactions between nucleic acids and g-C₃ N₄ . It is revealed that, toxicity is modulated through a subtle balance between electrostatic and van der Waals interactions. When the exposed nucleobases interact through predominantly short-ranged van der Waals and π-π stacking interactions, they get deviated from their native disposition and adsorb on the surface, leading to loss of self-stacking and intra-quartet H-bonding along with partial disruption of the native structure. In contrast, for the interaction with double-stranded structures of both DNA and RNA, long-range electrostatics govern the adsorption phenomena since the constituent nucleobases are relatively concealed and wrapped, thereby resulting in almost complete preservation of the nucleic acid structures. Construction of free energy landscapes for lateral translation of adsorbed nucleic acids suggests decent targeting specificity owing to their restricted movement on g-C₃ N₄ . The release times of nucleic acids adsorbed through predominant electrostatics are significantly less than those adsorbed through stacking with the surface. It is therefore proposed that g-C₃ N₄ would induce toxicity towards any biomolecule having bare residues available for strong van der Waals and π-π stacking interactions relative to those predominantly interacting through electrostatics.

Understanding The Role of Reline, a Natural DES, on Temperature-Induced Conformational Changes of C-Kit G-Quadruplex DNA: A Molecular Dynamics Study

The noncanonical guanine-rich DNAs have drawn particular attention to the scientific world due to their controllable diverse and polymorphic structures. Apart from biological and medical significance, G-quadruplex DNAs are widely used in various fields such as nanotechnology, nanomachine, biosensors, and biocatalyst. So far, the applications of the G-quadruplex DNA are mainly limited in the water medium. Recently, a new generation of solvent named deep eutectic solvent (DES) has become very popular and has been widely used as a reaction medium of biocatalytic reactions and long-term storage medium for nucleic acids, even at high temperature. Hence, it is essential to understand the role of DES on temperature-induced conformational changes of a G-quadruplex DNA. In this research work, we have explored the temperature-mediated conformational dynamics of c-kit oncogene promoter G-quadruplex DNA in reline medium in the temperature range of 300-500 K, using a total of 10 μs unbiased all-atom molecular dynamics simulation. Here, from RMSD, RMSF, R_g and principal component analyses, we notice that the c-kit G-quadruplex DNA is stable up to 450 K in reline medium. However, it unfolds in water medium at 450 K. It is found that the hydrogen bonding interactions between c-kit G-quadruplex DNA and reline play a key role in the stabilization of the G-quadruplex DNA even at high temperature. Furthermore, in this work we have observed a very interesting and distinctive phenomenon of the central cation of the G-quadruplex DNA. Its position was seen to fluctuate between the two tetrad cores, that is, the region between tetrad-1 and tetrad-2 and that between tetrad-2 and tetrad-3 and vice versa at 450 and 500 K in reline medium which is absent in water medium at 450 K. Moreover, the rate of its oscillation is increased when temperature is increased.

Structural transitions in the RNA 7SK 5' hairpin and their effect on HEXIM binding

7SK RNA, as part of the 7SK ribonucleoprotein complex, is crucial to the regulation of transcription by RNA-polymerase II, via its interaction with the positive transcription elongation factor P-TEFb. The interaction is induced by binding of the protein HEXIM to the 5′ hairpin (HP1) of 7SK RNA. Four distinct structural models have been obtained experimentally for HP1. Here, we employ computational methods to investigate the relative stability of these structures, transitions between them, and the effects of mutations on the observed structural ensembles. We further analyse the results with respect to mutational binding assays, and hypothesize a mechanism for HEXIM binding. Our results indicate that the dominant structure in the wild type exhibits a triplet involving the unpaired nucleotide U40 and the base pair A43-U66 in the GAUC/GAUC repeat. This conformation leads to an open major groove with enough potential binding sites for peptide recognition. Sequence mutations of the RNA change the relative stability of the different structural ensembles. Binding affinity is consequently lost if these changes alter the dominant structure.

Screening two dimensional materials for the transportation and delivery of diverse genetic materials

In spite of several reports of graphene and other 2D materials concerning their capacity for biomolecular adsorption and delivery, recent toxicity evaluations found them to be nanotoxic toward different biomolecules, especially nucleic acids. Therefore, there is urgent demand for the synthesis of 2D materials exhibiting biocompatible and non-nanotoxic features. In this article, employing classical molecular dynamics simulations, we provide a benchmarking of h2D-C2N, graphene and hexagonal boron nitride (h-BN) toward the adsorption, preservation, targeting and delivery of various classes of nucleic acids namely single stranded DNA, double stranded natural as well as unnatural base substituted DNA and two different types of human telomeric guanine quadruplexes, all comprising different secondary structures. Our simulations reveal that, while h2D-C2N preserves the structures of most of the nucleic acids, graphene and h-BN disrupt them through strong π-π stacking with aromatic nucleobases. Interestingly, for the first time we identified a 'quartet-by-quartet' disruption mechanism of guanine quadruplexes, but only on graphene and h-BN. The lateral diffusion of adsorbed nucleic acids over C2N is restricted unlike that over both graphene and h-BN, thereby increasing the targeting efficacy for C2N. Modeling of the delivery phenomena suggests orders of magnitude longer release times from graphene and h-BN compared to C2N, thereby demonstrating the preferential suitability of C2N for all the hierarchical steps of nucleic acid transportation.

Parallel G-triplexes and G-hairpins as potential transitory ensembles in the folding of parallel-stranded DNA G-Quadruplexes

Guanine quadruplexes (G4s) are non-canonical nucleic acids structures common in important genomic regions. Parallel-stranded G4 folds are the most abundant, but their folding mechanism is not fully understood. Recent research highlighted that G4 DNA molecules fold via kinetic partitioning mechanism dominated by competition amongst diverse long-living G4 folds. The role of other intermediate species such as parallel G-triplexes and G-hairpins in the folding process has been a matter of debate. Here, we use standard and enhanced-sampling molecular dynamics simulations (total length of ∼0.9 ms) to study these potential folding intermediates. We suggest that parallel G-triplex per se is rather an unstable species that is in local equilibrium with a broad ensemble of triplex-like structures. The equilibrium is shifted to well-structured G-triplex by stacked aromatic ligand and to a lesser extent by flanking duplexes or nucleotides. Next, we study propeller loop formation in GGGAGGGAGGG, GGGAGGG and GGGTTAGGG sequences. We identify multiple folding pathways from different unfolded and misfolded structures leading towards an ensemble of intermediates called cross-like structures (cross-hairpins), thus providing atomistic level of description of the single-molecule folding events. In summary, the parallel G-triplex is a possible, but not mandatory short-living (transitory) intermediate in the folding of parallel-stranded G4.

Coarse-grained dynamic RNA titration simulations

Electrostatic interactions play a pivotal role in many biomolecular processes. The molecular organization and function in biological systems are largely determined by these interactions. Owing to the highly negative charge of RNA, the effect is expected to be more pronounced in this system. Moreover, RNA base pairing is dependent on the charge of the base, giving rise to alternative secondary and tertiary structures. The equilibrium between uncharged and charged bases is regulated by the solution pH, which is therefore a key environmental condition influencing the molecule's structure and behaviour. By means of constant-pH Monte Carlo simulations based on a fast proton titration scheme, coupled with the coarse-grained model HiRE-RNA, molecular dynamic simulations of RNA molecules at constant pH enable us to explore the RNA conformational plasticity at different pH values as well as to compute electrostatic properties as local pK _a values for each nucleotide.

Fast-Folding Pathways of the Thrombin-Binding Aptamer G-Quadruplex Revealed by a Markov State Model

G-quadruplex structures participate in many important cellular processes. For a better understanding of their functions, knowledge of the mechanism by which they fold into the functional native structures is necessary. In this work, we studied the folding process of the thrombin-binding aptamer G-quadruplex. Enabled by a computational paradigm that couples an advanced sampling method and a Markov state model, four folding intermediates were identified, including an antiparallel G-hairpin, two G-triplex structures, and a double-hairpin conformation. Likewise, a misfolded structure with a nonnative distribution of syn/anti guanines was also observed. Based on these states, a transition path analysis revealed three fast-folding pathways, along which the thrombin-binding aptamer would fold to the native state directly, with no evidence of potential nonnative competing conformations. The results also showed that the TGT-loop plays an important role in the folding process. The findings of this research may provide general insight about the folding of other G-quadruplex structures.

Conformational deviation of Thrombin binding G-quadruplex aptamer (TBA) in presence of divalent cation Sr2+: A classical molecular dynamics simulation study

Thrombin binding TBA-G-quadruplex aptamer (TBA) plays a major role in blood coagulation cascade. The 15-mer TBA sequence tends to form four-stranded TBA-G-quadruplex structure. In this research work, a series of explicit solvent classical MD simulations of the TBA is carried out using different salt (SrCl₂) concentrations (0, 50, 100 and 200 mM). Here we have also testified the effect of salt concentration of divalent cation Sr²⁺ on the conformational change of quadruplex DNA. The structural deviations, fluctuations, torsional angles and the affinity of the ion are explored at different salt concentrations. It is found that the conformation of TBA-G-quadruplex at 0 mM and 50 mM salt concentrations, is very much different than the other salt concentrations (100 mM and 200 mM). Also observed are as follows: (i) no exchange of Sr²⁺ ion between inside and outside of the channel, (ii) an enhancement in the Sr²⁺ ion density around the phosphate region of the loop residues as salt concentration increases and (iii) the stacking of T3 and T4 residues of loop-1 that appears up to 50 mM concentration, vanishes as the salt concentration is increased further.

Structural dynamics of propeller loop: towards folding of RNA G-quadruplex

We have carried out an extended set of standard and enhanced-sampling MD simulations (for a cumulative simulation time of 620 μs) with the aim to study folding landscapes of the rGGGUUAGGG and rGGGAGGG parallel G-hairpins (PH) with propeller loop. We identify folding and unfolding pathways of the PH, which is bridged with the unfolded state via an ensemble of cross-like structures (CS) possessing mutually tilted or perpendicular G-strands interacting via guanine-guanine H-bonding. The oligonucleotides reach the PH conformation from the unfolded state via a conformational diffusion through the folding landscape, i.e. as a series of rearrangements of the H-bond interactions starting from compacted anti-parallel hairpin-like structures. Although isolated PHs do not appear to be thermodynamically stable we suggest that CS and PH-types of structures are sufficiently populated during RNA guanine quadruplex (GQ) folding within the context of complete GQ-forming sequences. These structures may participate in compact coil-like ensembles that involve all four G-strands and already some bound ions. Such ensembles can then rearrange into the fully folded parallel GQs via conformational diffusion. We propose that the basic atomistic folding mechanism of propeller loops suggested in this work may be common for their formation in RNA and DNA GQs.

RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview

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.

Insilico direct folding of thrombin-binding aptamer G-quadruplex at all-atom level

The reversible folding of the thrombin-binding DNA aptamer G-quadruplexes (GQs) (TBA-15) starting from fully unfolded states was demonstrated using a prolonged time scale (10-12 μs) parallel tempering metadynamics (PTMetaD) simulation method in conjunction with a modified version of the AMBER bsc1 force field. For unbiased descriptions of the folding free energy landscape of TBA-15, this force field was minimally modified. From this direct folding simulation using the modified bsc1 force field, reasonably converged free energy landscapes were obtained in K+-rich aqueous solution (150 mM), providing detailed atomistic pictures of GQ folding mechanisms for TBA-15. This study found that the TBA folding occurred via multiple folding pathways with two major free energy barriers of 13 and 15 kcal/mol in the presence of several intermediate states of G-triplex variants. The early formation of these intermediates was associated with a single K+ ion capturing. Interestingly, these intermediate states appear to undergo facile transitions among themselves through relatively small energy barriers.

Mesoscopic model for DNA G-quadruplex unfolding

Genomes contain rare guanine-rich sequences capable of assembling into four-stranded helical structures, termed G-quadruplexes, with potential roles in gene regulation and chromosome stability. Their mechanical unfolding has only been reported to date by all-atom simulations, which cannot dissect the major physical interactions responsible for their cohesion. Here, we propose a mesoscopic model to describe both the mechanical and thermal stability of DNA G-quadruplexes, where each nucleotide of the structure, as well as each central cation located at the inner channel, is mapped onto a single bead. In this framework we are able to simulate loading rates similar to the experimental ones, which are not reachable in simulations with atomistic resolution. In this regard, we present single-molecule force-induced unfolding experiments by a high-resolution optical tweezers on a DNA telomeric sequence capable of adopting a G-quadruplex conformation. Fitting the parameters of the model to the experiments we find a correct prediction of the rupture-force kinetics and a good agreement with previous near equilibrium measurements. Since G-quadruplex unfolding dynamics is halfway in complexity between secondary nucleic acids and tertiary protein structures, our model entails a nanoscale paradigm for non-equilibrium processes in the cell.

Multi-scale simulations of biological systems using the OPEP coarse-grained model

Biomolecules are complex machines that are optimized by evolution to properly fulfill or contribute to a variety of biochemical tasks in the cellular environment. Computer simulations based on quantum mechanics and atomistic force fields have been proven to be a powerful microscope for obtaining valuable insights into many biological, physical, and chemical processes. Many interesting phenomena involve, however, a time scale and a number of degrees of freedom, notably if crowding is considered, that cannot be explored at an atomistic resolution. To bridge the gap between reality and simulation, many different advanced computational techniques and coarse-grained (CG) models have been developed. Here, we report some applications of the CG OPEP protein model to amyloid fibril formation, the response of catch-bond proteins to two types of fluid flow, and interactive simulations to fold peptides with well-defined 3D structures or with intrinsic disorder.

Effect of Monovalent Ion Parameters on Molecular Dynamics Simulations of G-Quadruplexes

G-quadruplexes (GQs) are key noncanonical DNA and RNA architectures stabilized by desolvated monovalent cations present in their central channels. We analyze extended atomistic molecular dynamics simulations (∼580 μs in total) of GQs with 11 monovalent cation parametrizations, assessing GQ overall structural stability, dynamics of internal cations, and distortions of the G-tetrad geometries. Majority of simulations were executed with the SPC/E water model; however, test simulations with TIP3P and OPC water models are also reported. The identity and parametrization of ions strongly affect behavior of a tetramolecular d[GGG]₄ GQ, which is unstable with several ion parametrizations. The remaining studied RNA and DNA GQs are structurally stable, though the G-tetrad geometries are always deformed by bifurcated H-bonding in a parametrization-specific manner. Thus, basic 10-μs-scale simulations of fully folded GQs can be safely done with a number of cation parametrizations. However, there are parametrization-specific differences and basic force-field errors affecting the quantitative description of ion-tetrad interactions, which may significantly affect studies of the ion-binding processes and description of the GQ folding landscape. Our d[GGG]₄ simulations indirectly suggest that such studies will also be sensitive to the water models. During exchanges with bulk water, the Na⁺ ions move inside the GQs in a concerted manner, while larger relocations of the K⁺ ions are typically separated. We suggest that the Joung-Cheatham SPC/E K⁺ parameters represent a safe choice in simulation studies of GQs, though variation of ion parameters can be used for specific simulation goals.

What Can Human-Guided Simulations Bring to RNA Folding?

Inspired by the recent success of scientific-discovery games for predicting protein tertiary and RNA secondary structures, we have developed an open software for coarse-grained RNA folding simulations, guided by human intuition. To determine the extent to which interactive simulations can accurately predict 3D RNA structures of increasing complexity and lengths (four RNAs with 22-47 nucleotides), an interactive experiment was conducted with 141 participants who had very little knowledge of nucleic acids systems and computer simulations, and had received only a brief description of the important forces stabilizing RNA structures. Their structures and full trajectories have been analyzed statistically and compared to standard replica exchange molecular dynamics simulations. Our analyses show that participants gain easily chemical intelligence to fold simple and nontrivial topologies, with little computer time, and this result opens the door for the use of human-guided simulations to RNA folding. Our experiment shows that interactive simulations have better chances of success when the user widely explores the conformational space. Interestingly, providing on-the-fly feedback of the root mean square deviation with respect to the experimental structure did not improve the quality of the proposed models.

Martini Coarse-Grained Force Field: Extension to RNA

RNA has an important role not only as the messenger of genetic information but also as a regulator of gene expression. Given its central role in cell biology, there is significant interest in studying the structural and dynamic behavior of RNA in relation to other biomolecules. Coarse-grain molecular dynamics simulations are a key tool to that end. Here, we have extended the coarse-grain Martini force field to include RNA after our recent extension to DNA. In the same way DNA was modeled, the tertiary structure of RNA is constrained using an elastic network. This model, therefore, is not designed for applications involving RNA folding but rather offers a stable RNA structure for studying RNA interactions with other (bio)molecules. The RNA model is compatible with all other Martini models and opens the way to large-scale explicit-solvent molecular dynamics simulations of complex systems involving RNA.

Channeling through Two Stacked Guanine Quartets of One and Two Alkali Cations in the Li+, Na+, K+, and Rb+ Series. Assessment of the Accuracy of the SIBFA Anisotropic Polarizable Molecular Mechanics Potential

Stacking of guanine quartets (GQs) can trigger the formation of DNA or RNA quadruple helices, which play numerous biochemical roles. The GQs are stabilized by alkali cations, mainly K⁺ and Na⁺, which can reside in, or channel through, the central axis of the GQ stems. Further, ion conduction through GQ wires can be leveraged for nanochemistry applications. G-quadruplex systems have been extensively studied by classical molecular dynamics (MD) simulations using pair-additive force fields or by quantum-chemical (QC) calculations. However, the non-polarizable force fields are very approximate, while QC calculations lack the necessary sampling. Thus, ultimate description of GQ systems would require long-enough simulations using advanced polarizable molecular mechanics (MM). However, to perform such calculations, it is first mandatory to evaluate the method's accuracy using benchmark QC. We report such an evaluation for SIBFA polarizable MM, bearing on the channeling (movement) of an alkali cation (Li⁺, Na⁺, K⁺, or Rb⁺) along the axis of two stacked G quartets interacting with either one or two ions. The QC energy profiles display markedly different features depending upon the cation but can be retrieved in the majority of cases by the SIBFA profiles. An appropriate balance of first-order (electrostatic and short-range repulsion) and second-order (polarization, charge-transfer, and dispersion) contributions within ΔE is mandatory. With two cations in the channel, the relative weights of the second-order contributions increase steadily upon increasing the ion size. In the G8 complexes with two K⁺ or two Rb⁺ cations, the sum of polarization and charge-transfer exceeds the first order terms for all ion positions.

Folding of guanine quadruplex molecules-funnel-like mechanism or kinetic partitioning? An overview from MD simulation studies

BACKGROUND

Guanine quadruplexes (GQs) play vital roles in many cellular processes and are of much interest as drug targets. In contrast to the availability of many structural studies, there is still limited knowledge on GQ folding.

SCOPE OF REVIEW

We review recent molecular dynamics (MD) simulation studies of the folding of GQs, with an emphasis paid to the human telomeric DNA GQ. We explain the basic principles and limitations of all types of MD methods used to study unfolding and folding in a way accessible to non-specialists. We discuss the potential role of G-hairpin, G-triplex and alternative GQ intermediates in the folding process. We argue that, in general, folding of GQs is fundamentally different from funneled folding of small fast-folding proteins, and can be best described by a kinetic partitioning (KP) mechanism. KP is a competition between at least two (but often many) well-separated and structurally different conformational ensembles.

MAJOR CONCLUSIONS

The KP mechanism is the only plausible way to explain experiments reporting long time-scales of GQ folding and the existence of long-lived sub-states. A significant part of the natural partitioning of the free energy landscape of GQs comes from the ability of the GQ-forming sequences to populate a large number of syn-anti patterns in their G-tracts. The extreme complexity of the KP of GQs typically prevents an appropriate description of the folding landscape using just a few order parameters or collective variables.

GENERAL SIGNIFICANCE

We reconcile available computational and experimental studies of GQ folding and formulate basic principles characterizing GQ folding landscapes. This article is part of a Special Issue entitled "G-quadruplex" Guest Editor: Dr. Concetta Giancola and Dr. Daniela Montesarchio.