High-resolution structures of the SARS-CoV-2 2'-O-methyltransferase reveal strategies for structure-based inhibitor design

Natural evidence of coronaviral 2'-O-methyltransferase activity affecting viral pathogenesis via improved substrate RNA binding

Previous studies through targeted mutagenesis of K-D-K-E motif have demonstrated that 2'-O-MTase activity is essential for efficient viral replication and immune evasion. However, the K-D-K-E catalytic motif of 2'-O-MTase is highly conserved across numerous viruses, including flaviviruses, vaccinia viruses, coronaviruses, and extends even to mammals. Here, we observed a stronger 2'-O-MTase activity in SARS-CoV-2 compared to SARS-CoV, despite the presence of a consistently active catalytic center. We further identified critical residues (Leu-36, Asn-138 and Ile-153) which served as determinants of discrepancy in 2'-O-MTase activity between SARS-CoV-2 and SARS-CoV. These residues significantly enhanced the RNA binding affinity of 2'-O-MTase and boosted its versatility toward RNA substrates. Of interest, a triple substitution (Leu36 → Ile36, Asn138 → His138, Ile153 → Leu153, from SARS-CoV-2 to SARS-CoV) within nsp16 resulted in a proportional reduction in viral 2'-O-methylation and impaired viral replication. Furthermore, it led to a significant upregulation of type I interferon (IFN-I) and proinflammatory cytokines both in vitro and vivo, relying on the cooperative sensing of melanoma differentiation-associated protein 5 (MDA5) and laboratory of genetics and physiology 2 (LGP2). In conclusion, our findings demonstrated that alterations in residues other than K-D-K-E of 2'-O-MTase may affect viral replication and subsequently influence pathogenesis. Monitoring changes in nsp16 residues is crucial as it may aid in identifying and assessing future alteration in viral pathogenicity resulting from natural mutations occurring in nsp16.

Structure-Based Virtual Screening for Methyltransferase Inhibitors of SARS-CoV-2 nsp14 and nsp16

The ongoing COVID-19 pandemic still threatens human health around the world. The methyltransferases (MTases) of SARS-CoV-2, specifically nsp14 and nsp16, play crucial roles in the methylation of the N7 and 2'-O positions of viral RNA, making them promising targets for the development of antiviral drugs. In this work, we performed structure-based virtual screening for nsp14 and nsp16 using the screening workflow (HTVS, SP, XP) of Schrödinger 2019 software, and we carried out biochemical assays and molecular dynamics simulation for the identification of potential MTase inhibitors. For nsp14, we screened 239,000 molecules, leading to the identification of three hits A1-A3 showing N7-MTase inhibition rates greater than 60% under a concentration of 50 µM. For the SAM binding and nsp10-16 interface sites of nsp16, the screening of 210,000 and 237,000 molecules, respectively, from ZINC15 led to the discovery of three hit compounds B1-B3 exhibiting more than 45% of 2'-O-MTase inhibition under 50 µM. These six compounds with moderate MTase inhibitory activities could be used as novel candidates for the further development of anti-SARS-CoV-2 drugs.

Therapeutic implications of quercetin and its derived-products in COVID-19 protection and prophylactic

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), a novel human coronavirus, which has triggered a global pandemic of the coronavirus infectious disease 2019 (COVID-19). Outbreaks of emerging infectious diseases continue to challenge human health worldwide. The virus conquers human cells through the angiotensin-converting enzyme 2 receptor-driven pathway by mostly targeting the human respiratory tract. Quercetin is a natural flavonoid widely represented in the plant kingdom. Cumulative evidence has demonstrated that quercetin and its derivatives have various pharmacological properties including anti-cancer, anti-hypertension, anti-hyperlipidemia, anti-hyperglycemia, anti-microbial, antiviral, neuroprotective, and cardio-protective effects, because it is a potential treatment for severe inflammation and acute respiratory distress syndrome. Furthermore, it is the main life-threatening condition in patients with COVID-19. This article provides a comprehensive review of the primary literature on the predictable effectiveness of quercetin and its derivatives docked to multi-target of SARS-CoV-2 and host cells via in silico and some of validation through in vitro, in vivo, and clinically to fight SARS-CoV-2 infections, contribute to the reduction of inflammation, which suggests the preventive and therapeutic latency of quercetin and its derived-products against COVID-19 pandemic, multisystem inflammatory syndromes (MIS), and long-COVID.

Cheminformatics approach to identify andrographolide derivatives as dual inhibitors of methyltransferases (nsp14 and nsp16) of SARS-CoV-2

The Covid-19 pandemic outbreak has accelerated tremendous efforts to discover a therapeutic strategy that targets severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) to control viral infection. Various viral proteins have been identified as potential drug targets, however, to date, no specific therapeutic cure is available against the SARS-CoV-2. To address this issue, the present work reports a systematic cheminformatic approach to identify the potent andrographolide derivatives that can target methyltransferases of SARS-CoV-2, i.e. nsp14 and nsp16 which are crucial for the replication of the virus and host immune evasion. A consensus of cheminformatics methodologies including virtual screening, molecular docking, ADMET profiling, molecular dynamics simulations, free-energy landscape analysis, molecular mechanics generalized born surface area (MM-GBSA), and density functional theory (DFT) was utilized. Our study reveals two new andrographolide derivatives (PubChem CID: 2734589 and 138968421) as natural bioactive molecules that can form stable complexes with both proteins via hydrophobic interactions, hydrogen bonds and electrostatic interactions. The toxicity analysis predicts class four toxicity for both compounds with LD50 value in the range of 500-700 mg/kg. MD simulation reveals the stable formation of the complex for both the compounds and their average trajectory values were found to be lower than the control inhibitor and protein alone. MMGBSA analysis corroborates the MD simulation result and showed the lowest energy for the compounds 2734589 and 138968421. The DFT and MEP analysis also predicts the better reactivity and stability of both the hit compounds. Overall, both andrographolide derivatives exhibit good potential as potent inhibitors for both nsp14 and nsp16 proteins, however, in-vitro and in vivo assessment would be required to prove their efficacy and safety in clinical settings. Moreover, the drug discovery strategy aiming at the dual target approach might serve as a useful model for inventing novel drug molecules for various other diseases.

AI-guided pipeline for protein-protein interaction drug discovery identifies a SARS-CoV-2 inhibitor

Protein-protein interactions (PPIs) offer great opportunities to expand the druggable proteome and therapeutically tackle various diseases, but remain challenging targets for drug discovery. Here, we provide a comprehensive pipeline that combines experimental and computational tools to identify and validate PPI targets and perform early-stage drug discovery. We have developed a machine learning approach that prioritizes interactions by analyzing quantitative data from binary PPI assays or AlphaFold-Multimer predictions. Using the quantitative assay LuTHy together with our machine learning algorithm, we identified high-confidence interactions among SARS-CoV-2 proteins for which we predicted three-dimensional structures using AlphaFold-Multimer. We employed VirtualFlow to target the contact interface of the NSP10-NSP16 SARS-CoV-2 methyltransferase complex by ultra-large virtual drug screening. Thereby, we identified a compound that binds to NSP10 and inhibits its interaction with NSP16, while also disrupting the methyltransferase activity of the complex, and SARS-CoV-2 replication. Overall, this pipeline will help to prioritize PPI targets to accelerate the discovery of early-stage drug candidates targeting protein complexes and pathways.

Informatics and Computational Approaches for the Discovery and Optimization of Natural Product-Inspired Inhibitors of the SARS-CoV-2 2'-O-Methyltransferase

