Team:XMU-China/Content3

Hi, Oscillation!

Are you familiar with oscillation?

No? Oh, you must be kidding!

Let me present you its definition in Wikipedia to help you get a clue. Oscillation is the repetitive variation, typically in time, of some measure about a central value (often a point of equilibrium) or between two or more different states^[1] .

It seems that you have met it somewhere, right?

Have you ever used a GPS? Or lasers? Alright, you must have used AC power that supports nearly all electrical appliances in your houses. These three are typical applications of oscillation and amplified that signals in the frequency domain have obvious advantages over those steady-state design in terms of information gathering and procession. Take lasers for example, which are known for their intensity and can be focused to a tight spot over long distance. These characteristics all owe to their spatial coherent in the frequency of the light source.

Oscillation in bacteria

Scientists have proved oscillations also pervade biological systems at all scales as well, from gene expression to cell cycle progression, and these oscillations can incorporate the periodic variation in a parameter over time to generate an oscillatory output^[2]. As mentioned above, oscillations can lead to fantastic applications and benefit our everyday life. Since the output of bacteria oscillations can be detected as a frequency, people started their journey on designing a bacteria reporter.

The noise of gene expression in biological systems, however, slowed our pace in this process, because it will generate noisy or stochastic oscillation with varying amplitudes and frequencies. To deal with this problem, we have to unify the expression of the reporter gene in bacteria. And we found…

First generation of oscillation, the second and the third…Brief intro

Synchronized Oscillation

Yes, we found synchronized oscillation!

It was in 1670s that Christiaan Huygens first observed coupled oscillations: two of his pendulum clocks mounted next to each other on the same support often became synchronized3. To get cells communicate and oscillate in the same amplitude and frequency like the two pendulum clocks, we adopted quorum sensing, a cell-to-cell signaling mechanism that refers to the ability of bacteria to respond to chemical hormone-like molecules called autoinducers, from Vibrio fischeri into our host E. coli (MG1655) using our “hero” synthetic biology to realize synchronization among cells over tens of micrometers4.

QS now or later

As far as we concern, more consistent the oscillation can be if more colonies are synchronized. In order to enhance the communication range of bacteria, a gas-phase redox signal molecular H₂O₂ is introduced into our circuit. According to a published research5, through H₂O₂ a faster and long distance instantaneous communication can be achieved to strengthen the oscillation.

For more detailed information about our circuit, please refer to the Mechanism part of our project.

Introduction

In a synchronized oscillatory system, three important parts should be included: the oscillator, which is the biochemical machinery that generate the oscillatory output; the coupling pathway that ensure the connection among cells; and output pathway, which is also known as a reporter that reflect the state of the oscillator to downstream targets.

Oscillator

Quorum sensing (QS) is a cell-to-cell signaling mechanism that refers to the ability of bacteria to respond to chemical hormone-like molecules called autoinducers. When an autoinducer reaches a critical threshold, the bacteria detect and respond to this signal by altering their gene expression. In our circuit, QS (from Vibrio fischeri) is installed as a positive feedback while aiiA (from Bacillus Thurigensis) acts as a negative one to compose an oscillator together. In the quorum sensing part, the luxIgene is at low expression level and produces LuxI protein that synthesize a kind of acyl-homoserine lactone (AHL), which is a small molecule that can diffuse across the cell membrane and mediate intercellular coupling when it reaches the threshold as enough biomass accumulated. AHL will bind intracellular protein LuxR, which is also consecutively produced by luxR gene. The LuxR-AHL complex can activate the luxIpromoter, and the positive feedback loop is built. At the same time, the aiiA gene, which is under control of the same promoter is expressed and produce a lactonase enzyme known as aiiA that hydrolyzes the lactone ring of AHL. (Fig. 1) In this system, the activator enhances the expression of both activator and repressor, which shares the common motif of many synthetic oscillators.

Coupling

