In-silico media optimization for continuous cultures using genome scale metabolic networks: The case of CHO-K1.

The cell culture is the central piece of a biotechnological industrial process. It includes upstream (e.g. media preparation, fixed costs, etc.) and downstream steps (e.g. product purification, waste disposal, etc.). In the continuous mode of cell culture, a constant flow of fresh media replaces culture fluid until the system reaches a steady state. This steady state is the standard operation mode which, under very general conditions, is a function of the ratio between the cell density and the dilution rate and depends on the media supplied to the culture. To optimize the production process it is widely accepted that the concentration of the metabolites in this media should be carefully tuned. A poor media may not provide enough nutrients to the culture, while a media too rich in nutrients may be a waste of resources because, either the cells do not use all of the available nutrients, or worse, they over-consume them producing toxic byproducts. In this study, we show how an in-silico study of a genome scale metabolic network coupled to the dynamics of a chemostat could guide the strategy to optimize the media to be used in a continuous process. Given a known media we model the concentrations of the cells in a chemostat as a function of the dilution rate. Then, we cast the problem of optimizing the production process within a linear programming framework in which the goal is to minimize the cost of the media keeping fixed the cell concentration for a given dilution rate in the chemostat. We evaluate our results in two metabolic models: first a simplified model of mammalian cell metabolism, and then in a realistic genome-scale metabolic network of mammalian cells, the Chinese hamster ovary cell line. We explore the latter in more detail given specific meaning to the predictions of the concentrations of several metabolites.

Path-integral solution of MacArthur's resource-competition model for large ecosystems with random species-resources couplings.

We solve MacArthur's resource-competition model with random species-resource couplings in the "thermodynamic" limit of infinitely many species and resources using dynamical path integrals à la De Domincis. We analyze how the steady state picture changes upon modifying several parameters, including the degree of heterogeneity of metabolic strategies (encoding the preferences of species) and of maximal resource levels (carrying capacities), and discuss its stability. Ultimately, the scenario obtained by other approaches is recovered by analyzing an effective one-species-one-resource ecosystem that is fully equivalent to the original multi-species one. The technique used here can be applied for the analysis of other model ecosystems related to the version of MacArthur's model considered here.

Maximum entropy and population heterogeneity in continuous cell cultures.

Continuous cultures of mammalian cells are complex systems displaying hallmark phenomena of nonlinear dynamics, such as multi-stability, hysteresis, as well as sharp transitions between different metabolic states. In this context mathematical models may suggest control strategies to steer the system towards desired states. Although even clonal populations are known to exhibit cell-to-cell variability, most of the currently studied models assume that the population is homogeneous. To overcome this limitation, we use the maximum entropy principle to model the phenotypic distribution of cells in a chemostat as a function of the dilution rate. We consider the coupling between cell metabolism and extracellular variables describing the state of the bioreactor and take into account the impact of toxic byproduct accumulation on cell viability. We present a formal solution for the stationary state of the chemostat and show how to apply it in two examples. First, a simplified model of cell metabolism where the exact solution is tractable, and then a genome-scale metabolic network of the Chinese hamster ovary (CHO) cell line. Along the way we discuss several consequences of heterogeneity, such as: qualitative changes in the dynamical landscape of the system, increasing concentrations of byproducts that vanish in the homogeneous case, and larger population sizes.

Theory of Nonequilibrium Local Search on Random Satisfaction Problems.

We study local search algorithms to solve instances of the random k-satisfiability problem, equivalent to finding (if they exist) zero-energy ground states of statistical models with disorder on random hypergraphs. It is well known that the best such algorithms are akin to nonequilibrium processes in a high-dimensional space. In particular, algorithms known as focused, and which do not obey detailed balance, outperform simulated annealing and related methods in the task of finding the solution to a complex satisfiability problem, that is to find (exactly or approximately) the minimum in a complex energy landscape. A physical question of interest is if the dynamics of these processes can be well predicted by the well-developed theory of equilibrium Gibbs states. While it has been known empirically for some time that this is not the case, an alternative systematic theory that does so has been lacking. In this Letter we introduce such a theory based on the recently developed technique of cavity master equations and test it on the paradigmatic random 3-satisfiability problem. Our theory predicts the qualitative form of the phase boundary between the satisfiable (SAT) and unsatisfiable (UNSAT) region of the phase diagram where the numerics of a focused Metropolis search and cavity master equation cannot be distinguished.

