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

Oxidative lesions modulate G-quadruplex stability and structure in the human BCL2 promoter

Misregulation of BCL2 expression has been observed with many diseases and is associated with cellular exposure to reactive oxygen species. A region upstream of the P1 promoter in the human BCL2 gene plays a major role in regulating transcription. This G/C-rich region is highly polymorphic and capable of forming G-quadruplex structures. Herein we report that an oxidative event simulated with an 8-oxo-7,8-dihydroguanine (^oxoG) substitution within a long G-tract results in a reduction of structural polymorphism. Surprisingly, ^oxoG within a 25-nt construct boosts thermal stability of the resulting G-quadruplex. This is achieved by distinct hydrogen bonding properties of ^oxoG, which facilitate formation of an antiparallel basket-type G-quadruplex with a three G-quartet core and a G·^oxoG·C base triad. While ^oxoG has previously been considered detrimental for G-quadruplex formation, its stabilizing effect within a promoter described in this study suggests a potential novel regulatory role of oxidative stress in general and specifically in BCL2 gene transcription.

Highlighting the potential utility of MBP crystallization chaperone for Arabidopsis BIL1/BZR1 transcription factor-DNA complex

The maltose-binding protein (MBP) fusion tag is one of the most commonly utilized crystallization chaperones for proteins of interest. Recently, this MBP-mediated crystallization technique was adapted to Arabidopsis thaliana (At) BRZ-INSENSITIVE-LONG (BIL1)/BRASSINAZOLE-RESISTANT (BZR1), a member of the plant-specific BZR TFs, and revealed the first structure of AtBIL1/BZR1 in complex with target DNA. However, it is unclear how the fused MBP affects the structural features of the AtBIL1/BZR1-DNA complex. In the present study, we highlight the potential utility of the MBP crystallization chaperone by comparing it with the crystallization of unfused AtBIL1/BZR1 in complex with DNA. Furthermore, we assessed the validity of the MBP-fused AtBIL1/BZR1-DNA structure by performing detailed dissection of crystal packings and molecular dynamics (MD) simulations with the removal of the MBP chaperone. Our MD simulations define the structural basis underlying the AtBIL1/BZR1-DNA assembly and DNA binding specificity by AtBIL1/BZR1. The methodology employed in this study, the combination of MBP-mediated crystallization and MD simulation, demonstrates promising capabilities in deciphering the protein-DNA recognition code.

Base-pair resolution analysis of the effect of supercoiling on DNA flexibility and major groove recognition by triplex-forming oligonucleotides

In the cell, DNA is arranged into highly-organised and topologically-constrained (supercoiled) structures. It remains unclear how this supercoiling affects the detailed double-helical structure of DNA, largely because of limitations in spatial resolution of the available biophysical tools. Here, we overcome these limitations, by a combination of atomic force microscopy (AFM) and atomistic molecular dynamics (MD) simulations, to resolve structures of negatively-supercoiled DNA minicircles at base-pair resolution. We observe that negative superhelical stress induces local variation in the canonical B-form DNA structure by introducing kinks and defects that affect global minicircle structure and flexibility. We probe how these local and global conformational changes affect DNA interactions through the binding of triplex-forming oligonucleotides to DNA minicircles. We show that the energetics of triplex formation is governed by a delicate balance between electrostatics and bonding interactions. Our results provide mechanistic insight into how DNA supercoiling can affect molecular recognition, that may have broader implications for DNA interactions with other molecular species.

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.

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.

Phosphorothioate Substitutions in RNA Structure Studied by Molecular Dynamics Simulations, QM/MM Calculations, and NMR Experiments

Phosphorothioates (PTs) are important chemical modifications of the RNA backbone where a single nonbridging oxygen of the phosphate is replaced with a sulfur atom. PT can stabilize RNAs by protecting them from hydrolysis and is commonly used as a tool to explore their function. It is, however, unclear what basic physical effects PT has on RNA stability and electronic structure. Here, we present molecular dynamics (MD) simulations, quantum mechanical (QM) calculations, and NMR spectroscopy measurements, exploring the effects of PT modifications in the structural context of the neomycin-sensing riboswitch (NSR). The NSR is the smallest biologically functional riboswitch with a well-defined structure stabilized by a U-turn motif. Three of the signature interactions of the U-turn: an H-bond, an anion-π interaction, and a potassium binding site; are formed by RNA phosphates, making the NSR an ideal model for studying how PT affects RNA structure and dynamics. By comparing with high-level QM calculations, we reveal the distinct physical properties of the individual interactions facilitated by the PT. The sulfur substitution, besides weakening the direct H-bond interaction, reduces the directionality of H-bonding while increasing its dispersion and induction components. It also reduces the induction and increases the dispersion component of the anion-π stacking. The sulfur force-field parameters commonly employed in the literature do not reflect these distinctions, leading to the unsatisfactory description of PT in simulations of the NSR. We show that it is not possible to accurately describe the PT interactions using one universal set of van der Waals sulfur parameters and provide suggestions for improving the force-field performance.

WASP: a software package for correctly characterizing the topological development of ribbon structures

