The role of conserved water molecules in the catalytic domain of protein kinases

Protein kinases are essential signaling molecules with a characteristic bilobal shape that has been studied for over 15 years. Despite the number of crystal structures available, little study has been directed away from the prototypical functional elements of the kinase domain. We have performed a structural alignment of 13 active‐conformation kinases and discovered the presence of six water molecules that occur in conserved locations across this group of diverse kinases. Molecular dynamics simulations demonstrated that these waters confer a great deal of stability to their local environment and to a key catalytic residue. Our results highlight the importance of novel elements within the greater kinase family and suggest that conserved water molecules are necessary for efficient kinase function. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

Hydration energy landscape of the active site cavity in cytochrome P450cam

Hydration of protein cavities influences protein stability, dynamics, and function. Protein active sites usually contain water molecules that, upon ligand binding, are either displaced into bulk solvent or retained to mediate protein–ligand interactions. The contribution of water molecules to ligand binding must be accounted for to compute accurate values of binding affinities. This requires estimation of the extent of hydration of the binding site. However, it is often difficult to identify the water molecules involved in the binding process when ligands bind on the surface of a protein. Cytochrome P450cam is, therefore, an ideal model system because its substrate binds in a buried active site, displacing partially disordered solvent, and the protein is well characterized experimentally. We calculated the free energy differences for having five to eight water molecules in the active site cavity of the unliganded enzyme from molecular dynamics simulations by thermodynamic integration employing a three-stage perturbation scheme. The computed free energy differences between the hydration states are small (within 12 kJ mol⁻¹) but distinct. Consistent with the crystallographic determination and studies employing hydrostatic pressure, we calculated that, although ten water molecules could in principle occupy the volume of the active site, occupation by five to six water molecules is thermodynamically most favorable. Proteins 32:381–396, 1998. © 1998 Wiley-Liss, Inc.     

Evolutionary variation and adaptation in a conserved protein kinase allosteric network: Implications for inhibitor design

The activation of protein kinases involves conformational changes in key functional regions of the kinase domain, a detailed understanding of which is essential for the design of selective protein kinase inhibitors. Through statistical analysis of protein kinase sequences and crystal structures from diverse organisms, we recently proposed that the activation of protein kinases involves a hidden strain switch in the catalytic loop. Specifically, we demonstrated that the backbone torsion-angles of residues in the catalytic loop switch from a “relaxed” to “strained” conformation upon kinase activation and the strained geometry results in a network of hydrogen bonds involving conserved non-catalytic residues in the ATP and substrate binding lobes. Here, we further explore this activation mechanism by analyzing families that lack the canonical hydrogen bonding interactions with the strained backbone. We find that alternative mechanisms have evolved to maintain catalytic loop strain. In PIM kinase, for example, two water molecules account for the lack of a conserved aspartate in the substrate binding by hydrogen bonds to the strained backbone. We discuss the relevance of these findings in the design of family-specific allosteric inhibitors, and in predicting the structural and functional impact of cancer mutations that alter the strain associated hydrogen bonding network. This article is part of a Special Issue entitled: Inhibitors of Protein Kinases (2012).     

Dynamics of water molecules buried in cavities of apolipoprotein E studied by molecular dynamics simulations and continuum electrostatic calculations

Molecular dynamics (MD) simulations of several nanoseconds each were used to monitor the dynamic behavior of the five crystal water molecules buried in the interior of the N-terminal domain of apolipoprotein E. These crystal water molecules are fairly well conserved in several apolipoprotein E structures, suggesting that they are not an artifact of the crystal and that they may have a structural and/or functional role for the protein. All five buried crystal water molecules leave the protein interior in the course of the longest simulations and exchange with water molecules from the bulk. The free energies of binding evaluated from the electrostatic binding free energy computed using a continuum model and estimates of the binding entropy changes represent shallow minima. The corresponding calculated residence times of the buried water molecules range from tens of picoseconds to hundreds of nanoseconds, which denote rather short times as for buried water molecules. Several water exchanges monitored in the simulations show that water molecules along the exit/entrance pathway use a relay of H bonds primarily formed with charged residues which helps either the exit or the entrance from or into the buried site. The exit/entrance of water molecules from/into the sites is permitted essentially by local motions of, at most, two side chains, indicating that, in these cases, complex correlated atomic motions are not needed to open the buried site toward the surface of the protein. This provides a possible explanation for the short residence times.     

Catalytic independent functions of a protein kinase as revealed by a kinase-dead mutant: study of the Lys72His mutant of cAMP-dependent kinase

