2D-QSAR and docking studies of 4-anilinoquinazoline derivatives as epidermal growth factor receptor tyrosine kinase inhibitors

authors:

avatar Mahtab Ghasemi Dogaheh , avatar Heshmat Ebrahimi , avatar Fatemeh Yousefbeyk , avatar saeed ghasemi , *

Koomesh: Vol.24, issue 3; 405-420
published online: September 01, 2022
article type: Research Article
received: June 01, 2021
accepted: September 01, 2021

how to cite: Ghasemi Dogaheh M, Ebrahimi H, Yousefbeyk F, ghasemi S. 2D-QSAR and docking studies of 4-anilinoquinazoline derivatives as epidermal growth factor receptor tyrosine kinase inhibitors. koomesh. 2022;24(3):e152750. 

Abstract

Introduction: Epidermal growth factor receptor (EGFR) tyrosine kinase inhibitor derivatives play an important role in the treatment of cancer. We aim to construct 2D-QSAR models using various chemometrics using 4-anilinoquinazoline-containing EGFR TKIs. In addition, the binding profile of these compounds was evaluated using a docking study. Materials and Methods: In this study, 122 compounds of seven different structural categories with a 4-anilinoquinazoline scaffold were obtained from the pieces of literature. 2D-QSAR models were prepared using the linear method, including multiple linear regression (MLR), alongside non-linear methods, including artificial neural networks (ANN) and support vector machines (SVM). The validation of suggested 2D-QSAR models was performed using internal and external validation techniques. For molecular docking, three compounds, including compounds 32, 75, and 98, which had the highest pIC50, were used. Molecular docking was performed using AutoDock 4.2. Results: Selected descriptors indicated that atomic electronegativity, hydrogen bonding ability, lipophilicity, and molecular shape and volume were the factors affecting the activity of the compounds. Statistical criteria, related to interpretation and validation of the models, were in the appropriate range. The model obtained by the ANN method, with the lowest mean absolute error of 0.365 for both training and testing sets, was the best. The compound 75 as the most potent compound indicated binding energies of -8.33 kcal/mol. Conclusion: Finally, it was found that the models could effectively predict the activity of the compounds. It was also discovered that the compounds were properly bound to the active site of the receptor. In addition, validation results showed that all processes were sufficiently valid.

