Human cortical bone: the SiNuPrOs model

Many models that have been developed for cortical bone oversimplify much of the architectural and physical complexity. With SiNuPrOs model, a more complete approach is investigated: it is multiscale because it contains five structural levels and multi physic because it takes into account simultaneously structure (with various properties: elasticity, piezoelectricity, porous medium), fluid and mineralization process modelization. The multiscale aspect is modeled by using 18 structural parameters in a specific application of the mathematical theory of homogenization and 10 other physical parameters are necessary for the multi physic aspect. The modelization of collagen as a piezoelectric medium has needed the development of a new behaviour law allowing a better simulation of the effect of a medium considered as evolving during a mineralization process. Then the main interest of SiNuPrOs deals with the possibility to study, at each level of the cortical architecture, either the elastic properties or the fluid motion or the piezoelectric effects or both of them. All these possibilities constitute a very large work and all this mass of information (fluid aspects, even at the nanoscopic scale, piezoelectric phenomena and simulations) will be presented in several papers. This first one is only devoted to the presentation of this model with an application to the computation of elastic properties at the macroscopic scale. The computational methods have been packed into software also called SiNuPrOs and allowing a large number of predictive simulations corresponding to various different configurations.     

Mechanisms of seizure propagation in a cortical model

We consider a mathematical model of mesoscopic human cortical ictal electrical activity. We compare the model results with ictal electrocortical data recorded from three human subjects and show how the two agree. We determine that, in the model system, seizures result from increased connectivity between excitatory and inhibitory cell populations, or from decreased connectivity within either excitatory or inhibitory cell populations. We compare the model results with the disinhibition and 4-AP models of epilepsy and suggest how the model may guide the development of new anticonvulsant therapies.     

On the mechanical characterization of compact bone structure using the homogenization theory

In a previous paper (Crolet et al., 1993, J. Biomechanics 26, 677–687), a modelling of the mechanical behavior of compact bone was presented, in which the homogenization theory was the basic tool of computation. In this simulation, approximations were used for the modelling of the lamellae and the osteons: the lamella and the osteon were divided into cylindrical sectors, each sector being approximated as a parallelepiped having a periodic structure (fibrous composite for the lamella, superimposition of plates for the osteon). The present study deals with a new model without these approximations. First, it can be proved that the homogenized elasticity tensor for a lamella, which has a non-periodic structure, is obtained at each geometrical point as a homogenized tensor of a periodic problem. Similarly, for the osteonal structure, the components of the homogenized tensor are determined at each point as the result of a periodic homogenization.

The software OSTEON, which is the computational method associated with this model, allows one to obtain a better understanding of the effects of many bony parameters. The obtained results are in accordance with experimental data.     

Human bone cell cultures: a new model for studying the mechanism of action of calcitonin

An investigation on cell cultures obtained from temporal human bone fragments showed that they provide a suitable model for studying the mechanism involved in calcitonin action on bone cells. Furthermore they demonstrated: a transitory increase in 45Ca uptake that returned to control values ten minutes after the hormone was added; a relation between 45Ca uptake and increased cAMP concentrations when these were measured at the same time intervals; a reproduction of the salmon calcitonin (sCT) effect after incubation of the cultures with either db-cAMP or db-cGMP and inhibition of 45Ca uptake and parallel decrease in cAMP levels with propanol. These results suggest that in human bone cell cultures, sCT acts as a temporary promoter of 45Ca uptake, probably by activating an adenylate-cyclase system through a beta-receptor.     

The effect of restriction of dietary calcium on trabecular and cortical bone mineral density in the rats

