Binding sites and dynamics of ammonium ions in a telomere repeat DNA quadruplex

Basic protein- and peptide-induced stabilization of long-loop DNA G-guadruplexes

The human genome contains many G-quadruplex-forming sequences, including sequences containing long single-stranded loops that are believed to be unfavorable for G-quadruplex formation. The intracellular environment of biological cells is crowded with proteins with charged surfaces. Understanding the effects of protein-rich environments is important for understanding the formation of G-quadruplexes in an intracellular environment. In this study, we investigated the structural stability of DNA G-quadruplexes in the presence of several types of globular proteins (lysozyme, cytochrome c, bovine serum albumin, myoglobin, histone proteins, and serum proteins), unstructured polypeptides (protamine and poly-l-lysine), and oligopeptides (RGG/RG-domain peptides and short repeated peptides). Thermal melting studies of G-quadruplex-forming oligonucleotides derived from the human telomeric repeat sequence revealed that environments containing high concentrations of proteins and peptides differently affected the G-quadruplex stability according to their loop lengths. We found that weak electrostatic interactions of G-quadruplex loops with basic proteins and peptides improved the stability of long-loop G-quadruplexes and the interactions were strengthened under crowded conditions simulated by dextran. The comparison of the effects of different types of proteins and peptides indicated that excluded volume interactions and structural flexibility of both DNA and polypeptide chains influenced the efficiency of their interactions. This study provides insights into long-loop G-quadruplex stability in a crowded intracellular environment and the recognition of G-quadruplexes by arginine-rich domains of G-quadruplex-binding proteins.

Formation and Nanomechanical Properties of Silver-Mediated Guanine DNA Duplexes in Aqueous Solution

Silver cations can mediate base pairing of guanine (G) DNA oligomers, yielding linear parallel G-Ag+-G duplexes with enhanced stabilities compared to those of canonical DNA duplexes. To enable their use in programmable DNA nanotechnologies, it is critical to understand solution-state formation and the nanomechanical stiffness of G-Ag+-G duplexes. Using temperature-controlled circular dichroism (CD) spectroscopy, we find that heating mixtures of G oligomers and silver salt above 50 °C fully destabilizes G-quadruplex structures and converts oligomers to G-Ag+-G duplexes. Electrospray ionization mass spectrometry supports that G-Ag+-G duplexes form at stoichiometries of 1 Ag+ per base pair, and CD spectroscopy suggests that as the Ag+/base stoichiometry increases further, G-Ag+-G duplexes undergo additional morphological changes. Using liquid-phase atomic force microscopy, we find that this excess Ag+ enables assembly of long fiberlike structures with ∼2.5 nm heights equivalent to a single DNA duplex but with lengths that far exceed a single duplex. Finally, using the conditions established to form single G-Ag+-G duplexes, we use a surface forces apparatus (SFA) to compare the solution-phase stiffness of single G-Ag+-G duplexes with dG-dC Watson-Crick-Franklin duplexes. SFA shows that G-Ag+-G duplexes are 1.3 times stiffer than dG-dC duplexes, confirming gas-phase ion mobility spectrometry measurements and computational predictions. These findings may guide the development of structural DNA nanotechnologies that rely on silver-mediated base pairing.

Insight into Tetramolecular DNA G-Quadruplexes Associated with ALS and FTLD: Cation Interactions and Formation of Higher-Ordered Structure

The G4C2 hexanucleotide repeat expansion in the c9orf72 gene is a major genetic cause of familial amyotrophic lateral sclerosis (ALS) and frontotemporal lobar degeneration (FTLD), with the formation of G-quadruplexes directly linked to the development of these diseases. Cations play a crucial role in the formation and structure of G-quadruplexes. In this study, we investigated the impact of biologically relevant potassium ions on G-quadruplex structures and utilized 15N-labeled ammonium cations as a substitute for K+ ions to gain further insights into cation binding and exchange dynamics. Through nuclear magnetic resonance spectroscopy and molecular dynamics simulations, we demonstrate that the single d(G4C2) repeat, in the presence of 15NH4+ ions, adopts a tetramolecular G-quadruplex with an all-syn quartet at the 5'-end. The movement of 15NH4+ ions through the central channel of the G-quadruplex, as well as to the bulk solution, is governed by the vacant cation binding site, in addition to the all-syn quartet at the 5'-end. Furthermore, the addition of K+ ions to G-quadruplexes folded in the presence of 15NH4+ ions induces stacking of G-quadruplexes via their 5'-end G-quartets, leading to the formation of stable higher-ordered species.

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.

Unprecedented hour-long residence time of a cation in a left-handed G-quadruplex

