Potent, Novel SARS-CoV-2 PLpro Inhibitors Block Viral Replication in Monkey and Human Cell Cultures

Targeting SARS-CoV-2 nonstructural protein 3: Function, structure, inhibition, and perspective in drug discovery

As a highly contagious human pathogen, severe acute respiratory syndrome-associated coronavirus-2 (SARS-CoV-2) has infected billions of people worldwide with more than 6 million deaths. With several effective vaccines and antiviral drugs now available, the SARS-CoV-2 pandemic been brought under control. However, a new pathogenic coronavirus could emerge in the future, given the zoonotic nature of this virus. Natural evolution and drug-induced mutations of SARS-CoV-2 also require continued efforts for new anti-coronavirus drugs. Nonstructural protein (nsp) 3 of CoVs is a large, multifunctional protein, containing a papain-like protease (PLpro) and a macrodomain (Mac1), which are essential for viral replication. Here, we provide a comprehensive review of the function, structure, and inhibition of SARS-CoV/-CoV-2 PLpro and Mac1. We also discuss advances in, and challenges to, the discovery of drugs against these targets.

Discovery of Chalcone-Based Hybrid Structures as High Affinity and Site-Specific Inhibitors against SARS-CoV-2: A Comprehensive Structural Analysis Based on Various Host-Based and Viral Targets

Previous studies indicated that natural-based chalcones have significant inhibitory effects on the coronavirus enzymes 3CLpro and PLpro as well as modulation of some host-based antiviral targets (HBATs). In this study, a comprehensive computational and structural study was performed to investigate the affinity of our compound library consisting of 757 chalcone-based structures (CHA-1 to CHA-757) for inhibiting the 3CLpro and PLpro enzymes and against twelve selected host-based targets. Our results indicated that CHA-12 (VUF 4819) is the most potent and multi-target inhibitor in our chemical library over all viral and host-based targets. Correspondingly, CHA-384 and its congeners containing ureide moieties were found to be potent and selective 3CLpro inhibitors, and benzotriazole moiety in CHA-37 was found to be a main fragment for inhibiting the 3CLpro and PLpro. Surprisingly, our results indicate that the ureide and sulfonamide moieties are integral fragments for the optimum 3CLpro inhibition while occupying the S1 and S3 subsites, which is fully consistent with recent reports on the site-specific 3CLpro inhibitors. Finding the multi-target inhibitor CHA-12, previously reported as an LTD4 antagonist for the treatment of inflammatory pulmonary diseases, prompted us to suggest it as a concomitant agent for relieving respiratory symptoms and suppressing COVID-19 infection.

Biophysical Interpretation of Evolutionary Consequences on the SARS-CoV2 Main Protease through Molecular Dynamics Simulations and Network Topology Analysis

In this study, we present a combined analysis procedure between atomistic molecular dynamics (MD) simulations and network topology to obtain more understanding on the evolutionary consequences on protein stability and substrate binding of the main protease enzyme of SARS-CoV2. Communicability matrices of the protein residue networks (PRNs) were extracted from MD trajectories of both Mpro enzymes in complex with the nsp8/9 peptide substrate to compare the local communicability within both proteases that would affect the enzyme function, along with biophysical details on global protein conformation, flexibility, and contribution of amino acid side chains to both intramolecular and intermolecular interactions. The analysis displayed the significance of the mutated residue 46 with the highest communicability gain to the binding pocket closure. Interestingly, the mutated residue 134 with the highest communicability loss corresponded to a local structural disruption of the adjacent peptide loop. The enhanced flexibility of the disrupted loop connecting to the catalytic residue Cys145 introduced an extra binding mode that brought the substrate in proximity and could facilitate the reaction. This understanding might provide further help in the drug development strategy against SARS-CoV2 and prove the capability of the combined techniques of MD simulations and network topology analysis as a "reverse" protein engineering tool.

Large-Scale Virtual Screening for the Discovery of SARS-CoV-2 Papain-like Protease (PLpro) Non-covalent Inhibitors

