The Functional Human C-Terminome

All translated proteins end with a carboxylic acid commonly called the C-terminus. Many short functional sequences (minimotifs) are located on or immediately proximal to the C-terminus. However, information about the function of protein C-termini has not been consolidated into a single source. Here, we built a new “C-terminome” database and web system focused on human proteins. Approximately 3,600 C-termini in the human proteome have a minimotif with an established molecular function. To help evaluate the function of the remaining C-termini in the human proteome, we inferred minimotifs identified by experimentation in rodent cells, predicted minimotifs based upon consensus sequence matches, and predicted novel highly repetitive sequences in C-termini. Predictions can be ranked by enrichment scores or Gene Evolutionary Rate Profiling (GERP) scores, a measurement of evolutionary constraint. By searching for new anchored sequences on the last 10 amino acids of proteins in the human proteome with lengths between 3–10 residues and up to 5 degenerate positions in the consensus sequences, we have identified new consensus sequences that predict instances in the majority of human genes. All of this information is consolidated into a database that can be accessed through a C-terminome web system with search and browse functions for minimotifs and human proteins. A known consensus sequence-based predicted function is assigned to nearly half the proteins in the human proteome. Weblink: http://cterminome.bio-toolkit.com.     

Role of the carboxy-termini of tubulin on its chaperone-like activity.

Mutational analysis and the enzymatic digestion of many chaperones indicate the importance of both hydrophobic and hydrophilic residues for their unique property. Thus, the chaperone activity of alpha-crystallin is lost due to the substitution of hydrophobic residues or upon enzymatic digestion of the negatively charged residues. Tubulin, an eukaryotic cytoskeletal protein, exhibits chaperone-like activity as demonstrated by prevention of DTT-induced aggregation of insulin, thermal aggregation of alcohol dehydrogenase, betagamma-crystallin, and other proteins. We have shown that the tubulin lost its chaperone-like activity upon digestion of its negatively charged C-termini. In this article, the role of the C-terminus of individual subunits has been investigated. We observe that the digestion of C-terminus of beta-subunit with subtilisin causes loss of chaperone-like activity of tubulin. The contribution of C-terminus of alpha-subunit is difficult to establish directly as subtilisin cleaves C-terminus of beta-subunit first. This has been ascertained indirectly using a 14-residue peptide P2 having the sequence corresponding to a conserved region of MHC class I molecules and that binds tightly to the C-terminus of alpha-subunit. We have shown that the binding of P2 peptide to alphabeta-tubulin causes complete loss of its chaperone-like activity. NMR and gel-electrophoresis studies indicate that the P2 peptide has a significant higher binding affinity for the C-terminus of alpha-subunit compared to that of beta-subunit. Thus, we conclude that both the C-termini are necessary for the chaperone-like activity of tubulin. Implications for the chaperone functions in vivo have been discussed.     

Binding of coatomer by the PEX11 C-terminus is not required for function

Microbodies are single membrane-bound organelles found in eukaryotes from trypanosomes to man. Although they have diverse roles in metabolism, the mechanisms and molecules involved in membrane biogenesis and matrix protein import are conserved. Similarly, the basic mechanisms and structures involved in vesicular transport are similar throughout eukaryotic evolution. The PEX11 proteins are required for the division of microbodies in trypanosomes, yeast and mammals, and a role of coatomer in this process has been suggested. We show here that the binding of trypanosome, yeast and bovine coatomers to selected peptides is identical. Coatomer binds to the C-termini of trypanosome PEX11 and rat Pex11alpha, but not yeast Pex11p or human Pex11beta. Mutations of the C-terminus of trypanosome PEX11 that eliminated coatomer binding did not affect function in yeast or trypanosomes. Thus binding of coatomer to the C-terminus of PEX11 is not required for PEX11 function.     

Eukaryotic genomes contain a [2Fez.sbnd;2S] ferredoxin isoform with a conserved C-terminal sequence motif

Apicomplexan protists contain a single mitochondrial [2Fe-2S] ferredoxin sequence (mtFd) with a highly conserved C-terminal motif, VDGxxpxPH, that distinguishes it from other mtFd, which have heterogeneous C-termini. This isoform of mtFd, called 'type II ferredoxin', is widespread in eukaryotes, some species having two isoforms and others possessing only one. Because of the known modulating role of the C-terminus of type I mtFd during association with itself and other interacting proteins, the presence of a conserved C-terminus in type II mtFd suggests it evolved either as a means for optimized homodimerization or to allow interaction with a highly conserved partner(s) that is yet to be defined.     

Secondary Structure,a Missing Component of Sequence-Based Minimotif Definitions

