The use of XFEM to assess the influence of intra-cortical porosity on crack propagation

This study aimed at using eXtended finite element method (XFEM) to characterize crack growth through bone's intra-cortical pores. Two techniques were compared using Abaqus: (1) void material properties were assigned to pores; (2) multiple enrichment regions with independent crack-growth possibilities were employed. Both were applied to 2D models of transverse images of mouse bone with differing porous structures. Results revealed that assigning multiple enrichment regions allows for multiple cracks to be initiated progressively, which cannot be captured when the voids are filled. Therefore, filling pores with one enrichment region in the model will not create realistic fracture patterns in Abaqus-XFEM

Insight into differences in nanoindentation properties of bone

Nanoindentation provides the ideal framework to determine mechanical properties of bone at the tissue scale without being affected by the size, shape, and porosity of the bone. However, the values of tissue level mechanical properties vary significantly between studies. Since the differences in the bone sample, hydration state, and test parameters complicate direct comparisons across the various studies, these discrepancies in values cannot be compared directly. The objective of the current study is to evaluate and compare mechanical properties of the same bones using a broad range of testing parameters. Wild type C56BL6 mice tibiae were embedded following different processes and tested in dry and rehydrated conditions. Spherical and Berkovich indenter probes were used, and data analysis was considered within the elasto-plastic (Oliver-Pharr), viscoelastic and visco-elastic-plastic frameworks. The mean values of plane strain modulus varied significantly depending on the hydration state, probe geometry and analysis method. Indentations in dry bone analyzed using a visco-elastic-plastic approach gave values of 34 GPa. After rehydrating the same bones and indenting them with a spherical tip and utilizing a viscoelastic analysis, the mean modulus value was 4 GPa, nearly an order of magnitude smaller. Results suggest that the hydration state, probe geometry and the limitations and assumptions of each analysis method influence significantly the measured mechanical properties. This is the first time that such a systematic study has been carried out and it has been concluded that the discrepancies in the mechanical properties of bone measured by nanoindentation found in the literature should not be attributed only to the differences between the bones themselves, but also to the testing and analysis protocols

Age-related changes in mouse bone permeability

The determination of lacunar-canalicular permeability is essential for understanding local fluid flow in bone, which may indicate how bone senses changes in the mechanical environment to regulate mechano-adaptation. The estimates of lacunar-canalicular permeability found in the literature vary by up to eight orders of magnitude, and age-related permeability changes have not been measured in non-osteonal mouse bone. The objective of this study is to use a poroelastic approach based on nanoindentation data to characterize lacunar-canalicular permeability in murine bone as a function of age. Nine wild type C57BL/6 mice of different ages (2, 7 and 12 months) were used. Three tibiae from each age group were embedded in epoxy resin, cut in half and indented in the longitudinal direction in the mid-cortex using two spherical fluid indenter tips (R=238 μm and 500 μm). Results suggest that the lacunar-canalicular intrinsic permeability of mouse bone decreases from 2 to 7 months, with no significant changes from 7 to 12 months. The large indenter tip imposed larger contact sizes and sampled larger ranges of permeabilities, particularly for the old bone. This age-related difference in the distribution was not seen for indents with the smaller radius tip. We conclude that the small tip effectively measured lacunar-canalicular permeability, while larger tip indents were influenced by vascular permeability. Exploring the age-related changes in permeability of bone measured by nanoindentation will lead to a better understanding of the role of fluid flow in mechano-transduction. This understanding may help indicate alterations in bone adaptation and remodeling

Analytical model for the prediction of permeability of triply periodic minimal surfaces

Triply periodic minimal surfaces (TPMS) are mathematically defined cellular structures whose geometry can be quickly adapted to target desired mechanical response (structural and fluid). This has made them desirable for a wide range of bioengineering applications; especially as bioinspired materials for bone replacement. The main objective of this study was to develop a novel analytical framework which would enable calculating permeability of TPMS structures based on the desired architecture, pore size and porosity. To achieve this, computer-aided designs of three TPMS structures (Fisher-Koch S, Gyroid and Schwarz P) were generated with varying cell size and porosity levels. Computational Fluid Dynamics (CFD) was used to calculate permeability for all models under laminar flow conditions. Permeability values were then used to fit an analytical model dependent on geometry parameters only. Results showed that permeability of the three architectures increased with porosity at different rates, highlighting the importance of pore distribution and architecture. The computed values of permeability fitted well with the suggested analytical model (R2>0.99, p<0.001). In conclusion, the novel analytical framework presented in the current study enables predicting permeability values of TPMS structures based on geometrical parameters within a difference <5%. This model, which could be combined with existing structural analytical models, could open new possibilities for the smart optimisation of TPMS structures for biomedical applications where structural and fluid flow properties need to be optimised

