A mathematical model of P53 gene regulatory networks under radiotherapy

Control of apoptosis by SMAR1

The nuclear matrix associated protein SMAR1 is sensitive to p53 and acts as a stress inducer as well as a regulator in the p53 regulatory network. Depending on the amount of stress SMAR1 stimulates, it can drive the p53 dynamics in the system to various dynamical states which correspond to various cellular states. The behavior of p53 in these dynamical states is found to be multifractal, due to the mostly long range correlations and large scale fluctuations imparted by stress. This fractal behavior is exhibited in the topological properties of the networks constructed from these dynamical states, and is a signature of self-organization to optimize information flow in the dynamics. The assortativity found in these networks is due to perturbation induced by stress, and indicates that the hubs in the time series play a significant role in stress management. SMAR1 can also regulate apoptosis in the presence of HDAC1, depending on the stress induced by it.

Kinetic aspects of enzyme-mediated repair of DNA single-strand breaks

In cells and bacteria, DNA can be damaged in different ways. The efficient damage repair, mediated by various enzymes, is crucial for their survival. Most frequently, the damage is reduced to single-strand breaks. In human cells, according to the experiments, the repair of such breaks can mechanistically be divided into four steps including (i) the break detection, (ii) processing of damaged ends, (iii) gap filling, and (iv) ligation of unbound ends of the broken strand. The first and second steps run in parallel while the third and fourth steps are sequential. The author proposes a kinetic model describing these steps. It allows one to understand the likely dependence of the number of breaks in different states on enzyme concentrations. The dependence of these concentrations on the rate of the formation of breaks can be understood as well. In addition, the likely role of unzipping and zipping of the fragments of broken ends of the strand in the ligation step has been scrutinized taking the specifics of binding of DNA stands into account.

Non-linear feedback control of the p53 protein-mdm2 inhibitor system using the derivative-free non-linear Kalman filter

It is proven that the model of the p53-mdm2 protein synthesis loop is a differentially flat one and using a diffeomorphism (change of state variables) that is proposed by differential flatness theory it is shown that the protein synthesis model can be transformed into the canonical (Brunovsky) form. This enables the design of a feedback control law that maintains the concentration of the p53 protein at the desirable levels. To estimate the non-measurable elements of the state vector describing the p53-mdm2 system dynamics, the derivative-free non-linear Kalman filter is used. Moreover, to compensate for modelling uncertainties and external disturbances that affect the p53-mdm2 system, the derivative-free non-linear Kalman filter is re-designed as a disturbance observer. The derivative-free non-linear Kalman filter consists of the Kalman filter recursion applied on the linearised equivalent of the protein synthesis model together with an inverse transformation based on differential flatness theory that enables to retrieve estimates for the state variables of the initial non-linear model. The proposed non-linear feedback control and perturbations compensation method for the p53-mdm2 system can result in more efficient chemotherapy schemes where the infusion of medication will be better administered.

Change detection in the dynamics of an intracellular protein synthesis model using nonlinear Kalman filtering

A method for early diagnosis of parametric changes in intracellular protein synthesis models (e.g. the p53 protein - mdm2 inhibitor model) is developed with the use of a nonlinear Kalman Filtering approach (Derivative-free nonlinear Kalman Filter) and of statistical change detection methods. The intracellular protein synthesis dynamic model is described by a set of coupled nonlinear differential equations. It is shown that such a dynamical system satisfies differential flatness properties and this allows to transform it, through a change of variables (diffeomorphism), to the so-called linear canonical form. For the linearized equivalent of the dynamical system, state estimation can be performed using the Kalman Filter recursion. Moreover, by applying an inverse transformation based on the previous diffeomorphism it becomes also possible to obtain estimates of the state variables of the initial nonlinear model. By comparing the output of the Kalman Filter (which is assumed to correspond to the undistorted dynamical model) with measurements obtained from the monitored protein synthesis system, a sequence of differences (residuals) is obtained. The statistical processing of the residuals with the use of x2 change detection tests, can provide indication within specific confidence intervals about parametric changes in the considered biological system and consequently indications about the appearance of specific diseases (e.g. malignancies).

Kinetic theory approach to modeling of cellular repair mechanisms under genome stress

Under acute perturbations from outer environment, a normal cell can trigger cellular self-defense mechanism in response to genome stress. To investigate the kinetics of cellular self-repair process at single cell level further, a model of DNA damage generating and repair is proposed under acute Ion Radiation (IR) by using mathematical framework of kinetic theory of active particles (KTAP). Firstly, we focus on illustrating the profile of Cellular Repair System (CRS) instituted by two sub-populations, each of which is made up of the active particles with different discrete states. Then, we implement the mathematical framework of cellular self-repair mechanism, and illustrate the dynamic processes of Double Strand Breaks (DSBs) and Repair Protein (RP) generating, DSB-protein complexes (DSBCs) synthesizing, and toxins accumulating. Finally, we roughly analyze the capability of cellular self-repair mechanism, cellular activity of transferring DNA damage, and genome stability, especially the different fates of a certain cell before and after the time thresholds of IR perturbations that a cell can tolerate maximally under different IR perturbation circumstances.

