Part:BBa_K4387004

Nitric Oxide Sensing Genetic Circuit with Promoter BBa_K2116002

Usage and Biology

This composite part consists of the ETH promoter pNorV, a superfolder GFP preceded by one strong ribosomal binding site (BBa_B0034), the NorR regulator, and a double forward terminator. We chose a high-copy backbone from Twist Bioscience for this part. We wanted to compare this ETH NorV promoter to the pNorVβ promoter and see which one was better suited for sensing nitric oxide at lower concentration ranges.

This construct was tested in the bacterial strain E.coli Nissle 1917.

Characterization

Measuring parts with different approaches and comparing them to provide a more insightful and multilayered characterization is essential in Synthetic Biology. Here, we focused on two methods:

(i) time-lapse plate reader assays to measure the sensitivity of our circuit to NorR in a dynamic manner and under different concentrations of inducer; and

(ii) endpoint flow cytometry assays to measure the behavior of our circuits at the single-cell scale.

With the first assay, we uncovered essential kinetic information about the circuits on the populational level (every measurement is an average of the individual expression patterns in the sample). With the second assay, we delved deeper into the cell populations to characterize other essential properties of our system, such as expression noise and dose-dependent responses to different inducer concentrations.

We performed all analyses using in-house R scripts.

Time-Lapse Plate Reader Assay

Figure 1a: Response of the Genetic Circuit pNorV (2016 ETH Team) including 1 RBS to DETA/NO induction. We observe that the GFP expression is ~2300 GFP/OD for NO-concentration of 2mM. Comparing this construct with the others in figure 1b and 1c, we could prove that this genetic circuit with the 2016 ETH promoter is the least effective out of the other constructs with pNorVβ.

To make our experiments reproducible, during plate reader assays (PHERAstar FSX - λEx: 485 nm, λEm: 530 nm), we measured each sample for 16 hours at 37°C and constant orbital shaking, using three biological replicates (three individual colonies per circuit) and three technical replicates (three wells per biological replicate).

We performed the data analysis as follows:

1. Subtracted blanks from raw data
2. Normalize GFP by OD600
3. Calculate technical means of GFP/OD600 normalized data
4. Calculate the biological means of GFP/OD600 normalized data
5. Calculate the biological standard deviation of GFP/OD600 normalized data
6. When necessary, perform imputation. Usually, the first normalized measurements are noisy and unreliable as OD600 values can be very low and significantly impact normalization. Thus, when individual normalized values are extremely high or low (sometimes negative due to blank correction), imputation was used following a na_kalman() function from the ImputeTS R package.
7. Additional transformations, such as log transformations.

Hence, our plots (Figure 1a and Figure 1b) show the averages and standard deviations for the biological replicates for each sample for each time point.

Figure 1b: Response of Different Gene Circuits to DETA/NO Induction. Induction of plasmids BBa_K4387003 to BBa_K4387007. For each plasmid, it plots the emission of GFP over time, after induction with 6 different DETA/NO concentrations. Construct pNorVβ including 2 RBS has the highest GFP emission, whereas pNorV of the 2016 ETH team has the lowest GFP emission.

Figure 1c: Dose Response of all Constructs to DETA/NO. Plotted are the highest GFP expression levels under different concentrations of DETA/NO

Because the standard deviations overlap, we thought that we might be able to reduce the standard deviations and get a clearer result if we had more samples. To get more samples, we also performed a flow cytometry.

Endpoint Flow Cytometry Assay

For the flow cytometry experiment, cell cultures were grown overnight in LB medium supplemented with antibiotic, diluted in 2mL of M9 (supplemented with glucose, cas amino acids and an antibiotic) in a 1:10 ratio (v/v), induced with different NO concentrations and grown for 7 hours in a shaker (37C, 220 RPM). Samples were then chilled in ice to halt cell growth and diluted in 1mL of cold PBS (1:500 v/v ratio). A total of 100,000 cells per sample was measured in a BD FACSCanto II flow cytometer (FSC: 625V, SSC: 420V, FITC: 650V, Event threshold: FSC & SSC > 200, Channel: FITC (λEx 488 nm / λEm. 530/30 nm, High flow rate: ~ 10,000 events/s).

Figure 2a: Violin Plots and Boxplots for the Genetic Circuit Containing 2022 UZurich's pNorVβ Promoter Including 1 RBS and the Genetic Circuit Containing 2016 ETH's pNorV Promoter Including 1 RBS. These plots represent the distribution of FITC-H (green fluorescence) values for the two different promoters at 4 different concentrations of DETA/NO. Each measurement consisted of 100'000 cells.

Figure 2b: Noise Levels of 2022 UZurich's pNorVβ Promoter and 2016 ETH's pNorV Promoter.: Noise levels are measured as the coefficient of variation for construct sample and NO concentration. The constructs show low noise levels at high concentrations.

The results showed that our construct with pNorVβ has a higher response to induction with DETA/NO than the construct with pNorV (figure 2a). It further revealed that the constructs with pNorVβ and pNorV have relatively low noise compared to the negative control (figure 2b). Due to these reasons and the large sample size (100'000 sample points per plasmid), we can infer that pNorVβ is more responsive to DETA/NO induction than pNorV.

Sequence and Features

Assembly Compatibility:

• 10
  COMPATIBLE WITH RFC[10]
• 12
  COMPATIBLE WITH RFC[12]
• 21
  INCOMPATIBLE WITH RFC[21]

  Illegal XhoI site found at 742
  
• 23
  COMPATIBLE WITH RFC[23]
• 25
  COMPATIBLE WITH RFC[25]
• 1000
  COMPATIBLE WITH RFC[1000]