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

Computer Folding of Parallel DNA G-Quadruplex: Hitchhiker's Guide to the Conformational Space

Guanine quadruplexes (GQs) play crucial roles in various biological processes, and understanding their folding pathways provides insight into their stability, dynamics, and functions. This knowledge aids in designing therapeutic strategies, as GQs are potential targets for anticancer drugs and other therapeutics. Although experimental and theoretical techniques have provided valuable insights into different stages of the GQ folding, the structural complexity of GQs poses significant challenges, and our understanding remains incomplete. This study introduces a novel computational protocol for folding an entire GQ from single-strand conformation to its native state. By combining two complementary enhanced sampling techniques, we were able to model folding pathways, encompassing a diverse range of intermediates. Although our investigation of the GQ free energy surface (FES) is focused solely on the folding of the all-anti parallel GQ topology, this protocol has the potential to be adapted for the folding of systems with more complex folding landscapes.

Complex Biophysical and Computational Analyses of G-Quadruplex Ligands: The Porphyrin Stacks Back

G-quadruplexes (G4 s), as non-canonical DNA structures, attract a great deal of research interest in the molecular biology as well as in the material science fields. The use of small molecules as ligands for G-quadruplexes has emerged as a tool to regulate gene expression and telomeres maintenance. Meso-tetrakis-(N-methyl-4-pyridyl) porphyrin (TMPyP4) was shown as one of the first ligands for G-quadruplexes and it is still widely used. We report an investigation comprising molecular docking and dynamics, synthesis and multiple spectroscopic and spectrometric determinations on simple cationic porphyrins and their interaction with different DNA sequences. This study enabled the synthesis of tetracationic porphyrin derivatives that exhibited binding and stabilizing capacity against G-quadruplex structures; the detailed characterization has shown that the presence of amide groups at the periphery improves selectivity for parallel G4 s binding over other structures. Taking into account the ease of synthesis, 5,10,15,20-tetrakis-(1-acetamido-4-pyridyl) porphyrin bromide could be considered a better alternative to TMPyP4 in studies involving G4 binding.

Insights into the Molecular Structure, Stability, and Biological Significance of Non-Canonical DNA Forms, with a Focus on G-Quadruplexes and i-Motifs

This article provides a comprehensive examination of non-canonical DNA structures, particularly focusing on G-quadruplexes (G4s) and i-motifs. G-quadruplexes, four-stranded structures formed by guanine-rich sequences, are stabilized by Hoogsteen hydrogen bonds and monovalent cations like potassium. These structures exhibit diverse topologies and are implicated in critical genomic regions such as telomeres and promoter regions of oncogenes, playing significant roles in gene expression regulation, genome stability, and cellular aging. I-motifs, formed by cytosine-rich sequences under acidic conditions and stabilized by hemiprotonated cytosine-cytosine (C:C+) base pairs, also contribute to gene regulation despite being less prevalent than G4s. This review highlights the factors influencing the stability and dynamics of these structures, including sequence composition, ionic conditions, and environmental pH. Molecular dynamics simulations and high-resolution structural techniques have been pivotal in advancing our understanding of their folding and unfolding mechanisms. Additionally, the article discusses the therapeutic potential of small molecules designed to selectively bind and stabilize G4s and i-motifs, with promising implications for cancer treatment. Furthermore, the structural properties of these DNA forms are explored for applications in nanotechnology and molecular devices. Despite significant progress, challenges remain in observing these structures in vivo and fully elucidating their biological functions. The review underscores the importance of continued research to uncover new insights into the genomic roles of G4s and i-motifs and their potential applications in medicine and technology. This ongoing research promises exciting developments in both basic science and applied fields, emphasizing the relevance and future prospects of these intriguing DNA structures.

Interlocking G-Quadruplexes Using a G-Triad•G Connection: Implications for G-Wire Assembly

G-quadruplexes are noncanonical structures of nucleic acids formed mainly by G-rich sequences and play crucial roles in important cellular processes. They are also increasingly used in nanotechnology for their valuable properties. Various unexpected structures of G-quadruplexes have been solved recently, including a stable G-quadruplex lacking one guanine in the G-tetrad core, harboring a vacant site. In this study, we demonstrate the interlocking of two intramolecular G-quadruplexes: one containing a vacant site (4n - 1) and the other with an unbound guanine (4n + 1). These G-quadruplexes interact through a G-triad-G connection with unprecedented 5'-3' stacking. Using these interconnection properties, we have identified a sequence capable of self-assembling into G-wires in K+ solutions with potential nanotechnological applications.

Early Events in G-quadruplex Folding Captured by Time-Resolved Small-Angle X-Ray Scattering

Time-resolved small-angle X-ray experiments (TR-SAXS) are reported here that capture and quantify a previously unknown rapid collapse of the unfolded oligonucleotide as an early step in G4 folding of hybrid 1 and hybrid 2 telomeric G-quadruplex structures. The rapid collapse, initiated by a pH jump, is characterized by an exponential decrease in the radius of gyration from 20.6 to 12.6 Å. The collapse is monophasic and is complete in less than 600 ms. Additional hand-mixing pH-jump kinetic studies show that slower kinetic steps follow the collapse. The folded and unfolded states at equilibrium were further characterized by SAXS studies and other biophysical tools, to show that G4 unfolding was complete at alkaline pH, but not in LiCl solution as is often claimed. The SAXS Ensemble Optimization Method (EOM) analysis reveals models of the unfolded state as a dynamic ensemble of flexible oligonucleotide chains with a variety of transient hairpin structures. These results suggest a G4 folding pathway in which a rapid collapse, analogous to molten globule formation seen in proteins, is followed by a confined conformational search within the collapsed particle to form the native contacts ultimately found in the stable folded form.

Structural Differences at Quadruplex-Duplex Interfaces Enable Ligand-Induced Topological Transitions

