DiSCUS: A Simulation Platform for Conjugation Computing

Is it selfish to be filamentous in biofilms? Individual-based modeling links microbial growth strategies with morphology using the new and modular iDynoMiCS 2.0

Microbial communities are found in all habitable environments and often occur in assemblages with self-organized spatial structures developing over time. This complexity can only be understood, predicted, and managed by combining experiments with mathematical modeling. Individual-based models are particularly suited if individual heterogeneity, local interactions, and adaptive behavior are of interest. Here we present the completely overhauled software platform, the individual-based Dynamics of Microbial Communities Simulator, iDynoMiCS 2.0, which enables researchers to specify a range of different models without having to program. Key new features and improvements are: (1) Substantially enhanced ease of use (graphical user interface, editor for model specification, unit conversions, data analysis and visualization and more). (2) Increased performance and scalability enabling simulations of up to 10 million agents in 3D biofilms. (3) Kinetics can be specified with any arithmetic function. (4) Agent properties can be assembled from orthogonal modules for pick and mix flexibility. (5) Force-based mechanical interaction framework enabling attractive forces and non-spherical agent morphologies as an alternative to the shoving algorithm. The new iDynoMiCS 2.0 has undergone intensive testing, from unit tests to a suite of increasingly complex numerical tests and the standard Benchmark 3 based on nitrifying biofilms. A second test case was based on the "biofilms promote altruism" study previously implemented in BacSim because competition outcomes are highly sensitive to the developing spatial structures due to positive feedback between cooperative individuals. We extended this case study by adding morphology to find that (i) filamentous bacteria outcompete spherical bacteria regardless of growth strategy and (ii) non-cooperating filaments outcompete cooperating filaments because filaments can escape the stronger competition between themselves. In conclusion, the new substantially improved iDynoMiCS 2.0 joins a growing number of platforms for individual-based modeling of microbial communities with specific advantages and disadvantages that we discuss, giving users a wider choice.

Pressure-dependent growth controls 3D architecture of Pseudomonas putida microcolonies

Colony formation is key to many ecological and biotechnological processes. In its early stages, colony formation involves the concourse of a number of physical and biological parameters for generation of a distinct 3D structure-the specific influence of which remains unclear. We focused on a thus far neglected aspect of the process, specifically the consequences of the differential pressure experienced by cells in the middle of a colony versus that endured by bacteria located in the growing periphery. This feature was characterized experimentally in the soil bacterium Pseudomonas putida. Using an agent-based model we recreated the growth of microcolonies in a scenario in which pressure was the only parameter affecting proliferation of cells. Simulations exposed that, due to constant collisions with other growing bacteria, cells have virtually no free space to move sideways, thereby delaying growth and boosting chances of overlapping on top of each other. This scenario was tested experimentally on agar surfaces. Comparison between experiments and simulations suggested that the inside/outside differential pressure determines growth, both timewise and in terms of spatial directions, eventually moulding colony shape. We thus argue that-at least in the case studied-mere physical pressure of growing cells suffices to explain key dynamics of colony formation.

Calibrating spatiotemporal models of microbial communities to microscopy data: A review

Spatiotemporal models that account for heterogeneity within microbial communities rely on single-cell data for calibration and validation. Such data, commonly collected via microscopy and flow cytometry, have been made more accessible by recent advances in microfluidics platforms and data processing pipelines. However, validating models against such data poses significant challenges. Validation practices vary widely between modelling studies; systematic and rigorous methods have not been widely adopted. Similar challenges are faced by the (macrobial) ecology community, in which systematic calibration approaches are often employed to improve quantitative predictions from computational models. Here, we review single-cell observation techniques that are being applied to study microbial communities and the calibration strategies that are being employed for accompanying spatiotemporal models. To facilitate future calibration efforts, we have compiled a list of summary statistics relevant for quantifying spatiotemporal patterns in microbial communities. Finally, we highlight some recently developed techniques that hold promise for improved model calibration, including algorithmic guidance of summary statistic selection and machine learning approaches for efficient model simulation.

Modelling co-translational dimerization for programmable nonlinearity in synthetic biology

Nonlinearity plays a fundamental role in the performance of both natural and synthetic biological networks. Key functional motifs in living microbial systems, such as the emergence of bistability or oscillations, rely on nonlinear molecular dynamics. Despite its core importance, the rational design of nonlinearity remains an unmet challenge. This is largely due to a lack of mathematical modelling that accounts for the mechanistic basis of nonlinearity. We introduce a model for gene regulatory circuits that explicitly simulates protein dimerization-a well-known source of nonlinear dynamics. Specifically, our approach focuses on modelling co-translational dimerization: the formation of protein dimers during-and not after-translation. This is in contrast to the prevailing assumption that dimer generation is only viable between freely diffusing monomers (i.e. post-translational dimerization). We provide a method for fine-tuning nonlinearity on demand by balancing the impact of co- versus post-translational dimerization. Furthermore, we suggest design rules, such as protein length or physical separation between genes, that may be used to adjust dimerization dynamics in vivo. The design, build and test of genetic circuits with on-demand nonlinear dynamics will greatly improve the programmability of synthetic biological systems.

