Appearance and location of secondary ossification centres may be explained by a reaction–diffusion mechanism

Computational model of articular cartilage regeneration induced by scaffold implantation in vivo

Computational models allow to explain phenomena that cannot be observed through an animal model, such as the strain and stress states which can highly influence regeneration of the tissue. For this purpose, we have developed a simulation tool to determine the mechanical conditions provided by the polymeric scaffold. The computational model considered the articular cartilage, the subchondral bone, and the scaffold. All materials were modeled as poroelastic, and the cartilage had linear-elastic oriented collagen fibers. This model was able to explain the remodeling process that subchondral bone goes through, and how the scaffold allowed the conditions for cartilage regeneration. These results suggest that the use of scaffolds might lead the cartilaginous tissue growth in vivo by providing a better mechanical environment. Moreover, the developed computational model demonstrated to be useful as a tool prior experimental in vivo studies, by predicting the possible outcome of newly proposed treatments allowing to discard approaches that might not bring good results.

Computational model of endochondral ossification: Simulating growth of a long bone

Mechanical loading is a crucial factor in joint and bone development. Using a computational model, we investigated the role of mechanics on cartilage growth rate, ossification of the secondary center, formation of the growth plate, and overall bone shape. A computational algorithm was developed and implemented into finite element models to simulate the endochondral ossification for symmetric and asymmetric motion in a generic diarthrodial joint. Under asymmetric loading condition the secondary center ossifies asymmetrically leaning toward the external load and results in tilted growth plate. Also the mechanics seems to have greater influence in the early onset of the ossification of the secondary center rather than later progression of the center. While previous models have simulated select stages of skeletal development, our model can simulate growth and ossification during the entirety of post-natal development. Such computational models of skeletal development may provide insight into specific loading conditions that cause bone and joint deformities, and the required timing for rehabilitative repair.

Cellular scale model of growth plate: An in silico model of chondrocyte hypertrophy

The growth plate is the responsible for longitudinal bone growth. It is a cartilaginous structure formed by chondrocytes that are continuously undergoing a differentiation process that starts with a highly proliferative state, followed by cellular hypertrophy, and finally tissue ossification. Within the growth plate chondrocytes display a characteristic columnar organization that potentiates longitudinal growth. Both chondrocyte organization and hypertrophy are highly regulated processes influenced by biochemical and mechanical stimuli. These processes have been studied mainly using in vivo models, although there are few computational approaches focused on the rate of ossification rather than events at cellular level. Here, we developed a model of cellular behavior integrating biochemical and structural factors in a single column of cells in the growth plate. In our model proliferation and hypertrophy were controlled by biochemical regulatory loop formed between Ihh and PTHrP (modeled as a set of reaction-diffusion equations), while cell growth was controlled by mechanical loading. We also examined the effects of static loading. The model reproduced the proliferation and hypertrophy of chondrocytes in organized columns. This model constitutes a first step towards the development of mechanobiological models that can be used to study biochemical interactions during endochondral ossification.

Evolution and functional significance of derived sternal ossification patterns in ornithothoracine birds

The midline pattern of sternal ossification characteristic of the Cretaceous enantiornithine birds is unique among the Ornithodira, the group containing birds, nonavian dinosaurs and pterosaurs. This has been suggested to indicate that Enantiornithes is not the sister group of Ornithuromorpha, the clade that includes living birds and their close relatives, which would imply rampant convergence in many nonsternal features between enantiornithines and ornithuromorphs. However, detailed comparisons reveal greater similarity between neornithine (i.e. crown group bird) and enantiornithine modes of sternal ossification than previously recognized. Furthermore, a new subadult enantiornithine specimen demonstrates that sternal ossification followed a more typically ornithodiran pattern in basal members of the clade. This new specimen, referable to the Pengornithidae, indicates that the unique ossification pattern observed in other juvenile enantiornithines is derived within Enantiornithes. A similar but clearly distinct pattern appears to have evolved in parallel in the ornithuromorph lineage. The atypical mode of sternal ossification in some derived enantiornithines should be regarded as an autapomorphic condition rather than an indication that enantiornithines are not close relatives of ornithuromorphs. Based on what is known about molecular mechanisms for morphogenesis and the possible selective advantages, the parallel shifts to midline ossification that took place in derived enantiornithines and living neognathous birds appear to have been related to the development of a large ventral keel, which is only present in ornithuromorphs and enantiornithines. Midline ossification can serve to medially reinforce the sternum at a relatively early ontogenetic stage, which would have been especially beneficial during the protracted development of the superprecocial Cretaceous enantiornithines.

Growth plate stress distribution implications during bone development: a simple framework computational approach

