TY - JOUR A1 - Barlow, Axel A1 - Hartmann, Stefanie A1 - Gonzalez, Javier A1 - Hofreiter, Michael A1 - Paijmans, Johanna L. A. T1 - Consensify BT - a method for generating pseudohaploid genome sequences from palaeogenomic datasets with reduced error rates JF - Genes / Molecular Diversity Preservation International N2 - A standard practise in palaeogenome analysis is the conversion of mapped short read data into pseudohaploid sequences, frequently by selecting a single high-quality nucleotide at random from the stack of mapped reads. This controls for biases due to differential sequencing coverage, but it does not control for differential rates and types of sequencing error, which are frequently large and variable in datasets obtained from ancient samples. These errors have the potential to distort phylogenetic and population clustering analyses, and to mislead tests of admixture using D statistics. We introduce Consensify, a method for generating pseudohaploid sequences, which controls for biases resulting from differential sequencing coverage while greatly reducing error rates. The error correction is derived directly from the data itself, without the requirement for additional genomic resources or simplifying assumptions such as contemporaneous sampling. For phylogenetic and population clustering analysis, we find that Consensify is less affected by artefacts than methods based on single read sampling. For D statistics, Consensify is more resistant to false positives and appears to be less affected by biases resulting from different laboratory protocols than other frequently used methods. Although Consensify is developed with palaeogenomic data in mind, it is applicable for any low to medium coverage short read datasets. We predict that Consensify will be a useful tool for future studies of palaeogenomes. KW - palaeogenomics KW - ancient DNA KW - sequencing error KW - error reduction KW - D statistics KW - bioinformatics Y1 - 2020 U6 - https://doi.org/10.3390/genes11010050 SN - 2073-4425 VL - 11 IS - 1 PB - MDPI CY - Basel ER - TY - JOUR A1 - Christopher Ashwood, Wout Bittremieux A1 - Bittremieux, Wout A1 - Deutsch, Eric W. A1 - Doncheva, Nadezhda T. A1 - Dorfer, Viktoria A1 - Gabriels, Ralf A1 - Gorshkov, Vladimir A1 - Gupta, Surya A1 - Jones, Andrew R. A1 - Käll, Lukas A1 - Kopczynski, Dominik A1 - Lane, Lydie A1 - Lautenbacher, Ludwig A1 - Legeay, Marc A1 - Locard-Paulet, Marie A1 - Mesuere, Bart A1 - Sachsenberg, Timo A1 - Salz, Renee A1 - Samaras, Patroklos A1 - Schiebenhoefer, Henning A1 - Schmidt, Tobias A1 - Schwämmle, Veit A1 - Soggiu, Alessio A1 - Uszkoreit, Julian A1 - Van Den Bossche, Tim A1 - Van Puyvelde, Bart A1 - Van Strien, Joeri A1 - Verschaffelt, Pieter A1 - Webel, Henry A1 - Willems, Sander A1 - Perez-Riverolab, Yasset A1 - Netz, Eugen A1 - Pfeuffer, Julianus T1 - Proceedings of the EuBIC-MS 2020 Developers’ Meeting JF - EuPA Open Proteomics N2 - The 2020 European Bioinformatics Community for Mass Spectrometry (EuBIC-MS) Developers’ meeting was held from January 13th to January 17th 2020 in Nyborg, Denmark. Among the participants were scientists as well as developers working in the field of computational mass spectrometry (MS) and proteomics. The 4-day program was split between introductory keynote lectures and parallel hackathon sessions. During the latter, the participants developed bioinformatics tools and resources addressing outstanding needs in the community. The hackathons allowed less experienced participants to learn from more advanced computational MS experts, and to actively contribute to highly relevant research projects. We successfully produced several new tools that will be useful to the proteomics community by improving data analysis as well as facilitating future research. All keynote recordings are available on https://doi.org/10.5281/zenodo.3890181. KW - computational mass spectrometry KW - proteomics KW - bioinformatics KW - spectrum clustering KW - phosphoproteomics KW - XIC extraction KW - proteomics graph networks KW - predicted spectra Y1 - 2020 U6 - https://doi.org/10.1016/j.euprot.2020.11.001 SN - 2212-9685 VL - 24 SP - 1 EP - 6 PB - Elsevier CY - Amsterdam ER - TY - JOUR A1 - Luo, Ting A1 - Chen, Xiaoyi A1 - Zeng, Shufei A1 - Guan, Baozhang A1 - Hu, Bo A1 - Meng, Yu A1 - Liu, Fanna A1 - Wong, Taksui A1 - Lu, Yongpin A1 - Yun, Chen A1 - Hocher, Berthold A1 - Yin, Lianghong T1 - Bioinformatic identification of key genes and analysis of prognostic values in clear cell renal cell carcinoma JF - Oncology Letters N2 - The present study aimed to identify new key genes as potential biomarkers for the diagnosis, prognosis or targeted therapy of clear cell renal cell carcinoma (ccRCC). Three expression profiles (GSE36895, GSE46699 and GSE71963) were collected from Gene Expression Omnibus. GEO2R was used to identify differentially expressed genes (DEGs) in ccRCC tissues and normal samples. The Database for Annotation, Visualization and Integrated Discovery was utilized for functional and pathway enrichment analysis. STRING v10.5 and Molecular Complex Detection were used for protein-protein interaction (PPI) network construction and module analysis, respectively. Regulation network analyses were performed with the WebGestal tool. UALCAN web-portal was used for expression validation and survival analysis of hub genes in ccRCC patients from The Cancer Genome Atlas (TCGA). A total of 65 up- and 164 downregulated genes were identified as DEGs. DEGs were enriched with functional terms and pathways compactly related to ccRCC pathogenesis. Seventeen hub genes and one significant module were filtered out and selected from the PPI network. The differential expression of hub genes was verified in TCGA patients. Kaplan-Meier plot showed that high mRNA expression of enolase 2 (ENO2) was associated with short overall survival in ccRCC patients (P=0.023). High mRNA expression of cyclin D1 (CCND1) (P<0.001), fms related tyrosine kinase 1 (FLT1) (P=0.004), plasminogen (PLG) (P<0.001) and von Willebrand factor (VWF) (P=0.008) appeared to serve as favorable factors in survival. These findings indicate that the DEGs may be key genes in ccRCC pathogenesis and five genes, including ENO2, CCND1, PLT1, PLG and VWF, may serve as potential prognostic biomarkers in ccRCC. KW - clear cell renal cell carcinoma KW - bioinformatics KW - differentially expressed genes KW - biomarkers KW - Kaplan-Meier plot Y1 - 2018 U6 - https://doi.org/10.3892/ol.2018.8842 SN - 1792-1074 SN - 1792-1082 VL - 16 IS - 2 SP - 1747 EP - 1757 PB - Spandidos publ LTD CY - Athens ER - TY - JOUR A1 - Frioux, Clémence A1 - Schaub, Torsten H. A1 - Schellhorn, Sebastian A1 - Siegel, Anne A1 - Wanko, Philipp T1 - Hybrid metabolic network completion JF - Theory and practice of logic programming N2 - Metabolic networks play a crucial role in biology since they capture all chemical reactions in an organism. While there are networks of high quality for many model organisms, networks for less studied organisms are often of poor quality and suffer from incompleteness. To this end, we introduced in previous work an answer set programming (ASP)-based approach to metabolic network completion. Although this qualitative approach allows for restoring moderately degraded networks, it fails to restore highly degraded ones. This is because it ignores quantitative constraints capturing reaction rates. To address this problem, we propose a hybrid approach to metabolic network completion that integrates our qualitative ASP approach with quantitative means for capturing reaction rates. We begin by formally reconciling existing stoichiometric and topological approaches to network completion in a unified formalism. With it, we develop a hybrid ASP encoding and rely upon the theory reasoning capacities of the ASP system dingo for solving the resulting logic program with linear constraints over reals. We empirically evaluate our approach by means of the metabolic network of Escherichia coli. Our analysis shows that our novel approach yields greatly superior results than obtainable from purely qualitative or quantitative approaches. KW - answer set programming KW - metabolic network KW - gap-filling KW - linear programming KW - hybrid solving KW - bioinformatics Y1 - 2018 U6 - https://doi.org/10.1017/S1471068418000455 SN - 1471-0684 SN - 1475-3081 VL - 19 IS - 1 SP - 83 EP - 108 PB - Cambridge University Press CY - New York ER - TY - JOUR A1 - Kuckelkorn, Ulrike A1 - Stübler, Sabine A1 - Textoris-Taube, Kathrin A1 - Kilian, Christiane A1 - Niewienda, Agathe A1 - Henklein, Petra A1 - Janek, Katharina A1 - Stumpf, Michael P. H. A1 - Mishto, Michele A1 - Liepe, Juliane T1 - Proteolytic dynamics of human 20S thymoproteasome JF - The journal of biological chemistry N2 - An efficient immunosurveillance of CD8(+) T cells in the periphery depends on positive/negative selection of thymocytes and thus on the dynamics of antigen degradation and epitope production by thymoproteasome and immunoproteasome in the thymus. Although studies in mouse systems have shown how thymoproteasome activity differs from that of immunoproteasome and strongly impacts the T cell repertoire, the proteolytic dynamics and the regulation of human thymoproteasome are unknown. By combining biochemical and computational modeling approaches, we show here that human 20S thymoproteasome and immunoproteasome differ not only in the proteolytic activity of the catalytic sites but also in the peptide transport. These differences impinge upon the quantity of peptide products rather than where the substrates are cleaved. The comparison of the two human 20S proteasome isoforms depicts different processing of antigens that are associated to tumors and autoimmune diseases. KW - proteasome KW - protein degradation KW - antigen processing KW - computational biology KW - bioinformatics KW - thymoproteasome KW - thymus KW - proteolysis Y1 - 2019 U6 - https://doi.org/10.1074/jbc.RA118.007347 SN - 1083-351X VL - 294 IS - 19 SP - 7740 EP - 7754 PB - American Society for Biochemistry and Molecular Biology CY - Bethesda ER - TY - GEN A1 - Barlow, Axel A1 - Hartmann, Stefanie A1 - Gonzalez, Javier A1 - Hofreiter, Michael A1 - Paijmans, Johanna L. A. T1 - Consensify BT - a method for generating pseudohaploid genome sequences from palaeogenomic datasets with reduced error rates T2 - Postprints der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe N2 - A standard practise in palaeogenome analysis is the conversion of mapped short read data into pseudohaploid sequences, frequently by selecting a single high-quality nucleotide at random from the stack of mapped reads. This controls for biases due to differential sequencing coverage, but it does not control for differential rates and types of sequencing error, which are frequently large and variable in datasets obtained from ancient samples. These errors have the potential to distort phylogenetic and population clustering analyses, and to mislead tests of admixture using D statistics. We introduce Consensify, a method for generating pseudohaploid sequences, which controls for biases resulting from differential sequencing coverage while greatly reducing error rates. The error correction is derived directly from the data itself, without the requirement for additional genomic resources or simplifying assumptions such as contemporaneous sampling. For phylogenetic and population clustering analysis, we find that Consensify is less affected by artefacts than methods based on single read sampling. For D statistics, Consensify is more resistant to false positives and appears to be less affected by biases resulting from different laboratory protocols than other frequently used methods. Although Consensify is developed with palaeogenomic data in mind, it is applicable for any low to medium coverage short read datasets. We predict that Consensify will be a useful tool for future studies of palaeogenomes. T3 - Zweitveröffentlichungen der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe - 1033 KW - palaeogenomics KW - ancient DNA KW - sequencing error KW - error reduction KW - D statistics KW - bioinformatics Y1 - 2020 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-472521 SN - 1866-8372 IS - 1033 ER - TY - GEN A1 - Gebser, Martin A1 - Schaub, Torsten H. A1 - Thiele, Sven A1 - Veber, Philippe T1 - Detecting inconsistencies in large biological networks with answer set programming T2 - Postprints der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe N2 - We introduce an approach to detecting inconsistencies in large biological networks by using answer set programming. To this end, we build upon a recently proposed notion of consistency between biochemical/genetic reactions and high-throughput profiles of cell activity. We then present an approach based on answer set programming to check the consistency of large-scale data sets. Moreover, we extend this methodology to provide explanations for inconsistencies by determining minimal representations of conflicts. In practice, this can be used to identify unreliable data or to indicate missing reactions. T3 - Zweitveröffentlichungen der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe - 561 KW - answer set programming KW - bioinformatics KW - consistency KW - diagnosis Y1 - 2019 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-412467 SN - 1866-8372 IS - 561 ER - TY - THES A1 - Robaina Estevez, Semidan T1 - Context-specific metabolic predictions T1 - Kontextspezifische metabolische Vorhersagen BT - computational methods and applications BT - Berechnungsmethoden und Anwendungen N2 - 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. N2 - Alle lebenserhaltenden Prozesse werden durch tausende biochemische Reaktionen in der Zelle bestimmt, welche den Metabolismus charakterisieren. Diese Reaktionen bilden ein komplexes Netzwerk, welches alle notwendigen chemischen Verbindungen, die sogenannten Metabolite, aus einer bestimmten Menge an Ausgangsmolekülen produziert Zellen regulieren ihren Stoffwechsel kontextspezifisch, dies bedeutet, dass nur Reaktionen die in einem zellulären Kontext, zum Beispiel Zelltyp, Entwicklungsstadium oder verschiedenen Umwelteinflüssen, benötigt werden auch tatsächlich aktiv sind. Die übrigen Reaktionen werden als inaktiv betrachtet. Die Kontextspezifität des Metabolismus kann durch verschiedene experimentelle Daten, wie Gen- und Proteinexpressionen oder Metabolitprofile erfasst werden. Zusätzlich können diese Daten in Computersimulationen des Metabolismus integriert werden, um kontextspezifische (metabolische) Vorhersagen zu treffen. Diese Doktorarbeit besteht aus drei unabhängigen Studien, welche die Integration von kontextspezifischen experimentellen Daten in Computersimulationen des Metabolismus thematisieren. Die erste Studie beschreibt ein Konzept, basierend auf einem mathematischen Optimierungsproblem, welches es erlaubt kontextspezifische, metabolische Vorhersagen zu treffen. Dabei bietet diese vollautomatische Methode den Vorteil vom Nutzer unabhängige Parameter, zu verwenden. Die zweite Studie untersucht den Einfluss von alternativen optimalen Lösungen, welche bei kontextspezifischen metabolischen Vorhersagen generiert werden. Diese alternativen Lösungen stellen metabolische Modellvorhersagen da, welche die integrierten Daten gleichgut wiederspiegeln, sich aber grundlegend voneinander unterscheiden können. Diese Studie zeigt verschiedene Ansätze alternativen Lösungen zu analysieren und ihren Einfluss auf die Vorhersagen zu berücksichtigen. Schlussendlich, untersucht die dritte Studie die metabolische Spezialisierung der Schließzellen in Arabidopsis thaliana und vergleicht diese mit einer weiteren Zellart, den Mesophyllzellen. Zu diesem Zweck wurden die in dieser Doktorarbeit vorgestellten Methoden angewandt um metabolische Vorhersagen speziell für Schließzellen und Mesophyllzellen zu erhalten. Anschließend wurden die zellspezifischen Vorhersagen auf Unterschiede in der metabolischen Aktivität der Zelltypen, unter Berücksichtigung des Effekt von alternativen Optima, untersucht. Die Ergebnisse der Simulationen legen eine grundlegende Neuorganisation des Primärmetabolismus in Schließzellen verglichen mit Mesophyllzellen nahe. Diese Ergebnisse werden durch unabhängige 13C markierungs Experimente bestätigt. KW - systems biology KW - bioinformatics KW - metabolic networks KW - constraint-based modeling KW - data integration KW - Systemsbiologie KW - Bioinformatik KW - Stoffwechselnetze KW - Constraint-basierte Modellierung KW - Datenintegration Y1 - 2017 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-401365 ER - TY - JOUR A1 - Gebser, Martin A1 - Schaub, Torsten H. A1 - Thiele, Sven A1 - Veber, Philippe T1 - Detecting inconsistencies in large biological networks with answer set programming JF - Theory and practice of logic programming N2 - We introduce an approach to detecting inconsistencies in large biological networks by using answer set programming. To this end, we build upon a recently proposed notion of consistency between biochemical/genetic reactions and high-throughput profiles of cell activity. We then present an approach based on answer set programming to check the consistency of large-scale data sets. Moreover, we extend this methodology to provide explanations for inconsistencies by determining minimal representations of conflicts. In practice, this can be used to identify unreliable data or to indicate missing reactions. KW - answer set programming KW - bioinformatics KW - consistency KW - diagnosis Y1 - 2011 U6 - https://doi.org/10.1017/S1471068410000554 SN - 1471-0684 VL - 11 IS - 5-6 SP - 323 EP - 360 PB - Cambridge Univ. Press CY - New York ER - TY - THES A1 - Grimbs, Sergio T1 - Towards structure and dynamics of metabolic networks T1 - Struktur und Dynamik metabolischer Netzwerke N2 - 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. N2 - In dieser Arbeit werden mathematische und informatische Ansätze zur Behandlung diverser Probleme im Zusammenhang mit der Modellierung metabolischer Netzwerke vorgestellt, insbesondere unter Berücksichtigung der eingeschränkten Verfügbarkeit detaillierter Enzymkinetiken. Es wird gezeigt, dass präzise mathematische Formulierungen der Probleme notwendig sind, um erstens angemessene und, falls möglich, effiziente Algorithmen zur Lösung zu entwickeln. Und zweitens, um die Güte der so gefundenen Lösungen zu bewerten. Des weiteren werden Methoden zur Analyse dynamischer Eigenschaften metabolischer Netzwerke eingeführt, welche entweder nur auf der Struktur der Netzwerke basieren oder zusätzlich noch Informationen über stationäre Zustände mit berücksichtigen. Außerdem wird eine Strategie zur Bestimmung von Schlüsselreaktionen eines Netzwerkes vorgestellt, welche die Entwicklung kinetischer Modelle vereinfacht. Der Erfolg neuer Technologien ermöglicht eine immer billigere und schnellere Sequenzierung des Genoms. Dies wird in naher Zukunft die Analyse biologischer Netzwerke nicht nur für Spezies, sondern auch für einzelne Individuen ermöglichen. Die automatische Rekonstruktion metabolischer Netzwerke ist bestens dafür geeignet, diese großen Datenmengen auszuwerten. Eine mathematische Formulierung der Rekonstruktion als Optimierungsproblem wird vorgestellt, die sowohl bereits vorhandenes Wissen als auch theoretische Vorhersagen verschiedenster bioinformatischer Methoden berücksichtigt. Die rekonstruierten Netzwerke sind hinsichtlich möglichst großer und plausibler Zusammenhangskomponenten hin optimiert, um fragmentierte und isolierte Teilnetzwerke zu vermeiden. Als Beispiel dient die Rekonstruktion der Saccharosesynthese in Chlamydomonas reinhardtii. Es wird gezeigt, dass das Problem sehr rechenintensiv ist und somit Approximationsalgorithmen erforderlich macht. Das 'inverse scope' Problem hat als Optimierungsziel, für ein gegebenes metabolisches Netzwerk die minimale Menge notwendiger Metabolite zu bestimmen, um eine ebenfalls gegebene Menge von gewünschten Zielmetaboliten zu produzieren. Diese Zielmetabolite können entweder durch experimentellen Messungen festgelegt werden, oder sie sind die gewünschten Endprodukte einer biotechnologischen Anwendung. Es wird gezeigt, dass das 'inverse scope' Problem rechenintensiv ist. Allerdings wird angenommen, dass die Berechnungskomplexität stark von der Anzahl gerichteter Zyklen innerhalb des metabolischen Netzwerkes abhängt. Dies könnte die Entwicklung effizienter Approximationsalgorithmen ermöglichen. Unter der Annahme von Massenwirkungskinetiken erlaubt es die 'chemical reaction network theory' (CRNT), anhand der Struktur metabolischer Netzwerke Rückschlüsse auf Multistabilität zu ziehen. Auch weitere Kinetiken können durch Modellierung von Enzymmechanismen mit berücksichtigt werden. CRNT wird zum Vergleich von mehreren Modellen des Calvinzyklus, welche sich in Größe und Abstraktionsniveau unterscheiden, verwendet. Obwohl für kleinere Modelle Ergebnisse erzielt werden, erlauben es die verfügbaren Theoreme und Algorithmen der CRNT nicht, Aussagen für größere Modelle zu machen, da die gegenwärtigen Implementierungen der Algorithmen an ihre Berechnungsgrenzen stoßen. Sind sowohl die Stoichiometrie eines metabolischen Netzwerkes, als auch die Metabolitkonzentrationen und Flüsse im stationären Zustand bekannt, so kann 'structural kinetic modelling' angewandt werden, um das dynamische Verhalten des Netzwerkes zu analysieren, selbst wenn die expliziten Ratengleichung unbekannt sind. Dieser Ansatz wird verwendet, um den stabilisierenden Einfluss allosterischer Regulation in menschlichen Erythrozyten zu untersuchen. Des weiteren werden die Reaktionen anhand ihrer Bedeutung hinsichtlich Stabilität im stationären Zustand angeordnet. Die wichtigsten Reaktionen bezüglich dieser Ordnung sind Hexokinase, Phosphofructokinase und Pyruvatkinase, welche bekanntermaßen stark reguliert und irreversibel sind. Kinetische Modelle, die auf generischen Ratengleichung beruhen, werden mit detaillierten Referenzmodellen für Erythrozyten und Hepatozyten verglichen. Die generischen Modelle simulieren das Verhalten nur in der Nähe eines gegebenen stationären Zustandes recht gut. Der zuvor erwähnte Ansatz, wichtige Reaktionen bezüglich Stabilität zu identifizieren, wird zur Bestimmung von Schlüsselreaktionen genutzt. Diese Schlüsselreaktionen werden im Detail modelliert, während für alle anderen Reaktionen weiterhin generische Ratengleichung verwendet werden. Die so entstandenen Hybridmodelle können das Verhalten des Referenzmodells signifikant besser beschreiben. Die Hybridmodelle können als Ausgangspunkt zur Erstellung genomweiter kinetischer Modelle dienen. KW - metabolische Netzwerke KW - Modellierung KW - Struktur KW - Dynamik KW - Bioinformatik KW - metabolic networks KW - modelling KW - structure KW - dynamics KW - bioinformatics Y1 - 2009 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus-32397 ER -