In silico identification of hydantoin derivatives: a novel natural prolyl hydroxylase inhibitor

Molecular modeling and simulation approaches to characterize potential molecular targets for burdock inulin to instigate protection against autoimmune diseases

In the current study, we utilized molecular modeling and simulation approaches to define putative potential molecular targets for Burdock Inulin, including inflammatory proteins such as iNOS, COX-2, TNF-alpha, IL-6, and IL-1β. Molecular docking results revealed potential interactions and good binding affinity for these targets; however, IL-1β, COX-2, and iNOS were identified as the best targets for Inulin. Molecular simulation-based stability assessment demonstrated that inulin could primarily target iNOS and may also supplementarily target COX-2 and IL-1β during DSS-induced colitis to reduce the role of these inflammatory mechanisms. Furthermore, residual flexibility, hydrogen bonding, and structural packing were reported with uniform trajectories, showing no significant perturbation throughout the simulation. The protein motions within the simulation trajectories were clustered using principal component analysis (PCA). The IL-1β-Inulin complex, approximately 70% of the total motion was attributed to the first three eigenvectors, while the remaining motion was contributed by the remaining eigenvectors. In contrast, for the COX2-Inulin complex, 75% of the total motion was attributed to the eigenvectors. Furthermore, in the iNOS-Inulin complex, the first three eigenvectors contributed to 60% of the total motion. Furthermore, the iNOS-Inulin complex contributed 60% to the total motion through the first three eigenvectors. To explore thermodynamically favorable changes upon mutation, motion mode analysis was carried out. The Free Energy Landscape (FEL) results demonstrated that the IL-1β-Inulin achieved a single conformation with the lowest energy, while COX2-Inulin and iNOS-Inulin exhibited two lowest-energy conformations each. IL-1β-Inulin and COX2-Inulin displayed total binding free energies of - 27.76 kcal/mol and - 37.78 kcal/mol, respectively, while iNOS-Inulin demonstrated the best binding free energy results at - 45.89 kcal/mol. This indicates a stronger pharmacological potential of iNOS than the other two complexes. Thus, further experiments are needed to use inulin to target iNOS and reduce DSS-induced colitis and other autoimmune diseases.

Desmodin and isopongachromene as potential inhibitors of cyclin-dependent kinase 5: phytoconstituents targeting anticancer and neurological therapy

Cyclin-dependent kinase 5 (CDK5) is a proline-directed serine-threonine protein kinase vital for neuronal cell cycle arrest and differentiation. It activates by binding with p35 and p39 and is important for the functioning of the nervous system. A growing body of evidence suggests that CDK5 contributes to the onset and progression of neurodegeneration and tumorigenesis and represents itself as a potential therapeutic target. Our research illustrates virtual screening of phytochemicals from the IMPPAT (Indian Medicinal Plants, Phytochemistry and Therapeutics) library to search for potential inhibitors of CDK5. Initially, the compounds from the parent library were filtered out via their physicochemical properties following the Lipinski rule of five. Then sequentially, molecular docking-based virtual screening, PAINS filter, ADMET, PASS analysis, and molecular dynamics (MD) simulation were done using various computational tools to rule out adversities that can cause hindrances in the identification of potential inhibitors of CDK5. Finally, two compounds were selected via the extensive screening showing significant binding with CDK5 ATP-binding pocket and ultimately were selected as potent ATP-competitive inhibitors of CDK5. Finally, we propose that the elucidated compounds Desmodin and Isopongachromene can be used further in the drug discovery process and act as therapeutics in the medical industry to treat certain complex diseases, including cancer and neurodegeneration.Communicated by Ramaswamy H. Sarma.

Exploring the natural products chemical space through a molecular search to discover potential inhibitors that target the hypoxia-inducible factor (HIF) prolyl hydroxylase domain (PHD)

