Benchmarking the ability of novel compounds to inhibit SARS-CoV-2 main protease using steered molecular dynamics simulations

Dual Anti-Inflammatory and Anticancer Activity of Novel 1,5-Diaryl Pyrazole Derivatives: Molecular Modeling, Synthesis, In Vitro Activity, and Dynamics Study

A series of novel 1,5-diaryl pyrazole derivatives targeting the COX enzyme were designed by combined ligand and structure-based approach. The designed molecules were then further subjected to ADMET and molecular docking studies. Out of 34 designed compounds, the top-10 molecules from the computation studies were synthesized, characterized, and evaluated for COX-2 inhibition and anti-cancer activity. Initially, the target compounds were screened for the protein denaturation assay. The results of the top-five molecules T2, T3, T5, T6, and T9 were further subjected to in vitro COX-2 enzymatic assay and anti-cancer activity. As far as COX-2 inhibitory activity is considered, two compounds, T3 and T5, exhibited the half maximum inhibitory concentration (IC50) at 0.781 µM and 0.781 µM respectively. Further, the two compounds T3 and T5, when evaluated for COX-1 inhibition, exhibited excellent inhibitory activity with T3 IC50 of 4.655μM and T5 with IC50 of 5.596 μM. The compound T5 showed more significant human COX-2 inhibition, with a selectivity index of 7.16, when compared with T3, which had a selectivity index of 5.96. Further, in vitro anti-cancer activity was screened against two cancer cell lines in which compounds T2 and T3 were active against A549 cell lines and T6 was active against the HepG2 cell line. Stronger binding energy was found by comparing MM-PBSA simulations with molecular docking, which suggests that compounds T3 and T5 have a better possibility of being effective compounds, in which T5 showed higher binding affinity. The results suggest that these compounds have the potential to develop effective COX-2 inhibitors as anti-cancer agents.

Integration of ligand and structure-based pharmacophore screening for the identification of novel natural leads against Euchromatic histone lysine methyltransferase 2 (EHMT2/G9a)

Herein, we report a blended ligand and structure-based pharmacophore screening approach to identify new natural leads against the Protein Lysine Methyltransferase 2 (EHMT2/G9a). The EHMT2/G9a has been associated with Cancer, Alzheimer's, and aging and is considered an emerging drug target having no clinically passed inhibitor. Purposefully, we developed the ligand-based pharmacophore (Pharmacophore-L) based on the common features of known inhibitors and the structure-based pharmacophore (Pharmacophore-S) based on the interaction profile of available crystal structures. The Pharmacophore-L and Pharmacophore-S were subjected to multiple tiers of validations and utilized in combination for the screening of total 741543 compounds coming from multiple databases. Additional layers of stringency were applied in the screening process to test drug-likeness (using Lipinski's rule, Veber's rule, SMARTS and ADMET filtration), to rule out any toxicity (TOPKAT analysis). The interaction profiles, stabilities, and comparative analysis against the reference were carried out by flexible docking, MD simulation, and MM-GBSA analysis, which finally led to three leads as potential inhibitors of G9a.Communicated by Ramaswamy H. Sarma.

Harnessing Brazilian biodiversity database: identification of flavonoids as potential inhibitors of SARS-CoV-2 main protease using computational approaches and all-atom molecular dynamics simulation

SARS-CoV-2 (Severe Acute Respiratory Syndrome Coronavirus 2) is the etiological agent responsible for the global outbreak of COVID-19 (Coronavirus Disease 2019). The main protease of SARS-CoV-2, Mpro, is a key enzyme that plays a vital role in mediating viral replication and transcription. In this study, a comprehensive computational approach was employed to investigate the binding affinity, selectivity, and stability of natural product candidates as potential new antivirals acting on the viral polyprotein processing mediated by SARS-CoV-2 Mpro. A library of 288 flavonoids extracted from Brazilian biodiversity was screened to select potential Mpro inhibitors. An initial filter based on Lipinski's rule of five was applied, and 204 compounds that did not violate any of the Lipinski rules were selected. The compounds were then docked into the active site of Mpro using the GOLD program, and the poses were subsequently re-scored using MM-GBSA (Molecular Mechanics Generalized Born Surface Area) binding free energy calculations performed by AmberTools23. The top five flavonoids with the best MM-GBSA binding free energy values were selected for analysis of their interactions with the active site residues of the protein. Next, we conducted a toxicity and drug-likeness analysis, and non-toxic compounds were subjected to molecular dynamics simulation and free energy calculation using the MM-PBSA (Molecular Mechanics Poisson-Boltzmann Surface Area) method. It was observed that the five selected flavonoids had lower MM-GBSA binding free energy with Mpro than the co-crystal ligand. Furthermore, these compounds also formed hydrogen bonds with two important residues, Cys145 and Glu166, in the active site of Mpro. Two compounds that passed the drug-likeness filter showed stable conformations during the molecular dynamics simulations. Among these, NuBBE_867 exhibited the best MM-PBSA binding free energy value compared to the crystallographic inhibitor. Therefore, this study suggests that NuBBE_867 could be a potential inhibitor against the main protease of SARS-CoV-2 and may be further examined to confirm our results.

