Computational Exploration of Naturally Occurring Flavonoids as TGF-ß Inhibitors in Breast Cancer: Insights from Docking and Molecular Dynamics Simulations and In-vitro Cytotoxicity Study.

Breast cancer is a global health concern, demanding innovative treatments. Targeting the Transforming Growth Factor-beta (TGF-ß) signaling pathway, pivotal in breast cancer, is a promising approach. TGF-ß inhibits proliferation via G1 phase cell cycle arrest, acting as a suppressor initially, but in later stages, it promotes progression by enhancing motility, invasiveness, and metastasis formation. This study explores naturally occurring flavonoids' interactions with TGF-ß. Using molecular docking against the protein's crystal structure (PDB Id: 1PY5), Gossypin showed the highest docking score and underwent molecular dynamics simulation, revealing complex flexibility and explaining how flavonoids impede TGF-ß signaling in breast cancer. ADMET predictions adhered to Lipinski's rule of Five. Insights into flavonoid-TGF-ß binding offer a novel angle for breast cancer treatment. Flavonoids having a good docking score like gossypin, morin, luteolin and taxifolin shown potent cytotoxic effect on breast cancer cell line, MCF-7. Understanding these interactions could inspire flavonoid-based therapies targeting TGF-ß to halt breast cancer growth. These findings pave the way for personalized, targeted breast cancer therapies, offering hope against this formidable disease.

Rational design-aided discovery of novel 1,2,4-oxadiazole derivatives as potential EGFR inhibitors.

A molecular dynamics-based sampling of epidermal growth factor receptor tyrosine kinase (EGFR-TK) was carried out to search for energetically more stable protein, which was then used for molecular docking of a series of 1,2,4-oxadiazole derivatives previously reported from our laboratory. A total of 14 compounds were docked, where compounds 6a and 6b showed better binding to EGFR in silico. Further, physicochemical properties of all the compounds were calculated, which suggested that all the molecules obeyed Lipinski's rule of 5 and had favorable polar surface area and CaCO2 permeability along with the low potential for HERG inhibition. All the compounds were then screened for their ability to produce cytotoxicity in four different cell lines overexpressing EGFR (A549, HCT-116, HEPG2, MCF-7) and one EGFR negative cancer cell line (SW620); at three concentrations: 10, 1, and 0.1 µM. None of the compounds showed activity against SW620, which suggested that the compounds show cytotoxicity through inhibition of EGFR. Compounds that showed promise in this 3-concentration screen were further subjected to multiple dose-response curves to identify the IC50 values for the shortlisted eight compounds. It was encouraging to see 6a and 6b showing the best IC50 values against almost all the cell-lines which further suggests that our design protocol can be applied to optimize this lead (which are currently in the low micromolar range) to design the homologous compounds to achieve the desired potency in the nanomolar range and also to achieve selectivity across a range of kinases.

Structural Insights of PD-1/PD-L1 Axis: An In silico Approach.

BACKGROUND: Interaction of PD-1 protein (present on immune T-cell) with its ligand PD-L1 (over-expressed on cancerous cell) makes the cancerous cell survive and thrive. The association of PD-1/PD-L1 represents a classical protein-protein interaction (PPI), where receptor and ligand binding through a large flat surface. Blocking the PD-1/PDL-1 complex formation can restore the normal immune mechanism, thereby destroying cancerous cells. However, the PD-1/PDL1 interactions are only partially characterized. OBJECTIVE: We aim to comprehend the time-dependent behavior of PD-1 upon its binding with PD-L1. METHOD: The current work focuses on a molecular dynamics simulation (MDs) simulation study of apo and ligand bound PD-1. RESULTS: Our simulation reveals the flexible nature of the PD-1, both in apo and bound form. Moreover, the current study also differentiates the type of strong and weak interactions which could be targeted to overcome the complex formation. CONCLUSION: The current article could provide a valuable structural insight about the target protein (PD-1) and its ligand (PD-L1) which could open new opportunities in developing small molecule inhibitors (SMIs) targeting either PD-1 or PD-L1.

New molecular insights for 4H-1,2,4-triazole derivatives as inhibitors of tankyrase and Wnt-signaling antagonist: a molecular dynamics simulation study.

