key: cord-0805572-kakxmc14
authors: Achutha, A. S.; Pushpa, V. L.; Suchitra, Surendran
title: Theoretical Insights into the Anti-SARS-CoV-2 Activity of Chloroquine and Its Analogs and In Silico Screening of Main Protease Inhibitors
date: 2020-09-22
journal: J Proteome Res
DOI: 10.1021/acs.jproteome.0c00683
sha: 4ef32bacf98fffe587c960dba4451dfc533fdf83
doc_id: 805572
cord_uid: kakxmc14

[Image: see text] Corona virus disease (COVID-19) is a dangerous disease rapidly spreading all over the world today. Currently there are no treatment options for it. Drug repurposing studies explored the potency of antimalarial drugs, chloroquine and hydroxychloroquine, against SARS-CoV-2 virus. These drugs can inhibit the viral protease, called chymotrypsin-like cysteine protease, also known as Main protease (3CL(pro)); hence, we studied the binding efficiencies of 4-aminoquinoline and 8-aminoquinoline analogs of chloroquine. Six compounds furnished better binding energies than chloroquine and hydroxychloroquine. The interactions with the active site residues especially with Cys145 and His41, which are involved in catalytic diad for proteolysis, make these compounds potent main protease inhibitors. A regression model correlating binding energy and the molecular descriptors for chloroquine analogs was generated with R(2) = 0.9039 and Q(2) = 0.8848. This model was used to screen new analogs of primaquine and molecules from the Asinex compound library. The docking and regression analysis showed these analogs to be more potent inhibitors of 3CL(pro) than hydroxychloroquine and primaquine. The molecular dynamic simulations of the hits were carried out to determine the binding stabilities. Finally, we propose four compounds that show drug likeness toward SARS-CoV-2 that can be further validated through in vitro and in vivo studies.

In December 2019, patients were admitted in hospitals in Wuhan, China due to pneumonia with fever and other symptoms. 1, 2 In January 2020, it was discovered that a novel coronavirus, which is similar to SARS-CoV and MERS-CoV genomes, is responsible for this disease. 2 This virus is named as Severe Acute Respiratory Syndrome Corona Virus-2 (SARS-CoV-2) and the disease is called Corona Virus Disease-2019 . 3, 4 It was declared as pandemic on March 11 by the World Health Organization (WHO). 5 As of September 1, 2020, there are 25,327,098 confirmed cases with 848,255 deaths worldwide according to WHO reports. 6 Coronavirus is a single-stranded positive sense RNA with a membrane envelope composed of four types of viruses: α, β, γ, and δ coronaviruses. SARS-CoV-2 belongs to β-coronavirus family. 5, 7 Its genome encodes for structural proteins mainly spike glycoprotein (S), envelope glycoprotein (E), membrane glycoprotein (M), nucleocapsid protein (N), and nonstructural proteins consisting of chymotrypsin-like cysteine protease also known as Main protease (3CL pro ) and papain-like protease (PL pro ). 1, 8, 9 The viral RNA is released to the host cell, which is then translated to polyproteins inside the cell. These viral polyproteins are cleaved to functional proteins by proteolysis, which is catalyzed by 3CL pro and PL pro 1,8,9 . Since the main protease or 3CL pro cleaves at 11 sites of the polyprotein to produce smaller proteins for viral replication and because of its nonsimilarity with human proteins, 3CL pro is a potential target in anti-COVID-19 drug design. 8, 10 The main protease of SARS-CoV-2 has a Cys-His catalytic diad involving Cys145 and His41 and has four conserved binding subsites S1, S2, S4, and S1′. Main chains of Cys145, Phe140, and Leu141 and side chains of Asn142, Glu166, His163, and His172 residues are present in the S1 subsite. The side chains of His41, Met49, Tyr54, and Met165, along with the alkyl side chain of Asp187 forms the S2 subsite. The S4 subsite consists of Met165, Leu167, Phe185, and Gln192 side chain residues and the main chain of Gln189. S1′ includes Thr24 and Thr25. Therefore, a molecule to show 3CL pro inhibitory property should interact with one or more residues among these amino acids. Also, the inhibitory action can be assumed to enhance if they interact with Cys145 or His41, which is involved in the catalytic diad of Main protease 10 ( Figure 1 ). Cysteine acts as a nucleophile and Histidine as an acid/base catalyst in this catalytic reaction. 9, 11 Currently there are no specific approved treatment options for COVID-19. Researches are in search of vaccine development, druggable small molecule, monoclonal antibodies and cell-based therapies. 12 Several drug repurposing studies have been conducted so far, and many of them gave positive results. The 4-aminoquinolines, chloroquine (CQ), and its hydroxyl analog hydroxychloroquine (HCQ) have been found to disrupt the viral replication and infectivity. 13 Their broad spectrum of activity as antibacterial, antiviral, and antifungal infections has also been reported. 15, 16 HCQ has been found out to be less toxic than CQ, and it is currently being used to treat COVID-19 patients. 17 The present study deals with theoretical perspectives on the anti-SARS-CoV-2 activity of CQ and HCQ analogs by observing their main protease inhibitory property. Molecular docking studies with the 3CL pro protein were performed to analyze the drug likeness as well as to correlate the binding energy of the docked complex with various physicochemical properties of the active molecules, which will aid in the design of new anti-COVID-19 medicatives.

