Viscoelastic hydrogels regulate adipose-derived mesenchymal stem cells for nucleus pulposus regeneration.

Low back pain is a leading cause of disability worldwide, often attributed to intervertebral disc (IVD) degeneration with loss of the functional nucleus pulposus (NP). Regenerative strategies utilizing biomaterials and stem cells are promising for NP repair. Human NP tissue is highly viscoelastic, relaxing stress rapidly under deformation. However, the impact of tissue-specific viscoelasticity on the activities of adipose-derived stem cells (ASC) remains largely unexplored. Here, we investigated the role of matrix viscoelasticity in regulating ASC differentiation for IVD regeneration. Viscoelastic alginate hydrogels with stress relaxation time scales ranging from 100 s to 1000s were developed and used to culture human ASCs for 21 days. Our results demonstrated that the fast-relaxing hydrogel significantly enhanced ASCs long-term cell survival and NP-like extracellular matrix secretion of aggrecan and type-II collagen. Moreover, gene expression analysis revealed a substantial upregulation of the mechanosensitive ion channel marker TRPV4 and NP-specific markers such as SOX9, HIF-1α, KRT18, CDH2 and CD24 in ASCs cultured within the fast-relaxing hydrogel, compared to slower-relaxing hydrogels. These findings highlight the critical role of matrix viscoelasticity in regulating ASC behavior and suggest that viscoelasticity is a key parameter for novel biomaterials design to improve the efficacy of stem cell therapy for IVD regeneration. STATEMENT OF SIGNIFICANCE: Systematically characterized the influence of tissue-mimetic viscoelasticity on ASC. NP-mimetic hydrogels with tunable viscoelasticity and tissue-matched stiffness. Long-term survival and metabolic activity of ASCs are substantially improved in the fast-relaxing hydrogel. The fast-relaxing hydrogel allows higher rate of cell protrusions formation and matrix remodeling. ASC differentiation towards an NP-like cell phenotype is promoted in the fast-relaxing hydrogel, with more CD24 positive expression indicating NP committed cell fate. The expression of TRPV4, a molecular sensor of matrix viscoelasticity, is significantly enhanced in the fast-relaxing hydrogel, indicating ASC sensing matrix viscoelasticity during cell development. The NP-specific ECM secretion of ASC is considerably influenced by matrix viscoelasticity, where the deposition of aggrecan and type-II collagen are significantly enhanced in the fast-relaxing hydrogel.

Machine learning applications in spine biomechanics.

Spine biomechanics is at a transformation with the advent and integration of machine learning and computer vision technologies. These novel techniques facilitate the estimation of 3D body shapes, anthropometrics, and kinematics from as simple as a single-camera image, making them more accessible and practical for a diverse range of applications. This study introduces a framework that merges these methodologies with traditional musculoskeletal modeling, enabling comprehensive analysis of spinal biomechanics during complex activities from a single camera. Additionally, we aim to evaluate their performance and limitations in spine biomechanics applications. The real-world applications explored in this study include assessment in workplace lifting, evaluation of whiplash injuries in car accidents, and biomechanical analysis in professional sports. Our results demonstrate potential and limitations of various algorithms in estimating body shape, kinematics, and conducting in-field biomechanical analyses. In industrial settings, the potential to utilize these new technologies for biomechanical risk assessments offers a pathway for preventive measures against back injuries. In sports activities, the proposed framework provides new opportunities for performance optimization, injury prevention, and rehabilitation. The application in forensic domain further underscores the wide-reaching implications of this technology. While certain limitations were identified, particularly in accuracy of predictions, complex interactions, and external load estimation, this study demonstrates their potential for advancement in spine biomechanics, heralding an optimistic future in both research and practical applications.

Fibrin clot fracture under cyclic fatigue and variable rate loading.

