Development and validation of a spectral library searching method for peptide identification from MS/MS

A notable inefficiency of shotgun proteomics experiments is the repeated rediscovery of the same identifiable peptides by sequence database searching methods, which often are time-consuming and error-prone. A more precise and efficient method, in which previously observed and identified peptide MS/MS spectra are catalogued and condensed into searchable spectral libraries to allow new identifications by spectral matching, is seen as a promising alternative. To that end, an open-source, functionally complete, high-throughput and readily extensible MS/MS spectral searching tool, SpectraST, was developed. A high-quality spectral library was constructed by combining the high-confidence identifications of millions of spectra taken from various data repositories and searched using four sequence search engines. The resulting library consists of over 30,000 spectra for Saccharomyces cerevisiae. Using this library, SpectraST vastly outperforms the sequence search engine SEQUEST in terms of speed and the ability to discriminate good and bad hits. A unique advantage of SpectraST is its full integration into the popular Trans Proteomic Pipeline suite of software, which facilitates user adoption and provides important functionalities such as peptide and protein probability assignment, quantification, and data visualization. This method of spectral library searching is especially suited for targeted proteomics applications, offering superior performance to traditional sequence searching.     

Selective suppression of excipient signals in 2D <Superscript>1</Superscript>H–<Superscript>13</Superscript>C methyl spectra of biopharmaceutical products

While the use of ¹H–¹³C methyl correlated NMR spectroscopy at natural isotopic abundance has been demonstrated as feasible on protein therapeutics as large as monoclonal antibodies, spectral interference from aliphatic excipients remains a significant obstacle to its widespread application. These signals can cause large baseline artifacts, obscure protein resonances, and cause dynamic range suppression of weak peaks in non-uniform sampling applications, thus hampering both traditional peak-based spectral analyses as well as emerging chemometric methods of analysis. Here we detail modifications to the 2D ¹H–¹³C gradient-selected HSQC experiment that make use of selective pulsing techniques for targeted removal of interfering excipient signals in spectra of the NISTmAb prepared in several different formulations. This approach is demonstrated to selectively reduce interfering excipient signals while still yielding 2D spectra with only modest losses in protein signal. Furthermore, it is shown that spectral modeling based on the SMILE algorithm can be used to simulate and subtract any residual excipient signals and their attendant artifacts from the resulting 2D NMR spectra.     

A semi-empirical approach for predicting unobserved peptide MS/MS spectra from spectral libraries

Spectral library searching is a promising alternative to sequence database searching in peptide identification from MS/MS spectra. The key advantage of spectral library searching is the utilization of more spectral features to improve score discrimination between good and bad matches, and hence sensitivity. However, the coverage of reference spectral library is limited by current experimental and computational methods. We developed a computational approach to expand the coverage of spectral libraries with semi-empirical spectra predicted from perturbing known spectra of similar sequences, such as those with single amino acid substitutions. We hypothesized that the peptide of similar sequences should produce similar fragmentation patterns, at least in most cases. Our results confirm our hypothesis and specify when this approach can be applied. In actual spectral searching of real data sets, the sensitivity advantage of spectral library searching over sequence database searching can be mostly retained even when all real spectra are replaced by semi-empirical ones. We demonstrated the applicability of this approach by detecting several known non-synonymous single-nucleotide polymorphisms in three large human data sets by spectral searching.     

Spectrum-to-spectrum searching using a proteome-wide spectral library

The unambiguous assignment of tandem mass spectra (MS/MS) to peptide sequences remains a key unsolved problem in proteomics. Spectral library search strategies have emerged as a promising alternative for peptide identification, in which MS/MS spectra are directly compared against a reference library of confidently assigned spectra. Two problems relate to library size. First, reference spectral libraries are limited to rediscovery of previously identified peptides and are not applicable to new peptides, because of their incomplete coverage of the human proteome. Second, problems arise when searching a spectral library the size of the entire human proteome. We observed that traditional dot product scoring methods do not scale well with spectral library size, showing reduction in sensitivity when library size is increased. We show that this problem can be addressed by optimizing scoring metrics for spectrum-to-spectrum searches with large spectral libraries. MS/MS spectra for the 1.3 million predicted tryptic peptides in the human proteome are simulated using a kinetic fragmentation model (MassAnalyzer version2.1) to create a proteome-wide simulated spectral library. Searches of the simulated library increase MS/MS assignments by 24% compared with Mascot, when using probabilistic and rank based scoring methods. The proteome-wide coverage of the simulated library leads to 11% increase in unique peptide assignments, compared with parallel searches of a reference spectral library. Further improvement is attained when reference spectra and simulated spectra are combined into a hybrid spectral library, yielding 52% increased MS/MS assignments compared with Mascot searches. Our study demonstrates the advantages of using probabilistic and rank based scores to improve performance of spectrum-to-spectrum search strategies.     