We introduce the Writhe Application Software Package (WASP) which can be used to characterisze the topology of ribbon structures, the underlying mathematical model of DNA, Biopolymers, superfluid vorticies, elastic ropes and magnetic flux ropes. This characterization is achieved by the general twist–writhe decomposition of both open and closed ribbons, in particular through a quantity termed the polar writhe. We demonstrate how this decomposition is far more natural and straightforward than artificial closure methods commonly utilized in DNA modelling. In particular, we demonstrate how the decomposition of the polar writhe into local and non-local components distinctly characterizes the local helical structure and knotting/linking of the ribbon. This decomposition provides additional information not given by alternative approaches. As example applications, the WASP routines are used to characterise the evolving topology (writhe) of DNA minicircle and open ended plectoneme formation magnetic/optical tweezer simulations, and it is shown that the decomponsition into local and non-local components is particularly important for the detection of plectonemes. Finally it is demonstrated that a number of well known alternative writhe expressions are actually simplifications of the polar writhe measure.

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.

Locked nucleic acid building blocks as versatile tools for advanced G-quadruplex design

A hybrid-type G-quadruplex is modified with LNA (locked nucleic acid) and 2′-F-riboguanosine in various combinations at the two syn positions of its third antiparallel G-tract. LNA substitution in the central tetrad causes a complete rearrangement to either a V-loop or antiparallel structure, depending on further modifications at the 5′-neighboring site. In the two distinct structural contexts, LNA-induced stabilization is most effective compared to modifications with other G surrogates, highlighting a potential use of LNA residues for designing not only parallel but various more complex G4 structures. For instance, the conventional V-loop is a structural element strongly favored by an LNA modification at the V-loop 3′-end in contrast with an alternative V-loop, clearly distinguishable by altered conformational properties and base-backbone interactions as shown in a detailed analysis of V-loop structures.

Structure of a DNA G-Quadruplex Related to Osteoporosis with a G-A Bulge Forming a Pseudo-loop

Bone remodeling is a fine-tuned process principally regulated by a cascade triggered by interaction of receptor activator of NF-κB (RANK) and RANK ligand (RANKL). Excessive activity of the RANKL gene leads to increased bone resorption and can influence the incidence of osteoporosis. Although much has been learned about the intracellular signals activated by RANKL/RANK complex, significantly less is known about the molecular mechanisms of regulation of RANKL expression. Here, we report on the structure of an unprecedented DNA G-quadruplex, well-known secondary structure-mediated gene expression regulator, formed by a G-rich sequence found in the regulatory region of a RANKL gene. Solution-state NMR structural study reveals the formation of a three-layered parallel-type G-quadruplex characterized by an unique features, including a G-A bulge. Although a guanine within a G-tract occupies syn glycosidic conformation, bulge-forming residues arrange in a pseudo-loop conformation to facilitate partial 5/6-ring stacking, typical of G-quadruplex structures with parallel G-tracts orientation. Such distinctive structural features protruding from the core of the structure can represent a novel platform for design of highly specific ligands with anti-osteoporotic function. Additionally, our study suggests that the expression of RANKL gene may be regulated by putative folding of its G-rich region into non-B-DNA structure(s).

Molecular Dynamics Study on the Binding of an Anticancer DNA G-Quadruplex Stabilizer, CX-5461, to Human Telomeric, c-KIT1, and c-Myc G-Quadruplexes and a DNA Duplex

DNA G-quadruplex (G4) stabilizer, CX-5461, is in phase I/II clinical trials for advanced cancers with BRCA1/2 deficiencies. A FRET-melting temperature increase assay measured the stabilizing effects of CX-5461 to a DNA duplex (∼10 K), and three G4 forming sequences negatively implicated in the cancers upon its binding: human telomeric (∼30 K), c-KIT1 (∼27 K), and c-Myc (∼25 K). Without experimentally solved structures of these CX-5461-G4 complexes, CX-5461's interactions remain elusive. In this study, we performed a total of 73.5 μs free ligand molecular dynamics binding simulations of CX-5461 to the DNA duplex and three G4s. Three binding modes (top, bottom, and side) were identified for each system and their thermodynamic, kinetic, and structural nature were deciphered. The molecular mechanics/Poisson Boltzmann surface area binding energies of CX-5461 were calculated for the human telomeric (-28.6 kcal/mol), c-KIT1 (-23.9 kcal/mol), c-Myc (-22.0 kcal/mol) G4s, and DNA duplex (-15.0 kcal/mol) systems. These energetic differences coupled with structural differences at the 3' site explained the different melting temperatures between the G4s, while CX-5461's lack of intercalation to the duplex explained the difference between the G4s and duplex. Based on the interaction insight, CX-5461 derivatives were designed and docked, showing higher selectivity to the G4s over the duplex.

Transcription-facilitating histone chaperons interact with genomic and synthetic G4 structures

Affinity for G-quadruplex (G4) structures may be a common feature of transcription-facilitating histone chaperons (HCs). This assumption is based on previous unmatched studies of HCs FACT, nucleolin (NCL), BRD3, and ATRX. We verified this assumption and considered its implications for the therapeutic applications of synthetic (exogenous) G4s and the biological significance of genomic G4s. First, we questioned whether exogenous G4s that recognize cell-surface NCL and could trap other HCs in the nucleus are usable as anticancer agents. We performed in vitro binding assays and selected leading multi-targeted G4s. They exhibited minor effects on cell viability. The presumed NCL-regulated intracellular transport of G4s was inefficient or insufficient for tumor-specific G4 delivery. Next, to clarify whether G4s in the human genome could recruit HCs, we compared available HC ChIP-seq data with G4-seq/G4-ChIP-seq data. Several G4s, including the well-known c-Myc quadruplex structure, were found to be colocalized with HC occupancy sites in cancer cell lines. As evidenced by our molecular modeling data, c-Myc G4 might interfere with the HC function of BRD3 but is unlikely to prevent the BRD3-driven assembly of the chromatin remodeling complex. The c-Myc case illustrates the intricate role of genomic G4s in chromatin remodeling, nucleosome remodeling, and transcription.

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

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

