Stochastic regulation in early immune response

Trade-off between photosymbiosis and innate immunity influences cnidarian's response to pathogenic bacteria

Mutualistic relationships with photosynthetic organisms are common in cnidarians, which form an intracellular symbiosis with dinoflagellates in the family Symbiodiniaceae. The establishment and maintenance of these symbionts are associated with the suppression of key host immune factors. Because of this, there are potential trade-offs between the nutrition that cnidarian hosts gain from their symbionts and their ability to successfully defend themselves from pathogens. To investigate these potential trade-offs, we utilized the facultatively symbiotic polyps of the upside-down jellyfish Cassiopea xamachana and exposed aposymbiotic and symbiotic polyps to the pathogen Serratia marcescens. Symbiotic polyps had a lower probability of survival following S. marcescens exposure. Gene expression analyses 24 hours following pathogen exposure indicate that symbiotic animals mounted a more damaging immune response, with higher levels of inflammation and oxidative stress likely resulting in more severe disruptions to cellular homeostasis. Underlying this more damaging immune response may be differences in constitutive and pathogen-induced expression of immune transcription factors between aposymbiotic and symbiotic polyps rather than broadscale immune suppression during symbiosis. Our findings indicate that in facultatively symbiotic polyps, hosting symbionts limits C. xamachana's ability to survive pathogen exposure, indicating a trade-off between symbiosis and immunity that has potential implications for coral disease research.

Synergistic innate immune activation and anti-tumor immunity through combined STING and TLR4 stimulation

Previous work has shown that innate immune sensing of tumors involves the host STING pathway, which leads to IFN-β production, dendritic cell (DC) activation, and T cell priming against tumor antigens. This observation has led to the development of STING agonists as a potential cancer therapeutic. However, despite promising results in mouse studies using transplantable tumor models, clinical testing of STING agonists has shown activity in only a minority of patients. Thus, further study of innate immune pathways in anti-tumor immunity is paramount. Innate immune activation in response to a pathogen rarely occurs through stimulation of only one signaling pathway, and activating multiple innate immune pathways similar to a natural infection is one possible strategy to improve the efficacy of STING agonists. To test this, we performed experiments with the STING agonist DMXAA alone or in combination with several TLR agonists. We found that LPS + DMXAA induced significantly greater IFN-β transcription than the sum of either agonist alone. To explain this synergy, we assayed each step of STING pathway signaling. LPS did not increase STING protein aggregation, IRF3 phosphorylation, or IRF3 nuclear translocation beyond what occurred with DMXAA alone. However, since the IFN-β promoter also includes NF-κB binding sites, we additionally examined the NF-κB pathway. In fact, LPS increased the phosphorylation and nuclear translocation of the NF-κB subunit p65, and NF-κB signaling was required for the observed synergy. Intratumoral injection of suboptimal doses of LPS + DMXAA resulted in significantly improved tumor control of B16 melanoma in vivo compared to either agonist alone. Our results suggest that combinatorial signaling through TLR4 and STING results in optimal innate signaling via co-involvement of NF-κB and IRF3, and that combined engagement of these two pathways has therapeutic potential.

Endothelial junctional membrane protrusions serve as hotspots for neutrophil transmigration

Upon inflammation, leukocytes rapidly transmigrate across the endothelium to enter the inflamed tissue. Evidence accumulates that leukocytes use preferred exit sites, alhough it is not yet clear how these hotspots in the endothelium are defined and how they are recognized by the leukocyte. Using lattice light sheet microscopy, we discovered that leukocytes prefer endothelial membrane protrusions at cell junctions for transmigration. Phenotypically, these junctional membrane protrusions are present in an asymmetric manner, meaning that one endothelial cell shows the protrusion and the adjacent one does not. Consequently, leukocytes cross the junction by migrating underneath the protruding endothelial cell. These protrusions depend on Rac1 activity and by using a photo-activatable Rac1 probe, we could artificially generate local exit-sites for leukocytes. Overall, we have discovered a new mechanism that uses local induced junctional membrane protrusions to facilitate/steer the leukocyte escape/exit from inflamed vessel walls.

Novel Technologies in Studying Brain Immune Response

Over the past few decades, the immune system, including both the adaptive and innate immune systems, proved to be essential and critical to brain damage and recovery in the pathogenesis of several diseases, opening a new avenue for developing new immunomodulatory therapies and novel treatments for many neurological diseases. However, due to the specificity and structural complexity of the central nervous system (CNS), and the limit of the related technologies, the biology of the immune response in the brain is still poorly understood. Here, we discuss the application of novel technologies in studying the brain immune response, including single-cell RNA analysis, cytometry by time-of-flight, and whole-genome transcriptomic and proteomic analysis. We believe that advancements in technology related to immune research will provide an optimistic future for brain repair.

Of Keeping and Tipping the Balance: Host Regulation and Viral Modulation of IRF3-Dependent IFNB1 Expression

The type I interferon (IFN) response is a principal component of our immune system that allows to counter a viral attack immediately upon viral entry into host cells. Upon engagement of aberrantly localised nucleic acids, germline-encoded pattern recognition receptors convey their find via a signalling cascade to prompt kinase-mediated activation of a specific set of five transcription factors. Within the nucleus, the coordinated interaction of these dimeric transcription factors with coactivators and the basal RNA transcription machinery is required to access the gene encoding the type I IFN IFNβ (IFNB1). Virus-induced release of IFNβ then induces the antiviral state of the system and mediates further mechanisms for defence. Due to its key role during the induction of the initial IFN response, the activity of the transcription factor interferon regulatory factor 3 (IRF3) is tightly regulated by the host and fiercely targeted by viral proteins at all conceivable levels. In this review, we will revisit the steps enabling the trans-activating potential of IRF3 after its activation and the subsequent assembly of the multi-protein complex at the IFNβ enhancer that controls gene expression. Further, we will inspect the regulatory mechanisms of these steps imposed by the host cell and present the manifold strategies viruses have evolved to intervene with IFNβ transcription downstream of IRF3 activation in order to secure establishment of a productive infection.

The FLAME-accelerated signalling tool (FaST) for facile parallelisation of flexible agent-based models of cell signalling

Agent-based modelling is particularly adept at modelling complex features of cell signalling pathways, where heterogeneity, stochastic and spatial effects are important, thus increasing our understanding of decision processes in biology in such scenarios. However, agent-based modelling often is computationally prohibitive to implement. Parallel computing, either on central processing units (CPUs) or graphical processing units (GPUs), can provide a means to improve computational feasibility of agent-based applications but generally requires specialist coding knowledge and extensive optimisation. In this paper, we address these challenges through the development and implementation of the FLAME-accelerated signalling tool (FaST), a software that permits easy creation and parallelisation of agent-based models of cell signalling, on CPUs or GPUs. FaST incorporates validated new agent-based methods, for accurate modelling of reaction kinetics and, as proof of concept, successfully converted an ordinary differential equation (ODE) model of apoptosis execution into an agent-based model. We finally parallelised this model through FaST on CPUs and GPUs resulting in an increase in performance of 5.8× (16 CPUs) and 53.9×, respectively. The FaST takes advantage of the communicating X-machine approach used by FLAME and FLAME GPU to allow easy alteration or addition of functionality to parallel applications, but still includes inherent parallelisation optimisation. The FaST, therefore, represents a new and innovative tool to easily create and parallelise bespoke, robust, agent-based models of cell signalling.

