Hard-wired heterogeneity in blood stem cells revealed using a dynamic regulatory network model.

MOTIVATION: Combinatorial interactions of transcription factors with cis-regulatory elements control the dynamic progression through successive cellular states and thus underpin all metazoan development. The construction of network models of cis-regulatory elements, therefore, has the potential to generate fundamental insights into cellular fate and differentiation. Haematopoiesis has long served as a model system to study mammalian differentiation, yet modelling based on experimentally informed cis-regulatory interactions has so far been restricted to pairs of interacting factors. Here, we have generated a Boolean network model based on detailed cis-regulatory functional data connecting 11 haematopoietic stem/progenitor cell (HSPC) regulator genes. RESULTS: Despite its apparent simplicity, the model exhibits surprisingly complex behaviour that we charted using strongly connected components and shortest-path analysis in its Boolean state space. This analysis of our model predicts that HSPCs display heterogeneous expression patterns and possess many intermediate states that can act as 'stepping stones' for the HSPC to achieve a final differentiated state. Importantly, an external perturbation or 'trigger' is required to exit the stem cell state, with distinct triggers characterizing maturation into the various different lineages. By focusing on intermediate states occurring during erythrocyte differentiation, from our model we predicted a novel negative regulation of Fli1 by Gata1, which we confirmed experimentally thus validating our model. In conclusion, we demonstrate that an advanced mammalian regulatory network model based on experimentally validated cis-regulatory interactions has allowed us to make novel, experimentally testable hypotheses about transcriptional mechanisms that control differentiation of mammalian stem cells. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Transplantation tolerance induced by regulatory T cells: in vivo mechanisms and sites of action.

The mechanisms by which CD4(+)CD25(+)Foxp3(+) T (Treg) cells regulate effector T cells in a transplantation setting and their in vivo homeostasis still remain to be clarified. Using a mouse adoptive transfer model, we analyzed the in vivo expansion, trafficking, and effector function of alloreactive T cells and donor-specific Treg cells, in response to a full-thickness skin allograft. Fluorescent-labeled CD4(+)CD25(-) and antigen-specific Treg cells were transferred alone or co-injected into syngeneic BALB/c-Nude recipients transplanted with skins from (C57BL/6 x BALB/c) F1 donors. Treg cells divided in vivo, migrated and accumulated in the allograft draining lymph nodes as well as within the graft. The co-transfer of Treg cells did not modify the early activation and homing of CD4(+)CD25(-) T cells in secondary lymphoid organs. However, in the presence of Treg cells, alloreactive CD4(+)CD25(-) T cells produced significantly less IFN-gamma and were present in reduced numbers in the secondary lymphoid organs. Furthermore, time-course studies showed that Treg cells were recruited into the allograft at a very early stage after transplantation and effectively prevented the infiltration of effector T cells. In conclusion, suppression of rejection requires the early recruitment to the site of antigenic challenge of donor-specific Treg cells, which then mainly regulate the effector arm of T cell alloresponses.

Dynamical modeling and analysis of large cellular regulatory networks.

The dynamical analysis of large biological regulatory networks requires the development of scalable methods for mathematical modeling. Following the approach initially introduced by Thomas, we formalize the interactions between the components of a network in terms of discrete variables, functions, and parameters. Model simulations result in directed graphs, called state transition graphs. We are particularly interested in reachability properties and asymptotic behaviors, which correspond to terminal strongly connected components (or "attractors") in the state transition graph. A well-known problem is the exponential increase of the size of state transition graphs with the number of network components, in particular when using the biologically realistic asynchronous updating assumption. To address this problem, we have developed several complementary methods enabling the analysis of the behavior of large and complex logical models: (i) the definition of transition priority classes to simplify the dynamics; (ii) a model reduction method preserving essential dynamical properties, (iii) a novel algorithm to compact state transition graphs and directly generate compressed representations, emphasizing relevant transient and asymptotic dynamical properties. The power of an approach combining these different methods is demonstrated by applying them to a recent multilevel logical model for the network controlling CD4+ T helper cell response to antigen presentation and to a dozen cytokines. This model accounts for the differentiation of canonical Th1 and Th2 lymphocytes, as well as of inflammatory Th17 and regulatory T cells, along with many hybrid subtypes. All these methods have been implemented into the software GINsim, which enables the definition, the analysis, and the simulation of logical regulatory graphs.

