Journal of Cellular and Molecular Pharmacology
Open Access

Our Group organises 3000+ Global Conferenceseries Events every year across USA, Europe & Asia with support from 1000 more scientific Societies and Publishes 700+ Open Access Journals which contains over 50000 eminent personalities, reputed scientists as editorial board members.

Open Access Journals gaining more Readers and Citations
700 Journals and 15,000,000 Readers Each Journal is getting 25,000+ Readers

This Readership is 10 times more when compared to other Subscription Journals (Source: Google Analytics)
  • Review Article   
  • J Cell Mol Pharmacol 2022, Vol 6(5)Vol 6(5)
  • DOI: 10.4172/jcmp.1000135

Utilizing Virtual Screening and Structure-Based Drug Design to Repurpose Approved Medications as Possible Inhibitors of the 3CL-Protease of the SARS-Cov-2

Colin W G Fishwick*
School of Chemistry, University of Leeds, Leeds LS2 9JT, UK
*Corresponding Author: Colin W G Fishwick, School of Chemistry, University of Leeds, Leeds LS2 9JT, UK, Email: c.w.g.fiick@leeds.ac.uk

Received: 03-Sep-2022 / Manuscript No. jcmp-22-73438 / Editor assigned: 06-Sep-2022 / PreQC No. jcmp-22-73438 / Reviewed: 20-Sep-2022 / QC No. jcmp-22-73438 / Revised: 22-Sep-2022 / Published Date: 28-Sep-2022 DOI: 10.4172/jcmp.1000135

Abstract

Since there is currently no approved vaccine or small molecule therapeutic on the market, the recent global outbreak of SARS-CoV-2, which has nearly 15 million infected people and more than 600,000 fatalities (accessed on 20th July 2020), poses a significant challenge to all nations and societies. Nearly 1500 investigations are actively enrolling patients by invitation for clinical research against COVID-19 because to the urgent need for a causative therapy

The main protease of Coronaviruses, chymotrypsin-like protease (3CLpro), processes the large polyprotein 1ab releasing several further enzymes that are crucial for viral replication. Moreover, 3CLpro is unique for Coronaviruses and not found in higher organisms. This predestines 3CLpro as most attractive target for the development of antiinfective agents against SARS-CoV-2 and related Coronaviruses [1]. Consequently, several inhibitors of 3CLpro were developed mostly during the last 17 years right after the first wave of infection caused by the SARS-CoV-1. However, no experimental compound was developed further and reached the market. Most compounds were designed as covalent inactivators that react with the catalytic Cys145.

These guidelines, which are sometimes referred to as features, are made up of spheres where a heavy ligand atom with specific properties, such as hydrogen donor or acceptor properties, is meant to be positioned in order to make a specific interaction with a neighbouring amino acid. It is best to simulate the binding pocket using a sequence of overlapping excluded volume spheres in order to reduce the likelihood of receiving an excessive number of false positive hits for branched or bulky molecules that meet the requirements but would interfere with the receptor protein. These were produced utilising a radius of 2 from particular pocket amino acids [2]. Hits are not counted for molecules whose structures conflict with those of the excluded volume spheres. Using a pharmacophore model, one Very quick search for hits that satisfy all or some of the feature requirements across big 3D structure databases. The search process should very effectively enrich those chemical entities that provide high docking scores and have a higher likelihood of interacting with the relevant target protein if the pharmacophore model is adequately described. The MOE pharmacophore editor was used to launch the pharmacophore search [3]. All of the hits from the pharmacophore search were put through the previously mentioned follow-up virtual screening process using MOE to evaluate their binding affinity.

Introduction

