Conformational polymorphism in G-tetraplex structures: strand reversal by base flipover or sugar flipover

The Interactions of H2TMPyP, Analogues and Its Metal Complexes with DNA G-Quadruplexes-An Overview

The evidence that telomerase is overexpressed in almost 90% of human cancers justifies the proposal of this enzyme as a potential target for anticancer drug design. The inhibition of telomerase by quadruplex stabilizing ligands is being considered a useful approach in anticancer drug design proposals. Several aromatic ligands, including porphyrins, were exploited for telomerase inhibition by adduct formation with G-Quadruplex (GQ). 5,10,15,20-Tetrakis(N-methyl-4-pyridinium)porphyrin (H₂TMPyP) is one of the most studied porphyrins in this field, and although reported as presenting high affinity to GQ, its poor selectivity for GQ over duplex structures is recognized. To increase the desired selectivity, porphyrin modifications either at the peripheral positions or at the inner core through the coordination with different metals have been handled. Herein, studies involving the interactions of TMPyP and analogs with different DNA sequences able to form GQ and duplex structures using different experimental conditions and approaches are reviewed. Some considerations concerning the structural diversity and recognition modes of G-quadruplexes will be presented first to facilitate the comprehension of the studies reviewed. Additionally, considering the diversity of experimental conditions reported, we decided to complement this review with a screening where the behavior of H₂TMPyP and of some of the reviewed metal complexes were evaluated under the same experimental conditions and using the same DNA sequences. In this comparison under unified conditions, we also evaluated, for the first time, the behavior of the Ag^II complex of H₂TMPyP. In general, all derivatives showed good affinity for GQ DNA structures with binding constants in the range of 10⁶–10⁷ M⁻¹ and ligand-GQ stoichiometric ratios of 3:1 and 4:1. A promising pattern of selectivity was also identified for the new Ag^II derivative.

Targeting unimolecular G-quadruplex nucleic acids: a new paradigm for the drug discovery?

INTRODUCTION

G-quadruplexes (G4s) are targets of great interest because of their roles in crucial biological processes, such as aging and cancer. G4s are based on the formation of G-quartets, stabilised by Hoogsteen-type hydrogen bonds and by interaction with cations between the tetrads. These biologically relevant conformations were first discovered in eukaryotic chromosomal telomeric DNA, but have also been found in the proximal location of promoters in a number of human genes. Therefore, the extensive analysis of an intriguing target could move towards the rational drug design of new selective anticancer agents.

AREAS COVERED

The authors review G4 structural characterisation, with detailed insight related to the polymorphism issue. The authors describe the topologically distinct G4 structural forms and the factors involved in their interconversion mechanisms, such as the sequence of the oligonucleotides, the strand stoichiometry and orientation, the syn-anti conformation of the guanine glycosidic bonds and the G4 loop types and the environmental factors. Furthermore, the authors report several studies related to folding and unfolding kinetic profiles in order to understand the conformational view of monomolecular G4 formations.

EXPERT OPINION

G4 unimolecular nucleic acids can be considered as valid targets for the rational drug development of novel anticancer agents. Structural biology represents an essential link between the biology and medicinal chemistry knowledge in this field. In silico methods have already been demonstrated to be useful, especially if well integrated with biophysical tests. If this proves successful, the G4-targeting paradigm could also be extended to drug discovery beyond neoplastic pathologies.

IFI16 and NM23 bind to a common DNA fragment both in the P53 and the cMYC gene promoters

In the melanoma M14 cell line, we found that the antimetastatic protein NM23/nucleoside diphosphate kinase binds to the promoters of the oncogene cMYC and of P53, a gene often mutated in human cancer (Cervoni et al. [2006] J. Cell. Biochem. 98:421-428). In a further study, we find now that IFI16, a transcriptional repressor, in both promoters binds to the G-rich fragment that also binds NM23/NDPK. These fragments possess non-B DNA structures. Moreover, by sequential chromatin immunoprecipitation (re-ChIP) we show that the two proteins (IFI16 and NM23/NDPK) are simultaneously bound in vivo to the same DNA fragments. Since P53 stimulates apoptosis and inhibits cellular growth, and cMYC promotes cell growth and, in several instances, also apoptosis, the presence of NM23 and IFI16 on the same DNA fragments suggests their common involvement in the reduced development of some tumors.