Identification of synthetically tractable MERS-CoV main protease inhibitors using structure-based virtual screening and molecular dynamics potential of mean force (PMF) calculations

The Middle East Respiratory Syndrome Coronavirus (MERS-CoV) is a potentially lethal infection that presents a substantial threat to health, especially in Middle East nations. Given that no FDA-approved specific therapy for MERS infection exists, designing and discovering a potent antiviral therapy for MERS-CoV is crucial. One pivotal strategy for inhibiting MERS replication is to focus on the viral main protease (Mpro). In this study, we identify potential novel Mpro inhibitors employing structure-based virtual screening of our recently reported Ugi reaction-derived library (URDL) consisting of cherry-picked molecules from the literature. The key features of the URDL library include synthetic tractability (1-2 pot synthesis) of the molecules scaffold and unexplored chemical space. The hits were ranked based on the docking score, MM-GBSA free energy of binding, and the interaction pattern with the active site residues. A molecular dynamics (MD) simulation study was performed for the first two top-ranked compounds to analyze the stability and free binding energy based on the molecular mechanics Poisson-Boltzmann surface area. The potential mean force calculated from the steered molecular dynamics (SMD) simulations of the hits indicates improved H-bond potential, enhanced conformational stability, and binding affinity toward the target, compared to the cocrystallized ligand. The discovered hits represent novel synthetically tractable scaffolds as potential MERS-CoV Mpro inhibitors.Communicated by Ramaswamy H. Sarma.

Discovery of novel HDAC8 inhibitors from natural compounds by in silico high throughput screening

A class I histone deacetylase HDAC8 is associated with several diseases, including cancer, intellectual impairment and parasite infection. Most of the HDAC inhibitors that have so far been found to inhibit HDAC8 limit their efficacy in the clinic by producing toxicities. It is therefore very desirable to develop specific HDAC8 inhibitors. The emergence of HDAC inhibitors derived from natural sources has become quite popular. In recent decades, it has been shown that naturally occurring HDAC inhibitors have strong anticancer properties. A total of 0.2 million natural compounds were screened against HDAC8 from the Universal Natural Product Database (UNPD). Molecular docking was performed for these natural compounds and the top six hits were obtained. In addition, molecular dynamics (MD) simulations were used to evaluate the structural stability and binding affinity of the inhibitors, which showed that the protein-ligand complexes remained stable throughout the 100 ns simulation. MM-PBSA method demonstrated that the selected compounds have high affinity towards HDAC8. We infer from our findings that Hit-1 (-29.35 kcal mol-1), Hit-2 (-29.15 kcal mol-1) and Hit-6 (-30.28 kcal mol-1) have better binding affinity and adhesion to ADMET (absorption, distribution, metabolism, excretion and toxicity) characteristics against HDAC8. To compare our discussions and result in an effective way. We performed molecular docking, MD and MM-PBSA analysis for the FDA-approved drug romidepsin. The above results show that our hits show better binding affinity than the compound romidepsin (-12.03 ± 4.66 kcal mol-1). The important hotspot residues Asp29, Ile34, Trp141, Phe152, Asp267, Met274 and Tyr306 have significantly contributed to the protein-ligand interaction. These findings suggest that in vitro testing and additional optimization may lead to the development of HDAC8 inhibitors.Communicated by Ramaswamy H. Sarma.

Phytochemicals of Cocculus hirsutus deciphered SARS-CoV-2 inhibition by targeting main proteases in molecular docking, simulation, and pharmacological analyses

