Dynamic architecture of a protein kinase

Going beyond ATP binding site as a novel inhibitor design strategy for tau protein kinases in the treatment of Alzheimer's disease: A review

Alzheimer's disease (AD) is among the top mortality causing diseases worldwide. The presence of extracellular β-amyloidosis, as well as intraneuronal neurofibrillary aggregates of the abnormally hyperphosphorylated tau protein are two major characteristics of AD. Targeting protein kinases that are involved in the disease pathways has been a common approach in the fight against AD. Unfortunately, most kinase inhibitors currently available target the adenosine triphosphate (ATP)- binding site, which has proven unsuccessful due to issues with selectivity and resistance. As a result, a pressing need to find other alternative sites beyond the ATP- binding site has profoundly evolved. In this review, we will showcase some case studies of inhibitors of tau protein kinases acting beyond ATP binding site which have shown promising results in alleviating AD.

Role of Mutual Information Profile Shifts in Assessing the Pathogenicity of Mutations on Protein Functions: The Case of Pyrin Variants Associated With Familial Mediterranean Fever

This study presents a novel method to assess the pathogenicity of pyrin protein mutations by using mutual information (MI) as a measure to quantify the correlation between residue motions or fluctuations and associated changes affecting the phenotype. The concept of MI profile shift is presented to quantify changes in MI upon mutation, revealing insights into residue-residue interactions at critical positions. We apply this method to the pyrin protein variants, which are associated with an autosomal recessively inherited disease called familial Mediterranean fever (FMF) since the available tools do not help predict the pathogenicity of the most penetrant variants. We demonstrate the utility of MI profile shifts in assessing the effects of mutations on protein stability, function, and disease phenotype. The importance of MI shifts, particularly the negative shifts observed in the pyrin example, as indicators of severe functional effects is emphasized. Additionally, the exploration of potential compensatory mechanisms suggested by positive MI shifts, which are otherwise random and inconsequential, is highlighted. The study also discusses challenges in relating MI profile changes to disease severity and advocates for comprehensive analysis considering genetic, environmental, and stochastic factors. Overall, this study provides insights into the molecular mechanisms underlying the pathogenesis of FMF and offers a framework for identifying potential therapeutic targets based on MI profile changes induced by mutations.

Redox Regulation of cAMP-Dependent Protein Kinase and Its Role in Health and Disease

Protein kinase A (PKA) is a key regulator of cellular signaling that regulates key physiological processes such as metabolism, cell proliferation, and neuronal function. While its activation by the second messenger 3',5'-cyclic adenosine triphosphate (cAMP) is well characterized, recent research highlights additional regulatory mechanisms, particularly oxidative post-translational modifications, that influence PKA's structure, activity, and substrate specificity. Both the regulatory and catalytic subunits of PKA are susceptible to redox modifications, which have been shown to play important roles in the regulation of key cellular functions, including cardiac contractility, lipid metabolism, and the immune response. Likewise, redox-dependent modulation of PKA signaling has been implicated in numerous diseases, including cardiovascular disorders, diabetes, and neurodegenerative conditions, making it a potential therapeutic target. However, the mechanisms of crosstalk between redox- and PKA-dependent signaling remain poorly understood. This review examines the structural and functional regulation of PKA, with a focus on redox-dependent modifications and their impact on PKA-dependent signaling. A deeper understanding of these mechanisms may provide new strategies for targeting oxidative stress in disease and restoring balanced PKA signaling in cells.

Probing the Effect of Protein and Inhibitor Conformational Flexibility on the Reaction of Rocelitinib-Like Covalent Inhibitors of Epidermal Growth Factor Receptor. A Quantum Mechanics/Molecular Mechanics Study

Epidermal growth factor receptor (EGFR) is a tyrosine kinase and a validated target for non-small cell lung cancer (NSCLC). Drug discovery efforts on this target initially focused on traditional competitive, reversible ATP-binding site inhibitors; however, irreversible covalent binding EGFR inhibitors have become increasingly more popular. Covalent EGFR inhibitors have been developed using a range of different scaffolds, and unsurprisingly, the incorporation of an electrophilic acrylamide group can result in sizable orientation differences relative to the Cys797 nucleophile and the Asp800 general base. In this work, we report a QM/MM study aiming to better understand the aspects of covalent adduct formation, including the role of protein flexibility on chemical reactivity, the impact of electrophile location within the ATP binding site, and the impact of the acrylamide conformation (s-cis vs s-trans). We focus here on the diaminopyrimidine scaffold, as exemplified by Rocelitinib, where the electrophile is attached to its back pocket binding group. Our goal is to elucidate how electrophilic groups can be incorporated onto different inhibitor scaffolds targeting reactive active site residues. We find that irrespective of the EGFR MD conformation chosen, acrylamide, in both the s-cis or s-trans, can undergo reaction with rate-determining barriers of ∼20 kcal/mol. Interestingly, the nature of the rate-determining step for Rocelitinib-like inhibitors was found to be either direct nucleophilic attack or keto-enol tautomerization, depending on the precise protein and inhibitor conformation.

Distinct Allosteric Networks in CDK4 and CDK6 in the Cell Cycle and in Drug Resistance

Cyclin-dependent kinases 4 and 6 (CDK4 and CDK6) are key regulators of the G1-S phase transition in the cell cycle. In cancer cells, CDK6 overexpression often outcompetes CDK4 in driving cell cycle progression, contributing to resistance against CDK4/6 inhibitors (CDK4/6i). This suggests distinct functional and conformational differences between these two kinases, despite their striking structural and sequence similarities. Understanding the mechanisms that differentiate CDK4 and CDK6 is crucial, as resistance to CDK4/6i-frequently linked to CDK6 overexpression-remains a significant therapeutic challenge. Notably, CDK6 is often upregulated in CDK4/6i-resistant cancers and rapidly proliferating hematopoietic stem cells, underscoring its unique regulatory roles. We hypothesize that their distinct conformational dynamics explain their differences in phosphorylation of retinoblastoma protein, Rb, inhibitor efficacy, and cell cycle control. This leads us to question how their dissimilar conformational dynamics encode their distinct actions. To elucidate their differential activities, molecular mechanisms, and inhibitor binding, we combine biochemical assays and molecular dynamics (MD) simulations. We discover that CDK4 and CDK6 have distinct allosteric networks connecting the β3-αC loop and the G-loop. CDK6 exhibits stronger coupling and shorter path lengths between these regions, resulting in higher kinase activity upon cyclin binding and impacting inhibitor specificity. We also discover an unrecognized role of the unstructured CDK6 C-terminus, which allosterically connects and stabilizes the R-spine, facilitating slightly higher activity. Our findings bridge the gap between the structural similarity and functional divergence of CDK4 and CDK6, advancing the understanding of kinase regulation in cancer biology.

