Generative ß-hairpin design using a residue-based physicochemical property landscape.

De novo peptide design is a new frontier that has broad application potential in the biological and biomedical fields. Most existing models for de novo peptide design are largely based on sequence homology that can be restricted based on evolutionarily derived protein sequences and lack the physicochemical context essential in protein folding. Generative machine learning for de novo peptide design is a promising way to synthesize theoretical data that are based on, but unique from, the observable universe. In this study, we created and tested a custom peptide generative adversarial network intended to design peptide sequences that can fold into the ß-hairpin secondary structure. This deep neural network model is designed to establish a preliminary foundation of the generative approach based on physicochemical and conformational properties of 20 canonical amino acids, for example, hydrophobicity and residue volume, using extant structure-specific sequence data from the PDB. The beta generative adversarial network model robustly distinguishes secondary structures of ß hairpin from α helix and intrinsically disordered peptides with an accuracy of up to 96% and generates artificial ß-hairpin peptide sequences with minimum sequence identities around 31% and 50% when compared against the current NCBI PDB and nonredundant databases, respectively. These results highlight the potential of generative models specifically anchored by physicochemical and conformational property features of amino acids to expand the sequence-to-structure landscape of proteins beyond evolutionary limits.

Competition between inside-out unfolding and pathogenic aggregation in an amyloid-forming ß-propeller.

Studies of folded-to-misfolded transitions using model protein systems reveal a range of unfolding needed for exposure of amyloid-prone regions for subsequent fibrillization. Here, we probe the relationship between unfolding and aggregation for glaucoma-associated myocilin. Mutations within the olfactomedin domain of myocilin (OLF) cause a gain-of-function, namely cytotoxic intracellular aggregation, which hastens disease progression. Aggregation by wild-type OLF (OLFWT) competes with its chemical unfolding, but only below the threshold where OLF loses tertiary structure. Representative moderate (OLFD380A) and severe (OLFI499F) disease variants aggregate differently, with rates comparable to OLFWT in initial stages of unfolding, and variants adopt distinct partially folded structures seen along the OLFWT urea-unfolding pathway. Whether initiated with mutation or chemical perturbation, unfolding propagates outward to the propeller surface. In sum, for this large protein prone to amyloid formation, the requirement for a conformational change to promote amyloid fibrillization leads to direct competition between unfolding and aggregation.

Depletion assisted hemin affinity (DAsHA) proteomics reveals an expanded landscape of heme-binding proteins in the human proteome.

Heme b (iron protoporphyrin IX) plays important roles in biology as a metallocofactor and signaling molecule. However, the targets of heme signaling and the network of proteins that mediate the exchange of heme from sites of synthesis or uptake to heme dependent or regulated proteins are poorly understood. Herein, we describe a quantitative mass spectrometry (MS)-based chemoproteomics strategy to identify exchange labile hemoproteins in human embryonic kidney HEK293 cells that may be relevant to heme signaling and trafficking. The strategy involves depleting endogenous heme with the heme biosynthetic inhibitor succinylacetone (SA), leaving putative heme-binding proteins in their apo-state, followed by the capture of those proteins using hemin-agarose resin, and finally elution and identification by MS. By identifying only those proteins that interact with high specificity to hemin-agarose relative to control beaded agarose in an SA-dependent manner, we have expanded the number of proteins and ontologies that may be involved in binding and buffering labile heme or are targets of heme signaling. Notably, these include proteins involved in chromatin remodeling, DNA damage response, RNA splicing, cytoskeletal organization, and vesicular trafficking, many of which have been associated with heme through complementary studies published recently. Taken together, these results provide support for the emerging role of heme in an expanded set of cellular processes from genome integrity to protein trafficking and beyond.

Sod1 integrates oxygen availability to redox regulate NADPH production and the thiol redoxome.

