Class I Phospho-inositide-3-kinases (PI3Ks) Isoform-Specific Inhibition Study by the Combination of Docking and Molecular Dynamics Simulation

Isobavachalcone Activates Antitumor Immunity on Orthotopic Pancreatic Cancer Model: A Screening and Validation

Pancreatic cancer is accompanied by poor prognosis and accounts for a significant number of deaths every year. Since Psoralea corylifolia L. (PCL) possesses a broad spectrum of bioactivities, it is commonly used in traditional Chinese medicine. The study explored potential antitumor agents of PCL and underlying mechanisms in vitro and vivo. Based on network pharmacology, bioinformatics, and molecular docking, we considered isobavachalcone (IBC) as a valuable compound. The activity and potential mechanisms of IBC were investigated by RT-qPCR, immunohistochemistry, immunofluorescence, and flow cytometry. It was confirmed that IBC could inhibit Panc 02 cell proliferation and induce apoptosis via increasing the production of reactive oxygen species. IBC could attenuate the weight of solid tumors, increase CD8⁺ T cells, and reduce M2 macrophages in the tumor tissue and spleen. Another promising finding was that IBC alleviated the proportion of myeloid-derived suppressor cells (MDSCs) in the tumor tissue but had no change in the spleen. The study of pharmacological effects of IBC was carried out and suggested IBC restrained M2-like polarization of RAW 264.7 cells by inhibiting the expression of ARG1 and MRC1 and suppressed the expression of ARG1 and TGF-β in bone marrow-derived MDSC. In summary, this research screened IBC as an antineoplastic agent, which could attenuate the growth of pancreatic cancer via activating the immune activity and inducing cell apoptosis. It might be a reference for the antitumor ability of IBC and the treatment of the tumor microenvironment in pancreatic cancer.

Mechanism of drug resistance of BVDV induced by F224S mutation in RdRp: A case study of VP32947

Bovine viral diarrhea virus (BVDV) is an enveloped virus with an RNA genome, causing serious economic losses to the areas dominated by livestock industry. Currently, although several compounds with biological activities of inhibiting virus replication have been reported, amino acid mutations (especially F224S mutation) frequently occurring in the RNA-dependent RNA polymerase (RdRp) have greatly reduce their value of further research. In this study, we introduced an effective and rapid in silico strategy to explore the differences in the binding modes of VP32947 between the wild/mutant-type RdRp at the molecular level, and further explained the main reasons for the variations in the inhibitory activities of VP32947 against the two types of enzymes. Firstly, the binding site of VP32947 in the finger domain was determined based on the previously reported experimental data, and the initial conformation of VP32947 in the wild RdRp was constructed using molecular docking. Then, the mutant research system was obtained directly by artificial mutation strategy. Afterwards, the built research systems were subjected to microsecond-timescale molecular dynamic simulation, and the conformational and energic profile analyses were performed according to the simulation trajectories. It was found that after 1 μs simulation, VP32947 in the mutant system was transferred to the left side of Loop α, and its interactions with the residues in the loop region were weakened. However, VP32947 in the wild system remained at the right side of Loop α, and could have a good fit with the sub-pocket formed by F224, I261, P262, N264, S532, which was conducive to maintaining its stable binding conformation in the wild RdRp. The illustration of the difference in the binding mechanisms of VP32947 in the wild/mutant RdRp would provide a theoretical basis for the rational design of innovative inhibitors based on the enzyme.

Probing Allosteric Regulation Mechanism of W7.35 on Agonist-Induced Activity for μOR by Mutation Simulation

