Hydrophobic Core Variations Provide a Structural Framework for Tyrosine Kinase Evolution and Functional Specialization

CDK2 and CDK4: Cell Cycle Functions Evolve Distinct, Catalysis-Competent Conformations, Offering Drug Targets

Cyclin-dependent kinases (CDKs), particularly CDK4 and CDK2, are crucial for cell cycle progression from the Gap 1 (G1) to the Synthesis (S) phase by phosphorylating targets such as the Retinoblastoma Protein (Rb). CDK4, paired with cyclin-D, operates in the long G1 phase, while CDK2 with cyclin-E, manages the brief G1-to-S transition, enabling DNA replication. Aberrant CDK signaling leads to uncontrolled cell proliferation, which is a hallmark of cancer. Exactly how they accomplish their catalytic phosphorylation actions with distinct efficiencies poses the fundamental, albeit overlooked question. Here we combined available experimental data and modeling of the active complexes to establish their conformational functional landscapes to explain how the two cyclin/CDK complexes differentially populate their catalytically competent states for cell cycle progression. Our premise is that CDK catalytic efficiencies could be more important for cell cycle progression than the cyclin-CDK biochemical binding specificity and that efficiency is likely the prime determinant of cell cycle progression. We observe that CDK4 is more dynamic than CDK2 in the ATP binding site, the regulatory spine, and the interaction with its cyclin partner. The N-terminus of cyclin-D acts as an allosteric regulator of the activation loop and the ATP-binding site in CDK4. Integrated with a suite of experimental data, we suggest that the CDK4 complex is less capable of remaining in the active catalytically competent conformation, and may have a lower catalytic efficiency than CDK2, befitting their cell cycle time scales, and point to critical residues and motifs that drive their differences. Our mechanistic landscape may apply broadly to kinases, and we propose two drug design strategies: (i) allosteric Inhibition by conformational stabilization for targeting allosteric CDK4 regulation by cyclin-D, and (ii) dynamic entropy-optimized targeting which leverages the dynamic, entropic aspects of CDK4 to optimize drug binding efficacy.

Evolutionary sequence and structural basis for the distinct conformational landscapes of Tyr and Ser/Thr kinases

Protein kinases are molecular machines with rich sequence variation that distinguishes the two main evolutionary branches - tyrosine kinases (TKs) from serine/threonine kinases (STKs). Using a sequence co-variation Potts statistical energy model we previously concluded that TK catalytic domains are more likely than STKs to adopt an inactive conformation with the activation loop in an autoinhibitory "folded" conformation, due to intrinsic sequence effects. Here we investigated the structural basis for this phenomenon by integrating the sequence-based model with structure-based molecular dynamics (MD) to determine the effects of mutations on the free energy difference between active and inactive conformations, using a novel thermodynamic cycle involving many (n=108) protein-mutation free energy perturbation (FEP) simulations in the active and inactive conformations. The sequence and structure-based results are consistent and support the hypothesis that the inactive conformation "DFG-out Activation Loop Folded", is a functional regulatory state that has been stabilized in TKs relative to STKs over the course of their evolution via the accumulation of residue substitutions in the activation loop and catalytic loop that facilitate distinct substrate binding modes in trans and additional modes of regulation in cis for TKs.

Redox Regulation of Brain Selective Kinases BRSK1/2: Implications for Dynamic Control of the Eukaryotic AMPK family through Cys-based mechanisms

In eukaryotes, protein kinase signaling is regulated by a diverse array of post-translational modifications (PTMs), including phosphorylation of Ser/Thr residues and oxidation of cysteine (Cys) residues. While regulation by activation segment phosphorylation of Ser/Thr residues is well understood, relatively little is known about how oxidation of cysteine residues modulate catalysis. In this study, we investigate redox regulation of the AMPK-related Brain-selective kinases (BRSK) 1 and 2, and detail how broad catalytic activity is directly regulated through reversible oxidation and reduction of evolutionarily conserved Cys residues within the catalytic domain. We show that redox-dependent control of BRSKs is a dynamic and multilayered process involving oxidative modifications of several Cys residues, including the formation of intramolecular disulfide bonds involving a pair of Cys residues near the catalytic HRD motif and a highly conserved T-Loop Cys with a BRSK-specific Cys within an unusual CPE motif at the end of the activation segment. Consistently, mutation of the CPE-Cys increases catalytic activity in vitro and drives phosphorylation of the BRSK substrate Tau in cells. Molecular modeling and molecular dynamics simulations indicate that oxidation of the CPE-Cys destabilizes a conserved salt bridge network critical for allosteric activation. The occurrence of spatially proximal Cys amino acids in diverse Ser/Thr protein kinase families suggests that disulfide mediated control of catalytic activity may be a prevalent mechanism for regulation within the broader AMPK family.