Cu/Zn superoxide dismutase (Sod1) is a highly conserved and abundant antioxidant enzyme that detoxifies superoxide (O2•-) by catalyzing its conversion to dioxygen (O2) and hydrogen peroxide (H2O2). Using Saccharomyces cerevisiae and mammalian cells, we discovered that a major aspect of the antioxidant function of Sod1 is to integrate O2 availability to promote NADPH production. The mechanism involves Sod1-derived H2O2 oxidatively inactivating the glycolytic enzyme, GAPDH, which in turn reroutes carbohydrate flux to the oxidative phase of the pentose phosphate pathway (oxPPP) to generate NADPH. The aerobic oxidation of GAPDH is dependent on and rate-limited by Sod1. Thus, Sod1 senses O2 via O2•- to balance glycolytic and oxPPP flux, through control of GAPDH activity, for adaptation to life in air. Importantly, this mechanism for Sod1 antioxidant activity requires the bulk of cellular Sod1, unlike for its role in protection against O2•- toxicity, which only requires <1% of total Sod1. Using mass spectrometry, we identified proteome-wide targets of Sod1-dependent redox signaling, including numerous metabolic enzymes. Altogether, Sod1-derived H2O2 is important for antioxidant defense and a master regulator of metabolism and the thiol redoxome.

Combinatorial phosphorylation modulates the structure and function of the G protein γ subunit in yeast.

Intrinsically disordered regions (IDRs) in proteins are often targets of combinatorial posttranslational modifications, which serve to regulate protein structure and function. Emerging evidence suggests that the N-terminal tails of G protein γ subunits, which are essential components of heterotrimeric G proteins, are intrinsically disordered, phosphorylation-dependent determinants of G protein signaling. Here, we found that the yeast Gγ subunit Ste18 underwent combinatorial, multisite phosphorylation events within its N-terminal IDR. G protein-coupled receptor (GPCR) activation and osmotic stress induced phosphorylation at Ser7, whereas glucose and acid stress induced phosphorylation at Ser3, which was a quantitative indicator of intracellular pH. Each site was phosphorylated by a distinct set of kinases, and phosphorylation of one site affected phosphorylation of the other, as determined through exposure to serial stimuli and through phosphosite mutagenesis. Last, we showed that phosphorylation resulted in changes in IDR structure and that different combinations of phosphorylation events modulated the activation rate and amplitude of the downstream mitogen-activated protein kinase Fus3. These data place Gγ subunits among intrinsically disordered proteins that undergo combinatorial posttranslational modifications that govern signaling pathway output.

Different Grp94 components interact transiently with the myocilin olfactomedin domain in vitro to enhance or retard its amyloid aggregation.

The inherited form of open angle glaucoma arises due to a toxic gain-of-function intracellular misfolding event involving a mutated myocilin olfactomedin domain (OLF). Mutant myocilin is recognized by the endoplasmic reticulum (ER)-resident heat shock protein 90 paralog, glucose regulated protein 94 (Grp94), but their co-aggregation precludes mutant myocilin clearance by ER-associated degradation. When the Grp94-mutant myocilin interaction is abrogated by inhibitors or siRNA, mutant myocilin is efficiently degraded. Here we dissected Grp94 into component domains (N, NM, MC) to better understand the molecular factors governing its interaction with OLF. We show that the Grp94 N-terminal nucleotide-binding N domain is responsible for accelerating OLF aggregation in vitro. Upon inhibiting the isolated N domain pharmacologically or removing the Pre-N terminal 57 residues from full-length Grp94, OLF aggregation rates revert to those seen for OLF alone, but only pharmacological inhibition rescues co-aggregation. The Grp94-OLF interaction is below the detection limit of fluorescence polarization measurements, but chemical crosslinking paired with mass spectrometry analyses traps a reproducible interaction between OLF and the Grp94 N domain, as well as between OLF and the Grp94 M domain. The emerging molecular-level picture of quinary interactions between Grp94 and myocilin points to a role for the far N-terminal sequence of the Grp94 N domain and a cleft in the M domain. Our work further supports drug discovery efforts to inhibit these interactions as a strategy to treat myocilin-associated glaucoma.

Genetic heterogeneity within collective invasion packs drives leader and follower cell phenotypes.

