RNA-RNA and RNA-protein interactions in coronavirus replication and transcription

Mechanistic insights into the activity of SARS-CoV-2 RNA polymerase inhibitors using single-molecule FRET

The COVID-19 pandemic, caused by the SARS-CoV-2 virus, has resulted in significant global mortality, with over 7 million cases reported. Despite extensive research and high vaccination rates, highly mutated forms of the virus continue to circulate. It is therefore important to understand the viral lifecycle and the precise molecular mechanisms underlying SARS-CoV-2 replication. To address this, we developed a single-molecule Förster resonance energy transfer (smFRET) assay to directly visualize and analyse in vitro RNA synthesis by the SARS-CoV-2 RNA-dependent RNA polymerase (RdRp). We purified the minimal replication complex, comprising nsp12, nsp7, and nsp8, and combined it with fluorescently labelled RNA substrates, enabling real-time monitoring of RNA primer elongation at the single-molecule level. This platform allowed us to investigate the mechanisms of action of key inhibitors of SARS-CoV-2 replication. In particular, our data provides evidence for remdesivir's mechanism of action, which involves polymerase stalling and subsequent chain termination dependent on the concentration of competing nucleotide triphosphates. Our study demonstrates the power of smFRET to provide dynamic insights into SARS-CoV-2 replication, offering a valuable tool for antiviral screening and mechanistic studies of viral RdRp activity.

An intra-family conserved high-order RNA structure within the M ORF is important for arterivirus subgenomic RNA accumulation and infectious virus production

Synthesis of subgenomic RNAs is a strategy commonly used by polycistronic positive-sense single-stranded RNA viruses to express 3'-proximal genes. Members of the order Nidovirales, including coronaviruses and arteriviruses, use a unique discontinuous transcription strategy to synthesize subgenomic RNAs. In this study, in silico synonymous site conservation analysis and RNA structure folding predicted the existence of intra-family conserved high-order RNA structure within the M ORF of arteriviral genomes, which was further confirmed by RNA secondary structure probing. This RNA structure was determined to be important for the transcription/accumulation of subgenomic RNAs and the production of infectious viral particles. Mutations disrupting the stability of the RNA structures significantly decreased the accumulation of multiple subgenomic RNAs. In contrast, the impact of mutagenesis on full-length genomic RNA accumulation was limited. The degree to which wild-type levels of subgenomic RNA accumulation were maintained was found to correlate with the efficiency of infectious virus production. Moreover, the thermo-stability of stems within the high-order RNA structure is also well correlated with viral replication capacity and the maintenance of subgenomic RNA accumulation. This study is the first to report an intra-Arteriviridae conserved high-order RNA structure that is located in a protein-coding region and functions as an important cis-acting element to control the accumulation/transcription of arteriviral subgenomic RNAs. This work suggests a complex regulation mechanism between genome replication and discontinuous transcription in nidoviruses.IMPORTANCEArteriviruses are a group of RNA viruses that infect different animal species. They can cause diseases associated with respiratory/reproductive syndromes, abortion, or hemorrhagic fever. Among arteriviruses, porcine reproductive and respiratory syndrome virus (PRRSV) and equine arteritis virus (EAV) are economically important veterinary pathogens. The challenge in control of arterivirus infection reflects our limited knowledge of viral biology. In this study, we conducted a comprehensive analysis of arteriviral genomes and discovered intra-family conserved regions in the M ORF with a high-order RNA structure. The thermo-stability of the RNA structure influences sgRNA transcription/accumulation and correlates with the level of infectious virus production. Our studies provide new insight into arterivirus replication mechanisms, which may have implications for developing disease control and prevention strategies.

Structural insights into the RNA binding inhibitors of the C-terminal domain of the SARS-CoV-2 nucleocapsid

The SARS-CoV-2 nucleocapsid (N) protein is an essential structural element of the virion, playing a crucial role in enclosing the viral genome into a ribonucleoprotein (RNP) assembly, as well as viral replication and transmission. The C-terminal domain of the N-protein (N-CTD) is essential for encapsidation, contributing to the stabilization of the RNP complex. In a previous study, three inhibitors (ceftriaxone, cefuroxime, and ampicillin) were screened for their potential to disrupt the RNA packaging process by targeting the N-protein. However, the binding efficacy, mechanism of RNA binding inhibition, and molecular insights of binding with N-CTD remain unclear. In this study, we evaluated the binding efficacy of these inhibitors using isothermal titration calorimetry (ITC), revealing the affinity of ceftriaxone (18 ± 1.3 μM), cefuroxime (55 ± 4.2 μM), and ampicillin (28 ± 1.2 μM) with the N-CTD. Further inhibition assay and fluorescence polarisation assay demonstrated RNA binding inhibition, with IC50 ranging from ∼ 12 to 18 μM and KD values between 24 μM to 32 μM for the inhibitors, respectively. Additionally, we also determined the inhibitor-bound complex crystal structures of N-CTD-Ceftriaxone (2.0 Å) and N-CTD-Ampicillin (2.2 Å), along with the structure of apo N-CTD (1.4 Å). These crystal structures revealed previously unobserved interaction sites involving residues K261, K266, R293, Q294, and W301 at the oligomerization interface and the predicted RNA-binding region of N-CTD. These findings provide valuable molecular insights into the inhibition of N-CTD, highlighting its potential as an underexplored but promising target for the development of novel antiviral agents against coronaviruses.

Transcription Kinetics in the Coronavirus Life Cycle

Coronaviruses utilize a positive-sense single-strand RNA, functioning simultaneously as mRNA and the genome. An RNA-dependent RNA polymerase (RdRP) plays a dual role in transcribing genes and replicating the genome, making RdRP a critical target in therapies against coronaviruses. This review explores recent advancements in understanding the coronavirus transcription machinery, discusses it within virus infection context, and incorporates kinetic considerations on RdRP activity. We also address steric limitations in coronavirus replication, particularly during early infection phases, and outline hypothesis regarding translation-transcription conflicts, postulating the existence of mechanisms that resolve these issues. In cells infected by coronaviruses, abundant structural proteins are synthesized from subgenomic RNA fragments (sgRNAs) produced via discontinuous transcription. During elongation, RdRP can skip large sections of the viral genome, resulting in the creation of shorter sgRNAs that reflects the stoichiometry of viral structural proteins. Although the precise mechanism of discontinuous transcription remains unknown, we discuss recent hypotheses involving long-distance RNA-RNA interactions, helicase-mediated RdRP backtracking, dissociation and reassociation of RdRP, and RdRP dimerization.

Emergence of SARS-CoV-2 subgenomic RNAs that enhance viral fitness and immune evasion

Coronaviruses express their structural and accessory genes via a set of subgenomic RNAs, whose synthesis is directed by transcription regulatory sequences (TRSs) in the 5' genomic leader and upstream of each body open reading frame. In SARS-CoV-2, the TRS has the consensus AAACGAAC; upon searching for emergence of this motif in the global SARS-CoV-2 sequences, we find that it evolves frequently, especially in the 3' end of the genome. We show well-supported examples upstream of the Spike gene-within the nsp16 coding region of ORF1b-which is expressed during human infection, and upstream of the canonical Envelope gene TRS, both of which have evolved convergently in multiple lineages. The most frequent neo-TRS is within the coding region of the Nucleocapsid gene, and is present in virtually all viruses from the B.1.1 lineage, including the variants of concern Alpha, Gamma, Omicron and descendants thereof. Here, we demonstrate that this TRS leads to the expression of a novel subgenomic mRNA encoding a truncated C-terminal portion of Nucleocapsid, which is an antagonist of type I interferon production and contributes to viral fitness during infection. We observe distinct phenotypes when the Nucleocapsid coding sequence is mutated compared to when the TRS alone is ablated. Our findings demonstrate that SARS-CoV-2 is undergoing evolutionary changes at the functional RNA level in addition to the amino acid level.

Dissecting the Conformational Heterogeneity of Stem-Loop Substructures of the Fifth Element in the 5'-Untranslated Region of SARS-CoV-2

