Overview
In synthetic biology, Ordinary Differential Equation (ODE)-based models are essential for the precise modeling of gene regulatory circuits. These models allow for the mathematical description of dynamic interactions between genes, proteins, and molecules, facilitating insights into how gene expression changes over time1. Using ODE-based models, precise control over gene expression is achieved, ensuring the engineering of tailored biological functions in synthetic biology2.
For our project’s modeling, we have developed a comprehensive mathematical model tailored to the biological components of our project. The primary objective of our model is to enhance the design process in the wet lab by providing a structured, in silico representation that emphasizes modularity.
The System
The toxicity of oil mill wastewater (OMW) makes it a particularly destructive by-product capable of causing serious environmental complications. This is primarily associated with the high concentrations of phenolic compounds in OMW. Given these circumstances, the iGEM Thessaly wet lab team undertook the development of a system comprising two bacteria that collaborate in effectively degrading these phenolic compounds.
P. putida possesses the necessary mechanisms for the degradation of synthetic phenolic compounds, concurrently yielding a valuable by-product known as Polyhydroxyalkanoates (PHAs).
The use of E. coli is equally decisive, as it has an assistive role in the degradation of phenolic compounds. This is achieved through a mechanism based on the regulation of laccase production.
The mechanisms are based on a system of controlled protein expression by different regulatory factors. Utilizing activation and inhibition, the expression of each gene in the system can be regulated and controlled to produce the required output value.
The integration of P. Putida and E. Coli inside a microbial consortium provides a promising answer for minimizing the environmental effect of OMW and additionally provides a very valuable product (PHAs). For further information regarding this bacterial synergy, visit our Design page.
Our Aim
Using differential equations as the fundamental mathematical approach, the focus of the modeling is to define the dynamics of the two-bacteria system with sufficient accuracy to obtain a better estimation of protein and enzyme production rates. Some molecules with high monitoring value are PcaU, FadR and LacI as they can regulate the expression of other genes as inhibitors or activators. The simulation should provide useful information for the clarification of the dependencies regulating the kinetics of molecule productions. It should also help determine the best possible promoter to achieve a faster expression of the gene. The main objective is to develop a model that approximates reality to a limited but acceptable degree. Our modeling work serves the following key purposes:
Assumptions
The model developed does not fully correspond to reality as it is exceedingly complicated to represent multiple complex processes with a simple deterministic approach. In consideration of the project's optimal advancement, several assumptions have been defined.
Model Overview
Employing the MATLAB tool, SimBiology, we formulated a computational model of Pseudomonas putida composed of three distinct subsystems. These subsystems are:
Parameters & Equations
To analyze the mechanisms of our model, it is essential to specify and define the associated parameters. Careful selection of parameters is critical, as it is the foundation on which further evaluations are based. This involves a thorough examination of scientific papers, MATLAB Projects and model architectures.
In developing our model, we used a combination of ordinary differential equations (ODEs) to describe key processes such as translation, binding and transcription. This deterministic approach allows us to understand the complex interaction of molecules, enabling a better visualization of the regulation of gene expressions. While most of the equations are relatively simple, the repression of a gene’s expression from another transcriptional factor (enzyme) is considered and approached with the necessary equations.
Model Simulation
The plots generated by the model simulation provide a visual representation of the concentrations of each species in the subsystems. This visualization provides valuable insights into the dynamic behavior of the biochemical network. By communicating effectively with the wet lab team, we successfully conducted experiments to modify various parameters in the model, to meet specific requirements of the laboratory research. Main objective of the simulation was to uncover the regulatory mechanisms underlying the expression of the enzymes and the prior processes that controlled their concentration levels. The plots in figure 3 showed the concentration levels of the target enzymes that are mainly Phaj, PHAC1, PHAC2, AlkK and PhaG. We adjusted the parameters that controlled their expression such as the translation rates, binding and unbinding rates and the transcription rates. On example is the adjustments on the transcription rate of phaj mRNA where the repression equation was used.
Model Overview
Employing the MATLAB tool, SimBiology, we formulated a computational model of Escherichia Coli composed of three distinct subsystems. These subsystems are:
- The araC system. This system involves regulating the production of the enzyme laccase by utilizing the enzyme AraC while using Arabinose. AraC reacts with arabinose, and they form a complex which acts as an activator for the transcription of the lac gene.
- The TE system. This is a rather simple subsystem that can provide useful information about the transcription of the TE gene and its translation to produce Thioesterase.
- The FapR system. This system serves to regulate the expression of the enzyme carboxylase of acetyl-CoA (ACC), which is crucial for the catalytic conversion of acetyl-CoA to malonyl-CoA. This regulatory mechanism is accomplished through the regulated intervention of transcription factors, LacI and FapR, which act as repressors. Basically, the system incorporates a negative feedback loop, with FapR and LacI having inhibitory effects. The whole process depends on the presence of two critical molecules, Malonyl-CoA and isopropyl β-D-1-thiogalactopyranoside (IPTG). LacI forms a complex with IPTG, blocking its interaction with the promoter region of the ACC gene. Similarly, FapR forms a complex with Malonyl-CoA. Therefore, this prevents FapR from affecting the expression of the gene. The gene expression proceeds successfully. By adjusting the values of these molecules, we could extract useful information about how this oscillatory behavior affects enzyme production.
Parameters & Equations
To analyze the mechanisms of our model, it is essential to specify and define the associated parameters. Careful selection of parameters is critical, as it is the foundation on which further evaluations are based. This involves a thorough examination of scientific papers, MATLAB Projects and model architectures.
When developing this model, we focused on understanding how oscillators and feedback loops work. For the FapR system to operate, a series of rules and abbreviations needs to be followed. The negative feedback loop is rather sensitive, and a small mistake can alter results. With this system that we proposed, our aim is to provide a solid foundation for one to build upon and understand the biological mechanism of E. coli.
Model Simulation
In our simulation, the plots strongly display the concentrations of each species in all systems, providing a complete picture of their dynamic interactions. Our system features a negative feedback loop, which is visually evident in the diagrams through discrete oscillations. These oscillations, characteristic of dynamical systems, reflect periodic variations in the concentrations of species over time. By changing the hill coefficient (n), oscillations became more frequent with its increase and less frequent with its decrease. Experimenting with other parameters helped when analyzing the plots and the results with the wet lab team.
Model Overview
Employing the MATLAB tool, SimBiology, we formulated a computational model of the PHA recovery model that focuses on biosafety composed of one distinct subsystem. This model was used primarily as an experimental design rather than a tool to extract analytical data.
- The mekR system. This regulatory mechanism focuses on the controlled production of the Lys enzyme using a MekR enzyme in combination with the MEK substrate. The interaction between MekR and MEK leads to the formation of a complex, which then serves as a potent activator for the transcriptional process of the lys gene. This tuned regulatory cascade highlights a key mechanism in cellular homeostasis.
Parameters & Equations
The parameters used in the construction of this system's model are depicted below:
The equations used in this model are similar to the equations of araC system of P. putida Model.
Model Simulation
This simlutation helped our team observe the Lys enzyme regulation. Here, the plots display the concentrations of each species in this system, providing a complete picture of their dynamic interactions.
Discussion
Despite the limitations associated with the modeling process, the result is encouraging, as the resulting models provided useful information. The models demonstrated a good level of efficiency, providing powerful data that enriched our understanding of the project’s biological phenomena. By utilizing these data, we were able to provide valuable insights to our colleagues in the wet lab, aiding their experimental design. In particular, the implementation of the FapR negative feedback loop emerged as the most effective regulatory mechanism, showing effectiveness in understanding regulating system behavior. The results highlight the potential of our model as a fundamental structure for constructing more complex and integrated models for the microbial consortium. Eventually, our ambition is to reflect the complex reality of the oPHAelia microbial consortium within our simulations, trying to provide future iGEM teams with a strong basis in their attempts to work with similar complex systems.
You can access our complete modeling documentation here:
References
- Chandran, D., Copeland, W. B., Sleight, S. C., & Sauro, H. M. (2008). Mathematical modeling and synthetic biology. Drug discovery today. Disease models, 5(4), 299–309. https://doi.org/10.1016/j.ddmod.2009.07.002
- Zheng, Y., & Sriram, G. (2010). Mathematical modeling: bridging the gap between concept and realization in synthetic biology. Journal of biomedicine & biotechnology, 2010, 541609. https://doi.org/10.1155/2010/541609
- Liu, D., & Zhang, F. (2018). Metabolic Feedback Circuits Provide Rapid Control of Metabolite Dynamics. ACS Synthetic Biology, 7(2), 347–356. doi:10.1021/acssynbio.7b00342
- Selinger, Saxena, Cheung, Church, & Rosenow. Global RNA half-life analysis in Escherichia coli reveals positional patterns of transcript degradation. Retrieved October 11, 2023, from https://arep.med.harvard.edu/rna_decay/
- Sato, S., Kanazawa, H. & Tsuge, T. Expression and characterization of (R)-specific enoyl coenzyme A hydratases making a channeling route to polyhydroxyalkanoate biosynthesis in Pseudomonas putida . Appl Microbiol Biotechnol 90, 951–959 (2011). https://doi.org/10.1007/s00253-011-3150-5
- Trautwein, G., & Gerischer, U. (2001, February). Effects Exerted by Transcriptional Regulator PcaU from Acinetobacter sp. Strain ADP1. Journal of Bacteriology, 183(3), 873–881. https://doi.org/10.1128/jb.183.3.873-881.2001
- Zimmermann, T., Sorg, T., Siehler, S. Y., & Gerischer, U. (2009, April 15). Role of Acinetobacter baylyi Crc in Catabolite Repression of Enzymes for Aromatic Compound Catabolism. Journal of Bacteriology, 191(8), 2834–2842. https://doi.org/10.1128/jb.00817-08
- Wang, X., Lin, L., Dong, J., Ling, J., Wang, W., Wang, H., … Yu, X. (2018). Simultaneous improvements of cell growth and polyhydroxyalkanoate production in Pseudomonas from a lignin derivative for lignin consolidated bioprocessing. Applied and Environmental Microbiology, AEM.01469–18. doi:10.1128/aem.01469-18
- J. Rawlings, B. (1998). Biosynthesis of fatty acids and related metabolites. Natural Product Reports, 15(3), 275. doi:10.1039/a815275y
- Wang, Q., Tappel, R. C., Zhu, C., & Nomura, C. T. (2011). Development of a New Strategy for Production of Medium-Chain-Length Polyhydroxyalkanoates by Recombinant Escherichia coli via Inexpensive Non-Fatty Acid Feedstocks. Applied and Environmental Microbiology, 78(2), 519–527. doi:10.1128/aem.07020-11
- Thomas, B. S., Okamoto, K., Bankowski, M. J., & Seto, T. B. (2013). A Lethal Case of Pseudomonas putida Bacteremia Due to Soft Tissue Infection. Infectious diseases in clinical practice (Baltimore, Md.), 21(3), 147–213. https://doi.org/10.1097/IPC.0b013e318276956b
- Ballerstedt, H., Volkers, R. J., Mars, A. E., Hallsworth, J. E., dos Santos, V. A., Puchalka, J., van Duuren, J., Eggink, G., Timmis, K. N., de Bont, J. A., & Wery, J. (2007). Genomotyping of Pseudomonas putida strains using P. putida KT2440-based high-density DNA microarrays: implications for transcriptomics studies. Applied microbiology and biotechnology, 75(5), 1133–1142. https://doi.org/10.1007/s00253-007-0914-z
- Loeschcke, A., & Thies, S. (2015). Pseudomonas putida-a versatile host for the production of natural products. Applied microbiology and biotechnology, 99(15), 6197–6214. https://doi.org/10.1007/s00253-015-6745-4
- Bonner, W. M., & Bloch, K. (1972). Purification and properties of fatty acyl thioesterase I from Escherichia coli. Journal of Biological Chemistry, 247(10), 3123-3133.
- Yang, J., Wang, G., Ng, T. B., Lin, J., & Ye, X. (2016). Laccase Production and Differential Transcription of Laccase Genes in Cerrena sp. in Response to Metal Ions, Aromatic Compounds, and Nutrients. Frontiers in microbiology, 6, 1558. https://doi.org/10.3389/fmicb.2015.01558
- Janusz, G., Pawlik, A., Świderska-Burek, U., Polak, J., Sulej, J., Jarosz-Wilkołazka, A., & Paszczyński, A. (2020). Laccase Properties, Physiological Functions, and Evolution. International journal of molecular sciences, 21(3), 966. https://doi.org/10.3390/ijms21030966
- Stricker, J., Cookson, S., Bennett, M. et al. A fast, robust and tunable synthetic gene oscillator. Nature 456, 516–519 (2008). https://doi.org/10.1038/nature07389
- Daniel, R., Rubens, J., Sarpeshkar, R. et al. Synthetic analog computation in living cells. Nature 497, 619–623 (2013). https://doi.org/10.1038/nature12148
- Megerle, J. A., Fritz, G., Gerland, U., Jung, K., & Rädler, J. O. (2008). Timing and dynamics of single cell gene expression in the arabinose utilization system. Biophysical journal, 95(4), 2103–2115. https://doi.org/10.1529/biophysj.107.127191
- Shraddha, Shekher, R., Sehgal, S., Kamthania, M., & Kumar, A. (2011). Laccase: microbial sources, production, purification, and potential biotechnological applications. Enzyme research, 2011, 217861. https://doi.org/10.4061/2011/217861
- Uri, M. (n.d.). Half-life of AraC protein in vivo. Half-life of AraC protein in vivo - Bacteria Escherichia coli - BNID 105138. https://bionumbers.hms.harvard.edu/bionumber.aspx?id=105138&ver=1&trm=araC%2BmRNA&org=
- Leon Y Chan, Christopher F Mugler, Stephanie Heinrich, Pascal Vallotton, Karsten Weis (2018) Non-invasive measurement of mRNA decay reveals translation initiation as the major determinant of mRNA stability eLife, https://doi.org/10.7554/eLife.32536
- Schujman, G. E., Paoletti, L., Grossman, A. D., & de Mendoza, D. (2003). FapR, a Bacterial Transcription Factor Involved in Global Regulation of Membrane Lipid Biosynthesis. Developmental Cell, 4(5), 663–672. doi:10.1016/s1534-5807(03)00123-0
- Heath, R. J., White, S. W., & Rock, C. O. (2001). Lipid biosynthesis as a target for antibacterial agents. Progress in Lipid Research, 40(6), 467–497. doi:10.1016/s0163-7827(01)00012-1
- Parameters of the lactose repressor according to mathematical model. BioNumbers. Retrieved October 11, 2023, from https://bionumbers.hms.harvard.edu/bionumber.aspx?id=106217&ver=2&trm=LacI&org=
- Graf, N., Altenbuchner, J. Functional characterization and application of a tightly regulated MekR/P mekA expression system in Escherichia coli and Pseudomonas putida . Appl Microbiol Biotechnol 97, 8239–8251 (2013). https://doi.org/10.1007/s00253-013-5030-7
- Hoffmann, J., & Altenbuchner, J. (2015). Functional Characterization of the Mannitol Promoter of Pseudomonas fluorescens DSM 50106 and Its Application for a Mannitol-Inducible Expression System for Pseudomonas putida KT2440. PLOS ONE, 10(7), e0133248. doi:10.1371/journal.pone.0133248
- Inoue, I., & Rechsteiner, M. (1994). On the relationship between the metabolic and thermodynamic stabilities of T4 lysozymes. Measurements in Escherichia coli. The Journal of biological chemistry, 269(46), 29241–29246. https://doi.org/10.1016/S0021-9258(19)62036-3
- Bergamo, A., Gerdol, M., Pallavicini, A., Greco, S., Schepens, I., Hamelin, R., Armand, F., et al. (2019). Lysozyme-Induced Transcriptional Regulation of TNF-α Pathway Genes in Cells of the Monocyte Lineage. International Journal of Molecular Sciences, 20(21), 5502. MDPI AG. Retrieved from http://dx.doi.org/10.3390/ijms20215502
- Lin, J., & Amir, A. (2018, October 29). Homeostasis of protein and mRNA concentrations in growing cells. Nature Communications; Nature Portfolio. https://doi.org/10.1038/s41467-018-06714-zBergamo, A., Gerdol, M., Pallavicini, A., Greco, S., Schepens, I., Hamelin, R., Armand, F., et al. (2019). Lysozyme-Induced Transcriptional Regulation of TNF-α Pathway Genes in Cells of the Monocyte Lineage. International Journal of Molecular Sciences, 20(21), 5502. MDPI AG. Retrieved from http://dx.doi.org/10.3390/ijms20215502
- MIS for BCMP 309 (drug delivery and clearance example based on SimBiology Intro file by JW). Modeling with SimBiology. Quantitative Methods Boot Camp.
- Semsey, S., Jauffred, L., Csiszovszki, Z., Erdossy, J., Stéger, V., Hansen, S., & Krishna, S. (2013). The effect of LacI autoregulation on the performance of the lactose utilization system in Escherichia coli. Nucleic acids research, 41(13), 6381–6390. https://doi.org/10.1093/nar/gkt351