Comprehensive survey of conserved RNA secondary structures in full-genome alignment of Hepatitis C virus.

Hepatitis C virus (HCV) is a plus-stranded RNA virus that often chronically infects liver hepatocytes and causes liver cirrhosis and cancer. These viruses replicate their genomes employing error-prone replicases. Thereby, they routinely generate a large 'cloud' of RNA genomes (quasispecies) which-by trial and error-comprehensively explore the sequence space available for functional RNA genomes that maintain the ability for efficient replication and immune escape. In this context, it is important to identify which RNA secondary structures in the sequence space of the HCV genome are conserved, likely due to functional requirements. Here, we provide the first genome-wide multiple sequence alignment (MSA) with the prediction of RNA secondary structures throughout all representative full-length HCV genomes. We selected 57 representative genomes by clustering all complete HCV genomes from the BV-BRC database based on k-mer distributions and dimension reduction and adding RefSeq sequences. We include annotations of previously recognized features for easy comparison to other studies. Our results indicate that mainly the core coding region, the C-terminal NS5A region, and the NS5B region contain secondary structure elements that are conserved beyond coding sequence requirements, indicating functionality on the RNA level. In contrast, the genome regions in between contain less highly conserved structures. The results provide a complete description of all conserved RNA secondary structures and make clear that functionally important RNA secondary structures are present in certain HCV genome regions but are largely absent from other regions. Full-genome alignments of all branches of Hepacivirus C are provided in the supplement.

Exploring the parameter space of complex self-assembly through virus capsid models.

We use discrete event stochastic simulations to characterize the parameter space of a model of icosahedral viral capsid assembly as functions of monomer-monomer binding rates. The simulations reveal a parameter space characterized by three major assembly mechanisms, a standard nucleation-limited monomer-accretion pathway and two distinct hierarchical assembly pathways, as well as unproductive regions characterized by kinetically trapped species. Much of the productive parameter space also consists of border regions between these domains where hybrid pathways are likely to operate. A simpler octamer system studied for comparison reveals three analogous pathways, but is characterized by much lesser sensitivity to parameter variations in contrast to the sharp changes visible in the icosahedral model. The model suggests that modest changes in assembly conditions, consistent with expected differences between in vitro and in vivo assembly environments, could produce substantial shifts in assembly pathways. These results suggest that we must be cautious in drawing conclusions about in vivo capsid self-assembly dynamics from theoretical or in vitro models, as the nature of the basic assembly mechanisms accessible to a system can substantially differ between simple and complex model systems, between theoretical models and simulation results, and between in vitro and in vivo assembly conditions.