An Adaptive Pedaling Assistive Device for Asymmetric Torque Assistant in Cycling.

Dynamic loads have short and long-term effects in the rehabilitation of lower limb joints. However, an effective exercise program for lower limb rehabilitation has been debated for a long time. Cycling ergometers were instrumented and used as a tool to mechanically load the lower limbs and track the joint mechano-physiological response in rehabilitation programs. Current cycling ergometers apply symmetrical loading to the limbs, which may not reflect the actual load-bearing capacity of each limb, as in Parkinson's and Multiple Sclerosis diseases. Therefore, the present study aimed to develop a new cycling ergometer capable of applying asymmetric loads to the limbs and validate its function using human tests. The instrumented force sensor and crank position sensing system recorded the kinetics and kinematics of pedaling. This information was used to apply an asymmetric assistive torque only to the target leg using an electric motor. The performance of the proposed cycling ergometer was studied during a cycling task at three different intensities. It was shown that the proposed device reduced the pedaling force of the target leg by 19% to 40%, depending on the exercise intensity. This reduction in pedal force caused a significant reduction in the muscle activity of the target leg (p < 0.001), without affecting the muscle activity of the non-target leg. These results demonstrated that the proposed cycling ergometer device is capable of applying asymmetric loading to lower limbs, and thus has the potential to improve the outcome of exercise interventions in patients with asymmetric function in lower limbs.

Gait Phase Detection in Walking and Stairs Using Machine Learning.

Machine learning-based activity and gait phase recognition algorithms are used in powered motion assistive devices to inform control of motorized components. The objective of this study was to develop a supervised multiclass classifier to simultaneously detect activity and gait phase (stance, swing) in real-world walking, stair ascent, and stair descent using inertial measurement data from the thigh and shank. The intended use of this algorithm was for control of a motion assistive device local to the knee. Using data from 80 participants, two decision trees and five long short-term memory (LSTM) models that each used different feature sets were initially tested and evaluated using a novel performance metric: proportion of perfectly classified strides (PPCS). Based on the PPCS of these initial models, five additional posthoc LSTM models were tested. Separate models were developed to classify (i) both activity and gait phase simultaneously (one model predicting six states), and (ii) activity-specific models (three individual binary classifiers predicting stance/swing phases). The superior activity-specific model had an accuracy of 98.0% and PPCS of 55.7%. The superior six-phase model used filtered inertial measurement data as its features and a median filter on its predictions and had an accuracy of 92.1% and PPCS of 22.9%. Pooling stance and swing phases from all activities and treating this model as a binary classifier, this model had an accuracy of 97.1%, which may be acceptable for real-world lower limb exoskeleton control if only stance and swing gait phases must be detected. Keywords: machine learning, deep learning, inertial measurement unit, activity recognition, gait.

Effect of Simulated Mass-Tunable Auxetic Midsole on Vertical Ground Reaction Force.

When runners impact the ground, they experience a sudden peak ground reaction force (GRF), which may be up to 4× greater than their bodyweight. Increased GRF impact peak magnitude has been associated with lower limb injuries in runners. Yet, shoe midsoles are capable of cushioning the impact between the runner and the ground to reduce GRF. It has been proposed that midsoles should be tunable with subject mass to minimize GRF and reduce risk of injury. Auxetic metamaterials, structures designed to achieve negative Poisson's ratios, demonstrate superior impact properties and are highly tunable. Recently, auxetic structures have been introduced in footwear, but their effects on GRF are not documented in literature. This work investigates the viability of a three-dimensional auxetic impact structure with a tunable force plateau as a midsole through mass-spring-damper simulation. An mass-spring-damper model was used to perform 315 simulations considering combinations of seven subject masses (45-90 kg), 15 auxetic plateau forces (72-1080 N), and three auxetic damping conditions (450, 725, and 1000 Ns/m) and regression analysis was used to determine their influence on GRF impact peak, energy, instantaneous, and average loading rate. Simulations showed that tuning auxetic plateau force and damping based on subject mass may reduce GRF impact and loading rate versus simulated conventional midsoles. Auxetic plateau force and damping conditions of 450 Ns/m and ∼1 bodyweight (BW), respectively, minimized peak impact GRF. This work demonstrates the need for tunable auxetic midsoles and may inform future work involving midsole testing.

