In silico analysis of ORF1ab in coronavirus HKU1 genome reveals a unique putative cleavage site of coronavirus HKU1 3C-like protease

Proteolytic Processing of the Coronavirus Replicase Nonstructural Protein 14 Exonuclease Is Not Required for Virus Replication but Alters RNA Synthesis and Viral Fitness

Coronaviruses (CoVs) initiate replication by translation of the positive-sense RNA genome into the replicase polyproteins connecting 16 nonstructural protein domains (nsp1-16), which are subsequently processed by viral proteases to yield mature nsp. For the betacoronavirus murine hepatitis virus (MHV), total inhibition of translation or proteolytic processing of replicase polyproteins results in rapid cessation of RNA synthesis. The nsp5-3CLpro (Mpro) processes nsps7-16, which assemble into functional replication-transcription complexes (RTCs), including the enzymatic nsp12-RdRp and nsp14-exoribonuclease (ExoN)/N7-methyltransferase. The nsp14-ExoN activity mediates RNA-dependent RNA proofreading, high-fidelity RNA synthesis, and replication. To date, the solved partial RTC structures, biochemistry, and models use or assume completely processed, mature nsp. Here, we demonstrate that in MHV, engineered deletion of the cleavage sites between nsp13-14 and nsp14-15 allowed recovery of replication-competent virus. Compared to wild-type (WT) MHV, the nsp13-14 and nsp14-15 cleavage deletion mutants demonstrated delayed replication kinetics, impaired genome production, altered abundance and patterns of recombination, and impaired competitive fitness. Further, the nsp13-14 and nsp14-15 mutant viruses demonstrated mutation frequencies that were significantly higher than with the WT. The results demonstrate that cleavage of nsp13-14 or nsp14-15 is not required for MHV viability and that functions of the RTC/nsp14-ExoN are impaired when assembled with noncleaved intermediates. These data will inform future genetic, structural, biochemical, and modeling studies of coronavirus RTCs and nsp 13, 14, and 15 and may reveal new approaches for inhibition or attenuation of CoV infection.

IMPORTANCE Coronavirus replication requires proteolytic maturation of the nonstructural replicase proteins to form the replication-transcription complex. Coronavirus replication-transcription complex models assume mature subunits; however, mechanisms of coronavirus maturation and replicase complex formation have yet to be defined. Here, we show that for the coronavirus murine hepatitis virus, cleavage between the nonstructural replicase proteins nsp13-14 and nsp14-15 is not required for replication but does alter RNA synthesis and recombination. These results shed new light on the requirements for coronavirus maturation and replication-transcription complex assembly, and they may reveal novel therapeutic targets and strategies for attenuation.

Correlated sequence signatures are present within the genomic 5'UTR RNA and NSP1 protein in coronaviruses

The 5'UTR part of coronavirus genomes plays key roles in the viral replication cycle and translation of viral mRNAs. The first 75-80 nt, also called the leader sequence, are identical for genomic mRNA and subgenomic mRNAs. Recently, it was shown that cooperative actions of a 5'UTR segment and the nonstructural protein NSP1 are essential for both the inhibition of host mRNAs and for specific translation of viral mRNAs. Here, sequence analyses of both the 5'UTR RNA segment and the NSP1 protein have been done for several coronaviruses, with special attention to the betacoronaviruses. The conclusions are: (i) precise specific molecular signatures can be found in both the RNA and the NSP1 protein; (ii) both types of signatures correlate between each other. Indeed, definite sequence motifs in the RNA correlate with sequence motifs in the protein, indicating a coevolution between the 5'UTR and NSP1 in betacoronaviruses. Experimental mutational data on 5'UTR and NSP1 from SARS-CoV-2 using cell-free translation extracts support these conclusions and show that some conserved key residues in the amino-terminal half of the NSP1 protein are essential for evasion to the inhibitory effect of NSP1 on translation.

SARS-CoV-2 3CLpro whole human proteome cleavage prediction and enrichment/depletion analysis

A novel coronavirus (SARS-CoV-2) has devastated the globe as a pandemic that has killed millions of people. Widespread vaccination is still uncertain, so many scientific efforts have been directed toward discovering antiviral treatments. Many drugs are being investigated to inhibit the coronavirus main protease, 3CLpro, from cleaving its viral polyprotein, but few publications have addressed this protease’s interactions with the host proteome or their probable contribution to virulence. Too few host protein cleavages have been experimentally verified to fully understand 3CLpro’s global effects on relevant cellular pathways and tissues. Here, I set out to determine this protease’s targets and corresponding potential drug targets. Using a neural network trained on cleavages from 392 coronavirus proteomes with a Matthews correlation coefficient of 0.985, I predict that a large proportion of the human proteome is vulnerable to 3CLpro, with 4898 out of approximately 20,000 human proteins containing at least one putative cleavage site. These cleavages are nonrandomly distributed and are enriched in the epithelium along the respiratory tract, brain, testis, plasma, and immune tissues and depleted in olfactory and gustatory receptors despite the prevalence of anosmia and ageusia in COVID-19 patients. Affected cellular pathways include cytoskeleton/motor/cell adhesion proteins, nuclear condensation and other epigenetics, host transcription and RNAi, ribosomal stoichiometry and nascent-chain detection and degradation, ubiquitination, pattern recognition receptors, coagulation, lipoproteins, redox, and apoptosis. This whole proteome cleavage prediction demonstrates the importance of 3CLpro in expected and nontrivial pathways affecting virulence, lead me to propose more than a dozen potential therapeutic targets against coronaviruses, and should therefore be applied to all viral proteases and subsequently experimentally verified.