Analogs of CQ and HCQ as well as primaquine derivatives were obtained through literature search, from the database of small molecules Pubchem as well as created manually. 18−22 All the molecules were drawn in Marvin Sketch and minimized their energy using Avogadro software. 23, 24 The structure of 3CL pro protein having PDB ID 6LU7 was retrieved from RCSB Protein Data Bank. 10, 25 Characterization studies of the protein were done by using Expasy-protparam, an online server, which provides all the details regarding the protein.

The energy minimization of the retrieved protein was done by employing Swiss PBD Viewer. 26, 27 Molecular docking was done using Autodock 4.2.6. 28 Both protein and ligands were prepared in pdbqt format. Polar hydrogens and Gasteiger charges were added to the receptor. Grid was prepared with 60 × 60 × 60 Å 3 with spacing 0.375 Å and centered at −18.617, 13.208, and 60.061 Å along x, y, and z axes, respectively. Grids were prepared for every atom present in the ligand data set. The genetic algorithm was employed as search parameter with 50 runs, 300 and 27,000 population size and number of generations, respectively.

A multiple linear regression (MLR) model was generated to find the correlation with binding energy obtained from docking and the physicochemical properties using QSARINS. 29, 30 The molecular descriptors for the ligands were calculated utilizing PaDEL-Descriptor. 31 To build the model, 80% of the data set were randomly divided as training set and test set in 80:20 ratio. Thirty-five compounds (28 training and 7 test) including CQ and HCQ were employed for modeling purpose (Supplementary Tables S2 and S3 ). The models were generated using training set with MLR analysis of 1−5 variables. The generated models were examined by internal validation and external validation. Internal validation was done by cross-validated leave-one-out (LOO) method, which involves iteratively leaving one compound from the training set and generating regression model with the remaining molecules and predicting the value of response for the excluded one. For a good model, the regression coefficient (R 2 ) as well as the cross validated R 2 , that is, Q 2 has to be close and >0. 6. 29,32 The model was validated using test set, which was not included in model generation. 29, 33 The best model was used to predict the binding energy of novel set of compounds.

