Comparing genome-scale metabolic models of the non-resistant Enterococcus faecalis ATCC 19433 and the multi-resistant Enterococcus faecalis V583.

Enterococcus faecalis is a versatile lactic acid bacterium with a large variety of implications for humans. While some strains of this species are pathobionts being resistant against most of the common antibiotics, other strains are regarded as biological protectants or even probiotics. Accordingly, E. faecalis strains largely differ in the size and content of their accessory genome. In this study, we describe the genome-scale metabolic network reconstruction of E. faecalis ATCC 19433, a non-resistant human-associated strain. A comparison of the genome-scale metabolic model (GSM) of E. faecalis ATCC 19433 with a previously published GSM of the multi-resistant pathobiontic E. faecalis V583 reveals high similarities in the central metabolic abilities of these two human associated strains. This is reflected, e.g., in the identical amino acid auxotrophies. The ATCC 19433 strain, however, has a 14.1 % smaller genome than V583 and lacks the multiple antibiotic resistance genes and genes involved in capsule formation. Based on the measured metabolic fluxes at different growth rates, the energy demand at zero growth was calculated to be about 40 % lower for the ATCC 19433 strain compared to V583. Furthermore, the ATCC 19433 strain seems less prone to the depletion of amino acids utilizable for energy metabolism. This might hint at a lower overall energy demand of the ATCC 19433 strain as compared to V583.

The pH-dependent lactose metabolism of Lactobacillus delbrueckii subsp. bulgaricus: An integrative view through a mechanistic computational model.

The fermentation process of milk to yoghurt using Lactobacillus delbrueckii subsp. bulgaricus in co-culture with Streptococcus thermophilus is hallmarked by the breakdown of lactose to organic acids such as lactate. This leads to a substantial decrease in pH - both in the medium, as well as cytosolic. The latter impairs metabolic activities due to the pH-dependence of enzymes, which compromises microbial growth. To quantitatively elucidate the impact of the acidification on metabolism of L. bulgaricus in an integrated way, we have developed a proton-dependent computational model of lactose metabolism and casein degradation based on experimental data. The model accounts for the influence of pH on enzyme activities as well as cellular growth and proliferation of the bacterial population. We used a machine learning approach to quantify the cell volume throughout fermentation. Simulation results show a decrease in metabolic flux with acidification of the cytosol. Additionally, the validated model predicts a similar metabolic behaviour within a wide range of non-limiting substrate concentrations. This computational model provides a deeper understanding of the intricate relationships between metabolic activity and acidification and paves the way for further optimization of yoghurt production under industrial settings.

Spatial modeling reveals nuclear phosphorylation and subcellular shuttling of YAP upon drug-induced liver injury.

The Hippo signaling pathway controls cell proliferation and tissue regeneration via its transcriptional effectors yes-associated protein (YAP) and transcriptional coactivator with PDZ-binding motif (TAZ). The canonical pathway topology is characterized by sequential phosphorylation of kinases in the cytoplasm that defines the subcellular localization of YAP and TAZ. However, the molecular mechanisms controlling the nuclear/cytoplasmic shuttling dynamics of both factors under physiological and tissue-damaging conditions are poorly understood. By implementing experimental in vitro data, partial differential equation modeling, as well as automated image analysis, we demonstrate that nuclear phosphorylation contributes to differences between YAP and TAZ localization in the nucleus and cytoplasm. Treatment of hepatocyte-derived cells with hepatotoxic acetaminophen (APAP) induces a biphasic protein phosphorylation eventually leading to nuclear protein enrichment of YAP but not TAZ. APAP-dependent regulation of nuclear/cytoplasmic YAP shuttling is not an unspecific cellular response but relies on the sequential induction of reactive oxygen species (ROS), RAC-alpha serine/threonine-protein kinase (AKT, synonym: protein kinase B), as well as elevated nuclear interaction between YAP and AKT. Mouse experiments confirm this sequence of events illustrated by the expression of ROS-, AKT-, and YAP-specific gene signatures upon APAP administration. In summary, our data illustrate the importance of nuclear processes in the regulation of Hippo pathway activity. YAP and TAZ exhibit different shuttling dynamics, which explains distinct cellular responses of both factors under physiological and tissue-damaging conditions.

Stochastic dynamics of Type-I interferon responses.