Fibrin clot is a vital class of fibrous materials, governing the mechanical response of blood clots. Fracture behavior of fibrin clots under complex physiological load is relevant for hemostasis and thrombosis. But how they fracture under cyclic and variable rate loading has not been reported. Here we conduct cyclic fatigue and monotonic variable rate loading tests on fibrin clots to characterize their fracture properties in terms of fatigue threshold and rate-dependent fracture toughness. We demonstrate that the fracture behavior of fibrin clots is sensitive to the amplitude of cyclic load and the loading rate. The cyclic fatigue tests show the fatigue threshold of fibrin clots at 1.66 J/m2, compared to the overall fracture toughness 15.8 J/m2. Furthermore, we rationalize the fatigue threshold using a semi-empirical model parameterized by 3D morphometric quantification to account for the hierarchical molecular structure of fibrin fibers. The variable loading tests reveal rate dependence of the overall fracture toughness of fibrin clots. Our analysis with a viscoelastic fracture model suggests the viscoelastic origin of the rate-dependent fracture toughness. The toughening mechanism of fibrin clots is further compared with biological tissues and hydrogels. This study advances the understanding and modeling of fatigue and fracture of blood clots and would motivate further investigation on the mechanics of fibrous materials. STATEMENT OF SIGNIFICANCE: Fibrin clot is a soft fibrous gel, exhibiting nonlinear mechanical responses under complex physiological loads. It is the main load-bearing constituent of blood clots where red blood cells, platelets and other cells are trapped. How the fibrin clot fractures under complex mechanical loads is critical for hemostasis and thrombosis. We study the fracture behavior of fibrin clots under cyclic fatigue and monotonic variable rate loads. We characterize the fatigue-threshold and viscous energy dissipation of fibrin clots. We compare the toughness enhancement of fibrin clots with hydrogels. The findings offer new insights into the fatigue and fracture of blood clots and fibrous materials, which could improve design guidelines for bioengineered materials.

Comparative evaluation of different spinal stability metrics.

Direct in vivo measurements of spinal stability are not possible, leaving computational estimations (such as dynamic time series and structural analyses) as the feasible option. However, differences between different stability assessment approaches and metrics remain unclear. To explore this, we asked 32 participants to perform 35 cycles of repetitive lifts with and without load (4/2.6 kg for males/females). EMG signals and 3D kinematics were collected via 12 surface electrodes and 17 inertial sensors, and three dynamical stability measures were computed: short and long temporal and conventional maximum Lyapunov exponents (LyE) and maximum Floquet multipliers (FM). A dynamic subject-specific EMG-assisted musculoskeletal model computed four structural stability measures (critical muscle stiffness coefficient at which spine becomes unstable, average spine stiffness, minimum and geometric average of Hessian matrix eigenvalues). Across cycles, dynamical and structural stability outcomes varied noticeably. Temporal short-term LyE and all structural stability measures were more influenced by the cycle percentage (posture factor) than by phase (lifting, lowering) or load factor. The effect of all factors were non-significant for FM and long LyE, except for the posture on LyE-L with a small effect size. Pearson's correlations revealed a weak to moderate, or non-existent, correlation between structural and dynamical stability metrics, with small shared variances, underscoring their distinct and independent nature and theoretical foundations. Moreover, the low sensitivity of dynamic measures to posture and load factors, found in this study, calls for further examination. Considering the limitations and shortcomings of both dynamical and structural stability assessment approaches, there is a need for the development of improved musculoskeletal stability evaluation techniques.

Tissue-mimetic hybrid bioadhesives for intervertebral disc repair.

Intervertebral disc (IVD) degeneration and herniation often necessitate surgical interventions including a discectomy with or without a nucleotomy, which results in a loss of the normal nucleus pulposus (NP) and a defect in the annulus fibrosus (AF). Due to the limited regenerative capacity of the IVD tissue, the annular tear may remain a persistent defect and result in recurrent herniation post-surgery. Bioadhesives are promising alternatives but show limited adhesion performance, low regenerative capacity, and inability to prevent re-herniation. Here, we report hybrid bioadhesives that combine an injectable glue and a tough sealant to simultaneously repair and regenerate IVD post-nucleotomy. The glue fills the NP cavity while the sealant seals the AF defect. Strong adhesion occurs with the IVD tissues and survives extreme disc loading. Furthermore, the glue can match native NP mechanically, and support the viability and matrix deposition of encapsulated cells, serving as a suitable cell delivery vehicle to promote NP regeneration. Besides, biomechanical tests with bovine IVD motion segments demonstrate the capacity of the hybrid bioadhesives to restore the biomechanics of bovine discs under cyclic loading and to prevent permanent herniation under extreme loading. This work highlights the synergy of bioadhesive and tissue-engineering approaches. Future works are expected to further improve the tissue specificity of bioadhesives and prove their efficacy for tissue repair and regeneration.

