The intrinsic substrate specificity of the human tyrosine kinome.

Phosphorylation of proteins on tyrosine (Tyr) residues evolved in metazoan organisms as a mechanism of coordinating tissue growth1. Multicellular eukaryotes typically have more than 50 distinct protein Tyr kinases that catalyse the phosphorylation of thousands of Tyr residues throughout the proteome1-3. How a given Tyr kinase can phosphorylate a specific subset of proteins at unique Tyr sites is only partially understood4-7. Here we used combinatorial peptide arrays to profile the substrate sequence specificity of all human Tyr kinases. Globally, the Tyr kinases demonstrate considerable diversity in optimal patterns of residues surrounding the site of phosphorylation, revealing the functional organization of the human Tyr kinome by substrate motif preference. Using this information, Tyr kinases that are most compatible with phosphorylating any Tyr site can be identified. Analysis of mass spectrometry phosphoproteomic datasets using this compendium of kinase specificities accurately identifies specific Tyr kinases that are dysregulated in cells after stimulation with growth factors, treatment with anti-cancer drugs or expression of oncogenic variants. Furthermore, the topology of known Tyr signalling networks naturally emerged from a comparison of the sequence specificities of the Tyr kinases and the SH2 phosphotyrosine (pTyr)-binding domains. Finally we show that the intrinsic substrate specificity of Tyr kinases has remained fundamentally unchanged from worms to humans, suggesting that the fidelity between Tyr kinases and their protein substrate sequences has been maintained across hundreds of millions of years of evolution.

Proteomic and Phosphoproteomic Reprogramming in Epithelial Ovarian Cancer Metastasis.

Epithelial ovarian cancer (EOC) is a high-risk cancer presenting with heterogeneous tumors. The high incidence of EOC metastasis from primary tumors to nearby tissues and organs is a major driver of EOC lethality. We used cellular models of spheroid formation and readherence to investigate cellular signaling dynamics in each step toward EOC metastasis. In our system, adherent cells model primary tumors, spheroid formation represents the initiation of metastatic spread, and readherent spheroid cells represent secondary tumors. Proteomic and phosphoproteomic analyses show that spheroid cells are hypoxic and show markers for cell cycle arrest. Aurora kinase B abundance and downstream substrate phosphorylation are significantly reduced in spheroids and readherent cells, explaining their cell cycle arrest phenotype. The proteome of readherent cells is most similar to spheroids, yet greater changes in the phosphoproteome show that spheroid cells stimulate Rho-associated kinase 1 (ROCK1)-mediated signaling, which controls cytoskeletal organization. In spheroids, we found significant phosphorylation of ROCK1 substrates that were reduced in both adherent and readherent cells. Application of the ROCK1-specific inhibitor Y-27632 to spheroids increased the rate of readherence and altered spheroid density. The data suggest ROCK1 inhibition increases EOC metastatic potential. We identified novel pathways controlled by Aurora kinase B and ROCK1 as major drivers of metastatic behavior in EOC cells. Our data show that phosphoproteomic reprogramming precedes proteomic changes that characterize spheroid readherence in EOC metastasis.

Immune-based classification of HPV-associated oropharyngeal cancer with implications for biomarker-driven treatment de-intensification.

BACKGROUND: There is significant interest in treatment de-escalation for human papillomavirus-associated (HPV+) oropharyngeal squamous cell carcinoma (OPSCC) patients given the generally favourable prognosis. However, 15-30% of patients recur after primary treatment, reflecting a need for improved risk-stratification tools. We sought to develop a molecular test to risk stratify HPV+ OPSCC patients. METHODS: We created an immune score (UWO3) associated with survival outcomes in six independent cohorts comprising 906 patients, including blinded retrospective and prospective external validations. Two aggressive radiation de-escalation cohorts were used to assess the ability of UWO3 to identify patients who recur. Multivariate Cox models were used to assess the associations between the UWO3 immune class and outcomes. FINDINGS: A three-gene immune score classified patients into three immune classes (immune rich, mixed, or immune desert) and was strongly associated with disease-free survival in six datasets, including large retrospective and prospective datasets. Pooled analysis demonstrated that the immune rich group had superior disease-free survival compared to the immune desert (HR = 9.0, 95% CI: 3.2-25.5, P = 3.6 × 10-5) and mixed (HR = 6.4, 95% CI: 2.2-18.7, P = 0.006) groups after adjusting for age, sex, smoking status, and AJCC8 clinical stage. Finally, UWO3 was able to identify patients from two small treatment de-escalation cohorts who remain disease-free after aggressive de-escalation to 30 Gy radiation. INTERPRETATION: With additional prospective validation, the UWO3 score could enable biomarker-driven clinical decision-making for patients with HPV+ OPSCC based on robust outcome prediction across six independent cohorts. Prospective de-escalation and intensification clinical trials are currently being planned. FUNDING: CIHR, European Union, and the NIH.

