Peptide mapping by capillary zone electrophoresis: how close is theoretical simulation to experimental determination

A multi-variable computer model is presented for the prediction of the electrophoretic mobilities of peptides at pH 2.5 from known physico-chemical constants of their amino acid residues. The model is empirical and does not claim any theoretical dependencies; however, the results suggest that, at least at this pH, peptides may be theoretically represented as classical polymers of freely joined amino acid residues of unequal sizes. The model assumes that the electrophoretic mobility can be represented by a product of three functions that return the contributions of peptide charge, length and width, respectively to the mobility. The model relies on accurate experimental determination of the electrophoretic mobilities of a diverse set of peptides, by capillary zone electrophoresis (CZE), at 22 degrees C, with a 50 mM phosphate buffer, at pH 2.5. The electrophoretic mobilities of a basis set of 102 peptides that varied in charge from 0.65 to 16 and in size from two to 42 amino acid residues were accurately measured at these fixed experimental conditions using a stable 10% linear polyacrylamide-coated column. Data from this basis set was used to derive the peptide charge, length, and width functions respectively. The main purpose of this endeavor is to use the model for the prediction of peptide mobilities at pH 2.5, and for simulation of CZE peptide maps of protein digests. Excellent agreement was obtained between predicted and experimental electrophoretic mobilities for all categories of peptides, including the highly charged and the hydrophobic. To illustrate the utility of this model in protein studies it was used to simulate theoretical peptide maps of the digests of glucagon and horse cytochrome c. The resulting maps were compared and contrasted with their experimental counterparts. The potential of this approach and its limitations are discussed.     

Artificial neural network modeling of peptide mobility and peptide mapping in capillary zone electrophoresis

Recently, we have developed an artificial neural network model, which was able to predict accurately the electrophoretic mobilities of relatively small peptides. To examine the robustness of this methodology, a 3-3-1 back-propagation artificial neural network (BP-ANN) model was developed using the same inputs as the previous model, which were the Offord's charge over mass term (Q/M(2/3)), corrected steric substituent constant (E(s,c)) and molar refractivity (MR). The data set consisted of 102 peptides with a larger range of size than that of our earlier report - up to 42 amino acid residues as compared to 13 amino acids in the initial study - that also included highly charged and hydrophobic peptides. The entire data set was obtained from the published result by Janini and co-workers. The results of this model are compared with those obtained using multiple linear regressions (MLR) model developed in this work and the multi-variable model released by Janini et al. Better predictive ability of the BP-ANN model over the MLR indicates the non-linear characteristics of the electrophoretic mobility of peptides. The present model exhibits better robustness than the MLR models in predicting CZE mobilities of a diverse data set at different experimental conditions. To explore the utility of the ANN model in simulation of the CZE peptide maps, the profiles for the endoproteinase digests of melittin, glucagon and horse cytochrome C is studied in the present work.     

Peptide mapping through the coupling of capillary electrophoresis and high-performance liquid chromatography: map prediction of the tryptic digest of myoglobin.

The tryptic map of horse myoglobin was analysed through capillary electrophoresis using capillaries modified by a monolayer of acrylamide. The results were reproducible and the map was obtained in less than 30 min from ca. 8 pmol of tryptic digest. The peptide identification was performed using peptides previously identified by high-performance liquid chromatography. The peak areas measured using the two techniques are closely related, and the comparison of elution and migration times shows that the two techniques provide different maps. Furthermore, using the semiempirical relationship suggested by Grossman et al. [Anal. Biochem., 179 (1989) 28], which links the electrophoretic mobility to the charge of the peptide and its number of amino acids, a good agreement between predicted and experimental mobilities was observed.     

Peptide maps at picomolar levels obtained by reversed-phase high-performance liquid chromatography and pre-column derivatization with phenyl isothiocyanate. Microsequencing of Phenylthiocarbamyl Peptides

