Biochemical and genetic studies of the initiation of human rhinovirus 2 RNA replication: identification of a cis-replicating element in the coding sequence of 2A(pro)

AAACH is a conserved motif in a cis-acting replication element that is artificially inserted into Senecavirus A genome

Cis-acting replication element (cre) is required for generating a diuridylylated VPg that acts as a protein primer to initiate the synthesis of picornaviral genome or antigenome. The cre is a stem-loop structure, dependent of different picornaviruses, located in different genomic regions. The AAACA motif is highly conserved in the apical loop of cre among several picornaviral members, and plays a key role in synthesizing a diuridylylated VPg. We previously demonstrated that senecavirus A (SVA) also possesses an AAACA-containing cre in its genome. Its natural cre (Nc), if functionally inactivated through site-directed mutagenesis (SDM), would confer a lethal impact on virus recovery, whereas an artificial cre (Ac) is able to compensate for the Nc-caused functional inactivation, leading to successful rescue of a viable SVA. In this study, we constructed a set of SVA cDNA clones. Each of them contained one functionally inactivated Nc, and an extra SDM-modified Ac. Every cDNA clone had a unique SDM-modified Ac. The test of virus recovery showed that only two SVAs were rescued from their individual cDNA clones. They were AAACU- and AAACC-containing Ac genotypes. Both viruses were serially passaged in vitro for analyzing their viral characteristics. The results showed that both AAACU and AAACC genotypes were genetically stable during twenty passages, implying when the Nc was functionally inactivated, SVA could still use an AAACH-containing Ac to complete its own replication cycle.

Structural and nonstructural proteins of Senecavirus A: Recent research advances, and lessons learned from those of other picornaviruses

Senecavirus A (SVA) is an emerging virus, causing vesicular disease in swine. SVA is a single-stranded, positive-sense RNA virus, which is the only member of the genus Senecavirus in the family Picornaviridae. SVA genome encodes 12 proteins: L, VP4, VP2, VP3, VP1, 2A, 2B, 2C, 3A, 3B, 3C and 3D. The VP1 to VP4 are structural proteins, and the others are nonstructural proteins. The replication of SVA in host cells is a complex process coordinated by an elaborate interplay between the structural and nonstructural proteins. Structural proteins are primarily involved in the invasion and assembly of virions. Nonstructural proteins modulate viral RNA translation and replication, and also take part in antagonizing the antiviral host response and in disrupting some cellular processes to allow virus replication. Here, we systematically reviewed the molecular functions of SVA structural and nonstructural proteins by reference to literatures of SVA itself and other picornaviruses.

Identification of cis-acting replication element in VP2-encoding region of Senecavirus A genome

Picornavirus possesses one positive-sense, single-stranded RNA genome, in which a cis-acting replication element (cre) is located. The cre is a stem-loop structure that harbors a conserved AAACA motif within its loop region. This motif functions as a template for adding two U residues to the viral VPg, therefore generating a VPg-pUpU that is required for viral RNA synthesis. Senecavirus A (SVA) is an emerging picornavirus. Its cre has not been identified as yet. In the present study, one putative cre containing a typical AAACA motif was computationally predicted to exist within the VP2-encoding sequence of SVA. To test the role of this putative cre, 22 SVA cDNA clones with different point mutations in their cre-formed sequences were constructed in an attempt to rescue replication-competent SVAs. A total of 11 viruses were rescued from their individual cDNA clones, implying that some mutated cres exerted lethal impacts on SVA replication. To eliminate these impacts, an intact cre was artificially inserted into those SVA cDNA clones without ability of recovering virus. The artificial cre was proven to be able of compensating for some, but not all, defects caused by mutated cres, leading to successful recovery of SVAs. These results indicated that the putative cre of SVA was functionally similar to those of other picornaviruses, perhaps involved in the uridylylation of VPg.

A single mutation in the cis-acting replication element identified within the EV-A71 2C-coding region causes defects in virus production in cell culture

ABSTRACTEnterovirus A71 (EV-A71) can cause hand, foot and mouth disease with neurological and systemic complications, most frequently affecting children and infants. We describe a cis-acting replication element (cre) with a conserved stem-loop structure within the EV-A71 2C-coding region. By site-directed mutagenesis and reverse genetics using the EV-A71 full-length genome and the EV-A71 replicon containing the firefly luciferase reporter gene in place of the P1 region, the stem-loop structure and the AAACA in the loop of the cre were confirmed to be required for the EV-A71 replication phenotype. EV-A71 genomes containing a mutation at the first or third A residue of AAACA could not be recovered. Insertion of a wild-type cre from EV-A71 or poliovirus in the 5'UTR led to successful recovery of the replication of nonviable mutants. Furthermore, the cre mutants showed lower binding capacity with the host cellular factor IGF2BP2, knockdown of which resulted in a significant decrease in EV-A71 production. All the available evidence shows the location independence but functional importance of the interaction of the cre with the cellular host for efficient production of EV-A71, contributing to the growing body of knowledge regarding picornavirus cres.

Protein Nucleotidylylation in +ssRNA Viruses

Nucleotidylylation is a post-transcriptional modification important for replication in the picornavirus supergroup of RNA viruses, including members of the Caliciviridae, Coronaviridae, Picornaviridae and Potyviridae virus families. This modification occurs when the RNA-dependent RNA polymerase (RdRp) attaches one or more nucleotides to a target protein through a nucleotidyl-transferase reaction. The most characterized nucleotidylylation target is VPg (viral protein genome-linked), a protein linked to the 5' end of the genome in Caliciviridae, Picornaviridae and Potyviridae. The nucleotidylylation of VPg by RdRp is a critical step for the VPg protein to act as a primer for genome replication and, in Caliciviridae and Potyviridae, for the initiation of translation. In contrast, Coronaviridae do not express a VPg protein, but the nucleotidylylation of proteins involved in replication initiation is critical for genome replication. Furthermore, the RdRp proteins of the viruses that perform nucleotidylylation are themselves nucleotidylylated, and in the case of coronavirus, this has been shown to be essential for viral replication. This review focuses on nucleotidylylation within the picornavirus supergroup of viruses, including the proteins that are modified, what is known about the nucleotidylylation process and the roles that these modifications have in the viral life cycle.

Dengue virus 2 capsid protein chaperones the strand displacement of 5'-3' cyclization sequences

By virtue of its chaperone activity, the capsid protein of dengue virus strain 2 (DENV2C) promotes nucleic acid structural rearrangements. However, the role of DENV2C during the interaction of RNA elements involved in stabilizing the 5′-3′ panhandle structure of DENV RNA is still unclear. Therefore, we determined how DENV2C affects structural functionality of the capsid-coding region hairpin element (cHP) during annealing and strand displacement of the 9-nt cyclization sequence (5CS) and its complementary 3CS. cHP has two distinct functions: a role in translation start codon selection and a role in RNA synthesis. Our results showed that cHP impedes annealing between 5CS and 3CS. Although DENV2C does not modulate structural functionality of cHP, it accelerates annealing and specifically promotes strand displacement of 3CS during 5′-3′ panhandle formation. Furthermore, DENV2C exerts its chaperone activity by favouring one of the active conformations of cHP. Based on our results, we propose mechanisms for annealing and strand displacement involving cHP. Thus, our results provide mechanistic insights into how DENV2C regulates RNA synthesis by modulating essential RNA elements in the capsid-coding region, that in turn allow for DENV replication.

Effects of TDP2/VPg Unlinkase Activity on Picornavirus Infections Downstream of Virus Translation

