The eukaryotic linear motif resource ELM: 10 years and counting

Plasticity of an ultrafast interaction between nucleoporins and nuclear transport receptors

The mechanisms by which intrinsically disordered proteins engage in rapid and highly selective binding is a subject of considerable interest and represents a central paradigm to nuclear pore complex (NPC) function, where nuclear transport receptors (NTRs) move through the NPC by binding disordered phenylalanine-glycine-rich nucleoporins (FG-Nups). Combining single-molecule fluorescence, molecular simulations, and nuclear magnetic resonance, we show that a rapidly fluctuating FG-Nup populates an ensemble of conformations that are prone to bind NTRs with near diffusion-limited on rates, as shown by stopped-flow kinetic measurements. This is achieved using multiple, minimalistic, low-affinity binding motifs that are in rapid exchange when engaging with the NTR, allowing the FG-Nup to maintain an unexpectedly high plasticity in its bound state. We propose that these exceptional physical characteristics enable a rapid and specific transport mechanism in the physiological context, a notion supported by single molecule in-cell assays on intact NPCs.

•

Integrative structural biology reveals the basis of rapid nuclear transport

•

Transient binding of disordered nucleoporins leaves their plasticity unaffected

•

Multiple minimalistic low-affinity binding motifs create a polyvalent complex

•

A highly reactive and dynamic surface permits an ultrafast binding mechanism

Intrinsically disordered proteins: emerging interaction specialists

Intrinsically disordered proteins or regions of proteins (IDPs/IDRs) most often function through protein-protein interactions, when they permanently or transiently bind partner molecules with diverse functional consequences. There is a rapid advance in our understanding of the ensuing functional modes, obtained from describing atomic details of individual complexes, proteome-wide studies of interactomes and characterizing loosely assembled hydrogels and tightly packed amyloids. Here we briefly survey the most important recent methodological developments and structural-functional observations, with the aim of increasing the general appreciation of IDPs/IDRs as 'interaction specialists'.

Extensive identification and analysis of conserved small ORFs in animals

Background

There is increasing evidence that transcripts or transcript regions annotated as non-coding can harbor functional short open reading frames (sORFs). Loss-of-function experiments have identified essential developmental or physiological roles for a few of the encoded peptides (micropeptides), but genome-wide experimental or computational identification of functional sORFs remains challenging.

Results

Here, we expand our previously developed method and present results of an integrated computational pipeline for the identification of conserved sORFs in human, mouse, zebrafish, fruit fly, and the nematode C. elegans. Isolating specific conservation signatures indicative of purifying selection on amino acid (rather than nucleotide) sequence, we identify about 2,000 novel small ORFs located in the untranslated regions of canonical mRNAs or on transcripts annotated as non-coding. Predicted sORFs show stronger conservation signatures than those identified in previous studies and are sometimes conserved over large evolutionary distances. The encoded peptides have little homology to known proteins and are enriched in disordered regions and short linear interaction motifs. Published ribosome profiling data indicate translation of more than 100 novel sORFs, and mass spectrometry data provide evidence for more than 70 novel candidates.

Conclusions

Taken together, we identify hundreds of previously unknown conserved sORFs in major model organisms. Our computational analyses and integration with experimental data show that these sORFs are expressed, often translated, and sometimes widely conserved, in some cases even between vertebrates and invertebrates. We thus provide an integrated resource of putatively functional micropeptides for functional validation in vivo.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0742-x) contains supplementary material, which is available to authorized users.

Hammock: a hidden Markov model-based peptide clustering algorithm to identify protein-interaction consensus motifs in large datasets

Motivation: Proteins often recognize their interaction partners on the basis of short linear motifs located in disordered regions on proteins’ surface. Experimental techniques that study such motifs use short peptides to mimic the structural properties of interacting proteins. Continued development of these methods allows for large-scale screening, resulting in vast amounts of peptide sequences, potentially containing information on multiple protein-protein interactions. Processing of such datasets is a complex but essential task for large-scale studies investigating protein-protein interactions.

Results: The software tool presented in this article is able to rapidly identify multiple clusters of sequences carrying shared specificity motifs in massive datasets from various sources and generate multiple sequence alignments of identified clusters. The method was applied on a previously published smaller dataset containing distinct classes of ligands for SH3 domains, as well as on a new, an order of magnitude larger dataset containing epitopes for several monoclonal antibodies. The software successfully identified clusters of sequences mimicking epitopes of antibody targets, as well as secondary clusters revealing that the antibodies accept some deviations from original epitope sequences. Another test indicates that processing of even much larger datasets is computationally feasible.