Vancouver B and C periprosthetic fractures around the cemented Exeter Stem: sex is associate with fracture pattern.

INTRODUCTION: The aim of this study was to identify factors associated with the level of periprosthetic fracture involving a cemented polished tapered stem: Vancouver B or Vancouver C. METHODS: A retrospective cohort study of 181 unilateral periprosthetic fractures involving Exeter stems was assessed by three observers (mean age 78.5, range 39-103; mean BMI 27.1, 17-39; 97 (54%) male). Patient demographics, deprivation scores, BMI and time since primary prosthesis were recorded. Femoral diameter, femoral cortical thickness, Dorr classification and distal cement mantle length were measured from calibrated radiographs. Interobserver reliability was calculated using intraclass correlation coefficients (ICCs). Univariate and multivariate analysis was performed to identify associations with Vancouver B or C fractures. RESULTS: 160/181 (88%) Vancouver B and 21/181 (12%) Vancouver C-level fractures occurred at a mean of 5.9 ± 5.4 years (0.2-26.5) following primary surgery. Radiographic measurements demonstrated excellent agreement (ICC > 0.8, p < 0.001). Mortality was significantly higher following Vancouver C compared to B fractures: 90 day 14/160 Vs 5/21 (p = 0.05); 1 year 29/160 Vs 8/21 (p = 0.03). Univariate analysis demonstrated that Vancouver C fractures were associated with female sex, bisphosphonate use, cortical bone thickness, and distal cement mantle length (p < 0.05). On multivariate analysis, only female sex was an independent predictor of Vancouver C-level fractures (R2 =0.354, p = 0.005). CONCLUSION: Most PFFs involving the Exeter stem design are Vancouver B-type fractures and appear to be independent of osteoporosis. In contrast, Vancouver C periprosthetic fractures display typical fragility fracture characteristics and are associated with female sex, thinner femoral cortices, longer distal cement mantles and high mortality.

Lumbar spine loads are reduced for activities of daily living when using a braced arm-to-thigh technique.

PURPOSE: To evaluate the effect of the braced arm-to-thigh technique (BATT) (versus self-selected techniques) on three-dimensional trunk kinematics and spinal loads for three common activities of daily living (ADLs) simulated in the laboratory: weeding (gardening), reaching for an object in a low cupboard, and car egress using the two-legs out technique. METHODS: Ten young healthy males performed each task using a self-selected technique, and then using the BATT. The pulling action of weeding was simulated using a magnet placed on a steel plate. Cupboard and car egress tasks were simulated using custom apparatus representing the dimensions of a kitchen cabinet and a medium-sized Australian car, respectively. Three-dimensional trunk kinematics and L4/L5 spinal loads were estimated using the Lifting Full-Body OpenSim model and compared between techniques. Paired t-tests were used to compare peak values between methods (self-selected vs BATT). RESULTS: The BATT significantly reduced peak extension moments (13-51%), and both compression (27-45%) and shear forces (31-62%) at L4/L5, compared to self-selected techniques for all three tasks (p < 0.05). Lateral bending angles increased with the BATT for weeding and cupboard tasks, but these changes were expected as the BATT inherently introduces asymmetric trunk motion. CONCLUSION: The BATT substantially reduced L4/L5 extension moments, and L4/L5 compression and shear forces, compared to self-selected methods, for three ADLs, in a small cohort of ten young healthy males without prior history of back pain. These study findings can be used to inform safe procedures for these three ADLs, as the results are considered representative of a mature population.

Prediction of Knee Joint Contact Forces From External Measures Using Principal Component Prediction and Reconstruction.

