Molecular dynamics simulations reveal the balance of forces governing the formation of a guanine tetrad—a common structural unit of G-quadruplex DNA

G-quadruplexes (G4) are nucleic acid conformations of guanine-rich sequences, in which guanines are arranged in the square-planar G-tetrads, stacked on one another. G4 motifs form in vivo and are implicated in regulation of such processes as gene expression and chromosome maintenance. The structure and stability of various G4 topologies were determined experimentally; however, the driving forces for their formation are not fully understood at the molecular level. Here, we used all-atom molecular dynamics to probe the microscopic origin of the G4 motif stability. By computing the free energy profiles governing the dissociation of the 3′-terminal G-tetrad in the telomeric parallel-stranded G4, we examined the thermodynamic and kinetic stability of a single G-tetrad, as a common structural unit of G4 DNA. Our results indicate that the energetics of guanine association alone does not explain the overall stability of the G-tetrad and that interactions involving sugar–phosphate backbone, in particular, the constrained minimization of the phosphate–phosphate repulsion energy, are crucial in providing the observed enthalpic stabilization. This enthalpic gain is largely compensated by the unfavorable entropy change due to guanine association and optimization of the backbone topology.     

Human telomere, oncogenic promoter and 5'-UTR G-quadruplexes: diverse higher order DNA and RNA targets for cancer therapeutics

Guanine-rich DNA sequences can form G-quadruplexes stabilized by stacked G–G–G–G tetrads in monovalent cation-containing solution. The length and number of individual G-tracts and the length and sequence context of linker residues define the diverse topologies adopted by G-quadruplexes. The review highlights recent solution NMR-based G-quadruplex structures formed by the four-repeat human telomere in K⁺ solution and the guanine-rich strands of c-myc, c-kit and variant bcl-2 oncogenic promoters, as well as a bimolecular G-quadruplex that targets HIV-1 integrase. Such structure determinations have helped to identify unanticipated scaffolds such as interlocked G-quadruplexes, as well as novel topologies represented by double-chain-reversal and V-shaped loops, triads, mixed tetrads, adenine-mediated pentads and hexads and snap-back G-tetrad alignments. The review also highlights the recent identification of guanine-rich sequences positioned adjacent to translation start sites in 5′-untranslated regions (5′-UTRs) of RNA oncogenic sequences. The activity of the enzyme telomerase, which maintains telomere length, can be negatively regulated through G-quadruplex formation at telomeric ends. The review evaluates progress related to ongoing efforts to identify small molecule drugs that bind and stabilize distinct G-quadruplex scaffolds associated with telomeric and oncogenic sequences, and outlines progress towards identifying recognition principles based on several X-ray-based structures of ligand–G-quadruplex complexes.     

Effect of Coordinated Ions on Structure and Flexibility of Parallel G-quadruplexes: A Molecular Dynamics Study

Abstract Single tract guanine residues can associate to form stable parallel quadruplex structures in the presence of certain cations. Nanosecond scale molecular dynamics simulations have been performed on fully solvated fibre model of parallel d(G(7)) quadruplex structures with Na(+) or K(+) ions coordinated in the cavity formed by the O6 atoms of the guanine bases. The AMBER 4.1 force field and Particle Mesh Ewald technique for electrostatic interactions have been used in all simulations. These quadruplex structures are stable during the simulation, with the middle four base tetrads showing root mean square deviation values between 0.5 to 0.8 ? from the initial structure as well the high resolution crystal structure. Even in the absence of any coordinated ion in the initial structure, the G-quadruplex structure remains intact throughout the simulation. During the 1.1 ns MD simulation, one Na(+) counter ion from the solvent as well as several water molecules enter the central cavity to occupy the empty coordination sites within the parallel quadruplex and help stabilize the structure. Hydrogen bonding pattern depends on the nature of the coordinated ion, with the G-tetrad undergoing local structural variation to accommodate cations of different sizes. In the absence of any coordinated ion, due to strong mutual repulsion, O6 atoms within G-tetrad are forced farther apart from each other, which leads to a considerably different hydrogen bonding scheme within the G-tetrads and very favourable interaction energy between the guanine bases constituting a G-tetrad. However, a coordinated ion between G-tetrads provides extra stacking energy for the G-tetrads and makes the quadruplex structure more rigid. Na(+) ions, within the quadruplex cavity, are more mobile than coordinated K(+) ions. A number of hydrogen bonded water molecules are observed within the grooves of all quadruplex structures.     