The residue located at 15 positions before the most conserved residue in TM7 (7.35 of Ballesteros-Weinstein number) plays important roles in ligand binding and the receptor activity for class A GPCRs. Nevertheless, its regulation mechanism has not been clearly clarified in experiments, and some controversies also exist for its impact on μ-opioid receptors (μOR) bound by agonists. Thus, we chose the μ-opioid receptor (μOR) of class A GPCRs as a representative and utilized a microsecond accelerated molecular dynamics simulation (aMD) coupled with a protein structure network (PSN) to explore the effect of W318^7.35 on its functional activity induced by the agonist endomorphin2 mainly by a comparison of the wild system and its W7.35A mutant. When endomorphin2 binds to the wild-type μOR, TM6 in μOR moves outward to form an open intracellular conformation that is beneficial to accommodating the β-arrestin transducer, rather than the G-protein transducer due to the clash with the α5 helix of G-protein, thus acting as a β-arrestin biased agonist. However, the W318A mutation induces the intracellular part of μOR to form a closed state, which disfavors coupling with either G-protein or β-arrestin. The allosteric pathway analysis further reveals that the binding of endomorphin2 to the wild-type μOR transmits more activation signals to the β-arrestin binding site while the W318A mutation induces more structural signals to transmit to common binding residues of the G protein and β-arrestin. More interestingly, the residue at the 7.35 position regulates the shortest allosteric pathway in indirect ways by influencing the interactions between other ligand-binding residues and endomorphin2. W293^6.48 and F289^6.44 are important for regulating the different activities of μOR induced either by the agonist or by the mutation. Y336^7.53, F343^8.50, and D340^8.47 play crucial roles in activating the β-arrestin biased signal induced by the agonist endomorphin2, while L158^3.43 and V286^6.41 devote important contributions to the change in the activity of endomorphin2 from the β-arrestin biased agonist to the antagonist upon the W318A mutation.

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.

The Chinese herbal prescription JZ-1 induces autophagy to protect against herpes simplex Virus-2 in human vaginal epithelial cells by inhibiting the PI3K/Akt/mTOR pathway

ETHNOPHAMACOLOGICAL RELEVANCE

The Chinese herbal prescription JieZe-1 (JZ-1) is based on the modification of Yihuang Tang, which was first described in Fu Qingzhu Nvke by the famous Qing Dynasty doctor Shan Fu as a treatment for leukorrheal diseases. As an in-hospital preparation, JZ-1 has been used in Tongji Hospital for many years to treat various infectious diseases of the lower female genital tract, including cervicitis, vaginitis, genital herpes and condyloma acuminatum. Our previous studies have shown that JZ-1 has curative effects on Candida albicans, Trichomonas vaginalis and Ureaplasma urealyticum infections.

AIM OF THE STUDY

Genital herpes is among the most common sexually transmitted diseases (STDs) worldwide and is mainly caused by herpes simplex virus type-2 (HSV-2). Current therapies can relieve symptoms in patients but do not cure or prevent the spread of the virus. This study was designed to investigate the effect of JZ-1 on HSV-2 infection and its mechanism, which is based on autophagy induction, to provide new ideas and a basis for the study of antiviral drugs.

MATERIALS AND METHODS

Evaluation of the antiviral activity of JZ-1 was conducted by MTT assay and western blotting. Then, Western blot and immunofluorescence analyses, observations through transmission electron microscopy and experiments with the recombinant lentivirus vector mRFP-GFP-LC3B were used to monitor autophagic flux in VK2/E6E7 cells. To explore the mechanism by which JZ-1 regulates autophagy, western blotting and real-time quantitative PCR (qRT-PCR) were used to determine the expression of phosphoinositide 3-kinase (PI3K)/Akt/mTOR pathway proteins and to detect changes in critical molecules in the pathway after the application of a PI3K inhibitor. Additionally, the mRNA expression levels of inflammatory cytokines, namely, IL-6, IFN-α, IFN-β and TNF-α, were measured with qRT-PCR.

RESULTS

HSV-2 infection inhibited autophagy in the VK2/E6E7 cells. Further study revealed that the activation of the PI3K/Akt/mTOR pathway induced by HSV-2 infection may result in the blocked autophagic flux and inhibited autophagosome and autolysosome formation. JZ-1 exhibited significant antiviral activity in the VK2/E6E7 cells, which showed increased cell vitality and reduced viral protein expression, namely, earliest virus-specific infected cell polypeptides 5 (ICP5) and glycoprotein D (gD). We found that JZ-1 treatment inhibited the upregulation of the PI3K/Akt/mTOR pathway proteins and promoted autophagy to combat HSV-2 infection, while PI3K inhibitor pretreatment prevented the enhanced autophagy induced by JZ-1. Moreover, JZ-1 attenuated the increase in inflammatory cytokines that had been induced HSV-2 infection.

