Mathematical models disentangle the role of IL-10 feedbacks in human monocytes upon proinflammatory activation.

Inflammation is one of the vital mechanisms through which the immune system responds to harmful stimuli. During inflammation, proinflammatory and anti-inflammatory cytokines interplay to orchestrate fine-tuned and dynamic immune responses. The cytokine interplay governs switches in the inflammatory response and dictates the propagation and development of the inflammatory response. Molecular pathways underlying the interplay are complex, and time-resolved monitoring of mediators and cytokines is necessary as a basis to study them in detail. Our understanding can be advanced by mathematical models that enable to analyze the system of interactions and their dynamical interplay in detail. We, therefore, used a mathematical modeling approach to study the interplay between prominent proinflammatory and anti-inflammatory cytokines with a focus on tumor necrosis factor and interleukin 10 (IL-10) in lipopolysaccharide-primed primary human monocytes. Relevant time-resolved data were generated by experimentally adding or blocking IL-10 at different time points. The model was successfully trained and could predict independent validation data and was further used to perform simulations to disentangle the role of IL-10 feedbacks during an acute inflammatory event. We used the insight to obtain a reduced predictive model including only the necessary IL-10-mediated feedbacks. Finally, the validated reduced model was used to predict early IL-10-tumor necrosis factor switches in the inflammatory response. Overall, we gained detailed insights into fine-tuning of inflammatory responses in human monocytes and present a model for further use in studying the complex and dynamic process of cytokine-regulated acute inflammation.

A systems biology analysis of adrenergically stimulated adiponectin exocytosis in white adipocytes.

Circulating levels of the adipocyte hormone adiponectin are typically reduced in obesity, and this deficiency has been linked to metabolic diseases. It is thus important to understand the mechanisms controlling adiponectin exocytosis. This understanding is hindered by the high complexity of both the available data and the underlying signaling network. To deal with this complexity, we have previously investigated how different intracellular concentrations of Ca2+, cAMP, and ATP affect adiponectin exocytosis, using both patch-clamp recordings and systems biology mathematical modeling. Recent work has shown that adiponectin exocytosis is physiologically triggered via signaling pathways involving adrenergic ß3 receptors (ß3ARs). Therefore, we developed a mathematical model that also includes adiponectin exocytosis stimulated by extracellular epinephrine or the ß3AR agonist CL 316243. Our new model is consistent with all previous patch-clamp data as well as new data (collected from stimulations with a combination of the intracellular mediators and extracellular adrenergic stimuli) and can predict independent validation data. We used this model to perform new in silico experiments where corresponding wet lab experiments would be difficult to perform. We simulated adiponectin exocytosis in single cells in response to the reduction of ß3ARs that is observed in adipocytes from animals with obesity-induced diabetes. Finally, we used our model to investigate intracellular dynamics and to predict both cAMP levels and adiponectin release by scaling the model from single-cell to a population of cells-predictions corroborated by experimental data. Our work brings us one step closer to understanding the intricate regulation of adiponectin exocytosis.

Perturbation biology links temporal protein changes to drug responses in a melanoma cell line.

Cancer cells have genetic alterations that often directly affect intracellular protein signaling processes allowing them to bypass control mechanisms for cell death, growth and division. Cancer drugs targeting these alterations often work initially, but resistance is common. Combinations of targeted drugs may overcome or prevent resistance, but their selection requires context-specific knowledge of signaling pathways including complex interactions such as feedback loops and crosstalk. To infer quantitative pathway models, we collected a rich dataset on a melanoma cell line: Following perturbation with 54 drug combinations, we measured 124 (phospho-)protein levels and phenotypic response (cell growth, apoptosis) in a time series from 10 minutes to 67 hours. From these data, we trained time-resolved mathematical models that capture molecular interactions and the coupling of molecular levels to cellular phenotype, which in turn reveal the main direct or indirect molecular responses to each drug. Systematic model simulations identified novel combinations of drugs predicted to reduce the survival of melanoma cells, with partial experimental verification. This particular application of perturbation biology demonstrates the potential impact of combining time-resolved data with modeling for the discovery of new combinations of cancer drugs.

Inhibition of FOXO1 transcription factor in primary human adipocytes mimics the insulin-resistant state of type 2 diabetes.

Type 2 diabetes is characterized by insulin resistance in the expanding adipose tissue of obesity. The insulin resistance manifests in human adipocytes as system-wide impairment of insulin signalling. An exception is the regulation of transcription factor FOXO1 (forkhead box protein O1), which is phosphorylated downstream of mTORC2 (mammalian/mechanistic target of rapamycin in complex with raptor) and is therefore not exhibiting impaired response to insulin. However, the abundance, and activity, of FOXO1 is reduced by half in adipocytes from patients with diabetes. To elucidate the effect of reduced FOXO1 activity, we here transduced human adipocytes with a dominant-negative construct of FOXO1 (DN-FOXO1). Inhibition of FOXO1 reduced the abundance of insulin receptor, glucose transporter-4, ribosomal protein S6, mTOR and raptor. Functionally, inhibition of FOXO1 induced an insulin-resistant state network-wide, a state that qualitatively and quantitatively mimicked adipocytes from patients with type 2 diabetes. In contrast, and in accordance with these effects of DN-FOXO1, overexpression of wild-type FOXO1 appeared to augment insulin signalling. We combined experimental data with mathematical modelling to show that the impaired insulin signalling in FOXO1-inhibited cells to a large extent can be explained by reduced mTORC1 activity - a mechanism that defines much of the diabetic state in human adipocytes. Our findings demonstrate that FOXO1 is critical for maintaining normal insulin signalling of human adipocytes.