Potts Hamiltonian Models and Molecular Dynamics Free Energy Simulations for Predicting the Impact of Mutations on Protein Kinase Stability

Single-point mutations in kinase proteins can affect their stability and fitness, and computational analysis of these effects can provide insights into the relationships among protein sequence, structure, and function for this enzyme family. To assess the impact of mutations on protein stability, we used a sequence-based Potts Hamiltonian model trained on a kinase family multiple-sequence alignment (MSA) to calculate the statistical energy (fitness) effects of mutations and compared these against relative folding free energies (ΔΔGs) calculated from all-atom molecular dynamics free energy perturbation (FEP) simulations in explicit solvent. The fitness effects of mutations in the Potts model (ΔEs) showed good agreement with experimental thermostability data (Pearson r = 0.68), similar to the correlation we observed with ΔΔGs predicted from structure-based relative FEP simulations. Recognizing the possible advantages of using Potts models to rapidly estimate protein stability effects of kinase mutations seen in cancer genomics data, we used the Potts statistical energy model to estimate the stability effects of 65 conservative and nonconservative mutations across three distinct kinases (Wee1, Abl1, and Cdc7) with somatic mutations reported in the Genomic Data Commons (GDC) database. The ΔEs of these mutations calculated from the Potts model are consistent with the corresponding ΔΔGs from FEP simulations (Pearson ratio of 0.72). The agreement between these methods suggests that the Potts model may be used as a sequence-based tool for high-throughput screening of mutational effects as part of a computational pipeline for predicting the stability effects of mutations. We also demonstrate how the scalability of the fitness-based Potts model calculations permits analyses that are not easily accessed using FEP simulations. To this end, we employed site-saturation mutagenesis in the Potts model in order to investigate the relative stability effects of mutations seen in different cancer evolutionary scenarios. We used this approach to analyze the effects of drug pressure in Abl kinase by contrasting the relative fitness penalties of somatic mutations seen in miscellaneous cancer types with those calculated for mutations associated with cancer drug resistance. We observed that, in contrast to somatic mutations of Abl seen in various tumors that appear to have evolved neutrally, cancer mutations that evolved under drug pressure in Abl-targeted therapies tend to preserve enzyme stability.

Dark kinase annotation, mining, and visualization using the Protein Kinase Ontology

The Protein Kinase Ontology (ProKinO) is an integrated knowledge graph that conceptualizes the complex relationships among protein kinase sequence, structure, function, and disease in a human and machine-readable format. In this study, we have significantly expanded ProKinO by incorporating additional data on expression patterns and drug interactions. Furthermore, we have developed a completely new browser from the ground up to render the knowledge graph visible and interactive on the web. We have enriched ProKinO with new classes and relationships that capture information on kinase ligand binding sites, expression patterns, and functional features. These additions extend ProKinO's capabilities as a discovery tool, enabling it to uncover novel insights about understudied members of the protein kinase family. We next demonstrate the application of ProKinO. Specifically, through graph mining and aggregate SPARQL queries, we identify the p21-activated protein kinase 5 (PAK5) as one of the most frequently mutated dark kinases in human cancers with abnormal expression in multiple cancers, including a previously unappreciated role in acute myeloid leukemia. We have identified recurrent oncogenic mutations in the PAK5 activation loop predicted to alter substrate binding and phosphorylation. Additionally, we have identified common ligand/drug binding residues in PAK family kinases, underscoring ProKinO's potential application in drug discovery. The updated ontology browser and the addition of a web component, ProtVista, which enables interactive mining of kinase sequence annotations in 3D structures and Alphafold models, provide a valuable resource for the signaling community. The updated ProKinO database is accessible at https://prokino.uga.edu.

Custom Workflow for the Confident Identification of Sulfotyrosine-Containing Peptides and Their Discrimination from Phosphopeptides