All the MD simulations were done using the NAMD package (Version 2.14) developed by the Theoretical and Computational Biophysics Group in the Beckman Institute for Advanced Science and Technology at the University of Illinois at Urbana−Champaign 34 with the CHARMM 36 force field. 35, 36 The docked structures were used as the initial coordinates for MD simulations. TIP3P water box was used to solvate the complex with 5 Å padding. The solvated system was minimized for 100 ps with a time step of 2 fs. The simulations were done with Particle mesh Ewald electrostatics and periodic boundary conditions. Atom-based cutoff of 12 Å was applied for nonbonding interactions. The temperature and pressure were kept constant at 310 K and 1 atm, respectively, using Langevin dynamics. The system was equilibrated for 100 ps followed by a production simulation for 20 ns using the NPT ensemble. All the preparation and analysis steps were done using the VMD package 1.9.3 (available at http://www. ks.uiuc.edu/Research/vmd/) (see Figure 2 ). 37 

Being an inevitable enzyme for replication and transcription of the virus, the 3-chymotrypsin-like-cysteine protease (3CL pro ) having PDB ID 6LU7 was retrieved from RCSB PDB and characterized by primary and secondary structure analysis as shown in Table 1 .

The pI value (isoelectric point) of 5.95 shows the slightly acidic nature of the protein. The aliphatic index (AI) is a measure of the relative volume occupied by the side chains Ala, Leu, Ile, and Val, which are the aliphatic residues in a protein.

High AI ranging from 66.5 to 84.33 implies high thermal stability and hydrophobicity, which help it for membrane penetration in biological system. 39, 40 The value for grand average of hydropathy (GRAVY) indicates the hydrophilicity of the protein. 40 The GRAVY value of −0.019, which is close to zero, shows that the protein is hydrophobic. Instability index <40 shows that in vivo half-life is greater than 60 h, which indicates the stability of the protein. 41, 39 All the parameters show that the 6LU7 is a stable, hydrophobic protein.

The 6LU7 protein consists of three domains in which the catalytically active site is present between Domain I and II. Figure S1 (Supporting Information) shows the domains and active site of 6LU7 protein bound with its cocrystallized inhibitor N3. 10,42

The antimalarial drugs chloroquine (CQ) and its hydroxy analog hydroxychloroquine (HCQ) are currently used as drugs for the treatment of COVID-19. 13, 17 To compare its efficiency as an anti-COVID-19 drug to that of an antimalarial drug, docking studies were carried out for both targets, 3CL pro and PfDHFR-TS. Plasmodium falciparum dihydrofolate reductase thymidylate synthase (PfDHFR-TS) is one of the important targets for antimalarial drugs 43, 44 (Supplementary Table S1 ).

It was identified that for both the targets HCQ has lower binding energy owing to its higher binding affinity to the protein. Both the drugs are more active on PfDHFR-TS than 3CL pro . The 2D interaction diagrams and docked poses of CQ and HCQ with the target 3CL pro are given in Figures S2 and S3, respectively. HCQ fits more perfectly to the active site of the target than CQ, and hence, HCQ has high negative binding energy, which imparts more stability to the docked complex (Supplementary Figure S3) .

CQ forms three hydrogen bonds with Gly143, Cys145, and His164; van der Waals interaction with Asn142 and Gln189; πalkyl interaction with His163; and two π-donor hydrogen bond interactions with the −SH group of Cys145 at distances 3.78 Å, 5.16 Å, 5.31 Å, 4.30 Å, 5.10 Å, 6.63 Å, 4.93 Å, and 4.95 Å, respectively (Supplementary Figure S2a , Table 2 ). The Cys145-His41 diad is blocked by these interactions. The interaction of alkyl group of the CQ with the π-cloud of imidazole ring of histidine may disturb the delocalization of π electrons. The basicity of the imidazole ring is due to the resonance stabilization of the positive imidazolium ion ( Figure  3 ). 45 The alkyl−π interaction of the alkyl donor and aromatic acceptor increases the stability of the docked complex as well as reduces the ability of the imidazole ring to accept H + from the Cys residue, which is the initial step in the catalytic cycle in proteolysis, thereby reducing the activity. 46, 11 HCQ forms two hydrogen bonds with Phe140, and one each with Asn142, Ser144, and Glu166 at distances 5.15 Å, 5.0 Å, 3.94 Å, 4.31 Å, and 5.04 Å, respectively. It forms two van der Waals interactions with Phe140, and one each with Glu166 and His172 at distances 6.59 Å, 7.12 Å, 5.30 Å, and 5.63 Å, respectively. It also forms a π−σ interaction with Asn142 and two π−alkyl interactions with Cys145 (Supplementary Figure  S2b , Table 2 ). The π−alkyl interaction with Cys145 may disturb the conformation required for the interaction with polypeptide in the catalytic cycle.

There have been reports about the entry of SARS coronavirus as acidic pH dependent in which the activation occurs by the fusion of the membranes of virus and cellular endosomes and viral genome enters the cytoplasm. 47, 48, 15 The presence of amino side chain in CQ and HCQ makes them basic and increases the pH of cell organelles like endosomes and lysosomes, which interferes with the replication and inhibits the viral activation. 17,15

To find more potent inhibitors of 3CL pro , analogs of chloroquine were subjected to molecular docking, which were obtained from literature, Pubchem database, by replacement with isosteres and functional group variations 18−20 of chloroquine (Supplementary Tables S2 and  S3 ). The binding energy thus obtained from docking with 6LU7 protein is given in Table 2 .

All the chloroquine analogs show some of the key interactions with the active site residues. The negative value shows a release of energy while forming a protein−ligand docked complex, which imparts stability. The more negative the binding energy, the higher will be the stability and binding affinity. The compounds 29, 25, 22, 8, 24 , and 10 have more negative binding energy than CQ and HCQ. Their 2D ligand interaction diagrams are shown in Supplementary Figure S4 . From the structures of these six compounds, we can infer that if the side chain has an extra −NH group, it can introduce additional hydrogen bonds, which results in a decrease in binding energy. In compound 24, the tail of side chain has become a part of a ring, which can interact with Met165 residue via π−alkyl interaction. Because of the presence of an additional −NH group, the basicity of these compounds is higher than CQ and HCQ, and thus, they can act as better inhibitors against SARS-CoV-2.

Compounds 29, 25 22, 8, 24 , and 1 form π-donor hydrogen bonds with Cys145, which will make the −SH hydrogen less available for catalytic cycle. Compounds 29, 25, 22, 24, and 1 have a π−alkyl interaction with Cys145. Overall, these interactions may affect the catalysis of proteolysis reaction Journal of Proteome Research pubs.acs.org/jpr Article and hence block viral replication. Compound 29 is primaquine (an 8-aminoquinoline), which is an antimalarial drug. 49, 50 The effect of the side chain in the chloroquinoline scaffold was studied by changing the chain length and methyl substitution of the side chain. Compounds 12 and 23 are obtained by the removal of the branching methyl group from CQ and HCQ structures, which results in an increase in the binding energy to −5.23 and −5.80 kcal/mol, respectively. Therefore, it is understood that the branching next to the −NH− group in 4-aminopentyl side chain is important. Isosteric substitution of −CH 3 group in compounds 21, 22, and 26 with −OH, −NH 2 , and −CF 3 results in binding energies −5.75, −7.21, and −5.25 kcal/mol, respectively. Substitution with NH 2 resulted in a decrease in binding energy compared to HCQ by a factor of −0.63 kcal/mol, which means that an extra 0.63 kcal/mol of energy is released while forming a docked complex, which is more stable than that of HCQ.

An increase in binding energy due to the variation of side chain length in CQ and HCQ structures of compounds 6, 9, 11, 14, 15, 16, 17, 18 , and 19 makes them less stable complexes with the protein than with CQ and HCQ. Therefore, the chain length of 4-aminopentyl group is optimal for showing biological activity. When an additional −OH group is added to the ethyl chain of HCQ resulting in 2-[(4aminopentyl)(2-hydroxyethyl)amino]ethan-1-ol side chain as in compound 13, as a result, the binding energy is decreased to −6.0 kcal/mol. Therefore, the −OH group at one terminal ethyl chain only is required for their inhibitory activity. Rearrangement of the amino side chain from the fourth position to fifth position of quinoline ring in HCQ (molecule 33) increased the binding energy by 0.53 kcal/mol. Replacement of quinoline scaffold by isoquinoline scaffold also increased the binding energy by an amount of 0.61 kcal/ mol.

The CQ, HCQ, and its 33 analogs taken for this study were subjected to MLR analysis. Regression models were generated with binding energy as the response and molecular descriptors as the variables. Many models were generated by varying 1−5 variables at a time from a set of 1602 descriptors. The fitness criteria of a model include R 2 ≥ 0.6, R 2 adj ≥ 0.6, Q 2 loo ≥ 0.6, high value for Fisher ratio (F), low regression standard deviation (s), lower correlation between descriptors (K xx ), positive value for δ k , which is the difference in correlation between descriptors and the response (y) value and correlation among descriptors, and smaller root−mean−squared error for training calculation (RMSE tr ). 51, 52 The parameters of the best five models calculated are given in Table 3 .

Considering all the models generated, model 1 has the highest fitness score. However, upon inspecting its external validation parameters RMSE ext = 0.9003 and R 2 ext = 0.3245, the values did not meet the cutoff and hence show poor predictability. For a good prediction, RMSE ext should not vary much from RMSE tr (difference should be less) and R 2 ext should be greater than 0.6.

Model 2 was selected for further studies since it was the best fit model with next highest R 2 and Q 2 loo , which represents the fitness, robustness of the model, and the proximity of actual and predicted binding energy values. The model was built by incorporating the compounds 2, 5, 11, 21, 26, 28, and 30 in the test set and all others in the training set. Thus, the binding energy, the response value (y), can be expressed in terms of the descriptors considered for the generation of model 2 in the form of a linear relationship y = Σm i x i + c where m i is the coefficient obtained for each descriptor in MLR model and x i is the value for the selected descriptor. The internal and external validation parameters determine the feasibility and predictivity of the model 29, 32 (Table S4 ). The linear equation correlating binding energy and descriptors generated by model 2 can be written as shown in eq 1: The root−mean−square error for cross-validation RMSE cv = 0.2296 is greater than RMSE tr = 0.2008, making the model stable and predictive. The criteria for Y scrambling R 2 y scr < 0.2, Q 2 y scr < 0.2, and R 2 y scr > Q 2 y scr implies that there is no correlation between descriptors, and hence, the model is not simply correlated by chance. Model 2 was assessed with external validation parameters R 2 ext > 0.6 and small difference between RMSE ext and RMSE tr and was hence acceptable due to its high predictivity. The predicted values of binding energy for both training and test set of molecules are given in Table  S5 . The plot of observed versus predicted binding energy of training and test set molecules is given in Figure 4 . Substitution of −Cl with −NH 2 in quinoline scaffold in compound 30 led to inaccurate prediction and hence was detected as an outlier.

Model 2 shows that the binding energy is a linear combination of the descriptors SCH-7, nHsNH2, minHBint6, FP402, and KRFP402. SCH-7 is simple seventh order chi chain that is a topological descriptor based on interatomic distances calculated by the bonds between them representing molecular connectivity as a chemical graph like a bond-line formula of chemical structures. It considers the specificity of the structures at a fragment level rather than the whole molecule. 53 The order 7 represents the number of edges in the graph, which indicates the branching. 54, 55 The negative value of its coefficient indicates its negative impact on biological activity. Thus, a low value for SCH-7 is required for a compound to show low binding energy, which makes a restriction on branching. nHsNH2 represents the number of NH 2 groups and amine hydrogens. 56 The positive coefficient of 1.0801 indicates that as the number of NH 2 groups increases, the binding energy decreases, which is in perfect agreement 57 It is an atom-type electro-topological state descriptor in which both electronic and topological functions are considered. 53 There is a positive dependence of minHBint6 with the binding energy because as the chances of H-bond increases, the binding energy will be more negative. FP402 and KRFPC476 are fingerprint descriptors, which account for structure, specific properties, connectivity, or pharmacophores. 58 KRFPC476 accounts for the count of substructures also. 56 Both of these fingerprints are inversely related to binding energy. Designing molecules in such a way that having high values for nHsNH2 and minHBint6 and low values for SCH-7, FP402, and KRFPC476 will help to reduce the binding energy considerably.

Compound 29, primaquine, shows the least binding energy among the CQ analogs. It is the only antimalarial drug that can interfere with a majority of the life cycle stages of the malarial protozoan species. 59, 60 The combined usage of chloroquine and primaquine may mask chloroquine resistance. 61 However, its toxic effect is mainly due to the hemolytic lesions caused by methemoglobin production. 59, 60 It is reported that introduction of a tert-butyl group at C2 position (R1 in scaffold given in Table 4 ) will increase the efficiency and is completely drained of from methemoglobin production as in compound pq1. 59,60

By using the formulated regression Model 2, we predicted the binding energy of some primaquine analogs obtained from the literature and PubChem database and then carried out their molecular docking studies on 3CL pro target to check the inhibitory potency of the ligands, given in Table 4 . Also screened, 10,584 Asinex compounds having antiviral properties using the regression model 2. It was observed that the predicted binding energy of some compounds obtained from Journal of Proteome Research pubs.acs.org/jpr Article the Asinex database differed largely from its actual value, which was obtained from docking analysis. This may be due to their difference in scaffold from that of the compounds used to build the regression model. Even though the predicted binding energy was less accurate, they have furnished better binding energies on docking with the target protein implying the reliability of the model. Thirty molecules that showed lower binding energies were subjected to molecular docking analysis to identify the potent 3CL pro inhibitors (Supplementary Figure  S5) . Binding energies obtained through both the methods and their interactions with the receptor protein are given in Table  5 .

All the primaquine analogs showed lower binding energies than CQ and HCQ. All the compounds except pq5 furnished binding energies lower than that of primaquine. The binding energies calculated using the regression Model 2 are in good agreement with the value obtained from docking analysis. All the compounds exhibited interactions with some of the active site residues including the catalytic diad residues either Cys145 or His41 or both, which make them potent 3CL pro inhibitors. The compounds pq8, pq9, and pq10 exhibited lower binding energies of −8.58, −8.48, and −8.32 kcal/mol. The 2D interaction diagrams of these ligands with the protein are given in Figure 5 . In accordance with the findings based on docking studies, the regression analysis of CQ analogs reveals the presence of one extra −NH 2 bond in their structure.

Inspecting the 2D interaction diagrams of compounds A6, A9, A10, A15, A16, A23, A29, A30 it can be seen that there are unfavorable donor−donor or acceptor-acceptor interactions present between the compounds and target protein which may negatively affect the stability of the docked complexes (Supplementary Figure S6 ). Hence these can be considered as poor inhibitors of 3CL pro protein even though they have binding energies less than HCQ. From these results it is clear that binding energy as well as interactions play decisive roles in determining the druggability of compounds.

Asinex database compounds A12, A5, A13, A11, A20, A19, A18, A4, A14, and A3 obtained lower binding energies among all ( Figure 5 and Figure S7 ). These compounds form interactions with Cys145, His 41, or both along with other active site residues, which are responsible for their inhibitory property. The compounds A12, A5, A13, A11, A19, and A14 form hydrogen bonding interactions with Cys145, which will obstruct the initial deprotonation step in the catalytic cycle of 3CL pro as shown in Figure 1 . The π−π interaction of His 41 with π systems of the compounds A12 and A5 will affect the proton accepting ability of His41. The π−sulfur interactions are present in aromatic rings with −SH group of the compounds A13, A20, A18, A4, A14, and A3 with Cys145 residue. It involves the interaction with the aromatic π orbitals and σ* acceptor orbital of SH bond of Cys145. S−H/π interactions are the important contributors to binding energy. 62 The thiol group of Cys145 is engaged in S−H/π Comparing the structures of the hit compounds pq8, pq9, pq10, and A12 with the rest, it was clear that the high flexibility of the structures owing to more single bonds and less compact rings makes them able to undergo conformational changes at the active site of the target leading to more interactions with Cys145 and His41 along with other active site residues. Hence, these four compounds were selected for further analysis.

Even though docking studies give the best binding poses, they do not account for the conformational changes taking place upon ligand binding as the docking analysis preferably treats protein as rigid. The conformational changes affecting the binding affinity could be analyzed using molecular dynamic simulations. 63 Four top ranked compounds via docking (pq8, pq9, pq10, and A12) were selected for 20 ns molecular dynamics simulation studies to understand the conformational changes in the protein on binding of the ligand. RMSD, RMSF, and Radius of gyration plots were used to evaluate the stability of the complexes. The binding conformations of the compounds at the last frame of simulation are given in Figure  6 .

The root−mean−square deviation (RMSD) plot of apoprotein 3CL pro , pq8, pq9, pq10, and A12 with average RMSDs 1.75, 1.65, 1.93, 2.02, and 1.86 Å, respectively, is given in Figure 7a . It can be inferred that for the apo structure as well as the docked complexes except pq10, RMSD values rose to 2.54 Å and then formed a plateau after 12 ns of simulation, which indicates that the system has reached stable equilibrium state and has been confined in the active site. The RMSD of pq10 requires an extension of stabilization time to reach The root−mean−square fluctuation (RMSF) plot gives the fluctuations of individual residues in the protein backbone. 64 Higher RMSF denotes higher flexibility and vice versa. 65 The RMSF plot (Figure 7b ) showed greater residue flexibility for pq9 and pq10 compared to the unbound form. pq8 and A12 have comparable RMSF values with that of the apo form. The Cys145 and His41 residues involved in the catalytic cycle have fluctuations with greater than 0.7 Å in apoprotein. In pq8 and pq9, it fluctuates more than 3 Å, while in pq10 and A12, there is fluctuation less than 0.6 Å. In all the systems, the active site residues from 140 to 190 have the highest fluctuations.

The structural flexibility or compactness of the protein molecule is analyzed by the radius of gyration. 66 The protein is supposed to maintain a relatively steady value for R g . From inspection of Figure 7c , it can be seen that R g has almost no significant variance with the average R g values of 22 Å. The ligand bound form has lower R g values compared to apoprotein, which indicates that the bound structures are more compact.

Molecular docking analysis of chloroquine and hydroxychloroquine on main protease (3CL pro) protein has shown that HCQ is more efficient than CQ. More potent inhibitors of 3CL pro have been identified through docking studies. The interactions with the active site residues, especially with Cys145 and His41, which are the key residues involved in catalytic diad for proteolysis in the enzyme, are responsible for the activity of these compounds. A regression model with R 2 = 0.9039 and Q 2 = 0.8848 was developed correlating binding energy with the molecular descriptors. This model was used to predict the binding energies of novel molecules. The primaquine derivatives and antiviral compounds obtained from Asinex compound library were screened using the regression model, and the hit compounds thus obtained were further docked with 3CL pro protein. They showed less binding energies and proper interactions with the target, main protease (3CL pro ). Molecular dynamics studies on these top four compounds pq8, pq9, pq10, and A12 revealed their binding stability and conformational changes associated with the protein−ligand complexes. Hence, these four compounds are proposed for validation and further studies as more efficient medication for COVID-19.

The Supporting Information is available free of charge at https://pubs.acs.org/doi/10.1021/acs.jproteome.0c00683.

Binding energy obtained for docking CQ and HCQ with 3CL pro and PfDHFR-TS; structure of CQ analogs with 4-aminoquinoline scaffold; structure of chloroquine analogs with quinoline scaffold, values for internal and external validation criteria, binding energy values of training and test molecules predicted with Model 2; domains and active site of 6LU7 bound with inhibitor N3; 2D interaction diagrams of docked complexes of CQ, HCQ with 3CL pro ; docked poses of CQ and HCQ in active site of 6LU7; 2D interaction diagrams of 29, 25, 22, 8, 24, 1, and 10 with 6LU7; structures of screened compounds from Asinex compound library; 2D interaction diagram of compounds having unfavorable interaction with target; 2D interaction diagrams of Asinex compounds with 6LU7 protein (PDF)

Research and Development on Therapeutic Agents and Vaccines for COVID-19 and Related Human Coronavirus Diseases

Diagnosing COVID-19: The Disease and Tools for Detection

A Review of the 2019 Novel Coronavirus (COVID-19) Based on Current Evidence

Coronavirus Disease 2019: Coronaviruses and Blood Safety

Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2): An Overview of Viral Structure and Host Response

