A critical appreciation of pathway analysis in atherosclerotic disease. Cellular phenotypic plasticity as an illustrative example

The rapid advancements in genome-scale (omics) techniques has created significant opportunities to investigate complex disease mechanisms in tissues and cells. Nevertheless, interpreting -omics data can be challenging, and pathway enrichment analysis is a frequently used method to identify candidate molecular pathways that drive gene expression changes. With a growing number of -omics studies dedicated to atherosclerosis, there has been a significant increase in studies and hypotheses relying on enrichment analysis. This brief review discusses the benefits and limitations of pathway enrichment analysis within atherosclerosis research. We highlight the challenges of identifying complex biological processes, such as cell phenotypic switching, within -omics data. Additionally, we emphasize the need for more comprehensive and curated gene sets that reflect the biological complexity of atherosclerosis. Pathway enrichment analysis is a valuable tool for gaining insights into the molecular mechanisms of atherosclerosis. Nevertheless, it is crucial to remain aware of the intrinsic limitations of this approach. By addressing these weaknesses, enrichment analysis in atherosclerosis can lead to breakthroughs in identifying the mechanisms of disease progresses, the identification of key driver genes, and consequently, advance personalized patient care.

Musculoskeletal modeling-based definition of load cases and worst-case fracture orientation for the design of clavicle fixation plates

Mechanical performance of clavicle fracture fixation plates is often evaluated using finite element (FE) analysis. Typically, these studies use simplified loading conditions and assume a transversal fracture orientation. However, the loading conditions and fracture orientation influence how the fracture site and thus fixation plate is loaded. In this study, a musculoskeletal model that included the clavicle muscles and scapulohumeral rhythm was defined based on previously published models. The standard OpenSim workflow (inverse kinematics, inverse dynamics, static optimization, and joint reaction analysis) was used to calculate muscle and joint contact forces based on 3D marker data collected in three subjects during seven activities of daily living (ADL). These loading conditions were then applied to a 3D clavicle model with three different fracture orientations and the mean resulting moments on both fragments were calculated to assess fracture stability. Magnitude of glenohumeral contact forces showed good agreement with instrumented shoulder prosthesis data, whereas simulated muscle activations were comparable to experimental EMG data. An oblique fracture orienting from superomedial to inferolateral was the least self-stabilizing. The loading to which the clavicle is exposed during ADL tasks is more complex than the simplified loading conditions typically used as boundary conditions in FE analyses of clavicle fracture fixation plates. Additionally, transversal fractures did not represent the least self-stabilizing fracture orientation, and thus calculated stresses in the plate could be underestimated. Therefore, more complex loading conditions and evaluation of a midshaft fracture running from superomedial to inferolateral is more relevant in FE analyses of midshaft clavicle fracture fixation plates.

Thin patient-specific clavicle fracture fixation plates can mechanically outperform commercial plates: An in silico approach

Current fixation plates used to operatively stabilize clavicular fractures are suboptimal, leading to reoperation rates up to 53%. Plate irritation, which can be caused by a poor geometric fit and plate thickness, has been found to be an important factor for reoperation. Moreover, muscle attachment sites (MAS) have to be detached occasionally. To improve current surgical clavicle fracture treatment with plate osteosynthesis, a change in plate design is required. The goal of this study was to design a patient-specific clavicle fracture fixation plate that includes geometrical optimization and stiffness determination. Its biomechanical performance has been compared with a commercial plate by examining the geometric fit, anatomical outline, stresses and interfragmentary motion using a finite element analysis with physiological loading and boundary conditions. Evaluation showed a better geometrical fit of the patient-specific plate as well as an improved fracture reduction. Displacements between fracture fragments were lower in case of the patient-specific plate, both when a fracture gap and no fracture gap were present. Results indicate a superior mechanical performance in terms of all investigated outcomes of the patient-specific plate compared to the commercial plate, while better aligning with the patient-specific geometry and without the need for MAS release. Due to the patient-specific geometry and reduced thickness, these fixation plates are expected to decrease the operation time, as intraoperative contouring will become irrelevant, and to decrease reoperation rates as implant irritation will be minimized.

Time-dependent in silico modelling of orthognathic surgery to support the design of biodegradable bone plates

Orthognathic surgery is performed to realign the jaws of a patient through several osteotomies. The state-of-the-art bone plates used to maintain the bone fragments in place are made of titanium. The presence of these non-degradable plates can have unwanted side effects on the long term (e.g. higher infection risk) if they are not removed. Using a biodegradable material such as magnesium may be a possible solution to this problem. However, biodegradation leads to a decrease of mechanical strength, therefore a time-dependent computational approach can help to evaluate the performance of such plates. In the present work, a computational framework has been developed to include biodegradation and bone healing algorithms in a finite element model. It includes bone plates and the mandible, which are submitted to masticatory loads during the early healing period (two months following the surgery). Two different bone plate designs with different stiffnesses have been tested. The stiff design exhibited good mechanical stability, with maximum Von Mises stress being less than 40% of the yield strength throughout the simulation. The flexible design shows high stresses when the bone healing has not started in the fracture gaps, indicating possible failure of the plate. However, this design leads to a higher bone healing quality after two months, as more cartilage is formed due to higher strains exerted in fracture gaps. We therefore conclude that in silico modelling can support tuning of the design parameters to ensure mechanical stability and while promoting bone healing.

Computationally Efficient Optimization Method to Quantify the Required Surgical Accuracy for a Ligament Balanced TKA

OBJECTIVE

This study proposes a computationally efficient method to quantify the effect of surgical inaccuracies on ligament strain in total knee arthroplasty (TKA). More specifically, this study describes a framework to determine the implant position and required surgical accuracy that results in a ligament balanced post-operative outcome with a probability of 90%.

METHODS

The response surface method is used to translate uncertainty in the implant position parameters to uncertainty in the ligament strain. The designed uncertainty quantification technique allows for an optimization with feasible computational cost towards the planned implant position and the tolerated surgical error for each of the twelve degrees of freedom of the implant position.

RESULTS

It is shown that the error does not allow for a ligament balanced TKA with a probability of 90% using preoperative planning. Six critical implant position parameters can be identified, namely AP translation, PD translation, VV rotation, IE rotation for the femoral component and PD translation, VV rotation for the tibial component.

CONCLUSION

