Crippled Coronavirus: 5'-PolyU targeted Oligo prevents development of infectious Virions

Current RNA viral pandemic of COVID-19 has been worsened by rapidly spreading viral variants. To inhibit mutation-based development of new escape variants, elements that are indispensable for the virus may be targeted. The 5-polyU tract of the antigenome offers one such target. Host cells do not harbor 5-polyU tracts on any of their transcripts, making the tract an attractive virus-specific target. We hypothesize that inhibiting the 5-polyU by complementary oligonucleotide can limit the use of the tract as template for virus to generate 3 polyA tails of RNA viruses. Here, we used a frameshift-inducing DNA oligonucleotide with 3 polyAs to target the 5-polyU tract of mouse coronavirus (MHV-A59). Results from assays for double stranded RNA (dsRNA) synthesis, infectivity of released virions, and syncytium formation indicate that the oligonucleotide treatment prevented generation of infectious virions. Our results show a unique mode of action of the designed 3-polyA oligonucleotide against mouse coronavirus which leaves host cells unaffected. This strategy can be adopted for the development of novel classes of oligonucleotide-based drugs that inhibit the production of infectious RNA viruses, including the coronaviruses. Since the 5-polyU tract is conserved and is essential for variants of coronaviruses, this strategy can potentially address coronavirus variant emergence as well.

SARS-CoV-2 triggers Golgi fragmentation via downregulation of GRASP55 to facilitate viral trafficking

The ongoing COVID-19 pandemic is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), an enveloped RNA virus. Despite the high economic and life losses caused by SARS-CoV-2, the detailed viral cycle, especially how it assembles and traffics in the secretory pathway, remains largely unknown. Here, we show that SARS-CoV-2 infection induces global alterations of the host endomembrane system, including dramatic Golgi fragmentation. Disrupting Golgi function with small molecules strongly inhibits viral infection. Furthermore, expression of several SARS-CoV-2 proteins individually is sufficient to trigger Golgi fragmentation. Significantly, SARS-CoV-2 infection down-regulates GRASP55 but up-regulates TGN46 expression, while expression of GRASP55 or knockdown of TGN46 reduces the infection rate of both USA-WA1 and Delta variants of SARS-CoV-2. Our study reveals that SARS-CoV-2 modulates Golgi structure and function via altering GRASP55 and TGN46 expression to facilitate viral trafficking, indicating the Golgi as a novel therapeutic target to block SARS-CoV-2 infection.

Genome-wide characterization of SARS-CoV-2 cytopathogenic proteins in the search of antiviral targets

Therapeutic inhibition of critical viral functions is important for curtailing coronavirus disease-2019 (COVID-19). We sought to identify antiviral targets through genome-wide characterization of SARS-CoV-2 proteins that are crucial for viral pathogenesis and that cause harmful cytopathic effects. All twenty-nine viral proteins were tested in a fission yeast cell-based system using inducible gene expression. Twelve proteins including eight non-structural proteins (NSP1, NSP3, NSP4, NSP5, NSP6, NSP13, NSP14 and NSP15) and four accessory proteins (ORF3a, ORF6, ORF7a and ORF7b) were identified that altered cellular proliferation and integrity, and induced cell death. Cell death correlated with the activation of cellular oxidative stress. Of the twelve proteins, ORF3a was chosen for further study in mammalian cells. In human pulmonary and kidney epithelial cells, ORF3a induced cellular oxidative stress associated with apoptosis and necrosis, and caused activation of pro-inflammatory response with production of the cytokines TNF-, IL-6, and IFN-{beta}1, possibly through the activation of NF-{kappa}B. To further characterize the mechanism, we tested a natural ORF3a Beta variant, Q57H, and a mutant with deletion of the highly conserved residue, {Delta}G188. Compared to wild type ORF3a, the {Delta}G188 variant yielded more robust activation of cellular oxidative stress, cell death, and innate immune response. Since cellular oxidative stress and inflammation contribute to cell death and tissue damage linked to the severity of COVID-19, our findings suggest that ORF3a is a promising, novel therapeutic target against COVID-19. SignificanceThe ongoing SARS-CoV-2 pandemic has claimed over 5 million lives with more than 250 million people infected world-wide. While vaccines are effective, the emergence of new viral variants could jeopardize vaccine protection. Antiviral drugs provide an alternative to battle against COVID-19. Our goal was to identify viral therapeutic targets that can be used in antiviral drug discovery. Utilizing a genome-wide functional analysis in a fission yeast cell-based system, we identified twelve viral candidates, including ORF3a, which cause cellular oxidative stress, inflammation and apoptosis and necrosis that contribute to COVID-19. Our findings indicate that antiviral agents targeting ORF3a could greatly impact COVID-19.

The Achilles' heel of coronaviruses: targeting the 5' Polyuridines tract of the antigenome to inhibit Mouse coronavirus virus-induced cell death

