Difference between revisions of "Team:UCSC/Model"
| Line 851: | Line 851: | ||
</div> | </div> | ||
| − | <p>The quantity of DMB B<sub>12</sub> produced depends on a successful production and integration of the active B<sub>12</sub> lower ligand, composed of 5,6-dimethylbenzimidazole (5,6-DMB) and α-ribazole 5’-phosphate. For cyanobacteria in the wild, cobalt is often the limiting factor for growth and production of B<sub>12</sub><sup>[11,12]</sup>. Both <i>ssuE</i> and <i>bluB</i> genes are engineered to be regulated using a strong PtrC promoter, creating an abundance of the activating portion of the lower ligand, 5,6-DMB. Due to the abundant synthesis of 5,6-DMB, it should not act as the limiting factor for DMB B<sub>12</sub><sup>[25]</sup>. <i>Synechococcus elongatus</i> PCC 7942 and <i>Arthrospira platensis</i> C1 both have the <i>cobS</i> gene which codes for a protein that binds the α-ribazole within the 5,6-dimethylbenzimidazolyl nucleotide moiety to the cobalt-containing adenosylcobinamide-GDP complex<sup>[14,15]</sup>. Previous research involving the DMB B<sub>12</sub> pathway suggests that a lower ligand constructed with 5,6-DMB has a 100 fold higher binding affinity for cobalt than that of one constructed with adenine<sup>[16, 17, 10]</sup>. Due to the high binding affinity and abundance of 5,6-DMB as a lower ligand when compared to adenine, we assume a 100:1 ratio for DMB B<sub>12</sub> to adenine B<sub>12</sub>; for every 100 active DMB B<sub>12</sub> synthesized, there will be 1 inactive adenine B<sub>12</sub> synthesized | + | <p>The quantity of DMB B<sub>12</sub> produced depends on a successful production and integration of the active B<sub>12</sub> lower ligand, composed of 5,6-dimethylbenzimidazole (5,6-DMB) and α-ribazole 5’-phosphate. For cyanobacteria in the wild, cobalt is often the limiting factor for growth and production of B<sub>12</sub><sup>[11,12]</sup>. Both <i>ssuE</i> and <i>bluB</i> genes are engineered to be regulated using a strong PtrC promoter, creating an abundance of the activating portion of the lower ligand, 5,6-DMB. Due to the abundant synthesis of 5,6-DMB, it should not act as the limiting factor for DMB B<sub>12</sub><sup>[25]</sup>. <i>Synechococcus elongatus</i> PCC 7942 and <i>Arthrospira platensis</i> C1 both have the <i>cobS</i> gene which codes for a protein that binds the α-ribazole within the 5,6-dimethylbenzimidazolyl nucleotide moiety to the cobalt-containing adenosylcobinamide-GDP complex<sup>[14,15]</sup>. Previous research involving the DMB B<sub>12</sub> pathway suggests that a lower ligand constructed with 5,6-DMB has a 100 fold higher binding affinity for cobalt than that of one constructed with adenine<sup>[16, 17, 10]</sup>. Due to the high binding affinity and abundance of 5,6-DMB as a lower ligand when compared to adenine, we assume a 100:1 ratio for DMB B<sub>12</sub> to adenine B<sub>12</sub>; for every 100 active DMB B<sub>12</sub> synthesized, there will be 1 inactive adenine B<sub>12</sub> synthesized.</p> |
| + | |||
<br> | <br> | ||
<br> | <br> | ||
| Line 877: | Line 878: | ||
<br> | <br> | ||
<br> | <br> | ||
| + | <h2 style="text-align: left; font-weight: 400;">B<sub>12</sub> Conclusion</h2> | ||
| + | <p>Published HPLC results show that <i>Arthrospira platensis</i> produces between 1.5-2.5µg B<sub>12</sub> analogs per gram dry weight<sup>[18]</sup>. Considering the assumed ratio for conversion of B<sub>12</sub> analogs to DMB B<sub>12</sub>, we predict 1.49-2.47µg DMB B<sub>12</sub> per gram dry weight of our engineered <i>A. platensis</i>, thus meeting the USDA’s recommended daily value of 6µg in one 2-4g serving.</p> | ||
| − | |||
<!-- <div class="container"> | <!-- <div class="container"> | ||
<div class="row"> | <div class="row"> | ||
Revision as of 02:31, 2 November 2017
MODELING
Predict and optimize yield
Background
The purpose of modeling is to carefully examine the pathways of each intended biosynthetic product, 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
Acetaminophen Overview
To predict acetaminophen biosynthesis, we analyzed the abundance of the acetaminophen's precursors anthranilate and PABA. Anthranilate comes from chorismate, which is primarily used by the cell to produce the aromatic amino acids phenylalanine, tyrosine, tryptophan, and folate[1]. We used publicly available tryptophan and folate concentrations, estimations from genomic sequences, and ribosomal protein sequences to estimate the amount of precursor available to be made into acetaminophen. Combining this information with enzyme kinetics, we were able to model how much product would be produced per gram of biomass of Arthrospira platensis. Our results indicate that a backyard culture could produce 200 doses of acetaminophen every 10 days.
Method using Published Tryptophan Data
We used published tryptophan and folate concentration data for A. platensis to calculate the moles of the acetaminophen precursors per gram dry mass. Tryptophan and folate precursors anthranilate and PABA are the substrates for our enzyme 4ABH, converting those precursors to 4-aminophenol before nhoA finished the conversion to acetaminophen. We show several different calculations below using different sources of data.
Assumptions
- We assumed an even split of precursor at the branch points of PABA and anthranilate even though enzyme Km’s will affect the amount of precursor going down each pathway.
- We assumed the moles of anthranilate are equal to the moles of tryptophan, because at least as many moles of anthranilate must be present to make that much Tryptophan.
- 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.442\ µmol\ anth}{1g\ biomass}\times\frac{1\ mol\ acet}{3\ moles\ anth}\times\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, which can be used to estimate the abundance of the acetaminophen precursor.
Assumptions
- Using the entire genome sequence to approximate the cellular amino acid composition assumes each codon is expressed in equal concentrations. This does not account for some proteins being less expressed and the resulting skew in amino acid composition.
- Using the ribosomal sequences assumes that the genes encoding ribosomal proteins are representative of the total cellular amino acids ratio.
- 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 precursors.
- Folate is present in such small amounts, at 0.21nM of PABA, that its effect on acetaminophen production is neglegible[7].
Sequence Analysis Details
We ran both the genome and all 55 listed ribosomal protein sequences through custom Python program, which counted the instances of each amino acid in the sequence, totaled the mass, and calculated the ratio of tryptophan to the total mass. This gave a ratio of tryptophan mass relative to total amino acid mass, 0.9% and 0.6% tryptophan for the genome and ribosomal proteins respectively. Multiplying the mass percentage tryptophan by the stated 60% of Arthrospira platensis is protein by mass, we can calculate the grams of tryptophan per gram biomass and convert that to acetaminophen production as previously described.
$$\frac{0.054\ 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}$$
Enzyme Competition
Besides precursor concentrations, the next limiting factor would be how effective our enzyme 4ABH is at converting precursors from the tryptophan and folate pathways into 4-aminophenol. We can either assume all three competing enzymes will have similar precursor affinity and produce all three products in equal quantities or use ratios of each competing enzyme's Km for the limiting intermediates. To make up for unknown enzyme rates and quantities, we made a chorismate metabolic simulation which ultimately was not used because of unreliable Km ratios. For that reason, all of the previous calulations were computed using a one third precursor to acetaminophen conversion rate.
Assumptions
- We assume that our genes were designed, inserted, and translated successfully and in reasonable quantities.
- We assume Km values for Arthrospira platensis are the same as organisms' enzyme Kms. Ortholog sequences were BLAST compared in each case, resulting in a range of similarities between 52% and 63% alignment, suggesting our organism's Km 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.
Enzyme Competition Details
Since few cyanobacterial proteins have been isolated and tested for kinetic data, many of the Kms used are from dramatically different species with different enzymes. Our inserted 4ABH has a Km of 20.4µM for p-amino benzoate(PABA), while the folate 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 Km ratio of 20.4:0.37 (55:1) will represent how many moles of PABA goes to folate versus acetaminophen. For the tryptophan pathway, 4ABH activity for antranilate is assayed as having a rate of 24% of PABA, equal to 34µM for anthranilate. This compares to TrpD's Km for anthranilate of 40nM, meaning we'd have 850:1 tryptophan to acetaminophen product ratio. Since these ratios are based on Kms from completely different organisms, this enzyme kinetics comparison will be left for later troubleshooting while we use a one third precursor conversion rate for simplicity.
Acetaminophen Conclusion
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 tryptophan, we can assume there would be at least that many moles of anthranilate precursor for our inserted pathway to push towards acetaminophen. The three calculations above can be averaged to finally predict 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 3 by 12 foot 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
The quantity of DMB B12 produced depends on a successful production and integration of the active B12 lower ligand, composed of 5,6-dimethylbenzimidazole (5,6-DMB) and α-ribazole 5’-phosphate. For cyanobacteria in the wild, cobalt is often the limiting factor for growth and production of B12[11,12]. Both ssuE and bluB genes are engineered to be regulated using a strong PtrC promoter, creating an abundance of the activating portion of the lower ligand, 5,6-DMB. Due to the abundant synthesis of 5,6-DMB, it should not act as the limiting factor for DMB B12[25]. Synechococcus elongatus PCC 7942 and Arthrospira platensis C1 both have the cobS gene which codes for a protein that binds the α-ribazole within the 5,6-dimethylbenzimidazolyl nucleotide moiety to the cobalt-containing adenosylcobinamide-GDP complex[14,15]. Previous research involving the DMB B12 pathway suggests that a lower ligand constructed with 5,6-DMB has a 100 fold higher binding affinity for cobalt than that of one constructed with adenine[16, 17, 10]. Due to the high binding affinity and abundance of 5,6-DMB as a lower ligand when compared to adenine, we assume a 100:1 ratio for DMB B12 to adenine B12; for every 100 active DMB B12 synthesized, there will be 1 inactive adenine B12 synthesized.
$$\frac{2.5µg\ B_{12}\ analog}{1g\ drymass} * \frac{100\ DMB\ bindings}{100\ DMB\ bindings+1\ adenine\ binding} = \frac{2.47µg\ DMB\ B_{12}}{1g\ biomass}$$
Assumptions
- Cobalt is the primary limiting factor of DMB B12 synthesis[11,12].
- The PtrC promoter is a strong promoter in S. elongatus 7942[25].
- Gene inserts will be expressed, oxidizing FMN to FMNH2, then catalyzing the synthesis of 5,6-DMB in excess[13].
- The cobS protein in A. platensis will attach the α-ribazole part of the 5,6-DMB lower ligand to cobalt-containing adenosylcobinamide-GDP complex 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].
- There will be a 100:1 conversion of inactive adenine B12 analogs to active DMB B12 given the 100 fold higher binding affinity for 5,6-DMB as a lower ligand[16, 17, 10].
B12 Conclusion
Published HPLC results show that Arthrospira platensis produces between 1.5-2.5µg B12 analogs per gram dry weight[18]. Considering the assumed ratio for conversion of B12 analogs to DMB B12, we predict 1.49-2.47µg DMB B12 per gram dry weight of our engineered A. platensis, thus meeting the USDA’s recommended daily value of 6µg in one 2-4g serving.
