|Home | About | Journals | Submit | Contact Us | Français|
Coronavirus (CoV) 3C-like proteinase (3CLpro), located in nonstructural protein 5 (nsp5), processes the replicase polyproteins 1a and 1ab (pp1a and pp1ab) at 11 specific sites to produce 12 mature nonstructural proteins (nsp5 to nsp16). Structural and biochemical studies suggest that a conserved Gln residue at the P1 position is absolutely required for efficient cleavage. Here, we investigate the effects of amino acid substitution at the P1 position of 3CLpro cleavage sites of infectious bronchitis virus (IBV) on the cleavage efficiency and viral replication by in vitro cleavage assays and reverse genetic approaches. Our results demonstrated that a P1-Asn substitution at the nsp4-5/Q2779, nsp5-6/Q3086, nsp7-8/Q3462, nsp8-9/Q3672, and nsp9-10/Q3783 sites, a P1-Glu substitution at the nsp8-9/Q3672 site, and a P1-His substitution at the nsp15-16/Q6327 site were tolerated and allowed recovery of infectious mutant viruses, albeit with variable degrees of growth defects. In contrast, a P1-Asn substitution at the nsp6-7/Q3379, nsp12-13/Q4868, nsp13-14/Q5468, and nsp14-15/Q5989 sites, as well as a P1-Pro substitution at the nsp15-16/Q6327 site, abolished 3CLpro-mediated cleavage at the corresponding position and blocked the recovery of infectious viruses. Analysis of the effects of these lethal mutations on RNA synthesis suggested that processing intermediates, such as the nsp6-7, nsp12-13, nsp13-14, nsp14-15, and nsp15-16 precursors, may function in negative-stranded genomic RNA replication, whereas mature proteins may be required for subgenomic RNA (sgRNA) transcription. More interestingly, a mutant 3CLpro with either a P166S or P166L mutation was selected when an IBV infectious cDNA clone carrying the Q6327N mutation at the nsp15-16 site was introduced into cells. Either of the two mutations was proved to enhance significantly the 3CLpro-mediated cleavage efficiency at the nsp15-16 site with a P1-Asn substitution and compensate for the detrimental effects on recovery of infectious virus.
Coronavirus (CoV)-encoded 3C-like proteinase (3CLpro), together with one or two papain-like proteinases (PLpro), mediates the extensive proteolytic processing of two large replicase polyproteins (pp), pp1a and pp1ab, and yields multiple mature nonstructural proteins (nsp) essential for the assembly and function of the viral replication complex. PLpro cleaves the N-terminal regions of the polyproteins at two or three sites, and 3CLpro is responsible for the processing of all downstream parts of the replicase polyproteins at 11 conserved cleavage sites (50, 55, 72, 74). Like its homologs in other coronaviruses, 3CLpro of avian infectious bronchitis virus (IBV) is encoded by open reading frame 1a (ORF1a) and resides in nsp5. This proteinase specifically cleaves polyproteins 1a and 1ab at 11 sites to produce 12 mature products (nsp5 to nsp16) (24, 36-39, 41-43, 64) (Fig. (Fig.11).
The structures of many nonstructural proteins have been determined, and their functions are partially elucidated. nsp7 and nsp8 form a hexadecameric structure able to encircle and bind RNA (70); nsp8 alone is described as a second RdRp of severe acute respiratory syndrome CoV (SARS-CoV) (28). nsp9 is a single-stranded RNA-binding protein that may stabilize nascent and template RNA strands during replication and transcription and may also be involved in RNA processing (9, 16, 53). nsp10 forms a dodecamer structure, and two zinc fingers have been identified in the monomer, implying that it may function in the RNA synthesis machinery (32, 52). More recently, murine hepatitis virus (MHV) nsp10 has been shown to be a critical regulator of viral RNA synthesis and to play an important role in polyprotein processing (13, 14). nsp14 is a bifunctional protein: a 3′-to-5′ exonuclease (ExoN) (40) with a role in maintaining the fidelity of RNA transcription (15) and a cap N7 methyltransferase (10) involved in RNA cap formation. nsp15 is a poly(U)-specific endoribonuclease (NendoU) (6, 8, 25, 33). nsp16, an S-adenosylmethionine-dependent ribo-2′-O-methyltransferase (2′O-MTase), has also been thought to be involved in cap formation (11, 50, 60).
Because of its pivotal role in the viral life cycle and the conservation in structure, 3CLpro becomes the most promising target for drug design against coronaviruses (23, 63, 67, 69). The substrate specificity of coronavirus 3CLpro is determined mainly by the P1, P2, and P1′ positions. Comparative analyses of 77 3CLpro cleavage sites from seven coronaviruses suggest amino acid conservation at the P2 (Leu, Val, Ile, Phe, and Met), P1 (Gln), and P1′ (Ser, Ala, Gly, Asn, and Cys) positions. Among these, the P1 position is exclusively occupied by Gln, and the bulky hydrophobic residues (mainly Leu) are dominant at the P2 position (22, 26, 56, 74). Insights into the catalytic activity and substrate specificity of 3CLpros obtained by the recently determined crystal structures of five 3CLpros from transmissible gastroenteritis virus (TGEV), human coronavirus 229E (HCoV-229E), SARS-CoV, IBV, and HCoV HKU1 further confirm the necessity of a conserved Gln residue at the P1 position (2, 3, 66, 68, 71).
The importance of coronavirus 3CLpro-mediated processing of replicase polyproteins in viral replication was verified by proteinase inhibitors (35, 66, 67), but data on the necessity of coronavirus 3CLpro-mediated individual cleavages are very limited. Reverse genetics studies indicated recently that proteolytic processing at the nsp7-8, nsp8-9, and nsp10-11/12 cleavage sites of MHV was essential for viral replication, while cleavage at the nsp9-10 site of MHV and at the nsp10-11/12 site of IBV was dispensable for viral replication (12, 20). In this report, we investigate the effects of amino acid substitutions at the P1 position of 10 IBV 3CLpro cleavage sites, except the nsp10-11/12 site, on the cleavage efficiency and viral replication by in vitro cleavage assays and IBV reverse genetics. Our results demonstrated that mutation of Gln to Asn at the P1 position of the nsp4-5/Q2779, nsp5-6/Q3086, nsp7-8/Q3462, nsp8-9/Q3672, and nsp9-10/Q3783 sites, P1-Glu substitution at the nsp8-9/Q3672 site, and P1-His substitution at the nsp15-16/Q6327 site were tolerated, since each of these mutations allowed the recovery of infectious virus, albeit with various degrees of defects in viral growth. Our data also showed that proteolytic processing at the nsp6-7/Q3379, nsp12-13/Q4868, nsp13-14/Q5468, nsp14-15/Q5989, and nsp15-16/Q6327 cleavage sites was required for IBV replication. In addition, a distal mutation, P166S or P166L, in 3CLpro was identified when a point mutation, nsp15-16/Q6327N, representing a P1-Asn substitution at the nsp15-16 site, was introduced into the IBV infectious cDNA clone. Either of the two mutations was proved to enhance significantly the 3CLpro-mediated cleavage efficiency at the nsp15-16 site with a P1-Asn substitution and compensate for the detrimental effects on recovery of infectious virus. This study provides further insights into the substrate specificity of coronavirus 3CLpro and the regulatory mechanism of viral RNA synthesis as well as helpful information on the drug design against coronaviruses.
Vero cells and H1299 cells were cultured at 37°C, respectively, in minimal essential medium (MEM) and RPMI supplemented with 10% fetal bovine serum (FBS), penicillin (100 units/ml), and streptomycin (100 μg/ml). A recombinant IBV (rIBV) generated from an infectious clone (19) was used as the wild-type control. All mutant IBVs were propagated in Vero cells in FBS-free MEM.
ORFs placed under the control of a T7 promoter were expressed transiently in mammalian cells using a vaccinia virus-T7 expression system (21). Briefly, confluent H1299 cells in 6-well plates were inoculated with a recombinant vaccinia virus (vTF7-3) expressing T7 RNA polymerase. At 1 h postinfection, medium was removed and plasmids carrying wild-type or mutant cleavage sites of IBV 3CLpro together with pIBV3C containing the full-length IBV 3CLpro sequence (39) were transfected into cells with Lipofectamine 2000 transfection reagent according to the manufacturer's instructions (Invitrogen). Cells were harvested after incubation at 37°C for 16 to 20 h.
Transfected cells were lysed with 2× SDS loading buffer containing 100 mM dithiothreitol (DTT) and 5 mM iodoacetamide, boiled at 100°C for 5 min, and clarified. The proteins were separated by 15% SDS-PAGE and transferred to a polyvinylidene difluoride (PVDF) membrane (Stratagene). The membrane was blocked overnight at 4°C or for 2 h at room temperature in blocking buffer (5% fat-free milk powder in phosphate-buffered saline [PBS] buffer containing 0.1% Tween 20 [PBST]) and then was incubated with diluted primary antibodies in blocking buffer for 2 h at room temperature. After the membrane was washed three times with PBST, it was incubated with 1:2,000 diluted anti-mouse or anti-rabbit IgG antibodies conjugated with horseradish peroxidase (Dako) in blocking buffer for 1 h at room temperature. After the membrane was washed three times with PBST, the polypeptides were detected with a chemiluminescence detection kit (ECL kit; Amersham Biosciences) according to the manufacturer's instructions. The films were exposed and developed.
The DNA fragment containing the nsp4-5/Q2779, nsp5-6/Q3086, nsp6-7/Q3379, nsp7-8/Q3462, nsp8-9/Q3672, nsp9-10/Q3783, and nsp12-13/Q4868 cleavage sites was mapped to plasmid pXL-IBVC, and the fragment containing the nsp13-14/Q5468, nsp14-15/Q5989, and nsp15-16/Q6327 cleavage sites to pGEM-IBVD (19, 54). Constructs containing a single Gln-to-Asn/Glu/Pro/His substitution at the P1 position were produced by using a QuikChange site-directed mutagenesis kit (Stratagene). The full-length infectious clone was assembled as previously described (19, 54) by replacing the corresponding fragment with the mutant fragment. The mutations were verified by automated nucleotide sequencing.
Full-length transcripts were generated in vitro using the mMessage mMachine T7 kit (Ambion, Austin, TX) and electroporated into Vero cells with one pulse at 450 V and 50 μF with a Bio-Rad Gene Pulser II electroporator. The transfected Vero cells were cultured overnight in 1% FBS-containing MEM and further cultured in MEM without FBS (19, 54). The transfected cells were monitored daily for formation of cytopathic effect (CPE). Recovered viruses were plaque purified and passaged on Vero cells. For the mutants that failed to produce observable CPE, three independent experiments were conducted.
Vero cells were infected with rIBV and recovered mutant viruses at a multiplicity of infection (MOI) of ~1 PFU/cell, and 10 μg of total RNA extracted from the infected cells was added to a mixture of 1× MOPS, 37% formaldehyde, and formamide and incubated at 65°C for 20 min prior to gel electrophoresis. The separated RNA bands were transferred onto a Hybond N+ membrane (Amersham Biosciences) via capillary action overnight and fixed by UV cross-linking (Stratalinker). Hybridization of digoxigenin (Dig)-labeled DNA probes corresponding to the IBV 3′ untranslated region (nucleotides [nt] 27104 to 27510) was carried out at 50°C in a hybridization oven overnight. Membranes were washed three times for 15 min each with the stringency buffers, and signals were detected with CDPStar (Roche) according to the manufacturer's instructions.
A confluent monolayer of Vero cells on six-well plates was infected with recovered viruses at the indicated MOI and harvested at different times postinfection. Viral stocks were prepared by freeze-thawing the cells three times. The 50% tissue culture infection dose (TCID50) of each sample was determined by infecting five wells of Vero cells on 96-well plates with a 10-fold serial dilution of each viral stock and using Reed-Muench calculation.
Total RNAs were prepared from the electroporated or infected Vero cells with TRI reagent (Molecular Research Center, Inc.). For the mutations that did not produce the infectious virus, the in vitro transcripts prior to electroporation and the extracted RNAs were treated with DNase I and subjected to phenol-chloroform extraction and ethanol precipitation to completely rule out possible DNA contamination. Reverse transcription (RT) was performed with Expand reverse transcriptase (Roche) according to the manufacturer's instructions using the sense primer IBV leader (5′-26CTATTACACTAGCCTTGCGCT46-3′) for the detection of negative-stranded subgenomic RNA (sgRNA) and the antisense primer IBV24803-R (5′-24803CTCTGGATCC AATAACCTAC24784-3′) for the detection of positive-stranded sgRNA. The two primers were then used for PCR. If transcription of subgenomic mRNAs did occur, a 415-bp PCR product corresponding to the 5′-terminal region of subgenomic mRNA 4 and a 1,010-bp fragment corresponding to the 5′-terminal region of subgenomic mRNA 3 would be expected. Similarly, RT was carried out with the sense primer IBV14931-F (5′-14931GCTTATCCACTAGTACATC14949-3′) for the detection of negative-stranded genomic RNA, and IBV14931-F and the antisense primer IBV15600-R (5′-15600CTTCTCGCACTTCTGCACTAGCA15578-3′) were used for PCR. If replication of viral RNA occurred, a 670-bp PCR fragment would be expected.
Plasmid pXL-IBV-C (nt 8689 to 15532) or pGEM-IBV-D (nt 15511 to 20930) (19) was used to construct a series of T7-driven expression plasmids. A DNA fragment covering each of 10 wild-type cleavage sites on IBV polyproteins was amplified with Turbo Pfu DNA polymerase (Stratagene), and the purified PstI-digested PCR product was inserted in frame into EcoRV/PstI-digested vector pFlag. The resulting plasmid allowed the expression of an N-terminal Flag-tagged protein under the control of the T7 promoter. As templates, these plasmids were used to make the corresponding mutant plasmids containing different single mutations in place of the conserved Gln residue at the P1 position of every cleavage site by using the QuikChange site-directed mutagenesis kit (Stratagene). Plasmids and the primers used for PCR are listed in Table Table11 .
A DNA fragment encoding full-length 3CLpro of the SARS-CoV Sin2774 strain (GenBank accession no. AY283798) was amplified by PCR with primers 5′-CGGGATCCACCATG9969AGTGGTTTTAGGAAAATGGCA-3′ and 5′-CGGAATTCTTA10886TTGGAAGGTAACACCAGAGCA-3′. PCR products were digested with BamHI and EcoRI and cloned into BglII- and EcoRI-digested vector pKTO to get plasmid pSARS-3C.
To examine whether IBV 3CLpro can tolerate other amino acids instead of a Q residue at the P1 position, a set of expression plasmids, containing IBV sequences flanking either wild-type or mutant cleavage sites with mutation of P1-Q to N or E, was constructed (Fig. (Fig.22 A). These plasmids and the corresponding primers used are listed in Table Table1.1. To facilitate detection of the full-length and cleavage products, a Flag tag was fused to the N termini of the expressed IBV peptides. The calculated molecular masses of the Flag-tagged uncleaved (U) and cleaved (C) products are illustrated in Fig. Fig.2A.2A. Wild-type and mutant substrates were coexpressed with IBV 3CLpro, and processing was analyzed by Western blotting with an anti-Flag antibody. The effect of the mutations on the cleavage efficiency was evaluated by measuring the relative density of the U and C bands using the QUANTITY ONE program (Bio-Rad) and calculating the relative cleavage efficiency (RCE) [RCE = C/(U + C) × 100%]. In this expression system, IBV 3CLpro effectively cleaved the substrates at all 10 wild-type cleavage sites, albeit with different cleavage efficiencies (Fig. 2B to D). Notably, nearly complete cleavage at the nsp8-9 site, a so-called noncanonical site with an N residue instead of S, A, or G at the P1′ position, was observed (Fig. (Fig.2C).2C). In contrast, its counterparts in HCoV-229E, TGEV, MHV, and SARS-CoV were found to be hydrolyzed with far less efficiency (18, 26). The reasons for this discrepancy remain unclear but are likely due to different experimental approaches and substrates used.
Compared to wild-type sites, mutation of P1-Q to N showed variable effects on the cleavage efficiency, depending on the location of the site. Complete abolishment of cleavage was observed for the nsp6-7/Q3379N, nsp12-13/Q4868N, nsp13-14/Q5468N, and nsp14-15/Q5989N P1-N mutants, while no inhibitory effect for the nsp4-5/Q2279N mutant was observed (Fig. 2B to D). A minor or moderate reduction in the cleavage efficiency for the nsp7-8/Q3462N, nsp5-6/Q3086N, and nsp8-9/Q3672N mutants (from 9% to 31%), as well as a more dramatic decrease for the nsp9-10/Q3783N and nsp15-16/Q6327N mutants (58% and 69%, respectively), was detected (Fig. 2B to D). Mutation of Q to E showed effects similar to those of the Q-to-N substitution in most cases, except at the nsp13-14 and nsp15-16 sites. Cleavage at the nsp13-14/Q5468E site showed low efficiency (17%), but no cleavage was observed at the nsp13-14/Q5468N site (0%) (Fig. (Fig.2D).2D). To the contrary, higher cleavage efficiency was observed at the nsp15-16/Q6327E site (36%) than at the nsp15-16/Q6327N site (5%) (Fig. (Fig.2D).2D). Replacement of the P1-Q residues with H and P at the nsp15-16 site was also analyzed. Mutation of P1-Q to H at the nsp15-16 site maintained a relatively high level of cleavage efficiency (44%), but no cleavage was observed when P1-Q was replaced by P (0%) (Fig. (Fig.2D).2D). Taken together, these findings show that IBV 3CLpro was able to tolerate mutation of Q to other amino acids, such as N, E, or H, at several sites. The requirement for a conserved P1-Q residue at cleavage sites in the 1b region was more stringent than at those in the 1a region, except the nsp6-7 site.
In previous in vitro cleavage assays, synthetic peptides flanking the N-terminal autocleavage site of a certain coronavirus 3CLpro could be cleaved by 3CLpros from all three groups of coronaviruses (3, 26, 67), probably due to the structural similarity among these proteinases (3, 68). To test if this low-stringency requirement for enzyme substrate specificity could be extended to all other sites, cleavage of substrates with either wild-type or mutant sites by a heterogeneous 3CLpro from SARS-CoV was analyzed. Cleavage with various efficiencies (from 45% to 100%) was observed at the nsp4-5, nsp5-6, nsp7-8, nsp8-9, nsp9-10, and nsp15-16 wild-type sites, whereas no detectable cleavage was observed at other sites (Fig. 2B to D). Moreover, SARS-CoV 3CLpro could also tolerate the mutation of P1-Q to N or E at the two sites (nsp4-5/Q2279 and nsp5-6/Q3086) flanking IBV 3CLpro, although the processing efficiency of the mutant N-terminal site was far higher than that of the mutant C-terminal site (Fig. (Fig.2B),2B), in contrast with the previous results showing that cleavage was completely abolished when the P1-N or P1-E substitution was introduced into the N-terminal self-cleavage site of SARS-CoV (17). In addition, a low level (8%) of cleavage at the nsp15-16 site with a P1-H substitution was also detected (Fig. (Fig.2D2D).
To evaluate the impact of Q-N point mutations at the P1 positions on virus replication, full-length cDNA clones containing each of the mutant cleavage sites were assembled. In vitro transcripts were synthesized and electroporated into Vero cells. For mutations that allowed efficient proteolytic processing, including nsp4-5/Q2779N, nsp5-6/Q3086N, nsp7-8/Q3462N, nsp8-9/Q3672N, and nsp9-10/Q3783N, the corresponding mutant virus was recovered at 2 or 3 days postelectroporation. In contrast, mutations which completely abolished the cleavage activity, including nsp6-7/Q3379N, nsp12-13/Q4868N, nsp13-14/Q5468N, and nsp14-15/Q5989N, blocked the generation of infectious viruses. In cells electroporated with nsp15-16/Q6327N mutant transcripts, formation of several small syncytial cells was observed at 3 days postelectroporation. These syncytia disappeared after prolonged incubation. Following 6 rounds of successive reinfection of fresh Vero cells with freeze-thawed cell preparations, sgRNA was still detectable by RT-PCR, but no syncytium formation was observed (data not shown).
As cleavage was also observed at sites between nsp4 and nsp5, nsp5 and nsp6, nsp8 and nsp9, nsp9 and nsp10, and nsp15 and nsp16 when Q was replaced with E, this point mutation at the site between nsp8 and nsp9 (nsp8-9/Q3672E) was introduced into the infectious IBV clone, and infectious virus was recovered from cells electroporated with the mutant transcripts.
The five recovered mutant viruses, the nsp4-5/Q2779N, nsp5-6/Q3086N, nsp7-8/Q3462N, nsp8-9/Q3672N, and nsp9-10/Q3783N mutants, were plaque purified, and their growth properties were characterized by TCID50 analysis of passage 2 of each mutant virus. The results showed that the nsp7-8/Q3462N, nsp5-6/Q3086N, and nsp9-10/Q3783N mutants exhibited reduced growth phenotypes, decreasing about 2, 1, and 0.5 log10 of the peak TCID50 value at 16 h postinfection compared to wild-type rIBV, while the nsp4-5/Q2779N and nsp8-9/Q3672N mutants displayed delayed and slightly reduced growth phenotypes (Fig. (Fig.33 A). The TCID50 values of these two viruses reached the peak at 24 h postinfection, which was 8 h later than that for rIBV (Fig. (Fig.3A).3A). Analysis of the growth properties of the nsp8-9/Q3672E mutant virus by a TCID50 assay showed that the virus exhibited a delayed and reduced growth phenotype compared to rIBV, reaching the peak titer at 24 h postinfection, which was 8 h later than for the wild-type virus (Fig. (Fig.3A).3A). Moreover, the peak titer was approximately 10-fold less than that for rIBV (Fig. (Fig.3A).3A). Compared to the nsp8-9/Q3672N mutant virus, the nsp8-9/Q3672E mutant virus showed similar growth kinetics, but its peak titer was approximately 5-fold lower (Fig. (Fig.3A3A).
Northern blot analysis was performed to analyze the total RNA extracted from cells infected with the wild type and the recovered mutant viruses. The results showed that most mutant viruses displayed patterns similar to those of the wild-type virus, except that an additional, relatively strong band migrating below sgRNA6 and a relatively reduced sgRNA3 were simultaneously detected in cells infected with the Q3672E mutant virus (Fig. (Fig.3B).3B). As the complete sequence of this mutant virus was not determined, it was not clear at the moment if the detection of these aberrant bands was due to the presence of additional mutations in the genome. Further studies are required to clarify this issue.
The genetic stability of these mutant viruses was investigated by successive passage of plaque-purified mutant viruses on fresh Vero cells up to at least passage 5. Total RNA was extracted from infected cells, and the specific regions flanking cleavage sites were amplified by RT-PCR. Sequence analysis of the RT-PCR products confirmed the presence of the introduced mutations in the nsp4-5/Q2779N, nsp5-6/Q3086N, nsp8-9/Q3672N, nsp8-9/Q3672E, and nsp9-10/Q3783N mutant viruses, suggesting that these mutant viruses were stable in cultured cells. No amino acid change in the nsp5 coding region was identified by sequencing the entire region. Since the recovered mutant viruses were not fully sequenced, additional mutations in other regions could not be ruled out. Nevertheless, these results demonstrated that no other mutations in 3CLpro may be involved in the recovery of these mutant viruses. Interestingly, an increasing trend of reverting to the wild-type sequence was observed during passage of the nsp7-8/Q3462N mutant virus. This reversion (AAC→CAA) was initially detected in passage 3 despite the predominance of the mutant AAC codon (Fig. (Fig.3C).3C). Subsequently, the wild-type codon (CAA) became more dominant (Fig. (Fig.3C).3C). The RT-PCR products from passage 5 were cloned and sequenced. Eight of ten clones contained the wild-type codon (CAA), while the other two clones contained the mutant codon (AAC). Similarly, no mutations were found in 3CLpro.
The lethality of the nsp6-7/Q3379N, nsp12-13/Q4868N, nsp13-14/Q5468N, and nsp14-15/Q5989N substitutions to the recovery of infectious virus may be due to the detrimental effects of these mutations on RNA synthesis. To examine this possibility, RT-PCR was carried out to amplify the genomic and subgenomic RNAs. Vero cells were electroporated with either wild-type or mutant RNA transcripts and cultured at 37°C for 48 h. Total RNA was extracted from the electroporated cells for reverse transcription. RT-PCR products from nucleotides 14931 to 15600, representing both positive and negative strands of the genomic RNA, were detected from cells electroporated with all four full-length mutant transcripts (Fig. (Fig.44 A, lanes 4, 5, 7, 8, 10, 11, 13, and 14). Analysis of the same RNA extracts for the detection of subgenomic RNA species 3 and 4 (sgRNA3 and 4) by RT-PCR showed no detectable RT-PCR products (Fig. (Fig.4B,4B, lanes 3 to 10). These results suggested that mutations at these positions allowed the replication of the negative-stranded genomic RNA but disrupted the transcription of the subgenomic RNA.
Interestingly, substitution of nsp15-16/Q6327 with an N significantly decreased the cleavage efficiency but allowed the synthesis of both negative genomic RNA and sgRNA in cells electroporated with the nsp15-16/Q6327N mutant transcripts (Fig. (Fig.4A,4A, lanes 16 and 17; Fig. Fig.4B,4B, lanes 11 and 12). These results suggest that minimal levels of viral RNA replication and transcription occurred in these cells and prompt more detailed characterization of the site.
As sgRNA was still detected after six rounds of successive incubation of fresh cells with the nsp15-16/Q6327N mutant, but no apparent syncytium formation was observed, the possibility of additional mutations occurring in the spike protein was checked. Total RNA was extracted from the infected cells, and the RT-PCR products covering the full-length spike protein gene were obtained and sequenced. Sequence comparison illustrated that no amino acid change was found. Also, no additional mutation was detected in sequences flanking nsp15-16, except the introduced mutation. Western blotting showed no detectable expression of the structural proteins (data not shown). These data support the argument that only a minimal amount of infectious virus was produced in cells infected with this mutant virus.
To investigate whether this phenotype is associated with a nucleotide sequence across this site that may affect viral RNA packaging, a mutant construct containing six silent mutations at nucleotides 19499 (A→G), 19502 (G→A), 19505 (T→C), 19508 (A→G), 19511 (A→G), and 19514 (A→T) was made and introduced into the infectious clone. Infectious virus carrying six silent mutations was recovered, suggesting that this region is not part of the IBV packaging signal.
Mutation of P1-Q to H allowed the production of viable virus (nsp15-16/Q6327H) with a delayed growth phenotype (Fig. (Fig.55 A). Sequencing the RT-PCR products from infected Vero cells confirmed the presence of this mutation. It was reported that a putative cleavage site in which the P1 position is occupied by the H residue instead of the regular Q was predicted in human coronavirus NL63 (HCoV-NL63) (58), HCoV-HKU1 (61), and bat coronavirus HKU5 (62) based on sequence alignment. In addition, Goetz et al. reported that IBV, MHV, and SARS-CoV 3CLpros exhibit a strong preference for P1-H-containing substrates and both H and Q can adopt the same binding mode in the S1 pocket (23). Furthermore, a recently resolved structure of HCoV-HKU1 and enzyme activity assays confirmed the existence of such a cleavage event in vitro (71). This result provided the first biological evidence that IBV 3CLpro can recognize the P1-H-containing nsp15-16 site.
Replacement of nsp15-16/Q6327 with P abolished completely proteolytic processing as well as virus recovery. Analysis of RNA synthesis in cells electroporated with the mutant transcripts showed that nsp15-16/Q6327P allowed the synthesis of negative-strand genomic RNA but disrupted the sgRNA transcription (Fig. (Fig.5B),5B), suggesting that proteolytic processing at the nsp15-16 site is essential for viral replication.
Coincidently, in a repeated experiment using double amounts of the in vitro synthesized full-length nsp15-16/Q6327N transcripts for electroporation, mutant virus was recovered. After propagation of the recovered virus for 9 passages, five RT-PCR products covering the entire IBV genome were obtained and directly sequenced, showing mixed peaks at two positions (C-to-T changes at nucleotide positions 9361 and 9362) (Fig. (Fig.6A),6A), besides the presence of the introduced Q6327N mutation. To trace the source of the mutations, the corresponding regions of the early virus passages P1 and P3 were sequenced, revealing the presence of the two mutations (Fig. (Fig.6A).6A). In contrast, only the wild-type sequence was identified in passage 6 of the Q6327N virus which did not produce CPE (Fig. (Fig.6A,6A, top). Similarly, only wild-type sequences were detected by sequencing the RT-PCR products from in vitro transcripts prior to electroporation, ruling out the possibility that the mutations were introduced during in vitro transcription. To further confirm the mutations, 11 plaque-purified isolates were randomly selected from passage 9 for sequence analysis, showing that 7 out of 11 isolates carried C9362T, while the other four isolates possessed C9361T. These mutations result in amino acid changes from P to L and from P to S, respectively, at amino acid position 166 of 3CLpro. No wild-type sequences were identified in the same region.
It was speculated that P166S and P166L mutations in 3CLpro may compensate for the detrimental effects on viral replication caused by the nsp15-16/Q6327N mutation. To confirm this possibility, P166S, P166L, P166S/Q6327N, and P166L/Q6327N were introduced individually into the IBV genome. As expected, viable mutant viruses were recovered from all four constructs, and TCID50 assays showed that all four mutant viruses (P166L, P166S, P166L/Q6327N, and P166S/Q6327N) exhibited a slightly reduced growth phenotype, decreasing about 0.4 to 1.0 log unit of the peak TCID50 value at 16 h postinfection, compared to rIBV (Fig. (Fig.6B6B).
To determine whether 3CLpro carrying the P166S or P166L mutation affects cleavage at the mutant nsp15-16 site, processing of wild-type and mutant substrates by wild-type or mutant 3CLpro was tested. Two expression plasmids containing P166S or P166L in 3CLpro were constructed. Upon cotransfection of wild-type substrate with wild-type 3CLpro and each of the mutant 3CLpros, nearly complete cleavage was observed (Fig. (Fig.6C).6C). Upon cotransfection of nsp15-16/Q6327N mutant substrate with wild-type 3CLpro, a trace amount of the cleaved product was detected besides the uncleaved peptide (Fig. (Fig.6C).6C). In contrast, much more efficient cleavage of the mutant substrate was observed when either P166S or P166L mutant 3CLpro was coexpressed (Fig. (Fig.6C),6C), confirming that the mutant 3CLpros could strongly enhance the cleavage efficiency of the mutant substrates flanking the nsp15-16 site.
Attempts to recover mutant viruses carrying double mutations P166S/Q5468N, P166L/Q5468N, P166S/Q5989N, and P166L/Q5989N failed. It was also found that neither 3CLpro with P166S nor 3CLpro with P166L was able to process P1-N-containing substrates flanking the nsp13-14 and nsp14-15 sites, suggesting that this enhancement effect is specific for the P1-N-containing nsp15-16 site.
Among nidoviral 3CLpros, arterivirus and torovirus 3CLpros are serine proteinases, while ronivirus and coronavirus 3CLpros are cysteine proteinases (49, 51, 73, 74). With respect to substrate specificity, arterivirus and ronivirus 3CLpros prefer the Glu residue at the P1 position of the substrates, while coronavirus and torovirus 3CLpros require the Gln residue (49, 51, 73, 74). This differential preference at the P1 position has been thought to be determined by the S1 substrate-binding pocket (2, 3, 5, 66, 68). Among all identified 3CLpro cleavage sites of arterivirus, one site contains P1-Gln, which can also tolerate Gln-to-Asn/Glu substitution at the P1 position (59). The catalytic site Cys-to-Ser substitution in both IBV and SARS-CoV 3CLpro yields an active enzyme (27, 37, 57); conversely, Ser-to-Cys substitution in equine torovirus retains the enzymatic activity (49). In this report, we show that replacement of P1-Gln with Asn at the nsp4-5, nsp5-6, nsp7-8, nsp8-9, and nsp9-10 sites, Glu at the nsp8-9 site, and His at the nsp15-16 site of IBV allowed the recovery of the corresponding infectious mutant virus and the 3CLpro-mediated cleavage of these mutant substrates, suggesting that IBV 3CLpro is able to accept Asn, Glu, and His instead of the conserved Gln residue at the P1 position of several cleavage sites (Table (Table2).2). Further confirmation of the effect of these mutations on cleavage at the individual mutant sites would be required for cells infected with the corresponding recovered mutant viruses. However, due to the lack of specific antibodies, these studies cannot be performed. On the other hand, although the nsp4-5/Q2779N, nsp5-6/Q3086N, nsp7-8/Q3462N, nsp8-9/Q3672N, nsp8-9/Q3672E, and nsp9-10/Q3783N mutant viruses were recovered, all of them displayed reduced and/or delayed growth to a certain extent. Moreover, a reversion of Asn to Gln was observed in the nsp7-8/Q3462N mutant virus, indicating that IBV 3CLpro has a strong preference for the Gln residue at the P1 position, though it can tolerate the P1-Asn or -Glu substitution in these sites. The flexibility for substrate specificity displayed by IBV and arterivirus 3CLpros may reflect the evolutionary relationships between virus-encoded 3C-like serine and 3C-like cysteine proteinases, implying that the two proteinases may derive from a common ancestor.
Complete abolishment of proteolysis and failure to recover infectious virus caused by mutations nsp6-7/Q3379N, nsp12-13/Q4868N, nsp13-14/Q5468N, nsp14-15/Q5989N, and nsp15-16/Q6327P showed that cleavage at the nsp6-7, nsp12-13, nsp13-14, nsp14-15, and nsp15-16 sites is essential for IBV replication. nsp6 contains multiple transmembrane regions and is believed to be involved in the assembly of the replicase complex (4, 44). nsp7 and nsp8 form a hexadecameric structure with RNA binding activity (70). It is possible that the unprocessed nsp6-7 intermediate may affect the formation of a functional replicase complex as well as the nsp7-nsp8 supercomplex and thus impede the initiation of sgRNA synthesis.
Considering the significance and requirement of the processed products of the ORF1b-encoded polyprotein, including RdRp, helicase, and three highly conserved RNA-processing enzymes (NendoU, ExoN, and 2′-O-MTase) for virus replication (1, 7, 10, 11, 15, 29-31, 34, 40, 45, 55, 56), the lethality caused by these mutations may be due to the loss of related enzymatic and other accessory activities. It is likely that the nsp12-nsp13, nsp13-nsp14, nsp14-nsp15, or nsp15-nsp16 fusion protein results in the change in the conformation or structure of proteins, or disruption of the protein-protein interactions associated with sgRNA transcription. Coronavirus nsp14 contains an ExoN domain with the conserved residues DE-D-D (50). Alanine substitutions for the DE-D-D residues of SARS-CoV nsp14 nearly abolished its nucleolytic activities in vitro (40). Interestingly, alanine substitutions for these active-site residues were lethal to HCoV-229E replication (40), while viable mutant MHV viruses with growth and RNA synthesis defects were recovered in the context of an infectious cDNA system (15). The reason for this marked difference is unclear, but it would be interesting to investigate whether ExoN mutants of IBV are viable and whether the nsp13-14 or nsp14-15 intermediate affects ExoN activity. nsp15 exhibits endoribonuclease activity (6, 29, 34, 65), and its critical role in coronavirus replication was established with a SARS-CoV replicon (1) and reverse-genetics-based mutagenesis studies of HCoV-229E and MHV (29, 34). The biochemical and structural data reveal an active site with a highly conserved H-H-K triad as putative catalytic residues and the requirement of nsp15 hexamerization for optimal enzymatic activity (25, 29, 33, 34, 45, 65). Mutation of either or both of the two catalytic H residues in MHV nsp15 resulted in greatly reduced but detectable endoribonuclease activity in vitro and the recovery of mutant viruses with only modest defects in viral RNA synthesis, suggesting that the nsp15 endoribonuclease is needed for optimal MHV replication (34). In contrast, alanine substitutions for a noncatalytic residue, D324, in MHV nsp15 and the corresponding D6408 in HCoV-229E nsp15 were lethal to virus production (29, 34). Expression of MHV nsp15 containing the D324A mutation and corresponding SARS-CoV nsp15 carrying D296A in bacteria resulted in formation of insoluble proteins (25, 34). Structurally, this Asp residue is highly conserved in the nsp15 orthologs of all coronaviruses and is involved in an extensive hydrogen bond network (45). The D296A mutation was believed to have structural effects on nsp15 that indirectly impair the function of this protein, suggesting that the structural integrity of nsp15 is essential for viral replication (45). In the present study, the nsp14-15 and nsp15-16 intermediates resulting from nsp14-15/Q5989N and nsp15-16/Q6327P may block the hexamer formation of IBV nsp15, completely abolishing the enzymatic activity.
Some stable intermediate cleavage products, such as the IBV nsp12-13 (160 kDa) and nsp14-15-16 (132 kDa) intermediates and MHV nsp4-10/11 (150 kDa) intermediates, were identified in virus-infected cells and were thought to play a role in virus replication (48, 64). A recent model suggests that unprocessed intermediates, such as MHV nsp4-10/11 (P150), function as a polyprotein involved in negative-stranded RNA synthesis, while mature proteins are required for positive-stranded RNA synthesis (46, 47). The results from this report revealed an intimate correlation between 3CLpro-mediated cleavage and regulation of RNA synthesis (Table (Table2).2). Mutation of P1-Gln to Asn at the nsp6-7, nsp12-13, nsp13-14, and nsp14-15 sites, as well as mutation of P1-Gln to Pro at the nsp15-16 site, which resulted in the complete abolishment of proteolytic processing, allowed negative-stranded genomic RNA replication but blocked sgRNA transcription. In contrast, mutation of P1-Gln to Asn at the nsp15-16 site allowed the synthesis of both genomic and sgRNA, even though the amount of cleaved products was extremely low due to significantly reduced processing efficiency. These findings suggested that unprocessed intermediates, such as the nsp6-7, nsp12-13, nsp13-14, nsp14-15, and nsp15-16 precursors, may function in genomic RNA replication, while mature proteins may be involved in sgRNA transcription. More extensive studies of cells infected with the recovered mutant viruses and transfected with full-length transcripts carrying the mutations would be required to support this assumption.
Coronavirus nsp15 is an endoribonuclease involved in RNA processing (6, 8, 29, 65), and nsp16 possesses 2′O-methyltransferase activity, thought to be responsible for RNA capping (11, 50, 60). Both proteins were proven to be essential for RNA synthesis (1, 34). In this study, low-level, noncytopathic infection was observed with cells infected with Q6327N mutant virus, due to inefficient transcription of positive-stranded sgRNA. It suggests that one of the main functions of mature nsp15 and nsp16 is at this step of the viral replication cycle. However, attempts to rescue the phenotype by coexpression of individual nsp15 and nsp16 with the full-length mutant transcripts in cells were unsuccessful. Further biochemical and functional characterization of the two proteins would be required to address this issue.
Perhaps the most intriguing result is the selection of 3CLpro carrying the P166S or P166L mutation. Both the P166S and P166L mutants specifically enhance the cleavage efficiency at the nsp15-16 site containing the P1-Asn substitution. Structurally, Pro-166 in IBV 3CLpro, corresponding to Pro-168 in SARS-CoV 3CLpro, is not conserved in the 3CLpros of MHV, TGEV, and HCoV-229E where a Ser or a Gly residue is located (66). Residues 164 to 168 form an antiparallel sheet with residues P6 to P1 of a substrate, but Pro-166 is not a required determinant for substrate specificity (66). It is possible that the P166S or P166L mutation may alter the interaction between amino acid residues in the substrate-binding pocket, leading to the increased proteolytic cleavage activity at the mutant nsp15-16 site and promoting virus replication. Alternatively, interaction between the mutant 3CLpro and the mutant nsp15-16 intermediate may induce conformational change of the substrate to facilitate 3CLpro-mediated cleavage. As this enhancement effect of the mutant 3CLpros is observed only with the mutant nsp15-16 site containing the P1-Asn substitution, the substrate that induces the emergence of the mutant enzyme, further structural studies would reveal novel insights into the enzyme substrate specificity and catalytic mechanisms of this important proteinase.
Published ahead of print on 5 May 2010.