Protein-protein interactions in the allosteric regulation of protein kinases

Perspectives on Computational Enzyme Modeling: From Mechanisms to Design and Drug Development

Understanding enzyme mechanisms is essential for unraveling the complex molecular machinery of life. In this review, we survey the field of computational enzymology, highlighting key principles governing enzyme mechanisms and discussing ongoing challenges and promising advances. Over the years, computer simulations have become indispensable in the study of enzyme mechanisms, with the integration of experimental and computational exploration now established as a holistic approach to gain deep insights into enzymatic catalysis. Numerous studies have demonstrated the power of computer simulations in characterizing reaction pathways, transition states, substrate selectivity, product distribution, and dynamic conformational changes for various enzymes. Nevertheless, significant challenges remain in investigating the mechanisms of complex multistep reactions, large-scale conformational changes, and allosteric regulation. Beyond mechanistic studies, computational enzyme modeling has emerged as an essential tool for computer-aided enzyme design and the rational discovery of covalent drugs for targeted therapies. Overall, enzyme design/engineering and covalent drug development can greatly benefit from our understanding of the detailed mechanisms of enzymes, such as protein dynamics, entropy contributions, and allostery, as revealed by computational studies. Such a convergence of different research approaches is expected to continue, creating synergies in enzyme research. This review, by outlining the ever-expanding field of enzyme research, aims to provide guidance for future research directions and facilitate new developments in this important and evolving field.

Fenretinide induces apoptosis and synergises the apoptosis inducing effect of gemcitabine through inhibition of key signalling molecules involved in A549 cell survival in in silico and in vitro analyses

Fenretinide is a synthetic retinoid compound, which induces apoptosis via generating reactive oxygen species (ROS) and modulating PI3K/Akt/mTOR signalling pathway. We hypothesise that fenretinide's mechanism of action in triggering apoptosis may involve other targets, beside mTOR signalling pathway and it may augment apoptosis inducing effects of chemotherapeutic drugs in lung cancer. Time-lapse microscopy and Western blotting were used to evaluate apoptosis and apoptotic marker cleaved-Caspase 3 in A549 cells. Relative levels of protein phosphorylation and ROS were quantified by Human Phospho-Kinase Array Kit and CellROX® Green Reagent, respectively. Docking and simulation analyses of proteins and fenretinide interactions were identified and visualised by Discovery Studio Visualizer and AutoDock Vina software. Our results showed that fenretinide induced apoptosis in a dose dependant manner and combinations of fenretinide (5 μg/mL) and gemcitabine (1, 2, 4, 8 and 16 μg/mL) synergistically enhanced apoptosis in A549 cells. Fenretinide caused significant increase of cleaved-Caspase 3, de-phosphorylated p-S473 of Akt and failed to inhibit mTORC1 downstream targets. In silico results revealed that Akt required the lowest energy (-10.2 kcal/mol) to interact with fenretinide in comparison with other proteins. In conclusion, Akt may be exploited as a good target for induction of apoptosis in A549 cells and fenretinide has great potentials to fulfil this task. The mechanism by which fenretinide boosts the apoptosis inducing effects of gemcitabine, which is likely expected to be via inhibiting mTORC2 downstream targets. However, docking investigation revealed that fenretinide lacks specificity as it may also interact with several secondary targets beside Akt.

Insights from the protein interaction Universe of the multifunctional "Goldilocks" kinase DYRK1A

Human Dual specificity tyrosine (Y)-Regulated Kinase 1A (DYRK1A) is encoded by a dosage-dependent gene located in the Down syndrome critical region of human chromosome 21. The known substrates of DYRK1A include proteins involved in transcription, cell cycle control, DNA repair and other processes. However, the function and regulation of this kinase is not fully understood, and the current knowledge does not fully explain the dosage-dependent function of this kinase. Several recent proteomic studies identified DYRK1A interacting proteins in several human cell lines. Interestingly, several of known protein substrates of DYRK1A were undetectable in these studies, likely due to a transient nature of the kinase-substrate interaction. It is possible that the stronger-binding DYRK1A interacting proteins, many of which are poorly characterized, are involved in regulatory functions by recruiting DYRK1A to the specific subcellular compartments or distinct signaling pathways. Better understanding of these DYRK1A-interacting proteins could help to decode the cellular processes regulated by this important protein kinase during embryonic development and in the adult organism. Here, we review the current knowledge of the biochemical and functional characterization of the DYRK1A protein-protein interaction network and discuss its involvement in human disease.

The fitness cost of spurious phosphorylation

The fidelity of signal transduction requires the binding of regulatory molecules to their cognate targets. However, the crowded cell interior risks off-target interactions between proteins that are functionally unrelated. How such off-target interactions impact fitness is not generally known, but quantifying this is required to understand the constraints faced by cell systems as they evolve. Here, we use the model organism S. cerevisiae to inducibly express tyrosine kinases. Because yeast lacks bona fide tyrosine kinases, most of the resulting tyrosine phosphorylation is spurious. This provides a suitable system to measure the impact of artificial protein interactions on fitness. We engineered 44 yeast strains each expressing a tyrosine kinase, and quantitatively analysed their phosphoproteomes. This analysis resulted in ~30,000 phosphosites mapping to ~3,500 proteins. Examination of the fitness costs in each strain revealed a strong correlation between the number of spurious pY sites and decreased growth. Moreover, the analysis of pY effects on protein structure and on protein function revealed over 1000 pY events that we predict to be deleterious. However, we also find that a large number of the spurious pY sites have a negligible effect on fitness, possibly because of their low stoichiometry. This result is consistent with our evolutionary analyses demonstrating a lack of phosphotyrosine counter-selection in species with bona fide tyrosine kinases. Taken together, our results suggest that, alongside the risk for toxicity, the cell can tolerate a large degree of non-functional crosstalk as interaction networks evolve.

Computational analysis of regulatory regions in human protein kinases

