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

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.

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

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

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.

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 (K_a ) in the range of 2-5×10⁵  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 ΔT_m >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.

Dynamic interaction of BRCA2 with telomeric G-quadruplexes underlies telomere replication homeostasis

BRCA2-deficient cells precipitate telomere shortening upon collapse of stalled replication forks. Here, we report that the dynamic interaction between BRCA2 and telomeric G-quadruplex (G4), the non-canonical four-stranded secondary structure, underlies telomere replication homeostasis. We find that the OB-folds of BRCA2 binds to telomeric G4, which can be an obstacle during replication. We further demonstrate that BRCA2 associates with G-triplex (G3)-derived intermediates, which are likely to form during direct interconversion between parallel and non-parallel G4. Intriguingly, BRCA2 binding to G3 intermediates promoted RAD51 recruitment to the telomere G4. Furthermore, MRE11 resected G4-telomere, which was inhibited by BRCA2. Pathogenic mutations at the OB-folds abrogated the binding with telomere G4, indicating that the way BRCA2 associates with telomere is innate to its tumor suppressor activity. Collectively, we propose that BRCA2 binding to telomeric G4 remodels it and allows RAD51-mediated restart of the G4-driven replication fork stalling, simultaneously preventing MRE11-mediated breakdown of telomere.

G-quadruplex (G4) can be formed in telomeric DNA. Here the authors show that BRCA2 interacts with telomere G4 structure generated during telomere replication, protecting telomere from nuclease attack.

Never Cared for What They Do: High Structural Stability of Guanine-Quadruplexes in the Presence of Strand-Break Damage

DNA integrity is an important factor that assures genome stability and, more generally, the viability of cells and organisms. In the presence of DNA damage, the normal cell cycle is perturbed when cells activate their repair processes. Although efficient, the repair system is not always able to ensure complete restoration of gene integrity. In these cases, mutations not only may occur, but the accumulation of lesions can either lead to carcinogenesis or reach a threshold that induces apoptosis and programmed cell death. Among the different types of DNA lesions, strand breaks produced by ionizing radiation are the most toxic due to the inherent difficultly of repair, which may lead to genomic instability. In this article we show, by using classical molecular simulation techniques, that compared to canonical double-helical B-DNA, guanine-quadruplex (G4) arrangements show remarkable structural stability, even in the presence of two strand breaks. Since G4-DNA is recognized for its regulatory roles in cell senescence and gene expression, including oncogenes, this stability may be related to an evolutionary cellular response aimed at minimizing the effects of ionizing radiation.

G4-quadruplex-binding proteins: review and insights into selectivity

There are over 700,000 putative G4-quadruplexes (G4Qs) in the human genome, found largely in promoter regions, telomeres, and other regions of high regulation. Growing evidence links their presence to functionality in various cellular processes, where cellular proteins interact with them, either stabilizing and/or anchoring upon them, or unwinding them to allow a process to proceed. Interest in understanding and manipulating the plethora of processes regulated by these G4Qs has spawned a new area of small-molecule binder development, with attempts to mimic and block the associated G4-binding protein (G4BP). Despite the growing interest and focus on these G4Qs, there is limited data (in particular, high-resolution structural information), on the nature of these G4Q-G4BP interactions and what makes a G4BP selective to certain G4Qs, if in fact they are at all. This review summarizes the current literature on G4BPs with regards to their interactions with G4Qs, providing groupings for binding mode, drawing conclusions around commonalities and highlighting information on specific interactions where available.

Multimeric G-quadruplexes: A review on their biological roles and targeting