The communication range of quorum sensing, however, is limited by the diffusion rate of AHL, could only reach cells over tens of micrometers. So, we introduced the coding gene for NADH dehydrogenase II (ndh) and put it under the control of the same luxIpromoter, which means ndh also has a periodical expression in accordance with the oscillator. NDH-2 is a membrane-bound respiratory enzyme that generates low level H₂O₂ and superoxide (O₂^-) and H₂O₂ will permeate to neighboring colonies. Periodic production of H₂O₂ changes the redox state of a cell immediately and interacts with the synthetic circuit through the native aerobic response control systems, including ArcAB, which has a binding site in the lux promoter region. Before the oxidizing condition is triggered, ArcAB is partially expressed, so lux is partially repressed. When the concentration of H₂O₂ is increasing, ArcAB is gradually inactivated, and the repression of lux is relieved. With the periodically produced vapor phase H₂O₂ that can diffuse quickly among colonies, the oscillation is not only strengthened but also expanded to millimeter scales.

When compared with the strong but short ranged coupling by QS, H₂O₂ might be weaker but long ranged because its disperse characteristic. These two levels of communication between cells formed the basic oscillatory system inside our host. [the host choosing]

Reporter

A good reporter in an oscillatory circuit should be steadily expressed and can be easily detected by regular equipment in laboratory, so broadly used reporter green fluorescent protein (GFP) came into our mind.

First, we built regular gfp gene into our circuit and it works well. Its fluorescence can be observed by fluorescence microscopy under 1s of exposure time on the microfluidic array 1.

However, there is a newly found GFP variant Superfolder GFP, which has shown outstanding performance in many ways. It has two advantages in observing:

1) Fold 3.5 to 4 times faster than regular GFP;

2) Yield up to four times more total Fluorescence than regular GFP.

The greater the fluorescence strength is, the shorter exposure time will be needed, thus can decrease the photobleaching (GFP will be eventually destroyed by the light exposure that tries to stimulate it into fluorescing) that will lead into lapses in data processing.

Besides, we also want to make a comparison between two different GFPs’ performances in our oscillatory circuit. For example, peaks, troughs and periods in oscillatory curves.

[Data/Pictures] P. S. Thanks the Peking iGEM team for providing us the sfGFP part.

Host

After building the whole circuit in DH5α, we have to find a suitable host for it. It is as challenging as casting for a right leading actor for our blockbuster Synchronized Oscillatory.

At first, we chose BL21 (DE3) to be this lucky guy for the following reasons:
1) ndh gene is originally expressed in BL21 (DE3), and that’s where we get this part by PCR. So we think BL21 must have a full developed ndh expression system that can be a credit for the normal function our coupling part.

2) BL21 is always known for its competent in transformation and protein expression among a variety of E. coli cells.

However, BL21with three plasmids let us down after some characterizations.

The DNA gel electrophoresis for digestion confirmation results showed that there is something wrong with the gfp plasmid.

The single digestion by EcoR I (Band 2) should only generate a single band which runs a little faster than the plasmid itself (Band 1). From the picture above, however, an unknown band is found (Circled in red), and we couldn’t find a solution for it.

To be continued…

Note:

Quorum sensing (QS) is a cell-to-cell signaling mechanism that refers to the ability of bacteria to respond to chemical hormone-like molecules called autoinducers. When an autoinducer reaches a critical threshold, the bacteria detect and respond to this signal by altering their gene expression.

The autoinducer in our circuit is acyl-homoserine lactone (AHL).

http://en.wikipedia.org/wiki/Photobleaching

Photobleaching is the photochemical destruction of a dye or a fluorophore. In microscopy, photobleaching may complicate the observation of fluorescent molecules, since they will eventually be destroyed by the light exposure necessary to stimulate them into fluorescing. This is especially problematic in time-lapse microscopy.

Loss of activity caused by photobleaching can be controlled by reducing the intensity or time-span of light exposure, by increasing the concentration of fluorophores, by reducing the frequency and thus the photon energy of the input light, or by employing more robust fluorophores that are less prone to bleaching (e.g. Alexa Fluors or DyLight Fluors). To a reasonable approximation, a given molecule will be destroyed after a constant exposure (intensity of emission X emission time X number of cycles) because, in a constant environment, each absorption-emission cycle has an equal probability of causing photobleaching.

We designed three plasmids to function as the three important pA rts above according to a published research, and all of them all in the charge of the same promoter – lux promoter to make sure that all target genes share the same oscillatory period.

