Papain-like protease regulates SARS-CoV-2 viral spread and innate immunity.

The papain-like protease PLpro is an essential coronavirus enzyme that is required for processing viral polyproteins to generate a functional replicase complex and enable viral spread1,2. PLpro is also implicated in cleaving proteinaceous post-translational modifications on host proteins as an evasion mechanism against host antiviral immune responses3-5. Here we perform biochemical, structural and functional characterization of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) PLpro (SCoV2-PLpro) and outline differences with SARS-CoV PLpro (SCoV-PLpro) in regulation of host interferon and NF-κB pathways. SCoV2-PLpro and SCoV-PLpro share 83% sequence identity but exhibit different host substrate preferences; SCoV2-PLpro preferentially cleaves the ubiquitin-like interferon-stimulated gene 15 protein (ISG15), whereas SCoV-PLpro predominantly targets ubiquitin chains. The crystal structure of SCoV2-PLpro in complex with ISG15 reveals distinctive interactions with the amino-terminal ubiquitin-like domain of ISG15, highlighting the high affinity and specificity of these interactions. Furthermore, upon infection, SCoV2-PLpro contributes to the cleavage of ISG15 from interferon responsive factor 3 (IRF3) and attenuates type I interferon responses. Notably, inhibition of SCoV2-PLpro with GRL-0617 impairs the virus-induced cytopathogenic effect, maintains the antiviral interferon pathway and reduces viral replication in infected cells. These results highlight a potential dual therapeutic strategy in which targeting of SCoV2-PLpro can suppress SARS-CoV-2 infection and promote antiviral immunity.

The SARS-unique domain (SUD) of SARS-CoV and SARS-CoV-2 interacts with human Paip1 to enhance viral RNA translation.

The ongoing outbreak of severe acute respiratory syndrome (SARS) coronavirus 2 (SARS-CoV-2) demonstrates the continuous threat of emerging coronaviruses (CoVs) to public health. SARS-CoV-2 and SARS-CoV share an otherwise non-conserved part of non-structural protein 3 (Nsp3), therefore named as "SARS-unique domain" (SUD). We previously found a yeast-2-hybrid screen interaction of the SARS-CoV SUD with human poly(A)-binding protein (PABP)-interacting protein 1 (Paip1), a stimulator of protein translation. Here, we validate SARS-CoV SUD:Paip1 interaction by size-exclusion chromatography, split-yellow fluorescent protein, and co-immunoprecipitation assays, and confirm such interaction also between the corresponding domain of SARS-CoV-2 and Paip1. The three-dimensional structure of the N-terminal domain of SARS-CoV SUD ("macrodomain II", Mac2) in complex with the middle domain of Paip1, determined by X-ray crystallography and small-angle X-ray scattering, provides insights into the structural determinants of the complex formation. In cellulo, SUD enhances synthesis of viral but not host proteins via binding to Paip1 in pBAC-SARS-CoV replicon-transfected cells. We propose a possible mechanism for stimulation of viral translation by the SUD of SARS-CoV and SARS-CoV-2.

Oligomeric assembly of the C-terminal and transmembrane region of SARS-CoV-2 nsp3.