Interferon (IFN) activates the transcription of several hundred of IFN stimulated genes (ISGs) that constitute a highly effective antiviral defense program. Cell-to-cell variability in the induction of ISGs is well documented, but its source and effects are not completely understood. The molecular mechanisms behind this heterogeneity have been related to randomness in molecular events taking place during the JAK-STAT signaling pathway. Here, we study the sources of variability in the induction of the IFN-alpha response by using MxA and IFIT1 activation as read-out. To this end, we integrate time-resolved flow cytometry data and stochastic modeling of the JAK-STAT signaling pathway. The complexity of the IFN response was matched by fitting probability distributions to time-course flow cytometry snapshots. Both, experimental data and simulations confirmed that the MxA and IFIT1 induction circuits generate graded responses rather than all-or-none responses. Subsequently, we quantify the size of the intrinsic variability at different steps in the pathway. We found that stochastic effects are transiently strong during the ligand-receptor activation steps and the formation of the ISGF3 complex, but negligible for the final induction of the studied ISGs. We conclude that the JAK-STAT signaling pathway is a robust biological circuit that efficiently transmits information under stochastic environments.

Computational modeling and quantitative physiology reveal central parameters for brassinosteroid-regulated early cell physiological processes linked to elongation growth of the Arabidopsis root.

Brassinosteroids (BR) are key hormonal regulators of plant development. However, whereas the individual components of BR perception and signaling are well characterized experimentally, the question of how they can act and whether they are sufficient to carry out the critical function of cellular elongation remains open. Here, we combined computational modeling with quantitative cell physiology to understand the dynamics of the plasma membrane (PM)-localized BR response pathway during the initiation of cellular responses in the epidermis of the Arabidopsis root tip that are be linked to cell elongation. The model, consisting of ordinary differential equations, comprises the BR-induced hyperpolarization of the PM, the acidification of the apoplast and subsequent cell wall swelling. We demonstrate that the competence of the root epidermal cells for the BR response predominantly depends on the amount and activity of H+-ATPases in the PM. The model further predicts that an influx of cations is required to compensate for the shift of positive charges caused by the apoplastic acidification. A potassium channel was subsequently identified and experimentally characterized, fulfilling this function. Thus, we established the landscape of components and parameters for physiological processes potentially linked to cell elongation, a central process in plant development.

All Driven by Energy Demand? Integrative Comparison of Metabolism of Enterococcus faecalis Wildtype and a Glutamine Synthase Mutant.

Lactic acid bacteria (LAB) play a significant role in biotechnology, e.g., food industry and also in human health. Many LAB genera have developed a multidrug resistance in the past few years, causing a serious problem in controlling hospital germs worldwide. Enterococcus faecalis accounts for a large part of the human infections caused by LABs. Therefore, studying its adaptive metabolism under various environmental conditions is particularly important to promote the development of new therapeutic approaches. In this study, we investigated the effect of glutamine auxotrophy (ΔglnA mutant) on metabolic and proteomic adaptations of E. faecalis in response to a changing pH in its environment. Changing pH values are part of the organism's natural environment in the human body and play a role in the food industry. We compared the results with those of the wildtype. Using a genome-scale metabolic model constrained by metabolic and proteomic data, our integrative method allows us to understand the bigger picture of the adaptation strategies of this bacterium. The study showed that energy demand is the decisive factor in adapting to a new environmental pH. The energy demand of the mutant was higher at all conditions. It has been reported that ΔglnA mutants of bacteria are energetically less effective. With the aid of our data and model we are able to explain this phenomenon as a consequence of a failure to regulate glutamine uptake and the costs for the import of glutamine and the export of ammonium. Methodologically, it became apparent that taking into account the nonspecificity of amino acid transporters is important for reproducing metabolic changes with genome-scale models because it affects energy balance. IMPORTANCE The integration of new pH-dependent experimental data on metabolic uptake and release fluxes, as well as of proteome data with a genome-scale computational model of a glutamine synthetase mutant of E. faecalis is used and compared with those of the wildtype to understand why glutamine auxotrophy results in a less efficient metabolism and how-in comparison with the wildtype-the glutamine synthetase knockout impacts metabolic adjustments during acidification or simply exposure to lower pH. We show that forced glutamine auxotrophy causes more energy demand and that this is likely due to a disregulated glutamine uptake. Proteome changes during acidification observed for the mutant resemble those of the wildtype with the exception of glycolysis-related genes, as the mutant is already energetically stressed at a higher pH and the respective proteome changes were in effect.

Inspecting the Solution Space of Genome-Scale Metabolic Models.

