Carbon labeling for Metabolic Flux Analysis. Nicholas Wayne Henderson Computational and Applied Mathematics. Abstract

Size: px
Start display at page:

Download "Carbon labeling for Metabolic Flux Analysis. Nicholas Wayne Henderson Computational and Applied Mathematics. Abstract"

Transcription

1 Carbon labeling for Metabolic Flux Analysis Nicholas Wayne Henderson Computational and Applied Mathematics Abstract Metabolic engineering is rapidly growing as a discipline. Applications in pharmaceuticals, energy, and materials have the potential for a large commercial impact. This paper covers the basics of a mathematical toolset that uses carbon labeling information to assist in metabolic flux analysis. The purpose is to introduce users to the inner workings of a software tool developed in Matlab. It is hoped that this work makes the process and mathematics of Metabolic Flux Analysis easier to understand and implement.

2 Background Through studies in molecular biology and genetics, researchers have amassed significant knowledge of cellular metabolism and the controlling genetic regulation systems. Metabolic engineers strive to use this knowledge to manipulate metabolic reaction networks in order to achieve some objective. The text Metabolic Engineering (Stephanopoulos 1998) presents the following definition: metabolic engineering [is] the directed improvement of product formation or cellular properties through the modification of specific biochemical reactions or the introduction of new ones with the use of recombinant DNA technology. The key work to be carried out in metabolic engineering is selecting the reactions for modification. Once the reactions are targeted, established recombinant techniques are used to inhibit or amplify the genes corresponding to the necessary enzymes. Metabolic engineers use Metabolic Flux Analysis (MFA) in order to target the specific reactions. The flux of a metabolic reaction is the rate of change from input metabolites to output metabolites. As an example, one might write the following reaction: A + B" v C A,B,C # metabolites v # flux The flux, v, measures the rate at which metabolites A and B are consumed and C is produced. The purpose of MFA is to learn the fluxes of as many reactions in a metabolic pathway as possible. This knowledge gives the engineer an idea of which reactions are more active in a pathway. The engineer then has the ability to better predict the outcome of manipulations and can chose ones that help achieve his objective. Uncovering the intracellular fluxes is a complicated process. We want to know the flux of reactions deep inside a pathway, where the pools of intermediate metabolites are too small to measure. Basic stoichiometry, combined with knowledge of extracellular conditions, can be used to build a system of equations in order to solve for the flux values. The catch is that these systems are often underdetermined because there are more fluxes than metabolites in a pathway. Methods to solve this problem involve the use of 13 C labeling and knowledge of carbon skeleton rearrangements [1]. If an engineer knows the labeling patterns of input metabolites and is able to observe the labeling patterns of products then he is able to determine intracellular fluxes with greater certainty. The Metabolic Engineering research group in the Computational and Applied Mathematics department is currently interested in ways that 13 C labeling experiments can augment the underdetermined stoichiometric system in order to fully solve for the intracellular fluxes. 2

3 Mathematics of Metabolic Flux Analysis The literature on carbon labeling experiments for metabolic flux analysis is quite complex. This stems from a confusing notation that must handle various fluxes, metabolites, and other constraints. This section will attempt to convey the mathematics behind this project in order to educate users to the inner workings of the Matlab toolset. Care will be taken to keep the notation simple. Readers are encouraged to explore the literature for analysis with greater depth. Example System This section will use an example from Wiechert s paper, titled Bidirectional Reaction Steps in Metabolic Networks: III. Explicit Solution and Analysis of Isotopomer Labeling System. It is a vastly simplified version of the citric acid cycle: A f " 1 # B f B 2,b $ " 2 # E f B + E 3,b $ " 3 # C E f " 4 # H C f 6 "# D + F D f 7 "# E + G Figure 1: The example network. Reactions 2 and 3 are bidirectional. Figure 2: The reaction list corresponding to the example metabolic network in Figure 1. Flux variables are denoted f for forward and b for backward. Here A is the input metabolite. B, C, D, and E are intracellular metabolites. K, F, and H are output metabolites. Experimenters are typically able to measure only the fluxes corresponding to input and output metabolites. Thus, we must use modeling to gain information regarding the intracellular metabolites. It should be noted that the reaction fluxes are simply the rate at which educts are consumed and products are formed. Flux values should not be considered constants in the sense of reaction kinetics. Yet, metabolic systems with constant inputs will lead to constant flux values. If certain components of the system are modified, such as input metabolite or enzyme concentrations, the flux values will change. Stationarity of Intracellular Metabolites The key assumption in metabolic flux analysis is that the concentrations of intracellular metabolites are in a steady state. Thus, build up or depletion of these 3

4 metabolites does not occur with in a cell. They are created as fast as they are consumed. This is called a dynamic equilibrium and can be expressed by the following equations: d[b], d[c], d[d], d[e] [1] This allows the modeling of the system in terms of stoichiometric balance equations. Stoichiometry Each reaction is associated with a flux. Bidirectional reactions have 2 fluxes: one for the forward direction and one for the reverse direction. The flux describes the rate at which educts are turned into products. Referring back to the example, metabolite B is created at a rate of f 1 + b 2 + b 3 and consumed at a rate of f 2 + f 3 + f 5. The following balance equations results when considering all intracellular metabolites: d[b] = f 1 " f 2 + b 2 " f 3 + b 3 " f 5 d[c] = f 3 " b 3 " f 6 d[d] = f 6 " f 7 d[e] = f 2 " b 2 " f 3 + b 3 " f 4 + f 7 [2] Which can be written as a linear system: or, # 1 "1 1 "1 1 0 "1 0 0& # 0& % ( % ( % "1 0 0 "1 0 ( v = % 0 ( % "1( % 0( % ( % ( $ 0 1 "1 "1 1 " ' $ 0' Sv [3] where v = [ f 1, f 2,b 2, f 3,b 3, f 4, f 5, f 6, f 7 ] " is the vector of fluxes. This system is underdetermined, which is a common feature of realistic networks. There are typically many more fluxes than metabolites. The overall goal is to uncover the flux vector v. The stoichiometric balance system imposes constraints on v, but cannot be used to obtain a unique solution. We plan on using carbon labeling information to add more constraints to this system. 4