A highly conserved lysine in subdomain II is required for high catalytic activity among the protein kinases. This lysine interacts directly with ATP and mutation of this residue leads to a classical "kinase-dead" mutant. This study describes the biophysical and functional properties of a kinase-dead mutant of cAMP-dependent kinase where Lys72 was replaced with His. Although the mutant protein is less stable than the wild-type catalytic subunit, it is fully capable of binding ATP. The results highlight the effect of the mutation on stability and overall organization of the protein, especially the small lobe. Phosphorylation of the activation loop by a heterologous kinase, 3-phosphoinositide-dependent protein kinase-1 (PDK-1) also contributes dramatically to the global organization of the entire active site region. Deuterium-exchange mass spectrometry (DXMS) indicates a concerted stabilization of the entire active site following the addition of this single phosphate to the activation loop. Furthermore the mutant C-subunit is capable of binding both the type I and II regulatory subunits, but only after phosphorylation of the activation loop. This highlights the role of the large lobe as a scaffold for the regulatory subunits independent of catalytic competency and suggests that kinase dead members of the protein kinase superfamily may still have other important biological roles although they lack catalytic activity.     

DRoP: Automated detection of conserved solvent-binding sites on proteins

Water and ligand binding play critical roles in the structure and function of proteins, yet their binding sites and significance are difficult to predict a priori. Multiple solvent crystal structures (MSCS) is a method where several X-ray crystal structures are solved, each in a unique solvent environment, with organic molecules that serve as probes of the protein surface for sites evolved to bind ligands, while the first hydration shell is essentially maintained. When superimposed, these structures contain a vast amount of information regarding hot spots of protein-protein or protein-ligand interactions, as well as conserved water-binding sites retained with the change in solvent properties. Optimized mining of this information requires reliable structural data and a consistent, objective analysis tool. Detection of related solvent positions (DRoP) was developed to automatically organize and rank the water or small organic molecule binding sites within a given set of structures. It is a flexible tool that can also be used in conserved water analysis given multiple structures of any protein independent of the MSCS method. The DRoP output is an HTML format list of the solvent sites ordered by conservation rank in its population within the set of structures, along with renumbered and recolored PDB files for visualization and facile analysis. Here, we present a previously unpublished set of MSCS structures of bovine pancreatic ribonuclease A (RNase A) and use it together with published structures to illustrate the capabilities of DRoP.     

Thermodynamic analysis of water molecules at the surface of proteins and applications to binding site prediction and characterization

Water plays an essential role in determining the structure and function of all biological systems. Recent methodological advances allow for an accurate and efficient estimation of the thermodynamic properties of water molecules at the surface of proteins. In this work, we characterize these thermodynamic properties and relate them to various structural and functional characteristics of the protein. We find that high-energy hydration sites often exist near protein motifs typically characterized as hydrophilic, such as backbone amide groups. We also find that waters around alpha helices and beta sheets tend to be less stable than waters around loops. Furthermore, we find no significant correlation between the hydration site-free energy and the solvent accessible surface area of the site. In addition, we find that the distribution of high-energy hydration sites on the protein surface can be used to identify the location of binding sites and that binding sites of druggable targets tend to have a greater density of thermodynamically unstable hydration sites. Using this information, we characterize the FKBP12 protein and show good agreement between fragment screening hit rates from NMR spectroscopy and hydration site energetics. Finally, we show that water molecules observed in crystal structures are less stable on average than bulk water as a consequence of the high degree of spatial localization, thereby resulting in a significant loss in entropy. These findings should help to better understand the characteristics of waters at the surface of proteins and are expected to lead to insights that can guide structure-based drug design efforts.     

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.     

Evolutionary constraints associated with functional specificity of the CMGC protein kinases MAPK, CDK, GSK, SRPK, DYRK, and CK2alpha

Amino acid residues associated with functional specificity of cyclin-dependent kinases (CDKs), mitogen-activated protein kinases (MAPKs), glycogen synthase kinases (GSKs), and CDK-like kinases (CLKs), which are collectively termed the CMGC group, were identified by categorizing and quantifying the selective constraints acting upon these proteins during evolution. Many constraints specific to CMGC kinases correspond to residues between the N-terminal end of the activation segment and a CMGC-conserved insert segment associated with coprotein binding. The strongest such constraint is imposed on a "CMGC-arginine" near the substrate phosphorylation site with a side chain that plays a role both in substrate recognition and in kinase activation. Two nearby buried waters, which are also present in non-CMGC kinases, typically position the main chain of this arginine relative to the catalytic loop. These and other CMGC-specific features suggest a structural linkage between coprotein binding, substrate recognition, and kinase activation. Constraints specific to individual subfamilies point to mechanisms for CMGC kinase specialization. Within casein kinase 2alpha (CK2alpha), for example, the binding of one of the buried waters appears prohibited by the side chain of a leucine that is highly conserved within CK2alpha and that, along with substitution of lysine for the CMGC-arginine, may contribute to the broad substrate specificity of CK2alpha by relaxing characteristically conserved, precise interactions near the active site. This leucine is replaced by a conserved isoleucine or valine in other CMGC kinases, thereby illustrating the potential functional significance of subtle amino acid substitutions. Analysis of other CMGC kinases similarly suggests candidate family-specific residues for experimental follow-up.     

