Secondary structure of the central region of bacteriophage MS2 RNA. Conservation and biological significance

Effect of particulate matters on inactivation of bacteriophage MS2 under irradiation above 320 nm

The inactivation of bacteriophage MS2 under irradiation above 320 nm was investigated, focusing on different solution pH, ionic strength, and Suwannee River natural organic matter (SRNOM) concentrations when solutions contained organic or inorganic particle matters. Inorganic and organic particles were modeled using kaolinite (KAO) and Microcystis aeruginosa (MA), respectively. The results showed that the two types of particles influenced on MS2 inactivation under different conditions. The lower pH contributed to the greater MS2 aggregation within pH range of 3.0 to 8.0, leading to an increasing inactivation rate. The presence of KAO induced reactive oxygen species (ROS) under the action of irradiation above 320 nm, which promoted the inactivation of MS2. By comparison, the [Formula: see text] produced by MA after irradiation promoted the inactivation at pH < 6, whereas when the pH is ≥ 6, the inactivation effect of MS2 was lower than that of particle-free solution because MS2 was no longer aggregated and MA has a shading effect. In the presence of Na⁺ or Ca²⁺ cation, irradiation above 320 nm could not effectively inactivate the MS2 under particle-free solution. By comparison, KAO increased the inactivation efficiency as a photosensitizer. With the increase of Ca²⁺ concentration, MS2 was more easily adsorbed to MA than aggregation. Until Ca²⁺ concentration reached 20 mM, the inactivation effect in MA solution was enhanced. In the presence of SRNOM, the inactivation effect increased with the increase of SRNOM concentration. When the SRNOM was 20 mM, the inactivation increased in the particle-free solution due to the greater production of [Formula: see text]. Compared with the particle-free solution, the KAO and MA inactivation efficiency was lower.

A combination of mRNA features influence the efficiency of leaderless mRNA translation initiation

Bacterial translation is thought to initiate by base pairing of the 16S rRNA and the Shine-Dalgarno sequence in the mRNA's 5' untranslated region (UTR). However, transcriptomics has revealed that leaderless mRNAs, which completely lack any 5' UTR, are broadly distributed across bacteria and can initiate translation in the absence of the Shine-Dalgarno sequence. To investigate the mechanism of leaderless mRNA translation initiation, synthetic in vivo translation reporters were designed that systematically tested the effects of start codon accessibility, leader length, and start codon identity on leaderless mRNA translation initiation. Using these data, a simple computational model was built based on the combinatorial relationship of these mRNA features that can accurately classify leaderless mRNAs and predict the translation initiation efficiency of leaderless mRNAs. Thus, start codon accessibility, leader length, and start codon identity combine to define leaderless mRNA translation initiation in bacteria.

Functional resilience of activated sludge exposed to Bacillus globigii and bacteriophage MS2

This study evaluated the effect of two biocontaminant surrogates on the activity and performance of activated sludge. In the presence of bacteriophage MS2 at 3.2 × 10⁸ PFU/ml, the peak respiration rates varied between 9.1 and 10.5 mg O₂/gVSS/hr, generally similar to the rates observed in negative controls. MS2 did not alter the molar CO₂-produced-to-O₂ consumed stoichiometry observed during respiration. Similar results were observed for Bacillus globigii (BG). Ethanol, a potential co-contaminant associated with biocontamination incidents, can inhibit initial oxygen uptake. MS2 and BG both adsorbed to activated sludge; post-exposure viability was confirmed for BG but not for MS2. This study is the first to evaluate the effects of BG and MS2 on activated sludge, and it presents a protocol that can be used in operational situations.

ssRNA Phages: Life Cycle, Structure and Applications

ssRNA phages belonging to the family Leviviridae are among the tiniest viruses, infecting various Gram-negative bacteria by adsorption to their pilus structures. Due to their simplicity, they have been intensively studied as models for understanding various problems in molecular biology and virology. Several of the studied ssRNA characteristics, such as coat protein–RNA interactions and the ability to readily form virus-like particles in recombinant expression systems, have fueled many practical applications such as RNA labeling and tracking systems and vaccine development. In this chapter, we review the life cycle, structure and applications of these small yet fascinating viruses.

The p53 mRNA: an integral part of the cellular stress response

A large number of signalling pathways converge on p53 to induce different cellular stress responses that aim to promote cell cycle arrest and repair or, if the damage is too severe, to induce irreversible senescence or apoptosis. The differentiation of p53 activity towards specific cellular outcomes is tightly regulated via a hierarchical order of post-translational modifications and regulated protein-protein interactions. The mechanisms governing these processes provide a model for how cells optimize the genetic information for maximal diversity. The p53 mRNA also plays a role in this process and this review aims to illustrate how protein and RNA interactions throughout the p53 mRNA in response to different signalling pathways control RNA stability, translation efficiency or alternative initiation of translation. We also describe how a p53 mRNA platform shows riboswitch-like features and controls the rate of p53 synthesis, protein stability and modifications of the nascent p53 protein. A single cancer-derived synonymous mutation disrupts the folding of this platform and prevents p53 activation following DNA damage. The role of the p53 mRNA as a target for signalling pathways illustrates how mRNA sequences have co-evolved with the function of the encoded protein and sheds new light on the information hidden within mRNAs.

Pervasive allele-specific regulation on RNA decay in hybrid mice

