- Open Access
Predicted coronavirus Nsp5 protease cleavage sites in the human proteome
BMC Genomic Data volume 23, Article number: 25 (2022)
The coronavirus nonstructural protein 5 (Nsp5) is a cysteine protease required for processing the viral polyprotein and is therefore crucial for viral replication. Nsp5 from several coronaviruses have also been found to cleave host proteins, disrupting molecular pathways involved in innate immunity. Nsp5 from the recently emerged SARS-CoV-2 virus interacts with and can cleave human proteins, which may be relevant to the pathogenesis of COVID-19. Based on the continuing global pandemic, and emerging understanding of coronavirus Nsp5-human protein interactions, we set out to predict what human proteins are cleaved by the coronavirus Nsp5 protease using a bioinformatics approach.
Using a previously developed neural network trained on coronavirus Nsp5 cleavage sites (NetCorona), we made predictions of Nsp5 cleavage sites in all human proteins. Structures of human proteins in the Protein Data Bank containing a predicted Nsp5 cleavage site were then examined, generating a list of 92 human proteins with a highly predicted and accessible cleavage site. Of those, 48 are expected to be found in the same cellular compartment as Nsp5. Analysis of this targeted list of proteins revealed molecular pathways susceptible to Nsp5 cleavage and therefore relevant to coronavirus infection, including pathways involved in mRNA processing, cytokine response, cytoskeleton organization, and apoptosis.
This study combines predictions of Nsp5 cleavage sites in human proteins with protein structure information and protein network analysis. We predicted cleavage sites in proteins recently shown to be cleaved in vitro by SARS-CoV-2 Nsp5, and we discuss how other potentially cleaved proteins may be relevant to coronavirus mediated immune dysregulation. The data presented here will assist in the design of more targeted experiments, to determine the role of coronavirus Nsp5 cleavage of host proteins, which is relevant to understanding the molecular pathology of coronavirus infection.
Coronaviruses are major human and livestock pathogens, and are the current focus of international attention due to an ongoing global pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). This recently emerged coronavirus likely originated in bats in China, before passing to humans in late 2019 through a secondary animal vector [1, 2]. Although 79% identical at the nucleotide level to SARS-CoV , differences in the infectious period and community spread of SARS-CoV-2 has caused a greater number of cases and deaths worldwide [3, 4]. Individuals infected by SARS-CoV-2 can develop COVID-19 disease which primarily affects the lungs, but can also cause kidney damage, coagulopathy, liver damage, and neuropathy [5,6,7,8,9,10]. Hyperinflammation, resulting from dysregulation of the immune response to SARS-CoV-2 infection, has emerged as a leading hypothesis regarding severe COVID-19 cases, which may also explain the diverse and systemic symptoms observed [11,12,13,14].
Similar to other coronaviruses, once a cell is infected, the 5′ portion of the SARS-CoV-2 (+)ssRNA genome is translated into nonstructural proteins (Nsps) required for viral replication, which are expressed covalently linked to one another (Fig. 1a) . This polyprotein must therefore be cleaved to free the individual Nsps, which is performed by two virally encoded proteases: Nsp3/papain-like protease (PLpro) and Nsp5/Main Protease (Mpro)/3C-like protease (3CLpro). Nsp5 is responsible for the majority of polyprotein cleavages and its function is conserved across coronaviruses [16, 17], making it a key drug target as its inhibition impedes viral replication (reviewed by ). Notably, the recently developed SARS-CoV-2 Nsp5 inhibitor Paxlovid reduced COVID-19 related hospital admission or death by 89% in clinical trials .
All coronavirus Nsp5 proteases identified to date are cysteine proteases in the chymotrypsin family, which primarily cleave peptides at P2-P1-P1’ residues leucine-glutamine-alanine/serine [16, 17, 20, 21], where the cleavage occurs between the P1 and P1’ residues. Nsp5 forms a homodimer for optimal catalytic function but may function as a monomer when processing its own excision from the polyprotein [22,23,24]. SARS-CoV-2 Nsp5 shares 96.1% sequence identity with SARS-CoV Nsp5 and has similar substrate specificity in vitro, but SARS-CoV-2 Nsp5 accommodates more diverse residues at substrate position P2 and may have a higher catalytic efficiency [23,24,25,26].
Coronavirus proteases also manipulate the cellular environment of infected cells to favor viral replication [27, 28], and disrupt host interferon (IFN) signaling pathways to suppress the anti-viral response of the innate immune system (reviewed by [29,30,31]). The role of coronavirus protease Nsp3 as an IFN antagonist has been well documented, including SARS-CoV-2 Nsp3 [32, 33]. Although Nsp3 proteolytic activity contributes to IFN antagonism, it is the deubiquitinating and deISGylating activities of Nsp3 that are primarily responsible [33,34,35,36,37,38,39,40]. In contrast, fewer examples of Nsp5 mediated disruption of host molecular pathways have been identified, and all are a result of its proteolytic activity [41,42,43,44,45,46,47].
Coronavirus Nsp5 antagonism of IFN is not yet clear [32, 48,49,50], but in vitro evidence supports SARS-CoV-2 Nsp5 mediated cleavage of TAB1, NLRP12, RIG-I, and RNF20 which are involved in innate immunity [51,52,53]. Hundreds of potentially cleaved peptides containing the Nsp5 consensus sequence appeared when lysate from human cells were incubated with recombinant Nsp5 from SARS-CoV, SARS-CoV-2, or hCoV-NL63, indicating a significant potential for Nsp5 mediated disruption of host proteins [45, 54]. Similarly, the abundance of potentially cleaved peptides containing a Nsp5 consensus sequence was increased in cells infected in vitro with SARS-CoV-2, which was dependent on the cell type studied . Knock down or inhibition of some of these human proteins likely cleaved by Nsp5, suppressed SARS-CoV-2 replication in vitro, suggesting that targeted host protein proteolysis is involved in viral replication . Many other SARS-CoV-2 Nsp5-host protein interactions have been predicted using proximity labeling and co-immunoprecipitation [55,56,57,58,59,60], but it is unknown if these interactions lead to Nsp5 mediated cleavage. Indeed, in vitro studies may miss Nsp5-host protein interactions due to cleavage of the host protein upon Nsp5 binding , and because individual cell types only express a limited set of human proteins. A proteome-wide prediction of coronavirus Nsp5 mediated cleavage of human proteins is therefore relevant to understanding COVID-19 pathogenesis, and how coronaviruses in general disrupt host biology.
The neural network NetCorona was previously developed in 2004, and was trained on a dataset of Nsp5 cleavage sites from seven coronaviruses including SARS-CoV . NetCorona outperforms traditional consensus motif-based approaches for identifying cleavage sites, and based on the similar specificities of SARS-CoV and SARS-CoV-2 Nsp5, we believed it could be applied to the study of SARS-CoV-2 Nsp5 interactions with human proteins. However, NetCorona only analyzes the primary amino acid sequence to predict cleavage sites, which lacks information about the 3D structure of the folded protein, and therefore how exposed a predicted cleavage site is to a protease. In particular, the solvent accessibility of a peptide motif is closely related to proteolytic susceptibility [62, 63], and in silico measurement of solvent accessibility has previously been used to help predict proteolysis [64,65,66].
In this study we used NetCorona to make predictions of Nsp5 cleavage sites across the entire human proteome, and additionally analyzed available protein structures in silico to identify highly predicted cleavage sites. We extended this analysis to examine subcellular and tissue expression patterns of the proteins predicted to be cleaved, and applied protein network analysis to identify potential key pathways disrupted by Nsp5 cleavage. Predicted Nsp5 cleavage sites in human proteins were similar to those recently identified in vitro, and human proteins predicted to be cleaved by Nsp5 were found to be involved in molecular pathways that may be relevant to the pathogenesis of COVID-19 and other coronavirus diseases.
Evaluating NetCorona performance with the SARS-CoV-2 Polyprotein
As we sought to utilize the NetCorona neural network, which had not been trained on the SARS-CoV-2 polyprotein sequence (Fig. 1a), we examined if the 11 polyprotein cleavage sites homologous to SARS-CoV would be correctly scored as cleaved (NetCorona score > 0.5). Due to the high polyprotein pp1ab sequence similarity between SARS-CoV and SARS-CoV-2, there were only three cleavage sites containing different residues (Fig. 1b). The mean NetCorona score for 10 out of the 11 SARS-CoV-2 Nsp5 cleavage sites was 0.859 (SD = 0.08), indicating highly predicted cleavages. Additional Nsp5 cleavage sites have not been identified in the SARS-CoV-2 polyprotein, and no others were predicted by NetCorona. The cleavage site at Nsp5-Nsp6 was classified as uncleaved, with a score of 0.458. SARS-CoV contains the same unique phenylalanine at position P2 of Nsp5-Nsp6, but with different P1’-P3’ residues, and received a marginal score of 0.607 in the original NetCorona paper . Phenylalanine at P2 is not found in other coronaviruses that infect humans [21, 67], nor in the other viruses used to train NetCorona, which contributed to these low scores. A P2 phenylalanine may be intentionally unfavorable at the Nsp5-Nsp6 cleavage site, to assist in its autoprocessing from the polypeptide, by limiting the ability of the cleaved peptide’s C-terminus to bind the Nsp5 active site . Interestingly, the swapped identity of the P3 and P2 residues at the SARS-CoV-2 Nsp10-Nsp11 cleavage site resulted in a higher score versus SARS-CoV (0.865 vs 0.65), due to leucine being more common at P2 versus methionine. This mutation may result in a more rapid cleavage at this site in SARS-CoV-2 versus SARS-CoV, as Nsp5 favors leucine above all other residues at P2 [17, 26].
To investigate if NetCorona can distinguish between cleaved and uncleaved motifs, NetCorona scores for all glutamine motifs in the SARS-CoV-2 pp1ab polyprotein were also determined. To gather context from the ongoing pandemic and to investigate glutamine motifs across different viral variants, 8017 SARS-CoV-2 pp1ab polyprotein sequences obtained from patient samples were scored with NetCorona (Fig. 1c, Additional file 1: Table S1). Apart from two motifs present in only 40 sequences, all glutamine motifs not naturally processed by Nsp5 received a NetCorona score < 0.5, indicating they were correctly predicted not to be cleaved. Mutations at native Nsp5 cleavage sites were also rare, with only 28 such mutated cleavage sites present in 63 sequences. Except for three mutations present in one sequence each, mutations at native Nsp5 cleavage sites were conservative and only modestly changed the NetCorona score. One sequence contained a histidine at Nsp8-Nsp9 P1 (QIO04366), resulting in NetCorona not scoring the motif. SARS-CoV and SARS-CoV-2 Nsp5 may be able to cleave motifs with histidine at P1, albeit with reduced efficiency [17, 54].
These combined results indicate that despite NetCorona not being trained on the SARS-CoV-2 sequence, it was able to correctly distinguish between cleaved versus uncleaved motifs in the pp1ab polyprotein, except for Nsp5-Nsp6. The rarity of mutated canonical cleavage sites and mutations introducing new cleavage sites (0.8 and 0.5% of sequences respectively), indicates stabilizing selection for a distinction between Nsp5 cleavage sites and all other glutamine motifs.
NetCorona predictions of Nsp5 cleavage sites in the human proteome
To generate a global view of Nsp5 cleavage sites in the human proteome, datasets were batch analyzed using NetCorona (Fig. 2). Every 9-residue motif flanking a glutamine was scored, where glutamine acts as P1 and four resides were analyzed on either side (P5-P4’). Using a NetCorona score cutoff of > 0.5, 15,057 proteins (~ 20%) in the “All Human Proteins” dataset contained a predicted cleavage site, 6056 (~ 29%) proteins in the “One Protein Per Gene”, and 2167 (~ 32%) proteins in the “Proteins With PDB” dataset (Additional file 1: Table S2-S4, raw data sets in Additional file 2, 3 and 4).
To help interpret these results, we compared the output from “One Protein Per Gene” to proteins that have been directly tested in vitro for cleavage by a coronavirus Nsp5 protease (Additional file 1: Table S5). There are 18 human proteins where cleavage sites have been mapped to the protein sequence and confirmed using an in vitro cleavage assay (CDH6, CDH20, CREB1, F2, GOLGA3, LGALS8, MAP4K5, NEMO, NLRP12, NOTCH1, OBSCN, PAICS, PNN, PTBP1, RIG-I, RNF20, RPAP1, TAB1) [41, 45,46,47, 51, 53], and also two proteins from pigs (NEMO, STAT2) [42, 44], and one from cats (NEMO) . NetCorona accurately scored 14 out of the 25 unique cleavage sites mapped in these proteins, where a glutamine was at P1. NetCorona struggled with an identical cleavage motif at Q231 in NEMO from cats, pigs, and humans, which contains an uncommon valine at P1’. Interestingly, NetCorona predicted a cleavage site in PNN at Q495, which was not identified in the original study but matches the size of a reported secondary cleavage product .
Instances where NetCorona predicted cleavages but they are not observed in vitro are also relevant to interpreting the full proteome results. NetCorona predicted cleavage sites in 22 of the 71 proteins Moustaqil et al. studied, however only TAB1 and NLRP12 were observed to be cleaved by SARS-CoV-2 Nsp5 . NetCorona predicted three cleavage sites in TAB1 and two in NLRP12, but just one predicted site in each protein matched the mapped cleavage sites.
Many other potential cleavage sites have been identified by Koudelka et al. and Pablos et al., where N-terminomics was used to identify possible cleavage sites, after cell lysate was incubated with various coronavirus Nsp5 proteases [45, 54]. Out of the 383 unique peptides identified by Koudelka et al. where a glutamine was at P1, NetCorona predicted that 167 (44%) of them would be cleaved (Additional file 1: Table S6). Similarly, out of the 155 unique peptides identified by Pablos et al. where a glutamine was at P1, NetCorona predicted that 73 (47%) of them would be cleaved (Additional file 1: Table S7). Meyer et al. also used N-terminomics to study potential Nsp5 cleavage events, following in vitro infection with SARS-CoV-2 . They identified 12 motifs in human proteins that were likely cleaved by Nsp5, of which NetCorona predicted 8 of these to be cleaved (Additional file 1: Table S8).
Several SARS-CoV-2 human protein interactomes have been made available [55,56,57,58,59,60], where predicted interactions between Nsp5 and human proteins have been reported. Interactions predicted by Samavarchi-Therani et al. were the most numerous, and the data was well annotated , which enabled a comparison to our results. These interaction scores, which varied depending on where the BioID tag was located on Nsp5 (Nsp5 C-term, N-term, or N-term on the C145A catalytically inactive mutant), were plotted against the NetCorona score from our study, which is illustrated in Additional file 5: Fig. S1 (raw data in Additional file 1: Table S9). Although statistically significant, the negative correlation between the strength of the Nsp5-human protein interaction and the maximum NetCorona score was small: ρ ranged from − 0.18 to − 0.29, r2 ranged from 0.03 to 0.08, depending on where the BioID tag was located on Nsp5. When examining only the human proteins with a positive interaction score, the mean NetCorona score ranged from 0.35 to 0.38 (SD = 0.25). Thus, Nsp5-human protein interactions predicted in vitro by Samavarchi-Therani et al. did not reflect an increased likelihood of cleavage predicted by NetCorona.
Structural characterization of predicted Nsp5 cleavage sites
We next sought to incorporate available structural information of potential protein substrates into our analysis, to address the discrepancy between the cleavage events predicted by NetCorona and mapped cleavage sites observed in vitro. The “Proteins With PDB” dataset contains only human proteins that have a solved structure available in the RCSB Protein Data Bank (PDB), however technical limitations for solving protein structures means that certain protein domains, such as transmembrane and disordered regions, may be underrepresented . To investigate if the available PDB structures contained a biased distribution of NetCorona scores, similarity between the distribution of NetCorona scores for “Proteins With PDB” and proteins in the other two datasets was assessed through the non-parametric KS test (Fig. 3a). There was insufficient evidence to reject the null hypothesis that the distribution of scores for “Proteins With PDB” proteins was equivalent to scores for “All Human Proteins” and “One Protein Per Gene” (p = 0.121 and p = 0.856, respectively), indicating that there was not significant bias in the distribution of NetCorona scores.
NetCorona scores are derived from the primary amino acid sequence, but targeted proteolysis is also dependent on the 3D structural context of the potential substrate peptide within a protein [62, 63]. Many methods have been developed to quantify this structural context in silico, and solvent accessibility has been shown to be a strong predictor of proteolysis . Accessible surface area (ASA) is commonly used to measure solvent accessibility, where a probe that approximates a water molecule is rolled around the surface of the protein, and the path traced out is the accessible surface . Thin slices are then cut through this path, to calculate the accessible surface of individual atoms. After obtaining PDB files containing motifs predicted to be cleaved by NetCorona, the total ASA of each 9 amino acid motif was calculated using Protein Structure and Interaction Analyzer (PSAIA) . This ASA was then multiplied by the motif’s NetCorona score to provide a “Nsp5 access score”, which represents both the solvent accessibility and substrate sequence preference. A Nsp5 access score was obtained for 914 glutamine motifs in 794 unique human proteins (Additional file 1: Table S10), with the process for selecting PDB files to analyze listed in Additional file 6.
Specific examples are presented to illustrate the utility of the Nsp5 access score (Fig. 3b-e). Acetylcholinesterase (ACHE) contains a motif at Q259 that was highly scored by NetCorona (0.890), but due to its presence in a tightly packed beta sheet in the core of the protein, the low ASA (38.4 Å2) results in a similarly low Nsp5 access score (34.1) and is therefore unlikely to be cleaved by Nsp5 (Fig. 3b). TGF-beta-activated kinase 1 (TAB1) is one of the few human proteins with a structure and experimental evidence of SARS-CoV-2 cleavage at specific sites (Q132 and Q444) . As illustrated in Fig. 3c, the nearby motif at Q108 was scored higher than Q132 by NetCorona, but the greater ASA of the Q132 motif contributes to a higher Nsp5 access score, which matches the experimental evidence. The human protein with the highest Nsp5 access score was DEAH box protein 15 (DHX15), as the motif surrounding Q788 was both highly scored by NetCorona and its location proximal to the C-terminus of the protein makes it highly solvent exposed (Fig. 3d).
Rationale for Nsp5 access score cut-off
To focus analysis on human proteins most likely to be cleaved by Nsp5, we determined a relevant cut-off to the Nsp5 access score. Using available structures and homology models, the Nsp5 access score of SARS-CoV-2 native cleavage sites was calculated, which ranged from 487 (Nsp15-Nsp16) to 923 (Nsp4-Nsp5) (Additional file 1: Table S11). The Nsp15-Nsp16 site (Fig. 3e) had a surprisingly low ASA (542 Å2) versus the other SARS-CoV-2 cleavage sites analyzed (mean of others 890 Å2, SD = 102 Å2), and as compared to the P5-P4’ residues of known substrates of other proteases in the chymotrypsin family (mean 678 Å2, SD = 297 Å2) (Additional file 1: Table S12).
As previously noted, NetCorona predicted cleavage sites in 22 of the 71 proteins Moustaqil et al. studied, but cleavages were only observed in vitro in two proteins . Based on available protein structures, Nsp5 access scores could be assigned to 8 unique motifs from the 22 proteins NetCorona incorrectly predicted to be cleaved, the mean of which was 332 (SD = 143). The sum of this mean and one standard deviation gives a Nsp5 access score of 475. As these were incorrectly predicted to be cleaved, this number set a lower bound for the Nsp5 access score cut-off. The score cut-off was further informed by cleavage sites recently identified by Koudelka et al. and Pablos et al. that could be assigned a Nsp5 access score (Table 1). Only a single site identified as cleaved from Moustaquil et al. (TAB1, Q132) and Yucel et al. (F2, Q494) could be assigned Nsp5 access scores, at 375 and 532 respectively.
Based on these comparisons to available experimental data, a Nsp5 access score cut-off of 500 was selected, which is further illustrated in Additional file 7: Fig. S2 (full data in Additional file 1: Table S13). This cut-off accommodates motifs with marginal NetCorona scores (~ 0.5) but maximally observed ASA (~ 1000 Å2), and the opposite scenario where a low ASA comparable to Nsp15-Nsp16 (~ 500 Å2) is matched with a high NetCorona score (~ 0.9). Ninety-two motifs in ninety-two human proteins were found to have a Nsp5 access score > 500 (Fig. 3f), which were forwarded to the next rounds of analysis.
Analysis of tissue expression and subcellular localization of predicted cleaved proteins
Proteins with a Nsp5 access score above 500 were imputed in STRING within the Cytoscape environment [71,72,73]. The STRING app computes protein network interaction by integrating information from publicly available databases, such as Reactome and Uniprot. Through textmining of the articles reported in those databases, it also compiles scores for multiple tissues and cellular compartment. The nucleus and cytosol were the top locations for human proteins with a highly predicted Nsp5 cleavage site (Fig. 4a), and the highest expression was in the nervous system and liver (Fig. 4b). The mean or summed expression score did not correlate with the Nsp5 access score (ρ = 0.03 and 0.05 respectively), nor was there a correlation between the Nsp5 access score and subcellular localization scores (ρ = − 0.08 for mean and − 0.17 for sum).
Studies of the subcellular localization of coronavirus Nsps provide insight into where Nsp5 may exist in infected cells, and thus what human proteins it may be exposed to. Flanked by transmembrane proteins Nsp4 and Nsp6 in the polyprotein, Nsp5 is exposed to the cytosol when first expressed, where it colocalizes with Nsp3 once released [74,75,76]. Recent studies have indicated that SARS-CoV-2 Nsp5 activity can be detected throughout the cytosol of a patient’s cells ex vivo , and Nsp5 is also found in the nucleus and ER [57, 77].
Through the Human Protein Atlas (HPA), we obtained information on protein expression in tissue by immunohistochemistry (IHC) together with intracellular localization obtained by confocal imaging for most of the proteins in our dataset . Proteins that are not found in the same cellular compartment as Nsp5 (nucleus, cytoplasm, endoplasmic reticulum), or where intracellular localization was unknown, were filtered out. Out of the initial 92 proteins with a Nsp5 access score over 500 and based on current knowledge, only 48 proteins were likely to be found in the same cellular compartment as Nsp5 (Fig. 5, Additional file 1: Table S14–15), indicating the greatest potential for interacting with and being cleaved by the protease. Proteins involved in apoptosis, such as CASP2, E2F1, and FNTA, had both a high Nsp5 access score and an above average expression.
Network analysis and pathways of interest
Imputation in STRING of these 48 human proteins with a Nsp5 access score over 500 and plausible colocalization, revealed multiple pathways of interest (Fig. 6, Additional file 1: Table S16). The pathway containing the most proteins that may be targeted by and colocalize with Nsp5 was mRNA processing (DHX15, ELAVL1, LTV1, PABPC3, RPL10, RPUSD1, SKIV2L2, SMG7, TDRD7). Another prominent pathway was apoptosis, with multiple proteins involved directly in apoptosis or its regulation (CASP2, E2F1, FNTA, MAPT, PTPN13). DNA damage response, mediated through ATF2, NEIL1, PARP2, and RAD50 may also be targeted by Nsp5. PARP2 had the second highest Nsp5 access score in our analysis, and the predicted cleavage site at Q352 is located between the DNA-binding domain and the catalytic domain .
Proteins involved in membrane trafficking (RAB27B and SNX10), or in microtubule organization (DNM1, HTT, MAPRE3, TSC1) were also enriched in this focused dataset, which were grouped together under the descriptor “vesicle trafficking”. Two proteins related to ubiquitination (UBA1 and USP4) were also amongst these potential Nsp5 targets. Finally, a group of proteins implicated in cytokine response was also strongly predicted to be cleaved (AIMP1, MAPK12, and PTPN2), which are involved in downstream signaling of multiple cytokines [80,81,82,83].
To provide context to the growing list of coronavirus-host protein-protein interactions, and to aid in the interpretation of experiments focused on human proteins cleaved by coronavirus Nsp5, we applied a bioinformatics approach to predict human proteins cleaved by Nsp5. Our proteome-wide investigation complements in vitro experiments, which are limited to only a subset of potential human protein substrates based on what proteins are expressed by the cell type chosen, resulting in different proteins appearing to be cleaved by [46, 54], or interact with Nsp5 [55,56,57,58,59,60].
The NetCorona neural network generated long lists of potentially cleaved human proteins, but mismatches between these predictions and the in vitro mapping of Nsp5 cleavage sites indicated that NetCorona scores alone were insufficient for accurate predictions. We added to these NetCorona predictions, which are based on primary sequence alone, by calculating solvent accessibility of the predicted cleaved motifs, which is closely related to proteolytic susceptibility [62, 63]. We focused this analysis to high quality protein structures, and avoided homology models and predicted structures, to connect our predictions to real protein structures. This was made possible thanks to the PSAIA tool which automated the measurement of motif solvent accessibility with an easy-to-use GUI that handled batch input of PDB files .
Human proteins predicted to be cleaved by Nsp5 did not correlate with Nsp5-human protein-protein interactions predicted in vitro, and Nsp5 overall appears to interact with fewer human proteins compared to other Nsps and structural proteins . This may be because the proteolytic activity of Nsp5 reduces the efficiency of proximity labeling/affinity purification, whereby Nsp5 may cleave proteins it interacts with most favorably, reducing the appearance of host protein interactions. The small but statistically significant negative correlation between the strength of the Nsp5-human protein interaction and the human protein’s maximum NetCorona score may be evidence of this. Indeed, different sets of interacting proteins are obtained when using the catalytically inactive Nsp5 mutant C145A versus the wildtype Nsp5 [55, 57, 60]. These protein-protein interaction studies also rely on the overexpression of viral proteins in a non-native context. We therefore hypothesize that the interactions observed by proximity labeling/affinity purification do not reflect Nsp5 mediated proteolysis and instead represent non-proteolytic protein-protein interactions, which may still be important to understanding Nsp5’s role in modulating host protein networks.
N-terminomics based approaches have identified many potential Nsp5 cleavage sites in human proteins [45, 46, 54], but they have some limitations that a bioinformatics approach can complement. Trypsin is used in the preparation of samples for mass spectrometry, which generates cleavages at lysine and arginine residues that are not N-terminal to a proline. Lysine and arginine appear in many cleavage sites predicted by NetCorona, meaning that cleavage by trypsin may mask true cleavage sites by artificially generating a N-terminus proximal to a P1 glutamine residue. Only 38 cleavage sites were commonly identified by both Koudelka et al. and Pablos et al. using similar N-terminomics approaches, out of the hundreds of potentially cleaved peptides that each study identified [45, 54], likely as these studies used different cell lines and thus different proteins will be expressed. Meyer et al. point out that the lysate-based method used by Koudelka et al. and Pablos et al. strips proteins of their subcellular context, which may lead to observed cleavage events that are not possible in vivo during infection . Even so, the SARS-CoV-2 cellular infection-based method Meyer et al. used, paired with N-terminomics, resulted in cell-type dependent differences . Recently, Yucel et al. presented evidence that proteins predicted to be cleaved by SARS-CoV-2 Nsp5 are rapidly degraded by the proteosome, meaning that N-terminomics approaches may miss potentially cleaved proteins due to fragments being cleared.
In contrast, our bioinformatics analysis is cell-type and methodology agnostic as it examined the entire human proteome, and it is not affected by the potential for the proteosome to degrade the fragments cleaved by Nsp5. The cleavage sites predicted in silico, combined with analysis of solved protein structures, and knowledge of Nsp5 subcellular localization and protein networks, identified several interesting human cellular pathways and proteins.
mRNA processing, vesicle trafficking, cytokine signaling, ubiquitination, DNA damage recognition, and apoptosis were the cellular pathways found to be most affected by the predicted Nsp5 cleavage of human proteins. mRNA processing was also predicted by Pablos et al. to be affected by SARS-CoV-2 Nsp5 , and coronaviruses are known to manipulate mRNA processing both for suppressing antiviral responses , and for increasing viral takeover of the cell by degrading host mRNA [15, 85].
Vesicle trafficking is known to be manipulated by coronaviruses to favor the formation of complex double membrane vesicles during viral replication . RAB27B was predicted to be cleaved by Nsp5, and has roles in the trafficking of vesicles from the trans-Golgi network , which coronavirus proteins accumulate in . RAB27B also plays a role in the exocytosis of neutrophil granules , suggesting that cleavage of RAB27B could also disrupt the immune response to infection.
AIMP1, MAPK12, PTPN2 have roles in cytokine response and antiviral defense, and were identified as highly predicted targets of Nsp5 cleavage. AIMP1 is crucial to antiviral defense in mice , and it is targeted for degradation by hepatitis C virus envelope protein E2 . MAPK12 function appears to be important for SARS-CoV-2 replication, as knockdown by siRNA results in lower virus titers . This phenotype is similar to the knockdown of Nsp5 targeted proteins resulting in lower virus replication, as observed by Meyer et al. . T-cell protein tyrosine phosphatase, PTPN2, negatively regulates the antiviral response of MITA  and the JAK-STAT pathway of the innate immune system . Knockout of PTPN2 results in systemic inflammatory responses in mice resulting in premature death , and a genetic polymorphism resulting in PTPN2 loss of function increases ACE2 expression, resulting in greater susceptibility to SARS-CoV-2 infection .
DHX15 contained a predicted cleavage site with the highest Nsp5 access score, and the protein may co-localize with Nsp5 in the nuclei of infected cells, making it a significant protein of interest. DHX15 is a DExD/H-box helicase, a family of proteins that serves to detect foreign RNA, triggering an antiviral response (reviewed by ). The role of DHX15 in anti-viral defense is diverse, including by binding viral RNA with NLRP6, which activates Type I/III interferons and IFN-stimulated genes in the intestine of mice . DHX15 mediated sensing of viral RNA activates MAPK and NK-κB innate immune signaling , and also acts as a coreceptor of viral RNA with RIG-I that increases antiviral response and cytokine production . Interestingly, RIG-I was recently identified as being cleaved by SARS-CoV-2 Nsp5 in vitro .
DHX15 is not predicted to be a significant interactor with SARS-CoV-2 proteins , however it is capable of binding both dsRNA and ssRNA , suggesting it could bind coronavirus ssRNA. The location of the predicted Nsp5 cleavage site (Q788) is very close to its C-terminus (Y795), so cleavage by Nsp5 would remove only seven amino acids. However, there is a SUMOylation site (K786) at P3 of the cleavage motif, and the de-SUMOylation of DHX15 results in increased antiviral signaling . Therefore, it is possible that Nsp5 cleavage at Q788 may disrupt SUMOylation, by reducing the length of the peptide accessible by the SUMOylation complex, contributing to the general dysregulation of a coordinated innate immune response to viral infection.
Nsp3 mediated modulation of ubiquitination has been shown to be important for IFN antagonism [33,34,35,36,37,38,39,40], and there is also evidence for Nsp5 mediated reduction of ubiquitination [49, 54]. USP4, which was predicted to be cleaved by Nsp5, has a stabilizing effect on RIG-I, resulting in an elevated innate immune response . Thus, the predicted cleavage of USP4 would attenuate this response. UBA1 catalyzes the first step in ubiquitin conjugation, thus its potential cleavage by Nsp5 would have a myriad of disruptive effects in an infected cell.
PARP2 contained a predicted Nsp5 cleavage site with the second highest Nsp5 access score, at Q352. PARP2 is involved in DNA damage recognition and repair, and its correct functioning prevents apoptosis in the event of a double stranded break (DSB) . PARP2 also plays a role in the adaptive immune system, as it helps prevent the accumulation of DSBs during TCRα gene rearrangements in thymocytes, promoting T-cell maturation . The predicted Nsp5 site occurs between the DNA binding and catalytic domains of PARP2 , meaning the cleaved protein would be unable to recognize damaged DNA, contributing to apoptosis. Interestingly, this is similar to the native activity of human caspase-8, which cleaves PARP2 in its DNA binding domain during apoptosis . SARS-CoV-2 infection of lung epithelial cells was found to increase caspase-8 activity, resulting in cleaved PARP1, a homolog of PARP2 . SARS-CoV Nsp5 activity is known to be pro-apoptotic, via the activation of caspase-3 and caspase-9 . Overall, if coronavirus Nsp5 cleaves PARP2, it could contribute to the pro-apoptotic cell state observed in infected cells.
It is therefore possible that Nsp5 mediated cleavage of human proteins manipulates a variety of host cellular pathways, which may bias infected cells to favor viral replication while also disrupting antiviral responses, contributing to the molecular pathogenesis of diseases caused by coronaviruses. These pathways and specific proteins mentioned above represent interesting targets for further analysis in vitro.
The predictions presented in this study are not specific to SARS-CoV-2 Nsp5, but rather are predictions of the Nsp5 of coronaviruses in general, due to their shared mechanism and similar substrate specificities. Recent studies comparing Nsp5 proteases from various coronaviruses have indicated that despite sharing significant sequence and structural similarity, they cleave and interact with different human proteins in vitro [54, 60]. In general, a pan-coronavirus predictor of Nsp5 cleavage sites may not be feasible. For example, SARS-CoV-2 Nsp5 accommodates more diversity at P2 than SARS-CoV Nsp5 , which would influence the human proteins that could be cleaved. Refinements to the NetCorona neural network to improve its predictive accuracy, or make them virus-specific, would be beneficial and have recently been attempted [47, 104]. Indeed, the significant amount of research into SARS-CoV-2 Nsp5 mediated cleavage of human proteins has provided large datasets that could be used to create a neural network specific to SARS-CoV-2 Nsp5. Nsp5 also appears to accommodate more diversity at P1’ and P3 than what NetCorona was trained on, based on the naturally occurring SARS-CoV-2 pp1ab variants we report here. That histidine and methionine can be accommodated at P1, and phenylalanine at P2, albeit unfavorably, further adds complexity to what human proteins may be cleaved by Nsp5 [17, 26, 45, 54, 67]. The results of this study are instead meant to provide guidance for in vitro experimental design and interpretation of experimental results, in addition to suggesting proteome-wide trends in molecular pathways that Nsp5 from coronaviruses may disrupt to evade the immune system, or augment to favor viral replication through targeted cleavage of human proteins.
The large volume of recent coronavirus research and data requires proteome-wide views of interpretation. The results of this study are intended to complement the various in vitro approaches that have been used to identify coronavirus Nsp5-human protein interactions, and to map specific Nsp5 cleavage sites in human proteins. As Meyer et al. discuss, specific targeting of proteins by Nsp5 appears likely, as the knockdown of certain Nsp5-targeted proteins reduces viral reproduction . We have built upon the original NetCorona study by performing detailed structural analysis of predicted cleavage sites, and protein network and pathway analysis. Nsp5 was predicted to play a role in the targeted disruption or augmentation of several key molecular pathways that relate to recent studies of SARS-CoV-2, and thus are interesting targets for future analysis and characterization both in the context of the COVID-19 pandemic and to improve our understanding of coronaviruses overall. We hope that our analysis and the proteome-wide datasets generated will aid in the interpretation and design of additional experiments towards understanding Nsp5’s role in coronavirus molecular pathology. Specifically, the 48 proteins highly predicted to be cleaved, and also expected to be found in the same cellular compartment as Nsp5, should be investigated in vitro for Nsp5 mediated cleavage.
Amino acid sequence datasets
Three datasets of amino acid sequences of human proteins were downloaded from the UniProt reference human proteome on April 14 2020: “All Human Proteins” contains 74,811 human protein sequences including splice variants and predicted proteins; “One Protein Per Gene” contains 20,595 human protein sequences where only one sequence is provided for each gene; “Proteins With PDB” is the “All Human Proteins” dataset filtered in UniProt for associated entries in the RCSB Protein Data Bank, generating a dataset of 6806 human protein sequences. 9404 amino acid sequences of the SARS-CoV-2 pp1ab polyprotein were obtained from NCBI Virus on August 4 2020. Sequences with inconclusive “X” residues were filtered out as they were not correctly handled by NetCorona, leaving 8017 SARS-CoV-2 pp1ab polyprotein sequences to be analyzed.
The command line version of NetCorona was used to predict Nsp5 cleavage site scores for human and viral protein sequences . To overcome the input file limit of 50,000 amino acids per submission and handle sequences with non-standard amino acids, a Python script was developed. This script partitions the input data, runs the NetCorona neural network on each subset, and parses and concatenates the output data. The output file includes sequence accession number, position of P1-Glutamine (Q) residue, Netcorona score (0.000–1.000) and a 10 amino acid sequence motif of positions P5-P5’ (Additional file 2, 3 and 4). Note that the neural network itself uses positions P5-P4’ (9 residues) for calculating the score. Gene names and other identifiers associated with each UniProt ID containing a NetCorona score > 0.5 were collected in Microsoft Excel (Additional file 1: Table S2-S4). Scores from NetCorona run on each dataset of proteins were parsed and compared with a Kolmogorov-Smirnov (KS) test, to assess the null hypothesis that the scores for each population are drawn from the same distribution. Unique glutamine motifs from pp1ab polyprotein sequences were identified using Microsoft Excel (Additional file 1: Table S1). Statistical analysis and the generation of graphs was performed using GraphPad Prism (version 9.1.0).
PDB metadata associated with proteins in the “Proteins With PDB” dataset that also contained a predicted Nsp5 cleavage (NetCorona score > 0.5), were downloaded from the RCSB PDB website by generating a custom report in .csv format. Homology models, and structures with a resolution greater than 8 Å or where resolution was not reported, were removed. Nsp5 cleavage sites predicted by NetCorona were matched with one PDB file per cleavage site, by searching the PDB metadata for the predicted 9 amino acid cleavage motif using Microsoft Excel (Additional file 6). The entire predicted 9 amino acid motif must appear in the PDB file to be considered a match. Matches between a PDB file and predicted cleavage motif were manually corrected when the motif sequence appeared by chance in a PDB containing the incorrect protein.
PDB files containing a predicted Nsp5 cleavage site were then batch downloaded from the RCSB PDB, and analyzed 100 at a time using the Protein Structure and Interaction Analyzer (PSAIA) tool using default settings , with chains in each PDB analyzed independently. The total accessible surface area (ASA) of each residue was calculated using a Z slice of 0.25 Å and a probe radius of 1.4 Å. XML files output by PSAIA were combined in Microsoft Excel, to create searchable datasets for each 9 amino acid motif predicted to be cleaved by NetCorona, and the total ASA of all atoms in each 9 amino acid motif were summed. The motif’s ASA was then multiplied by the NetCorona score to provide a Nsp5 access score.
Proteins known to be cleaved by mammalian chymotrypsin-like proteases were independently obtained from the RCSB PDB, and the known cleaved motifs were analyzed as above. Protein structures and homology models of SARS-CoV-2 proteins were obtained the RCSB PDB and from SWISS-MODEL  and were analyzed as above. Publication quality figures were generated using PyMOL 2.3.0.
Tissue expression and subcellular localization analysis
Proteins with Nsp5 access score above 500 were loaded into the STRING app  within Cytoscape  (version 1.6.0 and 3.8.2 respectively) using Uniprot ID, Homo sapiens background, 0.80 confidence score cut-off and no additional interactor for the pathway enrichment analysis (Additional file 1: Table S14-S16). The node table (including tissue expression scores and compartments score for each protein) was exported to R for wrangling and data visualization using the tidyverse and ggrepel packages [106,107,108].
To increase confidence, tissue expression and subcellular localization data were obtained from the Human Protein Atlas which are all based on immunohistochemistry (tissue expression) or confocal microscopy (subcellular localization) [78, 109]. Each entry was then matched in R, table joining was done using Uniprot IDs. Expression levels noted as “Not detected”, “Low”, “Medium” or “High” were replaced by numeric values ranging from 0 to 3. Mean expression was calculated as the mean expression across all tissues, removing missing values from the analysis.
The following intracellular locations were used to encompass the nucleus, cytoplasm, and endoplasmic reticulum: “Cytosol”, “Nucleoplasm”, “Endoplasmic reticulum”, “Microtubules”, “Nuclear speckles”, “Intermediate filaments”, “Nucleoli”, “Nuclear bodies”. All the proteins that did not include one or more of these locations in the HPA database were excluded from further analysis.
Protein network analysis
The 48 proteins with a Nsp5 access score > 500 and that had the potential to be found in the same cellular compartment as Nsp5 were imported into the STRING app (again within Cytoscape) while allowing a maximum of 5 additional interactor for the network generation instead of none. All the other parameters were left unchanged. Individual nodes that had no protein-protein interactions with other proteins in the network were manually moved closer to other nodes presenting the same or similar pathway. When proteins could interact in multiple pathways represented here, a “main pathway” was assigned based on literature search. Node color was a gradient based on Nsp5 access score. Node size increased with the mean expression. Edges represent protein-protein interaction (confidence > 0.80). Gene name labels were colored based on the Nsp5 access acore for readability only.
Availability of data and materials
The datasets generated and/or analyzed during the current study are available in the Figshare repository, https://doi.org/10.6084/m9.figshare.19158689
PSAIA was obtained here: http://complex.zesoi.fer.hr/index.php/en/10-category-en-gb/tools-en/19-psaia-en
Andersen KG, Rambaut A, Lipkin WI, Holmes EC, Garry RF. The proximal origin of SARS-CoV-2. Nat Med. 2020;26(4):450–2.
Zhou H, Chen X, Hu T, Li J, Song H, Liu Y, et al. A novel bat coronavirus closely related to SARS-CoV-2 contains natural insertions at the S1/S2 cleavage site of the spike protein. Curr Biol. 2020;30(11):2196–203 e3.
Wilder-Smith A, Chiew CJ, Lee VJ. Can we contain the COVID-19 outbreak with the same measures as for SARS? Lancet Infect Dis. 2020;20(5):e102–e7.
Petersen E, Koopmans M, Go U, Hamer DH, Petrosillo N, Castelli F, et al. Comparing SARS-CoV-2 with SARS-CoV and influenza pandemics. Lancet Infect Dis. 2020;20(9):e238–e44.
Yang X, Yu Y, Xu J, Shu H, Xia J, Liu H, et al. Clinical course and outcomes of critically ill patients with SARS-CoV-2 pneumonia in Wuhan, China: a single-centered, retrospective, observational study. Lancet Respir Med. 2020;8(5):475–81.
Robba C, Battaglini D, Pelosi P, Rocco PRM. Multiple organ dysfunction in SARS-CoV-2: MODS-CoV-2. Expert rev. Respir Med. 2020;14(9):865-8.
Puelles VG, Lutgehetmann M, Lindenmeyer MT, Sperhake JP, Wong MN, Allweiss L, et al. Multiorgan and renal tropism of SARS-CoV-2. N Engl J Med. 2020.
Helms J, Tacquard C, Severac F, Leonard-Lorant I, Ohana M, Delabranche X, et al. High risk of thrombosis in patients with severe SARS-CoV-2 infection: a multicenter prospective cohort study. Intensive Care Med. 2020;46(6):1089–98.
Xu L, Liu J, Lu M, Yang D, Zheng X. Liver injury during highly pathogenic human coronavirus infections. Liver Int. 2020;40(5):998–1004.
Whittaker A, Anson M, Harky A. Neurological manifestations of COVID-19: a systematic review and current update. Acta Neurol Scand. 2020;142(1):14–22.
Giamarellos-Bourboulis EJ, Netea MG, Rovina N, Akinosoglou K, Antoniadou A, Antonakos N, et al. Complex immune Dysregulation in COVID-19 patients with severe respiratory failure. Cell Host Microbe. 2020;27(6):992–1000 e3.
Lucas C, Wong P, Klein J, Castro TBR, Silva J, Sundaram M, et al. Longitudinal analyses reveal immunological misfiring in severe COVID-19. Nature. 2020;584(7821):463–9.
Kuri-Cervantes L, Pampena MB, Meng W, Rosenfeld AM, Ittner CAG, Weisman AR, et al. Comprehensive mapping of immune perturbations associated with severe COVID-19. Sci Immunol. 2020;5(49):eabd7114.
Mathew D, Giles JR, Baxter AE, Oldridge DA, Greenplate AR, Wu JE, et al. Deep immune profiling of COVID-19 patients reveals distinct immunotypes with therapeutic implications. Science. 2020;369(6508):eabc8511.
Hartenian E, Nandakumar D, Lari A, Ly M, Tucker JM, Glaunsinger BA. The molecular virology of coronaviruses. J Biol Chem. 2020;295(37):12910–34.
Ziebuhr J, Snijder EJ, Gorbalenya AE. Virus-encoded proteinases and proteolytic processing in the Nidovirales. J Gen Virol. 2000;81(Pt 4):853–79.
Chuck CP, Chow HF, Wan DC, Wong KB. Profiling of substrate specificities of 3C-like proteases from group 1, 2a, 2b, and 3 coronaviruses. PLoS One. 2011;6(11):e27228.
Ullrich S, Nitsche C. The SARS-CoV-2 main protease as drug target. Bioorg Med Chem Lett. 2020;30(17):127377.
Mahase E. Covid-19: Pfizer’s paxlovid is 89% effective in patients at risk of serious illness, company reports. BMJ. 2021;375:n2713.
Hegyi A, Ziebuhr J. Conservation of substrate specificities among coronavirus main proteases. J Gen Virol. 2002;83(Pt 3):595–9.
Tomar S, Johnston ML, St John SE, Osswald HL, Nyalapatla PR, Paul LN, et al. Ligand-induced dimerization of Middle East respiratory syndrome (MERS) coronavirus nsp5 protease (3CLpro): implications for nsp5 regulation and the development of antivirals. J Biol Chem. 2015;290(32):19403–22.
Muramatsu T, Kim YT, Nishii W, Terada T, Shirouzu M, Yokoyama S. Autoprocessing mechanism of severe acute respiratory syndrome coronavirus 3C-like protease (SARS-CoV 3CLpro) from its polyproteins. FEBS J. 2013;280(9):2002–13.
Zhang L, Lin D, Sun X, Curth U, Drosten C, Sauerhering L, et al. Crystal structure of SARS-CoV-2 main protease provides a basis for design of improved alpha-ketoamide inhibitors. Science. 2020;368(6489):409–12.
Jin Z, Du X, Xu Y, Deng Y, Liu M, Zhao Y, et al. Structure of M(pro) from SARS-CoV-2 and discovery of its inhibitors. Nature. 2020;582(7811):289–93.
Yoshimoto FK. The proteins of severe acute respiratory syndrome Coronavirus-2 (SARS CoV-2 or n-COV19), the cause of COVID-19. Protein J. 2020;39(3):198–216.
Rut W, Groborz K, Zhang L, Sun X, Zmudzinski M, Pawlik B, et al. SARS-CoV-2 M(pro) inhibitors and activity-based probes for patient-sample imaging. Nat Chem Biol. 2021;17(2):222–8.
Angelini MM, Akhlaghpour M, Neuman BW, Buchmeier MJ. Severe acute respiratory syndrome coronavirus nonstructural proteins 3, 4, and 6 induce double-membrane vesicles. mBio. 2013;4(4):e00524–13.
Wolff G, Limpens R, Zevenhoven-Dobbe JC, Laugks U, Zheng S, de Jong AWM, et al. A molecular pore spans the double membrane of the coronavirus replication organelle. Science. 2020;369(6509):1395–8.
Totura AL, Baric RS. SARS coronavirus pathogenesis: host innate immune responses and viral antagonism of interferon. Curr Opin Virol. 2012;2(3):264–75.
Lei J, Hilgenfeld R. RNA-virus proteases counteracting host innate immunity. FEBS Lett. 2017;591(20):3190–210.
de Wilde AH, Snijder EJ, Kikkert M, van Hemert MJ. Host factors in coronavirus replication. Curr Top Microbiol Immunol. 2018;419:1–42.
Lei X, Dong X, Ma R, Wang W, Xiao X, Tian Z, et al. Activation and evasion of type I interferon responses by SARS-CoV-2. Nat Commun. 2020;11(1):3810.
Shin D, Mukherjee R, Grewe D, Bojkova D, Baek K, Bhattacharya A, et al. Papain-like protease regulates SARS-CoV-2 viral spread and innate immunity. Nature. 2020;587(7835):657–62.
Frieman M, Ratia K, Johnston RE, Mesecar AD, Baric RS. Severe acute respiratory syndrome coronavirus papain-like protease ubiquitin-like domain and catalytic domain regulate antagonism of IRF3 and NF-kappaB signaling. J Virol. 2009;83(13):6689–705.
Yang X, Chen X, Bian G, Tu J, Xing Y, Wang Y, et al. Proteolytic processing, deubiquitinase and interferon antagonist activities of Middle East respiratory syndrome coronavirus papain-like protease. J Gen Virol. 2014;95(Pt 3):614–26.
Mielech AM, Kilianski A, Baez-Santos YM, Mesecar AD, Baker SC. MERS-CoV papain-like protease has deISGylating and deubiquitinating activities. Virology. 2014;450-451:64–70.
Chen X, Yang X, Zheng Y, Yang Y, Xing Y, Chen Z. SARS coronavirus papain-like protease inhibits the type I interferon signaling pathway through interaction with the STING-TRAF3-TBK1 complex. Protein Cell. 2014;5(5):369–81.
Li SW, Wang CY, Jou YJ, Huang SH, Hsiao LH, Wan L, et al. SARS coronavirus papain-like protease inhibits the TLR7 signaling pathway through removing Lys63-linked Polyubiquitination of TRAF3 and TRAF6. Int J Mol Sci. 2016;17(5):678.
Knaap RCM, Fernández-Delgado R, Dalebout TJ, Oreshkova N, Bredenbeek PJ, Enjuanes L, et al. The deubiquitinating activity of Middle East respiratory syndrome coronavirus papain-like protease delays the innate immune response and enhances virulence in a mouse model. bioRxiv. 2019:751578.
Freitas BT, Durie IA, Murray J, Longo JE, Miller HC, Crich D, et al. Characterization and noncovalent inhibition of the Deubiquitinase and deISGylase activity of SARS-CoV-2 papain-like protease. ACS Infect Dis. 2020;6(8):2099–109.
Wang D, Fang L, Shi Y, Zhang H, Gao L, Peng G, et al. Porcine epidemic diarrhea virus 3C-like protease regulates its interferon antagonism by cleaving NEMO. J Virol. 2016;90(4):2090–101.
Zhu X, Fang L, Wang D, Yang Y, Chen J, Ye X, et al. Porcine deltacoronavirus nsp5 inhibits interferon-beta production through the cleavage of NEMO. Virology. 2017;502:33–8.
Chen S, Tian J, Li Z, Kang H, Zhang J, Huang J, et al. Feline infectious peritonitis virus Nsp5 inhibits type I interferon production by cleaving NEMO at multiple sites. Viruses. 2019;12(1):43.
Zhu X, Wang D, Zhou J, Pan T, Chen J, Yang Y, et al. Porcine Deltacoronavirus nsp5 antagonizes type I interferon signaling by cleaving STAT2. J Virol. 2017;91(10):e00003–17.
Pablos I, Machado Y, de Jesus HCR, Mohamud Y, Kappelhoff R, Lindskog C, et al. Mechanistic insights into COVID-19 by global analysis of the SARS-CoV-2 3CL(pro) substrate degradome. Cell Rep. 2021;37(4):109892.
Meyer B, Chiaravalli J, Gellenoncourt S, Brownridge P, Bryne DP, Daly LA, et al. Characterising proteolysis during SARS-CoV-2 infection identifies viral cleavage sites and cellular targets with therapeutic potential. Nat Commun. 2021;12(1):5553.
Yucel N, Marchiano S, Tchelepi E, Paterlini G, McAfee Q, Nimmagadda N, et al. Prediction and validation of host cleavage targets of SARS-CoV-2 3C-like protease. bioRxiv. 2022; 2022.01.17.476677.
Xia H, Cao Z, Xie X, Zhang X, Chen JY, Wang H, et al. Evasion of type I interferon by SARS-CoV-2. Cell Rep. 2020;33:108234.
Wu Y, Ma L, Zhuang Z, Cai S, Zhao Z, Zhou L, et al. Main protease of SARS-CoV-2 serves as a bifunctional molecule in restricting type I interferon antiviral signaling. Signal Transduct Target Ther. 2020;5(1):221.
Zheng Y, Deng J, Han L, Zhuang MW, Xu Y, Zhang J, et al. SARS-CoV-2 NSP5 and N protein counteract the RIG-I signaling pathway by suppressing the formation of stress granules. Signal Transduct Target Ther. 2022;7(1):22.
Moustaqil M, Ollivier E, Chiu HP, Van Tol S, Rudolffi-Soto P, Stevens C, et al. SARS-CoV-2 proteases PLpro and 3CLpro cleave IRF3 and critical modulators of inflammatory pathways (NLRP12 and TAB1): implications for disease presentation across species. Emerg Microbes Infect. 2021;10(1):178–95.
Liu Y, Qin C, Rao Y, Ngo C, Feng JJ, Zhao J, et al. SARS-CoV-2 Nsp5 demonstrates two distinct mechanisms targeting RIG-I and MAVS to evade the innate immune response. mBio. 2021;12(5):e0233521.
Zhang S, Wang J, Cheng G. Protease cleavage of RNF20 facilitates coronavirus replication via stabilization of SREBP1. Proc Natl Acad Sci U S A. 2021;118(37).
Koudelka T, Boger J, Henkel A, Schonherr R, Krantz S, Fuchs S, et al. N-Terminomics for the identification of in vitro substrates and cleavage site specificity of the SARS-CoV-2 Main protease. Proteomics. 2021;21(2):e2000246.
Gordon DE, Jang GM, Bouhaddou M, Xu J, Obernier K, White KM, et al. A SARS-CoV-2 protein interaction map reveals targets for drug repurposing. Nature. 2020;583(7816):459–68.
Bojkova D, Klann K, Koch B, Widera M, Krause D, Ciesek S, et al. Proteomics of SARS-CoV-2-infected host cells reveals therapy targets. Nature. 2020;583(7816):469–72.
Samavarchi-Tehrani P, Abdouni H, Knight JDR, Astori A, Samson R, Lin Z-Y, et al. A SARS-CoV-2 – host proximity interactome. bioRxiv. 2020; 2020.09.03.282103.
Laurent EMN, Sofianatos Y, Komarova A, Gimeno J-P, Tehrani PS, Kim D-K, et al. Global BioID-based SARS-CoV-2 proteins proximal interactome unveils novel ties between viral polypeptides and host factors involved in multiple COVID19-associated mechanisms. bioRxiv. 2020; 2020.08.28.272955.
Li J, Guo M, Tian X, Wang X, Yang X, Wu P, et al. Virus-host Interactome and proteomic survey reveal potential virulence factors influencing SARS-CoV-2 pathogenesis. Med (N Y). 2021;2(1):99–112 e7.
Gordon DE, Hiatt J, Bouhaddou M, Rezelj VV, Ulferts S, Braberg H, et al. Comparative host-coronavirus protein interaction networks reveal pan-viral disease mechanisms. Science. 2020;370(6521):eabe9403.
Kiemer L, Lund O, Brunak S, Blom N. Coronavirus 3CLpro proteinase cleavage sites: possible relevance to SARS virus pathology. BMC Bioinformatics. 2004;5:72.
Hubbard SJ, Campbell SF, Thornton JM. Molecular recognition. Conformational analysis of limited proteolytic sites and serine proteinase protein inhibitors. J Mol Biol. 1991;220(2):507–30.
Kazanov MD, Igarashi Y, Eroshkin AM, Cieplak P, Ratnikov B, Zhang Y, et al. Structural determinants of limited proteolysis. J Proteome Res. 2011;10(8):3642–51.
Boyd SE, Pike RN, Rudy GB, Whisstock JC, Garcia de la Banda M. PoPS: a computational tool for modeling and predicting protease specificity. J Bioinforma Comput Biol. 2005;3(3):551–85.
Verspurten J, Gevaert K, Declercq W, Vandenabeele P. SitePredicting the cleavage of proteinase substrates. Trends Biochem Sci. 2009;34(7):319–23.
Song J, Tan H, Perry AJ, Akutsu T, Webb GI, Whisstock JC, et al. PROSPER: an integrated feature-based tool for predicting protease substrate cleavage sites. PLoS One. 2012;7(11):e50300.
Muramatsu T, Takemoto C, Kim YT, Wang H, Nishii W, Terada T, et al. SARS-CoV 3CL protease cleaves its C-terminal autoprocessing site by novel subsite cooperativity. Proc Natl Acad Sci U S A. 2016;113(46):12997–3002.
Peng K, Obradovic Z, Vucetic S. Exploring bias in the protein data Bank using contrast classifiers. Pac Symp Biocomput. 2004:435–46.
Lee B, Richards FM. The interpretation of protein structures: estimation of static accessibility. J Mol Biol. 1971;55(3):379–400.
Mihel J, Sikic M, Tomic S, Jeren B, Vlahovicek K. PSAIA - protein structure and interaction analyzer. BMC Struct Biol. 2008;8:21.
Doncheva NT, Morris JH, Gorodkin J, Jensen LJ. Cytoscape StringApp: network analysis and visualization of proteomics data. J Proteome Res. 2019;18(2):623–32.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504.
Su G, Morris JH, Demchak B, Bader GD. Biological network exploration with Cytoscape 3. Curr Protoc Bioinformatics. 2014;47(1):8–13.
Oostra M, Hagemeijer MC, van Gent M, Bekker CP, te Lintelo EG, Rottier PJ, et al. Topology and membrane anchoring of the coronavirus replication complex: not all hydrophobic domains of nsp3 and nsp6 are membrane spanning. J Virol. 2008;82(24):12392–405.
Snijder EJ, van der Meer Y, Zevenhoven-Dobbe J, Onderwater JJ, van der Meulen J, Koerten HK, et al. Ultrastructure and origin of membrane vesicles associated with the severe acute respiratory syndrome coronavirus replication complex. J Virol. 2006;80(12):5927–40.
Knoops K, Kikkert M, Worm SH, Zevenhoven-Dobbe JC, van der Meer Y, Koster AJ, et al. SARS-coronavirus replication is supported by a reticulovesicular network of modified endoplasmic reticulum. PLoS Biol. 2008;6(9):e226.
Zhang J, Cruz-Cosme R, Zhuang MW, Liu D, Liu Y, Teng S, et al. A systemic and molecular study of subcellular localization of SARS-CoV-2 proteins. Signal Transduct Target Ther. 2020;5(1):269.
Thul PJ, Akesson L, Wiking M, Mahdessian D, Geladaki A, Ait Blal H, et al. A subcellular map of the human proteome. Science. 2017;356(6340).
Riccio AA, Cingolani G, Pascal JM. PARP-2 domain requirements for DNA damage-dependent activation and localization to sites of DNA damage. Nucleic Acids Res. 2016;44(4):1691–702.
Liang D, Tian L, You R, Halpert MM, Konduri V, Baig YC, et al. AIMp1 potentiates TH1 polarization and is critical for effective antitumor and antiviral immunity. Front Immunol. 2017;8:1801.
Kim MS, Kim S, Myung H. Degradation of AIMP1/p43 induced by hepatitis C virus E2 leads to upregulation of TGF-beta signaling and increase in surface expression of gp96. PLoS One. 2014;9(5):e96302.
Bouhaddou M, Memon D, Meyer B, White KM, Rezelj VV, Correa Marrero M, et al. The global phosphorylation landscape of SARS-CoV-2 infection. Cell. 2020;182(3):685–712 e19.
Spalinger MR, Hai R, Li J, Santos AN, Nordgren TM, Tremblay ML, et al. Identification of a novel susceptibility marker for SARS-CoV-2 infection in human subjects and risk mitigation with a clinically approved JAK inhibitor in human/mouse cells. bioRxiv. 2020.
Banerjee AK, Blanco MR, Bruce EA, Honson DD, Chen LM, Chow A, et al. SARS-CoV-2 disrupts splicing, translation, and protein trafficking to suppress host defenses. Cell. 2020;183(5):1325–39 e21.
Finkel Y, Gluck A, Nachshon A, Winkler R, Fisher T, Rozman B, et al. SARS-CoV-2 uses a multipronged strategy to impede host protein synthesis. Nature. 2021;594(7862):240–5.
Ostrowski M, Carmo NB, Krumeich S, Fanget I, Raposo G, Savina A, et al. Rab27a and Rab27b control different steps of the exosome secretion pathway. Nat Cell Biol. 2010;12(1):19–30.
Hutchison JM, Capone R, Luu DD, Shah KH, Hadziselimovic A, Van Horn WD, et al. Recombinant SARS-CoV-2 envelope protein traffics to the trans-Golgi network following amphipol-mediated delivery into human cells. J Biol Chem. 2021;297(2):100940.
Johnson JL, Brzezinska AA, Tolmachova T, Munafo DB, Ellis BA, Seabra MC, et al. Rab27a and Rab27b regulate neutrophil azurophilic granule exocytosis and NADPH oxidase activity by independent mechanisms. Traffic. 2010;11(4):533–47.
Xia T, Yi XM, Wu X, Shang J, Shu HB. PTPN1/2-mediated dephosphorylation of MITA/STING promotes its 20S proteasomal degradation and attenuates innate antiviral response. Proc Natl Acad Sci U S A. 2019;116(40):20063–9.
Nishiyama-Fujita Y, Shimizu T, Sagawa M, Uchida H, Kizaki M. The role of TC-PTP (PTPN2) in modulating sensitivity to imatinib and interferon-alpha in CML cell line, KT-1 cells. Leuk Res. 2013;37(9):1150–5.
Heinonen KM, Nestel FP, Newell EW, Charette G, Seemayer TA, Tremblay ML, et al. T-cell protein tyrosine phosphatase deletion results in progressive systemic inflammatory disease. Blood. 2004;103(9):3457–64.
Taschuk F, Cherry S. DEAD-box helicases: sensors, regulators, and effectors for antiviral defense. Viruses. 2020;12(2):181.
Wang P, Zhu S, Yang L, Cui S, Pan W, Jackson R, et al. Nlrp6 regulates intestinal antiviral innate immunity. Science. 2015;350(6262):826–30.
Mosallanejad K, Sekine Y, Ishikura-Kinoshita S, Kumagai K, Nagano T, Matsuzawa A, et al. The DEAH-box RNA helicase DHX15 activates NF-kappaB and MAPK signaling downstream of MAVS during antiviral responses. Sci Signal. 2014;7(323):ra40.
Pattabhi S, Knoll ML, Gale M Jr, Loo YM. DHX15 is a Coreceptor for RLR signaling that promotes antiviral defense against RNA virus infection. J Interf Cytokine Res. 2019;39(6):331–46.
Wang Y, Wu X, Ge R, Song L, Li K, Tian S, et al. Global Screening of Sentrin-Specific Protease Family Substrates in SUMOylation. bioRxiv. 2020; 2020.02.25.964072.
Wang L, Zhao W, Zhang M, Wang P, Zhao K, Zhao X, et al. USP4 positively regulates RIG-I-mediated antiviral response through deubiquitination and stabilization of RIG-I. J Virol. 2013;87(8):4507–15.
Lord CJ, Ashworth A. PARP inhibitors: synthetic lethality in the clinic. Science. 2017;355(6330):1152–8.
Yelamos J, Monreal Y, Saenz L, Aguado E, Schreiber V, Mota R, et al. PARP-2 deficiency affects the survival of CD4+CD8+ double-positive thymocytes. EMBO J. 2006;25(18):4350–60.
Obaji E, Haikarainen T, Lehtio L. Characterization of the DNA dependent activation of human ARTD2/PARP2. Sci Rep. 2016;6:34487.
Benchoua A, Couriaud C, Guegan C, Tartier L, Couvert P, Friocourt G, et al. Active caspase-8 translocates into the nucleus of apoptotic cells to inactivate poly(ADP-ribose) polymerase-2. J Biol Chem. 2002;277(37):34217–22.
Li S, Zhang Y, Guan Z, Li H, Ye M, Chen X, et al. SARS-CoV-2 triggers inflammatory responses and cell death through caspase-8 activation. Signal Transduct Target Ther. 2020;5(1):235.
Lin CW, Lin KH, Hsieh TH, Shiu SY, Li JY. Severe acute respiratory syndrome coronavirus 3C-like protease-induced apoptosis. FEMS Immunol Med Microbiol. 2006;46(3):375–80.
Prescott L. SARS-CoV-2 3CLpro whole human proteome cleavage prediction and enrichment/depletion analysis. bioRxiv. 2020; 2020.08.24.265645.
Bienert S, Waterhouse A, de Beer TA, Tauriello G, Studer G, Bordoli L, et al. The SWISS-MODEL repository-new features and functionality. Nucleic Acids Res. 2017;45(D1):D313–D9.
R Core Team. R: A language and environment for statistical computing Vienna, Austria. 2020. Available from: https://www.R-project.org/.
Wickham H, Averick M, Bryan J, Chang W, McGowan LD, François R, et al. Welcome to the Tidyverse. J Open Source Softw. 2019;4(43):1686.
Slowikowski K. ggrepel: Automatically Position Non-Overlapping Text Labels with 'ggplot2'. R package version 0.8.2. 2020. Available from: https://CRAN.R-project.org/package=ggrepel.
Uhlen M, Fagerberg L, Hallstrom BM, Lindskog C, Oksvold P, Mardinoglu A, et al. Proteomics. Tissue-based map of the human proteome. Science. 2015;347(6220):1260419.
We wish to thank the Crowdfight COVID (crowdfightcovid19.org) for facilitating this collaboration, and Ian M. Boyce for assistance editing figures. The National Institute of Standards and Technology (NIST) notes that certain software and materials are identified in this article to explain a procedure as completely as possible. This does not imply a recommendation or endorsement by NIST, nor does it imply that the particular materials are necessarily the best available for the purpose. The opinions expressed in this article are the authors’ own and do not necessarily represent the views of NIST.
BMS was an International Associate of NIST, supported through the Professional Research Experience Program agreement with the University of Maryland. PDT was supported by a National Research Council Postdoctoral Research Associateship. This study was not funded by any specific grant.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1. Supplementary Tables S1-S16
This is the raw data output by NetCorona following analysis of the “All Human Proteins” dataset.
This is the raw data output by NetCorona following analysis of the “One Protein Per Gene” dataset.
This is the raw data output by NetCorona following analysis of the “Proteins With PDB” dataset.
Interaction scores vs max NetCorona score.
Raw data displaying how predicted cleaved motifs were matched to a PDB file.
ASA vs NetCorona score.
About this article
Cite this article
Scott, B.M., Lacasse, V., Blom, D.G. et al. Predicted coronavirus Nsp5 protease cleavage sites in the human proteome. BMC Genom Data 23, 25 (2022). https://doi.org/10.1186/s12863-022-01044-y
- Human proteins
- Human proteome