Bridging the gap between mechanistic biological models and machine learning surrogates.

Mechanistic models have been used for centuries to describe complex interconnected processes, including biological ones. As the scope of these models has widened, so have their computational demands. This complexity can limit their suitability when running many simulations or when real-time results are required. Surrogate machine learning (ML) models can be used to approximate the behaviour of complex mechanistic models, and once built, their computational demands are several orders of magnitude lower. This paper provides an overview of the relevant literature, both from an applicability and a theoretical perspective. For the latter, the paper focuses on the design and training of the underlying ML models. Application-wise, we show how ML surrogates have been used to approximate different mechanistic models. We present a perspective on how these approaches can be applied to models representing biological processes with potential industrial applications (e.g., metabolism and whole-cell modelling) and show why surrogate ML models may hold the key to making the simulation of complex biological systems possible using a typical desktop computer.

Design and Assembly of Multilevel Transcriptional and Translational Regulators for Stringent Control of Gene Expression.

Precise control of gene expression is crucial when reprogramming the behavior of living cells. However, common inducible systems often lack the ability to stringently control gene expression due to the use of a single type of regulator that can be susceptible to unavoidable biomolecular fluctuations. In contrast, multilevel controllers (MLCs) employ several forms of regulation simultaneously to overcome this issue, ensuring a reduced basal expression while minimally affecting the maximum induced expression level that can be achieved. Here, we show how our publicly available genetic toolkit can be used to simplify the assembly of MLCs for the stringent control of gene expression. We demonstrate how new compatible parts can be designed and explain the rapid end-to-end assembly procedure.

Towards an engineering theory of evolution.

Biological technologies are fundamentally unlike any other because biology evolves. Bioengineering therefore requires novel design methodologies with evolution at their core. Knowledge about evolution is currently applied to the design of biosystems ad hoc. Unless we have an engineering theory of evolution, we will neither be able to meet evolution's potential as an engineering tool, nor understand or limit its unintended consequences for our biological designs. Here, we propose the evotype as a helpful concept for engineering the evolutionary potential of biosystems, or other self-adaptive technologies, potentially beyond the realm of biology.

Cheetah: A Computational Toolkit for Cybergenetic Control.

Advances in microscopy, microfluidics, and optogenetics enable single-cell monitoring and environmental regulation and offer the means to control cellular phenotypes. The development of such systems is challenging and often results in bespoke setups that hinder reproducibility. To address this, we introduce Cheetah, a flexible computational toolkit that simplifies the integration of real-time microscopy analysis with algorithms for cellular control. Central to the platform is an image segmentation system based on the versatile U-Net convolutional neural network. This is supplemented with functionality to robustly count, characterize, and control cells over time. We demonstrate Cheetah's core capabilities by analyzing long-term bacterial and mammalian cell growth and by dynamically controlling protein expression in mammalian cells. In all cases, Cheetah's segmentation accuracy exceeds that of a commonly used thresholding-based method, allowing for more accurate control signals to be generated. Availability of this easy-to-use platform will make control engineering techniques more accessible and offer new ways to probe and manipulate living cells.

A Centrifuge-Based Method for Identifying Novel Genetic Traits That Affect Root-Substrate Adhesion in Arabidopsis thaliana.

The physical presence of roots and the compounds they release affect the cohesion between roots and their environment. However, the plant traits that are important for these interactions are unknown and most methods that quantify the contributions of these traits are time-intensive and require specialist equipment and complex substrates. Our lab developed an inexpensive, high-throughput phenotyping assay that quantifies root-substrate adhesion in Arabidopsis thaliana. We now report that this method has high sensitivity and versatility for identifying different types of traits affecting root-substrate adhesion including root hair morphology, vesicle trafficking pathways, and root exudate composition. We describe a practical protocol for conducting this assay and introduce its use in a forward genetic screen to identify novel genes affecting root-substrate interactions. This assay is a powerful tool for identifying and quantifying genetic contributions to cohesion between roots and their environment.