Eukaryotic proteins often feature modular domain structures comprising globular domains that are connected by linker regions and intrinsically disordered regions that may contain important functional motifs. The intramolecular interactions of globular domains and nonglobular regions can play critical roles in different aspects of protein function. However, studying these interactions and their regulatory roles can be challenging due to the flexibility of nonglobular regions, the long insertions separating interacting modules, and the transient nature of some interactions. Obtaining the experimental structures of multiple domains and functional regions is more difficult than determining the structures of individual globular domains. High-quality structural models generated by AlphaFold offer a unique opportunity to study intramolecular interactions in eukaryotic proteins. In this study, we systematically explored intramolecular interactions between human protein kinase domains (KDs) and potential regulatory regions, including globular domains, N- and C-terminal tails, long insertions, and distal nonglobular regions. Our analysis identified intramolecular interactions between human KDs and 35 different types of globular domains, exhibiting a variety of interaction modes that could contribute to orthosteric or allosteric regulation of kinase activity. We also identified prevalent interactions between human KDs and their flanking regions (N- and C-terminal tails). These interactions exhibit group-specific characteristics and can vary within each specific kinase group. Although long-range interactions between KDs and nonglobular regions are relatively rare, structural details of these interactions offer new insights into the regulation mechanisms of several kinases, such as HASPIN, MAPK7, MAPK15, and SIK1B.

Engineered Allosteric Regulation of Protein Function

Allosteric regulation of proteins has been utilized to study various aspects of cell signaling, from unicellular events to organism-wide phenotypes. However, traditional methods of allosteric regulation, such as constitutively active mutants and inhibitors, lack tight spatiotemporal control. This often leads to unintended signaling consequences that interfere with data interpretation. To overcome these obstacles, researchers employed protein engineering approaches that enable tight control of protein function through allosteric mechanisms. These methods provide high specificity as well as spatial and temporal precision in regulation of protein activity in vitro and in vivo. In this review, we focus on the recent advancements in engineered allosteric regulation and discuss the various bioengineered allosteric techniques available now, from chimeric GPCRs to chemogenetic and optogenetic switches. We highlight the benefits and pitfalls of each of these techniques as well as areas in which future improvements can be made. Additionally, we provide a brief discussion on implementation of engineered allosteric regulation approaches, demonstrating that these tools can shed light on elusive biological events and have the potential to be utilized in precision medicine.

Parallel Chemoselective Profiling for Mapping Protein Structure

Solution-based structural techniques complement high-resolution structural data by providing insight into the oft-missed links between protein structure and dynamics. Here, we present Parallel Chemoselective Profiling, a solution-based structural method for characterizing protein structure and dynamics. Our method utilizes deep mutational scanning saturation mutagenesis data to install amino acid residues with specific chemistries at defined positions on the solvent-exposed surface of a protein. Differences in the extent of labeling of installed mutant residues are quantified using targeted mass spectrometry, reporting on each residue's local environment and structural dynamics. Using our method, we studied how conformation-selective, ATP-competitive inhibitors affect the local and global structure and dynamics of full-length Src kinase. Our results highlight how parallel chemoselective profiling can be used to study a dynamic multi-domain protein, and suggest that our method will be a useful addition to the relatively small toolkit of existing protein footprinting techniques.

How ATP-Competitive Inhibitors Allosterically Modulate Tyrosine Kinases That Contain a Src-like Regulatory Architecture

Small molecule kinase inhibitors that stabilize distinct ATP binding site conformations can differentially modulate the global conformation of Src-family kinases (SFKs). However, it is unclear which specific ATP binding site contacts are responsible for modulating the global conformation of SFKs and whether these inhibitor-mediated allosteric effects generalize to other tyrosine kinases. Here, we describe the development of chemical probes that allow us to deconvolute which features in the ATP binding site are responsible for the allosteric modulation of the global conformation of Src. We find that the ability of an inhibitor to modulate the global conformation of Src's regulatory domain-catalytic domain module relies mainly on the influence it has on the conformation of a structural element called helix αC. Furthermore, by developing a set of orthogonal probes that target a drug-sensitized Src variant, we show that stabilizing Src's helix αC in an active conformation is sufficient to promote a Src-mediated, phosphotransferase-independent alteration in cell morphology. Finally, we report that ATP-competitive, conformation-selective inhibitors can influence the global conformation of tyrosine kinases beyond the SFKs, suggesting that the allosteric networks we observe in Src are conserved in kinases that have a similar regulatory architecture. Our study highlights that an ATP-competitive inhibitor's interactions with helix αC can have a major influence on the global conformation of some tyrosine kinases.

Multi-state recognition pathway of the intrinsically disordered protein kinase inhibitor by protein kinase A

In the nucleus, the spatiotemporal regulation of the catalytic subunit of cAMP-dependent protein kinase A (PKA-C) is orchestrated by an intrinsically disordered protein kinase inhibitor, PKI, which recruits the CRM1/RanGTP nuclear exporting complex. How the PKA-C/PKI complex assembles and recognizes CRM1/RanGTP is not well understood. Using NMR, SAXS, fluorescence, metadynamics, and Markov model analysis, we determined the multi-state recognition pathway for PKI. After a fast binding step in which PKA-C selects PKI's most competent conformations, PKI folds upon binding through a slow conformational rearrangement within the enzyme's binding pocket. The high-affinity and pseudo-substrate regions of PKI become more structured and the transient interactions with the kinase augment the helical content of the nuclear export sequence, which is then poised to recruit the CRM1/RanGTP complex for nuclear translocation. The multistate binding mechanism featured by PKA-C/PKI complex represents a paradigm on how disordered, ancillary proteins (or protein domains) are able to operate multiple functions such as inhibiting the kinase while recruiting other regulatory proteins for nuclear export.

Molecular basis for receptor tyrosine kinase A-loop tyrosine transphosphorylation

A long-standing mystery shrouds the mechanism by which catalytically repressed receptor tyrosine kinase domains accomplish transphosphorylation of activation loop (A-loop) tyrosines. Here we show that this reaction proceeds via an asymmetric complex that is thermodynamically disadvantaged because of an electrostatic repulsion between enzyme and substrate kinases. Under physiological conditions, the energetic gain resulting from ligand-induced dimerization of extracellular domains overcomes this opposing clash, stabilizing the A-loop-transphosphorylating dimer. A unique pathogenic fibroblast growth factor receptor gain-of-function mutation promotes formation of the complex responsible for phosphorylation of A-loop tyrosines by eliminating this repulsive force. We show that asymmetric complex formation induces a more phosphorylatable A-loop conformation in the substrate kinase, which in turn promotes the active state of the enzyme kinase. This explains how quantitative differences in the stability of ligand-induced extracellular dimerization promotes formation of the intracellular A-loop-transphosphorylating asymmetric complex to varying extents, thereby modulating intracellular kinase activity and signaling intensity.

