Preparing to Work with Big Data in Proteomics – A Report on the HUPO‐PSI Spring Workshop

The Human Proteome Organisation Proteomics Standards Initiative (HUPO‐PSI) was established in 2002 with the aim of defining community standards for data representation in proteomics and facilitating data comparison, exchange and verification. The 2013 annual spring workshop was hosted by the University of Liverpool, UK and concentrated on updating and refining the existing standards in the light of new methodologies and technologies. To control the inflation of file sizes, strategies for file compression, particularly for mzML files, were explored. Best practices for encoding information such as protein grouping and PTM localisation were refined and documented. Additional example files for the mzQuantML format were designed to provide support for selected reaction monitoring techniques. Enhancements to the PSI Common Query Interface (PSICQUIC) and PSI‐MI XML were discussed. Finally, the group engaged in discussion on how the existing work of the HUPO‐PSI can be leveraged by the Metabolomics Standards Initiative to improve the capture of metabolite data.     

HTAPP: High‐throughput autonomous proteomic pipeline

Recent advances in the speed and sensitivity of mass spectrometers and in analytical methods, the exponential acceleration of computer processing speeds, and the availability of genomic databases from an array of species and protein information databases have led to a deluge of proteomic data. The development of a lab‐based automated proteomic software platform for the automated collection, processing, storage, and visualization of expansive proteomic data sets is critically important. The high‐throughput autonomous proteomic pipeline described here is designed from the ground up to provide critically important flexibility for diverse proteomic workflows and to streamline the total analysis of a complex proteomic sample. This tool is composed of a software that controls the acquisition of mass spectral data along with automation of post‐acquisition tasks such as peptide quantification, clustered MS/MS spectral database searching, statistical validation, and data exploration within a user‐configurable lab‐based relational database. The software design of high‐throughput autonomous proteomic pipeline focuses on accommodating diverse workflows and providing missing software functionality to a wide range of proteomic researchers to accelerate the extraction of biological meaning from immense proteomic data sets. Although individual software modules in our integrated technology platform may have some similarities to existing tools, the true novelty of the approach described here is in the synergistic and flexible combination of these tools to provide an integrated and efficient analysis of proteomic samples.     

Meeting New Challenges: The 2014 HUPO‐PSI/COSMOS Workshop

The Annual 2014 Spring Workshop of the Proteomics Standards Initiative (PSI) of the Human Proteome Organization (HUPO) was held this year jointly with the metabolomics COordination of Standards in MetabOlomicS (COSMOS) group. The range of existing MS standards (mzML, mzIdentML, mzQuantML, mzTab, TraML) was reviewed and updated in the light of new methodologies and advances in technologies. Adaptations to meet the needs of the metabolomics community were incorporated and a new data format for NMR, nmrML, was presented. The molecular interactions workgroup began work on a new version of the existing XML data interchange format. PSI‐MI XML3.0 will enable the capture of more abstract data types such as protein complex topology derived from experimental data, allosteric binding, and dynamic interactions. Further information about the work of the HUPO‐PSI can be found at http://www.psidev.info .     

Specific substitutions of light‐harvesting complex I proteins associated with photosystem I are required for supercomplex formation with chloroplast NADH dehydrogenase‐like complex

In Arabidopsis, the chloroplast NADH‐dehydrogenase‐like (NDH) complex is sandwiched between two copies of photosystem I (PSI) supercomplex, consisting of a PSI core and four light‐harvesting complex I (LHCI) proteins (PSI‐LHCI) to form the NDH–PSI supercomplex. Two minor LHCI proteins, Lhca5 and Lhca6, contribute to the interaction of each PSI–LHCI copy with the NDH complex. Here, large‐pore blue‐native gel electrophoresis revealed that, in addition to this complex, there were at least two types of higher‐order association of more LHCI copies with the NDH complex. In single‐particle images, this higher‐order association of PSI–LHCI preferentially occurs at the left side of the NDH complex when viewed from the stromal side, placing subcomplex A at the top (Yadav et al., Biochim. Biophys. Acta ‐ Bioenerg., 1858, 2017, 12). The association was impaired in the lhca6 mutant but not in the lhca5 mutant, suggesting that the left copy of PSI–LHCI was linked to the NDH complex via Lhca6. From an analysis of subunit compositions of the NDH–PSI supercomplex in lhca5 and lhca6 mutants, we propose that Lhca6 substitutes for Lhca2 in the left copy of PSI–LHCI, whereas Lhca5 substitutes for Lhca4 in the right copy. In the lhca2 mutant, Lhca3 was specifically stabilized in the NDH–PSI supercomplex through heterodimer formation with Lhca6. In the left copy of PSI–LHCI, subcomplex B, Lhca6 and NdhD likely formed the core of the supercomplex interaction. In contrast, a larger protein complex, including at least subcomplexes B and L and NdhB, was needed to form the contact site with Lhca5 in the right copy of PSI–LHCI.     