Quadruplex-duplex (QD) junctions, which represent unique structural motifs of both biological and technological significance, have been shown to constitute high-affinity binding sites for various ligands. A QD hybrid construct based on a human telomeric sequence, which harbors a duplex stem-loop in place of a short lateral loop, is structurally characterized by NMR. It folds into two major species with a (3+1) hybrid and a chair-type (2+2) antiparallel quadruplex domain coexisting in a K+ buffer solution. The antiparallel species is stabilized by an unusual capping structure involving a thymine and protonated adenine base AH+ of the lateral loop facing the hairpin duplex to form a T·AH+·G·C quartet with the interfacial G·C base pair at neutral pH. Addition and binding of Phen-DC3 to the QD hybrid mixture by its partial intercalation at corresponding QD junctions leads to a topological transition with exclusive formation of the (3+1) hybrid fold. In agreement with the available experimental data, such an unprecedented discrimination of QD junctions by a ligand can be rationalized following an induced fit mechanism.

Mechanical Stability and Unfolding Pathways of Parallel Tetrameric G-Quadruplexes Probed by Pulling Simulations

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.

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.

A General Framework to Interpret Hydrogen-Deuterium Exchange Native Mass Spectrometry of G-Quadruplex DNA

G-quadruplexes (G4s) are secondary structures formed by guanine-rich oligonucleotides involved in various biological processes. However, characterizing G4s is challenging, because of their structural polymorphism. Here, we establish how hydrogen-deuterium exchange native mass spectrometry (HDX/MS) can help to characterize the G4 structures and dynamics in solution. We correlated the time range of G4 exchange to the number of guanines involved in the inner and outer tetrads. We also established relationships among exchange rates, numbers of tetrads and bound cations, and stability. The use of HDX/native MS allows for the determination of tetrads formed and assessment of G4 stability at a constant temperature. A key finding is that stable G4s exchange through local fluctuations (EX2 exchange), whereas less stable G4s also undergo exchange through partial or complete unfolding (EX1 exchange). Deconvolution of the bimodal isotope distributions resulting from EX1 exchange provides valuable insight into the kinetics of folding and unfolding processes and allows one to detect and characterize transiently unfolded intermediates, even if scarcely populated. HDX/native MS thus represents a powerful tool for a more comprehensive exploration of the folding landscapes of G4s.

Computational Probing of the Folding Mechanism of Human Telomeric G-Quadruplex DNA

The human telomeric (htel) sequences in the terminal regions of human telomeres form diverse G-quadruplex (GQ) structures. Despite much experimental efforts to elucidate the folding pathways of htel GQ, no comprehensive model of htel GQ folding has been presented. Here, we describe folding pathways of the htel GQ determined by state-of-the-art enhanced sampling molecular dynamics simulation at the all-atom level. Briefly, GQ folding is initiated by the formation of a single-hairpin and then followed by the formation of double-hairpins, which then branch via distinct folding pathways to produce different GQ topologies (antiparallel chair, antiparallel basket, hybrids 1 and 2, and parallel propeller). In addition to these double-hairpin states, three-triad and two-tetrad structures in antiparallel backbone alignment serve as key intermediates that connect the GQ folding and transition between two different GQs.

Complexity of Guanine Quadruplex Unfolding Pathways Revealed by Atomistic Pulling Simulations

Guanine quadruplexes (GQs) are non-canonical nucleic acid structures involved in many biological processes. GQs formed in single-stranded regions often need to be unwound by cellular machinery, so their mechanochemical properties are important. Here, we performed steered molecular dynamics simulations of human telomeric GQs to study their unfolding. We examined four pulling regimes, including a very slow setup with pulling velocity and force load accessible to high-speed atomic force microscopy. We identified multiple factors affecting the unfolding mechanism, i.e.,: (i) the more the direction of force was perpendicular to the GQ channel axis (determined by GQ topology), the more the base unzipping mechanism happened, (ii) the more parallel the direction of force was, GQ opening and cross-like GQs were more likely to occur, (iii) strand slippage mechanism was possible for GQs with an all-anti pattern in a strand, and (iv) slower pulling velocity led to richer structural dynamics with sampling of more intermediates and partial refolding events. We also identified that a GQ may eventually unfold after a force drop under forces smaller than those that the GQ withstood before the drop. Finally, we found out that different unfolding intermediates could have very similar chain end-to-end distances, which reveals some limitations of structural interpretations of single-molecule spectroscopic data.

Inhibited complete folding of consecutive human telomeric G-quadruplexes

Noncanonical DNA structures, termed G-quadruplexes, are present in human genomic DNA and are important elements in many DNA metabolic processes. Multiple sites in the human genome have G-rich DNA stretches able to support formation of several consecutive G-quadruplexes. One of those sites is the telomeric overhang region that has multiple repeats of TTAGGG and is tightly associated with both cancer and aging. We investigated the folding of consecutive G-quadruplexes in both potassium- and sodium-containing solutions using single-molecule FRET spectroscopy, circular dichroism, thermal melting and molecular dynamics simulations. Our observations show coexistence of partially and fully folded DNA, the latter consisting of consecutive G-quadruplexes. Following the folding process over hours in sodium-containing buffers revealed fast G-quadruplex folding but slow establishment of thermodynamic equilibrium. We find that full consecutive G-quadruplex formation is inhibited by the many DNA structures randomly nucleating on the DNA, some of which are off-path conformations that need to unfold to allow full folding. Our study allows describing consecutive G-quadruplex formation in both nonequilibrium and equilibrium conditions by a unified picture, where, due to the many possible DNA conformations, full folding with consecutive G-quadruplexes as beads on a string is not necessarily achieved.

Phen-DC3 Induces Refolding of Human Telomeric DNA into a Chair-Type Antiparallel G-Quadruplex through Ligand Intercalation

Human telomeric G-quadruplex DNA structures are attractive anticancer drug targets, but the target's polymorphism complicates the drug design: different ligands prefer different folds, and very few complexes have been solved at high resolution. Here we report that Phen-DC3 , one of the most prominent G-quadruplex ligands in terms of high binding affinity and selectivity, causes dTAGGG(TTAGGG)3 to completely change its fold in KCl solution from a hybrid-1 to an antiparallel chair-type structure, wherein the ligand intercalates between a two-quartet unit and a pseudo-quartet, thereby ejecting one potassium ion. This unprecedented high-resolution NMR structure shows for the first time a true ligand intercalation into an intramolecular G-quadruplex.

