Towards a widespread adoption of metabolic modeling tools in biopharmaceutical industry: a process systems biology engineering perspective

Modeling of astaxanthin biosynthesis via machine learning, mathematical and metabolic network modeling

Natural astaxanthin is synthesized by diverse organisms including: bacteria, fungi, microalgae, and plants involving complex cellular processes, which depend on numerous interrelated parameters. Nonetheless, existing knowledge regarding astaxanthin biosynthesis and the conditions influencing astaxanthin accumulation is fairly limited. Thus, manipulation of the growth conditions to achieve desired biomass and astaxanthin yields can be a complicated process requiring cost-intensive and time-consuming experiment-based research. As a potential solution, modeling and simulation of biological systems have recently emerged, allowing researchers to predict/estimate astaxanthin production dynamics in selected organisms. Moreover, mathematical modeling techniques would enable further optimization of astaxanthin synthesis in a shorter period of time, ultimately contributing to a notable reduction in production costs. Thus, the present review comprehensively discusses existing mathematical modeling techniques which simulate the bioaccumulation of astaxanthin in diverse organisms. Associated challenges, solutions, and future perspectives are critically analyzed and presented.

Logistic PCA explains differences between genome-scale metabolic models in terms of metabolic pathways

Genome-scale metabolic models (GSMMs) offer a holistic view of biochemical reaction networks, enabling in-depth analyses of metabolism across species and tissues in multiple conditions. However, comparing GSMMs Against each other poses challenges as current dimensionality reduction algorithms or clustering methods lack mechanistic interpretability, and often rely on subjective assumptions. Here, we propose a new approach utilizing logisitic principal component analysis (LPCA) that efficiently clusters GSMMs while singling out mechanistic differences in terms of reactions and pathways that drive the categorization. We applied LPCA to multiple diverse datasets, including GSMMs of 222 Escherichia-strains, 343 budding yeasts (Saccharomycotina), 80 human tissues, and 2943 Firmicutes strains. Our findings demonstrate LPCA's effectiveness in preserving microbial phylogenetic relationships and discerning human tissue-specific metabolic profiles, exhibiting comparable performance to traditional methods like t-distributed stochastic neighborhood embedding (t-SNE) and Jaccard coefficients. Moreover, the subsystems and associated reactions identified by LPCA align with existing knowledge, underscoring its reliability in dissecting GSMMs and uncovering the underlying drivers of separation.

Giving the cells what they need when they need it: Biosensor-based feeding control

"Giving the cells exactly what they need, when they need it" is the core idea behind the proposed bioprocess control strategy: operating bioprocess based on the physiological behavior of the microbial population rather than exclusive monitoring of environmental parameters. We are envisioning to achieve this through the use of genetically encoded biosensors combined with online flow cytometry (FCM) to obtain a time-dependent "physiological fingerprint" of the population. We developed a biosensor based on the glnA promoter (glnAp) and applied it for monitoring the nitrogen-related nutritional state of Escherichia coli. The functionality of the biosensor was demonstrated through multiple cultivation runs performed at various scales-from microplate to 20 L bioreactor. We also developed a fully automated bioreactor-FCM interface for on-line monitoring of the microbial population. Finally, we validated the proposed strategy by performing a fed-batch experiment where the biosensor signal is used as the actuator for a nitrogen feeding feedback control. This new generation of process control, -based on the specific needs of the cells, -opens the possibility of improving process development on a short timescale and therewith, the robustness and performance of fermentation processes.

Driving towards digital biomanufacturing by CHO genome-scale models

Genome-scale metabolic models (GEMs) of Chinese hamster ovary (CHO) cells are valuable for gaining mechanistic understanding of mammalian cell metabolism and cultures. We provide a comprehensive overview of past and present developments of CHO-GEMs and in silico methods within the flux balance analysis (FBA) framework, focusing on their practical utility in rational cell line development and bioprocess improvements. There are many opportunities for further augmenting the model coverage and establishing integrative models that account for different cellular processes and data for future applications. With supportive collaborative efforts by the research community, we envisage that CHO-GEMs will be crucial for the increasingly digitized and dynamically controlled bioprocessing pipelines, especially because they can be successfully deployed in conjunction with artificial intelligence (AI) and systems engineering algorithms.

