Binding Selectivity Studies of Phosphoinositide 3-Kinases Using Free Energy Calculations

Targeting the PI3K/AKT signaling pathway in anticancer research: a recent update on inhibitor design and clinical trials (2020-2023)

INTRODUCTION

Recent years have witnessed great achievements in drug design and development targeting the phosphatidylinositol 3-kinase/protein kinase-B (PI3K/AKT) signaling pathway, a pathway central to cell growth and proliferation. The nearest neighbor protein-protein interaction networks for PI3K and AKT show the interplays between these target proteins which can be harnessed for drug discovery. In this review, we discuss the drug design and clinical development of inhibitors of PI3K/AKT in the past three years. We review in detail the structures, selectivity, efficacy, and combination therapy of 35 inhibitors targeting these proteins, classified based on the target proteins. Approaches to overcoming drug resistance and to minimizing toxicities are discussed. Future research directions for developing combinational therapy and PROTACs of PI3K and AKT inhibitors are also discussed.

AREA COVERED

This review covers clinical trial reports and patent literature on inhibitors of PI3K and AKT published between 2020 and 2023.

EXPERT OPINION

To address drug resistance and drug toxicity of inhibitors of PI3K and AKT, it is highly desirable to design and develop subtype-selective PI3K inhibitors or subtype-selective AKT1 inhibitors to minimize toxicity or to develop allosteric drugs that can form covalent bonds. The development of PROTACs of PI3Kα or AKT helps to reduce off-target toxicities.

Binding and selectivity studies of phosphatidylinositol 3-kinase (PI3K) inhibitors

Overexpression of the Phosphatidylinositol 3-kinase (PI3K) proteins have been observed in cancer cells. Targeting the phosphatidylinositol 3-kinase (PI3K) signaling transduction pathway by inhibition of the PI3K substrate recognition sites has been proved to be an effective approach to block cancer progression. Many PI3K inhibitors have been developed. Seven drugs have been approved by the US FDA with a mechanism of targeting the phosphatidylinositol 3-kinase/protein kinase-B/mammalian target of rapamycin (PI3K/AKT/mTOR) signaling pathway. In this study, we used docking tools to investigate selective binding of ligands toward four different subtypes of PI3Ks (PI3Kα, PI3Kβ, PI3Kγ and PI3Kδ). The affinity predicted from both the Glide dock and the Movable-Type (MT)-based free energy calculations agreed well with the experimental data. The validation of our predicted methods with a large dataset of 147 ligands showed very small mean errors. We identified residues that may dictate the subtype-specific binding. Particularly, residues Asp964, Ser806, Lys890 and Thr886 of PI3Kγ might be utilized for PI3Kγ-selective inhibitor design. Residues Val828, Trp760, Glu826 and Tyr813 may be important for PI3Kδ-selective inhibitor binding.

Mechanistic Investigation of Lysine-Targeted Covalent Inhibition of PI3Kδ via ONIOM QM:QM Computations

Phosphoinositide 3-kinase (PI3K) enzymes are important drug targets, especially in oncology, and several inhibitors are currently under investigation in clinical trials for the treatment of lymphocytic leukemia, follicular lymphoma, breast, thyroid, colorectal, and lung cancer. Targeted covalent inhibitors hold significant promise for drug discovery research especially for kinases. Targeting the lysine residues attracts attention as a new strategy in designing targeted covalent inhibitors, since the lysine residue provides several advantages over the traditional cysteine residue. Recently, new highly selective covalent inhibitors of PI3Kδ with activated ester warheads, targeting the conserved Lys779 residue, were reported. Based on the observed kinetics, a covalent inhibition mechanism was proposed, but the atomistic details of the reaction are still not understood. Therefore, in the present work, we have conducted quantum chemical ONIOM M06-2X/6-31+G(d,p):PM6 calculations on the active site cluster structure of PI3Kδ to elucidate the microscopic details of the mechanism of the aminolysis reaction between Lys779 and the ester inhibitors. Our calculations clearly discriminate the noncovalent methyl ester inhibitor and the covalent inhibitors with activated phenolic esters. For the representative p-NO₂, p-F, p-H, and p-OCH₃ phenolic esters, the Gibbs free energy profiles of alternative mechanistic paths through either Asp782 or Asp911 demonstrate the modulatory role of active site aspartate residues. The most plausible path alters depending on the electron-withdrawing/donating nature of the p-substituted phenolate leaving group. Inhibitors with sufficiently strong electron-withdrawing group prefer direct dissociation of the leaving group from the tetrahedral zwitterion intermediate, while the ones with electron-donating group favor the formation of a neutral tetrahedral intermediate prior to the dissociation. The relative Gibbs free energy barriers of p-NO₂ < p-F < p-H < p-OCH₃ substituted phenyl esters display the same qualitative trend as the experimentally measured k_inact/K_I values. Our results provide in depth insight into the mechanism, which can pave the way for optimizing the inhibitor efficiency.