A new reversed-phase high-performance liquid chromatography approach to the production of analytical peptide maps by pre-column derivatization using phenylisothiocyanate is described. Tryptic peptide digests were derivatized with phenyl isothiocyanate to form the phenylthiocarbamyl peptides followed by reversed-phase high-performance liquid chromatographic analysis. The phenylthiocarbamyl peptides were separated by reversed-phase high-performance liquid chromatography with the conventional gradient elution system of water-acetonitrile containing trifluoroacetic acid. The sensitivity of detection of these peptide derivatives was within the range 5-10 pmol with a constant baseline at 254-260 nm. The isolated phenylthiocarbamyl peptides can be subjected to automatic Edman degradation. The effectiveness of this method was exemplified by microsequencing of phenylthiocarbamyl peptides isolated from tryptic digests of three different proteins: alpha-lactalbumin, beta-lactoglobulin and a lambda light-chain immunoglobulin.     

Modeling the electrophoresis of peptides and proteins: improvements in the "bead method" to include ion relaxation and "finite size effects"

A bead model methodology developed in our lab (Xin et al. J. Phys. Chem. B 2006, 110, 1038) and applicable to modeling the free solution electrophoretic mobility of peptides and proteins is generalized in two significant ways. First, an approximate account is taken of the relaxation effect, which makes the methodology applicable to more highly charged peptides and proteins than was previously possible. Second, a more accurate account is taken of the finite size of the beads making up the model structure. This improvement makes the method applicable at higher salt concentrations and/or to models consisting of larger sized subunits. The relaxation effect is accounted for by correcting "unrelaxed" mobilities on the basis of model size and average electrostatic surface, or zeta potential. Correction factors are estimated using those of spheres with the same hydrodynamic radius and zeta potential as the model structure. The correction factors of spheres are readily determined. The more general methodology is first applied to two sets of peptides (74 different peptides total) varying in size from 2 to 42 amino acids. The sets also cover a wide range of net charges. It is shown that accounting for finite bead size results in a small change in model mobilities under the conditions of the experiments (35 mM monovalent salt). The correction for ion relaxation, however, can be significant for highly charged peptides and improves agreement between model and experimental mobilities. Our correction procedure is also tested by examining the electrophoretic mobility of a particular protein "charge ladder" (Carbeck et al. J. Am. Chem. Soc. 1999, 121, 10,671), where the protein charge is varied over a wide range yet the conformation remains essentially constant. In summary, the effects of ion relaxation can be significant if the absolute electrophoretic mobility of a peptide exceeds approximately 0.20 cm2/(kV s).     

Applicability of predictive models to the peptide mobility analysis by capillary electrophoresis-electrospray mass spectrometry

The prediction of peptide mobility by capillary electrophoresis (CE) coupled to electrospray mass spectrometry (MS) is studied in order to verify the validity of the semi-empirical models developed in classical CE. This work relies on the experimental determination of the electrophoretic mobilities of 68 peptides, different in charge and in size. The results indicate that the prediction is possible in CE-MS experiments, in spite of the restraints inherent in the coupling conditions. The best fit of experimental data was obtained with the Offord's model. The efficiency of the model was confirmed by the analysis of a peptide mixture in CE-MS.     

Electrophoretic mobility for peptides with post-translational modifications in capillary electrophoresis

Several semiempirical models for peptide electrophoretic mobility have been tested for capillary electrophoresis (CE) separation with a positively charged capillary using on-line CE combined with electrospray ionization-mass spectrometry (ESI-MS). In the current system with pH 2.7, the expression q/M(0.56) provided the best correlation with the electrophoretic mobility in the analysis of a set of 18 standard samples, where q is the calculated net charge and M is the molecular weight. The peptides resulting from various digests of horse heart myoglobin or bovine hemoglobin were used to demonstrate the validity of this correlation. Post-translationally modified peptides from tryptic digest of human myelin basic protein were also investigated and were found to provide excellent correlation with the linear plot when the total charge of the peptide was correctly calculated. If the total charge was not properly calculated then the post-translationally modified peptides fell off the linear plot. Using this method five arginine residues (residues 5, 49, 54, 97 and 130) were found to be partially citrullinated, four glutamine residues (residues 8, 103, 121 and 147) were found to be partially deamidated and both methionines at residues 21 and 167 were found to be partially oxidized. Three peptides were found with phosphorylation; TPPPSQGK (residue 98 to 105), FSWGAEGQR (residue 114 to 122), and SGSPMAR (residue 163 to 169) and Arginine at 107 was found to be partially monomethylated or dimethylated. The method may provide an excellent means of identifying the presence of peptides with post-translational modifications in conjunction with MS.     