Cellular RNA abundance is determined by both RNA transcription and decay. Therefore, change in RNA abundance, which can drive phenotypic diversity between different species, could arise from genetic variants affecting either process. However, previous studies in the evolution of RNA expression have been largely focused on transcription. Here, to globally investigate the effects of cis-regulatory divergence on RNA decay in mammals for the first time, we quantified allele-specific differences in RNA decay rates (ASD) in an F1 hybrid mouse. Out of 8,815 genes with sufficient data, we identified 621 genes exhibiting significant cis-divergence. Systematic analysis of these genes revealed that the genetic variants affecting microRNA binding and RNA secondary structures contribute to the observed divergences. Finally, we demonstrated that although the divergences in RNA abundance were predominantly determined by allelic differences in RNA transcription, most genes with significant ASD did not exhibit significant difference in RNA abundance. For these genes, the apparently compensatory effect between the allelic differences in RNA transcription and ASD suggests that changes in RNA decay could serve as important means to stabilize RNA abundances during mammalian evolution.

Protein-RNA Interactions in the Single-Stranded RNA Bacteriophages

Bacteriophages of the Leviviridae family are small viruses with short single-stranded RNA (ssRNA) genomes. Protein-RNA interactions play a key role throughout the phage life cycle, and all of the conserved phage proteins - the maturation protein, the coat protein and the replicase - are able to recognize specific structures in the RNA genome. The phage-coded replicase subunit associates with several host proteins to form a catalytically active complex. Recognition of the genomic RNA by the replicase complex is achieved in a remarkably complex manner that exploits the RNA-binding properties of host proteins and the particular three-dimensional structure of the phage genome. The coat protein recognizes a hairpin structure at the beginning of the replicase gene. The binding interaction serves to regulate the expression of the replicase gene and can be remarkably different in various ssRNA phages. The maturation protein is a minor structural component of the virion that binds to the genome, mediates attachment to the host and guides the genome into the cell. The maturation protein has two distinct RNA-binding surfaces that are in contact with different regions of the genome. The maturation and coat proteins also work together to ensure the encapsidation of the phage genome in new virus particles. In this chapter, the different ssRNA phage protein-RNA interactions, as well as some of their practical applications, are discussed in detail.

Direct Evidence for Packaging Signal-Mediated Assembly of Bacteriophage MS2

Using cross-linking coupled to matrix-assisted laser desorption/ionization mass spectrometry and CLIP-Seq sequencing, we determined the peptide and oligonucleotide sequences at the interfaces between the capsid proteins and the genomic RNA of bacteriophage MS2. The results suggest that the same coat protein (CP)-RNA and maturation protein (MP)-RNA interfaces are used in every viral particle. The portions of the viral RNA in contact with CP subunits span the genome, consistent with a large number of discrete and similar contacts within each particle. Many of these sites match previous predictions of the locations of multiple, dispersed and degenerate RNA sites with cognate CP affinity termed packaging signals (PSs). Chemical RNA footprinting was used to compare the secondary structures of protein-free genomic fragments and the RNA in the virion. Some PSs are partially present in protein-free RNA but others would need to refold from their dominant solution conformations to form the contacts identified in the virion. The RNA-binding peptides within the MP map to two sections of the N-terminal half of the protein. Comparison of MP sequences from related phages suggests a similar arrangement of RNA-binding sites, although these N-terminal regions have only limited sequence conservation. In contrast, the sequences of the C-termini are highly conserved, consistent with them encompassing pilin-binding domains required for initial contact with host cells. These results provide independent and unambiguous support for the assembly of MS2 virions via a PS-mediated mechanism involving a series of induced-fit viral protein interactions with RNA.

p53 mRNA and p53 Protein Structures Have Evolved Independently to Interact with MDM2

The p53 tumor suppressor and its key regulator MDM2 play essential roles in development, ageing, cancer, and cellular stress responses in mammals. Following DNA damage, MDM2 interacts with p53 mRNA in an ATM kinase-dependent fashion and stimulates p53 synthesis, whereas under normal conditions, MDM2 targets the p53 protein for degradation. The peptide- and RNA motifs that interact with MDM2 are encoded by the same conserved BOX-I sequence, but how these interactions have evolved is unknown. Here, we show that a temperature-sensitive structure in the invertebrate Ciona intestinalis (Ci) p53 mRNA controls its interaction with MDM2. We also show that a nonconserved flanking region of Ci-BOX-I domain prevents the p53-MDM2 protein-protein interaction. These results indicate that the temperature-regulated p53 mRNA-MDM2 interaction evolved to become kinase regulated in the mammalian DNA damage response. The data also suggest that the negative regulation of p53 by MDM2 via protein-protein interaction evolved in vertebrates following changes in the BOX-I flanking sequence.

Four RNA families with functional transient structures

Protein-coding and non-coding RNA transcripts perform a wide variety of cellular functions in diverse organisms. Several of their functional roles are expressed and modulated via RNA structure. A given transcript, however, can have more than a single functional RNA structure throughout its life, a fact which has been previously overlooked. Transient RNA structures, for example, are only present during specific time intervals and cellular conditions. We here introduce four RNA families with transient RNA structures that play distinct and diverse functional roles. Moreover, we show that these transient RNA structures are structurally well-defined and evolutionarily conserved. Since Rfam annotates one structure for each family, there is either no annotation for these transient structures or no such family. Thus, our alignments either significantly update and extend the existing Rfam families or introduce a new RNA family to Rfam. For each of the four RNA families, we compile a multiple-sequence alignment based on experimentally verified transient and dominant (dominant in terms of either the thermodynamic stability and/or attention received so far) RNA secondary structures using a combination of automated search via covariance model and manual curation. The first alignment is the Trp operon leader which regulates the operon transcription in response to tryptophan abundance through alternative structures. The second alignment is the HDV ribozyme which we extend to the 5' flanking sequence. This flanking sequence is involved in the regulation of the transcript's self-cleavage activity. The third alignment is the 5' UTR of the maturation protein from Levivirus which contains a transient structure that temporarily postpones the formation of the final inhibitory structure to allow translation of maturation protein. The fourth and last alignment is the SAM riboswitch which regulates the downstream gene expression by assuming alternative structures upon binding of SAM. All transient and dominant structures are mapped to our new alignments introduced here.