Macrophages employ quorum licensing to regulate collective activation

Macrophage-initiated inflammation is tightly regulated to eliminate threats such as infections while suppressing harmful immune activation. However, individual cells’ signaling responses to pro-inflammatory cues are heterogeneous, with subpopulations emerging with high or low activation states. Here, we use single-cell tracking and dynamical modeling to develop and validate a revised model for lipopolysaccharide (LPS)-induced macrophage activation that invokes a mechanism we term quorum licensing. The results show that bimodal phenotypic partitioning of macrophages is primed during the resting state, dependent on cumulative history of cell density, predicted by extrinsic noise in transcription factor expression, and independent of canonical LPS-induced intercellular feedback in the tumor necrosis factor (TNF) response. Our analysis shows how this density-dependent coupling produces a nonlinear effect on collective TNF production. We speculate that by linking macrophage density to activation, this mechanism could amplify local responses to threats and prevent false alarms.

Macrophage activation is tightly regulated to maintain immune homeostasis, yet activation is also heterogeneous. Here, the authors show that macrophages coordinate activation by partitioning into two phenotypes that can nonlinearly amplify collective inflammatory cytokine production as a function of cell density.

Sensitivity analysis methods in the biomedical sciences

Sensitivity analysis is an important part of a mathematical modeller's toolbox for model analysis. In this review paper, we describe the most frequently used sensitivity techniques, discussing their advantages and limitations, before applying each method to a simple model. Also included is a summary of current software packages, as well as a modeller's guide for carrying out sensitivity analyses. Finally, we apply the popular Morris and Sobol methods to two models with biomedical applications, with the intention of providing a deeper understanding behind both the principles of these methods and the presentation of their results.

An ecosystem framework for understanding and treating disease

Pathogens and cancers are pervasive health risks in the human population. I argue that if we are to better understand disease and its treatment, then we need to take an ecological perspective of disease itself. I generalize and extend an emerging framework that views disease as an ecosystem and many of its components as interacting in a community. I develop the framework for biological etiological agents (BEAs) that multiply within humans—focusing on bacterial pathogens and cancers—but the framework could be extended to include other host and parasite species. I begin by describing why we need an ecosystem framework to understand disease, and the main components and interactions in bacterial and cancer disease ecosystems. Focus is then given to the BEA and how it may proceed through characteristic states, including emergence, growth, spread and regression. The framework is then applied to therapeutic interventions. Central to success is preventing BEA evasion, the best known being antibiotic resistance and chemotherapeutic resistance in cancers. With risks of evasion in mind, I propose six measures that either introduce new components into the disease ecosystem or manipulate existing ones. An ecosystem framework promises to enhance our understanding of disease, BEA and host (co)evolution, and how we can improve therapeutic outcomes.

Quantitative analysis reveals crosstalk mechanisms of heat shock-induced attenuation of NF-κB signaling at the single cell level

Elevated temperature induces the heat shock (HS) response, which modulates cell proliferation, apoptosis, the immune and inflammatory responses. However, specific mechanisms linking the HS response pathways to major cellular signaling systems are not fully understood. Here we used integrated computational and experimental approaches to quantitatively analyze the crosstalk mechanisms between the HS-response and a master regulator of inflammation, cell proliferation, and apoptosis the Nuclear Factor κB (NF-κB) system. We found that populations of human osteosarcoma cells, exposed to a clinically relevant 43°C HS had an attenuated NF-κB p65 response to Tumor Necrosis Factor α (TNFα) treatment. The degree of inhibition of the NF-κB response depended on the HS exposure time. Mathematical modeling of single cells indicated that individual crosstalk mechanisms differentially encode HS-mediated NF-κB responses while being consistent with the observed population-level responses. In particular "all-or-nothing" encoding mechanisms were involved in the HS-dependent regulation of the IKK activity and IκBα phosphorylation, while others involving transport were "analogue". In order to discriminate between these mechanisms, we used live-cell imaging of nuclear translocations of the NF-κB p65 subunit. The single cell responses exhibited "all-or-nothing" encoding. While most cells did not respond to TNFα stimulation after a 60 min HS, 27% showed responses similar to those not receiving HS. We further demonstrated experimentally and theoretically that the predicted inhibition of IKK activity was consistent with the observed HS-dependent depletion of the IKKα and IKKβ subunits in whole cell lysates. However, a combination of "all-or-nothing" crosstalk mechanisms was required to completely recapitulate the single cell data. We postulate therefore that the heterogeneity of the single cell responses might be explained by the cell-intrinsic variability of HS-modulated IKK signaling. In summary, we show that high temperature modulates NF-κB responses in single cells in a complex and unintuitive manner, which needs to be considered in hyperthermia-based treatment strategies.

Modeling and Hemofiltration Treatment of Acute Inflammation

The body responds to endotoxins by triggering the acute inflammatory response system to eliminate the threat posed by gram-negative bacteria (endotoxin) and restore health. However, an uncontrolled inflammatory response can lead to tissue damage, organ failure, and ultimately death; this is clinically known as sepsis. Mathematical models of acute inflammatory disease have the potential to guide treatment decisions in critically ill patients. In this work, an 8-state (8-D) differential equation model of the acute inflammatory response system to endotoxin challenge was developed. Endotoxin challenges at 3 and 12 mg/kg were administered to rats, and dynamic cytokine data for interleukin (IL)-6, tumor necrosis factor (TNF), and IL-10 were obtained and used to calibrate the model. Evaluation of competing model structures was performed by analyzing model predictions at 3, 6, and 12 mg/kg endotoxin challenges with respect to experimental data from rats. Subsequently, a model predictive control (MPC) algorithm was synthesized to control a hemoadsorption (HA) device, a blood purification treatment for acute inflammation. A particle filter (PF) algorithm was implemented to estimate the full state vector of the endotoxemic rat based on time series cytokine measurements. Treatment simulations show that: (i) the apparent primary mechanism of HA efficacy is white blood cell (WBC) capture, with cytokine capture a secondary benefit; and (ii) differential filtering of cytokines and WBC does not provide substantial improvement in treatment outcomes vs. existing HA devices.

Signal transduction controls heterogeneous NF-κB dynamics and target gene expression through cytokine-specific refractory states