MAM kinases: physiological roles, related diseases, and therapeutic perspectives-a systematic review

Mitochondria-associated membranes (MAMs) are tethering regions amid the membranes of the endoplasmic reticulum (ER) and mitochondria. They are a lipid raft-like structure occupied by various proteins that facilitates signal transduction between the two organelles. The MAM proteome participates in cellular functions such as calcium (Ca2+) homeostasis, lipid synthesis, ER stress, inflammation, autophagy, mitophagy, and apoptosis. The human kinome is a superfamily of homologous proteins consisting of 538 kinases. MAM-associated kinases participate in the aforementioned cellular functions and act as cell fate executors. Studies have proved the dysregulated kinase interactions in MAM as an etiology for various diseases including cancer, diabetes mellitus, neurodegenerative diseases, cardiovascular diseases (CVDs), and obesity. Several small kinase inhibitory molecules have been well explored as promising drug candidates in clinical trials with an accelerating impact in the field of precision medicine. This review narrates the physiological actions, pathophysiology, and therapeutic potential of MAM-associated kinases with recent updates in the field.

Deciphering the structural consequences of R83 and R152 methylation on DNA polymerase β using molecular modeling

DNA polymerase β, a member of the X-family of DNA polymerases, undergoes complex regulations both in vitro and in vivo through various posttranslational modifications, including phosphorylation and methylation. The impact of these modifications varies depending on the specific amino acid undergoing alterations. In vitro, methylation of DNA polymerase β with the enzyme protein arginine methyltransferase 6 (PRMT6) at R83 and R152 enhances polymerase activity by improving DNA binding and processivity. Although these studies have shown that methylation improves DNA binding, the underlying mechanism of enhancement of polymerase activity in terms of structure and dynamics remains poorly understood. To address this gap, we modeled the methylated enzyme/DNA complex and conducted a microsecond-long simulation in the presence of Mg ions. Our results revealed significant structural changes induced by methylating both R83 and R152 sites in the enzyme. Specifically, these changes caused the DNA fragment to move closer to the C- and N-subdomains, forming additional hydrogen bonds. Furthermore, the cross-correlation map demonstrated that methylation enhanced long-range correlations within the domains/subdomains of DNA polymerase β, along with an increase in the linear mutual information value between the domains/subdomains and DNA fragments. The graph connectivity network also illustrated that methylation modulates the information pathway and identifies residues exhibiting long-distance coupling with the methylated sites. Our results provide an atomic-level understanding of the structural transition induced by methylation, shedding light on the mechanisms underlying the methylation-induced enhancement of activity in DNA polymerase β.

Distinct allosteric networks in CDK4 and CDK6 in the cell cycle and in drug resistance

Cyclin-dependent kinases 4 and 6 (CDK4 and CDK6) are key regulators of the G1-S phase transition in the cell cycle. In cancer cells, CDK6 overexpression often outcompetes CDK4 in driving cell cycle progression, contributing to resistance against CDK4/6 inhibitors (CDK4/6i). This suggests distinct functional and conformational differences between these two kinases, despite their striking structural and sequence similarities. Understanding the mechanisms that differentiate CDK4 and CDK6 is crucial, as resistance to CDK4/6i-frequently linked to CDK6 overexpression-remains a significant therapeutic challenge. Notably, CDK6 is often upregulated in CDK4/6i-resistant cancers and rapidly proliferating hematopoietic stem cells, underscoring its unique regulatory roles. We hypothesize that their distinct conformational dynamics explain their differences in phosphorylation of retinoblastoma protein, Rb, inhibitor efficacy, and cell cycle control. This leads us to question how their dissimilar conformational dynamics encode their distinct actions. To elucidate their differential activities, molecular mechanisms, and inhibitor binding, we combine biochemical assays and molecular dynamics (MD) simulations. We discover that CDK4 and CDK6 have distinct allosteric networks connecting the β3-αC loop and the G-loop. CDK6 exhibits stronger coupling and shorter path lengths between these regions, resulting in higher kinase activity upon cyclin binding and impacting inhibitor specificity. We also discover an unrecognized role of the unstructured CDK6 C-terminus, which allosterically connects and stabilizes the R-spine, facilitating slightly higher activity. Our findings bridge the gap between the structural similarity and functional divergence of CDK4 and CDK6, advancing the understanding of kinase regulation in cancer biology.

PIM-1L Kinase Binds to and Inactivates SRPK1: A Biochemical and Molecular Dynamics Study

SR/RS dipeptide repeats vary in both length and position, and are phosphorylated by SR protein kinases (SRPKs). PIM-1L, the long isoform of PIM-1 kinase, the splicing of which has been implicated in acute myeloid leukemia, contains a domain that consists largely of repeating SR/RS and SH/HS dipeptides (SR/SH-rich). In order to extend our knowledge on the specificity and cellular functions of SRPK1, here we investigate whether PIM-1L could act as substrate of SRPK1 by a combination of biochemical and computational approaches. Our biochemical data showed that the SR/SH-rich domain of PIM-1L was able to associate with SRPK1, yet it could not act as a substrate but, instead, inactivated the kinase. In line with our biochemical data, molecular modeling followed by a microsecond-scale all-atom molecular dynamics (MD) simulation suggests that the SR/SH-rich domain acts as a pseudo-docking peptide that binds to the same acidic docking-groove used in other SRPK1 interactions and induces inactive SRPK1 conformations. Comparative community network analysis of the MD trajectories, unraveled the dynamic architecture of apo SRPK1 and notable alterations of allosteric communications upon PIM-1L peptide binding. This analysis also allowed us to identify key SRPK1 residues, including unique ones, with a pivotal role in mediating allosteric signal propagation within the kinase core. Interestingly, most of the identified amino acids correspond to cancer-associated amino acid changes, validating our results. In total, this work provides insights not only on the details of SRPK1 inhibition by the PIM-1L SR/SH-domain, but also contributes to an in-depth understanding of SRPK1 regulation.

Exploring the conformational space of ROS1 kinase domain and the impact of allosteric mutations

