Search

In cluster analysis, a common first step is to scale the data aiming to better partition them into clusters. Even though many different techniques have throughout many years been introduced to this end, it is probably fair to say that the workhorse in this preprocessing phase has been to divide the data by the standard deviation along each dimension. Like division by the standard deviation, the great majority of scaling techniques can be said to have roots in some sort of statistical take on the data. Here we explore the use of multidimensional shapes of data, aiming to obtain scaling factors for use prior to clustering by some method, like k-means, that makes explicit use of distances between samples. We borrow from the field of cosmology and related areas the recently introduced notion of shape complexity, which in the variant we use is a relatively simple, data-dependent nonlinear function that we show can be used to help with the determination of appropriate scaling factors. Focusing on what might be called "midrange" distances, we formulate a constrained nonlinear programming problem and use it to produce candidate scaling-factor sets that can be sifted on the basis of further considerations of the data, say via expert knowledge. We give results on some iconic data sets, highlighting the strengths and potential weaknesses of the new approach. These results are generally positive across all the data sets used., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2023 Aguilar, Barbosa. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Motivation: There are several well-established paradigms for identifying and pinpointing discriminative peptides/proteins using shotgun proteomic data; examples are peptide-spectrum matching, de novo sequencing, open searches, and even hybrid approaches. Such an arsenal of complementary paradigms can provide deep data coverage, albeit some unidentified discriminative peptides remain., Results: We present DiagnoMass, software tool that groups similar spectra into spectral clusters and then shortlists those clusters that are discriminative for biological conditions. DiagnoMass then communicates with proteomic tools to attempt the identification of such clusters. We demonstrate the effectiveness of DiagnoMass by analyzing proteomic data from Escherichia coli, Salmonella, and Shigella, listing many high-quality discriminative spectral clusters that had thus far remained unidentified by widely adopted proteomic tools. DiagnoMass can also classify proteomic profiles. We anticipate the use of DiagnoMass as a vital tool for pinpointing biomarkers., Availability: DiagnoMass and related documentation, including a usage protocol, are available at http://www.diagnomass.com., Competing Interests: Declaration of Competing Interest None declared., (Copyright © 2023 Elsevier B.V. All rights reserved.)

Complex protein mixtures typically generate many tandem mass spectra produced by different peptides coisolated in the gas phase. Widely adopted proteomic data analysis environments usually fail to identify most of these spectra, succeeding at best in identifying only one of the multiple cofragmenting peptides. We present PatternLab V (PLV), an updated version of PatternLab that integrates the YADA 3 deconvolution algorithm to handle such cases efficiently. In general, we expect an increase of 10% in spectral identifications when dealing with complex proteomic samples. PLV is freely available at http://patternlabforproteomics.org.

Motivation: Confident deconvolution of proteomic spectra is critical for several applications such as de novo sequencing, cross-linking mass spectrometry and handling chimeric mass spectra., Results: In general, all deconvolution algorithms may eventually report mass peaks that are not compatible with the chemical formula of any peptide. We show how to remove these artifacts by considering their mass defects. We introduce Y.A.D.A. 3.0, a fast deconvolution algorithm that can remove peaks with unacceptable mass defects. Our approach is effective for polypeptides with less than 10 kDa, and its essence can be easily incorporated into any deconvolution algorithm., Availability and Implementation: Y.A.D.A. 3.0 is freely available for academic use at http://patternlabforproteomics.org/yada3., Supplementary Information: Supplementary information is available at Bioinformatics online., (© The Author(s) 2022. Published by Oxford University Press. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Shotgun proteomics aims to identify and quantify the thousands of proteins in complex mixtures such as cell and tissue lysates and biological fluids. This approach uses liquid chromatography coupled with tandem mass spectrometry and typically generates hundreds of thousands of mass spectra that require specialized computational environments for data analysis. PatternLab for proteomics is a unified computational environment for analyzing shotgun proteomic data. PatternLab V (PLV) is the most comprehensive and crucial update so far, the result of intensive interaction with the proteomics community over several years. All PLV modules have been optimized and its graphical user interface has been completely updated for improved user experience. Major improvements were made to all aspects of the software, ranging from boosting the number of protein identifications to faster extraction of ion chromatograms. PLV provides modules for preparing sequence databases, protein identification, statistical filtering and in-depth result browsing for both labeled and label-free quantitation. The PepExplorer module can even pinpoint de novo sequenced peptides not already present in the database. PLV is of broad applicability and therefore suitable for challenging experimental setups, such as time-course experiments and data handling from unsequenced organisms. PLV interfaces with widely adopted software and community initiatives, e.g., Comet, Skyline, PEAKS and PRIDE. It is freely available at http://www.patternlabforproteomics.org ., (© 2022. Springer Nature Limited.)