CONCLUSION

Our results showed that JZ-1 protects against HSV-2 infection, and this beneficial effect may be mediated by inducing autophagy via inhibition of the PI3K/Akt/mTOR signaling axis.

Probing the Druggablility on the Interface of the Protein-Protein Interaction and Its Allosteric Regulation Mechanism on the Drug Screening for the CXCR4 Homodimer

Modulating protein–protein interactions (PPIs) with small drug-like molecules targeting it exhibits great promise in modern drug discovery. G protein-coupled receptors (GPCRs) are the largest family of targeted proteins and could form dimers in living biological cells through PPIs. However, compared to drug development of the orthosteric site, there has been lack of investigations on the druggability of the PPI interface for GPCRs and its functional implication on experiments. Thus, in order to address these issues, we constructed a novel computational strategy, which involved in molecular dynamics simulation, virtual screening and protein structure network (PSN), to study one representative GPCR homodimer (CXCR4). One druggable pocket was identified in the PPI interface and one small molecule targeting it was screened, which could strengthen PPI mainly through hydrophobic interaction between the benzene rings of the PPI molecule and TM4 of the receptor. The PSN results further reveals that the PPI molecule could increase the number of the allosteric regulation pathways between the druggable pocket of the dimer interface to the orthostatic site for the subunit A but only play minor role for the other subunit B, leading to the asymmetric change in the volume of the binding pockets for the two subunits (increase for the subunit A and minor change for the subunit B). Consequently, the screening performance of the subunit A to the antagonists is enhanced while the subunit B is unchanged nearly, implying that the PPI molecule may be beneficial to enhance the drug efficacies of the antagonists. In addition, one main regulation pathway with the highest frequency was identified for the subunit A, which consists of Trp195^5.34–Tyr190^ECL2–Val196^5.35–Gln200^5.39–Asp262^6.58–Cys28^N-term, revealing their importance in the allosteric regulation from the PPI molecule. The observations from the work could provide valuable information for the development of the PPI drug-like molecule for GPCRs.

Joint action of miR-126 and MAPK/PI3K inhibitors against metastatic melanoma

Emerging data support the rationale of combined therapies in advanced melanoma. Specifically, the combined use of drugs with different mechanisms of action can reduce the probability of selecting resistant clones. To identify agents active against melanoma cells, we screened a library of 349 anti‐cancer compounds, currently in clinical use or trials, and selected PIK‐75, an inhibitor of the phosphatidylinositol 3‐kinase/protein kinase B (PI3K/AKT) pathway, as the ‘top active’ drug. PIK‐75 was then used alone or in combination with vemurafenib, the first BRAF inhibitor approved for patients with melanoma harboring BRAF mutations. We identified a combined dose of PIK‐75 and vemurafenib that inhibited both the PI3K/AKT and mitogen‐activated protein kinase pathways, thereby overcoming any compensatory activation. In view of the important tumor suppressor function induced by restoring expression of microRNA (miR)‐126 in metastatic melanoma cells, we examined whether miR‐126 has a synergistic role when included in a triple combination alongside PIK‐75 and vemurafenib. We found that enforced expression of miR‐126 (which alone can reduce tumorigenicity) significantly increased PIK‐75 activity when used as either a single agent or in combination with vemurafenib. Interestingly, PIK‐75 proved to be effective against early passage cell lines derived from patients’ biopsies and on melanoma cell lines resistant to either vemurafenib or dabrafenib, thus suggesting that it potentially has the capability to overcome drug resistance. Finally, the synergistic role played by miR‐126 in combination with vemurafenib and/or PIK‐75 was demonstrated in vivo in mouse xenograft models, in which tumor growth inhibition was associated with increased apoptosis. These results not only show the efficacy of PIK‐75 and vemurafenib co‐treatment but also indicate that restoration of miR‐126 expression in advanced melanoma can enhance their antitumor activity, which may possibly allow dose reduction to decrease adverse events without reducing the therapeutic benefits.

Structural Determinants of Isoform Selectivity in PI3K Inhibitors