In this study, we characterized the role of host cell protein tyrosyl-DNA phosphodiesterase 2 (TDP2) activity, also known as VPg unlinkase, in picornavirus infections in a human cell model of infection. TDP2/VPg unlinkase is used by picornaviruses to remove the small polypeptide, VPg (Virus Protein genome-linked, the primer for viral RNA synthesis), from virus genomic RNA. We utilized a CRISPR/Cas-9-generated TDP2 knock out (KO) human retinal pigment epithelial-1 (hRPE-1) cell line, in addition to the wild type (WT) counterpart for our studies. We determined that in the absence of TDP2, virus growth kinetics for two enteroviruses (poliovirus and coxsackievirus B3) were delayed by about 2 h. Virus titers were reduced by ~2 log10 units for poliovirus and 0.5 log10 units for coxsackievirus at 4 hours post-infection (hpi), and by ~1 log10 unit at 6 hpi for poliovirus. However, virus titers were nearly indistinguishable from those of control cells by the end of the infectious cycle. We determined that this was not the result of an alternative source of VPg unlinkase activity being activated in the absence of TPD2 at late times of infection. Viral protein production in TDP2 KO cells was also substantially reduced at 4 hpi for poliovirus infection, consistent with the observed growth kinetics delay, but reached normal levels by 6 hpi. Interestingly, this result differs somewhat from what has been reported previously for the TDP2 KO mouse cell model, suggesting that either cell type or species-specific differences might be playing a role in the observed phenotype. We also determined that catalytically inactive TDP2 does not rescue the growth defect, confirming that TDP2 5' phosphodiesterase activity is required for efficient virus replication. Importantly, we show for the first time that polysomes can assemble efficiently on VPg-linked RNA after the initial round of translation in a cell culture model, but both positive and negative strand RNA production is impaired in the absence of TDP2 at mid-times of infection, indicating that the presence of VPg on the viral RNA affects a step in the replication cycle downstream of translation (e.g., RNA synthesis). In agreement with this conclusion, we found that double-stranded RNA production (a marker of viral RNA synthesis) is delayed in TDP2 KO RPE-1 cells. Moreover, we show that premature encapsidation of nascent, VPg-linked RNA is not responsible for the observed virus growth defect. Our studies provide the first lines of evidence to suggest that either negative- or positive-strand RNA synthesis (or both) is a likely candidate for the step that requires the removal of VPg from the RNA for an enterovirus infection to proceed efficiently.

Modification of betanodavirus virulence by substitutions in the 3' terminal region of RNA2

Betanodaviruses have bi-segmented positive-sense RNA genomes, consisting of RNAs 1 and 2. For some members of the related genus alphanodavirus, the 3' terminal 50 nucleotides (nt) of RNA2, including a predicted stem-loop structure (3'SL), are essential for replication. We investigate the possible existence and role of a similar structure in a reassortant betanodavirus strain (RGNNV/SJNNV). In this study, we developed three recombinant strains containing nucleotide changes at positions 1408 and 1412. Predictive models showed stem-loop structures involving nt 1398-1421 of the natural reassortant whereas this structure is modified in the recombinant viruses harbouring point mutations r1408 and r1408-1412, but not in r1412. Results obtained from infectivity assays showed differences between the reference strains and the mutants in both RNA1 and RNA2 synthesis. Moreover, an imbalance between the synthesis of both segments was demonstrated, mainly with the double mutant. All these results suggest an interaction between RNA1 and the 3' non-coding regions (3'NCR) of RNA2. In addition, the significant attenuation of the virulence for Senegalese sole and the delayed replication of r1408-1412 in brain tissues may point to an interaction of RNA2 with host cellular proteins.

Biological function of Foot-and-mouth disease virus non-structural proteins and non-coding elements

Foot-and-mouth disease virus (FMDV) represses host translation machinery, blocks protein secretion, and cleaves cellular proteins associated with signal transduction and the innate immune response to infection. Non-structural proteins (NSPs) and non-coding elements (NCEs) of FMDV play a critical role in these biological processes. The FMDV virion consists of capsid and nucleic acid. The virus genome is a positive single stranded RNA and encodes a single long open reading frame (ORF) flanked by a long structured 5ʹ-untranslated region (5ʹ-UTR) and a short 3ʹ-UTR. The ORF is translated into a polypeptide chain and processed into four structural proteins (VP1, VP2, VP3, and VP4), 10 NSPs (L^pro, 2A, 2B, 2C, 3A, 3B_1–3, 3C^pro, and 3D^pol), and some cleavage intermediates. In the past decade, an increasing number of studies have begun to focus on the molecular pathogenesis of FMDV NSPs and NCEs. This review collected recent research progress on the biological functions of these NSPs and NCEs on the replication and host cellular regulation of FMDV to understand the molecular mechanism of host–FMDV interactions and provide perspectives for antiviral strategy and development of novel vaccines.

Mutational Disruption of cis-Acting Replication Element 2C in Coxsackievirus B3 Leads to 5'-Terminal Genomic Deletions

Following natural human or experimental murine infections and in cell culture, coxsackievirus B (CVB) RNA can persist for weeks in the absence of a cytopathic effect, yet viral RNA remains detectable. Our earlier studies demonstrated that this persistence produced viral RNA with up to 49 nucleotide deletions at the genomic 5' terminus which partially degraded the cloverleaf (or domain I), an RNA structure required for efficient viral replication. A cis-acting replication element (CRE) in the 2C protein-coding region [CRE(2C)] templates the addition of two uridine residues to the virus genome-encoded RNA replication primer VPg prior to positive-strand synthesis. Because our previous work also demonstrated that the genomes of CVB with a 5'-terminal deletion (CVB-TD) have VPg covalently linked, even though they rarely terminate in the canonical UU donated by CRE(2C)-mediated uridylylation of VPg, we hypothesized that a functional (uridylylating) CRE(2C) would be unnecessary for CVB-TD replication. Using the same 16 mutations in the CVB3 CRE(2C) structure that were considered lethal for this virus by others, we demonstrate here both in infected cell cultures and in mice that wild-type (wt) and CVB3-TD strains carrying these mutations with a nonuridylylating CRE(2C) are viable. While the wt genome with the mutated CRE(2C) displays suppressed replication levels similar to those observed in a CVB3-TD strain, mutation of the CRE(2C) function in a CVB3-TD strain does not further decrease replication. Finally, we show that replication of the parental CVB3 strain containing the mutated CRE(2C) drives the de novo generation of genomic deletions at the 5' terminus.

IMPORTANCE

In this report, we demonstrate that while CVB can replicate without a uridylylating CRE(2C), the replication rate suffers significantly. Further, deletions at the 5' terminus of the genome are generated in this virus population, with this virus population supplanting the wild-type population. This demonstrates that VPg can prime without being specifically uridylylated and that this priming is error prone, resulting in the loss of sequence information from the 5' terminus. These findings have significance when considering the replication of human enteroviruses, and we believe that these data are unattainable in a cell-free system due to the poor replication of these CRE-deficient viruses.

Replication and Inhibitors of Enteroviruses and Parechoviruses

The Enterovirus (EV) and Parechovirus genera of the picornavirus family include many important human pathogens, including poliovirus, rhinovirus, EV-A71, EV-D68, and human parechoviruses (HPeV). They cause a wide variety of diseases, ranging from a simple common cold to life-threatening diseases such as encephalitis and myocarditis. At the moment, no antiviral therapy is available against these viruses and it is not feasible to develop vaccines against all EVs and HPeVs due to the great number of serotypes. Therefore, a lot of effort is being invested in the development of antiviral drugs. Both viral proteins and host proteins essential for virus replication can be used as targets for virus inhibitors. As such, a good understanding of the complex process of virus replication is pivotal in the design of antiviral strategies goes hand in hand with a good understanding of the complex process of virus replication. In this review, we will give an overview of the current state of knowledge of EV and HPeV replication and how this can be inhibited by small-molecule inhibitors.

Chimeric rhinoviruses obtained via genetic engineering or artificially induced recombination are viable only if the polyprotein coding sequence derives from the same species

Recombination is a widespread phenomenon that ensures both the stability and variation of RNA viruses. This phenomenon occurs with different frequencies within species of the Enterovirus genus. Intraspecies recombination is described frequently among non-rhinovirus enteroviruses but appears to be sporadic in rhinoviruses. Interspecies recombination is even rarer for rhinoviruses and mostly is related to ancient events which contributed to the speciation of these viruses. We reported that artificially engineered 5' untranslated region (UTR) interspecies rhinovirus/rhinovirus or rhinovirus/non-rhinovirus enterovirus recombinants are fully viable. Using a similar approach, we demonstrated in this study that exchanges of the P1-2A polyprotein region between members of the same rhinovirus species, but not between members of different species, give rise to competent chimeras. To further assess the rhinovirus intra- and interspecies recombination potential, we used artificially induced recombination by cotransfection of 5'-end-deleted and 3'-end-deleted and replication-deficient genomes. In this system, intraspecies recombination also resulted in viable viruses with high frequency, whereas no interspecies rhinovirus recombinants could be recovered. Mapping intraspecies recombination sites within the polyprotein highlighted recombinant hotspots in nonstructural genes and at gene boundaries. Notably, all recombinants occurring at gene junctions presented in-frame sequence duplications, whereas most intragenic recombinants were homologous. Taken together, our results suggest that only intraspecies recombination gives rise to viable rhinovirus chimeras in the polyprotein coding region and that recombination hotspots map to nonstructural genes with in-frame duplications at gene boundaries. These data provide new insights regarding the mechanism and limitations of rhinovirus recombination.