Strategies for Engineering and Rewiring Kinase Regulation

Eukaryotic protein kinases (EPKs) catalyze the transfer of a phosphate group onto another protein in response to appropriate regulatory cues. In doing so, they provide a primary means for cellular information transfer. Consequently, EPKs play crucial roles in cell differentiation and cell-cycle progression, and kinase dysregulation is associated with numerous disease phenotypes including cancer. Nonnative cues for synthetically regulating kinases are thus much sought after, both for dissecting cell signaling pathways and for pharmaceutical development. In recent years advances in protein engineering and sequence analysis have led to new approaches for manipulating kinase activity, localization, and in some instances specificity. These tools have revealed fundamental principles of intracellular signaling and suggest paths forward for the design of therapeutic allosteric kinase regulators.

Engineering allosteric regulation in protein kinases

Phosphoregulation, in which the addition of a negatively charged phosphate group modulates protein activity, enables dynamic cellular responses. To understand how new phosphoregulation might be acquired, we mutationally scanned the surface of a prototypical yeast kinase (Kss1) to identify potential regulatory sites. The data revealed a set of spatially distributed "hotspots" that might have coevolved with the active site and preferentially modulated kinase activity. By engineering simple consensus phosphorylation sites at these hotspots, we rewired cell signaling in yeast. Using the same approach with a homolog yeast mitogen-activated protein kinase, Hog1, we introduced new phosphoregulation that modified its localization and signaling dynamics. Beyond revealing potential use in synthetic biology, our findings suggest that the identified hotspots contribute to the diversity of natural allosteric regulatory mechanisms in the eukaryotic kinome and, given that some are mutated in cancers, understanding these hotspots may have clinical relevance to human disease.

Aurora-PLK1 cascades as key signaling modules in the regulation of mitosis

Mitosis is controlled by reversible protein phosphorylation involving specific kinases and phosphatases. A handful of major mitotic protein kinases, such as the cyclin B-CDK1 complex, the Aurora kinases, and Polo-like kinase 1 (PLK1), cooperatively regulate distinct mitotic processes. Research has identified proteins and mechanisms that integrate these kinases into signaling cascades that guide essential mitotic events. These findings have important implications for our understanding of the mechanisms of mitotic regulation and may advance the development of novel antimitotic drugs. We review collected evidence that in vertebrates, the Aurora kinases serve as catalytic subunits of distinct complexes formed with the four scaffold proteins Bora, CEP192, INCENP, and TPX2, which we deem "core" Aurora cofactors. These complexes and the Aurora-PLK1 cascades organized by Bora, CEP192, and INCENP control crucial aspects of mitosis and all pathways of spindle assembly. We compare the mechanisms of Aurora activation in relation to the different spindle assembly pathways and draw a functional analogy between the CEP192 complex and the chromosomal passenger complex that may reflect the coevolution of centrosomes, kinetochores, and the actomyosin cleavage apparatus. We also analyze the roles and mechanisms of Aurora-PLK1 signaling in the cell and centrosome cycles and in the DNA damage response.

Structure and evolution of the Fam20 kinases

The Fam20 proteins are novel kinases that phosphorylate secreted proteins and proteoglycans. Fam20C phosphorylates hundreds of secreted proteins and is activated by the pseudokinase Fam20A. Fam20B phosphorylates a xylose residue to regulate proteoglycan synthesis. Despite these wide-ranging and important functions, the molecular and structural basis for the regulation and substrate specificity of these kinases are unknown. Here we report molecular characterizations of all three Fam20 kinases, and show that Fam20C is activated by the formation of an evolutionarily conserved homodimer or heterodimer with Fam20A. Fam20B has a unique active site for recognizing Galβ1-4Xylβ1, the initiator disaccharide within the tetrasaccharide linker region of proteoglycans. We further show that in animals the monomeric Fam20B preceded the appearance of the dimeric Fam20C, and the dimerization trait of Fam20C emerged concomitantly with a change in substrate specificity. Our results provide comprehensive structural, biochemical, and evolutionary insights into the function of the Fam20 kinases.

Two FGF Receptor Kinase Molecules Act in Concert to Recruit and Transphosphorylate Phospholipase Cγ

The molecular basis by which receptor tyrosine kinases (RTKs) recruit and phosphorylate Src Homology 2 (SH2) domain-containing substrates has remained elusive. We used X-ray crystallography, NMR spectroscopy, and cell-based assays to demonstrate that recruitment and phosphorylation of Phospholipase Cγ (PLCγ), a prototypical SH2 containing substrate, by FGF receptors (FGFR) entails formation of an allosteric 2:1 FGFR-PLCγ complex. We show that the engagement of pTyr-binding pocket of the cSH2 domain of PLCγ by the phosphorylated tail of an FGFR kinase induces a conformational change at the region past the cSH2 core domain encompassing Tyr-771 and Tyr-783 to facilitate the binding/phosphorylation of these tyrosines by another FGFR kinase in trans. Our data overturn the current paradigm that recruitment and phosphorylation of substrates are carried out by the same RTK monomer in cis and disclose an obligatory role for receptor dimerization in substrate phosphorylation in addition to its canonical role in kinase activation.

A structural atlas of kinases inhibited by clinically approved drugs

The aberrant activation of protein kinases is associated with many human diseases, most notably cancer. Due to this link between kinase deregulation and disease progression, kinases are one of the most targeted protein families for small-molecule inhibition. Within the last 15 years, the U.S. Food and Drug Administration has approved over 20 small-molecule inhibitors of protein kinases for use in the clinic. These inhibitors target the kinase active site and represent the successful hurdling by medicinal chemists of the formidable challenge posed by the high similarity among the active sites of the approximately 500 human kinases. We review the conserved structural features of kinases that are important for inhibitor binding as well as for catalysis. Many clinically approved drugs elicit selectivity by exploiting subtle variation within the kinase active site. We highlight some of the crystallographic studies on the kinase-inhibitor complexes that have provided valuable guidance for the development of these drugs as well as for future drug design efforts.