Protein tyrosine sulfation (sY) is a post-translational modification (PTM) catalyzed by Golgi-resident tyrosyl protein sulfo transferases (TPSTs). Information on sY in humans is currently limited to ∼50 proteins, with only a handful having verified sites of sulfation. As such, the contribution of sulfation to the regulation of biological processes remains poorly defined. Mass spectrometry (MS)-based proteomics is the method of choice for PTM analysis but has yet to be applied for systematic investigation of the "sulfome", primarily due to issues associated with discrimination of sY-containing from phosphotyrosine (pY)-containing peptides. In this study, we developed an MS-based workflow for sY-peptide characterization, incorporating optimized Zr4+ immobilized metal-ion affinity chromatography (IMAC) and TiO2 enrichment strategies. Extensive characterization of a panel of sY- and pY-peptides using an array of fragmentation regimes (CID, HCD, EThcD, ETciD, UVPD) highlighted differences in the generation of site-determining product ions and allowed us to develop a strategy for differentiating sulfated peptides from nominally isobaric phosphopeptides based on low collision energy-induced neutral loss. Application of our "sulfomics" workflow to a HEK-293 cell extracellular secretome facilitated identification of 21 new sulfotyrosine-containing proteins, several of which we validate enzymatically, and reveals new interplay between enzymes relevant to both protein and glycan sulfation.

Dynamic conformational states of apo, ATP and cabozantinib bound TAM kinases to differentiate active-inactive kinetic models

The dynamically active and inactive conformations of kinases play a crucial role in the activation of intracellular downstream signaling pathways. The all-atom molecular dynamics (MD) simulations at microsecond (µs) timescale and longer provide robust insights into the structural details of conformational alterations in kinases that contribute to their cellular metabolic activities and signaling pathways. Tyro3, Axl and Mer (TAM) receptor tyrosine kinases (RTKs) are overexpressed in several types of human cancers. Cabozantinib, a small molecule inhibitor constrains the activity of TAM kinases at nanomolar concentrations. The apo, complexes of ATP (active state) and cabozantinib (active and inactive states) with TAM RTKs were studied by 1 µs MD simulations followed by trajectory analyses. The dynamic mechanistic pathways intrinsic to the kinase activity and protein conformational landscape in the cabozantinib bound TAM kinases are revealed due to the alterations in the P-loop, α-helix and activation loop that result in breaking the regulatory (R) and catalytic (C) spines, while the active states of ATP bound TAM kinases are retained. The co-existence of dynamical states when bound to cabozantinib was observed and the long-lived kinetic transition states of distinct active and inactive structural models were deciphered from MD simulation trajectories that have not been revealed so far.Communicated by Ramaswamy H. Sarma.

Evolution of functional diversity in the holozoan tyrosine kinome

The emergence of multicellularity is strongly correlated with the expansion of tyrosine kinases, a conserved family of signaling enzymes that regulates pathways essential for cell-to-cell communication. Although tyrosine kinases have been classified from several model organisms, a molecular-level understanding of tyrosine kinase evolution across all holozoans is currently lacking. Using a hierarchical sequence constraint-based classification of diverse holozoan tyrosine kinases, we construct a new phylogenetic tree that identifies two ancient clades of cytoplasmic and receptor tyrosine kinases separated by the presence of an extended insert segment in the kinase domain connecting the D and E-helices. Present in nearly all receptor tyrosine kinases, this fast-evolving insertion imparts diverse functionalities, such as post-translational modification sites and regulatory interactions. Eph and EGFR receptor tyrosine kinases are two exceptions which lack this insert, each forming an independent lineage characterized by unique functional features. We also identify common constraints shared across multiple tyrosine kinase families which warrant the designation of three new subgroups: Src module (SrcM), insulin receptor kinase-like (IRKL), and fibroblast, platelet-derived, vascular, and growth factor receptors (FPVR). Subgroup-specific constraints reflect shared autoinhibitory interactions involved in kinase conformational regulation. Conservation analyses describe how diverse tyrosine kinase signaling functions arose through the addition of family-specific motifs upon subgroup-specific features and coevolving protein domains. We propose the oldest tyrosine kinases, IRKL, SrcM, and Csk, originated from unicellular premetazoans and were coopted for complex multicellular functions. The increased frequency of oncogenic variants in more recent tyrosine kinases suggests that lineage-specific functionalities are selectively altered in human cancers.

Cataloguing the dead: breathing new life into pseudokinase research