Tracking Topological and Electronic Effects on the Folding and Stability of Guanine-Deficient RNA G-Quadruplexes, Engineered with a New Computational Tool for De Novo Quadruplex Folding

The initial aim of this work was to elucidate the mutual influence of different single-stranded segments (loops and caps) on the thermodynamic stability of RNA G-quadruplexes. To this end, we used a new NAB-GQ-builder software program, to construct dozens of two-tetrad G-quadruplex topologies, based on a designed library of sequences. Then, to probe the sequence-morphology-stability relationships of the designed topologies, we performed molecular dynamics simulations. Their results provide guidance for the design of G-quadruplexes with balanced structures, and in turn programmable physicochemical properties for applications as biomaterials. Moreover, by comparative examinations of the single-stranded segments of three oncogene promoter G-quadruplexes, we assess their druggability potential for future therapeutic strategies. Finally, on the basis of a thorough analysis at the quantum mechanical level of theory on a series of guanine assemblies, we demonstrate how a valence tautomerism, triggered by a coordination of cations, initiates the process of G-quadruplex folding, and we propose a sequential folding mechanism, otherwise dictated by the cancellation of the dipole moments on guanines.

Stability and Existence of Noncanonical I-motif DNA Structures in Computer Simulations Based on Atomistic and Coarse-Grained Force Fields

Cytosine-rich DNA sequences are able to fold into noncanonical structures, in which semi-protonated cytosine pairs develop extra hydrogen bonds, and these bonds are responsible for the overall stability of a structure called the i-motif. The i-motif can be formed in many regions of the genome, but the most representative is the telomeric region in which the CCCTAA sequences are repeated thousands of times. The ability to reverse folding/unfolding in response to pH change makes the above sequence and i-motif very promising components of nanomachines, extended DNA structures, and drug carriers. Molecular dynamics analysis of such structures is highly beneficial due to direct insights into the microscopic structure of the considered systems. We show that Amber force fields for DNA predict the stability of the i-motif over a long timescale; however, these force fields are not able to predict folding of the cytosine-rich sequences into the i-motif. The reason is the kinetic partitioning of the folding process, which makes the transitions between various intermediates too time-consuming in atomistic force field representation. Application of coarse-grained force fields usually highly accelerates complex structural transitions. We, however, found that three of the most popular coarse-grained force fields for DNA (oxDNA, 3SPN, and Martini) were not able to predict the stability of the i-motif structure. Obviously, they were not able to accelerate the folding of unfolded states into an i-motif. This observation must be strongly highlighted, and the need to develop suitable extensions of coarse-grained force fields for DNA is pointed out. However, it will take a great deal of effort to successfully solve these problems.

Emerging accessibility patterns in long telomeric overhangs

We present single-molecule experimental and computational modeling studies investigating the accessibility of human telomeric overhangs of physiologically relevant lengths. We studied 25 different overhangs that contain 4-28 repeats of GGGTTA (G-Tract) sequence and accommodate one to seven tandem G-quadruplex (GQ) structures. Using the FRET-PAINT method, we probed the distribution of accessible sites via a short imager strand, which is complementary to a G-Tract and transiently binds to available sites. We report accessibility patterns that periodically change with overhang length and interpret these patterns in terms of the underlying folding landscape and folding frustration. Overhangs that have [4n]G-Tracts, (12, 16, 20…) demonstrate the broadest accessibility patterns where the peptide nucleic acid probe accesses G-Tracts throughout the overhang. On the other hand, constructs with [4n+2]G-Tracts, (14, 18, 22…) have narrower patterns where the neighborhood of the junction between single- and double-stranded telomeres is most accessible. We interpret these results as the folding frustration being higher in [4n]G-Tract constructs compared to [4n+2]G-Tract constructs. We also developed a computational model that tests the consistency of different folding stabilities and cooperativities between neighboring GQs with the observed accessibility patterns. Our experimental and computational studies suggest the neighborhood of the junction between single- and double-stranded telomeres is least stable and most accessible, which is significant as this is a potential site where the connection between POT1/TPP1 (bound to single-stranded telomere) and other shelterin proteins (localized on double-stranded telomere) is established.

Ion-Dependent Conformational Plasticity of Telomeric G-Hairpins and G-Quadruplexes

Telomeric DNA is guanine-rich and can adopt structures such as G-quadruplexes (GQs) and G-hairpins. Telomeric GQs influence genome stability and telomerase activity, making understanding of enzyme-GQ interactions and dynamics important for potential drug design. GQs have a characteristic tetrad core, which is connected by loop regions. Within this architecture are G-hairpins, fold-back motifs that are thought to represent the first intermediate in GQ folding. To better understand the relationship between G-hairpin motifs and GQs, we performed polarizable simulations of a two-tetrad telomeric GQ and an isolated SC11 telomeric G-hairpin. The telomeric GQ contains a G-triad, which functions as part of the tetrad core or linker regions, depending on local conformational change. This triad and another motif below the tetrad core frequently bound ions and may represent druggable sites. Further, we observed the unbiased formation of a G-triad and a G-tetrad in simulations of the SC11 G-hairpin and found that cations can be partially hydrated while facilitating the formation of these motifs. Finally, we demonstrated that K⁺ ions form specific interactions with guanine bases, while Na⁺ ions interact nonspecifically with bases in the structure. Together, these simulations provide new insights into the influence of ions on GQs, G-hairpins, and G-triad motifs.

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.

Insights into the free energy landscape and salt-controlled mechanism of the conformational conversions between human telomeric G-quadruplex structures