We introduced an optimization process that allows for the computation of the required surgical accuracy for a ligament balanced postoperative outcome using preoperative planning with feasible computational cost.

SIGNIFICANCE

Towards the research society, the proposed method allows for a computationally efficient uncertainty quantification on a complex model. Towards surgical technique developers, six critical implant position parameters were identified, which should be the focus when refining surgical accuracy of TKA, leveraging better patient satisfaction.

Mechano-driven regeneration predicts response variations in large animal model based on scaffold implantation site and individual mechano-sensitivity

It is well founded that the mechanical environment may regulate bone regeneration in orthopedic applications. The purpose of this study is to investigate the mechanical contributions of the scaffold and the host to bone regeneration, in terms of subject specificity, implantation site and sensitivity to the mechanical environment. Using a computational approach to model mechano-driven regeneration, bone ingrowth in porous titanium scaffolds was simulated in the distal femur and proximal tibia of three goats and compared to experimental results. The results showed that bone ingrowth shifted from a homogeneous distribution pattern, when scaffolds were in contact with trabecular bone (max local ingrowth 12.47%), to a localized bone ingrowth when scaffolds were implanted in a diaphyseal location (max local ingrowth 20.64%). The bone formation dynamics revealed an apposition rate of 0.37±0.28%/day in the first three weeks after implantation, followed by limited increase in bone ingrowth until the end of the experiment (12 weeks). According to in vivo data, we identified one animal whose sensitivity to mechanical stimulation was higher than the other two. Moreover, we found that the stimulus initiating bone formation was consistently higher in the femur than in the tibia for all the individuals. Overall, the dependence of the osteogenic response on the host biomechanics means that, from a mechanical perspective, the regenerative potential depends on both the scaffold and the host environment. Therefore, this work provides insights on how the mechanical conditions of both the recipient and the scaffold contribute to meet patient and location-specific characteristics.

The effect of hip muscle weakness and femoral bony deformities on gait performance

BACKGROUND

Children with cerebral palsy (CP) present with a pathological gait pattern due to musculoskeletal impairments, such as muscle weakness and altered bony geometry. However, the effect of these impairments on gait performance is still unknown. Research aim:This study aimed to explore the effect of hip muscle weakness and femoral deformities on the gait performance of CP and typical developing (TD) subjects.

METHODS

6400 musculoskeletal models were created by weakening the hip extensors, abductors, adductors and flexors from 0% to 75 % and increasing the femoral anteversion angle (FAA) and neck shaft angle (NSA) from 20° to 60° and 120° to 160°, respectively. One TD and five CP gait patterns were imposed to each model and muscle forces were calculated. The effect of weakness and bony deformities on the capability gap (CG) at the hip, i.e. the lack in hip moment generating capacity to perform the gait pattern, was investigated using regression analysis.

RESULTS

The CG of apparent equinus, stiff knee gait, TD gait, jump gait and true equinus increased with 0.080, 0.038, 0.015, 0.023 and 0.005 Nm/kg per 10 percent hip abductor weakness increase, with 0.211, 0.130, 0.056, 0.045 and 0.011 Nm/kg per 10 degrees FAA increase and with 0.163, 0.080, 0.036, 0.043 and 0.011 Nm/kg per 10 degrees NSA increase, respectively. Combined weakness and bony deformities explained 96 %, 85 %, 82 %, 65 %, 40 % and 35 % of the variance in the CG of apparent equinus, TD gait, stiff knee gait, jump gait, true equinus and crouch gait, respectively.

SIGNIFICANCE

The results suggest that surgical correction of femoral deformities is more likely to be effective than strength training of hip muscles in enhancing CP gait performance. Jump gait, true equinus and especially crouch were more robust, while apparent equinus and stiff knee gait were limited by hip weakness and femoral deformities.

Hip Muscle Forces and Contact Loading During Squatting After Cam-Type FAI Surgery

BACKGROUND

The purpose of this study was to compare muscle forces and hip contact forces (HCFs) during squatting in patients with cam-type femoroacetabular impingement (cam-FAI) before and after hip corrective surgery and with healthy control participants.

METHODS

Ten symptomatic male patients with cam-FAI performed deep squatting preoperatively and at 2 years postoperatively. Patients were matched by age and body mass index to 10 male control participants. Full-body kinematics and kinetics were computed, and muscle forces and HCFs were estimated using a musculoskeletal model and static optimization. Normalized squat cycle (%SC) trials were compared using statistical nonparametric mapping (SnPM).

RESULTS

Postoperatively, patients with cam-FAI squatted down with higher anterior pelvic tilt, higher hip flexion, and greater hip extension moments than preoperatively. Preoperative patients demonstrated lower anterior pelvic tilt and lower hip flexion compared with the participants in the control group. Postoperative patients showed increased semimembranosus force compared with their preoperative values. Preoperative forces were lower than the control group for the adductor magnus, the psoas major, and the semimembranosus; however, the preoperative patients showed greater inferior gluteus maximus forces than the patients in the control group, whereas the postoperative patients did not differ from the control patients. Higher posterior, superior, and resultant HCF magnitudes were identified postoperatively in comparison with the preoperative values. Preoperative posterior HCF was lower than in the control group, whereas the postoperative posterior HCF did not differ from those in the control group.

CONCLUSIONS

Higher postoperative anterior pelvic tilt was associated with an indication of return to closer to normal pelvic motion, which resembled data from the control group. Lower preoperative anterior pelvic tilt was associated with muscle force imbalance, indicated by decreased semimembranosus and increased gluteus maximus forces. The overall increased postoperative muscle forces were associated with improved pelvic mobility and increased HCFs that were comparable with the control-group standards.

CLINICAL RELEVANCE

Muscle forces and HCFs may be indicative of postoperative joint health restoration and alleviated symptoms.

Botulinum toxin injections minimally affect modelled muscle forces during gait in children with cerebral palsy

BACKGROUND

Children with cerebral palsy (CP) present altered gait patterns and electromyography (EMG) activity compared to typically developing children. To temporarily reduce muscular activity and to correct the abnormal muscle force balance, Botulinum Toxin type A (BTX-A) injections are used.

RESEARCH QUESTION

What is the effect of BTX-A injections on dynamic muscle forces during gait, when calculated using an EMG-constrained approach?.

METHODS

