Defective bone mineralization has serious clinical manifestations, including deformities and fractures, but the regulation of this extracellular process is not fully understood. We have developed a mathematical model consisting of ordinary differential equations that describe collagen maturation, production and degradation of inhibitors, and mineral nucleation and growth. We examined the roles of individual processes in generating normal and abnormal mineralization patterns characterized using two outcome measures: mineralization lag time and degree of mineralization. Model parameters describing the formation of hydroxyapatite mineral on the nucleating centers most potently affected the degree of mineralization, while the parameters describing inhibitor homeostasis most effectively changed the mineralization lag time. Of interest, a parameter describing the rate of matrix maturation emerged as being capable of counter-intuitively increasing both the mineralization lag time and the degree of mineralization. We validated the accuracy of model predictions using known diseases of bone mineralization such as osteogenesis imperfecta and X-linked hypophosphatemia. The model successfully describes the highly nonlinear mineralization dynamics, which includes an initial lag phase when osteoid is present but no mineralization is evident, then fast primary mineralization, followed by secondary mineralization characterized by a continuous slow increase in bone mineral content. The developed model can potentially predict the function for a mutated protein based on the histology of pathologic bone samples from mineralization disorders of unknown etiology.

2016-12-23T18:49:38

We present a generalized Keller-Segel model where an arbitrary number of chemical compounds react, some of which are produced by a species, and one of which is a chemoattractant for the species. To investigate the stability of homogeneous stationary states of this generalized model, we consider the eigenvalues of a linearized system. We are able to reduce this infinite dimensional eigenproblem to a parametrized finite dimensional eigenproblem. By matrix theoretic tools, we then provide easily verifiable sufficient conditions for destabilizing the homogeneous stationary states. In particular, one of the sufficient conditions is that the chemotactic feedback is sufficiently strong. Although this mechanism was already known to exist in the original Keller-Segel model, here we show that it is more generally applicable by significantly enlarging the class of models exhibiting this instability phenomenon which may lead to pattern formation.