Exploring the architecture of viral RNA genomes

The genomes of RNA viruses contain local structural elements and long-range interactions that control various steps in virus replication. While many individual RNA elements have been characterized, it remains less clear how the structure and activity of such elements are integrated and regulated within the complex context of complete viral genomes. Recent technical advances, particularly the development of high-throughput solution structure mapping methods, have made secondary structural analysis of entire viral RNA genomes feasible. As a consequence, whole-genome structural models have been deduced for a number of plus-strand RNA viruses and retroviruses and these structures have provided intriguing functional and evolutionary insights into global genome architecture.

Insights into RNA structure and function from genome-wide studies

A comprehensive understanding of RNA structure will provide fundamental insights into the cellular function of both coding and non-coding RNAs. Although many RNA structures have been analysed by traditional biophysical and biochemical methods, the low-throughput nature of these approaches has prevented investigation of the vast majority of cellular transcripts. Triggered by advances in sequencing technology, genome-wide approaches for probing the transcriptome are beginning to reveal how RNA structure affects each step of protein expression and RNA stability. In this Review, we discuss the emerging relationships between RNA structure and the regulation of gene expression.

On the role of four small hairpins in the HIV-1 RNA genome

An RNA secondary structure model for the complete HIV-1 genome has recently been published based on SHAPE technology. Several well-known RNA motifs such as TAR and RRE were confirmed and numerous new structured motifs were described that may play important roles in virus replication. The 9 kb viral RNA genome is densely packed with many RNA hairpin motifs and the collective fold may play an important role in HIV-1 biology. We initially focused on 16 RNA hairpin motifs scattered along the viral genome. We considered conservation of these structures, despite sequence variation among virus isolates, as a first indication for a significant function. Four relatively small hairpins exhibited considerable structural conservation and were selected for experimental validation in virus replication assays. Mutations were introduced into the HIV-1 RNA genome to destabilize individual RNA structures without affecting the protein-coding properties (silent codon changes). No major virus replication defects were scored, suggesting that these four hairpin structures do not play essential roles in HIV-1 replication.

A two-stage mechanism of viral RNA compaction revealed by single molecule fluorescence

Long RNAs often exist as multiple conformers in equilibrium. For the genomes of single-stranded RNA viruses, one of these conformers must include a compacted state allowing the RNA to be confined within the virion. We have used single molecule fluorescence correlation spectroscopy to monitor the conformations of viral genomes and sub-fragments in the absence and presence of coat proteins. Cognate RNA-coat protein interactions in two model viruses cause a rapid collapse in the hydrodynamic radii of their respective RNAs. This is caused by protein binding at multiple sites on the RNA that facilitate additional protein-protein contacts. The collapsed species recruit further coat proteins to complete capsid assembly with great efficiency and fidelity. The specificity in RNA-coat protein interactions seen at single-molecule concentrations reflects the packaging selectivity seen for such viruses in vivo. This contrasts with many in vitro reassembly measurements performed at much higher concentrations. RNA compaction by coat protein or polycation binding are distinct processes, implying that defined RNA-coat protein contacts are required for assembly.

Incorporating global features of RNA motifs in predictions for an ensemble of secondary structures for encapsidated MS2 bacteriophage RNA

The secondary structure of encapsidated MS2 genomic RNA poses an interesting RNA folding challenge. Cryoelectron microscopy has demonstrated that encapsidated MS2 RNA is well-ordered. Models of MS2 assembly suggest that the RNA hairpin-protein interactions and the appropriate placement of hairpins in the MS2 RNA secondary structure can guide the formation of the correct icosahedral particle. The RNA hairpin motif that is recognized by the MS2 capsid protein dimers, however, is energetically unfavorable, and thus free energy predictions are biased against this motif. Computer programs called Crumple, Sliding Windows, and Assembly provide useful tools for prediction of viral RNA secondary structures when the traditional assumptions of RNA structure prediction by free energy minimization may not apply. These methods allow incorporation of global features of the RNA fold and motifs that are difficult to include directly in minimum free energy predictions. For example, with MS2 RNA the experimental data from SELEX experiments, crystallography, and theoretical calculations of the path for the series of hairpins can be incorporated in the RNA structure prediction, and thus the influence of free energy considerations can be modulated. This approach thoroughly explores conformational space and generates an ensemble of secondary structures. The predictions from this new approach can test hypotheses and models of viral assembly and guide construction of complete three-dimensional models of virus particles.

Selection for thermostability can lead to the emergence of mutational robustness in an RNA virus

