Stability and cations coordination of DNA and RNA 14-mer G-quadruplexes: a multiscale computational approach.

Molecular dynamics simulations have been used to study the differences between two DNA and RNA 14-mer quadruplexes of analogous sequences. Their structures present a completely different fold: DNA forms a bimolecular quadruplex containing antiparallel strands and diagonal loops; RNA forms an intrastrand parallel quadruplex containing a G-tetrad and an hexad, which dimerizes by hexad stacking. We used a multiscale computational approach combining classical Molecular dynamics simulations and density functional theory calculations to elucidate the difference in stability of the 2-folds and their ability in coordinating cations. The presence of 2'-OH groups in the RNA promotes the formation of a large number of intramolecular hydrogen bonds that account for the difference in fold and stability of the two 14-mers. We observe that the adenines in the RNA quadruplex play a key role in conserving the geometry of the hexad. We predict the cation coordination mode of the two quadruplexes, not yet observed experimentally, and we offer a rationale for the corresponding binding energies involved.

Structure, interactions and effects on activity of the 5'-terminal region of human telomerase RNA.

Telomerase is an enzyme that catalyzes addition of telomeric repeat sequences to the 3'-termini of eukaryotic chromosome DNA. The catalytic core of telomerase consists of a protein component, telomerase reverse transcriptase (TERT), for the catalysis and an RNA component, telomerase RNA (TR), containing the template for the sequence. Human telomerase RNA (hTR) consists of 451 nucleotides (nt) and contains consecutive G-stretches in the 5'-terminal region. We examined the effects of the 5'-terminal sequence (nt 1-17) in hTR, which is assumed to be a single-stranded region (region 1), on interaction and telomerase activity in vitro. Mutation and binding experiments for hTR and its variants suggest that region 1 has repressive effects on telomerase activity by interaction with the region(s) in the 3'-half part. We prepared various hTR variants with mutations in region 1 and two possible target regions (region 2: nt 229-244; region 3: nt 284-297). Studies on these variants showed that region 1 can interact with regions 2 and 3 and the interactions between regions 1 and 3 may contribute to the repressive effects of region 1. We found that a mutation in region 2 markedly enhances telomerase activity. We also found that some deletion and sequence mutations in region 1 enhance the activity.

The RNA ligands for mouse proline-rich RNA-binding protein (mouse Prrp) contain two consensus sequences in separate loop structure.

Mouse proline-rich RNA-binding protein (mPrrp) is a mouse ortholog of Xenopus Prrp, which binds to a vegetal localization element (VLE) in the 3'-untranslated region (3'-UTR) of Vg1 mRNA and is expected to be involved in the transport and/or localization of Vg1 mRNA to the vegetal cortex of oocytes. In mouse testis, mPrrp protein is abundantly expressed in the nuclei of pachytene spermatocytes and round spermatids, and shifts to the cytoplasm in elongating spermatids. To gain an insight into the function of mPrrp in male germ cells, we performed in vitro RNA selection (SELEX) to determine the RNA ligand sequence of mPrrp. This analysis revealed that many of the selected clones contained both of two conserved elements, AAAUAG and GU1-3AG. RNA-binding study on deletion mutants and secondary structure analyses of the selected RNA revealed that a two-loop structure containing the conserved elements is required for high-affinity binding to mPrrp. Furthermore, we found that the target mRNAs of Xenopus Prrp contain intact AAAUAG and GU1-3AG sequences in the 3'-UTR, suggesting that these binding sequences are shared by Prrps of Xenopus and mouse.

Method for direct discrimination of intra- and intermolecular hydrogen bonds, and characterization of the G(:A):G(:A):G(:A):G heptad, with scalar couplings across hydrogen bonds.