Modeling the electrophoretic mobility and diffusion of weakly charged peptides

A bead model to determine the electrophoretic mobilities and translational diffusion constants of weakly charged peptides is developed that is based on a approximate structural model of peptides and is also grounded in electrohydrodynamic theory. A peptide made up of X amino acids is modeled as N=2X beads with 2 beads representing each amino acid in the chain. For the two beads representing a particular amino acid in a peptide, the radius of one bead is set to one-half the nearest neighbor Calpha-Calpha distance, and the radius of the other bead is chosen on the basis of the diffusion constant of the free amino acid. Peptide conformations, which are defined by a set of psi-phi dihedral angles, are randomly generated by using the transformation matrix approach of Flory (Flory, P. Statistical Mechanics of Chain Molecules; John Wiley: New York, 1969) and rejecting conformations which result in bead overlap. The mobility and diffusion constants are computed for each conformation and at least 100 independent conformations are examined for each peptide. In general, the mobility is found to depend only weakly on peptide conformation. Model and experimental mobilities are compared by examining the data of Janini and co-workers (Janini, G.; et al. J. Chromatogr. 1999, 848, 417-433). A total of 58 peptides consisting of from 2 to 39 amino acids are considered. The average relative error between experimental and model mobilities is found to be 1.0% and the rms relative error 7.7%. In specific cases, the discrepancy can be substantial and possible reasons for this are discussed. It should be emphasized that the input parameters of the peptide model are totally independent of experimental mobilities. It is hoped that the peptide model developed here will be useful in the prediction of peptide mobility as well as in using peptide mobilities to extract information about peptide structure, conformation, and charge. Finally, we show how simultaneous measurements of translational diffusion and mobility can be used to estimate peptide charge.     

Capillary electrophoresis-mass spectrometry of peptides from enzymatic protein hydrolysis: simulation and optimization

Two important limitations still exist for the characterization of protein digests by capillary electrophoresis-mass spectrometry (CE-MS): (i) the buffer choice (i.e., the buffer must provide an adequate CE separation without ruining the MS signal), and (ii) the frequent generation of "unexpected" peptidic fragments during the enzymatic protein hydrolysis. In this work, a new approach is used to solve these difficulties, namely a theoretical model that relates the electrophoretic behavior of peptides to their sequence. The effectiveness of this procedure is demonstrated by the fast attainment of good CE-MS conditions for analyzing the peptides obtained from an enzymatic protein hydrolysate in a single run. This strategy can provide useful information for helping to characterize "unexpected" fragments from protein digests.     

Studies on the relationship between structure and electrophoretic mobility of alpha-helical and beta-sheet peptides using capillary zone electrophoresis.

The electrophoretic behaviour of a series of 33 different synthetic peptides has been investigated using free solution high-performance capillary zonal electrophoretic (HPCZE) methods. The dependency of the electrophoretic mobility, mu(em), on the peptide charge, q, and on the charge-to-size ratio parameter, zeta, determined according to several electromobility models, have been examined. Significant divergences from linearity in the mu(em) vs. q or the mu(em) vs. zeta plots were noted for several peptides, possibly due to the proclivity of specific arrangements of their amino acid sequences to assume preferred alpha-helical or beta-sheet conformational features rather than random coil structures under the HPCZE conditions. These results provide further demonstration of the facility of HPCZE procedures to probe the effects of compositional, sequential and conformational differences of closely-related peptides and their consequences on their physicochemical behaviour in solution.     

An automated method for the analysis of stable isotope labeling data in proteomics

An algorithm is presented for the generation of a reliable peptide component peak table from liquid chromatography-mass spectrometry (LC-MS) and subsequent quantitative analysis of stable isotope coded peptide samples. The method uses chemical noise filtering, charge state fitting, and deisotoping toward improved analysis of complex peptide samples. Overlapping peptide signals in mass spectra were deconvoluted by correlation with modeled peptide isotopic peak profiles. Isotopic peak profiles for peptides were generated in silico from a protein database producing reference model distributions. Doublets of heavy and light labeled peak clusters were identified and compared to provide differential quantification of pairs of stable isotope coded peptides. Algorithms were evaluated using peptides from digests of a single protein and a seven-protein mixture that had been differentially coded with stable isotope labeling agents and mixed in known ratios. The experimental results correlated well with known mixing ratios.     