The rapid global spread of the SARS-CoV-2 virus facilitated the development of novel direct-acting antiviral agents (DAAs). The papain-like protease (PLpro) has been proposed as one of the major SARS-CoV-2 targets for DAAs due to its dual role in processing viral proteins and facilitating the host's immune suppression. This dual role makes identifying small molecules that can effectively neutralize SARS-CoV-2 PLpro activity a high-priority task. However, PLpro drug discovery faces a significant challenge due to the high mobility and induced-fit effects in the protease's active site. Herein, we virtually screened the ZINC20 database with Deep Docking (DD) to identify prospective noncovalent PLpro binders and combined ultra-large consensus docking with two pharmacophore (ph4)-filtering strategies. The analysis of active compounds revealed their somewhat-limited diversity, likely attributed to the induced-fit nature of PLpro's active site in the crystal structures, and therefore, the use of rigid docking protocols poses inherited limitations. The top hits were assessed against recombinant viral proteins and live viruses, demonstrating desirable inhibitory activities. The best compound VPC-300195 (IC₅₀: 15 μM) ranks among the top noncovalent PLpro inhibitors discovered through in silico methodologies. In the search for novel SARS-CoV-2 PLpro-specific chemotypes, the identified inhibitors could serve as diverse templates for the development of effective noncovalent PLpro inhibitors.

Language models for the prediction of SARS-CoV-2 inhibitors

The COVID-19 pandemic highlights the need for computational tools to automate and accelerate drug design for novel protein targets. We leverage deep learning language models to generate and score drug candidates based on predicted protein binding affinity. We pre-trained a deep learning language model (BERT) on ∼9.6 billion molecules and achieved peak performance of 603 petaflops in mixed precision. Our work reduces pre-training time from days to hours, compared to previous efforts with this architecture, while also increasing the dataset size by nearly an order of magnitude. For scoring, we fine-tuned the language model using an assembled set of thousands of protein targets with binding affinity data and searched for inhibitors of specific protein targets, SARS-CoV-2 Mpro and PLpro. We utilized a genetic algorithm approach for finding optimal candidates using the generation and scoring capabilities of the language model. Our generalizable models accelerate the identification of inhibitors for emerging therapeutic targets.

Hydrogen bonding penalty used for virtual screening to discover potent inhibitors for Papain-Like cysteine proteases of SARS-CoV-2

The Papain-Like proteases (PLpro) of SARS-CoV-2 play a crucial role in viral replication and the formation of nonstructural proteins. To find available inhibitors, the 3D structure of PLpro of SARS2 was obtained by homologous modelling, and we used this structure as a target to search for inhibitors through molecular docking and MM/GBSA binding free energy rescoring. A novel hydrogen bonding penalty was applied to the screening process, which meanwhile took desolvation into account. Finally, 61 compounds were acquired and 4 of them with IC50 at micromolar level tested in vitro enzyme activity assay, which includes clinical drugs tegaserod. Considering the importance of crystal water molecules, the 4 compounds were re-docked and considered bound waters in the active site as a part of PLpro. The binding modes of these 4 compounds were further explored with metadynamics simulations. The hits will provide a starting point for future key interactions identified and lead optimization targetting PLpro.

Parthenolide reveals an allosteric mode to inhibit the deISGylation activity of SARS-CoV‑2 papain-like protease

The coronavirus papain-like protease (PLpro) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is responsible for viral polypeptide cleavage and the deISGylation of interferon-stimulated gene 15 (ISG15), which enable it to participate in virus replication and host innate immune pathways. Therefore, PLpro is considered an attractive antiviral drug target. Here, we show that parthenolide, a germacrane sesquiterpene lactone, has SARS-CoV-2 PLpro inhibitory activity. Parthenolide covalently binds to Cys-191 or Cys-194 of the PLpro protein, but not the Cys-111 at the PLpro catalytic site. Mutation of Cys-191 or Cys-194 reduces the activity of PLpro. Molecular docking studies show that parthenolide may also form hydrogen bonds with Lys-192, Thr-193, and Gln-231. Furthermore, parthenolide inhibits the deISGylation but not the deubiquitinating activity of PLpro in vitro. These results reveal that parthenolide inhibits PLpro activity by allosteric regulation.

Inhibitor induced conformational changes in SARS-COV-2 papain-like protease

SARS-CoV-2's papain-like protease (PLpro) interaction with ligands has recently been explored with a myriad of crystal structures. We used molecular dynamics (MD) simulations to study different PLpro-ligand complexes, their ligand-induced conformational changes, and interactions. We focused on inhibitors reported with known IC50 against PLpro, namely GRL-0617, XR8-89, PLP_Snyder530, and Sander's recently published compound 7 (CPD7), and compared these trajectories against the apostructure (Apo), with a total of around 60 µs worth simulation data. We aimed to study the conformational changes using molecular dynamics simulations for the inhibitors in the PLpro. PCA analyses and the MSM models revealed distinct conformations of PLpro in the absence/presence of ligands and proposed that BL2-loop contributes to the accessibility of these inhibitors. Further, bulkier substituents closer to Tyr268 and Gln269 could improve inhibition of SARS-CoV-2 PLpro by occupying the region between BL2-groove and BL2-loop, but we also expand on the relevance of exploring multiple PLpro sub-pockets to improve inhibition.

