Free Achilles tendon strain during selected rehabilitation, locomotor, jumping, and landing tasks

A better understanding of the strains experienced by the Achilles tendon during commonly prescribed exercises and locomotor tasks is needed to improve efficacy of Achilles tendon training and rehabilitation programs. The aim of this study was to estimate in vivo free Achilles tendon strain during selected rehabilitation, locomotor, jumping, and landing tasks. Sixteen trained runners with no symptoms of Achilles tendinopathy participated in this study. Personalised free Achilles tendon moment arm and force-strain curve were obtained from imaging data and used in conjunction with motion capture and surface electromyography to estimate free Achilles tendon strain using electromyogram-informed neuromusculoskeletal modelling. There was a strong correspondence between Achilles tendon force estimates from the present study and experimental data reported in the literature (R² > 0.85). The average tendon strain was highest for maximal hop landing (8.8±1.6%), lowest for walking at 1.4 m/s (3.1±0.8%) and increased with locomotor speed during running (run 3.0 m/s: 6.5±1.6%; run 5.0 m/s: 7.9±1.7%) and during heel rise exercise with added mass (BW: 5.8±1.3%; 1.2 BW: 6.9±1.7%). The peak tendon strain was highest during running (5 m/s: 13.7±2.5%) and lowest during walking (1.4 m/s: 7±1.8%). Overall findings provide a preliminary evidence base for exercise selection to maximise anabolic tendon remodelling during training and rehabilitation of the Achilles tendon.

Targeted Achilles Tendon Training and Rehabilitation Using Personalized and Real-Time Multiscale Models of the Neuromusculoskeletal System

Musculoskeletal tissues, including tendons, are sensitive to their mechanical environment, with both excessive and insufficient loading resulting in reduced tissue strength. Tendons appear to be particularly sensitive to mechanical strain magnitude, and there appears to be an optimal range of tendon strain that results in the greatest positive tendon adaptation. At present, there are no tools that allow localized tendon strain to be measured or estimated in training or a clinical environment. In this paper, we first review the current literature regarding Achilles tendon adaptation, providing an overview of the individual technologies that so far have been used in isolation to understand in vivo Achilles tendon mechanics, including 3D tendon imaging, motion capture, personalized neuromusculoskeletal rigid body models, and finite element models. We then describe how these technologies can be integrated in a novel framework to provide real-time feedback of localized Achilles tendon strain during dynamic motor tasks. In a proof of concept application, Achilles tendon localized strains were calculated in real-time for a single subject during walking, single leg hopping, and eccentric heel drop. Data was processed at 250 Hz and streamed on a smartphone for visualization. Achilles tendon peak localized strains ranged from ∼3 to ∼11% for walking, ∼5 to ∼15% during single leg hop, and ∼2 to ∼9% during single eccentric leg heel drop, overall showing large strain variation within the tendon. Our integrated framework connects, across size scales, knowledge from isolated tendons and whole-body biomechanics, and offers a new approach to Achilles tendon rehabilitation and training. A key feature is personalization of model components, such as tendon geometry, material properties, muscle geometry, muscle-tendon paths, moment arms, muscle activation, and movement patterns, all of which have the potential to affect tendon strain estimates. Model personalization is important because tendon strain can differ substantially between individuals performing the same exercise due to inter-individual differences in these model components.

Applying Physiologically Relevant Strains to Tenocytes in an In Vitro Cell Device Induces In Vivo Like Behaviors

We have developed a novel cell stretching device (called Cell Gym) capable of applying physiologically relevant low magnitude strains to tenocytes on a collagen type I coated membrane. We validated our device thoroughly on two levels: (1) substrate strains, (2) cell level strains. Our cell level strain results showed that the applied stretches were transferred to cells accurately (∼90%). Our gene expression data showed that mechanically stimulated tenocytes (4%) expressed a lower level of COL I gene. COX2 gene was increased but did not reach statistical significance. Our device was then tested to see if it could reproduce results from an in vivo study that measured time-dependent changes in collagen synthesis. Our results showed that collagen synthesis peaked at 24 hrs after exercise and then decreased, which matched the results from the in vivo study. Our study demonstrated that it is important to incorporate physiologically relevant low strain magnitudes in in vitro cell mechanical studies and the need to validate the device thoroughly to operate the device at small strains. This device will be used in designing novel tendon tissue engineering scaffolds in the future.

