- Research article
- Open Access
2D and 3D similarity landscape analysis identifies PARP as a novel off-target for the drug Vatalanib
© Gohlke et al. 2015
- Received: 20 May 2015
- Accepted: 8 September 2015
- Published: 24 September 2015
Searching for two-dimensional (2D) structural similarities is a useful tool to identify new active compounds in drug-discovery programs. However, as 2D similarity measures neglect important structural and functional features, similarity by 2D might be underestimated. In the present study, we used combined 2D and three-dimensional (3D) similarity comparisons to reveal possible new functions and/or side-effects of known bioactive compounds.
We utilised more than 10,000 compounds from the SuperTarget database with known inhibition values for twelve different anti-cancer targets. We performed all-against-all comparisons resulting in 2D similarity landscapes. Among the regions with low 2D similarity scores are inhibitors of vascular endothelial growth factor receptor (VEGFR) and inhibitors of poly ADP-ribose polymerase (PARP). To demonstrate that 3D landscape comparison can identify similarities, which are untraceable in 2D similarity comparisons, we analysed this region in more detail. This 3D analysis showed the unexpected structural similarity between inhibitors of VEGFR and inhibitors of PARP. Among the VEGFR inhibitors that show similarities to PARP inhibitors was Vatalanib, an oral “multi-targeted” small molecule protein kinase inhibitor being studied in phase-III clinical trials in cancer therapy. An in silico docking simulation and an in vitro HT universal colorimetric PARP assay confirmed that the VEGFR inhibitor Vatalanib exhibits off-target activity as a PARP inhibitor, broadening its mode of action.
In contrast to the 2D-similarity search, the 3D-similarity landscape comparison identifies new functions and side effects of the known VEGFR inhibitor Vatalanib.
- Drug action
- Drug discovery
- Vascular endothelial growth factor (VEGF)
- 3D similarity landscapes
Drugs often not only interact with their intended target but also with so-called off-targets, thereby causing side-effects . Prediction of side-effects is still a big challenge during drug design and studies have shown the potential of computational methods for target and off-target analysis [2, 3]. These studies deal with pathway- and network-based approaches combined with the chemical structure of small molecule compounds regarding their binding site at the target protein [4–6].
To identify compound similarities it is important to take a detailed look at 2D- and 3D-similarities . 2D structural similarity algorithms were generated to predict and create a drug-target adverse drug reactions (ADR) network [3, 8]. These 2D-fingerprints represent the structure and properties of small molecules by a bit or integer string. Although several methods exist to measure the similarity between 2D fingerprints, the Tanimoto coefficient has been proven to be reliable [9, 10]. However, several problems can occur while working with fingerprints: size of compounds as well as functional groups or side chains have an impact on the similarity calculations. Hence, functional and structural features of compounds can be neglected. These problems can be overcome by using 3D similarity search methods.
Non-commercial drug- or target-related databases, which have been established in the last decade, can be used for 2D and 3D comparisons. Millions of compounds can be found in databases like ChEMBL  or PubChem  and their availability can be verified via the ZINC database .
We recently established our SuperTarget database, which was developed with the intention to accentuate drug–target interactions and to provide references to other resources for more elaborate analysis . The SuperTarget database contains a core dataset of about 330,000 drug-target interactions, of which about 310,000 interactions have calculated binding affinity data  and were used to compare 2D and 3D structures of promising anticancer drugs.
Among these drugs are inhibitors of the poly ADP-ribose polymerase (PARP). PARP binds to single-strand DNA breaks and plays a critical role in cell recovery from DNA damage. PARP inhibitors show activities not only in cancer therapy but are also being evaluated for the treatment of stroke, myocardial infarction and other diseases. Additional promising anticancer drugs, which can be found in the SuperTarget database, are inhibitors of the vascular endothelial growth factor receptor (VEGFR). The approved VEGFR inhibitor Vatalanib (PTK787 or PTK/ZK) is currently studied in several phases of clinical trials for different cancer therapies [16–18]. Vatalanib is an oral “multi-targeted” small molecule protein kinase inhibitor that binds to the intracellular kinase domain of all VEGF receptor subtypes, thereby inhibiting angiogenesis . In addition, it binds to c-KIT and platelet-derived growth factor receptor (PDGFR) but with lower affinity.
While applying a 3D similarity landscape analysis on inhibitors for different cancer targets by using the SuperTarget database, we found unexpected similarities between PARP and VEGFR inhibitors, which could not be detected by 2D similarity searches. As a proof of concept of our similarity landscape analysis, both in silico and in vitro assays confirmed Vatalanib’s off-target activity as a PARP inhibitor. In this paper we provide a combined approach of 2D and 3D similarity landscapes for target and off-target analysis, which can be applied to a larger number of targeted anti-cancer therapeutics.
where AB is the number of bits set to one in both molecules, A is the number of bits set to one in molecule A and B is the number of bits set to one in molecule B.
Another method to calculate the Tanimoto coefficient are the extended-connectivity fingerprints (ECFP) . These are used to cover the calculated 2D-similarity by OpenBabel fingerprints, which belong to the class of radial fingerprints and are based on the Morgan algorithm . To calculate the extended-connectivity fingerprints the cheminformatics toolkit of ChemAxon was used (JChem compr (18.104.22.168), 201n (2014), ChemAxon (http://www.chemaxon.com)).
Every molecule’s conformation was compared with each conformation of the second molecule, resulting in up to 2,500 separately calculated rmsd values. Here, only the smallest rmsd value, i.e. the best superposition of the compounds, was stored.
Ligand Docking - The docking study was performed by using LibDock, a high-throughput docking algorithm for library design and library prioritisation. This docking program was provided by Accelrys Discovery Studio (http://accelrys.com). The algorithm positioned ligands in the protein’s active site based on polar and non-polar interaction sites.
MCF-7 cell lines - Breast cancer cell lines MCF-7 were cultured in RPMI-1640 medium supplemented with 10 % inactivated FBS, 100 U/ml penicillin and 0.1 mg/ml streptomycin. Cells were cultured at 37 °C with 5 % CO2 in a fully humidified atmosphere.
IC 50 values of PARP inhibitors - For the determination of IC 50 values of Vatalanib and Compound 1 we used the HT universal colorimetric PARP assay kit with histone-coated strip wells (Trevigen, USA). Absorbance was measured in a Sunrise microplate reader (Tecan, Switzerland) at 450 nm.
γH2AX foci analysis - For immunofluorescence microscopic analyses, MCF-7 cells were grown on coverslips. 24 h post treatment with 0 (control), 1, 10, and 100 μM Compound 1 or Vatalanib, cells were washed in PBS, fixed in 3 % paraformaldehyde/PBS (15 min), permeabilised with 0.5 % Triton-X 100/PBS (2 min) and blocked in 5 % fetal bovine serum for 60 min at room temperature. After incubation with anti-phospho-Histone H2A.X (Ser139) clone JBW301 (mouse monoclonal IgG from Millipore, Billerica, MA, USA) overnight at 4 °C, cells were incubated with Alexa Fluor 488-labelled chicken anti-mouse IgG secondary antibody (Molecular Probes, The Netherlands) for 2 h at room temperature and counterstained with DAPI. Images of γH2AX foci and DAPI-labelled nuclei were acquired with a fluorescent microscope (BX50; Olympus, Germany) equipped with a 40×/0.75 objective lens (UPlanFL; Olympus, Germany) and a camera (micropublisher 5.0 RTV; QImaging, Canada) with Openlab software (Perkin Elmer, Germany).
Because 2D-similarity analyses often neglect important structural and functional features, we expanded our comparison to 3D-superpositions measured by the root-mean-square deviation (rmsd) based on the Kabsch algorithm. Although current computers calculate 3D-comparisons of compounds relatively fast, it would still take months to compare all inhibitors with each other. We therefore focused on the 3D-structural comparison of VEGFR and PARP inhibitors, which showed only little structural similarity in the 2D-similarity analyses. For this comparison up to 50 conformers were calculated by using Accelrys Discovery studio 3.5 (Accelrys Software Inc., Discovery Studio Modeling Environment, Release 3.5, San Diego: Accelrys Software Inc., 2012). To create diverse ligand conformations, the ‘fast’ search method was used to generate multiple low-energy conformations. The rmsd calculations are based on overlaying the anchor points of both conformers and to rotate at single bonds (degrees of freedom) to minimise the rmsd. We arbitrarily chose a 5 % quantile (rmsd of 0.215) to evaluate the similarity of related 3D structures and analysed only compounds with both low rmsd (high similarity) and low Tanimoto scores (low 2D-similarity).
The 3D screening identifies Vatalanib as a potential inhibitor of PARP (rmsd: 0.194) whereas OpenBabel fingerprints calculated a low Tanimoto score of 0.4 and therefore failed to recognise the similarity of these two compounds. To analyse in more detail if other 2D fingerprints uncover the similarity between Vatalanib and Compound 1, Tanimoto scores were likewise calculated by extended connectivity fingerprints (ECFP), which result in a score of 0.32 for Compound 1 and Vatalanib (data matrix not shown). In addition, MACCS and FP3 fingerprints also computed low Tanimoto scores of 0.37 and 0.33, respectively. This confirmed that 2D analyses are unsuccessful in identifying the similarity between Vatalanib and Compound 1.
To confirm inhibition of PARP and to rule out in vivo vs. in vitro activity discrepancies, i.e. differences of PARP expression in cells vs. an isolated PARP enzyme, we analysed accumulation of DNA damage in a human breast cancer cell line upon treatment with Compound 1 and Vatalanib. Because of PARP's involvement in DNA strand break repair, its inhibition has been proposed to lead to double-strand break (DSB) formation . These DSBs induce phosphorylation of histone H2AX on Ser-139 at sites flanking the breakage [27, 28]. Therefore, we analysed whether treated cells accumulate phosphorylated H2AX, denoted as γH2AX, which provides a common marker for DNA damage in vitro [26, 29].
2D fingerprints similarity search methods are widely used approaches in the discovery of novel molecules with high affinity to specific targets and, despite the fact that molecules are active in three dimensions, surprisingly powerful . In this study we used the freely available software packages OpenBabel and ChemAxon to analyse the 2D-similarity of about 10,000 inhibitors against twelve promising anti-cancer targets. Among these VEGFR and PARP inhibitors showed only little structural similarity, however, similarity by 2D might be underestimated .
Accurate target prediction can be achieved by combining different measures of chemical similarity based on both chemical structure and molecular shape . Furthermore it has been shown recently that the combination of a 2D similarity search and a 3D shape/flexibility-based similarity search led to an increased hit rate . Therefore 3D-similarity of VEGFR and PARP inhibitors was then analysed by a proprietary 3D-superpostion algorithm, which produces reproducible results because of pre-calculated conformers for every compound. The 3D-similarity method in combination with 2D-similarity comparison performs quite well by applying a 5 % quantile threshold (corresponding to an rmsd-value of 0.215) for early discovery detection. This is possible as our data follows a normal distribution. Nevertheless, using a threshold means missing out on compounds with larger rmsd-values (>0.215), which could have the same inhibitory function. By taking 50 conformers into account to simulate the flexibility and to cover the conformational space, 2,500 superpositions of the two compounds were calculated and the best superposition with a minimal rmsd-value was taken for the similarity measurement. This makes the method robust with respect to conformation changes. Although the landscape of the 3D-screening shows an overall reduced similarity of inhibitors compared to the 2D-landscape, selected PARP inhibitors display high similarity when compared to the VEGFR inhibitors, confirming that it is important to take 3D in addition to 2D similarity into account to increase the hit rate. Among the inhibitors with similarity we found the anticancer agent Vatalanib, and Compound 1, a similarity not identified by different 2D fingerprint algorithms. It might be that the 2D algorithms failed to identify this similarity, because these fingerprints are based on atom labels whereas rmsd does not consider these labels. Structural comparisons as performed by the 3D algorithm used here might perform better generally, as they compensate for fragments of the molecule and their connections of different atom types. Still, there might be clusters where 2D similarity algorithms might be faster and better.
Compound 1 was identified as a direct PARP1 inhibitor in a yeast screen assay with an EC50 of approximately 60 μM . By using the HT universal colorimetric PARP assay kit we calculated an IC50 value of approximately 3,000 μM. The differences in the IC50 values of Compound 1 can be explained by differences in the assays used as well as by differences of PARP expressed in yeast compared to an isolated PARP enzyme used in our assay. According to the colorimetric PARP assay, Vatalanib, with an IC50 value of 200 μM is fifteen-fold more effective than Compound 1 in inhibiting PARP. The in silico docking simulation indicates that Vatalanib’s additional chloride atom, which is missing in Compound 1, is in close proximity to the arginine 204 at the bottom of the binding site of PARP. A halogen bond between the nitrogen of arginine 204 might be formed, stabilising the binding position, which results in a more effective drug target inhibition. Both Vatalanib and Compound 1 have higher LibDock docking scores than 3-AB, which might be attributed to the relative small molecule size of 3-AB.
Despite better LibDock docking scores, both Vatalanib and Compound 1, are less potent than the prototype PARP inhibitor 3-AB (3-aminobenzamide), which has an IC50 of about 30 μM [34, 35]. Despite these comparatively high IC50 values of Compound 1 and Vatalanib, 10 μM of either of the drugs was able to induce γH2AX foci formation in human breast cancer cells, demonstrating DNA damage and PARP inhibition. In regard to this activity, these results again point to an in vivo vs. in vitro discrepancy with a higher bioactivity in cells compared to the enzyme assay. Due to the relative high IC50 of Vatalanib, Vatalanib might be of interest as a new PARP inhibitor or for drug design. In addition a putative positive side effect of Vatalanib against PARP might be important for the use of Vatalanib as a chemotherapeutical. Vatalanib doses up to 1,000 mg twice-a-day are well tolerated reaching plasma concentration in patients in the μM range. Accordingly, the 10 μM of Vatalanib, able to induce γH2AX foci, is in the rage of cmax plasma concentrations achieved in patients [36–39].
Vatalanib (PTK787/ZK222584) was initially described as a selective tyrosine kinase inhibitor (TKI) of VEGFR1-3. TKIs commonly have additional activity against other tyrosine kinases. Likewise Vatalanib, which at higher concentrations also inhibits other protein tyrosine kinases of the same family, such as the platelet-derived growth factor receptor beta tyrosine kinase and the c-Kit protein tyrosine kinase . Interestingly, activity across other classes of drug targets have also been documented for Vatalanib. It has been shown that Vatalanib significantly inhibits aromatase and thus might cross-inhibit two important classes of targets in breast cancer [40, 41]. This “multi-targeting” activity, which might also include PARP as a target, could potentially contribute to the antitumor effect of Vatalanib and indicates that a drug’s efficacy often might not only be based on the inhibition of one target but of multiple targets [42, 43].
3D similarity landscape comparison, as shown in this study, has the potential to identify new targets of known drugs. As a proof of principle, we identified Vatalanib’s additional ability to target PARP, which was demonstrated in vitro and in vivo. Thus, combined 2D and 3D similarity landscape comparison analysis can identify new functions and/or side effects of known bioactive compounds that are untraceable with 2D similarity searching alone.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Medina-Franco JL, Giulianotti MA, Welmaker GS, Houghten RA. Shifting from the single to the multitarget paradigm in drug discovery. Drug Discov Today. 2013;18(9-10):495–501.View ArticlePubMedPubMed CentralGoogle Scholar
- Gregori-Puigjane E, Mestres J. A ligand-based approach to mining the chemogenomic space of drugs. Com Chem High t Scr. 2008;11(8):669–76.Google Scholar
- Lagunin A, Stepanchikova A, Filimonov D, Poroikov V. PASS: prediction of activity spectra for biologically active substances. Bioinformatics. 2000;16(8):747–8.View ArticlePubMedGoogle Scholar
- von Eichborn J, Murgueitio MS, Dunkel M, Koerner S, Bourne PE, Preissner R. PROMISCUOUS: a database for network-based drug-repositioning. Nucleic Acids Res. 2011;39(Database issue):D1060–1066.View ArticleGoogle Scholar
- Brouwers L, Iskar M, Zeller G, van Noort V, Bork P. Network neighbors of drug targets contribute to drug side-effect similarity. PLoS One. 2011;6(7):e22187.View ArticlePubMedPubMed CentralGoogle Scholar
- Yamanishi Y, Pauwels E, Kotera M. Drug side-effect prediction based on the integration of chemical and biological spaces. J Chem Inf Model. 2012;52(12):3284–92.View ArticlePubMedGoogle Scholar
- Willett P, Barnard JM, Downs GM. Chemical similarity searching. J Chem Inf Comp Sci. 1998;38(6):983–96.View ArticleGoogle Scholar
- Godden JW, Xue L, Stahura FL, Bajorath J: Searching for molecules with similar biological activity: analysis by fingerprint profiling. Pac Symp Biocomput 2000;8:566-575Google Scholar
- Holliday JD, Salim N, Whittle M, Willett P. Analysis and display of the size dependence of chemical similarity coefficients. J Chem Inf Comp Sci. 2003;43(3):819–28.View ArticleGoogle Scholar
- Rademacher C, Paulson JC. Glycan fingerprints: calculating diversity in glycan libraries. ACS Chem Biol. 2012;7(5):829–34.View ArticlePubMedPubMed CentralGoogle Scholar
- Gaulton A, Bellis LJ, Bento AP, Chambers J, Davies M, Hersey A, et al. ChEMBL: a large-scale bioactivity database for drug discovery. Nucleic Acids Res. 2012;40(Database issue):D1100–1107.View ArticlePubMedGoogle Scholar
- Wang Y, Xiao J, Suzek TO, Zhang J, Wang J, Bryant SH. PubChem: a public information system for analyzing bioactivities of small molecules. Nucleic Acids Res. 2009;37(Web Server issue):W623–633.View ArticlePubMedPubMed CentralGoogle Scholar
- Irwin JJ, Shoichet BK. ZINC--a free database of commercially available compounds for virtual screening. J Chem Inf Model. 2005;45(1):177–82.View ArticlePubMedPubMed CentralGoogle Scholar
- Gunther S, Kuhn M, Dunkel M, Campillos M, Senger C, Petsalaki E, et al. SuperTarget and Matador: resources for exploring drug-target relationships. Nucleic Acids Res. 2008;36(Database issue):D919–922.PubMedGoogle Scholar
- Hecker N, Ahmed J, von Eichborn J, Dunkel M, Macha K, Eckert A, et al. SuperTarget goes quantitative: update on drug-target interactions. Nucleic Acids Res. 2012;40(Database issue):D1113–1117.View ArticlePubMedGoogle Scholar
- Los M, Roodhart JM, Voest EE. Target practice: lessons from phase III trials with bevacizumab and vatalanib in the treatment of advanced colorectal cancer. Oncologist. 2007;12(4):443–50.View ArticlePubMedGoogle Scholar
- Raizer JJ, Grimm SA, Rademaker A, Chandler JP, Muro K, Helenowski I, et al. A phase II trial of PTK787/ZK 222584 in recurrent or progressive radiation and surgery refractory meningiomas. J Neurooncol. 2014;117(1):93–101.View ArticlePubMedGoogle Scholar
- Bitting RL, Healy P, Creel PA, Turnbull J, Morris K, Wood SY, et al. A phase Ib study of combined VEGFR and mTOR inhibition with vatalanib and everolimus in patients with advanced renal cell carcinoma. Clin Genitourin Cancer. 2014;12(4):241–50.View ArticlePubMedGoogle Scholar
- Wood JM, Bold G, Buchdunger E, Cozens R, Ferrari S, Frei J, et al. PTK787/ZK 222584, a novel and potent inhibitor of vascular endothelial growth factor receptor tyrosine kinases, impairs vascular endothelial growth factor-induced responses and tumor growth after oral administration. Cancer Res. 2000;60(8):2178–89.PubMedGoogle Scholar
- Willett P, Winterman V. A comparison of some measures for the determination of intermolecular structural similarity measures of intermolecular structural similarity. Quant Struct-Act Rel. 1986;5(1):18–25.View ArticleGoogle Scholar
- Rogers D, Hahn M. Extended-connectivity fingerprints. J Chem Inf Model. 2010;50(5):742–54.View ArticlePubMedGoogle Scholar
- Morgan HL. Generation of a unique machine description for chemical structures-a technique developed at chemical abstracts service. J Chem Doc. 1965;5(2):107–13.View ArticleGoogle Scholar
- Kabsch W. Solution for best rotation to relate 2 sets of vectors. Acta Crystallogr A. 1976;32(Sep1):922–3.View ArticleGoogle Scholar
- Perkins E, Sun D, Nguyen A, Tulac S, Francesco M, Tavana H, et al. Novel inhibitors of poly(ADP-ribose) polymerase/PARP1 and PARP2 identified using a cell-based screen in yeast. Cancer Res. 2001;61(10):4175–83.PubMedGoogle Scholar
- Diller DJ, Merz Jr KM. High throughput docking for library design and library prioritization. Proteins. 2001;43(2):113–24.View ArticlePubMedGoogle Scholar
- Redon CE, Nakamura AJ, Zhang YW, Ji JJ, Bonner WM, Kinders RJ, et al. Histone gammaH2AX and poly(ADP-ribose) as clinical pharmacodynamic biomarkers. Clin Cancer Res. 2010;16(18):4532–42.View ArticlePubMedPubMed CentralGoogle Scholar
- West MH, Bonner WM. Histone 2A, a heteromorphous family of eight protein species. Biochemistry. 1980;19(14):3238–45.View ArticlePubMedGoogle Scholar
- Rogakou EP, Pilch DR, Orr AH, Ivanova VS, Bonner WM. DNA double-stranded breaks induce histone H2AX phosphorylation on serine 139. J Biol Chem. 1998;273(10):5858–68.View ArticlePubMedGoogle Scholar
- Huang X, Darzynkiewicz Z. Cytometric assessment of histone H2AX phosphorylation: a reporter of DNA damage. Methods Mol Biol. 2006;314:73–80.View ArticlePubMedPubMed CentralGoogle Scholar
- Bajorath J. Selected concepts and investigations in compound classification, molecular descriptor analysis, and virtual screening. J Chem Inf Comp Sci. 2001;41(2):233–45.View ArticleGoogle Scholar
- Thimm M, Goede A, Hougardy S, Preissner R. Comparison of 2D similarity and 3D superposition. Application to searching a conformational drug database. J Chem Inf Comp Sci. 2004;44(5):1816–22.View ArticleGoogle Scholar
- Gfeller D, Michielin O, Zoete V. Shaping the interaction landscape of bioactive molecules. Bioinformatics. 2013;29(23):3073–9.View ArticlePubMedGoogle Scholar
- Dobi K, Hajdu I, Flachner B, Fabo G, Szaszko M, Bognar M, et al. Combination of 2D/3D ligand-based similarity search in rapid virtual screening from multimillion compound repositories. Selection and biological evaluation of potential PDE4 and PDE5 inhibitors. Molecules. 2014;19(6):7008–39.View ArticlePubMedGoogle Scholar
- Costantino G, Macchiarulo A, Camaioni E, Pellicciari R. Modeling of poly(ADP-ribose)polymerase (PARP) inhibitors. Docking of ligands and quantitative structure-activity relationship analysis. J Med Chem. 2001;44(23):3786–94.View ArticlePubMedGoogle Scholar
- Dillon KJ, Smith GC, Martin NM. A FlashPlate assay for the identification of PARP-1 inhibitors. J Biomol Screen. 2003;8(3):347–52.View ArticlePubMedGoogle Scholar
- Reardon DA, Egorin MJ, Desjardins A, Vredenburgh JJ, Beumer JH, Lagattuta TF, et al. Phase I pharmacokinetic study of the vascular endothelial growth factor receptor tyrosine kinase inhibitor vatalanib (PTK787) plus imatinib and hydroxyurea for malignant glioma. Cancer. 2009;115(10):2188–98.View ArticlePubMedPubMed CentralGoogle Scholar
- Thomas AL, Morgan B, Horsfield MA, Higginson A, Kay A, Lee L, et al. Phase I study of the safety, tolerability, pharmacokinetics, and pharmacodynamics of PTK787/ZK 222584 administered twice daily in patients with advanced cancer. J Clin Oncol. 2005;23(18):4162–71.View ArticlePubMedGoogle Scholar
- Morgan B, Thomas AL, Drevs J, Hennig J, Buchert M, Jivan A, et al. Dynamic contrast-enhanced magnetic resonance imaging as a biomarker for the pharmacological response of PTK787/ZK 222584, an inhibitor of the vascular endothelial growth factor receptor tyrosine kinases, in patients with advanced colorectal cancer and liver metastases: results from two phase I studies. J Clin Oncol. 2003;21(21):3955–64.View ArticlePubMedGoogle Scholar
- Jost LM, Gschwind HP, Jalava T, Wang Y, Guenther C, Souppart C, et al. Metabolism and disposition of vatalanib (PTK787/ZK-222584) in cancer patients. Drug Metab Dispos. 2006;34(11):1817–28.View ArticlePubMedGoogle Scholar
- Banerjee S, Zvelebil M, Furet P, Mueller-Vieira U, Evans DB, Dowsett M, et al. The vascular endothelial growth factor receptor inhibitor PTK787/ZK222584 inhibits aromatase. Cancer Res. 2009;69(11):4716–23.View ArticlePubMedGoogle Scholar
- Banerjee S, A'Hern R, Detre S, Littlewood-Evans AJ, Evans DB, Dowsett M, et al. Biological evidence for dual antiangiogenic-antiaromatase activity of the VEGFR inhibitor PTK787/ZK222584 in vivo. Clin Cancer Res. 2010;16(16):4178–87.View ArticlePubMedGoogle Scholar
- Hopkins AL, Mason JS, Overington JP. Can we rationally design promiscuous drugs? Curr Opin Struct Biol. 2006;16(1):127–36.View ArticlePubMedGoogle Scholar
- Petrelli A, Giordano S. From single- to multi-target drugs in cancer therapy: when aspecificity becomes an advantage. Curr Med Chem. 2008;15(5):422–32.View ArticlePubMedGoogle Scholar