IMPORTANCE

Recombination represents a means to ensure both the stability and the variation of RNA viruses. While intraspecies recombination is described frequently among non-rhinovirus enteroviruses, it seems to occur more rarely in rhinoviruses. Interspecies recombination is even rarer in this virus group and is mostly related to ancient events, which contributed to its speciation. We used engineered chimeric genomes and artificially induced RNA recombination to study experimentally the recombination potential of rhinoviruses and analyze recombination sites. Our results suggest that only intraspecies recombination gives rise to viable chimeras in the polyprotein coding region. Furthermore, characterization of intraspecies chimeras provides new insight into putative recombination hotspots within the polyprotein. In summary, we applied two powerful and complementary experimental approaches to improve current knowledge on rhinovirus recombination.

Initiation of protein-primed picornavirus RNA synthesis

Plus strand RNA viruses use different mechanisms to initiate the synthesis of their RNA chains. The Picornaviridae family constitutes a large group of plus strand RNA viruses that possess a small terminal protein (VPg) covalently linked to the 5'-end of their genomes. The RNA polymerases of these viruses use VPg as primer for both minus and plus strand RNA synthesis. In the first step of the initiation reaction the RNA polymerase links a UMP to the hydroxyl group of a tyrosine in VPg using as template a cis-replicating element (cre) positioned in different regions of the viral genome. In this review we will summarize what is known about the initiation reaction of protein-primed RNA synthesis by the RNA polymerases of the Picornaviridae. As an example we will use the RNA polymerase of poliovirus, the prototype of Picornaviridae. We will also discuss models of how these nucleotidylylated protein primers might be used, together with viral and cellular replication proteins and other cis-replicating RNA elements, during minus and plus strand RNA synthesis.

Full-length genomic analysis of Korean porcine Sapelovirus strains

Porcine sapelovirus (PSV), a species of the genus Sapelovirus within the family Picornaviridae, is associated with diarrhea, pneumonia, severe neurological disorders, and reproductive failure in pigs. However, the structural features of the complete PSV genome remain largely unknown. To analyze the structural features of PSV genomes, the full-length nucleotide sequences of three Korean PSV strains were determined and analyzed using bioinformatic techniques in comparison with other known PSV strains. The Korean PSV genomes ranged from 7,542 to 7,566 nucleotides excluding the 3′ poly(A) tail, and showed the typical picornavirus genome organization; 5′untranslated region (UTR)-L-VP4-VP2-VP3-VP1-2A-2B-2C-3A-3B-3C-3D-3′UTR. Three distinct cis-active RNA elements, the internal ribosome entry site (IRES) in the 5′UTR, a cis-replication element (CRE) in the 2C coding region and 3′UTR were identified and their structures were predicted. Interestingly, the structural features of the CRE and 3′UTR were different between PSV strains. The availability of these first complete genome sequences for PSV strains will facilitate future investigations of the molecular pathogenesis and evolutionary characteristics of PSV.

Amiloride inhibits the initiation of Coxsackievirus and poliovirus RNA replication by inhibiting VPg uridylylation

The mechanism of amiloride inhibition of Coxsackievirus B3 (CVB3) and poliovirus type 1 (PV1) RNA replication was investigated using membrane-associated RNA replication complexes. Amiloride was shown to inhibit viral RNA replication and VPgpUpU synthesis. However, the drug had no effect on polymerase elongation activity during either (-) strand or (+) strand synthesis. These findings indicated that amiloride inhibited the initiation of RNA synthesis by inhibiting VPg uridylylation. In addition, in silico binding studies showed that amiloride docks in the VPg binding site on the back of the viral RNA polymerase, 3D(pol). Since VPg binding at this site on PV1 3D(pol) was previously shown to be required for VPg uridylylation, our results suggest that amiloride inhibits VPg binding to 3D(pol). In summary, our findings are consistent with a model in which amiloride inhibits VPgpUpU synthesis and viral RNA replication by competing with VPg for binding to 3D(pol).

The influence of viral RNA secondary structure on interactions with innate host cell defences

RNA viruses infecting vertebrates differ fundamentally in their ability to establish persistent infections with markedly different patterns of transmission, disease mechanisms and evolutionary relationships with their hosts. Although interactions with host innate and adaptive responses are complex and persistence mechanisms likely multi-factorial, we previously observed associations between bioinformatically predicted RNA secondary formation in genomes of positive-stranded RNA viruses with their in vivo fitness and persistence. To analyse this interactions functionally, we transfected fibroblasts with non-replicating, non-translated RNA transcripts from RNA viral genomes with differing degrees of genome-scale ordered RNA structure (GORS). Single-stranded RNA transcripts induced interferon-β mediated though RIG-I and PKR activation, the latter associated with rapid induction of antiviral stress granules. A striking inverse correlation was observed between induction of both cellular responses with transcript RNA structure formation that was independent of both nucleotide composition and sequence length. The consistent inability of cells to recognize RNA transcripts possessing GORS extended to downstream differences from unstructured transcripts in expression of TNF-α, other interferon-stimulated genes and induction of apoptosis. This functional association provides novel insights into interactions between virus and host early after infection and provides evidence for a novel mechanism for evading intrinsic and innate immune responses.

Human rhinovirus VPg uridylylation AlphaScreen for high-throughput screening

As an obligate step for picornaviruses to replicate their genome, the small viral peptide VPg must first be specifically conjugated with uridine nucleotides at a conserved tyrosine hydroxyl group. The resulting VPg-pUpU serves as the primer for genome replication. The uridylylation reaction requires the coordinated activity of many components, including the viral polymerase, a conserved internal RNA stem loop structure, and additional viral proteins. Formation of this complex and the resulting conjugation reaction catalyzed by the polymerase, offers a number of biochemical targets for inhibition of an essential process in the viral life cycle. Therefore, an assay recapitulating uridylylation would provide multiple opportunities for discovering potential antiviral agents. Our goal was to identify inhibitors of human rhinovirus (HRV) VPg uridylylation, which might ultimately be useful to reduce or prevent HRV-induced lower airway immunologic inflammatory responses, a major cause of asthma and chronic obstructive pulmonary disease exacerbations. We have reconstituted the complex uridylylation reaction in an AlphaScreen suitable for high-throughput screening, in which a rabbit polyclonal antiserum specific for uridylylated VPg serves as a key reagent. Assay results were validated by quantitative mass spectrometric detection of uridylylation.

Non-template functions of viral RNA in picornavirus replication

The genomic RNA of poliovirus and closely related picornaviruses perform template and non-template functions during viral RNA replication. The non-template functions are mediated by cis-active RNA sequences that bind viral and cellular proteins to form RNP complexes. The RNP complexes mediate temporally dynamic, long-range interactions in the viral genome and ensure the specificity of replication. The 5' cloverleaf (5' CL)-RNP complex serves as a key cis-active element in all of the non-template functions of viral RNA. The 5'CL-RNP complex is proposed to interact with the cre-RNP complex during VPgpUpU synthesis, the 3'NTR-poly(A) RNP complex during negative-strand initiation and the 30 end negative-strand-RNP complex during positive-strand initiation. Co-ordinating these long-range interactions is important in regulating each step in the replication cycle.

The encephalomyocarditis virus