Molecular insight into isoform specific inhibition of PI3K-α and PKC-η with dietary agents through an ensemble pharmacophore and docking studies

Dietary compounds play an important role in the prevention and treatment of many cancers, although their specific molecular mechanism is not yet known. In the present study, thirty dietary agents were analyzed on nine drug targets through in silico studies. However, nine dietary scaffolds, such as silibinin, flavopiridol, oleandrin, ursolic acid, α-boswellic acid, β-boswellic acid, triterpenoid, guggulsterone, and oleanolic acid potentially bound to the cavity of PI3K-α, PKC-η, H-Ras, and Ras with the highest binding energy. Particularly, the compounds silibinin and flavopiridol have been shown to have broad spectrum anticancer activity. Interestingly, flavopiridol was embedded in the pockets of PI3K-α and PKC-η as bound crystal inhibitors in two different conformations and showed significant interactions with ATP binding pocket residues. However, complex-based pharmacophore modeling achieved two vital pharmacophoric features namely, two H-bond acceptors for PI3K-α, while three are hydrophobic, one cat-donor and one H-bond donor and acceptor for PKC-η, respectively. The database screening with the ChemBridge core library explored potential hits on a valid pharmacophore query. Therefore, to optimize perspective lead compounds from the hits, which were subjected to various constraints such as docking, MM/GBVI, Lipinski rule of five, ADMET and toxicity properties. Henceforth, the top ligands were sorted out and examined for vital interactions with key residues, arguably the top three promising lead compounds for PI3K-α, while seven for PKC-η, exhibiting binding energy from - 11.5 to - 8.5 kcal mol^-1. Therefore, these scaffolds could be helpful in the development of novel class of effective anticancer agents.

Phosphatidylinositol 3-kinase (PI3K) inhibitors: a recent update on inhibitor design and clinical trials (2016-2020)

Introduction: The phosphatidylinositol 3-kinase/protein kinase-B/mammalian target of rapamycin (PI3K/AKT/mTOR) signaling pathway plays a central role in regulating cell growth and proliferation and thus has been considered as effective anticancer drug targets. Many PI3K inhibitors have been developed and progressed to various stages of clinical trials, and some have been approved as anticancer treatment. In this review, we discuss the drug design and clinical development of PI3K inhibitors over the past 4 years. We review the selectivity and potency of 47 PI3K inhibitors. Structural determinants for increasing selectivity toward PI3K subtype-selectivity or mutant selectivity are discussed. Future research direction and current clinical development in combination therapy of inhibitors involved in PI3Ks are also discussed.Area covered: This review covers clinical trial reports and patent literature on PI3K inhibitors and their selectivity published between 2016 and 2020.Expert opinion: To PI3Kα mutants (E542K, E545K, and H1047R), it is highly desirable to design and develop mutant-specific PI3K inhibitors. It is also necessary to develop subtype-selective PI3Kα inhibitors to minimize toxicity. To reduce drug resistance and to improve efficacy, future studies should include combination therapy of PI3K inhibitors with existing anticancer drugs from different pathways.

Docking Studies and Antiproliferative Activities of 6-(3-aryl-2-propenoyl)-2(3H)- benzoxazolone Derivatives as Novel Inhibitors of Phosphatidylinositol 3-Kinase (PI3Kα)

BACKGROUND

Cancer is a life-threatening group of diseases and universally, the second main cause of death. The design and development of new scaffolds targeting selective cancer cells are considered a promising goal for cancer treatment.

AIMS AND OBJECTIVE

Chalcone derivatives; 6-(3-aryl-2-propenoyl)-2(3H)-benzoxazolone, were previously prepared and evaluated against the oral cavity squamous cell carcinoma cell line, HSC-2, and were reported to have remarkably high tumor selectivity. The aim of this study was to further investigate the anticancer activities of the chalcone derivatives against human colon cancer cells with a possible elucidation of their mechanism of action.