Mutational robustness has important evolutionary implications, yet the mechanisms leading to its emergence remain poorly understood. One possibility is selection acting on a correlated trait, as for instance thermostability (plastogenetic congruence). Here, we examine the correlation between mutational robustness and thermostability in experimental populations of the RNA bacteriophage Qβ. Thermostable viruses evolved after only six serial passages in the presence of heat shocks, and genome sequencing suggested that thermostability can be conferred by several alternative mutations. To test whether thermostable viruses have increased mutational robustness, we performed additional passages in the presence of nitrous acid. Whereas in control lines this treatment produced the expected reduction in growth rate caused by the accumulation of deleterious mutations, thermostable viruses showed no such reduction, indicating that they are more resistant to mutagenesis. Our results suggest that selection for thermostability can lead to the emergence of mutational robustness driven by plastogenetic congruence. As temperature is a widespread selective pressure in nature, the mechanism described here may be relevant to the evolution of mutational robustness.

Mutually-induced conformational switching of RNA and coat protein underpins efficient assembly of a viral capsid

Single-stranded RNA viruses package their genomes into capsids enclosing fixed volumes. We assayed the ability of bacteriophage MS2 coat protein to package large, defined fragments of its genomic, single-stranded RNA. We show that the efficiency of packaging into a T=3 capsid in vitro is inversely proportional to RNA length, implying that there is a free-energy barrier to be overcome during assembly. All the RNAs examined have greater solution persistence lengths than the internal diameter of the capsid into which they become packaged, suggesting that protein-mediated RNA compaction must occur during assembly. Binding ethidium bromide to one of these RNA fragments, which would be expected to reduce its flexibility, severely inhibited packaging, consistent with this idea. Cryo-EM structures of the capsids assembled in these experiments with the sub-genomic RNAs show a layer of RNA density beneath the coat protein shell but lack density for the inner RNA shell seen in the wild-type virion. The inner layer is restored when full-length virion RNA is used in the assembly reaction, implying that it becomes ordered only when the capsid is filled, presumably because of the effects of steric and/or electrostatic repulsions. The cryo-EM results explain the length dependence of packaging. In addition, they show that for the sub-genomic fragments the strongest ordered RNA density occurs below the coat protein dimers forming the icosahedral 5-fold axes of the capsid. There is little such density beneath the proteins at the 2-fold axes, consistent with our model in which coat protein dimers binding to RNA stem-loops located at sites throughout the genome leads to switching of their preferred conformations, thus regulating the placement of the quasi-conformers needed to build the T=3 capsid. The data are consistent with mutual chaperoning of both RNA and coat protein conformations, partially explaining the ability of such viruses to assemble so rapidly and accurately.

Folding kinetics of large RNAs

We introduce here a heuristic approach to kinetic RNA folding that constructs secondary structures by stepwise combination of building blocks. These blocks correspond to subsequences and their thermodynamically optimal structures. These are determined by the standard dynamic programming approach to RNA folding. Folding trajectories are modeled at base-pair resolution using the Morgan-Higgs heuristic and a barrier tree-based heuristic to connect combinations of the local building blocks. Implemented in the program Kinwalker, the algorithm allows co-transcriptional folding and can be used to fold sequences of up to about 1500 nucleotides in length. A detailed comparison with several well-studied examples from the literature, including the delayed folding of bacteriophage cloverleaf structures, the adenine sensing riboswitch, and the hok RNA, shows an excellent agreement of predicted trajectories and experimental evidence. The software is available as part of the ViennaRNA Package.

Complete genome sequence of the broad host range single-stranded RNA phage PRR1 places it in the Levivirus genus with characteristics shared with Alloleviviruses

Single-stranded RNA (ssRNA) bacteriophages of the family Leviviridae infect gram-negative bacteria. They are restricted to a single host genus. Phage PRR1 is an exception, having a broad host range due to the promiscuity of the receptor encoded by the IncP plasmid. Here we report the complete genome sequence of PRR1. Three proteins homologous with those of other ssRNA phages, i.e., maturation, coat, and replicase proteins, were identified. A fourth protein has a lysis function. Comparison of PRR1 with other members of the Leviviridae family places PRR1 in the genus Levivirus with some characteristics more similar to those of members of the genus Allolevivirus.

Site-specific cleavage of MS2 RNA by a thermostable DNA-linked RNase H

A series of DNA-linked RNases H, in which the 15-mer DNA is cross-linked to the Thermus thermophilus RNase HI (TRNH) variants at positions 135, 136, 137 and 138, were constructed and analyzed for their abilities to cleave the complementary 15-mer RNA. Of these, that with the DNA adduct at position 135 most efficiently cleaved the RNA substrate, indicating that position 135 is the most appropriate cross-linking site among those examined. To examine whether DNA-linked RNase H also site-specifically cleaves a highly structured natural RNA, DNA-linked TRNHs with a series of DNA adducts varying in size at position 135 were constructed and analyzed for their abilities to cleave MS2 RNA. These DNA adducts were designed such that DNA-linked enzymes cleave MS2 RNA at a loop around residue 2790. Of the four DNA-linked TRNHs with the 8-, 12-, 16- and 20-mer DNA adducts, only that with the 16-mer DNA adduct efficiently and site-specifically cleaved MS2 RNA. Primer extension revealed that this DNA-linked TRNH cleaved MS2 RNA within the target sequence.

Peptide display on live MS2 phage: restrictions at the RNA genome level