Cells respond dynamically to pulsatile cytokine stimulation. Here we report that single, or well-spaced pulses of TNFα (>100 min apart) give a high probability of NF-κB activation. However, fewer cells respond to shorter pulse intervals (<100 min) suggesting a heterogeneous refractory state. This refractory state is established in the signal transduction network downstream of TNFR and upstream of IKK, and depends on the level of the NF-κB system negative feedback protein A20. If a second pulse within the refractory phase is IL-1β instead of TNFα, all of the cells respond. This suggests a mechanism by which two cytokines can synergistically activate an inflammatory response. Gene expression analyses show strong correlation between the cellular dynamic response and NF-κB-dependent target gene activation. These data suggest that refractory states in the NF-κB system constitute an inherent design motif of the inflammatory response and we suggest that this may avoid harmful homogenous cellular activation.

Regulation of signal transduction by spatial parameters: a case in NF-κB oscillation

NF-κB is a transcription factor regulating expression of more than 500 genes, and its dysfunction leads to the autoimmune and inflammatory diseases. In malignant cancer cells, NF-κB is constitutively activated. Thus the elucidation of mechanisms for NF-κB regulation is important for the establishment of therapeutic treatment caused by incorrect NF-κB responses. Cytoplasmic NF-κB translocates to the nucleus by the application of extracellular stimuli such as cytokines. Nuclear NF-κB is known to oscillate with the cycle of 1.5-4.5 h, and it is thought that the oscillation pattern regulates the expression profiles of genes. In this review, first we briefly describe regulation mechanisms of NF-κB. Next, published computational simulations on the oscillation of NF-κB are summarised. There are at least 60 reports on the computational simulation and analysis of NF-κB oscillation. Third, the importance of a 'space' for the regulation of oscillation pattern of NF-κB is discussed, showing altered oscillation pattern by the change in spatial parameters such as diffusion coefficient, nuclear to cytoplasmic volume ratio (N/C ratio), and transport through nuclear membrane. Finally, simulations in a true intracellular space (TiCS), which is an intracellular 3D space reconstructed in a computer with organelles such as nucleus and mitochondria are discussed.

NF-κB oscillations translate into functionally related patterns of gene expression

Several transcription factors (TFs) oscillate, periodically relocating between the cytoplasm and the nucleus. NF-κB, which plays key roles in inflammation and cancer, displays oscillations whose biological advantage remains unclear. Recent work indicated that NF-κB displays sustained oscillations that can be entrained, that is, reach a persistent synchronized state through small periodic perturbations. We show here that for our GFP-p65 knock-in cells NF-κB behaves as a damped oscillator able to synchronize to a variety of periodic external perturbations with no memory. We imposed synchronous dynamics to prove that transcription of NF-κB-controlled genes also oscillates, but mature transcript levels follow three distinct patterns. Two sets of transcripts accumulate fast or slowly, respectively. Another set, comprising chemokine and chemokine receptor mRNAs, oscillates and resets at each new stimulus, with no memory of the past. We propose that TF oscillatory dynamics is a means of segmenting time to provide renewing opportunity windows for decision.

A Neutrophil Phenotype Model for Extracorporeal Treatment of Sepsis

Neutrophils play a central role in eliminating bacterial pathogens, but may also contribute to end-organ damage in sepsis. Interleukin-8 (IL-8), a key modulator of neutrophil function, signals through neutrophil specific surface receptors CXCR-1 and CXCR-2. In this study a mechanistic computational model was used to evaluate and deploy an extracorporeal sepsis treatment which modulates CXCR-1/2 levels. First, a simplified mechanistic computational model of IL-8 mediated activation of CXCR-1/2 receptors was developed, containing 16 ODEs and 43 parameters. Receptor level dynamics and systemic parameters were coupled with multiple neutrophil phenotypes to generate dynamic populations of activated neutrophils which reduce pathogen load, and/or primed neutrophils which cause adverse tissue damage when misdirected. The mathematical model was calibrated using experimental data from baboons administered a two-hour infusion of E coli and followed for a maximum of 28 days. Ensembles of parameters were generated using a Bayesian parallel tempering approach to produce model fits that could recreate experimental outcomes. Stepwise logistic regression identified seven model parameters as key determinants of mortality. Sensitivity analysis showed that parameters controlling the level of killer cell neutrophils affected the overall systemic damage of individuals. To evaluate rescue strategies and provide probabilistic predictions of their impact on mortality, time of onset, duration, and capture efficacy of an extracorporeal device that modulated neutrophil phenotype were explored. Our findings suggest that interventions aiming to modulate phenotypic composition are time sensitive. When introduced between 3–6 hours of infection for a 72 hour duration, the survivor population increased from 31% to 40–80%. Treatment efficacy quickly diminishes if not introduced within 15 hours of infection. Significant harm is possible with treatment durations ranging from 5–24 hours, which may reduce survival to 13%. In severe sepsis, an extracorporeal treatment which modulates CXCR-1/2 levels has therapeutic potential, but also potential for harm. Further development of the computational model will help guide optimal device development and determine which patient populations should be targeted by treatment.

Sepsis occurs when a patient develops a whole body immune response due to infection. In this condition, white blood cells called neutrophils circulate in an active state, seeking and eliminating invading bacteria. However, when neutrophils are activated, healthy tissue is inadvertently targeted, leading to organ damage and potentially death. Even though sepsis kills millions worldwide, there are still no specific treatments approved in the United States. This may be due to the complexity and diversity of the body’s immune response, which can be managed well using computational modeling. We have developed a computational model to predict how different levels of neutrophil activation impact survival in an overactive inflammatory conditions. The model was utilized to assess the effectiveness of a simulated experimental sepsis treatment which modulates neutrophil populations and activity. This evaluation determined that treatment timing plays a critical role in therapeutic effectiveness. When utilized properly the treatment drastically improves survival, but there is also risk of causing patient harm when introduced at the wrong time. We intend for this computational model to support and guide further development of sepsis treatments and help translate these preliminary results from bench to bedside.

Roles of factorial noise in inducing bimodal gene expression

Some gene regulatory systems can exhibit bimodal distributions of mRNA or protein although the deterministic counterparts are monostable. This noise-induced bimodality is an interesting phenomenon and has important biological implications, but it is unclear how different sources of expression noise (each source creates so-called factorial noise that is defined as a component of the total noise) contribute separately to this stochastic bimodality. Here we consider a minimal model of gene regulation, which is monostable in the deterministic case. Although simple, this system contains factorial noise of two main kinds: promoter noise due to switching between gene states and transcriptional (or translational) noise due to synthesis and degradation of mRNA (or protein). To better trace the roles of factorial noise in inducing bimodality, we also analyze two limit models, continuous and adiabatic approximations, apart from the exact model. We show that in the case of slow gene switching, the continuous model where only promoter noise is considered can exhibit bimodality; in the case of fast switching, the adiabatic model where only transcriptional or translational noise is considered can also exhibit bimodality but the exact model cannot; and in other cases, both promoter noise and transcriptional or translational noise can cooperatively induce bimodality. Since slow gene switching and large protein copy numbers are characteristics of eukaryotic cells, whereas fast gene switching and small protein copy numbers are characteristics of prokaryotic cells, we infer that eukaryotic stochastic bimodality is induced mainly by promoter noise, whereas prokaryotic stochastic bimodality is induced primarily by transcriptional or translational noise.