Effect of Coordinated Ions on Structure and Flexibility of Parallel G-quadruplexes: A Molecular Dynamics Study

Abstract

Single tract guanine residues can associate to form stable parallel quadruplex structures in the presence of certain cations. Nanosecond scale molecular dynamics simulations have been performed on fully solvated fibre model of parallel d(G₇) quadruplex structures with Na⁺ or K⁺ ions coordinated in the cavity formed by the O6 atoms of the guanine bases. The AMBER 4.1 force field and Particle Mesh Ewald technique for electrostatic interactions have been used in all simulations. These quadruplex structures are stable during the simulation, with the middle four base tetrads showing root mean square deviation values between 0.5 to 0.8 Å from the initial structure as well the high resolution crystal structure. Even in the absence of any coordinated ion in the initial structure, the G-quadruplex structure remains intact throughout the simulation. During the 1.1 ns MD simulation, one Na⁺ counter ion from the solvent as well as several water molecules enter the central cavity to occupy the empty coordination sites within the parallel quadruplex and help stabilize the structure. Hydrogen bonding pattern depends on the nature of the coordinated ion, with the G-tetrad undergoing local structural variation to accommodate cations of different sizes. In the absence of any coordinated ion, due to strong mutual repulsion, O6 atoms within G-tetrad are forced farther apart from each other, which leads to a considerably different hydrogen bonding scheme within the G-tetrads and very favourable interaction energy between the guanine bases constituting a G-tetrad. However, a coordinated ion between G-tetrads provides extra stacking energy for the G-tetrads and makes the quadruplex structure more rigid. Na⁺ ions, within the quadruplex cavity, are more mobile than coordinated K⁺ ions. A number of hydrogen bonded water molecules are observed within the grooves of all quadruplex structures.     

Effect of coordinated ions on structure and flexiblity of parallel G-quandruplexes: a molecular dynamics study

Single tract guanine residues can associate to form stable parallel quadruplex structures in the presence of certain cations. Nanosecond scale molecular dynamics simulations have been performed on fully solvated fibre model of parallel d(G7) quadruplex structures with Na+ or K+ ions coordinated in the cavity formed by the 06 atoms of the guanine bases. The AMBER 4.1 force field and Particle Mesh Ewald technique for electrostatic interactions have been used in all simulations. These quadruplex structures are stable during the simulation, with the middle four base tetrads showing root mean square deviation values between 0.5 to 0.8 A from the initial structure as well the high resolution crystal structure. Even in the absence of any coordinated ion in the initial structure, the G-quadruplex structure remains intact throughout the simulation. During the 1.1 ns MD simulation, one Na+ counter ion from the solvent as well as several water molecules enter the central cavity to occupy the empty coordination sites within the parallel quadruplex and help stabilize the structure. Hydrogen bonding pattern depends on the nature of the coordinated ion, with the G-tetrad undergoing local structural variation to accommodate cations of different sizes. In the absence of any coordinated ion, due to strong mutual repulsion, 06 atoms within G-tetrad are forced farther apart from each other, which leads to a considerably different hydrogen bonding scheme within the G-tetrads and very favourable interaction energy between the guanine bases constituting a G-tetrad. However, a coordinated ion between G-tetrads provides extra stacking energy for the G-tetrads and makes the quadruplex structure more rigid. Na+ ions, within the quadruplex cavity, are more mobile than coordinated K+ ions. A number of hydrogen bonded water molecules are observed within the grooves of all quadruplex structures.     

Xanthine and 8-oxoguanine in G-quadruplexes: formation of a G·G·X·O tetrad