Characterizing steady states of genome-scale metabolic networks in continuous cell cultures.

In the continuous mode of cell culture, a constant flow carrying fresh media replaces culture fluid, cells, nutrients and secreted metabolites. Here we present a model for continuous cell culture coupling intra-cellular metabolism to extracellular variables describing the state of the bioreactor, taking into account the growth capacity of the cell and the impact of toxic byproduct accumulation. We provide a method to determine the steady states of this system that is tractable for metabolic networks of arbitrary complexity. We demonstrate our approach in a toy model first, and then in a genome-scale metabolic network of the Chinese hamster ovary cell line, obtaining results that are in qualitative agreement with experimental observations. We derive a number of consequences from the model that are independent of parameter values. The ratio between cell density and dilution rate is an ideal control parameter to fix a steady state with desired metabolic properties. This conclusion is robust even in the presence of multi-stability, which is explained in our model by a negative feedback loop due to toxic byproduct accumulation. A complex landscape of steady states emerges from our simulations, including multiple metabolic switches, which also explain why cell-line and media benchmarks carried out in batch culture cannot be extrapolated to perfusion. On the other hand, we predict invariance laws between continuous cell cultures with different parameters. A practical consequence is that the chemostat is an ideal experimental model for large-scale high-density perfusion cultures, where the complex landscape of metabolic transitions is faithfully reproduced.

Improving triplet-triplet-annihilation based upconversion systems by tuning their topological structure.

Materials capable to perform upconversion of light transform the photon spectrum and can be used to increase the efficiency of solar cells by upconverting sub-bandgap photons, increasing the density of photons able to generate an electron-hole pair in the cell. Incoherent solar radiation suffices to activate upconverters based on sensitized triplet-triplet annihilation, which makes them particularly suited for this task. This process requires two molecular species, sensitizers absorbing low energy photons, and emitters generating higher frequency photons. Successful implementations exist in solutions and solids. However, solid upconverters exhibit lower efficiency than those in solution, which poses a serious problem for real applications. In the present work, we suggest a new strategy to increase the efficiency of sensitized upconverters that exploits the solid nature of the material. We show that an upconversion model system with molecules distributed as clusters outperforms a system with a random distribution of molecules, as used in current upconverters. Our simulations reveal a high potential for improvement of upconverter systems by exploring different structural configurations of the molecules. The implementation of advanced structures can push the performance of solid upconverters further towards the theoretical limit and a step closer to technological application of low power upconversion.

Sodium acetate increases the productivity of HEK293 cells expressing the ECD-Her1 protein in batch cultures: experimental results and metabolic flux analysis.

Human Embryonic Kidney cells (HEK293) are a popular host for recombinant protein expression and production in the biotechnological industry. This has driven within both, the scientific and the engineering communities, the search for strategies to increase their protein productivity. The present work is inserted into this search exploring the impact of adding sodium acetate (NaAc) into a batch culture of HEK293 cells. We monitored, as a function of time, the cell density, many external metabolites, and the supernatant concentration of the heterologous extra-cellular domain ECD-Her1 protein, a protein used to produce a candidate prostate cancer vaccine. We observed that by adding different concentrations of NaAc (0, 4, 6 and 8 mM), the production of ECD-Her1 protein increases consistently with increasing concentration, whereas the carrying capacity of the medium decreases. To understand these results we exploited a combination of experimental and computational techniques. Metabolic Flux Analysis (MFA) was used to infer intracellular metabolic fluxes from the concentration of external metabolites. Moreover, we measured independently the extracellular acidification rate and oxygen consumption rate of the cells. Both approaches support the idea that the addition of NaAc to the culture has a significant impact on the metabolism of the HEK293 cells and that, if properly tuned, enhances the productivity of the heterologous ECD-Her1 protein.

Optimally designed quantum transport across disordered networks.

