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

A top-down approach for studying the in-silico effect of the novel phytocompound tribulusamide B on the inhibition of Nipah virus transmission through targeting fusion glycoprotein and matrix protein

The proteins of Nipah virus ascribe to its lifecycle and are crucial to infections caused by the virus. In the absence of approved therapeutics, these proteins can be considered as drug targets. This study examined the potential of fifty-three (53) natural compounds to inhibit Nipah virus fusion glycoprotein (NiV F) and matrix protein (NiV M) in silico. The molecular docking experiment, supported by the principal component analysis (PCA), showed that out of all the phytochemicals considered, Tribulusamide B had the highest inhibitory potential against the target proteins NiV F and NiV M (-9.21 and -8.66 kcal mol-1, respectively), when compared to the control drug, Ribavirin (-7.01 and -6.52 kcal mol-1, respectively). Furthermore, it was found that Tribulusamide B pharmacophores, namely, hydrogen donors, acceptors, aromatic and hydrophobic groups, contributed towards the effective residual interactions with the target proteins. The molecular dynamic simulation further validated the results of the docking studies and concluded that Tribulusamide B formed a stable complex with the target proteins. The data obtained from MM-PBSA study further explained that the phytochemical could strongly bind with NiV F (-31.26 kJ mol-1) and NiV M (-40.26 kJ mol-1) proteins in comparison with the control drug Ribavirin (-13.12 and -13.94 kJ mol-1, respectively). Finally, the results indicated that Tribulusamide B, a common inhibitor effective against multiple proteins, can be considered a potential therapeutic entity in treating the Nipah virus infection.

In silico screening of potential inhibitors from Cordyceps species against SARS-CoV-2 main protease

Coronavirus disease 2019 (COVID-19) is a result of a retroviral infection of SARS-CoV-2. Due to its virulence and high infection rate, it is a matter of serious concern and a global health emergency. Currently available COVID-19 vaccines approved by regulatory bodies around the world have been shown to provide significant protection against COVID-19. But no vaccine is 100% effective at preventing infection, also they have varying efficacy rates and different side effects. However, the main protease (Mpro) of SARS-CoV-2 has been identified as a key drug target due to its essential role in viral infection and its minimal similarity with human proteases. Cordyceps mushrooms have been found to have various therapeutic properties that could effectively combat SARS-CoV-2, including improve lung functioning, anti-viral, immunomodulators, anti-infectious, and anti-inflammatory. The present study aims to screen and evaluate the inhibitory potential of the bioactive molecules from the Cordyceps species against the Mpro of SARS-CoV-2. The bioactive molecules were screened based on their docking score, molecular interactions in the binding pocket, ADME properties, toxicity, carcinogenicity, and mutagenicity. Among all the molecules that were tested, cordycepic acid was the most effective and promising candidate, with a binding affinity of -8.10 kcal/mol against Mpro. The molecular dynamics (MD) simulation and free binding energy calculations revealed that the cordycepic acid-Mpro complex was highly stable and showed fewer conformational fluctuations. These findings need to be investigated further through in-vitro and in-vivo studies for additional validation.Communicated by Ramaswamy H. Sarma.

Abrogation of ORF8-IRF3 binding interface with Carbon nanotube derivatives to rescue the host immune system against SARS-CoV-2 by using molecular screening and simulation approaches

The COVID-19 pandemic, caused by the SARS-CoV-2 virus, has led to over six million deaths worldwide. In human immune system, the type 1 interferon (IFN) pathway plays a crucial role in fighting viral infections. However, the ORF8 protein of the virus evade the immune system by interacting with IRF3, hindering its nuclear translocation and consequently downregulate the type I IFN signaling pathway. To block the binding of ORF8-IRF3 and inhibit viral pathogenesis a quick discovery of an inhibitor molecule is needed. Therefore, in the present study, the interface between the ORF8 and IRF3 was targeted on a high-affinity carbon nanotube by using computational tools. After analysis of 62 carbon nanotubes by multiple docking with the induced fit model, the top five compounds with high docking scores of - 7.94 kcal/mol, - 7.92 kcal/mol, - 7.28 kcal/mol, - 7.19 kcal/mol and - 7.09 kcal/mol (top hit1-5) were found to have inhibitory activity against the ORF8-IRF3 complex. Molecular dynamics analysis of the complexes revealed the high compactness of residues, stable binding, and strong hydrogen binding network among the ORF8-nanotubes complexes. Moreover, the total binding free energy for top hit1-5 was calculated to be - 43.21 ± 0.90 kcal/mol, - 41.17 ± 0.99 kcal/mol, - 48.85 ± 0.62 kcal/mol, - 43.49 ± 0.77 kcal/mol, and - 31.18 ± 0.78 kcal/mol respectively. These results strongly suggest that the identified top five nanotubes (hit1-5) possess significant potential for advancing and exploring innovative drug therapies. This underscores their suitability for subsequent in vivo and in vitro experiments, marking them as promising candidates worthy of further investigation.

