Late stages of mineralization and their signature on the bone mineral density distribution

Experimental measurements of bone mineral density distributions (BMDDs) enable a determination of secondary mineralization kinetics in bone, but the maximum degree of mineralization and how this maximum is approached remain uncertain. We thus test computationally different hypotheses on late stages of bone mineralization by simulating BMDDs in low-turnover conditions.

An established computational model of the BMDD that accounts for mineralization and remodeling processes was extended to limit mineralization to various maximum calcium capacities of bone. Simulated BMDDs obtained by reducing turnover rate from the reference trabecular BMDD under different assumptions on late stage mineralization kinetics were compared with experimental BMDDs of low-turnover bone.

Simulations show that an abrupt stopping of mineralization near a maximum calcium capacity induces a pile-up of minerals in the BMDD statistics that is not observed experimentally. With a smooth decrease of mineralization rate, imposing low maximum calcium capacities helps to match peak location and width of simulated low-turnover BMDDs with peak location and width of experimental BMDDs, but results in a distinctive asymmetric peak shape. No tuning of turnover rate and maximum calcium capacity was able to explain the differences found in experimental BMDDs between trabecular bone (high turnover) and femoral cortical bone (low turnover).

Secondary mineralization in human bone does not stop abruptly, but continues slowly up to a calcium content greater than 30 wt% Ca. The similar mineral heterogeneity seen in trabecular and femoral cortical bones at different peak locations was unexplained by the turnover differences tested.