Difference between revisions of "BME494s2013 Project Team2"

From OpenWetWare
Jump to: navigation, search
(Testing: Modeling and GFP Imaging)
(Testing: Modeling and GFP Imaging)
Line 283: Line 283:
[[Image:mRNA06.jpg |thumb|noframe|300px|left|Figure 4: mRNA Concentration vs. Time, I = 0.06]]
[[Image:mRNA06.jpg |thumb|noframe|300px|left|Figure 4: mRNA Concentration vs. Time, I = 0.06]]

Revision as of 23:17, 29 April 2013

BME494Sp2013 banner.png
BME494 Asu logo.png

Home        People        Course Projects        Course Materials        Schedule        Photos        Wiki Editing Help


Overview & Purpose

Effective food testing represents a very large portion of the culinary economy.

The food packaging industry is worth $100 billion dollars annually. Among that, about $77 billion dollars represents that meat packaging industry. Despite representing such a high economic value, the food industry still has recalls and contamination scandals every year. Contaminants often include pathogenic bacteria and allergens such as gluten. Furthermore, food waste equates 1.3 billion tons annually. Lack of effective processing results in greater environmental damage through waste and economic damage through business scandals.

Synthetic biology offers a promising solution. A biological assay could be constructed to produce a fluorescent output depending on the presence of various contaminants. This could be done through a modification of the promoter of the GFP to detect a contaminant such as gluten. A quick swab of a food sample could be applied to the bacteria, and if gluten is present, then a green fluorescence could be identified.

Our project is a proof-of-concept of this potential assay through an IPTG inducible modified Lac switch. When exposed to IPTG, a lactose mimic, GFP instead of β-galactosidase is produced, therefore indicating that various input substrates can yield fluorescent outputs. Future applications of the switch could be modified to incorporate a promoter that detects allergens such as gluten or even bacterial contaminants to help mitigate overall health risks.

Background: The Lac Operon

The Lac Operon is a gene specific to E. Coli that controls the cell's digestion of lactose. It consists of a promoter, an operator, three structural genes, and a terminator. It is both positively and negatively regulated, allowing expression to be contingent on the concentrations of glucose and lactose in the cell.

Structure of the Lac Operon [1]

The Lac Operon encodes three structural genes:

  • LacZ: The Lac Z structural region, or β-galactosidase, hydrolyzes the disaccharide lactose into glucose and galactose, sugars that are smaller and easier for the cell to digest. However, in low concentrations of lactose, β-galactosidase cleaves and rearranges lactose into allolactose, which acts as an inducer for the LacI repressor (see Negative Regulation).
  • LacY: LacY, or lactose permease, is a transmembrane protein that transports lactose into the cell.
  • LacA: LacA is a transacetylase. While it has functionality, it has little effect on the output of our design, so it will not be discussed.

In addition to the structural genes, the Lac Operon includes a promoter and an operator region. The promoter region is the area to which the Lac I repressor and the CAP-cAMP complex bind, the mechanics of which will be discussed later (see Positive Regulation and Negative Regulation).

PURPOSE: Efficiency
Expression of the Lac Operon is determined jointly by the levels of glucose and lactose in the cell. Being a monosaccharide, glucose is easier (i.e., takes less energy) to digest; therefore, if glucose is present, the cell will prefer to use it as an energy source. However, if glucose is not available as an energy source, the cell will use lactose instead. A table describing this relationship is below:

Table 1: Glucose/Lactose Relationship to Lac Operon Transcription (2)
Carbohydrates CAP-cAMP Complex LacI Repressor RNA Polymerase Transcription of Lac Operon
+ Glucose, + Lactose Not bound to DNA Lifted off operator site Keeps falling off promoter site Very low transcription
+ Glucose, - Lactose Not bound to DNA Bound to operator site Blocked by repressor No transcription
- Glucose, - Lactose Bound to DNA Bound to operator site Blocked by the repressor No transcription
- Glucose, + Lactose Bound to DNA Lifted off operator site Sits on promoter site TRANSCRIPTION