Nucleotide Loading Modes of Human RNA Polymerase II as Deciphered by Molecular Simulations

Mapping the route of nucleoside triphosphate (NTP) entry into the sequestered active site of RNA polymerase (RNAP) has major implications for elucidating the complete nucleotide addition cycle. Constituting a dichotomy that remains to be resolved, two alternatives, direct NTP delivery via the secondary channel (CH2) or selection to downstream sites in the main channel (CH1) prior to catalysis, have been proposed. In this study, accelerated molecular dynamics simulations of freely diffusing NTPs about RNAPII were applied to refine the CH2 model and uncover atomic details on the CH1 model that previously lacked a persuasive structural framework to illustrate its mechanism of action. Diffusion and binding of NTPs to downstream DNA, and the transfer of a preselected NTP to the active site, are simulated for the first time. All-atom simulations further support that CH1 loading is transcription factor IIF (TFIIF) dependent and impacts catalytic isomerization. Altogether, the alternative nucleotide loading systems may allow distinct transcriptional landscapes to be expressed.

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.

GC ends control topology of DNA G-quadruplexes and their cation-dependent assembly

GCn and GCnCG, where n = (G2AG4AG2), fold into well-defined, dimeric G-quadruplexes with unprecedented folding topologies in the presence of Na+ ions as revealed by nuclear magnetic resonance spectroscopy. Both G-quadruplexes exhibit unique combination of structural elements among which are two G-quartets, A(GGGG)A hexad and GCGC-quartet. Detailed structural characterization uncovered the crucial role of 5'-GC ends in formation of GCn and GCnCG G-quadruplexes. Folding in the presence of 15NH4+ and K+ ions leads to 3'-3' stacking of terminal G-quartets of GCn G-quadruplexes, while 3'-GC overhangs in GCnCG prevent dimerization. Results of the present study expand repertoire of possible G-quadruplex structures. This knowledge will be useful in DNA sequence design for nanotechnological applications that may require specific folding topology and multimerization properties.

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

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

PDGFR-β Promoter Forms a Vacancy G-Quadruplex that Can Be Filled in by dGMP: Solution Structure and Molecular Recognition of Guanine Metabolites and Drugs

Aberrant expression of PDGFR-β is associated with a number of diseases. The G-quadruplexes (G4s) formed in PDGFR-β gene promoter are transcriptional modulators and amenable to small molecule targeting. The major G4 formed in the PDGFR-β gene promoter was previously shown to have a broken G-strand. Herein, we report that the PDGFR-β gene promoter sequence forms a vacancy G-quadruplex (vG4) which can be filled in and stabilized by physiologically relevant guanine metabolites, such as dGMP, GMP, and cGMP, as well as guanine-derivative drugs. We determined the NMR structure of the dGMP-fill-in PDGFR-β vG4 in K⁺ solution. This is the first structure of a guanine-metabolite-fill-in vG4 based on a human gene promoter sequence. Our structure and systematic analysis elucidate the contributions of Hoogsten hydrogen bonds, sugar, and phosphate moieties to the specific G-vacancy fill-in. Intriguingly, an equilibrium of 3'- and 5'-end vG4s is present in the PDGFR-β promoter sequence, and dGMP favors the 5'-end fill-in. Guanine metabolites and drugs were tested and showed a conserved selectivity for the 5'-vacancy, except for cGMP. cGMP binds both the 3'- and 5'-end vG4s and forms two fill-in G4s with similar population. Significantly, guanine metabolites are involved in many physiological and pathological processes in human cells; thus, our results provide a structural basis to understand their potential regulatory functions by interaction with promoter vG4s. Moreover, the NMR structure can guide rational design of ligands that target the PDGFR-β vG4.

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.

Truncated tetrahedral RNA nanostructures exhibit enhanced features for delivery of RNAi substrates

Using RNA as a material for nanoparticle construction provides control over particle size and shape at the nano-scale. RNA nano-architectures have shown promise as delivery vehicles for RNA interference (RNAi) substrates, allowing multiple functional entities to be combined on a single particle in a programmable fashion. Rather than employing a completely bottom-up approach to scaffold design, here multiple copies of an existing synthetic supramolecular RNA nano-architecture serve as building blocks along with additional motifs for the design of a novel truncated tetrahedral RNA scaffold, demonstrating that rationally designed RNA assemblies can themselves serve as modular pieces in the construction of larger rationally designed structures. The resulting tetrahedral scaffold displays enhanced characteristics for RNAi-substrate delivery in comparison to similar RNA-based scaffolds, as evidenced by its increased functional capacity, increased cellular uptake and ultimately an increased RNAi efficacy of its adorned Dicer substrate siRNAs. The unique truncated tetrahedral shape of the nanoparticle core appears to contribute to this particle's enhanced function, indicating the physical characteristics of RNA scaffolds merit significant consideration when designing platforms for delivery of functional RNAs via RNA nanoparticles.

