Noncanonical G(syn)-G(anti) base pairs stabilized by sulphate anions in two X-ray structures of the (GUGGUCUGAUGAGGCC) RNA duplex

The structures of two crystal forms of the RNA 16-mer with the sequence GUGGUCUGAUGAGGCC, grown in the presence of a high concentration of sulphate ions, have been determined using synchrotron radiation at 1.4- and 2.0-Å resolution. RNA with this sequence is known as one of the two strands of the noncleavable form of the hammerhead ribozyme. In both crystal structures, two G(syn)–G(anti) noncanonical base pairs are observed in the middle of a 14 base-pair (bp) duplex having 5′-dangling GU residues. Both structures contain sulphate anions interacting with the G–G bp stabilizing G in its syn conformation and bridging the two RNA strands. In both cases the interactions take place in the major groove, although the anions are accommodated within different helix geometries, most pronounced in the changing width of the major groove. In one structure, where a single sulphate spans both G–G pairs, the major groove is closed around the anion, while in the other structure, where each of the two G–G pairs is associated with a separate sulphate, the groove is open. This work provides the first examples of a G–G pair in syn-anti conformation, which minimizes the purine–purine clash in the center of the duplex, while utilizing its residual hydrogen bonding potential in specific interactions with sulphate anions.     

Database of non-canonical base pairs found in known RNA structures

Atomic resolution RNA structures are being published at an increasing rate. It is common to find a modest number of non-canonical base pairs in these structures in addition to the usual Watson-Crick pairs. This database summarizes the occurrence of these rare base pairs in accordance with standard nomenclature. The database, http://prion.bchs.uh.edu/, contains information such as sequence context, sugar pucker conformation, anti / syn base conformations, chemical shift, p K (a)values, melting temperature and free energy. Of the 29 anticipated pairs with two or more hydrogen bonds, 20 have been encountered to date. In addition, four unexpected pairs with two hydrogen bonds have been reported bringing the total to 24. Single hydrogen bond versions of five of the expected geometries have been encountered among the single hydrogen bond interactions. In addition, 18 different types of base triplets have been encountered, each of which involves three to six hydrogen bonds. The vast majority of the rare base pairs are antiparallel with the bases in the anti configuration relative to the ribose. The most common are the GU wobble, the Sheared GA pair, the Reverse Hoogsteen pair and the GA imino pair.     

ClaRNA: a classifier of contacts in RNA 3D structures based on a comparative analysis of various classification schemes

The understanding of folding and function of RNA molecules depends on the identification and classification of interactions between ribonucleotide residues. We developed a new method named ClaRNA for computational classification of contacts in RNA 3D structures. Unique features of the program are the ability to identify imperfect contacts and to process coarse-grained models. Each doublet of spatially close ribonucleotide residues in a query structure is compared to clusters of reference doublets obtained by analysis of a large number of experimentally determined RNA structures, and assigned a score that describes its similarity to one or more known types of contacts, including pairing, stacking, base–phosphate and base–ribose interactions. The accuracy of ClaRNA is 0.997 for canonical base pairs, 0.983 for non-canonical pairs and 0.961 for stacking interactions. The generalized squared correlation coefficient (GC²) for ClaRNA is 0.969 for canonical base pairs, 0.638 for non-canonical pairs and 0.824 for stacking interactions. The classifier can be easily extended to include new types of spatial relationships between pairs or larger assemblies of nucleotide residues. ClaRNA is freely available via a web server that includes an extensive set of tools for processing and visualizing structural information about RNA molecules.     

RNA structure and dynamics: A base pairing perspective