G-quadruplexes have become attractive drug targets in cancer therapy. However, due to the polymorphism of G-quadruplex structures, it is difficult to experimentally verify the relevant structures of multiple intermediates and transition states in dynamic equilibrium. Hence, understanding the mechanism by which structural conversions of G-quadruplexes occur is still challenging. We conducted targeted molecular dynamics simulation with umbrella sampling to investigate how salt affects the conformational conversion of human telomeric G-quadruplex. Our results explore a unique view into the structures and energy barrier of the intermediates and transition states in the interconversion process. The pathway of G-quadruplex conformational interconversion was mapped out by a free energy landscape, consisting of branched parallel pathways with multiple energy basins. We propose a salt-controlled mechanism that as the salt concentration increases, the conformational conversion mechanism switches from multi-pathway folding to sequential folding pathways. The hybrid-I and hybrid-II structures are intermediates in the basket-propeller transformation. In high-salt solutions, the conformational conversion upon K⁺ binding is more feasible than upon Na⁺ binding. The free energy barrier for conformational conversions ranges from 1.6 to 4.6 kcal/mol. Our work will be beneficial in developing anticancer agents.

Folding dynamics of polymorphic G-quadruplex structures

G-quadruplexes (G4), found in numerous places within the human genome, are involved in essential processes of cell regulation. Chromosomal DNA G4s are involved for example, in replication and transcription as first steps of gene expression. Hence, they influence a plethora of downstream processes. G4s possess an intricate structure that differs from canonical B-form DNA. Identical DNA G4 sequences can adopt multiple long-lived conformations, a phenomenon known as G4 polymorphism. A detailed understanding of the molecular mechanisms that drive G4 folding is essential to understand their ambivalent regulatory roles. Disentangling the inherent dynamic and polymorphic nature of G4 structures thus is key to unravel their biological functions and make them amenable as molecular targets in novel therapeutic approaches. We here review recent experimental approaches to monitor G4 folding and discuss structural aspects for possible folding pathways. Substantial progress in the understanding of G4 folding within the recent years now allows drawing comprehensive models of the complex folding energy landscape of G4s that we herein evaluate based on computational and experimental evidence.

Thermodynamic Stability of G-Quadruplexes: Impact of Sequence and Environment

G-quadruplexes have attracted growing interest in recent years due to their occurrence in vivo and their possible biological functions. In addition to being promising targets for drug design, these four-stranded nucleic acid structures have also been recognized as versatile tools for various technological applications. Whereas a large number of studies have yielded insight into their remarkable structural diversity, our current knowledge on G-quadruplex stabilities as a function of sequence and environmental factors only gradually emerges with an expanding collection of thermodynamic data. This minireview provides an overview of general rules that may be used to better evaluate quadruplex thermodynamic stabilities but also discusses present challenges in predicting most stable folds for a given sequence and environment.

Mass Spectrometry of Nucleic Acid Noncovalent Complexes

Nucleic acids have been among the first targets for antitumor drugs and antibiotics. With the unveiling of new biological roles in regulation of gene expression, specific DNA and RNA structures have become very attractive targets, especially when the corresponding proteins are undruggable. Biophysical assays to assess target structure as well as ligand binding stoichiometry, affinity, specificity, and binding modes are part of the drug development process. Mass spectrometry offers unique advantages as a biophysical method owing to its ability to distinguish each stoichiometry present in a mixture. In addition, advanced mass spectrometry approaches (reactive probing, fragmentation techniques, ion mobility spectrometry, ion spectroscopy) provide more detailed information on the complexes. Here, we review the fundamentals of mass spectrometry and all its particularities when studying noncovalent nucleic acid structures, and then review what has been learned thanks to mass spectrometry on nucleic acid structures, self-assemblies (e.g., duplexes or G-quadruplexes), and their complexes with ligands.

New perspectives in cancer drug development: computational advances with an eye to design

Computational chemistry has come of age in drug discovery. Indeed, most pharmaceutical development programs rely on computer-based data and results at some point. Herein, we discuss recent applications of advanced simulation techniques to difficult challenges in drug discovery. These entail the characterization of allosteric mechanisms and the identification of allosteric sites or cryptic pockets determined by protein motions, which are not immediately evident in the experimental structure of the target; the study of ligand binding mechanisms and their kinetic profiles; and the evaluation of drug-target affinities. We analyze different approaches to tackle challenging and emerging biological targets. Finally, we discuss the possible perspectives of future application of computation in drug discovery.

G-Quadruplex Formation by DNA Sequences Deficient in Guanines: Two Tetrad Parallel Quadruplexes Do Not Fold Intramolecularly

Guanine quadruplexes (G4s) are noncanonical forms of nucleic acids that are frequently found in genomes. The stability of G4s depends, among other factors, on the number of G-tetrads. Three- or four-tetrad G4s and antiparallel two-tetrad G4s have been characterized experimentally; however, the existence of an intramolecular (i. e., not dimeric or multimeric) two-tetrad parallel-stranded DNA G4 has never been experimentally observed. Many sequences compatible with two-tetrad G4 can be found in important genomic regions, such as promoters, for which parallel G4s predominate. Using experimental and theoretical approaches, the propensity of the model sequence AATGGGTGGGTTTGGGTGGGTAA to form an intramolecular parallel-stranded G4 upon increasing the number of GGG-to-GG substitutions has been studied. Deletion of a single G leads to the formation of intramolecular G4s with a stacked G-triad, whose topology depends on the location of the deletion. Removal of another guanine from another G-tract leads to di- or multimeric G4s. Further deletions mostly prevent the formation of any stable G4. Thus, a solitary two-tetrad parallel DNA G4 is not thermodynamically stable and requires additional interactions through capping residues. However, transiently populated metastable two-tetrad species can associate to form stable dimers, the dynamic formation of which might play additional delicate roles in gene regulation. These findings provide essential information for bioinformatics studies searching for potential G4s in genomes.

Mechanical diversity and folding intermediates of parallel-stranded G-quadruplexes with a bulge

A significant number of sequences in the human genome form noncanonical G-quadruplexes (G4s) with bulges or a guanine vacancy. Here, we systematically characterized the mechanical stability of parallel-stranded G4s with a one to seven nucleotides bulge at various positions. Our results show that G4-forming sequences with a bulge form multiple conformations, including fully-folded G4 with high mechanical stability (unfolding forces > 40 pN), partially-folded intermediates (unfolding forces < 40 pN). The folding probability and folded populations strongly depend on the positions and lengths of the bulge. By combining a single-molecule unfolding assay, dimethyl sulfate (DMS) footprinting, and a guanine-peptide conjugate that selectively stabilizes guanine-vacancy-bearing G-quadruplexes (GVBQs), we identified that GVBQs are the major intermediates of G4s with a bulge near the 5′ or 3′ ends. The existence of multiple structures may induce different regulatory functions in many biological processes. This study also demonstrates a new strategy for selectively stabilizing the intermediates of bulged G4s to modulate their functions.

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.

