De novo design of potential peptide analogs against the main protease of Omicron variant using in silico studies.

SARS-CoV-2 and its variants are crossing the immunity barrier induced through vaccination. Recent Omicron sub-variants are highly transmissible and have a low mortality rate. Despite the low severity of Omicron variants, these new variants are known to cause acute post-infectious syndromes. Nowadays, novel strategies to develop new potential inhibitors for SARS-CoV-2 and other Omicron variants have gained prominence. For viral replication and survival the main protease of SARS-CoV-2 plays a vital role. Peptide-like inhibitors that mimic the substrate peptide have already proved to be effective in inhibiting the Mpro of SARS-CoV-2 variants. Our systematic canonical amino acid point mutation analysis on the native peptide has revealed various ways to improve the native peptide of the main protease. Multi mutation analysis has led us to identify and design potent peptide-analog inhibitors that act against the Mpro of the Omicron sub-variants. Our in-depth analysis of all-atom molecular dynamics studies has paved the way to characterize the atomistic behavior of Mpro in Omicron variants. Our goal is to develop potent peptide-analogs that could be therapeutically effective against Omicron and its sub-variants.

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.

Exploring the Efficacy of Noncovalent SARS-CoV-2 Main Protease Inhibitors: A Computational Simulation Analysis Study.

The SARS-CoV-2 main protease, as a key target for antiviral therapeutics, is instrumental in maintaining virus stability, facilitating translation, and enabling the virus to evade innate immunity. Our research focused on designing non-covalent inhibitors to counteract the action of this protease. Utilizing a 3D-QSAR model and contour map, we successfully engineered eight novel non-covalent inhibitors. Further evaluation and comparison of these novel compounds through methodologies including molecular docking, ADMET analysis, frontier molecular orbital studies, molecular dynamics simulations, and binding free energy revealed that the inhibitors N02 and N03 demonstrated superior research performance (N02 ΔGbind=-206.648 kJ/mol, N03 ΔGbind=-185.602 kJ/mol). These findings offer insightful guidance for the further refinement of molecular structures and the development of more efficacious inhibitors. Consequently, future investigations can draw upon these findings to unearth more potent inhibitors, thereby amplifying their impact in the treatment and prevention of associated diseases.

Design of a SARS-CoV-2 papain-like protease inhibitor with antiviral efficacy in a mouse model.

The emergence of SARS-CoV-2 variants and drug-resistant mutants calls for additional oral antivirals. The SARS-CoV-2 papain-like protease (PLpro) is a promising but challenging drug target. We designed and synthesized 85 noncovalent PLpro inhibitors that bind to a recently discovered ubiquitin binding site and the known BL2 groove pocket near the S4 subsite. Leads inhibited PLpro with the inhibitory constant Ki values from 13.2 to 88.2 nanomolar. The co-crystal structures of PLpro with eight leads revealed their interaction modes. The in vivo lead Jun12682 inhibited SARS-CoV-2 and its variants, including nirmatrelvir-resistant strains with EC50 from 0.44 to 2.02 micromolar. Oral treatment with Jun12682 improved survival and reduced lung viral loads and lesions in a SARS-CoV-2 infection mouse model, suggesting that PLpro inhibitors are promising oral SARS-CoV-2 antiviral candidates.

Discovery of Covalent Lead Compounds Targeting 3CL Protease with a Lateral Interactions Spiking Neural Network.