Introduction: Hypoxia-inducible factor (HIF) prolyl hydroxylase domain (PHD) enzymes are major therapeutic targets of anemia and ischemic/hypoxia diseases. To overcome safety issues, liver failure, and problems associated with on-/off-targets, natural products due to their novel and unique structures offer promising alternatives as drug targets. Methods: In the current study, the Marine Natural Products, North African, South African, East African, and North-East African chemical space was explored for HIF-PHD inhibitors discovery through molecular search, and the final hits were validated using molecular simulation and free energy calculation approaches. Results: Our results revealed that CMNPD13808 with a docking score of -8.690 kcal/mol, CID15081178 with a docking score of -8.027 kcal/mol, CID71496944 with a docking score of -8.48 kcal/mol and CID11821407 with a docking score of -7.78 kcal/mol possess stronger activity than the control N-[(4-hydroxy-8-iodoisoquinolin-3-yl)carbonyl]glycine, 4HG (-6.87 kcal/mol). Interaction analysis revealed that the target compounds interact with Gln239, Tyr310, Tyr329, Arg383 and Trp389 residues, and chelate the active site iron in a bidentate manner in PHD2. Molecular simulation revealed that these target hits robustly block the PHD2 active site by demonstrating stable dynamics. Furthermore, the half-life of the Arg383 hydrogen bond with the target ligands, which is an important factor for PHD2 inhibition, remained almost constant in all the complexes during the simulation. Finally, the total binding free energy of each complex was calculated as CMNPD13808-PHD2 -72.91 kcal/mol, CID15081178-PHD2 -65.55 kcal/mol, CID71496944-PHD2 -68.47 kcal/mol, and CID11821407-PHD2 -62.06 kcal/mol, respectively. Conclusion: The results show the compounds possess good activity in contrast to the control drug (4HG) and need further in vitro and in vivo validation for possible usage as potential drugs against HIF-PHD2-associated diseases.

Inhibiting Cyclin-Dependent Kinase 6 by Taurine: Implications in Anticancer Therapeutics

Cyclin-dependent kinase 6 (CDK6) is linked with a cyclin partner and plays a crucial role in the early stages of cancer development. It is currently a potential drug target for developing therapeutic molecules targeting cancer therapy. Here, we have identified taurine as an inhibitor of CDK6 using combined in silico and experimental studies. We performed various experiments to find the binding affinity of taurine with CDK6. Molecular docking analysis revealed critical residues of CDK6 that are involved in taurine binding. Fluorescence measurement studies showed that taurine binds to CDK6 with a significant binding affinity, with a binding constant of K = 0.7 × 10⁷ M^-1 for the CDK6-taurine complex. Enzyme inhibition assay suggested taurine as a good inhibitor of CDK6 possessing an IC₅₀ value of 4.44 μM. Isothermal titration calorimetry analysis further confirmed a spontaneous binding of taurine with CDK6 and delineated the thermodynamic parameters for the CDK6-taurine system. Altogether, this study established taurine as a CDK6 inhibitor, providing a base for using taurine and its derivatives in CDK6-associated cancer and other diseases.

Discovering Tuberosin and Villosol as Potent and Selective Inhibitors of AKT1 for Therapeutic Targeting of Oral Squamous Cell Carcinoma

Oral squamous cell carcinoma (OSCC) is a major cause of death in developing countries because of high tobacco consumption. RAC-alpha serine-threonine kinase (AKT1) is considered as an attractive drug target because its prolonged activation and overexpression are associated with cancer progression and metastasis. In addition, several AKT1 inhibitors are being developed to control OSCC and other associated forms of cancers. We performed a screening of the IMPPAT (Indian Medicinal Plants, Phytochemistry and Therapeutics) database to discover promising AKT1 inhibitors which pass through various important filters such as ADMET (absorption, distribution, metabolism, excretion, and toxicity) properties, physicochemical properties, PAINS (pan-assay interference compounds) filters, PASS (prediction of activity spectra for substances) analysis, and specific interactions with AKT1. Molecules bearing admirable binding affinity and specificity towards AKT1 were selected for further analysis. Initially, we identified 30 natural compounds bearing appreciable affinity and specific interaction with AKT1. Finally, tuberosin and villosol were selected as potent and selective AKT1 inhibitors. To obtain deeper insights into binding mechanism and selectivity, we performed an all-atom molecular dynamics (MD) simulation and principal component analysis (PCA). We observed that both tuberosin and villosol strongly bind to AKT1, and their complexes were stable throughout the simulation trajectories. Our in-depth structure analysis suggested that tuberosin and villosol could be further exploited in the therapeutic targeting of OSCC and other cancers after further clinical validations.

Bioactive phytoconstituents as potent inhibitors of casein kinase-2: dual implications in cancer and COVID-19 therapeutics

