A mathematical model to derive N-glycan structures and cellular enzyme activities from mass spectrometric data

Effective representation and characterization of biosynthetic pathways of glycosylation can be facilitated by mathematical modeling. This paper describes the expansion of a previously developed detailed model for N-linked glycosylation with the further application of the model to analyze MALDI-TOF mass spectra of human N-glycans in terms of underlying cellular enzyme activities. The glycosylation reaction network is automatically generated by the model, based on the reaction specificities of the glycosylation enzymes. The use of a molecular mass cutoff and a network pruning method typically limits the model size to about 10,000 glycan structures. This allows prediction of the complete glycan profile and its abundances for any set of assumed enzyme concentrations and reaction rate parameters. A synthetic mass spectrum from model-calculated glycan profiles is obtained and enzyme concentrations are adjusted to bring the theoretically calculated mass spectrum into agreement with experiment. The result of this process is a complete characterization of a measured glycan mass spectrum containing hundreds of masses in terms of the activities of 19 enzymes. In addition, a complete annotation of the mass spectrum in terms of glycan structure is produced, including the proportions of isomers within each peak. The method was applied to mass spectrometric data of normal human monocytes and monocytic leukemia (THP1) cells to derive glycosyltransferase activity changes underlying the differences in glycan structure between the normal and diseased cells. Model predictions could lead to a better understanding of the changes associated with disease states, identification of disease-associated biomarkers, and bioengineered glycan modifications.