Harnessing the central dogma for stringent multi-level control of gene expression.

Strictly controlled inducible gene expression is crucial when engineering biological systems where even tiny amounts of a protein have a large impact on function or host cell viability. In these cases, leaky protein production must be avoided, but without affecting the achievable range of expression. Here, we demonstrate how the central dogma offers a simple solution to this challenge. By simultaneously regulating transcription and translation, we show how basal expression of an inducible system can be reduced, with little impact on the maximum expression rate. Using this approach, we create several stringent expression systems displaying >1000-fold change in their output after induction and show how multi-level regulation can suppress transcriptional noise and create digital-like switches between 'on' and 'off' states. These tools will aid those working with toxic genes or requiring precise regulation and propagation of cellular signals, plus illustrate the value of more diverse regulatory designs for synthetic biology.

In Vivo Feedback Control of an Antithetic Molecular-Titration Motif in Escherichia coli Using Microfluidics.

We study both in silico and in vivo the real-time feedback control of a molecular titration motif that has been earmarked as a fundamental component of antithetic and multicellular feedback control schemes in E. coli. We show that an external feedback control strategy can successfully regulate the average fluorescence output of a bacterial cell population to a desired constant level in real-time. We also provide in silico evidence that the same strategy can be used to track a time-varying reference signal where the set-point is switched to a different value halfway through the experiment. We use the experimental data to refine and parametrize an in silico model of the motif that can be used as an error computation module in future embedded or multicellular control experiments.

Micro-scale interactions between Arabidopsis root hairs and soil particles influence soil erosion.

Soil is essential for sustaining life on land. Plant roots play a crucial role in stabilising soil and minimising erosion, although these mechanisms are still not completely understood. Consequently, identifying and breeding for plant traits to enhance erosion resistance is challenging. Root hair mutants in Arabidopsis thaliana were studied using three different quantitative methods to isolate their effect on root-soil cohesion. We present compelling evidence that micro-scale interactions of root hairs with surrounding soil increase soil cohesion and reduce erosion. Arabidopsis seedlings with root hairs were more difficult to detach from soil, compost and sterile gel media than those with hairless roots, and it was 10-times harder to erode soil from roots with than without hairs. We also developed a model that can consistently predict the impact root hairs make to soil erosion resistance. Our study thus provides new insight into the mechanisms by which roots maintain soil stability.

Organization of feed-forward loop motifs reveals architectural principles in natural and engineered networks.

Network motifs are significantly overrepresented subgraphs that have been proposed as building blocks for natural and engineered networks. Detailed functional analysis has been performed for many types of motif in isolation, but less is known about how motifs work together to perform complex tasks. To address this issue, we measure the aggregation of network motifs via methods that extract precisely how these structures are connected. Applying this approach to a broad spectrum of networked systems and focusing on the widespread feed-forward loop motif, we uncover striking differences in motif organization. The types of connection are often highly constrained, differ between domains, and clearly capture architectural principles. We show how this information can be used to effectively predict functionally important nodes in the metabolic network of Escherichia coli. Our findings have implications for understanding how networked systems are constructed from motif parts and elucidate constraints that guide their evolution.

An Orthogonal Multi-input Integration System to Control Gene Expression in Escherichia coli.

In many biotechnological applications, it is useful for gene expression to be regulated by multiple signals, as this allows the programming of complex behavior. Here we implement, in Escherichia coli, a system that compares the concentration of two signal molecules, and tunes GFP expression proportionally to their relative abundance. The computation is performed via molecular titration between an orthogonal σ factor and its cognate anti-σ factor. We use mathematical modeling and experiments to show that the computation system is predictable and able to adapt GFP expression dynamically to a wide range of combinations of the two signals, and our model qualitatively captures most of these behaviors. We also demonstrate in silico the practical applicability of the system as a reference-comparator, which compares an intrinsic signal (reflecting the state of the system) with an extrinsic signal (reflecting the desired reference state) in a multicellular feedback control strategy.

