Next Article in Journal
Characterization of Different Functionalized Lipidic Nanocapsules as Potential Drug Carriers
Previous Article in Journal
Impaired Iron Status in Aging Research
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Pyrrolo[3,2-d]pyrimidine Derivatives as Type II Kinase Insert Domain Receptor (KDR) Inhibitors: CoMFA and CoMSIA Studies

School of Pharmaceutical Sciences, Southern Medical University, Guangzhou 510515, China
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2012, 13(2), 2387-2404; https://doi.org/10.3390/ijms13022387
Submission received: 21 December 2011 / Revised: 13 February 2012 / Accepted: 16 February 2012 / Published: 22 February 2012
(This article belongs to the Section Physical Chemistry, Theoretical and Computational Chemistry)

Abstract

:
Kinase insert domain receptor (KDR) inhibitors have been proved to be very effective anticancer agents. Molecular docking, 3D-QSAR methods, CoMFA and CoMSIA were performed on pyrrolo[3,2-d]pyrimidine derivatives as non-ATP competitive KDR inhibitors (type II). The bioactive conformation was explored by docking one potent compound 20 into the active site of KDR in its DFG-out inactive conformation. The constructed CoMFA and CoMSIA models produced statistically significant results with the cross-validated correlation coefficients q2 of 0.542 and 0.552, non-cross-validated correlation coefficients r2 of 0.912 and 0.955, and predicted correction coefficients r2pred of 0.913 and 0.897, respectively. These results ensure the CoMFA and CoMSIA models as a tool to guide the design of a series of new potent KDR inhibitors.

1. Introduction

Angiogenesis, the form of new blood vessels by capillary sprouting from pre-existing vasculatures, is a normal process for organ development during embryogenesis, wound healing and female reproductive cycling [1,2]. Abnormal regulation of angiogenesis has been shown to get involved in many diseases, such as diabetic retinopathy, psoriasis, rheumatoid arthritis, and cancer. In particular, it is widely recognized that the growth and metastasis of solid tumors is dependent on angiogenesis [3,4]. Out of the many factors that are involved in angiogenesis, vascular endothelial growth factors (VEGFs) are of particular interests [5,6]. The VEGFs are required for vasculogenesis and angiogenic sprouting, and act through receptor tyrosine kinase (VEGFR-1, -2, and -3). Among them, VEGF receptor-2 (VEGFR-2) or kinase insert domain receptor (KDR) play crucial roles in vessel sprouting and new vessel initiation in early stage of angiogenesis.
It is well known that inhibiting of KDR leads to suppression of angiogenesis and tumor growth. A number of preclinical and clinical studies have shown that many small-molecule KDR inhibitors are capable of inhibiting angiogenesis, tumor progression, and dissemination [711]. The vast majority of the KDR inhibitors known to date, such as Gefitinib, are ATP-competitive and classified as type I inhibitors. Such inhibitors target the ATP binding pocket in its active conformation of the activation loop. This conformation is normally referred to as DFG “in” based on the position of the conserved triad aspartate-phenylalanine-glycine (DFG) at the entrance of the activation loop. Type I inhibitors typically function in the DFG “in” conformation of KDR through hydrogen bonding with the backbone residues of the hinge region as well as hydrophobic interactions in and around the adenine region. These features, however, make it quite difficult to design highly selective type I inhibitors [12,13]. In addition, ATP-competitive inhibitors have to compete with high levels of intracellular ATP, leading to a significant difference between the in vitro and in vivo activities. In response to these issues, non-ATP competitive kinase inhibitors, such as Imatinib have been identified. These inhibitors, i.e., type II inhibitors bind to and stabilize an inactive kinase form that features the DFG motif in an “out” conformation. The different position of the DFG residues in the “out” form creates a new hydrophobic binding pocket that is adjacent to the ATP-binding site. This pocket, also known as the allosteric site, is characteristic of kinase in an inactive conformation. Type II inhibitors predominantly occupy the ATP binding site, but they also exploit unique hydrogen bonding and hydrophobic interactions with the allosteric site. Compared with type I kinase inhibitors, type II inhibitors have several advantages, including great cellular potency and improved kinase selectivity. In addition, type II inhibitors, because of their interactions with both the ATP pocket and the allosteric site, may provide an avenue to overcome the mutations that induce resistance to the other types of inhibitors [1418].
Recently, Oguro et al. reported a set of pyrrolo[3,2-d]pyrimidine derivatives 152 (Table 1) as potent and selective KDR inhibitors in biochemical and cellular assays [19,20]. These compounds are structurally distinct from the ATP-competitive inhibitors, and the X-ray crystallographic analysis has shown that compound 27 induces an inactive conformation, in which the diphenylurea moiety occupies the hydrophobic pocket created by the conformation change of DFG motif (DFG-out) [19,20]. In the present study, to gain further insight into the interactions of these compounds with the KDR, we employed the molecular docking-guided three dimensional quantitative structure-activity relationship (3D-QSAR) study, comparative molecular field analysis (CoMFA) [21] and comparative molecular similarity indices analysis (CoMSIA) [22] to address how steric, electrostatic, hydrophobic and hydrogen-bonding interactions modulate the inhibitory activities.

2. Materials and Methods

2.1. General