Availability and implementation: Hammock is published under GNU GPL v. 3 license and is freely available as a standalone program (from http://www.recamo.cz/en/software/hammock-cluster-peptides/) or as a tool for the Galaxy toolbox (from https://toolshed.g2.bx.psu.edu/view/hammock/hammock). The source code can be downloaded from https://github.com/hammock-dev/hammock/releases.

Contact:muller@mou.cz

Supplementaryinformation:Supplementary data are available at Bioinformatics online.

The Arabidopsis mitogen-activated protein kinase 6 is associated with γ-tubulin on microtubules, phosphorylates EB1c and maintains spindle orientation under nitrosative stress

Stress-activated plant mitogen-activated protein (MAP) kinase pathways play roles in growth adaptation to the environment by modulating cell division through cytoskeletal regulation, but the mechanisms are poorly understood. We performed protein interaction and phosphorylation experiments with cytoskeletal proteins, mass spectrometric identification of MPK6 complexes and immunofluorescence analyses of the microtubular cytoskeleton of mitotic cells using wild-type, mpk6-2 mutant and plants overexpressing the MAP kinase-inactivating phosphatase, AP2C3. We showed that MPK6 interacted with γ-tubulin and co-sedimented with plant microtubules polymerized in vitro. It was the active form of MAP kinase that was enriched with microtubules and followed similar dynamics to γ-tubulin, moving from poles to midzone during the anaphase-to-telophase transition. We found a novel substrate for MPK6, the microtubule plus end protein, EB1c. The mpk6-2 mutant was sensitive to 3-nitro-l-tyrosine (NO2 -Tyr) treatment with respect to mitotic abnormalities, and root cells overexpressing AP2C3 showed defects in chromosome segregation and spindle orientation. Our data suggest that the active form of MAP kinase interacts with γ-tubulin on specific subsets of mitotic microtubules during late mitosis. MPK6 phosphorylates EB1c, but not EB1a, and has a role in maintaining regular planes of cell division under stress conditions.

Identifying and characterising key alternative splicing events in Drosophila development

BACKGROUND

In complex Metazoans a given gene frequently codes for multiple protein isoforms, through processes such as alternative splicing. Large scale functional annotation of these isoforms is a key challenge for functional genomics. This annotation gap is increasing with the large numbers of multi transcript genes being identified by technologies such as RNASeq. Furthermore attempts to characterise the functions of splicing in an organism are complicated by the difficulty in distinguishing functional isoforms from those produced by splicing errors or transcription noise. Tools to help prioritise candidate isoforms for testing are largely absent.

RESULTS

In this study we implement a Time-course Switch (TS) score for ranking isoforms by their likelihood of producing additional functions based on their developmental expression profiles, as reported by modENCODE. The TS score allows us to better investigate functional roles of different isoforms expressed in multi transcript genes. From this analysis, we find that isoforms with high TS scores have sequence feature changes consistent with more deterministic splicing and functional changes and tend to gain domains or whole exons which could carry additional functions. Furthermore these functions appear to be particularly important for essential regulatory roles, establishing functional isoform switching as key for regulatory processes. Based on the TS score we develop a Transcript Annotations Pipeline for Alternative Splicing (TAPAS) that identifies functional neighbourhoods of potentially interesting isoforms.

CONCLUSIONS

We have identified a subset of protein isoforms which appear to have high functional significance, particularly in regulation. This has been made possible through the development of novel methods that make use of transcript expression profiles. The methods and analyses we present here represent important first steps in the development of tools to address the near complete lack of isoform specific function annotation. In turn the tools allow us to better characterise the regulatory functions of alternative splicing in more detail.

MIEN1, a novel interactor of Annexin A2, promotes tumor cell migration by enhancing AnxA2 cell surface expression

Background

Migration and invasion enhancer 1 (MIEN1) is a novel gene found to be abundantly expressed in breast tumor tissues and functions as a critical regulator of tumor cell migration and invasion to promote systemic metastases. Previous studies have identified post-translational modifications by isoprenylation at the C-terminal tail of MIEN1 to favor its translocation to the inner leaflet of plasma membrane and its function as a membrane-bound adapter molecule. However, the exact molecular events at the membrane interface activating the MIEN1-driven tumor cell motility are vaguely understood.

Methods

MIEN1 was first studied using in-silico analysis on available RNA sequencing data of human breast tissues and its expression was ascertained in breast cells. We performed several assays including co-immunoprecipitation, wound healing, western blotting and immunofluorescence to decipher the molecular events involved in MIEN1-mediated tumor cell migration.

Results

Clinically, MIEN1 is predominantly overexpressed in Her-2 and luminal B subtypes of breast tumors, and its increased expression correlates with poor disease free survival. Molecular studies identified a phosphorylation-dependent activation signal in the immunoreceptor tyrosine based activation motif (ITAM) of MIEN1 and the phosphorylation-deficient MIEN1-mutants (Y39F/50 F) to regulate filopodia generation, migration and invasion. We found that ITAM-phosphorylation of MIEN1 is significantly impaired in isoprenylation-deficient MIEN1 mutants indicating that prenylation of MIEN1 and membrane association is required for cross-phosphorylation of tyrosine residues. Furthermore, we identified MIEN1 as a novel interactor of Annexin A2 (AnxA2), a Ca²⁺ -dependent phospholipid binding protein, which serves as an extracellular proteolytic center regulating plasmin generation. Fluorescence resonance energy transfer (FRET) confirmed that MIEN1 physically interacts with AnxA2 and functional studies revealed that they mutually cooperate to accentuate tumor cell motility. Interestingly, our study identified that ectopic overexpression of MIEN1 significantly enhances Tyr23-phosphorylation on AnxA2, thereby stimulating cell surface translocation of AnxA2 and catalyzing the activation of its proteolytic activity.

Conclusion

Our data show that the presence and interaction of both MIEN1 and AnxA2 in breast tumors are crucial drivers of cell motility. Our study has now deciphered a novel regulatory network governing the vicious process of breast tumor cell invasion-metastasis, and findings suggest MIEN1-AnxA2 as prospective targets to counter the deadly disease.

Electronic supplementary material

The online version of this article (doi:10.1186/s12943-015-0428-8) contains supplementary material, which is available to authorized users.

Physiological, pathological, and structural implications of non-enzymatic protein-protein interactions of the multifunctional human transglutaminase 2

Transglutaminase 2 (TG2) is a ubiquitously expressed member of an enzyme family catalyzing Ca(2+)-dependent transamidation of proteins. It is a multifunctional protein having several well-defined enzymatic (GTP binding and hydrolysis, protein disulfide isomerase, and protein kinase activities) and non-enzymatic (multiple interactions in protein scaffolds) functions. Unlike its enzymatic interactions, the significance of TG2's non-enzymatic regulation of its activities has recently gained importance. In this review, we summarize all the partners that directly interact with TG2 in a non-enzymatic manner and analyze how these interactions could modulate the crosslinking activity and cellular functions of TG2 in different cell compartments. We have found that TG2 mostly acts as a scaffold to bridge various proteins, leading to different functional outcomes. We have also studied how specific structural features, such as intrinsically disordered regions and embedded short linear motifs contribute to multifunctionality of TG2. Conformational diversity of intrinsically disordered regions enables them to interact with multiple partners, which can result in different biological outcomes. Indeed, ID regions in TG2 were identified in functionally relevant locations, indicating that they could facilitate conformational transitions towards the catalytically competent form. We reason that these structural features contribute to modulating the physiological and pathological functions of TG2 and could provide a new direction for detecting unique regulatory partners. Additionally, we have assembled all known anti-TG2 antibodies and have discussed their significance as a toolbox for identifying and confirming novel TG2 regulatory functions.

Computational and statistical methods for high-throughput analysis of post-translational modifications of proteins

The investigation of post-translational modifications (PTMs) represents one of the main research focuses for the study of protein function and cell signaling. Mass spectrometry instrumentation with increasing sensitivity improved protocols for PTM enrichment and recently established pipelines for high-throughput experiments allow large-scale identification and quantification of several PTM types. This review addresses the concurrently emerging challenges for the computational analysis of the resulting data and presents PTM-centered approaches for spectra identification, statistical analysis, multivariate analysis and data interpretation. We furthermore discuss the potential of future developments that will help to gain deep insight into the PTM-ome and its biological role in cells. This article is part of a Special Issue entitled: Computational Proteomics.

Discovering short linear protein motif based on selective training of profile hidden Markov models

Short linear motifs (SLiMs) in proteins are relatively conservative sequence patterns within disordered regions of proteins, typically 3-10 amino acids in length. They play an important role in mediating protein-protein interactions. Discovering SLiMs by computational methods has attracted more and more attention, most of which were based on regular expressions and profiles. In this paper, a de novo motif discovery method was proposed based on profile hidden Markov models (HMMs), which can not only provide the emission probabilities of amino acids in the defined positions of SLiMs, but also model the undefined positions. We adopted the ordered region masking and the relative local conservation (RLC) masking to improve the signal to noise ratio of the query sequences while applying evolutionary weighting to make the important sequences in evolutionary process get more attention by the selective training of profile HMMs. The experimental results show that our method and the profile-based method returned different subsets within a SLiMs dataset, and the performance of the two approaches are equivalent on a more realistic discovery dataset. Profile HMM-based motif discovery methods complement the existing methods and provide another way for SLiMs analysis.

Globular and disordered-the non-identical twins in protein-protein interactions

In biology proteins from different structural classes interact across and within classes in ways that are optimized to achieve balanced functional outputs. The interactions between intrinsically disordered proteins (IDPs) and other proteins rely on changes in flexibility and this is seen as a strong determinant for their function. This has fostered the notion that IDP's bind with low affinity but high specificity. Here we have analyzed available detailed thermodynamic data for protein-protein interactions to put to the test if the thermodynamic profiles of IDP interactions differ from those of other protein-protein interactions. We find that ordered proteins and the disordered ones act as non-identical twins operating by similar principles but where the disordered proteins complexes are on average less stable by 2.5 kcal mol(-1).

A Dynein Light Chain 1 Binding Motif in Rabies Virus Polymerase L Protein Plays a Role in Microtubule Reorganization and Viral Primary Transcription

Rabies virus (RABV) polymerase L together with phosphoprotein P forms the PL polymerase complex that is essential for replication and transcription. However, its exact mechanism of action, interactions with cellular factors, and intracellular distribution are yet to be understood. Here by imaging a fluorescently tagged polymerase (mCherry-RABV-L), we show that L accumulates at acetylated and reorganized microtubules (MT). In silico analysis revealed a dynein light chain 1 (DLC1) binding motif in L that could mediate MT binding through dynein motors. As DLC1 binding by polymerase cofactor P is known, we compared the impact of the DLC1-binding motifs in P and L. Viruses with mutations in the respective motifs revealed that both motifs are required for efficient primary transcription, indicating that DLC1 acts as a transcription enhancer by binding to both P and L. Notably, also the levels of cellular DLC1 protein were regulated by both motifs, suggesting regulation of the DLC1 gene expression by both P and L. Finally, disruption of the motif in L resulted in a cell-type-specific loss of MT localization, demonstrating that DLC1 is involved in L-mediated cytoskeleton reorganization. Overall, we conclude that DLC1 acts as a transcription factor that stimulates primary RABV transcription by binding to both P and L. We further conclude that L influences MT organization and posttranslational modification, suggesting a model in which MT manipulation by L contributes to efficient intracellular transport of virus components and thus may serve as an important step in virus replication.

IMPORTANCE

Regulation of rabies virus polymerase complex by viral and cellular factors thus far has not been fully understood. Although cellular dynein light chain 1 (DLC1) has been reported to increase primary transcription by binding to polymerase cofactor phosphoprotein P, the detailed mechanism is unknown, and it is also not known whether the large enzymatic polymerase subunit L is involved. By fluorescence microscopy analysis of fluorescence-tagged rabies virus L, in silico identification of a potential DLC1 binding site in L, and characterization of recombinant rabies virus mutants, we show that a DLC1 binding motif in L is involved in cytoskeleton localization and reorganization, primary transcription regulation by DLC1, and regulation of cellular DLC1 gene expression. By providing evidence for a direct contribution of a DLC1 binding motif in L, our data significantly increase the understanding of rabies virus polymerase regulation and host manipulation by the virus as well.

Gibson assembly: an easy way to clone potyviral full-length infectious cDNA clones expressing an ectopic VPg

BACKGROUND

Approaches to simplify and accelerate the construction of full-length infectious cDNA clones for plant potyviruses have been described, based on cloning strategies involving in vitro ligation or homologous recombination in yeast. In the present study, we developed a faster and more efficient in vitro recombination system using Gibson assembly (GA), to engineer a Lettuce mosaic virus (LMV) infectious clone expressing an ectopic mcherry-tagged VPg (Viral protein genome-linked) for in planta subcellular localization of the viral protein in an infection context.

METHODS

Three overlapping long distance PCR fragments were amplified and assembled in a single-step process based on in vitro recombination (Gibson assembly). The resulting 17.5 kbp recombinant plasmids (LMVmchVPg_Ec) were inoculated by biolistic on lettuce plants and then propagated mechanically on Nicotiana benthamiana. Confocal microscopy was used to analyze the subcellular localization of the ectopically expressed mcherry-VPg fusion protein.

RESULTS

The Gibson assembly allowed the cloning of the expected plasmids without any deletion. All the inoculated plants displayed symptoms characteristic of LMV infection. The majority of the mcherry fluorescent signal observed using confocal microscopy was located in the nucleus and nucleolus as expected for a potyviral VPg.

CONCLUSIONS

This is the first report of the use of the Gibson assembly method to construct full-length infectious cDNA clones of a potyvirus genome. This is also the first description of the ectopic expression of a tagged version of a potyviral VPg without affecting the viability of the recombinant potyvirus.

Conserved Functional Motifs and Homology Modeling to Predict Hidden Moonlighting Functional Sites

Moonlighting functional centers within proteins can provide them with hitherto unrecognized functions. Here, we review how hidden moonlighting functional centers, which we define as binding sites that have catalytic activity or regulate protein function in a novel manner, can be identified using targeted bioinformatic searches. Functional motifs used in such searches include amino acid residues that are conserved across species and many of which have been assigned functional roles based on experimental evidence. Molecules that were identified in this manner seeking cyclic mononucleotide cyclases in plants are used as examples. The strength of this computational approach is enhanced when good homology models can be developed to test the functionality of the predicted centers in silico, which, in turn, increases confidence in the ability of the identified candidates to perform the predicted functions. Computational characterization of moonlighting functional centers is not diagnostic for catalysis but serves as a rapid screening method, and highlights testable targets from a potentially large pool of candidates for subsequent in vitro and in vivo experiments required to confirm the functionality of the predicted moonlighting centers.

Looking below the surface of nicotinic acetylcholine receptors

The amino acid sequences of nicotinic acetylcholine receptors (nAChRs) from diverse species can be compared across extracellular, transmembrane, and intracellular domains. The intracellular domains are most divergent among subtypes, yet relatively consistent among species. The diversity indicates that each nAChR subtype has a unique language for communication with its host cell. The conservation across species also suggests that the intracellular domains have defining functional roles for each subtype. Secondary structure prediction indicates two relatively conserved alpha helices within the intracellular domains of all nAChRs. Among all subtypes, the intracellular domain of α7 nAChR is one of the most well conserved, and α7 nAChRs have effects in non-neuronal cells independent of generating ion currents, making it likely that the α7 intracellular domain directly mediates signal transduction. There are potential phosphorylation and protein-binding sites in the α7 intracellular domain, which are conserved and may be the basis for α7-mediated signal transduction.

The knottin-like Blufensin family regulates genes involved in nuclear import and the secretory pathway in barley-powdery mildew interactions

Plants have evolved complex regulatory mechanisms to control a multi-layered defense response to microbial attack. Both temporal and spatial gene expression are tightly regulated in response to pathogen ingress, modulating both positive and negative control of defense. BLUFENSINs, small knottin-like peptides in barley, wheat, and rice, are highly induced by attack from fungal pathogens, in particular, the obligate biotrophic fungus, Blumeria graminis f. sp. hordei (Bgh), causal agent of barley powdery mildew. Previous research indicated that Blufensin1 (Bln1) functions as a negative regulator of basal defense mechanisms. In the current report, we show that BLN1 and BLN2 can both be secreted to the apoplast and Barley stripe mosaic virus (BSMV)-mediated overexpression of Bln2 increases susceptibility of barley to Bgh. Bimolecular fluorescence complementation (BiFC) assays signify that BLN1 and BLN2 can interact with each other, and with calmodulin. We then used BSMV-induced gene silencing to knock down Bln1, followed by Barley1 GeneChip transcriptome analysis, to identify additional host genes influenced by Bln1. Analysis of differential expression revealed a gene set enriched for those encoding proteins annotated to nuclear import and the secretory pathway, particularly Importin α1-b and Sec61 γ subunits. Further functional analysis of these two affected genes showed that when silenced, they also reduced susceptibility to Bgh. Taken together, we postulate that Bln1 is co-opted by Bgh to facilitate transport of disease-related host proteins or effectors, influencing the establishment of Bgh compatibility on its barley host.

Design and evaluation of antimalarial peptides derived from prediction of short linear motifs in proteins related to erythrocyte invasion

The purpose of this study was to investigate the blood stage of the malaria causing parasite, Plasmodium falciparum, to predict potential protein interactions between the parasite merozoite and the host erythrocyte and design peptides that could interrupt these predicted interactions. We screened the P. falciparum and human proteomes for computationally predicted short linear motifs (SLiMs) in cytoplasmic portions of transmembrane proteins that could play roles in the invasion of the erythrocyte by the merozoite, an essential step in malarial pathogenesis. We tested thirteen peptides predicted to contain SLiMs, twelve of them palmitoylated to enhance membrane targeting, and found three that blocked parasite growth in culture by inhibiting the initiation of new infections in erythrocytes. Scrambled peptides for two of the most promising peptides suggested that their activity may be reflective of amino acid properties, in particular, positive charge. However, one peptide showed effects which were stronger than those of scrambled peptides. This was derived from human red blood cell glycophorin-B. We concluded that proteome-wide computational screening of the intracellular regions of both host and pathogen adhesion proteins provides potential lead peptides for the development of anti-malarial compounds.

Genomic and single nucleotide polymorphism analysis of infectious bronchitis coronavirus

Infectious bronchitis virus (IBV) is a Gammacoronavirus that causes a highly contagious respiratory disease in chickens. A QX-like strain was analysed by high-throughput Illumina sequencing and genetic variation across the entire viral genome was explored at the sub-consensus level by single nucleotide polymorphism (SNP) analysis. Thirteen open reading frames (ORFs) in the order 5'-UTR-1a-1ab-S-3a-3b-E-M-4b-4c-5a-5b-N-6b-3'UTR were predicted. The relative frequencies of missense: silent SNPs were calculated to obtain a comparative measure of variability in specific genes. The most variable ORFs in descending order were E, 3b, 5'UTR, N, 1a, S, 1ab, M, 4c, 5a, 6b. The E and 3b protein products play key roles in coronavirus virulence, and RNA folding demonstrated that the mutations in the 5'UTR did not alter the predicted secondary structure. The frequency of SNPs in the Spike (S) protein ORF of 0.67% was below the genomic average of 0.76%. Only three SNPS were identified in the S1 subunit, none of which were located in hypervariable region (HVR) 1 or HVR2. The S2 subunit was considerably more variable containing 87% of the polymorphisms detected across the entire S protein. The S2 subunit also contained a previously unreported multi-A insertion site and a stretch of four consecutive mutated amino acids, which mapped to the stalk region of the spike protein. Template-based protein structure modelling produced the first theoretical model of the IBV spike monomer. Given the lack of diversity observed at the sub-consensus level, the tenet that the HVRs in the S1 subunit are very tolerant of amino acid changes produced by genetic drift is questioned.

The FAD-dependent glycerol-3-phosphate dehydrogenase of Giardia duodenalis: an unconventional enzyme that interacts with the g14-3-3 and it is a target of the antitumoral compound NBDHEX

The flagellated protozoan Giardia duodenalis is a worldwide parasite causing giardiasis, an acute and chronic diarrheal disease. Metabolism in G. duodenalis has a limited complexity thus making metabolic enzymes ideal targets for drug development. However, only few metabolic pathways (i.e., carbohydrates) have been described so far. Recently, the parasite homolog of the mitochondrial-like glycerol-3-phosphate dehydrogenase (gG3PD) has been identified among the interactors of the g14-3-3 protein. G3PD is involved in glycolysis, electron transport, glycerophospholipids metabolism, and hyperosmotic stress response, and is emerging as promising target in tumor treatment. In this work, we demonstrate that gG3PD is a functional flavoenzyme able to convert glycerol-3-phosphate into dihydroxyacetone phosphate and that its activity and the intracellular glycerol level increase during encystation. Taking advantage of co-immunoprecipitation assays and deletion mutants, we provide evidence that gG3PD and g14-3-3 interact at the trophozoite stage, the intracellular localization of gG3PD is stage dependent and it partially co-localizes with mitosomes during cyst development. Finally, we demonstrate that the gG3PD activity is affected by the antitumoral compound 6-(7-nitro-2,1,3-benzoxadiazol-4-ylthio)hexanol, that results more effective in vitro at killing G. duodenalis trophozoites than the reference drug metronidazole. Overall, our results highlight the involvement of gG3PD in processes crucial for the parasite survival thus proposing this enzyme as target for novel antigiardial interventions.

Proteomic Validation of Transcript Isoforms, Including Those Assembled from RNA-Seq Data

Human proteome analysis now requires an understanding of protein isoforms. We recently published the PG Nexus pipeline, which facilitates high confidence validation of exons and splice junctions by integrating genomics and proteomics data. Here we comprehensively explore how RNA-seq transcriptomics data, and proteomic analysis of the same sample, can identify protein isoforms. RNA-seq data from human mesenchymal (hMSC) stem cells were analyzed with our new TranscriptCoder tool to generate a database of protein isoform sequences. MS/MS data from matching hMSC samples were then matched against the TranscriptCoder-derived database, along with Ensembl and the neXtProt database. Querying the TranscriptCoder-derived or Ensembl database could unambiguously identify ∼450 protein isoforms, with isoform-specific proteotypic peptides, including candidate hMSC-specific isoforms for the genes DPYSL2 and FXR1. Where isoform-specific peptides did not exist, groups of nonisoform-specific proteotypic peptides could specifically identify many isoforms. In both the above cases, isoforms will be detectable with targeted MS/MS assays. Unfortunately, our analysis also revealed that some isoforms will be difficult to identify unambiguously as they do not have peptides that are sufficiently distinguishing. We covisualize mRNA isoforms and peptides in a genome browser to illustrate the above situations. Mass spectrometry data is available via ProteomeXchange (PXD001449).

DisCons: a novel tool to quantify and classify evolutionary conservation of intrinsic protein disorder

Background

Analyzing the amino acid sequence of an intrinsically disordered protein (IDP) in an evolutionary context can yield novel insights on the functional role of disordered regions and sequence element(s). However, in the case of many IDPs, the lack of evolutionary conservation of the primary sequence can hamper the study of functionality, because the conservation of their disorder profile and ensuing function(s) may not appear in a traditional analysis of the evolutionary history of the protein.

Results

Here we present DisCons (Disorder Conservation), a novel pipelined tool that combines the quantification of sequence- and disorder conservation to classify disordered residue positions. According to this scheme, the most interesting categories (for functional purposes) are constrained disordered residues and flexible disordered residues. The former residues show conservation of both the sequence and the property of disorder and are associated mainly with specific binding functionalities (e.g., short, linear motifs, SLiMs), whereas the latter class correspond to segments where disorder as a feature is important for function as opposed to the identity of the underlying sequence (e.g., entropic chains and linkers). DisCons therefore helps with elucidating the function(s) arising from the disordered state by analyzing individual proteins as well as large-scale proteomics datasets.

Conclusions

DisCons is an openly accessible sequence analysis tool that identifies and highlights structurally disordered segments of proteins where the conformational flexibility is conserved across homologs, and therefore potentially functional. The tool is freely available both as a web application and as stand-alone source code hosted at http://pedb.vib.be/discons.

DisCons: a novel tool to quantify and classify evolutionary conservation of intrinsic protein disorder

BACKGROUND

Analyzing the amino acid sequence of an intrinsically disordered protein (IDP) in an evolutionary context can yield novel insights on the functional role of disordered regions and sequence element(s). However, in the case of many IDPs, the lack of evolutionary conservation of the primary sequence can hamper the study of functionality, because the conservation of their disorder profile and ensuing function(s) may not appear in a traditional analysis of the evolutionary history of the protein.

RESULTS

Here we present DisCons (Disorder Conservation), a novel pipelined tool that combines the quantification of sequence- and disorder conservation to classify disordered residue positions. According to this scheme, the most interesting categories (for functional purposes) are constrained disordered residues and flexible disordered residues. The former residues show conservation of both the sequence and the property of disorder and are associated mainly with specific binding functionalities (e.g., short, linear motifs, SLiMs), whereas the latter class correspond to segments where disorder as a feature is important for function as opposed to the identity of the underlying sequence (e.g., entropic chains and linkers). DisCons therefore helps with elucidating the function(s) arising from the disordered state by analyzing individual proteins as well as large-scale proteomics datasets.

CONCLUSIONS

DisCons is an openly accessible sequence analysis tool that identifies and highlights structurally disordered segments of proteins where the conformational flexibility is conserved across homologs, and therefore potentially functional. The tool is freely available both as a web application and as stand-alone source code hosted at http://pedb.vib.be/discons .

RBFOX2 protein domains and cellular activities

RBFOX2 (RNA-binding protein, Fox-1 homologue 2)/RBM9 (RNA-binding-motif protein 9)/RTA (repressor of tamoxifen action)/HNRBP2 (hexaribonucleotide-binding protein 2) encodes an RNA-binding protein involved in tissue specific alternative splicing regulation and steroid receptors transcriptional activity. Its ability to regulate specific splicing profiles depending on context has been related to different expression levels of the RBFOX2 protein itself and that of other splicing regulatory proteins involved in the shared modulation of specific genes splicing. However, this cannot be the sole explanation as to why RBFOX2 plays a widespread role in numerous cellular mechanisms from development to cell survival dependent on cell/tissue type. RBFOX2 isoforms with altered protein domains exist. In the present article, we describe the main RBFOX2 protein domains, their importance in the context of splicing and transcriptional regulation and we propose that RBFOX2 isoform distribution may play a fundamental role in RBFOX2-specific cellular effects.

SELPHI: correlation-based identification of kinase-associated networks from global phospho-proteomics data sets

While phospho-proteomics studies have shed light on the dynamics of cellular signaling, they mainly describe global effects and rarely explore mechanistic details, such as kinase/substrate relationships. Tools and databases, such as NetworKIN and PhosphoSitePlus, provide valuable regulatory details on signaling networks but rely on prior knowledge. They therefore provide limited information on less studied kinases and fewer unexpected relationships given that better studied signaling events can mask condition- or cell-specific ‘network wiring’.

SELPHI is a web-based tool providing in-depth analysis of phospho-proteomics data that is intuitive and accessible to non-bioinformatics experts. It uses correlation analysis of phospho-sites to extract kinase/phosphatase and phospho-peptide associations, and highlights the potential flow of signaling in the system under study. We illustrate SELPHI via analysis of phospho-proteomics data acquired in the presence of erlotinib—a tyrosine kinase inhibitor (TKI)—in cancer cells expressing TKI-resistant and -sensitive variants of the Epidermal Growth Factor Receptor. In this data set, SELPHI revealed information overlooked by the reporting study, including the known role of MET and EPHA2 kinases in conferring resistance to erlotinib in TKI sensitive strains. SELPHI can significantly enhance the analysis of phospho-proteomics data contributing to improved understanding of sample-specific signaling networks. SELPHI is freely available via http://llama.mshri.on.ca/SELPHI.

Predicting peptide-mediated interactions on a genome-wide scale

We describe a method to predict protein-protein interactions (PPIs) formed between structured domains and short peptide motifs. We take an integrative approach based on consensus patterns of known motifs in databases, structures of domain-motif complexes from the PDB and various sources of non-structural evidence. We combine this set of clues using a Bayesian classifier that reports the likelihood of an interaction and obtain significantly improved prediction performance when compared to individual sources of evidence and to previously reported algorithms. Our Bayesian approach was integrated into PrePPI, a structure-based PPI prediction method that, so far, has been limited to interactions formed between two structured domains. Around 80,000 new domain-motif mediated interactions were predicted, thus enhancing PrePPI’s coverage of the human protein interactome.

Complexes formed between a structured domain on one protein and an unstructured peptide on another are ubiquitous. However, they are often quite difficult to detect experimentally. The development of computational approaches to predict domain-motif interactions is therefore an important goal. We report a method to predict domain-motif interactions using a Bayesian approach to integrate evidence from a variety of sources, including three-dimensional structural and non-structural information. The method was applied to the entire human proteome and showed significant improvement over existing methods. The method was incorporated into PrePPI, a computational pipeline for the prediction of protein-protein interactions that relies heavily on structural information. Approximately 80,000 new interactions were detected. The new PrePPI database provides easy access to about 400,000 human protein-protein interactions and should thus constitute a valuable resource in a variety of biological applications including the characterization of molecular interaction networks and, more generally, in the study of interactions mediated by proteins in families that may not be extensively studied experimentally.

Discovery of a novel ligand that modulates the protein-protein interactions of the AAA+ superfamily oncoprotein reptin

Developing approaches to discover protein-protein interactions (PPIs) remains a fundamental challenge. A chemical biology platform is applied here to identify novel PPIs for the AAA+ superfamily oncoprotein reptin. An in silico screen coupled with chemical optimization provided Liddean, a nucleotide-mimetic which modulates reptin's oligomerization status, protein-binding activity and global conformation. Combinatorial peptide phage library screening of Liddean-bound reptin with next generation sequencing identified interaction motifs including a novel reptin docking site on the p53 tumor suppressor protein. Proximity ligation assays demonstrated that endogenous reptin forms a predominantly cytoplasmic complex with its paralog pontin in cancer cells and Liddean promotes a shift of this complex to the nucleus. An emerging view of PPIs in higher eukaryotes is that they occur through a striking diversity of linear peptide motifs. The discovery of a compound that alters reptin's protein interaction landscape potentially leads to novel avenues for therapeutic development.

Co-expressed Cyclin D variants cooperate to regulate proliferation of germline nuclei in a syncytium

The role of the G1-phase Cyclin D-CDK 4/6 regulatory module in linking germline stem cell (GSC) proliferation to nutrition is evolutionarily variable. In invertebrate Drosophila and C. elegans GSC models, G1 is nearly absent and Cyclin E is expressed throughout the cell cycle, whereas vertebrate spermatogonial stem cells have a distinct G1 and Cyclin D1 plays an important role in GSC renewal. In the invertebrate, chordate, Oikopleura, where germline nuclei proliferate asynchronously in a syncytium, we show a distinct G1-phase in which 2 Cyclin D variants are co-expressed. Cyclin Dd, present in both somatic endocycling cells and the germline, localized to germline nuclei during G1 before declining at G1/S. Cyclin Db, restricted to the germline, remained cytoplasmic, co-localizing in foci with the Cyclin-dependent Kinase Inhibitor, CKIa. These foci showed a preferential spatial distribution adjacent to syncytial germline nuclei at G1/S. During nutrient-restricted growth arrest, upregulated CKIa accumulated in arrested somatic endoreduplicative nuclei but did not do so in germline nuclei. In the latter context, Cyclin Dd levels gradually decreased. In contrast, the Cyclin Dbβ splice variant, lacking the Rb-interaction domain and phosphodegron, was specifically upregulated and the number of cytoplasmic foci containing this variant increased. This upregulation was dependent on stress response MAPK p38 signaling. We conclude that under favorable conditions, Cyclin Dbβ-CDK6 sequesters CKIa in the cytoplasm to cooperate with Cyclin Dd-CDK6 in promoting germline nuclear proliferation. Under nutrient-restriction, this sequestration function is enhanced to permit continued, though reduced, cycling of the germline during somatic growth arrest.

CD2v Interacts with Adaptor Protein AP-1 during African Swine Fever Infection

African swine fever virus (ASFV) CD2v protein is believed to be involved in virulence enhancement, viral hemadsorption, and pathogenesis, although the molecular mechanisms of the function of this viral protein are still not fully understood. Here we describe that CD2v localized around viral factories during ASFV infection, suggesting a role in the generation and/or dynamics of these viral structures and hence in disturbing cellular traffic. We show that CD2v targeted the regulatory trans-Golgi network (TGN) protein complex AP-1, a key element in cellular traffic. This interaction was disrupted by brefeldin A even though the location of CD2v around the viral factory remained unchanged. CD2v-AP-1 binding was independent of CD2v glycosylation and occurred on the carboxy-terminal part of CD2v, where a canonical di-Leu motif previously reported to mediate AP-1 binding in eukaryotic cells, was identified. This motif was shown to be functionally interchangeable with the di-Leu motif present in HIV-Nef protein in an AP-1 binding assay. However, we demonstrated that it was not involved either in CD2v cellular distribution or in CD2v-AP-1 binding. Taken together, these findings shed light on CD2v function during ASFV infection by identifying AP-1 as a cellular factor targeted by CD2v and hence elucidate the cellular pathways used by the virus to enhance infectivity.

Inhibitory activities of short linear motifs underlie Hox interactome specificity in vivo

Hox proteins are well-established developmental regulators that coordinate cell fate and morphogenesis throughout embryogenesis. In contrast, our knowledge of their specific molecular modes of action is limited to the interaction with few cofactors. Here, we show that Hox proteins are able to interact with a wide range of transcription factors in the live Drosophila embryo. In this context, specificity relies on a versatile usage of conserved short linear motifs (SLiMs), which, surprisingly, often restrains the interaction potential of Hox proteins. This novel buffering activity of SLiMs was observed in different tissues and found in Hox proteins from cnidarian to mouse species. Although these interactions remain to be analysed in the context of endogenous Hox regulatory activities, our observations challenge the traditional role assigned to SLiMs and provide an alternative concept to explain how Hox interactome specificity could be achieved during the embryonic development.

DOI:http://dx.doi.org/10.7554/eLife.06034.001

In all animals, it is important that cells are correctly organised into tissues and organs. This organisation starts in the embryo, and cells are instructed to perform different roles depending on their position within the body.

A family of proteins called the Hox proteins coordinates the organisation of the cells in the animal embryo by binding to and controlling the expression of specific genes. To properly control their target genes, Hox proteins need to interact with other proteins called transcription factors that can also bind to the genes. However, only a few of these transcription factors have been identified so far, and it is not clear how Hox proteins are able to interact with them.

Here, Baëza, Viala, Heim et al. identified several more transcription factors that can bind to the Hox proteins in fruit fly embryos. The experiments show that Hox proteins are able to bind to many transcription factors that are very different from each other. Baëza, Viala, Heim et al. also show that two short sections within the Hox proteins known as short linear motifs are important for controlling these interactions. A fly Hox protein that was missing these motifs was able to interact with new transcription factors. This inhibitory role was found in Hox proteins from mice and sea anemones, suggesting that these motifs may play the same role in all animals.

Baëza, Viala, Heim et al.'s findings challenge the traditional view of the role of the short linear motifs in interactions between proteins. Also, the findings provide an alternative explanation for how the Hox proteins are only able to interact with particular transcription factors in animal embryos. The next step will be to find out whether the inhibitory role of short linear motifs could more generally apply to many other protein families.

DOI:http://dx.doi.org/10.7554/eLife.06034.002

A modular toolkit to inhibit proline-rich motif-mediated protein-protein interactions

Small-molecule competitors of protein-protein interactions are urgently needed for functional analysis of large-scale genomics and proteomics data. Particularly abundant, yet so far undruggable, targets include domains specialized in recognizing proline-rich segments, including Src-homology 3 (SH3), WW, GYF, and Drosophila enabled (Ena)/vasodilator-stimulated phosphoprotein (VASP) homology 1 (EVH1) domains. Here, we present a modular strategy to obtain an extendable toolkit of chemical fragments (ProMs) designed to replace pairs of conserved prolines in recognition motifs. As proof-of-principle, we developed a small, selective, peptidomimetic inhibitor of Ena/VASP EVH1 domain interactions. Highly invasive MDA MB 231 breast-cancer cells treated with this ligand showed displacement of VASP from focal adhesions, as well as from the front of lamellipodia, and strongly reduced cell invasion. General applicability of our strategy is illustrated by the design of an ErbB4-derived ligand containing two ProM-1 fragments, targeting the yes-associated protein 1 (YAP1)-WW domain with a fivefold higher affinity.

Phosphorylation regulates the ubiquitin-independent degradation of yeast Pah1 phosphatidate phosphatase by the 20S proteasome

Saccharomyces cerevisiae Pah1 phosphatidate phosphatase, which catalyzes the conversion of phosphatidate to diacylglycerol for triacylglycerol synthesis and simultaneously controls phosphatidate levels for phospholipid synthesis, is subject to the proteasome-mediated degradation in the stationary phase of growth. In this study, we examined the mechanism for its degradation using purified Pah1 and isolated proteasomes. Pah1 expressed in S. cerevisiae or Escherichia coli was not degraded by the 26S proteasome, but by its catalytic 20S core particle, indicating that its degradation is ubiquitin-independent. The degradation of Pah1 by the 20S proteasome was dependent on time and proteasome concentration at the pH optimum of 7.0. The 20S proteasomal degradation was conserved for human lipin 1 phosphatidate phosphatase. The degradation analysis using Pah1 truncations and its fusion with GFP indicated that proteolysis initiates at the N- and C-terminal unfolded regions. The folded region of Pah1, in particular the haloacid dehalogenase-like domain containing the DIDGT catalytic sequence, was resistant to the proteasomal degradation. The structural change of Pah1, as reflected by electrophoretic mobility shift, occurs through its phosphorylation by Pho85-Pho80, and the phosphorylation sites are located within its N- and C-terminal unfolded regions. Phosphorylation of Pah1 by Pho85-Pho80 inhibited its degradation, extending its half-life by ∼2-fold. The dephosphorylation of endogenously phosphorylated Pah1 by the Nem1-Spo7 protein phosphatase, which is highly specific for the sites phosphorylated by Pho85-Pho80, stimulated the 20S proteasomal degradation and reduced its half-life by 2.6-fold. These results indicate that the proteolysis of Pah1 by the 20S proteasome is controlled by its phosphorylation state.

A single PXXP motif in the C-terminal region of srGAP3 mediates binding to multiple SH3 domains

The Slit-Robo GTPase-activating protein 3 (srGAP3) has been implicated in different critical aspects of neuronal development. These findings have mainly been based on the characterisation of the three conserved globular N-terminal domains, while the function of the C-terminal region (CTR) is still unknown. We show that this predicted unstructured region acts as an adaptor by binding to the endocytic proteins Amphiphysin, Endophilin-A2, Endophilin-A1, as well as the Ras signalling protein Grb2. All these interactions depend on a single proline-rich motif in the CTR and the Src-homology 3 domains of the binding partners. Via these interactions srGAP3 could link receptor signalling events to the endocytic machinery.

Characterization of the avian Trojan gene family reveals contrasting evolutionary constraints

"Trojan" is a leukocyte-specific, cell surface protein originally identified in the chicken. Its molecular function has been hypothesized to be related to anti-apoptosis and the proliferation of immune cells. The Trojan gene has been localized onto the Z sex chromosome. The adjacent two genes also show significant homology to Trojan, suggesting the existence of a novel gene/protein family. Here, we characterize this Trojan family, identify homologues in other species and predict evolutionary constraints on these genes. The two Trojan-related proteins in chicken were predicted as a receptor-type tyrosine phosphatase and a transmembrane protein, bearing a cytoplasmic immuno-receptor tyrosine-based activation motif. We identified the Trojan gene family in ten other bird species and found related genes in three reptiles and a fish species. The phylogenetic analysis of the homologues revealed a gradual diversification among the family members. Evolutionary analyzes of the avian genes predicted that the extracellular regions of the proteins have been subjected to positive selection. Such selection was possibly a response to evolving interacting partners or to pathogen challenges. We also observed an almost complete lack of intracellular positively selected sites, suggesting a conserved signaling mechanism of the molecules. Therefore, the contrasting patterns of selection likely correlate with the interaction and signaling potential of the molecules.

Phosphorylation of ORF1p is required for L1 retrotransposition

Although members of the L1 (LINE-1) clade of non-LTR retrotransposons can be deleterious, the L1 clade has remained active in most mammals for ∼100 million years and generated almost 40% of the human genome. The details of L1-host interaction are largely unknown, however. Here we report that L1 activity requires phosphorylation of the protein encoded by the L1 ORF1 (ORF1p). Critical phospho-acceptor residues (two serines and two threonines) reside in four conserved proline-directed protein kinase (PDPK) target sites. The PDPK family includes mitogen-activated protein kinases and cyclin-dependent kinases. Mutation of any PDPK phospho-acceptor inhibits L1 retrotransposition. The phosphomimetic aspartic acid can restore activity at the two serine sites, but not at either threonine site, where it is strongly inhibitory. ORF1p also contains conserved PDPK docking sites, which promote specific interaction of PDPKs with their targets. As expected, mutations in these sites also inhibit L1 activity. PDPK mutations in ORF1p that inactivate L1 have no significant effect on the ability of ORF1p to anneal RNA in vitro, an important biochemical property of the protein. We show that phosphorylated PDPK sites in ORF1p are required for an interaction with the peptidyl prolyl isomerase 1 (Pin1), a critical component of PDPK-mediated regulation. Pin1 acts via isomerization of proline side chains at phosphorylated PDPK motifs, thereby affecting substrate conformation and activity. Our demonstration that L1 activity is dependent on and integrated with cellular phosphorylation regulatory cascades significantly increases our understanding of interactions between L1 and its host.

QSLiMFinder: improved short linear motif prediction using specific query protein data

Motivation: The sensitivity of de novo short linear motif (SLiM) prediction is limited by the number of patterns (the motif space) being assessed for enrichment. QSLiMFinder uses specific query protein information to restrict the motif space and thereby increase the sensitivity and specificity of predictions.

Results: QSLiMFinder was extensively benchmarked using known SLiM-containing proteins and simulated protein interaction datasets of real human proteins. Exploiting prior knowledge of a query protein likely to be involved in a SLiM-mediated interaction increased the proportion of true positives correctly returned and reduced the proportion of datasets returning a false positive prediction. The biggest improvement was seen if a short region of the query protein flanking the interaction site was known.

Availability and implementation: All the tools and data used in this study, including QSLiMFinder and the SLiMBench benchmarking software, are freely available under a GNU license as part of SLiMSuite, at: http://bioware.soton.ac.uk.

Contact:richard.edwards@unsw.edu.au

Supplementary information:Supplementary data are available at Bioinformatics online.

LMPID: a manually curated database of linear motifs mediating protein-protein interactions

Linear motifs (LMs), used by a subset of all protein-protein interactions (PPIs), bind to globular receptors or domains and play an important role in signaling networks. LMPID (Linear Motif mediated Protein Interaction Database) is a manually curated database which provides comprehensive experimentally validated information about the LMs mediating PPIs from all organisms on a single platform. About 2200 entries have been compiled by detailed manual curation of PubMed abstracts, of which about 1000 LM entries were being annotated for the first time, as compared with the Eukaryotic LM resource. The users can submit their query through a user-friendly search page and browse the data in the alphabetical order of the bait gene names and according to the domains interacting with the LM. LMPID is freely accessible at http://bicresources.jcbose. ac.in/ssaha4/lmpid and contains 1750 unique LM instances found within 1181 baits interacting with 552 prey proteins. In summary, LMPID is an attempt to enrich the existing repertoire of resources available for studying the LMs implicated in PPIs and may help in understanding the patterns of LMs binding to a specific domain and develop prediction model to identify novel LMs specific to a domain and further able to predict inhibitors/modulators of PPI of interest.

Embryonic expression of endogenous retroviral RNAs in somatic tissues adjacent to the Oikopleura germline

Selective pressure to maintain small genome size implies control of transposable elements, and most old classes of retrotransposons are indeed absent from the very compact genome of the tunicate Oikopleura dioica. Nonetheless, two families of retrotransposons are present, including the Tor elements. The gene organization within Tor elements is similar to that of LTR retrotransposons and retroviruses. In addition to gag and pol, many Tor elements carry a third gene encoding viral envelope-like proteins (Env) that may mediate infection. We show that the Tor family contains distinct classes of elements. In some classes, env mRNA is transcribed from the 5′LTR as in retroviruses. In others, env is transcribed from an additional promoter located downstream of the 5′LTR. Tor Env proteins are membrane-associated glycoproteins which exhibit some features of viral membrane fusion proteins. Whereas some elements are expressed in the adult testis, many others are specifically expressed in embryonic somatic cells adjacent to primordial germ cells. Such embryonic expression depends on determinants present in the Tor elements and not on their surrounding genomic environment. Our study shows that unusual modes of transcription and expression close to the germline may contribute to the proliferation of Tor elements.

The plant Polycomb repressive complex 1 (PRC1) existed in the ancestor of seed plants and has a complex duplication history

Background

Polycomb repressive complex 1 (PRC1) is an essential protein complex for plant development. It catalyzes ubiquitination of histone H2A that is an important part of the transcription repression machinery. Absence of PRC1 subunits in Arabidopsis thaliana plants causes severe developmental defects. Many aspects of the plant PRC1 are elusive, including its origin and phylogenetic distribution.

Results

We established the evolutionary history of the plant PRC1 subunits (LHP1, Ring1a-b, Bmi1a-c, EMF1, and VRN1), enabled by sensitive phylogenetic methods and newly sequenced plant genomes from previously unsampled taxonomic groups.

We showed that all PRC1 core subunits exist in gymnosperms, earlier than previously thought, and that VRN1 is a recent addition, found exclusively in eudicots. The retention of individual subunits in chlorophytes, mosses, lycophytes and monilophytes indicates that they can moonlight as part of other complexes or processes. Moreover, we showed that most PRC1 subunits underwent a complex, duplication-rich history that differs significantly between Brassicaceae and other eudicots.

Conclusions

PRC1 existed in the last common ancestor of seed plants where it likely played an important regulatory role, aiding their radiation. The presence of LHP1, Ring1 and Bmi1 in mosses, lycophytes and monilophytes also suggests the presence of a primitive yet functional PRC1.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0319-z) contains supplementary material, which is available to authorized users.

