Protein kinases: evolution of dynamic regulatory proteins

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 investigate 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 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.

The bacterial defense system MADS interacts with CRISPR-Cas to limit phage infection and escape

The constant arms race between bacteria and their parasites has resulted in a large diversity of bacterial defenses, with many bacteria carrying multiple systems. Here, we report the discovery of a phylogenetically widespread defense system, coined methylation-associated defense system (MADS), which is distributed across gram-positive and gram-negative bacteria. MADS interacts with a CRISPR-Cas system in its native host to provide robust and durable resistance against phages. While phages can acquire epigenetic-mediated resistance against MADS, co-existence of MADS and a CRISPR-Cas system limits escape emergence. MADS comprises eight genes with predicted nuclease, ATPase, kinase, and methyltransferase domains, most of which are essential for either self/non-self discrimination, DNA restriction, or both. The complex genetic architecture of MADS and MADS-like systems, relative to other prokaryotic defenses, points toward highly elaborate mechanisms of sensing infections, defense activation, and/or interference.

Molecular Dynamics Simulation on the Suppression Mechanism of Phosphorylation to Ser222 by Allosteric Inhibitors Targeting MEK1/2 Kinase

Allosteric inhibitors of mitogen-activated protein kinase 1 (MEK1) reveal distinct interactions with MEK1 activation loop residues. The structural analyses will determine whether, and how, distinct inhibitors suppress the phosphorylation of MEK1 and may guide future therapeutic development. In this study, we explored the suppression mechanism of the phosphorylation process in the presence of MEK allosteric inhibitors, such as selumetinib, trametinib, cobimetinib, and CH5126766, by employing molecular dynamics simulations accompanied by principal component analysis. The simulations of wildtype MEK1 show that Ser222 can come close to γ-phosphate but not Ser218. We have found the conformation where Ser222 is within 5 Å of distance, which makes Ser222 accessible for γ-phosphate. The conformation analysis from the simulations of MEK1 in the presence of allosteric inhibitors reveals that the inhibitor restricts the flexibility of Ser222 through strong interactions with the activation loop, Lys97, and water mediates interactions with amino acids in the vicinity. The results reveal that all the inhibitors act as screeners between the activation loop and Mg-ATP and restricting the flexibility of the activation loop through strong interaction causes the suppression of the phosphorylation process of MEK1. The results conclude that a strong interaction of allosteric inhibitors with the activation loop restricts the movement of Ser222 toward Mg-ATP, which could be the dominant factor for the suppression of phosphorylation in MEK1. This research will provide novel insights to design effective anticancer therapeutics for targeting MEK1 in the future.

DEAD box RNA helicases are pervasive protein kinase interactors and activators

DEAD box (DDX) RNA helicases are a large family of ATPases, many of which have unknown functions. There is emerging evidence that besides their role in RNA biology, DDX proteins may stimulate protein kinases. To investigate if protein kinase-DDX interaction is a more widespread phenomenon, we conducted three orthogonal large-scale screens, including proteomics analysis with 32 RNA helicases, protein array profiling, and kinome-wide in vitro kinase assays. We retrieved Ser/Thr protein kinases as prominent interactors of RNA helicases and report hundreds of binary interactions. We identified members of ten protein kinase families, which bind to, and are stimulated by, DDX proteins, including CDK, CK1, CK2, DYRK, MARK, NEK, PRKC, SRPK, STE7/MAP2K, and STE20/PAK family members. We identified MARK1 in all screens and validated that DDX proteins accelerate the MARK1 catalytic rate. These findings indicate pervasive interactions between protein kinases and DEAD box RNA helicases, and provide a rich resource to explore their regulatory relationships.

Expression and clinical significance of U2AF homology motif kinase 1 in oral squamous cell carcinoma

OBJECTIVE

U2AF homology motif kinase 1 (UHMK1) is a newly discovered molecule that may have multiple functions. Recent studies have revealed that UHMK1 had aberrant expression in many tumors and was associated with tumor progression. However, UHMK1 was rarely reported in oral squamous cell carcinoma (OSCC).

STUDY DESIGN

In this study, Western blot, quantitative real-time polymerase chain reaction (PCR), and immunohistochemistry were used to detect the expression of UHMK1 in OSCC and peritumoral non-neoplastic tissues. Then, its relationship with clinicopathologic parameters was analyzed. The Kaplan-Meier method and Cox regression model were used to analyze the effects of UHMK1 expression on the prognosis and survival of OSCC patients.

RESULTS

Our results showed that UHMK1 had higher expression in OSCC tissues compared with in peritumoral non-neoplastic tissues, and its high expression was associated with high TNM stage and lymph node metastasis. High UHMK1 expression was related to short overall and disease-free survival times. Moreover, UHMK1 expression was identified as an independent prognostic factor that influences overall and disease-free survival of OSCC patients.

CONCLUSIONS

High expression of UHMK1 is associated with the poor prognosis of patients, and it can be used as a potential prognostic molecule for OSCC.

Allosteric activation of the co-receptor BAK1 by the EFR receptor kinase initiates immune signaling

Transmembrane signaling by plant receptor kinases (RKs) has long been thought to involve reciprocal trans-phosphorylation of their intracellular kinase domains. The fact that many of these are pseudokinase domains, however, suggests that additional mechanisms must govern RK signaling activation. Non-catalytic signaling mechanisms of protein kinase domains have been described in metazoans, but information is scarce for plants. Recently, a non-catalytic function was reported for the leucine-rich repeat (LRR)-RK subfamily XIIa member EFR (elongation factor Tu receptor) and phosphorylation-dependent conformational changes were proposed to regulate signaling of RKs with non-RD kinase domains. Here, using EFR as a model, we describe a non-catalytic activation mechanism for LRR-RKs with non-RD kinase domains. EFR is an active kinase, but a kinase-dead variant retains the ability to enhance catalytic activity of its co-receptor kinase BAK1/SERK3 (brassinosteroid insensitive 1-associated kinase 1/somatic embryogenesis receptor kinase 3). Applying hydrogen-deuterium exchange mass spectrometry (HDX-MS) analysis and designing homology-based intragenic suppressor mutations, we provide evidence that the EFR kinase domain must adopt its active conformation in order to activate BAK1 allosterically, likely by supporting αC-helix positioning in BAK1. Our results suggest a conformational toggle model for signaling, in which BAK1 first phosphorylates EFR in the activation loop to stabilize its active conformation, allowing EFR in turn to allosterically activate BAK1.

Expression, Purification and Biophysical Characterisation of Klebsiella Pneumoniae Protein Adenylyltransferase: A Systematic Integration of Empirical and Computational Modelling Approaches

Infections that are acquired due to a prolonged hospital stay and manifest 2 days following the admission of a patient to a health-care institution can be classified as hospital-acquired infections. Klebsiella pneumoniae (K. pneumoniae) has become a critical pathogen, posing serious concern globally due to the rising incidences of hypervirulent and carbapenem-resistant strains. Glutaredoxin is a redox protein that protects cells from oxidative stress as it associates with glutathione to reduce mixed disulfides. Protein adenylyltransferase (PrAT) is a pseudokinase with a proposed mechanism of transferring an AMP group from ATP to glutaredoxin. Inducing oxidative stress to the bacterium by inhibiting the activity of PrAT is a promising approach to combating its contribution to hospital-acquired infections. Thus, this study aims to overexpress, purify, and analyse the effects of ATP and Mg2+ binding to Klebsiella pneumoniae PrAT (KpPrAT). The pET expression system and nickel affinity chromatography were effective in expressing and purifying KpPrAT. Far-UV CD spectroscopy demonstrates that the protein is predominantly α-helical, even in the presence of Mg2+. Extrinsic fluorescence spectroscopy with ANS indicates the presence of a hydrophobic pocket in the presence of ATP and Mg2+, while mant-ATP studies allude to the potential nucleotide binding ability of KpPrAT. The presence of Mg2+ increases the thermostability of the protein. Isothermal titration calorimetry provides insight into the binding affinity and thermodynamic parameters associated with the binding of ATP to KpPrAT, with or without Mg2+. Conclusively, the presence of Mg2+ induces a conformation in KpPrAT that favours nucleotide binding.

