Protein engineering of the colony-stimulating factor-1 receptor kinase domain for structural studies

A parallel approach to designing crystallization constructs for the c-FMS kinase domain was implemented, resulting in proteins suitable for structural studies. Sequence alignment and limited proteolysis were used to identify and eliminate unstructured and surface-exposed domains. A small library of chimeras was prepared in which the kinase insert domain of FMS was replaced with the kinase insert domain of previously crystallized receptor-tyrosine kinases. Characterization of the newly generated FMS constructs by enzymology and thermoshift assays demonstrated similar activities and compound binding to the FMS full-length cytoplasmic domain. Two chimeras were evaluated for crystallization in the presence and absence of a variety of ligands resulting in crystal structures, and leading to a successful structure-based drug design project for this important inflammation target.     

Inhibition of protein kinase CK2 by anthraquinone-related compounds. A structural insight

Protein kinases play key roles in signal transduction and therefore are among the most attractive targets for drug design. The pharmacological aptitude of protein kinase inhibitors is highlighted by the observation that various diseases with special reference to cancer are because of the abnormal expression/activity of individual kinases. The resolution of the three-dimensional structure of the target kinase in complex with inhibitors is often the starting point for the rational design of this kind of drugs, some of which are already in advanced clinical trial or even in clinical practice. Here we present and discuss three new crystal structures of ATP site-directed inhibitors in complex with "casein kinase-2" (CK2), a constitutively active protein kinase implicated in a variety of cellular functions and misfunctions. With the help of theoretical calculations, we disclose some key features underlying the inhibitory efficiency of anthraquinone derivatives, outlining three different binding modes into the active site. In particular, we show that a nitro group in a hydroxyanthraquinone scaffold decreases the inhibitory constants K(i) because of electron-withdrawing and resonance effects that enhance the polarization of hydroxylic substituents in paraposition.     

Effective handling of induced-fit motion in flexible docking

For structure-based drug design, where various ligand structures need to be docked to a target protein structure, a docking method that can handle conformational flexibility of not only the ligand, but also the protein, is indispensable. We have developed a simple and effective approach for dealing with the local induced-fit motion of the target protein, and implemented it in our docking tool, ADAM. Our approach efficiently combines the following two strategies: a vdW-offset grid in which the protein cavity is enlarged uniformly, and structure optimization allowing the motion of ligand and protein atoms. To examine the effectiveness of our approach, we performed docking validation studies, including redocking in 18 test cases and foreign-docking, in which various ligands from foreign crystal structures of complexes are docked into a target protein structure, in 22 cases (on five target proteins). With the original ADAM, the correct docking modes (RMSD < 2.0 A) were not present among the top 20 models in one case of redocking and four cases of foreign-docking. When the handling of induced-fit motion was implemented, the correct solutions were acquired in all 40 test cases. In foreign-docking on thymidine kinase, the correct docking modes were obtained as the top-ranked solutions for all 10 test ligands by our combinatorial approach, and this appears to be the best result ever reported with any docking tool. The results of docking validation have thus confirmed the effectiveness of our approach, which can provide reliable docking models even in the case of foreign-docking, where conformational change of the target protein cannot be ignored. We expect that this approach will contribute substantially to actual drug design, including virtual screening.     

Identifying protein construct variants with increased crystallization propensity--a case study

This study describes an efficient multiparallel automated workflow of cloning, expression, purification, and crystallization of a large set of construct variants for isolated protein domains aimed at structure determination by X-ray crystallography. This methodology is applied to MAPKAP kinase 2, a key enzyme in the inflammation pathway and thus an attractive drug target. The study reveals a distinct subset of truncation variants with improved crystallization properties. These constructs distinguish themselves by increased solubility and stability during a parallel automated multistep purification process including removal of the recombinant tag. High-throughput protein melting point analysis characterizes this subset of constructs as particularly thermostable. Both parallel purification screening and melting point determination clearly identify residue 364 as the optimal C terminus for the kinase domain. Moreover, all three constructs that ultimately crystallized feature this C terminus. At the N terminus, only three amino acids differentiate a noncrystallizing from a crystallizing construct. This study addresses the very common issues associated with difficult to crystallize proteins, those of solubility and stability, and the crucial importance of particular residues in the formation of crystal contacts. A methodology is suggested that includes biophysical measurements to efficiently identify and produce construct variants of isolated protein domains which exhibit higher crystallization propensity.     