Analyzing bone remodeling patterns after total hip arthroplasty using quantitative computed tomography and patient-specific 3D computational models

BACKGROUND

Computational models in the form of finite element analysis technique that incorporates bone remodeling theories along with DEXA scans has been extensively used in predicting bone remodeling patterns around the implant. However, majority of such studies used generic models. Therefore, the aim of this study is to develop patient-specific finite element models of total hip replacement patients using their quantitative computed tomography (QCT) scans and accurately analyse bone remodelling patterns after total hip arthroplasty (THA).

METHODS

Patient-specific finite element models have been generated using the patients' QCT scans from a previous clinical follow-up study. The femur was divided into five regions in proximal-distal direction and then further divided into four quadrants for detailed analysis of bone remodeling patterns. Two types of analysis were performed-inter-patient and intra patient to compare them and then the resulting bone remodeling patterns were quantitatively analyzed.

RESULTS

Our results show that cortical bone density decrease is higher in diaphyseal region over time and the cancellous bone density decreases significantly in metaphyseal region over time. In metaphyseal region, posterior-medial (P-M) quadrant showed high bone loss while diaphyseal regions show high bone loss in anterior-lateral (A-L) quadrant.

CONCLUSIONS

Our study demonstrated that combining QCT with 3D patient-specific models has the ability of monitoring bone density change patterns after THA in much finer details. Future studies include using these findings for the development of a bone remodelling algorithm capable of predicting surgical outcomes for THA patients.

Altered load transfer in the pelvis in the presence of periprosthetic osteolysis

Periprosthetic osteolysis in the retroacetabular region with cancellous bone loss is a recognized phenomenon in the long-term follow-up of total hip replacement. The effects on load transfer in the presence of defects are less well known. A validated, patient-specific, 3D finite element (FE) model of the pelvis was used to assess changes in load transfer associated with periprosthetic osteolysis adjacent to a cementless total hip arthroplasty (THA) component. The presence of a cancellous defect significantly increased (p < 0.05) von Mises stress in the cortical bone of the pelvis during walking and a fall onto the side. At loads consistent with single leg stance, this was still less than the predicted yield stress for cortical bone. During higher loads associated with a fall onto the side, highest stress concentrations occurred in the superior and inferior pubic rami and in the anterior column of the acetabulum with larger cancellous defects.

Integrating degenerative mechanisms in bone and cartilage: a multiscale approach

At the whole organ level, degenerative mechanisms in bone and cartilage are primarily attributed to modifications in loading pattern. Either a change in magnitude or location can initiate a degenerative path. At the micro scale we often see changes in structure such as porosity increase in bone and fibrillation in cartilage. These changes contribute to a reduced structural integrity that weakens the bulk strength of tissue. Finally, at the cell level we have modeling and remodeling pathways that may be disrupted through disease, drugs and altered stimulus from the micro and macro scales. In order to understand this entire process and the roles each level plays a multiscale modeling framework is necessary. This framework can take whole body loadings and pass information through finer spatial scales in order to understand how everyday dynamic movements influence micro and cellular response. In a similar manner, cellular and microstructural processes regulate whole bulk properties and modify whole organ strength. In this study we highlight the multiscale links developed as part of the open-source ontologies for the Physiome Project using the lower limb as an example. We consider the influence of remodeling in (i) anabolic treatments in cortical bone; and (ii) subchondral bone and cartilage degeneration.

Quantitative CT with finite element analysis: towards a predictive tool for bone remodelling around an uncemented tapered stem