Covalent drugs exhibit advantages in that noncovalent drugs cannot match, and covalent docking is an important method for screening covalent lead compounds. However, it is difficult for covalent docking to screen covalent compounds on a large scale because covalent docking requires determination of the covalent reaction type of the compound. Here, we propose to use deep learning of a lateral interactions spiking neural network to construct a covalent lead compound screening model to quickly screen covalent lead compounds. We used the 3CL protease (3CL Pro) of SARS-CoV-2 as the screen target and constructed two classification models based on LISNN to predict the covalent binding and inhibitory activity of compounds. The two classification models were trained on the covalent complex data set targeting cysteine (Cys) and the compound inhibitory activity data set targeting 3CL Pro, respected, with good prediction accuracy (ACC > 0.9). We then screened the screening compound library with 6 covalent binding screening models and 12 inhibitory activity screening models. We tested the inhibitory activity of the 32 compounds, and the best compound inhibited SARS-CoV-2 3CL Pro with an IC50 value of 369.5 nM. Further assay implied that dithiothreitol can affect the inhibitory activity of the compound to 3CL Pro, indicating that the compound may covalently bind 3CL Pro. The selectivity test showed that the compound had good target selectivity to 3CL Pro over cathepsin L. These correlation assays can prove the rationality of the covalent lead compound screening model. Finally, covalent docking was performed to demonstrate the binding conformation of the compound with 3CL Pro. The source code can be obtained from the GitHub repository (https://github.com/guzh970630/Screen_Covalent_Compound_by_LISNN).

The PKA-CREB1 axis regulates coronavirus proliferation by viral helicase nsp13 association.

The COVID-19 pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has posed a worldwide threat in the past 3 years. Although it has been widely and intensively investigated, the mechanism underlying the coronavirus-host interaction requires further elucidation, which may contribute to the development of new antiviral strategies. Here, we demonstrated that the host cAMP-responsive element-binding protein (CREB1) interacts with the non-structural protein 13 (nsp13) of SARS-CoV-2, a conserved helicase for coronavirus replication, both in cells and in lung tissues subjected to SARS-CoV-2 infection. The ATPase and helicase activity of viral nsp13 were shown to be potentiated by CREB1 association, as well as by Protein kinase A (PKA)-mediated CREB1 activation. SARS-CoV-2 replication is significantly suppressed by PKA Cα, cAMP-activated protein kinase catalytic subunit alpha (PRKACA), and CREB1 knockdown or inhibition. Consistently, the CREB1 inhibitor 666-15 has shown significant antiviral effects against both the WIV04 strain and the Omicron strain of the SARS-CoV-2. Our findings indicate that the PKA-CREB1 signaling axis may serve as a novel therapeutic target against coronavirus infection. IMPORTANCE: In this study, we provide solid evidence that host transcription factor cAMP-responsive element-binding protein (CREB1) interacts directly with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) helicase non-structural protein 13 (nsp13) and potentiate its ATPase and helicase activity. And by live SARS-CoV-2 virus infection, the inhibition of CREB1 dramatically impairs SARS-CoV-2 replication in vivo. Notably, the IC50 of CREB1 inhibitor 666-15 is comparable to that of remdesivir. These results may extend to all highly pathogenic coronaviruses due to the conserved nsp13 sequences in the virus.

Structural Basis for Coronaviral Main Proteases Inhibition by the 3CLpro Inhibitor GC376.

The main protease (Mpro) of coronaviruses participates in viral replication, serving as a hot target for drug design. GC376 is able to effectively inhibit the activity of Mpro, which is due to nucleophilic addition of GC376 by binding covalently with Cys145 in Mpro active site. Here, we used fluorescence resonance energy transfer (FRET) assay to analyze the IC50 values of GC376 against Mpros from six different coronaviruses (SARS-CoV-2, HCoV-229E, HCoV-HUK1, MERS-CoV, SARS-CoV, HCoV-NL63) and five Mpro mutants (G15S, M49I, K90R, P132H, S46F) from SARS-CoV-2 variants. The results showed that GC376 displays effective inhibition to various coronaviral Mpros and SARS-CoV-2 Mpro mutants. In addition, the crystal structures of SARS-CoV-2 Mpro (wide type)-GC376, SARS-CoV Mpro-GC376, MERS-CoV Mpro-GC376, and SARS-CoV-2 Mpro mutants (G15S, M49I, S46F, K90R, and P132H)-GC376 complexes were solved. We found that GC376 is able to fit into the active site of Mpros from different coronaviruses and different SARS-CoV-2 variants properly. Detailed structural analysis revealed key molecular determinants necessary for inhibition and illustrated the binding patterns of GC376 to these different Mpros. In conclusion, we not only proved the inhibitory activity of GC376 against different Mpros including SARS-CoV-2 Mpro mutants, but also revealed the molecular mechanism of inhibition by GC376, which will provide scientific guidance for the development of broad-spectrum drugs against SARS-CoV-2 as well as other coronaviruses.

