Generalized open-source workflows for atomistic molecular dynamics simulations of viral helicases.

Viral helicases are promising targets for the development of antiviral therapies. Given their vital function of unwinding double-stranded nucleic acids, inhibiting them blocks the viral replication cycle. Previous studies have elucidated key structural details of these helicases, including the location of substrate binding sites, flexible domains, and the discovery of potential inhibitors. Here we present a series of new Galaxy tools and workflows for performing and analyzing molecular dynamics simulations of viral helicases. We first validate them by demonstrating recapitulation of data from previous simulations of Zika (NS3) and SARS-CoV-2 (NSP13) helicases in apo and complex with inhibitors. We further demonstrate the utility and generalizability of these Galaxy workflows by applying them to new cases, proving their usefulness as a widely accessible method for exploring antiviral activity.

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

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

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.

What Strengthens Protein-Protein Interactions: Analysis and Applications of Residue Correlation Networks.

Identifying residues critical to protein-protein binding and efficient design of stable and specific protein binders are challenging tasks. Extending beyond the direct contacts in a protein-protein binding interface, our study employs computational modeling to reveal the essential network of residue interactions and dihedral angle correlations critical in protein-protein recognition. We hypothesized that mutating residues exhibiting highly correlated dynamic motion within the interaction network could efficiently optimize protein-protein interactions to create tight and selective protein binders. We tested this hypothesis using the ubiquitin (Ub) and MERS coronaviral papain-like protease (PLpro) complex, since Ub is a central player in multiple cellular functions and PLpro is an antiviral drug target. Our designed ubiquitin variant (UbV) hosting three mutated residues displayed a ∼3,500-fold increase in functional inhibition relative to wild-type Ub. Further optimization of two C-terminal residues within the Ub network resulted in a KD of 1.5 nM and IC50 of 9.7 nM for the five-point Ub mutant, eliciting 27,500-fold and 5,500-fold enhancements in affinity and potency, respectively, as well as improved selectivity, without destabilizing the UbV structure. Our study highlights residue correlation and interaction networks in protein-protein interactions, and introduces an effective approach to design high-affinity protein binders for cell biology research and future therapeutics.

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.

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.

Two Binding Sites of SARS-CoV-2 Macrodomain 3 Probed by Oxaprozin and Meclomen.

Severe acute respiratory syndrome-coronavirus-1/2 (SARS-CoV-1/2) macrodomain 3 (Mac3) is critical for replication and transcription of the viral genome and is therefore a potential therapeutic target. Here, we solved the crystal structure of SARS-CoV-2 Mac3, which reveals a small-molecule binding pocket. Two low-molecular-weight drugs, oxaprozin and meclomen, induced different patterns of nuclear magnetic resonance (NMR) chemical shift perturbations (CSPs). Meclomen binds to site I of SARS-CoV-2 Mac3 with binding pose determined by NMR CSP and transferred paramagnetic relaxation enhancement, while oxaprozin binds to site II as revealed by the crystal structure. Interestingly, oxaprozin and meclomen both perturb residues in site I of SARS-CoV Mac3. Fluorescence polarization experiments further demonstrated that oxaprozin and meclomen inhibited the binding of DNA-G4s to SARS-CoV-2 Mac3. Our work identified two adjacent ligand-binding sites of SARS-CoV-2 Mac3 that shall facilitate structure-guided fragment linking of these compounds for more potent inhibitors.

Insights into Coronavirus Papain-like Protease Structure, Function and Inhibitors.

The coronavirus family consists of pathogens that seriously affect human and animal health. They mostly cause respiratory or enteric diseases, which can be severe and life-threatening, such as coronavirus disease 2019 (COVID-19), severe acute respiratory syndrome (SARS), and Middle East Respiratory Syndrome (MERS) in humans. The conserved coronaviral papain-like protease is an attractive antiviral drug target because it is essential for coronaviral replication, and it also inhibits host innate immune responses. This review focuses on the latest research progress relating to the mechanism of coronavirus infection, the structural and functional characteristics of coronavirus papain-like protease, and the potent inhibitors of the protease.

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.

Identification of SARS-CoV-2 Papain-like Protease (PLpro) Inhibitors Using Combined Computational Approach.