Retrospective data of ten typically developing (TD) and fourteen children with spastic diplegic CP were used for musculoskeletal modeling and dynamic simulations of gait, before and after BTX-A treatment. Individual muscle forces were calculated using an EMG-constrained optimization, in which EMG of eight muscles was used as muscle excitation signal to constrain the muscle activation patterns. Paired t-tests were used to compare average modelled muscle forces in different phases of the gait cycle pre- and post-BTX-A, summarized in the muscle profile score. Two-sample t-tests were used to determine significant differences between TD and pre- and post-BTX-A modelled muscle forces.

RESULTS

For most muscles, the force was decreased in CP compared to TD children in all phases of the gait cycle, both before and after BTX-A treatment. Differences in muscle forces before and after BTX-A treatment were limited, with only few significant differences between pre- and post-BTX-A. Compared to a standard static optimization approach, imposing the EMG activity increased modelled muscle forces for most muscles.

SIGNIFICANCE

Our findings indicate that BTX-A treatment has a limited effect on the muscle balance in CP children. Besides that, the use of EMG-constrained optimization is recommended when studying muscle balance in children with CP.

Knee Load Distribution in Hip Osteoarthritis Patients After Total Hip Replacement

Reduced external knee adduction moments in the second half of stance after total hip replacement have been reported in hip osteoarthritis patients. This reduction is thought to shift the load from the medial to the lateral knee compartment and as such increase the risk for knee osteoarthritis. The knee adduction moment is a surrogate for the load distribution between the medial and lateral compartments of the knee and not a valid measure for the tibiofemoral contact forces which are the result of externally applied forces and muscle forces. The purpose of this study was to investigate whether the distribution of the tibiofemoral contact forces over the knee compartments in unilateral hip osteoarthritis patients 1 year after receiving a primary total hip replacement differs from healthy controls. Musculoskeletal modeling on gait was performed in OpenSim using the detailed knee model of Lerner et al. (2015) for 19 patients as well as for 15 healthy controls of similar age. Knee adduction moments were calculated by the inverse dynamics analysis, medial and lateral tibiofemoral contact forces with the joint reaction force analysis. Moments and contact forces of patients and controls were compared using Statistical Parametric Mapping two-sample t-tests. Knee adduction moments and medial tibiofemoral contact forces of both the ipsi- and contralateral leg were not significantly different compared to healthy controls. The contralateral leg showed 14% higher medial tibiofemoral contact forces compared to the ipsilateral (operated) leg during the second half of stance. During the first half of stance, the lateral tibiofemoral contact force of the contralateral leg was 39% lower and the ratio 32% lower compared to healthy controls. In contrast, during the second half of stance the forces were significantly higher (39 and 26%, respectively) compared to healthy controls. The higher ratio indicates a changed distribution whereas the increased lateral tibiofemoral contact forces indicate a higher lateral knee joint loading in the contralateral leg in OA patients after total hip replacement (THR). Musculoskeletal modeling using a detailed knee model can be useful to detect differences in the load distribution between the medial and lateral knee compartment which cannot be verified with the knee adduction moment.

Functional movement assessment by means of inertial sensor technology to discriminate between movement behaviour of healthy controls and persons with knee osteoarthritis

BACKGROUND

Apart from biomechanical alterations in movement patterns, it is known that movement limitations in persons with knee osteoarthritis (PwKOA) are related to an individual's perception and belief regarding pain and disability. To gain more insights into the functional movement behaviour of PwKOA in a clinical setting, inertial sensor technology can be applied. This study first aims to evaluate the ability of inertial sensors to discriminate between healthy controls (HC) and PwKOA. Secondly, this study aims to determine the relationship between movement behaviour, pain-related factors and disability scores.

METHODS

Twelve HC and 19 PwKOA were included. Five repetitions of six functional movement tasks (walking, forward lunge, sideward lunge, ascent and descent stairs, single leg squat and sit-to-stand) were simultaneously recorded by the inertial sensor system and a camera-based motion analysis system. Statistically significant differences in angular waveforms of the trunk, pelvis and lower limb joints between HC and PwKOA were determined using one-dimensional statistical parametric mapping (SPM1D). The Knee injury and Osteoarthritis Outcome Score and TAMPA scale for Kinesiophobia were used to evaluate the relationship between discriminating joint motion, pain-related factors and disability using spearman's correlation coefficients.

RESULTS

PwKOA had significantly less trunk rotation, internal pelvis rotation and knee flexion ROM during walking. Additionally, the reduced knee flexion (i.e. at the end of the stance phase and swing phase) was related to increased level of perceived pain. During the sideward lunge, PwKOA had significantly less knee flexion, ankle plantarflexion and hip abduction. This decreased hip abduction (i.e. during stance) was related to higher fear of movement. Finally, PwKOA had significantly less knee flexion during the forward lunge, single leg squat and during ascent and descent stairs. No significant correlations were observed with disability.

CONCLUSIONS

Inertial sensors were able to discriminate between movement characteristics of PwKOA and HC. Additionally, significant relationships were found between joint motion, perceived pain and fear of movement. Since inertial sensors can be used outside the laboratory setting, these results are promising as they indicate the ability to evaluate movement deviations. Further research is required to enable measurements of small movement deviations in clinically relevant tasks.

Discriminant validity of 3D joint kinematics and centre of mass displacement measured by inertial sensor technology during the unipodal stance task

Background

The unipodal stance task is a clinical task that quantifies postural stability and alignment of the lower limb joints, while weight bearing on one leg. As persons with knee osteoarthritis (KOA) have poor postural and knee joint stability, objective assessment of this task might be useful.

Objective

To investigate the discriminant validity of three-dimensional joint kinematics and centre of mass displacement (COM) between healthy controls and persons with knee KOA, during unipodal stance using inertial sensors. Additionally, the reliability, agreement and construct validity are assessed to determine the reproducibility and accuracy of the discriminating parameters.

Methods

Twenty healthy controls and 19 persons with unilateral severe KOA were included. Five repetitions of the unipodal stance task were simultaneously recorded by an inertial sensor system and a camera-based system (gold standard). Statistical significant differences in kinematic waveforms between healthy controls and persons with severe knee KOA were determined using one-dimensional statistical parametric mapping (SPM1D).

Results

