A New Improved and Extended Version of the Multicell Bacterial Simulator gro

Biology System Description Language (BiSDL): a modeling language for the design of multicellular synthetic biological systems

BACKGROUND

The Biology System Description Language (BiSDL) is an accessible, easy-to-use computational language for multicellular synthetic biology. It allows synthetic biologists to represent spatiality and multi-level cellular dynamics inherent to multicellular designs, filling a gap in the state of the art. Developed for designing and simulating spatial, multicellular synthetic biological systems, BiSDL integrates high-level conceptual design with detailed low-level modeling, fostering collaboration in the Design-Build-Test-Learn cycle. BiSDL descriptions directly compile into Nets-Within-Nets (NWNs) models, offering a unique approach to spatial and hierarchical modeling in biological systems.

RESULTS

BiSDL's effectiveness is showcased through three case studies on complex multicellular systems: a bacterial consortium, a synthetic morphogen system and a conjugative plasmid transfer process. These studies highlight the BiSDL proficiency in representing spatial interactions and multi-level cellular dynamics. The language facilitates the compilation of conceptual designs into detailed, simulatable models, leveraging the NWNs formalism. This enables intuitive modeling of complex biological systems, making advanced computational tools more accessible to a broader range of researchers.

CONCLUSIONS

BiSDL represents a significant step forward in computational languages for synthetic biology, providing a sophisticated yet user-friendly tool for designing and simulating complex biological systems with an emphasis on spatiality and cellular dynamics. Its introduction has the potential to transform research and development in synthetic biology, allowing for deeper insights and novel applications in understanding and manipulating multicellular systems.

Synthetic biological neural networks: From current implementations to future perspectives

Artificial neural networks, inspired by the biological networks of the human brain, have become game-changing computing models in modern computer science. Inspired by their wide scope of applications, synthetic biology strives to create their biological counterparts, which we denote synthetic biological neural networks (SYNBIONNs). Their use in the fields of medicine, biosensors, biotechnology, and many more shows great potential and presents exciting possibilities. So far, many different synthetic biological networks have been successfully constructed, however, SYNBIONN implementations have been sparse. The latter are mostly based on neural networks pretrained in silico and being heavily dependent on extensive human input. In this paper, we review current implementations and models of SYNBIONNs. We briefly present the biological platforms that show potential for designing and constructing perceptrons and/or multilayer SYNBIONNs. We explore their future possibilities along with the challenges that must be overcome to successfully implement a scalable in vivo biological neural network capable of online learning.

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.

pyBrick-DNA: A Python-Based Environment for Automated Genetic Component Assembly

Genetic component assembly is key in the simulation and implementation of genetic circuits. Automating this process, thus accelerating prototyping, is a necessity. We present pyBrick-DNA, a software written in Python, that assembles components for the construction of genetic circuits. pyBrick-DNA (colab.pyBrick.com) is a user-friendly environment where scientists can select genetic sequences or input custom sequences to build genetic assemblies. All components are modularly fused to generate a ready-to-go single DNA fragment. It includes Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR) and plant gene-editing components. Hence, pyBrick-DNA can generate a functional CRISPR construct composed of a single-guided RNA integrated with Cas9, promoters, and terminator elements. The outcome is a DNA sequence, along with a graphical representation, composed of user-selected genetic parts, ready to be synthesized and cloned in vivo. Moreover, the sequence can be exported as a GenBank file allowing its use with other synthetic biology tools.

A roadmap for the creation of synthetic lichen

Lichens represent a charismatic corner of biology that has a rich history of scientific exploration, but to which modern biological techniques have been sparsely applied. This has limited our understanding of phenomena unique to lichen, such as the emergent development of physically coupled microbial consortia or distributed metabolisms. The experimental intractability of natural lichens has prevented studies of the mechanistic underpinnings of their biology. Creating synthetic lichen from experimentally tractable, free-living microbes has the potential to overcome these challenges. They could also serve as powerful new chassis for sustainable biotechnology. In this review we will first briefly introduce what lichen are, what remains mysterious about their biology, and why. We will then articulate the scientific insights that creating a synthetic lichen will generate and lay out a roadmap for how this could be achieved using synthetic biology. Finally, we will explore the translational applications of synthetic lichen and detail what is needed to advance the pursuit of their creation.

