Virtual-Experimental 2DE Approach in Chromosome-Centric Human Proteome Project

Puzzle of Proteoform Variety-Where Is a Key?

One of the human proteome puzzles is an imbalance between the theoretically calculated and experimentally measured amounts of proteoforms. Considering the possibility of combinations of different post-translational modifications (PTMs), the quantity of possible proteoforms is huge. An estimation gives more than a million different proteoforms in each cell type. But, it seems that there is strict control over the production and maintenance of PTMs. Although the potential complexity of proteoforms due to PTMs is tremendous, available information indicates that only a small part of it is being implemented. As a result, a protein could have many proteoforms according to the number of modification sites, but because of different systems of personal regulation, the profile of PTMs for a given protein in each organism is slightly different.

Variability of haptoglobin beta-chain proteoforms

Existing knowledge on changes of the haptoglobin (Hp) molecule suggests that it may exist in multiple proteoforms, which obviously exhibit different functions. Using two-dimensional electrophoresis (2DE) in combination with mass spectrometry and immunodetection, we have analyzed blood plasma samples from both healthy donors and patients with primary grade IV glioblastoma (GBM), and obtained a detailed composite 2DE distribution map of β-chain proteoforms, as well as the full-length form of Hp (zonulin). Although the total level of plasma Hp exceeded normal values in cancer patients (especially patients with GBM), the presence of particuar proteoforms, detected by their position on the 2DE map, was very individual. Variability was found in both zonulin and the Hp β-chain. The presence of an alkaline form of zonulin in plasma can be considered a conditional, but insufficient, GBM biomarker. In other words, we found that at the level of minor proteoforms of Hp, even in normal conditions, there was a high individual variability. On the one hand, this raises questions about the reasons for such variability, if it is present not only in Hp, but also in other proteins. On the other hand, this may explain the discrepancy between the number of experimentally detected proteoforms and the theoretically possible ones not only in Hp, but also in other proteins.

Quantitative Aspects of the Human Cell Proteome

The number and identity of proteins and proteoforms presented in a single human cell (a cellular proteome) are fundamental biological questions. The answers can be found with sophisticated and sensitive proteomics methods, including advanced mass spectrometry (MS) coupled with separation by gel electrophoresis and chromatography. So far, bioinformatics and experimental approaches have been applied to quantitate the complexity of the human proteome. This review analyzed the quantitative information obtained from several large-scale panoramic experiments in which high-resolution mass spectrometry-based proteomics in combination with liquid chromatography or two-dimensional gel electrophoresis (2DE) were used to evaluate the cellular proteome. It is important that even though all these experiments were performed in different labs using different equipment and calculation algorithms, the main conclusion about the distribution of proteome components (proteins or proteoforms) was basically the same for all human tissues or cells. It follows Zipf's law and has a formula N = A/x, where N is the number of proteoforms, A is a coefficient, and x is the limit of proteoform detection in terms of abundance.

Proteomic Markers in the Muscles and Brain of Pigs Recovered from Hemorrhagic Stroke

(1) Background: Stroke is the leading cause of serious long-term disability. Walking dysfunction and paresis of the upper extremities occurs in more than 80% of people who have had a stroke. (2) Methods: We studied post-genomic markers in biosamples of muscle and brain tissue from animals that underwent intracerebral hematoma and recovered after 42 days. Our purpose was to understand the biological mechanisms associated with recovery from hemorrhagic stroke. We analyzed the peptides formed after trypsinolysis of samples by HPLC-MS, and the results were processed by bioinformatics methods, including the establishment of biochemical relationships (gene to gene) using topological omics databases such as Reactome and KEGG. (3) Results: In the pig brain, unique compounds were identified which are expressed during the recovery period after traumatic injury. These are molecular factors of activated microglia, and they contribute to the functional recovery of neurons and reduce instances of hematoma, edema, and oxidative stress. Complexes of the main binding factors of the neurotrophins involved in the differentiation and survival of nerve cells were found in muscles. (4) Conclusions: A network of gene interactions has been constructed for proteins involved in the regulation of synaptic transmission, in particular presynaptic vesicular and endocytic processes. The presence of transmitters and transporters associated with stimulation of NMDA receptors at neuromuscular junctions shows the relationship between upper motor neurons and neuromuscular junctions.

Construction of 2DE Patterns of Plasma Proteins: Aspect of Potential Tumor Markers

