Glycosylation flux analysis reveals dynamic changes of intracellular glycosylation flux distribution in Chinese hamster ovary fed-batch cultures.

N-linked glycosylation of proteins has both functional and structural significance. Importantly, the glycan structure of a therapeutic protein influences its efficacy, pharmacokinetics, pharmacodynamics and immunogenicity. In this work, we developed glycosylation flux analysis (GFA) for predicting intracellular production and consumption rates (fluxes) of glycoforms, and applied this analysis to CHO fed-batch immunoglobulin G (IgG) production using two different media compositions, with and without additional manganese feeding. The GFA is based on a constraint-based modeling of the glycosylation network, employing a pseudo steady state assumption. While the glycosylation fluxes in the network are balanced at each time point, the GFA allows the fluxes to vary with time by way of two scaling factors: (1) an enzyme-specific factor that captures the temporal changes among glycosylation reactions catalysed by the same enzyme, and (2) the cell specific productivity factor that accounts for the dynamic changes in the IgG production rate. The GFA of the CHO fed-batch cultivations showed that regardless of the media composition, galactosylation fluxes decreased with the cultivation time more significantly than the other glycosylation reactions. Furthermore, the GFA showed that the addition of Mn, a cofactor of galactosyltransferase, has the effect of increasing the galactosylation fluxes but only during the beginning of the cultivation period. The results thus demonstrated the power of the GFA in delineating the dynamic alterations of the glycosylation fluxes by local (enzyme-specific) and global (cell specific productivity) factors.

Assessing and Resolving Model Misspecifications in Metabolic Flux Analysis.

Metabolic flux analysis (MFA) is an indispensable tool in metabolic engineering. The simplest variant of MFA relies on an overdetermined stoichiometric model of the cell's metabolism under the pseudo-steady state assumption to evaluate the intracellular flux distribution. Despite its long history, the issue of model error in overdetermined MFA, particularly misspecifications of the stoichiometric matrix, has not received much attention. We evaluated the performance of statistical tests from linear least square regressions, namely Ramsey's Regression Equation Specification Error Test (RESET), the F-test, and the Lagrange multiplier test, in detecting model misspecifications in the overdetermined MFA, particularly missing reactions. We further proposed an iterative procedure using the F-test to correct such an issue. Using Chinese hamster ovary and random metabolic networks, we demonstrated that: (1) a statistically significant regression does not guarantee high accuracy of the flux estimates; (2) the removal of a reaction with a low flux magnitude can cause disproportionately large biases in the flux estimates; (3) the F-test could efficiently detect missing reactions; and (4) the proposed iterative procedure could robustly resolve the omission of reactions. Our work demonstrated that statistical analysis and tests could be used to systematically assess, detect, and resolve model misspecifications in the overdetermined MFA.