Going for broke: targeting the human cancer pseudokinome

Protein phosphorylation lies at the heart of cell signalling, and somatic mutation(s) in kinases drives and sustains a multitude of human diseases, including cancer. The human protein kinase superfamily (the kinome) encodes approximately 50 'pseudokinases', which were initially predicted to be incapable of dynamic cell signalling when compared with canonical enzymatically active kinases. This assumption was supported by bioinformatics, which showed that amino acid changes at one or more key loci, making up the nucleotide-binding site or phosphotransferase machinery, were conserved in multiple vertebrate and non-vertebrate pseudokinase homologues. Protein kinases are highly attractive targets for drug discovery, as evidenced by the approval of almost 30 kinase inhibitors in oncology, and the successful development of the dual JAK1/2 (Janus kinase 1/2) inhibitor ruxolitinib for inflammatory indications. However, for such a large (>550) protein family, a remarkable number have still not been analysed at the molecular level, and only a surprisingly small percentage of kinases have been successfully targeted clinically. This is despite evidence that many are potential candidates for the development of new therapeutics. Indeed, several recent reports confirm that disease-associated pseudokinases can bind to nucleotide co-factors at concentrations achievable in the cell. Together, these findings suggest that drug targeting using either ATP-site or unbiased ligand-discovery approaches should now be attempted using the validation technology currently employed to evaluate their classic protein kinase counterparts. In the present review, we discuss members of the human pseudokinome repertoire, and catalogue somatic amino acid pseudokinase mutations that are emerging as the depth and clinical coverage of the human cancer pseudokinome expand.

Specificity rendering 'hot-spots' for aurora kinase inhibitor design: the role of non-covalent interactions and conformational transitions

The present study examines the conformational transitions occurring among the major structural motifs of Aurora kinase (AK) concomitant with the DFG-flip and deciphers the role of non-covalent interactions in rendering specificity. Multiple sequence alignment, docking and structural analysis of a repertoire of 56 crystal structures of AK from Protein Data Bank (PDB) has been carried out. The crystal structures were systematically categorized based on the conformational disposition of the DFG-loop [in (DI) 42, out (DO) 5 and out-up (DOU) 9], G-loop [extended (GE) 53 and folded (GF) 3] and αC-helix [in (CI) 42 and out (CO) 14]. The overlapping subsets on categorization show the inter-dependency among structural motifs. Therefore, the four distinct possibilities a) 2W1C (DI, CI, GE) b) 3E5A (DI, CI, GF) c) 3DJ6 (DI, CO, GF) d) 3UNZ (DOU, CO, GF) along with their co-crystals and apo-forms were subjected to molecular dynamics simulations of 40 ns each to evaluate the variations of individual residues and their impact on forming interactions. The non-covalent interactions formed by the 157 AK co-crystals with different regions of the binding site were initially studied with the docked complexes and structure interaction fingerprints. The frequency of the most prominent interactions was gauged in the AK inhibitors from PDB and the four representative conformations during 40 ns. Based on this study, seven major non-covalent interactions and their complementary sites in AK capable of rendering specificity have been prioritized for the design of different classes of inhibitors.

Dimerization-induced allostery in protein kinase regulation

The ability of protein kinases to switch between inactive and active states is critical to control the outputs of cellular signaling pathways. In several protein kinases, the conformation of helix αC is a key hub on which regulatory inputs converge to induce catalytic switching. An emerging mechanism involved in regulating helix αC orientation is the allosteric coupling with kinase domain surfaces involved in homo- or heterodimerization. In this review, we discuss dimerization-mediated regulation of the rapidly accelerated fibrosarcoma (RAF) and eIF2α kinase families and draw parallels with the analogous behavior of the epidermal growth factor receptor (EGFR) and serine/threonine-protein kinase endoribonuclease 1 (IRE1)/ribonuclease L (RNAse L) kinase families. Given that resistance to RAF-targeted therapeutics often stems from dimerization-dependent mechanisms, we suggest that a better understanding of dimerization-induced allostery may assist in developing alternate therapeutic strategies.

IRAK4 dimerization and trans-autophosphorylation are induced by Myddosome assembly

Trans-autophosphorylation is among the most prevalent means of protein kinase activation, yet its molecular basis is poorly defined. In Toll-like receptor and interleukin-1 receptor signaling pathways, the kinase IRAK4 is recruited to the membrane-proximal adaptor MyD88 through death domain (DD) interactions, forming the oligomeric Myddosome and mediating NF-κB activation. Here we show that unphosphorylated IRAK4 dimerizes in solution with a KD of 2.5 μM and that Myddosome assembly greatly enhances IRAK4 kinase domain (KD) autophosphorylation at sub-KD concentrations. The crystal structure of the unphosphorylated IRAK4(KD) dimer captures a conformation that appears to represent the actual trans-autophosphorylation reaction, with the activation loop phosphosite of one IRAK4 monomer precisely positioned for phosphotransfer by its partner. We show that dimerization is crucial for IRAK4 autophosphorylation in vitro and ligand-dependent signaling in cells. These studies identify a mechanism for oligomerization-driven allosteric autoactivation of IRAK4 that may be general to other kinases activated by autophosphorylation.

Divergent modulation of Src-family kinase regulatory interactions with ATP-competitive inhibitors

Multidomain protein kinases, central controllers of signal transduction, use regulatory domains to modulate catalytic activity in a complex cellular environment. Additionally, these domains regulate noncatalytic functions, including cellular localization and protein–protein interactions. Src-family kinases (SFKs) are promising therapeutic targets for a number of diseases and are an excellent model for studying the regulation of multidomain kinases. Here, we demonstrate that the regulatory domains of the SFKs Src and Hck are divergently affected by ligands that stabilize two distinct inactive ATP-binding site conformations. Conformation-selective, ATP-competitive inhibitors differentially modulate the ability of the SH3 and SH2 domains of Src and Hck to engage in intermolecular interactions and the ability of the kinase–inhibitor complex to undergo post-translational modification by effector enzymes. This surprising divergence in regulatory domain behavior by two classes of inhibitors that each stabilize inactive ATP-binding site conformations is found to occur through perturbation or stabilization of the αC helix. These studies provide insight into how conformation-selective, ATP-competitive inhibitors can be designed to modulate domain interactions and post-translational modifications distal to the ATP-binding site of kinases.