Identification of putative antiviral bioactive compounds derived from family Asteraceae: An in silico approach

This computational study investigates 21 bioactive compounds from the Asteraceae family as potential inhibitors targeting the Spike protein (S protein) of SARS-CoV-2. Employing in silico methods and simulations, particularly CDOCKER and MM-GBSA, the study identifies two standout compounds, pterodontic acid and cichoric acid, demonstrating robust binding affinities (-46.1973 and -39.4265 kcal/mol) against the S protein. Comparative analysis with Favipiravir underscores their potential as promising inhibitors. Remarkably, these bioactives exhibit favorable ADMET properties, suggesting safety and efficacy. Molecular dynamics simulations validate their stability and interactions, signifying their potential as effective SARS-CoV-2 inhibitors.

Taming the cytokine storm: small molecule inhibitors targeting IL-6/IL-6α receptor

Given the increasing effectiveness of immune-based therapies, management of their associated toxicities is of utmost importance. Cytokine release syndrome (CRS), characterized by elevated levels of cytokine, poses a significant challenge following the administration of antibodies and CAR-T cell therapies. CRS also contributes to multiple organ dysfunction in severe viral infections, notably in COVID-19. Given the pivotal role of IL-6 cytokine in initiating CRS, it has been considered a most potential therapeutic target to mitigate hyperactivated immune responses. While monoclonal antibodies of IL-6 show promise in mitigating cytokine storm, concerns about immunotoxicity persist, and small molecule IL-6 antagonists remain unavailable. The present study employed sophisticated computational techniques to identify potential hit compounds as IL-6 inhibitors, with the aim of inhibiting IL-6/IL-6R protein-protein interactions. Through ligand-based pharmacophore mapping and shape similarity in combination with docking-based screening, we identified nine hit compounds with diverse chemical scaffolds as potential binders of IL-6. Further, the MD simulation of 300 ns of five virtual hits in a complex with IL-6 was employed to study the dynamic behavior. To provide a more precise prediction, binding free energy was also estimated. The identified compounds persistently interacted with the residues lining the binding site of the IL-6 protein. These compounds displayed low binding energy during MMPBSA calculations, substantiating their strong association with IL-6. This study suggests promising scaffolds as potential inhibitors of IL-6/IL-6R protein-protein interactions and provides direction for lead optimization.

In silico anti-viral assessment of phytoconstituents in a traditional (Siddha Medicine) polyherbal formulation - Targeting Mpro and pan-coronavirus post-fusion Spike protein

Background and aim

Novel nature of the viral pathogen SARS-CoV-2 and the absence of standard drugs for treatment, have been a major challenge to combat this deadly infection. Natural products offer safe and effective remedy, for which traditional ethnic medicine can provide leads. An indigenous poly-herbal formulation, Kabasura Kudineer from Siddha system of medicine was evaluated here using a combination of computational approaches, to identify potential inhibitors against two anti-SARS-CoV-2 targets - post-fusion Spike protein (structural protein) and main protease (Mpro, non-structural protein).

Experimental procedure

We docked 32 phytochemicals from the poly-herbal formulation against viral post-fusion Spike glycoprotein and Mpro followed by molecular dynamics using Schrodinger software. Drug-likeness analysis was performed using machine learning (ML) approach and pkCSM.

Results

The binding affinity of the phytochemicals in Kabasura Kudineer revealed the following top-five bioactives: Quercetin > Luteolin > Chrysoeriol > 5-Hydroxy-7,8-Dimethoxyflavone > Scutellarein against Mpro target, and Gallic acid > Piperlonguminine > Chrysoeriol > Elemol > Piperine against post-fusion Spike protein target. Quercetin and Gallic acid exhibited binding stability in complexation with their respective viral-targets and favourable free energy change as revealed by the molecular dynamics simulations and MM-PBSA analysis. In silico predicted pharmacokinetic profiling of these ligands revealed appropriate drug-likeness properties.