The encephalomyocarditis virus (EMCV) is a small non-enveloped single-strand RNA virus, the causative agent of not only myocarditis and encephalitis, but also neurological diseases, reproductive disorders and diabetes in many mammalian species. EMCV pathogenesis appears to be viral strain- and host-specific, and a better understanding of EMCV virulence factors is increasingly required. Indeed, EMCV is often used as a model for diabetes and viral myocarditis, and is also widely used in immunology as a double-stranded RNA stimulus in the study of Toll-like as well as cytosolic receptors. However, EMCV virulence and properties have often been neglected. Moreover, EMCV is able to infect humans albeit with a low morbidity. Progress on xenografts, such as pig heart transplantation in humans, has raised safety concerns that need to be explored. In this review we will highlight the biology of EMCV and all known and potential virulence factors.

Toward genetics-based virus taxonomy: comparative analysis of a genetics-based classification and the taxonomy of picornaviruses

Virus taxonomy has received little attention from the research community despite its broad relevance. In an accompanying paper (C. Lauber and A. E. Gorbalenya, J. Virol. 86:3890-3904, 2012), we have introduced a quantitative approach to hierarchically classify viruses of a family using pairwise evolutionary distances (PEDs) as a measure of genetic divergence. When applied to the six most conserved proteins of the Picornaviridae, it clustered 1,234 genome sequences in groups at three hierarchical levels (to which we refer as the "GENETIC classification"). In this study, we compare the GENETIC classification with the expert-based picornavirus taxonomy and outline differences in the underlying frameworks regarding the relation of virus groups and genetic diversity that represent, respectively, the structure and content of a classification. To facilitate the analysis, we introduce two novel diagrams. The first connects the genetic diversity of taxa to both the PED distribution and the phylogeny of picornaviruses. The second depicts a classification and the accommodated genetic diversity in a standardized manner. Generally, we found striking agreement between the two classifications on species and genus taxa. A few disagreements concern the species Human rhinovirus A and Human rhinovirus C and the genus Aphthovirus, which were split in the GENETIC classification. Furthermore, we propose a new supergenus level and universal, level-specific PED thresholds, not reached yet by many taxa. Since the species threshold is approached mostly by taxa with large sampling sizes and those infecting multiple hosts, it may represent an upper limit on divergence, beyond which homologous recombination in the six most conserved genes between two picornaviruses might not give viable progeny.

Replication enhancer elements within the open reading frame of tick-borne encephalitis virus and their evolution within the Flavivirus genus

We provide experimental evidence of a replication enhancer element (REE) within the capsid gene of tick-borne encephalitis virus (TBEV, genus Flavivirus). Thermodynamic and phylogenetic analyses predicted that the REE folds as a long stable stem–loop (designated SL6), conserved among all tick-borne flaviviruses (TBFV). Homologous sequences and potential base pairing were found in the corresponding regions of mosquito-borne flaviviruses, but not in more genetically distant flaviviruses. To investigate the role of SL6, nucleotide substitutions were introduced which changed a conserved hexanucleotide motif, the conformation of the terminal loop and the base-paired dsRNA stacking. Substitutions were made within a TBEV reverse genetic system and recovered mutants were compared for plaque morphology, single-step replication kinetics and cytopathic effect. The greatest phenotypic changes were observed in mutants with a destabilized stem. Point mutations in the conserved hexanucleotide motif of the terminal loop caused moderate virus attenuation. However, all mutants eventually reached the titre of wild-type virus late post-infection. Thus, although not essential for growth in tissue culture, the SL6 REE acts to up-regulate virus replication. We hypothesize that this modulatory role may be important for TBEV survival in nature, where the virus circulates by non-viraemic transmission between infected and non-infected ticks, during co-feeding on local rodents.

Synthesis of nucleotidylated poliovirus VPg proteins

Phosphitylation of the side chain hydroxyl function of Fmoc protected tyrosine with 5'-phosphoramidites of suitably protected cytidine, adenosine, and guanosine, followed by oxidation gave three novel nucleotidylated amino acid building blocks. After protective group manipulation, these building blocks were used in a solid phase peptide synthesis to afford the nucleotidylated poliovirus proteins VPgpC, VPgpA, and VPgpG.

Analysis of genetic diversity and sites of recombination in human rhinovirus species C

Human rhinoviruses (HRVs) are a highly prevalent and diverse group of respiratory viruses. Although HRV-A and HRV-B are traditionally detected by virus isolation, a series of unculturable HRV variants have recently been described and assigned as a new species (HRV-C) within the picornavirus Enterovirus genus. To investigate their genetic diversity and occurrence of recombination, we have performed comprehensive phylogenetic analysis of sequences from the 5' untranslated region (5' UTR), VP4/VP2, VP1, and 3Dpol regions amplified from 89 HRV-C-positive respiratory samples and available published sequences. Branching orders of VP4/VP2, VP1, and 3Dpol trees were identical, consistent with the absence of intraspecies recombination in the coding regions. However, numerous tree topology changes were apparent in the 5' UTR, where >60% of analyzed HRV-C variants showed recombination with species A sequences. Two recombination hot spots in stem-loop 5 and the polypyrimidine tract in the 5' UTR were mapped using the program GroupingScan. Available HRV-C sequences showed evidence for additional interspecies recombination with HRV-A in the 2A gene, with breakpoints mapping precisely to the boundaries of the C-terminal domain of the encoded proteinase. Pairwise distances between HRV-C variants in VP1 and VP4/VP2 regions fell into two separate distributions, resembling inter- and intraserotype distances of species A and B. These observations suggest that, without serological cross-neutralization data, HRV-C genetic groups may be equivalently classified into types using divergence thresholds derived from distance distributions. The extensive sequence data from multiple genome regions of HRV-C and analyses of recombination in the current study will assist future formulation of consensus criteria for HRV-C type assignment and identification.

Proposals for the classification of human rhinovirus species C into genotypically assigned types

Human rhinoviruses (HRVs) are common respiratory pathogens associated with mild upper respiratory tract infections, but also increasingly recognized in the aetiology of severe lower respiratory tract disease. Wider use of molecular diagnostics has led to a recent reappraisal of HRV genetic diversity, including the discovery of HRV species C (HRV-C), which is refractory to traditional virus isolation procedures. Although it is heterogeneous genetically, there has to date been no attempt to classify HRV-C into types analogous to the multiple serotypes identified for HRV-A and -B and among human enteroviruses. Direct investigation of cross-neutralization properties of HRV-C is precluded by the lack of methods for in vitro culture, but sequences from the capsid genes (VP1 and partial VP4/VP2) show evidence for marked phylogenetic clustering, suggesting the possibility of a genetically based system comparable to that used for the assignment of new enterovirus types. We propose a threshold of 13% divergence for VP1 nucleotide sequences for type assignment, a level that classifies the current dataset of 86 HRV-C VP1 sequences into a total of 33 types. We recognize, however, that most HRV-C sequence data have been collected in the VP4/VP2 region (currently 701 sequences between positions 615 and 1043). We propose a subsidiary classification of variants showing > 10% divergence in VP4/VP2, but lacking VP1 sequences, to 28 provisionally assigned types (subject to confirmation once VP1 sequences are determined). These proposals will assist in future epidemiological and clinical studies of HRV-C conducted by different groups worldwide, and provide the foundation for future exploration of type-associated differences in clinical presentations and biological properties.

Rhinovirus genome evolution during experimental human infection

Human rhinoviruses (HRVs) evolve rapidly due in part to their error-prone RNA polymerase. Knowledge of the diversity of HRV populations emerging during the course of a natural infection is essential and represents a basis for the design of future potential vaccines and antiviral drugs. To evaluate HRV evolution in humans, nasal wash samples were collected daily for five days from 15 immunocompetent volunteers experimentally infected with a reference stock of HRV-39. In parallel, HeLa-OH cells were inoculated to compare HRV evolution in vitro. Nasal wash in vivo assessed by real-time PCR showed a viral load that peaked at 48–72 h. Ultra-deep sequencing was used to compare the low-frequency mutation populations present in the HRV-39 inoculum in two human subjects and one HeLa-OH supernatant collected 5 days post-infection. The analysis revealed hypervariable mutation locations in VP2, VP3, VP1, 2C and 3C genes and conserved regions in VP4, 2A, 2B, 3A, 3B and 3D genes. These results were confirmed by classical sequencing of additional samples, both from inoculated volunteers and independent cell infections, and suggest that HRV inter-host transmission is not associated with a strong bottleneck effect. A specific analysis of the VP1 capsid gene of 15 human cases confirmed the high mutation incidence in this capsid region, but not in the antiviral drug-binding pocket. We could also estimate a mutation frequency in vivo of 3.4×10⁻⁴ mutations/nucleotides and 3.1×10⁻⁴ over the entire ORF and VP1 gene, respectively. In vivo, HRV generate new variants rapidly during the course of an acute infection due to mutations that accumulate in hot spot regions located at the capsid level, as well as in 2C and 3C genes.