This study aimed to investigate effects of restricted calcium intake on cortical and trabecular bone density in white rats. Low Ca diet was fed for six weeks, and bone density and bone metabolism parameters were assessed in blood. This study was carried out on 12 male white rats aged 12 weeks (Sprague-Dawley; SD). These rats were bred for 1 week and randomly assigned to the standard calcium diet group (SCa group, n = 6) and the low calcium diet group (LCa group; n = 6). The SCa group was given a modified AIN-93M mineral mix (with 0.5% Ca), which was made by adding calcium to a standard AIN93 diet, and the LCa Group was fed a modified AIN-93 Mineral mix (with 0.1% Ca). Femoral BMD and BMC were measured by DEXA in each rat. After trabecular bone was separated from cortical bone, volumetric bone mineral density (vBMD) was measured using pQCT. Serum Ca and P levels were measured as parameters of bone metabolism, and S-ALP, S-TrACP and-Dpd levels were also measured. The results revealed no significant differences in weight, growth rate, feed consumption and feed efficiency between the two groups before and after calcium-restricted diet (p > .05). No significant differences were also observed in bone length and bone mass between the two groups (p > .05). Although bilateral femoral BMDs were not significantly different between the two groups, bilateral femoral BMCs significantly decreased in the LCa group, compared with the SCa group (p = .023, p = .047). Bilateral cortical MDs were not significantly different between the two groups, either. However, trabecular BMD significantly decreased in the LCa group, compared with the SCa group (p = .041). U-Dpd and S-TrACP levels significantly declined in the LCa group, compared to the SCa group (p = .039, p = .010). There were no significant differences in serum Ca and P levels between the two groups (p > .05). However, a significant decrease in urinary Ca level (p = .001) and a significant increase in urinary P (p = .001) were observed in the LCa group, compared to the Sca group. These findings described that six-week low calcium diet led to decreased trabecular bone density, reduced urinary excretion of Ca and increased urinary excretion of P. As a result, Ca hemeostasis can be maintained.     

Collagen's role in the cortical bone's behaviour: a numerical approach

Literature devoted to experimental measurements of the elastic properties of the human cortical bone gives us a relatively wide spectrum of values. This proves that the result depends on the bone itself and on the area from where the sample has been taken. Ultrasonic measurement, which is a fine technology, points out complex maps. The reason of this very strong heterogeneity is not completely explained. The present study is based on a numerical model of the human cortical bone, the SiNuPrOs model and aims to suggest an explanation. If one admits that mineral apposition occurs around collagen fibres, the spatial orientation of these fibres would have an important consequence on the elastic properties of the medium. On the basis of homogenisation theory allowing to compute all the components of the elasticity tensor, this study quantifies the main influence of this architectural orientation and its effect on the anisotropy of the cortical bone.     

Micromechanical modeling of elastic properties of cortical bone accounting for anisotropy of dense tissue

The paper analyzes the connection between microstructure of the osteonal cortical bone and its overall elastic properties. The existing models either neglect anisotropy of the dense tissue or simplify cortical bone microstructure (accounting for Haversian canals only). These simplifications (related mostly to insufficient mathematical apparatus) complicate quantitative analysis of the effect of microstructural changes – produced by age, microgravity, or some diseases – on the overall mechanical performance of cortical bone. The present analysis fills this gap; it accounts for anisotropy of the dense tissue and uses realistic model of the porous microstructure. The approach is based on recent results of Sevostianov et al. (2005) and Saadat et al. (2012) on inhomogeneities in a transversely-isotropic material. Bone?s microstructure is modeled according to books of Martin and Burr (1989), Currey (2002), and Fung (1993) and includes four main families of pores. The calculated elastic constants for porous cortical bone are in agreement with available experimental data. The influence of each of the pore types on the overall moduli is examined.     

On the role of bone damage in calcium homeostasis

Bone serves as the reservoir of some minerals including calcium. If calcium is needed anywhere in the body, it can be removed from the bone matrix by resorption and put back into the blood flow. During bone remodelling the resorbed tissue is replaced by osteoid which gets mineralized very slowly. Then, calcium homeostasis is controlled by bone remodelling, among other processes: the more intense is the remodelling activity, the lower is the mineral content of bone matrix. Bone remodelling is initiated by the presence of microstructural damage. Some experimental evidences show that the fatigue properties of bone are degraded and more microdamage is accumulated due to the external load as the mineral content increases. That damage initiates bone remodelling and the mineral content is so reduced. Therefore, this process prevents the mineral content of bone matrix to reach very high (non-physiological) values. A bone remodelling model has been used to simulate this regulatory process. In this model, damage is an initiation factor for bone remodelling and is estimated through a fatigue algorithm, depending on the macroscopic strain level. Mineral content depends on bone remodelling and mineralization rate. Finally, the bone fatigue properties are defined as dependent on the mineral content, closing the interconnection between damage and mineral content. The remodelling model was applied to a simplified example consisting of a bar under tension with an initially heterogeneous mineral distribution. Considering the fatigue properties as dependent on the mineral content, the mineral distribution tends to be homogeneous with an ash fraction within the physiological range. If such dependance is not considered and fatigue properties are assumed constant, the homogenization is not always achieved and the mineral content may rise up to high non-physiological values. Thus, the interconnection between mineral content and fatigue properties is essential for the maintenance of bone's structural integrity as well as for the calcium homeostasis.     