NMR mapping of protein conformational landscapes using coordinated behavior of chemical shifts upon ligand binding

Proteins exist as an ensemble of conformers that are distributed on free energy landscapes resembling folding funnels. While the most stable conformers populate low energy basins, protein function is often carried out through low-populated conformational states that occupy high energy basins. Ligand binding shifts the populations of these states, changing the distribution of these conformers. Understanding how the equilibrium among the states is altered upon ligand binding, interaction with other binding partners, and/or mutations and post-translational modifications is of critical importance for explaining allosteric signaling in proteins. Here, we propose a statistical analysis of the linear trajectories of NMR chemical shifts (CONCISE, COordiNated ChemIcal Shifts bEhavior) for the interpretation of protein conformational equilibria. CONCISE enables one to quantitatively measure the population shifts associated with ligand titrations and estimate the degree of collectiveness of the protein residues' response to ligand binding, giving a concise view of the structural transitions. The combination of CONCISE with thermocalorimetric and kinetic data allows one to depict a protein's approximate conformational energy landscape. We tested this method with the catalytic subunit of cAMP-dependent protein kinase A, a ubiquitous enzyme that undergoes conformational transitions upon both nucleotide and pseudo-substrate binding. When complemented with chemical shift covariance analysis (CHESCA), this new method offers both collective response and residue-specific correlations for ligand binding to proteins.

Cracking the molecular origin of intrinsic tyrosine kinase activity through analysis of pathogenic gain-of-function mutations

The basal (ligand-independent) kinase activity of receptor tyrosine kinases (RTKs) promotes trans-phosphorylation on activation loop tyrosines upon ligand-induced receptor dimerization, thus upregulating intrinsic kinase activity and triggering intracellular signaling. To understand the molecular determinants of intrinsic kinase activity, we used X-ray crystallography and NMR spectroscopy to analyze pathogenic FGF receptor mutants with gradations in gain-of-function activity. These structural analyses revealed a "two-state" dynamic equilibrium model whereby the kinase toggles between an "inhibited," structurally rigid ground state and a more dynamic and heterogeneous active state. The pathogenic mutations have different abilities to shift this equilibrium toward the active state. The increase in the fractional population of FGF receptors in the active state correlates with the degree of gain-of-function activity and clinical severity. Our data demonstrate that the fractional population of RTKs in the active state determines intrinsic kinase activity and underscore how a slight increase in the active population of kinases can have grave consequences for human health.

Structural and functional diversity in the activity and regulation of DAPK-related protein kinases

Within the large group of calcium/calmodulin-dependent protein kinases (CAMKs) of the human kinome, there is a distinct branch of highly related kinases that includes three families: death-associated protein-related kinases, myosin light-chain-related kinases and triple functional domain protein-related kinases. In this review, we refer to these collectively as DMT kinases. There are several functional features that span the three families, such as a broad involvement in apoptotic processes, cytoskeletal association and cellular plasticity. Other CAMKs contain a highly conserved HRD motif, which is a prerequisite for kinase regulation through activation-loop phosphorylation, but in all 16 members of the DMT branch, this is replaced by an HF/LD motif. This DMT kinase signature motif substitutes phosphorylation-dependent active-site interactions with a local hydrophobic core that maintains an active kinase conformation. Only about half of the DMT kinases have an additional autoregulatory domain, C-terminal to the kinase domain that binds calcium/calmodulin in order to regulate kinase activity. Protein substrates have been identified for some of the DMT kinases, but little is known about the mechanism of recognition. Substrate conformation could be an equally important parameter in substrate recognition as specific preferences in sequence position. Taking the data together, this kinase branch encapsulates a treasure trove of features that renders it distinct from many other protein kinases and calls for future research activities in this field.

Regulation of JAK2 activation by Janus homology 2: evidence from molecular dynamics simulations

Janus kinase 2 (JAK2) is a protein tyrosine kinase implicated in signaling by specific members of the cytokine receptor family. Although it has been established that the JAK2 tyrosine kinase is negatively regulated by the JAK homology 2 (JH2) pseudokinase domain, the underlying mechanism of JH2 mediated regulation remains elusive. To elucidate the regulation of JAK2 kinase, we have built a structural model for the kinase and pseudokinase domains of JAK2. An asymmetric dimer is proposed, in which the kinase domain JH1 occupies a position where it could not be activated. We investigate the dynamic and energetic properties of the dimer by molecular dynamics simulation. JAK2 activation requires the two domains to be dissociated and rearranged in a form such that the JH1 kinase domain can adopt an active conformation. The significance of the above mechanism is emphasized by the finding that the activating V617F mutation destabilizes JH1-JH2 association in the proposed asymmetric dimer. Thus abrogation of the domain-domain interaction seems to be a possible first step for the structural rearrangement of the two domains, resulting in constitutive activation of JAK2 by the V617F mutation.

Response of rotation-translation blocked proteins using Langevin dynamics on a locally harmonic landscape