Throughout the family of coronaviruses, structured RNA elements within the 5' region of the genome are highly conserved. The fifth stem-loop element from SARS-CoV-2 (5_SL5) represents an example of an RNA structural element, repeatedly occurring in coronaviruses. It contains a conserved, repetitive fold within its substructures SL5a and SL5b. We herein report the detailed characterization of the structure and dynamics of elements SL5a and SL5b that are located immediately upstream of the SARS-CoV-2 ORF1a/b start codon. Exploiting the unique ability of solution NMR methods, we show that the structures of both apical loops are modulated by structural differences in the remote parts located in their stem regions. We further integrated our high-resolution models of SL5a/b into the context of full-length 5_SL5 structures by combining different structural biology methods. Finally, we evaluated the impact of the two most common VoC mutations within 5_SL5 with respect to individual base-pair stability.

Regulation viral RNA transcription and replication by higher-order RNA structures within the nsp1 coding region of MERS coronavirus

Coronavirus (CoV) possesses numerous functional cis-acting elements in its positive-strand genomic RNA. Although most of these RNA structures participate in viral replication, the functions of RNA structures in the genomic RNA of CoV in viral replication remain unclear. In this study, we investigated the functions of the higher-order RNA stem-loop (SL) structures SL5B, SL5C, and SL5D in the ORF1a coding region of Middle East respiratory syndrome coronavirus (MERS-CoV) in viral replication. Our approach, using reverse genetics of a bacterial artificial chromosome system, revealed that SL5B and SL5C play essential roles in the discontinuous transcription of MERS-CoV. In silico analyses predicted that SL5C interacts with a bulged stem-loop (BSL) in the 3' untranslated region, suggesting that the RNA structure of SL5C is important for viral RNA transcription. Conversely, SL5D did not affect transcription, but mediated the synthesis of positive-strand genomic RNA. Additionally, the RNA secondary structure of SL5 in the revertant virus of the SL5D mutant was similar to that of the wild-type, indicating that the RNA structure of SL5D can finely tune RNA replication in MERS-CoV. Our data indicate novel regulatory mechanisms of viral RNA transcription and replication by higher-order RNA structures in the MERS-CoV genomic RNA.

m6A Methylation of Transcription Leader Sequence of SARS-CoV-2 Impacts Discontinuous Transcription of Subgenomic mRNAs

The SARS-CoV-2 genome has been shown to be m6A methylated at several positions in vivo. Strikingly, a DRACH motif, the recognition motif for adenosine methylation, resides in the core of the transcriptional regulatory leader sequence (TRS-L) at position A74, which is highly conserved and essential for viral discontinuous transcription. Methylation at position A74 correlates with viral pathogenicity. Discontinuous transcription produces a set of subgenomic mRNAs that function as templates for translation of all structural and accessory proteins. A74 is base-paired in the short stem-loop structure 5'SL3 that opens during discontinuous transcription to form long-range RNA-RNA interactions with nascent (-)-strand transcripts at complementary TRS-body sequences. A74 can be methylated by the human METTL3/METTL14 complex in vitro. Here, we investigate its impact on the structural stability of 5'SL3 and the long-range TRS-leader:TRS-body duplex formation necessary for synthesis of subgenomic mRNAs of all four viral structural proteins. Methylation uniformly destabilizes 5'SL3 and long-range duplexes and alters their relative equilibrium populations, suggesting that the m6A74 modification acts as a regulator for the abundance of viral structural proteins due to this destabilization.

Assembly of SARS-CoV-2 nucleocapsid protein with nucleic acid

The viral genome of SARS-CoV-2 is packaged by the nucleocapsid (N-)protein into ribonucleoprotein particles (RNPs), 38 ± 10 of which are contained in each virion. Their architecture has remained unclear due to the pleomorphism of RNPs, the high flexibility of N-protein intrinsically disordered regions, and highly multivalent interactions between viral RNA and N-protein binding sites in both N-terminal (NTD) and C-terminal domain (CTD). Here we explore critical interaction motifs of RNPs by applying a combination of biophysical techniques to ancestral and mutant proteins binding different nucleic acids in an in vitro assay for RNP formation, and by examining nucleocapsid protein variants in a viral assembly assay. We find that nucleic acid-bound N-protein dimers oligomerize via a recently described protein-protein interface presented by a transient helix in its long disordered linker region between NTD and CTD. The resulting hexameric complexes are stabilized by multivalent protein-nucleic acid interactions that establish crosslinks between dimeric subunits. Assemblies are stabilized by the dimeric CTD of N-protein offering more than one binding site for stem-loop RNA. Our study suggests a model for RNP assembly where N-protein scaffolding at high density on viral RNA is followed by cooperative multimerization through protein-protein interactions in the disordered linker.

Bioinformatics Insights on Viral Gene Expression Transactivation: From HIV-1 to SARS-CoV-2

Viruses provide vital insights into gene expression control. Viral transactivators, with other viral and cellular proteins, regulate expression of self, other viruses, and host genes with profound effects on infected cells, underlying inflammation, control of immune responses, and pathogenesis. The multifunctional Tat proteins of lentiviruses (HIV-1, HIV-2, and SIV) transactivate gene expression by recruiting host proteins and binding to transacting responsive regions (TARs) in viral and host RNAs. SARS-CoV-2 nucleocapsid participates in early viral transcription, recruits similar cellular proteins, and shares intracellular, surface, and extracellular distribution with Tat. SARS-CoV-2 nucleocapsid interacting with the replication-transcription complex might, therefore, transactivate viral and cellular RNAs in the transcription and reactivation of self and other viruses, acute and chronic pathogenesis, immune evasion, and viral evolution. Here, we show, by using primary and secondary structural comparisons, that the leaders of SARS-CoV-2 and other coronaviruses contain TAR-like sequences in stem-loops 2 and 3. The coronaviral nucleocapsid C-terminal domains harbor a region of similarity to TAR-binding regions of lentiviral Tat proteins, and coronaviral nonstructural protein 12 has a cysteine-rich metal binding, dimerization domain, as do lentiviral Tat proteins. Although SARS-CoV-1 nucleocapsid transactivated gene expression in a replicon-based study, further experimental evidence for coronaviral transactivation and its possible implications is warranted.

Classification, replication, and transcription of Nidovirales

Nidovirales is one order of RNA virus, with the largest single-stranded positive sense RNA genome enwrapped with membrane envelope. It comprises four families (Arterividae, Mesoniviridae, Roniviridae, and Coronaviridae) and has been circulating in humans and animals for almost one century, posing great threat to livestock and poultry，as well as to public health. Nidovirales shares similar life cycle: attachment to cell surface, entry, primary translation of replicases, viral RNA replication in cytoplasm, translation of viral proteins, virion assembly, budding, and release. The viral RNA synthesis is the critical step during infection, including genomic RNA (gRNA) replication and subgenomic mRNAs (sg mRNAs) transcription. gRNA replication requires the synthesis of a negative sense full-length RNA intermediate, while the sg mRNAs transcription involves the synthesis of a nested set of negative sense subgenomic intermediates by a discontinuous strategy. This RNA synthesis process is mediated by the viral replication/transcription complex (RTC), which consists of several enzymatic replicases derived from the polyprotein 1a and polyprotein 1ab and several cellular proteins. These replicases and host factors represent the optimal potential therapeutic targets. Hereby, we summarize the Nidovirales classification, associated diseases, "replication organelle," replication and transcription mechanisms, as well as related regulatory factors.

Identification of subgenomic mRNAs derived from the coronavirus 1a/1b protein gene: Implications for coronavirus transcription

Synthesis of coronavirus subgenomic mRNA (sgmRNA) is guided by the transcription regulatory sequence (TRS). sgmRNA derived from the body TRS (TRS-B) located at the 1a/1b protein gene is designated 1ab/sgmRNA. In the current study, we comprehensively identified the 1ab/sgmRNAs synthesized from TRS-Bs located at the 1a/1b protein genes of different coronavirus genera both in vitro and in vivo by RT‒PCR and sequencing. The results suggested that the degree of sequence homology between the leader TRS (TRS-L) and TRS-B may not be a decisive factor for 1ab/sgmRNA synthesis. This observation led us to revisit the coronavirus transcription mechanism and to propose that the disassociation of coronavirus polymerase from the viral genome may be a prerequisite for sgmRNA synthesis. Once the polymerase can disassociate at TRS-B, the sequence homology between TRS-L and TRS-B is important for sgmRNA synthesis. The study therefore extends our understanding of transcription mechanisms.

Prediction of LncRNA-protein Interactions Using Auto-Encoder, SE-ResNet Models and Transfer Learning

BACKGROUND

