A method for the generation of standardized qualitative dynamical systems of regulatory networks
- Luis Mendoza^{1}Email author and
- Ioannis Xenarios^{1}
DOI: 10.1186/1742-4682-3-13
© Mendoza and Xenarios; licensee BioMed Central Ltd. 2006
Received: 12 December 2005
Accepted: 16 March 2006
Published: 16 March 2006
Abstract
Background
Modeling of molecular networks is necessary to understand their dynamical properties. While a wealth of information on molecular connectivity is available, there are still relatively few data regarding the precise stoichiometry and kinetics of the biochemical reactions underlying most molecular networks. This imbalance has limited the development of dynamical models of biological networks to a small number of well-characterized systems. To overcome this problem, we wanted to develop a methodology that would systematically create dynamical models of regulatory networks where the flow of information is known but the biochemical reactions are not. There are already diverse methodologies for modeling regulatory networks, but we aimed to create a method that could be completely standardized, i.e. independent of the network under study, so as to use it systematically.
Results
We developed a set of equations that can be used to translate the graph of any regulatory network into a continuous dynamical system. Furthermore, it is also possible to locate its stable steady states. The method is based on the construction of two dynamical systems for a given network, one discrete and one continuous. The stable steady states of the discrete system can be found analytically, so they are used to locate the stable steady states of the continuous system numerically. To provide an example of the applicability of the method, we used it to model the regulatory network controlling T helper cell differentiation.
Conclusion
The proposed equations have a form that permit any regulatory network to be translated into a continuous dynamical system, and also find its steady stable states. We showed that by applying the method to the T helper regulatory network it is possible to find its known states of activation, which correspond the molecular profiles observed in the precursor and effector cell types.
Background
The increasing use of high throughput technologies in different areas of biology has generated vast amounts of molecular data. This has, in turn, fueled the drive to incorporate such data into pathways and networks of interactions, so as to provide a context within which molecules operate. As a result, a wealth of connectivity information is available for multiple biological systems, and this has been used to understand some global properties of biological networks, including connectivity distribution [1], recurring motifs [2] and modularity [3]. Such information, while valuable, provides only a static snapshot of a network. For a better understanding of the functionality of a given network it is important to study its dynamical properties. The consideration of dynamics allows us to answer questions related to the number, nature and stability of the possible patterns of activation, the contribution of individual molecules or interactions to establishing such patterns, and the possibility of simulating the effects of loss- or gain-of-function mutations, for example.
Mathematical modeling of metabolic networks requires specification of the biochemical reactions involved. Each reaction has to incorporate the appropriate stoichiometric coefficients to account for the principle of mass conservation. This characteristic simplifies modeling, because it implies that at equilibrium every node of the metabolic network has a total mass flux of zero [4, 5]. There are cases, however, where the underlying biochemical reactions are not known for large parts of a pathway, but the direction of the flow of information is known, which is the case for so-called regulatory networks (see for example [6, 7]). In these cases, the directionality of signaling is sufficient for developing mathematical models of how the patterns of activation and inhibition determine the state of activation of the network (for a review, see [8]).
When cells receive external stimuli such as hormones, mechanical forces, changes in osmolarity, membrane potential etc., there is an internal response in the form of multiple intracellular signals that may be buffered or may eventually be integrated to trigger a global cellular response, such as growth, cell division, differentiation, apoptosis, secretion etc. Modeling the underlying molecular networks as dynamical systems can capture this channeling of signals into coherent and clearly identifiable stable cellular behaviors, or cellular states. Indeed, qualitative and semi-quantitative dynamical models provide valuable information about the global properties of regulatory networks. The stable steady states of a dynamical system can be interpreted as the set of all possible stable patterns of expression that can be attained within the particular biological network that is being modeled. The advantages of focusing the modeling on the stable steady states of the network are two-fold. First, it reduces the quantity of experimental data required to construct a model, e.g. kinetic and rate limiting step constants, because there is no need to describe the transitory response of the network under different conditions, only the final states. Second, it is easier to verify the predictions of the model experimentally, since it requires stable cellular states to be identified; that is, long-term patterns of activation and not short-lived transitory states of activation that may be difficult to determine experimentally.
Results and discussion
By applying Equation 1 to the network in Figure 2, we obtained Equation 2, which constitutes the discrete version of the dynamical system representing the Th network. Similarly, the continuous version of the Th network was obtained by applying Equation 3 to the network in Figure 2. In this case, however, some of the resulting equations are too large to be presented inside the main text, so we included them as the Additional file 1. Moreover, instead of just typing the equations, we decided to present them in a format that might be used directly to run simulations. The continuous dynamical system of the Th network is included as a plain text file that is able to run on the numerical computation software package GNU Octave http://www.octave.org.
Stable steady states of the dynamical systems. ^{a}
DISCRETE SYSTEM | CONTINUOUS SYSTEM | |||||
---|---|---|---|---|---|---|
Th0 | Th1 | Th2 | Th0 | Th1 | Th2 | |
GATA3 | 0 | 0 | 1 | 0 | 0 | 1 |
IFN-β | 0 | 0 | 0 | 0 | 0 | 0 |
IFN-βR | 0 | 0 | 0 | 0 | 0 | 0 |
IFN-γ | 0 | 1 | 0 | 0 | 0.71443 | 0 |
IFN-γR | 0 | 1 | 0 | 0 | 0.9719 | 0 |
IL-10 | 0 | 0 | 1 | 0 | 0 | 1 |
IL-10R | 0 | 0 | 1 | 0 | 0 | 1 |
IL-12 | 0 | 0 | 0 | 0 | 0 | 0 |
IL-12R | 0 | 0 | 0 | 0 | 0 | 0 |
IL-18 | 0 | 0 | 0 | 0 | 0 | 0 |
IL-18R | 0 | 0 | 0 | 0 | 0 | 0 |
IL-4 | 0 | 0 | 1 | 0 | 0 | 1 |
IL-4R | 0 | 0 | 1 | 0 | 0 | 1 |
IRAK | 0 | 0 | 0 | 0 | 0 | 0 |
JAK1 | 0 | 0 | 0 | 0 | 0.00489 | 0 |
NFAT | 0 | 0 | 0 | 0 | 0 | 0 |
SOCS1 | 0 | 1 | 0 | 0 | 0.89479 | 0 |
STAT1 | 0 | 0 | 0 | 0 | 0.00051 | 0 |
STAT3 | 0 | 0 | 1 | 0 | 0 | 1 |
STAT4 | 0 | 0 | 0 | 0 | 0 | 0 |
STAT6 | 0 | 0 | 1 | 0 | 0 | 1 |
T-bet | 0 | 1 | 0 | 0 | 0.89479 | 0 |
TCR | 0 | 0 | 0 | 0 | 0 | 0 |
Equation 3 defines a highly non-linear continuous dynamical system. In contrast with the discrete system, these continuous equations have to be studied numerically. Numerical methods for solving differential equations require the specification of an initial state, since they proceed via iterations. In our method, we propose to use the stable steady states of the discrete system as the initial states to solve the continuous system that results from application of equation 3 to a given network. We used a standard numerical simulation method to solve the continuous version of the Th model (see Methods). Starting alternatively from each of the three stable steady states found in the discrete model, i.e. the Th0, Th1 and Th2 states, the continuous system was solved numerically until it converged. The continuous system converged to values that could be compared directly with the stable steady states of the discrete system (Table 1). Note that the Th0 and Th2 stable steady states fall in exactly the same position for both the discrete and the continuous dynamical systems, and in close proximity for the Th1 state. This finding highlights the similarity in qualitative behavior of the two models constructed using equations 1 and 3, despite their different mathematical frameworks.
Regions of the state space reached by the continuous version of the Th model, as revealed by a large number of simulations starting from a random initial state. ^{a}
Th0 | Th1 | Th2 | ||||
---|---|---|---|---|---|---|
Avrg. | Std. Dev. | Avrg. | Std. Dev. | Avrg. | Std. Dev. | |
GATA3 | 0.00003 | 0.00008 | 0.00000 | 0.00000 | 0.99997 | 0.00007 |
IFN-β | 0.00000 | 0.00000 | 0.00000 | 0.00000 | 0.00000 | 0.00000 |
IFN-βR | 0.00000 | 0.00001 | 0.00000 | 0.00001 | 0.00000 | 0.00001 |
IFN-γ | 0.00005 | 0.00013 | 0.71438 | 0.00059 | 0.00000 | 0.00001 |
IFN-γR | 0.00004 | 0.00011 | 0.97169 | 0.00040 | 0.00001 | 0.00004 |
IL-10 | 0.00003 | 0.00007 | 0.00000 | 0.00001 | 0.99999 | 0.00004 |
IL-10R | 0.00005 | 0.00010 | 0.00000 | 0.00001 | 0.99999 | 0.00002 |
IL-12 | 0.00000 | 0.00001 | 0.00000 | 0.00000 | 0.00000 | 0.00001 |
IL-12R | 0.00000 | 0.00002 | 0.00000 | 0.00001 | 0.00000 | 0.00001 |
IL-18 | 0.00000 | 0.00001 | 0.00000 | 0.00000 | 0.00000 | 0.00001 |
IL-18R | 0.00000 | 0.00002 | 0.00000 | 0.00001 | 0.00000 | 0.00001 |
IL-4 | 0.00002 | 0.00006 | 0.00000 | 0.00001 | 0.99995 | 0.00011 |
IL-4R | 0.00002 | 0.00004 | 0.00000 | 0.00001 | 0.99990 | 0.00022 |
IRAK | 0.00001 | 0.00005 | 0.00000 | 0.00003 | 0.00001 | 0.00004 |
JAK1 | 0.00002 | 0.00008 | 0.00487 | 0.00005 | 0.00001 | 0.00005 |
NFAT | 0.00001 | 0.00003 | 0.00000 | 0.00002 | 0.00001 | 0.00003 |
SOCS1 | 0.00009 | 0.00022 | 0.89486 | 0.00037 | 0.00002 | 0.00006 |
STAT1 | 0.00001 | 0.00005 | 0.00051 | 0.00003 | 0.00002 | 0.00005 |
STAT3 | 0.00012 | 0.00023 | 0.00001 | 0.00002 | 1.00000 | 0.00002 |
STAT4 | 0.00001 | 0.00003 | 0.00000 | 0.00003 | 0.00000 | 0.00001 |
STAT6 | 0.00001 | 0.00004 | 0.00000 | 0.00002 | 0.99990 | 0.00023 |
T-bet | 0.00007 | 0.00018 | 0.89485 | 0.00036 | 0.00000 | 0.00000 |
TCR | 0.00000 | 0.00001 | 0.00000 | 0.00000 | 0.00000 | 0.00001 |
The whole process resulted in the creation of a model with qualitative characteristics fully comparable to those found in the experimental Th system. Notably, the model used default values for all parameters. Indeed, the continuous dynamical system of the Th network has a total of 58 parameters, all of which were set to the default value of 1, and one parameter (the gain of the sigmoids) with a default value of 10. This set of default values sufficed to capture the correct qualitative behavior of the biological system, namely, the existence of three stable steady states that represent Th0, Th1 and Th2 cells. Readers can run simulations on the model by using the equations provided in the "Th_continuous_model.octave.txt" file. The file was written to allow easy modification of the initial states for the simulations, as well as the values of all parameters.
Analysis of previously published regulatory networks related to Th cell differentiation
Stable steady states of the signaling network in Figure 4
Discrete state 1 | Discrete state 2 | Continuous state 1 | Continuous state 2 | |
---|---|---|---|---|
IFN-γ | 0 | 0 | 0 | 0 |
IL-10 | 0 | 1 | 0 | 0.78995 |
IL-12 | 0 | 0 | 0 | 0 |
IL-4 | 0 | 1 | 0 | 0.89469 |
IL-5 | 0 | 0 | 0 | 0.01343 |
Inf. Resp. | 0 | 0 | 0 | 0.00737 |
Steroids | 0 | 0 | 0 | 0.00105 |
Stable steady states of the signaling network in Figure 5
Discrete state 1 | Discrete state 2 | Discrete state 3 | Discrete state 4 | Discrete state 5 | Discrete state 6 | Discrete state 7 | |
---|---|---|---|---|---|---|---|
CSIF | 0 | 0 | 1 | 0 | 0.5 | 0.5 | 0 |
IFN-γ | 0 | 1 | 0 | 0.5 | 0 | 0 | 0. 5 |
IL-2 | 0 | 1 | 0 | 0.5 | 0.5 | 0.5 | 0 |
IL-4 | 0 | 0 | 1 | 0.5 | 0 | 0.5 | 0.5 |
Continuous state 1 | Continuous state 2 | Continuous state 3 | Continuous state 4 | Continuous state 5 | Continuous state 6 | Continuous state 7 | |
CSIF | 0 | 0.0034416 | 0.8888881 | 0.0034999 | 4.9132E-5 | 0.8881746 | 4.3001E-5 |
IFN-γ | 0 | 0.8888881 | 0.0034416 | 0.8881746 | 4.300E-5 | 0.0034999 | 4.9132E-5 |
IL-2 | 0 | 0.8888881 | 0.0034416 | 0.8881746 | 4.3154E-5 | 0.0035227 | 4.8979E-5 |
IL-4 | 0 | 0.0034416 | 0.8888881 | 0.0035227 | 4.8979E-5 | 0.8881746 | 4.3154E-5 |
Stable steady states of the signaling network in Figure 6
Discrete state 1 | Discrete state 2 | Discrete state 3 | Discrete state 4 | Continuous state 1 | Continuous state 2 | Continuous state 3 | Continuous state 4 | |
---|---|---|---|---|---|---|---|---|
GATA3 | 0 | 0 | 11 | 1 | 0 | 0 | 0.93037 | 0.93037 |
IFN-γ | 0 | 1 | 0 | 1 | 0 | 0.99914 | 0 | 0.90967 |
IFN-γR | 0 | 1 | 0 | 1 | 0 | 0.99997 | 0 | 0.99617 |
IL-12 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
IL-12R | 0 | 1 | 0 | 0 | 0 | 0.9096 | 0 | 0.00193 |
IL-13 | 0 | 0 | 1 | 1 | 0 | 0 | 0.99719 | 0.99719 |
IL-4 | 0 | 0 | 1 | 1 | 0 | 0 | 0.99719 | 0.99719 |
IL-4R | 0 | 0 | 1 | 1 | 0 | 0 | 0.99991 | 0.99991 |
IL-5 | 0 | 0 | 1 | 1 | 0 | 0 | 0.99719 | 0.99719 |
STAT1 | 0 | 1 | 0 | 1 | 0 | 1 | 0 | 0.99988 |
STAT4 | 0 | 1 | 0 | 0 | 0 | 0.99617 | 0 | 2.4E-4 |
STAT6 | 0 | 0 | 1 | 1 | 0 | 0 | 1 | 1 |
T-bet | 0 | 1 | 0 | 1 | 0 | 0.93037 | 0 | 0.93034 |
TCR | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
Stable steady states of the signaling network in Figure 7
Discrete state 1 | Discrete state 2 | Continuous state 1 | Continuous state 2 | |
---|---|---|---|---|
Ag/MHC | 0 | 0 | 0 | 0 |
ATF2 | 0 | 0 | 0 | 0 |
c-Maf | 0 | 0 | 0 | 0 |
CD4 | 0 | 0 | 0 | 0 |
GATA3 | 0 | 1 | 0 | 0.99999 |
IFN-γ | 0 | 0 | 0 | 0 |
IL-12 | 0 | 0 | 0 | 0 |
IL-12R | 0 | 0 | 0 | 0 |
IL-13 | 0 | 1 | 0 | 0.8468 |
IL-18 | 0 | 0 | 0 | 0 |
IL-18R | 0 | 0 | 0 | 0 |
IL-4 | 0 | 1 | 0 | 0.8468 |
IL-4R | 0 | 1 | 0 | 0.99176 |
IL-5 | 0 | 1 | 0 | 0.8469 |
IRAK | 0 | 0 | 0 | 0 |
Itk | 0 | 0 | 0 | 0 |
JNK | 0 | 0 | 0 | 0 |
JNK2 | 0 | 0 | 0 | 0 |
Lck | 0 | 0 | 0 | 0 |
MKK3 | 0 | 0 | 0 | 0 |
NFAT | 0 | 0 | 0 | 0 |
NFkB | 0 | 0 | 0 | 0 |
p38/MAPK | 0 | 0 | 0 | 0 |
STAT4 | 0 | 0 | 0 | 0 |
STAT6 | 0 | 1 | 0 | 0.99975 |
T-bet | 0 | 0 | 0 | 0 |
TCR | 0 | 0 | 0 | 0 |
TRAF6 | 0 | 0 | 0 | 0 |
It is important to note that the failure of these four alternative networks to capture the three states representing Th cells is not attributable to the use of very simplistic and/or outdated data. Indeed, the network in Figure 6 comes from a relatively recent review, while that in Figure 7 is rather complex and contains five more nodes than our own proposed network (Figure 2). All this stresses the importance of using a correctly reconstructed network to develop dynamical models, either with our approach or any other.
Conclusion
There is a great deal of interest in the reconstruction and analysis of regulatory networks. Unfortunately, kinetic information about the elements that constitute a network or pathway is not easily gathered, and hence the analysis of its dynamical properties (via simulation packages such as [13]) is severely restricted to a small set of well-characterized systems. Moreover, the translation from a static to a dynamical representation normally requires the use of a network-specific set of equations to represent the expression or concentration of every molecule in the system.
We herein propose a method for generating a system of ordinary differential equations to construct a model of a regulatory network. Since the equations can be unambiguously applied to any signaling or regulatory network, the construction and analysis of the model can be carried out systematically. Moreover, the process of finding the stable steady states is based on the application of an analytical method (generalized logical analysis [14, 15] on a discrete version of the model), followed by a numerical method (on the continuous version) starting from specific initial states (the results obtained from the logical analysis). This characteristic allows a fully automated implementation of our methodology for modeling. In order to construct the equations of the continuous dynamical system with the exclusive use of the topological information from the network, the equations have to incorporate a set of default values for all the parameters. Therefore, the resulting model is not optimized in any sense. However, the advantage of using Equation 3 is that the user can later modify the parameters so as to refine the performance of the model, approximating it to the known behavior of the biological system under study. In this way, the user has a range of possibilities, from a purely qualitative model to one that is highly quantitative.
There are studies that compare the dynamical behavior of discrete and continuous dynamical systems. Hence, it is known that while the steady state of a Boolean model will correspond qualitatively to an analogous steady state in a continuous approach, the reverse is not necessarily true. Moreover, periodic solutions in one representation may be absent in the other [16]. This discrepancy between the discrete and continuous models is more evident for steady states where at least one of the nodes has an activation state precisely at, or near, its threshold of activation. Because of this characteristic, discrete and continuous models for a given regulatory network differ in the total number of steady states [17]. For this reason, our method focuses on the study of only one type of steady state; namely, the regular stationary points [18]. These points do not have variables near an activation threshold, and they are always stable steady states. Moreover, it has been shown that this type of stable steady state can be found in discrete models, and then used to locate their analogous states in continuous models of a given genetic regulatory network [19].
It is beyond the scope of this paper to present a detailed mathematical analysis of the dynamical system described by Equation 3. Instead, we present a framework that can help to speed up the analysis of the qualitative behavior of signaling networks. Under this perspective, the usefulness of our method will ultimately be determined through building and analyzing concrete models. To show the capabilities of our proposed methodology, we applied it to analysis of the regulatory network that controls differentiation in T helper cells. This biological system was well suited to evaluating our methodology because the network contains several known components, and it has three alternative stable patterns of activation. Moreover, it is of great interest to understand the behavior of this network, given the role of T helper cell subsets in immunity and pathology [20]. Our method applied to the Th network generated a model with the same qualitative behavior as the biological system. Specifically, the model has three stable states of activation, which can be interpreted as the states of activation found in Th0, Th1 and Th2 cells. In addition, the system is capable of being moved from the Th0 state to either the Th1 or Th2 states, given a sufficiently large IFN-γ or IL-4 signal, respectively. This characteristic reflects the known qualitative properties of IFN-γ and IL-4 as key cytokines that control the fate of T helper cell differentiation.
Regarding the numerical values returned by the model, it is not possible yet to evaluate their accuracy, given that (to our knowledge) no quantitative experimental data are available for this biological system. The resulting model, then, should be considered as a qualitative representation of the system. However, representing the nodes in the network as normalized continuous variables will eventually permit an easy comparison with quantitative experimental data whenever they become available. Towards this end, the equations in our methodology define a sigmoid function, with values ranging from 0 to 1, regardless of the values of assigned to the parameters in the equations. This characteristic has been used before to represent and model the response of signaling pathways [21, 22]. It is important to note, however, that the modification of the parameters allow the model to be fitted against experimental data.
One benefit of a mathematical model of a particular biological network is the possibility of predicting the behavior of complex experimental setups. Therefore, it is important to be aware of its limitations beforehand, to avoid generating experimental data that cannot be handled by the model. The method we present in this paper has been developed to obtain the number and relative position of the stable steady states of a regulatory network. Equations 1 and 3 include a number of parameters that allow the response of the model to be fine-tuned, but the equations were not designed to describe the transitory responses of molecules with great detail. Therefore, failure to predict a stable steady state with high numerical accuracy should not be interpreted as a failure of the approach presented here. By contrast, failure to describe and/or predict the number and approximate location of stable steady states under a wide range of values for the parameters would call the validity of the reconstruction of a particular network into question. Here, however, it is essential to establish the validity of the network used as input. Indeed, we applied our method to four alternative forms of the network that regulates Th cell differentiation. The alternative networks (Figures 4 through 7) were taken from previously published attempts to discover the molecular basis of this differentiation process. Originally, such networks were not developed with the idea of studying dynamical properties. It is not surprising, then, that these networks do not reflect the existence of three stable steady states, representing the molecular states of Th0, Th1 and Th2 cells, respectively. In these cases, the failure to find the correct stable steady states is not a problem in the modeling methodology, but a problem in the inference of the regulatory network.
In conclusion, we have shown that the creation of a dynamical model of a regulatory network can be considerably simplified with the aid of a standardized set of equations, where the feature that distinguishes one molecule from another is the number of regulatory inputs. Such standardization permits a continuous dynamical system to be systematically and analytically constructed together with a basic analysis of its global properties, based exclusively on the information provided by the connectivity of the network. While the use of a standardized set of functions to model a network may severely restrict the capability to fit specific datasets, we believe that the loss in flexibility is balanced by the possibility of rapidly developing models and gaining knowledge of the dynamical behavior of a network, especially in those cases where few kinetic data are available. Thus, we provide a method for incorporating the dynamical perspective in the analysis of regulatory networks, using the topological information of a network, without the need to collect extensive time-series or kinetic data.
Methods
Molecular basis of the Th network topology
The following paragraphs detail the evidence used to infer the topology of the Th regulatory network, updating the data summarized in [11]. Th1 cells are producers of IFN-γ [10, 23], which acts on its target cells by binding to a cell-membrane receptor [24–26] to start a signaling cascade, which involves JAK1 and STAT-1 [27–29]. STAT-1 can be activated by a number of ligands besides IFN-γ, but importantly, it cannot be activated by IL-4 [30], which is a major Th2 signal. In contrast, STAT-1 plays a role in modulating IL-4, being an intermediate in the negative regulation of IFN-γ exerted on IL-4 expression [31]. Different signals converge in STAT-1, among them that of IFN-β/IFN-βR [32]. The IFN-γ signaling continues downstream to activate SOCS-1 in a STAT-1-dependent pathway [33, 34]. SOCS-1, in turn, influences both the IFN-γ and IL-4 pathways. On the one hand, SOCS-1 is a negative regulator of IFN-γ signaling, blocking the interaction of IFN-γR and STAT-1 [35] due to direct inhibition of JAK1 [29, 36]. On the other hand, SOCS-1 blocks the IL-4R/STAT-6 pathway [37]. SOCS-1 is, therefore, a key element for the inhibition from the IFN-γ to the IL-4 pathway. Th1 cells express high levels of SOCS-1 mRNA, while it is barely detectable in Th0 and Th2 cells [38]. Finally, another key molecule is T-bet, which is a transcription factor detected in Th1 but not Th0 or Th2 cells. T-bet expression is upregulated by IFN-γ in a STAT-1-dependent mechanism [39]. Importantly, T-bet is an inhibitor of GATA-3 [40], an activator of IFN-γ [40] and activator of T-bet itself [41, 42].
Th2 cells express IL-4, which is the major known determinant of the Th2 phenotype itself [43]. IL-4 binds to its receptor, IL-4R, which is preferentially expressed in Th2 cells [23, 44]. The IL-4R signaling is transduced by STAT-6, which in turn activates GATA-3 [10]. GATA-3, in turn, is capable of inducing IL-4 [45], thus establishing a feedback loop. The influence from the IL-4 pathway on the IFN-γ pathway seems to be mediated by GATA-3 via STAT-4 [46]. Like T-bet, GATA-3 also presents a self-activation loop [47–49].
IL-12 and IL-18 are two molecules that affect the IFN-γ pathway. IL-12 is a cytokine produced by monocytes and dendritic cells and promotes the development of Th1 cells [50]. The IL-12 receptor is present in its functional form in Th0 and Th1 but not Th2 cells [51]. IL-12R signaling is mediated by STAT-4 [52], which is able to activate IFN-γ [41, 46, 53]. The IL-12 signaling pathway can be blocked by IL-4 by the STAT-6 dependent down-regulation of one subunit of IL-12R [54]. IL-18 is a cytokine produced by many cell types and promotes IFN-γ production in Th cells [55]. It acts upon binding to its receptor, IL-18R, which acts through IRAK [56]. IL-12 and IL-18 act synergistically to increase IFN-γ production, but using different pathways [57, 58]. Finally, IL-4 is able to block IL-18 signaling in a STAT-6 dependent manner [59].
IL-10 is a cytokine actively produced by Th2 cells, and it inhibits cytokine production by Th1 cells. As with the other cytokines mentioned above, IL-10 acts upon binding to a cell surface receptor, IL-10R, which in turn activates the STAT signaling system [60]. In particular, it has been shown that the functioning of IL-10 signaling is dependent upon the presence of STAT-3 [61]. As for the signals affecting IL-10 expression, it has been shown that IL-4 enhances IL-10 gene expression in Th2 but not Th1 cells [62]. This requirement implies that the intracellular signaling from IL-4 to IL-10 should pass through a Th2 specific molecule, which from the molecules considered here can only be GATA-3. Finally, IL-10 has been shown to be a very powerful inhibitor of IFN-γ production [60, 63].
Cytokine gene expression in T cells is induced by the activation of the T cell receptor (TCR) by ligand binding. Different signaling pathways are activated by the TCR [64]. Among these is the pathway including the NFAT family of transcription factors, which are implicated in the T cell activation-dependent regulation of numerous cytokines. A constitutively active form of one of the NFAT proteins, specifically NFATc1, increases the expression of IFN-γ [65]. Importantly, the same experimental procedure does not affect the expression of IL-4. All this indicates that the NFAT family members play a central role in the TCR-induced expression of cytokines during Th cell differentiation, especially in the Th1 pathway.
The discrete dynamical system
The discrete system represents the network as a series of interconnected elements that have only two possible states of activation, 0 (or inactive) and 1 (or active). Given this property, the network is completely described by the following set of Boolean equations:
Equation 1.
A node x in the network can have only one of three possible forms depending on whether it has activator and inhibitor input nodes, or only activators, or only inhibitors. In the first case, i.e. form § in Eqn.1, the Boolean function can be read as: x will be active in the next time step if at this time any of its activators and none of its inhibitors are acting upon it. Similarly, form §§ can be translated as: x will be active if any of its activators is acting upon it. And finally, form §§§ reads as: x will be active if none of its inhibitors are acting upon it. Note than in all cases inhibitors are strong enough to change the state of a node from 1 to 0, while activators are strong enough to change the state of a node from 0 to 1 if no inhibitor is acting on the node of reference. The three alternative forms of representing a node in Equation 1 imply two possible default states of activation, i.e. the state of a node when there are neither activators nor inhibitors acting upon it. If the connectivity of the node includes either only positive inputs, or both positive and negative inputs, then the node has an inactive state by default. Alternatively, if the connectivity of a node has only negative inputs, then the node has an active state by default.
The Th network (Figure 2) can be converted into a discrete dynamical system using Equation 1. The resulting system of equations is as follows:
Equation 2.
GATA 3(t + 1) = (GATA 3(t) ∨ STAT 6(t)) ∧ ¬(T - bet(t))
IFN - βR(t + 1) = IFN - β(t)
IFN - γ(t + 1) = (IRAK(t) ∨ NFAT(t) ∨ STAT - 4(t) ∨ T - bet(t)) ∧ ¬(STAT 3(t))
IFN - γR(t + 1) = IFN - γ(t)
IL - 10(t + 1) = GATA 3(t)
IL - 10R(t + 1) = IL - 10(t)
IL - 12R(t + 1) = IL - 12(t)
IL - 18R(t + 1) = IL - 18(t) ∧ ¬(STAT6(t))
IL - 4(t + 1) = GATA 3(t) ∧ ¬(STAT 1(t))
IL - 4R(t + 1) = IL - 4(t) ∧ ¬(SOCS 1(t))
IRAK(t + 1) = IL - 18R(t)
JAK 1(t + 1) = IFN - γR(t) ∧ ¬(SOCS 1(t))
NFAT(t + 1) = TCR(t)
SOCS 1(t + 1) = STAT 1(t) ∨ T - bet(t)
STAT 1(t + 1) = IFN - βR(t) ∨ JAK 1(t)
STAT 3(t + 1) = IL - 10R(t)
STAT 4(t + 1) = IL - 12R(t) ∧ ¬(GATA 3(t))
STAT 6(t + 1) = IL - 4R(t)
T - bet(t + 1) = (STAT 1(t) ∨ T - bet(t)) ∧ ¬(GATA 3(t))
Notice that there are only 19 equations out of a total of 23 elements in the Th network. The reason is that four elements, namely IFN-β, IL-12, IL-18 and TCR, do not have inputs. These four elements are thus treated as constants, since there are no interactions that regulate their behavior. Throughout the text, these four elements are considered as having a value of 0.
Stable steady states of the discrete system
Circuits of the Th network ^{a}
1 | IFNγ→IFNγR→JAK1→STAT1¬IL4→IL4R→STAT6¬IL18R→IRAK→ |
---|---|
2 | IFNγ→IFNγR→JAK1→STAT1¬IL4→IL4R→STAT6¬IL12R→STAT4→ |
3 | IFNγ→IFNγR→JAK1→STAT1¬IL4→IL4R→STAT6→GATA3→IL10→IL10R→STAT3¬ |
4 | IFNγ→IFNγR→JAK1→STAT1¬IL4→IL4R→STAT6→GATA3¬STAT4→ |
5 | IFNγ→IFNγR→JAK1→STAT1¬IL4→IL4R→STAT6→GATA3¬Tbet→ |
6 | IFNγ→IFNγR→JAK1→STAT1→SOCS1¬IL4R→STAT6¬IL18R→IRAK→ |
7 | IFNγ→IFNγR→JAK1→STAT1→SOCS1¬IL4R→STAT6¬IL12R→STAT4→ |
8 | IFNγ→IFNγR→JAK1→STAT1→SOCS1¬IL4R→STAT6→GATA3→IL10→IL10R→STAT3¬ |
9 | IFNγ→IFNγR→JAK1→STAT1→SOCS1¬IL4R→STAT6→GATA3¬STAT4→ |
10 | IFNγ→IFNγR→JAK1→STAT1→SOCS1¬IL4R→STAT6→GATA3¬Tbet→ |
11 | IFNγ→IFNγR→JAK1→STAT1→Tbet→ |
12 | IFNγ→IFNγR→JAK1→STAT1→Tbet→SOCS1¬IL4R→STAT6¬IL18R→IRAK→ |
13 | IFNγ→IFNγR→JAK1→STAT1→Tbet→SOCS1¬IL4R→STAT6¬IL12R→STAT4→ |
14 | IFNγ→IFNγR→JAK1→STAT1→Tbet→SOCS1¬IL4R→STAT6→GATA3→IL10→IL10R→STAT3¬ |
15 | IFNγ→IFNγR→JAK1→STAT1→Tbet→SOCS1¬IL4R→STAT6→GATA3¬STAT4→ |
16 | IFNγ→IFNγR→JAK1→STAT1→Tbet¬GATA3→IL4→IL4R→STAT6¬IL18R→IRAK→ |
17 | IFNγ→IFNγR→JAK1→STAT1→Tbet¬GATA3→IL4→IL4R→STAT6¬IL12R→STAT4→ |
18 | IFNγ→IFNγR→JAK1→STAT1→Tbet¬GATA3→IL10→IL10R→STAT3¬ |
19 | IFNγ→IFNγR→JAK1→STAT1→Tbet¬GATA3¬STAT4→ |
20 | IL4→IL4R→STAT6→GATA3→ |
21 | IL4R→STAT6→GATA3¬ Tbet→SOCS1¬ |
22 | Tbet→ |
23 | Tbet¬GATA3¬ |
24 | GATA3→ |
25 | IL4→IL4R→STAT6→GATA3¬Tbet→SOCS1¬JAK1→STAT1¬ |
26 | JAK1→STAT1→SOCS1¬ |
27 | JAK1→STAT1→Tbet→ SOCS1¬ |
The continuous dynamical system
To describe the network as a continuous dynamical system, we use the following set of ordinary differential equations:
Equation 3.
Once a network is translated to a dynamical system using Equation 3, it is necessary to specify values for all parameters. For a system with n nodes and m interactions, there are m+2n parameters. However, there are usually insufficient experimental data to assign realistic values for each and every one of the parameters. Nevertheless, it is possible to use a series of default values for all the parameters in Equation 3. The reason is that, as we showed in the previous paragraph, the equations have the same qualitative shape for any value assigned to the parameters. Hence, for the sake of simplicity, it is possible to assign the same values to most of the parameters, as a first approach. For the present study on the Th model, we use a value of 1 for all α s, β s and γ s; and we use h = 10, since we currently lack quantitative data to estimate more realistic values. Moreover, the use of default values ensures the possibility of creating the dynamical system in a fully automated way. Nonetheless, after the initial construction and analysis of the resulting system, the modeler may modify the values of the parameters so as to fine-tune the dynamical behavior of the equations, whenever more experimental quantitative data become available. The continuous dynamical system of the Th model, constructed with the use of Equation 3, yields a system of 23 equations, which is included in the file "Th_continuous_model.octave.txt".
Stable steady states of the continuous system
Nonlinear systems of ordinary differential equations are studied numerically. Hence the continuous dynamical system defined by Equation 3 poses the problem of how to find all its stable steady states without using very time-consuming and computing-intensive methods. This is where the creation of two dynamical systems of the same network, one discrete and one continuous, bears fruit. Since a Boolean (step) function is a limiting case of a very steep sigmoid curve, networks made of binary elements share many qualitative features with systems modeled using continuous functions [68]. Indeed, it has been shown [19] that the qualitative information resulted from generalized logical analysis can be directly used to find the number, nature and approximate location of the steady states of a system of differential equations representing the same network. We therefore decided to use this characteristic to speed up the process of finding all the stable steady states in the continuous dynamical system. Specifically, the stable steady states of the discrete system are used as initial states to solve the differential equations, running them until the system converges to its own stable steady states. Calculating the convergence of a system of ordinary differential equations from a given initial state is a straightforward procedure using any numerical solver. For our simulations we used the lsode function of the GNU Octave package http://www.octave.org, stopping the numerical integration when all the variables of the system changed by less than 10^{-4} for at least 10 consecutive steps of the procedure. The final values of the variables in the system are considered to be the stable steady states of the continuous model of the network.
Implementation
The methodology was fully implemented in a java program, and it has been tested under a linux environment using java version 1.5.0 (JRE 5.0), as well as octave version 2.1.34. The bytecode version of the program is included as Additional file 2.
Declarations
Acknowledgements
We want to thank Massimo de Francesco, Mark Ibberson, Caroline Johnson-Leger, Maria Karmirantzou, Lukasz Salwinski, François Talabot and Francisca Zanoguera for their valuable comments and suggestions.
Authors’ Affiliations
References
- Jeong H, Tombor B, Albert R, Oltvai ZN, Barabási AL: The large-scale organization of metabolic networks. Nature. 2000, 407: 651-654. 10.1038/35036627.View ArticlePubMedGoogle Scholar
- Milo R, Shen-Orr S, Itzkovitz S, Kashtan N, Chklovskii D, Alon U: Network motifs: Simple building blocks of complex networks. Science. 2002, 298: 824-827. 10.1126/science.298.5594.824.View ArticlePubMedGoogle Scholar
- Ravasz E, Somera AL, Mongru DA, Oltvai ZN, Barabási AL: Hierarchical organization of modularity in metabolic networks. Science. 2002, 297: 1551-1555. 10.1126/science.1073374.View ArticlePubMedGoogle Scholar
- Covert MW, Schilling CH, Famili I, Edwards JS, Goryanin II, Selkov E, Palsson BO: Metabolic modeling of microbial strains in silico. Trends Biochem Sci. 2001, 26: 179-186. 10.1016/S0968-0004(00)01754-0.View ArticlePubMedGoogle Scholar
- Herrgård MJ, Covert MW, Palsson BØ: Reconstruction of microbial transcriptional regulatory networks. Curr Opin Biotechnol. 2004, 15: 70-77. 10.1016/j.copbio.2003.11.002.View ArticlePubMedGoogle Scholar
- Mendoza L, Thieffry D, Alvarez-Buylla ER: Genetic control of flower morphogenesis in Arabidopsis thaliana: a logical analysis. Bioinformatics. 1999, 15: 593-606. 10.1093/bioinformatics/15.7.593.View ArticlePubMedGoogle Scholar
- Sánchez L, Thieffry D: Segmenting the fly embryo: a logical analysis of the pair-rule cross-regulatory module. J Theor Biol. 2003, 224: 517-537. 10.1016/S0022-5193(03)00201-7.View ArticlePubMedGoogle Scholar
- de Jong H: Modeling and simulation of genetic regulatory systems: a literature review. J Comp Biol. 2002, 9: 67-103. 10.1089/10665270252833208.View ArticleGoogle Scholar
- Lok L, Brent R: Automatic generation of cellular reaction networks with Moleculizer 1.0. Nature Biotechnol. 2005, 23: 131-136. 10.1038/nbt1054.View ArticleGoogle Scholar
- Murphy KM, Reiner SL: The lineage decisions on helper T cells. Nat Rev Immunol. 2002, 2: 933-944. 10.1038/nri954.View ArticlePubMedGoogle Scholar
- Mendoza L: A network model for the control of the differentiation process in Th cells. BioSystems.Google Scholar
- Remy E, Ruet P, Mendoza L, Thieffry D, Chaouiya C: From Logical Regulatory Graphs to Standard Petri Nets: Dynamic Roles and Functionality of Feedback Circuits. Transactions on Computational Systems Biology.Google Scholar
- Mendes P: Biochemistry by numbers: simulation of biochemical pathways with Gepasi 3. Trends Biochem Sci. 1997, 22: 361-363. 10.1016/S0968-0004(97)01103-1.View ArticlePubMedGoogle Scholar
- Thomas R: Regulatory networks seen as asynchronous automata: a logical description. J Theor Biol. 1991, 153: 1-23.View ArticleGoogle Scholar
- Thomas R, Thieffry D, Kaufman M: Dynamical behaviour of biological regulatory networks-I. Biological role of feedback loops and practical use of the concept of the loop-characteristic state. Bull Math Biol. 1995, 57: 247-276. 10.1016/0092-8240(94)00036-C.View ArticlePubMedGoogle Scholar
- Glass L, Kaufman M: The logical analysis of continuous, non-linear biochemical control networks. J Theor Biol. 1973, 39: 103-129. 10.1016/0022-5193(73)90208-7.View ArticlePubMedGoogle Scholar
- Mochizuki A: An analytical study of the number of steady states in gene regulatory networks. J Theor Biol. 2005, 236: 291-310. 10.1016/j.jtbi.2005.03.015.View ArticlePubMedGoogle Scholar
- Snoussi EH, Thomas R: Logical identification of all steady states: the concept of feedback loop characteristic states. Bull Math Biol. 1993, 55: 973-991. 10.1016/S0092-8240(05)80199-5.View ArticleGoogle Scholar
- Muraille E, Thieffry D, Leo O, Kaufman M: Toxicity and neuroendocrine regulation of the immune response: a model analysis. J Theor Biol. 1996, 183: 285-305. 10.1006/jtbi.1996.0221.View ArticlePubMedGoogle Scholar
- Singh VK, Mehrotra S, Agarwal SS: The paradigm of Th1 and Th2 cytokines: its relevance to autoimmunity and allergy. Immunol Res. 1999, 20: 147-161.View ArticlePubMedGoogle Scholar
- Hautaniemi S, Kharait S, Iwabu A, Wells A, Lauffenburger DA: Modeling of signal-response cascades using decision tree analysis. Bioinformatics. 2005, 21: 2027-2035. 10.1093/bioinformatics/bti278.View ArticlePubMedGoogle Scholar
- Sachs K, Perez O, Pe'er D, Lauffenburger DA, Nolan GP: Causal protein-signaling networks derived from multiparameter single-cell data. Science. 2005, 308: 523-529. 10.1126/science.1105809.View ArticlePubMedGoogle Scholar
- Hamalainen H, Zhou H, Chou W, Hashizume H, Heller R, Lahesmaa R: Distinct gene expression profiles of human type 1 and type 2 T helper cells. Genome Biology. 2001, 2: 1-0022. 10.1186/gb-2001-2-7-research0022.View ArticleGoogle Scholar
- Groux H, Sornasse T, Cottrez F, de Vries JE, Coffman RL, Roncarolo MG, Yssel H: Induction of human T helper cell type 1 differentiation results in loss of IFN-γ receptor β-chain expression. J Immunol. 1997, 158: 5627-5631.PubMedGoogle Scholar
- Novelli F, D'Elios MM, Bernabei P, Ozmen L, Rigamonti L, Almerigogna F, Forni G, Del Prete G: Expression and role in apoptosis of the α- and β-chains of the IFN-γ receptor in human Th1 and Th2 clones. J Immunol. 1997, 159: 206-213.PubMedGoogle Scholar
- Rigamonti L, Ariotti S, Losana G, Gradini R, Russo MA, Jouanguy E, Casanova JL, Forni G, Novelli F: Surface expression of the IFN-γR2 chain is regulated by intracellular trafficking in human T lymphocytes. J Immunol. 2000, 164: 201-207.View ArticlePubMedGoogle Scholar
- Kotenko SV, Pestka S: Jak-Stat signal transduction pathway through the eyes of cytokine class II receptor complexes. Oncogene. 2000, 19: 2557-2565. 10.1038/sj.onc.1203524.View ArticlePubMedGoogle Scholar
- Kerr IM, Costa-Pereira AP, Lillemeier BF, Strobl B: Of JAKs, STATs, blind watchmakers, jeeps and trains. FEBS Lett. 2003, 546: 1-5. 10.1016/S0014-5793(03)00411-3.View ArticlePubMedGoogle Scholar
- Krebs DL, Hilton DJ: SOCS proteins: negative regulators of cytokine signaling. Stem Cells. 2001, 19: 378-387. 10.1634/stemcells.19-5-378.View ArticlePubMedGoogle Scholar
- Moriggl R, Kristofic C, Kinzel B, Volarevic S, Groner B, Brinkmann V: Activation of STAT proteins and cytokine genes in human Th1 and Th2 cells generated in the absence of IL-12 and IL-4. J Immunol. 1998, 160: 3385-3392.PubMedGoogle Scholar
- Elser B, Lohoff M, Kock S, Giaisi M, Kirchhoff S, Krammer PH, Li-Weber M: IFN-γ represses IL-4 expression via IRF-1 and IRF-2. Immunity. 2002, 17: 703-712. 10.1016/S1074-7613(02)00471-5.View ArticlePubMedGoogle Scholar
- Goodbourn S, Didcock L, Randal RE: Interferons: cell signalling, immune modulation, antiviral responses and virus countermeasures. J Gen Virol. 2000, 81: 2341-2364.View ArticlePubMedGoogle Scholar
- Chen XP, Losman JA, Rothman P: SOCS proteins, regulators of intracellular signaling. Immunity. 2000, 13: 287-290. 10.1016/S1074-7613(00)00028-5.View ArticlePubMedGoogle Scholar
- Saito H, Morita Y, Fujimoto M, Narazaki M, Naka T, Kishimoto T: IFN regulatory factor-1-mediated transcriptional activation of mouse STAT-induced STAT inhibitor-1 gene promoter by IFN-γ. J Immunol. 2000, 164: 5833-5843.View ArticlePubMedGoogle Scholar
- Diehl S, Anguita J, Hoffmeyer A, Zapton T, Ihle JN, Fikrig E, Rincón M: Inhibition of Th1 differentiation by IL-6 is mediated by SOCS1. Immunity. 2000, 13: 805-815. 10.1016/S1074-7613(00)00078-9.View ArticlePubMedGoogle Scholar
- Zhang JG, Metcalf D, Rakar S, Asimakis M, Greenhalgh CJ, Willson TA, Starr R, Nicholson SE, Carter W, Alexander WS, Hilton J: The SOCS box of suppressor of cytokine signaling-1 is important for inhibition of cytokine action in vivo. Proc Natl Acad Sci USA. 2001, 98: 13261-13265. 10.1073/pnas.231486498.PubMed CentralView ArticlePubMedGoogle Scholar
- Losman JA, Chen XP, Hilton D, Rothman P: Cutting edge: SOCS-1 is a potent inhibitor of IL-4 signal transduction. J Immunol. 1999, 162: 3770-3774.PubMed CentralPubMedGoogle Scholar
- Egwuagu CE, Yu CR, Zhang M, Mahdi RM, Kim SJ, Gery I: Suppressors of cytokine signaling proteins are differentially expresses in Th1 and Th2 cells: implications for the Th cell lineage commitment and maintenance. J Immunol. 2002, 168: 3181-3187.View ArticlePubMedGoogle Scholar
- Lighvani A., Frucht DM, Jankovic D, Yamane H, Aliberti J, Hissong BD, Nguyen BV, Gadina M, Sher A, Paul WE, O'Shea JJ: T-bet is rapidly induced by interferon-γ in lymphoid and myeloid cells. Proc Natl Acad Sci USA. 2001, 98: 15137-15142. 10.1073/pnas.261570598.PubMed CentralView ArticlePubMedGoogle Scholar
- Szabo SJ, Kim ST, Costa GL, Zhang X, Fathman CG, Glimcher LH: A novel transcription factor, T-bet, directs Th1 lineage commitment. Cell. 2000, 100: 655-669. 10.1016/S0092-8674(00)80702-3.View ArticlePubMedGoogle Scholar
- Mullen AC, High FA, Hutchins AS, Lee HW, Villarino AV, Livingston DM, Kung AL, Cereb N, Yao TP, Yang SY, Reiner SL: Role of T-bet in commitment of TH1 cells before IL-12-dependent selection. Science. 2001, 292: 1907-1910. 10.1126/science.1059835.View ArticlePubMedGoogle Scholar
- Zhang Y, Apilado R, Coleman J, Ben-Sasson S, Tsang S, Hu-Li J, Paul WE, Huang H: Interferon γ stabilizes the T helper cell type 1 phenotype. J Exp Med. 2001, 194: 165-172. 10.1084/jem.194.2.165.PubMed CentralView ArticlePubMedGoogle Scholar
- Agnello D, Lankford CSR, Bream J, Morinobu A, Gadina M, O'Shea J, Frucht DM: Cytokines and transcription factors that regulate T helper cell differentiation: new players and new insights. J Clin Immunol. 2003, 23: 147-161. 10.1023/A:1023381027062.View ArticlePubMedGoogle Scholar
- Nelms K, Keegan AD, Zamorano J, Ryan JJ, Paul WE: The IL-4 receptor: signaling mechanisms and biologic functions. Annu Rev Immunol. 1999, 17: 701-738. 10.1146/annurev.immunol.17.1.701.View ArticlePubMedGoogle Scholar
- Ouyang W, Löhning M, Gao Z, Assenmacher M, Ranganath S, Radbruch A, Murphy KM: Stat6-independent GATA-3 autoactivation directs IL-4-independent Th2 development and commitment. Immunity. 2000, 12: 27-37. 10.1016/S1074-7613(00)80156-9.View ArticlePubMedGoogle Scholar
- Usui T, Nishikomori R, Kitani A, Strober W: GATA-3 suppresses Th1 development by downregulation of Stat4 and not through effects on IL-12Rβ2 chain or T-bet. Immunity. 2003, 18: 415-428. 10.1016/S1074-7613(03)00057-8.View ArticlePubMedGoogle Scholar
- Zhou M, Ouyang W, Gong Q, Katz SG, White JM, Orkin SH, Murphy KM: Friend of GATA-1 Represses GATA-3-dependent activity in CD4+ cells. J Exp Med. 2001, 194: 1461-1471. 10.1084/jem.194.10.1461.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhou M, Ouyang W: The function role of GATA-3 in Th1 and Th2 differentiation. Immunol Res. 2003, 28: 25-37. 10.1385/IR:28:1:25.View ArticlePubMedGoogle Scholar
- Höfer T, Nathansen H, Löhning M, Radbruch A, Heinrich R: GATA-3 transcriptional imprinting in Th2 lymphocytes: A mathematical model. Proc Natl Acad Sci USA. 2002, 99: 9364-9368. 10.1073/pnas.142284699.PubMed CentralView ArticlePubMedGoogle Scholar
- Trinchieri G: Interleukin-12: a proinflammatory cytokine with immunoregulatory functions that bridge innate resistance and antigen-specific adaptive immunity. Annu Rev Immunol. 1995, 13: 251-276.View ArticlePubMedGoogle Scholar
- Szabo SJ, Jacobson NG, Dighe AS, Gubler U, Murphy KM: Developmental commitment to the Th2 lineage by extinction of IL-12 signaling. Immunity. 1995, 2: 665-675. 10.1016/1074-7613(95)90011-X.View ArticlePubMedGoogle Scholar
- Thierfelder WE, van Deursen JM, Yamamoto K, Tripp RA, Sarawar SR, Carson RT, Sangster MY, Vignali DA, Doherty PC, Grosveld GC, Ihle JN: Requirement for Stat4 in interleukin-12-mediated responses of natural killer and T cells. Nature. 1996, 382: 171-174. 10.1038/382171a0.View ArticlePubMedGoogle Scholar
- Kaplan MH, Sun YL, Hoey T, Grusby MJ: Impaired IL-12 responses and enhanced development of Th2 cells in Stat4-deficient mice. Nature. 1996, 382: 174-177. 10.1038/382174a0.View ArticlePubMedGoogle Scholar
- Szabo SJ, Dighe AS, Gubler U, Murphy KM: Regulation of the interleukin (IL)-12R β2 subunit expression in developing T helper 1 (Th1) and Th2 cells. J Exp Med. 1997, 185: 817-824. 10.1084/jem.185.5.817.PubMed CentralView ArticlePubMedGoogle Scholar
- Swain SL: Interleukin 18: tipping the balance towards a T helper cell 1 response. J Exp Med. 2001, 194: F11-F14. 10.1084/jem.194.3.F11.PubMed CentralView ArticlePubMedGoogle Scholar
- Chang JT, Segal BM, Nakanishi K, Okamura H, Shevach EM: The costimulatory effect of IL-18 on the induction of antigen-specific IFN-gamma production by resting T cells is IL-12 dependent and is mediated by up-regulation of the IL-12 receptor beta2 subunit. Eur J Immunol. 2000, 30: 1113-1119. 10.1002/(SICI)1521-4141(200004)30:4<1113::AID-IMMU1113>3.0.CO;2-P.View ArticlePubMedGoogle Scholar
- Akira S: The role of IL-18 in innate immunity. Curr Opin Immunol. 2000, 12: 59-63. 10.1016/S0952-7915(99)00051-5.View ArticlePubMedGoogle Scholar
- Kanakaraj P, Ngo K, Wu Y, Angulo A, Ghazal P, Harris CA, Siekierka JJ, Peterson PA, Fung-Leung WP: Defective interleukin (IL)-18-mediated natural killer and T helper cell type 1 response in IL-1 receptor-associated kinase (IRAK)-deficient mice. J Exp Med. 1999, 189: 1129-1138. 10.1084/jem.189.7.1129.PubMed CentralView ArticlePubMedGoogle Scholar
- Smeltz RB, Chen J, Hu-Li J, Shevach EM: Regulation of interleukin (IL)-18 receptor α chain expression on CD4+ T cells during T helper (Th)1/Th2 differentiation: critical downregulatory role of IL-4. J Exp Med. 2001, 194: 143-153. 10.1084/jem.194.2.143.PubMed CentralView ArticlePubMedGoogle Scholar
- Moore KW, de Waal Malefyt R, Coffman R, O'Garra A: Interleukin-10 and the interleukin-10 receptor. Annu Rev Immunol. 2001, 19: 683-765. 10.1146/annurev.immunol.19.1.683.View ArticlePubMedGoogle Scholar
- Riley JK, Takeda K, Akira S, Schreiber RD: Interleukin-10 receptor signaling through the JAK-STAT pathway. J Biol Chem. 1999, 274: 16513-16521. 10.1074/jbc.274.23.16513.View ArticlePubMedGoogle Scholar
- Schmidt-Weber CB, Alexander SI, Henault LE, James L, Lichtman AH: IL-4 enhances IL-10 gene expression in murine Th2 cells in the absence of TCR engagement. J Immunol. 1999, 162: 238-244.PubMedGoogle Scholar
- Skapenko A, Niedobitek GU, Kalden JR, Lipsky PE, Schulze-Koops H: Generation and regulation of human Th1-biased immune response in vivo: A critical role for IL-4 and IL-10. J Immunol. 2004, 172: 6427-6434.View ArticlePubMedGoogle Scholar
- Huang Y, Wange RL: T cell receptor signaling: beyond complex complexes. J Biol Chem. 2004, 279: 28827-28830. 10.1074/jbc.R400012200.View ArticlePubMedGoogle Scholar
- Porter CM, Clipstone NA: Sustained NFAT signaling promotes a Th1-like pattern of gene expression in primary murine CD4^{+} T cells. J Immunol. 2002, 168: 4936-4945.View ArticlePubMedGoogle Scholar
- Kauffman SA: Metabolic stability and epigenesis in randomly constructed genetic nets. J Theor Biol. 1969, 22: 437-467. 10.1016/0022-5193(69)90015-0.View ArticlePubMedGoogle Scholar
- Kauffman SA: Antichaos and adaptation. Sci Am. 1991, 265: 78-84.View ArticlePubMedGoogle Scholar
- Thomas R: Laws for the dynamics of regulatory networks. Int J Dev Biol. 1998, 42: 479-485.PubMedGoogle Scholar
- Muraille E, Leo O: Revisiting the Th1/Th2 paradigm. Scand J Immunol. 1998, 47: 1-9.View ArticlePubMedGoogle Scholar
- Street NE, Mosmann TM: Functional diversity of T lymphocytes due to secretion of different cytokine patterns. FASEB J. 1991, 5: 171-177.PubMedGoogle Scholar
- Murphy KM, Ouyang W, Farrar JD, Yang J, Ranganath S, Asnagli H, Afkarian M, Murphy TL: Signaling and transcription in T helper development. Annu Rev Immunol. 2000, 18: 451-494. 10.1146/annurev.immunol.18.1.451.View ArticlePubMedGoogle Scholar
Copyright
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.