Emerging Roles of DHHC-mediated Protein S-palmitoylation in Physiological and Pathophysiological Context

Protein S-palmitoylation refers to a post-translational modification (PTM) wherein palmitic acid, a 16-carbon long saturated fatty acid gets covalently attached to Cys sidechain of a protein. It has been known to the literature for almost 50 years and in general, this PTM is believed to facilitate membrane attachments of proteins for the obvious hydrophobicity of the palmitoyl group. But after the discovery of the protein palmitoyl acyltransferases (PATs, also known as DHHC-PATs), a major paradigm shift has been observed in the field of protein S-palmitoylation. A family of 23 mammalian DHHC-PATs has been identified and the majority of them are associated with many human diseases spanning from neuropsychiatric diseases to cancers. Novel unique and essential role of DHHC-mediated protein S-palmitoylation has been revealed apart from its membrane trafficking role. Biomedical importance of DHHCs has also been reiterated with small molecule inhibitors for DHHCs as well as in DHHC-knockout mice or mouse Xenograft models. In this review, we present recent advances in the field of protein S-palmitoylation and the involvement of individual DHHC isoforms in human diseases. In addition, the recent development of the analytical tools to study S-palmitoylation and their inhibitors are discussed in detail. We also highlight the issues that need to be addressed in detail to further develop our understanding on protein S-palmitoylation and strongly believe that pharmacological modulation of DHHC-mediated protein S-palmitoylation has a massive potential to emerge as a novel therapeutic strategy for human diseases. It will not be surprising if reversible protein S-palmitoylation prove to be an indispensable PTM that regulates a host of cellular processes, just like protein phosphorylation or ubiquitination.     

A recent advance in the intracellular and extracellular redox post‐translational modification of proteins in plants

Plants, as sessile organisms, have acquired through evolution sophisticated regulatory signal pathways to overcome external variable factors during each stage of the life cycle. Among these regulatory signals, two pathways in particular, reactive oxygen species and reactive nitrogen species, have become of significant interest in several aspects of plant biology, underpinning these molecules as critical regulators during development, cellular differentiation, and plant‐pathogen interaction. Recently, redox posttranslational modifications (PTM), such as S‐nitrosylation on cysteine residues and tyrosine nitration, have shed light on multiple protein targets, as they are associated with signal networks/downstream metabolic pathways, capable of transducing the imbalance of redox hemostasis and consequently redirecting the biochemical status under stress conditions. However, most of the redox PTM have been studied only in the intracellular compartment, providing limited information concerning redox PTM in the extracellular matrix of plant cells. Nevertheless, recent studies have indicated the plausibility of redox PTM in extracellular proteins, including cell wall associated proteins. Accordingly, in this review, we endeavor to examine evidence of redox PTM supported by mass spectrometry data in the intracellular and extracellular space in plant cells. As a further example, we focus the last section of this review on illustrating, using molecular dynamics simulation, the effect of S‐nitrosylation on the structural conformation of well‐known cell wall‐associated proteins including pectin methylesterase and xyloglucan endo‐transglycosylases.     

ESNOQ,Proteomic Quantification of Endogenous S-Nitrosation