5 It is possible to measure some extracellular fluxes, those that enter or exit the system with no other involvement. In the example, let s say fluxes f 1 and f 5 are directly measurable. This information can be used to write a flux measurement system: or, " % " $ ' v = w % 1 $ ' # & # & w 2 Hv = w. [4] The combination of the stoichiometric [eq. 3] and flux [eq. 4] measurement systems reduces the dimensionality of the flux identification problem. First, an optimization method is used to find an initial solution (v 0 ) to the underdetermined system: " S% $ # H& ' v = " 0 % $ # w ' & that forces all of the fluxes to be positive (v 0 > 0). Second, Gaussian elimination or the singular value decomposition is used to find the null space of the augmented matrix ([S;H]). The result is the expression of all fluxes in terms of a few free fluxes: v = Nx + v 0 [5] where v is the desired vector of fluxes, the columns of N span the null space, x is the vector of free fluxes, and v 0 is the initial solution. In the example there are 9 fluxes in total, but only 3 free fluxes. Because of stoichiometry and flux measurement, only 3 flux values need to be uncovered in order to solve the entire system. In order to move forward, we must use information gained from carbon labeling experiments. Isotopomers Carbons atoms are naturally found in 3 forms: 12 C with 6 protons and 6 neutrons has a natural relative abundance of 98.9%, 13 C with 7 neutrons has a natural relative abundance of 1.1%, and 14 C with 8 neutrons exists in trace amounts ( 14 C is the radioactive isotope used in carbon dating). 13 C is useful in MFA, because it can be distinguished from 12 C through techniques such as NMR and mass spectrometry. 5

6 The metabolites of interest are typically molecules based on a carbon chain with 1 to 6 atoms of length. The word isotopomer describes a particular isotope isomer, or a carbon based molecule with a specific labeling pattern. For an example let s look at a metabolite with 3 carbons: Unlabeled 1 13 C atom 2 13 C atoms Fully Labeled Figure 3: Light circles represent C 12 atoms. Dark circles represent C 13 atoms. There are a total of 8 isotopomers (2 3 ) with distinct labeling patterns in a molecule with a 3-carbon chain (assuming the molecule is not symmetric). Isotopomer measurements combined with knowledge of the carbon transitions in a metabolic network can help us gain information about the flux values. To this extent, a system of differential equations is formed to model the network. The problem is that these equations are non-linear and cannot be solved analytically. Nevertheless, they are helpful. Wiechert analyzed the balance equations and found that isotopomer systems are stable. After a period of time, the isotopomer fractions reach steady values. This allows us to think of isotopomers in terms of balance equations. Wiechert did more work to show that the isotopomer balance equations could be transformed into a compartmental linear system. To understand this, one must first grasp Wiechert s conceptual invention of cumomers. Cumomers The word cumomer describes a cumulative isotopomer. It is not a distinct physical entity, like isotopomers, but rather a collection of isotopomers. Cumomers are marked just as isotopomers are labeled. The difference is that with cumomers specific carbons are marked if they have a potential to be labeled in the 13 C sense. For example, look at the hypothetical 5-carbon molecule in figure 4. 6

7 Figure 4: Dark carbons are marked or labeled. The example cumomer in figure 4 is marked in the 2 and 4 positions. This cumomer represents the collection of the 4 isotopomers on the right. A completely unmarked cumomer represents the unlabeled isotopomer. A completely marked cumomer represents all possible isotopomers. Thinking of the labeled system in terms of cumomers reduces the complexity of the mathematics at very little cost. Values for isotopomer fractions and cumomer fractions are interchangeable through a simple linear transformation. This allows for the analytical solution of all intracellular cumomer fractions in terms of input cumomer fractions and flux values. This is possible because the cumomer balance equations are a compartmental linear system. Central to all of this is the notion of cumomer weight, which is assigned based on the number of marked carbons. For example, the cumomer in figure 4 has a weight of 2, because it has 2 marked carbons. Cumomer balance systems are compartmentalized in terms of weight. For each weight there exists a square linear system. Say we have a system with maximum weight of 3. The linear systems could be written as follows: F 1 y 1 = r 1 F 2 y 2 = r 2 (y 1 ) F 3 y 3 = r 3 (y 2,y 1 ) where the subscripts denote weight. In this notation y i symbolizes the fraction vector for cumomers of weight i. F i symbolizes the flux matrix for weight i. It is called a flux matrix because its elements are expressions of only flux values. r i symbolizes the right hand side vector. An important property of the r i vector is that its elements are expressions dependent on flux values and lower weight cumomer fractions. A cumomer of a given weight will never depend on a cumomer of higher weight. In the cumomer balance equations above this is shown by r i s dependence on lower weight cumomers. Solving the systems from low to high weight leads to the complete solution for the cumomer fractions: 7

8 y 1 = F "1 1 r 1 y 2 = F "1 2 r 2 (y 1 ). y 3 = F "1 3 r 3 (y 2,y 1 ) [6] The cumomer balance equations can also be written in another form: Cv. [7] This is closer to the stoichiometric balance equations. C is called a cumomer matrix, because its elements are expressions involving cumomer fractions. v is the vector of fluxes. Each row of C corresponds to the balance equation for a single isotopomer. This form is useful in determining if labeling experiments will help to fully uncover the intracellular fluxes. Information Analysis This section will describe how to use the stoichiometric and cumomer balance systems to determine whether labeling experiments will completely uncover the intracellular fluxes. This process is fairly straightforward. First, we must obtain the flux information by substituting the constrained flux equation [eq. 5] into the cumomer balance [eq. 7]: Cv C(Nx + v 0 ) CNx = "Cv 0 Mx = "m Here, M is called the flux information matrix and x is the vector of free fluxes gained from stoichiometry. M " R m#n, where m (number of matrix rows) is the number of intracellular cumomers and n (number of matrix columns) is the number of free fluxes. For realistic networks, m > n so the flux information system is over determined. Next, we must remove redundancies from the elements of M. This is done by replacing all of the cumomer fraction values in M with their free representation [eq. 6] obtained by solving the linear compartmental system. Finally, Gaussian elimination is used to analyze the rank of M. If the rank of M equals the dimension of x, the system is completely solvable. Even if M is rank deficient (less than the dimension of x) it is able to further constrain the flux values. The metabolic engineer should use this information to decide whether pursuit of carbon labeling data is worthwhile. 8