Pseudoenzymes are present within many, but not all, known enzyme families and lack one or more conserved canonical amino acids that help define their catalytically active counterparts. Recent findings in the pseudokinase field confirm that evolutionary repurposing of the structurally defined bilobal protein kinase fold permits distinct biological functions to emerge, many of which rely on conformational switching, as opposed to canonical catalysis. In this analysis, we evaluate progress in evaluating several members of the 'dark' pseudokinome that are pertinent to help drive this expanding field. Initially, we discuss how adaptions in erythropoietin-producing hepatocellular carcinoma (Eph) receptor tyrosine kinase domains resulted in two vertebrate pseudokinases, EphA10 and EphB6, in which co-evolving sequences generate new motifs that are likely to be important for both nucleotide binding and catalysis-independent signalling. Secondly, we discuss how conformationally flexible Tribbles pseudokinases, which have radiated in the complex vertebrates, control fundamental aspects of cell signalling that may be targetable with covalent small molecules. Finally, we show how species-level adaptions in the duplicated canonical kinase protein serine kinase histone (PSKH)1 sequence have led to the appearance of the pseudokinase PSKH2, whose physiological role remains mysterious. In conclusion, we show how the patterns we discover are selectively conserved within specific pseudokinases, and that when they are modelled alongside closely related canonical kinases, many are found to be located in functionally important regions of the conserved kinase fold. Interrogation of these patterns will be useful for future evaluation of these, and other, members of the unstudied human kinome.

Aurora A regulation by reversible cysteine oxidation reveals evolutionarily conserved redox control of Ser/Thr protein kinase activity

Reactive oxygen species (ROS) are physiological mediators of cellular signaling and play potentially damaging roles in human diseases. In this study, we found that the catalytic activity of the Ser/Thr kinase Aurora A was inhibited by the oxidation of a conserved cysteine residue (Cys290) that lies adjacent to Thr288, a critical phosphorylation site in the activation segment. Cys is present at the equivalent position in ~100 human Ser/Thr kinases, a residue that we found was important not only for the activity of human Aurora A but also for that of fission yeast MAPK-activated kinase (Srk1) and PKA (Pka1). Moreover, the presence of this conserved Cys predicted biochemical redox sensitivity among a cohort of human CAMK, AGC, and AGC-like kinases. Thus, we predict that redox modulation of the conserved Cys290 of Aurora A may be an underappreciated regulatory mechanism that is widespread in eukaryotic Ser/Thr kinases. Given the key biological roles of these enzymes, these findings have implications for understanding physiological and pathological responses to ROS and highlight the importance of protein kinase regulation through multivalent modification of the activation segment.

Analysis of 1- and 3-Phosphohistidine (pHis) Protein Modification Using Model Enzymes Expressed in Bacteria

Despite the discovery of protein histidine (His) phosphorylation nearly six decades ago, difficulties in measuring and quantifying this unstable post-translational modification (PTM) have limited its mechanistic analysis in prokaryotic and eukaryotic signaling. Here, we describe reliable procedures for affinity purification, cofactor-binding analysis and antibody-based detection of phosphohistidine (pHis), on the putative human His kinases NME1 (NDPK-A) and NME2 (NDPK-B) and the glycolytic phosphoglycerate mutase PGAM1. By exploiting isomer-specific monoclonal N1-pHis and N3-pHis antibodies, we describe robust protocols for immunological detection and isomer discrimination of site-specific pHis, including N3-pHis on His 11 of PGAM1.

Lipid-targeting pleckstrin homology domain turns its autoinhibitory face toward the TEC kinases

Bruton’s tyrosine kinase (BTK) is targeted in treatment of immune cancers. As patients experience drug resistance, there is a need for alternative approaches to inhibit BTK. Other recently published findings clarify the role of the BTK pleckstrin homology (PH) domain in mediating activation via dimerization and sensing of ligand concentration at the membrane. Work presented here provides insight into the autoinhibitory BTK structure that has so far been elusive via crystallographic methods. In the resting state, the BTK PH domain binds to the activation loop face of the kinase domain and allosterically alters key sites within the kinase domain. The findings define a new regulatory site, the PH/kinase interface, that can be exploited in drug discovery efforts.