Modeling dynamics of cell-to-cell variability in TRAIL-induced apoptosis explains fractional killing and predicts reversible resistance

Isogenic cells sensing identical external signals can take markedly different decisions. Such decisions often correlate with pre-existing cell-to-cell differences in protein levels. When not neglected in signal transduction models, these differences are accounted for in a static manner, by assuming randomly distributed initial protein levels. However, this approach ignores the a priori non-trivial interplay between signal transduction and the source of this cell-to-cell variability: temporal fluctuations of protein levels in individual cells, driven by noisy synthesis and degradation. Thus, modeling protein fluctuations, rather than their consequences on the initial population heterogeneity, would set the quantitative analysis of signal transduction on firmer grounds. Adopting this dynamical view on cell-to-cell differences amounts to recast extrinsic variability into intrinsic noise. Here, we propose a generic approach to merge, in a systematic and principled manner, signal transduction models with stochastic protein turnover models. When applied to an established kinetic model of TRAIL-induced apoptosis, our approach markedly increased model prediction capabilities. One obtains a mechanistic explanation of yet-unexplained observations on fractional killing and non-trivial robust predictions of the temporal evolution of cell resistance to TRAIL in HeLa cells. Our results provide an alternative explanation to survival via induction of survival pathways since no TRAIL-induced regulations are needed and suggest that short-lived anti-apoptotic protein Mcl1 exhibit large and rare fluctuations. More generally, our results highlight the importance of accounting for stochastic protein turnover to quantitatively understand signal transduction over extended durations, and imply that fluctuations of short-lived proteins deserve particular attention.

Agent-based model of fecal microbial transplant effect on bile acid metabolism on suppressing Clostridium difficile infection: an example of agent-based modeling of intestinal bacterial infection

Agent-based modeling is a computational modeling method that represents system-level behavior as arising from multiple interactions between the multiple components that make up a system. Biological systems are thus readily described using agent-based models (ABMs), as multi-cellular organisms can be viewed as populations of interacting cells, and microbial systems manifest as colonies of individual microbes. Intersections between these two domains underlie an increasing number of pathophysiological processes, and the intestinal tract represents one of the most significant locations for these inter-domain interactions, so much so that it can be considered an internal ecology of varying robustness and function. Intestinal infections represent significant disturbances of this internal ecology, and one of the most clinically relevant intestinal infections is Clostridium difficile infection (CDI). CDI is precipitated by the use of broad-spectrum antibiotics, involves the depletion of commensal microbiota, and alterations in bile acid composition in the intestinal lumen. We present an example ABM of CDI (the C. difficile Infection ABM, or CDIABM) to examine fundamental dynamics of the pathogenesis of CDI and its response to treatment with anti-CDI antibiotics and a newer treatment therapy, fecal microbial transplant. The CDIABM focuses on one specific mechanism of potential CDI suppression: commensal modulation of bile acid composition. Even given its abstraction, the CDIABM reproduces essential dynamics of CDI and its response to therapy, and identifies a paradoxical zone of behavior that provides insight into the role of intestinal nutritional status and the efficacy of anti-CDI therapies. It is hoped that this use case example of the CDIABM can demonstrate the usefulness of both agent-based modeling and the application of abstract functional representation as the biomedical community seeks to address the challenges of increasingly complex diseases with the goal of personalized medicine.

Spontaneous NF-κB activation by autocrine TNFα signaling: a computational analysis

NF-κB is a key transcription factor that regulates innate immune response. Its activity is tightly controlled by numerous feedback loops, including two negative loops mediated by NF-κB inducible inhibitors, IκBα and A20, which assure oscillatory responses, and by positive feedback loops arising due to the paracrine and autocrine regulation via TNFα, IL-1 and other cytokines. We study the NF-κB system of interlinked negative and positive feedback loops, combining bifurcation analysis of the deterministic approximation with stochastic numerical modeling. Positive feedback assures the existence of limit cycle oscillations in unstimulated wild-type cells and introduces bistability in A20-deficient cells. We demonstrated that cells of significant autocrine potential, i.e., cells characterized by high secretion of TNFα and its receptor TNFR1, may exhibit sustained cytoplasmic-nuclear NF-κB oscillations which start spontaneously due to stochastic fluctuations. In A20-deficient cells even a small TNFα expression rate qualitatively influences system kinetics, leading to long-lasting NF-κB activation in response to a short-pulsed TNFα stimulation. As a consequence, cells with impaired A20 expression or increased TNFα secretion rate are expected to have elevated NF-κB activity even in the absence of stimulation. This may lead to chronic inflammation and promote cancer due to the persistent activation of antiapoptotic genes induced by NF-κB. There is growing evidence that A20 mutations correlate with several types of lymphomas and elevated TNFα secretion is characteristic of many cancers. Interestingly, A20 loss or dysfunction also leaves the organism vulnerable to septic shock and massive apoptosis triggered by the uncontrolled TNFα secretion, which at high levels overcomes the antiapoptotic action of NF-κB. It is thus tempting to speculate that some cancers of deregulated NF-κB signaling may be prone to the pathogen-induced apoptosis.

RECEPTOR-DEPENDENT SENSITIVITY OF NF-κB TO LOW PHYSIOLOGICAL LEVEL

In the NFκB signaling pathway, cells respond to different concentrations of the TNFα signal by means of NFκB transcription factors. Previous studies showed that most cells are activated under high-dose stimulations and NFκB activation results in oscillations in nuclear NFκB abundance. Here, by analyzing sensitivity gain for the response of the nuclear NFκB to the number of cell-surface receptors under low-dose stimulations, we show that changes in the receptor number can give rise to significant changes in the nonsaturation part of the dose–response curve, where the receptor activation rates are very sensitive to stimulations. In addition, the number of the activated receptors tends to increase in a large range of stimulation dose and can significantly influence the expression of the downstream genes. These results imply that the number of cell-surface receptors plays a role of information encoding like frequency or amplitude encoding described in previous studies.

Statistical ensemble analysis for simulating extrinsic noise-driven response in NF-κB signaling networks

Background

Gene expression profiles and protein dynamics in single cells have a large cell-to-cell variability due to intracellular noise. Intracellular fluctuations originate from two sources: intrinsic noise due to the probabilistic nature of biochemical reactions and extrinsic noise due to randomized interactions of the cell with other cellular systems or its environment. Presently, there is no systematic parameterization and modeling scheme to simulate cellular response at the single cell level in the presence of extrinsic noise.

Results

In this paper, we propose a novel statistical ensemble method to simulate the distribution of heterogeneous cellular responses in single cells. We capture the effects of extrinsic noise by randomizing values of the model parameters. In this context, a statistical ensemble is a large number of system replicates, each with randomly sampled model parameters from biologically feasible intervals. We apply this statistical ensemble approach to the well-studied NF-κB signaling system. We predict several characteristic dynamic features of NF-κB response distributions; one of them is the dosage-dependent distribution of the first translocation time of NF-κB.

