Polymodal allosteric regulation of Type 1 Serine/Threonine Kinase Receptors via a conserved electrostatic lock

Type 1 Serine/Threonine Kinase Receptors (STKR1) transduce a wide spectrum of biological signals mediated by TGF-β superfamily members. The STKR1 activity is tightly controlled by their regulatory glycine-serine rich (GS) domain adjacent to the kinase domain. Despite decades of studies, it remains unknown how physiological or pathological GS domain modifications are coupled to STKR1 kinase activity. Here, by performing molecular dynamics simulations and free energy calculation of Activin-Like Kinase 2 (ALK2), we found that GS domain phosphorylation, FKBP12 dissociation, and disease mutations all destabilize a D354-R375 salt-bridge, which normally acts as an electrostatic lock to prevent coordination of adenosine triphosphate (ATP) to the catalytic site. We developed a WAFEX-guided principal analysis and unraveled how phosphorylation destabilizes this highly conserved salt-bridge in temporal and physical space. Using current-flow betweenness scores, we identified an allosteric network of residue-residue contacts between the GS domain and the catalytic site that controls the formation and disruption of this salt bridge. Importantly, our novel network analysis approach revealed how certain disease-causing mutations bypass FKBP12-mediated kinase inhibition to produce leaky signaling. We further provide experimental evidence that this salt-bridge lock exists in other STKR1s, and acts as a general safety mechanism in STKR1 to prevent pathological leaky signaling. In summary, our study provides a compelling and unifying allosteric activation mechanism in STKR1 kinases that reconciles a large number of experimental studies and sheds light on a novel therapeutic avenue to target disease-related STKR1 mutants.

[1]  Stefan Knapp,et al.  Structure of the Bone Morphogenetic Protein Receptor ALK2 and Implications for Fibrodysplasia Ossificans Progressiva* , 2012, The Journal of Biological Chemistry.

[2]  Duncan Poole,et al.  Routine Microsecond Molecular Dynamics Simulations with AMBER on GPUs. 1. Generalized Born , 2012, Journal of chemical theory and computation.

[3]  Chris Jones,et al.  ACVR1 mutations in DIPG: lessons learned from FOP. , 2014, Cancer research.

[4]  Jan H. Jensen,et al.  Improved Treatment of Ligands and Coupling Effects in Empirical Calculation and Rationalization of pKa Values. , 2011, Journal of chemical theory and computation.

[5]  J. Massagué,et al.  Mechanisms of TGF-β Signaling from Cell Membrane to the Nucleus , 2003, Cell.

[6]  Ivet Bahar,et al.  ProDy: Protein Dynamics Inferred from Theory and Experiments , 2011, Bioinform..

[7]  K Schulten,et al.  VMD: visual molecular dynamics. , 1996, Journal of molecular graphics.

[8]  Robert I Cukier,et al.  Molecular dynamics of apo-adenylate kinase: a principal component analysis. , 2006, The journal of physical chemistry. B.

[9]  Enrico Bozzo,et al.  Effective and efficient approximations of the generalized inverse of the graph Laplacian matrix with an application to current-flow betweenness centrality , 2012, ArXiv.

[10]  Laxmikant V. Kalé,et al.  Scalable molecular dynamics with NAMD , 2005, J. Comput. Chem..

[11]  Adrian E Roitberg,et al.  Optimization of Umbrella Sampling Replica Exchange Molecular Dynamics by Replica Positioning. , 2013, Journal of chemical theory and computation.

[12]  J. Massagué TGFβ signalling in context , 2012, Nature Reviews Molecular Cell Biology.

[13]  Susan S. Taylor,et al.  Surface comparison of active and inactive protein kinases identifies a conserved activation mechanism , 2006, Proceedings of the National Academy of Sciences.

[14]  In Ho Choi,et al.  A recurrent mutation in the BMP type I receptor ACVR1 causes inherited and sporadic fibrodysplasia ossificans progressiva , 2006, Nature Genetics.

[15]  D. Case,et al.  Constant pH molecular dynamics in generalized Born implicit solvent , 2004, J. Comput. Chem..

[16]  Mark E. J. Newman A measure of betweenness centrality based on random walks , 2005, Soc. Networks.

[17]  Zaida Luthey-Schulten,et al.  NetworkView: 3D display and analysis of protein·RNA interaction networks , 2012, Bioinform..

