Molecular screening of glycyrrhizin-based inhibitors against ACE2 host receptor of SARS-CoV-2.
Sajjad AhmadYasir WaheedAsma AbroSumra Wajid AbbasiSaba IsmailPublished in: Journal of molecular modeling (2021)
The interaction between SARS-CoV-2 Spike protein and angiotensin-converting enzyme 2 (ACE2) is essential to viral attachment and the subsequent fusion process. Interfering with this event represents an attractive avenue for the development of therapeutics and vaccine development. Here, a hybrid approach of ligand- and structure-based virtual screening techniques were employed to disclose similar analogues of a reported antiviral phytochemical, glycyrrhizin, targeting the blockade of ACE2 interaction with the SARS-CoV-2 Spike. A ligand-based similarity search using a stringent cut-off revealed 40 FDA-approved compounds in DrugBank. These filtered hits were screened against ACE2 using a blind docking approach to determine the natural binding tendency of the compounds with ACE2. Three compounds, deslanoside, digitoxin, and digoxin, were reported to show strong binding with ACE2. These compounds bind at the H1-H2 binding pocket, in a manner similar to that of glycyrrhizin which was used as a control. To achieve consistency in the docking results, docking calculations were performed via two sets of docking software that predicted binding energy as ACE2-Deslanoside (AutoDock, -10.3 kcal/mol and DockThor, -9.53 kcal/mol), ACE2-Digitoxin (AutoDock, -10.6 kcal/mol and DockThor, -8.84 kcal/mol), and ACE2-Digoxin (AutoDock, -10.6 kcal/mol and DockThor, -8.81 kcal/mol). The docking results were validated by running molecular simulations in aqueous solution that demonstrated the stability of ACE2 with no major conformational changes in the ligand original binding mode (~ 2 Å average RMSD). Binding interactions remained quite stable with an increased potential for getting stronger as the simulation proceeded. MMGB/PBSA binding free energies were also estimated and these supported the high stability of the complexes compared to the control (~ -50 kcal/mol net MMGB/PBSA binding energy versus ~ -30 kcal/mol). Collectively, the data demonstrated that the compounds shortlisted in this study might be subjected to experimental evaluation to uncover their real blockade capacity of SARS-CoV-2 host ACE2 receptor.
Keyphrases
- angiotensin converting enzyme
- sars cov
- angiotensin ii
- molecular dynamics
- molecular dynamics simulations
- binding protein
- protein protein
- dna binding
- density functional theory
- respiratory syndrome coronavirus
- small molecule
- single molecule
- magnetic resonance
- drug delivery
- risk assessment
- single cell
- computed tomography
- aqueous solution
- climate change
- coronavirus disease
- monte carlo