Screening of Severe Acute Respiratory Syndrome Coronavirus 2 RNA-Dependent RNA Polymerase Inhibitors Using Computational Approach

The detrimental effect of coronavirus disease 2019 (COVID-19) pandemic has manifested itself as a global crisis. Currently, no specific treatment options are available for COVID-19, so therapeutic interventions to tackle the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection must be urgently established. Therefore, cohesive and multidimensional efforts are required to identify new therapies or investigate the efficacy of small molecules and existing drugs against SARS-CoV-2. Since the RNA-dependent RNA Polymerase (RdRP) of SARS-CoV-2 is a promising therapeutic target, this study addresses the identification of antiviral molecules that can specifically target SARS-CoV-2 RdRP. The computational approach of drug development was used to screen the antiviral molecules from two antiviral libraries (Life Chemicals [LC] and ASINEX) against RdRP. Here, we report six antiviral molecules (F3407-4105, F6523-2250, F6559-0746 from LC and BDG 33693278, BDG 33693315, LAS 34156196 from ASINEX), which show substantial interactions with key amino acid residues of the active site of SARS-CoV-2 RdRP and exhibit higher binding affinity (>7.5 kcalmol^-1) than Galidesivir, an Food and Drug Administration-approved inhibitor of the same. Further, molecular dynamics simulation and Molecular Mechanics Poisson-Boltzmann Surface Area results confirmed that identified molecules with RdRP formed higher stable RdRP-inhibitor(s) complex than RdRP-Galidesvir complex. Our findings suggest that these molecules could be potential inhibitors of SARS-CoV-2 RdRP. However, further in vitro and preclinical experiments would be required to validate these potential inhibitors of SARS-CoV-2 protein.

High-throughput analysis of the interactions between viral proteins and host cell RNAs

RNA-protein interactions of a virus play a major role in the replication of RNA viruses. The replication and transcription of these viruses take place in the cytoplasm of the host cell; hence, there is a probability for the host RNA-viral protein and viral RNA-host protein interactions. The current study applies a high-throughput computational approach, including feature extraction and machine learning methods, to predict the affinity of protein sequences of ten viruses to three categories of RNA sequences. These categories include RNAs involved in the protein-RNA complexes stored in the RCSB database, the human miRNAs deposited at the mirBase database, and the lncRNA deposited in the LNCipedia database. The results show that evolution not only tries to conserve key viral proteins involved in the replication and transcription but also prunes their interaction capability. These proteins with specific interactions do not perturb the host cell through undesired interactions. On the other hand, the hypermutation rate of NSP3 is related to its affinity to host cell RNAs. The Gene Ontology (GO) analysis of the miRNA with affiliation to NSP3 suggests that these miRNAs show strongly significantly enriched GO terms related to the known symptoms of COVID-19. Docking and MD simulation study of the obtained miRNA through high-throughput analysis suggest a non-coding RNA (an RNA antitoxin, ToxI) as a natural aptamer drug candidate for NSP5 inhibition. Finally, a significant interplay of the host RNA-viral protein in the host cell can disrupt the host cell's system by influencing the RNA-dependent processes of the host cells, such as a differential expression in RNA. Furthermore, our results are useful to identify the side effects of mRNA-based vaccines, many of which are caused by the off-label interactions with the human lncRNAs.

FDA-approved thiol-reacting drugs that potentially bind into the SARS-CoV-2 main protease, essential for viral replication

Emergent novel SARS-CoV-2 is responsible for the current pandemic outbreak of severe acute respiratory syndrome with high mortality among the symptomatic population worldwide. Given the absence of a current vaccine or specific antiviral treatment, it is urgent to search for FDA-approved drugs that can potentially inhibit essential viral enzymes. The inhibition of 3CLpro has potential medical application, due to the fact that it is required for processing of the first translated replicase polyproteins into a series of native proteins, which are essential for viral replication in the host cell. We employed an in silico approach to test if disulfiram, as well as its metabolites, and captopril could be used as potential antiviral drugs against COVID-19. We provide data on the potential covalent interaction of disulfiram and its metabolites with the substrate binding subsite of 3CLpro and propose a possible mechanism for the irreversible protease inactivation thought the reaction of the aforementioned compounds with the Cys145. Although, captopril is shown to be a potential ligand of 3CLpro, it is not recommended anti-COVID-19 therapy, due to the fact that it can induce the expression of the viral cellular receptor such as, angiotensin-converting enzyme ACE-2, and thus, making the patient potentially more susceptible to infection. On the other hand, disulfiram, an alcoholism-averting drug, has been previously proposed as an antimicrobial and anti-SARS and MERS agent, safe to use even at higher doses with low side effects, it is recommended to be tested for control of SARS-CoV-2 infection.Communicated by Ramaswamy H. Sarma.

