Structural and thermodynamic studies on mutant RNA motifs that impair the specificity between a viral replicase and its promoter

The 3'-end region of the genomic RNA of brome mosaic virus forms a tRNA-like structure that is critical for its replication. Previous studies have shown that in this region, a stem-loop structure, called SLC, is necessary and sufficient for the binding of the RNA replicase, and for RNA replication. Recently, we determined the high-resolution NMR structure of SLC, which demonstrated that a 5'-AUA-3' triloop region is an important structural element for the enzymatic recognition. We proposed that the 5'-adenine of the triloop, which is rigidly fixed ("clamped") to the stem, is a key recognition element for the replicase. To elucidate the role of this "clamped base motif" for the enzymatic recognition, we have now investigated the solution conformations of several stem-loop molecules with mutant triloops, 5'-UUA-3', 5'-GUA-3', 5'-CUA-3' and 5'-UUU-3', that destroy the enzymatic recognition. For the GUA and UUA mutants, we have obtained high-resolution solution structures using 2D NMR. All four mutants have very similar thermodynamic stabilities, and all have the same secondary structures, a triloop with a five base-paired stem helix. In addition, they have quite similar sugar puckering patterns in the triloop region. The NMR structures of the GUA and UUA show that the 5' nucleotide of the triloop (G6 in GUA or U6 in UUA) lacks the strong interactions that hold its base in a fixed position. In particular, the U6 of UUA is found in two different conformations. Neither of these two mutants has the clamped base motif that was observed in the wild-type. While UUA also shows global change in the overall triloop conformation, GUA shows a very similar triloop conformation to the wild-type except for the lack of this motif. The absence of the clamped base motif is the only common structural difference between these two mutants and the wild-type. These results clearly indicate that the loss of function of the UUA and GUA mutants comes mainly from the destruction of a small key recognition motif rather than from global changes in their triloop conformations. Based on this study, we conclude that the key structural motif in the triloop recognized by the replicase is a solution-exposed, 5'-adenine base in the triloop that is clamped to the stem helix, which is called a clamped adenine motif.     

On the structural features of hairpin triloops in rRNA: from nucleotide to global conformational change upon ligand binding

RNA structure can be viewed as both a construct composed of various structural motifs and a flexible polymer that is substantially influenced by its environment. In this light, the present paper represents an attempt to reconcile the two standpoints. By using the 3D structures both of four (16S and 23S) portions of unbound 50S, H50S, and T30S ribosomal subunits and of 38 large ribonucleoligand complexes as the starting point, the behavior, which is induced by ligand binding, of 73 hairpin triloops with closing g-c and c-g base pairs was investigated using root-mean-square deviation (RMSD) approach and pseudotorsional (eta,theta) convention at the nucleotide-by-nucleotide level. Triloops were annotated in accordance with a recent proposal of geometric nomenclature. A simple measure for the determination of the strain of a triloop is introduced. It is believed that a possible classification of the interior triloops, based on the 2D eta-theta unique path, will aid to conceive their local behavior upon ligand binding. All rRNA residues in contact with ligands as well as regions of considerable conformational changes upon complex formation were identified. The analysis offers the answer to: how proximal to and how far from the actual ligand-binding sites the structural changes occur?     

The lonepair triloop: a new motif in RNA structure