Genome-scale metabolic models are frequently used in computational biology. They offer an integrative view on the metabolic network of an organism without the need to know kinetic information in detail. However, the huge solution space which comes with the analysis of genome-scale models by using, e.g., Flux Balance Analysis (FBA) poses a problem, since it is hard to thoroughly investigate and often only an arbitrarily selected individual flux distribution is discussed as an outcome of FBA. Here, we introduce a new approach to inspect the solution space and we compare it with other approaches, namely Flux Variability Analysis (FVA) and CoPE-FBA, using several different genome-scale models of lactic acid bacteria. We examine the extent to which different types of experimental data limit the solution space and how the robustness of the system increases as a result. We find that our new approach to inspect the solution space is a good complementary method that offers additional insights into the variance of biological phenotypes and can help to prevent wrong conclusions in the analysis of FBA results.

An integrative view on vacuolar pH homeostasis in Arabidopsis thaliana: Combining mathematical modeling and experimentation.

The acidification of plant vacuoles is of great importance for various physiological processes, as a multitude of secondary active transporters utilize the proton gradient established across the vacuolar membrane. Vacuolar-type H+ -translocating ATPases and a pyrophosphatase are thought to enable vacuoles to accumulate protons against their electrochemical potential. However, recent studies pointed to the ATPase located at the trans-Golgi network/early endosome (TGN/EE) to contribute to vacuolar acidification in a manner not understood as of now. Here, we combined experimental data and computational modeling to test different hypotheses for vacuolar acidification mechanisms. For this, we analyzed different models with respect to their ability to describe existing experimental data. To better differentiate between alternative acidification mechanisms, new experimental data have been generated. By fitting the models to the experimental data, we were able to prioritize the hypothesis in which vesicular trafficking of Ca2+ /H+ -antiporters from the TGN/EE to the vacuolar membrane and the activity of ATP-dependent Ca2+ -pumps at the tonoplast might explain the residual acidification observed in Arabidopsis mutants defective in vacuolar proton pump activity. The presented modeling approach provides an integrative perspective on vacuolar pH regulation in Arabidopsis and holds potential to guide further experimental work.

Impact of explicit area scaling on kinetic models involving multiple compartments.

BACKGROUND: Computational modelling of cell biological processes is a frequently used technique to analyse the underlying mechanisms and to generally understand the behaviour of these processes in the context of a pathway, network or even the whole cell. The most common technique in this context is the usage of ordinary differential equations that describe the kinetics of the relevant processes in mechanistic detail. Here, it is usually assumed that the content of the cell is well-stirred and thus homogeneous - which is of course an over-simplification, but often worked in the past. However, many processes happen at membranes and thus not in 3D, but in 2D. The scaling of the rates of these processes poses a special problem, if volumes of compartments are changed. They will typically scale with an area, but not with the volume of the involved compartment. However, commonly, this is neglected when setting up models and/or volume scaling also sometimes automatically happens when using modelling software in the field. RESULTS: Here, we investigate generic as well as specific, realistic cases to find out, how strong the impact of the wrong scaling is for the outcome of simulations. We show that the importance of correct area scaling depends on the architecture of the reaction site and its changes upon volume alterations and it is hard to foresee, if it has a significant impact or not just by looking at the original model set-up. Moreover, scaled rates might exhibit more or less control over the behaviour of the system and therefore, accordingly, incorrect scaling will have more or less influence. CONCLUSIONS: Working with multi-compartment reactions requires a careful consideration of the correct scaling of the rates when changing the volumes of the involved compartments. The error following incorrect scaling - often done by scaling with the volume of the respective compartments can lead to significant aberrations of model behaviour.

Protein cost allocation explains metabolic strategies in Escherichia coli.

In-depth understanding of microbial growth is crucial for the development of new advances in biotechnology and for combating microbial pathogens. Condition-specific proteome expression is central to microbial physiology and growth. A multitude of processes are dependent on the protein expression, thus, whole-cell analysis of microbial metabolism using genome-scale metabolic models is an attractive toolset to investigate the behaviour of microorganisms and their communities. However, genome-scale models that incorporate macromolecular expression are still inhibitory complex: the conceptual and computational complexity of these models severely limits their potential applications. In the need for alternatives, here we revisit some of the previous attempts to create genome-scale models of metabolism and macromolecular expression to develop a novel framework for integrating protein abundance and turnover costs to conventional genome-scale models. We show that such a model of Escherichia coli successfully reproduces experimentally determined adaptations of metabolism in a growth condition-dependent manner. Moreover, the model can be used as means of investigating underutilization of the protein machinery among different growth settings. Notably, we obtained strongly improved predictions of flux distributions, considering the costs of protein translation explicitly. This finding in turn suggests protein translation being the main regulation hub for cellular growth.

