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.

Testing Theoretical Minimal Genomes Using Whole-Cell Models.

The minimal gene set for life has often been theorized, with at least ten produced for Mycoplasma genitalium (M. genitalium). Due to the difficulty of using M. genitalium in the lab, combined with its long replication time of 12-15 h, none of these theoretical minimal genomes have been tested, even with modern techniques. The publication of the M. genitalium whole-cell model provided the first opportunity to test them, simulating the genome edits in silico. We simulated minimal gene sets from the literature, finding that they produced in silico cells that did not divide. Using knowledge from previous research, we reintroduced specific essential and low essential genes in silico; enabling cellular division. This reinforces the need to identify species-specific low essential genes and their interactions. Any genome designs created using the currently incomplete and fragmented gene essentiality information will very likely require in vivo reintroductions to correct issues and produce dividing cells.

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.

ChipSeg: An Automatic Tool to Segment Bacterial and Mammalian Cells Cultured in Microfluidic Devices.

Extracting quantitative measurements from time-lapse images is necessary in external feedback control applications, where segmentation results are used to inform control algorithms. We describe ChipSeg, a computational tool that segments bacterial and mammalian cells cultured in microfluidic devices and imaged by time-lapse microscopy, which can be used also in the context of external feedback control. The method is based on thresholding and uses the same core functions for both cell types. It allows us to segment individual cells in high cell density microfluidic devices, to quantify fluorescent protein expression over a time-lapse experiment, and to track individual mammalian cells. ChipSeg enables robust segmentation in external feedback control experiments and can be easily customized for other experimental settings and research aims.

Understanding Metabolic Flux Behaviour in Whole-Cell Model Output.

Whole-cell modelling is a newly expanding field that has many applications in lab experiment design and predictive drug testing. Although whole-cell model output contains a wealth of information, it is complex and high dimensional and thus hard to interpret. Here, we present an analysis pipeline that combines machine learning, dimensionality reduction, and network analysis to interpret and visualise metabolic reaction fluxes from a set of single gene knockouts simulated in the Mycoplasma genitalium whole-cell model. We found that the reaction behaviours show trends that correlate with phenotypic classes of the simulation output, highlighting particular cellular subsystems that malfunction after gene knockouts. From a graphical representation of the metabolic network, we saw that there is a set of reactions that can be used as markers of a phenotypic class, showing their importance within the network. Our analysis pipeline can support the understanding of the complexity of in silico cells without detailed knowledge of the constituent parts, which can help to understand the effects of gene knockouts and, as whole-cell models become more widely built and used, aid genome design.

Minimal Genome Design Algorithms Using Whole-Cell Models.

Synthetic biologists engineer cells and cellular functions using design-build-test cycles; when the task is to extensively engineer entire genomes, the lack of appropriate design tools and biological knowledge about each gene in a cell can lengthen the process, requiring time-consuming and expensive experimental iterations.Whole-cell models represent a new avenue for genome design; the bacteria Mycoplasma genitalium has the first (and currently only published) whole-cell model which combines 28 cellular submodels and represents the integrated functions of every gene and molecule in a cell.We created two minimal genome design algorithms, GAMA and Minesweeper, that produced 1000s of in silico minimal genomes by running simulations on multiple supercomputers. Here we describe the steps to produce in silico cells with reduced genomes, combining minimisation algorithms with whole-cell model simulations.We foresee that the combination of similar algorithms and whole-cell models could later be used for a broad spectrum of genome design applications across cellular species when appropriate models become available.

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.

Computer-Aided Whole-Cell Design: Taking a Holistic Approach by Integrating Synthetic With Systems Biology.

Computer-aided design (CAD) for synthetic biology promises to accelerate the rational and robust engineering of biological systems. It requires both detailed and quantitative mathematical and experimental models of the processes to (re)design biology, and software and tools for genetic engineering and DNA assembly. Ultimately, the increased precision in the design phase will have a dramatic impact on the production of designer cells and organisms with bespoke functions and increased modularity. CAD strategies require quantitative models of cells that can capture multiscale processes and link genotypes to phenotypes. Here, we present a perspective on how whole-cell, multiscale models could transform design-build-test-learn cycles in synthetic biology. We show how these models could significantly aid in the design and learn phases while reducing experimental testing by presenting case studies spanning from genome minimization to cell-free systems. We also discuss several challenges for the realization of our vision. The possibility to describe and build whole-cells in silico offers an opportunity to develop increasingly automatized, precise and accessible CAD tools and strategies.

