Ensemble of secondary structures for encapsidated satellite tobacco mosaic virus RNA consistent with chemical probing and crystallography constraints

Measuring intramolecular connectivity in long RNA molecules using two-dimensional DNA patch-probe arrays

We describe a simple method to infer intramolecular connections in a population of long RNA molecules in vitro. First we add DNA oligonucleotide "patches" that perturb the RNA connections, then we use a microarray containing a complete set of DNA oligonucleotide "probes" to record where perturbations occur. The pattern of perturbations reveals couplings between different regions of the RNA sequence, from which we infer connections as well as their prevalences in the population. We validate this patch-probe method using the 1,058-nucleotide RNA genome of satellite tobacco mosaic virus (STMV), which has previously been shown to have multiple long-range connections. Our results not only indicate long duplexes that agree with previous structures but also reveal the prevalence of competing connections. Together, these results suggest that globally-folded and locally-folded structures coexist in solution. We show that the prevalence of connections changes when pseudouridine, an important component of natural and synthetic RNA molecules, is substituted for uridine in STMV RNA.

Discovery of RNA-targeted small molecules through the merging of experimental and computational technologies

INTRODUCTION

The field of RNA-targeted small molecules is rapidly evolving, owing to the advances in experimental and computational technologies. With the identification of several bioactive small molecules that target RNA, including the FDA-approved risdiplam, the biopharmaceutical industry is gaining confidence in the field. This review, based on the literature obtained from PubMed, aims to disseminate information about the various technologies developed for targeting RNA with small molecules and propose areas for improvement to develop drugs more efficiently, particularly those linked to diseases with unmet medical needs.

AREAS COVERED

The technologies for the identification of RNA targets, screening of chemical libraries against RNA, assessing the bioactivity and target engagement of the hit compounds, structure determination, and hit-to-lead optimization are reviewed. Along with the description of the technologies, their strengths, limitations, and examples of how they can impact drug discovery are provided.

EXPERT OPINION

Many existing technologies employed for protein targets have been repurposed for use in the discovery of RNA-targeted small molecules. In addition, technologies tailored for RNA targets have been developed. Nevertheless, more improvements are necessary, such as artificial intelligence to dissect important RNA structures and RNA-small-molecule interactions and more powerful chemical probing and structure prediction techniques.

Genome-regulated Assembly of a ssRNA Virus May Also Prepare It for Infection

Many single-stranded, positive-sense RNA viruses regulate assembly of their infectious virions by forming multiple, cognate coat protein (CP)-genome contacts at sites termed Packaging Signals (PSs). We have determined the secondary structures of the bacteriophage MS2 ssRNA genome (gRNA) frozen in defined states using constraints from X-ray synchrotron footprinting (XRF). Comparison of the footprints from phage and transcript confirms the presence of multiple PSs in contact with CP dimers in the former. This is also true for a virus-like particle (VLP) assembled around the gRNA in vitro in the absence of the single-copy Maturation Protein (MP) found in phage. Since PS folds are present at many sites across gRNA transcripts, it appears that this genome has evolved to facilitate this mechanism of assembly regulation. There are striking differences between the gRNA-CP contacts seen in phage and the VLP, suggesting that the latter are inappropriate surrogates for aspects of phage structure/function. Roughly 50% of potential PS sites in the gRNA are not in contact with the protein shell of phage. However, many of these sit adjacent to, albeit not in contact with, PS-binding sites on CP dimers. We hypothesize that these act as PSs transiently during assembly but subsequently dissociate. Combining the XRF data with PS locations from an asymmetric cryo-EM reconstruction suggests that the genome positions of such dissociations are non-random and may facilitate infection. The loss of many PS-CP interactions towards the 3' end of the gRNA would allow this part of the genome to transit more easily through the narrow basal body of the pilus extruding machinery. This is the known first step in phage infection. In addition, each PS-CP dissociation event leaves the protein partner trapped in a non-lowest free-energy conformation. This destabilizes the protein shell which must disassemble during infection, further facilitating this stage of the life-cycle.

How and why RNA genomes are (partially) ordered in viral capsids

There is a long and productive progression of X-ray crystallographic and electron microscopy studies establishing the structures of the spherical/icosahedral and cylindrical/helical capsids of a wide range of virus particles. This is because of the high degree of order - down to the Angstrom scale - in the secondary/tertiary/quaternary structure of the proteins making up the capsids. In stark contradistinction, very little is known about the structure of DNA or RNA genomes inside these capsids. This is because of the relatively large extent of disorder in the confined DNA or RNA, due to several fundamental reasons: topological defects in the DNA case, and secondary/tertiary structural disorder in the RNA case. In this article we discuss the range of partial order associated with the encapsidated genomes of single-stranded RNA viruses, focusing on the contrast between mono-partite and multi-partite viruses and on the effects of sequence-specific and non-specific interactions between RNA and capsid proteins.

