Control of RelB during dendritic cell activation integrates canonical and noncanonical NF-κB pathways.

The NF-κB protein RelB controls dendritic cell (DC) maturation and may be targeted therapeutically to manipulate T cell responses in disease. Here we report that RelB promoted DC activation not as the expected RelB-p52 effector of the noncanonical NF-κB pathway, but as a RelB-p50 dimer regulated by canonical IκBs, IκBα and IκBɛ. IκB control of RelB minimized spontaneous maturation but enabled rapid pathogen-responsive maturation. Computational modeling of the NF-κB signaling module identified control points of this unexpected cell type-specific regulation. Fibroblasts that we engineered accordingly showed DC-like RelB control. Canonical pathway control of RelB regulated pathogen-responsive gene expression programs. This work illustrates the potential utility of systems analyses in guiding the development of combination therapeutics for modulating DC-dependent T cell responses

Estimation of kinetic rates of MAP kinase activation from experimental data

Mathematical model is an important tool in systems biology to study the dynamics of biological systems inside the cell. One of the significant challenges in systems biology is the lack of kinetic rates that should be measured in experiments or estimated from experimental data. This work addresses this issue by using a genetic algorithm to estimate reaction rates related to the phosphorylation and dephosphorylation of MAP kinase (ERK) in the mitogen-activated protein (MAP) kinase pathway from biological measurements. In addition, we discuss the robustness of the mathematical model with regards to the variation of kinetic rates together with external noise due to environmental fluctuations. This has been proposed as an additional criterion to choose the estimate from the candidate parameter sets that are obtained from different implementations of the genetic algorithm

Perturbation analysis analyzed—mathematical modeling of intact and perturbed gene regulatory circuits for animal development

Gene regulatory networks for animal development are the underlying mechanisms controlling cell fate specification and differentiation. The architecture of gene regulatory circuits determines their information processing properties and their developmental function. It is a major task to derive realistic network models from exceedingly advanced high throughput experimental data. Here we use mathematical modeling to study the dynamics of gene regulatory circuits to advance the ability to infer regulatory connections and logic function from experimental data. This study is guided by experimental methodologies that are commonly used to study gene regulatory networks that control cell fate specification. We study the effect of a perturbation of an input on the level of its downstream genes and compare between the cis-regulatory execution of OR and AND logics. Circuits that initiate gene activation and circuits that lock on the expression of genes are analyzed. The model improves our ability to analyze experimental data and construct from it the network topology. The model also illuminates information processing properties of gene regulatory circuits for animal development

Approximations and their consequences for dynamic modelling of signal transduction pathways

Signal transduction is the process by which the cell converts one kind of signal or stimulus into another. This involves a sequence of biochemical reactions, carried out by proteins. The dynamic response of complex cell signalling networks can be modelled and simulated in the framework of chemical kinetics. The mathematical formulation of chemical kinetics results in a system of coupled differential equations. Simplifications can arise through assumptions and approximations. The paper provides a critical discussion of frequently employed approximations in dynamic modelling of signal transduction pathways. We discuss the requirements for conservation laws, steady state approximations, and the neglect of components. We show how these approximations simplify the mathematical treatment of biochemical networks but we also demonstrate differences between the complete system and its approximations with respect to the transient and steady state behavior

Novel decay dynamics revealed for virus-mediated drug activation in cytomegalovirus infection

