Molecular dynamics simulations and their application to four-stranded DNA

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.

Differences in Conformational Sampling and Intrinsic Electric Fields Drive Ion Binding in Telomeric and TERRA G-Quadruplexes

The formation of G-quadruplexes (GQs) occurs in guanine-rich sequences of DNA and RNA, producing highly stable and structurally diverse noncanonical nucleic acid structures. GQs play crucial roles in regulating transcription, translation, and replication and maintaining the genome, among others; thus, changes to their structures can lead to diseases such as cancer. Previous studies using polarizable molecular dynamics simulations have shown differences in ion binding properties between telomeric and telomeric repeat-containing RNA GQs despite architectural similarities. Here, we used volume-based metadynamics and repulsive potential simulations in conjunction with polarizable force fields to quantify the impact of ion binding on the GQ dynamics and ion binding free energies. Furthermore, we describe how GQs exert electric fields on their surroundings to link dynamics with variations in the electronic structure. Our findings provide new insights into the energetic, physical, and conformational properties of GQs and expose subtle but important differences between DNA and RNA GQs with the same fold.

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.

Differences in Conformational Sampling and Intrinsic Electric Fields Drive Ion Binding in Telomeric and TERRA G-Quadruplexes

The formation of G-quadruplexes (GQs) occurs in guanine-rich sequences of DNA and RNA, producing highly stable and structurally diverse noncanonical nucleic acid structures. GQs play crucial roles in regulating transcription, translation, and replication; and maintaining the genome, among others, thus changes to their structures can lead to diseases such as cancer. Previous studies using polarizable molecular dynamics simulations have shown differences in ion binding properties between telomeric and TERRA GQs despite architectural similarities. Here, we used volume-based metadynamics and repulsive potential simulations in conjunction with polarizable force fields to quantify the impact of ion binding on GQ dynamics and ion binding free energies. Furthermore, we describe how GQs exert electric fields on their surroundings to link dynamics with variations in electronic structure. Our findings provide new insights into the energetic, physical, and conformational properties of GQs and expose subtle, but important, differences between DNA and RNA GQs with the same fold.

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.

Early steps of oxidative damage in DNA quadruplexes are position-dependent: Quantum mechanical and molecular dynamics analysis of human telomeric sequence containing ionized guanine

Guanine radical cation (G^•+) is a key intermediate in many oxidative processes occurring in nucleic acids. Here, by combining mixed Quantum Mechanical/Molecular Mechanics calculations and Molecular Dynamics (MD) simulations, we study how the structural behaviour of a tract GGG(TTAGGG)₃ (hereafter Tel21) of the human telomeric sequence, folded in an antiparallel quadruple helix, changes when one of the G bases is ionized to G^•+ (Tel21⁺). Once assessed that the electron-hole is localized on a single G, we perform MD simulations of twelve Tel21⁺ systems, differing in the position of G^•+ in the sequence. When G^•+ is located in the tetrad adjacent to the diagonal loop, we observe substantial structural rearrangements, which can decrease the electrostatic repulsion with the inner Na⁺ ions and increase the solvent exposed surface of G^•+. Analysis of solvation patterns of G^•+ provides new insights on the main reactions of G^•+, i.e. the deprotonation at two different sites and hydration at the C8 atom, the first steps of the processes producing 8oxo-Guanine. We suggest the main structural determinants of the relative reactivity of each position and our conclusions, consistent with the available experimental trends, can help rationalizing the reactivity of other G-quadruplex topologies.

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.

Benchmark Force Fields for the Molecular Dynamic Simulation of G-Quadruplexes

G-quadruplexes have drawn widespread attention for serving as a potential anti-cancer target and their application in material science. Molecular dynamics (MD) simulation is the key theoretical tool in the study of GQ's structure-function relationship. In this article, we systematically benchmarked the five force fields of parmbsc0, parmbsc1, OL15, AMOEBA, and Drude2017 on the MD simulation of G-quadruplex from four aspects: structural stability, central ion channel stability, description of Hoogsteen hydrogen bond network, and description of the main chain dihedral angle. The results show that the overall performance of the Drude force field is the best. Although there may be a certain over-polarization effect, it is still the best choice for the MD simulation of G-quadruplexes.

dUMP/F-dUMP Binding to Thymidylate Synthase: Human Versus Mycobacterium tuberculosis