Abnormal loading of the knee joint contributes to the pathogenesis of knee osteoarthritis. Gait retraining is a noninvasive intervention that aims to reduce knee loads by providing audible, visual, or haptic feedback of gait parameters. The computational expense of joint contact force prediction has limited real-time feedback to surrogate measures of the contact force, such as the knee adduction moment. We developed a method to predict knee joint contact forces using motion analysis and a statistical regression model that can be implemented in near real-time. Gait waveform variables were deconstructed using principal component analysis, and a linear regression was used to predict the principal component scores of the contact force waveforms. Knee joint contact force waveforms were reconstructed using the predicted scores. We tested our method using a heterogenous population of asymptomatic controls and subjects with knee osteoarthritis. The reconstructed contact force waveforms had mean (SD) root mean square differences of 0.17 (0.05) bodyweight compared with the contact forces predicted by a musculoskeletal model. Our method successfully predicted subject-specific shape features of contact force waveforms and is a potentially powerful tool in biofeedback and clinical gait analysis.

Impact on short-lived climate forcers (SLCFs) from a realistic land-use change scenario via changes in biogenic emissions.

More than one quarter of natural forests have been cleared by humans to make way for other land-uses, with changes to forest cover projected to continue. The climate impact of land-use change (LUC) is dependent upon the relative strength of several biogeophysical and biogeochemical effects. In addition to affecting the surface albedo and exchanging carbon dioxide (CO2) and moisture with the atmosphere, vegetation emits biogenic volatile organic compounds (BVOCs), altering the formation of short-lived climate forcers (SLCFs) including aerosol, ozone (O3) and methane (CH4). Once emitted, BVOCs are rapidly oxidised by O3, and the hydroxyl (OH) and nitrate (NO3) radicals. These oxidation reactions yield secondary organic products which are implicated in the formation and growth of aerosol particles and are estimated to have a negative radiative effect on the climate (i.e. a cooling). These reactions also deplete OH, increasing the atmospheric lifetime of CH4, and directly affect concentrations of O3; the latter two being greenhouse gases which impose a positive radiative effect (i.e. a warming) on the climate. Our previous work assessing idealised deforestation scenarios found a positive radiative effect due to changes in SLCFs; however, since the radiative effects associated with changes to SLCFs result from a combination of non-linear processes it may not be appropriate to scale radiative effects from complete deforestation scenarios according to the deforestation extent. Here we combine a land-surface model, a chemical transport model, a global aerosol model, and a radiative transfer model to assess the net radiative effect of changes in SLCFs due to historical LUC between the years 1850 and 2000.

Predicting sagittal plane biomechanics that minimize the axial knee joint contact force during walking.

Both development and progression of knee osteoarthritis have been associated with the loading of the knee joint during walking. We are, therefore, interested in developing strategies for changing walking biomechanics to offload the knee joint without resorting to surgery. In this study, simulations of human walking were performed using a 2D bipedal forward dynamics model. A simulation generated by minimizing the metabolic cost of transport (CoT) resembled data measured from normal human walking. Three simulations targeted at minimizing the peak axial knee joint contact force instead of the CoT reduced the peak force by 12-25% and increased the CoT by 11-14%. The strategies used by the simulations were (1) reduction in gastrocnemius muscle force, (2) avoidance of knee flexion during stance, and (3) reduced stride length. Reduced gastrocnemius force resulted from a combination of changes in activation and changes in the gastrocnemius contractile component kinematics. The simulations that reduced the peak contact force avoided flexing the knee during stance when knee motion was unrestricted and adopted a shorter stride length when the simulated knee motion was penalized if it deviated from the measured human knee motion. A higher metabolic cost in an offloading gait would be detrimental for covering a long distance without fatigue but beneficial for exercise and weight loss. The predicted changes in the peak axial knee joint contact force from the simulations were consistent with estimates of the joint contact force in a human subject who emulated the predicted kinematics. The results demonstrate the potential of using muscle-actuated forward dynamics simulations to predict novel joint offloading interventions.

The DNA sequence and biological annotation of human chromosome 1.