The COVID-19 pandemic is spreading rapidly due to the outbreak of novel coronavirus SARS-CoV-2 across the globe. Anti-COVID-19 drugs are urgently required in this situation. In this regard, the discovery of promising new anti-COVID-19 moieties is expected from traditional medicine. The study is aimed to discover phytochemicals of Cocculus hirsutus having anti-COVID-19 activity via inhibiting main proteases of SARS-CoV-2. Main proteases (Mpro) of SARS-CoV-2 serve as a protuberant target for anti-COVID-19 drug discovery because it is a key enzyme of coronaviruses and has a pivotal role in mediating viral replication and transcription that makes it an attractive drug target. Recent studies indicated the utility of C. hirsutus in the treatment of viral disorders like Dengue. Phytochemicals from C. hirsutus were docked against SARS-CoV-2 main proteases (6LU7, 5R7Y, 5R7Z, 5R80, 5R81, 5R82) using the PyRx virtual screen tool and discovery studio visualizer. Further, molecular dynamics simulations were performed (for 100 ns) to see conformational stability for all complexes. Pharmacokinetic properties and drug-likeness prediction of selected C. hirsutus phytoconstituents were also performed. Betulin, coclaurine, and quinic acid of C. hirsutus were found promising with significant binding affinity to SARS-CoV-2 Mpro in comparison to control. They have shown stable interactions with the amino acid residues present on the active site of most of the SARS-CoV-2 Mpro and were found as promising anti-COVID-19 candidates. These compounds could be potential leads for the development of target-specific anti-COVID-19 therapeutics while ethnomedicinal uses of this herb could further needed for its detailed antiviral therapeutic exploration.Communicated by Ramaswamy H. Sarma.

Viral target and metabolism-based rationale for combined use of recently authorized small molecule COVID-19 medicines: Molnupiravir, nirmatrelvir, and remdesivir

The COVID-19 pandemic remains a major health concern worldwide, and SARS-CoV-2 is continuously evolving. There is an urgent need to identify new antiviral drugs and develop novel therapeutic strategies. Combined use of newly authorized COVID-19 medicines including molnupiravir, nirmatrelvir, and remdesivir has been actively pursued. Mechanistically, nirmatrelvir inhibits SARS-CoV-2 replication by targeting the viral main protease (Mpro ), a critical enzyme in the processing of the immediately translated coronavirus polyproteins for viral replication. Molnupiravir and remdesivir, on the other hand, inhibit SARS-CoV-2 replication by targeting RNA-dependent RNA-polymerase (RdRp), which is directly responsible for genome replication and production of subgenomic RNAs. Molnupiravir targets RdRp and induces severe viral RNA mutations (genome), commonly referred to as error catastrophe. Remdesivir, in contrast, targets RdRp and causes chain termination and arrests RNA synthesis of the viral genome. In addition, all three medicines undergo extensive metabolism with strong therapeutic significance. Molnupiravir is hydrolytically activated by carboxylesterase-2 (CES2), nirmatrelvir is inactivated by cytochrome P450-based oxidation (e.g., CYP3A4), and remdesivir is hydrolytically activated by CES1 but covalently inhibits CES2. Additionally, remdesivir and nirmatrelvir are oxidized by the same CYP enzymes. The distinct mechanisms of action provide strong rationale for their combined use. On the other hand, these drugs undergo extensive metabolism that determines their therapeutic potential. This review discusses how metabolism pathways and enzymes involved should be carefully considered during their combined use for therapeutic synergy.

A novel monospecific tetravalent IgG1-(scFv)2 version shown enhanced neutralizing and Fc-mediated effector functions against SARS-CoV-2

Neutralizing monoclonal antibodies (nMABs) have been proved to be effective therapeutics in treating coronavirus disease (COVID-19). To enhance the potency of nMAB 553-15, we generated a novel monospecific tetravalent IgG1-(scFv)₂ version. This was achieved by covalently fusing two forms of 553-15-derived single chain variable fragments (scFv) to the C-terminus of the hIgG1 (human Immunoglobulin G1) Fc fragment. We found that the Fc-fused VL-linker-VH format achieved similar binding affinity and neutralizing behavior as 553-15. The tetravalent versions were constructed by fusing the scFv domains to the C-terminus of nMAB 553-15. As a result, the tetravalent version 55,315-VLVH exhibited significantly higher binding activity to target spike protein variants and enhanced neutralization against VOCs (variants of concern) pseudovirus compared to 553-15. We also measured the Fc effector responses of candidates using wild-type Spike-expressing CHOK1 cells. The 55,315-VLVH enhanced the function of ADCP (antibody-dependent cellular phagocytosis) but had similar IL-6 release levels compared to the bivalent 553-15. It seemed that the novel tetravalent version avoids the pro-inflammatory effect induced by macrophage activation. However, the 55,315-VLVH displayed slightly increased potency in ADCC (antibody-dependent cell-mediated cytotoxicity) and CDC (complement-dependent cytotoxicity), which might contribute to higher systemic inflammation. Further investigation is necessary to determine whether the tetravalent version is beneficial to balance efficiency and safety against COVID-19.

Repurposing and computational design of PARP inhibitors as SARS-CoV-2 inhibitors

