Peri-implant Bone Remodeling after Total Hip Replacement Combined with Systemic Alendronate Treatment: A Finite Element Analysis

In order to decrease the peri-implant bone loss during the life-time of the implant, oral use of anti-osteoporosis drugs (like bisphosphonates) has been suggested.

In this study, bone remodeling parameters identified from clinical trials of alendronate were used to simulate the effect of those drugs used after total hip arthroplasty on the peri-implant bone density. Results of the simulation show that the oral administrated drugs increase bone density around the implant and decreases, at the same time, the micromovements between the implant and the surrounding bone tissue.

Incorporation of drug effect in numerical studies of bone remodeling is a promising tool especially to predetermine safe bisphosphonate doses that could be used with orthopedic implants.     

Prediction of bone density around orthopedic implants delivering bisphosphonate

The fixation of an orthopedic implant depends strongly upon its initial stability. Peri-implant bone may resorb shortly after the surgery. This resorption is directly followed by new bone formation and implants fixation strengthening, the so-called secondary fixation. If the initial stability is not reached, the resorption continues and the implant fixation weakens, which leads to implant loosening. Studies with rats and dogs have shown that a solution to prevent peri-implant resorption is to deliver bisphosphonate from the implant surface.The aims of the study were, first, to develop a model of bone remodeling around an implant delivering bisphosphonate, second, to predict the bisphosphonate dose that would induce the maximal peri-implant bone density, and third to verify in vivo that peri-implant bone density is maximal with the calculated dose.The model consists of a bone remodeling equation and a drug diffusion equation. The change in bone density is driven by a mechanical stimulus and a drug stimulus. The drug stimulus function and the other numerical parameters were identified from experimental data. The model predicted that a dose of 0.3 μg of zoledronate on the implant would induce a maximal bone density. Implants with 0.3 μg of zoledronate were then implanted in rat femurs for 3, 6 and 9 weeks. We measured that peri-implant bone density was 4% greater with the calculated dose compared to the dose empirically described as best.The approach presented in this paper could be used in the design and analysis processes of experiments in local delivery of drug such as bisphosphonate.     

Orthopaedic Implant as Drug Delivery System: a Numerical Approach

Abstract

There is a current trend to propose cementless total joint arthroplasty (TJA) to younger patients. These patients have more demanding physical activity resulting in an increased failure rate of the implants. In particular for these type of patients, the desired service life of the implant should be extended. The actual implant used do not fulfil this requirement.

In this study, a new concept of orthopaedic implant is presented where the implant is not only a structural support but also a local drug delivery system. The delivered drug is meant to influence the bone remodeling in a way so as to compensate the effects of peri-implant osteolysis. To test this concept, we extended an existing bone remodeling model to include the effect of a drug. The results show that a more homogeneous bone density distribution can be obtained around the implant. Implants used as drug delivery systems could then be an alternative way to increase implant service life.     

Predicting bone remodeling around tissue- and bone-level dental implants used in reduced bone width

The objective of this study was to predict time-dependent bone remodeling around tissue- and bone-level dental implants used in patients with reduced bone width. The remodeling of bone around titanium tissue-level, and titanium and titanium–zirconium alloy bone-level implants was studied under 100 N oblique load for one month by implementing the Stanford theory into three-dimensional finite element models. Maximum principal stress, minimum principal stress, and strain energy density in peri-implant bone and displacement in x- and y- axes of the implant were evaluated. Maximum and minimum principal stresses around tissue-level implant were higher than bone-level implants and both bone-level implants experienced comparable stresses. Total strain energy density in bone around titanium implants slightly decreased during the first two weeks of loading followed by a recovery, and the titanium–zirconium implant showed minor changes in the axial plane. Total strain energy density changes in the loading and contralateral sides were higher in tissue-level implant than other implants in the cortical bone at the horizontal plane. The displacement values of the implants were almost constant over time. Tissue-level implants were associated with higher stresses than bone-level implants. The time-dependent biomechanical outcome of titanium–zirconium alloy bone-level implant was comparable to the titanium implant.     

Predictions of bone remodeling around dental implant systems

