MM-PBSA Analysis Research Articles

Rational Approach toward COVID-19's Main Protease Inhibitors: A Hierarchical Biochemoinformatics Analysis.

This study investigated the potential of selected compounds as inhibitors of SARS-CoV-2 Mpro through pharmacokinetic and toxicological analyses, molecular docking, and molecular dynamics simulations. In silico molecular docking simulations revealed promising ligands with favorable binding affinities for Mpro, ranging from -6.2 to -9.5 kcal/mol. Moreover, molecular dynamics simulations demonstrated the stability of protein-ligand complexes over 200 ns, maintaining protein secondary structures. MM-PBSA analysis revealed favorable interactions between ligands and Mpro, with negative binding energy values. Hydrogen bond formation capacity during molecular dynamics was confirmed, indicating consistent interactions with Mpro catalytic residues. Based on these findings, selected ligands show promise for future studies in developing COVID-19 treatments.

Natural products can be potential inhibitors of metalloproteinase II from Bacteroides fragilis to intervene colorectal cancer

Bacteroides fragilis, a gram negative and obligate anaerobe bacterium, is a member of normal gut microbiota and facilitates many essential roles being performed in human body in normal circumstances specifically in Gastrointestinal or GI tract. Sometimes, due to genetics, epigenetics, and environmental factors, Bacteroides fragilis and their protein(s) start interacting with intestinal epithelium thus damaging the lining leading to colorectal cancers (CRC). To identify these protein(s), we incorporated a novel subtractive proteomics approach in the study. Metalloproteinase II (MPII), a Bacteroides fragilis toxin (bft), was investigated for its virulence and unique pathways to demonstrate its specificity and uniqueness in pathogenicity followed by molecular docking against a set of small drug-like natural molecules to discover potential inhibitors against the toxin. All these identified inhibitor-like molecules were analyzed for their ADMET calculations and detailed physiochemical properties to predict their druggability, GI absorption, blood brain barrier and skin permeation, and others. Resultantly, a total of ten compounds with the least binding energies were obtained and were subjected to protein-compound interaction analysis. Interaction analysis revealed the most common ligand-interacting residues in MPII are His 345, Glu 346, His 339, Gly 310, Tyr 341, Pro 340, Asp 187, Phe 309, Lys 307, Ile 185, Thr 308, and Pro 184. Therefore, top three compounds complexed with MPII having best binding energies were selected in order to analyze their trajectories. RMSD, RMSF, Rg and MMPBSA analysis revealed that all compounds showed good binding and keeping the complex stable and compact throughout the simulation time in addition to all properties and qualities of being a potential inhibitor against MPII.

Exploring Quinoxalinone Derivatives as Promising Epidermal Growth Factor Receptor (EGFR) Kinase Inhibitors for Cancer Therapy

E GFR plays an extensive role in the signaling pathway such as proliferation, migration, invasion. However, the malfunctioning of the epidermal growth factor receptor (EGFR) is a crucial factor in cancer progression. The search for new kinase inhibitors is necessary due to the evolving nature of diseases and the emergence of resistant strains, prompting the need for novel compounds that can effectively target and inhibit specific kinases for therapeutic purposes. This study explores quinoxalinone derivatives, comprised of benzene and pyrazine aromatic rings, as novel potential EGFR kinase inhibitors for cancer therapy, fulfilling an essential need in targeted cancer treatments. We utilized molecular docking to examine 46 synthesized quinoxalinone derivatives for their interaction with the tyrosine kinase domain of EGFR. Fifteen compounds exhibited strong binding affinities. Notably, LS3c, MN333, and MN343 displayed significant inhibitory effects on A431 cells with half-maximum inhibitory concentration (IC50) values comparable to gefitinib, erlotinib, and afatinib in enzymatic kinase assays. Subsequent experiments revealed that LS3c and MN343 effectively inhibit EGFR-expressing cell lines with low hepatotoxicity. MN343, the most potent compound, was selected for molecular dynamic simulation to understand its interaction with the tyrosine kinase of EGFR. To investigate the interaction between MN343 and EGFR’s tyrosine kinase, MM-PBSA analysis identified key residues (L694, V702, L768, M769, and L820) crucial for interaction, highlighting the potential of MN343 as a clinical EGFR tyrosine kinase inhibitor. The need for discovering new kinase inhibitors arises from the dynamic nature of diseases and the appearance of resistant strains, emphasizing the requirement for innovative compounds capable of targeting and inhibiting specific kinases effectively for therapeutic applications.