The use of tumor markers aids in the early detection of cancer recurrence and prognosis. There is a hope that they might also be useful in screening tests for the early detection of cancer. Here, the question of finding ideal tumor markers, which should be sensitive, specific, and reliable, is an acute issue. Human plasma is one of the most popular samples as it is commonly collected in the clinic and provides noninvasive, rapid analysis for any type of disease including cancer. Many efforts have been applied in searching for “ideal” tumor markers, digging very deep into plasma proteomes. The situation in this area can be improved in two ways—by attempting to find an ideal single tumor marker or by generating panels of different markers. In both cases, proteomics certainly plays a major role. There is a line of evidence that the most abundant, so-called “classical plasma proteins”, may be used to generate a tumor biomarker profile. To be comprehensive these profiles should have information not only about protein levels but also proteoform distribution for each protein. Initially, the profile of these proteins in norm should be generated. In our work, we collected bibliographic information about the connection of cancers with levels of “classical plasma proteins”. Additionally, we presented the proteoform profiles (2DE patterns) of these proteins in norm generated by two-dimensional electrophoresis with mass spectrometry and immunodetection. As a next step, similar profiles representing protein perturbations in plasma produced in the case of different cancers will be generated. Additionally, based on this information, different test systems can be developed.

Evaluation of Haptoglobin and Its Proteoforms as Glioblastoma Markers

Haptoglobin (Hp) is a blood plasma glycoprotein that plays a critical role in tissue protection and the prevention of oxidative damage. Haptoglobin is an acute-phase protein, its concentration in plasma changes in pathology, and the test for its concentration is part of normal clinical practice. Haptoglobin is a conservative protein and is the subject of research as a potential biomarker of many diseases, including malignant neoplasms. The Human Hp gene is polymorphic and controls the synthesis of three major phenotypes-homozygous Hp1-1 and Hp2-2, and heterozygous Hp2-1, determined by a combination of allelic variants that are inherited. Numerous studies indicate that the phenotype of haptoglobin can be used to judge the individual's predisposition to various diseases. In addition, Hp undergoes various post-translational modifications (PTMs). Glioblastoma multiform (GBM) is the most malignant primary brain tumor. In our study, we have analyzed the state of Hp proteoforms in plasma and cells using 1D (SDS-PAGE) and 2D electrophoresis (2DE) with the following mass spectrometry (LC ES-MS/MS) or Western blotting. We found that the levels of α2- and β-chain proteoforms are up-regulated in the plasma of GBM patients. An unprocessed form of Hp2-2 (PreHp2-2, zonulin) with unusual biophysical parameters (pI/Mw) was also detected in the plasma of GBM patients and glioblastoma cells. Altogether, this data shows the possibility to use proteoforms of haptoglobin as a potential GBM-specific plasma biomarker.

Empowering Shotgun Mass Spectrometry with 2DE: A HepG2 Study

One of the major goals of the Chromosome-Centric Human Proteome Project (C-HPP) is to catalog and annotate a myriad of heterogeneous proteoforms, produced by ca. 20 thousand genes. To achieve a detailed and personalized understanding into proteomes, we suggest using a customized RNA-seq library of potential proteoforms, which includes aberrant variants specific to certain biological samples. Two-dimensional electrophoresis coupled with high-performance liquid chromatography allowed us to downgrade the difficulty of biological mixing following shotgun mass spectrometry. To benchmark the proposed pipeline, we examined heterogeneity of the HepG2 hepatoblastoma cell line proteome. Data are available via ProteomeXchange with identifier PXD018450.

A database for inventory of proteoform profiles: "2DE-pattern"

The human proteome is composed of a diverse and heterogeneous range of gene products/proteoforms/protein species. Because of the growing amount of information about proteoforms generated by different methods, we need a convenient approach to make an inventory of the data. Here, we present a database of proteoforms that is based on information obtained by separation of proteoforms using 2DE followed by shotgun ESI-LC-MS/MS. The database's principles and structure are described. The database is called "2DE-pattern" as it contains multiple isoform-centric patterns of proteoforms separated according to 2DE principles. The database can be freely used at http://2de-pattern.pnpi.nrcki.ru.

Innovating the Concept and Practice of Two-Dimensional Gel Electrophoresis in the Analysis of Proteomes at the Proteoform Level