miRNA target identification and prediction as a function of time in gene expression data.

The understanding of miRNA target interactions is still limited due to conflicting data and the fact that high-quality validation of targets is a time-consuming process. Faster methods like high-throughput screens and bioinformatics predictions are employed but suffer from several problems. One of these, namely the potential occurrence of downstream (i.e. secondary) effects in high-throughput screens has been only little discussed so far. However, such effects limit usage for both the identification of interactions and for the training of bioinformatics tools. In order to analyse this problem more closely, we performed time-dependent microarray screening experiments overexpressing human miR-517a-3p, and, together with published time-dependent datasets of human miR-17-5p, miR-135b and miR-124 overexpression, we analysed the dynamics of deregulated genes. We show that the number of deregulated targets increases over time, whereas seed sequence content and performance of several miRNA target prediction algorithms actually decrease over time. Bioinformatics recognition success of validated miR-17 targets was comparable to that of data gained only 12 h post-transfection. We therefore argue that the timing of microarray experiments is of critical importance for detecting direct targets with high confidence and for the usability of these data for the training of bioinformatics prediction tools.

Specifying the role of BAK1-interacting receptor-like kinase 3 in brassinosteroid signaling.

Brassinosteroids (BR) are involved in the control of several developmental processes ranging from root elongation to senescence and adaptation to environmental cues. Thus, BR perception and signaling have to be precisely regulated. One regulator is BRI1-associated kinase 1 (BAK1)-interacting receptor-like kinase 3 (BIR3). In the absence of BR, BIR3 forms complexes with BR insensitive 1 (BRI1) and BAK1. However, the biophysical and energetic requirements for complex formation in the absence of the ligand have yet to be determined. Using computational modeling, we simulated the potential complexes between the cytoplasmic domains of BAK1, BRI1 and BIR3. Our calculations and experimental data confirm the interaction of BIR3 with BAK1 and BRI1, with the BAK1 BIR3 interaction clearly favored. Furthermore, we demonstrate that BIR3 and BRI1 share the same interaction site with BAK1. This suggests a competition between BIR3 and BRI1 for binding to BAK1, which results in preferential binding of BIR3 to BAK1 in the absence of the ligand thereby preventing the active participation of BAK1 in BR signaling. Our model also suggests that BAK1 and BRI1 can interact even while BAK1 is in complex with BIR3 at an additional binding site of BAK1 that does not allow active BR signaling.

Computational systems biology of cellular processes in Arabidopsis thaliana: an overview.

Systems biology strives for gaining an understanding of biological phenomena by studying the interactions of different parts of a system and integrating the knowledge obtained into the current view of the underlying processes. This is achieved by a tight combination of quantitative experimentation and computational modeling. While there is already a large quantity of systems biology studies describing human, animal and especially microbial cell biological systems, plant biology has been lagging behind for many years. However, in the case of the model plant Arabidopsis thaliana, the steadily increasing amount of information on the levels of its genome, proteome and on a variety of its metabolic and signalling pathways is progressively enabling more researchers to construct models for cellular processes for the plant, which in turn encourages more experimental data to be generated, showing also for plant sciences how fruitful systems biology research can be. In this review, we provide an overview over some of these recent studies which use different systems biological approaches to get a better understanding of the cell biology of A. thaliana. The approaches used in these are genome-scale metabolic modeling, as well as kinetic modeling of metabolic and signalling pathways. Furthermore, we selected several cases to exemplify how the modeling approaches have led to significant advances or new perspectives in the field.

Quantitative systems pharmacology of interferon alpha administration: A multi-scale approach.

