Modeling the role of p53 pulses in DNA damage- induced cell death decision

Comprehensive Bioinformatic Investigation of TP53 Dysregulation in Diverse Cancer Landscapes

P53 overexpression plays a critical role in cancer pathogenesis by disrupting the intricate regulation of cellular proliferation. Despite its firmly established function as a tumor suppressor, elevated p53 levels can paradoxically contribute to tumorigenesis, influenced by factors such as exposure to carcinogens, genetic mutations, and viral infections. This phenomenon is observed across a spectrum of cancer types, including bladder (BLCA), ovarian (OV), cervical (CESC), cholangiocarcinoma (CHOL), colon adenocarcinoma (COAD), diffuse large B-cell lymphoma (DLBC), esophageal carcinoma (ESCA), head and neck squamous cell carcinoma (HNSC), kidney chromophobe (KICH), kidney renal clear cell carcinoma (KIRC), liver hepatocellular carcinoma (LIHC), lung adenocarcinoma (LUAD), lung squamous cell carcinoma (LUSC), and uterine corpus endometrial carcinoma (UCEC). This broad spectrum of cancers is often associated with increased aggressiveness and recurrence risk. Effective therapeutic strategies targeting tumors with p53 overexpression require a comprehensive approach, integrating targeted interventions aimed at the p53 gene with conventional modalities such as chemotherapy, radiation therapy, and targeted drugs. In this extensive study, we present a detailed analysis shedding light on the multifaceted role of TP53 across various cancers, with a specific emphasis on its impact on disease-free survival (DFS). Leveraging data from the TCGA database and the GTEx dataset, along with GEPIA, UALCAN, and STRING, we identify TP53 overexpression as a significant prognostic indicator, notably pronounced in prostate adenocarcinoma (PRAD). Supported by compelling statistical significance (p < 0.05), our analysis reveals the distinct influence of TP53 overexpression on DFS outcomes in PRAD. Additionally, graphical representations of overall survival (OS) underscore the notable disparity in OS duration between tumors exhibiting elevated TP53 expression (depicted by the red line) and those with lower TP53 levels (indicated by the blue line). The hazard ratio (HR) further emphasizes the profound impact of TP53 on overall survival. Moreover, our investigation delves into the intricate TP53 protein network, unveiling genes exhibiting robust positive correlations with TP53 expression across 13 out of 27 cancers. Remarkably, negative correlations emerge with pivotal tumor suppressor genes. This network analysis elucidates critical proteins, including SIRT1, CBP, p300, ATM, DAXX, HSP 90-alpha, Mdm2, RPA70, 14-3-3 protein sigma, p53, and ASPP2, pivotal in regulating cell cycle dynamics, DNA damage response, and transcriptional regulation. Our study underscores the paramount importance of deciphering TP53 dynamics in cancer, providing invaluable insights into tumor behavior, disease-free survival, and potential therapeutic avenues.

Coordination of MAPK and p53 dynamics in the cellular responses to DNA damage and oxidative stress

In response to different cellular stresses, the transcription factor p53 undergoes different dynamics. p53 dynamics, in turn, control cell fate. However, distinct stresses can generate the same p53 dynamics but different cell fate outcomes, suggesting integration of dynamic information from other pathways is important for cell fate regulation. To determine how MAPK activities affect p53-mediated responses to DNA breaks and oxidative stress, we simultaneously tracked p53 and either ERK, JNK, or p38 activities in single cells. While p53 dynamics were comparable between the stresses, cell fate outcomes were distinct. Combining MAPK dynamics with p53 dynamics was important for distinguishing between the stresses and for generating temporal ordering of cell fate pathways. Furthermore, cross-talk between MAPKs and p53 controlled the balance between proliferation and cell death. These findings provide insight into how cells integrate signaling pathways with distinct temporal patterns of activity to encode stress specificity and drive different cell fate decisions.

A common pathway to cancer: Oncogenic mutations abolish p53 oscillations

The tumor suppressor p53 oscillates in response to DNA double-strand breaks, a behavior that has been suggested to be essential to its anti-cancer function. Nearly all human cancers have genetic alterations in the p53 pathway; a number of these alterations have been shown to be oncogenic by experiment. These alterations include somatic mutations and copy number variations as well as germline polymorphisms. Intriguingly, they exhibit a mixed pattern of interactions in tumors, such as co-occurrence, mutual exclusivity, and paradoxically, mutual antagonism. Using a differential equation model of p53-Mdm2 dynamics, we employ Hopf bifurcation analysis to show that these alterations have a common mode of action, to abolish the oscillatory competence of p53, thereby, we suggest, impairing its tumor suppressive function. In this analysis, diverse genetic alterations, widely associated with human cancers clinically, have a unified mechanistic explanation of their role in oncogenesis.

