Motif refinement of the peroxisomal targeting signal 1 and evaluation of taxon-specific differences

A 27.368 kDa retinal reductase in New Zealand white rabbit liver cytosol encoded by the peroxisomal retinol dehydrogenase-reductase cDNA: purification and characterization of the enzyme

We obtained a full-length cDNA based on a sequence deposited in GenBank (accession No. AB045133), annotated as rabbit peroxisomal NADP(H)-dependent retinol dehydrogenase-reductase (NDRD). The rabbit NDRD gene, like its mouse and human homologs, harbors 2 initiation sites, one of which theoretically encodes a 29.6 kDa protein with 279 amino acids, and the other encodes a 27.4 kDa protein with 260 amino acids. The purification of a rabbit cytosolic retinol oxidoreductase with a subunit molecular mass of 34 kDa and an N terminus that is not completely identical to that of NDRD, has been reported. An enzyme responsible for the all-trans retinal reductase activity in the liver cytosol of New Zealand white rabbit was purified to homogeneity using differential centrifugation and successive chromatographic analyses. The subunit molecular mass of the purified enzyme, revealed by SDS-PAGE, was approximately 27 kDa. The intact molecular mass, measured by MALDI-TOF mass spectrometry, was 27.368 kDa. The 60 kDa relative mobility observed in size-exclusion chromatography indicates that the native protein probably exists as a dimer. The purified enzyme was positively confirmed to be the product of NDRD by peptide mass fingerprinting, tandem mass spectrometry, and N-terminal sequencing. Taken together, the results suggested that the native protein is truncated at the N terminus.

Quantitative analysis of peroxisomal targeting signal type-1 binding to wild-type and pathogenic mutants of Pex5p supports an affinity threshold for peroxisomal protein targeting

Peroxisomal biogenesis disorders (PBDs) are caused by mutations in 12 distinct genes that encode the components of the peroxisome assembly machinery. Three mutations in the gene encoding Pex5p, the peroxisomal targeting signal type-1 (PTS1) receptor, have been reported, each associated with a disorder of the Zellweger spectrum of different severity. Here, we report studies of the affinities of mutated forms of Pex5p for a series of PTS1 peptides and conclude that PTS1-affinity reductions are correlated with disease severity and cell biological phenotype. A quantitative model has been developed that allows estimation of the dissociation constants for complexes with a wide range of PTS1 sequences bound to wild-type and mutant Pex5p. In the context of this model, the binding measurements suggest that no PTS1-containing proteins are targeted by Pex5p(N489K) and only a relatively small subset of PTS1-containing proteins with the highest affinity for Pex5p are targeted to peroxisomes by Pex5p(S563W). Furthermore, the results of the analysis are consistent with an approximate dissociation constant threshold near 500 nM required for efficient protein targeting to peroxisomes.

pkaPS: prediction of protein kinase A phosphorylation sites with the simplified kinase-substrate binding model

Background

Protein kinase A (cAMP-dependent kinase, PKA) is a serine/threonine kinase, for which ca. 150 substrate proteins are known. Based on a refinement of the recognition motif using the available experimental data, we wished to apply the simplified substrate protein binding model for accurate prediction of PKA phosphorylation sites, an approach that was previously successful for the prediction of lipid posttranslational modifications and of the PTS1 peroxisomal translocation signal.

Results

Approximately 20 sequence positions flanking the phosphorylated residue on both sides have been found to be restricted in their sequence variability (region -18...+23 with the site at position 0). The conserved physical pattern can be rationalized in terms of a qualitative binding model with the catalytic cleft of the protein kinase A. Positions -6...+4 surrounding the phosphorylation site are influenced by direct interaction with the kinase in a varying degree. This sequence stretch is embedded in an intrinsically disordered region composed preferentially of hydrophilic residues with flexible backbone and small side chain. This knowledge has been incorporated into a simplified analytical model of productive binding of substrate proteins with PKA.

Conclusion

The scoring function of the pkaPS predictor can confidently discriminate PKA phosphorylation sites from serines/threonines with non-permissive sequence environments (sensitivity of ~96% at a specificity of ~94%). The tool "pkaPS" has been applied on the whole human proteome. Among new predicted PKA targets, there are entirely uncharacterized protein groups as well as apparently well-known families such as those of the ribosomal proteins L21e, L22 and L6.

Availability

The supplementary data as well as the prediction tool as WWW server are available at .

Reviewers

Erik van Nimwegen (Biozentrum, University of Basel, Switzerland), Sandor Pongor (International Centre for Genetic Engineering and Biotechnology, Trieste, Italy), Igor Zhulin (University of Tennessee, Oak Ridge National Laboratory, USA).

Recognition of a functional peroxisome type 1 target by the dynamic import receptor pex5p

Peroxisomes require the translocation of folded and functional target proteins of various sizes across the peroxisomal membrane. We have investigated the structure and function of the principal import receptor Pex5p, which recognizes targets bearing a C-terminal peroxisomal targeting signal type 1. Crystal structures of the receptor in the presence and absence of a peroxisomal target, sterol carrier protein 2, reveal major structural changes from an open, snail-like conformation into a closed, circular conformation. These changes are caused by a long loop C terminal to the 7-fold tetratricopeptide repeat segments. Mutations in residues of this loop lead to defects in peroxisomal import in human fibroblasts. The structure of the receptor/cargo complex demonstrates that the primary receptor-binding site of the cargo is structurally and topologically autonomous, enabling the cargo to retain its native structure and function.

Computational prediction of subcellular localization

It is widely recognized that much of the information for determining the final subcellular localization of proteins is found in their amino acid sequences. Thus the prediction of protein localization sites is of both theoretical and practical interest. In most cases, the prediction has been attempted in two ways: one is based on the knowledge of experimentally characterized targeting signals, while the other utilizes the statistical differences of general sequence characteristics, such as amino acid composition, between localization sites. Both approaches have limitations, and it is recommended to check the results of various prediction methods based on different principles as well as training data. Recently, increased proteomic analyses of localization sites have provided new data to assess the current status of predictive methods. In this chapter we discuss these issues and close with an example illustrating the use of the WoLF PSORT web server for localization prediction.