Conclusion

The distributions of heterogeneous cellular responses that our statistical ensemble formulation generates reveal the effect of different cellular conditions, e.g., effects due to wild type versus mutant cells or between different dosages of external stimulants. Distributions generated in the presence of extrinsic noise yield valuable insight into underlying regulatory mechanisms, which are sometimes otherwise hidden.

An agent-based model of cellular dynamics and circadian variability in human endotoxemia

As cellular variability and circadian rhythmicity play critical roles in immune and inflammatory responses, we present in this study an agent-based model of human endotoxemia to examine the interplay between circadian controls, cellular variability and stochastic dynamics of inflammatory cytokines. The model is qualitatively validated by its ability to reproduce circadian dynamics of inflammatory mediators and critical inflammatory responses after endotoxin administration in vivo. Novel computational concepts are proposed to characterize the cellular variability and synchronization of inflammatory cytokines in a population of heterogeneous leukocytes. Our results suggest that there is a decrease in cell-to-cell variability of inflammatory cytokines while their synchronization is increased after endotoxin challenge. Model parameters that are responsible for IκB production stimulated by NFκB activation and for the production of anti-inflammatory cytokines have large impacts on system behaviors. Additionally, examining time-dependent systemic responses revealed that the system is least vulnerable to endotoxin in the early morning and most vulnerable around midnight. Although much remains to be explored, proposed computational concepts and the model we have pioneered will provide important insights for future investigations and extensions, especially for single-cell studies to discover how cellular variability contributes to clinical implications.

Heterogeneous kinetics of AKT signaling in individual cells are accounted for by variable protein concentration

In most solid cancers, cells harboring oncogenic mutations represent only a sub-fraction of the entire population. Within this sub-fraction the expression level of mutated proteins can vary significantly due to cellular variability limiting the efficiency of targeted therapy. To address the causes of the heterogeneity, we performed a systematic analysis of one of the most frequently mutated pathways in cancer cells, the phosphatidylinositol 3 kinase (PI3K) signaling pathway. Among others PI3K signaling is activated by the hepatocyte growth factor (HGF) that regulates proliferation of hepatocytes during liver regeneration but also fosters tumor cell proliferation. HGF-mediated responses of PI3K signaling were monitored both at the single cell and cell population level in primary mouse hepatocytes and in the hepatoma cell line Hepa1_6. Interestingly, we observed that the HGF-mediated AKT responses at the level of individual cells is rather heterogeneous. However, the overall average behavior of the single cells strongly resembled the dynamics of AKT activation determined at the cell population level. To gain insights into the molecular cause for the observed heterogeneous behavior of individual cells, we employed dynamic mathematical modeling in a stochastic framework. Our analysis demonstrated that intrinsic noise was not sufficient to explain the observed kinetic behavior, but rather the importance of extrinsic noise has to be considered. Thus, distinct from gene expression in the examined signaling pathway fluctuations of the reaction rates has only a minor impact whereas variability in the concentration of the various signaling components even in a clonal cell population is a key determinant for the kinetic behavior.

Roles of spatial parameters on the oscillation of nuclear NF-κB: computer simulations of a 3D spherical cell

Transcription factor NF-κB resides in the cytoplasm and translocates to the nucleus by application of extracellular stimuli. It is known that the nuclear NF-κB oscillates and different oscillation patterns lead to different gene expression. Nearly forty reports on modeling and simulation of nuclear NF-κB have been published to date. The computational models reported so far are temporal or two-dimensional, and the discussions on spatial parameters have not been involved or limited. Since spatial parameters in cancer cells such as nuclear to cytoplasmic volume (N/C) ratio are different from normal cells, it is important to understand the relationship between oscillation patterns and spatial parameters. Here we report simulations of a 3D computational model for the oscillation of nuclear NF-κB using A-Cell software. First, we found that the default biochemical kinetic constants used in the temporal model cannot replicate the experimentally observed oscillation in the 3D model. Thus, the default parameters should be changed in the 3D model. Second, spatial parameters such as N/C ratio, nuclear transport, diffusion coefficients, and the location of IκB synthesis were found to alter the oscillation pattern. Third, among them, larger N/C ratios resulted in persistent oscillation of nuclear NF-κB, and larger nuclear transport resulted in faster oscillation frequency. Our simulation results suggest that the changes in spatial parameters seen in cancer cells is one possible mechanism for alteration in the oscillation pattern of nuclear NF-κB and lead to the altered gene expression in these cells.

Lessons from mathematically modeling the NF-κB pathway

Mathematical modeling has proved to be a critically important approach in the study of many complex networks and dynamic systems in physics, engineering, chemistry, and biology. The nuclear factor κB (NF-κB) system consists of more than 50 proteins and protein complexes and is both a highly networked and dynamic system. To date, mathematical modeling has only addressed a small fraction of the molecular species and their regulation, but when employed in conjunction with experimental analysis has already led to important insights. Here, we provide a personal account of studying how the NF-κB signaling system functions using mathematical descriptions of the molecular mechanisms. We focus on the insights gained about some of the key regulatory components: the control of the steady state, the signaling dynamics, and signaling crosstalk. We also discuss the biological relevance of these regulatory systems properties.

Agent-based dynamic knowledge representation of Pseudomonas aeruginosa virulence activation in the stressed gut: Towards characterizing host-pathogen interactions in gut-derived sepsis

Background

There is a growing realization that alterations in host-pathogen interactions (HPI) can generate disease phenotypes without pathogen invasion. The gut represents a prime region where such HPI can arise and manifest. Under normal conditions intestinal microbial communities maintain a stable, mutually beneficial ecosystem. However, host stress can lead to changes in environmental conditions that shift the nature of the host-microbe dialogue, resulting in escalation of virulence expression, immune activation and ultimately systemic disease. Effective modulation of these dynamics requires the ability to characterize the complexity of the HPI, and dynamic computational modeling can aid in this task. Agent-based modeling is a computational method that is suited to representing spatially diverse, dynamical systems. We propose that dynamic knowledge representation of gut HPI with agent-based modeling will aid in the investigation of the pathogenesis of gut-derived sepsis.

Methodology/Principal Findings

An agent-based model (ABM) of virulence regulation in Pseudomonas aeruginosa was developed by translating bacterial and host cell sense-and-response mechanisms into behavioral rules for computational agents and integrated into a virtual environment representing the host-microbe interface in the gut. The resulting gut milieu ABM (GMABM) was used to: 1) investigate a potential clinically relevant laboratory experimental condition not yet developed - i.e. non-lethal transient segmental intestinal ischemia, 2) examine the sufficiency of existing hypotheses to explain experimental data - i.e. lethality in a model of major surgical insult and stress, and 3) produce behavior to potentially guide future experimental design - i.e. suggested sample points for a potential laboratory model of non-lethal transient intestinal ischemia. Furthermore, hypotheses were generated to explain certain discrepancies between the behaviors of the GMABM and biological experiments, and new investigatory avenues proposed to test those hypotheses.