Persons with severe knee KOA had more lateral trunk lean towards the contralateral leg, more hip flexion throughout the performance of the unipodal stance task, more pelvic obliquity and COM displacement towards the contralateral side. However, for the latter two parameters the minimum detectable change was greater than the difference between healthy controls and persons with severe knee KOA. The construct validity was good (coefficient of multiple correlation 0.75, 0.83 respectively) and the root mean squared error (RMSE) was low (RMSE <1.5°) for the discriminant parameters.

Conclusion

Inertial sensor based movement analysis can discriminate between healthy controls and persons with severe knee KOA for lateral trunk lean and hip flexion, but unfortunately not for the knee angles. Further research is required to improve the reproducibility and accuracy of the inertial sensor measurements before they can be used to assess differences in tasks with a small range of motion.

Physics-Based Simulations to Predict the Differential Effects of Motor Control and Musculoskeletal Deficits on Gait Dysfunction in Cerebral Palsy: A Retrospective Case Study

Physics-based simulations of walking have the theoretical potential to support clinical decision-making by predicting the functional outcome of treatments in terms of walking performance. Yet before using such simulations in clinical practice, their ability to identify the main treatment targets in specific patients needs to be demonstrated. In this study, we generated predictive simulations of walking with a medical imaging based neuro-musculoskeletal model of a child with cerebral palsy presenting crouch gait. We explored the influence of altered muscle-tendon properties, reduced neuromuscular control complexity, and spasticity on gait dysfunction in terms of joint kinematics, kinetics, muscle activity, and metabolic cost of transport. We modeled altered muscle-tendon properties by personalizing Hill-type muscle-tendon parameters based on data collected during functional movements, simpler neuromuscular control by reducing the number of independent muscle synergies, and spasticity through delayed muscle activity feedback from muscle force and force rate. Our simulations revealed that, in the presence of aberrant musculoskeletal geometries, altered muscle-tendon properties rather than reduced neuromuscular control complexity and spasticity were the primary cause of the crouch gait pattern observed for this child, which is in agreement with the clinical examination. These results suggest that muscle-tendon properties should be the primary target of interventions aiming to restore an upright gait pattern for this child. This suggestion is in line with the gait analysis following muscle-tendon property and bone deformity corrections. Future work should extend this single case analysis to more patients in order to validate the ability of our physics-based simulations to capture the gait patterns of individual patients pre- and post-treatment. Such validation would open the door for identifying targeted treatment strategies with the aim of designing optimized interventions for neuro-musculoskeletal disorders.

Subject-specific geometry affects acetabular contact pressure during gait more than subject-specific loading patterns

Finite element modeling (FEM) can predict hip cartilage contact mechanics. This study investigated how subject-specific boundary conditions and joint geometry affect acetabular cartilage contact mechanics using a multi-scale workflow. For two healthy subjects, musculoskeletal models calculated subject-specific hip kinematics and loading, which were used as boundary conditions for FEM. Cartilage contact mechanics were predicted using either generic or subject-specific FEM and boundary conditions. A subject-specific mesh resulted in a more lateral contact. Effects of subject-specific boundary conditions varied between both subjects. Results highlight the complex interplay between loading and kinematics and their effect on cartilage contact mechanics.

Comparison of modelling and tracking methods for analysing elbow and forearm kinematics

This study aimed to find an optimal measurement protocol of elbow and forearm kinematics using different modelling and tracking methods. Kinematic data of elbow flexion/extension and forearm pronation/supination was acquired using optical motion capture from 12 healthy male volunteers. Segment coordinate systems for humerus, forearm, radius, ulna, and hand were defined. Different tracking methods, using anatomical markers or rigid or point maker clusters, were used to compute the three-dimensional rotations. Marker placement errors were assessed to evaluate the rigid body assumption. Multiple comparisons demonstrated statistical differences between tracking methods: compared to using only anatomical markers, tracking using clusters reduced the estimated range of pronation/supination by 14.9%–43.2%, while it estimated increased flexion/extension by 5.3%–9.1%. The study suggests using only anatomical markers exerts the optimal estimation of elbow and forearm kinematics. Modelling using the coordinate systems of the humerus and forearm and of the humerus and ulna, respectively, demonstrated good consistency with literature and are correspondingly regarded as the most appropriate approach for measuring pronation/supination and flexion/extension. The results are valuable in establishing a measurement protocol for analysing elbow and forearm kinematics, avoiding confusions and misinterpretations in communicating results from different methodologies.

The influence of knee joint geometry and alignment on the tibiofemoral load distribution: A computational study

BACKGROUND

Deviations in knee joint geometry and alignment were previously related to an increased risk for knee OA. These were hypothesized to influence the load distribution over the articular cartilage. Therefore, this study evaluated the effect of altered knee joint geometry and alignment in the coronal and transverse plane on the medial-lateral load distribution and ligament strain using a musculoskeletal modeling approach.

METHODS

Joint kinematics during gait were measured in 15 healthy adults. Using different musculoskeletal models with altered geometry of the tibia plateau or knee joint malalignment in the coronal and transverse plane, the resulting muscle, ligament and contact forces were calculated. Next, the distribution of the load over the medial and lateral condyle was analyzed and compared to the reference loading distribution, with neutral geometry and alignment, using repeated-measures ANOVA and individual t-tests, with a Bonferroni-corrected alpha level.

RESULTS

Coronal plane malalignment significantly affected the load distribution. Small changes in coronal tibial slope had less pronounced effects on the load distribution, but increased ligament strains. Transverse plane malalignment only minimally affected the load distribution.

CONCLUSION

Coronal plane knee malalignment affected knee loading, with increased varus alignment resulting in increased medial loading. This confirms a causal relation between coronal malalignment and increased medial compartment loading and suggests a potential role of aberrant coronal plane alignment on OA initiation. Altered coronal tibial slope induced increased ligament strains, potentially contributing to a cascade of knee laxity and subsequently more extreme knee malalignment.

SimCP: A Simulation Platform to Predict Gait Performance Following Orthopedic Intervention in Children With Cerebral Palsy

