John S. Tokarski

The Structure of Dasatinib (BMS-354825) Bound to Activated ABL Kinase Domain Elucidates Its Inhibitory Activity against Imatinib-Resistant ABL Mutants

Chronic myeloid leukemia (CML) is caused by the constitutively activated tyrosine kinase breakpoint cluster (BCR)-ABL. Current frontline therapy for CML is imatinib, an inhibitor of BCR-ABL. Although imatinib has a high rate of clinical success in early phase CML, treatment resistance is problematic, particularly in later stages of the disease, and is frequently mediated by mutations in BCR-ABL. Dasatinib (BMS-354825) is a multitargeted tyrosine kinase inhibitor that targets oncogenic pathways and is a more potent inhibitor than imatinib against wild-type BCR-ABL. It has also shown preclinical activity against all but one of the imatinib-resistant BCR-ABL mutants tested to date. Analysis of the crystal structure of dasatinib-bound ABL kinase suggests that the increased binding affinity of dasatinib over imatinib is at least partially due to its ability to recognize multiple states of BCR-ABL. The structure also provides an explanation for the activity of dasatinib against imatinib-resistant BCR-ABL mutants.

Discovery of N-(4-(2-Amino-3-chloropyridin-4-yloxy)-3-fluorophenyl)-4-ethoxy-1-(4-fluorophenyl)-2-oxo-1,2-dihydropyridine-3-carboxamide (BMS-777607), a Selective and Orally Efficacious Inhibitor of the Met Kinase Superfamily

Substituted N-(4-(2-aminopyridin-4-yloxy)-3-fluoro-phenyl)-1-(4-fluorophenyl)-2-oxo-1,2-dihydropyridine-3-carboxamides were identified as potent and selective Met kinase inhibitors. Substitution of the pyridine 3-position gave improved enzyme potency, while substitution of the pyridone 4-position led to improved aqueous solubility and kinase selectivity. Analogue 10 demonstrated complete tumor stasis in a Met-dependent GTL-16 human gastric carcinoma xenograft model following oral administration. Because of its excellent in vivo efficacy and favorable pharmacokinetic and preclinical safety profiles, 10 has been advanced into phase I clinical trials.

Discovery of Pyrrolopyridine-Pyridone Based Inhibitors of Met Kinase : Synthesis, X-ray Crystallographic Analysis, and Biological Activities

Conformationally constrained 2-pyridone analogue 2 is a potent Met kinase inhibitor with an IC50 value of 1.8 nM. Further SAR of the 2-pyridone based inhibitors of Met kinase led to potent 4-pyridone and pyridine N-oxide inhibitors such as 3 and 4. The X-ray crystallographic data of the inhibitor 2 bound to the ATP binding site of Met kinase protein provided insight into the binding modes of these inhibitors, and the SAR of this series of analogues was rationalized. Many of these analogues showed potent antiproliferative activities against the Met dependent GTL-16 gastric carcinoma cell line. Compound 2 also inhibited Flt-3 and VEGFR-2 kinases with IC50 values of 4 and 27 nM, respectively. It possesses a favorable pharmacokinetic profile in mice and demonstrates significant in vivo antitumor activity in the GTL-16 human gastric carcinoma xenograft model.

Discovery of orally active pyrrolopyridine- and aminopyridine-based Met kinase inhibitors

A series of acylurea analogs derived from pyrrolopyridine and aminopyridine scaffolds were identified as potent inhibitors of Met kinase activity. The SAR at various positions of the two kinase scaffolds was investigated. These studies led to the discovery of compounds 3b and 20b, which demonstrated favorable pharmacokinetic properties in mice and significant antitumor activity in a human gastric carcinoma xenograft model.

1H-Pyrazolo[3,4-b]pyridine inhibitors of cyclin-dependent kinases: highly potent 2,6-Difluorophenacyl analogues

Structure-activity studies of 1H-pyrazolo[3,4-b]pyridine 1 have resulted in the discovery of potent CDK1/CDK2 selective inhibitor 21h, BMS-265246 (CDK1/cycB IC(50)=6 nM, CDK2/cycE IC(50)=9 nM). The 2,6-difluorophenyl substitution was critical for potent inhibitory activity. A solid state structure of 21j, a close di-fluoro analogue, bound to CDK2 shows the inhibitor resides coincident with the ATP purine binding site and forms important H-bonds with Leu83 on the protein backbone.

1H-Pyrazolo[3,4-b]pyridine Inhibitors of Cyclin-Dependent Kinases

