Difference between revisions of "Team:Harvard BioDesign/Model"

(Removed hyperlink in page title)
 
(17 intermediate revisions by 3 users not shown)
Line 1: Line 1:
{{Harvard_BioDesign}}
+
<!--
 +
Helios by HTML5 UP
 +
html5up.net | @ajlkn
 +
Free for personal and commercial use under the CCA 3.0 license (html5up.net/license)
 +
-->
 
<html>
 
<html>
 +
<head>
 +
<title>Harvard BioDesign 2016</title>
 +
<meta charset="utf-8" />
 +
<meta name="viewport" content="width=device-width, initial-scale=1" />
  
 +
<!--[if lte IE 8]><script src="assets/js/ie/html5shiv.js"></script><![endif]-->
 +
 +
<link rel="stylesheet" type="text/css" href="https://2016.igem.org/Template:Harvard_BioDesign/assets_css_main?action=raw&ctype=text/css" />
 +
<!--[if lte IE 8]><link rel="stylesheet" href="assets/css/ie8.css" /><![endif]-->
  
 +
 +
<!-- Scripts -->
 +
<script type="text/javascript" src="https://2016.igem.org/Template:Harvard_BioDesign/assets_js_jquery_min?
 +
action=raw&ctype=text/javascript"></script>
 +
<script type="text/javascript" src="https://2016.igem.org/Template:Harvard_BioDesign/assets_js_jquery_dropotron_min?
 +
action=raw&ctype=text/javascript"></script>
 +
<script type="text/javascript" src="https://2016.igem.org/Template:Harvard_BioDesign/assets_js_jquery_scrolly_min?
 +
action=raw&ctype=text/javascript"></script>
 +
<script type="text/javascript" src="https://2016.igem.org/Template:Harvard_BioDesign/assets_js_jquery_onvisible_min?
 +
action=raw&ctype=text/javascript"></script>
 +
<script type="text/javascript" src="https://2016.igem.org/Template:Harvard_BioDesign/assets_js_skel_min?
 +
action=raw&ctype=text/javascript"></script>
 +
<script type="text/javascript" src="https://2016.igem.org/Template:Harvard_BioDesign/assets_js_util?
 +
action=raw&ctype=text/javascript"></script>
 +
<script type="text/javascript" src="https://2016.igem.org/Template:Harvard_BioDesign/assets_js_main?
 +
action=raw&ctype=text/javascript"></script>
 +
</head>
 +
<body class="homepage">
 +
<div id="page-wrapper">
  
 +
<!-- Header -->
 +
<div id="header">
  
 +
<!-- Inner -->
 +
<!--
 +
<div class="inner">
 +
<header>
 +
<h1><a href="https://2016.igem.org/Team:Harvard_BioDesign" id="logo">Harvard BioDesign</a></h1>
 +
<hr />
 +
<p>Description</p>
 +
</header>
 +
<footer>
 +
<a href="#banner" class="button circled scrolly">Start</a>
 +
</footer>
 +
</div>
 +
-->
  
  
 +
<!-- Nav -->
 +
<nav id="nav" >
 +
<ul>
 +
<li><a href="https://2016.igem.org/Team:Harvard_BioDesign">Home</a></li>
 +
<li>
 +
<a href="#">Team</a>
 +
<ul>
 +
<li><a href="https://2016.igem.org/Team:Harvard_BioDesign/Team">Team</a></li>
 +
<li><a href="https://2016.igem.org/Team:Harvard_BioDesign/Collaborations">Collaborations</a></li>
 +
</ul>
 +
</li>
 +
<li>
 +
<a href="#">Project</a>
 +
<ul>
 +
<li><a href="https://2016.igem.org/Team:Harvard_BioDesign/Description">Description</a></li>
 +
<li><a href="https://2016.igem.org/Team:Harvard_BioDesign/Design">Design</a></li>
 +
<li><a href="https://2016.igem.org/Team:Harvard_BioDesign/Experiments">Experiments</a></li>
 +
<li><a href="https://2016.igem.org/Team:Harvard_BioDesign/Proof">Proof of Concept</a></li>
 +