The urgent need for new classes of orally available, safe, and effective antivirals─covering a breadth of emerging viruses─is evidenced by the loss of life and economic challenges created by the HIV-1 and SARS-CoV-2 pandemics. As frontline interventions, small-molecule antivirals can be deployed prophylactically or postinfection to control the initial spread of outbreaks by reducing transmissibility and symptom severity. Natural products have an impressive track record of success as prototypic antivirals and continue to provide new drugs through synthesis, medicinal chemistry, and optimization decades after discovery. Here, we demonstrate an approach using computational analysis typically used for rational drug design to identify and develop natural product-inspired antivirals. This was done with the goal of identifying natural product prototypes to aid the effort of progressing toward safe, effective, and affordable broad-spectrum inhibitors of Betacoronavirus replication by targeting the highly conserved RNA 2'-O-methyltransferase (2'-O-MTase). Machaeriols RS-1 (7) and RS-2 (8) were identified using a previously outlined informatics approach to first screen for natural product prototypes, followed by in silico-guided synthesis. Both molecules are based on a rare natural product group. The machaeriols (3-6), isolated from the genus Machaerium, endemic to Amazonia, inhibited the SARS-CoV-2 2'-O-MTase more potently than the positive control, Sinefungin (2), and in silico modeling suggests distinct molecular interactions. This report highlights the potential of computationally driven screening to leverage natural product libraries and improve the efficiency of isolation or synthetic analog development.

5'-cap RNA/SAM mimetic conjugates as bisubstrate inhibitors of viral RNA cap 2'-O-methyltransferases

Viral RNA cap 2'-O-methyltransferases are considered promising therapeutic targets for antiviral treatments, as they play a key role in the formation of viral RNA cap-1 structures to escape the host immune system. A better understanding of how they interact with their natural substrates (RNA and the methyl donor SAM) would enable the rational development of potent inhibitors. However, as few structures of 2'-O-MTases in complex with RNA have been described, little is known about substrate recognition by these MTases. For this, chemical tools mimicking the state in which the cap RNA substrate and SAM cofactor are bound in the enzyme's catalytic pocket may prove useful. In this work, we designed and synthesized over 30 RNA conjugates that contain a short oligoribonucleotide (ORN with 4 or 6 nucleotides) with the first nucleotide 2'-O-attached to an adenosine by linkers of different lengths and containing S or N-heteroatoms, or a 1,2,3-triazole ring. These ORN conjugates bearing or not a cap structure at 5'-extremity mimic the methylation transition state with RNA substrate/SAM complex as bisubstrates of 2'-O-MTases. The ORN conjugates were synthesized either by the incorporation of a dinucleoside phosphoramidite during RNA elongation or by click chemistry performed on solid-phase post-RNA elongation. Their ability to inhibit the activity of the nsp16/nsp10 complex of SARS-CoV-2 and the NS5 protein of dengue and Zika viruses was assessed. Significant submicromolar IC50 values and Kd values in the µM range were found, suggesting a possible interaction of some ORN conjugates with these viral 2'-O-MTases.

Integrated molecular and quantum mechanical approach to identify novel potent natural bioactive compound against 2'-O-methyltransferase (nsp16) of SARS-CoV-2

With the advent of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) outbreak, efforts are still in progress to find out a functional cure for the infection. Among the various protein targets, nsp16 capping protein is one of the vital targets for drug development as it protects the virus against the host cell nucleases and evading innate immunity. The nsp16 protein forms a heterodimer with a co-factor nsp10 and triggers 2'-O-methyltransferase activity which catalyzes the conversion of S-adenosyl methionine into S-adenosyl homocysteine. The free methyl group is transferred to the 2'-O position on ribose sugar at the 5' end of mRNA to form the cap-1 structure which is essential for replication of the virus and evading the innate immunity of the host. In this study, we identify a potential lead natural bioactive compound against nsp16 protein by systematic cheminformatic analysis of more than 144k natural compounds. Virtual screening, molecular docking interactions, ADMET profiling, molecular dynamics (MD) simulations, molecular mechanics-generalized born surface area (MM-GBSA), free energy analysis and density functional theory analysis were used to discover the potential lead compound. Our investigation revealed that ZINC8952607 (methyl-[(6-methyl-2,3,4,9-tetrahydro-1H-carbazol-1-yl)aminomethyl]BLAHone) has the greatest binding affinity and best pharmacokinetic parameters due to presence of carbazol and BLAHone (biaryl moiety). Further, time-dependent MD simulation analysis substantiates the stability and rigidness of nsp16 protein even after interaction with the lead compound. We believe that the compound ZINC8952607 might establish as a novel natural drug candidate against CoVID-19 infection.Communicated by Ramaswamy H. Sarma.

Targeting cap1 RNA methyltransferases as an antiviral strategy

Methylation is one of the critical modifications that regulates numerous biological processes. Guanine capping and methylation at the 7th position (m7G) have been shown to mature mRNA for increased RNA stability and translational efficiency. The m7G capped cap0 RNA remains immature and requires additional methylation at the first nucleotide (N1-2'-O-Me), designated as cap1, to achieve full maturation. This cap1 RNA with N1-2'-O-Me prevents its recognition by innate immune sensors as non-self. Viruses have also evolved various strategies to produce self-like capped RNAs with the N1-2'-O-Me that potentially evades the antiviral response and establishes an efficient replication. In this review, we focus on the importance of the presence of N1-2'-O-Me in viral RNAs and discuss the potential for drug development by targeting host and viral N1-2'-O-methyltransferases.

Applying the bioisosterism strategy to obtain lead compounds against SARS-CoV-2 cysteine proteases: An in-silico approach

SARS-CoV-2 cysteine proteases are essential nonstructural proteins due to their role in the formation of the virus multiple enzyme replication-transcription complex. As a result, those functional proteins are extremely relevant targets in the development of a new drug candidate to fight COVID-19. Based on this fact and guided by the bioisosterism strategy, the present work has selected 126 out of 1050 ligands from DrugBank website. Subsequently, 831 chemical analogs containing bioisosteres, some of which became structurally simplified, were created using the MB-Isoster software, and molecular docking simulations were performed using AutoDock Vina. Finally, a study of physicochemical properties, along with pharmacokinetic profiles, was carried out through SwissADME and ADMETlab 2.0 platforms. The promising results obtained with the molecules encoded as DB00549_BI_005, DB04868_BI_003, DB11984_BI_002, DB12364_BI_006 and DB12805_BI_004 must be confirmed by molecular dynamics studies, followed by in vitro and in vivo empirical tests that ratify the advocated in-silico results.

Emerging variants of SARS-CoV-2 NSP10 highlight strong functional conservation of its binding to two non-structural proteins, NSP14 and NSP16

The coronavirus SARS-CoV-2 protects its RNA from being recognized by host immune responses by methylation of its 5' end, also known as capping. This process is carried out by two enzymes, non-structural protein 16 (NSP16) containing 2'-O-methyltransferase and NSP14 through its N7 methyltransferase activity, which are essential for the replication of the viral genome as well as evading the host's innate immunity. NSP10 acts as a crucial cofactor and stimulator of NSP14 and NSP16. To further understand the role of NSP10, we carried out a comprehensive analysis of >13 million globally collected whole-genome sequences (WGS) of SARS-CoV-2 obtained from the Global Initiative Sharing All Influenza Data (GISAID) and compared it with the reference genome Wuhan/WIV04/2019 to identify all currently known variants in NSP10. T12I, T102I, and A104V in NSP10 have been identified as the three most frequent variants and characterized using X-ray crystallography, biophysical assays, and enhanced sampling simulations. In contrast to other proteins such as spike and NSP6, NSP10 is significantly less prone to mutation due to its crucial role in replication. The functional effects of the variants were examined for their impact on the binding affinity and stability of both NSP14-NSP10 and NSP16-NSP10 complexes. These results highlight the limited changes induced by variant evolution in NSP10 and reflect on the critical roles NSP10 plays during the SARS-CoV-2 life cycle. These results also indicate that there is limited capacity for the virus to overcome inhibitors targeting NSP10 via the generation of variants in inhibitor binding pockets.

The role of post-transcriptional regulation in SARS-CoV-2 infection and pathogenicity

The COVID-19 pandemic caused by the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) virus has had a significant impact on global social and economic stability. To combat this, researchers have turned to omics approaches, particularly epitranscriptomics, to limit infection and develop effective therapeutic strategies. Multi-omics can provide the host response dynamics during multiple disease phases to reveal the molecular and cellular landscapes. Epitranscriptomics focuses on the mechanisms of gene transcription in cells and tissues and the relationship between genetic material and epigenetic regulation. This review highlights the role of post-transcriptional regulation in SARS-CoV-2, which affect various processes such as virus infection, replication, immunogenicity, and pathogenicity. The review also explains the formation mechanism of post-transcriptional modifications and how they can be regulated to combat viral infection and pathogenicity.