The potential of the RNA phage MS2 to accommodate extra amino acids in its major coat protein has been examined. Accordingly, a pentapeptide was encoded in the genome as an N-terminal extension. In the MS2 crystal structure, this part of the coat protein forms a loop that extends from the outer surface of the icosahedral virion. At the RNA level, the insert forms a large loop at the top of an existing hairpin. This study shows that it is possible to maintain inserts in the coat protein of live phages. However, not all inserts were genetically stable. Some suffer deletions, while others underwent adaptation by base substitutions. Whether or not an insert is stable appears to be determined by the choice of the nucleic acid sequence used to encode the extra peptide. This effect was not caused by differential translation, because coat-protein synthesis was equal in wild-type and mutants. We conclude that the stability of the insert depends on the structure of the large RNA hairpin loop, as demonstrated by the fact that a single substitution can convert an unstable loop into a stable one.

The forces driving molecular evolution

Competitive replication among RNA or DNA molecules at linear and non-linear rates of propagation has been reviewed from the perspective of a recent physicochemical model of molecular evolution and the findings are applied to pre-replication, prebiotic and biological evolution. A system of competitively replicating molecules was seen to follow a path of least action on both its thermodynamic and kinetic branch, in evolving toward steady state kinetics and equilibrium for the nucleotide condensation reaction. Stable and unstable states of coexistence, between competing molecular species, arise at nonlinear rates of propagation, and they derive from an equilibrium between kinetic forces. The de novo formation of self-replicating RNA molecules involves damping of these scalar forces, error tolerance and RNA driven strand separation. Increases in sequence complexity in the transition to self-replication does not exceed the free energy dissipated in RNA synthesis. Retrodiction of metabolic pathways and phylogenetic evidence point to the occurrence of three pre-replication metabolic systems, driven by autocatalytic C-fixation cycles. Thermodynamic and kinetic factors led to the replication take over. Biological evolution was found to involve resource capture, in addition to competition for a shared resource.

Long-range translational coupling in single-stranded RNA bacteriophages: an evolutionary analysis

In coliphage MS2 RNA a long-distance interaction (LDI) between an internal segment of the upstream coat gene and the start region of the replicase gene prevents initiation of replicase synthesis in the absence of coat gene translation. Elongating ribosomes break up the repressor LDI and thus activate the hidden initiation site. Expression studies on partial MS2 cDNA clones identified base pairing between 1427-1433 and 1738-1744, the so-called Min Jou (MJ) interaction, as the molecular basis for the long-range coupling mechanism. Here, we examine the biological significance of this interaction for the control of replicase gene translation. The LDI was disrupted by mutations in the 3'-side and the evolutionary adaptation was monitored upon phage passaging. Two categories of pseudorevertants emerged. The first type had restored the MJ interaction but not necessarily the native sequence. The pseudorevertants of the second type acquired a compensatory substitution some 80 nt downstream of the MJ interaction that stabilizes an adjacent LDI. In one examined case we confirmed that the second site mutations had restored coat-replicase translational coupling. Our results show the importance of translational control for fitness of the phage. They also reveal that the structure that buries the replicase start extends to structure elements bordering the MJ interaction.

Rescue of the RNA phage genome from RNase III cleavage

The secondary structure of the RNA from the single-stranded RNA bacteriophages, like MS2 and Qb, has evolved to serve a variety of functions such as controlling gene expression, exposing binding sites for the replicase and capsid proteins, allowing strand separation and so forth. On the other hand, all of these foldings have to perform in bacterial cells in which various RNA splitting enzymes are present. We therefore examined whether phage RNA structure is under selective pressure by host RNases. Here we show this to be true for RNase III. A fully double-stranded hairpin of 17 bp, which is an RNase III target, was inserted into a non-coding region of the MS2 RNA genome. In an RNase III-host these phages survived but in wild-type bacteria they did not. Here the stem underwent Darwinian evolution to a structure that was no longer a substrate for RNase III. This was achieved in three different ways: (i) the perfect stem was maintained but shortened by removing all or most of the insert; (ii) the stem acquired suppressor mutations that replaced Watson-Crick base pairs by mismatches; (iii) the stem acquired small deletions or insertions that created bulges. These insertions consist of short stretches of non-templated A or U residues. Their origin is ascribed to polyadenylation at the site of the RNase III cut (in the + or - strand) either by Escherichia coli poly(A) polymerase or by idling MS2 replicase.

Evolutionary reconstruction of a hairpin deleted from the genome of an RNA virus

The intercistronic region between the maturation and coat-protein genes of RNA phage MS2 contains important regulatory and structural information. The sequence participates in two adjacent stem-loop structures, one of which, the coat-initiator hairpin, controls coat-gene translation and is thus under strong selection pressure. We have removed 19 out of the 23 nucleotides constituting the intercistronic region, thereby destroying the capacity of the phage to build the two hairpins. The deletion lowered coat-protein yield more than 1000-fold, and the titer of the infectious clone carrying the deletion dropped 10 orders of magnitude as compared with the wild type. Two types of revertants were recovered. One had, in two steps, recruited 18 new nucleotides that served to rebuild the two hairpins and the lost Shine-Dalgarno sequence. The other type had deleted an additional six nucleotides, which allowed the reconstruction of the Shine-Dalgarno sequence and the initiator hairpin, albeit by sacrificing the remnants of the other stem-loop. The results visualize the immense genetic repertoire created by, what appears as, random RNA recombination. It would seem that in this genetic ensemble every possible new RNA combination is represented.

Random removal of inserts from an RNA genome: selection against single-stranded RNA