Chromosomal rearrangements are common oncogenic events in Non-Small Cell Lung Cancer. An example is the fusion of the ROS1 kinase domain with extracellular receptors. Although the fusion leads to a target that is druggable with multi-kinase inhibitors, several reports indicate the emergence of point mutations leading to drug resistance. Although these mutations are often located in the ATP binding pocket, a subset of them is neighboring the pocket without a direct effect on drug binding. Due to the clinical impact of these allosteric mutations, there is an urge to identify the mechanism of resistance and characterize the pocket for further drug design studies. This study aimed to unravel the resistance mechanism of L1982F and S1986F/Y mutations. The variants were modeled and simulated using classical Molecular Dynamics simulations and accessed for their conformational flexibility. Our results indicate a direct effect of these allosteric mutants in the binding pocket volume with an indication of the G-loop playing a central role.

A comprehensive review of targeting RAF kinase in cancer

RAF kinases, particularly the BRAF isoform, play a crucial role in the MAPK/ERK signaling pathway, regulating key cellular processes such as proliferation, differentiation, and survival. Dysregulation of this pathway often caused by mutations in the BRAF gene or alterations in upstream regulators like Ras and receptor tyrosine kinases contributes significantly to cancer development. Mutations, such as BRAF-V600E, are present in a variety of malignancies, with the highest prevalence in melanoma. Targeted therapies against RAF kinases have achieved substantial success, especially in BRAF-V600E-mutant melanomas, where inhibitors like vemurafenib and dabrafenib have demonstrated remarkable efficacy, leading to improved patient outcomes. These inhibitors have also shown clinical benefits in cancers such as thyroid and colorectal carcinoma, although to a lesser extent. Despite these successes, therapeutic resistance remains a major hurdle. Resistance mechanisms, including RAF dimerization, feedback reactivation of the MAPK pathway, and paradoxical activation of ERK signaling, often lead to diminished efficacy over time, resulting in disease progression or even secondary malignancies. In response, current research is focusing on novel therapeutic strategies, including combination therapies that target multiple components of the pathway simultaneously, such as MEK inhibitors used in tandem with RAF inhibitors. Additionally, next-generation RAF inhibitors are being developed to address resistance and enhance therapeutic specificity. This review discusses the clinical advancements in RAF-targeted therapies, with a focus on ongoing efforts to overcome therapeutic resistance and enhance outcomes for cancer patients. It also underscores the persistent challenges in effectively targeting RAF kinase in oncology.

The Role of Protonation in the PfMATE Transporter Protein Structural Transitions

Multi-antimicrobial extrusion (MATE) transporter membrane proteins provide drug and toxin resistivity by expelling compounds from cells. MATE proteins can be pictured as V-shaped. To regulate its functioning, the protein structure can switch between outward-facing (OF) and inward-facing (IF). Pyrococcus furiosus MATE (PfMATE) is the only member of the multidrug/oligosaccharidyl-lipid/polysaccharide (MOP) superfamily that has available both the IF and OF crystal structures. With the availability of both the IF and OF structures, we are able to perform computational investigations to determine how protonation of specific amino acids causes a cascade of changes in the protein conformation that allow PfMATE to change its state from OF to IF in order to regulate its antiporter function. Using a variety of computational and theoretical techniques, we investigated four different systems of IF and OF PfMATE along with the native archaeal lipid bilayer, without or with protonation at the experimentally determined locations within the protein. We performed molecular dynamics (MD) simulations to investigate the flexibility of the four different PfMATE structures and also performed targeted molecular dynamics (TMD) simulations, during which we observed occluded conformations. Our analysis of hydrogen bond changes, potential of mean force, dynamic network analysis, and transfer entropy analysis provides information on how protonation can induce cascading structural changes responsible for the transition between the IF and OF configurations.

Role of the αC-β4 loop in protein kinase structure and dynamics

Although the αC-β4 loop is a stable feature of all protein kinases, the importance of this motif as a conserved element of secondary structure, as well as its links to the hydrophobic architecture of the kinase core, has been underappreciated. We first review the motif and then describe how it is linked to the hydrophobic spine architecture of the kinase core, which we first discovered using a computational tool, local spatial Pattern (LSP) alignment. Based on NMR predictions that a mutation in this motif abolishes the synergistic high-affinity binding of ATP and a pseudo substrate inhibitor, we used LSP to interrogate the F100A mutant. This comparison highlights the importance of the αC-β4 loop and key residues at the interface between the N- and C-lobes. In addition, we delved more deeply into the structure of the apo C-subunit, which lacks ATP. While apo C-subunit showed no significant changes in backbone dynamics of the αC-β4 loop, we found significant differences in the side chain dynamics of K105. The LSP analysis suggests disruption of communication between the N- and C-lobes in the F100A mutant, which would be consistent with the structural changes predicted by the NMR spectroscopy.

The expanding landscape of canonical and non-canonical protein phosphorylation

Protein phosphorylation is a crucial regulatory mechanism in cell signaling, acting as a molecular switch that modulates protein function. Catalyzed by protein kinases and reversed by phosphoprotein phosphatases, it is essential in both normal physiological and pathological states. Recent advances have uncovered a vast and intricate landscape of protein phosphorylation that include histidine phosphorylation and more unconventional events, such as pyrophosphorylation and polyphosphorylation. Many questions remain about the true size of the phosphoproteome and, more importantly, its site-specific functional relevance. The involvement of unconventional actors such as pseudokinases and pseudophosphatases adds further complexity to be resolved. This review explores recent discoveries and ongoing challenges, highlighting the need for continued research to fully elucidate the roles and regulation of protein phosphorylation.

Disruptions to protein kinase A localization in adrenal pathology

Cell signaling fidelity requires specificity in protein-protein interactions and precise subcellular localization of signaling molecules. In the case of protein phosphorylation, many kinases and phosphatases exhibit promiscuous substrate pairing and therefore require targeting interactions to modify the appropriate substrates and avoid cross-talk among different pathways. In the past 10 years, researchers have discovered and investigated how loss of specific interactions and subcellular targeting for the protein kinase A catalytic subunit (PKAc) lead to cortisol-producing adenomas and the debilitating stress disorder adrenal Cushing's syndrome. This article reviews classical studies regarding PKA localization in glucocorticoid-producing adrenal cells and synthesizes recent evidence of disrupted PKA localization and selective regulatory interactions in adrenal pathology.

