Combining High-Resolution and Exact Calibration To Boost Statistical Power: A Well-Calibrated Score Function for High-Resolution MS2 Data.

To achieve accurate assignment of peptide sequences to observed fragmentation spectra, a shotgun proteomics database search tool must make good use of the very high-resolution information produced by state-of-the-art mass spectrometers. However, making use of this information while also ensuring that the search engine's scores are well calibrated, that is, that the score assigned to one spectrum can be meaningfully compared to the score assigned to a different spectrum, has proven to be challenging. Here we describe a database search score function, the "residue evidence" (res-ev) score, that achieves both of these goals simultaneously. We also demonstrate how to combine calibrated res-ev scores with calibrated XCorr scores to produce a "combined p value" score function. We provide a benchmark consisting of four mass spectrometry data sets, which we use to compare the combined p value to the score functions used by several existing search engines. Our results suggest that the combined p value achieves state-of-the-art performance, generally outperforming MS Amanda and Morpheus and performing comparably to MS-GF+. The res-ev and combined p-value score functions are freely available as part of the Tide search engine in the Crux mass spectrometry toolkit ( http://crux.ms ).

Computing exact p-values for a cross-correlation shotgun proteomics score function.

The core of every protein mass spectrometry analysis pipeline is a function that assesses the quality of a match between an observed spectrum and a candidate peptide. We describe a procedure for computing exact p-values for the oldest and still widely used score function, SEQUEST XCorr. The procedure uses dynamic programming to enumerate efficiently the full distribution of scores for all possible peptides whose masses are close to that of the spectrum precursor mass. Ranking identified spectra by p-value rather than XCorr significantly reduces variance because of spectrum-specific effects on the score. In combination with the Percolator postprocessor, the XCorr p-value yields more spectrum and peptide identifications at a fixed false discovery rate than Mascot, X!Tandem, Comet, and MS-GF+ across a variety of data sets.

Crux: rapid open source protein tandem mass spectrometry analysis.