Cations are critical for the folding and assembly of nucleic acids. In G-quadruplex structures, cations can bind between stacked G-tetrads and coordinate with negatively charged guanine carbonyl oxygens. They usually exchange between binding sites and with the bulk in solution with time constants ranging from sub-millisecond to seconds. Here we report the first observation of extremely long-lived K+ and NH4 + ions, with an exchange time constant on the order of an hour, when coordinated at the center of a left-handed G-quadruplex DNA. A single-base mutation, that switched one half of the structure from left- to right-handed conformation resulting in a right-left hybrid G-quadruplex, was shown to remove this long-lived behaviour of the central cation.

Experimental approaches for investigating ion atmospheres around nucleic acids and proteins

Ionic interactions are crucial to biological functions of DNA, RNA, and proteins. Experimental research on how ions behave around biological macromolecules has lagged behind corresponding theoretical and computational research. In the 21st century, quantitative experimental approaches for investigating ionic interactions of biomolecules have become available and greatly facilitated examinations of theoretical electrostatic models. These approaches utilize anomalous small-angle X-ray scattering, atomic emission spectroscopy, mass spectrometry, or nuclear magnetic resonance (NMR) spectroscopy. We provide an overview on the experimental methodologies that can quantify and characterize ions within the ion atmospheres around nucleic acids, proteins, and their complexes.

A novel minimal motif for left-handed G-quadruplex formation

A recent study on the left-handed G-quadruplex (LHG4) DNA revealed a 12-nt minimal motif GTGGTGGTGGTG with the ability to independently form an LHG4 and to drive an adjacent sequence to LHG4 formation. Here we have identified a second LHG4-forming motif, GGTGGTGGTGTG, and determined the X-ray crystal structure of an LHG4 involving this motif. Our structural analysis indicated the role of split guanines and single thymine loops in promoting LHG4 formation.

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.

Dynamics of Ionic Interactions at Protein-Nucleic Acid Interfaces

Molecular association of proteins with nucleic acids is required for many biological processes essential to life. Electrostatic interactions via ion pairs (salt bridges) of nucleic acid phosphates and protein side chains are crucial for proteins to bind to DNA or RNA. Counterions around the macromolecules are also key constituents for the thermodynamics of protein–nucleic acid association. Until recently, there had been only a limited amount of experiment-based information about how ions and ionic moieties behave in biological macromolecular processes. In the past decade, there has been significant progress in quantitative experimental research on ionic interactions with nucleic acids and their complexes with proteins. The highly negatively charged surfaces of DNA and RNA electrostatically attract and condense cations, creating a zone called the ion atmosphere. Recent experimental studies were able to examine and validate theoretical models on ions and their mobility and interactions with macromolecules. The ionic interactions are highly dynamic. The counterions rapidly diffuse within the ion atmosphere. Some of the ions are released from the ion atmosphere when proteins bind to nucleic acids, balancing the charge via intermolecular ion pairs of positively charged side chains and negatively charged backbone phosphates. Previously, the release of counterions had been implicated indirectly by the salt-concentration dependence of the association constant.

Recently, direct detection of counterion release by NMR spectroscopy has become possible and enabled more accurate and quantitative analysis of the counterion release and its entropic impact on the thermodynamics of protein–nucleic acid association. Recent studies also revealed the dynamic nature of ion pairs of protein side chains and nucleic acid phosphates. These ion pairs undergo transitions between two major states. In one of the major states, the cation and the anion are in direct contact and form hydrogen bonds. In the other major state, the cation and the anion are separated by water. Transitions between these states rapidly occur on a picosecond to nanosecond time scale. When proteins interact with nucleic acids, interfacial arginine (Arg) and lysine (Lys) side chains exhibit considerably different behaviors. Arg side chains show a higher propensity to form rigid contacts with nucleotide bases, whereas Lys side chains tend to be more mobile at the molecular interfaces. The dynamic ionic interactions may facilitate adaptive molecular recognition and play both thermodynamic and kinetic roles in protein–nucleic acid interactions.

Exploration of the formation and structure characteristics of a miR-92a promoter G-quadruplex by ESI-MS and CD

Combining electrospray ionization mass spectrometry (ESI-MS) with circular dichroism (CD) spectroscopy, a G-rich sequence from miR-92a promoter region was discovered to form a parallel G-quadruplex structures in KCl or NH₄OAc solution. In case of high concentration of NH₄OAc, the ESI mass spectra showed peaks of a dimeric G-quadruplex structure with 4 ammonium ions. Meanwhile, palmatine, a natural alkaloid, was screened by ESI-MS to bind with the miR-92a G-quadruplex affinitively. The variable temperature experiment of CD also proved that high concentration of NH₄OAc or palmatine could promote the stability of the dimeric G-quadruplex structure. To get the specific characteristics of the miR-92a G-quadruplex structure, systematic mutations of guanine were tested. Based on the number of NH₄⁺ or ligands, the important guanines involved in the G-quadruplex could be determined. Considering the importance of involved guanines and the number of G-quartets, we speculated an interlocked dimeric parallel G-quadruplex as a possible conformation of the miR-92a promoter G-quadruplex. All results obtained from ESI-MS and CD illustrate structure characteristics of the miR-92a G-quadruplex, which is a promising method for preliminary structural analysis of G-quadruplexes. Besides, this study also provides a strategy for regulating the functions of microRNA by exploring and targeting higher-order structures of miRNAs.