Localization of the pre-squalene segment of the isoprenoid biosynthetic pathway in mammalian peroxisomes

Previous studies have indicated that the early steps in the isoprenoid/cholesterol biosynthetic pathway occur in peroxisomes. However, the role of peroxisomes in cholesterol biosynthesis has recently been questioned in several reports concluding that three of the peroxisomal cholesterol biosynthetic enzymes, namely mevalonate kinase, phosphomevalonate kinase, and mevalonate diphosphate decarboxylase, do not localize to peroxisomes in human cells even though they contain consensus peroxisomal targeting signals. We re-investigated the subcellular localization of the cholesterol biosynthetic enzymes of the pre-squalene segment in human cells by using new stable isotopic techniques and data computations with isotopomer spectral analysis, in combination with immunofluorescence and cell permeabilization techniques. Our present findings clearly show and confirm previous studies that the pre-squalene segment of the cholesterol biosynthetic pathway is localized to peroxisomes. In addition, our data are consistent with the hypothesis that acetyl-CoA derived from peroxisomal beta-oxidation of very long-chain fatty acids and medium-chain dicarboxylic acids is preferentially channeled to cholesterol synthesis inside the peroxisomes without mixing with the cytosolic acetyl-CoA pool.

Analysis of the targeting sequences of an iron-containing superoxide dismutase (SOD) of the dinoflagellate Lingulodinium polyedrum suggests function in multiple cellular compartments

One of the proteins targeted to the peridinin plastid of the dinoflagellate Lingulodinium polyedrum is the iron-containing superoxide dismutase (LpSOD). Like dinoflagellate plastid proteins of class II, LpSOD carries a bipartite presequence comprising a signal peptide followed by a transit peptide. Our bioinformatic studies suggest that its signal peptide is atypical, however, and that the entire presequence may function as a mitochondrial targeting signal. It is possible that LpSOD represents a new class of proteins in algae with complex plastids, which are co-targeted to the plastid and mitochondrion. In addition to the ambiguous N-terminal targeting signal, LpSOD contains a potential type-1 peroxisome-targeting signal (PTS1) located at its C-terminus. In accordance with a peroxisome localization of this dismutase, its mRNA has two in-frame AUG codons. Our bioinformatic analyses indicate that the first start codon resides in a much weaker oligonucleotide context than the second one. This suggests that synthesis of the plastid/mitochondrion-targeted and peroxisome-targeted isoforms could proceed through so-called leaky scanning. Moreover, our results show that expression of the two isoforms could be regulated by a 'hairpin' structure located between the first and second start codons.

PeroxisomeDB: a database for the peroxisomal proteome, functional genomics and disease

Peroxisomes are essential organelles of eukaryotic origin, ubiquitously distributed in cells and organisms, playing key roles in lipid and antioxidant metabolism. Loss or malfunction of peroxisomes causes more than 20 fatal inherited conditions. We have created a peroxisomal database () that includes the complete peroxisomal proteome of Homo sapiens and Saccharomyces cerevisiae, by gathering, updating and integrating the available genetic and functional information on peroxisomal genes. PeroxisomeDB is structured in interrelated sections ‘Genes’, ‘Functions’, ‘Metabolic pathways’ and ‘Diseases’, that include hyperlinks to selected features of NCBI, ENSEMBL and UCSC databases. We have designed graphical depictions of the main peroxisomal metabolic routes and have included updated flow charts for diagnosis. Precomputed BLAST, PSI-BLAST, multiple sequence alignment (MUSCLE) and phylogenetic trees are provided to assist in direct multispecies comparison to study evolutionary conserved functions and pathways. Highlights of the PeroxisomeDB include new tools developed for facilitating (i) identification of novel peroxisomal proteins, by means of identifying proteins carrying peroxisome targeting signal (PTS) motifs, (ii) detection of peroxisomes in silico, particularly useful for screening the deluge of newly sequenced genomes. PeroxisomeDB should contribute to the systematic characterization of the peroxisomal proteome and facilitate system biology approaches on the organelle.

Mechanisms of copper loading on the Schizosaccharomyces pombe copper amine oxidase 1 expressed in Saccharomyces cerevisiae

Copper amine oxidases (CAOs) are found in almost every living kingdom. Although Saccharomyces cerevisiae is one of the few yeast species that lacks an endogenous CAO, heterologous gene expression of CAOs from other organisms produces a functional enzyme. To begin to characterize their function and mechanisms of copper acquisition, two putative cao(+) genes from Schizosaccharomyces pombe were expressed in S. cerevisiae. Expression of spao1(+) resulted in the production of an active enzyme capable of catalysing the oxidative deamination of primary amines. On the other hand, expression of spao2(+) failed to produce an active CAO. Using a functional spao1(+)-GFP fusion allele, the SPAO1 protein was localized in the cytosol. Under copper-limiting conditions, yeast cells harbouring deletions of the MAC1, CTR1 and CTR3 genes were defective in amine oxidase activity. Likewise, atx1Delta null cells exhibited no CAO activity, while ccc2Delta mutant cells exhibited decreased levels of amine oxidase activity, and mutations in cox17Delta and ccs1Delta did not cause any defects in this activity. Copper-deprived S. cerevisiae cells expressing spao1(+) required a functional atx1(+) gene for growth on minimal medium containing ethylamine as the sole nitrogen source. Under these conditions, the inability of the atx1Delta cells to utilize ethylamine correlated with the lack of SPAO1 activity, in spite of the efficient expression of the protein. Cells carrying a disrupted ccc2Delta allele exhibited only weak growth on ethylamine medium containing a copper chelator. The results of these studies reveal that expression of the heterologous spao1(+) gene in S. cerevisiae is required for its growth in medium containing ethylamine as the sole nitrogen source, and that expression of an active Schiz. pombe SPAO1 protein in S. cerevisiae depends on the acquisition of copper through the high-affinity copper transporters Ctr1 and Ctr3, and the copper chaperone Atx1.