Modeling of DNA Damage Repair and Cell Response in Relation to p53 System Exposed to Ionizing Radiation

Repair of DNA damage induced by ionizing radiation plays an important role in the cell response to ionizing radiation. Radiation-induced DNA damage also activates the p53 system, which determines the fate of cells. The kinetics of repair, which is affected by the cell itself and the complexity of DNA damage, influences the cell response and fate via affecting the p53 system. To mechanistically study the influences of the cell response to different LET radiations, we introduce a new repair module and a p53 system model with NASIC, a Monte Carlo track structure code. The factors determining the kinetics of the double-strand break (DSB) repair are modeled, including the chromosome environment and complexity of DSB. The kinetics of DSB repair is modeled considering the resection-dependent and resection-independent compartments. The p53 system is modeled by simulating the interactions among genes and proteins. With this model, the cell responses to low- and high-LET irradiation are simulated, respectively. It is found that the kinetics of DSB repair greatly affects the cell fate and later biological effects. A large number of DSBs and a slow repair process lead to severe biological consequences. High-LET radiation induces more complex DSBs, which can be repaired by slow processes, subsequently resulting in a longer cycle arrest and, furthermore, apoptosis and more secreting of TGFβ. The Monte Carlo track structure simulation with a more realistic repair module and the p53 system model developed in this study can expand the functions of the NASIC code in simulating mechanical radiobiological effects.

Dual roles of SIRT1 in the BAX switch through the P53 module: A mathematical modeling study

SIRT1 is a multifunctional deacetylase that participates in a variety of cellular physiological processes to cope with stress. The anticancer protein P53 is an important target of SIRT1. It has been found that SIRT1 is involved in apoptosis by regulating the activity and intracellular location of P53. Moreover, P53-dependent apoptosis is inseparable from the BCL-2 protein family. Among the members of this family, BAX’s switching dynamics may play a key role in apoptosis. Therefore, a challenging question arises: what effect does SIRT1 have on the BAX switch? To answer this question, we built a small-scale protein network model. Through computer simulation, the properties of SIRT1 that on the one hand promote and on the other inhibit apoptosis are revealed. We found that the opening time of the BAX switch will be delayed in the case of either SIRT1 excess or deficiency. Similarly, the stimulus threshold required for apoptosis will also increase in the above two scenarios. Thereby, we proposed that SIRT1 has an optimal content at which the probability of apoptosis is greatest. In addition, P53 oscillation requires the concentration of SIRT1 to be higher than the optimal value. This work may be helpful both experimentally and clinically.

Modeling the neuro-protection of theaflavic acid from black tea and its synergy with nimodipine via mitochondria apoptotic pathway

Ischemic stroke presents a leading cause of mortality and morbidity worldwide. Theaflavic acid (TFA) is a theaflavin isolated from black tea that exerts a potentially neuro-protective effect. However, the dynamic properties of TFA-mediated protection remain largely unknown. In the current study, we evaluated the function of TFA in the mitochondria apoptotic pathway using mathematical modeling. We found that TFA-enhanced B-cell lymphoma 2 (Bcl-2) overexpression can theoretically give rise to bistability. The bistability is highly robust against parametric stochasticity while also conferring considerable variability in survival threshold. Stochastic simulations faithfully match the TFA dose response pattern seen in experimental studies. In addition, we identified a dose- and time-dependent synergy between TFA and nimodipine, a clinically used neuro-protective drug. This synergistic effect was enhanced by bistability independent of temporal factors. Precise application of pulsed doses of TFA can also promote survival compared with sustained TFA treatment. These data collectively demonstrate that TFA treatment can give rise to bistability and that synergy between TFA and nimodipine may offer a promising strategy for developing therapeutic neuro-protection against ischemic stroke.

Protein stability of p53 targets determines their temporal expression dynamics in response to p53 pulsing

Oscillations in p53 expression are critical for regulating the cellular response to DNA damage. Hanson et al. show that the relationship between p53 pulse frequency and target mRNA and protein decay rates regulates stress response pathway dynamics and function.