Mechanical stimuli play a significant role in the process of long bone development as evidenced by clinical observations and in vivo studies. Up to now approaches to understand stimuli characteristics have been limited to the first stages of epiphyseal development. Furthermore, growth plate mechanical behavior has not been widely studied. In order to better understand mechanical influences on bone growth, we used Carter and Wong biomechanical approximation to analyze growth plate mechanical behavior, and explore stress patterns for different morphological stages of the growth plate. To the best of our knowledge this work is the first attempt to study stress distribution on growth plate during different possible stages of bone development, from gestation to adolescence. Stress distribution analysis on the epiphysis and growth plate was performed using axisymmetric (3D) finite element analysis in a simplified generic epiphyseal geometry using a linear elastic model as the first approximation. We took into account different growth plate locations, morphologies and widths, as well as different epiphyseal developmental stages. We found stress distribution during bone development established osteogenic index patterns that seem to influence locally epiphyseal structures growth and coincide with growth plate histological arrangement.

Growth of the flat bones of the membranous neurocranium: a computational model

This article assumes two stages in the formation of the bones in the calvaria, the first one takes into account the formation of the primary centers of ossification. This step counts on the differentiation from mesenchymal cells into osteoblasts. A molecular mechanism is used based on a system of reaction-diffusion between two antagonistic molecules, which are BMP2 and Noggin. To this effect we used equations whose behavior allows finding Turing patterns that determine the location of the primary centers. In the second step of the model we used a molecule that is expressed by osteoblasts, called Dxl5 and that is expressed from the osteoblasts of each flat bone. This molecule allows bone growth through its borders through cell differentiation adjacent to each bone of the skull. The model has been implemented numerically using the finite element method. The results allow us to observe a good approximation of the formation of flat bones of the membranous skull as well as the formation of fontanelles and sutures.

Spongiosa primary development: a biochemical hypothesis by Turing patterns formations

We propose a biochemical model describing the formation of primary spongiosa architecture through a bioregulatory model by metalloproteinase 13 (MMP13) and vascular endothelial growth factor (VEGF). It is assumed that MMP13 regulates cartilage degradation and the VEGF allows vascularization and advances in the ossification front through the presence of osteoblasts. The coupling of this set of molecules is represented by reaction-diffusion equations with parameters in the Turing space, creating a stable spatiotemporal pattern that leads to the formation of the trabeculae present in the spongy tissue. Experimental evidence has shown that the MMP13 regulates VEGF formation, and it is assumed that VEGF negatively regulates MMP13 formation. Thus, the patterns obtained by ossification may represent the primary spongiosa formation during endochondral ossification. Moreover, for the numerical solution, we used the finite element method with the Newton-Raphson method to approximate partial differential nonlinear equations. Ossification patterns obtained may represent the primary spongiosa formation during endochondral ossification.

A mathematical model of the process of ligament repair: effect of cold therapy and mechanical stress

This article proposes a mathematical model that predicts the wound healing process of the ligament after a sprain, grade II. The model describes the swelling, expression of the platelet-derived growth factor (PDGF), formation and migration of fibroblasts into the injury area and the expression of collagen fibers. Additionally, the model can predict the effect of ice treatment in reducing inflammation and the action of mechanical stress in the process of remodeling of collagen fibers. The results obtained from computer simulation show a high concordance with the clinical data previously reported by other authors.

Relating the chondrocyte gene network to growth plate morphology: from genes to phenotype

During endochondral ossification, chondrocyte growth and differentiation is controlled by many local signalling pathways. Due to crosstalks and feedback mechanisms, these interwoven pathways display a network like structure. In this study, a large-scale literature based logical model of the growth plate network was developed. The network is able to capture the different states (resting, proliferating and hypertrophic) that chondrocytes go through as they progress within the growth plate. In a first corroboration step, the effect of mutations in various signalling pathways of the growth plate network was investigated.

A MATHEMATICAL MODEL OF THE GROWTH PLATE

The growth plate is a structure formed of cells called chondrocytes; these are arranged in columns and provide the elongation of bone due to their proliferation and hypertrophy. In each column, we can see chondrocytes in their proliferating state, which are constantly dividing, and in hypertrophic state, which grow in a nearly spherical shape. These cells express different proteins and molecules throughout their half-life and exhibit a special behavior depending on their local mechanical and biochemical environments. This article develops a mathematical model that describes the relationship of geometry, growth by proliferation and hypertrophy, and vascular invasion with biochemical and mechanical factors present during endochondral ossification.

A model of cerebral cortex formation during fetal development using reaction-diffusion-convection equations with Turing space parameters

The cerebral cortex is a gray lamina formed by bodies of neurons covering the cerebral hemispheres, varying in thickness from 1.25 mm in the occipital lobe to 4mm in the anterior lobe. The brain's surface is about 30 times greater that of the skull because of its many folds; such folds form the gyri, sulci and fissures and mark out areas having specific functions, divided into five lobes. Convolution formation may vary between individuals and is an important feature of brain formation; such patterns can be mathematically represented as Turing patterns. This article describes how a phenomenological model was developed by describing the formation pattern for the gyri occurring in the cerebral cortex by reaction diffusion equations with Turing space parameters. Numerical examples for simplified geometries of a brain were solved to study pattern formation. The finite element method was used for the numerical solution, in conjunction with the Newton-Raphson method. The numerical examples showed that the model can represent cerebral cortex fold formation and reproduce pathologies related to gyri formation, such as polymicrogyria and lissencephaly.