Collective invasion, the coordinated movement of cohesive packs of cells, has become recognized as a major mode of metastasis for solid tumors. These packs are phenotypically heterogeneous and include specialized cells that lead the invasive pack and others that follow behind. To better understand how these unique cell types cooperate to facilitate collective invasion, we analyzed transcriptomic sequence variation between leader and follower populations isolated from the H1299 non-small cell lung cancer cell line using an image-guided selection technique. We now identify 14 expressed mutations that are selectively enriched in leader or follower cells, suggesting a novel link between genomic and phenotypic heterogeneity within a collectively invading tumor cell population. Functional characterization of two phenotype-specific candidate mutations showed that ARP3 enhances collective invasion by promoting the leader cell phenotype and that wild-type KDM5B suppresses chain-like cooperative behavior. These results demonstrate an important role for distinct genetic variants in establishing leader and follower phenotypes and highlight the necessity of maintaining a capacity for phenotypic plasticity during collective cancer invasion.

Heme bioavailability and signaling in response to stress in yeast cells.

Protoheme (hereafter referred to as heme) is an essential cellular cofactor and signaling molecule that is also potentially cytotoxic. To mitigate heme toxicity, heme synthesis and degradation are tightly coupled to heme utilization in order to limit the intracellular concentration of "free" heme. Such a model, however, would suggest that a readily accessible steady-state, bioavailable labile heme (LH) pool is not required for supporting heme-dependent processes. Using the yeast Saccharomyces cerevisiae as a model and fluorescent heme sensors, site-specific heme chelators, and molecular genetic approaches, we found here that 1) yeast cells preferentially use LH in heme-depleted conditions; 2) sequestration of cytosolic LH suppresses heme signaling; and 3) lead (Pb2+) stress contributes to a decrease in total heme, but an increase in LH, which correlates with increased heme signaling. We also observed that the proteasome is involved in the regulation of the LH pool and that loss of proteasomal activity sensitizes cells to Pb2+ effects on heme homeostasis. Overall, these findings suggest an important role for LH in supporting heme-dependent functions in yeast physiology.

Genetic Analysis of Rare Human Variants of Regulators of G Protein Signaling Proteins and Their Role in Human Physiology and Disease.

Regulators of G protein signaling (RGS) proteins modulate the physiologic actions of many neurotransmitters, hormones, and other signaling molecules. Human RGS proteins comprise a family of 20 canonical proteins that bind directly to G protein-coupled receptors/G protein complexes to limit the lifetime of their signaling events, which regulate all aspects of cell and organ physiology. Genetic variations account for diverse human traits and individual predispositions to disease. RGS proteins contribute to many complex polygenic human traits and pathologies such as hypertension, atherosclerosis, schizophrenia, depression, addiction, cancers, and many others. Recent analysis indicates that most human diseases are due to extremely rare genetic variants. In this study, we summarize physiologic roles for RGS proteins and links to human diseases/traits and report rare variants found within each human RGS protein exome sequence derived from global population studies. Each RGS sequence is analyzed using recently described bioinformatics and proteomic tools for measures of missense tolerance ratio paired with combined annotation-dependent depletion scores, and protein post-translational modification (PTM) alignment cluster analysis. We highlight selected variants within the well-studied RGS domain that likely disrupt RGS protein functions and provide comprehensive variant and PTM data for each RGS protein for future study. We propose that rare variants in functionally sensitive regions of RGS proteins confer profound change-of-function phenotypes that may contribute, in newly appreciated ways, to complex human diseases and/or traits. This information provides investigators with a valuable database to explore variation in RGS protein function, and for targeting RGS proteins as future therapeutic targets.

Negative Feedback Phosphorylation of Gγ Subunit Ste18 and the Ste5 Scaffold Synergistically Regulates MAPK Activation in Yeast.

