Prof. Michał Komorowski, PhD, DSc

Department of Biosystems and Soft Matter (ZBiMM)
Division of Modelling in Biology and Medicine (PMBM)
position: Associate Professor
telephone: (+48) 22 826 12 81 ext.: 449
room: 324
e-mail: mkomor
personal site: http://bluebox.ippt.pan.pl/~mkomor/

Doctoral thesis
2009-11-30Statistical methods for estimation of biochemical kinetic parameters  (DS,UW)
supervisor -- Bärbel Finkenstädt, DS,UW
supervisor -- Prof. David A. Rand, DS,UW
1242
 
Habilitation thesis
2018-10-23Rozwój metod analizy stochastycznych modeli reakcji biochemicznych 
Professor
2023-07-24Title of professor
Supervision of doctoral theses
1.2023-12-01Nienałtowski Karol  Parametric and non-parametric methods to address complexity of cellular signaling pathways1468
 
2.2023-06-29Topolewski Piotr  Sources of the cell-to-cell heterogeneity in interferon gamma and oncostatin M signaling responses1487
 
3.2020-10-06Jetka Tomasz Quantitative methods to describe information flow in biochemical signalling pathways 

Recent publications
1.Bartak M., Krahel W., Chodkowski M., Grel H., Walczak J., Pallepati A., Komorowski M., Cymerys J., ATPase Valosin-Containing Protein (VCP) Is Involved During the Replication and Egress of Sialodacryoadenitis Virus (SDAV) in Neurons, International Journal of Molecular Sciences, ISSN: 1422-0067, DOI: 10.3390/ijms252111633, Vol.25, No.21, pp.11633-1-23, 2024
Abstract:

Sialodacryoadenitis virus (SDAV) has been identified as the etiological agent responsible for the respiratory system and salivary gland infections in rats. The existing literature on SDAV infections is insufficient to address the topic adequately, particularly in relation to the central nervous system. In order to ascertain how SDAV gains access to neuronal cells and subsequently exits, our attention was focused on the small molecule valosin-containing protein (VCP), which is an ATPase. VCP is acknowledged for its function in the ubiquitin-mediated proteasomal degradation of proteins, including those of viral origin. To ascertain the potential influence of VCP on SDAV replication and egress, high-content screening was employed to determine the viral titer and protein content. Western blot analysis was employed to ascertain the relative expression of VCP. Real-time imaging of SDAV-infected cells and confocal imaging for qualitative morphological analysis were conducted. The Eeyarestatin I (EerI) inhibitor was employed to disrupt VCP involvement in the endoplasmic reticulum-associated protein degradation pathway (ERAD) in both pre- and post-incubation systems, with concentrations of 5 μM/mL and 25 μM/mL, respectively. We demonstrated for the first time that SDAV productively replicates in cultured primary neurons. VCP expression is markedly elevated during SDAV infection. The application of 5 μM/mL EerI in the post-treatment system yielded a statistically significant inhibition of the SDAV yield. It is likely that this modulates the efficacy of virion assembly by arresting viral proteins in the submembrane area.

Keywords:

SDAV,VCP,primary neurons,virion assembly,ERAD,eeyarestatin I

Affiliations:
Bartak M.-other affiliation
Krahel W.-other affiliation
Chodkowski M.-other affiliation
Grel H.-other affiliation
Walczak J.-IPPT PAN
Pallepati A.-IPPT PAN
Komorowski M.-IPPT PAN
Cymerys J.-other affiliation
2.Braniewska A., Skorzynski M., Sas Z., Dlugolecka M., Marszalek I., Kurpiel D., Marcel B., Strzemecki D., Magiera A., Bialasek M., Walczak J., Cheda Ł., Komorowski M., Tobias W., Czystowska-Kuzmicz M., Kwapiszewska K., Alberto B., Krol M., Rygiel Tomasz P., A novel process for transcellular hemoglobin transport from macrophages to cancer cells, Cell Communication and Signaling, ISSN: 1478-811X, DOI: 10.1186/s12964-024-01929-8, Vol.22, pp.570-1-21, 2024
Abstract:

Hemoglobin (Hb) performs its physiological function within the erythrocyte. Extracellular Hb has prooxidative and proinflammatory properties and is therefore sequestered by haptoglobin and bound by the CD163 receptor on macrophages. In the present study, we demonstrate a novel process of Hb uptake by macrophages independent of haptoglobin and CD163. Unexpectedly, macrophages do not degrade the entire Hb, but instead transfer it to neighboring cells. We have shown that the phenomenon of Hb transfer from macrophages to other cells is mainly mediated by extracellular vesicles. In contrast to the canonical Hb degradation pathway by macrophages, Hb transfer has not been reported before. In addition, we have used the process of Hb transfer in anticancer therapy, where macrophages are loaded with a Hb-anticancer drug conjugate and act as cellular drug carriers. Both mouse and human macrophages loaded with Hb-monomethyl auristatin E (MMAE) effectively killed cancer cells when co-cultured in vitro.

Keywords:

Hemoglobin,Macrophages,CD163,Extracellular vesicles,Monomethyl auristatin E

Affiliations:
Braniewska A.-other affiliation
Skorzynski M.-other affiliation
Sas Z.-other affiliation
Dlugolecka M.-other affiliation
Marszalek I.-other affiliation
Kurpiel D.-other affiliation
Marcel B.-other affiliation
Strzemecki D.-other affiliation
Magiera A.-other affiliation
Bialasek M.-other affiliation
Walczak J.-IPPT PAN
Cheda Ł.-other affiliation
Komorowski M.-IPPT PAN
Tobias W.-other affiliation
Czystowska-Kuzmicz M.-other affiliation
Kwapiszewska K.-other affiliation
Alberto B.-other affiliation
Krol M.-other affiliation
Rygiel Tomasz P.-other affiliation
3.Komorowski M., Making sense of BMP signaling complexity, Cell Systems, ISSN: 2405-4712, DOI: 10.1016/j.cels.2022.04.002, Vol.13, No.5, pp.349-351, 2022
Abstract:

