All iGEM projects involve modified organisms. When we work with those organisms, the question of confinement is essential to prevent their spreading out of the lab. Even if each team thinks about the best tool to answer this question, our team has decided to think about the worst situations.

What happens when a bacterial population escapes from our test tubes ?

To answer this question, we decided to create a computational simulation model in order to see :

• the number of bacteria that could proliferate on a lab bench
• what percentage of these bacteria would able to maintain their plasmid

The presentation of our work will be done in different sections. First, we are going to explain our approach in choosing a simulation model and our reasons for our choices. In the next section, we will describe how we have created our model and explain its initialization. After this, we will explain our mathematical choices to modelize bacterial growth and plasmid loss. Finally, we will make an assessment of our model, explain how we have validated it and give our perspectives for this project.

1. Multi-agent modelling

The definition

In informatics, a multi-agent system aims to represent intelligent agents which interact with one another and with a specific environment. In our case, our agents are bacteria, dispersed randomly on a grid (our environment). The bacteria are submitted to successive actions such as growth or division, and have specific attributes such as a an individual cell mass, or more importantly a plasmid. Our model aims to provide friendly users with information about bacterial growth. If you wish to evaluate a risk of contamination, or to experiment with parameters to characterise bacterial plasmid transfer or plasmid maintenance, you are welcome to test it out!

Why have we chosen this model?

As stated before, our modelling choice was to implement a multi-agent systems. Multi-agent modelling are frequently used to model biological phenomena. It stands out against continuous mathematical modelling, which is often used for predictions at a population level. Multi-agent modelling offers different advantages :

• Stochasticity : contrary to mathematical deterministic equations, it is much easier to induce randomness in multi-agents models.
• Heterogeneity within a population.
• BLA BLA BLA

2. Overview

Our programming choices

Our main concern in this project was the rapidity of the simulation. Indeed, bacteria can quickly develop, and our processor may have to face a horde of ever-growing agents. Hence, we chose to implement our simulation in C++.

As for the interface, we chose a much simpler implementation in Python 2.7. Our program's dependencies are the library gtkm for C++ and for Python Matplotlib.

Creating the agents and environment

Initially, an empty lattice is created, containing n number of squares. This number is determined by the height and length of the lattice provided by the user, and set by default at 1000x1000. A square is identified by its id in the lattice and can contain bacteria. If there is more than one bacterium, the bacteria are piled up in the square, with of course a limit to the piling,also set by the user.

The bacterial agents are then created, and randomly set in different squares. Each bacterial cell possesses certain characteristics such as a mass, a volume, an identifier of the square the cell is in etc.

3. How does it work?

Bacterial growth

This model can de divided into two main parts : the growth of the bacterial population and the transmission and maintenance of plasmids in the population.

For the growth of our bacterium agents, we tried two different implementations. The first one is an adaptation of bacterial growth implemented for the first time in Bacsim and based on Donachie and Robison model for of cell division. The second is adapted from the INDISIM methodology for bacterial growth simulation.

Bacterial growth based on Donachie's model

This model follows for each bacterium the following steps : INSERT KRISTINA IMAGE

Donachie's model is based on the dependence of the median cell's volume on growth rate during exponential growth. Each cell volume depends on the cell's growth rate, which itself depends on substrate uptake and metabolism.

Hence, after substrate uptake and metabolism, the growth rate of each bacterium is updated, and its new volume is calculated.

If the bacterium volume becomes superior to the volume at division, the bacterium divides in its "division neighbourhood" » but only if at least one square is unoccupied.

If the bacterium volume falls below a minimum volume, the bacterium dies and the agent is deleted from the lattice.

Plasmid loss