Long non-coding RNA (lncRNA) plays a crucial role in various biological processes, and mutations or imbalances of lncRNAs can lead to several diseases, including cancer, Prader-Willi syndrome, autism, Alzheimer's disease, cartilage-hair hypoplasia, and hearing loss. Understanding lncRNA-protein interactions (LPIs) is vital for elucidating basic cellular processes, human diseases, viral replication, transcription, and plant pathogen resistance. Despite the development of several LPI calculation methods, predicting LPI remains challenging, with the selection of variables and deep learning structure being the focus of LPI research.

METHODS

We propose a deep learning framework called AR-LPI, which extracts sequence and secondary structure features of proteins and lncRNAs. The framework utilizes an auto-encoder for feature extraction and employs SE-ResNet for prediction. Additionally, we apply transfer learning to the deep neural network SE-ResNet for predicting small-sample datasets.

RESULTS

Through comprehensive experimental comparison, we demonstrate that the AR-LPI architecture performs better in LPI prediction. Specifically, the accuracy of AR-LPI increases by 2.86% to 94.52%, while the F-value of AR-LPI increases by 2.71% to 94.73%.

CONCLUSION

Our experimental results show that the overall performance of AR-LPI is better than that of other LPI prediction tools.

3'UTR of ALV-J can affect viral replication through promoting transcription and mRNA nuclear export

IMPORTANCE

3'UTRs can affect gene transcription and post-transcriptional regulation in multiple ways, further influencing the function of proteins in a unique manner. Recently, ALV-J has been mutating and evolving rapidly, especially the 3'UTR of viral genome. Meanwhile, clinical symptoms caused by ALV-J have changed significantly. In this study, we found that the ALV-J strains containing △-r-TM-type 3'UTR are the most abundant. By constructing ALV-J infectious clones and subgenomic vectors containing different 3'UTRs, we prove that 3'UTRs directly affect viral tissue preference and can promote virus replication as an enhancer. ALV-J strain containing 3'UTR of △-r-TM proliferated fastest in primary cells. All five forms of 3'UTRs can assist intron-containing viral mRNA nuclear export, with similar efficiency. ALV-J mRNA half-life is not influenced by different 3'UTRs. Our results dissect the roles of 3'UTR on regulating viral replication and pathogenicity, providing novel insights into potential anti-viral strategies.

RNA structure and multiple weak interactions balance the interplay between RNA binding and phase separation of SARS-CoV-2 nucleocapsid

The nucleocapsid (N) protein of SARS-CoV-2 binds viral RNA, condensing it inside the virion, and phase separating with RNA to form liquid-liquid condensates. There is little consensus on what differentiates sequence-independent N-RNA interactions in the virion or in liquid droplets from those with specific genomic RNA (gRNA) motifs necessary for viral function inside infected cells. To identify the RNA structures and the N domains responsible for specific interactions and phase separation, we use the first 1,000 nt of viral RNA and short RNA segments designed as models for single-stranded and paired RNA. Binding affinities estimated from fluorescence anisotropy of these RNAs to the two-folded domains of N (the NTD and CTD) and comparison to full-length N demonstrate that the NTD binds preferentially to single-stranded RNA, and while it is the primary RNA-binding site, it is not essential to phase separation. Nuclear magnetic resonance spectroscopy identifies two RNA-binding sites on the NTD: a previously characterized site and an additional although weaker RNA-binding face that becomes prominent when binding to the primary site is weak, such as with dsRNA or a binding-impaired mutant. Phase separation assays of nucleocapsid domains with double-stranded and single-stranded RNA structures support a model where multiple weak interactions, such as with the CTD or the NTD's secondary face promote phase separation, while strong, specific interactions do not. These studies indicate that both strong and multivalent weak N-RNA interactions underlie the multifunctional abilities of N.

Inhibiting mechanism of Alpiniae oxyphyllae fructus polysaccharide 3 against the replication of porcine epidemic diarrhea virus

Porcine epidemic diarrhea virus (PEDV) causes diarrhea, vomiting, and death in piglets. Our previous study has revealed the anti-PEDV activity of Alpiniae oxyphyllae fructus polysaccharide 3 (AOFP3). However, it is still unknown whether AOFP3 can inhibit the replication of PEDV. Therefore, the effect of AOFP3 on PEDV replication was investigated in the present study, along with analysis of viral RdRp activity and expression of hnRNP A1 by RNA polymerase activity assay in vitro, RIP assay, and Western blotting. The results showed that both the PEDV gene and protein levels in IPEC-J2 cells decreased with AOFP3 treatment. In addition, AOFP3 significantly reduced PEDV's replication by down-regulating the activity of PEDV RdRp and reducing the expression of hnRNP A1, whereas only the bind of RdRp to PEDV 3'UTR was inhibited in AOFP3 treated cells.

The International Virus Bioinformatics Meeting 2023

The 2023 International Virus Bioinformatics Meeting was held in Valencia, Spain, from 24-26 May 2023, attracting approximately 180 participants worldwide. The primary objective of the conference was to establish a dynamic scientific environment conducive to discussion, collaboration, and the generation of novel research ideas. As the first in-person event following the SARS-CoV-2 pandemic, the meeting facilitated highly interactive exchanges among attendees. It served as a pivotal gathering for gaining insights into the current status of virus bioinformatics research and engaging with leading researchers and emerging scientists. The event comprised eight invited talks, 19 contributed talks, and 74 poster presentations across eleven sessions spanning three days. Topics covered included machine learning, bacteriophages, virus discovery, virus classification, virus visualization, viral infection, viromics, molecular epidemiology, phylodynamic analysis, RNA viruses, viral sequence analysis, viral surveillance, and metagenomics. This report provides rewritten abstracts of the presentations, a summary of the key research findings, and highlights shared during the meeting.

Recombination between coronaviruses and synthetic RNAs and biorisk implications motivated by a SARS-CoV-2 FCS origin controversy

The urgent need for improved policy, regulation, and oversight of research with potential pandemic pathogens (PPPs) has been widely acknowledged. A 2022 article in Frontiers in Virology raises questions, reporting on a 100% sequence homology between the SARS-CoV-2 furin cleavage site (FCS) and the negative strand of a 2017 patented sequence. Even though Ambati and collaborators suspect a possible inadvertent or intentional cause leading to the FCS insert, the related underpinnings have not been studied from the perspective of potential biorisk policy gaps. A commentary on their article contests the low coincidence likelihood that was calculated by Ambati et al., arguing that the sequence match could have been a chance occurrence alone. Additionally, it has been suggested that the odds of the recombination event may be low. These considerations seem to have put many speculations related to any implied viral beginnings, notably from a research setting likely outside the Wuhan Institute of Virology, to rest. However, potential implications for future disasters in terms of biosafety and biosecurity have not been addressed. To demonstrate the feasibility of the Ambati et al. postulate, a theoretical framework is developed that substantially extends the research orientations implicated by these authors and the related patent. It is argued that specific experimental conditions, in combination, could significantly increase the implied recombination profile between coronaviruses and synthetic RNAs. Consequently, this article scrutinizes these largely unrecognized vulnerabilities to discuss implications across the spectrum of the biological risk landscape, with special attention to a potential "crime harvest." Focusing on insufficiently understood features of interaction between the natural and man-made world, vulnerabilities related to contaminants, camouflaging, and various misuse potentials fostered by the digitization and computerization of synthetic biology, it highlights novel biorisk gaps not covered by existing PPP policy. Even though this work does not aim to provide proof of the viral origin, it will make the point that, in theory, a convergence of under-appreciated lab experiments and technologies could have led to the SARS-CoV-2 FCS insert, which analogously could be exploited by various threat actors for the clandestine genesis of similar or even worse pathogens.

Mutagenesis and structural studies reveal the basis for the specific binding of SARS-CoV-2 SL3 RNA element with human TIA1 protein

Viral RNA-host protein interactions are indispensable during RNA virus transcription and replication, but their detailed structural and dynamical features remain largely elusive. Here, we characterize the binding interface for the SARS-CoV-2 stem-loop 3 (SL3) cis-acting element to human TIA1 protein with a combined theoretical and experimental approaches. The highly structured SARS-CoV-2 SL3 has a high binding affinity to TIA1 protein, in which the aromatic stacking, hydrogen bonds, and hydrophobic interactions collectively direct this specific binding. Further mutagenesis studies validate our proposed 3D binding model and reveal two SL3 variants have enhanced binding affinities to TIA1. And disruptions of the identified RNA-protein interactions with designed antisense oligonucleotides dramatically reduce SARS-CoV-2 infection in cells. Finally, TIA1 protein could interact with conserved SL3 RNA elements within other betacoronavirus lineages. These findings open an avenue to explore the viral RNA-host protein interactions and provide a pioneering structural basis for RNA-targeting antiviral drug design.