Thymidylate synthase is an enzyme that catalyzes deoxythymidine monophosphate (dTMP) synthesis from substrate deoxyuridine monophosphate (dUMP). Thymidylate synthase of Mycobacterium tuberculosis (MtbThyX) is structurally distinct from its human analogue human thymidylate synthase (hThyA), thus drawing attention as an attractive drug target for combating tuberculosis. Fluorodeoxyuridylate (F-dUMP) is a successful inhibitor of both MtbThyX and hThyA, thus limited by poor selectivity. Understanding the dynamics and energetics associated with substrate/inhibitor binding to thymidylate synthase in atomic details remains a fundamental unsolved problem, which is necessary for a new selective inhibitor design. Structural studies of MtbThyX and hThyA bound substrate/inhibitor complexes not only revealed the extensive specific interaction network between protein and ligands but also opened up the possibility of directly computing the energetics of the substrate versus inhibitor recognition. Using experimentally determined structures as a template, we report extensive computer simulations (∼4.5 μs) that allow us to quantitatively estimate ligand selectivity (dUMP vs F-dUMP) by MtbThyX and hThyA. We show that MtbThyX prefers deprotonated dUMP (enolate form) as the substrate, whereas hThyA binds to the keto form of dUMP. Computed energetics clearly show that MtbThyX is less selective between dUMP and F-dUMP, favoring the latter, relative to hThyA. The simulations reveal the role of tyrosine at position 135 (Y135) of hThyA in amplifying the selectivity. The protonation state of the pyrimidine base of the ligand (i.e., keto or enolate) seems to have no role in MtbThyX ligand selectivity. A molecular gate (consists of Y108, K165, H203, and a water molecule) restricts water accessibility and offers a desolvated dry ligand-binding pocket for MtbThyX. The ligand-binding pocket of hThyA is relatively wet and exposed to bulk water.

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.

To probe the binding pathway of a selective compound (D089-0563) to c-MYC Pu24 G-quadruplex using free ligand binding simulations and Markov state model analysis

D089-0563 is a highly promising anti-cancer compound that selectively binds the transcription-silencing G-quadruplex element (Pu27) at the promoter region of the human c-MYC oncogene; however, its binding mechanism remains elusive.

Studying the Potassium-Induced G-Quadruplex DNA Folding Process Using Microscale Thermophoresis

Guanine (G) quadruplexes (G4s) can be formed by G-rich sequences when stabilized by the binding of cations (typically K⁺ or Na⁺) and play an essential role in replication, recombination, transcription, and telomere maintenance. Understanding of the G4 folding process is crucial for determining their cellular functions. However, G4-K⁺ interactions and folding pathways are still not well understood. By using human telomeric G4 (hTG4) as an example, two binding states corresponding to two K⁺ cations binding to hTG4 were distinguished clearly and fitted precisely. The basic binding parameters during G4-K⁺ interactions were measured and calculated by taking advantage of microscale thermophoresis (MST), which monitors the changes in charge and size at the same time. The G-hairpin and G-triplex have been suggested as intermediates during G4 folding and unfolding. We further analyzed the equilibrium dissociation constants of 10 possible folding intermediates using MST; thus, the energetically favorable folding/unfolding pathways were proposed. The results might not only shed new light on G4-K⁺ interactions and G4 folding pathways but also provide an example for experimentally studying DNA-ion interactions.

Human-specific features of the G-quadruplex in the androgen receptor gene promoter: A comparative structural and dynamics study

The androgen receptor (AR) promoter contains guanine-rich regions that are able to fold into polymorphic G-quadruplex (GQ) structures, and whose deletion decreases AR gene transcription. Our attention was focused on this region because of the frequent termination of sequencing reactions during promoter methylation studies. UV and circular dichroism (CD) spectroscopy of synthetic oligonucleotides encompassing these guanine-rich regions suggested a parallel quadruplex topology with three guanine quartets and three side loops in the three cases. Melting curves revealed a lower thermostability of the human GQ compared to the rat/mouse QG structures, which is attributed to the presence of a longer central loop in the former. One molecular model is proposed for the highly similar sequences in the rat/mouse. Due to the polymorphism resulting from possible arrangements of the guanine tracts, two models were derived for the human GQ. Molecular dynamics (MD) simulations determined that both models for the human GQ had higher flexibility and lower stability than the rodent GQ models. These properties result from the presence of a longer central loop in the human GQ models, which contains 11 and 13 nucleotides, in comparison to the 2-nucleotide long loop in the rat/mouse GQ. Overall, the unveiled structural and dynamics features provide sufficient detail for the intelligent design of drugs targeting the human AR promoter.

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

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

A coarse-grained model for DNA origami

Modeling tools provide a valuable support for DNA origami design. However, current solutions have limited application for conformational analysis of the designs. In this work we present a tool for a thorough study of DNA origami structure and dynamics. The tool is based on a novel coarse-grained model dedicated to geometry optimization and conformational analysis of DNA origami. We explored the ability of the model to predict dynamic behavior, global shapes, and fine details of two single-layer systems designed in hexagonal and square lattices using atomic force microscopy, Förster resonance energy transfer spectroscopy, and all-atom molecular dynamic simulations for validation of the results. We also examined the performance of the model for multilayer systems by simulation of DNA origami with published cryo-electron microscopy and atomic force microscopy structures. A good agreement between the simulated and experimental data makes the model suitable for conformational analysis of DNA origami objects. The tool is available at http://vsb.fbb.msu.ru/cosm as a web-service and as a standalone version.

Probing the Limits of Supramolecular G-Quadruplexes Using Atomistic Molecular Dynamics Simulations

