Stochastic simulations of the tetracycline operon

BackgroundThe tetracycline operon is a self-regulated system. It is found naturally in bacteria where it confers resistance to antibiotic tetracycline. Because of the performance of the molecular elements of the tetracycline operon, these elements are widely used as parts of synthetic gene networks where the protein production can be efficiently turned on and off in response to the presence or the absence of tetracycline. In this paper, we investigate the dynamics of the tetracycline operon. To this end, we develop a mathematical model guided by experimental findings. Our model consists of biochemical reactions that capture the biomolecular interactions of this intriguing system. Having in mind that small biological systems are subjects to stochasticity, we use a stochastic algorithm to simulate the tetracycline operon behavior. A sensitivity analysis of two critical parameters embodied this system is also performed providing a useful understanding of the function of this system.ResultsSimulations generate a timeline of biomolecular events that confer resistance to bacteria against tetracycline. We monitor the amounts of intracellular TetR2 and TetA proteins, the two important regulatory and resistance molecules, as a function of intrecellular tetracycline. We find that lack of one of the promoters of the tetracycline operon has no influence on the total behavior of this system inferring that this promoter is not essential for Escherichia coli. Sensitivity analysis with respect to the binding strength of tetracycline to repressor and of repressor to operators suggests that these two parameters play a predominant role in the behavior of the system. The results of the simulations agree well with experimental observations such as tight repression, fast gene expression, induction with tetracycline, and small intracellular TetR2 amounts.ConclusionsComputer simulations of the tetracycline operon afford augmented insight into the interplay between its molecular components. They provide useful explanations of how the components and their interactions have evolved to best serve bacteria carrying this operon. Therefore, simulations may assist in designing novel gene network architectures consisting of tetracycline operon components.

[1]  Yiannis N. Kaznessis,et al.  Models for synthetic biology , 2007, BMC Systems Biology.

[2]  D. W. Daniels,et al.  Promoter mutations affecting divergent transcription in the Tn10 tetracycline resistance determinant. , 1985, Journal of molecular biology.

[3]  W. Hillen,et al.  Kinetic and equilibrium characterization of the Tet repressor-tetracycline complex by fluorescence measurements. Evidence for divalent metal ion requirement and energy transfer. , 1986, Journal of Molecular Biology.

[4]  A. Arkin,et al.  Stochastic mechanisms in gene expression. , 1997, Proceedings of the National Academy of Sciences of the United States of America.

[5]  J. Davies,et al.  Molecular Biology of the Cell , 1983, Bristol Medico-Chirurgical Journal.

[6]  Robert J. Brooker,et al.  Genetics: Analysis and Principles , 1998 .

[7]  W. Reznikoff,et al.  Construction of a single-copy promoter vector and its use in analysis of regulation of the transposon Tn10 tetracycline resistance determinant , 1984, Journal of bacteriology.

[8]  Howard M. Salis,et al.  Numerical simulation of stochastic gene circuits , 2005, Comput. Chem. Eng..

[9]  T J FRANKLIN,et al.  RESISTANCE OF ESCHERICHIA COLI TO TETRACYCLINES. , 1965, The Biochemical journal.

[10]  W Hillen,et al.  Differential regulation of the Tn10‐encoded tetracycline resistance genes tetA and tetR by the tandem tet operators O1 and O2. , 1988, The EMBO journal.

[11]  J. Rawlings,et al.  Approximate simulation of coupled fast and slow reactions for stochastic chemical kinetics , 2002 .

[12]  Vassilios Sotiropoulos,et al.  Multiscale models for synthetic biology , 2009, 2009 Annual International Conference of the IEEE Engineering in Medicine and Biology Society.

[13]  A. Yamaguchi,et al.  Transport of divalent cations with tetracycline as mediated by the transposon Tn10-encoded tetracycline resistance protein. , 1990, The Journal of biological chemistry.

[14]  Yiannis N. Kaznessis,et al.  Forward engineering of synthetic bio-logical and gates , 2008 .

[15]  W. Hillen,et al.  Determination of the equilibrium association constant between Tet repressor and tetracycline at limiting Mg2+ concentrations: a generally applicable method for effector-dependent high-affinity complexes. , 1991, Analytical biochemistry.

[16]  W. Hillen,et al.  Dynamics of repressor-operator recognition: the Tn10-encoded tetracycline resistance control. , 1988, Biochemistry.