Phosphatidylinositol 3-kinases (PI3Ks) are important therapeutic targets for the treatment of cancer, thrombosis, and inflammatory and immune diseases. The four highly homologous Class I isoforms, PI3K, PI3K, PI3K and PI3K have unique, non-redundant physiological roles and as such, isoform selectivity has been a key consideration driving inhibitor design and development. In this review, we discuss the structural biology of PI3Ks and how our growing knowledge of structure has influenced the medicinal chemistry of PI3K inhibitors. We present an analysis of the available structure-selectivity-activity relationship data to highlight key insights into how the various regions of the PI3K binding site influence isoform selectivity. The picture that emerges is one that is far from simple and emphasizes the complex nature of protein-inhibitor binding, involving protein flexibility, energetics, water networks and interactions with non-conserved residues.

Structural Features and Ligand Selectivity for 10 Intermediates in the Activation Process of β2-Adrenergic Receptor

It has already been suggested by researchers that there should be multiple intermediate states in the activation process for G-protein-coupled receptors (GPCRs). However, the intermediate states are very short-lived and hardly captured by the experiments, leading to very limited understanding of their structural features and drug efficacies. In this work, a novel joint strategy of targeted molecular dynamics simulation, conventional molecular dynamics simulation, and virtual screening is developed to address the problems. The results from 10 intermediate conformations obtained from the work reveal that the ligand pocket is very unstable and fluctuates between the inactive state and the active one in the case of ligand-free, in particular for ECL2 as a gate-keeper of the ligand-binding. The ligand-binding site could be stable in the active state with a small volume and a completely closed ECL2, only when the G-protein-binding region is fully activated. In addition, the activations of the ligand-binding pocket and G-protein-binding site are relatively independent and exhibit a loose allosteric coupling, which contributes to the existence of multiple intermediate conformations. Interestingly, the screening performance of the agonists does not increase on increasing the overall activity of the intermediate state, but is dependent on the activated extent of the ligand pocket. The receptor is prone to bind the agonist when closing ECL2 and reducing the ligand-binding pocket volume, whereas it is more favorable for binding the antagonist when opening ECL2 and increasing the pocket volume. These observations added to previous studies could help us better understand the activation mechanism of GPCRs and provide valuable information for drug design.

Definition of the binding mode of phosphoinositide 3-kinase α-selective inhibitor A-66S through molecular dynamics simulation

Activation of the phosphatidylinositol 3-kinase α (PI3Kα) is commonly observed in human cancer and is critical for tumor progression, which has made PI3Kα an attractive target for anticancer drug discovery. To systematically investigate the binding mode of A-66S, a new selective PI3Kα inhibitor for PI3Kα, molecular docking, molecular dynamics simulation and ensuing energetic analysis were performed. The binding free energy between PI3Kα and A-66S is -11.27 kcal•mol⁻¹ using MMPBSA method, while -14.67 kcal•mol⁻¹ using MMGBSA method, which is beneficial for the binding, and the van der Waals/hydrophobic and electrostatic interactions are critical for the binding. The conserved hydrophobic adenine region of PI3Kα made up of Met772, Pro778, Ile800, Tyr836, Ile848, Val850, Val851, Met922, Phe930 and Ile932 accommodates the flat 2-tert-butyl-4'-methyl-4,5'-bithiazol moiety of A-66S, and the NH of Val851 forms a hydrogen with the nitrogen atom embedded in the aminothiazole ring of A-66S. The (S)-pyrrolidine carboxamide urea moiety especially extends toward the region of the binding site wall (Ser854-Gln859) defined by the C-terminal lobe, and has three hydrogen-bond arms with the backbone of Ser854 and the side chain of Gln859. Notably the interaction between the non-conserved residue Gln859 and A-66S is responsible for the selectivity profile of A-66S. The binding mode of A-66S for PI3Kα presented in this study should aid in the design of a new highly selective PI3Kα inhibitor.

Structure, function and inhibition of the phosphoinositide 3-kinase p110α enzyme