In the current pandemic, finding an effective drug to prevent or treat the infection is the highest priority. A rapid and safe approach to counteract COVID-19 is in silico drug repurposing. The SARS-CoV-2 PLpro promotes viral replication and modulates the host immune system, resulting in inhibition of the host antiviral innate immune response, and therefore is an attractive drug target. In this study, we used a combined in silico virtual screening for candidates for SARS-CoV-2 PLpro protease inhibitors. We used the Informational spectrum method applied for Small Molecules for searching the Drugbank database followed by molecular docking. After in silico screening of drug space, we identified 44 drugs as potential SARS-CoV-2 PLpro inhibitors that we propose for further experimental testing.

Design and Evaluation of a Novel Peptide-Drug Conjugate Covalently Targeting SARS-CoV-2 Papain-like Protease.

Coronavirus disease 2019 (COVID-19) pandemic, a global health threat, was caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The SARS-CoV-2 papain-like cysteine protease (PLpro) was recognized as a promising drug target because of multiple functions in virus maturation and antiviral immune responses. Inhibitor GRL0617 occupied the interferon-stimulated gene 15 (ISG15) C-terminus-binding pocket and showed an effective antiviral inhibition. Here, we described a novel peptide-drug conjugate (PDC), in which GRL0617 was linked to a sulfonium-tethered peptide derived from PLpro-specific substrate LRGG. The EM-C and EC-M PDCs showed a promising in vitro IC50 of 7.40 ± 0.37 and 8.63 ± 0.55 µM, respectively. EC-M could covalently label PLpro active site C111 and display anti-ISGylation activities in cellular assays. The results represent the first attempt to design PDCs composed of stabilized peptide inhibitors and GRL0617 to inhibit PLpro. These novel PDCs provide promising opportunities for antiviral drug design.

Computational investigation of binding of chloroquinone and hydroxychloroquinone against PLPro of SARS-CoV-2.

Novel coronavirus SARS-CoV-2 has infected 18 million people with 700,000+ mortalities worldwide and this deadly numeric figure is rapidly rising. With very few success stories, the therapeutic targeting of this epidemic has been mainly attributed to main protease (Mpro), whilst Papain-like proteases (PLpro) also plays a vital role in the processing of replicase polyprotein. Multifunctional roles of PLpro such as viral polypeptide cleavage, de-ISGlyation and immune suppression have made it a promising drug target for therapeutic interventions. Whilst there have been a number of studies and others are on-going on repurposing and new-small molecule screening, albeit previously FDA approved drugs viz. Chloroquine (CQ) and Hydroxychloroquine (HCQ) have only been found effective against this pandemic. Inspired by this fact, we have carried out molecular docking and dynamics simulation studies of FDA approved CQ and HCQ against SARS-CoV-2 PLpro. The end aim is to characterise the binding mode of CQ and HCQ and identify the key amino acid residues involved in the mechanism of action. Further, molecular dynamics simulations (MDS) were carried out with the docked complex to search for the conformational space and for understanding the integrity of binding mode. We showed that the CQ and HCQ can bind with better binding affinity with PLpro as compared to reference known PLpro inhibitor. Based on the presented findings, it can be anticipated that the SARS-CoV-2 PLpro may act as molecular target of CQ and HCQ, and can be projected for further exploration to design potent inhibitors of SARS-CoV-2 PLpro in the near future.

Evaluating Stability and Activity of SARS-CoV-2 PLpro for High-throughput Screening of Inhibitors.

Because of the essential roles of SARS-CoV-2 papain-like protease (PLpro) in the viral polyprotein processing and suppression of host immune responses, it is a crucial target for drug discovery against COVID-19. To develop robust biochemical methodologies for inhibitor screening against PLpro, extensive characterization of recombinant protein is important. Here we report cloning, expression, and purification of the recombinant SARS-CoV-2 PLpro, and explore various parameters affecting its stability and the catalytic activity. We also report the optimum conditions which should be used for high-throughput inhibitor screening using a fluorogenic tetrapeptide substrate.

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.

Carbon fullerene and nanotube are probable binders to multiple targets of SARS-CoV-2: Insights from computational modeling and molecular dynamic simulation studies.