Metatranscriptomics analysis reveals a novel transcriptional and translational landscape during Middle East respiratory syndrome coronavirus infection

Among all RNA viruses, coronavirus RNA transcription is the most complex and involves a process termed "discontinuous transcription" that results in the production of a set of 3'-nested, co-terminal genomic and subgenomic RNAs during infection. While the expression of the classic canonical set of subgenomic RNAs depends on the recognition of a 6- to 7-nt transcription regulatory core sequence (TRS), here, we use deep sequence and metagenomics analysis strategies and show that the coronavirus transcriptome is even more vast and more complex than previously appreciated and involves the production of leader-containing transcripts that have canonical and noncanonical leader-body junctions. Moreover, by ribosome protection and proteomics analyses, we show that both positive- and negative-sense transcripts are translationally active. The data support the hypothesis that the coronavirus proteome is much vaster than previously noted in the literature.

The preference signature of the SARS-CoV-2 Nucleocapsid NTD for its 5'-genomic RNA elements

The nucleocapsid protein (N) of SARS-CoV-2 plays a pivotal role during the viral life cycle. It is involved in RNA transcription and accounts for packaging of the large genome into virus particles. N manages the enigmatic balance of bulk RNA-coating versus precise RNA-binding to designated cis-regulatory elements. Numerous studies report the involvement of its disordered segments in non-selective RNA-recognition, but how N organizes the inevitable recognition of specific motifs remains unanswered. We here use NMR spectroscopy to systematically analyze the interactions of N's N-terminal RNA-binding domain (NTD) with individual cis RNA elements clustering in the SARS-CoV-2 regulatory 5'-genomic end. Supported by broad solution-based biophysical data, we unravel the NTD RNA-binding preferences in the natural genome context. We show that the domain's flexible regions read the intrinsic signature of preferred RNA elements for selective and stable complex formation within the large pool of available motifs.

Host microRNA interactions with the SARS-CoV-2 viral genome 3'-untranslated region

The 2019 pandemic, caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has marked the spread of a novel human coronavirus. While the viral life cycle is well understood, most of the interactions at the virus-host interface remain elusive. Furthermore, the molecular mechanisms behind disease severity and immune evasion are still largely unknown. Conserved elements of the viral genome such as secondary structures within the 5'- and 3'-untranslated regions (UTRs) serve as attractive targets of interest and could prove crucial in furthering our understanding of virus-host interactions. It has been proposed that microRNA (miR) interactions with viral components could be used by both the virus and host for their own benefit. Analysis of the SARS-CoV-2 viral genome 3'-UTR has revealed the potential for host cellular miR binding sites, providing sites for specific interactions with the virus. In this study, we demonstrate that the SARS-CoV-2 genome 3'-UTR binds the host cellular miRNAs miR-760-3p, miR-34a-5p, and miR-34b-5p, which have been shown to influence translation of interleukin-6 (IL-6), the IL-6 receptor (IL-6R), as well as progranulin (PGRN), respectively, proteins that have roles in the host immune response and inflammatory pathways. Furthermore, recent work suggests the potential of miR-34a-5p and miR-34b-5p to target and inhibit translation of viral proteins. Native gel electrophoresis and steady-state fluorescence spectroscopy were utilized to characterize the binding of these miRs to their predicted sites within the SARS-CoV-2 genome 3'-UTR. Additionally, we investigated 2'-fluoro-D-arabinonucleic acid (FANA) analogs of these miRNAs as competitive binding inhibitors for these miR binding interactions. The mechanisms detailed in this study have the potential to drive the development of antiviral treatments for SARS-CoV-2 infection, and provide a potential molecular basis for cytokine release syndrome and immune evasion which could implicate the host-virus interface.

The 5'UTR of HCoV-OC43 adopts a topologically constrained structure to intrinsically repress translation

The emergence of SARS-CoV-2, which is responsible for the COVID-19 pandemic, has highlighted the need for rapid characterization of viral mechanisms associated with cellular pathogenesis. Viral UTRs represent conserved genomic elements that contribute to such mechanisms. Structural details of most CoV UTRs are not available, however. Experimental approaches are needed to allow for the facile generation of high-quality viral RNA tertiary structural models, which can facilitate comparative mechanistic efforts. By integrating experimental and computational techniques, we herein report the efficient characterization of conserved RNA structures within the 5'UTR of the HCoV-OC43 genome, a lab-tractable model coronavirus. We provide evidence that the 5'UTR folds into a structure with well-defined stem-loops (SLs) as determined by chemical probing and direct detection of hydrogen bonds by NMR. We combine experimental base-pair restraints with global structural information from SAXS to generate a 3D model that reveals that SL1-4 adopts a topologically constrained structure wherein SLs 3 and 4 coaxially stack. Coaxial stacking is mediated by short linker nucleotides and allows SLs 1 to 2 to sample different cojoint orientations by pivoting about the SL3,4 helical axis. To evaluate the functional relevance of the SL3,4 coaxial helix, we engineered luciferase reporter constructs harboring the HCoV-OC43 5'UTR with mutations designed to abrogate coaxial stacking. Our results reveal that the SL3,4 helix intrinsically represses translation efficiency since the destabilizing mutations correlate with increased luciferase expression relative to wildtype without affecting reporter mRNA levels, thus highlighting how the 5'UTR structure contributes to the viral mechanism.

Identification of host factors that bind to the 5′ end of the MERS-CoV RNA genome

Aim: The aim of this study was to identify host factors that interact with the 5′ end of the MERS-CoV RNA genome. Materials & methods: RNA affinity chromatography followed by mass spectrometry analysis was used to identify the binding of host factors in Vero E6 cells. Results: A total of 59 host factors that bound the MERS-CoV RNA genome in non-infected Vero E6 cells were identified. Most of the identified cellular proteins were previously reported to interact with the genome of other RNA viruses. We validated our mass spectrometry results using western blotting. Conclusion: These data enhance our knowledge about the RNA–host interactions of coronaviruses, which could serve as targets for developing antiviral therapeutics against MERS-CoV.

Deletion of the s2m RNA Structure in the Avian Coronavirus Infectious Bronchitis Virus and Human Astrovirus Results in Sequence Insertions

Coronaviruses infect a wide variety of host species, resulting in a range of diseases in both humans and animals. The coronavirus genome consists of a large positive-sense single-stranded molecule of RNA containing many RNA structures. One structure, denoted s2m and consisting of 41 nucleotides, is located within the 3' untranslated region (3' UTR) and is shared between some coronavirus species, including infectious bronchitis virus (IBV), severe acute respiratory syndrome coronavirus (SARS-CoV), and SARS-CoV-2, as well as other pathogens, including human astrovirus. Using a reverse genetic system to generate recombinant viruses, we investigated the requirement of the s2m structure in the replication of IBV, a globally distributed economically important Gammacoronavirus that infects poultry causing respiratory disease. Deletion of three nucleotides predicted to destabilize the canonical structure of the s2m or the deletion of the nucleotides corresponding to s2m impacted viral replication in vitro. In vitro passaging of the recombinant IBV with the s2m sequence deleted resulted in a 36-nucleotide insertion in place of the deletion, which was identified to be composed of a duplication of flanking sequences. A similar result was observed following serial passage of human astrovirus with a deleted s2m sequence. RNA modeling indicated that deletion of the nucleotides corresponding to the s2m impacted other RNA structures present in the IBV 3' UTR. Our results indicated for both IBV and human astrovirus a preference for nucleotide occupation in the genome location corresponding to the s2m, which is independent of the specific s2m sequence. IMPORTANCE Coronaviruses infect many species, including humans and animals, with substantial effects on livestock, particularly with respect to poultry. The coronavirus RNA genome consists of structural elements involved in viral replication whose roles are poorly understood. We investigated the requirement of the RNA structural element s2m in the replication of the Gammacoronavirus infectious bronchitis virus, an economically important viral pathogen of poultry. Using reverse genetics to generate recombinant IBVs with either a disrupted or deleted s2m, we showed that the s2m is not required for viral replication in cell culture; however, replication is decreased in tracheal tissue, suggesting a role for the s2m in the natural host. Passaging of these viruses as well as human astrovirus lacking the s2m sequence demonstrated a preference for nucleotide occupation, independent of the s2m sequence. RNA modeling suggested deletion of the s2m may negatively impact other essential RNA structures.