LC-MS characterization and purity assessment of a prototype bispecific antibody

Bispecific IgG asymmetric (heterodimeric) antibodies offer enhanced therapeutic efficacy, but present unique challenges for drug development. These challenges are related to the proper assembly of heavy and light chains. Impurities such as symmetric (homodimeric) antibodies can arise with improper assembly. A new method to assess heterodimer purity of such bispecific antibody products is needed because traditional separation-based purity assays are unable to separate or quantify homodimer impurities. This paper presents a liquid chromatography-mass spectrometry (LC-MS)-based method for evaluating heterodimeric purity of a prototype asymmetric antibody containing two different heavy chains and two identical light chains. The heterodimer and independently expressed homodimeric standards were characterized by two complementary LC-MS techniques: Intact protein mass measurement of deglycosylated antibody and peptide map analyses. Intact protein mass analysis was used to check molecular integrity and composition. LC-MS^E peptide mapping of Lys-C digests was used to verify protein sequences and characterize post-translational modifications, including C-terminal truncation species. Guided by the characterization results, a heterodimer purity assay was demonstrated by intact protein mass analysis of pure deglycosylated heterodimer spiked with each deglycosylated homodimeric standard. The assay was capable of detecting low levels (2%) of spiked homodimers in conjunction with co-eluting half antibodies and multiple mass species present in the homodimer standards and providing relative purity differences between samples. Detection of minor homodimer and half-antibody C-terminal truncation species at levels as low as 0.6% demonstrates the sensitivity of the method. This method is suitable for purity assessment of heterodimer samples during process and purification development of bispecific antibodies, e.g., clone selection.     

Characterization of glycation in an IgG1 by capillary electrophoresis sodium dodecyl sulfate and mass spectrometry

We report a case study of characterization of a non-enzymatically glycated IgG1 using reducing capillary electrophoresis sodium dodecyl sulfate (CE–SDS) and mass spectrometry (MS). Glycation was found to occur nonspecifically at multiple sites in both the light and heavy chains. The glycated light and heavy chains result in wider peaks eluting late in the reducing CE–SDS profile; in particular, the glycated light chain behaved as a shoulder peak detected by either ultraviolet (UV) or laser-induced fluorescence (LIF) signals. The glycated species can be enriched by boronate affinity chromatography. Analyzing the enriched samples by reversed phase high-performance liquid chromatography in line with time-of-flight MS (RP–HPLC–TOF/MS) revealed adducts of +162 and +324 Da to both the light and heavy chains, suggesting the presence of multiple glycation sites. Tryptic peptide mapping and tandem mass sequencing were used to identify two glycation sites on each of the light and heavy chains.     

Fast protein sequencing of monoclonal antibody by real-time digestion on emitter during nanoelectrospray

Growth in the pharmaceutical industry has led to an increasing demand for rapid characterization of therapeutic monoclonal antibodies. The current methods for antibody sequence confirmation (e.g., N-terminal Edman sequencing and traditional peptide mapping methods) are not sufficient; thus, we developed a fast method for sequencing recombinant monoclonal antibodies using a novel digestion-on-emitter technology. Using this method, a monoclonal antibody can be denatured, reduced, digested, and sequenced in less than an hour. High throughput and satisfactory protein sequence coverage were achieved by using a non-specific protease from Aspergillus saitoi, protease XIII, to digest the denatured and reduced monoclonal antibody on an electrospray emitter, while electrospray high voltage was applied to the digestion mixture through the emitter. Tandem mass spectrometry data was acquired over the course of enzyme digestion, generating similar information compared to standard peptide mapping experiments in much less time. We demonstrated that this fast protein sequencing method provided sufficient sequence information for bovine serum albumin and two commercially available monoclonal antibodies, mouse IgG1 MOPC21 and humanized IgG1 NISTmAb. For two monoclonal antibodies, we obtained sequence coverage of 90.5–95.1% for the heavy chains and 98.6–99.1% for the light chains. We found that on-emitter digestion by protease XIII generated peptides of various lengths during the digestion process, which was critical for achieving sufficient sequence coverage. Moreover, we discovered that the enzyme-to-substrate ratio was an important parameter that affects protein sequence coverage. Due to its highly automatable and efficient design, our method offers a major advantage over N-terminal Edman sequencing and traditional peptide mapping methods in the identification of protein sequence, and is capable of meeting an ever-increasing demand for monoclonal antibody sequence confirmation in the biopharmaceutical industry.     