14-3-3-Pred: improved methods to predict 14-3-3-binding phosphopeptides

MOTIVATION

The 14-3-3 family of phosphoprotein-binding proteins regulates many cellular processes by docking onto pairs of phosphorylated Ser and Thr residues in a constellation of intracellular targets. Therefore, there is a pressing need to develop new prediction methods that use an updated set of 14-3-3-binding motifs for the identification of new 14-3-3 targets and to prioritize the downstream analysis of >2000 potential interactors identified in high-throughput experiments.

RESULTS

Here, a comprehensive set of 14-3-3-binding targets from the literature was used to develop 14-3-3-binding phosphosite predictors. Position-specific scoring matrix, support vector machines (SVM) and artificial neural network (ANN) classification methods were trained to discriminate experimentally determined 14-3-3-binding motifs from non-binding phosphopeptides. ANN, position-specific scoring matrix and SVM methods showed best performance for a motif window spanning from -6 to +4 around the binding phosphosite, achieving Matthews correlation coefficient of up to 0.60. Blind prediction showed that all three methods outperform two popular 14-3-3-binding site predictors, Scansite and ELM. The new methods were used for prediction of 14-3-3-binding phosphosites in the human proteome. Experimental analysis of high-scoring predictions in the FAM122A and FAM122B proteins confirms the predictions and suggests the new 14-3-3-predictors will be generally useful.

