RNA helical packing in solution: NMR structure of a 30 kDa GAAA tetraloop-receptor complex

Tertiary interactions are critical for proper RNA folding and ribozyme catalysis. RNA tertiary structure is often condensed through long-range helical packing interactions mediated by loop-receptor motifs. RNA structures displaying helical packing by loop-receptor interactions have been solved by X-ray crystallography, but not by NMR. Here, we report the NMR structure of a 30 kDa GAAA tetraloop-receptor RNA complex. In order to stabilize the complex, we used a modular design in which the RNA was engineered to form a homodimer, with each subunit containing a GAAA tetraloop phased one helical turn apart from its cognate 11-nucleotide receptor domain. The structure determination utilized specific isotopic labeling patterns (2H, 13C and 15N) and refinement against residual dipolar couplings. We observe a unique and highly unusual chemical shift pattern for an adenosine platform interaction that reveals a spectroscopic fingerprint for this motif. The structure of the GAAA tetraloop-receptor interaction is well defined solely from experimental NMR data, shows minor deviations from previously solved crystal structures, and verifies the previously inferred hydrogen bonding patterns within this motif. This work demonstrates the feasibility of using engineered homodimers as modular systems for the determination of RNA tertiary interactions by NMR.     

Crystal structure of an RNA aptamer bound to thrombin

Aptamers, an emerging class of therapeutics, are DNA or RNA molecules that are selected to bind molecular targets that range from small organic compounds to large proteins. All of the determined structures of aptamers in complex with small molecule targets show that aptamers cage such ligands. In structures of aptamers in complex with proteins that naturally bind nucleic acid, the aptamers occupy the nucleic acid binding site and often mimic the natural interactions. Here we present a crystal structure of an RNA aptamer bound to human thrombin, a protein that does not naturally bind nucleic acid, at 1.9 A resolution. The aptamer, which adheres to thrombin at the binding site for heparin, presents an extended molecular surface that is complementary to the protein. Protein recognition involves the stacking of single-stranded adenine bases at the core of the tertiary fold with arginine side chains. These results exemplify how RNA aptamers can fold into intricate conformations that allow them to interact closely with extended surfaces on non-RNA binding proteins.     

The adenosine wedge: A new structural motif in ribosomal RNA

Here, we present a new recurrent RNA arrangement, the so-called adenosine wedge (A-wedge), which is found in three places of the ribosomal RNA in both ribosomal subunits. The arrangement has a hierarchical structure, consisting of elements previously described as recurrent motifs, namely, the along-groove packing motif, the A-minor and the hook-turn. Within the A-wedge, these elements are involved in different types of cause–effect relationships, providing together for the particular tertiary structure of the motif.     

Determining the minimum number of types necessary to represent the sizes of protein atoms.

MOTIVATION: Traditionally, for packing calculations people have collected atoms together into a number of distinct 'types'. These, in fact, often represent a heavy atom and its associated hydrogens (i.e. a united atom). Also, atom typing is usually done according to basic chemistry, giving rise to 20-30 protein atom types, such as carbonyl carbons, methyl groups, and hydroxyl groups. No one has yet investigated how similar in packing these chemically derived types are. Here we address this question in detail, using Voronoi volume calculations on a set of high-resolution crystal structures. RESULTS: We perform a rigorous clustering analysis with cross-validation on tens of thousands of atom volumes and attempt to compile them into types based purely on packing. From our analysis, we are able to determine a 'minimal' set of 18 atom types that most efficiently represent the spectrum of packing in proteins. Furthermore, we are able to uncover a number of inconsistencies in traditional chemical typing schemes, where differently typed atoms have almost the same effective size. In particular, we find that tetrahedral carbons with two hydrogens are almost identical in size to many aromatic carbons with a single hydrogen. AVAILABILITY: Programs available from http://geometry.molmovdb.org. CONTACT: JerryTsai@TAMU.edu; neil.voss@yale.edu; Mark.Gerstein@yale.edu SUPPLEMENTARY INFORMATION: Available at http://geometry.molmovdb.org.     

