Identification of natural compounds as SARS-CoV-2 inhibitors via molecular docking and molecular dynamic simulation

Background

Base mutations increase the contagiousness and transmissibility of the Delta and Lambda strains and lead to the severity of the COVID-19 pandemic. Molecular docking and molecular dynamics (MD) simulations are frequently used for drug discovery and relocation. Small molecular compounds from Chinese herbs have an inhibitory effect on the virus. Therefore, this study used computational simulations to investigate the effects of small molecular compounds on the spike (S) protein and the binding between them and angiotensin-converting enzyme 2 (ACE2) receptors.

Methods

In this study, molecular docking, MD simulation, and protein-protein analysis were used to explore the medicinal target inhibition of Chinese herbal medicinal plant chemicals on SARS-CoV-2. 12,978 phytochemicals were screened against S proteins of SARS-CoV-2 Lambda and Delta mutants.

Results

Molecular docking showed that 65.61% and 65.28% of the compounds had the relatively stable binding ability to the S protein of Lambda and Delta mutants (docking score ≤ -6). The top five compounds with binding energy with Lambda and Delta mutants were clematichinenoside AR2 (-9.7), atratoglaucoside,b (-9.5), physalin b (-9.5), atratoglaucoside, a (-9.4), Ochnaflavone (-9.3) and neo-przewaquinone a (-10), Wikstrosin (-9.7), xilingsaponin A (-9.6), ardisianoside G (-9.6), and 23-epi-26-deoxyactein (-9.6), respectively. Four compounds (Casuarictin, Heterophylliin D, Protohypericin, and Glansrin B) could interact with S protein mutation sites of Lambda and Delta mutants, respectively, and MD simulation results showed that four plant chemicals and spike protein have good energy stable complex formation ability. In addition, protein-protein docking was carried out to evaluate the changes in ACE2 binding ability caused by the formation of four plant chemicals and S protein complexes. The analysis showed that the binding of four plant chemicals to the S protein could reduce the stability of the binding to ACE2, thereby reducing the replication ability of the virus.

Conclusion

To sum up, the study concluded that four phytochemicals (Casuarictin, Heterophylliin D, Protohypericin, and Glansrin B) had significant effects on the binding sites of the SARS-CoV-2 S protein. This study needs further in vitro and in vivo experimental validation of these major phytochemicals to assess their potential anti-SARS-CoV-2. Graphical abstract.

The Development of Pharmacophore Models for the Search of New Natural Inhibitors of SARS-CoV-2 Spike RBD-ACE2 Binding Interface

To date, some succeeding variants of SARS-CoV-2 have become more contagious. This virus is known to enter human cells by binding the receptor-binding domain (RBD) of spike protein with the angiotensin-converting enzyme 2 (ACE2), the latter being a membrane protein that regulates the renin-angiotensin system. Since the host cell receptor plays a critical role in viral entry, inhibition of the RBD-ACE2 complex is a promising strategy for preventing COVID-19 infection. In the present communication, we propose and utilize an approach based on the generation of a complex of pharmacophore models and subsequent Induced Fit Docking (IFD) to identify potential inhibitors of the main binding sites of the Omicron SARS-CoV-2 RBD(S1)-ACE2 complex (PDB ID: 7T9L) among a number of natural products of various types and origins. Several natural compounds have been found to provide a high affinity for the receptor of interest. It is expected that the present results will stimulate further research aimed at the development of specialized drugs against this virus.

A multi-targeting drug design strategy for identifying potent anti-SARS-CoV-2 inhibitors

The COVID-19, caused by SARS-CoV-2, is threatening public health, and there is no effective treatment. In this study, we have implemented a multi-targeted anti-viral drug design strategy to discover highly potent SARS-CoV-2 inhibitors, which simultaneously act on the host ribosome, viral RNA as well as RNA-dependent RNA polymerases, and nucleocapsid protein of the virus, to impair viral translation, frameshifting, replication, and assembly. Driven by this strategy, three alkaloids, including lycorine, emetine, and cephaeline, were discovered to inhibit SARS-CoV-2 with EC50 values of low nanomolar levels potently. The findings in this work demonstrate the feasibility of this multi-targeting drug design strategy and provide a rationale for designing more potent anti-virus drugs.

