Long-range molecular dynamics show that inactive forms of Protein Kinase A are more dynamic than active forms

Dynamic Interplay of Loop Motions Governs the Molecular Level Regulatory Dynamics in Spleen Tyrosine Kinase: Insights from Molecular Dynamics Simulations

The spleen tyrosine kinase (Syk) is a key regulator in immune cell signaling and is linked to various mechanisms in cancer and neurodegenerative diseases. Although most computational research on Syk focuses on novel drug design, the molecular-level regulatory dynamics remain unexplored. In this study, we utilized 5 × 1 μs all-atom molecular dynamics simulations of the Syk kinase domain, examining it in combinations of activation segment phosphorylated/unphosphorylated (at Tyr525, Tyr526) and the "DFG"-Asp protonated/deprotonated (at Asp512) states to investigate conformational variations and regulatory dynamics of various loops and motifs within the kinase domain. Our findings revealed that the formation and disruption of several electrostatic interactions among residues within and near the activation segment likely influenced its dynamics. The protein structure network analysis indicated that the N-terminal and C-terminal anchors were stabilized by connections with the nearby stable helical regions. The P-loop showed conformational variation characterized by movements toward and away from the conserved "HRD"-motif. Additionally, there was a significant correlation between the movement of the β3-αC loop and the P-loop, which controls the dimensions of the adenine-binding cavity of the C-spine region. Overall, understanding these significant motions of the Syk kinase domain enhances our knowledge of its functional regulatory mechanism and can guide future research.

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.

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

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

Allosteric regulation of the inactive to active state conformational transition in CDPK1 protein of Plasmodium falciparum

The aim of the current study is to investigate the role of the CAD domain in the activation mechanism of calcium dependent protein kinase-1 of Plasmodium falciparum (PfCDPK1) and explore the possibility of allosteric inhibition of this kinase. PfCDPK1 belongs to CDPK family of apicomplexan kinases which have a C-terminal CAD domain. Microsecond scale MD simulations were performed on modeled structures of complete PfCDPK1 and its kinase domain alone. The simulations revealed that in absence of CAD the salt bridge between Glu116 in αC-helix and Lys85 in β3-sheet of kinase breaks after 200 ns resulting in inactive conformation of the kinase, but the salt bridge stays intact in the complete protein stabilizing it in active conformation. These results highlight the novel CAD mediated allosteric stabilization of the crucial salt bridge which is a hallmark of active conformation of kinase domains. The mechanistic details of the allosteric activation revealed by our study, opens up the possibility for design of allosteric inhibitors of PfCDPK1 kinase by disrupting the kinase:CAD interactions. Using a combination of machine learning and structure-based in silico screening, we have identified novel PPI modulators for allosteric inactivation of PfCDPK1 kinase.

Phosphorylation-Induced Conformational Dynamics and Inhibition of Janus Kinase 1 by Suppressors of Cytokine Signaling 1

The dysfunction of the JAK/STAT (Janus kinase/signal transducers and activators of transcription) pathway results in several pathophysiological conditions, including autoimmune disorders. The negative feedback regulators of the JAK/STAT signaling pathway, suppressors of cytokine signaling (SOCS), act as a natural inhibitor of JAK and inhibit aberrant activity. SOCS1 is the most potent member of the SOCS family, whose kinase inhibitory region targets the substrate-binding groove of JAK with high affinity and blocks the phosphorylation of JAK kinases. Overall, we performed an aggregate of 13 μs molecular dynamics simulations on the activation loop's three different phosphorylation (double and single) states. Results from our simulations show that the single Tyr¹⁰³⁴ phosphorylation could stabilize the JAK1/SOCS1 complex as well as the flexible activation segment. The phosphate-binding loop (P-loop) shows conformational variability at dual and single phosphorylated states. Principal component analysis and protein structure network (PSN) analysis reveal that the different phosphorylation states and SOCS1 binding induce intermediate inactive conformations of JAK1, which could be a better target for future JAK1 selective drug design. PSN analysis suggests that the com-pY1034 system is stabilized due to higher values of network hubs than the other two complex systems. Moreover, the binding free energy calculations suggest that pTyr¹⁰³⁴ states show a higher affinity toward SOCS1 than the dual and pTyr¹⁰³⁵ states. We believe that the mechanistic understanding of JAK1/SOCS1 complexation will aid future studies related to peptide inhibitors based on SOCS1.

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.