A yeast sterol carrier protein with fatty-acid and fatty-acyl-CoA binding activity

The 14-kDa sterol carrier protein 2 (SCP2) domain is present in Eukaria, Bacteria and Archaea, and has been implicated in the transport and metabolism of lipids. We report the cloning, expression, purification and physicochemical characterization of a SCP2 from the yeast Yarrowia lipolytica (YLSCP2). Analytical size-exclusion chromatography, circular dichroism and fluorescence spectra, indicate that recombinant YLSCP2 is a well-folded monomer. Thermal unfolding experiments show that SCP2 maximal stability is at pH 7.0-9.0. YLSCP2 binds cis-parinaric acid and palmitoyl-CoA with KD values of 81+/-40 nM and 73+/-33 nM, respectively, sustaining for the first time the binding of fatty acids and their CoA esters to a nonanimal SCP2. The role of yeast SCP2 and other lipid binding proteins in transport, storage and peroxisomal oxidation of fatty acids is discussed.

Peroxisome targeting signal 1: is it really a simple tripeptide?

Originally, the peroxisomal targeting signal 1 (PTS1) was defined as a tripeptide at the C-terminus of proteins prone to be imported into the peroxisomal matrix. The corresponding receptor PEX5 initiates the translocation of proteins by identifying potential substrates via their C-termini and trapping PTS1s through remodeling of its TPR domain. Thorough studies on the interaction between PEX5 and PTS1 as well as sequence-analytic tools revealed the influence of amino acid residues further upstream of the ultimate tripeptide. Altogether, PTS1s should be defined as dodecamer sequences at the C-terminal ends of proteins. These sequences accommodate physical contacts with both the surface and the binding cavity of PEX5 and ensure accessibility of the extreme C-terminus. Knowledge-based approaches in applied Bioinformatics provide reliable tools to accurately predict the peroxisomal location of proteins not yet determined experimentally.

In silico prediction of the glycosomal enzymes of Leishmania major and trypanosomes

In total, 37080 protein sequences of the three trypanosomatids Leishmania major, Trypanosoma brucei and Trypanosoma cruzi, were used to predict the trypanosomatid glycosomal proteome. All protein sequences were analyzed for the presence of either a C-terminal (PTS1) or an N-terminal (PTS2) peroxisomal targeting sequence. For L. major 191 potential PTS1-containing proteins and 68 potential PTS2-containing proteins with homologues in T. brucei and T. cruzi were identified. About 50% of them were hypothetical proteins to which no function was attributed. From those proteins with known function it appears that the predicted glycosomal proteome of L. major strongly resembles that of T. brucei and T. cruzi with respect to enzyme content. Glycosomes are not only involved in glycolysis, but are predicted to carry out also gluconeogenesis, reactions of the hexose-monophosphate pathway, purine salvage and pyrimidine biosynthesis, beta-oxidation of fatty acids, fatty acid elongation and the biosynthesis of ether lipids. In addition, they seem to catalyze several reactions of isoprenoid synthesis and are involved in oxidant stress protection.

Identification of paralogous genes of firefly luciferase in the Japanese firefly, Luciola cruciata

Two homologous genes of firefly luciferase, LcLL1 and LcLL2, were cloned from the Japanese firefly Luciola cruciata, and were expressed and characterized. The gene product of LcLL1 had long-chain fatty acyl-CoA synthetic activity, but not luciferase activity. The other gene product of LcLL2 did not show enzymatic activities of acyl-CoA synthetase and luciferase. RT-PCR analysis showed that the transcript of LcLL1 was abundant in larva but very low in adult, while LcLL2 was expressed in both larva and adult. Phylogenetic analysis indicated that LcLL1 and LcLL2 are paralogous genes of firefly luciferase. Recently, we found that CG6178 in Drosophila melanogaster is an orthologue of firefly luciferase and shows fatty acyl-CoA synthetic activity, but not luciferase activity. These results suggest that firefly luciferase might be evolved from a fatty acyl-CoA synthetase by gene duplication in insects.

Sequence-based discovery of the human and rodent peroxisomal proteome

BACKGROUND

Peroxisomes are metabolic organelles present in virtually all eukaryotic cells. They contain enzymes involved in hydrogen peroxide-based respiration and lipid metabolism. At present, only a small number of peroxisomal enzymes that are associated with oxidative stress response and metabolic disorders have been characterised biochemically. Therefore, we devised a sequence-based, multistep knowledge discovery strategy to identify potential novel peroxisomal protein candidates in small rodent model organisms and human.

METHODS

Screening of 130,629 putative translations of GenBank rodent and primate mRNA sequences was limited to the classical type-1 peroxisomal targeting signal [SA]-K-L. This motif is over-represented among peroxisomal proteins and has a high targeting efficiency. Subsequent steps of identifying co-occurring motifs, secondary structure properties, orthologues and variants, in combination with literature searching and visual inspection by domain experts, aimed at reduction of both false positive and negative validation targets.

RESULTS

Our method yielded 117 known peroxisome-targeted proteins and 29 novel candidate proteins. Of special interest were the mouse C530046K17Rik and 1300019N10Rik protein sequences that contain domains associated with enzymatic functions. C530046K17Rik showed no similarity to any known sequence of the animal kingdom, but weak similarity to the possible Leishmania quinone oxidoreductase and a putative cyanobacterium nicotinamide adenine dinucleotide phosphate (NADP)-dependent oxidoreductase. 1300019N10Rik contains two protease-related domains, glutamyl endopeptidase I and trypsin-like serine and cysteine proteases, which may have unique specificities to achieve efficient breakdown of proteins in the peroxisomes.