Spectroscopic and In Silico Studies on the Interaction of Substituted Pyrazolo[1,2-a]benzo[1,2,3,4]tetrazine-3-one Derivatives with c-Myc G4-DNA

Herein we describe a combined experimental and in silico study of the interaction of a series of pyrazolo[1,2-a]benzo[1,2,3,4]tetrazin-3-one derivatives (PBTs) with parallel G-quadruplex (GQ) DNA aimed at correlating their previously reported anticancer activities and the stabilizing effects observed by us on c-myc oncogene promoter GQ structure. Circular dichroism (CD) melting experiments were performed to characterize the effect of the studied PBTs on the GQ thermal stability. CD measurements indicate that two out of the eight compounds under investigation induced a slight stabilizing effect (2-4 °C) on GQ depending on the nature and position of the substituents. Molecular docking results allowed us to verify the modes of interaction of the ligands with the GQ and estimate the binding affinities. The highest binding affinity was observed for ligands with the experimental melting temperatures (Tms). However, both stabilizing and destabilizing ligands showed similar scores, whilst Molecular Dynamics (MD) simulations, performed across a wide range of temperatures on the GQ in water solution, either unliganded or complexed with two model PBT ligands with the opposite effect on the Tms, consistently confirmed their stabilizing or destabilizing ability ascertained by CD. Clues about a relation between the reported anticancer activity of some PBTs and their ability to stabilize the GQ structure of c-myc emerged from our study. Furthermore, Molecular Dynamics simulations at high temperatures are herein proposed for the first time as a means to verify the stabilizing or destabilizing effect of ligands on the GQ, also disclosing predictive potential in GQ-targeting drug discovery.

The Folding Landscapes of Human Telomeric RNA and DNA G-Quadruplexes are Markedly Different

We investigated the folding kinetics of G-quadruplex (G4) structures by comparing the K+ -induced folding of an RNA G4 derived from the human telomeric repeat-containing RNA (TERRA25) with a sequence homologous DNA G4 (wtTel25) using CD spectroscopy and real-time NMR spectroscopy. While DNA G4 folding is biphasic, reveals kinetic partitioning and involves kinetically favoured off-pathway intermediates, RNA G4 folding is faster and monophasic. The differences in kinetics are correlated to the differences in the folded conformations of RNA vs. DNA G4s, in particular with regard to the conformation around the glycosidic torsion angle χ that uniformly adopts anti conformations for RNA G4s and both, syn and anti conformation for DNA G4s. Modified DNA G4s with 19 F bound to C2' in arabino configuration adopt exclusively anti conformations for χ. These fluoro-modified DNA (antiTel25) reveal faster folding kinetics and monomorphic conformations similar to RNA G4s, suggesting the correlation between folding kinetics and pathways with differences in χ angle preferences in DNA and RNA, respectively.

Unraveling the Kinetics of Spare-Tire DNA G-Quadruplex Folding

The folding of DNA G-quadruplexes (G4) is essential to regulate expression of oncogenes and involves polymorphic long-lived intermediate states. G4 formation requires four G-tracts, but human gene-promoters often contain multiple G-tracts that act as spare-tires. These additional G-tracts are highly conserved and add multiple layers of functional complexity, as they are crucial to maintain G4 function after oxidative damage. Herein, we unravel the folding dynamics of the G4 sequence containing five G-tracts from cMYC, the major proliferation-driving oncogene. We devise a general method to induce folding at constant experimental conditions using a photochemical trapping strategy. Our data dissect the individual kinetics and thermodynamics of the spare-tire mechanism of cMYC-G4.

Insight into formation propensity of pseudocircular DNA G-hairpins

We recently showed that Saccharomyces cerevisiae telomeric DNA can fold into an unprecedented pseudocircular G-hairpin (PGH) structure. However, the formation of PGHs in the context of extended sequences, which is a prerequisite for their function in vivo and their applications in biotechnology, has not been elucidated. Here, we show that despite its ‘circular’ nature, PGHs tolerate single-stranded (ss) protrusions. High-resolution NMR structure of a novel member of PGH family reveals the atomistic details on a junction between ssDNA and PGH unit. Identification of new sequences capable of folding into one of the two forms of PGH helped in defining minimal sequence requirements for their formation. Our time-resolved NMR data indicate a possibility that PGHs fold via a complex kinetic partitioning mechanism and suggests the existence of K⁺ ion-dependent PGH folding intermediates. The data not only provide an explanation of cation-type-dependent formation of PGHs, but also explain the unusually large hysteresis between PGH melting and annealing noted in our previous study. Our findings have important implications for DNA biology and nanotechnology. Overrepresentation of sequences able to form PGHs in the evolutionary-conserved regions of the human genome implies their functionally important biological role(s).

Ligand Binding to Dynamically Populated G-Quadruplex DNA

Several small-molecule ligands specifically bind and stabilize G-quadruplex (G4) nucleic acid structures, which are considered to be promising therapeutic targets. G4s are polymorphic structures of varying stability, and their formation is dynamic. Here, we investigate the mechanisms of ligand binding to dynamically populated human telomere G4 DNA by using the bisquinolinium based ligand Phen-DC3 and a combination of single-molecule FRET microscopy, ensemble FRET and CD spectroscopies. Different cations are used to tune G4 polymorphism and folding dynamics. We find that ligand binding occurs to pre-folded G4 structures and that Phen-DC3 also induces G4 formation in unfolded single strands. Following ligand binding to dynamically populated G4s, the DNA undergoes pronounced conformational redistributions that do not involve direct ligand-induced G4 conformational interconversion. On the contrary, the redistribution is driven by ligand-induced G4 folding and trapping of dynamically populated short-lived conformation states. Thus, ligand-induced stabilization does not necessarily require the initial presence of stably folded G4s.