Cheminformatics-Based Identification of Potential Novel Anti-SARS-CoV-2 Natural Compounds of African Origin

The coronavirus disease 2019 (COVID-19) pandemic caused by the severe acute respiratory syndrome virus 2 (SARS-CoV-2) has impacted negatively on public health and socioeconomic status, globally. Although, there are currently no specific drugs approved, several existing drugs are being repurposed, but their successful outcomes are not guaranteed. Therefore, the search for novel therapeutics remains a priority. We screened for inhibitors of the SARS-CoV-2 main protease and the receptor-binding domain of the spike protein from an integrated library of African natural products, compounds generated from machine learning studies and antiviral drugs using AutoDock Vina. The binding mechanisms between the compounds and the proteins were characterized using LigPlot+ and molecular dynamics simulations techniques. The biological activities of the hit compounds were also predicted using a Bayesian-based approach. Six potential bioactive molecules NANPDB2245, NANPDB2403, fusidic acid, ZINC000095486008, ZINC0000556656943 and ZINC001645993538 were identified, all of which had plausible binding mechanisms with both viral receptors. Molecular dynamics simulations, including molecular mechanics Poisson-Boltzmann surface area (MM/PBSA) computations revealed stable protein-ligand complexes with all the compounds having acceptable free binding energies <-15 kJ/mol with each receptor. NANPDB2245, NANPDB2403 and ZINC000095486008 were predicted as antivirals; ZINC000095486008 as a membrane permeability inhibitor; NANPDB2403 as a cell adhesion inhibitor and RNA-directed RNA polymerase inhibitor; and NANPDB2245 as a membrane integrity antagonist. Therefore, they have the potential to inhibit viral entry and replication. These drug-like molecules were predicted to possess attractive pharmacological profiles with negligible toxicity. Novel critical residues identified for both targets could aid in a better understanding of the binding mechanisms and design of fragment-based de novo inhibitors. The compounds are proposed as worthy of further in vitro assaying and as scaffolds for the development of novel SARS-CoV-2 therapeutic molecules.

Druggable targets of SARS-CoV-2 and treatment opportunities for COVID-19

COVID-19 caused by the novel SARS-CoV-2 has been declared a pandemic by the WHO is causing havoc across the entire world. As of May end, about 6 million people have been affected, and 367 166 have died from COVID-19. Recent studies suggest that the SARS-CoV-2 genome shares about 80% similarity with the SARS-CoV-1 while their protein RNA dependent RNA polymerase (RdRp) shares 96% sequence similarity. Remdesivir, an RdRp inhibitor, exhibited potent activity against SARS-CoV-2 in vitro. 3-Chymotrypsin like protease (also known as Mpro) and papain-like protease, have emerged as the potential therapeutic targets for drug discovery against coronaviruses owing to their crucial role in viral entry and host-cell invasion. Crystal structures of therapeutically important SARS-CoV-2 target proteins, namely, RdRp, Mpro, endoribonuclease Nsp15/NendoU and receptor binding domain of CoV-2 spike protein has been resolved, which have facilitated the structure-based design and discovery of new inhibitors. Furthermore, studies have indicated that the spike proteins of SARS-CoV-2 use the Angiotensin Converting Enzyme-2 (ACE-2) receptor for its attachment similar to SARS-CoV-1, which is followed by priming of spike protein by Transmembrane protease serine 2 (TMPRSS2) which can be targeted by a proven inhibitor of TMPRSS2, camostat. The current treatment strategy includes repurposing of existing drugs that were found to be effective against other RNA viruses like SARS, MERS, and Ebola. This review presents a critical analysis of druggable targets of SARS CoV-2, new drug discovery, development, and treatment opportunities for COVID-19.

Dual inhibitors of SARS-CoV-2 proteases: pharmacophore and molecular dynamics based drug repositioning and phytochemical leads

