Systematic identification of phosphorylation-mediated protein interaction switches

Waves of regulated protein expression and phosphorylation rewire the proteome to drive gametogenesis in budding yeast

Sexually reproducing eukaryotes employ a developmentally regulated cell division program-meiosis-to generate haploid gametes from diploid germ cells. To understand how gametes arise, we generated a proteomic census encompassing the entire meiotic program of budding yeast. We found that concerted waves of protein expression and phosphorylation modify nearly all cellular pathways to support meiotic entry, meiotic progression, and gamete morphogenesis. Leveraging this comprehensive resource, we pinpointed dynamic changes in mitochondrial components and showed that phosphorylation of the FoF1-ATP synthase complex is required for efficient gametogenesis. Furthermore, using cryoET as an orthogonal approach to visualize mitochondria, we uncovered highly ordered filament arrays of Ald4ALDH2, a conserved aldehyde dehydrogenase that is highly expressed and phosphorylated during meiosis. Notably, phosphorylation-resistant mutants failed to accumulate filaments, suggesting that phosphorylation regulates context-specific Ald4ALDH2 polymerization. Overall, this proteomic census constitutes a broad resource to guide the exploration of the unique sequence of events underpinning gametogenesis.

Phosphorylation of disordered proteins tunes local and global intramolecular interactions

Protein post-translational modifications, such as phosphorylation, are important regulatory signals for diverse cellular functions. In particular, intrinsically disordered protein regions (IDRs) are subject to phosphorylation as a means to modulate their interactions and functions. Toward understanding the relationship between phosphorylation in IDRs and specific functional outcomes, we must consider how phosphorylation affects the IDR conformational ensemble. Various experimental techniques are suited to interrogate the features of IDR ensembles; molecular simulations can provide complementary insights and even illuminate ensemble features that may be experimentally inaccessible. Therefore, we sought to expand the tools available to study phosphorylated IDRs by all-atom Monte Carlo simulations. To this end, we implemented parameters for phosphoserine (pSer) and phosphothreonine (pThr) into the OPLS version of the continuum solvent model, ABSINTH, and assessed their performance in all-atom simulations compared to published findings. We simulated short (< 20 residues) and long (> 80 residues) phospho-IDRs that, collectively, survey both local and global phosphorylation-induced changes to the ensemble. Our simulations of four well-studied phospho-IDRs show near-quantitative agreement with published findings for these systems via metrics including changes to radius of gyration, transient helicity, and persistence length. We also leveraged the inherent advantage of sequence control in molecular simulations to explore the conformational effects of diverse combinations of phospho-sites in two multi-phosphorylated IDRs. Our results support and expand on prior observations that connect phosphorylation to changes in the IDR conformational ensemble. Herein, we describe phosphorylation as a means to alter sequence chemistry, net charge and charge patterning, and intramolecular interactions, which can collectively modulate the local and global IDR ensemble features.

Large-scale phosphomimetic screening identifies phospho-modulated motif-based protein interactions

Phosphorylation is a ubiquitous post-translation modification that regulates protein function by promoting, inhibiting or modulating protein-protein interactions. Hundreds of thousands of phosphosites have been identified but the vast majority have not been functionally characterised and it remains a challenge to decipher phosphorylation events modulating interactions. We generated a phosphomimetic proteomic peptide-phage display library to screen for phosphosites that modulate short linear motif-based interactions. The peptidome covers ~13,500 phospho-serine/threonine sites found in the intrinsically disordered regions of the human proteome. Each phosphosite is represented as wild-type and phosphomimetic variant. We screened 71 protein domains to identify 248 phosphosites that modulate motif-mediated interactions. Affinity measurements confirmed the phospho-modulation of 14 out of 18 tested interactions. We performed a detailed follow-up on a phospho-dependent interaction between clathrin and the mitotic spindle protein hepatoma-upregulated protein (HURP), demonstrating the essentiality of the phospho-dependency to the mitotic function of HURP. Structural characterisation of the clathrin-HURP complex elucidated the molecular basis for the phospho-dependency. Our work showcases the power of phosphomimetic ProP-PD to discover novel phospho-modulated interactions required for cellular function.

Sequence-based machine learning method for predicting the effects of phosphorylation on protein-protein interactions

Protein phosphorylation, catalyzed by kinases, is an important biochemical process, which plays an essential role in multiple cell signaling pathways. Meanwhile, protein-protein interactions (PPI) constitute the signaling pathways. Abnormal phosphorylation status on protein can regulate protein functions through PPI to evoke severe diseases, such as Cancer and Alzheimer's disease. Due to the limited experimental evidence and high costs to experimentally identify novel evidence of phosphorylation regulation on PPI, it is necessary to develop a high-accuracy and user-friendly artificial intelligence method to predict phosphorylation effect on PPI. Here, we proposed a novel sequence-based machine learning method named PhosPPI, which achieved better identification performance (Accuracy and AUC) than other competing predictive methods of Betts, HawkDock and FoldX. PhosPPI is now freely available in web server (https://phosppi.sjtu.edu.cn/). This tool can help the user to identify functional phosphorylation sites affecting PPI and explore phosphorylation-associated disease mechanism and drug development.