Gait deficits in cerebral palsy (CP) are often treated with a single-event multi-level surgery (SEMLS). Selecting the treatment options (combination of bony and soft tissue corrections) for a specific patient is a complex endeavor and very often treatment outcome is not satisfying. A deterioration in 22.8% of the parameters describing gait performance has been reported and there is need for additional surgery in 11% of the patients. Computational simulations based on musculoskeletal models that allow clinicians to test the effects of different treatment options before surgery have the potential to drastically improve treatment outcome. However, to date, no such simulation and modeling method is available. Two important challenges are the development of methods to include patient-specific neuromechanical impairments into the models and to simulate the effect of different surgical procedures on post-operative gait performance. Therefore, we developed the SimCP framework that allows the evaluation of the effect of different simulated surgeries on gait performance of a specific patient and includes a graphical user interface (GUI) that enables performing virtual surgery on the models. We demonstrated the potential of our framework for two case studies. Models reflecting the patient-specific musculoskeletal geometry and muscle properties are generated based solely on data collected before the treatment. The patient's motor control is described based on muscle synergies derived from pre-operative EMG. The GUI is then used to modify the musculoskeletal properties according to the surgical plan. Since SEMLS does not affect motor control, the same motor control model is used to define gait performance pre- and post-operative. We use the capability gap (CG), i.e., the difference between the joint moments needed to perform healthy walking and the joint moments the personalized model can generate, to quantify gait performance. In both cases, the CG was smaller post- then pre-operative and this was in accordance with the measured change in gait kinematics after treatment.

Selective dorsal rhizotomy improves muscle forces during walking in children with spastic cerebral palsy

BACKGROUND

Selective dorsal rhizotomy aims to reduce spasticity in children with cerebral palsy. Early investigations indicated postoperative weakness, whereas more recent studies showed that selective dorsal rhizotomy either does not change or improves muscle strength. All previous studies assessed muscle strength in a static position, which did not represent the walking situation. The aim of this study was to analyze the influence of selective dorsal rhizotomy on muscle forces during gait.

METHODS

Motion capture data of 25 children with spastic cerebral palsy and 10 typically developing participants were collected. A musculoskeletal OpenSim model was used to calculate joint kinematics, joint kinetics and muscle forces during gait. Static optimization and an electromyography-informed approach to calculate muscle forces were compared. A Muscle-Force-Profile was introduced and used to compare the muscle forces during walking before and after a selective dorsal rhizotomy.

FINDINGS

Independent of the approach used (electromyography-informed versus static optimization), selective dorsal rhizotomy significantly normalized forces in spastic muscles during walking and did not reduce the contribution of non-spastic muscles.

INTERPRETATION

This study showed that selective dorsal rhizotomy improves dynamic muscle forces in children with cerebral palsy and leads to less gait pathology, as shown in the improvement in joint kinematics and joint kinetics. Individual muscle force analyses using the Muscle-Force-Profile extend standard joint kinematics and joint moment analyses, which might improve clinical-decision making in children with cerebral palsy in the future. The reference data of our participants and MATLAB code for the Muscle-Force-Profile are publicly available on simtk.org/projects/muscleprofile.

Topographical Variation of Human Femoral Articular Cartilage Thickness, T1rho and T2 Relaxation Times Is Related to Local Loading during Walking

OBJECTIVE

Early detection of degenerative changes in the cartilage matrix composition is essential for evaluating early interventions that slow down osteoarthritis (OA) initiation. T1rho and T2 relaxation times were found to be effective for detecting early changes in proteoglycan and collagen content. To use these magnetic resonance imaging (MRI) methods, it is important to document the topographical variation in cartilage thickness, T1rho and T2 relaxation times in a healthy population. As OA is partially mechanically driven, the relation between these MRI-based parameters and localized mechanical loading during walking was investigated.

DESIGN

MR images were acquired in 14 healthy adults and cartilage thickness and T1rho and T2 relaxation times were determined. Experimental gait data was collected and processed using musculoskeletal modeling to identify weight-bearing zones and estimate the contact force impulse during gait. Variation of the cartilage properties (i.e., thickness, T1rho, and T2) over the femoral cartilage was analyzed and compared between the weight-bearing and non-weight-bearing zone of the medial and lateral condyle as well as the trochlea.

RESULTS

Medial condyle cartilage thickness was correlated to the contact force impulse ( r = 0.78). Lower T1rho, indicating increased proteoglycan content, was found in the medial weight-bearing zone. T2 was higher in all weight-bearing zones compared with the non-weight-bearing zones, indicating lower relative collagen content.

CONCLUSIONS

The current results suggest that medial condyle cartilage is adapted as a long-term protective response to localized loading during a frequently performed task and that the weight-bearing zone of the medial condyle has superior weight bearing capacities compared with the non-weight-bearing zones.

Longitudinal joint loading in patients before and up to one year after unilateral total hip arthroplasty

Abnormal kinematics and kinetics have been reported in hip osteoarthritis (OA) patients before and after total hip arthroplasty (THA). These changes can affect the loading of the ipsilateral hip, as well as the contralateral hip and knee joint. As it is not clear how hip and knee loading evolves in THA patients during the first year after surgery, the goal of this study is to define how joint loading changes in patients before and at three evaluation times after THA surgery. Musculoskeletal modelling in combination with gait analysis data was used to calculate hip and knee contact forces in 14 patients before and 3-, 6- and 12-months after unilateral THA, as well as in 18 healthy controls. Results showed that bilateral hip and knee loading were decreased compared to controls, both before and after THA surgery. Loading symmetry was altered compared to controls at 3-months post-surgery for the hip and at all evaluation times, except for 6-months post-surgery, for the knee, with ipsilateral joint loading decreased compared to the contralateral side. To conclude, 12-months after THA joint loading was not normalized, with both hip and knee loading in patients decreased compared to controls. Therefore, no overloading of the ipsi- or contralateral hip and knee joint was found before and up to one year after unilateral THA.

The morphological and molecular mechanisms of epithelial/endothelial-to-mesenchymal transition and its involvement in atherosclerosis

Cell transdifferentiation occurs during cardiovascular development or remodeling either as a pathologic feature in the progression of disease or as a response to injury. Endothelial-to-Mesenchymal Transition (EndMT) is a process that is classified as a specialized form of Epithelial-to-Mesenchymal Transition (EMT), in which epithelial cells lose their epithelial characteristics and gain a mesenchymal phenotype. During transdifferentiation, cells lose both cell-cell contacts and their attachment to the basement membrane. Subsequently, the shape of the cells changes from a cuboidal to an elongated shape. A rearrangement of actin filaments facilitates the cells to become motile and prime their migration into the underlying tissue. EMT is a key process during embryonic development, wound healing and tissue regeneration, but has also been implicated in pathophysiological processes, such organ fibrosis and tumor metastases. EndMT has been associated with additional pathophysiological processes in cardiovascular related diseases, including atherosclerosis. Recent studies prove a significant role for EndMT in the progression and destabilization of atherosclerotic plaques, as a consequence of EndMT-derived fibroblast infiltration and the increased secretion of matrix metalloproteinase respectively. In this review we will discuss the essential molecular and morphological mechanisms of EMT and EndMT, along with their common denominators and key differences. Finally, we will discuss the role of EMT/EndMT in developmental and pathophysiological processes, focusing on the potential role of EndMT in atherosclerosis in more depth.

