Next Article in Journal
Inactivation of Dairy Bacteriophages by Thermal and Chemical Treatments
Previous Article in Journal
Animals as Reservoir for Human Norovirus
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Characterization of Mutational Tolerance of a Viral RNA–Protein Interaction

by
Maria A. Prostova
1,2,
Elena Smertina
1,3,
Denis V. Bakhmutov
1,†,
Anna A. Gasparyan
1,4,
Elena V. Khitrina
1,5,
Marina S. Kolesnikova
1,†,
Anna A. Shishova
1,
Anatoly P. Gmyl
1,4,6,† and
Vadim I. Agol
1,5,*
1
Institute of Poliomyelitis, M. P. Chumakov Center for Research and Development of Immunobiological Products, Russian Academy of Sciences, 108819 Moscow, Russia
2
Institute of Molecular Genetics, Russian Academy of Sciences, 123182 Moscow, Russia
3
Faculty of Fundamental Medicine, M. V. Lomonosov Moscow State University, 117192 Moscow, Russia
4
Faculty of Biology, M. V. Lomonosov Moscow State University, 119991 Moscow, Russia
5
A. N. Belozersky Institute of Physical-Chemical Biology, M. V. Lomonosov Moscow State University, 119899 Moscow, Russia
6
Sechenov First Moscow State Medical University, 119991 Moscow, Russia
*
Author to whom correspondence should be addressed.
Deceased.
Viruses 2019, 11(5), 479; https://doi.org/10.3390/v11050479
Submission received: 26 March 2019 / Revised: 20 May 2019 / Accepted: 22 May 2019 / Published: 25 May 2019
(This article belongs to the Section Animal Viruses)

Abstract

:
Replication of RNA viruses is generally markedly error-prone. Nevertheless, these viruses usually retain their identity under more or less constant conditions due to different mechanisms of mutation tolerance. However, there exists only limited information on quantitative aspects of the mutational tolerance of distinct viral functions. To address this problem, we used here as a model the interaction between a replicative cis-acting RNA element (oriL) of poliovirus and its ligand (viral protein 3CD). The mutational tolerance of a conserved tripeptide of 3CD, directly involved in this interaction, was investigated. Randomization of the relevant codons and reverse genetics were used to define the space of viability-compatible sequences. Surprisingly, at least 11 different amino acid substitutions in this tripeptide were not lethal. Several altered viruses exhibited wild-type-like phenotypes, whereas debilitated (but viable) genomes could increase their fitness by the acquisition of reversions or compensatory mutations. Together with our study on the tolerance of oriL (Prostova et al., 2015), the results demonstrate that at least 42 out of 51 possible nucleotide replacements within the two relevant genomic regions are viability-compatible. These results provide new insights into structural aspects of an important viral function as well as into the general problems of viral mutational robustness and evolution.

1. Introduction

Genome replication of RNA-containing viruses is generally characterized by a high level of mutations due to the infidelity of their RNA-dependent RNA polymerases and the overwhelming absence of proofreading. The error rate varies in the range of ~10−3–10−6 nt/site/copying act [1,2,3]. Thus, newly generated molecules of poliovirus RNA (~7.4 × 103 nt) differ from their templates averagely at one nucleotide position [4]. Consequently, even clonal populations of such viruses exist as highly heterogeneous swarms of variants referred to as quasispecies [5,6]. Nevertheless, under constant conditions, RNA viruses demonstrate a marked conservation of their genetic structure owing mainly to the selective disadvantage of debilitated mutants. However, viral intra- and inter-host transmissions often involve bottlenecks, upon which haphazardly picked-up genomes establish new lineages. Although such situations may contribute to the viral evolvability, they are also associated, especially if occurring repeatedly, with high risks of severe fitness cost or even the dead-end of lineages, a phenomenon known as the Muller ratchet [7,8]. To counteract or minimize effects of such disadvantageous developments, viruses possess mechanisms ensuring their relative genetic robustness, which include such well-known phenomena as degeneracy of the genetic code or of spatial structures of functional elements of RNA or protein molecules [9]. However, there exists only limited information concerning detailed and quantitative aspects of the mutational tolerance of distinct viral functions.
Interactions between elements of viral RNA and proteins required for efficient viral reproduction are attractive objects with which to study this problem. Quantitative evaluation of functional effects of mutational alterations of these partners can provide deeper insights into not only the nature of robustness of these interactions but also of their mechanistic details. As a model, we used here the poliovirus, a member of the Enterovirus genus of the Picornaviridae family of animal viruses. Its positive-strand RNA genome harbors a single open reading frame (ORF) encoding the polyprotein, which is eventually processed into a dozen mature proteins [10,11]. Among others, this set includes the RNA-dependent RNA polymerase 3Dpol, the viral primer protein VPg (also known as 3B), and the protease 3Cpro responsible for the polyprotein processing. The ORF is flanked by the 5′- and 3′-untranslated regions, containing, in particular, two cis-elements required for efficient RNA replication: oriL, also known as the cloverleaf because of its secondary structure [12,13], and oriR, respectively (Figure 1).
In particular, oriL is involved in initiation of the synthesis of both the viral (positive) and complementary (negative) RNA strands [13,14,15,16]. An essential step in accomplishing these functions is by recruiting oriL of several viral and host proteins to generate a ribonucleoprotein (RNP) complex [13,17,18]. A key reaction in the formation of the relevant RNP is the interaction of the hairpin-like domain d of oriL with the viral protein 3CD, the uncleaved precursor of 3Cpro and 3Dpol [17,18,19,20] (Figure 1). Specifically, the apical tetraloop of domain d is believed to interact with the TGK motif of the 3C moiety of 3CD, since alterations in the structure of the former may be functionally compensated by the mutational changes in the latter [20,21,22].
In a previous study [20] we assayed the level of mutational robustness of the domain d of oriL by the randomization (in the context of the full genome) of the sequence of its tetraloop and two adjoining base pairs, followed by selection and investigation of the viable progeny (the SELEX procedure). In addition to such in vitro selection experiments, certain altered genomes were engineered and investigated. Unexpectedly, it was found that the overwhelming majority, if not all, of the possible point mutations of this highly conserved octanucleotide were compatible with the viral viability. Some of these mutations did not change the viral phenotype, whereas others inflicted functional defects of variable strength, as evidenced by the impairment of viral reproduction, RNA synthesis, and oriL/3CD binding, these three alterations well correlating with each other. The debilitated mutants demonstrated a good potential for phenotypic recovery by either acquisition of advantageous mutations in the octanucleotide or by changes in the TGK motif of the 3CD protein, and in particular its conversion into IGK or TGR.
To further characterize the mutational tolerance of the poliovirus oriL/3CD interaction, the TGK-encoding nonanucleotide of 3C was now subjected to randomization. After deriving some preliminary conclusions from the selection and sequencing of viable viruses, an additional set of mutant genomes was engineered and investigated. At least 11 amino acid replacements in TGK were found to be compatible with the virus viability. The results markedly enlarged the space of 3CD sequences compatible with the virus viability and in combination with our previous study on mutational tolerance of oriL [20], they demonstrated that in the two genome segments controlling the tetraloop/TGK interaction (encompassing the apical octanucleotide of domain d of oriL and the TGK-encoding nonanucleotide, together 17 nt-long), at least 42 out of 51 theoretically possible nucleotide replacements are not lethal. These results not only provide further insights into the mechanistic aspects of an important step of enterovirus reproduction but also give a deeper understanding of the mutational tolerance and evolvability of RNA viruses.

2. Materials and Methods

2.1. Generation of Plasmids Encoding Poliovirus Genomes with Randomized Segments