This study presents the implementation of a mathematical bone remodeling algorithm to bone adaptation in the premolar area of the mandible around various dental implant systems, and thus sheds a new perspective to the complex interactions in dental implant mechanics. A two-dimensional, plane strain model of the bone was built from a CT-scan. The effect of implant contour on internal bone remodeling was investigated by considering four dental implant systems with contours similar to commercially available ones and another four with cylindrical and conical cross-sections. The remodeling algorithm predicts non-homogeneous density/elastic modulus distribution; and, implant contour has some effect on how this is distributed. Bone density is predicted to increase on the tips of the threads of the implants, but to decrease inside the grooves. Threadless implants favor to develop a softer bone around their periphery, compared to implant systems that have threads. The overall contour (dimensions and the shape) of an implant affect the bone density redistribution, but the differences between different implant systems are relatively small.     

Mechanisms of reduced implant stability in osteoporotic bone

The determining factors for the fixation of uncemented screws in bone are the bone-implant interface and the peri-implant bone. The goal of this work was to explore the role of the peri-implant bone architecture on the mechanics of the bone-implant system. In particular, the specific aims of the study were to investigate: (i) the impact of the different architectural parameters, (ii) the effects of disorder, and (iii) the deformations in the peri-implant region. A three-dimensional beam lattice model to describe trabecular bone was developed. Various microstructural features of the lattice were varied in a systematic way. Implant pull-out tests were simulated, and the stiffness and strength of the bone-implant system were computed. The results indicated that the strongest decrease in pull-out strength was obtained by trabecular thinning, whereas pull-out stiffness was mostly affected by trabecular removal. These findings could be explained by investigating the peri-implant deformation field. For small implant displacements, a large amount of trabeculae in the peri-implant region were involved in the load transfer from implant to bone. Therefore, trabecular removal in this region had a strong negative effect on pull-out stiffness. Conversely, at higher displacements, deformations mainly localized in the trabeculae in contact with the implant; hence, thinning those trabeculae produced the strongest decrease in the strength of the system. Although idealized, the current approach is helpful for a mechanical understanding of the role played by peri-implant bone.     

The discrete nature of trabecular bone microarchitecture affects implant stability

Small endosseous implants, such as screws, are important components of modern orthopedics and dentistry. Hence they have to reliably fulfill a variety of requirements, which makes the development of such implants challenging. Finite element analysis is a widely used computational tool used to analyze and optimize implant stability in bone. For these purposes, bone is generally modeled as a continuum material. However, bone failure and bone adaptation processes are occurring at the discrete level of individual trabeculae; hence the assessment of stresses and strains at this level is relevant. Therefore, the aim of the present study was to investigate how peri-implant strain distribution and load transfer between implant and bone are affected by the continuum assumption. We performed a computational study in which cancellous screws were inserted in continuum and discrete models of trabecular bone; axial loading was simulated. We found strong differences in bone-implant stiffness between the discrete and continuum bone model. They depended on bone density and applied boundary conditions. Furthermore, load transfer from the screw to the surrounding bone differed strongly between the continuum and discrete models, especially for low-density bone. Based on our findings we conclude that continuum bone models are of limited use for finite element analysis of peri-implant mechanical loading in trabecular bone when a precise quantification of peri-implant stresses and strains is required. Therefore, for the assessment and improvement of trabecular bone implants, finite element models which accurately represent trabecular microarchitecture should be used.     

Evaluation of bone remodeling around single dental implants of different lengths: a mechanobiological numerical simulation and validation using clinical data

Algorithmic models have been proposed to explain adaptive behavior of bone to loading; however, these models have not been applied to explain the biomechanics of short dental implants. Purpose of present study was to simulate bone remodeling around single implants of different lengths using mechanoregulatory tissue differentiation model derived from the Stanford theory, using finite elements analysis (FEA) and to validate the theoretical prediction with the clinical findings of crestal bone loss. Loading cycles were applied on 7-, 10-, or 13-mm-long dental implants to simulate daily mastication and bone remodeling was assessed by changes in the strain energy density of bone after a 3, 6, and 12 months of function. Moreover, clinical findings of marginal bone loss in 45 patients rehabilitated with same implant designs used in the simulation (n = 15) were computed to validate the theoretical results. FEA analysis showed that although the bone density values reduced over time in the cortical bone for all groups, bone remodeling was independent of implant length. Clinical data showed a similar pattern of bone resorption compared with the data generated from mathematical analyses, independent of implant length. The results of this study showed that the mechanoregulatory tissue model could be employed in monitoring the morphological changes in bone that is subjected to biomechanical loads. In addition, the implant length did not influence the bone remodeling around single dental implants during the first year of loading.     

Mandibular bone remodeling induced by dental implant