As for all single-stranded, positive-sense RNA (+RNA) viruses, intracellular RNA synthesis relies on extensive remodeling of host cell membranes that leads to the formation of specialized structures. In the case of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) coronavirus causing COVID-19, endoplasmic reticulum membranes are modified, resulting in the formation of double-membrane vesicles (DMVs), which contain the viral dsRNA intermediate and constitute membrane-bound replication organelles. The non-structural and transmembrane protein nsp3 is a key player in the biogenesis of DMVs and, therefore, represents an interesting antiviral target. However, as an integral transmembrane protein, it is challenging to express for structural biology. The C-terminus of nsp3 encompasses all the membrane-spanning, -interacting, and -remodeling elements. By using a cell-free expression system, we successfully produced the C-terminal region of nsp3 (nsp3C) and reconstituted purified nsp3C into phospholipid nanodiscs, opening the way for structural studies. Negative-stain transmission electron microscopy revealed the presence of nsp3C oligomers very similar to the region abutting and spanning the membrane on the cytosolic side of DMVs in a recent subtomogram average of the SARS-CoV-2 nsp3-4 pore (1). AlphaFold-predicted structural models fit particularly well with our experimental data and support a pore-forming hexameric assembly. Altogether, our data give unprecedented clues to understand the structural organization of nsp3, the principal component that shapes the molecular pore that spans the DMVs and is required for the export of RNA in vivo. IMPORTANCE: Membrane remodeling is at the heart of intracellular replication for single-stranded, positive-sense RNA viruses. In the case of coronaviruses, including severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), this leads to the formation of a network of double-membrane vesicles (DMVs). Targeting DMV biogenesis offers promising prospects for antiviral therapies. This requires a better understanding of the molecular mechanisms and proteins involved. Three non-structural proteins (nsp3, nsp4, and nsp6) direct the intracellular membrane rearrangements upon SARS-CoV-2 infection. All of them contain transmembrane helices. The nsp3 component, the largest and multi-functional protein of the virus, plays an essential role in this process. Aiming to understand its structural organization, we used a cell-free protein synthesis assay to produce and reconstitute the C-terminal part of nsp3 (nsp3C) including transmembrane domains into phospholipid nanodiscs. Our work reveals the oligomeric organization of one key player in the biogenesis of SARS-CoV-2 DMVs, providing basis for the design of future antiviral strategies.

A single inactivating amino acid change in the SARS-CoV-2 NSP3 Mac1 domain attenuates viral replication in vivo.

Despite unprecedented efforts, our therapeutic arsenal against SARS-CoV-2 remains limited. The conserved macrodomain 1 (Mac1) in NSP3 is an enzyme exhibiting ADP-ribosylhydrolase activity and a possible drug target. To determine the role of Mac1 catalytic activity in viral replication, we generated recombinant viruses and replicons encoding a catalytically inactive NSP3 Mac1 domain by mutating a critical asparagine in the active site. While substitution to alanine (N40A) reduced catalytic activity by ~10-fold, mutations to aspartic acid (N40D) reduced activity by ~100-fold relative to wild-type. Importantly, the N40A mutation rendered Mac1 unstable in vitro and lowered expression levels in bacterial and mammalian cells. When incorporated into SARS-CoV-2 molecular clones, the N40D mutant only modestly affected viral fitness in immortalized cell lines, but reduced viral replication in human airway organoids by 10-fold. In mice, the N40D mutant replicated at >1000-fold lower levels compared to the wild-type virus while inducing a robust interferon response; all animals infected with the mutant virus survived infection. Our data validate the critical role of SARS-CoV-2 NSP3 Mac1 catalytic activity in viral replication and as a promising therapeutic target to develop antivirals.

Design of inhibitors of SARS-CoV-2 papain-like protease deriving from GRL0617: Structure-activity relationships.

The unique and complex structure of papain-like protease (PLpro) of the SARS-CoV-2 virus represents a difficult challenge for antiviral development, yet it offers a compelling validated target for effective therapy of COVID-19. The surge in scientific interest in inhibiting this cysteine protease emerged after its demonstrated connection to the cytokine storm in patients with COVID-19 disease. Furthermore, the development of new inhibitors against PLpro may also be beneficial for the treatment of respiratory infections caused by emerging coronavirus variants of concern. This review article provides a comprehensive overview of PLpro inhibitors, focusing on the structural framework of the known inhibitor GRL0617 and its analogs. We categorize PLpro inhibitors on the basis of their structures and binding site: Glu167 containing site, BL2 groove, Val70Ub site, and Cys111 containing catalytic site. We summarize and evaluate the majority of GRL0617-like inhibitors synthesized so far, highlighting their published biochemical parameters, which reflect their efficacy. Published research has shown that strategic modifications to GRL0617, such as decorating the naphthalene ring, extending the aromatic amino group or the orthomethyl group, can substantially decrease the IC50 from micromolar up to nanomolar concentration range. Some advantageous modifications significantly enhance inhibitory activity, paving the way for the development of new potent compounds. Our review places special emphasis on structures that involve direct modifications to the GRL0617 scaffold, including piperidine carboxamides and modified benzylmethylnaphthylethanamines (Jun9 scaffold). All these compounds are believed to inhibit the proteolytic, deubiquitination, and deISGylation activity of PLpro, biochemical processes linked to the severe progression of COVID-19. Finally, we summarize the development efforts for SARS-CoV-2 PLpro inhibitors, in detailed structure-activity relationships diagrams. This aims to inform and inspire future research in the search for potent antiviral agents against PLpro of current and emerging coronavirus threats.