Major epigenetic factors associated with the novel coronavirus disease-2019 (COVID-19) severity

The worldwide spread and high rate of viral transmission and related morbidity and mortality of Coronavirus disease-19 (COVID-19) is a crisis. Some epigenetic determinants predispose individuals to severe infection. Patients with prior chronic medical illnesses (hypertension, diabetes, lupus, and chronic obstructive lung disease) are highly susceptible to the infection. The aging and diabetes pandemic possibly exacerbate the COVID-19 or SARS-CoV-2 pandemic by enhancing COVID-19 associated comorbidities. COVID-19 utilizes several proteins for tackling the host immune response associated with enhancing comorbidities. The angiotensin-converting enzyme (ACE) is a significant receptor for SARS-CoV-2, which significantly expresses higher among individuals with comorbidities and under stress conditions. Patients with systemic lupus erythematosus are also prone to be susceptible to the disease. Viral infections cause a defect in the DNA methylation in lupus, causing further ACE2 hypomethylation and overexpression, leading to viral binding and cytokine storm and tissue damage during COVID-19 infection. The microRNAs (miRNAs) epigenetics regulations also play a critical role in the suppression of immune responses. Meanwhile, viral proteins interplays with the host cell are conferred primarily through TGF-β and HIF-1 signaling, endocytosis, autophagy, and Toll-like receptor signaling RIG-I signaling, Il-17 signaling, and fatty acid oxidation/degradation. Furthermore, the COVID19 patient's metabolic states determine the infection severity. Noticeably, ten human metabolic proteins, including SGTA, SPECC1, FGL2, PHB, STAT3, BCL2L1, CAV1, JUN, PPP1CA, and XPO1, interact with the SARSE-CoV-2. Interactions between SARSCoV's spike protein-containing lipid-rich membrane compartments and epigenetic modulations are considered targets to inhibit the viral infection. Therefore, it seems that epigenetics plays a substantial role in the COVID-19 severity. Future in-depth studies will be promising. Vaccine design, particularly regarding ACE viral receptor monoclonal antibodies, is a proposal alongside adhering to personal hygiene.

Application of Tensor Decomposition to Gene Expression of Infection of Mouse Hepatitis Virus Can Identify Critical Human Genes and Efffective Drugs for SARS-CoV-2 Infection

To better understand the genes with altered expression caused by infection with the novel coronavirus strain SARS-CoV-2 causing COVID-19 infectious disease, a tensor decomposition (TD)-based unsupervised feature extraction (FE) approach was applied to a gene expression profile dataset of the mouse liver and spleen with experimental infection of mouse hepatitis virus, which is regarded as a suitable model of human coronavirus infection. TD-based unsupervised FE selected 134 altered genes, which were enriched in protein-protein interactions with orf1ab, polyprotein, and 3C-like protease that are well known to play critical roles in coronavirus infection, suggesting that these 134 genes can represent the coronavirus infectious process. We then selected compounds targeting the expression of the 134 selected genes based on a public domain database. The identified drug compounds were mainly related to known antiviral drugs, several of which were also included in those previously screened with an in silico method to identify candidate drugs for treating COVID-19.

Understanding genomic diversity, pan-genome, and evolution of SARS-CoV-2

Coronovirus disease 2019 (COVID-19) infection, which originated from Wuhan, China, has seized the whole world in its grasp and created a huge pandemic situation before humanity. Since December 2019, genomes of numerous isolates have been sequenced and analyzed for testing confirmation, epidemiology, and evolutionary studies. In the first half of this article, we provide a detailed review of the history and origin of COVID-19, followed by the taxonomy, nomenclature and genome organization of its causative agent Severe Acute Respiratory Syndrome-related Coronavirus-2 (SARS-CoV-2). In the latter half, we analyze subgenus Sarbecovirus (167 SARS-CoV-2, 312 SARS-CoV, and 5 Pangolin CoV) genomes to understand their diversity, origin, and evolution, along with pan-genome analysis of genus Betacoronavirus members. Whole-genome sequence-based phylogeny of subgenus Sarbecovirus genomes reasserted the fact that SARS-CoV-2 strains evolved from their common ancestors putatively residing in bat or pangolin hosts. We predicted a few country-specific patterns of relatedness and identified mutational hotspots with high, medium and low probability based on genome alignment of 167 SARS-CoV-2 strains. A total of 100-nucleotide segment-based homology studies revealed that the majority of the SARS-CoV-2 genome segments are close to Bat CoV, followed by some to Pangolin CoV, and some are unique ones. Open pan-genome of genus Betacoronavirus members indicates the diversity contributed by the novel viruses emerging in this group. Overall, the exploration of the diversity of these isolates, mutational hotspots and pan-genome will shed light on the evolution and pathogenicity of SARS-CoV-2 and help in developing putative methods of diagnosis and treatment.