Development of multi epitope subunit vaccines against emerging carp viruses Cyprinid herpesvirus 1 and 3 using immunoinformatics approach

Cyprinid herpesvirus is a causative agent of a destructive disease in common and koi carp (Cyprinus carpio), which leads to substantial global financial losses in aquaculture industries. Among the strains of C. herpesvirus, C. herpesvirus 1 (CyHV-1) and C. herpesvirus 3 (CyHV-3) are known as highly pathogenic to carp fishes in Europe, Asia, and Africa. To date, no effective vaccine has been developed to combat these viruses. This study aimed to develop unique multi-epitope subunit vaccines targeting the CyHV-1 and CyHV-3 using a reverse vaccinology approach. The study began with a comprehensive literature review to identify the most critical proteins, which were then subjected to in silico analyses to predict highly antigenic epitopes. These analyses involved assessing antigenicity, transmembrane topology screening, allergenecity, toxicity, and molecular docking approaches. We constructed two multi-epitope-based vaccines incorporating a suitable adjuvant and appropriate linkers. It revealed that both the vaccines are non-toxic and immunogenic. The tertiary structures of the vaccine proteins were generated, refined, and validated to ensure their suitability. The binding affinity between the vaccine constructs and TLR3 and TLR5 receptors were assessed by molecular docking studies. Molecular dynamics simulations indicated that vaccine construct V1 exhibited greater stability with both TLR3 and TLR5 based on RMSD analysis. Hydrogen bond analysis revealed a stronger binding affinity between the vaccine constructs and TLR5 compared to TLR3. Furthermore, MM-PBSA analysis suggested that both vaccine constructs exhibited a better affinity for TLR5. Considering all aspects, the results suggest that in silico development of CyHV vaccines incorporating multiple epitopes holds promise for management of diseases caused by CyHV-1 and CyHV-3. However, further in vivo trials are highly recommended to validate the efficacies of these vaccines.

Investigating the Mechanistic of Danhong Injection in Brain Damage Caused by Cardiac I/R Injury via Bioinformatics, Computer Simulation, and Experimental Validation.

Cardiac ischemia-reperfusion (I/R) injury has negative effects on the brain and can even lead to the occurrence of ischemic stroke. Clinical evidence shows that Danhong injection (DHI) protects the heart and brain following ischemic events. This study investigated the mechanisms and key active compounds underlying the therapeutic effect of DHI against brain damage induced by cardiac I/R injury. The gene expression omnibus database provided GSE66360 and GSE22255 data sets. The R programming language was used to identify the common differentially expressed genes (cDEGs). Gene ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analysis were performed, and protein-protein interaction network was constructed. Active compounds of DHI were collected from the Traditional Chinese Medicine Systems Pharmacology database. Molecular docking and molecular dynamics simulations were performed. The MMPBSA method was used to calculate the binding-free energy. The pkCSM server and DruLiTo software were used for Absorption, Distribution, metabolism, excretion, and toxicity (ADMET) analysis and drug-likeness analysis. Finally, in vitro experiments were conducted to validate the results. A total of 27 cDEGs had been identified. The PPI and enrichment results indicated that TNF-α was considered to be the core target. A total of 80 active compounds were retrieved. The molecular docking results indicated that tanshinone I (TSI), tanshinone IIA (TSIIA), and hydroxyl safflower yellow A (HSYA) were selected as core active compounds. Molecular dynamics verification revealed that the conformations were relatively stable without significant fluctuations. MMPBSA analysis revealed that the binding energies of TSI, TSIIA, and HSYA with TNF-α were -36.01, -21.71, and -14.80 kcal/mol, respectively. LEU57 residue of TNF-α has the highest contribution. TSI and TSIIA passed both the ADMET analysis and drug-likeness screening, whereas HSYA did not. Experimental verification confirmed that DHI and TSIIA reduced the expression of TNF-α, NLRP3, and IL-1β in the injured H9C2 and rat brain microvascular endothelial cells. TNF-α can be considered to be a key target for BD-CI/R. TSIIA in DHI exerts a significant inhibitory effect on the inflammatory damage of BD-CI/R, providing new insights for future drug development.