The crystallographic coordinates of KDR in complex with small-molecule inhibitors were obtained from the Brookheaven Protein Databank as entries 2OH4 [23]. All the molecular modeling and calculations were performed using Sybyl 7.3 molecular modeling package [24].

2.2. Data Sets

In this study, 52 pyrrolo[3,2-d]pyrimidine derivatives were taken from the work of Oguro et al. [19,20]. Their structures and inhibitory activities are listed in Tables 1 and 2. The IC50 values (M) were converted to the corresponding pIC50 (= −logIC50). On the basis of the diversity in the structures and activities, these 52 compounds were divided into two groups: 44 compounds were used as the training set to build the 3D-QSAR models, and 8 compounds that are marked with an asterisk in Table 1 were used as the test set to evaluate the predictive power of the developed CoMFA and CoMSIA models.

2.3. Molecular Docking

Because the crystal structure of KDR in complex with pyrrolo[3,2-d]pyrimidine was not available in the Brookheaven Protein Databank (PDB), the bioactive conformation was simulated by docking using Surflex-dock program. The crystallographic coordinates of KDR in complex with its inhibitor, which was reported to be in the inactive DFG-out conformation of KDR, were obtained from the PDB as entries 2OH4 [23]. Surflex-Dock program [25,26] has been widely used to calculate the protein–ligand interactions, and to efficiently predict the active conformations [2735]. Surflex-Dock uses a Protomol-based method and an empirical scoring function to dock a ligand into the binding site of a receptor. The Protomol is an idealized representation of a ligand that forms every potential interaction with the binding site. Surflex-Dock’s scoring function contains the factors that play crucial roles in the ligand-receptor interaction, including hydrophobic, polar, repulsive, entropic and solvation terms. In this study, the Protomol was generated using a ligand-based approach. During the Protomol generation process, two particular parameters, Protomol_bloat and Protomol_threshold, were specified to form the appropriate binding pocket. The former determines how far the site should extend from a potential ligand, whereas the latter determines how deep the atomic probes that are used to define the Protomol can penetrate into the protein. In the present work, Protomol_bloat and Protomol_threshold default values (0 and 0.5, respectively) were used when a reasonable binding pocket was obtained. During the docking process, the default values of all the other parameters were assigned. The highest-scored conformation of a potent compound 20 based on the Surflex-Dock scoring functions, was selected as the final bioactive conformation.

2.4. Molecular Modeling

In the 3D-QSAR study, the selection of active conformations is a key step for CoMFA and CoMSIA studies. The bioactive conformation of compound 20 was simulated using Surflex-Dock. The docked conformation with the highest total score was used as the template to construct the 3D structures of the rest compounds in the data set. Structural energy minimization process was performed using the Tripos force field with a distance-dependent dielectric and Powell gradient algorithm with a convergence criterion of 0.001 kcal/mol. Partial atomic charges were calculated using Gasteiger-Hückel method.

2.5. Molecular Alignment

In the 3D-QSAR study, the alignment rule is also a key step. The predictive accuracy of the CoMFA and CoMSIA models and the reliability of the contour maps are directly dependent on the structural alignment rule. The compounds were aligned by the atomfit to the template 20. The aligned compounds are shown in Figure 1.

2.6. CoMFA and CoMSIA Studies

Standard CoMFA and CoMSIA procedures were performed. A 3D cubic lattice was created automatically by extending at least 4 Å beyond all the aligned molecules in X, Y and Z directions with 2.0 Å grid spacing. The CoMFA steric (Lennard-Jones potential) and electrostatic (Coulomb potential) fields at each lattice were calculated using the standard Tripos force field method. A distance dependent dielectric constant of 1.0 was used, and an sp3 hybridized carbon atom with one positive charge and a radius of 1.52 Å served as a probe atom to calculate the steric and electrostatic fields. The default cutoff value of 30.0 kcal/mol was adopted.
Compared with CoMFA, CoMSIA methodology has the advantage of exploring the impacts of more fields. In addition to the steric (S) and electrostatic (E) fields used in CoMFA, the CoMSIA method defines hydrophobic (H), hydrogen bond donor (D), and hydrogen bond acceptor (A) descriptors. The CoMSIA fields were derived, according to Klebe et al. [22], from the same lattice box that was used in the CoMFA calculations, with a grid spacing of 2 Å and a probe carbon atom with one positive charge and a radius of 1.0 Å as implemented in Sybyl. Arbitrary definition of cutoff limits was not required in CoMSIA method, wherein the abrupt changes of potential energy near the molecular surface were taken into account in the distance dependent Gaussian type functional form. The default value of 0.3 was used as the attenuation factor.

2.7. PLS Regression Analysis and Validation of QSAR Models

Partial least squares (PLS) approach was used to derive the 3D QSAR models. The CoMFA and CoMSIA descriptors were used as independent variables and the pIC50 values were used as dependent variables. CoMFA and CoMSIA column filtering was set to 2.0 kcal/mol to improve the signal-to-noise ratio. The leave-one-out (LOO) cross-validation was carried out to obtain the optimal number of components (N) and the correlation coefficient q2. The obtained N was then used to derive the final QSAR model and to obtain the non-cross-validation correlation coefficient r2, standard error of estimate (SEE), and Fischer (F) ratio value.
To assess the predictive power of the derived 3D-models, a set of test compounds that had known biological activities, was used to validate the obtained models. The predictive correlation r2pred value was calculated using
r 2 pred = ( SD - PRESS ) / SD
Wherein SD is the sum of the squared deviations between the biological activities of the test compounds and the mean activities of the training compounds, and PRESS is the sum of the squared deviations between the experimental and the predicted activities of the test compounds.