Cucurbit[7]uril-driven modulation of ligand-DNA interactions by ternary assembly

The design of small organic molecules with a predictable and desirable DNA-binding mechanism is a topical research task for biomedicine application. Herein, we demonstrate an attractive supramolecular strategy for controlling the non-covalent ligand-DNA interaction by binding with cucurbituril as a synthetic receptor. With a combination of UV/vis, CD and NMR experiments, we demonstrate that the bis-styryl dye with two suitable binding sites can involve double stranded DNA and cucurbituril in the formation of the supramolecular triad. The ternary assembly is formed as a result of the interaction of macrocyclic cucurbituril with one pyridinium fragment of the bis-styryl dye, while the second pyridinium fragment of the dye is effectively associated with DNA backbones, which leads to a change in the ligand-DNA binding mode from aggregation to a minor groove. This exciting outcome was supported by molecular docking studies that help to understand the molecular orientation of the supramolecular triad and elucidate the destruction of dye aggregates caused by cucurbituril. These studies provide valuable information on the mechanisms of DNA binding to small molecules and recognition processes in bioorganic supramolecular assemblies constructed from multiple non-covalent interactions.

High-resolution EPR distance measurements on RNA and DNA with the non-covalent Ǵ spin label

Pulsed electron paramagnetic resonance (EPR) experiments, among them most prominently pulsed electron-electron double resonance experiments (PELDOR/DEER), resolve the conformational dynamics of nucleic acids with high resolution. The wide application of these powerful experiments is limited by the synthetic complexity of some of the best-performing spin labels. The recently developed $\bf\acute{G}$ (G-spin) label, an isoindoline-nitroxide derivative of guanine, can be incorporated non-covalently into DNA and RNA duplexes via Watson-Crick base pairing in an abasic site. We used PELDOR and molecular dynamics (MD) simulations to characterize $\bf\acute{G}$, obtaining excellent agreement between experiments and time traces calculated from MD simulations of RNA and DNA double helices with explicitly modeled $\bf\acute{G}$ bound in two abasic sites. The MD simulations reveal stable hydrogen bonds between the spin labels and the paired cytosines. The abasic sites do not significantly perturb the helical structure. $\bf\acute{G}$ remains rigidly bound to helical RNA and DNA. The distance distributions between the two bound $\bf\acute{G}$ labels are not substantially broadened by spin-label motions in the abasic site and agree well between experiment and MD. $\bf\acute{G}$ and similar non-covalently attached spin labels promise high-quality distance and orientation information, also of complexes of nucleic acids and proteins.

Sugar Puckering Drives G-Quadruplex Refolding: Implications for V-Shaped Loops

A DNA G-quadruplex adopting a (3+1) hybrid structure was modified in two adjacent syn positions of the antiparallel strand with anti-favoring 2'-deoxy-2'-fluoro-riboguanosine (F rG) analogues. The two substitutions promoted a structural rearrangement to a topology with the 5'-terminal G residue located in the central tetrad and the two modified residues linked by a V-shaped zero-nucleotide loop. Strikingly, whereas a sugar pucker in the preferred north domain is found for both modified nucleotides, the F rG analogue preceding the V-loop is forced to adopt the unfavored syn conformation in the new quadruplex fold. Apparently, a preferred C3'-endo sugar pucker within the V-loop architecture outweighs the propensity of the F rG analogue to adopt an anti glycosidic conformation. Refolding into a V-loop topology is likewise observed for a sequence modified at corresponding positions with two riboguanosine substitutions. In contrast, 2'-F-arabinoguanosine analogues with their favored south-east sugar conformation do not support formation of the V-loop topology. Examination of known G-quadruplexes with a V-shaped loop highlights the critical role of the sugar conformation for this distinct structural motif.

Three-Dimensional Structure of RNA Monomeric G-Quadruplex Containing ALS and FTD Related G4C2 Repeat and Its Binding with TMPyP4 Probed by Homology Modeling based on Experimental Constraints and Molecular Dynamics Simulations

The G-quadruplex-forming hexanucleotide repeat expansion (HRE), d(G4C2)_n, within the human C9orf72 gene is the root cause for familial amyotrophic lateral sclerosis-frontotemporal dementia (ALS-FTD). A recent study has shown that TMPyP4 has good potential to work as a RNA G-quadruplex binder in treating ALS and FTD. Although the high-resolution structure of the monomeric DNA antiparallel G-quadruplex form of the monomeric hexanucleotide repeat was recently solved, the RNA parallel G-quadruplex structure and its complex with TMPyP4 are not available yet. In this study, we first constructed the homology model for the parallel monomeric RNA G-quadruplex of r(G4C2)₃G4 based on experimental constraints and the parallel monomeric G-quadruplex DNA crystal structure. Although the G-tetra core of the homology model was stable observed in 15 μs molecular dynamics (MD) simulations, we observed that the loops adopt additional conformations besides the initial crystal conformation, where TMPyP4 binding was found to reduce the loop fluctuation of the RNA monomeric G-quadruplex. Next, we probed the elusive binding behavior of TMPyP4 to the RNA monomeric G-quadruplex. Encouragingly, the binding modes observed are similar to the modes observed in two experimental complexes of a parallel DNA G-quadruplex with TMPyP4. We also constructed a Markov state model to provide insights into the binding pathways. Together, the findings from our study may assist future development of G-quadruplex-specific ligands in the treatment of neurodegenerative diseases like ALS and FTD.

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.

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.