Trunk muscle forces and spinal loads during heavy deadlift: Effects of personalization, muscle wrapping, muscle lever arm, and lumbopelvic rhythm.

Heavy deadlift is used as a physical fitness screening tool in the U.S. Army. Despite the relevance of such a screening tool to military tasks performed by Service Members, the biomechanical impact of heavy deadlift and its risk of low-back injury remain unknown. A kinematics-driven musculoskeletal model of spine was implemented to investigate biomechanics of the lower back in a volunteer (23 years old, height of 1.82 m, and body mass of 98.8 kg) during a 68 kg deadlift. In search of protective mechanisms, effects of model personalization and variations in trunk musculature and lumbopelvic rhythm were also investigated. The net moment, compression and shear forces at the L5-S1 reached peaks of 684 Nm, 17.2 and 4.2 kN, respectively. Geometrical personalization and changes in lumbopelvic rhythm had the least effects on predictions while increases in muscle moment arms (40%) had the largest effects that caused, respectively, 32% and 36% decrease in the maximum compressive and shearing forces. Initiating wrapping of back muscles at farther distances from the spine had opposing effects on spinal loads; peak compression at the L5-S1 decreased by 12% whereas shear increased by 19%. Despite mechanisms considered, spinal loads during heavy deadlift exceed the existing evidence concerning the threshold of injury for spinal segments, suggesting the vulnerability to injury. Chronic exposure to such high-spinal loads may lead to (micro) fractures, degeneration, pathoanatomical changes and finally low-back pain.

Validation of an EMG submaximal method to calibrate a novel dynamic EMG-driven musculoskeletal model of the trunk: Effects on model estimates.

BACKGROUND: Multijoint EMG-assisted optimization models are reliable tools to predict muscle forces as they account for inter- and intra-individual variations in activation. However, the conventional method of normalizing EMG signals using maximum voluntary contractions (MVCs) is problematic and introduces major limitations. The sub-maximal voluntary contraction (SVC) approaches have been proposed as a remedy, but their performance against the MVC approach needs further validation particularly during dynamic tasks. METHODS: To compare model outcomes between MVC and SVC approaches, nineteen healthy subjects performed a dynamic lifting task with two loading conditions. RESULTS: Results demonstrated that these two approaches produced highly correlated results with relatively small absolute and relative differences (<10 %) when considering highly-aggregated model outcomes (e.g. compression forces, stability indices). Larger differences were, however, observed in estimated muscle forces. Although some model outcomes, e.g. force of abdominal muscles, were statistically different, their effect sizes remained mostly small (ηG2 ≤ 0.13) and in a few cases moderate (ηG2 ≤ 0.165). CONCLUSION: The findings highlight that the MVC calibration approach can reliably be replaced by the SVC approach when the true MVC exertion is not accessible due to pain, kinesiophobia and/or the lack of proper training.

Computational biomechanics of human knee joint in stair ascent: Muscle-ligament-contact forces and comparison with level walking.

