Robust Analytical Methods for the Accurate Quantification of the Total Biomass Composition of Mammalian Cells

Protein production dynamics and physiological adaptation of recombinant Komagataella phaffii at near-zero growth rates

BACKGROUND

Specific productivity (qP) in yeast correlates with growth, typically peaking at intermediate or maximum specific growth rates (μ). Understanding the factors limiting productivity at extremely low μ might reveal decoupling strategies, but knowledge of production dynamics and physiology in such conditions is scarce. Retentostats, a type of continuous cultivation, enable the well-controlled transition to near-zero µ through the combined retention of biomass and limited substrate supply. Recombinant Komagataella phaffii (syn Pichia pastoris) secreting a bivalent single domain antibody (VHH) was cultivated in aerobic, glucose-limited retentostats to investigate recombinant protein production dynamics and broaden our understanding of relevant physiological adaptations at near-zero growth conditions.

RESULTS

By the end of the retentostat cultivation, doubling times of approx. two months were reached, corresponding to µ = 0.00047 h-1. Despite these extremely slow growth rates, the proportion of viable cells remained high, and de novo synthesis and secretion of the VHH were observed. The average qP at the end of the retentostat was estimated at 0.019 mg g-1 h-1. Transcriptomics indicated that genes involved in protein biosynthesis were only moderately downregulated towards zero growth, while secretory pathway genes were mostly regulated in a manner seemingly detrimental to protein secretion. Adaptation to near-zero growth conditions of recombinant K. phaffii resulted in significant changes in the total protein, RNA, DNA and lipid content, and lipidomics revealed a complex adaptation pattern regarding the lipid class composition. The higher abundance of storage lipids as well as storage carbohydrates indicates that the cells are preparing for long-term survival.

CONCLUSIONS

In conclusion, retentostat cultivation proved to be a valuable tool to identify potential engineering targets to decouple growth and protein production and gain important insights into the physiological adaptation of K. phaffii to near-zero growth conditions.

Thermodynamic Genome-Scale Metabolic Modeling of Metallodrug Resistance in Colorectal Cancer

Simple Summary

Cancer, but also its treatment, can lead to a reprogramming of cellular metabolism. These changes are observable in metabolite abundances, which can be unbiasedly measured via mass spectrometry metabolomics. However, even when the metabolome changes strongly, a (mechanistic) interpretation is difficult as metabolite levels do not necessarily directly correspond to pathway activities. Here we measure the changes of the cellular metabolome in colorectal cancer cell lines sensitive and resistant to the ruthenium-based drug BOLD-100/KP1339 and the platinum-based drug oxaliplatin. We map these changes onto a cancer-specific genome-scale metabolic model, which allows us not only to compute intracellular flux distributions, but also to disentangle drug-specific effects from growth differences from differences in metabolic adaptations due to resistance. Specifically, we find that resistance to BOLD-100/KP1339 induces more extensive reprogramming than oxaliplatin, especially with respect to fatty acid and amino acid metabolism.

Abstract

Background: Mass spectrometry-based metabolomics approaches provide an immense opportunity to enhance our understanding of the mechanisms that underpin the cellular reprogramming of cancers. Accurate comparative metabolic profiling of heterogeneous conditions, however, is still a challenge. Methods: Measuring both intracellular and extracellular metabolite concentrations, we constrain four instances of a thermodynamic genome-scale metabolic model of the HCT116 colorectal carcinoma cell line to compare the metabolic flux profiles of cells that are either sensitive or resistant to ruthenium- or platinum-based treatments with BOLD-100/KP1339 and oxaliplatin, respectively. Results: Normalizing according to growth rate and normalizing resistant cells according to their respective sensitive controls, we are able to dissect metabolic responses specific to the drug and to the resistance states. We find the normalization steps to be crucial in the interpretation of the metabolomics data and show that the metabolic reprogramming in resistant cells is limited to a select number of pathways. Conclusions: Here, we elucidate the key importance of normalization steps in the interpretation of metabolomics data, allowing us to uncover drug-specific metabolic reprogramming during acquired metal-drug resistance.

Inclusion of maintenance energy improves the intracellular flux predictions of CHO