COVID-19) Weekly Epidemiological Update and Weekly Operational Update; WHO

Potential Interventions for Novel Coronavirus in China: A Systematic Review

Structural Basis of SARS-CoV-2 3CLpro and Anti-COVID-19 Drug Discovery from Medicinal Plants

An Overview of Severe Acute Respiratory Syndrome-Coronavirus (SARS-CoV) 3CL Protease Inhibitors: Peptidomimetics and Small Molecule Chemotherapy Department of Medicinal Chemistry, Tokyo University of Pharmacy and Life Sciences

Structure of Mpro from COVID-19 Virus and Discovery of Its Inhibitors

Distinct Roles of Catalytic Cysteine and Histidine in the Protease and Ligase Mechanisms of Human Legumain As Revealed by DFT-Based QM/MM Simulations

Ongoing Clinical Trials for the Management of the COVID-19 Pandemic

Chloroquine and Hydroxychloroquine in the Treatment of COVID-19 with or without Diabetes: A Systematic Search and a Narrative Review with a Special Reference to India and Other Developing Countries

Review: Hydroxychloroquine and Chloroquine for Treatment of SARS-CoV-2 (COVID-19)

New Insights on the Antiviral Effects of Chloroquine against Coronavirus: What to Expect for COVID-19?

Recycling of Chloroquine and Its Hydroxyl Analogue to Face Bacterial, Fungal and Viral Infections in the 21st Century