Heterotrimeric G proteins (Gαßγ) are essential transducers in G protein signaling systems in all eukaryotes. In yeast, G protein signaling differentially activates mitogen-activated protein kinases (MAPKs)-Fus3 and Kss1-a phenomenon controlled by plasma membrane (PM) association of the scaffold protein Ste5. Here, we show that phosphorylation of the yeast Gγ subunit (Ste18), together with Fus3 docking on Ste5, controls the rate and stability of Ste5/PM association. Disruption of either element alone by point mutation has mild but reciprocal effects on MAPK activation. Disabling both elements results in ultra-fast and stable bulk Ste5/PM localization and Fus3 activation that is 6 times faster and 4 times more amplified compared to wild-type cells. These results further resolve the mechanism by which MAPK negative feedback phosphorylation controls pathway activation and provides compelling evidence that Gγ subunits can serve as intrinsic regulators of G protein signaling.

Tyrosine phosphorylation switching of a G protein.

Heterotrimeric G protein complexes are molecular switches relaying extracellular signals sensed by G protein-coupled receptors (GPCRs) to downstream targets in the cytoplasm, which effect cellular responses. In the plant heterotrimeric GTPase cycle, GTP hydrolysis, rather than nucleotide exchange, is the rate-limiting reaction and is accelerated by a receptor-like regulator of G signaling (RGS) protein. We hypothesized that posttranslational modification of the Gα subunit in the G protein complex regulates the RGS-dependent GTPase cycle. Our structural analyses identified an invariant phosphorylated tyrosine residue (Tyr166 in the Arabidopsis Gα subunit AtGPA1) located in the intramolecular domain interface where nucleotide binding and hydrolysis occur. We also identified a receptor-like kinase that phosphorylates AtGPA1 in a Tyr166-dependent manner. Discrete molecular dynamics simulations predicted that phosphorylated Tyr166 forms a salt bridge in this interface and potentially affects the RGS protein-accelerated GTPase cycle. Using a Tyr166 phosphomimetic substitution, we found that the cognate RGS protein binds more tightly to the GDP-bound Gα substrate, consequently reducing its ability to accelerate GTPase activity. In conclusion, we propose that phosphorylation of Tyr166 in AtGPA1 changes the binding pattern with AtRGS1 and thereby attenuates the steady-state rate of the GTPase cycle. We coin this newly identified mechanism "substrate phosphoswitching."

Both positional and chemical variables control in vitro proteolytic cleavage of a presenilin ortholog.

Mechanistic details of intramembrane aspartyl protease (IAP) chemistry, which is central to many biological and pathogenic processes, remain largely obscure. Here, we investigated the in vitro kinetics of a microbial intramembrane aspartyl protease (mIAP) fortuitously acting on the renin substrate angiotensinogen and the C-terminal transmembrane segment of amyloid precursor protein (C100), which is cleaved by the presenilin subunit of γ-secretase, an Alzheimer disease (AD)-associated IAP. mIAP variants with substitutions in active-site and putative substrate-gating residues generally exhibit impaired, but not abolished, activity toward angiotensinogen and retain the predominant cleavage site (His-Thr). The aromatic ring, but not the hydroxyl substituent, within Tyr of the catalytic Tyr-Asp (YD) motif plays a catalytic role, and the hydrolysis reaction incorporates bulk water as in soluble aspartyl proteases. mIAP hydrolyzes the transmembrane region of C100 at two major presenilin cleavage sites, one corresponding to the AD-associated Aß42 peptide (Ala-Thr) and the other to the non-pathogenic Aß48 (Thr-Leu). For the former site, we observed more favorable kinetics in lipid bilayer-mimicking bicelles than in detergent solution, indicating that substrate-lipid and substrate-enzyme interactions both contribute to catalytic rates. High-resolution MS analyses across four substrates support a preference for threonine at the scissile bond. However, results from threonine-scanning mutagenesis of angiotensinogen demonstrate a competing positional preference for cleavage. Our results indicate that IAP cleavage is controlled by both positional and chemical factors, opening up new avenues for selective IAP inhibition for therapeutic interventions.

Predicted Functional Implications of Phosphorylation of Regulator of G Protein Signaling Protein in Plants.