Coronavirus disease 2019 (COVID-19) is a recent pandemic that caused serious global emergency. To identify new and effective therapeutics, we employed a drug repurposing approach. The poly (ADP ribose) polymerase inhibitors were used for this purpose and were repurposed against the main protease (Mpro) target of severe acute respiratory syndrome Coronavirus 2 (SARS-CoV-2). The results from these studies were used to design compounds using the 'Grow Scaffold' modules available on Discovery Studio v2018. The three designed compounds, olaparib 1826 and olaparib 1885, and rucaparib 184 demonstrated better CDOCKER docking scores for Mpro than their parent compounds. Moreover, the compounds adhered to Lipinski's rule of five and demonstrated a synthetic accessibility score of 3.55, 3.63, and 4.30 for olaparib 1826, olaparib 1885, and rucaparib 184, respectively. The short-range Coulombic and Lennard-Jones potentials also support the potential binding of the modified compounds to Mpro. Therefore, we propose these three compounds as novel SARS-CoV-2 inhibitors.

In-silico study: docking simulation and molecular dynamics of peptidomimetic fullerene-based derivatives against SARS-CoV-2 Mpro

COVID-19 is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2, has become a global pandemic resulting in significant morbidity and mortality. This study presents 12 new peptidomimetic fullerene-based derivatives in three groups that are investigated theoretically as SARS-CoV-2 M^pro inhibitors to increase the chance of treating COVID-19. Studied compounds are designed and optimized at B88-LYP/DZVP method. Molecular descriptors results show the stability and reactivity of the compounds with M^pro, especially in the 3rd group (Ser compounds). However, Lipinski's Rule of Five values indicates that the compounds are not suitable as oral drugs. Furthermore, molecular docking simulations are carried out to investigate the binding affinity and interaction modes of the top five compounds (compounds 1, 9, 11, 2, and 10) with the M^pro protein, which have the lowest binding energy. Molecular dynamics simulations are also performed to evaluate the stability of the protein-ligand complexes with compounds 1 and 9 and compare them with natural substrate interaction. The analysis of RMSD, H-bonds, Rg, and SASA indicates that both compounds 1 (Gly-α acid) and 9 (Ser-α acid) have good stability and strong binding affinity with the M^pro protein. However, compound 9 shows slightly better stability and binding affinity compared to compound 1.

Discovery of Chalcone-Based Hybrid Structures as High Affinity and Site-Specific Inhibitors against SARS-CoV-2: A Comprehensive Structural Analysis Based on Various Host-Based and Viral Targets

Previous studies indicated that natural-based chalcones have significant inhibitory effects on the coronavirus enzymes 3CLpro and PLpro as well as modulation of some host-based antiviral targets (HBATs). In this study, a comprehensive computational and structural study was performed to investigate the affinity of our compound library consisting of 757 chalcone-based structures (CHA-1 to CHA-757) for inhibiting the 3CLpro and PLpro enzymes and against twelve selected host-based targets. Our results indicated that CHA-12 (VUF 4819) is the most potent and multi-target inhibitor in our chemical library over all viral and host-based targets. Correspondingly, CHA-384 and its congeners containing ureide moieties were found to be potent and selective 3CLpro inhibitors, and benzotriazole moiety in CHA-37 was found to be a main fragment for inhibiting the 3CLpro and PLpro. Surprisingly, our results indicate that the ureide and sulfonamide moieties are integral fragments for the optimum 3CLpro inhibition while occupying the S1 and S3 subsites, which is fully consistent with recent reports on the site-specific 3CLpro inhibitors. Finding the multi-target inhibitor CHA-12, previously reported as an LTD4 antagonist for the treatment of inflammatory pulmonary diseases, prompted us to suggest it as a concomitant agent for relieving respiratory symptoms and suppressing COVID-19 infection.

Deciphering molecular mechanisms of SARS-CoV-2 pathogenesis and drug repurposing through GRN motifs: a comprehensive systems biology study