Conclusions/Significance

Agent-based modeling can account for the spatio-temporal dynamics of an HPI, and, even when carried out with a relatively high degree of abstraction, can be useful in the investigation of system-level consequences of putative mechanisms operating at the individual agent level. We suggest that an integrated and iterative heuristic relationship between computational modeling and more traditional laboratory and clinical investigations, with a focus on identifying useful and sufficient degrees of abstraction, will enhance the efficiency and translational productivity of biomedical research.

Sources of cell-to-cell variability in canonical nuclear factor-κB (NF-κB) signaling pathway inferred from single cell dynamic images

The canonical nuclear factor-κB (NF-κB) signaling pathway controls a gene network important in the cellular inflammatory response. Upon activation, NF-κB/RelA is released from cytoplasmic inhibitors, from where it translocates into the nucleus, subsequently activating negative feedback loops producing either monophasic or damped oscillatory nucleo-cytoplasmic dynamics. Although the population behavior of the NF-κB pathway has been extensively modeled, the sources of cell-to-cell variability are not well understood. We describe an integrated experimental-computational analysis of NF-κB/RelA translocation in a validated cell model exhibiting monophasic dynamics. Quantitative measures of cellular geometry and total cytoplasmic concentration and translocated RelA amounts were used as priors in Bayesian inference to estimate biophysically realistic parameter values based on dynamic live cell imaging studies of enhanced GFP-tagged RelA in stable transfectants. Bayesian inference was performed on multiple cells simultaneously, assuming identical reaction rate parameters, whereas cellular geometry and initial and total NF-κB concentration-related parameters were cell-specific. A subpopulation of cells exhibiting distinct kinetic profiles was identified that corresponded to differences in the IκBα translation rate. We conclude that cellular geometry, initial and total NF-κB concentration, IκBα translation, and IκBα degradation rates account for distinct cell-to-cell differences in canonical NF-κB translocation dynamics.

Crosstalk between p53 and nuclear factor-B systems: pro- and anti-apoptotic functions of NF-B

Nuclear factors p53 and NF-B control many physiological processes including cell cycle arrest, DNA repair, apoptosis, death, innate and adaptive immune responses, and inflammation. There are numerous pathways linking these systems and there is a bulk of evidence for cooperation as well as for antagonisms between p53 and NF-B. In this theoretical study, the authors use earlier models of p53 and NF-B systems and construct a crosstalk model of p53-NF-B network in order to explore the consequences of the two-way coupling, in which NF-B upregulates the transcription of p53, whereas in turn p53 attenuates transcription of NF-B inhibitors IB and A20. We consider a number of protocols in which cells are stimulated by tumour necrosis factor- (TNF) (that activates NF-B pathway) and/or gamma irradiation (that activates p53 pathway). The authors demonstrate that NF-B may have both anti- and pro-apoptotic roles. TNF stimulation, preceding DNA damaging irradiation, makes cells more resistant to irradiation-induced apoptosis, whereas the same TNF stimulation, when preceded by irradiation, increases the apoptotic cell fraction. The finding suggests that diverse roles of NF-B in apoptosis and cancer could be related to the dynamical context of activation of p53 and NF-B pathways. [Includes supplementary material].

Exploring mechanisms of oscillations in p53 and nuclear factor-B systems

A number of regulatory networks have the potential to generate sustained oscillations of irregular amplitude, but well conserved period. Single-cell experiments revealed that in p53 and nuclear factor (NF)-B systems the oscillation period is homogenous in cell populations, insensitive to the strength of the stimulation, and is not influenced by the overexpression of p53 or NF-B transcription factors. We propose a novel computational method of validation of molecular pathways models, based on the analysis of sensitivity of the oscillation period to the particular gene(s) copy number and the level of stimulation. Using this method, the authors demonstrate that existing p53 models, in which oscillations are borne at a saddle-node-on-invariant-circle or subcritical Hopf bifurcations (characteristic for systems with interlinked positive and negative feedbacks), are highly sensitive to gene copy number. Hence, these models cannot explain existing experiments. Analysing NF-B system, the authors show the importance of saturation in transcription of the NF-B inhibitor IB. Models without saturation predict that the oscillation period is a rapidly growing function of total NF-B level, which is in disagreement with experiments. The study supports the hypothesis that oscillations of robust period are based on supercritical Hopf bifurcation, characteristic for dynamical systems involving negative feedback and time delay. We hypothesise that in the p53 system, the role of positive feedback is not sustaining oscillations, but terminating them in severely damaged cells in which the apoptotic programme should be initiated.

Agent-based models in translational systems biology

Effective translational methodologies for knowledge representation are needed in order to make strides against the constellation of diseases that affect the world today. These diseases are defined by their mechanistic complexity, redundancy, and nonlinearity. Translational systems biology aims to harness the power of computational simulation to streamline drug/device design, simulate clinical trials, and eventually to predict the effects of drugs on individuals. The ability of agent-based modeling to encompass multiple scales of biological process as well as spatial considerations, coupled with an intuitive modeling paradigm, suggests that this modeling framework is well suited for translational systems biology. This review describes agent-based modeling and gives examples of its translational applications in the context of acute inflammation and wound healing.

Understanding the temporal codes of intra-cellular signals

The health of organisms and cells depends on appropriate responses to diverse internal and external cues, stimuli, or challenges, such as changes in hormone or cytokine levels, or exposure to a pathogen. Cellular responses must be tailored to the identity and intensity of the stimulus and therefore intra-cellular signals must carry information about both. However, signaling mediators often form intricate networks that react to multiple stimuli yet manage to produce stimulus-specific responses. The multi-functionality ('functional pleiotropism') of signaling nodes suggests that biological networks have evolved ways of passing physiologically relevant stimulus information through shared channels. Increasing evidence supports the notion that this is achieved in part through temporal regulation of signaling mediators' activities. The present challenge is to identify the features of temporal activity profile that represent information about a given stimulus and understand how cells read the temporal codes to control their responses.

Nuclear translocation kinetics of NF-kappaB in macrophages challenged with pathogens in a microfluidic platform

We have developed a microfluidic platform for real-time imaging of host-pathogen interactions and cellular signaling events. Host cells are immobilized in a controlled environment for optical interrogation of the kinetics and stochasticity of immune response to pathogenic challenges. Here, we have quantitatively measured activation of the toll-like receptor 4 (TLR4) pathway in RAW264.7 murine macrophage-like cells. This was achieved by measuring the cytoplasm-to-nucleus translocation kinetics of a green fluorescent protein fusion construct to the NF-kappaB transcription factor subunit RelA (GFP-RelA). Translocation kinetics in response to live bacteria and purified lipopolysaccharide (LPS) challenges were measured, and this work presents the first demonstration of live imaging of host cell infection on a microfluidic platform with quantitative analysis of an early (<0.5 h from infection) immune signaling event. Our data show that a 1,000x increase in the LPS dose led to a ~10x increase in a host cell activation metric we developed in order to describe NF-kappaB translocation kinetics. Using this metric, live bacteria challenges were assigned an equivalent LPS dose as a first step towards comparing NF-kappaB translocation kinetics between TLR4-only pathway signaling (activated by LPS) and multiple pathway signaling (activated by whole bacteria). The device also contains a unique architecture for capturing and fluidically isolating single host cells for the purpose of differentiating between primary and secondary immune signaling.