RNA is now known to possess various structural, regulatory and enzymatic functions for survival of cellular organisms. Functional RNA structures are generally created by three-dimensional organization of small structural motifs, formed by base pairing between self-complementary sequences from different parts of the RNA chain. In addition to the canonical Watson–Crick or wobble base pairs, several non-canonical base pairs are found to be crucial to the structural organization of RNA molecules. They appear within different structural motifs and are found to stabilize the molecule through long-range intra-molecular interactions between basic structural motifs like double helices and loops. These base pairs also impart functional variation to the minor groove of A-form RNA helices, thus forming anchoring site for metabolites and ligands. Non-canonical base pairs are formed by edge-to-edge hydrogen bonding interactions between the bases. A large number of theoretical studies have been done to detect and analyze these non-canonical base pairs within crystal or NMR derived structures of different functional RNA. Theoretical studies of these isolated base pairs using ab initio quantum chemical methods as well as molecular dynamics simulations of larger fragments have also established that many of these non-canonical base pairs are as stable as the canonical Watson–Crick base pairs. This review focuses on the various structural aspects of non-canonical base pairs in the organization of RNA molecules and the possible applications of these base pairs in predicting RNA structures with more accuracy.     

Inferring the conformation of RNA base pairs and triples from patterns of sequence variation.

The success of comparative analysis in resolving RNA secondary structure and numerous tertiary interactions relies on the presence of base covariations. Although the majority of base covariations in aligned sequences is associated to Watson-Crick base pairs, many involve non-canonical or restricted base pair exchanges (e.g. only G:C/A:U), reflecting more specific structural constraints. We have developed a computer program that determines potential base pairing conformations for a given set of paired nucleotides in a sequence alignment. This program (ISOPAIR) assumes that the base pair conformation is maintained through sequence variation without significantly affecting the path of the sugar-phosphate backbone. ISOPAIR identifies such 'isomorphic' structures for any set of input base pair or base triple sequences. The program was applied to base pairs and triples with known structures and sequence exchanges. In several instances, isomorphic structures were correctly identified with ISOPAIR. Thus, ISOPAIR is useful when assessing non-canonical base pair conformations in comparative analysis. ISOPAIR applications are limited to those cases where unusual base pair exchanges indeed reflect a non-canonical conformation.     

Structural characterization of naturally occurring RNA single mismatches

RNA is known to be involved in several cellular processes; however, it is only active when it is folded into its correct 3D conformation. The folding, bending and twisting of an RNA molecule is dependent upon the multitude of canonical and non-canonical secondary structure motifs. These motifs contribute to the structural complexity of RNA but also serve important integral biological functions, such as serving as recognition and binding sites for other biomolecules or small ligands. One of the most prevalent types of RNA secondary structure motifs are single mismatches, which occur when two canonical pairs are separated by a single non-canonical pair. To determine sequence–structure relationships and to identify structural patterns, we have systematically located, annotated and compared all available occurrences of the 30 most frequently occurring single mismatch-nearest neighbor sequence combinations found in experimentally determined 3D structures of RNA-containing molecules deposited into the Protein Data Bank. Hydrogen bonding, stacking and interaction of nucleotide edges for the mismatched and nearest neighbor base pairs are described and compared, allowing for the identification of several structural patterns. Such a database and comparison will allow researchers to gain insight into the structural features of unstudied sequences and to quickly look-up studied sequences.     

Loss of G-A base pairs is insufficient for achieving a large opening of U4 snRNA K-turn motif

Upon binding to the 15.5K protein, two tandem-sheared G–A base pairs are formed in the internal loop of the kink-turn motif of U4 snRNA (Kt-U4). We have reported that the folding of Kt-U4 is assisted by protein binding. Unstable interactions that contribute to a large opening of the free RNA (‘k–e motion’) were identified using locally enhanced sampling molecular dynamics simulations, results that agree with experiments. A detailed analysis of the simulations reveals that the k–e motion in Kt-U4 is triggered both by loss of G–A base pairs in the internal loop and backbone flexibility in the stems. Essential dynamics show that the loss of G–A base pairs is correlated along the first mode but anti-correlated along the third mode with the k–e motion. Moreover, when enhanced sampling was confined to the internal loop, the RNA adopted an alternative conformation characterized by a sharper kink, opening of G–A base pairs and modified stacking interactions. Thus, loss of G–A base pairs is insufficient for achieving a large opening of the free RNA. These findings, supported by previously published RNA structure probing experiments, suggest that G–A base pair formation occurs upon protein binding, thereby stabilizing a selective orientation of the stems.     