In response to DNA damage, the transcription factor p53 accumulates in a series of pulses. While p53 dynamics play a critical role in regulating stress responses, how p53 pulsing affects target protein expression is not well understood. Recently, we showed that p53 pulses generate diversity in target mRNA expression dynamics; however, given that mRNA and protein expression are not necessarily well correlated, it remains to be determined how p53 pulses impact target protein expression. Using computational and experimental approaches, we show that target protein decay rates filter p53 pulses: Distinct target protein expression dynamics are generated depending on the relationship between p53 pulse frequency and target mRNA and protein stability. Furthermore, by mutating the targets MDM2 and PUMA to alter their stabilities, we show that downstream pathways are sensitive to target protein decay rates. This study delineates the mechanisms by which p53 dynamics play a crucial role in orchestrating the timing of events in the DNA damage response network.

Recent progress and open challenges in modeling p53 dynamics in single cells

In mammalian cells, the tumor suppressor p53 is activated upon a variety of cellular stresses and ensures an appropriate response ranging from arrest and repair to the induction of senescence and apoptosis. Quantitative measurements in individual living cells showed stimulus-dependent dynamics of p53 accumulation upon stress induction. Due to the complexity of the underlying biochemical interactions, mathematical models were indispensable for understanding the topology of the network regulating p53 dynamics. Recent work provides furhter insights into the causes of heterogeneous responses in individual cells, the rewiring of the network in response to different inputs and the role of the downstream processes in determining the cellular fate upon stress.

Modeling the heterogeneity of p53 dynamics in DNA damage response

The tumor suppressor p53 can be activated by DNA damage and exhibits undamped pulses. Recent reports have demonstrated a non-threshold mechanism for p53 dynamics. However, no related theoretical studies have been proposed. Here, we constructed a refined DNA damage repair model that incorporated both intrinsic and extrinsic DNA lesions. We proposed that the basal DNA damage may trigger significant fractions of p53 pulses. We also reproduced the heterogeneity of p53 dynamics in experiments. The number of p53 pulses showed no correlations with DNA damage. We also replicated the linear correlation between DNA damage and the probability of igniting a pulse. Our model has unraveled the heterogeneous p53 responses.

Bifurcation in Cell Cycle Dynamics Regulated by p53

We study the regulating mechanism of p53 on the properties of cell cycle dynamics in the light of the proposed model of interacting p53 and cell cycle networks via p53. Irradiation (IR) introduce to p53 compel p53 dynamics to suffer different phases, namely oscillating and oscillation death (stabilized) phases. The IR induced p53 dynamics undergo collapse of oscillation with collapse time Δt which depends on IR strength. The stress p53 via IR drive cell cycle molecular species MPF and cyclin dynamics to different states, namely, oscillation death, oscillations of periods, chaotic and sustain oscillation in their bifurcation diagram. We predict that there could be a critical Δtc induced by p53 via IRc, where, if Δt〈Δtc the cell cycle may come back to normal state, otherwise it will go to cell cycle arrest (apoptosis).

Investigation of oscillation accumulation triggered genetic switch in gene regulatory networks

Genetic oscillator motifs and genetic switch motifs are blocks of biochemical reaction networks, which are involved in the regulation of rhythms, cell cycle progression, signal processing and cell fate decision. These motifs often interact to constitute complex signal processing systems. There widely exists the oscillation accumulation triggered mechanism in gene regulatory networks, i.e., an oscillator promotes the accumulation of a signaling protein over a threshold value, and activates a switch. The structure can be fund in some important biological function switches, such as apoptosis and DNA repair. We propose the structure as an oscillation accumulation triggered genetic switch (OATGS). Through mathematical modeling and analysis, results show the OATGS with features of robustness to noise and triggered mode. In addition, we show the existence of OATGS features and triggered manner in p53 gene regulatory networks, and explain some of the p53 regulation process, such as counting mechanism and pulse shape. We speculate that OATGS with oscillation accumulation triggered manner is a new important biological function switch.

A plausible model for bimodal p53 switch in DNA damage response

p53 is a tumor suppressor and the p53 dynamics displays stimulus dependent patterns. Recent evidence suggests a bimodal p53 switch in cell fate decision. However, no theoretical studies have been proposed to investigate bimodal p53 induction. Here we constructed a model and showed that MDM2-p53 mRNA binding might contribute to bimodal p53 switch through an intrinsic positive feedback loop. Lower damage favored pulsing while monotonic increasing was generated with higher damage. Bimodal p53 dynamics was largely influenced by cellular MDM2 and elevated p53/MDM2 ratios with increasing etoposide favor mono-ubiquitination. Our model replicated recent experiments and provided potential insights into dynamic mechanisms of bimodal switch.