About a third of knee joint disorders originate from the patellofemoral (PF) site that makes stair ascent a difficult activity for patients. A detailed finite element model of the knee joint is coupled to a lower extremity musculoskeletal model to simulate the stance phase of stair ascent. It is driven by the mean of measurements on the hip-knee-ankle moments-angles as well as ground reaction forces reported in healthy individuals. Predicted muscle activities compare well to the recorded electromyography data. Peak forces in quadriceps (3.87 BW, body weight, at 20% instance in our 607 N subject), medial hamstrings (0.77 BW at 20%), and gastrocnemii (1.21 BW at 80%) are estimated. Due to much greater flexion angles-moments in the first half of stance, large PF contact forces (peak of 3.1 BW at 20% stance) and stresses (peak of 4.83 MPa at 20% stance) are estimated that exceed their peaks in level walking by fourfold and twofold, respectively. Compared with level walking, ACL forces diminish in the first half of stance but substantially increase later in the second half (peak of 0.76 BW at 75% stance). Under nearly similar contact forces at 20% of stance, the contact stress on the tibiofemoral (TF) medial plateau reaches a peak (9.68 MPa) twice that on the PF joint suggesting the vulnerability of both joints. Compared with walking, stair ascent increases peak ACL force and both peak TF and PF contact stresses. Reductions in the knee flexion moment and/or angle appear as a viable strategy to mitigate internal loads and pain.

Submaximal electromyography-driven musculoskeletal modeling of the human trunk during static tasks: Equilibrium and stability analyses.

Conventional electromyography-driven (EMG) musculoskeletal models are calibrated during maximum voluntary contraction (MVC) tasks, but individuals with low back pain cannot perform unbiased MVCs. To address this issue, EMG-driven models can be calibrated in submaximal tasks. However, the effects of maximal (when data points include the maximum contraction) and submaximal calibration techniques on model outputs (e.g., muscle forces, spinal loads) remain yet unknown. We calibrated a subject-specific EMG-driven model, using maximal/submaximal isometric contractions, and simulated different independent tasks. Both approaches satisfactorily predicted external moments (Pearson's correlation ∼ 0.75; relative error = 44%), and removing calibration tasks under axial torques markedly improved the model performance (Pearson's correlation âˆ¼ 0.92; relative error âˆ¼ 28%). Unlike individual muscle forces, gross (aggregate) model outputs (i.e., spinal loads, stability index, and sum of abdominal/back muscle forces) estimated from maximal and submaximal calibration techniques were highly correlated (r > 0.78). Submaximal calibration method overestimated spinal loads (6% in average) and abdominal muscle forces (11% in average). Individual muscle forces estimated from maximal and submaximal approaches were substantially different; however, gross model outputs (especially internal loads and stability index) remained highly correlated with small to moderate relative differences; therefore, the submaximal calibration technique can be considered as an alternative to the conventional maximal calibration approach.

Blood clot behaves as a poro-visco-elastic material.

Mechanical properties of blood clots play a crucial role in hemostasis and embolization. They are time-dependent and often described with viscoelastic models. But blood clots also exhibit some hallmark features of poroelasticity, as most biological tissues exhibit concurrent viscoelasticity and poroelasticity. In this study, we characterized the time-dependent behavior of blood clots, and developed an experimental-computational framework to decouple and model poroelastic and viscoelastic responses. Compression stress relaxation tests were conducted on bovine blood clots with different diameters to evaluate the influence of sample size on the relaxation time. In the compression tests, the mass of blood clots was measured to estimate fluid migration. To capture pure viscoelastic responses, rheological shear stress relaxation tests were carried out. A poroviscoelastic model was also proposed and calibrated to capture the complex multiaxial (compression and shear) relaxation behavior of blood clots. In unconfined compression tests, stresses relaxed markedly (average: 83%; range: 76-90%), and samples with larger diameters showed longer relaxation time. Blood clots lost about 24% of their initial masses, and the mass transport took place gradually in compression tests. Under shear deformation, blood clots relaxed in average 37% (range: 32-39%) which was much less than those under unconfined compression tests (in average 37% versus 83%). Unlike poroelastic and viscoelastic theories, the poroviscoelastic model accurately predicted multiaxial responses of blood clots under compression and shear; additionally, the estimated Darcy's coefficient (4.4×10-9 cm2) was found within the reported physiological range (0.1×10-9 to 36 × 10-9 cm2). The combination of size-dependent stress relaxation and mass loss under compression (due to poroelasticity) as well as substantial stress relaxation under shear deformation (due to viscoelasticity) demonstrate that blood clot behaves as a poroviscoelastic material; therefore, accurate interpretation of transient responses of blood clots requires a validated poroviscoelastic model. This work provides understanding and methodology on blood clot mechanics and will further motivate the development of clot-like biomaterials.

