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.

A last-in first-out stack data structure implemented in DNA.

DNA-based memory systems are being reported with increasing frequency. However, dynamic DNA data structures able to store and recall information in an ordered way, and able to be interfaced with external nucleic acid computing circuits, have so far received little attention. Here we present an in vitro implementation of a stack data structure using DNA polymers. The stack is able to record combinations of two different DNA signals, release the signals into solution in reverse order, and then re-record. We explore the accuracy limits of the stack data structure through a stochastic rule-based model of the underlying polymerisation chemistry. We derive how the performance of the stack increases with the efficiency of washing steps between successive reaction stages, and report how stack performance depends on the history of stack operations under inefficient washing. Finally, we discuss refinements to improve molecular synchronisation and future open problems in implementing an autonomous chemical data structure.

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.

Easybiotics: a GUI for 3D physical modelling of multi-species bacterial populations.

MOTIVATION: 3D physical modelling is a powerful computational technique that allows for the simulation of complex systems such as consortia of mixed bacterial species. The complexities in physical modelling reside in the knowledge intensive model building process and the computational expense in calculating their numerical solutions. These models can offer insights into microbiology, both in understanding natural systems and as design tools for developing novel synthetic bacterial systems. Developing a robust synthetic system typically requires multiple iterations around the specify→design→build→test cycle to meet specifications. This process is laborious and expensive for both the computational and laboratory aspects, hence any improvement in any of the workflow steps would be welcomed. We have previously introduced Simbiotics, a powerful and flexible platform for designing and analyzing 3D simulations of mixed species bacterial populations. Simbiotics requires programming experience to use which creates barriers to entry for use of the tool. RESULTS: In the spirit of enabling biologists who may not have programming skills to install and utilize Simbiotics, we present in this application note Easybiotics, a user-friendly graphical user interface for Simbiotics. Users may design, simulate and analyze models from within the graphical user interface, with features such as live graph plotting and parameter sweeps. Easybiotics provides full access to all of Simbiotics simulation features, such as cell growth, motility and gene regulation. AVAILABILITY AND IMPLEMENTATION: Easybiotics and Simbiotics are free to use under the GPL3.0 licence, and can be found at: http://ico2s.org/software/simbiotics.html. We also provide readily downloadable virtual machine sandboxes to facilitate rapid installation.

Population Dynamics of Autocatalytic Sets in a Compartmentalized Spatial World.

Autocatalytic sets are self-sustaining and collectively catalytic chemical reaction networks which are believed to have played an important role in the origin of life. Simulation studies have shown that autocatalytic sets are, in principle, evolvable if multiple autocatalytic subsets can exist in different combinations within compartments, i.e., so-called protocells. However, these previous studies have so far not explicitly modeled the emergence and dynamics of autocatalytic sets in populations of compartments in a spatial environment. Here, we use a recently developed software tool to simulate exactly this scenario, as an important first step towards more realistic simulations and experiments on autocatalytic sets in protocells.

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.

Designing Uniquely Addressable Bio-orthogonal Synthetic Scaffolds for DNA and RNA Origami.

Nanotechnology and synthetic biology are rapidly converging, with DNA origami being one of the leading bridging technologies. DNA origami was shown to work well in a wide array of biotic environments. However, the large majority of extant DNA origami scaffolds utilize bacteriophages or plasmid sequences thus severely limiting its future applicability as a bio-orthogonal nanotechnology platform. In this paper we present the design of biologically inert (i.e., "bio-orthogonal") origami scaffolds. The synthetic scaffolds have the additional advantage of being uniquely addressable (unlike biologically derived ones) and hence are better optimized for high-yield folding. We demonstrate our fully synthetic scaffold design with both DNA and RNA origamis and describe a protocol to produce these bio-orthogonal and uniquely addressable origami scaffolds.

Sequence selection by dynamical symmetry breaking in an autocatalytic binary polymer model.

Template-directed replication of nucleic acids is at the essence of all living beings and a major milestone for any origin of life scenario. We present an idealized model of prebiotic sequence replication, where binary polymers act as templates for their autocatalytic replication, thereby serving as each others reactants and products in an intertwined molecular ecology. Our model demonstrates how autocatalysis alters the qualitative and quantitative system dynamics in counterintuitive ways. Most notably, numerical simulations reveal a very strong intrinsic selection mechanism that favors the appearance of a few population structures with highly ordered and repetitive sequence patterns when starting from a pool of monomers. We demonstrate both analytically and through simulation how this "selection of the dullest" is caused by continued symmetry breaking through random fluctuations in the transient dynamics that are amplified by autocatalysis and eventually propagate to the population level. The impact of these observations on related prebiotic mathematical models is discussed.

Programming chemistry in DNA-addressable bioreactors.