Structural and thermodynamic properties of conserved water molecules in Mpro native: A combined approach by MD simulation and Grid Inhomogeneous Solvation Theory.

The new viral strains of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) are continuously rising, becoming more virulent, and transmissible. Therefore, the development of new antiviral drugs is essential. Due to its significant role in the viral life cycle of SARS-CoV-2, the main protease (Mpro) enzyme is a leading target for antiviral drug design. The Mpro monomer consists of domain DI, DII, and DI-DII interface. Twenty-one conserved water molecules (W4-W24) are occupied at these domains according to multiple crystal structure analyses. The crystal and MD structures reveal the presence of eight conserved water sites in domain DI, DII and remaining in the DI-DII interface. Grid-based inhomogeneous fluid solvation theory (GIST) was employed on MD structures of Mpro native to predict structural and thermodynamic properties of each conserved water site for focusing to identify the specific conserved water molecules that can easily be displaced by proposed ligands. Finally, MD water W13 is emerged as a promising candidate for water mimic drug design due to its low mean interaction energy, loose binding character with the protein, and its involvement in a water-mediated H-bond with catalytic His41 via the interaction Thr25(OG)---W13---W---His41(NE2). In this context, water occupancy, relative interaction energy, entropy, and topologies of W13 are thermodynamically acceptable for the water displacement method. Therefore, the strategic use of W13's geometrical position in the DI domain may be implemented for drug discovery against COVID disease by designing new ligands with appropriately oriented chemical groups to mimic its structural, electronic, and thermodynamic properties.

SARS-CoV-2 and oncolytic EV-D68-encoded proteases differentially regulate pyroptosis.

Pyroptosis, a pro-inflammatory programmed cell death, has been implicated in the pathogenesis of coronavirus disease 2019 and other viral diseases. Gasdermin family proteins (GSDMs), including GSDMD and GSDME, are key regulators of pyroptotic cell death. However, the mechanisms by which virus infection modulates pyroptosis remain unclear. Here, we employed a mCherry-GSDMD fluorescent reporter assay to screen for viral proteins that impede the localization and function of GSDMD in living cells. Our data indicated that the main protease NSP5 of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) blocked GSDMD-mediated pyroptosis via cleaving residues Q29 and Q193 of GSDMD. While another SARS-CoV-2 protease, NSP3, cleaved GSDME at residue G370 but activated GSDME-mediated pyroptosis. Interestingly, respiratory enterovirus EV-D68-encoded proteases 3C and 2A also exhibit similar differential regulation on the functions of GSDMs by inactivating GSDMD but initiating GSDME-mediated pyroptosis. EV-D68 infection exerted oncolytic effects on human cancer cells by inducing pyroptotic cell death. Our findings provide insights into how respiratory viruses manipulate host cell pyroptosis and suggest potential targets for antiviral therapy as well as cancer treatment.IMPORTANCEPyroptosis plays a crucial role in the pathogenesis of coronavirus disease 2019, and comprehending its function may facilitate the development of novel therapeutic strategies. This study aims to explore how viral-encoded proteases modulate pyroptosis. We investigated the impact of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and respiratory enterovirus D68 (EV-D68) proteases on host cell pyroptosis. We found that SARS-CoV-2-encoded proteases NSP5 and NSP3 inactivate gasdermin D (GSDMD) but initiate gasdermin E (GSDME)-mediated pyroptosis, respectively. We also discovered that another respiratory virus EV-D68 encodes two distinct proteases 2A and 3C that selectively trigger GSDME-mediated pyroptosis while suppressing the function of GSDMD. Based on these findings, we further noted that EV-D68 infection triggers pyroptosis and produces oncolytic effects in human carcinoma cells. Our study provides new insights into the molecular mechanisms underlying virus-modulated pyroptosis and identifies potential targets for the development of antiviral and cancer therapeutics.