PURPOSE

We used quantitative CT in conjunction with finite element analysis to provide a new tool for assessment of bone quality after total hip arthroplasty in vivo. The hypothesis of this prospective five-year study is that the combination of the two modalities allows 3D patient-specific imaging of cortical and cancellous bone changes and stress shielding.

METHOD

We tested quantitative CT in conjunction with finite elements on a cohort of 29 patients (31 hips) who have been scanned postoperatively and at one year, two years and five years follow-up. The method uses cubic Hermite finite element interpolation for efficient mesh generation directly from qCT datasets. The element Gauss points that are used for the geometric interpolation functions are also used for interpolation of osteodensitometry data.

RESULTS

The study showed changes of bone density suggestive of proximal femur diaphysis load transfer with osteointegration and moderate metaphyseal stress shielding. Our model revealed that cortical bone initially became porous in the greater trochanter, but this phenomenon progressed to the cortex of the lesser trochanter and the posterior aspect of the metaphysis. The diaphyseal area did not experience major change in bone density for either cortical or cancellous bone.

CONCLUSION

The combination of quantitative CT with finite element analysis allows visualization of changes to bone density and architecture. It also provides correlation of bone density/architectural changes with stress patterns enabling the study of the effects of stress shielding on bone remodelling in vivo. This technology can be useful in predicting bone remodeling and the quality of implant fixation using prostheses with different design and/or biomaterials.

An efficient and accurate prediction of the stability of percutaneous fixation of acetabular fractures with finite element simulation

Posterior wall fracture is one of the most common fracture types of the acetabulum and a conventional approach is to perform open reduction and internal fixation with a plate and screws. Percutaneous screw fixations, on the other hand, have recently gained attention due to their benefits such as less exposure and minimization of blood loss. However their biomechanical stability, especially in terms interfragmentary movement, has not been investigated thoroughly. The aims of this study are twofold: (1) to measure the interfragmentary movements in the conventional open approach with plate fixations and the percutaneous screw fixations in the acetabular fractures and compare them; and (2) to develop and validate a fast and efficient way of predicting the interfragmentary movement in percutaneous fixation of posterior wall fractures of the acetabulum using a 3D finite element (FE) model of the pelvis. Our results indicate that in single fragment fractures of the posterior wall of the acetabulum, plate fixations give superior stability to screw fixations. However screw fixations also give reasonable stability as the average gap between fragment and the bone remained less than 1 mm when the maximum load was applied. Our finite element model predicted the stability of screw fixation with good accuracy. Moreover, when the screw positions were optimized, the stability predicted by our FE model was comparable to the stability obtained by plate fixations. Our study has shown that FE modeling can be useful in examining biomechanical stability of osteosynthesis and can potentially be used in surgical planning of osteosynthesis.

A multiscale framework based on the physiome markup languages for exploring the initiation of osteoarthritis at the bone-cartilage interface

The initiation of osteoarthritis (OA) has been linked to the onset and progression of pathologic mechanisms at the cartilage-bone interface. Most importantly, this degenerative disease involves cross-talk between the cartilage and subchondral bone environments, so an informative model should contain the complete complex. In order to evaluate this process, we have developed a multiscale model using the open-source ontologies developed for the Physiome Project with cartilage and bone descriptions at the cellular, micro, and macro levels. In this way, we can effectively model the influence of whole body loadings at the macro level and the influence of bone organization and architecture at the micro level, and have cell level processes that determine bone and cartilage remodeling. Cell information is then passed up the spatial scales to modify micro architecture and provide a macro spatial characterization of cartilage inflammation. We evaluate the framework by linking a common knee injury (anterior cruciate ligament deficiency) to proinflammatory mediators as a possible pathway to initiate OA. This framework provides a "virtual bone-cartilage" tool for evaluating hypotheses, treatment effects, and disease onset to inform and strengthen clinical studies.

Femoral bone density changes after total hip arthroplasty with uncemented taper-design stem: a five year follow-up study