Non-canonical base pairs and higher order structures in nucleic acids: crystal structure database analysis

Non-canonical base pairs, mostly present in the RNA, often play a prominent role towards maintaining their structural diversity. Higher order structures like base triples are also important in defining and stabilizing the tertiary folded structure of RNA. We have developed a new program BPFIND to analyze different types of canonical and non-canonical base pairs and base triples involving at least two direct hydrogen bonds formed between polar atoms of the bases or sugar O2' only. We considered 104 possible types of base pairs, out of which examples of 87 base pair types are found to occur in the available RNA crystal structures. Analysis indicates that approximately 32.7% base pairs in the functional RNA structures are non-canonical, which include different types of GA and GU Wobble base pairs apart from a wide range of base pair possibilities. We further noticed that more than 10.4% of these base pairs are involved in triplet formation, most of which play important role in maintaining long-range tertiary contacts in the three-dimensional folded structure of RNA. Apart from detection, the program also gives a quantitative estimate of the conformational deformation of detected base pairs in comparison to an ideal planar base pair. This helps us to gain insight into the extent of their structural variations and thus assists in understanding their specific role towards structural and functional diversity.     

The first crystal structures of RNA–PNA duplexes and a PNA-PNA duplex containing mismatches—toward anti-sense therapy against TREDs

PNA is a promising molecule for antisense therapy of trinucleotide repeat disorders. We present the first crystal structures of RNA–PNA duplexes. They contain CUG repeats, relevant to myotonic dystrophy type I, and CAG repeats associated with poly-glutamine diseases. We also report the first PNA–PNA duplex containing mismatches. A comparison of the PNA homoduplex and the PNA–RNA heteroduplexes reveals PNA''s intrinsic structural properties, shedding light on its reported sequence selectivity or intolerance of mismatches when it interacts with nucleic acids. PNA has a much lower helical twist than RNA and the resulting duplex has an intermediate conformation. PNA retains its overall conformation while locally there is much disorder, especially peptide bond flipping. In addition to the Watson–Crick pairing, the structures contain interesting interactions between the RNA''s phosphate groups and the Π electrons of the peptide bonds in PNA.     

Insights into the DNA stabilizing contributions of a bicyclic cytosine analogue: crystal structures of DNA duplexes containing 7,8-dihydropyrido [2,3-d]pyrimidin-2-one

The incorporation of the bicyclic cytosine analogue 7,8-dihydropyrido[2,3-d]pyrimidin-2-one (X) into DNA duplexes results in a significant enhancement of their stability (3–4 K per modification). To establish the effects of X on the local hydrogen-bonding and base stacking interactions and the overall DNA conformation, and to obtain insights into the correlation between the structure and stability of X-containing DNA duplexes, the crystal structures of [d(CGCGAATT-X-GCG)]₂ and [d(CGCGAAT-X-CGCG)]₂ have been determined at 1.9–2.9 Å resolutions. In all of the structures, the analogue X base pairs with the purine bases on the opposite strands through Watson–Crick and/or wobble type hydrogen bonds. The additional ring of the X base is stacked on the thymine bases at the 5′-side and overall exhibits greatly enhanced stacking interactions suggesting that this is a major contribution to duplex stabilization.     

Helical coherence of DNA in crystals and solution

