To elucidate the regulatory dynamics of the gene expression activation and inactivation, an in silico biochemical model of the lac circuit in Escherichia coli was used to evaluate the transcription rates that yield the steady-state mRNA production in active and inactive states of the lac circuit. This result can be used in synthetic biology applications to understand the limits of the genetic synthesis. Since most genetic networks involve many interconnected components with positive and negative feedback control, intuitive understanding of their dynamics is often difficult to obtain. Although the kinetic model of the lac circuit considered involves only a single positive feedback, the developed computational framework can be used to evaluate supported ranges of other reaction rates in genetic circuits with more complex regulatory networks. More specifically, the inducible lac gene switch in E. coli is regulated by unbinding and binding of the inducer–repressor complexes to or from the DNA operator to switch the gene expression on and off. The dependency of mRNA production at steady state on different transcription rates and the repressor complexes has been studied by computer simulations in the Lattice Microbe software. Provided that the lac circuit is in active state, the transcription rate is independent of the inducer–repressor complexes present in the cell. In inactive state, the transcription rate is dependent on the specific inducer–repressor complex bound to the operator that inactivates the gene expression. We found that the repressor complex with the largest affinity to the operator yields the smallest range of the feasible transcription rates to yield the steady state while the lac circuit is in inactive state. In contrast, the steady state in active state can be obtained for any value of the transcription rate.
Get full access to this article
View all access options for this article.
References
1.
AbhyudaiS., VargasC.A., and KarmakarR. (2013). Stochastic analysis and inference of a two-state genetic promoter model. In American Control Conference (ACC). IEEE, Washington, DC.
2.
AnanyaR., BellecourtM.J., and LandickR. (2016). Mechanisms of bacterial transcription termination: All good things must end. Annu. Rev. Biochem., 85, 319–347.
3.
AulanderS., and FusseneggerM. (2013). From gene switches to mammalian designer cells: Present and future prospects. Trends Biotechnol. 31, 155–168.
4.
BakshiS., SiryapornA., GoulianM., et al. (2012). Superresolution imaging of ribosomes and RNA polymerase in live Escherichia coli cells. Mol. Microbiol., 85, 21–38.
5.
BasanM., HuiS., OkanoH., et al. (2015). Overflow metabolism in Escherichia coli results from efficient proteome allocation. Nature, 528, 99–104.
6.
BernsteinJ.A., LinP., CohenS.N., et al. (2004). Global analysis of Escherichia coli RNA degradosome function using DNA microarrays. Proc. Natl. Acad. Sci. U. S. A., 101, 2758–2763.
7.
BridgesD., FraserM.E., and MoorheadG.B. (2005). Cyclic nucleotide binding proteins in the Arabidopsis thaliana and Oryza sativa genomes. BMC Bioinformatics, 6, 6.
8.
BurgessR.R. (1969). Separation and Characterization of the Subunits of Ribonucleic Acid Polymerase. J. Biol. Chem., 244, 6168–6176.
ChappellJ., TakahashiM.K., MeyerS., et al. (2013). The centrality of RNA for engineering gene expression. Biotechnol. J., 8, 1379–1395.
11.
ChappellJ., WestbrookA., VerosloffM., et al. (2017). Computational design of small transcription activating RNAs for versatile and dynamic gene regulation. Nat. Commun., 8, 1051.
12.
ChenJ., DarstS.A., and ThirumalaiD. (2010). Promoter melting triggered by bacterial RNA polymerase occurs in three steps. Proc. Natl. Acad. Sci. U. S. A., 107, 12523–12528.
13.
ChenY., HoJ.M., ShisD.L., et al. (2017). Tuning the dynamic range of bacterial promoters regulated by ligand-inducible transcription factors. Nat. Commun., 9, 64.
14.
ChungS.Y., LernerE., JinY., et al. (2017). The effect of macromolecular crowding on singleround transcription by E. coli RNA polymerase. Bio archive J.https://doi.org/10.1038/s41467-017-02473-5 (bioRxiv).
15.
CoxC.D., McCollumJ.M., AllenM.S., et al. (2008). Using noise to probe and Characterize gene circuits. Proc. Natl. Acad. Sci. U. S. A., 105, 10809–10814.
16.
DaberR., StayrookS., RosenbergA., et al. (2007). Structural analysis of lac repressor bound to allosteric effectors. J. Mol. Biol., 370, 609–619.
17.
DavisC.A., CappM.W., RecordM.T., et al. (2005). The effects of upstream DNA on open complex formation by Escherichia coli RNA polymerase. Proc. Natl. Acad. Sci. U. S. A., 102, 285–290.
18.
DeutscherM.P. (2006). Degradation of RNA in bacteria: Comparison of mRNA and stable RNA. Nucleic Acids Res. 34, 659–666.
19.
ErbanR., ChapmanJ., and MainiP. (2007). A Practical Guide to Stochastic Simulations of Reaction-Diffusion Processes. Cornell University Library, pp. 1–35. (http://arxiv.org/abs/0704-1908).
20.
EsmaeiliA., DavisonT., WuA., et al. (2015). PROKARYO: An illustrative and interactive computational model of the lactose operon in the bacterium Escherichia coli. BMC Bioinformatics, 16, 311.
21.
EsyuninaD., TurtolaM., PupovD., et al. (2016). Lineage-specific variations in the trigger loop modulate RNA proofreading by bacterial RNA polymerases. Nucleic Acids Res. 44, 1298–1308.
22.
FuchsG., VoichekY., BenjaminS., et al. (2014). 4sUDRB-seq: Measuring genomewide transcriptional elongation rates and initiation frequencies within cells. BMC Genome Biol. 15, R69.
23.
GardnerT.S., CantorC.R., and CollinsJ.J. (2000). Construction of a genetic toggle switch in Escherichia coli. Nature, 403, 339–342.
24.
JacobF., and MonodJ. (1961). Genetic regulatoty mechanism in the synthesis of proteins. J. Mol. Biol., 3, 318–356.
25.
JainC. (2002). Degradation of mRNA in Escherichia coli. IUBMB Life, 54, 315–321.
26.
KaliskyT., DekelE., and AlonU. (2007). Cost-benefit theory and optimal design of gene regulation functions. Phys. Biol., 4, 229–245.
27.
KhalilA.S., and OllinsJ.J. (2010). Synthetic biology: Applications come of age. Nat. Rev. Genet., 11, 367–379.
28.
KobayashiH., KaernM., ArakiM., et al. (2004). Programmable cells: Interfacing natural and engineered gene networks. Proc. Natl. Acad. Sci. U. S. A., 101, 8414–8419.
29.
KondoT., StrayerC.A., KulkarniR.D., et al. (1993). Circadian rhythms in prokaryotes: Luciferase as a reporter of circadian gene expression in cyanobacteria. Proc. Natl. Acad. Sci. U. S. A., 90, 5672–5676.
30.
KushnerS.R. (2002). mRNA decay in Escherichia coli comes of age. J. Bacteriol., 187, 4658–4665.
31.
LeeD.H., BaeJ.E., LeeJ.H., et al. (2010). Quantitative detection of residual E. coli host cell DNA by real-time PCR. J. Microbiol. Biotechnol., 20, 1463–1470.
32.
LeeD.J., MinchinS.D., and BusbyS.J. (2012). Activating transcription in bacteria. Annu. Rev. Microbiol., 66, 125–152.
33.
LeiS.Q., and ArkinA.P. (2014). A versatile framework for microbial engineering using synthetic noncoding RNAs. Nat. Rev. Microbiol., 12, 341–354.
MaH-W., BuerJ., and ZengA.-P. (2004). Hierarchical structure and modules in the Escherichia coli transcriptional regulatory network revealed by a new top-down approach. BMC Bioinformatics, 5, 199.
36.
McClureW.R. (1980). Rate-limiting steps in RNA chain initiation. Proc. Natl. Acad. Sci. U. S. A., 77, 5634–5638.
37.
McDowellJ.C., RobertsJ.W., JinD.J., et al. (1994). Determination of intrinsic transcription termination efficiency by RNA polymerase elongation rate. Science, 266, 822–825.
38.
McKeagueM., WongR.S., and SmolkeC.D. (2016). Opportunities in the design and application of RNA for gene expression control. Nucleic Acids Res. 44, 2987–2999.
39.
MejiaY.X., NudlerE., and BustamanteaC. (2014). Trigger loop folding determines transcription rate of Escherichia coli's RNA polymerase. Proc. Natl. Acad. Sci. U. S. A., 112, 743–748.
40.
MeklerV., KortkhonjiaE., MukhopadhyayJ., et al. (2002). Structural organization of bacterial RNA polymerase holoenzyme and the RNA polymerase-promoter open complex. Cell, 108, 599–614.
41.
MettetalJ., MuzzeyD., PedrazaJ.M., et al. (2005). Predicting stochastic gene expression dynamics in single cells. Proc. Natl. Acad. Sci. U. S. A., 103, 7304–7309.
42.
MunroP.D., AckersG.K., and ShearwinK.E. (2016). Aspects of protein–DNA interactions: A review of quantitative thermodynamic theory for modelling synthetic circuits utilising LacI and CI repressors, IPTG and the reporter gene lacZ. Biophys. Rev., 8, 331–345.
43.
NathK., and KochA.L. (1970). Protein degradation in Escherichia coli. J. Biol. Chem., 245, 2889–2900.
44.
OehlerS., AlbertiS., and Muller-HillB. (2006). Induction of the lac promoter in the absence of DNA loops and the stoichiometry of induction. Nucleic Acids Res. 34, 606–612.
45.
OehlerS., EismannE.R., KramerH., et al. (1990). The three operators of the lac operon cooperate in repression. EMBO J. 9, 973–979.
46.
ProshkinS., RahmouniA.R., MironovA., et al. (2010). Cooperation between translating ribosomes and RNA polymerase in transcription elongation. Science, 328, 504–508.
47.
PtashneM., and GannA. (1997). Transcriptional activation by recruitment. Nature, 386, 569–577.
48.
QayyumM.Z., DeyD., and SenR. (2016). Transcription elongation factor NusA is a general antagonist of Rho-dependent termination in Escherichia coli. J. Biol. Chem., 291, 8090–8108.
49.
RayJ.C., WckersheimM.L., JalihalA.P., et al. (2016). Cellular growth arrest and persistence from enzyme saturation. PLoS Comput. Biol. 12, e1004825.
50.
RobertsE., MagisA., OrtizJ.O., et al. (2011). Noise contributions in an inducible genetic switch: A whole-cell simulation study. PLoS Comput. Biol., 7, 1–21.
51.
RobertsE., StoneJ.E., and Luthey-SchultenZ. (2013). Lattice microbes: High-performance stochastic simulation method for the reaction-diffusion master equation. J. Comput. Chem., 34, 245–255.
52.
ReznikoffW.S., WinterR.B., and HurleyC.K. (1974). The location of the repressor binding sites in the lac operon. Proc. Nat. Acad. U. S. A., 71, 2314–2318.
53.
RussellP.J., et al. (2008). Biology: The Dynamic Science, 1st ed. Yolanda Cossio. Pacific Grove, CA, USA.
54.
SchlaxP.J., CappM.W., and RecordM.T. (1995). Inhibition of transcription initiation by lac repressor. J. Mol. Biol., 245, 331–350.
55.
SharabianiM.T., SiermalM., LehtinenT., et al. (2005). Dynamic covariation between gene expression and proteome characteristics. BMC Bioinformatics, 6, 215.
56.
SinghJ., and PadgettR.A. (2009). Rates of in situ transcription and splicing in large human genes. Nat. Struct. Mol. Biol., 16, 1128–1133.
57.
SkinnerS.O., SepulvedaL.A., XuH., et al. (2013). Measuring mRNA copy-number in individual Escherichia coli cells using single-molecule fluorescent in situ hybridization (smFISH). Nat. Protoc., 8, 1100–1113.
58.
SmithK.C., and HanawaltP.C. (1969). Molecular Photobiology: Inactivation and Recovery. Academic Press, Inc. (London) LTD, United Kingdom.
59.
SoL., GhoshA., ZongC., et al. (2011). General properties of transcriptional time series in Escherichia coli. Nat. Genet., 43, 554–560.
60.
SolomonE.P., BergL.R., MartinD.W. (2005). Biology, 7th ed. Thomson Brooks/Cole. Pacific Grove, CA, USA.
61.
StamatakisM., and MantzarisN.V. (2009). Comparison of deterministic and stochastic models of lac operon genetic network. Biophys. J., 96, 887–906.
62.
StraneyD.C., and CrothersD.M. (1985). Intermediates in transcription initiation from the E. coli lac UV5 Promoter. Cell, 43, 449–459.
63.
Swint-KruseL. (2004). Using networks to identify fine structural differences between functionally distinct protein states. J. Biochem., 43, 10886–10895.
64.
TangG-Q., RoyR., BandwarR.P., et al. (2009). Real-time observation of the transition from transcription initiation to elongation of the RNA polymerase. Proc. Natl. Acad. Sci. U. S. A., 106, 22175–22180.
65.
Tolic-NorrelykkeS.F., EnghA.M., LandickR., et al. (2004). Diversity in the rates of transcript elongation by single RNA polymerase molecules. J. Biol. Chem., 279, 3292–3299.
66.
TroppB.E. (2010). Molecular Biology: Genes to Proteins, 4th ed. Jones and Bartlett Learning, Canada.
67.
TroppB.E. (2011). Molecular Biology: Genes to Proteins, 3rd ed. Jones and Bartlett Learning, Canada.
68.
VanziF., BroggioC., SacconiL., et al. (2006). Lac repressor hinge flexibility andDNAlooping: Single molecule kinetics by tethered particle motion. Nucleic Acids Res. 34, 3409–3420.
69.
WangR., JinG., ZhangX., et al. (2009). Modeling post-transcriptional regulation activity of small non-coding RNAs in Escherichia coli. BMC Bioinformatics, 10(Suppl 4), S6.
70.
WashburnR.S., and GottesmanM.E. (2015). Regulation of transcription elongation and termination. Biomolecules, 5, 1063–1078.
71.
XuJ., LiuS., ChenM., et al. (2011). Altering residues N125 and D149 impacts sugar effector binding and allosteric parameters in Escherichia coli lactose repressor. Biochem. J., 50, 9002–9013.
72.
XueQ., and YeungE.S. (1995). Differences in the chemical reactivity of individual molecules of an enzyme. Nature, 373, 681–683.
73.
YagilG., and YagilE. (1971). On the relation between effector concentration and the rate of induced enzyme synthesis. Biophys. J., 11, 11–27.
74.
YangX., ChenB., CaiY., et al. (2015). Understanding the lac operon with GeneAct. J. Comput. Biol. Drug Des., 8, 168–188.
75.
ZhuangX., BartleyL.E., BabcockH.P., et al. (2000). A single-molecule study of RNA catalysis and folding. Science, 288, 2048–2051.
76.
ZuoY., and SteitzT.A. (2016). A structure-based kinetic model of transcription. Transcription, 8, 1–8.
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
