SHAPE analysis of long-range interactions reveals extensive and thermodynamically preferred misfolding in a fragile group I intron RNA

Most functional RNAs require proteins to facilitate formation of their active structures. In the case of the yeast bI3 group I intron, splicing requires binding by two proteins, the intron-encoded bI3 maturase and the nuclear encoded Mrs1. Here, we use selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) chemistry coupled with analysis of point mutants to map long-range interactions in this RNA. This analysis reveals two critical features of the free RNA state. First, the catalytic intron is separated from the flanking exons via a stable anchoring helix. This anchoring helix creates an autonomous structural domain for the intron and functions to prevent misfolding with the flanking exons. Second, the thermodynamically most stable structure for the free RNA is not consistent with the catalytically active conformation as phylogenetically conserved elements form stable, non-native structures. These results highlight a fragile bI3 RNA for which binding by protein cofactors functions to promote extensive secondary structure rearrangements that are an obligatory prerequisite for forming the catalytically active tertiary structure.     

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/.     

HIV-1 TAR RNA subverts RNA interference in transfected cells through sequestration of TAR RNA-binding protein, TRBP

TAR RNA-binding protein, TRBP, was recently discovered to be an essential partner for Dicer and a crucial component of the RNA-induced silencing complex (RISC), a critical element of the RNA interference (RNAi) of the cell apparatus. Human TRBP was originally characterized and cloned 15 years ago based on its high affinity for binding the HIV-1 encoded leader RNA, TAR. RNAi is used, in part, by cells to defend against infection by viruses. Here, we report that transfected TAR RNA can attenuate the RNAi machinery in human cells. Our data suggest that TAR RNA sequesters TRBP rendering it unavailable for downstream Dicer-RISC complexes. TAR-induced inhibition of Dicer-RISC activity in transfected cells was partially relieved by exogenous expression of TRBP.     

Cardioviral RNA structure logo analysis: entropy, correlations, and prediction

In recent years, there has been an increased number of sequenced RNAs leading to the development of new RNA databases. Thus, predicting RNA structure from multiple alignments is an important issue to understand its function. Since RNA secondary structures are often conserved in evolution, developing methods to identify covariate sites in an alignment can be essential for discovering structural elements. Structure Logo is a technique established on the basis of entropy and mutual information measured to analyze RNA sequences from an alignment. We proposed an efficient Structure Logo approach to analyze conservations and correlations in a set of Cardioviral RNA sequences. The entropy and mutual information content were measured to examine the conservations and correlations, respectively. The conserved secondary structure motifs were predicted on the basis of the conservation and correlation analyses. Our predictive motifs were similar to the ones observed in the viral RNA structure database, and the correlations between bases also corresponded to the secondary structure in the database.     

TAR cloning: insights into gene function, long-range haplotypes and genome structure and evolution

The structural and functional analysis of mammalian genomes would benefit from the ability to isolate from multiple DNA samples any targeted chromosomal segment that is the size of an average human gene. A cloning technique that is based on transformation-associated recombination (TAR) in the yeast Saccharomyces cerevisiae satisfies this need. It is a unique tool to selectively recover chromosome segments that are up to 250 kb in length from complex genomes. In addition, TAR cloning can be used to characterize gene function and genome variation, including polymorphic structural rearrangements, mutations and the evolution of gene families, and for long-range haplotyping.     

Structural features in TAR RNA of human and simian immunodeficiency viruses: a phylogenetic analysis.

A comparative analysis of TAR RNA structures in human and simian immunodeficiency viruses reveals the conservation of certain structural features despite the divergence in sequence. Both the TAR elements of HIV-1 and SIV-chimpanzee can be folded into relatively simple one-stem hairpin structures. Chemical and RNAase probes were used to analyze the more complex structure of HIV-2 TAR RNA, which folds into a branched hairpin structure. A surprisingly similar RNA conformation can be proposed for SIV-mandrill, despite considerable divergence in nucleotide sequence. A third structural presentation of TAR sequences is seen for SIV-african green monkey. These results are generally consistent with the classification of HIV-SIV viruses in four subgroups based on sequence analyses (both nucleotide- and amino acid-sequences). However, some conserved TAR structures were detected for members of different virus subgroups. It is therefore proposed that RNA structure analysis might provide an additional tool for determining phylogenetic relationships among the HIV-SIV viruses.     

Three-dimensional analysis of a viral RNA replication complex reveals a virus-induced mini-organelle