Improving docking and virtual screening performance using AlphaFold2 multi-state modeling for kinases

Structure-based virtual screening (SBVS) is a crucial computational approach in drug discovery, but its performance is sensitive to structural variations. Kinases, which are major drug targets, exemplify this challenge due to active site conformational changes caused by different inhibitor types. Most experimentally determined kinase structures have the DFGin state, potentially biasing SBVS towards type I inhibitors and limiting the discovery of diverse scaffolds. We introduce a multi-state modeling (MSM) protocol for AlphaFold2 (AF2) kinase structures using state-specific templates to address these challenges. Our comprehensive benchmarks evaluate predicted model qualities, binding pose prediction accuracy, and hit compound identification through ensemble SBVS. Results demonstrate that MSM models exhibit comparable or improved structural accuracy compared to standard AF2 models, enhancing pose prediction accuracy and effectively capturing kinase-ligand interactions. In virtual screening experiments, our MSM approach consistently outperforms standard AF2 and AF3 modeling, particularly in identifying diverse hit compounds. This study highlights the potential of MSM in broadening kinase inhibitor discovery by facilitating the identification of chemically diverse inhibitors, offering a promising solution to the structural bias problem in kinase-targeted drug discovery.

Dimerization Promotes PKR Activation by Modulating Energetics of αC Helix Conversion between Active and Inactive Conformations

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 the 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 critical 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.

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. Here, we use Saccharomyces cerevisiae to inducibly express tyrosine kinases. Because yeast lacks bona fide tyrosine kinases, the resulting tyrosine phosphorylation is biologically spurious. We engineered 44 yeast strains each expressing a tyrosine kinase, and quantitatively analysed their phosphoproteomes. This analysis resulted in ~30,000 phosphosites mapping to ~3500 proteins. The number of spurious pY sites generated correlates strongly with decreased growth, and we predict over 1000 pY events to be deleterious. However, we also find that many 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 tyrosine kinases. Our results suggest that, alongside the risk for toxicity, the cell can tolerate a large degree of non-functional crosstalk as interaction networks evolve.

Deciphering Internal Regulatory Patterns within the p53 Core Tetramer: Insights from Community Network Analysis

Gene therapy is one of the most effective strategies for cancer treatment. The p53 protein, commonly known as the "guardian of the genome", plays a critical role in gene activation and tumor suppression. Tetramerization of the p53 core domain is an essential allosteric process that supports its suppression functions. This letter presents a framework to analyze the structure, function, and dynamic connectivity of the p53 tetramer, using community network analysis based on all-atom molecular dynamics simulations. The communities within the p53 monomer exhibit distinct functional roles, while interactions between molecules establish a symmetrical network structure. We identified direct evidence of single, double, and multiple pathway regulations within the p53 tetramer and crucial residue pairs involved in these connections. Our study provides a comprehensive framework to understand the community network of the p53 tetramer, offering new insights into the stable formation of the p53 core tetramer.

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.

Reused Protein Segments Linked to Functional Dynamics

Protein space is characterized by extensive recurrence, or "reuse," of parts, suggesting that new proteins and domains can evolve by mixing-and-matching of existing segments. From an evolutionary perspective, for a given combination to persist, the protein segments should presumably not only match geometrically but also dynamically communicate with each other to allow concerted motions that are key to function. Evidence from protein space supports the premise that domains indeed combine in this manner; we explore whether a similar phenomenon can be observed at the sub-domain level. To this end, we use Gaussian Network Models (GNMs) to calculate the so-called soft modes, or low-frequency modes of motion for a dataset of 150 protein domains. Modes of motion can be used to decompose a domain into segments of consecutive amino acids that we call "dynamic elements", each of which belongs to one of two parts that move in opposite senses. We find that, in many cases, the dynamic elements, detected based on GNM analysis, correspond to established "themes": Sub-domain-level segments that have been shown to recur in protein space, and which were detected in previous research using sequence similarity alone (i.e. completely independently of the GNM analysis). This statistically significant correlation hints at the importance of dynamics in evolution. Overall, the results are consistent with an evolutionary scenario where proteins have emerged from themes that need to match each other both geometrically and dynamically, e.g. to facilitate allosteric regulation.

Accurate estimation of the normalized mutual information of multidimensional data

While the linear Pearson correlation coefficient represents a well-established normalized measure to quantify the inter-relation of two stochastic variables X and Y, it fails for multidimensional variables, such as Cartesian coordinates. Avoiding any assumption about the underlying data, the mutual information I(X, Y) does account for multidimensional correlations. However, unlike the normalized Pearson correlation, it has no upper bound (I ∈ [0, ∞)), i.e., it is not clear if say, I = 0.4 corresponds to a low or a high correlation. Moreover, the mutual information (MI) involves the estimation of high-dimensional probability densities (e.g., six-dimensional for Cartesian coordinates), which requires a k nearest-neighbor algorithm, such as the estimator by Kraskov et al. [Phys. Rev. E 69, 066138 (2004)]. As existing methods to normalize the MI cannot be used in connection with this estimator, a new approach is presented, which uses an entropy estimation method that is invariant under variable transformations. The algorithm is numerically efficient and does not require more effort than the calculation of the (un-normalized) MI. After validating the method by applying it to various toy models, the normalized MI between the Cα-coordinates of T4 lysozyme is considered and compared to a correlation analysis of inter-residue contacts.

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 αC-β4 loop controls the allosteric cooperativity between nucleotide and substrate in the catalytic subunit of protein kinase A

Allosteric cooperativity between ATP and substrates is a prominent characteristic of the cAMP-dependent catalytic subunit of protein kinase A (PKA-C). This long-range synergistic action is involved in substrate recognition and fidelity, and it may also regulate PKA's association with regulatory subunits and other binding partners. To date, a complete understanding of this intramolecular mechanism is still lacking. Here, we integrated NMR(Nuclear Magnetic Resonance)-restrained molecular dynamics simulations and a Markov State Model to characterize the free energy landscape and conformational transitions of PKA-C. We found that the apoenzyme populates a broad free energy basin featuring a conformational ensemble of the active state of PKA-C (ground state) and other basins with lower populations (excited states). The first excited state corresponds to a previously characterized inactive state of PKA-C with the αC helix swinging outward. The second excited state displays a disrupted hydrophobic packing around the regulatory (R) spine, with a flipped configuration of the F100 and F102 residues at the αC-β4 loop. We validated the second excited state by analyzing the F100A mutant of PKA-C, assessing its structural response to ATP and substrate binding. While PKA-CF100A preserves its catalytic efficiency with Kemptide, this mutation rearranges the αC-β4 loop conformation, interrupting the coupling of the two lobes and abolishing the allosteric binding cooperativity. The highly conserved αC-β4 loop emerges as a pivotal element to control the synergistic binding of nucleotide and substrate, explaining how mutations or insertions near or within this motif affect the function and drug sensitivity in homologous kinases.

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.

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.

