A theoretical exploration of birhythmicity in the p53-Mdm2 network.
Bottom Line:
To do so, we reduced this model, in a first step, into a 3-dimensional piecewise linear differential model where the Hill functions have been approximated by step functions, and, in a second step, into a 2-dimensional piecewise linear differential model by setting one autonomous variable as a constant in each domain of the phase space.We find that two features related to the phase space structure of the system are at the origin of the birhythmic behavior: the existence of two embedded cycles in the transition graph of the reduced models; the presence of a bypass in the orbit of the large amplitude oscillatory regime of low frequency.From a methodological point of view, this approach greatly facilitates the computational analysis of complex oscillatory behavior and could represent a valuable tool to explore mathematical models of biological rhythms showing sufficiently steep nonlinearities.
View Article:
PubMed Central - PubMed
Affiliation: COMORE Project-team, INRIA Sophia Antipolis, Sophia Antipolis, France. wabou@sophia.inria.fr
ABSTRACT
Show MeSH
Experimental observations performed in the p53-Mdm2 network, one of the key protein modules involved in the control of proliferation of abnormal cells in mammals, revealed the existence of two frequencies of oscillations of p53 and Mdm2 in irradiated cells depending on the irradiation dose. These observations raised the question of the existence of birhythmicity, i.e. the coexistence of two oscillatory regimes for the same external conditions, in the p53-Mdm2 network which would be at the origin of these two distinct frequencies. A theoretical answer has been recently suggested by Ouattara, Abou-Jaoudé and Kaufman who proposed a 3-dimensional differential model showing birhythmicity to reproduce the two frequencies experimentally observed. The aim of this work is to analyze the mechanisms at the origin of the birhythmic behavior through a theoretical analysis of this differential model. To do so, we reduced this model, in a first step, into a 3-dimensional piecewise linear differential model where the Hill functions have been approximated by step functions, and, in a second step, into a 2-dimensional piecewise linear differential model by setting one autonomous variable as a constant in each domain of the phase space. We find that two features related to the phase space structure of the system are at the origin of the birhythmic behavior: the existence of two embedded cycles in the transition graph of the reduced models; the presence of a bypass in the orbit of the large amplitude oscillatory regime of low frequency. Based on this analysis, an experimental strategy is proposed to test the existence of birhythmicity in the p53-Mdm2 network. From a methodological point of view, this approach greatly facilitates the computational analysis of complex oscillatory behavior and could represent a valuable tool to explore mathematical models of biological rhythms showing sufficiently steep nonlinearities. |
Related In:
Results -
Collection
getmorefigures.php?uid=PMC3038873&req=5
Mentions: Projection of the two oscillatory regimes composing the birhythmic behavior of Model 2 in the plane (Mn,P). The dashed lines represent the thresholds of the step functions: KMn and KMc for P, KP for Mn. The parameter values are the same as for Figure 3 (KMn = 0.1 nM, KP = 0.2 nM) except KMc = 0.6 nM. The period of the large amplitude oscillatory regime is significantly longer than the period of the small amplitude oscillatory regime (see Figure 8A and 8C). |
View Article: PubMed Central - PubMed
Affiliation: COMORE Project-team, INRIA Sophia Antipolis, Sophia Antipolis, France. wabou@sophia.inria.fr