Conclusion

These outcomes provide: (a) potential mechanism for the anti-viral efficacy of the indigenous Siddha formulation, targeting Mpro and post-fusion Spike protein (b) top bioactive lead-molecules that may be developed as natural product-based anti-viral pharmacotherapy and their pleiotropic protective effects may be leveraged to manage co-morbidities associated with COVID-19.

In silico guided screening of active components of C. lanceolata as 3-chymotrypsin-like protease inhibitors of novel coronavirus

Despite the intense worldwide efforts towards the identification of potential anti-CoV therapeutics, no antiviral drugs have yet been discovered. Numerous vaccines are now approved for use, but they all serve as preventative measures. To effectively treat viral infections, it is crucial to find new antiviral drugs that are derived from natural sources. Various compounds with potential activity against 3 chymotrypsin-like protease (3CLpro) were reported and some are validated by bioassay studies. Therefore, we performed the computational screening of phytoconstituents of Codonopsis lanceolata to search for potential antiviral hit candidates. The curated compounds of the plant C. lanceolata were collected and downloaded from the literature. The binding affinity of the curated datasets was predicted for the target 3CLpro. Stigmasterol exhibits the highest docking score for the 3CLpro target. In addition, molecular dynamics (MD) simulations were conducted for the validation of docking results using root mean square deviation and root mean square fluctuation plots. The MD results indicated that the docked complex was stable and retained hydrogen bonding and non-bonding interactions. Furthermore, the calculation of pharmacokinetic parameters and Lipinski's rule of five suggest that C. lanceolata has the potential for drug-likeness. In order to develop new medicines for this debilitating disease, we will focus on the primary virus-based and host-based targets that can direct medicinal chemists to identify novel treatments to produce new drugs for it.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03745-2.

De novo design of bioactive phenol and chromone derivatives for inhibitors of Spike glycoprotein of SARS-CoV-2 in silico

This work presents the synthesis of 12 phenol and chromone derivatives, prepared by the analogs, and the possibility of conducting an in silico study of its derivatives as a therapeutic alternative to combat the SARS-CoV-2, pathogen responsible for COVID-19 pandemic, using its S-glycoprotein as a macromolecular target. After the initial screening for the ranking of the products, it was chosen which structure presented the best energy bond with the target. As a result, derivative 4 was submitted to a molecular growth study using artificial intelligence, where 8436 initial structures were obtained that passed through the interaction filters and similarity to the active glycoprotein pocket through the MolAICal computational package. Thus, 557 Hits with active configuration were generated, which is very promising compared to the BLA reference link for inhibiting the biological target. Molecular dynamics also simulated these compounds to verify their stability within the active protein site to seek new therapeutic propositions to fight against the pandemic. The Hit 48 and 250 are the most active compounds against SARS-CoV-2. In summary, the results show that the Hit 250 would be more active than the natural compound, which could be further developed for further testing against SARS-CoV-2. The study employs the de novo approach to design new drugs, combining artificial intelligence and molecular dynamics simulations to create efficient molecular structures. This research aims to contribute to the development of effective therapeutic strategies against the pandemic.

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.

In silico identification of natural compounds against SARS-CoV-2 main protease from Chinese herbal medicines

Aims

To determine natural compounds with inhibitory effects toward SARS-CoV-2 Mpro from Chinese herbal medicines.

Materials & methods

∼1200 natural compounds from 19 Chinese herbal medicines were collected. Computational methods including molecular docking, drug-likeness assessment, molecular dynamics simulation and molecular mechanics Poisson-Boltzmann surface area analysis were combined to obtain potent inhibitors against SARS-CoV-2 Mpro.

Results

Top 20 compounds mainly originated from Ranunculus ternatus and Picrasma quassioides exhibited low binding free energies which below -9.0 kcal/mol. Compounds Japonicone G and Picrasidine T were obtained with favorable drug-likeness. Moreover, the complex of Japonicone G and Mpro had prominent stability.

Conclusion

Natural compound Japonicone G is highly promising as a potent inhibitor against SARS-CoV-2 for further study.