Hydroxychloroquine, a Less Toxic Derivative of Chloroquine, Is Effective in Inhibiting SARS-CoV-2 Infection in Vitro

Introducing New Antimalarial Analogues of Chloroquine and Amodiaquine: A Narrative Review. Iran

Chloroquine Analogs as Antimalarial Candidates with Potent in Vitro and in Vivo Activity

Chloroquine Analogues in Drug Discovery: New Directions of Uses, Mechanisms of Actions and Toxic Manifestations from Malaria to Multifarious Diseases

Update: Improved Access to Chemical Data

Avogadro: An Advanced Semantic Chemical Editor, Visualization, and Analysis Platform

Avogadro: An Open-Source Molecular Builder and Visualization Tool. Version 1.XX; Avogadro Chemistry

The Protein Data Bank

An Environment for Comparative Protein Modeling. Electrophoresis

Software News and Updates AutoDock4 and AutoDockTools4: Automated Docking with Selective Receptor Flexibility

A New Software for the Development, Analysis, and Validation of QSAR MLR Models

Insubria Datasets and New QSAR/QSPR Models for Environmental Pollutants in QSARINS

PaDEL-Descriptor: An Open Source Software to Calculate Molecular Descriptors and Fingerprints

Comparative QSAR Model Generation Using Pyrazole Derivatives for Screening Janus Kinase-1 Inhibitors