Langevin dynamics is used to compute the time evolution of the nonequilibrium motion of the atomic coordinates of a protein in response to ligand dissociation. The protein potential energy surface (PES) is approximated by a harmonic basin about the minimum of the unliganded state. Upon ligand dissociation, the protein undergoes relaxation from the bound to the unbound state. A coarse graining scheme based on rotation translation blocks (RTB) is applied to the relaxation of the two domain iron transport protein, ferric binding protein. This scheme provides a natural and efficient way to freeze out the small amplitude, high frequency motions within each rigid fragment, thereby allowing for the number of dynamical degrees of freedom to be reduced. The results obtained from all flexible atom (constraint free) dynamics are compared to those obtained using RTB-Langevin dynamics. To assess the impact of the assumed rigid fragment clustering on the temporal relaxation dynamics of the protein molecule, three distinct rigid block decompositions were generated and their responses compared. Each of the decompositions was a variant of the one-block-per-residue grouping, with their force and friction matrices being derived from their fully flexible counterpart. Monitoring the time evolution of the distance separating a selected pair of amino acids, the response curves of the blocked decompositions were similar in shape to each other and to the control system in which all atomic degrees of freedom are fully independent. The similar shape of the blocked responses showed that the variations in grouping had only a minor impact on the kinematics. Compared with the all atom responses, however, the blocked responses were faster as a result of the instantaneous transmission of force throughout each rigid block. This occurred because rigid blocking does not permit any intrablock deformation that could store or divert energy. It was found, however, that this accelerated response could be successfully corrected by scaling each eigenvalue in the appropriate propagation matrix by the least-squares fitted slope of the blocked vs nonblocked eigenvalue spectra. The RTB responses for each test system were dominated by small eigenvalue overdamped Langevin modes. The large eigenvalue members of each response dissipated within the first 5 ps, after which the long time response was dominated by a modest set of low energy, overdamped normal modes, that were characterized by highly cooperative, functionally relevant displacements. The response assuming that the system is in the overdamped limit was compared to the full phase space Langevin dynamics results. The responses after the first 5 ps were nearly identical, confirming that the inertial components were significant only in the initial stages of the relaxation. Since the propagator matrix in the overdamped formulation is real-symmetric and does not require the inertial component in the propagator, the computation time and memory footprint was reduced by 1 order of magnitude.

Effect of multimeric structure of CaMKII in the GluN2B-mediated modulation of kinetic parameters of ATP

Interaction of GluN2B subunit of N-methyl-D-aspartate receptor with calcium/calmodulin dependent protein kinase II (CaMKII) is critical for the induction of long term potentiation at hippocampal CA3-CA1 synapses. We have previously reported that CaMKII binding to GluN2B increases its affinity but abolishes the cooperativity for ATP. In the present study, we demonstrate that the reduction in S_0.5 for ATP of an individual CaMKII subunit seems to be directly induced by the binding of GluN2B to the same subunit, while any GluN2B induced effects on the cooperativity and maximal velocity would additionally require the CaMKII holoenzyme structure. We measured the apparent kinetic parameters for ATP using an association domain truncated monomeric CaMKII and a heteromultimeric CaMKII (having subunits that are either GluN2B binding defective or ATP binding defective), in the presence of GluN2A or GluN2B substrates. The S_0.5 value for ATP of monomeric CaMKII is reduced ∼ 3 fold by the presence of GluN2B suggesting that the induced change in affinity for ATP is independent of the holoenzyme structure. The heteromultimeric mutant of CaMKII, did not exhibit cooperativity of ATP binding probably because of the interspersing of ATP binding defective subunits in the holoenzyme. In contrast to the wild type holoenzyme, presence of GluN2B increased the V_max of monomeric CaMKII which resulted in an approximately 4.0 fold increase in the apparent catalytic constant (V_max/S_0.5) as compared to GluN2A. The kinetic parameter values of the heteromultimeric CaMKII for ATP, on the other hand, did not show any significant difference between the phosphorylation of GluN2B and GluN2A suggesting that modulation requires binding of GluN2B to the same subunit. Overall, our present study provides insights into the role of multimeric structure of CaMKII in GluN2B-mediated regulation.

Regulation of protein kinases by lipids

Membranes are sites of intense signaling activity within the cell, serving as dynamic scaffolds for the recruitment of signaling molecules and their substrates. The specific and reversible localization of these signaling molecules to membranes is critical for the appropriate activation of downstream signaling pathways. Phospholipid-binding domains, including C1, C2, PH, and PX domains, play critical roles in the membrane targeting of protein kinases. Recent structural studies have identified a new membrane association domain, the Kinase Associated 1 (KA1) domain, which targets a number of yeast and mammalian protein kinases to membranes containing acidic phospholipids. Despite an abundance of localization studies on lipid-binding proteins and structural studies of the isolated lipid-binding domains, the question of how membrane binding is coupled to the activation of the kinase catalytic domain has been virtually untouched. Recently, structural studies on protein kinase C (PKC) have provided some of the first structural insights into the allosteric regulation of protein kinases by lipid second messengers.

Computational modeling of allosteric communication reveals organizing principles of mutation-induced signaling in ABL and EGFR kinases

The emerging structural information about allosteric kinase complexes and the growing number of allosteric inhibitors call for a systematic strategy to delineate and classify mechanisms of allosteric regulation and long-range communication that control kinase activity. In this work, we have investigated mechanistic aspects of long-range communications in ABL and EGFR kinases based on the results of multiscale simulations of regulatory complexes and computational modeling of signal propagation in proteins. These approaches have been systematically employed to elucidate organizing molecular principles of allosteric signaling in the ABL and EGFR multi-domain regulatory complexes and analyze allosteric signatures of the gate-keeper cancer mutations. We have presented evidence that mechanisms of allosteric activation may have universally evolved in the ABL and EGFR regulatory complexes as a product of a functional cross-talk between the organizing αF-helix and conformationally adaptive αI-helix and αC-helix. These structural elements form a dynamic network of efficiently communicated clusters that may control the long-range interdomain coupling and allosteric activation. The results of this study have unveiled a unifying effect of the gate-keeper cancer mutations as catalysts of kinase activation, leading to the enhanced long-range communication among allosterically coupled segments and stabilization of the active kinase form. The results of this study can reconcile recent experimental studies of allosteric inhibition and long-range cooperativity between binding sites in protein kinases. The presented study offers a novel molecular insight into mechanistic aspects of allosteric kinase signaling and provides a quantitative picture of activation mechanisms in protein kinases at the atomic level.