Studies of chain substitution caused sub-fibril level differences in stiffness and ultrastructure of wildtype and oim/oim collagen fibers using multifrequency-AFM and molecular modeling.

Molecular alteration in type I collagen, i.e., substituting the α2 chain with α1 chain in tropocollagen molecule, can cause osteogenesis imperfecta (OI), a brittle bone disease, which can be represented by a mouse model (oim/oim). In this work, we use dual-frequency Atomic Force Microscopy (AFM) and incorporated with molecular modeling to quantify the ultrastructure and stiffness of the individual native collagen fibers from wildtype (+/+) and oim/oim diseased mice humeri. Our work presents direct experimental evidences that the +/+ fibers have highly organized and compact ultrastructure and corresponding ordered stiffness distribution. In contrast, oim/oim fibers have ordered but loosely packed ultrastructure with uncorrelated stiffness distribution, as well as local defects. The molecular model also demonstrates the structural and molecular packing differences between +/+ and oim/oim collagens. The molecular mutation significantly altered sub-fibril structure and mechanical property of collagen fibers. This study can give the new insight for the mechanisms and treatment of the brittle bone disease

Predictors of diabetes risk in urban and rural areas in Colombia

Background: Nutritional habits low in fruits and vegetables and sedentary lifestyle are associated with a higher risk of developing Type 2 Diabetes (T2D). However, it is important to assess differences between urban and rural areas. This study aimed to analyze the associations between the risk of developing T2D and setting in the Colombian north coast in 2017. Methods: This cross-sectional study included 1,005 subjects. Data was collected by interviewing self-identified members of an urban community and a rural-indigenous population. The interaction terms were evaluated as well as the confounders. Then, adjusted binary logistic regressions were used to estimate the odds ratio (OR) and 95% Confidence Intervals (CI). Results: subjects with a high risk of T2D are more likely to belong to the urban setting (OR = 1.908; 95%CI = 1.201-2.01) compared with those with lower T2D after adjusting for age, Body Mass Index (BMI), physical activity, history of high levels of glycemia, and diabetes in relatives. Conclusions: Urban communities are more likely to have T2D compared with rural-indigenous populations. These populations have differences from the cultural context, including personal, and lifestyle factors.Peer reviewe

In silico assessment of the bone regeneration potential of complex porous scaffolds.

Mechanical environment plays a crucial role in regulating bone regeneration in bone defects. Assessing the mechanobiological behavior of patient-specific orthopedic scaffolds in-silico could help guide optimal scaffold designs, as well as intra- and post-operative strategies to enhance bone regeneration and improve implant longevity. Additively manufactured porous scaffolds, and specifically triply periodic minimal surfaces (TPMS), have shown promising structural properties to act as bone substitutes, yet their ability to induce mechanobiologially-driven bone regeneration has not been elucidated. The aim of this study is to i) explore the bone regeneration potential of TPMS scaffolds made of different stiffness biocompatible materials, to ii) analyze the influence of pre-seeding the scaffolds and increasing the post-operative resting period, and to iii) assess the influence of patient-specific parameters, such as age and mechanosensitivity, on outcomes. To perform this study, an in silico model of a goat tibia is used. The bone ingrowth within the scaffold pores was simulated with a mechano-driven model of bone regeneration. Results showed that the scaffold's architectural properties affect cellular diffusion and strain distribution, resulting in variations in the regenerated bone volume and distribution. The softer material improved the bone ingrowth. An initial resting period improved the bone ingrowth but not enough to reach the scaffold's core. However, this was achieved with the implantation of a pre-seeded scaffold. Physiological parameters like age and health of the patient also influence the bone regeneration outcome, though to a lesser extent than the scaffold design. This analysis demonstrates the importance of the scaffold's geometry and its material, and highlights the potential of using mechanobiological patient-specific models in the design process for bone substitutes