1H-Pyrazolo[3,4-b]pyridine 3 (SQ-67563) has been shown to be a potent, selective inhibitor of CDK1/CDK2 in vitro. In cells 3 acts as a cytotoxic agent with the ability to block cell cycle progression and/or induce apoptosis. The solid state structure of 3 bound to CDK2 shows 3 resides coincident with the ATP purine binding site and forms important H-bonding interactions with Leu83 on the protein backbone.

Prediction of Ligand−Receptor Binding Thermodynamics by Free Energy Force Field (FEFF) 3D-QSAR Analysis: Application to a Set of Peptidometic Renin Inhibitors

A methodology is presented and applied in which the accurate estimation of ligand-receptor binding thermodynamics is achieved by formulating the calculation as a QSAR problem. When the receptor geometry is known, the free energy force field (FEFF) ligand-receptor binding energy terms can be calculated and used as independent variables in constructing FEFF 3D-QSARs. The FEFF 3D-QSAR analysis of a series of transition state inhibitors of renin was carried out. From a statistical analysis of the free energy contributions to the binding process, FEFF 3D-QSARs were constructed that reveal the change in solvation free energy upon binding and the intramolecular vacuum internal energy of the ligand in the unbound state are the most significant FEFF terms in determining the binding free energy, delta G. Other terms, such as ligand stretching, bending, and torsion energy changes, the intermolecular van der Waals interaction energy, and change in ligand conformational entropy upon binding, are also found to make significant contributions in some FEFF 3D-QSAR delta G models and in delta H and delta S binding models. Overall, a relatively small number of the thermodynamic contributions to the ligand-receptor binding process dominates the thermodynamics of binding in a given model.

Identification of pyrrolo[2,1-f][1,2,4]triazine-based inhibitors of Met kinase.

An amide library derived from the pyrrolo[2,1-f][1,2,4]triazine scaffold led to the identification of modest inhibitors of Met kinase activity. Introduction of polar side chains at C-6 of the pyrrolotriazine core provided significant improvements in in vitro potency. The amide moiety could be replaced with acylurea and malonamide substituents to give compounds with improved potency in the Met-driven GTL-16 human gastric carcinoma cell line. Acylurea pyrrolotriazines with substitution at C-5 demonstrated single digit nanomolar kinase activity. X-ray crystallography revealed that the C-5 substituted pyrrolotriazines bind to the Met kinase domain in an ATP-competitive manner.

Analysis and optimization of structure-based virtual screening protocols. 2. Examination of docked ligand orientation sampling methodology: mapping a pharmacophore for success.

An important element of any structure-based virtual screening (SVS) technique is the method used to orient the ligands in the target active site. This has been a somewhat overlooked issue in recent SVS validation studies, with the assumption being made that the performance of an algorithm for a given set of orientation sampling settings will be representative for the general behavior of said technique. Here, we analyze five different SVS targets using a variety of sampling paradigms within the DOCK, GOLD and PROMETHEUS programs over a data set of approximately 10,000 noise compounds, combined with data sets containing multiple active compounds. These sets have been broken down by chemotype, with chemotype hit rate used to provide a measure of enrichment with a potentially improved relevance to real world SVS experiments. The variability in enrichment results produced by different sampling paradigms is illustrated, as is the utility of using pharmacophores to constrain sampling to regions that reflect known structural biology. The difference in results when comparing chemotype with compound hit rates is also highlighted.

CONSTRUCTING PROTEIN MODELS FOR LIGAND-RECEPTOR BINDING THERMODYNAMIC SIMULATIONS : AN APPLICATION TO A SET OF PEPTIDOMETIC RENIN INHIBITORS

Structure-based design is the application of ligand-receptor modeling to predict the activity of a series of molecules that bind to a common receptor for which the molecular geometry is available. Successful structure-based design requires an accurate receptor model which can be economically employed in the design calculations. One goal of the work reported here has been to reduce the size of a model structure of a macromolecular receptor to allow multiple ligand-receptor molecular dynamic (MD) simulations to be computationally economical yet still provide meaningful binding thermodynamic data. A scaled-down 10 A receptor model of the enzyme renin, when subjected to an alternate atomic mass constraint, maintains the structural integrity of the composite parent crystal structure. A second goal of the work has been to develop schemes to explore and characterize the protonation states of receptors and ligand-receptor systems. Application of the charge state characterization schemes to the hydroxyethylene and statine transition state inhibitors of renin in the training set suggests a monoprotonation state of the two active-site aspartate residues, where the lone proton resides on the outer carboxylate oxygen of Asp226 is most likely. For the reduced amide transition state inhibitors an active site consisting of both aspartates in the totally ionized state, and the ligand carrying a net +1.0 charge, is most stable and consistent with experimental data.