Team:INSA-UPS France/Model
iGEM UPS-INSA Toulouse 2017
Our general strategy involves i) the detection of Vibrio cholerae using as input signal a quorum sensing molecule, ii) a molecular communication between two organisms thanks to the production of a signaling metabolite, and eventually iii) the detection of this signal to activate the production of antimicrobial peptides as output, at a lethal dose for V. cholerae. Compared to synthetic biology approaches based on a single microorganism, multi-organisms synthetic biology implies a more complex global behaviour.
Would the quorum sensing molecule (CAI-1) induce a sufficient response to activate the sensor (Vibrio harveyi)? Would the receptor be able to produce enough molecular message (diacetyl) to transmit the signal to the effector Pichia pastoris? Would the effector produce enough antimicrobial peptides to deliver the expected output, which is the death of V. cholerae to reach a non-toxic concentration? Modeling was vital to address these questions since this cannot be done intuitively.
Therefore, to understand, predict and ultimately control the behaviour of the synthetic microbial consortium, we have developed a mechanistic, dynamic model. This model has been exploited to test the feasibility of the system, characterize its emerging properties (e.g. robustness), and ultimately optimize its behaviour by adapting some of the important parameters that we can control (e.g. the initial concentrations of each microbial specie).
A model was also crucial regarding the entrepreneurship and the integrated human practices parts of our project: we needed to show to clients and investors, but also to citizens, how our system works, how we will dimensionate our device, and how long we have to wait before drinking a non-contaminated water.
Working with a complex biological system involving three microorganisms and several molecules, the Systems Biology Graphical Notation (SBGN) was a perfect way to ordinate all the elements and represent their interactions. A system of Ordinary Differential Equations (ODEs) was then written to describe and integrate the individual physical and biological processes to simulate the dynamics of our microbial consortium. Many data gathered from publications and experiments have been used to build the resulting kinetic model. This model has been used to assess the feasibility of our system and understand its properties by analyzing its sensitivity and robustness using different mathematical tools (global sensitivity analyses and an original extension of the metabolic control analysis framework). Following an iterative systems biology approach, it has been updated with experimental data. In turn, it has driven both the wet lab strategy and the device design, which is summed up in the optimization section. Finally, to discuss with non-mathematicians, a user-friendly interface has been created.
Systems biology is a new approach to understand and to manipulate biological material. It proposes to represent biological systems such as networks transforming inputs into outputs that could be modeled mathematically, similarly to technological and electronical devices.(1) It requires communication between biochemists, mathematicians and computer scientists. The SBGN (Systems Biology Graphical Notation) project, launched in 2005 to facilitate this communication, proposes a standard notation to represent biochemical and cellular processes, and to easily share biological systems to the community.(2) The standard notation can be found there.
Our system is composed of microorganisms interacting by responding to stimulations (inputs) and producing a molecular response (outputs) through a network composed of 12 entities and 21 reactions. This is a perfect illustration of a systematic use of biological material, at the scale of both a unique microorganism and a microbial consortium.
Abbreviations:
Our SBGN representation allows to present our system to different interlocutors (biologist, bioinformatician, mathematician), and to make it reusable and adaptable. This SBGN representation is also an rational way to visualize and inventory all the elements we had to consider to build the equation system.
The dynamics of our microbial consortium was summed up in 12 differential equations determining the concentrations of all the components (from genes to populations) and their dynamics. Every biological or physical process was described mathematically which were gathered to constitute our system of ODEs. Data from publications and from preliminary experiments (growth kinetics essay) were used to define the 63 parameters values. An ODEs solver from MATLAB R2017a was used to simulate the system. Under realistic conditions, simulation showed that the response time of the system (i.e. the time to reach a non-pathogenic concentration) remains below 1 hour (53.6 minutes), hence demonstrating the feasibility of the project. Different approaches used in metabolic modeling and systems biology were implemented to analyze emerging properties of the synthetic system, such as robustness and sensitivity, and understand how we can control it: Results indicated that most parameters have low control (coefficients close to 0, black), and highlighted the parameters having the stongest influence, most of them being related to the antimicrobial peptide: initial concentration of the producer Pichia pastoris, and AMP efficiency (Vc_death_rate and MIC50). Initial level of contamination of water (Vc_0) as well as the volume of the device were also identified as controlling parameter. Further simulations were thus carried out to evaluate the system efficiency for different concentrations in Pichia pastoris, under a wide range of contamination levels (concentrations in Vibrio cholerae). Simulation results, represented in the 3D graph below, indicate that the response time remains low (< 120 min) even at very high contamination levels. These two parameters, which exerts most of the control on the system, were selected as tunable parameters when creating a user-friendly simulation interface that can show their influence on the system to non-specialists. We wanted this model to be accessible to non-specialists. This required a platform to perform simulations of the behaviour of our synthetic system without the use of complex numerical computing environment such as Matlab. Thus, we developed a visual interface with Javascript to show the kinetics of Vibrio cholerae death depending of the two major parameters determined previously. Our mathematical model was useful for different parts of our project, and was also implemented by other parts:Model overview
Modeling a microbial consortium from genes to populations
Modeling approaches
SBGN: a standard representation of models in systems biology
“Systems biology is based on the understanding that the whole is greater than the sum of the parts.”
(1)
Simulating the system's dynamics
under realistic conditions (device volume, initial concentrations of each microorganism in the device, etc)→ The data, model and simulation results can be found in our simulation page
Understanding and optimizing the system
→ Sensitivity analysis methods and their results are described in our model analysis page
An intuitive modeling interface for non-specialists
→ Our visual interface is available on our interface page
Model driven optimization
References