The reference sequence for each human chromosome provides the framework for understanding genome function, variation and evolution. Here we report the finished sequence and biological annotation of human chromosome 1. Chromosome 1 is gene-dense, with 3,141 genes and 991 pseudogenes, and many coding sequences overlap. Rearrangements and mutations of chromosome 1 are prevalent in cancer and many other diseases. Patterns of sequence variation reveal signals of recent selection in specific genes that may contribute to human fitness, and also in regions where no function is evident. Fine-scale recombination occurs in hotspots of varying intensity along the sequence, and is enriched near genes. These and other studies of human biology and disease encoded within chromosome 1 are made possible with the highly accurate annotated sequence, as part of the completed set of chromosome sequences that comprise the reference human genome.

Investigating the effects of activation state and location on lower limb tissue stiffness.

Lower limb tissue stiffness is contingent on various factors, including location, tissue composition, loading rates, and the geometry of the indenting object. Previous studies demonstrated that tissue stiffness varies greatly between individuals and between locations on an individual. Additionally, some studies have shown that activation of underlying muscle tissue increases bulk soft tissue stiffness. Yet, few studies have simultaneously considered both location and activation; this could be particularly important for measuring and predicting the function of devices such as prostheses and exoskeletons that interact with limbs at various locations during dynamic movement. In the present study, a custom handheld indentation device was used to explore changes in bulk leg tissue stiffness at rest and during isometric contractions. The indentation force-displacement curves were modelled using a Hertz model. At each level of activation (active/inactive), the shank had dramatically (∼150%) greater tissue stiffness than the thigh (p < 0.001). However, results suggested location independence for stiffness ratio (active/inactive, p = 0.42); for either location, stiffness was approximately 2x greater for active vs inactive muscle. These results should be considered during the development of biomechanical models to simulate human tissue indentation stiffness across a range of activation states and locations.

DNA sequence and analysis of human chromosome 9.

Chromosome 9 is highly structurally polymorphic. It contains the largest autosomal block of heterochromatin, which is heteromorphic in 6-8% of humans, whereas pericentric inversions occur in more than 1% of the population. The finished euchromatic sequence of chromosome 9 comprises 109,044,351 base pairs and represents >99.6% of the region. Analysis of the sequence reveals many intra- and interchromosomal duplications, including segmental duplications adjacent to both the centromere and the large heterochromatic block. We have annotated 1,149 genes, including genes implicated in male-to-female sex reversal, cancer and neurodegenerative disease, and 426 pseudogenes. The chromosome contains the largest interferon gene cluster in the human genome. There is also a region of exceptionally high gene and G + C content including genes paralogous to those in the major histocompatibility complex. We have also detected recently duplicated genes that exhibit different rates of sequence divergence, presumably reflecting natural selection.

The DNA sequence and analysis of human chromosome 6.

Chromosome 6 is a metacentric chromosome that constitutes about 6% of the human genome. The finished sequence comprises 166,880,988 base pairs, representing the largest chromosome sequenced so far. The entire sequence has been subjected to high-quality manual annotation, resulting in the evidence-supported identification of 1,557 genes and 633 pseudogenes. Here we report that at least 96% of the protein-coding genes have been identified, as assessed by multi-species comparative sequence analysis, and provide evidence for the presence of further, otherwise unsupported exons/genes. Among these are genes directly implicated in cancer, schizophrenia, autoimmunity and many other diseases. Chromosome 6 harbours the largest transfer RNA gene cluster in the genome; we show that this cluster co-localizes with a region of high transcriptional activity. Within the essential immune loci of the major histocompatibility complex, we find HLA-B to be the most polymorphic gene on chromosome 6 and in the human genome.

The DNA sequence and analysis of human chromosome 13.

Chromosome 13 is the largest acrocentric human chromosome. It carries genes involved in cancer including the breast cancer type 2 (BRCA2) and retinoblastoma (RB1) genes, is frequently rearranged in B-cell chronic lymphocytic leukaemia, and contains the DAOA locus associated with bipolar disorder and schizophrenia. We describe completion and analysis of 95.5 megabases (Mb) of sequence from chromosome 13, which contains 633 genes and 296 pseudogenes. We estimate that more than 95.4% of the protein-coding genes of this chromosome have been identified, on the basis of comparison with other vertebrate genome sequences. Additionally, 105 putative non-coding RNA genes were found. Chromosome 13 has one of the lowest gene densities (6.5 genes per Mb) among human chromosomes, and contains a central region of 38 Mb where the gene density drops to only 3.1 genes per Mb.