Annotation of tertiary interactions in RNA structures reveals variations and correlations

RNA tertiary motifs play an important role in RNA folding and biochemical functions. To help interpret the complex organization of RNA tertiary interactions, we comprehensively analyze a data set of 54 high-resolution RNA crystal structures for motif occurrence and correlations. Specifically, we search seven recognized categories of RNA tertiary motifs (coaxial helix, A-minor, ribose zipper, pseudoknot, kissing hairpin, tRNA D-loop/T-loop, and tetraloop-tetraloop receptor) by various computer programs. For the nonredundant RNA data set, we find 613 RNA tertiary interactions, most of which occur in the 16S and 23S rRNAs. An analysis of these motifs reveals the diversity and variety of A-minor motif interactions and the various possible loop-loop receptor interactions that expand upon the tetraloop-tetraloop receptor. Correlations between motifs, such as pseudoknot or coaxial helix with A-minor, reveal higher-order patterns. These findings may ultimately help define tertiary structure restraints for RNA tertiary structure prediction. A complete annotation of the RNA diagrams for our data set is available at http://www.biomath.nyu.edu/motifs/.     

Nanometer scale pores similar in size to the entrance of the ribosomal exit cavity are a common feature of large RNAs

The highly conserved peptidyl transferase center (PTC) of the ribosome contains an RNA pore that serves as the entrance to the exit tunnel. Analysis of available ribosome crystal structures has revealed the presence of multiple additional well-defined pores of comparable size in the ribosomal (rRNA) RNAs. These typically have dimensions of 1–2 nm, with a total area of ∼100 Ų or more, and most are associated with one or more ribosomal proteins. The PTC example and the other rRNA pores result from the packing of helices. However, in the non-PTC cases the nitrogenous bases do not protrude into the pore, thereby limiting the potential for hydrogen bonding within the pore. Instead, it is the RNA backbone that largely defines the pore likely resulting in a negatively charged environment. In many but not all cases, ribosomal proteins are associated with the pores to a greater or lesser extent. With the exception of the PTC case, the large subunit pores are not found in what are thought to be the evolutionarily oldest regions of the 23S rRNA. The unusual nature of the PTC pore may reflect a history of being created by hybridization between two or more RNAs early in evolution rather than simple folding of a single RNA. An initial survey of nonribosomal RNA crystal structures revealed additional pores, thereby showing that they are likely a general feature of RNA tertiary structure.     

Crystal structure of the 30 S ribosomal subunit from Thermus thermophilus: structure of the proteins and their interactions with 16 S RNA

We present a detailed analysis of the protein structures in the 30 S ribosomal subunit from Thermus thermophilus, and their interactions with 16 S RNA based on a crystal structure at 3.05 A resolution. With 20 different polypeptide chains, the 30 S subunit adds significantly to our data base of RNA structure and protein-RNA interactions. In addition to globular domains, many of the proteins have long, extended regions, either in the termini or in internal loops, which make extensive contact to the RNA component and are involved in stabilizing RNA tertiary structure. Many ribosomal proteins share similar alpha+beta sandwich folds, but we show that the topology of this domain varies considerably, as do the ways in which the proteins interact with RNA. Analysis of the protein-RNA interactions in the context of ribosomal assembly shows that the primary binders are globular proteins that bind at RNA multihelix junctions, whereas proteins with long extensions assemble later. We attempt to correlate the structure with a large body of biochemical and genetic data on the 30 S subunit.     

X‐ray vs. NMR structures as templates for computational protein design