Generic Plug-and-Play Strategy for High-Throughput Analysis of PTM-Mediated Protein Complexes.

Protein complexes mediated by various post-translational modifications (PTMs) play important roles in almost every aspect of biological processes. PTM-mediated protein complexes often have weak and transient binding properties, which limit their unbiased profiling especially in complex biological samples. Here, we developed a plug-and-play chemical proteomic approach for high-throughput analyis of PTM-mediated protein complexes. Taking advantage of the glutathione-S-transferase (GST) tag, which is the gold standard for protein purification and has wide access to a variety of proteins of interest (POIs), a glutathione (GSH) group- and photo-cross-linking group-containing trifunctional chemical probe was developed to tag POIs and assembled onto a streptavidin-coated 96-well plate for affinity purification, photo-cross-linking, and proteomics sample preparation in a fully integrated manner. Compared with the previously developed photo-pTyr-scaffold strategy, by assembling the tyrosine phosphorylation (pTyr) binding domain through covalent NHS chemistry, the new plug-and-play strategy using a noncovalent GST-GSH interaction has comparable enrichment efficiency for EGF stimulation-dependent pTyr protein complexes. To further prove its feasibility, we additionally assembled four pTyr-binding domains in the 96-well plate and selectively identified their pTyr-dependent interacting proteins. Importantly, we systematically optimized and applied the plug-and-play approach for exploring protein methylation-mediated protein complexes, which are difficult to be characterized due to their weak binding affinity and the lack of efficient enrichment strategies. We explored a comprehensive protein methylation-mediated interaction network assembled by five protein methylation binding domains including the chromo domain of MPP8, tandem tudor domain of KDM4A, full-length CBX1, PHD domain of RAG2, and tandem tudor domain of TP53BP1 and validated the chromo domain- and tudor domain-mediated interaction with histone H3. Collectively, this plug-and-play approach provides a convenient and generic strategy for exploring PTM-dependent protein complexes for any POIs with the GST tag.

Discovery of an exosite on the SOCS2-SH2 domain that enhances SH2 binding to phosphorylated ligands.

Suppressor of cytokine signaling (SOCS)2 protein is a key negative regulator of the growth hormone (GH) and Janus kinase (JAK)-Signal Transducers and Activators of Transcription (STAT) signaling cascade. The central SOCS2-Src homology 2 (SH2) domain is characteristic of the SOCS family proteins and is an important module that facilitates recognition of targets bearing phosphorylated tyrosine (pTyr) residues. Here we identify an exosite on the SOCS2-SH2 domain which, when bound to a non-phosphorylated peptide (F3), enhances SH2 affinity for canonical phosphorylated ligands. Solution of the SOCS2/F3 crystal structure reveals F3 as an α-helix which binds on the opposite side of the SH2 domain to the phosphopeptide binding site. F3:exosite binding appears to stabilise the SOCS2-SH2 domain, resulting in slower dissociation of phosphorylated ligands and consequently, enhances binding affinity. This biophysical enhancement of SH2:pTyr binding affinity translates to increase SOCS2 inhibition of GH signaling.

Dynamic interplay of two molecular switches enabled by the MEK1/2-ERK1/2 and IL-6-STAT3 signaling axes controls epithelial cell migration in response to growth factors.