AVAILABILITY AND IMPLEMENTATION

A standalone prediction web server is available at http://www.compbio.dundee.ac.uk/1433pred. Human candidate 14-3-3-binding phosphosites were integrated in ANIA: ANnotation and Integrated Analysis of the 14-3-3 interactome database.

Comparative analysis of the conserved functions of Arabidopsis DRL1 and yeast KTI12

Patterning of the polar axis during the early leaf developmental stage is established by cell-to-cell communication between the shoot apical meristem (SAM) and the leaf primordia. In a previous study, we showed that the DRL1 gene, which encodes a homolog of the Elongator-associated protein KTI12 of yeast, acts as a positive regulator of adaxial leaf patterning and shoot meristem activity. To determine the evolutionally conserved functions of DRL1, we performed a comparison of the deduced amino acid sequence of DRL1 and its yeast homolog, KTI12, and found that while overall homology was low, well-conserved domains were presented. DRL1 contained two conserved plant-specific domains. Expression of the DRL1 gene in a yeast KTI12-deficient yeast mutant suppressed the growth retardation phenotype, but did not rescue the caffeine sensitivity, indicating that the role of Arabidopsis Elongator-associated protein is partially conserved with yeast KTI12, but may have changed between yeast and plants in response to caffeine during the course of evolution. In addition, elevated expression of DRL1 gene triggered zymocin sensitivity, while overexpression of KTI12 maintained zymocin resistance, indicating that the function of Arabidopsis DRL1 may not overlap with yeast KTI12 with regards to toxin sensitivity. In this study, expression analysis showed that class-I KNOX genes were downregulated in the shoot apex, and that YAB and KAN were upregulated in leaves of the Arabidopsis drl1-101 mutant. Our results provide insight into the communication network between the SAM and leaf primordia required for the establishment of leaf polarity by mediating histone acetylation or through other mechanisms.