Lining the pockets of kinases and phosphatases

The regulation of the activity of kinases and phosphatases is an essential aspect of intracellular signal transduction. Recently determined structures of AGC protein kinases, including isoforms of PKB, PKC, GRK and ROCK, indicate that occupancy of a hydrophobic pocket in the kinase N-lobe by a segment of the protein immediately C terminal to the kinase domain provides a mechanism for regulating kinase activity. In addition, crystal structures of Aurora-A and Aurora-B, which are closely related to AGC family kinases, in complex with their activators, TPX2 and INCENP, respectively, show how allosteric kinase activation is achieved by the binding of the activator protein to an equivalent hydrophobic pocket. Hence, regulation of kinase activity by analogous interactions is a shared regulatory mechanism of these kinases. Two crystal structures have explained the molecular basis of PKA anchoring through its regulatory subunits by members of the AKAP family of scaffold proteins. AKAPs can also interact directly with protein kinase and phosphatase catalytic domains. The crystal structure of the PP1 catalytic subunit in complex with the targeting subunit MYPT1 indicates that there is also scope for intimate phosphatase regulation by scaffold proteins.     

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.     

Analysis of water patterns in protein kinase binding sites

Deregulation of protein kinases is associated with numerous diseases, making them important targets for drug discovery. The majority of drugs target the catalytic site of these proteins, but due to the high level of similarity within the ATP binding sites of protein kinases, it is often difficult to achieve the required pharmacological selectivity. In this study, we describe the identification and subsequent analysis of water patterns in the ATP binding sites of 171 protein kinase structures, comprising 19 different kinases from various branches of the kinome, and demonstrate that structurally similar binding sites often have significantly different water patterns. We show that the observed variations in water patterns of different, but structurally similar kinases can be exploited in the structure-based design of potent and selective kinase inhibitors.     

Chaperone-assisted crystallography with DARPins

The structure of proteins that are difficult to crystallize can often be solved by forming a noncovalent complex with a helper protein--a crystallization "chaperone." Although several such applications have been described to date, their handling usually is still very laborious. A valuable addition to the present repertoire of binding proteins is the recently developed designed ankyrin repeat protein (DARPin) technology. DARPins are built based on the natural ankyrin repeat protein fold with randomized surface residue positions allowing specific binding to virtually any target protein. The broad potential of these binding proteins for X-ray crystallography is illustrated by five cocrystal structures that have been determined recently comprising target proteins from distinct families, namely a sugar binding protein, two kinases, a caspase, and a membrane protein. This article reviews the opportunities of this technology for structural biology and the structural aspects of the DARPin-protein complexes.     

Conservation of structural fluctuations in homologous protein kinases and its implications on functional sites

Our aim is to explore the similarities in structural fluctuations of homologous kinases. Gaussian Network Model based Normal Mode Analysis was performed on 73 active conformation structures in Ser/Thr/Tyr kinase superfamily. Categories of kinases with progressive evolutionary divergence, viz. (i) Same kinase with many crystal structures, (ii) Within‐Subfamily, (iii) Within‐Family, (iv) Within‐Group, and (v) Across‐Group, were analyzed. We identified a flexibility signature conserved in all kinases involving residues in and around the catalytic loop with consistent low‐magnitude fluctuations. However, the overall structural fluctuation profiles are conserved better in closely related kinases (Within‐Subfamily and Within‐family) than in distant ones (Within‐Group and Across‐Group). A substantial 65.4% of variation in flexibility was not accounted by variation in sequences or structures. Interestingly, we identified substructural residue‐wise fluctuation patterns characteristic of kinases of different categories. Specifically, we recognized statistically significant fluctuations unique to families of protein kinase A, cyclin‐dependent kinases, and nonreceptor tyrosine kinases. These fluctuation signatures localized to sites known to participate in protein‐protein interactions typical of these kinase families. We report for the first time that residues characterized by fluctuations unique to the group/family are involved in interactions specific to the group/family. As highlighted for Src family, local regions with differential fluctuations are proposed as attractive targets for drug design. Overall, our study underscores the importance of consideration of fluctuations, over and above sequence and structural features, in understanding the roles of sites characteristic of kinases. Proteins 2016; 84:957–978. © 2016 Wiley Periodicals, Inc.     