We present a formal calculus, termed the chemtainer calculus, able to capture the complexity of compartmentalized reaction systems such as populations of possibly nested vesicular compartments. Compartments contain molecular cargo as well as surface markers in the form of DNA single strands. These markers serve as compartment addresses and allow for their targeted transport and fusion, thereby enabling reactions of previously separated chemicals. The overall system organization allows for the set-up of programmable chemistry in microfluidic or other automated environments. We introduce a simple sequential programming language whose instructions are motivated by state-of-the-art microfluidic technology. Our approach integrates electronic control, chemical computing and material production in a unified formal framework that is able to mimic the integrated computational and constructive capabilities of the subcellular matrix. We provide a non-deterministic semantics of our programming language that enables us to analytically derive the computational and constructive power of our machinery. This semantics is used to derive the sets of all constructable chemicals and supermolecular structures that emerge from different underlying instruction sets. Because our proofs are constructive, they can be used to automatically infer control programs for the construction of target structures from a limited set of resource molecules. Finally, we present an example of our framework from the area of oligosaccharide synthesis.

The MATCHIT automaton: exploiting compartmentalization for the synthesis of branched polymers.

We propose an automaton, a theoretical framework that demonstrates how to improve the yield of the synthesis of branched chemical polymer reactions. This is achieved by separating substeps of the path of synthesis into compartments. We use chemical containers (chemtainers) to carry the substances through a sequence of fixed successive compartments. We describe the automaton in mathematical terms and show how it can be configured automatically in order to synthesize a given branched polymer target. The algorithm we present finds an optimal path of synthesis in linear time. We discuss how the automaton models compartmentalized structures found in cells, such as the endoplasmic reticulum and the Golgi apparatus, and we show how this compartmentalization can be exploited for the synthesis of branched polymers such as oligosaccharides. Lastly, we show examples of artificial branched polymers and discuss how the automaton can be configured to synthesize them with maximal yield.

Specific and reversible DNA-directed self-assembly of oil-in-water emulsion droplets.

Higher-order structures that originate from the specific and reversible DNA-directed self-assembly of microscopic building blocks hold great promise for future technologies. Here, we functionalized biotinylated soft colloid oil-in-water emulsion droplets with biotinylated single-stranded DNA oligonucleotides using streptavidin as an intermediary linker. We show the components of this modular linking system to be stable and to induce sequence-specific aggregation of binary mixtures of emulsion droplets. Three length scales were thereby involved: nanoscale DNA base pairing linking microscopic building blocks resulted in macroscopic aggregates visible to the naked eye. The aggregation process was reversible by changing the temperature and electrolyte concentration and by the addition of competing oligonucleotides. The system was reset and reused by subsequent refunctionalization of the emulsion droplets. DNA-directed self-assembly of oil-in-water emulsion droplets, therefore, offers a solid basis for programmable and recyclable soft materials that undergo structural rearrangements on demand and that range in application from information technology to medicine.

Coarse graining and scaling in dissipative particle dynamics.

Dissipative particle dynamics (DPD) is now a well-established method for simulating soft matter systems. However, its applicability was recently questioned because some investigations showed an upper coarse-graining limit that would prevent the applicability of the method to the whole mesoscopic range. This article aims to re-establish DPD as a truly mesoscopic method by analyzing the problems reported by other authors and by presenting a scaling scheme that allows one to apply DPD simulations directly to any desired length scale.

Life cycle of a minimal protocell--a dissipative particle dynamics study.

Cross-reactions and other systematic difficulties generated by the coupling of functional chemical subsystems pose the largest challenge for assembling a viable protocell in the laboratory. Our current work seeks to identify and clarify such key issues as we represent and analyze in simulation a full implementation of a minimal protocell. Using a 3D dissipative particle dynamics simulation method, we are able to address the coupled diffusion, self-assembly, and chemical reaction processes required to model a full life cycle of a protocell composed of coupled genetic, metabolic, and container subsystems. Utilizing this minimal structural and functional representation of the constituent molecules, their interactions, and their reactions, we identify and explore the nature of the many linked processes for the full protocellular system. Obviously the simplicity of this simulation method combined with the inherent system complexity prevents us from expecting quantitative simulation predictions from these investigations. However, we report important findings on systemic processes, some previously predicted and some newly discovered, as we couple the protocellular self-assembly processes and chemical reactions.

Minimal model of self-replicating nanocells: a physically embodied information-free scenario.

The building of minimal self-reproducing systems with a physical embodiment (generically called protocells) is a great challenge, with implications for both theory and applied sciences. Although the classical view of a living protocell assumes that it includes information-carrying molecules as an essential ingredient, a dividing cell-like structure can be built from a metabolism-container coupled system only. An example of such a system, modelled with dissipative particle dynamics, is presented here. This article demonstrates how a simple coupling between a precursor molecule and surfactant molecules forming micelles can experience a growth-division cycle in a predictable manner, and analyses the influence of crucial parameters on this replication cycle. Implications of these results for origins of cellular life and living technology are outlined.