Introduction
Several circuits have been proposed, constructed, and implemented, leading to landmark discoveries in synthetic biology. These include systems such as the bi-stable toggle switch, and the repressilator, which brought about a paradigm shift in the field. Since then, several systems have been constructed to employ memory modules, create counters, adders, digital biosensors, and a whole wide range of other products.
Figure – A brief timeline of major notable events in the creation and development of synthetic biology (Source: Del Vecchio, Domitilla et al, Journal of The Royal Society Interface 13.120 (2016): 20160380.)
However, there are several limitations that still need to be overcome, as the field continues to make strides in every area. These involve the fact that biological systems have a lot of noises that cannot be modeled accurately to date, and the fact that metabolic burden is a major issue. Along these lines, one of the central issues is the distinct lack of digital responses in synthetic biology.
Thus, as described in the project overview, we wished to use the high cooperativity TetR homologs in such a manner so as to generate a square wave oscillator circuit. Such a system could have a whole multitude of applications; some of which were mentioned briefly in the overview, and the same are also discussed below,
Thus, as described in the project overview, we wished to use the high cooperativity TetR homologs in such a manner so as to generate a square wave oscillator circuit. Such a system could have a whole multitude of applications; some of which were mentioned briefly in the overview, and the same are also discussed below,
Applications
The potential of genetic clock lies in its role to triggering logic reaction for sequential biological circuits. A square wave generator could be used as a genetic clock, since square waves lie at the heart of clocks. Further, these clocks could be used in any cellular system to time particular events.
This could be used to study correlation between two genes, by coupling one of the genes to the oscillator, then observe the dynamics of the second gene. In this manner, the effect of one gene on the others could be studied.
Like diabetic patients, insulin needs to be provided externally through injections or other means. Our oscillator system could be used as a patch containing bacteria that are oscillating to produce square levels of insulin to the patient once. This then would deliver insulin automatically at intervals, guided by our system .Further when the requirement for insulin would be high, the amount of insulin being delivered to them could be changed through some source which could dive a change in frequency of the oscillations (future applications could be focused on engineering frequency modulation).
A bacterial species could be engineered to produce levels of the permease for a particular sugar in the form of a square wave. Thus, at varying intervals of time, the permease (say lac permease) would be expressed, which would cause the bacteria to start metabolizing lactose. When it goes off, the bacteria would not express the lac permease and consume glucose. In this manner, we could tune the frequency of the oscillations to ensure metabolic switching and activation of pathway shunts in the manner that we want.
If we can control the frequency of the oscillations and modulate the time spent by the wave in the ON and OFF state (which is basically a function of the frequency itself), we could generate a combination of 0’s and 1’s, in order to generate a bar code. This bar code would be read in time, and therefore would be a temporal bar code. A typical example of this would be to encode the word “iGEM” by 11001001, where 1 represents an ON state for say, 20 minutes. Thus, an 11 response, which represents the letter i, would then be read if the fluorescence stays on for 40 minutes. In this fashion, our device could be used for encryption of data.
The system could be used as the bacterial analogue of memory storage. The base of the oscillations could be called the 0 bit, and the high point could be called the 1 bit, and these combinations of 0 and 1 bits could be used to store short term memory in biological cultures, performing the functions that the RAM (random access memory) does in computers.
From Sinusoids to Square Waves
Sinusoidal waves have been generated and characterized in several circuits employing various topologies so far, including single cell oscillators, small colony synchronized quorum sensing based oscillators, wide range synchronized oscillators (that synchronize over different channels inside a microfluidic chamber), and even mammalian oscillators. Thus, we used these as a starting point in our project, and thought of ways to transform this sinusoid, so as to generate a square wave.
After brainstorming, basic modeling and simulations, and discussing several possible topologies (which have been linked below), we finally found the 5 node ring topology with our design modifications to be the best solution to the problem.
The basic topology for the system was inspired by the 3 node repressilator, which basically had repressors arranged in a cyclic fashion, where one represses the second, the second represses the third, and the third represses the first. The same schematic, with an even number addition to the number of nodes can generate oscillations, as has been shown by several groups in the literature. Thus, we employed a novel 5 node system in the same fashion, with our design modifications, in order to get the square waves to work.
For this transformation from “y = sin (t)” to “y = square(t)” (as MATLAB calls it), the major strategies for conversion, and how our 5 node oscillator employs these to generate square waves have been outlined below.
Graded to Digital Response – High Cooperativity
The first issue that we identified was that sine waves go up and come down in a slow fashion; slower than what we would want in our square waves. Thus, the first thing to do was to increase the rate of the response, by employing a mechanism wherein the repressors act quickly, and only after a certain point. This would make rise and fall times faster.
As mentioned in the earlier section as well, cooperativity is the property of a repressor that represents the number of molecules of the repressor that need to combine and polymerize before they can repress the promoter of their corresponding gene. Repressor cooperativities bring down the rate of production of the protein in a multiplicative fashion via the hill function, which can be represented as (P represents protein concentration, n represents cooperativity) –
And again, as the cooperativity of the repressor changes, the level of the protein changes as follows –
Therefore, the high cooperativity promotes faster degradation. This fact was exploited by us for this purpose, employing repressors with the highest cooperativity from the 73 TetR homologs that were reported in Voigt’s paper (see references below). Thus, we were able to take care of the rise and fall times using these high cooperativity TetR homologs. The cooperativities and names of the repressors used in our system are as follows –
Node
Repressor Used
Cooperativity Reported
A
Orf2
6.3
B
TetR
3
C
SrpR
3.2
D
PhlF
4.5
E
Bm3RI
4.5
Node | Repressor Used | Cooperativity Reported |
---|---|---|
A | Orf2 | 6.3 |
B | TetR | 3 |
C | SrpR | 3.2 |
D | PhlF | 4.5 |
E | Bm3RI | 4.5 |
Holding the Constant ON State – Time Delay
Once we were ensured that the rates would be faster, we next turned to holding the constant ON and OFF states, a property that square wave oscillators exhibit, and sinusoidal ones don’t.
In order to maintain its state once it is high, we required a time delay between when the level of a node rises and when it falls. This is taken care of by our square wave oscillator by having the two extra repressor nodes. Basically, the two nodes do not qualitatively change the state of the oscillations, since if A has to repress D, then now in our circuit, A represses B, which in turn now cannot repress C, and therefore C is produced, which represses D. Therefore, the effect remains the same, but there is a delay introduced. This is shown diagrammatically below –
Noise Reduction – Low Copy Reporter and Removal of Degradation Tags
Having catered to the issues of rise time and fall time, as well as holding the ON and OFF states for a long time similar to how square waves behave, we then also wanted to have the oscillator be noise free. Noise in biological systems can hamper the oscillations by driving the cells out of phase due to copy number variations of the plasmid, leaky expression, and other factors.
How our square wave oscillator caters to the noise is as follows –
- High cooperativity repressors – High cooperativity ensures a near digital response, and therefore lower concentrations of the repressor would not be able to repress their corresponding node well enough, due to their requirement to dimerize in high numbers.
- Low copy reporter – Demonstrated by Paulsson et al, one of the main sources of noise in the repressilator was that the reporter was placed in a high copy number plasmid. Such high copy numbers (containing pUC19 or pMB1 origins of replication) can have large cell to cell variations in copy number, thereby reporting in a faulty, noisy manner. Therefore, as was done by their group, our square wave oscillator contains the reporter on a low copy backbone (p15A ori).
- Removing degradation Tags – Paulsson et al, in the same paper, reported that the ssrA degradation tags that were used at the end of the repressor genes in the repressilator also employed machinery that was noisy. Removing these tags from his system brought down the variance significantly. Our square wave oscillator employs the same strategy for noise reduction.
References