METHODS

Computational studies were conducted to explore the potential interaction of the synthesized molecules with the phosphatidylinositol-4,5-bisphosphate 3-kinaseα (PI3Kα). Biological evaluation of the antiproliferative activities associated with compounds 1-23 was carried out against the colon cancer cell line, HCT116. Lactate Dehydrogenase (LDH) activity was measured to study necrosis, while the caspase-3 activation and DNA measurements were used to evaluate apoptosis in the treated cells.

RESULTS

Glide studies against PI3Kα kinase domain demonstrated that the 6-(3-aryl-2-propenoyl)-2(3H)- benzoxazolone scaffold forms H-bond with K802, Y836, E849, V851, N853, Q859, and D933, and it fits the fingerprint of PI3Kα active inhibitors. Biological evaluation of the reported compounds in HCT116 cell line confirmed that the series inhibited PI3Kα activity and induced apoptosis via activation of caspase-3 and reduction of DNA content.

CONCLUSION

The recently developed compounds might be employed as lead structures for the design of new antitumor drugs targeting PI3Kα.

N-Phenyl-6-Chloro-4-Hydroxy-2-Quinolone-3-CarboxAmides: Molecular Docking, Synthesis, and Biological Investigation as Anticancer Agents

Cancer is a multifactorial disease and the second leading cause of death worldwide. Diverse factors induce carcinogenesis, such as diet, smoking, radiation, and genetic defects. The phosphatidylinositol 3-kinase (PI3Kα) has emerged as an attractive target for anticancer drug design. Eighteen derivatives of N-phenyl-6-chloro-4-hydroxy-2-quinolone-3-carboxamide were synthesized and characterized using FT-IR, NMR (¹H and ¹³C), and high-resolution mass spectra (HRMS). The series exhibited distinct antiproliferative activity (IC₅₀ µM) against human epithelial colorectal adenocarcinoma (Caco-2) and colon carcinoma (HCT-116) cell lines, respectively: compounds 16 (37.4, 8.9 µM), 18 (50.9, 3.3 µM), 19 (17.0, 5.3 µM), and 21 (18.9, 4.9 µM). The induced-fit docking (IFD) studies against PI3Kαs showed that the derivatives occupy the PI3Kα binding site and engage with key binding residues.

Molecular Modeling, Synthesis and Biological Evaluation of N-Phenyl-4-Hydroxy-6-Methyl-2-Quinolone-3-CarboxAmides as Anticancer Agents

The emergence of phosphatidylinositol 3-kinase (PI3Kα) in cancer development has accentuated its significance as a potential target for anticancer drug design. Twenty one derivatives of N-phenyl-4-hydroxy-6-methyl-2-quinolone-3-carboxamide were synthesized and characterized using NMR (¹H and ¹³C) and HRMS. The derivatives displayed inhibitory activity against human epithelial colorectal adenocarcinoma (Caco-2) and human colon cancer (HCT-116) cell lines: compounds 8 (IC₅₀ Caco-2 = 98 µM, IC₅₀ HCT-116 = 337 µM) and 16 (IC₅₀ Caco-2 = 13 µM, IC₅₀ HCT-116 = 240.2 µM). Results showed that compound 16 significantly affected the gene encoding AKT, BAD, and PI3K. The induced-fit docking (IFD) studies against PI3Kα demonstrated that the scaffold accommodates the kinase domains and forms H-bonds with significant binding residues.

Discovery of a novel phosphoinositide 3-kinase gamma (PI3Kγ) inhibitor against hematologic malignancies and theoretical studies on its PI3Kγ-specific binding mechanisms