A highly conserved program of neuronal microexons is misregulated in autistic brains

Alternative splicing (AS) generates vast transcriptomic and proteomic complexity. However, which of the myriad of detected AS events provide important biological functions is not well understood. Here, we define the largest program of functionally coordinated, neural-regulated AS described to date in mammals. Relative to all other types of AS within this program, 3-15 nucleotide "microexons" display the most striking evolutionary conservation and switch-like regulation. These microexons modulate the function of interaction domains of proteins involved in neurogenesis. Most neural microexons are regulated by the neuronal-specific splicing factor nSR100/SRRM4, through its binding to adjacent intronic enhancer motifs. Neural microexons are frequently misregulated in the brains of individuals with autism spectrum disorder, and this misregulation is associated with reduced levels of nSR100. The results thus reveal a highly conserved program of dynamic microexon regulation associated with the remodeling of protein-interaction networks during neurogenesis, the misregulation of which is linked to autism.

Versatile communication strategies among tandem WW domain repeats

Interactions mediated by short linear motifs in proteins play major roles in regulation of cellular homeostasis since their transient nature allows for easy modulation. We are still far from a full understanding and appreciation of the complex regulation patterns that can be, and are, achieved by this type of interaction. The fact that many linear-motif-binding domains occur in tandem repeats in proteins indicates that their mutual communication is used extensively to obtain complex integration of information toward regulatory decisions. This review is an attempt to overview, and classify, different ways by which two and more tandem repeats cooperate in binding to their targets, in the well-characterized family of WW domains and their corresponding polyproline ligands.