BSim 2.0: An Advanced Agent-Based Cell Simulator.

Agent-based models (ABMs) provide a number of advantages relative to traditional continuum modeling approaches, permitting incorporation of great detail and realism into simulations, allowing in silico tracking of single-cell behaviors and correlation of these with emergent effects at the macroscopic level. In this study we present BSim 2.0, a radically new version of BSim, a computational ABM framework for modeling dynamics of bacteria in typical experimental environments including microfluidic chemostats. This is facilitated through the implementation of new methods including cells with capsular geometry that are able to physically and chemically interact with one another, a realistic model of cellular growth, a delay differential equation solver, and realistic environmental geometries.

A proteomic approach identifies many novel palmitoylated proteins in Arabidopsis.

S-acylation (palmitoylation) is a poorly understood post-translational modification of proteins involving the addition of acyl lipids to cysteine residues. S-acylation promotes the association of proteins with membranes and influences protein stability, microdomain partitioning, membrane targeting and activation state. No consensus motif for S-acylation exists and it therefore requires empirical identification. Here, we describe a biotin switch isobaric tagging for relative and absolute quantification (iTRAQ)-based method to identify S-acylated proteins from Arabidopsis. We use these data to predict and confirm S-acylation of proteins not in our dataset. We identified c. 600 putative S-acylated proteins affecting diverse cellular processes. These included proteins involved in pathogen perception and response, mitogen-activated protein kinases (MAPKs), leucine-rich repeat receptor-like kinases (LRR-RLKs) and RLK superfamily members, integral membrane transporters, ATPases, soluble N-ethylmaleimide-sensitive factor-activating protein receptors (SNAREs) and heterotrimeric G-proteins. The prediction of S-acylation of related proteins was demonstrated by the identification and confirmation of S-acylation sites within the SNARE and LRR-RLK families. We showed that S-acylation of the LRR-RLK FLS2 is required for a full response to elicitation by the flagellin derived peptide flg22, but is not required for localization to the plasma membrane. Arabidopsis contains many more S-acylated proteins than previously thought. These data can be used to identify S-acylation sites in related proteins. We also demonstrated that S-acylation is required for full LRR-RLK function.

BSim: an agent-based tool for modeling bacterial populations in systems and synthetic biology.

Large-scale collective behaviors such as synchronization and coordination spontaneously arise in many bacterial populations. With systems biology attempting to understand these phenomena, and synthetic biology opening up the possibility of engineering them for our own benefit, there is growing interest in how bacterial populations are best modeled. Here we introduce BSim, a highly flexible agent-based computational tool for analyzing the relationships between single-cell dynamics and population level features. BSim includes reference implementations of many bacterial traits to enable the quick development of new models partially built from existing ones. Unlike existing modeling tools, BSim fully considers spatial aspects of a model allowing for the description of intricate micro-scale structures, enabling the modeling of bacterial behavior in more realistic three-dimensional, complex environments. The new opportunities that BSim opens are illustrated through several diverse examples covering: spatial multicellular computing, modeling complex environments, population dynamics of the lac operon, and the synchronization of genetic oscillators. BSim is open source software that is freely available from http://bsim-bccs.sf.net and distributed under the Open Source Initiative (OSI) recognized MIT license. Developer documentation and a wide range of example simulations are also available from the website. BSim requires Java version 1.6 or higher.

Temperature dependence of ssrA-tag mediated protein degradation.

: Building synthetic gene networks with highly transient dynamics requires rapid protein degradation. We show that the degradation conferred by two commonly used ssrA tags is highly temperature dependent. Synthetic gene networks are being used increasingly in real-world applications where they may be subjected to variable conditions, and be required to display precise, quantitative dynamics, which will be more susceptible to environmental changes than the general qualitative dynamics focussed on so far.

Using aging to visually uncover evolutionary processes on networks.

