key: cord-0685624-s95i50eh
authors: Amin, Sk. Abdul; Banerjee, Suvankar; Singh, Samayaditya; Qureshi, Insaf Ahmed; Gayen, Shovanlal; Jha, Tarun
title: First structure–activity relationship analysis of SARS-CoV-2 virus main protease (Mpro) inhibitors: an endeavor on COVID-19 drug discovery
date: 2021-01-05
journal: Mol Divers
DOI: 10.1007/s11030-020-10166-3
sha: ebcd682cf505d01f67653ea98c3cd19074d4cca1
doc_id: 685624
cord_uid: s95i50eh

ABSTRACT: Main protease (Mpro) of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) intervenes in the replication and transcription processes of the virus. Hence, it is a lucrative target for anti-viral drug development. In this study, molecular modeling analyses were performed on the structure activity data of recently reported diverse SARS-CoV-2 Mpro inhibitors to understand the structural requirements for higher inhibitory activity. The classification-based quantitative structure–activity relationship (QSAR) models were generated between SARS-CoV-2 Mpro inhibitory activities and different descriptors. Identification of structural fingerprints to increase or decrease in the inhibitory activity was mapped for possible inclusion/exclusion of these fingerprints in the lead optimization process. Challenges in ADME properties of protease inhibitors were also discussed to overcome the problems of oral bioavailability. Further, depending on the modeling results, we have proposed novel as well as potent SARS-CoV-2 Mpro inhibitors. GRAPHIC ABSTRACT: [Image: see text] ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (10.1007/s11030-020-10166-3) contains supplementary material, which is available to authorized users.

At the end of 2019, a new coronavirus-related infection namely severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) had spread its wings across the globe [1, 2] . For its worldwide impact, this coronavirus Extended author information available on the last page of the article disease-2019 (COVID-19) was declared as a global pandemic by the World Health Organization (WHO) and to date over million confirmed cases along with million COVID-19-related mortalities have been reported [1, 3] .

Belonging to the betacoronavirus genus, SARS-CoV-2 is responsible for lower respiratory tract infections similar to severe acute respiratory syndrome coronavirus (SARS-CoV) and Middle-East respiratory syndrome coronavirus (MERS-CoV) [1] . Ongoing research highlighted some important druggable targets like spike (S) protein, papainlike protease (PLpro), RNA-dependent RNA polymerase (RdRp) and SARS-CoV-2 main protease/3C-like protease (Mpro/3CLpro). These possess potentiality to become important targets for achieving the most desirable goal that humanity craves in the current situation [1, 2, 4] . The open reading frame 1ab (ORF 1a/b) of coronaviruses translates polyprotein 1a and polyprotein 1ab. The Mpro and PLpro enzymes produce non-structural proteins by processing these polyproteins which in term aids the production of viral structural proteins [5, 6] . Thus, SARS-CoV-2 Mpro enzyme can be a valuable target as it intervenes in the replication and transcription processes of the virus [2] . It possesses high structural similarity (96% sequential resemblance) to SARS-CoV Mpro [5] .

Additionally, targeting proteases were successful to provide anti-viral agents for the treatment of viral infections like human immunodeficiency virus (HIV) and hepatitis C virus (HCV) [7, 8] . Thus, small molecule-mediated blocking of Mpro activity is a feasible option for SARS-CoV-2 anti-viral drug development [9] [10] [11] [12] [13] [14] [15] [16] [17] [18] . The computer-aided drug design (CADD) and virtual screenings (VS) are viable options. These techniques may be useful to identify promising hit that can aid the design and development of potent anti-viral agents [4] . Meanwhile, drug repurposing was employed as an instant weapon against coronavirus [19] . However, the ongoing rampage of COVID-19 has employed researches in an assignment to discover a permanent solution for this pandemic. In this panorama, the small molecule inhibitors carefully designed by different modeling approaches are one of the most promising tools to achieve success.

