Linear discriminant analysis-based estimation of the false discovery rate for phosphopeptide identifications

The development of liquid chromatography coupled with tandem mass spectrometry (LC-MS/MS) has made it possible to characterize phosphopeptides in an increasingly large-scale and high-throughput fashion. However, extracting confident phosphopeptide identifications from the resulting large data sets in a similar high-throughput fashion remains difficult, as does rigorously estimating the false discovery rate (FDR) of a set of phosphopeptide identifications. This article describes a data analysis pipeline designed to address these issues. The first step is to reanalyze phosphopeptide identifications that contain ambiguous assignments for the incorporated phosphate(s) to determine the most likely arrangement of the phosphate(s). The next step is to employ an expectation maximization algorithm to estimate the joint distribution of the peptide scores. A linear discriminant analysis is then performed to determine how to optimally combine peptide scores (in this case, from SEQUEST) into a discriminant score that possesses the maximum discriminating power. Based on this discriminant score, the p- and q-values for each phosphopeptide identification are calculated, and the phosphopeptide identification FDR is then estimated. This data analysis approach was applied to data from a study of irradiated human skin fibroblasts to provide a robust estimate of FDR for phosphopeptides. The Phosphopeptide FDR Estimator software is freely available for download at http://ncrr.pnl.gov/software/.     

Comparative proteomics of human monkeypox and vaccinia intracellular mature and extracellular enveloped virions.

Orthopoxviruses are among the largest and most complex of the animal viruses. In response to the recent emergence of monkeypox in Africa and the threat of smallpox bioterrorism, two orthopoxviruses with different pathogenic potentials, human monkeypox virus and vaccinia virus, were proteomically compared with the goal of identifying proteins required for pathogenesis. Orthopoxviruses were grown in HeLa cells to two different viral forms (intracellular mature virus and extracellular enveloped virus), purified by sucrose gradient ultracentrifugation, denatured using RapiGest surfactant, and digested with trypsin. Unfractionated samples and strong cation exchange HPLC fractions were analyzed by high-resolution reversed-phase nano-LC-MS/MS, and analyses of the MS/MS spectra using SEQUEST and X! Tandem resulted in the confident identification of hundreds of monkeypox, vaccinia, and copurified host-cell proteins. The unfractionated samples were additionally analyzed by LC-MS using an LTQ-Orbitrap, and the accurate mass and elution time tag approach was used to perform quantitative comparisons. Possible pathophysiological roles of differentially abundant Orthopoxvirus proteins are discussed. Data, processed results, and protocols are available at http://www.proteomicsresource.org/.     

Modelling the transition from simple to complex Ca2+oscillations in pancreatic acinar cells

A mathematical model is proposed which systematically investigates complex calcium oscillations in pancreatic acinar cells. This model is based on calcium-induced calcium release via inositol trisphosphate receptors (IPR) and ryanodine receptors (RyR) and includes calcium modulation of inositol (1,4,5) trisphosphate (IP₃) levels through feedback regulation of degradation and production. In our model, the apical and the basal regions are separated by a region containing mitochondria, which is capable of restricting Ca²⁺ responses to the apical region. We were able to reproduce the observed oscillatory patterns, from baseline spikes to sinusoidal oscillations. The model predicts that calcium-dependent production and degradation of IP₃ is a key mechanism for complex calcium oscillations in pancreatic acinar cells. A partial bifurcation analysis is performed which explores the dynamic behaviour of the model in both apical and basal regions.