The pleckstrin homology (PH) domain is well known for its phospholipid targeting function. The PH-TEC homology (PHTH) domain within the TEC family of tyrosine kinases is also a crucial component of the autoinhibitory apparatus. The autoinhibitory surface on the PHTH domain has been previously defined, and biochemical investigations have shown that PHTH-mediated inhibition is mutually exclusive with phosphatidylinositol binding. Here we use hydrogen/deuterium exchange mass spectrometry, nuclear magnetic resonance (NMR), and evolutionary sequence comparisons to map where and how the PHTH domain affects the Bruton’s tyrosine kinase (BTK) domain. The data map a PHTH-binding site on the activation loop face of the kinase C lobe, suggesting that the PHTH domain masks the activation loop and the substrate-docking site. Moreover, localized NMR spectral changes are observed for non–surface-exposed residues in the active site and on the distal side of the kinase domain. These data suggest that the association of PHTH induces allosteric conformational shifts in regions of the kinase domain that are critical for catalysis. Through statistical comparisons of diverse tyrosine kinase sequences, we identify residues unique to BTK that coincide with the experimentally determined PHTH-binding surface on the kinase domain. Our data provide a more complete picture of the autoinhibitory conformation adopted by full-length TEC kinases, creating opportunities to target the regulatory domains to control the function of these kinases in a biological setting.

Tracing the origin and evolution of pseudokinases across the tree of life

Protein phosphorylation by eukaryotic protein kinases (ePKs) is a fundamental mechanism of cell signaling in all organisms. In model vertebrates, ~10% of ePKs are classified as pseudokinases, which have amino acid changes within the catalytic machinery of the kinase domain that distinguish them from their canonical kinase counterparts. However, pseudokinases still regulate various signaling pathways, usually doing so in the absence of their own catalytic output. To investigate the prevalence, evolutionary relationships, and biological diversity of these pseudoenzymes, we performed a comprehensive analysis of putative pseudokinase sequences in available eukaryotic, bacterial, and archaeal proteomes. We found that pseudokinases are present across all domains of life, and we classified nearly 30,000 eukaryotic, 1500 bacterial, and 20 archaeal pseudokinase sequences into 86 pseudokinase families, including ~30 families that were previously unknown. We uncovered a rich variety of pseudokinases with notable expansions not only in animals but also in plants, fungi, and bacteria, where pseudokinases have previously received cursory attention. These expansions are accompanied by domain shuffling, which suggests roles for pseudokinases in plant innate immunity, plant-fungal interactions, and bacterial signaling. Mechanistically, the ancestral kinase fold has diverged in many distinct ways through the enrichment of unique sequence motifs to generate new families of pseudokinases in which the kinase domain is repurposed for noncanonical nucleotide binding or to stabilize unique, inactive kinase conformations. We further provide a collection of annotated pseudokinase sequences in the Protein Kinase Ontology (ProKinO) as a new mineable resource for the signaling community.

New tools for carbohydrate sulfation analysis: heparan sulfate 2-O-sulfotransferase (HS2ST) is a target for small-molecule protein kinase inhibitors

Sulfation of carbohydrate residues occurs on a variety of glycans destined for secretion, and this modification is essential for efficient matrix-based signal transduction. Heparan sulfate (HS) glycosaminoglycans control physiological functions ranging from blood coagulation to cell proliferation. HS biosynthesis involves membrane-bound Golgi sulfotransferases, including HS 2-O-sulfotransferase (HS2ST), which transfers sulfate from the cofactor PAPS (3'-phosphoadenosine 5'-phosphosulfate) to the 2-O position of α-l-iduronate in the maturing polysaccharide chain. The current lack of simple non-radioactive enzyme assays that can be used to quantify the levels of carbohydrate sulfation hampers kinetic analysis of this process and the discovery of HS2ST inhibitors. In the present paper, we describe a new procedure for thermal shift analysis of purified HS2ST. Using this approach, we quantify HS2ST-catalysed oligosaccharide sulfation using a novel synthetic fluorescent substrate and screen the Published Kinase Inhibitor Set, to evaluate compounds that inhibit catalysis. We report the susceptibility of HS2ST to a variety of cell-permeable compounds in vitro, including polyanionic polar molecules, the protein kinase inhibitor rottlerin and oxindole-based RAF kinase inhibitors. In a related study, published back-to-back with the present study, we demonstrated that tyrosyl protein sulfotranferases are also inhibited by a variety of protein kinase inhibitors. We propose that appropriately validated small-molecule compounds could become new tools for rapid inhibition of glycan (and protein) sulfation in cells, and that protein kinase inhibitors might be repurposed or redesigned for the specific inhibition of HS2ST.