G-quadruplexes are four-stranded structures built from stacked G-tetrads (G·G·G·G), which are planar cyclical assemblies of four guanine bases interacting through Hoogsteen hydrogen bonds. A G-quadruplex containing a single guanine analog substitution, such as 8-oxoguanine (O) or xanthine (X), would suffer from a loss of a Hoogsteen hydrogen bond within a G-tetrad and/or potential steric hindrance. We show that a proper arrangement of O and X bases can reestablish the hydrogen-bond pattern within a G·G·X·O tetrad. Rational incorporation of G·G·X·O tetrads in a (3+1) G-quadruplex demonstrated a similar folding topology and thermal stability to that of the unmodified G-quadruplex. pH titration conducted on X·O-modified G-quadruplexes indicated a protonation-deprotonation equilibrium of X with a pKa ∼6.7. The solution structure of a G-quadruplex containing a G·G·X·O tetrad was determined, displaying the same folding topology in both the protonated and deprotonated states. A G-quadruplex containing a deprotonated X·O pair was shown to exhibit a more electronegative groove compared to that of the unmodified one. These differences are likely to manifest in the electronic properties of G-quadruplexes and may have important implications for drug targeting and DNA-protein interactions.     

Stacking of G-quadruplexes: NMR structure of a G-rich oligonucleotide with potential anti-HIV and anticancer activity

G-rich oligonucleotides T30695 (or T30923), with the sequence of (GGGT)(4), and T40214, with the sequence of (GGGC)(4), have been reported to exhibit anti-HIV and anticancer activity. Here we report on the structure of a dimeric G-quadruplex adopted by a derivative of these sequences in K(+) solution. It comprises two identical propeller-type parallel-stranded G-quadruplex subunits each containing three G-tetrad layers that are stacked via the 5'-5' interface. We demonstrated control over the stacking of the two monomeric subunits by sequence modifications. Our analysis of possible structures at the stacking interface provides a general principle for stacking of G-quadruplexes, which could have implications for the assembly and recognition of higher-order G-quadruplex structures.     

The crystal structure of human telomeric DNA complexed with berberine: an interesting case of stacked ligand to G-tetrad ratio higher than 1:1

The first crystal structure of human telomeric DNA in complex with the natural alkaloid berberine, produced by different plant families and used in folk medicine for millennia, was solved by X-ray diffraction method. The G-quadruplex unit features all-parallel strands. The overall folding assumed by DNA is the same found in previously reported crystal structures. Similarly to previously reported structures the ligand molecules were found to be stacked onto the external 5′ and 3′-end G-tetrads. However, the present crystal structure highlighted for the first time, the presence of two berberine molecules in the two binding sites, directly interacting with each tetrad. As a consequence, our structural data point out a 2:1 ligand to G-tetrad molar ratio, which has never been reported before in a telomeric intramolecular quadruplex structure.     

Highly fluorescent guanosine mimics for folding and energy transfer studies

Guanosines with substituents at the 8-position can provide useful fluorescent probes that effectively mimic guanine residues even in highly demanding model systems such as polymorphic G-quadruplexes and duplex DNA. Here, we report the synthesis and photophysical properties of a small family of 8-substituted-2′-deoxyguanosines that have been incorporated into the human telomeric repeat sequence using phosphoramidite chemistry. These include 8-(2-pyridyl)-2′-deoxyguanosine (2PyG), 8-(2-phenylethenyl)-2′-deoxyguanosine (StG) and 8-[2-(pyrid-4-yl)-ethenyl]-2′-deoxyguanosine (4PVG). On DNA folding and stability, 8-substituted guanosines can exhibit context-dependent effects but were better tolerated by G-quadruplex and duplex structures than pyrimidine mismatches. In contrast to previously reported fluorescent guanine analogs, 8-substituted guanosines exhibit similar or even higher quantum yields upon their incorporation into nucleic acids (Φ = 0.02–0.45). We have used these highly emissive probes to quantify energy transfer efficiencies from unmodified DNA nucleobases to 8-substituted guanosines. The resulting DNA-to-probe energy transfer efficiencies (η_t) are highly structure selective, with η_t(duplex) < η_t(single-strand) < η_t(G-quadruplex). These trends were independent of the exact structural features and thermal stabilities of the G-quadruplexes or duplexes containing them. The combination of efficient energy transfer, high probe quantum yield, and high molar extinction coefficient of the DNA provides a highly sensitive and reliable readout of G-quadruplex formation even in highly diluted sample solutions of 0.25 nM.     