CONCLUSION

One mouse C57BL/6J strain-specific isocitrate dehydrogenase 1 isoform might be suitable to investigate potential phenotypes associated with the deficit of the intraperoxisomal reduced form of NADP (NADPH) and 2-oxoglutarate. Our biological knowledge discovery strategy enabled not only the identification of peroxisomal enzymes already described in the literature, but also the prediction of several novel proteins with possible roles in peroxisomal biochemistry and metabolism that are currently under experimental validation.

Hansenula polymorpha and Saccharomyces cerevisiae Pex5p's recognize different, independent peroxisomal targeting signals in alcohol oxidase

Peroxisomal alcohol oxidase (AO) from Hansenula polymorpha is inactive and partially mislocalized to the cytosol upon synthesis in Saccharomyces cerevisiae. Co-production with H. polymorpha pyruvate carboxylase (HpPyc1p) resulted in AO activation, but did not improve import into peroxisomes. We show that import of AO mediated by S. cerevisiae Pex5p is strictly dependent on the peroxisomal targeting signal 1 (PTS1) of AO and independent of HpPyc1p. In contrast, HpPex5p-mediated sorting of AO into S. cerevisiae peroxisomes is independent of the PTS1, but requires an alternative PTS that is only formed when HpPyc1p is co-produced and most likely involves folding and co-factor binding to AO.

6-Phosphofructo-2-kinase and fructose-2,6-bisphosphatase in Trypanosomatidae. Molecular characterization, database searches, modelling studies and evolutionary analysis

Fructose 2,6-bisphosphate is a potent allosteric activator of trypanosomatid pyruvate kinase and thus represents an important regulator of energy metabolism in these protozoan parasites. A 6-phosphofructo-2-kinase, responsible for the synthesis of this regulator, was highly purified from the bloodstream form of Trypanosoma brucei and kinetically characterized. By searching trypanosomatid genome databases, four genes encoding proteins homologous to the mammalian bifunctional enzyme 6-phosphofructo-2-kinase/fructose-2,6-bisphosphatase (PFK-2/FBPase-2) were found for both T. brucei and the related parasite Leishmania major and four pairs in Trypanosoma cruzi. These genes were predicted to each encode a protein in which, at most, only a single domain would be active. Two of the T. brucei proteins showed most conservation in the PFK-2 domain, although one of them was predicted to be inactive due to substitution of residues responsible for ligating the catalytically essential divalent metal cation; the two other proteins were most conserved in the FBPase-2 domain. The two PFK-2-like proteins were expressed in Escherichia coli. Indeed, the first displayed PFK-2 activity with similar kinetic properties to that of the enzyme purified from T. brucei, whereas no activity was found for the second. Interestingly, several of the predicted trypanosomatid PFK-2/FBPase-2 proteins have long N-terminal extensions. The N-terminal domains of the two polypeptides with most similarity to mammalian PFK-2s contain a series of tandem repeat ankyrin motifs. In other proteins such motifs are known to mediate protein-protein interactions. Phylogenetic analysis suggests that the four different PFK-2/FBPase-2 isoenzymes found in Trypanosoma and Leishmania evolved from a single ancestral bifunctional enzyme within the trypanosomatid lineage. A possible explanation for the evolution of multiple monofunctional enzymes and for the presence of the ankyrin-motif repeats in the PFK-2 isoenzymes is presented.

Peroxisomal Import of Human Alanine:glyoxylate Aminotransferase Requires Ancillary Targeting Information Remote from Its C Terminus

Although human alanine:glyoxylate aminotransferase (AGT) is imported into peroxisomes by a Pex5p-dependent pathway, the properties of its C-terminal tripeptide (KKL) are unlike those of any other type 1 peroxisomal targeting sequence (PTS1). We have previously suggested that AGT might possess ancillary targeting information that enables its unusual PTS1 to work. In this study, we have attempted to locate this information and to determine whether or not it is a characteristic of all vertebrate AGTs. Using the two-hybrid system, we show that human AGT interacts with human Pex5p in mammalian cells, but not yeast cells. Using (immuno)fluorescence microscopic analysis of the distribution of various constructs expressed in COS cells, we show the following. 1) The putative ancillary peroxisomal targeting information (PTS1A) in human AGT is located entirely within the smaller C-terminal structural domain of 110 amino acids, with the sequence between Val-324 and Ile-345 being the most likely candidate region. 2) The PTS1A is present in all mammalian AGTs studied (human, rat, guinea pig, rabbit, and cat), but not amphibian AGT (Xenopus). 3) The PTS1A is necessary for peroxisomal import of human, rabbit, and cat AGTs, but not rat and guinea pig AGTs. We speculate that the internal PTS1A of human AGT works in concert with the C-terminal PTS1 by interacting with Pex5p indirectly with the aid of a yet-to-be-identified mammal-specific adaptor molecule. This interaction might reshape the tetratricopeptide repeat domain allosterically, enabling it to accept KKL as a functional PTS1.

Refinement and prediction of protein prenylation motifs

Three prenylation motif predictors are presented that allow discrimination between proteins that are unique substrates of farnesyltransferase (FT) and those that can be alternatively processed by geranylgeranyltransferase I (GGT1).

We refined the motifs for carboxy-terminal protein prenylation by analysis of known substrates for farnesyltransferase (FT), geranylgeranyltransferase I (GGT1) and geranylgeranyltransferase II (GGT2). In addition to the CaaX box for the first two enzymes, we identify a preceding linker region that appears constrained in physicochemical properties, requiring small or flexible, preferably hydrophilic, amino acids. Predictors were constructed on the basis of sequence and physical property profiles, including interpositional correlations, and are available as the Prenylation Prediction Suite (PrePS, ) which also allows evaluation of evolutionary motif conservation. PrePS can predict partially overlapping substrate specificities, which is of medical importance in the case of understanding cellular action of FT inhibitors as anticancer and anti-parasite agents.