Cellular signaling systems are immensely complex. Dedicated experimental and theoretical approaches are therefore required to decipher how they function. In this issue of Cell Systems, two studies systematically interrogate the Bone Morphogenetic Protein (BMP) pathway, uncovering mechanisms and consequences of distinct responses to combinations of BMP ligands.

Affiliations:
Komorowski M.-IPPT PAN
4.Topolewski P., Zakrzewska K.E., Walczak J., Nienałtowski K., Müller-Newen G., Singh A., Komorowski M., Phenotypic variability, not noise, accounts for most of the cell-to-cell heterogeneity in IFN-γ and oncostatin M signaling responses, Science Signaling, ISSN: 1945-0877, DOI: 10.1126/scisignal.abd9303, Vol.15, No.721, pp.eabd9303-1-16, 2022
Abstract:

Cellular signaling responses show substantial cell-to-cell heterogeneity, which is often ascribed to the inherent randomness of biochemical reactions, termed molecular noise, wherein high noise implies low signaling fidelity. Alternatively, heterogeneity could arise from differences in molecular content between cells, termed molecular phenotypic variability, which does not necessarily imply imprecise signaling. The contribution of these two processes to signaling heterogeneity is unclear. Here, we fused fibroblasts to produce binuclear syncytia to distinguish noise from phenotypic variability in the analysis of cytokine signaling. We reasoned that the responses of the two nuclei within one syncytium could approximate the signaling outcomes of two cells with the same molecular content, thereby disclosing noise contribution, whereas comparison of different syncytia should reveal contribution of phenotypic variability. We found that ~90% of the variance in the primary response (which was the abundance of phosphorylated, nuclear STAT) to stimulation with the cytokines interferon-γ and oncostatin M resulted from differences in the molecular content of individual cells. Thus, our data reveal that cytokine signaling in the system used here operates in a reproducible, high-fidelity manner.

Affiliations:
Topolewski P.-IPPT PAN
Zakrzewska K.E.-IPPT PAN
Walczak J.-IPPT PAN
Nienałtowski K.-IPPT PAN
Müller-Newen G.-RWTH Aachen University (DE)
Singh A.-University of Delaware (US)
Komorowski M.-IPPT PAN
5.Nienałtowski K., Rigby R.E., Walczak J., Zakrzewska K.E., Głów E., Rehwinkel J., Komorowski M., Fractional response analysis reveals logarithmic cytokine responses in cellular populations, Nature Communications, ISSN: 2041-1723, DOI: 10.1038/s41467-021-24449-2, Vol.12, pp.4175-1-10, 2021
Abstract:

Although we can now measure single-cell signaling responses with multivariate, high-throughput techniques our ability to interpret such measurements is still limited. Even interpretation of dose–response based on single-cell data is not straightforward: signaling responses can differ significantly between cells, encompass multiple signaling effectors, and have dynamic character. Here, we use probabilistic modeling and information-theory to introduce fractional response analysis (FRA), which quantifies changes in fractions of cells with given response levels. FRA can be universally performed for heterogeneous, multivariate, and dynamic measurements and, as we demonstrate, quantifies otherwise hidden patterns in single-cell data. In particular, we show that fractional responses to type I interferon in human peripheral blood mononuclear cells are very similar across different cell types, despite significant differences in mean or median responses and degrees of cell-to-cell heterogeneity. Further, we demonstrate that fractional responses to cytokines scale linearly with the log of the cytokine dose, which uncovers that heterogeneous cellular populations are sensitive to fold-changes in the dose, as opposed to additive changes.

Affiliations:
Nienałtowski K.-other affiliation
Rigby R.E.-University of Oxford (GB)
Walczak J.-IPPT PAN
Zakrzewska K.E.-IPPT PAN
Głów E.-IPPT PAN
Rehwinkel J.-University of Oxford (GB)
Komorowski M.-IPPT PAN
6.Fiebelkow J., Guendel A., Guendel B., Mehwald N., Jetka T., Komorowski M., Waldherr S., Schaper F., Dittrich A., The tyrosine phosphatase SHP2 increases robustness and information transfer within IL-6-induced JAK/STAT signalling, Cell Communication and Signaling, ISSN: 1478-811X, DOI: 10.1186/s12964-021-00770-7, Vol.19, No.1, pp.94-1-19, 2021
Abstract:

Background - Cell-to-cell heterogeneity is an inherent feature of multicellular organisms and is central in all physiological and pathophysiological processes including cellular signal transduction. The cytokine IL-6 is an essential mediator of pro- and anti-inflammatory processes. Dysregulated IL-6-induced intracellular JAK/STAT signalling is associated with severe inflammatory and proliferative diseases. Under physiological conditions JAK/STAT signalling is rigorously controlled and timely orchestrated by regulatory mechanisms such as expression of the feedback-inhibitor SOCS3 and activation of the protein-tyrosine phosphatase SHP2 (PTPN11). Interestingly, the function of negative regulators seems not to be restricted to controlling the strength and timely orchestration of IL-6-induced STAT3 activation. Exemplarily, SOCS3 increases robustness of late IL-6-induced STAT3 activation against heterogenous STAT3 expression and reduces the amount of information transferred through JAK/STAT signalling. Methods - Here we use multiplexed single-cell analyses and information theoretic approaches to clarify whether also SHP2 contributes to robustness of STAT3 activation and whether SHP2 affects the amount of information transferred through IL-6-induced JAK/STAT signalling. Results - SHP2 increases robustness of both basal, cytokine-independent STAT3 activation and early IL-6-induced STAT3 activation against differential STAT3 expression. However, SHP2 does not affect robustness of late IL-6-induced STAT3 activation. In contrast to SOCS3, SHP2 increases the amount of information transferred through IL-6-induced JAK/STAT signalling, probably by reducing cytokine-independent STAT3 activation and thereby increasing sensitivity of the cells. These effects are independent of SHP2-dependent MAPK activation. Conclusion - In summary, the results of this study extend our knowledge of the functions of SHP2 in IL-6-induced JAK/STAT signalling. SHP2 is not only a repressor of basal and cytokine-induced STAT3 activity, but also ensures robustness and transmission of information. Plain English summary - Cells within a multicellular organism communicate with each other to exchange information about the environment. Communication between cells is facilitated by soluble molecules that transmit information from one cell to the other. Cytokines such as interleukin-6 are important soluble mediators that are secreted when an organism is faced with infections or inflammation. Secreted cytokines bind to receptors within the membrane of their target cells. This binding induces activation of an intracellular cascade of reactions called signal transduction, which leads to cellular responses. An important example of intracellular signal transduction is JAK/STAT signalling. In healthy organisms signalling is controlled and timed by regulatory mechanisms, whose activation results in a controlled shutdown of signalling pathways. Interestingly, not all cells within an organism are identical. They differ in the amount of proteins involved in signal transduction, such as STAT3. These differences shape cellular communication and responses to intracellular signalling. Here, we show that an important negative regulatory protein called SHP2 (or PTPN11) is not only responsible for shutting down signalling, but also for steering signalling in heterogeneous cell populations. SHP2 increases robustness of STAT3 activation against variable STAT3 amounts in individual cells. Additionally, it increases the amount of information transferred through JAK/STAT signalling by increasing the dynamic range of pathway activation in heterogeneous cell populations. This is an amazing new function of negative regulatory proteins that contributes to communication in heterogeneous multicellular organisms in health and disease.

Keywords:

signal transduction, SHP2, PTPN11, JAK/STAT, MAPK, information theory, channel capacity, mutual information

Affiliations:
Fiebelkow J.-other affiliation
Guendel A.-other affiliation
Guendel B.-other affiliation
Mehwald N.-other affiliation
Jetka T.-other affiliation
Komorowski M.-IPPT PAN
Waldherr S.-Katholieke Universiteit Leuven (BE)
Schaper F.-Otto-von-Guericke University (DE)
Dittrich A.-Otto-von-Guericke University (DE)
7.Topolewski P., Komorowski M., Information-theoretic analyses of cellular strategies for achieving high signaling capacity—dynamics, cross-wiring, and heterogeneity of cellular states, Current Opinion in Systems Biology, ISSN: 2452-3100, DOI: 10.1016/j.coisb.2021.06.003, Vol.27, pp.100352-1-9, 2021
Abstract:

An individual eukaryotic cell senses identity and quantity of ligands through molecular receptors and signaling pathways, dynamically activating signaling effectors. A distinct ligand often activates multiple different effectors, and a distinct effector is activated by numerous different ligands, which results in cross-wired signaling. In apparently identical cells, the activity of signaling effectors can vary considerably, raising questions about the accuracy of cellular signaling and the interpretation of heterogeneous responses, as either functional or simply noise. Cell-to-cell variability of signaling outcomes, signaling dynamics, and cross-wiring all give rise to signaling complexity, complicating the analysis of signaling mechanisms. Here, we consider a simple input–output modeling approach of information theory that is suitable to analyze signaling complexity and highlight recent studies that have advanced our understanding of the role different components of signaling complexity play in achieving effective information transfer along cellular signaling pathways.

Keywords:

signaling pathways, hormones, growth factors or cytokines, signaling dynamics, cross-wired signaling, Shannon information, Fisher information

Affiliations:
Topolewski P.-IPPT PAN
Komorowski M.-IPPT PAN
8.Komorowski M., Tawfik D.S., The limited information capacity of cross-reactive sensors drives the evolutionary expansion of signaling, Cell Systems, ISSN: 2405-4712, DOI: 10.1016/j.cels.2018.12.006, Vol.8, No.1, pp.76-85.e6, 2019
Abstract:

Signaling systems expand by duplications of various components, be it receptors or downstream effectors. However, whether and how duplicated components contribute to higher signaling capacity is unclear, especially because in most cases, their specificities overlap. Using information theory, we found that augmentation of capacity by an increase in the copy number is strongly limited by logarithmic diminishing returns. Moreover, counter to conventional biochemical wisdom, refinements of the response mechanism, e.g., by cooperativity or allostery, do not increase the overall signaling capacity. However, signaling capacity nearly doubles when a promiscuous, non-cognate ligand becomes explicitly recognized via duplication and partial divergence of signaling components. Our findings suggest that expansion of signaling components via duplication and enlistment of promiscuously acting cues is virtually the only accessible evolutionary strategy to achieve overall high-signaling capacity despite overlapping specificities and molecular noise. This mode of expansion also explains the highly cross-wired architecture of signaling pathways.

Keywords:

paralog expansion, gene duplication, allostery, cooperativity, biochemical signal processing, information capacity

Affiliations:
Komorowski M.-IPPT PAN
Tawfik D.S.-Weizmann Institute of Science (IL)
9.Jetka T., Nienałtowski K., Winarski T., Błoński S., Komorowski M., Information-theoretic analysis of multivariate single-cell signaling responses, PLOS COMPUTATIONAL BIOLOGY, ISSN: 1553-7358, DOI: 10.1371/journal.pcbi.1007132, Vol.15, No.7, pp.e1007132-1-23, 2019
Abstract:

Mathematical methods of information theory appear to provide a useful language to describe how stimuli are encoded in activities of signaling effectors. Exploring the information-theoretic perspective, however, remains conceptually, experimentally and computationally challenging. Specifically, existing computational tools enable efficient analysis of relatively simple systems, usually with one input and output only. Moreover, their robust and readily applicable implementations are missing. Here, we propose a novel algorithm, SLEMI—statistical learning based estimation of mutual information, to analyze signaling systems with high-dimensional outputs and a large number of input values. Our approach is efficient in terms of computational time as well as sample size needed for accurate estimation. Analysis of the NF-κB single—cell signaling responses to TNF-α reveals that NF-κB signaling dynamics improves discrimination of high concentrations of TNF-α with a relatively modest impact on discrimination of low concentrations. Provided R-package allows the approach to be used by computational biologists with only elementary knowledge of information theory.

Affiliations:
Jetka T.-other affiliation
Nienałtowski K.-other affiliation
Winarski T.-IPPT PAN
Błoński S.-IPPT PAN
Komorowski M.-IPPT PAN
10.Billing U., Jetka T., Nortmann L., Wundrack N., Komorowski M., Waldherr S., Schaper F., Dittrich A., Robustness and information transfer within IL-6-induced JAK/STAT signalling, Communications Biology, ISSN: 2399-3642, DOI: 10.1038/s42003-018-0259-4, Vol.2, pp.27-1-14, 2019
Abstract:

Cellular communication via intracellular signalling pathways is crucial. Expression and activation of signalling proteins is heterogenous between isogenic cells of the same cell-type. However, mechanisms evolved to enable sufficient communication and to ensure cellular functions. We use information theory to clarify mechanisms facilitating IL-6-induced JAK/STAT signalling despite cell-to-cell variability. We show that different mechanisms enabling robustness against variability complement each other. Early STAT3 activation is robust as long as cytokine concentrations are low. Robustness at high cytokine concentrations is ensured by high STAT3 expression or serine phosphorylation. Later the feedback-inhibitor SOCS3 increases robustness. Channel Capacity of JAK/STAT signalling is limited by cell-to-cell variability in STAT3 expression and is affected by the same mechanisms governing robustness. Increasing STAT3 amount increases Channel Capacity and robustness, whereas increasing STAT3 tyrosine phosphorylation reduces robustness but increases Channel Capacity. In summary, we elucidate mechanisms preventing dysregulated signalling by enabling reliable JAK/STAT signalling despite cell-to-cell heterogeneity.

Affiliations:
Billing U.-Otto-von-Guericke University (DE)
Jetka T.-other affiliation
Nortmann L.-Otto-von-Guericke University (DE)
Wundrack N.-Otto-von-Guericke University (DE)
Komorowski M.-IPPT PAN
Waldherr S.-Katholieke Universiteit Leuven (BE)
Schaper F.-Otto-von-Guericke University (DE)
Dittrich A.-Otto-von-Guericke University (DE)
11.Jetka T., Nienałtowski K., Filippi S., Stumpf M.P.H., Komorowski M., An information-theoretic framework for deciphering pleiotropic and noisy biochemical signaling, Nature Communications, ISSN: 2041-1723, DOI: 10.1038/s41467-018-07085-1, Vol.9, pp.4591-1-9, 2018
Abstract:

Many components of signaling pathways are functionally pleiotropic, and signaling responses are marked with substantial cell-to-cell heterogeneity. Therefore, biochemical descriptions of signaling require quantitative support to explain how complex stimuli (inputs) are encoded in distinct activities of pathways effectors (outputs). A unique perspective of information theory cannot be fully utilized due to lack of modeling tools that account for the complexity of biochemical signaling, specifically for multiple inputs and outputs. Here, we develop a modeling framework of information theory that allows for efficient analysis of models with multiple inputs and outputs; accounts for temporal dynamics of signaling; enables analysis of how signals flow through shared network components; and is not restricted by limited variability of responses. The framework allows us to explain how identity and quantity of type I and type III interferon variants could be recognized by cells despite activating the same signaling effectors.

Affiliations:
Jetka T.-other affiliation
Nienałtowski K.-other affiliation
Filippi S.-Imperial College London (GB)
Stumpf M.P.H.-Imperial College London (GB)
Komorowski M.-IPPT PAN
12.Martincuks A., Andryka K., Küster A., Schmitz-Van de Leur H., Komorowski M., Müller-Newen G., Nuclear translocation of STAT3 and NF-κB are independent of each other but NF-κB supports expression and activation of STAT3, Cellular Signalling, ISSN: 0898-6568, DOI: 10.1016/j.cellsig.2017.01.006, Vol.32, pp.36-47, 2017
Abstract:

NF-κB and STAT3 are essential transcription factors in immunity and act at the interface of the transition from chronic inflammation to cancer. Different functional crosstalks between NF-κB and STAT3 have been recently described arguing for a direct interaction of both proteins. During a systematic analysis of NF-κB/STAT3 crosstalk we observed that appearance of the subcellular distribution of NF-κB and STAT3 in immunofluorescence heavily depends on the fixation procedure. Therefore, we established an optimized fixation protocol for the reliable simultaneous analysis of the subcellular distributions of both transcription factors. Using this protocol we found that cytokine-induced nuclear accumulation of NF-κB or STAT3 did not alter the subcellular distribution of the other transcription factor. Both knockout and overexpression of STAT3 does not have any major effect on canonical TNFα-NF-κB signalling in MEF or HeLa cells. Similarly, knockout of p65 did not alter nuclear accumulation of STAT3 in response to IL-6. However, p65 expression correlates with elevated total cellular levels of STAT3 and STAT1 and supports activation of these transcription factors. Our findings in MEF cells argue against a direct physical interaction of free cellular NF-κB and STAT3 but point to more intricate functional interactions.