9 Matlab Tool A set of Matlab tools has been developed to carry out the mathematics described above. Documentation and examples can be found online (Code is available on Dr. Steve Cox s Metabolic Engineering website). The user must specify a metabolic network with carbon transitions and feed it into the code. The tool creates all of the necessary data structures and carries out symbolic full information analysis. At the end, it notifies the user of the rank of the information matrix in comparison to the dimension of the free flux vector. Conclusions This project has completed two goals: to absorb the information spread throughout Wiechert s papers and implement the ideas in Matlab. This paper has covered the basics of mathematics behind metabolic flux analysis using carbon labeling experiments. Future work should be done to make the code handle measurement and statistical concerns. References 1. Nicole Isermann and Wolfgang Wiechert. Metabolic isotopomer labeling systems. Part II: structural flux identifiability analysis. Mathematical Biosciences, Volume 183, Issue 2, June 2003, Pages ( 2. Wolfgang Wiechert and Michael Wurzel. Metabolic isotopomer labeling systems: Part I: global dynamic behavior. Mathematical Biosciences, Volume 169, Issue 2, February 2001, Pages ( 4/2/9d6dbbd409d4e2ebf821021e3e4dbb48) 3. Wolfgang Wiechert et al. Bidirectional Reaction Steps in Metabolic Networks: III. Explicit Solution and Analysis of Isotopomer Labeling Systems. Biotechnology and Bioengineering, Volume 66, Number 2, July 1999, Pages Gregory N. Stephanopoulos, Aristos A. Aristidou, and Jens Nielsen. Metabolic Engineering: Principles and Methodologies. San Diego: Academic Press, Thomas Szyperski. 13 C-NMR, MS and metabolic flux balancing in biotechnology research. Quarterly Reviews of Biophysics, 31, 1 (1998), Pages

Systems II. Metabolic Cycles

Systems II. Metabolic Cycles Systems II Metabolic Cycles Overview Metabolism is central to cellular functioning Maintenance of mass and energy requirements for the cell Breakdown of toxic compounds Consists of a number of major pathways

More information

Kleene s Theorem and the Solution of Metabolic Carbon Labeling Systems

Kleene s Theorem and the Solution of Metabolic Carbon Labeling Systems Kleene s Theorem and the Solution of Metabolic Carbon Labeling Systems Nicole Isermann Michael Weitzel Wolfgang Wiechert Department of Simulation, University of Siegen Paul-Bonatz-Str. 9-11, D-57068 Siegen,

More information

Metabolic Network Analysis

Metabolic Network Analysis Metabolic Network Analysis Christoph Flamm Institute for Theoretical Chemistry University Vienna Harvest Network Course, Brno, April 1-6, 2011 The Systems Biology Modelling Space Network Models Level of

More information

New Results on Energy Balance Analysis of Metabolic Networks

New Results on Energy Balance Analysis of Metabolic Networks New Results on Energy Balance Analysis of Metabolic Networks Qinghua Zhou 1, Simon C.K. Shiu 2,SankarK.Pal 3,andYanLi 1 1 College of Mathematics and Computer Science, Hebei University, Baoding City 071002,

More information

Equivalence of metabolite fragments and flow analysis of isotopomer distributions for flux estimation

Equivalence of metabolite fragments and flow analysis of isotopomer distributions for flux estimation Equivalence of metabolite fragments and flow analysis of isotopomer distributions for flux estimation Ari Rantanen, Hannu Maaheimo 1, Esa Pitkänen, Juho Rousu, Esko Ukkonen, Firstname.Lastname@cs.Helsinki.Fi,

More information

Module 6 : Reaction Kinetics and Dynamics Lecture 28 : Elementary Reactions and Reaction Mechanisms

Module 6 : Reaction Kinetics and Dynamics Lecture 28 : Elementary Reactions and Reaction Mechanisms Module 6 : Reaction Kinetics and Dynamics Lecture 28 : Elementary Reactions and Reaction Mechanisms Objectives In this Lecture you will learn to do the following Define what is an elementary reaction.

More information

Regulation of metabolism

Regulation of metabolism Regulation of metabolism So far in this course we have assumed that the metabolic system is in steady state For the rest of the course, we will abandon this assumption, and look at techniques for analyzing

More information

Basic modeling approaches for biological systems. Mahesh Bule

Basic modeling approaches for biological systems. Mahesh Bule Basic modeling approaches for biological systems Mahesh Bule The hierarchy of life from atoms to living organisms Modeling biological processes often requires accounting for action and feedback involving

More information

General Chemistry Notes Name

General Chemistry Notes Name Bio Honors General Chemistry Notes Name Directions: Carefully read the following information. Look for the ** directions in italics** for prompts where you can do some work. Use the information you have

More information

SPA for quantitative analysis: Lecture 6 Modelling Biological Processes

SPA for quantitative analysis: Lecture 6 Modelling Biological Processes 1/ 223 SPA for quantitative analysis: Lecture 6 Modelling Biological Processes Jane Hillston LFCS, School of Informatics The University of Edinburgh Scotland 7th March 2013 Outline 2/ 223 1 Introduction

More information

Theoretical aspects of C13 metabolic flux analysis with sole quantification of carbon dioxide labeling. Guangquan Shi 04/28/06

Theoretical aspects of C13 metabolic flux analysis with sole quantification of carbon dioxide labeling. Guangquan Shi 04/28/06 Theoretical aspects of C13 metabolic flux analysis with sole quantification of carbon dioxide labeling Guangquan Shi 04/28/06 Omes? One June 26, 2000 President Clinton, with J. Craig Venter, left, and

More information

Overview of Kinetics

Overview of Kinetics Overview of Kinetics [P] t = ν = k[s] Velocity of reaction Conc. of reactant(s) Rate of reaction M/sec Rate constant sec -1, M -1 sec -1 1 st order reaction-rate depends on concentration of one reactant

More information

Unit 1. Unit 2. Unit 6. Unit 3. Unit 5. Unit 4. Solutions

Unit 1. Unit 2. Unit 6. Unit 3. Unit 5. Unit 4. Solutions Chemistry Chemistry Standards The Cobb Teaching and Learning Standards (CT & LS) for science are designed to provide foundational knowledge and skills for all students to develop proficiency in science.

More information

Stoichiometric network analysis

Stoichiometric network analysis Stoichiometric network analysis In stoichiometric analysis of metabolic networks, one concerns the effect of the network structure on the behaviour and capabilities of metabolism. Questions that can be

More information

The Chemistry department approved by the American Chemical Society offers a Chemistry degree in the following concentrations:

The Chemistry department approved by the American Chemical Society offers a Chemistry degree in the following concentrations: Chemistry 1 Chemistry 203-C Materials Science Building Telephone: 256.824.6153 Email: chem.admin@uah.edu (chem@uah.edu) The Chemistry department approved by the American Chemical Society offers a Chemistry

More information

Chapter 6- An Introduction to Metabolism*

Chapter 6- An Introduction to Metabolism* Chapter 6- An Introduction to Metabolism* *Lecture notes are to be used as a study guide only and do not represent the comprehensive information you will need to know for the exams. The Energy of Life

More information

Unit 3. 4 weeks BL/8 weeks YR

Unit 3. 4 weeks BL/8 weeks YR Chemistry Chemistry Standards The Cobb Teaching and Learning Standards (CT & LS) for science are designed to provide foundational knowledge and skills for all students to develop proficiency in science.

More information

V15 Flux Balance Analysis Extreme Pathways

V15 Flux Balance Analysis Extreme Pathways V15 Flux Balance Analysis Extreme Pathways Stoichiometric matrix S: m n matrix with stochiometries of the n reactions as columns and participations of m metabolites as rows. The stochiometric matrix is

More information

Introduction to Bioinformatics

Introduction to Bioinformatics Systems biology Introduction to Bioinformatics Systems biology: modeling biological p Study of whole biological systems p Wholeness : Organization of dynamic interactions Different behaviour of the individual

More information

A Correlation of. to the. Georgia Standards of Excellence Chemistry

A Correlation of. to the. Georgia Standards of Excellence Chemistry A Correlation of to the A Correlation of Pearson Introduction The following document demonstrates how Pearson supports the Georgia Standards for Excellence in. Correlation references are to the Student

More information

V19 Metabolic Networks - Overview

V19 Metabolic Networks - Overview V19 Metabolic Networks - Overview There exist different levels of computational methods for describing metabolic networks: - stoichiometry/kinetics of classical biochemical pathways (glycolysis, TCA cycle,...

More information

Lab 2 Worksheet. Problems. Problem 1: Geometry and Linear Equations

Lab 2 Worksheet. Problems. Problem 1: Geometry and Linear Equations Lab 2 Worksheet Problems Problem : Geometry and Linear Equations Linear algebra is, first and foremost, the study of systems of linear equations. You are going to encounter linear systems frequently in

More information

Grundlagen der Systembiologie und der Modellierung epigenetischer Prozesse

Grundlagen der Systembiologie und der Modellierung epigenetischer Prozesse Grundlagen der Systembiologie und der Modellierung epigenetischer Prozesse Sonja J. Prohaska Bioinformatics Group Institute of Computer Science University of Leipzig October 25, 2010 Genome-scale in silico

More information

Stoichiometric analysis of metabolic networks

Stoichiometric analysis of metabolic networks 1 SGN-6156 Computational systems biology II Antti Larjo antti.larjo@tut.fi Computational Systems Biology group 30.4.2008 30.4.2008 2 Contents Networks in cells Metabolism Metabolic networks and pathways

More information

Program for the rest of the course

Program for the rest of the course Program for the rest of the course 16.4 Enzyme kinetics 17.4 Metabolic Control Analysis 19.4. Exercise session 5 23.4. Metabolic Control Analysis, cont. 24.4 Recap 27.4 Exercise session 6 etabolic Modelling

More information

MATH 3511 Lecture 1. Solving Linear Systems 1

MATH 3511 Lecture 1. Solving Linear Systems 1 MATH 3511 Lecture 1 Solving Linear Systems 1 Dmitriy Leykekhman Spring 2012 Goals Review of basic linear algebra Solution of simple linear systems Gaussian elimination D Leykekhman - MATH 3511 Introduction

More information

International Journal of Scientific & Engineering Research, Volume 6, Issue 2, February ISSN

International Journal of Scientific & Engineering Research, Volume 6, Issue 2, February ISSN International Journal of Scientific & Engineering Research, Volume 6, Issue 2, February-2015 273 Analogizing And Investigating Some Applications of Metabolic Pathway Analysis Methods Gourav Mukherjee 1,

More information

Computer Simulation and Data Analysis in Molecular Biology and Biophysics

Computer Simulation and Data Analysis in Molecular Biology and Biophysics Victor Bloomfield Computer Simulation and Data Analysis in Molecular Biology and Biophysics An Introduction Using R Springer Contents Part I The Basics of R 1 Calculating with R 3 1.1 Installing R 3 1.1.1

More information

Designing Information Devices and Systems I Fall 2018 Lecture Notes Note Introduction to Linear Algebra the EECS Way

Designing Information Devices and Systems I Fall 2018 Lecture Notes Note Introduction to Linear Algebra the EECS Way EECS 16A Designing Information Devices and Systems I Fall 018 Lecture Notes Note 1 1.1 Introduction to Linear Algebra the EECS Way In this note, we will teach the basics of linear algebra and relate it

More information

Crosswalk of Georgia Performance Standards & Georgia Standards of Excellence GSE Implementation in Chemistry

Crosswalk of Georgia Performance Standards & Georgia Standards of Excellence GSE Implementation in Chemistry Implementation in 2017-18 SC1. Students will analyze the nature of matter and its classifications. a. Relate the role of nuclear fusion in producing essentially all elements heavier than helium. b. Identify

More information

Flux Balance Analysis

Flux Balance Analysis Lecture 10: Flux Balance Analysis Tue 7 March 2006 with the collaboration of Luna De Ferrari 1 Images and text from: E. Klipp, Systems Biology in Practice, Wiley-VCH, 2005 Ch 5 Edwards JS, Palsson BO,

More information

Objectives INTRODUCTION TO METABOLISM. Metabolism. Catabolic Pathways. Anabolic Pathways 3/6/2011. How to Read a Chemical Equation

Objectives INTRODUCTION TO METABOLISM. Metabolism. Catabolic Pathways. Anabolic Pathways 3/6/2011. How to Read a Chemical Equation Objectives INTRODUCTION TO METABOLISM. Chapter 8 Metabolism, Energy, and Life Explain the role of catabolic and anabolic pathways in cell metabolism Distinguish between kinetic and potential energy Distinguish

More information

Designing Information Devices and Systems I Spring 2018 Lecture Notes Note Introduction to Linear Algebra the EECS Way

Designing Information Devices and Systems I Spring 2018 Lecture Notes Note Introduction to Linear Algebra the EECS Way EECS 16A Designing Information Devices and Systems I Spring 018 Lecture Notes Note 1 1.1 Introduction to Linear Algebra the EECS Way In this note, we will teach the basics of linear algebra and relate

More information

MTH Linear Algebra. Study Guide. Dr. Tony Yee Department of Mathematics and Information Technology The Hong Kong Institute of Education

MTH Linear Algebra. Study Guide. Dr. Tony Yee Department of Mathematics and Information Technology The Hong Kong Institute of Education MTH 3 Linear Algebra Study Guide Dr. Tony Yee Department of Mathematics and Information Technology The Hong Kong Institute of Education June 3, ii Contents Table of Contents iii Matrix Algebra. Real Life

More information

Linear Algebra. The analysis of many models in the social sciences reduces to the study of systems of equations.

Linear Algebra. The analysis of many models in the social sciences reduces to the study of systems of equations. POLI 7 - Mathematical and Statistical Foundations Prof S Saiegh Fall Lecture Notes - Class 4 October 4, Linear Algebra The analysis of many models in the social sciences reduces to the study of systems

More information

Cellular Metabolic Models

Cellular Metabolic Models Cellular Metabolic Models. Cellular metabolism. Modeling cellular metabolism. Flux balance model of yeast glycolysis 4. Kinetic model of yeast glycolysis Cellular Metabolic Models Cellular Metabolism Basic

More information

A. One-Substrate Reactions (1) Kinetic concepts

A. One-Substrate Reactions (1) Kinetic concepts A. One-Substrate Reactions (1) Kinetic concepts (2) Kinetic analysis (a) Briggs-Haldane steady-state treatment (b) Michaelis constant (K m ) (c) Specificity constant (3) Graphical analysis (4) Practical

More information

Exam I Solutions Chem 6, 9 Section, Spring 2002

Exam I Solutions Chem 6, 9 Section, Spring 2002 1. (a) Two researchers at the University of Nebraska recently published a paper on the rate of the disappearance of World Wide Web links, a phenomenon called link rot. They asked the question, If I place

More information

Predicting rice (Oryza sativa) metabolism

Predicting rice (Oryza sativa) metabolism Predicting rice (Oryza sativa) metabolism Sudip Kundu Department of Biophysics, Molecular Biology & Bioinformatics, University of Calcutta, WB, India. skbmbg@caluniv.ac.in Collaborators: Mark G Poolman

More information

Field 045: Science Life Science Assessment Blueprint

Field 045: Science Life Science Assessment Blueprint Field 045: Science Life Science Assessment Blueprint Domain I Foundations of Science 0001 The Nature and Processes of Science (Standard 1) 0002 Central Concepts and Connections in Science (Standard 2)

More information

V14 extreme pathways

V14 extreme pathways V14 extreme pathways A torch is directed at an open door and shines into a dark room... What area is lighted? Instead of marking all lighted points individually, it would be sufficient to characterize

More information

CHEM Chemical Kinetics. Reaction Mechanisms

CHEM Chemical Kinetics. Reaction Mechanisms Chemical Kinetics Deri ed Rate La s from Derived Rate Laws from Reaction Mechanisms Reaction Mechanism Determine the rate law by experiment Devise a reaction mechanism If the predicted and experimental

More information

Stoichiometric vector and matrix

Stoichiometric vector and matrix Stoichiometric vector and matrix The stoichiometric coefficients of a reaction are collected to a vector s r In s r there is a one position for each metabolite in the metabolic system, and the stoichiometric

More information

Metabolic modelling. Metabolic networks, reconstruction and analysis. Esa Pitkänen Computational Methods for Systems Biology 1 December 2009

Metabolic modelling. Metabolic networks, reconstruction and analysis. Esa Pitkänen Computational Methods for Systems Biology 1 December 2009 Metabolic modelling Metabolic networks, reconstruction and analysis Esa Pitkänen Computational Methods for Systems Biology 1 December 2009 Department of Computer Science, University of Helsinki Metabolic

More information

Introduction to Systems Biology. James Gomes KSBS, IITD

Introduction to Systems Biology. James Gomes KSBS, IITD Introduction to Systems Biology James Gomes KSBS, IITD Inner Life of a Cell Cells Exhibit a High Degree of Complexity What are the different types of cells? Is Life an emergent property? What is Systems

More information

Chapter 7: Metabolic Networks

Chapter 7: Metabolic Networks Chapter 7: Metabolic Networks 7.1 Introduction Prof. Yechiam Yemini (YY) Computer Science epartment Columbia University Introduction Metabolic flux analysis Applications Overview 2 1 Introduction 3 Metabolism:

More information

Lecture 6 Asymptotic Structure for Four-Step Premixed Stoichiometric Methane Flames

Lecture 6 Asymptotic Structure for Four-Step Premixed Stoichiometric Methane Flames Lecture 6 Asymptotic Structure for Four-Step Premixed Stoichiometric Methane Flames 6.-1 Previous lecture: Asymptotic description of premixed flames based on an assumed one-step reaction. basic understanding

More information

Computational Systems Biology Exam

Computational Systems Biology Exam Computational Systems Biology Exam Dr. Jürgen Pahle Aleksandr Andreychenko, M.Sc. 31 July, 2012 Name Matriculation Number Do not open this exam booklet before we ask you to. Do read this page carefully.

More information

Advanced Chemical Reaction Engineering Prof. H. S. Shankar Department of Chemical Engineering IIT Bombay. Lecture - 03 Design Equations-1

Advanced Chemical Reaction Engineering Prof. H. S. Shankar Department of Chemical Engineering IIT Bombay. Lecture - 03 Design Equations-1 (Refer Slide Time: 00:19) Advanced Chemical Reaction Engineering Prof. H. S. Shankar Department of Chemical Engineering IIT Bombay Lecture - 03 Design Equations-1 We are looking at advanced reaction engineering;

More information

Chemistry. Faculty Kent Davis, chair; Pablo Hilario, Denise Lee-Haye, Marie Pak, Robert Wilson Departmental Office: 355 Chan Shun Hall; (707)

Chemistry. Faculty Kent Davis, chair; Pablo Hilario, Denise Lee-Haye, Marie Pak, Robert Wilson Departmental Office: 355 Chan Shun Hall; (707) Faculty Kent Davis, chair; Pablo Hilario, Denise Lee-Haye, Marie Pak, Robert Wilson Departmental Office: 355 Chan Shun Hall; (707) 965-7597 Degrees and Programs Chemistry, B.S....77 Chemistry, B.A....78

More information

Section 1.1: Systems of Linear Equations. A linear equation: a 1 x 1 a 2 x 2 a n x n b. EXAMPLE: 4x 1 5x 2 2 x 1 and x x 1 x 3

Section 1.1: Systems of Linear Equations. A linear equation: a 1 x 1 a 2 x 2 a n x n b. EXAMPLE: 4x 1 5x 2 2 x 1 and x x 1 x 3 Section 1.1: Systems of Linear Equations A linear equation: a 1 x 1 a 2 x 2 a n x n b EXAMPLE: 4x 1 5x 2 2 x 1 and x 2 2 6 x 1 x 3 rearranged rearranged 3x 1 5x 2 2 2x 1 x 2 x 3 2 6 Not linear: 4x 1 6x

More information

Lecture 18: The Rank of a Matrix and Consistency of Linear Systems

Lecture 18: The Rank of a Matrix and Consistency of Linear Systems Lecture 18: The Rank of a Matrix and Consistency of Linear Systems Winfried Just Department of Mathematics, Ohio University February 28, 218 Review: The linear span Definition Let { v 1, v 2,..., v n }

More information

Getting Started with Communications Engineering

Getting Started with Communications Engineering 1 Linear algebra is the algebra of linear equations: the term linear being used in the same sense as in linear functions, such as: which is the equation of a straight line. y ax c (0.1) Of course, if we

More information

Analyze the roles of enzymes in biochemical reactions

Analyze the roles of enzymes in biochemical reactions ENZYMES and METABOLISM Elements: Cell Biology (Enzymes) Estimated Time: 6 7 hours By the end of this course, students will have an understanding of the role of enzymes in biochemical reactions. Vocabulary

More information

Chemical Kinetics. Goal. Objectives

Chemical Kinetics. Goal. Objectives 4 Chemical Kinetics To understand the physical Goal factors that determine reaction rates. bjectives After this chapter, you should be able to: describe the factors that determine reaction rates. identify

More information

Kinetic isotope effects significantly influence intracellular metabolite [superscript 13]C labeling patterns and flux determination

Kinetic isotope effects significantly influence intracellular metabolite [superscript 13]C labeling patterns and flux determination Kinetic isotope effects significantly influence intracellular metabolite [superscript 13]C labeling patterns and flux determination The MIT Faculty has made this article openly available. Please share

More information

MTH 2032 Semester II

MTH 2032 Semester II MTH 232 Semester II 2-2 Linear Algebra Reference Notes Dr. Tony Yee Department of Mathematics and Information Technology The Hong Kong Institute of Education December 28, 2 ii Contents Table of Contents

More information

Updated: 10/11/2018 Page 1 of 5

Updated: 10/11/2018 Page 1 of 5 A. Academic Division: Health Sciences B. Discipline: Biology C. Course Number and Title: BIOL1230 Biology I MASTER SYLLABUS 2018-2019 D. Course Coordinator: Justin Tickhill Assistant Dean: Melinda Roepke,

More information

Simplicity is Complexity in Masquerade. Michael A. Savageau The University of California, Davis July 2004

Simplicity is Complexity in Masquerade. Michael A. Savageau The University of California, Davis July 2004 Simplicity is Complexity in Masquerade Michael A. Savageau The University of California, Davis July 2004 Complexity is Not Simplicity in Masquerade -- E. Yates Simplicity is Complexity in Masquerade One

More information

Time-Dependent Statistical Mechanics 5. The classical atomic fluid, classical mechanics, and classical equilibrium statistical mechanics

Time-Dependent Statistical Mechanics 5. The classical atomic fluid, classical mechanics, and classical equilibrium statistical mechanics Time-Dependent Statistical Mechanics 5. The classical atomic fluid, classical mechanics, and classical equilibrium statistical mechanics c Hans C. Andersen October 1, 2009 While we know that in principle

More information

Introduction on metabolism & refresher in enzymology

Introduction on metabolism & refresher in enzymology Introduction on metabolism & refresher in enzymology Daniel Kahn Laboratoire de Biométrie & Biologie Evolutive Lyon 1 University & INRA MIA Department Daniel.Kahn@univ-lyon1.fr General objectives of the

More information

Biological Pathways Representation by Petri Nets and extension

Biological Pathways Representation by Petri Nets and extension Biological Pathways Representation by and extensions December 6, 2006 Biological Pathways Representation by and extension 1 The cell Pathways 2 Definitions 3 4 Biological Pathways Representation by and

More information

Bioinformatics: Network Analysis

Bioinformatics: Network Analysis Bioinformatics: Network Analysis Reaction Kinetics COMP 572 (BIOS 572 / BIOE 564) - Fall 2013 Luay Nakhleh, Rice University 1 Reaction kinetics is the study of how fast chemical reactions take place, what

More information

METABOLIC PATHWAY PREDICTION/ALIGNMENT

METABOLIC PATHWAY PREDICTION/ALIGNMENT COMPUTATIONAL SYSTEMIC BIOLOGY METABOLIC PATHWAY PREDICTION/ALIGNMENT Hofestaedt R*, Chen M Bioinformatics / Medical Informatics, Technische Fakultaet, Universitaet Bielefeld Postfach 10 01 31, D-33501

More information

Unit 7 Part I: Introductions to Biochemistry

Unit 7 Part I: Introductions to Biochemistry Unit 7 Part I: Introductions to Biochemistry Chemical Reactions, Enzymes and ATP 19-Mar-14 Averett 1 Chemical Reactions Chemical Reactions Process by which one set of chemicals is changed into another

More information

CORRELATIONS FOR THE GEORGIA SCIENCE STANDARDS OF EXCELLENCE 2016 GRADES CHEMISTRY

CORRELATIONS FOR THE GEORGIA SCIENCE STANDARDS OF EXCELLENCE 2016 GRADES CHEMISTRY ! CORRELATIONS FOR THE GEORGIA SCIENCE STANDARDS OF EXCELLENCE 2016 GRADES 10-12 CHEMISTRY A Natural Approach to Chemistry (NAC) is written by Hsu, Chaniotakis, Carlisle, and Damelin, and is published

More information

Satisfying Flux Balance and Mass-Action Kinetics in a Network of Biochemical Reactions

Satisfying Flux Balance and Mass-Action Kinetics in a Network of Biochemical Reactions Satisfying Flux Balance and Mass-Action Kinetics in a Network of Biochemical Reactions Speakers Senior Personnel Collaborators Michael Saunders and Ines Thiele Ronan Fleming, Bernhard Palsson, Yinyu Ye

More information

Math Camp II. Basic Linear Algebra. Yiqing Xu. Aug 26, 2014 MIT

Math Camp II. Basic Linear Algebra. Yiqing Xu. Aug 26, 2014 MIT Math Camp II Basic Linear Algebra Yiqing Xu MIT Aug 26, 2014 1 Solving Systems of Linear Equations 2 Vectors and Vector Spaces 3 Matrices 4 Least Squares Systems of Linear Equations Definition A linear

More information

How fast reactants turn into products. Usually measured in Molarity per second units. Kinetics

How fast reactants turn into products. Usually measured in Molarity per second units. Kinetics How fast reactants turn into products. Usually measured in Molarity per second units. Kinetics Reaction rated are fractions of a second for fireworks to explode. Reaction Rates takes years for a metal

More information

Integrated Knowledge-based Reverse Engineering of Metabolic Pathways

Integrated Knowledge-based Reverse Engineering of Metabolic Pathways Integrated Knowledge-based Reverse Engineering of Metabolic Pathways Shuo-Huan Hsu, Priyan R. Patkar, Santhoi Katare, John A. Morgan and Venkat Venkatasubramanian School of Chemical Engineering, Purdue

More information

Canadian Advanced Senior High

Canadian Advanced Senior High Canadian Advanced Senior High Department: Science Course Development Date: November 2017 Course Title: Biology Grade: 12 Course Type: Ministry Course Code: University SBI4U Credit Value: 1 Hours: 110 Ministry

More information

Chemical Kinetics I: The Dry Lab. Up until this point in our study of physical chemistry we have been interested in

Chemical Kinetics I: The Dry Lab. Up until this point in our study of physical chemistry we have been interested in Chemical Kinetics I: The Dry Lab Up until this point in our study of physical chemistry we have been interested in equilibrium properties; now we will begin to investigate non-equilibrium properties and

More information

FUNDAMENTALS of SYSTEMS BIOLOGY From Synthetic Circuits to Whole-cell Models

FUNDAMENTALS of SYSTEMS BIOLOGY From Synthetic Circuits to Whole-cell Models FUNDAMENTALS of SYSTEMS BIOLOGY From Synthetic Circuits to Whole-cell Models Markus W. Covert Stanford University 0 CRC Press Taylor & Francis Group Boca Raton London New York Contents /... Preface, xi

More information

Chapter 8: An Introduction to Metabolism

Chapter 8: An Introduction to Metabolism AP Biology Reading Guide Name Chapter 8: An Introduction to Metabolism Concept 8.1 An organism s metabolism transforms matter and energy, subject to the laws of thermodynamics 1. Define metabolism. 2.

More information

CHAPTER 2 THE CHEMICAL BASIS OF LIFE

CHAPTER 2 THE CHEMICAL BASIS OF LIFE CHAPTER 2 THE CHEMICAL BASIS OF LIFE CHAPTER OVERVIEW This chapter introduces very basic concepts of chemistry, emphasizing the structure of atoms and how they combine (form bonds). The types of bonds,

More information

ANALYSIS OF BIOLOGICAL NETWORKS USING HYBRID SYSTEMS THEORY. Nael H. El-Farra, Adiwinata Gani & Panagiotis D. Christofides

ANALYSIS OF BIOLOGICAL NETWORKS USING HYBRID SYSTEMS THEORY. Nael H. El-Farra, Adiwinata Gani & Panagiotis D. Christofides ANALYSIS OF BIOLOGICAL NETWORKS USING HYBRID SYSTEMS THEORY Nael H El-Farra, Adiwinata Gani & Panagiotis D Christofides Department of Chemical Engineering University of California, Los Angeles 2003 AIChE

More information

Some Notes on Linear Algebra

Some Notes on Linear Algebra Some Notes on Linear Algebra prepared for a first course in differential equations Thomas L Scofield Department of Mathematics and Statistics Calvin College 1998 1 The purpose of these notes is to present

More information

ENZYME KINETICS AND INHIBITION

ENZYME KINETICS AND INHIBITION ENZYME KINETICS AND INHIBITION The kinetics of reactions involving enzymes are a little bit different from other reactions. First of all, there are sometimes lots of steps involved. Also, the reaction

More information

SECTION II: KINETICS AND BIOREACTOR DESIGN: JAVIER CALZADA FUNES

SECTION II: KINETICS AND BIOREACTOR DESIGN: JAVIER CALZADA FUNES SECTION II: KINETICS AND BIOREACTOR DESIGN: LESSON 9.1. - Enzymatic kinetics, microbial kinetics and metabolic stoichiometry - Brief review on enzymatic reaction kinetics JAVIER CALZADA FUNES Biotechnology

More information

UNIT 1: CHEMISTRY FOUNDATIONS

UNIT 1: CHEMISTRY FOUNDATIONS Advanced Placement AP Chemistry builds students' understanding of the nature and reactivity of matter. After studying chemical reactions and electrochemistry, students move on to understand how the chemical

More information

Prentice Hall. Chemistry, (Wilbraham) 2008, National SE, Georgia TE. Grades 9-12

Prentice Hall. Chemistry, (Wilbraham) 2008, National SE, Georgia TE. Grades 9-12 Prentice Hall Chemistry, (Wilbraham) 2008, National SE, Georgia TE Grades 9-12 C O R R E L A T E D T O Georgia Performance s for Chemistry (Course Number 40.05100 Chemistry I) Grades 9-12 FORMAT FOR CORRELATION

More information

BASIC NOTIONS. x + y = 1 3, 3x 5y + z = A + 3B,C + 2D, DC are not defined. A + C =

BASIC NOTIONS. x + y = 1 3, 3x 5y + z = A + 3B,C + 2D, DC are not defined. A + C = CHAPTER I BASIC NOTIONS (a) 8666 and 8833 (b) a =6,a =4 will work in the first case, but there are no possible such weightings to produce the second case, since Student and Student 3 have to end up with

More information

Georgia Standards of Excellence Curriculum Map High School Chemistry

Georgia Standards of Excellence Curriculum Map High School Chemistry Vocabulary Standards and elements with main topics Unit Name Focus Georgia Standards of Excellence Curriculum Map 2017-18 High School Chemistry Unit # Unit 1 Unit 2 Unit 3 Unit 4 Pace 26 days 12 days 30

More information

Chapter 8: Photosynthesis. Name Block

Chapter 8: Photosynthesis. Name Block Fred and Theresa Holtzclaw Updated by Chris Chou for Campbell Biology in Focus, 2nd Ed. (Oct. 2017) Name Block This chapter is as challenging as the one you just finished on cellular respiration. However,

More information

Chapter 8: An Introduction to Metabolism

Chapter 8: An Introduction to Metabolism Name Period Concept 8.1 An organism s metabolism transforms matter and energy, subject to the laws of thermodynamics 1. Define metabolism. 2. There are two types of reactions in metabolic pathways: anabolic

More information

6. MESH ANALYSIS 6.1 INTRODUCTION

6. MESH ANALYSIS 6.1 INTRODUCTION 6. MESH ANALYSIS INTRODUCTION PASSIVE SIGN CONVENTION PLANAR CIRCUITS FORMATION OF MESHES ANALYSIS OF A SIMPLE CIRCUIT DETERMINANT OF A MATRIX CRAMER S RULE GAUSSIAN ELIMINATION METHOD EXAMPLES FOR MESH

More information

Mathematical Biology - Lecture 1 - general formulation

Mathematical Biology - Lecture 1 - general formulation Mathematical Biology - Lecture 1 - general formulation course description Learning Outcomes This course is aimed to be accessible both to masters students of biology who have a good understanding of the

More information

BMC Systems Biology. Open Access. Abstract

BMC Systems Biology. Open Access. Abstract BMC Systems Biology BioMed Central Research article Hybrid optimization for 13 C metabolic flux analysis using systems parametrized by compactification Tae Hoon Yang* 1,, Oliver Frick and Elmar Heinzle

More information

Concept 10.1 Photosynthesis converts light energy to the chemical energy of food

Concept 10.1 Photosynthesis converts light energy to the chemical energy of food Name Period Chapter 10: Photosynthesis This chapter is as challenging as the one you just finished on cellular respiration. However, conceptually it will be a little easier because the concepts learned

More information

General Chemistry 201 Section ABC Harry S. Truman College Spring Semester 2014

General Chemistry 201 Section ABC Harry S. Truman College Spring Semester 2014 Instructor: Michael Davis Office: 3226 Phone: 773 907 4718 Office Hours: Tues 9:00 12:00 Wed 1:00 3:00 Thurs 9:00 12:00 Email: mdavis@ccc.edu Website: http://faradaysclub.com http://ccc.blackboard.com

More information

AP BIOLOGY: READING ASSIGNMENT FOR CHAPTER 2. Particle Charge Mass Location

AP BIOLOGY: READING ASSIGNMENT FOR CHAPTER 2. Particle Charge Mass Location 1) Fill in the names beside the symbols of the following elements commonly found in living matter: a. Ca b. P c. K d. S e. Na f. Cl g. Mg 2) The different between the mass number and the atomic number