SARS-CoV-2 sneezing explained.

Viral protein prods nerve cells, may be treatment target.

Reversible and irreversible inhibitors of coronavirus Nsp15 endoribonuclease.

The emergence of severe acute respiratory syndrome coronavirus 2, the causative agent of coronavirus disease 2019, has resulted in the largest pandemic in recent history. Current therapeutic strategies to mitigate this disease have focused on the development of vaccines and on drugs that inhibit the viral 3CL protease or RNA-dependent RNA polymerase enzymes. A less-explored and potentially complementary drug target is Nsp15, a uracil-specific RNA endonuclease that shields coronaviruses and other nidoviruses from mammalian innate immune defenses. Here, we perform a high-throughput screen of over 100,000 small molecules to identify Nsp15 inhibitors. We characterize the potency, mechanism, selectivity, and predicted binding mode of five lead compounds. We show that one of these, IPA-3, is an irreversible inhibitor that might act via covalent modification of Cys residues within Nsp15. Moreover, we demonstrate that three of these inhibitors (hexachlorophene, IPA-3, and CID5675221) block severe acute respiratory syndrome coronavirus 2 replication in cells at subtoxic doses. This study provides a pipeline for the identification of Nsp15 inhibitors and pinpoints lead compounds for further development against coronavirus disease 2019 and related coronavirus infections.

SARS-CoV-2 papain-like protease plays multiple roles in regulating cellular proteins in the endoplasmic reticulum.

Nsp3s are the largest nonstructural proteins of coronaviruses. These transmembrane proteins include papain-like proteases (PLpro) that play essential roles in cleaving viral polyproteins into their mature units. The PLpro of SARS-CoV viruses also have deubiquitinating and deISGylating activities. As Nsp3 is an endoplasmic reticulum (ER)-localized protein, we asked if the deubiquitinating activity of SARS-CoV-2 PLpro affects proteins that are substrates for ER-associated degradation (ERAD). Using full-length Nsp3 as well as a truncated transmembrane form we interrogated, by coexpression, three potential ERAD substrates, all of which play roles in regulating lipid biosynthesis. Transmembrane PLpro increases the level of INSIG-1 and decreases its ubiquitination. However, different effects were seen with SREBP-1 and SREBP-2. Transmembrane PLpro cleaves SREBP-1 at three sites, including two noncanonical sites in the N-terminal half of the protein, resulting in a decrease in precursors of the active transcription factor. Conversely, cleavage of SREBP-2 occurs at a single canonical site that disrupts a C-terminal degron, resulting in increased SREBP-2 levels. When this site is mutated and the degron can no longer be interrupted, SREBP-2 is still stabilized by transmembrane PLpro, which correlates with a decrease in SREBP-2 ubiquitination. All of these observations are dependent on PLpro catalytic activity. Our findings demonstrate that, when anchored to the ER membrane, SARS-CoV-2 Nsp3 PLpro can function as a deubiquitinating enzyme to stabilize ERAD substrates. Additionally, SARS-CoV-2 Nsp3 PLpro can cleave ER-resident proteins, including at sites that could escape analyses based on the established consensus sequence.

Computational study on the mechanisms of inhibition of SARS-CoV-2 Mpro by aldehyde warheads based on DFT.