Multiple roles of the NF-kappaB signaling pathway regulated by coupled negative feedback circuits

The NF-kappaB signaling pathway can perform multiple functional roles depending on specific cellular environments and cell types. Even in the same cell clones, the pathway can show different kinetic and phenotypic properties. It is believed that the complex networks controlling the NF-kappaB signaling pathway can generate these diverse and sometimes ambiguous phenomena. We noted, however, that the dynamics of NF-kappaB signaling pathway is highly stochastic and that the NF-kappaB signaling pathway contains multiple negative feedback circuits formed by IkappaB isoform proteins, IkappaBalpha and IkappaBepsilon in particular. Considering the topological similarity, their functional roles seem to be redundant, raising the question why different types of IkappaB isoforms need to exist. From extensive stochastic simulations of the NF-kappaB signaling pathway, we found that each IkappaB isoform actually conducts a different regulatory role through its own negative feedback. Specifically, our data suggest that IkappaBalpha controls the dynamic patterns of nuclear NF-kappaB, while IkappaBepsilon induces cellular heterogeneity of the NF-kappaB activities. These results may provide an answer to the question of how a single NF-kappaB signaling pathway can perform multiple biological functions even in the same clonal populations.

Pulsatile stimulation determines timing and specificity of NF-kappaB-dependent transcription

The nuclear factor kappaB (NF-kappaB) transcription factor regulates cellular stress responses and the immune response to infection. NF-kappaB activation results in oscillations in nuclear NF-kappaB abundance. To define the function of these oscillations, we treated cells with repeated short pulses of tumor necrosis factor-alpha at various intervals to mimic pulsatile inflammatory signals. At all pulse intervals that were analyzed, we observed synchronous cycles of NF-kappaB nuclear translocation. Lower frequency stimulations gave repeated full-amplitude translocations, whereas higher frequency pulses gave reduced translocation, indicating a failure to reset. Deterministic and stochastic mathematical models predicted how negative feedback loops regulate both the resetting of the system and cellular heterogeneity. Altering the stimulation intervals gave different patterns of NF-kappaB-dependent gene expression, which supports the idea that oscillation frequency has a functional role.

Floral morphogenesis: stochastic explorations of a gene network epigenetic landscape

In contrast to the classical view of development as a preprogrammed and deterministic process, recent studies have demonstrated that stochastic perturbations of highly non-linear systems may underlie the emergence and stability of biological patterns. Herein, we address the question of whether noise contributes to the generation of the stereotypical temporal pattern in gene expression during flower development. We modeled the regulatory network of organ identity genes in the Arabidopsis thaliana flower as a stochastic system. This network has previously been shown to converge to ten fixed-point attractors, each with gene expression arrays that characterize inflorescence cells and primordial cells of sepals, petals, stamens, and carpels. The network used is binary, and the logical rules that govern its dynamics are grounded in experimental evidence. We introduced different levels of uncertainty in the updating rules of the network. Interestingly, for a level of noise of around 0.5–10%, the system exhibited a sequence of transitions among attractors that mimics the sequence of gene activation configurations observed in real flowers. We also implemented the gene regulatory network as a continuous system using the Glass model of differential equations, that can be considered as a first approximation of kinetic-reaction equations, but which are not necessarily equivalent to the Boolean model. Interestingly, the Glass dynamics recover a temporal sequence of attractors, that is qualitatively similar, although not identical, to that obtained using the Boolean model. Thus, time ordering in the emergence of cell-fate patterns is not an artifact of synchronous updating in the Boolean model. Therefore, our model provides a novel explanation for the emergence and robustness of the ubiquitous temporal pattern of floral organ specification. It also constitutes a new approach to understanding morphogenesis, providing predictions on the population dynamics of cells with different genetic configurations during development.

Detailed qualitative dynamic knowledge representation using a BioNetGen model of TLR-4 signaling and preconditioning

INTRODUCTION

Intracellular signaling/synthetic pathways are being increasingly extensively characterized. However, while these pathways can be displayed in static diagrams, in reality they exist with a degree of dynamic complexity that is responsible for heterogeneous cellular behavior. Multiple parallel pathways exist and interact concurrently, limiting the ability to integrate the various identified mechanisms into a cohesive whole. Computational methods have been suggested as a means of concatenating this knowledge to aid in the understanding of overall system dynamics. Since the eventual goal of biomedical research is the identification and development of therapeutic modalities, computational representation must have sufficient detail to facilitate this 'engineering' process. Adding to the challenge, this type of representation must occur in a perpetual state of incomplete knowledge. We present a modeling approach to address this challenge that is both detailed and qualitative. This approach is termed 'dynamic knowledge representation,' and is intended to be an integrated component of the iterative cycle of scientific discovery.

METHODS

BioNetGen (BNG), a software platform for modeling intracellular signaling pathways, was used to model the toll-like receptor 4 (TLR-4) signal transduction cascade. The informational basis of the model was a series of reference papers on modulation of (TLR-4) signaling, and some specific primary research papers to aid in the characterization of specific mechanistic steps in the pathway. This model was detailed with respect to the components of the pathway represented, but qualitative with respect to the specific reaction coefficients utilized to execute the reactions. Responsiveness to simulated lipopolysaccharide (LPS) administration was measured by tumor necrosis factor (TNF) production. Simulation runs included evaluation of initial dose-dependent response to LPS administration at 10, 100, 1000 and 10,000, and a subsequent examination of preconditioning behavior with increasing LPS at 10, 100, 1000 and 10,000 and a secondary dose of LPS at 10,000 administered at approximately 27h of simulated time. Simulations of 'knockout' versions of the model allowed further examination of the interactions within the signaling cascade.

RESULTS

The model demonstrated a dose-dependent TNF response curve to increasing stimulus by LPS. Preconditioning simulations demonstrated a similar dose-dependency of preconditioning doses leading to attenuation of response to subsequent LPS challenge - a 'tolerance' dynamic. These responses match dynamics reported in the literature. Furthermore, the simulated 'knockout' results suggested the existence and need for dual negative feedback control mechanisms, represented by the zinc ring-finger protein A20 and inhibitor kappa B proteins (IkappaB), in order for both effective attenuation of the initial stimulus signal and subsequent preconditioned 'tolerant' behavior.

CONCLUSIONS