3. Results and Discussion

3.1. Binding Modes of Pyrrolo[3,2-d]pyrimidine Derivatives

To determine the probable binding conformations of these compounds, Surflex-Dock was used to dock one potent compound 20 into the active site of KDR (PDB code: 2OH4). First, the docking reliability was validated by a known inhibitor 53 (Figure 2) that was reported to bind in the DFG-out inactive conformation of KDR [23]. The co-crystallized 53 was re-docked into the binding site, and the docked conformation having the highest total score was selected as the most probable binding conformation (Figure 3). The low root mean-square deviation (RMSD) of 0.58 Å between the docked and the crystal conformations demonstrated the high reliability of Surflex-dock in reproducing the experimentally observed binding mode for these KDR inhibitors. As shown in Figure 3a, redocked 53 was almost in the same orientation with co-crystallized 53 at the active site of KDR. Therefore, Surflex-Dock docking protocol and the used parameters were extended in search for the binding conformations of KDR inhibitors.
The binding mode of compound 20 at the active site of KDR in its inactive conformation is shown in Figure 3b. It is clear that compound 20 adopts an overall conformation that is very similar to that of compound 53, and forms multiple hydrogen bonds with KDR. Specifically, the 1-nitrogen of the pyrrolo[3,2-d]pyrimidine subunit forms one hydrogen bond with the backbone-NH of Cys917 in the hinge region of KDR at the angle and distance of 130.76° and 3.253 Å, respectively. In the allosteric site, the two NH groups of the urea moiety form two hydrogen bonds with the side chain carboxylate of Glu883 at the angles of 99.12° and 97.85°, and the distances of 2.566 Å and 2.873 Å, respectively, whereas the CO moiety interacts with the backbone-NH of Asp1044 in the DFG motif through hydrogen bonding at the angle and distance of 151.02° and of 2.666 Å, respectively. These interactions are typical characteristics of the interactions between the inhibitors and the inactive conformation of the kinase. In addition, the terminal phenyl moiety occupies the hydrophobic pocket that is formed from residues Ile886, Leu887, Ile890, Val896 and Leu1017.

3.2. CoMFA and CoMSIA Results

The CoMFA and CoMSIA 3D-QSAR models were derived from a training set of 44 compounds. The statistical results of the CoMFA and CoMSIA 3D-QSAR models are presented in Table 3. The CoMFA model gave a cross-validated correlation coefficient q2 of 0.542, an optimal number of principal components (N) of 4 and a non-cross-validated correlation coefficient r2 of 0.912. The corresponding contributions of steric and electrostatic fields were 52.5% and 47.5%, respectively. The CoMSIA model gave a cross-validated correlation coefficient q2 of 0.552, an optimal number of principal components of 5 and a non-cross-validated correlation coefficient r2 of 0.955. The corresponding contributions of steric, electrostatic, hydrophobic, hydrogen bond donor and acceptor fields were 18.4%, 22.8%, 34.3%, 6.3% and 18.2%, respectively. Both the CoMFA and CoMSIA models were satisfactory from the viewpoint of statistical significance. The activities of the 44 training compounds were predicted with the constructed CoMFA and CoMSIA models. The predicted pIC50 values are shown in Table 2 and Figure 4. It can be seen that the predicted pIC50 values were in good agreement with the experimental values, indicating that the obtained CoMFA and CoMSIA models had strong predictive ability.

3.3. Validation of the 3D-QSAR Models

The predictive powers of the CoMFA and CoMSIA models were validated by the eight test compounds. The predicted pIC50 values were found to be in good agreement with the experimental data within an acceptable error range (Table 2 and Figure 4). The predictive correction coefficients of the CoMFA and CoMSIA models were 0.913 and 0.897, respectively. This result indicates that the CoMFA and CoMSIA models may be used to predict the inhibitory activities of novel pyrrolo[3,2-d]pyrimidine derivatives as type-II KDR inhibitors.

3.4. Contour Analysis

To visualize the results of the CoMFA and CoMSIA models, the 3D coefficient contour maps were generated. The CoMFA and CoMSIA results were graphically interpreted by the field contribution maps using the STDEV*COEFF field type. The contour maps of CoMFA (steric and electrostatic) and CoMSIA (steric, electrostatic, hydrophobic, hydrogen bond donor and acceptor fields) are shown in Figures 5 and 6, respectively. Compound 20 was displayed in the map in aid of visualization. All the contours represented the default 80% and 20% level contributions for favorable and unfavorable regions, respectively, except 90% and 10% level contributions in Figure 6c, respectively.

3.4.1. CoMFA Contour Maps