The current coronavirus pandemic situation is worsened by the rapidly-spreading SARS-CoV-2 virus variants. Identification of viral targets that are indispensable for the virus can be targeted to inhibit mutation-based new escape variant development. The 5-polyU tract of the antigenome offers such a target. Host cells do not harbor 5-polyU tracts on any of their transcripts, making the tract an attractive, virus-specific target. Inhibiting the 5-polyU can limit the use of the tract as template to generate 3 polyA tails of +RNAs of coronaviruses. Here, a modified DNA oligo with 3 polyAs is used to target the 5-polyU tract in mouse coronavirus (MHV-A59). The oligo treatment in mouse 17CL-1 cells infected with MHV-A59 significantly prevented virus-induced cell deaths. This proof-of-concept result shows a unique mode of action against mouse coronavirus without affecting host cells, and can be used for the development of novel classes of drugs that inhibit coronavirus infection in host cells, specifically by the COVID-19-causing virus SARS-CoV-2. In addition, as the 5-polyU tract is immediately generated upon infection, the tag can also be targeted for reliable early detection of viral infection.

Computational saturation mutagenesis of SARS-CoV-1 spike glycoprotein: stability, binding affinity, and comparison with SARS-CoV-2

Severe Acute respiratory syndrome coronavirus (SARS-CoV-1) attaches to the host cell surface to initiate the interaction between the receptor-binding domain (RBD) of its spike glycoprotein (S) and the human Angiotensin-converting enzyme (hACE2) receptor. SARS-CoV-1 mutates frequently because of its RNA genome, which challenges the antiviral development. Here, we performed computational saturation mutagenesis of the S protein of SARS-CoV-1 to identify the residues crucial for its functions. We used the structure-based energy calculations to analyze the effects of the missense mutations on the SARS-CoV-1 S stability and the binding affinity with hACE2. The sequence and structure alignment showed similarities between the S proteins of SARS-CoV-1 and SARS-CoV-2. Interestingly, we found that target mutations of S protein amino acids generate similar effects on their stabilities between SARS-CoV-1 and SARS-CoV-2. For example, G839W of SARS-CoV-1 corresponds to G857W of SARS-CoV-2, which decrease the stability of their S glycoproteins. The viral mutation analysis of the two different SARS-CoV-1 isolates showed that mutations, T487S and L472P, weakened the S-hACE2 binding of the 2003-2004 SARS-CoV-1 isolate. In addition, the mutations of L472P and F360S destabilized the 2003-2004 viral isolate. We further predicted that many mutations on N-linked glycosylation sites would increase the stability of the S glycoprotein. Our results can be of therapeutic importance in the design of antivirals or vaccines against SARS-CoV-1 and SARS-CoV-2. Author SummarySevere acute respiratory syndrome coronavirus (SARS-CoV-1) is an RNA virus that undergoes frequent mutations, which may result in more virulent SARS-CoV-1 variants. To prevent another pandemic in the future, scientists must understand the mechanisms of viral mutations and predict if any variants could become a dominant. The infection of SARS-CoV-1 in cells is largely depending on the interactions of the viral Spike (S) and human angiotensin-converting enzyme 2 (hACE2). We applied a computational method to predict S missense mutations that will make SARS-CoV-1 more virulent. We are interested in the variants that can change SARS-CoV-1 spike protein stability and/or change the virus-receptor interactions. We mutated each residue of SARS-CoV-1 spike to all possible amino acids; we calculated the differences between the folding energy and binding energy of each variant and the wildtype and identified the target S mutations with significant effects on protein stability and protein-protein interaction. We found some viral mutations could destabilize S and weaken S-hACE2 binding of SARS-CoV-1 isolate. Our results show that the computational saturation mutagenesis is a reliable approach in the analysis and prediction of missense mutations.

Circular RNA profiling reveals abundant and diverse circRNAs of SARS-CoV-2, SARS-CoV and MERS-CoV origin

Circular RNAs (circRNAs) encoded by DNA genomes have been identified across host and pathogen species as parts of the transcriptome. Accumulating evidences indicate that circRNAs play critical roles in autoimmune diseases and viral pathogenesis. Here we report that RNA viruses of the Betacoronavirus genus of Coronaviridae, SARS-CoV-2, SARS-CoV and MERS-CoV, encode a novel type of circRNAs. Through de novo circRNA analyses of publicly available coronavirus-infection related deep RNA-Sequencing data, we identified 351, 224 and 2,764 circRNAs derived from SARS-CoV-2, SARS-CoV and MERS-CoV, respectively, and characterized two major back-splice events shared by these viruses. Coronavirus-derived circRNAs are more abundant and longer compared to host genome-derived circRNAs. Using a systematic strategy to amplify and identify back-splice junction sequences, we experimentally identified over 100 viral circRNAs from SARS-CoV-2 infected Vero E6 cells. This collection of circRNAs provided the first line of evidence for the abundance and diversity of coronavirus-derived circRNAs and suggested possible mechanisms driving circRNA biogenesis from RNA genomes. Our findings highlight circRNAs as an important component of the coronavirus transcriptome. SummaryWe report for the first time that abundant and diverse circRNAs are generated by SARS-CoV-2, SARS-CoV and MERS-CoV and represent a novel type of circRNAs that differ from circRNAs encoded by DNA genomes.