Heterotrimeric G proteins function in development, biotic, and abiotic stress responses, hormone signaling as well as sugar sensing. We previously proposed that discrimination of these various external signals in the G protein pathway is accomplished in plants by membrane-localized receptor-like kinases (RLKs) rather than G-protein-coupled receptors. Arabidopsis thaliana Regulator of G Signaling protein 1 (AtRGS1) modulates G protein activation and is phosphorylated by several RLKs and by WITH-NO-LYSINE kinases (WNKs). Here, a combination of in vitro kinase assays, mass spectrometry, and computational bioinformatics identified and functionally prioritized phosphorylation sites in AtRGS1. Phosphosites for two more RLKs (BRL3 and PEPR1) were identified and added to the AtRGS1 phosphorylation profile. Bioinformatics analyses revealed that RLKs and WNK kinases phosphorylate plant RGS proteins within regions that are conserved across eukaryotes and at a high frequency. Four phospho-sites among 14 identified are proximal to equivalent mammalian phosphosites that impact RGS function, including: pS437 and pT267 in GmRGS2, and pS339 and pS436 in AtRGS1. Based on these analyses, we propose that pS437 and pS436 regulate GmRGS2 and AtRGS1 protein interactions and/or localization, whereas pT267 is important for modulation of GmRGS2 GAP activity and localization. Moreover, pS339 most likely affects AtRGS1 activation.

Acylation of Superoxide Dismutase 1 (SOD1) at K122 Governs SOD1-Mediated Inhibition of Mitochondrial Respiration.

In this study, we employed proteomics to identify mechanisms of posttranslational regulation on cell survival signaling proteins. We focused on Cu-Zn superoxide dismutase (SOD1), which protects cells from oxidative stress. We found that acylation of K122 on SOD1, while not impacting SOD1 catalytic activity, suppressed the ability of SOD1 to inhibit mitochondrial metabolism at respiratory complex I. We found that deacylase depletion increased K122 acylation on SOD1, which blocked the suppression of respiration in a K122-dependent manner. In addition, we found that acyl-mimicking mutations at K122 decreased SOD1 accumulation in mitochondria, initially hinting that SOD1 may inhibit respiration directly within the intermembrane space (IMS). However, surprisingly, we found that forcing the K122 acyl mutants into the mitochondria with an IMS-targeting tag did not recover their ability to suppress respiration. Moreover, we found that suppressing or boosting respiration levels toggled SOD1 in or out of the mitochondria, respectively. These findings place SOD1-mediated inhibition of respiration upstream of its mitochondrial localization. Lastly, deletion-rescue experiments show that a respiration-defective mutant of SOD1 is also impaired in its ability to rescue cells from toxicity caused by SOD1 deletion. Together, these data suggest a previously unknown interplay between SOD1 acylation, metabolic regulation, and SOD1-mediated cell survival.

DRILL: An Electrospray Ionization-Mass Spectrometry Interface for Improved Sensitivity via Inertial Droplet Sorting and Electrohydrodynamic Focusing in a Swirling Flow.

We describe the DRILL (dry ion localization and locomotion) device, which is an interface for electrospray ionization (ESI)-mass spectrometry (MS) that exploits a swirling flow to enable the use of inertial separation to prescribe different fates for electrosprayed droplets based on their size. This source adds a new approach to charged droplet trajectory manipulation which, when combined with hydrodynamic drag forces and electric field forces, provides a rich range of possible DRILL operational modes. Here, we experimentally demonstrate sensitivity improvement obtained via vortex-induced inertial sorting of electrosprayed droplets/ions: one possible mode of DRILL operation. In this mode, DRILL removes larger droplets while accelerating the remainder of the ESI plume, producing a high velocity stream of gas-enriched spray with small, highly charged droplets and ions and directing it toward the MS inlet. The improved signal-to-noise ratio (10-fold enhancement) in the detection of angiotensin I is demonstrated using the DRILL interface coupled to ESI-MS along with an improved limit of detection (10-fold enhancement, 100 picomole) in the detection of angiotensin II. The utility of DRILL has also been demonstrated by liquid chromatography (LC)-MS: a stable isotope labeled peptide cocktail was spiked into a complex native tissue extract and quantified by unscheduled multiple reaction monitoring on a TSQ Vantage. DRILL demonstrated improved signal strength (up to a 700-fold) for 8 out of 9 peptides and had no effects on the peak shape of the transitions.