Exploring skyline for both MSE‐based label‐free proteomics and HRMS quantitation of small molecules

The MS^E (where MS^E is low energy (MS) and elevated energy (E) mode of acquisition) acquisition method commercialized by Waters on its Q‐TOF instruments is regarded as a unique data‐independent fragmentation approach that improves the accuracy and dynamic range of label‐free proteomic quantitation. Due to its special format, MS^E acquisition files cannot be independently analyzed with most widely used open‐source proteomic software specialized for processing data‐dependent acquisition files. In this study, we established a workflow integrating Skyline, a popular and versatile peptide‐centric quantitation program, and a statistical tool DiffProt to fulfill MS^E‐based proteomic quantitation. Comparison with the vendor software package for analyzing targeted phosphopeptides and global proteomic datasets reveals distinct advantages of Skyline in MS^E data mining, including sensitive peak detection, flexible peptide filtering, and transparent step‐by‐step workflow. Moreover, we developed a new procedure such that Skyline MS1 filtering was extended to small molecule quantitation for the first time. This new utility of Skyline was examined in a protein–ligand interaction experiment to identify multiple chemical compounds specifically bound to NDM‐1 (where NDM is New Delhi metallo‐β‐lactamase 1), an antibiotics‐resistance target. Further improvement of the current weaknesses in Skyline MS1 filtering is expected to enhance the reliability of this powerful program in full scan‐based quantitation of both peptides and small molecules.     

The composition and structure of photosystem I‐associated antenna from Cyanidioschyzon merolae

Red algae contain two types of light‐harvesting antenna systems, the phycobilisomes and chlorophyll a binding polypeptides (termed Lhcr), which expand the light‐harvesting capacity of the photosynthetic reaction centers. In this study, photosystem I (PSI) and its associated light‐harvesting proteins were isolated from the red alga Cyanidioschyzon merolae. The structural and functional properties of the largest PSI particles observed were investigated by biochemical characterization, mass spectrometry, fluorescence emission and excitation spectroscopy, and transmission electron microscopy. Our data provide strong evidence for a stable PSI complex in red algae that possesses two distinct types of functional peripheral light‐harvesting antenna complex, comprising both Lhcr and a PSI‐linked phycobilisome sub‐complex. We conclude that the PSI antennae system of red algae represents an evolutionary intermediate between the prokaryotic cyanobacteria and other eukaryotes, such as green algae and vascular plants.     

PROTEOMER: A workflow‐optimized laboratory information management system for 2‐D electrophoresis‐centered proteomics

In recent years proteomics became increasingly important to functional genomics. Although a large amount of data is generated by high throughput large‐scale techniques, a connection of these mostly heterogeneous data from different analytical platforms and of different experiments is limited. Data mining procedures and algorithms are often insufficient to extract meaningful results from large datasets and therefore limit the exploitation of the generated biological information. In our proteomic core facility, which almost exclusively focuses on 2‐DE/MS‐based proteomics, we developed a proteomic database custom tailored to our needs aiming at connecting MS protein identification information to 2‐DE derived protein expression profiles. The tools developed should not only enable an automatic evaluation of single experiments, but also link multiple 2‐DE experiments with MS‐data on different levels and thereby helping to create a comprehensive network of our proteomics data. Therefore the key feature of our “PROTEOMER” database is its high cross‐referencing capacity, enabling integration of a wide range of experimental data. To illustrate the workflow and utility of the system, two practical examples are provided to demonstrate that proper data cross‐referencing can transform information into biological knowledge.     

