A novel deconvolution method for modeling UDP-N-acetyl-D-glucosamine biosynthetic pathways based on ¹³C mass isotopologue profiles under non-steady-state conditions

Background: Stable isotope tracing is a powerful technique for following the fate of individual atoms through metabolic pathways. Measuring isotopic enrichment in metabolites provides quantitative insights into the biosynthetic network and enables flux analysis as a function of external perturbations. NMR and mass spectrometry are the techniques of choice for global profiling of stable isotope labeling patterns in cellular metabolites. However, meaningful biochemical interpretation of the labeling data requires both quantitative analysis and complex modeling. Here, we demonstrate a novel approach that involved acquiring and modeling the timecourses of ¹³C isotopologue data for UDP-N-acetyl-D-glucosamine (UDP-GlcNAc) synthesized from [U-¹³C]-glucose in human prostate cancer LnCaP-LN3 cells. UDP-GlcNAc is an activated building block for protein glycosylation, which is an important regulatory mechanism in the development of many prominent human diseases including cancer and diabetes.Results: We utilized a stable isotope resolved metabolomics (SIRM) approach to determine the timecourse of ¹³C incorporation from [U-¹³C]-glucose into UDP-GlcNAc in LnCaP-LN3 cells. ¹³C Positional isotopomers and isotopologues of UDP-GlcNAc were determined by high resolution NMR and Fourier transform-ion cyclotron resonance-mass spectrometry. A novel simulated annealing/genetic algorithm, called 'Genetic Algorithm for Isotopologues in Metabolic Systems' (GAIMS) was developed to find the optimal solutions to a set of simultaneous equations that represent the isotopologue compositions, which is a mixture of isotopomer species. The best model was selected based on information theory. The output comprises the timecourse of the individual labeled species, which was deconvoluted into labeled metabolic units, namely glucose, ribose, acetyl and uracil. The performance of the algorithm was demonstrated by validating the computed fractional ¹³C enrichment in these subunits against experimental data. The reproducibility and robustness of the deconvolution were verified by replicate experiments, extensive statistical analyses, and cross-validation against NMR data.Conclusions: This computational approach revealed the relative fluxes through the different biosynthetic pathways of UDP-GlcNAc, which comprises simultaneous sequential and parallel reactions, providing new insight into the regulation of UDP-GlcNAc levels and O-linked protein glycosylation. This is the first such analysis of UDP-GlcNAc dynamics, and the approach is generally applicable to other complex metabolites comprising distinct metabolic subunits, where sufficient numbers of isotopologues can be unambiguously resolved and accurately measured.