The p53 response in single cells is linearly correlated to the number of DNA breaks without a distinct threshold

BACKGROUND

The tumor suppressor protein p53 is activated by cellular stress. DNA double strand breaks (DSBs) induce the activation of the kinase ATM, which stabilizes p53 and activates its transcriptional activity. Single cell analysis revealed that DSBs induced by gamma irradiation trigger p53 accumulation in a series of pulses that vary in number from cell to cell. Higher levels of irradiation increase the number of p53 pulses suggesting that they arise from periodic examination of the damage by ATM. If damage persists, additional pulses of p53 are triggered. The threshold of damage required for activating a p53 pulse is unclear. Previous studies that averaged the response across cell populations suggested that one or two DNA breaks are sufficient for activating ATM and p53. However, it is possible that by averaging over a population of cells important features of the dependency between DNA breaks and p53 dynamics are missed.

RESULTS

Using fluorescent reporters we developed a system for following in individual cells the number of DSBs, the kinetics of repair and the p53 response. We found a large variation in the initial number of DSBs and the rate of repair between individual cells. Cells with higher number of DSBs had higher probability of showing a p53 pulse. However, there was no distinct threshold number of breaks for inducing a p53 pulse. We present evidence that the decision to activate p53 given a specific number of breaks is not entirely stochastic, but instead is influenced by both cell-intrinsic factors and previous exposure to DNA damage. We also show that the natural variations in the initial amount of p53, rate of DSB repair and cell cycle phase do not affect the probability of activating p53 in response to DNA damage.

CONCLUSIONS

The use of fluorescent reporters to quantify DNA damage and p53 levels in live cells provided a quantitative analysis of the complex interrelationships between both processes. Our study shows that p53 activation differs even between cells that have a similar number of DNA breaks. Understanding the origin and consequences of such variability in normal and cancerous cells is crucial for developing efficient and selective therapeutic interventions.

Bayesian parameter inference by Markov chain Monte Carlo with hybrid fitness measures: theory and test in apoptosis signal transduction network

When exact values of model parameters in systems biology are not available from experiments, they need to be inferred so that the resulting simulation reproduces the experimentally known phenomena. For the purpose, Bayesian statistics with Markov chain Monte Carlo (MCMC) is a useful method. Biological experiments are often performed with cell population, and the results are represented by histograms. On another front, experiments sometimes indicate the existence of a specific bifurcation pattern. In this study, to deal with both type of such experimental results and information for parameter inference, we introduced functions to evaluate fitness to both type of experimental results, named quantitative and qualitative fitness measures respectively. We formulated Bayesian formula for those hybrid fitness measures (HFM), and implemented it to MCMC (MCMC-HFM). We tested MCMC-HFM first for a kinetic toy model with a positive feedback. Inferring kinetic parameters mainly related to the positive feedback, we found that MCMC-HFM reliably infer them with both qualitative and quantitative fitness measures. Then, we applied the MCMC-HFM to an apoptosis signal transduction network previously proposed. For kinetic parameters related to implicit positive feedbacks, which are important for bistability and irreversibility of the output, the MCMC-HFM reliably inferred these kinetic parameters. In particular, some kinetic parameters that have the experimental estimates were inferred without these data and the results were consistent with the experiments. Moreover, for some parameters, the mixed use of quantitative and qualitative fitness measures narrowed down the acceptable range of parameters. Taken together, our approach could reliably infer the kinetic parameters of the target systems.

Mechanisms that enhance sustainability of p53 pulses

The tumor suppressor p53 protein shows various dynamic responses depending on the types and extent of cellular stresses. In particular, in response to DNA damage induced by γ-irradiation, cells generate a series of p53 pulses. Recent research has shown the importance of sustaining repeated p53 pulses for recovery from DNA damage. However, far too little attention has been paid to understanding how cells can sustain p53 pulses given the complexities of genetic heterogeneity and intrinsic noise. Here, we explore potential molecular mechanisms that enhance the sustainability of p53 pulses by developing a new mathematical model of the p53 regulatory system. This model can reproduce many experimental results that describe the dynamics of p53 pulses. By simulating the model both deterministically and stochastically, we found three potential mechanisms that improve the sustainability of p53 pulses: 1) the recently identified positive feedback loop between p53 and Rorα allows cells to sustain p53 pulses with high amplitude over a wide range of conditions, 2) intrinsic noise can often prevent the dampening of p53 pulses even after mutations, and 3) coupling of p53 pulses in neighboring cells via cytochrome-c significantly reduces the chance of failure in sustaining p53 pulses in the presence of heterogeneity among cells. Finally, in light of these results, we propose testable experiments that can reveal important mechanisms underlying p53 dynamics.