A systems biology analysis connects insulin receptor signaling with glucose transporter translocation in rat adipocytes.

Type 2 diabetes is characterized by insulin resistance, which arises from malfunctions in the intracellular insulin signaling network. Knowledge of the insulin signaling network is fragmented, and because of the complexity of this network, little consensus has emerged for the structure and importance of the different branches of the network. To help overcome this complexity, systems biology mathematical models have been generated for predicting both the activation of the insulin receptor (IR) and the redistribution of glucose transporter 4 (GLUT4) to the plasma membrane. Although the insulin signal transduction between IR and GLUT4 has been thoroughly studied with modeling and time-resolved data in human cells, comparable analyses in cells from commonly used model organisms such as rats and mice are lacking. Here, we combined existing data and models for rat adipocytes with new data collected for the signaling network between IR and GLUT4 to create a model also for their interconnections. To describe all data (>140 data points), the model needed three distinct pathways from IR to GLUT4: (i) via protein kinase B (PKB) and Akt substrate of 160 kDa (AS160), (ii) via an AS160-independent pathway from PKB, and (iii) via an additional pathway from IR, e.g. affecting the membrane constitution. The developed combined model could describe data not used for training the model and was used to generate predictions of the relative contributions of the pathways from IR to translocation of GLUT4. The combined model provides a systems-level understanding of insulin signaling in rat adipocytes, which, when combined with corresponding models for human adipocytes, may contribute to model-based drug development for diabetes.

Mathematical modeling of white adipocyte exocytosis predicts adiponectin secretion and quantifies the rates of vesicle exo- and endocytosis.

Adiponectin is a hormone secreted from white adipocytes and takes part in the regulation of several metabolic processes. Although the pathophysiological importance of adiponectin has been thoroughly investigated, the mechanisms controlling its release are only partly understood. We have recently shown that adiponectin is secreted via regulated exocytosis of adiponectin-containing vesicles, that adiponectin exocytosis is stimulated by cAMP-dependent mechanisms, and that Ca2+ and ATP augment the cAMP-triggered secretion. However, much remains to be discovered regarding the molecular and cellular regulation of adiponectin release. Here, we have used mathematical modeling to extract detailed information contained within our previously obtained high-resolution patch-clamp time-resolved capacitance recordings to produce the first model of adiponectin exocytosis/secretion that combines all mechanistic knowledge deduced from electrophysiological experimental series. This model demonstrates that our previous understanding of the role of intracellular ATP in the control of adiponectin exocytosis needs to be revised to include an additional ATP-dependent step. Validation of the model by introduction of data of secreted adiponectin yielded a very close resemblance between the simulations and experimental results. Moreover, we could show that Ca2+-dependent adiponectin endocytosis contributes to the measured capacitance signal, and we were able to predict the contribution of endocytosis to the measured exocytotic rate under different experimental conditions. In conclusion, using mathematical modeling of published and newly generated data, we have obtained estimates of adiponectin exo- and endocytosis rates, and we have predicted adiponectin secretion. We believe that our model should have multiple applications in the study of metabolic processes and hormonal control thereof.

LASSIM-A network inference toolbox for genome-wide mechanistic modeling.

Recent technological advancements have made time-resolved, quantitative, multi-omics data available for many model systems, which could be integrated for systems pharmacokinetic use. Here, we present large-scale simulation modeling (LASSIM), which is a novel mathematical tool for performing large-scale inference using mechanistically defined ordinary differential equations (ODE) for gene regulatory networks (GRNs). LASSIM integrates structural knowledge about regulatory interactions and non-linear equations with multiple steady state and dynamic response expression datasets. The rationale behind LASSIM is that biological GRNs can be simplified using a limited subset of core genes that are assumed to regulate all other gene transcription events in the network. The LASSIM method is implemented as a general-purpose toolbox using the PyGMO Python package to make the most of multicore computers and high performance clusters, and is available at https://gitlab.com/Gustafsson-lab/lassim. As a method, LASSIM works in two steps, where it first infers a non-linear ODE system of the pre-specified core gene expression. Second, LASSIM in parallel optimizes the parameters that model the regulation of peripheral genes by core system genes. We showed the usefulness of this method by applying LASSIM to infer a large-scale non-linear model of naïve Th2 cell differentiation, made possible by integrating Th2 specific bindings, time-series together with six public and six novel siRNA-mediated knock-down experiments. ChIP-seq showed significant overlap for all tested transcription factors. Next, we performed novel time-series measurements of total T-cells during differentiation towards Th2 and verified that our LASSIM model could monitor those data significantly better than comparable models that used the same Th2 bindings. In summary, the LASSIM toolbox opens the door to a new type of model-based data analysis that combines the strengths of reliable mechanistic models with truly systems-level data. We demonstrate the power of this approach by inferring a mechanistically motivated, genome-wide model of the Th2 transcription regulatory system, which plays an important role in several immune related diseases.