Human cytomegalovirus (CMV) infection is a substantial cause of morbidity and mortality in immunocompromised hosts and globally is one of the most important congenital infections. The nucleoside analogue ganciclovir (GCV), which requires initial phosphorylation by the viral UL97 kinase, is the mainstay for treatment. To date, CMV decay kinetics during GCV therapy have not been extensively investigated and its clinical implications not fully appreciated. We measured CMV DNA levels in the blood of 92 solid organ transplant recipients with CMV disease over the initial 21 days of ganciclovir therapy and identified four distinct decay patterns, including a new pattern exhibiting a transient viral rebound (Hump) following initial decline. Since current viral dynamics models were unable to account for this Hump profile, we developed a novel multi-level model, which includes the intracellular role of UL97 in the continued activation of ganciclovir, that successfully described all the decline patterns observed. Fitting the data allowed us to estimate ganciclovir effectiveness in vivo (mean 92%), infected cell half-life (mean 0.7 days), and other viral dynamics parameters that determine which of the four kinetic patterns will ensue. An important clinical implication of our results is that the virological efficacy of GCV operates over a broad dose range. The model also raises the possibility that GCV can drive replication to a new lower steady state but ultimately cannot fully eradicate it. This model is likely to be generalizable to other anti-CMV nucleoside analogs that require activation by viral enzymes such as UL97 or its homologues

The case for absolute ligand discrimination : modeling information processing and decision by immune T cells

Some cells have to take decision based on the quality of surroundings ligands, almost irrespective of their quantity, a problem we name "absolute discrimination". An example of absolute discrimination is recognition of not-self by immune T Cells. We show how the problem of absolute discrimination can be solved by a process called "adaptive sorting". We review several implementations of adaptive sorting, as well as its generic properties such as antagonism. We show how kinetic proofreading with negative feedback implements an approximate version of adaptive sorting in the immune context. Finally, we revisit the decision problem at the cell population level, showing how phenotypic variability and feedbacks between population and single cells are crucial for proper decision

Boolean network model predicts cell cycle sequence of fission yeast

A Boolean network model of the cell-cycle regulatory network of fission yeast (Schizosaccharomyces Pombe) is constructed solely on the basis of the known biochemical interaction topology. Simulating the model in the computer, faithfully reproduces the known sequence of regulatory activity patterns along the cell cycle of the living cell. Contrary to existing differential equation models, no parameters enter the model except the structure of the regulatory circuitry. The dynamical properties of the model indicate that the biological dynamical sequence is robustly implemented in the regulatory network, with the biological stationary state G1 corresponding to the dominant attractor in state space, and with the biological regulatory sequence being a strongly attractive trajectory. Comparing the fission yeast cell-cycle model to a similar model of the corresponding network in S. cerevisiae, a remarkable difference in circuitry, as well as dynamics is observed. While the latter operates in a strongly damped mode, driven by external excitation, the S. pombe network represents an auto-excited system with external damping.Comment: 10 pages, 3 figure

Simulating non-small cell lung cancer with a multiscale agent-based model

Background The epidermal growth factor receptor (EGFR) is frequently overexpressed in many cancers, including non-small cell lung cancer (NSCLC). In silcio modeling is considered to be an increasingly promising tool to add useful insights into the dynamics of the EGFR signal transduction pathway. However, most of the previous modeling work focused on the molecular or the cellular level only, neglecting the crucial feedback between these scales as well as the interaction with the heterogeneous biochemical microenvironment. Results We developed a multiscale model for investigating expansion dynamics of NSCLC within a two-dimensional in silico microenvironment. At the molecular level, a specific EGFR-ERK intracellular signal transduction pathway was implemented. Dynamical alterations of these molecules were used to trigger phenotypic changes at the cellular level. Examining the relationship between extrinsic ligand concentrations, intrinsic molecular profiles and microscopic patterns, the results confirmed that increasing the amount of available growth factor leads to a spatially more aggressive cancer system. Moreover, for the cell closest to nutrient abundance, a phase-transition emerges where a minimal increase in extrinsic ligand abolishes the proliferative phenotype altogether. Conclusions Our in silico results indicate that, in NSCLC, in the presence of a strong extrinsic chemotactic stimulus, and depending on the cell's location, downstream EGFR-ERK signaling may be processed more efficiently, thereby yielding a migration-dominant cell phenotype and overall, an accelerated spatio-temporal expansion rate.Comment: 37 pages, 7 figure