Pursuing origins of (poly)ethylene glycol-induced G-quadruplex structural modulations

Molecular crowding conditions provided by high concentration of cosolutes are utilized for characterization of biomolecules in cell-mimicking environment and development of drug-delivery systems. In this context, (poly)ethylene glycols are often used for studying non-canonical DNA structures termed G-quadruplexes, which came into focus by emerging structural biology findings and new therapeutic drug design approaches. Recently, several reports were made arguing against using (poly)ethylene glycols in role of molecular crowding agents due to their direct impact on DNA G-quadruplex stability and topology. However, the available data on structural details underlying DNA interaction is very scarce and thus limits in-depth comprehension. Herein, structural and thermodynamic analyses were strategically combined to assess G-quadruplex-cosolute interactions and address previously reported variances regarding the driving forces of G-rich DNA structural transformations under molecular crowding conditions. With the use of complementary (CD, NMR and UV) spectroscopic methods and model approach we characterized DNA G-quadruplex in the presence of the smallest and one of the largest typically used (poly)ethylene glycols. Dehydration effect is the key contributor to ethylene-glycol-induced increased stability of the G-quadruplex, which is in the case of the large cosolute mainly guided by the subtle direct interactions between PEG 8000 and the outer G-quartet regions.

Redox-dependent control of i-Motif DNA structure using copper cations

Previous computational studies have shown that Cu⁺ can act as a substitute for H⁺ to support formation of cytosine (C) dimers with similar conformation to the hemi-protonated base pair found in i-motif DNA. Through a range of biophysical methods, we provide experimental evidence to support the hypothesis that Cu⁺ can mediate C–C base pairing in i-motif DNA and preserve i-motif structure. These effects can be reversed using a metal chelator, or exposure to ambient oxygen in the air that drives oxidation of Cu⁺ to Cu²⁺, a comparatively weak ligand. Herein, we present a dynamic and redox-sensitive system for conformational control of an i-motif forming DNA sequence in response to copper cations.

Principles of chemical geometry underlying chiral selectivity in RNA minihelix aminoacylation

The origin of homochirality in L-amino acid in proteins is one of the mysteries of the evolution of life. Experimental studies show that a non-enzymatic aminoacylation reaction of an RNA minihelix has a preference for L-amino acid over D-amino acid. The reaction initiates by approaching of a 3'-oxygen of the RNA minihelix to the carbonyl carbon of an aminoacyl phosphate oligonucleotide. Here, employing molecular dynamics simulations, we examined the possible mechanisms that determine this chiral selectivity. The simulation system adopted a geometry required for the chemical reaction to occur more frequently with L-alanine than that with D-alanine. For L-alanine, the structure with this geometry was formed by a combination of stable dihedral angles along alanyl phosphate backbone with a canonical RNA structure, where the methyl group of alanine was placed on the opposite side of the approaching 3'-hydroxyl group with respect to the carbonyl plane. For D-alanine, the methyl group and the 3'-hydroxyl group were placed on the same side with respect to the carbonyl plane, which significantly decreased its ability to approach 3'-oxygen close to the carbonyl carbon compared to L-alanine. The mechanism suggested herein can explain experimentally observed chiral preferences.

Sulfur-substitution-induced base flipping in the DNA duplex

Base flipping is widely observed in a number of important biological processes. The genetic codes deposited inside the DNA duplex become accessible to external agents upon base flipping. The sulfur substitution of guanine leads to thioguanine, which alters the thermodynamic stability of the GC base pairs and the GT mismatches. Experimental studies conclude that the sulfur substitution decreases the lifetime of the GC base pair. In this work, under three AMBER force fields for nucleotide systems, we firstly performed equilibrium and nonequilibrium free energy simulations to investigate the variation of the thermodynamic profiles in base flipping upon sulfur substitution. It is found that the bsc0 modification, the bsc1 modification and the OL15 modification of AMBER force fields are able to qualitatively describe the sulfur-substitution dependent behavior of the thermodynamics. However, only the two last-generation AMBER force fields are able to provide quantitatively correct predictions. The second computational study on the sulfur substitutions focused on the relative stability of the S6G-C base pair and the S6G-T mismatch. Two conflicting experimental observations were reported by the same authors. One suggested that the S6G-C base pair was more stable, while the other concludes that the S6G-T mismatch was more stable. We answered this question by constructing the free energy profiles along the base flipping pathway computationally.

NMR structure of a G-quadruplex formed by four d(G4C2) repeats: insights into structural polymorphism