The high-resolution crystal structure of a parallel intermolecular DNA G-4 quadruplex/drug complex employing syn glycosyl linkages

We have determined the X-ray structure of the complex between the DNA quadruplex d(5′-GGGG-3′)₄ and daunomycin, as a potential model for studying drug–telomere interactions. The structure was solved at 1.08 Å by direct methods in space group I4. The asymmetric unit comprises a linear arrangement of one d(GGGG) strand, four daunomycin molecules, a second d(GGGG) strand facing in the opposite direction to the first, and Na and Mg cations. The crystallographic 4-fold axis generates the biological unit, which is a 12-layered structure comprising two sets of four guanine layers, with four layers each of four daunomycins stacked between the 5′ faces of the two quadruplexes. The daunomycin layers fall into two groups which are novel in their mode of self assembly. The only contacts between daunomycin molecules within any one of these layers are van der Waals interactions, however there is substantial π–π stacking between successive daunomycin layers and also with adjacent guanine layers. The structure differs significantly from all other parallel d(TGGGGT)₄ quadruplexes in that the 5′ guanine adopts the unusual syn glycosyl linkage, refuting the widespread belief that such conformations should all be anti. In contrast to the related d(TGGGGT)/daunomycin complex, there are no ligand–quadruplex groove insertion interactions.     

G-quadruplexes with (4n?- 1) guanines in the G-tetrad core: formation of a G-triad·water complex and implication for small-molecule binding

G-quadruplexes are non-canonical structures of nucleic acids, in which guanine bases form planar G-tetrads (G·G·G·G) that stack on each other in the core of the structure. G-quadruplexes generally contain multiple times of four (4n) guanines in the core. Here, we study the structure of G-quadruplexes with only (4n - 1) guanines in the core. The solution structure of a DNA sequence containing 11 guanines showed the formation of a parallel G-quadruplex involving two G-tetrads and one G-triad with a vacant site. Molecular dynamics simulation established the formation of a stable G-triad·water complex, where water molecules mimic the position of the missing guanine in the vacant site. The concept of forming G-quadruplexes with missing guanines in the core broadens the current definition of G-quadruplex-forming sequences. The potential ability of such structures to bind different metabolites, including guanine, guanosine and GTP, in the vacant site, could have biological implications in regulatory functions. Formation of this unique binding pocket in the G-triad could be used as a specific target in drug design.     

Two-repeat Tetrahymena telomeric d(TGGGGTTGGGGT) Sequence interconverts between asymmetric dimeric G-quadruplexes in solution

Recently, the two-repeat human telomeric d(TAGGGTTAGGGT) sequence has been shown to form interconverting parallel and antiparallel G-quadruplex structures in solution. Here, we examine the structures formed by the two-repeat Tetrahymena telomeric d(TGGGGTTGGGGT) sequence, which differs from the human sequence only by one G-for-A replacement in each repeat. We show by NMR that this sequence forms two novel G-quadruplex structures in Na+-containing solution. Both structures are asymmetric, dimeric G-quadruplexes involving a core of four stacked G-tetrads and two edgewise loops. The adjacent strands of the G-tetrad core are alternately parallel and antiparallel. All G-tetrads adopt syn.syn.anti.anti alignments, which occur with 5'-syn-anti-syn-anti-3' alternations along G-tracks. In the first structure (head-to-head), two loops are at one end of the G-tetrad core; in the second structure (head-to-tail), two loops are located on opposite ends of the G-tetrad core. In contrast to the human telomere counterpart, the proportions of the two forms here are similar for a wide range of temperatures; their unfolding rates are also similar, with an activation enthalpy of 153 kJ/mol.     

Interaction of cyclic cytosine-, guanine-, thymine-, uracil- and mixed guanine-cytosine base tetrads with K+, Na+ and Li+ ions -- a density functional study