A basic phosphoproteomic-DIA workflow integrating precise quantification of phosphosites in systems biology

Phosphorylation is one of the most important post-translational modifications (PTMs) of proteins, governing critical protein functions. Most human proteins have been shown to undergo phosphorylation, and phosphoproteomic studies have been widely applied due to recent advancements in high-resolution mass spectrometry technology. Although the experimental workflow for phosphoproteomics has been well-established, it would be useful to optimize and summarize a detailed, feasible protocol that combines phosphoproteomics and data-independent acquisition (DIA), along with follow-up data analysis procedures due to the recent instrumental and bioinformatic advances in measuring and understanding tens of thousands of site-specific phosphorylation events in a single experiment. Here, we describe an optimized Phos-DIA protocol, from sample preparation to bioinformatic analysis, along with practical considerations and experimental configurations for each step. The protocol is designed to be robust and applicable for both small-scale phosphoproteomic analysis and large-scale quantification of hundreds of samples for studies in systems biology and systems medicine.

Unified access to up-to-date residue-level annotations from UniProtKB and other biological databases for PDB data

More than 61,000 proteins have up-to-date correspondence between their amino acid sequence (UniProtKB) and their 3D structures (PDB), enabled by the Structure Integration with Function, Taxonomy and Sequences (SIFTS) resource. SIFTS incorporates residue-level annotations from many other biological resources. SIFTS data is available in various formats like XML, CSV and TSV format or also accessible via the PDBe REST API but always maintained separately from the structure data (PDBx/mmCIF file) in the PDB archive. Here, we extended the wwPDB PDBx/mmCIF data dictionary with additional categories to accommodate SIFTS data and added the UniProtKB, Pfam, SCOP2, and CATH residue-level annotations directly into the PDBx/mmCIF files from the PDB archive. With the integrated UniProtKB annotations, these files now provide consistent numbering of residues in different PDB entries allowing easy comparison of structure models. The extended dictionary yields a more consistent, standardised metadata description without altering the core PDB information. This development enables up-to-date cross-reference information at the residue level resulting in better data interoperability, supporting improved data analysis and visualisation.

Towards a structurally resolved human protein interaction network

Cellular functions are governed by molecular machines that assemble through protein-protein interactions. Their atomic details are critical to studying their molecular mechanisms. However, fewer than 5% of hundreds of thousands of human protein interactions have been structurally characterized. Here we test the potential and limitations of recent progress in deep-learning methods using AlphaFold2 to predict structures for 65,484 human protein interactions. We show that experiments can orthogonally confirm higher-confidence models. We identify 3,137 high-confidence models, of which 1,371 have no homology to a known structure. We identify interface residues harboring disease mutations, suggesting potential mechanisms for pathogenic variants. Groups of interface phosphorylation sites show patterns of co-regulation across conditions, suggestive of coordinated tuning of multiple protein interactions as signaling responses. Finally, we provide examples of how the predicted binary complexes can be used to build larger assemblies helping to expand our understanding of human cell biology.

Computational systems approach towards phosphoproteomics and their downstream regulation

Protein phosphorylation plays an essential role in modulating cell signalling and its downstream transcriptional and translational regulations. Until recently, protein phosphorylation has been studied mostly using low-throughput biochemical assays. The advancement of mass spectrometry (MS)-based phosphoproteomics transformed the field by enabling measurement of proteome-wide phosphorylation events, where tens of thousands of phosphosites are routinely identified and quantified in an experiment. This has brought a significant challenge in analysing large-scale phosphoproteomic data, making computational methods and systems approaches integral parts of phosphoproteomics. Previous works have primarily focused on reviewing the experimental techniques in MS-based phosphoproteomics, yet a systematic survey of the computational landscape in this field is still missing. Here, we review computational methods and tools, and systems approaches that have been developed for phosphoproteomics data analysis. We categorise them into four aspects including data processing, functional analysis, phosphoproteome annotation and their integration with other omics, and in each aspect, we discuss the key methods and example studies. Lastly, we highlight some of the potential research directions on which future work would make a significant contribution to this fast-growing field. We hope this review provides a useful snapshot of the field of computational systems phosphoproteomics and stimulates new research that drives future development.

PTMint database of experimentally verified PTM regulation on protein-protein interaction

MOTIVATION

Post-translational modification (PTM) is an important biochemical process. which includes six most well-studied types: phosphorylation, acetylation, methylation, sumoylation, ubiquitylation and glycosylation. PTM is involved in various cell signaling pathways and biological processes. Abnormal PTM status is closely associated with severe diseases (such as cancer and neurologic diseases) by regulating protein functions, such as protein-protein interactions (PPIs). A set of databases was constructed separately for PTM sites and PPI; however, the resource of regulation for PTM on PPI is still unsolved.

RESULTS