Certain protein‐design calculations involve using an experimentally determined high‐resolution structure as a template to identify new sequences that can adopt the same fold. This approach has led to the successful design of many novel, well‐folded, native‐like proteins. Although any atomic‐resolution structure can serve as a template in such calculations, most successful designs have used high‐resolution crystal structures. Because there are many proteins for which crystal structures are not available, it is of interest whether nuclear magnetic resonance (NMR) templates are also appropriate. We have analyzed differences between using X‐ray and NMR templates in side‐chain repacking and design calculations. We assembled a database of 29 proteins for which both a high‐resolution X‐ray structure and an ensemble of NMR structures are available. Using these pairs, we compared the rotamericity, χ₁‐angle recovery, and native‐sequence recovery of X‐ray and NMR templates. We carried out design using RosettaDesign on both types of templates, and compared the energies and packing qualities of the resulting structures. Overall, the X‐ray structures were better templates for use with Rosetta. However, for ～20% of proteins, a member of the reported NMR ensemble gave rise to designs with similar properties. Re‐evaluating RosettaDesign structures with other energy functions indicated much smaller differences between the two types of templates. Ultimately, experiments are required to confirm the utility of particular X‐ray and NMR templates. But our data suggest that the lack of a high‐resolution X‐ray structure should not preclude attempts at computational design if an NMR ensemble is available. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

Sequence and structural conservation in RNA ribose zippers

The "ribose zipper", an important element of RNA tertiary structure, is characterized by consecutive hydrogen-bonding interactions between ribose 2'-hydroxyls from different regions of an RNA chain or between RNA chains. These tertiary contacts have previously been observed to also involve base-backbone and base-base interactions (A-minor type). We searched for ribose zipper tertiary interactions in the crystal structures of the large ribosomal subunit RNAs of Haloarcula marismortui and Deinococcus radiodurans, and the small ribosomal subunit RNA of Thermus thermophilus and identified a total of 97 ribose zippers. Of these, 20 were found in T. thermophilus 16 S rRNA, 44 in H. marismortui 23 S rRNA (plus 2 bridging 5 S and 23 S rRNAs) and 30 in D. radiodurans 23 S rRNA (plus 1 bridging 5 S and 23 S rRNAs). These were analyzed in terms of sequence conservation, structural conservation and stability, location in secondary structure, and phylogenetic conservation. Eleven types of ribose zippers were defined based on ribose-base interactions. Of these 11, seven were observed in the ribosomal RNAs. The most common of these is the canonical ribose zipper, originally observed in the P4-P6 group I intron fragment. All ribose zippers were formed by antiparallel chain interactions and only a single example extended beyond two residues, forming an overlapping ribose zipper of three consecutive residues near the small subunit A-site. Almost all ribose zippers link stem (Watson-Crick duplex) or stem-like (base-paired), with loop (external, internal, or junction) chain segments. About two-thirds of the observed ribose zippers interact with ribosomal proteins. Most of these ribosomal proteins bridge the ribose zipper chain segments with basic amino acid residues hydrogen bonding to the RNA backbone. Proteins involved in crucial ribosome function and in early stages of ribosomal assembly also stabilize ribose zipper interactions. All ribose zippers show strong sequence conservation both within these three ribosomal RNA structures and in a large database of aligned prokaryotic sequences. The physical basis of the sequence conservation is stacked base triples formed between consecutive base-pairs on the stem or stem-like segment with bases (often adenines) from the loop-side segment. These triples have previously been characterized as Type I and Type II A-minor motifs and are stabilized by base-base and base-ribose hydrogen bonds. The sequence and structure conservation of ribose zippers can be directly used in tertiary structure prediction and may have applications in molecular modeling and design.     

Evidence for ligandable sites in structured RNA throughout the Protein Data Bank

RNA has attracted considerable attention as a target for small molecules. However, methods to identify, study, and characterize suitable RNA targets have lagged behind strategies for protein targets. One approach that has received considerable attention for protein targets has been to utilize computational analysis to investigate ligandable “pockets” on proteins that are amenable to small molecule binding. These studies have shown that selected physical properties of pockets are important parameters that govern the ability of a structure to bind to small molecules. This work describes a similar analysis to study pockets on all RNAs in the Protein Data Bank (PDB). Using parameters such as buriedness, hydrophobicity, volume, and other properties, the set of all RNAs is analyzed and compared to all proteins. Considerable overlap is observed between the properties of pockets on RNAs and proteins. Thus, many RNAs are capable of populating conformations with pockets that are likely suitable for small molecule binding. Further, principal moment of inertia (PMI) calculations reveal that liganded RNAs exist in diverse structural space, much of which overlaps with protein structural space. Taken together, these results suggest that complex folded RNAs adopt unique structures with pockets that may represent viable opportunities for small molecule targeting.     