Novel Ground-Up 3D Multicellular Simulators for Synthetic Biology CAD Integrating Stochastic Gillespie Simulations Benchmarked with Topologically Variable SBML Models

The elevation of Synthetic Biology from single cells to multicellular simulations would be a significant scale-up. The spatiotemporal behavior of cellular populations has the potential to be prototyped in silico for computer assisted design through ergonomic interfaces. Such a platform would have great practical potential across medicine, industry, research, education and accessible archiving in bioinformatics. Existing Synthetic Biology CAD systems are considered limited regarding population level behavior, and this work explored the in silico challenges posed from biological and computational perspectives. Retaining the connection to Synthetic Biology CAD, an extension of the Infobiotics Workbench Suite was considered, with potential for the integration of genetic regulatory models and/or chemical reaction networks through Next Generation Stochastic Simulator (NGSS) Gillespie algorithms. These were executed using SBML models generated by in-house SBML-Constructor over numerous topologies and benchmarked in association with multicellular simulation layers. Regarding multicellularity, two ground-up multicellular solutions were developed, including the use of Unreal Engine 4 contrasted with CPU multithreading and Blender visualization, resulting in a comparison of real-time versus batch-processed simulations. In conclusion, high-performance computing and client-server architectures could be considered for future works, along with the inclusion of numerous biologically and physically informed features, whilst still pursuing ergonomic solutions.

Agent-Based Modeling of Microbial Communities

Microbial communities are complex living systems that populate the planet with diverse functions and are increasingly harnessed for practical human needs. To deepen the fundamental understanding of their organization and functioning as well as to facilitate their engineering for applications, mathematical modeling has played an increasingly important role. Agent-based models represent a class of powerful quantitative frameworks for investigating microbial communities because of their individualistic nature in describing cells, mechanistic characterization of molecular and cellular processes, and intrinsic ability to produce emergent system properties. This review presents a comprehensive overview of recent advances in agent-based modeling of microbial communities. It surveys the state-of-the-art algorithms employed to simulate intracellular biomolecular events, single-cell behaviors, intercellular interactions, and interactions between cells and their environments that collectively serve as the driving forces of community behaviors. It also highlights three lines of applications of agent-based modeling, namely, the elucidation of microbial range expansion and colony ecology, the design of synthetic gene circuits and microbial populations for desired behaviors, and the characterization of biofilm formation and dispersal. The review concludes with a discussion of existing challenges, including the computational cost of the modeling, and potential mitigation strategies.

CellSium: versatile cell simulator for microcolony ground truth generation

Summary

To train deep learning-based segmentation models, large ground truth datasets are needed. To address this need in microfluidic live-cell imaging, we present CellSium, a flexibly configurable cell simulator built to synthesize realistic image sequences of bacterial microcolonies growing in monolayers. We illustrate that the simulated images are suitable for training neural networks. Synthetic time-lapse videos with and without fluorescence, using programmable cell growth models, and simulation-ready 3D colony geometries for computational fluid dynamics are also supported.

Availability and implementation

CellSium is free and open source software under the BSD license, implemented in Python, available at github.com/modsim/cellsium (DOI: 10.5281/zenodo.6193033), along with documentation, usage examples and Docker images.

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

Ecological modelling approaches for predicting emergent properties in microbial communities

Recent studies have brought forward the critical role of emergent properties in shaping microbial communities and the ecosystems of which they are a part. Emergent properties-patterns or functions that cannot be deduced linearly from the properties of the constituent parts-underlie important ecological characteristics such as resilience, niche expansion and spatial self-organization. While it is clear that emergent properties are a consequence of interactions within the community, their non-linear nature makes mathematical modelling imperative for establishing the quantitative link between community structure and function. As the need for conservation and rational modulation of microbial ecosystems is increasingly apparent, so is the consideration of the benefits and limitations of the approaches to model emergent properties. Here we review ecosystem modelling approaches from the viewpoint of emergent properties. We consider the scope, advantages and limitations of Lotka-Volterra, consumer-resource, trait-based, individual-based and genome-scale metabolic models. Future efforts in this research area would benefit from capitalizing on the complementarity between these approaches towards enabling rational modulation of complex microbial ecosystems.

