Host cell proteins (HCPs) are process-related impurities generated during biotherapeutic protein production. HCPs can be problematic if they pose a significant metabolic demand, degrade product quality, or contaminate the final product. Here, we present an effort to create a “clean” Chinese hamster ovary (CHO) cell by disrupting multiple genes to eliminate HCPs. Using a model of CHO cell protein secretion, we predict that the elimination of unnecessary HCPs could have a non-negligible impact on protein production. We analyze the HCP content of 6-protein, 11-protein, and 14-protein knockout clones. These cell lines exhibit a substantial reduction in total HCP content (40%-70%). We also observe higher productivity and improved growth characteristics in specific clones. The reduced HCP content facilitates purification of a monoclonal antibody. Thus, substantial improvements can be made in protein titer and purity through large-scale HCP deletion, providing an avenue to increased quality and affordability of high-value biopharmaceuticals.

Genome-scale metabolic models provide a valuable context for analyzing data from diverse high-throughput experimental techniques. Models can quantify the activities of diverse pathways and cellular functions. Since some metabolic reactions are only catalyzed in specific environments, several algorithms exist that build context-specific models. However, these methods make differing assumptions that influence the content and associated predictive capacity of resulting models, such that model content varies more due to methods used than cell types. Here we overcome this problem with a novel framework for inferring the metabolic functions of a cell before model construction. For this, we curated a list of metabolic tasks and developed a framework to infer the activity of these functionalities from transcriptomic data. We protected the data-inferred tasks during the implementation of diverse context-specific model extraction algorithms for 44 cancer cell lines. We show that the protection of data-inferred metabolic tasks decreases the variability of models across extraction methods. Furthermore, resulting models better capture the actual biological variability across cell lines. This study highlights the potential of using biological knowledge, inferred from omics data, to obtain a better consensus between existing extraction algorithms. It further provides guidelines for the development of the next-generation of data contextualization methods. Genome-scale models of human metabolism have facilitated numerous exciting discoveries regarding human physiology and therapeutics. The accuracy of results from such studies requires that models capture the tissue or cell-type specific metabolism. In hopes to obtain accurate models, several algorithms have been developed to extract cell- or tissue-specific metabolic models. Each algorithm has provided useful insights into the metabolism of specific cell and tissue types. However, since each of these methods use different assumptions to guide reaction inclusion and removal, they result in considerable differences in size, functionality, accuracy, and ultimate biological interpretation, even when using the same data set. To overcome this, the enclosed research proposes an approach to infer the functionalities of a cell or tissue from omics data, and then protect these functions to guide the construction of a context-specific model. Through this study, we highlight the value of using experimental data to help infer the set of metabolic functions that should be included in a model, in an effort to obtain greater consensus across existing extraction algorithms. This study further provides guidelines for the development of the next-generation of data contextualization methods.

Human Milk Oligosaccharides (HMOs) are abundant carbohydrates fundamental to infant health and development. Although these oligosaccharides were discovered more than half a century ago, their biosynthesis in the mammary gland remains largely uncharacterized. Here, we use a systems biology framework that integrates glycan and RNA expression data to construct an HMO biosynthetic network and predict glycosyltransferases involved. To accomplish this, we construct models describing the most likely pathways for the synthesis of the oligosaccharides accounting for >95% of the HMO content in human milk. Through our models, we propose candidate genes for elongation, branching, fucosylation, and sialylation of HMOs. Our model aggregation approach recovers 2 of 2 previously known gene-enzyme relations and 2 of 3 empirically confirmed gene-enzyme relations. The top genes we propose for the remaining 5 linkage reactions are consistent with previously published literature. These results provide the molecular basis of HMO biosynthesis necessary to guide progress in HMO research and application with the goal of understanding and improving infant health and development.

Glycans are fundamental cellular building blocks, involved in many organismal functions. Advances in glycomics are elucidating the essential roles of glycans. Still, it remains challenging to properly analyze large glycomics datasets, since the abundance of each glycan is dependent on many other glycans that share many intermediate biosynthetic steps. Furthermore, the overlap of measured glycans can be low across samples. We address these challenges with GlyCompare, a glycomic data analysis approach that accounts for shared biosynthetic steps for all measured glycans to correct for sparsity and non-independence in glycomics, which enables direct comparison of different glycoprofiles and increases statistical power. Using GlyCompare, we study diverse N-glycan profiles from glycoengineered erythropoietin. We obtain biologically meaningful clustering of mutant cell glycoprofiles and identify knockout-specific effects of fucosyltransferase mutants on tetra-antennary structures. We further analyze human milk oligosaccharide profiles and find mother’s fucosyltransferase-dependent secretor-status indirectly impact the sialylation. Finally, we apply our method on mucin-type O-glycans, gangliosides, and site-specific compositional glycosylation data to reveal tissues and disease-specific glycan presentations. Our substructure-oriented approach will enable researchers to take full advantage of the growing power and size of glycomics data.