Cited by 734Advanced Proteomics Techniques and ApplicationsProtein Kinase Regulation and GTPase SignalingProtein Tyrosine Phosphatases

The DOI matches its Crossref record; 120 of the 120 references it cites were found in a registry or an index.

An atlas of substrate specificities for the human serine/threonine kinome

Open access

Published version of a preprint

Abstract

Protein phosphorylation is one of the most widespread post-translational modifications in biology 1,2 . With advances in mass-spectrometry-based phosphoproteomics, 90,000 sites of serine and threonine phosphorylation have so far been identified, and several thousand have been associated with human diseases and biological processes 3,4 . For the vast majority of phosphorylation events, it is not yet known which of the more than 300 protein serine/threonine (Ser/Thr) kinases encoded in the human genome are responsible 3 . Here we used synthetic peptide libraries to profile the substrate sequence specificity of 303 Ser/Thr kinases, comprising more than 84% of those predicted to be active in humans. Viewed in its entirety, the substrate specificity of the kinome was substantially more diverse than expected and was driven extensively by negative selectivity. We used our kinome-wide dataset to computationally annotate and identify the kinases capable of phosphorylating every reported phosphorylation site in the human Ser/Thr phosphoproteome. For the small minority of phosphosites for which the putative protein kinases involved have been previously reported, our predictions were in excellent agreement. When this approach was applied to examine the signalling response of tissues and cell lines to hormones, growth factors, targeted inhibitors and environmental or genetic perturbations, it revealed unexpected insights into pathway complexity and compensation. Overall, these studies reveal the intrinsic substrate specificity of the human Ser/Thr kinome, illuminate cellular signalling responses and provide a resource to link phosphorylation events to biological pathways.

Overview

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Integrity signals

Registry

Crossref · 10.1038/s41586-022-05575-3

Cites

120 of 120 references verified0 citation errors · 0 not checkable

Statistics

Figures

