A MULTISCALE COMPUTATIONAL MODEL FOR THE GROWTH OF THE CRANIAL VAULT IN CRANIOSYNOSTOSIS

A COMPUTATIONAL ANALYSIS OF BONE FORMATION IN THE CRANIAL VAULT USING A COUPLED REACTION-DIFFUSION-STRAIN MODEL

Bones of the murine cranial vault are formed by differentiation of mesenchymal cells into osteoblasts, a process that is primarily understood to be controlled by a cascade of reactions between extracellular molecules and cells. We assume that the process can be modeled using Turing's reaction-diffusion equations, a mathematical model describing the pattern formation controlled by two interacting molecules (activator and inhibitor). In addition to the processes modeled by reaction-diffusion equations, we hypothesize that mechanical stimuli of the cells due to growth of the underlying brain contribute significantly to the process of cell differentiation in cranial vault development. Structural analysis of the surface of the brain was conducted to explore the effects of the mechanical strain on bone formation. We propose a mechanobiological model for the formation of cranial vault bones by coupling the reaction-diffusion model with structural mechanics. The mathematical formulation was solved using the finite volume method. The computational domain and model parameters are determined using a large collection of experimental data that provide precise three dimensional (3D) measures of murine cranial geometry and cranial vault bone formation for specific embryonic time points. The results of this study suggest that mechanical strain contributes information to specific aspects of bone formation. Our mechanobiological model predicts some key features of cranial vault bone formation that were verified by experimental observations including the relative location of ossification centers of individual vault bones, the pattern of cranial vault bone growth over time, and the position of cranial vault sutures.

A computational analysis of bone formation in the cranial vault in the mouse

Bones of the cranial vault are formed by the differentiation of mesenchymal cells into osteoblasts on a surface that surrounds the brain, eventually forming mineralized bone. Signaling pathways causative for cell differentiation include the actions of extracellular proteins driven by information from genes. We assume that the interaction of cells and extracellular molecules, which are associated with cell differentiation, can be modeled using Turing's reaction-diffusion model, a mathematical model for pattern formation controlled by two interacting molecules (activator and inhibitor). In this study, we hypothesize that regions of high concentration of an activator develop into primary centers of ossification, the earliest sites of cranial vault bone. In addition to the Turing model, we use another diffusion equation to model a morphogen (potentially the same as the morphogen associated with formation of ossification centers) associated with bone growth. These mathematical models were solved using the finite volume method. The computational domain and model parameters are determined using a large collection of experimental data showing skull bone formation in mouse at different embryonic days in mice carrying disease causing mutations and their unaffected littermates. The results show that the relative locations of the five ossification centers that form in our model occur at the same position as those identified in experimental data. As bone grows from these ossification centers, sutures form between the bones.