Modulating the interaction between CDK2 and cyclin A with a quinoline-based inhibitor.

A new class of quinoline-based kinase inhibitors has been discovered that both disrupt cyclin dependent 2 (CDK2) interaction with its cyclin A subunit and act as ATP competitive inhibitors. The key strategy for discovering this class of protein-protein disrupter compounds was to screen the monomer CDK2 in an affinity-selection/mass spectrometry-based technique and to perform secondary assays that identified compounds that bound only to the inactive CDK2 monomer and not the active CDK2/cyclin A heterodimer. Through a series of chemical modifications the affinity (Kd) of the original hit improved from 1 to 0.005µM.

Application of fragment-based NMR screening, X-ray crystallography, structure-based design, and focused chemical library design to identify novel microM leads for the development of nM BACE-1 (beta-site APP cleaving enzyme 1) inhibitors.

Fragment-based NMR screening, X-ray crystallography, structure-based design, and focused chemical library design were used to identify novel inhibitors for BACE-1. A rapid optimization of an initial NMR hit was achieved by a combination of NMR and a functional assay, resulting in the identification of an isothiourea hit with a K(d) of 15 microM for BACE-1. NMR data and the crystal structure revealed that this hit makes H-bond interactions with the two catalytic aspartates, occupies the nonprime side region of the active site of BACE-1, and extends toward the S3 subpocket (S3sp). A focused NMR-based search for heterocyclic isothiourea isosteres resulted in several distinct classes of BACE-1 active site directed compounds with improved chemical stability and physicochemical properties. The strategy for optimization of the 2-aminopyridine lead series to potent inhibitors of BACE-1 was demonstrated. The structure-based design of a cyclic acylguanidine lead series and its optimization into nanomolar BACE-1 inhibitors are the subject of the companion paper

Crystal structures of MEK1 binary and ternary complexes with nucleotides and inhibitors.

MEK1 is a member of the MAPK signal transduction pathway that responds to growth factors and cytokines. We have determined that the kinase domain spans residues 35-382 by proteolytic cleavage. The complete kinase domain has been crystallized and its X-ray crystal structure as a complex with magnesium and ATP-gammaS determined at 2.1 A. Unlike crystals of a truncated kinase domain previously published, the crystals of the intact domain can be grown either as a binary complex with a nucleotide or as a ternary complex with a nucleotide and one of a multitude of allosteric inhibitors. Further, the crystals allow for the determination of costructures with ATP competitive inhibitors. We describe the structures of nonphosphorylated MEK1 (npMEK1) binary complexes with ADP and K252a, an ATP-competitive inhibitor (see Table 1), at 1.9 and 2.7 A resolution, respectively. Ternary complexes have also been solved between npMEK1, a nucleotide, and an allosteric non-ATP competitive inhibitor: ATP-gammaS with compound 1 and ADP with either U0126 or the MEK1 clinical candidate PD325089 at 1.8, 2.0, and 2.5 A, respectively. Compound 1 is structurally similar to PD325901. These structures illustrate fundamental differences among various mechanisms of inhibition at the molecular level. Residues 44-51 have previously been shown to play a negative regulatory role in MEK1 activity. The crystal structure of the integral kinase domain provides a structural rationale for the role of these residues. They form helix A and repress enzymatic activity by stabilizing an inactive conformation in which helix C is displaced from its active state position. Finally, the structure provides for the first time a molecular rationale that explains how mutations in MEK may lead to the cardio-facio-cutaneous syndrome.

Induced-fit docking of mometasone furoate and further evidence for glucocorticoid receptor 17alpha pocket flexibility.

An induced-fit docking method was used to characterize the interactions of the glucocorticoid receptor binding-site with mometasone furoate, a glucocorticoid with a lipophilic ester at the C17alpha position. Two validation studies demonstrated that the protocol can reproduce crystal structures of nuclear receptors, and is appropriate for modeling ligand binding to the glucocorticoid receptor. Key hydrogen bonding interactions between mometasone furoate and the glucocorticoid receptor, as well as favorable hydrophobic interactions between the furoate group and the 17alpha pocket, contribute to high affinity and specificity of this ligand for the receptor. Using the glucocorticoid des-ciclesonide, which has an even larger moiety at the 16,17alpha position, induced-fit docking demonstrates the ability of the 17alpha pocket of the receptor to expand even further to accommodate the ligand.

Cross-docking of inhibitors into CDK2 structures. 1.

Predicting protein/ligand binding affinity is one of the most challenging computational chemistry tasks. Numerous methods have been developed to address this challenge, but they all have limitations. Failure to account for protein flexibility has been a shortcoming of many methods. In this cross-docking study the data set comprised 150 inhibitor complexes of the protein kinase CDK2. Gold and Glide performed well in terms of docking accuracy. The chance of cross-docking a ligand within a 2 A RMSD of its experimental pose was found to be 50%. Relative binding potency was not properly predicted from scoring functions, even though cross-docking of each inhibitor into each protein structure was performed and only scores of correctly docked ligands were considered. An accompanying paper (Voigt, J. H.; Elkin, C.; Madison, V. S. Duca, J. S. J. Chem. Inf. Model. 2008, 48, 669-678) covers cross-docking and docking accuracy from the perspective of using multiple protein structures.

Cross-docking of inhibitors into CDK2 structures. 2.