Here, we have explored SARS-CoV-2 Mpro inhibitors by different molecular modeling strategies with four main mottos-(i) development of a mathematical relationship between the derivatives and SARS-CoV-2 Mpro enzyme (ii) identification of important fingerprints that module the SARS-CoV-2 Mpro inhibition, (iii) scope of these derivatives to address ADME properties, (iv) design of potent SARS-CoV-2 Mpro inhibitors with significant ADME properties. The current study, a part of our rational drug design and discovery program, [4, [19] [20] [21] may offer an initiative to explore the possibility of potent inhibitor design against the Mpro enzyme of SARS-CoV-2.

Dataset A number of 33 derivatives, represented by SARS-CoV-2 Mpro inhibitory activity IC 50 (µM), were obtained from the published data [5, 6, 9, 14, 15] . The SARS-CoV-2 Mpro inhibitory activity values of the inhibitors are presented in Supplementary Table S1 . The pIC 50 (i.e., -log IC 50 ) values were used to derive QSAR models [22] [23] [24] .

The classification modeling assists to classify the active and inactive molecules in terms of their biological data [25] [26] [27] [28] [29] [30] . Here, we employed Bayesian classification approach [31] [32] [33] .

Performing Bayesian classification study by the aid of Discovery Studio (DS) software [34] enables graphical visualization of critical chemical sub-structural features (fingerprint or fragments) attributed to enhance or decrease the SARS-CoV-2 Mpro inhibitory activity. Additionally, as to conduct this classification-based study, on the basis of their SARS-CoV-2 Mpro inhibitory activity, the dataset molecules were grouped into active (SARS-CoV-2 Mpro pIC 50 > 5.0) and inactive (SARS-CoV-2 Mpro pIC 50 < 5.0) molecules (e.g., active = 1, inactive = 0) [23] .

The selection of the training and test sets was done by using Generate training and test data tool in DS [34] . The whole data were divided into 20 clusters by maximum dissimilarity approach on the basis of Predefined Set properties including ALogP, Molecular_Weight, Num_H_ Donors, Num_H_Acceptors, Num_RotatableBonds, Num_ Atoms, Num_Rings, Num_AromaticRings, Num_Fragments, Molecular_PolarSurfaceArea. The whole data set compounds were separated into two groups, a training set SARS-CoV-2 Mpro inhibitors, a test set SARS-CoV-2 Mpro inhibitors (Supplementary Table S1 ).

Further, to ensure whether the selected test set compounds truly represent the training set or not, principal components analysis (PCA) was performed by Calculate principal components tool in DS [34] . The DS default properties such as ALogP, Molecular_Weight, Num_H_ Donors, Num_H_Acceptors, Num_RotatableBonds, Num_ Rings, Num_AromaticRings, Molecular_FractionalPolar-SurfaceArea were considered for the PCA calculation. The uniform distribution of the test set SARS-CoV-2 Mpro inhibitors in the PCA three-dimensional plot (as given in Supplementary Figure S1 ) referred a proper division of the training and the test sets.

Finally, the Bayesian classification model was constructed on the training set and was cross-validated by using the test set. Before conducting this Bayesian classification study, several fundamental molecular features namely, ALogP, Molecular_Weight, Num_H_Donors, Num_H_Acceptors, Num_RotatableBonds, Num_Rings, Num_AromaticRings, Molecular_FractionalPolarSurfaceArea of the dataset molecules have been calculated [34] . Alongside those molecular properties, a topological fingerprint descriptor namely extended connectivity fingerprint of diameter 6 (ECFP_6) [35] was also considered for this study. The quality of this classification model was evaluated using the Receiver operating characteristics (ROC) plot [36] , sensitivity (Se), specificity (Sp) and accuracy (Acc) for both the training and the test sets [23] .

The derivatives with no activity and without definite SARS-CoV-2 Mpro inhibitory activity were not considered for the multiple linear regression (MLR) analysis [23] . Hence, only 25 molecules were recognized for the regression-based QSAR study (Supplementary Table S1 ).

Meanwhile, a number of 2D and fingerprint descriptors were calculated [37] . Then, the descriptors with constant values were removed from the data matrix [23] . Next, the highly inter-correlated variables were stocked out depending on the specified variance of 0.001 and correlation coefficient cutoff values of 0.99 [38, 39] . Then, several genetic function approximation (GFA) runs were employed to collect a bunch of important descriptors [39] . Finally, stepwise multiple linear regression (S-MLR) model was developed to identify the linear correlation between the structure of SARS-CoV-2 Mpro inhibitors and their respective Mpro inhibitory activities. The robustness of the constructed model was justified by correlation coefficient (R), adjusted R 2 (R A 2 ), variance ratio (F) at specified degrees of freedom (df), cross-validated R 2 (Q 2 ), standard error of estimate (SEE), and other validation metrics [23] . In addition, Euclidean distance-based applicability domain was also constructed [23, 38] to check the applicability of the MLR model.

For the docking studies, the SARS-CoV-2 Mpro structure was obtained from Protein Data Bank (PDB ID: 6LZE). Subsequently, the compounds were docked in the active site of the Mpro protein using Auto Dock Vina v1.1.2 [40] , wherein a grid box of size 16, 14, and 14 with spacing of 1 Å were set around the active site of SARS-CoV-2 Mpro.

Later, the molecular dynamics simulation was performed by the GROMACS 5.1.4 version [41] using the GRO-MOS43A2 force field and SPC/E water model. To neutralize the charges on each simulating system, an appropriate number of ions (Na + ) were added. The energy of each system was minimized by the steepest descent algorithm followed by NVT (at 300 K) and NPT (at 1 bar) ensemble equilibrations for 100 ps. Subsequently, each of the equilibrated system was carried on for the production simulation of 20 ns. The trajectory data of the production simulations were further used for the calculation of root mean square deviation (RMSD), root mean square fluctuations (RMSF), and radius of gyration (Rg) data of each system. The binding energy of each compound in the complex was calculated by g_MMPBSA package of GROMACS, for every 0.1 ns frame of each 20 ns simulation [42] .

SARS-CoV-2 Mpro is a homodimeric protein. Each subunit is termed as protomer. A number of 306 amino acid residue is found in each protomer. It is constructed by three domains [4, 9] . The domain I is 8 to 100 residues long followed by domain II (101 to 184 residues) and domain III (from 199 to 306 residues). Besides, domains II and III are bridged by a long loop (from 185 to 198 residues) [9] .

Domains I and II allocated the same fold i.e., an antiparallel six stranded β-barrel structure, while domain III is semblance by five α-helices arranged into a largely antiparallel globular cluster. Meanwhile, the domain III helps in the regulation of Mpro dimerization through an inter subunit salt-bridge between E290 from one protomer and R4 from the other protomer. Notably, the substrate-binding site or catalytic site of SARS-CoV-2 Mpro is located at a cleft between domains I and II. The N-terminal amino acid residue of a protomer namely S1 interacts with the E166 of another to form the S1 subsite of the substrate-binding pocket. Hence, the dimerization is essential for protease activity [17] .

The research on SARS-CoV-2 Mpro has moved at a much faster after delivering several ligand bound crystal structures. Those have provided useful information for developing inhibitors, but it seems that it is not enough. Analysis of different crystal structures has shown that there is an intrinsic flexibility in the catalytic site. In order to explore the detail binding interactions, few contour maps of the binding site of the SARS-CoV-2 Mpro (PDB: 6WTT) were determined by Display receptor surfaces tool of DS [34] . Six structure-based contour maps for hydrophobic, hydrogen bond, charge, aromatic, ionizability and solvent accessible surface (SAS) are provided in Fig. 1 . (Fig. 1a) . Analyzing Fig. 1a , b suggests that hydrophobic aromatic substitution may be allowable in the binding pocket. The hydrogen-bond site map (Fig. 1c) shows that acceptor feature exists close to the straight chain amide residues. The S1 cavity is acceptor specific. Near S1 site, H172 endorses acidic ionizability (Fig. 1d) , and it is slight negatively charged (Fig. 1e) . Significantly, ionizability and interpolated charge contours are more or less consistent. From the SAS contour (Fig. 1f) , it may be suggested that a significant part of the catalytic site is solvent exposed. In order to explore the details contribution of fragments/ fingerprints of the inhibitors, we moved forward to quantitative structure-activity relationship (QSAR) studies and design of specific SARS-CoV-2 Mpro inhibitors.

The Bayesian classification modeling is a classification QSAR technique based on the Bayes' theorem which utilizes data to predict the probability of specific events [43] [44] [45] . Additionally, another advantage of this Bayesian classification study with fingerprint descriptor is its capability to recognize important structural fragments of molecules while indicating their positive or negative influence on the activity [23] .

In order to describe the statistical quality of the generated Bayesian classification model, different statistical parameters like sensitivity (Se), specificity (Sp), and accuracy (Acc) were calculated. The results were found to be statistically significant as all the parameters were having decent scores to consider the model as robust and predictive as specified in Table 1 .

Further, the ROC (Receiver operating characteristic) curve for training and test sets are found to be 0.747 and Figure S2 . The mechanistic interpretation of the Bayesian classification study is performed using a fingerprint descriptor ECFP_6. A set of 20 good and 20 bad molecular sub-structural features has been procured with positive and negative influences, respectively, on SARS-CoV-2 Mpro inhibition of the compounds. Twenty good (G1-G20) and twenty bad sub-structural fragments (B1-B20) as constructed from the ECFP_6 fingerprint descriptor are shown in Supplementary Figure S3 and S4, respectively.

Upon observation, the set of 20 good molecular substructures can be clustered into four major groups namely: bi-acetyl amine and acetamido group containing 2-oxo pyrrolidine moiety (G1, G7 and G9-G15), cyclohexyl and cyclohexyl methyl groups (G2-G3, G8, G16-G17 and G19-G20), and acetamido methylene (iso-butyl) acetamide moiety (G4-G6). Beside these frequent sub-structures, the oxyanion function (G18) is upheld as positive influencers for the Mpro inhibition as shown in Supplementary Figure S3 .

In contrast, among the proposed negatively influencing features, 4-fluoro phenyl and 4-fluoro benzyl moieties are the most commonly displayed bad features (B11-B15 and B17). The branched alkyl (B3-B4, B6, and B8) and amino-alkyl (B1-B2, B5, B10, and B16) groups are indicated to be detrimental for activity. Moreover, oxymethylene carbonyl (B7) and acetate (B9) functions are also suggested as negative regulators of Mpro inhibitory activity ( Supplementary Figure S4 ).

Further analysis of the fragments and the dataset molecules, it is found that the most active M027 having an acetamido methylene (iso-butyl) acetamide function and a negatively charged oxygen ion similar to the sub-structures G4-G6 and G18, respectively (Fig. 2) .

From the crystal structure analysis of M027 with SARS-CoV-2 Mpro (PDB: 6WTT), the mentioned fragments are found to involve in several interactions at the enzyme active site. The iso-butyl group of M027 enters into the S2 pocket of the enzyme while the carboxamide function interacts with Q189 and E166 amino acids (Fig. 3) [15] .

Meanwhile, fragments G2-G3, G8, G16-G17 and G19-G20 exhibit the importance of cyclohexyl moiety. From the SARs, it may be observed that cyclohexyl function is important for the activity. The cyclohexyl methyl moiety is found in active molecules like M009, M011, and M012. The cyclohexyl function embeds itself in the hydrophobic S2 site of SARS-CoV-2 Mpro [6] . Therefore, hydrophobic interactions are essential in these regions. Similarly, the (S)γ-lactam ring is directed to a hydrophobic S1 pocket (Fig. 3) . The cyclohexyl methyl group of M009 is found to be important for entering into the S2 pocket of the enzyme where its indole ring enters into the S4 pocket (PDB: 6M0K). Meanwhile, the indole moiety of compound M010 also shows identical binding to that of M009 (PDB: 6LZE) [6] .

A substituted 2-pyridinone moiety is present in both compound M012 and M013 whereas the 3-amino-Boc substituted 2-pyridinone moiety of M013 forms more than one interaction at the active site of SARS-CoV-2 Mpro (PDB: 6Y2F) as shown in Fig. 3 . Also, the 2-carbonyl and the 3-amino functions of the moiety interact with E166 through hydrogen bond formation [5] . The presence of 2-phenyl-4-chromenone moiety can be observed in both M033 and M032. The 6-and 7-hydroxyl group of M033 interacts with L141 and G143 (PDB: 6M2N), respectively [14] .

Regarding the bad molecular fragments, compound M021 possessing a 4-fluorophenyl group is inactive (Fig. 4) . The acetate function containing M029 exhibits inactivity against Mpro. The oxymethylene carbonyl moiety containing M030 also shows inactivity against Mpro (Fig. 4) .

An effective drug candidate/drug-like ligands having promising biological responses should possess the ability to reach its desire domain in sufficient concentration. Drug design and discovery obviously depends on the assessment of absorption, distribution, metabolism and excretion (ADME) characteristics.

In order to check the drug-likeliness of the investigated derivatives Filters ligands using Lipinski and Veber rules protocol of DS was employed [34] . It selects drug-like [46] and Veber [47] . The default settings for Lipinski Rule of Five (Hydrogen Bond Donors: 5, Hydrogen Bond Acceptors: 10, Molecular Weight: 500, AlogP: 5, Number of Violations Allowed: 1) and Veber Rule (Rotatable Bonds 10, Polar Surface Area 140, Hydrogen Bond Donors, and Acceptors 12) were considered for this study.

Notably, 19 compounds (Fig. 5) fail to pass the Lipinski and Veber rules. Only 14 compounds (Fig. 5) pass these two rules, therefore, those have a higher probability of good oral bioavailability.

The protease targeted peptidomimetic inhibitors design is very challenging due to their undesirable pharmacokinetic properties. In contrast, compounds with low molecular weight or non-peptidomimetics exhibit good druglikeliness. However, non-peptidomimetics/low molecular weight derivatives fail to effectively block the proteolytic activity of SARS-CoV-2 Mpro. In these circumstances, the structure of SARS-CoV-2 Mpro in complex with a small molecule baicalein (M033) may be a good option for baicaleinderived lead optimization. The binding mode of baicalein at the active site of SARS-CoV-2 Mpro facilitated a unique protein-ligand interaction pattern.

Since baicalein (M033) possesses a molecular weight of 270.24 Da, it encourages us to anticipate new derivative design by keeping the baicalein core. Lead optimization of baicalein with good molecular fingerprint (as suggested by the Bayesian modeling study) may render new derivatives directed toward S1 and/or S4 site(s). It may effectively block the proteolytic activity of SARS-CoV-2 Mpro.

Taken together these modeling efforts may give rise in a new candidate with broad-spectrum anti-viral properties.

However, substitution at the wrong position (as in case of baicalin) resulted in ~ sevenfold loss in SARS-CoV-2 Mpro inhibition (baicalein vs baicalin) [14] .

Considering the finding of the performed QSAR studies, we have designed a set of four chromenone-based molecules (Fig. 6 ).

Primarily, the Bayesian classification model was used to predict the Mpro inhibitory activity of these molecules. The designed compounds (D1-D4) predicted as active. Hence, these compounds may serve as promising molecules against SARS-CoV-2 Mpro.

To further revalidate prediction credibility, a stepwise multiple linear regression (S-MLR) model has been constructed on the available data. At first, a pool of 2D and fingerprint descriptors for these derivatives was calculated [37] . Then, dataset thinning was introduced followed by several genetic function approximation (GFA) runs [38, 39] . The best model (Eq. 1) through the S-MLR analysis (the stepping criterion of F = 4 for inclusion and F = 3.99 for exclusion) is as follows Table 2 . Other details are given in the Supplementary files (Table S2-S4) .

Additionally, Euclidean distance-based applicability domain was constructed [23, 38] as illustrated in Fig. 7 . It justifies that all the compounds are within the boundary of the hypothetical domain (Fig. 7) . Hence, there is no outlier for this dataset. The designed molecules (D1-D4) predicted pIC 50 more than 7.523 as depicted in Table 3 . This result supports the potential of these designed molecules to become promising Mpro inhibitors.

Since the drug-likeliness is one of the major challenges in drug discovery. The drug-likeliness of these designed molecules (D1-D4) was investigated using the DruLiTo software [48] . The drug-like properties of the designed molecules are also given in Table 3 . Fingerprint O=C-C-C-C-C-N Negative

To understand the structural basis of inhibition by compounds (D1-D4), protein-ligand docking studies were performed using AutoDockVina [40] . The outcome of molecular docking shows that all compounds are found to be docked into the active site of SARS-CoV-2 Mpro (Fig. 8) .

The binding energies of the selected conformation of the compounds are depicted in the Supplementary Table S5 , which indicates that the complex with compound D4 has higher binding energy in comparison with the other compounds. Moreover, interacting residues of the docked complexes also reveals similarity with the interacting residues of the reported protein co-crystal structure (PDB: 6LZE).

The respective average RMSD of apo, prt_D1, prt_D2, prt_D3, and prt_D4 are enumerated as 0.284, 0.292, 0.213, 0.258, and 0.276 nm, respectively. Wherein, protein shows lesser deviations in the structure with compound D2 in comparison with the other compounds as well as apo form (Fig. 8) . The analysis infers an increase in the stabilization of protein backbone structure after interaction with D2 during the simulation. Simultaneously, the fluctuations in the backbone atom of the protein residues in each system was analyzed by RMSF, which shows apo, prt_D1, prt_D2, prt_D3, and prt_D4 system to possess average RMSF of 0.135, 0.137, 0.117, 0.138, and 0.126 nm, respectively. However, protein residues presented lower fluctuation after interacting with compound D2 than the other complex and Further, the comparative analysis of Rg data was performed to determine the protein compactness after interaction with compounds. The respective average Rg value of 2.13, 2.10, 2.13, 2.10, and 2.12 nm are computed for apo, prt_D1, prt_D2, prt_D3, and prt_D4, respectively. The obtained data reveals that all protein complex forms attain level of compactness similar to the apo form, which indicates that each compound interact with protein without disturbing its structural folding in the dynamic environment (Fig. 8) . Altogether, our study highlighted that the compound D2 has shown much stable interaction along with the induction of low deviations and fluctuations in the protein structure as compared to apo form and other compounds during MD simulation.

The affinity of the compounds with protein was also analyzed by the binding energy calculation. The average binding energy calculated for each protein-ligand complex is presented in Table 4 , which shows that the compound D4 has more binding affinity with Mpro protein during the simulation.

Wherein, the van der Waals energy plays a major role in the binding of compound D4 at the active site of the Mpro in comparison to the other free energies. Moreover, the binding energy analysis shows corroboration with docking studies of the compounds. It infers that the compound D4 has more affinity for the static and dynamic SARS-CoV-2 Mpro structure.

COVID-19 shows worldwide impact as a global pandemic. Till date over million confirmed cases have been reported worldwide. In this communication, QSAR analyses were performed on recently reported structurally diverse Mpro inhibitors to understand structural requirements for higher activity. The study is able to extract the significant molecular attributes of these SARS-CoV-2 Mpro inhibitors.

The main problems for design of SARS-CoV-2 Mpro inhibitors are the perfect binding of susbtituents in putative binding site and the ADME properties. To overcome these problems, we suggest baicalein-derived design as well as lead optimization. Since different ligands induce different conformational changes, the conformation of binding pocket residues could not be easily predicted for different inhibitors. Nonetheless, the S1, S1', S2, and S4 pockets bear intrinsic flexibility where hydrophobic susbtitutents may trigger the SARS-CoV-2 Mpro inhibition. Our structure-based contours result suggests that Mpro binding pockets should be analyzed carefully to design inhibitors with such flexibilities.

In our previous study, the Monte Carlo optimizationbased QSAR, structural and physico-chemical interpretation (SPCI) analysis were successful to deliver several important molecular features from the SARS-CoV Mpro inhibitors [21] . This can be useful to develop effective inhibitors against SARS-CoV-2 Mpro. Additionally, compared to the recent attempts to identify the promising attributes for previous coronavirus inhibitors (Table 5) , the current study deals with the existing SARS-CoV-2 Mpro inhibitors. In summary, the modeling results provide useful quantitative and qualitative information about the structural requirements of an effective Mpro inhibitor against SARS-CoV-2.

The Trinity of COVID-19: immunity, inflammation and intervention

Structure-based drug designing and immunoinformatics approach for SARS-CoV-2

WHO Coronavirus Disease

Fight against novel coronavirus: a perspective of medicinal chemists

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

Structure-based design of antiviral drug candidates targeting the SARS-CoV-2 main protease

Targeting proteases: successes, failures and future prospects

Emerging principles in proteasebased drug discovery

Structure of Mpro from SARS-CoV-2 and discovery of its inhibitors

2020) α-Ketoamides as broad-spectrum inhibitors of coronavirus and enterovirus replication: structure-based design, synthesis, and activity assessment

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