Capillary electrophoresis/tandem mass spectrometry for analysis of proteins from two-dimensional sodium dodecyl sulfate polyacrylamide gel electrophoresis

Capillary electrophoresis/time-of-flight mass spectrometry(CE/TOFMS) has been used for analysis of in-gel digests of protein spots excised from two-dimensional sodium dodecyl sulfate polyacrylamide gel electrophoresis (2-D SDS-PAGE). An off-line purification and preconcentration procedure with a Zip Tip is used before CE/TOFMS analysis which allows for detection of protein spots with <1 picomole of material from 2-D gels. The off-line procedure provides sufficient purification for analysis while maintaining the quality of the CE separation. Using this procedure, several proteins from Coomassie Blue and zinc negatively stained gels are identified by the peptide maps generated and database searching. CE/TOF tandem mass spectrometry is used for the confirmation of database searching results and structural analysis of peptides that do not match the expected peptide maps obtained from the database in order to identify structural modifications. Several modifications were pinpointed and identified by this method.     

Use of an integrated MS--multiplexed MS/MS data acquisition strategy for high-coverage peptide mapping studies

Liquid chromatography/mass spectrometry (LC/MS) peptide maps have become a basic tool for characterizing proteins of biological and pharmaceutical interest. The ability to generate reproducible maps with high protein sequence coverage is a central goal of methods development. We have applied a recently developed analytical approach (termed LC/MS(E)) to LC/MS peptide mapping. Using the LC/MS(E) approach, the mass detector alternates between a low-energy scanning mode (MS) for accurate mass peptide precursor identification, and an elevated-energy mode (MS(E)) for generation of accurate mass multiplex peptide fragmentation data. In this paper, we evaluate this analytical approach against a tryptic digest of yeast enolase. From the low-energy data, high peptide map coverage (98% of sequence from peptides >3 amino acids) was reproducibly obtained. The MS signal for essentially equimolar peptides varied over 2 orders of magnitude in intensity, and peptide intensities could be precisely and reproducibly measured. Using the temporal constraint that MS(E) peptide fragment ions exhibit chromatographic profiles that parallel the precursor ions that generated them, we were able to produce accurate mass time-resolved MS/MS information for all enolase peptides with sufficient abundance to produce a detectable fragment ion.     

Evaluation of DNA bending models in their capacity to predict electrophoretic migration anomalies of satellite DNA sequences

DNA containing a sequence that generates a local curvature exhibits a pronounced retardation in electrophoretic mobility. Various theoretical models have been proposed to explain relationship between DNA structural features and migration anomaly. Here, we studied the capacity of 15 static wedge‐bending models to predict electrophoretic behavior of 69 satellite monomers derived from four divergent families. All monomers exhibited retarded mobility in PAGE corresponding to retardation factors ranging 1.02–1.54. The curvature varied both within and across the groups and correlated with the number, position, and lengths of A‐tracts. Two dinucleotide models provided strong correlation between gel mobility and curvature prediction; two trinucleotide models were satisfactory while remaining dinucleotide models provided intermediate results with reliable prediction for subsets of sequences only. In some cases, similarly shaped molecules exhibited relatively large differences in mobility and vice versa. Generally less accurate predictions were obtained in groups containing less homogeneous sequences possessing distinct structural features. In conclusion, relatively universal theoretical models were identified suitable for the analysis of natural sequences known to harbor relatively moderate curvature. These models could be potentially applied to genome wide studies. However, in silico predictions should be viewed in context of experimental measurement of intrinsic DNA curvature.     

Use of the arginine-specific butanedione/phenylboronic acid tag for analysis of peptides and protein digests using matrix-assisted laser desorption/ionization mass spectrometry

We have applied an arginine-specific labeling technique to the study of peptides by matrix-assisted laser desorption/ionization mass spectrometry (MALDI-MS). The reaction converts the guanidine group of the arginine side chain by reacting it with 2,3-butanedione and an arylboronic acid. Despite the general chemical lability of the tag under acidic conditions, it was possible to employ acidic matrices like alpha-cyano-4-hydroxycinnamic acid without adverse effects, using the thin-layer technique for preparation. After optimizing the method using arginine-containing model peptides--for which sensitivity down to the low fmol range was demonstrated--the procedure was applied to enzymatic digests of several model proteins in solution and to protein spots in gels obtained by two-dimensional electrophoretic separation of cell lysate samples. Information on the presence of arginine in peptides can be easily obtained from the mass spectra by the characteristic mass shift and the isotope pattern resulting from the incorporation of boron. This information might serve as a valuable additional search constraint for achieving a higher degree of confidence for protein identification by peptide mass fingerprinting.     