Minimotifs are short contiguous segments of proteins that have a known biological function. The hundreds of thousands of minimotifs discovered thus far are an important part of the theoretical understanding of the specificity of protein-protein interactions, posttranslational modifications, and signal transduction that occur in cells. However, a longstanding problem is that the different abstractions of the sequence definitions do not accurately capture the specificity, despite decades of effort by many labs. We present evidence that structure is an essential component of minimotif specificity, yet is not used in minimotif definitions. Our analysis of several known minimotifs as case studies, analysis of occurrences of minimotifs in structured and disordered regions of proteins, and review of the literature support a new model for minimotif definitions that includes sequence, structure, and function.     

An isotope coding strategy for proteomics involving both amine and carboxyl group labeling

This paper describes a heavy isotope coding strategy for the analysis of all types of tryptic peptides, including those that are N-terminally blocked and from the C-terminus of proteins. The method exploits differential derivatization of amine and carboxyl groups generated during proteolysis as a means of coding. Carboxyl groups produced during proteolysis incorporate 18O from H218O. Peptides from the C-terminus of proteins were not labeled with 18O unless they contained a basic C-terminal amino acid. Primary amines from control and experimental samples were differentially acylated after proteolysis with either 1H3- or 2H3-N-acetoxysuccinamide. When these two types of labeling were combined, unique coding patterns were achieved for peptides arising from the C-termini and blocked N-termini of proteins. This method was used to (1) distinguish C-terminal peptides in model proteins, (2) recognize N-terminal peptides from proteins in which the amino terminus is acylated, and (3) identify primary structure variations between proteins from different sources.     

Unique gating properties of C. elegans ClC anion channel splice variants are determined by altered CBS domain conformation and the R-helix linker

All eukaryotic and some prokaryotic ClC anion transport proteins have extensive cytoplasmic C-termini containing two cystathionine-β-synthase (CBS) domains. CBS domain secondary structure is highly conserved and consists of two a-helices and three b-strands arranged as b1-a1-b2-b3-a2. ClC CBS domain mutations cause muscle and bone disease and alter ClC gating. However, the precise functional roles of CBS domains and the structural bases by which they regulate ClC function are poorly understood. CLH-3a and CLH-3b are C. elegans ClC anion channel splice variants with strikingly different biophysical properties. Splice variation occurs at cytoplasmic N- and C-termini and includes several amino acids that form a2 of the second CBS domain (CBS2). We demonstrate that interchanging a2 between CLH-3a and CLH-3b interchanges their gating properties. The "R-helix" of ClC proteins forms part of the ion conducting pore and selectivity filter and is connected to the cytoplasmic C-terminus via a short stretch of cytoplasmic amino acids termed the "R-helix linker". C-terminus conformation changes could cause R-helix structural rearrangements via this linker. X-ray structures of three ClC protein cytoplasmic C-termini suggest that a2 of CBS2 and the R-helix linker could be closely apposed and may therefore interact. We found that mutating apposing amino acids in a2 and the R-helix linker of CLH-3b was sufficient to give rise to CLH-3a-like gating. We postulate that the R-helix linker interacts with CBS2 a2, and that this putative interaction provides a pathway by which cytoplasmic C-terminus conformational changes induce conformational changes in membrane domains that in turn modulate ClC function.     

A Nonsynonymous/Synonymous Substitution Analysis of the B56 Gene Family Aids in Understanding B56 Isoform Diversity

Gene duplication leads to the formation of gene families, wherein purifying or neutral selection maintains the original gene function, while diversifying selection confers new functions onto duplicated genes. The B56 gene family is highly conserved; it is encoded by one gene in protists and fungi, and five genes in vertebrates. B56 regulates protein phosphatase 2A (PP2A), an abundant heterotrimeric serine/threonine phosphatase that functions as a tumor suppressor and consists of a scaffolding “A” and catalytic “C” subunit heterodimer bound to a regulatory “B” subunit. Individual regulatory B56 subunits confer disparate functions onto PP2A in various cell-cell signaling pathways. B56 proteins share a conserved central core domain, but have divergent N- and C-termini which play a role in isoform specificity. We carried out a nonsynonymous/synonymous substitution analysis to better understand the divergence of vertebrate B56 genes. When five B56 paralogs from ten vertebrate species were analyzed, the gene family displayed purifying selection; stronger purifying selection was revealed when individual B56 isoforms were analyzed separately. The B56 core experienced stronger purifying selection than the N- and C-termini, which correlates with the presence of several contacts between the core and the AC heterodimer. Indeed, the majority of the contact points that we analyzed between B56 and the AC heterodimer experienced strong purifying selection. B56 subfamilies showed distinct patterns of selection in their N- and C-termini. The C-terminus of the B56-1 subfamily and the N-terminus of the B56-2 subfamily exhibited strong purifying selection, suggesting that these termini carry out subfamily-specific functions, while the opposite termini exhibited diversifying selection and likely carry out isoform-specific functions. We also found reduced synonymous substitutions at the N- and C-termini when grouping B56 genes by species but not by isoform, suggesting species-specific codon bias may have a role in regulating B56 gene expression.     