We present an example of detailed, qualitative dynamic knowledge representation using the TLR-4 signaling pathway, its control mechanisms and overall behavior with respect to preconditioning. The intent of this approach is to demonstrate a method of translating the extensive mechanistic knowledge being generated at the basic science level into an executable framework that can provide a means of 'conceptual model verification.' This allows for both the 'checking' of the dynamic consequences of a mechanistic hypothesis and the creation of a modular component of an overall model directed at the engineering goal of biomedical research. It is hoped that this paper will increase the use of knowledge representation and communication in this fashion, and facilitate the concatenation and integration of community-wide knowledge.

Mechanistic simulations of inflammation: current state and future prospects

Inflammation is a normal, robust physiological process. It can also be viewed as a complex system that senses and attempts to resolve homeostatic perturbations initiated from within the body (for example, in autoimmune disease) or from the outside (for example, in infections). Virtually all acute and chronic diseases are either driven or modulated by inflammation. The complex interplay between beneficial and harmful arms of the inflammatory response may underlie the lack of fully effective therapies for many diseases. Mathematical modeling is emerging as a frontline tool for understanding the complexity of the inflammatory response. A series of articles in this issue highlights various modeling approaches to inflammation in the larger context of health and disease, from intracellular signaling to whole-animal physiology. Here we discuss the state of this emerging field. We note several common features of inflammation models, as well as challenges and prospects for future studies.

Virus Infection Induces NF-kappaB-dependent interchromosomal associations mediating monoallelic IFN-beta gene expression

Transcriptional activation of the IFN-beta gene by virus infection requires the cooperative assembly of an enhanceosome. We report that the stochastic and monoallelic expression of the IFN-beta gene depends on interchromosomal associations with three identified distinct genetic loci that could mediate binding of the limiting transcription factor NF-kappaB to the IFN-beta enhancer, thus triggering enhanceosome assembly and activation of transcription from this allele. The probability of a cell to express IFN-beta is dramatically increased when the cell is transfected with any of these loci. The secreted IFN-beta protein induces high-level expression of the enhanceosome factor IRF-7, which in turn promotes enhanceosome assembly and IFN-beta transcription from the remaining alleles and in other initially nonexpressing cells. Thus, the IFN-beta enhancer functions in a nonlinear fashion by working as a signal amplifier.

Modulating temporal control of NF-kappaB activation: implications for therapeutic and assay selection

The activation of transcription factor NF-kappaB (nuclear factor-kappaB) plays a central role in the induction of many inflammatory response genes. This process is characterized by either oscillations or stable induction of NF-kappaB nuclear binding. Changes in dynamics of binding result in the expression of distinct subsets of genes leading to different physiological outcomes. We examined NF-kappaB DNA binding activity in lipopolysaccharide (LPS)-stimulated IC-21 cells by electromobility shift assay and nonradioactive transcription factor assay and interpreted the results using a kinetic model of NF-kappaB activation. Both assays detected damped oscillatory behavior of NF-kappaB with differences in sensitivity and reproducibility. 3,4-Dichloropropionaniline (DCPA) was used to modulate the oscillatory behavior of NF-kappaB after LPS stimulation. DCPA is known to inhibit the production of two NF-kappaB-inducible cytokines, IL-6 and tumor necrosis factor alpha, by reducing but not completely abrogating NF-kappaB-induced transcription. DCPA treatment resulted in a potentiation of early LPS-induced NF-kappaB activation. The nonradioactive transcription factor assay, which has a higher signal/noise ratio than the electromobility shift assay, combined with in silico modeling, produced results that revealed changes in NF-kappaB dynamics which, to the best of our knowledge, have never been previously reported. These results highlight the importance of cell type and stimulus specificity in transcription factor activity assessment. In addition, assay selection has important implications for network inference and drug discovery.

Introduction of an agent-based multi-scale modular architecture for dynamic knowledge representation of acute inflammation

Background

One of the greatest challenges facing biomedical research is the integration and sharing of vast amounts of information, not only for individual researchers, but also for the community at large. Agent Based Modeling (ABM) can provide a means of addressing this challenge via a unifying translational architecture for dynamic knowledge representation. This paper presents a series of linked ABMs representing multiple levels of biological organization. They are intended to translate the knowledge derived from in vitro models of acute inflammation to clinically relevant phenomenon such as multiple organ failure.

Results and Discussion

ABM development followed a sequence starting with relatively direct translation from in-vitro derived rules into a cell-as-agent level ABM, leading on to concatenated ABMs into multi-tissue models, eventually resulting in topologically linked aggregate multi-tissue ABMs modeling organ-organ crosstalk. As an underlying design principle organs were considered to be functionally composed of an epithelial surface, which determined organ integrity, and an endothelial/blood interface, representing the reaction surface for the initiation and propagation of inflammation. The development of the epithelial ABM derived from an in-vitro model of gut epithelial permeability is described. Next, the epithelial ABM was concatenated with the endothelial/inflammatory cell ABM to produce an organ model of the gut. This model was validated against in-vivo models of the inflammatory response of the gut to ischemia. Finally, the gut ABM was linked to a similarly constructed pulmonary ABM to simulate the gut-pulmonary axis in the pathogenesis of multiple organ failure. The behavior of this model was validated against in-vivo and clinical observations on the cross-talk between these two organ systems

Conclusion

A series of ABMs are presented extending from the level of intracellular mechanism to clinically observed behavior in the intensive care setting. The ABMs all utilize cell-level agents that encapsulate specific mechanistic knowledge extracted from in vitro experiments. The execution of the ABMs results in a dynamic representation of the multi-scale conceptual models derived from those experiments. These models represent a qualitative means of integrating basic scientific information on acute inflammation in a multi-scale, modular architecture as a means of conceptual model verification that can potentially be used to concatenate, communicate and advance community-wide knowledge.

Understanding NF-kappaB signaling via mathematical modeling

Mammalian inflammatory signaling, for which NF-kappaB is a principal transcription factor, is an exquisite example of how cellular signaling pathways can be regulated to produce different yet specific responses to different inflammatory insults. Mathematical models, tightly linked to experiment, have been instrumental in unraveling the forms of regulation in NF-kappaB signaling and their underlying molecular mechanisms. Our initial model of the IkappaB-NF-kappaB signaling module highlighted the role of negative feedback in the control of NF-kappaB temporal dynamics and gene expression. Subsequent studies sparked by this work have helped to characterize additional feedback loops, the input-output behavior of the module, crosstalk between multiple NF-kappaB-activating pathways, and NF-kappaB oscillations. We anticipate that computational techniques will enable further progress in the NF-kappaB field, and the signal transduction field in general, and we discuss potential upcoming developments.

Wires in the soup: quantitative models of cell signaling

Living cells are capable of extracting information from their environments and mounting appropriate responses to a variety of challenges. The underlying signal transduction networks enabling this response can be quite complex, and so sophisticated computational modeling coupled with precise experimentation is required to unravel them. Although we are still at the beginning of this process, some recent examples of integrative analysis of cell signaling are very encouraging. Quantitative models of signaling pathways (e.g. NF-kappaB) can be gradually constructed through continuous experimental validation, and important lessons can be learnt from such exercises.