Time-Course In Vivo Trafficking Pattern and Effector Function of Donor-Specific Regulatory T Cells in Response to an Allograft.

Background: Experimental data have suggested that adoptive transfer of CD4+CD25+Foxp3+ regulatory T cells (Tregs), capable of controlling immune responses to specifi c auto- or alloantigens, could be used as a therapeutic strategy to promote specifi c tolerance in T-cell mediated diseases and in organ transplantation (Tx). However, before advocating the application of immunotherapy with Tregs in Tx, we need to improve our understanding of their in vivo homeostasis, traffi cking pattern and effector function in response to alloantigens. Methods : Donor-antigen specifi c murine Tregs were generated and characterized in vitro following our described protocols. Using an adoptive transfer and skin allotransplantation model, we have analyzed the in vivo expansion and homing of fl uorescent-labeled effector T cells (Teff) and Tregs, at different time-points after Tx, using fl ow-cytometry as well as fl uorescence microscopy techniques. Results: Tregs expressed CD62L, CCR7 and CD103 allowing their homing into lymphoid and non-lymphoid tissues (gut, skin) after intravenous injection. While hyporesponsive to TCR stimulation in vitro, transferred Tregs survived, migrated to secondary lymphoid organs and preferentially expanded within the allograft draining lymph nodes. Furthermore, Foxp3+ cells could be detected inside the allograft as early as day 3-5 after Tx. At a much later time-point (day 60 after Tx), graft-infi ltrating Foxp3+ cells were also detectable in tolerant recipients. When transferred alone, CD4+CD25- Teff cells expanded within secondary lymphoid organs and infi ltrated the allograft by day 3-5 after Tx. The co-transfer of Tregs limited the expansion of alloreactive Teff cells as well as their recruitment into the allograft. The promotion of graft survival observed in the presence of Tregs was in part mediated by the inhibition of the production of effector cytokines by CD4+CD25- T cells. Conclusion: Taken together, our results suggest that the suppression of allograft rejection and the induction of Tx tolerance are in part dependant on the alloantigendriven homing and expansion of Tregs. Thus, the appropriate localization of Tregs may be critical for their suppressive function in vivo.

CD4+CD25-mTGFbeta+ T cells induced by nasal application of ovalbumin transfer tolerance in a therapeutic model of asthma.

Background: Intranasal administration of high amount of allergen was shown to induce tolerance and to reverse the allergic phenotype. However, mechanisms of tolerance induction via the mucosal route are still unclear. Objectives: To characterize the therapeutic effects of intranasal application of ovalbumin (OVA) in a mouse model of bronchial inflammation as well as the cellular and molecular mechanisms leading to protection upon re-exposure to allergen. Methods: After induction of bronchial inflammation, mice were treated intranasally with OVA and re-exposed to OVA aerosols 10 days later. Bronchoalveolar lavage fluid (BALF), T cell proliferation and cytokine secretion were examined. The respective role of CD4(+)CD25(+) and CD4(+)CD25(-) T cells in the induction of tolerance was analysed. Results: Intranasal treatment with OVA drastically reduced inflammatory cell recruitment into BALF and bronchial hyperresponsiveness upon re-exposure to allergen. Both OVA- specific-proliferation of T cells, T(h)1 and T(h)2 cytokine production from lung and bronchial lymph nodes were inhibited. Transfer of CD4(+)CD25(-) T cells, which strongly expressed membrane-bound transforming growth factor beta (mTGF beta), from tolerized mice protected asthmatic recipient mice from subsequent aerosol challenges. The presence of CD4(+)CD25(+)(Foxp3(+)) T cells during the process of tolerization was indispensable to CD4(+)CD25(-) T cells to acquire regulatory properties. Whereas the presence of IL-10 appeared dispensable in this model, the suppression of CD4(+)CD25(-)mTGF beta(+) T cells in transfer experiments significantly impaired the down-regulation of airways inflammation. Conclusion: Nasal application of OVA in established asthma led to the induction of CD4(+)CD25(-)mTGF beta(+) T cells with regulatory properties, able to confer protection upon allergen re-exposure.