The CoMFA contour maps of the steric and electrostatic fields are shown in Figure 5. In the map of steric field, the green contours represent the regions in which bulky groups confer an increase in the activity, whereas the yellow ones represent the regions where bulky groups may lead to a decrease in the activity. Similarly, in the map of electrostatic field, the blue contours indicate the regions where electropositive substitution increases the inhibitory activity, whereas the red contours indicate the regions where electronegative substitution increases the activity.
In the CoMFA steric contour map (Figure 5a), a large green contour near the 3′-position of the terminal 3′-trifluoromethylphenyl group of compound 20 suggests that introducing of bulky groups at this position would increase the activity. In consistent with this, compounds bearing bulky groups at this position, for example, compounds 18, 20, 21 and 23 showed high activities, whereas the ones bearing small groups at the same position, for example compounds 16 and 22, showed low activities. Large yellow contours near the 4′- and 5′-positions of the terminal 3′-trifluoromethylphenyl group suggest that introducing of bulky groups at these positions would decrease the activity. This is in agreement with the fact that compounds 4452 with bulky groups at 4′- or 5′-positions of the terminal 3′-trifluoromethylphenyl showed decreased activities. In addition, big yellow contours near the pyrrolo[3,2-d]pyrimidine subunit suggest that steric bulkiness is unfavorable by the model. When pyrrolo[3,2-d]pyrimidine subunit is conjugated to the central phenyl group via a sulfur or nitrogen atom (X), the subunit falls into the region of the yellow contour, which suggests that a decrease in the inhibitory activity will be observed. For example, sulfur-linked derivative 12 (IC50 = 110 nM) and amine-linked derivative 13 (IC50 = 1400 nM) were 3- and 40-fold less active than oxygen-linked derivative 16 (IC50 = 33 nM), respectively.
In the CoMFA electrostatic contour map, a small red contour near the 3′-position of the 3′-trifluoromethylphenyl group of compound 20 indicates that introduction of electronegative groups around this position would increase the inhibitory activity. This, together with the green contour discussed above, suggests that electronegative and bulky groups near the 3′-position of the terminal phenyl group are favored by the CoMFA model. A blue contour near the 4′-position of the 3′-trifluoromethylphenyl group indicates that introducing of electropositive groups around this position would increase the inhibitory activity. For example, compound 24 with electropositive hydrogen at the 4′-position showed higher activity than the corresponding compounds 4548 with electronegative substituents.

3.4.2. CoMSIA Contour Maps

The CoMSIA steric and electrostatic contour maps of compound 20 are shown in Figure 6a,b, respectively. These contours are quite similar to those of CoMFA. Therefore, our following discussion will focus on the hydrophobic, hydrogen bond donor and acceptor fields.
Figure 6c shows the hydrophobic contour maps in which yellow and gray contours indicate the regions where hydrophobic and hydrophilic groups are favored by the model, respectively. A yellow contour near the 3′-position of the terminal 3′-trifluoromethylphenyl group indicates that hydrophobic substituent at this position would increase the activity. This has been noticed in compounds 18, 20, 21 and 23 bearing hydrophobic Cl, CF3, Br and CH3, respectively. These compounds exhibited increased activities. This hydrophobic interaction may play a crucial role in improving of the binding affinity, since it is also observed in the CoMFA and CoMSIA steric contour maps. The large gray contour near the urea indicates that hydrophilic urea at this position is favorable.
The CoMSIA hydrogen bond donor and acceptor contour plots are shown in Figure 6d,e respectively. The cyan contours represent the regions where hydrogen bond-donating groups increase the activity, whereas the purple contours represent the regions where hydrogen bond-donating groups decrease the activity. Similarly, the magenta contours indicate the regions where hydrogen bond-accepting groups increase the inhibitory activity, whereas the red contours indicate the regions where hydrogen bond-accepting groups decrease the activity.
The cyan contours near the NH of urea indicate that hydrogen bond-donating groups are favored. This is well consistent with the observations that NH group in this region forms stable hydrogen bonds with the residue of Glu883 as hydrogen bond donor, and that the amide derivatives 6 and 7 having a CH2 to replace the NH of the urea 16 led to more than 300- and 100-fold decreases in the activities, respectively. A large magenta contour located on the carbonyl oxygen of the urea moiety suggests that hydrogen bond-accepting groups are favored in this region. This is evident from the fact that thiourea derivative 8 was about 100-fold less active than the corresponding urea derivative 16, and that the carbonyl oxygen of the urea moiety was engaged in hydrogen bonding interaction with the Asp1044 in the DFG motif.

3.5. Design of New Inhibitors

As shown above, the CoMFA and CoMSIA have provided detailed insight into the key structural requirements for potent activities of the inhibitors of this class. Specifically, the urea plays a crucial role in the inhibitory activity—its replacement with thiourea or with an acetamide leads to a complete loss of the activity. Substituting of the oxygen linker with sulfur or nitrogen atoms affords less active compounds. Introducing of appropriately bulky and strongly hydrophobic groups at the 3′-position of the terminal phenyl group may significantly increase the activity. To demonstrate the practical values of these structure-activity relationships, a series of new inhibitors were designed and their pIC50 values were predicted with the established CoMFA and CoMSIA models (Table 4).

4. Concluding Remarks

In this study, 3D-QSAR analyses, CoMFA and CoMSIA, have been applied to a set of recently synthesized pyrrolo[3,2-d]pyrimidine derivatives as type II KDR inhibitors. The binding mode of the template molecule 20 was clarified by Surflex-dock. The CoMFA and CoMSIA models showed statistically significant results in terms of cross-validated coefficients and conventional coefficients. Their predictive capabilities were verified by the test compounds. The derived CoMFA and CoMSIA models showed predictive cross-validated coefficients of 0.913 and 0.897, respectively, and the activities of the training and test compounds were predicted with good accuracy. Based on the obtained structure-activity relationships, a series of new inhibitors were designed to have excellent activities that were predicted with the developed CoMFA and CoMSIA models. Thus, these models may be expected to serve as a tool to guide the future rational design of pyrrolo[3,2-d]pyrimidine-based novel type II KDR inhibitors with potent activities.