New tools for evaluating protein tyrosine sulfation: tyrosylprotein sulfotransferases (TPSTs) are novel targets for RAF protein kinase inhibitors

Protein tyrosine sulfation is a post-translational modification best known for regulating extracellular protein–protein interactions. Tyrosine sulfation is catalysed by two Golgi-resident enzymes termed tyrosylprotein sulfotransferases (TPSTs) 1 and 2, which transfer sulfate from the cofactor PAPS (3′-phosphoadenosine 5′-phosphosulfate) to a context-dependent tyrosine in a protein substrate. A lack of quantitative tyrosine sulfation assays has hampered the development of chemical biology approaches for the identification of small-molecule inhibitors of tyrosine sulfation. In the present paper, we describe the development of a non-radioactive mobility-based enzymatic assay for TPST1 and TPST2, through which the tyrosine sulfation of synthetic fluorescent peptides can be rapidly quantified. We exploit ligand binding and inhibitor screens to uncover a susceptibility of TPST1 and TPST2 to different classes of small molecules, including the anti-angiogenic compound suramin and the kinase inhibitor rottlerin. By screening the Published Kinase Inhibitor Set, we identified oxindole-based inhibitors of the Ser/Thr kinase RAF (rapidly accelerated fibrosarcoma) as low-micromolar inhibitors of TPST1 and TPST2. Interestingly, unrelated RAF inhibitors, exemplified by the dual BRAF/VEGFR2 inhibitor RAF265, were also TPST inhibitors in vitro. We propose that target-validated protein kinase inhibitors could be repurposed, or redesigned, as more-specific TPST inhibitors to help evaluate the sulfotyrosyl proteome. Finally, we speculate that mechanistic inhibition of cellular tyrosine sulfation might be relevant to some of the phenotypes observed in cells exposed to anionic TPST ligands and RAF protein kinase inhibitors.

Genome-Wide Identification and Characterization of Tyrosine Kinases in the Silkworm, Bombyx mori

The tyrosine kinases (TKs) are important parts of metazoan signaling pathways and play significant roles in cell growth, development, apoptosis and disease. Genome-wide characterization of TKs has been conducted in many metazoans, however, systematic information about this family in Lepidoptera is still lacking. We retrieved 33 TK-encoding genes in silkworm and classified them into 25 subfamilies by sequence analysis, without members in AXL, FRK, PDGFR, STYK1 and TIE subfamilies. Although domain sequences in each subfamily are conserved, TKs in vertebrates tend to be remarkably conserved and stable. Our results of phylogenetic analysis supported the previous conclusion for the second major expansion of TK family. Gene-Ontology (GO) analysis revealed that a higher proportion of BmTKs played roles in binding, catalysis, signal transduction, metabolism, biological regulation and response to stimulus, compared to all silkworm genes annotated in GO. Moreover, the expression profile analysis of BmTKs among multiple tissues and developmental stages demonstrated that many genes exhibited stage-specific and/or sex-related expression during embryogenesis, molting and metamorphosis, and that 8 BmTKs presented tissue-specific high expression. Our study provides systematic description of silkworm tyrosine kinases, and may also provide further insights into metazoan TKs and assist future studies addressing their functions.

Coupled regulation by the juxtamembrane and sterile α motif (SAM) linker is a hallmark of ephrin tyrosine kinase evolution

Ephrin (Eph) receptor tyrosine kinases have evolutionarily diverged from other tyrosine kinases to respond to specific activation and regulatory signals that require close coupling of kinase catalytic and regulatory functions. However, the evolutionary basis for such functional coupling is not fully understood. We employed an evolutionary systems approach involving statistical mining of large sequence and structural data sets to define the hallmarks of Eph kinase evolution and functional specialization. We found that some of the most distinguishing Eph-specific residues structurally tether the flanking juxtamembrane and sterile α motif (SAM) linker regions to the kinase domain, and substitutions of these residues in EphA3 resulted in faster kinase activation. We report for the first time that the SAM domain linker is functionally coupled to the juxtamembrane through co-conserved residues in the kinase domain and that together these residues provide a structural framework for coupling catalytic and regulatory functions. The unique organization of Eph-specific tethering networks and the identification of other Eph-specific sequence features of unknown functions provide new hypotheses for future functional studies and new clues to disease mutations altering Eph kinase-specific functions.