Systematic analysis of non-structural protein features for the prediction of PTM function potential by artificial neural networks.

Post-translational modifications (PTMs) provide an extensible framework for regulation of protein behavior beyond the diversity represented within the genome alone. While the rate of identification of PTMs has rapidly increased in recent years, our knowledge of PTM functionality encompasses less than 5% of this data. We previously developed SAPH-ire (Structural Analysis of PTM Hotspots) for the prioritization of eukaryotic PTMs based on function potential of discrete modified alignment positions (MAPs) in a set of 8 protein families. A proteome-wide expansion of the dataset to all families of PTM-bearing, eukaryotic proteins with a representational crystal structure and the application of artificial neural network (ANN) models demonstrated the broader applicability of this approach. Although structural features of proteins have been repeatedly demonstrated to be predictive of PTM functionality, the availability of adequately resolved 3D structures in the Protein Data Bank (PDB) limits the scope of these methods. In order to bridge this gap and capture the larger set of PTM-bearing proteins without an available, homologous structure, we explored all available MAP features as ANN inputs to identify predictive models that do not rely on 3D protein structural data. This systematic, algorithmic approach explores 8 available input features in exhaustive combinations (247 models; size 2-8). To control for potential bias in random sampling for holdback in training sets, we iterated each model across 100 randomized, sample training and testing sets-yielding 24,700 individual ANNs. The size of the analyzed dataset and iterative generation of ANNs represents the largest and most thorough investigation of predictive models for PTM functionality to date. Comparison of input layer combinations allows us to quantify ANN performance with a high degree of confidence and subsequently select a top-ranked, robust fit model which highlights 3,687 MAPs, including 10,933 PTMs with a high probability of biological impact but without a currently known functional role.

Proteome-wide Structural Analysis of PTM Hotspots Reveals Regulatory Elements Predicted to Impact Biological Function and Disease.

Post-translational modifications (PTMs) regulate protein behavior through modulation of protein-protein interactions, enzymatic activity, and protein stability essential in the translation of genotype to phenotype in eukaryotes. Currently, less than 4% of all eukaryotic PTMs are reported to have biological function - a statistic that continues to decrease with an increasing rate of PTM detection. Previously, we developed SAPH-ire (Structural Analysis of PTM Hotspots) - a method for the prioritization of PTM function potential that has been used effectively to reveal novel PTM regulatory elements in discrete protein families (Dewhurst et al., 2015). Here, we apply SAPH-ire to the set of eukaryotic protein families containing experimental PTM and 3D structure data - capturing 1,325 protein families with 50,839 unique PTM sites organized into 31,747 modified alignment positions (MAPs), of which 2010 (∼6%) possess known biological function. Here, we show that using an artificial neural network model (SAPH-ire NN) trained to identify MAP hotspots with biological function results in prediction outcomes that far surpass the use of single hotspot features, including nearest neighbor PTM clustering methods. We find the greatest enhancement in prediction for positions with PTM counts of five or less, which represent 98% of all MAPs in the eukaryotic proteome and 90% of all MAPs found to have biological function. Analysis of the top 1092 MAP hotspots revealed 267 of truly unknown function (containing 5443 distinct PTMs). Of these, 165 hotspots could be mapped to human KEGG pathways for normal and/or disease physiology. Many high-ranking hotspots were also found to be disease-associated pathogenic sites of amino acid substitution despite the lack of observable PTM in the human protein family member. Taken together, these experiments demonstrate that the functional relevance of a PTM can be predicted very effectively by neural network models, revealing a large but testable body of potential regulatory elements that impact hundreds of different biological processes important in eukaryotic biology and human health.

Catalytic Properties of Intramembrane Aspartyl Protease Substrate Hydrolysis Evaluated Using a FRET Peptide Cleavage Assay.