We establish a general mechanism for highly efficient quantum transport through finite, disordered 3D networks. It relies on the interplay of disorder with centrosymmetry and a dominant doublet spectral structure and can be controlled by the proper tuning of only coarse-grained quantities. Photosynthetic light harvesting complexes are discussed as potential biological incarnations of this design principle.

Efficiency scaling of non-coherent upconversion in a one-dimensional model system.

A very promising approach to obtain efficient upconversion of light is the use of triplet-triplet annihilation of excitations in molecular systems. In real materials, besides upconversion, many other physical processes take place--fluorescence, phosphorescence, non-radiative decay, annihilation, diffusion--and compete with upconversion. The main objective of this work is to design a proof of principle model that can be used to shed light on the interplay between these processes. Ultimately, we want to establish general principles that may guide experimentalists toward the design of solid state materials with maximum efficiency. Here we show, in a one-dimensional model system, that upconversion can be optimized by varying the ratio between the two molecular species used in triplet-triplet-annihilation based upconversion systems, even in the presence of undesired losses through phosphorescence, non-radiative decay, or annihilation. We derive scaling laws for this ratio and for the maximum efficiency of upconversion, as a function of the diffusion rate J, as well as of the creation and of the decay rate of the excitations.

HEK293 producing the extracellular domain HER1: Full datasets of continuous fermentation process and metabolites analysis.

The data for provide evidences of the multi steady state of the human cell line HEK 293 was obtained from 2 L bioreactor continuous culture. A HEK 293 cell line transfected to produce soluble HER1 receptor was used. The bioreactor was operated at three different dilution rates in sequential manner. Daily samples of culture broth were collected, a total of 85 samples were processed. Viable cell concentration and culture viability was addressing by trypan blue exclusion method using a hemocytometer. Heterologous HER1 supernatant concentration was quantified by a specific ELISA and the metabolites by mass spectrometry coupled to a liquid chromatography. The primary data were collected in excel files, where it was calculated the kinetic and other variables by using mass balance and mathematical principles. It was compared the steady states behavior each other's to find out the existence of steady states' multiplicity, taking into account the stationary phase with respect to the cell density (which means its coefficient of variation is less than 20 %). From the metabolic measurements by using Liquid Chromatography coupled to mass spectrometry (LC-MS), it was also built the data matrix with the specific rates of the 76 metabolites obtained. The data were processed and analyzed, using multivariate data asssnalysis (MVDA) to reduce the complexity and to find the main patterns present in the data. We describe also the full data of the metabolites not only for steady states but also in the time evolution, which could help others in terms of modeling and deep understanding of HEK293 metabolism, especially under different culture conditions.

Stochastic approximation to the T cell mediated specific response of the immune system.

We develop a stochastic model to study the specific response of the immune system. The model is based on the dynamical interaction between Regulatory and Effector CD4+ T cells in the presence of Antigen Presenting Cells inside a lymphatic node. At a mean field level the model predicts the existence of different regimes where active, tolerant, or cyclic immune responses are possible. To study the model beyond mean field and to understand the specific responses of the immune system we use the Linear Noise Approximation and show that fluctuations due to finite size effects may strongly alter the mean field scenario. Moreover, it was found that the existence of a certain characteristic frequency for the fluctuations. All the analytical predictions were compared with simulations using Gillespie's algorithm.

Inference of metabolic fluxes in nutrient-limited continuous cultures: A Maximum Entropy approach with the minimum information.

The study of cellular metabolism is limited by the amount of experimental data available. Formulations able to extract relevant predictions from accessible measurements are needed. Maximum Entropy (ME) inference has been successfully applied to genome-scale models of cellular metabolism, and recent data-driven studies have suggested that in chemostat cultures of Escherichia coli (E. coli), the growth rate and uptake rates of limiting nutrients are the most informative observables. We propose the thesis that this can be explained by the chemostat dynamics, which typically drives nutrient-limited cultures toward observable metabolic states maximally restricted in the dimensions of those fluxes. A practical consequence is that relevant flux observables can now be replaced by culture parameters usually controlled. We test our model by using simulations, and then we apply it to E. coli experimental data where we evaluate the quality of the inference, comparing it to alternative formulations that rest on convex optimization.

Statistical mechanics of interacting metabolic networks.