Class IB phosphoinositide 3-kinase gamma (PI3Kγ) is vital for regulating intracellular signaling pathways and has become an attractive drug target for the treatment of malignant tumors. In the present study, one potent PI3Kγ inhibitor (JN-PK1) with a novel scaffold against hematologic malignancies was identified based on a series of biological experiments, and then the selective mechanism of PI3Kγ inhibition was explored by a systematic computational method. JN-PK1 shows an effective antiproliferative activity on several cancer cell lines, especially blood cancer cells. Cell-free enzymatic studies demonstrated that JN-PK1 specifically inhibits PI3Kγ at low micromolar concentrations without affecting other isoforms of PI3K. In the cellular context, JN-PK1 potently inhibits PI3K/Akt/mTOR signaling pathway in a time- and concentration-dependent manner, which leads to the apoptosis of cancer cells. Further, the specific binding mode of JN-PK1 with PI3Kγ was illustrated by molecular docking, and the selective inhibition mechanism of PI3Kγ by JN-PK1 was revealed by molecular dynamics simulation. Finally, some key residues of PI3Kγ required for specificity and activity were identified. Taken together, JN-PK1 may be developed as a promising therapeutic agent for the treatment of hematologic malignancies.

Ligand-Based Drug Design: Synthesis and Biological Evaluation of Substituted Benzoin Derivatives as Potential Antitumor Agents

Background:

Phosphoinositide 3-kinase α (PI3Kα) has emerged as a promising target for anticancer drug design.

Objectives:

Target compounds were designed to investigate the effect of the p-OCH3 motifs on ligand/PI3Kα complex interaction and antiproliferative activity.

Methods:

Synthesis of the proposed compounds, biological examination tests against human colon adenocarcinoma (HCT-116), breast adenocarcinoma (MCF-7), and breast carcinoma (T47D) cell lines, along with Glide docking studies.

Results:

A series of 1,2-bis(4-methoxyphenyl)-2-oxoethyl benzoates was synthesized and characterized by means of FT-IR, 1H and 13C NMR, and by elemental analysis. Biological investigation demonstrated that the newly synthesized compounds exhibit antiproliferative activity in human colon adenocarcinoma (HCT-116), breast adenocarcinoma (MCF-7), and breast carcinoma (T47D) cell lines possibly via inhibition of PI3Kα and estrogen receptor alpha (ERα). Additionally, results revealed that these compounds exert selective inhibitory activity, induce apoptosis, and suppress VEGF production. Compound 3c exhibited promising antiproliferative activity in HCT-116 interrogating that hydrogen bond-acceptor mediates ligand/PI3Kα complex formation on m- position. Compounds 3e and 3i displayed high inhibitory activity in MCF-7 and T47D implying a wide cleft discloses the o-attachment. Furthermore, compound 3g exerted selective inhibitory activity against T47D. Glide docking studies against PI3Kα and ERα demonstrated that the series accommodate binding to PI3Kα and/or ERα.

Conclusion:

The series exhibited a potential antitumor activity in human carcinoma cell lines encoding PI3Kα and/or ERα.

Cholesterol-coated gold nanorods as an efficient nano-carrier for chemotherapeutic delivery and potential treatment of breast cancer: in vitro studies using the MCF-7 cell line

Gold nanorods (GNRs) have a recognized role in treatment of cancers as efficient nanocarriers for chemotherapeutic drug delivery. In this study, GNRs modified with cholesterol-PEG were employed as a nanocarrier for a hydrophobic compound having a promising phosphatidylinositol 3-kinase (PI3Kα) inhibitory activity. The acquired nanocomplex was characterized by optical and infra-red (IR) absorption spectroscopies, in addition to hydrodynamic size and zeta potential. Glide docking and superposing of docked poses of the hydrophobic ligand and cholesterol moiety demonstrated that hydrophobic interactions drive the conjugation and attachment of the ligand to the cholesterol moiety of the nanocarrier. In vitro release study under a cellular environment indicates that the presence of cells has enhanced the release and the cellular uptake of the conjugated ligand. Furthermore, the anti-proliferative assay of the nanocomplex revealed potent cytotoxicity over a low concentration range of the nanocomplex against MCF-7 breast cancer cells compared to the free compound or the nanocarrier alone. Analysis of cellular death modality by flow cytometry showed that the nanocomplex has a rapid effect on cell death, as cells went toward the late apoptotic/necrotic stage rapidly and proportionally to the increase of the nanocomplex concentration. The overall results propose that cholesterol-decorated GNRs could be considered as a promising nanocarrier for hydrophobic drugs to achieve efficient delivery and potential therapy against breast cancer cells.

Gold nanorods demonstrate a recognized role in the treatment of breast cancer cell lines as an efficient nanocarrier for chemotherapeutic drug delivery.

Novel amphiphilic pyridinium ionic liquids-supported Schiff bases: ultrasound assisted synthesis, molecular docking and anticancer evaluation