THE CYANOBACTERIAL CHLOROPHYLL‐BINDING‐PROTEIN ISIA ACTS TO INCREASE THE IN VIVO EFFECTIVE ABSORPTION CROSS‐SECTION OF PSI UNDER IRON LIMITATION1

Iron availability limits primary production in >30% of the world’s oceans; hence phytoplankton have developed acclimation strategies. In particular, cyanobacteria express IsiA (iron‐stress‐induced) under iron stress, which can become the most abundant chl‐binding protein in the cell. Within iron‐limited oceanic regions with significant cyanobacterial biomass, IsiA may represent a significant fraction of the total chl. We spectroscopically measured the effective cross‐section of the photosynthetic reaction center PSI (σ_PSI) in vivo and biochemically quantified the absolute abundance of PSI, PSII, and IsiA in the model cyanobacterium Synechocystis sp. PCC 6803. We demonstrate that accumulation of IsiA results in a ～60% increase in σ_PSI, in agreement with the theoretical increase in cross‐section based on the structure of the biochemically isolated IsiA‐PSI supercomplex from cyanobacteria. Deriving a chl budget, we suggest that IsiA plays a primary role as a light‐harvesting antenna for PSI. On progressive iron‐stress in culture, IsiA continues to accumulate without a concomitant increase in σ_PSI, suggesting that there may be a secondary role for IsiA. In natural populations, the potential physiological significance of the uncoupled pool of IsiA remains to be established. However, the functional role as a PSI antenna suggests that a large fraction of IsiA‐bound chl is directly involved in photosynthetic electron transport.     

Reconstructing a puzzle: existence of cyanophages containing both photosystem‐I and photosystem‐II gene suites inferred from oceanic metagenomic datasets

Cyanobacteria play a key role in marine photosynthesis, which contributes to the global carbon cycle and to the world oxygen supply. Genes encoding the photosystem‐II (PSII) reaction centre are found in many cyanophage genomes, and it was suggested that the horizontal transfer of these genes might be involved in increasing phage fitness. Recently, evidence for the existence of phages carrying Photosystem‐I (PSI) genes was also reported. Here, using a combination of different marine metagenomic datasets and a unique crossing of the datasets, we now describe the finding of phages that, as in plants and cyanobacteria, contain both PSII and PSI genes. In addition, these phages also contain NADH dehydrogenase genes. The presence of modified PSII and PSI genes in the same viral entities in combination with electron transfer proteins like NAD(P)H dehydrogenase (NDH‐1) strongly points to a role in perturbation of the cyanobacterial host photosynthetic electron flow. We therefore suggest that, depending on the physiological condition of the infected cyanobacterial host, the viruses may use different options to maximize survival. The modified PSI may alternate between functioning with PSII in linear electron transfer and contributing to the production of both NADPH and ATP or functioning independently of PSII in cyclic mode via the NDH‐1 complex and thus producing only ATP.     

Far‐red light enhances photochemical efficiency in a wavelength‐dependent manner

Linear electron transport depends on balanced excitation of photosystem I and II. Far‐red light preferentially excites photosystem I (PSI) and can enhance the photosynthetic efficiency when combined with light that over‐excites photosystem II (PSII). The efficiency of different wavelengths of far‐red light exciting PSI was quantified by measuring the change in quantum yield of PSII (Φ_PSII) of lettuce (Lactuca sativa) under red/blue light with narrowband far‐red light added (from 678 to 752 nm, obtained using laser diodes). The Φ_PSII of lettuce increased with increasing wavelengths of added light from 678 to 703 nm, indicating longer wavelengths within this region are increasingly used more efficiently by PSI than by PSII. Adding 721 nm light resulted in similar Φ_PSII as adding 703 nm light, but Φ_PSII tended to decrease as wavelength increased from 721 to 731 nm, likely due to decreasing absorptance and low photon energy. Adding 752 nm light did not affect Φ_PSII. Leaf chlorophyll fluorescence light response measurements showed lettuce had higher Φ_PSII under halogen light (rich in far‐red) than under red/blue light (which over‐excites PSII). Far‐red light is more photosynthetically active than commonly believed, because of its synergistic interaction with light of shorter wavelengths.     