Most frequent genetic cause of amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD), is a largely increased number of d(G4C2)n•(G2C4)n repeats located in the non-coding region of C9orf72 gene. Non-canonical structures, including G-quadruplexes, formed within expanded repeats have been proposed to drive repeat expansion and pathogenesis of ALS and FTD. Oligonucleotide d[(G4C2)3G4], which represents the shortest oligonucleotide model of d(G4C2) repeats with the ability to form a unimolecular G-quadruplex, forms two major G-quadruplex structures in addition to several minor species which coexist in solution with K+ ions. Herein, we used solution-state NMR to determine the high-resolution structure of one of the major G-quadruplex species adopted by d[(G4C2)3G4]. Structural characterization of the G-quadruplex named AQU was facilitated by a single substitution of dG with 8Br-dG at position 21 and revealed an antiparallel fold composed of four G-quartets and three lateral C-C loops. The G-quadruplex exhibits high thermal stability and is favored kinetically and under slightly acidic conditions. An unusual structural element distinct from a C-quartet is observed in the structure. Two C•C base pairs are stacked on the nearby G-quartet and are involved in a dynamic equilibrium between symmetric N3-amino and carbonyl-amino geometries and protonated C+•C state.

Disease-associated missense variants in ZBTB18 disrupt DNA binding and impair the development of neurons within the embryonic cerebral cortex

The activities of DNA-binding transcription factors, such as the multi-zinc-finger protein ZBTB18 (also known as RP58, or ZNF238), are essential to coordinate mammalian neurodevelopment, including the birth and radial migration of newborn neurons within the fetal brain. In humans, the majority of disease-associated missense mutations in ZBTB18 lie within the DNA-binding zinc-finger domain and are associated with brain developmental disorder, yet the molecular mechanisms explaining their role in disease remain unclear. To address this, we developed in silico models of ZBTB18, bound to DNA, and discovered that half of the missense variants map to residues (Asn461, Arg464, Glu486) predicted to be essential to sequence-specific DNA contact, whereas others map to residues (Leu434, Tyr447, Arg495) with limited contributions to DNA binding. We studied pathogenic variants to residues with close (N461S) and limited (R495G) DNA contact and found that each bound DNA promiscuously, displayed altered transcriptional regulatory activity in vitro, and influenced the radial migration of newborn neurons in vivo in different ways. Taken together, our results suggest that altered transcriptional regulation could represent an important pathological mechanism for ZBTB18 missense variants in brain developmental disease.

Determination of Base-Flipping Free-Energy Landscapes from Nonequilibrium Stratification

Correct calculation of the variation of free energy upon base flipping is crucial in understanding the dynamics of DNA systems. The free-energy landscape along the flipping pathway gives the thermodynamic stability and the flexibility of base-paired states. Although numerous free-energy simulations are performed in the base flipping cases, no theoretically rigorous nonequilibrium techniques are devised and employed to investigate the thermodynamics of base flipping. In the current work, we report a general nonequilibrium stratification scheme for the efficient calculation of the free-energy landscape of base flipping in DNA duplex. We carefully monitor the convergence behavior of the equilibrium sampling based free-energy simulation and the nonequilibrium stratification and determine the empirical length of time blocks required for converged sampling. Comparison between the performances of the equilibrium umbrella sampling and the nonequilibrium stratification is given. The results show that nonequilibrium free-energy simulation achieves similar accuracy and efficiency compared with the equilibrium enhanced sampling technique in the base flipping cases. We further test a convergence criterion we previously proposed and it comes out that the convergence determined by this criterion agrees with those given by the time-invariant behavior of PMF and the nonlinear dependence of standard deviation on the sample size.

Benzothiazole-based cyanines as fluorescent "light-up" probes for duplex and quadruplex DNA

Analogs of benzothiazole orange (BO) with one, two or three methylbenzothiazolylmethylidene substituents in the 1-methylpyridinium ring were obtained from the respective picolinium, lutidinium or collidinium salts. Fluorescence parameters of the known and new dyes in complexes with various DNA structures, including G-quadruplexes (G4s) and i-motifs (IMs), were analyzed. All dyes efficiently distinguished G4s and ss-DNA. The bi- and tri-substituted derivatives had basically similar distributions of relative fluorescence intensities. The mono-substituted derivatives exhibited enhanced sensitivity to parallel G4s. All dyes were particularly sensitive to a G4 structure with an additional duplex module (the thrombin-binding aptamer TBA31), presumably due to a distinctive binding mode (interaction with the junction between the two modules). In particular, BO showed a strong (160-fold) enhancement in fluorescence quantum yield in complex with TBA31 compared to the free dye. The fluorescence quantum yields of the 2,4-bisubstituted derivative in complex with well-characterized G4s from oncogene promoters were in the range of 0.04-0.28, i.e. comparable to those of ThT. The mono/bi-substituted derivatives should be considered as possible light-up probes for G4 formation.

Nanomechanical DNA resonators for sensing and structural analysis of DNA-ligand complexes

The effect of direct or indirect binding of intercalant molecules on DNA structure is of fundamental importance in understanding the biological functioning of DNA. Here we report on self-suspended DNA nanobundles as ultrasensitive nanomechanical resonators for structural studies of DNA-ligand complexes. Such vibrating nanostructures represent the smallest mechanical resonator entirely composed of DNA. A correlative analysis between the mechanical and structural properties is exploited to study the intrinsic changes of double strand DNA, when interacting with different intercalant molecules (YOYO-1 and GelRed) and a chemotherapeutic drug (Cisplatin), at different concentrations. Possible implications of our findings are related to the study of interaction mechanism of a wide category of molecules with DNA, and to further applications in medicine, such as optimal titration of chemotherapeutic drugs and environmental studies for the detection of heavy metals in human serum.

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.

Towards Understanding of Polymorphism of the G-rich Region of Human Papillomavirus Type 52