Background

Pyridinium Schiff bases and ionic liquids have attracted increasing interest in medicinal chemistry.

Results

A library of 32 cationic fluorinated pyridinium hydrazone-based amphiphiles tethering fluorinated counteranions was synthesized by alkylation of 4-fluoropyridine hydrazone with various long alkyl iodide exploiting lead quaternization and metathesis strategies. All compounds were assessed for their anticancer inhibition activity towards different cancer cell lines and the results revealed that increasing the length of the hydrophobic chain of the synthesized analogues appears to significantly enhance their anticancer activities. Substantial increase in caspase-3 activity was demonstrated upon treatment with the most potent compounds, namely 8, 28, 29 and 32 suggesting an apoptotic cellular death pathway.

Conclusions

Quantum-polarized ligand docking studies against phosphoinositide 3-kinase α displayed that compounds 2–6 bind to the kinase site and form H-bond with S774, K802, H917 and D933.

Electronic supplementary material

The online version of this article (10.1186/s13065-018-0489-z) contains supplementary material, which is available to authorized users.

Molecular Dynamics Simulations to Investigate the Binding Mode of the Natural Product Liphagal with Phosphoinositide 3-Kinase α

Phosphatidylinositol 3-kinase α (PI3Kα) is an attractive target for anticancer drug design. Liphagal, isolated from the marine sponge Aka coralliphaga, possesses the special “liphagane” meroterpenoid carbon skeleton and has been demonstrated as a PI3Kα inhibitor. Molecular docking and molecular dynamics simulations were performed to explore the dynamic behaviors of PI3Kα binding with liphagal, and free energy calculations and energy decomposition analysis were carried out by use of molecular mechanics/Poisson-Boltzmann (generalized Born) surface area (MM/PB(GB)SA) methods. The results reveal that the heteroatom rich aromatic D-ring of liphagal extends towards the polar region of the binding site, and the D-ring 15-hydroxyl and 16-hydroxyl form three hydrogen bonds with Asp810 and Tyr836. The cyclohexyl A-ring projects up into the upper pocket of the lipophilic region, and the hydrophobic/van der Waals interactions with the residues Met772, Trp780, Ile800, Ile848, Val850, Met922, Phe930, Ile932 could be the key interactions for the affinity of liphagal to PI3Kα. Thus, a new strategy for the rational design of more potent analogs of liphagal against PI3Kα is provided. Our proposed PI3Kα/liphagal binding mode would be beneficial for the discovery of new active analogs of liphagal against PI3Kα.

An Algorithm to Identify Target-Selective Ligands - A Case Study of 5-HT7/5-HT1A Receptor Selectivity

A computational procedure to search for selective ligands for structurally related protein targets was developed and verified for serotonergic 5-HT7/5-HT1A receptor ligands. Starting from a set of compounds with annotated activity at both targets (grouped into four classes according to their activity: selective toward each target, not-selective and not-selective but active) and with an additional set of decoys (prepared using DUD methodology), the SVM (Support Vector Machines) models were constructed using a selective subset as positive examples and four remaining classes as negative training examples. Based on these four component models, the consensus classifier was then constructed using a data fusion approach. The combination of two approaches of data representation (molecular fingerprints vs. structural interaction fingerprints), different training set sizes and selection of the best SVM component models for consensus model generation, were evaluated to determine the optimal settings for the developed algorithm. The results showed that consensus models with molecular fingerprints, a larger training set and the selection of component models based on MCC maximization provided the best predictive performance.

Computer-aided design, synthesis, and biological evaluation of new indole-2-carboxamide derivatives as PI3Kα/EGFR inhibitors