重金属结合肽（蛋白）的结构分析

金属离子与金属结合肽（蛋白）的相互作用与应用研究,一直是生物无机化学的重点和热点,也是分子间相互作用研究领域的难点。本研究利用ClustalX、BLAST等生物信息技术与方法对大量已知的重金属结合肽进行分析与数据挖掘。确定筛选获得的重金属结合肽常富含His,无Cys,无金属结合肽模式序列,进化不保守;部分氨基酸序列结构（如六肽）可在蛋白数据库中找到相似序列。序列特征主要为Zn^2＋相关的转录因子。本研究为重金属结合蛋白-重金属离子的相互作用分析简化为重金属结合肽-重金属离子的结构模拟与分析提供了重要的理论基础和研究手段。     

Vibrational circular dichroism and IR spectral analysis as a test of theoretical conformational modeling for a cyclic hexapeptide

A model cyclohexapeptide, cyclo-(Phe-(D)Pro-Gly-Arg-Gly-Asp) was synthesized and its IR and VCD spectra were used as a test of density functional theory (DFT) level predictions of spectral intensities for a peptide with a nonrepeating but partially constricted conformation. Peptide structure and flexibility was estimated by molecular dynamics (MD) simulations and the spectra were simulated using full quantum mechanical (QM) approaches for the complete peptide and for simplified models with truncated side chains. After simulated annealing, the backbone conformation of the ring structure is relatively stable, consisting of a normal beta-turn and a tight loop (no H-bond) which does not vary over short trajectories. Only in quite long MD runs at high temperatures do other conformations appear. MD simulations were carried out for the cyclic peptide in water and in TFE, which match experimental solvents, as well as with and without protonation of the Asp carboxyl group. DFT spectral simulations were made using the annealed structure and were extended to include basis set variation, to determine an optimal computational approach, and solvent simulation with a polarized continuum model (PCM). Stepwise full DFT simulation of spectra was done for various sequences with the same backbone geometry but based on (1) solely Gly residues, (2) Ala substitution except Gly and Pro, and (3) complete sequences with side chains. Additionally, a selection of structures was used to compute IR and VCD spectra with the optimal method to determine structural variation effects. The side chains, especially the Asp-COOH and Arg-NH(2) transitions, had an impact on the computed amide frequencies, IR intensities and VCD pattern. Since experimentally these groups would have little chirality, due to conformational variation, they do not impact the observed VCD spectra. Correcting for frequency shifts, the Ala model for the cyclopeptide gives the clearest representation of the amide VCD. The experimental sign pattern for the amide I' band in D(2)O and also the sharper, more intense amide I VCD band in TFE was seen to some degree in one conformer with Type II' turns, but the data favor a mix of structures.     

Intrachain disulfide bond in the core hinge region of human IgG4.