Molecular dynamics simulations of Guanine quadruplex loops: advances and force field limitations

A computational analysis of d(GGGGTTTTGGGG)(2) guanine quadruplexes containing either lateral or diagonal four-thymidine loops was carried out using molecular dynamics (MD) simulations in explicit solvent, locally enhanced sampling (LES) simulations, systematic conformational search, and free energy molecular-mechanics, Poisson Boltzmann, surface area (MM-PBSA) calculations with explicit inclusion of structural monovalent cations. The study provides, within the approximations of the applied all-atom additive force field, a qualitatively complete analysis of the available loop conformational space. The results are independent of the starting structures. Major conformational transitions not seen in conventional MD simulations are observed when LES is applied. The favored LES structures consistently provide lower free energies (as estimated by molecular-mechanics, Poisson Boltzmann, surface area) than other structures. Unfortunately, the predicted optimal structure for the diagonal loop arrangement differs substantially from the atomic resolution experiments. This result is attributed to force field deficiencies, such as the potential misbalance between solute-cation and solvent-cation terms. The MD simulations are unable to maintain the stable coordination of the monovalent cations inside the diagonal loops as reported in a recent x-ray study. The optimal diagonal and lateral loop arrangements appear to be close in energy although a proper inclusion of the loop monovalent cations could stabilize the diagonal architecture.

Crystal structure of a bulged RNA tetraplex at 1.1 a resolution: implications for a novel binding site in RNA tetraplex

Bulges are an important structural motif in RNA and can be used as recognition and interaction sites in RNA-protein interaction and RNA-RNA interaction. Here we report the first crystal structure of a bulged RNA tetraplex at 1.1 A resolution. The hexamer r(U)(BrdG)r(UGGU) forms a parallel tetraplex with the uridine sandwiched by guanines bulging out. The bulged uridine adopts the syn glycosidic conformation and its O2 and N3 atoms face outwards, serving as an effective recognition and interaction site. The bulge formation both widens the groove width and changes the groove hydrogen-bonding pattern on its 5' side. However, the bulge does not make any bends or kinks in the tetraplex structure. The present study demonstrates the dramatic difference between uridine and guanine in forming tetraplex structure. In addition, both G(syn) tetrad and G(anti) tetrad have been observed. They display the same base-pairing pattern and similar C1'-C1' distance but different hydrogen-bonding patterns in the groove.

An eight-stranded helical fragment in RNA crystal structure: implications for tetraplex interaction

Multistranded helical structures in nucleic acids play various functions in biological processes. Here we report the crystal structure of a hexamer, rU(BrdG)r(AGGU),at 1.5 A resolution containing a structural complex of an alternating antiparallel eight-stranded helical fragment that is sandwiched in two tetraplexes. The octaplex is formed by groove binding interaction and base tetrad intercalation between two tetraplexes. Two different forms of octaplexes have been proposed, which display different properties in interaction with proteins and nucleic acids. Adenines form a base tetrad in the novel N6-H em leader N3 conformation and further interact with uridines to form an adenine-uridine octad in the reverse Hoogsteen pairing scheme. The conformational flexibility of adenine tetrad indicates that it can optimize its conformation in different interactions.

Modelling studies on neurodegenerative disease-causing triplet repeat sequences d(GGC/GCC)n and d(CAG/CTG)n

Model building and molecular mechanics studies have been carried out to examine the potential structures for d(GGC/GCC)5 and d(CAG/CTG)5 that might relate to their biological function and association with triplet repeat expansion diseases. Model building studies suggested that hairpin and quadruplex structures could be formed with these repeat sequences. Molecular mechanics studies have demonstrated that the hairpin and hairpin dimer structures of triplet repeat sequences formed by looping out of the two strands are as favourable as the corresponding B-DNA type hetero duplex structures. Further, at high salt condition, Greek key type quadruplex structures are energetically comparable with hairpin dimer and B-DNA type duplex structures. All tetrads in the quadruplex structures are well stacked and provide favourable stacking energy values. Interestingly, in the energy minimized hairpin dimer and Greek key type quadruplex structures, all the bases even in the non-G tetrads are cyclically hydrogen bonded, even though the A, C and T-tetrads were not hydrogen bonded in the starting structures.