Viral and host proteins involved in picornavirus life cycle

Picornaviruses cause several diseases, not only in humans but also in various animal hosts. For instance, human enteroviruses can cause hand-foot-and-mouth disease, herpangina, myocarditis, acute flaccid paralysis, acute hemorrhagic conjunctivitis, severe neurological complications, including brainstem encephalitis, meningitis and poliomyelitis, and even death. The interaction between the virus and the host is important for viral replication, virulence and pathogenicity. This article reviews studies of the functions of viral and host factors that are involved in the life cycle of picornavirus. The interactions of viral capsid proteins with host cell receptors is discussed first, and the mechanisms by which the viral and host cell factors are involved in viral replication, viral translation and the switch from translation to RNA replication are then addressed. Understanding how cellular proteins interact with viral RNA or viral proteins, as well as the roles of each in viral infection, will provide insights for the design of novel antiviral agents based on these interactions.

Cis-acting RNA elements in human and animal plus-strand RNA viruses

The RNA genomes of plus-strand RNA viruses have the ability to form secondary and higher-order structures that contribute to their stability and to their participation in inter- and intramolecular interactions. Those structures that are functionally important are called cis-acting RNA elements because their functions cannot be complemented in trans. They can be involved not only in RNA/RNA interactions but also in binding of viral and cellular proteins during the complex processes of translation, RNA replication and encapsidation. Most viral cis-acting RNA elements are located in the highly structured 5'- and 3'-nontranslated regions of the genomes but sometimes they also extend into the adjacent coding sequences. In addition, some cis-acting RNA elements are embedded within the coding sequences far away from the genomic ends. Although the functional importance of many of these structures has been confirmed by genetic and biochemical analyses, their precise roles are not yet fully understood. In this review we have summarized what is known about cis-acting RNA elements in nine families of human and animal plus-strand RNA viruses with an emphasis on the most thoroughly characterized virus families, the Picornaviridae and Flaviviridae.

Phylogenetic prediction of cis-acting elements: a cre-like sequence in Norovirus genome?

BACKGROUND

Discrete RNA structures such as cis-acting replication elements (cre) in the coding region of RNA virus genomes create characteristic suppression of synonymous site variability (SSSV). Different phylogenetic methods have been developed to predict secondary structures in RNA viruses, for high-resolution thermodynamic scanning and for detecting SSSV. These approaches have been successfully in predicting cis-acting signals in different members of the family Picornaviridae and Caliciviridae. In order to gain insight into the identification of cis-acting signals in viruses whose mechanisms of replication are currently unknown, we performed a phylogenetic analysis of complete genome sequences from 49 Human Norovirus (NoV) strains.

FINDINGS

The complete coding sequences of NoV ORF1 were obtained from the DDBJ database and aligned. Shannon entropy calculations and RNAalifold consensus RNA structure prediction identified a discrete, conserved, invariant sequence region with a characteristic AAACG cre motif at positions 240 through 291 of the RNA dependant RNA polymerase (RdRp) sequence (relative to strain [EMBL:EU794713]). This sequence region has a high probability to conform a stem-loop.

CONCLUSION

A new predicted stem-loop has been identified near the 5' end of the RdRp of Human NoV genome. This is the same location recently reported for Hepatovirus cre stem-loop.

Novel antiviral agent DTriP-22 targets RNA-dependent RNA polymerase of enterovirus 71

Enterovirus 71 (EV71) has emerged as an important virulent neurotropic enterovirus in young children. DTriP-22 (4{4-[(2-bromo-phenyl)-(3-methyl-thiophen-2-yl)-methyl]-piperazin-1-yl}-1-pheny-1H-pyrazolo[3,4-d]pyrimidine) was found to be a novel and potent inhibitor of EV71. The molecular target of this compound was identified by analyzing DTriP-22-resistant viruses. A substitution of lysine for Arg163 in EV71 3D polymerase rendered the virus drug resistant. DTriP-22 exhibited the ability to inhibit viral replication by reducing viral RNA accumulation. The compound suppressed the accumulated levels of both positive- and negative-stranded viral RNA during virus infection. An in vitro polymerase assay indicated that DTriP-22 inhibited the poly(U) elongation activity, but not the VPg uridylylation activity, of EV71 polymerase. These findings demonstrate that the nonnucleoside analogue DTriP-22 acts as a novel inhibitor of EV71 polymerase. DTriP-22 also exhibited a broad spectrum of antiviral activity against other picornaviruses, which highlights its potential in the development of antiviral agents.

Cis-active RNA elements (CREs) and picornavirus RNA replication

Our understanding of picornavirus RNA replication has improved over the past 10 years, due in large part to the discovery of cis-active RNA elements (CREs) within picornavirus RNA genomes. CREs function as templates for the conversion of VPg, the Viral Protein of the genome, into VPgpUpU(OH). These so called CREs are different from the previously recognized cis-active RNA sequences and structures within the 5' and 3' NTRs of picornavirus genomes. Two adenosine residues in the loop of the CRE RNA structures allow the viral RNA-dependent RNA polymerase 3D(Pol) to add two uridine residues to the tyrosine residue of VPg. Because VPg and/or VPgpUpU(OH) prime the initiation of viral RNA replication, the asymmetric replication of viral RNA could not be explained without an understanding of the viral RNA template involved in the conversion of VPg into VPgpUpU(OH) primers. We review the growing body of knowledge regarding picornavirus CREs and discuss how CRE RNAs work coordinately with viral replication proteins and other cis-active RNAs in the 5' and 3' NTRs during RNA replication.

Picornavirus genome replication: roles of precursor proteins and rate-limiting steps in oriI-dependent VPg uridylylation

The 5' ends of all picornaviral RNAs are linked covalently to the genome-encoded peptide, VPg (or 3B). VPg linkage is thought to occur in two steps. First, VPg serves as a primer for production of diuridylylated VPg (VPg-pUpU) in a reaction catalyzed by the viral polymerase that is templated by an RNA element (oriI). It is currently thought that the viral 3AB protein is the source of VPg in vivo. Second, VPg-pUpU is transferred to the 3' end of plus- and/or minus-strand RNA and serves as primer for production of full-length RNA. Nothing is known about the mechanism of transfer. We present biochemical and biological evidence refuting the use of 3AB as the donor for VPg uridylylation. Our data are consistent with precursors 3BC and/or 3BCD being employed for uridylylation. This conclusion is supported by in vitro uridylylation of these proteins, the ability of a mutant replicon incapable of producing processed VPg to replicate in HeLa cells and cell-free extracts and corresponding precursor processing profiles, and the demonstration of 3BC-linked RNA in mutant replicon-transfected cells. These data permit elaboration of our model for VPg uridylylation to include the use of precursor proteins and invoke a possible mechanism for location of the diuridylylated, VPg-containing precursor at the 3' end of plus- or minus-strand RNA for production of full-length RNA. Finally, determinants of VPg uridylylation efficiency suggest formation and/or collapse or release of the uridylylated product as the rate-limiting step in vitro depending upon the VPg donor employed.

Identification of a conserved RNA replication element (cre) within the 3Dpol-coding sequence of hepatoviruses