Resource-aware whole-cell model of division of labour in a microbial consortium for complex-substrate degradation

BACKGROUND

Low-cost sustainable feedstocks are essential for commercially viable biotechnologies. These feedstocks, often derived from plant or food waste, contain a multitude of different complex biomolecules which require multiple enzymes to hydrolyse and metabolise. Current standard biotechnology uses monocultures in which a single host expresses all the proteins required for the consolidated bioprocess. However, these hosts have limited capacity for expressing proteins before growth is impacted. This limitation may be overcome by utilising division of labour (DOL) in a consortium, where each member expresses a single protein of a longer degradation pathway.

RESULTS

Here, we model a two-strain consortium, with one strain expressing an endohydrolase and a second strain expressing an exohydrolase, for cooperative degradation of a complex substrate. Our results suggest that there is a balance between increasing expression to enhance degradation versus the burden that higher expression causes. Once a threshold of burden is reached, the consortium will consistently perform better than an equivalent single-cell monoculture.

CONCLUSIONS

We demonstrate that resource-aware whole-cell models can be used to predict the benefits and limitations of using consortia systems to overcome burden. Our model predicts the region of expression where DOL would be beneficial for growth on starch, which will assist in making informed design choices for this, and other, complex-substrate degradation pathways.

Substrate availability and toxicity shape the structure of microbial communities engaged in metabolic division of labor

Metabolic division of labor (MDOL) represents a widespread natural phenomenon, whereby a complex metabolic pathway is shared between different strains within a community in a mutually beneficial manner. However, little is known about how the composition of such a microbial community is regulated. We hypothesized that when degradation of an organic compound is carried out via MDOL, the concentration and toxicity of the substrate modulate the benefit allocation between the two microbial populations, thus affecting the structure of this community. We tested this hypothesis by combining modeling with experiments using a synthetic consortium. Our modeling analysis suggests that the proportion of the population executing the first metabolic step can be simply estimated by Monod-like formulas governed by substrate concentration and toxicity. Our model and the proposed formula were able to quantitatively predict the structure of our synthetic consortium. Further analysis demonstrates that our rule is also applicable in estimating community structures in spatially structured environments. Together, our work clearly demonstrates that the structure of MDOL communities can be quantitatively predicted using available information on environmental factors, thus providing novel insights into how to manage artificial microbial systems for the wide application of the bioindustry.

A metabolic modeling platform for the computation of microbial ecosystems in time and space (COMETS)

Genome-scale stoichiometric modeling of metabolism has become a standard systems biology tool for modeling cellular physiology and growth. Extensions of this approach are emerging as a valuable avenue for predicting, understanding and designing microbial communities. Computation of microbial ecosystems in time and space (COMETS) extends dynamic flux balance analysis to generate simulations of multiple microbial species in molecularly complex and spatially structured environments. Here we describe how to best use and apply the most recent version of COMETS, which incorporates a more accurate biophysical model of microbial biomass expansion upon growth, evolutionary dynamics and extracellular enzyme activity modules. In addition to a command-line option, COMETS includes user-friendly Python and MATLAB interfaces compatible with the well-established COBRA models and methods, as well as comprehensive documentation and tutorials. This protocol provides a detailed guideline for installing, testing and applying COMETS to different scenarios, generating simulations that take from a few minutes to several days to run, with broad applicability to microbial communities across biomes and scales.

Toward Full-Stack In Silico Synthetic Biology: Integrating Model Specification, Simulation, Verification, and Biological Compilation