Fig. 1 Profiling the substrate specificity of the human serine/threonine kinome. a, Experimental workflow for the PSPA analysis and representative results. The schematic was created using BioRender. Z denotes fixed positions containing one of the 20 natural amino acids, or either phosphorylated Thr (pThr) or phosphorylated Tyr (pTyr). X denotes unfixed positions containing randomized mixtures of all natural amino acids except Ser, Thr and Cys. Darker spots indicate preferred residues. b, Dendrogram of the human protein kinome, highlighting the Ser/Thr kinases analysed in this study.
Fig. 1 Profiling the substrate specificity of the human serine/threonine kinome. a, Experimental workflow for the PSPA analysis and representative results. The schematic was created using BioRender. Z denotes fixed positions containing one of the 20 natural amino acids, or either phosphorylated Thr (pThr) or phosphorylated Tyr (pTyr). X denotes unfixed positions containing randomized mixtures of all natural amino acids except Ser, Thr and Cys. Darker spots indicate preferred residues. b, Dendrogram of the human protein kinome, highlighting the Ser/Thr kinases analysed in this study.CC BY · PubMed Central
Fig. 2 Phosphorylation-site motif tree of the human Ser/Thr kinome. Hierarchical clustering of 303 Ser/Thr kinases on the basis of their amino acid motif selectivity (PSSMs). Kinase names are colour labelled according to their phylogenetic relationships (top right)2.
Fig. 2 Phosphorylation-site motif tree of the human Ser/Thr kinome. Hierarchical clustering of 303 Ser/Thr kinases on the basis of their amino acid motif selectivity (PSSMs). Kinase names are colour labelled according to their phylogenetic relationships (top right)2.CC BY · PubMed Central
Fig. 3 Phosphorylation motifs for the human Ser/Thr kinome enable comprehensive scoring and annotation of the human phosphoproteome. a, Schematic of the substrate-scoring process4. b, Results for Ser15 on glycogen phosphorylase alongside PSSM and the substrate motif logo of its established kinase glycogen phosphorylase kinase. c, The results for Ser15 of p53 alongside its established kinase ATM. d, Annotation of the human Ser and Thr phosphoproteome by percentile scores from 303 Ser/Thr kinases performed as shown in a. A total of 89,752 phosphorylation sites that were identified using high-throughput approaches4 and/or low-throughput approaches3 were sorted along the x-axis by their numbers of kinases with percentile scores higher than 90. On the y-axis, kinase percentile scores were sorted by rank separately for each site and represented in the heat map. Examples of well-studied kinase–substrate relationships are highlighted (yellow squares). Inset: phosphorylation sites on the left end of the plot scored favourably for many kinases, whereas sites on the right end scored favourably for fewer kinases.
Fig. 3 Phosphorylation motifs for the human Ser/Thr kinome enable comprehensive scoring and annotation of the human phosphoproteome. a, Schematic of the substrate-scoring process4. b, Results for Ser15 on glycogen phosphorylase alongside PSSM and the substrate motif logo of its established kinase glycogen phosphorylase kinase. c, The results for Ser15 of p53 alongside its established kinase ATM. d, Annotation of the human Ser and Thr phosphoproteome by percentile scores from 303 Ser/Thr kinases performed as shown in a. A total of 89,752 phosphorylation sites that were identified using high-throughput approaches4 and/or low-throughput approaches3 were sorted along the x-axis by their numbers of kinases with percentile scores higher than 90. On the y-axis, kinase percentile scores were sorted by rank separately for each site and represented in the heat map. Examples of well-studied kinase–substrate relationships are highlighted (yellow squares). Inset: phosphorylation sites on the left end of the plot scored favourably for many kinases, whereas sites on the right end scored favourably for fewer kinases.CC BY · PubMed Central
Fig. 4 Global motif analysis reveals how kinase perturbations and pathway rewiring reshape the phosphoproteome. a, Workflow of the motif enrichment analysis of phosphoproteomics data. The schematic was created using BioRender. b–g, Results from published datasets. b, Conditioned medium of HepG2 cells after genetic deletion of FAM20C44. c, Cultured myotubes after 30 min treatment with 2 μM isoproterenol45. d, HeLa cells after mitotic arrest by treatment for 45 min with 0.1 μM PLK1 inhibitor BI 2536 (ref. 46). e, A549 cells 2 h after exposure to 6 Gy of ionizing radiation49. f, 3T3-L1 adipocytes after serum starvation and then 1 min and 60 min treatment with 100 nM insulin54. g, C57BL/6J mouse bone-marrow-derived dendritic cells after 30 min and 4 h treatment with 100 ng ml−1 lipopolysaccharide (LPS)55. The enrichments in b–g were determined using one-sided exact Fisher’s tests and corrected for multiple hypotheses using the Benjamini–Hochberg method. Fully annotated versions of these plots are presented in Supplementary Fig. 2.
Fig. 4 Global motif analysis reveals how kinase perturbations and pathway rewiring reshape the phosphoproteome. a, Workflow of the motif enrichment analysis of phosphoproteomics data. The schematic was created using BioRender. b–g, Results from published datasets. b, Conditioned medium of HepG2 cells after genetic deletion of FAM20C44. c, Cultured myotubes after 30 min treatment with 2 μM isoproterenol45. d, HeLa cells after mitotic arrest by treatment for 45 min with 0.1 μM PLK1 inhibitor BI 2536 (ref. 46). e, A549 cells 2 h after exposure to 6 Gy of ionizing radiation49. f, 3T3-L1 adipocytes after serum starvation and then 1 min and 60 min treatment with 100 nM insulin54. g, C57BL/6J mouse bone-marrow-derived dendritic cells after 30 min and 4 h treatment with 100 ng ml−1 lipopolysaccharide (LPS)55. The enrichments in b–g were determined using one-sided exact Fisher’s tests and corrected for multiple hypotheses using the Benjamini–Hochberg method. Fully annotated versions of these plots are presented in Supplementary Fig. 2.CC BY · PubMed Central
Extended Data Fig. 1 Representation of phosphorylation site motifs in the human serine and threonine phosphoproteome. Venn diagram representation of the percentages of three prominent Ser/Thr kinase motif features, pertaining to Clusters 1, 2, and 3 in Fig. 2, across 82,735 human serine and threonine phosphorylation sites confidently identified in mass spectrometry experiments4. The phosphorylated residues in the logos are represented as S/T.
Extended Data Fig. 1 Representation of phosphorylation site motifs in the human serine and threonine phosphoproteome. Venn diagram representation of the percentages of three prominent Ser/Thr kinase motif features, pertaining to Clusters 1, 2, and 3 in Fig. 2, across 82,735 human serine and threonine phosphorylation sites confidently identified in mass spectrometry experiments4. The phosphorylated residues in the logos are represented as S/T.CC BY · PubMed Central
Extended Data Fig. 2 Subcategorization of the basophilic kinases of Cluster 1. Subcategorization of Cluster 1 from Fig. 2 into 11 motif classes.
Extended Data Fig. 2 Subcategorization of the basophilic kinases of Cluster 1. Subcategorization of Cluster 1 from Fig. 2 into 11 motif classes.CC BY · PubMed Central
Extended Data Fig. 3 Subcategorization of the proline-directed kinases of Cluster 2. Subcategorization of Cluster 2 from Fig. 2 into 5 motif classes.
Extended Data Fig. 3 Subcategorization of the proline-directed kinases of Cluster 2. Subcategorization of Cluster 2 from Fig. 2 into 5 motif classes.CC BY · PubMed Central
Extended Data Fig. 4 Subcategorization of the acidophilic kinases of Cluster 3. Subcategorization of Cluster 3 from Fig. 2 into 8 motif classes.
Extended Data Fig. 4 Subcategorization of the acidophilic kinases of Cluster 3. Subcategorization of Cluster 3 from Fig. 2 into 8 motif classes.CC BY · PubMed Central
Extended Data Fig. 5 Structural models of kinase–substrate complexes. a, Synthetic peptide from its complex with PAK4 (PDB: 2Q0N) modelled onto WNK3 (PDB: 5O26). Dotted circle highlights a shallow hydrophobic pocket accommodating a +3 Phe residue. b, GSK3 peptide from its complex with AKT2 (PDB: 1O6L) modelled onto CAMKK2 (PDB: 2ZV2). Circle indicates a hydrophobic pocket that could accommodate a −2 aliphatic residue. c, Monophosphorylated peptide from p63 bound to CK1δ (PDB: 6RU6) modelled onto GRK2 (PDB: 1YM7). Circle shows positive surface potential in the vicinity of the −2 and −3 pSer residues. d, p63 peptide bound to CK1δ (PDB: 6RU8) was modelled onto YANK1 (PDB: 4FR4) showing potential binding sites for −3 and +2 phosphorylated residues. Surface electrostatics are represented with Coulombic potential values were computed in ChimeraX and represented by scale bars (kcal/mol·e).
Extended Data Fig. 5 Structural models of kinase–substrate complexes. a, Synthetic peptide from its complex with PAK4 (PDB: 2Q0N) modelled onto WNK3 (PDB: 5O26). Dotted circle highlights a shallow hydrophobic pocket accommodating a +3 Phe residue. b, GSK3 peptide from its complex with AKT2 (PDB: 1O6L) modelled onto CAMKK2 (PDB: 2ZV2). Circle indicates a hydrophobic pocket that could accommodate a −2 aliphatic residue. c, Monophosphorylated peptide from p63 bound to CK1δ (PDB: 6RU6) modelled onto GRK2 (PDB: 1YM7). Circle shows positive surface potential in the vicinity of the −2 and −3 pSer residues. d, p63 peptide bound to CK1δ (PDB: 6RU8) was modelled onto YANK1 (PDB: 4FR4) showing potential binding sites for −3 and +2 phosphorylated residues. Surface electrostatics are represented with Coulombic potential values were computed in ChimeraX and represented by scale bars (kcal/mol·e).CC BY · PubMed Central
Extended Data Fig. 6 Profiling phospho-acceptor specificity. a, in vitro phosphorylation assays with recombinant kinases and substrate peptides containing either serine or threonine phospho-acceptors. Results shown for 208 recombinant Ser/Thr kinases. b, Correlation plot of the experimental results in (a) with the position sum approach applied in this study to score Ser/Thr phospho-acceptor preference.
Extended Data Fig. 6 Profiling phospho-acceptor specificity. a, in vitro phosphorylation assays with recombinant kinases and substrate peptides containing either serine or threonine phospho-acceptors. Results shown for 208 recombinant Ser/Thr kinases. b, Correlation plot of the experimental results in (a) with the position sum approach applied in this study to score Ser/Thr phospho-acceptor preference.CC BY · PubMed Central
Extended Data Fig. 7 Global analysis of the relationship between the DFG+1 amino acid and preference for the serine versus threonine phospho-acceptor. a, Bottom, relative preferences for Ser or Thr phospho-acceptor residues for each kinase, arranged in order of decreasing Ser/Thr selectivity. Top, frequency of amino acids at the DGF+1 positions of corresponding kinases (bin size: 15 kinases). b, Structural illustration of the proximity between the DFG+1 residue and substrate phospho-acceptor residue, shown with the AKT1 kinase domain bound to substrate (GSK3β) peptide (pdb 1O6K).
Extended Data Fig. 7 Global analysis of the relationship between the DFG+1 amino acid and preference for the serine versus threonine phospho-acceptor. a, Bottom, relative preferences for Ser or Thr phospho-acceptor residues for each kinase, arranged in order of decreasing Ser/Thr selectivity. Top, frequency of amino acids at the DGF+1 positions of corresponding kinases (bin size: 15 kinases). b, Structural illustration of the proximity between the DFG+1 residue and substrate phospho-acceptor residue, shown with the AKT1 kinase domain bound to substrate (GSK3β) peptide (pdb 1O6K).CC BY · PubMed Central
Extended Data Fig. 8 Global performance analysis of substrate percentile scores for their literature-annotated kinases. a, Percentile-score distributions of substrates for their literature-annotated kinases (AUCDF= area under the cumulative distribution function). b, Percentile-score of literature-annotated kinase–substrate pairs as a function of number of reports. Higher number of reports correlates with more favourable percentile-scores between the reported kinase and its substrate. n = 9,073, n = 3,945, n = 544, n = 224, and n = 201 for kinase–substrate relationships with 1, 2, 3, 4, and 5 or more reports, respectively. Statistical analyses were performed using Double-sided Mann-Whitney U-test. Box minima=25th percentile, centre=50th percentile, maxima=75th percentile. Whiskers extend from the box maxima or minima to the largest or smallest value no further than 1.5 x interquartile range. (ns p > 0.05, * p ≤ 0.05, ** p ≤ 10−3, *** p ≤ 10−4, **** p ≤ 10−5).
Extended Data Fig. 8 Global performance analysis of substrate percentile scores for their literature-annotated kinases. a, Percentile-score distributions of substrates for their literature-annotated kinases (AUCDF= area under the cumulative distribution function). b, Percentile-score of literature-annotated kinase–substrate pairs as a function of number of reports. Higher number of reports correlates with more favourable percentile-scores between the reported kinase and its substrate. n = 9,073, n = 3,945, n = 544, n = 224, and n = 201 for kinase–substrate relationships with 1, 2, 3, 4, and 5 or more reports, respectively. Statistical analyses were performed using Double-sided Mann-Whitney U-test. Box minima=25th percentile, centre=50th percentile, maxima=75th percentile. Whiskers extend from the box maxima or minima to the largest or smallest value no further than 1.5 x interquartile range. (ns p > 0.05, * p ≤ 0.05, ** p ≤ 10−3, *** p ≤ 10−4, **** p ≤ 10−5).CC BY · PubMed Central
Extended Data Fig. 9 Global performance analysis of kinase ranks for their literature-annotated substrates. a, Rank distributions of kinases for their literature-annotated substrates (AUCDF= area under the cumulative distribution function). b, Rank of the literature-annotated kinase–substrate pairs, as a function of number of reports. Higher number of reports correlates with more favourable ranking of reported kinase for its substrate. n = 9,073, n = 3,945, n = 544, n = 224, and n = 201 for kinase–substrate relationships with 1, 2, 3, 4, and 5 or more reports, respectively. Statistical analyses were performed using Double-sided Mann-Whitney U-test. Box minima=25th percentile, centre=50th percentile, maxima=75th percentile. Whiskers extend from the box maxima or minima to the largest or smallest value no further than 1.5 x interquartile range. (ns p > 0.05, * p ≤ 0.05, ** p ≤ 10−3, *** p ≤ 10−4, **** p ≤ 10−5).
Extended Data Fig. 9 Global performance analysis of kinase ranks for their literature-annotated substrates. a, Rank distributions of kinases for their literature-annotated substrates (AUCDF= area under the cumulative distribution function). b, Rank of the literature-annotated kinase–substrate pairs, as a function of number of reports. Higher number of reports correlates with more favourable ranking of reported kinase for its substrate. n = 9,073, n = 3,945, n = 544, n = 224, and n = 201 for kinase–substrate relationships with 1, 2, 3, 4, and 5 or more reports, respectively. Statistical analyses were performed using Double-sided Mann-Whitney U-test. Box minima=25th percentile, centre=50th percentile, maxima=75th percentile. Whiskers extend from the box maxima or minima to the largest or smallest value no further than 1.5 x interquartile range. (ns p > 0.05, * p ≤ 0.05, ** p ≤ 10−3, *** p ≤ 10−4, **** p ≤ 10−5).CC BY · PubMed Central
Extended Data Fig. 10 Motif-based scoring results for phosphorylation events facilitated by subcellular localization or docking. a, Illustration of the mitochondrial-localized regulation of the pyruvate dehydrogenase complex through phosphorylation by the PDHKs. Scoring results for PDHA1 Ser293, highlighting the PDHK family kinases. b, Illustration of docking-driven phosphorylation of ERK1/2 by MEK1/2. Scoring results for ERK1 Thr202 / ERK2 Thr203 (identical sequences), highlighting MEK1 and MEK2.
Extended Data Fig. 10 Motif-based scoring results for phosphorylation events facilitated by subcellular localization or docking. a, Illustration of the mitochondrial-localized regulation of the pyruvate dehydrogenase complex through phosphorylation by the PDHKs. Scoring results for PDHA1 Ser293, highlighting the PDHK family kinases. b, Illustration of docking-driven phosphorylation of ERK1/2 by MEK1/2. Scoring results for ERK1 Thr202 / ERK2 Thr203 (identical sequences), highlighting MEK1 and MEK2.CC BY · PubMed Central
Extended Data Fig. 11 Scoring comparison of CDK subfamilies. Illustration of the phosphoregulation of RNA Polymerase II (POLR2A) CTD and Retinoblastoma protein (Rb) by their respective canonical CDKs, the transcriptional CDKs (purple) and the cell cycle progression CDKs (green) (left). Links between kinases and substrates correspond to favourable scores between motifs and phosphorylation sites (right).
Extended Data Fig. 11 Scoring comparison of CDK subfamilies. Illustration of the phosphoregulation of RNA Polymerase II (POLR2A) CTD and Retinoblastoma protein (Rb) by their respective canonical CDKs, the transcriptional CDKs (purple) and the cell cycle progression CDKs (green) (left). Links between kinases and substrates correspond to favourable scores between motifs and phosphorylation sites (right).CC BY · PubMed Central

