@phdthesis{Arnold2014, author = {Arnold, Anne}, title = {Modeling photosynthesis and related metabolic processes : from detailed examination to consideration of the metabolic context}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus-72277}, school = {Universit{\"a}t Potsdam}, year = {2014}, abstract = {Mathematical modeling of biological systems is a powerful tool to systematically investigate the functions of biological processes and their relationship with the environment. To obtain accurate and biologically interpretable predictions, a modeling framework has to be devised whose assumptions best approximate the examined scenario and which copes with the trade-off of complexity of the underlying mathematical description: with attention to detail or high coverage. Correspondingly, the system can be examined in detail on a smaller scale or in a simplified manner on a larger scale. In this thesis, the role of photosynthesis and its related biochemical processes in the context of plant metabolism was dissected by employing modeling approaches ranging from kinetic to stoichiometric models. The Calvin-Benson cycle, as primary pathway of carbon fixation in C3 plants, is the initial step for producing starch and sucrose, necessary for plant growth. Based on an integrative analysis for model ranking applied on the largest compendium of (kinetic) models for the Calvin-Benson cycle, those suitable for development of metabolic engineering strategies were identified. Driven by the question why starch rather than sucrose is the predominant transitory carbon storage in higher plants, the metabolic costs for their synthesis were examined. The incorporation of the maintenance costs for the involved enzymes provided a model-based support for the preference of starch as transitory carbon storage, by only exploiting the stoichiometry of synthesis pathways. Many photosynthetic organisms have to cope with processes which compete with carbon fixation, such as photorespiration whose impact on plant metabolism is still controversial. A systematic model-oriented review provided a detailed assessment for the role of this pathway in inhibiting the rate of carbon fixation, bridging carbon and nitrogen metabolism, shaping the C1 metabolism, and influencing redox signal transduction. The demand of understanding photosynthesis in its metabolic context calls for the examination of the related processes of the primary carbon metabolism. To this end, the Arabidopsis core model was assembled via a bottom-up approach. This large-scale model can be used to simulate photoautotrophic biomass production, as an indicator for plant growth, under so-called optimal, carbon-limiting and nitrogen-limiting growth conditions. Finally, the introduced model was employed to investigate the effects of the environment, in particular, nitrogen, carbon and energy sources, on the metabolic behavior. This resulted in a purely stoichiometry-based explanation for the experimental evidence for preferred simultaneous acquisition of nitrogen in both forms, as nitrate and ammonium, for optimal growth in various plant species. The findings presented in this thesis provide new insights into plant system's behavior, further support existing opinions for which mounting experimental evidences arise, and posit novel hypotheses for further directed large-scale experiments.}, language = {en} } @phdthesis{Basler2012, author = {Basler, Georg}, title = {Mass-balanced randomization : a significance measure for metabolic networks}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus-62037}, school = {Universit{\"a}t Potsdam}, year = {2012}, abstract = {Complex networks have been successfully employed to represent different levels of biological systems, ranging from gene regulation to protein-protein interactions and metabolism. Network-based research has mainly focused on identifying unifying structural properties, including small average path length, large clustering coefficient, heavy-tail degree distribution, and hierarchical organization, viewed as requirements for efficient and robust system architectures. Existing studies estimate the significance of network properties using a generic randomization scheme - a Markov-chain switching algorithm - which generates unrealistic reactions in metabolic networks, as it does not account for the physical principles underlying metabolism. Therefore, it is unclear whether the properties identified with this generic approach are related to the functions of metabolic networks. Within this doctoral thesis, I have developed an algorithm for mass-balanced randomization of metabolic networks, which runs in polynomial time and samples networks almost uniformly at random. The properties of biological systems result from two fundamental origins: ubiquitous physical principles and a complex history of evolutionary pressure. The latter determines the cellular functions and abilities required for an organism's survival. Consequently, the functionally important properties of biological systems result from evolutionary pressure. By employing randomization under physical constraints, the salient structural properties, i.e., the smallworld property, degree distributions, and biosynthetic capabilities of six metabolic networks from all kingdoms of life are shown to be independent of physical constraints, and thus likely to be related to evolution and functional organization of metabolism. This stands in stark contrast to the results obtained from the commonly applied switching algorithm. In addition, a novel network property is devised to quantify the importance of reactions by simulating the impact of their knockout. The relevance of the identified reactions is verified by the findings of existing experimental studies demonstrating the severity of the respective knockouts. The results suggest that the novel property may be used to determine the reactions important for viability of organisms. Next, the algorithm is employed to analyze the dependence between mass balance and thermodynamic properties of Escherichia coli metabolism. The thermodynamic landscape in the vicinity of the metabolic network reveals two regimes of randomized networks: those with thermodynamically favorable reactions, similar to the original network, and those with less favorable reactions. The results suggest that there is an intrinsic dependency between thermodynamic favorability and evolutionary optimization. The method is further extended to optimizing metabolic pathways by introducing novel chemically feasibly reactions. The results suggest that, in three organisms of biotechnological importance, introduction of the identified reactions may allow for optimizing their growth. The approach is general and allows identifying chemical reactions which modulate the performance with respect to any given objective function, such as the production of valuable compounds or the targeted suppression of pathway activity. These theoretical developments can find applications in metabolic engineering or disease treatment. The developed randomization method proposes a novel approach to measuring the significance of biological network properties, and establishes a connection between large-scale approaches and biological function. The results may provide important insights into the functional principles of metabolic networks, and open up new possibilities for their engineering.}, language = {en} } @article{BaslerGrimbsEbenhoehetal.2012, author = {Basler, Georg and Grimbs, Sergio and Ebenh{\"o}h, Oliver and Selbig, Joachim and Nikoloski, Zoran}, title = {Evolutionary significance of metabolic network properties}, series = {Interface : journal of the Royal Society}, volume = {9}, journal = {Interface : journal of the Royal Society}, number = {71}, publisher = {Royal Society}, address = {London}, issn = {1742-5689}, doi = {10.1098/rsif.2011.0652}, pages = {1168 -- 1176}, year = {2012}, abstract = {Complex networks have been successfully employed to represent different levels of biological systems, ranging from gene regulation to protein-protein interactions and metabolism. Network-based research has mainly focused on identifying unifying structural properties, such as small average path length, large clustering coefficient, heavy-tail degree distribution and hierarchical organization, viewed as requirements for efficient and robust system architectures. However, for biological networks, it is unclear to what extent these properties reflect the evolutionary history of the represented systems. Here, we show that the salient structural properties of six metabolic networks from all kingdoms of life may be inherently related to the evolution and functional organization of metabolism by employing network randomization under mass balance constraints. Contrary to the results from the common Markov-chain switching algorithm, our findings suggest the evolutionary importance of the small-world hypothesis as a fundamental design principle of complex networks. The approach may help us to determine the biologically meaningful properties that result from evolutionary pressure imposed on metabolism, such as the global impact of local reaction knockouts. Moreover, the approach can be applied to test to what extent novel structural properties can be used to draw biologically meaningful hypothesis or predictions from structure alone.}, language = {en} } @phdthesis{Grimbs2009, author = {Grimbs, Sergio}, title = {Towards structure and dynamics of metabolic networks}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus-32397}, school = {Universit{\"a}t Potsdam}, year = {2009}, abstract = {This work presents mathematical and computational approaches to cover various aspects of metabolic network modelling, especially regarding the limited availability of detailed kinetic knowledge on reaction rates. It is shown that precise mathematical formulations of problems are needed i) to find appropriate and, if possible, efficient algorithms to solve them, and ii) to determine the quality of the found approximate solutions. Furthermore, some means are introduced to gain insights on dynamic properties of metabolic networks either directly from the network structure or by additionally incorporating steady-state information. Finally, an approach to identify key reactions in a metabolic networks is introduced, which helps to develop simple yet useful kinetic models. The rise of novel techniques renders genome sequencing increasingly fast and cheap. In the near future, this will allow to analyze biological networks not only for species but also for individuals. Hence, automatic reconstruction of metabolic networks provides itself as a means for evaluating this huge amount of experimental data. A mathematical formulation as an optimization problem is presented, taking into account existing knowledge and experimental data as well as the probabilistic predictions of various bioinformatical methods. The reconstructed networks are optimized for having large connected components of high accuracy, hence avoiding fragmentation into small isolated subnetworks. The usefulness of this formalism is exemplified on the reconstruction of the sucrose biosynthesis pathway in Chlamydomonas reinhardtii. The problem is shown to be computationally demanding and therefore necessitates efficient approximation algorithms. The problem of minimal nutrient requirements for genome-scale metabolic networks is analyzed. Given a metabolic network and a set of target metabolites, the inverse scope problem has as it objective determining a minimal set of metabolites that have to be provided in order to produce the target metabolites. These target metabolites might stem from experimental measurements and therefore are known to be produced by the metabolic network under study, or are given as the desired end-products of a biotechological application. The inverse scope problem is shown to be computationally hard to solve. However, I assume that the complexity strongly depends on the number of directed cycles within the metabolic network. This might guide the development of efficient approximation algorithms. Assuming mass-action kinetics, chemical reaction network theory (CRNT) allows for eliciting conclusions about multistability directly from the structure of metabolic networks. Although CRNT is based on mass-action kinetics originally, it is shown how to incorporate further reaction schemes by emulating molecular enzyme mechanisms. CRNT is used to compare several models of the Calvin cycle, which differ in size and level of abstraction. Definite results are obtained for small models, but the available set of theorems and algorithms provided by CRNT can not be applied to larger models due to the computational limitations of the currently available implementations of the provided algorithms. Given the stoichiometry of a metabolic network together with steady-state fluxes and concentrations, structural kinetic modelling allows to analyze the dynamic behavior of the metabolic network, even if the explicit rate equations are not known. In particular, this sampling approach is used to study the stabilizing effects of allosteric regulation in a model of human erythrocytes. Furthermore, the reactions of that model can be ranked according to their impact on stability of the steady state. The most important reactions in that respect are identified as hexokinase, phosphofructokinase and pyruvate kinase, which are known to be highly regulated and almost irreversible. Kinetic modelling approaches using standard rate equations are compared and evaluated against reference models for erythrocytes and hepatocytes. The results from this simplified kinetic models can simulate acceptably the temporal behavior for small changes around a given steady state, but fail to capture important characteristics for larger changes. The aforementioned approach to rank reactions according to their influence on stability is used to identify a small number of key reactions. These reactions are modelled in detail, including knowledge about allosteric regulation, while all other reactions were still described by simplified reaction rates. These so-called hybrid models can capture the characteristics of the reference models significantly better than the simplified models alone. The resulting hybrid models might serve as a good starting point for kinetic modelling of genome-scale metabolic networks, as they provide reasonable results in the absence of experimental data, regarding, for instance, allosteric regulations, for a vast majority of enzymatic reactions.}, language = {en} } @phdthesis{Kueken2020, author = {K{\"u}ken, Anika}, title = {Predictions from constraint-based approaches including enzyme kinetics}, school = {Universit{\"a}t Potsdam}, pages = {116, A-16, B-7, C-8}, year = {2020}, abstract = {The metabolic state of an organism reflects the entire phenotype that is jointly affected by genetic and environmental changes. Due to the complexity of metabolism, system-level modelling approaches have become indispensable tools to obtain new insights into biological functions. In particular, simulation and analysis of metabolic networks using constraint-based modelling approaches have helped the analysis of metabolic fluxes. However, despite ongoing improvements in prediction of reaction flux through a system, approaches to directly predict metabolite concentrations from large-scale metabolic networks remain elusive. In this thesis, we present a computational approach for inferring concentration ranges from genome-scale metabolic models endowed with mass action kinetics. The findings specify a molecular mechanism underling facile control of concentration ranges for components in large-scale metabolic networks. Most importantly, an extended version of the approach can be used to predict concentration ranges without knowledge of kinetic parameters, provided measurements of concentrations in a reference state. We show that the approach is applicable with large-scale kinetic and stoichiometric metabolic models of organisms from different kingdoms of life. By challenging the predictions of concentration ranges in the genome-scale metabolic network of Escherichia coli with real-world data sets, we further demonstrate the prediction power and limitations of the approach. To predict concentration ranges in other species, e.g. model plant species Arabidopsis thaliana, we would rely on estimates of kinetic parameters (i.e. enzyme catalytic rates) since plant-specific enzyme catalytic rates are poorly documented. Using the constraint-based approach of Davidi et al. for estimation of enzyme catalytic rates, we obtain values for 168 plant enzymes. The approach depends on quantitative proteomics data and flux estimates obtained from constraint-based model of plant leaf metabolism integrating maximal rates of selected enzymes, plant-specific constraints on fluxes through canonical pathways, and growth measurements from Arabidopsis thaliana rosette under ten conditions. We demonstrate a low degree of plant enzyme saturation, supported by the agreement between concentrations of nicotinamide adenine dinucleotide, adenosine triphosphate, and glyceraldehyde 3-phosphate, based on our maximal in vivo catalytic rates, and available quantitative metabolomics data. Hence, our results show genome-wide estimation for plant-specific enzyme catalytic rates is feasible. These can now be readily employed to study resource allocation, to predict enzyme and metabolite concentrations using recent constrained-based modelling approaches. Constraint-based methods do not directly account for kinetic mechanisms and corresponding parameters. Therefore, a number of workflows have already been proposed to approximate reaction kinetics and to parameterize genome-scale kinetic models. We present a systems biology strategy to build a fully parameterized large-scale model of Chlamydomonas reinhardtii accounting for microcompartmentalization in the chloroplast stroma. Eukaryotic algae comprise a microcompartment, the pyrenoid, essential for the carbon concentrating mechanism (CCM) that improves their photosynthetic performance. Since the experimental study of the effects of microcompartmentation on metabolic pathways is challenging, we employ our model to investigate compartmentation of fluxes through the Calvin-Benson cycle between pyrenoid and stroma. Our model predicts that ribulose-1,5-bisphosphate, the substrate of Rubisco, and 3-phosphoglycerate, its product, diffuse in and out of the pyrenoid. We also find that there is no major diffusional barrier to metabolic flux between the pyrenoid and stroma. Therefore, our computational approach represents a stepping stone towards understanding of microcompartmentalized CCM in other organisms. This thesis provides novel strategies to use genome-scale metabolic networks to predict and integrate metabolite concentrations. Therefore, the presented approaches represent an important step in broadening the applicability of large-scale metabolic models to a range of biotechnological and medical applications.}, language = {en} } @misc{OmranianKleessenTohgeetal.2015, author = {Omranian, Nooshin and Kleessen, Sabrina and Tohge, Takayuki and Klie, Sebastian and Basler, Georg and M{\"u}ller-R{\"o}ber, Bernd and Fernie, Alisdair R. and Nikoloski, Zoran}, title = {Differential metabolic and coexpression networks of plant metabolism}, series = {Trends in plant science}, volume = {20}, journal = {Trends in plant science}, number = {5}, publisher = {Elsevier}, address = {London}, issn = {1360-1385}, doi = {10.1016/j.tplants.2015.02.002}, pages = {266 -- 268}, year = {2015}, abstract = {Recent analyses have demonstrated that plant metabolic networks do not differ in their structural properties and that genes involved in basic metabolic processes show smaller coexpression than genes involved in specialized metabolism. By contrast, our analysis reveals differences in the structure of plant metabolic networks and patterns of coexpression for genes in (non)specialized metabolism. Here we caution that conclusions concerning the organization of plant metabolism based on network-driven analyses strongly depend on the computational approaches used.}, 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} }