We present the Infobiotics Workbench (IBW), a user-friendly, scalable, and integrated computational environment for the computer-aided design of synthetic biological systems. It supports an iterative workflow that begins with specification of the desired synthetic system, followed by simulation and verification of the system in high-performance environments and ending with the eventual compilation of the system specification into suitable genetic constructs. IBW integrates modeling, simulation, verification, and biocompilation features into a single software suite. This integration is achieved through a new domain-specific biological programming language, the Infobiotics Language (IBL), which tightly combines these different aspects of in silico synthetic biology into a full-stack integrated development environment. Unlike existing synthetic biology modeling or specification languages, IBL uniquely blends modeling, verification, and biocompilation statements into a single file. This allows biologists to incorporate design constraints within the specification file rather than using decoupled and independent formalisms for different in silico analyses. This novel approach offers seamless interoperability across different tools as well as compatibility with SBOL and SBML frameworks and removes the burden of doing manual translations for standalone applications. We demonstrate the features, usability, and effectiveness of IBW and IBL using well-established synthetic biological circuits.

Control engineering and synthetic biology: working in synergy for the analysis and control of microbial systems

The implementation of novel functionalities in living cells is a key aspect of synthetic biology. In the last decade, the field of synthetic biology has made progress working in synergy with control engineering, whose solid framework has provided concepts and tools to analyse biological systems and guide their design. In this review, we briefly highlight recent work focused on the application of control theoretical concepts and tools for the analysis and design of synthetic biology systems in microbial cells.

A Framework for Implementing Metaheuristic Algorithms Using Intercellular Communication

Metaheuristics (MH) are Artificial Intelligence procedures that frequently rely on evolution. MH approximate difficult problem solutions, but are computationally costly as they explore large solution spaces. This work pursues to lay the foundations of general mappings for implementing MH using Synthetic Biology constructs in cell colonies. Two advantages of this approach are: harnessing large scale parallelism capability of cell colonies and, using existing cell processes to implement basic dynamics defined in computational versions. We propose a framework that maps MH elements to synthetic circuits in growing cell colonies to replicate MH behavior in cell colonies. Cell-cell communication mechanisms such as quorum sensing (QS), bacterial conjugation, and environmental signals map to evolution operators in MH techniques to adapt to growing colonies. As a proof-of-concept, we implemented the workflow associated to the framework: automated MH simulation generators for the gro simulator and two classes of algorithms (Simple Genetic Algorithms and Simulated Annealing) encoded as synthetic circuits. Implementation tests show that synthetic counterparts mimicking MH are automatically produced, but also that cell colony parallelism speeds up the execution in terms of generations. Furthermore, we show an example of how our framework is extended by implementing a different computational model: The Cellular Automaton.

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.

A Microbial Screening in Silico Method for the Fitness Step Evaluation in Evolutionary Algorithms

One of the most delicate stages of an evolutionary algorithm is the evaluation of the goodness of the solutions by some procedure providing a fitness value. However, although there are general rules, it is not always easy to find an appropriate evaluation function for a given problem. In the biological realm, today, there is a variety of experimental methods under the name of microbial screening to identify and select bacteria from their traits, as well as to obtain their fitness. In this paper, we show how given an optimization problem, a colony of synthetic bacteria or bacterial agents is able to evaluate the fitness of candidate solutions by building an evaluation function. The evaluation function is obtained simulating, in silico, a bacterial colony conducting the laboratory methods used in microbiology, biotechnology and synthetic biology to measure microbial fitness. Once the evaluation function is built, it is included in the code of the genetic algorithm as part of the fitness routine. The practical use of this approach is illustrated in two classic optimization problems. In silico routines have been programmed in Gro, a cell programming language oriented to synthetic biology, and can easily be customized to many other optimization problems.

Hybrid Model of Bacterial Biofilm Growth

Bacterial biofilms play a critical role in environmental processes, water treatment, human health, and food processing. They exhibit highly complex dynamics due to the interactions between the bacteria and the extracellular polymeric substance (EPS), water, and nutrients and minerals that make up the biofilm. We present a hybrid computational model in which the dynamics of discrete bacterial cells are simulated within a multiphase continuum, consisting of EPS and water as separate interacting phases, through which nutrients and minerals diffuse. Bacterial cells in our model consume water and nutrients in order to grow, divide, and produce EPS. Consequently, EPS flows outward from the bacterial colony, while water flows inward. The model predicts bacterial colony formation as a treelike structure. The distribution of bacterial growth and EPS production is found to be sensitive to the pore spacing between bacteria and the consumption of nutrients within the bacterial colony. Forces that are sometimes neglected in biofilm simulations, such as lubrication force between nearby bacterial cells and osmotic (swelling) pressure force resulting from gradients in EPS concentration, are observed to have an important effect on biofilm growth via their influence on bacteria pore spacing and associated water/nutrient percolation into the bacterial colony.