Discrimination of intra- and intermolecular hydrogen bonds in a symmetric multimer has not been accomplished yet, although such discrimination would provide a crucial basis for construction of the multimeric architecture of nucleic acids by NMR. We have developed a direct and unambiguous method for such discrimination involving the use of scalar couplings across hydrogen bonds. The method has been validated with a symmetric dimer of d(GGGCTTTTGGGC), for which the structure including both intra- and intermolecular hydrogen bonds was already reported. This has demonstrated that our method can clearly discriminate these two kinds of hydrogen bonds. Then, the method was applied to a symmetric dimer of d(GGAGGAGGAGGA) and has provided decisive information on its multimeric architecture. Additionally, the values for scalar couplings across hydrogen bonds for G:G and G:A base pairs in the G(:A):G(:A):G(:A):G heptad formed by d(GGAGGAGGAGGA) were determined for the first time. This determination has provided an insight into the nature of the heptad.

Imino proton NMR analysis of HDV ribozymes: nested double pseudoknot structure and Mg2+ ion-binding site close to the catalytic core in solution.

Minimized trans-acting HDV ribozyme systems consisting of three (Rz-3) and two (Rz-2) RNA strands were prepared and their folding conformations were analyzed by NMR spectroscopy. The guanosine residues in one of the enzyme components of Rz-3 were labeled with 13C and 15N. Imino proton signals were assigned by analysis of NOESY and HSQC spectra. The results are consistent with the nested double pseudoknot model, which contains novel base pairs (P1.1), as observed in the crystal structure of a genomic HDV ribozyme. The NOE connectivities suggest an additional G:G pair at the bottom of P1.1 and at the top of P4. The effects of temperature and Mg2+ ions on base pairs for Rz-3 were examined. The temperature variation experiment on Rz-3 showed that P3 is the most stable and that P1.1 is as stable as P1 and P2. The imino proton signals of the G:U pair at the bottom of P1 and the top of P1.1, which are close to the cleavage site, showed the largest changes upon Mg2+ titration of Rz-3. The results suggest that the catalytic Mg2+ ion binds to the pocket formed by P1 and L3.

Structural basis of the highly efficient trapping of the HIV Tat protein by an RNA aptamer.

An RNA aptamer containing two binding sites exhibits extremely high affinity to the HIV Tat protein. We have determined the structure of the aptamer complexed with two argininamide molecules. Two adjacent U:A:U base triples were formed, which widens the major groove to make space for the two argininamide molecules. The argininamide molecules bind to the G bases through hydrogen bonds. The binding is stabilized through stacking interactions. The structure of the aptamer complexed with a Tat-derived arginine-rich peptide was also characterized. It was suggested that the aptamer structure is similar for both complexes and that the aptamer interacts with two different arginine residues of the peptide simultaneously at the two binding sites, which could explain the high affinity to Tat.

Addition of an extra substrate binding site and partial destabilization of stem structures in HDV ribozyme give rise to high sequence-specificity for its target RNA.

Because the substrate binding site (P1) of HDV ribozyme consists of only seven nucleotides, cleavage of undesired RNA is likely to occur when applied for a specific long RNA target such as mRNA. To overcome this problem, we designed modified trans-acting HDV ribozymes with an extra substrate-binding site (P5) in addition to the original binding site (P1). By inserting an additional seven base-pair stem (P5 stem) into the J1/2 single-stranded region of the ribozyme core system and partial destabilization of the P2 or P4 stem, we succeeded in preparation of new HDV ribozymes that can cleave the target RNA depending on the formation of P5 stem. Moreover, the ribozyme with a six-nucleotide P1 site was able to distinguish the substrate RNA with a complete match from that with a single mismatch in the P1 region. These results suggest that the HDV ribozyme system is useful for the application in vivo.

Structural property of DNA that migrates faster in gel electrophoresis, as deduced by CD spectroscopy.

Bent DNAs are known to migrate slower than ordinary DNA in non-denaturing polyacrylamide gel electrophoresis. In contrast, several satellite DNAs have been shown to migrate fast. The structural property that causes the fast migration, however, is not clarified so far on molecular basis. We have investigated the structural property of a satellite DNA, which contains consecutive purine sequences and migrates faster in gel, by CD spectroscopy. Partial formation of an A-form-like structure has been suggested. Reduction in DNA length due to the formation of the A-form-like structure may be responsible for the fast migration. The pronounced rigidity of DNA may also contribute to the behavior.