Insights into G-Quadruplex-Hemin Dynamics Using Atomistic Simulations: Implications for Reactivity and Folding

Guanine quadruplex nucleic acids (G4s) are involved in key biological processes such as replication or transcription. Beyond their biological relevance, G4s find applications as biotechnological tools since they readily bind hemin and enhance its peroxidase activity, creating a G4-DNAzyme. The biocatalytic properties of G4-DNAzymes have been thoroughly studied and used for biosensing purposes. Despite hundreds of applications and massive experimental efforts, the atomistic details of the reaction mechanism remain unclear. To help select between the different hypotheses currently under investigation, we use extended explicit-solvent molecular dynamics (MD) simulations to scrutinize the G4/hemin interaction. We find that besides the dominant conformation in which hemin is stacked atop the external G-quartets, hemin can also transiently bind to the loops and be brought to the external G-quartets through diverse delivery mechanisms. The simulations do not support the catalytic mechanism relying on a wobbling guanine. Similarly, the catalytic role of the iron-bound water molecule is not in line with our results; however, given the simulation limitations, this observation should be considered with some caution. The simulations rather suggest tentative mechanisms in which the external G-quartet itself could be responsible for the unique H₂O₂-promoted biocatalytic properties of the G4/hemin complexes. Once stacked atop a terminal G-quartet, hemin rotates about its vertical axis while readily sampling shifted geometries where the iron transiently contacts oxygen atoms of the adjacent G-quartet. This dynamics is not apparent from the ensemble-averaged structure. We also visualize transient interactions between the stacked hemin and the G4 loops. Finally, we investigated interactions between hemin and on-pathway folding intermediates of the parallel-stranded G4 fold. The simulations suggest that hemin drives the folding of parallel-stranded G4s from slip-stranded intermediates, acting as a G4 chaperone. Limitations of the MD technique are briefly discussed.

Ion Binding Properties and Dynamics of the bcl-2 G-Quadruplex Using a Polarizable Force Field

G-quadruplexes (GQs) are topologically diverse, highly thermostable noncanonical nucleic acid structures that form in guanine-rich sequences in DNA and RNA. GQs are implicated in transcriptional and translational regulation and genome maintenance, and deleterious alterations to their structures contribute to diseases such as cancer. The expression of the B-cell lymphoma 2 (Bcl-2) antiapoptotic protein, for example, is under transcriptional control of a GQ in the promoter of the bcl-2 gene. Modulation of the bcl-2 GQ by small molecules is of interest for chemotherapeutic development but doing so requires knowledge of the factors driving GQ folding and stabilization. To develop a greater understanding of the electrostatic properties of the bcl-2 promoter GQ, we performed molecular dynamics simulations using the Drude-2017 polarizable force field and compared relevant outcomes to the nonpolarizable CHARMM36 force field. Our simulation outcomes highlight the importance of dipole-dipole interactions in the bcl-2 GQ, particularly during the recruitment of a bulk K⁺ ion to the solvent-exposed face of the tetrad stem. We also predict and characterize an "electronegative pocket" at the tetrad-long loop junction that induces local backbone conformational change and may induce local conformational changes at cellular concentrations of K⁺. These outcomes suggest that moieties within the bcl-2 GQ can be targeted by small molecules to modulate bcl-2 GQ stability.

Dynamics and Barrier of Movements of Sodium and Potassium Ions Across the Oxytricha nova G-Quadruplex Core

G-quadruplexes (GQs) are highly stable noncanonical forms of nucleic acids that are present in important genomic regions. The central core of the GQ is lined up by four closely spaced carbonyl groups from the G-quartets, and the resulting electrostatic repulsion is neutralized by the coordinating cations. In spite of several reports on GQ structure and cation-GQ interactions, the atomic- to molecular-level understanding of the ion dynamics and ion exchange in the GQ core is quite poor. Here, we attempt to elucidate the mechanism of Na⁺ and K⁺ binding to the GQ core and trace the exchange of these ions with the ions in bulk by means of all-atomic molecular dynamics (MD) simulations. One of the most studied GQs, Oxytricha nova telomeric G-quadruplex (OxyGQ), is taken as the representative GQ. Subsequently, umbrella sampling MD simulations were performed to elucidate the energetics of ion translocation from one end to the other end of the GQ central core. Our study highlights the importance of ion hydration for the uptake and correct positioning of the cations in the core. The free-energy landscape of ion transport has shown favorable in-plane binding of Na⁺ ions with GQ quartets, which matches very well with the crystal structure. The binding of K⁺ ions, on the other hand, was out-of-plane and its translocation required a larger barrier to cross.

Targeting Telomeres: Molecular Dynamics and Free Energy Simulation of Gold-Carbene Binding to DNA

DNA sequences in regulatory regions and in telomers at the ends of chromosomes frequently contain tandem repeats of guanine nucleotides that can form stacked structures stabilized by Hoogsten pairing and centrally bound monovalent cations. The replication and elongation of telomeres requires the disruption of these G-quadruplex structures. Hence, drug molecules such as gold (Au)-carbene that stabilize G-quadruplexes may also interfere with the elongation of telomeres and, in turn, could be used to control cell replication and growth. To better understand the molecular mechanism of Au-carbene binding to G-quadruplexes, we employed molecular dynamics simulations and free energy simulations. Whereas very restricted mobility of two Au-carbene ligands was found upon binding as a doublet to one side of the G-quadruplex, much larger translational and orientational mobility was observed for a single Au-carbene binding at the second G-quadruplex surface. Comparative simulations on duplex DNA in the presence of Au-carbene ligands indicates a preference for the minor groove and weaker unspecific and more salt-dependent binding than to the G-quadruplex surface. Analysis of energetic contributions reveals a dominance of nonpolar and van der Waals interactions to drive binding. The simulations can also be helpful for proposing possible modifications that could improve Au-carbene affinity and specificity for G-quadruplex binding.

Recent developments in empirical atomistic force fields for nucleic acids and applications to studies of folding and dynamics