Chinese hamster ovary (CHO) cells are the leading platform for the production of biopharmaceuticals with human-like glycosylation. The standard practice for cell line generation relies on trial and error approaches such as adaptive evolution and high-throughput screening, which typically take several months. Metabolic modeling could aid in designing better producer cell lines and thus shorten development times. The genome-scale metabolic model (GSMM) of CHO can accurately predict growth rates. However, in order to predict rational engineering strategies it also needs to accurately predict intracellular fluxes. In this work we evaluated the agreement between the fluxes predicted by parsimonious flux balance analysis (pFBA) using the CHO GSMM and a wide range of ¹³C metabolic flux data from literature. While glycolytic fluxes were predicted relatively well, the fluxes of tricarboxylic acid (TCA) cycle were vastly underestimated due to too low energy demand. Inclusion of computationally estimated maintenance energy significantly improved the overall accuracy of intracellular flux predictions. Maintenance energy was therefore determined experimentally by running continuous cultures at different growth rates and evaluating their respective energy consumption. The experimentally and computationally determined maintenance energy were in good agreement. Additionally, we compared alternative objective functions (minimization of uptake rates of seven nonessential metabolites) to the biomass objective. While the predictions of the uptake rates were quite inaccurate for most objectives, the predictions of the intracellular fluxes were comparable to the biomass objective function.

There is an increasing demand for protein pharmaceuticals, especially monoclonal antibodies. Chinese Hamster Ovary (CHO) are currently the leading production host due to their ability to perform human-like post-translational modifications. However, it typically takes several months of trial-and-error approaches to develop a high-producer cell line. Metabolic modelling has the potential to make cell line and process development faster and cheaper by predicting targeted modifications to the cell line genome, cultivation medium or bioprocess. In fact, genome-scale metabolic reconstructions of CHO are already available, and ready for use in cell line development. However, in order to successfully use these models, we need to make sure that they are able to accurately predict metabolic phenotypes. Here we use genome-scale metabolic models of CHO to evaluate the models’ ability to correctly predict intracellular flux distributions. We find that a crucial key ingredient for the correct estimation of central carbon fluxes is the non-growth associated maintenance energy (mATP). We estimated mATP computationally and confirmed it experimentally. Adding this single constraint leads to significantly better predictions of intracellular fluxes, especially in glycolysis and the tricarboxylic acid cycle.

Genome-scale metabolic modelling when changes in environmental conditions affect biomass composition

Genome-scale metabolic modeling is an important tool in the study of metabolism by enhancing the collation of knowledge, interpretation of data, and prediction of metabolic capabilities. A frequent assumption in the use of genome-scale models is that the in vivo organism is evolved for optimal growth, where growth is represented by flux through a biomass objective function (BOF). While the specific composition of the BOF is crucial, its formulation is often inherited from similar organisms due to the experimental challenges associated with its proper determination.

A cell’s macro-molecular composition is not fixed and it responds to changes in environmental conditions. As a consequence, initiatives for the high-fidelity determination of cellular biomass composition have been launched. Thus, there is a need for a mathematical and computational framework capable of using multiple measurements of cellular biomass composition in different environments. Here, we propose two different computational approaches for directly addressing this challenge: Biomass Trade-off Weighting (BTW) and Higher-dimensional-plane InterPolation (HIP).

In lieu of experimental data on biomass composition-variation in response to changing nutrient environment, we assess the properties of BTW and HIP using three hypothetical, yet biologically plausible, BOFs for the Escherichia coli genome-scale metabolic model iML1515. We find that the BTW and HIP formulations have a significant impact on model performance and phenotypes. Furthermore, the BTW method generates larger growth rates in all environments when compared to HIP. Using acetate secretion and the respiratory quotient as proxies for phenotypic changes, we find marked differences between the methods as HIP generates BOFs more similar to a reference BOF than BTW. We conclude that the presented methods constitute a conceptual step in developing genome-scale metabolic modelling approaches capable of addressing the inherent dependence of cellular biomass composition on nutrient environments.

Changes in the environment promote changes in an organism’s metabolism. To achieve balanced growth states for near-optimal function, cells respond through metabolic rearrangements, which may influence the biosynthesis of metabolic precursors for building a cell’s molecular constituents. Therefore, it is necessary to take the dependence of biomass composition on environmental conditions into consideration. While measuring the biomass composition for some environments is possible, and should be done, it cannot be completed for all possible environments.

In this work, we propose two main approaches, BTW and HIP, for addressing the challenge of estimating biomass composition in response to environmental changes. We evaluate the phenotypic consequences of BTW and HIP by characterizing their effect on growth, secretion potential, respiratory efficiency, and gene essentiality of a cell.

Our work constitutes a first conceptual step in accounting for the influence of growth conditions on biomass composition, and in turn the biomass composition’s effect on metabolic phenotypic traits, within constraint-based modelling. As such, we believe it will improve the relevance of constraint-based methods in metabolic engineering and drug discovery, since the biosynthetic potential of microbes for generating industrially relevant products or drugs often is closely linked to their biomass composition.