RESUMO
We propose and solve analytically a stochastic model for the dynamics of a binary biological switch, defined as a DNA unit with two mutually exclusive configurations, each one triggering the expression of a different gene. Such a device has the potential to be used as a memory unit for biological computing systems designed to operate in noisy environments. We discuss a recent implementation of this switch in living cells, the recombinase addressable data (RAD) module. In order to understand the behavior of a RAD module we compute the exact time-dependent joint distribution of the two expressed genes starting in one state and evolving to another asymptotic state. We consider two operating regimes of the RAD module, a fast and a slow stochastic switching regime. The fast regime is aggregative and produces unimodal distributions, whereas the slow regime is separative and produces bimodal distributions. Both regimes can serve to prepare pure memory states when all cells are expressing the same gene. The slow regime can also separate mixed states by producing two subpopulations, each one expressing a different gene. Compared to the genetic toggle switch based on positive feedback, the RAD module ensures more rapid memory operations for the same quality of the separation between binary states. Our model provides a simplified phenomenological framework for studying RAD memory devices and our analytic solution can be further used to clarify theoretical concepts in biocomputation and for optimal design in synthetic biology.
Assuntos
Modelos Genéticos , Expressão Gênica/genética , Processos Estocásticos , Biologia SintéticaRESUMO
In this manuscript, we propose a mathematical framework to couple transcription and translation in which mRNA production is described by a set of master equations, while the dynamics of protein density is governed by a random differential equation. The coupling between the two processes is given by a stochastic perturbation whose statistics satisfies the master equations. In this approach, from the knowledge of the analytical time-dependent distribution of mRNA number, we are able to calculate the dynamics of the probability density of the protein population.