IgG is a tetrameric protein composed of two copies each of the light and heavy chains. The four-chain structure is maintained by strong noncovalent interactions between the amino-terminal half of pairs of heavy-light chains and between the carboxyl-terminal regions of the two heavy chains. In addition, interchain disulfide bonds link each heavy-light chain and also link the paired heavy chains. An engineered human IgG4 specific for human tumor necrosis factor-alpha (CDP571) is similar to human myeloma IgG4 in that it is secreted as both disulfide bonded tetramers (approximately 75% of the total amount of IgG) and as tetramers composed of nondisulfide bonded half-IgG4 (heavy chain disulfide bonded to light chain) molecules. However, when CDP571 was genetically engineered with a proline at residue 229 of the core hinge region rather than serine, CDP571 (S229P), or with an IgG1 rather than IgG4 hinge region, CDP571(gamma 1), only trace amounts of nondisulfide bonded half-IgG tetramers were observed. Trypsin digest reversephase HPLC peptide mapping studies of CDP571 and CDP571(gamma 1) with on-line electrospray ionization mass spectroscopy supplemented with Edman sequencing identified the chemical factor preventing inter-heavy chain disulfide bond formation between half-IgG molecules: the two cysteines in the IgG4 and IgG1 core hinge region (CPSCP and CPPCP, respectively) are capable of forming an intrachain disulfide bond. Conformational modeling studies on cyclic disulfide bonded CPSCP and CPPCP peptides yielded energy ranges for the low-energy conformations of 31-33 kcal/mol and 40-42 kcal/mol, respectively. In addition, higher torsion and angle bending energies were observed for the CPPCP peptide due to backbone constraints caused by the extra proline. These modeling results suggest a reason why a larger fraction of intrachain bonds are observed in IgG4 rather than IgG1 molecules: the serine in the core hinge region of IgG4 allows more hinge region flexibility than the proline of IgG1 and thus may permit formation of a stable intrachain disulfide bond more readily.     

Quantitative analysis of glycation and its impact on antigen binding

Glycation has been observed in antibody therapeutics manufactured by the fed-batch fermentation process. It not only increases the heterogeneity of antibodies, but also potentially affects product safety and efficacy. In this study, non-glycated and glycated fractions enriched from a monoclonal antibody (mAb1) as well as glucose-stressed mAb1 were characterized using a variety of biochemical, biophysical and biological assays to determine the effects of glycation on the structure and function of mAb1. Glycation was detected at multiple lysine residues and reduced the antigen binding activity of mAb1. Heavy chain Lys100, which is located in the complementary-determining region of mAb1, had the highest levels of glycation in both stressed and unstressed samples, and glycation of this residue was likely responsible for the loss of antigen binding based on hydrogen/deuterium exchange mass spectrometry analysis. Peptide mapping and intact liquid chromatography-mass spectrometry (LC-MS) can both be used to monitor the glycation levels. Peptide mapping provides site specific glycation results, while intact LC-MS is a quicker and simpler method to quantitate the total glycation levels and is more useful for routine testing. Capillary isoelectric focusing (cIEF) can also be used to monitor glycation because glycation induces an acidic shift in the cIEF profile. As expected, total glycation measured by intact LC-MS correlated very well with the percentage of total acidic peaks or main peak measured by cIEF. In summary, we demonstrated that glycation can affect the function of a representative IgG1 mAb. The analytical characterization, as described here, should be generally applicable for other therapeutic mAbs.     

Rates and impact of human antibody glycation in vivo

Glycation of immunoglobulin G (IgG) can result from incubation with a reducing sugar in vitro or during circulation in vivo. Upon injection of a recombinantly produced human therapeutic IgG into humans, changes in the glycation levels could be observed as a function of circulation time. Mass changes on the individual IgG polypeptide chains as the results of glycation were determined using reversed-phase liquid chromatography/mass spectrometry. Changes to the light and heavy chains were low but easily detectable at 0.00092 and 0.0021 glucose (Glc) additions per chain per day, respectively. Levels of glycation found on the Fc portion of IgG isolated from healthy subjects, using a similar analytical approach, were on average 0.045 Glc molecules per fragment. In vivo glycation rates could be approximated in vitro by modeling the physiological glycation reaction with a simplified incubation containing physiological Glc concentrations, pH and temperature but with a high concentration of a single purified IgG. To test the impact of glycation on IgG function, highly glycated IgG1 and IgG2 were prepared containing on average 42-49 Glc molecules per IgG. Binding to FcγIIIa receptors, neonatal Fc receptor or protein A was similar or identical to the non-glycated IgG controls. Although the modifications were well distributed throughout the protein sequence, and at high enough levels to affect the elution position by size-exclusion chromatography, no changes in the tested Fc functions were observed.     

Validation of peptide MS/MS spectra using metabolic isotope labeling for spectral matching-based shotgun proteome analysis