Evaluating stability of human spine in static tasks: a combined in vivo-computational study.

Various interpretations and parameters have been proposed to assess spinal stability such as antagonist muscle coactivity, trunk stiffness and spinal buckling load; however, the correlation between these parameters remains unknown. We evaluated spinal stability during different tasks while changing the external moment and load height and investigated likely relationships between different EMG- and model-based parameters (e.g., EMG coactivity ratio, trunk stiffness, force coactivity ratio) and stability margins. EMG and kinematics of 40 young healthy subjects were recorded during various quasi-static tasks. Muscle forces, trunk stiffness and stability margins were calculated by a nonlinear subject-specific EMG-assisted-optimization musculoskeletal model of the trunk. The load elevation and external moment increased muscle activities and trunk stiffness while all stability margins (i.e., buckling loads) decreased. The force coactivity ratio was strongly correlated with the hand-load stability margin (i.e., additional weight in hands to initiate instability; R2 = 0.54) demonstrating the stabilizing role of abdominal muscles. The total trunk stiffness (Pearson's r = 0.96) and the sum of EMGs of back muscles (Pearson's r = 0.65) contributed the most to the T1 stability margin (i.e., additional required load at T1 for instability/buckling). Force coactivity ratio and trunk stiffness can be used as alternative spinal stability metrics.

Modeling of human intervertebral disc annulus fibrosus with complex multi-fiber networks.

Collagen fibers within the annulus fibrosus (AF) lamellae are unidirectionally aligned with alternating orientations between adjacent layers. AF constitutive models often combine two adjacent lamellae into a single equivalent layer containing two fiber networks with a crisscross pattern. Additionally, AF models overlook the inter-lamellar matrix (ILM) as well as elastic fiber networks in between lamellae. We developed a nonhomogenous micromechanical model as well as two coarser homogenous hyperelastic and microplane models of the human AF, and compared their performances against measurements (tissue level uniaxial and biaxial tests as well as whole disc experiments) and seven published hyperelastic models. The micromechanical model had a realistic non-homogenous distribution of collagen fiber networks within each lamella and elastic fiber network in the ILM. For small matrix linear moduli (<0.2 MPa), the ILM showed substantial anisotropy (>10%) due to the elastic fiber network. However, at moduli >0.2 MPa, the effects of the elastic fiber network on differences in stress-strain responses at different directions disappeared (<10%). Variations in sample geometry and boundary conditions (due to uncertainty) markedly affected stress-strain responses of the tissue in uniaxial and biaxial tests (up to 16 times). In tissue level tests, therefore, simulations should represent testing conditions (e.g., boundary conditions, specimen geometry, preloads) as closely as possible. Stress/strain fields estimated from the single equivalent layer approach (conventional method) yielded different results from those predicted by the anatomically more accurate apparoach (i.e., layerwise). In addition, in a disc under a compressive force (symmetric loading), asymmetric stress-strain distributions were computed when using a layerwise simulation. Although all developed and selected published AF models predicted gross compression-displacement responses of the whole disc within the range of measured data, some showed excessively stiff or compliant responses under tissue-level uniaxial/biaxial tests. This study emphasizes, when constructing and validating constitutive models of AF, the importance of the proper simulation of individual lamellae as distinct layers, and testing parameters (sample geometric dimensions/loading/boundary conditions).

Comparison of different lifting analysis tools in estimating lower spinal loads - Evaluation of NIOSH criterion.