We cast the metabolism of interacting cells within a statistical mechanics framework considering both the actual phenotypic capacities of each cell and its interaction with its neighbors. Reaction fluxes will be the components of high-dimensional spin vectors, whose values will be constrained by the stochiometry and the energy requirements of the metabolism. Within this picture, finding the phenotypic states of the population turns out to be equivalent to searching for the equilibrium states of a disordered spin model. We provide a general solution of this problem for arbitrary metabolic networks and interactions. We apply this solution to a simplified model of metabolism and to a complex metabolic network, the central core of Escherichia coli, and demonstrate that the combination of selective pressure and interactions defines a complex phenotypic space. We also present numerical results for cells fixed in a grid. These results reproduce the qualitative picture discussed for the mean-field model. Cells may specialize in producing or consuming metabolites complementing each other, and this is described by an equilibrium phase space with multiple minima, like in a spin-glass model.

Modeling functional resting-state brain networks through neural message passing on the human connectome.

In this work, we propose a natural model for information flow in the brain through a neural message-passing dynamics on a structural network of macroscopic regions, such as the human connectome (HC). In our model, each brain region is assumed to have a binary behavior (active or not), the strengths of interactions among them are encoded in the anatomical connectivity matrix defined by the HC, and the dynamics of the system is defined by the Belief Propagation (BP) algorithm, working near the critical point of the network. We show that in the absence of direct external stimuli the BP algorithm converges to a spatial map of activations that is similar to the Default Mode Network (DMN) of the brain, which has been defined from the analysis of functional MRI data. Moreover, we use Susceptibility Propagation (SP) to compute the matrix of long-range correlations between the different regions and show that the modules defined by a clustering of this matrix resemble several Resting State Networks (RSN) determined experimentally. Both results suggest that the functional DMN and RSNs can be seen as simple consequences of the anatomical structure of the brain and a neural message-passing dynamics between macroscopic regions. With the new model, we explore predictions on how functional maps change when the anatomical brain network suffers structural alterations, like in Alzheimer's disease and in lesions of the Corpus Callosum. The implications and novel interpretations suggested by the model, as well as the role of criticality, are discussed.

Cell population heterogeneity driven by stochastic partition and growth optimality.

A fundamental question in biology is how cell populations evolve into different subtypes based on homogeneous processes at the single cell level. Here we show that population bimodality can emerge even when biological processes are homogenous at the cell level and the environment is kept constant. Our model is based on the stochastic partitioning of a cell component with an optimal copy number. We show that the existence of unimodal or bimodal distributions depends on the variance of partition errors and the growth rate tolerance around the optimal copy number. In particular, our theory provides a consistent explanation for the maintenance of aneuploid states in a population. The proposed model can also be relevant for other cell components such as mitochondria and plasmids, whose abundances affect the growth rate and are subject to stochastic partition at cell division.

Estimating the size of the solution space of metabolic networks.

BACKGROUND: Cellular metabolism is one of the most investigated system of biological interactions. While the topological nature of individual reactions and pathways in the network is quite well understood there is still a lack of comprehension regarding the global functional behavior of the system. In the last few years flux-balance analysis (FBA) has been the most successful and widely used technique for studying metabolism at system level. This method strongly relies on the hypothesis that the organism maximizes an objective function. However only under very specific biological conditions (e.g. maximization of biomass for E. coli in reach nutrient medium) the cell seems to obey such optimization law. A more refined analysis not assuming extremization remains an elusive task for large metabolic systems due to algorithmic limitations. RESULTS: In this work we propose a novel algorithmic strategy that provides an efficient characterization of the whole set of stable fluxes compatible with the metabolic constraints. Using a technique derived from the fields of statistical physics and information theory we designed a message-passing algorithm to estimate the size of the affine space containing all possible steady-state flux distributions of metabolic networks. The algorithm, based on the well known Bethe approximation, can be used to approximately compute the volume of a non full-dimensional convex polytope in high dimensions. We first compare the accuracy of the predictions with an exact algorithm on small random metabolic networks. We also verify that the predictions of the algorithm match closely those of Monte Carlo based methods in the case of the Red Blood Cell metabolic network. Then we test the effect of gene knock-outs on the size of the solution space in the case of E. coli central metabolism. Finally we analyze the statistical properties of the average fluxes of the reactions in the E. coli metabolic network. CONCLUSION: We propose a novel efficient distributed algorithmic strategy to estimate the size and shape of the affine space of a non full-dimensional convex polytope in high dimensions. The method is shown to obtain, quantitatively and qualitatively compatible results with the ones of standard algorithms (where this comparison is possible) being still efficient on the analysis of large biological systems, where exact deterministic methods experience an explosion in algorithmic time. The algorithm we propose can be considered as an alternative to Monte Carlo sampling methods.