We report an isotope labeling shotgun proteome analysis strategy to validate the spectrum-to-sequence assignments generated by using sequence-database searching for the construction of a more reliable MS/MS spectral library. This strategy is demonstrated in the analysis of the E. coli K12 proteome. In the workflow, E. coli cells were cultured in normal and (15)N-enriched media. The differentially labeled proteins from the cell extracts were subjected to trypsin digestion and two-dimensional liquid chromatography quadrupole time-of-flight tandem mass spectrometry (2D-LC QTOF MS/MS) analysis. The MS/MS spectra of the two samples were individually searched using Mascot against the E. coli proteome database to generate lists of peptide sequence matches. The two data sets were compared by overlaying the spectra of unlabeled and labeled matches of the same peptide sequence for validation. Two cutoff filters, one based on the number of common fragment ions and another one on the similarity of intensity patterns among the common ions, were developed and applied to the overlaid spectral pairs to reject the low quality or incorrectly assigned spectra. By examining 257,907 and 245,156 spectra acquired from the unlabeled and (15)N-labeled samples, respectively, an experimentally validated MS/MS spectral library of tryptic peptides was constructed for E. coli K12 that consisted of 9,302 unique spectra with unique sequence and charge state, representing 7,763 unique peptide sequences. This E. coli spectral library could be readily expanded, and the overall strategy should be applicable to other organisms. Even with this relatively small library, it was shown that more peptides could be identified with higher confidence using the spectral search method than by sequence-database searching.     

Quantitative Proteomics Using Reductive Dimethylation for Stable Isotope Labeling

Stable isotope labeling of peptides by reductive dimethylation (ReDi labeling) is a method to accurately quantify protein expression differences between samples using mass spectrometry. ReDi labeling is performed using either regular (light) or deuterated (heavy) forms of formaldehyde and sodium cyanoborohydride to add two methyl groups to each free amine. Here we demonstrate a robust protocol for ReDi labeling and quantitative comparison of complex protein mixtures. Protein samples for comparison are digested into peptides, labeled to carry either light or heavy methyl tags, mixed, and co-analyzed by LC-MS/MS. Relative protein abundances are quantified by comparing the ion chromatogram peak areas of heavy and light labeled versions of the constituent peptide extracted from the full MS spectra. The method described here includes sample preparation by reversed-phase solid phase extraction, on-column ReDi labeling of peptides, peptide fractionation by basic pH reversed-phase (BPRP) chromatography, and StageTip peptide purification. We discuss advantages and limitations of ReDi labeling with respect to other methods for stable isotope incorporation. We highlight novel applications using ReDi labeling as a fast, inexpensive, and accurate method to compare protein abundances in nearly any type of sample.     

Building and searching tandem mass (MS/MS) spectral libraries for peptide identification in proteomics

Spectral library searching is an emerging approach in peptide identifications from tandem mass spectra, a critical step in proteomic data analysis. In spectral library searching, a spectral library is first meticulously compiled from a large collection of previously observed peptide MS/MS spectra that are conclusively assigned to their corresponding amino acid sequence. An unknown spectrum is then identified by comparing it to all the candidates in the spectral library for the most similar match. This review discusses the basic principles of spectral library building and searching, describes its advantages and limitations, and provides a primer for researchers interested in adopting this new approach in their data analysis. It will also discuss the future outlook on the evolution and utility of spectral libraries in the field of proteomics.     

Improving SRM assay development: a global comparison between triple quadrupole, ion trap, and higher energy CID peptide fragmentation spectra

In proteomics, selected reaction monitoring (SRM) is rapidly gaining importance for targeted protein quantification. The triple quadrupole mass analyzers used in SRM assays allow for levels of specificity and sensitivity hard to accomplish by more standard shotgun proteomics experiments. Often, an SRM assay is built by in silico prediction of transitions and/or extraction of peptide precursor and fragment ions from a spectral library. Spectral libraries are typically generated from nonideal ion trap based shotgun proteomics experiments or synthetic peptide libraries, consuming considerable time and effort. Here, we investigate the usability of beam type CID (or "higher energy CID" (HCD)) peptide fragmentation spectra, as acquired using an Orbitrap Velos, to facilitate SRM assay development. Therefore, peptide fragmentation spectra, obtained by ion-trap CID, triple-quadrupole CID (QqQ-CID) and Orbitrap HCD, originating from digested cellular lysates, were compared. Spectral comparison and a dedicated correlation algorithm indicated significantly higher similarity between QqQ-CID and HCD fragmentation spectra than between QqQ-CID and ion trap-CID spectra. SRM transitions generated using a constructed HCD spectral library increased SRM assay sensitivity up to 2-fold, when compared to the use of a library created from more conventionally used ion trap-CID spectra, showing that HCD spectra can assist SRM assay development.     