Coronavirus 2'-O-methyltransferase: A promising therapeutic target

Coronaviruses (CoVs) have been the source of multiple epidemics and a global pandemic since the start of century, and there is an urgent need to understand CoV biology and develop better therapeutics. Here, we review the role of NSP16 in CoV replication, specifically its importance to 2'-O-methylation and CoV RNA capping. We describe the attenuation phenotypes of NSP16-mutant CoVs, the roles of MDA5 and IFITs in sensing and antagonizing viral RNA lacking 2'O methylation, and the dependence on 2'-O-methylation in other virus families. We also detail the growing body of research into targeting 2'-O-methylation for therapeutics or as a platform for live attenuated vaccines. Beyond its role in RNA capping, NSP16 may have yet uncharacterized importance to CoV replication, highlighting the need for continued studies into NSP16 functions. Understanding the full contribution of NSP16 to the replicative fitness of CoVs will better inform the development of treatments against future CoV outbreaks.

Discovery of a Druggable, Cryptic Pocket in SARS-CoV-2 nsp16 Using Allosteric Inhibitors

A collaborative, open-science team undertook discovery of novel small molecule inhibitors of the SARS-CoV-2 nsp16-nsp10 2'-O-methyltransferase using a high throughput screening approach with the potential to reveal new inhibition strategies. This screen yielded compound 5a, a ligand possessing an electron-deficient double bond, as an inhibitor of SARS-CoV-2 nsp16 activity. Surprisingly, X-ray crystal structures revealed that 5a covalently binds within a previously unrecognized cryptic pocket near the S-adenosylmethionine binding cleft in a manner that prevents occupation by S-adenosylmethionine. Using a multidisciplinary approach, we examined the mechanism of binding of compound 5a to the nsp16 cryptic pocket and developed 5a derivatives that inhibited nsp16 activity and murine hepatitis virus replication in rat lung epithelial cells but proved cytotoxic to cell lines canonically used to examine SARS-CoV-2 infection. Our study reveals the druggability of this newly discovered SARS-CoV-2 nsp16 cryptic pocket, provides novel tool compounds to explore the site, and suggests a new approach for discovery of nsp16 inhibition-based pan-coronavirus therapeutics through structure-guided drug design.

Potent multi-target natural inhibitors against SARS-CoV-2 from medicinal plants of the Himalaya: a discovery from hybrid machine learning, chemoinformatics, and simulation assisted screening

The emergence and immune evasion ability of SARS-CoV-2 Omicron strains, mainly BA.5.2 and BF.7 and other variants of concern have raised global apprehensions. With this context, the discovery of multitarget inhibitors may be proven more comprehensive paradigm than its one-drug-to-one target counterpart. In the current study, a library of 271 phytochemicals from 25 medicinal plants from the Indian Himalayan Region has been virtually screened against SARS-CoV-2 by targeting nine virus proteins, viz., papain-like protease, main protease, nsp12, helicase, nsp14, nsp15, nsp16, envelope, and nucleocapsid for screening of a multi-target inhibitor against the viral replication. Initially, 94 phytochemicals were screened by a hybrid machine learning model constructed by combining 6 confirmatory bioassays against SARS-CoV-2 replication using an instance-based learner lazy k-nearest neighbour classifier. Further, 25 screened compounds with excellent drug-like properties were subjected to molecular docking. The phytochemical Cepharadione A from the plant Piper longum showed binding potential against four proteins with the highest binding energy of -10.90 kcal/mol. The compound has acceptable absorption, distribution, metabolism, excretion, and toxicity properties and exhibits stable binding behaviour in terms of root mean square deviation (0.068 ± 0.05 nm), root-mean-square fluctuation, hydrogen bonds, solvent accessible surface area (83.88-161.89 nm2), and molecular mechanics Poisson-Boltzmann surface area during molecular dynamics simulation of 200 ns with selected target proteins. Concerning the utility of natural compounds in the therapeutics formulation, Cepharadione A could be further investigated as a remarkable lead candidate for the development of therapeutic drugs against SARS-CoV-2.Communicated by Ramaswamy H. Sarma.

What do we know about the function of SARS-CoV-2 proteins?

The COVID-19 pandemic has highlighted the importance in the understanding of the biology of SARS-CoV-2. After more than two years since the first report of COVID-19, it remains crucial to continue studying how SARS-CoV-2 proteins interact with the host metabolism to cause COVID-19. In this review, we summarize the findings regarding the functions of the 16 non-structural, 6 accessory and 4 structural SARS-CoV-2 proteins. We place less emphasis on the spike protein, which has been the subject of several recent reviews. Furthermore, comprehensive reviews about COVID-19 therapeutic have been also published. Therefore, we do not delve into details on these topics; instead we direct the readers to those other reviews. To avoid confusions with what we know about proteins from other coronaviruses, we exclusively report findings that have been experimentally confirmed in SARS-CoV-2. We have identified host mechanisms that appear to be the primary targets of SARS-CoV-2 proteins, including gene expression and immune response pathways such as ribosome translation, JAK/STAT, RIG-1/MDA5 and NF-kβ pathways. Additionally, we emphasize the multiple functions exhibited by SARS-CoV-2 proteins, along with the limited information available for some of these proteins. Our aim with this review is to assist researchers and contribute to the ongoing comprehension of SARS-CoV-2's pathogenesis.

Identification of motif-based interactions between SARS-CoV-2 protein domains and human peptide ligands pinpoint antiviral targets

The virus life cycle depends on host-virus protein-protein interactions, which often involve a disordered protein region binding to a folded protein domain. Here, we used proteomic peptide phage display (ProP-PD) to identify peptides from the intrinsically disordered regions of the human proteome that bind to folded protein domains encoded by the SARS-CoV-2 genome. Eleven folded domains of SARS-CoV-2 proteins were found to bind 281 peptides from human proteins, and affinities of 31 interactions involving eight SARS-CoV-2 protein domains were determined (KD ∼ 7-300 μM). Key specificity residues of the peptides were established for six of the interactions. Two of the peptides, binding Nsp9 and Nsp16, respectively, inhibited viral replication. Our findings demonstrate how high-throughput peptide binding screens simultaneously identify potential host-virus interactions and peptides with antiviral properties. Furthermore, the high number of low-affinity interactions suggest that overexpression of viral proteins during infection may perturb multiple cellular pathways.

Identifying inhibitors of NSP16-NSP10 of SARS-CoV-2 from large databases

The COVID-19 pandemic, which has already claimed millions of lives, continues to pose a serious threat to human health, requiring the development of new effective drugs. Non-structural proteins of SARS-CoV-2 play an important role in viral replication and infection. Among them, NSP16 (non-structured protein 16) and its cofactor NSP10 (non-structured protein 10) perform C2'-O methylation at the 5' end of the viral RNA, which promotes efficient virus replication. Therefore, the NSP16-NSP10 complex becomes an attractive target for drug development. Using a multi-step virtual screening protocol which includes Lipinski's rule, docking, steered molecular dynamics and umbrella sampling, we searched for potential inhibitors from the PubChem and anti-HIV databases. It has been shown that CID 135566620 compound from PubChem is the best candidate with an inhibition constant in the sub-μM range. The Van der Waals interaction was found to be more important than the electrostatic interaction in the binding affinity of this compound to NSP16-NSP10. Further in vitro and in vivo studies are needed to test the activity of the identified compound against COVID-19.Communicated by Ramaswamy H. Sarma.

Structural basis for RNA-cap recognition and methylation by the mpox methyltransferase VP39

