BlockMaster: partitioning protein kinase structures using normal-mode analysis

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.

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.

Molecular Dynamics and Machine Learning Give Insights on the Flexibility-Activity Relationships in Tyrosine Kinome

Tyrosine kinases are a subfamily of kinases with critical roles in cellular machinery. Dysregulation of their active or inactive forms is associated with diseases like cancer. This study aimed to holistically understand their flexibility-activity relationships, focusing on pockets and fluctuations. We studied 43 different tyrosine kinases by collecting 120 μs of molecular dynamics simulations, pocket and residue fluctuation analysis, and a complementary machine learning approach. We found that the inactive forms often have increased flexibility, particularly at the DFG motif level. Noteworthy, thanks to these long simulations combined with a decision tree, we identified a semiquantitative fluctuation threshold of the DGF+3 residue over which the kinase has a higher probability to be in the inactive form.

Protein kinases in cardiovascular diseases

ABSTRACT

Cardiovascular disease (CVD) remains the leading cause of death worldwide. Therefore, exploring the mechanism of CVDs and critical regulatory factors is of great significance for promoting heart repair, reversing cardiac remodeling, and reducing adverse cardiovascular events. Recently, significant progress has been made in understanding the function of protein kinases and their interactions with other regulatory proteins in myocardial biology. Protein kinases are positioned as critical regulators at the intersection of multiple signals and coordinate nearly every aspect of myocardial responses, regulating contractility, metabolism, transcription, and cellular death. Equally, reconstructing the disrupted protein kinases regulatory network will help reverse pathological progress and stimulate cardiac repair. This review summarizes recent researches concerning the function of protein kinases in CVDs, discusses their promising clinical applications, and explores potential targets for future treatments.

Pseudokinases repurpose flexibility signatures associated with the protein kinase fold for noncatalytic roles

The bilobal protein kinase-like fold in pseudokinases lack one or more catalytic residues, conserved in canonical protein kinases, and are considered enzymatically deficient. Tertiary structures of pseudokinases reveal that their loops topologically equivalent to activation segments of kinases adopt contracted configurations, which is typically extended in active conformation of kinases. Herein, anisotropic network model based normal mode analysis (NMA) was conducted on 51 active conformation structures of protein kinases and 26 crystal structures of pseudokinases. Our observations indicate that although backbone fluctuation profiles are similar for individual kinase-pseudokinase families, low intensity mean square fluctuations in pseudo-activation segment and other sub-structures impart rigidity to pseudokinases. Analyses of collective motions from functional modes reveal that pseudokinases, compared to active kinases, undergo distinct conformational transitions using the same structural fold. All-atom NMA of protein kinase-pseudokinase pairs from each family, sharing high amino acid sequence identities, yielded distinct community clusters, partitioned by residues exhibiting highly correlated fluctuations. It appears that atomic fluctuations from equivalent activation segments guide community membership and network topologies for respective kinase and pseudokinase. Our findings indicate that such adaptations in backbone and side-chain fluctuations render pseudokinases competent for catalysis-independent roles.

Extracting Dynamical Correlations and Identifying Key Residues for Allosteric Communication in Proteins by correlationplus

Extracting dynamical pairwise correlations and identifying key residues from large molecular dynamics trajectories or normal-mode analysis of coarse-grained models are important for explaining various processes like ligand binding, mutational effects, and long-distance interactions. Efficient and flexible tools to perform this task can provide new insights about residues involved in allosteric regulation and protein function. In addition, combining and comparing dynamical coupling information with sequence coevolution data can help to understand better protein function. To this aim, we developed a Python package called correlationplus to calculate, visualize, and analyze pairwise correlations. In this way, the package aids to identify key residues and interactions in proteins. The source code of correlationplus is available under LGPL version 3 at https://github.com/tekpinar/correlationplus. The current version of the package (0.2.0) can be installed with common installation methods like conda or pip in addition to source code installation. Moreover, docker images are also available for usage of the code without installation.

Tuning the "violin" of protein kinases: The role of dynamics-based allostery