The PI3K (phosphoinositide 3-kinase) p110α isoform is activated by oncogenic mutations in many cancers. This has stimulated intense interest in identifying inhibitors of the PI3K pathway as well as p110α-selective inhibitors, and understanding the mechanisms underlying activation by the oncogenic mutations. In the present article, we review recent progress in the structure and function of the p110α enzyme and two of its most common oncogenic mutations, the development of isoform-selective inhibitors, and p110α pharmacology.

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.

Synthesis and cancer stem cell-based activity of substituted 5-morpholino-7H-thieno[3,2-b]pyran-7-ones designed as next generation PI3K inhibitors

Dysregulation of the phosphatidylinositol-3-kinase (PI3K) pathway in a wide range of tumors has made PI3K a consensus target to inhibit as illustrated by more than 15 inhibitors now in clinical trials. Our previous work, built on the early pioneering multikinase inhibitor LY294002, resulted in the only PI3K vascular-targeted PI3K inhibitor prodrug, SF1126, which has now completed Phase I clinical trials. This inhibitor has properties that impart more in vivo activity than should be warranted by its enzymatic potency, which in general is much lower than other clinical stage PI3K inhibitors. We embarked on the exploration of scaffolds that retained such properties while simultaneously exhibiting an increased potency toward PI3K. This work resulted in the discovery of the 5-morpholino-7H-thieno[3,2-b]pyran-7-one system as the foundation of a new compound class of potential PI3K inhibitors having improved potency toward PI3K. The synthesis and cancer stem cell-based activity of these compounds are reported herein.

Molecular simulation of a series of benzothiazole PI3Kα inhibitors: probing the relationship between structural features, anti-tumor potency and selectivity

The phosphatidylinositol 3-kinase α (PI3Kα) was genetically validated as a promising therapeutic target for developing novel anticancer drugs. In order to explore the structure-activity correlation of benzothiazole series as inhibitors of PI3Kα, comparative molecular field analysis (CoMFA), comparative molecular similarity indices analysis (CoMSIA) were performed on 61 promising molecules to build 3D-QSAR models based on both the ligand- and receptor-based methods. The best CoMFA and CoMSIA models had a cross-validated coefficient r(cv)(2) of 0.618 and 0.621, predicted correlation coefficient r(pred) (2) of 0.812 and 0.83, respectively, proving their high correlative and predictive abilities on both the training and test sets. In addition, docking analysis and molecular dynamics simulation (MD) were also applied to elucidate the probable binding modes of these inhibitors at the ATP binding pocket. Based on the contour maps and MD results, some key structural factors responsible for the activity of this series of compounds were revealed as follows: (1) Ring-A has a strong preference for bulky hydrophobic or aromatic groups; (2) Electron-withdrawing groups at the para position of ring-B and hydrophilic substituents in ring-B region may benefit the potency; (3) A polar substituent like -NHSO(2)- between ring-A and ring-B can enhance the activity of the drug by providing hydrogen bonding interaction with the protein target. The satisfactory results obtained from this work strongly suggest that the developed 3D-QSAR models and the obtained PI3Kα inhibitor binding structures are reasonable for the prediction of the activity of new inhibitors and be helpful in future PI3Kα inhibitor design.

Characterization of molecular recognition of phosphoinositide-3-kinase α inhibitor through molecular dynamics simulation

Phosphatidylinositol 3-kinase α (PI3Kα) is a promising target for anticancer drug discovery due to its overactivation in tumor cells. To systematically investigate the interactions between PI3Kα and PIK75 which is the most selective PI3Kα inhibitor reported to date, molecular docking, molecular dynamics simulation, and ensuing energetic analysis were utilized. The binding free energy between PI3Kα and PIK75 is -10.04 kcal•mol(-1) using MMPBSA method, while -13.88 kcal•mol(-1) using MMGBSA method, which is beneficial for the binding. The van der Waals/hydrophobic and electrostatic interactions play critical roles for the binding. The binding mode of PIK75 for PI3Kα is predicted. The conserved hydrophobic adenine region of PI3Kα made up of Ile800, Ile848, Val850, Val851, Met922, Phe930, and Ile932 accommodates the flat 6-bromine imidazo[1,2-a]pyridine ring of PIK75. The 2-methyl-5-nitrophenyl group of PIK75 extends to the P-loop region, and has four hydrogen-bond arms with the backbone and side chain of Ser773 and Ser774. And the distinct conformation of the P-loop induced by PIK75 is speculated to be responsible for the selectivity profile of PIK75. The predicted binding mode of PIK75 for PI3Kα presented in this study may help design high affinity and selective compounds to target PI3Kα.