SARS-CoV-2 main protease, Mpro, plays a crucial role in the virus replication cycle, making it an important target for antiviral research. In this study, a simplified model obtained through truncation is used to explore the reaction mechanism of aldehyde warhead compounds inhibiting Mpro at the level of density functional theory. According to the calculation results, proton transfer (P_T)-nucleophilic attack (N_A) is the rate-determining step in the entire reaction pathway. The water molecule that plays a catalytic role occupies the oxyanion hole, which is unfavorable for the aldehyde warhead to approach the Cys145 SH. Through a hypothetical study of substituting the main chain NH with methylene, it is further confirmed that the P_T-N_A is a proton transfer-dominated process accompanied by a nucleophilic attack reaction. In this process, the oxyanion hole serves only to stabilize the aldehyde oxygen anion and therefore does not have a significant impact on the activation free energy barrier of the step. Our research results provide a unique perspective for understanding the covalent inhibition reaction of the Mpro active site. This study also offers theoretical guidance for the design of new Mpro covalent inhibitors.

Molecular mechanisms of SARS-CoV-2 resistance to nirmatrelvir.

Nirmatrelvir is a specific antiviral drug that targets the main protease (Mpro) of SARS-CoV-2 and has been approved to treat COVID-191,2. As an RNA virus characterized by high mutation rates, whether SARS-CoV-2 will develop resistance to nirmatrelvir is a question of concern. Our previous studies have shown that several mutational pathways confer resistance to nirmatrelvir, but some result in a loss of viral replicative fitness, which is then compensated for by additional alterations3. The molecular mechanisms for this observed resistance are unknown. Here we combined biochemical and structural methods to demonstrate that alterations at the substrate-binding pocket of Mpro can allow SARS-CoV-2 to develop resistance to nirmatrelvir in two distinct ways. Comprehensive studies of the structures of 14 Mpro mutants in complex with drugs or substrate revealed that alterations at the S1 and S4 subsites substantially decreased the level of inhibitor binding, whereas alterations at the S2 and S4' subsites unexpectedly increased protease activity. Both mechanisms contributed to nirmatrelvir resistance, with the latter compensating for the loss in enzymatic activity of the former, which in turn accounted for the restoration of viral replicative fitness, as observed previously3. Such a profile was also observed for ensitrelvir, another clinically relevant Mpro inhibitor. These results shed light on the mechanisms by which SARS-CoV-2 evolves to develop resistance to the current generation of protease inhibitors and provide the basis for the design of next-generation Mpro inhibitors.

Observing inhibition of the SARS-CoV-2 helicase at single-nucleotide resolution.

The genome of SARS-CoV-2 encodes for a helicase (nsp13) that is essential for viral replication and highly conserved across related viruses, making it an attractive antiviral target. Here we use nanopore tweezers, a high-resolution single-molecule technique, to gain detailed insight into how nsp13 turns ATP-hydrolysis into directed motion along nucleic acid strands. We measured nsp13 both as it translocates along single-stranded DNA or unwinds double-stranded DNA. Our data reveal nsp13's single-nucleotide steps, translocating at ∼1000 nt/s or unwinding at ∼100 bp/s. Nanopore tweezers' high spatiotemporal resolution enables detailed kinetic analysis of nsp13 motion. As a proof-of-principle for inhibition studies, we observed nsp13's motion in the presence of the ATPase inhibitor ATPγS. We construct a detailed picture of inhibition in which ATPγS has multiple mechanisms of inhibition. The dominant mechanism of inhibition depends on the application of assisting force. This lays the groundwork for future single-molecule inhibition studies with viral helicases.

Validation of X-ray Crystal Structure Ensemble Representations of SARS-CoV-2 Main Protease by Solution NMR Residual Dipolar Couplings.