More information

Chapter 8 Educational Goals

Chapter 8 Educational Goals Chapter 8 Educational Goals 1. Given a chemical equation, write the law of mass action. 2. Given the equilibrium constant (K eq ) for a reaction, predict whether the reactants or products are predominant.

More information

Lecture: Computational Systems Biology Universität des Saarlandes, SS Introduction. Dr. Jürgen Pahle

Lecture: Computational Systems Biology Universität des Saarlandes, SS Introduction. Dr. Jürgen Pahle Lecture: Computational Systems Biology Universität des Saarlandes, SS 2012 01 Introduction Dr. Jürgen Pahle 24.4.2012 Who am I? Dr. Jürgen Pahle 2009-2012 Manchester Interdisciplinary Biocentre, The University

More information

Lassen Community College Course Outline

Lassen Community College Course Outline Lassen Community College Course Outline BIOL-1 Principles of Molecular and Cellular Biology 4.0 Units I. Catalog Description A course in principles of biology, with special emphasis given to molecular

More information

Sample Questions for the Chemistry of Life Topic Test

Sample Questions for the Chemistry of Life Topic Test Sample Questions for the Chemistry of Life Topic Test 1. Enzymes play a crucial role in biology by serving as biological catalysts, increasing the rates of biochemical reactions by decreasing their activation