In-cell bioluminescence resonance energy transfer (BRET)-based assay uncovers ceritinib and CA-074 as SARS-CoV-2 papain-like protease inhibitors.

Papain-like protease (PLpro) is an attractive anti-coronavirus target. The development of PLpro inhibitors, however, is hampered by the limitations of the existing PLpro assay and the scarcity of validated active compounds. We developed a novel in-cell PLpro assay based on BRET and used it to evaluate and discover SARS-CoV-2 PLpro inhibitors. The developed assay demonstrated remarkable sensitivity for detecting the reduction of intracellular PLpro activity while presenting high reliability and performance for inhibitor evaluation and high-throughput screening. Using this assay, three protease inhibitors were identified as novel PLpro inhibitors that are structurally disparate from those previously known. Subsequent enzymatic assays and ligand-protein interaction analysis based on molecular docking revealed that ceritinib directly inhibited PLpro, showing high geometric complementarity with the substrate-binding pocket in PLpro, whereas CA-074 methyl ester underwent intracellular hydrolysis, exposing a free carboxyhydroxyl group essential for hydrogen bonding with G266 in the BL2 groove, resulting in PLpro inhibition.

Postinfection treatment with a protease inhibitor increases survival of mice with a fatal SARS-CoV-2 infection.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection continues to be a serious global public health threat. The 3C-like protease (3CLpro) is a virus protease encoded by SARS-CoV-2, which is essential for virus replication. We have previously reported a series of small-molecule 3CLpro inhibitors effective for inhibiting replication of human coronaviruses including SARS-CoV-2 in cell culture and in animal models. Here we generated a series of deuterated variants of a 3CLpro inhibitor, GC376, and evaluated the antiviral effect against SARS-CoV-2. The deuterated GC376 displayed potent inhibitory activity against SARS-CoV-2 in the enzyme- and the cell-based assays. The K18-hACE2 mice develop mild to lethal infection commensurate with SARS-CoV-2 challenge doses and were proposed as a model for efficacy testing of antiviral agents. We treated lethally infected mice with a deuterated derivative of GC376. Treatment of K18-hACE2 mice at 24 h postinfection with a derivative (compound 2) resulted in increased survival of mice compared to vehicle-treated mice. Lung virus titers were decreased, and histopathological changes were ameliorated in compound 2-treated mice compared to vehicle-treated mice. Structural investigation using high-resolution crystallography illuminated binding interactions of 3CLpro of SARS-CoV-2 and SARS-CoV with deuterated variants of GC376. Taken together, deuterated GC376 variants have excellent potential as antiviral agents against SARS-CoV-2.

Structural and Biochemical Characterization of Porcine Epidemic Diarrhea Virus Papain-Like Protease 2.

Coronaviral papain-like proteases (PLpros) are essential enzymes that mediate not only the proteolytic processes of viral polyproteins during virus replication but also the deubiquitination and deISGylation of cellular proteins that attenuate host innate immune responses. Therefore, PLpros are attractive targets for antiviral drug development. Here, we report the crystal structure of papain-like protease 2 (PLP2) of porcine epidemic diarrhea virus (PEDV) in complex with ubiquitin (Ub). The X-ray structural analyses reveal that PEDV PLP2 interacts with the Ub substrate mainly through the Ub core region and C-terminal tail. Mutations of Ub-interacting residues resulted in a moderately or completely abolished deubiquitinylating function of PEDV PLP2. In addition, our analyses also indicate that 2-residue-extended blocking loop 2 at the S4 subsite contributes to the substrate selectivity and binding affinity of PEDV PLP2. Furthermore, the PEDV PLP2 Glu99 residue, conserved in alphacoronavirus PLpros, was found to govern the preference of a positively charged P4 residue of peptidyl substrates. Collectively, our data provided structure-based information for the substrate binding and selectivity of PEDV PLP2. These findings may help us gain insights into the deubiquitinating (DUB) and proteolytic functions of PEDV PLP2 from a structural perspective. IMPORTANCE Current challenges in coronaviruses (CoVs) include a comprehensive understanding of the mechanistic effects of associated enzymes, including the 3C-like and papain-like proteases. We have previously reported that the PEDV PLP2 exhibits a broader substrate preference, superior DUB function, and inferior peptidase activity. However, the structural basis for these functions remains largely unclear. Here, we show the high-resolution X-ray crystal structure of PEDV PLP2 in complex with Ub. Integrated structural and biochemical analyses revealed that (i) three Ub core-interacting residues are essential for DUB function, (ii) 2-residue-elongated blocking loop 2 regulates substrate selectivity, and (iii) a conserved glutamate residue governs the substrate specificity of PEDV PLP2. Collectively, our findings provide not only structural insights into the catalytic mechanism of PEDV PLP2 but also a model for developing antiviral strategies.