Regulation of RIPK1 Phosphorylation: Implications for Inflammation, Cell Death, and Therapeutic Interventions

Receptor-interacting protein kinase 1 (RIPK1) plays a crucial role in controlling inflammation and cell death. Its function is tightly controlled through post-translational modifications, enabling its dynamic switch between promoting cell survival and triggering cell death. Phosphorylation of RIPK1 at various sites serves as a critical mechanism for regulating its activity, exerting either activating or inhibitory effects. Perturbations in RIPK1 phosphorylation status have profound implications for the development of severe inflammatory diseases in humans. This review explores the intricate regulation of RIPK1 phosphorylation and dephosphorylation and highlights the potential of targeting RIPK1 phosphorylation as a promising therapeutic strategy for mitigating human diseases.

Investigating Pyruvate Dehydrogenase Kinase 3 Inhibitory Potential of Myricetin Using Integrated Computational and Spectroscopic Approaches

Protein kinases are involved in various diseases and currently represent potential targets for drug discovery. These kinases play major roles in regulating the cellular machinery and control growth, homeostasis, and cell signaling. Dysregulation of kinase expression is associated with various disorders such as cancer and neurodegeneration. Pyruvate dehydrogenase kinase 3 (PDK3) is implicated in cancer therapeutics as a potential drug target. In this current study, a molecular docking exhibited a strong binding affinity of myricetin to PDK3. Further, a 100 ns all-atom molecular dynamics (MD) simulation study provided insights into the structural dynamics and stability of the PDK3-myricetin complex, revealing the formation of a stable complex with minimal structural alterations upon ligand binding. Additionally, the actual affinity was ascertained by fluorescence binding studies, and myricetin showed appreciable binding affinity to PDK3. Further, the kinase inhibition assay suggested significant inhibition of PDK3 by myricetin, revealing an excellent inhibitory potential with an IC50 value of 3.3 μM. In conclusion, this study establishes myricetin as a potent PDK3 inhibitor that can be implicated in therapeutic targeting cancer and PDK3-associated diseases. In addition, this study underscores the efficacy of myricetin as a potential lead to drug discovery and provides valuable insights into the inhibition mechanism, enabling advancements in cancer therapeutics.

Functional and Structural Characterization of Clinical-Stage Janus Kinase 2 Inhibitors Identifies Determinants for Drug Selectivity

Janus kinase 2 (JAK2) plays a critical role in orchestrating hematopoiesis, and its deregulation leads to various blood disorders, most importantly myeloproliferative neoplasms (MPNs). Ruxolitinib, fedratinib, momelotinib, and pacritinib are FDA-/EMA-approved JAK inhibitors effective in relieving symptoms in MPN patients but show variable clinical profiles due to poor JAK selectivity. The development of next-generation JAK2 inhibitors is hampered by the lack of comparative functional analysis and knowledge of the molecular basis of their selectivity. Here, we provide mechanistic profiling of the four approved and six clinical-stage JAK2 inhibitors and connect selectivity data with high-resolution structural and thermodynamic analyses. All of the JAK inhibitors potently inhibited JAK2 activity. Inhibitors differed in their JAK isoform selectivity and potency for erythropoietin signaling, but their general cytokine inhibition signatures in blood cells were comparable. Structural data indicate that high potency and moderate JAK2 selectivity can be obtained by targeting the front pocket of the adenosine 5'-triphosphate-binding site.

Integration of a Randomized Sequence Scanning Approach in AlphaFold2 and Local Frustration Profiling of Conformational States Enable Interpretable Atomistic Characterization of Conformational Ensembles and Detection of Hidden Allosteric States in the ABL1 Protein Kinase

Despite the success of AlphaFold methods in predicting single protein structures, these methods showed intrinsic limitations in the characterization of multiple functional conformations of allosteric proteins. The recent NMR-based structural determination of the unbound ABL kinase in the active state and discovery of the inactive low-populated functional conformations that are unique for ABL kinase present an ideal challenge for the AlphaFold2 approaches. In the current study, we employ several adaptations of the AlphaFold2 methodology to predict protein conformational ensembles and allosteric states of the ABL kinase including randomized alanine sequence scanning combined with the multiple sequence alignment subsampling proposed in this study. We show that the proposed new AlphaFold2 adaptation combined with local frustration profiling of conformational states enables accurate prediction of the protein kinase structures and conformational ensembles, also offering a robust approach for interpretable characterization of the AlphaFold2 predictions and detection of hidden allosteric states. We found that the large high frustration residue clusters are uniquely characteristic of the low-populated, fully inactive ABL form and can define energetically frustrated cracking sites of conformational transitions, presenting difficult targets for AlphaFold2. The results of this study uncovered previously unappreciated fundamental connections between local frustration profiles of the functional allosteric states and the ability of AlphaFold2 methods to predict protein structural ensembles of the active and inactive states. This study showed that integration of the randomized sequence scanning adaptation of AlphaFold2 with a robust landscape-based analysis allows for interpretable atomistic predictions and characterization of protein conformational ensembles, providing a physical basis for the successes and limitations of current AlphaFold2 methods in detecting functional allosteric states that play a significant role in protein kinase regulation.

Gene-edited protein kinases and phosphatases in molecular plant breeding

Protein phosphorylation, the most common and essential post-translational modification, belongs to crucial regulatory mechanisms in plants, affecting their metabolism, intracellular transport, cytoarchitecture, cell division, growth, development, and interactions with the environment. Protein kinases and phosphatases, two important families of enzymes optimally regulating phosphorylation, have now become important targets for gene editing in crops. We review progress on gene-edited protein kinases and phosphatases in crops using clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (Cas9). We also provide guidance for computational prediction of alterations and/or changes in function, activity, and binding of protein kinases and phosphatases as consequences of CRISPR/Cas9-based gene editing with its possible application in modern crop molecular breeding towards sustainable agriculture.

Prediction of Conformational Ensembles and Structural Effects of State-Switching Allosteric Mutants in the Protein Kinases Using Comparative Analysis of AlphaFold2 Adaptations with Sequence Masking and Shallow Subsampling