COVID-19: drug targets and potential treatments

Activity profiling and structures of inhibitor-bound SARS-CoV-2-PLpro protease provides a framework for anti-COVID-19 drug design

Discovery of baicalin and baicalein as novel, natural product inhibitors of SARS-CoV-2 3CL protease in vitro

Boceprevir, GC-376, and calpain inhibitors II, XII inhibit SARS-CoV-2 viral replication by targeting the viral main protease

Targeting the dimerization of the main protease of coronaviruses: a potential broad-spectrum therapeutic strategy

Prediction of novel inhibitors of the main protease (M-pro) of SARS-CoV-2 through consensus docking and drug reposition

AI-aided design of novel targeted covalent inhibitors against SARS-CoV-2

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

Gupta SP (ed) Viral proteases and their 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

Quantitative prediction of imprinting factor of molecularly imprinted polymers by artificial neural network

First report on the structural exploration and prediction of new BPTES analogs as glutaminase inhibitors

In silico studies of novel scaffold of thiazolidin-4-one derivatives as anti-Toxoplasma gondii agents by 2D/3D-QSAR, molecular docking, and molecular dynamics simulations

QSAR modeling of endpoints for peptides which is based on representation of the molecular structure by a sequence of amino acids

Classical and 3D QSAR studies of larvicidal monoterpenes against Aedes aegypti: new molecular insights for the rational design of more active compounds