Quantitative Structure Activity Relationship Study of Anticonvulsant Activity of α _ Substituted Acetamido-N-Benzylacetamide Derivatives

Scalable Molecular Dynamics on CPU and GPU Architectures with NAMD Scalable Molecular Dynamics on CPU and GPU Architectures with NAMD

Ligand Reader & Modeler for CHARMM Force Field Generation of Small Molecules

VMD: Visual Molecular Dynamics

Research Paper State-of-the-Art Tools Unveil Potent Drug Targets amongst Clinically Approved Drugs to Inhibit Helicase in SARS-CoV-2

Physicochemical Characterization and Functional Analysis of Some Snake Venom Toxin Proteins and Related Non-Toxin Proteins of Other Chordates Abstract: Background: Bioinformation

Genetic Involvement of Interleukin 4 for Asthma and Identification of Potential Phytochemical Scaffold Through Molecular Docking Studies

Understanding the Relationship between the Primary Structure of Proteins and Its Propensity to Be Soluble on Overexpression in Escherichia Coli

Structure of Mpro from COVID-19 Virus and Discovery of Its Inhibitors

Structure -Activity Relationship and Comparative Docking Studies for Cycloguanil Analogs as PfDHFR-TS Inhibitors

Chalcone Analogue as Potent Anti-Malarial Compounds against Plasmodium Falciparum: Synthesis, Biological Evaluation, and Docking Simulation Study. Asian Pac