In proteomics, the identification of peptides from mass spectral data can be mathematically described as the partitioning of mass spectra into clusters (i.e., groups of spectra derived from the same peptide). The way partitions are validated is just as important, having evolved side by side with the clustering algorithms themselves and given rise to many partition assessment measures. An assessment measure is said to have a selection bias if, and only if, the probability that a randomly chosen partition scoring a high value depends on the number of clusters in the partition. In the context of clustering mass spectra, this might mislead the validation process to favor clustering algorithms that generate too many (or few) spectral clusters, regardless of the underlying peptide sequence. A selection bias toward the number of peptides is desirable for proteomics as it estimates the number of peptides in a complex protein mixture. Here, we introduce an assessment measure that is purposely biased toward the number of peptide ion species. We also introduce a partition assessment framework for proteomics, called the Partition Assessment Tool, and demonstrate its importance by evaluating the performance of eight clustering algorithms on seven proteomics datasets while discussing the trade-offs involved. SIGNIFICANCE: Clustering algorithms are widely adopted in proteomics for undertaking several tasks such as speeding up search engines, generating consensus mass spectra, and to aid in the classification of proteomic profiles. Choosing which algorithm is most fit for the task at hand is not simple as each algorithm has advantages and disadvantages; furthermore, specifying clustering parameters is also a necessary and fundamental step. For example, deciding on whether to generate "pure clusters" or fewer clusters but accepting noise. With this as motivation, we verify the performance of several widely adopted algorithms on proteomic datasets and introduce a theoretical framework for drawing conclusions on which approach is suitable for the task at hand., (Copyright © 2021 Elsevier B.V. All rights reserved.)

Bacterial quorum sensing is the communication that takes place between bacteria as they secrete certain molecules into the intercellular medium that later get absorbed by the secreting cells themselves and by others. Depending on cell density, this uptake has the potential to alter gene expression and thereby affect global properties of the community. We consider the case of multiple bacterial species coexisting, referring to each one of them as a genotype and adopting the usual denomination of the molecules they collectively secrete as public goods. A crucial problem in this setting is characterizing the coevolution of genotypes as some of them secrete public goods (and pay the associated metabolic costs) while others do not but may nevertheless benefit from the available public goods. We introduce a network model to describe genotype interaction and evolution when genotype fitness depends on the production and uptake of public goods. The model comprises a random graph to summarize the possible evolutionary pathways the genotypes may take as they interact genetically with one another, and a system of coupled differential equations to characterize the behavior of genotype abundance in time. We study some simple variations of the model analytically and more complex variations computationally. Our results point to a simple trade-off affecting the long-term survival of those genotypes that do produce public goods. This trade-off involves, on the producer side, the impact of producing and that of absorbing the public good. On the nonproducer side, it involves the impact of absorbing the public good as well, now compounded by the molecular compatibility between the producer and the nonproducer. Depending on how these factors turn out, producers may or may not survive.

We present the Mixed-Data Acquisition (MDA) strategy for mass spectrometry data acquisition. MDA combines Data-Dependent Acquisition (DDA) and Data-Independent Acquisition (DIA) in the same run, thus doing away with the requirements for separate DDA spectral libraries. MDA is a natural result from advances in mass spectrometry, such as high scan rates and multiple analyzers, and is tailored toward exploiting these features. We demonstrate MDA's effectiveness on a yeast proteome analysis by overcoming a common bottleneck for XIC-based label-free quantitation; namely, the coelution of precursors when m/z values cannot be distinguished. We anticipate that MDA will become the next mainstream data generation approach for proteomics. MDA can also serve as an orthogonal validation approach for DDA experiments. Specialized software for MDA data analysis is made available on the project's website., (Copyright © 2020 Elsevier B.V. All rights reserved.)

Motivation: We present the first tool for unbiased quality control of top-down proteomics datasets. Our tool can select high-quality top-down proteomics spectra, serve as a gateway for building top-down spectral libraries and, ultimately, improve identification rates., Results: We demonstrate that a twofold rate increase for two E. coli top-down proteomics datasets may be achievable., Availability and Implementation: http://patternlabforproteomics.org/tdgc, freely available for academic use., Supplementary Information: Supplementary data are available at Bioinformatics online., (© The Author(s) 2019. Published by Oxford University Press. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

We present a new module integrated into the widely adopted PatternLab for proteomics to enable analysis of isotope-labeled peptides produced using dimethyl or SILAC. The accurate quantitation of proteins lies within the heart of proteomics; dimethylation has shown to be reliable, inexpensive, and applicable to any sample type. We validate our algorithm using an M. tuberculosis dataset obtained from two biological conditions; we used three dimethyl labels, one serving as an internal control for labeling a mixture of samples from both biological conditions. This internal control certified the proper functioning of our software. Availability: http://patternlabforproteomics.org, freely available for academic use., (Copyright © 2019. Published by Elsevier B.V.)

Background: Worldwide, breast cancer is the main cause of cancer mortality in women. Most cases originate in mammary ductal cells that produce the nipple aspirate fluid (NAF). In cancer patients, this secretome contains proteins associated with the tumor microenvironment. NAF studies are challenging because of inter-individual variability. We introduced a paired-proteomic shotgun strategy that relies on NAF analysis from both breasts of patients with unilateral breast cancer and extended PatternLab for Proteomics software to take advantage of this setup., Methods: The software is based on a peptide-centric approach and uses the binomial distribution to attribute a probability for each peptide as being linked to the disease; these probabilities are propagated to a final protein p-value according to the Stouffer's Z-score method., Results: A total of 1227 proteins were identified and quantified, of which 87 were differentially abundant, being mainly involved in glycolysis (Warburg effect) and immune system activation (activated stroma). Additionally, in the estrogen receptor-positive subgroup, proteins related to the regulation of insulin-like growth factor transport and platelet degranulation displayed higher abundance, confirming the presence of a proliferative microenvironment., Conclusions: We debuted a differential bioinformatics workflow for the proteomic analysis of NAF, validating this secretome as a treasure-trove for studying a paired-organ cancer type.