Binding properties of SUMO-interacting motifs (SIMs) in yeast

Small ubiquitin-like modifier (SUMO) conjugation and interaction play an essential role in many cellular processes. A large number of yeast proteins is known to interact non-covalently with SUMO via short SUMO-interacting motifs (SIMs), but the structural details of this interaction are yet poorly characterized. In the present work, sequence analysis of a large dataset of 148 yeast SIMs revealed the existence of a hydrophobic core binding motif and a preference for acidic residues either within or adjacent to the core motif. Thus the sequence properties of yeast SIMs are highly similar to those described for human. Molecular dynamics simulations were performed to investigate the binding preferences for four representative SIM peptides differing in the number and distribution of acidic residues. Furthermore, the relative stability of two previously observed alternative binding orientations (parallel, antiparallel) was assessed. For all SIMs investigated, the antiparallel binding mode remained stable in the simulations and the SIMs were tightly bound via their hydrophobic core residues supplemented by polar interactions of the acidic residues. In contrary, the stability of the parallel binding mode is more dependent on the sequence features of the SIM motif like the number and position of acidic residues or the presence of additional adjacent interaction motifs. This information should be helpful to enhance the prediction of SIMs and their binding properties in different organisms to facilitate the reconstruction of the SUMO interactome.

Human FAM154A (SAXO1) is a microtubule-stabilizing protein specific to cilia and related structures