Positive-strand RNA viruses are the largest genetic class of viruses and include many serious human pathogens. All positive-strand RNA viruses replicate their genomes in association with intracellular membrane rearrangements such as single- or double-membrane vesicles. However, the exact sites of RNA synthesis and crucial topological relationships between relevant membranes, vesicle interiors, surrounding lumens, and cytoplasm generally are poorly defined. We applied electron microscope tomography and complementary approaches to flock house virus (FHV)-infected Drosophila cells to provide the first 3-D analysis of such replication complexes. The sole FHV RNA replication factor, protein A, and FHV-specific 5-bromouridine 5'-triphosphate incorporation localized between inner and outer mitochondrial membranes inside approximately 50-nm vesicles (spherules), which thus are FHV-induced compartments for viral RNA synthesis. All such FHV spherules were outer mitochondrial membrane invaginations with interiors connected to the cytoplasm by a necked channel of approximately 10-nm diameter, which is sufficient for ribonucleotide import and product RNA export. Tomographic, biochemical, and other results imply that FHV spherules contain, on average, three RNA replication intermediates and an interior shell of approximately 100 membrane-spanning, self-interacting protein As. The results identify spherules as the site of protein A and nascent RNA accumulation and define spherule topology, dimensions, and stoichiometry to reveal the nature and many details of the organization and function of the FHV RNA replication complex. The resulting insights appear relevant to many other positive-strand RNA viruses and support recently proposed structural and likely evolutionary parallels with retrovirus and double-stranded RNA virus virions.     

Tensor decomposition reveals concurrent evolutionary convergences and divergences and correlations with structural motifs in ribosomal RNA

Evolutionary relationships among organisms are commonly described by using a hierarchy derived from comparisons of ribosomal RNA (rRNA) sequences. We propose that even on the level of a single rRNA molecule, an organism's evolution is composed of multiple pathways due to concurrent forces that act independently upon different rRNA degrees of freedom. Relationships among organisms are then compositions of coexisting pathway-dependent similarities and dissimilarities, which cannot be described by a single hierarchy. We computationally test this hypothesis in comparative analyses of 16S and 23S rRNA sequence alignments by using a tensor decomposition, i.e., a framework for modeling composite data. Each alignment is encoded in a cuboid, i.e., a third-order tensor, where nucleotides, positions and organisms, each represent a degree of freedom. A tensor mode-1 higher-order singular value decomposition (HOSVD) is formulated such that it separates each cuboid into combinations of patterns of nucleotide frequency variation across organisms and positions, i.e., "eigenpositions" and corresponding nucleotide-specific segments of "eigenorganisms," respectively, independent of a-priori knowledge of the taxonomic groups or rRNA structures. We find, in support of our hypothesis that, first, the significant eigenpositions reveal multiple similarities and dissimilarities among the taxonomic groups. Second, the corresponding eigenorganisms identify insertions or deletions of nucleotides exclusively conserved within the corresponding groups, that map out entire substructures and are enriched in adenosines, unpaired in the rRNA secondary structure, that participate in tertiary structure interactions. This demonstrates that structural motifs involved in rRNA folding and function are evolutionary degrees of freedom. Third, two previously unknown coexisting subgenic relationships between Microsporidia and Archaea are revealed in both the 16S and 23S rRNA alignments, a convergence and a divergence, conferred by insertions and deletions of these motifs, which cannot be described by a single hierarchy. This shows that mode-1 HOSVD modeling of rRNA alignments might be used to computationally predict evolutionary mechanisms.     

Pulse Dipolar ESR of Doubly Labeled Mini TAR DNA and Its Annealing to Mini TAR RNA

Pulse dipolar electron-spin resonance in the form of double electron electron resonance was applied to strategically placed, site-specifically attached pairs of nitroxide spin labels to monitor changes in the mini TAR DNA stem-loop structure brought on by the HIV-1 nucleocapsid protein NCp7. The biophysical structural evidence was at Ångstrom-level resolution under solution conditions not amenable to crystallography or NMR. In the absence of complementary TAR RNA, double labels located in both the upper and the lower stem of mini TAR DNA showed in the presence of NCp7 a broadened distance distribution between the points of attachment, and there was evidence for several conformers. Next, when equimolar amounts of mini TAR DNA and complementary mini TAR RNA were present, NCp7 enhanced the annealing of their stem-loop structures to form duplex DNA-RNA. When duplex TAR DNA-TAR RNA formed, double labels initially located 27.5 Å apart at the 3′- and 5′-termini of the 27-base mini TAR DNA relocated to opposite ends of a 27 bp RNA-DNA duplex with 76.5 Å between labels, a distance which was consistent with the distance between the two labels in a thermally annealed 27-bp TAR DNA-TAR RNA duplex. Different sets of double labels initially located 26–27 Å apart in the mini TAR DNA upper stem, appropriately altered their interlabel distance to ∼35 Å when a 27 bp TAR DNA-TAR RNA duplex formed, where the formation was caused either through NCp7-induced annealing or by thermal annealing. In summary, clear structural evidence was obtained for the fraying and destabilization brought on by NCp7 in its biochemical function as an annealing agent and for the detailed structural change from stem-loop to duplex RNA-DNA when complementary RNA was present.     

Binding of a cyclic BIV beta-Tat peptide with its TAR RNA construct

The ability of RNA structures to adopt diverse yet complex tertiary structures has resulted in numerous fascinating RNA-protein recognition events. It was recently reported that a close relative of the HIV Rev peptide, namely a 17 residue Tat peptide from bovine immuno-deficiency virus (BIV), is able to bind to the 28 nucleotide BIV TAR RNA construct. Here we report that by simply converting the 17 residue beta-ribbon peptide structure to a 19 residue cyclopeptide, the binding affinity (Kd) of the resulting cyclopeptide to the TAR RNA target, observed by fluorescence binding study, was enhanced approximately 5-fold.     