Exploring Noncovalent Protease Inhibitors for the Treatment of Severe Acute Respiratory Syndrome and Severe Acute Respiratory Syndrome-Like Coronaviruses

Over the last 20 years, both severe acute respiratory syndrome coronavirus-1 and severe acute respiratory syndrome coronavirus-2 have transmitted from animal hosts to humans causing zoonotic outbreaks of severe disease. Both viruses originate from a group of betacoronaviruses known as subgroup 2b. The emergence of two dangerous human pathogens from this group along with previous studies illustrating the potential of other subgroup 2b members to transmit to humans has underscored the need for antiviral development against them. Coronaviruses modify the host innate immune response in part through the reversal of ubiquitination and ISGylation with their papain-like protease (PLpro). To identify unique or overarching subgroup 2b structural features or enzymatic biases, the PLpro from a subgroup 2b bat coronavirus, BtSCoV-Rf1.2004, was biochemically and structurally evaluated. This evaluation revealed that PLpros from subgroup 2b coronaviruses have narrow substrate specificity for K48 polyubiquitin and ISG15 originating from certain species. The PLpro of BtSCoV-Rf1.2004 was used as a tool alongside PLpro of CoV-1 and CoV-2 to design 30 novel noncovalent drug-like pan subgroup 2b PLpro inhibitors that included determining the effects of using previously unexplored core linkers within these compounds. Two crystal structures of BtSCoV-Rf1.2004 PLpro bound to these inhibitors aided in compound design as well as shared structural features among subgroup 2b proteases. Screening of these three subgroup 2b PLpros against this novel set of inhibitors along with cytotoxicity studies provide new directions for pan-coronavirus subgroup 2b antiviral development of PLpro inhibitors.

Targeting SARS-CoV-2 Proteases for COVID-19 Antiviral Development

The emergence of severe acute respiratory syndrome (SARS-CoV-2) in 2019 marked the third occurrence of a highly pathogenic coronavirus in the human population since 2003. As the death toll surpasses 5 million globally and economic losses continue, designing drugs that could curtail infection and disease progression is critical. In the US, three highly effective Food and Drug Administration (FDA)-authorized vaccines are currently available, and Remdesivir is approved for the treatment of hospitalized patients. However, moderate vaccination rates and the sustained evolution of new viral variants necessitate the ongoing search for new antivirals. Several viral proteins have been prioritized as SARS-CoV-2 antiviral drug targets, among them the papain-like protease (PLpro) and the main protease (Mpro). Inhibition of these proteases would target viral replication, viral maturation, and suppression of host innate immune responses. Knowledge of inhibitors and assays for viruses were quickly adopted for SARS-CoV-2 protease research. Potential candidates have been identified to show inhibitory effects against PLpro and Mpro, both in biochemical assays and viral replication in cells. These results encourage further optimizations to improve prophylactic and therapeutic efficacy. In this review, we examine the latest developments of potential small-molecule inhibitors and peptide inhibitors for PLpro and Mpro, and how structural biology greatly facilitates this process.

Structure-based inhibitor design and repurposing clinical drugs to target SARS-CoV-2 proteases

SARS-CoV-2, the coronavirus responsible for the current COVID-19 pandemic, encodes two proteases, 3CLpro and PLpro, two of the main antiviral research targets. Here we provide an overview of the structures and functions of 3CLpro and PLpro and examine strategies of structure-based drug designing and drug repurposing against these proteases. Rational structure-based drug design enables the generation of potent and target-specific antivirals. Drug repurposing offers an attractive prospect with an accelerated turnaround. Thus far, several protease inhibitors have been identified, and some candidates are undergoing trials that may well prove to be effective antivirals against SARS-CoV-2.

Natural biflavones are potent inhibitors against SARS-CoV-2 papain-like protease

Papain-like protease (PL^pro) is a key enzyme encoded by SARS-CoV-2 that is essential for viral replication and immune evasion. Significant suppression of viral spread and promotion of antiviral immunity can be achieved by inhibition of PL^pro, revealing an inspiring strategy for COVID-19 treatment. This study aimed to discover PL^pro inhibitors by investigating the national compound library of traditional Chinese medicines (NCLTCMs), a phytochemical library comprising over 9000 TCM-derived compounds. Through virtual screening and enzymatic evaluations, nine natural biflavones were confirmed to be effective PL^pro inhibitors with IC₅₀ values ranging from 9.5 to 43.2 μM. Pro-ISG15 cleavage assays further demonstrated that several biflavones exhibited potent inhibitory effects against PL^pro-mediated deISGylation, a key process involved in viral immune evasion. Herein, we report the discovery, antiviral evaluation, structure-activity relationship elucidation and molecular docking investigation of biflavones as potent inhibitors of SARS-CoV-2 PL^pro.

