Predicting the Conformational Variability of Abl Tyrosine Kinase using Molecular Dynamics Simulations and Markov State Models.

Understanding protein conformational variability remains a challenge in drug discovery. The issue arises in protein kinases, whose multiple conformational states can affect the binding of small-molecule inhibitors. To overcome this challenge, we propose a comprehensive computational framework based on Markov state models (MSMs). Our framework integrates the information from explicit-solvent molecular dynamics simulations to accurately rank-order the accessible conformational variants of a target protein. We tested the methodology using Abl kinase with a reference and blind-test set. Only half of the Abl conformational variants discovered by our approach are present in the disclosed X-ray structures. The approach successfully identified a protein conformational state not previously observed in public structures but evident in a retrospective analysis of Lilly in-house structures: the X-ray structure of Abl with WHI-P154. Using a MSM-derived model, the free energy landscape and kinetic profile of Abl was analyzed in detail highlighting opportunities for targeting the unique metastable states.

Utilization of peptide phage display to investigate hotspots on IL-17A and what it means for drug discovery.

To date, IL-17A antibodies remain the only therapeutic approach to correct the abnormal activation of the IL-17A/IL-17R signaling complex. Why is it that despite the remarkable success of IL-17 antibodies, there is no small molecule antagonist of IL-17A in the clinic? Here we offer a unique approach to address this question. In order to understand the interaction of IL-17A with its receptor, we combined peptide discovery using phage display with HDX, crystallography, and functional assays to map and characterize hot regions that contribute to most of the energetics of the IL-17A/IL-17R interaction. These functional maps are proposed to serve as a guide to aid in the development of small molecules that bind to IL-17A and block its interaction with IL-17RA.

Synthesis and Pharmacological Characterization of C4ß-Amide-Substituted 2-Aminobicyclo[3.1.0]hexane-2,6-dicarboxylates. Identification of (1 S,2 S,4 S,5 R,6 S)-2-Amino-4-[(3-methoxybenzoyl)amino]bicyclo[3.1.0]hexane-2,6-dicarboxylic Acid (LY2794193), a Highly Potent and Selective mGlu3 Receptor Agonist.

Multiple therapeutic opportunities have been suggested for compounds capable of selective activation of metabotropic glutamate 3 (mGlu3) receptors, but small molecule tools are lacking. As part of our ongoing efforts to identify potent, selective, and systemically bioavailable agonists for mGlu2 and mGlu3 receptor subtypes, a series of C4ß-N-linked variants of (1 S,2 S,5 R,6 S)-2-amino-bicyclo[3.1.0]hexane-2,6-dicarboxylic acid 1 (LY354740) were prepared and evaluated for both mGlu2 and mGlu3 receptor binding affinity and functional cellular responses. From this investigation we identified (1 S,2 S,4 S,5 R,6 S)-2-amino-4-[(3-methoxybenzoyl)amino]bicyclo[3.1.0]hexane-2,6-dicarboxylic acid 8p (LY2794193), a molecule that demonstrates remarkable mGlu3 receptor selectivity. Crystallization of 8p with the amino terminal domain of hmGlu3 revealed critical binding interactions for this ligand with residues adjacent to the glutamate binding site, while pharmacokinetic assessment of 8p combined with its effect in an mGlu2 receptor-dependent behavioral model provides estimates for doses of this compound that would be expected to selectively engage and activate central mGlu3 receptors in vivo.

Synthesis and Pharmacological Characterization of C4-(Thiotriazolyl)-substituted-2-aminobicyclo[3.1.0]hexane-2,6-dicarboxylates. Identification of (1R,2S,4R,5R,6R)-2-Amino-4-(1H-1,2,4-triazol-3-ylsulfanyl)bicyclo[3.1.0]hexane-2,6-dicarboxylic Acid (LY2812223), a Highly Potent, Functionally Selective mGlu2 Receptor Agonist.

Identification of orthosteric mGlu(2/3) receptor agonists capable of discriminating between individual mGlu2 and mGlu3 subtypes has been highly challenging owing to the glutamate-site sequence homology between these proteins. Herein we detail the preparation and characterization of a series of molecules related to (1S,2S,5R,6S)-2-aminobicyclo[3.1.0]hexane-2,6-dicarboxylate 1 (LY354740) bearing C4-thiotriazole substituents. On the basis of second messenger responses in cells expressing other recombinant human mGlu2/3 subtypes, a number of high potency and efficacy mGlu2 receptor agonists exhibiting low potency mGlu3 partial agonist/antagonist activity were identified. From this, (1R,2S,4R,5R,6R)-2-amino-4-(1H-1,2,4-triazol-3-ylsulfanyl)bicyclo[3.1.0]hexane-2,6-dicarboxylic acid 14a (LY2812223) was further characterized. Cocrystallization of 14a with the amino terminal domains of hmGlu2 and hmGlu3 combined with site-directed mutation studies has clarified the underlying molecular basis of this unique pharmacology. Evaluation of 14a in a rat model responsive to mGlu2 receptor activation coupled with a measure of central drug disposition provides evidence that this molecule engages and activates central mGlu2 receptors in vivo.

Crystal Structures of mPGES-1 Inhibitor Complexes Form a Basis for the Rational Design of Potent Analgesic and Anti-Inflammatory Therapeutics.

Microsomal prostaglandin E synthase 1 (mPGES-1) is an α-helical homotrimeric integral membrane inducible enzyme that catalyzes the formation of prostaglandin E2 (PGE2) from prostaglandin H2 (PGH2). Inhibition of mPGES-1 has been proposed as a therapeutic strategy for the treatment of pain, inflammation, and some cancers. Interest in mPGES-1 inhibition can, in part, be attributed to the potential circumvention of cardiovascular risks associated with anti-inflammatory cyclooxygenase 2 inhibitors (coxibs) by targeting the prostaglandin pathway downstream of PGH2 synthesis and avoiding suppression of antithrombotic prostacyclin production. We determined the crystal structure of mPGES-1 bound to four potent inhibitors in order to understand their structure-activity relationships and provide a framework for the rational design of improved molecules. In addition, we developed a light-scattering-based thermal stability assay to identify molecules for crystallographic studies.

Synthesis and pharmacological characterization of C4-disubstituted analogs of 1S,2S,5R,6S-2-aminobicyclo[3.1.0]hexane-2,6-dicarboxylate: identification of a potent, selective metabotropic glutamate receptor agonist and determination of agonist-bound human mGlu2 and mGlu3 amino terminal domain structures.

As part of our ongoing research to identify novel agents acting at metabotropic glutamate 2 (mGlu2) and 3 (mGlu3) receptors, we have previously reported the identification of the C4α-methyl analog of mGlu2/3 receptor agonist 1 (LY354740). This molecule, 1S,2S,4R,5R,6S-2-amino-4-methylbicyclo[3.1.0]hexane-2,6-dicarboxylate 2 (LY541850), exhibited an unexpected mGlu2 agonist/mGlu3 antagonist pharmacological profile, whereas the C4ß-methyl diastereomer (3) possessed dual mGlu2/3 receptor agonist activity. We have now further explored this structure-activity relationship through the preparation of cyclic and acyclic C4-disubstituted analogs of 1, leading to the identification of C4-spirocyclopropane 5 (LY2934747), a novel, potent, and systemically bioavailable mGlu2/3 receptor agonist which exhibits both antipsychotic and analgesic properties in vivo. In addition, through the combined use of protein-ligand X-ray crystallography employing recombinant human mGlu2/3 receptor amino terminal domains, molecular modeling, and site-directed mutagenesis, a molecular basis for the observed pharmacological profile of compound 2 is proposed.

Structural context of disease-associated mutations and putative mechanism of autoinhibition revealed by X-ray crystallographic analysis of the EZH2-SET domain.

The enhancer-of-zeste homolog 2 (EZH2) gene product is an 87 kDa polycomb group (PcG) protein containing a C-terminal methyltransferase SET domain. EZH2, along with binding partners, i.e., EED and SUZ12, upon which it is dependent for activity forms the core of the polycomb repressive complex 2 (PRC2). PRC2 regulates gene silencing by catalyzing the methylation of histone H3 at lysine 27. Both overexpression and mutation of EZH2 are associated with the incidence and aggressiveness of various cancers. The novel crystal structure of the SET domain was determined in order to understand disease-associated EZH2 mutations and derive an explanation for its inactivity independent of complex formation. The 2.00 Å crystal structure reveals that, in its uncomplexed form, the EZH2 C-terminus folds back into the active site blocking engagement with substrate. Furthermore, the S-adenosyl-L-methionine (SAM) binding pocket observed in the crystal structure of homologous SET domains is notably absent. This suggests that a conformational change in the EZH2 SET domain, dependent upon complex formation, must take place for cofactor and substrate binding activities to be recapitulated. In addition, the data provide a structural context for clinically significant mutations found in the EZH2 SET domain.

The 2.5 Å crystal structure of the SIRT1 catalytic domain bound to nicotinamide adenine dinucleotide (NAD+) and an indole (EX527 analogue) reveals a novel mechanism of histone deacetylase inhibition.

The sirtuin SIRT1 is a NAD(+)-dependent histone deacetylase, a Sir2 family member, and one of seven human sirtuins. Sirtuins are conserved from archaea to mammals and regulate transcription, genome stability, longevity, and metabolism. SIRT1 regulates transcription via deacetylation of transcription factors such as PPARγ, NFκB, and the tumor suppressor protein p53. EX527 (27) is a nanomolar SIRT1 inhibitor and a micromolar SIRT2 inhibitor. To elucidate the mechanism of SIRT inhibition by 27, we determined the 2.5 Å crystal structure of the SIRT1 catalytic domain (residues 241-516) bound to NAD(+) and the 27 analogue compound 35. 35 binds deep in the catalytic cleft, displacing the NAD(+) nicotinamide and forcing the cofactor into an extended conformation. The extended NAD(+) conformation sterically prevents substrate binding. The SIRT1/NAD(+)/35 crystal structure defines a novel mechanism of histone deacetylase inhibition and provides a basis for understanding, and rationally improving, inhibition of this therapeutically important target by drug-like molecules.

Crystal structure of the human PRMT5:MEP50 complex.

Protein arginine methyltransferases (PRMTs) play important roles in several cellular processes, including signaling, gene regulation, and transport of proteins and nucleic acids, to impact growth, differentiation, proliferation, and development. PRMT5 symmetrically di-methylates the two-terminal ω-guanidino nitrogens of arginine residues on substrate proteins. PRMT5 acts as part of a multimeric complex in concert with a variety of partner proteins that regulate its function and specificity. A core component of these complexes is the WD40 protein MEP50/WDR77/p44, which mediates interactions with binding partners and substrates. We have determined the crystal structure of human PRMT5 in complex with MEP50 (methylosome protein 50), bound to an S-adenosylmethionine analog and a peptide substrate derived from histone H4. The structure of the surprising hetero-octameric complex reveals the close interaction between the seven-bladed ß-propeller MEP50 and the N-terminal domain of PRMT5, and delineates the structural elements of substrate recognition.

Structural bases of PAS domain-regulated kinase (PASK) activation in the absence of activation loop phosphorylation.

Per-Arnt-Sim (PAS) domain-containing protein kinase (PASK) is an evolutionary conserved protein kinase that coordinates cellular metabolism with metabolic demand in yeast and mammals. The molecular mechanisms underlying PASK regulation, however, remain unknown. Herein, we describe a crystal structure of the kinase domain of human PASK, which provides insights into the regulatory mechanisms governing catalysis. We show that the kinase domain adopts an active conformation and has catalytic activity in vivo and in vitro in the absence of activation loop phosphorylation. Using site-directed mutagenesis and structural comparison with active and inactive kinases, we identified several key structural features in PASK that enable activation loop phosphorylation-independent activity. Finally, we used combinatorial peptide library screening to determine that PASK prefers basic residues at the P-3 and P-5 positions in substrate peptides. Our results describe the key features of the PASK structure and how those features are important for PASK activity and substrate selection.

SGX523 is an exquisitely selective, ATP-competitive inhibitor of the MET receptor tyrosine kinase with antitumor activity in vivo.