Genome Structure, Life Cycle, and Taxonomy of Coronaviruses and the Evolution of SARS-CoV-2

Coronaviruses have a broad host range and exhibit high zoonotic potential. In this chapter, we describe their genomic organization in terms of encoded proteins and provide an introduction to the peculiar discontinuous transcription mechanism. Further, we present evolutionary conserved genomic RNA secondary structure features, which are involved in the complex replication mechanism. With a focus on computational methods, we review the emergence of SARS-CoV-2 starting with the 2019 strains. In that context, we also discuss the debated hypothesis of whether SARS-CoV-2 was created in a laboratory. We focus on the molecular evolution and the epidemiological dynamics of this recently emerged pathogen and we explain how variants of concern are detected and characterised. COVID-19, the disease caused by SARS-CoV-2, can spread through different transmission routes and also depends on a number of risk factors. We describe how current computational models of viral epidemiology, or more specifically, phylodynamics, have facilitated and will continue to enable a better understanding of the epidemic dynamics of SARS-CoV-2.

Antiviral Activity of Oligonucleotides Targeting the SARS-CoV-2 Genomic RNA Stem-Loop Sequences within the 3'-End of the ORF1b

Increased evidence shows vaccines against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) exhibited no long-term efficacy and limited worldwide availability, while existing antivirals and treatment options have only limited efficacy. In this study, the main objective was the development of antiviral strategies using nucleic acid-based molecules. To this purpose, partially overlapped 6-19-mer phosphorothioate deoxyoligonucleotides (S-ONs) designed on the SARS-CoV-2 genomic RNA stem-loop packaging sequences within the 3' end of the ORF1b were synthetized using the direct and complementary sequence. Among the S-ONs tested, several oligonucleotides exhibited a fifty percent inhibitory concentration antiviral activity ranging from 0.27 to 34 μM, in the absence of cytotoxicity. The S-ON with a scrambled sequence used in the same conditions was not active. Moreover, selected 10-mer S-ONs were tested using different infectious doses and against different SARS-CoV-2 variants, showing comparable antiviral activity that was abrogated when the central sequence was mutated. Experiments to evaluate the intracellular functional target localization of the S-ON inhibitory activity were also performed. Collectively the data indicate that the SARS-CoV-2 packaging region in the 3' end of the ORF1b may be a promising target candidate for further investigation to develop innovative nucleic-acid-based antiviral therapy.

CovInter: interaction data between coronavirus RNAs and host proteins

Coronavirus has brought about three massive outbreaks in the past two decades. Each step of its life cycle invariably depends on the interactions among virus and host molecules. The interaction between virus RNA and host protein (IVRHP) is unique compared to other virus-host molecular interactions and represents not only an attempt by viruses to promote their translation/replication, but also the host's endeavor to combat viral pathogenicity. In other words, there is an urgent need to develop a database for providing such IVRHP data. In this study, a new database was therefore constructed to describe the interactions between coronavirus RNAs and host proteins (CovInter). This database is unique in (a) unambiguously characterizing the interactions between virus RNA and host protein, (b) comprehensively providing experimentally validated biological function for hundreds of host proteins key in viral infection and (c) systematically quantifying the differential expression patterns (before and after infection) of these key proteins. Given the devastating and persistent threat of coronaviruses, CovInter is highly expected to fill the gap in the whole process of the 'molecular arms race' between viruses and their hosts, which will then aid in the discovery of new antiviral therapies. It's now free and publicly accessible at: https://idrblab.org/covinter/.

Nanopore ReCappable sequencing maps SARS-CoV-2 5' capping sites and provides new insights into the structure of sgRNAs

The SARS-CoV-2 virus has a complex transcriptome characterised by multiple, nested subgenomic RNAsused to express structural and accessory proteins. Long-read sequencing technologies such as nanopore direct RNA sequencing can recover full-length transcripts, greatly simplifying the assembly of structurally complex RNAs. However, these techniques do not detect the 5' cap, thus preventing reliable identification and quantification of full-length, coding transcript models. Here we used Nanopore ReCappable Sequencing (NRCeq), a new technique that can identify capped full-length RNAs, to assemble a complete annotation of SARS-CoV-2 sgRNAs and annotate the location of capping sites across the viral genome. We obtained robust estimates of sgRNA expression across cell lines and viral isolates and identified novel canonical and non-canonical sgRNAs, including one that uses a previously un-annotated leader-to-body junction site. The data generated in this work constitute a useful resource for the scientific community and provide important insights into the mechanisms that regulate the transcription of SARS-CoV-2 sgRNAs.

The Neighborhood of the Spike Gene Is a Hotspot for Modular Intertypic Homologous and Nonhomologous Recombination in Coronavirus Genomes

Coronaviruses (CoVs) have very large RNA viral genomes with a distinct genomic architecture of core and accessory open reading frames (ORFs). It is of utmost importance to understand their patterns and limits of homologous and nonhomologous recombination, because such events may affect the emergence of novel CoV strains, alter their host range, infection rate, tissue tropism pathogenicity, and their ability to escape vaccination programs. Intratypic recombination among closely related CoVs of the same subgenus has often been reported; however, the patterns and limits of genomic exchange between more distantly related CoV lineages (intertypic recombination) need further investigation. Here, we report computational/evolutionary analyses that clearly demonstrate a substantial ability for CoVs of different subgenera to recombine. Furthermore, we show that CoVs can obtain-through nonhomologous recombination-accessory ORFs from core ORFs, exchange accessory ORFs with different CoV genera, with other viruses (i.e., toroviruses, influenza C/D, reoviruses, rotaviruses, astroviruses) and even with hosts. Intriguingly, most of these radical events result from double crossovers surrounding the Spike ORF, thus highlighting both the instability and mobile nature of this genomic region. Although many such events have often occurred during the evolution of various CoVs, the genomic architecture of the relatively young SARS-CoV/SARS-CoV-2 lineage so far appears to be stable.

The Remarkable Evolutionary Plasticity of Coronaviruses by Mutation and Recombination: Insights for the COVID-19 Pandemic and the Future Evolutionary Paths of SARS-CoV-2

Coronaviruses (CoVs) constitute a large and diverse subfamily of positive-sense single-stranded RNA viruses. They are found in many mammals and birds and have great importance for the health of humans and farm animals. The current SARS-CoV-2 pandemic, as well as many previous epidemics in humans that were of zoonotic origin, highlights the importance of studying the evolution of the entire CoV subfamily in order to understand how novel strains emerge and which molecular processes affect their adaptation, transmissibility, host/tissue tropism, and patho non-homologous genicity. In this review, we focus on studies over the last two years that reveal the impact of point mutations, insertions/deletions, and intratypic/intertypic homologous and non-homologous recombination events on the evolution of CoVs. We discuss whether the next generations of CoV vaccines should be directed against other CoV proteins in addition to or instead of spike. Based on the observed patterns of molecular evolution for the entire subfamily, we discuss five scenarios for the future evolutionary path of SARS-CoV-2 and the COVID-19 pandemic. Finally, within this evolutionary context, we discuss the recently emerged Omicron (B.1.1.529) VoC.

Amilorides inhibit SARS-CoV-2 replication in vitro by targeting RNA structures

The SARS-CoV-2 pandemic, and the likelihood of future coronavirus pandemics, emphasized the urgent need for development of novel antivirals. Small-molecule chemical probes offer both to reveal aspects of virus replication and to serve as leads for antiviral therapeutic development. Here, we report on the identification of amiloride-based small molecules that potently inhibit OC43 and SARS-CoV-2 replication through targeting of conserved structured elements within the viral 5′-end. Nuclear magnetic resonance–based structural studies revealed specific amiloride interactions with stem loops containing bulge like structures and were predicted to be strongly bound by the lead amilorides in retrospective docking studies. Amilorides represent the first antiviral small molecules that target RNA structures within the 5′ untranslated regions and proximal region of the CoV genomes. These molecules will serve as chemical probes to further understand CoV RNA biology and can pave the way for the development of specific CoV RNA–targeted antivirals.