The de facto independence of regulatory agencies and its consequences for policy making and regulatory outcomes

This dissertation focuses on the practice of regulatory governance, throughout the study of the functioning of formally independent regulatory agencies (IRAs), with special attention to their de facto independence. The research goals are grounded on a "neo-positivist" (or "reconstructed positivist") position (Hawkesworth 1992; Radaelli 2000b; Sabatier 2000). This perspective starts from the ontological assumption that even if subjective perceptions are constitutive elements of political phenomena, a real world exists beyond any social construction and can, however imperfectly, become the object of scientific inquiry. Epistemologically, it follows that hypothetical-deductive theories with explanatory aims can be tested by employing a proper methodology and set of analytical techniques. It is thus possible to make scientific inferences and general conclusions to a certain extent, according to a Bayesian conception of knowledge, in order to update the prior scientific beliefs in the truth of the related hypotheses (Howson 1998), while acknowledging the fact that the conditions of truth are at least partially subjective and historically determined (Foucault 1988; Kuhn 1970). At the same time, a sceptical position is adopted towards the supposed disjunction between facts and values and the possibility of discovering abstract universal laws in social science. It has been observed that the current version of capitalism corresponds to the golden age of regulation, and that since the 1980s no government activity in OECD countries has grown faster than regulatory functions (Jacobs 1999). Following an apparent paradox, the ongoing dynamics of liberalisation, privatisation, decartelisation, internationalisation, and regional integration hardly led to the crumbling of the state, but instead promoted a wave of regulatory growth in the face of new risks and new opportunities (Vogel 1996). Accordingly, a new order of regulatory capitalism is rising, implying a new division of labour between state and society and entailing the expansion and intensification of regulation (Levi-Faur 2005). The previous order, relying on public ownership and public intervention and/or on sectoral self-regulation by private actors, is being replaced by a more formalised, expert-based, open, and independently regulated model of governance. Independent regulation agencies (IRAs), that is, formally independent administrative agencies with regulatory powers that benefit from public authority delegated from political decision makers, represent the main institutional feature of regulatory governance (Gilardi 2008). IRAs constitute a relatively new technology of regulation in western Europe, at least for certain domains, but they are increasingly widespread across countries and sectors. For instance, independent regulators have been set up for regulating very diverse issues, such as general competition, banking and finance, telecommunications, civil aviation, railway services, food safety, the pharmaceutical industry, electricity, environmental protection, and personal data privacy. Two attributes of IRAs deserve a special mention. On the one hand, they are formally separated from democratic institutions and elected politicians, thus raising normative and empirical concerns about their accountability and legitimacy. On the other hand, some hard questions about their role as political actors are still unaddressed, though, together with regulatory competencies, IRAs often accumulate executive, (quasi-)legislative, and adjudicatory functions, as well as about their performance.

Synchronous versus asynchronous modeling of gene regulatory networks.