Analyzing the information content of DNA, though holding the promise to help quantify how the processes of evolution have led to information gain throughout the ages, has remained an elusive goal. Paradoxically, one of the main reasons for this has been precisely the great diversity of life on the planet: if on the one hand this diversity is a rich source of data for information-content analysis, on the other hand there is so much variation as to make the task unmanageable. During the past decade or so, however, succinct fragments of the COI mitochondrial gene, which is present in all animal phyla and in a few others, have been shown to be useful for species identification through DNA barcoding. A few million such fragments are now publicly available through the BOLD systems initiative, thus providing an unprecedented opportunity for relatively comprehensive information-theoretic analyses of DNA to be attempted. Here we show how a generalized form of total correlation can yield distinctive information-theoretic descriptors of the phyla represented in those fragments. In order to illustrate the potential of this analysis to provide new insight into the evolution of species, we performed principal component analysis on standardized versions of the said descriptors for 23 phyla. Surprisingly, we found that, though based solely on the species represented in the data, the first principal component correlates strongly with the natural logarithm of the number of all known living species for those phyla. The new descriptors thus constitute clear information-theoretic signatures of the processes whereby evolution has given rise to current biodiversity, which suggests their potential usefulness in further related studies., (Copyright © 2018 Elsevier Ltd. All rights reserved.)

Cross-linking coupled with mass spectrometry (XL-MS) has emerged as a powerful strategy for the identification of protein-protein interactions, characterization of interaction regions, and obtainment of structural information on proteins and protein complexes. In XL-MS, proteins or complexes are covalently stabilized with cross-linkers and digested, followed by identification of the cross-linked peptides by tandem mass spectrometry (MS/MS). This provides spatial constraints that enable modeling of protein (complex) structures and regions of interaction. However, most XL-MS approaches are not capable of differentiating intramolecular from intermolecular links in multimeric complexes, and therefore they cannot be used to study homodimer interfaces. We have recently developed an approach that overcomes this limitation by stable isotope-labeling of one of the two monomers, thereby creating a homodimer with one 'light' and one 'heavy' monomer. Here, we describe a step-by-step protocol for stable isotope-labeling, followed by controlled denaturation and refolding in the presence of the wild-type protein. The resulting light-heavy dimers are cross-linked, digested, and analyzed by mass spectrometry. We show how to quantitatively analyze the corresponding data with SIM-XL, an XL-MS software with a module tailored toward the MS/MS data from homodimers. In addition, we provide a video tutorial of the data analysis with this protocol. This protocol can be performed in ∼14 d, and requires basic biochemical and mass spectrometry skills.

Venoms are a rich source for the discovery of molecules with biotechnological applications, but their analysis is challenging even for state-of-the-art proteomics. Here we report on a large-scale proteomic assessment of the venom of Loxosceles intermedia, the so-called brown spider. Venom was extracted from 200 spiders and fractioned into two aliquots relative to a 10 kDa cutoff mass. Each of these was further fractioned and digested with trypsin (4 h), trypsin (18 h), pepsin (18 h), and chymotrypsin (18 h), then analyzed by MudPIT on an LTQ-Orbitrap XL ETD mass spectrometer fragmenting precursors by CID, HCD, and ETD. Aliquots of undigested samples were also analyzed. Our experimental design allowed us to apply spectral networks, thus enabling us to obtain meta-contig assemblies, and consequently de novo sequencing of practically complete proteins, culminating in a deep proteome assessment of the venom. Data are available via ProteomeXchange, with identifier PXD005523.

Motivation: Around 75% of all mass spectra remain unidentified by widely adopted proteomic strategies. We present DiagnoProt, an integrated computational environment that can efficiently cluster millions of spectra and use machine learning to shortlist high-quality unidentified mass spectra that are discriminative of different biological conditions., Results: We exemplify the use of DiagnoProt by shortlisting 4366 high-quality unidentified tandem mass spectra that are discriminative of different types of the Aspergillus fungus., Availability and Implementation: DiagnoProt, a demonstration video and a user tutorial are available at http://patternlabforproteomics.org/diagnoprot ., Contact: andrerfsilva@gmail.com or paulo@pcarvalho.com., Supplementary Information: Supplementary data are available at Bioinformatics online., (© The Author 2017. Published by Oxford University Press. All rights reserved. For Permissions, please e-mail: journals.permissions@oup.com)

Background: Promoting physical activity (PA) in low- and middle-income countries is an important public health topic as well as a challenge for practice. This study aimed to assess the effect of a school-based intervention on different PA-related variables among students., Methods: This cluster-randomized-controlled trial included 548 students in the intervention group and 537 in the control group (11-18 years-old) from 6 schools in neighborhoods with low Human Development Index (0.170-0.491) in Fortaleza, Brazil. The intervention included strategies focused on training teachers, opportunities for PA in the school environment and health education. Variables measured at baseline and again at the 4-months follow-up included the weekly time in different types of moderate-to-vigorous PA (MVPA), preference for PA during leisure-time, PA behavioral change stage and active commuting to school. Generalized linear models and binary logistic regressions were used., Results: An intervention effect was found by increasing the weekly time in MVPA (effect size = 0.17), popular games (effect size = 0.35), and the amount of PA per week (effect size = 0.27) among students (all P < .05)., Conclusions: The intervention was effective in promoting improvements in some PA outcomes, but the changes were not sufficient to increase the proportion of those meeting PA recommendations.