A nanosecond molecular dynamics study of antiparallel d(G)7 quadruplex structures: effect of the coordinated cations

Nanosecond scale molecular dynamics simulations have been performed on antiparallel Greek key type d(G7) quadruplex structures with different coordinated ions, namely Na+ and K+ ion, water and Na+ counter ions, using the AMBER force field and Particle Mesh Ewald technique for electrostatic interactions. Antiparallel structures are stable during the simulation, with root mean square deviation values of approximately 1.5 A from the initial structures. Hydrogen bonding patterns within the G-tetrads depend on the nature of the coordinated ion, with the G-tetrad undergoing local structural variation to accommodate different cations. However, alternating syn-anti arrangement of bases along a chain as well as in a quartet is maintained through out the MD simulation. Coordinated Na+ ions, within the quadruplex cavity are quite mobile within the central channel and can even enter or exit from the quadruplex core, whereas coordinated K+ ions are quite immobile. MD studies at 400K indicate that K+ ion cannot come out from the quadruplex core without breaking the terminal G-tetrads. Smaller grooves in antiparallel structures are better binding sites for hydrated counter ions, while a string of hydrogen bonded water molecules are observed within both the small and large grooves. The hydration free energy for the K+ ion coordinated structure is more favourable than that for the Na+ ion coordinated antiparallel quadruplex structure.

G-quadruplex DNA structures--variations on a theme

To be functional, nucleic acids need to adopt particular three-dimensional structures. For a long time DNA was regarded as a rigid and passive molecule with the sole purpose to store genetic information, but experimental data has now accumulated that indicates the full dynamic repertoire of this macromolecule. During the last decade, four-stranded DNA structures known as G-quadruplexes, or DNA tetraplexes, have emerged as a three-dimensional structure of special interest. Motifs for the formation of G-quadruplex DNA structures are widely dispersed in eukaryotic genomes, and are abundant in regions of biological significance, for example, at telomeres, in the promoters of many important genes, and at recombination hotspots, to name but a few in man. Here I explore the plethora of G-quadruplex DNA structures, and discuss their possible biological functions as well as the proteins that interact with them.

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.

Molecular dynamics simulations on parallel and antiparallel C.G*G triplexes

Molecular dynamics (MD) studies have been carried out on the Hoogsteen hydrogen bonded parallel and the reverse Hoogsteen hydrogen bonded antiparallel C.G*G triplexes. Earlier, the molecular mechanics studies had shown that the parallel structure was energetically more favourable than the antiparallel structure. To characterize the structural stability of the two triplexes and to investigate whether the antiparallel structure can transit to an energetically more favourable structure, due to the local fluctuations in the structure during the MD simulation, the two structures were subjected to 200ps of constant temperature vacuum MD simulations at 300K. Initially no constraints were applied to the structures and it was observed that for the antiparallel triplex, the structure showed a large root mean square deviation from the starting structure within the first 12ps and the N4-H41--O6 hydrogen bond in the WC duplex got distorted due to a high propeller twist and a moderate increase in the opening angle in the basepairs. Starting from an initial value of 30 degrees , helical twist of the average structure from this simulation had a value of 36 degrees , while the parallel structure stabilized at a twist of 33 degrees. In spite of the hydrogen bond distortions in the antiparallel triplex, it was energetically comparable to the parallel triplex. To examine the structural characteristics of an undistorted structure, another MD simulation was performed on the antiparallel triplex by constraining all the hydrogen bonds. This structure stabilized at an average twist of 33 degrees. In the course of the dynamics though the energy of the molecule - compared to the initial structure - improved, it did not become comparable to the parallel structure. Energy minimization studies performed in the presence of explicit water and counterions also showed the two structures to be equally favourable energetically. Together these results indicate that the parallel C.G*G triplex with Hoogsteen hydrogen bonds also represents a stereochemically and energetically favourable structure for this class of triplexes.

Sequence-independent recombination triple helices: a molecular dynamics study