Crystal structure of the kinase domain of human vascular endothelial growth factor receptor 2: a key enzyme in angiogenesis.

BACKGROUND: Angiogenesis is involved in tumor growth, macular degeneration, retinopathy and other diseases. Vascular endothelial growth factor (VEGF) stimulates angiogenesis by binding to specific receptors (VEGFRs) on the surface of vascular endothelial cells. VEGFRs are receptor tyrosine kinases that, like the platelet-derived growth factor receptors (PDGFRs), contain a large insert within the kinase domain. RESULTS: We report here the generation, kinetic characterization, and 2.4 A crystal structure of the catalytic kinase domain of VEGF receptor 2 (VEGFR2). This protein construct, which lacks 50 central residues of the 68-residue kinase insert domain (KID), has comparable kinase activity to constructs containing the entire KID. The crystal structure, determined in an unliganded phosphorylated state, reveals an overall fold and catalytic residue positions similar to those observed in other tyrosine-kinase structures. The kinase activation loop, autophosphorylated on Y1059 prior to crystallization, is mostly disordered; however, a portion of it occupies a position inhibitory to substrate binding. The ends of the KID form a beta-like structure, not observed in other known tyrosine kinase structures, that packs near to the kinase C terminus. CONCLUSIONS: The majority of the VEGFR2 KID residues are not necessary for kinase activity. The unique structure observed for the ends of the KID may also occur in other PDGFR family members and may serve to properly orient the KID for signal transduction. This VEGFR2 kinase structure provides a target for design of selective anti-angiogenic therapeutic agents.     

Imprint of evolutionary conservation and protein structure variation on the binding function of protein tyrosine kinases

MOTIVATION: According to the models of divergent molecular evolution, the evolvability of new protein function may depend on the induction of new phenotypic traits by a small number of mutations of the binding site residues. Evolutionary relationships between protein kinases are often employed to infer inhibitor binding profiles from sequence analysis. However, protein kinases binding profiles may display inhibitor selectivity within a given kinase subfamily, while exhibiting cross-activity between kinases that are phylogenetically remote from the prime target. The emerging insights into kinase function and evolution combined with a rapidly growing number of publically available crystal structures of protein kinases complexes have motivated structural bioinformatics analysis of sequence-structure relationships in determining the binding function of protein tyrosine kinases. RESULTS: In silico profiling of Imatinib mesylate and PD-173955 kinase inhibitors with protein tyrosine kinases is conducted on kinome scale by using evolutionary analysis and fingerprinting inhibitor-protein interactions with the panel of all publically available protein tyrosine kinases crystal structures. We have found that sequence plasticity of the binding site residues alone may not be sufficient to enable protein tyrosine kinases to readily evolve novel binding activities with inhibitors. While evolutionary signal derived solely from the tyrosine kinase sequence conservation can not be readily translated into the ligand binding phenotype, the proposed structural bioinformatics analysis can discriminate a functionally relevant kinase binding signal from a simple phylogenetic relationship. The results of this work reveal that protein conformational diversity is intimately linked with sequence plasticity of the binding site residues in achieving functional adaptability of protein kinases towards specific drug binding. This study offers a plausible molecular rationale to the experimental binding profiles of the studied kinase inhibitors and provides a theoretical basis for constructing functionally relevant kinase binding trees.     

Exploring c-Met kinase flexibility by sampling and clustering its conformational space