Nucleic acids play critical roles in carrying genetic information, participating in catalysis, and preserving chromosomal structure. Despite over a century of study, efforts to understand the dynamics and structure-function relationships of DNA and RNA at the atomic level are still ongoing. Molecular dynamics (MD) simulations augment experiments by providing atomistic resolution and quantitative relationships between structure and conformational energy. Steady advancements in computer hardware, software, and atomistic force fields (FFs) over 40 years have facilitated new discoveries. Here, we review nucleic acid FF development with emphasis on recent refinements that have improved descriptions of important nucleic acid properties. We then discuss several key examples of successes and challenges in modeling nucleic acid structure and dynamics using the latest FFs.

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.

Diversity of Parallel Guanine Quadruplexes Induced by Guanine Substitutions

Recently, we reported an inhibitory effect of guanine substitutions on the conformational switch from antiparallel to parallel quadruplexes (G4) induced by dehydrating agents. As a possible cause, we proposed a difference in the sensitivity of parallel and antiparallel quadruplexes to the guanine substitutions in the resulting thermodynamic stability. Reports on the influence of guanine substitutions on the biophysical properties of intramolecular parallel quadruplexes are rare. Moreover, such reports are often complicated by the multimerisation tendencies of parallel quadruplexes. To address this incomplete knowledge, we employed circular dichroism spectroscopy (CD), both as stopped-flow-assisted fast kinetics measurements and end-point measurements, accompanied by thermodynamic analyses, based on UV absorption melting profiles, and electrophoretic methods. We showed that parallel quadruplexes are significantly more sensitive towards guanine substitutions than antiparallel ones. Furthermore, guanine-substituted variants, which in principle might correspond to native genomic sequences, distinctly differ in their biophysical properties, indicating that the four guanines in each tetrad of parallel quadruplexes are not equal. In addition, we were able to distinguish by CD an intramolecular G4 from intermolecular ones resulting from multimerisation mediated by terminal tetrad association, but not from intermolecular G4s formed due to inter-strand Hoogsteen hydrogen bond formation. In conclusion, our study indicates significant variability in parallel quadruplex structures, otherwise disregarded without detailed experimental analysis.

Folding topology, structural polymorphism, and dimerization of intramolecular DNA G-quadruplexes with inverted polarity strands and non-natural loops

Synthetically modified DNA G-quadruplexes (GQs) have great potential in the development of designer molecules for a wide range of applications. Identification of the role of various structural elements in the folding and final topology of artificial GQs is necessary to predict their secondary structure. We report here the results of spectroscopic and electrophoretic studies of GQ scaffolds formed by G-rich sequences comprising four G3-tracts of different polarity connected by either a single-nucleotide thymine loop or a non-natural tetraethyleneglycol loop. Depending on G-strand polarities, loop arrangement and the presence of extra 5'-base G-rich oligonucleotides form monomeric, dimeric, or multimeric species of different topologies. In most cases, oligonucleotides were able to fold into stable parallel or hybrid GQs. However, certain specific arrangements of loops and G-tracts resulted in a diverse mixture of low stable structures. Comparative analysis of topology, stability, and structural heterogeneity of different G-rich sequences suggests the important role of loop type and arrangement, G3-tract polarities, and the presence of 5'-capping residues in the outcome of the folding process. The results also imply that the formation of anti-parallel G-hairpin intermediates is a key event in major favourable folding pathways.

Affinity-Selected Bicyclic Peptide G-Quadruplex Ligands Mimic a Protein-like Binding Mechanism

The study of G-quadruplexes (G4s) in a cellular context has demonstrated links between these nucleic acid secondary structures, gene expression, and DNA replication. Ligands that bind to the G4 structure therefore present an excellent opportunity for influencing gene expression through the targeting of a nucleic acid structure rather than sequence. Here, we explore cyclic peptides as an alternative class of G4 ligands. Specifically, we describe the development of de novo G4-binding bicyclic peptides selected by phage display. Selected bicyclic peptides display submicromolar affinity to G4 structures and high selectivity over double helix DNA. Molecular simulations of the bicyclic peptide-G4 complexes corroborate the experimental binding strengths and reveal molecular insights into G4 recognition by bicyclic peptides via the precise positioning of amino acid side chains, a binding mechanism reminiscent of endogenous G4-binding proteins. Overall, our results demonstrate that selection of (bi)cyclic peptides unlocks a valuable chemical space for targeting nucleic acid structures.

Stability of Two-Quartet G-Quadruplexes and Their Dimers in Atomistic Simulations

G-quadruplexes (GQs) are four-stranded noncanonical DNA and RNA architectures that can be formed by guanine-rich sequences. The stability of GQs increases with the number of G-quartets, and three G-quartets generally form stable GQs. However, the stability of two-quartet GQs is an open issue. To understand the intrinsic stability of two-quartet GQ stems, we have carried out a series of unbiased molecular dynamics (MD) simulations (505 μs in total) of two- and four-quartet DNA and RNA GQs, with attention paid mainly to parallel-stranded arrangements. We used AMBER DNA parmOL15 and RNA parmOL3 force fields and tested different ion and water models. Two-quartet parallel-stranded DNA GQs unfolded in all the simulations, while the equivalent RNA GQ was stable in most of the simulations. GQs composed of two stacked units of two-quartet GQs were stable for both DNA and RNA. The simulations suggest that a minimum of three quartets are needed to form an intrinsically stable all-anti parallel-stranded DNA GQ. Parallel two-quartet DNA GQ may exist if substantially stabilized by another molecule or structural element, including multimerization. On the other hand, we predict that isolated RNA two-quartet parallel GQs may form, albeit being weakly stable. We also show that ionic parameters and water models should be chosen with caution because some parameter combinations can cause spurious instability of GQ stems. Some in-so-far unnoticed limitations of force-field description of multiple ions inside the GQs are discussed, which compromise the capability of simulations to fully capture the effect of increase in the number of quartets on the GQ stability.

Polarizable Molecular Dynamics Simulations of Two c-kit Oncogene Promoter G-Quadruplexes: Effect of Primary and Secondary Structure on Loop and Ion Sampling