[17]  An-Ping Zeng,et al.  Nonlinear Dynamics of Regulation of Bacterial trpOperon: Model Analysis of Integrated Effects of Repression, Feedback Inhibition, and Attenuation , 2002, Biotechnology progress.

[18]  A. Salyers,et al.  Bacterial resistance to tetracycline: mechanisms, transfer, and clinical significance , 1992, Clinical Microbiology Reviews.

[19]  Yiannis N. Kaznessis,et al.  SynBioSS designer: a web-based tool for the automated generation of kinetic models for synthetic biological constructs , 2010, Briefings Bioinform..

[20]  Julian I. Rood,et al.  Nomenclature for New Tetracycline Resistance Determinants , 1999, Antimicrobial Agents and Chemotherapy.

[21]  W. Hillen,et al.  Computer simulations and experimental studies of gel mobility patterns for weak and strong non-cooperative protein binding to two targets on the same DNA: application to binding of tet repressor variants to multiple and single tet operator sites. , 1991, Nucleic acids research.

[22]  A. Arkin,et al.  It's a noisy business! Genetic regulation at the nanomolar scale. , 1999, Trends in genetics : TIG.

[23]  M. Elowitz,et al.  A synthetic oscillatory network of transcriptional regulators , 2000, Nature.

[24]  C. Beck,et al.  Overproduction of transposon Tn10-encoded tetracycline resistance protein results in cell death and loss of membrane potential , 1989, Journal of bacteriology.

[25]  W. Hillen,et al.  Quantitative analysis of Tn10 Tet repressor binding to a complete set of tet operator mutants. , 1990, Nucleic acids research.

[26]  W. Saenger,et al.  Structural basis of gene regulation by the tetracycline inducible Tet repressor–operator system , 2000, Nature Structural Biology.

[27]  L. McMurry,et al.  Tetracycline resistance: Efflux, Mutations and, Other Mechanisms , 2005 .

[28]  Vassilios Sotiropoulos,et al.  Multiscale Hy3S: Hybrid stochastic simulation for supercomputers , 2006, BMC Bioinformatics.

[29]  Vassilios Sotiropoulos,et al.  Synthetic tetracycline-inducible regulatory networks: computer-aided design of dynamic phenotypes , 2007, BMC Systems Biology.

[30]  C. Gatz,et al.  Control of expression of the Tn10-encoded tetracycline resistance operon. II. Interaction of RNA polymerase and TET repressor with the tet operon regulatory region. , 1984, Journal of molecular biology.

[31]  W. Hillen,et al.  Regulated expression of heterologous genes inBacillus subtilis using the Tn10 encodedtet regulatory elements , 1990, Applied Microbiology and Biotechnology.

[32]  Marilyn Roberts,et al.  Tetracycline Antibiotics: Mode of Action, Applications, Molecular Biology, and Epidemiology of Bacterial Resistance , 2001, Microbiology and Molecular Biology Reviews.

[33]  Daniel T Gillespie,et al.  Stochastic simulation of chemical kinetics. , 2007, Annual review of physical chemistry.

[34]  D. Gillespie.  Exact Stochastic Simulation of Coupled Chemical Reactions , 1977 .

[35]  W. Hillen,et al.  The Role of the Variable Region in Tet Repressor for Inducibility by Tetracycline* , 1997, The Journal of Biological Chemistry.

[36]  D. A. Mcquarrie.  Stochastic approach to chemical kinetics , 1967, Journal of Applied Probability.

[37]  Vincent A. Fischetti,et al.  Gram-positive pathogens , 2006 .

[38]  Hermann Bujard,et al.  Tetracycline-regulated gene expression in the brain , 2000, Current Opinion in Neurobiology.

[39]  A. Arkin,et al.  Stochastic kinetic analysis of developmental pathway bifurcation in phage lambda-infected Escherichia coli cells. , 1998, Genetics.

[40]  B. Cooperman,et al.  Photoincorporation of tetracycline into Escherichia coli ribosomes , 1980, FEBS letters.

[41]  M. Gossen,et al.  Control of gene activity in higher eukaryotic cells by prokaryotic regulatory elements. , 1994, Trends in biochemical sciences.

[42]  Michail Stamatakis,et al.  Comparison of deterministic and stochastic models of the lac operon genetic network. , 2009, Biophysical journal.

