James B. Procter

Jalview Version 2--a multiple sequence alignment editor and analysis workbench.

Summary: Jalview Version 2 is a system for interactive WYSIWYG editing, analysis and annotation of multiple sequence alignments. Core features include keyboard and mouse-based editing, multiple views and alignment overviews, and linked structure display with Jmol. Jalview 2 is available in two forms: a lightweight Java applet for use in web applications, and a powerful desktop application that employs web services for sequence alignment, secondary structure prediction and the retrieval of alignments, sequences, annotation and structures from public databases and any DAS 1.53 compliant sequence or annotation server. Availability: The Jalview 2 Desktop application and JalviewLite applet are made freely available under the GPL, and can be downloaded from www.jalview.org Contact: g.j.barton@dundee.ac.uk

JPred4: a protein secondary structure prediction server

JPred4 (http://www.compbio.dundee.ac.uk/jpred4) is the latest version of the popular JPred protein secondary structure prediction server which provides predictions by the JNet algorithm, one of the most accurate methods for secondary structure prediction. In addition to protein secondary structure, JPred also makes predictions of solvent accessibility and coiled-coil regions. The JPred service runs up to 94 000 jobs per month and has carried out over 1.5 million predictions in total for users in 179 countries. The JPred4 web server has been re-implemented in the Bootstrap framework and JavaScript to improve its design, usability and accessibility from mobile devices. JPred4 features higher accuracy, with a blind three-state (α-helix, β-strand and coil) secondary structure prediction accuracy of 82.0% while solvent accessibility prediction accuracy has been raised to 90% for residues <5% accessible. Reporting of results is enhanced both on the website and through the optional email summaries and batch submission results. Predictions are now presented in SVG format with options to view full multiple sequence alignments with and without gaps and insertions. Finally, the help-pages have been updated and tool-tips added as well as step-by-step tutorials.

MSAViewer: interactive JavaScript visualization of multiple sequence alignments

Summary: The MSAViewer is a quick and easy visualization and analysis JavaScript component for Multiple Sequence Alignment data of any size. Core features include interactive navigation through the alignment, application of popular color schemes, sorting, selecting and filtering. The MSAViewer is ‘web ready’: written entirely in JavaScript, compatible with modern web browsers and does not require any specialized software. The MSAViewer is part of the BioJS collection of components. Availability and Implementation: The MSAViewer is released as open source software under the Boost Software License 1.0. Documentation, source code and the viewer are available at http://msa.biojs.net/. Supplementary information: Supplementary data are available at Bioinformatics online. Contact: msa@bio.sh

Distinct donor and acceptor specificities of Trypanosoma brucei oligosaccharyltransferases

Asparagine‐linked glycosylation is catalysed by oligosaccharyltransferase (OTase). In Trypanosoma brucei OTase activity is catalysed by single‐subunit enzymes encoded by three paralogous genes of which TbSTT3B and TbSTT3C can complement a yeast Δstt3 mutant. The two enzymes have overlapping but distinct peptide acceptor specificities, with TbSTT3C displaying an enhanced ability to glycosylate sites flanked by acidic residues. TbSTT3A and TbSTT3B, but not TbSTT3C, are transcribed in the bloodstream and procyclic life cycle stages of T. brucei. Selective knockdown and analysis of parasite protein N‐glycosylation showed that TbSTT3A selectively transfers biantennary Man5GlcNAc2 to specific glycosylation sites whereas TbSTT3B selectively transfers triantennary Man9GlcNAc2 to others. Analysis of T. brucei glycosylation site occupancy showed that TbSTT3A and TbSTT3B glycosylate sites in acidic to neutral and neutral to basic regions of polypeptide, respectively. This embodiment of distinct specificities in single‐subunit OTases may have implications for recombinant glycoprotein engineering. TbSTT3A and TbSTT3B could be knocked down individually, but not collectively, in tissue culture. However, both were independently essential for parasite growth in mice, suggesting that inhibiting protein N‐glycosylation could have therapeutic potential against trypanosomiasis.

Visualizing biological data-now and in the future.

Methods and tools for visualizing biological data have improved considerably over the last decades, but they are still inadequate for some high-throughput data sets. For most users, a key challenge is to benefit from the deluge of data without being overwhelmed by it. This challenge is still largely unfulfilled and will require the development of truly integrated and highly useable tools.

Phosphoproteomic screening identifies Rab GTPases as novel downstream targets of PINK1

Mutations in the PTEN‐induced kinase 1 (PINK1) are causative of autosomal recessive Parkinsons disease (PD). We have previously reported that PINK1 is activated by mitochondrial depolarisation and phosphorylates serine 65 (Ser65) of the ubiquitin ligase Parkin and ubiquitin to stimulate Parkin E3 ligase activity. Here, we have employed quantitative phosphoproteomics to search for novel PINK1‐dependent phosphorylation targets in HEK (human embryonic kidney) 293 cells stimulated by mitochondrial depolarisation. This led to the identification of 14,213 phosphosites from 4,499 gene products. Whilst most phosphosites were unaffected, we strikingly observed three members of a sub‐family of Rab GTPases namely Rab8A, 8B and 13 that are all phosphorylated at the highly conserved residue of serine 111 (Ser111) in response to PINK1 activation. Using phospho‐specific antibodies raised against Ser111 of each of the Rabs, we demonstrate that Rab Ser111 phosphorylation occurs specifically in response to PINK1 activation and is abolished in HeLa PINK1 knockout cells and mutant PINK1 PD patient‐derived fibroblasts stimulated by mitochondrial depolarisation. We provide evidence that Rab8A GTPase Ser111 phosphorylation is not directly regulated by PINK1 in vitro and demonstrate in cells the time course of Ser111 phosphorylation of Rab8A, 8B and 13 is markedly delayed compared to phosphorylation of Parkin at Ser65. We further show mechanistically that phosphorylation at Ser111 significantly impairs Rab8A activation by its cognate guanine nucleotide exchange factor (GEF), Rabin8 (by using the Ser111Glu phosphorylation mimic). These findings provide the first evidence that PINK1 is able to regulate the phosphorylation of Rab GTPases and indicate that monitoring phosphorylation of Rab8A/8B/13 at Ser111 may represent novel biomarkers of PINK1 activity in vivo. Our findings also suggest that disruption of Rab GTPase‐mediated signalling may represent a major mechanism in the neurodegenerative cascade of Parkinsons disease.

Visualization of multiple alignments, phylogenies and gene family evolution

Software for visualizing sequence alignments and trees are essential tools for life scientists. In this review, we describe the major features and capabilities of a selection of stand-alone and web-based applications useful when investigating the function and evolution of a gene family. These range from simple viewers, to systems that provide sophisticated editing and analysis functions. We conclude with a discussion of the challenges that these tools now face due to the flood of next generation sequence data and the increasingly complex network of bioinformatics information sources.

Functional analysis of the methylmalonyl-CoA epimerase from Caenorhabditis elegans.

Methylmalonyl‐CoA epimerase (MCE) is an enzyme involved in the propionyl‐CoA metabolism that is responsible for the degradation of branched amino acids and odd‐chain fatty acids. This pathway typically functions in the reversible conversion of propionyl‐CoA to succinyl‐CoA. The Caenorhabditis elegans genome contains a single gene encoding MCE (mce‐1) corresponding to a 15 kDa protein. This was expressed in Escherichia coli and the enzymatic activity was determined. Analysis of the protein expression pattern at both the tissue and subcellular level by microinjection of green fluorescent protein constructs revealed expression in the pharynx, hypodermis and, most prominently in body wall muscles. The subcellular pattern agrees with predictions of mitochondrial localization. The sequence similarity to an MCE of known structure was high enough to permit a three‐dimensional model to be built, suggesting conservation of ligand and metal binding sites. Comparison with corresponding sequences from a variety of organisms shows more than 1/6 of the sequence is completely conserved. Mutants allelic to mce‐1 showed no obvious phenotypic alterations, demonstrating that the enzyme is not essential for normal worm development under laboratory conditions. However, survival of the knockout mutants was altered when exposed to stress conditions, with mutants surprisingly showing an increased resistance to oxidative stress.

MACSIMS : multiple alignment of complete sequences information management system

BackgroundIn the post-genomic era, systems-level studies are being performed that seek to explain complex biological systems by integrating diverse resources from fields such as genomics, proteomics or transcriptomics. New information management systems are now needed for the collection, validation and analysis of the vast amount of heterogeneous data available. Multiple alignments of complete sequences provide an ideal environment for the integration of this information in the context of the protein family.ResultsMACSIMS is a multiple alignment-based information management program that combines the advantages of both knowledge-based and ab initio sequence analysis methods. Structural and functional information is retrieved automatically from the public databases. In the multiple alignment, homologous regions are identified and the retrieved data is evaluated and propagated from known to unknown sequences with these reliable regions. In a large-scale evaluation, the specificity of the propagated sequence features is estimated to be >99%, i.e. very few false positive predictions are made. MACSIMS is then used to characterise mutations in a test set of 100 proteins that are known to be involved in human genetic diseases. The number of sequence features associated with these proteins was increased by 60%, compared to the features available in the public databases. An XML format output file allows automatic parsing of the MACSIM results, while a graphical display using the JalView program allows manual analysis.ConclusionMACSIMS is a new information management system that incorporates detailed analyses of protein families at the structural, functional and evolutionary levels. MACSIMS thus provides a unique environment that facilitates knowledge extraction and the presentation of the most pertinent information to the biologist. A web server and the source code are available at http://bips.u-strasbg.fr/MACSIMS/.

Wurst: a protein threading server with a structural scoring function, sequence profiles and optimized substitution matrices

Wurst is a protein threading program with an emphasis on high quality sequence to structure alignments (http://www.zbh.uni-hamburg.de/wurst). Submitted sequences are aligned to each of about 3000 templates with a conventional dynamic programming algorithm, but using a score function with sophisticated structure and sequence terms. The structure terms are a log-odds probability of sequence to structure fragment compatibility, obtained from a Bayesian classification procedure. A simplex optimization was used to optimize the sequence-based terms for the goal of alignment and model quality and to balance the sequence and structural contributions against each other. Both sequence and structural terms operate with sequence profiles.