Universal and strain specific structure features of segment 8 genomic RNA of influenza A virus-application of 4-thiouridine photocrosslinking

RNA structure in the influenza A virus (IAV) has been the focus of several studies that have shown connections between conserved secondary structure motifs and their biological function in the virus replication cycle. Questions have arisen on how to best recognize and understand the pandemic properties of IAV strains from an RNA perspective, but determination of the RNA secondary structure has been challenging. Herein, we used chemical mapping to determine the secondary structure of segment 8 viral RNA (vRNA) of the pandemic A/California/04/2009 (H1N1) strain of IAV. Additionally, this long, naturally occurring RNA served as a model to evaluate RNA mapping with 4-thiouridine (4sU) crosslinking. We explored 4-thiouridine as a probe of nucleotides in close proximity, through its incorporation into newly transcribed RNA and subsequent photoactivation. RNA secondary structural features both universal to type A strains and unique to the A/California/04/2009 (H1N1) strain were recognized. 4sU mapping confirmed and facilitated RNA structure prediction, according to several rules: 4sU photocross-linking forms efficiently in the double-stranded region of RNA with some flexibility, in the ends of helices, and across bulges and loops when their structural mobility is permitted. This method highlighted three-dimensional properties of segment 8 vRNA secondary structure motifs and allowed to propose several long-range three-dimensional interactions. 4sU mapping combined with chemical mapping and bioinformatic analysis could be used to enhance the RNA structure determination as well as recognition of target regions for antisense strategies or viral RNA detection.

Perspectives on Viral RNA Genomes and the RNA Folding Problem

Viral RNA genomes change shape as virus particles disassemble, form replication complexes, attach to ribosomes for translation, evade host defense mechanisms, and assemble new virus particles. These structurally dynamic RNA shapeshifters present a challenging RNA folding problem, because the RNA sequence adopts multiple structures and may sometimes contain regions of partial disorder. Recent advances in high resolution asymmetric cryoelectron microscopy and chemical probing provide new ways to probe the degree of structure and disorder, and have identified more than one conformation in dynamic equilibrium in viral RNA. Chemical probing and the Detection of RNA Folding Ensembles using Expectation Maximization (DREEM) algorithm has been applied to studies of the dynamic equilibrium conformations in HIV RNA in vitro, in virio, and in vivo. This new type of data provides insight into important questions about virus assembly mechanisms and the fundamental physical forces driving virus particle assembly.

Challenges and approaches to predicting RNA with multiple functional structures

The revolution in sequencing technology demands new tools to interpret the genetic code. As in vivo transcriptome-wide chemical probing techniques advance, new challenges emerge in the RNA folding problem. The emphasis on one sequence folding into a single minimum free energy structure is fading as a new focus develops on generating RNA structural ensembles and identifying functional structural features in ensembles. This review describes an efficient combinatorially complete method and three free energy minimization approaches to predicting RNA structures with more than one functional fold, as well as two methods for analysis of a thermodynamics-based Boltzmann ensemble of structures. The review then highlights two examples of viral RNA 3'-UTR regions that fold into more than one conformation and have been characterized by single molecule fluorescence energy resonance transfer or NMR spectroscopy. These examples highlight the different approaches and challenges in predicting structure and function from sequence for RNA with multiple biological roles and folds. More well-defined examples and new metrics for measuring differences in RNA structures will guide future improvements in prediction of RNA structure and function from sequence.

Swellix: a computational tool to explore RNA conformational space

Background

The sequence of nucleotides in an RNA determines the possible base pairs for an RNA fold and thus also determines the overall shape and function of an RNA. The Swellix program presented here combines a helix abstraction with a combinatorial approach to the RNA folding problem in order to compute all possible non-pseudoknotted RNA structures for RNA sequences. The Swellix program builds on the Crumple program and can include experimental constraints on global RNA structures such as the minimum number and lengths of helices from crystallography, cryoelectron microscopy, or in vivo crosslinking and chemical probing methods.

Results

The conceptual advance in Swellix is to count helices and generate all possible combinations of helices rather than counting and combining base pairs. Swellix bundles similar helices and includes improvements in memory use and efficient parallelization. Biological applications of Swellix are demonstrated by computing the reduction in conformational space and entropy due to naturally modified nucleotides in tRNA sequences and by motif searches in Human Endogenous Retroviral (HERV) RNA sequences. The Swellix motif search reveals occurrences of protein and drug binding motifs in the HERV RNA ensemble that do not occur in minimum free energy or centroid predicted structures.

Conclusions

Swellix presents significant improvements over Crumple in terms of efficiency and memory use. The efficient parallelization of Swellix enables the computation of sequences as long as 418 nucleotides with sufficient experimental constraints. Thus, Swellix provides a practical alternative to free energy minimization tools when multiple structures, kinetically determined structures, or complex RNA-RNA and RNA-protein interactions are present in an RNA folding problem.

