Computational investigation of the selectivity mechanisms of PI3Kδ inhibition with marketed idelalisib: combined molecular dynamics simulation and free energy calculation

Chemical analogue based drug design for cancer treatment targeting PI3K: integrating machine learning and molecular modeling

Cancer is a generic term for a group of disorders defined by uncontrolled cell growth and the potential to invade or spread to other parts of the body. Gene and epigenetic alterations disrupt normal cellular control, leading to abnormal cell proliferation, resistance to cell death, blood vessel development, and metastasis (spread to other organs). One of the several routes that play an important role in the development and progression of cancer is the phosphoinositide 3-kinase (PI3K) signaling pathway. Moreover, the gene PIK3CG encodes the catalytic subunit gamma (p110γ) of phosphoinositide 3-kinase (PI3Kγ), a member of the PI3K family. Therefore, in this study, PIK3CG was targeted to inhibit cancer by identifying a novel inhibitor through computational methods. The study screened 1015 chemical fragments against PIK3CG using machine learning-based binding estimation and docking to select the potential compounds. Later, the analogues were generated from the selected hits, and 414 analogues were selected, which were further screened, and as most potential candidates, three compounds were obtained: (a) 84,332, 190,213, and 885,387. The protein-ligand complex's stability and flexibility were then investigated by dynamic modeling. The 100 ns simulation revealed that 885,387 exhibited the steadiest deviation and constant creation of hydrogen bonds. Compared to the other compounds, 885,387 demonstrated a superior binding free energy (ΔG = -18.80 kcal/mol) with the protein when the MM/GBSA technique was used. The study determined that 885,387 showed significant therapeutic potential and justifies further experimental investigation as a possible inhibitor of the PIK3CG target implicated in cancer.

Evaluation of the anti-inflammatory effects of PI3Kδ/γ inhibitors for treating acute lung injury

Phosphatidylinositol 3-kinase delta (PI3Kδ) and gamma (PI3Kγ) are predominantly located in immune and hematopoietic cells. It is well-established that PI3Kδ/γ plays important roles in the immune system and participates in inflammation; hence, it could be a potential target for anti-inflammatory therapy. Currently, several PI3K inhibitors are used clinically to treat cancers with aberrant PI3K signaling; however, their role in treating acute respiratory inflammatory diseases has rarely been explored. Herein, we investigated the potential anti-inflammatory activities of several pharmacological PI3K inhibitors, including marketed drugs idelalisib (PI3Kδ), duvelisib (PI3Kδ/γ), and copanlisib (pan-PI3K with preferential α/δ) and the clinical drug eganelisib (PI3Kγ), for treating acute lung injury (ALI). In the lipopolysaccharide-induced RAW264.7 macrophage inflammatory model, the four inhibitors significantly suppressed proinflammatory cytokine expression by inhibiting the PI3K signaling pathway. Oral administration of PI3K inhibitors markedly improved lung injury in a murine model of ALI. PI3K pathway inhibition decreased inflammatory cell infiltration and totalprotein levels, as well as reduced the expression of associated lung inflammatory factors. Collectively, all four representative PI3K inhibitors exerted prominent anti-inflammatory properties, indicating that PI3K δ and/or γ inhibition could be ideal targets to treat respiratory inflammatory diseases by reducing the inflammatory response. The findings of the current study provide a new basis for utilizing PI3K inhibitors to treat acute respiratory inflammatory diseases.

Molecular modeling strategy for detailing the primary mechanism of action of copanlisib to PI3K: combined ligand-based and target-based approach

Since dysregulation of the phosphatidylinositol 3-kinase (PI3K) signaling pathway is associated with the pathogenesis of cancer, inflammation, and autoimmunity, PI3K has emerged as an attractive target for drug development. Although copanlisib is the first pan-PI3K inhibitor to be approved for clinical use, the precise mechanism by which it acts on PI3K has not been fully elucidated. To reveal the binding mechanisms and structure-activity relationship between PI3K and copanlisib, a comprehensive modeling approach that combines 3D-quantitative structure-activity relationship (3D-QSAR), pharmacophore model, and molecular dynamics (MD) simulation was utilized. Initially, the structure-activity relationship of copanlisib and its derivatives were explored by constructing a 3D-QSAR. Then, the key chemical characteristics were identified by building common feature pharmacophore models. Finally, MD simulations were performed to elucidate the important interactions between copanlisib and different PI3K subtypes, and highlight the key residues for tight-binding inhibitors. The present study uncovered the principal mechanism of copanlisib's action on PI3K at the theoretical level, and these findings might provide guidance for the rational design of pan-PI3K inhibitors.Communicated by Ramaswamy H. Sarma.

Deciphering Selectivity Mechanism of BRD9 and TAF1(2) toward Inhibitors Based on Multiple Short Molecular Dynamics Simulations and MM-GBSA Calculations