Despite the success of AlphaFold2 approaches in predicting single protein structures, these methods showed intrinsic limitations in predicting multiple functional conformations of allosteric proteins and have been challenged to accurately capture of the effects of single point mutations that induced significant structural changes. We systematically examined several implementations of AlphaFold2 methods to predict conformational ensembles for state-switching mutants of the ABL kinase. The results revealed that a combination of randomized alanine sequence masking with shallow multiple sequence alignment subsampling can significantly expand the conformational diversity of the predicted structural ensembles and capture shifts in populations of the active and inactive ABL states. Consistent with the NMR experiments, the predicted conformational ensembles for M309L/L320I and M309L/H415P ABL mutants that perturb the regulatory spine networks featured the increased population of the fully closed inactive state. On the other hand, the predicted conformational ensembles for the G269E/M309L/T334I and M309L/L320I/T334I triple ABL mutants that share activating T334I gate-keeper substitution are dominated by the active ABL form. The proposed adaptation of AlphaFold can reproduce the experimentally observed mutation-induced redistributions in the relative populations of the active and inactive ABL states and capture the effects of regulatory mutations on allosteric structural rearrangements of the kinase domain. The ensemble-based network analysis complemented AlphaFold predictions by revealing allosteric mediating centers that often directly correspond to state-switching mutational sites or reside in their immediate local structural proximity, which may explain the global effect of regulatory mutations on structural changes between the ABL states. This study suggested that attention-based learning of long-range dependencies between sequence positions in homologous folds and deciphering patterns of allosteric interactions may further augment the predictive abilities of AlphaFold methods for modeling of alternative protein sates, conformational ensembles and mutation-induced structural transformations.

Highlighting the Major Role of Cyclin C in Cyclin-Dependent Kinase 8 Activity through Molecular Dynamics Simulations

Dysregulation of cyclin-dependent kinase 8 (CDK8) activity has been associated with many diseases, including colorectal and breast cancer. As usual in the CDK family, the activity of CDK8 is controlled by a regulatory protein called cyclin C (CycC). But, while human CDK family members are generally activated in two steps, that is, the binding of the cyclin to CDK and the phosphorylation of a residue in the CDK activation loop, CDK8 does not require the phosphorylation step to be active. Another peculiarity of CDK8 is its ability to be associated with CycC while adopting an inactive form. These specificities raise the question of the role of CycC in the complex CDK8-CycC, which appears to be more complex than the other members of the CDK family. Through molecular dynamics (MD) simulations and binding free energy calculations, we investigated the effect of CycC on the structure and dynamics of CDK8. In a second step, we particularly focused our investigation on the structural and molecular basis of the protein-protein interaction between the two partners by finely analyzing the energetic contribution of residues and simulating the transition between the active and the inactive form. We found that CycC has a stabilizing effect on CDK8, and we identified specific interaction hotspots within its interaction surface compared to other human CDK/Cyc pairs. Targeting these specific interaction hotspots could be a promising approach in terms of specificity to effectively disrupt the interaction between CDK8. The simulation of the conformational transition from the inactive to the active form of CDK8 suggests that the residue Glu99 of CycC is involved in the orientation of three conserved arginines of CDK8. Thus, this residue may assume the role of the missing phosphorylation step in the activation mechanism of CDK8. In a more general view, these results point to the importance of keeping the CycC in computational studies when studying the human CDK8 protein in both the active and the inactive form.

Computational Analysis of Activation of Dimerized Epidermal Growth Factor Receptor Kinase Using the String Method and Markov State Model

Epidermal growth factor receptor (EGFR) activation is accompanied by dimerization. During the activation of the intracellular kinase domain, two EGFR kinases form an asymmetric dimer, and one side of the dimer (receiver) is activated. Using the string method and Markov state model (MSM), we performed a computational analysis of the structural changes in the activation of the EGFR dimer in this study. The string method reveals the minimum free-energy pathway (MFEP) from the inactive to active structure. The MSM was constructed from numerous trajectories of molecular dynamics simulations around the MFEP, which revealed the free-energy map of structural changes. In the activation of the receiver kinase, the unfolding of the activation loop (A-loop) is followed by the rearrangement of the C-helix, as observed in other kinases. However, unlike other kinases, the free-energy map of EGFR at the asymmetric dimer showed that the active state yielded the highest stability and revealed how interactions at the dimer interface induced receiver activation. As the H-helix of the activator approaches the C-helix of the receiver during activation, the A-loop unfolds. Subsequently, L782 of the receiver enters the pocket between the G- and H-helices of the activator, leading to a rearrangement of the hydrophobic residues around L782 of the receiver, which constitutes a structural rearrangement of the C-helix of the receiver from an outward to an inner position. The MSM analysis revealed long-time scale trajectories via kinetic Monte Carlo.

Allosteric activation of the co-receptor BAK1 by the EFR receptor kinase initiates immune signaling

Transmembrane signaling by plant receptor kinases (RKs) has long been thought to involve reciprocal trans-phosphorylation of their intracellular kinase domains. The fact that many of these are pseudokinase domains, however, suggests that additional mechanisms must govern RK signaling activation. Non-catalytic (pseudo)kinase signaling mechanisms have been described in metazoans, but information is scarce for plants. Recently, a non-catalytic function was reported for the leucine-rich repeat (LRR)-RK subfamily XIIa member EFR (ELONGATION FACTOR TU RECEPTOR) and phosphorylation-dependent conformational changes were proposed to regulate signaling of RKs with non-RD kinase domains. Here, using EFR as a model, we describe a non-catalytic activation mechanism for LRR-RKs with non-RD kinase domains. EFR is an active kinase, but a kinase-dead variant retains the ability to enhance catalytic activity of its co-receptor kinase BAK1/SERK3 (BRASSINOSTEROID INSENSITIVE 1-ASSOCIATED KINASE 1/SOMATIC EMBRYOGENESIS RECEPTOR KINASE 3). Applying hydrogen-deuterium exchange mass spectrometry (HDX-MS) analysis and designing homology-based intragenic suppressor mutations, we provide evidence that the EFR kinase domain must adopt its active conformation in order to activate BAK1 allosterically, likely by supporting αC-helix positioning in BAK1. Our results suggest a conformational toggle model for signaling, in which BAK1 first phosphorylates EFR in the activation loop to stabilize its active conformation, allowing EFR in turn to allosterically activate BAK1.

Conserved regulatory motifs in the juxtamembrane domain and kinase N-lobe revealed through deep mutational scanning of the MET receptor tyrosine kinase domain

MET is a receptor tyrosine kinase (RTK) responsible for initiating signaling pathways involved in development and wound repair. MET activation relies on ligand binding to the extracellular receptor, which prompts dimerization, intracellular phosphorylation, and recruitment of associated signaling proteins. Mutations, which are predominantly observed clinically in the intracellular juxtamembrane and kinase domains, can disrupt typical MET regulatory mechanisms. Understanding how juxtamembrane variants, such as exon 14 skipping (METΔEx14), and rare kinase domain mutations can increase signaling, often leading to cancer, remains a challenge. Here, we perform a parallel deep mutational scan (DMS) of the MET intracellular kinase domain in two fusion protein backgrounds: wild type and METΔEx14. Our comparative approach has revealed a critical hydrophobic interaction between a juxtamembrane segment and the kinase αC-helix, pointing to potential differences in regulatory mechanisms between MET and other RTKs. Additionally, we have uncovered a β5 motif that acts as a structural pivot for the kinase domain in MET and other TAM family of kinases. We also describe a number of previously unknown activating mutations, aiding the effort to annotate driver, passenger, and drug resistance mutations in the MET kinase domain.

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.

The α C helix is a central regulator of PKR activation