New enantioenriched β-indolyl ketones as aromatase inhibitors: Unraveling heme-ligand interactions by MD simulationand MMPBSA analysis.

A series of enantioenriched β-indolyl ketones as aromatase inhibitors (AI) is synthesized through the Michael-type Friedel-Crafts alkylation of indole. A highly efficient bifunctionalized amino catalyst is developed to access structurally diverse β-indolyl ketones in high yields (up to 91%) and excellent enantioselectivity (enantiomeric ratio up to 98:2). All the synthesized compounds demonstrated promising aromatase inhibitory potential, where ortho-substituted analogs (3c and 3e) were found most active with IC50 values of 0.68 and 0.90 µM, respectively. Both of these compounds exhibited significant cytotoxicity (IC50 = 0.34 and 0.37 µM) against the MCF-7 breast cancer cell line in the (3-[4,5-dimethylthiazol-2-yl]-2,5-diphenyl tetrazolium bromide) assay. Molecular docking studies of the synthesized compounds demonstrate favorable binding interactions with the estrogens controlling CYP19A1 (3EQM) and metabolizing CYP3A4 (5VCC) enzymes. Molecular dynamic (MD) simulation analysis revealed the essentiality of heme-ligand interactions to build a stable protein-ligand complex. An average root mean square deviation of 0.35 nm observed during a 100-ns MD simulation and binding free energy in the range of -190 to -227 kJ/mol calculated by g_mmpbsa analysis authenticated the stability of the 3c-3EQM complex. ADMET and drug-likeness parameters supported the suitability of these indole derivatives as the drug lead to develop potent inhibitors for estrogen-dependent breast cancer.

Deciphering the interaction mechanism between soy protein isolate and fat-soluble anthocyanin on experiments and molecular simulations

In this work, the acylated anthocyanin (Ca-An) was prepared by enzymatic modification of black rice anthocyanin with caffeic acid, and the binding mechanism of Ca-An to soybean protein isolate (SPI) was investigated by experiments and computer simulation to expand the potential application of anthocyanin in food industry. Multi-spectroscopic studies revealed that the stable binding of Ca-An to SPI induced the folding of protein polypeptide chain, which transformed the secondary structure of SPI trended to be flexible. The microenvironment of protein was transformed from hydrophobic to hydrophilic, while tyrosine played dominant role in quenching process. The binding sites and forces of the complexes were determined by computer simulation for further explored. The protein conformation of the 7S and 11S binding regions to Ca-An changed, and the amino acid microenvironment shifted to hydrophilic after binding. The results showed that more non-polar amino acids existed in the binding sites, while in binding process van der Waals forces and hydrogen bonding played a major role hydrophobicity played a minor role. Based on MM-PBSA analysis, the binding constants of 7S-Ca-An and 11S-Ca-An were 0.518 × 106 mol−1 and 5.437 × 10−3 mol−1, respectively. This information provides theoretical guidance for further studying the interaction between modified anthocyanins and biomacromolecules.

Comparative study on interactions of phillyrin and phillygenol with lysozyme: Spectroscopy, differential scanning calorimetry and molecular modeling approaches

Phillyrin (PR) and its derivative phillygenol (PG) are major active extracts of the Chinese medicine Forsythiae Fructus. Lysozyme, an important carrier protein, has a variety of physiological along with pharmaceutical functions. Herein, the interactions of PR and PG with hen egg white lysozyme (HEWL) are investigated in details by spectroscopic analyses, differential scanning calorimetry (DSC) tests, and molecular modeling approaches. Fluorescence measurement and UV–visible absorption difference spectroscopy reveal that an “untypical static quenching” takes place pointing towards the formation of ground state complex. The binding constants of HEWL-PG are greater than that of HEWL-PR. Moreover, the major driving force for the binding process is elucidated as hydrophobic interactions for both cases. Circular dichroism spectroscopy and three-dimensional fluorescence indicate that PR and PG decrease the helical structure of the HEWL and hydrophobicity surrounding Trp and Tyr residues. The thermal stability of HEWL decreases in the presence of PR and PG by DSC. Docking study infers that both PR and PG could bind to the hydrophobic cavity of HEWL and supports the higher binding of PG than PR. Molecular dynamics study illustrates that PR has higher relative position change, which may attribute to its lower binding ability. In terms of the radius of gyration (Rg), PG causes the protein structure to swell even more. The average hydrogen bonds number of HEWL-PG is less than that of HEWL-PR and more changes of the residual fluctuation are observed in HEWL-PG system. Furthermore, MM-PBSA analysis exhibits that HEWL-PG has lower binding energy. Collectively, these results support that both PR and PG could bind to HEWL, but PG possesses higher affinity which leads to stronger secondary structural changes of HEWL than PR. These findings are beneficial to comprehensively understand the binding characteristic of PR and PG with HEWL, providing experimental basis for pharmacokenitics of PR and PG.