Haloferax volcanii archaeosortase is required for motility,mating, and C‐terminal processing of the S‐layer glycoprotein

Cell surfaces are decorated by a variety of proteins that facilitate interactions with their environments and support cell stability. These secreted proteins are anchored to the cell by mechanisms that are diverse, and, in archaea, poorly understood. Recently published in silico data suggest that in some species a subset of secreted euryarchaeal proteins, which includes the S‐layer glycoprotein, is processed and covalently linked to the cell membrane by enzymes referred to as archaeosortases. In silico work led to the proposal that an independent, sortase‐like system for proteolysis‐coupled, carboxy‐terminal lipid modification exists in bacteria (exosortase) and archaea (archaeosortase). Here, we provide the first in vivo characterization of an archaeosortase in the haloarchaeal model organism Haloferax volcanii. Deletion of the artA gene (HVO_0915) resulted in multiple biological phenotypes: (a) poor growth, especially under low‐salt conditions, (b) alterations in cell shape and the S‐layer, (c) impaired motility, suppressors of which still exhibit poor growth, and (d) impaired conjugation. We studied one of the ArtA substrates, the S‐layer glycoprotein, using detailed proteomic analysis. While the carboxy‐terminal region of S‐layer glycoproteins, consisting of a putative threonine‐rich O‐glycosylated region followed by a hydrophobic transmembrane helix, has been notoriously resistant to any proteomic peptide identification, we were able to identify two overlapping peptides from the transmembrane domain present in the ΔartA strain but not in the wild‐type strain. This clearly shows that ArtA is involved in carboxy‐terminal post‐translational processing of the S‐layer glycoprotein. As it is known from previous studies that a lipid is covalently attached to the carboxy‐terminal region of the S‐layer glycoprotein, our data strongly support the conclusion that archaeosortase functions analogously to sortase, mediating proteolysis‐coupled, covalent cell surface attachment.     

Integrated platform for manual and high‐throughput statistical validation of tandem mass spectra

As proteomic data sets increase in size and complexity, the necessity for database‐centric software systems able to organize, compare, and visualize all the proteomic experiments in a lab grows. We recently developed an integrated platform called high‐throughput autonomous proteomic pipeline (HTAPP) for the automated acquisition and processing of quantitative proteomic data, and integration of proteomic results with existing external protein information resources within a lab‐based relational database called PeptideDepot. Here, we introduce the peptide validation software component of this system, which combines relational database‐integrated electronic manual spectral annotation in Java with a new software tool in the R programming language for the generation of logistic regression spectral models from user‐supplied validated data sets and flexible application of these user‐generated models in automated proteomic workflows. This logistic regression spectral model uses both variables computed directly from SEQUEST output in addition to deterministic variables based on expert manual validation criteria of spectral quality. In the case of linear quadrupole ion trap (LTQ) or LTQ‐FTICR LC/MS data, our logistic spectral model outperformed both XCorr (242% more peptides identified on average) and the X!Tandem E‐value (87% more peptides identified on average) at a 1% false discovery rate estimated by decoy database approach.     

Analysis of fast chlorophyll fluorescence rise (O‐K‐J‐I‐P) curves in green fruits indicates electron flow limitations at the donor side of PSII and the acceptor sides of both photosystems

Limited evidence up to now indicates low linear photosynthetic electron flow and CO₂ assimilation rates in non‐foliar chloroplasts. In this investigation, we used chlorophyll fluorescence techniques to locate possible limiting steps in photosystem function in exposed, non‐stressed green fruits (both pericarps and seeds) of three species, while corresponding leaves served as controls. Compared with leaves, fruit photosynthesis was characterized by less photon trapping and less quantum yields of electron flow, while the non‐photochemical quenching was higher and potentially linked to enhanced carotenoid/chlorophyll ratios. Analysis of fast chlorophyll fluorescence rise curves revealed possible limitations both in the donor (oxygen evolving complex) and the acceptor (Q_A^?→ intermediate carriers) sides of photosystem II (PSII) indicating innately low PSII photochemical activity. On the other hand, PSI was characterized by faster reduction of its final electron acceptors and their small pool sizes. We argue that the fast reductive saturation of final PSI electron acceptors may divert electrons back to intermediate carriers facilitating a cyclic flow around PSI, while the partial inactivation of linear flow precludes strong reduction of plastoquinone. As such, the photosynthetic attributes of fruit chloroplasts may act to replenish the ATP lost because of hypoxia usually encountered in sink organs with high diffusive resistance to gas exchange.     

