Dissecting the functional program of Escherichia coli promoters: the combined mode of action of Lac repressor and AraC activator

Determining noisy attractors of delayed stochastic gene regulatory networks from multiple data sources

MOTIVATION

Gene regulatory networks (GRNs) are stochastic, thus, do not have attractors, but can remain in confined regions of the state space, i.e. the 'noisy attractors', which define the cell type and phenotype.

RESULTS

We propose a gamma-Bernoulli mixture model clustering algorithm (GammaBMM), tailored for quantizing states from gamma and Bernoulli distributed data, to determine the noisy attractors of stochastic GRN. GammaBMM uses multiple data sources, naturally selects the number of states and can be extended to other parametric distributions according to the number and type of data sources available. We apply it to protein and RNA levels, and promoter occupancy state of a toggle switch and show that it can be bistable, tristable or monostable depending on its internal noise level. We show that these results are in agreement with the patterns of differentiation of model cells whose pathway choice is driven by the switch. We further apply GammaBMM to a model of the MeKS module of Bacillus subtilis, and the results match experimental data, demonstrating the usability of GammaBMM.

AVAILABILITY

Implementation software is available upon request.

Delayed stochastic model of transcription at the single nucleotide level

We present a delayed stochastic model of transcription at the single nucleotide level. The model accounts for the promoter open complex formation and includes alternative pathways to elongation, namely pausing, arrest, misincorporation and editing, pyrophosphorolysis, and premature termination. We confront the dynamics of this detailed model with a single-step multi-delayed stochastic model and with measurements of expression of a repressed gene at the single molecule level. At low expression rates both models match the experiments but, at higher rates the two models differ significantly, with consequences to cell-to-cell phenotypic variability. The alternative pathway reactions, due to, for example, causing polymerases to collide more often on the template, are the cause for the difference in dynamical behaviors. Next, we confront the model with measurements of the transcriptional dynamics at the single RNA level of an induced gene and show that RNA production, besides its bursting dynamics, also exhibits pulses (2 or more RNAs produced in intervals smaller than the smallest interval between initiations). The distribution of occurrences and amplitudes of pulses match the experimental measurements. This pulsing and the noise at the elongation stage are shown to play a role in the dynamics of a genetic switch.

Characterization of ExsA and of ExsA-dependent promoters required for expression of the Pseudomonas aeruginosa type III secretion system

Expression of the Pseudomonas aeruginosa type III secretion system (T3SS) is activated by ExsA, a member of the AraC/XylS family of transcriptional regulators. In the present study we examine the DNA-binding properties of ExsA. ExsA was purified as a histidine-tagged fusion protein (ExsA(His)) and found to be monomeric in solution. ExsA(His) specifically bound T3SS promoters with high affinity as determined by electrophoretic mobility shift assays (EMSA). For each promoter tested two distinct ExsA-DNA complexes were detected. Biochemical analyses indicate that the higher-mobility complex consists of a single ExsA(His) molecule bound to DNA while the lower-mobility complex results from the binding of two ExsA(His) molecules. DNase I protection assays demonstrate that the ExsA(His) binding site overlaps the -35 RNA polymerase binding site and extends upstream an additional approximately 34 bp. An alignment of all 10 ExsA-dependent promoters revealed a number of highly conserved nucleotides within the footprinted region. We find that most of the highly conserved nucleotides are required for transcription in vivo; EMSA-binding assays confirm that several of these nucleotides are essential determinants of ExsA(His) binding. The combined data support a model in which two ExsA(His) molecules bind adjacent sites on the promoter to activate T3SS gene transcription.

Genetic parts to program bacteria

Genetic engineering is entering a new era, where microorganisms can be programmed using synthetic constructs of DNA encoding logic and operational commands. A toolbox of modular genetic parts is being developed, comprised of cell-based environmental sensors and genetic circuits. Systems have already been designed to be interconnected with each other and interfaced with the control of cellular processes. Engineering theory will provide a predictive framework to design operational multicomponent systems. On the basis of these developments, increasingly complex cellular machines are being constructed to build specialty chemicals, weave biomaterials, and to deliver therapeutics.

In-cell NMR for protein-protein interactions (STINT-NMR)

We describe an in-cell NMR-based method for mapping the structural interactions (STINT-NMR) that underlie protein-protein complex formation. This method entails sequentially expressing two (or more) proteins within a single bacterial cell in a time-controlled manner and monitoring their interactions using in-cell NMR spectroscopy. The resulting NMR data provide a complete titration of the interaction and define structural details of the interacting surfaces at atomic resolution. Unlike the case where interacting proteins are simultaneously overexpressed in the labeled medium, in STINT-NMR the spectral complexity is minimized because only the target protein is labeled with NMR-active nuclei, which leaves the interactor protein(s) cryptic. This method can be combined with genetic and molecular screens to provide a structural foundation for proteomic studies. The protocol takes 4 d from the initial transformation of the bacterial cells to the acquisition of the NMR spectra.

