Modeling, predicting and mining metabolism at atom-level resolution

Living organisms can catalyze many thousands of biochemical reactions that they use to convert energy and matter, which provides them with the essentials for life. The sum of these chemical reactions happening in an organism is called metabo-lism. Understanding metabolism is crucial to elucidate the fundamental principles of biology, and further to enable us to redesign it for the sustainable, biosynthetic pro-duction of bio-based fuels, commodity chemicals and medicines. Mathematical mod-els are essential to organize and understand the complexity of metabolism. They usu-ally represent metabolism as a network of reactions, but they tend to neglect the exact molecular structure of the metabolites. To redesign metabolic reactions, how-ever, a mechanistic understanding of metabolic reactions and their catalysts, proteins called enzymes, is essential. In this work, a mathematical description of enzymatic reaction mechanism, called generalized reaction rules, is applied to computationally simulate and predict meta-bolic processes at the level of atoms. Each reaction rule describes the catalytic activi-ty of an enzyme, or a group of enzymes, at the mechanistic level by encoding the re-arrangement of atoms in the reaction. The reaction rules are called âgeneralizedâ, because they mimic the ability of a single enzyme to catalyze multiple reactions by acting on a range of substrates. Using these reaction rules, we first developed a computational representation of me-tabolism that allowed tracking single atoms throughout complex metabolic reaction networks. The principle of atom-tracking was then used to develop a graph-theory based method to represent and analyze metabolic networks, and to reliably identify of metabolic pathways for the biosynthesis of chemicals. Next, we applied the gener-alized reaction rules to predict all possible novel, hypothetical reactions from known biological compounds, and we stored the five million generated novel reactions them in a database called ATLAS. Finally, the developed tools and resources were applied to specific engineering and research problems, such as the biosynthetic pathway design for the biofuel bisabolene and the plastic precursor 1,4-butanediol. We further pre-dicted a biosynthesis route for the pharmaceutical tetrahydropalmatine and engi-neered a yeast strain to produce it. Finally, we show that our tools can be used to mine available genome sequences to find organisms that can degrade xenobiotics. Our findings suggest that the atom-level representation of metabolism can greatly contribute to its understanding, exploration and prediction. Given the complexity of atom-level modeling of metabolic processes, we propose metrics that can approxi-mate the atom-level information to conserve the information at the level of big, hy-pothetical metabolic networks like ATLAS. This database plus the developed pathway search techniques form a valuable resource for scientist to help characterizing un-known biosynthesis pathways towards secondary metabolites, and for metabolic engi-neers for designing novel bioproduction pathways for chemicals. Hopefully, these considerations will contribute to a better understanding of metabolism, advance the exploration of the bioproduction of drugs and other valuable molecules, and acceler-ate metabolic engineering efforts to realize the switch from a petroleum-based chem-ical industry towards a more sustainable, bio-based production of societyâs chemical needs.

Algorithms and flowchart for the design of synthetic biochemical networks

Homa Mohammadi Peyhani

New drugs are needed to assure effective therapies for previously untreated diseases, emerging diseases, and personalized medicine, but the process of drug development is complex, costly, and time-consuming. This is especially problematic considering that 90% of drug candidates in clinical trials are discarded due to unexpected toxicity or other secondary effects. This inefficiency threatens our health care system and economy. Despite the advances in the cellular metabolism, our knowledge of the mechanisms governing enzymatic biotransformations in cells is far from complete, in particular regarding degradation pathways, mode of action, or side effects of drugs. Examining the mechanisms of enzymatic reactions at the cellular scale could improve our fundamental understanding of their catalytic capability, and facilitate identifying and filling the knowledge gaps. The scale and the complexity of metabolic data is ever-expanding, requiring scientists to apply more advanced computational methods to systematically store, explore, and interpret the enzymatic potential of cells. The first step toward simulating enzymes in silico is to learn from their biochemical reactions in nature. To do this, we use distilled knowledge of known biochemistry in the form of generalized enzymatic reaction rules. Enzymatic rules are mathematical representations of enzymatic action mimicking the catalytic function of enzymes. They are formulated in a less specific manner (more promiscuous) to act on a broad range of substrates. In addition to reconstructing known biochemistry, the application of these reaction rules paves the way toward the discovery of novel enzymatic interactions. In this thesis, I developed computational models, tools, and methodologies to facilitate the study of metabolism and catalytic action of enzymes. We analyzed different aspects of metabolism through five distinct studies: In a first study, in order to provide a holistic view of currently known biochemistry, we gathered biochemical data from 14 sources, covering the known metabolic networks of all species. We integrated all biological data into a high-performance database based on ontology, named LCSB DB. We further expanded the scope of LCSB DB to cover all bioactive and chemicals. LCSB DB offers fast and efficient searching of biochemical data and serves as a platform for sharing, storing, and analyzing biochemical data. In a second study, we used enzymatic reaction rules to predict all theoretically possible metabolic reactions between biological and bioactive compounds in LCSB DB. In a third study, we developed a method to find enzymes are able to catalyze orphan and predicted reactions, called BridgIT. BridgIT uses the knowledge of reactive sites on substrates to find the most similar, known biochemical reactions. We then validated the utility of BridgIT in enzyme discovery for the design of de novo synthetic pathways producing tetrahydropalmatine and adipic acid. In the last study, we propose a workflow for rational drug design and systems-level analysis of drug metabolism, called NICEdrug.ch. NICEdrug.ch allows large-scale computational analysis of drug biochemistry (metabolic precursors or prodrugs and metabolic fate or degradation), enzymatic targets, and toxicity in the context of cellular metabolism. Finally, in the conclusion chapter, we discuss the contribution and the potential further applications of the computational tools that were developed in this thesis.