SARS-CoV-2 Nsp3 unique domain SUD interacts with guanine quadruplexes and G4-ligands inhibit this interaction.

The multidomain non-structural protein 3 (Nsp3) is the largest protein encoded by coronavirus (CoV) genomes and several regions of this protein are essential for viral replication. Of note, SARS-CoV Nsp3 contains a SARS-Unique Domain (SUD), which can bind Guanine-rich non-canonical nucleic acid structures called G-quadruplexes (G4) and is essential for SARS-CoV replication. We show herein that the SARS-CoV-2 Nsp3 protein also contains a SUD domain that interacts with G4s. Indeed, interactions between SUD proteins and both DNA and RNA G4s were evidenced by G4 pull-down, Surface Plasmon Resonance and Homogenous Time Resolved Fluorescence. These interactions can be disrupted by mutations that prevent oligonucleotides from folding into G4 structures and, interestingly, by molecules known as specific ligands of these G4s. Structural models for these interactions are proposed and reveal significant differences with the crystallographic and modeled 3D structures of the SARS-CoV SUD-NM/G4 interaction. Altogether, our results pave the way for further studies on the role of SUD/G4 interactions during SARS-CoV-2 replication and the use of inhibitors of these interactions as potential antiviral compounds.

A Novel Ambroxol-Derived Tetrahydroquinazoline with a Potency against SARS-CoV-2 Proteins.

We report synthesis of a novel 1,2,3,4-tetrahydroquinazoline derivative, named 2-(6,8-dibromo-3-(4-hydroxycyclohexyl)-1,2,3,4-tetrahydroquinazolin-2-yl)phenol (1), which was obtained from the hydrochloride of 4-((2-amino-3,5-dibromobenzyl)amino)cyclohexan-1-ol (ambroxol hydrochloride) and salicylaldehyde in EtOH. The resulting compound was produced in the form of colorless crystals of the composition 1∙0.5EtOH. The formation of the single product was confirmed by the IR and 1H spectroscopy, single-crystal and powder X-ray diffraction, and elemental analysis. The molecule of 1 contains a chiral tertiary carbon of the 1,2,3,4-tetrahydropyrimidine fragment and the crystal structure of 1∙0.5EtOH is a racemate. Optical properties of 1∙0.5EtOH were revealed by UV-vis spectroscopy in MeOH and it was established that the compound absorbs exclusively in the UV region up to about 350 nm. 1∙0.5EtOH in MeOH exhibits dual emission and the emission spectra contains bands at about 340 and 446 nm upon excitation at 300 and 360 nm, respectively. The DFT calculations were performed to verify the structure as well as electronic and optical properties of 1. ADMET properties of the R-isomer of 1 were evaluated using the SwissADME, BOILED-Egg, and ProTox-II tools. As evidenced from the blue dot position in the BOILED-Egg plot, both human blood-brain barrier penetration and gastrointestinal absorption properties are positive with the positive PGP effect on the molecule. Molecular docking was applied to examine the influence of the structures of both R-isomer and S-isomer of 1 on a series of the SARS-CoV-2 proteins. According to the docking analysis results, both isomers of 1 were found to be active against all the applied SARS-CoV-2 proteins with the best binding affinities with Papain-like protease (PLpro) and nonstructural protein 3 (Nsp3_range 207-379-AMP). Ligand efficiency scores for both isomers of 1 inside the binding sites of the applied proteins were also revealed and compared with the initial ligands. Molecular dynamics simulations were also applied to evaluate the stability of complexes of both isomers with Papain-like protease (PLpro) and nonstructural protein 3 (Nsp3_range 207-379-AMP). The complex of the S-isomer with Papain-like protease (PLpro) was found to be highly unstable, while the other complexes are stable.