Detecting Counterion Dynamics in DNA-Protein Association

Due to a high density of negative charges on its surface, DNA condenses cations as counterions, forming the so-called "ion atmosphere". Although the release of counterions upon DNA-protein association has been postulated to have a major contribution to the binding thermodynamics, this release remains to be confirmed through a direct observation of the ions. Herein, we report the characterization of the ion atmosphere around DNA using NMR spectroscopy and directly detect the release of counterions upon DNA-protein association. NMR-based diffusion data reveal the highly dynamic nature of counterions within the ion atmosphere around DNA. Counterion release is observed as an increase in the apparent ionic diffusion coefficient, which directly provides the number of counterions released upon DNA-protein association.

Stabilization of DNA Loop Structures by Large Cations

The DNA-binding properties of large cations differ from those of metal ions due to steric exclusion from base-paired regions. In this study, the thermal stability of DNA secondary structures, including duplexes, internal loops, bulge loops, hairpin loops, dangling ends, and G-quadruplexes, was investigated in the presence of cations of different sizes. Large cations, such as tetrabutylammonium and tetrapentylammonium ions, reduced the stability of fully matched duplexes but increased the stability of duplexes with a long loop. The cations also increased the stability of G-quadruplexes with a long loop, and the degree of stabilization was greater for low-stability G-quadruplexes. Analysis of the salt concentration dependence indicates that large cations bind to the loop nucleotides, leading to counteracting the destabilization effect on base pairing. It is likely that binding occurs when loop nucleotides are sufficiently flexible to allow for greater accessibility for large cations. These results provide insight into nucleic acid interactions with large cationic molecules and suggest a potential method for stabilizing noncanonical DNA structures under intracellular conditions.

Predicting Site-Binding Modes of Ions and Water to Nucleic Acids Using Molecular Solvation Theory

Site binding of ions and water shapes nucleic acids folding, dynamics, and biological function, complementing the more diffuse, nonspecific "territorial" ion binding. Unlike territorial binding, prediction of site-specific binding to nucleic acids remains an unsolved challenge in computational biophysics. This work presents a new toolset based on the 3D-RISM molecular solvation theory and topological analysis that predicts cation and water site binding to nucleic acids. 3D-RISM is shown to accurately capture alkali cations and water binding to the central channel, transversal loops, and grooves of the Oxytricha nova's telomeres' G-quadruplex ( Oxy-GQ), in agreement with high-resolution crystallographic data. To improve the computed cation occupancy along the Oxy-GQ central channel, it was necessary to refine and validate new cation-oxygen parameters using structural and thermodynamic data available for crown ethers and ion channels. This single set of parameters that describes both localized and delocalized binding to various biological systems is used to gain insight into cation occupancy along the Oxy-GQ channel under various salt conditions. The paper concludes with prospects for extending the method to predict divalent cation binding to nucleic acids. This work advances the forefront of theoretical methods able to provide predictive insight into ion atmosphere effects on nucleic acids function.

Conformational deviation of Thrombin binding G-quadruplex aptamer (TBA) in presence of divalent cation Sr2+: A classical molecular dynamics simulation study

Thrombin binding TBA-G-quadruplex aptamer (TBA) plays a major role in blood coagulation cascade. The 15-mer TBA sequence tends to form four-stranded TBA-G-quadruplex structure. In this research work, a series of explicit solvent classical MD simulations of the TBA is carried out using different salt (SrCl₂) concentrations (0, 50, 100 and 200 mM). Here we have also testified the effect of salt concentration of divalent cation Sr²⁺ on the conformational change of quadruplex DNA. The structural deviations, fluctuations, torsional angles and the affinity of the ion are explored at different salt concentrations. It is found that the conformation of TBA-G-quadruplex at 0 mM and 50 mM salt concentrations, is very much different than the other salt concentrations (100 mM and 200 mM). Also observed are as follows: (i) no exchange of Sr²⁺ ion between inside and outside of the channel, (ii) an enhancement in the Sr²⁺ ion density around the phosphate region of the loop residues as salt concentration increases and (iii) the stacking of T3 and T4 residues of loop-1 that appears up to 50 mM concentration, vanishes as the salt concentration is increased further.