Acknowledgments

The authors are gratefully acknowledged financial support from the National Science Foundation of China (No.81173097) and the Science and Technology Planning Project of Guangzhou (2009A1- E011-8).

References

  1. Risau, W. Mechanisms of angiogenesis. Nature 1997, 386, 671–674. [Google Scholar]
  2. Ribatti, D.; Vacca, A.; Nico, B.; Roncali, L.; Dammacco, F. Postnatal vasculogenesis. Mech. Dev 2001, 100, 157–163. [Google Scholar]
  3. Folkman, J. Anti-angiogenesis: New concept for therapy of solid tumors. Ann. Surg 1972, 175, 409–416. [Google Scholar]
  4. Liotta, L.A.; Steeg, P.S.; Stetler-Stevenson, W.G. Cancer metastasis and angiogenesis: An imbalance of positive and negative regulation. Cell 1991, 64, 327–336. [Google Scholar]
  5. Yancopoulos, G.D.; Davis, S.; Gale, N.W.; Rudge, J.S.; Wiegand, S.J.; Holash, J. Vascular-specific growth factors and blood vessel formation. Nature 2000, 407, 242–248. [Google Scholar]
  6. Hartman, G.D.; Fraley, M.E.; Bilodeau, M.T. Kinase insert domain-containing receptor kinase inhibitors as antiangiogenic agents. Expert. Opin. Invest. Drugs 2002, 11, 737–745. [Google Scholar]
  7. Hicklin, D.J.; Ellis, L.M. Role of the vascular endothelial growth factor pathway in tumor growth and angiogenesis. J. Clin. Oncol 2005, 23, 1011–1027. [Google Scholar]
  8. Ferrara, N. Vascular endothelial growth factor: Basic science and clinical progress. Endocr. Rev 2004, 25, 581–611. [Google Scholar]
  9. Carmeliet, P.; Ferreira, V.; Breier, G.; Pollefeyt, S.; Kieckens, L.; Gertsenstein, M.; Fahrig, M.; Vandenhoeck, A.; Harpal, K.; Eberhardt, C.; et al. Abnormal blood vessel development and lethality in embryos lacking a single VEGF allele. Nature 1996, 380, 435–439. [Google Scholar]
  10. Ferrara, N.; Davis-Smyth, T. The biology of vascular endothelial growth factor. Endocr. Rev 1997, 18, 4–25. [Google Scholar]
  11. Holmes, K.; Roberts, O.L.; Thomas, A.M.; Cross, M.J. Vascular endothelial growth factor receptor-2: Structure, function, intracellular signaling and therapeutic inhibition. Cell Signal 2007, 19, 2003–2012. [Google Scholar]
  12. Fischer, P.M. The design of drug candidate molecules as selective inhibitors of therapeutically relevant protein kinases. Curr. Med. Chem 2004, 11, 1563–1583. [Google Scholar]
  13. Bogoyevitch, M.A.; Fairlie, D. A new paradigm for protein kinase inhibition: Blocking phosphorylation without directly targeting ATP binding. Drug Discov. Today 2007, 12, 622–633. [Google Scholar]
  14. Liu, Y.; Gray, N.S. Rational design of inhibitors that bind to inactive kinase conformations. Nat. Chem. Biol 2006, 2, 358–364. [Google Scholar]
  15. Garuti, L.; Roberti, M.; Bottegoni, G. Non-ATP competitive protein kinase inhibitors. Curr. Med. Chem 2010, 17, 2804–2821. [Google Scholar]
  16. Okram, B.; Nagle, A.; Adrian, F.J.; Lee, C.; Ren, P.; Wang, X.; Sim, T.; Xie, Y.; Wang, X.; Xia, G.; et al. A general strategy for creating “inactive-conformation” Abl inhibitors. Chem. Biol 2006, 13, 779–786. [Google Scholar]
  17. Kirkland, L.O.; McInnes, C. Non-ATP competitive protein kinase inhibitors as anti-tumor therapeutics. Biochem. Pharmacol 2009, 77, 1561–1571. [Google Scholar]
  18. Kiselyov, A.; Balakin, K.V.; Tkachenko, S.E.; Savchuk, N.P. Recent progress in development of non-ATP competitive small-molecule inhibitors of protein kinesis. Mini Rev. Med. Chem 2006, 6, 711–717. [Google Scholar]
  19. Oguro, Y.; Miyamoto, N.; Okada, K.; Takagi, T.; Iwata, H.; Awazu, Y.; Miki, H.; Hori, A.; Kamiyama, K.; Imamura, S. Design, synthesis, and evaluation of 5-methyl-4-phenoxy-5H-pyrrolo-[ 3,2-d]pyrimidine derivatives: Novel VEGFR2 kinase inhibitors binding to inactive kinase conformation. Bioorg. Med. Chem 2010, 18, 7260–7273. [Google Scholar]
  20. Oguro, Y.; Miyamoto, N.; Takagi, T.; Okada, K.; Awazu, Y.; Miki, H.; Hori, A.; Kamiyama, K.; Imamura, S. N-Phenyl-N′-[4-(5H-pyrrolo[3,2-d]pyrimidin-4-yloxy) phenyl] ureas as novel inhibitors of VEGFR and FGFR kinases. Bioorg. Med. Chem 2010, 18, 7150–7163. [Google Scholar]
  21. Cramer, R.D., III; Patterson, D.E.; Bunce, J.D. Comparative molecular field analysis (CoMFA). 1. Effect of shape on binding of steroids to carrier proteins. J. Am. Chem. Soc 1988, 110, 5959–5967. [Google Scholar]
  22. Klebe, G.; Abraham, U.; Mietzner, T. Molecular similarity indices in a comparative analysis (CoMSIA) of drug molecules to correlate and predict their biological activity. J. Med. Chem 1994, 37, 4130–4146. [Google Scholar]
  23. Hasegawa, M.; Nishigaki, N.; Washio, Y.; Kano, K.; Harris, P.A.; Sato, H.; Mori, I.; West, R.I.; Shibahara, M.; Toyoda, H.; et al. Discovery of novel benzimidazoles as potent inhibitors of TIE-2 and VEGFR-2 tyrosine kinase receptors. J. Med. Chem 2007, 50, 4453–4470. [Google Scholar]
  24. Sybyl 7.3; Tripos Associates Inc.: St. Louis, MO, USA, 2009.
  25. Jain, A.N. Surflex: Fully automatic flexible molecular docking using a molecular similarity-based search engine. J. Med. Chem 2003, 46, 499–511. [Google Scholar]
  26. Jain, A.N. Surflex-Dock 2.1: Robzust performance from ligand energetic modeling, ring flexibility, and knowledge-based search. J. Comput. Aided Mol. Des 2007, 21, 281–306. [Google Scholar]
  27. Pan, J.; Liu, G.Y.; Cheng, J.; Chen, X.J.; Ju, X.L. CoMFA and molecular docking studies of benzoxazoles and benzothiazoles as CYP450 1A1 inhibitors. Eur. J. Med. Chem 2010, 45, 967–972. [Google Scholar]
  28. Sun, J.Y.; Cai, S.; Yan, N.; Mei, H. Docking and 3D-QSAR studies of influenza neuraminidase inhibitors using three-dimensional holographic vector of atomic interaction field analysis. Eur. J. Med. Chem 2010, 45, 1008–1014. [Google Scholar]
  29. Lan, P.; Chen, W.N.; Sun, P.H.; Chen, W.M. 3D-QSAR and molecular docking studies of azaindole derivatives as Aurora B kinase inhibitors. J. Mol. Model 2011, 17, 1191–1205. [Google Scholar]
  30. Ai, Y.; Wang, S.-T.; Sun, P.-H.; Song, F.-J. Molecular modeling studies of 4,5-Dihydro-1H-pyrazolo[ 4,3-h] quinazoline derivatives as potent CDK2/Cyclin a inhibitors using 3D-QSAR and docking. Int. J. Mol. Sci 2010, 11, 3705–3724. [Google Scholar]
  31. Ai, Y.; Wang, S.-T.; Sun, P.-H.; Song, F.-J. Combined 3D-QSAR modeling and molecular docking studies on Pyrrole-Indolin-2-ones as Aurora A Kinase inhibitors. Int. J. Mol. Sci 2011, 12, 1605–1624. [Google Scholar]
  32. Zhang, B.-D.; Li, Y.; Zhang, H.-X.; Ai, C.-Z. 3D-QSAR and molecular docking studies on derivatives of MK-0457, GSK1070916 and SNS-314 as inhibitors against Aurora B Kinase. Int. J. Mol. Sci 2010, 11, 4326–4347. [Google Scholar]
  33. Liu, J.-L.; Wang, F-F.; Ma, Z.; Wang, X.; Wang, Y.H. Structural determination of three different series of compounds as Hsp90 inhibitors using 3D-QSAR modeling, molecular docking and molecular dynamics methods. Int. J. Mol. Sci. 2011, 12, 946–970. [Google Scholar]
  34. Wu, X.-Y; Wu, S.-G; Chen, W.-H. Molecular docking and 3D-QSAR study on 4-(1H-indazol-4- yl)phenylamino and aminopyrazolopyridine urea derivatives as kinase insert domain receptor (KDR) inhibitors. J. Mol. Model 2011. [Google Scholar] [CrossRef]
  35. Sun, J.; Cai, S.; Mei, H.; Li, J.; Yan, N.; Wang, Y. Docking and 3D QSAR study of thiourea analogs as potent inhibitors of influenza virus neuraminidase. J. Mol. Model 2010, 16, 1809–1818. [Google Scholar]