Systems-wide Experimental and Modeling Analysis of Insulin Signaling through Forkhead Box Protein O1 (FOXO1) in Human Adipocytes, Normally and in Type 2 Diabetes.

Insulin resistance is a major aspect of type 2 diabetes (T2D), which results from impaired insulin signaling in target cells. Signaling to regulate forkhead box protein O1 (FOXO1) may be the most important mechanism for insulin to control transcription. Despite this, little is known about how insulin regulates FOXO1 and how FOXO1 may contribute to insulin resistance in adipocytes, which are the most critical cell type in the development of insulin resistance. We report a detailed mechanistic analysis of insulin control of FOXO1 in human adipocytes obtained from non-diabetic subjects and from patients with T2D. We show that FOXO1 is mainly phosphorylated through mTORC2-mediated phosphorylation of protein kinase B at Ser(473) and that this mechanism is unperturbed in T2D. We also demonstrate a cross-talk from the MAPK branch of insulin signaling to stimulate phosphorylation of FOXO1. The cellular abundance and consequently activity of FOXO1 are halved in T2D. Interestingly, inhibition of mTORC1 with rapamycin reduces the abundance of FOXO1 to the levels in T2D. This suggests that the reduction of the concentration of FOXO1 is a consequence of attenuation of mTORC1, which defines much of the diabetic state in human adipocytes. We integrate insulin control of FOXO1 in a network-wide mathematical model of insulin signaling dynamics based on compatible data from human adipocytes. The diabetic state is network-wide explained by attenuation of an mTORC1-to-insulin receptor substrate-1 (IRS1) feedback and reduced abundances of insulin receptor, GLUT4, AS160, ribosomal protein S6, and FOXO1. The model demonstrates that attenuation of the mTORC1-to-IRS1 feedback is a major mechanism of insulin resistance in the diabetic state.

A single mechanism can explain network-wide insulin resistance in adipocytes from obese patients with type 2 diabetes.

The response to insulin is impaired in type 2 diabetes. Much information is available about insulin signaling, but understanding of the cellular mechanisms causing impaired signaling and insulin resistance is hampered by fragmented data, mainly obtained from different cell lines and animals. We have collected quantitative and systems-wide dynamic data on insulin signaling in primary adipocytes and compared cells isolated from healthy and diabetic individuals. Mathematical modeling and experimental verification identified mechanisms of insulin control of the MAPKs ERK1/2. We found that in human adipocytes, insulin stimulates phosphorylation of the ribosomal protein S6 and hence protein synthesis about equally via ERK1/2 and mTORC1. Using mathematical modeling, we examined the signaling network as a whole and show that a single mechanism can explain the insulin resistance of type 2 diabetes throughout the network, involving signaling both through IRS1, PKB, and mTOR and via ERK1/2 to the nuclear transcription factor Elk1. The most important part of the insulin resistance mechanism is an attenuated feedback from the protein kinase mTORC1 to IRS1, which spreads signal attenuation to all parts of the insulin signaling network. Experimental inhibition of mTORC1 using rapamycin in adipocytes from non-diabetic individuals induced and thus confirmed the predicted network-wide insulin resistance.

Insulin signaling in type 2 diabetes: experimental and modeling analyses reveal mechanisms of insulin resistance in human adipocytes.

Type 2 diabetes originates in an expanding adipose tissue that for unknown reasons becomes insulin resistant. Insulin resistance reflects impairments in insulin signaling, but mechanisms involved are unclear because current research is fragmented. We report a systems level mechanistic understanding of insulin resistance, using systems wide and internally consistent data from human adipocytes. Based on quantitative steady-state and dynamic time course data on signaling intermediaries, normally and in diabetes, we developed a dynamic mathematical model of insulin signaling. The model structure and parameters are identical in the normal and diabetic states of the model, except for three parameters that change in diabetes: (i) reduced concentration of insulin receptor, (ii) reduced concentration of insulin-regulated glucose transporter GLUT4, and (iii) changed feedback from mammalian target of rapamycin in complex with raptor (mTORC1). Modeling reveals that at the core of insulin resistance in human adipocytes is attenuation of a positive feedback from mTORC1 to the insulin receptor substrate-1, which explains reduced sensitivity and signal strength throughout the signaling network. Model simulations with inhibition of mTORC1 are comparable with experimental data on inhibition of mTORC1 using rapamycin in human adipocytes. We demonstrate the potential of the model for identification of drug targets, e.g. increasing the feedback restores insulin signaling, both at the cellular level and, using a multilevel model, at the whole body level. Our findings suggest that insulin resistance in an expanded adipose tissue results from cell growth restriction to prevent cell necrosis.