Umbrella sampling molecular dynamics simulations reveal concerted ion movement through G-quadruplex DNA channels

We have applied the umbrella sampling (US) method in all-atom molecular dynamics (MD) simulations to obtain potential of mean force (PMF) profiles for ion transport through three representative G-quadruplex DNA channels: [d(TG₄T)]₄, [d(G₃T₄G₄)]₂, and d[G₄(T₄G₄)₃]. The US MD results are in excellent agreement with those obtained previously with the adaptive biasing force (ABF) method. We then utilized the unique features in the US MD method to investigate multi-ion effects in [d(G₃T₄G₄)]₂ and discovered that the concerted ion movement is crucial for fully explaining the unusual experimental results on ion movement in this particular G-quadruplex system. We anticipate that these modern free-energy methods will be useful tools in evaluating ion transport properties of other G-quadruplex DNA channels.

Isothermal folding of a light-up bio-orthogonal RNA origami nanoribbon

RNA presents intringuing roles in many cellular processes and its versatility underpins many different applications in synthetic biology. Nonetheless, RNA origami as a method for nanofabrication is not yet fully explored and the majority of RNA nanostructures are based on natural pre-folded RNA. Here we describe a biologically inert and uniquely addressable RNA origami scaffold that self-assembles into a nanoribbon by seven staple strands. An algorithm is applied to generate a synthetic De Bruijn scaffold sequence that is characterized by the lack of biologically active sites and repetitions larger than a predetermined design parameter. This RNA scaffold and the complementary staples fold in a physiologically compatible isothermal condition. In order to monitor the folding, we designed a new split Broccoli aptamer system. The aptamer is divided into two nonfunctional sequences each of which is integrated into the 5' or 3' end of two staple strands complementary to the RNA scaffold. Using fluorescence measurements and in-gel imaging, we demonstrate that once RNA origami assembly occurs, the split aptamer sequences are brought into close proximity forming the aptamer and turning on the fluorescence. This light-up 'bio-orthogonal' RNA origami provides a prototype that can have potential for in vivo origami applications.

Hsa-miR-1587 G-quadruplex formation and dimerization induced by NH4+, molecular crowding environment and jatrorrhizine derivatives

A guanine-rich human mature microRNA, miR-1587, was discovered to form stable intramolecular G-quadruplexes in the presence of K⁺, Na⁺ and low concentration of NH₄⁺ (25mM) by electrospray ionization mass spectrometry (ESI-MS) combined with circular dichroism (CD) spectroscopy. Furthermore, under high concentration of NH₄⁺ (100mM) or molecular crowding environments, miR-1587 formed a dimeric G-quadruplex through 3'-to-3' stacking of two monomeric G-quadruplex subunits with one ammonium ion sandwiched between the interfaces. Specifically, two synthesized jatrorrhizine derivatives with terminal amine groups could also induce the dimerization of miR-1587 G-quadruplex and formed 1:1 and 2:1 complexes with the dimeric G-quadruplex. In contrast, jatrorrhizine could bind with the dimeric miR-1587 G-quadruplex, but could not induce dimerization of miR-1587 G-quadruplex. These results provide a new strategy to regulate the functions of miR-1587 through induction of G-quadruplex formation and dimerization.

Measurement of 15N longitudinal relaxation rates in 15NH4+ spin systems to characterise rotational correlation times and chemical exchange

Many chemical and biological processes rely on the movement of monovalent cations and an understanding of such processes can therefore only be achieved by characterising the dynamics of the involved ions. It has recently been shown that 15N-ammonium can be used as a proxy for potassium to probe potassium binding in bio-molecules such as DNA quadruplexes and enzymes. Moreover, equations have been derived to describe the time-evolution of 15N-based spin density operator elements of 15NH4+ spin systems. Herein NMR pulse sequences are derived to select specific spin density matrix elements of the 15NH4+ spin system and to measure their longitudinal relaxation in order to characterise the rotational correlation time of the 15NH4+ ion as well as report on chemical exchange events of the 15NH4+ ion. Applications to 15NH4+ in acidic aqueous solutions are used to cross-validate the developed pulse sequence while measurements of spin-relaxation rates of 15NH4+ bound to a 41kDa domain of the bacterial Hsp70 homologue DnaK are presented to show the general applicability of the derived pulse sequence. The rotational correlation time obtained for 15N-ammonium bound to DnaK is similar to the correlation time that describes the rotation about the threefold axis of a methyl group. The methodology presented here provides, together with the previous theoretical framework, an important step towards characterising the motional properties of cations in macromolecular systems.

The role of loops and cation on the volume of unfolding of G-quadruplexes related to HTel