<li><a href="https://2016.igem.org/Team:Harvard_BioDesign/Demonstrate">Demonstrate</a></li>
 +
<li><a href="https://2016.igem.org/Team:Harvard_BioDesign/Results">Results</a></li>
 +
<li><a href="https://2016.igem.org/Team:Harvard_BioDesign/Notebook">Notebook</a></li>
 +
</ul>
 +
</li>
 +
<li>
 +
<a href="#">Parts</a>
 +
<ul>
 +
<li><a href="https://2016.igem.org/Team:Harvard_BioDesign/Parts">Parts</a></li>
 +
<li><a href="https://2016.igem.org/Team:Harvard_BioDesign/Basic_Part">Basic Parts</a></li>
 +
<li><a href="https://2016.igem.org/Team:Harvard_BioDesign/Composite_Part">Composite Parts</a></li>
 +
<li><a href="https://2016.igem.org/Team:Harvard_BioDesign/Part_Collection">Part Collection</a></li>
 +
</ul>
 +
</li>
  
 +
<!-- R edits -->
  
<div class="column full_size">
+
<li><a href="https://2016.igem.org/Team:Harvard_BioDesign/Safety">Safety</a></li>
<h2> Modeling</h2>
+
<li><a href="https://2016.igem.org/Team:Harvard_BioDesign/Attributions">Attributions</a></li>
<p>Mathematical models and computer simulations provide a great way to describe the function and operation of BioBrick Parts and Devices. Synthetic Biology is an engineering discipline, and part of engineering is simulation and modeling to determine the behavior of your design before you build it. Designing and simulating can be iterated many times in a computer before moving to the lab. This award is for teams who build a model of their system and use it to inform system design or simulate expected behavior in conjunction with experiments in the wetlab.</p>
+
<li>
 +
<a href="#">Human Practices</a>
 +
<ul>
 +
<li><a href="https://2016.igem.org/Team:Harvard_BioDesign/Human_Practices">Human Practices</a></li>
 +
<li><a href="https://2016.igem.org/Team:Harvard_BioDesign/HP/Silver">Silver</a></li>
 +
<li><a href="https://2016.igem.org/Team:Harvard_BioDesign/HP/Gold">Gold</a></li>
 +
<li><a href="https://2016.igem.org/Team:Harvard_BioDesign/Integrated_Practices">Integrated Practices</a></li>
 +
<li><a href="https://2016.igem.org/Team:Harvard_BioDesign/Engagement">Engagement</a></li>
 +
</ul>
 +
</li>
 +
<li>
 +
<a href="#">Awards</a>
 +
<ul>
 +
<li><a href="https://2016.igem.org/Team:Harvard_BioDesign/Model">Model</a></li>
 +
<li><a href="https://2016.igem.org/Team:Harvard_BioDesign/Achievements">Achievements</a></li>
 +
</ul>
 +
</li>
  
 +
</ul>
 +
</nav>
 +
 +
</div>
 +
 +
 +
<!-- Main -->
 +
<div class="wrapper style1">
 +
 +
<div class="container">
 +
<article id="main" class="special">
 +
<header>
 +
<h2>Model</h2>
 +
</header>
 +
 +
 +
 +
 +
<!--Rebekah - START - insert contents here-->
 +
<section>
 +
<h3>Informing our Design: Modelling Michaelis-Menten Kinetics</h3>
 +
 +
<p>As we were refining the design of our plastic-sensing device, we were having a difficult time accounting for the lag between plastic uptake and plastic degradation. Recall that plastic in our bioreactor needs to be broken down by PETase in order to produce the terephthalic acid necessary to power our sensor. Though PETase is the most efficient PET-degrading enzyme discovered to date,<sup>1</sup> it by no means degrades plastic instantaneously. This means that our devices gives off a signal at a later time and later location from where plastic was originally collected. </p>
 +
 +
<p>The team decided to turn to modelling to find quantitative solutions to this design challenge. We realized that understanding degradation rates as determined by enzyme kinetics would allow us to adjust our reactor design to minimize the lag in signal. </p>
 +
 +