We have monitored the evolution of insertions in two MS2 RNA regions of known secondary structure where coding pressure is negligible or absent. Base changes and shortening of the inserts proceed until the excessive nucleotides can be accommodated in the original structure. The stems of hairpins can be dramatically extended but the loops cannot, revealing natural selection against single-stranded RNA. The 3' end of the MS2 A-protein gene forms a small hairpin with an XbaI sequence in the loop. This site was used to insert XbaI fragments of various sizes. Phages produced by these MS2 cDNA clones were not wild type, nor had they retained the full insert. Instead, every revertant phage had trimmed the insert in a different way to leave a four- to seven-membered loop to the now extended stem. Similar results were obtained with inserts in the 5' untranslated region. The great number of different revertants obtained from a single starting mutant as well as sequence inspection of the crossover points suggest that the removal of redundant RNA occurs randomly. The only common feature among all revertants appears the potential to form a hairpin with a short loop, suggesting that single-stranded RNA negatively affects the viability of the phage. To test this hypothesis, we introduced XbaI fragments of 34 nucleotides that could form either a long stem with a small loop or a short stem with a large loop (26 nucleotides). The base-paired inserts were perfectly maintained for many generations, whereas the unpaired versions were quickly trimmed back to reduce the size of the loop. These data confirm that single-stranded RNA adversely affects phage fitness and is strongly selected against. The repair of the RNA genome that we describe here appears as the result of random recombination. Of the plethora of recombinants, only those able to adopt a base-paired structure survive. The frequency with which our inserts are removed seems higher than measured by others for small inserts in a reading frame in Q beta RNA. To account for this higher frequency, we suggest models in which the single-stranded nature of our inserts induces random recombination at the site of the insertion.

Psoralen crosslinking between human immunodeficiency virus type 1 RNA and primer tRNA3(Lys)

Initiation of reverse transcription is a crucial step of retroviral infection. In HIV-1, it involves hybridization of the 18 3'-terminal nucleotides of the primer tRNA3(Lys) to the primer binding site (PBS) of the viral RNA. Moreover, additional interactions between the two RNAs were recently evidenced [Isel et al. (1995) J. Mol. Biol. 247, 25269-25272]. To get further information on the topology of the viral RNA/tRNA3(Lys) complex, we used psoralen to induce RNA-RNA crosslinking. A defined intermolecular crosslinked complex was obtained. The crosslinked regions were characterized by RNase T1 digestion followed by bi-dimensional gel electrophoresis. The crosslinked residues (nucleotide mcm5S2U34 and U35 in the anticodon loop of tRNA3(Lys) and UCU154 in the viral RNA upstream of the PBS) were mapped using a retardation method coupled with random hydrolysis. The formation of this crosslink depends on the same elements that are required for the formation of the extended interactions between primer and template RNAs, i.e., the modified bases of the tRNA and a conserved A-rich loop located upstream of the PBS in the genomic RNA. Therefore, the present crosslinking data provide relevant information on the topology of the template/primer binary complex.

Secondary structure model of the Mason-Pfizer monkey virus 5' leader sequence: identification of a structural motif common to a variety of retroviruses

A stable secondary structure model is presented for the region 3' of the primer-binding site to 130 bases into the gag sequence of the prototype type D retrovirus Mason-Pfizer monkey virus. Using biochemical probing of RNA from this region in association with free energy minimization, we have identified a stem-loop structure in the region, which from other studies has been shown to be important for genomic RNA encapsidation. The structure involves a highly stable stem of five G-C pairs terminating in a heptaloop. Comparison of the Mason-Pfizer monkey virus structure with one predicted for squirrel monkey retrovirus demonstrates an identical stem and a common ACC motif in the loop. Free energy studies of the secondary structure of the 5' regions of eight other retroviruses predict stem loops which have similar GAYC motifs. We believe this may represent a common structural and sequence motif which among other functions may be involved in genomic RNA packaging in these viruses.

Nucleotide sequence of a single-stranded RNA phage from Pseudomonas aeruginosa: kinship to coliphages and conservation of regulatory RNA structures

We report the complete nucleotide sequence of the single-stranded RNA phage PP7 from Pseudomonas aeruginosa. There are three open reading frames which code for apparent protein homologues of the single-stranded RNA coliphages, i.e., maturation protein, coat protein, and replicase. A fourth overlapping reading frame exists that probably encodes a lysis protein, similar to what has been found in the group A coliphages such as MS2. The genetic map of PP7 is colinear with group A coliphages and we accordingly classify the phage as a levivirus. There is, generally speaking, no significant nucleotide sequence identity between PP7 and the coliphages except for a few regions where homologous parts of proteins are encoded, most notable in the replicase gene. In these regions the nucleotide sequence similarity between PP7 and MS2 is no greater than between PP7 and the group B coliphages such as Q beta. Surprisingly, Q beta and MS2 are no closer to each other than they are to PP7. Several regulatory RNA secondary structure features that are present in the coliphages were identified also in PP7 RNA although the sequences involved cannot be aligned. Among these are the coat protein binding helix at the start of the replicase gene, structures at the 5' and 3' terminus of the RNA, a replicase binding site, and the structure of the coat protein cistron start. Some of these features resemble MS2 type coliphages but others the Q beta type. These findings suggest that PP7 is related to the coliphages but branched off before the coliphages diverged into separate groups.

Translational initiation at the coat-protein gene of phage MS2: native upstream RNA relieves inhibition by local secondary structure