Analysis of SARS-CoV-2 RNA-dependent RNA polymerase as a potential therapeutic drug target using a computational approach

Background

The Severe acute respiratory syndrome-related coronavirus 2 (SARS-CoV-2) outbreak originating in Wuhan, China, has raised global health concerns and the pandemic has now been reported on all inhabited continents. Hitherto, no antiviral drug is available to combat this viral outbreak.

Methods

Keeping in mind the urgency of the situation, the current study was designed to devise new strategies for drug discovery and/or repositioning against SARS-CoV-2. In the current study, RNA-dependent RNA polymerase (RdRp), which regulates viral replication, is proposed as a potential therapeutic target to inhibit viral infection.

Results

Evolutionary studies of whole-genome sequences of SARS-CoV-2 represent high similarity (> 90%) with other SARS viruses. Targeting the RdRp active sites, ASP760 and ASP761, by antiviral drugs could be a potential therapeutic option for inhibition of coronavirus RdRp, and thus viral replication. Target-based virtual screening and molecular docking results show that the antiviral Galidesivir and its structurally similar compounds have shown promise against SARS-CoV-2.

Conclusions

The anti-polymerase drugs predicted here—CID123624208 and CID11687749—may be considered for in vitro and in vivo clinical trials.

Analysis of SARS-CoV-2 RNA-dependent RNA polymerase as a potential therapeutic drug target using a computational approach

BACKGROUND

The Severe acute respiratory syndrome-related coronavirus 2 (SARS-CoV-2) outbreak originating in Wuhan, China, has raised global health concerns and the pandemic has now been reported on all inhabited continents. Hitherto, no antiviral drug is available to combat this viral outbreak.

METHODS

Keeping in mind the urgency of the situation, the current study was designed to devise new strategies for drug discovery and/or repositioning against SARS-CoV-2. In the current study, RNA-dependent RNA polymerase (RdRp), which regulates viral replication, is proposed as a potential therapeutic target to inhibit viral infection.

RESULTS

Evolutionary studies of whole-genome sequences of SARS-CoV-2 represent high similarity (> 90%) with other SARS viruses. Targeting the RdRp active sites, ASP760 and ASP761, by antiviral drugs could be a potential therapeutic option for inhibition of coronavirus RdRp, and thus viral replication. Target-based virtual screening and molecular docking results show that the antiviral Galidesivir and its structurally similar compounds have shown promise against SARS-CoV-2.

CONCLUSIONS

The anti-polymerase drugs predicted here-CID123624208 and CID11687749-may be considered for in vitro and in vivo clinical trials.

Proteolytic processing of mesonivirus replicase polyproteins by the viral 3C-like protease

Mesoniviridae are a family of insect RNA viruses that diverged profoundly from other families of the Nidovirales. Mesonivirus replicative proteins are produced from large polyprotein (pp) precursors (pp1a and pp1ab) through proteolytic cleavage by the viral 3C-like protease (3CLpro) and, possibly, other proteases. Using recombinant forms of the Cavally virus 3CLpro and pp1a/pp1ab-derived substrates, we characterized 3CLpro cleavage sites in mesonivirus polyproteins. Our data lead us to suggest that 3CLpro cleaves the central and C-proximal regions of mesonivirus pp1a/pp1ab at 12 conserved sites. Compared to other nidovirus homologues, the mesonivirus 3CLpro features a distinct substrate specificity, with asparagine at P2 being a major specificity determinant. Furthermore, we provide evidence that expression of the ORF1b-encoded part of pp1ab involves a -1 ribosomal frameshift at a conserved GGAUUUU heptanucleotide sequence in the ORF1a/1b overlap region. Taken together, the study identifies critical steps in the expression and maturation of mesonivirus replicative proteins.

Middle East respiratory syndrome coronavirus: another zoonotic betacoronavirus causing SARS-like disease

The source of the severe acute respiratory syndrome (SARS) epidemic was traced to wildlife market civets and ultimately to bats. Subsequent hunting for novel coronaviruses (CoVs) led to the discovery of two additional human and over 40 animal CoVs, including the prototype lineage C betacoronaviruses, Tylonycteris bat CoV HKU4 and Pipistrellus bat CoV HKU5; these are phylogenetically closely related to the Middle East respiratory syndrome (MERS) CoV, which has affected more than 1,000 patients with over 35% fatality since its emergence in 2012. All primary cases of MERS are epidemiologically linked to the Middle East. Some of these patients had contacted camels which shed virus and/or had positive serology. Most secondary cases are related to health care-associated clusters. The disease is especially severe in elderly men with comorbidities. Clinical severity may be related to MERS-CoV's ability to infect a broad range of cells with DPP4 expression, evade the host innate immune response, and induce cytokine dysregulation. Reverse transcription-PCR on respiratory and/or extrapulmonary specimens rapidly establishes diagnosis. Supportive treatment with extracorporeal membrane oxygenation and dialysis is often required in patients with organ failure. Antivirals with potent in vitro activities include neutralizing monoclonal antibodies, antiviral peptides, interferons, mycophenolic acid, and lopinavir. They should be evaluated in suitable animal models before clinical trials. Developing an effective camel MERS-CoV vaccine and implementing appropriate infection control measures may control the continuing epidemic.

