A mathematical model for radiation-induced myelopoiesis.

ABSTRACT

A model for damage, repair, killing, and repopulation of myelopoietic marrow is presented. Evaluation produces time and dose-rate profiles during and following any complex irradiation. Equations model variable dose rates, multiple exposures, different sources, and arbitrary intervals between treatments. If factors which dominate the control of biological processes can be demonstrated, an option is to set biological rate constants to experimentally determined values. Previously, knowledge did not permit identification of dominating biological processes and their temporal rates. But a unique feature of this study is that unspecified lesions for killing and injury of cells are evaluated from mortality data on the animal species of choice. "Unspecified" is used to indicate a condition of assumption-free modeling of molecular processes, whereby rate constants for cellular effects are simply computed directly from animal mortality data. Coefficients (estimated by maximum-likelihood methods for nonspecific processes) are compared with experimental values for specific processes. The model has many uses, including modeling of the myelopoietic potential as a function of time. Another option is to calculate the whole-body survival curve for cells that control myelopoiesis as a result of the treatment schedule. Also through simple extensions of the model, an extremely complex protocol can be identified with an equivalent prompt dose value--even for partial-body, fractionated exposures.