Virtual Cell Virtual Cell WebsiteSchool of Health Sciences WebsitePurdue University Website
About the Virtual CellDownload VCExamplesTheoryUser ManualHelpExecute VC Simulation
Multiscale modeling in radiation biology
Author: Robert Stewart, Ph.D.
Theory
Cell-to-cell interactions and the tissue microenvironment

Low-dose experiments with light ions and microbeam and medium-transfer experiments demonstrate that radiation-damaged cells can induce genomic instability, mutations, and enhance cell killing and neoplastic transformation in nearby undamaged cells (reviewed in Ballarini et al. 2002, Smith et al. 2003, and Morgan 2003a, 2003b). Communication through gap junctions and through secreted (diffusible) messengers both contribute to the observed detrimental bystander effects. Radiation damage to the extracellular matrix (ECM) is another possible mode of action. For example, reactive oxygen species (ROSs) produced in the ECM by radiation may convert latent transforming growth factor β (TGF-β to an active complex that fosters neoplastic growth (Park et al. 2000, Barcellos-Hoff 2001). Persistent disruption of the tissue microenvironment by radiation or other means (e.g., disease, age, and wounds) may give unstable cells a selective growth advantage over normal cells and contribute to tumor formation (Barcellos-Hoff 2001, Barcellos-Hoff and Brooks 2001, Nunney 1999, Radisky 2001, Robert 2000, 2002).

On the other hand, many converging lines of evidence also indicate that multi-cellular communities (i.e., tissues) take an active role in suppressing neoplastic growth.  One of the ways in which a tissue defends against neoplastic growth is through the selective ablation of transformed cells through apoptosis (Bauer 1996, 2002, Dormann et al. 1998, Häufel et al. 1999).  Normal cells have also been shown to inhibit aberrant-cell growth by limiting proliferation (Bignami 1988a, 1988b, Stoker 1964, 1967, Stoker et al. 1966) and by inducing the terminal differentiation of transformed cells (La Rocca et al. 1989).  The effective inhibition of neoplastic growth for all of these modes of action depends critically on the density and the spatial arrangement (size) of the normal and transformed cell populations (reviewed in Bauer 1996).

The evolution of multi-cellular organisms required mechanisms to coordinate cell function and behavior, and the observed positive and negative effects that arise from cell-to-cell interactions may be due to these same homeostatic controls (Barcellos-Hoff and Brooks 2001).  Regardless, cell signaling and tissue microenvironmental factors cast doubt on cell-autonomous modes of action, a central tenet of the current radiobiological paradigm, and pose many new challenges for radiation-response modeling at the cell and tissue levels.

Multiscale Modeling

The physical, chemical and biological processes that ultimately lead to cancer span at least ten orders of magnitude on the time scale (< 10-3 s to > 107 s) and ten orders of magnitude on the spatial scale (~ 10-10 m to 1 m).

Multi-stage carcinogenesis models, such as the two-mutation clonal expansion model (Moolgavkar et al. 1988, Moolgavkar and Luebeck 1990), capture many of the putative key steps associated with the development of cancer, and they do so using as few adjustable parameters as possible. However, the desire to keep the number of adjustable model inputs to a minimum means that the effects of many molecular and cellular processes must be lumped together. That is, the details of the specific biological mechanisms of action involved in the development of cancer must be replaced by nonspecific (average) mechanisms that capture events and processes occurring over very large time and spatial scales. However, some in vitro experiments suggest that the risks of low dose radiation exposures may be larger (Sawant et al. 2001, Mothersill et al. 2002) or smaller (Azzam et al. 1996, Redpath et al. 2001, 2003) than the risks expected using simple linear high-to-low-dose extrapolation procedures. Studies such as these indicate that we need to find practical strategies to introduce additional (more specific) biological details into the cancer modeling process.

Multi-scale modeling provides one avenue to better integrate cell signaling and microenvironmental factors into radiation response models. The multi-scale modeling approach begins by formulating relatively detailed mathematical models for specific molecular and cellular processes. For example, one model might capture key elements of the base excision repair (BER) pathway while another model relates the overall mutation rate to the onset of genomic instability and to the eventual neoplastic transformation of a cell or one of its progeny. As a first test of the postulated mechanisms of action, these lower-level models (“sub-models”) are tested against measured data from in vitro systems. To model the emergent behavior (response) of a group of cells or a tissue, these sub-models are linked together to form a “supermodel.” The links (information pathways) among the sub-models represent specific hypotheses about the key events and processes involved in the development of cancer. As a second test of the postulated mechanisms, results from the multi-scale supermodel can be compared to data from in vivo systems.

An attractive feature of the multi-scale modeling approach is that it provides a systematic method to balance the need for greater biological detail against the need for simplicity (fewer adjustable parameters).  That is, each sub-model uses a minimum number of adjustable parameters to capture the putative mechanisms of action underlying a particular endpoint of interest (e.g., mutagenesis, genomic instability or apoptosis).  The potential usefulness of multi-scale biological models is driven by the hypothesis that mechanisms of action operate in approximately the same manner in vitro and in vivo.  If mechanisms of action operate in approximately the same manner, in vitro data can be used to identify appropriate mathematical models for specific mechanisms of action.  The in vitro data may also be used help define bounds on the range of parameter values (model inputs) appropriate for in vivo systems.

The hypothesis that mechanisms and parameters are approximately the same in vitro and in vivo is reasonably well-justified for endpoints such as DNA damage and, perhaps to a lesser extent, the outcome from a pathway-specific DNA repair event (probability a lesion is correctly or incorrectly repaired).  As we progress up through the biological hierarchy to processes such as the induction of apoptosis or the regulation of cell growth or differentiation, the use of model inputs derived from in vitro data becomes increasingly questionable because of differences in the in vitro and in vivo microenvironment.  However, setting some model inputs a priori to an approximate in vitro value may still be preferable to treating all model inputs as purely adjustable.  If all model inputs are treated as adjustable, parameter identifiability issues arise (Hanin and Yakovlev 1996, Hanin and Boucher 1999).  Parameter estimation issues are a, if not the, major challenge associated with the effort to construct mechanism-based radiation response models.


References

  • E.I. Azzam, S.M. de Toledo, and G.P. Raaphorst, Mitchel R.E. Low-dose ionizing radiation decreases the frequency of neoplastic transformation to a level below the spontaneous rate in C3H 10T1/2 cells. Radiat. Res. 146(4), 369-73 (1996).
  • F. Ballarini, M. Biaggi, A. Ottolenghi and O. Sapora. Cellular communication and bystander effects: a critical review for modelling low-dose radiation action. Mutat. Res. 501(1-2), 1-12 (2002).
  • M.H. Barcellos-Hoff It takes a tissue to make a tumor: epigenetics, cancer and the microenvironment. J. Mammary Gland Biol Neoplasia 6(2), 213-221 (2001).
  • M.H. Barcellos-Hoff and A. L. Brooks.  Extracellular signaling through the microenvironment: a hypothesis relating carcinogenesis, bystander effects, and genomic instability.  Radiat. Res. 156, 618-627 (2001).
  • G. Bauer. Elimination of transformed cells by normal cells: a novel concept for the control of carcinogenesis. Histol. Histopathol. 11(1), 237-255 (1996).
  • G. Bauer. Signaling and proapoptotic functions of transformed cell-derived reactive oxygen species. Prostaglandins Leukot. Essent. Fatty Acids. 66(1), 41-56 (2002).
  • M. Bignami, S. Rosa, S.A. La Rocca, G. Falcone and F. Tato. Differential influence of adjacent normal cells on the proliferation of mammalian cells transformed by the viral oncogenes myc, ras and src. Oncogene 2(5), 509-514 (1988a).
  • M. Bignami, S. Rosa, G. Falcone, F. Tato, F. Katoh, H. Yamasaki. Specific viral oncogenes cause differential effects on cell-to-cell communication, relevant to the suppression of the transformed phenotype by normal cells. Mol. Carcinog. 1(1), 67-75 (1988b).
  • S. Dormann and G. Bauer. TGF-beta and FGF-trigger intercellular induction of apoptosis: analogous activity on non-transformed but differential activity on transformed cells. Int. J. Oncol. 13(6), 1247-1252 (1998).
  • L.G. Hanin, A.Y. Yakovlev. A nonidentifiability aspect of the two-stage model of carcinogenesis. Risk Anal. 16(5), 711-715 (1996).
  • L.G. Hanin, K.M. Boucher. Identifiability of parameters in the Yakovlev-Polig model of carcinogenesis. Math Biosci. 160(1), 1-24 (1999).
  • T. Häufel, S. Dormann, J. Hanusch, A. Schwieger and G. Bauer. Three distinct roles for TGF-beta during intercellular induction of apoptosis: a review. Anticancer Res. 19(1A), 105-111 (1999).
  • La Rocca S.A., M. Grossi, G. Falcone, S. Alema and F. Tato. Interaction with normal cells suppresses the transformed phenotype of v-myc-transformed quail muscle cells. Cell 58(1), 123-31 (1989).
  • S.H. Moolgavkar, A. Dewanji and D.J. Venzon. A stochastic two-stage model for cancer risk assessment. I. The hazard function and the probability of tumor. Risk Anal. 8(3), 383-392 (1988).
  • S.H. Moolgavkar and G. Luebeck. Two-event model for carcinogenesis: biological, mathematical, and statistical considerations. Risk Anal. 10(2), 323-341 (1990).
  • W.F. Morgan. Non-targeted and delayed effects of exposure to ionizing radiation: I. Radiation-induced genomic instability and bystander effects In Vitro. Radiat. Res. 159(5), 567-580 (2003a).
  • W.F. Morgan. Non-targeted and delayed effects of exposure to ionizing radiation: II. Radiation-induced genomic instability and bystander effects in vivo, clastogenic factors and transgenerational effects. Radiat. Res. 159(5), 581-596 (2003b).
  • C. Mothersill, C. B. Seymour and M.C. Joiner. Relationship between radiation-induced low-dose hypersensitivity and the bystander effect. Radiat. Res. 157(5), 526-32 (2002).
  • L. Nunney. Lineage selection and the evolution of multistage carcinogenesis. Proc. R. Soc. Lond. B Biol. Sci. 266(1418), 493-498 (1999).
  • C.C. Park, M.J. Bissell and M.H. Barcellos-Hoff. The influence of the microenvironment on the malignant phenotype. Mol. Med. Today 6(8), 324-329 (2000).
  • D. Radisky, C. Hagios and M.J. Bissell. Tumors are unique organs defined by abnormal signaling and context. Semin. Cancer Biol. 11(2), 87-95 (2001).
  • J.L. Redpath, D. Liang, T.H. Taylor, C.  Christie and E. Elmore. The shape of the dose-response curve for radiation-induced neoplastic transformation in vitro: evidence for an adaptive response against neoplastic transformation at low doses of low-LET radiation. Radiat Res. 156(6), 700-707 (2001).
  • J.L. Redpath, S.C. Short, M. Woodcock and P.J. Johnston. Low-dose reduction in transformation frequency compared to unirradiated controls: the role of hyper-radiosensitivity to cell death. Radiat Res. 159(3), 433-6 (2003).
  • L. Robert. Cellular and molecular mechanisms of aging and age related diseases. Pathol. Oncol. Res. 6(1), 3-9 (2000).
  • L. Robert.  Cell-matrix interactions in cancer spreading--effect of aging. Semin. Cancer Biol. 12(3), 157-63 (2002).
  • S.G. Sawant, G. Randers-Pehrson, C. R. Geard, D.J. Brenner and E.J. Hall. The bystander effect in radiation oncogenesis: I. Transformation in C3H 10T1/2 cells in vitro can be initiated in the unirradiated neighbors of irradiated cells. Radiat. Res. 155(3), 397-401 (2001).
  • L.E. Smith, S. Nagar, G.J. Kim and W.F. Morgan. Radiation-induced genomic instability: radiation quality and dose response. Health Phys. 85(1), 23-29 (2003).
  • M.G.P. Stoker. Regulation of growth and orientation in hamster cells transformed by polyoma virus. Virology 24, 165-174 (1964).
  • M.G.P. Stoker. Transfer of growth inhibition between normal and virus-transformed cells: autoradiographic studies using marked cells. J. Cell Sci. 2, 293-304 (1967).
  • M.G.P. Stoker, M. Shearer M and C. O’ Neil. Growth inhibition of polyoma-transformed cells by contact with static normal fibroblasts. J. Cell Sci. 1, 297-310 (1966).

School of Health Sciences
Purdue University
Disclaimer

Last updated: 10 June, 2011
Webmaster