Recognition and long-range interactions of a minimal nanos RNA localization signal element

Localization of nanos (nos) mRNA to the germ plasm at the posterior pole of the Drosophila embryo is essential to activate nos translation and thereby generate abdominal segments. nos RNA localization is mediated by a large cis-acting localization signal composed of multiple, partially redundant elements within the nos 3' untranslated region. We identify a protein of approximately 75 kDa (p75) that interacts specifically with the nos +2' localization signal element. We show that the function of this element can be delimited to a 41 nucleotide domain that is conserved between D. melanogaster and D. virilis, and confers near wild-type localization when present in three copies. Two small mutations within this domain eliminate both +2' element localization function and p75 binding, consistent with a role for p75 in nos RNA localization. In the intact localization signal, the +2' element collaborates with adjacent localization elements. We show that different +2' element mutations not only abolish collaboration between the +2' and adjacent +1 element but also produce long-range deleterious effects on localization signal function. Our results suggest that higher order structural interactions within the localization signal, which requires factors such as p75, are necessary for association of nos mRNA with the germ plasm.     

Identification of a novel HIV-1 TAR RNA bulge binding protein.

The Tat protein binds to TAR RNA to stimulate the expression of the human immunodeficiency virus type 1 (HIV-1) genome. Tat is an 86 amino acid protein that contains a short region of basic residues (aa49-aa57) that are required for RNA binding and TAR is a 59 nucleotide stem-loop with a tripyrimidine bulge in the upper stem. TAR is located at the 5' end of all viral RNAs. In vitro, Tat specifically interacts with TAR by recognising the sequence of the bulge and upper stem, with no requirement for the loop. However, in vivo the loop sequence is critical for activation, implying a requirement for accessory cellular TAR RNA binding factors. A number of TAR binding cellular factors have been identified in cell extracts and various models for the function of these factors have been suggested, including roles as coactivators and inhibitors. We have now identified a novel 38 kD cellular factor that has little general, single-stranded or double-stranded RNA binding activity, but that specifically recognises the bulge and upper stem region of TAR. The protein, referred to as BBP (bulge binding protein), is conserved in mammalian and amphibian cells and in Schizosaccharomyces pombe but is not found in Saccharomyces cerevisiae. BBP is an effective competitive inhibitor of Tat binding to TAR in vitro. Our data suggest that the bulge-stem recognition motif in TAR is used to mediate cellular factor/RNA interactions and indicates that Tat action might be inhibited by such competing reactions in vivo.     

Probing the mechanical folding kinetics of TAR RNA by hopping, force-jump, and force-ramp methods

Mechanical unfolding and refolding of single RNA molecules have previously been observed in optical traps as sudden changes in molecular extension. Two methods have been traditionally used: "force-ramp", with the applied force continuously changing, and "hopping". In hopping experiments the force is held constant and the molecule jumps spontaneously between two different states. Unfolding/refolding rates are measured directly, but only over a very narrow range of forces. We have now developed a force-jump method to measure the unfolding and refolding rates independently over a wider range of forces. In this method, the applied force is rapidly stepped to a new value and either the unfolding or refolding event is monitored through changes in the molecular extension. The force-jump technique is compared to the force-ramp and hopping methods by using a 52-nucleotide RNA hairpin with a three-nucleotide bulge, i.e., the transactivation response region RNA from the human immunodeficiency virus. We find the unfolding kinetics and Gibbs free energies obtained from all three methods to be in good agreement. The transactivation response region RNA hairpin unfolds in an all-or-none two-state reaction at any loading rate with the force-ramp method. The unfolding reaction is reversible at small loading rates, but shows hysteresis at higher loading rates. Although the RNA unfolds and refolds without detectable intermediates in constant-force conditions (hopping and force-jump), it shows partially folded intermediates in force-ramp experiments at higher unloading rates. Thus, we find that folding of RNA hairpins can be more complex than a simple single-step reaction, and that application of several methods can improve understanding of reaction mechanisms.     

The dynamical implications of disease interference: correlations and coexistence

Ecological interference between unrelated diseases, caused by the temporary or permanent removal of individuals susceptible to one disease following infection with another, might be an important mechanism underlying epidemics. In this paper, we explore the potential dynamic consequences of interference by analyzing a two-disease model. By studying the stability domain of the model's equilibria, we find that the stable region of the two-disease endemic state becomes increasingly smaller as the strength of interference (largely determined by the disease-induced mortality) increases. When seasonal changes are included in the transmission rates, the bifurcation structure of the model's periodic cycles reveals that when the two diseases have similar mean transmission rates, multiple attractors in which the two diseases are strongly correlated can coexist, and that when the two diseases have very different mean transmission rates, the one with higher mean transmission rate may determine the dynamics of the system, with the other infection mimicking the behavior. We conclude that ecological interference can have important effects on the dynamical pattern of interacting diseases, the extent of which is determined by the epidemiological features of the diseases, their mean transmission rates in particular.