Large-scale analysis of human alternative protein isoforms: pattern classification and correlation with subcellular localization signals

We investigated human alternative protein isoforms of >2600 genes based on full-length cDNA clones and SwissProt. We classified the isoforms and examined their co-occurrence for each gene. Further, we investigated potential relationships between these changes and differential subcellular localization. The two most abundant patterns were the one with different C-terminal regions and the one with an internal insertion, which together account for 43% of the total. Although changes of the N-terminal region are less common than those of the C-terminal region, extension of the C-terminal region is much less common than that of the N-terminal region, probably because of the difficulty of removing stop codons in one isoform. We also found that there are some frequently used combinations of co-occurrence in alternative isoforms. We interpret this as evidence that there is some structural relationship which produces a repertoire of isoformal patterns. Finally, many terminal changes are predicted to cause differential subcellular localization, especially in targeting either peroxisomes or mitochondria. Our study sheds new light on the enrichment of the human proteome through alternative splicing and related events. Our database of alternative protein isoforms is available through the internet.

Biogenesis of peroxisomes and glycosomes: trypanosomatid glycosome assembly is a promising new drug target

In trypanosomatids (Trypanosoma and Leishmania), protozoa responsible for serious diseases of mankind in tropical and subtropical countries, core carbohydrate metabolism including glycolysis is compartmentalized in peculiar peroxisomes called glycosomes. Proper biogenesis of these organelles and the correct sequestering of glycolytic enzymes are essential to these parasites. Biogenesis of glycosomes in trypanosomatids and that of peroxisomes in other eukaryotes, including the human host, occur via homologous processes involving proteins called peroxins, which exert their function through multiple, transient interactions with each other. Decreased expression of peroxins leads to death of trypanosomes. Peroxins show only a low level of sequence conservation. Therefore, it seems feasible to design compounds that will prevent interactions of proteins involved in biogenesis of trypanosomatid glycosomes without interfering with peroxisome formation in the human host cells. Such compounds would be suitable as lead drugs against trypanosomatid-borne diseases.

Plant peroxisomes

Peroxisomes, one of single membrane-bound organelles, are present ubiquitously in eukaryotic cells. They were originally identified as organelles for production of hydrogen peroxide, the degradation of its hydrogen peroxide, and metabolism of fatty acids, which are functions common to almost all the organisms. Meanwhile, photorespiration and assimilation of symbiotically induced nitrogen are plant-specific functions. Recent postgenetic approaches such as transcriptome and proteome showed that plant peroxisomes are differentiated in various tissues, and revealed that peroxisomes have more important roles in various metabolic processes including biosynthesis of plant hormones than we speculated. All peroxisomal proteins, including metabolic enzymes in the matrix, membrane proteins, and factors responsible for peroxisome biogenesis, are nuclear encoded, and are provided from the outside of peroxisomes. Peroxisome biogenesis, such as protein transport, division, and enlargement, requires various complicated steps and is one of the most intriguing topics. Analyses using peroxisome biogenesis mutants and the whole-scale sequencing projects among several organisms revealed the existence of essential factors responsible for peroxisome biogenesis such as peroxins. This review addresses a comprehensive issue relating to function and biogenesis of plant peroxisomes and Arabidopsis mutants that have been accelerating our understanding of peroxisomes in planta.

Functional similarity between the peroxisomal PTS2 receptor binding protein Pex18p and the N-terminal half of the PTS1 receptor Pex5p

Within the extended receptor cycle of peroxisomal matrix import, the function of the import receptor Pex5p comprises cargo recognition and transport. While the C-terminal half (Pex5p-C) is responsible for PTS1 binding, the contribution of the N-terminal half of Pex5p (Pex5p-N) to the receptor cycle has been less clear. Here we demonstrate, using different techniques, that in Saccharomyces cerevisiae Pex5p-N alone facilitates the import of the major matrix protein Fox1p. This finding suggests that Pex5p-N is sufficient for receptor docking and cargo transport into peroxisomes. Moreover, we found that Pex5p-N can be functionally replaced by Pex18p, one of two auxiliary proteins of the PTS2 import pathway. A chimeric protein consisting of Pex18p (without its Pex7p binding site) fused to Pex5p-C is able to partially restore PTS1 protein import in a PEX5 deletion strain. On the basis of these results, we propose that the auxiliary proteins of the PTS2 import pathway fulfill roles similar to those of the N-terminal half of Pex5p in the PTS1 import pathway.

Hidden localization motifs: naturally occurring peroxisomal targeting signals in non-peroxisomal proteins

Functional but silent peroxisomal targeting signals have been found in non- peroxisomal proteins. This discovery has important implications for sequence-based signal prediction and for evolution.

Background

Can sequence segments coding for subcellular targeting or for posttranslational modifications occur in proteins that are not substrates in either of these processes? Although considerable effort has been invested in achieving low false-positive prediction rates, even accurate sequence-analysis tools for the recognition of these motifs generate a small but noticeable number of protein hits that lack the appropriate biological context but cannot be rationalized as false positives.

Results

We show that the carboxyl termini of a set of definitely non-peroxisomal proteins with predicted peroxisomal targeting signals interact with the peroxisomal matrix protein receptor peroxin 5 (PEX5) in a yeast two-hybrid test. Moreover, we show that examples of these proteins - chicken lysozyme, human tyrosinase and the yeast mitochondrial ribosomal protein L2 (encoded by MRP7) - are imported into peroxisomes in vivo if their original sorting signals are disguised. We also show that even prokaryotic proteins can contain peroxisomal targeting sequences.

Conclusions