Cell migration is an essential physiological process, and aberrant migration of epithelial cells underlies many pathological conditions. However, the molecular mechanisms governing cell migration are not fully understood. We report here that growth factor-induced epithelial cell migration is critically dependent on the crosstalk of two molecular switches, namely phosphorylation switch (P-switch) and transcriptional switch (T-switch). P-switch refers to dynamic interactions of deleted in liver cancer 1 (DLC1) and PI3K with tensin-3 (TNS3), phosphatase and tensin homolog (PTEN), C-terminal tension, and vav guanine nucleotide exchange factor 2 (VAV2) that are dictated by mitogen-activated protein kinase kinase 1/2-extracellular signal-regulated protein kinase 1/2-dependent phosphorylation of TNS3, PTEN, and VAV2. Phosphorylation of TNS3 and PTEN on specific Thr residues led to the switch of DLC1-TNS3 and PI3K-PTEN complexes to DLC1-PTEN and PI3K-TNS3 complexes, whereas Ser phosphorylation of VAV2 promotes the transition of the PI3K-TNS3/PTEN complexes to PI3K-VAV2 complex. T-switch denotes an increase in C-terminal tension transcription/expression regulated by both extracellular signal-regulated protein kinase 1/2 and signal transducer and activator of transcription 3 (STAT3) via interleukin-6-Janus kinase-STAT3 signaling pathway. We have found that, the P-switch is indispensable for both the initiation and continuation of cell migration induced by growth factors, whereas the T-switch is only required to sustain cell migration. The interplay of the two switches facilitated by the interleukin-6-Janus kinase-STAT3 pathway governs a sequence of dynamic protein-protein interactions for sustained cell migration. That a similar mechanism is employed by both normal and tumorigenic epithelial cells to drive their respective migration suggests that the P-switch and T-switch are general regulators of epithelial cell migration and potential therapeutic targets.

Rapid and accurate agglutination-based testing for SARS-CoV-2 antibodies.

We have developed a rapid, accurate, and cost-effective serologic test for SARS-CoV-2 virus, which caused the COVID-19 pandemic, on the basis of antibody-dependent agglutination of antigen-coated latex particles. When validated using plasma samples that are positive or negative for SARS-CoV-2, the agglutination assay detected antibodies against the receptor-binding domain of the spike (S-RBD) or the nucleocapsid protein of SARS-CoV-2 with 100% specificity and ∼98% sensitivity. Furthermore, we found that the strength of the S-RBD antibody response measured by the agglutination assay correlated with the efficiency of the plasma in blocking RBD binding to the angiotensin-converting enzyme 2 in a surrogate neutralization assay, suggesting that the agglutination assay might be used to identify individuals with virus-neutralizing antibodies. Intriguingly, we found that >92% of patients had detectable antibodies on the day of a positive viral RNA test, suggesting that the agglutination antibody test might complement RNA testing for the diagnosis of SARS-CoV-2 infection.

HDAC8 Activates AKT through Upregulating PLCB1 and Suppressing DESC1 Expression in MEK1/2 Inhibition-Resistant Cells.

Inhibition of the RAF-MEK1/2-ERK signaling pathway is an ideal strategy for treating cancers with NRAS or BRAF mutations. However, the development of resistance due to incomplete inhibition of the pathway and activation of compensatory cell proliferation pathways is a major impediment of the targeted therapy. The anthrax lethal toxin (LT), which cleaves and inactivates MEKs, is a modifiable biomolecule that can be delivered selectively to tumor cells and potently kills various tumor cells. However, resistance to LT and the mechanism involved are yet to be explored. Here, we show that LT, through inhibiting MEK1/2-ERK activation, inhibits the proliferation of cancer cells with NRAS/BRAF mutations. Among them, the human colorectal tumor HT-29 and murine melanoma B16-BL6 cells developed resistance to LT in 2 to 3 days of treatment. These resistant cells activated AKT through a histone deacetylase (HDAC) 8-dependent pathway. Using an Affymetrix microarray, followed by qPCR validation, we identified that the differential expression of the phospholipase C-ß1 (PLCB1) and squamous cell carcinoma-1 (DESC1) played an important role in HDAC8-mediated AKT activation and resistance to MEK1/2-ERK inhibition. By using inhibitors, small interference RNAs and/or expression vectors, we found that the inhibition of HDAC8 suppressed PLCB1 expression and induced DESC1 expression in the resistant cells, which led to the inhibition of AKT and re-sensitization to LT and MEK1/2 inhibition. These results suggest that targeting PLCB1 and DESC1 is a novel strategy for inhibiting the resistance to MEK1/2 inhibition.

Single-Cell Analysis of Signaling Proteins Provides Insights into Proapoptotic Properties of Anticancer Drugs.