The twist, rise, slide, shift, tilt and roll between adjoining base pairs in DNA depend on the identity of the bases. The resulting dependence of the double helix conformation on the nucleotide sequence is important for DNA recognition by proteins, packaging and maintenance of genetic material, and other interactions involving DNA. This dependence, however, is obscured by poorly understood variations in the stacking geometry of the same adjoining base pairs within different sequence contexts. In this article, we approach the problem of sequence-dependent DNA conformation by statistical analysis of X-ray and NMR structures of DNA oligomers. We evaluate the corresponding helical coherence length—a cumulative parameter quantifying sequence-dependent deviations from the ideal double helix geometry. We find, e.g. that the solution structure of synthetic oligomers is characterized by 100–200 Å coherence length, which is similar to ~150 Å coherence length of natural, salmon-sperm DNA. Packing of oligomers in crystals dramatically alters their helical coherence. The coherence length increases to 800–1200 Å, consistent with its theoretically predicted role in interactions between DNA at close separations.     

Ion-induced folding of a kink turn that departs from the conventional sequence

Kink turns (k-turns) are important structural motifs that create a sharp axial bend in RNA. Most conform to a consensus in which a three-nucleotide bulge is followed by consecutive G•A and A•G base pairs, and when these G•A pairs are modified in vitro this generally leads to a failure to adopt the k-turn conformation. Kt-23 in the 30S ribosomal subunit of Thermus thermophilus is a rare exception in which the bulge-distal A•G pair is replaced by a non-Watson–Crick A•U pair. In the context of the ribosome, Kt-23 adopts a completely conventional k-turn geometry. We show here that this sequence is induced to fold into a k-turn structure in an isolated RNA duplex by Mg²⁺ or Na⁺ ions. Therefore, the Kt-23 is intrinsically stable despite lacking the key A•G pair; its formation requires neither tertiary interactions nor protein binding. Moreover, the Kt-23 k-turn is stabilized by the same critical hydrogen-bonding interactions within the core of the structure that are found in more conventional sequences such as the near-consensus Kt-7. T. thermophilus Kt-23 has two further non-Watson–Crick base pairs within the non-canonical helix, three and four nucleotides from the bulge, and we find that the nature of these pairs influences the ability of the RNA to adopt k-turn conformation, although the base pair adjacent to the A•U pair is more important than the other.     

NMR structure of the DNA decamer duplex containing double T*G mismatches of cis-syn cyclobutane pyrimidine dimer: implications for DNA damage recognition by the XPC-hHR23B complex

The cis–syn cyclobutane pyrimidine dimer (CPD) is a cytotoxic, mutagenic and carcinogenic DNA photoproduct and is repaired by the nucleotide excision repair (NER) pathway in mammalian cells. The XPC–hHR23B complex as the initiator of global genomic NER binds to sites of certain kinds of DNA damage. Although CPDs are rarely recognized by the XPC–hHR23B complex, the presence of mismatched bases opposite a CPD significantly increased the binding affinity of the XPC–hHR23B complex to the CPD. In order to decipher the properties of the DNA structures that determine the binding affinity for XPC–hHR23B to DNA, we carried out structural analyses of the various types of CPDs by NMR spectroscopy. The DNA duplex which contains a single 3′ T·G wobble pair in a CPD (CPD/GA duplex) induces little conformational distortion. However, severe distortion of the helical conformation occurs when a CPD contains double T·G wobble pairs (CPD/GG duplex) even though the T residues of the CPD form stable hydrogen bonds with the opposite G residues. The helical bending angle of the CPD/GG duplex was larger than those of the CPD/GA duplex and properly matched CPD/AA duplex. The fluctuation of the backbone conformation and significant changes in the widths of the major and minor grooves at the double T·G wobble paired site were also observed in the CPD/GG duplex. These structural features were also found in a duplex that contains the (6–4) adduct, which is efficiently recognized by the XPC–hHR23B complex. Thus, we suggest that the unique structural features of the DNA double helix (that is, helical bending, flexible backbone conformation, and significant changes of the major and/or minor grooves) might be important factors in determining the binding affinity of the XPC–hHR23B complex to DNA.     

Analysis of RNA flexibility by scanning force spectroscopy