References

  • 1.

    Garofalo A, Goossens L, Lemoine A, Ravez S, Six P, Howsam M, et al. [4-(6, 7-disubstituted quinazolin-4-ylamino) phenyl] carbamic acid esters: A novel series of dual egfr/vegfr-2 tyrosine kinase inhibitors. Med Chem Comm 2011; 2: 65-72.

  • 2.

    Zhang HQ, Gong FH, Ye JQ, Zhang C, Yue XH, Li CG, et al. Design and discovery of 4-anilinoquinazoline-urea derivatives as dual tk inhibitors of egfr and vegfr-2. Eur J Med Chem 2017; 125: 245-254.

  • 3.

    Ahmadzadeh M, Farshdari F, Nematollahi L, Maleknia S, Mohit E. Optimization of her2-based and cell-based elisa for detection of trastuzumab biosimilar. Koomesh 2020; 22: 341-350. (Persian).

  • 4.

    Osmani F, Rasekhi A, Hajizadeh E, Akbari M E. Simultaneous modeling of multiple recurrences in breast cancer patients. Koomesh 2020; 22: 359-364. (Persian).

  • 5.

    Yu H, Li Y, Ge Y, Song Z, Wang C, Huang S, et al. Novel 4-anilinoquinazoline derivatives featuring an 1-adamantyl moiety as potent egfr inhibitors with enhanced activity against nsclc cell lines. Eur J Med Chem 2016; 110: 195-203.

  • 6.

    Cohen RB. Epidermal growth factor receptor as a therapeutic target in colorectal cancer. Clin Colorectal Cancer 2003; 2: 246-251.

  • 7.

    Chandregowda V, Kush A, Reddy GC. Synthesis and in vitro antitumor activities of novel 4-anilinoquinazoline derivatives. Eur J Med Chem 2009; 44: 3046-3055.

  • 8.

    Hennequin LF, Stokes ES, Thomas AP, Johnstone C, Pl PA, Ogilvie DJ, et al. Novel 4-anilinoquinazolines with c-7 basic side chains: Design and structure activity relationship of a series of potent, orally active, vegf receptor tyrosine kinase inhibitors. Eur J Med Chem 2002; 45: 1300-1312.

  • 9.

    Hansch C, Hoekman D, Gao H. Comparative qsar: Toward a deeper understanding of chemicobiological interactions. Chem Rev 1996; 96: 1045-1076.

  • 10.

    Shayanfar A, Ghasemi S, Soltani S, Asadpour-Zeynali K, J Doerksen R, Jouyban A. Quantitative structure-activity relationships of imidazole-containing farnesyltransferase inhibitors using different chemometric methods. Med Chem 2013; 9: 434-448.

  • 11.

    Kharkar PS. Two-dimensional (2d) in silico models for absorption, distribution, metabolism, excretion and toxicity (adme/t) in drug discovery. Current topics in medicinal chemistry 2010; 10: 116-126.

  • 12.

    Vapnik V. The nature of statistical learning theory: Springer science business media; 2013.

  • 13.

    Verma J, Khedkar VM, Coutinho EC. 3d-qsar in drug design-a review. Current topics in medicinal chemistry 2010; 10: 95-115.

  • 14.

    Shahlaei M, Fassihi A, Saghaie L. Application of pc-ann and pc-ls-svm in qsar of ccr1 antagonist compounds: A comparative study. Eur J Med Chem 2010; 45: 1572-1582.

  • 15.

    Cheng Z, Zhang Y, Fu W. Qsar study of carboxylic acid derivatives as hiv-1 integrase inhibitors. Eur J Med Chem 2010; 45: 3970-3980.

  • 16.

    Darnag R, Mazouz EM, Schmitzer A, Villemin D, Jarid A, Cherqaoui D. Support vector machines: Development of qsar models for predicting anti-hiv-1 activity of tibo derivatives. Eur J Med Chem 2010; 45: 1590-1597.

  • 17.

    Bridges AJ, Zhou H, Cody DR, Rewcastle GW, McMichael A, Showalter HD, et al. Tyrosine kinase inhibitors. 8. An unusually steep structureactivity relationship for analogues of 4-(3-bromoanilino)-6,7-dimethoxyquinazoline (pd 153035), a potent inhibitor of the epidermal growth factor receptor. Eur J Med Chem 1996; 39: 267-276.

  • 18.

    Palmer BD, Trumpp-Kallmeyer S, Fry DW, Nelson JM, Showalter HD, Denny WA. Tyrosine kinase inhibitors. 11. Soluble analogues of pyrrolo- and pyrazoloquinazolines as epidermal growth factor receptor inhibitors: Synthesis, biological evaluation, and modeling of the mode of binding. Eur J Med Chem 1997; 40: 1519-1529.

  • 19.

    Rewcastle GW, Denny WA, Bridges AJ, Zhou H, Cody DR, McMichael A, et al. Tyrosine kinase inhibitors. 5. Synthesis and structure-activity relationships for 4-[(phenylmethyl)amino]- and 4-(phenylamino)quinazolines as potent adenosine 5'-triphosphate binding site inhibitors of the tyrosine kinase domain of the epidermal growth factor receptor. Eur J Med Chem 1995; 38: 3482-3487.

  • 20.

    Rewcastle GW, Palmer BD, Bridges AJ, Showalter HD, Sun L, Nelson J, et al. Tyrosine kinase inhibitors. 9. Synthesis and evaluation of fused tricyclic quinazoline analogues as atp site inhibitors of the tyrosine kinase activity of the epidermal growth factor receptor. Eur J Med Chem 1996; 39: 918-928.

  • 21.

    Mauri A, Consonni V, Pavan M, Todeschini R. Dragon software: An easy approach to molecular descriptor calculations. Match 2006; 56: 237-248.

  • 22.

    Teppen BJ. Hyperchem, release 2: Molecular modeling for the personal computer. J Chem Inform Comput Sci 1992; 32: 757-759.

  • 23.

    Ghanbarzadeh S, Ghasemi S, Shayanfar A, Ebrahimi-Najafabadi H. 2d-qsar study of some 2, 5-diaminobenzophenone farnesyltransferase inhibitors by different chemometric methods. EXCLI J 2015; 14: 484.

  • 24.

    Nunes CA, Freitas MP. Aug-mia-qspr on the modeling of sweetness values of disaccharide derivatives. LWT Food Sci Technol 2013; 51: 405-408.

  • 25.

    Habibi-Yangjeh A. Qsar study of the 5-ht 1a receptor affinities of arylpiperazines using a genetic algorithm-artificial neural network model. Monatsh Chem 2009; 140: 523-530.

  • 26.

    Leardi R, Seasholtz MB, Pell RJ. Variable selection for multivariate calibration using a genetic algorithm: Prediction of additive concentrations in polymer films from fourier transform-infrared spectral data. Anal Chim Acta 2002; 461: 189-200.

  • 27.

    Soltani S, Abolhasani H, Zarghi A, Jouyban A. Qsar analysis of diaryl cox-2 inhibitors: Comparison of feature selection and train-test data selection methods. Eur J Med Chem 2010; 45: 2753-2760.

  • 28.

    Dastmalchi S, Hamzeh-Mivehroud M, Ghafourian T, Hamzeiy H. Molecular modeling of histamine h3 receptor and qsar studies on arylbenzofuran derived h3 antagonists. J Mol Graph Model 2008; 26: 834-844.

  • 29.

    Jalali-Heravi M, Asadollahi-Baboli M, Shahbazikhah P. Qsar study of heparanase inhibitors activity using artificial neural networks and levenberg-marquardt algorithm. Eur J Med Chem 2008; 43: 548-556.

  • 30.

    Ghandadi M, Shayanfar A, Hamzeh-Mivehroud M, Jouyban A. Quantitative structure activity relationship and docking studies of imidazole-based derivatives as p-glycoprotein inhibitors. Med Chem Res 2014; 23: 4700-4712.

  • 31.

    Golbraikh A, Tropsha A. Beware of q2! J Mol Graph Model 2002; 20: 269-276.

  • 32.

    Roy PP, Roy K. On some aspects of variable selection for partial least squares regression models. QSAR Comb Sci 2008; 27: 302-313.

  • 33.

    Moghadam FA, Evazalipour M, Kefayati H, Ghasemi S. 6, 7-disubstituted-4-anilinoquinazoline: Design, synthesis and anticancer activity as a novel series of potent anticancer agents. Pharm Sci 2021; 27: 209-218.

  • 34.

    Caballero J, Fernndez M, Saavedra M, Gonzlez-Nilo F D. 2d autocorrelation, comfa, and comsia modeling of protein tyrosine kinases' inhibition by substituted pyrido [2, 3-d] pyrimidine derivatives. Bioorg Med Chem 2008; 16: 810-821.

  • 35.

    Gozalbes R, Doucet J, Derouin F. Application of topological descriptors in qsar and drug design: History and new trends. Curr Drug Targets Infect Disord 2002; 2: 93-102.

  • 36.

    Khan AU. Descriptors and their selection methods in qsar analysis: Paradigm for drug design. Drug Discov Today 2016; 21: 1291-302.

  • 37.

    Seedher N, Bhatia S, Singh B. Quantitative correlation between theoretical molecular descriptors and drug-hsa binding affinities for various cox-2 inhibitors. Chem Biol Drug Des 2008; 72: 297-302.

  • 38.

    Gosav S, Praisler M, Dorohoi D. Ann expert system screening for illicit amphetamines using molecular descriptors. J Mol Struct 2007; 834: 188-194.

  • 39.

    Gramatica P, Corradi M, Consonni V. Modelling and prediction of soil sorption coefficients of non-ionic organic pesticides by molecular descriptors. Chemosphere 2000; 41: 763-777.

  • 40.

    Devinyak O, Havrylyuk D, Lesyk R. 3d-morse descriptors explained. J Mol Graph Model 2014; 54: 194-203.

  • 41.

    Gramatica P. Whim descriptors of shape. QSAR Comb Sci 2006; 25: 327-332.

  • 42.

    D'Archivio AA, Maggi MA, Mazzeo P, Ruggieri F. Quantitative structure-retention relationships of pesticides in reversed-phase high-performance liquid chromatography based on whim and getaway molecular descriptors. Anal Chim Acta 2008; 628: 162-172.

  • 43.

    Gonzlez MP, Tern C, Teijeira M, Helguera AM. Radial distribution function descriptors: An alternative for predicting a2 a adenosine receptors agonists. Eur J Med Chem 2006; 41: 56-62.##.