Project description:There is a vast amount of molecular information regarding the differentiation of T lymphocytes, in particular regarding in vitro experimental treatments that modify their differentiation process. This publicly available information was used to infer the regulatory network that controls the differentiation of T lymphocytes into CD4+ and CD8+ cells. Hereby we present a network that consists of 50 nodes and 97 regulatory interactions, representing the main signaling circuits established among molecules and molecular complexes regulating the differentiation of T cells. The network was converted into a continuous dynamical system in the form of a set of coupled ordinary differential equations, and its dynamical behavior was studied. With the aid of numerical methods, nine fixed point attractors were found for the dynamical system. These attractors correspond to the activation patterns observed experimentally for the following cell types: CD4?CD8?, CD4+CD8+, CD4+ naive, Th1, Th2, Th17, Treg, CD8+ naive, and CTL. Furthermore, the model is able to describe the differentiation process from the precursor CD4?CD8? to any of the effector types due to a specific series of extracellular signals.
Project description:Systems models of biological networks show promise for informing drug target selection/qualification, identifying lead compounds and factors regulating disease progression, rationalizing combinatorial regimens, and explaining sources of intersubject variability and adverse drug reactions. However, most models of biological systems are qualitative and are not easily coupled with dynamical models of drug exposure-response relationships. In this proof-of-concept study, logic-based modeling of signal transduction pathways in U266 multiple myeloma (MM) cells is used to guide the development of a simple dynamical model linking bortezomib exposure to cellular outcomes. Bortezomib is a commonly used first-line agent in MM treatment; however, knowledge of the signal transduction pathways regulating bortezomib-mediated cell cytotoxicity is incomplete. A Boolean network model of 66 nodes was constructed that includes major survival and apoptotic pathways and was updated using responses to several chemical probes. Simulated responses to bortezomib were in good agreement with experimental data, and a reduction algorithm was used to identify key signaling proteins. Bortezomib-mediated apoptosis was not associated with suppression of nuclear factor B (NFB) protein inhibition in this cell line, which contradicts a major hypothesis of bortezomib pharmacodynamics. A pharmacodynamic model was developed that included three critical proteins (phospho-NFB, BclxL, and cleaved poly (ADP ribose) polymerase). Model-fitted protein dynamics and cell proliferation profiles agreed with experimental data, and the model-predicted IC50 (3.5 nM) is comparable to the experimental value (1.5 nM). The cell-based pharmacodynamic model successfully links bortezomib exposure to MM cellular proliferation via protein dynamics, and this model may show utility in exploring bortezomib-based combination regimens.
Project description:Due to the large number of diseases associated to a malfunction of the hematopoietic system, there is an interest in knowing the molecular mechanisms controlling the differentiation of blood cell lineages. However, the structure and dynamical properties of the underlying regulatory network controlling this process is not well understood. This manuscript presents a regulatory network of 81 nodes, representing several types of molecules that regulate each other during the process of lymphopoiesis. The regulatory interactions were inferred mostly from published experimental data. However, 15 out of 159 regulatory interactions are predictions arising from the present study. The network is modelled as a continuous dynamical system, in the form of a set of differential equations. The dynamical behaviour of the model describes the differentiation process from the common lymphocyte precursor (CLP) to several mature B and T cell types; namely, plasma cell (PC), cytotoxic T lymphocyte (CTL), T helper 1 (Th1), Th2, Th17, and T regulatory (Treg) cells. The model qualitatively recapitulates key cellular differentiation events, being able to represent the directional and branched nature of lymphopoiesis, going from a multipotent progenitor to fully differentiated cell types.
Project description:Co-infections alter the host immune response but how the systemic and local processes at the site of infection interact is still unclear. The majority of studies on co-infections concentrate on one of the infecting species, an immune function or group of cells and often focus on the initial phase of the infection. Here, we used a combination of experiments and mathematical modelling to investigate the network of immune responses against single and co-infections with the respiratory bacterium Bordetella bronchiseptica and the gastrointestinal helminth Trichostrongylus retortaeformis. Our goal was to identify representative mediators and functions that could capture the essence of the host immune response as a whole, and to assess how their relative contribution dynamically changed over time and between single and co-infected individuals. Network-based discrete dynamic models of single infections were built using current knowledge of bacterial and helminth immunology; the two single infection models were combined into a co-infection model that was then verified by our empirical findings. Simulations showed that a T helper cell mediated antibody and neutrophil response led to phagocytosis and clearance of B. bronchiseptica from the lungs. This was consistent in single and co-infection with no significant delay induced by the helminth. In contrast, T. retortaeformis intensity decreased faster when co-infected with the bacterium. Simulations suggested that the robust recruitment of neutrophils in the co-infection, added to the activation of IgG and eosinophil driven reduction of larvae, which also played an important role in single infection, contributed to this fast clearance. Perturbation analysis of the models, through the knockout of individual nodes (immune cells), identified the cells critical to parasite persistence and clearance both in single and co-infections. Our integrated approach captured the within-host immuno-dynamics of bacteria-helminth infection and identified key components that can be crucial for explaining individual variability between single and co-infections in natural populations.
Project description:Analyzing the long-term behaviors (attractors) of dynamic models of biological systems can provide valuable insight into biological phenotypes and their stability. In this paper we identify the allowed long-term behaviors of a multi-level, 70-node dynamic model of the stomatal opening process in plants. We start by reducing the models huge state space. We first reduce unregulated nodes and simple mediator nodes, then simplify the regulatory functions of selected nodes while keeping the model consistent with experimental observations. We perform attractor analysis on the resulting 32-node reduced model by two methods: 1. converting it into a Boolean model, then applying two attractor-finding algorithms; 2. theoretical analysis of the regulatory functions. We further demonstrate the robustness of signal propagation by showing that a large percentage of single-node knockouts does not affect the stomatal opening level. Combining both methods with analysis of perturbation scenarios, we conclude that all nodes except two in the reduced model have a single attractor; and only two nodes can admit oscillations. The multistability or oscillations of these four nodes do not affect the stomatal opening level in any situation. This conclusion applies to the original model as well in all the biologically meaningful cases. In addition, the stomatal opening level is resilient against single-node knockouts. Thus, we conclude that the complex structure of this signal transduction network provides multiple information propagation pathways while not allowing extensive multistability or oscillations, resulting in robust signal propagation. Our innovative combination of methods offers a promising way to analyze multi-level models.
Project description:Two types of distinct cardiac progenitor cell populations can be identified during early heart development: the first heart field (FHF) and second heart field (SHF) lineage that later form the mature heart. They can be characterized by differential expression of transcription and signaling factors. These regulatory factors influence each other forming a gene regulatory network. Here, we present a core gene regulatory network for early cardiac development based on published temporal and spatial expression data of genes and their interactions. This gene regulatory network was implemented in a Boolean computational model. Simulations reveal stable states within the network model, which correspond to the regulatory states of the FHF and the SHF lineages. Furthermore, we are able to reproduce the expected temporal expression patterns of early cardiac factors mimicking developmental progression. Additionally, simulations of knock-down experiments within our model resemble published phenotypes of mutant mice. Consequently, this gene regulatory network retraces the early steps and requirements of cardiogenic mesoderm determination in a way appropriate to enhance the understanding of heart development.