Excessive loads on the human spine is recognized as a risk factor for back injuries/pain. Various lifting analysis tools such as musculoskeletal models, regression equations and NIOSH (National Institute for Occupational Safety and Health) lifting equation (NLE) have been proposed to evaluate and mitigate associated risks during manual material handling activities. Present study aims to compare predicted spinal loads from 5 different lifting analysis tools as well as to critically evaluate the NIOSH recommended weight limit (RWL). Spinal loads were estimated under different symmetric/asymmetric lifting tasks in which hand-load mass at each task was set based on RWL from NLE. Estimated intradiscal pressures (IDPs) of various tools were also compared with in vivo measurements. We compared RWL by NLE versus our estimations of RWL calculated from our regression equations using biomechanical criteria (compression <3400 N with/without shear <1000, 1250 or 1500 N). Our regression equations followed by OpenSim, AnyBody, simple polynomial and 3DSSPP satisfactorily predicted L4-L5 IDPs. Lifting analysis tools estimated comparable spinal compression forces (mean Pearson's r = 0.80; standard deviation of relative difference = 26%) while in shear, differences were greater (mean Pearson's r = 0.68; standard deviation of relative difference = 56%). NLE estimations of RWL were conservative in comparison with our estimations for lean individuals (BMI < 25 kg/m2) when compression <3400 N and shear <1250 N were considered as the biomechanical criteria. For heavier individuals, however, NLE estimations of RWL generated spinal compression >3400 N (NIOSH biomechanical safety threshold) as well as shear >1000 N. Although RWLs estimated by NLE was body weight independent, body weight substantially altered RWLs estimated from our regression equations. For improved estimation of the risk of injury, more accurate failure criteria for spinal segments are essential.

Internal load-sharing in the human passive lumbar spine: Review of in vitro and finite element model studies.

Human lumbar motion segment is composed of various components with distinct contributions to its gross mechanical response. By employing experimental and computational approaches, many studies have investigated the relative role of each component as well as effects of various factors such as boundary-initial conditions, load magnitude-combination-direction, load temporal regime, preload, posture, degeneration, failures and surgical interventions on load-sharing. This paper reviews and critically discusses the relevant findings of in vitro and finite element model studies on load-sharing in healthy, aged, degenerate and damaged human lumbar motion segments. Two systematic searches were performed in PubMed (October 2018 - March 2019) using three sets of concepts ("lumbar spine", "load-sharing" and "motion segment components") followed by a complementary generic search. The segment overall response as well as the relative role of its constituents are markedly influenced by alterations in resection sequence, boundary conditions, geometry, loading characteristics (rate, magnitude, combinations and preloads), disc hydration, bone quality, posture and time (creep and cyclic). Structural transection order affects both findings and conclusions not only in force-control protocols but also in displacement-control loading regimes. Disc degeneration, endplate fracture and surgical resections significantly alter load transmission in the lumbar spine. In summary, in vitro and finite element model studies have together substantially improved our understanding of functional biomechanics (load-sharing) of human lumbar spine in normal and perturbed conditions acting as invaluable complementary tools in clinical applications.

On the modeling of human intervertebral disc annulus fibrosus: Elastic, permanent deformation and failure responses.

As a primary load-resisting component, annulus fibrosus (AF) maintains structural integrity of the entire intervertebral disc. Experiments have demonstrated that permanent deformation and damage take place in the tissue under mechanical loads. Development of an accurate model to capture the complex behaviour of AF tissue is hence crucial in disc model studies. We, therefore, aimed to develop a non-homogenous model to capture elastic, inelastic and failure responses of the AF tissue and the entire disc model under axial load. Our model estimations satisfactorily agreed with results of existing uniaxial (along fiber, circumferential and axial directions) and biaxial tissue-level tests. The model accurately predicted the failure of the tissue in various directions in uniaxial extension. Collagen fiber content, type and orientation substantially altered AF tissue responses in uni- and bi-axial tests. Although collagen fiber content and type mostly affected failure stress, fiber orientation significantly influenced the tissue failure strain. The entire L2-L3 disc model accurately replicated load-displacement as well as loading-unloading responses of the disc under compression-tension forces. Preconditioning of the disc-body unit substantially stiffened response. Poisson's ratio of both AF and nucleus considerably affected compression-displacement responses of the disc (173% increase in compression at 1.49 mm displacement when it was changed from 0.499 to 0.49999). Any AF constitutive model should be calibrated under various tissue-level loads and directions as well as the entire disc model responses since using a single tissue-level loading (e.g. uniaxial) for calibration can lead to unrealistic responses in other tests (e.g., biaxial). Special attentions should be given to the choice of Poisson's ratio and the realistic consideration of preconditioning load.