Ensemble-based modeling and rigidity decomposition of allosteric interaction networks and communication pathways in cyclin-dependent kinases: Differentiating kinase clients of the Hsp90-Cdc37 chaperone

The overarching goal of delineating molecular principles underlying differentiation of protein kinase clients and chaperone-based modulation of kinase activity is fundamental to understanding activity of many oncogenic kinases that require chaperoning of Hsp70 and Hsp90 systems to attain a functionally competent active form. Despite structural similarities and common activation mechanisms shared by cyclin-dependent kinase (CDK) proteins, members of this family can exhibit vastly different chaperone preferences. The molecular determinants underlying chaperone dependencies of protein kinases are not fully understood as structurally similar kinases may often elicit distinct regulatory responses to the chaperone. The regulatory divergences observed for members of CDK family are of particular interest as functional diversification among these kinases may be related to variations in chaperone dependencies and can be exploited in drug discovery of personalized therapeutic agents. In this work, we report the results of a computational investigation of several members of CDK family (CDK5, CDK6, CDK9) that represented a broad repertoire of chaperone dependencies—from nonclient CDK5, to weak client CDK6, and strong client CDK9. By using molecular simulations of multiple crystal structures we characterized conformational ensembles and collective dynamics of CDK proteins. We found that the elevated dynamics of CDK9 can trigger imbalances in cooperative collective motions and reduce stability of the active fold, thus creating a cascade of favorable conditions for chaperone intervention. The ensemble-based modeling of residue interaction networks and community analysis determined how differences in modularity of allosteric networks and topography of communication pathways can be linked with the client status of CDK proteins. This analysis unveiled depleted modularity of the allosteric network in CDK9 that alters distribution of communication pathways and leads to impaired signaling in the client kinase. According to our results, these network features may uniquely define chaperone dependencies of CDK clients. The perturbation response scanning and rigidity decomposition approaches identified regulatory hotspots that mediate differences in stability and cooperativity of allosteric interaction networks in the CDK structures. By combining these synergistic approaches, our study revealed dynamic and network signatures that can differentiate kinase clients and rationalize subtle divergences in the activation mechanisms of CDK family members. The therapeutic implications of these results are illustrated by identifying structural hotspots of pathogenic mutations that preferentially target regions of the increased flexibility to enable modulation of activation changes. Our study offers a network-based perspective on dynamic kinase mechanisms and drug design by unravelling relationships between protein kinase dynamics, allosteric communications and chaperone dependencies.

Identification and characterization of tyrosine kinases in anole lizard indicate the conserved tyrosine kinase repertoire in vertebrates

The tyrosine kinases (TKs) play principal roles in regulation of multicellular aspects of the organism and are implicated in many cancer types and congenital disorders. The anole lizard has recently been introduced as a model organism for laboratory-based studies of organismal function and field studies of ecology and evolution. However, the TK family of anole lizard has not been systematically identified and characterized yet. In this study, we identified 82 TK-encoding genes in the anole lizard genome and classified them into 28 subfamilies through phylogenetic analysis, with no member from ROS and STYK1 subfamilies identified. Although TK domain sequences and domain organization in each subfamily were conserved, the total number of TKs in different species was much variable. In addition, extensive evolutionary analysis in metazoans indicated that TK repertoire in vertebrates tends to be remarkably stable. Phylogenetic analysis of Eph subfamily indicated that the divergence of EphA and EphB occurred prior to the whole genome duplication (WGD) but after the split of Urochordates and vertebrates. Moreover, the expression pattern analysis of lizard TK genes among 9 different tissues showed that 14 TK genes exhibited tissue-specific expression and 6 TK genes were widely expressed. Comparative analysis of TK expression suggested that the tissue specifically expressed genes showed different expression pattern but the widely expressed genes showed similar pattern between anole lizard and human. These results may provide insights into the evolutionary diversification of animal TK genes and would aid future studies on TK protein regulation of key growth and developmental processes.

DRP-1 is required for BH3 mimetic-mediated mitochondrial fragmentation and apoptosis