Despite recent progress in computational and experimental studies of dynamic regulation in protein kinases, a mechanistic understanding of long-range communication and mechanisms of mutation-induced signaling controlling kinase activity remains largely qualitative. In this study, we have performed a systematic modeling and analysis of allosteric activation in ABL and EGFR kinases at the increasing level of complexity - from catalytic domain to multi-domain regulatory complexes. The results of this study have revealed organizing structural and mechanistic principles of allosteric signaling in protein kinases. Although activation mechanisms in ABL and EGFR kinases have evolved through acquisition of structurally different regulatory complexes, we have found that long-range interdomain communication between common functional segments (αF-helix and αC-helix) may be important for allosteric activation. The results of study have revealed molecular signatures of activating cancer mutations and have shed the light on general mechanistic aspects of mutation-induced signaling in protein kinases. An advanced understanding and further characterization of molecular signatures of kinase mutations may aid in a better rationalization of mutational effects on clinical outcomes and facilitate molecular-based therapeutic strategies to combat kinase mutation-dependent tumorigenesis.

The energy landscape analysis of cancer mutations in protein kinases

The growing interest in quantifying the molecular basis of protein kinase activation and allosteric regulation by cancer mutations has fueled computational studies of allosteric signaling in protein kinases. In the present study, we combined computer simulations and the energy landscape analysis of protein kinases to characterize the interplay between oncogenic mutations and locally frustrated sites as important catalysts of allostetric kinase activation. While structurally rigid kinase core constitutes a minimally frustrated hub of the catalytic domain, locally frustrated residue clusters, whose interaction networks are not energetically optimized, are prone to dynamic modulation and could enable allosteric conformational transitions. The results of this study have shown that the energy landscape effect of oncogenic mutations may be allosteric eliciting global changes in the spatial distribution of highly frustrated residues. We have found that mutation-induced allosteric signaling may involve a dynamic coupling between structurally rigid (minimally frustrated) and plastic (locally frustrated) clusters of residues. The presented study has demonstrated that activation cancer mutations may affect the thermodynamic equilibrium between kinase states by allosterically altering the distribution of locally frustrated sites and increasing the local frustration in the inactive form, while eliminating locally frustrated sites and restoring structural rigidity of the active form. The energy landsape analysis of protein kinases and the proposed role of locally frustrated sites in activation mechanisms may have useful implications for bioinformatics-based screening and detection of functional sites critical for allosteric regulation in complex biomolecular systems.

A conserved non-canonical motif in the pseudoactive site of the ROP5 pseudokinase domain mediates its effect on Toxoplasma virulence

The ROP5 family is a closely related set of polymorphic pseudokinases that are critical to the ability of Toxoplasma to cause disease. Polymorphisms in ROP5 also make it a major determinant of strain-specific differences in virulence. ROP5 possesses all of the major kinase motifs required for catalysis except for a substitution at the catalytic Asp. We show that this substitution in the catalytic loop of ROP5 is part of a motif conserved in other pseudokinases of both Toxoplasma and human origin, and that this motif is required for the full activity in vivo of ROP5. This suggests evolutionary selection at this site for a biochemical function, rather than simple drift away from catalysis. We present the crystal structures of a virulent isoform of ROP5 both in its ATP-bound and -unbound states and have demonstrated that despite maintaining the canonical ATP-binding motifs, ROP5 binds ATP in a distorted conformation mediated by unusual magnesium coordination sites that would not be predicted from the primary sequence. In addition, we have mapped the polymorphisms spread throughout the primary sequence of ROP5 to two major surfaces, including the activation segment of ROP5. This suggests that the pseudoactive site of this class of pseudokinases may have evolved to use the canonical ATP-binding motifs for non-catalytic signaling through allostery.

Structural insights into the architecture and allostery of full-length AMP-activated protein kinase

AMP-activated protein kinase (AMPK) is a heterotrimeric complex composed of α catalytic subunit, β scaffolding subunit, and γ regulatory subunit with critical roles in maintaining cellular energy homeostasis. However, the molecular architecture of the intact complex and the allostery associated with the adenosine binding-induced regulation of kinase activity remain unclear. Here, we determine the three-dimensional reconstruction and subunit organization of the full-length rat AMPK (α1β1γ1) through single-particle electron-microscopy. By comparing the structures of AMPK in ATP- and AMP-bound states, we are able to visualize the sequential conformational changes underlying kinase activation that transmits from the adenosine binding sites in the γ subunit to the kinase domain of the α subunit. These results not only make substantial revision to the current model of AMPK assembly, but also highlight a central role of the linker sequence of the α subunit in mediating the allostery of AMPK.

Protein kinase biochemistry and drug discovery

Protein kinases are fascinating biological catalysts with a rapidly expanding knowledge base, a growing appreciation in cell regulatory control, and an ascendant role in successful therapeutic intervention. To better understand protein kinases, the molecular underpinnings of phosphoryl group transfer, protein phosphorylation, and inhibitor interactions are examined. This analysis begins with a survey of phosphate group and phosphoprotein properties which provide context to the evolutionary selection of phosphorylation as a central mechanism for biological regulation of most cellular processes. Next, the kinetic and catalytic mechanisms of protein kinases are examined with respect to model aqueous systems to define the elements of catalysis. A brief structural biology overview further delves into the molecular basis of catalysis and regulation of catalytic activity. Concomitant with a prominent role in normal physiology, protein kinases have important roles in the disease state. To facilitate effective kinase drug discovery, classic and emerging approaches for characterizing kinase inhibitors are evaluated including biochemical assay design, inhibitor mechanism of action analysis, and proper kinetic treatment of irreversible inhibitors. As the resulting protein kinase inhibitors can modulate intended and unintended targets, profiling methods are discussed which can illuminate a more complete range of an inhibitor's biological activities to enable more meaningful cellular studies and more effective clinical studies. Taken as a whole, a wealth of protein kinase biochemistry knowledge is available, yet it is clear that a substantial extent of our understanding in this field remains to be discovered which should yield many new opportunities for therapeutic intervention.

An Atg13 protein-mediated self-association of the Atg1 protein kinase is important for the induction of autophagy