Verification

Claims the paper makes

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Claims other papers make about this one

Each is checked against this paper’s own text: does it say what the citing paper says it does?

Self-check

Percentages, sums and reported differences re-computed from the paper’s own numbers.

The self-check needs the paper’s full text.

Citations

Cites

Show all 120 references

Cited by

Sorted by citations

See all 734

Metadata

Authors
Johnson, J. L., Yaron, T. M., Huntsman, E. M., Kerelsky, A., Song, J., Regev, A., Lin, T.-Y., Liberatore, K., Cizin, D. M., Cohen, B. M., Vasan, N., Ma, Y., Krismer, K., Robles, J. T., van de Kooij, B., van Vlimmeren, A. E., Andrée-Busch, N., Käufer, N. F., Dorovkov, M. V., Ryazanov, A. G., Takagi, Y., Kastenhuber, E. R., Goncalves, M. D., Hopkins, B. D., Elemento, O., Taatjes, D. J., Maucuer, A., Yamashita, A., Degterev, A., Uduman, M., Lu, J., Landry, S. D., Zhang, B., Cossentino, I., Linding, R., Blenis, J., Hornbeck, P. V., Turk, B. E., Yaffe, M. B., & Cantley, L. C.
Affiliations
  • Weill Cornell Medicine (US)
  • Memorial Sloan Kettering Cancer Center (US)
  • Cornell University (US)
  • Columbia University Irving Medical Center (US)
  • MIT Computer Science and Artificial Intelligence Laboratory (US)
  • Yale University (US)
  • Massachusetts Institute of Technology (US)
  • Technische Universität Braunschweig (DE)
  • Rutgers, The State University of New Jersey (US)
  • Indiana University School of Medicine (US)
  • Icahn School of Medicine at Mount Sinai (US)
  • University of Colorado Boulder (US)
  • Inserm (FR)
  • University of the Ryukyus (JP)
  • Tufts University (US)
  • Cell Signaling Technology (United States) (US)
  • Humboldt-Universität zu Berlin (DE)
  • National Institutes of Health (US)