The lonepair triloop (LPTL) is an RNA structural motif that contains a single ("lone") base-pair capped by a hairpin loop containing three nucleotides. The two nucleotides immediately outside of this motif (5' and 3' to the lonepair) are not base-paired to one another, restricting the length of this helix to a single base-pair. Four examples of this motif, along with three tentative examples, were initially identified in the 16S and 23S rRNAs with covariation analysis. An evaluation of the recently determined crystal structures of the Thermus thermophilus 30S and Haloarcula marismortui 50S ribosomal subunits revealed the authenticity for all of these proposed interactions and identified 16 more LPTLs in the 5S, 16S and 23S rRNAs. This motif is found in the T loop in the tRNA crystal structures. The lonepairs are positioned, in nearly all examples, immediately 3' to a regular secondary structure helix and are stabilized by coaxial stacking onto this flanking helix. In all but two cases, the nucleotides in the triloop are involved in a tertiary interaction with another section of the rRNA, establishing an overall three-dimensional function for this motif. Of these 24 examples, 14 occur in multi-stem loops, seven in hairpin loops and three in internal loops. While the most common lonepair, U:A, occurs in ten of the 24 LPTLs, the remaining 14 LPTLs contain seven different base-pair types. Only a few of these lonepairs adopt the standard Watson-Crick base-pair conformations, while the majority of the base-pairs have non-standard conformations. While the general three-dimensional conformation is similar for all examples of this motif, characteristic differences lead to several subtypes present in different structural environments. At least one triloop nucleotide in 22 of the 24 LPTLs in the rRNAs and tRNAs forms a tertiary interaction with another part of the RNA. When a LPTL containing the GNR or UYR triloop sequence forms a tertiary interaction with the first (and second) triloop nucleotide, it recruits a fourth nucleotide to mediate stacking and mimic the tetraloop conformation. Approximately half of the LPTL motifs are in close association with proteins. The majority of these LPTLs are positioned at sites in rRNAs that are conserved in the three phylogenetic domains; a few of these occur in regions of the rRNA associated with ribosomal function, including the presumed site of peptidyl transferase activity in the 23S rRNA.     

Sequence–structure relationships in RNA loops: establishing the basis for loop homology modeling

The specific function of RNA molecules frequently resides in their seemingly unstructured loop regions. We performed a systematic analysis of RNA loops extracted from experimentally determined three-dimensional structures of RNA molecules. A comprehensive loop-structure data set was created and organized into distinct clusters based on structural and sequence similarity. We detected clear evidence of the hallmark of homology present in the sequence–structure relationships in loops. Loops differing by <25% in sequence identity fold into very similar structures. Thus, our results support the application of homology modeling for RNA loop model building. We established a threshold that may guide the sequence divergence-based selection of template structures for RNA loop homology modeling. Of all possible sequences that are, under the assumption of isosteric relationships, theoretically compatible with actual sequences observed in RNA structures, only a small fraction is contained in the Rfam database of RNA sequences and classes implying that the actual RNA loop space may consist of a limited number of unique loop structures and conserved sequences. The loop-structure data sets are made available via an online database, RLooM. RLooM also offers functionalities for the modeling of RNA loop structures in support of RNA engineering and design efforts.     

Analyzing the sequence-structure relationship of a library of local structural prototypes

We present a thorough analysis of the relation between amino acid sequence and local three-dimensional structure in proteins. A library of overlapping local structural prototypes was built using an unsupervised clustering approach called “hybrid protein model” (HPM). The HPM carries out a multiple structural alignment of local folds from a non-redundant protein structure databank encoded into a structural alphabet composed of 16 protein blocks (PBs). Following previous research focusing on the HPM protocol, we have considered gaps in the local structure prototype. This methodology allows to have variable length fragments. Hence, 120 local structure prototypes were obtained. Twenty-five percent of the protein fragments learnt by HPM had gaps.An investigation of tight turns suggested that they are mainly derived from three PB series with precise locations in the HPM. The amino acid information content of the whole conformational classes was tackled by multivariate methods, e.g., canonical correlation analysis. It points out the presence of seven amino acid equivalence classes showing high propensities for preferential local structures. In the same way, definition of “contrast factors” based on sequence-structure properties underline the specificity of certain structural prototypes, e.g., the dependence of Gly or Asn-rich turns to a limited number of PBs, or, the opposition between Pro-rich coils to those enriched in Ser, Thr, Asn and Glu. These results are so useful to analyze the sequence-structure relationships, but could also be used to improve fragment-based method for protein structure prediction from sequence.     

A molecular model of the carbohydrate recognition domain of a rat macrophage lectin and analysis of its binding site

A three-dimensional model of the carbohydrate recognition domain of a rat macrophage C-type lectin has been constructed by comparative modeling and assessed by inverse folding analysis. Comparative modeling in the presence of low sequence similarity was based on information provided by comparison of X-ray structures and sequence-structure alignments. The sequence-structure compatibility of the model was sound. Its binding site was analyzed in comparison to the X-ray structure of a galactose-specific mutant of the mannose-binding protein. The specificity of the macrophage lectin was discussed in light of mutagenesis data on asialoglycoprotein receptors.     