Isoform-selective inhibition of phosphoinositide 3-kinase: identification of a new region of nonconserved amino acids critical for p110α inhibition

The combination of molecular modeling and X-ray crystallography has failed to yield a consensus model of the mechanism for selective binding of inhibitors to the phosphoinositide 3-kinase (PI3K) p110 α-isoform. Here we have used kinetic analysis to determine that the p110α-selective inhibitor 2-methyl-5-nitro-2-[(6-bromoimidazo[1,2-α]pyridin-3-yl)methylene]-1-methylhydrazide-benzenesulfonic acid (PIK-75) is a competitive inhibitor with respect to a substrate, phosphatidylinositol (PI) in contrast to most other PI3K inhibitors, which bind at or near the ATP site. Using sequence analysis and the existing crystal structures of inhibitor complexes with the p110γ and -δ isoforms, we have identified a new region of nonconserved amino acids (region 2) that was postulated to be involved in PIK-75 p110α selectivity. Analysis of region 2, using in vitro mutation of identified nonconserved amino acids to alanine, showed that Ser773 was a critical amino acid involved in PIK-75 binding, with an 8-fold-increase in the IC(50) compared with wild-type. Kinetic analysis showed that, with respect to PI, the PIK-75 K(i) for the isoform mutant S773D increased 64-fold compared with wild-type enzyme. In addition, a nonconserved amino acid, His855, from the previously identified region 1 of nonconserved amino acids, was found to be involved in PIK-75 binding. These results show that these two regions of nonconserved amino acids that are close to the substrate binding site could be targeted to produce p110α isoform-selective inhibitors.

Docking studies on isoform-specific inhibition of phosphoinositide-3-kinases

Phosphatidylinositol 3-kinase α (PI3Kα) is a promising target for anticancer drug design. Oncogenic mutation H1047R in the catalytic domain is observed in many tumors and may enhance PI3Kα kinase activity by affecting loop confirmations as well as membrane binding. We applied docking methods to 33 PI3K inhibitors against the wild type (wt) PI3Kα, the H1047R mutant of PI3Kα and the γ isoform of PI3K (PI3Kγ). We also investigated the effect of protein flexibility on ligand binding by docking the same set of ligands to conformations of the wt and mutant PI3Kα generated by molecular dynamics simulations. Our data suggests that conformational differences in Gln859, Ser854, Tyr836, and Ser774 between the PI3Kα wt and H1047R mutant may be used to design ligands that are active against both the wt and H1047R mutant isoforms. Gln859, Ser854 and Ser774 may play critical roles in ligand binding to the α isoform H1047R mutant while formation of H-bonds with Ser806 of PI3Kγ may enhance γ-isoform-specific inhibition. In addition to H-bond interactions, structural and size differences in the activation and hydrophobic domains of PI3Kα, PI3Kγ, and the PI3Kα H1047R mutant could be exploited to direct the design of isoform- and/or mutant-specific PI3K inhibitors. Our data provide a reasonable explanation for the activity and selectivity of small molecular PI3K inhibitors and are in good agreement with available experimental and computational data.

Thiazolidinedione-based PI3Kα inhibitors: an analysis of biochemical and virtual screening methods

A series of synthesized and commercially available compounds were assessed against PI3Kα for in vitro inhibitory activity and the results compared to binding calculated in silico. Using published crystal structures of PI3Kγ and PI3Kδ co-crystallized with inhibitors as a template, docking was able to identify the majority of potent inhibitors from a decoy set of 1000 compounds. On the other hand, PI3Kα in the apo-form, modeled by induced fit docking, or built as a homology model gave only poor results. A PI3Kα homology model derived from a ligand-bound PI3Kδ crystal structure was developed that has a good ability to identify active compounds. The docking results identified binding poses for active compounds that differ from those identified to date and can contribute to our understanding of structure-activity relationships for PI3K inhibitors.