The potential to affect gene expression via G-quadruplex stabilization has been extended to all domains of life, including viruses. Here, we investigate the polymorphism and structures of G-quadruplexes of the human papillomavirus type 52 with UV, CD and NMR spectroscopy and gel electrophoresis. We show that oligonucleotide with five G-tracts folds into several structures and that naturally occurring single nucleotide polymorphisms (SNPs) have profound effects on the structural polymorphism in the context of G-quadruplex forming propensity, conformational heterogeneity and folding stability. With help of SNP analysis, we were able to select one of the predominant forms, formed by G-rich sequence d(G₃TAG₃CAG₄ACACAG₃T). This oligonucleotide termed HPV52_(1-4) adopts a three G-quartet snap back (3 + 1) type scaffold with four syn guanine residues, two edgewise loops spanning the same groove, a no-residue V loop and a propeller type loop. The first guanine residue is incorporated in the central G-quartet and all four-guanine residues from G4 stretch are included in the three quartet G-quadruplex core. Modification studies identified several structural elements that are important for stabilization of the described G-quadruplex fold. Our results expand set of G-rich targets in viral genomes and address the fundamental questions regarding folding of G-rich sequences.

Selective recognition of c-MYC Pu22 G-quadruplex by a fluorescent probe

Nucleic acid mimics of fluorescent proteins can be valuable tools to locate and image functional biomolecules in cells. Stacking between the internal G-quartet, formed in the mimics, and the exogenous fluorophore probes constitutes the basis for fluorescence emission. The precision of recognition depends upon probes selectively targeting the specific G-quadruplex in the mimics. However, the design of probes recognizing a G-quadruplex with high selectivity in vitro and in vivo remains a challenge. Through structure-based screening and optimization, we identified a light-up fluorescent probe, 9CI that selectively recognizes c-MYC Pu22 G-quadruplex both in vitro and ex vivo. Upon binding, the biocompatible probe emits both blue and green fluorescence with the excitation at 405 nm. With 9CI and c-MYC Pu22 G-quadruplex complex as the fluorescent response core, a DNA mimic of fluorescent proteins was constructed, which succeeded in locating a functional aptamer on the cellular periphery. The recognition mechanism analysis suggested the high selectivity and strong fluorescence response was attributed to the entire recognition process consisting of the kinetic match, dynamic interaction, and the final stacking. This study implies both the single stacking state and the dynamic recognition process are crucial for designing fluorescent probes or ligands with high selectivity for a specific G-quadruplex structure.

Binding of BRACO19 to a Telomeric G-Quadruplex DNA Probed by All-Atom Molecular Dynamics Simulations with Explicit Solvent

Although BRACO19 is a potent G-quadruplex binder, its potential for clinical usage is hindered by its low selectivity towards DNA G-quadruplex over duplex. High-resolution structures of BRACO19 in complex with neither single-stranded telomeric DNA G-quadruplexes nor B-DNA duplex are available. In this study, the binding pathway of BRACO19 was probed by 27.5 µs molecular dynamics binding simulations with a free ligand (BRACO19) to a DNA duplex and three different topological folds of the human telomeric DNA G-quadruplex (parallel, anti-parallel and hybrid). The most stable binding modes were identified as end stacking and groove binding for the DNA G-quadruplexes and duplex, respectively. Among the three G-quadruplex topologies, the MM-GBSA binding energy analysis suggested that BRACO19′s binding to the parallel scaffold was most energetically favorable. The two lines of conflicting evidence plus our binding energy data suggest conformation-selection mechanism: the relative population shift of three scaffolds upon BRACO19 binding (i.e., an increase of population of parallel scaffold, a decrease of populations of antiparallel and/or hybrid scaffold). This hypothesis appears to be consistent with the fact that BRACO19 was specifically designed based on the structural requirements of the parallel scaffold and has since proven effective against a variety of cancer cell lines as well as toward a number of scaffolds. In addition, this binding mode is only slightly more favorable than BRACO19s binding to the duplex, explaining the low binding selectivity of BRACO19 to G-quadruplexes over duplex DNA. Our detailed analysis suggests that BRACO19′s groove binding mode may not be stable enough to maintain a prolonged binding event and that the groove binding mode may function as an intermediate state preceding a more energetically favorable end stacking pose; base flipping played an important role in enhancing binding interactions, an integral feature of an induced fit binding mechanism.

Solution Structure of a MYC Promoter G-Quadruplex with 1:6:1 Loop Length

The important MYC oncogene is deregulated in many cancer cells and comprises one of the most prominent G-quadruplex (G4) forming sequences in its promoter regions, the NHE III₁ motif. Formation of G4s suppresses MYC transcription and can be modulated by drug binding, establishing these DNA structures as promising targets in cancer therapy. The NHE III₁ motif can fold into more than one parallel G4s, including 1:2:1 and 1:6:1 loop length conformers, with the 1:2:1 conformer shown as the major species under physiological conditions in solution. However, additional factors such as protein interactions may affect the cellular folding equilibrium. Nucleolin, a protein shown to bind MYC G4 and repress MYC transcription, is reported herein to preferably bind to the 1:6:1 loop length conformer suggesting a physiological significance of this species. The high-resolution NMR solution structure of the 1:6:1 conformer is determined, which reveals a 5'-capping structure distinctive from the 1:2:1 form, with the 6 nt central loop playing an essential role for this specific capping structure. This suggests that each parallel G-quadruplex likely adopts unique capping and loop structures determined by the specific central loop and flanking sequences. The resulting structural information at the molecular level will help to understand protein recognition of different G4s, contribution of G4 polymorphism to gene regulation, and to rationally design small molecules selectively targeting the 1:6:1 MYC G4.