PDCD5-regulated cell fate decision after ultraviolet-irradiation-induced DNA damage

Programmed cell death 5 (PDCD5) is a human apoptosis-related molecule that is involved in both the cytoplasmic caspase-3 activity pathway (by regulating Bax translocation from cytoplasm to mitochondria) and the nuclear pathway (by interacting with Tip60). In this study, we developed a mathematical model of the PDCD5-regulated switching of the cell response from DNA repair to apoptosis after ultraviolet irradiation-induced DNA damage. We established the model by combining several hypotheses with experimental observations. Our simulations indicate that the ultimate cell response to DNA damage is dependent on a signal threshold mechanism, and the PDCD5 promotion of Bax translocation plays an essential role in PDCD5-regulated cell apoptosis. Furthermore, the model simulations revealed that PDCD5 nuclear translocation can attenuate cell apoptosis, and PDCD5 interactions with Tip60 can accelerate DNA damage-induced apoptosis, but the final cell fate decision is insensitive to the PDCD5-Tip60 interaction. These results are consistent with experimental observations. The effect of recombinant human PDCD5 was also investigated and shown to sensitize cells to DNA damage by promoting caspase-3 activity.

Rigor of cell fate decision by variable p53 pulses and roles of cooperative gene expression by p53

Upon DNA damage, the cell fate decision between survival and apoptosis is largely regulated by p53-related networks. Recent experiments found a series of discrete p53 pulses in individual cells, which led to the hypothesis that the cell fate decision upon DNA damage is controlled by counting the number of p53 pulses. Under this hypothesis, Sun et al. (2009) modeled the Bax activation switch in the apoptosis signal transduction pathway that can rigorously "count" the number of uniform p53 pulses. Based on experimental evidence, here we use variable p53 pulses with Sun et al.'s model to investigate how the variability in p53 pulses affects the rigor of the cell fate decision by the pulse number. Our calculations showed that the experimentally anticipated variability in the pulse sizes reduces the rigor of the cell fate decision. In addition, we tested the roles of the cooperativity in PUMA expression by p53, finding that lower cooperativity is plausible for more rigorous cell fate decision. This is because the variability in the p53 pulse height is more amplified in PUMA expressions with more cooperative cases.

Modeling the basal dynamics of p53 system

BACKGROUND

The tumor suppressor p53 has become one of most investigated genes. Once activated by stress, p53 leads to cellular responses such as cell cycle arrest and apoptosis.

METHODOLOGY/PRINCIPAL FINDINGS

Most previous models have ignored the basal dynamics of p53 under nonstressed conditions. To explore the basal dynamics of p53, we constructed a stochastic delay model by incorporating two negative feedback loops. We found that protein distribution of p53 under nonstressed condition is highly skewed with a fraction of cells showing high p53 levels comparable to those observed under stressed conditions. Under nonstressed conditions, asynchronous and spontaneous p53 pulses are triggered by basal DNA double strand breaks produced during normal cell cycle progression. The first peaking times show a predominant G1 distribution while the second ones are more widely distributed. The spontaneous pulses are triggered by an excitable mechanism. Once initiated, the amplitude and duration of pulses remain unchanged. Furthermore, the spontaneous pulses are filtered by ataxia telangiectasia mutated protein mediated posttranslational modifications and do not result in substantial p21 transcription. If challenged by externally severe DNA damage, cells generate synchronous p53 pulses and induce significantly high levels of p21. The high expression of p21 can also be partially induced by lowering the deacetylation rate.

CONCLUSIONS

Our results demonstrated that the dynamics of p53 under nonstressed conditions is initiated by an excitable mechanism and cells become fully responsive only when cells are confronted with severe damage. These findings advance our understanding of the mechanism of p53 pulses and unlock many opportunities to p53-based therapy.

Noise amplification in human tumor suppression following gamma irradiation