Electronic supplementary material

The online version of this article (10.1186/s12859-017-1910-7) contains supplementary material, which is available to authorized users.

Influenza virus segment 5 (+)RNA - secondary structure and new targets for antiviral strategies

Influenza A virus is a threat for humans due to seasonal epidemics and occasional pandemics. This virus can generate new strains that are dangerous through nucleotide/amino acid changes or through segmental recombination of the viral RNA genome. It is important to gain wider knowledge about influenza virus RNA to create new strategies for drugs that will inhibit its spread. Here, we present the experimentally determined secondary structure of the influenza segment 5 (+)RNA. Two RNAs were studied: the full-length segment 5 (+)RNA and a shorter construct containing only the coding region. Chemical mapping data combined with thermodynamic energy minimization were used in secondary structure prediction. Sequence/structure analysis showed that the determined secondary structure of segment 5 (+)RNA is mostly conserved between influenza virus type A strains. Microarray mapping and RNase H cleavage identified accessible sites for oligonucleotides in the revealed secondary structure of segment 5 (+)RNA. Antisense oligonucleotides were designed based on the secondary structure model and tested against influenza virus in cell culture. Inhibition of influenza virus proliferation was noticed, identifying good targets for antisense strategies. Effective target sites fall within two domains, which are conserved in sequence/structure indicating their importance to the virus.

Packaged and Free Satellite Tobacco Mosaic Virus (STMV) RNA Genomes Adopt Distinct Conformational States

The RNA genomes of viruses likely undergo multiple functionally important conformational changes during their replication cycles, changes that are poorly understood at present. We used two complementary in-solution RNA structure probing strategies (SHAPE-MaP and RING-MaP) to examine the structure of the RNA genome of satellite tobacco mosaic virus inside authentic virions and in a capsid-free state. Both RNA states feature similar three-domain architectures in which each major replicative function-translation, capsid coding, and genome synthesis-fall into distinct domains. There are, however, large conformational differences between the in-virion and capsid-free states, primarily in one arm of the central T domain. These data support a model in which the packaged capsid-bound RNA is constrained in a local high-energy conformation by the native capsid shell. The removal of the viral capsid then allows the RNA genome to relax into a more thermodynamically stable conformation. The RNA architecture of the central T domain thus likely changes during capsid assembly and disassembly and may play a role in genome packaging.

Secondary structure model of the naked segment 7 influenza A virus genomic RNA

The influenza A virus (IAV) genome comprises eight negative-sense viral (v)RNA segments. The seventh segment of the genome encodes two essential viral proteins and is specifically packaged alongside the other seven vRNAs. To gain insights into the possible roles of RNA structure both within and without virions, a secondary structure model of a naked (protein-free) segment 7 vRNA (vRNA7) has been determined using chemical mapping and thermodynamic energy minimization. The proposed structure model was validated using microarray mapping, RNase H cleavage and comparative sequence analysis. Additionally, the detailed structures of three vRNA7 fragment constructs - comprising independently folded subdomains - were determined. Much of the proposed vRNA7 structure is preserved between IAV strains, suggesting their importance in the influenza replication cycle. Possible structure rearrangements, which allow or preclude long-range RNA interactions, are also proposed.

Self-Folding of Naked Segment 8 Genomic RNA of Influenza A Virus

Influenza A is a negative sense RNA virus that kills hundreds of thousands of humans each year. Base pairing in RNA is very favorable, but possibilities for RNA secondary structure of the influenza genomic RNA have not been investigated. This work presents the first experimentally-derived exploration of potential secondary structure in an influenza A naked (protein-free) genomic segment. Favorable folding regions are revealed by in vitro chemical structure mapping, thermodynamics, bioinformatics, and binding to isoenergetic microarrays of an entire natural sequence of the 875 nt segment 8 vRNA and of a smaller fragment. Segment 8 has thermodynamically stable and evolutionarily conserved RNA structure and encodes essential viral proteins NEP and NS1. This suggests that vRNA self-folding may generate helixes and loops that are important at one or more stages of the influenza life cycle.

Crumple: An Efficient Tool to Explore Thoroughly the RNA Folding Landscape

The folding landscape for an RNA sequence contains many diverse structures and motifs, which are often sampled rather than completely explored. Today's supercomputers make the complete enumeration of all possible folds for an RNA and a detailed description of the RNA folding landscape a more feasible task. This chapter provides protocols for using the Crumple folding algorithm, an efficient tool to generate all possible non-pseudoknotted folds for an RNA sequence. Crumple in conjunction with Sliding Windows and Assembly can incorporate experimental constraints on the global features of an RNA, such as the minimum number and lengths of helices, which may be determined by crystallography or cryo-electron microscopy. This complete enumeration method is independent of free-energy minimization and allows the user to incorporate experimental data such as chemical probing, SELEX data on RNA-protein binding motifs, and phylogenetic covariation.

