Key Takeaways
- Fragment ion calculators: Predict theoretical b-ion and y-ion masses from a peptide sequence for CID/HCD, or c/z-ion masses for ETD, enabling database-based peptide identification in mass spectrometry.
- Isoelectric point calculators: Estimate the pH at which a peptide carries zero net charge, using pKa values assigned to ionizable residues and termini; the result predicts gel migration position or off-gel fraction assignment.
- pKa scale dependence: Different pI calculators use different pKa scales (IPC, EMBOSS, ExPASy); the IPC scale has the best empirical validation for 2-DE migration prediction as of 2021.
- Limitations: Both calculators assume linear, unmodified peptides with standard amino acids; cyclic peptides, glycopeptides, and derivatized termini require extended calculation approaches.
- Pharmaceutical context: Fragment ion analysis and pI prediction are commonly part of structural characterization workflows for pharmaceutical-grade peptide active ingredients.
Peptide identification in mass spectrometry rests on a straightforward calculation: if you know the amino acid sequence of a peptide, you can predict the masses of the fragment ions it will produce when fragmented. Match those predictions against the observed spectrum, and you have an identification. The same principle underlies the major proteomics database search engines in common use today, including Mascot, SEQUEST, and OMSSA. The isoelectric point calculation operates on a different but equally foundational principle: given the ionizable residues in a sequence, predict the pH at which the peptide carries no net charge. Both computations are deterministic arithmetic — the value of understanding them lies in knowing when the arithmetic applies and when it does not.
The Fragment Ion Series Formula
Fragment ion calculation is based on the peptide backbone cleavage pattern generated by the fragmentation method used.
For CID and HCD (b-ions and y-ions):
b-ion mass = sum of N-terminal residue masses + 1 (proton)
y-ion mass = sum of C-terminal residue masses + 18.011 (water) + 1 (proton)
For ETD (c-ions and z-ions):
c-ion mass = b-ion mass + 17.027 (NH₃)
z-ion mass = y-ion mass − 16.019 (NH₂)
These formulas generate a series of masses for every possible cleavage position along the peptide backbone. The b-series reads from the N-terminus; the y-series reads from the C-terminus; together they produce overlapping ladders that together specify the full sequence. Chi and colleagues, writing in the Journal of Proteome Research in 2010, described fragment ion mass calculation from amino acid residue masses as the foundational arithmetic underlying de novo peptide sequencing using HCD spectra. The residue mass of each amino acid is the molecular weight of that amino acid minus water (lost during peptide bond formation): glycine contributes 57.021 Da, alanine 71.037 Da, leucine/isoleucine 113.084 Da, and so forth across all 20 standard residues. Shen and colleagues, in a 2012 paper in the Journal of Proteome Research, described how the b-ion and y-ion series are predicted computationally from the peptide sequence in identification tools that use complementary spectra to improve coverage.
The Isoelectric Point Formula
The isoelectric point (pI) is calculated by finding the pH at which the sum of charges from all ionizable groups equals zero.
Net charge = Σ (charge contribution of each ionizable group at pH x)
pI = pH at which net charge = 0
Each ionizable group (N-terminus, C-terminus, Asp, Glu, His, Cys, Tyr, Lys, Arg) contributes a fractional charge at any given pH, calculated from the Henderson-Hasselbalch equation using the pKa for that group. The algorithm iterates across pH values to find the crossing point. Kozlowski's 2021 paper describing IPC 2.0 in Nucleic Acids Research establishes this as the definitive algorithmic reference for pI and pKa prediction, with validation across multiple pKa scales and peptide datasets. The critical practical issue is which pKa scale the calculator uses. The 2016 IPC validation paper in Biology Direct demonstrated that the IPC scale outperforms EMBOSS and ExPASy scales specifically for 2-DE migration prediction, with lower root-mean-square deviation against observed gel positions.
How to Apply the Fragment Ion Formula
Run through these steps for any peptide sequence using CID or HCD fragmentation.
- Write out the amino acid sequence from N-terminus to C-terminus. Example sequence: GLSDGEWQ (8 residues).
- List the residue masses from the standard monoisotopic residue mass table. G = 57.021, L = 113.084, S = 87.032, D = 115.027, G = 57.021, E = 129.043, W = 186.079, Q = 128.059.
- Calculate b-ions by summing residue masses from the N-terminus, adding 1 for the proton at each position. b1 = 57.021 + 1 = 58.029; b2 = 57.021 + 113.084 + 1 = 171.113; b3 = 57.021 + 113.084 + 87.032 + 1 = 258.145; and so on.
- Calculate y-ions by summing residue masses from the C-terminus, adding 18.011 (water) and 1 (proton) at each position. y1 = 128.059 + 18.011 + 1 = 147.076; y2 = 128.059 + 186.079 + 18.011 + 1 = 333.155; and so on.
- The complete b-series and y-series together span all cleavage positions. For an 8-residue peptide: 7 b-ions and 7 y-ions. The sum of any complementary b/y pair (b_n + y_(length-n)) equals the full precursor mass plus water plus 2 protons.
- For modifications: Add the modification mass delta to every ion in the series that includes the modified residue. Phosphoserine (+79.966 Da) shifts every b-ion from the phosphorylated position onward and every y-ion that includes that position.
Bern and colleagues, writing in the Journal of Computational Biology in 2006, established that mass differences between adjacent peaks in a spectrum correspond to residue masses, which is the computational basis for graph-theory approaches to de novo sequencing.
How to Apply the Isoelectric Point Formula
The manual calculation of pI follows a systematic titration algorithm. For most research purposes, a validated calculator implementing the IPC scale will give the best empirical accuracy for gel-based applications.
- Identify all ionizable groups in the sequence. These are: the N-terminus (pKa approximately 8.0), the C-terminus (pKa approximately 3.1), and the side chains of Asp (3.86), Glu (4.07), His (6.04), Cys (8.14), Tyr (10.46), Lys (10.54), and Arg (12.48) per the IPC scale.
- For each pH value (iterating from 0 to 14): calculate the fractional charge of each ionizable group using the Henderson-Hasselbalch equation: charge = 1 / (1 + 10^(pH - pKa)) for basic groups; charge = -1 / (1 + 10^(pKa - pH)) for acidic groups.
- Sum all fractional charges at each pH to get net charge.
- The pI is the pH at which net charge = 0. Narrow the search using bisection: find the interval where sign of net charge changes, halve it, and repeat to the desired precision.
- For practical gel-based work: the pI predicts where the peptide will focus on an IPG strip. Gorg and colleagues, writing in Proteomics in 2004, established that isoelectric focusing on immobilized pH gradient strips separates proteins and peptides by pI, with the focused position determined by the exact pI relative to the gradient endpoints. Stastna and colleagues, writing in Electrophoresis in 2005, showed that pI precision directly affects band localization accuracy in 2-D gel IEF separations.
Worked Example: Fragment Ion Prediction
The following example uses a hypothetical 5-residue peptide for illustration. It does not represent any clinical compound.
Sequence: ACDEF (5 residues, linear, unmodified)
Monoisotopic residue masses: A = 71.037, C = 103.009, D = 115.027, E = 129.043, F = 147.068.
B-ions (N-terminal series):
- b1 = 71.037 + 1 = 72.044
- b2 = 71.037 + 103.009 + 1 = 175.046
- b3 = 71.037 + 103.009 + 115.027 + 1 = 290.073
- b4 = 71.037 + 103.009 + 115.027 + 129.043 + 1 = 419.116
Y-ions (C-terminal series):
- y1 = 147.068 + 18.011 + 1 = 166.079
- y2 = 147.068 + 129.043 + 18.011 + 1 = 295.122
- y3 = 147.068 + 129.043 + 115.027 + 18.011 + 1 = 410.149
- y4 = 147.068 + 129.043 + 115.027 + 103.009 + 18.011 + 1 = 513.158
Full precursor neutral mass: 71.037 + 103.009 + 115.027 + 129.043 + 147.068 + 18.011 = 583.195 Da.
This example is illustrative only. In practice, database search tools such as Mascot, SEQUEST, and OMSSA embed this arithmetic and match the theoretical series against observed spectra. Good and colleagues, writing in Proteomics in 2010, confirmed that all three major search engines use predicted fragment ion masses as the core matching criterion, with SEQUEST, Mascot, and OMSSA all performing this calculation internally.
Worked Example: Isoelectric Point Calculation
The following uses a hypothetical 4-residue peptide for illustration.
Sequence: GKDE (Gly-Lys-Asp-Glu)
Ionizable groups: N-terminus (pKa 8.0, basic), C-terminus (pKa 3.1, acidic), Lys side chain (pKa 10.54, basic), Asp side chain (pKa 3.86, acidic), Glu side chain (pKa 4.07, acidic).
At pH 4.5: basic groups (N-terminus, Lys) contribute approximately +1.97 total; acidic groups (C-terminus, Asp, Glu) contribute approximately -0.94 total; net charge approximately +1.03.
At pH 7.0: basic groups contribute approximately +1.47; acidic groups contribute approximately -1.98; net charge approximately -0.51.
The sign change between pH 4.5 and pH 7.0 places the pI in that interval. Iterating further narrows it to approximately pH 5.8. The exact value depends on the pKa scale used. Kozlowski's 2016 IPC validation paper showed that pKa scale selection is the single most important variable in pI accuracy for gel migration prediction — the IPC scale reduces mean absolute error versus EMBOSS by approximately 0.3 pI units for peptide datasets.
Fragmentation Method Reference
The fragment ion series produced depends on the fragmentation method. Frese and colleagues' 2011 comparison across CID, HCD, and ETD established the ion-type dependence on fragmentation method with direct empirical data on an LTQ-Orbitrap. Pejchinovski and colleagues, writing in Proteomics: Clinical Applications in 2015, compared HCD and CID specifically for identification of naturally occurring urinary peptides, providing practical guidance on which method performs better in biofluid-based peptidome analysis.
- CID (collision-induced dissociation): Generates primarily b-ions and y-ions. Standard for ion trap and triple-quadrupole instruments. b-ions can be weak or absent for short sequences.
- HCD (higher-energy collisional dissociation): Also generates b-ions and y-ions, with better representation of low-mass ions. Used on Orbitrap instruments. Generally provides better sequence coverage than CID for short peptides.
- ETD (electron-transfer dissociation): Generates c-ions and z-ions instead of b/y. Better for phosphopeptides and large peptides. Fragment ion calculators must use ETD-specific ion formulas.
- ECD (electron-capture dissociation): Similar to ETD; also produces c/z-ions. Less common in modern LC-MS/MS workflows but relevant for top-down proteomics.
Limitations of These Formulas
Both fragment ion and pI formulas apply cleanly to linear, unmodified peptides with standard amino acids. Several classes of peptides require extended handling.
Cyclic and constrained peptides
Cyclic peptides have no N-terminus or C-terminus. Standard b/y-ion prediction does not apply: cleavage of any bond in the ring produces a linear fragment with termini at different positions, generating a complex spectrum that cannot be matched against a standard linear fragment ion series. Mohimani and colleagues, describing Cycloquest in 2011, addressed this limitation with a cycle-specific database search approach. Banerjee and colleagues, writing in Organic and Biomolecular Chemistry in 2011, showed that fragmentation of cyclodepsipeptides in the presence of metal ions produces patterns further diverging from linear peptide predictions. McGee and colleagues, writing in the Journal of Mass Spectrometry in 2013, described facile cleavage C-terminal to ornithine as a non-standard fragmentation event that challenges standard ion series prediction even in linear peptides containing this residue.
Glycopeptides and glycan-bearing modifications
Glycan modifications add substantial mass to fragment ions and generate characteristic oxonium ions (e.g., m/z 204 for HexNAc) that are not part of the standard b/y series. Nilsson, writing in the Glycoconjugate Journal in 2016, described the extended fragmentation analysis required for glycopeptide identification by LC-ESI-MS/MS, including the oxonium ion series and the glycan-modified fragment masses that a standard calculator cannot predict without glycan mass inputs.
Derivatized termini and isobaric labeling
Isobaric termini labeling (IPTL) and other N- or C-terminal derivatizations shift b- and y-ion masses by the mass of the reagent. Koehler and colleagues, writing in Methods in Molecular Biology in 2012, described how isotopic IPTL labels shift both termini-associated ions, requiring the calculator to accept modification masses for both termini simultaneously. Chacon and colleagues, writing in Bioorganic and Medicinal Chemistry in 2006, illustrated how N-terminal amino acid side-chain cleavage in chemically modified peptides generates specific fragmentation patterns used in N-terminal sequencing, where the expected ion series differs from unmodified peptides.
pI limitations for modified or non-standard sequences
The pI formula uses the pKa values of the 20 standard amino acids. Non-standard residues (hydroxyproline, selenocysteine, D-amino acids) have different pKa values not present in standard tables. Chemical modifications that alter ionizable groups (phosphorylation, acetylation of lysine, methylation of arginine) change the effective pKa of those residues. A calculator that does not accept these inputs will return an incorrect pI for modified peptides. For pharmaceutical peptide characterization, pI calculations must use modification-aware algorithms. Qian Cutrone and colleagues' 2017 work on macrocyclic peptide characterization illustrates how full structural characterization of pharmaceutical peptides integrates multiple analytical tools, of which pI calculation is one component.
Application to Pharmaceutical Peptide Characterization
Fragment ion analysis and pI prediction are commonly part of structural characterization workflows for pharmaceutical-grade peptide active pharmaceutical ingredients. In drug development quality control, fragment ion calculators provide the expected spectral signature against which observed MS/MS data is compared to confirm compound identity and sequence integrity. The isoelectric point informs formulation decisions: peptides with pI near the intended formulation pH may have solubility challenges, and buffer selection for reconstitution sometimes accounts for pI. Kozlowski's IPC 2.0 paper noted that pI prediction is relevant not only for gel-based proteomics but also for peptide drug solubility estimation in aqueous formulations.
IGF-1 is the primary downstream biomarker of GH axis activity, and for GH-axis pharmaceutical peptides used under a licensed provider, IGF-1 measurement is a standard monitoring parameter — separate from the research and QC use of the formulas described above.
IMPORTANT SAFETY INFORMATION
This page discusses laboratory and computational tools for peptide analysis in research and pharmaceutical development contexts. The formulas and calculations described are for informational and educational purposes only. This content does not constitute medical advice, a clinical protocol, or a dosing recommendation. Superpower Health provides biomarker testing services through licensed healthcare providers; it does not prescribe or facilitate access to unregulated peptide compounds.
For information about FDA-approved peptide medications, prescribing information is available at dailymed.nlm.nih.gov.
Scientifically reviewed by the Superpower Medical Advisory Team.