Figure 1. Superimposition of compounds in the training and test set.
Figure 1. Superimposition of compounds in the training and test set.
Ijms 13 02387f1
Figure 2. Compound 53 from 2OH4 and atom-numbered compound 20.
Figure 2. Compound 53 from 2OH4 and atom-numbered compound 20.
Ijms 13 02387f2
Figure 3. Binding conformations of (a) co-crystallized (magenta) and re-docked (green) 53 and (b) docked compound 20 at the active site of kinase insert domain receptor (KDR) in the inactive conserved triad aspartate-phenylalanine-glycine (DFG)-out conformation. Key residues are displayed and hydrogen bonds are displayed in dotted lines.
Figure 3. Binding conformations of (a) co-crystallized (magenta) and re-docked (green) 53 and (b) docked compound 20 at the active site of kinase insert domain receptor (KDR) in the inactive conserved triad aspartate-phenylalanine-glycine (DFG)-out conformation. Key residues are displayed and hydrogen bonds are displayed in dotted lines.
Ijms 13 02387f3
Figure 4. Graphs of the experimental versus predicted pIC50 values of the training (■) and test (▴) compounds from the CoMFA and CoMSIA models.
Figure 4. Graphs of the experimental versus predicted pIC50 values of the training (■) and test (▴) compounds from the CoMFA and CoMSIA models.
Ijms 13 02387f4
Figure 5. CoMFA STDEV*COEFF contour maps. (a) Favorable (green) and unfavorable (yellow) steric fields. (b) Electropositive (blue) and electronegative (red) fields. Compound 20 was overlaid in each map.
Figure 5. CoMFA STDEV*COEFF contour maps. (a) Favorable (green) and unfavorable (yellow) steric fields. (b) Electropositive (blue) and electronegative (red) fields. Compound 20 was overlaid in each map.
Ijms 13 02387f5
Figure 6. STDEV*COEFF contour maps. (a) Favorable (green) and unfavorable (yellow) steric fields. (b) Electropositive (blue) and electronegative (red) fields. (c) Favorable (yellow) and unfavorable (gray) hydrophobic fields. (d) Favorable (cyan) and unfavorable (purple) hydrogen bond donor fields. (e) Favorable (magenta) and unfavorable (red) hydrogen bond acceptor fields. Compound 20 was overlaid in each plot.
Figure 6. STDEV*COEFF contour maps. (a) Favorable (green) and unfavorable (yellow) steric fields. (b) Electropositive (blue) and electronegative (red) fields. (c) Favorable (yellow) and unfavorable (gray) hydrophobic fields. (d) Favorable (cyan) and unfavorable (purple) hydrogen bond donor fields. (e) Favorable (magenta) and unfavorable (red) hydrogen bond acceptor fields. Compound 20 was overlaid in each plot.
Ijms 13 02387f6
Table 1. Structures of compounds 152.
Table 1. Structures of compounds 152.
StructureCompoundSubstituent