Real-time single-molecule imaging of CaMKII-calmodulin interactions

The binding of calcium/calmodulin (CAM) to calcium/calmodulin-dependent protein kinase II (CaMKII) initiates an ATP-driven cascade that triggers CaMKII autophosphorylation. The autophosphorylation in turn increases the CaMKII affinity for CAM. Here, we studied the ATP dependence of CAM association with the actin-binding CaMKIIβ isoform using single-molecule total internal reflection fluorescence microscopy. Rhodamine-CAM associations/dissociations to surface-immobilized Venus-CaMKIIβ were resolved with 0.5 s resolution from video records, batch-processed with a custom algorithm. CAM occupancy was determined simultaneously with spot-photobleaching measurement of CaMKII holoenzyme stoichiometry. We show the ATP-dependent increase of the CAM association requires dimer formation for both the α and β isoforms. The study of mutant β holoenzymes revealed that the ATP-dependent increase in CAM affinity results in two distinct states. The phosphorylation-defective (T287.306-307A) holoenzyme resides only in the low-affinity state. CAM association is further reduced in the T287A holoenzyme relative to T287.306-307A. In the absence of ATP, the affinity of CAM for the T287.306-307A mutant and the wild-type monomer are comparable. The affinity of the ATP-binding impaired (K43R) mutant is even weaker. In ATP, the K43R holoenzyme resides in the low-affinity state. The phosphomimetic mutant (T287D) resides only in a 1000-fold higher-affinity state, with mean CAM occupancy of more than half of the 14-mer holoenzyme stoichiometry in picomolar CAM. ATP promotes T287D holoenzyme disassembly but does not elevate CAM occupancy. Single Poisson distributions characterized the ATP-dependent CAM occupancy of mutant holoenzymes. In contrast, the CAM occupancy of the wild-type population had a two-state distribution with both low- and high-affinity states represented. The low-affinity state was the dominant state, a result different from published in vitro assays. Differences in assay conditions can alter the balance between activating and inhibitory autophosphorylation. Bound ATP could be sufficient for CaMKII structural function, while antagonistic autophosphorylations may tune CaMKII kinase-regulated action-potential frequency decoding in vivo.

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.

Community analysis of large-scale molecular dynamics simulations elucidated dynamics-driven allostery in tyrosine kinase 2

TYK2 is a nonreceptor tyrosine kinase, member of the Janus kinases (JAK), with a central role in several diseases, including cancer. The JAKs' catalytic domains (KD) are highly conserved, yet the isolated TYK2-KD exhibits unique specificities. In a previous work, using molecular dynamics (MD) simulations of a catalytically impaired TYK2-KD variant (P1104A) we found that this amino acid change of its JAK-characteristic insert (αFG), acts at the dynamics level. Given that structural dynamics is key to the allosteric activation of protein kinases, in this study we applied a long-scale MD simulation and investigated an active TYK2-KD form in the presence of adenosine 5'-triphosphate and one magnesium ion that represents a dynamic and crucial step of the catalytic cycle, in other protein kinases. Community analysis of the MD trajectory shed light, for the first time, on the dynamic profile and dynamics-driven allosteric communications within the TYK2-KD during activation and revealed that αFG and amino acids P1104, P1105, and I1112 in particular, hold a pivotal role and act synergistically with a dynamically coupled communication network of amino acids serving intra-KD signaling for allosteric regulation of TYK2 activity. Corroborating our findings, most of the identified amino acids are associated with cancer-related missense/splice-site mutations of the Tyk2 gene. We propose that the conformational dynamics at this step of the catalytic cycle, coordinated by αFG, underlie TYK2-unique substrate recognition and account for its distinct specificity. In total, this work adds to knowledge towards an in-depth understanding of TYK2 activation and may be valuable towards a rational design of allosteric TYK2-specific inhibitors.

New pan-ALK inhibitor-resistant EML4::ALK mutations detected by liquid biopsy in lung cancer patients

ALK and ROS1 fusions are effectively targeted by tyrosine kinase inhibitors (TKIs), however patients inevitably relapse after an initial response, often due to kinase domain mutations. We investigated circulating DNA from TKI-relapsed NSCLC patients by deep-sequencing. New EML4::ALK substitutions, L1198R, C1237Y and L1196P, were identified in the plasma of NSCLC ALK patients and characterized in a Ba/F3 cell model. Variants C1237Y and L1196P demonstrated pan-inhibitor resistance across 5 clinical and 2 investigational TKIs.

Allostery in Protein Tyrosine Phosphatases is Enabled by Divergent Dynamics

Dynamics-driven allostery provides important insights into the working mechanics of proteins, especially enzymes. In this study, we employ this paradigm to answer a basic question: in enzyme superfamilies, where the catalytic mechanism, active sites, and protein fold are conserved, what accounts for the difference in the catalytic prowess of the individual members? We show that when subtle changes in sequence do not translate to changes in structure, they do translate to changes in dynamics. We use sequentially diverse PTP1B, TbPTP1, and YopH as representatives of the conserved protein tyrosine phosphatase (PTP) superfamily. Using amino acid network analysis of group behavior (community analysis) and influential node dominance on networks (eigenvector centrality), we explain the dynamic basis of the catalytic variations seen between the three proteins. Importantly, we explain how a dynamics-based blueprint makes PTP1B amenable to allosteric control and how the same is abstracted in TbPTP1 and YopH.

Dynamic correlations: exact and approximate methods for mutual information

MOTIVATION

Proteins are dynamic entities that undergo conformational changes critical for their functions. Understanding the communication pathways and information transfer within proteins is crucial for elucidating allosteric interactions in their mechanisms. This study utilizes mutual information (MI) analysis to probe dynamic allostery. Using two cases, Ubiquitin and PLpro, we have evaluated the accuracy and limitations of different approximations including the exact anisotropic and isotropic models, multivariate Gaussian model, isotropic Gaussian model, and the Gaussian Network Model (GNM) in revealing allosteric interactions.