SARS-CoV-2 Spike Protein Interaction Space

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a +sense single-strand RNA virus. The virus has four major surface proteins: spike (S), envelope (E), membrane (M), and nucleocapsid (N), respectively. The constitutive proteins present a high grade of symmetry. Identifying a binding site is difficult. The virion is approximately 50-200 nm in diameter. Angiotensin-converting enzyme 2 (ACE2) acts as the cell receptor for the virus. SARS-CoV-2 has an increased affinity to human ACE2 compared with the original SAR strain. Topological space, and its symmetry, is a critical component in molecular interactions. By exploring this space, a suitable ligand space can be characterized accordingly. A spike protein (S) computational model in a complex with ACE 2 was generated using silica methods. Topological spaces were probed using high computational throughput screening techniques to identify and characterize the topological space of both SARS and SARS-CoV-2 spike protein and its ligand space. In order to identify the symmetry clusters, computational analysis techniques, together with statistical analysis, were utilized. The computations are based on crystallographic protein data bank PDB-based models of constitutive proteins. Cartesian coordinates of component atoms and some cluster maps were generated and analyzed. Dihedral angles were used in order to compute a topological receptor space. This computational study uses a multimodal representation of spike protein interactions with some fragment proteins. The chemical space of the receptors (a dimensional volume) suggests the relevance of the receptor as a drug target. The spike protein S of SARS and SARS-CoV-2 is analyzed and compared. The results suggest a mirror symmetry of SARS and SARS-CoV-2 spike proteins. The results show thatSARS-CoV-2 space is variable and has a distinct topology. In conclusion, surface proteins grant virion variability and symmetry in interactions with a potential complementary target (protein, antibody, ligand). The mirror symmetry of dihedral angle clusters determines a high specificity of the receptor space.

In silico screening, ADMET analysis and MD simulations of phytochemicals of Onosma bracteata Wall. as SARS CoV-2 inhibitors

Being attracted with their cardiotonic, antidiabetic, cough relieving activity, treatment of fever, absorbent, anti-asthmatic, etc. activities reported in ancient Ayurvedic literature, phytochemicals of Onosma bracteata wall should be evaluated for their activity against SARS-CoV-2 virus. The main objective of this study is to identify a hit molecule for the inhibition of entry, replication, and protein synthesis of SARS CoV-2 virus into the host. To achieve given objective, computational virtual screening of phytochemicals of Onosma bracteata wall has been performed against three main viral targets: spike, RdRp, and M^pro. Further, the analysis of Lipinski's Ro5 and their estimation of ADMET profiles were performed using computational tools. The MD simulations studies of top hits against each viral target have also been performed for 20 ns to ensure their stability. The analysis of results revealed that Pulmonarioside C (9) and other plant compounds showed better binding affinity towards targets than existing antiviral compounds, making them probable lead compounds against SARS-CoV-2. Structural modifications and studies through in silico analysis provided the founding stone for the establishment of SARS CoV-2 inhibitory potential of phytoconstitutents of Onosma bracteata wall.

Identification of active compounds of traditional chinese medicine derived from maxing shigan decoction for COVID-19 treatment: a meta-analysis and in silico study

BACKGROUND

Coronavirus 2019 (COVID-19) poses a serious threat to human health. In China, traditional Chinese medicine (TCM), mainly based on the Maxing Shigan decoction (MXSGD), is used in conjunction with western medicine to treat COVID-19.

RESEARCH DESIGN AND METHODS

We conducted a network meta-analysis to investigate whether MXSGD-related TCM combined with western medicine is more effective in treating COVID-19 compared to western medicine alone. Additionally, using network pharmacology, cross-docking, and molecular dynamics (MD) simulation to explore the potential active compounds and possible targets underlying the therapeutic effects of MXSGD-related TCM.

RESULTS

MXSGD-related TCM combined with western medicine was better for treating COVID-19 compared to western medicine alone. Network pharmacological analysis identified 43 shared ingredients in the MXSGD-related TCM prescriptions and 599 common target genes. Cross-docking of the 43 compounds with 154 proteins that matched these genes led to the identification of 60 proteins. Pathway profiling revealed that the active ingredients participated in multiple signaling pathways that contribute to their efficacy. Molecular docking and MD simulation demonstrated that MOL007214, the most promising molecule, could stably bind to the active site of SARS-CoV-2 3CLpro.

CONCLUSION

This study demonstrates the important role of MXSGD-related TCM in the treatment of COVID-19.

Molecular modeling and simulations of some antiviral drugs, benzylisoquinoline alkaloid, and coumarin molecules to investigate the effects on Mpro main viral protease inhibition