Metabolomics Signatures of SARS-CoV-2 Infection

For a very long time, viral infections have been considered as one of the most important causes of death and disability around the world. Through the viral infection, viruses as small pathogens enter the host cells and use hosts' biosynthesis machinery to replicate and collect infectious lineages. Moreover, they can modify hosts' metabolic pathways in order to their own purposes. Nowadays (in 2019-2020), the most famous type of viral infection which was caused by a novel type of coronavirus is called COVID-19 disease. It has claimed the lives of many people around the world and is a very serious threat to health. Since investigations of the effects of viruses on host metabolism using metabolomics tools may have given focuses on novel appropriate treatments, in the current review the authors highlighted the virus-host metabolic interactions and metabolomics perspective in COVID-19.

CoV-er all the bases: Structural perspectives of SARS-CoV-2 RNA synthesis

The ongoing Covid-19 pandemic has spurred research in the biology of the nidovirus severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2). Much focus has been on the viral RNA synthesis machinery due to its fundamental role in viral propagation. The central and essential enzyme of the RNA synthesis process, the RNA-dependent RNA polymerase (RdRp), functions in conjunction with a coterie of viral-encoded enzymes that mediate crucial nucleic acid transactions. Some of these enzymes share common features with other RNA viruses, while others play roles unique to nidoviruses or CoVs. The RdRps are proven targets for viral pathogens, and many of the other nucleic acid processing enzymes are promising targets. The purpose of this review is to summarize recent advances in our understanding of the mechanisms of RNA synthesis in CoVs. By reflecting on these studies, we hope to emphasize the remaining gaps in our knowledge. The recent onslaught of structural information related to SARS-CoV-2 RNA synthesis, in combination with previous structural, genetic and biochemical studies, have vastly improved our understanding of how CoVs replicate and process their genomic RNA. Structural biology not only provides a blueprint for understanding the function of the enzymes and cofactors in molecular detail, but also provides a basis for drug design and optimization. The concerted efforts of researchers around the world, in combination with the renewed urgency toward understanding this deadly family of viruses, may eventually yield new and improved antivirals that provide relief to the current global devastation.

In-Silico analysis reveals lower transcription efficiency of C241T variant of SARS-CoV-2 with host replication factors MADP1 and hnRNP-1

Novel severe acute respiratory syndrome-coronavirus-2 (SARS-CoV-2) has claimed more than 3.3 million lives worldwide and still counting. As per the GISAID database, the genomics of SARS-CoV-2 has been extensively studied, with more than 500 genome submissions per day. Out of several hotspot mutations within the SARS-CoV-2 genome, recent research has focused mainly on the missense variants. Moreover, significantly less attention has been accorded to delineate the role of the untranslated regions (UTRs) of the SARS-CoV-2 genome in the disease progression and etiology. One of the most frequent 5' UTR variants in the SARS-CoV-2 genome is the C241T, with a global frequency of more than 95 %. In the present study, the effect of the C241T mutation has been studied with respect to the changes in RNA structure and its interaction with the host replication factors MADP1 Zinc finger CCHC-type and RNA-binding motif 1 (hnRNP1). The results obtained from molecular docking and molecular dynamics simulation indicated weaker interaction of C241T mutant stem-loops with the host transcription factor MADP1, indicating a reduced replication efficiency. The results are also correlated with increased recovery rates and decreased death rates of global SARS-CoV-2 cases.

RNA-Protein Interaction Analysis of SARS-CoV-2 5' and 3' Untranslated Regions Reveals a Role of Lysosome-Associated Membrane Protein-2a during Viral Infection

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a positive-strand RNA virus. The viral genome is capped at the 5′ end, followed by an untranslated region (UTR). There is a poly(A) tail at the 3′ end, preceded by a UTR. The self-interaction between the RNA regulatory elements present within the 5′ and 3′ UTRs and their interaction with host/virus-encoded proteins mediate the function of the 5′ and 3′ UTRs. Using an RNA-protein interaction detection (RaPID) assay coupled to liquid chromatography with tandem mass spectrometry, we identified host interaction partners of SARS-CoV-2 5′ and 3′ UTRs and generated an RNA-protein interaction network. By combining these data with the previously known protein-protein interaction data proposed to be involved in virus replication, we generated the RNA-protein-protein interaction (RPPI) network, likely to be essential for controlling SARS-CoV-2 replication. Notably, bioinformatics analysis of the RPPI network revealed the enrichment of factors involved in translation initiation and RNA metabolism. Lysosome-associated membrane protein-2a (Lamp2a), the receptor for chaperone-mediated autophagy, is one of the host proteins that interact with the 5′ UTR. Further studies showed that the Lamp2 level is upregulated in SARS-CoV-2-infected cells and that the absence of the Lamp2a isoform enhanced the viral RNA level whereas its overexpression significantly reduced the viral RNA level. Lamp2a and viral RNA colocalize in the infected cells, and there is an increased autophagic flux in infected cells, although there is no change in the formation of autophagolysosomes. In summary, our study provides a useful resource of SARS-CoV-2 5′ and 3′ UTR binding proteins and reveals the role of Lamp2a protein during SARS-CoV-2 infection.

IMPORTANCE Replication of a positive-strand RNA virus involves an RNA-protein complex consisting of viral genomic RNA, host RNA(s), virus-encoded proteins, and host proteins. Dissecting out individual components of the replication complex will help decode the mechanism of viral replication. 5′ and 3′ UTRs in positive-strand RNA viruses play essential regulatory roles in virus replication. Here, we identified the host proteins that associate with the UTRs of SARS-CoV-2, combined those data with the previously known protein-protein interaction data (expected to be involved in virus replication), and generated the RNA-protein-protein interaction (RPPI) network. Analysis of the RPPI network revealed the enrichment of factors involved in translation initiation and RNA metabolism, which are important for virus replication. Analysis of one of the interaction partners of the 5′-UTR (Lamp2a) demonstrated its role in reducing the viral RNA level in SARS-CoV-2-infected cells. Collectively, our study provides a resource of SARS-CoV-2 UTR-binding proteins and identifies an important role for host Lamp2a protein during viral infection.

Lantern: an integrative repository of functional annotations for lncRNAs in the human genome

BACKGROUND

With advancements in omics technologies, the range of biological processes where long non-coding RNAs (lncRNAs) are involved, is expanding extensively, thereby generating the need to develop lncRNA annotation resources. Although, there are a plethora of resources for annotating genes, despite the extensive corpus of lncRNA literature, the available resources with lncRNA ontology annotations are rare.

RESULTS

We present a lncRNA annotation extractor and repository (Lantern), developed using PubMed's abstract retrieval engine and NCBO's recommender annotation system. Lantern's annotations were benchmarked against lncRNAdb's manually curated free text. Benchmarking analysis suggested that Lantern has a recall of 0.62 against lncRNAdb for 182 lncRNAs and precision of 0.8. Additionally, we also annotated lncRNAs with multiple omics annotations, including predicted cis-regulatory TFs, interactions with RBPs, tissue-specific expression profiles, protein co-expression networks, coding potential, sub-cellular localization, and SNPs for ~ 11,000 lncRNAs in the human genome, providing a one-stop dynamic visualization platform.

CONCLUSIONS

Lantern integrates a novel, accurate semi-automatic ontology annotation engine derived annotations combined with a variety of multi-omics annotations for lncRNAs, to provide a central web resource for dissecting the functional dynamics of long non-coding RNAs and to facilitate future hypothesis-driven experiments. The annotation pipeline and a web resource with current annotations for human lncRNAs are freely available on sysbio.lab.iupui.edu/lantern.

The SARS-CoV-2 subgenome landscape and its novel regulatory features