multiplierz v2.0: A Python‐based ecosystem for shared access and analysis of native mass spectrometry data

The continued evolution of modern mass spectrometry instrumentation and associated methods represents a critical component in efforts to decipher the molecular mechanisms which underlie normal physiology and understand how dysregulation of biological pathways contributes to human disease. The increasing scale of these experiments combined with the technological diversity of mass spectrometers presents several challenges for community‐wide data access, analysis, and distribution. Here we detail a redesigned version of multiplierz, our Python software library which leverages our common application programming interface (mzAPI) for analysis and distribution of proteomic data. New features include support for a wider range of native mass spectrometry file types, interfaces to additional database search engines, compatibility with new reporting formats, and high‐level tools to perform post‐search proteomic analyses. A GUI desktop environment, mzDesktop, provides access to multiplierz functionality through a user friendly interface. multiplierz is available for download from: https://github.com/BlaisProteomics/multiplierz ; and mzDesktop is available for download from: https://sourceforge.net/projects/multiplierz/     

Addressing the Challenges of High‐Throughput Cancer Tissue Proteomics for Clinical Application: ProCan

The cancer tissue proteome has enormous potential as a source of novel predictive biomarkers in oncology. Progress in the development of mass spectrometry (MS)‐based tissue proteomics now presents an opportunity to exploit this by applying the strategies of comprehensive molecular profiling and big‐data analytics that are refined in other fields of ‘omics research. ProCan (ProCan is a registered trademark) is a program aiming to generate high‐quality tissue proteomic data across a broad spectrum of cancer types. It is based on data‐independent acquisition–MS proteomic analysis of annotated tissue samples sourced through collaboration with expert clinical and cancer research groups. The practical requirements of a high‐throughput translational research program have shaped the approach that ProCan is taking to address challenges in study design, sample preparation, raw data acquisition, and data analysis. The ultimate goal is to establish a large proteomics knowledge‐base that, in combination with other cancer ‘omics data, will accelerate cancer research.     

2‐D PAGE and MS analysis of proteins from formalin‐fixed,paraffin‐embedded tissues

In the past decade, encouraging results have been obtained in extraction and analysis of proteins from formalin‐fixed, paraffin‐embedded (FFPE) tissues. However, 2‐D PAGE protein maps with satisfactory proteomic information and comparability to fresh tissues have never been described to date. In the present study, we report 2‐D PAGE separation and MS identification of full‐length proteins extracted from FFPE skeletal muscle tissue. The 2‐D protein profiles obtained from FFPE tissues could be matched to those achieved from frozen tissues replicates. Up to 250 spots were clearly detected in 2‐D maps of proteins from FFPE tissue following standard mass‐compatible silver staining. Protein spots from both FFPE and frozen tissue 2‐D gels were excised, subjected to in situ hydrolysis, and identified by MS analysis. Matched spots produced matched protein identifications. Moreover, 2‐D protein maps from FFPE tissues were successfully subjected to Western immunoblotting, producing comparable results to fresh‐frozen tissues. In conclusion, this study provides evidence that, when adequately extracted, full‐length proteins from FFPE tissues might be suitable to 2‐D PAGE‐MS analysis, allowing differential proteomic studies on the vast existing archives of healthy and pathological‐fixed tissues.     

Possible target‐related proteins and signal network of bufalin in A549 cells suggested by both iTRAQ‐based and label‐free proteomic analysis