In human cells, nucleic acids adopt several non-canonical structures that regulate key cellular processes. Among them, G-quadruplexes (G4s) are stable structures that form in guanine-rich regions in vitro and in cells. G4 folded/unfolded state shapes numerous cellular processes, including genome replication, transcription, and translation. Moreover, G4 folding is involved in genomic instability. G4s have been described to multimerize, forming high-order structures in both DNA and/or RNA strands. Multimeric G4s can be formed by adjacent intramolecular G4s joined by stacking interactions or connected by short loops. Multimeric G4s can also originate from the assembly of guanines embedded on independent DNA or RNA strands. Notably, crucial regions of the human genome, such as the 3'-terminal overhang of the telomeric DNA as well as the open reading frame of genes involved in the preservation of neuron viability in the human central and peripheral nervous system are prone to form multimeric G4s. The biological importance of such structures has been recently described, with multimeric G4s playing potentially protective or deleterious effects in the pathogenic cascade of various diseases. Here, we portray the multifaceted scenario of multimeric G4s, in terms of structural properties, biological roles, and targeting strategies.

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.

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.

Frontiers in G-Quadruplex Therapeutics in Cancer: Selection of Small Molecules, Peptides and Aptamers

G-quadruplex, a unique secondary structure in nucleic acids found throughout human genome, elicited widespread interest in the field of therapeutic research. Being present in key regulatory regions of oncogenes, RNAs and telomere, G-quadruplex structure regulates transcription, translation, splicing etc. Changes in its structure and stability leads to differential expression of oncogenes causing cancer. Thus, targeting G-Quadruplex structures with small molecules/other biologics has shown elevated research interest. Covering previous reports, in this review we try to enlighten the facts on the structural diversity in G-quadruplex ligands aiming to provide newer insights to design first-in-class drugs for the next generation cancer treatment.

Unraveling the binding characteristics of small ligands to telomeric DNA by pressure modulation

Recently, non-canonical DNA structures, such as G-quadruplexes (GQs), were found to be highly pressure sensitive, suggesting that pressure modulation studies can provide additional mechanistic details of such biomolecular systems. Using FRET and CD spectroscopy as well as binding equilibrium measurements, we investigated the effect of pressure on the binding reaction of the ligand ThT to the quadruplex 22AG in solutions containing different ionic species and a crowding agent mimicking the intracellular milieu. Pressure modulation helped us to identify the different conformational substates adopted by the quadruplex at the different solution conditions and to determine the volumetric changes during complex formation and the conformational transitions involved. The magnitudes of the binding volumes are a hallmark of packing defects and hydrational changes upon ligand binding. The conformational substates of the GQ as well as the binding strength and the stoichiometry of complex formation depend strongly on the solution conditions as well as on pressure. High hydrostatic pressure can also impact GQs inside living cells and thus affect expression of genetic information in deep sea organisms. We show that sub-kbar pressures do not only affect the conformational dynamics and structures of GQs, but also their ligand binding reactions.

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.

In Silico Food-Drug Interaction: A Case Study of Eluxadoline and Fatty Meal

Food-drug interaction is an infrequently considered aspect in clinical practice. Usually, drugs are taken together with meals and what follows may adversely affect pharmacokinetic and pharmacodynamic properties, and hence, the therapeutic effects. In this study, a computational protocol was proposed to explain the different assimilations of two µ-receptors agonists, eluxadoline and loperamide, with a peculiar pharmacokinetic profile. Compared to loperamide, eluxadoline is absorbed less after the intake of a fatty meal, and the LogP values do not explain this event. Firstly, keeping in mind the different pH in the intestinal tract, the protonation states of both compounds were calculated. Then, all structures were subjected to a conformational search by using MonteCarlo and Molecular Dynamics methods, with solvation terms mimicking the water and weak polar solvent (octanol). Both computational results showed that eluxadoline has less conformational freedom in octanol, unlike loperamide, which exhibits constant behavior in both solvents. Therefore, we hypothesize that fatty meal causes the "closure" of the eluxadoline molecule to prevent the exposure of the polar groups and their interaction with water, necessary for the drug absorption. Based on our results, this work could be a reasonable "case study", useful for future investigation of the drug pharmacokinetic profile.

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.