From Microbial Communities to Distributed Computing Systems

A distributed biological system can be defined as a system whose components are located in different subpopulations, which communicate and coordinate their actions through interpopulation messages and interactions. We see that distributed systems are pervasive in nature, performing computation across all scales, from microbial communities to a flock of birds. We often observe that information processing within communities exhibits a complexity far greater than any single organism. Synthetic biology is an area of research which aims to design and build synthetic biological machines from biological parts to perform a defined function, in a manner similar to the engineering disciplines. However, the field has reached a bottleneck in the complexity of the genetic networks that we can implement using monocultures, facing constraints from metabolic burden and genetic interference. This makes building distributed biological systems an attractive prospect for synthetic biology that would alleviate these constraints and allow us to expand the applications of our systems into areas including complex biosensing and diagnostic tools, bioprocess control and the monitoring of industrial processes. In this review we will discuss the fundamental limitations we face when engineering functionality with a monoculture, and the key areas where distributed systems can provide an advantage. We cite evidence from natural systems that support arguments in favor of distributed systems to overcome the limitations of monocultures. Following this we conduct a comprehensive overview of the synthetic communities that have been built to date, and the components that have been used. The potential computational capabilities of communities are discussed, along with some of the applications that these will be useful for. We discuss some of the challenges with building co-cultures, including the problem of competitive exclusion and maintenance of desired community composition. Finally, we assess computational frameworks currently available to aide in the design of microbial communities and identify areas where we lack the necessary tools.

Toward Engineering Biosystems With Emergent Collective Functions

Many complex behaviors in biological systems emerge from large populations of interacting molecules or cells, generating functions that go beyond the capabilities of the individual parts. Such collective phenomena are of great interest to bioengineers due to their robustness and scalability. However, engineering emergent collective functions is difficult because they arise as a consequence of complex multi-level feedback, which often spans many length-scales. Here, we present a perspective on how some of these challenges could be overcome by using multi-agent modeling as a design framework within synthetic biology. Using case studies covering the construction of synthetic ecologies to biological computation and synthetic cellularity, we show how multi-agent modeling can capture the core features of complex multi-scale systems and provide novel insights into the underlying mechanisms which guide emergent functionalities across scales. The ability to unravel design rules underpinning these behaviors offers a means to take synthetic biology beyond single molecules or cells and toward the creation of systems with functions that can only emerge from collectives at multiple scales.

Pathways to cellular supremacy in biocomputing

Synthetic biology uses living cells as the substrate for performing human-defined computations. Many current implementations of cellular computing are based on the “genetic circuit” metaphor, an approximation of the operation of silicon-based computers. Although this conceptual mapping has been relatively successful, we argue that it fundamentally limits the types of computation that may be engineered inside the cell, and fails to exploit the rich and diverse functionality available in natural living systems. We propose the notion of “cellular supremacy” to focus attention on domains in which biocomputing might offer superior performance over traditional computers. We consider potential pathways toward cellular supremacy, and suggest application areas in which it may be found.

Synthetic biology uses cells as its computing substrate, often based on the genetic circuit concept. In this Perspective, the authors argue that existing synthetic biology approaches based on classical models of computation limit the potential of biocomputing, and propose that living organisms have under-exploited capabilities.

A Model for the Spatiotemporal Design of Gene Regulatory Circuits †

Mathematical modeling assists the design of synthetic regulatory networks by providing a detailed mechanistic understanding of biological systems. Models that can predict the performance of a design are fundamental for synthetic biology since they minimize iterations along the design-build-test lifecycle. Such predictability depends crucially on what assumptions (i.e., biological simplifications) the model considers. Here, we challenge a common assumption when it comes to the modeling of bacterial-based gene regulation: considering negligible the effects of intracellular physical space. It is commonly assumed that molecules, such as transcription factors (TF), are homogeneously distributed inside a cell, so there is no need to model their diffusion. We describe a mathematical model that accounts for molecular diffusion and show how simulations of network performance are decisively affected by the distance between its components. Specifically, the model focuses on the search by a TF for its target promoter. The combination of local searches, via one-dimensional sliding along the chromosome, and global searches, via three-dimensional diffusion through the cytoplasm, determine TF-promoter interplay. Previous experimental results with engineered bacteria in which the distance between TF source and target was minimized or enlarged were successfully reproduced by the spatially resolved model we introduce here. This suggests that the spatial specification of the circuit alone can be exploited as a design parameter in synthetic biology to select programmable output levels.