Internally located, cis-acting RNA replication elements (cre) have been identified within the genomes of viruses representing each of the major picornavirus genera (Enterovirus, Rhinovirus, Aphthovirus, and Cardiovirus) except Hepatovirus. Previous efforts to identify a stem-loop structure with cre function in hepatitis A virus (HAV), the type species of this genus, by phylogenetic analyses or thermodynamic predictions have not succeeded. However, a region of markedly suppressed synonymous codon variability was identified in alignments of HAV sequences near the 5' end of the 3D(pol)-coding sequence of HAV, consistent with noncoding constraints imposed by an underlying RNA secondary structure. Subsequent MFOLD predictions identified a 110-nucleotide (nt) complex stem-loop in this region with a typical AAACA/G cre motif in its top loop. A potentially homologous RNA structure was identified in this region of the avian encephalitis virus genome, despite little nucleotide sequence relatedness between it and HAV. Mutations that disrupted secondary RNA structure or the AAACA/G motif, without altering the amino acid sequence of 3D(pol), ablated replication of a subgenomic HAV replicon in transfected human hepatoma cells. Replication competence could be rescued by reinsertion of the native 110-nt stem-loop structure (but not an abbreviated 45-nt stem-loop) upstream of the HAV coding sequence in the replicon. These results suggest that this stem-loop is functionally similar to cre elements of other picornaviruses and likely involved in templating VPg uridylylation as in other picornaviruses, despite its significantly larger size and lower free folding energy.

The capsid-coding region hairpin element (cHP) is a critical determinant of dengue virus and West Nile virus RNA synthesis

Dengue virus (DENV) and West Nile virus (WNV) are members of the Flavivirus genus of positive-strand RNA viruses. RNA sequences and structures, primarily in the untranslated regions, have been shown to modulate flaviviral gene expression and genome replication. Previously, we demonstrated that a structure in the DENV coding region (cHP) enhances translation start codon selection and is required for viral replication. Here we further characterize the role of the cHP in the DENV life cycle. We demonstrate that the cHP is required for efficient viral RNA synthesis in a sequence-independent manner. Viruses with a disrupted cHP are rescued by a spontaneous compensatory mutation that restabilizes the structure. Furthermore, the cHP, which is predicted to be conserved among arthropod-borne flaviviruses, is required for WNV replication. We propose that the cHP is a multifunctional determinant of flavivirus replication, functioning in both translation and RNA synthesis.

The cis-acting replication elements define human enterovirus and rhinovirus species

Replication of picornaviruses is dependent on VPg uridylylation, which is linked to the presence of the internal cis-acting replication element (cre). Cre are located within the sequence encoding polyprotein, yet at distinct positions as demonstrated for poliovirus and coxsackievirus-B3, cardiovirus, and human rhinovirus (HRV-A and HRV-B), overlapping proteins 2C, VP2, 2A, and VP1, respectively. Here we report a novel distinct cre element located in the VP2 region of the recently reported HRV-A2 species and provide evolutionary evidence of its functionality. We also experimentally interrogated functionality of recently identified HRV-B cre in the 2C region that is orthologous to the human enterovirus (HEV) cre and show that it is dispensable for replication and appears to be a nonfunctional evolutionary relic. In addition, our mutational analysis highlights two amino acids in the 2C protein that are crucial for replication. Remarkably, we conclude that each genetic clade of HRV and HEV is characterized by a unique functional cre element, where evolutionary success of a new genetic lineage seems to be associated with an invention of a novel cre motif and decay of the ancestral one. Therefore, we propose that cre element could be considered as an additional criterion for human rhinovirus and enterovirus classification.

Poliovirus cis-acting replication element-dependent VPg Uridylylation lowers the Km of the initiating nucleoside triphosphate for viral RNA replication

Initiation of RNA synthesis by RNA-dependent RNA polymerases occurs when a phosphodiester bond is formed between the first two nucleotides in the 5' terminus of product RNA. The concentration of initiating nucleoside triphosphates (NTPi) required for RNA synthesis is typically greater than the concentration of NTPs required for elongation. VPg, a small viral protein, is covalently attached to the 5' end of picornavirus negative- and positive-strand RNAs. A cis-acting replication element (CRE) within picornavirus RNAs serves as a template for the uridylylation of VPg, resulting in the synthesis of VPgpUpU(OH). Mutations within the CRE RNA structure prevent VPg uridylylation. While the tyrosine hydroxyl of VPg can prime negative-strand RNA synthesis in a CRE- and VPgpUpU(OH)-independent manner, CRE-dependent VPgpUpU(OH) synthesis is absolutely required for positive-strand RNA synthesis. As reported herein, low concentrations of UTP did not support negative-strand RNA synthesis when CRE-disrupting mutations prevented VPg uridylylation, whereas correspondingly low concentrations of CTP or GTP had no negative effects on the magnitude of CRE-independent negative-strand RNA synthesis. The experimental data indicate that CRE-dependent VPg uridylylation lowers the K(m) of UTP required for viral RNA replication and that CRE-dependent VPgpUpU(OH) synthesis was required for efficient negative-strand RNA synthesis, especially when UTP concentrations were limiting. By lowering the concentration of UTP needed for the initiation of RNA replication, CRE-dependent VPg uridylylation provides a mechanism for a more robust initiation of RNA replication.

Picornavirus genome replication. Identification of the surface of the poliovirus (PV) 3C dimer that interacts with PV 3Dpol during VPg uridylylation and construction of a structural model for the PV 3C2-3Dpol complex

Picornaviruses have a peptide termed VPg covalently linked to the 5'-end of the genome. Attachment of VPg to the genome occurs in at least two steps. First, Tyr-3 of VPg, or some precursor thereof, is used as a primer by the viral RNA-dependent RNA polymerase, 3Dpol, to produce VPg-pUpU. Second, VPg-pUpU is used as a primer to produce full-length genomic RNA. Production of VPg-pUpU is templated by a single adenylate residue located in the loop of an RNA stem-loop structure termed oriI by using a slide-back mechanism. Recruitment of 3Dpol to and its stability on oriI have been suggested to require an interaction between the back of the thumb subdomain of 3Dpol and an undefined region of the 3C domain of viral protein 3CD. We have performed surface acidic-to-alanine-scanning mutagenesis of 3C to identify the surface of 3C with which 3Dpol interacts. This analysis identified numerous viable poliovirus mutants with reduced growth kinetics that correlated to reduced kinetics of RNA synthesis that was attributable to a change in VPg-pUpU production. Importantly, these 3C derivatives were all capable of binding to oriI as well as wild-type 3C. Synthetic lethality was observed for these mutants when placed in the context of a poliovirus mutant containing 3Dpol-R455A, a residue on the back of the thumb required for VPg uridylylation. These data were used to guide molecular docking of the structures for a poliovirus 3C dimer and 3Dpol, leading to a structural model for the 3C(2)-3Dpol complex that extrapolates well to all picornaviruses.

Human rhinovirus type 14 gain-of-function mutants for oriI utilization define residues of 3C(D) and 3Dpol that contribute to assembly and stability of the picornavirus VPg uridylylation complex

VPg linkage to the 5' ends of picornavirus RNAs requires production of VPg-pUpU. VPg-pUpU is templated by an RNA stem-loop (the cre or oriI) found at different locations in picornavirus genomes. At least one adaptive mutation is required for human rhinovirus type 14 (HRV-14) to use poliovirus type 3 (PV-3) or PV-1 oriI efficiently. One mutation changes Leu-94 of 3C to Pro; the other changes Asp-406 of 3Dpol to Asn. By using an in vitro VPg uridylylation system for HRV-14 that recapitulates biological phenotypes, we show that the 3C adaptive mutation functions at the level of 3C(D) and the 3D adaptive mutation functions at the level of 3Dpol. Pro-94 3C(D) has an expanded specificity and enhanced stability relative to wild-type 3C(D) that leads to production of more processive uridylylation complexes. PV-1/HRV-14 oriI chimeras reveal sequence specificity in 3C(D) recognition of oriI that resides in the upper stem. Asn-406 3Dpol is as active as wild-type 3Dpol in RNA-primed reactions but exhibits greater VPg uridylylation activity due to more efficient recruitment to and retention in the VPg uridylylation complex. Asn-406 3Dpol from PV-1 exhibits identical behavior. These studies suggest a two-step binding mechanism in the assembly of the 3C(D)-oriI complex that leads to unwinding of at least the upper stem of oriI and provide additional support for a direct interaction between the back of the thumb of 3Dpol and 3C that is required for 3Dpol recruitment to and retention in the uridylylation complex.

Universal and mutation-resistant anti-enteroviral activity: potency of small interfering RNA complementary to the conserved cis-acting replication element within the enterovirus coding region