Structures and Energetics of Four Adjacent G·U Pairs That Stabilize an RNA Helix

Consecutive G·U base pairs inside RNA helices can be destabilizing, while those at the ends of helices are thermodynamically stabilizing. To determine if this paradox could be explained by differences in base stacking, we determined the high-resolution (1.32 Å) crystal structure of (5'-GGUGGCUGUU-3')2 and studied three sequences with four consecutive terminal G·U pairs by NMR spectroscopy. In the crystal structure of (5'-GGUGGCUGUU-3')2, the helix is overwound but retains the overall features of A-form RNA. The penultimate base steps at each end of the helix have high base overlap and contribute to the unexpectedly favorable energetic contribution for the 5'-GU-3'/3'-UG-5' motif in this helix position. The balance of base stacking and helical twist contributes to the positional dependence of G·U pair stabilities. The energetic stabilities and similarity to A-form RNA helices suggest that consecutive G·U pairs would be recognized by RNA helix binding proteins, such as Dicer and Ago. Thus, these results will aid future searches for target sites of small RNAs in gene regulation.

Biologically-supported structural model for a viral satellite RNA

Satellite RNAs (satRNAs) are a class of small parasitic RNA replicon that associate with different viruses, including plus-strand RNA viruses. Because satRNAs do not encode a polymerase or capsid subunit, they rely on a companion virus to provide these proteins for their RNA replication and packaging. SatRNAs recruit these and other required factors via their RNA sequences and structures. Here, through a combination of chemical probing analysis of RNA structure, phylogenetic structural comparisons, and viability assays of satRNA mutants in infected cells, the biological importance of a deduced higher-order structure for a 619 nt long tombusvirus satRNA was assessed. Functionally-relevant secondary and tertiary RNA structures were identified throughout the length of the satRNA. Notably, a 3′-terminal segment was found to adopt two mutually-exclusive RNA secondary structures, both of which were required for efficient satRNA accumulation. Accordingly, these alternative conformations likely function as a type of RNA switch. The RNA switch was also found to engage in a required long-range kissing-loop interaction with an upstream sequence. Collectively, these results establish a high level of conformational complexity within this small parasitic RNA and provide a valuable structural framework for detailed mechanistic studies.

Properties of satellite tobacco mosaic virus phenotypes expressed in the presence and absence of helper virus

In this study, we assembled an Agrobacterium-based transient expression system for the ectopic expression of Satellite tobacco mosaic virus (STMV) (+) or (-) transcripts and their biological activity was confirmed when Nicotiana benthamiana plants were co-expressed with helper Tobacco mosaic virus replicase. Characterization of STMV in the presence and absence of its HV revealed: (i) HV-dependent expression of STMV (+) in N. benthamiana, but not in N. tabacum, generated a replication-deficient but translation and encapsidation competent variant lacking the highly conserved 3' 150 nucleotides (nt) (STMVΔ150); (ii) mutational analysis demonstrated that a conserved 3' stem-loop structure in wild type and STMVΔ150 located between nt 874 and 897 is essential for translation of CP; (iii) helper virus-independent expression of CP from wt STMV was competent for the assembly of empty aberrant virion-like particles; whereas, CP translated from STMVΔ150 resulted in disorganized CP aggregates suggesting a role for the 3'tRNA-like structure in STMV assembly.

The Role of Packaging Sites in Efficient and Specific Virus Assembly

During the life cycle of many single-stranded RNA viruses, including many human pathogens, a protein shell called the capsid spontaneously assembles around the viral genome. Understanding the mechanisms by which capsid proteins selectively assemble around the viral RNA amidst diverse host RNAs is a key question in virology. In one proposed mechanism, short sequences (packaging sites) within the genomic RNA promote rapid and efficient assembly through specific interactions with the capsid proteins. In this work, we develop a coarse-grained particle-based computational model for capsid proteins and RNA that represents protein-RNA interactions arising both from nonspecific electrostatics and from specific packaging site interactions. Using Brownian dynamics simulations, we explore how the efficiency and specificity of assembly depend on solution conditions (which control protein-protein and nonspecific protein-RNA interactions) and the strength and number of packaging sites. We identify distinct regions in parameter space in which packaging sites lead to highly specific assembly via different mechanisms and others in which packaging sites lead to kinetic traps. We relate these computational predictions to in vitro assays for specificity in which cognate viral RNAs compete against non-cognate RNAs for assembly by capsid proteins.

Visualizing the global secondary structure of a viral RNA genome with cryo-electron microscopy

