William Mather

A fast, robust and tunable synthetic gene oscillator

One defining goal of synthetic biology is the development of engineering-based approaches that enable the construction of gene-regulatory networks according to ‘design specifications’ generated from computational modelling. This approach provides a systematic framework for exploring how a given regulatory network generates a particular phenotypic behaviour. Several fundamental gene circuits have been developed using this approach, including toggle switches and oscillators, and these have been applied in new contexts such as triggered biofilm development and cellular population control. Here we describe an engineered genetic oscillator in Escherichia coli that is fast, robust and persistent, with tunable oscillatory periods as fast as 13 min. The oscillator was designed using a previously modelled network architecture comprising linked positive and negative feedback loops. Using a microfluidic platform tailored for single-cell microscopy, we precisely control environmental conditions and monitor oscillations in individual cells through multiple cycles. Experiments reveal remarkable robustness and persistence of oscillations in the designed circuit; almost every cell exhibited large-amplitude fluorescence oscillations throughout observation runs. The oscillatory period can be tuned by altering inducer levels, temperature and the media source. Computational modelling demonstrates that the key design principle for constructing a robust oscillator is a time delay in the negative feedback loop, which can mechanistically arise from the cascade of cellular processes involved in forming a functional transcription factor. The positive feedback loop increases the robustness of the oscillations and allows for greater tunability. Examination of our refined model suggested the existence of a simplified oscillator design without positive feedback, and we construct an oscillator strain confirming this computational prediction.

Queueing up for enzymatic processing: correlated signaling through coupled degradation.

High‐throughput technologies have led to the generation of complex wiring diagrams as a post‐sequencing paradigm for depicting the interactions between vast and diverse cellular species. While these diagrams are useful for analyzing biological systems on a large scale, a detailed understanding of the molecular mechanisms that underlie the observed network connections is critical for the further development of systems and synthetic biology. Here, we use queueing theory to investigate how ‘waiting lines’ can lead to correlations between protein ‘customers’ that are coupled solely through a downstream set of enzymatic ‘servers’. Using the E. coli ClpXP degradation machine as a model processing system, we observe significant cross‐talk between two networks that are indirectly coupled through a common set of processors. We further illustrate the implications of enzymatic queueing using a synthetic biology application, in which two independent synthetic networks demonstrate synchronized behavior when common ClpXP machinery is overburdened. Our results demonstrate that such post‐translational processes can lead to dynamic connections in cellular networks and may provide a mechanistic understanding of existing but currently inexplicable links.

Translational Cross Talk in Gene Networks

It has been shown experimentally that competition for limited translational resources by upstream mRNAs can lead to an anticorrelation between protein counts. Here, we investigate a stochastic model for this phenomenon, in which gene transcripts of different types compete for a finite pool of ribosomes. Throughout, we utilize concepts from the theory of multiclass queues to describe a qualitative shift in protein count statistics as the system transitions from being underloaded (ribosomes exceed transcripts in number) to being overloaded (transcripts exceed ribosomes in number). The exact analytical solution of a simplified stochastic model, in which the numbers of competing mRNAs and ribosomes are fixed, exhibits weak positive correlations between steady-state protein counts when total transcript count slightly exceeds ribosome count, whereas the solution can exhibit strong negative correlations when total transcript count significantly exceeds ribosome count. Extending this analysis, we find approximate but reasonably accurate solutions for a more realistic model, in which abundances of mRNAs and ribosomes are allowed to fluctuate randomly. Here, ribosomal fluctuations contribute positively and mRNA fluctuations contribute negatively to correlations, and when mRNA fluctuations dominate ribosomal fluctuations, a strong anticorrelation extremum reliably occurs near the transition from the underloaded to the overloaded regime.

Correlation resonance generated by coupled enzymatic processing.

A major challenge for systems biology is to deduce the molecular interactions that underlie correlations observed between concentrations of different intracellular molecules. Although direct explanations such as coupled transcription or direct protein-protein interactions are often considered, potential indirect sources of coupling have received much less attention. Here we show how correlations can arise generically from a posttranslational coupling mechanism involving the processing of multiple protein species by a common enzyme. By observing a connection between a stochastic model and a multiclass queue, we obtain a closed form expression for the steady-state distribution of the numbers of molecules of each protein species. Upon deriving explicit analytic expressions for moments and correlations associated with this distribution, we discover a striking phenomenon that we call correlation resonance: for small dilution rate, correlations peak near the balance-point where the total rate of influx of proteins into the system is equal to the maximum processing capacity of the enzyme. Given the limited number of many important catalytic molecules, our results may lead to new insights into the origin of correlated behavior on a global scale.

Buckling instability in ordered bacterial colonies

Bacterial colonies often exhibit complex spatio-temporal organization. This collective behavior is affected by a multitude of factors ranging from the properties of individual cells (shape, motility, membrane structure) to chemotaxis and other means of cell-cell communication. One of the important but often overlooked mechanisms of spatio-temporal organization is direct mechanical contact among cells in dense colonies such as biofilms. While in natural habitats all these different mechanisms and factors act in concert, one can use laboratory cell cultures to study certain mechanisms in isolation. Recent work demonstrated that growth and ensuing expansion flow of rod-like bacteria Escherichia coli in confined environments leads to orientation of cells along the flow direction and thus to ordering of cells. However, the cell orientational ordering remained imperfect. In this paper we study one mechanism responsible for the persistence of disorder in growing cell populations. We demonstrate experimentally that a growing colony of nematically ordered cells is prone to the buckling instability. Our theoretical analysis and discrete-element simulations suggest that the nature of this instability is related to the anisotropy of the stress tensor in the ordered cell colony.