Protein kinase R (PKR) functions in the eukaryotic innate immune system as a first-line defense against viral infections. PKR binds viral dsRNA, leading to autophosphorylation and activation. In its active state, PKR can phosphorylate its primary substrate, eIF2 α , which blocks initiation of translation in the infected cell. It has been established that PKR activation occurs when the kinase domain dimerizes in a back-to-back configuration. However, the mechanism by which dimerization leads to enzymatic activation is not fully understood. Here, we investigate the structural mechanistic basis and energy landscape for PKR activation, with a focus on the α C helix - a kinase activation and signal integration hub - using all-atom equilibrium and enhanced sampling molecular dynamics simulations. By employing window-exchange umbrella sampling, we compute free energy profiles of activation which show that back-to-back dimerization stabilizes a catalytically competent conformation of PKR. Key hydrophobic residues in the homodimer interface contribute to stabilization of the α C helix in an active conformation and the position of its glutamate residue. Using linear mutual information analysis, we analyze allosteric communication connecting the protomers' N-lobes and the α C helix dimer interface with the α C helix.

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.

Dabrafenib for Pilocytic Astrocytoma With BRAF V599ins

This report highlights the first pediatric case of pilocytic astrocytoma with BRAF V599ins mutation, successfully treated with dabrafenib.

RedOx regulation of LRRK2 kinase activity by active site cysteines

Mutations of the human leucine-rich repeat kinase 2 (LRRK2) have been associated with both, idiopathic and familial Parkinson's disease (PD). Most of these pathogenic mutations are located in the kinase domain (KD) or GTPase domain of LRRK2. In this study we describe a mechanism in which protein kinase activity can be modulated by reversible oxidation or reduction, involving a unique pair of adjacent cysteines, the "CC" motif. Among all human protein kinases, only LRRK2 contains this "CC" motif (C2024 and C2025) in the Activation Segment (AS) of the kinase domain. In an approach combining site-directed mutagenesis, biochemical analyses, cell-based assays, and Gaussian accelerated Molecular Dynamics (GaMD) simulations we could attribute a role for each of those cysteines. We employed reducing and oxidizing agents with potential clinical relevance to investigate effects on kinase activity and microtubule docking. We find that each cysteine gives a distinct contribution: the first cysteine, C2024, is essential for LRRK2 protein kinase activity, while the adjacent cysteine, C2025, contributes significantly to redox sensitivity. Implementing thiolates (R-S-) in GaMD simulations allowed us to analyse how each of the cysteines in the "CC" motif interacts with its surrounding residues depending on its oxidation state. From our studies we conclude that oxidizing agents can downregulate kinase activity of hyperactive LRRK2 PD mutations and may provide promising tools for therapeutic strategies.

Phosphorylation of the compartmentalized PKA substrate TAF15 regulates RNA-protein interactions

Spatiotemporal-controlled second messengers alter molecular interactions of central signaling nodes for ensuring physiological signal transmission. One prototypical second messenger molecule which modulates kinase signal transmission is the cyclic-adenosine monophosphate (cAMP). The main proteinogenic cellular effectors of cAMP are compartmentalized protein kinase A (PKA) complexes. Their cell-type specific compositions precisely coordinate substrate phosphorylation and proper signal propagation which is indispensable for numerous cell-type specific functions. Here we present evidence that TAF15, which is implicated in the etiology of amyotrophic lateral sclerosis, represents a novel nuclear PKA substrate. In cross-linking and immunoprecipitation experiments (iCLIP) we showed that TAF15 phosphorylation alters the binding to target transcripts related to mRNA maturation, splicing and protein-binding related functions. TAF15 appears to be one of multiple PKA substrates that undergo RNA-binding dynamics upon phosphorylation. We observed that the activation of the cAMP-PKA signaling axis caused a change in the composition of a collection of RNA species that interact with TAF15. This observation appears to be a broader principle in the regulation of molecular interactions, as we identified a significant enrichment of RNA-binding proteins within endogenous PKA complexes. We assume that phosphorylation of RNA-binding domains adds another layer of regulation to binary protein-RNAs interactions with consequences to RNA features including binding specificities, localization, abundance and composition.

Molecular dynamics simulations reveal phosphorylation-induced conformational dynamics of the fibroblast growth factor receptor 1 kinase

The Fibroblast Growth Factor Receptor1 (FGFR1) kinase wields exquisite control on cell fate, proliferation, differentiation, and homeostasis. An imbalance of FGFR1 signaling leads to several pathogeneses of diseases ranging from multiple cancers to allergic and neurodegenerative disorders. In this study, we investigated the phosphorylation-induced conformational dynamics of FGFR1 in apo and ATP-bound states via all-atom molecular dynamics simulations. All simulations were performed for 2 × 2 µs. We have also investigated the energetics of the binding of ATP to FGFR1 using the molecular mechanics Poisson-Boltzmann scheme. Our study reveals that the FGFR1 kinase can reach a fully active configuration through phosphorylation and ATP binding. A 3-10 helix formation in the activation loop signifies its rearrangement leading to stability upon ATP binding. The interaction of phosphorylated tyrosine (pTyr654) with positively charged residues forms strong salt-bridge interactions, driving the compactness of the structure. The dynamic cross-correlation map reveals phosphorylation enhances correlated motions and reduces anti-correlated motions between different domains. We believe that the mechanistic understanding of large-conformational changes upon the activation of the FGFR1 kinase will aid the development of novel targeted therapeutics.Communicated by Ramaswamy H. Sarma.

Architecture and activation of human muscle phosphorylase kinase

The study of phosphorylase kinase (PhK)-regulated glycogen metabolism has contributed to the fundamental understanding of protein phosphorylation; however, the molecular mechanism of PhK remains poorly understood. Here we present the high-resolution cryo-electron microscopy structures of human muscle PhK. The 1.3-megadalton PhK α4β4γ4δ4 hexadecamer consists of a tetramer of tetramer, wherein four αβγδ modules are connected by the central β4 scaffold. The α- and β-subunits possess glucoamylase-like domains, but exhibit no detectable enzyme activities. The α-subunit serves as a bridge between the β-subunit and the γδ subcomplex, and facilitates the γ-subunit to adopt an autoinhibited state. Ca2+-free calmodulin (δ-subunit) binds to the γ-subunit in a compact conformation. Upon binding of Ca2+, a conformational change occurs, allowing for the de-inhibition of the γ-subunit through a spring-loaded mechanism. We also reveal an ADP-binding pocket in the β-subunit, which plays a role in allosterically enhancing PhK activity. These results provide molecular insights of this important kinase complex.

Conflicting Interfacial Electrostatic Interactions as a Design Principle to Modulate Long-Range Interdomain Communication

The extent and molecular basis of interdomain communication in multidomain proteins, central to understanding allostery and function, is an open question. One simple evolutionary strategy could involve the selection of either conflicting or favorable electrostatic interactions across the interface of two closely spaced domains to tune the magnitude of interdomain connectivity. Here, we study a bilobed domain FF34 from the eukaryotic p190A RhoGAP protein to explore one such design principle that mediates interdomain communication. We find that while the individual structural units in wild-type FF34 are marginally coupled, they exhibit distinct intrinsic stabilities and low cooperativity, manifesting as slow folding. The FF3-FF4 interface harbors a frustrated network of highly conserved electrostatic interactions-a charge troika-that promotes the population of multiple, decoupled, and non-native structural modes on a rugged native landscape. Perturbing this network via a charge-reversal mutation not only enhances stability and cooperativity but also dampens the fluctuations globally and speeds up the folding rate by at least an order of magnitude. Our work highlights how a conserved but nonoptimal network of interfacial electrostatic interactions shapes the native ensemble of a bilobed protein, a feature that could be exploited in designing molecular systems with long-range connectivity and enhanced cooperativity.

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.