Cilia and flagella are microtubule-based organelles present at the surface of most cells, ranging from protozoa to vertebrates, in which these structures are implicated in processes from morphogenesis to cell motility. In vertebrate neurons, microtubule-associated MAP6 proteins stabilize cold-resistant microtubules through their Mn and Mc modules, and play a role in synaptic plasticity. Although centrioles, cilia and flagella have cold-stable microtubules, MAP6 proteins have not been identified in these organelles, suggesting that additional proteins support this role in these structures. Here, we characterize human FAM154A (hereafter referred to as hSAXO1) as the first human member of a widely conserved family of MAP6-related proteins specific to centrioles and cilium microtubules. Our data demonstrate that hSAXO1 binds specifically to centriole and cilium microtubules. We identify, in vivo and in vitro, hSAXO1 Mn modules as responsible for microtubule binding and stabilization as well as being necessary for ciliary localization. Finally, overexpression and knockdown studies show that hSAXO1 modulates axoneme length. Taken together, our findings suggest a fine regulation of hSAXO1 localization and important roles in cilium biogenesis and function.

The ABBA motif binds APC/C activators and is shared by APC/C substrates and regulators

The anaphase-promoting complex or cyclosome (APC/C) is the ubiquitin ligase that regulates mitosis by targeting specific proteins for degradation at specific times under the control of the spindle assembly checkpoint (SAC). How the APC/C recognizes its different substrates is a key problem in the control of cell division. Here, we have identified the ABBA motif in cyclin A, BUBR1, BUB1, and Acm1, and we show that it binds to the APC/C coactivator CDC20. The ABBA motif in cyclin A is required for its proper degradation in prometaphase through competing with BUBR1 for the same site on CDC20. Moreover, the ABBA motifs in BUBR1 and BUB1 are necessary for the SAC to work at full strength and to recruit CDC20 to kinetochores. Thus, we have identified a conserved motif integral to the proper control of mitosis that connects APC/C substrate recognition with the SAC.