Author Correction: Designing minimal genomes using whole-cell models.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

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.

Designing minimal genomes using whole-cell models.

In the future, entire genomes tailored to specific functions and environments could be designed using computational tools. However, computational tools for genome design are currently scarce. Here we present algorithms that enable the use of design-simulate-test cycles for genome design, using genome minimisation as a proof-of-concept. Minimal genomes are ideal for this purpose as they have a simple functional assay whether the cell replicates or not. We used the first (and currently only published) whole-cell model for the bacterium Mycoplasma genitalium. Our computational design-simulate-test cycles discovered novel in silico minimal genomes which, if biologically correct, predict in vivo genomes smaller than JCVI-Syn3.0; a bacterium with, currently, the smallest genome that can be grown in pure culture. In the process, we identified 10 low essential genes and produced evidence for at least two Mycoplasma genitalium in silico minimal genomes. This work brings combined computational and laboratory genome engineering a step closer.

Genome-driven cell engineering review: in vivo and in silico metabolic and genome engineering.

Producing 'designer cells' with specific functions is potentially feasible in the near future. Recent developments, including whole-cell models, genome design algorithms and gene editing tools, have advanced the possibility of combining biological research and mathematical modelling to further understand and better design cellular processes. In this review, we will explore computational and experimental approaches used for metabolic and genome design. We will highlight the relevance of modelling in this process, and challenges associated with the generation of quantitative predictions about cell behaviour as a whole: although many cellular processes are well understood at the subsystem level, it has proved a hugely complex task to integrate separate components together to model and study an entire cell. We explore these developments, highlighting where computational design algorithms compensate for missing cellular information and underlining where computational models can complement and reduce lab experimentation. We will examine issues and illuminate the next steps for genome engineering.

Developing a graduate training program in Synthetic Biology: SynBioCDT.

This article presents the experience of a team of students and academics in developing a post-graduate training program in the new field of Synthetic Biology. Our Centre for Doctoral Training in Synthetic Biology (SynBioCDT) is an initiative funded by the United Kingdom's Research Councils of Engineering and Physical Sciences (EPSRC), and Biotechnology and Biological Sciences (BBSRC). SynBioCDT is a collaboration between the Universities of Oxford, Bristol and Warwick, and has been successfully running since 2014, training 78 students in this field. In this work, we discuss the organization of the taught, research and career development training. We also address the challenges faced when offering an interdisciplinary program. The article concludes with future directions to continue the development of the SynBioCDT.

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.

The Auxin-Regulated CrRLK1L Kinase ERULUS Controls Cell Wall Composition during Root Hair Tip Growth.

Root hairs facilitate a plant's ability to acquire soil anchorage and nutrients. Root hair growth is regulated by the plant hormone auxin and dependent on localized synthesis, secretion, and modification of the root hair tip cell wall. However, the exact cell wall regulators in root hairs controlled by auxin have yet to be determined. In this study, we describe the characterization of ERULUS (ERU), an auxin-induced Arabidopsis receptor-like kinase, whose expression is directly regulated by ARF7 and ARF19 transcription factors. ERU belongs to the Catharanthus roseus RECEPTOR-LIKE KINASE 1-LIKE (CrRLK1L) subfamily of putative cell wall sensor proteins. Imaging of a fluorescent fusion protein revealed that ERU is localized to the apical root hair plasma membrane. ERU regulates cell wall composition in root hairs and modulates pectin dynamics through negative control of pectin methylesterase (PME) activity. Mutant eru (-/-) root hairs accumulate de-esterified homogalacturonan and exhibit aberrant pectin Ca2+-binding site oscillations and increased PME activity. Up to 80% of the eru root hair phenotype is rescued by pharmacological supplementation with a PME-inhibiting catechin extract. ERU transcription is altered in specific cell wall-related root hair mutants, suggesting that it is a target for feedback regulation. Loss of ERU alters the phosphorylation status of FERONIA and H+-ATPases 1/2, regulators of apoplastic pH. Furthermore, H+-ATPases 1/2 and ERU are differentially phosphorylated in response to auxin. We conclude that ERULUS is a key auxin-controlled regulator of cell wall composition and pectin dynamics during root hair tip growth.