Targeting human inosine 5’ monophosphate dehydrogenase type 2 for anti-dengue lead identification – a computational approach

Dengue is a rapidly evolving arboviral disease that mainly affects tropical and subtropical regions of the world. The lack of therapeutic drugs and effective vaccines suggests that further resources need to be investigated. The effectiveness of the existing dengue vaccine is improbable as its efficacy depends on prior exposure to the dengue virus(DENV). Although the mechanism underlying the action of bioactive compounds to limit viral replication is less studied and still needs to be further explored, medicinal plants are excellent alternatives to combat DENV infection. In the current study, an in silico screening of phytochemicals from Annona reticulata Linn. against human Impdh2 was performed using Autodock Vina. Daucosterol (−9.0 kcal/mol) and Kaurenoic acid (−8.5 kcal/mol) were chosen as the top hits based on molecular interaction analysis. The hits were further exposed to pharmacokinetics and toxicity properties to determine their drug-like parameters. Molecular dynamics simulation studies of the Impdh2-top hits were carried out to investigate their kinetic behaviour and structural stabilities. The binding free energies of the Impdh2-hit complexes were determined using MM-PBSA analysis. According to the overall conclusions of the study, Daucosterol showed good binding affinity and high structural stability to the binding site residues of the target, therefore it is recommended as a lead compound against dengue. Communicated by Ramaswamy H. Sarma

Guideline directed medical therapy induced nephrotoxicity in HFrEF patients; an insight to its mechanism

Guideline Directed Medical Therapy (GDMT) has been the standard pharmacotherapy for the treatment of Heart Failure patients with reduced Ejection Fraction (HFrEF) recommended by the European Society of Cardiology (ESC). However, patients on GDMT are likely to possess nephrotoxicity as an adverse effect. We utilized multiple system biology tools like ADVER-Pred, gene enrichment analysis, molecular docking, molecular dynamic simulations, and MMPBSA analysis to predict a possible molecular mechanism of how selected combinations of GDMT may cause nephrotoxicity. As per the ACC/AHA/ESC guidelines, we categorized the drugs as category 1 including β-blockers (BB), angiotensin receptor blockers (ARB), and sodium-glucose cotransporter-2 inhibitors (SGLT2I), category 2 includes BB’s, SGLT2I, and angiotensin receptor-neprilysin inhibitors (ARNI), and category 3 includes BB’s, SGLT2I, and angiotensin-converting enzyme (ACE) inhibitors. Enrichment analysis predicted category 2 drugs to possess the highest number of proteins to be involved in the development of nephrotoxicity i.e. 79.41%. The targets HBA1, CBR1, ATG5, and SLC6A3 were the top hub genes with an edge count of 7 followed by GPX1 with an edge count of 6. Molecular docking studies revealed candesartan-SLC6A3 to possess the highest binding affinity of −10.2 kcal/mol. In addition, simulation studies displayed empagliflozin-CBR1 to possess the highest stability followed by candesartan-ATG5. A combination of β-blockers, ARBs, and SGLT2I are predicted to likely possess nephrotoxicity which may be due to the modulation of HBA1, CBR1, ATG5, and GPX1. In conclusion, candesartan and empagliflozin are most likely to cause nephrotoxicity via the modulation of HBA1, CBR1, ATG5, and GPX1. Communicated by Ramaswamy H. Sarma

Identification of bioactive compounds of Zanthoxylum armatum as potential inhibitor of pyruvate kinase M2 (PKM2): Computational and virtual screening approaches