Scanning force spectroscopy was used to measure the mechanical properties of double stranded RNA molecules in comparison with DNA. We find that, similar to the B–S transition in DNA, RNA molecules are stretched from the assumed A′ conformation to a stretched conformation by applying a defined force (plateau force). The force depends on the G + C content of the RNA and is distinct from that required for the B–S transition of a homologous DNA molecule. After the conformational change, DNA can be further extended by a factor of 0.7 ± 0.2 (S-factor) before melting occurs and the binding of the molecule to the cantilever is finally disrupted. For RNA, the S-factor was higher (1.0 ± 0.2) and more variable. Experiments to measure secondary structures in single stranded RNA yielded a large number of different force-distance curves, suggesting disruption and stretching of various secondary structures. Oriented attachment of the molecules to the substrate, a defined pick-up point and an increased resolution of the instrument could provide the means to analyse RNA secondary structures by scanning force spectroscopy.     

Structural characterization of cationic lipid–tRNA complexes

Despite considerable interest and investigations on cationic lipid–DNA complexes, reports on lipid–RNA interaction are very limited. In contrast to lipid–DNA complexes where lipid binding induces partial B to A and B to C conformational changes, lipid–tRNA complexation preserves tRNA folded state. This study is the first attempt to investigate the binding of cationic lipid with transfer RNA and the effect of lipid complexation on tRNA aggregation and condensation. We examine the interaction of tRNA with cholesterol (Chol), 1,2-dioleoyl-3-trimethylammonium-propane (DOTAP), dioctadecyldimethylammoniumbromide (DDAB) and dioleoylphosphatidylethanolamine (DOPE), at physiological condition, using constant tRNA concentration and various lipid contents. FTIR, UV-visible, CD spectroscopic methods and atomic force microscopy (AFM) were used to analyze lipid binding site, the binding constant and the effects of lipid interaction on tRNA stability, conformation and condensation. Structural analysis showed lipid–tRNA interactions with G–C and A–U base pairs as well as the backbone phosphate group with overall binding constants of K_Chol = 5.94 (± 0.8) × 10⁴ M^–1, K_DDAB = 8.33 (± 0.90) × 10⁵ M^–1, K_DOTAP = 1.05 (± 0.30) × 10⁵ M^–1 and K_DOPE = 2.75 (± 0.50) × 10⁴ M^–1. The order of stability of lipid–tRNA complexation is DDAB > DOTAP > Chol > DOPE. Hydrophobic interactions between lipid aliphatic tails and tRNA were observed. RNA remains in A-family structure, while biopolymer aggregation and condensation occurred at high lipid concentrations.     

Determination of thermodynamic parameters for HIV DIS type loop-loop kissing complexes

The HIV-1 type dimerization initiation signal (DIS) loop was used as a starting point for the analysis of the stability of Watson–Crick (WC) base pairs in a tertiary structure context. We used ultraviolet melting to determine thermodynamic parameters for loop–loop tertiary interactions and compared them with regular secondary structure RNA helices of the same sequences. In 1 M Na⁺ the loop–loop interaction of a HIV-1 DIS type pairing is 4 kcal/mol more stable than its sequence in an equivalent regular and isolated RNA helix. This difference is constant and sequence independent, suggesting that the rules governing the stability of WC base pairs in the secondary structure context are also valid for WC base pairs in the tertiary structure context. Moreover, the effect of ion concentration on the stability of loop–loop tertiary interactions differs considerably from that of regular RNA helices. The stabilization by Na⁺ and Mg²⁺ is significantly greater if the base pairing occurs within the context of a loop–loop interaction. The dependence of the structural stability on salt concentration was defined via the slope of a T_m/log [ion] plot. The short base-paired helices are stabilized by 8°C/log [Mg²⁺] or 11°C/log [Na⁺], whereas base-paired helices forming tertiary loop–loop interactions are stabilized by 16°C/log [Mg²⁺] and 26°C/log [Na⁺]. The different dependence on ionic strength that is observed might reflect the contribution of specific divalent ion binding to the preformation of the hairpin loops poised for the tertiary kissing loop–loop contacts.     