The concept of using BH3 mimetics as anticancer agents has been substantiated by the efficacy of selective drugs, such as Navitoclax and Venetoclax, in treating BCL-2-dependent haematological malignancies. However, most solid tumours depend on MCL-1 for survival, which is highly amplified in multiple cancers and a major factor determining chemoresistance. Most MCL-1 inhibitors that have been generated so far, while demonstrating early promise in vitro, fail to exhibit specificity and potency in a cellular context. To address the lack of standardised assays for benchmarking the in vitro binding of putative inhibitors before analysis of their cellular effects, we developed a rapid differential scanning fluorimetry (DSF)-based assay, and used it to screen a panel of BH3 mimetics. We next contrasted their binding signatures with their ability to induce apoptosis in a MCL-1 dependent cell line. Of all the MCL-1 inhibitors tested, only A-1210477 induced rapid, concentration-dependent apoptosis, which strongly correlated with a thermal protective effect on MCL-1 in the DSF assay. In cells that depend on both MCL-1 and BCL-XL, A-1210477 exhibited marked synergy with A-1331852, a BCL-XL specific inhibitor, to induce cell death. Despite this selectivity and potency, A-1210477 induced profound structural changes in the mitochondrial network in several cell lines that were not phenocopied following MCL-1 RNA interference or transcriptional repression, suggesting that A-1210477 induces mitochondrial fragmentation in an MCL-1-independent manner. However, A-1210477-induced mitochondrial fragmentation was dependent upon DRP-1, and silencing expression levels of DRP-1 diminished not just mitochondrial fragmentation but also BH3 mimetic-mediated apoptosis. These findings provide new insights into MCL-1 ligands, and the interplay between DRP-1 and the anti-apoptotic BCL-2 family members in the regulation of apoptosis.

Tribbles in the 21st Century: The Evolving Roles of Tribbles Pseudokinases in Biology and Disease

The Tribbles (TRIB) pseudokinases control multiple aspects of eukaryotic cell biology and evolved unique features distinguishing them from all other protein kinases. The atypical pseudokinase domain retains a regulated binding platform for substrates, which are ubiquitinated by context-specific E3 ligases. This plastic configuration has also been exploited as a scaffold to support the modulation of canonical MAPK and AKT modules. In this review, we discuss the evolution of TRIBs and their roles in vertebrate cell biology. TRIB2 is the most ancestral member of the family, whereas the emergence of TRIB3 homologs in mammals supports additional biological roles, many of which are currently being dissected. Given their pleiotropic role in diseases, the unusual TRIB pseudokinase conformation provides a highly attractive opportunity for drug design.

Pseudoenzymes are inactive counterparts of classical enzymes and have evolved in all kingdoms of life, where they regulate a vast array of biological processes. The pseudokinases are one of the best-studied families of human pseudoenzymes.

Eukaryotic TRIB pseudokinases evolved from a common ancestor (the human TRIB2 homolog), and contain a highly atypical pseudokinase domain fused to a unique docking site in an extended C tail that binds to ubiquitin E3 ligases.

TRIB evolution has led to the appearance of three mammalian TRIB pseudokinases, termed TRIB1, TRIB2, and TRIB3, which contain both unique and shared features.

In cells, TRIB pseudokinases act as modulators of substrate ubiquitination and as molecular scaffolds for the assembly and regulation of signaling modules, including the C/EBPα transcription factor and AKT and ERK networks.

TRIB1 and TRIB2 have potent oncogenic activities in vertebrate cells, and recent evidence also suggests that TRIB2 acts as a tumour suppressor, consistent with the requirement for balanced TRIB signaling in the regulation of transcription, differentiation, proliferation, and apoptosis.

cAMP-dependent protein kinase (PKA) complexes probed by complementary differential scanning fluorimetry and ion mobility-mass spectrometry

cAMP-dependent protein kinase (PKA) is an archetypal biological signaling module and a model for understanding the regulation of protein kinases. In the present study, we combine biochemistry with differential scanning fluorimetry (DSF) and ion mobility–mass spectrometry (IM–MS) to evaluate effects of phosphorylation and structure on the ligand binding, dynamics and stability of components of heteromeric PKA protein complexes in vitro. We uncover dynamic, conformationally distinct populations of the PKA catalytic subunit with distinct structural stability and susceptibility to the physiological protein inhibitor PKI. Native MS of reconstituted PKA R₂C₂ holoenzymes reveals variable subunit stoichiometry and holoenzyme ablation by PKI binding. Finally, we find that although a ‘kinase-dead’ PKA catalytic domain cannot bind to ATP in solution, it interacts with several prominent chemical kinase inhibitors. These data demonstrate the combined power of IM–MS and DSF to probe PKA dynamics and regulation, techniques that can be employed to evaluate other protein-ligand complexes, with broad implications for cellular signaling.