PKM2 (Pyruvate kinase M2) is the isoform of pyruvate kinase which is known to catalyse the last step of glycolysis that is responsible for energy production. This specific isoform is known to be highly expressed in certain cancerous conditions. Considering the role of this protein in various cancer conditions, we used PKM2 as a target protein to identify the potential compounds against this target. In this study, we have examined 96 compounds of Zanthoxylum armatum using an array of computational and in silico tools. The compounds were assessed for toxicity then their anticancer potential was predicted. The virtual screening was done with molecular docking followed by a detailed examination using molecular dynamics simulation. The majority of the compounds showed a higher probability of being antineoplastic. Based on toxicity, predicted anticancer potential, binding affinity, and binding site, three compounds (nevadensin, asarinin, and kaempferol) were selected as hit compounds. The binding energy of these compounds with PKM2 ranged from −7.7 to −8.3 kcal/mol and all hit compounds interact at the active site of the protein. The selected hit compounds formed a stable complex with PKM2 when simulated under physiological conditions. The dynamic analysis showed that these compounds remained attached to the active site till the completion of molecular simulation. MM-PBSA analysis showed that nevadensin exhibited a higher affinity towards PKM2 compared to asarinin and kaempferol. These compounds need to be assessed properties in vivo and in vitro to validate their efficacy.

Identification of alpha-glucosidase enzyme inhibitors from phytochemicals via integrated deep learning, molecular docking, molecular dynamics simulation, and MMPBSA analysis

Alpha-glucosidase is a crucial enzyme involved in carbohydrate metabolism, and inhibiting its activity holds promise for treating Type 2 Diabetes Mellitus (T2DM). This study aimed to identify potential drugs that can selectively target alpha-glucosidase's essential enzyme, Maltase-glucoamylase, to discover potent antidiabetic compounds. A phytochemical library containing 2000 compounds from various plants was subjected to screening using a deep learning approach, molecular docking, ADMET, and DruliTo against the Maltase-glucoamylase enzyme. The compounds exhibiting strong binding affinity with the enzyme were further analyzed through Molecular Dynamics Simulation (MDS). MDS analysis, including (RMSD) Root Mean Square Deviation, Root Mean Square Fluctuation (RMSF), Radius of Gyration, Hydrogen-bond, and PCA, indicates that Carpaine, Pseudocarpaine, Pollinastanol, and Annosquamosin C have favorable interactions with alpha-glucosidase. Moreover, analysis of free binding energy substantiates the stability of these interactions, as evidenced by the binding energies of α-Glucosidase-Reference (-4.6 kJ/mol), α-Glucosidase-Pollinastaol (-4.7 kJ/mol), α-Glucosidase-Annosquamosoin C (-5.6 kJ/mol), α-Glucosidase-Pseudocarpaine (-8.0 kJ/mol), and α-Glucosidase-Carpaine (-10.2 kJ/mol). Based on these results, we conclude that these phytochemicals possess significant potential for inhibiting the alpha-glucosidase enzyme. The findings of this study suggest that these phytochemicals may exhibit promising inhibitory activity against α-Glucosidase and should be further investigated in vitro and in vivo experiments as potential natural resources for developing effective antidiabetic drugs.

Unveiling the azo-reductase mechanism in Pseudomonas putida for efficient decolorization of textile Reactive dyes: an in-silico study

The textile industry utilizing affordable azo dyes is a high threat to aquatic life and causes environmental problems due to their toxicity. Biodegradation of azo dyes employing microbes and enzymes has proved to be an efficient method for treating industrial effluent. This study used the novel microbial consortium to decolorize reactive azo dyes (Reactive Red 120; Reactive Black 5 and Reactive Blue 13), and its azo-reductase activity was evaluated. The metagenomic analysis of the consortium identified azo-reductase-producing bacterial species. The molecular docking revealed that PpAzoR from Pseudomonas putida had the highest binding affinities for all the three dyes such as Reactive Black 5 (−9.3 kcal/mol), Reactive Blue 13 (−9.8 kcal/mol) and Reactive Red 120 (−10.7 kcal/mol). The structural rigidity and stability of the docked complex were confirmed through MD simulations evaluated across multiple descriptors from the simulation trajectories. Further, MMPBSA analysis validated the results that binding of the ligands, i.e. dye molecules Reactive Black (RB5), Reactive Blue (RB13) and Reactive Red (RR120) binding with the Azoreductase (PpAzoR) to the screened Azo-dyes was spontaneous. Based on molecular dynamics simulations for 100 ns, RR 120 showed the highest binding affinity (−411.336 ± 46.799 KJ/mol), followed by RB5 (−288.012 ± 33.371 KJ/mol). The dyes (RR120 and RB5) exhibited stable interactions with the target azoreductase (PpAzoR). The present study provides insights that PpAzoR shows the highest decolorization potency, which could be interpreted as a potential dye-degrading protein based on dye-degrading assay findings. Communicated by Ramaswamy H. Sarma