<p>The task required us to understand enzyme kinetics under two broad sets of circumstances: </p>
 +
<ul style="list-style-type:disc; font-size: 20px; color: #545454">
 +
<li>When the pieces of plastic were small enough in size so that they were effectively suspended in a solution. </li>
 +
<li>When the pieces of plastic were too large to be considered in suspension. This would require us to apply models geared towards heterogeneous reactions that involve interactions between a solid substrate and an enzyme in solution. </li>
 +
</ul>
 +
<br>
 +
<p>To understand conditions involving small pieces of plastic, our team relied on Michaelis-Menten, a widely-applied model for enzyme kinetics. Prior to beginning our work, we ensured that the assumptions inherent in the model were reasonable for modelling PETase activity:<sup>2</sup></p>
 +
<ul style="list-style-type:disc; font-size: 20px; color: #545454">
 +
<li>The model assumes that the product is not converted back into the substrate. This assumption holds when the concentration of product is very low. It is valid at the beginning of reactions and is especially compatible with our reactor design considering that the products of the reaction are funneled to another chamber where none of the original substrate or enzyme are present. The assumption also makes sense considering that our enzyme functions by breaking down polymers: after the reaction, the newly formed monomers no longer fit in the enzyme's active site. </li>
 +
<li>Another assumption is that the enzyme must be either free or bound to the substrate. That is, PETase should not at any time be bound to a molecule that does not consist of PET. The pnpb assays of Yoshida et al. (2016)<sup>1</sup> during their characterization of different esterases demonstrate how PETase is significantly less likely that other cutinases to catalyze the breakdown of esters. This suggests that PETase has a more specific binding site and comes closer than alternative enzymes to satisfying the assumption. </li>
 +
</ul>
 +
<br>
 +
<p>An overview of Michaelis-Menten:<sup>3</sup> The Michaelis-Menten constant (K<sub>m</sub>) corresponds to the substrate concentration at which the speed of the reaction is half of the maximum possible velocity: </p>
 +
 +
 +
<center><img src="https://static.igem.org/mediawiki/2016/0/02/T--Harvard_BioDesign--Modelling_MichaelisGeneric.png"/></center>
 +
 +
<p>Our team initially intended to compare the Michaelis-Menten constant associated with small plastics suspended in solutions to that of large pieces of plastics. However, we soon realized that the Michaelis-Menten model, with its references to substrate concentration, assumes that the substrate is in solution and would therefore be unsuitable for describing the kinetics associated with large solid substrates. </p>
 +
 +
<p>As our team and mentors investigated more sophisticated variations of Michaelis-Menten, we uncovered a highly relevant paper, "Surface Enzyme Kinetics for Biopolymer Microarrays: a Combination of Langmuir and Michaelis-Menten Concepts"<sup>4</sup>. This model incorporates additional parameters - adsorption and desorption coefficients - to help explain surface kinetics. Using the paper as a blueprint, our team constructed a model in Matlab (code available at the bottom of this page) that enables its users to enter values for adsorption, desorption, and catalysis constants and values for initial substrate and enzyme concentrations to immediately receive values for the Michaelis-Menten constant and the percentage of substrate that is bound to the enzyme at any point in time. </p>
 +
 +
<p>Since our assessment of enzyme kinetics was purely comparative (that is, we judged the efficiency of surface kinetics by comparing them to kinetics in solution), we focused our investigation on the parameters that are unique to surfaces: adsorption and desorption.
 +
</p>
 +
 +
</section>
 +
 +
<section>
 +
<h3>Supporting Future Experiments: Calculating Extinction Coefficients</h3>
 +
 +
<p>During our wet lab work with PETase this summer, we <a href="https://2016.igem.org/Team:Harvard_BioDesign/Results">managed to successfully isolate the enzyme</a>. In future quantitative experiments involving PETase, it would be useful to have a means by which we could calculate the concentration of this enzyme in a given solution. Although we have the necessary tools to measure the optical density (OD) of a solution that contains PETase, OD measurements are not sufficient on their own to determine concentration because we do not know how much light the enzyme absorbs.</p>
 +
 +