Isolation and In Silico Anti-SARS-CoV-2 Papain-Like Protease Potentialities of Two Rare 2-Phenoxychromone Derivatives from Artemisia spp.

Two rare 2-phenoxychromone derivatives, 6-demethoxy-4`-O-capillarsine (1) and tenuflorin C (2), were isolated from the areal parts of Artemisia commutata and A. glauca, respectively, for the first time. Being rare in nature, the inhibition potentialities of 1 and 2 against SARS-CoV-2 was investigated using multistage in silico techniques. At first, molecular similarity and fingerprint studies were conducted for 1 and 2 against co-crystallized ligands of eight different COVID-19 enzymes. The carried-out studies indicated the similarity of 1 and 2 with TTT, the co-crystallized ligand of COVID-19 Papain-Like Protease (PLP), (PDB ID: 3E9S). Therefore, molecular docking studies of 1 and 2 against the PLP were carried out and revealed correct binding inside the active site exhibiting binding energies of -18.86 and -18.37 Kcal/mol, respectively. Further, in silico ADMET in addition to toxicity evaluation of 1 and 2 against seven models indicated the general safety and the likeness of 1 and 2 to be drugs. Lastly, to authenticate the binding and to investigate the thermodynamic characters, molecular dynamics (MD) simulation studies were conducted on 1 and PLP.

Predicting cryptic ligand binding sites based on normal modes guided conformational sampling.

To greatly expand the druggable genome, fast and accurate predictions of cryptic sites for small molecules binding in target proteins are in high demand. In this study, we have developed a fast and simple conformational sampling scheme guided by normal modes solved from the coarse-grained elastic models followed by atomistic backbone refinement and side-chain repacking. Despite the observations of complex and diverse conformational changes associated with ligand binding, we found that simply sampling along each of the lowest 30 modes is near optimal for adequately restructuring cryptic sites so they can be detected by existing pocket finding programs like fpocket and concavity. We further trained machine-learning protocols to optimize the combination of the sampling-enhanced pocket scores with other dynamic and conservation scores, which only slightly improved the performance. As assessed based on a training set of 84 known cryptic sites and a test set of 14 proteins, our method achieved high accuracy of prediction (with area under the receiver operating characteristic curve >0.8) comparable to the CryptoSite server. Compared with CryptoSite and other methods based on extensive molecular dynamics simulation, our method is much faster (1-2 hours for an average-size protein) and simpler (using only pocket scores), so it is suitable for high-throughput processing of large datasets of protein structures at the genome scale.

Site-Specific Labeling of Endogenous Proteins Using CoLDR Chemistry.

Chemical modifications of native proteins can affect their stability, activity, interactions, localization, and more. However, there are few nongenetic methods for the installation of chemical modifications at a specific protein site in cells. Here we report a covalent ligand directed release (CoLDR) site-specific labeling strategy, which enables the installation of a variety of functional tags on a target protein while releasing the directing ligand. Using this approach, we were able to label various proteins such as BTK, K-RasG12C, and SARS-CoV-2 PLpro with different tags. For BTK we have shown selective labeling in cells of both alkyne and fluorophores tags. Protein labeling by traditional affinity methods often inhibits protein activity since the directing ligand permanently occupies the target binding pocket. We have shown that using CoLDR chemistry, modification of BTK by these probes in cells preserves its activity. We demonstrated several applications for this approach including determining the half-life of BTK in its native environment with minimal perturbation, as well as quantification of BTK degradation by a noncovalent proteolysis targeting chimera (PROTAC) by in-gel fluorescence. Using an environment-sensitive "turn-on" fluorescent probe, we were able to monitor ligand binding to the active site of BTK. Finally, we have demonstrated efficient CoLDR-based BTK PROTACs (DC50 < 100 nM), which installed a CRBN binder onto BTK. This approach joins very few available labeling strategies that maintain the target protein activity and thus makes an important addition to the toolbox of chemical biology.

Molecular docking of potential SARS-CoV-2 papain-like protease inhibitors.

SARS-CoV-2 papain-like protease is considered as an important potential target for anti-SARS-CoV-2 drug discovery due to its crucial roles in viral spread and innate immunity. Here, we have utilized an in silico molecular docking approach to identify the possible inhibitors of the SARS-CoV-2 papain-like protease, by screening 21 antiviral, antifungal and anticancer compounds. Among them, Neobavaisoflavone has the highest binding energy for SARS-CoV-2 papain-like protease. These molecules could bind near the SARS-CoV-2 papain-like protease crucial catalytic triad, ubiquitination and ISGylation residues: Trp106, Asn109, Cys111, Met208, Lys232, Pro247, Tyr268, Gln269, His272, Asp286 and Thr301. Because blocking the papain-like protease is an important strategy in fighting against viruses, these compounds might be promising candidates for therapeutic intervention against COVID-19.

Targeting SARS-CoV-2 viral proteases as a therapeutic strategy to treat COVID-19.

The 21st century has witnessed three outbreaks of coronavirus (CoVs) infections caused by severe acute respiratory syndrome (SARS)-CoV, Middle East respiratory syndrome (MERS)-CoV, and SARS-CoV-2. Coronavirus disease 2019 (COVID-19), caused by SARS-CoV-2, spreads rapidly and since the discovery of the first COVID-19 infection in December 2019, has caused 1.2 million deaths worldwide and 226,777 deaths in the United States alone. The high amino acid similarity between SARS-CoV and SARS-CoV-2 viral proteins supports testing therapeutic molecules that were designed to treat SARS infections during the 2003 epidemic. In this review, we provide information on possible COVID-19 treatment strategies that act via inhibition of the two essential proteins of the virus, 3C-like protease (3CLpro ) or papain-like protease (PLpro ).

Can molecular mimicry explain the cytokine storm of SARS-CoV-2?: An in silico approach.

PARP14 and PARP9 play a key role in macrophage immune regulation. SARS-CoV-2 is an emerging viral disease that triggers hyper-inflammation known as a cytokine storm. In this study, using in silico tools, we hypothesize about the immunological phenomena of molecular mimicry between SARS-CoV-2 Nsp3 and the human PARP14 and PARP9. The results showed an epitope of SARS-CoV-2 Nsp3 protein that contains consensus sequences for both human PARP14 and PARP9 that are antigens for MHC Classes 1 and 2, which can potentially induce an immune response against human PARP14 and PARP9; while its depletion causes a hyper-inflammatory state in SARS-CoV-2 patients.

Microscopic interactions between ivermectin and key human and viral proteins involved in SARS-CoV-2 infection.

The identification of chemical compounds able to bind specific sites of the human/viral proteins involved in the SARS-CoV-2 infection cycle is a prerequisite to design effective antiviral drugs. Here we conduct a molecular dynamics study with the aim to assess the interactions of ivermectin, an antiparasitic drug with broad-spectrum antiviral activity, with the human Angiotensin-Converting Enzyme 2 (ACE2), the viral 3CLpro and PLpro proteases, and the viral SARS Unique Domain (SUD). The drug/target interactions have been characterized in silico by describing the nature of the non-covalent interactions found and by measuring the extent of their time duration along the MD simulation. Results reveal that the ACE2 protein and the ACE2/RBD aggregates form the most persistent interactions with ivermectin, while the binding with the remaining viral proteins is more limited and unspecific.

Identification of potential plant-based inhibitor against viral proteases of SARS-CoV-2 through molecular docking, MM-PBSA binding energy calculations and molecular dynamics simulation.

The Coronavirus disease 2019 (COVID-19), caused by the novel coronavirus, SARS-CoV-2, has recently emerged as a pandemic. Here, an attempt has been made through in-silico high throughput screening to explore the antiviral compounds from traditionally used plants for antiviral treatments in India namely, Tea, Neem and Turmeric, as potential inhibitors of two widely studied viral proteases, main protease (Mpro) and papain-like protease (PLpro) of the SARS-CoV-2. Molecular docking study using BIOVIA Discovery Studio 2018 revealed, (-)-epicatechin-3-O-gallate (ECG), a tea polyphenol has a binding affinity toward both the selected receptors, with the lowest CDocker energy - 46.22 kcal mol-1 for SARS-CoV-2 Mpro and CDocker energy - 44.72 kcal mol-1 for SARS-CoV-2 PLpro, respectively. The SARS-CoV-2 Mpro complexed with (-)-epicatechin-3-O-gallate, which had shown the best binding affinity was subjected to molecular dynamics simulations to validate its binding affinity, during which, the root-mean-square-deviation values of SARS-CoV-2 Mpro-Co-crystal ligand (N3) and SARS-CoV-2 Mpro- (-)-epicatechin-3-O-gallate systems were found to be more stable than SARS-CoV-2 Mpro system. Further, (-)-epicatechin-3-O-gallate was subjected to QSAR analysis which predicted IC50 of 0.3281 nM against SARS-CoV-2 Mpro. Overall, (-)-epicatechin-3-O-gallate showed a potential binding affinity with SARS-CoV-2 Mpro and could be proposed as a potential natural compound for COVID-19 treatment.

Molecular docking studies, molecular dynamics and ADME/tox reveal therapeutic potentials of STOCK1N-69160 against papain-like protease of SARS-CoV-2.

SARS-CoV-2 is a new strain of Coronavirus that caused the pneumonia outbreak in Wuhan, China and has spread to over 200 countries of the world. It has received worldwide attention due to its virulence and high rate of infection. So far, several drugs have experimented against SARS-CoV-2, but the failure of these drugs to specifically interact with the viral protease necessitates urgent measure to boost up researches for the development of effective therapeutics against SARS-CoV-2. Papain-like protease (PLpro) of the viral polyproteins is essential for maturation and infectivity of the virus, making it one of the prime targets explored for SARS-CoV-2 drug design. This study was conducted to evaluate the efficacy of ~ 50,000 natural compounds retrieved from IBS database against COVID-19 PLpro using computer-aided drug design. Based on molecular dock scores, molecular interaction with active catalytic residues and molecular dynamics (MD) simulations studies, STOCK1N-69160 [(S)-2-((R)-4-((R)-2-amino-3-methylbutanamido)-3-(4-chlorophenyl) butanamido) propanoic acid hydrochloride] has been proposed as a novel inhibitor against COVID-19 PLpro. It demonstrated favourable docking score, the free energy of binding, interacted with key amino acid residues necessary for PLpro inhibition and also showed significant moderation for parameters investigated for ADME/tox (Adsorption, distribution, metabolism, excretion and toxicological) properties. The edge of the compound was further established by its stability in MD simulation conducted for 30 ns employing GROMACS software. We propose that STOCK1N-69160 is worth further investigation for preventing SARS-CoV-2.

Metal-Bound Methisazone; Novel Drugs Targeting Prophylaxis and Treatment of SARS-CoV-2, a Molecular Docking Study.

SARS-CoV-2 currently lacks effective first-line drug treatment. We present promising data from in silico docking studies of new Methisazone compounds (modified with calcium, Ca; iron, Fe; magnesium, Mg; manganese, Mn; or zinc, Zn) designed to bind more strongly to key proteins involved in replication of SARS-CoV-2. In this in silico molecular docking study, we investigated the inhibiting role of Methisazone and the modified drugs against SARS-CoV-2 proteins: ribonucleic acid (RNA)-dependent RNA polymerase (RdRp), spike protein, papain-like protease (PlPr), and main protease (MPro). We found that the highest binding interactions were found with the spike protein (6VYB), with the highest overall binding being observed with Mn-bound Methisazone at -8.3 kcal/mol, followed by Zn and Ca at -8.0 kcal/mol, and Fe and Mg at -7.9 kcal/mol. We also found that the metal-modified Methisazone had higher affinity for PlPr and MPro. In addition, we identified multiple binding pockets that could be singly or multiply occupied on all proteins tested. The best binding energy was with Mn-Methisazone versus spike protein, and the largest cumulative increases in binding energies were found with PlPr. We suggest that further studies are warranted to identify whether these compounds may be effective for treatment and/or prophylaxis.