The promising potential of RNA interference-based antiviral therapies has been well established. However, the antiviral efficacy is largely limited by genomic diversity and genetic instability of various viruses, including human enterovirus B (HEB). In this work, the first evidence supporting the anti-HEB activity of the small interfering RNA (siRNA) targeting the highly conserved cis-acting replication element (CRE) within virus coding region 2C is presented. HeLa cells pre-treated with siRNA complementary to the conserved sequence of the loop region of CRE(2C) were effectively rescued from the cytopathic effects of HEBs. Downregulation of virus replication and attenuation of cytotoxicity were consistently observed in various reference strains and clinical isolates. Cells treated with this siRNA were resistant to the emergence of viable escape mutants and showed sustained antiviral ability. Collectively, the data suggest that the siRNA based on the disordered structure within the highly conserved cis-acting coding region has potential as a universal, persistent anti-HEB agent. The same strategy can be successfully applied to the development of siRNA with consistent antiviral effects in other virus groups possessing similar RNA elements.

New complete genome sequences of human rhinoviruses shed light on their phylogeny and genomic features

Background

Human rhinoviruses (HRV), the most frequent cause of respiratory infections, include 99 different serotypes segregating into two species, A and B. Rhinoviruses share extensive genomic sequence similarity with enteroviruses and both are part of the picornavirus family. Nevertheless they differ significantly at the phenotypic level. The lack of HRV full-length genome sequences and the absence of analysis comparing picornaviruses at the whole genome level limit our knowledge of the genomic features supporting these differences.

Results

Here we report complete genome sequences of 12 HRV-A and HRV-B serotypes, more than doubling the current number of available HRV sequences. The whole-genome maximum-likelihood phylogenetic analysis suggests that HRV-B and human enteroviruses (HEV) diverged from the last common ancestor after their separation from HRV-A. On the other hand, compared to HEV, HRV-B are more related to HRV-A in the capsid and 3B-C regions. We also identified the presence of a 2C cis-acting replication element (cre) in HRV-B that is not present in HRV-A, and that had been previously characterized only in HEV. In contrast to HEV viruses, HRV-A and HRV-B share also markedly lower GC content along the whole genome length.

Conclusion

Our findings provide basis to speculate about both the biological similarities and the differences (e.g. tissue tropism, temperature adaptation or acid lability) of these three groups of viruses.

Picornavirus genome replication: assembly and organization of the VPg uridylylation ribonucleoprotein (initiation) complex

All picornaviruses have a protein, VPg, covalently linked to the 5'-ends of their genomes. Uridylylated VPg (VPg-pUpU) is thought to serve as the protein primer for RNA synthesis. VPg-pUpU can be produced in vitro by the viral polymerase, 3Dpol, in a reaction in which a single adenylate residue of a stem-loop structure, termed oriI, templates processive incorporation of UMP into VPg by using a "slide-back" mechanism. This reaction is greatly stimulated by viral precursor protein 3CD or its processed derivative, 3C; both contain RNA-binding and protease activities. We show that the 3C domain encodes specificity for oriI, and the 3D domain enhances the overall affinity for oriI. Thus, 3C(D) stimulation exhibits an RNA length dependence. By using a minimal system to evaluate the mechanism of VPg uridylylation, we show that the active complex contains polymerase, oriI, and 3C(D) at stoichiometry of 1:1:2. Dimerization of 3C(D) is supported by physical and structural data. Polymerase recruitment to and retention in this complex require a protein-protein interaction between the polymerase and 3C(D). Physical and functional data for this interaction are provided for three picornaviruses. VPg association with this complex is weak, suggesting that formation of a complex containing all necessary components of the reaction is rate-limiting for the reaction. We suggest that assembly of this complex in vivo would be facilitated by use of precursor proteins instead of processed proteins. These data provide a glimpse into the organization of the ribonucleoprotein complex that catalyzes this key step in picornavirus genome replication.

Genome-wide diversity and selective pressure in the human rhinovirus

Background

The human rhinoviruses (HRV) are one of the most common and diverse respiratory pathogens of humans. Over 100 distinct HRV serotypes are known, yet only 6 genomes are available. Due to the paucity of HRV genome sequence, little is known about the genetic diversity within HRV or the forces driving this diversity. Previous comparative genome sequence analyses indicate that recombination drives diversification in multiple genera of the picornavirus family, yet it remains unclear if this holds for HRV.

Results

To resolve this and gain insight into the forces driving diversification in HRV, we generated a representative set of 34 fully sequenced HRVs. Analysis of these genomes shows consistent phylogenies across the genome, conserved non-coding elements, and only limited recombination. However, spikes of genetic diversity at both the nucleotide and amino acid level are detectable within every locus of the genome. Despite this, the HRV genome as a whole is under purifying selective pressure, with islands of diversifying pressure in the VP1, VP2, and VP3 structural genes and two non-structural genes, the 3C protease and 3D polymerase. Mapping diversifying residues in these factors onto available 3-dimensional structures revealed the diversifying capsid residues partition to the external surface of the viral particle in statistically significant proximity to antigenic sites. Diversifying pressure in the pleconaril binding site is confined to a single residue known to confer drug resistance (VP1 191). In contrast, diversifying pressure in the non-structural genes is less clear, mapping both nearby and beyond characterized functional domains of these factors.

Conclusion

This work provides a foundation for understanding HRV genetic diversity and insight into the underlying biology driving evolution in HRV. It expands our knowledge of the genome sequence space that HRV reference serotypes occupy and how the pattern of genetic diversity across HRV genomes differs from other picornaviruses. It also reveals evidence of diversifying selective pressure in both structural genes known to interact with the host immune system and in domains of unassigned function in the non-structural 3C and 3D genes, raising the possibility that diversification of undiscovered functions in these essential factors may influence HRV fitness and evolution.

Crystal structure of poliovirus 3CD protein: virally encoded protease and precursor to the RNA-dependent RNA polymerase

Poliovirus 3CD is a multifunctional protein that serves as a precursor to the protease 3C(pro) and the viral polymerase 3D(pol) and also plays a role in the control of viral replication. Although 3CD is a fully functional protease, it lacks polymerase activity. We have solved the crystal structures of 3CD at a 3.4-A resolution and the G64S fidelity mutant of 3D(pol) at a 3.0-A resolution. In the 3CD structure, the 3C and 3D domains are joined by a poorly ordered polypeptide linker, possibly to facilitate its cleavage, in an arrangement that precludes intramolecular proteolysis. The polymerase active site is intact in both the 3CD and the 3D(pol) G64S structures, despite the disruption of a network proposed to position key residues in the active site. Therefore, changes in molecular flexibility may be responsible for the differences in fidelity and polymerase activities. Extensive packing contacts between symmetry-related 3CD molecules and the approach of the 3C domain's N terminus to the VPg binding site suggest how 3D(pol) makes biologically relevant interactions with the 3C, 3CD, and 3BCD proteins that control the uridylylation of VPg during the initiation of viral replication. Indeed, mutations designed to disrupt these interfaces have pronounced effects on the uridylylation reaction in vitro.

Role of RNA structure and RNA binding activity of foot-and-mouth disease virus 3C protein in VPg uridylylation and virus replication

The uridylylation of the VPg peptide primer is the first stage in the replication of picornavirus RNA. This process can be achieved in vitro using purified components, including 3B (VPg) with the RNA dependent RNA polymerase (3Dpol), the precursor 3CD, and an RNA template containing the cre/bus. We show that certain RNA sequences within the foot-and-mouth disease virus (FMDV) 5' untranslated region but outside of the cre/bus can enhance VPg uridylylation activity. Furthermore, we have shown that the FMDV 3C protein alone can substitute for 3CD, albeit less efficiently. In addition, the VPg precursors, 3B(3)3C and 3B(123)3C, can function as substrates for uridylylation in the absence of added 3C or 3CD. Residues within the FMDV 3C protein involved in interaction with the cre/bus RNA have been identified and are located on the face of the protein opposite from the catalytic site. These residues within 3C are also essential for VPg uridylylation activity and efficient virus replication.

Intramolecular and intermolecular uridylylation by poliovirus RNA-dependent RNA polymerase