Hydration Structures of the Human Protein Kinase CK2α Clarified by Joint Neutron and X-ray Crystallography

Casein kinase 2 (CK2) has broad phosphorylation activity against various regulatory proteins, which are important survival factors in eukaryotic cells. To clarify the hydration structure and catalytic mechanism of CK2, we determined the crystal structure of the alpha subunit of human CK2 containing hydrogen and deuterium atoms using joint neutron (1.9 Å resolution) and X-ray (1.1 Å resolution) crystallography. The analysis revealed the structure of conserved water molecules at the active site and a long potential hydrogen bonding network originating from the catalytic Asp156 that is well known to enhance the nucleophilicity of the substrate OH group to the γ-phospho group of ATP by proton elimination. His148 and Asp214 conserved in the protein kinase family are located in the middle of the network. The water molecule forming a hydrogen bond with Asp214 appears to be deformed. In addition, mutational analysis of His148 in CK2 showed significant reductions by 40%–75% in the catalytic efficiency with similar affinity for ATP. Likewise, remarkable reductions to less than 5% were shown by corresponding mutations on His131 in death-associated protein kinase 1, which belongs to a group different from that of CK2. These findings shed new light on the catalytic mechanism of protein kinases in which the hydrogen bond network through the C-terminal domain may assist the general base catalyst to extract a proton with a link to the bulk solvent via intermediates of a pair of residues.     

A "solvated rotamer" approach to modeling water-mediated hydrogen bonds at protein-protein interfaces

Water-mediated hydrogen bonds play critical roles at protein-protein and protein-nucleic acid interfaces, and the interactions formed by discrete water molecules cannot be captured using continuum solvent models. We describe a simple model for the energetics of water-mediated hydrogen bonds, and show that, together with knowledge of the positions of buried water molecules observed in X-ray crystal structures, the model improves the prediction of free-energy changes upon mutation at protein-protein interfaces, and the recovery of native amino acid sequences in protein interface design calculations. We then describe a "solvated rotamer" approach to efficiently predict the positions of water molecules, at protein-protein interfaces and in monomeric proteins, that is compatible with widely used rotamer-based side-chain packing and protein design algorithms. Finally, we examine the extent to which the predicted water molecules can be used to improve prediction of amino acid identities and protein-protein interface stability, and discuss avenues for overcoming current limitations of the approach.     

Evolved to be active: sulfate ions define substrate recognition sites of CK2alpha and emphasise its exceptional role within the CMGC family of eukaryotic protein kinases

CK2alpha is the catalytic subunit of protein kinase CK2 and a member of the CMGC family of eukaryotic protein kinases like the cyclin-dependent kinases, the MAP kinases and glycogen-synthase kinase 3. We present here a 1.6 A resolution crystal structure of a fully active C-terminal deletion mutant of human CK2alpha liganded by two sulfate ions, and we compare this structure systematically with representative structures of related CMGC kinases. The two sulfate anions occupy binding pockets at the activation segment and provide the structural basis of the acidic consensus sequence S/T-D/E-X-D/E that governs substrate recognition by CK2. The anion binding sites are conserved among those CMGC kinases. In most cases they are neutralized by phosphorylation of a neighbouring threonine or tyrosine side-chain, which triggers conformational changes for regulatory purposes. CK2alpha, however, lacks both phosphorylation sites at the activation segment and structural plasticity. Here the anion binding sites are functionally changed from regulation to substrate recognition. These findings underline the exceptional role of CK2alpha as a constitutively active enzyme within a family of strictly controlled protein kinases.     

A novel lipid binding site formed by the MAP kinase insert in p38 alpha

The p38 mitogen-activated protein (MAP) kinases function as signaling molecules essential for many cellular processes, particularly mediating stress response. The activity of p38 MAP kinases is meticulously regulated to reach the desired cellular phenotype. Several alternative activation and attenuation mechanisms have been characterized recently which include new phosphorylation sites. Here we present the crystal structure of p38α MAP kinase in complex with n-octyl-β-glucopyranoside detergent. The complex unveils a novel lipid-binding site formed by a local conformational change of the MAP kinase insert. This binding is the first attribution for a possible role of the MAP kinase insert in p38. The binding site can accommodate a large selection of lipidic molecules. In addition, we also show via biophysical methods that arachidonic acid and its derivatives bind p38α in vitro. Based on our analysis we propose that the binding of lipids could fine-tune p38α catalytic activity towards a preferred phenotype.     

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.     

Structural and dynamic roles of permanent water molecules in ligand molecular recognition by chicken liver bile acid binding protein