Thus, functional localization signals can evolve in unrelated protein sequences as a result of neutral mutations, and subcellular targeting is hierarchically organized, with signal accessibility playing a decisive role. The occurrence of silent functional motifs in unrelated proteins is important for the development of sequence-based function prediction tools and the interpretation of their results. Silent functional signals have the potential to acquire importance in future evolutionary scenarios and in pathological conditions.

The Arabidopsis peroxisomal targeting signal type 2 receptor PEX7 is necessary for peroxisome function and dependent on PEX5

Plant peroxisomal proteins catalyze key metabolic reactions. Several peroxisome biogenesis PEROXIN (PEX) genes encode proteins acting in the import of targeted proteins necessary for these processes into the peroxisomal matrix. Most peroxisomal matrix proteins bear characterized Peroxisomal Targeting Signals (PTS1 or PTS2), which are bound by the receptors PEX5 or PEX7, respectively, for import into peroxisomes. Here we describe the isolation and characterization of an Arabidopsis peroxin mutant, pex7-1, which displays peroxisome-defective phenotypes including reduced PTS2 protein import. We also demonstrate that the pex5-1 PTS1 receptor mutant, which contains a lesion in a domain conserved among PEX7-binding proteins from various organisms, is defective not in PTS1 protein import, but rather in PTS2 protein import. Combining these mutations in a pex7-1 pex5-1 double mutant abolishes detectable PTS2 protein import and yields seedlings that are entirely sucrose-dependent for establishment, suggesting a severe block in peroxisomal fatty acid beta-oxidation. Adult pex7-1 pex5-1 plants have reduced stature and bear abnormally shaped seeds, few of which are viable. The pex7-1 pex5-1 seedlings that germinate have dramatically fewer lateral roots and often display fused cotyledons, phenotypes associated with reduced auxin response. Thus PTS2-directed peroxisomal import is necessary for normal embryonic development, seedling establishment, and vegetative growth.

A new definition for the consensus sequence of the peroxisome targeting signal type 2

All organisms except the nematode Caenorhabditis elegans have been shown to possess an import system for peroxisomal proteins containing a peroxisome targeting signal type 2 (PTS2). The currently accepted consensus sequence for this amino-terminal nonapeptide is -(R/K)(L/V/I)X(5)(H/Q)(L/A)-. Some C.elegans proteins contain putative PTS2 motifs, including the ortholog (CeMeK) of human mevalonate kinase, an enzyme known to be targeted by PTS2 to mammalian peroxisomes. We cloned the gene for CeMeK (open reading frame Y42G9A.4) and examined the subcellular localization of CeMeK and of two other proteins with putative PTS2s at their amino termini encoded by the open reading frames D1053.2 and W10G11.11. All three proteins localized to the cytosol, confirming and extending the finding that C.elegans lacks PTS2-dependent peroxisomal protein import. The putative PTS2s of the proteins encoded by D1053.2 and W10G11.11 did not function in targeting to peroxisomes in yeast or mammalian cells, suggesting that the current PTS2 consensus sequence is too broad. Analysis of available experimental data on both functional and nonfunctional PTS2s led to two re-evaluated PTS2 consensus sequences: -R(L/V/I/Q)XX(L/V/I/H)(L/S/G/A)X(H/Q)(L/A)-, describes the most common variants of PTS2, while -(R/K)(L/V/I/Q)XX(L/V/I/H/Q)(L/S/G/A/K)X(H/Q)(L/A/F)-, describes essentially all variants of PTS2. These redefined PTS2 consensus sequences will facilitate the identification of proteins of unknown cellular localization as possible peroxisomal proteins.

Penicillium chrysogenum Pex5p mediates differential sorting of PTS1 proteins to microbodies of the methylotrophic yeast Hansenula polymorpha

We have isolated the Penicillium chrysogenum pex5 gene encoding the receptor for microbody matrix proteins containing a type 1 peroxisomal targeting signal (PTS1). Pc-pex5 contains 2 introns and encodes a protein of approximately 75 kDa. P. chrysogenum pex5 disruptants appear to be highly unstable, show poor growth, and are unable to sporulate asexually. Furthermore, pex5 cells mislocalize a fluorescent PTS1 reporter protein to the cytosol. Pc-pex5 was expressed in a PEX5 null mutant of the yeast Hansenula polymorpha. Detailed analysis demonstrated that the PTS1 proteins dihydroxyacetone synthase and catalase were almost fully imported into microbodies. Surprisingly, alcohol oxidase, which also depends on Pex5p for import into microbodies, remained mainly in the cytosol. Thus, P. chrysogenum Pex5p has a different specificity of cargo recognition than its H. polymorpha counterpart. This was also suggested by the observation that Pc-Pex5p sorted a reporter protein fused to various functional PTS1 signals with different efficiencies.

AraPerox. A database of putative Arabidopsis proteins from plant peroxisomes

To identify unknown proteins from plant peroxisomes, the Arabidopsis genome was screened for proteins with putative major or minor peroxisome targeting signals type 1 or 2 (PTS1 or PTS2), as defined previously (Reumann S [2004] Plant Physiol 135: 783-800). About 220 and 60 proteins were identified that carry a putative PTS1 or PTS2, respectively. To further support postulated targeting to peroxisomes, several prediction programs were applied and the putative targeting domains analyzed for properties conserved in peroxisomal proteins and for PTS conservation in homologous plant expressed sequence tags. The majority of proteins with a major PTS and medium to high overall probability of peroxisomal targeting represent novel nonhypothetical proteins and include several enzymes involved in beta-oxidation of unsaturated fatty acids and branched amino acids, and 2-hydroxy acid oxidases with a predicted function in fatty acid alpha-oxidation, as well as NADP-dependent dehydrogenases and reductases. In addition, large protein families with many putative peroxisomal isoforms were recognized, including acyl-activating enzymes, GDSL lipases, and small thioesterases. Several proteins are homologous to prokaryotic enzymes of a novel aerobic hybrid degradation pathway for aromatic compounds and proposed to be involved in peroxisomal biosynthesis of plant hormones like jasmonic acid, auxin, and salicylic acid. Putative regulatory proteins of plant peroxisomes include protein kinases, small heat shock proteins, and proteases. The information on subcellular targeting prediction, homology, and in silico expression analysis for these Arabidopsis proteins has been compiled in the public database AraPerox to accelerate discovery and experimental investigation of novel metabolic and regulatory pathways of plant peroxisomes.