Background

SARS-CoV-2 is a deadly viral disease and uncounted deaths occurs since its first appearance in the year 2019. The antiviral drugs, benzylisoquinoline alkaloids, and coumarin molecules were searched using different online engines for drug repurposing with SARS-CoV-2 and to investigate the effects on main viral protease (Mpro) upon their bindings.

Methods

A database composed of antiviral drugs, benzylisoquinoline alkaloids, and Coumarin molecules was screened through a molecular docking strategy to uncover the interactions of collected molecules with SARS-CoV-2 Mpro. Further, molecular dynamics simulations (MDS) were implemented for 100 ns to calculate the stability of the best complexed molecular scaffold with Mpro. The conformations of the simulated complexes were investigated by using principal component analysis (PCA) and Gibbs energy landscape (FEL) and DSSP together. Next, free binding energy (ΔG_bind) was calculated using the mmpbsa method.

Results

Molecular docking simulations demonstrate 17 molecules exhibited better binding affinity out of 99 molecules present in the database with the viral protease Mpro, followed ADMET properties and were documented. The Coumarin-EM04 molecular scaffold exhibited interactions with catalytical dyad HIS41, CYS145, and neighboring amino acids SER165 and GLN189 in the catalytical site. The crucial factor RMSD was calculated to determine the orientations of Coumarin-EM04. The Coumarin-EM04 complexed with Mpro was found stable in the binding site during MDS. Furthermore, the free energy binding ΔG_bind of Coumarin-EM04 was found to be −187.471 ± 2.230 kJ/mol, and for Remdesivir ΔG_bind was −171.926 ± 2.237 kJ/mol with SARS-CoV-2 Mpro^.

Conclusion

In this study, we identify potent molecules that exhibit interactions with catalytical dyad HIS41 and CYS145 amino acids and unravel Coumarin-EM04 exhibited ΔG_bind higher than Remdesivir against Mpro and thus may serve better antiviral agent against SARS-CoV-2.

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.

Prediction of human protein interactome of dengue virus non-structural protein 5 (NS5) and its downstream immunological implications

The non-structural protein 5 (NS5) is the most conserved protein among flaviviruses, a family that includes the dengue virus. It functions both as an RNA-dependent RNA polymerase and an RNA-methyltransferase and is therefore essential for the replication of viral RNA. The discovery that dengue virus NS5 protein (DENV-NS5) can also localize to the nucleus has resulted in renewed interest in its potential roles at the host-virus interface. In this study, we have used two complementary computational approaches in parallel - one based on linear motifs (ELM) and another based on tertiary structure of the protein (DALI) - to predict the host proteins that DENV-NS5 might interact with. Of the 42 human proteins predicted by both these methods, 34 are novel. Pathway analysis of these 42 human proteins shows that they are involved in key host cellular processes related to cell cycle regulation, proliferation, protein degradation, apoptosis, and immune responses. A focused analysis of transcription factors that directly interact with the predicted DENV-NS5 interacting proteins was performed, followed by the identification of downstream genes that are differentially expressed after dengue infection using previously published RNA-seq data. Our study provides unique insights into the DENV-NS5 interaction network and delineates mechanisms whereby DENV-NS5 could impact the host-virus interface. The novel interactors identified in this study could be potentially targeted by NS5 to modulate the host cellular environment in general, and the immune response in particular, thereby extending the role of DENV-NS5 beyond its known enzymatic functions.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03569-0.

Evaluation of therapeutic potentials of selected phytochemicals against Nipah virus, a multi-dimensional in silico study

The current study attempted to evaluate the potential of fifty-three (53) natural compounds as Nipah virus attachment glycoprotein (NiV G) inhibitors through in silico molecular docking study. Pharmacophore alignment of the four (4) selected compounds (Naringin, Mulberrofuran B, Rutin and Quercetin 3-galactoside) through Principal Component Analysis (PCA) revealed that common pharmacophores, namely four H bond acceptors, one H bond donor and two aromatic groups were responsible for the residual interaction with the target protein. Out of these four compounds, Naringin was found to have the highest inhibitory potential ( - 9.19 kcal mol^-1) against the target protein NiV G, when compared to the control drug, Ribavirin ( - 6.95 kcal mol^-1). The molecular dynamic simulation revealed that Naringin could make a stable complex with the target protein in the near-native physiological condition. Finally, MM-PBSA (Molecular Mechanics-Poisson-Boltzmann Solvent-Accessible Surface Area) analysis in agreement with our molecular docking result, showed that Naringin ( - 218.664 kJ mol^-1) could strongly bind with the target protein NiV G than the control drug Ribavirin ( - 83.812 kJ mol^-1).

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03595-y.