Tankyrase (TNKS) enzymes remained central biotargets to treat Wnt-driven colorectal cancers. The success of Olaparib posited the druggability of PARP family enzymes depending on their role in tumor proliferation. In this work, an MD-simulation-based comparative assessment of the protein-ligand interactions using the best-docked poses of three selected compounds (two of the designed and previously synthesized molecules obtained through molecular docking and one reported TNKS inhibitor) was performed for a 500 ns period. The PDB:ID-7KKP and 3U9H were selected for TNKS1 and TNKS2, respectively. The Molecular Mechanics Generalized Born Surface Area (MM-GBSA) based binding energy data exhibited stronger binding of compound-15 (average values of -102.92 and -104.32 kcal/mol for TNKS1 and TNKS2, respectively) as compared to compound-22 (average values of -82.99 and -85.68 kcal/mol for TNKS1 and TNKS2, respectively) and the reported compound-32 (average values of -81.89 and -74.43 kcal/mol for TNKS1 and TNKS2, respectively). Compound-15 and compound-22 exhibited comparable or superior binding to both receptors forming stable complexes when compared to that of compound-32 upon examining their MD trajectories. The key contributors were hydrophobic stacking and optimum hydrogen bonding allowing these molecules to occupy the adenosine pocket by interfacing D-loop residues. The results of bond distance analysis, radius of gyration, root mean square deviation, root mean square fluctuation, snapshots at different time intervals, LUMO-HUMO energy differences, electrostatic potential calculations, and binding free energy suggested better binding efficiency for compound-15 to TNKS enzymes. The computed physicochemical and ADMET properties of compound-15 were encouraging and could be explored further for drug development.Communicated by Ramaswamy H. Sarma.

Identification of 1,2,4-Oxadiazoles-Based Novel EGFR Inhibitors: Molecular Dynamics Simulation-Guided Identification and in vitro ADME Studies.

Background: In this work, we have identified heterocyclic derivatives with 1,2,4 oxadiazole scaffold mimicking the functions of tyrosine kinase inhibitors. Fourteen molecules that displayed the best fit were picked from the library of compounds and studied under in-silico and in-vitro conditions. Four compounds were selected for further cytotoxicity and ADME (Absorption, Distribution, Metabolism, Elimination) profiling showing IC50 (from 8-13 µM) values against EGFR positive cancer cell line (MCF7). Methods: A molecular dynamics simulation study was performed to understand the correlation of non-covalent binding energies with biological activity. The drug-like properties of the selected four compounds (7a, 7b, 7e, and 7m) were evaluated by in-vitro ADME studies. Compounds 7a, 7b, and 7m were the active compounds in the molecular dynamics simulations study. Further, EGFR binding activity was confirmed with EGFRWT and EGFRT790M kinase assay using a luminescence-based method. Results: These compounds (7a, 7b, and 7m) showed activity against EGFRWT and mutant EGFRT790M, exhibiting IC50 values of <10 and <50 micromolar, respectively. These compounds also possess moderate aqueous solubility in 40-70 µg/mL at pH 7.4 and 30-100 µg/mL at pH 4.0. Further, 7a, 7b, and 7m showed balanced lipophilicity with Log D values ranging from 1-3. They demonstrated a good correlation in Caco-2 permeability with Apparent permeability (Papp) 1 to 5 × 10-6 cm/s in comparison with 7e, which was found to be highly lipophilic (Log D >5) and showed high permeability (Papp 17 × 10-6 cm/s). Lastly, all these compounds were moderately stable in liver microsomes at alkaline pH with a half-life of 30-60 min, while at a highly acidic pH (2.0), the compounds were stable up to 15-20 min. Conclusion: Overall, in-vitro ADME results of these molecules showed good drug-like properties, which are well correlated with the in-silico ADME data, making them ideal for developing an oral drug delivery formulation.

Protein kinase C ßII in diabetic complications: survey of structural, biological and computational studies.

INTRODUCTION: PKC-ßII is a conventional isoform of PKC. It is overexpressed in hyperglycemic conditions and is known to trigger various diabetic complications, mainly cardiovascular complications and to a certain extent nephropathy, neuropathy, retinopathy etc. Selective inhibition of this enzyme will be one of the favorable approaches to treat diabetes-mellitus-related complications. Due to high sequence similarities among PKC isoforms, selective inhibition of PKC-ßII is difficult and yet to be achieved successfully. AREAS COVERED: This review discusses the studies carried out in various aspects of PKC-ßII. The biological aspects, crystal structure data, structure­activity relationship study (SAR) and in silico studies related to PKC-ßII such as homology modeling, molecular docking, molecular dynamics, quantitative structure­activity relationship (QSAR) studies and pharmacophore modeling etc. are summarized. EXPERT OPINION: PKC-ßII is a potential target for treating diabetes-related complications. Selective inhibitors of this enzyme are under clinical trials but to date, success has not been achieved. Thus, extensive research is essential in this direction; the contribution of in silico tools in designing and optimizing selective inhibitors of PKC-ßII is valuable.