Plasmids with the randomized nucleotides encoding amino acids 154–156 of the viral 3C protein were prepared as follows (Figure 2). A 68 nt-long synthetic DNA 3CTGK (Table 1) corresponding to the region with coordinates 5875–5942 of the poliovirus RNA but with the randomized relevant nonanucleotide (positions 5897–5905) and containing a synonymous marker mutation G5914C, was purchased from Syntol (Moscow Russia).
For its PCR-amplification, the reaction mixture containing 0.6 pmole 3CTGK and 20 pmoles each of 3CMluI and 3CSacII primers (Table 1) was subjected to 10 cycles of heating/cooling (20 s at 55 °C, 20 s at 72 °C, 15 s at 95 °C) and the product was treated with SacII and MluI endonucleases (Fermentas, Vilnius, Lithuania) and purified by electrophoresis in 2.5% agarose. The DNA fragment thus generated was inserted into poliovirus genome-harboring plasmids in a two-step procedure aimed to minimize possible contamination with non-mutated genomes. First, it was inserted at the proper position into the plasmid pT7PV1RibMS, which contained the full-genome copy of the type 1 poliovirus RNA with the artificially created unique cleavage sites for the restriction endonucleases MluI (at position 5861/5862) and SacII (5952/5953). To this end, the large fragment obtained after treatment of pT7PV1RibMS with MluI and SacII was purified by electrophoresis in 1% agarose and ligated with the above-described fragment containing the randomized nonanucleotide. The created plasmid was used for the transformation of E. coli TOP10 (Invitrogen, Carlsbad, CA, USA) and the relevant region of a clone (pT7PV1RibMSa, clone #18 in Table 3) was sequenced and found to lack G5906, resulting in the frameshifting of the open reading frame downstream of the randomized nonanucleotide. This intermediate plasmid was used as a donor of the large SacII/MluI fragment in the ligation reaction with the small SacII/MluI fragment harboring the randomized codons, as described above. Pools of the randomized plasmids were obtained upon transformation of E. coli TOP10 cells.
To generate plasmids with randomized octanucleotide of oriL and nonanucleotide of the 3C gene, clone #18 (pT7PV1RibMSa) was treated with ApaI and SplI and the larger ApaI/SplI fragment was ligated with the small ApaI/SplI fragment obtained by PCR with the previously described [20] synthetic DNA possessing randomized relevant octanucleotide in oriL. The plasmid thus generated was propagated in E. coli and pooled ~2250 plasmids (pT7PV1RibMSb) were further treated as pT7PV1RibMSa in the protocol for the generation of plasmids with the randomized 9 nt in 3C.

2.2. Construction of Genomes with Desired Mutations

Since the viruses generated by the SELEX procedure not infrequently contained mutations outside the randomized nonanucleotide, certain analogous genomes lacking these undesired mutations were constructed. To this end, RNA from the preparations of the relevant viruses generated from the randomized plasmids was reverse transcribed and PCR-amplified using 3CMluI and 3CSacII primers. The product was digested with MluI and SacII and used to generate constructs encoding the full-length viral genome, as described above.
To generate plasmids encoding viral genomes with novel mutations, two separate PCR amplifications were performed using pT7PV1RibMS as the template and pairs of primers, mut-s and 3EP4, in one reaction, and mut-a and B5594 in the other, where mut-s and mut-a harbored the desired mutation (Table 1). The products obtained in these reactions were fused together in an additional PCR with 3EP17 and B5594 primers, digested with SacII and MluI, and used to generate plasmids encoding the full-length viral genome as above.
For the experiments requiring purification of mutated 3CD proteins, the small MluI-SacII fragment with the desired mutation was generated by treatment with MluI and SacII of the PCR product of the reaction containing pT7PV1RibMS/mut as the template and 3CMluI and 3CSacII as primers. The gel-purified fragments were ligated to the large fragment obtained upon digestion with the same endonucleases of the expression plasmid pQE60-3CD [20], which harbored the 3CD sequence of wild-type poliovirus with the H40A mutation (inactivating the proteolytic activity of the protein) and the N-terminal His-tag.
The presence of introduced mutations in all the engineered plasmid clones was verified by sequencing.

2.3. E. coli Transformation and Preparation of DNA

E. coli TOP10 cells were transformed by the plasmids and the plasmid DNA was isolated from individual or pooled clones as described previously [20].

2.4. Transcription and Transfection

The procedures were performed as described in [20]. Briefly, DNA samples from 19 plasmid pools containing from 31 to ~2250 variants were linearized by digestion with EcoRI and transcribed by T7pol (Thermo Scientific, Carlsbad, CA, USA), and the transcription mixtures were used for the transfection of Vero cells (ATCC® CCL-81™) as described [23]. The concentration of RNA in the transcripts was evaluated by EtBr staining in 1% agarose gels. For the experiments with engineered genomes, the viral RNA was purified by centrifugation in 5–20% sucrose gradients and the concentration of RNA was determined spectrophotometrically.

2.5. Sequencing of the Viral Genomes

The material from a plaque or virus suspension was suspended in 0.3 mL of the nutrient medium. RNA was isolated using either a phenol-chloroform extraction, Trizol reagent (Invitrogen), or Qiagen RNAeasy kit (Qiagen, Hilden, Germany), and was reverse transcribed using random or 3EP4 primers (Table 1). The oriL-containing PCR products were prepared with primers Rib2 and DEN3, whereas primers B5594 and 3EP17 were used for the preparation of DNA fragments encoding the portions of 3C that included the relevant tripeptide. The PCR products were gel-purified and sequenced either manually using afmol® DNA Cycle Sequencing System (Promega, Fitchburg, WI, USA) or by automatic sequencers Beckman Coulter Seq 8000 or ABI 3130 Genetic Analyzer.

2.6. Time-Course of Viral RNA Replication

The experiments were performed as described in [20] with minor modifications. Vero cells monolayers grown in 12-well panels (Corning Incorporated, Corning, NY, USA), ~2.4 × 105 cells/well, were transfected with 50 ng RNA transcripts per well, and the total RNA was extracted with Trizol reagent at 0, 12, 16, 20, and 24 h post-transfection (p.t.). Three wells were used as parallels for each time point. One μg of the purified RNA was used for reverse transcription with a random hexamer primer and the Maxima Reverse Transcriptase (Thermo Fisher Scientific, Carlsbad, CA USA). The standard curve was generated by serial dilutions of the wild-type transcript supplemented with 1 μg RNA from mock-transfected Vero cells. Real-time PCR was carried out using a ABI 7500 Real Time PCR System analyzer with primers PVL1 and PVR1 and FAM-tagged oligonucleotide PVP1 as the probe (Table 1).

2.7. Expression and Purification of Modified 3СD Proteins

The procedures used have been described previously [20]. Briefly, pQE60-3CD with the desired substitutions was used for the transformation of E. coli JM109 cells. The transformed cells in 200 mL of SOB medium (2% tryptone, 0.5% yeast extract, 10 mM NaCl, 2.5 mM KCl) containing 20 mM glucose and 100 μg/mL ampicillin were grown in a rotary shaker for 4–5 h at 37 °C. The medium was changed to the one lacking glucose and containing 2 mM IPTG and incubated in a rotary shaker overnight at room temperature. The cells were subjected to centrifugation and lysed by sonication, and recombinant 3CD was purified by the Ni-chelating chromatography. The electrophoretic pattern of the purified proteins is shown in Figure 3.

2.8. Electrophoretic Mobility Shift Assay (EMSA)

Direct interaction between the oriL-containing fragments of viral RNA and variants of the 3CD proteins was investigated as described in [20]. Briefly, motilities of RNA fragments corresponding to purified preparations of the 5′-terminal 115 nt of viral RNAs and 3CD proteins were determined by native 6% PAGE separately and in mixtures.

2.9. 3C Structure Analysis

3C structure (PDB ID 1L1N, chain A) visualization and distance labeling was performed with VMD software [24]. Electrostatic potentials at the solvent accessible surfaces of 3C were calculated using a PyMOL APBS Electrostatics Plugin with default settings [25].

3. Results

3.1. Generation of Plasmids Encoding Full-Length Poliovirus Genomes with a Randomized Tripeptide Corresponding to the TGK-Motif of the 3C Protein

As described in Materials and Methods, a set of plasmids was engineered, in which nine nucleotides (positions 5897–5905 of the poliovirus genome) encoding tripeptide TGK (positions 154–156 of the viral 3C protein) were randomized. The randomization was intentionally incomplete: the triplets in the inserted nonanucleotide contained neither A nor T at the thid positions. Although the total number of potentially encoded codons decreased from 64 to 32, their capacity to translate into all 20 amino acids was retained (the UAA and UGA stop-codons were not encoded). The decrease in the number of codons made the possible nucleotide sequence space more compact, enhancing the probability of representation of codons for all amino acids in relatively small samples of the genomes investigated.
To assess the extent of randomization, the segment corresponding to positions 5800–5960 of the poliovirus genome was sequenced in 24 randomly selected plasmids. Table 2 shows that the distribution of nucleotides across the randomized region was rather uneven, but taking into account the relatively small size of the investigated set, the level of randomization could be considered satisfactory.
Sequencing demonstrated that six plasmids did not harbor non-intended alterations in the codons of interest and adjoining regions, two plasmids exhibited some heterogeneity in the randomized sequence, and one plasmid possessed a synonymous substitution outside the randomized region (Table 3). These nine plasmids could be regarded as consistent with the goal of the experiment. One plasmid has a non-synonymous mutation outside the randomized octanucleotide, whose phenotypic effect, if any, is unknown. However, three plasmids possessed termination codons within the randomized region, and 11 plasmids had deletions, or deletions coupled with insertions, resulting in full or partial frameshifting. Thus, less than a half of the randomized plasmids could be considered as fulfilling the requirements of the experiment.

3.2. Selection of Viable Mutated Viruses and Primary Structures of Their Presumptive oriL-Interacting Region of 3C Proteins

Nineteen pools were assembled from different numbers (31 to ~2250) of randomly selected plasmid clones. DNA preparations isolated from these pools were transcribed in vitro and serial dilutions of these transcripts were used to transfect Vero cells. The cell cultures were observed for at least six days. All these pooled transcripts generated plaque-forming viruses, however, with variable specific infectivity, i.e., the number of plaque-forming units (pfu) per µg RNA, which was several orders of magnitude lower than that of the wild type virus (Table 4), suggesting a very low abundance of the infectious genomes in the samples. The plaques generated by different samples appeared at different times and had different sizes, suggesting the possibility that the plaque-forming capacity of some pools could be due to quasi-infectious genomes, i.e., marginally replicating RNAs, which could obtain plaque-forming capacity after the post-transfectional acquisition of additional mutations.
To characterize the variety of sequences compatible with the viral viability, a portion of the 3C-encoding region (positions 5800–5960) of the RNA from 37 randomly selected mutant clones was sequenced. In 17 such clones, the primary plaques generated upon transfection contained sufficient material for sequencing, but in other cases the material from primary plaques was subjected to one or more bulk passages. The investigated set of genomes contained 22 unique nucleotide sequences.
The prevalence of different codons in the randomized RNA region of these viruses differed markedly from that in the plasmids used for their generation (compare Table 2 and Table 5). For example, the nucleotides A and C that were underrepresented at certain positions of codons 154 and 156 of the plasmids were abundant enough in the sequenced viral RNAs. Of note is that the presence of U and A at the 3rd positions of codons 154 and 156, respectively, may indicate that the genomes of some selected viruses had indeed acquired mutations after the transfection, since the paternal plasmids should not possess the relevant nucleotides at these positions. In other words, the relevant plasmids were likely quasi-infectious. Alternatively, the mutations could be due to errors arising upon in vitro transcription.
In the selected viruses, the tripeptide at positions 154–156 of 3C was represented by 13 unique variants of amino acid sequences, including the parental TGK (Table 6). Although this set did not necessarily denote the whole variety of acceptable sequences, it revealed several important features. The tripeptide of all isolates invariably had the central Gly residue. In the genetically stable variants, position 156 of 3C could be occupied not only by Lys but also by another positively charged residue, Arg, accompanied with either the original Thr or with Val, Ile, or Cys at position 154. The latter position could also be solitarily changed into Val or Cys. Viable but genetically unstable isolates could also contain Met or Leu at position 154 (in combination with Arg156), which were converted upon passages into more comfortable Thr or Val. Similarly, the heterogeneity at this position in one isolate (#17) was likely due to the conversion of poor but acceptable Met in the virus that initiated infection into Val soon after the transfection. Three unstable primary isolates (with VGM, CGC, and IGW) harbored changed amino acids at both positions 154 and 156; they increased their fitness by the acquisition of a positively charged residue, Lys or Arg, at the latter position. Somewhat distinctly, a primary isolate (#21) with Ser156 (and Thr154) retained this residue upon passages but acquired Arg at the position just preceding the randomized triplet. Another virus (#22), having uncomfortable Met156, already possessed Arg153 upon isolation from the primary passage and appeared to convert, during further passages, its acceptable RTGM tetrapeptide into an even better (wild-type) CTGK.
The above results suggest some preliminary regularity. Various amino acids at position 154 are compatible with viral viability but certain ones (Thr, Val, Ile, Cis) are obviously preferable over some others (Leu and Met), the latter tending to be replaced by the former upon passaging. Position 155 appeared to be invariably occupied by Gly. Position 156 requires either Lys or Arg but tolerated to some extent non-positively charged amino acids (e.g., Met, Ser, or Trp), especially if Cys153 was replaced by Arg.

3.3. Engineering and Properties of Mutants with Definite Alterations in the 3C Protein

To define more exactly the validity of the above-formulated preliminary regularity and to obtain some information on the phenotypic properties of viruses with different 3C sequences, a number of mutant genomes were engineered and their plaque-forming capacities were compared with those obtained through the in vivo selection of the randomized RNA. This was important because the latter procedure could be accompanied by the acquisition of undetected mutations outside the randomized sequence. The results obtained with the engineered genomes are presented in Figure 4 and Table 7.
Firstly, mutants with substitutions at position 154 of 3C were engineered. In line with the results obtained with the randomized plasmids, the variants having at this position Val, Ile, or Cys exhibited phenotypes (specific infectivity as well as the time of appearance and size of plaques) comparable to those of the wild-type genome and were genetically stable. The wild-type-like phenotype was exhibited also by the virus with the TGR-encoding RNA. The replacement of Thr154 by Met generated a less-fit, unstable genome, which increased its fitness upon passages by either true reversion, or changing this residue to Val or Leu, or by acquiring a mutation (Pro88Ser) in an upstream region of the 3C sequence. However, the genome with Leu154 also exhibited decreased specific infectivity and plaque sizes as well as genetic instability (manifested by either Leu154Val replacement or the above-mentioned mutation in the nonrandomized area). The SGK-encoding genome appeared to be quasi-infectious, generating plaques at day 8 p.t. in one out of three attempts. This genome should experience two transversions after transfection to acquire a more favorable C154 residue, explaining the poor reproducibility of the experiments.
The replacement of Lys156 by Met or Ser resulted in a marked decrease in both specific infectivity and fitness, as evidenced by the very late (days 7–8 p.t.) appearance of plaque-producing viruses with either the true reversion to TGK (in the case of TGM) or a variety of compensatory mutations in other parts of 3C (in both viruses) (Figure 4B, Table 7). The set of such mutations included the appearance of positively charged Arg at position 153 just preceding the relevant tripeptide as well as various mutations in the 3C sequence upstream (Pro88 to Leu or Thr; Val101 to lle) or downstream of it (Asp282 to Asn). The appearance of Arg153 in similar circumstances was observed in the above experiment with randomization as well. This prompted us to construct genomes with RTGM and RTGS sequences. Such viruses appeared to be markedly fitter compared to their Cys153-containing counterparts, as evidenced by earlier appearances of plaques, but still were genetically unstable, accumulating various mutations in other parts of 3C (Figure 4C, Table 7). These compensatory changes included mutations at already mentioned positions 88 and 101 as well as various other positions of 3С.
Notably, introduction of mutations P88S or I151L into the RTGS-containing RNA resulted in a significant gain of fitness (Figure 4C). The replacement of Lys156 by a negatively-charged Glu in 3C with Thr154 or Val154 generated non-infectious transcripts (Table 7).
Several genomes were reconstructed, in which alterations were introduced at positions 154 and 156 simultaneously (Figure 4A, Table 7). Some of the low-fit, genetically unstable genomes of the above set of in vivo selected viruses were reconstructed and additional variations were introduced at these positions. The MGR-containing virus exhibited a wild-type-like specific infectivity, but generated plaques later and was eventually transformed into VGR- or TGR-containing variants (similarly to the one selected from the randomized genomes). The genomes with IGW, CGC or VGM were unstable, producing late small plaques with a newly-acquired “good” Lys156 or Arg156.
It should be noted that the modifications of the TGK in the investigated cases were not accompanied with the acquisition of any alterations in its ligand, the domain d of oriL. In particular, this was demonstrated for the viruses generated by transcripts with TGM, RTGS, TGR, VGK, CGK, SGK, and IGK after a round of reproduction and viruses with MGR, MGK, IGW, LGK, and RTGM after one additional passage (not shown).

3.4. Checking Mutual Compatibility of Certain Altered 3СD Proteins with Altered oriL

As indicated in the introduction, previous studies demonstrated that adverse effects of certain alterations in the tetraloop of domain d could be, at least partially, compensated by the acquisition of some changes in the TGK motif of 3C, such as T154I and K156R [20,21,22]. Here, we investigated whether alterations of 3C could also increase fitness of the debilitated genomes with certain non-YNMG-like structures (Y = pyrimidine; N = any nucleotide; M = A or C) of the loop of domain d, for which such a compensatory effect was not observed earlier. The genomes were constructed with IGK- or TGR-containing 3C in combination with either agCUUGcu- or auGAGAgu-containing tetraloops of oriL. The former and the latter oriLs appeared to confer, to the TGK-harboring genomes, low-infectious and quasi-infectious phenotypes, respectively, and both could increase their fitness by some alterations of the mutated tetraloop [20]. As shown in Figure 5, both IGK and TGR variants of the genomes with the non-YNMG oriLs demonstrated a markedly higher fitness compared to their TGK-containing counterparts, with a somewhat lesser effect in the case of the auGAGAgu/TGR combination.
To obtain additional information about the mutual compatibility of structurally variable oriL/3CD ligands, plasmids were constructed that encoded full-length poliovirus genomes, in which interacting motifs of both partners (the tetraloop and two adjoining base pairs of domain d of oriL and 9 nt corresponding to codons 154–156 of the 3C gene) were randomized. Seven transcripts of pooled plasmids (~103–1.5 × 104 clones) were used for the transfection of Vero cells. RNA from one of the three pools with ~103 plasmid clones as well as all transcripts of larger pools proved to be infectious, although with a low specific infectivity. Sequencing 113 5′-terminal nucleotides (including oriL) and positions 5800–5960 (a region of 3C gene) of the genome of 12 viral isolates gave the results shown in Table 8. The relevant tripeptide in 3C was represented only by the above described TGR, VGR, and LGR, the latter being replaced by the stronger VGR between the third and fifth passages. The tetraloops of domain d in all but one of the primary isolates possessed sequences not belonging to any consensuses with known stable structures, but they clearly tended to acquire YNMG or GNUA sequences, which have been demonstrated to be optimal for the oriL/3CD interaction [20]. The results indicated that VGR (similarly to the previously investigated TGR and IGK) was able to ensure certain levels of functional interaction with “bad” tetraloops, but improvements of the latter were required to ensure a more efficient oriL/3CD interaction. It may be noted that the unintended mutations in the oriL of the transcripts in pools 1 and 4 (Table 8) could hardly markedly influence the efficiency of this interaction, though the former (insertion of G at position 14) should increase the size of the loop of the hairpin b.
Obviously, the set of available viable viruses was too small to demonstrate the whole space of acceptable combinations.

3.5. Efficiency of RNA Replication of Engineered Mutants

To ascertain whether the observed alterations in the phenotypic properties caused by the mutations in the conserved tripeptide TGK were due to the changed efficiency of the genome replication, the time-course of accumulation of viral RNA in Vero cells transfected with some engineered RNAs was assayed by the quantitative PCR. Generally, the efficiency of replication correlated reasonably well with the mutant’s plaque phenotypes. The genomes with Val, Ile, and Cys at position 154 displayed almost the same replication dynamics as their wild-type Thr-containing counterpart (Figure 6A). The occurrence of Leu (Figure 6A), Met, or Ser (Figure 6B) at this position resulted in a marked decrease in the replicative capacity. Substitution of Lys156 by Arg was without an appreciable effect on the RNA replication (Figure 6B). Non-positively charged amino acids (Met, Ser, or Glu) at this position allowed only a marginal level of accumulation of RNA (Figure 6C). Binary mutations at positions 154 and 156 in MGR, SGR (Figure 6B) and VGE (Figure 6C) endowed the transcripts with somewhat variable but decreased efficiency. However, the replacement of wild-type Cys by Arg at position 153 of the TGM- and TGS-encoding genomes, resulting in some fitness gain, was accompanied only with a very slight, statistically insignificant enhancement of RNA accumulation (Figure 6C). Interestingly, this group with a severely damaged replicative capacity contained both low-fit but viable (possibly quasi-infectious) (TGM and TGS) and apparently dead (TGE and VGE) genomes.
Importantly, the IGK-containing mutant was able to markedly, though not fully, rescue the profoundly inhibited replication of the RNA harboring a GAGA (i.e., GNRA-type) tetraloop in oriL (Figure 6D), congruently with the T154I mutation previously observed to take place upon transfection of the GAGA-containing quasi-infectious genome [20].

3.6. Interaction of Mutant 3CDs with oriL

To ascertain whether phenotypic changes caused by the TGK mutations were linked to the impaired oriL/3CD interaction, electrophoretic mobility shift assays (EMSA) were performed. As interacting partners, we used the recombinant 3CD proteins and 5′-terminal 113 nt-long fragment of the poliovirus RNA. The results were visualized by staining with ethidium bromide (EtBr). The VGK-, IGK-, and TGR-containing proteins demonstrated a level of oriL-binding comparable (in the case of IGK somewhat lower than in other cases) to that of the TGK-containing control (Figure 7) in sufficiently good agreement with their phenotypes. SGK-, MGR-, and TGM-containing 3CD, conferring markedly low viral fitness, exhibited undetectable or marginal affinity to oriL in EMSA. The presence of positively charged Arg153 enhanced, though slightly, binding of the “poor” ligands. Only the CGK–containing 3CD did not show correlation between the EMSA and other phenotypes: it failed to detectably bind oriL but endowed the virus with acceptable fitness. The reason for this discrepancy is unknown, but it might be caused by the presence of two neighboring Cys residues (at positions 153 and 154), which could change the protein conformation under the conditions of the EMSA assay (see also the Discussion section).

4. Discussion

4.1. The Problem

The replication of RNA-containing viruses is generally error-prone, and in the case of picornaviruses nearly each act of template copying may be associated with the acquisition of a mutation [1,2,3]. Such negligence is mainly due to a low fidelity of the viral RNA-dependent RNA polymerases (RdRP) [26,27,28] and the lack (with a very few exceptions) of proofreading mechanisms. This infidelity is not an inherently incorrigible property of RdRP, since their faithfulness could be markedly enhanced by various point mutations [29,30,31,32,33]. Somewhat counterintuitively, an increase in fidelity may result in a decreased viral fitness [29,34,35,36] and hence the frequency of RdRP-made errors appears to be evolutionally tuned.
To prevent or diminish the potential harm of replicative infidelity, viruses should possess a significant degree of mutational tolerance. This tolerance is largely due to the degeneracy of codons, phenotype-neutral character of many amino acid substitutions, the ability of diverse sequences in RNA regulatory elements to maintain analogous mutual orientations, and functional equality of certain nucleotides in these elements [9]. Though the general importance of these factors for the counteracting replicative infidelity is well appreciated, only rather limited information is available on the quantitative aspects of the mutational tolerance of distinct viral functions.

4.2. Mutational Tolerance of the TGK Motif and the oriL/3CD Interaction

The interaction between oriL and 3CD is an essential step of the poliovirus genome replication [13,17,18,20,37]. This interaction involves the tetraloop of domain d of oriL and TGK tripeptide of the 3C moiety of 3CD [16,18,19,20,21,22,38]. Notably, the TGK motif is highly conserved in 3C proteins of members of the Enterovirus C species, though Thr, being most abundant, could also be occupied by Val, Ile, and Met [38]. This study provides insights into quantitative and mechanistic aspects of the mutational tolerance of the genome regions controlling the structure of these ligands.
As summarized in Figure 8, at least 11 nucleotide mutations out of 27 possible in the TGK-encoding nonanuclotide are compatible with the viral viability. Taken together with our previous observation of the mutational robustness of domain d of oriL [20], this means that at least 34 point mutations out of 51 possible in the two-segmented 17 nt-long stretch of RNA (octanucleotide of domain d and nonanucleotide of the 3C gene) are not lethal. If a second point mutation in the TGK-encoding motif is allowed (such mutations could well be already present in the quasispecies populations), then the number of the viability-compatible substitutions in it would reach at least 19 (Figure 8) and the whole space of the permitted nucleotide replacements in the 17 nt-long stretch of RNA would rise to at least 42. Additionally, the tripeptide can sustain not less than 11 amino acid replacements, these being six and five at positions 154 and 156, respectively.
Being not lethal, the amino acid replacements in 3C exerted different fitness effects. A significant proportion of them did not demonstrate, in our in vitro experiments, any marked adverse effects. Other mutations negatively affected the oriL/3CD interaction to different degrees, with some of them bringing the virus on the verge of a catastrophe. However, even in the most debilitating cases, the surviving viruses have a resilience tool: the infidelity of RNA replication resulting in the acquisition of reversions or compensatory mutation.
If there exist such a variety of structures of the relevant tripeptide in protein 3C with apparently more or less equal phenotypic impacts, why is TGK so strictly conserved in wild-type polioviruses? It may be speculated that the laboratory assays do not completely reflect fitness of circulating viruses. It should also be taken into account that even in tissue culture experiments competitive capacity of the relevant mutant viruses has not been assayed.
Though we are focusing here on the direct interaction between the tetraloop of domain d of oriL with the TGK motif of 3CD, it should be kept in mind that the both RNA and protein partners of this interaction have several separate functions and that the formation of the oriL/3CD complex involves several other viral and cellular participants and is significant not only for the RNA replication but for its translation as well [39,40,41]. The viability-compatible mutations identified in this study may affect some of these activities but obviously such effects, if any, are not virus-killing.

4.3. Possible Mechanistic Features of the oriL/3CD Interaction

Although the significance of the oriL/3CD interaction for viral RNA replication is well established, detailed information about the mechanistic aspects underlying their mutual affinity is lacking. The results reported here and in our previous paper [20], though insufficient to suggest a specific molecular model of this interaction, may nevertheless contribute to the development of such models in future.
In particular, the requirements for distinct amino acids at positions 154–156 of poliovirus 3C became partially defined. Although all full wild-type poliovirus genomes in the NCBI database have TGK in the corresponding region, the tripeptide could endure numerous modifications (Figure 8) either without any appreciable loss of fitness or with some debilitating but still viability-compatible effects. Only the central Gly155 appeared to be indispensable, although, admittedly, no exhaustive attempts to prove this were undertaken. The strong requirement for this residue may be related to its position at the loop between the two β-strands [42]. Gly is frequently found in loops because it provides a high flexibility to peptide chains and is often conserved as a structure determinant [43]. “Good” residues at position 154, Val and Ile, share with the wild-type Thr a methyl group at the β carbon atom, hinting that this group may be involved in a hydrophobic interaction. The sebilitating effect of Ser154 is in line with this assumption. On the other hand, Cys154, which was also able to confer a stable wild-type phenotype, has an SH group at the β carbon. It is tempting to assume that this distinction is responsible for a weak interaction of the CGK-containing 3CD with domain d in the EMSA assay (Figure 7). The discrepancy between this inefficiency and functional competence in the RNA replication (Figure 6) may be due to the presence of two neighboring Cys residues (see above).
Discussing the phenotypic effects of 3C mutations, additional possibilities to accomplish the oriL/3C interaction, e.g., via another RNA-binding motif of 3C, KFRDI at positions 82–86 [22,44], should be taken into account. Adaptive changes of Pro88 into Ser, Thr, or Leu observed in several viruses with unfavorable tripeptides at positions 154–156 (Table 7) may presumably be linked to the proximity of position 88 to Tyr6 and His89, involved in oriL recognition [45]. Pro88, being located in a small helix, can affect the orientation of the neighboring His89, which is known to interact with Tyr6 of 3C [42], the distance between their aromatic rings being 3.31Å (Figure 9A), which is common for stacking. These two residues have been reported also to be involved in the oriL/3C interaction [45] and are highly conserved in polioviruses [38]. The close proximity of His89 to TGK (6.43 Å and 10.66 Å to Lys156 and Cys153, respectively) and Tyr6 to Gly155 (5.28 Å) points to possible effects of substitutions in CTGK to the mutual orientation of His89 and Tyr6, which could be compensated by substitution of Pro88 by more conformationally flexible Ser, retaining a His89/Tyr6 interaction in the RNA-recognition.
For full functionality, position 156 could be occupied not only by the wild-type Lys but also by positively charged Arg, whereas the negatively charged Glu at this position was lethal, suggesting an electrostatic interaction in the tetraloop/3CD affinity. It may also be noted that the Lys156Ala replacement was reported to inhibit the capacity of 3CD to stimulate uridylylation of VPg [42], which is known to depend on the oriL/3CD interaction [16]. The lack of a positive charge at position 156 (in mutants with TGS and TGM tripeptides) could be partially compensated by the appearance of such a charge (e.g., in Arg) at position 153, just preceding the relevant tripeptide. Of note is that Lys153 has almost the same steric potential to interact with RNA-ligands as Lys156, as follows from the comparison of crystal structures of TGK-containing (poliovirus) and KIGQ156-containing (rhinovirus A2) 3C proteins: Lys153 of rhinovirus exposes its positive charge to the same surface area as Lys156 of poliovirus, though this area in the former 3C has a somewhat lower overall positive charge, due to a lesser abundance of basic amino acids [42,46] (Figure 9, compare panels (B) and (C)).
It is not clear whether debilitating effects of certain “poor” residues in the relevant tripeptide were linked to the disappearance of distinct RNA-protein interactions directly involving these residues or to changes in the protein conformation and solubility. In the latter cases, the possibility of dynamic changes of this conformation to modulate its functionality should be considered. It may be worth remembering that the functionally optimal conformation of the tetraloop of domain d of oriL could be provided by different sequences of the YNMG consensus, and it has been proposed that certain non–YNMG sequences are able to temporarily acquire a YNMG-like conformation as a result of molecular dynamics, acquiring thereby some level of functionality [20].

5. Conclusions

This study provides new quantitative and mechanistic insights into how the fundamental conflict between the infidelity of genome replication and the need to retain identity may be solved in RNA viruses. This conflict becomes especially biologically relevant in the case of various bottlenecks, often accompanying viral life history. After such bottlenecks, new viral lineages are often initiated by individual or a small number of mutated genomes haphazardly picked up from always highly heterogeneous (quasispecies) viral populations. By using as the exemplary model an essential viral RNA/protein (oriL/3CD) interaction required for the efficient replication of poliovirus, it was demonstrated that at least 11 amino acid replacements in the involved highly conserved key tripeptide (TGK) of the protein ligand did not kill the virus. Combining the present results with those reported by us previously [20], we concluded that at least 42 out of 51 possible nucleotide replacements within the two relevant genomic regions controlling this interaction (octa- and nonanucleotides of the oriL and 3C gene, respectively) could be tolerated without the loss of viral viability. A significant proportion of these acceptable mutations did not markedly affect the viral phenotype as studied in vitro. The viruses adversely affected by other viability-compatible mutations exhibited a remarkable level of resilience, i.e., the capacity to regain, fully or partially, the fitness by alterations (reversions or compensatory changes) of either RNA or protein ligands. This resilience is based again on the infidelity of viral RNA replication as well as on selection.
The levels of mutational tolerance and robustness described here could be extrapolated (with some obvious limitations related to differences in the replicative mechanisms) to other RNA viruses and thus could contribute to better understanding of their conservation and evolvability.

Author Contributions

A.P.G. and V.I.A. conceived this study; M.A.P., D.V.B., A.P.G., and V.I.A. designed the experiments; M.A.P., E.I.S., D.V.B., A.P.G., A.A.G., E.V.K., M.S.K., and A.A.S. performed the experiments, M.A.P., D.V.B., A.A.S., A.P.G., and V.I.A. analyzed the data, and M.A.P., A.P.G., and V.I.A. wrote the paper.

Funding

This research was supported by the Russian Science Foundation under grant N15-15-00147.

Acknowledgments

We thank Alexander E. Gorbalenya for useful suggestions.

Conflicts of Interest

The authors declare no conflict of interest. The founding sponsors had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, and in the decision to publish the results.

References

  1. Domingo, E.; Holland, J.J. RNA virus mutations and fitness for survival. Annu. Rev. Microbiol. 1997, 51, 151–178. [Google Scholar] [CrossRef]
  2. Drake, J.W.; Holland, J.J. Mutation rates among RNA viruses. Proc. Natl. Acad. Sci. USA 1999, 96, 13910–13913. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Sanjuán, R.; Nebot, M.R.; Chirico, N.; Mansky, L.M.; Belshaw, R. Viral mutation rates. J. Virol. 2010, 84, 9733–9748. [Google Scholar] [CrossRef] [PubMed]
  4. Acevedo, A.; Brodsky, L.; Andino, R. Mutational and fitness landscapes of an RNA virus revealed through population sequencing. Nature 2014, 505, 686–690. [Google Scholar] [CrossRef] [PubMed]
  5. Andino, R.; Domingo, E. Viral quasispecies. Virology 2015, 479, 46–51. [Google Scholar]
  6. Domingo, E.; Perales, C. Quasispecies and virus. Euro. Biophys. J. 2018, 47, 443–457. [Google Scholar] [CrossRef]
  7. Escarmis, C.; Lazaro, E.; Manrubia, S.C. Population bottlenecks in quasispecies dynamics. In Quasispecies: Concept and Implications for Virology; Springer: Berlin/Heidelberg, Germany, 2006; pp. 141–170. [Google Scholar]
  8. Zwart, M.P.; Elena, S.F. Matters of Size: Genetic Bottlenecks in Virus Infection and Their Potential Impact on Evolution. Ann. Rev. Virol. 2015, 2, 161–179. [Google Scholar] [CrossRef]
  9. Agol, V.I.; Gmyl, A.P. Emergency Services of Viral RNAs: Repair and Remodeling. Microbiol. Mol. Biol. Rev. 2018, 82, e00067-17. [Google Scholar] [CrossRef]
  10. Agol, V.I. Picornavirus genome: An overview. In Molecular Biology of Picornaviruses; Semler, B.L., Wimmer, E., Eds.; ASM Press: Washington, DC, USA, 2002; pp. 127–148. [Google Scholar]
  11. Palmenberg, A.; Neubauer, D.; Skern, T. Genome organization and encoded proteins. In The Picornaviruses; Ehrenfeld, E., Domingo, E., Roos, R.P., Eds.; ASM Press: Washington, DC, USA, 2010; pp. 3–17. [Google Scholar]
  12. Rivera, V.M.; Welsh, J.D.; Maizel, J.V. Comparative Sequence-Analysis of the 5′ Noncoding Region of the Enteroviruses and Rhinoviruses. Virology 1988, 165, 42–50. [Google Scholar] [CrossRef]
  13. Andino, R.; Rieckhof, G.E.; Baltimore, D. A Functional Ribonucleoprotein Complex Forms around the 5′ End of Poliovirus RNA. Cell 1990, 63, 369–380. [Google Scholar] [CrossRef] [PubMed]
  14. Barton, D.J.; O’Donnell, B.J.; Flanegan, J.B. 5′ cloverleaf in poliovirus RNA is a cis-acting replication element required for negative-strand synthesis. EMBO J. 2001, 20, 1439–1448. [Google Scholar] [CrossRef] [PubMed]
  15. Lyons, T.; Murray, K.E.; Roberts, A.W.; Barton, D.J. Poliovirus 5′-Terminal Cloverleaf RNA Is Required in cis for VPg Uridylylation and the Initiation of Negative-Strand RNA Synthesis. J. Virol. 2001, 75, 10696–10708. [Google Scholar] [CrossRef] [PubMed]
  16. Vogt, D.A.; Andino, R. An RNA element at the 5′-end of the poliovirus genome functions as a general promoter for RNA synthesis. PLoS Pathogens 2010, 6, e1000936. [Google Scholar] [CrossRef]
  17. Harris, K.S.; Xiang, W.K.; Alexander, L.; Lane, W.S.; Paul, A.V.; Wimmer, E. Interaction of Poliovirus Polypeptide 3CD(pro) with the 5′-Termini and 3′-Termini of the Poliovirus Genome—Identification of Viral and Cellular Cofactors Needed for Efficient Binding. J. Biol. Chem. 1994, 269, 27004–27014. [Google Scholar]
  18. Xiang, W.K.; Harris, K.S.; Alexander, L.; Wimmer, E. Interaction between the 5′-Terminal Cloverleaf and 3AB/3CD(pro) of Poliovirus Is Essential for RNA Replication. J. Virol. 1995, 69, 3658–3667. [Google Scholar] [PubMed]
  19. Parsley, T.B.; Towner, J.S.; Blyn, L.B.; Ehrenfeld, E.; Semler, B.L. Poly (rC) binding protein 2 forms a ternary complex with the 5′-terminal sequences of poliovirus RNA and the viral 3CD proteinase. RNA 1997, 3, 1124–1134. [Google Scholar] [PubMed]
  20. Prostova, M.A.; Gmyl, A.P.; Bakhmutov, D.V.; Shishova, A.A.; Khitrina, E.V.; Kolesnikova, M.S.; Serebryakova, M.V.; Isaeva, O.V.; Agol, V.I. Mutational robustness and resilience of a replicative cis-element of RNA virus: promiscuity, limitations, relevance. RNA Biol. 2015, 12, 1338–1354. [Google Scholar] [CrossRef] [PubMed]
  21. Andino, R.; Rieckhof, G.E.; Trono, D.; Baltimore, D. Substitutions in the Protease (3Cpro) Gene of Poliovirus Can Suppress a Mutation in the 5′ Noncoding Region. J. Virol. 1990, 64, 607–612. [Google Scholar]
  22. Andino, R.; Rieckhof, G.E.; Achacoso, P.L.; Baltimore, D. Poliovirus RNA-Synthesis Utilizes an RNP Complex Formed around the 5′-End of Viral-RNA. EMBO J. 1993, 12, 3587–3598. [Google Scholar] [CrossRef]
  23. Pilipenko, E.V.; Gmyl, A.P.; Maslova, S.V.; Svitkin, Y.V.; Sinyakov, A.N.; Agol, V.I. Prokaryotic-Like Cis Elements in the Cap-Independent Internal Initiation of Translation on Picornavirus RNA. Cell 1992, 68, 119–131. [Google Scholar] [CrossRef]
  24. Humphrey, W.; Dalke, A.; Schulten, K. VMD: Visual molecular dynamics. J. Mol. Graph. 1996, 14, 33–38. [Google Scholar] [CrossRef]
  25. Jurrus, E.; Engel, D.; Star, K.; Monson, K.; Brandi, J.; Felberg, L.E.; Brookes, D.H.; Wilson, L.; Chen, J.; Liles, K.; et al. Improvements to the APBS biomolecular solvation software suite. Protein Sci. 2018, 27, 112–128. [Google Scholar] [CrossRef] [PubMed]
  26. Ward, C.D.; Stokes, M.A.M.; Flanegan, J.B. Direct Measurement of the Poliovirus RNA-Polymerase Error Frequency In vitro. J. Virol. 1988, 62, 558–562. [Google Scholar] [PubMed]
  27. Castro, C.; Arnold, J.J.; Cameron, C.E. Incorporation fidelity of the viral RNA-dependent RNA polymerase: A kinetic, thermodynamic and structural perspective. Virus Res. 2005, 107, 141–149. [Google Scholar] [CrossRef] [PubMed]
  28. Freistadt, M.S.; Vaccaro, J.A.; Eberle, K.E. Biochemical characterization of the fidelity of poliovirus RNA-dependent RNA polymerase. Virol. J. 2007, 4, 44. [Google Scholar] [CrossRef] [PubMed]
  29. Pfeiffer, J.K.; Kirkegaard, K. A single mutation in poliovirus RNA-dependent RNA polymerase confers resistance to mutagenic nucleotide analogs via increased fidelity. Proc. Natl. Acad. Sci. USA 2003, 100, 7289–7294. [Google Scholar] [CrossRef] [Green Version]
  30. Cameron, C.E.; Moustafa, I.M.; Arnold, J.J. Fidelity of Nucleotide Incorporation by the RNA-Dependent RNA Polymerase from Poliovirus. Enzymes 2016, 39, 293–323. [Google Scholar] [CrossRef] [Green Version]
  31. Borderia, A.V.; Rozen-Gagnon, K.; Vignuzzi, M. Fidelity Variants and RNA Quasispecies. Curr. Top. Microbiol. 2016, 392, 303–322. [Google Scholar] [CrossRef]
  32. Yang, X.R.; Liu, X.R.; Musser, D.M.; Moustafa, I.M.; Arnold, J.J.; Cameron, C.E.; Boehr, D.D. Triphosphate Reorientation of the Incoming Nucleotide as a Fidelity Checkpoint in Viral RNA-dependent RNA Polymerases. J. Biol. Chem. 2017, 292, 3810–3826. [Google Scholar] [CrossRef] [Green Version]
  33. Peersen, O.B. Picornaviral polymerase structure, function, and fidelity modulation. Virus Res. 2017, 234, 4–20. [Google Scholar] [CrossRef] [Green Version]
  34. Vignuzzi, M.; Wendt, E.; Andino, R. Engineering attenuated virus vaccines by controlling replication fidelity. Nat. Med. 2008, 14, 154–161. [Google Scholar] [CrossRef]
  35. Meng, T.; Kwang, J. Attenuation of Human Enterovirus 71 High-Replication-Fidelity Variants in AG129 Mice. J. Virol. 2014, 88, 5803–5815. [Google Scholar] [CrossRef] [Green Version]
  36. McDonald, S.; Block, A.; Beaucourt, S.; Moratorio, G.; Vignuzzi, M.; Peersen, O.B. Design of a Genetically Stable High Fidelity Coxsackievirus B3 Polymerase That Attenuates Virus Growth in Vivo. J Biol. Chem. 2016, 291, 13999–14011. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Gamarnik, A.V.; Andino, R. Switch from translation to RNA replication in a positive-stranded RNA virus. Gene Dev. 1998, 12, 2293–2304. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Prostova, M.A.; Deviatkin, A.A.; Tcelykh, I.O.; Lukashev, A.N.; Gmyl, A.P. Independent evolution of tetraloop in enterovirus oriL replicative element and its putative binding partners in protein 3C. PeerJ 2017, 5, e3896. [Google Scholar] [CrossRef] [PubMed]
  39. Rozovics, J.M.; Semler, B.L. Genome replication I: The players. In The Picornaviruses; Ehrenfeld, E., Domingo, E., Roos, R.P., Eds.; ASM Press: Wahington, DC, USA, 2010; pp. 107–125. [Google Scholar]
  40. Kirkegaard, K.; Semler, B.L. Genome replication II: The process. In The Picornaviruses; Ehrenfeld, E., Domingo, E., Roos, R.P., Eds.; ASM Press: Washington, DC, USA, 2010; pp. 127–140. [Google Scholar]
  41. Martinez-Salas, E.; Ryan, M.D. Translation and protein processing. In The Picornaviruses; Ehrenfeld, E., Domingo, E., Roos, R.P., Eds.; ASM Press: Washington, DC, USA, 2010; pp. 141–161. [Google Scholar]
  42. Mosimann, S.C.; Cherney, M.M.; Sia, S.; Plotch, S.; James, M.N.G. Refined x-ray crystallographic structure of the poliovirus 3C gene product. J. Mol. Biol. 1997, 273, 1032–1047. [Google Scholar] [CrossRef] [PubMed]
  43. Dahl, D.B.; Bohsnnan, Z.; Mo, Q.; Vannucci, M.; Tsai, J. Assessing side-chain perturbations of the protein backbone: A knowledge-based classification of residue Ramachandran space. J. Mol. Biol. 2008, 378, 749–758. [Google Scholar] [CrossRef] [PubMed]
  44. Hammerle, T.; Molla, A.; Wimmer, E. Mutational Analysis of the Proposed Fg Loop of Poliovirus Proteinase 3C Identifies Amino-Acids That Are Necessary for 3CD Cleavage and Might Be Determinants of a Function Distinct from Proteolytic Activity. J. Virol. 1992, 66, 6028–6034. [Google Scholar]
  45. Blair, W.S.; Parsley, T.B.; Bogerd, H.A.L.P.; Towner, J.S.; Semler, B.L.; Cullen, B.R. Utilization of a mammalian cell-based RNA binding assay to characterize the RNA binding properties of picornavirus 3C Utilization of a mammalian cell-based RNA binding assay to characterize the RNA binding properties of picornavirus 3C proteinases. RNA 1998, 4, 215–225. [Google Scholar]
  46. Matthews, D.A.; Dragovich, P.S.; Webber, S.E.; Fuhrman, S.A.; Patick, A.K.; Zalman, L.S.; Hendrickson, T.F.; Love, R.A.; Prins, T.J.; Marakovits, J.T.; et al. Structure-assisted design of mechanism-based irreversible inhibitors of human rhinovirus 3C protease with potent antiviral activity against multiple rhinovirus serotypes. Proc. Natl. Acad. Sci. USA 1999, 96, 11000–11007. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Schematic representation of the poliovirus genome, its replicative cis-elements, and formation of the oriL/3CD complex. The positions of only 3C- and 3D-coding sequences are shown in the viral polyprotein reading frame (rectangle). Locations of the major replicative cis-elements (oriL, cre, and oriR) are indicated.
Figure 1. Schematic representation of the poliovirus genome, its replicative cis-elements, and formation of the oriL/3CD complex. The positions of only 3C- and 3D-coding sequences are shown in the viral polyprotein reading frame (rectangle). Locations of the major replicative cis-elements (oriL, cre, and oriR) are indicated.
Viruses 11 00479 g001
Figure 2. Generation of plasmids encoding viral genomes with randomized segments (A) or distinct mutations (B).
Figure 2. Generation of plasmids encoding viral genomes with randomized segments (A) or distinct mutations (B).
Viruses 11 00479 g002
Figure 3. SDS PAGE profiles of the engineered mutant 3CD proteins.
Figure 3. SDS PAGE profiles of the engineered mutant 3CD proteins.
Viruses 11 00479 g003
Figure 4. Plaque phenotypes of the engineered viruses. Fifty ng of transcripts were used for the transfection with wild type genomes in all these experiments. (A,B) Altered tripeptide TGK. Twenty to 50 ng and ~1 µg of mutant transcripts were used for the transfection of a flask in (A) and (B), respectively. (C) Altered tetrapeptide CTGK. Fifty to 150 ng of transcript per flask. Days post-transfection are indicated. * the same flask is shown at days 3 and 6. The results presented in panels AC were obtained in three separate experiments, and some variability of the plaque sizes in the control (TGK-containing) samples could be due to different sensitivities of the cell cultures.
Figure 4. Plaque phenotypes of the engineered viruses. Fifty ng of transcripts were used for the transfection with wild type genomes in all these experiments. (A,B) Altered tripeptide TGK. Twenty to 50 ng and ~1 µg of mutant transcripts were used for the transfection of a flask in (A) and (B), respectively. (C) Altered tetrapeptide CTGK. Fifty to 150 ng of transcript per flask. Days post-transfection are indicated. * the same flask is shown at days 3 and 6. The results presented in panels AC were obtained in three separate experiments, and some variability of the plaque sizes in the control (TGK-containing) samples could be due to different sensitivities of the cell cultures.
Viruses 11 00479 g004
Figure 5. Plaque phenotypes (day 4 post transfection) of viruses in which both the apical loop of domain d of oriL and the TGK motif of 3C were modified.
Figure 5. Plaque phenotypes (day 4 post transfection) of viruses in which both the apical loop of domain d of oriL and the TGK motif of 3C were modified.
Viruses 11 00479 g005
Figure 6. Efficiency of replication of viral genomes with modifications in only the TGK motif of 3CD (AC) or in both oriL’s tetraloop and TGK (D). Monolayers of Vero cells were transfected with the RNAs modified as indicated and the accumulation of viral genomes was determined by quantitative PCR. Each point reflects the results obtained in triplicated samples.
Figure 6. Efficiency of replication of viral genomes with modifications in only the TGK motif of 3CD (AC) or in both oriL’s tetraloop and TGK (D). Monolayers of Vero cells were transfected with the RNAs modified as indicated and the accumulation of viral genomes was determined by quantitative PCR. Each point reflects the results obtained in triplicated samples.
Viruses 11 00479 g006
Figure 7. Electrophoretic mobility shift assay of the affinity between altered 3CD proteins and wild-type oriL. (A) Illustration of an experiment. The presence of 3CD in the complexes with the mobility indicated has been demonstrated previously by Western blotting [20]. (B) The optical densities of the oriL/3CD complex bands (relative to that of the TGK-containing controls) observed in two experiments (the interaction of oriL with the RTGM-containing 3C was investigated only once).
Figure 7. Electrophoretic mobility shift assay of the affinity between altered 3CD proteins and wild-type oriL. (A) Illustration of an experiment. The presence of 3CD in the complexes with the mobility indicated has been demonstrated previously by Western blotting [20]. (B) The optical densities of the oriL/3CD complex bands (relative to that of the TGK-containing controls) observed in two experiments (the interaction of oriL with the RTGM-containing 3C was investigated only once).
Viruses 11 00479 g007
Figure 8. The space of nucleotide (A) and amino acid (B) sequences corresponding to codons 154–156 of the poliovirus protein 3C. The nucleotides and amino acids in the wild-type poliovirus are given in black. The set of viability-compatible nucleotide substitutions caused by single point mutation and substitutions requiring an additional mutation are given in blue and orange, respectively. The viability-compatible amino acids are in green.
Figure 8. The space of nucleotide (A) and amino acid (B) sequences corresponding to codons 154–156 of the poliovirus protein 3C. The nucleotides and amino acids in the wild-type poliovirus are given in black. The set of viability-compatible nucleotide substitutions caused by single point mutation and substitutions requiring an additional mutation are given in blue and orange, respectively. The viability-compatible amino acids are in green.
Viruses 11 00479 g008
Figure 9. Some features of the crystal structure of 3C. Panel (A) (based on PDB ID 1L1N) demonstrates the proximity of the RNA-binding motif TGK to some other relevant amino acids. Oxygen, nitrogen, and carbon atoms are colored in red, blue, and cyan, respectively, and β-barrels and α-helices in yellow and blue/magenta, respectively. The distances between α carbons of His89, Gly155, Tyr6, β carbon of Cys153, ε carbon of Lys156, C4 of Tyr6, and C5 of His89 are indicated in Å. The electrostatic potentials of the solvent accessible surfaces of proteins 3C of poliovirus (based on PDB ID 4DCD) and rhinovirus (based on PDB ID 1CQQ), shown on panels (B) and (C), respectively, were calculated by using a PyMOL APBS Electrostatics Plugin. The potential range from negative −3 mV (blue) to positive +3 mV (red) is shown at the bottom of panel (C). Relevant residues (see the Discussion section) are highlighted.
Figure 9. Some features of the crystal structure of 3C. Panel (A) (based on PDB ID 1L1N) demonstrates the proximity of the RNA-binding motif TGK to some other relevant amino acids. Oxygen, nitrogen, and carbon atoms are colored in red, blue, and cyan, respectively, and β-barrels and α-helices in yellow and blue/magenta, respectively. The distances between α carbons of His89, Gly155, Tyr6, β carbon of Cys153, ε carbon of Lys156, C4 of Tyr6, and C5 of His89 are indicated in Å. The electrostatic potentials of the solvent accessible surfaces of proteins 3C of poliovirus (based on PDB ID 4DCD) and rhinovirus (based on PDB ID 1CQQ), shown on panels (B) and (C), respectively, were calculated by using a PyMOL APBS Electrostatics Plugin. The potential range from negative −3 mV (blue) to positive +3 mV (red) is shown at the bottom of panel (C). Relevant residues (see the Discussion section) are highlighted.
Viruses 11 00479 g009
Table 1. Oligonucleotides used.
Table 1. Oligonucleotides used.
No.NameSequence aPosition in the Genome
13CTGKGTGTGGTGGAGTCATCACATGTnn(g/c)nn(g/c)nn(g/c)GTCATCGGcATGCATGTTGGTGGGAACGGTTCACACG5875–5942
2 3CMluITTCCAACgcGtGCAGGACAGTGTGGTGGAGTCATCAC5856–5892
3 3CSacIICAGGGCcGCGGCAAACCCGTGTGAACCGTTCC5928–5959, complement
4Rib2gaggccgaaaggccgaaaagggcctatgggcccttcTTAAAACAGCTCTGG1–15
5DEN3GAAACAGAAGTGCTTGTTCG158–177, complement
63EP17CATTCTCCCATGTGAC7082–7097, complement
73EP4TTTTTTTTTTTTTTCTCCG7436–polyA, complement
8B5594GAAGTGGAGATCTTGGATGCC5594–5614
9mut-s bCATCACATGT(mut-s)GTCATC5892–5911
10mut-a b(mut-a)ACATGTGATGACTCC5902–5882, complement
11PVP1FAM-TTGATTCATGAATTTCCTTCATTGGCA-BHQ17161–7187 complement
12PVR1CGAACGTGATCCTGAGTGT7212–7230, complement
13PVL1GGCAGACGAGAAATACCCAT7123–7142
a The mutation-generating and nonviral nucleotides are in lower case. FAM, 6-carboxyfluorescein; BHQ1, black hole quencher 1. b mut-s and mut-a stand for various distinct mutagenic nonanucleotides of sense and anti-sense polarity, respectively.
Table 2. The prevalence (%) of different nucleotides in the codons corresponding to amino acids 154–156 of the protein 3C in 20 randomized plasmids a.
Table 2. The prevalence (%) of different nucleotides in the codons corresponding to amino acids 154–156 of the protein 3C in 20 randomized plasmids a.
NucleotideCodons of the 3C-Coding Region
154155156
N1N2N3N1N2N3N1N2N3
G252570104570604090
C20153030353054010
T453002010020150
A10300401001550
a Four plasmids (two with deletions and two with heterogeneity in the randomized region) were omitted from this analysis.
Table 3. The coding capacity of plasmids.
Table 3. The coding capacity of plasmids.
No.Sequence of the Randomized NonanucleotideEncoded TripeptideNon-Intended Alterations a
1ggg ggg ccgGly Gly Pronone
2tgg cgg gtgTrp Arg Valnone
3agg acg gcgArg Thr Alanone
4ctc cgg ggcLeu Arg Glyt5916a (Met160Lys)
5gcg ctg acgAla Leu Thrnone
6gcc acg tcgAla Thr Sernone
7ttg ggc ttgLeu Gly Leua5884g (synonymous)
8tgg aac tggTrp Asn TrpΔ5893, frameshift
9cag tcg gtgGln Ser ValΔ5860–5861, frameshift
10gcg agg acgAla Arg ThrΔ5893–5894, frameshift
11ttg agc gcgLeu Ser Alanone
12tag ctg gagStop Leu Glutermination of translation
13tac tсg gcgTyr Ser AlaΔ5893, frameshift
14gtc acg gggVal Thr GlyΔ5895, frameshift
15tag tcc gggStop Ser Glytermination of translation
16agc ccc aggSer Pro ArgInsertion of g at position 5882, Δ5925, local frameshift
17cac agc tggHis Ser TrpΔ5895, frameshift
18ctg cgg gggLeu Arg GlyΔ5906, frameshift
19ttg tag gcgLeu Stop AlaΔ5870, frameshift
20tag agg ggcStop Arg Glytermination of translation
21aa(c/g) tac (t/a)gg bAsn/Lys Tyr Trp/Argnone
22tg(c/g) (t/g)ac ccg bCys/Trp Tyr/Asp Pronone
23deletion in randomized region cdeletion in randomized regionframeshift
24deletion in randomized region cdeletion in randomized regionframeshift
a In the sequenced region (positions 5800–5960 of the virus RNA encoding amino acids 121–174 of 3C). b Nucleotide heterogeneity in the randomized region. с Deletion in the randomized region, as judged by its unchanged flanking sequences.
Table 4. Specific infectivity of pools with different numbers of variants.
Table 4. Specific infectivity of pools with different numbers of variants.
No.RNA, µgNumber of Variants in the PoolDay of Plaque Appearance and Their Size aNumber of PlaquesSpecific Infectivity, pfu/µg RNA b
11.5313, minute77, 12667.7
221003, small9, 207.25
31.51025, small102, 15084
41.51165, small100, 10066.7
51.51485, minute 20, 39 19.7
61.51855, minute20, 5023.4
722503, minute3, 6 2.25
81.52713, minute103, 13077.7
91.52973, minute89, 9059.7
1023503, small17, 23 10
1135253, small3, 183.5
121.56303, minute84, 130 71.4
131.56663, small34, 8539.7
1438503, small15, 307.5
153~17503, large5, 61.8
163~17503, large1, 20.5
172~22503, large14, 26 10
183~22503, large/minute16, 216.2
193~22503, large/minute43, 72 19.2
a The cultures were observed for plaque appearance at days 3, 5, and 6 p.t. The day at which the plaques were first noticed and their size are indicated. The sizes of large, small, and minute plaques were >3 mm, 1–3 mm, and <1 mm, respectively. b The specific infectivity of the transcript corresponding to the full genome of wild type poliovirus, strain Mahoney, was 7.1 × 104 pfu/µg.
Table 5. The prevalence (%) of different nucleotides in 22 unique sequences of 3C codons 154–156 of RNA of viable viruses.
Table 5. The prevalence (%) of different nucleotides in 22 unique sequences of 3C codons 154–156 of RNA of viable viruses.
NucleotidesCodons of the 3C-Coding Region
154155156
N1N2N3N1N2N3N1N2N3
G359541001008606886
C9274100144159
U136450001490
A430000045185
Table 6. Nucleotide and amino acid sequences in the randomized region of the selected viable viruses and their genetic stability a.
Table 6. Nucleotide and amino acid sequences in the randomized region of the selected viable viruses and their genetic stability a.
No.Primary Sequence bSequence after 5–6 Additional Passages
NucleotidesAmino AcidsNucleotidesAmino Acids
1acg ggg aaaTGKnot donenot done
2acg ggg cggTGRacc ggg cggTGR
3acc ggg cgcTGRnot donenot done
4acc ggg cggTGRnot donenot done
5gug ggg aagVGKgug ggg aagVGK
6guc ggg aagVGKnot donenot done
7guc ggg cggVGRguc ggg cggVGR
8gug ggg aggVGRnot donenot done
9gug ggg cggVGRnot donenot done
10guu ggg aggVGRnot donenot done
11auc ggg aggIGRauc ggg aggIGR
12ugc ggg aagCGKugc ggg aagCGK
13aug gg ccggMGRacg ggc cggTGR
14uug ggc cggLGRgug ggc cggVGR
15cuc ggg aggLGRnot donenot done
16cug ggg cggLGRnot donenot done
17(a/g)ug ggg cgg(М/V)GRgug ggg cggVGR
18gug ggg augVGMgug ggg aagVGK
19ugc ggg ugcCGCugc ggg cgcCGR
20auc ggg uggIGWauc ggg cggIGR
21ugu acg ggc ucg cС153TGScgu acg ggc ucgR153TGS
22cgu acg ggg aug cR153TGM(u/c)gu acg ggg a(u/a)g(C/R)153TG(K/M)
a The nucleotide sequences at positions 5897–5905 of the poliovirus genome and the corresponding amino acid sequences. In the case of changes after passages, the initially determined and changed nucleotides and amino acids are underlined. Parentheses indicate heterogeneous positions. b Initially determined sequences are given for the viruses isolated from the primary plaques (highlighted in bold) and for the viruses subjected to one or more passages to obtain sufficient material for the sequencing. c In addition to the changes in the randomized region, the codon at the preceding positions 5894–5896 (ugu) and the encoded amino acid (С153) were changed upon further passages of the primary isolate 21 and were already changed (and underwent further alterations) in the primary isolate 22.
Table 7. Phenotype and genetic stability of constructed genomes.
Table 7. Phenotype and genetic stability of constructed genomes.
No.Codons (153)154–156 of 3CAmino acid (153)154–156 of 3CDay of Plaque AppearancePlaque SizeRelative Specific Infectivity aChanges in Recovered Viruses b
NucleotideAmino Acids
1acu ggg aaaTGK3large1no changes
2guc ggg aagVGK3large0.8no changes
3auu ggc aaaIGK3large1.2no changes
4ugc ggg aagCGK3large1.4no changes
5acg gg gcggTGR3large1.6not done
6aug ggg aaaMGK3minute0.6a5897gM154V c
a5897uM154L
a5897(a/g)
u5898(u/c)
M154(V/T)
c5699uP88S
7uug ggg aaaLGK6small/heterogeneous0.2u5897gL154V
c5700uP88L
8ucg ggg aaaSGK8solitary plaques d0.0007c5898g
g5899c
S154C
9acg ggg augTGM7–8solitary plaques d0.07u5896cC153R c
u5904aM156K
10acg ggc ucgTGS7–8solitary plaques d0.003c5700uP88L
c5699aP88T
g5738aV101I
g6281aD282N
11(cgu) acg ggg aug(R)TGM3–5small/heterogeneous0.3g5921aV162Ic
a5846gM137V
a5775uY113F
a5625uE63V
g5777aV114I
12(cgu) acg ggc ucg(R)TGS3–5small/heterogeneous0.02c5699UP88Sc
g5798cE121Q
g5738aV101I
a5888cI151L
13acu ggg gaaTGEno plaques
14gug ggg gaaVGEno plaques
15aug ggc cggMGR4–6small1.3a5897gM154Vc
u5898cM154T
u5898(u/c)M154(T/M)
16auc ggg uggIGW6–8small/heterogeneous0.2u5903cW156Rc
17ugc ggg ugcCGC7–8solitary plaques d0.0015u5903cC156R
18gug ggg augVGM7–9solitary plaques d0.0015u5904(u/a)M156(K/M)
19(cgu) acg ggc ucg(R)TGS e3small1.2not done
20(cgu) acg ggc ucg(R)TGS f3small1.4not done
a Specific infectivity of the RNA-transcript of poliovirus strain Mahoney determined in the respective experiments, where (5.25–6.5) × 102 pfu/µg, is assumed to be 1.0. b Positions 5500–6500 of the poliovirus RNA were sequenced. c Such changes were found in several clones. d Late solitary plaques after transfection with a high dose (~1 μg) of the transcript. e Containing the engineered P88S mutation. f Containing the engineered I151L mutation.
Table 8. Partial nucleotide and amino acid sequences and genetic stability of viable viruses selected from the genomes with both randomized apical octanucleotide of domain d and nonanucleotide encoding the TGK motif of 3C.
Table 8. Partial nucleotide and amino acid sequences and genetic stability of viable viruses selected from the genomes with both randomized apical octanucleotide of domain d and nonanucleotide encoding the TGK motif of 3C.
Pool No. aNumber of IsolatesPrimary SequenceSequence after 2–6 Additional Passages
nt. 61–68 of oriLTetraloop ConsensusCodons 154–156 of 3CPassageNt 61–68 of oriLTetraloop ConsensusCodons 154–156 of 3C
NucleotidesAmino acidsNucleotidesAmino Acids
11ugAU(U/G)Uca + insertion G14no consensusgug ggg aggVGR5ugGUN(U/A)ca + insertion G14no consensus/GNUAgug ggg aggVGR
6ugGUUAca + insertion G14GNUAgug ggg aggVGR
21uaGCUCuano consensusno data3uaGCUCuano consensuscug ggg aggLGR
5uaGCU(C/A)uaGNUAgug ggg aggVGR
31uaUCAGugYNMGacg ggg cgcTGR2uaUCAGugYNMGacg ggg cgcTGR
5uaUCAGu(g/a)acg ggg cgcTGR
4 b2uaUCAGugYNMGacg ggg cgcTGR
5uaUCAGugacg ggg cgcTGR
41auUUAUgu + insertion C119no consensusgug ggg aggVGR2guUUAUgc + insertion C119no consensusgug ggg aggVGR
3guUUA(U/G)gc + insertion C119YNMGgug ggg aggVGR
5guUUAGgc + insertion C119gug ggg aggVGR
13auUUAGgu + insertion C119YNMGgug ggg aggVGR
3 b3auUUAUgu + insertion C119no consensusgug ggg aggVGR
a Pools contained transcripts from 1000–15,000 plasmid clones. b All isolates demonstrated the same genetic alterations upon passages.

Share and Cite

MDPI and ACS Style

Prostova, M.A.; Smertina, E.; Bakhmutov, D.V.; Gasparyan, A.A.; Khitrina, E.V.; Kolesnikova, M.S.; Shishova, A.A.; Gmyl, A.P.; Agol, V.I. Characterization of Mutational Tolerance of a Viral RNA–Protein Interaction. Viruses 2019, 11, 479. https://doi.org/10.3390/v11050479

AMA Style

Prostova MA, Smertina E, Bakhmutov DV, Gasparyan AA, Khitrina EV, Kolesnikova MS, Shishova AA, Gmyl AP, Agol VI. Characterization of Mutational Tolerance of a Viral RNA–Protein Interaction. Viruses. 2019; 11(5):479. https://doi.org/10.3390/v11050479

Chicago/Turabian Style

Prostova, Maria A., Elena Smertina, Denis V. Bakhmutov, Anna A. Gasparyan, Elena V. Khitrina, Marina S. Kolesnikova, Anna A. Shishova, Anatoly P. Gmyl, and Vadim I. Agol. 2019. "Characterization of Mutational Tolerance of a Viral RNA–Protein Interaction" Viruses 11, no. 5: 479. https://doi.org/10.3390/v11050479

APA Style

Prostova, M. A., Smertina, E., Bakhmutov, D. V., Gasparyan, A. A., Khitrina, E. V., Kolesnikova, M. S., Shishova, A. A., Gmyl, A. P., & Agol, V. I. (2019). Characterization of Mutational Tolerance of a Viral RNA–Protein Interaction. Viruses, 11(5), 479. https://doi.org/10.3390/v11050479

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop