Crystal structure of the tyrosine kinase domain of colony-stimulating factor-1 receptor (cFMS) in complex with two inhibitors

The cFMS proto-oncogene encodes for the colony-stimulating factor-1 receptor, a receptor-tyrosine kinase responsible for the differentiation and maturation of certain macrophages. Upon binding its ligand colony-stimulating factor-1 cFMS autophosphorylates, dimerizes, and induces phosphorylation of downstream targets. We report the novel crystal structure of unphosphorylated cFMS in complex with two members of different classes of drug-like protein kinase inhibitors. cFMS exhibits a typical bi-lobal kinase fold, and its activation loop and DFG motif are found to be in the canonical inactive conformation. Both ATP competitive inhibitors are bound in the active site and demonstrate a binding mode similar to that of STI-571 bound to cABL. The DFG motif is prevented from switching into the catalytically competent conformation through interactions with the inhibitors. Activation of cFMS is also inhibited by the juxtamembrane domain, which interacts with residues of the active site and prevents formation of the activated kinase. Together the structures of cFMS provide further insight into the autoinhibition of receptor-tyrosine kinases via their respective juxtamembrane domains; additionally the binding mode of two novel classes of kinase inhibitors will guide the design of novel molecules targeting macrophage-related diseases.     

Crystal structures of the N-terminal kinase domain of human RSK1 bound to three different ligands: Implications for the design of RSK1 specific inhibitors

The p90 ribosomal S6 kinases (RSKs) also known as MAPKAP-Ks are serine/threonine protein kinases that are activated by ERK or PDK1 and act as downstream effectors of mitogen-activated protein kinase (MAPK). RSK1, a member of the RSK family, contains two distinct kinase domains in a single polypeptide chain, the regulatory C-terminal kinase domain (CTKD) and the catalytic N-terminal kinase domain (NTKD). Autophosphorylation of the CTKD leads to activation of the NTKD that subsequently phosphorylates downstream substrates. Here we report the crystal structures of the unactivated RSK1 NTKD bound to different ligands at 2.0 A resolution. The activation loop and helix alphaC, key regulatory elements of kinase function, are disordered. The DFG motif of the inactive RSK1 adopts an "active-like" conformation. The beta-PO(4) group in the AMP-PCP complex adopts a unique conformation that may contribute to inactivity of the enzyme. Structures of RSK1 ligand complexes offer insights into the design of novel anticancer agents and into the regulation of the catalytic activity of RSKs.     

Cutting Edge: IL-1 receptor-associated kinase 4 structures reveal novel features and multiple conformations

IL-1R-associated kinase (IRAK)4 plays a central role in innate and adaptive immunity, and is a crucial component in IL-1/TLR signaling. We have determined the crystal structures of the apo and ligand-bound forms of human IRAK4 kinase domain. These structures reveal several features that provide opportunities for the design of selective IRAK4 inhibitors. The N-terminal lobe of the IRAK4 kinase domain is structurally distinctive due to a loop insertion after an extended N-terminal helix. The gatekeeper residue is a tyrosine, a unique feature of the IRAK family. The IRAK4 structures also provide insights into the regulation of its activity. In the apo structure, two conformations coexist, differing in the relative orientation of the two kinase lobes and the position of helix C. In the presence of an ATP analog only one conformation is observed, indicating that this is the active conformation.     

Structures of human Bruton's tyrosine kinase in active and inactive conformations suggest a mechanism of activation for TEC family kinases

Bruton's tyrosine kinase (BTK), a member of the TEC family of kinases, plays a crucial role in B‐cell maturation and mast cell activation. Although the structures of the unphosphorylated mouse BTK kinase domain and the unphosphorylated and phosphorylated kinase domains of human ITK are known, understanding the kinase selectivity profiles of BTK inhibitors has been hampered by the lack of availability of a high resolution, ligand‐bound BTK structure. Here, we report the crystal structures of the human BTK kinase domain bound to either Dasatinib (BMS‐354825) at 1.9 Å resolution or to 4‐amino‐5‐(4‐phenoxyphenyl)‐7H‐pyrrolospyrimidin‐ 7‐yl‐cyclopentane at 1.6 Å resolution. This data provides information relevant to the development of small molecule inhibitors targeting BTK and the TEC family of nonreceptor tyrosine kinases. Analysis of the structural differences between the TEC and Src families of kinases near the Trp‐Glu‐Ile motif in the N‐terminal region of the kinase domain suggests a mechanism of regulation of the TEC family members.     

Structural insights into the autoactivation mechanism of p21-activated protein kinase

p21-activated kinases (PAKs) play an important role in diverse cellular processes. Full activation of PAKs requires autophosphorylation of a critical threonine/serine located in the activation loop of the kinase domain. Here we report crystal structures of the phosphorylated and unphosphorylated PAK1 kinase domain. The phosphorylated PAK1 kinase domain has a conformation typical of all active protein kinases. Interestingly, the structure of the unphosphorylated PAK1 kinase domain reveals an unusual dimeric arrangement expected in an authentic enzyme-substrate complex, in which the activation loop of the putative "substrate" is projected into the active site of the "enzyme." The enzyme is bound to AMP-PNP and has an active conformation, whereas the substrate is empty and adopts an inactive conformation. Thus, the structure of the asymmetric homodimer mimics a trans-autophosphorylation complex, and suggests that unphosphorylated PAK1 could dynamically adopt both the active and inactive conformations in solution.     

A Src-like inactive conformation in the abl tyrosine kinase domain

The improper activation of the Abl tyrosine kinase results in chronic myeloid leukemia (CML). The recognition of an inactive conformation of Abl, in which a catalytically important Asp-Phe-Gly (DFG) motif is flipped by approximately 180 degrees with respect to the active conformation, underlies the specificity of the cancer drug imatinib, which is used to treat CML. The DFG motif is not flipped in crystal structures of inactive forms of the closely related Src kinases, and imatinib does not inhibit c-Src. We present a structure of the kinase domain of Abl, determined in complex with an ATP-peptide conjugate, in which the protein adopts an inactive conformation that resembles closely that of the Src kinases. An interesting aspect of the Src-like inactive structure, suggested by molecular dynamics simulations and additional crystal structures, is the presence of features that might facilitate the flip of the DFG motif by providing room for the phenylalanine to move and by coordinating the aspartate side chain as it leaves the active site. One class of mutations in BCR-Abl that confers resistance to imatinib appears more likely to destabilize the inactive Src-like conformation than the active or imatinib-bound conformations. Our results suggest that interconversion between distinctly different inactive conformations is a characteristic feature of the Abl kinase domain.     

Crystal structures of the Mnk2 kinase domain reveal an inhibitory conformation and a zinc binding site

Human mitogen-activated protein kinases (MAPK)-interacting kinases 1 and 2 (Mnk1 and Mnk2) target the translational machinery by phosphorylation of the eukaryotic initiation factor 4E (eIF4E). Here, we present the 2.1 A crystal structure of a nonphosphorylated Mnk2 fragment that encompasses the kinase domain. The results show Mnk-specific features such as a zinc binding motif and an atypical open conformation of the activation segment. In addition, the ATP binding pocket contains an Asp-Phe-Asp (DFD) in place of the canonical magnesium binding Asp-Phe-Gly (DFG) motif. The phenylalanine of this motif sticks into the ATP binding pocket and blocks ATP binding as observed with inhibitor bound and, thus, inactive p38 kinase. Replacement of the DFD by the canonical DFG motif affects the conformation of Mnk2, but not ATP binding and kinase activity. The results suggest that the ATP binding pocket and the activation segment of Mnk2 require conformational switches to provide kinase activity.     

A crystal structure of the human protein kinase Mps1 reveals an ordered conformation of the activation loop

Monopolar spindle 1 (Mps1) is a dual-specificity protein kinase, orchestrating faithful chromosome segregation during mitosis. All reported structures of the Mps1 kinase adopt the hallmarks of an inactive conformation, which includes a mostly disordered activation loop. Here, we present a 2.4 Å resolution crystal structure of an “extended” version of the Mps1 kinase domain, which shows an ordered activation loop. However, the other structural characteristics of an active kinase are not present. Our structure shows that the Mps1 activation loop can fit to the ATP binding pocket and interferes with ATP, but less so with inhibitors binding, partly explain the potency of various Mps1 inhibitors.     

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.     

Insights from free-energy calculations: protein conformational equilibrium, driving forces, and ligand-binding modes

Accurate free-energy calculations provide mechanistic insights into molecular recognition and conformational equilibrium. In this work, we performed free-energy calculations to study the thermodynamic properties of different states of molecular systems in their equilibrium basin, and obtained accurate absolute binding free-energy calculations for protein-ligand binding using a newly developed M2 algorithm. We used a range of Asp-Phe-Gly (DFG)-in/out p38α mitogen-activated protein kinase inhibitors as our test cases. We also focused on the flexible DFG motif, which is closely connected to kinase activation and inhibitor binding. Our calculations explain the coexistence of DFG-in and DFG-out states of the loop and reveal different components (e.g., configurational entropy and enthalpy) that stabilize the apo p38α conformations. To study novel ligand-binding modes and the key driving forces behind them, we computed the absolute binding free energies of 30 p38α inhibitors, including analogs with unavailable experimental structures. The calculations revealed multiple stable, complex conformations and changes in p38α and inhibitor conformations, as well as balance in several energetic terms and configurational entropy loss. The results provide relevant physics that can aid in designing inhibitors and understanding protein conformational equilibrium. Our approach is fast for use with proteins that contain flexible regions for structure-based drug design.     

Structure of the dimeric autoinhibited conformation of DAPK2, a pro-apoptotic protein kinase

The death-associated protein kinase (DAPK) family has been characterized as a group of pro-apoptotic serine/threonine kinases that share specific structural features in their catalytic kinase domain. Two of the DAPK family members, DAPK1 and DAPK2, are calmodulin-dependent protein kinases that are regulated by oligomerization, calmodulin binding, and autophosphorylation. In this study, we have determined the crystal and solution structures of murine DAPK2 in the presence of the autoinhibitory domain, with and without bound nucleotides in the active site. The crystal structure shows dimers of DAPK2 in a conformation that is not permissible for protein substrate binding. Two different conformations were seen in the active site upon the introduction of nucleotide ligands. The monomeric and dimeric forms of DAPK2 were further analyzed for solution structure, and the results indicate that the dimers of DAPK2 are indeed formed through the association of two apposed catalytic domains, as seen in the crystal structure. The structures can be further used to build a model for DAPK2 autophosphorylation and to compare with closely related kinases, of which especially DAPK1 is an actively studied drug target. Our structures also provide a model for both homodimerization and heterodimerization of the catalytic domain between members of the DAPK family. The fingerprint of the DAPK family, the basic loop, plays a central role in the dimerization of the kinase domain.     

Mutagenesis of p38alpha MAP kinase establishes key roles of Phe169 in function and structural dynamics and reveals a novel DFG-OUT state

In order to study the role of Phe169 in p38alpha MAP kinase structure and function, wild-type p38alpha and five p38alpha DFG motif mutants were examined in vitro for phosphorylation by MKK6, kinase activity toward ATF2 substrate, thermal stability, and X-ray crystal structure. All six p38alpha variants were efficiently phosphorylated by MKK6. However, only one activated p38alpha mutant (F169Y) possessed measurable kinase activity (1% compared to wild-type). The loss of kinase activity among the DFG mutants may result from an inability to correctly position Asp168 in the activated form of p38alpha. Two mutations significantly increased the thermal stability of p38alpha (F169A DeltaTm = 1.3 degrees C and D168G DeltaTm = 3.8 degrees C), and two mutations significantly decreased the stability of p38alpha (F169R DeltaTm = -3.2 degrees C and F169G DeltaTm = -4.7 degrees C). Interestingly, X-ray crystal structures of two thermally destabilized p38alpha-F169R and p38alpha-F169G mutants revealed a DFG-OUT conformation in the absence of an inhibitor molecule. This DFG-OUT conformation, termed alpha-DFG-OUT, is different from the ones previously identified in p38alpha crystal structures with bound inhibitors and postulated from high-temperature molecular dynamics simulations. Taken together, these results indicate that Phe169 is optimized for p38alpha functional activity and structural dynamics, rather than for structural stability. The alpha-DFG-OUT conformation observed for p38alpha-F169R and p38alpha-F169G may represent a naturally occurring intermediate state of p38alpha that provides access for binding of allosteric inhibitors. A model of the local forces driving the DFG IN-OUT transition in p38alpha is proposed.     

Crystal Structures of the FAK Kinase in Complex with TAE226 and Related Bis-Anilino Pyrimidine Inhibitors Reveal a Helical DFG Conformation

Focal Adhesion Kinase (FAK) is a non-receptor tyrosine kinase required for cell migration, proliferation and survival. FAK overexpression has been documented in diverse human cancers and is associated with a poor clinical outcome. Recently, a novel bis-anilino pyrimidine inhibitor, TAE226, was reported to efficiently inhibit FAK signaling, arrest tumor growth and invasion and prolong the life of mice with glioma or ovarian tumor implants. Here we describe the crystal structures of the FAK kinase bound to TAE226 and three related bis-anilino pyrimidine compounds. TAE226 induces a conformation of the N-terminal portion of the kinase activation loop that is only observed in FAK, but is distinct from the conformation in both the active and inactive states of the kinase. This conformation appears to require a glycine immediately N-terminal to the “DFG motif”, which adopts a helical conformation stabilized by interactions with TAE226. The presence of a glycine residue in this position contributes to the specificity of TAE226 and related compounds for FAK. Our work highlights the fact that kinases can access conformational space that is not necessarily utilized for their native catalytic regulation, and that such conformations can explain and be exploited for inhibitor specificity.     

Allosteric activation mechanism of the Mycobacterium tuberculosis receptor Ser/Thr protein kinase, PknB

The essential Mycobacterium tuberculosis Ser/Thr protein kinase (STPK), PknB, plays a key role in regulating growth and division, but the structural basis of activation has not been defined. Here, we provide biochemical and structural evidence that dimerization through the kinase-domain (KD) N-lobe activates PknB by an allosteric mechanism. Promoting KD pairing using a small-molecule dimerizer stimulates the unphosphorylated kinase, and substitutions that disrupt N-lobe pairing decrease phosphorylation activity in vitro and in vivo. Multiple crystal structures of two monomeric PknB KD mutants in complex with nucleotide reveal diverse inactive conformations that contain large active-site distortions that propagate > 30 ? from the mutation site. These results define flexible, inactive structures of a monomeric bacterial receptor KD and show how "back-to-back" N-lobe dimerization stabilizes the active KD conformation. This general mechanism of bacterial receptor STPK activation affords insights into the regulation of homologous eukaryotic kinases that form structurally similar dimers.     

Structural propensities of kinase family proteins from a Potts model of residue co‐variation

Understanding the conformational propensities of proteins is key to solving many problems in structural biology and biophysics. The co‐variation of pairs of mutations contained in multiple sequence alignments of protein families can be used to build a Potts Hamiltonian model of the sequence patterns which accurately predicts structural contacts. This observation paves the way to develop deeper connections between evolutionary fitness landscapes of entire protein families and the corresponding free energy landscapes which determine the conformational propensities of individual proteins. Using statistical energies determined from the Potts model and an alignment of 2896 PDB structures, we predict the propensity for particular kinase family proteins to assume a “DFG‐out” conformation implicated in the susceptibility of some kinases to type‐II inhibitors, and validate the predictions by comparison with the observed structural propensities of the corresponding proteins and experimental binding affinity data. We decompose the statistical energies to investigate which interactions contribute the most to the conformational preference for particular sequences and the corresponding proteins. We find that interactions involving the activation loop and the C‐helix and HRD motif are primarily responsible for stabilizing the DFG‐in state. This work illustrates how structural free energy landscapes and fitness landscapes of proteins can be used in an integrated way, and in the context of kinase family proteins, can potentially impact therapeutic design strategies.     

The R1275Q Neuroblastoma Mutant and Certain ATP-competitive Inhibitors Stabilize Alternative Activation Loop Conformations of Anaplastic Lymphoma Kinase

Anaplastic lymphoma kinase (ALK) is a receptor tyrosine kinase that, when genetically altered by mutation, amplification, chromosomal translocation or inversion, has been shown to play an oncogenic role in certain cancers. Small molecule inhibitors targeting the kinase activity of ALK have proven to be effective therapies in certain ALK-driven malignancies and one such inhibitor, crizotinib, is now approved for the treatment of EML4-ALK-driven, non-small cell lung cancer. In neuroblastoma, activating point mutations in the ALK kinase domain can drive disease progression, with the two most common mutations being F1174L and R1275Q. We report here crystal structures of the ALK kinase domain containing the F1174L and R1275Q mutations. Also included are crystal structures of ALK in complex with novel small molecule ALK inhibitors, including a classic type II inhibitor, that stabilize previously unobserved conformations of the ALK activation loop. Collectively, these structures illustrate a different series of activation loop conformations than has been observed in previous ALK crystal structures and provide insight into the activating nature of the R1275Q mutation. The novel active site topologies presented here may also aid the structure-based drug design of a new generation of ALK inhibitors.     

Crystallographic and solution studies of an activation loop mutant of the insulin receptor tyrosine kinase: insights into kinase mechanism

The tyrosine kinase domain of the insulin receptor is subject to autoinhibition in the unphosphorylated basal state via steric interactions involving the activation loop. A mutation in the activation loop designed to relieve autoinhibition, Asp-1161 --> Ala, substantially increases the ability of the unphosphorylated kinase to bind ATP. The crystal structure of this mutant in complex with an ATP analog has been determined at 2.4-A resolution. The structure shows that the active site is unobstructed, but the end of the activation loop is disordered and therefore the binding site for peptide substrates is not fully formed. In addition, Phe-1151 of the protein kinase-conserved DFG motif, at the beginning of the activation loop, hinders closure of the catalytic cleft and proper positioning of alpha-helix C for catalysis. These results, together with viscometric kinetic measurements, suggest that peptide substrate binding induces a reconfiguration of the unphosphorylated activation loop prior to the catalytic step. The crystallographic and solution studies provide new insights into the mechanism by which the activation loop controls phosphoryl transfer as catalyzed by the insulin receptor.     

A flexible-protein molecular docking study of the binding of ruthenium complex compounds to PIM1, GSK-3β, and CDK2/Cyclin A protein kinases

We employ ensemble docking simulations to characterize the interactions of two enantiomeric forms of a Ru-complex compound (1-R and 1-S) with three protein kinases, namely PIM1, GSK-3β, and CDK2/cyclin A. We show that our ensemble docking computational protocol adequately models the structural features of these interactions and discriminates between competing conformational clusters of ligand-bound protein structures. Using the determined X-ray crystal structure of PIM1 complexed to the compound 1-R as a control, we discuss the importance of including the protein flexibility inherent in the ensemble docking protocol, for the accuracy of the structure prediction of the bound state. A comparison of our ensemble docking results suggests that PIM1 and GSK-3β bind the two enantiomers in similar fashion, through two primary binding modes: conformation I, which is very similar to the conformation presented in the existing PIM1/compound 1-R crystal structure; conformation II, which represents a 180° flip about an axis through the NH group of the pyridocarbazole moiety, relative to conformation I. In contrast, the binding of the enantiomers to CDK2 is found to have a different structural profile including a suggested bound conformation, which lacks the conserved hydrogen bond between the kinase and the ligand (i.e., ATP, staurosporine, Ru-complex compound). The top scoring conformation of the inhibitor bound to CDK2 is not present among the top-scoring conformations of the inhibitor bound to either PIM1 or GSK-3β and vice-versa. Collectively, our results help provide atomic-level insights into inhibitor selectivity among the three kinases.

A Highly Selective Dual Insulin Receptor (IR)/Insulin-like Growth Factor 1 Receptor (IGF-1R) Inhibitor Derived from an Extracellular Signal-regulated Kinase (ERK) Inhibitor

Dual inhibitors of the closely related receptor tyrosine kinases insulin-like growth factor 1 receptor (IGF-1R) and insulin receptor (IR) are promising therapeutic agents in cancer. Here, we report an unusually selective class of dual inhibitors of IGF-1R and IR identified in a parallel screen of known kinase inhibitors against a panel of 300 human protein kinases. Biochemical and structural studies indicate that this class achieves its high selectivity by binding to the ATP-binding pocket of inactive, unphosphorylated IGF-1R/IR and stabilizing the activation loop in a native-like inactive conformation. One member of this compound family was originally reported as an inhibitor of the serine/threonine kinase ERK, a kinase that is distinct in the structure of its unphosphorylated/inactive form from IR/IGF-1R. Remarkably, this compound binds to the ATP-binding pocket of ERK in an entirely different conformation to that of IGF-1R/IR, explaining the potency against these two structurally distinct kinase families. These findings suggest a novel approach to polypharmacology in which two or more unrelated kinases are inhibited by a single compound that targets different conformations of each target kinase.     

Conformational dynamics of human protein kinase CK2α and its effect on function and inhibition

Protein kinase, casein kinase II (CK2), is ubiquitously expressed and highly conserved protein kinase that shows constitutive activity. It phosphorylates a diverse set of proteins and plays crucial role in several cellular processes. The catalytic subunit of this enzyme (CK2α) shows remarkable flexibility as evidenced in numerous crystal structures determined till now. Here, using analysis of multiple crystal structures and long timescale molecular dynamics simulations, we explore the conformational flexibility of CK2α. The enzyme shows considerably higher flexibility in the solution as compared to that observed in crystal structure ensemble. Multiple conformations of hinge region, located near the active site, were observed during the dynamics. We further observed that among these multiple conformations, the most populated conformational state was inadequately represented in the crystal structure ensemble. The catalytic spine, was found to be less dismantled in this state as compared to the “open” hinge/αD state crystal structures. The comparison of dynamics in unbound (Apo) state and inhibitor (CX4945) bound state exhibits inhibitor induced suppression in the overall dynamics of the enzyme. This is especially true for functionally important glycine‐rich loop above the active site. Together, this work gives novel insights into the dynamics of CK2α in solution and relates it to the function. This work also explains the effect of inhibitor on the dynamics of CK2α and paves way for development of better inhibitors.