<p>This is where modelling comes into the picture. In order to calculate the concentration of enzyme present, our team also needs to know PETase's extinction coefficient. The extinction coefficient of an enzyme describes how much light the enzyme absorbs at a given wavelength. The paper, "Calculation of Protein Extinction Coefficients from Amino Acid Sequence Data",<sup>5</sup> guided us through the process of finding this crucial coefficient for a wavelength of 280 nm and allowed us to use the number of tryptophan, tyrosine, and cystine amino acids present in PETase to arrive at an extinction coefficient of 39170 M<sup>-1</sup> cm<sup>-1</sup>. </p>
 +
 +
<p>Below is the equation from the paper that enabled us to determine the extinction coefficient: </p>
 +
 +
<img src="https://static.igem.org/mediawiki/2016/0/0e/T--Harvard_BioDesign--Modelling_ExtinctionEq1.png"/><br>
 +
<img src="https://static.igem.org/mediawiki/2016/8/86/T--Harvard_BioDesign--Modelling_subeq1.png"/><br>
 +
<img src="https://static.igem.org/mediawiki/2016/2/27/T--Harvard_BioDesign--Modelling_subeq2.png"/><br>
 +
<img src="https://static.igem.org/mediawiki/2016/8/88/T--Harvard_BioDesign--Modelling_subeq3.png"/>
 +
<p>Although this technically the coefficient for the denatured protein, the paper states explicitly that there is a "relatively small difference" between this coefficient and the coefficient of the native protein at 280 nm. It states that it is reasonably accurate to calculate it "directly... from amino acid composition information".</p>
 +
 +
</section>
 +
 +
<section>
 +
<h3>MATLAB Code</h3>
 +
<!--This is the code for the Matlab model: could either include in body of text in a different font/color/size, or could include as a downloadable matlab file-->
  
<h5> Inspiration </h5>
 
 
<p>
 
<p>
Here are a few examples from previous teams:
+
<br>%% PET Surface Enzyme Kinetics
 +
 +
<br>% Define variables
 +
<br>h = 0.001; % step size
 +
<br>t = 1:h:100; % time frame
 +
 +
<br>E = 0*t; %concentration of enzyme
 +
 +
<br>k_cat = 0; % enzyme catalysis
 +
<br>k_a = 0; % enzyme adsorption
 +
<br>k_d = 0; % enzyme desorption
 +
 +
<br>gamma_S = 0*t; % surface coverage of substrate
 +
<br>gamma_ES = 0*t; % surface coverage of enzyme-substrate complex
 +
<br>gamma_P = 0*t; % surface coverage of product
 +
<br>gamma_total = gamma_S + gamma_ES + gamma_P; % total number of surface sites
 +
 +
<br>d_gamma_S = diff(gamma_S)/h % derivative of surface coverage of substrate wrt t
 +
<br>d_gamma_ES = diff(gamma_ES)/h % derivative of surface coverage of complex wrt t
 +
<br>d_gamma_P = diff(gamma_P)/h % derivative of surface coverage of product wrt t
 +
 +
<br>theta_S = gamma_S/gamma_total; % relative surface coverage of substrate
 +
<br>theta_ES = gamma_ES/gamma_total; % relative surface coverage of complex
 +
<br>theta_P = gamma_P/gamma_total; % relative surface coverage of product
 +
 +
<br>d_theta_S = diff(theta_S)/h % derivative of relative coverage of substrate wrt t
 +
<br>d_theta_ES = diff(theta_ES)/h % derivative of relative coverage of complex wrt t
 +
<br>d_theta_P = diff(theta_P)/h % derivative of relative coverage of product wrt t
 +
 +
<br>% Set-up
 +
<br>d_gamma_ES = k_a*gamma_S - k_d*gamma_ES - k_cat*gamma_ES; % rxn rate for production of complex
 +
<br>d_gamma_P = k_cat*gamma_ES; % rxn rate for production of product
 +
<br>% Rates of change of surface coverage
 +
<br>d_theta_ES = k_a*theta_S - k_d*theta_ES - k_cat*theta_ES;
 +