Hypervariability of accessible and inaccessible conformational space of proteins

Proteins undergo motions in a range of amplitudes, from domain motions to backbone rotations, leading to changes in (φ,ψ) torsion angles and small-scale bond vibrations and angle bending. Here, we study the extent of variations in (φ,ψ) values in proteins and the effects of bond geometry variations due to vibrational motions in a protein on the accessible, (steric clash-free) (φ,ψ) space. We perform 1-fs timestep unconstrained molecular dynamics simulations on super-high-resolution protein structures. Extent of variations in bond geometry during the simulation is within acceptable ranges of bond lengths and angles. However, the steric clash-free (φ,ψ) space continuously changes as seen in bond geometry-specific (φ,ψ) steric maps at the residue level during simulations. (φ,ψ) regions that have steric clash at one timepoint can become steric clash-free at a different timepoint through minor adjustments to backbone bond lengths and angles. Also instances of (φ,ψ) transitions from the left to right half of the (φ,ψ) map in consecutive snapshots of the trajectory are seen. Although the two quadrants are separated by a steric clash-prone region, corresponding to a high-energy barrier, height of this barrier is lowered by adjusting the bond geometry such that a bridging region of steric clash-free, low-energy (φ,ψ) values is formed. We demonstrate the idea of dynamically varying nature of acceptable and accessible (φ,ψ) steric space in proteins, which has implications for protein folding; proteins could sample (φ,ψ) space which is originally considered to be inaccessible, during folding, through minor adjustments to their backbone bond geometry.

Characterizing an allosteric inhibitor-induced inactive state in with-no-lysine kinase 1 using Gaussian accelerated molecular dynamics simulations

The binding of an allosteric inhibitor in WNK1 leads to the inactive state.

Activation of RSK by phosphomimetic substitution in the activation loop is prevented by structural constraints

The activation of the majority of AGC kinases is regulated by two phosphorylation events on two conserved serine/threonine residues located on the activation loop and on the hydrophobic motif, respectively. In AGC kinase family, phosphomimetic substitutions with aspartate or glutamate, leading to constitutive activation, have frequently occurred at the hydrophobic motif site. On the contrary, phosphomimetic substitutions in the activation loop are absent across the evolution of AGC kinases. This observation is explained by the failure of aspartate and glutamate to mimic phosphorylatable serine/threonine in this regulatory site. By detailed 3D structural simulations of RSK2 and further biochemical evaluation in cells, we show that the phosphomimetic residue on the activation loop fails to form a critical salt bridge with R114, necessary to reorient the αC-helix and to activate the protein. By a phylogenetic analysis, we point at a possible coevolution of a phosphorylatable activation loop and the presence of a conserved positively charged amino acid on the αC-helix. In sum, our analysis leads to the unfeasibility of phosphomimetic substitution in the activation loop of RSK and, at the same time, highlights the peculiar structural role of activation loop phosphorylation.

ATP Controls the Aggregation of Aβ16-22 Peptides

The oligomerization of Aβ_16-22 peptide, which is the hydrophobic core region of full-length Aβ_1-42, causes Alzheimer's disease (AD). This progressive neurodegenerative disease affects over 44 million people worldwide. However, very few synthesized drug molecules are available to inhibit the aggregation of Aβ. Recently, experimental studies have shown that the biological ATP molecule prevents Aβ fibrillation at the millimolar scale; however, the significance of ATP molecules on Aβ fibrillation and the mechanism behind it remain elusive. We have carried out a total of 7.5 μs extensive all-atom molecular dynamics and 8.82 μs of umbrella sampling in explicit water using AMBER14SB, AMBER99SB-ILDN, and AMBER-FB15 force fields for Aβ_16-22 peptide, to investigate the role of ATP on the disruption of Aβ_16-22 prefibrils. From various analyses, such as secondary structure analysis, residue-wise contact map, SASA, and interaction energies, we have observed that, in the presence of ATP, the aggregation of Aβ_16-22 peptide is very unfavorable. Moreover, the biological molecule ATP interacts with the Aβ_16-22 peptide via hydrogen bonding, π-π stacking, and NH-π interactions which, ultimately, prevent the aggregation of Aβ_16-22 peptide. Hence, we assume that the deficiency of ATP may cause Alzheimer's disease (AD).