Casein kinase 2 (CK2) is a conserved serine/threonine-protein kinase involved in hematopoietic cell survival, cell cycle control, DNA repair, and other cellular processes. It plays a significant role in cancer progression and viral infection. CK2 is considered a potential drug target in cancers and COVID-19 therapy. In this study, we have performed a virtual screening of phytoconstituents from the IMPPAT database to identify some potential inhibitors of CK2. The initial filter was the physicochemical properties of the molecules following the Lipinski rule of five. Then binding affinity calculation, PAINS filter, ADMET, and PASS analyses followed by interaction analysis were carried out to discover nontoxic and better hits. Finally, two compounds, stylopine and dehydroevodiamines with appreciable affinity and specific interaction towards CK2, were identified. Their time-evolution analyses were carried out using all-atom molecular dynamics simulation, principal component analysis and free energy landscape. Altogether, we propose that stylopine and dehydroevodiamines can be further explored in in vitro and in vivo settings to develop anticancer and antiviral therapeutics.

Showing protein–ligands interactions, electrostatic potential of CK2 bound to selected compounds, free energy landscapes of CK2-stylopine, and CK2-dehydroevodiamines complexes.

Identification of Phytoconstituents as Potent Inhibitors of Casein Kinase-1 Alpha Using Virtual Screening and Molecular Dynamics Simulations

Casein kinase-1 alpha (CK1α) is a multifunctional protein kinase that belongs to the serine/threonine kinases of the CK1α family. It is involved in various signaling pathways associated with chromosome segregation, cell metabolism, cell cycle progression, apoptosis, autophagy, etc. It has been known to involve in the progression of many diseases, including cancer, neurodegeneration, obesity, and behavioral disorders. The elevated expression of CK1α in diseased conditions facilitates its selective targeting for therapeutic management. Here, we have performed virtual screening of phytoconstituents from the IMPPAT database seeking potential inhibitors of CK1α. First, a cluster of compounds was retrieved based on physicochemical parameters following Lipinski's rules and PAINS filter. Further, high-affinity hits against CK1α were obtained based on their binding affinity score. Furthermore, the ADMET, PAINS, and PASS evaluation was carried out to select more potent hits. Finally, following the interaction analysis, we elucidated three phytoconstituents, Semiglabrinol, Curcusone_A, and Liriodenine, posturing considerable affinity and specificity towards the CK1α binding pocket. The result was further evaluated by molecular dynamics (MD) simulations, dynamical cross-correlation matrix (DCCM), and principal components analysis (PCA), which revealed that binding of the selected compounds, especially Semiglabrinol, stabilizes CK1α and leads to fewer conformational fluctuations. The MM-PBSA analysis suggested an appreciable binding affinity of all three compounds toward CK1α.

In silico identification of prolyl hydroxylase inhibitor by per-residue energy decomposition-based pharmacophore approach

Hypoxia is an effective preconditioning stimulus and many cellular responses to hypoxia are mediated through a transcription control complex termed the hypoxia-inducible factor (HIF). The stability and activation of HIF are governed by HIF prolyl-4-hydroxylases 2 (PHD2). Hence, the development of a small molecule inhibitor for prolyl hydroxylase has been suggested as a potentially useful therapeutic strategy for the treatment of oxidative/ischemic stress conditions. Thus, to unveil a novel human PHD2 inhibitor, a custom-based virtual screening was carried out to identify the potential inhibitors against PHD2 based on; (1) the per-residue energy decomposition (PRED)-based pharmacophore model, (2) molecular docking, and (3) MD approaches. The PRED analysis was performed to identify the common interaction pattern of HIF fragment (5L9B) and crystallized ligand (4JZR) to develop a relevant accurate allosteric pharmacophore model. The custom pharmacophore model (AAARR) was developed and further used to screen multiple databases. The docking was performed as a secondary strategy for screening the pharmacophore hits. Furthermore, the docked complexes were screened by molecular dynamics (MD) simulation and molecular mechanics/generalized Born surface area (MM-GBSA) based binding free energy calculations to determine the binding energy of the inhibitors and to identify crucial interaction energy fingerprint. One hit has demonstrated good binding free energy and a better binding affinity for PHD2 compared to the other four selected ligands. Thus, the results obtained from pharmacophore, docking, and MD simulations depicted that linker length and metal binding in the scaffold could be effectively used as a potent inhibitor toward human PHD2 in AD therapeutics.