<br>d_theta_P = k_cat*theta_ES;
 +
 +
<br>% Assuming k_cat << k_a, k_d
 +
<br>theta_ES = k_a*theta_S*E/(k_d + k_cat); % equilibrium surface coverage of complex
 +
<br>K_m = k_d + k_cat/(k_a); % Michaelis-Menten constant
 
</p>
 
</p>
<ul>
+
</section>
<li><a href="https://2014.igem.org/Team:ETH_Zurich/modeling/overview">ETH Zurich 2014</a></li>
+
 
<li><a href="https://2014.igem.org/Team:Waterloo/Math_Book">Waterloo 2014</a></li>
+
<section>
</ul>
+
<h3>Footnotes</h3>
 +
 
 +
<p>
 +
<sup>1</sup>Yoshida, Shosuke et al. “A Bacterium That Degrades and Assimilates Poly(ethylene Terephthalate).” Science 351.6278 (2016): 1196–1199. science.sciencemag.org. Web.<br>
 +
 
 +
<sup>2</sup>Montagud, A., E. Navarro, P. Córdoba, and J. Urchueguía. Kinetic Chemistry. N.p.: InterTech, n.d. Openwetware.org. Web. 6 July 2016.<br>
 +
 
 +
<sup>3</sup>"Michaelis-Menten Equation - Interactive Graph." PhysiologyWeb. N.p., 1 Sept. 2014. Web. 14 July 2016.<br>
 +
 
 +
<sup>4</sup>Lee, Hye, Alastair Wark, Terry Goodrich, Shiping Fang, and Robert Corn. "Surface Enzyme Kinetics for Biopolymer Microarrays: A Combination of Langmuir and Michaelis-Menten Concepts." American Chemical Society (2005): n. pag. Web. 26 July 2016.
 +
 
 +
<sup>5</sup> Gill, Stanley, and Peter Hippel. "Calculation of Protein Extinction Coefficients from Amino Acid Sequence Data." Analytical Biochemistry (1989): n. pag. Web.
 +
</p>
 +
 
 +
</section>
 +
 
 +
<!-- Rebekah - STOP -->
 +
 
 +
 
 +
 
 +
 
 +
 
  
  
</div>
+
</article>
 +
</div>
  
 +
</div>
 +
<!--y Does it start here?-->
 +
<!-- Footer -->
 +
            <div id="footer">
 +
                <div class="container">
 +
                    <div class="row">
 +
                        <div class="12u">
 +
                           
 +
                            <!-- Contact -->
 +
                            <section class="contact">
 +
                               
 +
                                <ul class="icons">
 +
                                   
 +
                                    <li>
 +
                                        <a href="https://www.facebook.com/Harvard-iGEM-269794024286" target="_blank">
 +
                                            <img src="https://static.igem.org/mediawiki/2016/2/2f/T--Harvard_BioDesign--images_facebook01.png"alt="Facebook Logo" style="width:51px;height:51px;">
 +
                                                </a>
 +
                                    </li>
 +
                                   
 +
                                    <li>
 +
                                        <a href="mailto:crimsonigem@gmail.com">
 +
                                            <img src="https://static.igem.org/mediawiki/2016/e/e2/T--Harvard_BioDesign--images_gmail01.png" alt="Email Logo" style="width:51px;height:51px;">
 +
                                                </a>
 +
                                    </li>
 +
                                   
 +
                                   
 +
                                    <li>
 +
                                        <a href="https://twitter.com/harvardigem" target="_blank">
 +
                                            <img src="https://static.igem.org/mediawiki/2016/4/4e/T--Harvard_BioDesign--images_twitter01.png"alt="Twitter Logo" style="width:51px;height:51px;">
 +
                                                </a>
 +
                                    </li>
 +
                                   
 +
                                   
 +
                                </ul>
 +
                            </section>
 +
                           
 +
                            <!-- Copyright -->
 +
                            <div class="copyright">
 +
                                <ul class="menu">
 +
                                    <li>&copy; 2016 Harvard iGEM</li><li>Design: <a href="http://html5up.net">HTML5 UP</a></li>
 +
                                </ul>
 +
                               
 +
                               
 +
                            </div>
 +
                           
 +
                        </div>
 +
                    </div>
 +
                </div>
 +
               
 +