Two-dimensional gel electrophoresis (2DE) is an important and well-established technical platform enabling extensive top-down proteomic analysis. However, the long-held but now largely outdated conventional concepts of 2DE have clearly impacted its application to in-depth investigations of proteomes at the level of protein species/proteoforms. It is time to popularize a new concept of 2DE for proteomics. With the development and enrichment of the proteome concept, any given “protein” is now recognized to consist of a series of proteoforms. Thus, it is the proteoform, rather than the canonical protein, that is the basic unit of a proteome, and each proteoform has a specific isoelectric point (pI) and relative mass (M_r). Accordingly, using 2DE, each proteoform can routinely be resolved and arrayed according to its different pI and M_r. Each detectable spot contains multiple proteoforms derived from the same gene, as well as from different genes. Proteoforms derived from the same gene are distributed into different spots in a 2DE pattern. High-resolution 2DE is thus actually an initial level of separation to address proteome complexity and is effectively a pre-fractionation method prior to analysis using mass spectrometry (MS). Furthermore, stable isotope-labeled 2DE coupled with high-sensitivity liquid chromatography-tandem MS (LC-MS/MS) has tremendous potential for the large-scale detection, identification, and quantification of the proteoforms that constitute proteomes.

Comprehensive Analysis of Carbohydrate-Active Enzymes from the Filamentous Fungus Scytalidium candidum 3C

Complete enzymatic degradation of plant polysaccharides is a result of combined action of various carbohydrate-active enzymes (CAZymes). In this paper, we demonstrate the potential of the filamentous fungus Scytalidium candidum 3C for processing of plant biomass. Structural annotation of the improved assembly of S. candidum 3C genome and functional annotation of CAZymes revealed putative gene sequences encoding such proteins. A total of 190 CAZyme-encoding genes were identified, including 104 glycoside hydrolases, 52 glycosyltransferases, 28 oxidative enzymes, and 6 carbohydrate esterases. In addition, 14 carbohydrate-binding modules were found. Glycoside hydrolases secreted during the growth of S. candidum 3C in three media were analyzed with a variety of substrates. Mass spectrometry analysis of the fungal culture liquid revealed the presence of peptides identical to 36 glycoside hydrolases, three proteins without known enzymatic function belonging to the same group of families, and 11 oxidative enzymes. The activity of endo-hemicellulases was determined using specially synthesized substrates in which the glycosidic bond between monosaccharide residues was replaced by a thio-linkage. During analysis of the CAZyme profile of S. candidum 3C, four β-xylanases from the GH10 family and two β-glucanases from the GH7 and GH55 families were detected, partially purified, and identified.

Next Steps on in Silico 2DE Analyses of Chromosome 18 Proteoforms

In the boundaries of the chromosome-centric Human Proteome Project (c-HPP) to obtain information about proteoforms coded by chromosome 18, several cell lines (HepG2, glioblastoma, LEH), normal liver, and plasma were analyzed. In our study, we have been using proteoform separation by two-dimensional electrophoresis (2DE) (a sectional analysis) and a semivirtual 2DE with following shotgun mass spectrometry using LC-ESI-MS/MS. Previously, we published a first draft of this research, where only HepG2 cells were tested. Here, we present the next step using more detailed analysis and more samples. Altogether, confident (2 significant sequences minimum) information about proteoforms of 117 isoforms coded by 104 genes of chromosome 18 was obtained. The 3D-graphs showing distribution of different proteoforms from the same gene in the 2D map were generated. Additionally, a semivirtual 2DE approach has allowed for detecting more proteoforms and estimating their pI more precisely. Data are available via ProteomeXchange with identifier PXD010142.

Coagulation cascade and complement system in systemic lupus erythematosus

This study was conducted to (1) characterize coagulation cascade and complement system in systemic lupus erythematosus (SLE); (2) evaluate the associations between coagulation cascade, complement system, inflammatory response and SLE disease severity; (3) test the diagnostic value of a combination of D-dimer and C4 for lupus activity. Transcriptomics, proteomics and metabolomics were performed in 24 SLE patients and 24 healthy controls. The levels of ten coagulations, seven complements and three cytokines were measured in 112 SLE patients. Clinical data were collected from 2025 SLE patients. The analysis of multi-omics data revealed the common links for the components of coagulation cascade and complement system. The results of ELISA showed coagulation cascade and complement system had an interaction effect on SLE disease severity, this effect was pronounced among patients with excess inflammation. The analysis of clinical data revealed a combination of D-dimer and C4 provided good diagnostic performance for lupus activity. This study suggested that coagulation cascade and complement system become 'partners in crime', contributing to SLE disease severity and identified the diagnostic value of D-dimer combined with C4for lupus activity.

Variety and Dynamics of Proteoforms in the Human Proteome: Aspects of Markers for Hepatocellular Carcinoma

