Evan pepper (Talk | contribs) |
Evan pepper (Talk | contribs) |
||
Line 742: | Line 742: | ||
<br> | <br> | ||
<h3>Background</h3> | <h3>Background</h3> | ||
− | <p>The purpose of modeling is to carefully examine the pathways of each intended biosynthetic products, look for ways to optimize production, understand limiting factors | + | <p>The purpose of modeling is to carefully examine the pathways of each intended biosynthetic products, look for ways to optimize production, and understand limiting factors. To accomplish these goals, we used available metabolic pathways for our target organism, and evaluated several different methods to model production of acetaminophen and B<sub>12</sub> in cyanobacteria. Each of these modeling methods has different assumptions which allow these data to be averaged; providing reasonable quantitative estimates of our biosynthetic products.</p> |
<div class="container"> | <div class="container"> | ||
<div class="row"> | <div class="row"> | ||
Line 756: | Line 756: | ||
<br> | <br> | ||
<h2 style="text-align: left; font-weight: 500;">Overview</h2> | <h2 style="text-align: left; font-weight: 500;">Overview</h2> | ||
− | <p>To predict acetaminophen biosynthesis, we analyzed the abundance of the acetaminophen's precursor, chorismate. | + | <p>To predict acetaminophen biosynthesis, we analyzed the abundance of the acetaminophen's precursor, chorismate. Chorismate is primarily used by the cell to produce the three aromatic amino acids, but is also used to synthesize salicylic acid, folate, vitamins, and many alkaloids<sup>[1]</sup>. We used public data on these chorismate products and gene transcription of amino acids for multiple species of cyanobacteria to estimate the chorismate pool available for acetaminophen production. Using enzyme binding affinity K<sub>m</sub>'s we created a simulation of chorismate metabolism in Python to approximate enzyme rates. Insufficient data on cyanobacterial K<sub>m</sub> values necessitated using data from other related bacterial species. We assumed the enzymes with the highest matching sequence identity to <i>S. elongatus</i> PCC 7942 and <i>A. platensis</i> would have K<sub>m</sub>'s close to the enzyme rates of <i>S. elongatus</i> PCC 7942 and <i>A. platensis</i>. With this information, we made a quantitative metabolic model which estimates how much of our precursor goes down each pathway.</p> |
<br> | <br> | ||
<img src="https://static.igem.org/mediawiki/2017/7/7a/Chorismate-metabolism.png" style="width:70%;"> | <img src="https://static.igem.org/mediawiki/2017/7/7a/Chorismate-metabolism.png" style="width:70%;"> | ||
− | <figcaption>Chorismate is processed into the aromatic amino acids, phenylalanine, tyrosine, tryptophan, and folate. Our inserted enzyme <i>4ABH</i> metabolize PABA and anthranilate to make 4-aminophenol which is then processed by <i> | + | <figcaption>Chorismate is processed into the aromatic amino acids, phenylalanine, tyrosine, tryptophan, and folate. Our inserted enzyme <i>4ABH</i> metabolize PABA and anthranilate to make 4-aminophenol which is then processed by <i>nhoA</i> to make acetaminophen, top right. Note that unmarked arrows signify unmodified enzymatic reactions.</figcaption> |
<h2 style="text-align: left; font-weight: 500;">Enzyme Competition</h2> | <h2 style="text-align: left; font-weight: 500;">Enzyme Competition</h2> | ||
<p>Precursor concentrations would be the first limit to how much acetaminophen could be produced. The second would be how effective our enzyme <i>4ABH</i> was at converting precursors from the tryptophan and folate pathways into 4-aminophenol. To answer the competition question, we can either assume all three competing enzymes will have similar precursor affinity and produce each product in equal quantities or use ratios of each competing enzyme's K<sub>m</sub> for the limiting intermediates. This enzyme kinetics based method is difficult because there is so little data on rates, quantities, and affinities for our particular enzyme sequences. For that reason, all of the following calulations will be computed using a 33% assumed precursor to product conversion rate.</p> | <p>Precursor concentrations would be the first limit to how much acetaminophen could be produced. The second would be how effective our enzyme <i>4ABH</i> was at converting precursors from the tryptophan and folate pathways into 4-aminophenol. To answer the competition question, we can either assume all three competing enzymes will have similar precursor affinity and produce each product in equal quantities or use ratios of each competing enzyme's K<sub>m</sub> for the limiting intermediates. This enzyme kinetics based method is difficult because there is so little data on rates, quantities, and affinities for our particular enzyme sequences. For that reason, all of the following calulations will be computed using a 33% assumed precursor to product conversion rate.</p> | ||
+ | |||
+ | <br> | ||
<p>Since few cyanobacterial proteins have been isolated and tested for kinetic data, many of the K<sub>m</sub>s compared are from other species. <i>4ABH</i> has a K<sub>m</sub> of 20.4µM for p-amino benzoate(PABA), while the folate processing enzyme folp1 has a K<sub>m</sub> of 0.37µM<sup>[5,6]</sup>. Assuming the rate and quantity are the same, the folate enzyme will be fully saturated at a much lesser concentration of PABA than <i>4ABH</i>. You can assume that the logrithmic K<sub>m</sub> ratio of 20.4:0.37 (or 55:1) will represent how many moles of PABA goes to folate versus to acetaminophen. For the tryptophan intermediate anthranilate, no K<sub>m</sub> is given but <i>4ABH</i> is assayed as having a rate of 24% of PABA, equal to 34µM for anthranilate, this compares to <i>TrpD's</i> K<sub>m</sub> for chorismate of 40nM, meaning we'd have 850:1 tryptophan to acetaminophen. Since these ratios are based on poorly matching enzymes and make acetaminophen production undetectably small, we will keep this enzyme competition in mind while using 33% conversion for simplicity.</p> | <p>Since few cyanobacterial proteins have been isolated and tested for kinetic data, many of the K<sub>m</sub>s compared are from other species. <i>4ABH</i> has a K<sub>m</sub> of 20.4µM for p-amino benzoate(PABA), while the folate processing enzyme folp1 has a K<sub>m</sub> of 0.37µM<sup>[5,6]</sup>. Assuming the rate and quantity are the same, the folate enzyme will be fully saturated at a much lesser concentration of PABA than <i>4ABH</i>. You can assume that the logrithmic K<sub>m</sub> ratio of 20.4:0.37 (or 55:1) will represent how many moles of PABA goes to folate versus to acetaminophen. For the tryptophan intermediate anthranilate, no K<sub>m</sub> is given but <i>4ABH</i> is assayed as having a rate of 24% of PABA, equal to 34µM for anthranilate, this compares to <i>TrpD's</i> K<sub>m</sub> for chorismate of 40nM, meaning we'd have 850:1 tryptophan to acetaminophen. Since these ratios are based on poorly matching enzymes and make acetaminophen production undetectably small, we will keep this enzyme competition in mind while using 33% conversion for simplicity.</p> | ||
Line 778: | Line 780: | ||
<h2 style="text-align: left; font-weight: 500;">Amino Acid Method</h2> | <h2 style="text-align: left; font-weight: 500;">Amino Acid Method</h2> | ||
− | <p> | + | <p>Because chorismate is used as a precursor for the aromatic amino acids, we used the amino acid concentration as a proxy for the chorismate concentration. The amino acid composition for <i>A. platensis</i> is publicly available, and we used that to calculate the moles per cell of the amino acid precursors anthranilate and PABA. These molecules are the direct substrate for our enzyme <i>4ABH</i>, which converts anthranilate and PABA to the intermediate 4-aminophenol before <i>nhoA</i> converts that to acetaminophen. We show several different calculations below using different sources of data.</p> |
+ | |||
+ | <br> | ||
<h2 style="text-align: left; font-weight: 400;">Assumptions</h2> | <h2 style="text-align: left; font-weight: 400;">Assumptions</h2> | ||
<ul style="font-family: 'objektiv-mk1'; font-size: inherit; text-align: left;"> | <ul style="font-family: 'objektiv-mk1'; font-size: inherit; text-align: left;"> | ||
− | <li><i>Synechococcus</i> and <i>Arthrospira platensis</i> have similar amino acid ratios. Since there was no available amino acid data for | + | <li><i>Synechococcus</i> and <i>Arthrospira platensis</i> have similar amino acid ratios. Since there was no available amino acid data for <i>Synechococcus</i>, we assumed it has a similar ratio to the more well described <i>Arthrospira platensis</i>.</li> |
<br> | <br> | ||
<li>The amount of anthranilate and PABA precursors are equal to that of their products tryptophan and folate.</li> | <li>The amount of anthranilate and PABA precursors are equal to that of their products tryptophan and folate.</li> | ||
<br> | <br> | ||
− | <li> | + | <li>Based on the K<sub>m</sub> ratios of <i>4ABH</i> and its competitor <i>TrpD</i>, 30% of the available precursors will go down the acetaminophen pathway.</li> |
<br> | <br> | ||
<li>Folate is present in such small amounts, at 0.21nM of PABA, that its effect on acetaminophen production is neglegible<sup>[7]</sup>.</li> | <li>Folate is present in such small amounts, at 0.21nM of PABA, that its effect on acetaminophen production is neglegible<sup>[7]</sup>.</li> | ||
Line 828: | Line 832: | ||
<p>These numbers show that there will probably be enough precursor to produce a useful, detectable quantity of acetaminophen. Based on literature and sequence estimates of aromatic amino acids, we can assume there would be at least that many moles of chorismate from which our added pathway pushes towards acetaminophen. The three calculations above can be averaged to finally predict 8.26mg ± 2.77mg acetaminophen per gram of <i>Arthrospira platensis</i> biomass.</p> | <p>These numbers show that there will probably be enough precursor to produce a useful, detectable quantity of acetaminophen. Based on literature and sequence estimates of aromatic amino acids, we can assume there would be at least that many moles of chorismate from which our added pathway pushes towards acetaminophen. The three calculations above can be averaged to finally predict 8.26mg ± 2.77mg acetaminophen per gram of <i>Arthrospira platensis</i> biomass.</p> | ||
+ | |||
+ | <br> | ||
<p>We used different tryptophan estimates to reach several different predictions for acetaminophen, averaging 8.26mg ± 2.77mg acetaminophen per gram of <i>Arthrospira platensis</i> biomass or 8.26µg per mL. This would be significantly above the limit of detection for our HPLC, at 50ng per ml, and serve as a starting point for optimizing production. This means that one 325mg dose of acetaminophen could be obtained in ~39g of biomass, meaning a 12 by 3 feet round pool could produce enough acetaminophen for more than 200 people every 10 days. While 39 grams isn't an ideal amount of medicine to consume, it does show that <i>Arthrospira platensis</i> has significant potential as a molecular factory for acetaminophen.</p> | <p>We used different tryptophan estimates to reach several different predictions for acetaminophen, averaging 8.26mg ± 2.77mg acetaminophen per gram of <i>Arthrospira platensis</i> biomass or 8.26µg per mL. This would be significantly above the limit of detection for our HPLC, at 50ng per ml, and serve as a starting point for optimizing production. This means that one 325mg dose of acetaminophen could be obtained in ~39g of biomass, meaning a 12 by 3 feet round pool could produce enough acetaminophen for more than 200 people every 10 days. While 39 grams isn't an ideal amount of medicine to consume, it does show that <i>Arthrospira platensis</i> has significant potential as a molecular factory for acetaminophen.</p> | ||
Line 888: | Line 894: | ||
<br> | <br> | ||
<br> | <br> | ||
− | <p>To understand the production capacity of our organisms, we aggregated growth data from published papers<sup>[22,23,24]</sup> and our | + | <p>To understand the production capacity of our organisms, we aggregated growth data from published papers<sup>[22,23,24]</sup> and our own daily optical density growth data. Using carrying-capacity-limited logistic growth curves to fit our data to an equation, we modelled dried biomass with respect to time. We also added additional dependent variables of temperature and light intensity to our equation so that we can better understand how minimalist growth conditions might affect the maximum culture density and the maximum growth rate.</p> |
<br> | <br> |
Revision as of 18:20, 1 November 2017
![](https://static.igem.org/mediawiki/2017/f/fd/UCSCmodeling.png)
MODELING
Predict and optimize yield
Background
The purpose of modeling is to carefully examine the pathways of each intended biosynthetic products, look for ways to optimize production, and understand limiting factors. To accomplish these goals, we used available metabolic pathways for our target organism, and evaluated several different methods to model production of acetaminophen and B12 in cyanobacteria. Each of these modeling methods has different assumptions which allow these data to be averaged; providing reasonable quantitative estimates of our biosynthetic products.
ACETAMINOPHEN
![](https://upload.wikimedia.org/wikipedia/commons/thumb/2/29/Paracetamol-skeletal.svg/1200px-Paracetamol-skeletal.svg.png)
Overview
To predict acetaminophen biosynthesis, we analyzed the abundance of the acetaminophen's precursor, chorismate. Chorismate is primarily used by the cell to produce the three aromatic amino acids, but is also used to synthesize salicylic acid, folate, vitamins, and many alkaloids[1]. We used public data on these chorismate products and gene transcription of amino acids for multiple species of cyanobacteria to estimate the chorismate pool available for acetaminophen production. Using enzyme binding affinity Km's we created a simulation of chorismate metabolism in Python to approximate enzyme rates. Insufficient data on cyanobacterial Km values necessitated using data from other related bacterial species. We assumed the enzymes with the highest matching sequence identity to S. elongatus PCC 7942 and A. platensis would have Km's close to the enzyme rates of S. elongatus PCC 7942 and A. platensis. With this information, we made a quantitative metabolic model which estimates how much of our precursor goes down each pathway.
![](https://static.igem.org/mediawiki/2017/7/7a/Chorismate-metabolism.png)
Enzyme Competition
Precursor concentrations would be the first limit to how much acetaminophen could be produced. The second would be how effective our enzyme 4ABH was at converting precursors from the tryptophan and folate pathways into 4-aminophenol. To answer the competition question, we can either assume all three competing enzymes will have similar precursor affinity and produce each product in equal quantities or use ratios of each competing enzyme's Km for the limiting intermediates. This enzyme kinetics based method is difficult because there is so little data on rates, quantities, and affinities for our particular enzyme sequences. For that reason, all of the following calulations will be computed using a 33% assumed precursor to product conversion rate.
Since few cyanobacterial proteins have been isolated and tested for kinetic data, many of the Kms compared are from other species. 4ABH has a Km of 20.4µM for p-amino benzoate(PABA), while the folate processing enzyme folp1 has a Km of 0.37µM[5,6]. Assuming the rate and quantity are the same, the folate enzyme will be fully saturated at a much lesser concentration of PABA than 4ABH. You can assume that the logrithmic Km ratio of 20.4:0.37 (or 55:1) will represent how many moles of PABA goes to folate versus to acetaminophen. For the tryptophan intermediate anthranilate, no Km is given but 4ABH is assayed as having a rate of 24% of PABA, equal to 34µM for anthranilate, this compares to TrpD's Km for chorismate of 40nM, meaning we'd have 850:1 tryptophan to acetaminophen. Since these ratios are based on poorly matching enzymes and make acetaminophen production undetectably small, we will keep this enzyme competition in mind while using 33% conversion for simplicity.
Assumptions
- We assume that our genes were designed, inserted, and translated successfully and in reasonable quantities.
- Km values for Arthrospira platensis can be approximated by using other organisms' enzyme rates. Sequences were BLAST compared in each case, resulting in a range of similarity between 52% and 63% positive alignment, suggesting that enzyme orthologs might be significantly different.
- We assume that enzyme rate and quantity is the same for each 4ABH, folp1, and TrpD, else we could not compare Km ratios and calculate approximate the Km of 4ABH for anthranilate.
Amino Acid Method
Because chorismate is used as a precursor for the aromatic amino acids, we used the amino acid concentration as a proxy for the chorismate concentration. The amino acid composition for A. platensis is publicly available, and we used that to calculate the moles per cell of the amino acid precursors anthranilate and PABA. These molecules are the direct substrate for our enzyme 4ABH, which converts anthranilate and PABA to the intermediate 4-aminophenol before nhoA converts that to acetaminophen. We show several different calculations below using different sources of data.
Assumptions
- Synechococcus and Arthrospira platensis have similar amino acid ratios. Since there was no available amino acid data for Synechococcus, we assumed it has a similar ratio to the more well described Arthrospira platensis.
- The amount of anthranilate and PABA precursors are equal to that of their products tryptophan and folate.
- Based on the Km ratios of 4ABH and its competitor TrpD, 30% of the available precursors will go down the acetaminophen pathway.
- Folate is present in such small amounts, at 0.21nM of PABA, that its effect on acetaminophen production is neglegible[7].
$$\frac{0.442\ µmol\ Trp}{1g\ biomass}\approx \frac{0.455\ µmol\ anth}{1g\ biomass}\times\frac{1\ mol\ acet}{3\ moles\ anth}\time\frac{151g\ acet.}{1 mol acet.}=\frac{2.3mg\ acet.}{1g\ biomass}$$
Sequence Analysis Method
To validate our organism's quantity of tryptophan precursor, we used a custom Python program[9] to convert DNA sequences to amino acids and calculate molar and mass percentages of tryptophan. We ran both the genome and all 55 listed ribosomal protein sequences through our program, which resulted in 0.9% and 0.6% tryptophan by moles. Knowing that about 60% of Arthrospira platensis is protein by mass, we can predict acetaminophen production.
$$\frac{0.056\ g\ Trp}{1\ biomass}\rightarrow\frac{0.27\ mmol\ Trp}{1\ g\ biomass}\times\frac{1\ mol\ acet}{3\ mol\ Trp}\times\frac{151.163g\ acet}{1\ mol\ acet.}=\frac{13.7mg\ acet.}{1g\ biomass}$$
$$\frac{0.036 g\ Trp}{1 g\ biomass}\rightarrow \frac{0.0.175\ mmol\ Trp}{1\ g\ biomass} * \frac{1\ mol\ acet}{3\ mol\ chor} *\frac{151.163\ g}{1 mol\ acet} = \frac{8.8\ mg\ acet}{1 g\ biomass}$$
Assumptions
- Codon composition from the genome and ribosomal proteins alone apprimates the amino acid composition of the cyanobacteria.
- The amount of anthranilate and PABA precursors are equal to that of their products tryptophan and folate.
- Of the available precursors, 33% will go down our pathway. This is based off the idea that the enzymes may have similar affinities for the precursor.
- Folate is present in such small amounts, at 0.21nM of PABA, that its effect on acetaminophen production is neglegible[7].
These numbers show that there will probably be enough precursor to produce a useful, detectable quantity of acetaminophen. Based on literature and sequence estimates of aromatic amino acids, we can assume there would be at least that many moles of chorismate from which our added pathway pushes towards acetaminophen. The three calculations above can be averaged to finally predict 8.26mg ± 2.77mg acetaminophen per gram of Arthrospira platensis biomass.
We used different tryptophan estimates to reach several different predictions for acetaminophen, averaging 8.26mg ± 2.77mg acetaminophen per gram of Arthrospira platensis biomass or 8.26µg per mL. This would be significantly above the limit of detection for our HPLC, at 50ng per ml, and serve as a starting point for optimizing production. This means that one 325mg dose of acetaminophen could be obtained in ~39g of biomass, meaning a 12 by 3 feet round pool could produce enough acetaminophen for more than 200 people every 10 days. While 39 grams isn't an ideal amount of medicine to consume, it does show that Arthrospira platensis has significant potential as a molecular factory for acetaminophen.
VITAMIN B12
![](https://static.igem.org/mediawiki/2017/6/62/B12structure.png)
The quantity of DMB B12 produced depends on a successful production and integration of the active B12 lower ligand, 5,6-dimethyl-benzimidazole (5,6-DMB). For phytoplankton in the wild, cobalt is often the limiting factor for growth and production of B12[11,12], while B12 production is limited by growth need in optimal media[12]. With ssuE and bluB genes inserted and regulated using a strong PrtC promoter, the activating lower ligand 5,6-DMB will be created in abundance[13]. Synechococcus and Arthrospira platensis both have CobS, bluB, and pGam genes that code for proteins which bind 5,6-DMB to the cobalt[14,15]. If these proteins work as well as in their origin organism, then assays report that 5,6-DMB has at least 100 times higher affinity for cobalt than the B12 analog ligand, adenine[16,17], meaning the DMB B12 to B12 analog ratio would be 100:1. Published HPLC results show that that Arthrospira platensis produces between 1.5-2.5µg B12 analogs per gram dry weight[18].
$$\frac{2.5µg\ B_{12}\ analog}{1g\ drymass} * \frac{100\ DMB\ bindings}{101\ DMB+adenine\ binding} = \frac{2.47µg\ DMB\ B_{12}}{1g\ biomass}$$
An additional paper assayed Synechococcus elongatus sp WH 7803 at 10-18 moles per cell[16], and at a reported density 109 cells per liter (about a gram)[16,17], Synechococcus would produce 1.3µg B12 per liter dry mass. Another older paper used microbiological assays along with TLC and HPLC, finding maximums of 2.4, 1.47, and 1.27µg B12 per gram dry mass. Averaging the 6 data points and multiplying by a 100:1 conversion ratio results in a predicted production of 1.74 ± 0.23µg DMB B12 per gram of Arthrospira platensis drymass, meaning the USDA’s recommended daily value of 6µg could be obtained in one 3.5 gram serving.
Assumptions
- Gene inserts will be expressed, converting riboflavin-5′-phosphate to 5,6-DMB in excess[13].
- The bluB/CobS protein complex in Arthrospira platensis will attach 5,6-DMB to cobalt at rates similar to those assayed in Propionibacterium freudenreichii[14,15].
- Cobalt will be provided in excess of 0.3mM according the BG-11 recipe, ensuring maximum precursor availability[16,17]
- Synechococcus 7942, 7803 and Arthrospira platensis will have similar rates of B12 production.
Future B12 projects might use chemo-trophic bacteria such as Methanosarcina barkeri which produces more than 1000 times more B12 and could be converted to active form using a similar bioengineering proccess[21]. One exceptional use of cyanobacterial B12 is growing cyanobacteria in the water used to grow rice, increasing the carbon fixation, nitrogen fixation, and bioenriching the rice with B12[18].
BIOMASS
![](https://static.igem.org/mediawiki/2017/5/5c/Spirulina_powder.png)
To understand the production capacity of our organisms, we aggregated growth data from published papers[22,23,24] and our own daily optical density growth data. Using carrying-capacity-limited logistic growth curves to fit our data to an equation, we modelled dried biomass with respect to time. We also added additional dependent variables of temperature and light intensity to our equation so that we can better understand how minimalist growth conditions might affect the maximum culture density and the maximum growth rate.
Light Intensity: μE m-2 s-1
Temperature: ℃