</div>
 +
</body>
 
</html>
 
</html>

Latest revision as of 12:51, 19 January 2017

Harvard BioDesign 2016

Model

Informing our Design: Modelling Michaelis-Menten Kinetics

As we were refining the design of our plastic-sensing device, we were having a difficult time accounting for the lag between plastic uptake and plastic degradation. Recall that plastic in our bioreactor needs to be broken down by PETase in order to produce the terephthalic acid necessary to power our sensor. Though PETase is the most efficient PET-degrading enzyme discovered to date,1 it by no means degrades plastic instantaneously. This means that our devices gives off a signal at a later time and later location from where plastic was originally collected.

The team decided to turn to modelling to find quantitative solutions to this design challenge. We realized that understanding degradation rates as determined by enzyme kinetics would allow us to adjust our reactor design to minimize the lag in signal.

The task required us to understand enzyme kinetics under two broad sets of circumstances:

  • When the pieces of plastic were small enough in size so that they were effectively suspended in a solution.
  • When the pieces of plastic were too large to be considered in suspension. This would require us to apply models geared towards heterogeneous reactions that involve interactions between a solid substrate and an enzyme in solution.

To understand conditions involving small pieces of plastic, our team relied on Michaelis-Menten, a widely-applied model for enzyme kinetics. Prior to beginning our work, we ensured that the assumptions inherent in the model were reasonable for modelling PETase activity:2

  • The model assumes that the product is not converted back into the substrate. This assumption holds when the concentration of product is very low. It is valid at the beginning of reactions and is especially compatible with our reactor design considering that the products of the reaction are funneled to another chamber where none of the original substrate or enzyme are present. The assumption also makes sense considering that our enzyme functions by breaking down polymers: after the reaction, the newly formed monomers no longer fit in the enzyme's active site.
  • Another assumption is that the enzyme must be either free or bound to the substrate. That is, PETase should not at any time be bound to a molecule that does not consist of PET. The pnpb assays of Yoshida et al. (2016)1 during their characterization of different esterases demonstrate how PETase is significantly less likely that other cutinases to catalyze the breakdown of esters. This suggests that PETase has a more specific binding site and comes closer than alternative enzymes to satisfying the assumption.

An overview of Michaelis-Menten:3 The Michaelis-Menten constant (Km) corresponds to the substrate concentration at which the speed of the reaction is half of the maximum possible velocity:

Our team initially intended to compare the Michaelis-Menten constant associated with small plastics suspended in solutions to that of large pieces of plastics. However, we soon realized that the Michaelis-Menten model, with its references to substrate concentration, assumes that the substrate is in solution and would therefore be unsuitable for describing the kinetics associated with large solid substrates.

As our team and mentors investigated more sophisticated variations of Michaelis-Menten, we uncovered a highly relevant paper, "Surface Enzyme Kinetics for Biopolymer Microarrays: a Combination of Langmuir and Michaelis-Menten Concepts"4. This model incorporates additional parameters - adsorption and desorption coefficients - to help explain surface kinetics. Using the paper as a blueprint, our team constructed a model in Matlab (code available at the bottom of this page) that enables its users to enter values for adsorption, desorption, and catalysis constants and values for initial substrate and enzyme concentrations to immediately receive values for the Michaelis-Menten constant and the percentage of substrate that is bound to the enzyme at any point in time.

Since our assessment of enzyme kinetics was purely comparative (that is, we judged the efficiency of surface kinetics by comparing them to kinetics in solution), we focused our investigation on the parameters that are unique to surfaces: adsorption and desorption.

Supporting Future Experiments: Calculating Extinction Coefficients

During our wet lab work with PETase this summer, we managed to successfully isolate the enzyme. In future quantitative experiments involving PETase, it would be useful to have a means by which we could calculate the concentration of this enzyme in a given solution. Although we have the necessary tools to measure the optical density (OD) of a solution that contains PETase, OD measurements are not sufficient on their own to determine concentration because we do not know how much light the enzyme absorbs.