The lifecycle, and therefore the virulence, of single-stranded (ss)-RNA viruses is regulated not only by their particular protein gene products, but also by the secondary and tertiary structure of their genomes. The secondary structure of the entire genomic RNA of satellite tobacco mosaic virus (STMV) was recently determined by selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE). The SHAPE analysis suggested a single highly extended secondary structure with much less branching than occurs in the ensemble of structures predicted by purely thermodynamic algorithms. Here we examine the solution-equilibrated STMV genome by direct visualization with cryo-electron microscopy (cryo-EM), using an RNA of similar length transcribed from the yeast genome as a control. The cryo-EM data reveal an ensemble of branching patterns that are collectively consistent with the SHAPE-derived secondary structure model. Thus, our results both elucidate the statistical nature of the secondary structure of large ss-RNAs and give visual support for modern RNA structure determination methods. Additionally, this work introduces cryo-EM as a means to distinguish between competing secondary structure models if the models differ significantly in terms of the number and/or length of branches. Furthermore, with the latest advances in cryo-EM technology, we suggest the possibility of developing methods that incorporate restraints from cryo-EM into the next generation of algorithms for the determination of RNA secondary and tertiary structures.

Mechanisms of virus assembly

Viruses are nanoscale entities containing a nucleic acid genome encased in a protein shell called a capsid and in some cases are surrounded by a lipid bilayer membrane. This review summarizes the physics that govern the processes by which capsids assemble within their host cells and in vitro. We describe the thermodynamics and kinetics for the assembly of protein subunits into icosahedral capsid shells and how these are modified in cases in which the capsid assembles around a nucleic acid or on a lipid bilayer. We present experimental and theoretical techniques used to characterize capsid assembly, and we highlight aspects of virus assembly that are likely to receive significant attention in the near future.

Asymmetric genome organization in an RNA virus revealed via graph-theoretical analysis of tomographic data

Cryo-electron microscopy permits 3-D structures of viral pathogens to be determined in remarkable detail. In particular, the protein containers encapsulating viral genomes have been determined to high resolution using symmetry averaging techniques that exploit the icosahedral architecture seen in many viruses. By contrast, structure determination of asymmetric components remains a challenge, and novel analysis methods are required to reveal such features and characterize their functional roles during infection. Motivated by the important, cooperative roles of viral genomes in the assembly of single-stranded RNA viruses, we have developed a new analysis method that reveals the asymmetric structural organization of viral genomes in proximity to the capsid in such viruses. The method uses geometric constraints on genome organization, formulated based on knowledge of icosahedrally-averaged reconstructions and the roles of the RNA-capsid protein contacts, to analyse cryo-electron tomographic data. We apply this method to the low-resolution tomographic data of a model virus and infer the unique asymmetric organization of its genome in contact with the protein shell of the capsid. This opens unprecedented opportunities to analyse viral genomes, revealing conserved structural features and mechanisms that can be targeted in antiviral drug design.

A parallel implementation of the Wuchty algorithm with additional experimental filters to more thoroughly explore RNA conformational space

We present new modifications to the Wuchty algorithm in order to better define and explore possible conformations for an RNA sequence. The new features, including parallelization, energy-independent lonely pair constraints, context-dependent chemical probing constraints, helix filters, and optional multibranch loops, provide useful tools for exploring the landscape of RNA folding. Chemical probing alone may not necessarily define a single unique structure. The helix filters and optional multibranch loops are global constraints on RNA structure that are an especially useful tool for generating models of encapsidated viral RNA for which cryoelectron microscopy or crystallography data may be available. The computations generate a combinatorially complete set of structures near a free energy minimum and thus provide data on the density and diversity of structures near the bottom of a folding funnel for an RNA sequence. The conformational landscapes for some RNA sequences may resemble a low, wide basin rather than a steep funnel that converges to a single structure.

Microarrays for identifying binding sites and probing structure of RNAs

Oligonucleotide microarrays are widely used in various biological studies. In this review, application of oligonucleotide microarrays for identifying binding sites and probing structure of RNAs is described. Deep sequencing allows fast determination of DNA and RNA sequence. High-throughput methods for determination of secondary structures of RNAs have also been developed. Those methods, however, do not reveal binding sites for oligonucleotides. In contrast, microarrays directly determine binding sites while also providing structural insights. Microarray mapping can be used over a wide range of experimental conditions, including temperature, pH, various cations at different concentrations and the presence of other molecules. Moreover, it is possible to make universal microarrays suitable for investigations of many different RNAs, and readout of results is rapid. Thus, microarrays are used to provide insight into oligonucleotide sequences potentially able to interfere with biological function. Better understanding of structure-function relationships of RNA can be facilitated by using microarrays to find RNA regions capable to bind oligonucleotides. That information is extremely important to design optimal sequences for antisense oligonucleotides and siRNA because both bind to single-stranded regions of target RNAs.

Probing viral genomic structure: alternative viewpoints and alternative structures for satellite tobacco mosaic virus RNA