An in silico investigation on the binding site preference of PD-1 and PD-L1 for designing antibodies for targeted cancer therapy.

Cancer control and treatment remain a significant challenge in cancer therapy and recently immune checkpoints has considered as a novel treatment strategy to develop anti-cancer drugs. Many cancer types use the immune checkpoints and its ligand, PD-1/PD-L1 pathway, to evade detection and destruction by the immune system, which is associated with altered effector function of PD-1 and PD-L1 overexpression on cancer cells to deactivate T cells. In recent years, mAbs have been employed to block immune checkpoints, therefore normalization of the anti-tumor response has enabled the scientists to develop novel biopharmaceuticals. In vivo affinity maturation of antibodies in targeted therapy has sometimes failed, and current experimental methods cannot accommodate the accurate structural details of protein-protein interactions. Therefore, determining favorable binding sites on the protein surface for modulator design of these interactions is a major challenge. In this study, we used the in silico methods to identify favorable binding sites on the PD-1 and PD-L1 and to optimize mAb variants on a large scale. At first, all the binding areas on PD-1 and PD-L1 have been identified. Then, using the RosettaDesign protocol, thousands of antibodies have been generated for 11 different regions on PD-1 and PD-L1 and then the designs with higher stability, affinity, and shape complementarity were selected. Next, molecular dynamics simulations and MM-PBSA analysis were employed to understand the dynamic, structural features of the complexes and measure the binding affinity of the final designs. Our results suggest that binding sites 1, 3 and 6 on PD-1 and binding sites 9 and 11 on PD-L1 can be regarded as the most appropriate sites for the inhibition of PD-1-PD-L1 interaction by the designed antibodies. This study provides comprehensive information regarding the potential binding epitopes on PD-1 which could be considered as hotspots for designing potential biopharmaceuticals. We also showed that mutations in the CDRs regions will rearrange the interaction pattern between the designed antibodies and targets (PD-1 and PD-L1) with improved affinity to effectively inhibit protein-protein interaction and block the immune checkpoint.

Targeting EGFR allosteric site with marine-natural products of Clathria Sp.: A computational approach

The EGFR-C797S resistance mutation to third-generation drugs has been overcome by fourth-generation inhibitors, allosteric inhibitors, namely EAI045 and has reached phase 3 clinical trials, so the Allosteric Site is currently an attractive target for development. In this study, researchers are interested in knowing the activity of metabolite compounds from marine natural ingredients Clathria Sp. against the Allosteric Site of EGFR computationally. The methods used include molecular docking using Autodock4 software and Molecular Dynamics simulation performed using GROMACS software. The research began with the preparation of metabolite samples from Clathria Sp. through the KnapSack database site and the preparation of EGFR receptors that have been complexed with allosteric inhibitors, namely proteins with PDB code 5D41. Each compound was docked to the Allosteric Site of the natural ligand and then molecular dynamics simulations were performed on the compound with the best docking energy compared to the natural ligand. From the docking results, the Clathrin_A compound showed the lowest binding energy compared to other metabolites, and the value was close to the natural ligand. Then from the molecular dynamics results, the clathrin_A compound shows good stability and resembles the natural ligand, which is analyzed through RMSD, RMSF, SASA, Rg, and PCA, and shows the binding free energy from MMPBSA analysis which is close to the natural ligand. It can be concluded, Clathrin_A compound has potential as an allosteric inhibitor.

Exploring Therapeutic Potential of Nutraceutical Compounds from Propolis on MAPK1 Protein Using Bioinformatics Approaches as Anti-Coronavirus Disease 2019 (COVID-19)