Considerable debate has focused on whether sampling of molecular dynamics trajectories restrained by crystallographic data can be used to develop realistic ensemble models for proteins in their natural, solution state. For the SARS-CoV-2 main protease, Mpro, we evaluated agreement between solution residual dipolar couplings (RDCs) and various recently reported multi-conformer and dynamic-ensemble crystallographic models. Although Phenix-derived ensemble models showed only small improvements in crystallographic Rfree, substantially improved RDC agreement over fits to a conventionally refined 1.2-Å X-ray structure was observed, in particular for residues with above average disorder in the ensemble. For a set of six lower resolution (1.55-2.19 Å) Mpro X-ray ensembles, obtained at temperatures ranging from 100 to 310 K, no significant improvement over conventional two-conformer representations was found. At the residue level, large differences in motions were observed among these ensembles, suggesting high uncertainties in the X-ray derived dynamics. Indeed, combining the six ensembles from the temperature series with the two 1.2-Å X-ray ensembles into a single 381-member "super ensemble" averaged these uncertainties and substantially improved agreement with RDCs. However, all ensembles showed excursions that were too large for the most dynamic fraction of residues. Our results suggest that further improvements to X-ray ensemble refinement are feasible, and that RDCs provide a sensitive benchmark in such endeavors. Remarkably, a weighted ensemble of 350 PDB Mpro X-ray structures provided slightly better cross-validated agreement with RDCs than any individual ensemble refinement, implying that differences in lattice confinement also limit the fit of RDCs to X-ray coordinates.

Contribution of the catalytic dyad of SARS-CoV-2 main protease to binding covalent and noncovalent inhibitors.

The effect of mutations of the catalytic dyad residues of SARS-CoV-2 main protease (MProWT) on the thermodynamics of binding of covalent inhibitors comprising nitrile [nirmatrelvir (NMV), NBH2], aldehyde (GC373), and ketone (BBH1) warheads to MPro is examined together with room temperature X-ray crystallography. When lacking the nucleophilic C145, NMV binding is ∼400-fold weaker corresponding to 3.5 kcal/mol and 13.3 °C decrease in free energy (ΔG) and thermal stability (Tm), respectively, relative to MProWT. The H41A mutation results in a 20-fold increase in the dissociation constant (Kd), and 1.7 kcal/mol and 1.4 °C decreases in ΔG and Tm, respectively. Increasing the pH from 7.2 to 8.2 enhances NMV binding to MProH41A, whereas no significant change is observed in binding to MProWT. Structures of the four inhibitor complexes with MPro1-304/C145A show that the active site geometries of the complexes are nearly identical to that of MProWT with the nucleophilic sulfur of C145 positioned to react with the nitrile or the carbonyl carbon. These results support a two-step mechanism for the formation of the covalent complex involving an initial non-covalent binding followed by a nucleophilic attack by the thiolate anion of C145 on the warhead carbon. Noncovalent inhibitor ensitrelvir (ESV) exhibits a binding affinity to MProWT that is similar to NMV but differs in its thermodynamic signature from NMV. The binding of ESV to MProC145A also results in a significant, but smaller, increase in Kd and decrease in ΔG and Tm, relative to NMV.

Kinetic analysis of RNA cleavage by coronavirus Nsp15 endonuclease: Evidence for acid-base catalysis and substrate-dependent metal ion activation.

Understanding the functional properties of severe acute respiratory syndrome coronavirus 2 nonstructural proteins is essential for defining their roles in the viral life cycle, developing improved therapeutics and diagnostics, and countering future variants. Coronavirus nonstructural protein Nsp15 is a hexameric U-specific endonuclease whose functions, substrate specificity, mechanism, and dynamics are not fully defined. Previous studies report that Nsp15 requires Mn2+ ions for optimal activity; however, the effects of divalent ions on Nsp15 reaction kinetics have not been investigated in detail. Here, we analyzed the single- and multiple-turnover kinetics for model ssRNA substrates. Our data confirm that divalent ions are dispensable for catalysis and show that Mn2+ activates Nsp15 cleavage of two different ssRNA oligonucleotide substrates but not a dinucleotide. Biphasic kinetics of ssRNA substrates demonstrates that Mn2+ stabilizes alternative enzyme states that have faster substrate cleavage on the enzyme. However, we did not detect Mn2+-induced conformational changes using CD and fluorescence spectroscopy. The pH-rate profiles in the presence and absence of Mn2+ reveal active-site ionizable groups with similar pKas of ca. 4.8 to 5.2. An Rp stereoisomer phosphorothioate modification at the scissile phosphate had minimal effect on catalysis supporting a mechanism involving an anionic transition state. However, the Sp stereoisomer is inactive because of weak binding, consistent with models that position the nonbridging phosphoryl oxygen deep in the active site. Together, these data demonstrate that Nsp15 employs a conventional acid-base catalytic mechanism passing through an anionic transition state, and that divalent ion activation is substrate dependent.