Medial knee loading is altered in subjects with early osteoarthritis during gait but not during step-up-and-over task

This study evaluates knee joint loading during gait and step-up-and-over tasks in control subjects, subjects with early knee OA and those with established knee OA. Thirty-seven subjects with varying degrees of medial compartment knee OA severity (eighteen with early OA and sixteen with established OA), and nineteen healthy controls performed gait and step-up-and-over tasks. Knee joint moments, contact forces (KCF), the magnitude of contact pressures and center of pressure (CoP) location were analyzed for the three groups for both activities using a multi-body knee model with articular cartilage contact, 14 ligaments, and six degrees of freedom tibiofemoral and patellofemoral joints. During gait, the first peak of the medial KCF was significantly higher for patients with early knee OA (p = 0.048) and established knee OA (p = 0.001) compared to control subjects. Furthermore, the medial contact pressure magnitudes and CoP location were significantly different in both groups of patients compared to controls. Knee rotation moments (KRMs) and external rotation angles were significantly higher during early stance in both patient groups (p < 0.0001) compared to controls. During step-up-and-over, there was a high variability between the participants and no significant differences in KCF were observed between the groups. Knee joint loading and kinematics were found to be altered in patients with early knee OA only during gait. This is an indication that an excessive medial KCF and altered loading location, observed in these patients, is a contributor to early progression of knee OA.

Gait alterations can reduce the risk of edge loading

Following metal-on-metal hip arthroplasty, edge loading (i.e., loading near the edge of a prosthesis cup) can increase wear and lead to early revision. The position and coverage angle of the prosthesis cup influence the risk of edge loading. This study investigates the effect of altered gait patterns, more specific hip, and pelvis kinematics, on the orientation of hip contact force and the consequent risk of antero-superior edge loading using muscle driven simulations of gait. With a cup orientation of 25° anteversion and 50° inclination and a coverage angle of 168°, many gait patterns presented risk of edge loading. Specifically at terminal double support, 189 out of 405 gait patterns indicated a risk of edge loading. At this time instant, the high hip contact forces and the proximity of the hip contact force to the edge of the cup indicated the likelihood of the occurrence of edge loading. Although the cup position contributed most to edge loading, altering kinematics considerably influenced the risk of edge loading. Increased hip abduction, resulting in decreasing hip contact force magnitude, and decreased hip extension, resulting in decreased risk on edge loading, are gait strategies that could prevent edge loading. © 2015 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 34:1069-1076, 2016.

Subject-specific geometrical detail rather than cost function formulation affects hip loading calculation

This study assessed the relative importance of introducing an increasing level of medical image-based subject-specific detail in bone and muscle geometry in the musculoskeletal model, on calculated hip contact forces during gait. These forces were compared to introducing minimization of hip contact forces in the optimization criterion. With an increasing level of subject-specific detail, specifically MRI-based geometry and wrapping surfaces representing the hip capsule, hip contact forces decreased and were more comparable to contact forces measured using instrumented prostheses (average difference of 0.69 BW at the first peak compared to 1.04 BW for the generic model). Inclusion of subject-specific wrapping surfaces in the model had a greater effect than altering the cost function definition.

The role of altered proximal femoral geometry in impaired pelvis stability and hip control during CP gait: A simulation study

Children with cerebral palsy (CP) often present aberrant hip geometry, more specifically increased femoral anteversion and neck-shaft angle. Furthermore, altered gait patterns are present within this population. This study analyzed the effect of aberrant femoral geometry, as present in subjects with CP, on the ability of muscles to control hip and knee joint kinematics. Given the specific gait deficits observed during crouch gait, increased ability to abduct, externally rotate the hip and extend the knee and hip were denoted as beneficial effects. We ran dynamic simulations of CP and normal gait using two musculoskeletal models, one reflecting normal femoral geometry and one reflecting proximal femoral deformities. The results show that the combination of aberrant bone geometry and CP-specific gait characteristics beneficially increased the ability of gluteus medius and maximus to extend the hip and knee. In contrast, the potentials of the hamstrings to extend the hip decreased whereas the potentials to flex the knee increased. These changes closely followed the observed changes in the muscle moment arm lengths. In conclusion, this study emphasizes the concomitant effect of the presence of proximal femoral deformity and CP gait characteristics on the muscle control of hip and knee joint kinematics during single stance. Not accounting for subject-specific geometry will affect the calculated muscles' potential during gait. Therefore, the use of generic models to assess muscle function in the presence of femoral deformity and CP gait should be treated with caution.

Does surgical approach or prosthesis type affect hip joint loading one year after surgery?

Several approaches may be used for hip replacement surgery either in combination with conventional total hip arthroplasty (THA) or resurfacing hip arthroplasty (RHA). This study investigates the differences in hip loading during gait one year or more after surgery in three cohorts presenting different surgical procedures, more specific RHA placed using the direct lateral (RHA-DLA, n=8) and posterolateral (RHA-PLA, n=14) approach as well as THA placed using the direct anterior (THA-DAA, n=12) approach. For the DAA and control subjects, hip loading was also evaluated during stair ascent and descent to evaluate whether these motions can better discriminate between patients and controls compared to gait. Musculoskeletal modelling in OpenSim was used to calculate in vivo joint loading. Results showed that for all operated patients, regardless the surgical procedure, hip loading was decreased compared to control subjects, while no differences were found between patient groups. This indicates that THA via DAA results in similar hip loading as a RHA via DLA or PLA. Stair climbing did not result in more distinct differences in hip contact force magnitude between patients and controls, although differences in orientation were more distinct. However, patients after hip surgery did adjust their motion pattern to decrease the magnitude of loading on the hip joint compared to control subjects.

Evaluation of stair motion contributes to new insights into hip osteoarthritis-related motion pathomechanics

