The competitive nature of signal transducer and activator of transcription complex formation drives phenotype switching of T cells.

Signal transducers and activators of transcription (STATs) are key molecular determinants of T-cell fate and effector function. Several inflammatory diseases are characterized by an altered balance of T-cell phenotypes and cytokine secretion. STATs, therefore, represent viable therapeutic targets in numerous pathologies. However, the underlying mechanisms by which the same STAT proteins regulate both the development of different T-cell phenotypes and their plasticity during changes in extracellular conditions remain unclear. In this study, we investigated the STAT-mediated regulation of T-cell phenotype formation and plasticity using mathematical modelling and experimental data for intracellular STAT signalling proteins. The close fit of our model predictions to the experimental data allows us to propose a potential mechanism for T-cell switching. According to this mechanism, T-cell phenotype switching is the result of the relative redistribution of STAT dimer complexes caused by the extracellular cytokine-dependent STAT competition effects. The developed model predicts that the balance between the intracellular STAT species defines the amount of the produced cytokines and thereby T-cell phenotypes. The model predictions are consistent with the experimentally observed interferon-γ to interleukin-10 switching that regulates human T helper type 1/type 1 regulatory T-cell responses. The proposed model is applicable to a number of STAT signalling circuits.

A systems model for immune cell interactions unravels the mechanism of inflammation in human skin.

Inflammation is characterized by altered cytokine levels produced by cell populations in a highly interdependent manner. To elucidate the mechanism of an inflammatory reaction, we have developed a mathematical model for immune cell interactions via the specific, dose-dependent cytokine production rates of cell populations. The model describes the criteria required for normal and pathological immune system responses and suggests that alterations in the cytokine production rates can lead to various stable levels which manifest themselves in different disease phenotypes. The model predicts that pairs of interacting immune cell populations can maintain homeostatic and elevated extracellular cytokine concentration levels, enabling them to operate as an immune system switch. The concept described here is developed in the context of psoriasis, an immune-mediated disease, but it can also offer mechanistic insights into other inflammatory pathologies as it explains how interactions between immune cell populations can lead to disease phenotypes.

Kinetic regulation of multi-ligand binding proteins.

BACKGROUND: Second messengers, such as calcium, regulate the activity of multisite binding proteins in a concentration-dependent manner. For example, calcium binding has been shown to induce conformational transitions in the calcium-dependent protein calmodulin, under steady state conditions. However, intracellular concentrations of these second messengers are often subject to rapid change. The mechanisms underlying dynamic ligand-dependent regulation of multisite proteins require further elucidation. RESULTS: In this study, a computational analysis of multisite protein kinetics in response to rapid changes in ligand concentrations is presented. Two major physiological scenarios are investigated: i) Ligand concentration is abundant and the ligand-multisite protein binding does not affect free ligand concentration, ii) Ligand concentration is of the same order of magnitude as the interacting multisite protein concentration and does not change. Therefore, buffering effects significantly influence the amounts of free ligands. For each of these scenarios the influence of the number of binding sites, the temporal effects on intermediate apo- and fully saturated conformations and the multisite regulatory effects on target proteins are investigated. CONCLUSIONS: The developed models allow for a novel and accurate interpretation of concentration and pressure jump-dependent kinetic experiments. The presented model makes predictions for the temporal distribution of multisite protein conformations in complex with variable numbers of ligands. Furthermore, it derives the characteristic time and the dynamics for the kinetic responses elicited by a ligand concentration change as a function of ligand concentration and the number of ligand binding sites. Effector proteins regulated by multisite ligand binding are shown to depend on ligand concentration in a highly nonlinear fashion.

Erythrocyte: A systems model of the control of aggregation and deformability.

Human erythrocytes are highly specialized enucleate cells that are involved in providing efficient gas transport. Erythrocytes have been extensively studied both experimentally and by mathematical modeling in recent years. However, understanding of how aggregation and deformability are regulated is limited. These properties of the erythrocyte are essential for the physiological functioning of the cell. In this work, we propose a novel mathematical model of the molecular system that controls the aggregation and deformability of the erythrocyte. This model is based on the experimental results of previously published studies. Our model suggests fundamentally new mechanisms that regulate aggregation and deformability in a latch-like manner. The results of this work could be used as a general explanation of how the erythrocytes regulate their aggregation and deformability, and are essential in understanding erythrocyte disorders and aging.

A systems model of phosphorylation for inflammatory signaling events.