Mpox is a zoonotic disease caused by the mpox virus (MPXV), which has gained attention due to its rapid and widespread transmission, with reports from more than 100 countries. The virus belongs to the Orthopoxvirus genus, which also includes variola virus and vaccinia virus. In poxviruses, the RNA cap is crucial for the translation and stability of viral mRNAs and also for immune evasion. This study presents the crystal structure of the mpox 2'-O-methyltransfarase VP39 in complex with a short cap-0 RNA. The RNA substrate binds to the protein without causing any significant changes to its overall fold and is held in place by a combination of electrostatic interactions, π-π stacking and hydrogen bonding. The structure also explains the mpox VP39 preference for a guanine base at the first position; it reveals that guanine forms a hydrogen bond that an adenine would not be able to form.

Accelerating antiviral drug discovery: lessons from COVID-19

During the coronavirus disease 2019 (COVID-19) pandemic, a wave of rapid and collaborative drug discovery efforts took place in academia and industry, culminating in several therapeutics being discovered, approved and deployed in a 2-year time frame. This article summarizes the collective experience of several pharmaceutical companies and academic collaborations that were active in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) antiviral discovery. We outline our opinions and experiences on key stages in the small-molecule drug discovery process: target selection, medicinal chemistry, antiviral assays, animal efficacy and attempts to pre-empt resistance. We propose strategies that could accelerate future efforts and argue that a key bottleneck is the lack of quality chemical probes around understudied viral targets, which would serve as a starting point for drug discovery. Considering the small size of the viral proteome, comprehensively building an arsenal of probes for proteins in viruses of pandemic concern is a worthwhile and tractable challenge for the community.

AI-guided pipeline for protein-protein interaction drug discovery identifies a SARS-CoV-2 inhibitor

Protein-protein interactions (PPIs) offer great opportunities to expand the druggable proteome and therapeutically tackle various diseases, but remain challenging targets for drug discovery. Here, we provide a comprehensive pipeline that combines experimental and computational tools to identify and validate PPI targets and perform early-stage drug discovery. We have developed a machine learning approach that prioritizes interactions by analyzing quantitative data from binary PPI assays and AlphaFold-Multimer predictions. Using the quantitative assay LuTHy together with our machine learning algorithm, we identified high-confidence interactions among SARS-CoV-2 proteins for which we predicted three-dimensional structures using AlphaFold Multimer. We employed VirtualFlow to target the contact interface of the NSP10-NSP16 SARS-CoV-2 methyltransferase complex by ultra-large virtual drug screening. Thereby, we identified a compound that binds to NSP10 and inhibits its interaction with NSP16, while also disrupting the methyltransferase activity of the complex, and SARS-CoV-2 replication. Overall, this pipeline will help to prioritize PPI targets to accelerate the discovery of early-stage drug candidates targeting protein complexes and pathways.

Design and synthesis of naturally-inspired SARS-CoV-2 inhibitors

A naturally inspired chemical library of 25 molecules was synthesised guided by 3-D dimensionality and natural product likeness factors to explore a new chemical space. The synthesised chemical library, consisting of fused-bridged dodecahydro-2a,6-epoxyazepino[3,4,5-c,d]indole skeletons, followed lead likeness factors in terms of molecular weight, C-sp3 fraction and Clog P. Screening of the 25 compounds against lung cells infected with SARS-CoV-2 led to the identification of 2 hits. Although the chemical library showed cytotoxicity, the two hits (3b, 9e) showed the highest antiviral activity (EC50 values of 3.7 and 1.4 μM, respectively) with an acceptable cytotoxicity difference. Computational analysis based on docking and molecular dynamics simulations against main protein targets in SARS-CoV-2 (main protease Mpro, nucleocapsid phosphoprotein, non-structural protein nsp10-nsp16 complex and RBD/ACE2 complex) were performed. The computational analysis proposed the possible binding targets to be either Mpro or the nsp10-nsp16 complex. Biological assays were performed to confirm this proposition. A cell-based assay for Mpro protease activity using a reverse-nanoluciferase (Rev-Nluc) reporter confirmed that 3b targets Mpro. These results open the way towards further hit-to-lead optimisations.

Molecular docking with SARS-CoV-2 and potential drug property of a bioactive novel Zn(II) polymer: A combined experimental and theoretical study

A new Zn(II) coordination polymer based on o-phthalato (Phth) and 2-aminopyridine (2-Ampy) viz. {[Zn(2-Ampy)₂(Phth)]∙(H₂O)]}_n (1) has been synthesized at room temperature and characterized by elemental analyses, electronic spectroscopy, FT-IR spectroscopy, thermal analysis (TGA/DSC), powder X-ray diffraction (PXRD) and single crystal X-ray diffraction. The basic trimeric units of 1 form a polymeric chain by N-H⋯O and π⋯π interactions. These polymeric chains interconnect through various non-covalent interactions in two perpendicular directions to ultimately give rise to a 3D architecture of 1. The interesting non-covalent interactions in 1, contributing to its stability in the solid state are studied by Hirshfeld surface analysis and other different theoretical tools. Molecular docking study of 1 is performed against six different proteins of SARS-CoV-2. The drug potential of the synthesized compound is evaluated by ADMET calculations.

Targeting viral proteins for restraining SARS-CoV-2: focusing lens on viral proteins beyond spike for discovering new drug targets

INTRODUCTION

Emergence of highly infectious SARS-CoV-2 variants are reducing protection provided by current vaccines, requiring constant updates in antiviral approaches. The virus encodes four structural and sixteen nonstructural proteins which play important roles in viral genome replication and transcription, virion assembly, release , entry into cells, and compromising host cellular defenses. As alien proteins to host cells, many viral proteins represent potential targets for combating the SARS-CoV-2.

AREAS COVERED

Based on literature from PubMed and Web of Science databases, the authors summarize the typical characteristics of SARS-CoV-2 from the whole viral particle to the individual viral proteins and their corresponding functions in virus life cycle. The authors also discuss the potential and emerging targeted interventions to curb virus replication and spread in detail to provide unique insights into SARS-CoV-2 infection and countermeasures against it.

EXPERT OPINION

Our comprehensive analysis highlights the rationale to focus on non-spike viral proteins that are less mutated but have important functions. Examples of this include: structural proteins (e.g. nucleocapsid protein, envelope protein) and extensively-concerned nonstructural proteins (e.g. NSP3, NSP5, NSP12) along with the ones with relatively less attention (e.g. NSP1, NSP10, NSP14 and NSP16), for developing novel drugs to overcome resistance of SARS-CoV-2 variants to preexisting vaccines and antibody-based treatments.

SARS-CoV-2 Uses Nonstructural Protein 16 To Evade Restriction by IFIT1 and IFIT3

Understanding the molecular basis of innate immune evasion by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is an important consideration for designing the next wave of therapeutics. Here, we investigate the role of the nonstructural protein 16 (NSP16) of SARS-CoV-2 in infection and pathogenesis. NSP16, a ribonucleoside 2'-O-methyltransferase (MTase), catalyzes the transfer of a methyl group to mRNA as part of the capping process. Based on observations with other CoVs, we hypothesized that NSP16 2'-O-MTase function protects SARS-CoV-2 from cap-sensing host restriction. Therefore, we engineered SARS-CoV-2 with a mutation that disrupts a conserved residue in the active site of NSP16. We subsequently show that this mutant is attenuated both in vitro and in vivo, using a hamster model of SARS-CoV-2 infection. Mechanistically, we confirm that the NSP16 mutant is more sensitive than wild-type SARS-CoV-2 to type I interferon (IFN-I) in vitro. Furthermore, silencing IFIT1 or IFIT3, IFN-stimulated genes that sense a lack of 2'-O-methylation, partially restores fitness to the NSP16 mutant. Finally, we demonstrate that sinefungin, an MTase inhibitor that binds the catalytic site of NSP16, sensitizes wild-type SARS-CoV-2 to IFN-I treatment and attenuates viral replication. Overall, our findings highlight the importance of SARS-CoV-2 NSP16 in evading host innate immunity and suggest a target for future antiviral therapies. IMPORTANCE Similar to other coronaviruses, disruption of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) NSP16 function attenuates viral replication in a type I interferon-dependent manner. In vivo, our results show reduced disease and viral replication at late times in the hamster lung, but an earlier titer deficit for the NSP16 mutant (dNSP16) in the upper airway. In addition, our results confirm a role for IFIT1 but also demonstrate the necessity of IFIT3 in mediating dNSP16 attenuation. Finally, we show that targeting NSP16 activity with a 2'-O-methyltransferase inhibitor in combination with type I interferon offers a novel avenue for antiviral development.