Exploring the diluted ferromagnetic p-spin model with a cavity master equation.

We introduce an alternative solution to Glauber multispin dynamics on random graphs. The solution is based on the recently introduced cavity master equation (CME), a time-closure turning the, in principle, exact dynamic cavity method into a practical method of analysis and of fast simulation. Running CME once is of comparable computational complexity as one Monte Carlo run on the same problem. We show that CME correctly models the ferromagnetic p-spin Glauber dynamics from high temperatures down to and below the spinoidal transition. We also show that CME allows an alternative exploration of the low-temperature spin-glass phase of the model.

Gauge-free cluster variational method by maximal messages and moment matching.

We present an implementation of the cluster variational method (CVM) as a message passing algorithm. The kind of message passing algorithm used for CVM, usually named generalized belief propagation (GBP), is a generalization of the belief propagation algorithm in the same way that CVM is a generalization of the Bethe approximation for estimating the partition function. However, the connection between fixed points of GBP and the extremal points of the CVM free energy is usually not a one-to-one correspondence because of the existence of a gauge transformation involving the GBP messages. Our contribution is twofold. First, we propose a way of defining messages (fields) in a generic CVM approximation, such that messages arrive on a given region from all its ancestors, and not only from its direct parents, as in the standard parent-to-child GBP. We call this approach maximal messages. Second, we focus on the case of binary variables, reinterpreting the messages as fields enforcing the consistency between the moments of the local (marginal) probability distributions. We provide a precise rule to enforce all consistencies, avoiding any redundancy, that would otherwise lead to a gauge transformation on the messages. This moment matching method is gauge free, i.e., it guarantees that the resulting GBP is not gauge invariant. We apply our maximal messages and moment matching GBP to obtain an analytical expression for the critical temperature of the Ising model in general dimensions at the level of plaquette CVM. The values obtained outperform Bethe estimates, and are comparable with loop corrected belief propagation equations. The method allows for a straightforward generalization to disordered systems.

Microenvironmental cooperation promotes early spread and bistability of a Warburg-like phenotype.

We introduce an in silico model for the initial spread of an aberrant phenotype with Warburg-like overflow metabolism within a healthy homeostatic tissue in contact with a nutrient reservoir (the blood), aimed at characterizing the role of the microenvironment for aberrant growth. Accounting for cellular metabolic activity, competition for nutrients, spatial diffusion and their feedbacks on aberrant replication and death rates, we obtain a phase portrait where distinct asymptotic whole-tissue states are found upon varying the tissue-blood turnover rate and the level of blood-borne primary nutrient. Over a broad range of parameters, the spreading dynamics is bistable as random fluctuations can impact the final state of the tissue. Such a behaviour turns out to be linked to the re-cycling of overflow products by non-aberrant cells. Quantitative insight on the overall emerging picture is provided by a spatially homogeneous version of the model.

Statistical theory of designed quantum transport across disordered networks.

We explain how centrosymmetry, together with a dominant doublet of energy eigenstates in the local density of states, can guarantee interference-assisted, strongly enhanced, strictly coherent quantum excitation transport between two predefined sites of a random network of two-level systems. Starting from a generalization of the chaos-assisted tunnelling mechanism, we formulate a random matrix theoretical framework for the analytical prediction of the transfer time distribution, of lower bounds of the transfer efficiency, and of the scaling behavior of characteristic statistical properties with the size of the network. We show that these analytical predictions compare well to numerical simulations, using Hamiltonians sampled from the Gaussian orthogonal ensemble.