G-quadruplexes (GQs) are highly ordered nucleic acid structures that play fundamental roles in regulating gene expression and maintaining genomic stability. GQs are topologically diverse and enriched in promoter sequences of growth regulatory genes and proto-oncogenes, suggesting that they may serve as attractive targets for drug design at the level of transcription rather than inhibiting the activity of the protein products of these genes. The c-kit promoter contains three adjacent GQ-forming sequences that have proposed antagonistic effects on gene expression and thus are promising drug targets for diseases such as gastrointestinal stromal tumors, mast cell disease, and leukemia. Because GQ stability is influenced by primary structure, secondary structure, and ion interactions, a greater understanding of GQ structure, dynamics, and ion binding properties is needed to develop novel, GQ-targeting therapeutics. Here, we performed molecular dynamics simulations to systematically study the c-kit2 and c-kit* GQs, evaluating nonpolarizable and polarizable force fields (FFs) and examining the effects of base substitutions and cation type (K⁺, Na⁺, and Li⁺) on the dynamics of their isolated and linked structures. We found that the Drude polarizable FF outperformed the additive CHARMM36 FF in two- and three-tetrad GQs and solutions of KCl, NaCl, and LiCl. Drude simulations with different cations agreed with the known GQ stabilization preference (K⁺ > Na⁺ > Li⁺) and illustrated that tetrad core-ion coordination differs as a function of cation type. Finally, we showed that differences in primary and secondary structure influence loop sampling, ion binding, and core-ion energetics of GQs.

Pre-folded structures govern folding pathways of human telomeric G-quadruplexes

Understanding the mechanism by which biological macromolecules fold into their functional native conformations represents a problem of fundamental interest. DNA oligonucleotides derived from human telomeric repeat d[TAGGG(TTAGGG)₃] and d[TAGGG(TTAGGG)₃TT] fold into G-quadruplexes through diverse steps. Varying the pH and temperature by the use of nuclear magnetic resonance and other methods enabled detection of pre-folded structures that exist in solution before completely formed G-quadruplexes upon addition of cations. Pre-folded structures are in general hard to detect, however their knowledge is crucial to set up folding pathways into final structure since they are believed to be a starting point. Unexpectedly well-defined pre-folded structures composed of base triples for both oligonucleotides were detected at certain pH and temperature. These kinds of structures were up to now only hypothesized as intermediates in the folding process. All revealed pre-folded structures irrespective of the pH and temperature exhibited one common structural feature that could govern folding process.

Conformational Dynamics of Strand Register Shifts in DNA G-Quadruplexes

The complex folding energy landscape of DNA G-quadruplexes leads to numerous conformations for this functionally important class of noncanonical DNA structures. A new layer of conformational heterogeneity comes from sequences with different numbers of G-nucleotides in each of the DNA G-strands that form the four-stranded G-quartet core. The mechanisms by which G-quadruplexes transition from one folded conformation to another are currently unknown. To address this question, we studied two different G-quadruplexes, selecting a single conformation by blocking hydrogen bonding with photolabile protection groups. Upon irradiation, the block can be released and the kinetics of re-equilibration to the native conformational equilibrium can be determined by time-resolved NMR. We compared the NMR-derived refolding kinetics with data derived from thermal hysteresis folding kinetic experiments and found excellent agreement. The outlined methodological approach allows separation of K⁺-induced G-quadruplex formation and subsequent refolding and provides key insight into rate-limiting steps of G-quadruplex conformational dynamics.

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.

The crystal structure of an antiparallel chair-type G-quadruplex formed by Bromo-substituted human telomeric DNA

Human telomeric guanine-rich DNA, which could adopt different G-quadruplex structures, plays important roles in protecting the cell from recombination and degradation. Although many of these structures were determined, the chair-type G-quadruplex structure remains elusive. Here, we present a crystal structure of the G-quadruplex composed of the human telomeric sequence d[GGGTTAGG8GTTAGGGTTAGG20G] with two dG to 8Br-dG substitutions at positions 8 and 20 with syn conformation in the K+ solution. It forms a novel three-layer chair-type G-quadruplex with two linking trinucleotide loops. Particularly, T5 and T17 are coplanar with two water molecules stacking on the G-tetrad layer in a sandwich-like mode through a coordinating K+ ion and an A6•A18 base pair. While a twisted Hoogsteen A12•T10 base pair caps on the top of G-tetrad core. The three linking TTA loops are edgewise and each DNA strand has two antiparallel adjacent strands. Our findings contribute to a deeper understanding and highlight the unique roles of loop and water molecule in the folding of the G-quadruplex.

The investigation of the G-quadruplex aptamer selectivity to Pb2+ ion: a joint molecular dynamics simulation and density functional theory study

The aptamers with the ability to form a G-quadruplex structure can be stable in the presence of some ions. Hence, study of the interactions between such aptamers and ions can be beneficial to determine the highest selective aptamer toward an ion. In this article, molecular dynamics (MD) simulations and quantum mechanics (QM) calculations have been applied to investigate the selectivity of the T30695 aptamer toward Pb²⁺ in comparison with some ions. The Free Energy Landscape (FEL) analysis indicates that Pb²⁺ has remained inside the aptamer during the MD simulation, while the other ions have left it. The Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) binding energies prove that the conformational stability of the aptamer is the highest in the presence of Pb²⁺. According to the compaction parameters, the greatest compressed ion-aptamer complex, and hence, the highest ion-aptamer interaction have been induced in the presence of Pb²⁺. The contact maps clarify the closer contacts between the nucleotides of the aptamer in the presence of Pb²⁺. The density functional theory (DFT) results show that Pb²⁺ forms the most stable complex with the aptamer, which is consistent with the MD results. The QM calculations reveal that the N-H bonds and the O…H distances are the longest and the shortest, respectively, in the presence of Pb²⁺. The obtained results verify that the strongest hydrogen bonds (HBs), and hence, the most compressed aptamer structure are induced by Pb²⁺. Besides, atoms in molecules (AIM) and natural bond orbital (NBO) analyses confirm the results.Communicated by Ramaswamy H. Sarma.