The bone marrow is necessary for renewal of all hematopoietic cells and
critical for maintenance of a wide range of physiologic functions. Multiple
human diseases result from bone marrow dysfunction. It is also the site in which
“liquid” tumors, including leukemia and multiple myeloma,
develop as well as a frequent site of metastases. Understanding the complex
cellular and microenvironmental interactions that govern normal bone marrow
function as well as diseases and cancers of the bone marrow would be a valuable
medical advance. Our goal is the development of a spatially-explicit in
silico model of the bone marrow to understand both its normal
function and the evolutionary dynamics that govern the emergence of bone marrow
malignancy.
Here we introduce a multiscale computational model of the bone marrow
that incorporates three distinct spatial scales, cell, hematopoietic subunit,
whole marrow. Implemented as a fixed lattice 3D cellular automaton, it
reproduces the spatial characteristics of the normal bone marrow and is
validated against data from the daily production of mature blood cells and
response of hematopoiesis after irradiation. The major mechanisms modeled in
this work are: (1) replication, specialization and migration of hematopoietic
cells, (2) optimized spatial configuration of sinuses and hematopoietic
compartments and, (3) intravasation of mature hematopoietic cells into
sinuses.
Our results, using parameter estimates from literature, recapitulates
normal bone marrow function and suggest an explanation for the fractal-like
structure of trabeculae and sinuses in the marrow, which would be an
optimization of the hematopoietic function in order to maximize the number of
mature blood cells produced daily within the volumetric restrictions of the
marrow.