Submitted: 04 June 2020 | Approved: 02 July 2020 | Published: 03 July 2020
How to cite this article: Kumari V, Kanwar SS. In silico comparative analysis of HIV protease inhibitors effect on 2019-nCoV coronavirus 3CLpro. Ann Proteom Bioinform. 2020; 4: 012-014.
Copyright: © 2020 Kumari V, et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Keywords: 2019-nCoV; Protease inhibitors; Pneumonia; 3CLpro; Interaction analysis; Ligands
The novel coronavirus 2019-nCoV has become a bane to mankind and spread worldwide and infected many people. Thus, there is an urgent need of a cure for the severe pneumonia disease caused by this virus. In this study, In silico comparative analysis has been done for HIV protease inhibitors on coronavirus 3CLpro protein which has shown the major interactions and common amino acid residues involved in interactions. The amino acid interaction analysis has revealed two amino acids ARG4, LYS5 to be the major amino acids targets among selected ligands. The binding energy analysis has also revealed Cobicistat as one of these best suited ligand for 3CLpro.
The outbreak of novel coronavirus has caused major loss of human lives and many lives are still at risk. 2019-nCoV is an enveloped, positive-sense, single-stranded RNA beta-coronavirus . Corona viruses can infect respiratory, gastrointestinal, hepatic and central nervous system of human , livestock, birds, bats , mouse, and many other wild animals . Similar to SARS and MERS, non-structural proteins encoded by 2019-nCoV genome are 3-chymotrypsin-like protease (3CLpro), papain-like protease (PLpro), RNA-dependent RNA polymerase (RdRp) and helicase, structural proteins spike glycoprotein . 3CLpro protein is considered as one of the main target for 2019-nCoV and the urgent need of cure has made it more prominent. In this study, the 3CLpro has been considered as a target and HIV protease inhibitors (ligands) were comparatively analyzed in order to reveal the interactions and important amino acids involved in such interactions, so as to make it easier to identify drug-targets.
The three dimensional crystal structure of 2019-nCoV coronavirus protein 3CLpro (Pdb id: 6lu7) was obtained from RCSB PDB (https://www.rcsb.org/structure/). The ligand groups and water molecules already present were removed out in order to make this structure suitable for further interactions with other ligands.
The HIV protease inhibitors has been scanned for the selected coronavirus protein structure. Three ligands Lopinavir , Cobicistat , Amprenavir, Atazanavir, Indinavir and Saquinavir  were selected from the literature review. The chemical structures for these ligands were obtained from PubChem in SDF format and were converted to PDB format for further interactions analysis.
HIV protease inhibitors and 3CLpro interactions
The selected ligands were docked to coronavirus protein 3CLpro using Autodock tool 1.5.6 . Autodock is an automated docking tool to predict how the small molecule or substrate bind to a specific target protein. The protein and ligand PDB formatted files were converted to PDBQT formats followed by docking using Lamarckian genetic algorithm and a total of ten conformations for each ligand docking were observed with their binding energies to find out the best lead compound among them. The protein ligand interactions were obtained followed by the identification of interacting amino acid residues of coronavirus protein 3CLpro to chemical groups with their respective ligands.
The docking analysis of 3CLpro and target ligands revealed the interaction types in three ligands. The docking studies involving 3CLpro with ligands such as Lopinavir, Cobicistat, Amprenavir, Atazanavir, Indinavir and Saquinavir were done to seek greater insight about the amino acid residues of MATE interacting with the specified ligands (Figure 1).
The interaction analysis revealed the occurrence of two amino acid residues ARG4 and LYS5 in interactions with selected ligands which conferred these two amino acid residues to be major drug target sites (Table 1).
The binding energies for ligands and 3CLpro interactions were observed and analyzed comparatively. The lowest binding energy among the selected drugs; Cobicistat, Lopinavir, Amprenavir, Atazanavir, Indinavir and Saquinavir was found to be associated with Cobicistat which conferred it as one of the best suited ligand for 3CLpro (Table 2). However, Cobicistat lacks the enzyme inducing activity and thus requires the co-administered drugs with close monitoring of dose adjustments while treating the patients.
The novel coronavirus 2019-nCoV has caused major loss to human lives, socially and economically as well. The urgent cure is the need of hour for this disease. The results observed from present studies has revealed the ligand with lowest binding energy ‘Cobicistat’ among the selected ligands which confers it to be one of the best suited ligand for 3CLpro. Cobicistat and Ritonavir are structural analogues of each other but Ritonavir has major limitations like poor water solubility which corresponds to difficulty in its co-formulation and drug interactions. Furthermore, the switching between Cobicistat-based drug and Ritonavir-based drug treatment may produce significant problems among the patients. In order to avoid these limitations, only cobicistat has been taken under consideration in this study. The common amino acids ARG4 & LYS5 in the docking studies confirmed these two amino acids as main targets. The amino acid residues involved in interaction with novel coronavirus protein 3CLpro also revealed the common interactions which shall further help scientists to find out the new target sites and lead compounds for novel coronavirus 2019-nCoV.
This work has been funded by the Department of Science & Technology, New Delhi, under a DST-INSPIRE Fellowship [IF180179] awarded to one of the authors (VK). The authors are thankful to DST and DBT, Ministry of Science & Technology, New Delhi for continuous financial support to the Bioinformatics Centre, Himachal Pradesh University, Shimla (India).
- Li G, De Clercq E. Therapeutic options for the 2019 novel coronavirus (2019-nCoV). Nat Rev Drug Discov. 2020; 19: 149-150. PubMed: https://pubmed.ncbi.nlm.nih.gov/32127666/
- Ge XY, Li JL, Yang XL, Chmura AA, Zhu G, et al. Isolation and characterization of a bat SARS‐like coronavirus that uses the ACE2 receptor. Nature. 2013; 503: 535‐538. PubMed: https://pubmed.ncbi.nlm.nih.gov/24172901/
- Wang LF, Shi Z, Zhang S, Field H, Daszak P, et al. Review of bats and SARS. Emerg Infect Dis. 2006; 12: 1834‐1840. PubMed: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3291347/
- Chen Y, Guo D. Molecular mechanisms of coronavirus RNA capping and methylation. Virol Sin. 2016; 31: 3‐11. PubMed: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC7091378/
- Kim UJ, Won EJ, Kee SJ, Jung S, Jang HC. Combination therapy with lopinavir/ritonavir, ribavirin and interferon-alpha for Middle East respiratory syndrome. Antivir Ther. 2016; 21: 455-459. PubMed: https://pubmed.ncbi.nlm.nih.gov/26492219
- Pang J, Wang MX, Ang IYH, Tan SHX, Lewis RF, et al, Potential Rapid Diagnostics, Vaccine and Therapeutics for 2019 Novel Coronavirus (2019-nCoV): A Systematic Review J Clin Med. 2016; 9: 623. PubMed: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC7141113/
- Lv Z, Chu Y, Wang Y. HIV protease inhibitors: a review of molecular selectivity and toxicity. HIV AIDS. 2015;7: 95–104. PubMed: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4396582/
- Morris GM, Huey R, Lindstrom W, Sanner MF, Belew RK, et al. AutoDock4 and AutoDockTools4: Automated Docking with Selective Receptor Flexibility. J Comput Chem. 2009; 30: 2785-2791. PubMed: https://pubmed.ncbi.nlm.nih.gov/19399780/