Analysis of the leakage of gene repression by an artificial TetR-regulated promoter in cyanobacteria
© Huang et al. 2015
Received: 29 June 2014
Accepted: 8 September 2015
Published: 19 September 2015
There is a need for strong and tightly regulated promoters to construct more reliable and predictable genetic modules for synthetic biology and metabolic engineering. For this reason we have previously constructed a TetR regulated L promoter library for the cyanobacterium Synechocystis PCC 6803. In addition to the L03 promoter showing wide dynamic range of transcriptional regulation, we observed the L09 promoter as unique in high leaky gene expression under repressed conditions. In the present study, we attempted to identify the cause of L09 promoter leakage. TetR binding to the promoter was studied by theoretical simulations of DNA breathing dynamics and by surface plasmon resonance (SPR) biosensor technology to analyze the kinetics of the DNA–protein interactions.
DNA breathing dynamics of a promoter was computed with the extended nonlinear Peyrard–Bishop–Dauxois mesoscopic model to yield a DNA opening probability profile at a single nucleotide resolution. The L09 promoter was compared to the L10, L11, and L12 promoters that were point-mutated and different in repressed promoter strength. The difference between DNA opening probability profiles is trivial on the TetR binding site. Furthermore, the kinetic rate constants of TetR binding, as measured by SPR biosensor technology, to the respective promoters are practically identical. This suggests that a trivial difference in probability as low as 1 × 10−4 cannot lead to detectable variations in the DNA–protein interactions. Higher probability at the downstream region of transcription start site of the L09 promoter compared to the L10, L11, and L12 promoters was observed. Having practically the same kinetics of binding to TetR, the leakage problem of the L09 promoter might be due to enhanced RNA Polymerase (RNAP)-promoter interactions in the downstream region.
Both theoretical and experimental analyses of the L09 promoter’s leakage problem exclude a mechanism of reduced TetR binding but instead suggest enhanced RNAP binding. These results assist in creating more tightly regulated promoters for realizing synthetic biology and metabolic engineering in biotechnological applications.
Transcription regulation plays a major role in controlling gene expression and consequently in altering biochemical reactions catalyzed by these gene products. Hence, it has high impact in synthetic biology and metabolic engineering [1–3]. The regulation relies on DNA–protein interactions between a transcription factor and a promoter [4, 5]. Strong and tightly regulated promoters are the prerequisite to realize modular control in synthetic biology , and genetic control in metabolic engineering .
Promoter data from previous in vivo measurements 
15.6 ± 0.1
2.88 ± 0.01
5.15 ± 0.01
17.6 ± 0.1
0.235 ± 0.003
71 ± 1
19.1 ± 0.1
0.236 ± 0.003
77 ± 1
0.043 ± 0.003
0.022 ± 0.003
1.9 ± 0.3
Inspired by the observed long-range effect of a flanking single-nucleotide polymorphism on changing the binding affinity of the eukaryotic YY1 transcription factor , we investigated whether a flanking point mutation would have a similar effect on the binding characteristics of the L09 promoter in comparison with the L10, L11, and L12 promoters. We re-simulated the DNA breathing dynamics of the four promoters by the Extended nonlinear Peyrard–Bishop–Dauxois (EPBD) mesoscopic model to reach a single nucleotide resolution that our previous study could not achieve . We also determined the kinetics of the interaction between TetR and each promoter using a surface plasmon resonance (SPR) biosensor assay.
DNA breathing refers to the transient opening of base pairs of double-stranded DNA molecules subjected to thermal motions fluctuating at physiological temperature. DNA breathing dynamics computes these spontaneous opening and re-closing events of a double stranded DNA and presents the results as a DNA opening probability profile . The profile shows characteristic patterns in different probability, length, amplitude, and lifetime of DNA local transient separation (i.e. DNA bubble) . A strong correlation between DNA breathing dynamics and a transcription factor binding was evidenced to predict the location of the TSS and the potential regulatory sites [12, 13]. Maintaining the integrity of the base-specific interactions in a regulatory site, base substitutions or modifications in the flanking bases might lead to suppressed or enhanced binding of a transcription factor [9, 13, 14].
Kinetic studies using SPR biosensor technology provide detailed kinetic and mechanistic insights into biomolecular interactions. DNA–protein interactions using SPR biosensor assays have been successfully applied , including studies of TetR–tetO interactions . The assay principle is to coat a sensor surface with streptavidin and then capture biotinylated DNA via the strong biotin–streptavidin interaction [17, 18]. Once the target DNA is captured at low surface densities (20–50 RU), real-time interaction studies with the ligands of choice can be performed.
Identifying what causes leaky gene expression under repressed conditions of transcriptional regulation is important for constructing tightly regulated promoters, and consequently, for developing and realizing synthetic biology and metabolic engineering for biotechnological applications. Reasons based on simulations and experimental data to explain the leakage problem of the L09 promoter are discussed and future directions are proposed.
The reverse primer is 5′-CGCGCAGGATGGGCACCACCCCGGTGAACAGCTCC-3′ and biotinylated on its 5′-end. The DNA fragments were purified by the Zymoclean Gel DNA Recovery Kit (Zymo Research, Irvine, USA) from a 1.2 % agarose gel after gel electrophoresis.
DNA breathing dynamic simulation
Monte Carlo simulations on the Extended nonlinear Peyrard–Bishop–Dauxois (EPBD) mesoscopic model and its parameters describing DNA breathing dynamics  were performed with MATLAB (MathWorks, Natick, USA) using the same set  of 2000 different seeds of random number generator and with parallel computing with its distributed computing toolbox. A 160 bp DNA sequence (Fig. 1) containing clamp sequences on each end of a strand was simulated at 303 K with periodic boundary conditions to prevent the end effect. Each realization takes 2.1 × 107 steps and as the 1 × 106 steps reaches the initial equilibrium and then records every 500 steps to have 40,000 snap shots of the displacements of a base pair. Every accepted configuration in an advanced step is determined by the standard Metropolis algorithm. From 40,000 recorded displacements of each base pair in the 160 bp DNA sequence, if the displacements at a base pair and its following consecutive 3–10 base pairs are larger than 2.8 Å, it counts one opening event at the first base in the defined DNA bubble length from 4 to 11 bp. The opening probability of a base pair is the ratio of summed opening events to 40,000 recorded displacements. The DNA opening probability profile is averaged from 2000 realizations.
SPR biosensor based interaction analysis
The SPR based interaction studies were performed on a Biacore 2000 instrument (GE Healthcare, Uppsala, Sweden) using CM5 sensor chips (GE Healthcare). Streptavidin (Sigma Aldrich, St. Louis, USA) was immobilized using amine coupling chemistry at a temperature of 25 °C and a flow rate of 5 μL min−1. The surface was activated by injecting 200 mM EDC and 50 mM NHS for 7 min. Streptavidin (100 μg mL−1 in 10 mM Na-acetate, pH 5.5) was injected for 10 min and the surface was deactivated by injecting 1 M ethanolamine, pH 8.5, for 7 min. Finally, the surface was washed with five injections of 50 mM NaOH (30 s, 30 μL min−1), essentially as described previously [15, 17]. The running buffer consisted of 10 mM Hepes, 150 mM NaCl, 0.05 % Tween-20, pH 7.4 (HBS-P).
Promoters (L09, L10, L11, and L12) were diluted into HBS-P to a final concentration of 10 nM and injected over streptavidin containing sensor surfaces using a flow rate of 5 μL min−1 until surface densities of approximately 40 RU were reached. The surface was then washed at a flow rate of 30 μL min−1 by injecting 1 M NaCl for 100 s followed by 0.05 % SDS for 30 s. A separate flow cell with immobilized streptavidin was used as a reference surface.
TetR (Imgen BioSciences, Fall River, USA) storage buffer was exchanged to HBS-P using Protein desalting columns (Thermo Scientific) according to the manufacturer’s instructions. For kinetic analysis, the temperature was increased to 30 °C. TetR was injected in twofold concentration series (62–4 nM) at a flow rate of 60 μL min−1 for 90 s. The dissociation time was 240 s. HBS-P was used as running buffer. After every TetR injection, the surface was regenerated by injecting a mixture of 1 M NaCl and 1 M MgCl2 for 30 s (60 μL min−1). For interaction studies of induced TetR, the protein was incubated with 2 uM aTc and 2 mM MgCl2. The dilution series was performed into HBS-P supplemented with 2 uM aTc and 2 mM MgCl2.
Trivial difference in DNA opening probability profiles between the L09 promoter and other L promoters
Identical kinetic characteristics of TetR-promoter interactions
Kinetic parameters for the interaction between TetR and promoters L09, L10, L11, and L12
k 1 (s−1 M−1)a (pk1)
k −1 (s−1)a (pk−1)
K D (M)b (pK D )
2 × 105 (−5.2 ± 0.1)
1 × 10−3 (3.0 ± 0.1)
6 × 10−9 (8.2 ± 0.1)
2 × 105 (−5.3 ± 0.1)
1 × 10−3 (3.0 ± 0.1)
5 × 10−9 (8.3 ± 0.1)
2 × 105 (−5.3 ± 0.3)
1 × 10−3 (3.0 ± 0.1)
6 × 10−9 (8.2 ± 0.2)
2 × 105 (−5.3 ± 0.2)
2 × 10−3 (2.8 ± 0.1)
8 × 10−9 (8.1 ± 0.1)
For developing promoters regulated in a wide dynamic range by TetR, the L promoter library was constructed previously for Synechocystis PCC 6803 (Synechocysits) . A point mutation with adenine, thymine, cytosine, and guanine at 2 bp downstream of the −10 element generated the unique L09 promoter. To understand why a point mutation makes L10 and L11 promoters tightly regulated, while the L09 promoter is leaky upon TetR binding (Table 1), potential differences in DNA–protein binding were investigated. The developed SPR biosenor assay allowed to reliably measure the interactions between TetR and the promoters L09, L10, L11, and L12. The approach of capturing biotinylated DNA via covalently immobilized streptavidin has been described in detail previously [17, 18]. For this type of assay, the DNA is biotinylated at a defined position (here 5′ end of template strand) leading to a uniform orientation when immobilized on the sensor surface. The high affinity of the streptavidin–biotin interaction leads to a very stable sensor surface with practically no baseline drift. Furthermore, the structural stability of both streptavidin and the captured DNA allow interaction studies under a long period of time without a considerable loss of surface functionality. The surface density of the immobilized DNA should be reduced to a minimum to increase the accessibility of the binding sites. Different DNA capture levels have been reported, ranging from 300 to 600 RU  down to 0.7–3.5 RU . Generally, DNA levels from 20 to 50 RU are recommended [17, 18] and these were the levels that were achieved in the present study.
The kinetic analysis of the TetR-promoter interactions shows that the association and dissociation constants and the resulting equilibrium dissociation constants are practically identical (Table 2). Kinetic studies of TetR interactions with a 40 bp tetO-containing DNA fragment have been performed previously . Although a direct comparison is not valid due to the different DNA sequences, the conditions in the SPR based interaction analysis were similar and the determined affinity is in the low nanomolar range as well (≈0.2 nM).
The kinetic model that was fitted to the experimental data, a reversible 1-step model (Scheme 1), might seem inappropriate to describe the TetR-promoter interaction. There are two TetR binding sites on the promoter that would make the use of a model describing the parallel and independent interactions to those more appropriate (2-site model). However, two independent and parallel interactions that have identical kinetics lead to an overall binding curve that resembles a 1-step interaction. This is supported by comparing the theoretical R max values with those determined by the global nonlinear regression analysis. The capture levels of the promoters were between 30 and 50 RU corresponding to theoretical R max values in the same range assuming a valency of 2. The experimental R max levels were in the range of 30–50 RU, thereby supporting the parallel interaction of TetR at both promoter binding sites and indicating that almost all DNA binding sites were accessible.
The dimeric structure of TetR would further suggest to use a model describing the multivalent interaction of each dimer. Previously, the mechanism of the interaction between TetR and the tetO-containing promoter was analyzed using stopped-flow measurements and appears to involve the formation of an initial complex with association rate constant of 3 × 108 s−1 M−1, and the subsequent formation of a second more stable complex . Whether a similar 2-step mechanism occurs during the interaction of TetR with the promoters L09, L10, L11, and L12 is not known because the DNA sequences employed are similar but not identical. Furthermore, association rate constants in the range of 1 × 108 s−1 M−1 cannot be measured with the instrument employed. The fact that the reversible 1-step model (Scheme 1) gave an adequate description of the experimental data might reflect that it describes the rate limiting step of the interaction, namely the specific formation of the stable TetR-promoter complex.
Because DNA breathing dynamics simulation has shown a strong correlation with the DNA–protein interactions [12, 13], the trivial difference between DNA opening probability profiles at two TetR cognate sites (Fig. 3) might have caused different TetR binding to the L09, L10, L11, and L12 promoters. However, kinetic rate constants obtained in the SPR-based analysis are identical (Table 2). The observed trivial difference in the computed DNA opening probability does not lead to a difference in the kinetic rate constants. On the contrary, for the YY1 transcription factor, an obviously reduced opening probability at its cognate binding site is consistent to the loss of binding in the ChIP experiments due to the flanking point mutation . These results together indicate that flanking point mutations in the TetR-regulated L promoters do not have long-range effects on the binding of TetR.
If TetR displays practically the same kinetic rate constants when interacting with the L10, L11, and L12 promoters, what could be the reason causing leaky gene expression by the L09 promoter? The repression in transcription regulation is due to the steric hindrance upon a repressor’s binding to its cognate site in the vicinity of the core promoter to prevent RNAP binding . Therefore, under the same repression, the leakage might result from the enhanced RNAP binding to the L09 promoter. Comparing the DNA opening probability profiles in Fig. 2, L09 has a larger DNA bubble at the +2 to +4 base pairs and at the +20 base pair. These regions have a critical role in downstream interactions between RNAP and promoter . Specifically, the contacts with double-stranded DNA at +2 to +4 determine the formation and stability of RNAP-downstream fork junctions complex and the length of contacts possibly exceeding to the +20 base pair also assists the formation of a promoter open complex. The binding difference in this region might cause the L09 promoter to open more easily and form a more stable RNAP-promoter open complex. This might then lead to enhanced binding of RNAP to the L09 promoter. Further SPR based studies could be performed with SigA, the major sigma factor in Synechocystis under normal growth conditions , to measure its kinetics to the L09, L10, L11, and L12 promoters and also perform competition experiments against TetR binding to the promoter. However, as discussed in  the RNAP is structurally different between E. coli and cyanobacteria as well as the sigma factors used in cyanobacteria compared to in E. coli making such experiments much more challenging. In addition, sigma factors in Synechocystis may have different selectivities to differerent nucleotide sequences in the −35 and −10 regions of specific promoter sequences .
Previous crystallization studies show that the conformation of the [aTc∙Mg]2 +-bound TetR dimer is unable to bind DNA [19, 23]. As expected, no interaction was observed when TetR was incubated and injected together with aTc and MgCl2 (Fig. 4). The promoter strength under induced conditions could reveal a promoter’s ability to perform transcription (Table 1). Additionally, from the discussion in our previous study , the repressed promoter strength of the L12 promoter reached the detection limit and indicated that this promoter is fully repressed. Comparing the induced promoter strength 0.043 to the detection limit 0.022 (Table 1), we further confirm that the L12 promoter is an extremely weak promoter.
A flanking point mutation in the L09, L10, L11, and L12 promoters neither causes significant difference in the DNA breathing dynamics at the TetR binding sites nor affects the kinetic rate constants of the interactions between TetR and the respective L promoter. The leakage problem of the L09 promoter may be due to enhanced RNAP binding. These results assist in creating more tightly regulated promoters for realizing synthetic biology and metabolic engineering in biotechnological applications.
HHH designed and did the simulations. HHH and CS designed and conducted the SPR experiments, with support from UHD. PL coordinated the project. All authors were involved in the analysis and interpretation of data, and have contributed to the manuscript. All authors (HHH, CS, UHD and PL) read and approved the final manuscript.
The research leading to these results has received funding from the Swedish Energy Agency; the Knut and Alice Wallenberg Foundation (project MoSE); the Swedish Research Council; and the European Union Seventh Framework Programme (FP7/2007-2013) under grant agreement 308518 (project CyanoFactory). The simulations were performed on resources (project snic2013-1-185) provided by the Swedish National Infrastructure for Computing (SNIC) at Uppsala Multidisciplinary Center for Advanced Computational Science (UPPMAX).
Compliance with ethical guidelines
Competing interests The authors declare that they have no competing interests.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Andrianantoandro E, Basu S, Karig DK, Weiss R. Synthetic biology: new engineering rules for an emerging discipline. Mol Syst Biol. 2006;2:0028.View ArticlePubMedGoogle Scholar
- Keasling JD. Synthetic biology for synthetic chemistry. ACS Chem Biol. 2008;3:64–76.View ArticlePubMedGoogle Scholar
- Nielsen J, Keasling JD. Synergies between synthetic biology and metabolic engineering. Nat Biotechnol. 2011;29:693–5.View ArticlePubMedGoogle Scholar
- Browning DF, Busby SJ. The regulation of bacterial transcription initiation. Nat Rev Microbiol. 2004;2:57–65.View ArticlePubMedGoogle Scholar
- Rohs R, Jin X, West SM, Joshi R, Honig B, Mann R. Origins of specificity in protein-DNA recognition. Annu Rev Biochem. 2010;79:233–69.PubMed CentralView ArticlePubMedGoogle Scholar
- Del Vecchio D, Ninfa AJ, Sontag ED. Modular cell biology: retroactivity and insulation. Mol Syst Biol. 2008;4:161.PubMed CentralPubMedGoogle Scholar
- Carothers JM, Goler JA, Keasling JD. Chemical synthesis using synthetic biology. Curr Opin Biotechnol. 2009;20:498–503.View ArticlePubMedGoogle Scholar
- Huang HH, Lindblad P. Wide-dynamic-range promoters engineered for cyanobacteria. J Biol Eng. 2013;7:10.PubMed CentralView ArticlePubMedGoogle Scholar
- Alexandrov BS, Fukuyo Y, Lange M, Horikoshi N, Gelev V, Rasmussen KO, Bishop AR, Usheva A. DNA breathing dynamics distinguish binding from nonbinding consensus sites for transcription factor YY1 in cells. Nucleic Acids Res. 2012;40:10116–23.PubMed CentralView ArticlePubMedGoogle Scholar
- Bishop AR, Rasmussen KO, Usheva A, Alexandrov BS. Entropy-driven conformations controlling DNA functions. In: Kakeshita T, Fukuda T, Saxena A, Planes A, editors. Disorder and strain-induced complexity in functional materials. Berlin, Heidelberg: Springer; 2012. p. 273–292.Google Scholar
- Alexandrov BS, Gelev V, Yoo SW, Bishop AR, Rasmussen KO, Usheva A. Toward a detailed description of the thermally induced dynamics of the core promoter. PLoS Comput Biol. 2009;5:e1000313.PubMed CentralView ArticlePubMedGoogle Scholar
- Apostolaki A, Kalosakas G. Targets of DNA-binding proteins in bacterial promoter regions present enhanced probabilities for spontaneous thermal openings. Phys Biol. 2011;8:26006.View ArticleGoogle Scholar
- Nowak-Lovato K, Alexandrov LB, Banisadr A, Bauer AL, Bishop AR, Usheva A, Mu F, Hong-Geller E, Rasmussen KO, Hlavacek WS, Alexandrov BS. Binding of nucleoid-associated protein Fis to DNA is regulated by DNA breathing dynamics. PLoS Comput Biol. 2013;9:e1002881.PubMed CentralView ArticlePubMedGoogle Scholar
- Alexandrov BS, Gelev V, Yoo SW, Alexandrov LB, Fukuyo Y, Bishop AR, Rasmussen KO, Usheva A. DNA dynamics play a role as a basal transcription factor in the positioning and regulation of gene transcription initiation. Nucleic Acids Res. 2010;38:1790–5.PubMed CentralView ArticlePubMedGoogle Scholar
- Vaidyanathan VG, Xu L, Cho B. Binary and ternary binding anities between exonuclease-deficient Klenow fragment (Kf-exo−) and various arylamine DNA lesions characterized by surface plasmon resonance. Chem Res Toxicol. 2012;25:1568–70.PubMed CentralView ArticlePubMedGoogle Scholar
- Kamionka A, Bogdanska-Urbaniak J, Scholz O, Hillen W. Two mutations in the tetracycline repressor change the inducer anhydrotetracycline to a corepressor. Nucleic Acids Res. 2004;32:842–7.PubMed CentralView ArticlePubMedGoogle Scholar
- Stockley PG, Persson B. Surface plasmon resonance assays of DNA–protein interactions. In: Bishop AR, Rasmussen KØ, Usheva A, Alexandrov BS, editors. DNA–protein interactions. New York: Humana Press; 2009. p. 653–669.Google Scholar
- Majka J, Speck C. Analysis of protein–DNA interactions using surface plasmon resonance. In: Seitz H, editor. Analytics of protein–DNA interactions. Berlin, Heidelberg: Springer; 2007. p. 13–36.Google Scholar
- Orth P, Schnappinger D, Hillen W, Saenger W, Hinrichs W. Structural basis of gene regulation by the tetracycline inducible Tet repressor-operator system. Nat Struct Biol. 2000;7:215–9.View ArticlePubMedGoogle Scholar
- Hook-Barnard IG, Hinton DM. Transcription initiation by mix and match elements: flexibility for polymerase binding to bacterial promoters. Gene Regul Syst Bio. 2007;1:275–93.PubMed CentralPubMedGoogle Scholar
- Feklistov A, Darst SA. Structural basis for promoter −10 element recognition by the bacterial RNA polymerase sigma subunit. Cell. 2011;147:1257–69.PubMed CentralView ArticlePubMedGoogle Scholar
- Haugen SP, Berkmen MB, Ross W, Gaal T, Ward C, Gourse RL. rRNA promoter regulation by nonoptimal binding of sigma region 1.2: an additional recognition element for RNA polymerase. Cell. 2006;125:1069–82.View ArticlePubMedGoogle Scholar
- Orth P, Cordes F, Schnappinger D, Hillen W, Saenger W, Hinrichs W. Conformational changes of the Tet repressor induced by tetracycline trapping. J Mol Biol. 1998;279:439–47.View ArticlePubMedGoogle Scholar
- Daley SM, Kappell AD, Carrick MJ, Burnap RL. Regulation of the cyanobacterial CO2-concentrating mechanism involves internal sensing of NADP+ and α-ketogutarate levels by transcription factor CcmR. PLoS One. 2012;7:e41286.PubMed CentralView ArticlePubMedGoogle Scholar
- Kleinschmidt C, Tovar K, Hillen W, Porschke D. Dynamics of repressor-operator recognition. Tn10-encoded tetracycline resistance control. Biochemistry. 1988;27:1094–104.View ArticlePubMedGoogle Scholar
- Mekler V, Minakhin L, Severinov K. A critical role of downstream RNA polymerase-promoter interactions in the formation of initiation complex. J Biol Chem. 2011;286:22600–8.PubMed CentralView ArticlePubMedGoogle Scholar
- Imamura S, Asayama M. Sigma factors for cyanobacterial transcription. Gene Regul Syst Bio. 2009;3:65–78.PubMed CentralPubMedGoogle Scholar
- Berla BM, Saha R, Immethun CM, Maranas CD, Seok Moon T, Pakrasi HB. Synthetic biology of cyanobacteria: unique challenges and opportunities. Front Microbiol. 2013;4:246.PubMed CentralView ArticlePubMedGoogle Scholar
- Guerrero F, Carbonell V, Cossu M, Correddu D, Jones PR. Ethylene synthesis and regulated expression of recombinant protein in Synechocystis sp. PCC 6803. PLoS One. 2012;7:e50470.PubMed CentralView ArticlePubMedGoogle Scholar