Maximal translation of the coat-protein gene from RNA bacteriophage MS2 requires a contiguous stretch of native MS2 RNA that extends hundreds of nucleotides upstream from the translational start site. Deletion of these upstream sequences from MS2 cDNA plasmids results in a 30-fold reduction of translational efficiency. By site-directed mutagenesis, we show that this low level of expression is caused by a hairpin structure centred around the initiation codon. When this hairpin is destabilized by the introduction of mismatches, expression from the truncated messenger increases 20-fold to almost the level of the full-length construct. Thus, the translational effect of hundreds of upstream nucleotides can be mimicked by a single substitution that destabilizes the structure. The same hairpin is also present in full-length MS2 RNA, but there it does not impair ribosome binding. Apparently, the upstream RNA somehow reduces the inhibitory effect of the structure on translational initiation. The upstream MS2 sequence does not stimulate translation when cloned in front of another gene, nor can unrelated RNA segments activate the coat-protein gene. Several possible mechanisms for the activation are discussed and a function in gene regulation of the phage is suggested.

Identification of a binding site for the human immunodeficiency virus type 1 nucleocapsid protein

The nucleocapsid (NC) protein NCp7 of human immunodeficiency virus type 1 (HIV-1) is important for encapsidation of the virus genome, RNA dimerization, and primer tRNA annealing in vitro. Here we present evidence from gel mobility-shift experiments indicating that NCp7 binds specifically to an RNA sequence. Two complexes were identified in native gels. The more slowly migrating complex contained two RNA molecules and one peptide, while the more rapidly migrating one is composed of one RNA and one peptide. Further, mutational analysis of the RNA shows that the predicted stem and loop structure of stem-loop 1 plays a critical role. Our results show that NCp7 binds to a unique RNA structure within the psi region; in addition, this structure is necessary for RNA dimerization. We propose that NCp7 binds to the RNA via a direct interaction of one zinc-binding motif to stem-loop 1 followed by binding of the other zinc-binding motif to stem-loop 1, stem-loop 2, or the linker region of the second RNA molecule, forming a bridge between the two RNAs.

Translational control by a long range RNA-RNA interaction; a basepair substitution analysis

One of the two mechanisms that regulate expression of the replicase cistron in the single stranded RNA coliphages is translational coupling. This mechanism prevents ribosomes from binding at the start of the replicase cistron unless the upstream coat cistron is being translated. Genetic analysis had identified a maximal region of 132 nucleotides in the coat gene over which ribosomes should pass to activate the replicase start. Subsequent deletion studies in our laboratory had further narrowed down the regulatory region to 12 nucleotides. Here, we identify a long-distance RNA-RNA interaction of 6 base pairs as the basis of the translational polarity. The 3' side of the complementarity region is located in the coat-replicase intercistronic region, some 20 nucleotides before the start codon of the replicase. The 5' side encodes amino acids 31 and 32 of the coat protein. Mutations that disrupt the long-range interaction abolish the translational coupling. Repair of basepairing by second site base substitutions restores translational coupling.

Secondary structure of the HIV-2 leader RNA comprising the tRNA-primer binding site

The initiation of reverse transcription of a retroviral RNA genome occurs by a tRNA primer bound near the 5' end of the genomic RNA at a position called the primer-binding site (PBS). To understand the molecular basis for this RNA-RNA interaction, the secondary structure of the leader RNA of the human immunodeficiency virus type 2 (HIV-2) RNA was analyzed. In vitro synthesized HIV-2 RNA was probed with various structure-specific enzymes and chemicals. A computer program was then used to predict the secondary structure consistent with these data. In addition, the nucleotide sequences of different HIV-2 isolates were used to screen for the occurrence of covariation among putative base pairs. The primary sequences have diverged rapidly in some HIV-2 isolates, however, some strikingly conserved secondary structure elements were identified. Most nucleotides in the leader region are involved in base pairing. An exception is the PBS sequence, of which 15 out of 18 nucleotides are exposed in an internal loop. These findings suggest that the overall structure of the HIV-2 genome has evolved to facilitate an optimal interaction with its tRNA primer.

Extraordinarily stable mini-hairpins: electrophoretical and thermal properties of the various sequence variants of d(GCGAAAGC) and their effect on DNA sequencing

A small DNA fragment having a characteristic sequence d(GCGAAAGC) has been shown to form an extraordinarily stable mini-hairpin structure and to have an unusually rapid mobility in polyacrylamide gel electrophoresis, even when containing 7M urea. Here, we have studied the stability of the various sequence variants of d(GCGAAAGC) and the corresponding RNA fragments. Many such sequence variants form stable mini-hairpins in a similar manner to the d(GCGAAAGC) sequence. The RNA fragment, r(GCGAAAGC) also forms a mini-hairpin structure with less stability. The DNA mini-hairpins with GAAA or GAA loop are much more stable than DNA and RNA mini-hairpins with other loop sequence so far as has been examined. The stability difference between DNA and RNA mini-hairpins may be deduced to the stem structures formed by DNA (B form) and RNA (A form). The stable hairpins consisting of the GCGAAAGC sequence cause strong band compression on the sequencing gel. This phenomenon should be carefully considered in DNA sequencing.

Effect of dimerization on the conformation of the encapsidation Psi domain of Moloney murine leukemia virus RNA