Keywords:

STAT3, NF-κB, Signal transduction, Nuclear translocation, Crosstalk

Affiliations:
Martincuks A.-RWTH Aachen University (DE)
Andryka K.-other affiliation
Küster A.-RWTH Aachen University (DE)
Schmitz-Van de Leur H.-RWTH Aachen University (DE)
Komorowski M.-IPPT PAN
Müller-Newen G.-RWTH Aachen University (DE)
13.Nienałtowski K., Włodarczyk M., Lipniacki T., Komorowski M., Clustering reveals limits of parameter identifiability in multi-parameter models of biochemical dynamics, BMC SYSTEMS BIOLOGY, ISSN: 1752-0509, DOI: 10.1186/s12918-015-0205-8, Vol.9, pp.65-1-9, 2015
Abstract:

Background
Compared to engineering or physics problems, dynamical models in quantitative biology typically depend on a relatively large number of parameters. Progress in developing mathematics to manipulate such multi-parameter models and so enable their efficient interplay with experiments has been slow. Existing solutions are significantly limited by model size.

Results
In order to simplify analysis of multi-parameter models a method for clustering of model parameters is proposed. It is based on a derived statistically meaningful measure of similarity between groups of parameters. The measure quantifies to what extend changes in values of some parameters can be compensated by changes in values of other parameters. The proposed methodology provides a natural mathematical language to precisely communicate and visualise effects resulting from compensatory changes in values of parameters. As a results, a relevant insight into identifiability analysis and experimental planning can be obtained. Analysis of NF- κB and MAPK pathway models shows that highly compensative parameters constitute clusters consistent with the network topology. The method applied to examine an exceptionally rich set of published experiments on the NF- κB dynamics reveals that the experiments jointly ensure identifiability of only 60 % of model parameters. The method indicates which further experiments should be performed in order to increase the number of identifiable parameters.

Conclusions
We currently lack methods that simplify broadly understood analysis of multi-parameter models. The introduced tools depict mutually compensative effects between parameters to provide insight regarding role of individual parameters, identifiability and experimental design. The method can also find applications in related methodological areas of model simplification and parameters estimation.

Affiliations:
Nienałtowski K.-IPPT PAN
Włodarczyk M.-other affiliation
Lipniacki T.-IPPT PAN
Komorowski M.-IPPT PAN
14.Jetka T., Charzyńska A., Gambin A., Stumpf M.P.H., Komorowski M., StochDecomp—Matlab package for noise decomposition in stochastic biochemical systems, BIOINFORMATICS, ISSN: 1367-4803, DOI: 10.1093/bioinformatics/btt631, Vol.30, No.1, pp.137-138, 2014
Abstract:

Motivation: Stochasticity is an indispensable aspect of biochemical processes at the cellular level. Studies on how the noise enters and propagates in biochemical systems provided us with non-trivial insights into the origins of stochasticity, in total, however, they constitute a patchwork of different theoretical analyses.

Results: Here we present a flexible and widely applicable noise decomposition tool that allows us to calculate contributions of individual reactions to the total variability of a system’s output. With the package it is, therefore, possible to quantify how the noise enters and propagates in biochemical systems. We also demonstrate and exemplify using the JAK-STAT signalling pathway that the noise contributions resulting from individual reactions can be inferred from data experimental data along with Bayesian parameter inference. The method is based on the linear noise approximation, which is assumed to provide a reasonable representation of analyzed systems.

Affiliations:
Jetka T.-other affiliation
Charzyńska A.-University of Warsaw (PL)
Gambin A.-other affiliation
Stumpf M.P.H.-Imperial College London (GB)
Komorowski M.-IPPT PAN
15.Woodcock D.J., Vance K.W., Komorowski M., Koentges G., Finkenstädt B., Rand D.A., A hierarchical model of transcriptional dynamics allows robust estimation of transcription rates in populations of single cells with variable gene copy number, BIOINFORMATICS, ISSN: 1367-4803, DOI: 10.1093/bioinformatics/btt201, Vol.29, pp.1519-1525, 2013
Abstract:

Motivation: cis-regulatory DNA sequence elements, such as enhancers and silencers, function to control the spatial and temporal expression of their target genes. Although the overall levels of gene expression in large cell populations seem to be precisely controlled, transcription of individual genes in single cells is extremely variable in real time. It is, therefore, important to understand how these cis-regulatory elements function to dynamically control transcription at single-cell resolution. Recently, statistical methods have been proposed to back calculate the rates involved in mRNA transcription using parameter estimation of a mathematical model of transcription and translation. However, a major complication in these approaches is that some of the parameters, particularly those corresponding to the gene copy number and transcription rate, cannot be distinguished; therefore, these methods cannot be used when the copy number is unknown.

Results: Here, we develop a hierarchical Bayesian model to estimate biokinetic parameters from live cell enhancer–promoter reporter measurements performed on a population of single cells. This allows us to investigate transcriptional dynamics when the copy number is variable across the population. We validate our method using synthetic data and then apply it to quantify the function of two known developmental enhancers in real time and in single cells.

Affiliations:
Woodcock D.J.-University of Warwick (GB)
Vance K.W.-University of Warwick (GB)
Komorowski M.-IPPT PAN
Koentges G.-University of Warwick (GB)
Finkenstädt B.-University of Warwick (GB)
Rand D.A.-University of Warwick (GB)
16.Liepe J., Filippi S., Komorowski M., Stumpf M.P.H., Maximising the information content of experiments in systems biology, PLOS COMPUTATIONAL BIOLOGY, ISSN: 1553-7358, DOI: 10.1371/journal.pcbi.1002888, Vol.9, No.1, pp.e1002888-1-13, 2013
Abstract:

Our understanding of most biological systems is in its infancy. Learning their structure and intricacies is fraught with challenges, and often side-stepped in favour of studying the function of different gene products in isolation from their physiological context. Constructing and inferring global mathematical models from experimental data is, however, central to systems biology. Different experimental setups provide different insights into such systems. Here we show how we can combine concepts from Bayesian inference and information theory in order to identify experiments that maximize the information content of the resulting data. This approach allows us to incorporate preliminary information; it is global and not constrained to some local neighbourhood in parameter space and it readily yields information on parameter robustness and confidence. Here we develop the theoretical framework and apply it to a range of exemplary problems that highlight how we can improve experimental investigations into the structure and dynamics of biological systems and their behavior.

Affiliations:
Liepe J.-Imperial College London (GB)
Filippi S.-Imperial College London (GB)
Komorowski M.-IPPT PAN
Stumpf M.P.H.-Imperial College London (GB)
17.Finkenstädt B., Woodcock D.J., Komorowski M., Harper C.V., Davis J.R.E., White M.R.H., Rand D.A., Quantifying intrinsic and extrinsic noise in gene transcription using the linear noise approximation: An application to single cell data, Annals of Applied Statistics, ISSN: 1932-6157, DOI: 10.1214/13-AOAS669, Vol.7, No.4, pp.1960-1982, 2013
Abstract:

A central challenge in computational modeling of dynamic biological systems is parameter inference from experimental time course measurements. However, one would not only like to infer kinetic parameters but also study their variability from cell to cell. Here we focus on the case where single-cell fluorescent protein imaging time series data are available for a population of cells. Based on van Kampen’s linear noise approximation, we derive a dynamic state space model for molecular populations which is then extended to a hierarchical model. This model has potential to address the sources of variability relevant to single-cell data, namely, intrinsic noise due to the stochastic nature of the birth and death processes involved in reactions and extrinsic noise arising from the cell-to-cell variation of kinetic parameters. In order to infer such a model from experimental data, one must also quantify the measurement process where one has to allow for nonmeasurable molecular species as well as measurement noise of unknown level and variance. The availability of multiple single-cell time series data here provides a unique testbed to fit such a model and quantify these different sources of variation from experimental data.

Keywords:

Linear noise approximation, kinetic parameter estimation, intrinsic and extrinsic noise, state space model and Kalman filter, Bayesian hierarchical modeling

Affiliations:
Finkenstädt B.-University of Warwick (GB)
Woodcock D.J.-University of Warwick (GB)
Komorowski M.-IPPT PAN
Harper C.V.-University of Manchester (GB)
Davis J.R.E.-University of Manchester (GB)
White M.R.H.-University of Manchester (GB)
Rand D.A.-University of Warwick (GB)
18.Schumacher J., Behrends V., Pan Z., Brown D.R., Heydenreich F., Lewis M.R., Bennett M.H., Razzaghi B., Komorowski M., Barahona M., Stumpf M.P.H., Wigneshweraraj S., Bundy J.G., Buck M., Nitrogen and carbon status are integrated at the transcriptional level by the nitrogen regulator NtrC in vivo, mBio, ISSN: 2150-7511, DOI: 10.1128/mBio.00881-13, Vol.4, pp.00881-1-13, 2013
Abstract:

Nitrogen regulation in Escherichia coli is a model system for gene regulation in bacteria. Growth on glutamine as a sole nitrogen source is assumed to be nitrogen limiting, inferred from slow growth and strong NtrB/NtrC-dependent gene activation. However, we show that under these conditions, the intracellular glutamine concentration is not limiting but 5.6-fold higher than in ammonium-replete conditions; in addition, α-ketoglutarate concentrations are elevated. We address this glutamine paradox from a systems perspective. We show that the dominant role of NtrC is to regulate glnA transcription and its own expression, indicating that the glutamine paradox is not due to NtrC-independent gene regulation. The absolute intracellular NtrC and GS concentrations reveal molecular control parameters, where NtrC-specific activities were highest in nitrogen-starved cells, while under glutamine growth, NtrC showed intermediate specific activity. We propose an in vivo model in which α-ketoglutarate can derepress nitrogen regulation despite nitrogen sufficiency.

Affiliations:
Schumacher J.-Imperial College London (GB)
Behrends V.-Imperial College London (GB)
Pan Z.-Imperial College London (GB)
Brown D.R.-Imperial College London (GB)
Heydenreich F.-Imperial College London (GB)
Lewis M.R.-Imperial College London (GB)
Bennett M.H.-Imperial College London (GB)
Razzaghi B.-Imperial College London (GB)
Komorowski M.-IPPT PAN
Barahona M.-Imperial College London (GB)
Stumpf M.P.H.-Imperial College London (GB)
Wigneshweraraj S.-Imperial College London (GB)
Bundy J.G.-Imperial College London (GB)
Buck M.-Imperial College London (GB)
19.Komorowski M., Miękisz J., Stumpf M.P.H., Decomposing Noise in Biochemical Signalling Systems Highlights the Role of Protein Degradation, BIOPHYSICAL JOURNAL, ISSN: 0006-3495, DOI: 10.1016/j.bpj.2013.02.027, Vol.104, pp.1783-1793, 2013
Abstract:

Stochasticity is an essential aspect of biochemical processes at the cellular level. We now know that living cells take advantage of stochasticity in some cases and counteract stochastic effects in others. Here we propose a method that allows us to calculate contributions of individual reactions to the total variability of a system’s output. We demonstrate that reactions differ significantly in their relative impact on the total noise and we illustrate the importance of protein degradation on the overall variability for a range of molecular processes and signaling systems. With our flexible and generally applicable noise decomposition method, we are able to shed new, to our knowledge, light on the sources and propagation of noise in biochemical reaction networks; in particular, we are able to show how regulated protein degradation can be employed to reduce the noise in biochemical systems.