Understanding the influence of surgical parameters on craniofacial surgery outcomes: a computational study

Sagittal craniosynostosis (SC) is a congenital condition whereby the newborn skull develops abnormally owing to the premature ossification of the sagittal suture. Spring-assisted cranioplasty (SAC) is a minimally invasive surgical technique to treat SC, where metallic distractors are used to reshape the newborn’s head. Although safe and effective, SAC outcomes remain uncertain owing to the limited understanding of skull−distractor interaction and the limited information provided by the analysis of single surgical cases. In this work, an SC population-averaged skull model was created and used to simulate spring insertion by means of the finite-element analysis using a previously developed modelling framework. Surgical parameters were varied to assess the effect of osteotomy and spring positioning, as well as distractor combinations, on the final skull dimensions. Simulation trends were compared with retrospective measurements from clinical imaging (X-ray and three-dimensional photogrammetry scans). It was found that the on-table post-implantation head shape change is more sensitive to spring stiffness than to the other surgical parameters. However, the overall end-of-treatment head shape is more sensitive to spring positioning and osteotomy size parameters. The results of this work suggest that SAC surgical planning should be performed in view of long-term results, rather than immediate on-table reshaping outcomes

Towards a radiation free numerical modelling framework to predict spring assisted correction of scaphocephaly

Sagittal Craniosynostosis (SC) is a congenital craniofacial malformation, involving premature sagittal suture ossification; spring-assisted cranioplasty (SAC) – insertion of metallic distractors for skull reshaping – is an established method for treating SC. Surgical outcomes are predictable using numerical modelling, however published methods rely on computed tomography (CT) scans availability, which are not routinely performed. We investigated a simplified method, based on radiation-free 3D stereophotogrammetry scans.Eight SAC patients (age 5.1 ± 0.4 months) with preoperative CT and 3D stereophotogrammetry scans were included. Information on osteotomies, spring model and post-operative spring opening were recorded. For each patient, two preoperative models (PREOP) were created: i) CT model and ii) S model, created by processing patient specific 3D surface scans using population averaged skin and skull thickness and suture locations. Each model was imported into ANSYS Mechanical (Analysis System Inc., Canonsburg, PA) to simulate spring expansion. Spring expansion and cranial index (CI - skull width over length) at times equivalent to immediate postop (POSTOP) and follow up (FU) were extracted and compared with in-vivo measurements.Overall expansion patterns were very similar for the 2 models at both POSTOP and FU. Both models had comparable outcomes when predicting spring expansion. Spring induced CI increase was similar, with a difference of 1.2%±0.8% for POSTOP and 1.6%±0.6% for FU.This work shows that a simplified model created from the head surface shape yields acceptable results in terms of spring expansion prediction. Further modelling refinements will allow the use of this predictive tool during preoperative planning

Three-Dimensional Handheld Scanning to Quantify Head-Shape Changes in Spring-Assisted Surgery for Sagittal Craniosynostosis

Three-dimensional (3D) imaging is an important tool for diagnostics, surgical planning, and evaluation of surgical outcomes in craniofacial procedures. Gold standard for acquiring 3D imaging is computed tomography that entails ionizing radiations and, in young children, a general anaesthesia. Three-dimensional photographic imaging is an alternative method to assess patients who have undergone calvarial reconstructive surgery. The aim of this study was to assess the utility of 3D handheld scanning photography in a cohort of patients who underwent spring-assisted correction surgery for scaphocephaly. Pre- and postoperative 3D scans acquired in theater and at the 3-week follow-up in clinic were postprocessed for 9 patients. Cephalic index (CI), head circumference, volume, sagittal length, and coronal width over the head at pre-op, post-op, and follow-up were measured from the 3D scans. Cephalic index from 3D scans was compared with measurements from planar x-rays. Statistical shape modeling (SSM) was used to calculate the 3D mean anatomical head shape of the 9 patients at the pre-op, post-op, and follow-up. No significant differences were observed in the CI between 3D and x-ray. Cephalic index, volume, and coronal width increased significantly over time. Mean shapes from SSM visualized the overall and regional 3D changes due to the expansion of the springs in situ. Three-dimensional handheld scanning followed by SSM proved to be an efficacious and practical method to evaluate 3D shape outcomes after spring-assisted cranioplasty in individual patients and the population