Manipal Journal of Pharmaceucal Sciences | September 2020 | Volume 6 | Issue 2 65 Abstract A series of diaminopyrimidine derivatives were collected from various literatures and subjected to in silico 3D-QSAR (Quantitative Structure Activity Relationship) study. Particularly, a field-based 3D-QSAR by using steric, H-bond donor and acceptor, electrostatic and hydrophobic fields was performed. Field-based 3D-QSAR models were generated with a data set of 62 diaminopyrimidine derivatives using its inhibitory activity against focal adhesion kinase (FAK) receptor. The data set were comprised of the training set (45 compounds; 72%) and test set (17 compounds; 28%), which were randomly assigned. The relationship of molecular field versus inhibitory activity was established using the partial least square (PLS) method. Assessment of the model for stability and predictability was done by the leave-one-out method. The QSAR model (PLS factor 6) with good predictivity showed R² value 0.98, value 0.6321, Q² value 0.6151, RMSE value 0.6, SD value 0.1527 based on different fields. Key words: Anticancer agents, Diaminopyrimidines, FAK inhibitors, Field-based 3D QSAR Mirza Sohanabibi R 1 , Azmin M. Mogal1, Uam A. More 1 , Malleshappa N Noolvi 1 , Patel Salman Ismail 1 , Vaja Mukesh Bharatbhai 1 , Mahirakhatun R Rana 1 , Payal S Jain 1,3 , Navdeep Singh Sethi 2 , Vandana Kharb 4 1 Department of Pharmaceucal Chemistry, Shree Dhanvantary Pharmacy College, Kim (Surat) - 394110, Gujarat, India. 2 Department of Pharmaceucal Chemistry, Doaba College of Pharmacy, Kharar (Mohali) - 140103, Punjab, India 3 Pharmacy, Gujarat Technological University, Chandkheda, Ahmedabad-382424, Gujarat, India., 4 Department of Pharmaceucs, Sachdeva College of Pharmacy, Ghrauan (Mohali) - 140413, Punjab, India., * Corresponding Author Date of Submission: 18 Jul 2020, Date of Revision: 31 Jul 2020, Date of Acceptance: 01 Aug 2020 How to cite this article: Sohanabibi R M, M. Mogal A, More UA, N Noolvi M, Ismail P S, BharatbhaiV M, Mahirakhatun R Rana, Jain P, Sethi N S, Kharb V: In Silico Study of FAK Inhibitors Containing Pyrimidine Fragment as Anticancer Agents. MJPS 2020; 6(2): 65-78. Research Article In Silico Study of FAK Inhibitors Containing Pyrimidine Fragment as Ancancer Agents Mirza Sohanabibi R, Azmin M. Mogal, Uam A. More*, Malleshappa N. Noolvi, Patel Salman Ismail, Vaja Mukesh Bharatbhai, Mahirakhatun R Rana, Payal S Jain, Navdeep Singh Sethi, Vandana Kharb Email: [email protected]1. Introduction “Cancer” comes from the term “karkinos” which means crab in Greek. The hallmark of cancer is the growth of dangerous, carcinogenic irregular cells within the body. 1,2 During tumour growth, the cells divide and multiply progressively, forming aggressive tumours, which may invade adjacent typical cells while continuing to spread to distant areas within the body through lymphatic circulation. 3 Even though cancer may originate within a group of cells, it is often accompanied by spread to other organs resulting in tumours in the breast, lungs, liver, stomach and colorectal. 4,5 Malignant growth is regarded as the main source of death in the world today. 6 Malignancy is the second most normal reason for mortality after coronary illness. 7 Malignancy is one of the several driving reasons for death before the age of 70 in the majority of the countries. 6 According to the GLOBCAN 2018 database, there will be an estimated 18 million new cases with almost 10 million deaths in the near future. 6 The different targets used to treat cancer are HER-2, Flk-1/KDR (VEGFR-2), Bcr-Abl, EGFR, Src kinase, focal adhesion kinase (FAK), farnesyl transferases (FTases), histone deacetylases (HDACs), telomerases; and signal transducers and
14
Embed
In Silico Study of FAK Inhibitors Containing Pyrimidine ......Mukesh Bharatbhai 1, Mahirakhatun R Rana, Payal S Jain1,3, Navdeep Singh Sethi2, Vandana Kharb4 1 Department of Pharmaceutical
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Sohanabibi RM, et al: In Silico Study of FAK Inhibitors Containing Pyrimidine...
Manipal Journal of Pharmaceutical Sciences | September 2020 | Volume 6 | Issue 2 65
AbstractA series of diaminopyrimidine derivatives were collected from various literatures and subjected to in silico 3D-QSAR (Quantitative Structure Activity Relationship) study. Particularly, a field-based 3D-QSAR by using steric, H-bond donor and acceptor, electrostatic and hydrophobic fields was performed. Field-based 3D-QSAR models were generated with a data set of 62 diaminopyrimidine derivatives using its inhibitory activity against focal adhesion kinase (FAK) receptor. The data set were comprised of the training set (45 compounds; 72%) and test set (17 compounds; 28%), which were randomly assigned. The relationship of molecular field versus inhibitory activity was established using the partial least square (PLS) method. Assessment of the model for stability and predictability was done by the leave-one-out method. The QSAR model (PLS factor 6) with good predictivity showed R² value 0.98, value 0.6321, Q² value 0.6151, RMSE value 0.6, SD value 0.1527 based on different fields.
Key words: Anticancer agents, Diaminopyrimidines, FAK inhibitors, Field-based 3D QSAR
Mirza Sohanabibi R1, Azmin M. Mogal1, Uttam A. More1, Malleshappa N Noolvi1, Patel Salman Ismail1, Vaja Mukesh Bharatbhai1, Mahirakhatun R Rana1, Payal S Jain1,3, Navdeep Singh Sethi2, Vandana Kharb4
1 Department of Pharmaceutical Chemistry, Shree Dhanvantary Pharmacy College, Kim (Surat) - 394110, Gujarat, India.
2 Department of Pharmaceutical Chemistry, Doaba College of Pharmacy, Kharar (Mohali) - 140103, Punjab, India
4 Department of Pharmaceutics, Sachdeva College of Pharmacy, Ghrauan (Mohali) - 140413, Punjab, India.,
* Corresponding Author
Date of Submission: 18 Jul 2020, Date of Revision: 31 Jul 2020, Date of Acceptance: 01 Aug 2020
How to cite this article: Sohanabibi R M, M. Mogal A, More UA, N Noolvi M, Ismail P S, BharatbhaiV M, Mahirakhatun R Rana, Jain P, Sethi N S, Kharb V: In Silico Study of FAK Inhibitors Containing Pyrimidine Fragment as Anticancer Agents. MJPS 2020; 6(2): 65-78.
Research Article
In Silico Study of FAK Inhibitors Containing Pyrimidine Fragment as Anticancer AgentsMirza Sohanabibi R, Azmin M. Mogal, Uttam A. More*, Malleshappa N. Noolvi, Patel Salman Ismail, Vaja Mukesh Bharatbhai, Mahirakhatun R Rana, Payal S Jain, Navdeep Singh Sethi, Vandana Kharb
1. Introduction“Cancer” comes from the term “karkinos” which means crab in Greek. The hallmark of cancer is the growth of dangerous, carcinogenic irregular cells within the body.1,2 During tumour growth, the cells divide and multiply progressively, forming
aggressive tumours, which may invade adjacent typical cells while continuing to spread to distant areas within the body through lymphatic circulation.3 Even though cancer may originate within a group of cells, it is often accompanied by spread to other organs resulting in tumours in the breast, lungs, liver, stomach and colorectal.4,5 Malignant growth is regarded as the main source of death in the world today.6
Malignancy is the second most normal reason for mortality after coronary illness.7 Malignancy is one of the several driving reasons for death before the age of 70 in the majority of the countries.6 According to the GLOBCAN 2018 database, there will be an estimated 18 million new cases with almost 10 million deaths in the near future.6
The different targets used to treat cancer are HER-2, Flk-1/KDR (VEGFR-2), Bcr-Abl, EGFR, Src kinase, focal adhesion kinase (FAK), farnesyl transferases (FTases), histone deacetylases (HDACs), telomerases; and signal transducers and
Sohanabibi RM, et al: In Silico Study of FAK Inhibitors Containing Pyrimidine...
66 Manipal Journal of Pharmaceutical Sciences | September 2020 | Volume 6 | Issue 2
activators of transcription (STATs). Among them, FAK is one of the promising targets to develop new anticancer agents.
FAK was found by two combining lines of research in the mid ‘90s. Initially, a roughly 120 kDa protein was found as a noteworthy integrin-subordinate tyrosine phosphorylated protein. This protein was named as FAK due to its prominent localization in the focal adhesions.8,9 FAK is a significant supporter of development factor receptor and integrin signal, playing an important role in ordinary and malignant growth cells like controlling cell development and survival.10 FAK articulation and Y397 phosphorylation are over-expressed in a few advanced organized, strong malignant growths.11 It is accepted that FAK signaling is compartmentalized through communications with a variety of effector proteins that confine FAK either to central grips, endocytic vesicles, cell–cell intersections or the cell nucleus.10 The schematic representation of FAK is given in Figure 1.
Figure 1: Schematic representation of FAK
Working as a non-receptor tyrosine kinase FAK is universally communicated, directing explicit signals starting in the regions of the attachment between integrin and extracellular matrix (ECM),12 just as those activated by initiated development factor receptors.13 Assessment of human tumours has found that upgraded articulation of FAK transcripts,14 protein and expanded FAK activity 15,16 are emphatically related with metastasis, and are frequently connected with less fortunate clinical outcomes.17 Since FAK is over-expressed in several malignancies such as the tumours of the pancreas, ovaries, cervix, prostate, breast, lung, and kidney, it shows an explicable remedial potential.18,19
Until now, compounds like PF-573228, PF-562271, TAE-226, CEP-37440, and Defactinib which demonstrated reversible ATP-competitive FAK inhibition have been developed to bind to the tumour and endothelial cells.20
The present study aims to build predictive field-based 3D-QSAR model in order to gain knowledge of the structure-activity relationship for pyrimidine-based compounds as FAK inhibitors. The 3D-QSAR method was used to produce the best predictive model and this method allowed us to use various descriptors like electrostatic, hydrophobic, or steric fields, for a set of aligned ligands.
2. Materials and method3D-QSAR study for collected data set was carried out using Schrodinger suite.
2.1 DatasetUsing a data set of 62 diaminopyrimidine derivatives (Table 1) and extracting their inhibitory activity against FAK receptor (as per literature reports), field-based 3D-QSAR models were created.20-26 A training set of 45 compounds (72%) and a test set of 17 compounds (28%) were created by random assignment. The half maximal inhibitory concentrations (IC₅₀) of these 62 derivatives were converted into negative logarithmic units and represented as pIC₅₀. The relationship between the observed activity and the predicted activity was constructed using a scatter plot.
2.2. Ligands preparation and molecular alignmentThe 2D structures of the selected compounds were converted to 3D structures and energy minimization was performed by applying a force field with the help of OPLS-2003. From which, one compound having higher bioactivity with lowest energy conformation was selected as a template for molecular alignment with other compounds from dataset.27
2.3 Field-based 3D QSAR studySix 3D models were generated using various Gaussian fields, among which the one that showed maximum statistical robustness was selected and was visualized using a 3D contour map.
Sohanabibi RM, et al: In Silico Study of FAK Inhibitors Containing Pyrimidine...
Manipal Journal of Pharmaceutical Sciences | September 2020 | Volume 6 | Issue 2 67
Ligand-targeted drug design strategies like QSAR and mapping of the pharmacophore are utilized in the process of drug discovery in several ways: activity prediction of unique molecules, searching the database for new possible molecules and important structural feature identification for activity.28 Comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) are used to design the QSAR models. These models are the execution of the CoMFA and CoMSIA strategies with a particular arrangement of parameters. Field-based model is an arrangement subordinate strategy in which the communication vitality terms are related to biological reactions, utilizing multivariate measurable investigations. Field-based QSAR enables you to fabricate 3D-QSAR models dependent on fields, for example, electrostatic, hydrophobic, or steric fields, for a lot of adjusted ligands. Using 3D-QSAR model, the field around the ligand entity was determined.27-29
In Gaussian 3D-QSAR model collaboration, energy calculations were performed using various fields making use of Gaussian settings for field calculations. While the various fields represent the autonomic factors, the key variable is the pIC50 of the model.27 The various Gaussian contour maps are explained below.
2.3.1 Steric contour mapsThe green and yellow contour regions represent the steric field. While the green coloured sections demonstrate increased activity due to expanded mass, the yellow sections propose regions where expanded steric mass is sans any activity.
2.3.2 Electrostatic contour mapsThe blue and red contour areas demonstrate the electrostatic field, where the electronic substituents resulted in a negative to activity (represented by blue colour) while the electronegative substituents were positive to activity (represented by red colour).
2.3.3 Hydrophobic contour mapsHere, the yellow and white regions indicate the hydrophobic fields, where yellow areas are good while white areas are troublesome for activity.
2.3.4 Hydrogen bond acceptor (HBA) contour mapsThe colours red and magenta indicate the HBA fields, where red areas are good and magenta areas are troublesome for activity, respectively.
2.4.5 Hydrogen bond donor (HBD) contour mapsHere, purple, and cyan regions exhibit the HBD fields, where purple regions are positive and cyan regions are negative towards the activity, respectively.
2.4.6 Gaussian aromatic ring contour mapsIt is represented as orange and gray regions, where orange regions are good for activity and gray regions are not suitable for the activity.27, 29
2.4 QSAR model validationA leave-one-out (LOO) cross-validation strategy was used to assess the stability and predictability of the developed models. In this approach, one compound is omitted from the dataset, and using the model which is derived from the remaining compounds in the dataset, the activity of the omitted compound is predicted.27 The criteria for model selection was the robustness of the model. A test (external) set comprising of 17 compounds was used to validate the models. The root mean square error (RMSE), Q², and Pearson’s r (indicative of the predictability of the model) values were the key statistical parameters that were evaluated. Contour maps were used to demonstrate the various Gaussian fields [29].
2.4.1 Partial least square (PLS) analysisThe relationship between the molecular fields and the compound inhibitory activity was established using the PLS method. pIC₅₀ values were utilized as needy factors to produce 3D-QSAR models based on regression analysis of the autonomous factors, which comprised of the field and Gaussian forces. LOO method was used to validate the model.27
#Factors: It indicates the number of factors in the model.SD: Standard deviation of the regressionR²: Correlation coefficient of activity (experimentally observed and predicted) of the training set: Cross-validated value (LOO)R² Scramble: Average value of R² from a series of models built using scrambled activitiesStability: Effect of changes in the training set composition on the stability of model predictions
Sohanabibi RM, et al: In Silico Study of FAK Inhibitors Containing Pyrimidine...
68 Manipal Journal of Pharmaceutical Sciences | September 2020 | Volume 6 | Issue 2
F: Ratio of varianceP: Statistical significanceQ²: Q² value of test setRMSE: Root mean square errorPearson’s r: Correlation coefficient of activity (experimentally observed and predicted) of the test setVariance in the information is evaluated by the correlation coefficient ‘r’, that estimates in what way the observed data tracks the fitted regression line. Mistakes in either the model or in the data will prompt a lack of fit. The predictive correlation coefficient is based on the test set’s molecules and determined by the following equation (1).30
7
omitted compound is predicted.27 The criteria for model selection was the robustness of the model.
A test (external) set comprising of 17 compounds was used to validate the models. The root mean
square error (RMSE), Q², and Pearson’s r (indicative of the predictability of the model) values
were the key statistical parameters that were evaluated. Contour maps were used to demonstrate
the various Gaussian fields [29].
2.4.1 Partial least square (PLS) analysis
The relationship between the molecular fields and the compound inhibitory activity was
established using the PLS method. pIC₅₀ values were utilized as needy factors to produce 3D-
QSAR models based on regression analysis of the autonomous factors, which comprised of the
field and Gaussian forces. LOO method was used to validate the model.27
#Factors: It indicates the number of factors in the model.
SD: Standard deviation of the regression
R²: Correlation coefficient of activity (experimentally observed and predicted) of the training set
𝑅𝑅𝑐𝑐𝑐𝑐2 : Cross-validated value (LOO)
R² Scramble: Average value of R² from a series of models built using scrambled activities
Stability: Effect of changes in the training set composition on the stability of model predictions
F: Ratio of variance
P: Statistical significance
Q²: Q² value of test set
RMSE: Root mean square error
Pearson’s r: Correlation coefficient of activity (experimentally observed and predicted) of the test
set
Variance in the information is evaluated by the correlation coefficient ‘r’, that estimates in what
way the observed data tracks the fitted regression line. Mistakes in either the model or in the data
will prompt a lack of fit. The predictive correlation coefficient is based on the test set’s molecules
and determined by the following equation (1).30
𝑟𝑟𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝2 = 𝑆𝑆𝑆𝑆−𝑃𝑃𝑃𝑃𝑃𝑃𝑆𝑆𝑆𝑆
𝑆𝑆𝑆𝑆 (1)
7
omitted compound is predicted.27 The criteria for model selection was the robustness of the model.
A test (external) set comprising of 17 compounds was used to validate the models. The root mean
square error (RMSE), Q², and Pearson’s r (indicative of the predictability of the model) values
were the key statistical parameters that were evaluated. Contour maps were used to demonstrate
the various Gaussian fields [29].
2.4.1 Partial least square (PLS) analysis
The relationship between the molecular fields and the compound inhibitory activity was
established using the PLS method. pIC₅₀ values were utilized as needy factors to produce 3D-
QSAR models based on regression analysis of the autonomous factors, which comprised of the
field and Gaussian forces. LOO method was used to validate the model.27
#Factors: It indicates the number of factors in the model.
SD: Standard deviation of the regression
R²: Correlation coefficient of activity (experimentally observed and predicted) of the training set
𝑅𝑅𝑐𝑐𝑐𝑐2 : Cross-validated value (LOO)
R² Scramble: Average value of R² from a series of models built using scrambled activities
Stability: Effect of changes in the training set composition on the stability of model predictions
F: Ratio of variance
P: Statistical significance
Q²: Q² value of test set
RMSE: Root mean square error
Pearson’s r: Correlation coefficient of activity (experimentally observed and predicted) of the test
set
Variance in the information is evaluated by the correlation coefficient ‘r’, that estimates in what
way the observed data tracks the fitted regression line. Mistakes in either the model or in the data
will prompt a lack of fit. The predictive correlation coefficient is based on the test set’s molecules
and determined by the following equation (1).30
𝑟𝑟𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝2 = 𝑆𝑆𝑆𝑆−𝑃𝑃𝑃𝑃𝑃𝑃𝑆𝑆𝑆𝑆
𝑆𝑆𝑆𝑆 (1)
Where SD is the sum of squared deviation between the biological activities of the test set and the training set molecules. A predictive sum of square (Press) is the sum of the squared deviation between the real and predicted activities of each compound in the test set.
A sign of the presentation of the model is acquired from the cross-validation (or prescient) r² which is characterized as equation (2).30
8
Where SD is the sum of squared deviation between the biological activities of the test set and the
training set molecules. A predictive sum of square (Press) is the sum of the squared deviation
between the real and predicted activities of each compound in the test set.
A sign of the presentation of the model is acquired from the cross-validation (or prescient) r² which
Sohanabibi RM, et al: In Silico Study of FAK Inhibitors Containing Pyrimidine...
74 Manipal Journal of Pharmaceutical Sciences | September 2020 | Volume 6 | Issue 2
3. Results and discussionCollected all data set molecules aligned with the lowest energy conformation and highest activity compound 32 as a template (see figure 2 for the alignment of the dataset). Further in subsections, the results which are the outcome of this database alignment have been discussed.
Figure 2: Database alignment by using compound no. 32 as template
3.1 Gaussian 3D QSAR analysisTable 2 shows the data obtained from the training and test set of compounds which were used to build the model. The best model (PLS factor 6) indicated good predictive power (Table 3). The predictability of the models was assessed using the LOO method. In Gaussian field-based QSAR PLS factor no. 6 showed the best model, which had the lowest standard deviation and best R² value 0.98 with the highest F value 464.3. In case of the test set, we got the q² value as 0.6151 and Pearson-r value of more than 0.7.
3.2 Validation of the modelsThe scatter plot for the predicted activity versus the actual activity of both training as well as test set compounds is shown in Figure 3, which is based on
Sohanabibi RM, et al: In Silico Study of FAK Inhibitors Containing Pyrimidine...
Manipal Journal of Pharmaceutical Sciences | September 2020 | Volume 6 | Issue 2 75
factor 6 of the 3D-QSAR model described above. The validated model showed strong predictive power and provided comprehensive information with regards to the structural characteristics of the compounds. Table 4 shows the Gaussian field contributions.
Figure 3: Scatter plot between predicted activity and actual activity
3.3 Contour map analysisThe correlation between the biological activity and the contour maps are shown in figures 4, 5, 6, 7, 8, and 9. The contour map for the steric fields of the compound revealed that the green colour adjacent to phenyl ring attached to pyrimidine, and nearer to the squaramide group indicated the presence of bulky groups or substituents that may be favorable for activity with the pIC₅₀ value of 9.22 in compound no 32. In compound 32, the electrostatic contours of the compound revealed that electronegative groups (red colour region) attached to squaramide group favour the activity, while blue contour near to pyrimidine ring and around the squaramide group suggest an enhanced activity due to electropositive charges. In general, a pyrimidine ring in all the compounds is highly responsible for establishing hydrogen bonding with the Cys 502. Red contour
adjacent to the N-methylbenzamide, squaramide group and around the pyrimidine ring indicates that the presence of hydrogen bond acceptor subunits increases the activity. The HBD contour map of the compound suggests that purple contour around the diphenylpyrimidine ring and N-methylbenzamide group favoured the activity, while the distant cyan contour indicates that they are highly active. The hydrophobic contour map of the compound suggests that yellow colour contour around the phenyl ring of N-methylbenzamide moiety and substituent attached to squaramide group has better activity. In compound 32, the aromatic ring contour map suggests that orange contour around the N-methylbenzamide favour the activity which is present in almost all the compounds.
Figure 4: Steric contour map of compound no. 32
Figure 5: Electrostatic contour map of compound no. 32
Table 4: Field contributions of Gaussian fields
# FactorsGaussian
steric
Gaussian
electrostatic Gaussian hydrophobic
Gaussian H-bond
acceptor
Gaussian H-bond
donor
1 0.322 0.095 0.227 0.191 0.165
2 0.285 0.112 0.233 0.207 0.164
3 0.266 0.116 0.242 0.221 0.155
4 0.266 0.115 0.249 0.213 0.157
5 0.266 0.115 0.253 0.212 0.154
6 0.269 0.115 0.256 0.213 0.147
Sohanabibi RM, et al: In Silico Study of FAK Inhibitors Containing Pyrimidine...
76 Manipal Journal of Pharmaceutical Sciences | September 2020 | Volume 6 | Issue 2
Figure 6: Hydrogen bond acceptor contour map of the compound no. 32
Figure 8: Hydrophobic contour map of compound no. 32
4. Structure activity relationship (SAR)
On the basis of 3D-QSAR contour map analysis, we came up with the SAR figure. The structure activity analysis showed that 2, 4-diaminopyrimidine is the basic scaffold for activity against focal adhesion kinase. Substitution on R1 and R2 position with alkyl chain may decrease the activity. Aromatic and heterocyclic ring substitution at R1 and R2 position may enhance the activity. Also, at C5 position of pyrimidine fragment, electron withdrawing group like Cl and NO2 are essential for activity.
Figure 10: Structure activity relationship of compound no. 32
Figure 7: Hydrogen bond donor contour map of compound no. 32
Figure 9: Aromatic ring contour map of compound no. 32
Sohanabibi RM, et al: In Silico Study of FAK Inhibitors Containing Pyrimidine...
Manipal Journal of Pharmaceutical Sciences | September 2020 | Volume 6 | Issue 2 77
5. ConclusionIn silico study was carried out for the database of compounds collected from the different literature sources. A Gaussian-based 3D-QSAR model was developed which demonstrated satisfactory predictability and indicated the potential structural characteristics which influenced the FAK inhibitory activity. The developed QSAR model was robust and showed R²₌ 0.98, ₌0.6321, Q²₌ 0.6151, RMSE₌ 0.6, SD₌ 0.1527 based on different fields. The contour map analysis revealed that the aromatic substituent at C2 and C4 positions of diaminopyrimidine enhances the activity, while alkyl groups are not that favorable for the activity. Similarly, an aromatic ring with electronegative substituent at C2 position of diaminopyrimidine and an aromatic ring with electropositive groups at C4 position of diaminopyrimidine may enhance the activity. Notably, all the diaminopyrimidine derivatives from the database was shown to form three strong hydrogen bonds. These hydrogen bonds were formed due to the bonding between, the nitrogen atom of the pyrimidine fragment and amino acid of Cys502, the amino acid of SP564 and the carbonyl group of N-methylbenzamide, and the nitrogen atom of N-methylbenzamide with Arg505, Gly506, and water. In addition, C-5 substituent of the pyrimidine fragment, including the NO2 and Cl groups produced equivalent interaction forces with the amino acid Gln438. At C-5 position, chlorine atom produces van der waal contacts with the amino acid Met499. Due to the similar binding forces, most of the inhibitors from the dataset displayed potent inhibitory potency against FAK kinase activity.
Acknowledgements
This study is financially supported by Research Promotion Scheme, All India Council for Technical Education, New Delhi, file no. 8-246/RIFD/RPS(POLICY-1)2018-19. The authors would like to thank Dr N D Jivani, President, Dr Dinesh L Sutaria, Secretary, Shree Dhanvantary Pharmacy College for providing necessary facilities to complete this work.
References1. Long ER. A History of Pathology. Dover
Publications; 1965.
2. Schlotterbeck ME. What Is Cancer? The American Cancer Society invites all nurses to see and use this new teaching film. Am. J. Nurse 1949; 49(5):300.
3. Rajpal S, Kumar A, Joe W. Economic burden of cancer in India: Evidence from cross-sectional nationally representative household survey, 2014. PloS one 2018; 13(2): e0193320.
4. Joseph F, Fraumeni J. Epidemiologic approaches to cancer etiology. Annu. Rev. Public Health 1982; 3:80-100.
5. World Health Organization. Cancer control: knowledge into action, WHO guide for effective programmes. WHO. 2006; Geneva. https://www.who.int/cancer/modules/en/
6. GLOBOCAN-2018, “International agency for Research on cancer”, http://www.iarc.fr/index.php
7. Jemal A, Siegel R, Ward E, Murray T, Xu J, Thun MJ. Cancer statistics. CA: Cancer J. Clin 2007; 57:43-66.
8. Guan JL, Trevithick JE, Hynes RO. Fibronectin/integrin interaction induces tyrosine phosphorylation of a 120-kDa protein. Cell Regul 1991; 2(11):951–964.
9. Kornberg LJ. Signal transduction by integrins: increased protein tyrosine phosphorylation caused by clustering of beta 1 integrins. Proc. Natl Acad. Sci. USA 1991; 88 (19):8392–8396.
11. Sulzmaier FJ, Jean C, Schlaepfer DD. FAK in cancer: mechanistic findings and clinical applications. Nat. Rev. Cancer 2014; 14:598–610.
12. Frame MC, Patel H, Serrels B, Lietha D, Eck MJ. The FERM domain: organizing the structure and function of FAK. Nat Rev Mol Cell Biol 2010; 11:802-814.
13. Chen HC, Chan PC, Tang MJ, Cheng CH, Chang TJ. Tyrosine phosphorylation of focal adhesion kinase stimulated by hepatocyte growth factor leads to mitogen-activated protein kinase activation. J Biol Chem 1998; 273:25777-25782.
14. Weiner TM, Liu ET, Craven RJ, Cance WG. Expression of focal adhesion kinase gene and invasive cancer. Lancet 1993; 342:1024-1025.
15. Owens LV, Xu L, Craven RJ, Dent GA, Weiner TM, Kornberg L, Liu ET, & Cance WG.
Sohanabibi RM, et al: In Silico Study of FAK Inhibitors Containing Pyrimidine...
78 Manipal Journal of Pharmaceutical Sciences | September 2020 | Volume 6 | Issue 2
Overexpression of the focal adhesion kinase (p125FAK) in invasive human tumors. Cancer Re, 1995; 55:2752-2755.
16. Hess AR, Postovit LM, Margaryan NV, Seftor EA, Schneider GB, Seftor RE, Nickoloff BJ, & Hendrix MJ. Focal adhesion kinase promotes the aggressive melanoma phenotype. Cancer Res 2005; 65:9851-9860.
17. Pylayeva Y, Gillen K M, Gerald W, Beggs HE, Reichardt LF, & Giancotti FG. Ras- and PI3K-dependent breast tumorigenesis in mice and humans requires focal adhesion kinase signaling. J Clin Invest 2009; 119:252-266
18. Zhang J, He DH, Zajac-Kaye M, Hochwald SN. A small molecule FAK kinase inhibitor, GSK2256098, inhibits growth and survival of pancreatic ductal adenocarcinoma cells. Cell cycle. Pharmacol. Ther 2014; 13(19): 3143-3149.
19. Golubovskaya VM. Focal adhesion kinase as a cancer therapy target. Anti-cancer Agents Med. Chem 2010; 10; 735-741.
20. Su Y, Li R, Ning X, Lin Z, Zhao X, Zhou J, Liu J, Jin Y, Yin Y, “Discovery of 2, 4-diarylaminopyrimidine derivatives bearing dithiocarbamate moiety as novel FAK inhibitors with antitumor and anti-angiogenesis activities.” Eur. J. Med. Chem., 2019, 177, 32-46.
21. Liu H, Wu B, Ge Y, Huang J, Song S, Wang C, Yao J, Liu K, Li Y, Li Y, Ma X. Phosphamide-containing diphenylpyrimidine analogues (PA-DPPYs) as potent focal adhesion kinase (FAK) inhibitors with enhanced activity against pancreatic cancer cell lines. Bioorg. Med. Chem. Lett 2017; 25(24):6313-21.
22. Wang L, Ai M, Yu J, Jin L, Wang C, Liu Z, Shu X, Tang Z, Liu K, Luo H, Guan W. Structure-based modification of carbonyl-diphenylpyrimidines (Car-DPPYs) as a novel focal adhesion kinase (FAK) inhibitor against various stubborn cancer cells. Eur. J. Med. Chem 2019; 172:154-62.
23. Qu M, Liu Z, Zhao D, Wang C, Zhang J, Tang Z, Liu K, Shu X, Yuan H and Ma X. Design, synthesis and biological evaluation of sulfonamide-substituted diphenylpyrimidine derivatives (Sul-DPPYs) as potent focal adhesion
24. Yen-Pon E, Li B, Acebron-Garcia-de-Eulate M, Tomkiewicz-Raulet C, Dawson J, Lietha D, Frame MC, Coumoul X, Garbay C, Etheve-Quelquejeu M, Chen H. Structure-based design, synthesis, and characterization of the first irreversible inhibitor of focal adhesion kinase. ACS Chem. Biol 2018; 13(8):2067-2073.
25. Ott GR, Cheng M, Learn KS, Wagner J, Gingrich DE, Lisko JG, Curry M, Mesaros EF, Ghose AK, Quail MR, Wan W. Discovery of clinical candidate CEP-37440, a selective inhibitor of focal adhesion kinase (FAK) and anaplastic lymphoma kinase (ALK). J. Med. Chem 2016; 59(16):7478-7496.
26. Farand J, Mai N, Chandrasekhar J, Newby ZE, Van Veldhuizen J, Loyer-Drew J, Venkataramani C, Guerrero J, Kwok A, Li N, Zherebina Y. Selectivity switch between FAK and Pyk2: Macrocyclization of FAK inhibitors improve Pyk2 potency. Bioorg. Med. Chem. Lett 2016; 26(24):5926-5930.
27. Yadav DK, Sharma P, Misra S, Singh H, Mancera RL, Kim K, Jang C, Kim MH, Pérez-Sánchez H, Choi EH, Kumar S. Studies of the benzopyran class of selective COX-2 inhibitors using 3D-QSAR and molecular docking. Arch. Pharmacal. Res 2018; 41(12);1178-1189.
28. Kalva S, Vinod D, Saleena LM, “Field-and Gaussian-based 3D-QSAR studies on barbiturate analogs as MMP-9 inhibitors. Med. Chem. Res 2013; 22(11):5303-5313.
29. Joshi SD, More UA, Dixit SR, Korat HH, Aminabhavi TM, Badiger AM. Synthesis, characterization, biological activity, and 3D-QSAR studies on some novel class of pyrrole derivatives as antitubercular agents. Med. Chem. Res 2014; 23(3):1123-1147.
30. Gade DR, Makkapati A, Yarlagadda RB, Peters GJ, Sastry BS, Prasad VR. Elucidation of chemosensitization effect of acridones in cancer cell lines: Combined pharmacophore modeling, 3D QSAR, and molecular dynamics studies. Comput. Biol. Chem 2018; 74:63-75.