Lessons Learnt from COVID-19: Computational Strategies for Facing Present and Future Pandemics

Since its outbreak in December 2019, the COVID-19 pandemic has caused the death of more than 6.5 million people around the world. The high transmissibility of its causative agent, the SARS-CoV-2 virus, coupled with its potentially lethal outcome, provoked a profound global economic and social crisis. The urgency of finding suitable pharmacological tools to tame the pandemic shed light on the ever-increasing importance of computer simulations in rationalizing and speeding up the design of new drugs, further stressing the need for developing quick and reliable methods to identify novel active molecules and characterize their mechanism of action. In the present work, we aim at providing the reader with a general overview of the COVID-19 pandemic, discussing the hallmarks in its management, from the initial attempts at drug repurposing to the commercialization of Paxlovid, the first orally available COVID-19 drug. Furthermore, we analyze and discuss the role of computer-aided drug discovery (CADD) techniques, especially those that fall in the structure-based drug design (SBDD) category, in facing present and future pandemics, by showcasing several successful examples of drug discovery campaigns where commonly used methods such as docking and molecular dynamics have been employed in the rational design of effective therapeutic entities against COVID-19.

Masked inverse folding with sequence transfer for protein representation learning

Self-supervised pretraining on protein sequences has led to state-of-the art performance on protein function and fitness prediction. However, sequence-only methods ignore the rich information contained in experimental and predicted protein structures. Meanwhile, inverse folding methods reconstruct a protein's amino-acid sequence given its structure, but do not take advantage of sequences that do not have known structures. In this study, we train a masked inverse folding protein masked language model parameterized as a structured graph neural network. During pretraining, this model learns to reconstruct corrupted sequences conditioned on the backbone structure. We then show that using the outputs from a pretrained sequence-only protein masked language model as input to the inverse folding model further improves pretraining perplexity. We evaluate both of these models on downstream protein engineering tasks and analyze the effect of using information from experimental or predicted structures on performance.

Analysis of critical protein-protein interactions of SARS-CoV-2 capping and proofreading molecular machineries towards designing dual target inhibitory peptides

