Background: During the last decade, a number of authors have shown that the genetic regulation of metabolic networks may follow optimality principles. Optimal control theory has been succesfully used to compute optimal enzyme profiles considering simple metabolic pathways. However, applying this optimal control framework to more general networks (e.g. branched networks, or networks incorporating enzyme production dynamics) yields problems that are analytically intractable and/or numerically very challenging. Further, these previous studies have only considered a single-objective framework. Results: In this work we consider a more general multi-objective formulation and we present solutions based on recent developments in global dynamic optimization techniques. We illustrate the performance and capabilities of these techniques considering two sets of problems. First, we consider a set of single-objective examples of increasing complexity taken from the recent literature. We analyze the multimodal character of the associated non linear optimization problems, and we also evaluate different global optimization approaches in terms of numerical robustness, efficiency and scalability. Second, we consider generalized multi-objective formulations for several examples, and we show how this framework results in more biologically meaningful results. Conclusions: The proposed strategy was used to solve a set of single-objective case studies related to unbranched and branched metabolic networks of different levels of complexity. All problems were successfully solved in reasonable computation times with our global dynamic optimization approach, reaching solutions which were comparable or better than those reported in previous literature. Further, we considered, for the first time, multi-objective formulations, illustrating how activation in metabolic pathways can be explained in terms of the best trade-offs between conflicting objectives. This new methodology can be applied to metabolic networks with arbitrary topologies, non-linear dynamics and constraints.

As interactions between the immune system and tumour cells are governed by a complex network of cell–cell interactions, knowing the specific immune cell composition of a solid tumour may be essential to predict a patient’s response to immunotherapy. Here, we analyse in depth how to derive the cellular composition of a solid tumour from bulk gene expression data by mathematical deconvolution, using indication-specific and cell type-specific reference gene expression profiles (RGEPs) from tumour-derived single-cell RNA sequencing data. We demonstrate that tumour-derived RGEPs are essential for the successful deconvolution and that RGEPs from peripheral blood are insufficient. We distinguish nine major cell types, as well as three T cell subtypes. Using the tumour-derived RGEPs, we can estimate the content of many tumours associated immune and stromal cell types, their therapeutically relevant ratios, as well as an improved gene expression profile of the malignant cells.

We provide an integrated dynamic view on a eukaryotic osmolyte system, linking signaling with regulation of gene expression, metabolic control and growth. Adaptation to osmotic changes enables cells to adjust cellular activity and turgor pressure to an altered environment. The yeast Saccharomyces cerevisiae adapts to hyperosmotic stress by activating the HOG signaling cascade, which controls glycerol accumulation. The Hog1 kinase stimulates transcription of genes encoding enzymes required for glycerol production (Gpd1, Gpp2) and glycerol import (Stl1) and activates a regulatory enzyme in glycolysis (Pfk26/27). In addition, glycerol outflow is prevented by closure of the Fps1 glycerol facilitator. In order to better understand the contributions to glycerol accumulation of these different mechanisms and how redox and energy metabolism as well as biomass production are maintained under such conditions we collected an extensive dataset. Over a period of 180 min after hyperosmotic shock we monitored in wild type and different mutant cells the concentrations of key metabolites and proteins relevant for osmoadaptation. The dataset was used to parameterize an ODE model that reproduces the generated data very well. A detailed computational analysis using time-dependent response coefficients showed that Pfk26/27 contributes to rerouting glycolytic flux towards lower glycolysis. The transient growth arrest following hyperosmotic shock further adds to redirecting almost all glycolytic flux from biomass towards glycerol production. Osmoadaptation is robust to loss of individual adaptation pathways because of the existence and upregulation of alternative routes of glycerol accumulation. For instance, the Stl1 glycerol importer contributes to glycerol accumulation in a mutant with diminished glycerol production capacity. In addition, our observations suggest a role for trehalose accumulation in osmoadaptation and that Hog1 probably directly contributes to the regulation of the Fps1 glycerol facilitator. Taken together, we elucidated how different metabolic adaptation mechanisms cooperate and provide hypotheses for further experimental studies. Osmotic changes are common environmental challenges for cells, even in multi-cellular organisms, having led to sophisticated adaptation mechanisms. In order to adapt to hyperosmotic stress, yeast cells accumulate glycerol. This is achieved by short-term responses involving metabolic and transmembrane transport changes as well as long-term transcriptional responses. By integrating experimentation and simulation of a mathematical model we resolve the quantitative and temporal characteristics of different processes contributing to glycerol accumulation. We show that osmoadaptation prioritizes the redox and energy balance in glycolysis while rerouting flux from biomass to glycerol production. We further show that the glycerol accumulation network provides osmoadaptation with robustness by compensating for the loss of certain nodes and with the flexibility necessary for responding to different stress situations. Finally we provide novel insight into the roles of transport processes in glycerol accumulation and evidence that trehalose may play a role in yeast osmoadaptation. The present work provides for the first time an integrated dynamic view on a eukaryotic osmolyte system and links signaling with regulation of gene expression and metabolic control.

Parameterized models of biophysical and mechanical cell properties are important for predictive mathematical modeling of cellular processes. The concepts of turgor, cell wall elasticity, osmotically active volume, and intracellular osmolarity have been investigated for decades, but a consistent rigorous parameterization of these concepts is lacking. Here, we subjected several data sets of minimum volume measurements in yeast obtained after hyper-osmotic shock to a thermodynamic modeling framework. We estimated parameters for several relevant biophysical cell properties and tested alternative hypotheses about these concepts using a model discrimination approach. In accordance with previous reports, we estimated an average initial turgor of 0.6 ± 0.2 MPa and found that turgor becomes negligible at a relative volume of 93.3 ± 6.3% corresponding to an osmotic shock of 0.4 ± 0.2 Osm/l. At high stress levels (4 Osm/l), plasmolysis may occur. We found that the volumetric elastic modulus, a measure of cell wall elasticity, is 14.3 ± 10.4 MPa. Our model discrimination analysis suggests that other thermodynamic quantities affecting the intracellular water potential, for example the matrix potential, can be neglected under physiological conditions. The parameterized turgor models showed that activation of the osmosensing high osmolarity glycerol (HOG) signaling pathway correlates with turgor loss in a 1:1 relationship. This finding suggests that mechanical properties of the membrane trigger HOG pathway activation, which can be represented and quantitatively modeled by turgor.

EMBO Press is an editorially independent publishing platform for the development of EMBO scientific publications.