Pex5p binding affinities for canonical and noncanonical PTS1 peptides

The majority of proteins targeted to the peroxisomal lumen contain a C-terminal peroxisomal targeting signal-1 (PTS1) that is bound by the peroxin Pex5p. The PTS1 is generally regarded as a C-terminal tripeptide that adheres to the consensus (S/A/C)(K/R/H)(L/M). Previously, we studied the binding affinity of peptides of the form YQX(-3)X(-2)X(-1) to the peptide-binding domain of human Pex5p (referred to as Pex5p-C). Optimal affinity was found for YQSKL, which bound with an affinity of 200 +/- 40 nM. To extend this work, we investigated the properties of a peptide containing the last 9 residues of acyl-CoA oxidase (RHYLKPLQSKL) and discovered that it binds to Pex5p-C with a dissociation constant of 1.4 +/- 0.4 nM, 180 times tighter than YQSKL. Further analysis revealed that the enhanced affinity is primarily due to the presence of leucine in the (-5) position. In addition, a peptide corresponding to the luciferase C-terminus (YKGGKSKL) was found to bind Pex5p-C about 20 times tighter than YQSKL. The majority of this effect results from having lysine in position (-4). Catalase contains a noncanonical PTS1 (-AREKANL). The affinity of YQANL was found to be 3600 +/- 400 nM. This relatively weak binding is consistent with previous unsuccessful attempts to direct chloramphenicol acetyltransferase to the peroxisome by fusing -ANL to its C-terminus (-GGA-ANL). The peptides YKANL, YEKANL, YREKANL, and YAREKANL all bound Pex5p-C with higher affinities than did YQANL, but the affinities are still lower than peptides that correspond to functional targeting signals in other contexts. Because both catalase and Pex5p are tetramers (as opposed to the monomeric Pex5p-C and the peptides used in our studies), multidentate effects on binding affinity between Pex5p and other oligomeric proteins should be considered. Our study provides direct thermodynamic data revealing that peptide binding to Pex5p-C binding is favored by lysine in the (-4) position and leucine in the (-5) position. Our results suggest that peptides or proteins with optimized residues in the (-4) and/or (-5) positions can bind to Pex5p with affinities that are at least two orders of magnitude greater than that of YQSKL, and that this stabilization can compensates for otherwise weakly binding PTS1s.

Specification of the peroxisome targeting signals type 1 and type 2 of plant peroxisomes by bioinformatics analyses

To specify the C-terminal peroxisome targeting signal type 1 (PTS1) and the N-terminal PTS2 for higher plants, a maximum number of plant cDNAs and expressed sequence tags that are homologous to PTS1- and PTS2-targeted plant proteins was retrieved from the public databases and the primary structure of their targeting domains was analyzed for conserved properties. According to their high overall frequency in the homologs and their widespread occurence in different orthologous groups, nine major PTS1 tripeptides ([SA][RK][LM]> without AKM> plus SRI> and PRL>) and two major PTS2 nonapeptides (R[LI]x5HL) were defined that are considered good indicators for peroxisomal localization if present in unknown proteins. A lower but significant number of homologs contained 1 of 11 minor PTS1 tripeptides or of 9 minor PTS2 nonapeptides, many of which have not been identified before in plant peroxisomal proteins. The region adjacent to the PTS peptides was characterized by specific conserved properties as well, such as a pronounced incidence of basic and Pro residues and a high positive net charge, which probably play an auxiliary role in peroxisomal targeting. By contrast, several peptides with assumed peroxisomal targeting properties were not found in any of the 550 homologs and hence play--if at all--only a minor role in peroxisomal targeting. Based on the definition of these major and minor PTS and on the recognition of additional conserved properties, the accuracy of predicting peroxisomal proteins can be raised and plant genomes can be screened for novel proteins of peroxisomes more successfully.

Unfolding intermediate in the peroxisomal flavoprotein D-amino acid oxidase

The flavoenzyme d-amino acid oxidase (DAAO) from Rhodotorula gracilis is a peroxisomal enzyme and a prototypical member of the glutathione reductase family of flavoproteins. DAAO is a stable homodimer with a FAD molecule tightly bound to each 40-kDa subunit. In this work, the urea-induced unfolding of dimeric DAAO was compared with that of a monomeric form of the same protein, a deleted dimerization loop mutant. By using circular dichroism spectroscopy, protein and flavin fluorescence, 1,8-anilinonaphtalene sulfonic acid binding and activity assays, we demonstrated that the urea-induced unfolding of DAAO is a three-state process, yielding an intermediate, and that this process is reversible. The intermediate species lacks the catalytic activity and the characteristic tertiary structure of native DAAO but has significant secondary structure and retains flavin binding. Unfolding of DAAO proceeds through formation of an expanded, partially unfolded inactive intermediate, characterized by low solubility, by increased exposure of hydrophobic surfaces, and by increased sensitivity to trypsin of the beta-strand F5 belonging to the FAD binding domain. The oligomeric state does not modify the inferred folding process. The strand F5 is in contact with the C-terminal alpha-helix containing the Ser-Lys-Leu sequence corresponding to the type 1 peroxisomal targeting signal, and this structural element interacts with the N-terminal betaalphabeta flavin binding motif (Rossmann fold). The expanded conformation of the folding intermediate (and in particular the higher disorder of the mentioned secondary structure elements) could match the structure of the inactive holoenzyme required for in vivo trafficking of DAAO through the peroxisomal membrane.