Autophagy pathways in eukaryotic cells mediate the turnover of a diverse set of cytoplasmic components, including damaged organelles and abnormal protein aggregates. Autophagy-mediated degradation is highly regulated, and defects in these pathways have been linked to a number of human disorders. The Atg1 protein kinase appears to be a key site of this control and is targeted by multiple signaling pathways to ensure the appropriate autophagic response to changing environmental conditions. Despite the importance of this kinase, relatively little is known about the molecular details of Atg1 activation. In this study we show that Atg13, an evolutionarily conserved regulator of Atg1, promotes the formation of a specific Atg1 self-interaction in the budding yeast, Saccharomyces cerevisiae. The appearance of this Atg1-Atg1 complex is correlated with the induction of autophagy, and conditions that disrupt this complex result in diminished levels of both autophagy and Atg1 kinase activity. Moreover, the addition of a heterologous dimerization domain to Atg1 resulted in elevated kinase activity both in vivo and in vitro. The formation of this complex appears to be an important prerequisite for the subsequent autophosphorylation of Thr-226 in the Atg1 activation loop. Previous work indicates that this modification is necessary and perhaps sufficient for Atg1 kinase activity. Interestingly, this Atg1 self-association does not require Atg17, suggesting that this second conserved regulator might activate Atg1 in a manner mechanistically distinct from that of Atg13. In all, this work suggests a model whereby this self-association stimulates the autophosphorylation of Atg1 within its activation loop.

Centrosomal protein of 192 kDa (Cep192) promotes centrosome-driven spindle assembly by engaging in organelle-specific Aurora A activation

Centrosomes are primary microtubule (MT)-organizing centers (MTOCs). During mitosis, they dramatically increase their size and MT-nucleating activity and participate in spindle assembly from spindle poles. These events require the serine/threonine kinase, Aurora A (AurA), and the centrosomal protein of 192 kDa (Cep192)/spindle defective 2 (Spd-2), but the underlying mechanism remains unclear. We have found that Cep192, unlike targeting protein for Xklp2 (TPX2), a known MT-localizing AurA activator, is an AurA cofactor in centrosome-driven spindle assembly. Cep192, through a direct interaction, targets AurA to mitotic centrosomes where the locally accumulating AurA forms homodimers or oligomers. The dimerization of endogenous AurA, in the presence of bound Cep192, triggers potent kinase activation that, in turn, drives MT assembly. Depletion of Cep192 or specific interference with AurA-Cep192 binding did not prevent AurA oligomerization on MTs but abrogated AurA recruitment to centrosomes and its activation by either sperm nuclei or anti-AurA antibody (αAurA)-induced dimerization. In these settings, MT assembly by both centrosomes and αAurA-coated beads was also abolished or severely compromised. Hence, Cep192 activates AurA by a mechanism different from that previously described for TPX2. The Cep192-mediated mechanism maximizes AurA activity at centrosomes and appears essential for the function of these organelles as MTOCs.

Exploring and exploiting allostery: Models, evolution, and drug targeting

The concept of allostery was elaborated almost 50years ago by Monod and coworkers to provide a framework for interpreting experimental studies on the regulation of protein function. In essence, binding of a ligand at an allosteric site affects the function at a distant site exploiting protein flexibility and reshaping protein energy landscape. Both monomeric and oligomeric proteins can be allosteric. In the past decades, the behavior of allosteric systems has been analyzed in many investigations while general theoretical models and variations thereof have been steadily proposed to interpret the experimental data. Allostery has been established as a fundamental mechanism of regulation in all organisms, governing a variety of processes that range from metabolic control to receptor function and from ligand transport to cell motility. A number of studies have shed light on how evolutionary pressures have favored and molded the development of allosteric features in specific macromolecular systems. The widespread occurrence of allostery has been recently exploited for the development and design of allosteric drugs that bind to either physiological or non-physiological allosteric sites leading to gain of function or loss of function. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.

Affinity reagents that target a specific inactive form of protein kinases

A number of small-molecule inhibitors have been developed that target the catalytic domains of protein kinases that are not in an active conformation. An inactive form that has been observed in several kinases is the DFG-out conformation. This conformation is characterized by an almost 180 degrees rotation of the conserved Asp-Phe-Gly (DFG) motif in the ATP-binding cleft relative to the active form. However, the sequence and structural determinants that allow a kinase to stably adopt the DFG-out conformation are not known. Here, we characterize a series of inhibitors based on a general pharmacophore for this inactive form. We demonstrate that modified versions of these inhibitors can be used to study the thermodynamics and kinetics of ligand binding to DFG-out-adopting kinases and for enriching these kinases from complex protein mixtures.

Unraveling the allosteric mechanism of serine protease inhibition by an antibody

Recent structural studies have outlined the mechanism of protease inhibition by active site-directed antibodies. However, the molecular basis of allosteric inhibition by antibodies has been elusive. Here we report the 2.35 A resolution structure of the trypsin-like serine protease hepatocyte growth factor activator (HGFA) in complex with the allosteric antibody Ab40, a potent inhibitor of HGFA catalytic activity. The antibody binds at the periphery of the substrate binding cleft and imposes a conformational change on the entire 99-loop (chymotrypsinogen numbering). The altered conformation of the 99-loop is incompatible with substrate binding due to the partial collapse of subsite S2 and the reorganization of subsite S4. Remarkably, a single residue deletion of Ab40 abolished inhibition of HGFA activity, commensurate with the reversal of the 99-loop conformation to its "competent" state. The results define an "allosteric switch" mechanism as the basis of protease inhibition by an allosteric antibody.

Molecular basis of the death-associated protein kinase-calcium/calmodulin regulator complex

Death-associated protein kinase (DAPK) provides a model for calcium-bound calmodulin (CaM)-dependent protein kinases (CaMKs). Here, we report the crystal structure of the binary DAPK-CaM complex, using a construct that includes the DAPK catalytic domain and adjacent autoregulatory domain. When DAPK was in a complex with CaM, the DAPK autoregulatory domain formed a long seven-turn helix. This DAPK-CaM module interacted with the DAPK catalytic domain through two separate domain-domain interfaces, which involved the upper and the lower lobe of the catalytic domain. When bound to DAPK, CaM adopted an extended conformation, which was different from that in CaM-CaMK peptide complexes. Complementary biochemical analysis showed that the ability of DAPK to bind CaM correlated with its catalytic activity. Because many features of CaM binding are conserved in other CaMKs, our findings likely provide a generally applicable model for regulation of CaMK activity.