One-year Oxford knee scores should be used in preference to 6-month scores when assessing the outcome of total knee arthroplasty.

PURPOSE: The primary aim of this study was to assess whether there was a clinically significant difference in the mean Oxford knee score (OKS) between 6 and 12 months after total knee arthroplasty (TKA). The secondary aim was to identify variables associated with a clinically significant change in the OKS between 6 and 12 months. METHODS: A retrospective cohort study was undertaken using an established arthroplasty database of 1574 primary TKA procedures. Patient demographics, body mass index (BMI), comorbidities, OKS and EuroQoL 5-domain (EQ-5D) score were collected preoperatively and at 6 and 12 months postoperatively. A clinically significant change in the OKS was defined as 5 points or more. RESULTS: There was a 1.1-point increase in the OKS between 6 and 12 months postoperatively, which was statistically significant (95% confidence (CI) 0.8-1.3, p < 0.0001). There were 381 (24.2%) patients who had a clinically significant improvement in their OKS from 6 to 12 months. After adjusting for confounding, patients with a lower BMI (p = 0.028), without diabetes mellitus (p < 0.001), a better preoperative OKS (p < 0.001) or a worse 6-month OKS (p < 0.001) were more likely to have a clinically significant improvement. A 6-month OKS < 36 points was a reliable predictor of a clinically significant improvement in the 6-month to 12-month OKS (area under the curve 0.73, 95% CI 0.70-0.75, p < 0.001). CONCLUSION: Overall, there was no clinically significant change in the OKS from 6 to 12 months; however, a clinically significant improvement was observed in approximately a quarter of patients and was more likely in those scoring less than 36 points at 6 months. LEVEL OF EVIDENCE: retrospective diagnostic study, level III.

A braced arm-to-thigh (BATT) lifting technique reduces lumbar spine loads in healthy and low back pain participants.

Despite the common use of one-handed lifting techniques for activities of daily living, these techniques have received little attention in the biomechanics literature. The braced arm-to-thigh technique (BATT) is a one-handed lifting method in which the dominant hand picks up objects, while the free hand braces the trunk on the ipsilateral thigh. The aim of this study was to compare the BATT to two-handed or unsupported one-handed lifting techniques with loads of 2 and 10 kg, by evaluating trunk motion and spine loading at L4/L5. Twenty healthy participants (30-70 years old) matched in age and sex to 18 participants with low back pain were recruited to the study. A three-axis load cell secured to the distal anterior thigh measured the bracing forces applied by the hand. The OpenSim Lifting Full-Body model was used to estimate trunk kinematics and spinal loading at L4/L5. Linear mixed-effects models were developed to compare trunk angles and L4/L5 moments and forces between lifting techniques. Trunk flexion angles were significantly reduced for the BATT lift compared to one-handed and two-handed stoop lifts (9-20%). However, the BATT also increased asymmetric trunk kinematics and moments at L4/L5. The BATT produced significantly lower moments (28-38%), and compressive (25-32%) and antero-posterior shear (25-45%) forces at L4/L5, compared to unsupported lifting techniques. Bracing the hand on the thigh to support the trunk can substantially reduce low back loading during lifting tasks of 2 to 10 kg.

The effect of malalignment on proximal tibial strain in fixed-bearing unicompartmental knee arthroplasty: A comparison between metal-backed and all-polyethylene components using a validated finite element model.