RESULTS

Our findings emphasize the required trajectory length for capturing accurate mutual information profiles. Long molecular dynamics trajectories, 1 ms for Ubiquitin and 100 µs for PLpro are used as benchmarks, assuming they represent the ground truth. Trajectory lengths of approximately 5 µs for Ubiquitin and 1 µs for PLpro marked the onset of convergence, while the multivariate Gaussian model accurately captured mutual information with trajectories of 5 ns for Ubiquitin and 350 ns for PLpro. However, the isotropic Gaussian model is less successful in representing the anisotropic nature of protein dynamics, particularly in the case of PLpro, highlighting its limitations. The GNM, however, provides reasonable approximations of long-range information exchange as a minimalist network model based on a single crystal structure. Overall, the optimum trajectory lengths for effective Gaussian approximations of long-time dynamic behavior depend on the inherent dynamics within the protein's topology. The GNM, by showcasing dynamics across relatively diverse time scales, can be used either as a standalone method or to gauge the adequacy of MD simulation lengths.

AVAILABILITY AND IMPLEMENTATION

Mutual information codes are available at https://github.com/kemaldemirtas/prc-MI.git.

Imaging of plant calcium-sensor kinase conformation monitors real time calcium-dependent decoding in planta

Changes in cytosolic calcium (Ca2+) concentration are among the earliest reactions to a multitude of stress cues. While a plethora of Ca2+-permeable channels may generate distinct Ca2+ signatures and contribute to response specificities, the mechanisms by which Ca2+ signatures are decoded are poorly understood. Here, we developed a genetically encoded Förster resonance energy transfer (FRET)-based reporter that visualizes the conformational changes in Ca2+-dependent protein kinases (CDPKs/CPKs). We focused on two CDPKs with distinct Ca2+-sensitivities, highly Ca2+-sensitive Arabidopsis (Arabidopsis thaliana) AtCPK21 and rather Ca2+-insensitive AtCPK23, to report conformational changes accompanying kinase activation. In tobacco (Nicotiana tabacum) pollen tubes, which naturally display coordinated spatial and temporal Ca2+ fluctuations, CPK21-FRET, but not CPK23-FRET, reported oscillatory emission ratio changes mirroring cytosolic Ca2+ changes, pointing to the isoform-specific Ca2+-sensitivity and reversibility of the conformational change. In Arabidopsis guard cells, CPK21-FRET-monitored conformational dynamics suggest that CPK21 serves as a decoder of signal-specific Ca2+ signatures in response to abscisic acid and the flagellin peptide flg22. Based on these data, CDPK-FRET is a powerful approach for tackling real-time live-cell Ca2+ decoding in a multitude of plant developmental and stress responses.

Development of MPS1 Inhibitors: Recent Advances and Perspectives

Monopolar spindle kinase 1 (MPS1) plays a pivotal role as a dual-specificity kinase governing spindle assembly checkpoint activation and sister chromatid separation in mitosis. Its overexpression has been observed in various human malignancies. MPS1 reduces spindle assembly checkpoint sensitivity, allowing tumor cells with a high degree of aneuploidy to complete mitosis and survive. Thus, MPS1 has emerged as a promising candidate for cancer therapy. Despite the identification of numerous MPS1 inhibitors, only five have advanced to clinical trials with none securing FDA approval for cancer treatment. In this perspective, we provide a concise overview of the structural and functional characteristics of MPS1 by highlighting its relevance to cancer. Additionally, we explore the structure-activity relationships, selectivity, and pharmacokinetics of MPS1 inhibitors featuring diverse scaffolds. Moreover, we review the reported work on enhancing MPS1 inhibitor selectivity, offering valuable insights into the discovery of novel, highly potent small-molecule MPS1 inhibitors.

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.

Protein Kinase Structure and Dynamics: Role of the αC-β4 Loop

Although the αC-β4 loop is a stable feature of all protein kinases, the importance of this motif as a conserved element of secondary structure, as well as its links to the hydrophobic architecture of the kinase core, has been underappreciated. We first review the motif and then describe how it is linked to the hydrophobic spine architecture of the kinase core, which we first discovered using a computational tool, Local Spatial Pattern (LSP) alignment. Based on NMR predictions that a mutation in this motif abolishes the synergistic high-affinity binding of ATP and a pseudo substrate inhibitor, we used LSP to interrogate the F100A mutant. This comparison highlights the importance of the αC-β4 loop and key residues at the interface between the N- and C-lobes. In addition, we delved more deeply into the structure of the apo C-subunit, which lacks ATP. While apo C-subunit showed no significant changes in backbone dynamics of the αC-β4 loop, we found significant differences in the side chain dynamics of K105. The LSP analysis suggests disruption of communication between the N- and C-lobes in the F100A mutant, which would be consistent with the structural changes predicted by the NMR spectroscopy.

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.

The αC-β4 loop controls the allosteric cooperativity between nucleotide and substrate in the catalytic subunit of protein kinase A

Allosteric cooperativity between ATP and substrates is a prominent characteristic of the cAMP-dependent catalytic (C) subunit of protein kinase A (PKA). Not only this long-range synergistic action is involved in substrate recognition and fidelity, but it is likely to regulate PKA association with regulatory subunits and other binding partners. To date, a complete understanding of the molecular determinants for this intramolecular mechanism is still lacking. Here, we used an integrated NMR-restrained molecular dynamics simulations and a Markov Model to characterize the free energy landscape and conformational transitions of the catalytic subunit of protein kinase A (PKA-C). We found that the apo-enzyme populates a broad free energy basin featuring a conformational ensemble of the active state of PKA-C (ground state) and other basins with lower populations (excited states). The first excited state corresponds to a previously characterized inactive state of PKA-C with the αC helix swinging outward. The second excited state displays a disrupted hydrophobic packing around the regulatory (R) spine, with a flipped configuration of the F100 and F102 residues at the tip of the αC-β4 loop. To experimentally validate the second excited state, we mutated F100 into alanine and used NMR spectroscopy to characterize the binding thermodynamics and structural response of ATP and a prototypical peptide substrate. While the activity of PKA-CF100A toward a prototypical peptide substrate is unaltered and the enzyme retains its affinity for ATP and substrate, this mutation rearranges the αC-β4 loop conformation interrupting the allosteric coupling between nucleotide and substrate. The highly conserved αC-β4 loop emerges as a pivotal element able to modulate the synergistic binding between nucleotide and substrate and may affect PKA signalosome. These results may explain how insertion mutations within this motif affect drug sensitivity in other homologous kinases.