Stair motion in the presence of hip osteoarthritis (OA) has received less attention than level walking. Its more strenuous aspect may shed the light on different locomotor strategies when compared to walking. We, therefore, aimed to define stair motion features associated to hip OA and to evaluate whether these specific features would differ from level walking and better characterize the hip pathological condition. Principal component and linear discriminant analyses were, respectively, used as data reduction and classification techniques. Our study highlighted that most of stair motion features associated to hip OA were similar to the ones of walking. Stair descent presented with the lowest misclassification error rate, ranging from 12% to 19% (estimated by cross-validation). But, features that may be considered as a mechanism to reduce demand on the hip abductors were found to be more important in the stair ascent condition. This was reflected by both, greater importance in the classification rule and variance compared with walking, that is, decreased hip internal rotation moment at mid-stance (72.50% vs. 57.63%) and increased trunk lateroflexion toward affected side (56.43% vs. 29.37%). This study emphasized the importance of investigating stair motion in hip osteoarthritic population by highlighting specific locomotor strategies.

Biomechanical gait features associated with hip osteoarthritis: Towards a better definition of clinical hallmarks

Critical appraisal of the literature highlights that the discriminative power of gait-related features in patients with hip osteoarthritis (OA) has not been fully explored. We aimed to reduce the number of gait-related features and define the most discriminative ones comparing the three-dimensional gait analysis of 20 patients with hip osteoarthritis (OA) with those of 17 healthy peers. First, principal component analysis was used to reduce the high-dimensional gait data into a reduced set of interpretable variables for further analysis, including tests for group differences. These differences were indicative for the selection of the top 10 variables to be included into linear discriminant analysis models (LDA). Our findings demonstrated the successful data reduction of hip osteoarthritic-related gait features with a high discriminatory power. The combination of the top variables into LDA models clearly separated groups, with a maximum misclassification error rate of 19%, estimated by cross-validation. Decreased hip/knee extension, hip flexion and internal rotation moment were gait features with the highest discriminatory power. This study listed the most clinically relevant gait features characteristics of hip OA. Moreover, it will help clinicians and physiotherapists understand the movement pathomechanics related to hip OA useful in the management and design of rehabilitation intervention.

Loading of Hip Measured by Hip Contact Forces at Different Speeds of Walking and Running

Exercise plays a pivotal role in maximizing peak bone mass in adulthood and maintaining it through aging, by imposing mechanical loading on the bone that can trigger bone mineralization and growth. The optimal type and intensity of exercise that best enhances bone strength remains, however, poorly characterized, partly because the exact peak loading of the bone produced by the diverse types of exercises is not known. By means of integrated motion capture as an input to dynamic simulations, contact forces acting on the hip of 20 young healthy adults were calculated during walking and running at different speeds. During walking, hip contact forces (HCFs) have a two-peak profile whereby the first peak increases from 4.22 body weight (BW) to 5.41 BW and the second from 4.37 BW to 5.74 BW, by increasing speed from 3 to 6 km/h. During running, there is only one peak HCF that increases from 7.49 BW to 10.01 BW, by increasing speed from 6 to 12 km/h. Speed related profiles of peak HCFs and ground reaction forces (GRFs) reveal a different progression of the two peaks during walking. Speed has a stronger impact on peak HCFs rather than on peak GRFs during walking and running, suggesting an increasing influence of muscle activity on peak HCF with increased speed. Moreover, results show that the first peak of HCF during walking can be predicted best by hip adduction moment, and the second peak of HCF by hip extension moment. During running, peak HCF can be best predicted by hip adduction moment. The present study contributes hereby to a better understanding of musculoskeletal loading during walking and running in a wide range of speeds, offering valuable information to clinicians and scientists exploring bone loading as a possible nonpharmacological osteogenic stimulus. © 2015 American Society for Bone and Mineral Research.

Gait alterations to effectively reduce hip contact forces

Patients with hip pathology present alterations in gait which have an effect on joint moments and loading. In knee osteoarthritic patients, the relation between medial knee contact forces and the knee adduction moment are currently being exploited to define gait retraining strategies to effectively reduce pain and disease progression. However, the relation between hip contact forces and joint moments has not been clearly established. Therefore, this study aims to investigate the effect of changes in hip and pelvis kinematics during gait on internal hip moments and contact forces which is calculated using muscle driven simulations. The results showed that frontal plane kinetics have the largest effect on hip contact forces. Given the high correlation between the change in hip adduction moment and contact force at initial stance (R(2)  = 0.87), this parameter can be used to alter kinematics and predict changes in contact force. At terminal stance the hip adduction and flexion moment can be used to predict changes in contact force (R(2)  = 0.76). Therefore, gait training that focuses on decreasing hip adduction moments, a wide base gait pattern, has the largest potential to reduce hip contact forces.

Muscle optimization techniques impact the magnitude of calculated hip joint contact forces

In musculoskeletal modelling, several optimization techniques are used to calculate muscle forces, which strongly influence resultant hip contact forces (HCF). The goal of this study was to calculate muscle forces using four different optimization techniques, i.e., two different static optimization techniques, computed muscle control (CMC) and the physiological inverse approach (PIA). We investigated their subsequent effects on HCFs during gait and sit to stand and found that at the first peak in gait at 15-20% of the gait cycle, CMC calculated the highest HCFs (median 3.9 times peak GRF (pGRF)). When comparing calculated HCFs to experimental HCFs reported in literature, the former were up to 238% larger. Both static optimization techniques produced lower HCFs (median 3.0 and 3.1 pGRF), while PIA included muscle dynamics without an excessive increase in HCF (median 3.2 pGRF). The increased HCFs in CMC were potentially caused by higher muscle forces resulting from co-contraction of agonists and antagonists around the hip. Alternatively, these higher HCFs may be caused by the slightly poorer tracking of the net joint moment by the muscle moments calculated by CMC. We conclude that the use of different optimization techniques affects calculated HCFs, and static optimization approached experimental values best.

Hip contact force in presence of aberrant bone geometry during normal and pathological gait