In Moloney murine leukemia virus, the encapsidation Psi element was shown to be necessary and sufficient to promote packaging of viral RNA, and to be required for dimerization. The conformation of the Psi domain (nucleotides 215 to 565) was investigated in solution by chemical probing. The four bases were monitored at one of their Watson-Crick positions with dimethylsulfate at cytosine N3 and adenosine N1, and with a carbodiimide derivative at guanosine N1 and uridine N3. Position N7 of adenine residues was probed with diethylpyrocarbonate. The analyses were conducted on in vitro transcribed fragments corresponding either to the isolated Psi domain or to the 5'-terminal 725 nucleotides. The RNA fragments were analyzed in their monomeric and dimeric forms. A secondary structure model was derived from probing data, computer prediction and sequence analysis of related murine retroviruses. One major result is that Psi forms an independent and highly structured domain. Dimerization induces an extensive reduction of reactivity in region 278 to 309 that can be interpreted as the result of intermolecular interactions and/or intramolecular conformational rearrangements. A second region (around position 215) was shown to display discrete reactivity changes upon dimerization. These two regions represent likely elements of dimerization. More unexpectedly, reactivity changes (essentially enhancement of reactivity) were also detected in another part of Psi (around position 480) not believed to contain elements of dimerization. These reactivity changes could be interpreted as dimerization-induced allosteric transitions.

Secondary structure of coliphage Q beta RNA. Analysis by electron microscopy

The secondary structure of genomic RNA from the coliphage Q beta has been examined by electron microscopy in the presence of varying concentrations of spermidine using the Kleinschmidt spreading technique. The size and position of structural features that cover 70% of the viral genome have been mapped. The structural features that are visualized by electron microscopy in Q beta RNA are large. They range in size from 170 to 1600 nucleotides. A loop containing approximately 450 nucleotides is located at the 5' end of the RNA. It includes the initiation region for the viral maturation protein. A large hairpin containing approximately 1600 nucleotides is located in the center of the molecule. It is multibranched and includes most of the viral coat gene, the readthrough region of the A1 gene, and approximately one third of the viral replicase gene. Within the central hairpin, the initiation region for the viral replicase gene pairs with a region within the distal third of the viral coat gene. This structure may participate in the regulation of translational initiation of the viral replicase gene. Two structural variants of the central hairpin were observed. One of them brings the internal S and M viral replicase binding regions into juxtaposition. These observations suggest that the central hairpin may also participate in the regulation of translation of the viral coat gene. The secondary structures that are observed in Q beta RNA differ significantly from structures that we described previously in the genomic RNA of coliphage MS2 but are similar to structures we observed by electron microscopy in the related group B coliphage SP.

Secondary structure of the ribosome binding site determines translational efficiency: a quantitative analysis

We have quantitatively analyzed the relationship between translational efficiency and the mRNA secondary structure in the initiation region. The stability of a defined hairpin structure containing a ribosome binding site was varied over 12 kcal/mol (1 cal = 4.184 J) by site-directed mutagenesis and the effects on protein yields were analyzed in vivo. The results reveal a strict correlation between translational efficiency and the stability of the helix. An increase in its delta G0 of -1.4 kcal/mol (i.e., less than the difference between an A.U and a G.C pair) corresponds to the reduction by a factor of 10 in initiation rate. Accordingly, a single nucleotide substitution led to the decrease by a factor of 500 in expression because it turned a mismatch in the helix into a match. We find no evidence that exposure of only the Shine-Dalgarno region or the start codon preferentially favors recognition. Translational efficiency is strictly correlated with the fraction of mRNA molecules in which the ribosome binding site is unfolded, indicating that initiation is completely dependent on spontaneous unfolding of the entire initiation region. Ribosomes appear not to recognize nucleotides outside the Shine-Dalgarno region and the initiation codon.

Complete nucleotide sequence of the group I RNA bacteriophage fr

We report the complete nucleotide sequence of the group I RNA bacteriophage fr. The entire genome consists of 3575 nucleotides, six nucleotides more than the only other sequenced group I representative, MS2. The greatest divergence between these phages occurs in the 5' terminal region of the A gene, while the lysis-replicase gene overlap, the coat gene and the central region of the replicase gene are highly conserved. Overall sequence homology between fr and MS2 is 77%. Here, we present a general comparison between the two phages. In the accompanying paper we use phylogenetic sequence comparison between MS2 and fr to deduce the secondary structure at the 3' untranslated region.

Secondary structure at the 3' terminal region of RNA coliphages: comparison with tRNA

Secondary structure models for the 3' non-coding region of the four groups of coliphage RNA are proposed based on comparative sequence analysis and on previously published data on the sensitivity of nucleotides in MS2 RNA to chemical modification and enzymes. We report the following observations. (1) In contrast to the coding regions, the structure at the 3' terminus is characterized by stable regular helices. We note the occurrence of the loop sequences 5'-GUUCGC and 5'-CGAAAG, that are reported to confer exceptional stability to stem structures. These features are probably present to promote the segregation of mother and daughter strands during replication. (2) Comparison of homologous helices indicates that only those base pair substitutions are allowed that maintain the thermodynamic stability. (3) We have compared the structure of phage RNA with tRNA. Overall similarity is low, but one common element may exist. It is a quasi-continuous helix of 12 basepairs that could be the equivalent of the 12 basepair long coaxially stacked helix, formed by the T psi C arm and the aminoacyl acceptor arm in tRNA. As in tRNA, this structure element starts after the fourth nucleotide from the 3' end. (4) Phage RNA contains a large variable region of about 35 nucleotides bulging out from the quasi-continuous helix. We speculate that the variable loop in present-day tRNA could be the remnant of the variable region found in phage RNA. The variable region contains overlapping binding sites for the replicase enzyme and the maturation protein. This common binding site may serve as a switch from replication to packaging.