How Molecular Should Your Molecular Model Be?: On the Level of Molecular Detail Required to Simulate Biological Networks in Systems and Synthetic Biology
Abstract
The recent advance of genetic studies and the rapid accumulation of molecular data, together with the increasing performance of computers, led researchers to design more and more detailed mathematical models of biological systems. Many modeling approaches rely on ordinary differential equations (ODE) which are based on standard enzyme kinetics. Michaelis-Menten and Hill functions are indeed commonly used in dynamical models in systems and synthetic biology because they provide the necessary nonlinearity to make the dynamics nontrivial (i.e., limit-cycle oscillations or multistability). For most of the systems modeled, the actual molecular mechanism is unknown, and the enzyme equations should be regarded as phenomenological. In this chapter, we discuss the validity and accuracy of these approximations. In particular, we focus on the validity of the Michaelis-Menten function for open systems and on the use of Hill kinetics to describe transcription rates of regulated genes. Our discussion is illustrated by numerical simulations of prototype systems, including the Repressilator (a genetic oscillator) and the Toggle Switch model (a bistable system). We systematically compare the results obtained with the compact version (based on Michaelis-Menten and Hill functions) with its corresponding developed versions (based on "elementary" reaction steps and mass action laws). We also discuss the use of compact approaches to perform stochastic simulations (Gillespie algorithm). On the basis of these results, we argue that using compact models is suitable to model qualitatively biological.