Guanosine and related derivatives self-assemble in the presence of cations like potassium into supramolecular G-quadruplexes (SGQs), where four guanine moieties form planar tetrads (T) that coaxially stack into columnar aggregates with broad size distributions. However, SGQs made from 8-aryl-2'-deoxyguanosine derivatives (8ArGs), form mostly octamers, or two-tetrad (2T)-SGQs, while some form dodecamers (3T-SGQs), or hexadecamers (4T-SGQs), and none reported to date form higher assemblies. A theoretical model that addresses the configurational space available for the multiple pathways available for 8ArGs to self-assemble into SGQs is used to frame a series of molecular dynamics simulations (MDS) with selected SGQs. Some key insights from this work include: (a) The predicted entropic costs are not significantly higher for SGQs with more subunits due to their hierarchical assembly pathways; (b) The multiple isomeric SGQs vary in the interfacial contacts between consecutive tetrads, due to their two distinct sides (head, h; tail, t), with the MDS supporting the predicted order of stability of hh > ht > tt for octamers. (c) Such order also applies to dodecamers and hexadecamers, but with context-dependent exceptions due to strong allosteric effects. (d) The main factor disfavoring the tt interface is the repulsive dipolar interactions between the O4' from ribose moieties on adjacent tetrads. (e) SGQs with 5 or more tetrads are disfavored because the attractive interactions are not large or strong enough to overcome the many repulsive forces resulting from the addition of further tetrads. We expect these findings provide some guidelines to enable the further development of SGQs into functional materials.

Rotation of Guanine Amino Groups in G-Quadruplexes: A Probe for Local Structure and Ligand Binding

Nucleic acids are dynamic molecules whose functions may depend on their conformational fluctuations and local motions. In particular, amino groups are dynamic components of nucleic acids that participate in the formation of various secondary structures such as G-quadruplexes. Here, we present a cost-efficient NMR method to quantify the rotational dynamics of guanine amino groups in G-quadruplex nucleic acids. An isolated spectrum of amino protons from a specific tetrad-bound guanine can be extracted from the nuclear Overhauser effect spectroscopy spectrum based on the close proximity between the intra-residue imino and amino protons. We apply the method in different structural contexts of G-quadruplexes and their complexes. Our results highlight the role of stacking and hydrogen-bond interactions in restraining amino-group rotation. The measurement of the rotation rate of individual amino groups could give insight into the dynamic processes occurring at specific locations within G-quadruplex nucleic acids, providing valuable probes for local structure, dynamics, and ligand binding.

Role of Alkali Metal Ions in G-Quadruplex Nucleic Acid Structure and Stability

G-quadruplexes are guanine-rich nucleic acids that fold by forming successive quartets of guanines (the G-tetrads), stabilized by intra-quartet hydrogen bonds, inter-quartet stacking, and cation coordination. This specific although highly polymorphic type of secondary structure deviates significantly from the classical B-DNA duplex. G-quadruplexes are detectable in human cells and are strongly suspected to be involved in a number of biological processes at the DNA and RNA levels. The vast structural polymorphism exhibited by G-quadruplexes, together with their putative biological relevance, makes them attractive therapeutic targets compared to canonical duplex DNA. This chapter focuses on the essential and specific coordination of alkali metal cations by G-quadruplex nucleic acids, and most notably on studies highlighting cation-dependent dissimilarities in their stability, structure, formation, and interconversion. Section 1 surveys G-quadruplex structures and their interactions with alkali metal ions while Section 2 presents analytical methods used to study G-quadruplexes. The influence of alkali cations on the stability, structure, and kinetics of formation of G-quadruplex structures of quadruplexes will be discussed in Sections 3 and 4. Section 5 focuses on the cation-induced interconversion of G-quadruplex structures. In Sections 3 to 5, we will particularly emphasize the comparisons between cations, most often K(+) and Na(+) because of their prevalence in the literature and in cells.

Structural Dynamics of Thrombin-Binding DNA Aptamer d(GGTTGGTGTGGTTGG) Quadruplex DNA Studied by Large-Scale Explicit Solvent Simulations

The thrombin-binding aptamer (15-TBA) is a 15-mer DNA oligonucleotide with sequence d(GGTTGGTGTGGTTGG). 15-TBA folds into a quadruplex DNA (G-DNA) structure with two planar G-quartets connected by three single-stranded loops. The arrangement of the 15-TBA-thrombin complex is unclear, particularly with respect to the precise 15-TBA residues that interact with the thrombin structure. Our present understanding suggests either the 15-TBA single stranded loops containing sequential thymidines (TT) or alternatively a single-stranded loop, containing a guanine flanked by 2 thymidines (TGT), physically associates with thrombin protein. In the present study, the explicit solvent molecular dynamics (MD) simulation method was utilized to further analyze the 15-TBA-thrombin three-dimensional structure. Functional annotation of the loop residues was made with long simulations in the parmbsc0 force field. In total, the elapsed time of simulations carried out in this study exceeds 12 microseconds, substantially surpassing previous G-DNA simulation reports. Our simulations suggest that the TGT-loop function is to stabilize the structure of the aptamer, while the TT-loops participate in direct binding to thrombin. The findings of the present report advance our understanding of the molecular structure of the 15-TBA-thrombin structure further enabling the construction of biosensors for aptamer bases and the development of anticoagulant agents.

Extended molecular dynamics of a c-kit promoter quadruplex