[43]  J. Collins,et al.  Construction of a genetic toggle switch in Escherichia coli , 2000, Nature.

[44]  Jeffrey W. Smith,et al.  Stochastic Gene Expression in a Single Cell , .

[45]  H. Bremer.  Modulation of Chemical Composition and Other Parameters of the Cell by Growth Rate , 1999 .

[46]  Frederic M. Sapunaric,et al.  Frontiers in Antimicrobial Resistance: A Tribute to Stuart B. Levy , 2005 .

[47]  A. Kierzek,et al.  Bridging the gap between stochastic and deterministic regimes in the kinetic simulations of the biochemical reaction networks. , 2004, Biophysical journal.

[48]  C. Daub,et al.  BMC Systems Biology , 2007 .

[49]  W. Hillen,et al.  Mechanisms underlying expression of Tn10 encoded tetracycline resistance. , 1994, Annual review of microbiology.

[50]  W. Hillen,et al.  The Tn10-encoded tetR mRNA has heterogeneous 5' ends in vivo and in vitro. , 1992, Gene.

[51]  C. Gatz,et al.  Control of expression of the Tn10-encoded tetracycline resistance genes. Equilibrium and kinetic investigation of the regulatory reactions. , 1983, Journal of molecular biology.

[52]  W. Saenger,et al.  The Tetracycline Repressor-A Paradigm for a Biological Switch. , 2000, Angewandte Chemie.

[53]  Yiannis Kaznessis,et al.  Accurate hybrid stochastic simulation of a system of coupled chemical or biochemical reactions. , 2005, The Journal of chemical physics.

[54]  J. Keasling,et al.  Mathematical Model of the lac Operon: Inducer Exclusion, Catabolite Repression, and Diauxic Growth on Glucose and Lactose , 1997, Biotechnology progress.

[55]  Yiannis N Kaznessis,et al.  Model-driven designs of an oscillating gene network. , 2005, Biophysical journal.

[56]  B. Cooperman,et al.  Photoincorporation of tetracycline into Escherichia coli ribosomes. Identification of the major proteins photolabeled by native tetracycline and tetracycline photoproducts and implications for the inhibitory action of tetracycline on protein synthesis. , 1983, Biochemistry.

[57]  Michael A. Gibson,et al.  Efficient Exact Stochastic Simulation of Chemical Systems with Many Species and Many Channels , 2000 .

[58]  Konstantinos Michalodimitrakis,et al.  Noise in transcription negative feedback loops: simulation and experimental analysis , 2006, Molecular systems biology.

[59]  W Hillen,et al.  Nomenclature for tetracycline resistance determinants. , 1989, Antimicrobial agents and chemotherapy.

[60]  S. Levy,et al.  Active efflux of tetracycline encoded by four genetically different tetracycline resistance determinants in Escherichia coli. , 1980, Proceedings of the National Academy of Sciences of the United States of America.

[61]  A. R. English,et al.  Tetracycline-microbiologic, pharmacologic, and clinical evaluation. , 1954, Antibiotics & chemotherapy.

[62]  M. Strauss,et al.  Book ReviewChance and Necessity: An essay on the natural philosophy of modern biology. , 1972 .

[63]  J. Lefèvre,et al.  A new experimental approach for studying the association between RNA polymerase and the tet promoter of pBR322. , 1984, Nucleic acids research.

[64]  W. Hillen,et al.  Tetracycline analogs affecting binding to Tn10-Encoded Tet repressor trigger the same mechanism of induction. , 1996, Biochemistry.

[65]  Vassilios Sotiropoulos,et al.  SynBioSS: the synthetic biology modeling suite , 2008, Bioinform..

[66]  W. Hillen,et al.  Permeation of tetracyclines through membranes of liposomes and Escherichia coli. , 2000, European journal of biochemistry.

[67]  Patrick F. McDermott,et al.  Frontiers in antimicrobial resistance: a tribute to Stuart B. Levy. , 2005 .

[68]  M. Karp,et al.  A recombinant Escherichia coli sensor strain for the detection of tetracyclines. , 1998, Analytical chemistry.

[69]  D. Lilley,et al.  Activation of the leu-500 Promoter by a Reversed Polarity tetA Gene , 1998, The Journal of Biological Chemistry.