We have carried out B3LYP hybrid density functional studies of complexes formed by cyclic cytosine-, guanine-, thymine-, uracil- and mixed guanine cytosine-tetrads with Li+, Na+ and K+ ions to determine their structures and interaction energies. The conformations studied have been restricted to a hydrogen bond pattern closely related to the tetrads observed in experimental nucleic acid structures. A comparison of the alkali metal ion/tetrad complexes with the tetrads without cations indicates that alkali metal ions modulate the tetrad structures significantly and that even the hydrogen bond pattern may change. Guanine-tetrad cation complexes show the strongest interaction energy compared to other tetrads that occur less frequently in experimental structures. The most stable G-tetrad/metal ion structure adopts a nearly planar geometry that is especially suitable for tetraplex formation, which requires approximately parallel tetrad planes. In the cytosine-tetrad there is a very large central cavity suitable for cation recognition, but the complexes adopt a non-planar structure unsuitable for stacking, except possibly for ions with very large radii. Uracil and thymine tetrads show a significant different characteristics which may contribute to the differences between DNA and RNA     

Overlapping but distinct: a new model for G-quadruplex biochemical specificity

G-quadruplexes are noncanonical nucleic acid structures formed by stacked guanine tetrads. They are capable of a range of functions and thought to play widespread biological roles. This diversity raises an important question: what determines the biochemical specificity of G-quadruplex structures? The answer is particularly important from the perspective of biological regulation because genomes can contain hundreds of thousands of G-quadruplexes with a range of functions. Here we analyze the specificity of each sequence in a 496-member library of variants of a reference G-quadruplex with respect to five functions. Our analysis shows that the sequence requirements of G-quadruplexes with these functions are different from one another, with some mutations altering biochemical specificity by orders of magnitude. Mutations in tetrads have larger effects than mutations in loops, and changes in specificity are correlated with changes in multimeric state. To complement our biochemical data we determined the solution structure of a monomeric G-quadruplex from the library. The stacked and accessible tetrads rationalize why monomers tend to promote a model peroxidase reaction and generate fluorescence. Our experiments support a model in which the sequence requirements of G-quadruplexes with different functions are overlapping but distinct. This has implications for biological regulation, bioinformatics, and drug design.     

Interaction of Cyclic Cytosine-, Guanine-, Thymine-, Uracil- and Mixed Guanine-Cytosine Base Tetrads with K+, Na+ and Li+ Ions—A Density Functional Study

Abstract

We have carried out B3LYP hybrid density functional studies of complexes formed by cyclic cytosine-, guanine-, thymine-, uracil- and mixed guanine cytosine-tetrads with Li⁺, Na⁺ and K⁺ ions to determine their structures and interaction energies. The conformations studied have been restricted to a hydrogen bond pattern closely related to the tetrads observed in experimental nucleic acid structures. A comparison of the alkali metal ion/tetrad complexes with the tetrads without cations indicates that alkali metal ions modulate the tetrad structures significantly and that even the hydrogen bond pattern may change. Guanine-tetrad cation complexes show the strongest interaction energy compared to other tetrads that occur less frequently in experimental structures. The most stable G-tetrad/metal ion structure adopts a nearly planar geometry that is especially suitable for tetraplex formation, which requires approximately parallel tetrad planes. In the cytosine-tetrad there is a very large central cavity suitable for cation recognition, but the complexes adopt a non-planar structure unsuitable for stacking, except possibly for ions with very large radii. Uracil and thymine tetrads show a significant different characteristics which may contribute to the differences between DNA and RNA.     

A comparison of the properties and the solution structure for RNA and DNA quadruplexes which contain two GGAGG sequences joined with a tetranucleotide linker

