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

Journal of Molecular Modeling - Tập 18 - Trang 2943-2958 - 2011
Jinan Wang1, Fangfang Wang1, Zhengtao Xiao1, Guowen Sheng2, Yan Li3, Yonghua Wang1,2
1Center of Bioinformatics, Northwest A&F University, Yangling, China
2College of Life Sciences, Northwest A&F University, Yangling, China
3Department of Materials Science and Chemical Engineering, Dalian University of Technology, Dalian, China

Tóm tắt

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 rcv 2 of 0.618 and 0.621, predicted correlation coefficient rpred 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 -NHSO2- 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.

Tài liệu tham khảo

Leevers SJ, Vanhaesebroeck B, Waterfield MD (1999) Signalling through phosphoinositide 3-kinases: the lipids take centre stage. Curr Opin Cell Biol 11:219–225 Stein RC, Waterfield MD (2000) PI3-kinase inhibition: a target for drug development? Mol Med Today 6:347–357 Chang HW, Aoki M, Fruman D, Auger KR, Bellacosa A, Tsichlis PN, Cantley LC, Roberts TM, Vogt PK (1997) Transformation of chicken cells by the gene encoding the catalytic subunit of PI 3-kinase. Science 276:1848–1850 Katso R, Okkenhaug K, Ahmadi K, White S, Timms J, Waterfield MD (2001) Cellular function of phosphoinositide 3-kinases: implications for development, homeostasis, and cancer. Annu Rev Cell Dev Biol 17:615–675 Jiang BH, Liu LZ (2008) PI3K/PTEN signaling in tumorigenesis and angiogenesis. BBA-Proteins Proteom 1784:150–158 Samuels Y, Wang Z, Bardelli A, Silliman N, Ptak J, Szabo S, Yan H, Gazdar A, Powell SM, Riggins GJ, Willson JK, Markowitz S, Kinzler KW, Vogelstein B, Velculescu VE (2004) High frequency of mutations of the PIK3CA gene in human cancers. Science 304:554 Graupera M, Guillermet-Guibert J, Foukas LC, Phng LK, Cain RJ, Salpekar A, Pearce W, Meek S, Millan J, Cutillas PR, Smith AJ, Ridley AJ, Ruhrberg C, Gerhardt H, Vanhaesebroeck B (2008) Angiogenesis selectively requires the p110alpha isoform of PI3K to control endothelial cell migration. Nature 453:662–666 Schultz RM, Merriman RL, Andis SL, Bonjouklian R, Grindey GB, Rutherford PG, Gallegos A, Massey K, Powis G (1995) In vitro and in vivo antitumor activity of the phosphatidylinositol-3-kinase inhibitor, wortmannin. Anticancer Res 15:1135–1139 Wymann MP, Bulgarelli-Leva G, Zvelebil MJ, Pirola L, Vanhaesebroeck B, Waterfield MD, Panayotou G (1996) Wortmannin inactivates phosphoinositide 3-kinase by covalent modification of Lys-802, a residue involved in the phosphate transfer reaction. Mol Cell Biol 16:1722–1733 Vlahos CJ, Matter WF, Hui KY, Brown RF (1994) A specific inhibitor of phosphatidylinositol 3-kinase, 2-(4-morpholinyl)-8-phenyl-4H-1-benzopyran-4-one (LY294002). J Biol Chem 269:5241–5248 Frojdo S, Cozzone D, Vidal H, Pirola L (2007) Resveratrol is a class IA phosphoinositide 3-kinase inhibitor. Biochem J 406:511–518 Marion F, Williams DE, Patrick BO, Hollander I, Mallon R, Kim SC, Roll DM, Feldberg L, Van Soest R, Andersen RJ (2006) Liphagal, a Selective inhibitor of PI3 kinase alpha isolated from the sponge akacoralliphaga: structure elucidation and biomimetic synthesis. Org Lett 8:321–324 Pat. WO2006081659, 2006-08-10 Hayakawa M, Kaizawa H, Moritomo H, Koizumi T, Ohishi T, Okada M, Ohta M, Tsukamoto S, Parker P, Workman P, Waterfield M (2006) Synthesis and biological evaluation of 4-morpholino-2-phenylquinazolines and related derivatives as novel PI3 kinase p110alpha inhibitors. Bioorg Med Chem 14:6847–6858 Hayakawa M, Kawaguchi K, Kaizawa H, Koizumi T, Ohishi T, Yamano M, Okada M, Ohta M, Tsukamoto S, Raynaud FI, Parker P, Workman P, Waterfield MD (2007) Synthesis and biological evaluation of sulfonylhydrazone-substituted imidazo[1,2-a]pyridines as novel PI3 kinase p110alpha inhibitors. Bioorg Med Chem 15:5837–5844 Frederick R, Denny WA (2008) Phosphoinositide-3-kinases (PI3Ks): combined comparative modeling and 3D-QSAR to rationalize the inhibition of p110alpha. J Chem Inf Model 48:629–638 Han M, Zhang JZ (2010) Class I phospho-inositide-3-kinases (PI3Ks) isoform-specific inhibition study by the combination of docking and molecular dynamics simulation. J Chem Inf Model 50:136–145 Li Y, Wang Y, Zhang F (2010) Pharmacophore modeling and 3D-QSAR analysis of phosphoinositide 3-kinase p110alpha inhibitors. J Mol Model 16:1449–1460 D'Angelo ND, Kim TS, Andrews K, Booker SK, Caenepeel S, Chen K, D'Amico D, Freeman D, Jiang JA, Liu LB, McCarter JD, Miguel TS, Mullady EL, Schrag M, Subramanian R, Tang J, Wahl RC, Wang L, Whittington DA, Wu TA, Xi N, Xu Y, Yakowec P, Yang K, Zalameda LP, Zhang N, Hughes P, Norman MH (2011) Discovery and Optimization of a Series of Benzothiazole Phosphoinositide 3-Kinase (PI3K)/Mammalian Target of Rapamycin (mTOR) Dual Inhibitors. J Med Chem 54:1789–1811 Cramer RD III, Bunce JD, Patterson DE, Frank IE (1988) Crossvalidation, bootstrapping, and partial least squares compared with multiple regression in conventional QSAR studies. Quant Struct Act Rel 7:18–25 Klebe G, Abraham U, Mietzner T (1994) Molecular similarity indices in a comparative analysis (CoMSIA) of drug molecules to correlate and predict their biological activity. J Med Chem 37:4130–4146 Matthew C, Richard DC III, Van Nicole O (1989) J Comput Chem 10:982–1012 Huang CH, Mandelker D, Schmidt-Kittler O, Samuels Y, Velculescu VE, Kinzler KW, Vogelstein B, Gabelli SB, Amzel LM (2007) The structure of a human p110α/p85α complex elucidates the effects of oncogenic PI3Kα mutations. Science 318:1744–1748 Mandelker D, Gabelli SB, Schmidt-Kittler O, Zhu J, Cheong I, Huang CH, Kinzler KW, Vogelstein B, Amzel LM (2009) A frequent kinase domain mutation that changes the interaction between PI3Kα and the membrane. Proc Natl Acad Sci USA 106:16996–17001 Holm L, Park J (2000) DaliLite workbench for protein structure comparison. Bioinformatics 16:566–567 Wang Y, Li Y, Ma Z, Yang W, Ai C (2010) Mechanism of microrna-target interaction: Molecular dynamics simulations and thermodynamics analysis. PLoS Comput Biol 6:e1000866 Kuntz I, Blaney J, Oatley S, Langridge R, Ferrin T (1982) A geometric approach to macromolecule-ligand recognitions. J Mol Biol 161:269–288 Jones G, Willett P, Glen RC, Leach AR, Taylor R (1997) Development and validation of a genetic algorithm for flexible docking1. J Mol Biol 267:727–748 Eldridge M, Murray C, Auton T, Paolini G, Mee R (1997) The development of a fast empirical scoring function to estimate the binding affinity of ligands in receptor complexes. J Comput Aided Mol Des 11:425–445 Muegge I, Martin YC (1999) A general and fast scoring function for protein-ligand interactions: a simplified potential approach. J Med Chem 42:791–804 AbdulHameed MDM, Hamza A, Liu J, Zhan CG (2008) Combined 3D-QSAR modeling and molecular docking study on indolinone derivatives as inhibitors of 3-phosphoinositide-dependent protein kinase-1. J Chem Inf Model 48:1760–1772 Cramer RD III, Patterson DE, Bunce JD (1988) Comparative molecular field analysis (CoMFA). 1. Effect of shape on binding of steroids to carrier proteins. J Am Chem Soc 110:5959–5967 Nilsson J (1998) Multiway Calibration in 3D QSAR. University of Groningen, Applications to Dopamine Receptor Ligands Uddin R, Yuan H, Petukhov PA, Choudhary MI, Madura JD (2008) Receptor-based modeling and 3D-QSAR for a quantitative production of the butyrylcholinesterase inhibitors based on genetic algorithm. J Chem Inf Model 48:1092–1103 Wold S, Ruhe A, Wold H, Dunn W III (1984) The collinearity problem in linear regression. The partial least squares (PLS) approach to generalized inverses. SIAM J Sci Stat Comput 5:735–743 Case D, Darden T, Cheatham Iii T, Simmerling C, Wang J, Duke R, Luo R, Crowley M, Walker RC, Zhang W (2008) AMBER 10. University of California, San Francisco Duan Y, Wu C, Chowdhury S, Lee MC, Xiong G, Zhang W, Yang R, Cieplak P, Luo R, Lee T (2003) A point-charge force field for molecular mechanics simulations of proteins based on condensed-phase quantum mechanical calculations. J Comput Chem 24:1999–2012 Jakalian A, Jack DB, Bayly CI (2002) Fast, efficient generation of high-quality atomic charges. AM1-BCC model: II. Parameterization and validation. J Comput Chem 23:1623–1641 Wang J, Wolf RM, Caldwell JW, Kollman PA, Case DA (2004) Development and testing of a general amber force field. J Comput Chem 25:1157–1174 Jorgensen WL, Chandrasekhar J, Madura JD, Impey RW, Klein ML (1983) Comparison of simple potential functions for simulating liquid water. J Chem Phys 79:926–935 Ryckaert JP, Ciccotti G, Berendsen HJC (1977) Numerical integration of the Cartesian equations of motion of a system with constraints: molecular dynamics of n-alkanes. J Comput Phys 23:327–341 Essmann U, Perera L, Berkowitz ML, Darden T, Lee H, Pedersen LG (1995) A smooth particle mesh Ewald method. J Chem Phys 103:8577–8593 Wang F, Ma Z, Li Y, Zhu S, Xiao Z, Zhang H, Wang Y Development of in silico models for pyrazoles and pyrimidine derivatives as cyclin-dependent kinase 2 inhibitors. J Mol Graph Model 30:67–81 Liu J, Zhang H, Xiao Z, Wang F, Wang X, Wang Y (2011) Combined 3D-QSAR, Molecular Docking and Molecular Dynamics Study on Derivatives of Peptide Epoxyketone and Tyropeptin-Boronic Acid as Inhibitors Against the β Subunit of Human 20S Proteasome. Int J Mol Sci 12:1807–1835 Vulpetti A, Bosotti R (2004) Sequence and structural analysis of kinase ATP pocket residues. II Farmaco 59:759–765