In Silico Elucidation of Potent Inhibitors and Rational Drug Design against SARS-CoV-2 Papain-like Protease

Global public health has been a critical problem by the sudden increase of the COVID-19 outbreak. The papain-like protease (PLpro) of SARS-CoV-2 is a key promising target for antiviral drug development since it plays a pivotal role in viral replication and innate immunity. Here, we employed the all-atom molecular dynamics (MD) simulations and binding free energy calculations based on MM-PB(GB)SA and SIE methods to elucidate and compare the binding behaviors of five inhibitors derived from peptidomimetic inhibitors (VIR250 and VIR251) and naphthalene-based inhibitors (GRL-0617, compound 3, and compound Y96) against SARS-CoV-2 PLpro. The obtained results revealed that all inhibitors interacting within the PLpro active site are mostly driven by vdW interactions, and the hydrogen bond formation in residues G163 and G271 with peptidomimetics and the Q269 residue with naphthalene-based inhibitors was essential for stabilizing the protein-ligand complexes. Among the five studied inhibitors, VIR250 exhibited the most binding efficiency with SARS-CoV-2 PLpro, and thus, it was chosen for the rational drug design. Based on the computationally designed ligand-protein complexes, the replacement of aromatic rings including heteroatoms (e.g., thiazolopyridine) at the P2 and P4 sites could help to improve the inhibitor-binding efficiency. Furthermore, the hydrophobic interactions with residues at P1-P3 sites can be increased by enlarging the nonpolar moieties (e.g., ethene) at the N-terminal of VIR250. We expect that the structural data obtained will contribute to the development of new PLpro inhibitors with more inhibitory potency for COVID-19.

Analyzing the Effect of Mutations in SARS-CoV2 Papain-Like Protease from Saudi Isolates on Protein Structure and Drug-Protein Binding: Molecular Modelling and Dynamics Studies

The continuous and rapid development of the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) virus remains a health concern especially with the emergence of numerous variants and mutations worldwide. As with other RNA viruses, SARS-CoV-2 has a genetically high mutation rate. These mutations have an impact on the virus characteristics, including transmissibility, antigenicity and development of drug and vaccine resistance. This work was pursued to identify the differences that exist in the papain-like protease (PL^Pro) from 58 Saudi isolates in comparison to the first reported sequence from Wuhan, China and determine their implications on protein structure and the inhibitor binding. PL^pro is a key protease enzyme for the host cells invasion and viral proteolytic cleavage, hence, it emerges as a valuable antiviral therapeutic target. Two mutations were identified including D108G and A249V and shown to increase the molecular flexibility of PL^Pro protein and alter the protein stability, particularly with D108G mutation. The effect of these mutations on the stability and dynamic behavior of PL^Pro structures as well as their effect on the binding of a known inhibitor; GRL0617 were further investigated by molecular docking and dynamic simulation.

SARS-CoV-2 Papain-Like Protease Potential Inhibitors-In Silico Quantitative Assessment

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) encodes the papain-like protease (PLpro). The protein not only plays an essential role in viral replication but also cleaves ubiquitin and ubiquitin-like interferon-stimulated gene 15 protein (ISG15) from host proteins, making it an important target for developing new antiviral drugs. In this study, we searched for novel, noncovalent potential PLpro inhibitors by employing a multistep in silico screening of a 15 million compound library. The selectivity of the best-scored compounds was evaluated by checking their binding affinity to the human ubiquitin carboxy-terminal hydrolase L1 (UCH-L1), which, as a deubiquitylating enzyme, exhibits structural and functional similarities to the PLpro. As a result, we identified 387 potential, selective PLpro inhibitors, from which we retrieved the 20 best compounds according to their IC50 values toward PLpro estimated by a multiple linear regression model. The selected candidates display potential activity against the protein with IC50 values in the nanomolar range from approximately 159 to 505 nM and mostly adopt a similar binding mode to the known, noncovalent SARS-CoV-2 PLpro inhibitors. We further propose the six most promising compounds for future in vitro evaluation. The results for the top potential PLpro inhibitors are deposited in the database prepared to facilitate research on anti-SARS-CoV-2 drugs.