A Systemic and Molecular Study of Subcellular Localization of SARS-CoV-2 Proteins

Coronavirus possesses the largest RNA genome among all the RNA viruses. Its genome encodes about 29 proteins. Most of the viral proteins are non-structural proteins (NSP) except envelop (E), membrane (M), nucleocapsid (N) and Spike (S) proteins that constitute the viral nucleocapsid, envelop and surface. We have recently cloned all the 29 SARS-CoV-2 genes into vectors for their expressions in mammalian cells except NSP11 that has only 14 amino acids (aa). We are able to express all the 28 cloned SARS-CoV-2 genes in human cells to characterize their subcellular distributions. The proteins of SARS-CoV-2 are mostly cytoplasmic but some are both cytoplasmic and nuclear. Those punctate staining proteins were further investigated by immunofluorescent assay (IFA) using specific antibodies or by co-transfection with an organelle marker-expressing plasmid. As a result, we found that NSP15, ORF6, M and ORF7a are related to Golgi apparatus, and that ORF7b, ORF8 and ORF10 colocalize with endoplasmic reticulum (ER). Interestingly, ORF3a distributes in cell membrane, early endosome, endosome, late endosome and lysosome, which suggests that ORF3a might help the infected virus to usurp endosome and lysosome for viral use. Furthermore, we revealed that NSP13 colocalized with SC35, a protein standing for splicing compartments in the nucleus. Our studies for the first time visualized the subcellular locations of SARS-CoV-2 proteins and might provide novel insights into the viral proteins biological functions.

Systemic Effects of Missense Mutations on SARS-CoV-2 Spike Glycoprotein Stability and Receptor Binding Affinity

The spike (S) glycoprotein of SARS-CoV-2 is responsible for the binding to the permissive cells. The receptor-binding domain (RBD) of SARS-CoV-2 S protein directly interacts with the human angiotensin-converting enzyme 2 (ACE2) on the host cell membrane. In this study, we used computational saturation mutagenesis approaches, including structure-based energy calculations and sequence-based pathogenicity predictions, to quantify the systemic effects of missense mutations on SARS-CoV-2 S protein structure and function. A total of 18,354 mutations in S protein were analyzed and we discovered that most of these mutations could destabilize the entire S protein and its RBD. Specifically, residues G431 and S514 in SARS-CoV-2 RBD are important for S protein stability. We analyzed 384 experimentally verified S missense variations and revealed that the dominant pandemic form, D614G, can stabilize the entire S protein. Moreover, many mutations in N-linked glycosylation sites can increase the stability of the S protein. In addition, we investigated 3,705 mutations in SARS-CoV-2 RBD and 11,324 mutations in human ACE2 and found that SARS-CoV-2 neighbor residues G496 and F497 and ACE2 residues D355 and Y41 are critical for the RBD-ACE2 interaction. The findings comprehensively provide potential target sites in the development of drugs and vaccines against COVID-19.

Protein interactions in the murine cytomegalovirus capsid revealed by cryoEM

Cytomegalovirus (CMV) is distinct among members of the Herpesviridae family for having the largest dsDNA genome (230 kb). Packaging of large dsDNA genome is known to give rise to a highly pressurized viral capsid, but molecular interactions conducive to the formation of CMV capsid resistant to pressurization have not been described. Here, we report a cryo electron microscopy (cryoEM) structure of the murine cytomegalovirus (MCMV) capsid at a 9.1 Å resolution and describe the molecular interactions among the ∼3000 protein molecules in the MCMV capsid at the secondary structure level. Secondary structural elements are resolved to provide landmarks for correlating with results from sequence-based prediction and for structure-based homology modeling. The major capsid protein (MCP) upper domain (MCPud) contains α-helices and β-sheets conserved with those in MCPud of herpes simplex virus type 1 (HSV-1), with the largest differences identified as a "saddle loop" region, located at the tip of MCPud and involved in interaction with the smallest capsid protein (SCP). Interactions among the bacteriophage HK97-like floor domain of MCP, the middle domain of MCP, the hook and clamp domains of the triplex proteins (hoop and clamp domains of TRI-1 and clamp domain of TRI-2) contribute to the formation of a mature capsid. These results offer a framework for understanding how cytomegalovirus uses various secondary structural elements of its capsid proteins to build a robust capsid for packaging its large dsDNA genome inside and for attaching unique functional tegument proteins outside.