Identification of SARS-CoV-2 PLpro and 3CLpro human proteome substrates using substrate phage display coupled with protein network analysis.

Viral proteases play key roles in viral replication, and they also facilitate immune escape by proteolyzing diverse target proteins. Deep profiling of viral protease substrates in host cells is beneficial for understanding viral pathogenesis and for antiviral drug discovery. Here, we utilized substrate phage display coupled with protein network analysis to identify human proteome substrates of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) viral proteases, including papain-like protease (PLpro) and 3C-like protease (3CLpro). We first performed peptide substrates selection of PLpro and 3CLpro, and we then used the top 24 preferred substrate sequences to identify a total of 290 putative protein substrates. Protein network analysis revealed that the top clusters of PLpro and 3CLpro substrate proteins contain ubiquitin-related proteins and cadherin-related proteins, respectively. We verified that cadherin-6 and cadherin-12 are novel substrates of 3CLpro, and CD177 is a novel substrate of PLpro using in vitro cleavage assays. We thus demonstrated that substrate phage display coupled with protein network analysis is a simple and high throughput method to identify human proteome substrates of SARS-CoV-2 viral proteases for further understanding of virus-host interactions.

Conserved Characteristics of NMPylation Activities of Alpha- and Betacoronavirus NiRAN Domains.

Coronavirus genome replication and expression are mediated by the viral replication-transcription complex (RTC) which is assembled from multiple nonstructural proteins (nsp). Among these, nsp12 represents the central functional subunit. It harbors the RNA-directed RNA polymerase (RdRp) domain and contains, at its N terminus, an additional domain called NiRAN which is widely conserved in coronaviruses and other nidoviruses. In this study, we produced bacterially expressed coronavirus nsp12s to investigate and compare NiRAN-mediated NMPylation activities from representative alpha- and betacoronaviruses. We found that the four coronavirus NiRAN domains characterized to date have a number of conserved properties, including (i) robust nsp9-specific NMPylation activities that appear to operate largely independently of the C-terminal RdRp domain, (ii) nucleotide substrate preference for UTP followed by ATP and other nucleotides, (iii) dependence on divalent metal ions, with Mn2+ being preferred over Mg2+, and (iv) a key role of N-terminal residues (particularly Asn2) of nsp9 for efficient formation of a covalent phosphoramidate bond between NMP and the N-terminal amino group of nsp9. In this context, a mutational analysis confirmed the conservation and critical role of Asn2 across different subfamilies of the family Coronaviridae, as shown by studies using chimeric coronavirus nsp9 variants in which six N-terminal residues were replaced with those from other corona-, pito- and letovirus nsp9 homologs. The combined data of this and previous studies reveal a remarkable degree of conservation among coronavirus NiRAN-mediated NMPylation activities, supporting a key role of this enzymatic activity in viral RNA synthesis and processing. IMPORTANCE There is strong evidence that coronaviruses and other large nidoviruses evolved a number of unique enzymatic activities, including an additional RdRp-associated NiRAN domain, that are conserved in nidoviruses but not in most other RNA viruses. Previous studies of the NiRAN domain mainly focused on severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and suggested different functions for this domain, such as NMPylation/RNAylation of nsp9, RNA guanylyltransferase activities involved in canonical and/or unconventional RNA capping pathways, and other functions. To help resolve partly conflicting information on substrate specificities and metal ion requirements reported previously for the SARS-CoV-2 NiRAN NMPylation activity, we extended these earlier studies by characterizing representative alpha- and betacoronavirus NiRAN domains. The study revealed that key features of NiRAN-mediated NMPylation activities, such as protein and nucleotide specificity and metal ion requirements, are very well conserved among genetically divergent coronaviruses, suggesting potential avenues for future antiviral drug development targeting this essential viral enzyme.