Corresponding author
  • Turk, Benjamin E.
  • Yaffe, Michael B.
  • Cantley, Lewis C.
Venue
Nature
Published
2023-01-11
Volume
613
Issue
7945
Pages
759-766
ISSN
  • 0028-0836 (ISSN-L)
  • 1476-4687
OpenAlex
W4315642941
Cited by
728
References
120
Citations by year
2026: 206 · 2025: 223 · 2024: 193 · 2023: 102
FWCI
67.03
Type
journal article
Language
en
Open access
Locations
  • Nature (publishedVersion, cc-by)
  • bioRxiv (Cold Spring Harbor Laboratory) (acceptedVersion, cc-by-nc)
  • PubMed Central (publishedVersion, cc-by)
  • PubMed (publishedVersion)
  • HAL (submittedVersion, cc-by)
  • RePEc: Research Papers in Economics (submittedVersion)
  • Institutional Repositories DataBase (IRDB) (submittedVersion, cc-by)
  • University of Groningen research database (University of Groningen / Centre for Information Technology) (submittedVersion, cc-by)
  • IUScholarWorks (Indiana University) (submittedVersion, cc-by)
  • University oh the Ryukyus Repository (submittedVersion, cc-by-nc-nd)
Funders
  • Cancer Research UK (A28596, C42454/A28596)
  • National Cancer Institute (CA120964, CA14051, T32CA203702, CA197588, CA226898, R01CA226898, Grant P30-CA14051, R35CA197588, P01CA120964, P30CA014051, P01CA117969, P30-CA-14051)
  • National Institute of General Medical Sciences (R01GM104047, R35GM139550)
  • National Institute of Diabetes and Digestive and Kidney Diseases (U24DK116204)
  • National Institute of Environmental Health Sciences (R35ES028374)
Topics
  • Physical Sciences › Chemistry › Spectroscopy › Advanced Proteomics Techniques and Applications
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Protein Kinase Regulation and GTPase Signaling
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Protein Tyrosine Phosphatases
Keywords
serine/threonine kinase, kinase substrate specificity, protein phosphorylation, phosphoproteome, synthetic peptide libraries, kinase profiling
MeSH terms
Cell Line, Humans, Phosphoproteins, Phosphorylation, Phosphoserine, Phosphothreonine, Serine, Substrate Specificity, Threonine, Protein Serine-Threonine Kinases, Proteome, Datasets as Topic

About these checks

Registry, retraction, Cites and statistics results are computed, not judged, with one exception: a printed reference that no index lists is read once by a model before it counts as a citation error. Claims are extracted by a model and kept only when their quote appears verbatim in the text; each cited source is then checked, with a model’s help, for whether it says what the claim relies on. Those verdicts can be wrong — read the claim’s quote and the cited source before relying on one.

Statistics and arithmetic are recomputed, not judged.

Confidence bands are shown only for fields calibrated on our annotated test set. None is calibrated yet, so no band is shown.

What counts as a claim ›How we measure confidence ›How we check this ›

Are you an author, or do you see a problem with this page? Request a review