It is now widely recognized that the flexibility of both partners has to be considered in molecular docking studies. However, the question how to handle the best the huge computational complexity of exploring the protein binding site landscape is still a matter of debate. Here we investigate the flexibility of c-Met kinase as a test case for comparing several simulation methods. The c-Met kinase catalytic site is an interesting target for anticancer drug design. In particular, it harbors an unusual plasticity compared with other kinases ATP binding sites. Exploiting this feature may eventually lead to the discovery of new anticancer agents with exquisite specificity. We present in this article an extensive investigation of c-Met kinase conformational space using large-scale computational simulations in order to extend the knowledge already gathered from available X-ray structures. In the process, we compare the relevance of different strategies for modeling and injecting receptor flexibility information into early stage in silico structure-based drug discovery pipeline. The results presented here are currently being exploited in on-going virtual screening investigations on c-Met.     

Improved amino acid flexibility parameters

Protein molecules exhibit varying degrees of flexibility throughout their three-dimensional structures, with some segments showing little mobility while others may be so disordered as to be unresolvable by techniques such as X-ray crystallography. Atomic displacement parameters, or B-factors, from X-ray crystallographic studies give an experimentally determined indication of the degree of mobility in a protein structure. To provide better estimators of amino acid flexibility, we have examined B-factors from a large set of high-resolution crystal structures. Because of the differences among structures, it is necessary to normalize the B-factors. However, many proteins have segments of unusually high mobility, which must be accounted for before normalization can be performed. Accordingly, a median-based method from quality control studies was used to identify outliers. After removal of outliers from, and normalization of, each protein chain, the B-factors were collected for each amino acid in the set. It was found that the distribution of normalized B-factors followed a Gumbel, or extreme value distribution, and the location parameter, or mode, of this distribution was used as an estimator of flexibility for the amino acid. These new parameters have a higher correlation with experimentally determined B-factors than parameters from earlier methods.     

Target flexibility in molecular recognition

Induced-fit effects are well known in the binding of small molecules to proteins and other macromolecular targets. Among other targets, protein kinases are particularly flexible proteins, so that such effects should be considered in attempts at structure-based inhibitor design for kinase targets. This paper outlines some recent progress in methods for including target flexibility in computational studies of molecular recognition. A focus is the "relaxed complex method," in which ligands are docked to an ensemble of conformations of the target, and the best complexes are re-scored to provide predictions of optimal binding geometries. Early applications of this method have suggested a new approach to the development of inhibitors of HIV-1 Integrase.     

Conformational Transitions upon Ligand Binding: Holo-Structure Prediction from Apo Conformations

Biological function of proteins is frequently associated with the formation of complexes with small-molecule ligands. Experimental structure determination of such complexes at atomic resolution, however, can be time-consuming and costly. Computational methods for structure prediction of protein/ligand complexes, particularly docking, are as yet restricted by their limited consideration of receptor flexibility, rendering them not applicable for predicting protein/ligand complexes if large conformational changes of the receptor upon ligand binding are involved. Accurate receptor models in the ligand-bound state (holo structures), however, are a prerequisite for successful structure-based drug design. Hence, if only an unbound (apo) structure is available distinct from the ligand-bound conformation, structure-based drug design is severely limited. We present a method to predict the structure of protein/ligand complexes based solely on the apo structure, the ligand and the radius of gyration of the holo structure. The method is applied to ten cases in which proteins undergo structural rearrangements of up to 7.1 Å backbone RMSD upon ligand binding. In all cases, receptor models within 1.6 Å backbone RMSD to the target were predicted and close-to-native ligand binding poses were obtained for 8 of 10 cases in the top-ranked complex models. A protocol is presented that is expected to enable structure modeling of protein/ligand complexes and structure-based drug design for cases where crystal structures of ligand-bound conformations are not available.     

Abl kinase constructs expressed in bacteria: facilitation of structural and functional studies including segmental labeling by expressed protein ligation