The intricacies of allosteric regulation of protein kinases continue to engage the research community. Allostery, or control from a distance, is seen as a fundamental biomolecular mechanism for proteins. From the traditional methods of conformational selection and induced fit, the field has grown to include the role of protein motions in defining a dynamics-based allosteric approach. Harnessing of these continuous motions in the protein to exert allosteric effects can be defined by a "violin" model that focuses on distributions of protein vibrations as opposed to concerted pathways. According to this model, binding of an allosteric modifier causes global redistribution of dynamics in the protein kinase domain that leads to changes in its catalytic properties. This model is consistent with the "entropy-driven allostery" mechanism proposed by Cooper and Dryden in 1984 and does not require, but does not exclude, any major structural changes. We provide an overview of practical implementation of the violin model and how it stands amidst the other known models of protein allostery. Protein kinases have been described as the biomolecules of interest. © 2019 IUBMB Life, 71(6):685-696, 2019.

Protein dynamic communities from elastic network models align closely to the communities defined by molecular dynamics

Dynamic communities in proteins comprise the cohesive structural units that individually exhibit rigid body motions. These can correspond to structural domains, but are usually smaller parts that move with respect to one another in a protein's internal motions, key to its functional dynamics. Previous studies emphasized their importance to understand the nature of ligand-induced allosteric regulation. These studies reported that mutations to key community residues can hinder transmission of allosteric signals among the communities. Usually molecular dynamic (MD) simulations (~ 100 ns or longer) have been used to identify the communities-a demanding task for larger proteins. In the present study, we propose that dynamic communities obtained from MD simulations can also be obtained alternatively with simpler models-the elastic network models (ENMs). To verify this premise, we compare the specific communities obtained from MD and ENMs for 44 proteins. We evaluate the correspondence in communities from the two methods and compute the extent of agreement in the dynamic cross-correlation data used for community detection. Our study reveals a strong correspondence between the communities from MD and ENM and also good agreement for the residue cross-correlations. Importantly, we observe that the dynamic communities from MD can be closely reproduced with ENMs. With ENMs, we also compare the community structures of stable and unstable mutant forms of T4 Lysozyme with its wild-type. We find that communities for unstable mutants show substantially poorer agreement with the wild-type communities than do stable mutants, suggesting such ENM-based community structures can serve as a means to rapidly identify deleterious mutants.

Deciphering conformational transitions of proteins by small angle X-ray scattering and normal mode analysis

Structural flexibility and conformational rearrangements are often related to important functions of biological macromolecules, but the experimental characterization of such transitions with high-resolution techniques is challenging. At a lower resolution, small angle X-ray scattering (SAXS) can be used to obtain information on biomolecular shapes and transitions in solution. Here, we present SREFLEX, a hybrid modeling approach that uses normal mode analysis (NMA) to explore the conformational space of high-resolution models and refine the structure guided by the agreement with the experimental SAXS data. The method starts from a given conformation of the protein (which does not agree with the SAXS data). The structure is partitioned into pseudo-domains either using structural classification databases or automatically from the protein dynamics as predicted by the NMA. The algorithm proceeds hierarchically employing NMA to first probe large rearrangements and progresses into smaller and more localized movements. At the large rearrangements stage the pseudo-domains stay as rigid bodies allowing one to avoid structural disruptions inherent to the earlier NMA-based algorithms. To validate the approach, we compiled a representative benchmark set of 88 conformational states known experimentally at high resolution. The performance of the algorithm is demonstrated in the simulated data on the benchmark set and also in a number of experimental examples. SREFLEX is included into the ATSAS program package freely available to the academic users, both for download and in the on-line mode.

Dynamics-Driven Allostery in Protein Kinases

Protein kinases have very dynamic structures and their functionality strongly depends on their dynamic state. Active kinases reveal a dynamic pattern with residues clustering into semirigid communities that move in μs-ms timescale. Previously detected hydrophobic spines serve as connectors between communities. Communities do not follow the traditional subdomain structure of the kinase core or its secondary structure elements. Instead they are organized around main functional units. Integration of the communities depends on the assembly of the hydrophobic spine and phosphorylation of the activation loop. Single mutations can significantly disrupt the dynamic infrastructure and thereby interfere with long-distance allosteric signaling that propagates throughout the whole molecule. Dynamics is proposed to be the underlying mechanism for allosteric regulation in protein kinases.

Dynamic architecture of a protein kinase