The ability to assess the effects of an implant on bone remodeling is of particular importance to prosthesis placement planning and associated treatment assurance. Prediction of on-going bone responses will enable us to improve the performance of a restoration. Although the bone remodeling for long bones had been extensively studied, there have been relatively few reports for dental scenarios despite its increasing significance with more and more dental implant placements. This paper aimed to develop a systematic protocol to assess mandibular bone remodeling induced by dental implantation, which extends the remodeling algorithms established for the long bones into dental settings. In this study, a 3D model for a segment of a human mandible was generated from in vivo CT scan images, together with a titanium implant embedded to the mandible. The results examined the changes in bone density and stiffness as a result of bone remodeling over a period of 48 months. Resonance frequency analysis was also performed to relate natural frequencies to bone remodeling. The density contours are qualitatively compared with clinical follow-up X-ray images, thereby providing validity for the bone remodeling algorithm presented in dental bone analysis.     

Finite element modelling of implant designs and cortical bone thickness on stress distribution in maxillary type IV bone

The aims of this study were to examine the effect of implant neck design and cortical bone thickness using 3D finite element analysis and to analyse the stability of clinical evidence based on micromotion and principal stress. Four commercial dental implants for a type IV bone and maxillary segments were created. Various parameters were considered, including the osseointegration condition, loading direction and cortical bone thickness. Micromotion and principal stresses were used to evaluate the failure of osseointegration and bone overloading, respectively. It was found that the maximum stress of the peri-implant bone decreased as cortical bone thickness increased. The micromotion level in full osseointegration is less than that in non-osseointegration and it also decreases as cortical bone thickness increases. The cortical bone thickness should be measured before surgery to help select a proper implant. In the early stage of implantation, the horizontal loading component induces stress concentration in bone around the implant neck more easily than does the vertical loading component, and this may result in crestal bone loss.     

A Femur-Implant Model for the Prediction of Bone Remodeling Behavior Induced by Cementless Stem

Bone remodeling simulation is an effective tool for the prediction of long-term effect of implant on the bone tissue, as well as the selection of an appropriate implant in terms of architecture and material. In this paper, a finite element model of proximal femur was developed to simulate the structures of internal trabecular and cortical bones by incorporating quantitative bone functional adaptation theory with finite element analysis. Cementless stems made of titanium, two types of Functionally Graded Material (FGM) and flexible ‘iso-elastic’ material as comparison were implanted in the structure of proximal femur respectively to simulate the bone remodeling behaviors of host bone. The distributions of bone density, von Mises stress, and interface shear stress were obtained. All the prosthetic stems had effects on the bone remodeling behaviors of proximal femur, but the degrees of stress shielding were different. The amount of bone loss caused by titanium implant was in agreement with the clinical observation. The FGM stems caused less bone loss than that of the titanium stem, in which FGM I stem (titanium richer at the top to more HAP/Col towards the bottom) could relieve stress shielding effectively, and the interface shear stresses were more evenly distributed in the model with FGM I stem in comparison with those in the models with FGM II (titanium and bioglass) and titanium stems. The numerical simulations in the present study provided theoretical basis for FGM as an appropriate material of femoral implant from a biomechanical point of view. The next steps are to fabricate FGM stem and to conduct animal experiments to investigate the effects of FGM stem on the remodeling behaviors using animal model.     

The effect of surface roughness on the stress adaptation of trabecular architecture around a cylindrical implant

The effect of implant-bone bonding and the effect of implant surface roughness on bone remodeling near the bone-implant interface were studied by using a surface remodeling theory and the boundary element method. The study has shown that implant attachment plays an important role in bone remodeling near the implant. It has been observed in animal experiments and in clinical situations that the remodeled trabecular bone architecture around a cylindrical implant could vary, on one hand, from a hub surrounding the implant with a set of external spokes to, on the other hand, a hubless situation in which a set of spokes attach directly to the implant. It is shown here that the difference in these structures may be attributed to differences in implant attachment. The results show that the bone with perfect bonding or roller boundary condition without a gap remodeled to a hubless spoke trabecular bone architecture. On the other hand, the roller boundary condition with a specified gap yielded a spoke trabecular architecture with a hub or ring surrounding the implant. These quantitative results mirror the experimental and clinical observations. It is concluded that the hub is a consequence of the gap and not a consequence of the lack of friction between the implant and the bone.     

Trabecular bone strains around a dental implant and associated micromotions—A micro-CT-based three-dimensional finite element study