From this table, we can observe a multitude of things:

  • Transcription only occurs when lactose, but not glucose, is present.
  • When glucose is not present, the CAP-cAMP complex (or the activator protein) is not bound to DNA.
    • Thus, the absence of glucose promotes transcription.
  • When lactose is not present, the LacI repressor is bound to the operator site. When lactose is absent, the repressor is NOT bound to the operator site.
    • Thus, the presence of lactose promotes transcription.
  • For the RNA Polymerase to properly attach to the Lac Operon, the CAP-cAMP complex must be <b>attached to the DNA, and the LacI repressor must not be attached to the operator site.
  • Therefore, transcription only occurs when lactose, but not glucose, is present.

Why does this phenomenon occur? Well, like stated before, lactose is the cell's last resort energy source because it requires more energy from the cell to digest than does glucose. The enzyme that digests lactose is β-galactosidase, which can only be produced by initiating transcription of the Lac Operon. Thus, to be able to digest lactose, the cell needs to initiate transcription of the Lac Operon.


Example of negative regulation of the LacI repressor [5]

The genes encoding the LacI repressor are actually located upstream of the Lac Operon. The LacI gene is not regulated; therefore, it is produced continuously. It binds to the Lac Operon in the promoter region; however, it does not bind if there is lactose in the cell. Why is this? Well, the cell produces very low levels of β-galactosidase even when not in the presence of lactose. In these very low lactose conditions, β-galactosidase has a different function: it cleaves lactose and recombines it to form allolactose, which acts as an inducer for LacI. It binds to LacI and causes a conformational change, which in turn makes LacI unable to bind to the promoter region of the Lac Operon.

Positive and negative regulation of the Lac Operon [4]

Remember from before, the absence of the LacI repressor is not the only factor that allows transcription to occur. There is also a form of positive regulation that occurs via the CAP-cAMP Complex, the formation of which is controlled by the levels of glucose within the cell. As glucose levels in the cell begin to decline, E. Coli responds by beginning to synthesize cyclic adenosine monophosphate, or cAMP. As cAMP concentration increases, it binds to a catabolite activator protein, or CAP. cAMP acts as an inducer for CAP, causing a conformational change that allows CAP to bind to the promoter region of the Lac Operon. This cAMP-CAP complex interacts with the RNA polymerase, increasing its affinity for the Lac promoter. Without attachment of the cAMP-CAP complex, or in high levels of glucose, affinity wouldn't be high enough to cause significant transcription.

A circuit diagram of the natural Lac Operon illustrating glucose and lactose as inputs.

So far, what do we know about the natural Lac Operon? Well, we know that it is a gene that produces a number of structural proteins, including β-galactosidase. We also know that it requires lactose to be present and glucose to be absent for transcription to occur. What we haven't discussed is the function of the operon from an engineering perspective.

Table 2: Glucose NOT Gate
Input (Glucose) Output
1 0
0 1
Table 3: Glucose and Lactose AND Gate
Input 1 (Table 2 Output) Input 2 (Lactose) Output
0 0 0
0 1 0
1 0 0
1 1 1

If we analyze it from a digital logic context, we can describe glucose and lactose as inputs, and the transcription of β-galactosidase as an output. Furthermore, we can build a logic circuit symbolizing the operon's functionality (illustrated in diagram on left). When glucose acts as an input, it produces a NOT gate functionality (See Table 2).

When lactose and the NOT gate output of glucose are incorporated as inputs to the system, they produce an AND gate functionality (see Table 3).

Furthermore, there are a couple of other other proteins that "mimic" the function of lactose as an input for the natural lac operon. Among these are IPTG (used for our switch), and the previously mentioned allolactose which is an isomer of lactose.