A dimeric RNA quadruplex architecture comprised of two G:G(:A):G:G(:A) hexads, G:G:G:G tetrads and UUUU loops.

Using CD and NMR, we determined the structure of an RNA oligomer, r(GGAGGUUUUGGAGG) (R14), comprising two GGAGG segments joined by a UUUU segment. A modified quadruplex structure was observed for r(GGAGGUUUUGGAGG) in solution even in the absence of K(+). An unusually stable dimeric RNA quadruplex architecture formed from two strands of r(GGAGGUUUUGGAGG) at low K(+) concentration is reported here. In each strand of r(GGAGGUUUUGGAGG), two sets of successive turns in the GGAGG segments and turns at both ends of the UUUU loops drive four G-G steps to align in a parallel manner, a core with two stacked G-tetrads being formed. Two adenine bases bind to two edges of one G:G:G:G tetrad through the sheared G:A mismatch augmenting the tetrad into a G:G(:A):G:G(:A) hexad. Thus, one molecule of r(GGAGGUUUUGGAGG) folds into a modified quadruplex comprising a G:G:G:G tetrad, a UUUU double-chain reversal loop and a G:G(:A):G:G(:A) hexad. Two such molecules further associate by stacking through the dimeric hexad-hexad interface with a rotational symmetry. The ribose rings of most nucleotides take S (close to C2'-endo) puckering, which is unusual for an RNA. K(+) can increase the stability of this quadruplex structure; the number of bound K(+) was estimated from the results of the titration experiment. Besides G:G and G:A mismatches, a network of hydrogen bonds including O4'-NH(2) and C-H..O hydrogen bonds, and the extensive base stacking contribute to the high thermodynamic stability of R14. Our results could provide the stereochemical and thermodynamic basis for elucidating the biological role of the GGAGG-containing RNA segments abundantly existing in various RNAs. Relevance to quadruplex-mediated mRNA-FMRP binding and HIV-1 genome RNA dimerization is discussed.

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.

Studies on RNA-protein interactions and enzyme activity by mutational analysis of telomerase.

The protein component of telomerase (TERT: telomerase reverse transcriptase) contains a domain (RID2: RNA interaction domain 2), which interacts with the RNA component of telomerase (TR: telomerase RNA). RID2 includes a telomerase-specific motif (T-motif), which is highly conserved among species. But, the role of T-motif is not clearly understood. To elucidate the role of T-motif in the telomerase activity and the TERT-TR interaction, we designed two systems including full-length human TERT (hTERT) proteins for telomerase activity and RID2 domain peptides for binding activity. Variants of hTERT, in which a highly conserved amino acid residue in T-motif was replaced with an alanine residue, were expressed in the in vitro translation system of rabbit reticulocyte lysate. We examined telomerase activity of the hTERT variants by the stretch-PCR method. T564A, E565A, and W581A variants did not show remarkably decreased activity suggesting that these amino acid residues are not important for the enzymatic activity, although they are conserved with 89-100% identity. Next, we tried construction of an experimental system for examination of RNA-binding activity of RID2 variants with the same mutations in T-motif. The wild-type RID2 peptide with His-tags was prepared by expression in an E. coli system and purified with an affinity column. The binding activity to human TR (hTR) was examined by the EMSA (electrophoretic mobility shift assay) method. The RID2 (aa 300-617) peptide showed binding to the 3'half molecule of hTR but not to the 5'-half molecule. The RID2 variants are being prepared and analyzed in this system.

Studies on activity and interactions of human telomerase.

It is reported that a single stranded region at 5'-terminus of human telomerase RNA (hTR) affects telomerase activity, though the sequences of the region are not conserved among vertebrate species. The mechanism of this phenomenon is not known. We examined binding affinity of an RNA oligomer (R17), which corresponds to the single stranded region (1-17) at 5'-terminus of hTR, to deletion variants of hTR. R17 showed higher affinity toward the 3'-half part of hTR than that for the 5'-half part. We chose two regions of the 3'-half part, where R17 is assumed to bind, and prepared variants of hTR in the single stranded region and the selected regions. The interactions and telomerase activities of these variants were examined. We also prepared an RNA interaction domain of the protein component (hTERT) of telomerase and its variants and tried to elucidate influence on interactions with hTR.

PSPC1, NONO, and SFPQ are expressed in mouse Sertoli cells and may function as coregulators of androgen receptor-mediated transcription.

In Sertoli cells of testis, androgen receptor-regulated gene transcription plays an indispensable role in maintaining spermatogenesis. Androgen receptor activity is modulated by a number of coregulators which are associated with the androgen receptor. Non-POU-domain-containing, octamer binding protein (NONO), a member of the DBHS-containing proteins, complexes with androgen receptor and functions as a coactivator for the receptor. Paraspeckle protein 1 alpha isoform (PSPC1, previously known as PSP1) and Splicing factor, proline- and glutamine-rich (SFPQ, previously known as PSF), other members of the DBHS-containing proteins, are also found in androgen receptor complexes, suggesting that these DBHS-containing proteins may cooperatively regulate androgen receptor-mediated gene transcription. We demonstrated that PSPC1, NONO, and SFPQ are coexpressed in Sertoli cell line TTE3 and interact reciprocally. The effect of the DBHS-containing proteins on the transcriptional activity was assessed using the construct containing androgen-responsive elements followed by a luciferase gene. The results showed that all the DBHS-containing proteins activate androgen receptor-mediated transcription, and PSPC1 is the most effective coactivator among them. Furthermore, we confirmed the presence of PSPC1, NONO, and SFPQ proteins in Sertoli cells of adult mouse testis sections. These observations suggest that PSPC1, NONO, and SFPQ form complexes with each other in Sertoli cells and may regulate androgen receptor-mediated transcriptional activity.

Effects on telomerase activity of the 5'-terminal region of human telomerase RNA.

Human telomerase RNA (hTR) is assumed to contain a 17-nt segment in a single-stranded state at the 5'-terminus. There are four stretches of consecutive G residues, which are apt to form a quadruplex structure, in the segment. It is reported that deletion of some part of this region enhances telomerase activity when the core telomerase enzyme is reconstituted by addition of telomerase reverse transcriptase. To elucidate the reason for such effects, we constructed hTR mutants with deletion of the segment (hTRdelta20) and with the segment containing four G-to-A displacement to interrupt the G-stretch sequences (hTR17A) and their activity was compared with that of the wild-type hTR (hTRW). hTRdelta20 showed much higher activity than the other while hTR17A showed activity similar to that of hTRW. This result suggests that the lower activity for full-length hTR is not due to putative quadruplex formation.

Structure of hnRNP D complexed with single-stranded telomere DNA and unfolding of the quadruplex by heterogeneous nuclear ribonucleoprotein D.

Heterogeneous nuclear ribonucleoprotein D, also known as AUF1, has two DNA/RNA-binding domains, each of which can specifically bind to single-stranded d(TTAGGG)n, the human telomeric repeat. Here, the structure of the C-terminal-binding domain (BD2) complexed with single-stranded d(TTAGGG) determined by NMR is presented. The structure has revealed that each residue of the d(TAG) segment is recognized by BD2 in a base-specific manner. The interactions deduced from the structure have been confirmed by gel retardation experiments with mutant BD2 and DNA. It is known that single-stranded DNA with the telomeric repeat tends to form a quadruplex and that the quadruplex has an inhibitory effect on telomere elongation by telomerase. This time it is revealed that BD2 unfolds the quadruplex of such DNA upon binding. Moreover, the effect of BD2 on the elongation by telomerase was examined in vitro. These results suggest the possible involvement of heterogeneous nuclear ribonucleoprotein D in maintenance of the telomere 3'-overhang either through protection of a single-stranded DNA or destabilization of the potentially deleterious quadruplex structure for the elongation by telomerase.

Intramolecular higher order packing of parallel quadruplexes comprising a G:G:G:G tetrad and a G(:A):G(:A):G(:A):G heptad of GGA triplet repeat DNA.

GGA triplet repeats are widely dispersed throughout eukaryotic genomes and are frequently located within biologically important regions such as gene regulatory regions and recombination hot spot sites. We determined the structure of d(GGA)4 (12-mer) under physiological conditions and founded the formation of an intramolecular parallel quadruplex for the first time. Later, a similar architecture to that of the intramolecular parallel quadruplex was found for a telomere DNA in the crystalline state. Here, we have determined the structure of d(GGA)8 (24-mer) under physiological conditions. Two intramolecular parallel quadruplexes comprising a G:G:G:G tetrad and a G(:A):G(:A):G(:A):G heptad are formed in d(GGA)8. These quadruplexes are packed in a tail-to-tail manner. This is the first demonstration of the intramolecular higher order packing of quadruplexes at atomic resolution. K+ ions, but not Na+ ones, are critically required for the formation of this unique structure. The elucidated structure suggests the mechanisms underlying the biological events related to the GGA triplet repeat. Furthermore, in the light of the structure, the mode of the higher order packing of the telomere DNA is discussed.

Unique quadruplex structures of d(GGA)4 (12-mer) and d(GGA)8 (24-mer)--direct evidence of the formation of non-canonical base pairs and structural comparison.

We have reported that d(GGA)4 (12-mer) folds into an intramolecular parallel quadruplex with a G:G:G:G tetrad and a G(:A):G(:A):G(:A):G heptad and that two quadruplexes form a dimer. Here we present the unique structure of d(GGA)8 (24-mer) under physiological conditions. G1-A12 and G13-A24 segments of d(GGA)8 fold into an intramolecular parallel quadruplex with the tetrad and heptad, respectively, and two quadruplexes stack each other in a head-to-head manner. The formation of non-canonical G:G and G:A base pairs was directly certified by observing spin-spin couplings across hydrogen bonds with the aid of 13C-, 15N-labelling of DNA. It was noted that the monomeric architecture of d(GGA)8 resembles the dimeric architecture of d(GGA)4. Interestingly, the recent X-ray study shows that the telomere DNA also folds into the similar intramolecular parallel quadruplex as we have reported for GGA-repeat DNAs.

Quadruplex structures of RNA 14-mer, r(GGAGGUUUUGGAGG) and DNA 14-mer, d(GGAGGTTTTGGAGG).

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. We also synthesized the corresponding DNA 14-mer, d (GGAGGTTTTGGAGG) (D14), and examined its properties and structure by CD, gel electrophoresis, and NMR. 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 containing antiparalle strands and diagonal loops. The NMR study also gave the results that are consistent with such structure: alternate glycosidic conformation, 5'G(syn)-G(anti)3', and characteristic chemical shift data.

RNA and DNA, which contain two GGAGG segments connected with UUUU or TTTT sequences, form entirely different quadruplex structures.

We determined solution structure of d(GGAGGTTTTGGAGG) (D14) in the presence of Na+ ions by NMR. Two molecules of D14 form a quadruplex strucure with parallel and antiparallel strand alignments. The quadruplex is formed by four helical GGAGG segments and two diagonal TTTT loops at the top and bottom of the helix. This quadruplex structure is entirely different from that of r(GGAGGUUUUGGAGG) (R14), which has been previously determined by NMR. In the case of RNA, R14 forms an intra-strand parallel quadruplex with a UUUU loop a t lateral position ofthe helix and two such molecules form a dimer by stacking. The factors determining the type of RNA and DNA quadruplexes are discussed.

Structure of RNA aptamer for HIV Tat complexed with Tat-derived peptide.

An RNA aptamer containing two binding sites exhibits extremely high affinity to the HIV Tat protein. We previously reported the structure of the aptamer complexed with argininamide as the simplest analogue of Tat. Here, we have analyzed the structure of the aptamer complexed with the partial peptide of Tat, RKKRR. The profile of chemical sift perturbations for the aptamer upon complex formation with RKKRR revealed that RKKRR can be a realistic analogue of Tat to address the interactions between the arginine-rich motif of Tat and the aptamer. It was suggested that the aptamer interacts with different arginine residues of RKKRR simultaneously at the two binding sites, which can explain the extremely high affinity to Tat.