The 22-amino-acid protein VPg can be uridylylated in solution by purified poliovirus 3D polymerase in a template-dependent reaction thought to mimic primer formation during RNA amplification in infected cells. In the cell, the template used for the reaction is a hairpin RNA termed 2C-cre and, possibly, the poly(A) at the 3' end of the viral genome. Here, we identify several additional substrates for uridylylation by poliovirus 3D polymerase. In the presence of a 15-nucleotide (nt) RNA template, the poliovirus polymerase uridylylates other polymerase molecules in an intermolecular reaction that occurs in a single step, as judged by the chirality of the resulting phosphodiester linkage. Phosphate chirality experiments also showed that VPg uridylylation can occur by a single step; therefore, there is no obligatory uridylylated intermediate in the formation of uridylylated VPg. Other poliovirus proteins that could be uridylylated by 3D polymerase in solution were viral 3CD and 3AB proteins. Strong effects of both RNA and protein ligands on the efficiency and the specificity of the uridylylation reaction were observed: uridylylation of 3D polymerase and 3CD protein was stimulated by the addition of viral protein 3AB, and, when the template was poly(A) instead of the 15-nt RNA, the uridylylation of 3D polymerase itself became intramolecular instead of intermolecular. Finally, an antiuridine antibody identified uridylylated viral 3D polymerase and 3CD protein, as well as a 65- to 70-kDa host protein, in lysates of virus-infected human cells.

Novel, structure-based mechanism for uridylylation of the genome-linked peptide (VPg) of picornaviruses

The VPg peptide, which is found in poliovirus infected cells either covalently bound to the 5'-end of both plus and minus strand viral RNA, or in a uridylylated free form, is essential for picornavirus replication. Combining experimental structure and mutation results with molecular modeling suggests a new mechanism for VPg uridylylation, which assigns an additional function, that of scaffold, to the polymerase. The polarity of the NMR structure of VPg is complementary to the binding site on the surface of poliovirus polymerase determined previously by mutagenesis. Docking VPg at this position places the reactive tyrosinate close to the 5'-end of Poly(A)7 RNA when this is bound with its 3'-end in the active site of the polymerase. The triphosphate tail of a UTP moiety, base paired with the 5'-end of the RNA, projects back over the Tyr3-OH and is held in position by conserved positively charged side-chains of VPg. Other conserved residues mediate binding to the polymerase surface and serve as ligands for metal ion catalyzed transphosphorylation. Additional viral proteins or a second polymerase molecule may aid in stabilizing the components of the reaction. In the model complex, VPg can direct its own uridylylation before entering the polymerase active site.

Polyadenylation of genomic RNA and initiation of antigenomic RNA in a positive-strand RNA virus are controlled by the same cis-element

Genomes and antigenomes of many positive-strand RNA viruses contain 3'-poly(A) and 5'-poly(U) tracts, respectively, serving as mutual templates. Mechanism(s) controlling the length of these homopolymeric stretches are not well understood. Here, we show that in coxsackievirus B3 (CVB3) and three other enteroviruses the poly(A) tract is approximately 80-90 and the poly(U) tract is approximately 20 nt-long. Mutagenesis analysis indicate that the length of the CVB3 3'-poly(A) is determined by the oriR, a cis-element in the 3'-noncoding region of viral RNA. In contrast, while mutations of the oriR inhibit initiation of (-) RNA synthesis, they do not affect the 5'-poly(U) length. Poly(A)-lacking genomes are able to acquire genetically unstable AU-rich poly(A)-terminated 3'-tails, which may be generated by a mechanism distinct from the cognate viral RNA polyadenylation. The aberrant tails ensure only inefficient replication. The possibility of RNA replication independent of oriR and poly(A) demonstrate that highly debilitated viruses are able to survive by utilizing 'emergence', perhaps atavistic, mechanisms.

Structural and functional characterization of the coxsackievirus B3 CRE(2C): role of CRE(2C) in negative- and positive-strand RNA synthesis

A stem-loop element located within the 2C-coding region of the coxsackievirus B3 (CVB3) genome has been proposed to function as a cis-acting replication element (CRE). It is shown here that disruption of this structure indeed interfered with viral RNA replication in vivo and abolished uridylylation of VPg in vitro. Site-directed mutagenesis demonstrated that the previously proposed enteroviral CRE consensus loop sequence, R(1)NNNAAR(2)NNNNNNR(3), is also applicable to CVB3 CRE(2C) and that a positive correlation exists between the ability of CRE(2C) mutants to serve as template in the uridylylation reaction and the capacity of these mutants to support viral RNA replication. To further investigate the effects of the mutations on negative-strand RNA synthesis, an in vitro translation/replication system containing HeLa S10 cell extracts was used. Similar to the results observed for poliovirus and rhinovirus, it was found that a complete disruption of the CRE(2C) structure interfered with positive-strand RNA synthesis, but not with negative-strand synthesis. All CRE(2C) point mutants affecting the enteroviral CRE consensus loop, however, showed a marked decrease in efficiency to induce negative-strand synthesis. Moreover, a transition (A(5)G) regarding the first templating adenosine residue in the loop was even unable to initiate complementary negative-strand synthesis above detectable levels. Taken together, these results indicate that the CVB3 CRE(2C) is not only required for the initiation of positive-strand RNA synthesis, but also plays an essential role in the efficient initiation of negative-strand RNA synthesis, a conclusion that has not been reached previously by using the cell-free system.

Dual modes of modification of hepatitis A virus 3C protease by a serine-derived beta-lactone: selective crystallization and formation of a functional catalytic triad in the active site

Hepatitis A virus (HAV) 3C proteinase is a member of the picornain cysteine proteases responsible for the processing of the viral polyprotein, a function essential for viral maturation and infectivity. This and its structural similarity to other 3C and 3C-like proteases make it an attractive target for the development of antiviral drugs. Previous solution NMR studies have shown that a Cys24Ser (C24S) variant of HAV 3C protein, which displays catalytic properties indistinguishable from the native enzyme, is irreversibly inactivated by N-benzyloxycarbonyl-l-serine-beta-lactone (1a) through alkylation of the sulfur atom at the active site Cys172. However, crystallization of an enzyme-inhibitor adduct from the reaction mixture followed by X-ray structural analysis shows only covalent modification of the epsilon2-nitrogen of the surface His102 by the beta-lactone with no reaction at Cys172. Re-examination of the heteronuclear multiple quantum coherence (HMQC) NMR spectra of the enzyme-inhibitor mixture indicates that dual modes of single covalent modification occur with a >/=3:1 ratio of S-alkylation of Cys172 to N-alkylation of His102. The latter product crystallizes readily, probably due to the interaction between the phenyl ring of the N-benzyloxycarbonyl (N-Cbz) moiety and a hydrophobic pocket of a neighboring protein molecule in the crystal. Furthermore, significant structural changes are observed in the active site of the 3C protease, which lead to the formation of a functional catalytic triad with Asp84 accepting one hydrogen bond from His44. Although the 3C protease modified at Cys172 is catalytically inactive, the singly modified His102 N(epsilon2)-alkylated protein displays a significant level of enzymatic activity, which can be further modified/inhibited by N-iodoacetyl-valine-phenylalanine-amide (IVF) (in solution and in crystal) or excessive amount of the same beta-lactone inhibitor (in solution). The success of soaking IVF into HAV 3C-1a crystals demonstrates the usefulness of this new crystal form in the study of enzyme-inhibitor interactions in the proteolytic active site.

Relationship between poliovirus negative-strand RNA synthesis and the length of the 3' poly(A) tail

The precise relationship between the length of the 3' poly(A) tail and the replication and infectivity of poliovirus RNA was examined in this study. With both poly(A)(11) and poly(A)(12) RNAs, negative-strand synthesis was 1-3% of the level observed with poly(A)(80) RNA. In contrast, increasing the length of the poly(A) tail from (A)(12) to (A)(13) resulted in about a ten-fold increase in negative-strand synthesis. This increase continued with each successive increase in poly(A) tail length. With poly(A)(20) RNA, RNA synthesis approached the level observed with poly(A)(80) RNA. A similar relationship was observed between poly(A) tail length and the infectivity of the viral RNA. A replication model is described which suggests that viral RNA replication is dependent on a poly(A) tail that is long enough to bind poly(A) binding protein and to act as a template for VPg uridylylation and negative-strand initiation.