CONTRABASS: exploiting flux constraints in genome-scale models for the detection of vulnerabilities

MOTIVATION

Despite the fact that antimicrobial resistance is an increasing health concern, the pace of production of new drugs is slow due to the high cost and uncertain success of the process. The development of high-throughput technologies has allowed the integration of biological data into detailed genome-scale models of multiple organisms. Such models can be exploited by means of computational methods to identify system vulnerabilities such as chokepoint reactions and essential reactions. These vulnerabilities are appealing drug targets that can lead to novel drug developments. However, the current approach to compute these vulnerabilities is only based on topological data and ignores the dynamic information of the model. This can lead to misidentified drug targets.

RESULTS

This work computes flux constraints that are consistent with a certain growth rate of the modelled organism, and integrates the computed flux constraints into the model to improve the detection of vulnerabilities. By exploiting these flux constraints, we are able to obtain a directionality of the reactions of metabolism consistent with a given growth rate of the model, and consequently, a more realistic detection of vulnerabilities can be performed. Several sets of reactions that are system vulnerabilities are defined and the relationships among them are studied. The approach for the detection of these vulnerabilities has been implemented in the Python tool CONTRABASS. Such tool, for which an online web server has also been implemented, computes flux constraints and generates a report with the detected vulnerabilities.

AVAILABILITY AND IMPLEMENTATION

CONTRABASS is available as an open source Python package at https://github.com/openCONTRABASS/CONTRABASS under GPL-3.0 License. An online web server is available at http://contrabass.unizar.es.

SUPPLEMENTARY INFORMATION

A glossary of terms are available at Bioinformatics online.

Machine Learning and Hybrid Methods for Metabolic Pathway Modeling

Computational cell metabolism models seek to provide metabolic explanations of cell behavior under different conditions or following genetic alterations, help in the optimization of in vitro cell growth environments, or predict cellular behavior in vivo and in vitro. In the extremes, mechanistic models can include highly detailed descriptions of a small number of metabolic reactions or an approximate representation of an entire metabolic network. To date, all mechanistic models have required details of individual metabolic reactions, either kinetic parameters or metabolic flux, as well as information about extracellular and intracellular metabolite concentrations. Despite the extensive efforts and the increasing availability of high-quality data, required in vivo data are not available for the majority of known metabolic reactions; thus, mechanistic models are based primarily on ex vivo kinetic measurements and limited flux information. Machine learning approaches provide an alternative for derivation of functional dependencies from existing data. The increasing availability of metabolomic and lipidomic data, with growing feature coverage as well as sample set size, is expected to provide new data options needed for derivation of machine learning models of cell metabolic processes. Moreover, machine learning analysis of longitudinal data can lead to predictive models of cell behaviors over time. Conversely, machine learning models trained on steady-state data can provide descriptive models for the comparison of metabolic states in different environments or disease conditions. Additionally, inclusion of metabolic network knowledge in these analyses can further help in the development of models with limited data.This chapter will explore the application of machine learning to the modeling of cell metabolism. We first provide a theoretical explanation of several machine learning and hybrid mechanistic machine learning methods currently being explored to model metabolism. Next, we introduce several avenues for improving these models with machine learning. Finally, we provide protocols for specific examples of the utilization of machine learning in the development of predictive cell metabolism models using metabolomic data. We describe data preprocessing, approaches for training of machine learning models for both descriptive and predictive models, and the utilization of these models in synthetic and systems biology. Detailed protocols provide a list of software tools and libraries used for these applications, step-by-step modeling protocols, troubleshooting, as well as an overview of existing limitations to these approaches.

The Use of Interdisciplinary Approaches to Understand the Biology of Campylobacter jejuni

