Structural dynamics of possible late-stage intermediates in folding of quadruplex DNA studied by molecular simulations

Enzymatic bypass of G-quadruplex structures containing oxidative lesions

The function of many DNA processing enzymes involves sliding along the double helix or individual DNA strands. Stable secondary structures in the form of G-quadruplexes are difficult for such enzymes to bypass. We used a polymerase stop assay to determine which structural features of the human telomeric and the BCL2 promoter G-quadruplexes can stall progression of the Klenow fragment. Primer extension profiles revealed that G-quartets are effective roadblocks for the Klenow fragment, while auxiliary base pairs can be easily bypassed. Furthermore, we utilized 8-oxoguanine to simulate oxidative damage in G-rich regions and determine the effects on enzyme bypass. In rare cases, oxidative lesions reduce the level of G-quadruplex bypass. In general, however, oxidative lesions reduce G-quadruplex stability and facilitate bypassing of such G-rich regions, especially if the lesion persists in unfolding intermediates. Our findings using Klenow fragment can be extrapolated to other G-quadruplex forming sequences and enzymes that utilise a clamp-like structure to slide along DNA and are involved in processes such as gene expression regulation and telomere maintenance.

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.

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.

BG4 antibody can recognize telomeric G-quadruplexes harboring destabilizing base modifications and lesions

BG4 is a single-chain variable fragment antibody shown to bind various G-quadruplex (GQ) topologies with high affinity and specificity, and to detect GQ in cells, including GQ structures formed within telomeric TTAGGG repeats. Here, we used ELISA and single-molecule pull-down (SiMPull) detection to test how various lengths and GQ destabilizing base modifications in telomeric DNA constructs alter BG4 binding. We observed high-affinity BG4 binding to telomeric GQ independent of telomere length, although three telomeric repeat constructs that cannot form stable intramolecular GQ showed reduced affinity. A single guanine substitution with 8-aza-7-deaza-G, T, A, or C reduced affinity to varying degrees depending on the location and base type, whereas two G substitutions in the telomeric construct dramatically reduced or abolished binding. Substitution with damaged bases 8-oxoguanine and O6-methylguanine failed to prevent BG4 binding although affinity was reduced depending on lesion location. SiMPull combined with FRET revealed that BG4 binding promotes folding of telomeric GQ harboring a G to T substitution or 8-oxoguanine. Atomic force microscopy revealed that BG4 binds telomeric GQ with a 1:1 stoichiometry. Collectively, our data suggest that BG4 can recognize partially folded telomeric GQ structures and promote telomeric GQ stability.

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.

Conversion to Trimolecular G-Quadruplex by Spontaneous Hoogsteen Pairing-Based Strand Displacement Reaction between Bimolecular G-Quadruplex and Double G-Rich Probes

Bimolecular or tetramolecular G-quadruplexes (GQs) are predominantly self-assembled by the same sequence-identical G-rich oligonucleotides and usually remain inert to the strand displacement reaction (SDR) with other short G-rich invading fragments of DNA or RNA. Appealingly, in this study, we demonstrate that a parallel homomeric bimolecular GQ target of Tub10 d(CAGGGAGGGT) as the starting reactant, although completely folded in K+ solution and sufficiently stable (melting temperature of 57.7 °C), can still spontaneously accept strand invasion by a pair of short G-rich invading probes of P1 d(TGGGA) near room temperature. The final SDR product is a novel parallel heteromeric trimolecular GQ (tri-GQ) of Tub10/2P1 reassembled between one Tub10 strand and two P1 strands. Here we present, to the best of our knowledge, the first NMR solution structure of such a discrete heteromeric tri-GQ and unveil a unique mode of two probes vs one target in mutual recognition among G-rich canonical DNA oligomers. As a model system, the short invading probe P1 can spontaneously trap G-rich target Tub10 from a Watson-Crick duplex completely hybridized between Tub10 and its fully complementary strand d(ACCCTCCCTG). The Tub10 sequence of d(CAGGGAGGGT) is a fragment from the G-rich promoter region of the human β2-tubulin gene. Our findings provide new insights into the Hoogsteen pairing-based SDR between a GQ target and double invading probes of short G-rich DNA fragments and are expected to grant access to increasingly complex architectures in GQ-based DNA nanotechnology.

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.

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.