Networks are widely used to describe many natural and technological systems. Understanding how these evolve over time poses a challenge for existing visualization techniques originally developed for fixed network structures. We describe a method of incorporating the concept of aging into evolving networks, where nodes and edges store information related to the amount of local evolutionary change they have experienced. This property is used to generate visualizations that ensure stable substructures maintain relatively fixed spatial positions, allowing them to act as visual markers and providing context for evolutionary change elsewhere. By further supplementing these visualizations with color cues, the resultant animations enable a clearer portrayal of the underlying evolutionary process.

The ankyrin repeats and DHHC S-acyl transferase domain of AKR1 act independently to regulate switching from vegetative to mating states in yeast.

Signal transduction from G-protein coupled receptors to MAPK cascades through heterotrimeric G-proteins has been described for many eukaryotic systems. One of the best-characterised examples is the yeast pheromone response pathway, which is negatively regulated by AKR1. AKR1-like proteins are present in all eukaryotes and contain a DHHC domain and six ankyrin repeats. Whilst the DHHC domain dependant S-acyl transferase (palmitoyl transferase) function of AKR1 is well documented it is not known whether the ankyrin repeats are also required for this activity. Here we show that the ankyrin repeats of AKR1 are required for full suppression of the yeast pheromone response pathway, by sequestration of the Gßγ dimer, and act independently of AKR1 S-acylation function. Importantly, the functions provided by the AKR1 ankyrin repeats and DHHC domain are not required on the same molecule to fully restore WT phenotypes and function. We also show that AKR1 molecules are S-acylated at locations other than the DHHC cysteine, increasing the abundance of AKR1 in the cell. Our results have important consequences for studies of AKR1 function, including recent attempts to characterise S-acylation enzymology and kinetics. Proteins similar to AKR1 are found in all eukaryotes and our results have broad implications for future work on these proteins and the control of switching between Gßγ regulated pathways.

A multi-functional synthetic gene network: a frequency multiplier, oscillator and switch.

We present the design and analysis of a synthetic gene network that performs frequency multiplication. It takes oscillatory transcription factor concentrations, such as those produced from the currently available genetic oscillators, as an input, and produces oscillations with half the input frequency as an output. Analysis of the bifurcation structure also reveals novel, programmable multi-functionality; in addition to functioning as a frequency multiplier, the network is able to function as a switch or an oscillator, depending on the temporal nature of the input. Multi-functionality is often observed in neuronal networks, where it is suggested to allow for the efficient coordination of different responses. This network represents a significant theoretical addition that extends the capabilities of synthetic gene networks.

Evolving enhanced topologies for the synchronization of dynamical complex networks.

Enhancing the synchronization of dynamical networks is of great interest to those designing and analyzing many man-made and natural systems. In this work, we investigate how network topology can be evolved to improve this property through the rewiring of edges. A computational tool called NETEVO performs this task using a simulated annealing metaheuristic. In contrast to other work which considers topological attributes when assessing current performance, we instead take a dynamical approach using simulated output from the system to direct the evolution of the network. Resultant topologies are analyzed using standard network measures, B matrices, and motif distributions. These uncover the convergence of many similar features for all our networks, highlighting also significant differences between those evolved using topological rather than dynamical performance measures.

A comparative analysis of synthetic genetic oscillators.

Synthetic biology is a rapidly expanding discipline at the interface between engineering and biology. Much research in this area has focused on gene regulatory networks that function as biological switches and oscillators. Here we review the state of the art in the design and construction of oscillators, comparing the features of each of the main networks published to date, the models used for in silico design and validation and, where available, relevant experimental data. Trends are apparent in the ways that network topology constrains oscillator characteristics and dynamics. Also, noise and time delay within the network can both have constructive and destructive roles in generating oscillations, and stochastic coherence is commonplace. This review can be used to inform future work to design and implement new types of synthetic oscillators or to incorporate existing oscillators into new designs.