The world has recently been plagued by a new coronavirus infection called SARS-CoV-2. This virus may lead to severe acute respiratory syndrome followed by multiple organ failure. SARS-CoV-2 has approximately 80-90% genetic similarity to SARS-CoV. Given the limited omics data available for host response to the viruses (more limited data for SARS-CoV-2), we attempted to unveil the crucial molecular mechanisms underlying the SARS-CoV-2 pathogenesis by comparing its regulatory network motifs with SARS-CoV. We also attempted to identify the non-shared crucial molecules and their functions to predict the specific mechanisms for each infection and the processes responsible for their different manifestations. Deciphering the crucial shared and non-shared mechanisms at the molecular level and signaling pathways underlying both diseases may help shed light on their pathogenesis and pave the way for other new drug repurposing against COVID-19. We constructed the GRNs for host response to SARS-CoV and SARS-CoV-2 pathogens (in vitro) and identified the significant 3-node regulatory motifs by analyzing them topologically and functionally. We attempted to identify the shared and non-shared regulatory elements and signaling pathways between their host responses. Interestingly, our findings indicated that NFKB1, JUN, STAT1, FOS, KLF4, and EGR1 were the critical shared TFs between motif-related subnetworks in both SARS and COVID-1, which are considered genes with specific functions in the immune response. Enrichment analysis revealed that the NOD-like receptor signaling, TNF signaling, and influenza A pathway were among the first significant pathways shared between SARS and COVID-19 up-regulated DEGs networks, and the term "metabolic pathways" (hsa01100) among the down-regulated DEGs networks. WEE1, PMAIP1, and TSC22D2 were identified as the top three hubs specific to SARS. However, MYPN, SPRY4, and APOL6 were the tops specific to COVID-19 in vitro. The term "Complement and coagulation cascades" pathway was identified as the first top non-shared pathway for COVID-19 and the MAPK signaling pathway for SARS. We used the identified crucial DEGs to construct a drug-gene interaction network to propose some drug candidates. Zinc chloride, Fostamatinib, Copper, Tirofiban, Tretinoin, and Levocarnitine were the six drugs with higher scores in our drug-gene network analysis.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03518-x.

Investigation of the effects of N-Acetylglucosamine on the stability of the spike protein in SARS-CoV-2 by molecular dynamics simulations

A lot of effort has been made in developing vaccine and therapeutic agents against the SARS-CoV-2, concentrating on the Spike protein that binds angiotensin-converting enzyme 2 on human cells. Nowadays, some researches study the role of the N-linked glycans as potential targets for vaccines and new agents. Due to the flexibility and diversity of the N-linked glycans, in this work, we focus on the N-Acetylglucosamine moiety, which is the precursor of nearly all eukaryotic glycans. We performed molecular dynamics simulations to study the effects of the N-Acetylglucosamine on the stability of the spike glycoprotein in SARS-CoV-2. After a 100 ns of simulation on the spike proteins without and with the N-Acetylglucosamine molecules, we found that the presence of N-Acetylglucosamine increases the local stability in their vicinity; even though their effect on the full structure is negligible. Thus; it can be inferred that the N-Acetylglucosamine moieties can potentially affect the interaction of the S protein with the ACE2 receptor. We also found that the S1 domain is more flexible than the S2 domain. We propose which of the experimentally observed glycans found on the spike may be more functional than the others. Detailed understanding of glycans is key for the development of new therapeutic strategies.

Triterpene Derivatives as Potential Inhibitors of the RBD Spike Protein from SARS-CoV-2: An In Silico Approach

The appearance of a new coronavirus, SARS-CoV-2, in 2019 kicked off an international public health emergency. Although rapid progress in vaccination has reduced the number of deaths, the development of alternative treatments to overcome the disease is still necessary. It is known that the infection begins with the interaction of the spike glycoprotein (at the virus surface) and the angiotensin-converting enzyme 2 cell receptor (ACE2). Therefore, a straightforward solution for promoting virus inhibition seems to be the search for molecules capable of abolishing such attachment. In this work, we tested 18 triterpene derivatives as potential inhibitors of SARS-CoV-2 against the receptor-binding domain (RBD) of the spike protein by means of molecular docking and molecular dynamics simulations, modeling the RBD S1 subunit from the X-ray structure of the RBD-ACE2 complex (PDB ID: 6M0J). Molecular docking revealed that at least three triterpene derivatives of each type (i.e., oleanolic, moronic and ursolic) present similar interaction energies as the reference molecule, i.e., glycyrrhizic acid. Molecular dynamics suggest that two compounds from oleanolic and ursolic acid, OA5 and UA2, can induce conformational changes capable of disrupting the RBD-ACE2 interaction. Finally, physicochemical and pharmacokinetic properties simulations revealed favorable biological activity as antivirals.

Targeting viral proteins for restraining SARS-CoV-2: focusing lens on viral proteins beyond spike for discovering new drug targets

INTRODUCTION

Emergence of highly infectious SARS-CoV-2 variants are reducing protection provided by current vaccines, requiring constant updates in antiviral approaches. The virus encodes four structural and sixteen nonstructural proteins which play important roles in viral genome replication and transcription, virion assembly, release , entry into cells, and compromising host cellular defenses. As alien proteins to host cells, many viral proteins represent potential targets for combating the SARS-CoV-2.

AREAS COVERED

Based on literature from PubMed and Web of Science databases, the authors summarize the typical characteristics of SARS-CoV-2 from the whole viral particle to the individual viral proteins and their corresponding functions in virus life cycle. The authors also discuss the potential and emerging targeted interventions to curb virus replication and spread in detail to provide unique insights into SARS-CoV-2 infection and countermeasures against it.