Pi-Interactions in Proteins

PH-Dependent Entry of Severe Acute Respiratory Syndrome Coronavirus Is Mediated by the Spike Glycoprotein and Enhanced by Dendritic Cell Transfer through DC-SIGN

Primaquine or Other 8-Aminoquinolines for Reducing Plasmodium Falciparum Transmission

Primaquine Therapy for Malaria

Analysis for Antioxidant Activity of Dipicolinic Acid Derivatives

QSAR Study of Anti-Human African Trypanosomiasis Activity for 2-Phenylimidazopyridines Derivatives Using DFT and Lipinski's Descriptors. Heliyon

QSAR with Electro Topological State Atom Index: Antialirenergic Activi Ty of N, N-Dimethyl-2-Bromo-2-Phenylethylamines

The Molecular Connectivity Chi Indexes and Kappa Shape Indexes in Structure-Property Modeling. Reviews in Computational Chemistry

) Yap, C. W. PaDEL-Descriptor: An open source software to calculate molecular descriptors and fingerprints

Investigation of 5,6-Dihydro-2-Pyrones Derivatives as Potent Anti-HIV Agents Inhibitors

Quantitative Regression Models for the Prediction of Chemical Properties by an Efficient Workflow

Synthesis and Blood-Schizontocidal Antimalarial Activities of 2-Substituted/2, 5-Disubstituted-8-Quinolinamines and Some of Their Amino Acid Conjugates