EPFL2020

Kinetic reconstruction and analysis of sphingolipid metabolism

Vassily Hatzimanikatis, Howard Riezman, Georgios Savoglidis

Lipids are major constituents of the cell. They are responsible for major properties of the cellular membranes: hydrophobicity, selective permeability and being the scaffold of signaling proteins. Many diseases are associated with alterations in the lipid distribution in the cell and the composition of membrane domains. Metabolic syndrome, obesity, atherosclerosis, as well as Alzheimer’s, Huntington’s diseases and cancer, have an impact in the levels of lipids, with observed alterations in their concentrations compared to the healthy state. Applying computational techniques along with systematic modeling of lipid metabolism can provide insights that can guide biomedical research and develop potential strategies for prevention and cure. In the present study we developed a comprehensive model of the sphingolipid biosynthesis in the yeast Saccharomyces cerevisiae. Sphingolipids are one of the four major lipid categories, along with (glycero)phospholipids, sterols and fatty acids synthesized in the yeast S. cerevisiae. The importance of sphingolipids present in any higher eukaryote has been demonstrated in many recent studies. For this study, S. cerevisiae has been chosen as a model organism due to its high homology of cellular processes with mammalian cells. The developed model will be an essential part towards the construction of a detailed kinetic model of the whole lipidome of the cell. We first constructed a stoichiometric model that contains all the currently known reactions for the biosynthesis of ceramides and complex sphingolipids in the yeast. Additionally, we have accounted for all five reported hydroxylation states, along with the reactions synthesizing the necessary precursor metabolites from other lipid pathways (i.e. palmitate-CoA from fatty acid synthesis and phosphatidylinositol from the phospholipids metabolism). We next developed a kinetic model and we used a large number of lipidomic measurements of wild type yeast to consistently calibrate our model. Curation of the kinetic information of the model came from the comprehensive mining of references for operation of the enzymes as well as ranges of kinetic parameters from online databases. These datasets created the pool for a sampling technique that accounts for the uncertainty in the parameters of the model. This lead to a robust dynamic model, containing mass balances for all the components of the biochemical network as well as terms that accounted for the dillution of these molecules in the cell due to growth. We performed a thorough kinetic analysis of the system by examining the impact of different assumptions in enzyme operation on the levels of ceramides and complex sphingolipids (e.g. substrate competition, (un)competitive inhibition). By applying the principles of Metabolic Control Analysis (MCA) we were able to quantify the effect of enzyme activities on the lipid profiles and we identified enzymes of the biochemical reaction network as efficient targets for metabolic engineering towards a desired state. We also applied a reverse engineering method and we were able to identify enzyme perturbations responsible for an observed altered state compared to the wild type cellular lipidomic profile. Such an approach could lead to the identification of genetic mutations by imploring the information containd within metabolomic measurements. Although we demonstrate the application of the method in models of lipid metabolism it is possible to be applied to biochemical systems of various parts of metabolism and sizes creating a modular platform for kinetic analysis of cellular operations.

2014

Discovery and Evaluation of Novel Pathways for Production of the Second Generation of Biofuels

Meriç Ataman, Noushin Hadadi, Vassily Hatzimanikatis, Ljubisa Miskovic, Milenko Tokic

In an effort of overcoming the limited availability of fossil energy resources the focus of the research and development in the area of biofuels has moved towards developing the second generation of fuels that should be produced via microbial fermentation. The idea is to use as a feedstock inexpensive and abundant waste materials such as lignocellulosic biomass. The second generation biofuels should satisfy several criteria such as lower emission, higher energy density and should be less corrosive to engines. However, currently used industrial workhorses such as E. coli and S. cerevisiae do not produce many of these molecules naturally. Furthermore, for majority of these molecules there are no known biochemical pathways. The need for discovery of novel biosynthetic pathways towards desired molecules sparked the development of computational tools that are capable of reconstructing feasible reaction steps between a given set of starting compounds (precursors) and a molecule of interest (a target compound). In this study, we performed the retrobiosynthesis analysis for all the candidates applying Biochemical Network Integrated Computational Explorer (BNICE.ch) and we reconstructed the metabolic network for 69 fuel compounds. We analyzed compounds (biological or chemical) and reactions one-step away from these target molecules and the appearance of the potential substrates in their reconstructed metabolic network. Based on the results of this analysis, coupled with other criteria such as the Gibes free energy of formation, combustion potential and expert opinion, we chose the 5 highest ranked fuel candidates for further pathway reconstruction studies. For these 5 fuel candidates we have created 1.8 million de novo pathways with different number of reaction steps for further evaluation. We have then embedded the reconstructed pathways in the genome-°©‐scale models of the two chassis organisms E. coli and P. putida and we performed the pathway evaluation with respect to their maximum theoretical yield using two methods (i) Flux Balance Analysis (FBA), and (ii) Thermodynamic-based Flux Analysis (TFA). All the pathways that did not satisfy the thermodynamic constraints were rejected. Our results indicated that in the evaluation and ranking of the pathways the thermodynamics is a necessary consideration. Specifically, we showed that the majority of the pathways that were determined as feasible based on FBA, they were not feasible based on TFA. Moreover, we also showed that some pathways that were determined as feasible in both organisms based on FBA, they were not feasible in both of them based on TFA. Also, the yield of the feasible pathways can differ depending on the organism. These results demonstrate how TFA can also guide the selection of a suitable chassis organism. For all discovered reactions in the feasible pathways we have also identified known enzymes from the KEGG database with the closest reaction similarity that could be engineered for the implementation of novel reactions. This study shows the full potential of computational tools for discovery and design of novel synthetic pathways and their relevance for future developments in the area of metabolic engineering and synthetic biology.

2016