Recent experimental studies have shown that the Rec-A mediated homologous recombination reaction involves a triple helical intermediate, in which the third strand base forms hydrogen bonds with both the bases in the major groove of the Watson-Crick duplex. Such 'mixed' hydrogen bonds allow formation of sequence independent triplexes. DNA triple helices involving 'mixed' hydrogen bonds have been studied, using model building, molecular mechanics (MM) and molecular dynamics (MD). Models were built for a triplex comprising all four possible triplets viz., G.C*C, C.G*G, A.T*T and T.A*A. To check the stability of all the 'mixed' hydrogen bonds in such triplexes and the conformational preferences of such triplex structures, MD studies were carried out starting from two structures with 30 degrees and 36 degrees twist between the basepairs. It was observed that though the two triplexes converged towards a similar structure, the various hydrogen bonds between the WC duplex and the third strand showed differential stabilities. An MD simulation with restrained hydrogen bonds showed that the resulting structure was stable and remained close to the starting structure. These studies help us in defining stable hydrogen bond geometries involving the third strand and the WC duplex. It was observed that in the C.G*G triplets the N7 atom of the second strand is always involved in hydrogen bonding. In the G.C*C triplets, either N3 or O2 in the third strand cytosine can interchangeably act as a hydrogen bond acceptor.

A vibrational spectroscopic study of the metastable form of associated polyinosinic acid

We have studied by Raman and ir spectroscopy the metastable complex formed by the self-association of polyinosinic acid in aqueous solution. The complex is easily prepared by quickly cooling to ca. 0 degrees C a warm solution of the polyribonucleotide to which a small amount of rubidium salt has been added. Upon heating, this metastable form melts cooperatively near 13 degrees C, well below the dissociation temperature of a stable four-stranded complex, which occurs at 47 degrees C in the same conditions. The presence of several components in the stretching-mode region of the carbonyl groups in the vibrational spectra of the metastable complex suggests that it also has a parallel four-stranded structure. The difference in structure between the two forms is believed to be caused by the presence of fewer metal ions in the central channel of the metastable complex, in agreement with conclusions reached in previous investigations. The Raman spectra further show that the ribose units in the metastable form have a C3'-endo conformation, in contrast with the stable form, for which we have previously suggested a mixed C2'-endo/C3'-endo conformation.

Structural polymorphism in d(T)12.d(A)12*d(T)12 triple helices

DNA triple helices containing two thymine strands and one adenine strand have been studied, using model building followed by energy minimisation, for different orientations of the third strand resulting from variation in the hydrogen bonding between the Watson-Crick duplex and the third strand and the glycosidic torsion angle in the third strand. Our results show that the structure with a parallel orientation of the third strand, in which the third strand base forms Hoogsteen hydrogen bonds with the adenine base in the Watson-Crick duplex, is energetically the most favourable. An antiparallel orientation of the third strand is also possible, in which the third strand base hydrogen bonds to both the bases in the Watson-Crick duplex. This structure is energetically comparable to the parallel structure. For the parallel triplex a 200ps molecular dynamics simulation starting from two different starting structures indicates that at 300K significant structural heterogeneity exists in this triplex structure. The results are compared with existing structural data on this class of triplexes derived from theoretical and NMR techniques.

Chain folding and A:T pairing in human telomeric DNA: a model-building and molecular dynamics study

The various types of chain folding and possible intraloop as well as interloop base pairing in human telomeric DNA containing d(TTAG3) repeats have been investigated by model-building, molecular mechanics, and molecular dynamics techniques. Model-building and molecular mechanics studies indicate that it is possible to build a variety of energetically favorable folded-back structures with the two TTA loops on same side and the 5' end thymines in the two loops forming TATA tetrads involving a number of different intraloop as well as interloop A:T pairing schemes. In these folded-back structures, although both intraloop and interloop Watson-Crick pairing is feasible, no structure is possible with interloop Hoogsteen pairing. MD studies of representative structures indicate that the guanine-tetraplex stem is very rigid and, while the loop regions are relatively much more flexible, most of the hydrogen bonds remain intact throughout the 350-ps in vacuo simulation. The various possible TTA loop structures, although they are energetically similar, have characteristic inter proton distances, which could give rise to unique cross-peaks in two-dimensional nuclear Overhauser effect spectroscopy (NOESY) experiments. These folded-back structures with A:T pairings in the loop region help in rationalizing the data from chemical probing and other biochemical studies on human telomeric DNA.