The 22-mer c-kit promoter sequence folds into a parallel-stranded quadruplex with a unique structure, which has been elucidated by crystallographic and NMR methods and shows a high degree of structural conservation. We have carried out a series of extended (up to 10 μs long, ∼50 μs in total) molecular dynamics simulations to explore conformational stability and loop dynamics of this quadruplex. Unfolding no-salt simulations are consistent with a multi-pathway model of quadruplex folding and identify the single-nucleotide propeller loops as the most fragile part of the quadruplex. Thus, formation of propeller loops represents a peculiar atomistic aspect of quadruplex folding. Unbiased simulations reveal μs-scale transitions in the loops, which emphasizes the need for extended simulations in studies of quadruplex loops. We identify ion binding in the loops which may contribute to quadruplex stability. The long lateral-propeller loop is internally very stable but extensively fluctuates as a rigid entity. It creates a size-adaptable cleft between the loop and the stem, which can facilitate ligand binding. The stability gain by forming the internal network of GA base pairs and stacks of this loop may be dictating which of the many possible quadruplex topologies is observed in the ground state by this promoter quadruplex.

Investigation of the nanoviscosity effect of a G-quadruplex and single-strand DNA using fluorescence correlation spectroscopy

Guanine (G)-quadruplexes are of interest because of their presence in the telomere sequence and the oncogene promoter region. Their diffusion and change of structure, especially in high viscosity solutions, are important for understanding their dynamics. G-quadruplexes may have less effective viscosity (nanoviscosity) when they are smaller than the solvent molecules. In this paper, we report the difference in the diffusion dynamics of the G-rich DNA sequences of single-strand DNA (ssDNA) and the G-quadruplex in aqueous, sucrose, and polyethylene glycol (PEG) solutions. From experiments with aqueous and sucrose solutions, we confirm that a simple diffusion model according to the viscosity is appropriate. In the PEG experiments, the nanoviscosity effect is observed according to PEG's molecular weight. In the PEG 200 solution, both the ssDNA and the G-quadruplex possess macroviscosity. In the PEG 10,000 solution, the G-quadruplex possesses nanoviscosity and the ssDNA possesses macroviscosity, whereas, in the PEG 35,000 solution, both ssDNA and the G-quadruplex possess nanoviscosity. The experimental results are consistent with the theoretical predictions.

Tuning supramolecular G-quadruplexes with mono- and divalent cations

Supramolecular G-quadruplexes (SGQs) are formed via the cation promoted self-assembly of guanine derivatives into stacks of planar hydrogen-bonded tetramers. Here, we present results on the formation of SGQs made from the 8-(m-acetylphenyl)-2'-deoxyguanosine (mAGi) derivative in the presence of various mono- and divalent cations. NMR and HR ESI-MS data indicate that varying the cation can efficiently tune the molecularity, the fidelity and stability (thermal and kinetic) of the resulting SGQs. The results show that, parallel to the previously reported potassium-templated hexadecamer (mAGi₁₆·3K⁺), Na⁺, Rb⁺ and [Formula: see text] also promote the formation of similar supramolecules with high fidelity and molecularity. In contrast, the divalent cations Pb²⁺, Sr²⁺ and Ba²⁺ template the formation of octamers (mAGi₈), with the latter two inducing higher thermal stabilities. Molecular dynamics simulations for the hexadecamers containing monovalent cations enabled critical insights that help explain the experimental observations.

Molecular modeling of nucleic Acid structure: electrostatics and solvation

This unit presents an overview of computer simulation techniques as applied to nucleic acid systems, ranging from simple in vacuo molecular modeling techniques to more complete all-atom molecular dynamics treatments that include an explicit representation of the environment. The third in a series of four units, this unit focuses on critical issues in solvation and the treatment of electrostatics. UNITS 7.5 & 7.8 introduced the modeling of nucleic acid structure at the molecular level. This included a discussion of how to generate an initial model, how to evaluate the utility or reliability of a given model, and ultimately how to manipulate this model to better understand its structure, dynamics, and interactions. Subject to an appropriate representation of the energy, such as a specifically parameterized empirical force field, the techniques of minimization and Monte Carlo simulation, as well as molecular dynamics (MD) methods, were introduced as a way of sampling conformational space for a better understanding of the relevance of a given model. This discussion highlighted the major limitations with modeling in general. When sampling conformational space effectively, difficult issues are encountered, such as multiple minima or conformational sampling problems, and accurately representing the underlying energy of interaction. In order to provide a realistic model of the underlying energetics for nucleic acids in their native environments, it is crucial to include some representation of solvation (by water) and also to properly treat the electrostatic interactions. These subjects are discussed in detail in this unit.

Multiscale simulations of human telomeric G-quadruplex DNA