High-Performance Biocomputing in Synthetic Biology-Integrated Transcriptional and Metabolic Circuits

Biocomputing uses molecular biology parts as the hardware to implement computational devices. By following pre-defined rules, often hard-coded into biological systems, these devices are able to process inputs and return outputs—thus computing information. Key to the success of any biocomputing endeavor is the availability of a wealth of molecular tools and biological motifs from which functional devices can be assembled. Synthetic biology is a fabulous playground for such purpose, offering numerous genetic parts that allow for the rational engineering of genetic circuits that mimic the behavior of electronic functions, such as logic gates. A grand challenge, as far as biocomputing is concerned, is to expand the molecular hardware available beyond the realm of genetic parts by tapping into the host metabolism. This objective requires the formalization of the interplay of genetic constructs with the rest of the cellular machinery. Furthermore, the field of metabolic engineering has had little intersection with biocomputing thus far, which has led to a lack of definition of metabolic dynamics as computing basics. In this perspective article, we advocate the conceptualization of metabolism and its motifs as the way forward to achieve whole-cell biocomputations. The design of merged transcriptional and metabolic circuits will not only increase the amount and type of information being processed by a synthetic construct, but will also provide fundamental control mechanisms for increased reliability.

Deconvolution of Gene Expression Noise into Spatial Dynamics of Transcription Factor-Promoter Interplay

Gene expression noise is not only the mere consequence of stochasticity, but also a signal that reflects the upstream physical dynamics of the cognate molecular machinery. Soil bacteria facing recalcitrant pollutants exploit noise of catabolic promoters to deploy beneficial phenotypes such as metabolic bet-hedging and/or division of biochemical labor. Although the role of upstream promoter-regulator interplay in the origin of this noise is little understood, its specifications are probably ciphered in flow cytometry data patterns. We studied Pm promoter activity of the environmental bacterium Pseudomonas putida and its cognate regulator XylS by following expression of Pm-gfp fusions in single cells. Using mathematical modeling and computational simulations, we determined the kinetic properties of the system and used them as a baseline code to interpret promoter activity in terms of upstream regulator dynamics. Transcriptional noise was predicted to depend on the intracellular physical distance between regulator source (where XylS is produced) and the target promoter. Experiments with engineered bacteria in which this distance is minimized or enlarged confirmed the predicted effects of source/target proximity on noise patterns. This approach allowed deconvolution of cytometry data into mechanistic information on gene expression flow. It also provided a basis for selecting programmable noise levels in synthetic regulatory circuits.

Simbiotics: A Multiscale Integrative Platform for 3D Modeling of Bacterial Populations

Simbiotics is a spatially explicit multiscale modeling platform for the design, simulation and analysis of bacterial populations. Systems ranging from planktonic cells and colonies, to biofilm formation and development may be modeled. Representation of biological systems in Simbiotics is flexible, and user-defined processes may be in a variety of forms depending on desired model abstraction. Simbiotics provides a library of modules such as cell geometries, physical force dynamics, genetic circuits, metabolic pathways, chemical diffusion and cell interactions. Model defined processes are integrated and scheduled for parallel multithread and multi-CPU execution. A virtual lab provides the modeler with analysis modules and some simulated lab equipment, enabling automation of sample interaction and data collection. An extendable and modular framework allows for the platform to be updated as novel models of bacteria are developed, coupled with an intuitive user interface to allow for model definitions with minimal programming experience. Simbiotics can integrate existing standards such as SBML, and process microscopy images to initialize the 3D spatial configuration of bacteria consortia. Two case studies, used to illustrate the platform flexibility, focus on the physical properties of the biosystems modeled. These pilot case studies demonstrate Simbiotics versatility in modeling and analysis of natural systems and as a CAD tool for synthetic biology.

Predictive Modeling of a Batch Filter Mating Process

Quantitative characterizations of horizontal gene transfer are needed to accurately describe gene transfer processes in natural and engineered systems. A number of approaches to the quantitative description of plasmid conjugation have appeared in the literature. In this study, we seek to extend that work, motivated by the question of whether a mathematical model can accurately predict growth and conjugation dynamics in a batch process. We used flow cytometry to make time-point observations of a filter-associated mating between two E. coli strains, and fit ordinary differential equation models to the data. A model comparison analysis identified the model formulation that is best supported by the data. Identifiability analysis revealed that the parameters were estimated with acceptable accuracy. The predictive power of the model was assessed by comparison with test data that demanded extrapolation from the training experiments. This study represents the first attempt to assess the quality of model predictions for plasmid conjugation. Our successful application of this approach lays a foundation for predictive modeling that can be used both in the study of natural plasmid transmission and in model-based design of engineering approaches that employ conjugation, such as plasmid-mediated bioaugmentation.