A comprehensive computational analysis of cathinone and its metabolites using quantum mechanical approaches and docking studies

First report on predictive chemometric modeling, 3D-toxicophore mapping and in silico screening of in vitro basal cytotoxicity of diverse organic chemicals

The application of machine learning techniques to innovative antibacterial discovery and development

Development of classification and regression based QSAR models to predict rodent carcinogenic potency using oral slope factor

Novel naïve Bayes classification models for predicting the chemical ames mutagenicity

Classification of kinase inhibitors using a bayesian model

Novel Bayesian classification models for predicting compounds blocking hERG potassium channels

Discovery studio. Biovia

Extended-connectivity fingerprints

An introduction to ROC analysis

PaDEL-descriptor: an open source software to calculate molecular descriptors and fingerprints

The simple, user-friendly and reliable online standalone tools

Nano-BRIDGES" software: open access tools to perform QSAR and nano-QSAR modelling

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

Open Source Drug Discovery Consortium

g_mmpbsa-a GROMACS tool for high-throughput MM-PBSA calculations

Chemometrics tools in QSAR/QSPR studies: a historical perspective

Ecotoxicological assessment of pharmaceuticals and personal care products using predictive toxicology approaches

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

Experimental and computational approaches to estimate solubility and permeability in drug discovery and development settings

Molecular properties that influence the oral bioavailability of drug candidates

Drug Likeness Tool (DruLiTo), software

SARS-CoV Mpro inhibitory activity of aromatic disulfide compounds: qSAR model

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

In silico modeling of small molecule carboxamides as inhibitors of SARS-CoV 3CL protease: an approach towards combating COVID-19

Development of a simple, interpretable and easily transferable QSAR model for quick screening antiviral databases in search of novel 3C-like protease (3CLpro) enzyme inhibitors against SARS-CoV diseases

Publisher's Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations

The authors have no conflict of interests.

Sk. Abdul Amin 1 · Suvankar Banerjee 1 · Samayaditya Singh 2 · Insaf Ahmed Qureshi 2 · Shovanlal Gayen 3