SELF-ASSEMBLED SCAFFOLDS USING REACTION–DIFFUSION SYSTEMS: A HYPOTHESIS FOR BONE REGENERATION

One area of tissue engineering concerns research into alternatives for new bone formation and replacing its function. Scaffolds have been developed to meet this requirement, allowing cell migration, bone tissue growth, transport of growth factors and nutrients, and the improvement of the mechanical properties of bone. Scaffolds are made from different biomaterials and manufactured using several techniques that, in some cases, do not allow full control over the size and orientation of the pores characterizing the scaffold. A novel hypothesis that a reaction–diffusion (RD) system can be used for designing the geometrical specifications of the bone matrix is thus presented here. The hypothesis was evaluated by making simulations in two- and three-dimensional RD systems in conjunction with the biomaterial scaffold. The results showed the methodology's effectiveness in controlling features such as the percentage of porosity, size, orientation, and interconnectivity of pores in an injectable bone matrix produced by the proposed hypothesis.

A mechanobiological model of epiphysis structures formation

Developing bone consists of epiphysis, metaphysis and diaphysis. The secondary ossification centre (SOC) appears and grows within the epiphysis, involving two histological stages. Firstly, cartilage canals appear; they carry hypertrophy factors towards the central area of the epiphysis. Canal growth and expansion is modulated by stress on the epiphysis. Secondly, the diffusion of hypertrophy factors causes SOC growth. Hypertrophy is regulated by biological and mechanical factors present within the epiphysis. The finite element method has been used for solving a coupled system of differential equations for modelling these histological stages of epiphyseal development. Cartilage canal spatial-temporal growth patterns were obtained as well as the SOC formation pattern. This model qualitatively agreed with experimental results reported by other authors.

A mathematical model for describing the metastasis of cancer in bone tissue

Metastasis is the rapid proliferation of cancer cells (secondary tumour) at a specific place, generally leading to death. This occurs at anatomical parts providing the necessary environment for vascularity, oxygen and food to hide their actions and trigger the rapid growth of cancer. Prostate and breast cancers, for example, use bone marrow for their proliferation. Bone-supporting cancer cells thus adapt to the environment, mimicking the behaviour of genetic and molecular bone cells. Evidence of this has been given in Cecchini et al. (2005, EAU Update Ser. 3:214-226), providing arguments such as how cancer cell growth is so active during bone reabsorption. This paper simulates metastasis activation in bone marrow. A mathematical model has been developed involving the activation of molecules from bone tissue cells, which are necessary for cancer to proliferate. Here, we simulate two forms of secondary tumour growth depending on the type of metastasis: osteosclerosis and osteolysis.

A mathematical model of medial collateral ligament repair: migration, fibroblast proliferation and collagen formation

The partial rupture of ligament fibres leads to an injury known as grade 2 sprain. Wound healing after injury consists of four general stages: swelling, release of platelet-derived growth factor (PDGF), fibroblast migration and proliferation and collagen production. The aim of this paper is to present a mathematical model based on reaction-diffusion equations for describing the repair of the medial collateral ligament when it has suffered a grade 2 sprain. We have used the finite element method to solve the equations of this. The results have simulated the tissue swelling at the time of injury, predicted PDGF influence, the concentration of fibroblasts migrating towards the place of injury and reproduced the random orientation of immature collagen fibres. These results agree with experimental data reported by other authors. The model describes wound healing during the 9 days following such injury.

A mathematical model of epiphyseal development: hypothesis on the cartilage canals growth

The role of cartilage canals is to transport nutrients and biological factors that cause the appearance of the secondary ossification centre (SOC). The SOC appears in the centre of the epiphysis of long bones. The canal development is a complex interaction between mechanical and biological factors that guide its expansion into the centre of the epiphysis. This article introduces the 'Hypothesis on the growth of cartilage canals'. Here, we have considered that the development of these canals is an essential event for the appearance of SOC. Moreover, it is also considered to be important for the transport of molecular factors (RUNX2 and MMP9) at the ends of such canals. Once the canals are merged in the centre of the epiphysis, these factors are released causing hypertrophy of adjacent cells. This RUNX2 and MMP9 release occurs due to the action of mechanical loads that supports the epiphysis. In order to test this hypothesis, we use a hybrid approach using the finite element method to simulate the mechanical stresses present in the epiphysis and the cellular automata to simulate the expansion of the canals and the hypertrophy factors pathway. By using this hybrid approach, we have obtained as a result the spatial-temporal patterns for the growth of cartilage canals and hypertrophy factors within the epiphysis. The model is in qualitative agreement with experimental results previously reported by other authors. Thus, we conclude that this model may be used as a methodological basis to present a complete mathematical model of the processes involved in epiphyseal development.