Molecular Dynamics Simulations of 2-Aminopurine-Labeled Dinucleoside Monophosphates Reveal Multiscale Stacking Kinetics

Molecular dynamics (MD) simulations of 2-aminopurine (2Ap)-labeled DNA dinucleoside monophosphates (DNMPs) were performed to investigate the hypothesis that base stacking dynamics occur on timescales sufficiently rapid to influence the emission signals measured in time-resolved fluorescence experiments. Analysis of multiple microsecond-length trajectories shows that the DNMPs sample all four coplanar stacking motifs. In addition, three metastable unstacked conformations are detected. A hidden Markov-state model (HMSM) was applied to the simulations to estimate transition rates between the stacked and unstacked states. Transitions between different stacked states generally occur at higher rates when the number of nucleobase faces requiring desolvation is minimized. Time constants for structural relaxation range between 1.6 and 25 ns, suggesting that emission from photoexcited 2Ap, which has an excited-state lifetime of 10 ns, is sensitive to base stacking kinetics. A master equation model for the excited-state population of 2Ap predicts multiexponential emission decays that reproduce the sub-10 ns emission decay lifetimes and amplitudes seen in experiments. Combining MD simulations with HMSM analysis is a powerful way to understand the dynamics that influence 2Ap excited-state relaxation and represents an important step toward using observed emission signals to validate MD simulations.

Impact of Oxidative Lesions on the Human Telomeric G-Quadruplex

Telomere attrition is closely associated with cell aging and exposure to reactive oxygen species (ROS). While oxidation products of nucleotides have been studied extensively in the past, the underlying secondary/tertiary structural changes in DNA remain poorly understood. In this work, we systematically probed guanine positions in the human telomeric oligonucleotide sequence (hTel) by substitutions with the major product of ROS, 8-oxo-7,8-dihydroguanine (^oxoG), and evaluated the G-quadruplex forming ability of such oligonucleotides. Due to reduced hydrogen-bonding capability caused by ^oxoG, a loss of G-quadruplex structure was observed for most oligonucleotides containing oxidative lesions. However, some positions in the hTel sequence were found to tolerate substitutions with ^oxoG. Due to ^oxoG’s preference for the syn conformation, distinct responses were observed when replacing guanines with different glycosidic conformations. Accommodation of ^oxoG at sites originally in syn or anti in nonsubstituted hTel G-quadruplex requires a minor structural rearrangement or a major conformational shift, respectively. The system responds by retaining or switching to a fold where ^oxoG is in syn conformation. Most importantly, these G-quadruplex structures are still stable at physiological temperatures and should be considered detrimental in higher-order telomere structures.

Small molecule binders recognize DNA microstructural variations via an induced fit mechanism

Regulation of gene-expression by specific targeting of protein-nucleic acid interactions has been a long-standing goal in medicinal chemistry. Transcription factors are considered "undruggable" because they lack binding sites well suited for binding small-molecules. In order to overcome this obstacle, we are interested in designing small molecules that bind to the corresponding promoter sequences and either prevent or modulate transcription factor association via an allosteric mechanism. To achieve this, we must design small molecules that are both sequence-specific and able to target G/C base pair sites. A thorough understanding of the relationship between binding affinity and the structural aspects of the small molecule-DNA complex would greatly aid in rational design of such compounds. Here we present a comprehensive analysis of sequence-specific DNA association of a synthetic minor groove binder using long timescale molecular dynamics. We show how binding selectivity arises from a combination of structural factors. Our results provide a framework for the rational design and optimization of synthetic small molecules in order to improve site-specific targeting of DNA for therapeutic uses in the design of selective DNA binders targeting transcription regulation.

Dynamic Behavior of RNA Nanoparticles Analyzed by AFM on a Mica/Air Interface

RNA is an attractive biopolymer for engineering self-assembling materials suitable for biomedical applications. Previously, programmable hexameric RNA rings were developed for the controlled delivery of up to six different functionalities. To increase the potential for functionalization with little impact on nanoparticle topology, we introduce gaps into the double-stranded regions of the RNA rings. Molecular dynamic simulations are used to assess the dynamic behavior and the changes in the flexibility of novel designs. The changes suggested by simulations, however, cannot be clearly confirmed by the conventional techniques such as nondenaturing polyacrylamide gel electrophoresis (native-PAGE) and dynamic light scattering (DLS). Also, an in vitro analysis in primary cultures of human peripheral blood mononuclear cells does not reveal any discrepancy in the immunological recognition of new assemblies. To address these deficiencies, we introduce a computer-assisted quantification strategy. This strategy is based on an algorithmic atomic force microscopy (AFM)-resolved deformation analysis of the RNA nanoparticles studied on a mica/air interface. We validate this computational method by manual image analysis and fitting it to the simulation-predicted results. The presented nanoparticle modification strategy and subsequent AFM-based analysis are anticipated to provide a broad spectrum approach for the future development of nucleic acid-based nanotechnology.