The first plasmid is the pSB1C3-gfp-luxI (pG-L) plasmid, which will express GFP reporter and LuxI proteins; the second plasmid is pSB3T5-aiiA (pA ), which will express protein aiiA to degrade AHL; the third plasmid is pSB4K5-ndh (pN ), which expresses ndh to generate H₂O₂ to communicate between colonies.

You may have noticed that these three plasmids have different resistant genes and replication origins on their backbone. Different resistant genes are used mainly for …. A different replication origin means a different copy number of target genes in this plasmid. Using different replication origins is to produce different target proteins in an appropriate ratio, and only in this way can oscillation be observed.

At first, we assume that the relationship between pG-L, pA and pN should be high-middle-low. [reason] So we chose pMB1, p15A and pSC101 because they seem to meet the requirement of copy number. However, the oscillation we observed on our microfluidic array is not …, because the light of green fluorescence protein is not dim enough at the trough. (We also discussed this problem with the author of this article, and his suggestion was to lower the copy number of the pG-L). Then we thought the explanation for this phenomenon is there was too much GFP but not enough aiiA to degrade them all in time, so GFP remained longer than we expected and weak the oscillatory effect. So we changed the backbone of plasmid pSB1C3--gfp-luxI into … with a … replication origin, and the copy number is ….

Why M&M?

If we want to capture gene expression dynamics and distributions reported by fluorescence, we need to create a chemostat that provide a near-constant environment for the growth of bacteria. The near-constant environment mainly includes two aspects: continuous fresh media to ensure the optimum growth state and constant small populations of bacteria to ensure single-cell resolution in observation. Steady growth state is needed because the oscillation curve is expected to be a near simple harmonic one. Besides, bacteria tend to form biofilm which will lead to the overlapping of cells and create a false impression that the fluorescence increases, and that’s why we have to prevent this to maintain a constant bacteria population with single-cell resolution.

What we mentioned above might have covered all aspects we need to consider for communication in a single colony, but if we still want to observe the coupling between colonies we have to create an array that contains many colonies and allows H₂O₂ to pass through freely.

Among many different techniques, microfluidic array seems to be our best choice with following advantages:

1) Precisely control of flowing rate and cell populations when equipped with syringe pumps;
2) Cells behavior including fluorescence strength can be captured by high resolution microscopy;
3) Trap numbers and dimensions can be adapted to different experimental need;
4) Productive: fewer resources and less time are required when compared with macroscopic cultures.

Design of Microfluidic Array

According to what we have discussed above, the microfluidic device should be able to provide flowing fresh media and maintain a constant bacteria population with single-cell resolution. So we designed our device as shown below.

About devices:

a) Material: PDMS

b) Incubation temperature: around 30℃

c) *Tween20 in media: 0.075%

d) Trap size: ？μm height 100*100 μm space between traps was 25μm channel wide 100 μm

e) Flow rate: around 0.01 ml/hr, 3mm radius, for LB

will be increased according to the cell density

This device is composed by two parts: the flowing channels and near 3,000 square trapping chambers, and their layout is like there are many passages with symmetry rooms distributed on both sides, so we call this array “passage” model. Firstly, cells are loaded from the cell port. When the loading stops, most cells would remain in channels, but if the loading cell density is dense enough, many cells would be pushed into chambers and form a monolayer of Escherichia coli growing in it. This is the time for running LB to get in. Fresh LB is loaded from the media port with Tween 20* and appropriate antibiotics. Flowing LB does not only permit a constant nutrition and exponentially growing cells for more than two days (that’s the longest time we have tried), but also washes away newly generated cells to sustain a certain cell density, which is crucial for the determination of fluorescence strength.

It has been stated before in the Mechanism part that most quorum sensing systems have a critical cell density requirement, so the trap size is designed neither too big nor too small to capture enough cells and allows for nutrients to permeate to every corner of the chamber.

You may also notice that the height of the chamber is extremely small, which is really demanding for the researcher’s skills, but it’s important for the observation. The oscillation could only be seen when the height of the chamber is in a particular range as our modeling predicted for many complicated reasons. For more details, please refer to the Modeling part.

Polydimethylsiloxane (PDMS), a transparent polymer material, is chosen to build the devise because of its excellent breathability, which is necessary for the permeation of H₂O₂. Once H₂O₂ diffused outside cells, it can travel from colony to colony and synchronize them.