Viral RNA structure prediction is a valuable tool for development of drugs against viral disease. This work discusses different approaches to predicting encapsidated viral RNA and highlights satellite tobacco mosaic virus (STMV) RNA as a model system with excellent crystallography data. Fundamentally important issues for debate include thermodynamic versus kinetic control of virus assembly and the possible consequences of quasi-species in the primary structure on RNA secondary structure prediction of a single structure or an ensemble of structures. Multiple computational tools and chemical reagents are now available for improved viral RNA structure prediction. Two different predicted structures for encapsidated STMV RNA result from differences in three main areas: a different approach and philosophy to studying encapsidated viral RNA, an emphasis on different RNA motifs, and technical differences in computational methods and chemical reagents. The experiments with traditional chemical probing and SHAPE reagents are compared in terms of chemistry, results, and interpretation for STMV RNA as well as other RNA protein assemblies, such as the 5'UTR of HIV and the ribosome. This discussion of the challenges of viral RNA structure prediction will lead to new experiments and improved future predictions for viral RNA.

To build a virus on a nucleic acid substrate

Many viruses package their genomes concomitant with assembly. Here, we show that this reaction can be described by three coefficients: association of capsid protein (CP) to nucleic acid (NA), KNA; CP-CP interaction, ω; and α, proportional to the work required to package NA. The value of α can vary as NA is packaged. A phase diagram of average lnα versus lnω identifies conditions where assembly is likely to fail or succeed. NA morphology can favor (lnα > 0) or impede (lnα < 0) assembly. As lnω becomes larger, capsids become more stable and assembly becomes more cooperative. Where (lnα + lnω) < 0, the CP is unable to contain the NA, so that assembly results in aberrant particles. This phase diagram is consistent with quantitative studies of cowpea chlorotic mottle virus, hepatitis B virus, and simian virus 40 assembling on ssRNA and dsDNA substrates. Thus, the formalism we develop is suitable for describing and predicting behavior of experimental studies of CP assembly on NA.

Secondary structure of a conserved domain in the intron of influenza A NS1 mRNA

Influenza A virus is a segmented single-stranded (−)RNA virus that causes substantial annual morbidity and mortality. The transcriptome of influenza A is predicted to have extensive RNA secondary structure. The smallest genome segment, segment 8, encodes two proteins, NS1 and NEP, via alternative splicing. A conserved RNA domain in the intron of segment 8 may be important for regulating production of NS1. Two different multi-branch loop structures have been proposed for this region. A combination of in vitro chemical mapping and isoenergetic microarray techniques demonstrate that the consensus sequence for this region folds into a hairpin conformation. These results provide an alternative folding for this region and a foundation for designing experiments to probe its functional role in the influenza life cycle.

Packaging signals in two single-stranded RNA viruses imply a conserved assembly mechanism and geometry of the packaged genome

The current paradigm for assembly of single-stranded RNA viruses is based on a mechanism involving non-sequence-specific packaging of genomic RNA driven by electrostatic interactions. Recent experiments, however, provide compelling evidence for sequence specificity in this process both in vitro and in vivo. The existence of multiple RNA packaging signals (PSs) within viral genomes has been proposed, which facilitates assembly by binding coat proteins in such a way that they promote the protein-protein contacts needed to build the capsid. The binding energy from these interactions enables the confinement or compaction of the genomic RNAs. Identifying the nature of such PSs is crucial for a full understanding of assembly, which is an as yet untapped potential drug target for this important class of pathogens. Here, for two related bacterial viruses, we determine the sequences and locations of their PSs using Hamiltonian paths, a concept from graph theory, in combination with bioinformatics and structural studies. Their PSs have a common secondary structure motif but distinct consensus sequences and positions within the respective genomes. Despite these differences, the distributions of PSs in both viruses imply defined conformations for the packaged RNA genomes in contact with the protein shell in the capsid, consistent with a recent asymmetric structure determination of the MS2 virion. The PS distributions identified moreover imply a preferred, evolutionarily conserved assembly pathway with respect to the RNA sequence with potentially profound implications for other single-stranded RNA viruses known to have RNA PSs, including many animal and human pathogens.

The icosahedral RNA virus as a grotto: organizing the genome into stalagmites and stalactites

There are two important problems in the assembly of small, icosahedral RNA viruses. First, how does the capsid protein select the viral RNA for packaging, when there are so many other candidate RNA molecules available? Second, what is the mechanism of assembly? With regard to the first question, there are a number of cases where a particular RNA sequence or structure--often one or more stem-loops--either promotes assembly or is required for assembly, but there are others where specific packaging signals are apparently not required. With regard to the assembly pathway, in those cases where stem-loops are involved, the first step is generally believed to be binding of the capsid proteins to these "fingers" of the RNA secondary structure. In the mature virus, the core of the RNA would then occupy the center of the viral particle, and the stem-loops would reach outward, towards the capsid, like stalagmites reaching up from the floor of a grotto towards the ceiling. Those viruses whose assembly does not depend on protein binding to stem-loops could have a different structure, with the core of the RNA lying just under the capsid, and the fingers reaching down into the interior of the virus, like stalactites. We review the literature on these alternative structures, focusing on RNA selectivity and the assembly mechanism, and we propose experiments aimed at determining, in a given virus, which of the two structures actually occurs.