A great portion of tyrosine kinases are involved in cell development and their structural alteration is intimately involved in associated pathologies of development and oncology. These kinases are one of the major groups of targets under investigation for molecular therapeutics. To carry out biochemical and structural biological studies on these kinases, economical production of their purified forms is highly desirable. However over-expressing tyrosine kinases as recombinant forms in bacterial systems and their purification is a significant challenge. Abelson kinase (Abl) has previously been expressed on a large scale to facilitate X-ray crystallography and NMR structure studies mainly in baculovirus infected insect cells. Even though success has been achieved in expression of soluble tyrosine kinases in E. coli with chaperones to improve correct folding, low expression levels of kinases are intrinsic in such systems because of diversion of resources to produce chaperones. Here we present a straightforward method to express and purify isolated Abl kinase domain and SH3-SH2-kinase multi-domain structures. The expressed Abl protein retains its correct folding and biological function. The yield of soluble protein is in a several mg L(-1) range in minimal media. Furthermore we demonstrate that segmental isotopic labelling using expressed protein ligation can be achieved using bacterial expressed Abl kinase domain constructs, which is especially useful in NMR structure-activity studies.     

Crystal structures of ASK1‐inhibtor complexes provide a platform for structure‐based drug design

ASK1, a member of the MAPK Kinase Kinase family of proteins has been shown to play a key role in cancer, neurodegeneration and cardiovascular diseases and is emerging as a possible drug target. Here we describe a ‘replacement‐soaking’ method that has enabled the high‐throughput X‐ray structure determination of ASK1/ligand complexes. Comparison of the X‐ray structures of five ASK1/ligand complexes from 3 different chemotypes illustrates that the ASK1 ATP binding site is able to accommodate a range of chemical diversity and different binding modes. The replacement‐soaking system is also able to tolerate some protein flexibility. This crystal system provides a robust platform for ASK1/ligand structure determination and future structure based drug design.     

Insights into the conformational flexibility of Bruton's tyrosine kinase from multiple ligand complex structures

Bruton's tyrosine kinase (BTK) plays a key role in B cell receptor signaling and is considered a promising drug target for lymphoma and inflammatory diseases. We have determined the X-ray crystal structures of BTK kinase domain in complex with six inhibitors from distinct chemical classes. Five different BTK protein conformations are stabilized by the bound inhibitors, providing insights into the structural flexibility of the Gly-rich loop, helix C, the DFG sequence, and activation loop. The conformational changes occur independent of activation loop phosphorylation and do not correlate with the structurally unchanged WEI motif in the Src homology 2-kinase domain linker. Two novel activation loop conformations and an atypical DFG conformation are observed representing unique inactive states of BTK. Two regions within the activation loop are shown to structurally transform between 3(10)- and α-helices, one of which collapses into the adenosine-5'-triphosphate binding pocket. The first crystal structure of a Tec kinase family member in the pharmacologically important DFG-out conformation and bound to a type II kinase inhibitor is described. The different protein conformations observed provide insights into the structural flexibility of BTK, the molecular basis of its regulation, and the structure-based design of specific inhibitors.     

A structural comparison of inhibitor binding to PKB, PKA and PKA-PKB chimera

Although the crystal structure of the anti-cancer target protein kinase B (PKBbeta/Akt-2) has been useful in guiding inhibitor design, the closely related kinase PKA has generally been used as a structural mimic due to its facile crystallization with a range of ligands. The use of PKB-inhibitor crystallography would bring important benefits, including a more rigorous understanding of factors dictating PKA/PKB selectivity, and the opportunity to validate the utility of PKA-based surrogates. We present a "back-soaking" method for obtaining PKBbeta-ligand crystal structures, and provide a structural comparison of inhibitor binding to PKB, PKA, and PKA-PKB chimera. One inhibitor presented here exhibits no PKB/PKA selectivity, and the compound adopts a similar binding mode in all three systems. By contrast, the PKB-selective inhibitor A-443654 adopts a conformation in PKB and PKA-PKB that differs from that with PKA. We provide a structural explanation for this difference, and highlight the ability of PKA-PKB to mimic the true PKB binding mode in this case.