The influence of noise on oscillatory motion is a subject of permanent interest, both for fundamental and practical reasons. Cells respond properly to external stimuli by using noisy systems. We have clarified the effect of intrinsic noise on the dynamics in the human cancer cells following gamma irradiation. It is shown that the large amplification and increasing mutual information with delay are due to coherence resonance. Furthermore, frequency domain analysis is used to study the mechanisms.

Coordination of the nuclear and cytoplasmic activities of p53 in response to DNA damage

The tumor suppressor p53 plays a key role in the cellular response to various stresses. Most previous studies have focused on either the nuclear or cytoplasmic proapoptotic functions of p53, ignoring the combination of both functions. To explore how the two functions of p53 are coordinated in the DNA damage response via computer simulation, we construct a model for the p53 network comprising coupled positive and negative feedback loops involving p53, Mdm2, and Akt, as well as PUMA and Bax. In our model p53 is stabilized and accumulates in the nucleus and cytoplasm upon DNA damage. Nuclear p53 induces expression of Mdm2, PTEN, PUMA, and Bax. Cytoplasmic p53 is then released from the p53.Bcl-xL complex by PUMA to activate Bax directly. We find that the switching between low and high protein levels underlies the decision between cell survival and death. Moreover, a balance between the nuclear and cytoplasmic p53 levels and appropriate levels of Akt and PUMA are required for reliable cell fate decision. Our results indicate that coordination of the transcription-dependent and -independent activities of p53 is important in determining cellular outcomes. These findings advance our understanding of the mechanism for p53-mediated cellular responses and provide clues to p53-based cancer therapy.

Intrinsic noise and Hill dynamics in the p53 system

The p53 feedback loop can induce cellular senescence, cell cycle arrest and apoptosis in response to various stresses, including DNA damage, hypoxia and nutrient deprivation. Using a stochastic model of the negative feedback circuit involving p53 and its inhibitor Mdm2, we present the different oscillatory dynamics at the single-cell and population-cell levels as described in the experiments, and the resonant nature of the oscillations is captured. The stationary distributions of protein populations are characterized by non-Gaussian nature which is due to the interplay between time delay and nonlinearity of reactions.

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.

Expression of recombinant human adenovirus-p53 gene delivered into rabbit VX2 tumors by interventional methods

AIM: To investigate the expression of recombinant human adenovirus-p53 gene (rAd-p53, 1×10¹² VP) delivered into rabbit VX2 tumors by interventional methods.

METHODS: Thirty New Zealand rabbits were used to establish an animal model of hepatocarcinoma by implantation of the VX2 tumors into the liver. Eleven days after implantation, magnetic resonance imaging (MRI) and computed tomography (CT) scans were performed. Successful implantation was achieved in 27 rabbits. These rabbits were randomly divided into three groups: A, B and C. Groups A and B were infused through the hepatic artery with rAd-p53 (5 × 10⁶ VP) alone and rAd-p53 (5 × 10⁶ VP) plus lipiodol ultra fluid (0.5 mL), respectively, while group C was injected with rAd-p53 (5 × 10⁶ VP) into tumors by surgery. All model rabbits underwent digital subtraction angiography (DSA) during operation. The expression of P53, Bax and Bcl-2 proteins in tumor and peritumoral tissue was detected by immunohistochemistry at 24 and 72 hours and on 6 days after the operation.

RESULTS: MRI and CT scans showed that successful implantation was achieved in 27 rabbits, and the success rate was 95%. Immunohistochemical analysis indicated that P53, Bax and Bcl-2 proteins were expressed in all the three groups. The expression levels of P53 and Bax reached the peak at 72 hours. The expression levels of P53 and Bax in group C at 72 hours were significantly higher than those in groups A and B (P53: 70.3% ± 3.4% vs 62.4% ± 3.2% and 65.4% ± 2.1%; Bax: 43.7% ± 2.1% vs 35.7% ± 1.3% and 37.6% ± 2.4%, respectively; all P < 0.05). In contrast, the expression level of Bcl-2 in group C at 72 hours was significantly lower than those in groups A and B (12.0% ± 2.6% vs 18.0% ± 4.3% and 16.2% ± 3.1%, respectively; both P < 0.05).

CONCLUSION: P53 protein is expressed in rabbit VX2 tumors in a time-dependent manner. The expression level of P53 protein depends on the delivery methods used. The highest expression of P53 protein is achieved by intratumoral injection, followed by infusion of rAd-p53 plus lipiodol ultra fluid and rAd-p53 alone through the hepatic artery.