Coronavirus disease 2019 (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is currently a global pandemic. CoVs are known to generate negative subgenomes (subgenomic RNAs [sgRNAs]) through transcription-regulating sequence (TRS)-dependent template switching, but the global dynamic landscapes of coronaviral subgenomes and regulatory rules remain unclear. Here, using next-generation sequencing (NGS) short-read and Nanopore long-read poly(A) RNA sequencing in two cell types at multiple time points after infection with SARS-CoV-2, we identified hundreds of template switches and constructed the dynamic landscapes of SARS-CoV-2 subgenomes. Interestingly, template switching could occur in a bidirectional manner, with diverse SARS-CoV-2 subgenomes generated from successive template-switching events. The majority of template switches result from RNA-RNA interactions, including seed and compensatory modes, with terminal pairing status as a key determinant. Two TRS-independent template switch modes are also responsible for subgenome biogenesis. Our findings reveal the subgenome landscape of SARS-CoV-2 and its regulatory features, providing a molecular basis for understanding subgenome biogenesis and developing novel anti-viral strategies.

SARS-CoV-2 main protease suppresses type I interferon production by preventing nuclear translocation of phosphorylated IRF3

Suppression of type I interferon (IFN) response is one pathological outcome of the infection of highly pathogenic human coronaviruses. To effect this, severe acute respiratory syndrome coronavirus (SARS-CoV) and SARS-CoV-2 encode multiple IFN antagonists. In this study, we reported on the IFN antagonism of SARS-CoV-2 main protease NSP5. NSP5 proteins of both SARS-CoV and SARS-CoV-2 counteracted Sendai virus-induced IFN production. NSP5 variants G15S and K90R commonly seen in circulating strains of SARS-CoV-2 retained the IFN-antagonizing property. The suppressive effect of NSP5 on IFN-β gene transcription induced by RIG-I, MAVS, TBK1 and IKKϵ suggested that NSP5 likely acts at a step downstream of IRF3 phosphorylation in the cytoplasm. NSP5 did not influence steady-state expression or phosphorylation of IRF3, suggesting that IRF3, regardless of its phosphorylation state, might not be the substrate of NSP5 protease. However, nuclear translocation of phosphorylated IRF3 was severely compromised in NSP5-expressing cells. Taken together, our work revealed a new mechanism by which NSP5 proteins encoded by SARS-CoV and SARS-CoV-2 antagonize IFN production by retaining phosphorylated IRF3 in the cytoplasm. Our findings have implications in rational design and development of antiviral agents against SARS-CoV-2.

A novel end-to-end method to predict RNA secondary structure profile based on bidirectional LSTM and residual neural network

Background

Studies have shown that RNA secondary structure, a planar structure formed by paired bases, plays diverse vital roles in fundamental life activities and complex diseases. RNA secondary structure profile can record whether each base is paired with others. Hence, accurate prediction of secondary structure profile can help to deduce the secondary structure and binding site of RNA. RNA secondary structure profile can be obtained through biological experiment and calculation methods. Of them, the biological experiment method involves two ways: chemical reagent and biological crystallization. The chemical reagent method can obtain a large number of prediction data, but its cost is high and always associated with high noise, making it difficult to get results of all bases on RNA due to the limited of sequencing coverage. By contrast, the biological crystallization method can lead to accurate results, yet heavy experimental work and high costs are required. On the other hand, the calculation method is CROSS, which comprises a three-layer fully connected neural network. However, CROSS can not completely learn the features of RNA secondary structure profile since its poor network structure, leading to its low performance.

Results

In this paper, a novel end-to-end method, named as “RPRes, was proposed to predict RNA secondary structure profile based on Bidirectional LSTM and Residual Neural Network.

Conclusions

RPRes utilizes data sets generated by multiple biological experiment methods as the training, validation, and test sets to predict profile, which can compatible with numerous prediction requirements. Compared with the biological experiment method, RPRes has reduced the costs and improved the prediction efficiency. Compared with the state-of-the-art calculation method CROSS, RPRes has significantly improved performance.

Computational strategies to combat COVID-19: useful tools to accelerate SARS-CoV-2 and coronavirus research

SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) is a novel virus of the family Coronaviridae. The virus causes the infectious disease COVID-19. The biology of coronaviruses has been studied for many years. However, bioinformatics tools designed explicitly for SARS-CoV-2 have only recently been developed as a rapid reaction to the need for fast detection, understanding and treatment of COVID-19. To control the ongoing COVID-19 pandemic, it is of utmost importance to get insight into the evolution and pathogenesis of the virus. In this review, we cover bioinformatics workflows and tools for the routine detection of SARS-CoV-2 infection, the reliable analysis of sequencing data, the tracking of the COVID-19 pandemic and evaluation of containment measures, the study of coronavirus evolution, the discovery of potential drug targets and development of therapeutic strategies. For each tool, we briefly describe its use case and how it advances research specifically for SARS-CoV-2. All tools are free to use and available online, either through web applications or public code repositories. Contact:evbc@unj-jena.de.

Novel perspectives for SARS-CoV-2 genome browsing

SARS-CoV-2 has spread worldwide and caused social, economic, and health turmoil. The first genome assembly of SARS-CoV-2 was produced in Wuhan, and it is widely used as a reference. Subsequently, more than a hundred additional SARS-CoV-2 genomes have been sequenced. While the genomes appear to be mostly identical, there are variations. Therefore, an alignment of all available genomes and the derived consensus sequence could be used as a reference, better serving the science community. Variations are significant, but representing them in a genome browser can become, especially if their sequences are largely identical. Here we summarize the variation in one track. Other information not currently found in genome browsers for SARS-CoV-2, such as predicted miRNAs and predicted TRS as well as secondary structure information, were also added as tracks to the consensus genome. We believe that a genome browser based on the consensus sequence is better suited when considering worldwide effects and can become a valuable resource in the combating of COVID-19. The genome browser is available at http://cov.iaba.online.

The global and local distribution of RNA structure throughout the SARS-CoV-2 genome

SARS-CoV-2 is the causative viral agent of COVID-19, the disease at the center of the current global pandemic. While knowledge of highly structured regions is integral for mechanistic insights into the viral infection cycle, very little is known about the location and folding stability of functional elements within the massive, ∼30kb SARS-CoV-2 RNA genome. In this study, we analyze the folding stability of this RNA genome relative to the structural landscape of other well-known viral RNAs. We present an in-silico pipeline to predict regions of high base pair content across long genomes and to pinpoint hotspots of well-defined RNA structures, a method that allows for direct comparisons of RNA structural complexity within the several domains in SARS-CoV-2 genome. We report that the SARS-CoV-2 genomic propensity for stable RNA folding is exceptional among RNA viruses, superseding even that of HCV, one of the most structured viral RNAs in nature. Furthermore, our analysis suggests varying levels of RNA structure across genomic functional regions, with accessory and structural ORFs containing the highest structural density in the viral genome. Finally, we take a step further to examine how individual RNA structures formed by these ORFs are affected by the differences in genomic and subgenomic contexts, which given the technical difficulty of experimentally separating cellular mixtures of sgRNA from gRNA, is a unique advantage of our in-silico pipeline. The resulting findings provide a useful roadmap for planning focused empirical studies of SARS-CoV-2 RNA biology, and a preliminary guide for exploring potential SARS-CoV-2 RNA drug targets.Importance The RNA genome of SARS-CoV-2 is among the largest and most complex viral genomes, and yet its RNA structural features remain relatively unexplored. Since RNA elements guide function in most RNA viruses, and they represent potential drug targets, it is essential to chart the architectural features of SARS-CoV-2 and pinpoint regions that merit focused study. Here we show that RNA folding stability of SARS-CoV-2 genome is exceptional among viral genomes and we develop a method to directly compare levels of predicted secondary structure across SARS-CoV-2 domains. Remarkably, we find that coding regions display the highest structural propensity in the genome, forming motifs that differ between the genomic and subgenomic contexts. Our approach provides an attractive strategy to rapidly screen for candidate structured regions based on base pairing potential and provides a readily interpretable roadmap to guide functional studies of RNA viruses and other pharmacologically relevant RNA transcripts.

Inhibition of anti-viral stress granule formation by coronavirus endoribonuclease nsp15 ensures efficient virus replication