Chemical details of intramembrane proteolysis remain elusive despite its prevalence throughout biology. We developed a FRET peptide assay for the intramembrane aspartyl protease (IAP) from Methanoculleus marisnigri JR1 in combination with quantitative mass spectrometry cleavage site analysis. IAP can hydrolyze the angiotensinogen sequence, a substrate for the soluble aspartyl protease renin, at a predominant cut site, His-Thr. Turnover is slow (min(-1) × 10(-3)), affinity and Michaelis constant (Km) values are in the low micromolar range, and both catalytic rates and cleavage sites are the same in detergent as reconstituted into bicelles. Three well-established, IAP-directed inhibitors were directly confirmed as competitive, albeit with modest inhibitor constant (Ki) values. Partial deletion of the first transmembrane helix results in a biophysically similar but less active enzyme than full-length IAP, indicating a catalytic role. Our study demonstrates previously unappreciated similarities with soluble aspartyl proteases, provides new biochemical features of IAP and inhibitors, and offers tools to study other intramembrane protease family members in molecular detail.

Structural Analysis of PTM Hotspots (SAPH-ire)--A Quantitative Informatics Method Enabling the Discovery of Novel Regulatory Elements in Protein Families.

Predicting the biological function potential of post-translational modifications (PTMs) is becoming increasingly important in light of the exponential increase in available PTM data from high-throughput proteomics. We developed structural analysis of PTM hotspots (SAPH-ire)--a quantitative PTM ranking method that integrates experimental PTM observations, sequence conservation, protein structure, and interaction data to allow rank order comparisons within or between protein families. Here, we applied SAPH-ire to the study of PTMs in diverse G protein families, a conserved and ubiquitous class of proteins essential for maintenance of intracellular structure (tubulins) and signal transduction (large and small Ras-like G proteins). A total of 1728 experimentally verified PTMs from eight unique G protein families were clustered into 451 unique hotspots, 51 of which have a known and cited biological function or response. Using customized software, the hotspots were analyzed in the context of 598 unique protein structures. By comparing distributions of hotspots with known versus unknown function, we show that SAPH-ire analysis is predictive for PTM biological function. Notably, SAPH-ire revealed high-ranking hotspots for which a functional impact has not yet been determined, including phosphorylation hotspots in the N-terminal tails of G protein gamma subunits--conserved protein structures never before reported as regulators of G protein coupled receptor signaling. To validate this prediction we used the yeast model system for G protein coupled receptor signaling, revealing that gamma subunit-N-terminal tail phosphorylation is activated in response to G protein coupled receptor stimulation and regulates protein stability in vivo. These results demonstrate the utility of integrating protein structural and sequence features into PTM prioritization schemes that can improve the analysis and functional power of modification-specific proteomics data.

Guanine nucleotide-binding protein (Gα) endocytosis by a cascade of ubiquitin binding domain proteins is required for sustained morphogenesis and proper mating in yeast.

Heterotrimeric G proteins are well known to transmit signals from cell surface receptors to intracellular effector proteins. There is growing appreciation that G proteins are also present at endomembrane compartments, where they can potentially interact with a distinct set of signaling proteins. Here, we examine the cellular trafficking function of the G protein α subunit in yeast, Gpa1. Gpa1 contains a unique 109-amino acid insert within the α-helical domain that undergoes a variety of posttranslational modifications. Among these is monoubiquitination, catalyzed by the NEDD4 family ubiquitin ligase Rsp5. Using a newly optimized method for G protein purification together with biophysical measures of structure and function, we show that the ubiquitination domain does not influence enzyme activity. By screening a panel of 39 gene deletion mutants, each lacking a different ubiquitin binding domain protein, we identify seven that are necessary to deliver Gpa1 to the vacuole compartment including four proteins (Ede1, Bul1, Ddi1, and Rup1) previously not known to be involved in this process. Finally, we show that proper endocytosis of the G protein is needed for sustained cellular morphogenesis and mating in response to pheromone stimulation. We conclude that a cascade of ubiquitin-binding proteins serves to deliver the G protein to its final destination within the cell. In this instance and in contrast to the previously characterized visual system, endocytosis from the plasma membrane is needed for proper signal transduction rather than for signal desensitization.