In recent years, the emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), as the cause of the coronavirus disease (COVID-19) global pandemic, and its variants, especially those with higher transmissibility and substantial immune evasion, have highlighted the imperative for developing novel therapeutics as sustainable solutions other than vaccination to combat coronaviruses (CoVs). Beside receptor recognition and virus entry, members of the SARS-CoV-2 replication/transcription complex are promising targets for designing antivirals. Here, the interacting residues that mediate protein-protein interactions (PPIs) of nsp10 with nsp16 and nsp14 were comprehensively analyzed, and the key residues' interaction maps, interaction energies, structural networks, and dynamics were investigated. Nsp10 stimulates both nsp14's exoribonuclease (ExoN) and nsp16's 2'O-methyltransferase (2'O-MTase). Nsp14 ExoN is an RNA proofreading enzyme that supports replication fidelity. Nsp16 2'O-MTase is responsible for the completion of RNA capping to ensure efficient replication and translation and escape from the host cell's innate immune system. The results of the PPIs analysis proposed crucial information with implications for designing SARS-CoV-2 antiviral drugs. Based on the predicted shared protein-protein interfaces of the nsp16-nsp10 and nsp14-nsp10 interactions, a set of dual-target peptide inhibitors was designed. The designed peptides were evaluated by molecular docking, peptide-protein interaction analysis, and free energy calculations, and then further optimized by in silico saturation mutagenesis. Based on the predicted evolutionary conservation of the interacted target residues among CoVs, the designed peptides have the potential to be developed as dual target pan-coronavirus inhibitors.

A computational study on the molecular mechanisms of panduratin A as a potential inhibitor on SARS-CoV-2 protein targets

Panduratin A from Boesebergia rotunda was recently reported as a potent anti-SARS-CoV-2 compound. However, the molecular mechanisms underlying the inhibition by Panduratin A and its target remained unclear. Molecular docking calculations were performed between panduratin A and five important proteins, i.e., main protease (Mpro), papain-like protease (PLpro), receptor binding domain (RBD) of spike proteins, RNA-dependent-RNA-polymerase (RdRp), and 2'-O-methyltransferase (MTase). The estimated binding free energy and the interaction networks extracted from the best docking mode for each complex suggested that MTase was the most probable target for panduratin A inhibition. To further validate the ability of panduratin A to inhibit MTase, molecular dynamics (MD) simulations and binding free energy calculations were performed for panduratin A-MTase complex, in comparison with another MTase complex with sinefungin as a positive control. Chemical features of panduratin A and sinefungin were compared for their contribution in MTase binding. It was found that both molecules could bind to the S-Adenosyl methionine (SAM) binding pocket and prevent the SAM entrance co-substrate, which could eventually halt the function of MTase. Despite a slightly weaker binding free energy, the equilibrated positional binding of panduratin A was found at a closer distance to the active sites. Therefore, this study proposed MTase as a possible target of panduratin A, along with the mechanisms of inhibition, prompting another future in vitro study as a verification.

Potential targets of severe acute respiratory syndrome coronavirus 2 of clinical drug fluvoxamine: Docking and molecular dynamics studies to elucidate viral action

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has continued evolving for survival and adaptation by mutating itself into different variants of concern, including omicron. Several studies and clinical trials found fluvoxamine, an Food and Drug Administration-approved antidepressant drug, to be effective at preventing mild coronavirus disease 2019 (COVID-19) from progressing to severe diseases. However, the mechanism of fluvoxamine's direct antiviral action against COVID-19 is still unknown. Fluvoxamine was docked with 11 SARS-CoV-2 targets and subjected to stability, conformational changes, and binding free energy analyses to explore its mode of action. Of the targets, nonstructural protein 14 (NSP14), main protease (Mpro), and papain-like protease (PLpro) had the best docking scores with fluvoxamine. Consistent with the docking results, it was confirmed by molecular dynamics simulations that the NSP14 N7-MTase ((N7-guanine)-methyltransferase)-fluvoxamine, Mpro-fluvoxamine, and PLpro-fluvoxamine complexes are stable, with the lowest binding free energies of -105.1, -82.7, and - 38.5 kJ/mol, respectively. A number of hotspot residues involved in the interaction were also identified. These include Glu166, Asp187, His41, and Cys145 in Mpro, Gly163 and Arg166 in PLpro, and Glu302, Gly333, and Phe426 in NSP14, which could aid in the development of better antivirals against SARS-CoV-2.

Structural screening into the recognition of a potent inhibitor against non-structural protein 16: a molecular simulation to inhibit SARS-CoV-2 infection

COVID-19 infection is caused by endemic crown infection (SARS-CoV-2) and is associated with lung damage and severe immune response. Non-Structural Proteins are the central components of coronaviral transcription and replication machinery in SARS-CoV-2 and also stimulate mRNA cap methylation to avoid the immune response. Non-Structural Protein 16 (NSP16) is one of the primary targets for the drug discovery of coronaviruses. Discovering an effective inhibitor against the NSP16 in comparison with Sinefungin was the main purpose of this investigation. Binding free-energy calculations, computational methods of molecular dynamics, docking, and virtual screening were utilized in this study. The ZINC and PubChem databases were applied to screen some chemical compounds regarding Sinefungin as a control inhibitor. Based on structural similarity to Sinefungin, 355 structures were obtained from the mentioned databases. Subsequently, this set of compounds were monitored by AutoDock Vina software, and ultimately the potent inhibitor (PUBCHEM512713) was chosen. At the next stage, molecular dynamics were carried out by GROMACS software to evaluate the potential elected compounds in a simulated environment and in a timescale of 100 nanoseconds. MM-PBSA investigation exhibited that the value of binding free energy for PUBCHEM512713 (-30.829 kJ.mol^-1) is more potent than Sinefungin (-11.941 kJ.mol^-1). Furthermore, the results of ADME analysis illustrated that the pharmacokinetics, drug-likeness, and lipophilicity parameters of PUBCHEM512713 are admissible for human utilization. Finally, our data suggested that PUBCHEM512713 is an effective drug candidate for inhibiting the NSP16 and is suitable for in vitro and in vivo studies.Communicated by Ramaswamy H. Sarma.

Crystal structure and cap binding analysis of the methyltransferase of langat virus

Tick-borne encephalitis virus (TBEV) is a major dangerous human pathogen, as TBEV infection can cause serious illness that can lead to irreversible neurological sequelae and even death. Langat virus (LGTV), a member of the tick-borne encephalitis virus (TBEV) serogroup, belongs to the family Flaviviridae, genus Flavivirus. Its nonstructural protein 5 (NS5) protein contains a methyltransferase (MTase) domain that can methylate RNA cap structures, which is critical for viral replication. We determined the structure of LGTV NS5 methyltransferase bound to S-adenosyl-L-homocysteine (SAH) at a 1.70 Å resolution. Sequence analysis and structural comparison of homologous MTases suggests that folds and structures are closely conserved throughout Flavivirus species and play important roles. This study provides the key structural information on LGTV MTase and the foundation for research on antiviral drugs targeting LGTV MTase.

Multi-target direct-acting SARS-CoV-2 antivirals against the nucleotide-binding pockets of virus-specific proteins

The nucleotide-binding pockets (NBPs) in virus-specific proteins have proven to be the most successful antiviral targets for several viral diseases. Functionally important NBPs are found in various structural and non-structural proteins of SARS-CoV-2. In this study, the first successful multi-targeting attempt to identify effective antivirals has been made against NBPs in nsp12, nsp13, nsp14, nsp15, nsp16, and nucleocapsid (N) proteins of SARS-CoV-2. A structure-based drug repurposing in silico screening approach with ADME analysis identified small molecules targeting NBPs in SARS-CoV-2 proteins. Further, isothermal titration calorimetry (ITC) experiments validated the binding of top hit molecules to the purified N-protein. Importantly, cell-based antiviral assays revealed antiviral potency for INCB28060, darglitazone, and columbianadin with EC50 values 15.71 μM, 5.36 μM, and 22.52 μM, respectively. These effective antivirals targeting multiple proteins are envisioned to direct the development of antiviral therapy against SARS-CoV-2 and its emerging variants.

Identification and characterization of alternative sites and molecular probes for SARS-CoV-2 target proteins

Three protein targets from SARS-CoV-2, the viral pathogen that causes COVID-19, are studied: the main protease, the 2′-O-RNA methyltransferase, and the nucleocapsid (N) protein. For the main protease, the nucleophilicity of the catalytic cysteine C145 is enabled by coupling to three histidine residues, H163 and H164 and catalytic dyad partner H41. These electrostatic couplings enable significant population of the deprotonated state of C145. For the RNA methyltransferase, the catalytic lysine K6968 that serves as a Brønsted base has significant population of its deprotonated state via strong coupling with K6844 and Y6845. For the main protease, Partial Order Optimum Likelihood (POOL) predicts two clusters of biochemically active residues; one includes the catalytic H41 and C145 and neighboring residues. The other surrounds a second pocket adjacent to the catalytic site and includes S1 residues F140, L141, H163, E166, and H172 and also S2 residue D187. This secondary recognition site could serve as an alternative target for the design of molecular probes. From in silico screening of library compounds, ligands with predicted affinity for the secondary site are reported. For the NSP16-NSP10 complex that comprises the RNA methyltransferase, three different sites are predicted. One is the catalytic core at the conserved K-D-K-E motif that includes catalytic residues D6928, K6968, and E7001 plus K6844. The second site surrounds the catalytic core and consists of Y6845, C6849, I6866, H6867, F6868, V6894, D6895, D6897, I6926, S6927, Y6930, and K6935. The third is located at the heterodimer interface. Ligands predicted to have high affinity for the first or second sites are reported. Three sites are also predicted for the nucleocapsid protein. This work uncovers key interactions that contribute to the function of the three viral proteins and also suggests alternative sites for ligand design.

Nsp16 shields SARS-CoV-2 from efficient MDA5 sensing and IFIT1-mediated restriction

Methylation of the mRNA 5' cap by cellular methyltransferases enables efficient translation and avoids recognition by innate immune factors. Coronaviruses encode viral 2'-O-methyltransferases to shield their RNA from host factors. Here, we generate recombinant SARS-CoV-2 harboring a catalytically inactive 2'-O-methyltransferase Nsp16, Nsp16mut, and analyze viral replication in human lung epithelial cells. Although replication is only slightly attenuated, we find SARS-CoV-2 Nsp16mut to be highly immunogenic, resulting in a strongly enhanced release of type I interferon upon infection. The elevated immunogenicity of Nsp16mut is absent in cells lacking the RNA sensor MDA5. In addition, we report that Nsp16mut is highly sensitive to type I IFN treatment and demonstrate that this strong antiviral effect of type I IFN is mediated by the restriction factor IFIT1. Together, we describe a dual role for the 2'-O-methyltransferase Nsp16 during SARS-CoV-2 replication in avoiding efficient recognition by MDA5 and in shielding its RNA from interferon-induced antiviral responses, thereby identifying Nsp16 as a promising target for generating attenuated and highly immunogenic SARS-CoV-2 strains and as a potential candidate for therapeutic intervention.

Recent Drug Development and Medicinal Chemistry Approaches for the Treatment of SARS-CoV-2 Infection and COVID-19

COVID-19, caused by SARS-CoV-2 infection, continues to be a major public health crisis around the globe. Development of vaccines and the first cluster of antiviral drugs has brought promise and hope for prevention and treatment of severe coronavirus disease. However, continued development of newer, safer, and more effective antiviral drugs are critically important to combat COVID-19 and counter the looming pathogenic variants. Studies of the coronavirus life cycle revealed several important biochemical targets for drug development. In the present review, we focus on recent drug design and medicinal chemistry efforts in small molecule drug discovery, including the development of nirmatrelvir that targets viral protein synthesis and remdesivir and molnupiravir that target viral RdRp. These are recent FDA approved drugs for the treatment of COVID-19.

SARS-CoV-2 Uses Nonstructural Protein 16 to Evade Restriction by IFIT1 and IFIT3

Understanding the molecular basis of innate immune evasion by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is an important consideration for designing the next wave of therapeutics. Here, we investigate the role of the nonstructural protein 16 (NSP16) of SARS-CoV-2 in infection and pathogenesis. NSP16, a ribonucleoside 2'- O methyltransferase (MTase), catalyzes the transfer of a methyl group to mRNA as part of the capping process. Based on observations with other CoVs, we hypothesized that NSP16 2'- O MTase function protects SARS-CoV-2 from cap-sensing host restriction. Therefore, we engineered SARS-CoV-2 with a mutation that disrupts a conserved residue in the active site of NSP16. We subsequently show that this mutant is attenuated both in vitro and in vivo , using a hamster model of SARS-CoV-2 infection. Mechanistically, we confirm that the NSP16 mutant is more sensitive to type I interferon (IFN-I) in vitro . Furthermore, silencing IFIT1 or IFIT3, IFN-stimulated genes that sense a lack of 2'- O methylation, partially restores fitness to the NSP16 mutant. Finally, we demonstrate that sinefungin, a methyltransferase inhibitor that binds the catalytic site of NSP16, sensitizes wild-type SARS-CoV-2 to IFN-I treatment. Overall, our findings highlight the importance of SARS-CoV-2 NSP16 in evading host innate immunity and suggest a possible target for future antiviral therapies.

Importance

Similar to other coronaviruses, disruption of SARS-CoV-2 NSP16 function attenuates viral replication in a type I interferon-dependent manner. In vivo , our results show reduced disease and viral replication at late times in the hamster lung, but an earlier titer deficit for the NSP16 mutant (dNSP16) in the upper airway. In addition, our results confirm a role for IFIT1, but also demonstrate the necessity of IFIT3 in mediating dNSP16 attenuation. Finally, we show that targeting NSP16 activity with a 2'- O methyltransferase inhibitor in combination with type I interferon offers a novel avenue for antiviral development.

Toll-like Receptor Mediation in SARS-CoV-2: A Therapeutic Approach

The innate immune system facilitates defense mechanisms against pathogen invasion and cell damage. Toll-like receptors (TLRs) assist in the activation of the innate immune system by binding to pathogenic ligands. This leads to the generation of intracellular signaling cascades including the biosynthesis of molecular mediators. TLRs on cell membranes are adept at recognizing viral components. Viruses can modulate the innate immune response with the help of proteins and RNAs that downregulate or upregulate the expression of various TLRs. In the case of COVID-19, molecular modulators such as type 1 interferons interfere with signaling pathways in the host cells, leading to an inflammatory response. Coronaviruses are responsible for an enhanced immune signature of inflammatory chemokines and cytokines. TLRs have been employed as therapeutic agents in viral infections as numerous antiviral Food and Drug Administration-approved drugs are TLR agonists. This review highlights the therapeutic approaches associated with SARS-CoV-2 and the TLRs involved in COVID-19 infection.

SARS-CoV-2 Non-Structural Proteins and Their Roles in Host Immune Evasion

Coronavirus disease 2019 (COVID-19) has caused an unprecedented global crisis and continues to threaten public health. The etiological agent of this devastating pandemic outbreak is the severe acute respiratory syndrome-coronavirus-2 (SARS-CoV-2). COVID-19 is characterized by delayed immune responses, followed by exaggerated inflammatory responses. It is well-established that the interferon (IFN) and JAK/STAT signaling pathways constitute the first line of defense against viral and bacterial infections. To achieve viral replication, numerous viruses are able to antagonize or hijack these signaling pathways to attain productive infection, including SARS-CoV-2. Multiple studies document the roles of several non-structural proteins (NSPs) of SARS-CoV-2 that facilitate the establishment of viral replication in host cells via immune escape. In this review, we summarize and highlight the functions and characteristics of SARS-CoV-2 NSPs that confer host immune evasion. The molecular mechanisms mediating immune evasion and the related potential therapeutic strategies for controlling the COVID-19 pandemic are also discussed.

Crystal structure of SARS‐CoV ‐2 nsp10–nsp16 in complex with small molecule inhibitors, SS148 and WZ16

SARS‐CoV‐2 nsp10–nsp16 complex is a 2′‐O‐methyltransferase (MTase) involved in viral RNA capping, enabling the virus to evade the immune system in humans. It has been considered a valuable target in the discovery of antiviral therapeutics, as the RNA cap formation is crucial for viral propagation. Through cross‐screening of the inhibitors that we previously reported for SARS‐CoV‐2 nsp14 MTase activity against nsp10–nsp16 complex, we identified two compounds (SS148 and WZ16) that also inhibited nsp16 MTase activity. To further enable the chemical optimization of these two compounds towards more potent and selective dual nsp14/nsp16 MTase inhibitors, we determined the crystal structure of nsp10–nsp16 in complex with each of SS148 and WZ16. As expected, the structures revealed the binding of both compounds to S‐adenosyl‐L‐methionine (SAM) binding pocket of nsp16. However, our structural data along with the biochemical mechanism of action determination revealed an RNA‐dependent SAM‐competitive pattern of inhibition for WZ16, clearly suggesting that binding of the RNA first may help the binding of some SAM competitive inhibitors. Both compounds also showed some degree of selectivity against human protein MTases, an indication of great potential for chemical optimization towards more potent and selective inhibitors of coronavirus MTases.

PDB Code(s): 7R1T and 7R1U

Virtual screening and molecular dynamics simulation for identification of natural antiviral agents targeting SARS-CoV-2 NSP10

New variations of SARS-CoV-2 continue to emerge in the global pandemic, which may be resistant to at least some vaccines in COVID-19, indicating that drug and vaccine development must be continuously strengthened. NSP10 plays an essential role in SARS-CoV-2 viral life cycle. It stimulates the enzymatic activities of NSP14-ExoN and NSP16-O-MTase by the formation of NSP10/NSP14 and NSP10/NSP16 complexes. Inhibiting NSP10 can block the binding of NSP10 to NSP14 and NSP16. This study has identified potential natural NSP10 inhibitors from ZINC database. The protein druggable pocket was identified for screening candidates. Molecular docking of the selected compounds was performed and MM-GBSA binding energy was calculated. After ADMET assessment, 4 hits were obtained for favorable druggability. The analysis of site interactions suggested that the hits all had excellent binding. Molecular dynamics studies revealed that selected natural compounds stably bind to NSP10. These compounds were identified as potential leads against NSP10 for the development of strategies to combat SARS-CoV-2 replication and could serve as the basis for further studies.

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

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

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

Interferon induction, evasion, and paradoxical roles during SARS-CoV-2 infection

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causative agent of coronavirus disease 2019 (COVID-19), has caused millions of deaths in the past two years. Although initially little was understood about this virus, recent research has significantly advanced and landed interferons (IFNs) in the spotlight. While Type I and III IFN have long been known as central to antiviral immunity, in the case of COVID-19 their role was initially controversial. However, the protective function of IFN is now well supported by the identification of human deficiencies in IFN responses as a predictor of disease severity. Here, we will review the cell types and pathways that lead to IFN production as well as the importance of IFN timing and location for disease outcome. We will further discuss the mechanisms that SARS-CoV-2 uses to evade IFN responses, and the current efforts to implement IFNs as therapeutics in the treatment of COVID-19. It is essential to understand the relationships between SARS-CoV-2 and IFN to better inform treatments that exploit IFN functions to alleviate COVID-19.

In-silico screening and in-vitro assay show the antiviral effect of Indomethacin against SARS-CoV-2

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has caused the worldwide spread of coronavirus disease 19 (COVID-19), and till now, it has caused death to more than 6.2 million people. Although various vaccines and drug candidates are being tested globally with limited to moderate success, a comprehensive therapeutic cure is yet to be achieved. In this study, we applied computational drug repurposing methods complemented with the analyses of the already existing gene expression data to find better therapeutics in treatment and recovery. Primarily, we identified the most crucial proteins of SARS-CoV-2 and host human cells responsible for viral infection and host response. An in-silico screening of the existing drugs was performed against the crucial proteins for SARS-CoV-2 infection, and a few existing drugs were shortlisted. Further, we analyzed the gene expression data of SARS-CoV-2 in human lung epithelial cells and investigated the molecules that can reverse the cellular mRNA expression profiles in the diseased state. LINCS L1000 and Comparative Toxicogenomics Database (CTD) were utilized to obtain two sets of compounds that can be used to counter SARS-CoV-2 infection from the gene expression perspective. Indomethacin, a nonsteroidal anti-inflammatory drug (NSAID), and Vitamin-A were found in two sets of compounds, and in the in-silico screening of existing drugs to treat SARS-CoV-2. Our in-silico findings on Indomethacin were further successfully validated by in-vitro testing in Vero CCL-81 cells with an IC₅₀ of 12 μM. Along with these findings, we briefly discuss the possible roles of Indomethacin and Vitamin-A to counter the SARS-CoV-2 infection in humans.

Attenuation of SARS-CoV-2 replication and associated inflammation by concomitant targeting of viral and host cap 2'-O-ribose methyltransferases

The SARS‐CoV‐2 infection cycle is a multistage process that relies on functional interactions between the host and the pathogen. Here, we repurposed antiviral drugs against both viral and host enzymes to pharmaceutically block methylation of the viral RNA 2'‐O‐ribose cap needed for viral immune escape. We find that the host cap 2'‐O‐ribose methyltransferase MTr1 can compensate for loss of viral NSP16 methyltransferase in facilitating virus replication. Concomitant inhibition of MTr1 and NSP16 efficiently suppresses SARS‐CoV‐2 replication. Using in silico target‐based drug screening, we identify a bispecific MTr1/NSP16 inhibitor with anti‐SARS‐CoV‐2 activity in vitro and in vivo but with unfavorable side effects. We further show antiviral activity of inhibitors that target independent stages of the host SAM cycle providing the methyltransferase co‐substrate. In particular, the adenosylhomocysteinase (AHCY) inhibitor DZNep is antiviral in in vitro, in ex vivo, and in a mouse infection model and synergizes with existing COVID‐19 treatments. Moreover, DZNep exhibits a strong immunomodulatory effect curbing infection‐induced hyperinflammation and reduces lung fibrosis markers ex vivo. Thus, multispecific and metabolic MTase inhibitors constitute yet unexplored treatment options against COVID‐19.

Methodology-Centered Review of Molecular Modeling, Simulation, and Prediction of SARS-CoV-2

Despite tremendous efforts in the past two years, our understanding of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), virus-host interactions, immune response, virulence, transmission, and evolution is still very limited. This limitation calls for further in-depth investigation. Computational studies have become an indispensable component in combating coronavirus disease 2019 (COVID-19) due to their low cost, their efficiency, and the fact that they are free from safety and ethical constraints. Additionally, the mechanism that governs the global evolution and transmission of SARS-CoV-2 cannot be revealed from individual experiments and was discovered by integrating genotyping of massive viral sequences, biophysical modeling of protein-protein interactions, deep mutational data, deep learning, and advanced mathematics. There exists a tsunami of literature on the molecular modeling, simulations, and predictions of SARS-CoV-2 and related developments of drugs, vaccines, antibodies, and diagnostics. To provide readers with a quick update about this literature, we present a comprehensive and systematic methodology-centered review. Aspects such as molecular biophysics, bioinformatics, cheminformatics, machine learning, and mathematics are discussed. This review will be beneficial to researchers who are looking for ways to contribute to SARS-CoV-2 studies and those who are interested in the status of the field.