Synthetic biology: new engineering rules for an emerging discipline

Synthetic biologists engineer complex artificial biological systems to investigate natural biological phenomena and for a variety of applications. We outline the basic features of synthetic biology as a new engineering discipline, covering examples from the latest literature and reflecting on the features that make it unique among all other existing engineering fields. We discuss methods for designing and constructing engineered cells with novel functions in a framework of an abstract hierarchy of biological devices, modules, cells, and multicellular systems. The classical engineering strategies of standardization, decoupling, and abstraction will have to be extended to take into account the inherent characteristics of biological devices and modules. To achieve predictability and reliability, strategies for engineering biology must include the notion of cellular context in the functional definition of devices and modules, use rational redesign and directed evolution for system optimization, and focus on accomplishing tasks using cell populations rather than individual cells. The discussion brings to light issues at the heart of designing complex living systems and provides a trajectory for future development.

Trends and challenges in enzyme technology

Several major developments took place in the field of biocatalysis over the past few years. These include the invention of directed evolution as an extremely useful method for biocatalyst improvement on the molecular level in combination with high-throughput screening systems, methods for accessing "nonculturable" biodiversity using metagenome approaches and progress in sequence-based biocatalyst discovery. In addition, new carriers and tools for immobilization of enzymes have been developed. For the synthesis of optically active compounds impressive examples using new enzymes and major progress in dynamic kinetic resolutions of racemates took place. These achievements are summarized in this review.

Model-driven designs of an oscillating gene network

The current rapid expansion of biological knowledge offers a great opportunity to rationally engineer biological systems that respond to signals such as light and chemical inducers by producing specific proteins. Turning on and off the production of proteins on demand holds great promise for creating significant biotechnological and biomedical applications. With successful stories already registered, the challenge still lies with rationally engineering gene regulatory networks which, like electronic circuits, sense inputs and generate desired outputs. From the literature, we have found kinetic and thermodynamic information describing the molecular components and interactions of the transcriptionally repressing lac, tet, and ara operons. Connecting these components in a model gene network, we determine how to change the kinetic parameters to make this normally nonperiodic system one which has well-defined oscillations. Simulating the designed lac-tet-ara gene network using a hybrid stochastic-discrete and stochastic-continuous algorithm, we seek to elucidate the relationship between the strength and type of specific connections in the gene network and the oscillatory nature of the protein product. Modeling the molecular components of the gene network allows the simulation to capture the dynamics of the real biological system. Analyzing the effect of modifications at this level provides the ability to predict how changes to experimental systems will alter the network behavior, while saving the time and expense of trial and error experimental modifications.

RNA polymerase structure and function at lac operon

Transcription of E. coli lac operon by RNA polymerase (RNAP) is a classic example of how the basic functions of this enzyme, specifically the ability to recognize/bind promoters, melt the DNA and initiate RNA synthesis, is positively regulated by transcription activators, such as cyclic AMP-receptor protein, CRP, and negatively regulated by lac-repressor, LacI. In this review, we discuss the recent progress in structural and biochemical studies of RNAP and its binary and ternary complexes with CRP and lac promoter. With structural information now available for RNAP and models of binary and ternary elongation complexes, the interaction between these factors and RNAP can be modeled, and possible molecular mechanisms of their action can be inferred.

Advances in synthetic biology: on the path from prototypes to applications

Synthetic biology combines knowledge from various disciplines including molecular biology, engineering, mathematics and physics to design and build novel proteins, genetic circuits and metabolic networks. Early efforts aimed at altering the behavior of individual elements have now evolved to focus on the construction of complex networks in single-cell and multicellular systems. Recent achievements include the development of sophisticated non-native behaviors such as bi-stability, oscillations, proteins customized for biosensing, optimized drug synthesis and programmed spatial pattern formation. The de novo construction of such systems offers valuable quantitative insight into naturally occurring information processing activities. Furthermore, as the techniques for system design, synthesis and optimization mature, we will witness a rapid growth in the capabilities of synthetic systems with a wide-range of applications.

RNA dynamics in live Escherichia coli cells