Synthesis and screening of support-bound combinatorial peptide libraries with free C-termini: determination of the sequence specificity of PDZ domains

Preparation of support-bound combinatorial peptide libraries with free C-termini has been challenging in the past because solid-phase peptide synthesis usually starts from the C-terminus, which must be covalently attached to the solid support. In this work, we have developed a general methodology to synthesize and screen one-bead-one-compound peptide libraries containing free C-termini. TentaGel microbeads (90 mum) were spatially segregated into outer and inner layers, and peptides were synthesized on the beads in the conventional C --> N manner, with their C-termini attached to the support through an ester linkage on the bead surface but through an amide bond in the bead interior. The surface peptides were cyclized between their N-terminal amine and a carboxyl group installed at a C-terminal linker sequence, while the internal peptides were kept in the linear form. Base hydrolysis of the ester linkage in the cyclic peptides regenerated linear peptides that contained a free alpha-carboxyl group at their C-termini but remained covalently attached to the resin via the N-termini ("inverted" peptides). An inverted peptide library containing five random residues (theoretical diversity of 3.2 x 10 (6)) was synthesized and screened for binding to four postsynaptic density-95/discs large/zona occluden-1 (PDZ) domains of sodium-hydrogen exchanger regulatory factor-1 (NHERF1) and channel-interacting PDZ domain protein (CIPP). The identity of the binding peptides was determined by sequencing the linear encoding peptides inside the bead by partial Edman degradation/mass spectrometry. Consensus recognition motifs were identified for the PDZ domains, and representative peptides were resynthesized and confirmed for binding to their cognate PDZ domains. This method should be generally applicable to all PDZ domains as well as other protein domains and enzymes that recognize the C-terminus of their target proteins.     

Distinct overlapping sequences at the carboxy-terminus of merlin regulate its tumour suppressor and morphogenic activity

The Neurofibromatosis 2 (NF2) gene product merlin is a tumour suppressor, which in addition to inhibiting cell proliferation regulates cell morphology. The morphogenic properties of merlin may play a role in tumour suppression, as patient-derived tumour cells demonstrate cytoskeletal abnormalities. However, it is still unclear how these functions are linked. The N-terminal FERM-domain of merlin is highly homologous to the oncogenic protein ezrin, while the C-termini are less conserved, suggesting that the opposite effect of the proteins on proliferation could be mediated by their distinct C-terminal regions. In this study we characterize the role of the most C-terminal residues of merlin in the regulation of proliferation, cytoskeletal organization, phosphorylation and intramolecular associations. In addition to the two full-length merlin isoforms and truncating mutations found in patients, we focused on the evolutionally conserved C-terminal residues 545-547, also harbouring disease-causing mutations. We demonstrate that merlin induces cell extensions, which result from impaired retraction of protrusions rather than from increased formation of filopodia. The residues 538-568 were found particularly important for this morphogenic activity. The results further show that both merlin isoforms are able to equally inhibit proliferation, whereas C-terminal mutants affecting residues 545-547 are less effective in growth suppression. This study demonstrates that the C-terminus contains distinct but overlapping functional domains important for regulation of the morphogenic activity, intramolecular associations and cell proliferation.     

Structural insights into the dynamics and function of the C-terminus of the E. coli RNA chaperone Hfq

The hexameric Escherichia coli RNA chaperone Hfq (Hfq(Ec)) is involved in riboregulation of target mRNAs by small trans-encoded RNAs. Hfq proteins of different bacteria comprise an evolutionarily conserved core, whereas the C-terminus is variable in length. Although the structure of the conserved core has been elucidated for several Hfq proteins, no structural information has yet been obtained for the C-terminus. Using bioinformatics, nuclear magnetic resonance spectroscopy, synchrotron radiation circular dichroism (SRCD) spectroscopy and small angle X-ray scattering we provide for the first time insights into the conformation and dynamic properties of the C-terminal extension of Hfq(Ec). These studies indicate that the C-termini are flexible and extend laterally away from the hexameric core, displaying in this way features typical of intrinsically disordered proteins that facilitate intermolecular interactions. We identified a minimal, intrinsically disordered region of the C-terminus supporting the interactions with longer RNA fragments. This minimal region together with rest of the C-terminal extension provides a flexible moiety capable of tethering long and structurally diverse RNA molecules. Furthermore, SRCD spectroscopy supported the hypothesis that RNA fragments exceeding a certain length interact with the C-termini of Hfq(Ec).     