We measured bone density (BD) changes to assess adaptive bone remodelling five years after uncemented total hip arthroplasty with taper-design femoral component using quantitative computed-tomography-assisted osteodensitometry (qCT). Nineteen consecutive patients (21 hips) with degenerative joint disease were enrolled in the study. A press-fit cup and a tapered uncemented stem ceramic-ceramic pairing were used in all patients. Serial clinical, radiological and qCT osteodensitometry assessments were performed after the index operation and at the one, two and five year follow-ups. At the latest follow-up, the clinical outcome was rated satisfactory in all hips. The radiological assessment showed signs of osteointegration with stable fixation of all cups and stems. Overall, there was evidence of a BD loss at year five (p = 0.004). We estimate that BD loss was between 2.2% and 12.1% in comparison with baseline postoperative values. Progressive loss of BD in the metaphyseal region was observed in all hips. We found unremarkable BD changes of diaphyseal cortical BD throughout the five year follow-up period. qCT osteodensitometry technology allows differentiation of cortical and cancellous BD changes over time. Periprosthetic BD changes at the five year follow-up are suggestive of stable stem osteointegration with proximal femoral diaphysis load transfer and metaphyseal stress shielding.

Development and Validation of Patient-Specific Finite Element Models of the Hemipelvis Generated From a Sparse CT Data Set

To produce a patient-specific finite element (FE) model of a bone such as the pelvis, a complete computer tomographic (CT) or magnetic resonance imaging (MRI) geometric data set is desirable. However, most patient data are limited to a specific region of interest such as the acetabulum. We have overcome this problem by providing a hybrid method that is capable of generating accurate FE models from sparse patient data sets. In this paper, we have validated our technique with mechanical experiments. Three cadaveric embalmed pelves were strain gauged and used in mechanical experiments. FE models were generated from the CT scans of the pelves. Material properties for cancellous bone were obtained from the CT scans and assigned to the FE mesh using a spatially varying field embedded inside the mesh while other materials used in the model were obtained from the literature. Although our FE meshes have large elements, the spatially varying field allowed them to have location dependent inhomogeneous material properties. For each pelvis, five different FE meshes with a varying number of patient CT slices (8–12) were generated to determine how many patient CT slices are needed for good accuracy. All five mesh types showed good agreement between the model and experimental strains. Meshes generated with incomplete data sets showed very similar stress distributions to those obtained from the FE mesh generated with complete data sets. Our modeling approach provides an important step in advancing the application of FE models from the research environment to the clinical setting.

The use of sparse CT datasets for auto-generating accurate FE models of the femur and pelvis

The finite element (FE) method when coupled with computed tomography (CT) is a powerful tool in orthopaedic biomechanics. However, substantial data is required for patient-specific modelling. Here we present a new method for generating a FE model with a minimum amount of patient data. Our method uses high order cubic Hermite basis functions for mesh generation and least-square fits the mesh to the dataset. We have tested our method on seven patient data sets obtained from CT assisted osteodensitometry of the proximal femur. Using only 12 CT slices we generated smooth and accurate meshes of the proximal femur with a geometric root mean square (RMS) error of less than 1 mm and peak errors less than 8 mm. To model the complex geometry of the pelvis we developed a hybrid method which supplements sparse patient data with data from the visible human data set. We tested this method on three patient data sets, generating FE meshes of the pelvis using only 10 CT slices with an overall RMS error less than 3 mm. Although we have peak errors about 12 mm in these meshes, they occur relatively far from the region of interest (the acetabulum) and will have minimal effects on the performance of the model. Considering that linear meshes usually require about 70-100 pelvic CT slices (in axial mode) to generate FE models, our method has brought a significant data reduction to the automatic mesh generation step. The method, that is fully automated except for a semi-automatic bone/tissue boundary extraction part, will bring the benefits of FE methods to the clinical environment with much reduced radiation risks and data requirement.