OBJECTIVES: Elevated proximal tibial bone strain may cause unexplained pain, an important cause of unicompartmental knee arthroplasty (UKA) revision. This study investigates the effect of tibial component alignment in metal-backed (MB) and all-polyethylene (AP) fixed-bearing medial UKAs on bone strain, using an experimentally validated finite element model (FEM). METHODS: A previously experimentally validated FEM of a composite tibia implanted with a cemented fixed-bearing UKA (MB and AP) was used. Standard alignment (medial proximal tibial angle 90°, 6° posterior slope), coronal malalignment (3°, 5°, 10° varus; 3°, 5° valgus), and sagittal malalignment (0°, 3°, 6°, 9°, 12°) were analyzed. The primary outcome measure was the volume of compressively overstrained cancellous bone (VOCB) < -3000 µÎµ. The secondary outcome measure was maximum von Mises stress in cortical bone (MSCB) over a medial region of interest. RESULTS: Varus malalignment decreased VOCB but increased MSCB in both implants, more so in the AP implant. Varus malalignment of 10° reduced the VOCB by 10% and 3% in AP and MB implants but increased the MSCB by 14% and 13%, respectively. Valgus malalignment of 5° increased the VOCB by 8% and 4% in AP and MB implants, with reductions in MSCB of 7% and 10%, respectively. Sagittal malalignment displayed negligible effects. Well-aligned AP implants displayed greater VOCB than malaligned MB implants. CONCLUSION: All-polyethylene implants are more sensitive to coronal plane malalignments than MB implants are; varus malalignment reduced cancellous bone strain but increased anteromedial cortical bone stress. Sagittal plane malalignment has a negligible effect on bone strain.Cite this article: I. Danese, P. Pankaj, C. E. H. Scott. The effect of malalignment on proximal tibial strain in fixed-bearing unicompartmental knee arthroplasty: A comparison between metal-backed and all-polyethylene components using a validated finite element model. Bone Joint Res 2019;8:55-64. DOI: 10.1302/2046-3758.82.BJR-2018-0186.R2.

'Worse than death' and waiting for a joint arthroplasty.

AIMS: The EuroQol five-dimension (EQ-5D) questionnaire is a widely used multiattribute general health questionnaire where an EQ-5D < 0 defines a state 'worse than death' (WTD). The aim of this study was to determine the proportion of patients awaiting total hip arthroplasty (THA) or total knee arthroplasty (TKA) in a health state WTD and to identify associations with this state. Secondary aims were to examine the effect of WTD status on one-year outcomes. PATIENTS AND METHODS: A cross-sectional analysis of 2073 patients undergoing 2073 THAs (mean age 67.4 years (sd 11.6; 14 to 95); mean body mass index (BMI) 28.5 kg/m2 (sd 5.7; 15 to 72); 1253 female (60%)) and 2168 patients undergoing 2168 TKAs (mean age 69.3 years (sd 9.6; 22 to 91); BMI 30.8 kg/m2 (sd 5.8; 13 to 57); 1244 female (57%)) were recorded. Univariate analysis was used to identify variables associated with an EQ-5D score < 0: age, BMI, sex, deprivation quintile, comorbidities, and joint-specific function measured using the Oxford Hip Score (OHS) or Oxford Knee Score (OKS). Multivariate logistic regression was performed. EQ-5D and OHS/OKS were repeated one year following surgery in 1555 THAs and 1700 TKAs. RESULTS: Preoperatively, 391 THA patients (19%) and 263 TKA patients (12%) were WTD. Multivariate analysis identified preoperative OHS, deprivation, and chronic obstructive pulmonary disease in THA, and OKS, peripheral arterial disease, and inflammatory arthropathy in TKA as independently associated with WTD status (p < 0.05). One year following arthroplasty EQ-5D scores improved significantly (p < 0.001) and WTD rates reduced to 35 (2%) following THA and 53 (3%) following TKA. Patients who were WTD preoperatively achieved significantly (p < 0.001) worse joint-specific Oxford scores and satisfaction rates one year following joint arthroplasty, compared with those not WTD preoperatively. CONCLUSION: In total, 19% of patients awaiting THA and 12% awaiting TKA for degenerative joint disease are in a health state WTD. Although specific comorbidities contribute to this, hip- or knee-specific function, mainly pain, appear key determinants and can be reliably reversed with an arthroplasty. Cite this article: Bone Joint J 2019;101-B:941-950.

Can altered neuromuscular coordination restore soft tissue loading patterns in anterior cruciate ligament and menisci deficient knees during walking?

