key: cord-0883808-alkvt58v
authors: Amin, Sk. Abdul; Ghosh, Kalyan; Singh, Samayaditya; Qureshi, Insaf Ahmed; Jha, Tarun; Gayen, Shovanlal
title: Exploring naphthyl derivatives as SARS-CoV papain-like protease (PLpro) inhibitors and its implications in COVID-19 drug discovery
date: 2021-03-06
journal: Mol Divers
DOI: 10.1007/s11030-021-10198-3
sha: 15175434edbbe173346c48692a8ce1f64a885ed3
doc_id: 883808
cord_uid: alkvt58v

ABSTRACT: Novel coronavirus disease 2019 (COVID-19) emerges as a serious threat to public health globally. The rapid spreading of COVID-19, caused by severe acute respiratory syndrome (SARS) coronavirus 2 (SARS-CoV-2), proclaimed the multitude of applied research needed not only to save the human health but also for the environmental safety. As per the recent World Health Organization reports, the novel corona virus may never be wiped out completely from the world. In this connection, the inhibitors already designed against different targets of previous human coronavirus (HCoV) infections will be a great starting point for further optimization. Pinpointing biochemical events censorious to the HCoV lifecycle has provided two proteases: a papain-like protease (PLpro) and a 3C-like protease (3CLpro) enzyme essential for viral replication. In this study, naphthyl derivatives inhibiting PLpro enzyme were subjected to robust molecular modelling approaches to understand different structural fingerprints important for the inhibition. Here, we cover two main aspects such as (a) exploration of naphthyl derivatives by classification QSAR analyses to find important fingerprints that module the SARS-CoV PLpro inhibition and (b) implications of naphthyl derivatives against SARS-CoV-2 PLpro enzyme through detailed ligand–receptor interaction analysis. The modelling insights will help in the speedy design of potent broad spectrum PLpro inhibitors against infectious SARS-CoV and SARS-CoV-2 in the future. GRAPHICAL ABSTRACT: [Image: see text] SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at(10.1007/s11030-021-10198-3) .

It is just alarming and ever surprising how rapidly a communicable novel coronavius (2019-nCoV), cause of novel coronavirus disease 2019 (COVID- 19) , has been spreading in the world in the twenty-first century [1] [2] [3] . 2019-nCoV is also known as severe acute respiratory syndrome (SARS) coronavirus 2 (SARS-CoV-2). World Health Organization (WHO) characterized COVID-19 as world pandemic [4] . So far, more than million deaths have been reported from 216 countries and territories [5] . Every death makes us painfully aware that our swords are blunt till date in the battle against this hazardous COVID-19.