Interpretable Atomistic Prediction and Functional Analysis of Conformational Ensembles and Allosteric States in Protein Kinases Using AlphaFold2 Adaptation with Randomized Sequence Scanning and Local Frustration Profiling

The groundbreaking achievements of AlphaFold2 (AF2) approaches in protein structure modeling marked a transformative era in structural biology. Despite the success of AF2 tools in predicting single protein structures, these methods showed intrinsic limitations in predicting multiple functional conformations of allosteric proteins and fold-switching systems. The recent NMR-based structural determination of the unbound ABL kinase in the active state and two inactive low-populated functional conformations that are unique for ABL kinase presents an ideal challenge for AF2 approaches. In the current study we employ several implementations of AF2 methods to predict protein conformational ensembles and allosteric states of the ABL kinase including (a) multiple sequence alignments (MSA) subsampling approach; (b) SPEACH_AF approach in which alanine scanning is performed on generated MSAs; and (c) introduced in this study randomized full sequence mutational scanning for manipulation of sequence variations combined with the MSA subsampling. We show that the proposed AF2 adaptation combined with local frustration mapping of conformational states enable accurate prediction of the ABL active and intermediate structures and conformational ensembles, also offering a robust approach for interpretable characterization of the AF2 predictions and limitations in detecting hidden allosteric states. We found that the large high frustration residue clusters are uniquely characteristic of the low-populated, fully inactive ABL form and can define energetically frustrated cracking sites of conformational transitions, presenting difficult targets for AF2 methods. This study uncovered previously unappreciated, fundamental connections between distinct patterns of local frustration in functional kinase states and AF2 successes/limitations in detecting low-populated frustrated conformations, providing a better understanding of benefits and limitations of current AF2-based adaptations in modeling of conformational ensembles.

The molecular interaction pattern of lenvatinib enables inhibition of wild-type or kinase-mutated FGFR2-driven cholangiocarcinoma

Fibroblast growth factor receptor (FGFR)-2 can be inhibited by FGFR-selective or non-selective tyrosine kinase inhibitors (TKIs). Selective TKIs are approved for cholangiocarcinoma (CCA) with FGFR2 fusions; however, their application is limited by a characteristic pattern of adverse events or evocation of kinase domain mutations. A comprehensive characterization of a patient cohort treated with the non-selective TKI lenvatinib reveals promising efficacy in FGFR2-driven CCA. In a bed-to-bench approach, we investigate FGFR2 fusion proteins bearing critical tumor-relevant point mutations. These mutations confer growth advantage of tumor cells and increased resistance to selective TKIs but remain intriguingly sensitive to lenvatinib. In line with clinical observations, in-silico analyses reveal a more favorable interaction pattern of lenvatinib with FGFR2, including an increased flexibility and ligand efficacy, compared to FGFR-selective TKIs. Finally, the treatment of a patient with progressive disease and a newly developed kinase mutation during therapy with a selective inhibitor results in a striking response to lenvatinib. Our in vitro, in silico, and clinical data suggest that lenvatinib is a promising treatment option for FGFR2-driven CCA, especially when insurmountable adverse reactions of selective TKIs or acquired kinase mutations occur.

Deciphering the molecular choreography of Janus kinase 2 inhibition via Gaussian accelerated molecular dynamics simulations: a dynamic odyssey

The Janus kinases (JAK) are crucial targets in drug development for several diseases. However, accounting for the impact of possible structural rearrangements on the binding of different kinase inhibitors is complicated by the extensive conformational variability of their catalytic kinase domain (KD). The dynamic KD contains mainly four prominent mobile structural motifs: the phosphate-binding loop (P-loop), the αC-helix within the N-lobe, the Asp-Phe-Gly (DFG) motif, and the activation loop (A-loop) within the C-lobe. These distinct structural orientations imply a complex signal transmission path for regulating the A-loop's flexibility and conformational preference for optimal JAK function. Nevertheless, the precise dynamical features of the JAK induced by different types of inhibitors still remain elusive. We performed comparative, microsecond-long, Gaussian accelerated molecular dynamics simulations in triplicate of three phosphorylated JAK2 systems: the KD alone, type-I ATP-competitive inhibitor (CI) bound KD in the catalytically active DFG-in conformation, and the type-II inhibitor (AI) bound KD in the catalytically inactive DFG-out conformation. Our results indicate significant conformational variations observed in the A-loop and αC helix motions upon inhibitor binding. Our studies also reveal that the DFG-out inactive conformation is characterized by the closed A-loop rearrangement, open catalytic cleft of N and C-lobe, the outward movement of the αC helix, and open P-loop states. Moreover, the outward positioning of the αC helix impacts the hallmark salt bridge formation between Lys882 and Glu898 in an inactive conformation. Finally, we compared their ligand binding poses and free energy by the MM/PBSA approach. The free energy calculations suggested that the AI's binding affinity is higher than CI against JAK2 due to an increased favorable contribution from the total non-polar interactions and the involvement of the αC helix. Overall, our study provides the structural and energetic insights crucial for developing more promising type I/II JAK2 inhibitors for treating JAK-related diseases.

Kinase-impaired BTK mutations are susceptible to clinical-stage BTK and IKZF1/3 degrader NX-2127

Increasing use of covalent and noncovalent inhibitors of Bruton's tyrosine kinase (BTK) has elucidated a series of acquired drug-resistant BTK mutations in patients with B cell malignancies. Here we identify inhibitor resistance mutations in BTK with distinct enzymatic activities, including some that impair BTK enzymatic activity while imparting novel protein-protein interactions that sustain B cell receptor (BCR) signaling. Furthermore, we describe a clinical-stage BTK and IKZF1/3 degrader, NX-2127, that can bind and proteasomally degrade each mutant BTK proteoform, resulting in potent blockade of BCR signaling. Treatment of chronic lymphocytic leukemia with NX-2127 achieves >80% degradation of BTK in patients and demonstrates proof-of-concept therapeutic benefit. These data reveal an oncogenic scaffold function of mutant BTK that confers resistance across clinically approved BTK inhibitors but is overcome by BTK degradation in patients.

Investigating the role of thymol as a promising inhibitor of pyruvate dehydrogenase kinase 3 for targeted cancer therapy

Protein kinases have emerged as major contributors to various diseases. They are currently exploited as a potential target in drug discovery because they play crucial roles in cell signaling, growth, and regulation. Their dysregulation is associated with inflammatory disorders, cancer, and neurodegenerative diseases. Pyruvate dehydrogenase kinase 3 (PDK3) has become an attractive drug target in cancer therapeutics. In the present study, we investigated the effective role of thymol in PDK3 inhibition due to the high affinity predicted through molecular docking studies. Hence, to better understand this inhibition mechanism, we carried out a 100 ns molecular dynamics (MD) simulation to analyse the dynamics and stability of the PDK3-thymol complex. The PDK3-thymol complex was stable and energetically favourable, with many intramolecular hydrogen bond interactions in the PDK3-thymol complex. Enzyme inhibition assay showed significant inhibition of PDK3 by thymol, revealing potential inhibitory action of thymol towards PDK3 (IC50 = 2.66 μM). In summary, we established thymol as one of the potential inhibitors of PDK3, proposing promising therapeutic implications for severe diseases associated with PDK3 dysregulation. This study further advances our understanding of thymol's therapeutic capabilities and potential role in cancer treatment.