EXPERT OPINION

Our comprehensive analysis highlights the rationale to focus on non-spike viral proteins that are less mutated but have important functions. Examples of this include: structural proteins (e.g. nucleocapsid protein, envelope protein) and extensively-concerned nonstructural proteins (e.g. NSP3, NSP5, NSP12) along with the ones with relatively less attention (e.g. NSP1, NSP10, NSP14 and NSP16), for developing novel drugs to overcome resistance of SARS-CoV-2 variants to preexisting vaccines and antibody-based treatments.

Ginkgolic acids inhibit SARS-CoV-2 and its variants by blocking the spike protein/ACE2 interplay

Targeting the interaction between the spike protein receptor binding domain (S-RBD) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and angiotensin-converting enzyme 2 (ACE2) is a potential therapeutic strategy for treating coronavirus disease 2019 (COVID-19). However, we still lack small-molecule drug candidates for this target due to the missing knowledge in the hot spots for the protein-protein interaction. Here, we used NanoBiT technology to identify three Ginkgolic acids from an in-house traditional Chinese medicine (TCM) library, and they interfere with the S-RBD/ACE2 interplay. Our pseudovirus assay showed that one of the compounds, Ginkgolic acid C17:1 (GA171), significantly inhibits the entry of original SARS-CoV-2 and its variants into the ACE2-overexpressed HEK293T cells. We investigated and proposed the binding sites of GA171 on S-RBD by combining molecular docking and molecular dynamics simulations. Site-directed mutagenesis and surface plasmon resonance revealed that GA171 specifically binds to the pocket near R403 and Y505, critical residues of S-RBD for S-RBD interacting with ACE2. Thus, we provide structural insights into developing new small-molecule inhibitors and vaccines against the proposed S-RBD binding site.

Characterizing the ligand-binding affinity toward SARS-CoV-2 Mpro via physics- and knowledge-based approaches

Computational approaches, including physics- and knowledge-based methods, have commonly been used to determine the ligand-binding affinity toward SARS-CoV-2 main protease (Mpro or 3CLpro). Strong binding ligands can thus be suggested as potential inhibitors for blocking the biological activity of the protease. In this context, this paper aims to provide a short review of computational approaches that have recently been applied in the search for inhibitor candidates of Mpro. In particular, molecular docking and molecular dynamics (MD) simulations are usually combined to predict the binding affinity of thousands of compounds. Quantitative structure-activity relationship (QSAR) is the least computationally demanding and therefore can be used for large chemical collections of ligands. However, its accuracy may not be high. Moreover, the quantum mechanics/molecular mechanics (QM/MM) method is most commonly used for covalently binding inhibitors, which also play an important role in inhibiting the activity of SARS-CoV-2. Furthermore, machine learning (ML) models can significantly increase the searching space of ligands with high accuracy for binding affinity prediction. Physical insights into the binding process can then be confirmed via physics-based calculations. Integration of ML models into computational chemistry provides many more benefits and can lead to new therapies sooner.

Evaluation of the anti-diabetic drug sitagliptin as a novel attenuate to SARS-CoV-2 evidence-based in silico: molecular docking and molecular dynamics

The current outbreak of COVID-19 cases worldwide has been responsible for a significant number of deaths, especially in hospitalized patients suffering from comorbidities, such as obesity, diabetes, hypertension. The disease not only has prompted an interest in the pathophysiology, but also it has propelled a massive race to find new anti-SARS-CoV-2 drugs. In this scenario, known drugs commonly used to treat other diseases have been suggested as alternative or complementary therapeutics. Herein we propose the use of sitagliptin, an inhibitor of dipeptidyl peptidase-4 (DPP₄) used to treat type-II diabetes, as an agent to block and inhibit the activity of two proteases, 3CL^pro and PL^pro, related to the processing of SARS-CoV-2 structural proteins. Inhibition of these proteases may possibly reduce the viral load and infection on the host by hampering the synthesis of new viruses, thus promoting a better outcome. In silico assays consisting in the modeling of the ligand sitagliptin and evaluation of its capacity to interact with 3CL^pro and PL^pro through the prediction of the ligand bioactivity, molecular docking, overlapping of crystal structures, and molecular dynamic simulations were conducted. The experiments indicate that sitagliptin can interact and bind to both targets. However, this interaction seems to be stronger and more stable to 3CL^pro (ΔG = −7.8 kcal mol⁻¹), when compared to PL^pro (ΔG = −7.5 kcal mol⁻¹). This study suggests that sitagliptin may be suitable to treat COVID-19 patients, beyond its common use as an anti-diabetic medication. In vivo studies may further support this hypothesis.