From a genetic engineering perspective, the structural genes encoding for β-galactosidase can be replaced with different genes encoding other proteins, therefore creating a designer system that produces proteins when induced by lactose. A few genetic regions that can be substituted are a gene coding for insulin to produce insulin for people with diabetes, a gene for heparin to use as a coating on implant device, or even a gene that codes for a fluorescent protein, therefore producing an indicator of the presence of lactose. Applications of engineering the Lac Operon are endless, only bounded by the human imagination.

Design: Our genetic circuit

As described above, the structural protein regions of the natural Lac Operon can be replaced by various other protein coding regions to alter the output of the Lac Operon. In the case of our gene switch, we chose to replace β-galactosidase with a gene coding for GFP, or green fluorescent protein. We also chose to use IPTG instead of lactose as the system's input; we were able to do this because IPTG acts as a lactose mimic due to its similar structure and active regions. (Note: IPTG is frequently used as a synthetic replacement for lactose in systems created in laboratory settings.) Therefore, our switch turns "on" in the presence of IPTG (our input) and produces a green fluorescent color (GFP) as its output.

Also, we chose a promoter that was not sensitive to the CAP-cAMP complex so that our switch would not be influenced by the presence of glucose. Therefore, IPTG would be the only input affecting our system, as it would not be positively regulated by glucose intake.

Device design. Image adapted from [3].

Our design incorporates BioBrick parts from the Registry of Standard Biological Parts (see diagram on left for specific part numbers). There are two main BioBricks used in our system:

Brick 1: IPTG-Inducible Lac Promoter Brick
The first BioBrick includes the genetic region that codes for the LacI repressor protein, as well a promoter adapted from the natural Lac Operon that is negatively regulated by it. This BioBrick consists of:

  • A consitutive promoter: Causes transcription to begin.
  • A ribosome binding site: Ribosomes will attach here during transcription.
  • The gene for the LacI repressor protein: This gene is what will be transcribed to create the protein of interest, the LacI repressor.
  • Terminators: These signal the end of the transcription process.
  • LacI regulated promoter: This promoter will cause the next stage of transcription to begin, and is negatively regulated (repressed) by the LacI protein. Therefore, transcription will occur only in the absence of the LacI protein.

Brick 2: GFP Production Brick
The second BioBrick includes the parts necessary to produce an output of Green Fluorescence Protein (GFP). This output is regulated by the parts from the previous stage. This BioBrick consists of:

  • A ribosome binding site: This stage involves a second round of transcription, so it needs its own site for ribosome binding.
  • The GFP gene: This is the gene that will be transcribed to produce GFP.
  • Terminators: These signal the end of the transcription process.

The switch response of this device is due to the relationships it creates between IPTG, the LacI protein, and the GFP output. Transcription of the GFP output depends on the activity of the stage 2 promoter, the Lac-I regulated promoter. If this promoter is active, GFP will be produced. This promoter is regulated by the LacI repressor protein. Presence of the LacI protein inhibits the promoter, which turns off GFP production. The LacI protein is created in stage 1 of the genetic circuit. In its default state, the mechanism would operate as follows:

  • LacI protein is created → LacI regulated promoter is inhibited → Transcription of GFP is inhibitedNo Output

On the other hand, when an IPTG input is added to the system, this results in the following:

  • IPTG is added → LacI protein is created → IPTG binds to LacI → LacI can no longer inhibit the GFP promoterGFP Is Produced
Circuit diagram showing input IPTG and output GFP. PLac refers to the promoter for the creation of LacI; PGFP refers to the promoter for the creation of GFP, which is regulated by LacI.

Building: Assembly Scheme


As stated above, the device circuit is made out of 3 parts: The IPTG inducible Lac promoter cassette, the GFP brick and the vector. Listed below, are all the details and specifics for each part.

IPTG Inducible Lac promoter cassette:

  • Includes: promoter for Lac-I, RBS, the Lac-I repressor gene, terminator for the Lac-I transcription, promoter that is Lac-I regulated and IPTG induced
  • BBa_K418000
  • Length: 1416 base pairs