PatternLab for proteomics is an integrated computational environment that unifies several previously published modules for the analysis of shotgun proteomic data. The contained modules allow for formatting of sequence databases, peptide spectrum matching, statistical filtering and data organization, extracting quantitative information from label-free and chemically labeled data, and analyzing statistics for differential proteomics. PatternLab also has modules to perform similarity-driven studies with de novo sequencing data, to evaluate time-course experiments and to highlight the biological significance of data with regard to the Gene Ontology database. The PatternLab for proteomics 4.0 package brings together all of these modules in a self-contained software environment, which allows for complete proteomic data analysis and the display of results in a variety of graphical formats. All updates to PatternLab, including new features, have been previously tested on millions of mass spectra. PatternLab is easy to install, and it is freely available from http://patternlabforproteomics.org.

The production of structurally significant product ions during the dissociation of phosphopeptides is a key to the successful determination of phosphorylation sites. These diagnostic ions can be generated using the widely adopted MS/MS approach, MS3 (Data Dependent Neutral Loss - DDNL), or by multistage activation (MSA). The main purpose of this work is to introduce a false-localization rate (FLR) probabilistic model to enable unbiased phosphoproteomics studies. Briefly, our algorithm infers a probabilistic function from the distribution of the identified phosphopeptides' XCorr Delta scores (XD-Scores) in the current experiment. Our module infers p-values by relying on Gaussian mixture models and a logistic function. We demonstrate the usefulness of our probabilistic model by revisiting the "to MSA, or not to MSA" dilemma. For this, we use human leukemia-derived cells (K562) as a study model and enriched for phosphopeptides using the hydroxyapatite (HAP) chromatography. The aliquots were analyzed with and without MSA on an Orbitrap-XL. Our XD-Scoring analysis revealed that the MS/MS approach provides more identifications because of its faster scan rate, but that for the same given scan rate higher-confidence spectra can be achieved with MSA. Our software is integrated into the PatternLab for proteomics freely available for academic community at http://www.patternlabforproteomics.org. Biological significance Assigning statistical confidence to phosphorylation sites is necessary for proper phosphoproteomic assessment. Here we present a rigorous statistical model, based on Gaussian mixture models and a logistic function, which overcomes shortcomings of previous tools. The algorithm described herein is made readily available to the scientific community by integrating it into the widely adopted PatternLab for proteomics. This article is part of a Special Issue entitled: Computational Proteomics., (Copyright © 2015 Elsevier B.V. All rights reserved.)

Chemical cross-linking has emerged as a powerful approach for the structural characterization of proteins and protein complexes. However, the correct identification of covalently linked (cross-linked or XL) peptides analyzed by tandem mass spectrometry is still an open challenge. Here we present SIM-XL, a software tool that can analyze data generated through commonly used cross-linkers (e.g., BS3/DSS). Our software introduces a new paradigm for search-space reduction, which ultimately accounts for its increase in speed and sensitivity. Moreover, our search engine is the first to capitalize on reporter ions for selecting tandem mass spectra derived from cross-linked peptides. It also makes available a 2D interaction map and a spectrum-annotation tool unmatched by any of its kind. We show SIM-XL to be more sensitive and faster than a competing tool when analyzing a data set obtained from the human HSP90. The software is freely available for academic use at http://patternlabforproteomics.org/sim-xl. A video demonstrating the tool is available at http://patternlabforproteomics.org/sim-xl/video. SIM-XL is the first tool to support XL data in the mzIdentML format; all data are thus available from the ProteomeXchange consortium (identifier PXD001677). This article is part of a Special Issue entitled: Computational Proteomics., (Copyright © 2015 Elsevier B.V. All rights reserved.)

PepExplorer aids in the biological interpretation of de novo sequencing results; this is accomplished by assembling a list of homolog proteins obtained by aligning results from widely adopted de novo sequencing tools against a target-decoy sequence database. Our tool relies on pattern recognition to ensure that the results satisfy a user-given false-discovery rate (FDR). For this, it employs a radial basis function neural network that considers the precursor charge states, de novo sequencing scores, the peptide lengths, and alignment scores. PepExplorer is recommended for studies addressing organisms with no genomic sequence available. PepExplorer is integrated into the PatternLab for proteomics environment, which makes available various tools for downstream data analysis, including the resources for quantitative and differential proteomics., (Copyright © 2015 John Wiley & Sons, Inc.)

This work introduces a new methodology for the early detection of epileptic seizure based on the WiSARD weightless neural network model and a new approach in terms of preprocessing the electroencephalogram (EEG) data. WiSARD has, among other advantages, the capacity of perform the training phase in a very fast way. This speed in training is due to the fact that WiSARD's neurons work like Random Access Memories (RAM) addressed by input patterns. Promising results were obtained in the anticipation of seizure onsets in four representative patients from the European Database on Epilepsy (EPILEPSIAE). The proposed seizure early detection WNN architecture was explored by varying the detection anticipation (δ) in the 2 to 30 seconds interval, and by adopting 2 and 3 seconds as the width of the Sliding Observation Window (SOW) input. While in the most challenging patient (A) one obtained accuracies from 99.57% (δ=2s; SOW=3s) to 72.56% (δ=30s; SOW=2s), patient D seizures could be detected in the 99.77% (δ=2s; SOW=2s) to 99.93% (δ=30s; SOW=3s) accuracy interval.

