The Biological Input Processors - The Senders
A perceptron algorithm includes the recognition of input patterns and the conversion of these patterns to an intermediate signal in a weighted manner. Namely, each bit in its activated (“1”) state has a different power to determine the final result.
In our bacterial perceptron the input processors are the senders, which comprise distinct subpopulations. Each one has a specific weight, meaning that it responds to the recognition of the input by producing a different amount of the LuxI synthase, the synthase of OC6, and, consequently, a different amount of OC6 lactone per unit of time. This graded response is achieved via the manipulation of the Ribosome Binding Site (RBS), a major determinant of translation initiation rate and, therefore, protein production.
The Ribosome Binding Site as a plug-and-play device to tune the weights.
Tuning the weights in a biological perceptron corresponds to finding ways to tweak the contribution of each activated sender to the weighted sum. Possible ways to achieve that is transcriptional regulation via promoter strength manipulation, like in Ximing Li’s work [12], variation of the translation rate by changing the RBS, by different mutants of the synthase possessing different enzymatic capability, synthase degradation rate deploying degradation tags, mRNA stability and more.
In our project, translational regulation via manipulation of the RBS sequence was chosen as the determining factor of the weights. RBS sequences are versatile, potentially plug-and-play devices that could be placed upstream of any gene of interest to provide the desired translation rate.
Synthetic RBS Library
Aiming to find the right combination of weights for our pattern recognition problem, a synthetic RBS library using the RBS Calculator [13,14] was generated by our Dry Lab. After guidance of Prof. Howard Salis, the RBS Library Calculator was used to design a degenerate RBS sequence that systematically varies the translation initiation rate. Target translation rates at logarithmically varied numbers were set and a library of synthetic RBS was generated. Then, RBS variants were clustered as “Really Strong”, “Strong”, “Medium” and “Weak” according to the table:
Classes |
Really Strong (>150.000) |
Strong (70.000-90.000) |
Medium (25.000-40.000) |
Weak (<15.000) |
Table 1.
This clustering was arbitrary aiming to create RBS groups that are more probable to differ significantly in their translation initiation rate when tested. From the library, 6 synthetic RBS sequences were randomly selected to be characterized (Table 2).
RBS_Name |
RBS_sequence [Sequence of the RBS variant] |
CDS Start Position (nt) |
ORF_number [Open Reading Frame] |
TIR [Translation Initiation Rate (au)] |
Strength |
synthetic 4 |
TTGGCAACGGGTTATGGGAGGGATGACA |
28 |
1 |
85298.89117 |
Strong |
synthetic 5 |
TTGGCAACGGCTTATGGGAGGTATGTCA |
28 |
1 |
78686.23401 |
Strong |
synthetic 10 |
TTGGCAACGGGTTATCGGAGGGATGACA |
28 |
1 |
31126.53073 |
Medium |
synthetic 11 |
TTGGCAACGGGTTATCGGAGGGATGACA |
28 |
1 |
30580.84267 |
Medium |
synthetic 12 |
TTGGCAACGGGTTATGGGACGGATGTCA |
28 |
1 |
11058.30802 |
Weak |
synthetic 13 |
TTGGCAACGGCTTATCGGAGGTATGCCA |
28 |
1 |
10752.07857 |
Weak |
Table 2.
Results from the RBS calculator regarding the translation initiation rate of the selected synthetic RBS sequences.
Reliable translation elements
The function of a RBS sequence relies heavily on its upstream [15] and downstream [16] genetic context, i.e. the neighboring nucleotide sequence. On the grounds that the above synthetic RBS translation rates were estimates, we wanted to ensure that we had an already tested reliable alternative to express LuxI in our arsenal. Our research led us to BCDs (Bicistronic Design) by Mutalik et.[17] Inspired by overlapping genetic elements, they include short leading peptides followed downstream by the Gene of Interest (GOI). Their design “encodes a 16-amino-acid leader peptide in a first cistron that overlaps by 1 base pair with a variable downstream coding sequence, encoding both a stop and start codon via a −1 frame shift. The leader peptide is synthesized by ribosomes that bind to an upstream SD core sequence (SD1); translation of the downstream GOI is thought to result, primarily, from SD1-directed ribosomes that recognize and reinitiate translation via a second SD site (SD2) that is encoded entirely within the coding sequence of the leader peptide” -Mutalik et. al. The intrinsic helicase activity of ribosomes arriving at the stop codon of an upstream cistron might eliminate inhibitory RNA structures that would otherwise disrupt translation initiation of the downstream GOI.
Given the fact that we are aiming to demonstrate the construction of a biological perceptron on the foundational level, we decided to include the more reliable BCDs to our characterisation, so as some of our calculations to be reproducible to other systems that might include RBS variants but with a different quorum sensing synthase and, therefore, a different genetic context near the RBS. BCDs promise fewer fluctuations when applied to different systems and are therefore closer to be considered plug-and-play devices.
The iGEM distribution kit includes three different BCDs, BCD1 (BBa_J428034), BCD8 (BBa_K2680561) and BCD12 (BBa_K2680529). We decided to deploy these BCDs and we ordered the BCD2 (BBa_K4294102), the strongest variant of all according to Mutalik et. al and BCD14 (BBa_K4294114).
As a preliminary evaluation of their strengths we took into account both the paper’s plots that compare the strength and variability of the BCDs and the RBS calculator.
Figure 1.
BCDs relative performance scores. Derived from Mutalik et. al “Precise and reliable gene expression via standard transcription and translation initiation elements.” Supplementary Material.
RBS_Name |
RBS_sequence [Sequence of the RBS variant] |
CDS Start Position (nt) |
ORF_number [Open Reading Frame] |
TIR [Translation Initiation Rate (au)] |
Strength |
BCD2 |
GGGCCCAAGTTCACTTAAAAAGGAGATCAACAATGAAAGCAATTTTCGTACTGAAACATCTTAATCATGCTAAGGAGGTTTTCTA |
85 |
1 |
287840.97 |
Very Strong |
BCD12 |
GGGCCCAAGTTCACTTAAAAAGGAGATCAACAATGAAAGCAATTTTCGTACTGAAACATCTTAATCATGCTGCGGAGGGTTTCTA |
85 |
1 |
12782.75 |
Weak |
BCD1 |
GGGCCCAAGTTCACTTAAAAAGGAGATCAACAATGAAAGCAATTTTCGTACTGAAACATCTTAATCATGCACAGGAGACTTTCTA |
85 |
1 |
9455.22 |
Weak |
BCD8 |
GGGCCCAAGTTCACTTAAAAAGGAGATCAACAATGAAAGCAATTTTCGTACTGAAACATCTTAATCATGCATCGGACCGTTTCTA |
85 |
1 |
2190.04 |
Weak |
BCD14 |
GGGCCCAAGTTCACTTAAAAAGGAGATCAACAATGAAAGCAATTTTCGTACTGAAACATCTTAATCATGCGGTGGAGGGTTTCTA |
85 |
1 |
2019.63 |
Weak |
Table 3.
Translation initiation rates predicted by the RBS calculator. The TIR of the second SD sequence was used for this analysis, since ,according to Mutalik et. al, it is the main determinant of the GOI's translation rate. The ribosomes from SD1 reinitiate translation from SD2.
It is obvious that the plots and the calculations contradict each other regarding some RBS sequences and therefore wet lab characterisation was necessary.
Senders’ main circuit
Defining the input for pattern generation - The TetR repressor system.
After a long way to decide the right inducer as the input, we chose the well known and characterized TetR repressor system. TetR forms a dimer and binds to tet operators in PTet promoter and blocks the assembly of the transcription machinery. PTet is a medium strength promoter that is constitutively “ON” in the absence of TetR. It contains two tet operators, one inside the core promoter sequence and the other after the -10 hexamer, leading to efficient repression. The system is efficiently induced via anhydrotetracycline (aTc), a tetracycline derivative which binds Tet R with an ~35-fold higher binding constant, being an effective inducer at very low concentrations. Additionally, its antibiotic activity is ~100-fold lower and has minimal toxic effect on E.coli cells in the concentration used for induction[18].
Main characteristics needed from the senders’ induction system were the low leakiness and dynamic range. Steepness would be an extra asset, because it would make our system more compatible with our perceptron ideation if senders responded in an ON/OFF fashion as well. More information about the concept of steep response curves can be found in the Receivers’ circuit design, where it is more thoroughly examined.
An extensive characterization and optimization of inducible systems in E.coli was made by Meyer et. al to create the E.coli Marionette Strains [19]. This comprehensive characterization provides useful information about repressible systems’ characteristics, such as the repressor’s binding and dimerization constants, max and minimum output (i.e. the leakiness). Among the several common systems that are characterized, we thought that the TetR repressor system suits our design goals best. It lies at the lower end of the spectrum regarding its leakiness when it is genome encoded and in the middle when it is plasmid regulated, its dynamic range lies in the middle of the spectrum, and, most importantly, it shows the steepest response curve of all (the highest n constant in the Hill’s function). Even though these parameters were extracted from measurements of E.coli Marionette strains, these strains derive from DH10B (Marionette-Clo) and BL21 (Marionette Pro). Therefore, we believe it was rational to choose the senders’ induction system based on those parameters.
Inducer |
Plasmid |
Regulator |
Promoter |
ymax (RPU) |
ymin (RPUx10-3) |
K (μM) |
n |
Dynamic range |
|
aTc |
pAJM.011 |
TetR |
PTet* |
2.4 |
4.9 |
0.013 |
3.8 |
500 |
Plasmid Encoded Regulator |
aTc |
pAJM.717 |
TetR |
PTet* |
3.2 |
3.6 |
0.012 |
4.4 |
890 |
Genome Encoded Regulator |
Table 4.
Parameters from the E.coli Marionette Strains characterization.
Figure 2.
Senders’ Circuit Overview. TetR is expressed by a constitutive promoter (BBa_J23106 and PlacI promoters were tested, please refer to Engineering Success for more regarding these trials) and the BBa_B0034 RBS. It binds to specific operators in the PTet promoter and blocks RNA polymerase binding and, thus, transcription.
The circuit at its “OFF” state (“0”)
The circuit at its “ON” state (“1”)
Figure 3: Plasmid map of a Level 1 construct with one of the synthetic RBS controlling LuxI translation.
Figure 4: Plasmid map of a Level 2 construct.
Final TetR regulation circuit
As it is explained in more detail in the Engineering Success page, the constructed plasmid-encoded TetR regulator systems did not show a normal induction response. As a result, DH5a-z1 cells, which encode TetR in their genome, were utilized for our perceptron system. These cells were transformed with pTU1 high copy number plasmid containing the Level 1 LuxI TU constructs. This was an important compensation on the grounds that pTU1 plasmid offers resistance to a different antibiotic than pTU2, making experiments with mixed senders and receivers populations more complex.Our initial thoughts about the potential shadowing of the relative RBS strength differences due to the high copy number vector variations were disproved. The selected RBS sequences provided clearly different translation rates, as it can be found in the Results page.
Characterisation of the Ribosome Binding Sites
Characterisation of RBS strength utilizing a LuxI-mNeonGreen bicistronic design with a translation Termination-Reinitiation (TeRe) strategy.
In order to characterize the strength of our senders’ RBS, we had to find a way to measure the relative LuxI synthesis of each subpopulation in the scale of time.
From a very early stage of our designing sessions, we agreed to steer clear of complex and time-consuming measurement methods, such as column-based purification procedures, or processes that required specialized equipment we did not have access to, e.g. Western blotting. Nevertheless, we added a C-terminal hexa-HisTag in our senders’ final constructs, to perform protein analysis through Western blotting, just in case our main scheme did not work as planned.
After excluding these strategies, we decided to adopt the chimeric gene approach by coupling the gene coding for LuxI with a gene producing a visible readout, namely a fluorescent protein. We ruled out placing a linker between the two CDS, as this would require additional research, and incorporated an alternative linking mechanism instead; to combine LuxI’s synthesis with the production of the fluorescent protein, we deployed a phenomenon called Termination-Reinitiation (or TeRe) observed during translation [20]. TeRe occurs when two neighboring genes have overlapping start and stop codons. More specifically, once translation of the upstream gene is terminated, the 30S ribosomal subunit lingers in the area of the stop codon; if it immediately encounters a new start codon, it re-recruits the 50S subunit and reinitiates translation, this time, of the downstream gene, as shown in the figure below. This way, the downstream gene is translated at the same rate as the upstream gene.
Figure 5
Termination-Reinitiation translation mechanism. Figure derived from Huber et. al., 2019 (modified)
Based on this mechanism, we enriched the TAA stop codon of the luxI gene with an extra TG, to convert its terminal A into the initial A of the fluorescent protein’s gene start codon, resulting in the coupling sequence 5’-TAATG-3’ for the characterization of senders’ RBS.
Figure 6. Overlapping stop (luxI) and start codon (mNeonGreen) forming the sequence 5’-TAATG-3’
We opted for mNeonGreen as the fluorescent protein we would utilize for our TeRe-based measurement system, as we had already researched a lot about it -since we had already selected it as a visible readout for our entire perceptron-resembling device (more information about this aspect of our project below)- and believed it would also fulfil the requirements for our senders’ characterization experiments. According to the TeRe mechanism and our predictions, mNeonGreen fluorescence measurements would constitute a relatively accurate indicator of the amount of LuxI synthesized.
Figure 7.
Overview of the translational coupling circuit for the characterization of the LuxI expression.
Circuit at its off state (“0”)
Circuit at its on state (“1”)
Figure 8: Plasmid map of the Level 1 construct with the TeRe design.
Figure 9: Plasmid map of the Level 2 construct with the TeRe design.
Unfortunately, the circuit did not work as expected and it had to be replaced with a different strategy.
Characterisation of the RBS strength using the 36nt luxI - sfGFP fusion.
As it was stated, the function of the Ribosome Binding Site is often unpredictable, since it depends on its interaction with the upstream and downstream genetic context. Therefore, a characterisation of our RBS strength regarding the translation of LuxI would probably not be accurate, if a fluorescent reporter was simply placed after them, since the genetic context (i.e. the nucleotide sequence and so its interaction with the RBS) would differ.
To prove the reliability and context independence of the BCDs, Mutalik et. al engineered several fusions of the first 36 nucleotides of commonly used regulators and enzymes with superfolder Green Fluorescent Protein (sfGFP). These first 36 nucleotides are sufficient enough to introduce a similar genetic context of the whole sequence they originated to the RBS [21]. sfGFP is an engineered version of the green fluorescent protein that can be fused with poorly folded peptides without misfolding and losing its functionality [22]., making these 36nt fusions a simple and probably reliable way to characterize translation rates provided by the same RBS for different coding sequences. A similar strategy was recently deployed to determine the effect of the spacer sequence between a RBS and the start codon in the context dependency of the RBS [23].
Based on the above information, we built a fusion of the first 36 nucleotides of the LuxI synthase with sfGFP to characterize the relative strength of our deployed RBS in the context of the LuxI coding sequence. sfGFP maintained its fluorescent properties and allowed us to conduct the characterisation. The characterisation results can be found in the Results section.
Figure 10.
Final RBS characterisation circuit.
FInal RBS characterization circuit overview
Circuit at its off state (“0”)
Circuit at its on state (“1”)
Figure 11 Plasmid map of a Level 1 plasmid with the 36ntLuxI-sfGFP fusion CDS.
Figure 12 Plasmid map of a Level 2 plasmid with the 36ntLuxI-sfGFP fusion TU. TetR is constitutively expressed by the promoter BBa_J23106.
Figure 13 Plasmid map of a Level 2 plasmid with the 36ntLuxI-sfGFP fusion TU. TetR is constitutively expressed by the promoter PlacI. This is an alteration of the original circuit that was designed during troubleshooting to explore whether a change in TetR dynamics would improve the system's response to aTc.
The 36nt LuxI-sfGFP fusion was proved to be a successful part in the characterization to determine the relative strength of the RBS variants in the context of the LuxI CDS. Please, refer to the Engineering Success tab and Results tab for more.
Cross-talk
After deciding that the input patterns will be generated by a single chemical inducer, a problem that emerged was how to deal with cross-talk among senders. For this aspect of our project, cross-talk is defined as the unintended stimulation of supposedly inactivated sender cells by residual inducer molecules when they are mixed with activated senders to produce a specific pattern based on the strengths of their RBS.
1. We came up with various plans to tackle this situation:
2. One idea revolved around initiating the production of a tTa transactivator upon mixing senders’ subpopulations. The expression of a molecule that would with TetR for the binding of aTc could contribute to sequestering the inducer, thus keeping uninduced senders inactivated [24].
3. Another idea pertained to designing a CRISPR-Cas9 mediated mechanism that would knock out the luxI gene in inactivated senders upon mixing them with their induced counterparts to prevent cross-talk phenomena.
4. Similarly to the previous thought, a recombinase could be induced to be produced when mixing sender cells, to misplace the luxI gene, thus deactivating it.
5. Certainly, a simple wash of the different sender cultures’ supernatants could also resolve the issue by removing the inducer altogether.
6. We even thought about harnessing the photosensitive nature of our inducer, aTc, and exposing it to blue light to photobleach it before mixing different senders’ subpopulations, so that it could not activate sender cells anymore [25].
All potential solutions mentioned above present a number of challenges, such as leading us to a vicious circle by perpetuating cross-talk issues or requiring specialized optogenetics equipment which we could not afford, so we eventually culminated in a completely different approach. During our investigation of potential solutions and after consulting with our advisors, we encountered a paper analyzing the tetracycline destructase family, a group of enzymes that degrade tetracycline and its derivatives, such as aTc [26]. We decided to incorporate such an enzyme, called TetX, in our cross-talk tackling strategy by co-transforming it in senders’ subpopulations we wish to keep uninduced. As already mentioned, we designed this construct with pTKEI-Dest as its backbone, because its ori (incompatibility group A) would not interfere with the replication of the senders’ final construct (incompatibility group B). TetX is placed downstream the Ptac Promoter, which contains a lac operator. pTKEI vector incorporates a constitutive LacI expression and, therefore, TetX production is induced by IPTG.
Figure 14.
pSend-TetX cross-talk construct (created with SnapGene)
Having concluded the design of our cross-talk system, the parameters we had to take into account during every induction experiment with the combination of induced and uninduced can be summarized by the following three points:
1. All sender’ subpopulations that would have to remain inactivated to create the desired pattern should be co-transformed with pSend-TetX as described above.
2. All sender’ subpopulations that would have to be induced, but still be co-cultured with inactivated senders, should be co-transformed with a second plasmid as well, to acquire resistance to Kan too (the antibiotic of pSend-TetX); to achieve that, we selected pJUMP27-1A from the iGEM DNA Distribution kit (plate 1, position 4E), which is Kan resistant, to be our “dummy” senders’ plasmid based on its origin of replication (senders’ main construct belongs to incompatibility group B, while pJUMP27-1A belongs to incompatibility group C).
3. Similarly to induced senders, receivers being co-cultured with senders in order to extract the desired pattern should also have resistance to Kan; since receivers’ main constructs are classified as either incompatibility group B plasmids (p15A ori) or incompatibility group A plasmids (colE1 ori), pJUMP27-1A (pSC101 ori) could be utilized as the “dummy” plasmid for receivers too.
The Biological Activation Function - The Receivers
In a common perceptron, the input patterns are converted into a weighted sum which is then processed by an activation function. This is the point of the classification of the initial pattern into a specific class. If the weighted sum surpasses a set threshold, the “neuron” gets activated and gives an output “1”. Otherwise, the output is “0”.
In order to program the receiver cells to function similarly to the activation function, a circuit that imitates it is needed. Namely, a circuit that performs like an ON/OFF switch.
Often inspired by computer science data processing modules, synthetic genetic and metabolic circuits have been engineered to function like logic-gates, oscillators and more[27]. One of the main principles of Synthetic Biology is the use of characterized genetic elements, i.e. parts, that can be combined to perform complex and desired tasks. In the relatively simple context of a system that responds to the concentration of an input by an output expression, the relation of the output signal to inducer concentration forms a sigmoid curve if plotted on a Xy axis. These plots are the “response curves” and the mathematical equation that describes this response to the input is called a “transfer function”. Such systems often need the suitable tuning to function as expected and be suitable for combination with other systems and synthetic biologists have established various possible ways to tune a systems’ response curve in terms of dynamic range, leakiness, steepness and more[28]. Before the receivers’ genetic circuits are described, it is crucial to define in a nutshell some important terms:
1. Dynamic Range: The difference in output between the “OFF” and the saturated “ON” state.
2. Leakiness: The output production in the “OFF” state.
3. Steepness: A steep response curve is described by a large change in output within a narrow input concentration range; small changes in input result in large changes in output. Such switch-like biochemical circuits are often referred to as ultrasensitive or bistable.
The perceptron’s activation function is a switch-like function that gives a positive (“1”) output if the weighted sum is surpassed. It becomes, therefore, clear that if the activation function is to be imitated in a biological system, an ultrasensitive circuit, i.e. with a steep activation function, that responds to the quorum sensing signal summation is needed.
Figure 1.
Response curves with different steepness. The typical response curve of a simple inducer-output system is a sigmoid curve and its steepness depends on the intrinsic properties of its components (e.g. the repressor’s binding cooperativity). With specific circuit architectures, steeper, switch-like response curves can be engineered. Such circuits are described as “ultrasensitive” or “bistable”.
A step-by-step approach to achieve ultrasensitivity
A number of possible ways at different levels of gene expression and protein function to build a bistable switch exist[28]. A well known example is the construction of the first toggle switch in E.coli; a circuit that responds in an ON/OFF fashion effectively[29]. Such bistable circuits are based on mutual inhibition relationships between transcriptional repressors in the circuit. At the level of transcription, binding cooperativity of transcription factors and the number and topology of operators has been proven to correlate with a response’s steepness. Cooperative binding of multiple transcription factors to the same promoter or a biochemical process mimicking this effect may promote ultrasensitivity as well [28]. Moreover, circuit architectures with proven enhanced steepness are transcriptional positive feedbacks [30] and transcriptional cascades [31]. Further strategies include a sequestering molecule that binds a circuit component and prevents its effect. Sequestration has been achieved using sRNAs [32] that bind to mRNA, proteins that bind to transcription factors[33] or translation machinery components [34] and decoy DNA operators [35,36] that titrate the transcription factor away from the output promoter. The above approaches are often coupled in such a way that their respective effects generate a greater impact on the final response [34].
Irrespective of the existing literature, every circuit that is built requires testing and tuning. On the grounds of our narrow time window to perform extensive characterisation experiments, we followed a step-by-step approach regarding the receivers’ genetic circuits. Starting from a simple open loop circuit, which was necessary to obtain vital parameters of our parts (refer to Model section), two variations of a positive feedback and two variations of a coherent feed forward loop coupled with a positive feedback at its first step were built. Progressively introducing new components to our circuits enables troubleshooting and better interpretation of the results, a process that was also facilitated by estimations by our Dry Lab (see Engineering Success for more).
Open Loop Circuit (OL or OpLo)
This is a circuit of the LuxR quorum sensing regulator at its simplest form. LuxR is constitutively expressed. When concentration of OC6 increases, it binds with LuxR, LuxR forms stable dimmers and binds to its binding site (lux box), upstream of the -35 region of the Plux promoter, activating the transcription of the downstream output.
Figure 2.
Flowchart of the Open Loop circuit
OL genetic circuit overview. LuxR is constitutively expressed under the control of the constitutive BBa_J23105 promoter and BBa_B0034 RBS.
Output (mNeonGreen) production is controlled by the Plux inducible promoter and BBa_B0034 RBS.
OL at is “OFF” state (‘“0”)
OL at its “ON” state (“1”)
We constructed this circuit in both a medium (ColE1 origin of replication) and a low (p15a origin of replication) copy number pTU2 plasmid vector so as to choose the copy number that provides the most suitable dynamic range, leakiness and steepness. The medium copy number plasmid was chosen and used for the rest of the constructs (see Results section for more details) and the parameter estimation for our Model.
Figure 3.
a. pTU2 (ColE1) vector with OL circuit.
b. pTU2 (p15a) vector with the OL circuit.
Positive Feedback Circuit (PF and PFc).
One already tested and widely utilized method to engineer a bistable circuit are transcriptional Positive Feedbacks. In such circuits, an activator of an output induces its own production as well after induction with its respective ligand. A similar circuit was built by Li et. al [12] to achieve a steeper activation function of the CinR quorum sensing activator system.
In our case, Plux was placed upstream of the LuxR coding sequence and the output coding sequence. Therefore, LuxR activates its own production after binding with OC6 and establishes a positive feedback.
Two possible topologies exist; LuxR being controlled only by the Plux promoter or providing an extra TU for constitutive LuxR production.
In the former circuit topology, Plux leakiness provides a low basal LuxR concentration in the cell. This concentration is sufficient for the LuxR dimer formation after induction with OC6. The activated LuxR dimers promote the expression more LuxR resulting in a positive feedback.
The latter provides a higher basal LuxR concentration inside the cell. These circuits were both constructed and tested in order to find the most suitable regulator and circuit dynamics for our system.
Unfortunately neither topology resulted in the desired results. For more information and engineering considerations, please refer to Results and Engineering Success tabs respectively.
Figure 4.
a.Flow chart of a Positive Feedback circuit
b.PF (Positive Feedback) genetic circuit overview. LuxR is placed downstream the Plux inducible promoter and BBa_B0034 RBS. The leakiness of Plux promoter provides a basal LuxR concentration in the cell that will be activated by OC6 and promote the production of more LuxR and of the output, resulting in a positive feedback loop. Output (mNeonGreen) is placed downstream of a Plux promoter and BBa_B0034 RBS.
c.PF (Positive Feedback) genetic circuit overview. LuxR is placed downstream of the promoter Plux inducible promoter and BBa_B0034 RBS. The d.leakiness of Plux promoter provides a basal LuxR concentration in the cell that will be activated by OC6 and promote the production of more LuxR and of the output, resulting in a positive feedback loop. Output (mNeonGreen) is placed downstream Plux promoter and BBa_B0034 RBS.
e.PFc (Positive Feedback with constitutive LuxR production) genetic circuit overview. Two LuxR TUs exist in this circuit. One expresses LuxR constitutively and provides a higher LuxR concentration in the cell. The other LuxR TU (with the Plux promoter) is responsible for the generation of the positive feedback after induction with OC6. Output (mNeonGreen) is placed downstream of the Plux promoter and BBa_B0034 RBS.
f.PFc at its “OFF” state (“0”)
g.PFc at its “ON” state (“1”)
Figure 5.
a. PF circuit in the pTU2 (ColE1) plasmid vector.
b. PFc circuit in the pTU2 (ColE1) plasmid vector.
A coherent feed forward loop coupled with a Positive Feedback (PFR, PFcR)
For transcriptionally regulated gene expression, one approach to achieve a switch-like response is an engineered promoter onto which multiple transcription regulators bind cooperatively or a biochemical mimic of this effect.
This led us to the design of the part BBa_K4294301, a hybrid promoter which includes a LuxR activator binding domain (luxbox) upstream the -35 hexamer and a PhlF operator (phlO) that spans from the core sequence and the -10 hexamer to the proximal promoter sequence.
Figure 6.
Plux-phlO hybrid promoter.
PhlF is under the transcriptional control of a LuxR repressible promoter (Plux_rep, BBa_J107103). However, in this promoter the LuxR binding domain is not placed upstream of the -35 hexamer, but inside the core sequence between the -35 and -10 regions. Due to this localization change, the effect of the active LuxR dimers is reversed; instead of activating the production of PhlF, they repress it [37]
Even though LuxR and PhlF do not bind cooperatively to their respective binding sites, their dynamics are indirectly connected. PhlF was chosen since it is an effective repressor with a low association constant K, indicating its ability to form functional dimers in lower concentrations (according to the Ecoli Marionette strains characterization [19]). Therefore, in “intermediate” OC6 concentrations, it could possibly still reinforce repression to the output irrespective of its lower production rate (due to the effect of LuxR on the Plux_rep promoter) and the activation effect of LuxR to the output’s promoter. LuxR establishes a positive feedback loop to propagate its effect and both positive feedback topologies (with and without constitutive LuxR expression) were designed and tested. A very efficient ssrA [38] degradation tag (NDENYALAA) was added to PhlF to ensure its rapid degradation after maximum repression of the phlF gene by LuxR and to avoid PhlF concentrations that might lead to a complete circuit block in every inducer concentration. The degradation tag was empirically selected and further analysis of different degradation tags will be needed for a potential optimization of the circuit’s response.
Note: These circuit architectures are referred to as Positive Feedback + Repressor (PF+R) and Positive Feedback with constitutive LuxR expression + Repressor (PFc+R) in other pages as well.
Figure 7.
a. Flow chart of PFR (Positive Feedback + Repressor) and PFcR (Positive Feedback with constitutive LuxR production + Repressor) circuits. This topology resembles the one of a coherent feedforward loop with an added positive feedback involving the first step.
b. PFR genetic circuit overview. The output is controlled by the hybrid promoter BBa_K4294301, which is activated by LuxR and repressed by PhlF. LuxR is placed downstream of the Plux promoter, establishing a positive feedback after induction with OC6 as previously described. LuxR also binds to a Plux repressible promoter (BBa_J107103), which controls the production of PhlF, blocking PhlF expression after induction with OC6. Every coding sequence is placed downstream a BBa_B0034 RBS. PhlF has a high efficiency NDENYALAA degradation tag.
c. PFR circuit at its “OFF” state (“0”)
d. PFR circuit at its “ON” state (“1”)
e. PFcR genetic circuit overview. This circuit is similar to the PFR circuit, but with an added constitutive expression of LuxR, which provides a higher basal LuxR concentration in the cell.
f. PFcR at is “OFF” state (‘“0”)
g. PFcR at its “ON” state (“1”)
Figure 8.
PFR circuit in the pTU2 (ColE1) plasmid vector.
PFcR circuit in the pTU2 (ColE1) plasmid vector.
After induction with OC6, PFR circuit response showed a significant difference in steepness in comparison with the original open loop circuit. For more information, please refer to the Results page and the Model page for the circuit’s model parameters.