More information

CHEMISTRY (CHEM) CHEM 5800 Principles Of Materials Chemistry. Tutorial in selected topics in materials chemistry. S/U grading only.

CHEMISTRY (CHEM) CHEM 5800 Principles Of Materials Chemistry. Tutorial in selected topics in materials chemistry. S/U grading only. Chemistry (CHEM) 1 CHEMISTRY (CHEM) CHEM 5100 Principles of Organic and Inorganic Chemistry Study of coordination compounds with a focus on ligand bonding, electron counting, molecular orbital theory,

More information

ENZYME SCIENCE AND ENGINEERING PROF. SUBHASH CHAND DEPARTMENT OF BIOCHEMICAL ENGINEERING AND BIOTECHNOLOGY IIT DELHI LECTURE 6

ENZYME SCIENCE AND ENGINEERING PROF. SUBHASH CHAND DEPARTMENT OF BIOCHEMICAL ENGINEERING AND BIOTECHNOLOGY IIT DELHI LECTURE 6 ENZYME SCIENCE AND ENGINEERING PROF. SUBHASH CHAND DEPARTMENT OF BIOCHEMICAL ENGINEERING AND BIOTECHNOLOGY IIT DELHI LECTURE 6 KINETICS OF ENZYME CATALYSED REACTIONS Having understood the chemical and

More information

Biology A: Chapter 5 Annotating Notes

Biology A: Chapter 5 Annotating Notes Name: Pd: Biology A: Chapter 5 Annotating Notes -As you read your textbook, please fill out these notes. -Read each paragraph state the big/main idea on the left side. - On the right side you should take

More information

AP BIOLOGY

AP BIOLOGY AP BIOLOGY 2016 2017 Welcome to Advanced Placement Biology. The following contains statements of the four Big Ideas that College Board has developed. 1. The process of evolution explains the diversity

More information

It s a Small World After All

It s a Small World After All It s a Small World After All Engage: Cities, factories, even your own home is a network of dependent and independent parts that make the whole function properly. Think of another network that has subunits

More information