Human coronaviruses: Clinical features and phylogenetic analysis

Strains of human coronavirus (HCoV), namely HCoV-OC43, HCoV-229E, HCoV-NL63, and HCoV-HKU1, primarily infect the upper respiratory and gastrointestinal tracts and are the most common cause of non-rhinovirus-induced common cold in humans. Although the manifestations of coronavirus infection (i.e., rhinorrhea, sneezing, cough, nasal obstruction, and bronchitis) are generally self-limiting in healthy adults, certain strains such as HCoV-NL63 and HCoV-HKU1 can cause severe lower respiratory tract infection and febrile seizure, especially in infants, people of advanced age, and immunocompromised hosts. In 2003, a novel HCoV strain was identified as the causative agent of the severe acute respiratory syndrome (SARS) epidemic that began in Asia in 2002. The strain has hence been referred to as SARS-CoV. In addition, as recently as September 2012, another novel HCoV, human betacoronavirus 2c EMC2012, was identified as being the cause of fever, renal failure, pneumonia, and severe respiratory distress in two patients in the Middle East. Phylogenetic analysis has revealed highly conserved sequences of ORF1ab, spike, nucleocapsid, and envelope protein genes, but not membrane protein genes, between human betacoronavirus 2c EMC2012 and SARS-CoV. This review focuses on the differences in the genomes of certain HCoV strains, the pathogenesis of said strains, and recent developments in the establishment of therapeutic agents that might aid in the treatment of patients with such infections.

Characterization of Bafinivirus main protease autoprocessing activities

The production of functional nidovirus replication-transcription complexes involves extensive proteolytic processing by virus-encoded proteases. In this study, we characterized the viral main protease (M(pro)) of the type species, White bream virus (WBV), of the newly established genus Bafinivirus (order Nidovirales, family Coronaviridae, subfamily Torovirinae). Comparative sequence analysis and mutagenesis data confirmed that the WBV M(pro) is a picornavirus 3C-like serine protease that uses a Ser-His-Asp catalytic triad embedded in a predicted two-β-barrel fold, which is extended by a third domain at its C terminus. Bacterially expressed WBV M(pro) autocatalytically released itself from flanking sequences and was able to mediate proteolytic processing in trans. Using N-terminal sequencing of autoproteolytic processing products we tentatively identified Gln↓(Ala, Thr) as a substrate consensus sequence. Mutagenesis data provided evidence to suggest that two conserved His and Thr residues are part of the S1 subsite of the enzyme's substrate-binding pocket. Interestingly, we observed two N-proximal and two C-proximal autoprocessing sites in the bacterial expression system. The detection of two major forms of M(pro), resulting from processing at two different N-proximal and one C-proximal site, in WBV-infected epithelioma papulosum cyprini cells confirmed the biological relevance of the biochemical data obtained in heterologous expression systems. To our knowledge, the use of alternative M(pro) autoprocessing sites has not been described previously for other nidovirus M(pro) domains. The data presented in this study lend further support to our previous conclusion that bafiniviruses represent a distinct group of viruses that significantly diverged from other phylogenetic clusters of the order Nidovirales.

Coronavirus genomics and bioinformatics analysis

The drastic increase in the number of coronaviruses discovered and coronavirus genomes being sequenced have given us an unprecedented opportunity to perform genomics and bioinformatics analysis on this family of viruses. Coronaviruses possess the largest genomes (26.4 to 31.7 kb) among all known RNA viruses, with G + C contents varying from 32% to 43%. Variable numbers of small ORFs are present between the various conserved genes (ORF1ab, spike, envelope, membrane and nucleocapsid) and downstream to nucleocapsid gene in different coronavirus lineages. Phylogenetically, three genera, Alphacoronavirus, Betacoronavirus and Gammacoronavirus, with Betacoronavirus consisting of subgroups A, B, C and D, exist. A fourth genus, Deltacoronavirus, which includes bulbul coronavirus HKU11, thrush coronavirus HKU12 and munia coronavirus HKU13, is emerging. Molecular clock analysis using various gene loci revealed that the time of most recent common ancestor of human/civet SARS related coronavirus to be 1999-2002, with estimated substitution rate of 4×10(-4) to 2×10(-2) substitutions per site per year. Recombination in coronaviruses was most notable between different strains of murine hepatitis virus (MHV), between different strains of infectious bronchitis virus, between MHV and bovine coronavirus, between feline coronavirus (FCoV) type I and canine coronavirus generating FCoV type II, and between the three genotypes of human coronavirus HKU1 (HCoV-HKU1). Codon usage bias in coronaviruses were observed, with HCoV-HKU1 showing the most extreme bias, and cytosine deamination and selection of CpG suppressed clones are the two major independent biological forces that shape such codon usage bias in coronaviruses.

