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

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.

Conformational Strain Indicated by Ramachandran Angles for the Protein Backbone Is Only Weakly Related to the Flexibility

Studies on energy associated with free dipeptides have shown that conformers with unfavorable (ϕ,ψ) torsion angles have higher energy compared to conformers with favorable (ϕ,ψ) angles. It is expected that higher energy confers higher dynamics and flexibility to that part of the protein. Here, we explore a potential relationship between conformational strain in a residue due to unfavorable (ϕ,ψ) angles and its flexibility and dynamics in the context of protein structures. We compared flexibility of strained and relaxed residues, which are recognized based on outlier/allowed and favorable (ϕ,ψ) angles respectively, using normal-mode analysis (NMA). We also performed in-depth analysis on flexibility and dynamics at catalytic residues in protein kinases, which exhibit different strain status in different kinase structures using NMA and molecular dynamics simulations. We underline that strain of a residue, as defined by backbone torsion angles, is almost unrelated to the flexibility and dynamics associated with it. Even the overall trend observed among all high-resolution structures in which relaxed residues tend to have slightly higher flexibility than strained residues is counterintuitive. Consequently, we propose that identifying strained residues based on (ϕ,ψ) values is not an effective way to recognize energetic strain in protein structures.

Coevolutionary Landscape of Kinase Family Proteins: Sequence Probabilities and Functional Motifs

The protein kinase catalytic domain is one of the most abundant domains across all branches of life. Although kinases share a common core function of phosphoryl-transfer, they also have wide functional diversity and play varied roles in cell signaling networks, and for this reason are implicated in a number of human diseases. This functional diversity is primarily achieved through sequence variation, and uncovering the sequence-function relationships for the kinase family is a major challenge. In this study we use a statistical inference technique inspired by statistical physics, which builds a coevolutionary "Potts" Hamiltonian model of sequence variation in a protein family. We show how this model has sufficient power to predict the probability of specific subsequences in the highly diverged kinase family, which we verify by comparing the model's predictions with experimental observations in the Uniprot database. We show that the pairwise (residue-residue) interaction terms of the statistical model are necessary and sufficient to capture higher-than-pairwise mutation patterns of natural kinase sequences. We observe that previously identified functional sets of residues have much stronger correlated interaction scores than are typical.

Computational and Experimental Characterization of Patient Derived Mutations Reveal an Unusual Mode of Regulatory Spine Assembly and Drug Sensitivity in EGFR Kinase

The catalytic activation of protein kinases requires precise positioning of key conserved catalytic and regulatory motifs in the kinase core. The Regulatory Spine (RS) is one such structural motif that is dynamically assembled upon kinase activation. The RS is also a mutational hotspot in cancers; however, the mechanisms by which cancer mutations impact RS assembly and kinase activity are not fully understood. In this study, through mutational analysis of patient derived mutations in the RS of EGFR kinase, we identify an activating mutation, M766T, at the RS3 position. RS3 is located in the regulatory αC-helix, and a series of mutations at the RS3 position suggest a strong correlation between the amino acid type present at the RS3 position and ligand (EGF) independent EGFR activation. Small polar amino acids increase ligand independent activity, while large aromatic amino acids decrease kinase activity. M766T relies on the canonical asymmetric dimer for full activation. Molecular modeling and molecular dynamics simulations of WT and mutant EGFR suggest a model in which M766T activates the kinase domain by disrupting conserved autoinhibitory interactions between M766 and hydrophobic residues in the activation segment. In addition, a water mediated hydrogen bond network between T766, the conserved K745-E762 salt bridge, and the backbone amide of the DFG motif is identified as a key determinant of M766T-mediated activation. M766T is resistant to FDA approved EGFR inhibitors such as gefitinib and erlotinib, and computational estimation of ligand binding free energy identifies key residues associated with drug sensitivity. In sum, our studies suggest an unusual mode of RS assembly and oncogenic EGFR activation, and provide new clues for the design of allosteric protein kinase inhibitors.

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