Single-cell DNA analysis technology has provided unprecedented insights into many physiological and pathological processes. In contrast, technologies that allow protein analysis in single cells have lagged behind. Herein, a method called single-cell Plasmonic ImmunoSandwich Assay (scPISA) that is capable of measuring signaling proteins and protein complexes in single living cells is described. scPISA is straightforward, comprising specific in-cell extraction and ultrasensitive plasmonic detection. It is applied to evaluate the efficacy and kinetics of cytotoxic drugs. It reveals that different drugs exhibit distinct proapoptotic properties at the single-cell level. A set of new parameters is thus proposed for comprehensive and quantitative evaluation of the efficacy of anticancer drugs. It discloses that metformin can dramatically enhance the overall anticancer efficacy when combined with actinomycin D, although it itself is significantly less effective. Furthermore, scPISA reveals that survivin interacts with cytochrome C and caspase-3 in a dynamic fashion in single cells during continuous drug treatment. As compared with conventional assays, scPISA exhibits several significant advantages, such as ultrahigh sensitivity, single-cell resolution, fast speed, and so on. Therefore, this approach may provide a powerful tool for wide, important applications from basic research to clinical applications, particularly precision medicine.

MMS2plot: An R Package for Visualizing Multiple MS/MS Spectra for Groups of Modified and Non-Modified Peptides.

A large number of post-translational modifications (PTMs) in proteins are buried in the unassigned mass spectrometric (MS) spectra in shot-gun proteomics datasets. Because the modified peptide fragments are low in abundance relative to the corresponding non-modified versions, it is critical to develop tools that allow facile evaluation of assignment of PTMs based on the MS/MS spectra. Such tools will preferably have the ability to allow comparison of fragment ion spectra and retention time between the modified and unmodified peptide pairs or group. Herein, MMS2plot, an R package for visualizing peptide-spectrum matches (PSMs) for multiple peptides, is described. MMS2plot features a batch mode and generates the output images in vector graphics file format that facilitate evaluation and publication of the PSM assignment. MMS2plot is expected to play an important role in PTM discovery from large-scale proteomics datasets generated by liquid chromatography-MS/MS. The MMS2plot package is freely available at https://github.com/lileir/MMS2plot under the GPL-3 license.

Proteome-wide Prediction of Lysine Methylation Leads to Identification of H2BK43 Methylation and Outlines the Potential Methyllysine Proteome.

Protein Lys methylation plays a critical role in numerous cellular processes, but it is challenging to identify Lys methylation in a systematic manner. Here we present an approach combining in silico prediction with targeted mass spectrometry (MS) to identify Lys methylation (Kme) sites at the proteome level. We develop MethylSight, a program that predicts Kme events solely on the physicochemical properties of residues surrounding the putative methylation sites, which then requires validation by targeted MS. Using this approach, we identify 70 new histone Kme marks with a 90% validation rate. H2BK43me2, which undergoes dynamic changes during stem cell differentiation, is found to be a substrate of KDM5b. Furthermore, MethylSight predicts that Lys methylation is a prevalent post-translational modification in the human proteome. Our work provides a useful resource for guiding systematic exploration of the role of Lys methylation in human health and disease.

DLC1 SAM domain-binding peptides inhibit cancer cell growth and migration by inactivating RhoA.

Deleted-in-liver cancer 1 (DLC1) exerts its tumor suppressive function mainly through the Rho-GTPase-activating protein (RhoGAP) domain. When activated, the domain promotes the hydrolysis of RhoA-GTP, leading to reduced cell migration. DLC1 is kept in an inactive state by an intramolecular interaction between its RhoGAP domain and the DLC1 sterile α motif (SAM) domain. We have shown previously that this autoinhibited state of DLC1 may be alleviated by tensin-3 (TNS3) or PTEN. We show here that the TNS3/PTEN-DLC1 interactions are mediated by the C2 domains of the former and the SAM domain of the latter. Intriguingly, the DLC1 SAM domain was capable of binding to specific peptide motifs within the C2 domains. Indeed, peptides containing the binding motifs were highly effective in blocking the C2-SAM domain-domain interaction. Importantly, when fused to the tat protein-transduction sequence and subsequently introduced into cells, the C2 peptides potently promoted the RhoGAP function in DLC1, leading to decreased RhoA activation and reduced tumor cell growth in soft agar and migration in response to growth factor stimulation. To facilitate the development of the C2 peptides as potential therapeutic agents, we created a cyclic version of the TNS3 C2 domain-derived peptide and showed that this peptide readily entered the MDA-MB-231 breast cancer cells and effectively inhibited their migration. Our work shows, for the first time, that the SAM domain is a peptide-binding module and establishes the framework on which to explore DLC1 SAM domain-binding peptides as potential therapeutic agents for cancer treatment.