A queueing approach to multi-site enzyme kinetics

Multi-site enzymes, defined as where multiple substrate molecules can bind simultaneously to the same enzyme molecule, play a key role in a number of biological networks, with the Escherichia coli protease ClpXP a well-studied example. These enzymes can form a low latency ‘waiting line’ of substrate to the enzymes catalytic core, such that the enzyme molecule can continue to collect substrate even when the catalytic core is occupied. To understand multi-site enzyme kinetics, we study a discrete stochastic model that includes a single catalytic core fed by a fixed number of substrate binding sites. A natural queueing systems analogy is found to provide substantial insight into the dynamics of the model. From this, we derive exact results for the probability distribution of the enzyme configuration and for the distribution of substrate departure times in the case of identical but distinguishable classes of substrate molecules. Comments are also provided for the case when different classes of substrate molecules are not processed identically.

Marching along to an Offbeat Drum: Entrainment of Synthetic Gene Oscillators by a Noisy Stimulus.

Modulation of biological oscillations by stimuli lies at the root of many phenomena, including maintenance of circadian rhythms, propagation of neural signals, and somitogenesis. While it is well established that regular periodic modulation can entrain an oscillator, an aperiodic (noisy) modulation can also robustly entrain oscillations. This latter scenario may describe, for instance, the effect of irregular weather patterns on circadian rhythms, or why irregular neural stimuli can still reliably transmit information. A synthetic gene oscillator approach has already proven to be useful in understanding the entrainment of biological oscillators by periodic signaling, mimicking the entrainment of a number of noisy oscillating systems. We similarly seek to use synthetic biology as a platform to understand how aperiodic signals can strongly correlate the behavior of cells. This study should lead to a deeper understanding of how fluctuations in our environment and even within our body may promote substantial synchrony among our cells. Specifically, we investigate experimentally and theoretically the entrainment of a synthetic gene oscillator in E. coli by a noisy stimulus. This phenomenon was experimentally studied and verified by a combination of microfluidics and microscopy using the real synthetic circuit. Stochastic simulation of an associated model further supports that the synthetic gene oscillator can be strongly entrained by aperiodic signals, especially telegraph noise. Finally, widespread applicability of aperiodic entrainment beyond the synthetic gene oscillator is supported by results derived from both a model for a natural oscillator in D. discoideum and a model for predator-prey oscillations.

Synchronized DNA cycling across a bacterial population

A defining goal of synthetic biology is to engineer cells to coordinate tasks that often require precise temporal modulation of gene expression. Although a variety of relatively small gene circuits have been constructed and characterized, their logical combination into larger networks remains a central challenge. This is due primarily to the lack of compatible and orthogonal elements for predictable dynamic control of gene expression. As an alternative approach to promoter-level regulation, we explored the use of DNA copy number as a circuit control element. We engineered colony-wide DNA cycling in Escherichia coli in the form of plasmid copy number oscillations via a modular design that can be readily adapted for use with other gene circuitry. Copy number modulation is a generalizable principle that adds a layer of control to synthetic gene circuits, allowing dynamic regulation of circuit elements without requiring specially engineered promoters.

The spatiotemporal system dynamics of acquired resistance in an engineered microecology

Great strides have been made in the understanding of complex networks; however, our understanding of natural microecologies is limited. Modelling of complex natural ecological systems has allowed for new findings, but these models typically ignore the constant evolution of species. Due to the complexity of natural systems, unanticipated interactions may lead to erroneous conclusions concerning the role of specific molecular components. To address this, we use a synthetic system to understand the spatiotemporal dynamics of growth and to study acquired resistance in vivo. Our system differs from earlier synthetic systems in that it focuses on the evolution of a microecology from a killer-prey relationship to coexistence using two different non-motile Escherichia coli strains. Using empirical data, we developed the first ecological model emphasising the concept of the constant evolution of species, where the survival of the prey species is dependent on location (distance from the killer) or the evolution of resistance. Our simple model, when expanded to complex microecological association studies under varied spatial and nutrient backgrounds may help to understand the complex relationships between multiple species in intricate natural ecological networks. This type of microecological study has become increasingly important, especially with the emergence of antibiotic-resistant pathogens.

Entrainment of a Bacterial Synthetic Gene Oscillator through Proteolytic Queueing.

Internal chemical oscillators (chemical clocks) direct the behavior of numerous biological systems, and maintenance of a given period and phase among many such oscillators may be important for their proper function. However, both environmental variability and fundamental molecular noise can cause biochemical oscillators to lose coherence. One solution to maintaining coherence is entrainment, where an external signal provides a cue that resets the phase of the oscillators. In this work, we study the entrainment of gene networks by a queueing interaction established by competition between proteins for a common proteolytic pathway. Principles of queueing entrainment are investigated for an established synthetic oscillator in Escherichia coli. We first explore this theoretically using a standard chemical reaction network model and a map-based model, both of which suggest that queueing entrainment can be achieved through pulsatile production of an additional protein competing for a common degradation pathway with the oscillator proteins. We then use a combination of microfluidics and fluorescence microscopy to verify that pulse trains modulating the production rate of a fluorescent protein targeted to the same protease (ClpXP) as the synthetic oscillator can entrain the oscillator.