We have determined solution structure of r(GGAGGUUUUGGAGG) (R14) by NMR; the RNA 14-mer forms an intra-strand parallel quadruplex with a G-tetrad and a hexad, in which a G-tetrad core is augmented by association of two A residues. The quadruplex further forms a dimer through stacking interaction between the hexads. In order to obtain insight into the difference between RNA and DNA quadruplexes, we synthesized the corresponding DNA 14-mer, d(GGAGGTTTTGGAGG) (D14), and examined its properties and structure by CD, gel electrophoresis, and NMR. K+ ions increased the thermal stability of both R14 and D14 structures. The binding affinity of K+ ions to R14 was much higher than that to D14. The CD and gel electrophoretic studies suggest that D14 forms a quadruplex entirely different from that of R14 in the presence of K+ ions; two molecules of D14 form a quadruplex with both antiparallel and parallel strand alignments and with diagonal loops at both ends of the stacked G-tetrads. The NMR study also gave results that are consistent with such structure: alternate glycosidic conformation, 5'G(syn)-G(anti)3', and characteristic chemical shift data observed for many quadruplexes containing diagonal TTTT loops.     

Guanine residues in d(T2AG3) and d(T2G4) form parallel-stranded potassium cation stabilized G-quadruplexes with anti glycosidic torsion angles in solution.

We report below on proton NMR studies of the G-quadruplex structure formed by the human telomere sequence d(T2AG3) and the tetrahymena telomere sequence d(T2G4) in K cation containing solution. We observe well-resolved proton NMR spectra corresponding to a G-quadruplex monomer conformation predominant at 50 mM K cation concentration and a G-quadruplex dimer conformation predominant at 300 mM K cation concentration. By contrast, d(T2AG3T) and d(T2G4T) form only the G-quadruplex monomer structures independent of K cation concentration as reported previously [Sen, D., & Gilbert, W. (1992) Biochemistry 31, 65-70]. We detect well-resolved resonances for the exchangeable guanine imino and amino protons involved in G-tetrad formation with the hydrogen-bonded and exposed amino protons separated by up to 3.5 ppm. The observed NOEs between the amino and H8 protons on adjacent guanines within individual G-tetrads support the Hoogsteen pairing alignment around the tetrad. The imino protons of the internal G-tetrads exchange very slowly with solvent H2O in the d(T2AG3) and d(T2G4) quadruplexes. The nature and intensity of the observed NOE patterns establish formation of parallel-stranded right-handed G-quadruplexes with all anti guanine glycosidic torsion angles. A model for the parallel-stranded G-quadruplex is proposed which is consistent with the experimental NOE data on the d(T2AG3) and d(T2G4) quadruplexes in solution.(ABSTRACT TRUNCATED AT 250 WORDS)     

The physical basis of nucleic acid base stacking in water

It has been argued that the stacking of adenyl groups in water must be driven primarily by electrostatic interactions, based upon NMR data showing stacking for two adenyl groups joined by a 3-atom linker but not for two naphthyl groups joined by the same linker. In contrast, theoretical work has suggested that adenine stacking is driven primarily by nonelectrostatic forces, and that electrostatic interactions actually produce a net repulsion between adenines stacking in water. The present study provides evidence that the experimental data for the 3-atom-linked bis-adenyl and bis-naphthyl compounds are consistent with the theory indicating that nonelectrostatic interactions drive adenine stacking. First, a theoretical conformational analysis is found to reproduce the observed ranking of the stacking tendencies of the compounds studied experimentally. A geometric analysis identifies two possible reasons, other than stronger electrostatic interactions, why the 3-atom-linked bis-adenyl compounds should stack more than the bis-naphthyl compounds. First, stacked naphthyl groups tend to lie further apart than stacked adenyl groups, based upon both quantum calculations and crystal structures. This may prevent the bis-naphthyl compound from stacking as extensively as the bis-adenyl compound. Second, geometric analysis shows that more stacked conformations are sterically accessible to the bis-adenyl compound than to the bis-naphthyl compound because the linker is attached to the sides of the adenyl groups, but to the ends of the naphthyl groups. Finally, ab initio quantum mechanics calculations and energy decompositions for relevant conformations of adenine and naphthalene dimers support the view that stacking in these compounds is driven primarily by nonelectrostatic interactions. The present analysis illustrates the importance of considering all aspects of a molecular system when interpreting experimental data, and the value of computer models as an adjunct to chemical intuition.