Visualization and characterization of RNA–protein interactions in living cells

RNA–protein interactions are the structural and functional basis of significant numbers of RNA molecules. RNA–protein interaction assays though, still mainly depend on biochemical tests in vitro. Here, we establish a convenient and reliable RNA fluorescent three-hybrid (rF3H) method to detect/interrogate the interactions between RNAs and proteins in cells. A GFP tagged highly specific RNA trap is constructed to anchor the RNA of interest to an artificial or natural subcellular structure, and RNA–protein interactions can be detected and visualized by the enrichment of RNA binding proteins (RBPs) at these structures. Different RNA trapping systems are developed and detection of RNA–protein complexes at multiple subcellular structures are assayed. With this new toolset, interactions between proteins and mRNA or noncoding RNAs are characterized, including the interaction between a long noncoding RNA and an epigenetic modulator. Our approach provides a flexible and reliable method for the characterization of RNA–protein interactions in living cells.     

Revealing A-T and G-C Hoogsteen base pairs in stressed protein-bound duplex DNA

Watson–Crick base pairs (bps) are the fundamental unit of genetic information and the building blocks of the DNA double helix. However, A-T and G-C can also form alternative ‘Hoogsteen’ bps, expanding the functional complexity of DNA. We developed ‘Hoog-finder’, which uses structural fingerprints to rapidly screen Hoogsteen bps, which may have been mismodeled as Watson–Crick in crystal structures of protein–DNA complexes. We uncovered 17 Hoogsteen bps, 7 of which were in complex with 6 proteins never before shown to bind Hoogsteen bps. The Hoogsteen bps occur near mismatches, nicks and lesions and some appear to participate in recognition and damage repair. Our results suggest a potentially broad role for Hoogsteen bps in stressed regions of the genome and call for a community-wide effort to identify these bps in current and future crystal structures of DNA and its complexes.     

Crystal Structure of a Four-Layer Aggregate of Engineered TMV CP Implies the Importance of Terminal Residues for Oligomer Assembly

Background

Crystal structures of the tobacco mosaic virus (TMV) coat protein (CP) in its helical and disk conformations have previously been determined at the atomic level. For the helical structure, interactions of proteins and nucleic acids in the main chains were clearly observed; however, the conformation of residues at the C-terminus was flexible and disordered. For the four-layer aggregate disk structure, interactions of the main chain residues could only be observed through water–mediated hydrogen bonding with protein residues. In this study, the effects of the C-terminal peptides on the interactions of TMV CP were investigated by crystal structure determination.

Methodology/Principal Findings

The crystal structure of a genetically engineered TMV CP was resolved at 3.06 Å. For the genetically engineered TMV CP, a six-histidine (His) tag was introduced at the N-terminus, and the C-terminal residues 155 to 158 were truncated (N-His-TMV CP¹⁹). Overall, N-His-TMV CP¹⁹ protein self-assembled into the four-layer aggregate form. The conformations of residues Gln36, Thr59, Asp115 and Arg134 were carefully analyzed in the high radius and low radius regions of N-His-TMV CP¹⁹, which were found to be significantly different from those observed previously for the helical and four-layer aggregate forms. In addition, the aggregation of the N-His-TMV CP¹⁹ layers was found to primarily be mediated through direct hydrogen-bonding. Notably, this engineered protein also can package RNA effectively and assemble into an infectious virus particle.

Conclusion

The terminal sequence of amino acids influences the conformation and interactions of the four-layer aggregate. Direct protein–protein interactions are observed in the major overlap region when residues Gly155 to Thr158 at the C-terminus are truncated. This engineered TMV CP is reassembled by direct protein–protein interaction and maintains the normal function of the four-layer aggregate of TMV CP in the presence of RNA.