Lskipk Lsatpase double mutants are necessary and sufficient for the compact plant architecture of butterhead lettuce

Lettuce, an important leafy vegetable crop worldwide, has rich variations in plant architecture. Butterhead lettuce, a popular horticultural type, has a unique plant architecture with loose leafy heads. The genetic and molecular mechanisms for such a compact plant architecture remain unclear. In this study we constructed a segregating population through crossing a butterhead cultivar and a stem lettuce cultivar. Genetic analysis identified the LsKIPK gene, which encodes a kinase, as the candidate gene controlling butterhead plant architecture. The Lskipk gene in the butterhead parent had a nonsense mutation, leading to a partial predicted protein. CRISPR/Cas9 and complementation tests verified its functions in plant architecture. We showed that the loss of function of LsKIPK is necessary but not sufficient for the butterhead plant architecture. To identify additional genes required for butterhead lettuce, we crossed a butterhead cultivar and a crisphead cultivar, both with the mutated Lskipk gene. Genetic mapping identified a new gene encoding an ATPase contributing to butterhead plant architecture. Knockout and complementation tests showed that loss of function of LsATPase is also required for the development of butterhead plant architecture. The Lskipk Lsatpase double mutation could reduce leaf size and leaf angle, leading to butterhead plant architecture. Expression and cytology analysis indicated that the loss of function of LsKIPK and LsATPase contributed to butterhead plant architecture by regulating cell wall development, a regulatory mechanism different from that for crisphead. This study provides new gene resources and theory for the breeding of the crop ideotype.

Small Molecule Tyrosine Kinase Inhibitors (TKIs) for Glioblastoma Treatment

In the last decade, many small molecules, usually characterized by heterocyclic scaffolds, have been designed and synthesized as tyrosine kinase inhibitors (TKIs). Among them, several compounds have been tested at preclinical and clinical levels to treat glioblastoma multiforme (GBM). GBM is the most common and aggressive type of cancer originating in the brain and has an unfavorable prognosis, with a median survival of 15-16 months and a 5-year survival rate of 5%. Despite recent advances in treating GBM, it represents an incurable disease associated with treatment resistance and high recurrence rates. For these reasons, there is an urgent need for the development of new pharmacological agents to fight this malignancy. In this review, we reported the compounds published in the last five years, which showed promising activity in GBM preclinical models acting as TKIs. We grouped the compounds based on the targeted kinase: first, we reported receptor TKIs and then, cytoplasmic and peculiar kinase inhibitors. For each small molecule, we included the chemical structure, and we schematized the interaction with the target for some representative compounds with the aim of elucidating the mechanism of action. Finally, we cited the most relevant clinical trials.

Interaction of PINK1 with nucleotides and kinetin

The ubiquitin kinase PINK1 accumulates on damaged mitochondria to trigger mitophagy, and PINK1 loss-of-function mutations cause early onset Parkinson's disease. Nucleotide analogs such as kinetin triphosphate (KTP) were reported to enhance PINK1 activity and may represent a therapeutic strategy for the treatment of Parkinson's disease. Here, we investigate the interaction of PINK1 with nucleotides, including KTP. We establish a cryo-EM platform exploiting the dodecamer assembly of Pediculus humanus corporis (Ph) PINK1 and determine PINK1 structures bound to AMP-PNP and ADP, revealing conformational changes in the kinase N-lobe that help establish PINK1's ubiquitin binding site. Notably, we find that KTP is unable to bind PhPINK1 or human (Hs) PINK1 due to a steric clash with the kinase "gatekeeper" methionine residue, and mutation to Ala or Gly is required for PINK1 to bind and use KTP as a phosphate donor in ubiquitin phosphorylation and mitophagy. HsPINK1 M318G can be used to conditionally uncouple PINK1 stabilization and activity on mitochondria.

Location, location, location: Protein kinase nanoclustering for optimised signalling output

Protein kinases (PKs) are proteins at the core of cellular signalling and are thereby responsible for most cellular physiological processes and their regulations. As for all intracellular proteins, PKs are subjected to Brownian thermal energy that tends to homogenise their distribution throughout the volume of the cell. To access their substrates and perform their critical functions, PK localisation is therefore tightly regulated in space and time, relying upon a range of clustering mechanisms. These include post-translational modifications, protein-protein and protein-lipid interactions, as well as liquid-liquid phase separation, allowing spatial restriction and ultimately regulating access to their substrates. In this review, we will focus on key mechanisms mediating PK nanoclustering in physiological and pathophysiological processes. We propose that PK nanoclusters act as a cellular quantal unit of signalling output capable of integration and regulation in space and time. We will specifically outline the various super-resolution microscopy approaches currently used to elucidate the composition and mechanisms driving PK nanoscale clustering and explore the pathological consequences of altered kinase clustering in the context of neurodegenerative disorders, inflammation, and cancer.

A palmitoyl transferase chemical-genetic system to map ZDHHC-specific S-acylation

The 23 human zinc finger Asp-His-His-Cys motif-containing (ZDHHC) S-acyltransferases catalyze long-chain S-acylation at cysteine residues across an extensive network of hundreds of proteins important for normal physiology or dysregulated in disease. Here we present a technology to directly map the protein substrates of a specific ZDHHC at the whole-proteome level, in intact cells. Structure-guided engineering of paired ZDHHC 'hole' mutants and 'bumped' chemically tagged fatty acid probes enabled probe transfer to specific protein substrates with excellent selectivity over wild-type ZDHHCs. Chemical-genetic systems were exemplified for five human ZDHHCs (3, 7, 11, 15 and 20) and applied to generate de novo ZDHHC substrate profiles, identifying >300 substrates and S-acylation sites for new functionally diverse proteins across multiple cell lines. We expect that this platform will elucidate S-acylation biology for a wide range of models and organisms.

Structural basis and selectivity of sulfatinib binding to FGFR and CSF-1R

Acquired drug resistance poses a challenge for single-target FGFR inhibitors, leading to the development of dual- or multi-target FGFR inhibitors. Sulfatinib is a multi-target kinase inhibitor for treating neuroendocrine tumors, selectively targeting FGFR1/CSF-1R. To elucidate the molecular mechanisms behind its binding and kinase selectivity, we determined the crystal structures of sulfatinib with FGFR1/CSF-1R. The results reveal common structural features and distinct conformational adaptability of sulfatinib in response to FGFR1/CSF-1R binding. Further biochemical and structural analyses disclose sensitivity of sulfatinib to FGFR/CSF-1R gatekeeper mutations. The insensitivity of sulfatinib to FGFR gatekeeper mutations highlights the indispensable interactions with the hydrophobic pocket for FGFR selectivity, whereas the rotatory flexibility may enable sulfatinib to overcome CSF-1RT663I. This study not only sheds light on the structural basis governing sulfatinib's FGFR/CSF-1R inhibition, but also provides valuable insights into the rational design of dual- or multi-target FGFR inhibitors with selectivity for CSF-1R and sensitivity to gatekeeper mutations.

Discovery of dual rho-associated protein kinase 1 (ROCK1)/apoptosis signal-regulating kinase 1 (ASK1) inhibitors as a novel approach for non-alcoholic steatohepatitis (NASH) treatment