In aqueous solutions containing sodium or potassium cations, oligodeoxyribonucleotides (ODNs) rich in guanine form four-stranded DNA structures called G-quadruplexes (G4s). These structures are destabilized by elevated hydrostatic pressure. Here, we use pressure to investigate the volumetric changes arising from the formation of G4 structures. G4s display a great deal of structural heterogeneity that depends on the stabilizing cation as well as the oligonucleotide sequence. Using UV thermal unfolding at different pressures, we have investigated the volume change of the helix-coil equilibrium of a series of ODNs whose sequences are related to the G-rich ODN HTel (d[A(GGGTTA)₃GGG]), which contains four repeats of the human telomeric sequence. The experiments are conducted in aqueous buffers containing either 100mM NaCl or KCl at pH7.4. The G4s stabilized by Na⁺ are less sensitive to pressure perturbation than those stabilized by K⁺. The overall molar volume changes (ΔV_tot) of the unfolding transition for all of the G4s are large and negative. A large fraction of the measured ΔV_tot value arises from the re-hydration of the cations released from the interior of the folded structure. However, the differences in the measured ΔV_tot values demonstrate that variations in the structure of G4s formed by each ODN, arising from differences in the sequence of the loops, contribute significantly to ΔV_tot and presumably the hydration of the folded structures. Depending on the sequence of the loops, the magnitude of the measured ΔV_tot can be larger or smaller than that of HTel in solutions containing sodium. However, the magnitude of ΔV_tot is smaller than HTel for the unfolding of all G4s that are stabilized by potassium ions.

An intermolecular G-quadruplex as the basis for GTP recognition in the class V-GTP aptamer

Many naturally occurring or artificially created RNAs are capable of binding to guanine or guanine derivatives with high affinity and selectivity. They bind their ligands using very different recognition modes involving a diverse set of hydrogen bonding and stacking interactions. Apparently, the potential structural diversity for guanine, guanosine, and guanine nucleotide binding motifs is far from being fully explored. Szostak and coworkers have derived a large set of different GTP-binding aptamer families differing widely in sequence, secondary structure, and ligand specificity. The so-called class V-GTP aptamer from this set binds GTP with very high affinity and has a complex secondary structure. Here we use solution NMR spectroscopy to demonstrate that the class V aptamer binds GTP through the formation of an intermolecular two-layered G-quadruplex structure that directly incorporates the ligand and folds only upon ligand addition. Ligand binding and G-quadruplex formation depend strongly on the identity of monovalent cations present with a clear preference for potassium ions. GTP binding through direct insertion into an intermolecular G-quadruplex is a previously unobserved structural variation for ligand-binding RNA motifs and rationalizes the previously observed specificity pattern of the class V aptamer for GTP analogs.

Metal Cations in G-Quadruplex Folding and Stability

This review is focused on the structural and physicochemical aspects of metal cation coordination to G-Quadruplexes (GQ) and their effects on GQ stability and conformation. G-quadruplex structures are non-canonical secondary structures formed by both DNA and RNA. G-quadruplexes regulate a wide range of important biochemical processes. Besides the sequence requirements, the coordination of monovalent cations in the GQ is essential for its formation and determines the stability and polymorphism of GQ structures. The nature, location, and dynamics of the cation coordination and their impact on the overall GQ stability are dependent on several factors such as the ionic radii, hydration energy, and the bonding strength to the O6 of guanines. The intracellular monovalent cation concentration and the localized ion concentrations determine the formation of GQs and can potentially dictate their regulatory roles. A wide range of biochemical and biophysical studies on an array of GQ enabling sequences have generated at a minimum the knowledge base that allows us to often predict the stability of GQs in the presence of the physiologically relevant metal ions, however, prediction of conformation of such GQs is still out of the realm.

Divalent cation-aided identification of physico-chemical properties of metal ions that stabilize RNA G-quadruplexes

DNA and RNA sequences rich in guanosines (G) can form a four-stranded secondary structure known as a G-quadruplex (GQ), which plays a role in regulation of gene expression at the transcription and translation level. Both DNA and RNA GQs typically use the monovalent K(+) ion for stabilization of the structures. However, the fundamental reasons for K(+) acting as the most stabilizing metal ion for RNA GQs are not known. To identify the properties of a metal ion that stabilizes an RNA GQ we investigated the effect of alkaline earth metal cations and a set of divalent transition metal ions on two previously identified highly stable RNA GQs. Our results based upon circular dichroism and RNase T1 structure mapping data reveal that the RNA GQs are destabilized in the presence of the tested divalent metal cations. The destabilizing effect of a divalent metal cation is reversible upon increasing K(+) concentration. Results show that ionic radius, hydration energy, and binding strength towards the hard ligand (guanine O(6)) are important factors that determine a metal ion's ability to stabilize an RNA GQ. Additionally, the tested set of divalent metal cations incongruously affects RNA and DNA GQs.