This is where modelling comes into the picture. In order to calculate the concentration of enzyme present, our team also needs to know PETase's extinction coefficient. The extinction coefficient of an enzyme describes how much light the enzyme absorbs at a given wavelength. The paper, "Calculation of Protein Extinction Coefficients from Amino Acid Sequence Data",5 guided us through the process of finding this crucial coefficient for a wavelength of 280 nm and allowed us to use the number of tryptophan, tyrosine, and cystine amino acids present in PETase to arrive at an extinction coefficient of 39170 M-1 cm-1.

Below is the equation from the paper that enabled us to determine the extinction coefficient:




Although this technically the coefficient for the denatured protein, the paper states explicitly that there is a "relatively small difference" between this coefficient and the coefficient of the native protein at 280 nm. It states that it is reasonably accurate to calculate it "directly... from amino acid composition information".

MATLAB Code


%% PET Surface Enzyme Kinetics
% Define variables
h = 0.001; % step size
t = 1:h:100; % time frame
E = 0*t; %concentration of enzyme
k_cat = 0; % enzyme catalysis
k_a = 0; % enzyme adsorption
k_d = 0; % enzyme desorption
gamma_S = 0*t; % surface coverage of substrate
gamma_ES = 0*t; % surface coverage of enzyme-substrate complex
gamma_P = 0*t; % surface coverage of product
gamma_total = gamma_S + gamma_ES + gamma_P; % total number of surface sites
d_gamma_S = diff(gamma_S)/h % derivative of surface coverage of substrate wrt t
d_gamma_ES = diff(gamma_ES)/h % derivative of surface coverage of complex wrt t
d_gamma_P = diff(gamma_P)/h % derivative of surface coverage of product wrt t
theta_S = gamma_S/gamma_total; % relative surface coverage of substrate
theta_ES = gamma_ES/gamma_total; % relative surface coverage of complex
theta_P = gamma_P/gamma_total; % relative surface coverage of product
d_theta_S = diff(theta_S)/h % derivative of relative coverage of substrate wrt t
d_theta_ES = diff(theta_ES)/h % derivative of relative coverage of complex wrt t
d_theta_P = diff(theta_P)/h % derivative of relative coverage of product wrt t
% Set-up
d_gamma_ES = k_a*gamma_S - k_d*gamma_ES - k_cat*gamma_ES; % rxn rate for production of complex
d_gamma_P = k_cat*gamma_ES; % rxn rate for production of product
% Rates of change of surface coverage
d_theta_ES = k_a*theta_S - k_d*theta_ES - k_cat*theta_ES;
d_theta_P = k_cat*theta_ES;
% Assuming k_cat << k_a, k_d
theta_ES = k_a*theta_S*E/(k_d + k_cat); % equilibrium surface coverage of complex
K_m = k_d + k_cat/(k_a); % Michaelis-Menten constant

Footnotes

1Yoshida, Shosuke et al. “A Bacterium That Degrades and Assimilates Poly(ethylene Terephthalate).” Science 351.6278 (2016): 1196–1199. science.sciencemag.org. Web.
2Montagud, A., E. Navarro, P. Córdoba, and J. Urchueguía. Kinetic Chemistry. N.p.: InterTech, n.d. Openwetware.org. Web. 6 July 2016.
3"Michaelis-Menten Equation - Interactive Graph." PhysiologyWeb. N.p., 1 Sept. 2014. Web. 14 July 2016.
4Lee, Hye, Alastair Wark, Terry Goodrich, Shiping Fang, and Robert Corn. "Surface Enzyme Kinetics for Biopolymer Microarrays: A Combination of Langmuir and Michaelis-Menten Concepts." American Chemical Society (2005): n. pag. Web. 26 July 2016. 5 Gill, Stanley, and Peter Hippel. "Calculation of Protein Extinction Coefficients from Amino Acid Sequence Data." Analytical Biochemistry (1989): n. pag. Web.