Our device is not only environmental friendly, but also a reusable one. After the flushing of a lysis buffer, most cells can be washed away, even those attached to walls. Then the device will be washed for three times by sterilized water before reuse.

The green fluorescence of bacteria in chambers is captured every 5 minutes by an inverted microscope (Eclipse Ti, Nikon) with a standard CCD camera manually and usually last more than 7 hrs. (This is really a big challenge for our team and very physically demanding since this is the only available fluorescence we could find. However, we managed to finish this data collecting process anyway. Good jobs, guys！)

Data are transferred into oscillation curves by ImageJ and software designed by our teammate for digital image batch processing. (More detail for our self-designed software, please refer to protocol.)

Experiments with different models

Before our experiments on the microfluidic array we described before, a series of pretests were done on another kind of chip.

“Swimming pool” chip in black field, those white squares (100μm×100μm×100μm) are chambers for bacteria

From the sketch we can see that the chamber is like a swimming pool where bacteria swim happily in the fresh LB that running from up above, and newly generated cells would be washed away.

A “swimming pool” chamber model would have a lager AHL concentration when compared with a “passage” chamber because the dilution effect is minor. According to a published research that the oscillation period will be shorter as the AHL concentration increases. So we assume that it will be more time-consuming to use this “swimming pool” model in pretests.

However, this chamber is about 100 μm deep and could grow E. coli far more than monolayer. As the time goes by, the oscillation in a single colony signaling by AHL could never go back to the original baseline, and the oscillation curves have a tendency like the following graph.

In spite of the effect of growing bacteria, a slight oscillation was still observed. However, this data was not solid enough since we only captured one period of oscillation because of the limited microscope using time.

The oscillation period is about 120 minutes with a tendency of increasing fluorescence strength. However, it still proved that our circuit could function as an oscillation.

BL21 strain with three plasmids (pSB1C3-gfp-luxI & pSB3T5-aiiA & pSB4K5-ndh) has also been tried several times on this array, but no obvious phenomenon could be captured. There are many factors involved, for example cell density, trap size and flowing rate. The strict requirement for the observation of oscillation has already discussed before.

DH5α strains have been tested on this array, too. They didn’t work well, either. The reasons could be found in Host part of our project.

Mathematical model

We modeled the oscillate gene circuit (without ndh gene) based on the paper “A synchronized quorum of genetic clocks”^[12]. In that, the author characterized the oscillator by using the four delay differential equations (DDEs), which are the basic of our modeling. DDE can perfectly simulate our circuit due to its powerful function about solving the dynamic change in biological processes.

In the first two equations, item is a Hill function equaling to , which describes the delayed reactions about the complex processes in the cells (transcription, translation, maturation, etc.). This Hill function takes the history of the system into account, i.e. the concentration of AHL at the time it binds to LuxR to form the activation complex.
This model does not include an equation for LuxR assuming that it is constitutively produced at a constant level. All parameters are normalized for the convenience of analysis.

Modeling in MATLAB

We used MATLAB to help us solve the DDEs and draw the graphs. In addition, we developed a GUI user interface after finishing the main DDE code. The user interface makes people easier to modify and adjust the parameters without recoding the .m files and any MATLAB language knowledge. We called this software “Gene OS” which is short for Gene oscillate simulator.

After travelling a long way of circuit construction, you may wonder what we are going to do with those sparkling bacteria. Of course much thought was given, too, to make our oscillation circuit both showy and substantial as a safety guard in daily life. Here we are going to show you the potential of those shinny soldiers as biosensor, cancer cell killer and gene sinusoidal signal generator.

Biosensor

To date there appears a great number of biosensors that are used for detecting certain factors, such as small organic molecules or metal ions. And we’ve noticed that each year many iGEM teams relate their projects to biosensors. Biological detection has strength in its low detection limit, high sensitivity and bioaffinity when compared with chemical detection methods. We can use the same gene sensing circuit to detect various metal ions or other determinant via changing different promoters which responds to a specific target molecule.

Compared with traditional single-cell biosensor, our oscillatory system biosensor has three advantages.

1) Synchronized colonies will be more sensitive to concentration changes of the target molecule and the changes in fluorescence strength are easier to observe than single cell biosensors.