Inhibition of nonstructural protein 15 of SARS-CoV-2 by golden spice: A computational insight

The quick widespread of the coronavirus and speedy upsurge in the tally of cases demand the fast development of effective drugs. The uridine-directed endoribonuclease activity of nonstructural protein 15 (Nsp15) of the coronavirus is responsible for the invasion of the host immune system. Therefore, developing potential inhibitors against Nsp15 is a promising strategy. In this concern, the in silico approach can play a significant role, as it is fast and cost-effective in comparison to the trial and error approaches of experimental investigations. In this study, six turmeric derivatives (curcuminoids) were chosen for in silico analysis. The molecular interactions, pharmacokinetics, and drug-likeness of all the curcuminoids were measured. Further, the stability of Nsp15-curcuminoids complexes was appraised by employing molecular dynamics (MD) simulations and MM-PBSA approaches. All the molecules were affirmed to have strong interactions and pharmacokinetic profile. The MD simulations data stated that the Nsp15-curcuminoids complexes were stable during simulations. All the curcuminoids showed stable and high binding affinity, and these curcuminoids could be admitted as potential modulators for Nsp15 inhibition.

Pharmacogenomics deliberations of 2-deoxy-d-glucose in the treatment of COVID-19 disease: an in silico approach

The outbreak of COVID-19 caused by the coronavirus (SARS-CoV-2) prompted number of computational and laboratory efforts to discover molecules against the virus entry or replication. Simultaneously, due to the availability of clinical information, drug-repurposing efforts led to the discovery of 2-deoxy-d-glucose (2-DG) for treating COVID-19 infection. 2-DG critically accumulates in the infected cells to prevent energy production and viral replication. As there is no clarity on the impact of genetic variations on the efficacy and adverse effects of 2-DG in treating COVID-19 using in silico approaches, we attempted to extract the genes associated with the 2-DG pathway using the Comparative Toxicogenomics Database. The interaction between selected genes was assessed using ClueGO, to identify the susceptible gene loci for SARS-CoV infections. Further, SNPs that were residing in the distinct genomic regions were retrieved from the Ensembl genome browser and characterized. A total of 80 SNPs were retrieved using diverse bioinformatics resources after assessing their (a) detrimental influence on the protein stability using Swiss-model, (b) miRNA regulation employing miRNASNP3, PolymiRTS, MirSNP databases, (c) binding of transcription factors by SNP2TFBS, SNPInspector, and (d) enhancers regulation using EnhancerDB and HaploReg reported A2M rs201769751, PARP1 rs193238922 destabilizes protein, six polymorphisms of XIAP effecting microRNA binding sites, EGFR rs712829 generates 15 TFBS, BECN1 rs60221525, CASP9 rs4645980, SLC2A2 rs5393 impairs 14 TFBS, STK11 rs3795063 altered 19 regulatory motifs. These data may provide the relationship between genetic variations and drug effects of 2-DG which may further assist in assigning the right individuals to benefit from the treatment.

Targeting the SARS-CoV-2 HR1 with Small Molecules as Inhibitors of the Fusion Process

The rapid and global propagation of the novel human coronavirus that causes severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has produced an immediate urgency to discover promising targets for the treatment of this virus. In this paper, we studied the spike protein S2 domain of SARS-CoV-2 as it is the most conserved component and controls the crucial fusion process of SARS-CoV-2 as a target for different databases of small organic compounds. Our in silico methodology, based on pharmacophore modeling, docking simulation and molecular dynamics simulations, was first validated with ADS-J1, a potent small-molecule HIV fusion inhibitor that has already proved effective in binding the HR1 domain and inhibiting the fusion core of SARS-CoV-1. It then focused on finding novel small molecules and new peptides as fusion inhibitors. Our methodology identified several small molecules and peptides as potential inhibitors of the fusion process. Among these, NF 023 hydrate (MolPort-006-822-583) is one of the best-scored compounds. Other compounds of interest are ZINC00097961973, Salvianolic acid, Thalassiolin A and marine_160925_88_2. Two interesting active peptides were also identified: AP00094 (Temporin A) and AVP1227 (GBVA5). The inhibition of the spike protein of SARS-CoV-2 is a valid target to inhibit the virus entry in human cells. The discussed compounds reported in this paper led to encouraging results for future in vitro tests against SARS-CoV-2.

Computational design of a broad-spectrum multi-epitope vaccine candidate against seven strains of human coronaviruses