Functional and genetic studies of the substrate specificity of coronavirus infectious bronchitis virus 3C-like proteinase

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.

More and More Coronaviruses: Human Coronavirus HKU1

After human coronaviruses OC43, 229E and NL63, human coronavirus HKU1 (HCoV-HKU1) is the fourth human coronavirus discovered. HCoV-HKU1 is a group 2a coronavirus that is still not cultivable. The G + C contents of HCoV-HKU1 genomes are 32%, the lowest among all known coronaviruses with complete genome sequences available. Among all coronaviruses, HCoV-HKU1 shows the most extreme codon usage bias, attributed most importantly to severe cytosine deamination. All HCoV-HKU1 genomes contain unique tandem copies of a 30-base acidic tandem repeat of unknown function at the N-terminus of nsp3 inside the acidic domain upstream of papain-like protease 1. Three genotypes, A, B and C, of HCoV-HKU1 and homologous recombination among their genomes, are observed. The incidence of HCoV-HKU1 infections is the highest in winter. Similar to other human coronaviruses, HCoV-HKU1 infections have been reported globally, with a median (range) incidence of 0.9 (0 – 4.4) %. HCoV-HKU1 is associated with both upper and lower respiratory tract infections that are mostly self-limiting. The most common method for diagnosing HCoV-HKU1 infection is RT-PCR or real-time RT-PCR using RNA extracted from respiratory tract samples such as nasopharyngeal aspirates (NPA). Both the pol and nucleocapsid genes have been used as the targets for amplification. Monoclonal antibodies have been generated for direct antigen detection in NPA. For antibody detection, Escherichia coli BL21 and baculovirus-expressed recombinant nucleocapsid of HCoV-HKU1 have been used for IgG and IgM detection in sera of patients and normal individuals, using Western blot and enzyme-linked immunoassay.

Expression, crystallization and preliminary crystallographic study of human coronavirus HKU1 nonstructural protein 9

Human coronavirus HKU1 (HCoV-HKU1) belongs to coronavirus group II and encodes 16 nonstructural proteins (nsps) which mediate genome replication and transcription. Among these nsps, nsp9 has been shown to possess single-stranded DNA/RNA-binding properties. The gene that encodes HCoV-HKU1 nsp9 was cloned and expressed in Escherichia coli and the protein was subjected to crystallization trials. The crystals diffracted to 2.7 A resolution and belonged to space group P2(1)2(1)2, with unit-cell parameters a = 83.5, b = 88.4, c = 31.2 A, alpha = beta = gamma = 90 degrees and two molecules per asymmetric unit.

Comparative analysis of complete genome sequences of three avian coronaviruses reveals a novel group 3c coronavirus

In this territory-wide molecular epidemiology study of coronaviruses (CoVs) in Hong Kong involving 1,541 dead wild birds, three novel CoVs were identified in three different bird families (bulbul CoV HKU11 [BuCoV HKU11], thrush CoV HKU12 [ThCoV HKU12], and munia CoV HKU13 [MuCoV HKU13]). Four complete genomes of the three novel CoVs were sequenced. Their genomes (26,396 to 26,552 bases) represent the smallest known CoV genomes. In phylogenetic trees constructed using chymotrypsin-like protease (3CL(pro)), RNA-dependent RNA polymerase (Pol), helicase, spike, and nucleocapsid proteins, BuCoV HKU11, ThCoV HKU12, and MuCoV HKU13 formed a cluster distantly related to infectious bronchitis virus and turkey CoV (group 3a CoVs). For helicase, spike, and nucleocapsid, they were also clustered with a CoV recently discovered in Asian leopard cats, for which the complete genome sequence was not available. The 3CL(pro), Pol, helicase, and nucleocapsid of the three CoVs possessed higher amino acid identities to those of group 3a CoVs than to those of group 1 and group 2 CoVs. Unique genomic features distinguishing them from other group 3 CoVs include a distinct transcription regulatory sequence and coding potential for small open reading frames. Based on these results, we propose a novel CoV subgroup, group 3c, to describe this distinct subgroup of CoVs under the group 3 CoVs. Avian CoVs are genetically more diverse than previously thought and may be closely related to some newly identified mammalian CoVs. Further studies would be important to delineate whether the Asian leopard cat CoV was a result of interspecies jumping from birds, a situation analogous to that of bat and civet severe acute respiratory syndrome CoVs.

Structure of the main protease from a global infectious human coronavirus, HCoV-HKU1