MOTIVATION: In silico modeling of gene regulatory networks has gained some momentum recently due to increased interest in analyzing the dynamics of biological systems. This has been further facilitated by the increasing availability of experimental data on gene-gene, protein-protein and gene-protein interactions. The two dynamical properties that are often experimentally testable are perturbations and stable steady states. Although a lot of work has been done on the identification of steady states, not much work has been reported on in silico modeling of cellular differentiation processes. RESULTS: In this manuscript, we provide algorithms based on reduced ordered binary decision diagrams (ROBDDs) for Boolean modeling of gene regulatory networks. Algorithms for synchronous and asynchronous transition models have been proposed and their corresponding computational properties have been analyzed. These algorithms allow users to compute cyclic attractors of large networks that are currently not feasible using existing software. Hereby we provide a framework to analyze the effect of multiple gene perturbation protocols, and their effect on cell differentiation processes. These algorithms were validated on the T-helper model showing the correct steady state identification and Th1-Th2 cellular differentiation process. AVAILABILITY: The software binaries for Windows and Linux platforms can be downloaded from http://si2.epfl.ch/~garg/genysis.html.

New insights into the regulatory mechanisms of the NALP3 inflammasome

Résumé : Les vertébrés ont recours au système immunitaire inné et adaptatif pour combattre les pathogènes. La découverte des récepteurs Toll, il y a dix ans, a fortement augmenté l'intérêt porté à l'immunité innée. Depuis lors, des récepteurs intracellulaires tels que les membres de la famille RIG-like helicase (RLHs) et NOD-like receptor (NLRs) ont été décrits pour leur rôle dans la détection des pathogènes. L'interleukine-1 beta (IL-1β) est une cytokine pro-inflammatoire qui est synthétisée sous forme de précurseur, la proIL-1β. La proIL-1β requiert d'être clivée par la caspase-1 pour devenir active. La caspase-1 est elle-même activée par un complexe appelé inflammasome qui peut être formé par divers membres de la famille NLR. Plusieurs inflammasomes ont été décrits tels que le NALP3 inflammasome ou l'IPAF inflammasome. Dans cette étude nous avons identifié la co-chaperone SGT1 et la chaperone HSP90 comme partenaires d'interaction de NALP3. Ces deux protéines sont bien connues chez les plantes pour leurs rôles dans la régulation des gènes de résistance (gène R) qui sont structurellement apparentés à la famille NLR. Nous avons pu montrer que SGT1 et HSP90 jouent un rôle similaire dans la régulation de NALP3 et des protéines R. En effet, nous avons démontré que les deux protéines sont nécessaires pour l'activité du NALP3 inflammasome. De plus, la HSP90 est également requise pour la stabilité de NALP3. En se basant sur ces observations, nous avons proposé un modèle dans lequel SGT1 et HSP90 maintiennent NALP3 inactif mais prêt à percevoir un ligand activateur qui initierait la cascade inflammatoire. Nous avons également montré une interaction entre SGT1 et HSP90 avec plusieurs NLRs. Cette observation suggère qu'un mécanisme similaire pourrait être impliqué dans la régulation des membres de la famille des NLRs. Ces dernières années, plusieurs PAMPs mais également des DAMPs ont été identifiés comme activateurs du NALP3 inflammasome. Dans la seconde partie de cette étude, nous avons identifié la réponse au stress du réticulum endoplasmique (RE) comme nouvel activateur du NALP3 inflammasome. Cette réponse est initiée lors de l'accumulation dans le réticulum endoplasmique de protéines ayant une mauvaise conformation ce qui conduit, en autre, à l'arrêt de la synthèse de nouvelles protéines ainsi qu'une augmentation de la dégradation des protéines. Les mécanismes par lesquels la réponse du réticulum endoplasmique induit l'activation du NALP3 inflammasome doivent encore être déterminés. Summary : Vertebrates rely on the adaptive and the innate immune systems to fight pathogens. Awarness of the importance of the innate system increased with the identification of Toll-like receptors a decade ago. Since then, intracellular receptors such as the RIG-like helicase (RLH) and the NOD-like receptor (NLR) families have been described for their role in the recognition of microbes. Interleukin- 1ß (IL-1ß) is a key mediator of inflammation. This proinflammatory cytokine is synthesised as an inactive precursor that requires processing by caspase-1 to become active. Caspase-1 is, itself, activated in a complex termed the inflammasome that can be formed by members of the NLR family. Various inflammasome complexes have been described such as the IPAF and the NALP3 inflammasome. In this study, we have identified the co-chaperone SGT1 and the chaperone HSP90 as interacting partners of NALP3. SGT1 and HSP90 are both known for their role in the activity of plant resistance proteins (R proteins) which are structurally related to the NLR family. We have shown that HSP90 and SGT1 play a similar role in the regulation of NALP3 and in the regulation of plant R proteins. Indeed, we demonstrated that both HSP90 and SGT1 are essential for the activity of the NALP3 inflammasome complex. In addition, HSP90 is required for the stability of NALP3. Based on these observations, we have proposed a model in which SGT1 and HSP90 maintain NALP3 in an inactive but signaling-competent state, ready to receive an activating ligand that induces the inflammatory cascade. An interaction between several NLR members, SGTI and HSP90 was also shown, suggesting that similar mechanisms could be involved in the regulation of other NLRs. Several pathogen-associated molecular patterns (PAMPs) but also danger associated molecular patterns (DAMPs) have been identified as NALP3 activators. In the second part of this study, we have identified the ER stress response as a new NALP3 activator. The ER stress response is activated upon the accumulation of unfolded protein in the endoplasmic reticulum and results in a block in protein synthesis and increased protein degradation. The mechanisms of ER stress-mediated NALP3 activation remain to be determined.