Identification of local conformational similarity in structurally variable regions of homologous proteins using protein blocks

Structure comparison tools can be used to align related protein structures to identify structurally conserved and variable regions and to infer functional and evolutionary relationships. While the conserved regions often superimpose well, the variable regions appear non superimposable. Differences in homologous protein structures are thought to be due to evolutionary plasticity to accommodate diverged sequences during evolution. One of the kinds of differences between 3-D structures of homologous proteins is rigid body displacement. A glaring example is not well superimposed equivalent regions of homologous proteins corresponding to α-helical conformation with different spatial orientations. In a rigid body superimposition, these regions would appear variable although they may contain local similarity. Also, due to high spatial deviation in the variable region, one-to-one correspondence at the residue level cannot be determined accurately. Another kind of difference is conformational variability and the most common example is topologically equivalent loops of two homologues but with different conformations. In the current study, we present a refined view of the "structurally variable" regions which may contain local similarity obscured in global alignment of homologous protein structures. As structural alphabet is able to describe local structures of proteins precisely through Protein Blocks approach, conformational similarity has been identified in a substantial number of 'variable' regions in a large data set of protein structural alignments; optimal residue-residue equivalences could be achieved on the basis of Protein Blocks which led to improved local alignments. Also, through an example, we have demonstrated how the additional information on local backbone structures through protein blocks can aid in comparative modeling of a loop region. In addition, understanding on sequence-structure relationships can be enhanced through our approach. This has been illustrated through examples where the equivalent regions in homologous protein structures share sequence similarity to varied extent but do not preserve local structure.     

A structural bioinformatics approach to the analysis of nonsynonymous single nucleotide polymorphisms (nsSNPs) and their relation to disease

The prediction of the effects of nonsynonymous single nucleotide polymorphisms (nsSNPs) on function depends critically on exploiting all information available on the three-dimensional structures of proteins. We describe software and databases for the analysis of nsSNPs that allow a user to move from SNP to sequence to structure to function. In both structure prediction and the analysis of the effects of nsSNPs, we exploit information about protein evolution, in particular, that derived from investigations on the relation of sequence to structure gained from the study of amino acid substitutions in divergent evolution. The techniques developed in our laboratory have allowed fast and automated sequence-structure homology recognition to identify templates and to perform comparative modeling; as well as simple, robust, and generally applicable algorithms to assess the likely impact of amino acid substitutions on structure and interactions. We describe our strategy for approaching the relationship between SNPs and disease, and the results of benchmarking our approach -- human proteins of known structure and recognized mutation.     

Modeling RNA tertiary structure motifs by graph-grammars

A new approach, graph-grammars, to encode RNA tertiary structure patterns is introduced and exemplified with the classical sarcin-ricin motif. The sarcin-ricin motif is found in the stem of the crucial ribosomal loop E (also referred to as the sarcin-ricin loop), which is sensitive to the alpha-sarcin and ricin toxins. Here, we generate a graph-grammar for the sarcin-ricin motif and apply it to derive putative sequences that would fold in this motif. The biological relevance of the derived sequences is confirmed by a comparison with those found in known sarcin-ricin sites in an alignment of over 800 bacterial 23S ribosomal RNAs. The comparison raised alternative alignments in few sarcin-ricin sites, which were assessed using tertiary structure predictions and 3D modeling. The sarcin-ricin motif graph-grammar was built with indivisible nucleotide interaction cycles that were recently observed in structured RNAs. A comparison of the sequences and 3D structures of each cycle that constitute the sarcin-ricin motif gave us additional insights about RNA sequence-structure relationships. In particular, this analysis revealed the sequence space of an RNA motif depends on a structural context that goes beyond the single base pairing and base-stacking interactions.     

RNA motifs that determine specificity between a viral replicase and its promoter