Spike (S) proteins are an attractive target as it mediates the binding of the SARS-CoV-2 to the host through ACE-2 receptors. We hypothesize that the screening of the S protein sequences of all the seven known HCoVs would result in the identification of potential multi-epitope vaccine candidates capable of conferring immunity against various HCoVs. In the present study, several machine learning-based in-silico tools were employed to design a broad-spectrum multi-epitope vaccine candidate targeting the S protein of seven known strains of human coronaviruses. Herein, multiple B-cell epitopes and T-cell epitopes (CTL and HTL) were predicted from the S protein sequences of all seven known HCoVs. Post-prediction they were linked together with an adjuvant to construct a potential broad-spectrum vaccine candidate. Secondary and tertiary structures were predicted and validated, and the refined 3D-model was docked with an immune receptor. The vaccine candidate was evaluated for antigenicity, allergenicity, solubility, and its ability to achieve high-level expression in bacterial hosts. Finally, the immune simulation was carried out to evaluate the immune response after three vaccine doses. The designed vaccine is antigenic (with or without the adjuvant), non-allergenic, binds well with TLR-3 receptor and might elicit a diverse and strong immune response.

Identification of antiviral peptide inhibitors for receptor binding domain of SARS-CoV-2 omicron and its sub-variants: an in-silico approach

Omicron, a variant of concern (VOC) of SARS-CoV-2, emerged in South Africa in November 2021. Omicron has been continuously acquiring a series of new mutations, especially in the spike (S) protein that led to high infectivity and transmissibility. Peptides targeting the receptor-binding domain (RBD) of the spike protein by which omicron and its variants attach to the host receptor, angiotensin-converting enzyme (ACE2) can block the viral infection at the first step. This study aims to identify antiviral peptides from the Antiviral peptide database (AVPdb) and HIV-inhibitory peptide database (HIPdb) against the RBD of omicron by using a molecular docking approach. The lead RBD binder peptides obtained through molecular docking were screened for allergenicity and physicochemical criteria (isoelectric point (pI) and net charge) required for peptide-based drugs. The binding affinity of the best five peptide inhibitors with the RBD of omicron was validated further by molecular dynamics (MD) simulation. Our result introduces five antiviral peptides, including AVP1056, AVP1059, AVP1225, AVP1801, and HIP755, that may effectively hinder omicron-host interactions. It is worth mentioning that all the three major sub-variants of omicron, BA.1 (B.1.1.529.1), BA.2 (B.1.1.529.2), and BA.3 (B.1.1.529.3), exhibits conserved ACE-2 interacting residues. Hence, the screened antiviral peptides with similar affinity can also interrupt the RBD-mediated invasion of different major sub-variants of omicron. Altogether, these peptides can be considered in the peptide-based therapeutics development for omicron treatment after further experimentation.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03258-4.

Pharmacological Profile of Nigella sativa Seeds in Combating COVID-19 through In-Vitro and Molecular Docking Studies

COVID-19 infection is associated with elevated oxidative stress, systemic hyper-inflammatory responses, endothelial dysfunction, and red blood cell membrane deformability. Nigella sativa extract is widely used in alternative and complementary medicine systems in a large population, due to its highly therapeutic, economic, natural, and safe nature. The aim of this study was to evaluate the effect of N. sativa extract on oxidative stress, hemolysis, proteolysis, and glycation through in vitro studies, as well as to find out its anti-viral potential against severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) using in silico studies. N. sativa seed extract (at 600 µg/mL) displayed 67.33% scavenging activity in the 2,2-diphenyl-1-picrylhydrazyl (DPPH) test, and 70.28% hydrogen peroxide reducing activity. N. sativa exhibited anti-proteolytic activity by decreasing heat-induced denaturation of bovine serum albumin (BSA) and egg albumin by 63.14% and 57.95%, respectively, and exhibited anti-proteinase potential of 66.28% at 600 μg/mL. In addition, heat-induced hemolysis and hypersalinity-induced hemolysis were inhibited by 57.86% and 61.7%, respectively, by the N. sativa seeds. N. sativa also inhibited browning intensity by 56.38%, and percent aggregation index by 51.38%, amyloid structure by 48.28%, and AGE-specific fluorescence by 52.18%, thereby protecting the native structure of BSA from glycation. The binding interactions between bioactive molecules of N. sativa seed with SARS-CoV-2 spike glycoprotein were proven by using in silico molecular docking tools. The functional amino acids involved in the interactions are Asp467, Thr108, Thr114, Ile468, Asn234, Gln155, Glu465, Arg466, Gly232, and Ile233, indicating the inhibiting property of N. sativa on SARS-CoV-2. Finally, we may infer that phytoconstituents of N. sativa seeds have the potential to protect against the spike protein of SARS-CoV-2. Studies on N. sativa seeds might act as a path to develop a potent alternative therapy against viral infections, especially COVID-19 infection, in the future. However, the limitations linked with the use of natural products are also needed to be considered in this regard.