Structure-based drug design and molecular modeling were employed to identify a new series of indole-2-carboxamides as potential anticancer agents. These compounds were synthesized and characterized with the aid of several spectroscopic techniques, such as FT-IR, NMR, and mass spectrometry as well as by elemental analysis. Molecular docking studies confirmed that the newly synthesized compounds accommodate PI3Kα and EGFR kinase catalytic sites and form H-bonding with the key binding residues. The antitumor activity of these new compounds against an array of cancer cell lines (human colon carcinoma (HCT116), leukemia (K562), and breast cancer (MDA231) was evaluated. Results revealed that these compounds were selective against the kinase domain, and none of them showed any inhibitory activity against K562. In addition, results showed that compound 13 exhibited high potency in HCT116 and MDA231 with IC50 values of 19 and 15μM, respectively. Our findings recommend that further optimization of this series would be beneficial for colon and breast cancer treatment.

Molecular modeling of exquisitely selective c-Met inhibitors through 3D-QSAR and molecular dynamics simulations

c-Met has been considered as an attractive target for developing antitumor agents. The highly selective c-Met inhibitors provide invaluable opportunities for the combination with other therapies safely to achieve the optimal efficacy. In this work, a series of triazolopyrazine c-Met inhibitors with exquisitely selectivity were investigated using a combination of molecular docking, three-dimensional quantitative structure-activity relationship (3D-QSAR), and molecular dynamics simulation. Comparative molecular field analysis (CoMFA) and comparative molecular similarity index analysis (CoMSIA) models were developed to reveal the structural determinants for c-Met inhibition. Both models were validated to have high reliability and predictability, and contour map analysis suggested feature requirements for different substituents on the scaffold. It is worth noting that an important hydrogen bond rich region was identified in the unique narrow channel, which is distinct from other kinases. Molecular dynamics simulations and binding free energy calculations provided further support that suitable groups in this hydrogen bond rich region made great contributions to the binding of ligands. Moreover, hydrogen bonds with residues of the narrow channel were also indicated to be essential to improve the activity and selectivity. This study will facilitate the discovery and optimization of novel c-Met inhibitors with higher activity and selectivity.

Discovery of nanomolar phosphoinositide 3-kinase gamma (PI3Kγ) inhibitors using ligand-based modeling and virtual screening followed by in vitro analysis

Phosphoinositide 3-kinase gamma (PI3Kγ) is member of a family of enzymes involved in cancer pathogenesis. Accordingly, considerable efforts have been carried out to develop new PI3Kγ inhibitors. Towards this end we explored the pharmacophoric space of PI3Kγ using three diverse sets of inhibitors. Subsequently, we employed genetic algorithm-based QSAR analysis to select optimal combination of pharmacophoric models and physicochemical descriptors that can explain bioactivity variation within training inhibitors. Interestingly, two successful pharmacophores were selected within two statistically consistent QSAR models. The close similarity among the two binding models prompted us to merge them in a hybrid pharmacophore. The resulting model showed superior receiver operator characteristic curve (ROC) and closely resembled binding interactions seen in crystallographic ligand-PI3Kγ complexes. The resulting model was employed to screen the national cancer institute (NCI) list of compounds to search for new PI3Kγ ligands. After testing captured hits in vitro, 19 compounds showed nanomolar IC50 values against PI3Kγ. The chemical structures and purities of most potent hits were validated using NMR and MS experiments.

Theoretical studies on beta and delta isoform-specific binding mechanisms of phosphoinositide 3-kinase inhibitors

Phosphoinositide 3-kinase (PI3K) is known to be closely related to tumorigenesis and cell proliferation, and controls a variety of cellular processes, including proliferation, growth, apoptosis, migration, metabolism, etc. The PI3K family comprises eight catalytic isoforms, which are subdivided into three classes. Recently, the discovery of inhibitors that block a single isoform of PI3K has continued to attract special attention because they may have higher selectivity for certain tumors and less toxicity for healthy cells. The PI3Kβ and PI3Kδ share fewer studies than α/γ, and therefore, in this work, the combination of molecular dynamics simulations and free energy calculations was employed to explore the binding of three isoform-specific PI3K inhibitors (COM8, IC87114, and GDC-0941) to PI3Kβ or PI3Kδ. The isoform specificities of the studied inhibitors derived from the predicted binding free energies are in good agreement with the experimental data. In addition, the key residues critical for PI3Kβ or PI3Kδ selectivity were highlighted by decomposing the binding free energies into the contributions from individual residues. It was observed that although PI3Kβ and PI3Kδ share the conserved ATP-binding pockets, individual residues do behave differently, particularly the residues critical for PI3Kβ or PI3Kδ selectivity. It can be concluded that the inhibitor specificity between PI3Kβ and PI3Kδ is determined by the additive contributions from multiple residues, not just a single one. This study provides valuable information for understanding the isoform-specific binding mechanisms of PI3K inhibitors, and should be useful for the rational design of novel and selective PI3K inhibitors.