Aromatic Base Stacking in DNA: From ab initio Calculations to Molecular Dynamics Simulations

Abstract Aromatic stacking of nucleic acid bases is one of the key players in determining the structure and dynamics of nucleic acids. The arrangement of nucleic acid bases with extensive overlap of their aromatic rings gave rise to numerous often contradictory suggestions about the physical origins of stacking and the possible role of delocalized electrons in stacked aromatic π systems, leading to some confusion about the issue. The recent advance of computer hardware and software finally allowed the application of state of the art quantum-mechanical approaches with inclusion of electron correlation effects to study aromatic base stacking, now providing an ultimitate qualitative description of the phenomenon. Base stacking is determined by an interplay of the three most commonly encountered molecular interactions: dispersion attraction, electrostatic interaction, and short-range repulsion. Unusual (aromatic- stacking specific) energy contributions were in fact not evidenced and are not necessary to describe stacking. The currently used simple empirical potential form, relying on atom-centered constant point charges and Lennard-Jones van der Waals terms, is entirely able to reproduce the essential features of base stacking. Thus, we can conclude that base stacking is in principle one of the best described interactions in current molecular modeling and it allows to study base stacking in DNA using large-scale classical molecular dynamics simulations. Neglect of cooperativity of stacking appears to be the most serious approximation of the currently used force field form. This review summarizes recent developments in the field. It is written for an audience that is not necessarily expert in computational quantum chemistry and follows up on our previous contribution (Sponer et. al., J. Biomol. Struct. Dyn. 14, 117, (1997)). First, the applied methodology, its accuracy, and the physical nature of base stacking is briefly overviewed, including a comment on the accuracy of other molecular orbital methods and force fields. Then, base stacking is contrasted with hydrogen bonding, the other dominant force in nucleic acid structure. The sequence dependence and cooperativity of base stacking is commented on, and finally a brief introduction into recent progress in large-scale molecular dynamics simulations of nucleic acids is provided. Using four stranded DNA assemblies as an example, we demonstrate the efficacy of current molecular dynamics techniques that utilize refined and verified force fields in the study of stacking in nucleic acid molecules.

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.

Higher-order human telomeric G-quadruplex DNA metalloenzyme catalyzed Diels–Alder reaction: an unexpected inversion of enantioselectivity modulated by K+ and NH4+ ions

K⁺ and NH₄⁺, bearing approximately equal ionic radius, present different allosteric activation for higher-order human telomeric G-quadruplex DNA metalloenzyme.

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.

DNA switches: from principles to applications

The base sequence of nucleic acid encodes structural and functional properties into the biopolymer. Structural information includes the formation of duplexes, G-quadruplexes, i-motif, and cooperatively stabilized assemblies. Functional information encoded in the base sequence involves the strand-displacement process, the recognition properties by aptamers, and the catalytic functions of DNAzymes. This Review addresses the implementation of the information encoded in nucleic acids to develop DNA switches. A DNA switch is a supramolecular nucleic acid assembly that undergoes cyclic, switchable, transitions between two distinct states in the presence of appropriate triggers and counter triggers, such as pH value, metal ions/ligands, photonic and electrical stimuli. Applications of switchable DNA systems to tailor switchable DNA hydrogels, for the controlled drug-release and for the activation of switchable enzyme cascades, are described, and future perspectives of the systems are addressed.

G-triplex structure and formation propensity

The occurrence of a G-triplex folding intermediate of thrombin binding aptamer (TBA) has been recently predicted by metadynamics calculations, and experimentally supported by Nuclear Magnetic Resonance (NMR), Circular Dichroism (CD) and Differential Scanning Calorimetry (DSC) data collected on a 3' end TBA-truncated 11-mer oligonucleotide (11-mer-3'-t-TBA). Here we present the solution structure of 11-mer-3'-t-TBA in the presence of potassium ions. This structure is the first experimental example of a G-triplex folding, where a network of Hoogsteen-like hydrogen bonds stabilizes six guanines to form two G:G:G triad planes. The G-triplex folding of 11-mer-3'-t-TBA is stabilized by the potassium ion and destabilized by increasing the temperature. The superimposition of the experimental structure with that predicted by metadynamics shows a great similarity, with only significant differences involving two loops. These new structural data show that 11-mer-3'-t-TBA assumes a G-triplex DNA conformation as its stable form, reinforcing the idea that G-triplex folding intermediates may occur in vivo in human guanine-rich sequences. NMR and CD screening of eight different constructs obtained by removing from one to four bases at either the 3' and the 5' ends show that only the 11-mer-3'-t-TBA yields a relatively stable G-triplex.

Heteronuclear transverse and longitudinal relaxation in AX4 spin systems: application to (15)N relaxations in (15)NH4(+)