Design of novel pyrimidine based remdesivir analogues with dual target specificity for SARS CoV-2: A computational approach

As the world undergone unpreceded time of tragedy with the corona virus, many researchers have raised to showcase their scientific contributions in terms of novel configured anti-viral drugs until now. Herein, we designed pyrimidine based nucleotides and assessed for the binding capability with SARS-CoV-2 viral replication targets of nsp12 RNA-dependent RNA polymerase and M^pro main protease. Molecular docking studies showed all the designed compounds to possess good binding affinity, with a few compounds which outperforms the control drug remdesivir GS-5743 and its active form GS-441524. Further molecular dynamics simulation studies confirmed their stability and preservation of the non-covalent interactions. Based on the present findings Ligand2-BzV_0Tyr, ligand3-BzV_0Ura, and ligand5-EeV_0Tyr showed good binding affinity with M^pro, whereas, ligand1-BzV_0Cys and Ligand2-BzV_0Tyr showed good binding affinity with RdRp, thus could act as potential lead compounds against SARS-CoV-2, which needs further validation studies. In particular, Ligand2-BzV_0Tyr could be more beneficial candidate with the dual target specificity for M^pro and RdRp.

Molecular insights into the inhibition mechanism of harringtonine against essential proteins associated with SARS-CoV-2 entry

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has recently posed a serious threat to global public health. Harringtonine (HT), as a small-molecule antagonist, has antiviral activity against a variety of viruses. There is evidence that HT can inhibit the SARS-CoV-2 entry into host cells by blocking the Spike protein and transmembrane protease serine 2 (TMPRSS2). However, the molecular mechanism underlying the inhibition effect of HT is largely elusive. Here, docking and all-atom molecular dynamic simulations were used to investigate the mechanism of HT against the receptor binding domain (RBD) of Spike, TMPRSS2, as well as the complex of RBD and angiotensin-converting enzyme 2 complex (RBD-ACE2). The results reveal that HT binds to all proteins primarily through hydrogen bond and hydrophobic interactions. Binding with HT influences the structural stability and dynamic motility processes of each protein. The interactions of HT with residues N33, H34 and K353 of ACE2, and residue K417 and Y453 of RBD contribute to disrupting the binding affinity between RBD and ACE2, which may hinder the virus entry into host cells. Our research provides molecular insights into the inhibition mechanism of HT against SARS-CoV-2 associated proteins, which will help for the novel antiviral drugs development.

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.

Modulatory effects of phytol on the antiemetic property of domperidone, possibly through the D2 receptor interaction pathway: in vivo and in silico studies

The current study is designed to evaluate the antiemetic effect of the diterpenoid phytol (PHY) using in vivo and in silico studies. For this, emesis was induced in 4-day-old chicks by the oral administration of copper sulfate (CuSO4.5H2O) at 50 mg/kg. To see the possible antiemetic mechanism of PHY, we used a number of reference drugs such as domperidone (80 mg/kg), ondansetron (24 mg/kg) and hyoscine (100 mg/kg) as positive controls, while the vehicle served as a negative control group. PHY was administered orally at the doses of 50 and 75 mg/kg. Both PHY and reference drugs were given alone or in combined groups to evaluate their synergistic or antagonistic effects on the chicks. Molecular docking of PHY and reference drugs was carried out against 5HT3, D2, D3, H1, NK1, and mAChRs (M1-M5) receptors for estimating binding affinity to the receptors. Drug-receptor interactions and active sites of the receptors were observed with the aid of different computational tools. The drug-likeness and pharmacokinetics of all the drugs were predicted through the SwissADME online database. The results suggest that PHY reduces the mean number of retches and increases latency dose-dependently in the birds. In the combination groups, PHY75 showed better antiemetic effects with domperidone and ondansetron. In addition, PHY exhibited the highest binding affinity with the D2 receptor (6CM4) (- 7.3 kcal/mol). In conclusion, PHY showed an antiemetic activity in chicks, possibly through the D2 receptor interaction pathway.

Identification of RdRp inhibitors against SARS-CoV-2 through E-pharmacophore-based virtual screening, molecular docking and MD simulations approaches