Enhancement of ionization efficiency and selective enrichment of phosphorylated peptides from complex protein mixtures using a reversible poly-histidine tag

To improve the detection of phosphorylated peptides/proteins, a combination of optimized MS-based strategies were used involving chemical derivatization with a polyhistidine-tag (His-tag) and affinity enrichment of the resulting His-tag peptides on a nanoscale Ni(2+)-IMAC column. The phosphoserine and phosphothreonine peptides were derivatized using a one-pot beta-elimination/Michael addition reaction with a reversible His-tag possessing a thiol-containing Cys residue. The His-tag peptides were enriched selectively by Ni(2+)-IMAC and released using either imidazole or cleavage with Factor Xa. This novel capture and enzyme-mediated release provided an additional element of selectivity and yielded phosphopeptide-specific modifications with enhanced MS ionization characteristics. The eluted peptides were mapped using MALDI-TOF MS and QTRAP ESI-MS/MS techniques. The results obtained for a model peptide and two tryptic protein digests show that the method is highly specific and allows selective enrichment of phosphorylated peptides at low concentrations of femtomoles per microliter.     

A novel multifunctional labeling reagent for enhanced protein characterization with mass spectrometry.

Individual peptides with lysine at the C-terminus as well as protein tryptic digests were reacted with 2-methoxy-4,5-dihydro-1H-imidazole, converting lysine residues to their 4,5-dihydro-1H-imidazol-2-yl derivatives. The mass spectra of derivatized digests exhibit a greater number of more intense features than their underivatized counterparts, thus increasing the information obtained in peptide mapping experiments. Additionally, MS/MS spectra of the derivatized peptides are greatly simplified in comparison to their native species, yielding primarily an easily interpretable series of y-ions. Finally, this novel label also enables differential quantitation studies, as a stable isotopic form containing four deuterium atoms can readily be produced.     

A database of 660 peptide ion cross sections: Use of intrinsic size parameters for bona fide predictions of cross sections

An ion trap/ion mobility/time-of-flight mass spectrometry technique has been used to measure collision cross sections for 660 peptide ions generated by tryptic digestion of 34 common proteins. Measured cross sections have been compiled into a database that contains peptide molecular weight and sequence information. The database is used to generate average intrinsic contributions to cross section (size parameters) for different amino acid residues by solving systems of equations that relate the unknown contributions of individual residues to the sequences and cross sections of database peptides. Size parameters are combined with information about amino acid composition to calculate cross sections for database peptides. Bona fide cross section predictions (made prior to measurement) for peptides observed in tryptic digests of sperm whale myoglobin and yeast enolase are made. Eight of 10 predicted cross sections are within 2% of the experimental values and all 10 are within 3.2%. The utility of size parameters for cross section prediction is explored and discussed.     

Analysis of tryptic peptides using desorption electrospray ionisation combined with ion mobility spectrometry/mass spectrometry

A novel method is reported for rapid protein identification by the analysis of tryptic peptides using desorption electrospray ionisation (DESI) coupled with hyphenated ion mobility spectrometry and quadrupole time-of-flight mass spectrometry (IMS/Q-ToF-MS). Confident protein identification is demonstrated for the analysis of tryptically digested bovine serum albumin (BSA), with no sample pre-treatment or clean-up. Electrophoretic ion mobility separation of ions generated by DESI allowed examination of charge-state and mobility distributions for tryptic peptide mixtures. Selective interrogation of singly charged ions allowed isobaric peptide responses to be distinguished, along with a reduction in spectral noise. The mobility-selected singly charged peptide responses were presented as a pseudo-peptide mass fingerprint (p-PMF) for protein database searching. Comparative data are shown for electrospray ionisation (ESI) of the BSA digest, without sample clean-up, from which confident protein identification could not be made. Implications for the robustness of the DESI method, together with potential insights into mechanisms for DESI of proteolytic digests, are discussed.