The newly emergent human coronavirus HKU1 (HCoV-HKU1) was first identified in Hong Kong in 2005. Infection by HCoV-HKU1 occurs worldwide and causes syndromes such as the common cold, bronchitis, and pneumonia. The CoV main protease (M(pro)), which is a key enzyme in viral replication via the proteolytic processing of the replicase polyproteins, has been recognized as an attractive target for rational drug design. In this study, we report the structure of HCoV-HKU1 M(pro) in complex with a Michael acceptor, inhibitor N3. The structure of HCoV-HKU1 provides a high-quality model for group 2A CoVs, which are distinct from group 2B CoVs such as severe acute respiratory syndrome CoV. The structure, together with activity assays, supports the relative conservation at the P1 position that was discovered by sequencing the HCoV-HKU1 genome. Combined with structural data from other CoV M(pro)s, the HCoV-HKU1 M(pro) structure reported here provides insights into both substrate preference and the design of antivirals targeting CoVs.

CoVDB: a comprehensive database for comparative analysis of coronavirus genes and genomes

The recent SARS epidemic has boosted interest in the discovery of novel human and animal coronaviruses. By July 2007, more than 3000 coronavirus sequence records, including 264 complete genomes, are available in GenBank. The number of coronavirus species with complete genomes available has increased from 9 in 2003 to 25 in 2007, of which six, including coronavirus HKU1, bat SARS coronavirus, group 1 bat coronavirus HKU2, groups 2c and 2d coronaviruses, were sequenced by our laboratory. To overcome the problems we encountered in the existing databases during comparative sequence analysis, we built a comprehensive database, CoVDB (http://covdb.microbiology.hku.hk), of annotated coronavirus genes and genomes. CoVDB provides a convenient platform for rapid and accurate batch sequence retrieval, the cornerstone and bottleneck for comparative gene or genome analysis. Sequences can be directly downloaded from the website in FASTA format. CoVDB also provides detailed annotation of all coronavirus sequences using a standardized nomenclature system, and overcomes the problems of duplicated and identical sequences in other databases. For complete genomes, a single representative sequence for each species is available for comparative analysis such as phylogenetic studies. With the annotated sequences in CoVDB, more specific blast search results can be generated for efficient downstream analysis.

Cytosine deamination and selection of CpG suppressed clones are the two major independent biological forces that shape codon usage bias in coronaviruses

Using the complete genome sequences of 19 coronavirus genomes, we analyzed the codon usage bias, dinucleotide relative abundance and cytosine deamination in coronavirus genomes. Of the eight codons that contain CpG, six were markedly suppressed. The mean NNU/NNC ratio of the six amino acids using either NNC or NNU as codon is 3.262, suggesting cytosine deamination. Among the 16 dinucleotides, CpG was most markedly suppressed (mean relative abundance 0.509). No correlation was observed between CpG abundance and mean NNU/NNC ratio. Among the 19 coronaviruses, CoV-HKU1 showed the most extreme codon usage bias and extremely high NNU/NNC ratio of 8.835. Cytosine deamination and selection of CpG suppressed clones by the immune system are the two major independent biochemical and biological selective forces that shape codon usage bias in coronavirus genomes. The underlying mechanism for the extreme codon usage bias, cytosine deamination and G + C content in CoV-HKU1 warrants further studies.

Complete genome sequence of bat coronavirus HKU2 from Chinese horseshoe bats revealed a much smaller spike gene with a different evolutionary lineage from the rest of the genome

Apart from bat-SARS-CoV, we have identified a novel group 1 coronavirus, bat-CoV HKU2, in Rhinolophus sinicus (Chinese horseshoe bats). Since it has been suggested that the receptor-binding motif (RBM) of SARS-CoV may have been acquired from a group 1 coronavirus, we conducted a surveillance study and identified bat-SARS-CoV and bat-CoV HKU2 in 8.7% and 7.5% respectively of R. sinicus in Hong Kong and Guangdong. Complete genome sequencing of four strains of bat-CoV HKU2 revealed the smallest coronavirus genome (27164 nucleotides) and a unique spike protein evolutionarily distinct from the rest of the genome. This spike protein, sharing similar deletions with other group 2 coronaviruses in its C-terminus, also contained a 15-amino acid peptide homologous to a corresponding peptide within the RBM of spike protein of SARS-CoV, which was absent in other coronaviruses except bat-SARS-CoV. These suggest a common evolutionary origin in the spike protein of bat-CoV HKU2, bat-SARS-CoV, and SARS-CoV.

Comparative analysis of twelve genomes of three novel group 2c and group 2d coronaviruses reveals unique group and subgroup features

Twelve complete genomes of three novel coronaviruses-bat coronavirus HKU4 (bat-CoV HKU4), bat-CoV HKU5 (putative group 2c), and bat-CoV HKU9 (putative group 2d)-were sequenced. Comparative genome analysis showed that the various open reading frames (ORFs) of the genomes of the three coronaviruses had significantly higher amino acid identities to those of other group 2 coronaviruses than group 1 and 3 coronaviruses. Phylogenetic trees constructed using chymotrypsin-like protease, RNA-dependent RNA polymerase, helicase, spike, and nucleocapsid all showed that the group 2a and 2b and putative group 2c and 2d coronaviruses are more closely related to each other than to group 1 and 3 coronaviruses. Unique genomic features distinguishing between these four subgroups, including the number of papain-like proteases, the presence or absence of hemagglutinin esterase, small ORFs between the membrane and nucleocapsid genes and ORFs (NS7a and NS7b), bulged stem-loop and pseudoknot structures downstream of the nucleocapsid gene, transcription regulatory sequence, and ribosomal recognition signal for the envelope gene, were also observed. This is the first time that NS7a and NS7b downstream of the nucleocapsid gene has been found in a group 2 coronavirus. The high Ka/Ks ratio of NS7a and NS7b in bat-CoV HKU9 implies that these two group 2d-specific genes are under high selective pressure and hence are rapidly evolving. The four subgroups of group 2 coronaviruses probably originated from a common ancestor. Further molecular epidemiological studies on coronaviruses in the bats of other countries, as well as in other animals, and complete genome sequencing will shed more light on coronavirus diversity and their evolutionary histories.

Human coronavirus NL63, a new respiratory virus

From the mid-1960s onwards, it was believed that only two human coronavirus species infect humans: HCoV-229E and HCoV-OC43. Then, in 2003, a novel member of the coronavirus family was introduced into the human population: SARS-CoV, causing an aggressive lung disease. Fortunately, this virus was soon expelled from the human population, but it quickly became clear that the human coronavirus group contains more members then previously assumed, with HCoV-NL63 identified in 2004. Despite its recent discovery, ample results from HCoV-NL63 research have been described. We present an overview of the publications on this novel coronavirus.

Comparative analysis of 22 coronavirus HKU1 genomes reveals a novel genotype and evidence of natural recombination in coronavirus HKU1

We sequenced and compared the complete genomes of 22 strains of coronavirus HKU1 (CoV HKU1) obtained from nasopharyngeal aspirates of patients with respiratory tract infections over a 2-year period. Phylogenetic analysis of 24 putative proteins and polypeptides showed that the 22 CoV HKU1 strains fell into three clusters (genotype A, 13 strains; genotype B, 3 strains and genotype C, 6 strains). However, different phylogenetic relationships among the three clusters were observed in different regions of their genomes. From nsp4 to nsp6, the genotype A strains were clustered with the genotype B strains. For nsp7 and nsp8 and from nsp10 to nsp16, the genotype A strains were clustered with the genotype C strains. From hemagglutinin esterase (HE) to nucleocapsid (N), the genotype B strains were clustered closely with the genotype C strains. Bootscan analysis showed possible recombination between genotypes B and C from nucleotide positions 11,500 to 13,000, corresponding to the nsp6-nsp7 junction, giving rise to genotype A, and between genotypes A and B from nucleotide positions 21,500 to 22,500, corresponding to the nsp16-HE junction, giving rise to genotype C. Multiple alignments further narrowed the sites of crossover to a 143-bp region between nucleotide positions 11,750 and 11,892 and a 29-bp region between nucleotide positions 21,502 and 21,530. Genome analysis also revealed various numbers of tandem copies of a perfect 30-base acidic tandem repeat (ATR) which encodes NDDEDVVTGD and various numbers and sequences of imperfect repeats in the N terminus of nsp3 inside the acidic domain upstream of papain-like protease 1 among the 22 genomes. All 10 CoV HKU1 strains with incomplete imperfect repeats (1.4 and 4.4) belonged to genotype A. The present study represents the first evidence for natural recombination in coronavirus associated with human infection. Analysis of a single gene is not sufficient for the genotyping of CoV HKU1 strains but requires amplification and sequencing of at least two gene loci, one from nsp10 to nsp16 (e.g., pol or helicase) and another from HE to N (e.g., spike or N). Further studies will delineate whether the ATR is useful for the molecular typing of CoV HKU1.

Molecular diversity of coronaviruses in bats

The existence of coronaviruses in bats is unknown until the recent discovery of bat-SARS-CoV in Chinese horseshoe bats and a novel group 1 coronavirus in other bat species. Among 309 bats of 13 species captured from 20 different locations in rural areas of Hong Kong over a 16-month period, coronaviruses were amplified from anal swabs of 37 (12%) bats by RT-PCR. Phylogenetic analysis of RNA-dependent-RNA-polymerase (pol) and helicase genes revealed six novel coronaviruses from six different bat species, in addition to the two previously described coronaviruses. Among the six novel coronaviruses, four were group 1 coronaviruses (bat-CoV HKU2 from Chinese horseshoe bat, bat-CoV HKU6 from rickett's big-footed bat, bat-CoV HKU7 from greater bent-winged bat and bat-CoV HKU8 from lesser bent-winged bat) and two were group 2 coronaviruses (bat-CoV HKU4 from lesser bamboo bats and bat-CoV HKU5 from Japanese pipistrelles). An astonishing diversity of coronaviruses was observed in bats.