Partitioning of Minimotifs Based on Function with Improved Prediction Accuracy

Background

Minimotifs are short contiguous peptide sequences in proteins that are known to have a function in at least one other protein. One of the principal limitations in minimotif prediction is that false positives limit the usefulness of this approach. As a step toward resolving this problem we have built, implemented, and tested a new data-driven algorithm that reduces false-positive predictions.

Methodology/Principal Findings

Certain domains and minimotifs are known to be strongly associated with a known cellular process or molecular function. Therefore, we hypothesized that by restricting minimotif predictions to those where the minimotif containing protein and target protein have a related cellular or molecular function, the prediction is more likely to be accurate. This filter was implemented in Minimotif Miner using function annotations from the Gene Ontology. We have also combined two filters that are based on entirely different principles and this combined filter has a better predictability than the individual components.

Conclusions/Significance

Testing these functional filters on known and random minimotifs has revealed that they are capable of separating true motifs from false positives. In particular, for the cellular function filter, the percentage of known minimotifs that are not removed by the filter is ∼4.6 times that of random minimotifs. For the molecular function filter this ratio is ∼2.9. These results, together with the comparison with the published frequency score filter, strongly suggest that the new filters differentiate true motifs from random background with good confidence. A combination of the function filters and the frequency score filter performs better than these two individual filters.     

Conserved C-termini of Spidroins are secreted by the major ampullate glands and retained in the silk thread

The C-termini of Spidroins produced in the major and minor ampullate glands of spiders are highly conserved. Despite this conservation, no corresponding peptides have been identified in the spinning dopes or the silk filaments so far. To prove their presence or absence, polyclonal antibodies derived against fusion proteins containing the conserved C-terminal regions of both Spidroin 1 and 2 from the spider Nephila clavipes were generated. The antibodies reacted with high molecular weight polypeptides of the corresponding gland extracts and solubilized major ampullate filament and in addition to filament cross-sections. This demonstrates the existence of C-terminal specific peptides in the spinning dope and the mature Spidroins. Both the fusion proteins as well as the proteins contained within the gland lumen showed a reduction in their size under reducing conditions indicating the presence of disulfide bonds. Their high conservation and the biochemical data suggest crucial roles the C-termini play in the formation and/or structure of the corresponding silk filaments.     

Characterization of the mouse pancreatic islet proteome and comparative analysis with other mouse tissues