Structural basis for SARS-CoV-2 nucleocapsid (N) protein recognition by 14-3-3 proteins

14-3-3 proteins are important dimeric scaffolds that regulate the function of hundreds of proteins in a phosphorylation-dependent manner. The SARS-CoV-2 nucleocapsid (N) protein forms a complex with human 14-3-3 proteins upon phosphorylation, which has also been described for other coronaviruses. Here, we report a high-resolution crystal structure of 14-3-3 bound to an N phosphopeptide bearing the phosphoserine 197 in the middle. The structure revealed two copies of the N phosphopeptide bound, each in the central binding groove of each 14-3-3 monomer. A complex network of hydrogen bonds and water bridges between the peptide and 14-3-3 was observed explaining the high affinity of the N protein for 14-3-3 proteins.

The Robustness of Cellular Immunity Determines the Fate of SARS-CoV-2 Infection

Two years after the appearance of the SARS-CoV-2 virus, the causal agent of the current global pandemic, it is time to analyze the evolution of the immune protection that infection and vaccination provide. Cellular immunity plays an important role in limiting disease severity and the resolution of infection. The early appearance, breadth and magnitude of SARS-CoV-2 specific T cell response has been correlated with disease severity and it has been thought that T cell responses may be sufficient to clear infection with minimal disease in COVID-19 patients with X-linked or autosomal recessive agammaglobulinemia. However, our knowledge of the phenotypic and functional diversity of CD8+ cytotoxic lymphocytes, CD4+ T helper cells, mucosal-associated invariant T (MAIT) cells and CD4+ T follicular helper (Tfh), which play a critical role in infection control as well as long-term protection, is still evolving. It has been described how CD8+ cytotoxic lymphocytes interrupt viral replication by secreting antiviral cytokines (IFN-γ and TNF-α) and directly killing infected cells, negatively correlating with stages of disease progression. In addition, CD4+ T helper cells have been reported to be key pieces, leading, coordinating and ultimately regulating antiviral immunity. For instance, in some more severe COVID-19 cases a dysregulated CD4+ T cell signature may contribute to the greater production of pro-inflammatory cytokines responsible for pathogenic inflammation. Here we discuss how cellular immunity is the axis around which the rest of the immune system components revolve, since it orchestrates and leads antiviral response by regulating the inflammatory cascade and, as a consequence, the innate immune system, as well as promoting a correct humoral response through CD4+ Tfh cells. This review also analyses the critical role of cellular immunity in modulating the development of high-affinity neutralizing antibodies and germinal center B cell differentiation in memory and long-lived antibody secreting cells. Finally, since there is currently a high percentage of vaccinated population and, in some cases, vaccine booster doses are even being administered in certain countries, we have also summarized newer approaches to long-lasting protective immunity and the cross-protection of cellular immune response against SARS-CoV-2.