This study explores the potential of propolis, a natural substance, as a gene therapy for treating COVID-19. Despite the advent of COVID-19 vaccines, their side effects pose new health challenges. Utilizing network pharmacology, this research identifies propolis compounds through various databases and assesses their ability to target proteins associated with COVID-19. MAPK1 emerges as a potential therapeutic target, and molecular docking reveals Broussoflavonol F, Glyasperin A, and Sulabiroins as promising compounds with strong binding affinities, i.e.,- 9.0, -9.0, and -8.8 kcal/mol, respectively, exceeding the native ligand (-7.2 kcal/mol). Molecular Dynamics displays stable complex behavior, with backbone RMSD values consistently below 4 Angstroms and RMSF simulations showing minimal fluctuations within ±2 Angstroms error. Moreover, MM-PBSA analysis further supports the strong binding of Broussoflavonol F, Glyasperin A, and Sulabiroins A, with relative binding energies of -122.82±89.65, 131.48±95.39, and -155.97±111,37 kJ/mol, respectively. These results indicate that propolis has potential as an anti-COVID-19 agent, primarily through inhibiting the MAPK1 pathway. However, further research is needed to validate these results and develop practical applications for COVID-19 therapy. This study underscores the significance of network pharmacology and computational models in understanding propolis mechanisms, offering potential directions for future research and treatment strategies against COVID-19.

Identification of novel inhibitors against VP40 protein of Marburg virus by integrating molecular modeling and dynamics approaches

Marburg virus (MV) is a highly etiological agent of haemorrhagic fever in humans and has spread across the world. Its outbreaks caused a 23–90% human death rate. However, there are currently no authorized preventive or curative measures yet. VP40 is the MV matrix protein, which builds protein shell underneath the viral envelope and confers hallmark filamentous. VP40 alone is able to induce assembly and budding of filamentous virus-like particles (VLPs), which resemble authentic virions. As a result, this research is credited with clarifying the function of VP40 and leading to the discovery of new therapeutic targets effective in combating MV disease (MVD). Virtual screening, molecular docking and molecular dynamics (MD) simulation were used to find the putative active chemicals based on a 3D pharmacophore model of the protein’s active site cavity. Initially, andrographidine-C, a potent inhibitor was selected for the development of the pharmacophore model. Later, a library of 30,000 compounds along with the andrographidine-C was docked against VP40 protein. Three best hits including avanafil, diuvaretin and macrourone were subjected to further MD simulation analysis, as these compounds had better binding affinities as compared to andrographidine-C. Furthermore, throughout the 100 ns simulations, the back bone of VP40 protein in presence of avanafil, diuvaretin and macrourone remained stable which was further validated by MM-PBSA analysis. Additionally, all of these compounds depict maximum drug-like properties. The predicted drugs based on the ligand, avanafil, diuvaretin and macrourone could be exploited and developed as an alternative or complementary therapy for the treatment of MVD. Communicated by Ramaswamy H. Sarma

Discovery of Potential Prolyl-tRNA Synthetase Allosteric Inhibitor Through Virtual Screening and In Vitro Assay against Plasmodium falciparum

Objectives: This study aimed to identify novel antimalarial compounds based on allosteric inhibitor of prolyl-tRNA synthetase using hierarchical virtual screening. Materials and Methods: Pharmacophore model was designed initially, based on the structure-activity relationships data between several pyrazole-urea analogues and their IC50 enzymatic value. The model obtained was applied to screen ZINC15 database, after which followed by drug-likeness, toxicophore, and PAINS filter. The hit compounds were docked against P. falciparum prolyl-tRNA synthetase enzyme, using validated docking method. The resulting docking poses were ranked based on the docking score and re-evaluated based on the pharmacophore criteria. Top five compounds were obtained from this step and then evaluated using molecular dynamics simulation to verify its stability and hydrogen bond dynamics over 50 nanoseconds. MM-PBSA analysis was also performed to estimate their binding free energy. Ultimately, their potential bioactivity as antimalarial candidates have been verified against 3D7 strain. Results: The results showed that all five compounds obtained from virtual screening possess micromolar potency in vitro. Two compounds (ZINC 1029449 and ZINC1029453), yield high antimalarial activity (0.44 and 0.72 μM, respectively) Conclusions: Overall, the virtual screening approach has successfully produced lead compounds which can be further optimized to be antimalarial agents.

