Comprehensive deep mutational scanning of three Src family kinases using phospho-PCA coupled with hierarchical Bayesian modeling reveals that kinase domain variants impact substrate specificity primarily through long-range allosteric networks rather than direct binding interactions, with most specificity-determining residues located distal to the substrate interface and operating through diverse mechanisms including electrostatic remodeling, conformational biasing, and metastable state stabilization.

Protein kinases relay information to various cellular processes, and their dysregulation underlies numerous human diseases. Despite their importance, our understanding of how kinase domains and their variants impact protein stability, catalytic activity, and substrate recognition is incomplete. In this work, we develop the phosphorylation protein complementation assay (phospho-PCA), which enables quantitative measurements of kinase-substrate interactions by coupling them to the growth of bud-ding yeast, thereby enabling deep mutational scanning of kinase domains. When combined with deep mutational scans targeting folding stability and Bayesian modeling, phospho-PCA can disentangle the relative impact of mutations on kinase domain stability, catalytic activity, and substrate specificity. We demonstrate the accuracy and breadth of phospho-PCA, showing its applicability to both tyrosine and serine/threonine kinase domains. We then apply our method to three closely related protein kinases with distinct substrate preferences, evaluating over 15,000 kinase variants against a panel of three substrates for both catalytic activity and substrate specificity. The resulting dataset constitutes the largest and most detailed variant-to-function map assembled for this enzyme family to date, revealing numerous mutations that alter kinase activity and substrate specificity. Physics-based modeling reveals how these mutations operate through diverse mechanisms, including long-range allosteric communication, to alter both activity and substrate specificity. Given its scalability, we believe phospho-PCA can measure the functional impact of variants across the entire human kinome.