We present a coarse-grain (CG) model of human telomeric G-quadruplex, obtained using the inverse Monte Carlo (IMC) and iterative Boltzmann inversion (IBI) techniques implemented within the software package called MagiC. As a starting point, the 2HY9 human telomeric [3 + 1] hybrid, a 26-nucleobase sequence, was modeled performing a 1 μs long atomistic molecular dynamics (MD) simulation. The chosen quadruplex includes two kinds of loops and all possible combinations of relative orientations of guanine strands that can be found in quadruplexes. The effective CG potential for a one bead per nucleotide model has been developed from the radial distribution functions of this reference system. The obtained potentials take into account explicitly the interaction with counterions, while the effect of the solvent is included implicitly. The structural properties of the obtained CG model of the quadruplex provided a perfect match to those resulting from the reference atomistic MD simulation. The same set of interaction potentials was then used to simulate at the CG level another quadruplex topology (PDB id 1KF1 ) that can be formed by the human telomeric DNA sequence. This quadruplex differs from 2HY9 in the loop topology and G-strand relative orientation. The results of the CG MD simulations of 1KF1 are very encouraging and suggest that the CG model based on 2HY9 can be used to simulate quadruplexes with different topologies. The CG model was further applied to a higher order human telomeric quadruplex formed by the repetition, 20 times, of the 1KF1 quadruplex structure. In all cases, the developed model, which to the best of our knowledge is the first model of quadruplexes at the CG level presented in the literature, reproduces the main structural features remarkably well.

Molecular dynamics simulations identify time scale of conformational changes responsible for conformational selection in molecular recognition of HIV-1 transactivation responsive RNA

The HIV-1 Tat protein and several small molecules bind to HIV-1 transactivation responsive RNA (TAR) by selecting sparsely populated but pre-existing conformations. Thus, a complete characterization of TAR conformational ensemble and dynamics is crucial to understand this paradigmatic system and could facilitate the discovery of new antivirals targeting this essential regulatory element. We show here that molecular dynamics simulations can be effectively used toward this goal by bridging the gap between functionally relevant time scales that are inaccessible to current experimental techniques. Specifically, we have performed several independent microsecond long molecular simulations of TAR based on one of the most advanced force fields available for RNA, the parmbsc0 AMBER. Our simulations are first validated against available experimental data, yielding an excellent agreement with measured residual dipolar couplings and order parameter S(2). This contrast with previous molecular dynamics simulations (Salmon et al., J. Am. Chem. Soc. 2013 135, 5457-5466) based on the CHARMM36 force field, which could achieve only modest accord with the experimental RDC values. Next, we direct the computation toward characterizing the internal dynamics of TAR over the microsecond time scale. We show that the conformational fluctuations observed over this previously elusive time scale have a strong functionally oriented character in that they are primed to sustain and assist ligand binding.

Quantum mechanical calculations unveil the structure and properties of the absorbing and emitting excited electronic states of guanine quadruplex

Herein, a full quantum mechanical study, in solution, of several models of guanine-quadruplex helices, both parallel and antiparallel, containing up to eight guanine residues, in their electronic excited state is reported. By exploiting TD-DFT calculations and including solvent effects by the polarizable continuum model, we provide the first atomistic description of the processes triggered by the absorption of UV light, reproducing and assigning the experimental optical and electronic circular dichroism spectra. The absorbing excited states are delocalized over multiple bases, whereas emission involves a stacked guanine dimer or a monomer. Several states, with a varying degree of localization and charge-transfer character, rule the photoexcited dynamics, which are deeply affected by the quadruplex topology. The lowest excited-state minimum for parallel quadruplex is an asymmetric excimer involving two stacked guanines, with a small charge transfer character, whereas for the anti-parallel structure, with the same topology of the thrombin binding aptamer, it is a fully symmetric excimer, characterized by a strong decrease of the stacking distance. A monomer-like decay path is the most relevant nonradiative decay pathway. Insights on the effect of the ions (K(+) or Na(+)) on the excited state decay are also provided.

Structure and conformational dynamics of a stacked dimeric G-quadruplex formed by the human CEB1 minisatellite

CEB1 is a highly polymorphic human minisatellite. In yeast, the size variation of CEB1 tandem arrays has been associated with the capacity of the motif to form G-quadruplexes. Here we report on the NMR solution structure of a G-quadruplex formed by the CEB1 DNA G-rich fragment d(AGGGGGGAGGGAGGGTGG), harboring several G-tracts including one with six continuous guanines. This sequence forms a dimeric G-quadruplex involving the stacking of two subunits, each being a unique snapback parallel-stranded scaffold with three G-tetrad layers, three double-chain-reversal loops, and a V-shaped loop. The two subunits are stacked at their 5'-end tetrads, and multiple stacking rotamers may be present due to a high symmetry at the stacking interface. There is a conformational exchange in the millisecond time scale involving a swapping motion between two bases of the six-guanine tract. Our results not only add to the understanding of how the G-quadruplex formation in human minisatellite leads to genetic instability but also address the fundamental questions regarding stacking of G-quadruplexes and how a long continuous G-tract participates in the structure and conformational dynamics of G-quadruplexes.

Unfolding and conformational variations of thrombin-binding DNA aptamers: synthesis, circular dichroism and molecular dynamics simulations