We have previously developed an approach, where two-dimensional gel electrophoresis (2DE) was used, followed by sectional analysis of the whole gel using high-resolution nano-liquid chromatography-mass spectrometry (ESI LC-MS/MS). In this study, we applied this approach on the panoramic analysis of proteins and their proteoforms from normal (liver) and cancer (HepG2) cells. This allowed us to detect, in a single proteome, about 20,000 proteoforms coded by more than 4000 genes. A set of 3D-graphs showing distribution of these proteoforms in 2DE maps (profiles) was generated. A comparative analysis of these profiles between normal and cancer cells showed high variability and dynamics of many proteins. Among these proteins, there are some well-known features like alpha-fetoprotein (FETA) or glypican-3 (GPC3) and potential hepatocellular carcinoma (HCC) markers. More detailed information about their proteoforms could be used for generation of panels of more specific biomarkers.

Proteomic Cinderella: Customized analysis of bulky MS/MS data in one night

Proteomic challenges, stirred up by the advent of high-throughput technologies, produce large amount of MS data. Nowadays, the routine manual search does not satisfy the "speed" of modern science any longer. In our work, the necessity of single-thread analysis of bulky data emerged during interpretation of HepG2 proteome profiling results for proteoforms searching. We compared the contribution of each of the eight search engines (X!Tandem, MS-GF[Formula: see text], MS Amanda, MyriMatch, Comet, Tide, Andromeda, and OMSSA) integrated in an open-source graphical user interface SearchGUI ( http://searchgui.googlecode.com ) into total result of proteoforms identification and optimized set of engines working simultaneously. We also compared the results of our search combination with Mascot results using protein kit UPS2, containing 48 human proteins. We selected combination of X!Tandem, MS-GF[Formula: see text] and OMMSA as the most time-efficient and productive combination of search. We added homemade java-script to automatize pipeline from file picking to report generation. These settings resulted in rise of the efficiency of our customized pipeline unobtainable by manual scouting: the analysis of 192 files searched against human proteome (42153 entries) downloaded from UniProt took 11[Formula: see text]h.

[Multiomics study of HepG2 cell line proteome]

Current proteomic studies are generally focused on the most abundant proteoforms encoded by canonical nucleic sequences. Transcriptomic and proteomic data, accumulated in a variety of postgenome sources and coupled with state-of-art analytical technologies, allow to start the identification of aberrant (non-canonical) proteoforms. The main sources of aberrant proteoforms are alternative splicing, single nucleotide polymorphism, and post-translational modifications. The aim of this work was to estimate the heterogeneity of HepG2 proteome. We suggested multiomics approach, which combines transcriptomic (RNAseq) and proteomic (2DE-MS/MS) methods, as a promising strategy to explore the proteome.

A semi-virtual two dimensional gel electrophoresis: IF-ESI LC-MS/MS

A method for increasing the productivity of ESI LC-MS/MS (electrospray ionization-liquid chromatography-tandem mass spectrometry) was proposed and applied.

After IF (isoelectric focusing) of the sample using IPG (immobilized pH gradient) strip, the strip was cut to sections, and every section was treated according to trypsinolysis protocol for MS/MS analysis. The peptides produced were further analyzed by ESI LC-MS/MS. The procedure allows to:

•

identify many more proteins and proteoforms compared to shotgun analysis of extracts.

•

build a semi-virtual 2DE map of identified proteins.

Dataset of protein species from human liver

This article contains data related to the research article entitled "Zipf׳s law in proteomics" (Naryzhny et al., 2017) [1]. The protein composition in the human liver or hepatocarcinoma (HepG2) cells extracts was estimated using a filter-aided sample preparation (FASP) protocol. The protein species/proteoform composition in the human liver was determined by two-dimensional electrophoresis (2-DE) followed by Electrospray Ionization Liquid Chromatography-Tandem Mass Spectrometry (ESI LC-MS/MS). In the case of two-dimensional electrophoresis (2-DE), the gel was stained with Coomassie Brilliant Blue R350, and image analysis was performed with ImageMaster 2D Platinum software (GE Healthcare). The 96 sections in the 2D gel were selected and cut for subsequent ESI LC-MS/MS and protein identification. If the same protein was detected in different sections, it was considered to exist as different protein species/proteoforms. A list of human liver proteoforms detected in this way is presented.

Towards the Full Realization of 2DE Power

Here, approaches that allow disclosure of the information hidden inside and outside of two-dimensional gel electrophoresis (2DE) are described. Experimental identification methods, such as mass spectrometry of high resolution and sensitivity (MALDI-TOF MS and ESI LC-MS/MS) and immunodetection (Western and Far-Western) in combination with bioinformatics (collection of all information about proteoforms), move 2DE to the next level of power. The integration of these technologies will promote 2DE as a powerful methodology of proteomics technology.