Long-range architecture in a viral RNA genome

We have developed a model for the secondary structure of the 1058-nucleotide plus-strand RNA genome of the icosahedral satellite tobacco mosaic virus (STMV) using nucleotide-resolution SHAPE chemical probing of the viral RNA isolated from virions and within the virion, perturbation of interactions distant in the primary sequence, and atomic force microscopy. These data are consistent with long-range base pairing interactions and a three-domain genome architecture. The compact domains of the STMV RNA have dimensions of 10-45 nm. Each of the three domains corresponds to a specific functional component of the virus: The central domain corresponds to the coding sequence of the single (capsid) protein encoded by the virus, whereas the 5' and 3' untranslated domains span signals essential for translation and replication, respectively. This three-domain architecture is compatible with interactions between the capsid protein and short RNA helices previously visualized by crystallography. STMV is among the simplest of the icosahedral viruses but, nonetheless, has an RNA genome with a complex higher-order structure that likely reflects high information content and an evolutionary relationship between RNA domain structure and essential replicative functions.

Packaging signals in single-stranded RNA viruses: nature's alternative to a purely electrostatic assembly mechanism

The formation of a protective protein container is an essential step in the life-cycle of most viruses. In the case of single-stranded (ss)RNA viruses, this step occurs in parallel with genome packaging in a co-assembly process. Previously, it had been thought that this process can be explained entirely by electrostatics. Inspired by recent single-molecule fluorescence experiments that recapitulate the RNA packaging specificity seen in vivo for two model viruses, we present an alternative theory, which recognizes the important cooperative roles played by RNA–coat protein interactions, at sites we have termed packaging signals. The hypothesis is that multiple copies of packaging signals, repeated according to capsid symmetry, aid formation of the required capsid protein conformers at defined positions, resulting in significantly enhanced assembly efficiency. The precise mechanistic roles of packaging signal interactions may vary between viruses, as we have demonstrated for MS2 and STNV. We quantify the impact of packaging signals on capsid assembly efficiency using a dodecahedral model system, showing that heterogeneous affinity distributions of packaging signals for capsid protein out-compete those of homogeneous affinities. These insights pave the way to a new anti-viral therapy, reducing capsid assembly efficiency by targeting of the vital roles of the packaging signals, and opens up new avenues for the efficient construction of protein nanocontainers in bionanotechnology.

Visualizing RNA base-pairing probabilities with RNAbow diagrams

There are many effective ways to represent a minimum free energy RNA secondary structure that make it easy to locate its helices and loops. It is a greater challenge to visualize the thermal average probabilities of all folds in a partition function sum; dot plot representations are often puzzling. Therefore, we introduce the RNAbows visualization tool for RNA base pair probabilities. RNAbows represent base pair probabilities with line thickness and shading, yielding intuitive diagrams. RNAbows aid in disentangling incompatible structures, allow comparisons between clusters of folds, highlight differences between wild-type and mutant folds, and are also rather beautiful.

In vitro secondary structure of the genomic RNA of satellite tobacco mosaic virus

Satellite tobacco mosaic virus (STMV) is a T = 1 icosahedral virus with a single-stranded RNA genome. It is widely accepted that the RNA genome plays an important structural role during assembly of the STMV virion. While the encapsidated form of the RNA has been extensively studied, less is known about the structure of the free RNA, aside from a purported tRNA-like structure at the 3' end. Here we use selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) analysis to examine the secondary structure of in vitro transcribed STMV RNA. The predicted secondary structure is unusual in the sense that it is highly extended, which could be significant for protecting the RNA from degradation. The SHAPE data are also consistent with the previously predicted tRNA-like fold at the 3' end of the molecule, which is also known to hinder degradation. Our data are not consistent with the secondary structure proposed for the encapsidated RNA by Schroeder et al., suggesting that, if the Schroeder structure is correct, either the RNA is packaged as it emerges from the replication complex, or the RNA undergoes extensive refolding upon encapsidation. We also consider the alternative, i.e., that the structure of the encapsidated STMV RNA might be the same as the in vitro structure presented here, and we examine how this structure might be organized in the virus. This possibility is not rigorously ruled out by the available data, so it remains open to examination by experiment.

Sequence-specific, RNA-protein interactions overcome electrostatic barriers preventing assembly of satellite tobacco necrosis virus coat protein