Diversity of base-pair conformations and their occurrence in rRNA structure and RNA structural motifs

In addition to the canonical base-pairs comprising the standard Watson-Crick (C:G and U:A) and wobble U:G conformations, an analysis of the base-pair types and conformations in the rRNAs in the high-resolution crystal structures of the Thermus thermophilus 30S and Haloarcula marismortui 50S ribosomal subunits has identified a wide variety of non-canonical base-pair types and conformations. However, the existing nomenclatures do not describe all of the observed non-canonical conformations or describe them with some ambiguity. Thus, a standardized system is required to classify all of these non-canonical conformations appropriately. Here, we propose a new, simple and systematic nomenclature that unambiguously classifies base-pair conformations occurring in base-pairs, base-triples and base-quadruples that are associated with secondary and tertiary interactions. This system is based on the topological arrangement of the two bases and glycosidic bonds in a given base-pair. Base-pairs in the internal positions of regular secondary structure helices usually form with canonical base-pair groups (C:G, U:A, and U:G) and canonical conformations (C:G WC, U:A WC, and U:G Wb). In contrast, non-helical base-pairs outside of regular structure helices usually have non-canonical base-pair groups and conformations. In addition, many non-helical base-pairs are involved in RNA motifs that form a defined set of non-canonical conformations. Thus, each rare non-canonical conformation may be functionally and structurally important. Finally, the topology-based isostericity of base-pair conformations can rationalize base-pair exchanges in the evolution of RNA molecules.     

RNA helicases: modulators of RNA structure

RNA molecules play an essential role in many cellular processes, often as components of ribonucleoprotein complexes. Like proteins, RNA molecules adopt sequence-specific secondary and tertiary structures that are essential for function; alteration of these structures therefore provides a means of regulating RNA function. The discovery of DEAD box proteins, a large family of proteins that share several highly conserved motifs and have known or putative ATP-dependent RNA helicase activity, has provoked growing interest in the concept that regulation of RNA function may occur through local unwinding of complex RNA structures.     

The polynucleotide ligase and RNA capping enzyme superfamily of covalent nucleotidyltransferases

ATP- and NAD(+)-dependent DNA ligases, ATP-dependent RNA ligases and GTP-dependent mRNA capping enzymes comprise a superfamily of proteins that catalyze nucleotidyl transfer to polynucleotide 5' ends via covalent enzyme-(lysyl-N)-NMP intermediates. The superfamily is defined by five peptide motifs that line the nucleotide-binding pocket and contribute amino acid sidechains essential for catalysis. Early crystal structures revealed a shared core tertiary structure for DNA ligases and capping enzymes, which are composed minimally of a nucleotidyltransferase domain fused to a distal OB-fold domain. Recent structures of viral and bacterial DNA ligases, and a fungal mRNA capping enzyme underscore how the substrate-binding and chemical steps of the ligation and capping pathways are coordinated with large rearrangements of the component protein domains and with remodeling of the atomic contacts between the enzyme and the nucleotide at the active site. The first crystal structure of an RNA ligase suggests that contemporary DNA ligases, RNA ligases and RNA capping enzymes evolved by fusion of ancillary effector domains to an ancestral catalytic module involved in RNA repair.     

Calculations of protein volumes: sensitivity analysis and parameter database