Here, we firstly constructed a public accessible database of PTMint (PTMs that are associated with PPIs) (https://ptmint.sjtu.edu.cn/) that contains manually curated complete experimental evidence of the PTM regulation on PPIs in multiple organisms, including Homo sapiens, Arabidopsis thaliana, Caenorhabditis elegans, Drosophila melanogaster, Saccharomyces cerevisiae and Schizosaccharomyces pombe. Currently, the first version of PTMint encompassed 2477 non-redundant PTM sites in 1169 proteins affecting 2371 protein-protein pairs involving 357 diseases. Various annotations were systematically integrated, such as protein sequence, structure properties and protein complex analysis. PTMint database can help to insight into disease mechanism, disease diagnosis and drug discovery associated with PTM and PPI.

AVAILABILITY AND IMPLEMENTATION

PTMint is freely available at: https://ptmint.sjtu.edu.cn/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Where protein structure and cell diversity meet

Protein-protein interaction networks - interactomes - are charted with the hope to understand how phenotypes emerge and how they are altered in disease states. Early efforts to map interactomes have focused on the assembly of context agnostic, reference networks. However, recent studies have mapped interactomes across different cell lines and tissues, finding highly variable interactomes due to the rewiring of protein-protein interactions in different contexts. Increasing evidence points to significant links between protein structure and interactome diversity seen across cell types and tissues. We discuss how recent findings support the key role of alternative splicing and phosphorylation, two well-established regulators of protein structural and functional diversity, in defining cell type- and tissue-specific interactomes. Moreover, we show that intrinsically disordered protein regions are most favorably equipped to support interactome rewiring by acting as hubs of protein structure and function regulation.

The evolution of post-translational modifications

Post-translational modifications (PTMs) are chemical modifications that can regulate the activity and function of proteins. From an evolutionary perspective, they also represent a fast mechanism for the generation of phenotypic diversity and divergence. Advances in mass spectrometry have now enabled the identification of over 600 distinct PTM classes collectively spanning an order of 10⁶ unique sites. However, the chemical detection of PTMs has lagged far behind their functional characterisation, and relatively little is still known about the selective constraints that govern PTM evolution. In particular, the true fraction of PTM sites that are functional - and thus subject to selection - remains an open question. Here, I review advances made in the past two years towards understanding the evolution of PTMs and their associated enzymes.

Systematic Discovery of FBXW7-Binding Phosphodegrons Highlights Mitogen-Activated Protein Kinases as Important Regulators of Intracellular Protein Levels

A FBXW7 is an F-box E3 ubiquitin-ligase affecting cell growth by controlling protein degradation. Mechanistically, its effect on its substrates depends on the phosphorylation of degron motifs, but the abundance of these phosphodegrons has not been systematically explored. We used a ratiometric protein degradation assay geared towards the identification of FBXW7-binding degron motifs phosphorylated by mitogen-activated protein kinases (MAPKs). Most of the known FBXW7 targets are localized in the nucleus and function as transcription factors. Here, in addition to more transcription affecting factors (ETV5, KLF4, SP5, JAZF1, and ZMIZ1 CAMTA2), we identified phosphodegrons located in proteins involved in chromatin regulation (ARID4B, KMT2E, KMT2D, and KAT6B) or cytoskeletal regulation (MAP2, Myozenin-2, SMTL2, and AKAP11), and some other proteins with miscellaneous functions (EIF4G3, CDT1, and CCAR2). We show that the protein level of full-length ARID4B, ETV5, JAZF1, and ZMIZ1 are affected by different MAPKs since their FBXW7-mediated degradation was diminished in the presence of MAPK-specific inhibitors. Our results suggest that MAPK and FBXW7 partnership plays an important cellular role by directly affecting the level of key regulatory proteins. The data also suggest that the p38α-controlled phosphodegron in JAZF1 may be responsible for the pathological regulation of the cancer-related JAZF1-SUZ12 fusion construct implicated in endometrial stromal sarcoma.

A peptidoform based proteomic strategy for studying functions of post-translational modifications

Protein post-translational modifications (PTMs) generate an enormous, but as yet undetermined, expansion of the produced proteoforms. In this Viewpoint, we firstly reviewed the concepts of proteoform and peptidoform. We show that many of the current PTM biological investigation and annotation studies largely follow a PTM site-specific rather than proteoform-specific approach. We further illustrate a potentially useful matching strategy in which a particular "modified peptidoform" is matched to the corresponding "unmodified peptidoform" as a reference for the quantitative analysis between samples and conditions. We suggest this strategy has the potential to provide more directly relevant information to learn the PTM site-specific biological functions. Accordingly, we advocate for the wider use of the nomenclature "peptidoform" in future bottom-up proteomic studies.

Molecular dynamics shows complex interplay and long-range effects of post-translational modifications in yeast protein interactions

Post-translational modifications (PTMs) play a vital, yet often overlooked role in the living cells through modulation of protein properties, such as localization and affinity towards their interactors, thereby enabling quick adaptation to changing environmental conditions. We have previously benchmarked a computational framework for the prediction of PTMs’ effects on the stability of protein-protein interactions, which has molecular dynamics simulations followed by free energy calculations at its core. In the present work, we apply this framework to publicly available data on Saccharomyces cerevisiae protein structures and PTM sites, identified in both normal and stress conditions. We predict proteome-wide effects of acetylations and phosphorylations on protein-protein interactions and find that acetylations more frequently have locally stabilizing roles in protein interactions, while the opposite is true for phosphorylations. However, the overall impact of PTMs on protein-protein interactions is more complex than a simple sum of local changes caused by the introduction of PTMs and adds to our understanding of PTM cross-talk. We further use the obtained data to calculate the conformational changes brought about by PTMs. Finally, conservation of the analyzed PTM residues in orthologues shows that some predictions for yeast proteins will be mirrored to other organisms, including human. This work, therefore, contributes to our overall understanding of the modulation of the cellular protein interaction networks in yeast and beyond.

Proteins are a diverse set of biological molecules responsible for numerous functions within cells, such as obtaining energy from food or transport of small molecules, and many processes rely on interactions of specific proteins. Moreover, a single protein may acquire different roles depending on cellular requirements and as a response to changes in the environment. A commonly used way to quickly change protein’s function or activity is by introducing small chemical modifications on specific locations within the protein. These modifications can cause the protein to interact in a more or less stable way with other proteins. We have previously developed a computational pipeline for predicting the effect of modifications on interactions of proteins, and in this work we apply it to all yeast proteins with known structures. We find differences in effects on the binding for different types of modifications. Importantly, we demonstrate that the modifications far from the interaction interface also significantly contribute to binding due to their impact on protein’s shape, which is often neglected by other methods. This work contributes to our understanding of the modulation of protein interactions in yeast due to modifications, while our widely applicable method will allow similar investigations in other organisms.

Proteomic characterization of the human lens and Cataractogenesis

INTRODUCTION

The goal of this review is to highlight the triumphs and frontiers in measurement of the lens proteome as it relates to onset of age-related nuclear cataract. As global life expectancy increases, so too does the frequency of age-related nuclear cataracts. Molecular therapeutics do not exist for delay or relief of cataract onset in humans. Since lens fiber cells are incapable of protein synthesis after initial maturation, age-related changes in proteome composition and post-translational modification accumulation can be measured with various techniques. Several of these modifications have been associated with cataract onset.

AREAS COVERED

We discuss the impact of long-lived proteins on the lens proteome and lens homeostasis as well as proteomic techniques that may be used to measure proteomes at various levels of proteomic specificity and spatial resolution.

EXPERT OPINION

There is clear evidence that several proteome modifications are correlated with cataract formation. Past studies should be enhanced with cutting-edge, spatially resolved mass spectrometry techniques to enhance the specificity and sensitivity of modification detection as it relates to cataract formation.

In-depth and 3-dimensional exploration of the budding yeast phosphoproteome

Phosphorylation is one of the most dynamic and widespread post-translational modifications regulating virtually every aspect of eukaryotic cell biology. Here, we assemble a dataset from 75 independent phosphoproteomic experiments performed in our laboratory using Saccharomyces cerevisiae. We report 30,902 phosphosites identified from cells cultured in a range of DNA damage conditions and/or arrested in distinct cell cycle stages. To generate a comprehensive resource for the budding yeast community, we aggregate our dataset with the Saccharomyces Genome Database and another recently published study, resulting in over 46,000 budding yeast phosphosites. With the goal of enhancing the identification of functional phosphorylation events, we perform computational positioning of phosphorylation sites on available 3D protein structures and systematically identify events predicted to regulate protein complex architecture. Results reveal hundreds of phosphorylation sites mapping to or near protein interaction interfaces, many of which result in steric or electrostatic "clashes" predicted to disrupt the interaction. With the advancement of Cryo-EM and the increasing number of available structures, our approach should help drive the functional and spatial exploration of the phosphoproteome.

Systematic Identification of Protein Phosphorylation-Mediated Interactions

Protein phosphorylation is a key regulatory mechanism involved in nearly every eukaryotic cellular process. Increasingly sensitive mass spectrometry approaches have identified hundreds of thousands of phosphorylation sites, but the functions of a vast majority of these sites remain unknown, with fewer than 5% of sites currently assigned a function. To increase our understanding of functional protein phosphorylation we developed an approach (phospho-DIFFRAC) for identifying the phosphorylation-dependence of protein assemblies in a systematic manner. A combination of nonspecific protein phosphatase treatment, size-exclusion chromatography, and mass spectrometry allowed us to identify changes in protein interactions after the removal of phosphate modifications. With this approach we were able to identify 316 proteins involved in phosphorylation-sensitive interactions. We recovered known phosphorylation-dependent interactors such as the FACT complex and spliceosome, as well as identified novel interactions such as the tripeptidyl peptidase TPP2 and the supraspliceosome component ZRANB2. More generally, we find phosphorylation-dependent interactors to be strongly enriched for RNA-binding proteins, providing new insight into the role of phosphorylation in RNA binding. By searching directly for phosphorylated amino acid residues in mass spectrometry data, we identified the likely regulatory phosphosites on ZRANB2 and FACT complex subunit SSRP1. This study provides both a method and resource for obtaining a better understanding of the role of phosphorylation in native macromolecular assemblies. All mass spectrometry data are available through PRIDE (accession #PXD021422).

Global and Site-Specific Effect of Phosphorylation on Protein Turnover

To date, the effects of specific modification types and sites on protein lifetime have not been systematically illustrated. Here, we describe a proteomic method, DeltaSILAC, to quantitatively assess the impact of site-specific phosphorylation on the turnover of thousands of proteins in live cells. Based on the accurate and reproducible mass spectrometry-based method, a pulse labeling approach using stable isotope-labeled amino acids in cells (pSILAC), phosphoproteomics, and a unique peptide-level matching strategy, our DeltaSILAC profiling revealed a global, unexpected delaying effect of many phosphosites on protein turnover. We further found that phosphorylated sites accelerating protein turnover are functionally selected for cell fitness, enriched in Cyclin-dependent kinase substrates, and evolutionarily conserved, whereas the glutamic acids surrounding phosphosites significantly delay protein turnover. Our method represents a generalizable approach and provides a rich resource for prioritizing the effects of phosphorylation sites on protein lifetime in the context of cell signaling and disease biology.

Co-regulation of the transcription controlling ATF2 phosphoswitch by JNK and p38

Transcription factor phosphorylation at specific sites often activates gene expression, but how environmental cues quantitatively control transcription is not well-understood. Activating protein 1 transcription factors are phosphorylated by mitogen-activated protein kinases (MAPK) in their transactivation domains (TAD) at so-called phosphoswitches, which are a hallmark in response to growth factors, cytokines or stress. We show that the ATF2 TAD is controlled by functionally distinct signaling pathways (JNK and p38) through structurally different MAPK binding sites. Moreover, JNK mediated phosphorylation at an evolutionarily more recent site diminishes p38 binding and made the phosphoswitch differently sensitive to JNK and p38 in vertebrates. Structures of MAPK-TAD complexes and mechanistic modeling of ATF2 TAD phosphorylation in cells suggest that kinase binding motifs and phosphorylation sites line up to maximize MAPK based co-regulation. This study shows how the activity of an ancient transcription controlling phosphoswitch became dependent on the relative flux of upstream signals.

Mass Spectrometry-Based Characterization of the Virion Proteome, Phosphoproteome, and Associated Kinase Activity of Human Cytomegalovirus

The assembly of human cytomegalovirus (HCMV) virions is an orchestrated process that requires, as an essential prerequisite, the complex crosstalk between viral structural proteins. Currently, however, the mechanisms governing the successive steps in the constitution of virion protein complexes remain elusive. Protein phosphorylation is a key regulator determining the sequential changes in the conformation, binding, dynamics, and stability of proteins in the course of multiprotein assembly. In this review, we present a comprehensive map of the HCMV virion proteome, including a refined view on the virion phosphoproteome, based on previous publications supplemented by new results. Thus, a novel dataset of viral and cellular proteins contained in HCMV virions is generated, providing a basis for future analyses of individual phosphorylation steps and sites involved in the orchestrated assembly of HCMV virion-specific multiprotein complexes. Finally, we present the current knowledge on the activity of pUL97, the HCMV-encoded and virion-associated kinase, in phosphorylating viral and host proteins.

Calcineurin Inhibition Causes α2δ-1-Mediated Tonic Activation of Synaptic NMDA Receptors and Pain Hypersensitivity

Calcineurin inhibitors, such as tacrolimus (FK506) and cyclosporine, are widely used as standard immunosuppressants in organ transplantation recipients. However, these drugs can cause severe pain in patients, commonly referred to as calcineurin inhibitor-induced pain syndrome (CIPS). Although calcineurin inhibition increases NMDAR activity in the spinal cord, the underlying mechanism remains enigmatic. Using an animal model of CIPS, we found that systemic administration of FK506 in male and female mice significantly increased the amount of α2δ-1-GluN1 complexes in the spinal cord and the level of α2δ-1-bound GluN1 proteins in spinal synaptosomes. Treatment with FK506 significantly increased the frequency of mEPSCs and the amplitudes of monosynaptic EPSCs evoked from the dorsal root and puff NMDAR currents in spinal dorsal horn neurons. Inhibiting α2δ-1 with gabapentin or disrupting the α2δ-1-NMDAR interaction with α2δ-1Tat peptide completely reversed the effects of FK506. In α2δ-1 gene KO mice, treatment with FK506 failed to increase the frequency of NMDAR-mediated mEPSCs and the amplitudes of evoked EPSCs and puff NMDAR currents in spinal dorsal horn neurons. Furthermore, systemic administration of gabapentin or intrathecal injection of α2δ-1Tat peptide reversed thermal and mechanical hypersensitivity in FK506-treated mice. In addition, genetically deleting GluN1 in dorsal root ganglion neurons or α2δ-1 genetic KO similarly attenuated FK506-induced thermal and mechanical hypersensitivity. Together, our findings indicate that α2δ-1-bound NMDARs mediate calcineurin inhibitor-induced tonic activation of presynaptic and postsynaptic NMDARs at the spinal cord level and that presynaptic NMDARs play a prominent role in the development of CIPS.SIGNIFICANCE STATEMENT Calcineurin inhibitors are immunosuppressants used to prevent rejection of transplanted organs and tissues. However, these drugs can cause severe, unexplained pain. We showed that calcineurin inhibition enhances physical interaction between α2δ-1 and NMDARs and their synaptic trafficking in the spinal cord. α2δ-1 is essential for calcineurin inhibitor-induced aberrant activation of presynaptic and postsynaptic NMDARs in the spinal cord. Furthermore, inhibiting α2δ-1 or disrupting α2δ-1-NMDAR interaction reduces calcineurin inhibitor-induced pain hypersensitivity. Eliminating NMDARs in primary sensory neurons or α2δ-1 KO also attenuates calcineurin inhibitor-induced pain hypersensitivity. This new information extends our mechanistic understanding of the role of endogenous calcineurin in regulating synaptic plasticity and nociceptive transmission and suggests new strategies for treating this painful condition.

The functional landscape of the human phosphoproteome

Protein phosphorylation is a key post-translational modification regulating protein function in almost all cellular processes. Although tens of thousands of phosphorylation sites have been identified in human cells, approaches to determine the functional importance of each phosphosite are lacking. Here, we manually curated 112 datasets of phospho-enriched proteins, generated from 104 different human cell types or tissues. We re-analyzed the 6,801 proteomics experiments that passed our quality control criteria, creating a reference phosphoproteome containing 119,809 human phosphosites. To prioritize functional sites, we used machine learning to identify 59 features indicative of proteomic, structural, regulatory or evolutionary relevance and integrate them into a single functional score. Our approach identifies regulatory phosphosites across different molecular mechanisms, processes and diseases, and reveals genetic susceptibilities at a genomic scale. Several regulatory phosphosites were experimentally validated, including identifying a role in neuronal differentiation for phosphosites in SMARCC2, a member of the SWI/SNF chromatin-remodeling complex.

Genetics' Piece of the PI: Inferring the Origin of Complex Traits and Diseases from Proteome-Wide Protein-Protein Interaction Dynamics

How do common and rare genetic polymorphisms contribute to quantitative traits or disease risk and progression? Multiple human traits have been extensively characterized at the genomic level, revealing their complex genetic architecture. However, it is difficult to resolve the mechanisms by which specific variants contribute to a phenotype. Recently, analyses of variant effects on molecular traits have uncovered intermediate mechanisms that link sequence variation to phenotypic changes. Yet, these methods only capture a fraction of genetic contributions to phenotype. Here, in reviewing the field, it is proposed that complex traits can be understood by characterizing the dynamics of biochemical networks within living cells, and that the effects of genetic variation can be captured on these networks by using protein-protein interaction (PPI) methodologies. This synergy between PPI methodologies and the genetics of complex traits opens new avenues to investigate the molecular etiology of human diseases and to facilitate their prevention or treatment.

Mass spectrometric analysis of PTM dynamics using stable isotope labeled metabolic precursors in cell culture

Biological organisms represent highly dynamic systems, which are continually exposed to environmental factors and always strive to restore steady-state homeostasis. Posttranslational modifications are key regulators with which biological systems respond to external stimuli. To understand how homeostasis is restored, it is important to study the kinetics of posttranslational modifications. In this review we discuss proteomic approaches using stable isotope labeled metabolic precursors to study dynamics of posttranslational modifications in cell culture.

Spatiotemporal changes in the human lens proteome: Critical insights into long-lived proteins

The ocular lens is a unique tissue that contains an age gradient of cells and proteins ranging from newly differentiated cells containing newly synthesized proteins to cells and proteins that are as old as the organism. Thus, the ocular lens is an excellent model for studying long-lived proteins (LLPs) and the effects of aging and post-translational modifications on protein structure and function. Given the architecture of the lens, with young fiber cells in the outer cortex and the oldest cells in the lens nucleus, spatially-resolved studies provide information on age-specific protein changes. In this review, experimental strategies and proteomic methods that have been used to examine age-related and cataract-specific changes to the human lens proteome are described. Measured spatio-temporal changes in the human lens proteome are summarized and reveal a highly consistent, time-dependent set of modifications observed in transparent human lenses. Such measurements have led to the discovery of cataract-specific modifications and the realization that many animal systems are unsuitable to study many of these modifications. Mechanisms of protein modifications such as deamidation, racemization, truncation, and protein-protein crosslinking are presented and the implications of such mechanisms for other long-lived proteins in other tissues are discussed in the context of age-related neurological diseases. A comprehensive understanding of LLP modifications will enhance our ability to develop new therapies for the delay, prevention or reversal of age-related diseases.

Multi-omics Characterization of Interaction-mediated Control of Human Protein Abundance levels

Proteogenomic studies of cancer samples have shown that copy-number variation can be attenuated at the protein level for a large fraction of the proteome, likely due to the degradation of unassembled protein complex subunits. Such interaction-mediated control of protein abundance remains poorly characterized. To study this, we compiled genomic, (phospho)proteomic and structural data for hundreds of cancer samples and find that up to 42% of 8,124 analyzed proteins show signs of post-transcriptional control. We find evidence of interaction-dependent control of protein abundance, correlated with interface size, for 516 protein pairs, with some interactions further controlled by phosphorylation. Finally, these findings in cancer were reflected in variation in protein levels in normal tissues. Importantly, expression differences due to natural genetic variation were increasingly buffered from phenotype differences for highly attenuated proteins. Altogether, this study further highlights the importance of posttranscriptional control of protein abundance in cancer and healthy cells.

FunFam protein families improve residue level molecular function prediction

BACKGROUND

The CATH database provides a hierarchical classification of protein domain structures including a sub-classification of superfamilies into functional families (FunFams). We analyzed the similarity of binding site annotations in these FunFams and incorporated FunFams into the prediction of protein binding residues.

RESULTS

FunFam members agreed, on average, in 36.9 ± 0.6% of their binding residue annotations. This constituted a 6.7-fold increase over randomly grouped proteins and a 1.2-fold increase (1.1-fold on the same dataset) over proteins with the same enzymatic function (identical Enzyme Commission, EC, number). Mapping de novo binding residue prediction methods (BindPredict-CCS, BindPredict-CC) onto FunFam resulted in consensus predictions for those residues that were aligned and predicted alike (binding/non-binding) within a FunFam. This simple consensus increased the F1-score (for binding) 1.5-fold over the original prediction method. Variation of the threshold for how many proteins in the consensus prediction had to agree provided a convenient control of accuracy/precision and coverage/recall, e.g. reaching a precision as high as 60.8 ± 0.4% for a stringent threshold.

CONCLUSIONS

The FunFams outperformed even the carefully curated EC numbers in terms of agreement of binding site residues. Additionally, we assume that our proof-of-principle through the prediction of protein binding residues will be relevant for many other solutions profiting from FunFams to infer functional information at the residue level.

Regulation of a plant aquaporin by a Casparian strip membrane domain protein-like

The absorption of soil water by roots allows plants to maintain their water status. At the endodermis, water transport can be affected by initial formation of a Casparian strip and further deposition of suberin lamellas and regulated by the function of aquaporins. Four Casparian strip membrane domain protein-like (CASPL; CASPL1B1, CASPL1B2, CASPL1D1, and CASPL1D2) were previously shown to interact with PIP2;1. The present work shows that CASPL1B1, CASPL1B2, and CASPL1D2 are exclusively expressed in suberized endodermal cells, suggesting a cell-specific role in suberization and/or water transport regulation. When compared with wild-type plants, and by contrast to caspl1b1*caspl1b2 double loss of function, caspl1d1*caspl1d2 double mutants showed, in some control or NaCl stress experiments and not upon abscisic acid (ABA) treatment, a weak enlargement of the continuous suberization zone. None of the mutants showed root hydraulic conductivity (Lp_r ) phenotype, whether in control, NaCl, or ABA treatment conditions. The data suggest a slight negative role for CASPL1D1 and CASPL1D2 in suberization under control or salt stress conditions, with no major impact on whole root transport functions. At the molecular level, CASPL1B1 was able to physically interact with PIP2;1 and potentially could influence the regulation of aquaporins by acting on their phosphorylated form.

Visualizing Post-Translational Modifications in Protein Interaction Networks Using PTMOracle

Post-translational modifications (PTMs) of proteins act as key regulators of protein activity, including the regulation of protein-protein interactions (PPIs). However, exploring functional links between PTMs and PPIs can be difficult. PTMOracle is a Cytoscape app that facilitates the co-visualization and co-analysis of PTMs in the context of PPI networks. PTMOracle also allows extensive data to be integrated and co-analyzed, allowing the role of domains, motifs, and disordered regions to be considered. Here, we describe several PTMOracle protocols investigating complex PTM-associated relationships and their role in PPIs. This is assisted by OraclePainter for coloring proteins by the modifications present and visualizing these in the context of networks, by OracleTools for cross-matching PTMs with sequence feature for all nodes in the network, and by OracleResults for exploring specific proteins and visualizing their PTMs in the context of protein sequences. This unit aims to demonstrate how PTMOracle can be used to systematically explore network visualizations and generate testable hypotheses regarding the functional role of PTMs in PPIs, and how the results can be analyzed to better understand the regulatory role of PTMs in PPIs. © 2019 by John Wiley & Sons, Inc.

Effects of Acetylation and Phosphorylation on Subunit Interactions in Three Large Eukaryotic Complexes

Protein post-translational modifications (PTMs) have an indispensable role in living cells as they expand chemical diversity of the proteome, providing a fine regulatory layer that can govern protein-protein interactions in changing environmental conditions. Here we investigated the effects of acetylation and phosphorylation on the stability of subunit interactions in purified Saccharomyces cerevisiae complexes, namely exosome, RNA polymerase II and proteasome. We propose a computational framework that consists of conformational sampling of the complexes by molecular dynamics simulations, followed by Gibbs energy calculation by MM/GBSA. After benchmarking against published tools such as FoldX and Mechismo, we could apply the framework for the first time on large protein assemblies with the aim of predicting the effects of PTMs located on interfaces of subunits on binding stability. We discovered that acetylation predominantly contributes to subunits' interactions in a locally stabilizing manner, while phosphorylation shows the opposite effect. Even though the local binding contributions of PTMs may be predictable to an extent, the long range effects and overall impact on subunits' binding were only captured because of our dynamical approach. Employing the developed, widely applicable workflow on other large systems will shed more light on the roles of PTMs in protein complex formation.

Proteome-wide analysis of phospho-regulated PDZ domain interactions

A key function of reversible protein phosphorylation is to regulate protein-protein interactions, many of which involve short linear motifs (3-12 amino acids). Motif-based interactions are difficult to capture because of their often low-to-moderate affinities. Here, we describe phosphomimetic proteomic peptide-phage display, a powerful method for simultaneously finding motif-based interaction and pinpointing phosphorylation switches. We computationally designed an oligonucleotide library encoding human C-terminal peptides containing known or predicted Ser/Thr phosphosites and phosphomimetic variants thereof. We incorporated these oligonucleotides into a phage library and screened the PDZ (PSD-95/Dlg/ZO-1) domains of Scribble and DLG1 for interactions potentially enabled or disabled by ligand phosphorylation. We identified known and novel binders and characterized selected interactions through microscale thermophoresis, isothermal titration calorimetry, and NMR We uncover site-specific phospho-regulation of PDZ domain interactions, provide a structural framework for how PDZ domains accomplish phosphopeptide binding, and discuss ligand phosphorylation as a switching mechanism of PDZ domain interactions. The approach is readily scalable and can be used to explore the potential phospho-regulation of motif-based interactions on a large scale.

Mcl-1 Phosphorylation without Degradation Mediates Sensitivity to HDAC Inhibitors by Liberating BH3-Only Proteins

Mcl-1, a prosurvival Bcl-2 family protein, is frequently overexpressed in cancer cells and plays a critical role in therapeutic resistance. It is well known that anticancer agents induce phosphorylation of Mcl-1, which promotes its binding to E3 ubiquitin ligases and subsequent proteasomal degradation and apoptosis. However, other functions of Mcl-1 phosphorylation in cancer cell death have not been well characterized. In this study, we show in colon cancer cells that histone deacetylase inhibitors (HDACi) induce GSK3β-dependent Mcl-1 phosphorylation, but not degradation or downregulation. The in vitro and in vivo anticancer effects of HDACi were dependent on Mcl-1 phosphorylation and were blocked by genetic knock-in of a Mcl-1 phosphorylation site mutant. Phosphorylation-dead Mcl-1 maintained cell survival by binding and sequestering BH3-only Bcl-2 family proteins PUMA, Bim, and Noxa, which were upregulated and necessary for apoptosis induction by HDACi. Resistance to HDACi mediated by phosphorylation-dead Mcl-1 was reversed by small-molecule Mcl-1 inhibitors that liberated BH3-only proteins. These results demonstrate a critical role of Mcl-1 phosphorylation in mediating HDACi sensitivity through a novel and degradation-independent mechanism. These results provide new mechanistic insights on how Mcl-1 maintains cancer cell survival and suggest that Mcl-1-targeting agents are broadly useful for overcoming therapeutic resistance in cancer cells.Significance: These findings present a novel degradation-independent function of Mcl-1 phosphorylation in anticancer therapy that could be useful for developing new Mcl-1-targeting agents to overcome therapeutic resistance. Cancer Res; 78(16); 4704-15. ©2018 AACR.

The effect of phosphorylation on the salt-tolerance-related functions of the soybean protein PM18, a member of the group-3 LEA protein family

Enzymatically driven post-translated modifications (PTMs) usually happen within the intrinsically disordered regions of a target protein and can modulate variety of protein functions. Late embryogenesis abundant (LEA) proteins are a family of the plant intrinsically disordered proteins (IDPs). Despite their important roles in plant stress response, there is currently limited knowledge on the presence and functional and structural effects of phosphorylation on LEA proteins. In this study, we identified three phosphorylation sites (Ser90, Tyr136, and Thr266) in the soybean PM18 protein that belongs to the group-3 LEA proteins. In yeast expression system, PM18 protein increased the salt tolerance of yeast, and the phosphorylation of this protein further enhanced its protective function. Further analysis revealed that Ser90 and Tyr136 are more important than Thr266, and these two sites might work cooperatively in regulating the salt resistance function of PM18. The circular dichroism analysis showed that PM18 protein was disordered in aqueous media, and phosphorylation did not affect the disordered status of this protein. However, phosphorylation promoted formation of more helical structure in the presence of sodium dodecyl sulfate (SDS) or trifluoroethanol (TFE). Furthermore, in dedicated in vitro experiments, phosphorylated PM18 protein was able to better protect lactate dehydrogenase (LDH) from the inactivation induced by the freeze-thaw cycles than its un- or dephosphorylated forms. All these data indicate that phosphorylation may have regulatory effects on the stress-tolerance-related function of LEA proteins. Therefore, further studies are needed to shed more light on functional and structural roles of phosphorylation in LEA proteins.