Affiliations:
Komorowski M.-IPPT PAN
Miękisz J.-University of Warsaw (PL)
Stumpf M.P.H.-Imperial College London (GB)
20.Komorowski M., Zurauskiene J., Stumpf M.P.H., StochSens - matlab package for sensitivity analysis of stochastic chemical systems, BIOINFORMATICS, ISSN: 1367-4803, DOI: 10.1093/bioinformatics/btr714, Vol.28, No.5, pp.731-733, 2012
Abstract:

Motivation: The growing interest in the role of stochasticity in biochemical systems drives the demand for tools to analyse stochastic dynamical models of chemical reactions. One powerful tool to elucidate performance of dynamical systems is sensitivity analysis. Traditionally, however, the concept of sensitivity has mainly been applied to deterministic systems, and the difficulty to generalize these concepts for stochastic systems results from necessity of extensive Monte Carlo simulations.

Results: Here we present a Matlab package, StochSens, that implements sensitivity analysis for stochastic chemical systems using the concept of the Fisher Information Matrix (FIM). It uses the linear noise approximation to represent the FIM in terms of solutions of ordinary differential equations. This is the first computational tool that allows for quick computation of the Information Matrix for stochastic systems without the need for Monte Carlo simulations.

Affiliations:
Komorowski M.-other affiliation
Zurauskiene J.-other affiliation
Stumpf M.P.H.-Imperial College London (GB)
21.Harrington H.A., Komorowski M., Beguerisse Díaz M., Ratto G.M., Stumpf M.P.H., Mathematical modeling reveals the functional implications of the different nuclear shuttling rates of Erk1 and Erk2, PHYSICAL BIOLOGY, ISSN: 1478-3967, DOI: 10.1088/1478-3975/9/3/036001, Vol.9, pp.036001-1-12, 2012
Abstract:

The mitogen-activated protein kinase (MAPK) family of proteins is involved in regulating cellular fates such as proliferation, differentiation and apoptosis. In particular, the dynamics of the Erk/Mek system, which has become the canonical example for MAPK signaling systems, have attracted considerable attention. Erk is encoded by two genes, Erk1 and Erk2, that until recently had been considered equivalent as they differ only subtly at the sequence level. However, these proteins exhibit radically different trafficking between cytoplasm and nucleus and this fact may have functional implications. Here we use spatially resolved data on Erk1/2 to develop and analyze spatio-temporal models of these cascades, and we discuss how sensitivity analysis can be used to discriminate between mechanisms. Our models elucidate some of the factors governing the interplay between signaling processes and the Erk1/2 localization in different cellular compartments, including competition between Erk1 and Erk2. Our approach is applicable to a wide range of signaling systems, such as activation cascades, where translocation of molecules occurs. Our study provides a first model of Erk1 and Erk2 activation and their nuclear shuttling dynamics, revealing a role in the regulation of the efficiency of nuclear signaling.

Affiliations:
Harrington H.A.-other affiliation
Komorowski M.-other affiliation
Beguerisse Díaz M.-other affiliation
Ratto G.M.-other affiliation
Stumpf M.P.H.-Imperial College London (GB)
22.Komorowski M., Costa M.J., Rand D.A., Stumpf M.P.H., Sensitivity, robustness, and identifiability in stochastic chemical kinetics models, PROCEEDINGS OF THE NATIONAL ACADEMY OF SCIENCES OF THE UNITED STATES OF AMERICA, ISSN: 0027-8424, DOI: 10.1073/pnas.1015814108, Vol.108, No.21, pp.8645-8650, 2011
Abstract:

We present a novel and simple method to numerically calculate Fisher information matrices for stochastic chemical kinetics models. The linear noise approximation is used to derive model equations and a likelihood function that leads to an efficient computational algorithm. Our approach reduces the problem of calculating the Fisher information matrix to solving a set of ordinary differential equations. This is the first method to compute Fisher information for stochastic chemical kinetics models without the need for Monte Carlo simulations. This methodology is then used to study sensitivity, robustness, and parameter identifiability in stochastic chemical kinetics models. We show that significant differences exist between stochastic and deterministic models as well as between stochastic models with time-series and time-point measurements. We demonstrate that these discrepancies arise from the variability in molecule numbers, correlations between species, and temporal correlations and show how this approach can be used in the analysis and design of experiments probing stochastic processes at the cellular level. The algorithm has been implemented as a Matlab package and is available from the authors upon request.

Affiliations:
Komorowski M.-other affiliation
Costa M.J.-University of Warwick (GB)
Rand D.A.-University of Warwick (GB)
Stumpf M.P.H.-Imperial College London (GB)
23.Komorowski M., Finkenstädt B., Rand D., Using a single fluorescent reporter gene to infer half-life of extrinsic noise and other parameters of gene expression, BIOPHYSICAL JOURNAL, ISSN: 0006-3495, DOI: 10.1016/j.bpj.2010.03.032, Vol.98, pp.2759-2769, 2010
Abstract:

Fluorescent and luminescent proteins are often used as reporters of transcriptional activity. Given the prevalence of noise in biochemical systems, the time-series data arising from these is of significant interest in efforts to calibrate stochastic models of gene expression and obtain information about sources of nongenetic variability. We present a statistical inference framework that can be used to estimate kinetic parameters of gene expression, as well as the strength and half-life of extrinsic noise from single fluorescent-reporter-gene time-series data. The method takes into account stochastic variability in a fluorescent signal resulting from intrinsic noise of gene expression, kinetics of fluorescent protein maturation, and extrinsic noise, which is assumed to arise at transcriptional level. We use the linear noise approximation and derive an explicit formula for the likelihood of observed fluorescent data. The method is embedded in a Bayesian paradigm, so that certain parameters can be informed from other experiments allowing portability of results across different studies. Inference is performed using Markov chain Monte Carlo. Fluorescent reporters are primary tools to observe dynamics of gene expression and the correct interpretation of fluorescent data is crucial to investigating these fundamental processes of cellular life. As both magnitude and frequency of the noise may have a dramatic effect on the cell fitness, the quantification of stochastic fluctuation is essential to the understanding of how genes are regulated. Our method provides a framework that addresses this important question.

Affiliations:
Komorowski M.-other affiliation
Finkenstädt B.-University of Warwick (GB)
Rand D.-University of Warwick (GB)

List of chapters in recent monographs
1.
587
Nienałtowski K., Jetka T., Komorowski M., Quantitative Biology. Theory, Computational Methods, and Models, rozdział: Sensitivity analysis, MIT Press, pp.293-319, 2018
2.
588
Jetka T., Nienałtowski K., Komorowski M., Quantitative Biology. Theory, Computational Methods, and Models, rozdział: Experimental design, MIT Press, pp.321-337, 2018

Conference papers
1.Vahdat Z., Nienałtowski K., Farooq Z., Komorowski M., Singh A., Information processing in unregulated and autoregulated gene expression, ECC20, European Control Conference, 2020-05-12/05-15, Saint Petersburg, virtual (RU), DOI: 10.23919/ECC51009.2020.9143689, pp.258-263, 2020
Abstract:

How living cells can reliably process biochemical cues in the presence of molecular noise is not fully understood. Here we investigate the fidelity of information transfer in the expression of a single gene. We use the established model of gene expression to examine how precisely the protein levels can be controlled by two distinct mechanisms: (i) the transcription rate of the gene, or (ii) the translation rate for the corresponding mRNA. The fidelity of gene expression is quantified with the information-theoretic notion of information capacity. Derived information capacity formulae reveal that transcriptional control generally provides a tangibly higher capacity as compared to the translational control. We next introduce negative feedback regulation in gene expression, where the protein directly inhibits its own transcription. While negative feedback reduces noise in the level of the protein for a given input signal, it also decreases the input-to-output sensitivity. Our results show that the combined effect of these two opposing forces is a reduced capacity in the presence of feedback. In summary, our analysis presents analytical quantification of information transfer in simple gene expression models, which provides insight into the fidelity of basic gene expression control mechanisms.

Affiliations:
Vahdat Z.-University of Delaware (US)
Nienałtowski K.-other affiliation
Farooq Z.-IPPT PAN
Komorowski M.-IPPT PAN
Singh A.-University of Delaware (US)

Conference abstracts
1.Topolewski P., Komorowski M., Cell cycle does not contribute to cell-to-cell heterogeneity of interferon responses, Cytokine, ISSN: 1043-4666, DOI: 10.1016/j.cyto.2017.09.011, Vol.100, pp.Mo-P7-11-100-100, 2017
2.Zakrzewska K.E., Jetka T., Nienałtowski K., Szymańska K., Andryka K., Topolewski P., Głów E., Komorowski M., Sensing and remembering IFNs concentrations, Cytokine, ISSN: 1043-4666, DOI: 10.1016/j.cyto.2017.09.011, Vol.100, pp.Mo-P7-12-100-100, 2017
3.Andryka K., Głów E., Nienałtowski K., Jetka T., Komorowski M., Sensing accuracy of interferons' concentrations, Cytokine, ISSN: 1043-4666, DOI: 10.1016/j.cyto.2015.08.238, Vol.76, pp.108, 2015
Abstract:

Interferons exhibit their key role of immune modulators through activation of the Jak-Stat signalling pathway. We know substantial amount of molecular details regarding functioning of the pathway. However, to what extend the action of the pathway is dose dependent at the single cell level remains largely unclear. Specifically it is not know if single cells respond in a digital fashion or their output is continuously dependent on the stimulant’s concentration. We have combined an information-theoretic framework with high-throughput confocal imaging of mouse embryonic fibroblasts to provide a thorough, single-cell analysis of the Jak-Stat signalling in response to interferon beta and interferon gamma. We showed that in a baseline state single cells have information hardly sufficient to distinguish between presence or absence of interferons. However they can be put in an alert state by an action of interferons, which allows them to respond more in an analogous fashion. Our results show that the accuracy with which signalling pathways transmit information is not fixed but can be modulated on the contextual basis.

Affiliations:
Andryka K.-other affiliation
Głów E.-IPPT PAN
Nienałtowski K.-IPPT PAN
Jetka T.-other affiliation
Komorowski M.-IPPT PAN
4.Jetka T., Komorowski M., How can we quantify ligand sensitivity for single-cell heterogenous dynamical responses?, FEBS Journal, ISSN: 1742-464X, Vol.281, No.Supplement: 1, Special Issue: SI, pp.625, 2014
Keywords:

Information Processing, Sensitivity Analysis, Signal Transduction

Affiliations:
Jetka T.-other affiliation
Komorowski M.-IPPT PAN
5.Głów E., Jetka T., Komorowski M., Sensitisation in the IFN-alpha/b, IFN-gamma crosstalk reveals mechanisms for enhanced information processing capacity of the STAT1, STAT2 signalling pathway, FEBS Journal, ISSN: 1742-464X, Vol.281, No.Supplement: 1, Special Issue: SI, pp.640, 2014
Keywords:

Interferon, Jak/Stat, Stimulation

Affiliations:
Głów E.-other affiliation
Jetka T.-other affiliation
Komorowski M.-IPPT PAN