The E3 ubiquitin ligase SINA1 and the protein kinase BIN2 cooperatively regulate PHR1 in apple anthocyanin biosynthesis

PHR1 (PHOSPHATE STARVATION RESPONSE1) plays key roles in the inorganic phosphate (Pi) starvation response and in Pi deficiency-induced anthocyanin biosynthesis in plants. However, the post-translational regulation of PHR1 is unclear, and the molecular basis of PHR1-mediated anthocyanin biosynthesis remains elusive. In this study, we determined that MdPHR1 was essential for Pi deficiency-induced anthocyanin accumulation in apple (Malus × domestica). MdPHR1 interacted with MdWRKY75, a positive regulator of anthocyanin biosynthesis, to enhance the MdWRKY75-activated transcription of MdMYB1, leading to anthocyanin accumulation. In addition, the E3 ubiquitin ligase SEVEN IN ABSENTIA1 (MdSINA1) negatively regulated MdPHR1-promoted anthocyanin biosynthesis via the ubiquitination-mediated degradation of MdPHR1. Moreover, the protein kinase apple BRASSINOSTEROID INSENSITIVE2 (MdBIN2) phosphorylated MdPHR1 and positively regulated MdPHR1-mediated anthocyanin accumulation by attenuating the MdSINA1-mediated ubiquitination degradation of MdPHR1. Taken together, these findings not only demonstrate the regulatory role of MdPHR1 in Pi starvation induced anthocyanin accumulation, but also provide an insight into the post-translational regulation of PHR1.

New Insights into the Cooperativity and Dynamics of Dimeric Enzymes

A survey of protein databases indicates that the majority of enzymes exist in oligomeric forms, with about half of those found in the UniProt database being homodimeric. Understanding why many enzymes are in their dimeric form is imperative. Recent developments in experimental and computational techniques have allowed for a deeper comprehension of the cooperative interactions between the subunits of dimeric enzymes. This review aims to succinctly summarize these recent advancements by providing an overview of experimental and theoretical methods, as well as an understanding of cooperativity in substrate binding and the molecular mechanisms of cooperative catalysis within homodimeric enzymes. Focus is set upon the beneficial effects of dimerization and cooperative catalysis. These advancements not only provide essential case studies and theoretical support for comprehending dimeric enzyme catalysis but also serve as a foundation for designing highly efficient catalysts, such as dimeric organic catalysts. Moreover, these developments have significant implications for drug design, as exemplified by Paxlovid, which was designed for the homodimeric main protease of SARS-CoV-2.

Allostery in Protein Tyrosine Phosphatases is Enabled by Divergent Dynamics

Dynamics-driven allostery provides important insights into the working mechanics of proteins, especially enzymes. In this study we employ this paradigm to answer a basic question: in enzyme superfamilies where the catalytic mechanism, active sites and protein fold are conserved, what accounts for the difference in the catalytic prowess of the individual members? We show that when subtle changes in sequence do not translate to changes in structure, they do translate to changes in dynamics. We use sequentially diverse PTP1B, TbPTP1, and YopH as the representatives of the conserved Protein Tyrosine Phosphatase (PTP) superfamily. Using amino acid network analysis of group behavior (community analysis) and influential node dominance on networks (eigenvector centrality), we explain the dynamic basis of catalytic variations seen between the three proteins. Importantly, we explain how a dynamics-based blueprint makes PTP1B amenable to allosteric control and how the same is abstracted in TbPTP1 and YopH.

Polypeptide Substrate Accessibility Hypothesis: Gain-of-Function R206H Mutation Allosterically Affects Activin Receptor-like Protein Kinase Activity

Although structurally similar to type II counterparts, type I or activin receptor-like kinases (ALKs) are set apart by a metastable helix-loop-helix (HLH) element preceding the protein kinase domain that, according to a longstanding paradigm, serves passive albeit critical roles as an inhibitor-to-substrate-binding switch. A single recurrent mutation in the codon of the penultimate residue, directly adjacent the position of a constitutively activating substitution, causes milder activation of ACVR1/ALK2 leading to sporadic heterotopic bone deposition in patients presenting with fibrodysplasia ossificans progressiva, or FOP. To determine the protein structural-functional basis for the gain of function, R206H mutant, Q207D (aspartate-substituted caALK2) and HLH subdomain-truncated (208 Ntrunc) forms were compared to one another and the wild-type enzyme through in vitro kinase and protein-protein interaction analyses that were complemented by signaling read-out (p-Smad) in primary mouse embryonic fibroblasts and Drosophila S2 cells. Contrary to the paradigm, the HLH subdomain actively suppressed the phosphotransferase activity of the enzyme, even in the absence of FKBP12. Unexpectedly, perturbation of the HLH subdomain elevated kinase activity at a distance, i.e., allosterically, at the ATP-binding and polypeptide-interacting active site cleft. Accessibility to polypeptide substrate (BMP Smad C-terminal tails) due to allosterically altered conformations of type I active sites within heterohexameric cytoplasmic signaling complexes-assembled noncanonically by activin-type II receptors extracellularly-is hypothesized to produce a gain of function of the R206H mutant protein responsible for episodic heterotopic ossification in FOP.

Thermodynamic understanding of flower pigmentation

We have reviewed and interpreted the thermodynamic principles for flower pigmentation. The basic thoughts are as follows: 1) any biological trait is associated with one thermodynamic system; 2) a thermodynamic system of biology cannot be physically isolated from complex thermal systems of biology but can be separately studied using thermodynamic methods; 3) a thermodynamic system of biology has all types of information, including volume, shape, and structure, unlike the traditional thermal system of gases; 4) a thermodynamic system of biology is associated with one type of biological structure that is not fully fixed but can change its conformation under different conditions; and 5) a thermodynamic system of biology shows a hierarchical structure. On the basis of these principles, several conclusions regarding flower pigmentation are obtained as follows: 1) processes of pigmentation formation can be divided into reversible and irreversible processes; 2) the reversible process is related to quantitative changes in pigments; 3) the irreversible process is related to the formation of stable pigmentation patterns that are physiologically inherited; 4) the spot pattern of color pigmentation represents an independent island of the physiological system; 5) many types of activators and inhibitors are involved in flower pigmentation production; 6) the patterns of flower pigmentation can be modulated; and 7) the evolution mechanism of organogenesis can be separated into several steps of independent thermodynamic processes. Our conclusion is that the thermodynamic system, rather than the dynamic system, is the essential and fundamental attribute of biological behaviors.