Objective: To describe the percentile values for body mass index (BMI), waist circumference (WC) and waist-to-height (WHtR) of children from Colombo, Brazil, and compare them with data of children from other countries., Methods: This was a cross-sectional study with a random sample of 2,035 children aged 6-11 years. Age-and sex-specific smoothed percentiles curves for BMI, WC and WHtR were created using the LMS method. Values of 10(th), 50(th) and 90(th) percentiles from Brazilian children were compared with data from other countries., Results: There was a trend of increasing BMI and WC with age in both sexes. WHtR remained constant with advancing age in boys and girls. Comparison of the growth pattern among countries showed clear differences. Southern Brazil boys and girls had elevated 90(th) percentile values for BMI, which was similar to German children and higher than the North American and World Health Organization percentile values. However, children from this study had intermediate values for WC and WHtR in comparison to children from other countries., Conclusions: Elevated BMI values were observed among southern Brazilian children, but WC and WHtR percentile values were lower in southern Brazilian children than in children from other countries. Interventions at different levels should be made to avoid a probable increase of nutritional disorders (especially general obesity) in the next years., (Copyright © 2014 Associação de Pediatria de São Paulo. Publicado por Elsevier Editora Ltda. All rights reserved.)

Peptide spectrum matching is the current gold standard for protein identification via mass-spectrometry-based proteomics. Peptide spectrum matching compares experimental mass spectra against theoretical spectra generated from a protein sequence database to perform identification, but protein sequences not present in a database cannot be identified unless their sequences are in part conserved. The alternative approach, de novo sequencing, can make it possible to infer a peptide sequence directly from a mass spectrum, but interpreting long lists of peptide sequences resulting from large-scale experiments is not trivial. With this as motivation, PepExplorer was developed to use rigorous pattern recognition to assemble a list of homologue proteins using de novo sequencing data coupled to sequence alignment to allow biological interpretation of the data. PepExplorer can read the output of various widely adopted de novo sequencing tools and converge to a list of proteins with a global false-discovery rate. To this end, it employs a radial basis function neural network that considers precursor charge states, de novo sequencing scores, peptide lengths, and alignment scores to select similar protein candidates, from a target-decoy database, usually obtained from phylogenetically related species. Alignments are performed using a modified Smith-Waterman algorithm tailored for the task at hand. We verified the effectiveness of our approach using a reference set of identifications generated by ProLuCID when searching for Pyrococcus furiosus mass spectra on the corresponding NCBI RefSeq database. We then modified the sequence database by swapping amino acids until ProLuCID was no longer capable of identifying any proteins. By searching the mass spectra using PepExplorer on the modified database, we were able to recover most of the identifications at a 1% false-discovery rate. Finally, we employed PepExplorer to disclose a comprehensive proteomic assessment of the Bothrops jararaca plasma, a known biological source of natural inhibitors of snake toxins. PepExplorer is integrated into the PatternLab for Proteomics environment, which makes available various tools for downstream data analysis, including resources for quantitative and differential proteomics., (© 2014 by The American Society for Biochemistry and Molecular Biology, Inc.)

Accessing localized proteomic profiles has emerged as a fundamental strategy to understand the biology of diseases, as recently demonstrated, for example, in the context of determining cancer resection margins with improved precision. Here, we analyze a gastric cancer biopsy sectioned into 10 parts, each one subjected to MudPIT analysis. We introduce a software tool, named Shotgun Imaging Analyzer and inspired in MALDI imaging, to enable the overlaying of a protein's expression heat map on a tissue picture. The software is tightly integrated with the NeXtProt database, so it enables the browsing of identified proteins according to chromosomes, quickly listing human proteins never identified by mass spectrometry (i.e., the so-called missing proteins), and the automatic search for proteins that are more expressed over a specific region of interest on the biopsy, all of which constitute goals that are clearly well-aligned with those of the C-HPP. Our software has been able to highlight an intense expression of proteins previously known to be correlated with cancers (e.g., glutathione S-transferase Mu 3), and in particular, we draw attention to Gastrokine-2, a "missing protein" identified in this work of which we were able to clearly delineate the tumoral region from the "healthy" with our approach. Data are available via ProteomeXchange with identifier PXD000584.

Mass-spectrometry-based shotgun proteomics has become a widespread technology for analyzing complex protein mixtures. Here we describe a new module integrated into PatternLab for Proteomics that allows the pinpointing of differentially expressed domains. This is accomplished by inferring functional domains through our cloud service, using HMMER3 and Pfam remotely, and then mapping the quantitation values into domains for downstream analysis. In all, spotting which functional domains are changing when comparing biological states serves as a complementary approach to facilitate the understanding of a system's biology. We exemplify the new module's use by reanalyzing a previously published MudPIT dataset of Cryptococcus gattii cultivated under iron-depleted and replete conditions. We show how the differential analysis of functional domains can facilitate the interpretation of proteomic data by providing further valuable insight., (Copyright © 2013 Elsevier B.V. All rights reserved.)

Summary: Protein identification by mass spectrometry is commonly accomplished using a peptide sequence matching search algorithm, whose sensitivity varies inversely with the size of the sequence database and the number of post-translational modifications considered. We present the Spectrum Identification Machine, a peptide sequence matching tool that capitalizes on the high-intensity b1-fragment ion of tandem mass spectra of peptides coupled in solution with phenylisotiocyanate to confidently sequence the first amino acid and ultimately reduce the search space. We demonstrate that in complex search spaces, a gain of some 120% in sensitivity can be achieved., Availability: All data generated and the software are freely available for academic use at http://proteomics.fiocruz.br/software/sim., Contact: paulo@pcarvalho.com, Supplementary Information: Supplementary data are available at Bioinformatics online.