The equations that describe the time-evolution of transverse and longitudinal (15)N magnetisations in tetrahedral ammonium ions, (15)NH4(+), are derived from the Bloch-Wangsness-Redfield density operator relaxation theory. It is assumed that the relaxation of the spin-states is dominated by (1) the intra-molecular (15)N-(1)H and (1)H-(1)H dipole-dipole interactions and (2) interactions of the ammonium protons with remote spins, which also include the contribution to the relaxations that arise from the exchange of the ammonium protons with the bulk solvent. The dipole-dipole cross-correlated relaxation mechanisms between each of the (15)N-(1)H and (1)H-(1)H interactions are explicitly taken into account in the derivations. An application to (15)N-ammonium bound to a 41kDa domain of the protein DnaK is presented, where a comparison between experiments and simulations show that the ammonium ion rotates rapidly within its binding site with a local correlation time shorter than approximately 1ns. The theoretical framework provided here forms the basis for further investigations of dynamics of AX4 spin systems, with ammonium ions in solution and bound to proteins of particular interest.

A G-quadruplex-containing RNA activates fluorescence in a GFP-like fluorophore

Spinach is an in vitro-selected RNA aptamer that binds a GFP-like ligand and activates its green fluorescence. Spinach is thus an RNA analog of GFP and has potentially widespread applications for in vivo labeling and imaging. We used antibody-assisted crystallography to determine the structures of Spinach both with and without bound fluorophore at 2.2-Å and 2.4-Å resolution, respectively. Spinach RNA has an elongated structure containing two helical domains separated by an internal bulge that folds into a G-quadruplex motif of unusual topology. The G-quadruplex motif and adjacent nucleotides comprise a partially preformed binding site for the fluorophore. The fluorophore binds in a planar conformation and makes extensive aromatic stacking and hydrogen bond interactions with the RNA. Our findings provide a foundation for structure-based engineering of new fluorophore-binding RNA aptamers.

Using ¹⁵N-ammonium to characterise and map potassium binding sites in proteins by NMR spectroscopy

A variety of enzymes are activated by the binding of potassium ions. The potassium binding sites of these enzymes are very specific, but ammonium ions can often replace potassium ions in vitro because of their similar ionic radii. In these cases, ammonium can be used as a proxy for potassium to characterise potassium binding sites in enzymes: the (1) H,(15) N spin-pair of enzyme-bound (15) NH4 (+) can be probed by (15) N-edited heteronuclear NMR experiments. Here, we demonstrate the use of NMR spectroscopy to characterise binding of ammonium ions to two different enzymes: human histone deacetylase 8 (HDAC8), which is activated allosterically by potassium, and the bacterial Hsp70 homologue DnaK, for which potassium is an integral part of the active site. Ammonium activates both enzymes in a similar way to potassium, thus supporting this non-invasive approach. Furthermore, we present an approach to map the observed binding site onto the structure of HDAC8. Our method for mapping the binding site is general and does not require chemical shift assignment of the enzyme resonances.

Cations form sequence selective motifs within DNA grooves via a combination of cation-pi and ion-dipole/hydrogen bond interactions

The fine conformational subtleties of DNA structure modulate many fundamental cellular processes including gene activation/repression, cellular division, and DNA repair. Most of these cellular processes rely on the conformational heterogeneity of specific DNA sequences. Factors including those structural characteristics inherent in the particular base sequence as well as those induced through interaction with solvent components combine to produce fine DNA structural variation including helical flexibility and conformation. Cation-pi interactions between solvent cations or their first hydration shell waters and the faces of DNA bases form sequence selectively and contribute to DNA structural heterogeneity. In this paper, we detect and characterize the binding patterns found in cation-pi interactions between solvent cations and DNA bases in a set of high resolution x-ray crystal structures. Specifically, we found that monovalent cations (Tl⁺) and the polarized first hydration shell waters of divalent cations (Mg²⁺, Ca²⁺) form cation-pi interactions with DNA bases stabilizing unstacked conformations. When these cation-pi interactions are combined with electrostatic interactions a pattern of specific binding motifs is formed within the grooves.

G-quadruplex formation of entirely hydrophobic DNA in organic solvents

We report herein a novel hydrophobic G-quadruplex DNA consisting of dodecyl phosphotriester linkages. The 6-mer hydrophobic DNA having a TG4T sequence binds to monovalent cations to form a tetramolecular G-quadruplex in low polarity organic solvents.

Ammonium ion binding to DNA G-quadruplexes: do electrospray mass spectra faithfully reflect the solution-phase species?