The impact of G-quadruplex dynamics on inter-tetrad electronic couplings: a hybrid computational study

The G-quadruplex is a fascinating nucleic acid motif with implications in biology, medicine, and nanotechnologies. G-quadruplexes can form in the telomeres at the edges of chromosomes and in other guanine-rich regions of the genome. They can also be engineered for exploitation as biological materials for nanodevices. Their higher stiffness and higher charge transfer rates make them better candidates in nanodevices than duplex DNA. For the development of molecular nanowires, it is important to optimize electron transport along the wire axis. One powerful basis to do so is by manipulating the structure, based on known effects that structural changes have on electron transport. Here, we investigate such effects, by a combination of classical simulations of the structure and dynamics and quantum calculations of electronic couplings. We find that this structure-function relationship is complex. A single helix shape parameter alone does not embody such complexity, but rather a combination of distances and angles between stacked bases influences charge transfer efficiency. By analyzing linear combinations of shape descriptors for different topologies, we identify the structural features that most affect charge transfer efficiency. We discuss the transferability of the proposed model and the limiting effects of inherent flexibility.

Gaining Insights on the Interactions of a Class of Decorated (2-([2,2'-Bipyridin]-6-yl)phenyl)platinum Compounds with c-Myc Oncogene Promoter G-Quadruplex and Other DNA Structures

Organometallic molecules offer some of the most promising scaffolds for interaction with G-quadruplex nucleic acids. We report the efficient synthesis of a family of organoplatinum(II) complexes, featuring a 2-([2,2'-bipyridin]-6-yl)phenyl tridentate (N∧ N∧ C) ligand, that incorporates peripheral side-chains aiming at enhancing and diversifying its interaction capabilities. These include a di-isopropyl carbamoyl amide, a morpholine ethylenamide, two enantiomeric proline imides and an oxazole. The binding affinities of the Pt-complexes were evaluated via UV-vis and fluorescence titrations, against 5 topologically-distinct DNA structures, including c-myc G-quadruplex, two telomeric (22AG) G-quadruplexes, a duplex (ds26) and a single-stranded (polyT) DNA. All compounds exhibited binding selectivity in favour of c-myc, with association constants (Ka ) in the range of 2-5×105  M-1 , lower affinity for both folds of 22AG and for ds26 and negligible affinity for polyT. Remarkable emission enhancements (up to 200-fold) upon addition of excess DNA were demonstrated by a subset of the compounds with c-myc, providing a basis for optical selectivity, since optical response to all other tested DNAs was low. A c-myc DNA-melting experiment showed significant stabilizing abilities for all compounds, with the most potent binder, the morpholine-Pt-complex, exhibiting a ΔTm >30 °C, at 1 : 5 DNA-to-ligand molar ratio. The same study implied contributions of the diverse side-chains to helix stabilization. To gain direct evidence of the nature of the interactions, mixtures of c-myc with the four most promising compounds were studied via UV Resonance Raman (UVRR) spectroscopy, which revealed end-stacking binding mode, combined with interactions of side-chains with loop nucleobase residues. Docking simulations were conducted to provide insights into the binding modes for the same four Pt-compounds, suggesting that the binding preference for two alternative orientations of the c-myc G-quadruplex thymine 'cap' ('open' vs. 'closed'), as well as the relative contributions to affinity from end-stacking and H-bonding, are highly dependent on the nature of the interacting Pt-complex side-chain.

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.

Computational and Experimental Characterization of rDNA and rRNA G-Quadruplexes

DNA G-quadruplexes in human telomeres and gene promoters are being extensively studied for their role in controlling the growth of cancer cells. G-quadruplexes have been unambiguously shown to exist both in vitro and in vivo, including in the guanine (G)-rich DNA genes encoding pre-ribosomal RNA (pre-rRNA), which is transcribed in the cell's nucleolus. Recent studies strongly suggest that these DNA sequences ("rDNA"), and the transcribed rRNA, are a potential anticancer target through the inhibition of RNA polymerase I (Pol I) in ribosome biogenesis, but the structures of ribosomal G-quadruplexes at atomic resolution are unknown and very little biophysical characterization has been performed on them to date. In the present study, circular dichroism (CD) spectroscopy is used to show that two putative rDNA G-quadruplex sequences, NUC 19P and NUC 23P and their counterpart rRNAs, predominantly adopt parallel topologies, reminiscent of the analogous telomeric quadruplex structures. Based on this information, we modeled parallel topology atomistic structures of the putative ribosomal G-quadruplexes. We then validated and refined the modeled ribosomal G-quadruplex structures using all-atom molecular dynamics (MD) simulations with the CHARMM36 force field in the presence and absence of stabilizing K⁺. Motivated by preliminary MD simulations of the telomeric parallel G-quadruplex (TEL 24P) in which the K⁺ ion is expelled, we used updated CHARMM36 force field K⁺ parameters that were optimized, targeting the data from quantum mechanical calculations and the polarizable Drude model force field. In subsequent MD simulations with optimized CHARMM36 parameters, the K⁺ ions are predominantly in the G-quadruplex channel and the rDNA G-quadruplexes have more well-defined, predominantly parallel-topology structures as compared to rRNA. In addition, NUC 19P is more structured than NUC 23P, which contains extended loops. Results from this study set the structural foundation for understanding G-quadruplex functions and the design of novel chemotherapeutics against these nucleolar targets and can be readily extended to other DNA and RNA G-quadruplexes.

DNA methylation is regulated by both the stability and topology of G-quadruplex

The methylation reaction was regulated by not only the stability of G4 but also the topology of G4.

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.

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.

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.

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.

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.

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.

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.

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.

Streamlining effects of extra telomeric repeat on telomeric DNA folding revealed by fluorescence-force spectroscopy

A human telomere ends in a single-stranded 3′ tail, composed of repeats of T₂AG₃. G-quadruplexes (GQs) formed from four consecutive repeats have been shown to possess high-structural and mechanical diversity. In principle, a GQ can form from any four repeats that are not necessarily consecutive. To understand the dynamics of GQs with positional multiplicity, we studied five and six repeats human telomeric sequence using a combination of single molecule FRET and optical tweezers. Our results suggest preferential formation of GQs at the 3′ end both in K⁺ and Na⁺ solutions, with minor populations of 5′-GQ or long-loop GQs. A vectorial folding assay which mimics the directional nature of telomere extension showed that the 3′ preference holds even when folding is allowed to begin from the 5′ side. In 100 mM K⁺, the unassociated T₂AG₃ segment has a streamlining effect in that one or two mechanically distinct species was observed at a single position instead of six or more observed without an unassociated repeat. We did not observe such streamlining effect in 100 mM Na⁺. Location of GQ and reduction in conformational diversity in the presence of extra repeats have implications in telomerase inhibition, T-loop formation and telomere end protection.

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.

Ligand binding to G-quadruplex DNA: new insights from ultraviolet resonance Raman spectroscopy

G-Quadruplexes (G4s) are noncanonical nucleic acid structures involved in the regulation of several biological processes of many organisms. The rational design of G4-targeting molecules developed as potential anticancer and antiviral therapeutics is a complex problem intrinsically due to the structural polymorphism of these peculiar DNA structures. The aim of the present work is to show how Ultraviolet Resonance Raman (UVRR) spectroscopy can complement other techniques in providing valuable information about ligand/G4 interactions in solution. Here, the binding of BRACO-19 and Pyridostatin - two of the most potent ligands - to selected biologically relevant G4s was investigated by polarized UVRR scattering at 266 nm. The results give new insights into the binding mode of these ligands to G4s having different sequences and topologies by performing an accurate analysis of peaks assigned to specific groups and their changes upon binding. Indeed, the UVRR data not only show that BRACO-19 and Pyridostatin interact with different G4 sites, but also shed light on the ligand and G4 chemical groups really involved in the interaction. In addition, UVRR results complemented by circular dichroism data clearly indicate that the binding mode of a ligand can also depend on the conformation(s) of the target G4. Overall, these findings demonstrate the utility of using UVRR spectroscopy in the investigation of G4s and G4-ligand interactions in solution.

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.

Folding intermediate states of the parallel human telomeric G-quadruplex DNA explored using Well-Tempered Metadynamics

An increasingly comprehension of the folding intermediate states of DNA G-quadruplexes (G4s) is currently an important scientific challenge, especially for the human telomeric (h-tel) G4s-forming sequences, characterized by a highly polymorphic nature. Despite the G-triplex conformation was proposed as one of the possible folding intermediates for the antiparallel and hybrid h-tel G4s, for the parallel h-tel topology with an all-anti guanine orientation, a vertical strand-slippage involving the G-triplets was proposed in previous works through microseconds-long standard molecular dynamics simulations (MDs). Here, in order to get further insights into the vertical strand-slippage and the folding intermediate states of the parallel h-tel G4s, we have carried out a Well-Tempered Metadynamics simulation (WT-MetaD), which allowed us to retrieve an ensemble of six G4s having two/G-tetrad conformations derived by the G-triplets vertical slippage. The insights highlighted in this work are aimed at rationalizing the mechanistic characterisation of the parallel h-tel G4 folding process.

Determination of a Structural Ensemble Representing the Dynamics of a G-Quadruplex DNA

It is increasingly recognized that the structures and dynamics of G-quadruplex DNA molecules are dictated by their sequences and greatly affected by environmental factors. The core guanine tetrads (G-tetrads) coordinate cations and display a strong conformational rigidity compared with that of the connecting loops. Although long loops linking the G-tetrads are typically disfavored, when present, they provide a striking illustration of the dynamics of short, single-stranded DNA regions. In addition to their role in determining the stability of the G-quadruplex state, these loops are also interesting as potential drug targets. To characterize accurately the dynamics of this DNA state, we apply here the principles of structural ensemble determination developed in the past two decades for protein molecules to DNA molecules. We thus perform extensive molecular dynamics simulations restrained with nuclear magnetic resonance residual dipolar couplings to determine a structural ensemble of the human CEB25 minisatellite G-quadruplex, which contains a connecting loop of nine nucleotides. This structural ensemble displays a wide set of arrangements for the loop and a compact, well-defined G-quadruplex core. Our results show the importance of stacking interactions in the loop and strengthen the ability of the closing base pairs to confer a large thermodynamic stability to the G-quadruplex structure.

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.

Modulation of DNA structure formation using small molecules

Genome integrity is essential for proper cell function such that genetic instability can result in cellular dysfunction and disease. Mutations in the human genome are not random, and occur more frequently at "hotspot" regions that often co-localize with sequences that have the capacity to adopt alternative (i.e. non-B) DNA structures. Non-B DNA-forming sequences are mutagenic, can stimulate the formation of DNA double-strand breaks, and are highly enriched at mutation hotspots in human cancer genomes. Thus, small molecules that can modulate the conformations of these structure-forming sequences may prove beneficial in the prevention and/or treatment of genetic diseases. Further, the development of molecular probes to interrogate the roles of non-B DNA structures in modulating DNA function, such as genetic instability in cancer etiology are warranted. Here, we discuss reported non-B DNA stabilizers, destabilizers, and probes, recent assays to identify ligands, and the potential biological applications of these DNA structure-modulating molecules.

Revealing conformational dynamics of 2'-O-methyl-RNA guanine modified G-quadruplex by replica exchange molecular dynamics

Thrombin-binding DNA aptamer (TBA) can fold into an antiparallel unimolecular G-quadruplex (G4) structure. Different types of modifications lead to various effects on the structure and stability of the G4 structure. Previous study has shown that a modified TBA (mTBA) that 2'-deoxy guanine (dG) at positions 10 and 11 in the TBA sequence were replaced by 2'-O-methyl-RNA guanine (2'OMe-G) can't fold into a well-defined G4 structure. In order to investigate the detailed structural information and probe the instability factors, we successfully employed the replica exchange molecular dynamics (REMD) to characterize the large conformational variations of the mTBA and systemically describe the influences of the 2'OMe-G on the mTBA in terms of conformation variations and the probability distributions of Hoogsteen hydrogen bonds, dihedral, sugar pucker and glycosyl torsion angle. Replacing position 10 with the 2'OMe-G (2'OMe-G10) induced a strong destabilization of the aptamer, while the 2'OMe-G at position 11(2'OMe-G11) was less destabilizing. More importantly, the glycosyl torsion angle and sugar pucker of 2'OMe-G10 were the most critical destabilization factors. These results were in good agreement with the theoretical and experimental results. Moreover, the structure information can be used as guidelines for the further design of modifications on G4 structure.

Structural and mechanistic insights into modified G-quadruplex thrombin-binding DNA aptamers

Thrombin-binding aptamer (TBA) can fold into a G-quadruplex structure necessary for interacting with thrombin. When one thymidine residue of the TGT loop at position 7 is replaced with unlocked uracil (UNA), d-isothymidine (D-isoT) or l-isothymidine (L-isoT), these modified sequences display different activities. To date, the mechanisms of how D/L-isoT and UNA influence the biological properties of TBA have not been illustrated in the literature. In this paper, we fill this gap by probing the structure variations and binding modes of these modified TBAs via molecular dynamics (MD) simulation and free energy calculation. Comparative structural analyses demonstrated that both D-IsoT and UNA changed the local conformation of TGT loop and formed stronger interactions with the target protein. Particularly, D-IsoT and UNA adopted similar conformation which can well explain their similar biological activities. In addition, the flexibility of the two TT loops were described clearly. In contrast, L-IsoT at position 7 led to an obvious tendency to unfold. Free energy calculation and the analysis of key residues energy contributions eventually provide a clear picture of interactions for further understanding of the structure-activity relationships. Collectively, our findings open the way for a rational design of modified aptamers.

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.

Extreme mechanical diversity of human telomeric DNA revealed by fluorescence-force spectroscopy

G-quadruplexes (GQs) can adopt diverse structures and are functionally implicated in transcription, replication, translation, and maintenance of telomere. Their conformational diversity under physiological levels of mechanical stress, however, is poorly understood. We used single-molecule fluorescence-force spectroscopy that combines fluorescence resonance energy transfer with optical tweezers to measure human telomeric sequences under tension. Abrupt GQ unfolding with K⁺ in solution occurred at as many as four discrete levels of force. Added to an ultrastable state and a gradually unfolding state, there were six mechanically distinct structures. Extreme mechanical diversity was also observed with Na⁺, although GQs were mechanically weaker. Our ability to detect small conformational changes at low forces enabled the determination of refolding forces of about 2 pN. Refolding was rapid and stochastically redistributed molecules to mechanically distinct states. A single guanine-to-thymine substitution mutant required much higher ion concentrations to display GQ-like unfolding and refolded via intermediates, contrary to the wild type. Contradicting an earlier proposal, truncation to three hexanucleotide repeats resulted in a single-stranded DNA-like mechanical behavior under all conditions, indicating that at least four repeats are required to form mechanically stable structures.

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.

Dynamics of physiologically relevant noncanonical DNA structures: an overview from experimental and theoretical studies

DNA is a complex molecule with phenomenal inherent plasticity and the ability to form different hydrogen bonding patterns of varying stabilities. These properties enable DNA to attain a variety of structural and conformational polymorphic forms. Structurally, DNA can exist in single-stranded form or as higher-order structures, which include the canonical double helix as well as the noncanonical duplex, triplex and quadruplex species. Each of these structural forms in turn encompasses an ensemble of dynamically heterogeneous conformers depending on the sequence composition and environmental context. In vivo, the widely populated canonical B-DNA attains these noncanonical polymorphs during important cellular processes. While several investigations have focused on the structure of these noncanonical DNA, studying their dynamics has remained nontrivial. Here, we outline findings from some recent advanced experimental and molecular simulation techniques that have significantly contributed toward understanding the complex dynamics of physiologically relevant noncanonical forms of DNA.

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.

Topological impact of noncanonical DNA structures on Klenow fragment of DNA polymerase

Noncanonical DNA structures that stall DNA replication can cause errors in genomic DNA. Here, we investigated how the noncanonical structures formed by sequences in genes associated with a number of diseases impacted DNA polymerization by the Klenow fragment of DNA polymerase. Replication of a DNA sequence forming an i-motif from a telomere, hypoxia-induced transcription factor, and an insulin-linked polymorphic region was effectively inhibited. On the other hand, replication of a mixed-type G-quadruplex (G4) from a telomere was less inhibited than that of the antiparallel type or parallel type. Interestingly, the i-motif was a better inhibitor of replication than were mixed-type G4s or hairpin structures, even though all had similar thermodynamic stabilities. These results indicate that both the stability and topology of structures formed in DNA templates impact the processivity of a DNA polymerase. This suggests that i-motif formation may trigger genomic instability by stalling the replication of DNA, causing intractable diseases.