Injuries to the anterior cruciate ligament (ACL) and menisci commonly lead to early onset osteoarthritis. Treatments that can restore normative cartilage loading patterns may mitigate the risk of osteoarthritis, though it is unclear whether such a goal is achievable through conservative rehabilitation. We used musculoskeletal simulation to predict cartilage and ligament loading patterns during walking in intact, ACL deficient, menisci deficient, and ACL-menisci deficient knees. Stochastic simulations with varying coordination strategies were then used to test whether neuromuscular control could be modulated to restore normative knee mechanics in the pathologic conditions. During early stance, a 3 mm increase in anterior tibial translation was predicted in the ACL deficient knee. Mean cartilage contact pressure increased by 18% and 24% on the medial and lateral plateaus, respectively, in the menisci deficient knee. Variations in neuromuscular coordination were insufficient to restore normative cartilage contact patterns in either the ACL or menisci deficient knees. Elevated cartilage contact pressures in the pathologic knees were observed in regions where cartilage wear patterns have previously been reported. These results suggest that altered cartilage tissue loading during gait may contribute to region-specific degeneration patterns, and that varying neuromuscular coordination in isolation is unlikely to restore normative knee mechanics.

Development of a novel MATLAB-based framework for implementing mechanical joint stability constraints within OpenSim musculoskeletal models.

The Static Optimization (SO) solver in OpenSim estimates muscle activations and forces that only equilibrate applied moments. In this study, SO was enhanced through an open-access MATLAB interface, where calculated muscle activations can additionally satisfy crucial mechanical stability requirements. This Stability-Constrained SO (SCSO) is applicable to many OpenSim models and can potentially produce more biofidelic results than SO alone, especially when antagonistic muscle co-contraction is required to stabilize body joints. This hypothesis was tested using existing models and experimental data in the literature. Muscle activations were calculated by SO and SCSO for a spine model during two series of static trials (i.e. simulation 1 and 2), and also for a lower limb model (supplementary material 2). In simulation 1, symmetric and asymmetric flexion postures were compared, while in simulation 2, various external load heights were compared, where increases in load height did not change the external lumbar flexion moment, but necessitated higher EMG activations. During the tasks in simulation 1, the predicted muscle activations by SCSO demonstrated less average deviation from the EMG data (6.8% -7.5%) compared to those from SO (10.2%). In simulation 2, SO predicts constant muscle activations and forces, while SCSO predicts increases in the average activations of back and abdominal muscles that better match experimental data. Although the SCSO results are sensitive to some parameters (e.g. musculotendon stiffness), when considering the strategy of the central nervous system in distributing muscle forces and in activating antagonistic muscles, the assigned activations by SCSO are more biofidelic than SO.

Contributions of muscles and external forces to medial knee load reduction due to osteoarthritis braces.

BACKGROUND: Braces for medial knee osteoarthritis can reduce medial joint loads through a combination of three mechanisms: application of an external brace abduction moment, alteration of gait dynamics, and reduced activation of antagonistic muscles. Although the effect of knee bracing has been reported independently for each of these parameters, no previous study has quantified their relative contributions to reducing medial knee loads. METHODS: In this study, we used a detailed musculoskeletal model to investigate immediate changes in medial and lateral loads caused by two different knee braces: OA Assist and OA Adjuster 3 (DJO Global). Seventeen osteoarthritis subjects and eighteen healthy controls performed overground gait trials in unbraced and braced conditions. RESULTS: Across all subjects, bracing reduced medial loads by 0.1 to 0.3 times bodyweight (BW), or roughly 10%, and increased lateral loads by 0.03 to 0.2 BW. Changes in gait kinematics due to bracing were subtle, and had little effect on medial and lateral joint loads. The knee adduction moment was unaltered unless the brace moment was included in its computation. Only one muscle, biceps femoris, showed a significant change in EMG with bracing, but this did not contribute to altered peak medial contact loads. CONCLUSIONS: Knee braces reduced medial tibiofemoral loads primarily by applying a direct, and substantial, abduction moment to each subject's knee. To further enhance brace effectiveness, future brace designs should seek to enhance the magnitude of this unloader moment, and possibly exploit additional kinematic or neuromuscular gait modifications.