@phdthesis{Stuebler2023,
  author    = {St{\"u}bler, Sabine},
  title     = {Mathematical model of the mucosal immune response to study inflammatory bowel diseases and their treatments},
  doi       = {10.25932/publishup-61230},
  url       = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-612301},
  school      = {Universit{\"a}t Potsdam},
  pages     = {xiv, 194},
  year      = {2023},
  abstract  = {Inflammatory bowel diseases (IBD), characterised by a chronic inflammation of the gut wall, develop as consequence of an overreacting immune response to commensal bacteria, caused by a combination of genetic and environmental conditions. Large inter-individual differences in the outcome of currently available therapies complicate the decision for the best option for an individual patient. Predicting the prospects of therapeutic success for an individual patient is currently only possible to a limited extent; for this, a better understanding of possible differences between responders and non-responders is needed. In this thesis, we have developed a mathematical model describing the most important processes of the gut mucosal immune system on the cellular level. The model is based on literature data, which were on the one hand used (qualitatively) to choose which cell types and processes to incorporate and to derive the model structure, and on the other hand (quantitatively) to derive the parameter values. Using ordinary differential equations, it describes the concentration-time course of neutrophils, macrophages, dendritic cells, T cells and bacteria, each subdivided into different cell types and activation states, in the lamina propria and mesenteric lymph nodes. We evaluate the model by means of simulations of the healthy immune response to salmonella infection and mucosal injury. A virtual population includes IBD patients, which we define through their initially asymptomatic, but after a trigger chronically inflamed gut wall. We demonstrate the model's usefulness in different analyses: (i) The comparison of virtual IBD patients with virtual healthy individuals shows that the disease is elicited by many small or fewer large changes, and allows to make hypotheses about dispositions relevant for development of the disease. (ii) We simulate the effects of different therapeutic targets and make predictions about the therapeutic outcome based on the pre-treatment state. (iii) From the analysis of differences between virtual responders and non-responders, we derive hypotheses about reasons for the inter-individual variability in treatment outcome. (iv) For the example of anti-TNF-alpha therapy, we analyse, which alternative therapies are most promising in case of therapeutic failure, and which therapies are most suited for combination therapies: For drugs also directly targeting the cytokine levels or inhibiting the recruitment of innate immune cells, we predict a low probability of success when used as alternative treatment, but a large gain when used in a combination treatment. For drugs with direct effects on T cells, via modulation of the sphingosine-1-phosphate receptor or inhibition of T cell proliferation, we predict a considerably larger probability of success when used as alternative treatment, but only a small additional gain when used in a combination therapy.},
  language  = {en}
}
@article{SteuerGrossSelbigetal.2006,
  author    = {Steuer, Ralf and Gross, Thilo and Selbig, Joachim and Blasius, Bernd},
  title     = {Structural kinetic modeling of metabolic networks},
  series = {Proceedings of the National Academy of Sciences of the United States of America},
  volume    = {103},
  journal   = {Proceedings of the National Academy of Sciences of the United States of America},
  number    = {32},
  publisher = {National Academy of Sciences},
  address   = {Washington},
  issn      = {0027-8424},
  doi       = {10.1073/pnas.0600013103},
  pages     = {11868 -- 11873},
  year      = {2006},
  abstract  = {To develop and investigate detailed mathematical models of metabolic processes is one of the primary challenges in systems biology. However, despite considerable advance in the topological analysis of metabolic networks, kinetic modeling is still often severely hampered by inadequate knowledge of the enzyme-kinetic rate laws and their associated parameter values. Here we propose a method that aims to give a quantitative account of the dynamical capabilities of a metabolic system, without requiring any explicit information about the functional form of the rate equations. Our approach is based on constructing a local linear model at each point in parameter space, such that each element of the model is either directly experimentally accessible or amenable to a straightforward biochemical interpretation. This ensemble of local linear models, encompassing all possible explicit kinetic models, then allows for a statistical exploration of the comprehensive parameter space. The method is exemplified on two paradigmatic metabolic systems: the glycolytic pathway of yeast and a realistic-scale representation of the photosynthetic Calvin cycle.},
  language  = {en}
}
@article{JueppnerMubeenLeisseetal.2017,
  author    = {J{\"u}ppner, Jessica and Mubeen, Umarah and Leisse, Andrea and Caldana, Camila and Brust, Henrike and Steup, Martin and Herrmann, Marion and Steinhauser, Dirk and Giavalisco, Patrick},
  title     = {Dynamics of lipids and metabolites during the cell cycle of Chlamydomonas reinhardtii},
  series = {The plant journal},
  volume    = {92},
  journal   = {The plant journal},
  publisher = {Wiley},
  address   = {Hoboken},
  issn      = {0960-7412},
  doi       = {10.1111/tpj.13642},
  pages     = {331 -- 343},
  year      = {2017},
  abstract  = {Metabolites and lipids are the final products of enzymatic processes, distinguishing the different cellular functions and activities of single cells or whole tissues. Understanding these cellular functions within a well-established model system requires a systemic collection of molecular and physiological information. In the current report, the green alga Chlamydomonas reinhardtii was selected to establish a comprehensive workflow for the detailed multi-omics analysis of a synchronously growing cell culture system. After implementation and benchmarking of the synchronous cell culture, a two-phase extraction method was adopted for the analysis of proteins, lipids, metabolites and starch from a single sample aliquot of as little as 10-15million Chlamydomonas cells. In a proof of concept study, primary metabolites and lipids were sampled throughout the diurnal cell cycle. The results of these time-resolved measurements showed that single compounds were not only coordinated with each other in different pathways, but that these complex metabolic signatures have the potential to be used as biomarkers of various cellular processes. Taken together, the developed workflow, including the synchronized growth of the photoautotrophic cell culture, in combination with comprehensive extraction methods and detailed metabolic phenotyping has the potential for use in in-depth analysis of complex cellular processes, providing essential information for the understanding of complex biological systems.},
  language  = {en}
}
@article{SchwahnBeleggiaOmranianetal.2017,
  author    = {Schwahn, Kevin and Beleggia, Romina and Omranian, Nooshin and Nikoloski, Zoran},
  title     = {Stoichiometric Correlation Analysis: Principles of Metabolic Functionality from Metabolomics Data},
  series = {Frontiers in plant science},
  volume    = {8},
  journal   = {Frontiers in plant science},
  publisher = {Frontiers Research Foundation},
  address   = {Lausanne},
  issn      = {1664-462X},
  doi       = {10.3389/fpls.2017.02152},
  pages     = {12},
  year      = {2017},
  abstract  = {Recent advances in metabolomics technologies have resulted in high-quality (time-resolved) metabolic profiles with an increasing coverage of metabolic pathways. These data profiles represent read-outs from often non-linear dynamics of metabolic networks. Yet, metabolic profiles have largely been explored with regression-based approaches that only capture linear relationships, rendering it difficult to determine the extent to which the data reflect the underlying reaction rates and their couplings. Here we propose an approach termed Stoichiometric Correlation Analysis (SCA) based on correlation between positive linear combinations of log-transformed metabolic profiles. The log-transformation is due to the evidence that metabolic networks can be modeled by mass action law and kinetics derived from it. Unlike the existing approaches which establish a relation between pairs of metabolites, SCA facilitates the discovery of higherorder dependence between more than two metabolites. By using a paradigmatic model of the tricarboxylic acid cycle we show that the higher-order dependence reflects the coupling of concentration of reactant complexes, capturing the subtle difference between the employed enzyme kinetics. Using time-resolved metabolic profiles from Arabidopsis thaliana and Escherichia coli, we show that SCA can be used to quantify the difference in coupling of reactant complexes, and hence, reaction rates, underlying the stringent response in these model organisms. By using SCA with data from natural variation of wild and domesticated wheat and tomato accession, we demonstrate that the domestication is accompanied by loss of such couplings, in these species. Therefore, application of SCA to metabolomics data from natural variation in wild and domesticated populations provides a mechanistic way to understanding domestication and its relation to metabolic networks.},
  language  = {en}
}
@article{GorochowskiAycilarKucukgozeBovenbergetal.2016,
  author    = {Gorochowski, Thomas E. and Aycilar-Kucukgoze, Irem and Bovenberg, Roel A. L. and Roubos, Johannes A. and Ignatova, Zoya},
  title     = {A Minimal Model of Ribosome Allocation Dynamics Captures Trade-offs in Expression between Endogenous and Synthetic Genes},
  series = {ACS synthetic biology},
  volume    = {5},
  journal   = {ACS synthetic biology},
  publisher = {American Chemical Society},
  address   = {Washington},
  issn      = {2161-5063},
  doi       = {10.1021/acssynbio.6b00040},
  pages     = {710 -- 720},
  year      = {2016},
  abstract  = {Cells contain a finite set of resources that must be distributed across many processes to ensure survival. Among them, the largest proportion of cellular resources is dedicated to protein translation. Synthetic biology often exploits these resources in executing orthogonal genetic circuits, yet the burden this places on the cell is rarely considered. Here, we develop a minimal model of ribosome allocation dynamics capturing the demands on translation when expressing a synthetic construct together with endogenous genes required for the maintenance of cell physiology. Critically, it contains three key variables related to design parameters of the synthetic construct covering transcript abundance, translation initiation rate, and elongation time. We show that model-predicted changes in ribosome allocation closely match experimental shifts in synthetic protein expression rate and cellular growth. Intriguingly, the model is also able to accurately infer transcript levels and translation times after further exposure to additional ambient stress. Our results demonstrate that a simple model of resource allocation faithfully captures the redistribution of protein synthesis resources when faced with the burden of synthetic gene expression and environmental stress. The tractable nature of the model makes it a versatile tool for exploring the guiding principles of efficient heterologous expression and the indirect interactions that can arise between synthetic circuits and their host chassis because of competition for shared translational resources.},
  language  = {en}
}
@misc{KaminskiSchaubSiegeletal.2013,
  author    = {Kaminski, Roland and Schaub, Torsten H. and Siegel, Anne and Videla, Santiago},
  title     = {Minimal intervention strategies in logical signaling networks with ASP},
  series = {Postprints der Universit{\"a}t Potsdam : Mathematisch Naturwissenschaftliche Reihe},
  journal   = {Postprints der Universit{\"a}t Potsdam : Mathematisch Naturwissenschaftliche Reihe},
  number    = {4-5},
  issn      = {1866-8372},
  doi       = {10.25932/publishup-41570},
  url       = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-415704},
  pages     = {675 -- 690},
  year      = {2013},
  abstract  = {Proposing relevant perturbations to biological signaling networks is central to many problems in biology and medicine because it allows for enabling or disabling certain biological outcomes. In contrast to quantitative methods that permit fine-grained (kinetic) analysis, qualitative approaches allow for addressing large-scale networks. This is accomplished by more abstract representations such as logical networks. We elaborate upon such a qualitative approach aiming at the computation of minimal interventions in logical signaling networks relying on Kleene's three-valued logic and fixpoint semantics. We address this problem within answer set programming and show that it greatly outperforms previous work using dedicated algorithms.},
  language  = {en}
}
@phdthesis{Schwahn2018,
  author    = {Schwahn, Kevin},
  title     = {Data driven approaches to infer the regulatory mechanism shaping and constraining levels of metabolites in metabolic networks},
  doi       = {10.25932/publishup-42324},
  url       = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-423240},
  school      = {Universit{\"a}t Potsdam},
  pages     = {109},
  year      = {2018},
  abstract  = {Systems biology aims at investigating biological systems in its entirety by gathering and analyzing large-scale data sets about the underlying components. Computational systems biology approaches use these large-scale data sets to create models at different scales and cellular levels. In addition, it is concerned with generating and testing hypotheses about biological processes. However, such approaches are inevitably leading to computational challenges due to the high dimensionality of the data and the differences in the dimension of data from different cellular layers. This thesis focuses on the investigation and development of computational approaches to analyze metabolite profiles in the context of cellular networks. This leads to determining what aspects of the network functionality are reflected in the metabolite levels. With these methods at hand, this thesis aims to answer three questions: (1) how observability of biological systems is manifested in metabolite profiles and if it can be used for phenotypical comparisons; (2) how to identify couplings of reaction rates from metabolic profiles alone; and (3) which regulatory mechanism that affect metabolite levels can be distinguished by integrating transcriptomics and metabolomics read-outs. I showed that sensor metabolites, identified by an approach from observability theory, are more correlated to each other than non-sensors. The greater correlations between sensor metabolites were detected both with publicly available metabolite profiles and synthetic data simulated from a medium-scale kinetic model. I demonstrated through robustness analysis that correlation was due to the position of the sensor metabolites in the network and persisted irrespectively of the experimental conditions. Sensor metabolites are therefore potential candidates for phenotypical comparisons between conditions through targeted metabolic analysis. Furthermore, I demonstrated that the coupling of metabolic reaction rates can be investigated from a purely data-driven perspective, assuming that metabolic reactions can be described by mass action kinetics. Employing metabolite profiles from domesticated and wild wheat and tomato species, I showed that the process of domestication is associated with a loss of regulatory control on the level of reaction rate coupling. I also found that the same metabolic pathways in Arabidopsis thaliana and Escherichia coli exhibit differences in the number of reaction rate couplings. I designed a novel method for the identification and categorization of transcriptional effects on metabolism by combining data on gene expression and metabolite levels. The approach determines the partial correlation of metabolites with control by the principal components of the transcript levels. The principle components contain the majority of the transcriptomic information allowing to partial out the effect of the transcriptional layer from the metabolite profiles. Depending whether the correlation between metabolites persists upon controlling for the effect of the transcriptional layer, the approach allows us to group metabolite pairs into being associated due to post-transcriptional or transcriptional regulation, respectively. I showed that the classification of metabolite pairs into those that are associated due to transcriptional or post-transcriptional regulation are in agreement with existing literature and findings from a Bayesian inference approach. The approaches developed, implemented, and investigated in this thesis open novel ways to jointly study metabolomics and transcriptomics data as well as to place metabolic profiles in the network context. The results from these approaches have the potential to provide further insights into the regulatory machinery in a biological system.},
  language  = {en}
}
@phdthesis{RobainaEstevez2017,
  author    = {Robaina Estevez, Semidan},
  title     = {Context-specific metabolic predictions},
  url       = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-401365},
  school      = {Universit{\"a}t Potsdam},
  pages     = {vi, 158},
  year      = {2017},
  abstract  = {All life-sustaining processes are ultimately driven by thousands of biochemical reactions occurring in the cells: the metabolism. These reactions form an intricate network which produces all required chemical compounds, i.e., metabolites, from a set of input molecules. Cells regulate the activity through metabolic reactions in a context-specific way; only reactions that are required in a cellular context, e.g., cell type, developmental stage or environmental condition, are usually active, while the rest remain inactive. The context-specificity of metabolism can be captured by several kinds of experimental data, such as by gene and protein expression or metabolite profiles. In addition, these context-specific data can be assimilated into computational models of metabolism, which then provide context-specific metabolic predictions. This thesis is composed of three individual studies focussing on context-specific experimental data integration into computational models of metabolism. The first study presents an optimization-based method to obtain context-specific metabolic predictions, and offers the advantage of being fully automated, i.e., free of user defined parameters. The second study explores the effects of alternative optimal solutions arising during the generation of context-specific metabolic predictions. These alternative optimal solutions are metabolic model predictions that represent equally well the integrated data, but that can markedly differ. This study proposes algorithms to analyze the space of alternative solutions, as well as some ways to cope with their impact in the predictions. Finally, the third study investigates the metabolic specialization of the guard cells of the plant Arabidopsis thaliana, and compares it with that of a different cell type, the mesophyll cells. To this end, the computational methods developed in this thesis are applied to obtain metabolic predictions specific to guard cell and mesophyll cells. These cell-specific predictions are then compared to explore the differences in metabolic activity between the two cell types. In addition, the effects of alternative optima are taken into consideration when comparing the two cell types. The computational results indicate a major reorganization of the primary metabolism in guard cells. These results are supported by an independent 13C labelling experiment.},
  language  = {en}
}
@phdthesis{Breuer2016,
  author    = {Breuer, David},
  title     = {The plant cytoskeleton as a transportation network},
  url       = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-93583},
  school      = {Universit{\"a}t Potsdam},
  pages     = {164},
  year      = {2016},
  abstract  = {The cytoskeleton is an essential component of living cells. It is composed of different types of protein filaments that form complex, dynamically rearranging, and interconnected networks. The cytoskeleton serves a multitude of cellular functions which further depend on the cell context. In animal cells, the cytoskeleton prominently shapes the cell's mechanical properties and movement. In plant cells, in contrast, the presence of a rigid cell wall as well as their larger sizes highlight the role of the cytoskeleton in long-distance intracellular transport. As it provides the basis for cell growth and biomass production, cytoskeletal transport in plant cells is of direct environmental and economical relevance. However, while knowledge about the molecular details of the cytoskeletal transport is growing rapidly, the organizational principles that shape these processes on a whole-cell level remain elusive. This thesis is devoted to the following question: How does the complex architecture of the plant cytoskeleton relate to its transport functionality? The answer requires a systems level perspective of plant cytoskeletal structure and transport. To this end, I combined state-of-the-art confocal microscopy, quantitative digital image analysis, and mathematically powerful, intuitively accessible graph-theoretical approaches. This thesis summarizes five of my publications that shed light on the plant cytoskeleton as a transportation network: (1) I developed network-based frameworks for accurate, automated quantification of cytoskeletal structures, applicable in, e.g., genetic or chemical screens; (2) I showed that the actin cytoskeleton displays properties of efficient transport networks, hinting at its biological design principles; (3) Using multi-objective optimization, I demonstrated that different plant cell types sustain cytoskeletal networks with cell-type specific and near-optimal organization; (4) By investigating actual transport of organelles through the cell, I showed that properties of the actin cytoskeleton are predictive of organelle flow and provided quantitative evidence for a coordination of transport at a cellular level; (5) I devised a robust, optimization-based method to identify individual cytoskeletal filaments from a given network representation, allowing the investigation of single filament properties in the network context. The developed methods were made publicly available as open-source software tools. Altogether, my findings and proposed frameworks provide quantitative, system-level insights into intracellular transport in living cells. Despite my focus on the plant cytoskeleton, the established combination of experimental and theoretical approaches is readily applicable to different organisms. Despite the necessity of detailed molecular studies, only a complementary, systemic perspective, as presented here, enables both understanding of cytoskeletal function in its evolutionary context as well as its future technological control and utilization.},
  language  = {en}
}
@phdthesis{Mauri2014,
  author    = {Mauri, Marco},
  title     = {A model for sigma factor competition in bacterial cells},
  url       = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-72098},
  school      = {Universit{\"a}t Potsdam},
  pages     = {167},
  year      = {2014},
  abstract  = {Bacteria respond to changing environmental conditions by switching the global pattern of expressed genes. In response to specific environmental stresses the cell activates several stress-specific molecules such as sigma factors. They reversibly bind the RNA polymerase to form the so-called holoenzyme and direct it towards the appropriate stress response genes. In exponentially growing E. coli cells, the majority of the transcriptional activity is carried out by the housekeeping sigma factor, while stress responses are often under the control of alternative sigma factors. Different sigma factors compete for binding to a limited pool of RNA polymerase (RNAP) core enzymes, providing a mechanism for cross talk between genes or gene classes via the sharing of expression machinery. To quantitatively analyze the contribution of sigma factor competition to global changes in gene expression, we develop a thermodynamic model that describes binding between sigma factors and core RNAP at equilibrium, transcription, non-specific binding to DNA and the modulation of the availability of the molecular components. Association of housekeeping sigma factor to RNAP is generally favored by its abundance and higher binding affinity to the core. In order to promote transcription by alternative sigma subunits, the bacterial cell modulates the transcriptional efficiency in a reversible manner through several strategies such as anti-sigma factors, 6S RNA and generally any kind of transcriptional regulators (e.g. activators or inhibitors). By shifting the outcome of sigma factor competition for the core, these modulators bias the transcriptional program of the cell. The model is validated by comparison with in vitro competition experiments, with which excellent agreement is found. We observe that transcription is affected via the modulation of the concentrations of the different types of holoenzymes, so saturated promoters are only weakly affected by sigma factor competition. However, in case of overlapping promoters or promoters recognized by two types of sigma factors, we find that even saturated promoters are strongly affected. Active transcription effectively lowers the affinity between the sigma factor driving it and the core RNAP, resulting in complex cross talk effects and raising the question of how their in vitro measure is relevant in the cell. We also estimate that sigma factor competition is not strongly affected by non-specific binding of core RNAPs, sigma factors, and holoenzymes to DNA. Finally, we analyze the role of increased core RNAP availability upon the shut-down of ribosomal RNA transcription during stringent response. We find that passive up-regulation of alternative sigma-dependent transcription is not only possible, but also displays hypersensitivity based on the sigma factor competition. Our theoretical analysis thus provides support for a significant role of passive control during that global switch of the gene expression program and gives new insights into RNAP partitioning in the cell.},
  language  = {en}
}