[18]  L. Mathews,et al.  Characterization of type II activin receptors. Binding, processing, and phosphorylation. , 1993, The Journal of biological chemistry.

[19]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[20]  Taehoon Kim,et al.  CHARMM‐GUI: A web‐based graphical user interface for CHARMM , 2008, J. Comput. Chem..

[21]  Anja Nohe,et al.  The Mode of Bone Morphogenetic Protein (BMP) Receptor Oligomerization Determines Different BMP-2 Signaling Pathways* , 2002, The Journal of Biological Chemistry.

[22]  Jan H. Jensen,et al.  PROPKA3: Consistent Treatment of Internal and Surface Residues in Empirical pKa Predictions. , 2011, Journal of chemical theory and computation.

[23]  Eileen M. Shore,et al.  In vitro Analyses of the Dysregulated R206H ALK2 Kinase-FKBP12 Interaction Associated with Heterotopic Ossification in FOP , 2011, Cells Tissues Organs.

[24]  Gábor Csárdi,et al.  The igraph software package for complex network research , 2006 .

[25]  Y. Sugita,et al.  Multidimensional replica-exchange method for free-energy calculations , 2000, cond-mat/0009120.

[26]  Daniel R Roe,et al.  PTRAJ and CPPTRAJ: Software for Processing and Analysis of Molecular Dynamics Trajectory Data. , 2013, Journal of chemical theory and computation.

[27]  J. Massagué,et al.  Mechanisms of TGF-beta signaling from cell membrane to the nucleus. , 2003, Cell.

[28]  Alexander D. MacKerell,et al.  Optimization of the additive CHARMM all-atom protein force field targeting improved sampling of the backbone φ, ψ and side-chain χ(1) and χ(2) dihedral angles. , 2012, Journal of chemical theory and computation.

[29]  Stephen Yip,et al.  Recurrent activating ACVR1 mutations in diffuse intrinsic pontine glioma , 2014, Nature Genetics.

[30]  J. Kuriyan,et al.  The Conformational Plasticity of Protein Kinases , 2002, Cell.

[31]  Jeffrey L. Wrana,et al.  Mechanism of activation of the TGF-β receptor , 1994, Nature.

[32]  Daniel R. Roe,et al.  Using Wavelet Analysis To Assist in Identification of Significant Events in Molecular Dynamics Simulations , 2016, J. Chem. Inf. Model..

[33]  J Kuriyan,et al.  The TGF beta receptor activation process: an inhibitor- to substrate-binding switch. , 2001, Molecular cell.

[34]  Morgan Huse,et al.  Crystal Structure of the Cytoplasmic Domain of the Type I TGF β Receptor in Complex with FKBP12 , 1999, Cell.

[35]  K. Miyazono,et al.  Signal transduction via serine/threonine kinase receptors. , 1994, Seminars in cell biology.

[36]  Enrico Bozzo,et al.  Approximations of the Generalized Inverse of the Graph Laplacian Matrix , 2012, Internet Math..

[37]  P. Donahoe,et al.  Specific interaction of type I receptors of the TGF-beta family with the immunophilin FKBP-12. , 1994, Science.

[38]  R Wieser,et al.  Mechanism of activation of the TGF-beta receptor. , 1994, Nature.

[39]  K. Miyazono,et al.  Characterization of the Interaction of FKBP12 with the Transforming Growth Factor-β Type I Receptor in Vivo* , 1996, The Journal of Biological Chemistry.

[40]  Benoît Roux,et al.  Calculation of Free Energy Landscape in Multi-Dimensions with Hamiltonian-Exchange Umbrella Sampling on Petascale Supercomputer. , 2012, Journal of chemical theory and computation.

[41]  Lei Huang,et al.  Generalized scalable multiple copy algorithms for molecular dynamics simulations in NAMD , 2014, Comput. Phys. Commun..

[42]  Duncan Poole,et al.  Routine Microsecond Molecular Dynamics Simulations with AMBER on GPUs. 2. Explicit Solvent Particle Mesh Ewald. , 2013, Journal of chemical theory and computation.

[43]  Amar Gajjar,et al.  The genomic landscape of diffuse intrinsic pontine glioma and pediatric non-brainstem high-grade glioma , 2014, Nature Genetics.

[44]  Kohei Miyazono,et al.  Bone morphogenetic protein receptors and signal transduction. , 2010, Journal of biochemistry.