Homogenization theory and digital imaging: A basis for studying the mechanics and design principles of bone tissue

Bone tissue is a complex multilevel composite which has the ability to sense ad respond to its mechanical environment. It is believed that bone cells called osteocytes within the bone matrix sense the mechanical environment and determine whether structural alterations are needed. At present it is not known, however, how loads are transferred from the whole bone level to cells. A computational procedure combining representative volume element (RVE) based homogenization theory with digital imaging is proposed to estimate strains at various levels of bone structure. Bone tissue structural organization and RVE based analysis are briefly reviewed. The digital image based computational procedure was applied to estimate strains in individual trabeculae (first-level microstructure). Homogenization analysis of an idealized model was used to estimate strains at one level of bone structure around osteocyte lacunae (second-level trabecular microstructure). The results showed that strain at one level of bone structure is amplified to a broad range at the next microstructural level. In one case, a zeor-level tensile principal strain of 495 muE engendered strains ranging between -1000 and 7000 muE in individual trabeculae (first-level microstructure). Subsequently, a first-level tensile principal strains of 1325 muE within an inidividual trabecula engendered strains ranging between 782 and 2530 muE around osteocyte lacunae. Lacunar orientation was found to influence strains around osteocyte lacunae much more than lacunar ellipticity. In conclusion, the computational procedure combining homogenization theory with digital imaging can proveide estimates of cell level strains within whole bones. Such results may be used to bridge experimental studies of bone adaptation at the whole bone and cell culture level. (c) 1994 John Wiley & Sons, Inc.     

Effect of mineral dissolution from bone specimens on the viscoelastic properties of cortical bone

Although physiological saline (0.15 M NaCl aqueous solution) has been used as a storage solution to prevent bone specimens from drying, there have been reports that Ca²⁺ ions dissolve from bone specimens during the storage in saline. In order to determine whether such storage has a marked effect on mechanical properties, the relaxation modulus of bovine cortical bone stored in physiological saline was compared with that stored in a buffer solution containing a sufficient amount of Ca²⁺ ions. After storage in saline, the modulus value of specimens was significantly reduced from that before storage. On the other hand, the modulus value of specimens soaked in the solution containing sufficient Ca²⁺ ions did not change after storage. The relaxation rate of a bone specimen stored in physiological saline was larger than that of a specimen stored in Ca²⁺-buffered saline solution and that of the control specimens. The results suggest that by the dissolution of Ca²⁺ ions from a bone specimen during storage in physiological saline, percolated paths of mineral phase and of reinforced matrix phase are disjoined, resulting in reduction in the elastic modulus and change in the viscoelastic properties of bone.     

Human temporal bone from the lower Paleolithic site of Castel di Guido, near Rome, Italy

The Castel di Guido site, with an estimated age of approximately 300,000 years, has yielded abundant animal remains, Acheulian stone and bone bifaces, and small tools. On the surface of the original deposit, turned over by agricultural activities, fragments of human remains were discovered between 1980 and 1986, including a temporal bone described here.     

How and why humans grow thin skulls: Experimental evidence for systemic cortical robusticity