Cytoplasmic stress granules (SGs) are generally triggered by stress-induced translation arrest for storing mRNAs. Recently, it has been shown that SGs exert anti-viral functions due to their involvement in protein synthesis shut off and recruitment of innate immune signaling intermediates. The largest RNA viruses, coronaviruses, impose great threat to public safety and animal health; however, the significance of SGs in coronavirus infection is largely unknown. Infectious Bronchitis Virus (IBV) is the first identified coronavirus in 1930s and has been prevalent in poultry farm for many years. In this study, we provided evidence that IBV overcomes the host antiviral response by inhibiting SGs formation via the virus-encoded endoribonuclease nsp15. By immunofluorescence analysis, we observed that IBV infection not only did not trigger SGs formation in approximately 80% of the infected cells, but also impaired the formation of SGs triggered by heat shock, sodium arsenite, or NaCl stimuli. We further demonstrated that the intrinsic endoribonuclease activity of nsp15 was responsible for the interference of SGs formation. In fact, nsp15-defective recombinant IBV (rIBV-nsp15-H238A) greatly induced the formation of SGs, along with accumulation of dsRNA and activation of PKR, whereas wild type IBV failed to do so. Consequently, infection with rIBV-nsp15-H238A strongly triggered transcription of IFN-β which in turn greatly affected rIBV-nsp15-H238A replication. Further analysis showed that SGs function as an antiviral hub, as demonstrated by the attenuated IRF3-IFN response and increased production of IBV in SG-defective cells. Additional evidence includes the aggregation of pattern recognition receptors (PRRs) and signaling intermediates to the IBV-induced SGs. Collectively, our data demonstrate that the endoribonuclease nsp15 of IBV interferes with the formation of antiviral hub SGs by regulating the accumulation of viral dsRNA and by antagonizing the activation of PKR, eventually ensuring productive virus replication. We further demonstrated that nsp15s from PEDV, TGEV, SARS-CoV, and SARS-CoV-2 harbor the conserved function to interfere with the formation of chemically-induced SGs. Thus, we speculate that coronaviruses employ similar nsp15-mediated mechanisms to antagonize the host anti-viral SGs formation to ensure efficient virus replication.

Coronavirus encodes the conserved endoribonuclease nsp15, which has been reported to antagonize IFN responses by mediating evasion of recognition by dsRNA sensors. SGs are part of the host cell anti-viral response; not surprisingly, viruses in turn produce an array of antagonists to counteract such host response. Here, we show that IBV prevents the formation of SGs via nsp15, by reducing the accumulation of viral dsRNA, thereby evading the activation of PKR, phosphorylation of eIF2α, and formation of SGs. Depletion of SG scaffold proteins G3BP1/2 decreases IRF3-IFN response and increases the production of IBV. When overexpressed alone, nsp15s from different coronaviruses (IBV, PEDV, TGEV, SARS-CoV, and SARS-CoV-2) interferes with chemically- and physically-induced SGs, probably by targeting essential SGs assembly factors. In this way, coronaviruses antagonize the formation of SGs by nsp15, via reducing the viral dsRNA accumulation and sequestering/depleting critical component of SGs. To our knowledge, this is the first report describing the role of coronavirus nsp15 in the suppression of integral stress response, in crosstalk with anti-innate immune response.

Rfam 14: expanded coverage of metagenomic, viral and microRNA families

Rfam is a database of RNA families where each of the 3444 families is represented by a multiple sequence alignment of known RNA sequences and a covariance model that can be used to search for additional members of the family. Recent developments have involved expert collaborations to improve the quality and coverage of Rfam data, focusing on microRNAs, viral and bacterial RNAs. We have completed the first phase of synchronising microRNA families in Rfam and miRBase, creating 356 new Rfam families and updating 40. We established a procedure for comprehensive annotation of viral RNA families starting with Flavivirus and Coronaviridae RNAs. We have also increased the coverage of bacterial and metagenome-based RNA families from the ZWD database. These developments have enabled a significant growth of the database, with the addition of 759 new families in Rfam 14. To facilitate further community contribution to Rfam, expert users are now able to build and submit new families using the newly developed Rfam Cloud family curation system. New Rfam website features include a new sequence similarity search powered by RNAcentral, as well as search and visualisation of families with pseudoknots. Rfam is freely available at https://rfam.org.

Structural Insight Into the SARS-CoV-2 Nucleocapsid Protein C-Terminal Domain Reveals a Novel Recognition Mechanism for Viral Transcriptional Regulatory Sequences

Coronavirus disease 2019 (COVID-19) has caused massive disruptions to society and the economy, and the transcriptional regulatory mechanisms behind the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) are poorly understood. Herein, we determined the crystal structure of the SARS-CoV-2 nucleocapsid protein C-terminal domain (CTD) at a resolution of 2.0 Å, and demonstrated that the CTD has a comparable distinct electrostatic potential surface to equivalent domains of other reported CoVs, suggesting that the CTD has novel roles in viral RNA binding and transcriptional regulation. Further in vitro biochemical assays demonstrated that the viral genomic intergenic transcriptional regulatory sequences (TRSs) interact with the SARS-CoV-2 nucleocapsid protein CTD with a flanking region. The unpaired adeno dinucleotide in the TRS stem-loop structure is a major determining factor for their interactions. Taken together, these results suggested that the nucleocapsid protein CTD is responsible for the discontinuous viral transcription mechanism by recognizing the different patterns of viral TRS during transcription.

Translational control of coronaviruses

Coronaviruses represent a large family of enveloped RNA viruses that infect a large spectrum of animals. In humans, the severe acute respiratory syndrome coronavirus type 2 (SARS-CoV-2) is responsible for the current COVID-19 pandemic and is genetically related to SARS-CoV and Middle East respiratory syndrome-related coronavirus (MERS-CoV), which caused outbreaks in 2002 and 2012, respectively. All viruses described to date entirely rely on the protein synthesis machinery of the host cells to produce proteins required for their replication and spread. As such, virus often need to control the cellular translational apparatus to avoid the first line of the cellular defense intended to limit the viral propagation. Thus, coronaviruses have developed remarkable strategies to hijack the host translational machinery in order to favor viral protein production. In this review, we will describe some of these strategies and will highlight the role of viral proteins and RNAs in this process.

Amilorides inhibit SARS-CoV-2 replication in vitro by targeting RNA structures

The SARS-CoV-2 pandemic, and the likelihood of future coronavirus pandemics, has rendered our understanding of coronavirus biology more essential than ever. Small molecule chemical probes offer to both reveal novel aspects of virus replication and to serve as leads for antiviral therapeutic development. The RNA-biased amiloride scaffold was recently tuned to target a viral RNA structure critical for translation in enterovirus 71, ultimately uncovering a novel mechanism to modulate positive-sense RNA viral translation and replication. Analysis of CoV RNA genomes reveal many conserved RNA structures in the 5'-UTR and proximal region critical for viral translation and replication, including several containing bulge-like secondary structures suitable for small molecule targeting. Following phylogenetic conservation analysis of this region, we screened an amiloride-based small molecule library against a less virulent human coronavirus, OC43, to identify lead ligands. Amilorides inhibited OC43 replication as seen in viral plaque assays. Select amilorides also potently inhibited replication competent SARS-CoV-2 as evident in the decreased levels of cell free virions in cell culture supernatants of treated cells. Reporter screens confirmed the importance of RNA structures in the 5'-end of the viral genome for small molecule activity. Finally, NMR chemical shift perturbation studies of the first six stem loops of the 5'-end revealed specific amiloride interactions with stem loops 4, 5a, and 6, all of which contain bulge like structures and were predicted to be strongly bound by the lead amilorides in retrospective docking studies. Taken together, the use of multiple orthogonal approaches allowed us to identify the first small molecules aimed at targeting RNA structures within the 5'-UTR and proximal region of the CoV genome. These molecules will serve as chemical probes to further understand CoV RNA biology and can pave the way for the development of specific CoV RNA-targeted antivirals.

SARS-CoV-2 Disrupts Splicing, Translation, and Protein Trafficking to Suppress Host Defenses

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a recently identified coronavirus that causes the respiratory disease known as coronavirus disease 2019 (COVID-19). Despite the urgent need, we still do not fully understand the molecular basis of SARS-CoV-2 pathogenesis. Here, we comprehensively define the interactions between SARS-CoV-2 proteins and human RNAs. NSP16 binds to the mRNA recognition domains of the U1 and U2 splicing RNAs and acts to suppress global mRNA splicing upon SARS-CoV-2 infection. NSP1 binds to 18S ribosomal RNA in the mRNA entry channel of the ribosome and leads to global inhibition of mRNA translation upon infection. Finally, NSP8 and NSP9 bind to the 7SL RNA in the signal recognition particle and interfere with protein trafficking to the cell membrane upon infection. Disruption of each of these essential cellular functions acts to suppress the interferon response to viral infection. Our results uncover a multipronged strategy utilized by SARS-CoV-2 to antagonize essential cellular processes to suppress host defenses.