Phosphorylation is a fundamental biochemical reaction that modulates protein activity in cells. While a single phosphorylation event is relatively easy to understand, multisite phosphorylation requires systems approaches for deeper elucidation of the underlying molecular mechanisms. In this paper we develop a mechanistic model for single- and multi-site phosphorylation. The proposed model is compared with previously reported studies. We compare the predictions of our model with experiments published in the literature in the context of inflammatory signaling events in order to provide a mechanistic description of the multisite phosphorylation-mediated regulation of Signal Transducer and Activator of Transcription 3 (STAT3) and Interferon Regulatory Factor 5 (IRF-5) proteins. The presented model makes crucial predictions for transcription factor phosphorylation events in the immune system. The model proposes potential mechanisms for T cell phenotype switching and production of cytokines. This study also provides a generic framework for the better understanding of a large number of multisite phosphorylation-regulated biochemical circuits.

Computational modelling elucidates the mechanism of ciliary regulation in health and disease.

BACKGROUND: Ciliary dysfunction leads to a number of human pathologies, including primary ciliary dyskinesia, nephronophthisis, situs inversus pathology or infertility. The mechanism of cilia beating regulation is complex and despite extensive experimental characterization remains poorly understood. We develop a detailed systems model for calcium, membrane potential and cyclic nucleotide-dependent ciliary motility regulation. RESULTS: The model describes the intimate relationship between calcium and potassium ionic concentrations inside and outside of cilia with membrane voltage and, for the first time, describes a novel type of ciliary excitability which plays the major role in ciliary movement regulation. Our model describes a mechanism that allows ciliary excitation to be robust over a wide physiological range of extracellular ionic concentrations. The model predicts the existence of several dynamic modes of ciliary regulation, such as the generation of intraciliary Ca2+ spike with amplitude proportional to the degree of membrane depolarization, the ability to maintain stable oscillations, monostable multivibrator regimes, all of which are initiated by variability in ionic concentrations that translate into altered membrane voltage. CONCLUSIONS: Computational investigation of the model offers several new insights into the underlying molecular mechanisms of ciliary pathologies. According to our analysis, the reported dynamic regulatory modes can be a physiological reaction to alterations in the extracellular environment. However, modification of the dynamic modes, as a result of genetic mutations or environmental conditions, can cause a life threatening pathology.

Systems biology of cell behavior.

Systems Biology approaches to drug discovery largely focus on the increasing understanding of intracellular and cellular circuits, by computational representation of a molecular system followed by parameter validation against experimental data. This chapter outlines a universal approach to systems biology that allows the linking of intracellular molecular machinery and cellular activity. This procedure is achieved by applying mathematical modeling to molecular modules of a cell in the light of systems biology techniques.

Crosstalk between G-protein and Ca2+ pathways switches intracellular cAMP levels.

Cyclic adenosine monophosphate and cyclic guanosine monophosphate are universal intracellular messengers whose concentrations are regulated by molecular networks comprised of different isoforms of the synthases adenylate cyclase or guanylate cyclase and the phosphodiesterases which degrade these compounds. In this paper, we employ a systems biology approach to develop mathematical models of these networks that, for the first time, take into account the different biochemical properties of the isoforms involved. To investigate the mechanisms underlying the joint regulation of cAMP and cGMP, we apply our models to analyse the regulation of cilia beat frequency in Paramecium by Ca(2+). Based on our analysis of these models, we propose that the diversity of isoform combinations that occurs in living cells provides an explanation for the huge variety of intracellular processes that are dependent on these networks. The inclusion of both G-protein receptor and Ca(2+)-dependent regulation of AC in our models allows us to propose a new explanation for the switching properties of G-protein subunits involved in nucleotide regulation. Analysis of the models suggests that, depending on whether the G-protein subunit is bound to AC, Ca(2+) can either activate or inhibit AC in a concentration-dependent manner. The resulting analysis provides an explanation for previous experimental results that showed that alterations in Ca(2+) concentrations can either increase or decrease cilia beat frequency over particular Ca(2+) concentration ranges.

Elucidating the mechanisms of cooperative calcium-calmodulin interactions: a structural systems biology approach.

BACKGROUND: Calmodulin is an important multifunctional molecule that regulates the activities of a large number of proteins in the cell. Calcium binding induces conformational transitions in calmodulin that make it specifically active to particular target proteins. The precise mechanisms underlying calcium binding to calmodulin are still, however, quite poorly understood. RESULTS: In this study, we adopt a structural systems biology approach and develop a mathematical model to investigate various types of cooperative calcium-calmodulin interactions. We compare the predictions of our analysis with physiological dose-response curves taken from the literature, in order to provide a quantitative comparison of the effects of different mechanisms of cooperativity on calcium-calmodulin interactions. The results of our analysis reduce the gap between current understanding of intracellular calmodulin function at the structural level and physiological calcium-dependent calmodulin target activation experiments. CONCLUSION: Our model predicts that the specificity and selectivity of CaM target regulation is likely to be due to the following factors: variations in the target-specific Ca2+ dissociation and cooperatively effected dissociation constants, and variations in the number of Ca2+ ions required to bind CaM for target activation.