MOTIVATION: The precise sizes of protein atoms in terms of occupied packing volume are of great importance. We have previously presented standard volumes for protein residues based on calculations with Voronoi-like polyhedra. To understand the applicability and limitations of our set, we investigated, in detail, the sensitivity of the volume calculations to a number of factors: (i) the van der Waals radii set, (ii) the criteria for including buried atoms in the calculations or atom selection, (iii) the method of positioning the dividing plane in polyhedra construction, and (iv) the set of structures used in the averaging. RESULTS: We find that different radii sets have only moderate affects to the distribution and mean of volumes. Atom selection and dividing plane methods cause larger changes in protein atoms volumes. More significantly, we show how the variation in volumes appears to be clearly related to the quality of the structures analyzed, with higher quality structures giving consistently smaller average volumes with less variance.     

RNA as a small molecule druggable target

Small molecule drugs have readily been developed against many proteins in the human proteome, but RNA has remained an elusive target for drug discovery. Increasingly, we see that RNA, and to a lesser extent DNA elements, show a persistent tertiary structure responsible for many diverse and complex cellular functions. In this digest, we have summarized recent advances in screening approaches for RNA targets and outlined the discovery of novel, drug-like small molecules against RNA targets from various classes and therapeutic areas. The link of structure, function, and small-molecule Druggability validates now for the first time that RNA can be the targets of therapeutic agents.     

Cavities in protein–DNA and protein–RNA interfaces

An analysis of cavities present in protein–DNA and protein–RNA complexes is presented. In terms of the number of cavities and their total volume, the interfaces formed in these complexes are akin to those in transient protein–protein heterocomplexes. With homodimeric proteins protein–DNA interfaces may contain cavities involving both the protein subunits and DNA, and these are more than twice as large as cavities involving a single protein subunit and DNA. A parameter, cavity index, measuring the degree of surface complementarity, indicates that the packing of atoms in protein–protein/DNA/RNA is very similar, but it is about two times less efficient in the permanent interfaces formed between subunits in homodimers. As within the tertiary structure and protein–protein interfaces, protein–DNA interfaces have a higher inclination to be lined by β-sheet residues; from the DNA side, base atoms, in particular those in minor grooves, have a higher tendency to be located in cavities. The larger cavities tend to be less spherical and solvated. A small fraction of water molecules are found to mediate hydrogen-bond interactions with both the components, suggesting their primary role is to fill in the void left due to the local non-complementary nature of the surface patches.     

The packing density in proteins: standard radii and volumes.

The sizes of atomic groups are a fundamental aspect of protein structure. They are usually expressed in terms of standard sets of radii for atomic groups and of volumes for both these groups and whole residues. Atomic groups, which subsume a heavy-atom and its covalently attached hydrogen atoms into one moiety, are used because the positions of hydrogen atoms in protein structures are generally not known. We have calculated new values for the radii of atomic groups and for the volumes of atomic groups. These values should prove useful in the analysis of protein packing, protein recognition and ligand design. Our radii for atomic groups were derived from intermolecular distance calculations on a large number (approximately 30,000) of crystal structures of small organic compounds that contain the same atomic groups to those found in proteins. Our radii show significant differences to previously reported values. We also use this new radii set to determine the packing efficiency in different regions of the protein interior. This analysis shows that, if the surface water molecules are included in the calculations, the overall packing efficiency throughout the protein interior is high and fairly uniform. However, if the water structure is removed, the packing efficiency in peripheral regions of the protein interior is underestimated, by approximately 3.5 %.     

A general edit distance between RNA structures.

Arc-annotated sequences are useful in representing the structural information of RNA sequences. In general, RNA secondary and tertiary structures can be represented as a set of nested arcs and a set of crossing arcs, respectively. Since RNA functions are largely determined by molecular confirmation and therefore secondary and tertiary structures, the comparison between RNA secondary and tertiary structures has received much attention recently. In this paper, we propose the notion of edit distance to measure the similarity between two RNA secondary and tertiary structures, by incorporating various edit operations performed on both bases and arcs (i.e., base-pairs). Several algorithms are presented to compute the edit distance between two RNA sequences with various arc structures and under various score schemes, either exactly or approximately, with provably good performance. Preliminary experimental tests confirm that our definition of edit distance and the computation model are among the most reasonable ones ever studied in the literature.