Single-Layer ß Sheets as Versatile Protein-Binding Modules.

In this issue of Structure, Li et al. (2019) describe the structure of the MORN4-Myo3a complex, revealing that the MORN repeats in MORN4 form a single-layer antiparallel ß sheet and employs an extensive surface area in binding Myo3a. Their findings suggest that single-layer ß sheets are versatile protein-binding modules.

Hybrid CuCoO-GO enables ultrasensitive detection of antibiotics with enhanced laser desorption/ionization at nano-interfaces.

The soaring concerns globally on antibiotic overuse have made calls for the development of rapid and sensitive detection methods urgent. Here we report that the hybrid CuCoO-GO matrix allows for sensitive detection of various antibiotics in combination with MALDI TOF MS. The new matrix is composed of few-layered GO nanosheets decorated with CuCoO nanoparticles with an average size of 10 nm, and exhibits excellent aqueous suspensibility. Accurate quantitation of the sulfonamide antibiotics in milk samples have been demonstrated using a CuCoO-GO matrix and a stable isotope (C13)-labeled analyte as the internal standard. Our experiments have achieved lower limits of detection (LOD) by several hundred fold for the detection of a panel of representative antibiotics, in comparison with the literature reports. Both intrabacterial and extrabacterial residual antibiotics can be sensitively detected with our method. We have further investigated the molecular mechanism of the enhanced desorption/ionization efficiency by the CuCoO-GO matrix with synchrotron radiation techniques for the first time. This work provides a sensitive matrix enabling MALDI-TOF MS to be applied in small molecular analysis, but also presents a distinct perspective on the mechanism behind the material functions.

N-acetylcysteine decreases malignant characteristics of glioblastoma cells by inhibiting Notch2 signaling.

BACKGROUND: Glioblastomas multiforme (GBM) is the most devastating primary intracranial malignancy lacking effective clinical treatments. Notch2 has been established to be a prognostic marker and probably involved in GBM malignant progression. N-acetylcysteine (NAC), a precursor of intracellular glutathione (GSH), has been widely implicated in prevention and therapy of several cancers. However, the role of NAC in GBM remains unclear and the property of NAC independent of its antioxidation is largely unknown. METHODS: The mRNA and protein levels of Notch family and other related factors were detected by RT-PCR and western blot, respectively. In addition, intracellular reactive oxygen species (ROS) was measured by flow cytometry-based DCFH-DA. Moreover, cell viability was assessed by CCK8 and cell cycle was analyzed by flow cytometry-based PI staining. The level of apoptosis was checked by flow cytometry-based Annexin V/PI. Cell migration and invasion were evaluated by wound healing and transwell invasion assays. At last, U87 Xenograft model was established to confirm whether NAC could restrain the growth of tumor. RESULTS: Our data showed that NAC could decrease the protein level of Notch2. Meanwhile, NAC had a decreasing effect on the mRNA and protein levels of its downstream targets Hes1 and Hey1. These effects caused by NAC were independent of cellular GSH and ROS levels. The mechanism of NAC-mediated Notch2 reduction was elucidated by promoting Notch2 degradation through Itch-dependent lysosome pathway. Furthermore, NAC could prevent proliferation, migration, and invasion and might induce apoptosis in GBM cells via targeting Notch2. Significantly, NAC could suppress the growth of tumor in vivo. CONCLUSIONS: NAC could facilitate Notch2 degradation through lysosomal pathway in an antioxidant-independent manner, thus attenuating Notch2 malignant signaling in GBM cells. The remarkable ability of NAC to inhibit cancer cell proliferation and tumor growth may implicate a novel application of NAC on GBM therapy.

Surface Loops in a Single SH2 Domain Are Capable of Encoding the Spectrum of Specificity of the SH2 Family.