Effects of immunosuppressive drugs on T helper cells subsets and potential to promote tolerance in a stringent experimental transplantation model.

To achieve the goal of sustained donor-specifi c transplantation (Tx) tolerance, research efforts are now focusing on therapies based on specifi c cell subsets with regulatory properties. We and others have previously highlighted the therapeutic potential of naturally occurring CD4+CD25+Foxp3+ regulatory T cells (nTreg) in promoting long-term graft acceptance. Using more stringent experimental Tx models, we were however confronted to limitations. Indeed, while the transfer of antigenspecifi c nTreg promoted long-term MHC-mismatched skin allograft acceptance in lymphopenic mice in the absence of any immunosuppressive drug, allograft survival was only slightly prolonged when nTreg were transferred alone into non-lymphopenic mice. This suggested that in more stringent conditions, adjuvant therapies may be needed to effectively control alloreactive T cells (Teff). Whether and how the expansion of the Treg pool could be best combined with current immunosuppressive regimens in clinical settings remains to be defi ned. In this study, we have used in vitro assays and an in vivo skin Tx model to investigate the effects of various immunosuppressive drugs on the survival, proliferation and effector function of Teff and nTreg in response to alloantigens. Teff proliferation was inhibited in a dose-dependent manner by rapamycin and cyclosporine A, while anti-CD154 mAb only marginally affected Teff survival, proliferation and effector fucntion in vitro. Rapamycin promoted apoptosis of Teff as compared to nTreg that were more resistant in the presence of IL-2. In vivo, the transfer and/or expansion of Treg could be advantageously combined with rapamycin and anti-CD154 mAb treatment to signifi cantly prolong MHC-mismatched skin allografts survival in non-lymphopenic recipients. Taken together our data indicate that immunosuppressive drugs differentially target T-cell subsets and that some regimens could promote Treg expansion while controlling the Teff pool in response to alloantigens.

Strength and limitations of regulatory T Cells for immunotherapy in transplantation.