S-nitrosation is a post-translational protein modification and is one of the most important mechanisms of NO signaling. Endogenous S-nitrosothiol (SNO) quantification is a challenge for detailed functional studies. Here we developed an ESNOQ (Endogenous SNO Quantification) method which combines the stable isotope labeling by amino acids in cell culture (SILAC) technique with the detergent-free biotin-switch assay and LC-MS/MS. After confirming the accuracy of quantification in this method, we obtained an endogenous S-nitrosation proteome for LPS/IFN-γ induced RAW264.7 cells. 27 S-nitrosated protein targets were confirmed and using our method we were able to obtain quantitative information on the level of S-nitrosation on each modified Cys. With this quantitative information, over 15 more S-nitrosated targets were identified than in previous studies. Based on the quantification results, we found that the S-nitrosation levels of different cysteines varied within one protein, providing direct evidence for differences in the sensitivity of cysteine residues to reactive nitrosative stress and that S-nitrosation is a site-specific modification. Gene ontology clustering shows that S-nitrosation targets in the LPS/IFN-γ induced RAW264.7 cell model were functionally enriched in protein translation and glycolysis, suggesting that S-nitrosation may function by regulating multiple pathways. The ESNOQ method described here thus provides a solution for quantification of multiple endogenous S-nitrosation events, and makes it possible to elucidate the network of relationships between endogenous S-nitrosation targets involved in different cellular processes.     

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.Post-translational modifications (PTMs)¹ are a rapidly expanding and important class of protein feature that broaden the functional diversity of proteins in a proteome. By definition, PTMs change protein structure and therefore have the potential to affect protein function by altering protein interactions, protein stability or catalytic activity (1, 2). As they have been found to occur on nearly every protein in the eukaryotic proteome, PTMs broadly impact nearly all known cellular processes. Over 300 different types of PTM are known, ranging from single atom modifications (e.g. oxide) to small protein modifiers (e.g. ubiquitin), which can occur on all but five amino acid residues resulting from enzymatic or nonenzymatic processes (3). Over 220,000 distinct PTM sites have been experimentally identified across ∼77,000 different proteins to date (dbPTM; http://dbptm.mbc.nctu.edu.tw/statistics.php) – numbers that continue to grow exponentially because of improved methods for high throughput detection by mass spectrometry (MS). By virtue of how they are detected, most PTM data are sequence-linked and lack structural context.The function of most PTMs is unknown because the rate of PTM detection far surpasses the rate at which any one modification can be studied empirically. Moreover, the functional impact of every PTM is likely not equivalent (4). For example, computational analysis of phosphorylation sites in yeast and human proteomes indicate that well-conserved phosphosites are more likely to have a functional consequence compared with poorly conserved sites, yet only a fraction of phosphosites are well conserved (5, 6). Consequently, the development of tools that provide functional prioritization of PTMs could have a broad impact on our understanding of protein regulation, biological mechanism, and molecular evolution.The emerging need for methods that predict the functional impact of a PTM has not yet been met. Longstanding methods capitalize predominantly on the sequence context of PTMs and have been used to predict sites of modification (expasy.org/proteomics/post-translational_modification) and to compare enzyme/substrate interactions (7–9). More recently, studies aimed at expanding the parameters associated with functional PTMs have emerged. In these cases, a set of common features correlated with functional importance are derived from the analysis of PTMs within and between organisms including: number of PTM observations at a multiple sequence alignment position (i.e. hotspots), measures of co-occurrence between different PTMs (e.g. distance between phosphorylation and ubiquitination sites), biological dynamics (up or down-regulation), and protein–protein interaction influence (7, 10–12). Recent efforts to provide structural context by linking individual PTMs to three-dimensional structures in the protein data bank (PDB) have also been described (13, 14). However, these resources are extensions of existing PTM databases that allow visualization of single instances of modification onto individual proteins, but do not provide quantitative or analytical value.In principle, combining PTM hotspot and structural analysis would offer multiple advantages over any one approach used in isolation. Sequence homology provides protein family membership—thereby clustering PTMs into hotspots for groups of proteins to provide information about: (1) the evolutionary conservation and (2) observation frequencies of PTMs within the family. A primary consequence of their sequence homology is that members of a protein family will exhibit similar structures and protein interactions—features that dictate the function of protein systems. A secondary consequence is that PTM hotspots generated by alignment can be projected onto family-representative protein structures, which places each PTM hotspot into a three-dimensional context that can be visualized for each family. The structural context enabled by this projection can also provide spatial information about the PTM site that can supplement the sequence characteristics of the hotspot, namely: (3) solvent accessibility, which provides an estimate of whether a modification could occur on the folded protein; and (4) protein interface residence, which indicates the potential of the PTM to disrupt protein–protein interactions. Despite the theoretical advantages, no single tool has been developed that exploits the quantitative output from both sequence and structural data to evaluate the function potential of PTMs.Here we describe a new analytical method – Structural Analysis of PTM Hotspots (SAPH-ire), which ranks PTM hotspots by their potential to impact biological function for distinct protein families (Fig. 1). We demonstrate the application of SAPH-ire to the complete set of PTMs for eight distinct protein families including large heterotrimeric G proteins—revealing high-ranking hotspots for which a biological function has not yet been determined. In particular, SAPH-ire revealed the N-terminal tail (Nt) of G protein gamma (Gγ) subunits as one of the highest ranking PTM hotspots for heterotrimeric G proteins (Gα, Gβ, and Gγ). We tested this prediction by monitoring the phosphorylation state and mutation effects of phosphorylation sites in the N terminus of the yeast Gγ subunit (Ste18). Consistent with SAPH-ire predictions, we found that phosphorylation of Ste18-Nt is biologically responsive to a GPCR stimulus and that phospho-null or phospho-mimic mutation of these sites controls protein abundance in an opposite manner in vivo. Thus, SAPH-ire is a powerful new method for predicting the function potential of PTM hotspots, which can guide empirical research toward the discovery of new protein regulatory elements based on high-throughput proteomics.Open in a separate windowFig. 1.Schematic diagram of the SAPH-ire method. A, SAPH-ire integrates InterPro, the Protein Data bank (PDB) and a customized database of experimentally validated PTMs. Uniprot entries with PTMs that belong to specific InterPro-classified protein families undergo multiple-sequence alignment (MSA) and PTM hotspot analysis (HSA), which layers all PTMs for a given alignment position in the MSA. The total PTMs observed in each hotspot and the conservation of a modifiable residue (e.g. conservation lysine at a ubiquitination hotspot) at the hotspot are quantified. B, PTM hotspots within the protein family are then projected onto all known crystal structures for the family using the Structural Projection of PTMs (SPoP) tool. From the structural topology of PTM hotspots generated by SPoP, the solvent accessible surface area (SASA) and protein interface residence is quantified for each hotspot. C, PTM Function Potential Calculator (FPC) integrates the output from HSA and SPoP, resulting in PTM function potential scores for each hotspot. The function potential score can be used to rank PTM hotspots within or between protein families – prioritizing hotspots with the greatest potential to be biologically regulated and/or effect a biological function for the protein family of interest.     

Cross-species comparison of site-specific evolutionary-rate variation in influenza haemagglutinin

We investigate the causes of site-specific evolutionary-rate variation in influenza haemagglutinin (HA) between human and avian influenza, for subtypes H1, H3, and H5. By calculating the evolutionary-rate ratio, ω = dN/dS as a function of a residue''s solvent accessibility in the three-dimensional protein structure, we show that solvent accessibility has a significant but relatively modest effect on site-specific rate variation. By comparing rates within HA subtypes among host species, we derive an upper limit to the amount of variation that can be explained by structural constraints of any kind. Protein structure explains only 20–40% of the variation in ω. Finally, by comparing ω at sites near the sialic-acid-binding region to ω at other sites, we show that ω near the sialic-acid-binding region is significantly elevated in both human and avian influenza, with the exception of avian H5. We conclude that protein structure, HA subtype, and host biology all impose distinct selection pressures on sites in influenza HA.     

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.Since the discovery of phosphorylation in 1954 (1), post-translational modifications (PTMs)¹ have emerged as a broad class of protein feature that expand the functional proteome in eukaryotes. Improvements in the detection of PTMs by mass spectrometry have resulted in an exponential increase in our knowledge of the number and type of PTMs that make up the landscape of a modified eukaryotic proteome. As a result, the rate at which PTMs are discovered now far surpasses the rate at which they can be experimentally tested for biological function - a characteristic that is specific for each PTM and likely not equivalent between all PTMs that have been observed (2–4). Thus, effective methods of prioritization are essential for quantifying the likelihood of a site to be regulatory and/or impactful on biological function, which we refer to as the function potential of a PTM.Several unique features have been identified as predictors of biological impact for any given PTM - the determination of which relies on placing each PTM in the context of a multiple sequence alignment for a discrete protein or domain family, which we refer to as a Modified Alignment Position (MAP). For example, MAPs that are evolutionarily well conserved are more likely to exhibit biological function (3, 4). Similarly, functional PTMs are more commonly found within MAPs that exhibit a higher PTM observation frequency, are dynamic with respect to biological condition, located at protein interaction interfaces, and more solvent-accessible within a folded protein structure (5–7). Although efforts to elucidate the features associated with functional PTMs are relatively longstanding, few if any have established an integrative approach to quantitatively prioritize the function potential of PTMs beyond the use of single features.Previous evidence from our lab first demonstrated that multiple feature integration can improve functional prioritization. To accomplish this, we built Structural Analysis of PTM Hotspots (SAPH-ire)—an algorithm through which multiple predictors of PTM function are integrated to produce a single, quantitative function potential (FP) score that rank orders each hotspot within or between protein families (6) (Fig. 1). Previously, we used SAPH-ire to predict novel PTM regulatory elements in G protein families—including heterotrimeric G proteins—for which we discovered and experimentally confirmed a novel PTM regulatory element that is critical for cell signaling (6, 8). We propose that similar analysis of PTMs across the entire eukaryotic proteome is likely to result in the discovery of several novel regulatory elements that have yet to be realized.Open in a separate windowFig. 1.Schematic diagram of SAPH-ire. A, A theoretical segment of the multiple sequence alignment for a protein family (IPR000276; G protein-coupled receptor, rhodopsin-like) used here for illustrating the concept of SAPH-ire. Circled amino acid residues represent PTM sites experimentally observed on respective protein family members. Circle and arrow color represents the PTM observation frequency at each aligned position, called a MAP (modified alignment position), where green indicates 1 observation, blue for 2, orange for 3, and red for 5 or more. B, Cartoon rendering of bovine rhodopsin ({"type":"entrez-protein","attrs":{"text":"P02699","term_id":"129204","term_text":"P02699"}}P02699, RHO; PDB 2PED, chain A) showing side chains with projected PTM hotspots colored according to the number of observations within the family at each position aligned with the structural sequence. PDB coordinate data from the structurally projected PTM hotspots is used for calculation of solvent accessible surface area (SASA) and determination of protein interface residence (PPI). C, Hotspot features derived from the sequence and structural data are extracted for each protein family, where each hotspot corresponds to a precise family alignment position containing at least one PTM observation. D, Comparison of the comprehensive and SAPH-ire datasets representing all known experimental PTM data versus PTM data included in this study, respectively. E, Values calculated and derived from extracted hotspot features are analyzed by logistic regression or neural network models to produce probability scores for each hotspot.Here we apply SAPH-ire to protein families for which PTMs and protein structure are currently available, resulting in function potential prediction for 50,839 experimental PTM sites distributed across 31,747 MAPs. Using a neural network model (SAPH-ire NN) trained to predict the identity of embedded known-function MAPs, we derived a probability score that allows rank ordering for the likelihood of function for all MAPs including those with unknown function. We show that the SAPH-ire NN model significantly outperforms all other single or multi-feature predictive models and exhibits a proportional increase in predictive power for known function hotspots that have been more frequently studied (and therefore published). Using a strictly conservative probability threshold, we characterized the top-ranked 1092 MAPs corresponding to “function potential hotspots,” revealing 267 with truly unknown function - a striking fraction of which are also found mutated in human disease irrespective of whether the human protein, specifically, contains an observed PTM.     

Roles of sirtuins in the regulation of antioxidant defense and bioenergetic function of mitochondria under oxidative stress

Abstract

In addition to serving as the power house of mammalian cells, mitochondria are crucial for the maintenance of cellular homeostasis in response to physiological or environmental changes. Several lines of evidence suggest that posttranslational modification (PTM) of proteins plays a pivotal role in the regulation of the bioenergetic function of mitochondria. Among them, reversible lysine acetylation of mitochondrial proteins has been established as one of the key mechanisms in cellular response to energy demand by modulating the flux of a number of key metabolic pathways. In this article, we focus on the role of Sirt3-mediated deacetylation in: (1) flexibility of energy metabolism, (2) activation of antioxidant defense, and (3) maintenance of cellular redox status in response to dietary challenge and oxidative stress. We suggest that oxidative stress-elicited down-regulation of Sirt3 plays a role in the pathophysiology of diabetes, cardiac hypotrophy, mitochondrial diseases, and age-related diseases. Besides, the physiological role of newly identified lysine acylation mediated by Sirt5 and its biochemical effects on oxidative metabolism are also discussed. Moreover, we have integrated the regulatory function of several protein kinases that are involved in the phosphorylation of mitochondrial enzymes during oxidative stress. Finally, the functional consequence of the synergistic regulation through diverse protein modifications is emphasized on the maintenance of the bioenergetic homeostasis and metabolic adaptation of the animal and human cells. Together, we have provided an updated review of PTM in mitochondrial biology and their implications in aging and human diseases through an intricate regulation of energy metabolism under oxidative stress.     

Chemical modifications of respiratory complex I for structural and functional studies

Studies on chemical modifications of bacterial and mitochondrial complex I by synthetic chemical probes as well as endogenous chemicals have provided useful information on the structural and functional aspects of this enzyme. We herein reviewed recent studies that investigated chemical modifications of complex I by endogenous chemicals (e.g. Cys-S-nitrosation, Cys-S-glutathionylation, and Ser-O-phosphorylation) and synthetic reagents (e.g. Cys-SH modification by SH-reagents and the cross-linking of nearby subunits by bifunctional cross-linkers). We also reviewed recent photoaffinity labeling studies using complex I inhibitors, which can be recognized as “site-specific modification” by synthetic chemicals. In addition, we discussed the possibility of site-specific modification by various functional probes via ligand-directed tosylate (LDT) chemistry as a promising approach for unique biophysical studies on complex I.     

Identification of the substrate recognition region in the Δ6-fatty acid and Δ8-sphingolipid desaturase by fusion mutagenesis

Δ⁸-sphingolipid desaturase and Δ⁶-fatty acid desaturase share high protein sequence identity. Thus, it has been hypothesized that Δ⁶-fatty acid desaturase is derived from Δ⁸-sphingolipid desaturase; however, there is no direct proof. The substrate recognition regions of Δ⁶-fatty acid desaturase and Δ⁸-sphingolipid desaturase, which aid in understanding the evolution of these two enzymes, have not been reported. A blackcurrant Δ⁶-fatty acid desaturase and a Δ⁸-sphingolipid desaturase gene, RnD6C and RnD8A, respectively, share more than 80 % identity in their coding protein sequences. In this study, a set of fusion genes of RnD6C and RnD8A were constructed and expressed in yeast. The Δ⁶- and Δ⁸-desaturase activities of the fusion proteins were characterized. Our results indicated that (1) the exchange of the C-terminal 172 amino acid residues can lead to a significant decrease in both desaturase activities; (2) amino acid residues 114–174, 206–257, and 258–276 played important roles in Δ⁶-substrate recognition, and the last two regions were crucial for Δ⁸-substrate recognition; and (3) amino acid residues 114–276 of Δ⁶-fatty acid desaturase contained the substrate recognition site(s) responsible for discrimination between ceramide (a substrate of Δ⁸-sphingolipid desaturase) and acyl-PC (a substrate of Δ⁶-fatty acid desaturase). Substituting the amino acid residues 114-276 of RnD8A with those of RnD6C resulted in a gain of Δ⁶-desaturase activity in the fusion protein but a loss in Δ⁸-sphingolipid desaturase activity. In conclusion, several regions important for the substrate recognition of Δ8-sphingolipid desaturase and Δ⁶-fatty acid desaturase were identified, which provide clues in understanding the relationship between the structure and function in desaturases.     

Gain-of-function mutation in Gnao1: A murine model of epileptiform encephalopathy (EIEE17)?

G protein-coupled receptors strongly modulate neuronal excitability but there has been little evidence for G protein mechanisms in genetic epilepsies. Recently, four patients with epileptic encephalopathy (EIEE17) were found to have mutations in GNAO1, the most abundant G protein in brain, but the mechanism of this effect is not known. The GNAO1 gene product, Gα_o, negatively regulates neurotransmitter release. Here, we report a dominant murine model of Gnao1-related seizures and sudden death. We introduced a genomic gain-of-function knock-in mutation (Gnao1 ^+/G184S) that prevents G_o turnoff by Regulators of G protein signaling proteins. This results in rare seizures, strain-dependent death between 15 and 40 weeks of age, and a markedly increased frequency of interictal epileptiform discharges. Mutants on a C57BL/6J background also have faster sensitization to pentylenetetrazol (PTZ) kindling. Both premature lethality and PTZ kindling effects are suppressed in the 129SvJ mouse strain. We have mapped a 129S-derived modifier locus on Chromosome 17 (within the region 41–70 MB) as a Modifer of G protein Seizures (Mogs1). Our mouse model suggests a novel gain-of-function mechanism for the newly defined subset of epileptic encephalopathy (EIEE17). Furthermore, it reveals a new epilepsy susceptibility modifier Mogs1 with implications for the complex genetics of human epilepsy as well as sudden death in epilepsy.     

Chemical tools to explore nutrient-driven O-GlcNAc cycling

Abstract

Posttranslational modifications (PTM) including glycosylation, phosphorylation, acetylation, methylation and ubiquitination dynamically alter the proteome. The evolutionarily conserved enzymes O-linked N-acetylglucosamine (O-GlcNAc) transferase (OGT) and O-GlcNAcase are responsible for the addition and removal, respectively, of the nutrient-sensitive PTM of protein serine and threonine residues with O-GlcNAc. Indeed, the O-GlcNAc modification acts at every step in the “central dogma” of molecular biology and alters signaling pathways leading to amplified or blunted biological responses. The cellular roles of OGT and the dynamic PTM O-GlcNAc have been clarified with recently developed chemical tools including high-throughput assays, structural and mechanistic studies and potent enzyme inhibitors. These evolving chemical tools complement genetic and biochemical approaches for exposing the underlying biological information conferred by O-GlcNAc cycling.     

An S-glutathiomimetic Provides Structural Insights into Stromal Interaction Molecule-1 Regulation

Stromal interaction molecule 1 (STIM1) is an endo/sarcoplasmic reticulum (ER/SR) calcium (Ca²⁺) sensing protein that regulates store-operated calcium entry (SOCE). In SOCE, STIM1 activates Orai1-composed Ca²⁺ channels in the plasma membrane (PM) after ER stored Ca²⁺ depletion. S-Glutathionylation of STIM1 at Cys56 evokes constitutive SOCE in DT40 cells; however, the structural and biophysical mechanisms underlying the regulation of STIM1 by this modification are poorly defined. By establishing a protocol for site-specific STIM1 S-glutathionylation using reduced glutathione and diamide, we have revealed that modification of STIM1 at either Cys49 or Cys56 induces thermodynamic destabilization and conformational changes that result in increased solvent-exposed hydrophobicity. Further, S-glutathionylation or point-mutation of Cys56 reduces Ca²⁺ binding affinity, as measured by intrinsic fluorescence and far-UV circular dichroism spectroscopies. Solution NMR showed S-glutathionylated-induced perturbations in STIM1 are localized to the α1 helix of the canonical EF-hand, the α3 and α4 helices of the non-canonical EF-hand and α6 and α8 helices of the SAM domain. Finally, we designed an S-glutathiomimetic mutation that strongly recapitulates the structural, biophysical and functional effects within the STIM1 luminal domain and we envision to be another tool for understanding the effects of protein S-glutathionylation in vitro, in cellulo and in vivo.     

Protein engineering for covalent immobilization and enhanced stability through incorporation of multiple noncanonical amino acids

In this study, we demonstrate the application of multiple functional properties of proteins generated through coupling of residue-specific and site-specific incorporation method. With green fluorescent protein (GFP) as a model protein, we constructed multifunctional GFP through sitespecific incorporation of L-3,4-dihydroxyphenylalanine (DOPA) and residue-specific incorporation of (2S, 4S)-4- fluoroproline (4S-FP) or L-homopropargylglycine (hpg). Fluorescence analysis revealed a conjugation efficiency of approximately 20% for conjugation of DOPA-containing variants GFPdopa, GFPdp[4S-FP], and GFPdphpg onto chitosan. While incorporation of 4S-FP improved protein folding and stability, hpg incorporation into GFP allowed conjugation with fluorescent dye/polyethylene glycol (PEG). In addition, the modification of GFPhpg and GFPdphpg with PEG through Cu(I)-catalyzed click reaction increased protein thermal stability by about two-fold of the wild-type GFP.     

Inhibitors of DHHC family proteins

Protein S-acylation is a prevalent post-translational protein lipidation that is dynamically regulated by ‘writer’ protein S-acyltransferases and ‘eraser’ acylprotein thioesterases. The protein S-acyltransferases comprise 23 aspartate–histidine–histidine–cysteine (DHHC)–containing proteins, which transfer fatty acid acyl groups from acyl-coenzyme A onto protein substrates. DHHC proteins are increasingly recognized as critical regulators of S-acylation–mediated cellular processes and pathology. As our understanding of the importance and breadth of DHHC-mediated biology and pathology expands, so too does the need for chemical inhibitors of this class of proteins. In this review, we discuss the challenges and progress in DHHC inhibitor development, focusing on 2-bromopalmitate, the most commonly used inhibitor in the field, and N-cyanomethyl-N-myracrylamide, a new broad-spectrum DHHC inhibitor. We believe that current and ongoing advances in structure elucidation, mechanistic interrogation, and novel inhibitor design around DHHC proteins will spark innovative strategies to modulate these critical proteins in living systems.     

Development and Validation of a Method for Profiling Post-Translational Modification Activities Using Protein Microarrays

Background

Post-translational modifications (PTMs) impact on the stability, cellular location, and function of a protein thereby achieving a greater functional diversity of the proteome. To fully appreciate how PTMs modulate signaling networks, proteome-wide studies are necessary. However, the evaluation of PTMs on a proteome-wide scale has proven to be technically difficult. To facilitate these analyses we have developed a protein microarray-based assay that is capable of profiling PTM activities in complex biological mixtures such as whole-cell extracts and pathological specimens.

Methodology/Principal Findings

In our assay, protein microarrays serve as a substrate platform for in vitro enzymatic reactions in which a recombinant ligase, or extracts prepared from whole cells or a pathological specimen is overlaid. The reactions include labeled modifiers (e.g., ubiquitin, SUMO1, or NEDD8), ATP regenerating system, and other required components (depending on the assay) that support the conjugation of the modifier. In this report, we apply this methodology to profile three molecularly complex PTMs (ubiquitylation, SUMOylation, and NEDDylation) using purified ligase enzymes and extracts prepared from cultured cell lines and pathological specimens. We further validate this approach by confirming the in vivo modification of several novel PTM substrates identified by our assay.

Conclusions/Significance

This methodology offers several advantages over currently used PTM detection methods including ease of use, rapidity, scale, and sample source diversity. Furthermore, by allowing for the intrinsic enzymatic activities of cell populations or pathological states to be directly compared, this methodology could have widespread applications for the study of PTMs in human diseases and has the potential to be directly applied to most, if not all, basic PTM research.     

Thioredoxin 1 regulation of protein S-desulfhydration

The importance of H₂S in biology and medicine has been widely recognized in recent years, and protein S-sulfhydration is proposed to mediate the direct actions of H₂S bioactivity in the body. Thioredoxin 1 (Trx1) is an important reducing enzyme that cleaves disulfides in proteins and acts as an S-denitrosylase. The regulation of Trx1 on protein S-sulfhydration is unclear. Here we showed that Trx1 facilitates protein S-desulfhydration. Overexpression of Trx1 attenuated the basal level and H₂S-induced protein S-sulfhydration by direct interaction with S-sulfhydrated proteins, i.e., glyceraldehyde 3-phosphate dehydrogenase and pyruvate carboxylase. In contrast, knockdown of Trx1 mRNA expression by short interfering RNA or blockage of Trx1 redox activity with PX12 or 2,4-dinitrochlorobenzene enhanced protein S-sulfhydration. Mutation of cysteine-32 but not cysteine-35 in the Trp–Cys³²–Gly–Pro–Cys³⁵ motif eliminated the binding of Trx1 with S-sulfhydrated proteins and abolished the S-desulfhydrating effect of Trx1. All these data suggest that Trx1 acts as an S-desulfhydrase.     

GSHSite: Exploiting an Iteratively Statistical Method to Identify S-Glutathionylation Sites with Substrate Specificity

S-glutathionylation, the covalent attachment of a glutathione (GSH) to the sulfur atom of cysteine, is a selective and reversible protein post-translational modification (PTM) that regulates protein activity, localization, and stability. Despite its implication in the regulation of protein functions and cell signaling, the substrate specificity of cysteine S-glutathionylation remains unknown. Based on a total of 1783 experimentally identified S-glutathionylation sites from mouse macrophages, this work presents an informatics investigation on S-glutathionylation sites including structural factors such as the flanking amino acids composition and the accessible surface area (ASA). TwoSampleLogo presents that positively charged amino acids flanking the S-glutathionylated cysteine may influence the formation of S-glutathionylation in closed three-dimensional environment. A statistical method is further applied to iteratively detect the conserved substrate motifs with statistical significance. Support vector machine (SVM) is then applied to generate predictive model considering the substrate motifs. According to five-fold cross-validation, the SVMs trained with substrate motifs could achieve an enhanced sensitivity, specificity, and accuracy, and provides a promising performance in an independent test set. The effectiveness of the proposed method is demonstrated by the correct identification of previously reported S-glutathionylation sites of mouse thioredoxin (TXN) and human protein tyrosine phosphatase 1b (PTP1B). Finally, the constructed models are adopted to implement an effective web-based tool, named GSHSite (http://csb.cse.yzu.edu.tw/GSHSite/), for identifying uncharacterized GSH substrate sites on the protein sequences.     

Arabidopsis PROTEIN S‐ACYL TRANSFERASE4 mediates root hair growth

Polar growth of root hairs is critical for plant survival and requires fine‐tuned Rho of plants (ROP) signaling. Multiple ROP regulators participate in root hair growth. However, protein S‐acyl transferases (PATs), mediating the S‐acylation and membrane partitioning of ROPs, are yet to be found. Using a reverse genetic approach, combining fluorescence probes, pharmacological drugs, site‐directed mutagenesis and genetic analysis with related root‐hair mutants, we have identified and characterized an Arabidopsis PAT, which may be responsible for ROP2 S‐acylation in root hairs. Specifically, functional loss of PAT4 resulted in reduced root hair elongation, which was rescued by a wild‐type but not an enzyme‐inactive PAT4. Membrane‐associated ROP2 was significantly reduced in pat4, similar to S‐acylation‐deficient ROP2 in the wild type. We further showed that PAT4 and SCN1, a ROP regulator, additively mediate the stability and targeting of ROP2. The results presented here indicate that PAT4‐mediated S‐acylation mediates the membrane association of ROP2 at the root hair apex and provide novel insights into dynamic ROP signaling during plant tip growth.     

Post-translational modifications of the polycystin proteins

Autosomal dominant polycystic kidney disease (ADPKD) is the most common inherited cause of kidney failure and affects up to 12 million people worldwide. Germline mutations in two genes, PKD1 or PKD2, account for almost all patients with ADPKD. The ADPKD proteins, polycystin-1 (PC1) and polycystin-2 (PC2), are regulated by post-translational modifications (PTM), with phosphorylation, glycosylation and proteolytic cleavage being the best described changes. A few PTMs have been shown to regulate polycystin trafficking, signalling, localisation or stability and thus their physiological function. A key challenge for the future will be to elucidate the functional significance of all the individual PTMs reported to date. Finally, it is possible that site-specific mutations that disrupt PTM could contribute to cystogenesis although in the majority of cases, confirmatory evidence is awaited.