Spatial propagation of electrical signals in circular biofilms: A combined experimental and agent-based fire-diffuse-fire study

Bacterial biofilms are a risk to human health, playing critical roles in persistent infections. Recent studies have observed electrical signaling in biofilms and thus biofilms represent a new class of active excitable matter in which cell division is the active process and the spiking of the individual bacterial cells is the excitable process. Electrophysiological models have predominantly been developed to describe eukaryotic systems, but we demonstrate their use in understanding bacterial biofilms. Our agent-based fire-diffuse-fire (ABFDF) model successfully simulates the propagation of both centrifugal (away from the center) and centripetal (toward the center) electrical signals through biofilms of Bacillus subtilis. Furthermore, the ABFDF model allows realistic spatial positioning of the bacteria in two dimensions to be included in the fire-diffuse-fire model and this is the crucial factor that improves agreement with experiments. The speed of propagation is not constant and depends on the radius of the propagating electrical wave front. Centripetal waves are observed to move faster than centrifugal waves, which is a curvature driven effect and is correctly captured by our simulations.

The ecology and evolution of microbial metabolic strategies

Free-living microbes are generally capable of growing on multiple different nutrients. Some of those nutrients are used simultaneously, while others are used sequentially. The pattern of nutrient preferences and co-utilization defines the metabolic strategy of a microorganism. Metabolic strategies can substantially affect ecological interactions between species, but their evolution and distribution across the tree of life remain poorly characterized. We discuss how the confluence of better computational models of genotype-phenotype maps and high-throughput experimental tools can help us fill gaps in our knowledge and incorporate metabolic strategies into quantitative predictive models of microbial consortia.

Synthetic Biology Tools to Engineer Microbial Communities for Biotechnology

Microbial consortia have been used in biotechnology processes, including fermentation, waste treatment, and agriculture, for millennia. Today, synthetic biologists are increasingly engineering microbial consortia for diverse applications, including the bioproduction of medicines, biofuels, and biomaterials from inexpensive carbon sources. An improved understanding of natural microbial ecosystems, and the development of new tools to construct synthetic consortia and program their behaviors, will vastly expand the functions that can be performed by communities of interacting microorganisms. Here, we review recent advancements in synthetic biology tools and approaches to engineer synthetic microbial consortia, discuss ongoing and emerging efforts to apply consortia for various biotechnological applications, and suggest future applications.

Small Universal Bacteria and Plasmid Computing Systems

Bacterial computing is a known candidate in natural computing, the aim being to construct "bacterial computers" for solving complex problems. In this paper, a new kind of bacterial computing system, named the bacteria and plasmid computing system (BP system), is proposed. We investigate the computational power of BP systems with finite numbers of bacteria and plasmids. Specifically, it is obtained in a constructive way that a BP system with 2 bacteria and 34 plasmids is Turing universal. The results provide a theoretical cornerstone to construct powerful bacterial computers and demonstrate a concept of paradigms using a "reasonable" number of bacteria and plasmids for such devices.

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.

Modeling mechanical interactions in growing populations of rod-shaped bacteria

Advances in synthetic biology allow us to engineer bacterial collectives with pre-specified characteristics. However, the behavior of these collectives is difficult to understand, as cellular growth and division as well as extra-cellular fluid flow lead to complex, changing arrangements of cells within the population. To rationally engineer and control the behavior of cell collectives we need theoretical and computational tools to understand their emergent spatiotemporal dynamics. Here, we present an agent-based model that allows growing cells to detect and respond to mechanical interactions. Crucially, our model couples the dynamics of cell growth to the cell's environment: Mechanical constraints can affect cellular growth rate and a cell may alter its behavior in response to these constraints. This coupling links the mechanical forces that influence cell growth and emergent behaviors in cell assemblies. We illustrate our approach by showing how mechanical interactions can impact the dynamics of bacterial collectives growing in microfluidic traps.