The present study aimed to predict the binding potential of carbon nanotube and nano fullerene towards multiple targets of SARS-CoV-2. Based on the virulent functions, the spike glycoprotein, RNA-dependent RNA polymerase, main protease, papain-like protease, and RNA binding domain of the nucleocapsid proteins of SARS-CoV-2 were prioritized as the molecular targets and their three-dimensional (3D) structures were retrieved from the Protein Data Bank. The 3D structures of carbon nanotubes and nano-fullerene were computationally modeled, and the binding potential of these nanoparticles to the selected molecular targets was predicted by molecular docking and molecular dynamic (MD) simulations. The drug-likeness and pharmacokinetic features of the lead molecules were computationally predicted. The current study suggested that carbon fullerene and nanotube demonstrated significant binding towards the prioritized multi-targets of SARS-CoV-2. Interestingly, carbon nanotube showed better interaction with these targets when compared to carbon fullerene. MD simulation studies clearly showed that the interaction of nanoparticles and selected targets possessed stability and conformational changes. This study revealed that carbon nanotubes and fullerene are probably used as effectual binders to multiple targets of SARS-CoV-2, and the study offers insights into the experimental validation and highlights the relevance of utilizing carbon nanomaterials as a therapeutic remedy against COVID-19.

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.

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.

Destruction mechanisms of ozone over SARS-CoV-2.

In this pandemic SARS-CoV-2 crisis, any attempt to contain and eliminate the virus will also stop its spread and consequently decrease the risk of severe illness and death. While ozone treatment has been suggested as an effective disinfection process, no precise mechanism of action has been previously reported. This study aimed to further investigate the effect of ozone treatment on SARS-CoV-2. Therefore, virus collected from nasopharyngeal and oropharyngeal swab and sputum samples from symptomatic patients was exposed to ozone for different exposure times. The virus morphology and structure were monitored and analyzed through Atomic Force Microscopy (AFM), Transmission Electron Microscopy (TEM), Atomic Absorption Spectroscopy (AAS), and ATR-FTIR. The obtained results showed that ozone treatment not only unsettles the virus morphology but also alters the virus proteins' structure and conformation through amino acid disturbance and Zn ion release from the virus non-structural proteins. These results could provide a clearer pathway for virus elimination and therapeutics preparation.

Hesperidin Is a Potential Inhibitor against SARS-CoV-2 Infection.

Hesperidin (HD) is a common flavanone glycoside isolated from citrus fruits and possesses great potential for cardiovascular protection. Hesperetin (HT) is an aglycone metabolite of HD with high bioavailability. Through the docking simulation, HD and HT have shown their potential to bind to two cellular proteins: transmembrane serine protease 2 (TMPRSS2) and angiotensin-converting enzyme 2 (ACE2), which are required for the cellular entry of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Our results further found that HT and HD suppressed the infection of VeroE6 cells using lentiviral-based pseudo-particles with wild types and variants of SARS-CoV-2 with spike (S) proteins, by blocking the interaction between the S protein and cellular receptor ACE2 and reducing ACE2 and TMPRSS2 expression. In summary, hesperidin is a potential TMPRSS2 inhibitor for the reduction of the SARS-CoV-2 infection.

Molecular modelling on SARS-CoV-2 papain-like protease: an integrated study with homology modelling, molecular docking, and molecular dynamics simulations.

SARS-CoV-2 PLpro was investigated as a therapeutic target for potent antiviral drugs due to its essential role in not only viral replication but also in regulating the inborn immune response. Several computational approaches, including homology modelling, molecular docking, and molecular dynamics (MD) studies, were employed to search for promising drugs in treating SARS-CoV-2. Eighty-one compounds, sub-structurally similar to the antiviral drug, were used as potential inhibitors of PLpro. From our results, three complexes containing the ligands with Pubchem IDs: 153012995, 12149203, and 123608715 showed lower binding energies than the control (Ritonavir), indicating that they may become promising inhibitors for PLpro. MD was performed in a water solvent to validate the stability of the three complexes. All complexes achieved stable structure during the simulation as no significant fluctuations were observed in the validation parameters. Moreover, the binding energy for each complex was estimated using the MM-GBSA method. Complex 1 was the most stable structure based on the lowest binding energy score and its structure remained in a similar cavity with the docket snapshot. Based on our studies, three ligands were assumed to be potential inhibitors. The ligand of complex 1 may become the most promising antiviral drug against SARS-CoV-2 targeting PLpro.