R1R2
Ijms 13 02387f71NHCONHPhH
2 *NHCOPhH
3CONHPhH
4HNHCOPh
5HCONHPh
6 *HNHCOCH2Ph
7HCH2CONHPh
8HNHCSNHPh
9HIjms 13 02387f8
10HNHCONHMe
11HNHCONHPr
Ijms 13 02387f9XR
12SH
13NHH
14N(Me)3-CF3
Ijms 13 02387f10
Table 1. Structures of compounds 152.
StructureCompoundSubstituent

RR1R2R3R4
Ijms 13 02387f1116 *PhHHMeH
172-ClPhHHMeH
183-ClPhHHMeH
194-ClPhHHMeH
20 *3-CF3PhHHMeH
213-BrPhHHMeH
223-FPhHHMeH
233-MePhHHMeH
24 *3-CF3PhClHMeH
253-CF3PhHClMeH
263-CF3PhOMeHMeH
273-CF3PhFHMeH
283-CF3PhClH(CH2)2OHH
29 *3-CF3PhClH(CH2)2OMeH
303-CF3PhClH(CH2)2O(CH2)2OHH
313-CF3PhClHMeCH2OH
323-CF3PhClHMeCH2OMe
333-CF3PhClHMeCMe2OH
342-C5H4NClHMeH
353-C5H4NClHMeH
364-C5H4NClHMeH
37Ijms 13 02387f12ClHMeH
38Ijms 13 02387f13ClHMeH
39Ijms 13 02387f14ClHMeH
40Ijms 13 02387f15ClHMeH
41Ijms 13 02387f16ClHMeH
42Ijms 13 02387f17ClHMeH
43Ijms 13 02387f18ClHMeH
R1R2
Ijms 13 02387f1944OMeH
45Ijms 13 02387f20H
46Ijms 13 02387f21H
47Ijms 13 02387f22H
48 *Ijms 13 02387f23H
49 *HOMe
50HIjms 13 02387f24
51HIjms 13 02387f25
52HIjms 13 02387f26
*Test set.
Table 2. The experimental pIC50, predicted pIC50 and their residuals of compounds 152.
Table 2. The experimental pIC50, predicted pIC50 and their residuals of compounds 152.
CompoundExp pIC50COMFACOMSIA