Discovery of a Bulky 2-Tert -Butyl Group Containing Primaquine Analogue That Exhibits Potent Blood-Schizontocidal Antimalarial Activities and Complete Elimination of Methemoglobin Toxicity

Blood Stage Antimalarial Efficacy of Primaquine in Plasmodium Vivax Malaria

Insights into Thiol-Aromatic Interactions: A Stereoelectronic Basis for S -H/π Interactions

Docking Studies and Molecular Dynamics Simulations of the Binding Characteristics of Waldiomycin and Its Methyl Ester Analog to Staphylococcus Aureus Histidine Kinase

Understanding the Structural Features of JAK2 Inhibitors: A Combined 3D-QSAR, DFT and Molecular Dynamics Study

Molecular Docking and Molecular Dynamics to Identify a Novel Human Immunodeficiency Virus Inhibitor from Alkaloids of Toddalia Asiatica

Molecular Docking and Molecular Dynamics Simulation Studies to Identify Potent AURKA Inhibitors: Assessing the Performance of Density Functional Theory, MM-GBSA and Mass Action Kinetics Calculations

The authors declare no competing financial interest.

The authors gratefully acknowledge generous funding from Kerala State Council for Science, Technology and Environment (KSCSTE) for the Junior Research Fellowship (No. 414/2020/KSCSTE).