Protein tyrosine kinases (PTKs) are a group of closely related enzymes that have evolutionarily diverged from serine/threonine kinases (STKs) to regulate pathways associated with multi-cellularity. Evolutionary divergence of PTKs from STKs has occurred through accumulation of mutations in the active site as well as in the commonly conserved hydrophobic core. While the functional significance of active site variations is well understood, relatively little is known about how hydrophobic core variations contribute to PTK evolutionary divergence. Here, using a combination of statistical sequence comparisons, molecular dynamics simulations, mutational analysis and in vitro thermostability and kinase assays, we investigate the structural and functional significance of key PTK-specific variations in the kinase core. We find that the nature of residues and interactions in the hydrophobic core of PTKs is strikingly different from other protein kinases, and PTK-specific variations in the core contribute to functional divergence by altering the stability and dynamics of the kinase domain. In particular, a functionally critical STK-conserved histidine that stabilizes the regulatory spine in STKs is selectively mutated to an alanine, serine or glutamate in PTKs, and this loss-of-function mutation is accommodated, in part, through compensatory PTK-specific interactions in the core. In particular, a PTK-conserved phenylalanine in the I-helix appears to structurally and functionally compensate for the loss of STK-histidine by interacting with the regulatory spine, which has far-reaching effects on enzyme activity, inhibitor sensing, and stability. We propose that hydrophobic core variations provide a selective advantage during PTK evolution by increasing the conformational flexibility, and therefore the allosteric potential of the kinase domain. Our studies also suggest that Tyrosine Kinase Like kinases such as RAF are intermediates in PTK evolutionary divergence inasmuch as they share features of both PTKs and STKs in the core. Finally, our studies provide an evolutionary framework for identifying and characterizing disease and drug resistance mutations in the kinase core.

Proteins evolve new functions through accumulation of mutations in the primary sequence. Understanding how naturally occurring mutations shape protein function can provide insights into how non-natural mutations contribute to disease. Here, we identify sequence variants associated with the functional specialization of tyrosine kinases, a large and medically important class of proteins associated with the evolution of complex multicellular functions and diseases such as cancer. We find that mutations distal from the active site contribute to functional specialization by altering the stability, activity and dynamics of the kinase core. Our findings have implications for understanding the evolution of allosteric regulation in tyrosine kinases, and in predicting the structural and functional impact of disease and drug resistance mutations in the kinase core.

Small-world networks of residue interactions in the Abl kinase complexes with cancer drugs: topology of allosteric communication pathways can determine drug resistance effects

Computational modelling of efficiency and robustness of the residue interaction networks and allosteric pathways in kinase structures can characterize protein kinase sensitivity to drug binding and drug resistance effects.

Structure-based network analysis of activation mechanisms in the ErbB family of receptor tyrosine kinases: the regulatory spine residues are global mediators of structural stability and allosteric interactions

The ErbB protein tyrosine kinases are among the most important cell signaling families and mutation-induced modulation of their activity is associated with diverse functions in biological networks and human disease. We have combined molecular dynamics simulations of the ErbB kinases with the protein structure network modeling to characterize the reorganization of the residue interaction networks during conformational equilibrium changes in the normal and oncogenic forms. Structural stability and network analyses have identified local communities integrated around high centrality sites that correspond to the regulatory spine residues. This analysis has provided a quantitative insight to the mechanism of mutation-induced “superacceptor” activity in oncogenic EGFR dimers. We have found that kinase activation may be determined by allosteric interactions between modules of structurally stable residues that synchronize the dynamics in the nucleotide binding site and the αC-helix with the collective motions of the integrating αF-helix and the substrate binding site. The results of this study have pointed to a central role of the conserved His-Arg-Asp (HRD) motif in the catalytic loop and the Asp-Phe-Gly (DFG) motif as key mediators of structural stability and allosteric communications in the ErbB kinases. We have determined that residues that are indispensable for kinase regulation and catalysis often corresponded to the high centrality nodes within the protein structure network and could be distinguished by their unique network signatures. The optimal communication pathways are also controlled by these nodes and may ensure efficient allosteric signaling in the functional kinase state. Structure-based network analysis has quantified subtle effects of ATP binding on conformational dynamics and stability of the EGFR structures. Consistent with the NMR studies, we have found that nucleotide-induced modulation of the residue interaction networks is not limited to the ATP site, and may enhance allosteric cooperativity with the substrate binding region by increasing communication capabilities of mediating residues.