The therapeutic effect of a drug is governed by its pharmacokinetics which determine the downstream pharmacodynamic response within the cellular network. A complete understanding of the drug-effect relationship therefore requires multi-scale models which integrate the properties of the different physiological scales. Computational modelling of these individual scales has been successfully established in the past. However, coupling of the scales remains challenging, although it will provide a unique possibility of mechanistic and holistic analyses of therapeutic outcomes for varied treatment scenarios. We present a methodology to combine whole-body physiologically-based pharmacokinetic (PBPK) models with mechanistic intracellular models of signal transduction in the liver for therapeutic proteins. To this end, we developed a whole-body distribution model of IFN-α in human and a detailed intracellular model of the JAK/STAT signalling cascade in hepatocytes and coupled them at the liver of the whole-body human model. This integrated model infers the time-resolved concentration of IFN-α arriving at the liver after intravenous injection while simultaneously estimates the effect of this dose on the intracellular signalling behaviour in the liver. In our multi-scale physiologically-based pharmacokinetic/pharmacodynamic (PBPK/PD) model, receptor saturation is seen at low doses, thus giving mechanistic insights into the pharmacodynamic (PD) response. This model suggests a fourfold lower intracellular response after administration of a typical IFN-α dose to an individual as compared to the experimentally observed responses in in vitro setups. In conclusion, this work highlights clear differences between the observed in vitro and in vivo drug effects and provides important suggestions for future model-based study design.

Distribution of control in the sulfur assimilation in Arabidopsis thaliana depends on environmental conditions.

Sulfur assimilation is central to the survival of plants and has been studied under different environmental conditions. Multiple studies have been published trying to determine rate-limiting or controlling steps in this pathway. However, the picture remains inconclusive with at least two different enzymes proposed to represent such rate-limiting steps. Here, we used computational modeling to gain an integrative understanding of the distribution of control in the sulfur assimilation pathway of Arabidopsis thaliana. For this purpose, we set up a new ordinary differential equation (ODE)-based, kinetic model of sulfur assimilation encompassing all biochemical reactions directly involved in this process. We fitted the model to published experimental data and produced a model ensemble to deal with parameter uncertainties. The ensemble was validated against additional published experimental data. We used the model ensemble to subsequently analyse the control pattern and robustly identified a set of processes that share the control in this pathway under standard conditions. Interestingly, the pattern of control is dynamic and not static, that is it changes with changing environmental conditions. Therefore, while adenosine-5'-phosphosulfate reductase (APR) and sulfite reductase (SiR) share control under standard laboratory conditions, APR takes over an even more dominant role under sulfur starvation conditions.

Robustness of frequency vs. amplitude coding of calcium oscillations during changing temperatures.

Intracellular calcium oscillations have been widely studied. It is assumed that information is conveyed in the frequency, amplitude and shape of these oscillations. In particular, calcium signalling in mammalian liver cells has repeatedly been reported to display frequency coding so that an increasing amount of stimulus is translated into an increasing frequency of the oscillations. However, recently, we have shown that calcium oscillations in fish liver cells rather exhibit amplitude coding with increasing stimuli being translated into increasing amplitudes. Practical consequences of this difference are unknown so far. Here we investigated advantages and disadvantages of frequency vs. amplitude coding, in particular in environments with substantially changing temperatures (e.g. 10-20 degrees). For this purpose, we use computational modelling and a new approach to generate a calcium model exactly displaying a specific frequency and/or amplitude. We conclude that despite the advantages in flexibility that frequencies might offer for the transmission of information in the cell, amplitude coding is obviously more robust with respect to changes in environmental temperatures. This potentially explains the observed differences between two classes of organisms, one operating at constant temperatures whereas the other is not.

Advances and current challenges in calcium signaling.

Content Summary 414 I. Introduction 415 II. Ca2+ importer and exporter in plants 415 III. The Ca2+ decoding toolkit in plants 415 IV. Mechanisms of Ca2+ signal decoding 417 V. Immediate Ca2+ signaling in the regulation of ion transport 418 VI. Ca2+ signal integration into long-term ABA responses 419 VII Integration of Ca2+ and hormone signaling through dynamic complex modulation of the CCaMK/CYCLOPS complex 420 VIII Ca2+ signaling in mitochondria and chloroplasts 422 IX A view beyond recent advances in Ca2+ imaging 423 X Modeling approaches in Ca2+ signaling 424 XI Conclusions: Ca2+ signaling a still young blooming field of plant research 424 Acknowledgements 425 ORCID 425 References 425 SUMMARY: Temporally and spatially defined changes in Ca2+ concentration in distinct compartments of cells represent a universal information code in plants. Recently, it has become evident that Ca2+ signals not only govern intracellular regulation but also appear to contribute to long distance or even organismic signal propagation and physiological response regulation. Ca2+ signals are shaped by an intimate interplay of channels and transporters, and during past years important contributing individual components have been identified and characterized. Ca2+ signals are translated by an elaborate toolkit of Ca2+ -binding proteins, many of which function as Ca2+ sensors, into defined downstream responses. Intriguing progress has been achieved in identifying specific modules that interconnect Ca2+ decoding proteins and protein kinases with downstream target effectors, and in characterizing molecular details of these processes. In this review, we reflect on recent major advances in our understanding of Ca2+ signaling and cover emerging concepts and existing open questions that should be informative also for scientists that are currently entering this field of ever-increasing breath and impact.