GFP brick (includes RBS, GFP gene, and terminator):

  • BBa_K283026
  • Length: 914 base pairs


  • pSB1A3
  • Length: 2155 base pairs


Type IIS Assembly is used to assemble all the DNA parts to form the final device. Type IIS Assembly is protocol that is used to seamlessly assemble DNA parts through the use of a Type IIS restriction enzyme.

Key features of Type IIS Assembly:

  • Implements "sticky overhangs"
  • Reduces ligation to single reaction
  • Uses the restriction enzyme BsmBI

Theoretical Gel Slide.
PCR Reagents.
Ligation and Digestion Reagents.


Before the parts can be assembled, PCR must be used to amplify the parts and to add the necessary primers. The copied DNA parts are then purified and ligated using the BsmBI/T4 ligase mediated assembly. To ensure a successful ligation step, the DNA parts were screened for BsmBI sites using APE. If any of the DNA parts contain the BsmBI sequence anywhere other then at their restriction sites, a separate primer is made to modify the sequence.


The assembly also requires the design and use of custom primers. The purpose of the primers is to add BsmBI sites to the ends of all the DNA parts. The following custom primers are used for this device:

Forward Lac-I primer: cacaccaCGTCTCaTAGAttgacggctagctc

Reverse Lac-I Primer: cacaccaCGTCTCaTAGAtgagctagccgtcaa

Forward GFP Brick Primer: cacaccaCGTCTCaaaagaggagaaata

Reverse GFP Brick Primer: cacaccaCGTCTCaTAGTtataaacgcagaaag

Forward Vector Primer: cacaccaCGTCTCaactagtagcggccgct

Reverse Vector Primer: cacaccaCGTCTCatctagatgcggccgcg

Testing: Modeling and GFP Imaging


Another example of a (very complex) mathematical model. [6]

We used a previously published synthetic switch, developed by Ceroni et al., to understand how our system could potentially be modeled and simulated. A mathematical model is a mathematical representation of system behaviors defined by the relationships between various system parameters. Parameters are simply different values that affect the behavior of the system. One could even use a simple algebraic equation to represent a mathematical model. In the following equation,

y = 3x - 7

the equation "3x - 7" would be a mathematical model of the system y. Because the value of x affects the ultimate value of y, x would be a parameter of this system.

Flowchart Diagram of Ceroni Data and Respective Values

The Ceroni model diagram is presented left. The Ceroni et. al model took advantage of a dual plasmid primer with one plasmid constitutively producing LacI and another coding for the GFP protein. Included in the model was binding, transcription, translation, and degradation rates. The diagram left is based on the 2011 Northwestern model combined with variables of the Ceroni experiment. The variables are as follows:

  • I - IPTG Concentration
  • G - Concentration of GFP
  • LF - Free, unbound LacI molecules in cell system
  • LI - LacI molecules bound to IPTG
  • MG - mRNA molecules of GFP
  • ML - mRNA molecules of LacI
  • aG - GFP rate of synthesis
  • AL - LacI rate of synthesis
  • aMG - GFP transcription rate
  • aML - LacI transcription rate

We used a model of the natural Lac operon to understand how changing the parameter values changes the behavior of the system. Some of the parameters that were used to describe its behavior are as follows:

  • Mu - Describes the dilution of the system input, or IPTG. A mathematical way of thinking about this would be to take IPTG concentration as a percentage of cell volume, or [IPTG]/cell volume.
  • GammaM - Every protein in a cell has a limited lifespan; at some point, chemical reactions will occur that degrade it or cause it to lose its functionality. GammaM represents the degradation rate of M, M being the concentration of mRNA for Bgal (β-galactosidase) in the cell.
  • GammaB - Represents the degradation rate of Bgal, or β-galactosidase.
  • K - Represents the half-max of the transfer function, or the point at which the output reaches half of its maximum output. It also represents the concentration value of the input at the point where the rate of increase of output is at a maximum.
  • AlphaM - Represents the production rate of M (mRNA), or the rate at which mRNA is transcribed from DNA.
  • AlphaB - Represents the production rate of Bgal (β-galactosidase), or the rate at which β-galactosidase is translated from mRNA.
  • Cell = 1 - "Cell" represents the bacterial cell volume. This statement sets this volume to an arbitrary value 1.