A chemical and computational approach to comprehensive glycation characterization on antibodies

Non-enzymatic glycation is a challenging post-translational modification to characterize due to the structural heterogeneity it generates in proteins. Glycation has become increasingly recognized as an important product quality attribute to monitor, particularly for the biotechnology sector, which produces recombinant proteins under conditions that are amenable to protein glycation. The elucidation of sites of glycation can be problematic using conventional collision-induced dissociation (CID)-based mass spectrometry because of the predominance of neutral loss ions. A method to characterize glycation using an IgG1 monoclonal antibody (mAb) as a model is reported here. The sugars present on this mAb were derivatized using sodium borohydride chemistry to stabilize the linkage and identified using CID-based MS² mass spectrometry and spectral search engines. Quantification of specific glycation sites was then done using a targeted MS¹ based approach, which allowed the identification of a glycation hot spot in the heavy chain complementarity-determining region 3 of the mAb. This targeted approach provided a path forward to developing a structural understanding of the propensity of sites to become glycated on mAbs. Through structural analysis we propose a model in which the number and 3-dimensional distances of carboxylic acid amino acyl residues create a favorable environment for glycation to occur.     

A mass spectrometry-based method to screen for α-amidated peptides

Amidation is a post-translational modification found at the C-terminus of ~50% of all neuropeptide hormones. Cleavage of the C(α)-N bond of a C-terminal glycine yields the α-amidated peptide in a reaction catalyzed by peptidylglycine α-amidating monooxygenase (PAM). The mass of an α-amidated peptide decreases by 58 Da relative to its precursor. The amino acid sequences of an α-amidated peptide and its precursor differ only by the C-terminal glycine meaning that the peptides exhibit similar RP-HPLC properties and tandem mass spectral (MS/MS) fragmentation patterns. Growth of cultured cells in the presence of a PAM inhibitor ensured the coexistence of α-amidated peptides and their precursors. A strategy was developed for precursor and α-amidated peptide pairing (PAPP): LC-MS/MS data of peptide extracts were scanned for peptide pairs that differed by 58 Da in mass, but had similar RP-HPLC retention times. The resulting peptide pairs were validated by checking for similar fragmentation patterns in their MS/MS data prior to identification by database searching or manual interpretation. This approach significantly reduced the number of spectra requiring interpretation, decreasing the computing time required for database searching and enabling manual interpretation of unidentified spectra. Reported here are the α-amidated peptides identified from AtT-20 cells using the PAPP method.     

The Generation of a Comprehensive Spectral Library for the Analysis of the Guinea Pig Proteome by SWATH‐MS

Advances in liquid chromatography‐mass spectrometry have facilitated the incorporation of proteomic studies to many biology experimental workflows. Data‐independent acquisition platforms, such as sequential window acquisition of all theoretical mass spectra (SWATH‐MS), offer several advantages for label‐free quantitative assessment of complex proteomes over data‐dependent acquisition (DDA) approaches. However, SWATH data interpretation requires spectral libraries as a detailed reference resource. The guinea pig (Cavia porcellus) is an excellent experimental model for translation to many aspects of human physiology and disease, yet there is limited experimental information regarding its proteome. To overcome this knowledge gap, a comprehensive spectral library of the guinea pig proteome is generated. Homogenates and tryptic digests are prepared from 16 tissues and subjected to >200 DDA runs. Analysis of >250 000 peptide‐spectrum matches resulted in a library of 73 594 peptides from 7666 proteins. Library validation is provided by i) analyzing externally derived SWATH files ( https://doi.org/10.1016/j.jprot.2018.03.023 ) and comparing peptide intensity quantifications; ii) merging of externally derived data to the base library. This furnishes the research community with a comprehensive proteomic resource that will facilitate future molecular‐phenotypic studies using (re‐engaging) the guinea pig as an experimental model of relevance to human biology. The spectral library and raw data are freely accessible in the MassIVE repository (MSV000083199).