We have examined the roles of RNA–coat protein (CP) interactions in the assembly of satellite tobacco necrosis virus (STNV). The viral genomic RNA encodes only the CP, which comprises a β-barrel domain connected to a positively charged N-terminal extension. In the previous crystal structures of this system, the first 11 residues of the protein are disordered. Using variants of an RNA aptamer sequence isolated against the CP, B3, we have studied the sequence specificity of RNA-induced assembly. B3 consists of a stem–loop presenting the tetra-loop sequence ACAA. There is a clear preference for RNAs encompassing this loop sequence, as measured by the yield of T = 1 capsids, which is indifferent to sequences within the stem. The B3-containing virus-like particle has been crystallised and its structure was determined to 2.3 Å. A lower-resolution map encompassing density for the RNA has also been calculated. The presence of B3 results in increased ordering of the N-terminal helices located at the particle 3-fold axes, which extend by roughly one and a half turns to encompass residues 8–11, including R8 and K9. Under assembly conditions, STNV CP in the absence of RNA is monomeric and does not self-assemble. These facts suggest that a plausible model for assembly initiation is the specific RNA-induced stabilisation of a trimeric capsomere. The basic nature of the helical extension suggests that electrostatic repulsion between CPs prevents assembly in the absence of RNA and that this barrier is overcome by correct placement of appropriately orientated helical RNA stems. Such a mechanism would be consistent with the data shown here for assembly with longer RNA fragments, including an STNV genome. The results are discussed in light of a first stage of assembly involving compaction of the genomic RNA driven by multiple RNA packaging signal–CP interactions.

Crumple: a method for complete enumeration of all possible pseudoknot-free RNA secondary structures

The diverse landscape of RNA conformational space includes many canyons and crevices that are distant from the lowest minimum free energy valley and remain unexplored by traditional RNA structure prediction methods. A complete description of the entire RNA folding landscape can facilitate identification of biologically important conformations. The Crumple algorithm rapidly enumerates all possible non-pseudoknotted structures for an RNA sequence without consideration of thermodynamics while filtering the output with experimental data. The Crumple algorithm provides an alternative approach to traditional free energy minimization programs for RNA secondary structure prediction. A complete computation of all non-pseudoknotted secondary structures can reveal structures that would not be predicted by methods that sample the RNA folding landscape based on thermodynamic predictions. The free energy minimization approach is often successful but is limited by not considering RNA tertiary and protein interactions and the possibility that kinetics rather than thermodynamics determines the functional RNA fold. Efficient parallel computing and filters based on experimental data make practical the complete enumeration of all non-pseudoknotted structures. Efficient parallel computing for Crumple is implemented in a ring graph approach. Filters for experimental data include constraints from chemical probing of solvent accessibility, enzymatic cleavage of paired or unpaired nucleotides, phylogenetic covariation, and the minimum number and lengths of helices determined from crystallography or cryo-electron microscopy. The minimum number and length of helices has a significant effect on reducing conformational space. Pairing constraints reduce conformational space more than single nucleotide constraints. Examples with Alfalfa Mosaic Virus RNA and Trypanosome brucei guide RNA demonstrate the importance of evaluating all possible structures when pseduoknots, RNA-protein interactions, and metastable structures are important for biological function. Crumple software is freely available at http://adenosine.chem.ou.edu/software.html.

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.

A model for the structure of satellite tobacco mosaic virus

Satellite tobacco mosaic virus (STMV) is an icosahedral T=1 single-stranded RNA virus with a genome containing 1058 nucleotides. X-ray crystallography revealed a structure containing 30 double-helical RNA segments, with each helix having nine base pairs and an unpaired nucleotide at the 3' end of each strand. Based on this structure, Larson and McPherson proposed a model of 30 hairpin-loop elements occupying the edges of the icosahedron and connected by single-stranded regions. More recently, Schroeder et al. have combined the results of chemical probing with a novel helix searching algorithm to propose a specific secondary structure for the STMV genome, compatible with the Larson-McPherson model. Here we report an all-atom model of STMV, using the complete protein and RNA sequences and the Schroeder RNA secondary structure. As far as we know, this is the first all-atom model for the complete structure of any virus (100% of the atoms) using the natural genomic sequence.

Visualizing large RNA molecules in solution

Single-stranded RNAs (ssRNAs) longer than a few hundred nucleotides do not have a unique structure in solution. Their equilibrium properties therefore reflect the average of an ensemble of structures. We use cryo-electron microscopy to image projections of individual long ssRNA molecules and characterize the anisotropy of their ensembles in solution. A flattened prolate volume is found to best represent the shapes of these ensembles. The measured sizes and anisotropies are in good agreement with complementary determinations using small-angle X-ray scattering and coarse-grained molecular dynamics simulations. A long viral ssRNA is compared with shorter noncoding transcripts to demonstrate that prolate geometry and flatness are generic properties independent of sequence length and origin. The anisotropy persists under physiological as well as low-ionic-strength conditions, revealing a direct correlation between secondary structure asymmetry and 3D shape and size. We discuss the physical origin of the generic anisotropy and its biological implications.