In silico characterization of the psilocybin biosynthesis pathway

Nearly all mushrooms of the Psilocybe genus contain the natural product psilocybin, which is a psychoactive alkaloid derived from l-tryptophan. Considering their use in ancient times, as well as their psychedelic properties, these mushrooms have re-emerged with psychotherapeutic potential for treating depression, which has triggered increased pharmaceutical interest. However, the psilocybin biosynthesis pathway was only recently defined and, as such, little exists in the way of structural data. Accordingly, the aim of this study was to structurally characterize this pathway by generating homology models for the four Psilocybe cubensis enzymes involved in psilocybin biosynthesis (PsiD, a decarboxylase; PsiH, a monooxygenase; PsiK, a phosphotransferase; PsiM, a methyltransferase). Following initial model generation and alignment with the identified structural templates, repeated refinement of the models was carried out using secondary structure prediction, geometry evaluation, energy minimization, and molecular dynamics simulations in water. The final models were then evaluated using molecular docking interactions with their substrates, i.e., psilocybin precursors (l-tryptophan, tryptamine, 4-hydroxytryptamine, and norbaeocystin/baeocystin), all of which generated feasible binding modes for the expected biotransformation. Further plausibility of the psilocybin → aeruginascin, 4-hydroxytryptamine → norpsilocin, and tryptamine → N,N-dimethyltryptamine conversions, all mediated by the generated model for PsiM, suggests valid routes of formation for these key secondary metabolites. The structural characterization of these enzymes and their binding modes which emerged from this study can lead to a better understanding of psilocybin synthesis, thereby paving the way for the development of novel substrates and selective inhibitors, as well as improved biotechnological manipulation and production of psilocybin in vitro.

Evaluation of Pathogenicity and Structural Alterations for the Mutations Identified in the Conserved Region of the C-Terminal Kinase Domain of Human-Ribosomal S6 Kinase 1

Human-ribosomal s6 kinase 1 (h-RSK1) is an effector kinase of the Ras/MAPK signaling pathway, which is involved in the regulation of the cell cycle, proliferation, and survival. RSKs comprise two functionally distinct kinase domains at the N-terminal (NTKD) and C-terminal (CTKD) separated by a linker region. The mutations in RSK1 may have the potential to provide an extra benefit to the cancer cell to proliferate, migrate, and survive. The present study focuses on evaluating the structural basis for the missense mutations identified at the C-terminal kinase domain of human-RSK1. A total of 139 mutations reported on RSK1 were retrieved from cBioPortal, where 62 were located at the CTKD region. Furthermore, 10 missense mutations Arg434Pro, Thr701Met, Ala704Thr, Arg725Trp, Arg726Gln, His533Asn, Pro613Leu, Ser720Cys, Arg725Gln, and Ser732Phe were predicted to be deleterious using in silico tools. To our observation, these mutations are located in the evolutionarily conserved region of RSK1 and shown to alter the inter- and intramolecular interactions and also the conformational stability of RSK1-CTKD. The molecular dynamics (MD) simulation study further revealed that the five mutations Arg434Pro, Thr701Met, Ala704Thr, Arg725Trp, and Arg726Gln showed maximum structural alterations in RSK1-CTKD. Thus, based on the in silico and MD simulation analysis, it can be concluded that the reported mutations may serve as potential candidates for further functional studies.

Decoding the Conformational Selective Mechanism of FGFR Isoforms: A Comparative Molecular Dynamics Simulation

Fibroblast growth factor receptors (FGFRs) play critical roles in the regulation of cell growth, differentiation, and proliferation. Specifically, FGFR2 gene amplification has been implicated in gastric and breast cancer. Pan-FGFR inhibitors often cause large toxic side effects, and the highly conserved ATP-binding pocket in the FGFR1/2/3 isoforms poses an immense challenge in designing selective FGFR2 inhibitors. Recently, an indazole-based inhibitor has been discovered that can selectively target FGFR2. However, the detailed mechanism involved in selective inhibition remains to be clarified. To this end, we performed extensive molecular dynamics simulations of the apo and inhibitor-bound systems along with multiple analyses, including Markov state models, principal component analysis, a cross-correlation matrix, binding free energy calculation, and community network analysis. Our results indicated that inhibitor binding induced the phosphate-binding loop (P-loop) of FGFR2 to switch from the open to the closed conformation. This effect enhanced extensive hydrophobic FGFR2-inhibitor contacts, contributing to inhibitor selectivity. Moreover, the key conformational intermediate states, dynamics, and driving forces of this transformation were uncovered. Overall, these findings not only provided a structural basis for understanding the closed P-loop conformation for therapeutic potential but also shed light on the design of selective inhibitors for treating specific types of cancer.

Strategy toward Kinase-Selective Drug Discovery

Kinase drug selectivity is the ground challenge in cancer research. Due to the structurally similar kinase drug pockets, off-target inhibitor toxicity has been a major cause for clinical trial failures. The pockets are similar but not identical. Here, we describe a transformation invariant protocol to identify distinct geometric features in the drug pocket that can distinguish one kinase from all others. We integrate available experimental structures with the artificial intelligence-based structural kinome, performing a kinome-wide structural bioinformatic analysis to establish the structural principles of kinase drug selectivity. We generate the structural landscape from the experimental kinase-ligand complexes and propose a binary network that encapsulates the information. The results show that all kinases contain binary units that are shared by less than seven other kinases in the kinome. 331 kinases contain unique binary units that may distinguish them from all others. The structural features encoded by these binary units in the network represent the inhibitor-accessible geometric space that may capture the kinome-wide selectivity. Our proposed binary network with the unsupervised clustering can serve as a general structural bioinformatic protocol for extracting the distinguishing structural features for any protein from their families. We apply the binary network to epidermal growth factor receptor tyrosine kinase inhibitor selectivity by targeting the gate area and the AKT1 serine/threonine kinase selectivity by binding to the αC-helix region and the allosteric pocket. Finally, we develop the cross-platform software, KDS (Kinase Drug Selectivity), for customized visualization and analysis of the binary networks in the human kinome (https://github.com/CBIIT/KDS).