Chicken liver bile acid binding protein (cL-BABP) crystallizes with water molecules in its binding site. To obtain insights on the role of internal water, we performed two 100 ns molecular dynamics (MD) simulations in explicit solvent for cL-BABP, as apo form and as a complex with two molecules of cholic acid, and analyzed in detail the dynamics properties of all water molecules. The diffusion coefficients of the more persistent internal water molecules are significantly different from the bulk, but similar between the two protein forms. A different number of molecules and a different organization are observed for apo- and holo-cL-BABP. Most water molecules identified in the binding site of the apo-crystal diffuse to the bulk during the simulation. In contrast, almost all the internal waters of the holo-crystal maintain the same interactions with internal sidechains and ligands, which suggests they have a relevant role in protein-ligand molecular recognition. Only in the presence of these water molecules we were able to reproduce, by a classical molecular docking approach, the structure of the complex cL-BABP::cholic acid with a low ligand root mean square deviation (RMSD) with respect to its reference positioning. Literature data reported a conserved pattern of hydrogen bonds between a single water molecule and three amino acid residues of the binding site in a series of crystallized FABP. In cL-BABP, the interactions between this conserved water molecule and the three residues are present in the crystal of both apo- and holo-cL-BABP but are lost immediately after the start of molecular dynamics. Copyright (c) 2008 John Wiley & Sons, Ltd.     

Buried water molecules in helical transmembrane proteins

Buried water molecules (having no contact with bulk solvent) in 30 helical transmembrane (TM) protein structures were identified. The average amount of buried water in helical TM proteins is about the same as for all water-soluble (WS) proteins, but it is greater than the average for helical WS proteins. Buried waters in TM proteins make more polar contacts, and are more frequently found contacting helices than in WS proteins. The distribution of the buried water binding sites across the membrane profile shows that the sites to some extent reflect protein function. There is also evidence for asymmetry of the sites, with more in the extracellular half of the membrane. Many of the buried water contact sites are conserved across families of proteins, including family members having different functions. This suggests that at least some buried waters play a role in structural stabilization. Disease-causing mutations, which are known to result in misfolded TM proteins, occur at buried water contact sites at a higher than random frequency, which also supports a stabilizing role for buried water molecules.     

Activation of cyclic AMP-dependent protein kinase isoenzymes: studies using specific antisera

A novel method for rapidly determining the amount and degree of association-dissociation of the Type I and Type II cAMP-dependent protein kinases has been developed and validated. Antibodies directed against the regulatory subunits of Type I and Type II cAMP-dependent protein kinases were used. The antibodies formed complexes with holoenzymes and regulatory subunits which were precipitated by goat anti-rabbit IgG (immunoglobulin G). These complexes bound [3H]cAMP with an apparent Kb of 20 nM for protein kinase I and 80 nM for protein kinase II. Immunoprecipitated protein kinases I and II were catalytically active when incubated with cAMP, [gamma-32P]ATP, and histone H2B. When mixtures of the two kinase isoenzymes or cytosol were incubated with various amounts of [3H]cAMP and the isoenzymes were separated by precipitation with antisera specific for each isoenzyme, the amount of [3H]cAMP associated with immunoprecipitates was proportional to the concentration of [3H]cAMP. In contrast, the catalytic activity that was immunoprecipitated varied inversely with the concentration of [3H]cAMP, showing that the activation of protein kinase could be assessed by the disappearance of catalytic activity from the immunoprecipitates. In the absence of MgATP protein kinase I was activated by a 10-fold lower concentration of cAMP than protein kinase II. However, when MgATP was added to the incubation, there was no significant difference in the binding of [3H]cAMP or dissociation of catalytic subunits of the two isoenzymes. The anti-R antibodies were also used to rapidly quantitate the concentration of regulatory subunits and the relative ratio of protein kinases I and II in tissue cytosols.     

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.     

Structural Characteristics in Protein Hydration Investigated by Cryogenic X-ray Crystal Structure Analyses

Cryogenic X-ray crystallography has heen applied to investigate thehydration structures of proteins. The amount of hydration water moleculesidentified at cryogenic temperature is more than twice those at ambienttemperature, and the structural models of proteins with a lot of hydrationwater molecules have provided much information to elucidate the static anddynamical characteristics of hydration structures of proteins. On proteinsurface, hydration water molecules distribute non-randomly and stillretain the tetrahedral hydrogen-bond geometry as well as in bulk solvent.In addition, water molecules form clathrate-like arrangements to cover thehydrophobic residues exposed to solvent. The standard interaction geometryenables the three-dimensional extension of hydrogen-bond networks aroundprotein molecules and, simultaneously, ensures the concerted reorganizationof hydration structures during the dynamical motion of proteins at work.The hydration structure analyses at cryogenic temperatures may contributeto understanding physical principles governing the dynamics of `molecularmachines' in aqueous environment.