When analyzing this model, we were concerned with three main things: how mRNA concentration changed over time, how IPTG concentration changed over time, and how the concentration of β-galactosidase changed over time. We ran a simulation setting the concentration of IPTG at an arbitrary value of 0.32 and produced the graphs of mRNA concentration over time (Figure 1) and IPTG concentration over time (Figure 2) as shown. The concentration of β-galactosidase over time for this input concentration is the top blue curve in Figure 3.

Figure 1: mRNA Concentration vs. Time
Figure 2: IPTG Concentration vs. Time

As we can see from these graphs, plotting mRNA concentration over time (Figure 1) results in a curve that oscillates around the value of 6.06x10^-4. Plotting IPTG concentration over time (Figure 2) results in a straight line, which is unsurprising considering that the model assumes that IPTG stays constant over time for simplification purposes, which is of course not true in real life. Plotting the concentration of β-galactosidase over time (blue curve in Figure 3) results in a curve that seems to max out at a value of 4.5x10^-4. Logically, this makes sense: at some point in time, the rate of production of β-galactosidase (AlphaB) must reach the degradation rate of β-galactosidase (GammaB), meaning that the concentration of β-galactosidase will no longer increase; it will simply stay the same no matter how much more IPTG is added to the system.

Figure 3: β-galactosidase Concentration vs. Time

To further analyze how IPTG concentration affects the concentration of β-galactosidase, we thought we would change the initial value for the concentration of IPTG. We set the concentration of IPTG to the arbitrary value of 0.25 and produced the second curve in Figure 3.

From this curve we can see that when the IPTG concentration is 0.25, the maximum output of β-galactosidase seems to be about 3.0x10^-4. This value is significantly lower than the value procured from Figure 3; thus, we can infer that as the system input of IPTG concentration decreases, so does the system output of β-galactosidase.

Changing the input parameters and observing the simulated response would allow us to decide what levels of input concentration would be beneficial for Bgal output. Therefore, we varied the parameters in search of a cutoff input concentration level, above which would produce a sustained Bgal output, and below which would produce a declining output over time. After simulating several different initial concentrations, we found that the turning point occurred at about I = 0.063. Figure 4 shows the resultant mRNA concentration in this case. The Bgal concentration over time was also graphed for this input concentration and appears as the bottom red curve in Figure 3. Therefore, transcription and the output Bgal decrease over time for an input concentration of I = 0.06, below the threshold we found.

Figure 4: mRNA Concentration vs. Time, I = 0.06

We explored how one technique, imaging via microscopy could be used to determine the production rate of an output protein, in this case GFP in yeast, could be used to determine a "real" value for maximum GFP production rate under our own laboratory conditions.

Frame 18 of GFP producing cell images.

Images of GFP producing cells were taken over time in order to track GFP output. There were 24 frames of images sorted chronologically. ImageJ software was used to track output data of one cell at a time across the 24 time points. The average intensity of the image across the chosen cell was recorded for each frame, since this brightness would correspond to GFP output level. Data of this nature was recorded for 8 different cells, allowing for the calculation of an average brightness or GFP production level at each frame or time point. This average production over time was then graphed using MatLab. Using MatLab's native fitting tools, a 6th order polynomial was fit to the plot and the equation for this polynomial was generated, as shown below.

6th degree polynomial fit to the raw GFP output data; K value and Hill Coefficient, n, shown.