G-quadruplex nucleic acids can bind ammonium ions in solution, and these complexes can be detected by electrospray mass spectrometry (ESI-MS). However, because ammonium ions are volatile, the extent to which ESI-MS quantitatively could provide an accurate reflection of such solution-phase equilibria is unclear. Here we studied five G-quadruplexes having known solution-phase structure and ammonium ion binding constants: the bimolecular G-quadruplexes (dG(4)T(4)G(4))(2), (dG(4)T(3)G(4))(2), and (dG(3)T(4)G(4))(2), and the intramolecular G-quadruplexes dG(4)(T(4)G(4))(3) and dG(2)T(2)G(2)TGTG(2)T(2)G(2) (thrombin binding aptamer). We found that not all mass spectrometers are equally suited to reflect the solution phase species. Ion activation can occur in the electrospray source, or in a high-pressure traveling wave ion mobility cell. When the softest instrumental conditions are used, ammonium ions bound between G-quartets, but also additional ammonium ions bound at specific sites outside the external G-quartets, can be observed. However, even specifically bound ammonium ions are in some instances too labile to be fully retained in the gas phase structures, and although the ammonium ion distribution observed by ESI-MS shows biases at specific stoichiometries, the relative abundances in solution are not always faithfully reflected. Ion mobility spectrometry results show that all inter-quartet ammonium ions are necessary to preserve the G-quadruplex fold in the gas phase. Ion mobility experiments, therefore, help assign the number of inner ammonium ions in the solution phase structure.

Strand directionality affects cation binding and movement within tetramolecular G-quadruplexes

Nuclear magnetic resonance study of G-quadruplex structures formed by d(TG₃T) and its modified analogs containing a 5′-5′ or 3′-3′ inversion of polarity sites, namely d(3′TG5′-5′G₂T3′), d(3′T5′-5′G₃T3′) and d(5′TG3′-3′G₂T5’) demonstrates formation of G-quadruplex structures with tetrameric topology and distinct cation-binding preferences. All oligonucleotides are able to form quadruplex structures with two binding sites, although the modified oligonucleotides also form, in variable amounts, quadruplex structures with only one bound cation. The inter-quartet cavities at the inversion of polarity sites bind ammonium ions less tightly than a naturally occurring 5′-3′ backbone. Exchange of ¹⁵ ions between G-quadruplex and bulk solution is faster at the 3′-end in comparison to the 5′-end. In addition to strand directionality, cation movement is influenced by formation of an all-syn G-quartet. Formation of such quartet has been observed also for the parent d(TG₃T) that besides the canonical quadruplex with only all-anti G-quartets, forms a tetramolecular parallel quadruplex containing one all-syn G-quartet, never observed before in unmodified quadruplex structures.

Investigating the interactions between cations, peroxidation substrates and G-quadruplex topology in DNAzyme peroxidation reactions using statistical testing

Cations affect the topology and enzymatic proficiency of most macromolecular catalysts but the role of cations in DNAzyme peroxidation reactions remains unresolved. Herein, we use statistical methods (ANOVA, t-test and Wilcoxon Mann-Whitney non-parametric test) to demonstrate that there are strong associations between cations, DNAzyme topology, peroxidation substrate and peroxidation rates of G-quadruplex peroxidises. Ammonium cation was found to be superior to all tested cations, including potassium. A t-test indicated that NH(4)(+) was better than K(+) with a p-value=0.05. Interestingly, the nature of the peroxidation substrate employed affected the dependence of peroxidation rate on the cation present and of the three substrates tested, 2,2'-azino-bis (3-ethylbenzothiazoline-6-sulfonic acid (ABTS), tyramine and 3,3',5,5'-tetramethylbenzidine (TMB), ABTS was the most sensitive to the nature of cation present.

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.

Mg2+ in the major groove modulates B-DNA structure and dynamics

This study investigates the effect of Mg(2+) bound to the DNA major groove on DNA structure and dynamics. The analysis of a comprehensive dataset of B-DNA crystallographic structures shows that divalent cations are preferentially located in the DNA major groove where they interact with successive bases of (A/G)pG and the phosphate group of 5'-CpA or TpG. Based on this knowledge, molecular dynamics simulations were carried out on a DNA oligomer without or with Mg(2+) close to an ApG step. These simulations showed that the hydrated Mg(2+) forms a stable intra-strand cross-link between the two purines in solution. ApG generates an electrostatic potential in the major groove that is particularly attractive for cations; its intrinsic conformation is well-adapted to the formation of water-mediated hydrogen bonds with Mg(2+). The binding of Mg(2+) modulates the behavior of the 5'-neighboring step by increasing the BII (ε-ζ>0°) population of its phosphate group. Additional electrostatic interactions between the 5'-phosphate group and Mg(2+) strengthen both the DNA-cation binding and the BII character of the 5'-step. Cation binding in the major groove may therefore locally influence the DNA conformational landscape, suggesting a possible avenue for better understanding how strong DNA distortions can be stabilized in protein-DNA complexes.