The outbreak of novel Coronavirus, an enduring pandemic declared by WHO, has consequences to an alarming ongoing public health menace which has already claimed several million human lives. In addition to numerous vaccinations and medications for mild to moderate COVID-19 infection, lack of promising medication or therapeutic pharmaceuticals remains a serious concern to counter the ongoing coronavirus infections and to hinder its dreadful spread. Global health emergencies have called for urgency for potential drug discovery and time is the biggest constraint apart from the financial and human resources required for the high throughput drug screening. However, computational screening or in-silico approaches appeared to be an effective and faster approach to discover potential molecules without sacrificing the model animals. Accumulated shreds of evidence on computational studies against viral diseases have revealed significance of in-silico drug discovery approaches especially in the time of urgency. The central role of RdRp in SARS-CoV-2 replication makes it promising drug target to curtain on going infection and its spread. The present study aimed to employ E-pharmacophore-based virtual screening to reveal potent inhibitors of RdRp as potential leads to block the viral replication. An energy-optimised pharmacophore model was generated to screen the Enamine REAL DataBase (RDB). Then, ADME/T profiles were determined to validate the pharmacokinetics and pharmacodynamics properties of the hit compounds. Moreover, High Throughput Virtual Screening (HTVS) and molecular docking (SP & XP) were employed to screen the top hits from pharmacophore-based virtual screening and ADME/T screen. The binding free energies of the top hits were calculated by conducting MM-GBSA analysis followed by MD simulations to determine the stability of molecular interactions between top hits and RdRp protein. These virtual investigations revealed six compounds having binding free energies of −57.498, −45.776, −46.248, −35.67, −25.15 and −24.90 kcal/mol respectively as calculated by the MM-GBSA method. The MD simulation studies confirmed the stability of protein ligand complexes, hence, indicating as potent RdRp inhibitors and are promising candidate drugs to be further validated and translated into clinics in future.

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.

Identification of potent COVID-19 main protease inhibitors by loading of favipiravir on Mg12O12 and Zn12O12 nanoclusters: an in silico strategy for COVID-19 treatment

Pandemic new severe acute respiratory syndrome coronavirus (SARS-CoV-2) virus has increased throughout the world. There is no effective treatment against this virus until now. Since its appearance in Wuhan, China in December 2019, SARS-CoV-2 becomes the largest challenge the world is opposite today, including the discovery of an antiviral drug for this virus. Several viral proteins have been prioritized as SARS-CoV-2 antiviral drug targets, among them the papain-like protease (PLpro) and the main protease (Mpro). Inhibition of these proteases would target viral replication, viral maturation and suppression of host innate immune responses. Potential candidates have been identified to show inhibitory effects against Mpro, both in biochemical assays and viral replication in cells. There are different molecules such as lopinavir and favipiravir considerably inhibit the activity of Mpro in vitro. Different studies have shown that structurally improved favipiravir and other similar compounds can inhibit SARS-CoV-2 main protease. In this work, we study the interactions between favipiravir with Mg12O12 and Zn12O12 nanoclusters by density functional theory (DFT) and quantum mechanics atoms in molecules (QMAIM) methods to summarize the ability to load favipiravir onto Mg12O12 and Zn12O12 nanoclusters. Favipiravir-Mg12O12 and favipiravir-Zn12O12 lowest structures complexes were chosen to dock inside the SARS-CoV-2 main protease by molecular docking study. The molecular docking analysis revealed that the binding affinity of Mg12O12 and Zn12O12 nanoclusters inside the Mpro receptor is larger than that of favipiravir. Also, the loading of favipiravir on the surface of Mg12O12 and Zn12O12 nanoclusters increased the binding affinity against the Mpro receptor. Subsequently, 100 ns molecular dynamics simulation of the favipiravir-Mg12O12, and favipiravir-Zn12O12 docked inside the Mpro complexes established that favipiravir-Mg12O12, forms the most stable complex with the Mpro. Further molecular mechanics Poisson Boltzmann surface area (MMPBSA) analyses using the MD trajectories also demonstrated the higher binding affinity of favipiravir-Mg12O12 inside the Mpro. In summary, this study demonstrates a new way to characterize leads for novel anti-viral drugs against SARS-CoV-2, by improving the drug ability of favipiravir via loading it on Mg12O12 and Zn12O12 nanoclusters.Communicated by Ramaswamy H. Sarma.