Children with cerebral palsy (CP) often present aberrant hip geometry, specifically increased femoral anteversion and neck-shaft angle. Furthermore, altered gait patterns are present within this population. We analyzed the effect of aberrant femoral geometry, as present in CP subjects, on hip contact force (HCF) during pathological and normal gait. We ran dynamic simulations of CP-specific and normal gait using two musculoskeletal models (MSMs), one reflecting normal femoral geometry and one reflecting proximal femoral deformities. The combination of aberrant bone geometry and CP-specific gait characteristics reduced HCF compared to normal gait on a CP subject-specific MSM, but drastically changed the orientation of the HCF vector. The HCF was orientated more vertically and anteriorly than compared to HCF orientation during normal gait. Furthermore, subjects with more pronounced bony deformities encountered larger differences in resultant HCF and HCF orientation. When bone deformities were not accounted for in MSMs of pathologic gait, the HCF orientation was more similar to normal children. Thus, our results support a relation between aberrant femoral geometry and joint loading during pathological/normal gait and confirm a compensatory effect of altered gait kinematics on joint loading.

Test-retest reliability of innovated strength tests for hip muscles

The burden of hip muscles weakness and its relation to other impairments has been well documented. It is therefore a pre-requisite to have a reliable method for clinical assessment of hip muscles function allowing the design and implementation of a proper strengthening program. Motor-driven dynamometry has been widely accepted as the gold-standard for lower limb muscle strength assessment but is mainly related to the knee joint. Studies focusing on the hip joint are less exhaustive and somewhat discrepant with regard to optimal participants position, consequently influencing outcome measures. Thus, we aimed to develop a standardized test setup for the assessment of hip muscles strength, i.e. flexors/extensors and abductors/adductors, with improved participant stability and to define its psychometric characteristics. Eighteen participants performed unilateral isokinetic and isometric contractions of the hip muscles in the sagittal and coronal plane at two separate occasions. Peak torque and normalized peak torque were measured for each contraction. Relative and absolute measures of reliability were calculated using the intraclass correlation coefficient and standard error of measurement, respectively. Results from this study revealed higher levels of between-day reliability of isokinetic/isometric hip abduction/flexion peak torque compared to existing literature. The least reliable measures were found for hip extension and adduction, which could be explained by a less efficient stabilization technique. Our study additionally provided a first set of reference normalized data which can be used in future research.

Influence of plasticization time, curing conditions, storage time, and core properties on the drug release from Aquacoat-coated pellets

Theophylline or chlorpheniramine maleate pellets were coated with an aqueous ethylcellulose dispersion, Aquacoat. The influence of the plasticization time, curing conditions, storage time, and core properties on the drug release were investigated. The plasticization time (time between plasticizer addition to the polymer dispersion and the spraying process) did not affect the drug release, when the water-soluble plasticizer triethyl citrate, was used because of its rapid uptake by the colloidal polymer particles. In contrast, with the water-insoluble plasticizer acetyltributyl citrate (ATBC), plasticization time (1/2 h vs 24 h) influenced the drug release, the longer plasticization time resulted in a slower drug release because of a more complete plasticizer uptake prior to the coating step. However a thermal aftertreatment of the coated pellets at eleylated temperatures (curing step) reduced/eliminated the effect of the plasticization time with ATBC. In general, curing reduced the drug release and resulted in stable drug release profiles. The time period between the coating and the curing step was not critical when the pellets were cured for a longer time. The structure of the pellet core (high dose matrix vs low dose layered pellet) strongly affected the drug release. A slow, zero-order drug release was obtained with high dose theophylline pellets, while a more rapid, first-order release pattern was obtained with low dose theophylline-layered nonpareil pellets.

Tackiness of acrylic and cellulosic polymer films used in the coating of solid dosage forms

The objective was to determine the tackiness of acrylic and cellulosic polymer films in order to make predictions on the tackiness (agglomeration) of coated dosage forms during coating and curing. Force-displacement curves of the detachment process of two polymeric films were used as a measure of tackiness. Various polymers (cellulosic (Aquacoat and acrylics (Eudragit RS 30D, L 30D, NE 30D)), plasticizers (triacetin, triethyl citrate, tributyl citrate, acetyltributyl citrate) and anti-tacking agents (talc and glyceryl monostearate) were investigated. The order of tackiness for films prepared from the different aqueous polymer dispersions was in order of Eudragit NE 30D > RS 30D > RL 30D > Aquacoat. The tackiness increased with increasing plasticizer concentration due to the softening of the polymer. A correlation between the minimum film formation temperature and the tackiness was observed, however, no correlation between the tackiness and the lipophilicity of the plasticizer was seen. Talc and glyceryl monostearate (GMS) reduced the tackiness of the films significantly, with GMS being effective at much lower concentrations. Curing of Eudragit RS 30D-coated theophylline beads at temperatures higher than 40 degrees C in an irreversible agglomeration of the beads and damage of the coating upon separation of the beads. This resulted in a faster release than with uncured beads. Blending the beads with talc just prior to the curing step eliminated the agglomeration and therefore film damage, even at a curing temperature of 60 degrees C.

Drug release from beads coated with an aqueous colloidal ethylcellulose dispersion, Aquacoat, or an organic ethylcellulose solution

The objective was to investigate several factors (composition of the coating formulation, the type and pH of the release medium and curing conditions), which influence the drug release from beads coated with either the aqueous ethylcellulose dispersion, Aquacoat or an organic ethylcellulose solution. The chlorpheniramine maleate release from Aquacoat-coated beads was faster in pH 7.4 buffer than in 0.1 N HCl. Increasing the curing time and curing temperature decreased the drug release in pH 7.4 buffer but did not affect the release in 0.1 N HCl. In contrast, the drug release from beads coated with the ethanolic ethylcellulose solution was not affected by the curing step, the release medium or the addition of sodium lauryl sulfate. Scanning electron microscopy and contact angle measurements explained the release data. The differences in the drug release behavior of aqueous--and organic solvent--ethylcellulose--coated beads could be attributed to the differences in the film formation process.

Multiprocessor and memory architecture of the neurocomputer SYNAPSE-1

A general purpose neurocomputer, SYNAPSE-1, which exhibits a multiprocessor and memory architecture is presented. It offers wide flexibility with respect to neural algorithms and a speed-up factor of several orders of magnitude--including learning. The computational power is provided by a 2-dimensional systolic array of neural signal processors. Since the weights are stored outside these NSPs, memory size and processing power can be adapted individually to the application needs. A neural algorithms programming language, embedded in C(+2) has been defined for the user to cope with the neurocomputer. In a benchmark test, the prototype of SYNAPSE-1 was 8000 times as fast as a standard workstation.