Campylobacter jejuni is a bacterial pathogen recognised as a major cause of foodborne illness worldwide. While Campylobacter jejuni generally does not grow outside its host, it can survive outside of the host long enough to pose a health concern. This review presents an up-to-date description and evaluation of biological, mathematical, and statistical approaches used to understand the behaviour of this foodborne pathogen and suggests future avenues which can be explored. Specifically, the incorporation of mathematical modelling may aid the understanding of C. jejuni biofilm formation both outside and inside the host. Predictive studies may be improved by the introduction of more standardised protocols for assessments of disinfection methods and by assessment of novel physical disinfection strategies as well as assessment of the efficiency of plant extracts on C. jejuni eradication. A full description of the metabolic pathways of C. jejuni, which is needed for the successful application of metabolic models, is yet to be achieved. Finally, a shift from animal models (except for those that are a source of human campylobacteriosis) to human-specific data may be made possible due to recent technological advancements, and this may lead to more accurate predictions of human infections.

Arduino Soft Sensor for Monitoring Schizochytrium sp. Fermentation, a Proof of Concept for the Industrial Application of Genome-Scale Metabolic Models in the Context of Pharma 4.0

Schizochytrium sp. is a microorganism cultured for producing docosahexaenoic acid (DHA). Genome-scale metabolic modeling (GEM) is a promising technique for describing gen-protein-reactions in cells, but with still limited industrial application due to its complexity and high computation requirements. In this work, we simplified GEM results regarding the relationship between the specific oxygen uptake rate (−rO2), the specific growth rate (µ), and the rate of lipid synthesis (rL) using an evolutionary algorithm for developing a model that can be used by a soft sensor for fermentation monitoring. The soft sensor estimated the concentration of active biomass (X), glutamate (N), lipids (L), and DHA in a Schizochytrium sp. fermentation using the dissolved oxygen tension (DO) and the oxygen mass transfer coefficient (kLa) as online input variables. The soft sensor model described the biomass concentration response of four reported experiments characterized by different kLa values. The average range normalized root-mean-square error for X, N, L, and DHA were equal to 1.1, 1.3, 1.1, and 3.2%, respectively, suggesting an acceptable generalization capacity. The feasibility of implementing the soft sensor over a low-cost electronic board was successfully tested using an Arduino UNO, showing a novel path for applying GEM-based soft sensors in the context of Pharma 4.0.

Leveraging rAAV bioprocess understanding and next generation bioanalytics development

Recombinant adeno-associated (rAAV) vector-based gene therapy has been the focus of intense research driven by the safety profile and several recent clinical breakthroughs. As of April 2021, there are two rAAV-based gene therapies approved and more than two-hundred active clinical trials (approximately thirty in Phase III). However, the expected increase in demand for rAAV vectors still poses several challenges. Discussed herein are key aspects related to R&D needs and Chemistry, Manufacturing and Control (CMC) efforts required to attend this growing demand. Authors provide their perspective on strategic topics for rAAV-based therapies success: scalability and productivity; improved safety; increased process understanding combined with development of orthogonal bioanalytics that are able to identify, monitor and control Critical Quality Attributes (CQAs) during bioprocessing.

Bacterial metabolism and pathogenesis intimate intertwining: time for metabolic modelling to come into action

We take a snapshot of the recent understanding of bacterial metabolism and the bacterial-host metabolic interplay during infection, and highlight key outcomes and challenges for the practical implementation of bacterial metabolic modelling computational tools in the pathogenesis field.

The Expanding Computational Toolbox for Engineering Microbial Phenotypes at the Genome Scale

Microbial strains are being engineered for an increasingly diverse array of applications, from chemical production to human health. While traditional engineering disciplines are driven by predictive design tools, these tools have been difficult to build for biological design due to the complexity of biological systems and many unknowns of their quantitative behavior. However, due to many recent advances, the gap between design in biology and other engineering fields is closing. In this work, we discuss promising areas of development of computational tools for engineering microbial strains. We define five frontiers of active research: (1) Constraint-based modeling and metabolic network reconstruction, (2) Kinetics and thermodynamic modeling, (3) Protein structure analysis, (4) Genome sequence analysis, and (5) Regulatory network analysis. Experimental and machine learning drivers have enabled these methods to improve by leaps and bounds in both scope and accuracy. Modern strain design projects will require these tools to be comprehensively applied to the entire cell and efficiently integrated within a single workflow. We expect that these frontiers, enabled by the ongoing revolution of big data science, will drive forward more advanced and powerful strain engineering strategies.