Given two subsets A and B of nodes in a directed graph, the conduciveness of the graph from A to B is the ratio representing how many of the edges outgoing from nodes in A are incoming to nodes in B. When the graph's nodes stand for the possible solutions to certain problems of combinatorial optimization, choosing its edges appropriately has been shown to lead to conduciveness properties that provide useful insight into the performance of algorithms to solve those problems. Here we study the conduciveness of CA-rule graphs, that is, graphs whose node set is the set of all CA rules given a cell's number of possible states and neighborhood size. We consider several different edge sets interconnecting these nodes, both deterministic and random ones, and derive analytical expressions for the resulting graph's conduciveness toward rules having a fixed number of non-quiescent entries. We demonstrate that one of the random edge sets, characterized by allowing nodes to be sparsely interconnected across any Hamming distance between the corresponding rules, has the potential of providing reasonable conduciveness toward the desired rules. We conjecture that this may lie at the bottom of the best strategies known to date for discovering complex rules to solve specific problems, all of an evolutionary nature.

A strategy for treating cancer is to surgically remove the tumor together with a portion of apparently healthy tissue surrounding it, the so-called "resection margin", to minimize recurrence. Here, we investigate whether the proteomic profiles from biopsies of gastric cancer resection margins are indeed more similar to those from healthy tissue than from cancer biopsies. To this end, we analyzed biopsies using an offline MudPIT shotgun proteomic approach and performed label-free quantitation through a distributed normalized spectral abundance factor approach adapted for extracted ion chromatograms (XICs). A multidimensional scaling analysis revealed that each of those tissue-types is very distinct from each other. The resection margin presented several proteins previously correlated with cancer, but also other overexpressed proteins that may be related to tumor nourishment and metastasis, such as collagen alpha-1, ceruloplasmin, calpastatin, and E-cadherin. We argue that the resection margin plays a key role in Paget's "soil to seed" hypothesis, that is, that cancer cells require a special microenvironment to nourish and that understanding it could ultimately lead to more effective treatments.

PatternLab for proteomics is a self-contained computational environment for analyzing shotgun proteomic data. Recent improvements incorporate modules to facilitate the computational analysis, such as FastaDBXtractor for sequence database preparation and ProLuCID runner for simplifying and managing the protein identification search engine; modules for pushing the limits on proteomics standards, such as SEPro, which relies on a semi-labeled decoy approach for increasing confidence in filtering and organizing peptide spectrum matches; and modules with novel features, such as SEProQ for enabling label-free quantitation by extracted ion chromatograms according to a distributed normalized ion abundance factor approach (dNIAF). Existing modules were also improved, such as the TFold module for pinpointing differentially expressed proteins. These new modules are integrated into the previously described arsenal of tools for further data analysis. Here we provide detailed instructions for operating and understanding them., (© 2012 by John Wiley & Sons, Inc.)

A quasispecies is a set of interrelated genotypes that have reached a stationary state while evolving according to the usual Darwinian principles of selection and mutation. Quasispecies studies invariably assume that it is possible for any genotype to mutate into any other, but recent finds indicate that this assumption is not necessarily true. Here we revisit the traditional quasispecies theory by adopting a network structure to constrain the occurrence of mutations. Such structure is governed by a random-graph model, whose single parameter (a probability p) controls both the graph's density and the dynamics of mutation. We contribute two further modifications to the theory, one to account for the fact that different loci in a genotype may be differently susceptible to the occurrence of mutations, the other to allow for a more plausible description of the transition from adaptation to degeneracy of the quasispecies as p is increased. We give analytical and simulation results for the usual case of binary genotypes, assuming the fitness landscape in which a genotype's fitness decays exponentially with its Hamming distance to the wild type. These results support the theory's assertions regarding the adaptation of the quasispecies to the fitness landscape and also its possible demise as a function of p.

Unlabelled: We present an updated version of the TFold software for pinpointing differentially expressed proteins in shotgun proteomics experiments. Given an FDR bound, the updated approach uses a theoretical FDR estimator to maximize the number of identifications that satisfy both a fold-change cutoff that varies with the t-test P-value as a power law and a stringency criterion that aims to detect lowly abundant proteins. The new version has yielded significant improvements in sensitivity over the previous one., Availability: Freely available for academic use at http://pcarvalho.com/patternlab.