Src homology 2 (SH2) domains play an essential role in cellular signal transduction by binding to proteins phosphorylated on Tyr residue. Although Tyr phosphorylation (pY) is a prerequisite for binding for essentially all SH2 domains characterized to date, different SH2 domains prefer specific sequence motifs C-terminal to the pY residue. Because all SH2 domains adopt the same structural fold, it is not well understood how different SH2 domains have acquired the ability to recognize distinct sequence motifs. We have shown previously that the EF and BG loops that connect the secondary structure elements on an SH2 domain dictate its specificity. In this study, we investigated if these surface loops could be engineered to encode diverse specificities. By characterizing a group of SH2 variants selected by different pY peptides from phage-displayed libraries, we show that the EF and BG loops of the Fyn SH2 domain can encode a wide spectrum of specificities, including all three major specificity classes (p + 2, p + 3 and p + 4) of the SH2 domain family. Furthermore, we found that the specificity of a given variant correlates with the sequence feature of the bait peptide used for its isolation, suggesting that an SH2 domain may acquire specificity by co-evolving with its ligand. Intriguingly, we found that the SH2 variants can employ a variety of different mechanisms to confer the same specificity, suggesting the EF and BG loops are highly flexible and adaptable. Our work provides a plausible mechanism for the SH2 domain to acquire the wide spectrum of specificity observed in nature through loop variation with minimal disturbance to the SH2 fold. It is likely that similar mechanisms may have been employed by other modular interaction domains to generate diversity in specificity.

Engineered SH2 domains with tailored specificities and enhanced affinities for phosphoproteome analysis.

Protein phosphorylation is the most abundant post-translational modification in cells. Src homology 2 (SH2) domains specifically recognize phosphorylated tyrosine (pTyr) residues to mediate signaling cascades. A conserved pocket in the SH2 domain binds the pTyr side chain and the EF and BG loops determine binding specificity. By using large phage-displayed libraries, we engineered the EF and BG loops of the Fyn SH2 domain to alter specificity. Engineered SH2 variants exhibited distinct specificity profiles and were able to bind pTyr sites on the epidermal growth factor receptor, which were not recognized by the wild-type Fyn SH2 domain. Furthermore, mass spectrometry showed that SH2 variants with additional mutations in the pTyr-binding pocket that enhanced affinity were highly effective for enrichment of diverse pTyr peptides within the human proteome. These results showed that engineering of the EF and BG loops could be used to tailor SH2 domain specificity, and SH2 variants with diverse specificities and high affinities for pTyr residues enabled more comprehensive analysis of the human phosphoproteome. STATEMENT: Src Homology 2 (SH2) domains are modular domains that recognize phosphorylated tyrosine embedded in proteins, transducing these post-translational modifications into cellular responses. Here we used phage display to engineer hundreds of SH2 domain variants with altered binding specificities and enhanced affinities, which enabled efficient and differential enrichment of the human phosphoproteome for analysis by mass spectrometry. These engineered SH2 domain variants will be useful tools for elucidating the molecular determinants governing SH2 domains binding specificity and for enhancing analysis and understanding of the human phosphoproteome.

Affinity Purification of Methyllysine Proteome by Site-Specific Covalent Conjugation.

A basic but critical step in targeted proteomics by mass spectrometry is the separation of the targeted proteins from the complex mixture of the whole proteome by affinity purification. The bait protein is usually immobilized on the surface of a solid support to enable affinity-based purification of the targeted proteome. Here, we developed a site-specific covalent immobilization of the bait protein through affinity-guided covalent coupling (AGCC) of a single cysteine residue of an SH2 domain (utilized as an affinity tag for the protein target) with an engineered ligand peptide. Site-specific covalent immobilization of a methyllysine-binding protein HP1ß chromodomain on the agarose resin was used to purify the methyllysine proteome from the whole-protein mixture. This new bait immobilization led to a notably low background in the affinity purification step, markedly outperforming the conventional (His)6 tag-nickel nitrilotriacetic acid (Ni-NTA) immobilization method. Subsequent analysis of the purified proteome identified 275 lysine methylated sites and 184 methylated proteins from 332 HP1ß CD-binding proteins, including 30 novel methylated proteins. This work demonstrates that a robust site-specific covalent protein immobilization method is well-suited for proteomic analysis of low-abundance proteins. This method also enables the identification of new methylated proteins and methylation sites in the methyllysine proteome.