System biology mediated assessment of molecular mechanism for sinapic acid against breast cancer: via network pharmacology and molecular dynamic simulation

Sinapic acid is a hydroxycinnamic acid widespread in the plant kingdom, known to be a potent anti-oxidant used for the treatment of cancer, infections, oxidative stress, and inflammation. However, the mode of action for its chemotherapeutic properties has yet not been unleashed. Hence, we aimed to identify potential targets to propose a possible molecular mechanism for sinapic acid against breast cancer. We utilized multiple system biology tools and databases like DisGeNET, DIGEP-Pred, Cytoscape, STRING, AutoDock 4.2, AutoDock vina, Schrodinger, and gromacs to predict a probable molecular mechanism for sinapic acid against breast cancer. Targets for the disease breast cancer, were identified via DisGeNET database which were further matched with proteins predicted to be modulated by sinapic acid. In addition, KEGG pathway analysis was used to identify pathways; a protein-pathway network was constructed via Cytoscape. Molecular docking was performed using three different algorithms followed by molecular dynamic simulations and MMPBSA analysis. Moreover, cluster analysis and gene ontology (GO) analysis were performed. A total of 6776 targets were identified for breast cancer; 95.38% of genes predicted to be modulated by sinapic acid were common with genes of breast cancer. The ‘Pathways in cancer’ was predicted to be modulated by most umber of proteins. Further, PRKCA, CASP8, and CTNNB1 were predicted to be the top 3 hub genes. In addition, molecular docking studies revealed CYP3A4, CYP1A1, and SIRT1 to be the lead proteins identified from AutoDock 4.2, AutoDock Vina, and Schrodinger suite Glide respectively. Molecular dynamic simulation and MMPBSA were performed for the complex of sinapic acid with above mentioned proteins which revealed a stable complex throughout simulation. The predictions revealed that the mechanism of sinapic acid in breast cancer may be due to regulation of multiple proteins like CTNNB1, PRKCA, CASP8, SIRT1, and cytochrome enzymes (CYP1A1 & CYP3A4); the majorly regulated pathway was predicted to be ‘Pathways in cancer’. This indicates the rationale for sinapic acid to be used in the treatment of breast cancer. However, these are predictions and need to be validated and looked upon in-depth to confirm the exact mechanism of sinapic acid in the treatment of breast cancer; this is future scope as well as a drawback of the current study.

An in-silico approach to identify bioactive phytochemicals from Houttuynia cordata Thunb. As potential inhibitors of human glutathione reductase

Cellular infections are central to the etiology of various diseases, notably cancer and malaria. Counteracting cellular oxidative stress via the inhibition of glutathione reductase (GR) has emerged as a promising therapeutic strategy. Houttuynia cordata, a medicinal plant known for its potent antioxidant properties, has been the focus of our investigation. In this study, we conducted comprehensive in silico analyses involving the phytochemical constituents of H. cordata to identify potential natural GR inhibitors. Our methodological approach encompassed multiple in silico techniques, including molecular docking, molecular dynamics simulations, MMPBSA analysis, and dynamic cross-correlation analysis. Out of 13 docked phytochemicals, Quercetin, Quercitrin, and Sesamin emerged as particularly noteworthy due to their exceptional binding affinities for GR. Notably, our investigation demonstrated that Quercetin and Sesamin exhibited promising outcomes compared to the well-established pharmaceutical agent N-acetylcysteine (NAC). Molecular dynamics analyses provided insights into the ability of these phytochemicals to induce structural compaction and stabilization of the GR protein, as evidenced by changes in radius of gyration and solvent-accessible surface area. Moreover, MMPBSA analysis highlighted the crucial roles of specific residues, namely Gly27, Gly28, Ser51, His52, and Val61, in mediating essential interactions with these phytochemicals. Furthermore, an assessment of Absorption, Distribution, Metabolism, Excretion, and Toxicity (ADME-Tox) profiles underscored the favourable drug-like attributes of these phytochemicals. Thus, the current findings underscore the immense potential of Houttuynia cordata phytochemicals as potent antioxidants with the capacity to combat a spectrum of maladies, including malaria and cancer. This study not only unveils novel therapeutic avenues but also underscores the distinctive outcomes and paramount significance of harnessing H. cordata phytochemicals for their efficacious antioxidant properties. Communicated by Ramaswamy H. Sarma