The pancreatic islets of Langerhans, and especially the insulin-producing beta cells, play a central role in the maintenance of glucose homeostasis. Alterations in the expression of multiple proteins in the islets that contribute to the maintenance of islet function are likely to underlie the pathogenesis of types 1 and 2 diabetes. To identify proteins that constitute the islet proteome, we provide the first comprehensive proteomic characterization of pancreatic islets for mouse, the most commonly used animal model in diabetes research. Using strong cation exchange fractionation coupled with reversed phase LC-MS/MS we report the confident identification of 17,350 different tryptic peptides covering 2612 proteins having at least two unique peptides per protein. The data set also identified approximately 60 post-translationally modified peptides including oxidative modifications and phosphorylation. While many of the identified phosphorylation sites corroborate those previously known, the oxidative modifications observed on cysteinyl residues reveal potentially novel information suggesting a role for oxidative stress in islet function. Comparative analysis with 15 available proteomic data sets from other mouse tissues and cells revealed a set of 133 proteins predominantly expressed in pancreatic islets. This unique set of proteins, in addition to those with known functions such as peptide hormones secreted from the islets, contains several proteins with as yet unknown functions. The mouse islet protein and peptide database accessible at (http://ncrr.pnl.gov), provides an important reference resource for the research community to facilitate research in the diabetes and metabolism fields.     

A computational tool for identifying minimotifs in protein-protein interactions and improving the accuracy of minimotif predictions

Protein-protein interactions are important to understanding cell functions; however, our theoretical understanding is limited. There is a general discontinuity between the well-accepted physical and chemical forces that drive protein-protein interactions and the large collections of identified protein-protein interactions in various databases. Minimotifs are short functional peptide sequences that provide a basis to bridge this gap in knowledge. However, there is no systematic way to study minimotifs in the context of protein-protein interactions or vice versa. Here we have engineered a set of algorithms that can be used to identify minimotifs in known protein-protein interactions and implemented this for use by scientists in Minimotif Miner. By globally testing these algorithms on verified data and on 100 individual proteins as test cases, we demonstrate the utility of these new computation tools. This tool also can be used to reduce false-positive predictions in the discovery of novel minimotifs. The statistical significance of these algorithms is demonstrated by an ROC analysis (P = 0.001).     

Conformational changes in bacteriophage P22 scaffolding protein induced by interaction with coat protein

Many prokaryotic and eukaryotic double-stranded DNA viruses use a scaffolding protein to assemble their capsid. Assembly of the double-stranded DNA bacteriophage P22 procapsids requires the interaction of 415 molecules of coat protein and 60-300 molecules of scaffolding protein. Although the 303-amino-acid scaffolding protein is essential for proper assembly of procapsids, little is known about its structure beyond an NMR structure of the extreme C-terminus, which is known to interact with coat protein. Deletion mutagenesis indicates that other regions of scaffolding protein are involved in interactions with coat protein and other capsid proteins. Single-cysteine and double-cysteine variants of scaffolding protein were generated for use in fluorescence resonance energy transfer and cross-linking experiments designed to probe the conformation of scaffolding protein in solution and within procapsids. We showed that the N-terminus and the C-terminus are proximate in solution, and that the middle of the protein is near the N-terminus but not accessible to the C-terminus. In procapsids, the N-terminus was no longer accessible to the C-terminus, indicating that there is a conformational change in scaffolding protein upon assembly. In addition, our data are consistent with a model where scaffolding protein dimers are positioned parallel with one another with the associated C-termini.     

Presentation of antigenic peptides by products of the major histocompatibility complex

Molecules encoded by the major histocompatibility complex (MHC) are polymorphic integral membrane proteins adapted to the presentation of peptide fragments of foreign antigens to antigen-specific T-cells. The diversity of infectious agents to which an immune response must be mounted poses a unique problem for receptor–ligand interactions; how can proteins whose polymorphism is necessarily limited bind an array of peptides almost infinite in its complexity? Both MHC class I and class II determinants have achieved this goal by harnessing a limited number of peptide side chains to anchor the epitope in place while exploiting conserved features of peptide structure, independent of their primary sequence. While class I molecules interact predominantly with the N- and C-termini of peptides, class II determinants form an extensive hydrogen bonding network along the length of the peptide backbone. Such a strategy ensures high-affinity binding, while selectively exposing the unique features of each ligand for recognition by the T-cell receptor. © 1998 European Peptide Society and John Wiley & Sons, Ltd.     

Snakin/GASA proteins: Involvement in hormone crosstalk and redox homeostasis

Snakin/GASA proteins are widely distributed among plant species. They are expressed in different plant organs with high tissue and temporal specificity, and their subcellular localization varies among the different members. Interestingly, all of them maintain 12 cysteines of the C-terminus in highly conserved positions of the aminoacid sequences that are essential for their biochemical activity and probably responsible for their protein structure. Despite their common features, their functions are not completely elucidated and little is known about their mode of action. This review focuses on the current knowledge about this intriguing family of peptides and advances comprising gene regulation analyses, expression pattern studies and phenotypic characterization of mutants and transgenic plants. Furthermore, we discuss the roles of Snakin/GASA proteins in several aspects of plant development, plant responses to biotic or abiotic stress and their participation in hormone crosstalk and redox homeostasis.     

Heat shock proteins: essential proteins for apoptosis regulation

Many different external and intrinsic apoptotic stimuli induce the accumulation in the cells of a set of proteins known as stress or heat shock proteins (HSPs). HSPs are conserved proteins present in both prokaryotes and eukaryotes. These proteins play an essential role as molecular chaperones by assisting the correct folding of nascent and stress-accumulated misfolded proteins, and by preventing their aggregation. HSPs have a protective function, that is they allow the cells to survive to otherwise lethal conditions. Various mechanisms have been proposed to account for the cytoprotective functions of HSPs. Several of these proteins have demonstrated to directly interact with components of the cell signalling pathways, for example those of the tightly regulated caspase-dependent programmed cell death machinery, upstream, downstream and at the mitochondrial level. HSPs can also affect caspase-independent apoptosis-like process by interacting with apoptogenic factors such as apoptosis-inducing factor (AIF) or by acting at the lysosome level. This review will describe the different key apoptotic proteins interacting with HSPs and the consequences of these interactions in cell survival, proliferation and apoptotic processes. Our purpose will be illustrated by emerging strategies in targeting these protective proteins to treat haematological malignancies.