The 3' end of brome mosaic virus RNA contains a tRNA-like sequence that directs its RNA synthesis. A stem loop structure in this sequence, stem loop C (SLC), was investigated using NMR, and correlated with its ability to direct RNA synthesis by its replicase. SLC consists of two discrete domains, a flexible stem with an internal loop and a rigid stem containing a 5'-AUA-3' triloop. Efficient RNA synthesis requires the sequence on only one side of the flexible stem and a specific compact conformation of the triloop. A high resolution structure of the triloop places the 5' adenine out in solution, and the 3' adenine within the triloop, held tightly through stacking and unusual hydrogen bonds. This high resolution structure of an RNA promoter from a (+)-strand RNA virus provides new insights into how the RNA-dependent RNA polymerase binds to the RNA to initiate synthesis.     

Solution structure of a truncated anti-MUC1 DNA aptamer determined by mesoscale modeling and NMR

Mucin 1 is a well-established target for the early diagnosis of epithelial cancers. The nucleotides of the S1.3/S2.2 DNA aptamer involved in binding to variable number tandem repeat mucin 1 peptides have been identified using footprinting experiments. The majority of these binding nucleotides are located in the 25-nucleotide variable region of the total aptamer. Imino proton and 2D NMR spectra of truncated and total aptamers in supercooled water reveal common hydrogen-bonding networks and point to a similar secondary structure for this 25-mer sequence alone or embedded within the total aptamer. NMR titration experiments confirm that the TTT triloop structure is the primary binding site and show that the initial structure of the truncated aptamers is conserved upon interaction with variable number tandem repeat peptides. The thermal dependence of the NMR chemical shift data shows that the base-paired nucleotides melt cooperatively at 47 ±?4°C. The structure of the 25-mer oligonucleotide was determined using a new combined mesoscale molecular modeling, molecular dynamics and NMR spectroscopy investigation. It contains three Watson-Crick pairs, three consecutive mispairs and four Watson-Crick pairs capped by a TTT triloop motif. The 3D model structures (PDB 2L5K) and biopolymer chain elasticity molecular models are consistent with both NMR and long unconstrained molecular dynamics (10 ns) in explicit water, respectively. Database Structural data are available in the Protein Data Bank and BioMagResBank databases under the accession numbers 2L5K and 17129, respectively.     

Practical tools for molecular modeling of complex carbohydrates and their interactions with proteins

Computer modeling has become a valuable component of studies of carbohydrate three-dimensional structures and their relationship to function and properties. In this paper we examine the methods required for conformational modeling of carbohydrates, and we present a series of tools that have been developed to this end. These tools can be integrated into three-dimensional real-time molecular modeling software. A data base of pre-optimized carbohydrate fragments has been established to be used further in the construction of much more complex molecules. In addition we describe some possible uses of a data base of three dimensional structures of the disaccharide fragments present in the glycan moiety ofN-glycoprotein. A molecular mechanical force field appropriate for the conformational analysis of oligosaccharides has been derived by the addition of new parameters to the Tripos force field and is compatible with protein simulations. The new parametrization has been assessed in three stages of increasing complexity: computations of potential energy surfaces, conformational refinement of relevant oligosaccharides, modeling at the atomic level of a protein/carbohydrate complex.     

RNA FRABASE 2.0: an advanced web-accessible database with the capacity to search the three-dimensional fragments within RNA structures

Background  

Recent discoveries concerning novel functions of RNA, such as RNA interference, have contributed towards the growing importance of the field. In this respect, a deeper knowledge of complex three-dimensional RNA structures is essential to understand their new biological functions. A number of bioinformatic tools have been proposed to explore two major structural databases (PDB, NDB) in order to analyze various aspects of RNA tertiary structures. One of these tools is RNA FRABASE 1.0, the first web-accessible database with an engine for automatic search of 3D fragments within PDB-derived RNA structures. This search is based upon the user-defined RNA secondary structure pattern. In this paper, we present and discuss RNA FRABASE 2.0. This second version of the system represents a major extension of this tool in terms of providing new data and a wide spectrum of novel functionalities. An intuitionally operated web server platform enables very fast user-tailored search of three-dimensional RNA fragments, their multi-parameter conformational analysis and visualization.     

Molecular Modeling of Immunoglobulin Superfamily Proteins: Predicting the Three Dimensional Structure of the Extracellular Domain of CTLA-4 (CD152)