Comparison of spinal stability following motor control and general exercises in nonspecific chronic low back pain patients.

BACKGROUND: Motor control exercise was claimed to improve spinal stability in patients with chronic non-specific back pain, but to investigate the effectiveness of this exercise, other outcome measures have been used rather than spinal stability itself. The aim of our study is to assess motor control exercise effects on spinal stability using a biomechanical model. METHODS: Fifty-one patients were assigned to either motor control or general exercises. Before and after trainings, participants were tested for spinal stability at seven isometric tasks. Electromyography signals were recorded from ten superficial muscles, and a hybrid EMG-driven musculoskeletal model estimated spinal stability indices at each task. FINDINGS: Pain and disability significantly decreased in both groups. After trainings, patients had both increase and decrease in stability depending on the task, and stability did not increase/decrease uniformly in all patients. In the motor control group, stability increased at all positions but reached to significance only at right lateral pulling. In the general exercise group, except for pulling the trunk backward, stability decreased at other positions and reached to statistical significance only at pulling the trunk forward. No significant difference between groups was found in changing stability after the intervention. INTERPRETATION: Interventions yielded no significant difference in disability, pain and stability index between two groups. Significant increase of stability in the motor control group at right lateral pulling may be attributed to more activity of abdominal muscles, and significant decrease of stability in the general exercise group at forward pulling may be attributed to more optimal activity of back muscles.

Obesity and Obesity Shape Markedly Influence Spine Biomechanics: A Subject-Specific Risk Assessment Model.

Underlying mechanisms of obesity-related back pain remain unexplored. Thus, we aim to determine the effect of obesity and its shapes on the spinal loads and the associated risks of injury. Obesity shapes were initially constructed by principal component analysis based on datasets on 5852 obese individuals. Spinal loads, cycles to vertebral failure and trunk stability margin were estimated in a subject-specific trunk model taking account of personalized musculature, passive ligamentous spine, obesity shapes, segmental weights, spine kinematics and bone mineral density. Three obesity shapes (mean and extreme abdominal circumferences) at three body weights (BWs) of 86, 98 and 109 kg were analyzed. Additional BW (12 kg) increased spinal loads by ~11.8%. Higher waist circumferences at identical BW increased spinal forces to the tune of ~20 kg additional BW and the risk of vertebral fatigue compression fracture by 3-7 times when compared with smaller waist circumferences. Forward flexion, greater BW and load in hands increased the trunk stability margin. Spinal loads markedly increased with BW, especially at greater waist circumferences. The risk of vertebral fatigue fracture also substantially increased at greater waist circumferences though not at smaller ones. Obesity and its shape should be considered in spine biomechanics.

Effect of intervertebral translational flexibilities on estimations of trunk muscle forces, kinematics, loads, and stability.

Due to the complexity of the human spinal motion segments, the intervertebral joints are often simulated in the musculoskeletal trunk models as pivots thus allowing no translational degrees of freedom (DOFs). This work aims to investigate, for the first time, the effect of such widely used assumption on trunk muscle forces, spinal loads, kinematics, and stability during a number of static activities. To address this, the shear deformable beam elements used in our nonlinear finite element (OFE) musculoskeletal model of the trunk were either substantially stiffened in translational directions (SFE model) or replaced by hinge joints interconnected through rotational springs (HFE model). Results indicated that ignoring intervertebral translational DOFs had in general low to moderate impact on model predictions. Compared with the OFE model, the SFE and HFE models predicted generally larger L4-L5 and L5-S1 compression and shear loads, especially for tasks with greater trunk angles; differences reached ~15% for the L4-L5 compression, ~36% for the L4-L5 shear and ~18% for the L5-S1 shear loads. Such differences increased, as location of the hinge joints in the HFE model moved from the mid-disc height to either the lower or upper endplates. Stability analyses of these models for some select activities revealed small changes in predicted margin of stability. Model studies dealing exclusively with the estimation of spinal loads and/or stability may, hence with small loss of accuracy, neglect intervertebral translational DOFs at smaller trunk flexion angles for the sake of computational simplicity.