Efficiently and accurately analyzing big protein tandem mass spectrometry data sets requires robust software that incorporates state-of-the-art computational, machine learning, and statistical methods. The Crux mass spectrometry analysis software toolkit ( http://cruxtoolkit.sourceforge.net ) is an open source project that aims to provide users with a cross-platform suite of analysis tools for interpreting protein mass spectrometry data.

MR-Tandem: parallel X!Tandem using Hadoop MapReduce on Amazon Web Services.

SUMMARY: MR-Tandem adapts the popular X!Tandem peptide search engine to work with Hadoop MapReduce for reliable parallel execution of large searches. MR-Tandem runs on any Hadoop cluster but offers special support for Amazon Web Services for creating inexpensive on-demand Hadoop clusters, enabling search volumes that might not otherwise be feasible with the compute resources a researcher has at hand. MR-Tandem is designed to drop in wherever X!Tandem is already in use and requires no modification to existing X!Tandem parameter files, and only minimal modification to X!Tandem-based workflows. AVAILABILITY AND IMPLEMENTATION: MR-Tandem is implemented as a lightly modified X!Tandem C++ executable and a Python script that drives Hadoop clusters including Amazon Web Services (AWS) Elastic Map Reduce (EMR), using the modified X!Tandem program as a Hadoop Streaming mapper and reducer. The modified X!Tandem C++ source code is Artistic licensed, supports pluggable scoring, and is available as part of the Sashimi project at http://sashimi.svn.sourceforge.net/viewvc/sashimi/trunk/trans_proteomic_pipeline/extern/xtandem/. The MR-Tandem Python script is Apache licensed and available as part of the Insilicos Cloud Army project at http://ica.svn.sourceforge.net/viewvc/ica/trunk/mr-tandem/. Full documentation and a windows installer that configures MR-Tandem, Python and all necessary packages are available at this same URL. CONTACT: brian.pratt@insilicos.com

PREDICTD PaRallel Epigenomics Data Imputation with Cloud-based Tensor Decomposition.

The Encyclopedia of DNA Elements (ENCODE) and the Roadmap Epigenomics Project seek to characterize the epigenome in diverse cell types using assays that identify, for example, genomic regions with modified histones or accessible chromatin. These efforts have produced thousands of datasets but cannot possibly measure each epigenomic factor in all cell types. To address this, we present a method, PaRallel Epigenomics Data Imputation with Cloud-based Tensor Decomposition (PREDICTD), to computationally impute missing experiments. PREDICTD leverages an elegant model called "tensor decomposition" to impute many experiments simultaneously. Compared with the current state-of-the-art method, ChromImpute, PREDICTD produces lower overall mean squared error, and combining the two methods yields further improvement. We show that PREDICTD data captures enhancer activity at noncoding human accelerated regions. PREDICTD provides reference imputed data and open-source software for investigating new cell types, and demonstrates the utility of tensor decomposition and cloud computing, both promising technologies for bioinformatics.

Correction: Forecasting Seizures Using Bivariate Intracranial EEG Measures and SVM in Naturally Occurring Canine Epilepsy.

[This corrects the article DOI: 10.1371/journal.pone.0133900.].

Forecasting Seizures Using Intracranial EEG Measures and SVM in Naturally Occurring Canine Epilepsy.

Management of drug resistant focal epilepsy would be greatly assisted by a reliable warning system capable of alerting patients prior to seizures to allow the patient to adjust activities or medication. Such a system requires successful identification of a preictal, or seizure-prone state. Identification of preictal states in continuous long- duration intracranial electroencephalographic (iEEG) recordings of dogs with naturally occurring epilepsy was investigated using a support vector machine (SVM) algorithm. The dogs studied were implanted with a 16-channel ambulatory iEEG recording device with average channel reference for a mean (st. dev.) of 380.4 (+87.5) days producing 220.2 (+104.1) days of intracranial EEG recorded at 400 Hz for analysis. The iEEG records had 51.6 (+52.8) seizures identified, of which 35.8 (+30.4) seizures were preceded by more than 4 hours of seizure-free data. Recorded iEEG data were stratified into 11 contiguous, non-overlapping frequency bands and binned into one-minute synchrony features for analysis. Performance of the SVM classifier was assessed using a 5-fold cross validation approach, where preictal training data were taken from 90 minute windows with a 5 minute pre-seizure offset. Analysis of the optimal preictal training time was performed by repeating the cross validation over a range of preictal windows and comparing results. We show that the optimization of feature selection varies for each subject, i.e. algorithms are subject specific, but achieve prediction performance significantly better than a time-matched Poisson random predictor (p<0.05) in 5/5 dogs analyzed.

Forecasting seizures in dogs with naturally occurring epilepsy.

Seizure forecasting has the potential to create new therapeutic strategies for epilepsy, such as providing patient warnings and delivering preemptive therapy. Progress on seizure forecasting, however, has been hindered by lack of sufficient data to rigorously evaluate the hypothesis that seizures are preceded by physiological changes, and are not simply random events. We investigated seizure forecasting in three dogs with naturally occurring focal epilepsy implanted with a device recording continuous intracranial EEG (iEEG). The iEEG spectral power in six frequency bands: delta (0.1-4 Hz), theta (4-8 Hz), alpha (8-12 Hz), beta (12-30 Hz), low-gamma (30-70 Hz), and high-gamma (70-180 Hz), were used as features. Logistic regression classifiers were trained to discriminate labeled pre-ictal and inter-ictal data segments using combinations of the band spectral power features. Performance was assessed on separate test data sets via 10-fold cross-validation. A total of 125 spontaneous seizures were detected in continuous iEEG recordings spanning 6.5 to 15 months from 3 dogs. When considering all seizures, the seizure forecasting algorithm performed significantly better than a Poisson-model chance predictor constrained to have the same time in warning for all 3 dogs over a range of total warning times. Seizure clusters were observed in all 3 dogs, and when the effect of seizure clusters was decreased by considering the subset of seizures separated by at least 4 hours, the forecasting performance remained better than chance for a subset of algorithm parameters. These results demonstrate that seizures in canine epilepsy are not randomly occurring events, and highlight the feasibility of long-term seizure forecasting using iEEG monitoring.

Identification of TNF-alpha inhibitors from a split-pool library based on a tyrosine-proline peptidomimetic scaffold.

The design and synthesis of a combinatorial library based on a 4-aryloxyproline scaffold with tyrosine as the aryl portion is described. The 1728 member library was prepared using the split-pool method to generate pools of compounds. Screening of the library components as mixtures followed by deconvolution led to the discovery of novel inhibitors of TNF-alpha induced apoptosis.

Benzobicyclooctanes as novel inhibitors of TNF-alpha signaling.

A novel series of TNF-alpha inhibitors based on a benzobicyclooctane scaffold is reported. The compounds display good potency in inhibiting TNF-alpha induced apoptosis and NF kappa B activation. Additionally, they are selective for TNF-alpha as they do not inhibit apoptosis induced by soluble Fas ligand. The compounds described here can act as leads for future medicinal chemistry efforts and may also be useful tools for elucidating the TNF-alpha signaling pathway.