This is reflected in the majority of the 3CLpro enzyme x-ray structures from different Coronaviruses. As opposed to covalent inhibitors, which are common, small reversible inhibitors are typically chosen since they have fewer toxicities and side effects. SARS-3CLpro CoV-2 recently showed two crystal structures, both of which were bound to covalent inhibitors). The crystal structures of the 3CLprocomplexes’ PDB-IDs 6LU7 and 6Y2F show total overlap with an RMSD-value of 0.48 over more than 300 amino acids. The most significant crystal structure of 3CLpro of SARS-CoV-1 with a noncovalent inhibitor is found in PDB-ID: 3V3M.It’s intriguing how many drugs that have been authorised for use in other ailments also demonstrated success against coronaviruses. The antiquated malaria drug hydroxychloroquine, which has shown effectiveness in clinical trials and appeared to be effective in inhibiting the multiplication of SARS-CoV-2 in vitro, is one of the active ingredients.Gao and colleagues 2020. These findings triggered a passionate debate about the hydroxylchloroquine’s possible benefits for COVID-19 sufferers [4-15]. A significant retrospective study, however, discovered that COVID-19 patients who took chloroquine or hydroxychloroquine developed severe dysrhythmias and received no benefit.

Subjective Heading

In addition, a very recent FDA review of the safety concerns associated with the use of these chloroquine derivates highlights the connection between treatment and serious heart rhythm problems as well as other safety concerns, such as blood and lymph system disorders, kidney injuries, and liver failure (https://www.fda.gov/ drugs/drug-safety-and-availability/fda-cautions-against-usehydroxychloroquine- or-chloroqu Because of this, the official NIH treatment recommendation (accessed on July 20, 2020) advises against using chloroquine or hydroxychloroquine to treat COVID-19, with the exception of in clinical setting up tests. A corticosteroid called dexamethasone is used to treat a wide range of inflammatory illnesses, including allergy disorders, ulcerative colitis, arthritis, lupus, psoriasis, and respiratory problems. Dexamethasone is suggested for the treatment of COVID-19 patients in the NIH’s actual COVID-19 treatment guideline due to lower mortality rates compared to patients who received conventional care.

In cell culture, other medications like mefloquine, selamectin, and cepharanthin also demonstrated positive benefits (Fan et al., 2020). The National Health Commission in China also advises using ribavirin and interferon together to treat COVID-19 because to its impact on MERS-CoV.Remdesivir, a broad-spectrum antiviral medication, is currently being studied in a clinical trial to treat COVID-19. Remdesivir was previously tested on humans who had the Ebola virus sickness. Initial findings in MERS-CoV-infected animal models were encouraging. Additionally, the clinically validated camostat mesilate has demonstrated that inhibition of the serine protease prevents SARSCoV- 2 from entering cells.

Discussion

The MOE 2019 programme (Chemical Computing Group, Montreal, Canada) and Autodock Vina (The Scripps Research Institute, La Jolla, USA) were used for two virtual screens.The latter was integrated into the PyRx environment’s free version, which is accessible at https://sourceforge.net/projects/pyrx/. The entire AMBER14:EHT force field was utilised for the MOE-based virtual screen. The RCSB Protein Data Bank provided the SARS-CoV-2 3CLpro crystal structure (PDB-ID: 6LU7). The covalently attached ligand to the catalytic Cys145 was eliminated, and the Cys145 was then made available as a free thiol. After that, the free target protein underwent MOE’s QuickPrep method, which included repairs for any missing atoms, misaligned geometries, or other 3D protonation, removing water molecules more than 4.5 from any receptor or ligand atom, and crystallographic aberrations. The binding pocket was established using the expelled covalent ligand. A MOE-specific database was imported and the SD-file containing the 3D structures of approved medications was chosen as the ligand in the docking setup for virtual screening. The stiff receptor and the triangular matcher were chosen as additional factors. Ten docking poses per ligand molecule were evaluated using the London dG score. The screening results were exported as SD files, and Osiris Datawarrior was used to integrate several findings for a single chemical structure to create a result file with distinct chemical structures. maintaining various docking scores for each molecule. The crystal structure of 3CLpro was likewise liberated from the covalent ligand for the virtual screen with Autodock Vina, and it was afterwards put through the Dockprep technique used by the UCSF Chimera programme (University of California, USA 2004). The procedure calls for the removal of solvent, the use of a rotamer library to replace missing side chains, protonation, and the insertion of charges.

PyRx was used to load an SD-file containing the 3D structures of the ligand library and convert it to the necessary pdbqt-format. The implemented Vina wizard was then used to pick the protein target and the compound library as numerous ligands in the following phase. An area cuboid was employed to specify the active site pocket before beginning the virtual screening process. The exhaustiveness, which measures docking precision, was set at 8 in the last stage. The Vina wizard then began the virtual screening process. Again, the data were acquired in pdbqt format, which needed to be changed for additional analysis into SD format.

Using MOE 2019 software, structural data visualisation and molecular docking were carried out (Chemical Computing Group, Montreal, Canada). RCSB Protein Data Bank provided the 3CLpro of SARS-CoV-2 crystal structure (PDB-ID: 6LU7). The structure file was imported into the computer and then underwent the previously described 3D protonation and structure preparation. The installed Amber14:EHT force field was used to determine the partial charges of all protein and ligand atoms. The ligand was placed in the binding site using the triangle matcher during molecular docking, and the London dG scoring function was used to rate the results. The best 50 poses were passed to the refinement and energy minimization in the pocket using the induced fit method and then rescored with the GBVI/ WSA dG scoring function. Best poses were further refined by energy minimization of all amino acids in a radius of 10 Å around the ligand.

In order to facilitate the investigation of protein-ligand interactions, the complex structures of 3CLpro with the highest scoring ligands were superposed and subjected to the Quickprep process implemented in MOE using the AMBER14:EHT force field once more. Several docked protein-ligand complexes with ideal binding postures were used as the foundation for the generation and analysis of protein ligand interaction fingerprinting (PLIF) utilising the corresponding MOE tool. A panel displaying the precise connections of each ligand to a variety of proteinligand complexes that the PLIFs have selected as the most significant relationships amino acids that are close to the binding pocket. Most crucially, a pharmacophore model encapsulating fundamental requirements that a possible ligand of the corresponding binding pocket must satisfy can be created using the PLIFs and frequency of proteinligand contacts. These guidelines, which are sometimes referred to as features, are made up of spheres where a heavy ligand atom with specific properties, such as hydrogen donor or acceptor properties, is meant to be positioned in order to make a specific interaction with a neighbouring amino acid. It is best to simulate the binding pocket using a sequence of overlapping excluded volume spheres in order to reduce the likelihood of receiving an excessive number of false positive hits for branched or bulky molecules that meet the requirements but would interfere with the receptor protein. These were produced utilising a radius of 2 from particular pocket amino acids. Hits are not counted for molecules whose structures conflict with those of the excluded volume spheres. Large 3D structure databases can be quickly searched for matches that satisfy all or some of the feature requirements using a pharmacophore model. The search process should very effectively enrich those chemical entities that provide high docking scores and have a higher likelihood of interacting with the relevant target protein if the pharmacophore model is adequately described. The MOE pharmacophore editor was used to launch the pharmacophore search. All of the hits from the pharmacophore search were put through the previously mentioned follow-up virtual screening process using MOE to evaluate their binding affinity.

Then, a comparison and activity Cliff Analysis was used to visualise the chemical landscape, find clusters and singletons with high docking scores, and group related molecules together on a 2D surface.There are four significant clusters that stand out, including sizable flavonoid, tetracycline, aminoglycoside, and anthracycline clusters (Fig. 2B). Quercetin, oxytetracycline, kanamycin, and doxorubicin serve as illustrative medications for these clusters, respectively. There are other singletons or groups of two with high scores, such as raloxifen. It is in great agreement with the very recent publication by Jo et al., who provide experimental proof that flavonoids are in fact inhibitors, that numerous flavonoids are among the hits with the best docking scores. that 3CLpro.It should be mentioned that 6.5% of the FDA-approved medications and 4.5% of the flavonoids quercetin, rutin, homoorientin, all of which are flavonoids, eltrombopag, and doxorubicin, struck the virtual screen PAINS patterns. In order to solve the issue of frequent hitters in experimental high throughput screening campaigns-often false positive hits-Baell and Holloway created the idea of PAINS. However, only six patented assays assessing protein-protein interaction (PPI) inhibition using the AlphaScreen detection technique were used to test the crucial substructural components of electronic PAINS filters. warn against using PAINS filters in a careless manner to identify and classify substances with liabilities and suggest that only orthogonal experiments should be used to reach such findings.All compounds were taken further to detailed docking analysis to elucidate the potential molecular interactions with 3Clpro, because a wide variety of orthogonal assays despite the fact that some of the approved drug molecules contain crucial substructures like labile ester (salvianolic acid B) or possibly redox active groups like electron rich scaffolds, tests have been conducted to demonstrate biological activity and safety of the molecules prior to drug approval (polyphenols). The virtual screen’s best hits were subjected to a more thorough docking technique that sampled 50 docking postures utilising a secondary GBVI/WSA dG score that was introduced in MOE. Table 1 summarises the findings of the compounds that showed the greatest promise. Curiously, most high Antioxidants called flavonoids are also used medicinally as anti-cancer or antibiotics. The two medicines with the greatest docking scores, kanamycin and salvianolic acid B, display comparable binding postures by occupying the S1- and S2-subpockets Both ligands form H-bonds with the catalytic Cys145 despite the fact that kanamycin interacts with the thiolate moiety and salvianolic acid B with the backbone nitrogen of Cys145. The most effective protein-ligand interactions with 3CLpro for typical ligands are shown in Fig. S2, and they are further discussed in the section.

The development of a pharmacophore model was made possible by the thorough analysis of the key chemical interactions involving ligands and amino acids that border the active site pocket. A simplified version of this model describes the chemical properties in terms of molecular interactions, such as H-bond donator or acceptor, and their 3D layout within the receptor protein’s active site pocket. By overlapping excluded volume spheres with the protein pocket shape, the pharmacophore model is finished. If any of a ligand’s atom centres cross across an excluded volume, the ligand will not match. A pharmacophore model for 3CLpro was created using the results of the PLIF-analysis mentioned above. A heavy atom’s position that can act as an H-bond donor, acceptor, or metal ligator is one of the model’s five essential characteristics, or features.The pharmacophore model outperforms conventional virtual screens for screening larger compound libraries for chemical entities that satisfy the requirements of the chemically and spatially defined model characteristics, which rely on quick rigid model docking techniques.

Conclusion

This study was prompted by the present SARS-CoV-2 epidemic to quickly find currently approved medications that could be modified to target 3CLpro, a cysteine protease crucial for Coronavirus replication. Several scaffolds with a high potential to inhibit 3CLpro were discovered by a virtual screening method, including flavonoids, which were recently characterised as experimentally proven inhibitors of this enzyme, antibiotics, and anticancer drugs. most optimistic Oxytetracycline, naringin, kanamycin, cefpiramide, salvianolic acid-B, teniposide, etoposide, and doxorubicin were among the drugs that were impacted. A further pharmacophore search that screened 7.2 million compounds from the ZINC15 collection greatly expanded the tested chemical space, to well-known chemical groups like flavonoids. Furthermore, the detailed examination of the critical interactions between ligands and the amino acids in the active site pocket opened up new avenues for creating and enhancing 3CLpr inhibitors.

Acknowledgement

I would like to thank my Professor for his support and encouragement.

Conflict of Interest

The authors declare that they are no conflict of interest.

References

  1. Baell JB,Holloway GA(2010) New substructure filters for removal of pan assay interference compounds (PAINS) from screening libraries and for their exclusion in bioassays.J Med Chem: 2719-2740.
  2. Indexed at, Google Scholar, Crossref

  3. Bajorath J, Peltason L,Wawer M(2009) Navigating structure-activity landscapes.Drug Discov Today 14 (13–14):698-705.
  4. Indexed at, Google Scholar, Crossref

  5.  Berry M, Fielding BC, Gamieldien J(2015) Potential broad Spectrum inhibitors of the coronavirus 3CLpro a virtual screening and structure-based drug design. studyViruses7 (12):6642-6660.
  6. Indexed at, Google Scholar, Crossref

  7.  Capuzzi SJ, . Muratov EN ,  TropshaPhantomA(2017) Problems with the Utility of Alerts for Pan-Assay INterference Compound. J Chem  Inf Model 57 (3):417-427.
  8. Indexed at, Google Scholar, Crossref

  9. Cortegiani A, Ingoglia G,  Ippolito M, Giarratano A(2020) A systematic review on the efficacy and safety of chloroquine for the treatment of COVID-19.J Crit. Care 57 (3):417-427.
  10. Indexed at, Google Scholar, Crossref

  11. Dong E ,  Du EL,  Gardner  L (2020) An interactive web-based dashboard to track COVID-19 in real time Lancet. Infect Dis 7 (12):6642-6660.
  12. Indexed at, Google Scholar, Crossref

  13.  Fan HH,  Wang LQ (2020) Repurposing of clinically approved drugs for treatment of coronavirus disease 2019 in a 2019-novel coronavirus.model Chin Med J.
  14. Indexed at, Google Scholar, Crossref

  15. Gao J, Tian Z, Yan X(2020) Breakthrough Chloroquine phosphate has shown apparent efficacy in treatment of COVID-19 associated pneumonia in clinical studies. Biosci Trends 14 (1):72-73.
  16. Indexed at, Google Scholar, Crossref

  17. Flexner C (1998) HIV-protease inhibitors N Engl J  Med 338:1281-1292.
  18. Indexed at, Google Scholar, Crossref

  19. Ghosh AK,  Osswald HL (2016) Prato Recent progress in the development of HIV-1 protease inhibitors for the treatment of HIV/AIDS. J Med Chem 59 (11):5172-5208.
  20. Indexed at, Google Scholar, Crossref

  21.  Lv Z,  Chu Y,Wang  Y(2015) HIV protease inhibitors a review of molecular selectivity and toxicity. HIV AIDS Res Palliat Care 7:95-104.
  22. Indexed at, Google Scholar, Crossref

  23. Wlodawer A, Vondrasek J  (1918) Inhibitors of HIV-1 protease  a major success of structure-assisted drug design. Annu Rev  Biophys  Biomol Struct 27:249-284.
  24. Indexed at, Google Scholar, Crossref

  25. Baell JB,Holloway GA(2010) New substructure filters for removal of pan assay interference compounds (PAINS) from screening libraries and for their exclusion in bioassays.J Med Chem 2719-2740.
  26. Indexed at, Google Scholar, Crossref

  27. Paterson DL, Swindells S, Mohr J, Brester M (2000) Adherence to protease inhibitor therapy and outcomes in patients with HIV infection.Ann Intern Med  133 (1):21-30.
  28. Indexed at, Google Scholar, Cross Ref

  29. Price GW, Gould PS, Mars A(2014)Use of freely available and open source tools for in silico screening in chemical biology. J Chem Educ 91 (4):602-604.
  30. Indexed at, Google Scholar, Crossref

Citation: Fishwick CWG (2022) Utilizing Virtual Screening and Structure-Based Drug Design to Repurpose Approved Medications as Possible Inhibitors of the 3CL-Protease of the SARS-Cov-2. J Cell Mol Pharmacol 6: 135. DOI: 10.4172/jcmp.1000135

Copyright: © 2022 Fishwick CWG. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Top