Bufalin (BF) exhibited antiproliferation and antimigration effects on human A549 lung cancer cells. To search its target‐related proteins, protein expression profiles of BF‐treated and control cells were compared using two quantitative proteomic methods, iTRAQ‐based and label‐free proteomic analysis. A total of 5428 proteins were identified in iTRAQ‐based analysis while 6632 proteins were identified in label‐free analysis. The number of common identified proteins of both methods was 4799 proteins. By application of 1.20‐fold for upregulated and 0.83‐fold for downregulated cutoff values, 273 and 802 differentially expressed proteins were found in iTRAQ‐based and label‐free analysis, respectively. The number of common differentially expressed proteins of both methods was 45 proteins. Results of bioinformational analysis using Metacore^TM showed that the two proteomic methods were complementary and both suggested the involvement of oxidative stress and regulation of gene expression in the effects of BF, and fibronectin‐related pathway was suggested to be an important pathway affected by BF. Western blotting assay results confirmed BF‐induced change in levels of fibronectin and other related proteins. Overexpression of fibronectin by plasmid transfection ameliorated antimigration effects of BF. Results of the present study provided information about possible target‐related proteins and signal network of BF.     

The secretion of the bacterial phytase PHY‐US417 by Arabidopsis roots reveals its potential for increasing phosphate acquisition and biomass production during co‐growth

Phytic acid (PA) is a major source of inorganic phosphate (Pi) in the soil; however, the plant lacks the capacity to utilize it for Pi nutrition and growth. Microbial phytases constitute a group of enzymes that are able to remobilize Pi from PA. Thus, the use of these phytases to increase the capacity of higher plants to remobilize Pi from PA is of agronomical interest. In the current study, we generate transgenic Arabidopsis lines (ePHY) overexpressing an extracellular form of the phytase PHY‐US417 of Bacillus subtilis, which are characterized by high levels of secreted phytase activity. In the presence of PA as sole source of Pi, while the wild‐type plants show hallmark of Pi deficiency phenotypes, including the induction of the expression of Pi starvation‐induced genes (PSI, e.g. PHT1;4) and the inhibition of growth capacity, the ePHY overexpressing lines show a higher biomass production and no PSI induction. Interestingly, when co‐cultured with ePHY overexpressors, wild‐type Arabidopsis plants (or tobacco) show repression of the PSI genes, improvement of Pi content and increases in biomass production. In line with these results, mutants in the high‐affinity Pi transporters, namely pht1;1 and pht1;1‐1;4, both fail to accumulate Pi and to grow when co‐cultured with ePHY overexpressors. Taken together, these data demonstrate the potential of secreted phytases in improving the Pi content and enhancing growth of not only the transgenic lines but also the neighbouring plants.     

Site‐specific inhibition of integrin αvβ3‐vitronectin association by a ser‐asp‐val sequence through an Arg‐Gly‐Asp‐binding site of the integrin

A functional proteomic technology using protein chip and molecular simulation was used to demonstrate a novel biomolecular interaction between P11, a peptide containing the Ser‐Asp‐Val (SDV) sequence and integrin αvβ3. P11 (HSDVHK) is a novel antagonistic peptide of integrin αvβ3 screened from hexapeptide library through protein chip system. An in silico docking study and competitive protein chip assay revealed that the SDV sequence of P11 is able to create a stable inhibitory complex onto the vitronectin‐binding site of integrin αvβ3. The Arg‐Gly‐Asp (RGD)‐binding site recognition by P11 was site specific because the P11 was inactive for the complex formation of a denatured form of integrin–vitronectin. P11 showed a strong antagonism against αvβ3‐GRGDSP interaction with an IC₅₀ value of 25.72±3.34 nM, whereas the value of GRGDSP peptide was 1968.73±444.32 nM. The binding‐free energies calculated from the docking simulations for each P11 and RGD peptide were ?3.99 and ?3.10 kcal/mol, respectively. The free energy difference between P11 and RGD corresponds to approximately a 4.5‐fold lower K_i value for the P11 than the RGD peptide. The binding orientation of the docked P11 was similar to the crystal structure of the RGD in αvβ3. The analyzed docked poses suggest that a divalent metal–ion coordination was a common driving force for the formation of both SDV/αvβ3 and RGD/αvβ3 complexes. This is the first report on the specific recognition of the RGD‐binding site of αvβ3 by a non‐RGD containing peptide using a computer‐assisted proteomic approach.