Thrombin-binding DNA aptamer (TBA), with a consensus 15-base sequence: d(GGTTGGTGTGGTTGG), can fold into an antiparallel unimolecular G-quadruplex structure that is necessary for its interaction with thrombin. For the first time, using steered molecular dynamics (SMD) simulations, we have successfully simulated the unfolding process of native TBA G-quadruplex. The unfolding pathway proposed is in agreement with previously reported experimental NMR data. Moreover, the critical intermediate structure in the unfolding pathway, predicted by the NMR results, was identified. The structural characteristics of several TBA oligonucleotides modified with locked nucleoside (LNA) or 2'-O-methyl-nucleoside (MNA) at different positions and number were also investigated by CD spectroscopy. An oligonucleotide substituted with either LNA or MNA at position 2 folds into a native-like G-quadruplex, while doubly substituted derivatives of TBA where LNA or MNA is incorporated at positions 11 and 14 are no longer able to form a G-quadruplex. Starting from the same initial intermediate structure, we successfully overcame sampling limitations, and simulated the large conformational variations in structures of native TBA and modified TBAs by classic MD. Analysis of the models showed that inversion of the glycosyl orientation at position 14 contributes significantly to the disruption of G-quadruplex formation in both of the di-substituted modified TBA systems. Our calculations provide a simple and reliable theoretical model that can be used to investigate and predict the effects of the modifications of an oligonucleotide on the resultant G-quadruplex structure. In addition, the computational protocol described can also be applied for other G-quadruplex systems.

Deprotonation mechanism of a single-stranded DNA i-motif

An atomistic Molecular Dynamics simulation to study the unfolding and deprotonation mechanism of a single-stranded and fully protonated DNA i-motif.

The critical effect of polarization on the dynamical structure of guanine quadruplex DNA

Guanine quadruplex DNA (G-DNA), found in eukaryotic telomeres and recently in non-telomeric genomic DNA, plays important biological roles and their structures are being explored as potential targets for therapeutic intervention. Since the quadruplex structure of G-DNA is stabilized by cations, electrostatic interaction is expected to play important roles in the dynamical structure of G-DNA. In current work, MD simulation was carried out to study the dynamical structure of a special G-DNA (with sequence d(G(4)T(4)G(4))) complexed with five K(+) ions. In order to properly include polarization in MD simulation, a new set of polarized nucleic acid specific charge based on fragment quantum chemistry calculation was developed for G-DNA. Our study showed that polarization of the nucleobases by K(+) enhanced electrostatic attraction between the base and ions. This increased attractive interaction is critical to stabilizing the stem-loop junction ions in G-DNA. Without this polarization effect, as in MD simulation using a standard (nonpolarizable) force field, the top and bottom cations escaped into the solvent within just a few nanoseconds. Furthermore, an incorrect bifurcated bonding geometry of G-DNA, found in previous MD simulation study under a standard force field but not observed in experiments, disappeared in the present stimulation using the new polarized force field. The current study bridged an important gap between the simulation study and experimental observation on the dynamical structure of G-DNA.

Structural dynamics of human telomeric G-quadruplex loops studied by molecular dynamics simulations

Loops which are linkers connecting G-strands and supporting the G-tetrad core in G-quadruplex are important for biological roles of G-quadruplexes. TTA loop is a common sequence which mainly resides in human telomeric DNA (hTel) G-quadruplex. A series of molecular dynamics (MD) simulations were carried out to investigate the structural dynamics of TTA loops. We found that (1) the TA base pair formed in TTA loops are very stable, the occupied of all hydrogen bonds are more than 0.95. (2) The TA base pair makes the adjacent G-quartet more stable than others. (3) For the edgewise loop and the diagonal loop, most loop bases are stacking with others, only few bases have considerable freedom. (4) The stabilities of these stacking structures are distinct. Part of the loops, especially TA base pairs, and bases stacking with the G-quartet, maintain certain stable conformations in the simulation, but other parts, like TT and TA stacking structures, are not stable enough. For the first time, spontaneous conformational switches of TTA edgewise loops were observed in our long time MD simulations. (5) For double chain reversal loop, it is really hard to maintain a stable conformation in the long time simulation under present force fields (parm99 and parmbsc0), as it has multiple conformations with similar free energies.

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

Explicit solvent molecular dynamics simulations have been used to complement preceding experimental and computational studies of folding of guanine quadruplexes (G-DNA). We initiate early stages of unfolding of several G-DNAs by simulating them under no-salt conditions and then try to fold them back using standard excess salt simulations. There is a significant difference between G-DNAs with all-anti parallel stranded stems and those with stems containing mixtures of syn and anti guanosines. The most natural rearrangement for all-anti stems is a vertical mutual slippage of the strands. This leads to stems with reduced numbers of tetrads during unfolding and a reduction of strand slippage during refolding. The presence of syn nucleotides prevents mutual strand slippage; therefore, the antiparallel and hybrid quadruplexes initiate unfolding via separation of the individual strands. The simulations confirm the capability of G-DNA molecules to adopt numerous stable locally and globally misfolded structures. The key point for a proper individual folding attempt appears to be correct prior distribution of syn and anti nucleotides in all four G-strands. The results suggest that at the level of individual molecules, G-DNA folding is an extremely multi-pathway process that is slowed by numerous misfolding arrangements stabilized on highly variable timescales.

Fragment-based drug discovery using a multidomain, parallel MD-MM/PBSA screening protocol