2) Single cell biosensor usually works, or “switches on” when the concentration of the specific molecules which can be responded by the promoter reaches to the lowest value. Sometimes this kind of detection is limited when we want to narrow the range of detection. Because the oscillation will happen only when concentration of certain molecule is between two value, the maximum and the minimum, it can dwindle the detection range of the determinant.

3) Oscillatory system is more robust and stable when used as a biosensor. These advantages will promote oscillatory gene circuits a better use in biosensing.

Jeff Hasty from UCSD used the same kind of oscillate circuit to engineer two kinds of arsenic-sensing macroscopic biosensors. These circuits are rewired by a native arsenite-responsive promoter, which is repressed by ArsR in the absence of arsenite. Due to the different mechanism, the sensor outputs are not the same. Figure 1 shows two constructed sensing modules and Figure 2 compares the outputs of two sensors^[1].

Figure 1.Two different constructed sensing modules.

Figure 2. Two different outputs of corresponding biosensors.

Due to the accuracy and specific condition in the field of bio-medicine, the oscillatory biosensor needs to meet the requirement such as synchronization. For example, all cells expected to twinkle at the same time only when cytotoxin existed. In addition, the period or amplitude of the oscillation needs to indicate the concentration of the cytotoxin. Our circuit has the potential to accomplish this kind of job. By further exploration in the relationship between oscillatory parameters such as period and amplitude with other environmental factors such as concentration of antibiotic or cytotoxin, we may create a oscillatory biosensor which can sense different kinds of biological toxin.

Cancer cell killer

Nowadays cancer becomes one of the greatest enemy of our human being. Our main methods used to diagnosis of cancer are based on tissue and organ level, such as MRI, CT, PET, etc. These are mainly physical diagnostic methods.

With the development of science and technology, we create various methods into molecule level concentrated on molecular recognition and signal transduction of tumor marker, such as AFP(alpha fetoprotein), CEA (carcino-embryonic antigen), PSA(Prostate Specific Antigen) and so on. Level of these special molecules in human bodies has become the main basis of diagnosis of cancer. These are mainly chemical diagnostic methods.

So how about using synthetic biology to do something for diagnosis of cancer? Our project can be a potential biological method to assist diagnosis of cancer.

ROS is reactive oxygen species in cells, which has been proved to be associated with

several neurological disorders including Alzheimer's disease (AD) and Parkinson's diseases^[2][3][4]. Concentration of ROS indicates the level of cell metabolism. Normally the stronger cell’ metabolism is, the more ROS it will produce. In the recent decade scientist found that the occurrence of cancer is closely related to the level of ROS. But it is not a simple positive correlation as we expect^[5].

From this picture we can get to know it clearly. As we know, we use oxidative stress to describe the level of ROS in our cell. There is a threshold of oxidative stress which is the critical value between normal cell and cancer cell. As a normal cell, its oxidative stress is lower than the critical value while the ROS production value is over the threshold. On the one hand, when the oxidative stress in our internal environment increases, normal cells may mutate to cancer cells if the ROS value is higher than the threshold value. On the other hand, cancer cells may die because of the accumulation of ROS. So the key point is, when ROS value higher than the maximum, we can kill the cancer cell while ROS value lower than the threshold we may push normal cells to change into cancer cells. So the level of ROS is in a specific range which we can measure using chemical methods.

Think about a proper range of ROS’s concentration. Is it familiar to you as we mentioned a similar problem when talking about advantages of oscillatory system as biosensors? H₂O₂ is a kind of ROS. The oscillation will happen only when the concentration of is H₂O₂ between the maximum and the minimum. By regulating our circuits to let it respond in the proper concentration range which matches the ROS range of cancer cells, we can use the oscillatory system as a biological cancer detector.

In our experiments, we found that we can regulate flow rate and some parameters of microfluidic chip such as the depth of the well, the width of the channel and the concentration of H₂O₂ to change the period and amplitude of the oscillation. So by changing parameters of microfluidic chip, we can get a proper concentration range of H₂O₂ if the oscillation works. From then on, we can extract some cancer tissue to detect the ROS concentration in the microfluidic chips or get some suspected tissue on chips and diagnose whether it contains cancer cell.

2) Periodical Killer

2.1 Drug delivery