In many experimental models, CD4+CD25+Foxp3+ regulatory T cells (nTreg) have been identifi ed as key players in promoting peripheral transplantation (Tx) tolerance. We have been focusing on therapies based on antigen-specifi c nTreg that can control effector T cells (Teff) and prevent allograft rejection. The use of nTreg in immunotherapeutic protocols for solid organ Tx is however limited by their overall low numbers as well as the low precursor frequency of alloantigen cross-reactive nTreg expected to be found in a normal individual. Moreover, although we previously described robust protocols to generate and expand antigen-specifi c nTreg in vitro, the process requires careful selection of highly pure nTreg and cumbersome ex-vivo manipulations, rendering this strategy not easily applicable in clinical solid organ Tx. In this study, we aimed to expand Treg directly in vivo and determine their suppressive function, effi cacy and stability in promoting donor-specifi c tolerance in a stringent murine Tx model. Our data suggest that IL-2-based therapies lead to a signifi cant increase of Treg in vivo. The expanded Treg suppressed Teff proliferation (albeit slightly less effi ciently than nTreg isolated from control mice) and allowed prolonged graft survival of major MHC-mismatched skin grafts in wild-type non-lymphopenic recipients. The expanded Treg alone were however not suffi cient to induce tolerance in stringent experimental conditions. Rapamycin reduced the frequency of Teff but did not impede expansion of Treg. Pro-infl ammatory stimuli hindered the expansion of Treg and resulted in an increase in the frequency of CD4+IFN-γ+ and CD4+IL17+ T cells. We propose that IL-2-based treatments would be an effi cient method for expanding functional Treg in vivo without affecting other immune cell populations, thereby favorably shifting the pool of alloreactive T cells towards regulation in response to an allograft. However, we also highlight some potential limitations of Treg expansion such as concomitant infl ammatory events.

Regulatory T cell defects are associated with autoimmunity in Wiskott-Aldrich Syndrome protein deficient mice

Abstract : The Wiskott-Aldrich Syndrome (WAS) is an X-linked recessive human primary immunodeficiency. It is caused by mutations in the gene encoding the hermatopoietic specific regulator of the actin cytoskeleton Wiskott-Aldrich Syndrome Protein (WASP). Importantly, a majority of affected patients develop autoimmunity including an inflammatory bowel disease (IBD)-like disease. WASP deficient mice share many similarities with the human WAS. One of these similarities is the spontaneous development of colitis. I have focused my dissertation studies on the pathogenesis of colitis in WASP deficient mice. Prior work from our laboratory had shown that lymphocytes were required and that CD4+ T cells sufficient for colitis development. This colitis was associated with a predominant Th2-cytokine skewing. I have contributed in exploring whether the Th2 cytokine IL-4 plays a role in disease maintenance. Using two approaches to neutralize IL-4, we found that this cytokine plays a role in disease maintenance. Natural CD4*CD25*Foxp3* regulatory T cells (nTreg cells) have been implicated in the pathogenesis of several autoimmune disorders. We found that WASP deficient mice have reduced nTreg cell numbers in peripheral lymphoid organs. This was associated with functional defects in suppressing T cell proliferation and preventing colitis induced by transfer of naïve T cells into SCID recipient, which lack lymphocytes. WASP deficiency affected homing of nTreg cells to lymphoid compartments, IL-2-mediated activation and secretion of the immunomodulatory cytokine IL-10. Finally, we could prevent colitis onset via adoptive transfer of WT nTreg cells prior to colitis development. This suggests that nTreg cells dysfunction is one of the mechanisms underlying colitis development in WASP deficient mice. Future directions will aim at deciphering the role of other immune cell types, the bacterial flora, and various cytokines in colitis development in this murine model of colitis. In addition, we believe that colitis in WASP deficient mice could serve as a useful tool to evaluate nTreg cells manipulation as novel therapeutic approach for IBD.

Modeling stochasticity and robustness in gene regulatory networks.