Conserved coronavirus proteins as targets of broad-spectrum antivirals

Coronaviruses are a class of single-stranded, positive-sense RNA viruses that have caused three major outbreaks over the past two decades: Middle East respiratory syndrome–related coronavirus (MERS-CoV), severe acute respiratory syndrome coronavirus (SARS-CoV), and severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). All outbreaks have been associated with significant morbidity and mortality. In this study, we have identified and explored conserved binding sites in the key coronavirus proteins for the development of broad-spectrum direct acting anti-coronaviral compounds and validated the significance of this conservation for drug discovery with existing experimental data. We have identified four coronaviral proteins with highly conserved binding site sequence and 3D structure similarity: PL^pro, M^pro, nsp10-nsp16 complex(methyltransferase), and nsp15 endoribonuclease. We have compiled all available experimental data for known antiviral medications inhibiting these targets and identified compounds active against multiple coronaviruses. The identified compounds representing potential broad-spectrum antivirals include: GC376, which is active against six viral M^pro (out of six tested, as described in research literature); mycophenolic acid, which is active against four viral PL^pro (out of four); and emetine, which is active against four viral RdRp (out of four). The approach described in this study for coronaviruses, which combines the assessment of sequence and structure conservation across a viral family with the analysis of accessible chemical structure – antiviral activity data, can be explored for the development of broad-spectrum drugs for multiple viral families.

Data-driven platform for identifying variants of interest in COVID-19 virus

•

We provide an automated way to identify emerging variants of concern using viral genome and patient-outcome data.

•

We assembled 10,000 sample-strong case-control dataset and identified 117 single nucleotide variants (SNV) associated with adverse patient outcomes.

•

We observe co-evolution of protective and pathogenic interactions between the spike, nsp14, and N region with either orf3a or nsp3.

•

Structural modelling reveals mutation clusters in the Zn binding domain of nsp14 suggesting ongoing adaptation to the human host.

•

Our approach identified Variants Being Monitored (VBM) a week before they were flagged by Health Organizations and offers a clade-independent function-orientated grouping.

New SARS-CoV-2 variants emerge as part of the virus’ adaptation to the human host. The Health Organizations are monitoring newly emerging variants with suspected impact on disease or vaccination efficacy as Variants Being Monitored (VBM), like Delta and Omicron. Genetic changes (SNVs) compared to the Wuhan variant characterize VBMs with current emphasis on the spike protein and lineage markers. However, monitoring VBMs in such a way might miss SNVs with functional effect on disease.

Here we introduce a lineage-agnostic genome-wide approach to identify SNVs associated with disease. We curated a case-control dataset of 10,520 samples and identified 117 SNVs significantly associated with adverse patient outcome. While 40% (47) SNV are already monitored and 36% (43) are in the spike protein, we also identified 70 new SNVs that are associated with disease outcome. 31 of these are disease-worsening and predominantly located in the 3′-5′ exonuclease (NSP14) with structural modelling revealing a concise cluster in the Zn binding domain that has known host-immune modulating function. Furthermore, we generate clade-independent VBM groupings by identifying interacting SNVs (epistasis). We find 37 sets of higher-order epistatic interactions joining 5 genomic regions (nsp3, nsp14, Spike S1, ORF3a, N). Structural modelling of these regions provides insights into potential mechanistic pathways of increased virulence as well as orthogonal methods of validation.

Clade-independent monitoring of functionally interacting (epistasis, co-evolution) SNVs detected emerging VBM a week before they were flagged by Health Organizations and in conjunction with structural modelling provides faster, mechanistic insight into emerging strains to guide public health interventions.