Researchers around the world have been trying different options to restrict the virus replications [6] [7] [8] [9] [10] [11] [12] [13] [14] [15] . However, there is still no effective drug/vaccine against this virus. In this situation, emphasis should also be given to the systematic rational drug discovery against different targets of the virus. Among the different targets, two proteases namely papain-like protease (PLpro) and a 3C-like protease (3CLpro) are very crucial for virus replication and are considered as important druggable targets [3, 6, 7, [16] [17] [18] [19] [20] [21] [22] . The PLpro enzyme also shows deubiquitinating (DUB) and deIS-Gylating activities [9, 16] . As a result, it is also responsible for host cell immune suppression due to the inactivation of NF-κB pathway (Fig. 1 ). In addition, the structures of different PLpro enzymes are very similar in different coronavirus, and therefore, it is considered as a target for broad-spectrum inhibitor development.

Naphthalene is the most straightforward member of the class of PLpro inhibitors [23] [24] [25] [26] , in which a couple of benzene rings are fused in the ortho positions. Numerous naphthalene-containing molecules have also been reported to boast significant antimicrobial property. A commonly used dye, β-naphthol, exhibits antimicrobial activity [27] . In addition, naphthyl-based drugs including naftifine, terbinafine, nafacillin, tolnaftate, etc., are found to possess antimicrobial property [27] [28] [29] . Ratia and collaborators first introduced naphthyl derivatives those were likely to act as non-covalent competitive inhibitors of PLpro [23] . Naphthyl derivative binds within the S4-S3 subsites of the enzyme, thereby inducing a loop closure which ultimately results in conformational change and manifests the PLpro active site as non-functional. Several other naphthyl derivatives were further reported as PLpro inhibitors [23] [24] [25] [26] . In recent times, the interest on naphthyl derivatives is tremendously increased as these derivatives have shown potential SARS-CoV-2 PLpro inhibition [30] .

Thus, in this work, we have focused our attention on naphthyl derivatives as SARS-CoV-2 PLpro inhibitors. Multiple modelling strategies were applied with these mottos: (a) identification of important fingerprints that modulate the SARS-CoV PLpro inhibition and (b) scope of naphthyl derivatives to target SARS-CoV-2 PLpro though ligand-receptor interaction analysis. The current study, a part of our rational antiviral drug design and discovery Sk. Abdul Amin and Kalyan Ghosh have equal contribution.

* Tarun Jha tjupharm@yahoo.com * Shovanlal Gayen shovanlal.gayen@gmail.com programme [31] [32] [33] [34] , may offer an initiative to explore the possibility of broad spectrum inhibitors against the PLpro enzyme in case of both SARS-CoV and SARS-CoV-2.

In order to accelerate the drug discovery effort against coronavirus, a set of diverse naphthyl derivatives belonging to a collection of SARS-CoV PLpro inhibitors were collected with inhibitory activities (IC 50 ) [23] [24] [25] [26] . The naphthyl derivatives with no activity and without definite activity were eliminated initially and followed by removal of duplicates from the study. Thus, remaining fifty-six molecules were considered for the further study (Table 1) .

The classification modelling assists to classify the active and inactive molecules in terms of their biological data. Here, we have employed structural and physico-chemical interpretation (SPCI) [35, 36] and Monte Carlo-based coral QSAR [37] [38] [39] [40] [41] studies. For these studies, the SARS-CoV PLpro pIC 50 value, 5.3, was considered as the threshold value. Compounds having the SARS-CoV PLpro pIC 50 value of 5.3 or more were classified as higher PLpro inhibitors or actives and those with less PLpro pIC 50 than threshold value were distinguished as inactives.

Structural and physico-chemical interpretation (SPCI) study [35, 36] was solicited to identify and estimate contributions of scaffolds and/or linkers and/or single substituent to these naphthyl-based SARS-CoV PLpro inhibitors. At first, descriptors were calculated by the aid of SiRMS tool followed by model development and validation. Four different classification QSAR models were undertaken by using machine learning approaches including gradient boosting classification (GBC), random forest (RF), support vector machine (SVM), and k-nearest neighbour (kNN). The developed models were evaluated by statistical parameters such as balanced accuracy, sensitivity, and specificity [36] . Furthermore, molecular fragmentation of the dataset was done to estimate the contribution from the developed models. Fragments consist of at most three attachment points were preferred, and subsequently, preferred fragments were counted by RDKit in combination with SMARTS pattern [! #1]! @! = ! # [! #1]. Lastly, the overall contribution of the different fragments and their median fragment contribution graphs were generated by using rspciR software package [42] .

Monte Carlo optimization method was used to identify the important structural fingerprints that are exclusively responsible for promoting or hindering of activity [37] [38] [39] [40] [41] . Here, the molecular structures of the different inhibitors were represented by SMILES (Simplified Molecular Input Line System) format. These symbolic notations were used to represent structural attributes such as atoms, bonds, etc., and to calculate In this method, the models were generated on the basis of three optimal descriptors like: SMILES, graph and hybrid. The SMILES-based optimal descriptors used in the study were calculated as:

where T represents threshold and N is number of epoch used for generation of models. Further, the constants like: a, b, c, d, α, β and γ represents different coefficients that are used to modify the descriptors. The CW stands for correlation weight of specific structural attributes for PLpro inhibition.

Global SMILE attributes include NOSP (presence or absence of nitrogen, oxygen, sulphur and phosphorus atom), HALO (presence or absence of halogen groups like: fluorine, chlorine, bromine and iodine), BOND (presence or absence of bond like: = is double bond, # is triple bond and @ is stereo-chemical bond), and ATOMPAIR (represents presence The graph-based descriptors comprise of three types they are: GAO (graph of atomic orbital), HSG (hydrogensuppressed graph) and HFG (hydrogen-filled graph). These graph-based descriptors are calculated by means of different molecular connectivity indices like: 0 EC k , 1 EC k and 3 EC k . The graph-based descriptors are represented as:

Here, the chemical atoms like C, N, O, etc., are represented by A k . The coefficients having value 0 and 1 are denoted by as α, β and γ. The coefficients having value 1 are generally used for model building, whereas 0 value are excluded. Further, the notations like: 0 EC k , 1 EC k and 3 EC k denotes different Morgan's connectivity indices used in the [41] . Thus, the combination of SMILES and graphbased descriptors forms hybrid descriptors which are represented as:

In our study, twenty-one robust classification models were generated from three different splits on the basis of optimum descriptors (SMILES, graph and hybrid) using the balance of correlation method. The dataset consisting of 56 compounds was distributed into the sub-training (25 compounds), calibration (21 compounds) and test (10 compounds). The training set compounds were used for model building, calibration set was used to check the predictive potential and to prevent overtraining of the model, and the test set is used as an estimator to validate the models developed using training and calibration set [43, 44] 

The molecular docking experiments were implemented with the help of AutoDock Vina [45] . Prior to docking study, polar hydrogen atoms were added by the aid of AutoDock Tools (ADT) [46] in order to relax the conformational strain. The docked poses of naphthyl derivatives were visualized by using the Discovery Studio 3.5 Visualizer [47] .

The MD simulations for apo and each complex form of the protein were conducted by the GROMACS v5.1.4 [48] using GROMOS96 45a3 force field [49] for the parameterization of the protein. The ligand topology was generated through GROMOS96 force field by the PRODRG 2 server [50]. Each Fig. 2 Schematic representation of current work design of the system was solved by the SPC/E water molecules, and the charge on simulating systems was neutralized by providing the equal number of counter ions [51] . Then, minimization of the systems was executed by the steepest descent algorithm in 50,000 steps. After the minimization, each system was subjected to the position restrained NVT (300 K) and NPT (1 atm) simulations of 100 ps for the equilibration according to the mentioned temperature and pressure. Equilibrated systems were then submitted for the non-position restrained production simulations of 20 ns each, and the generated data of MD were then used for the enumeration of the RMSD, RMSF and Rg data of the protein backbone, which helps in structural analysis during dynamics. The binding energy of each complex was calculated by g_MMPBSA package of the GROMACS considering 200 frames from 20 ns data of simulations [52] .

In this twenty-first century, health crisis posted by SARS-CoV-2 outbreak, drug repurposing and/or screening of SARS-CoV inhibitor databases is the fastest option in terms of strategic and economic way. However, systematic and rational drug discovery approaches to find potent inhibitors against different targets of SARS-CoV-2 should not be ignored. In this connection, the inhibitors already designed against SARS-CoV will be very helpful due to the sequence similarity between them. Different modelling approaches were applied on naphthyl derivatives as SARS-CoV PLpro inhibitors as outlined in Fig. 2 .

Performing classification-based structural and physicochemical interpretation (SPCI) and Monte Carlo-based coral QSAR methods enable interesting visualization of important chemical sub-structural features attributed to enhance or decrease the SARS-CoV PLpro inhibitory properties. Considering the potency threshold value, out of 56 compounds, 31 compounds were identified as lower and 25 compounds were denoted as higher SARS-CoV PLpro inhibitors.

In our study, four different models like gradient boosting machine (GBM), random forest (RF), support vector machine (SVM) and k-nearest neighbour (kNN) were developed. The different tuning parameters used to build the machine learning models for naphthyl-based PLpro inhibitors are given in Table S1 . The models generated by using SPCI analysis were found to have acceptable statistical parameters as shown in Table 2 .

Among four different models (GBM, RF, SVM and kNN), GBM and RF models are found to have similar statistical parameters. Further, a consensus model was generated in order to remove the biasness of the individual model; this was further considered for the interpretation. The contributions of different fragments obtained from the SPCI analysis are shown in Fig. 3 .

A limpid trend of structure-activity (SARS-CoV PLpro inhibition) relationship is recognized in Fig. 3 . Piperidine moiety, acyclic CONH, aromatic NH 2 , aromatic halogen and (Fig. 3) .

In Monte Carlo optimization, a total of twenty-one different models from three different splits were generated using SMILES and graph-based descriptors with a combination of different connectivity indices calculated for generation of different models (Table S2) . Each models were developed after the search for desirable T (threshold) and N (epoch) values. Among these models (Table S2) , the model M3 (SMILES and GAO with 1 EC k connectivity of split-1) was found to be the best model and was used for the physical interpretation of structural fingerprints for PLpro inhibition ( Table 3 ). The best model from each split is shown in bold face. The end point values calculated for M3 are shown below:

Different structural attributes of the best model M3 obtained from Monte Carlo optimization method are depicted in Table S3 .

Each compound of this dataset can be represented as consisting of three parts: head and tail connected by a linker moiety (Fig. 4) .

As the head region is structurally conserved in between 1 and 2-naphthyl, maximum influenced sub-structure/ Endpoint = 0.1381188 (± 0.0113900) + 0.0308263 (± 0.0006586) * DCW (4, 5) fingerprints were noticed from linker and tail portions. More interestingly, a limpid trend of structure-activity (SARS-CoV PLpro inhibition) relationship was recognized in Fig. 3 . From the SAR study, it can be conferred that the 1-naphthyl head was preferable over the 2-naphthyl prototypes. Thus, the polarizability of the 1-naphthyl function modulates the binding against SARS-CoV PLpro (Fig. 5a) . The importance of a single methyl substituent at the R 3 position when compared to unsubstituted compounds 43-45, 55-56. Notably, the stereo-chemical pattern of the methyl substituent was a critical factor to modulate PLpro binding affinity. This can be understood by observing the X-ray crystal structure of the compound 53-bounded SARS-CoV PLpro complex (Fig. 5b) where the (R)-methyl enantiomer extends into an interior of the PLpro enzyme between Tyr-265 and Thr-302. 24 Similarly, (R)-1-naphthylethylamide of compound 9 interacted with Tyr-265 and Tyr-269 aromatic rings and with side chains of Pro-248 and Pro-249 and consequently forms hydrophobic interactions. 23 Báez-Santos and collaborators 24 very nicely explained the loss of PLpro inhibitory activity with proportional to substituent size at R 3 position, due to the potential entropic gain by displacing the water molecules lead to a larger enthalpic penalty of breaking the H-bonds (between these water molecules and aminoacid residues D165, R167, Y274, T302 and D303). Thus, in all docking calculations as well as reported crystal structures with different inhibitors (PDB: 4OW0, 3MJ5, 4OVJ), the conserved water molecules were present in the binding site for ligand-receptor interactions. Piperidine moiety, acyclic CONH, NH 2 , aromatic halogen and phenyl unsubstituted functions exhibited positive effects on SARS-CoV PLpro inhibitory activities (Fig. 5c-f) . However, the methoxy group attached with an aromatic ring suggested clear negative influence on biological properties. Figure 3 predicted the positive contribution of acyclic CONH and aromatic NH 2 in the PLpro inhibition. This can be justified by the SAR observation. At the middle position, the amide NH was demonstrated to be an optimal feature for PLpro inhibition. Addition of a methyl group at amide nitrogen resulted in significantly differing activity levels (the N-methyl derivative 67, IC 50 = 22.6 μM vs compound 6, IC 50 = 2.3 μM). In addition, the aromatic NH 2 was found to be conducive for compounds 9 (IC 50 = 0.6 μM), (Table S3) . Table S3 demonstrated the importance of 'C…(…1…' that can be decoded as sp 3 carbons with branching attached to a ring. This is explained by the presence of methylamine fragments, another promoter of activity which contributes to increase of biological activity against SARS-CoV PLpro. This can be explained by compound 77 where incorporation of 5-methylamine substituent on the benzamide ring would result in escalation of the inhibitory activity (IC 50 = 0.46 μM). Thus, addition of donor groups (− NH 2 ) at the benzamide ring can lead to promising PLpro inhibitory activity (Fig. 5f) . Further introduction of a methyl function to the amine group of 77 yielded compound 78 with minutely narrowed PLpro inhibitory activity (IC 50 = 1.3 μM). However, compound 78 exhibits slight enhancement of anti-SARS-CoV potency (EC 50 = 5.2 μM) compared to compound 77 (EC 50 = 6 μM). Compounds 65 and 68 were poor PLpro inhibitors, both comprising p-NH 2 at R 2 position. Therefore, it was evidenced that p-NH 2 alone at R 2 position might hamper the binding with the PLpro enzyme.

The presence of 'O…(…1…' (as evidenced by the Monte Carlo-based QSAR analysis) can be decoded by comparing compounds 43-45, 50, 61 and 63 where OCH 3 substitution negatively influenced the PLpro inhibitory activity (IC 50 range 30 to 149 μM). Thus, increasing the level of negatively charged atom with hydrophobic characteristics of methoxy function at the terminal phenyl ring would result in fall of the biological activity against PLpro. This observation can be understood by the molecular interaction that substituents in 4-position of tail phenyl ring should be from a positively charged group with additional demand of hydrophilic effects.

As mentioned previously, Fig. 3 showed the contributions of the unsubstituted phenyl function, positively, to the SARS-CoV PLpro inhibitory activity (Fig. 5c) . This observation may clearly be explained by comparing the naphthyl derivatives 18-22 (bearing unsubstituted phenyl tail), whose IC 50 values were found in between the range of 1.9 and 18. However, the marked differences in the PLpro inhibitory activities of these compounds were due to the type of R 4 substituents. Moreover, the effect of stereochemistry of R 4 substituents 'C…@…….' (as evidenced by the Monte Carlo-based QSAR analysis) played critically (compounds 21 vs 22). A tenfold decrease in PLpro inhibitory activity was manifested for (R)-methoxymethyl containing 21 (IC 50 = 18 μM), compared to the corresponding S-isomer 22 (IC 50 = 1.9 μM). Overall, the minute observation of the SAR study of these naphthyl derivatives 18-22 suggested that R 4 substituents were not at all pleasurable towards biological activity because the most active analogue in this series, 22, was similar effective the R 4 unsubstituted prototype 18 (IC 50 = 2.2 μM).

More interestingly, methylcarboxamide (compounds 24-25) and more electronegative halogen (compounds 32-33) substitutions at the 3rd or 4th position of the terminal phenyl ring manifested improvement in potency compare to the corresponding unsubstituted analogue 18. Figure 3 clearly expressed the importance of aromatic halogen substitutions.

Additionally, dioxolane derivatives 53-54 and 4-ethyl prototype 23 (IC 50 = 0.47 μM) also showed promising inhibitory activity against PLpro. The dioxolane group contributed positively to the PLpro inhibitory property not because of the interaction with Gln270 only; there should be another effect which contributed largely to the substituted phenyl derivatives to improve the biological potency. As noted, the mono-fluoro substitution at the phenyl ring induced conspicuous polarization effects in the π-system of the associated terminal phenyl ring and consequently, showed improved binding affinity with the PLpro active site amino acids residues.

Notably 

PLpro inhibitors have the potential to be broad spectrum inhibitors due to high sequence similarities of PLpro enzyme in different CoVs. It is interesting to know whether these naphthyl derivatives are effective also in SARS-CoV-2 PLpro target along with SARS-CoV PLpro. We have done both molecular docking and molecular dynamics simulations of three naphthyl derivatives, compounds 23, 27 and 32, against SARS-CoV-2 PLpro target (PDB: 6WUU). The energy minimized geometry of naphthyl derivatives were considered for the docking experiments against SARS-CoV-2 PLpro by the aid of AutoDock Vina [45] . The docking results are depicted as Fig. 6 . Figure 6 highlighted the superimposition of docking poses of compounds 23, 27 and 32 in SARS-CoV-2 PLpro (PDB: 6WUU). The docking analysis clearly highlights that naphthyl derivatives can nicely bind with the binding pocket of SARS-CoV-2 PLpro.

After the selection of the lowest energy conformation from the docking output of compounds 23, 27 and 32 (Fig. 7a-c) , the stability of the protein structure with those ligand conformations was analysed by MD simulations. The calculated average RMSD of the apo, prt-23, prt-27, and prt-32 are 0.320, 0.315, 0.309, and 0.357 nm, respectively. Data of the average RMSD showed that each of the complex protein structure has almost similar structural deviations like the apo form during dynamics (Fig. 7d) . But as depicted in the Fig. 7d , the backbone deviation of protein structure is more stable with the compound 23 in comparison of other compounds. Average RMSF in the residues of the apo, prt-23, prt-27 and prt-32 are calculated to be 0.142, 0.132, 0.133 and 0.143 nm during the dynamics, respectively (Fig. 7e) . These revealed that the fluctuations in the residues of protein are more relaxed after binding of compounds 23 and 27 as compared to compound 32 that presented fluctuation similar to apo protein. To analyse the induced changes in the compactness of protein structure after binding of ligand, the Rg data of each complex is compared with that of the apo form (Fig. 7f) . The average Rg in the backbone structure of the apo, prt-23, prt-27, and prt-32, are enumerated as 2.338, 2.327, 2.330, and 2.305 nm, respectively, which showed that the compactness of the structure in the apo and complexes in similar during the dynamics. The plotted data (Fig. 7f) depicted that the compactness of the protein backbone is slightly better in the presence of compound 32 in comparison to other compounds and apo protein.

Hence, comparative MD data analysis delineated that the protein bound with compounds has retained apo-like backbone deviations and compactness along with low fluctuations in the residues throughout the dynamics. Additionally, stability of each ligand conformation in the active site of the PLpro substantiates the docking studies.

The g_MMPBSA tool of the GROMACS is also used to determine the affinity of the compounds with protein as well as the free energy terms responsible for the affinity ( Table 4) .

The binding energy analysis showed that the compound 32 has more affinity towards protein during the dynamics in comparison with other compounds. Throughout the simulation, both of the electrostatic and van der Waals energy terms contributed majorly for the affinity between the compound and protein in the complexes.

This above analysis proves that naphthyl derivatives have potential to be inhibitor against SARS-CoV-2 PLpro enzyme. All the structural attributes identified by our different modelling approaches may be valid for SARS-CoV-2 PLpro enzyme also. Therefore, it can be emphasized that the naphthyl derivatives have potential to use as a seed for ligand design as well as optimization against SARS-CoV-2 PLpro enzyme by taking different modelling insights performed in this study.

Human coronavirus infections had almost been forgotten, and it was not challenged until novel coronavirus outbreak in December 2019. As per the recent World Health Organization reports, the novel corona virus may never be wiped out completely from the world. There are various strategies including drug repurposing, vaccine and immunity approaches, etc., employed to quickly react in the situation. However, scientific community should also give serious thought to start finding inhibitors against the different targets of the virus through rational drug discovery approaches. In this connection, the inhibitors already designed against different targets of previous human coronavirus infections will be a great starting point for further optimization. In spite of few drugs being assessed clinically against COVID-19, there remains thirst for discovering new molecules with increased efficacy as well as safety. Our research unit previously demonstrated the quantitative structure-activity relationship studies on SARS-CoV protease inhibitors [31] [32] [33] [34] .

Here, we endorsed rational drug design efforts through computational drug discovery approaches including machine learning and molecular docking studies. The different molecular modelling techniques such as structural and physico-chemical interpretation (SPCI) analysis and Monte Carlo optimization-based QSAR study collectively deliver some crucial structural information modulating SARS-CoV PLpro inhibitory activities. By considering all these QSAR models, molecular docking and MD simulation studies, it can be concluded that:

(1) Presence of 1-naphthyl head affects the activity, since it modulates the interactions at the active site residues (Fig. 8 ). This conclusion is based on SPCI analysis, Monte Carlo optimization-based QSAR and docking studies.

(2) The presence of a piperidine moiety is important for binding interaction, since it plays a significant role in interaction with active site residues as suggested by the molecular docking study (Fig. 8) .

(3) The stereo-chemical pattern of the methyl substituent at R 3 position is a critical factor to modulate PLpro binding affinity (Fig. 8) . This can be understood by the docking study. Similarly, (R)-1-naphthylethylamide function is important since it interacts with amino acid residues Tyr-265 and Tyr-269 aromatic rings and with side chain of amino acid residues Pro-248 and Pro-249 and consequently, forming hydrophobic interactions.

(4) The unsubstituted R 4 position is more favourable for PLpro inhibition.

(5) The amino function at the linker area is important to form hydrogen bond interaction with PLpro active site residue as evidence by the molecular docking study.

(6) Presence of p-NH 2 alone at the R 2 position might hamper the binding with the PLpro enzyme.

(7) 5-Methylamine and halogen substituent at the R 2 position are beneficial for the PLpro inhibitory activity. (8) Methoxy substation at the R 2 position is unfavourable for biological activity towards PLpro enzyme. This conclusion is based on SPCI analysis and Monte Carlo optimization-based QSAR studies (Fig. 8) .

Significantly, the above requisite structural features (Fig. 8 ) enlighten the perspective of medicinal chemists to develop potent PLpro inhibitors in the future. Research funding agencies and industries should consider in vitro and in vivo studies of the investigated naphthyl derivatives as a seed which sustain significant hope against SARS-CoV-2. 

The global spread of 2019-nCoV: a molecular evolutionary analysis

Therapeutics for COVID-19: from computation to practices-where we are, where we are heading to

Fight against novel coronavirus: a perspective of medicinal chemists

Mesecar AD (2020) Drug development and medicinal chemistry efforts toward SARS-Coronavirus and Covid-19 therapeutics

The requisite structural features has to be considered to design new molecules with improved activity

In silico modeling for quick prediction of inhibitory activity against 3CLpro enzyme in SARS CoV diseases

Crystal structure of SARS-CoV-2 main protease provides a basis for design of improved α-ketoamide inhibitors

The SARS-coronavirus papain-like protease: structure, function and inhibition by designed antiviral compounds

Disulfiram can inhibit MERS and SARS coronavirus papain-like proteases via different modes

Tanshinones as selective and slow-binding inhibitors for SARS-CoV cysteine proteases

Thiopurine analogues inhibit papainlike protease of severe acute respiratory syndrome coronavirus

Thiopurine analogs and mycophenolic acid synergistically inhibit the papain-like protease of Middle East respiratory syndrome coronavirus

Severe acute respiratory syndrome coronavirus papain-like protease: structure of a viral deubiquitinating enzyme

Identification of severe acute respiratory syndrome coronavirus replicase products and characterization of papain-like protease activity

Thiopurine analogue inhibitors of severe acute respiratory syndrome-coronavirus papain-like protease, a deubiquitinating and deISGylating enzyme

Deubiquitinating and interferon antagonism activities of coronavirus papain-like proteases

Proteolytic processing, deubiquitinase and interferon antagonist activities of Middle East respiratory syndrome coronavirus papain-like protease

The papain-like protease determines a virulence trait that varies among members of the SARS-coronavirus species

Diarylheptanoids from Alnus japonica inhibit papain-like protease of severe acute respiratory syndrome coronavirus

Geranylated flavonoids displaying SARS-CoV papain-like protease inhibition from the fruits of Paulownia tomentosa

Inhibitor recognition specificity of MERS-CoV papain-like protease may differ from that of SARS-CoV

A noncovalent class of papain-like protease/deubiquitinase inhibitors blocks SARS virus replication

X-ray structural and biological evaluation of a series of potent and highly selective inhibitors of human coronavirus papain-like proteases

Severe acute respiratory syndrome coronavirus papain-like novel protease inhibitors: design, synthesis, protein−ligand X-ray structure and biological evaluation

Structure-based design, synthesis, and biological evaluation of a series of novel and reversible inhibitors for the severe acute respiratory syndrome−coronavirus papain-like protease

Textbook of organic medicinal and pharmaceutical chemistry. Lippincott, Williams and Wilkins, Philadelphia 28

Characterization and noncovalent inhibition of the deubiquitinase and deISGylase activity of SARS-CoV-2 papain-like protease

First structureactivity relationship analysis of SARS-CoV-2 virus main protease (Mpro) inhibitors: an endeavour on COVID-19 drug discovery

Chemical-informatics approach to COVID-19 drug discovery: monte carlo based QSAR, virtual screening and molecular docking study of some in-house molecules as papain-like protease (PLpro) inhibitors

Protease targeted COVID-19 drug discovery and its challenges: insight into viral main protease (Mpro) and papain-like protease (PLpro) inhibitors

Chemical-informatics approach to COVID-19 drug discovery: exploration of important fragments and data mining based prediction of some hits from natural origins as main protease (Mpro) inhibitors

Structural and physico-chemical interpretation (SPCI) of QSAR models and its comparison with matched molecular pair analysis

Interpretation of quantitative structureactivity relationship models: past, present, and future

In silico modelling, identification of crucial molecular fingerprints, and prediction of new possible substrates of human organic cationic transporters 1 and 2

CORAL: Building up the model for bioconcentration factor and defining it's applicability domain

A comparative QSAR on 1, 2, 5-thiadiazolidin-3-one 1, 1-dioxide compounds as selective inhibitors of human serine proteinases

SMILES-based optimal descriptors: QSAR modeling of carcinogenicity by balance of correlations with ideal slopes

QSPR modeling for enthalpies of formation of organometallic compounds by means of SMILES-based optimal descriptors

A quasi-QSPR modelling for the photocatalytic decolourization rate constants and cellular viability (CV%) of nanoparticles by CORAL

Large-scale QSAR study of aromatase inhibitors using SMILES-based descriptors

AutoDockVina: improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading

GROMACS: high performance molecular simulations through multi-level parallelism from laptops to supercomputers

Biomolecular simulation: the GRO-MOS96 manual and user guide

The missing term in effective pair potentials

Open Source Drug Discovery Consortium (2014) g_mmpbsa-a GROMACS tool for high-throughput MM-PBSA calculations

The authors have no conflict of interests.