MOTIVATION: Understanding gene regulation in biological processes and modeling the robustness of underlying regulatory networks is an important problem that is currently being addressed by computational systems biologists. Lately, there has been a renewed interest in Boolean modeling techniques for gene regulatory networks (GRNs). However, due to their deterministic nature, it is often difficult to identify whether these modeling approaches are robust to the addition of stochastic noise that is widespread in gene regulatory processes. Stochasticity in Boolean models of GRNs has been addressed relatively sparingly in the past, mainly by flipping the expression of genes between different expression levels with a predefined probability. This stochasticity in nodes (SIN) model leads to over representation of noise in GRNs and hence non-correspondence with biological observations. RESULTS: In this article, we introduce the stochasticity in functions (SIF) model for simulating stochasticity in Boolean models of GRNs. By providing biological motivation behind the use of the SIF model and applying it to the T-helper and T-cell activation networks, we show that the SIF model provides more biologically robust results than the existing SIN model of stochasticity in GRNs. AVAILABILITY: Algorithms are made available under our Boolean modeling toolbox, GenYsis. The software binaries can be downloaded from http://si2.epfl.ch/ approximately garg/genysis.html.

Lifelong dynamics of human CD4+CD25+ regulatory T cells: insights from in vivo data and mathematical modeling.

Despite their limited proliferation capacity, regulatory T cells (T(regs)) constitute a population maintained over the entire lifetime of a human organism. The means by which T(regs) sustain a stable pool in vivo are controversial. Using a mathematical model, we address this issue by evaluating several biological scenarios of the origins and the proliferation capacity of two subsets of T(regs): precursor CD4(+)CD25(+)CD45RO(-) and mature CD4(+)CD25(+)CD45RO(+) cells. The lifelong dynamics of T(regs) are described by a set of ordinary differential equations, driven by a stochastic process representing the major immune reactions involving these cells. The model dynamics are validated using data from human donors of different ages. Analysis of the data led to the identification of two properties of the dynamics: (1) the equilibrium in the CD4(+)CD25(+)FoxP3(+)T(regs) population is maintained over both precursor and mature T(regs) pools together, and (2) the ratio between precursor and mature T(regs) is inverted in the early years of adulthood. Then, using the model, we identified three biologically relevant scenarios that have the above properties: (1) the unique source of mature T(regs) is the antigen-driven differentiation of precursors that acquire the mature profile in the periphery and the proliferation of T(regs) is essential for the development and the maintenance of the pool; there exist other sources of mature T(regs), such as (2) a homeostatic density-dependent regulation or (3) thymus- or effector-derived T(regs), and in both cases, antigen-induced proliferation is not necessary for the development of a stable pool of T(regs). This is the first time that a mathematical model built to describe the in vivo dynamics of regulatory T cells is validated using human data. The application of this model provides an invaluable tool in estimating the amount of regulatory T cells as a function of time in the blood of patients that received a solid organ transplant or are suffering from an autoimmune disease.

Additive functions in boolean models of gene regulatory network modules.

Gene-on-gene regulations are key components of every living organism. Dynamical abstract models of genetic regulatory networks help explain the genome's evolvability and robustness. These properties can be attributed to the structural topology of the graph formed by genes, as vertices, and regulatory interactions, as edges. Moreover, the actual gene interaction of each gene is believed to play a key role in the stability of the structure. With advances in biology, some effort was deployed to develop update functions in Boolean models that include recent knowledge. We combine real-life gene interaction networks with novel update functions in a Boolean model. We use two sub-networks of biological organisms, the yeast cell-cycle and the mouse embryonic stem cell, as topological support for our system. On these structures, we substitute the original random update functions by a novel threshold-based dynamic function in which the promoting and repressing effect of each interaction is considered. We use a third real-life regulatory network, along with its inferred Boolean update functions to validate the proposed update function. Results of this validation hint to increased biological plausibility of the threshold-based function. To investigate the dynamical behavior of this new model, we visualized the phase transition between order and chaos into the critical regime using Derrida plots. We complement the qualitative nature of Derrida plots with an alternative measure, the criticality distance, that also allows to discriminate between regimes in a quantitative way. Simulation on both real-life genetic regulatory networks show that there exists a set of parameters that allows the systems to operate in the critical region. This new model includes experimentally derived biological information and recent discoveries, which makes it potentially useful to guide experimental research. The update function confers additional realism to the model, while reducing the complexity and solution space, thus making it easier to investigate.