Insect inhibitor-of-apoptosis (IAP) proteins are negatively regulated by signal-induced N-terminal degrons absent within viral IAP proteins

Inhibitor-of-apoptosis (IAP) proteins are key regulators of the innate antiviral response by virtue of their capacity to respond to signals affecting cell survival. In insects, wherein the host IAP provides a primary restriction to apoptosis, diverse viruses trigger rapid IAP depletion that initiates caspase-mediated apoptosis, thereby limiting virus multiplication. We report here that the N-terminal leader of two insect IAPs, Spodoptera frugiperda SfIAP and Drosophila melanogaster DIAP1, contain distinct instability motifs that regulate IAP turnover and apoptotic consequences. Functioning as a protein degron, the cellular IAP leader dramatically shortened the life span of a long-lived viral IAP (Op-IAP3) when fused to its N terminus. The SfIAP degron contains mitogen-activated kinase (MAPK)-like regulatory sites, responsible for MAPK inhibitor-sensitive phosphorylation of SfIAP. Hyperphosphorylation correlated with increased SfIAP turnover independent of the E3 ubiquitin-ligase activity of the SfIAP RING, which also regulated IAP stability. Together, our findings suggest that the SfIAP phospho-degron responds rapidly to a signal-activated kinase cascade, which regulates SfIAP levels and thus apoptosis. The N-terminal leader of dipteran DIAP1 also conferred virus-induced IAP depletion by a caspase-independent mechanism. DIAP1 instability mapped to previously unrecognized motifs that are not found in lepidopteran IAPs. Thus, the leaders of cellular IAPs from diverse insects carry unique signal-responsive degrons that control IAP turnover. Rapid response pathways that trigger IAP degradation and initiate apoptosis independent of canonical prodeath gene (Reaper-Grim-Hid) expression may provide important innate immune advantages. Furthermore, the elimination of these response motifs within viral IAPs, including those of baculoviruses, explains their unusual stability and their potent antiapoptotic activity.

IMPORTANCE

Apoptosis is an effective means by which a host controls virus infection. In insects, inhibitor-of-apoptosis (IAP) proteins act as regulatory sentinels by responding to cellular signals that determine the fate of infected cells. We discovered that lepidopteran (moth and butterfly) IAPs, which are degraded upon baculovirus infection, are controlled by a conserved phosphorylation-sensitive degron within the IAP N-terminal leader. The degron likely responds to virus-induced kinase-specific signals for degradation through SKP1/Cullin/F-box complex-mediated ubiquitination. Such signal-induced destruction of cellular IAPs is distinct from degradation caused by well-known IAP antagonists, which act to expel IAP-bound caspases. The major implication of this study is that insects have multiple signal-responsive mechanisms by which the sentinel IAPs are actively degraded to initiate host apoptosis. Such diversity of pathways likely provides insects with rapid and efficient strategies for pathogen control. Furthermore, the absence of analogous degrons in virus-encoded IAPs explains their relative stability and antiapoptotic potency.

Gamma-sarcoglycan is required for the response of archvillin to mechanical stimulation in skeletal muscle

Loss of gamma-sarcoglycan (γ-SG) induces muscle degeneration and signaling defects in response to mechanical load, and its absence is common to both Duchenne and limb girdle muscular dystrophies. Growing evidence suggests that aberrant signaling contributes to the disease pathology; however, the mechanisms of γ-SG-mediated mechanical signaling are poorly understood. To uncover γ-SG signaling pathway components, we performed yeast two-hybrid screens and identified the muscle-specific protein archvillin as a γ-SG and dystrophin interacting protein. Archvillin protein and message levels were significantly upregulated at the sarcolemma of murine γ-SG-null (gsg(-/-)) muscle but delocalized in dystrophin-deficient mdx muscle. Similar elevation of archvillin protein was observed in human quadriceps muscle lacking γ-SG. Reintroduction of γ-SG in gsg(-/-) muscle by rAAV injection restored archvillin levels to that of control C57 muscle. In situ eccentric contraction of tibialis anterior (TA) muscles from C57 mice caused ERK1/2 phosphorylation, nuclear activation of P-ERK1/2 and stimulus-dependent archvillin association with P-ERK1/2. In contrast, TA muscles from gsg(-/-) and mdx mice exhibited heightened P-ERK1/2 and increased nuclear P-ERK1/2 localization following eccentric contractions, but the archvillin-P-ERK1/2 association was completely ablated. These results position archvillin as a mechanically sensitive component of the dystrophin complex and demonstrate that signaling defects caused by loss of γ-SG occur both at the sarcolemma and in the nucleus.