skip to main content

Title: Integration of Time-Series Transcriptomic Data with Genome-Scale CHO Metabolic Models for mAb Engineering
Chinese hamster ovary (CHO) cells are the most commonly used cell lines in biopharmaceutical manufacturing. Genome-scale metabolic models have become a valuable tool to study cellular metabolism. Despite the presence of reference global genome-scale CHO model, context-specific metabolic models may still be required for specific cell lines (for example, CHO-K1, CHO-S, and CHO-DG44), and for specific process conditions. Many integration algorithms have been available to reconstruct specific genome-scale models. These methods are mainly based on integrating omics data (i.e., transcriptomics, proteomics, and metabolomics) into reference genome-scale models. In the present study, we aimed to investigate the impact of time points of transcriptomics integration on the genome-scale CHO model by assessing the prediction of growth rates with each reconstructed model. We also evaluated the feasibility of applying extracted models to different cell lines (generated from the same parental cell line). Our findings illustrate that gene expression at various stages of culture slightly impacts the reconstructed models. However, the prediction capability is robust enough on cell growth prediction not only across different growth phases but also in expansion to other cell lines.
Award ID(s):
Publication Date:
Journal Name:
Page Range or eLocation-ID:
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Constraint-based modeling has been applied to analyze metabolism of numerous organisms via flux balance analysis and genome-scale metabolic models, including mammalian cells such as the Chinese hamster ovary (CHO) cells—the principal cell factory platform for therapeutic protein production. Unfortunately, the application of genome-scale model methodologies using the conventional biomass objective function is challenged by the presence of overly-restrictive constraints, including essential amino acid exchange fluxes that can lead to improper predictions of growth rates and intracellular flux distributions. In this study, these constraints are found to be reliably predicted by an “essential nutrient minimization” approach. After modifying these constraints with the predicted minimal uptake values, a series of unconventional objective functions are applied to minimize each individual non-essential nutrient uptake rate, revealing useful insights about metabolic exchange rates and flows across different cell lines and culture conditions. This unconventional uptake-rate objective functions (UOFs) approach is able to distinguish metabolic differences between three distinct CHO cell lines (CHO-K1, -DG44, and -S) not directly observed using the conventional biomass growth maximization solutions. Further, a comparison of model predictions with experimental data from literature correctly correlates with the specific CHO-DG44-derived cell line used experimentally, and the corresponding dual prices provide fruitful informationmore »concerning coupling relationships between nutrients. The UOFs approach is likely to be particularly suited for mammalian cells and other complex organisms which contain multiple distinct essential nutrient inputs, and may offer enhanced applicability for characterizing cell metabolism and physiology as well as media optimization and biomanufacturing control.

    « less
  2. Hallam, Steven J. (Ed.)
    ABSTRACT Nitrite-oxidizing bacteria belonging to the genus Nitrospira mediate a key step in nitrification and play important roles in the biogeochemical nitrogen cycle and wastewater treatment. While these organisms have recently been shown to exhibit metabolic flexibility beyond their chemolithoautotrophic lifestyle, including the use of simple organic compounds to fuel their energy metabolism, the metabolic networks controlling their autotrophic and mixotrophic growth remain poorly understood. Here, we reconstructed a genome-scale metabolic model for Nitrospira moscoviensis ( i Nmo686) and used flux balance analysis to evaluate the metabolic networks controlling autotrophic and formatotrophic growth on nitrite and formate, respectively. Subsequently, proteomic analysis and [ 13 C]bicarbonate and [ 13 C]formate tracer experiments coupled to metabolomic analysis were performed to experimentally validate model predictions. Our findings corroborate that N. moscoviensis uses the reductive tricarboxylic acid cycle for CO 2 fixation, and we also show that N. moscoviensis can indirectly use formate as a carbon source by oxidizing it first to CO 2 followed by reassimilation, rather than direct incorporation via the reductive glycine pathway. Our study offers the first measurements of Nitrospira ’s in vivo central carbon metabolism and provides a quantitative tool that can be used for understanding and predicting theirmore »metabolic processes. IMPORTANCE Nitrospira spp. are globally abundant nitrifying bacteria in soil and aquatic ecosystems and in wastewater treatment plants, where they control the oxidation of nitrite to nitrate. Despite their critical contribution to nitrogen cycling across diverse environments, detailed understanding of their metabolic network and prediction of their function under different environmental conditions remains a major challenge. Here, we provide the first constraint-based metabolic model of Nitrospira moscoviensis representing the ubiquitous Nitrospira lineage II and subsequently validate this model using proteomics and 13 C-tracers combined with intracellular metabolomic analysis. The resulting genome-scale model will serve as a knowledge base of Nitrospira metabolism and lays the foundation for quantitative systems biology studies of these globally important nitrite-oxidizing bacteria.« less
  3. Optimization-based models have been used to predict cellular behavior for over 25 years. The constraints in these models are derived from genome annotations, measured macromolecular composition of cells, and by measuring the cell's growth rate and metabolism in different conditions. The cellular goal (the optimization problem that the cell is trying to solve) can be challenging to derive experimentally for many organisms, including human or mammalian cells, which have complex metabolic capabilities and are not well understood. Existing approaches to learning goals from data include (a) estimating a linear objective function, or (b) estimating linear constraints that model complex biochemical reactions and constrain the cell's operation. The latter approach is important because often the known reactions are not enough to explain observations; therefore, there is a need to extend automatically the model complexity by learning new reactions. However, this leads to nonconvex optimization problems, and existing tools cannot scale to realistically large metabolic models. Hence, constraint estimation is still used sparingly despite its benefits for modeling cell metabolism, which is important for developing novel antimicrobials against pathogens, discovering cancer drug targets, and producing value-added chemicals. Here, we develop the first approach to estimating constraint reactions from data that can scalemore »to realistically large metabolic models. Previous tools were used on problems having less than 75 reactions and 60 metabolites, which limits real-life-size applications. We perform extensive experiments using 75 large-scale metabolic network models for different organisms (including bacteria, yeasts, and mammals) and show that our algorithm can recover cellular constraint reactions. The recovered constraints enable accurate prediction of metabolic states in hundreds of growth environments not seen in training data, and we recover useful cellular goals even when some measurements are missing.« less
  4. Abstract MicroRNAs (miRNAs) play important roles in post-transcriptional gene regulation and phenotype development. Understanding the regulation of miRNA genes is critical to understand gene regulation. One of the challenges to study miRNA gene regulation is the lack of condition-specific annotation of miRNA transcription start sites (TSSs). Unlike protein-coding genes, miRNA TSSs can be tens of thousands of nucleotides away from the precursor miRNAs and they are hard to be detected by conventional RNA-Seq experiments. A number of studies have been attempted to computationally predict miRNA TSSs. However, high-resolution condition-specific miRNA TSS prediction remains a challenging problem. Recently, deep learning models have been successfully applied to various bioinformatics problems but have not been effectively created for condition-specific miRNA TSS prediction. Here we created a two-stream deep learning model called D-miRT for computational prediction of condition-specific miRNA TSSs ( ). D-miRT is a natural fit for the integration of low-resolution epigenetic features (DNase-Seq and histone modification data) and high-resolution sequence features. Compared with alternative computational models on different sets of training data, D-miRT outperformed all baseline models and demonstrated high accuracy for condition-specific miRNA TSS prediction tasks. Comparing with the most recent approaches on cell-specific miRNA TSS identification using cell linesmore »that were unseen to the model training processes, D-miRT also showed superior performance.« less
  5. Gralnick, Jeffrey A. (Ed.)
    ABSTRACT Rhodopseudomonas palustris CGA009 is a Gram-negative purple nonsulfur bacterium that grows phototrophically by fixing carbon dioxide and nitrogen or chemotrophically by fixing or catabolizing a wide array of substrates, including lignin breakdown products for its carbon and fixing nitrogen for its nitrogen requirements. It can grow aerobically or anaerobically and can use light, inorganic, and organic compounds for energy production. Due to its ability to convert different carbon sources into useful products during anaerobic growth, this study reconstructed a metabolic and expression (ME) model of R. palustris to investigate its anaerobic-photoheterotrophic growth. Unlike metabolic (M) models, ME models include transcription and translation reactions along with macromolecules synthesis and couple these reactions with growth rate. This unique feature of the ME model led to nonlinear growth curve predictions, which matched closely with experimental growth rate data. At the theoretical maximum growth rate, the ME model suggested a diminishing rate of carbon fixation and predicted malate dehydrogenase and glycerol-3 phosphate dehydrogenase as alternate electron sinks. Moreover, the ME model also identified ferredoxin as a key regulator in distributing electrons between major redox balancing pathways. Because ME models include the turnover rate for each metabolic reaction, it was used to successfully capturemore »experimentally observed temperature regulation of different nitrogenases. Overall, these unique features of the ME model demonstrated the influence of nitrogenases and rubiscos on R. palustris growth and predicted a key regulator in distributing electrons between major redox balancing pathways, thus establishing a platform for in silico investigation of R. palustris metabolism from a multiomics perspective. IMPORTANCE In this work, we reconstructed the first ME model for a purple nonsulfur bacterium (PNSB). Using the ME model, different aspects of R. palustris metabolism were examined. First, the ME model was used to analyze how reducing power entering the R. palustris cell through organic carbon sources gets partitioned into biomass, carbon dioxide fixation, and nitrogen fixation. Furthermore, the ME model predicted electron flux through ferredoxin as a major bottleneck in distributing electrons to nitrogenase enzymes. Next, the ME model characterized different nitrogenase enzymes and successfully recapitulated experimentally observed temperature regulations of those enzymes. Identifying the bottleneck responsible for transferring an electron to nitrogenase enzymes and recapitulating the temperature regulation of different nitrogenase enzymes can have profound implications in metabolic engineering, such as hydrogen production from R. palustris . Another interesting application of this ME model can be to take advantage of its redox balancing strategy to gain an understanding of the regulatory mechanism of biodegradable plastic production precursors, such as polyhydroxybutyrate (PHB).« less