The interactions between CD28/CTLA-4 (CD152) on T cells and their ligands CD80/CD86 on antigen presenting cells provide costimulatory signals critical for T cell activation. CD28/CTLA-4 and CD80/CD86 are members of the immunoglobulin superfamily (IgSF). CD28 and CTLA-4 both contain a single extracellular immunoglobulin (Ig) domain which binds CD80/CD86. Here we report modeling studies on the three-dimensional (3D) structure of the CTLA-4 binding domain. Since CTLA-4 displays only very weak sequence homology to proteins with known 3D structure, conventional modeling techniques were difficult to apply. Structure-oriented sequence comparison, consensus residue analysis, conformational searching, and inverse folding calculations were employed to aid in the generation of a comparative CTLA-4 model. Regions of high and low prediction confidence were identified, and the sequence-structure compatibility of the model was determined. Characteristics of the modeled structure, which resembles an Ig V domain, were analyzed, and the model was used to map N-linked glycosylation sites and residues critical for CTLA-4 function. The modeling approach described here can be applied to predict 3D structures of other IgSF proteins.     

Displaying 3D data on RNA secondary structures: coloRNA

RNA performs a variety of diverse functions and therefore must adopt many different three-dimensional conformations. The number and complexity of RNA structures that are currently available are steadily increasing, necessitating the generation of versatile structure visualization tools. Here, we describe a new RNA secondary and tertiary structure visualization tool, the display program coloRNA. This program colors each nucleotide in a secondary structure schematic according to the value of an assigned property of the corresponding backbone phosphate group, such as the distance between corresponding residues in two atomic models of the same RNA molecule. To assist in analyzing tertiary structure, coloRNA also colors nucleotides based on the three-dimensional distances between a user-selected nucleotide and all others. Minimum and maximum thresholds can be used to focus in on, or eliminate, a particular value range. coloRNA can display a user-specified group of nucleotides by outlining the structure in an automatically assigned, but user-changeable color. As an example, we have used coloRNA to analyze a pair of recently published structures of the Escherichia coli 70S ribosome. When coloRNA is used to display the conformational difference between the two structures, the large movement of the small subunit head stands visually out from the background changes in the remaining domains of the small subunit.     

Global mapping of nucleic acid conformational space: dinucleoside monophosphate conformations and transition pathways among conformational classes

A global conformational space of 6253 dinucleoside monophosphate (DMP) units consisting of RNA and DNA (free and protein/drug-bound) was 'mapped' using high resolution crystal structures cataloged in the Nucleic Acid Database (NDB). The torsion angles of each DMP were clustered in a reduced three-dimensional space using a classical multi-dimensional scaling method. The mapping of the conformational space reveals nine primary clusters which distinguish among the common A-, B- and Z-forms and their various substates, plus five secondary clusters for kinked or bent structures. Conformational relationships and possible transitional pathways among the substates are also examined using the conformational states of DNA and RNA bound with proteins or drugs as potential pathway intermediates.     

Programmed ribosomal frameshifting in SIV is induced by a highly structured RNA stem-loop

Simian immunodeficiency virus (SIV), like its human homologues (HIV-1, HIV-2), requires a -1 translational frameshift event to properly synthesize all of the proteins required for viral replication. The frameshift mechanism is dependent upon a seven-nucleotide slippery sequence and a downstream RNA structure. In SIV, the downstream RNA structure has been proposed to be either a stem-loop or a pseudoknot. Here, we report the functional, structural and thermodynamic characterization of the SIV frameshift site RNA. Translational frameshift assays indicate that a stem-loop structure is sufficient to promote efficient frameshifting in vitro. NMR and thermodynamic studies of SIV RNA constructs of varying length further support the absence of any pseudoknot interaction and indicate the presence of a stable stem-loop structure. We determined the structure of the SIV frameshift-inducing RNA by NMR. The structure reveals a highly ordered 12 nucleotide loop containing a sheared G-A pair, cross-strand adenine stacking, two G-C base-pairs, and a novel CCC triloop turn. The loop structure and its high thermostability preclude pseudoknot formation. Sequence conservation and modeling studies suggest that HIV-2 RNA forms the same structure. We conclude that, like the main sub-groups of HIV-1, SIV and HIV-2 utilize stable stem-loop structures to function as a thermodynamic barrier to translation, thereby inducing ribosomal pausing and frameshifting.