The first objective of this computational study was to assess the strain magnitude and distribution within the three-dimensional (3D) trabecular bone structure around an osseointegrated dental implant loaded axially. The second objective was to investigate the relative micromotions between the implant and the surrounding bone. The work hypothesis adopted was that these virtual measurements would be a useful indicator of bone adaptation (resorption, homeostasis, formation).In order to reach these objectives, a μCT-based finite element model of an oral implant implanted into a Berkshire pig mandible was developed along with a robust software methodology. The finite element mesh of the 3D trabecular bone architecture was generated from the segmentation of μCT scans. The implant was meshed independently from its CAD file obtained from the manufacturer. The meshes of the implant and the bone sample were registered together in an integrated software environment. A series of non-linear contact finite element (FE) analyses considering an axial load applied to the top of the implant in combination with three sets of mechanical properties for the trabecular bone tissue was devised. Complex strain distribution patterns are reported and discussed. It was found that considering the Young’s modulus of the trabecular bone tissue to be 5, 10 and 15 GPa resulted in maximum peri-implant bone microstrains of about 3000, 2100 and 1400. These results indicate that, for the three sets of mechanical properties considered, the magnitude of maximum strain lies within an homeostatic range known to be sufficient to maintain/form bone. The corresponding micro-motions of the implant with respect to the bone microstructure were shown to be sufficiently low to prevent fibrous tissue formation and to favour long-term osseointegration.     

Adaptive bone-remodeling theory applied to prosthetic-design analysis

The subject of this article is the development and application of computer-simulation methods to predict stress-related adaptive bone remodeling, in accordance with 'Wolff's Law'. These models are based on the Finite Element Method (FEM) in combination with numerical formulations of adaptive bone-remodeling theories. In the adaptive remodeling models presented, the Strain Energy Density (SED) is used as a feed-back control variable to determine shape or bone density adaptations to alternative functional requirements, whereby homeostatic SED distribution is assumed as the remodeling objective. These models are applied to investigate the relation between 'stress shielding' and bone resorption in the femoral cortex around intramedullary prostheses, such as used in Total Hip Arthroplasty (THA). It is shown that the amount of bone resorption depends mainly on the rigidity and the bonding characteristics of the implant. Homeostatic SED can be obtained when the resorption process occurs at the periosteal surface, rather than inside the cortex, provided that the stem is adequately flexible.     

Multi-factorial analysis of variables influencing the bone loss of an implant placed in the maxilla: Prediction using FEA and SED bone remodeling algorithm

The aim of this study was to investigate the interactions of implant position, implant–abutment connection and loading condition influencing bone loss of an implant placed in the maxilla using finite element (FE) analysis and mathematical bone remodeling theory. The maxilla section contours were acquired using CT images to construct FE models containing RS (internal retaining-screw) and the TIS (taper integrated screwed-in) implants placed in SC (along the axis of occlusal force) and RA (along the axis of residual ridge) positions. The adaptive strain energy density (SED) algorithm was combined with FE approach to study the preliminary bone remodeling around implant systems under different load conditions. The simulated results showed that the implant position obviously influenced the bone loss. An implant placed in the RA position resulted in substantially increased bone loss. Implant receiving a lateral load slightly increased bone loss compared with an axial load. The implant type did not significantly influence bone loss. It was found that buccal site suffered the most bone loss around the implant, followed by distal, lingual and mesial sites. The implant position primarily influenced bone loss and it was found most obviously at the buccal site. Implant placed along the axial load direction of a proposed prosthesis could obtain less bone loss around the implant. Attaining proper occlusal adjustments to reduce the lateral occlusal force is recommended in implant–bone–prosthesis system. Abutments of internal engagement with or without taper-fit did not affect the bone loss in the surrounding bone.     

Simulated bone remodeling around tilted dental implants in the anterior maxilla

Dental implants have to be placed with the long axis in different angulations due to the change in bone morphology. The objective of this study was to investigate the different bone remodeling response induced by the tilted dental implants and to assess whether it could lead to bone loss and implant failure. In this study, bone remodeling due to palato-labially inclined dental implants placed in the anterior maxillary incisor region was simulated. CT-based finite element models of a maxillary bone with dental implants were created herein. Five dental implants were placed at \(+10^{\circ }\), \(+5^{\circ }\), \(0^{\circ }\), \(-5^{\circ }\) and \(-10^{\circ }\), respectively. The remodeling progression was recorded and compared. Model \(-10^{\circ }\) (palatal side) shows the highest bone density values, but the inclined implant at \(+10^{\circ }\) (labial side) leads to significant bone loss. From a biomechanical perspective, it is speculated that a palatally inclined implant is more likely to enhance the bone density in the maxillary anterior region, but labial inclination of implant could jeopardize its stability.     