Harmonic oscillations in homeostatic controllers: Dynamics of the p53 regulatory system

Homeostatic mechanisms are essential for the protection and adaptation of organisms in a changing and challenging environment. Previously, we have described molecular mechanisms that lead to robust homeostasis/adaptation under inflow or outflow perturbations. Here we report that harmonic oscillations occur in models of such homeostatic controllers and that a close relationship exists between the control of the p53/Mdm2 system and that of a homeostatic inflow controller. This homeostatic control model of the p53 system provides an explanation why large fluctuations in the amplitude of p53/Mdm2 oscillations may arise as part of the homeostatic regulation of p53 by Mdm2 under DNA-damaging conditions. In the presence of DNA damage p53 is upregulated, but is subject to a tight control by Mdm2 and other factors to avoid a premature apoptotic response of the cell at low DNA damage levels. One of the regulatory steps is the Mdm2-mediated degradation of p53 by the proteasome. Oscillations in the p53/Mdm2 system are considered to be part of a mechanism by which a cell decides between cell cycle arrest/DNA repair and apoptosis. In the homeostatic inflow control model, harmonic oscillations in p53/Mdm2 levels arise when the binding strength of p53 to degradation complexes increases. Due to the harmonic character of the oscillations rapid fluctuating noise can lead, as experimentally observed, to large variations in the amplitude of the oscillation but not in their period, a behavior which has been difficult to simulate by deterministic limit-cycle models. In conclusion, the oscillatory response of homeostatic controllers may provide new insights into the origin and role of oscillations observed in homeostatically controlled molecular networks.

Dynamic modeling of cellular response to DNA damage based on p53 stress response networks

Under acute perturbations from the outside, cells can trigger self-defensive mechanisms to fight against genome stress. To investigate the cellular response to continuous ion radiation (IR), a dynamic model for p53 stress response networks at the cellular level is proposed. The model can successfully be used to simulate the dynamic processes of double-strand breaks (DSBs) generation and their repair, switch-like ataxia telangiectasia mutated (ATM) activation, oscillations occurring in the p53-MDM2 feedback loop, as well as toxins elimination triggered by p53 stress response networks. Especially, the model can predict the plausible outcomes of cellular response under different IR dose regimes.

A systems biology view of cancer

In order to understand how a cancer cell is functionally different from a normal cell it is necessary to assess the complex network of pathways involving gene regulation, signaling, and cell metabolism, and the alterations in its dynamics caused by the several different types of mutations leading to malignancy. Since the network is typically complex, with multiple connections between pathways and important feedback loops, it is crucial to represent it in the form of a computational model that can be used for a rigorous analysis. This is the approach of systems biology, made possible by new -omics data generation technologies. The goal of this review is to illustrate this approach and its utility for our understanding of cancer. After a discussion of recent progress using a network-centric approach, three case studies related to diagnostics, therapy, and drug development are presented in detail. They focus on breast cancer, B-cell lymphomas, and colorectal cancer. The discussion is centered on key mathematical and computational tools common to a systems biology approach.

Mathematical model of a network of interaction between p53 and Bcl-2 during genotoxic-induced apoptosis

Ionizing radiation like UV light, gamma and X rays, can produce genotoxic damage in fibroblast cells. This injury can be reverted by activation of specific nuclear molecules. However, intense genotoxic damage induces the activation of the p53 dependent apoptotic pathway. Activated nuclear p53 has the role of a transcription factor that switches on the transcription of the Puma protein, which once released into the cytoplasm leads to the activation of a network of chemical processes that produces cell death. This network is built up with the chemical interaction between pro-apoptotic p53, Puma and Bax proteins and the anti-apoptotic Bcl2 and Bcl-x(L) proteins. In this work we present a mathematical model of this modular network under different regimes of Puma release into the cytoplasm. In this model we use a set of coupled nonlinear differential equations, which is solved by numerical methods, to determine the nature of the equilibrium points of the dynamical system. With this information we construct the phase portrait associated with Puma induced apoptosis and we found that once the genotoxic injury is produced, Bax levels continuously increase. In this work we propose that a possible mechanism for the control of apoptosis is the Bax level in the cytoplasm, i.e., Bax is continuously released into the cytoplasm, even in absence of Puma, according to the kinetic properties of the set of chemical interactions in which the Bcl2/Bax complex takes part. If the damage is reverted before the Bax level reaches a threshold value the apoptotic process is stopped, otherwise the process goes on until cell death.