We have developed a rigorous computational screening protocol to identify novel fragment-like inhibitors of N(5)-CAIR mutase (PurE), a key enzyme involved in de novo purine synthesis that represents a novel target for the design of antibacterial agents. This computational screening protocol utilizes molecular docking, graphics processing unit (GPU)-accelerated molecular dynamics, and Molecular Mechanics/Poisson-Boltzmann Surface Area (MM/PBSA) free energy estimations to investigate the binding modes and energies of fragments in the active sites of PurE. PurE is a functional octamer comprised of identical subunits. The octameric structure, with its eight active sites, provided a distinct advantage in these studies because, for a given simulation length, we were able to place eight separate fragment compounds in the active sites to increase the throughput of the MM/PBSA analysis. To validate this protocol, we have screened an in-house fragment library consisting of 352 compounds. The theoretical results were then compared with the results of two experimental fragment screens, Nuclear Magnetic Resonance (NMR) and Surface Plasmon Resonance (SPR) binding analyses. In these validation studies, the protocol was able to effectively identify the competitive binders that had been independently identified by experimental testing, suggesting the potential utility of this method for the identification of novel fragments for future development as PurE inhibitors.

Conformational dynamics of the human propeller telomeric DNA quadruplex on a microsecond time scale

The human telomeric DNA sequence with four repeats can fold into a parallel-stranded propeller-type topology. NMR structures solved under molecular crowding experiments correlate with the crystal structures found with crystal-packing interactions that are effectively equivalent to molecular crowding. This topology has been used for rationalization of ligand design and occurs experimentally in a number of complexes with a diversity of ligands, at least in the crystalline state. Although G-quartet stems have been well characterized, the interactions of the TTA loop with the G-quartets are much less defined. To better understand the conformational variability and structural dynamics of the propeller-type topology, we performed molecular dynamics simulations in explicit solvent up to 1.5 μs. The analysis provides a detailed atomistic account of the dynamic nature of the TTA loops highlighting their interactions with the G-quartets including formation of an A:A base pair, triad, pentad and hexad. The results present a threshold in quadruplex simulations, with regards to understanding the flexible nature of the sugar-phosphate backbone in formation of unusual architecture within the topology. Furthermore, this study stresses the importance of simulation time in sampling conformational space for this topology.

Guanine base stacking in G-quadruplex nucleic acids

G-quadruplexes constitute a class of nucleic acid structures defined by stacked guanine tetrads (or G-tetrads) with guanine bases from neighboring tetrads stacking with one another within the G-tetrad core. Individual G-quadruplexes can also stack with one another at their G-tetrad interface leading to higher-order structures as observed in telomeric repeat-containing DNA and RNA. In this study, we investigate how guanine base stacking influences the stability of G-quadruplexes and their stacked higher-order structures. A structural survey of the Protein Data Bank is conducted to characterize experimentally observed guanine base stacking geometries within the core of G-quadruplexes and at the interface between stacked G-quadruplex structures. We couple this survey with a systematic computational examination of stacked G-tetrad energy landscapes using quantum mechanical computations. Energy calculations of stacked G-tetrads reveal large energy differences of up to 12 kcal/mol between experimentally observed geometries at the interface of stacked G-quadruplexes. Energy landscapes are also computed using an AMBER molecular mechanics description of stacking energy and are shown to agree quite well with quantum mechanical calculated landscapes. Molecular dynamics simulations provide a structural explanation for the experimentally observed preference of parallel G-quadruplexes to stack in a 5′–5′ manner based on different accessible tetrad stacking modes at the stacking interfaces of 5′–5′ and 3′–3′ stacked G-quadruplexes.

New insights from molecular dynamic simulation studies of the multiple binding modes of a ligand with G-quadruplex DNA

G-quadruplexes are higher-order DNA and RNA structures formed from guanine-rich sequences. These structures have recently emerged as a new class of potential molecular targets for anticancer drugs. An understanding of the three-dimensional interactions between small molecular ligands and their G-quadruplex targets in solution is crucial for rational drug design and the effective optimization of G-quadruplex ligands. Thus far, rational ligand design has been focused mainly on the G-quartet platform. It should be noted that small molecules can also bind to loop nucleotides, as observed in crystallography studies. Hence, it would be interesting to elucidate the mechanism underlying how ligands in distinct binding modes influence the flexibility of G-quadruplex. In the present study, based on a crystal structure analysis, the models of a tetra-substituted naphthalene diimide ligand bound to a telomeric G-quadruplex with different modes were built and simulated with a molecular dynamics simulation method. Based on a series of computational analyses, the structures, dynamics, and interactions of ligand-quadruplex complexes were studied. Our results suggest that the binding of the ligand to the loop is viable in aqueous solutions but dependent on the particular arrangement of the loop. The binding of the ligand to the loop enhances the flexibility of the G-quadruplex, while the binding of the ligand simultaneously to both the quartet and the loop diminishes its flexibility. These results add to our understanding of the effect of a ligand with different binding modes on G-quadruplex flexibility. Such an understanding will aid in the rational design of more selective and effective G-quadruplex binding ligands.

Molecular dynamics simulations to provide new insights into the asymmetrical ammonium ion movement inside of the [d(G3T4G4)]2 G-quadruplex DNA structure