Strain shielding in distal femur after patellofemoral arthroplasty under different activity conditions

Strain shielding, a mechanical effect occurring in structures combining stiff with more flexible materials, is considered to lead to a reduction of density in bone surrounding the implant. This effect can be related to the weakness of the implant fixation, which can promote implant loosening. Several studies describe a significant decrease in postoperative bone mineral density adjacent to joint implants, which can compromise their long-term fixation. The aim of the present study was to quantify the strain shielding effect on the distal femur after patellofemoral arthroplasty. For this purpose three activities of daily living were considered: level walking, stair climbing and deep bending at different angles of knee flexion. To determine the strain shielding effect, cortical bone strains were measured experimentally with triaxial strain gauges in synthetic femurs before and after patellofemoral arthroplasty for each of the different daily activities. The results showed that the patellofemoral arthroplasty in general reduced the strains in the medial and distal regions of the femur when deep bending activity occurred, consequently, strain shielding in these regions, with strain decreases of ?72.0% and ?67.5% were measured. On the other side, higher values of strain were found in the anterior region after patellofemoral replacement for this activity with an increase of +182.0%. The occurrence of strain shielding seems to be more significant when the angle of knee flexion and applied load increases. Strain shielding and over-loading may have relevant effects on bone remodeling surrounding the patellofemoral implant, suggesting a potential effect of later bone resorption in the medial and distal femur regions in case of regular deep bending activity.     

In vivo anti-osteoporotic activity of isotaxiresinol, a lignan from wood of Taxus yunnanensis

Isotaxiresinol, the main lignan isolated from the water extract of wood of Taxus yunnanensis, was investigated for its effect on bone loss, on serum biochemical markers for bone remodeling and on uterine tissue, using ovariectomized (OVX) rats as the model of postmenopausal osteoporosis. After oral administration of isotaxiresinol (50 and 100mg/kg/d) for 6 weeks, bone mineral content (BMC) and bone mineral density (BMD) in total and cortical bones were increased as compared to those of OVX control rats, and decreases of three bone strength indexes induced by OVX surgery were prevented. Serum biochemical markers for bone remodeling revealed that isotaxiresinol slightly increased bone formation and significantly inhibited bone resorption without side effect on uterine tissue. These results suggest that isotaxiresinol may be useful for treatment of postmenopausal osteoporosis, especially for prevention of bone fracture induced by estrogen deficiency.     

Modelling the early phases of bone regeneration around an endosseous oral implant

The objective of this study was to see whether a mathematical model of fracture healing was able to mimic bone formation around an unloaded screw-shaped titanium implant as it is well-believed that both processes exhibit many biological similarities. This model describes the spatio-temporal evolution of cellular activities, ranging from mesenchymal stem cell migration, proliferation, differentiation to bone formation, which are initiated and regulated by the growth factors present at the peri-implant site. For the simulations, a finite volume code was used and adequate initial and boundary conditions were applied. Two sets of analyses have been performed, in which either initial and boundary condition or model parameter values were changed with respect to the fracture healing model parameter values. For a number of combinations, the spatio-temporal evolution of bone density was well-predicted. However reducing cell proliferation rate and increasing osteoblast differentiation and osteogenic growth factor synthesis rates, the simulation results were in agreement with the experimental data.     

Importance of diameter-to-length ratio in selecting dental implants: a methodological finite element study

Implant dimensions greatly influence load transfer characteristics and the lifetime of a dental system. Excessive stresses at peri-implant area may result in bone failure. Finding the critical point at the implant–bone interface and evaluating the influence of implant diameter-to-length ratio on adjacent bone stresses makes it possible to select implant dimensions. For this, different cylindrical implants were numerically analysed using geometrical models generated from computed tomography images of mandible with osseointegrated implants. All materials were assumed to be linearly elastic and isotropic. Masticatory load was applied in its natural direction, oblique to occlusal plane. Maximum von Mises stresses were located around the implant neck at the critical point of its intersection with the plane of loading and were functions of implant diameter-to-length ratio. It was demonstrated that there exists a certain spectrum of diameter-to-length ratios, which will keep maximum bone stresses at a preset level chosen in accordance with patient's bone strength.