Our circuit also has a great potential application in periodic drug delivery. It is amazing when we realize that the oscillate circuit can control the drug to release in a specific period and the concentration of the drug is related to the amplitude. It is well-known that nowadays drug-release is a very hot research field. Many efforts are focusing on how to modify or delivery drug using materials such as synthetic polymers, liposome or nano shell which can respond to differences between different intracellular environments. Our oscillatory gene circuits provide a more biofriendly way for drug release. Our future plan is to make some deep exploration in this field.

2.2 Protein Synthesis

What else, a resonance may happen among the oscillators which results in a peak synthesis of protein under some conditions. This phenomenon can be used to increase the yield of some medical proteins such as antitoxic serum and antibody^[9][10].

2.3 Prevention pathogen

In the applications of controlling pathogens, quorum sensing system plays an important role in regulating growth and reputation of microorganism. So we could restrain quorum sensing system to control microorganism to make our life more healthy and nice. The most important part in the metabolism of pathogens controlled by quorum sensing system is the synthesis of toxin and the formation of membrane system. If we can restrain these two kinds of effect may we decrease some potential diseases resources and improve the effects of antitoxins. What’s more, it is a particularly attractive point to prevent pathogens through restraining their quorum sensing system instead of killing them directly as it offers the prospect of combating effectively infections resistant to all existing drugs. As we all know, just treat diseases via using antibotics or antitoxin can help these pathogens improve their resistance to various drugs, which is very dangerous for us human being. To date, although we’ve found that different fungus form different signal molecules, the quorum sensing system almost takes part in cell morphological transformation and information communication. We think our oscillatory gene circuits provide a very useful and meaningful model to promote more and more research in this field^[11].

Gene sinusoidal signal generator

With the development of science and technology, computer had been here and there with larger RAM, faster CPU and smaller volume. All these changes owe to VLSI （Very Large Scale Integrated Circuites）, which consists of different electronic component.

Just image what if we can control cell as a biocomputer. Cellular computer is a device consists of different kinds genetically engineered bacteria, which can do the calculation, such as add and subtract like a normal computer. To finish such a device, scientists need various gene circuits to complete different operations. These circuits are just like electronic component working in an electronic computer.

At the beginning of this century, biocomputer is still a crazy dream far more than a possible blueprint. But now scientists have had the technology to make this “impossible mission” comes true. In recent months, scientists from Stanford constructed a new “bio-crystal valve” by DNA and RNA, which can do the logical operations in the living cells^[6].

Two months later, scientists from MIT used less than 3 gene devices to do the digital operation in division, logarithmic and square root. Thus, constructing gene circuits as electronic devices have become possible. From their research, their cellular calculator which uses analog computation has more advantages in computing efficiency in part count, speed and energy consumption compared with digital computation^[7].

Figure: Different kinds of log-domain analog computation.

Due to the fact that our circuits can generate a steady oscillation, it is possible that using our circuits as a sinusoidal signal generator.

Signal generators, such as function generator (including sinusoidal signal generator), are electronic devices that generate repeating or non-repeating electronic signals. As an important part of basic circuit studies, they are generally used in designing, testing, troubleshooting, and repairing electronic or electroacoustic devices.

What’s more, Li Guang-bi from TUST reported that they have constructed two kinds of electronic circuits which well simulated the gene oscillator. This indicates that our circuit can develop to a new kind of gene sinusoidal signal generator^[8].

Figure 1.This picture shows a kind of electronic circuit generating the oscillation wave.

Figure 2.This graph shows the generated oscillation wave.

Future Version Trilogy

We blueprint our next step developing oscillatory gene circuits in following steps:

Step One

Complete our project by cleaning every obstacle on our way to oscillation.
Confirm that every part is functioning as expected;
Amplify the effect of different copy numbers;
Complete the comparison of different strains’ effect on oscillation;
Complete the comparison of different EGFP & sfGFP effect on oscillation;
Confirm most suitable microfluidic parameters for oscillation;
Optimize our model

Step Two

We will make our oscillation a controllable one by integrating the catalase gene into the circuit.