SARS-related coronaviruses poses continual threat to humanity by rapidly mutating and emerging as severe pandemic outbreaks, including the current nCoV-19 pandemic. Hence a rapid drug repositioning and lead identification strategy are required to mitigate these outbreaks. We report a pharmacophore and molecular dynamics-based approach for drug repositioning and lead identification against dual targets (3CLp and PLp) of SARS-CoV-2. The pharmacophore model of 3CLp inhibitors was apolar with two aromatic and two H-bond acceptors, whereas that of PLp was relatively polar, bearing one aromatic and three H-bond acceptors. Pharmacophore-based virtual screening yielded six existing FDA-approved drugs and twelve natural products with both the pharmacophoric features. Among them are nelfinavir, tipranavir and licochalcone-D, which has shown better binding characteristics with both the proteases compared to lopinavir. The molecular dynamics revealed that the connecting loop (residues 176–199) of 3CLp is highly flexible, and hence, inhibitors should avoid high-affinity interactions with it. Lopinavir, due to its high affinity with the loop region, exhibited unstable binding. Further, the van der Waals size of the 3CLp inhibitors positively correlated with their binding affinity with 3CLp. However, the van der Waals size of a ligand should not cross a threshold of 572ų, beyond which the ligands are likely to make high-affinity interaction with the loop and suffer unstable binding as observed in the case of lopinavir. Similarly, the total polar surface area of the ligands were found to be negatively correlated with their binding affinity with PLp.

Computational discovery of small drug-like compounds as potential inhibitors of SARS-CoV-2 main protease

A computational approach to in silico drug discovery was carried out to identify small drug-like compounds able to show structural and functional mimicry of the high affinity ligand X77, potent non-covalent inhibitor of SARS-COV-2 main protease (M^Pro). In doing so, the X77-mimetic candidates were predicted based on the crystal X77-M^Pro structure by a public web-oriented virtual screening platform Pharmit. Models of these candidates bound to SARS-COV-2 M^Pro were generated by molecular docking, quantum chemical calculations and molecular dynamics simulations. At the final point, analysis of the interaction modes of the identified compounds with M^Pro and prediction of their binding affinity were carried out. Calculation revealed 5 top-ranking compounds that exhibited a high affinity to the active site of SARS-CoV-2 M^Pro. Insights into the ligand - M^Pro models indicate that all identified compounds may effectively block the binding pocket of SARS-CoV-2 M^Pro, in line with the low values ​​of binding free energy and dissociation constant. Mechanism of binding of these compounds to M^Pro is mainly provided by van der Waals interactions with the functionally important residues of the enzyme, such as His-41, Met-49, Cys-145, Met-165, and Gln-189 that play a role of the binding hot spots assisting the predicted molecules to effectively interact with the M^Pro active site. The data obtained show that the identified X77-mimetic candidates may serve as good scaffolds for the design of novel antiviral agents able to target the active site of SARS-CoV-2 M^Pro.Communicated by Ramaswamy H. Sarma.

Ethnomedicines of Indian origin for combating COVID-19 infection by hampering the viral replication: using structure-based drug discovery approach

In the present study, we have explored the interaction of the active components from 10 different medicinal plants of Indian origin that are commonly used for treating cold and respiratory-related disorders, through molecular docking analysis. In the current scenario, COVID-19 patients experience severe respiratory syndromes, hence it is envisaged from our study that these traditional medicines are very likely to provide a favourable effect on COVID-19 infections. The active ingredients identified from these natural products are previously reported for antiviral activities against large group of viruses. Totally 47 bioactives identified from the medicinal plants were investigated against the structural targets of SARS-CoV-2 (Mpro and spike protein) and human ACE2 receptor. The top leads were identified based on interaction energies, number of hydrogen bond and other parameters that explain their potency to inhibit SARS-CoV-2. The bioactive ligands such as Cucurbitacin E, Orientin, Bis-andrographolide, Cucurbitacin B, Isocucurbitacin B, Vitexin, Berberine, Bryonolic acid, Piperine and Magnoflorine targeted the hotspot residues of SARS-CoV-2 main protease. In fact, this protease enzyme has an essential role in mediating the viral replication and therefore compounds targeting this key enzyme are expected to block the viral replication and transcription. The top scoring conformations identified through docking analysis were further demonstrated with molecular dynamics simulation. Besides, the stability of the conformation was studied in detail by investigating the binding free energy using MM-PBSA method. Overall, the study emphasized that the proposed hit Cucurbitacin E and orientin could serve as a promising scaffold for developing anti-COVID-19 drug.

Communicated by Ramaswamy H. Sarma