Multiple targeting modules on peroxisomal proteins are not redundant: discrete functions of targeting signals within Pmp47 and Pex8p

Several peroxisomal proteins have two nonoverlapping targeting signals. These signals have been termed "redundant" because targeting can still occur with only one signal. We now report that separate targeting motifs within both Pmp47 and Pex8 provide complementary function. Pmp47 is an ATP translocator that contains six transmembrane domains (TMDs). We had previously shown that the TMD2 region (termed TMD2R, consisting of TMD2 and a short adjacent segment of cytosolic loop) was required for targeting to proliferated peroxisomes in Saccharomyces cerevisiae. We now report that the analogous TMD4R, which cannot target to proliferated peroxisomes, targets at least as well, or much better (depending on strain and growth conditions) in cells containing only basal (i.e., nonproliferated) peroxisomes. These data suggest differences in the targeting pathway among peroxisome populations. Pex8p, a peripheral protein facing the matrix, contains a typical carboxy terminal targeting sequence (PTS1) that has been shown to be nonessential for targeting, indicating the existence of a second targeting domain (not yet defined in S. cerevisiae); thus, its function was unknown. We show that targeting to basal peroxisomes, but not to proliferated peroxisomes, is more efficient with the PTS1 than without it. Our results indicate that multiple targeting signals within peroxisomal proteins extend coverage among heterogeneous populations of peroxisomes and increase efficiency of targeting in some metabolic states.

TPR subunits of the anaphase-promoting complex mediate binding to the activator protein CDH1

BACKGROUND

Chromosome segregation and mitotic exit depend on activation of the anaphase-promoting complex (APC) by the substrate adaptor proteins CDC20 and CDH1. The APC is a ubiquitin ligase composed of at least 11 subunits. The interaction of APC2 and APC11 with E2 enzymes is sufficient for ubiquitination reactions, but the functions of most other subunits are unknown.

RESULTS

We have biochemically characterized subcomplexes of the human APC. One subcomplex, containing APC2/11, APC1, APC4, and APC5, can assemble multiubiquitin chains but is unable to bind CDH1 and to ubiquitinate substrates. The other subcomplex contains all known APC subunits except APC2/11. This subcomplex can recruit CDH1 but fails to support any ubiquitination reaction. In vitro, the C termini of CDC20 and CDH1 bind to the closely related TPR subunits APC3 and APC7. Homology modeling predicts that these proteins are similar in structure to the peroxisomal import receptor PEX5, which binds cargo proteins via their C termini. APC activation by CDH1 depends on a conserved C-terminal motif that is also found in CDC20 and APC10.

CONCLUSIONS

APC1, APC4, and APC5 may connect APC2/11 with TPR subunits. TPR domains in APC3 and APC7 recruit CDH1 to the APC and may thereby bring substrates into close proximity of APC2/11 and E2 enzymes. In analogy to PEX5, the different TPR subunits of the APC might function as receptors that interact with the C termini of regulatory proteins such as CDH1, CDC20, and APC10.

Prediction of lipid posttranslational modifications and localization signals from protein sequences: big-Pi, NMT and PTS1

Many posttranslational modifications (N-myristoylation or glycosylphosphatidylinositol (GPI) lipid anchoring) and localization signals (the peroxisomal targeting signal PTS1) are encoded in short, partly compositionally biased regions at the N- or C-terminus of the protein sequence. These sequence signals are not well defined in terms of amino acid type preferences but they have significant interpositional correlations. Although the number of verified protein examples is small, the quantification of several physical conditions necessary for productive protein binding with the enzyme complexes executing the respective transformations can lead to predictors that recognize the signals from the amino acid sequence of queries alone. Taxon-specific prediction functions are required due to the divergent evolution of the active complexes. The big-Pi tool for the prediction of the C-terminal signal for GPI lipid anchor attachment is available for metazoan, protozoan and plant sequences. The myristoyl transferase (NMT) predictor recognizes glycine N-myristoylation sites (at the N-terminus and for fragments after processing) of higher eukaryotes (including their viruses) and fungi. The PTS1 signal predictor finds proteins with a C-terminus appropriate for peroxisomal import (for metazoa and fungi). Guidelines for application of the three WWW-based predictors (http://mendel.imp.univie.ac.at/) and for the interpretation of their output are described.

Prediction of peroxisomal targeting signal 1 containing proteins from amino acid sequence

Peroxisomal matrix proteins have to be imported into their target organelle post-translationally. The major translocation pathway depends on a C-terminal targeting signal, termed PTS1. Our previous analysis of sequence variability in the PTS1 motif revealed that, in addition to the known C-terminal tripeptide, at least nine residues directly upstream are important for signal recognition in the PTS1-Pex5 receptor complex. The refined PTS1 motif description was implemented in a prediction tool composed of taxon-specific functions (metazoa, fungi, remaining taxa), capable of recognising potential PTS1s in query sequences. The composite score function consists of classical profile terms and additional terms penalising deviations from the derived physical property pattern over sequence segments. The prediction algorithm has been validated with a self-consistency and three different cross-validation tests. Additionally, we tested the tool on a large set of non-peroxisomal negatives, on mutation data, and compared the prediction rate to the PTS1 component of the PSORT2 program. The sensitivity of our predictor in recognising documented PTS1 signal containing proteins is close to 90% for reliable prediction. The predictor distinguishes even SKL-appended non-peroxisomally targeted proteins such as a mouse dihydrofolate reductase-SKL construct. The corresponding rate of false positives is not worse than 0.8%; thus, the tool can be applied for large-scale unsupervised sequence database annotation. A scan of public protein databases uncovered a number of yet uncharacterised proteins for which the PTS1 signal might be critical for biological function. The predicted presence of a PTS1 signal implies peroxisomal localisation in the absence of N-terminal targeting sequences such as the mitochondrial import signal.