We describe a method for tracking RNA molecules in Escherichia coli that is sensitive to single copies of mRNA, and, using the method, we find that individual molecules can be followed for many hours in living cells. We observe distinct characteristic dynamics of RNA molecules, all consistent with the known life history of RNA in prokaryotes: localized motion consistent with the Brownian motion of an RNA polymer tethered to its template DNA, free diffusion, and a few examples of polymer chain dynamics that appear to be a combination of chain fluctuation and chain elongation attributable to RNA transcription. We also quantify some of the dynamics, such as width of the displacement distribution, diffusion coefficient, chain elongation rate, and distribution of molecule numbers, and compare them with known biophysical parameters of the E. coli system.

Expression of the chitobiose operon of Escherichia coli is regulated by three transcription factors: NagC, ChbR and CAP

The chitobiose operon, chbBCARFG, encodes genes for the transport and degradation of the N-acetylglucosamine disaccharide, chitobiose. Chitobiose is transported by the phosphotransferase system (PTS) producing chitobiose-6P which is hydrolysed to GlcNAc-6P by the chbF gene product and then further degraded by the nagBA gene products. Expression of the chb operon is repressed by NagC, which regulates genes involved in amino sugar metabolism. The inducer for NagC is GlcNAc-6P. NagC binds to two sites separated by 115 bp and the transcription start point of the chb operon lies within the downstream NagC operator. In addition the chb operon encodes its own specific regulator, ChbR, an AraC-type dual repressor-activator, which binds to two direct repeats of 19 bp located between the two NagC sites. ChbR is necessary for transcription activation in the presence of chitobiose in vivo. Induction of the operon also requires CAP, which binds to a site upstream of the ChbR repeats. In the absence of chitobiose both NagC and ChbR act as repressors. Together these three factors cooperate in switching chb expression from the repressed to the activated state. The need for two specific inducing signals, one for ChbR to activate the expression of the operon and a second for NagC to relieve its repression, ensure that the chb operon is only induced when there is sufficient flux through the combined chb-nag metabolic pathway to activate expression of both the chb and nag operons.

The davDT operon of Pseudomonas putida, involved in lysine catabolism, is induced in response to the pathway intermediate delta-aminovaleric acid

Pseudomonas putida KT2440 is a soil microorganism that attaches to seeds and efficiently colonizes the plant's rhizosphere. Lysine is one of the major compounds in root exudates, and P. putida KT2440 uses this amino acid as a source of carbon, nitrogen, and energy. Lysine is channeled to delta-aminovaleric acid and then further degraded to glutaric acid via the action of the davDT gene products. We show that the davDT genes form an operon transcribed from a single sigma70-dependent promoter. The relatively high level of basal expression from the davD promoter increased about fourfold in response to the addition of exogenous lysine to the culture medium. However, the true inducer of this operon seems to be delta-aminovaleric acid because in a mutant unable to metabolize lysine to delta-aminovaleric acid, this compound, but not lysine, acted as an effector. Effective induction of the P. putida P(davD) promoter by exogenously added lysine requires efficient uptake of this amino acid, which seems to proceed by at least two uptake systems for basic amino acids that belong to the superfamily of ABC transporters. Mutants in these ABC uptake systems retained basal expression from the davD promoter but exhibited lower induction levels in response to exogenous lysine than the wild-type strain.

Direct stimulation of the lambdapaQ promoter by the transcription effector guanosine-3',5'-(bis)pyrophosphate in a defined in vitro system

The bacterial response to nutritional deprivation, called the stringent response, results in the introduction of the specific nucleotide guanosine-3',5'-(bis) pyrophosphate (ppGpp). This nucleotide interacts with RNA polymerase and alters its action so that transcription from certain promoters is inhibited, whereas transcription from others seems to be activated. The exact mechanism of transcriptional stimulation by ppGpp in vivo remains unknown. A passive control model has been proposed according to which transcription inhibition during the stringent response at several very active promoters, like those for rRNA and tRNA genes, makes more free RNA polymerase (RNAP) molecules available for transcription at promoters with weak binding affinities for RNAP, thus leading to their passive activation. Among promoters whose transcription is activated by ppGpp in vivo is the histidine operon promoter (hisGp). However, in vitro it is only possible to demonstrate this effect in a coupled transcription-translation system. Here we demonstrate, using another in vivo ppGpp-stimulated promoter, the phage lambdapaQ promoter, that activation by ppGpp in a defined in vitro system is direct. A systematic study of ppGpp effects on the stimulation of paQ revealed that, as in the case of promoters inhibited by this nucleotide, ppGpp decreases the half-life of paQ open complexes. Our results also indicate that the equilibrium binding affinity of RNA polymerase to paQ seems not to be affected in the presence of ppGpp. Our data indicate that the mechanism underlying ppGpp stimulation of paQ is due to an increased rate of productive open complex formation.