EpiProfile Quantifies Histone Peptides With Modifications by Extracting Retention Time and Intensity in High-resolution Mass Spectra

Histone post-translational modifications contribute to chromatin function through their chemical properties which influence chromatin structure and their ability to recruit chromatin interacting proteins. Nanoflow liquid chromatography coupled with high resolution tandem mass spectrometry (nanoLC-MS/MS) has emerged as the most suitable technology for global histone modification analysis because of the high sensitivity and the high mass accuracy of this approach that provides confident identification. However, analysis of histones with this method is even more challenging because of the large number and variety of isobaric histone peptides and the high dynamic range of histone peptide abundances. Here, we introduce EpiProfile, a software tool that discriminates isobaric histone peptides using the distinguishing fragment ions in their tandem mass spectra and extracts the chromatographic area under the curve using previous knowledge about peptide retention time. The accuracy of EpiProfile was evaluated by analysis of mixtures containing different ratios of synthetic histone peptides. In addition to label-free quantification of histone peptides, EpiProfile is flexible and can quantify different types of isotopically labeled histone peptides. EpiProfile is unique in generating layouts (i.e. relative retention time) of histone peptides when compared with manual quantification of the data and other programs (such as Skyline), filling the need of an automatic and freely available tool to quantify labeled and non-labeled modified histone peptides. In summary, EpiProfile is a valuable nanoflow liquid chromatography coupled with high resolution tandem mass spectrometry-based quantification tool for histone peptides, which can also be adapted to analyze nonhistone protein samples.The nucleosome, the basic unit of chromatin, consists of 147 base pairs of DNA wrapped around histone proteins (H2A, H2B, H3, and H4). Histones play vital roles in chromatin, interacting with many signaling proteins and chromatin-structural proteins through various post-translational modifications (PTMs)¹ (1–3). There are numerous PTMs on histones, including methylation (mono - me1, di - me2, tri - me3), acetylation (ac), phosphorylation (ph), ubiquitination, and SUMOylation (4). Histone PTMs can affect chromatin function, and therefore influence processes such as gene accessibility, DNA repair and chromosome condensation. Moreover, histone PTMs cross-talk in a synergistic manner to fine-tune gene expression (5). Therefore, quantification of histone PTMs has become a high priority to investigate cell regulation and epigenetics (6).Traditionally, antibody-based methods (e.g. Western blot) have been used to analyze histone modifications (7), which have multiple disadvantages. First, antibodies are not available for every new PTM discovered. Second, PTMs on neighboring amino acids (e.g. H3K9me1–3 and H3S10ph) may prevent antibody binding, a phenomenon called epitope occlusion. Third, the quantification of PTMs via antibody-based methods is not sensitive to small differences (e.g. <twofold). Mass spectrometry (MS) has emerged as a sensitive and efficient technique to detect known and novel PTMs (8). The high mass accuracy and the high speed of modern mass spectrometers allow for sensitive, confident, and accurate peptide quantification when coupled with nanoflow liquid chromatography (nanoLC).NanoLC-MS/MS analysis of protein digests (i.e. bottom-up MS) is nowadays a mature and widely applied technology. Data-dependent acquisition is the most commonly adopted MS acquisition method to identify peptides via bottom-up MS (9–12), generating MS1 and MS2 spectra. Nevertheless, histone proteins are particularly challenging to analyze by using the generalized bottom-up workflow. As histones are rich with lysines and arginines, tryptic digest of histones generates short peptides that are difficult to be retained on C18 columns. To improve histone peptide retention, the unmodified and mono-methylated lysines and peptide N terminus can be selectively chemically propionylated (13–16), preventing tryptic digest after lysine to generate longer peptides. Moreover, peptide identification through traditional database searches leads to a large number of false positives, as allowing several dynamic modifications (e.g. me1/me2/me3, ac, ph) dramatically increases the number of molecular candidates and thus the possibility to achieve a false hit (12). Therefore, software tools that quantify histone peptides require additional data to correctly map a given peptide, such as previous knowledge of peptide retention time.Quantification of histone peptides is particularly challenging because of presence of isobaric peptides, near isobaric PTMs such as tri-methylation (42.047 Da) and acetylation (42.011 Da), and low abundant species. Previous knowledge about relative peptide retention time (RT) enables differentiation between species close in mass and therefore selection of the correct peak for integration of the area of the chromatographic peak (i.e. area under curve or AUC). However, determination of peptide RT might be difficult because of their low abundance though acid extraction was performed to purify histones. This problem can be solved by using isotopically labeled synthetic histone peptides (17), or data independent approaches (18). When using relative retention time information to assign peak identities, reproducible nanoLC is crucial, especially because some isobaric peptides co-elute. In this case, the MS acquisition method must perform targeted MS2 for the co-eluting isobaric peptides at the specific time that they elute. These species can be discriminated and quantified based on the intensity of fragment ions unique to each species. For instance, the peptides KacSTGGKAPR (H3K9ac) and KSTGGKacAPR (H3K14ac) have the same mass and overlap at the nanoLC elution (the full protein sequence of human canonical histone H3 and H4 are shown in Fig. 1A). Thus, the co-eluting isobaric peptides could not be quantified separately based on the MS1 signal, but the unique fragment ions present in MS2 spectra allow them to be quantified individually.Open in a separate windowFig. 1.Histones are a challenge for quantitative mass spectrometry analyses. A, Human histone H3.1 and H4 protein sequences. B, Spline fitting to calculate AUC: blue lines are the original peaks and pink lines are the fitted peaks. C, An example of isobaric PTM modified peptides. The above MS2 is matched with H3K18ac, and the same MS2 is also matched with H3K23ac below. D, The workflow of EpiProfile: inputting precursor m/z and charge state, extracting elution profiles, selecting the correct chromatographic peak, calculating AUC, and outputting quantification tables and figures.There have been few computational investigations attempting to solve the problem of quantifying co-eluting isobaric peptides. DiMaggio et al. used a mixed integer linear optimization (MILP) framework to quantify partially co-eluting isobaric histone peptides from electron transfer dissociation (ETD) spectra (19). The framework is comprised of two MILP models: (1) enumerating the entire space of the modified forms that satisfy a given peptide mass and (2) determining the relative composition of the modified forms in the spectrum. Another study by Guan et al. identified isobaric peptides by searching ETD MS/MS spectra for ions representing all possible configurations of modified peptides using a visual assistance program. The relative abundances of these species were estimated by using a nonnegative least squares procedure (20). Other quantification programs can also perform accurate peak picking, but are commonly not as suitable for heavily modified and isobaric histone peptides (e.g. Skyline) (21). These software programs are unable to provide the layouts of histone peptides (i.e. relative RTs) or discriminate all isobaric modified peptides, two tasks that are vital for full characterization of a histone sample.In this study, we developed a new quantification program named EpiProfile. EpiProfile extracts ion chromatography for known histone peptides by using previous knowledge about their elution profiles. Moreover, it discriminates and quantifies the isobaric histone peptides by resolving the linear equations listed with the peak heights of unique fragment ions between the two modification sites in the MS2 spectra (e.g. ions between H3K9ac and H3K14ac). We evaluated the accuracy of EpiProfile by mixing different ratios of synthetic histone peptides, and then tested EpiProfile by analyzing nanoLC-MS/MS data sets of the following samples: purified histones from HeLa cells, a synthetic histone peptide library, and histone peptides labeled during cell growth with ¹³C-labeled glucose media or stable isotope labeling by amino acids in cell culture (SILAC) (22). We compared EpiProfile to manual quantification of the data, and also with the openly available program Skyline. We found that manual quantification is obviously time-consuming and that Skyline cannot generate the layouts of histone peptides and cannot discriminate four or six-component isobaric peptides, a common occurrence in histone data. Moreover, EpiProfile is highly flexible, and thus it can be used to analyze various protein samples, including isotopically labeled peptides and nonhistone data sets.