The MET receptor tyrosine kinase has emerged as an important target for the development of novel cancer therapeutics. Activation of MET by mutation or gene amplification has been linked to kidney, gastric, and lung cancers. In other cancers, such as glioblastoma, autocrine activation of MET has been demonstrated. Several classes of ATP-competitive inhibitor have been described, which inhibit MET but also other kinases. Here, we describe SGX523, a novel, ATP-competitive kinase inhibitor remarkable for its exquisite selectivity for MET. SGX523 potently inhibited MET with an IC50 of 4 nmol/L and is >1,000-fold selective versus the >200-fold selectivity of other protein kinases tested in biochemical assays. Crystallographic study revealed that SGX523 stabilizes MET in a unique inactive conformation that is inaccessible to other protein kinases, suggesting an explanation for the selectivity. SGX523 inhibited MET-mediated signaling, cell proliferation, and cell migration at nanomolar concentrations but had no effect on signaling dependent on other protein kinases, including the closely related RON, even at micromolar concentrations. SGX523 inhibition of MET in vivo was associated with the dose-dependent inhibition of growth of tumor xenografts derived from human glioblastoma and lung and gastric cancers, confirming the dependence of these tumors on MET catalytic activity. Our results show that SGX523 is the most selective inhibitor of MET catalytic activity described to date and is thus a useful tool to investigate the role of MET kinase in cancer without the confounding effects of promiscuous protein kinase inhibition.

Structure of a microsporidian methionine aminopeptidase type 2 complexed with fumagillin and TNP-470.

Microsporidia are protists that have been reported to cause infections in both vertebrates and invertebrates. They have emerged as human pathogens particularly in patients that are immunosuppressed and cases of gastrointestinal infection, encephalitis, keratitis, sinusitis, myositis and disseminated infection are well described in the literature. While benzimidazoles are active against many species of microsporidia, these drugs do not have significant activity against Enterocytozoon bieneusi. Fumagillin and its analogues have been demonstrated to have activity invitro and in animal models of microsporidiosis and human infections due to E. bieneusi. Fumagillin and its analogues inhibit methionine aminopeptidase type 2. Encephalitozoon cuniculi MetAP2 (EcMetAP2) was cloned and expressed as an active enzyme using a baculovirus system. The crystal structure of EcMetAP2 was determined with and without the bound inhibitors fumagillin and TNP-470. This structure classifies EcMetAP2 as a member of the MetAP2c family. The EcMetAP2 structure was used to generate a homology model of the E. bieneusi MetAP2. Comparison of microsporidian MetAP2 structures with human MetAP2 provides insights into the design of inhibitors that might exhibit specificity for microsporidian MetAP2.

A novel mode of Gleevec binding is revealed by the structure of spleen tyrosine kinase.

Spleen tyrosine kinase (Syk) is a non-receptor tyrosine kinase required for signaling from immunoreceptors in various hematopoietic cells. Phosphorylation of two tyrosine residues in the activation loop of the Syk kinase catalytic domain is necessary for signaling, a phenomenon typical of tyrosine kinase family members. Syk in vitro enzyme activity, however, does not depend on phosphorylation (activation loop tyrosine --> phenylalanine mutants retain catalytic activity). We have determined the x-ray structure of the unphosphorylated form of the kinase catalytic domain of Syk. The enzyme adopts a conformation of the activation loop typically seen only in activated, phosphorylated tyrosine kinases, explaining why Syk does not require phosphorylation for activation. We also demonstrate that Gleevec (STI-571, Imatinib) inhibits the isolated kinase domains of both unphosphorylated Syk and phosphorylated Abl with comparable potency. Gleevec binds Syk in a novel, compact cis-conformation that differs dramatically from the binding mode observed with unphosphorylated Abl, the more Gleevec-sensitive form of Abl. This finding suggests the existence of two distinct Gleevec binding modes: an extended, trans-conformation characteristic of tight binding to the inactive conformation of a protein kinase and a second compact, cis-conformation characteristic of weaker binding to the active conformation. Finally, the Syk-bound cis-conformation of Gleevec bears a striking resemblance to the rigid structure of the nonspecific, natural product kinase inhibitor staurosporine.