A20/TNFAIP3 Discriminates Tumor Necrosis Factor (TNF)-Induced NF-κB from JNK Pathway Activation in Hepatocytes.

In the liver tumor necrosis factor (TNF)-induced signaling critically regulates the immune response of non-parenchymal cells as well as proliferation and apoptosis of hepatocytes via activation of the NF-κB and JNK pathways. Especially, the induction of negative feedback regulators, such as IκBα and A20 is responsible for the dynamic and time-restricted response of these important pathways. However, the precise mechanisms responsible for different TNF-induced phenotypes under physiological stimulation conditions are not completely understood so far. In addition, it is not known if varying TNF concentrations may differentially affect the desensitization properties of both pathways. By using computational modeling, we first showed that TNF-induced activation and downstream signaling is qualitatively comparable between primary mouse hepatocytes and immortalized hepatocellular carcinoma (HCC) cells. In order to define physiologically relevant TNF levels, which allow for an adjustable and dynamic NF-κB/JNK pathway response in parenchymal liver cells, a range of cytokine concentrations was defined that led to gradual pathway responses in HCC cells (1-5 ng/ml). Repeated stimulations with low (1 ng/ml), medium (2.5 ng/ml) and high (5 ng/ml) TNF amounts demonstrated that JNK signaling was still active at cytokine concentrations, which led to dampened NF-κB signaling illustrating differential pathway responsiveness depending on TNF input dynamics. SiRNA-mediated inhibition of the negative feedback regulator A20 (syn. TNFAIP3) or its overexpression did not significantly affect the NF-κB response. In contrast, A20 silencing increased the JNK response, while its overexpression dampened JNK phosphorylation. In addition, the A20 knockdown sensitized hepatocellular cells to TNF-induced cleavage and activity of the effector caspase-3. In conclusion, a mathematical model-based approach shows that the TNF-induced pathway responses are qualitatively comparable in primary and immortalized mouse hepatocytes. The cytokine amount defines the pathway responsiveness under repeated treatment conditions with NF-κB signaling being dampened 'earlier' than JNK. A20 appears to be the molecular switch discriminating between NF-κB and JNK signaling when stimulating with varying physiological cytokine concentrations.

COPASI and its applications in biotechnology.

COPASI is software used for the creation, modification, simulation and computational analysis of kinetic models in various fields. It is open-source, available for all major platforms and provides a user-friendly graphical user interface, but is also controllable via the command line and scripting languages. These are likely reasons for its wide acceptance. We begin this review with a short introduction describing the general approaches and techniques used in computational modeling in the biosciences. Next we introduce the COPASI package, and its capabilities, before looking at typical applications of COPASI in biotechnology.

A new efficient approach to fit stochastic models on the basis of high-throughput experimental data using a model of IRF7 gene expression as case study.

BACKGROUND: Mathematical models are used to gain an integrative understanding of biochemical processes and networks. Commonly the models are based on deterministic ordinary differential equations. When molecular counts are low, stochastic formalisms like Monte Carlo simulations are more appropriate and well established. However, compared to the wealth of computational methods used to fit and analyze deterministic models, there is only little available to quantify the exactness of the fit of stochastic models compared to experimental data or to analyze different aspects of the modeling results. RESULTS: Here, we developed a method to fit stochastic simulations to experimental high-throughput data, meaning data that exhibits distributions. The method uses a comparison of the probability density functions that are computed based on Monte Carlo simulations and the experimental data. Multiple parameter values are iteratively evaluated using optimization routines. The method improves its performance by selecting parameters values after comparing the similitude between the deterministic stability of the system and the modes in the experimental data distribution. As a case study we fitted a model of the IRF7 gene expression circuit to time-course experimental data obtained by flow cytometry. IRF7 shows bimodal dynamics upon IFN stimulation. This dynamics occurs due to the switching between active and basal states of the IRF7 promoter. However, the exact molecular mechanisms responsible for the bimodality of IRF7 is not fully understood. CONCLUSIONS: Our results allow us to conclude that the activation of the IRF7 promoter by the combination of IRF7 and ISGF3 is sufficient to explain the observed bimodal dynamics.