Apoptosis is regulated by several signaling pathways which are extensively linked by crosstalks. Boolean or logical modeling has become a promising approach to capture the qualitative behavior of such complex networks. Here we built a large-scale literature-based Boolean model of the central intrinsic and extrinsic apoptosis pathways as well as pathways connected with them. The model responds to several external stimuli such as Fas ligand, TNF-α, UV-B irradiation, interleukin-1β and insulin. Timescales and multi-value node logic were used and turned out to be indispensable to reproduce the behavior of the apoptotic network. The coherence of the model was experimentally validated. Thereby an UV-B dose-effect is shown for the first time in mouse hepatocytes. Analysis of the model revealed a tight regulation emerging from high connectivity and spanning crosstalks and a particular importance of feedback loops. An unexpected feedback from Smac release to RIP could further increase complex II formation. The introduced Boolean model provides a comprehensive and coherent description of the apoptosis network behavior. It gives new insights into the complex interplay of pro- and antiapoptotic factors and can be easily expanded to other signaling pathways.

A mathematical model representing metabolite interconversions in the central carbohydrate metabolism of Arabidopsis (Arabidopsis thaliana) was developed to simulate the diurnal dynamics of primary carbon metabolism in a photosynthetically active plant leaf. The model groups enzymatic steps of central carbohydrate metabolism into blocks of interconverting reactions that link easily measurable quantities like CO 2 exchange and quasi-steady-state levels of soluble sugars and starch. When metabolite levels that fluctuate over diurnal cycles are used as a basic condition for simulation, turnover rates for the interconverting reactions can be calculated that approximate measured metabolite dynamics and yield kinetic parameters of interconverting reactions. We used experimental data for Arabidopsis wild-type plants, accession Columbia, and a mutant defective in vacuolar invertase, AtbFruct4, as input data. Reducing invertase activity to mutant levels in the wild-type model led to a correct prediction of increased sucrose levels. However, additional changes were needed to correctly simulate levels of hexoses and sugar phosphates, indicating that invertase knockout causes subsequent changes in other enzymatic parameters. Reduction of invertase activity caused a decline in photosynthesis and export of reduced carbon to associated metabolic pathways and sink organs (e.g. roots), which is in agreement with the reported contribution of vacuolar invertase to sink strength. According to model parameters, there is a role for invertase in leaves, where futile cycling of sucrose appears to have a buffering effect on the pools of sucrose, hexoses, and sugar phosphates. Our data demonstrate that modeling complex metabolic pathways is a useful tool to study the significance of single enzyme activities in complex, nonintuitive networks.

The dynamics of biological reaction networks are strongly constrained by thermodynamics. An holistic understanding of their behavior and regulation requires mathematical models that observe these constraints. However, kinetic models may easily violate the constraints imposed by the principle of detailed balance, if no special care is taken. Detailed balance demands that in thermodynamic equilibrium all fluxes vanish. We introduce a thermodynamic-kinetic modeling (TKM) formalism that adapts the concepts of potentials and forces from irreversible thermodynamics to kinetic modeling. In the proposed formalism, the thermokinetic potential of a compound is proportional to its concentration. The proportionality factor is a compound-specific parameter called capacity. The thermokinetic force of a reaction is a function of the potentials. Every reaction has a resistance that is the ratio of thermokinetic force and reaction rate. For mass-action type kinetics, the resistances are constant. Since it relies on the thermodynamic concept of potentials and forces, the TKM formalism structurally observes detailed balance for all values of capacities and resistances. Thus, it provides an easy way to formulate physically feasible, kinetic models of biological reaction networks. The TKM formalism is useful for modeling large biological networks that are subject to many detailed balance relations.