We have used both adaptive biasing force (ABF) and regular molecular dynamics (MD) simulations to investigate the asymmetrical NH(4)(+) ion movement inside of a bimolecular G-quadruplex DNA structure [d(G(3)T(4)G(4))](2). The free-energy landscapes obtained from ABF MD simulations suggest that the NH(4)(+) ion exiting the [d(G(3)T(4)G(4))](2) G-quadruplex stem in the direction toward the edge-type loop (denoted as the upper direction) experiences a lower free-energy barrier than that toward the diagonal loop (denoted as the lower direction) by approximately 3-4 kcal mol(-1). This result is in qualitative agreement with the previous discovery made by Šket and Plavec on the same G-quadruplex structure from (15)N NMR experiments (J. Am. Chem. Soc. 2007, 129, 8794). In the Na(+) form of the same G-quadruplex, Na(+) ion movement was found to be symmetrical, with a free-energy barrier of only 5-7 kcal mol(-1) to cross all three G-quartets, that is, [d(G(3)T(4)G(4))](2) still exhibits ion-channel-like behaviors for Na(+) ions. On the basis of the new computational results, we hypothesize that the stiffness of a G-quartet is primarily determined by the base stacking interactions within the G-quadruplex stem. Therefore, the structural origin for the asymmetrical NH(4)(+) ion movement in [d(G(3)T(4)G(4))](2) is the presence of two different modes of base stacking around the NH(4)(+) binding sites, a more stable 5'-syn-anti mode between lower and central G-quartets and a less stable 5'-anti-anti mode between upper and central G-quartets. Simulations also suggest that loop topology at the end of a G-quadruplex stem only controls the direction at which an exiting NH(4)(+) ion reaches bulk solution but does not impose significant free-energy barriers.

Reference simulations of noncanonical nucleic acids with different χ variants of the AMBER force field: quadruplex DNA, quadruplex RNA and Z-DNA

Refinement of empirical force fields for nucleic acids requires their extensive testing using as wide range of systems as possible. However, finding unambiguous reference data is not easy. In this paper, we analyze four systems which we suggest should be included in standard portfolio of molecules to test nucleic acids force fields, namely, parallel and antiparallel stranded DNA guanine quadruplex stems, RNA quadruplex stem, and Z-DNA. We highlight parameters that should be monitored to assess the force field performance. The work is primarily based on 8.4 μs of 100-250 ns trajectories analyzed in detail followed by 9.6 μs of additional selected back up trajectories that were monitored to verify that the results of the initial analyses are correct. Four versions of the Cornell et al. AMBER force field are tested, including an entirely new parmχ(OL4) variant with χ dihedral specifically reparametrized for DNA molecules containing syn nucleotides. We test also different water models and ion conditions. While improvement for DNA quadruplexes is visible, the force fields still do not fully represent the intricate Z-DNA backbone conformation.

Biocompatible Xanthine-Quadruplex Scaffold for Ion-Transporting DNA Channels

Molecular dynamics simulations and adaptive biasing force analysis of the quadruplex DNA dynamics in an explicit solvent reveal fundamentally different mechanisms of Na(+) transport in xanthine- and guanine-based DNA systems. The barrier to the transport of K(+) through the xanthine-based quadruplex is significantly lower than those reported for the guanine-based analogs.

Molecular dynamics simulations of G-DNA and perspectives on the simulation of nucleic acid structures

The article reviews the application of biomolecular simulation methods to understand the structure, dynamics and interactions of nucleic acids with a focus on explicit solvent molecular dynamics simulations of guanine quadruplex (G-DNA and G-RNA) molecules. While primarily dealing with these exciting and highly relevant four-stranded systems, where recent and past simulations have provided several interesting results and novel insight into G-DNA structure, the review provides some general perspectives on the applicability of the simulation techniques to nucleic acids.

Mass spectrometry and ion mobility spectrometry of G-quadruplexes. A study of solvent effects on dimer formation and structural transitions in the telomeric DNA sequence d(TAGGGTTAGGGT)

We survey here state of the art mass spectrometry methodologies for investigating G-quadruplexes, and will illustrate them with a new study on a simple model system: the dimeric G-quadruplex of the 12-mer telomeric DNA sequence d(TAGGGTTAGGGT), which can adopt either a parallel or an antiparallel structure. We will discuss the solution conditions compatible with electrospray ionisation, the quantification of complexes using ESI-MS, the interpretation of ammonium ion preservation in the complexes in the gas phase, and the use of ion mobility spectrometry to resolve ambiguities regarding the strand stoichiometry, or separate and characterise different structural isomers. We also describe that adding electrospray-compatible organic co-solvents (methanol, ethanol, isopropanol or acetonitrile) to aqueous ammonium acetate increases the stability and rate of formation of dimeric G-quadruplexes, and causes structural transitions to parallel structures. Structural changes were probed by circular dichroism and ion mobility spectrometry, and the excellent correlation between the two techniques validates the use of ion mobility to investigate G-quadruplex folding. We also demonstrate that parallel G-quadruplex structures are easier to preserve in the gas phase than antiparallel structures.