Project description:Cells and tissues are exposed to stress from numerous sources. Senescence is a protective mechanism that prevents malignant tissue changes and constitutes a fundamental mechanism of aging. It can be accompanied by a senescence associated secretory phenotype (SASP) that causes chronic inflammation. We present a Boolean network model-based gene regulatory network of the SASP, incorporating published gene interaction data. The simulation results describe current biological knowledge. The model predicts different in-silico knockouts that prevent key SASP-mediators, IL-6 and IL-8, from getting activated upon DNA damage. The NF-B Essential Modulator (NEMO) was the most promising in-silico knockout candidate and we were able to show its importance in the inhibition of IL-6 and IL-8 following DNA-damage in murine dermal fibroblasts in-vitro. We strengthen the speculated regulator function of the NF-B signaling pathway in the onset and maintenance of the SASP using in-silico and in-vitro approaches. We were able to mechanistically show, that DNA damage mediated SASP triggering of IL-6 and IL-8 is mainly relayed through NF-B, giving access to possible therapy targets for SASP-accompanied diseases.
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:The replication and life cycle of the influenza virus is governed by an intricate network of intracellular regulatory events during infection, including interactions with an even more complex system of biochemical interactions of the host cell. Computational modeling and systems biology have been successfully employed to further the understanding of various biological systems, however, computational studies of the complexity of intracellular interactions during influenza infection is lacking. In this work, we present the first large-scale dynamical model of the infection and replication cycle of influenza, as well as some of its interactions with the host's signaling machinery. Specifically, we focus on and visualize the dynamics of the internalization and endocytosis of the virus, replication and translation of its genomic components, as well as the assembly of progeny virions. Simulations and analyses of the models dynamics qualitatively reproduced numerous biological phenomena discovered in the laboratory. Finally, comparisons of the dynamics of existing and proposed drugs, our results suggest that a drug targeting PB1:PA would be more efficient than existing Amantadin/Rimantaine or Zanamivir/Oseltamivir.
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:The TOL system encoded by plasmid pWW0 of Pseudomonas putida mt-2 is able to sense a large number of both exogenous and endogenous signals as inputs for the genetic and metabolic circuit that determines the biodegradation of m-xylene. However, whether the enormous combinatorial space of inputs is translated into an equally variable response landscape or is processed into very few outcomes remains unclear. To address this question, we set out to define the number of states that can be obtained by a network of a given set of genes under the control of a specified regulatory circuit that is exposed to all possible combinations of inputs. To this end, the TOL network and its regulatory wiring were formalized as a synchronous logic Boolean circuit that had 10 signals (i.e. pathway substrates, temperature, sugars, amino acids, metabolic regimes and global regulators) as possible inputs. The analysis of the attractors of the circuit using a satisfiability (SAT) algorithm revealed that only eight transcriptional states are reached in response to the 1024 possible combinations of stimuli. The experimental validation resulted in a refinement of the model through the addition of a previously unknown interaction that controls the meta catabolic pathway. The full induction of the two xyl operons occurred with only 1.6% of the input combinations, which suggests that the architecture of the network allows the expression of the xyl genes only under a very narrow range of conditions. These data not only explain much of the unusual layout of the TOL circuit but also strengthen the view of the regulatory circuits of environmental bacteria as digital decision-making devices.
Project description:The blood cancer T cell large granular lymphocyte (T-LGL) leukemia is a chronic disease characterized by a clonal proliferation of cytotoxic T cells. As no curative therapy is yet known for this disease, identification of potential therapeutic targets is of immense importance. In this paper, we perform a comprehensive dynamical and structural analysis of a network model of this disease. By employing a network reduction technique, we identify the stationary states (fixed points) of the system, representing normal and diseased (T-LGL) behavior, and analyze their precursor states (basins of attraction) using an asynchronous Boolean dynamic framework. This analysis identifies the T-LGL states of 54 components of the network, out of which 36 (67%) are corroborated by previous experimental evidence and the rest are novel predictions. We further test and validate one of these newly identified states experimentally. Specifically, we verify the prediction that the node SMAD is over-active in leukemic T-LGL by demonstrating the predominant phosphorylation of the SMAD family members Smad2 and Smad3. Our systematic perturbation analysis using dynamical and structural methods leads to the identification of 19 potential therapeutic targets, 68% of which are corroborated by experimental evidence. The novel therapeutic targets provide valuable guidance for wet-bench experiments. In addition, we successfully identify two new candidates for engineering long-lived T cells necessary for the delivery of virus and cancer vaccines. Overall, this study provides a bird's-eye-view of the avenues available for identification of therapeutic targets for similar diseases through perturbation of the underlying signal transduction network.
Project description:This model represents a set of generic, commonly expressed receptor signaling pathways, including EGFR, G-protein-coupled receptor (alpha i, alpha q, alpha 12/13, and alpha s ligands), integrin, and stress pathways.
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.
Project description:Background In the yeast Saccharomyces cerevisiae, the presence of high levels of glucose leads to an array of down-regulatory effects known as glucose repression. This process is complex due to the presence of feedback loops and crosstalk between different pathways, complicating the use of intuitive approaches to analyze the system. Results We established a logical model of yeast glucose repression, formalized as a hypergraph. The model was constructed based on verified regulatory interactions and it includes 50 gene transcripts, 22 proteins, 5 metabolites and 118 hyperedges. We computed the logical steady states of all nodes in the network in order to simulate wildtype and deletion mutant responses to different sugar availabilities. Evaluation of the model predictive power was achieved by comparing changes in the logical state of gene nodes with transcriptome data. Overall, we observed 71% true predictions, and analyzed sources of errors and discrepancies for the remaining. Conclusion Though the binary nature of logical (Boolean) models entails inherent limitations, our model constitutes a primary tool for storing regulatory knowledge, searching for incoherencies in hypotheses and evaluating the effect of deleting regulatory elements involved in glucose repression.