As described above, fitting this data to a mathematical equation allows us to solve for parameter values of the system, from which we can infer information about production of the output. The Hill Coefficient, n, can be found from this fit. The Hill Coefficient occurs at the inflection point of the rising output curve, meaning it represents the maximum rate of GFP production. The parameter K described above is related to the Hill Coefficient, n. K represents the input concentration at which n occurs. In our mathematical model, K and n are the x and y values at which the inflection point of the GFP production curve occurs. The inflection point of the curve corresponds to the maximum of the derivative of the curve. Therefore, the maximum rate of GFP output production would occur at the maximum of the derivative curve of the polynomial fit equation. To find the maximum rate of GFP production relative to changing input concentration, this derivative was graphed and its maximum was found.

5th degree derivative of the polynomial fit; maximum GFP production rate shown.

Ideally, the GFP production rate measured by this method could be entered as a value for the GFP transcription rate, aMG in the Ceroni et al. model.

Human Practices

Gap Analysis of the GFP Assay

Food Waste
1.3 billion tons of food thrown away yearly. This creates massive wastes ranging from environmental problems to simple affluence. Use of a faster colorimetric assay would be effective in mitigating the overall waste generated by food. This problem goes beyond economic incentives for the businesses. It also supports the consumer--individuals who are exposed to various contaminants and allergens that should not be found otherwise. The impact of reducing food waste would encourage a healthier ecosystem of effective resource stewardship, less landfills, and less health complications on both corporate and social levels.

Our Team

Shay Ravacchioli

  • My name is Shay Ravacchioli, and I am a Junior majoring in Biomedical Engineering with minors in Biological Sciences and Psychology. I am taking BME 494 because I think Synthetic Biology is fascinating. An interesting fact about me is that I play piano and guitar.

Jenessa Lancaster

  • My name is Jenessa Lancaster, and I am a Junior majoring in Biomedical Engineering with a minor in Psychology. I am taking BME 494 because I have always wanted to learn more about Synthetic Biology and Genetic Engineering. An interesting fact about me is that I write songs.

Michael Rose

  • My name is Michael Rose, and I am a senior majoring in BME. I am taking BME 494 because I find synthetic biology very interesting and I want to learn more about it. An interesting fact about me is that I love to cook.

Kwanho Yun

  • My name is Kwanho Yun, and I am a Junior majoring in Molecular Biosciences and Biotechnology. I am taking BME 494 because I am a big fan of synthetic biology and hope to learn more about the industry and how the technology works. An interesting fact about me is that I grow red beard hairs.

Works Cited

[1] Heller, H. Craig., David M. Hillis, Gordon H. Orians, William K. Purves, and David Sadava. Life: The Science of Biology. Sunderland, MA,: Sinauer Ass., W.H. Freeman and, 2008. N. pag. Print.

[2] Escalante, Ananias. "Regulation I." Class Notes. University of Arizona. 20 February 2013.

[3] Registry of Standard Biological Parts. Web. 25 Apr 2013. <partsregistry.org>.

[4] Slonczewski, Joan, and John Watkins. Foster. Microbiology: An Evolving Science. New York: W.W. Norton &, 2009. Print.

[5] Schmidt, Markus. Synthetic Biology: Industrial and Environmental Applications. Weinheim, Germany: Wiley-Blackwell, 2012. Print.

[6] Filho, Ernesto R., Fransisco L. Nunes, Jr., and Sidney O. Nunes. "Synchronus Machine Field Current Calculation Taking Into Account the Magnetic Saturation." SciELO - Scientific Electronic Library Online. N.p., May 2002. Web. 26 Apr. 2013.

[7] "Haynes:TypeIIS Assembly." OpenWetWare, . 25 Feb 2013, 18:20 UTC. 18 Mar 2013, 03:49 <http://openwetware.org/index.php?title=Haynes:TypeIIS_Assembly&oldid=679110>.

[8] Engler C, Gruetzner R, Kandzia R, Marillonnet S. (2009) Golden Gate Shuffling: A One-Pot DNA Shuffling Method Based on Type IIs Restriction Enzymes. PLoS One. 4