In the preceding paper (Duca, J. S.; Madison, V. S.; Voigt, J. H. J. Chem. Inf. Model. 2008, 48, 659-668), the accuracy of docking and affinity predictions of the Gold and Glide programs were investigated using single protein conformations spanning 150 CDK2/inhibitor crystallographic complexes. High docking accuracy was observed with both methods; furthermore, Glide showed modest log(IC50)/score correlations. In this part of the study, the effect of combining docking results from multiple protein conformations in a consensus fashion was probed. This approach enhanced docking accuracy only for Glide, which was attributed to the nature of its scoring function. For log(IC50)/score correlations, particular emphasis was placed on considering only scores from correctly docked poses. Using multiple instead of single protein structures showed an improvement in the correlations. Validation sets and scrambling experiments were used to examine the statistical significance and predictivity of these correlations. Rather than actual improvements in scoring accuracy, docking to multiple protein conformations produced overfitting artifacts.

Identification of novel cannabinoid CB1 receptor antagonists by using virtual screening with a pharmacophore model.

CB1 receptor antagonists have proven to be clinically effective in treating obesity and related disorders. We report here the identification of a novel class of azetidinone CB1 antagonists by using virtual screening methods. For this purpose, we developed a pharmacophore model based on known representative CB1 antagonists and employed it to screen a database of about a half million Schering-Plough compounds. We applied a stepwise filtering protocol based on molecular weight, compound availability, and a modified rule-of-five to reduce the number of hits. We then combined Bayesian modeling and clustering techniques to select a final set of 420 compounds for in vitro testing. Five compounds were found to have >50% inhibition at 100 nM in a CB1 competitive binding assay and were further characterized by using both CB1 and CB2 assays. The most potent compound has a CB1 K i of 53 nM and >5-fold selectivity against the CB2 receptor.

Design and synthesis of depeptidized macrocyclic inhibitors of hepatitis C NS3-4A protease using structure-based drug design.

Hepatitis C virus (HCV) NS3, when bound to NS-4A cofactor, facilitates development of mature virons by catalyzing cleavage of a polyprotein to form functional and structural proteins of HCV. The enzyme has a shallow binding pocket at the catalytic site, making development of inhibitors difficult. We have designed, preorganized, and depeptidized macrocyclic inhibitors from P(4) to P(2)' and optimized binding to 0.1 microM. The structure of an inhibitor bound to the enzyme was also solved.

Insights from ab initio quantum chemical calculations into the preferred tautomeric forms and binding affinities to CDK2 of substituted pyrazolopyridines.

Ab initio calculations were employed to compute pKa values and tautomer properties of a series of substituted pyrazolopyridines. The results show that the neutral 1H tautomer predominates, but upon protonation this proton migrates to give the preferred charged [2H,7H] tautomer. The basicity of the pyrazolopyridines is correlated with the electron donating capability of the 4-substituent. Ab initio free energy calculations were also used to identify determinants of binding affinity for some recently published pyrazolopyridine inhibitors of CDK2. Hydrogen-bonding affinity may be one important component of binding strength.

NMR characterization of interleukin-2 in complexes with the IL-2Ralpha receptor component, and with low molecular weight compounds that inhibit the IL-2/IL-Ralpha interaction.

Nuclear magnetic resonance (NMR) methods were employed to study the interaction of the cytokine Interleukin-2 (IL-2) with the alpha-subunit of its receptor (IL-2Ralpha), and to help understand the behavior of small molecule inhibitors of this interaction. Heteronuclear (1)H-(15)N HSQC experiments were used to identify the interaction surface of (15)N-enriched Interleukin-2 ((15)N-IL-2) in complex with human IL-2Ralpha. In these experiments, chemical shift and line width changes in the heteronuclear single-quantum coherence (HSQC) spectra upon binding of (15)N-IL-2 enabled classification of NH atoms as either near to, or far from, the IL-2Ralpha interaction surface. These data were complemented by hydrogen/deuterium (H/D) exchange measurements, which illustrated enhanced protection of slowly-exchanging IL-2 NH protons near the site of interaction with IL-2Ralpha. The interaction surface defined by NMR compared well with the IL-2Ralpha binding site identified previously using mutagenesis of human and murine IL-2. Two low molecular weight inhibitors of the IL-2/IL-2Ralpha interaction were studied: one (a cyclic peptide derivative) was found to mimic a part of the cytokine and bind to IL-2Ralpha; the other (an acylphenylalanine derivative) was found to bind to IL-2. For the interaction between IL-2 and the acylphenylalanine, chemical shift perturbations of (15)N and (15)NH backbone resonances were tracked as a function of ligand concentration. The perturbation pattern observed for this complex revealed that the acylphenylalanine is a competitive inhibitor-it binds to the same site on IL-2 that interacts with IL-2Ralpha.

The active site of a zinc-dependent metalloproteinase influences the computed pK(a) of ligands coordinated to the catalytic zinc ion.

TNF-alpha converting enzyme (TACE) is a multidomain, membrane-anchored protein that includes a Zn-dependent protease domain. It releases the soluble form of cytokine tumor necrosis factor-alpha (TNF-alpha) from its membrane-bound precursor. TACE is a metalloprotease containing a catalytic glutamic acid, Glu-406, and a Zn(2+) ion ligated to three imidazoles. The protonation states of the active site glutamic acid and inhibitors are important factors in understanding the potency of inhibitors with acidic zinc-ligating groups such as hydroxamic and carboxylic acids. Density functional methods were utilized to compute pK(a) values using a model of the catalytic site of TACE and to predict a concomitant mechanism of binding, consistent with lowering the pK(a) of the bound ligand and raising the pK(a) of the active site Glu-406. Weak acids, such as hydroxamic acids, bind in their neutral form and then transfer an acidic proton to Glu-406. Stronger acids, such as carboxylic acids, bind in their anionic form and require preprotonation of Glu-406. Similar binding events would be expected for other zinc-dependent proteases.