To what extent is cranial vault thickness (CVT) a character that is strongly linked to the genome, or to what extent does it reflect the activity of an individual prior to skeletal maturity? Experimental data from pigs and armadillos indicate that CVT increases more rapidly in exercised juveniles than in genetically similar controls, despite the low levels of strain generated by chewing or locomotion in the neurocranium. CVT increases in these individuals appear to be a consequence of systemic cortical bone growth induced by exercise. In addition, an analysis of the variability in vault thickness in the genus Homo demonstrates that, until the Holocene, there has been only a slight, general decrease in vault thickness over time with no consistent significant differences between archaic and early anatomically modern humans from the Late Pleistocene. Although there may be some genetic component to variation in CVT, exercise-related, non-genetically heritable stimuli appear to account for most of the variance between individuals. The thick cranial vaults of most hunter-gatherers and early agriculturalists suggests that they may have experienced higher levels of sustained exercise relative to body mass than the majority of recent, post-industrial humans. © 1996 Wiley-Liss, Inc.     

Flow and Oxygen Transfer Characteristics of an Intravenous Membrane Oxygenator: A Computational Study

Abstract

Spectral element computational simulations of the conservation of mass, momentum and species equations are performed to investigate the flow and oxygen transfer characteristics of an Intravenous Membrane Oxygenator (IMO). The simulations consider a three-dimensional IMO computational model consisting of equally-spaced fibers, an elastic balloon with non-permeable walls positioned longitudinally within the vena cava, and a Newtonian and time-dependent incompressible flow. Flow characteristics and oxygen transfer parameters are determined for operating conditions of a stationary and a pulsating balloon. For the stationary balloon configuration the flow is two-dimensional, parallel, laminar and without secondary flows for the Reynolds number range of 5.7-455.2. Evaluations of the oxygen transfer characteristics for the stationary balloon indicate that the main transport mechanisms are diffusion and convection in the crosswise and streamwise directions, respectively. Additionally, evaluations of oxygen transfer rates and Sherwood numbers in this Reynolds number range indicate that the oxygen transfer rate reaches an asymptotic limit at relatively moderate Reynolds numbers. For the pulsating balloon, flow characteristic results demonstrate the existence of a strong secondary flow around the fiber, and between the balloon and the fiber. This secondary flow induces oscillatory crosswise and streamwise velocities and a seemingly random spanwise flow which enhances the flow mixing as well as the transport of oxygen from the fiber surface to the bulk flow.     

A novel in vitro retinal differentiation model by co-culturing adult human bone marrow stem cells with retinal pigmented epithelium cells

Human retinal pigment epithelium (HRPE) cells are important in maintaining the normal physiology within the neurosensory retina and photoreceptors. Recently, transplantation of HRPE has become a possible therapeutic approach for retinal degeneration. By negative immunoselection (CD45 and glycophorin A), in this study, we have isolated and cultivated adult human bone marrow stem cells (BMSCs) with multilineage differentiation potential. After a 2- to 4-week culture under chondrogenic, osteogenic, adipogenic, and hepatogenic induction medium, these BMSCs were found to differentiate into cartilage, bone, adipocyte, and hepatocyte-like cells, respectively. We also showed that these BMSCs could differentiate into neural precursor cells (nestin-positive) and mature neurons (MAP-2 and Tuj1-positive) following treatment of neural selection and induction medium for 1 month. Furthermore, the plasticity of BMSCs was confirmed by initiating their differentiation into retinal cells and photoreceptor lineages by co-culturing with HRPE cells. The latter system provides an ex vivo expansion model of culturing photoreceptors for the treatment of retinal degeneration diseases.     

Human demineralised bone matrix as a bone substitute for reconstruction of cystic defects of the lower jaw

Ina retrospective study validated by a standardized clinical and radiologicalexamination, the bone regeneration in 90 patients with cystic mandibulardefectswas examined. In 50 patients bony defect reconstructions with humandemineralised bone matrix (HDBM) were carried out, while in a comparable groupof 40 patients the hollow pockets were left to regenerate bone spontaneously.The bone regeneration after the implantation of human demineralised bone matrix(HDBM) was subjected to a comparative validation. Osteoinductive proteinspresent in HDBM (bone morphogenetic proteins) can diffuse into the implant seatand induce new bone formation (osteoinduction). A markedly faster and morethorough bone regeneration was demonstrated after the surgical therapy ofcysticmandibular lesions with HDBM than without. HDBM also proved to be exceptionallybiocompatible.