In the current study we suggest a novel approach to curb non-alcoholic steatohepatitis (NASH) progression, and we suggest privileged scaffolds for the design of novel compounds for this aim. NASH is an advanced form of non-alcoholic fatty liver disease that can further progress into fibrosis, cirrhosis, and hepatocellular carcinoma. It is a widely emerging disease affecting 25% of the global population and has no current approved treatments. Protein kinases are key regulators of cellular pathways, of which, Rho-associated protein kinase 1 (ROCK1) and apoptosis signal-regulating kinase 1 (ASK1) play an important role in the progression of NASH and they stand out as promising targets for NASH therapy. Interestingly, their kinase domains are found to be similar in sequence and topology; therefore, dual inhibition of ROCK1 and ASK1 is expected to be amenable and could achieve a more favourable outcome. To reach this goal, a training set of ROCK1 and ASK1 protein structures co-crystalized with type 1 (ATP-competitive) inhibitors was constructed to manually generate receptor-based pharmacophore models representing ROCK1 and ASK1 inhibitors' common pharmacophoric features. The models produced were assessed using a test set of both ROCK1 and ASK1 actives and decoys, and their performance was evaluated using different assessment metrics. The best pharmacophore model obtained, showing a Mathew's correlation coefficient (MCC) of 0.71, was then used to screen the ZINC purchasable database retrieving 6178 hits that were filtered accordingly using several medicinal chemistry and pharmacokinetics filters returning 407 promising compounds. To confirm that these compounds are capable of binding to the target kinases, they were subjected to molecular docking simulations at both protein structures. The results were then assessed individually and filtered, setting the spotlight on various privileged scaffolds that could be exploited as the nucleus for designing novel ROCK1/ASK1 dual inhibitors.

Targeting CRAF kinase in anti-cancer therapy: progress and opportunities

The RAS/mitogen-activated protein kinase (MAPK) signaling cascade is commonly dysregulated in human malignancies by processes driven by RAS or RAF oncogenes. Among the members of the RAF kinase family, CRAF plays an important role in the RAS-MAPK signaling pathway, as well as in the progression of cancer. Recent research has provided evidence implicating the role of CRAF in the physiological regulation and the resistance to BRAF inhibitors through MAPK-dependent and MAPK-independent mechanisms. Nevertheless, the effectiveness of solely targeting CRAF kinase activity remains controversial. Moreover, the kinase-independent function of CRAF may be essential for lung cancers with KRAS mutations. It is imperative to develop strategies to enhance efficacy and minimize toxicity in tumors driven by RAS or RAF oncogenes. The review investigates CRAF alterations observed in cancers and unravels the distinct roles of CRAF in cancers propelled by diverse oncogenes. This review also seeks to summarize CRAF-interacting proteins and delineate CRAF's regulation across various cancer hallmarks. Additionally, we discuss recent advances in pan-RAF inhibitors and their combination with other therapeutic approaches to improve treatment outcomes and minimize adverse effects in patients with RAF/RAS-mutant tumors. By providing a comprehensive understanding of the multifaceted role of CRAF in cancers and highlighting the latest developments in RAF inhibitor therapies, we endeavor to identify synergistic targets and elucidate resistance pathways, setting the stage for more robust and safer combination strategies for cancer treatment.

Frequent loss-of-function mutations in the AMPK-α2 catalytic subunit suggest a tumour suppressor role in human skin cancers

The AMP-activated protein kinase (AMPK) is a sensor of cellular energy status activated by increases in AMP or ADP relative to ATP. Once activated, it phosphorylates targets that promote ATP-generating catabolic pathways or inhibit ATP-consuming anabolic pathways, helping to restore cellular energy balance. Analysis of human cancer genome studies reveals that the PRKAA2 gene (encoding the α2 isoform of the catalytic subunit) is often subject to mis-sense mutations in cancer, particularly in melanoma and non-melanoma skin cancers, where up to 70 mis-sense mutations have been documented, often accompanied by loss of the tumour suppressor NF1. Recently it has been reported that knockout of PRKAA2 in NF1-deficient melanoma cells promoted anchorage-independent growth in vitro, as well as growth as xenografts in immunodeficient mice in vivo, suggesting that AMPK-α2 can act as a tumour suppressor in that context. However, very few of the mis-sense mutations in PRKAA2 that occur in human skin cancer and melanoma have been tested to see whether they cause loss-of-function. We have addressed this by making most of the reported mutations and testing their activity when expressed in AMPK knockout cells. Of 55 different mis-sense mutations (representing 75 cases), 9 (12%) appeared to cause a total loss of activity, 18 (24%) a partial loss, 11 (15%) an increase in phenformin-stimulated kinase activity, while just 37 (49%) had no clear effect on kinase activity. This supports the idea that AMPK-α2 acts as a tumour suppressor in the context of human skin cancer.

Allosteric regulation of kinase activity in living cells

The dysregulation of protein kinases is associated with multiple diseases due to the kinases' involvement in a variety of cell signaling pathways. Manipulating protein kinase function, by controlling the active site, is a promising therapeutic and investigative strategy to mitigate and study diseases. Kinase active sites share structural similarities, making it difficult to specifically target one kinase, and allosteric control allows specific regulation and study of kinase function without directly targeting the active site. Allosteric sites are distal to the active site but coupled via a dynamic network of inter-atomic interactions between residues in the protein. Establishing an allosteric control over a kinase requires understanding the allosteric wiring of the protein. Computational techniques offer effective and inexpensive mapping of the allosteric sites on a protein. Here, we discuss the methods to map and regulate allosteric communications in proteins, and strategies to establish control over kinase functions in live cells and organisms. Protein molecules, or 'sensors,' are engineered to function as tools to control allosteric activity of the protein as these sensors have high spatiotemporal resolution and help in understanding cell phenotypes after immediate activation or inactivation of a kinase. Traditional methods used to study protein functions, such as knockout, knockdown, or mutation, cannot offer a sufficiently high spatiotemporal resolution. We discuss the modern repertoire of tools to regulate protein kinases as we enter a new era in deciphering cellular signaling and developing novel approaches to treat diseases associated with signal dysregulation.

Role of the Residue Q1919 in Increasing Kinase Activity of G2019S LRRK2 Kinase: A Computational Study

Mutations in multi-domain leucine-rich repeat kinase 2 (LRRK2) have been an interest to researchers as these mutations are associated with Parkinson's disease. G2019S mutation in LRRK2 kinase domain leads to the formation of additional hydrogen bonds by S2019 which results in stabilization of the active state of the kinase, thereby increasing kinase activity. Two additional hydrogen bonds of S2019 are reported separately. Here, a mechanistic picture of the formation of additional hydrogen bonds of S2019 with Q1919 (also with E1920) is presented using 'active' Roco4 kinase as a homology model and its relationship with the stabilization of the 'active' G2019S LRRK2 kinase. A conformational flipping of residue Q1919 was found which helped to form stable hydrogen bond with S2019 and made 'active' state more stable in G2019S LRRK2. Two different states were found within the 'active' kinase with respect to the conformational change (flipping) in Q1919. Two doubly-mutated systems, G2019S/Q1919A and G2019S/E1920 K, were studied separately to check the effect of Q1919 and E1920. For both cases, the stable S2 state was not formed, leading to a decrease in kinase activity. These results indicate that both the additional hydrogen bonds of S2019 (with Q1919 and E1920) are necessary to stabilize the active G2019S LRRK2.

Roles of RIPK1 as a stress sentinel coordinating cell survival and immunogenic cell death