Conformational polymorphism in telomeric structures: loop orientation and interloop pairing in d(G4TnG4)

Sequence repeats constituting the telomeric regions of chromosomes are known to adopt a variety of unusual structures, consisting of a G tetraplex stem and short stretches of thymines or thymines and adenines forming loops over the stem. Detailed model building and molecular mechanics studies have been carried out for these telomeric sequences to elucidate different types of loop orientations and possible conformations of thymines in the loop. The model building studies indicate that a minimum of two thymines have to be interspersed between guanine stretches to form folded-back structures with loops across adjacent strands in a G tetraplex (both over the small as well as large groove), while the minimum number of thymines required to build a loop across the diagonal strands in a G tetraplex is three. For two repeat sequences, these hairpins, resulting from different types of folding, can dimerize in three distinct ways--i.e., with loops across adjacent strands and on same side, with loops across adjacent strands and on opposite sides, and with loops across diagonal strands and on opposite sides--to form hairpin dimer structures. Energy minimization studies indicate that all possible hairpin dimers have very similar total energy values, though different structures are stabilized by different types of interactions. When the two loops are on the same side, in the hairpin dimer structures of d(G4TnG4), the thymines form favorably stacked tetrads in the loop region and there is interloop hydrogen bonding involving two hydrogen bonds for each thymine-thymine pair. Our molecular mechanics calculations on various folded-back as well as parallel tetraplex structures of these telomeric sequences provide a theoretical rationale for the experimentally observed feature that the presence of intervening thymine stretches stabilizes folded-back structures, while isolated stretches of guanines adopt a parallel tetraplex structure.

Monoclonal antibodies targeted to alpha-oligonucleotides. Characterisation and application in nucleic acid detection

The aim of the present study was to test the antigenicity of alpha-deoxyribonucleotides in order to develop a new tool for the detection of nucleic acid sequences for use in diagnostic applications. We describe four monoclonal antibodies (Mabs) which recognize alpha-deoxyribonucleotides. Two were raised against a poly(alpha-dT) sequence and specifically recognized the alpha-dT nucleotide. Two were raised against a sequence containing all four common nucleotides as alpha-nucleotides and, surprisingly, only recognized the alpha-dG nucleotide. For all four Mabs, no cross reactivity was observed with beta-oligonucleotides. These Mabs were reactive with alpha-oligonucleotide sequences whether these sequences were single-stranded or hybridized to DNA or RNA. The four Mabs were tested in a sandwich hybridization assay that consisted of an alpha-oligonucleotide (for target sequence recognition), one of the four Mabs (for recognition of the hybridized alpha-oligonucleotide), and goat anti-mouse antibody conjugated to horse radish peroxidase (HRP) (for detection). One of the monoclonal antibodies, Mab 2E11D7, was directly conjugated to HRP and used in sandwich hybridization to detect PCR fragments of HPV 18 DNA. The sensitivity of this reaction was 1 pg of plasmid DNA containing the HPV 18 fragment. The specificity of the detection was demonstrated using HPV 6/11 and 16 DNA sequences.

A vibrational spectroscopic study of the self-association of polyinosinic acid and polyguanylic acid in aqueous solution

We have studied by Raman and ir spectroscopy the structure of self-associated polyinosinic acid and polyguanylic acid in aqueous solution. The results are consistent with the formation of a four-stranded complex, which melts cooperatively near 60 degrees C in the case of poly(I) in the presence of K+ ions. The conformation of the ribose in both systems is mixed C2'-endo/C3'-endo, giving a structure that is intermediate between the extremes proposed previously from x-ray diffraction studies. Characteristic Raman bands for the C2'-endo ribose conformation in polyribonucleotides are identified. The four-stranded structure of poly(I) appears to be very flexible, with approximately 15% of the tetrameric segments being disrupted and approximately 30% of the ribose units adopting a disordered conformation prior to melting. This disordering process increases to approximately 75% above the melting transition, with the remaining approximately 25% of the ribose units keeping an ordered C2'-endo or C3'-endo conformation.