PredResidPredResid
15.7455.1940.5515.834−0.089
2 *6.0275.5630.4646.327−0.300
36.0605.6210.4396.246−0.185
45.0274.7880.2394.8210.206
55.0005.192−0.1924.9800.020
6 *5.0004.8080.1925.074−0.074
75.4815.698−0.2165.2850.197
85.1746.042−0.8685.1290.045
95.3776.029−0.6525.1560.221
105.2014.8520.3495.1080.093
115.4095.3070.1025.420−0.011
126.9597.074−0.1166.8950.063
135.8545.7380.1165.987−0.133
145.0005.233−0.2335.341−0.341
158.1497.7290.4207.8100.339
16 *7.4817.1030.3797.1890.293
176.4696.863−0.3946.927−0.458
188.3877.9590.4288.1870.200
197.4957.2010.2937.503−0.008
20 *8.2768.2090.0668.341−0.066
218.3577.9690.3878.710−0.354
227.7217.4000.3217.3350.387
238.5698.3230.2468.1530.415
24 *8.4328.2550.1778.3240.108
257.5237.966−0.4437.987−0.464
267.8548.056−0.2028.171−0.317
278.2088.1760.0328.357−0.149
288.3988.812−0.4148.3860.012
29 *8.0468.717−0.6718.262−0.216
308.5698.977−0.4088.3080.261
318.0137.9870.0277.8510.162
327.8547.8350.0197.7740.080
337.7457.5540.1917.7010.044
346.6787.047−0.3697.042−0.364
356.6387.354−0.7167.090−0.451
366.1876.691−0.5046.342−0.155
378.5698.3070.2628.2650.303
388.5388.3570.1818.2880.249
398.4207.9830.4388.0840.336
408.3987.6890.7098.2730.125
418.1198.357−0.2388.282−0.163
428.0518.079−0.0287.9930.058
437.7218.019−0.2978.392−0.670
448.0278.0130.0137.8750.152
457.2157.1650.0507.2030.011
467.7217.6340.0887.805−0.084
478.0327.7210.3108.072−0.041
48 *7.9217.7060.2157.2260.695
49 *7.6997.2340.4656.9110.788
507.2227.1350.0877.0730.149
517.6387.689−0.0517.5530.086
527.9597.9140.0447.7350.223
*Test set.
Table 3. Statistical parameters for the CoMFA and CoMSIA models.
Table 3. Statistical parameters for the CoMFA and CoMSIA models.
Nq2r2SEEFr2predField contribution

SEHDA
CoMFA40.5420.9120.376100.4620.9130.5250.475---
CoMSIA50.5520.9550.272161.2450.8970.1840.2280.3430.0630.182
q2: Cross-validated correlation coefficient; r2: non-cross-validated correlation coefficient; r2pred: predictive correlation coefficient; SEE: standard error of estimate; F: Fischer ratio; N: optimal number of principal components; S: steric field; E: electrostatic field; H: hydrophobic field; D: hydrogen bond donor field; A: hydrogen bond acceptor field.
Table 4. Structures and predicted pIC50 values of newly designed derivatives.
Table 4. Structures and predicted pIC50 values of newly designed derivatives.
Ijms 13 02387f27

CompoundRR1Predicted pIC50

CoMFACoMSIA
D1ClIjms 13 02387f288.1628.052
D2ClIjms 13 02387f298.4638.049
D3ClIjms 13 02387f307.9727.876
D4ClIjms 13 02387f318.3739.102
D5ClIjms 13 02387f328.1248.696
D6ClIjms 13 02387f337.7978.405
D7H3-EtPh8.2168.507
D8H3-CH2CF3Ph8.3788.891
D9H3,4-2CH3Ph8.1938.609
D10H3-CF3,4-CH3Ph8.2348.837
D11Cl3-CH3Ph7.9938.096
D12Cl3,4-2CH3Ph8.1708.582
D13Cl3-CF3,4-CH3Ph8.2588.814
D14Cl3-ClPh7.8988.245
D15Cl3-Cl,4-CH3Ph8.1918.731

Share and Cite

MDPI and ACS Style

Wu, X.-Y.; Chen, W.-H.; Wu, S.-G.; Tian, Y.-X.; Zhang, J.-J. Pyrrolo[3,2-d]pyrimidine Derivatives as Type II Kinase Insert Domain Receptor (KDR) Inhibitors: CoMFA and CoMSIA Studies. Int. J. Mol. Sci. 2012, 13, 2387-2404. https://doi.org/10.3390/ijms13022387

AMA Style

Wu X-Y, Chen W-H, Wu S-G, Tian Y-X, Zhang J-J. Pyrrolo[3,2-d]pyrimidine Derivatives as Type II Kinase Insert Domain Receptor (KDR) Inhibitors: CoMFA and CoMSIA Studies. International Journal of Molecular Sciences. 2012; 13(2):2387-2404. https://doi.org/10.3390/ijms13022387

Chicago/Turabian Style

Wu, Xiao-Yun, Wen-Hua Chen, Shu-Guang Wu, Yuan-Xin Tian, and Jia-Jie Zhang. 2012. "Pyrrolo[3,2-d]pyrimidine Derivatives as Type II Kinase Insert Domain Receptor (KDR) Inhibitors: CoMFA and CoMSIA Studies" International Journal of Molecular Sciences 13, no. 2: 2387-2404. https://doi.org/10.3390/ijms13022387

APA Style

Wu, X. -Y., Chen, W. -H., Wu, S. -G., Tian, Y. -X., & Zhang, J. -J. (2012). Pyrrolo[3,2-d]pyrimidine Derivatives as Type II Kinase Insert Domain Receptor (KDR) Inhibitors: CoMFA and CoMSIA Studies. International Journal of Molecular Sciences, 13(2), 2387-2404. https://doi.org/10.3390/ijms13022387

Article Metrics

Back to TopTop