BRD9 and TAF1(2) have been regarded as significant targets of drug design for clinically treating acute myeloid leukemia, malignancies, and inflammatory diseases. In this study, multiple short molecular dynamics simulations combined with the molecular mechanics generalized Born surface area method were employed to investigate the binding selectivity of three ligands, 67B, 67C, and 69G, to BRD9/TAF1(2) with IC50 values of 230/59 nM, 1400/46 nM, and 160/410 nM, respectively. The computed binding free energies from the MM-GBSA method displayed good correlations with that provided by the experimental data. The results indicate that the enthalpic contributions played a critical factor in the selectivity recognition of inhibitors toward BRD9 and TAF1(2), indicating that 67B and 67C could more favorably bind to TAF1(2) than BRD9, while 69G had better selectivity toward BRD9 over TAF1(2). In addition, the residue-based free energy decomposition approach was adopted to calculate the inhibitor–residue interaction spectrum, and the results determined the gatekeeper (Y106 in BRD9 and Y1589 in TAF1(2)) and lipophilic shelf (G43, F44, and F45 in BRD9 and W1526, P1527, and F1528 in TAF1(2)), which could be identified as hotspots for designing efficient selective inhibitors toward BRD9 and TAF1(2). This work is also expected to provide significant theoretical guidance and insightful molecular mechanisms for the rational designs of efficient selective inhibitors targeting BRD9 and TAF1(2).

Theoretical study of myriocin-binding mechanism targeting serine palmitoyltransferase

Sphingolipids (SLs) are vital for cells as forming membrane and transducing signals. The first step for de novo biosynthesis of SLs is catalyzed by the pyridoxal-5'-phosphate (PLP)-dependent enzyme serine palmitoyltransferase (SPT), which has been proven to be a promising drug target for treating various diseases. However, there are few SPT-specific inhibitors have been identified so far. Myriocin, a natural fungal product, is confirmed as the most potent inhibitor of SPT and has been widely used, but studies of its molecular mechanism are still underway. Besides, there is no intact co-crystal structure of SPT-binding myriocin until now. Aiming to uncover the interaction mechanism between SPT- and PLP-binding myriocin at the molecular level, a systematic computational strategy was performed in this present study. Firstly, covalent docking was implemented to preliminarily predict the binding pose SPT/PLP-myriocin aldimine and its structurally similar intermediate SPT/PLP-β-ketoacid aldimine. Secondly, two binding complexes were treated as initial structures to perform molecular dynamics simulations and binding free energy calculations. The calculated docking scores and predicted binding energies were consistent with the reported bioactivities. Finally, the binding mechanism of myriocin binding with SPT was meticulously described, and the key residues making favorable contributions were highlighted. Taken together, the current study could provide some important information and valuable guidance for further rational screening, design, and modification of potent specific SPT inhibitors.

Tailored Parameterization of the LIE Method for Calculating the Binding Free Energy of Vps34-Inhibitor Complexes

Vps34 is the only isoform of the PI3K family in fungi, making this protein an attractive target to develop new treatments against pathogenic fungi. The high structural similarity between the active sites of the human and fungal Vps34 makes repurposing of human Vps34 inhibitors an appealing strategy. Nonetheless, while some of the cross-reactive inhibitors might have the potential to treat fungal infections, a safer approach to prevent undesired side effects would be to identify molecules that specifically inhibit the fungal Vps34. This study presents the parameterization of four LIE models for estimating the binding free energy of Vps34-inhibitor complexes. Two models are parameterized using a multiparametric linear regression leaving one or more free parameters, while the other two are based on the LIE-D model. All of the models show good predictive capacity (R ² > 0.7, r > 0.85) and a low mean absolute error (MAE < 0.71 kcal/mol). The current study highlights the advantages of LIE-D-derived models when predicting the weight of the different contributions to the binding free energy. It is expected that this study will provide researchers with a valuable tool to identify new Vps34 inhibitors for relevant applications such as cancer treatment and the development of new antimicrobial agents.

A multi-conformational virtual screening approach based on machine learning targeting PI3Kγ

Nowadays, more and more attention has been attracted to develop selective PI3Kγ inhibitors, but the unique structural features of PI3Kγ protein make it a very big challenge. In the present study, a virtual screening strategy based on machine learning with multiple PI3Kγ protein structures was developed to screen novel PI3Kγ inhibitors. First, six mainstream docking programs were chosen to evaluate their scoring power and screening power; CDOCKER and Glide show satisfactory reliability and accuracy against the PI3Kγ system. Next, virtual screening integrating multiple PI3Kγ protein structures was demonstrated to significantly improve the screening enrichment rate comparing to that with an individual protein structure. Last, a multi-conformational Naïve Bayesian Classification model with the optimal docking programs was constructed, and it performed a true capability in the screening of PI3Kγ inhibitors. Taken together, the current study could provide some guidance for the docking-based virtual screening to discover novel PI3Kγ inhibitors.

Targeting phosphatidylinositol 3-kinase gamma (PI3Kγ): Discovery and development of its selective inhibitors

Phosphatidylinositol 3-kinase gamma (PI3Kγ) has been regarded as a promising drug target for the treatment of advanced solid tumors, leukemia, lymphoma, and inflammatory and autoimmune diseases. However, the high level of structural conservation among the members of the PI3K family and the diverse physiological roles of Class I PI3K isoforms (α, β, δ, and γ) highlight the importance of isoform selectivity in the development of PI3Kγ inhibitors. In this review, we provide an overview of the structural features of PI3Kγ that influence γ-isoform selectivity and discuss the structure-selectivity-activity relationship of existing clinical PI3Kγ inhibitors. Additionally, we summarize the experimental and computational techniques utilized to identify PI3Kγ inhibitors. The insights gained so far could be used to overcome the main challenges in development and accelerate the discovery of PI3Kγ-selective inhibitors.