Multi-omic characterization of antibody-producing CHO cell lines elucidates metabolic reprogramming and nutrient uptake bottlenecks.

Characterizing the phenotypic diversity and metabolic capabilities of industrially relevant manufacturing cell lines is critical to bioprocess optimization and cell line development. Metabolic capabilities of production hosts limit nutrient and resource channeling into desired cellular processes and can have a profound impact on productivity. These limitations cannot be directly inferred from measured data such as spent media concentrations or transcriptomics. Here, we present an integrated multi-omic analysis pipeline combining exo-metabolomics, transcriptomics, and genome-scale metabolic network analysis and apply it to three antibody-producing Chinese Hamster Ovary cell lines to identify reprogramming features associated with high-producing clones and metabolic bottlenecks limiting product formation in an industrial bioprocess. Analysis of individual datatypes revealed a decreased nitrogenous byproduct secretion in high-producing clones and the topological changes in peripheral metabolic pathway expression associated with phase shifts. An integrated omics analysis in the context of the genome-scale metabolic model elucidated the differences in central metabolism and identified amino acid utilization bottlenecks limiting cell growth and antibody production that were not evident from exo-metabolomics or transcriptomics alone. Thus, we demonstrate the utility of a multi-omics characterization in providing an in-depth understanding of cellular metabolism, which is critical to efforts in cell engineering and bioprocess optimization.

COSMIC-dFBA: A novel multi-scale hybrid framework for bioprocess modeling.

Metabolism governs cell performance in biomanufacturing, as it fuels growth and productivity. However, even in well-controlled culture systems, metabolism is dynamic, with shifting objectives and resources, thus limiting the predictive capability of mechanistic models for process design and optimization. Here, we present Cellular Objectives and State Modulation In bioreaCtors (COSMIC)-dFBA, a hybrid multi-scale modeling paradigm that accurately predicts cell density, antibody titer, and bioreactor metabolite concentration profiles. Using machine-learning, COSMIC-dFBA decomposes the instantaneous metabolite uptake and secretion rates in a bioreactor into weighted contributions from each cell state (growth or antibody-producing state) and integrates these with a genome-scale metabolic model. A major strength of COSMIC-dFBA is that it can be parameterized with only metabolite concentrations from spent media, although constraining the metabolic model with other omics data can further improve its capabilities. Using COSMIC-dFBA, we can predict the final cell density and antibody titer to within 10% of the measured data, and compared to a standard dFBA model, we found the framework showed a 90% and 72% improvement in cell density and antibody titer prediction, respectively. Thus, we demonstrate our hybrid modeling framework effectively captures cellular metabolism and expands the applicability of dFBA to model the dynamic conditions in a bioreactor.

Advanced error modeling and Bayesian uncertainty quantification in mechanistic liquid chromatography modeling.

Current mechanistic chromatography process modeling methods lack the ability to account for the impact of experimental errors beyond detector noise (e.g. pump delays and variable feed composition) on the uncertainty in calibrated model parameters and the resulting model-predicted chromatograms. This paper presents an uncertainty quantification method that addresses this limitation by determining the probability distribution of parameters in calibrated models, taking into consideration multiple realistic sources of experimental error. The method, which is based on Bayes' theorem and utilizes Markov chain Monte Carlo with an ensemble sampler, is demonstrated to be robust and extensible using synthetic and industrial data. The corresponding software is freely available as open-source code at https://github.com/modsim/CADET-Match.

Guidelines for extracting biologically relevant context-specific metabolic models using gene expression data.

Genome-scale metabolic models comprehensively describe an organism's metabolism and can be tailored using omics data to model condition-specific physiology. The quality of context-specific models is impacted by (i) choice of algorithm and parameters and (ii) alternate context-specific models that equally explain the -omics data. Here we quantify the influence of alternate optima on microbial and mammalian model extraction using GIMME, iMAT, MBA, and mCADRE. We find that metabolic tasks defining an organism's phenotype must be explicitly and quantitatively protected. The scope of alternate models is strongly influenced by algorithm choice and the topological properties of the parent genome-scale model with fatty acid metabolism and intracellular metabolite transport contributing much to alternate solutions in all models. mCADRE extracted the most reproducible context-specific models and models generated using MBA had the most alternate solutions. There were fewer qualitatively different solutions generated by GIMME in E. coli, but these increased substantially in the mammalian models. Screening ensembles using a receiver operating characteristic plot identified the best-performing models. A comprehensive evaluation of models extracted using combinations of extraction methods and expression thresholds revealed that GIMME generated the best-performing models in E. coli, whereas mCADRE is better suited for complex mammalian models. These findings suggest guidelines for benchmarking -omics integration algorithms and motivate the development of a systematic workflow to enumerate alternate models and extract biologically relevant context-specific models.

Advanced score system and automated search strategies for parameter estimation in mechanistic chromatography modeling.

Least squares estimation of unknown parameters from measurement data is a well-established standard method in chromatography modeling but can suffer from critical disadvantages. The description of real-world systems is generally prone to unaccounted mechanisms, such as dispersion in external holdup volumes, and systematic measurement errors, such as caused by pump delays. In this scenario, matching the shape between simulated and measured chromatograms has been found to be more important than the exact peak positions. We have therefore developed a new score system that separately accounts for the shape, position and height of individual peaks. A genetic algorithm is used for optimizing these multiple objectives. Even for non-conflicting objectives, this approach shows superior convergence in comparison to single-objective gradient search, while conflicting objectives indicate incomplete models or inconsistent data. In the latter case, Pareto optima provide important information for understanding the system and improving experiments. The proposed method is demonstrated with synthetic and experimental case studies of increasing complexity. All software is freely available as open source code (https://github.com/modsim/CADET-Match).

Mechanistic model of pH and excipient concentration during ultrafiltration and diafiltration processes of therapeutic antibodies.

Ultrafiltration and diafiltration (UF/DF) unit operations are widely used for the manufacture of therapeutic antibodies to control drug substance protein concentration, pH, and excipient properties. During UF/DF, molecular interactions and volume exclusion effects often lead to substantial differences in pH and excipient concentrations between the diafiltration buffer and final UF/DF pool. These differences complicate the design process beyond simply specifying a buffer with the desired drug substance pH and excipient conditions. This article describes a UF/DF process model which dynamically and accurately simulates UF/DF retentate pool pH and excipient conditions throughout the UF/DF process. This multiscale model accounts for microscopic descriptions of ion-protein charge interactions using the Poisson-Boltzmann equation as well as macroscopic descriptions of volume exclusion and mass transfer. Model predictions of the final UF/DF pool properties were experimentally verified through comparisons to design of experiment (DoE) data from four monoclonal antibody (mAb) processes, each with differing formulations and UF/DF operating conditions. Additionally, model simulations of the retentate pool properties throughout the UF/DF process were verified for two mAb processes through comparisons to experimental data collected at intermediate process points. Model results were qualified, using statistical equivalence tests, against the outputs from large-scale GMP runs which confirmed that the model accurately captures large-scale process performance. Finally, the model was applied toward the simulation of process scenarios beyond those examined experimentally. These in-silico experiments demonstrate the model's capability as a tool for augmented process design and it is potential to reduce the extent of UF/DF laboratory experiments.