Protein kinases are dynamically regulated signaling proteins that act as switches in the cell by phosphorylating target proteins. To establish a framework for analyzing linkages between structure, function, dynamics, and allostery in protein kinases, we carried out multiple microsecond-scale molecular-dynamics simulations of protein kinase A (PKA), an exemplar active kinase. We identified residue-residue correlated motions based on the concept of mutual information and used the Girvan-Newman method to partition PKA into structurally contiguous "communities." Most of these communities included 40-60 residues and were associated with a particular protein kinase function or a regulatory mechanism, and well-known motifs based on sequence and secondary structure were often split into different communities. The observed community maps were sensitive to the presence of different ligands and provide a new framework for interpreting long-distance allosteric coupling. Communication between different communities was also in agreement with the previously defined architecture of the protein kinase core based on the "hydrophobic spine" network. This finding gives us confidence in suggesting that community analyses can be used for other protein kinases and will provide an efficient tool for structural biologists. The communities also allow us to think about allosteric consequences of mutations that are linked to disease.

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

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

Deciphering the Arginine-binding preferences at the substrate-binding groove of Ser/Thr kinases by computational surface mapping

Protein kinases are key signaling enzymes that catalyze the transfer of γ-phosphate from an ATP molecule to a phospho-accepting residue in the substrate. Unraveling the molecular features that govern the preference of kinases for particular residues flanking the phosphoacceptor is important for understanding kinase specificities toward their substrates and for designing substrate-like peptidic inhibitors. We applied ANCHORSmap, a new fragment-based computational approach for mapping amino acid side chains on protein surfaces, to predict and characterize the preference of kinases toward Arginine binding. We focus on positions P−2 and P−5, commonly occupied by Arginine (Arg) in substrates of basophilic Ser/Thr kinases. The method accurately identified all the P−2/P−5 Arg binding sites previously determined by X-ray crystallography and produced Arg preferences that corresponded to those experimentally found by peptide arrays. The predicted Arg-binding positions and their associated pockets were analyzed in terms of shape, physicochemical properties, amino acid composition, and in-silico mutagenesis, providing structural rationalization for previously unexplained trends in kinase preferences toward Arg moieties. This methodology sheds light on several kinases that were described in the literature as having non-trivial preferences for Arg, and provides some surprising departures from the prevailing views regarding residues that determine kinase specificity toward Arg. In particular, we found that the preference for a P−5 Arg is not necessarily governed by the 170/230 acidic pair, as was previously assumed, but by several different pairs of acidic residues, selected from positions 133, 169, and 230 (PKA numbering). The acidic residue at position 230 serves as a pivotal element in recognizing Arg from both the P−2 and P−5 positions.

Protein kinases are key signaling enzymes and major drug targets that catalyze the transfer of phosphate group to a phospho-accepting residue in the substrate. Unraveling molecular features that govern the preference of kinases for particular residues flanking the phosphoacceptor (substrate consensus sequence, SCS) is important for understanding kinase-substrates specificities and for designing peptidic inhibitors. Current methods used to predict this set of essential residues usually rely on linking between experimentally determined SCSs to kinase sequences. As such, these methods are less sensitive when specificity is dictated by subtle or kinase-unique sequence/structural features. In this study, we took a different approach for studying kinases specificities, by applying a new fragment-based method for mapping amino acid side chains on protein surfaces. We predicted and characterized the preference of Ser/Thr kinases toward Arginine binding, using the unbound kinase structures. The method produced high quality predictions and was able to provide novel insights and interesting departures from the prevailing views regarding the specificity-determining elements governing specificity toward Arginine. This work paves the way for studying the kinase binding preferences for other amino acids, for predicting protein-peptide structures, for facilitating the design of novel inhibitors, and for re-engineering of kinase specificities.

Protein kinases: evolution of dynamic regulatory proteins

Eukayotic protein kinases evolved as a family of highly dynamic molecules with strictly organized internal architecture. A single hydrophobic F-helix serves as a central scaffold for assembly of the entire molecule. Two non-consecutive hydrophobic structures termed "spines" anchor all the elements important for catalysis to the F-helix. They make firm, but flexible, connections within the molecule, providing a high level of internal dynamics of the protein kinase. During the course of evolution, protein kinases developed a universal regulatory mechanism associated with a large activation segment that can be dynamically folded and unfolded in the course of cell functioning. Protein kinases thus represent a unique, highly dynamic, and precisely regulated set of switches that control most biological events in eukaryotic cells.