The search engine processor (SEPro) is a tool for filtering, organizing, sharing, and displaying peptide spectrum matches. It employs a novel three-tier Bayesian approach that uses layers of spectrum, peptide, and protein logic to lead the data to converge to a single list of reliable protein identifications. SEPro is integrated into the PatternLab for proteomics environment, where an arsenal of tools for analyzing shotgun proteomic data is provided. By using the semi-labeled decoy approach for benchmarking, we show that SEPro significantly outperforms a commercially available competitor., (© 2012 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

The decoy-database approach is currently the gold standard for assessing the confidence of identifications in shotgun proteomic experiments. Here, we demonstrate that what might appear to be a good result under the decoy-database approach for a given false-discovery rate could be, in fact, the product of overfitting. This problem has been overlooked until now and could lead to obtaining boosted identification numbers whose reliability does not correspond to the expected false-discovery rate. To overcome this, we are introducing a modified version of the method, termed a semi-labeled decoy approach, which enables the statistical determination of an overfitted result., (Copyright © 2011 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

The breakup of alkaline glass and alumina plates due to planar impacts on one of their lateral sides is studied. Particular attention is given to investigating the spatial location of the cracks within the plates. Analysis based on a phenomenological model suggests that bifurcations along the cracks' paths are more likely to take place closer to the impact region than far away from it, i.e., the bifurcation probability seems to lower as the perpendicular distance from the impacted lateral increases. It is also found that many observables are not sensitive to the plate material used in this work, as long as the fragment multiplicities corresponding to the fragmentation of the plates are similar. This gives support to the universal properties of the fragmentation process reported in previous experiments. However, even under the just mentioned circumstances, some spatial observables are capable of distinguishing the material of which the plates are made, which therefore suggests that this universality should be carefully investigated.

An information-theoretic framework known as integrated information theory (IIT) has been introduced recently for the study of the emergence of consciousness in the brain [D. Balduzzi and G. Tononi, PLoS Comput. Biol. 4, e1000091 (2008)]. IIT purports that this phenomenon is to be equated with the generation of information by the brain surpassing the information that the brain's constituents already generate independently of one another. IIT is not fully plausible in its modeling assumptions nor is it testable due to severe combinatorial growth embedded in its key definitions. Here we introduce an alternative to IIT which, while inspired in similar information-theoretic principles, seeks to address some of IIT's shortcomings to some extent. Our alternative framework uses the same network-algorithmic cortical model we introduced earlier [A. Nathan and V. C. Barbosa, Phys. Rev. E 81, 021916 (2010)] and, to allow for somewhat improved testability relative to IIT, adopts the well-known notions of information gain and total correlation applied to a set of variables representing the reachability of neurons by messages in the model's dynamics. We argue that these two quantities relate to each other in such a way that can be used to quantify the system's efficiency in generating information beyond that which does not depend on integration. We give computational results on our cortical model and on variants thereof that are either structurally random in the sense of an Erdős-Rényi random directed graph or structurally deterministic. We have found that our cortical model stands out with respect to the others in the sense that many of its instances are capable of integrating information more efficiently than most of those others' instances.

Summary: We present an approach to statistically pinpoint differentially expressed proteins that have quantitation values near the quantitation threshold and are not identified in all replicates (marginal cases). Our method uses a Bayesian strategy to combine parametric statistics with an empirical distribution built from the reproducibility quality of the technical replicates., Availability: The software is freely available for academic use at http://pcarvalho.com/patternlab.

The Tri-Border Region has several social and health problems among young people, however, there are no data about the overweight between adolescents. This study investigated the prevalence of overweight and associated factors in students from Tri-Border Region: Argentina, Brazil and Paraguay. Study participants were 1.183 students, ages from 15 to 18 years. The overweight was identified according to body mass index cut-off points proposed by the World Health Organization. The Global School-Based Student Health Survey questionnaire was used to identify sociodemographic (gender and age) and behavioral factors (physical activity outside of school, commuting to school, time watching TV and fruits, vegetables, sweet, and salty snack consumption) associated with overweight in adolescents. It was used the descriptive statistics, analysis of variance, chi-square and binary logistic regression, adopting p < 0.05. Thirteen percent of students were overweight. Boys were approximately two times more likely to have overweight than girls, independently of nationality. Brazilian students that realized passive commuting to school and Argentineans students with low consumption of vegetables (< 1 time/day) were 2.2 and 2.9 times more likely to have overweight than their counterparts who performed active commuting to school and consumed vegetables daily, respectively. These results suggest that public policies to combat overweight should attention on promoting healthy lifestyle among young people from Tri-Border Region.

In evolutionary dynamics, the probability that a mutation spreads through the whole population, having arisen from a single individual, is known as the fixation probability. In general, it is not possible to find the fixation probability analytically given the mutant's fitness and the topological constraints that govern the spread of the mutation, so one resorts to simulations instead. Depending on the topology in use, a great number of evolutionary steps may be needed in each of the simulation events, particularly in those that end with the population containing mutants only. We introduce two techniques to accelerate the determination of the fixation probability. The first one skips all evolutionary steps in which the number of mutants does not change and thereby reduces the number of steps per simulation event considerably. This technique is computationally advantageous for some of the so-called layered networks. The second technique, which is not restricted to layered networks, consists of aborting any simulation event in which the number of mutants has grown beyond a certain threshold value and counting that event as having led to a total spread of the mutation. For advantageous mutations in large populations and regardless of the network's topology, we demonstrate, both analytically and by means of simulations, that using a threshold of about [N/(r-1)](1/4) mutants, where N is the number of simulation events and r is the ratio of the mutants' fitness to that of the remainder of the population, leads to an estimate of the fixation probability that deviates in no significant way from that obtained from the full-fledged simulations. We have observed speedups of two orders of magnitude for layered networks with 10,000 nodes.

Background: Given an undirected graph, we consider the two problems of combinatorial optimization, which ask that its chromatic and independence numbers be found. Although both problems are NP-hard, when either one is solved on the incrementally denser graphs of a random sequence, at certain critical values of the number of edges, it happens that the transition to the next value causes optimal solutions to be obtainable substantially more easily than right before it., Methodology/principal Findings: We introduce the notion of a network's conduciveness, a probabilistically interpretable measure of how the network's structure allows it to be conducive to roaming agents, in certain conditions, from one portion of the network to another. We demonstrate that the performance jumps of graph coloring and independent sets at the critical-value transitions in the number of edges can be understood by resorting to the network that represents the solution space of the problems for each graph and examining its conduciveness between the non-optimal solutions and the optimal ones. Right past each transition, this network becomes strikingly more conducive in the direction of the optimal solutions than it was just before it, while at the same time becoming less conducive in the opposite direction., Conclusions/significance: Network conduciveness provides a useful conceptual framework for explaining the performance jumps associated with graph coloring and independent sets. We believe it may also become instrumental in helping clarify further issues related to NP-hardness that remain poorly understood. Additionally, it may become useful also in other areas in which network theory has a role to play.

PatternLab for proteomics is a one-stop shop computational environment for analyzing shotgun proteomic data. Its modules provide means to pinpoint proteins/peptides that are differentially expressed and those that are unique to a state. It can also cluster the ones that share similar expression profiles in time-course experiments, as well as help in interpreting results according to Gene Ontology. PatternLab is user-friendly, simple, and provides a graphical user interface., ((c) 2010 by John Wiley & Sons, Inc.)

The fragmentation of alumina and glass plates due to lateral impact is studied. A few hundred plates have been fragmented at different impact velocities and the produced fragments are analyzed. The method employed in this work allows one to investigate some geometrical properties of the fragments, besides the traditional size distribution usually studied in former experiments. We found that, although both materials exhibit qualitative similar fragment size distribution function, their geometrical properties appear to be quite different. A schematic model for two-dimensional fragmentation is also presented and its predictions are compared to our experimental results. The comparison suggests that the analysis of the fragments' geometrical properties constitutes a more stringent test of the theoretical models' assumptions than the size distribution.

Summary: XDIA is a computational strategy for analyzing multiplexed spectra acquired using electron transfer dissociation and collision-activated dissociation; it significantly increases identified spectra (approximately 250%) and unique peptides (approximately 30%) when compared with the data-dependent ETCaD analysis on middle-down, single-phase shotgun proteomic analysis. Increasing identified spectra and peptides improves quantitation statistics confidence and protein coverage, respectively., Availability: The software and data produced in this work are freely available for academic use at http://fields.scripps.edu/XDIA CONTACT: paulo@pcarvalho.com, Supplementary Information: Supplementary data are available at Bioinformatics online.

Perillyl alcohol (POH) is a naturally occurring terpene and a promising chemotherapeutic agent for glioblastoma multiform; yet, little is known about its molecular effects. Here we present results of a semi-quantitative proteomic analysis of A172 cells exposed to POH for different time-periods (1', 10', 30', 60', 4h, and 24h). The analysis identified more than 4000 proteins; which were clustered using PatternLab for proteomics and then linked to Ras signaling, tissue homeostasis, induction of apoptosis, metallopeptidase activity, and ubiquitin-protein ligase activity. Our results make available one of the most complete protein repositories for the A172. Moreover, we detected the phosphorylation of GSK3beta (Glycogen synthase kinase) and the inhibition of ERK's (extracellular signal regulated kinase) phosphorylation after 10', which suggests a new mechanism of POH's activation for apoptosis., (Copyright (c) 2010 Elsevier B.V. All rights reserved.)

When a neuron fires and the resulting action potential travels down its axon toward other neurons' dendrites, the effect on each of those neurons is mediated by the strength of the synapse that separates it from the firing neuron. This strength, in turn, is affected by the postsynaptic neuron's response through a mechanism that is thought to underlie important processes such as learning and memory. Although of difficult quantification, cortical synaptic strengths have been found to obey a long-tailed unimodal distribution peaking near the lowest values (approximately lognormal), thus confirming some of the predictive models built previously. Most of these models are causally local, in the sense that they refer to the situation in which a number of neurons all fire directly at the same postsynaptic neuron. Consequently, they necessarily embody assumptions regarding the generation of action potentials by the presynaptic neurons that have little biological interpretability. We introduce a network model of large groups of interconnected neurons and demonstrate, making none of the assumptions that characterize the causally local models, that its long-term behavior gives rise to a distribution of synaptic weights (the mathematical surrogates of synaptic strengths) with the same properties that were experimentally observed. In our model, the action potentials that create a neuron's input are, ultimately, the product of network-wide causal chains relating what happens at a neuron to the firings of others. Our model is then of a causally global nature and predicates the emergence of the synaptic-weight distribution on network structure and function. As such, it has the potential to become instrumental also in the study of other emergent cortical phenomena.

Unlabelled: YADA can deisotope and decharge high-resolution mass spectra from large peptide molecules, link the precursor monoisotopic peak information to the corresponding tandem mass spectrum, and account for different co-fragmenting ion species (multiplexed spectra). We describe how YADA enables a pipeline consisting of ProLuCID and DTASelect for analyzing large-scale middle-down proteomics data., Availability: http://fields.scripps.edu/yada

Natural selection and random drift are competing phenomena for explaining the evolution of populations. Combining a highly fit mutant with a population structure that improves the odds that the mutation spreads through the whole population tips the balance in favor of natural selection. The probability that the spread occurs, known as the fixation probability, depends heavily on how the population is structured. Certain topologies, albeit highly artificially contrived, have been shown to exist that favor fixation. We present a randomized mechanism for network growth that is loosely inspired in some of these topologies' key properties and demonstrate, through simulations, that it is capable of giving rise to structured populations for which the fixation probability significantly surpasses that of an unstructured population. This discovery provides important support to the notion that natural selection can be enhanced over random drift in naturally occurring population structures.

Electron transfer dissociation (ETD) can dissociate highly charged ions. Efficient analysis of ions dissociated with ETD requires accurate determination of charge states for calculation of molecular weight. We created an algorithm to assign the charge state of ions often used for ETD. The program, Charge Prediction Machine (CPM), uses Bayesian decision theory to account for different charge reduction processes encountered in ETD and can also handle multiplex spectra. CPM correctly assigned charge states to 98% of the 13,097 MS2 spectra from a combined data set of four experiments. In a comparison between CPM and a competing program, Charger (ThermoFisher), CPM produced half the mistakes.