Cell death and inflammation are closely linked arms of the innate immune response to combat infection and tissue malfunction. Recent advancements in our understanding of the intricate signals originating from dying cells have revealed that cell death serves as more than just an end point. It facilitates the exchange of information between the dying cell and cells of the tissue microenvironment, particularly immune cells, alerting and recruiting them to the site of disturbance. Receptor-interacting serine/threonine-protein kinase 1 (RIPK1) is emerging as a critical stress sentinel that functions as a molecular switch, governing cellular survival, inflammatory responses and immunogenic cell death signalling. Its tight regulation involves multiple layers of post-translational modifications. In this Review, we discuss the molecular mechanisms that regulate RIPK1 to maintain homeostasis and cellular survival in healthy cells, yet drive cell death in a context-dependent manner. We address how RIPK1 mutations or aberrant regulation is associated with inflammatory and autoimmune disorders and cancer. Moreover, we tease apart what is known about catalytic and non-catalytic roles of RIPK1 and discuss the successes and pitfalls of current strategies that aim to target RIPK1 in the clinic.

AAK1-like: A putative pseudokinase with potential roles in cargo uptake in bloodstream form Trypanosoma brucei parasites

Selection and internalization of cargo via clathrin-mediated endocytosis requires adaptor protein complexes. One complex, AP-2, acts during cargo selection at the plasma membrane. African trypanosomes lack all components of the AP-2 complex, except for a recently identified orthologue of the AP-2-associated protein kinase 1, AAK1. In characterized eukaryotes, AAK1 phosphorylates the μ2 subunit of the AP-2 complex to enhance cargo recognition and uptake into clathrin-coated vesicles. Here, we show that kinetoplastids encode not one, but two AAK1 orthologues: one (AAK1L2) is absent from salivarian trypanosomes, while the other (AAK1L1) lacks important kinase-specific residues in a range of trypanosomes. These AAK1L1 and AAK1L2 novelties reinforce suggestions of functional divergence in endocytic uptake within salivarian trypanosomes. Despite this, we show that AAK1L1 null mutant Trypanosoma brucei, while viable, display slowed proliferation, morphological abnormalities including swelling of the flagellar pocket, and altered cargo uptake. In summary, our data suggest an unconventional role for a putative pseudokinase during endocytosis and/or vesicular trafficking in T. brucei, independent of AP-2.

Mechanistic and evolutionary insights into isoform-specific 'supercharging' in DCLK family kinases

Catalytic signaling outputs of protein kinases are dynamically regulated by an array of structural mechanisms, including allosteric interactions mediated by intrinsically disordered segments flanking the conserved catalytic domain. The doublecortin-like kinases (DCLKs) are a family of microtubule-associated proteins characterized by a flexible C-terminal autoregulatory 'tail' segment that varies in length across the various human DCLK isoforms. However, the mechanism whereby these isoform-specific variations contribute to unique modes of autoregulation is not well understood. Here, we employ a combination of statistical sequence analysis, molecular dynamics simulations, and in vitro mutational analysis to define hallmarks of DCLK family evolutionary divergence, including analysis of splice variants within the DCLK1 sub-family, which arise through alternative codon usage and serve to 'supercharge' the inhibitory potential of the DCLK1 C-tail. We identify co-conserved motifs that readily distinguish DCLKs from all other calcium calmodulin kinases (CAMKs), and a 'Swiss Army' assembly of distinct motifs that tether the C-terminal tail to conserved ATP and substrate-binding regions of the catalytic domain to generate a scaffold for autoregulation through C-tail dynamics. Consistently, deletions and mutations that alter C-terminal tail length or interfere with co-conserved interactions within the catalytic domain alter intrinsic protein stability, nucleotide/inhibitor binding, and catalytic activity, suggesting isoform-specific regulation of activity through alternative splicing. Our studies provide a detailed framework for investigating kinome-wide regulation of catalytic output through cis-regulatory events mediated by intrinsically disordered segments, opening new avenues for the design of mechanistically divergent DCLK1 modulators, stabilizers, or degraders.

Crystal Structures Reveal Hidden Domain Mechanics in Protein Kinase A (PKA)

Cyclic-AMP-dependent protein kinase A (PKA) is a critical enzyme involved in various signaling pathways that plays a crucial role in regulating cellular processes including metabolism, gene transcription, cell proliferation, and differentiation. In this study, the mechanisms of allostery in PKA were investigated by analyzing the vast repertoire of crystal structures available in the RCSB database. From existing structures of murine and human PKA, we elucidated the conformational ensembles and protein dynamics that are altered in a ligand-dependent manner. Distance metrics to analyze conformations of the G-loop were proposed to delineate different states of PKA and were compared to existing structural metrics. Furthermore, ligand-dependent flexibility was investigated through normalized B'-factors to better understand the inherent dynamics in PKA. The presented study provides a contemporary approach to traditional methods in engaging the use of crystal structures for understanding protein dynamics. Importantly, our studies provide a deeper understanding into the conformational ensemble of PKA as the enzyme progresses through its catalytic cycle. These studies provide insights into kinase regulation that can be applied to both PKA individually and protein kinases as a class.

Structure and Function of a Class III Metal-Independent Lanthipeptide Synthetase

In bacteria, Ser/Thr protein kinase-like sequences are found as part of large multidomain polypeptides that biosynthesize lanthipeptides, a class of natural products distinguished by the presence of thioether cross-links. The kinase domain phosphorylates Ser or Thr residues in the peptide substrates. Subsequent β-elimination by a lyase domain yields electrophilic dehydroamino acids, which can undergo cyclase domain-catalyzed cyclization to yield conformationally restricted, bioactive compounds. Here, we reconstitute the biosynthetic pathway for a class III lanthipeptide from Bacillus thuringiensis NRRL B-23139, including characterization of a two-component protease for leader peptide excision. We also describe the first crystal structures of a class III lanthipeptide synthetase, consisting of the lyase, kinase, and cyclase domains, in various states including complexes with its leader peptide and nucleotide. The structure shows interactions between all three domains that result in an active conformation of the kinase domain. Biochemical analysis demonstrates that the three domains undergo movement upon binding of the leader peptide to establish interdomain allosteric interactions that stabilize this active form. These studies inform on the regulatory mechanism of substrate recognition and provide a framework for engineering of variants of biotechnological interest.

An unconventional gatekeeper mutation sensitizes inositol hexakisphosphate kinases to an allosteric inhibitor

Inositol hexakisphosphate kinases (IP6Ks) are emerging as relevant pharmacological targets because a multitude of disease-related phenotypes has been associated with their function. While the development of potent IP6K inhibitors is gaining momentum, a pharmacological tool to distinguish the mammalian isozymes is still lacking. Here, we implemented an analog-sensitive approach for IP6Ks and performed a high-throughput screen to identify suitable lead compounds. The most promising hit, FMP-201300, exhibited high potency and selectivity toward the unique valine gatekeeper mutants of IP6K1 and IP6K2, compared to the respective wild-type (WT) kinases. Biochemical validation experiments revealed an allosteric mechanism of action that was corroborated by hydrogen deuterium exchange mass spectrometry measurements. The latter analysis suggested that displacement of the αC helix, caused by the gatekeeper mutation, facilitates the binding of FMP-201300 to an allosteric pocket adjacent to the ATP-binding site. FMP-201300 therefore serves as a valuable springboard for the further development of compounds that can selectively target the three mammalian IP6Ks; either as analog-sensitive kinase inhibitors or as an allosteric lead compound for the WT kinases.

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.