Catalase is a common enzyme found in nearly all living organisms exposed to oxygen. It catalyzes the decomposition of hydrogen peroxide to water and oxygen. It is a very important enzyme in protecting the cell from oxidative damage by reactive oxygen species (ROS). Since our lab is studying Shewanella oneidensis, and its catalase protein has showed outstanding performance in catalyzing H₂O₂. So we thought we could integrate the catalase gene into our circuit and under an artificial inducible promoter, by which we can decide when to start this catalase gene and stop the oscillation.

Step Three

Realize everything we imagined in APPLICATION!

Make our circuit function as biosensor, cancer cell killer and gene sinusoidal signal generator.

References:
1. Prindle, A., et al., A sensing array of radically coupled genetic 'biopixels'. Nature 481, 39-44 (2012);
2. Gaggelli, E. Kozlowski, H. Valensin, D. Valensin & G. Copper ,Homeostasis and Neurodegenerative Disorders (Alzheimer's,Prion, and Parkinson's Diseases and Amyotrophic Lateral Sclerosis). Chem. Rev. 106, 1995–2044 (2006);
3. Scott, L. E. Orvig, C. Medicinal Inorganic Chemistry Approaches to Passivation and Removal of Aberrant Metal Ions in Disease. Chem. Rev. 109, 4885–4910 (2009);
4. Zheng, Z. Q. White, C. Lee, J. Peterson, T. S. Bush, A. I. Sun, G. Y. Weisman, G. A. & Petris, M. J. Altered Microglial Copper Homeostasis in a Mouse Model of Alzheimer's Disease. J. Neurochem. 114, 1630–1638 (2010);
5. Dunyaporn T., Jerome A. & Peng H., Targeting cancer cells by ROS-mediated mechanisms: a radical therapeutic approach? Nature Reviews Drug Discovery 8, 579-591 (2009);
6. Bonnet J, Yin P, Ortiz ME, Subsoontorn P & Endy D, Amplifying Genetic Logic Gates. Science 340, 599-603 ,2013;
7. Ramiz D., Timothy K. Lu et al. Synthetic analog computation in living cells. Nature 497, 619-+ (2013);
8. Guang bi Li et al, Electronic Ciircuit Model Construction of A Synthetic Gene Network. Journal of Tianjin University of Science & Technology 26, 22-25 (2011);
9. Hasty, J., Dolnik, M., Rottschafer, V. & Collins, J. J. Synthetic gene network for entraining and amplifying cellular oscillations. Physical Review Letters 88, 148101, (2002);
10. Elowitz, M. B. & Leibler, S. A. Synthetic oscillatory network of transcriptional regulators. Nature 403, 335 (2000);
11. Yan C., Quorum--Sensing System in Microbes and Its Application,School of Resources and Environment，Yuxi Normal University，Yuxi，Yunnan 653100.
12. Danino, T., Mondragon-Palomino, O., Tsimring, L. & Hasty, J., A synchronized quorum of genetic clocks. Nature 463, 326-330 (2010)

First

We designed three plasmids to function as the three important pA rts above according to a published research, and all of them all in the charge of the same promoter – lux promoter to make sure that all target genes share the same oscillatory period.

The first plasmid is the pSB1C3-gfp-luxI (pG-L) plasmid, which will express GFP reporter and LuxI proteins; the second plasmid is pSB3T5-aiiA (pA ), which will express protein aiiA to degrade AHL; the third plasmid is pSB4K5-ndh (pN ), which expresses ndh to generate H₂O₂ to communicate between colonies.

You may have noticed that these three plasmids have different resistant genes and replication origins on their backbone. Different resistant genes are used mainly for …. A different replication origin means a different copy number of target genes in this plasmid. Using different replication origins is to produce different target proteins in an appropriate ratio, and only in this way can oscillation be observed.

At first, we assume that the relationship between pG-L, pA and pN should be high-middle-low. [reason] So we chose pMB1, p15A and pSC101 because they seem to meet the requirement of copy number. However, the oscillation we observed on our microfluidic array is not …, because the light of green fluorescence protein is not dim enough at the trough. (We also discussed this problem with the author of this article, and his suggestion was to lower the copy number of the pG-L). Then we thought the explanation for this phenomenon is there was too much GFP but not enough aiiA to degrade them all in time, so GFP remained longer than we expected and weak the oscillatory effect. So we changed the backbone of plasmid pSB1C3--gfp-luxI into … with a … replication origin, and the copy number is ….