On the relative relevance of subject-specific geometries and degeneration-specific mechanical properties for the study of cell death in human intervertebral disk models

Unveiling interactions between intervertebral disc morphologies and mechanical behavior through personalized finite element modeling

Introduction: Intervertebral Disc (IVD) Degeneration (IDD) is a significant health concern, potentially influenced by mechanotransduction. However, the relationship between the IVD phenotypes and mechanical behavior has not been thoroughly explored in local morphologies where IDD originates. This work unveils the interplays among morphological and mechanical features potentially relevant to IDD through Abaqus UMAT simulations. Methods: A groundbreaking automated method is introduced to transform a calibrated, structured IVD finite element (FE) model into 169 patient-personalized (PP) models through a mesh morphing process. Our approach accurately replicates the real shapes of the patient's Annulus Fibrosus (AF) and Nucleus Pulposus (NP) while maintaining the same topology for all models. Using segmented magnetic resonance images from the former project MySpine, 169 models with structured hexahedral meshes were created employing the Bayesian Coherent Point Drift++ technique, generating a unique cohort of PP FE models under the Disc4All initiative. Machine learning methods, including Linear Regression, Support Vector Regression, and eXtreme Gradient Boosting Regression, were used to explore correlations between IVD morphology and mechanics. Results: We achieved PP models with AF and NP similarity scores of 92.06\% and 92.10\% compared to the segmented images. The models maintained good quality and integrity of the mesh. The cartilage endplate (CEP) shape was represented at the IVD-vertebra interfaces, ensuring personalized meshes. Validation of the constitutive model against literature data showed a minor relative error of 5.20%. Discussion: Analysis revealed the influential impact of local morphologies on indirect mechanotransduction responses, highlighting the roles of heights, sagittal areas, and volumes. While the maximum principal stress was influenced by morphologies such as heights, the disc's ellipticity influenced the minimum principal stress. Results suggest the CEPs are not influenced by their local morphologies but by those of the AF and NP. The generated free-access repository of individual disc characteristics is anticipated to be a valuable resource for the scientific community with a broad application spectrum.

Comparative FEM study on intervertebral disc modeling: Holzapfel-Gasser-Ogden vs. structural rebars

Introduction: Numerical modeling of the intervertebral disc (IVD) is challenging due to its complex and heterogeneous structure, requiring careful selection of constitutive models and material properties. A critical aspect of such modeling is the representation of annulus fibers, which significantly impact IVD biomechanics. This study presents a comparative analysis of different methods for fiber reinforcement in the annulus fibrosus of a finite element (FE) model of the human IVD. Methods: We utilized a reconstructed L4-L5 IVD geometry to compare three fiber modeling approaches: the anisotropic Holzapfel-Gasser-Ogden (HGO) model (HGO fiber model) and two sets of structural rebar elements with linear-elastic (linear rebar model) and hyperelastic (nonlinear rebar model) material definitions, respectively. Prior to calibration, we conducted a sensitivity analysis to identify the most important model parameters to be calibrated and improve the efficiency of the calibration. Calibration was performed using a genetic algorithm and in vitro range of motion (RoM) data from a published study with eight specimens tested under four loading scenarios. For validation, intradiscal pressure (IDP) measurements from the same study were used, along with additional RoM data from a separate publication involving five specimens subjected to four different loading conditions. Results: The sensitivity analysis revealed that most parameters, except for the Poisson ratio of the annulus fibers and C01 from the nucleus, significantly affected the RoM and IDP outcomes. Upon calibration, the HGO fiber model demonstrated the highest accuracy (R2 = 0.95), followed by the linear (R2 = 0.89) and nonlinear rebar models (R2 = 0.87). During the validation phase, the HGO fiber model maintained its high accuracy (RoM R2 = 0.85; IDP R2 = 0.87), while the linear and nonlinear rebar models had lower validation scores (RoM R2 = 0.71 and 0.69; IDP R2 = 0.86 and 0.8, respectively). Discussion: The results of the study demonstrate a successful calibration process that established good agreement with experimental data. Based on our findings, the HGO fiber model appears to be a more suitable option for accurate IVD FE modeling considering its higher fidelity in simulation results and computational efficiency.

Dataset of Finite Element Models of Normal and Deformed Thoracolumbar Spine

Adult spine deformity (ASD) is prevalent and leads to a sagittal misalignment in the vertebral column. Computational methods, including Finite Element (FE) Models, have emerged as valuable tools for investigating the causes and treatment of ASD through biomechanical simulations. However, the process of generating personalised FE models is often complex and time-consuming. To address this challenge, we present a dataset of FE models with diverse spine morphologies that statistically represent real geometries from a cohort of patients. These models are generated using EOS images, which are utilized to reconstruct 3D surface spine models. Subsequently, a Statistical Shape Model (SSM) is constructed, enabling the adaptation of a FE hexahedral mesh template for both the bone and soft tissues of the spine through mesh morphing. The SSM deformation fields facilitate the personalization of the mean hexahedral FE model based on sagittal balance measurements. Ultimately, this new hexahedral SSM tool offers a means to generate a virtual cohort of 16807 thoracolumbar FE spine models, which are openly shared in a public repository.

Effects of geometric individualisation of a human spine model on load sharing: neuro-musculoskeletal simulation reveals significant differences in ligament and muscle contribution

In spine research, two possibilities to generate models exist: generic (population-based) models representing the average human and subject-specific representations of individuals. Despite the increasing interest in subject specificity, individualisation of spine models remains challenging. Neuro-musculoskeletal (NMS) models enable the analysis and prediction of dynamic motions by incorporating active muscles attaching to bones that are connected using articulating joints under the assumption of rigid body dynamics. In this study, we used forward-dynamic simulations to compare a generic NMS multibody model of the thoracolumbar spine including fully articulated vertebrae, detailed musculature, passive ligaments and linear intervertebral disc (IVD) models with an individualised model to assess the contribution of individual biological structures. Individualisation was achieved by integrating skeletal geometry from computed tomography and custom-selected muscle and ligament paths. Both models underwent a gravitational settling process and a forward flexion-to-extension movement. The model-specific load distribution in an equilibrated upright position and local stiffness in the L4/5 functional spinal unit (FSU) is compared. Load sharing between occurring internal forces generated by individual biological structures and their contribution to the FSU stiffness was computed. The main finding of our simulations is an apparent shift in load sharing with individualisation from an equally distributed element contribution of IVD, ligaments and muscles in the generic spine model to a predominant muscle contribution in the individualised model depending on the analysed spine level.

Patient-Specific Variations in Local Strain Patterns on the Surface of a Trussed Titanium Interbody Cage

Introduction: 3D printed trussed titanium interbody cages may deliver bone stimulating mechanobiological strains to cells attached at their surface. The exact size and distribution of these strains may depend on patient-specific factors, but the influence of these factors remains unknown. Therefore, this study aimed to determine patient-specific variations in local strain patterns on the surface of a trussed titanium interbody fusion cage.

Materials and Methods: Four patients eligible for spinal fusion surgery with the same cage size were selected from a larger database. For these cases, patient-specific finite element models of the lumbar spine including the same trussed titanium cage were made. Functional dynamics of the non-operated lumbar spinal segments, as well as local cage strains and caudal endplate stresses at the operated segment, were evaluated under physiological extension/flexion movement of the lumbar spine.

Results: All patient-specific models revealed physiologically realistic functional dynamics of the operated spine. In all patients, approximately 30% of the total cage surface experienced strain values relevant for preserving bone homeostasis and stimulating bone formation. Mean caudal endplate contact pressures varied up to 10 MPa. Both surface strains and endplate contact pressures varied more between loading conditions than between patients.

Conclusions: This study demonstrates the applicability of patient-specific finite element models to quantify the impact of patient-specific factors such as bone density, degenerative state of the spine, and spinal curvature on interbody cage loading. In the future, the same framework might be further developed in order to establish a pipeline for interbody cage design optimizations.

Patient-specific apparent diffusion maps used to model nutrient availability in degenerated intervertebral discs

INTRODUCTION

In this study, magnetic resonance imaging data was used to (1) model IVD-specific gradients of glucose, oxygen, lactate, and pH; and (2) investigate possible effects of covariate factors (i.e., disc geometry, and mean apparent diffusion coefficient values) on the IVD's microenvironment. Mathematical modeling of the patient's specific IVD microenvironment could be important when selecting patients for stem cell therapy due to the increased nutrient demand created by that treatment.

MATERIALS AND METHODS

Disc geometry and water diffusion coefficients were extracted from MRIs of 37 patients using sagittal T1-weighted images, T2-weighted images, and ADC Maps. A 2-D steady state finite element mathematical model was developed in COMSOL Multiphysics® 5.4 to compute concentration maps of glucose, oxygen, lactate and pH.

RESULTS

Concentration of nutrients (i.e., glucose, and oxygen) dropped with increasing distance from the cartilaginous endplates (CEP), whereas acidity levels increased. Most discs experienced poor nutrient levels along with high acidity values in the inner annulus fibrosus (AF). The disc's physiological microenvironment became more deficient as degeneration progressed. For example, minimum glucose concentration in grade 4 dropped by 31.1% compared to grade 3 (p < 0.0001). The model further suggested a strong effect of the following parameters: disc size, AF and CEP diffusivities, metabolic reactions, and cell density on solute concentrations in the disc (p < 0.05).

CONCLUSION

The significance of this work implies that the individual morphology and physiological conditions of each disc, even among discs of the same Pfirrmann grade, should be evaluated when modeling IVD solute concentrations.

Misaligned spinal rods can induce high internal forces consistent with those observed to cause screw pullout and disc degeneration

BACKGROUND CONTEXT

Manual contouring of spinal rods is often required intraoperatively for proper alignment of the rods within the pedicle screw heads. Residual misalignments are frequently reduced by using dedicated reduction devices. The forces exerted by these devices, however, are uncontrolled and may lead to excessive reaction forces. As a consequence, screw pullout might be provoked and surrounding tissue may experience unfavorable biomechanical loads. The corresponding loads and induced tissue deformations are however not well identified. Additionally, whether the forced reduction alters the biomechanical behavior of the lumbar spine during physiological movements postoperatively, remains unexplored.

PURPOSE

To predict whether the reduction of misaligned posterior instrumentation might result in clinical complications directly after reduction and during a subsequent physiological flexion movement.

STUDY DESIGN

Finite element analysis.

METHODS

A patient-specific, total lumbar (L1-S1) spine finite element model was available from previous research. The model consists of poro-elastic intervertebral discs with Pfirrmann grade-dependent material parameters, with linear elastic bone tissue with stiffness values related to the local bone density, and with the seven major ligaments per spinal motion segment described as nonlinear materials. Titanium instrumentation was implemented in this model to simulate a L4, L5, and S1 posterolateral fusion. Next, coronal and sagittal misalignments of 6 mm each were introduced between the rod and the screw head at L4. These misalignments were computationally reduced and a physiological flexion movement of 15° was prescribed. Non-instrumented and well-aligned instrumented models were used as control groups.

RESULTS

Pulling forces up to 1.0 kN were required to correct the induced misalignments of 6 mm. These forces affected the posture of the total lumbar spine, as motion segments were predicted to rotate up to 3 degrees and rotations propagated proximally to and even affect the L1-2 level. The facet contact pressures in the corrected misaligned models were asymmetrical suggesting non-physiological joint loading in the misaligned models. In addition, the discs and vertebrae experienced abnormally high forces as a result of the correction procedure. These effects were more pronounced after a 15° flexion movement following forced reduction.

CONCLUSIONS

The results of this study indicate that the correction of misaligned posterior instrumentation can result in high forces at the screws consistent with those reported to cause screw pullout, and may cause high-tissue strains in adjacent and downstream spinal segments.

CLINICAL SIGNIFICANCE

Proper alignment of spinal posterior instrumentation may reduce clinical complications secondary to unfavorable biomechanics.

Computational Challenges in Tissue Engineering for the Spine

This paper deals with a brief review of the recent developments in computational modelling applied to innovative treatments of spine diseases. Additionally, it provides a perspective on the research directions expected for the forthcoming years. The spine is composed of distinct and complex tissues that require specific modelling approaches. With the advent of additive manufacturing and increasing computational power, patient-specific treatments have moved from being a research trend to a reality in clinical practice, but there are many issues to be addressed before such approaches become universal. Here, it is identified that the major setback resides in validation of these computational techniques prior to approval by regulatory agencies. Nevertheless, there are very promising indicators in terms of optimised scaffold modelling for both disc arthroplasty and vertebroplasty, powered by a decisive contribution from imaging methods.

Multiscale Regulation of the Intervertebral Disc: Achievements in Experimental, In Silico, and Regenerative Research

Intervertebral disc (IVD) degeneration is a major risk factor of low back pain. It is defined by a progressive loss of the IVD structure and functionality, leading to severe impairments with restricted treatment options due to the highly demanding mechanical exposure of the IVD. Degenerative changes in the IVD usually increase with age but at an accelerated rate in some individuals. To understand the initiation and progression of this disease, it is crucial to identify key top-down and bottom-up regulations' processes, across the cell, tissue, and organ levels, in health and disease. Owing to unremitting investigation of experimental research, the comprehension of detailed cell signaling pathways and their effect on matrix turnover significantly rose. Likewise, in silico research substantially contributed to a holistic understanding of spatiotemporal effects and complex, multifactorial interactions within the IVD. Together with important achievements in the research of biomaterials, manifold promising approaches for regenerative treatment options were presented over the last years. This review provides an integrative analysis of the current knowledge about (1) the multiscale function and regulation of the IVD in health and disease, (2) the possible regenerative strategies, and (3) the in silico models that shall eventually support the development of advanced therapies.

Automated Grading of Lumbar Disc Degeneration Using a Push-Pull Regularization Network Based on MRI

BACKGROUND

Lower back pain is one of the most widely experienced health problems and is strongly associated with intervertebral disc (IVD) degeneration. The quantification of the degeneration is valuable for estimation of the material properties of the IVDs and then to perform biomechanical simulation of the spinal conditions and treatments. The MR image characteristics of relatively high intraclass variability and small interclass differences pose challenges for the classification algorithm to perform automatic grading of degenerated IVD.

PURPOSE

To assess the feasibility and improvement of a computer-assisted IVD degeneration grading method based on proposed push-pull regularization (PPR) strategy.

STUDY TYPE

Retrospective.

POPULATION

In total, 500 subjects (350 for training, 70 for validation, and 80 for test in a 10-time 10-fold cross validation setting) with varied lumbar disorders were included.

FIELD STRENGTH/SEQUENCE

3.0T; T₂ -weighted spin echo sequence.

ASSESSMENT

IVD degeneration grading was taken as a classification task of five classes according to the Pfirrmann grading system in this study. The classification results of deep-learning models with and without PPR were compared with the classifications made by three experienced spinal radiologists.

STATISTICAL TESTS

Paired t-tests.

RESULTS

The classification results show that in four classical CNN models of VGG-M, VGG-16, GoogleNet, and ResNet-34, by embedding a PPR strategy, the accuracies of grade II and III IVD classification were improved by more than 10% (P < 0.05), and that the overall accuracy (grades I to V) was improved by over 8% (P < 0.05).

DATA CONCLUSION

The embedded PPR significantly improved the classification performance, which enhanced CNN representation capability for IVD degeneration grading.

LEVEL OF EVIDENCE

3 TECHNICAL EFFICACY STAGE: 1.

Lumbar model generator: a tool for the automated generation of a parametric scalable model of the lumbar spine

Low back pain is a major cause of disability and requires the development of new devices to treat pathologies and improve prognosis following surgery. Understanding the effects of new devices on the biomechanics of the spine is crucial in the development of new effective and functional devices. The aim of this study was to develop a preliminary parametric, scalable and anatomically accurate finite-element model of the lumbar spine allowing for the evaluation of the performance of spinal devices. The principal anatomical surfaces of the lumbar spine were first identified, and then accurately fitted from a previous model supplied by S14 Implants (Bordeaux, France). Finally, the reconstructed model was defined according to 17 parameters which are used to scale the model according to patient dimensions. The developed model, available as a toolbox named the lumbar model generator, enables generating a population of models using subject-specific dimensions obtained from data scans or averaged dimensions evaluated from the correlation analysis. This toolbox allows patient-specific assessment, taking into account individual morphological variation. The models have applications in the design process of new devices, evaluating the biomechanics of the spine and helping clinicians when deciding on treatment strategies.

Theoretical Explorations Generate New Hypotheses About the Role of the Cartilage Endplate in Early Intervertebral Disk Degeneration

Altered cell nutrition in the intervertebral disk (IVD) is considered a main cause for disk degeneration (DD). The cartilage endplate (CEP) provides a major path for the diffusion of nutrients from the peripheral vasculature to the IVD nucleus pulposus (NP). In DD, sclerosis of the adjacent bony endplate is suggested to be responsible for decreased diffusion and disk cell nutrition. Yet, experimental evidence does not support this hypothesis. Hence, we evaluated how moderate CEP composition changes related to tissue degeneration can affect disk nutrition and cell viability. A novel composition-based permeability formulation was developed for the CEP, calibrated, validated, and used in a mechano-transport finite element IVD model. Fixed solute concentrations were applied at the outer surface of the annulus and the CEP, and three cycles of daily mechanical load were simulated. The CEP model indicated that CEP permeability increases with the degeneration/aging of the tissue, in accordance with recent measurements reported in the literature. Additionally, our results showed that CEP degeneration might be responsible for mechanical load-induced NP dehydration, which locally affects oxygen and lactate levels, and reduced glucose concentration by 16% in the NP-annulus transition zone. Remarkably, CEP degeneration was a condition sine-qua-non to provoke cell starvation and death, while simulating the effect of extracellular matrix depletion in DD. This theoretical study cast doubts about the paradigm that CEP calcification is needed to provoke cell starvation, and suggests an alternative path for DD whereby the early degradation of the CEP plays a key role.

Comparison of patient-specific computational models vs. clinical follow-up, for adjacent segment disc degeneration and bone remodelling after spinal fusion

Spinal fusion is a standard surgical treatment for patients suffering from low back pain attributed to disc degeneration. However, results are somewhat variable and unpredictable. With fusion the kinematic behaviour of the spine is altered. Fusion and/or stabilizing implants carrying considerable load and prevent rotation of the fused segments. Associated with these changes, a risk for accelerated disc degeneration at the adjacent levels to fusion has been demonstrated. However, there is yet no method to predict the effect of fusion surgery on the adjacent tissue levels, i.e. bone and disc. The aim of this study was to develop a coupled and patient-specific mechanoregulated model to predict disc generation and changes in bone density after spinal fusion and to validate the results relative to patient follow-up data. To do so, a multiscale disc mechanoregulation adaptation framework was developed and coupled with a previously developed bone remodelling algorithm. This made it possible to determine extra cellular matrix changes in the intervertebral disc and bone density changes simultaneously based on changes in loading due to fusion surgery. It was shown that for 10 cases the predicted change in bone density and degeneration grade conforms reasonable well to clinical follow-up data. This approach helps us to understand the effect of surgical intervention on the adjacent tissue remodelling. Thereby, providing the first insight for a spine surgeon as to which patient could potentially be treated successfully by spinal fusion and in which patient has a high risk for adjacent tissue changes.

A study of the damage behaviour of porcine intervertebral discs in a bioreactor environment

The intervertebral discs are cartilaginous, articulating structures that lie between vertebral bodies, allowing flexibility, transmission, modification, and also distribution of the forces to the spinal column. Disc degeneration is characterised by progressive loss of disc height and exaggerated radial bulging. Therefore, the spine becomes shorter, stiffer, and less mobile. In the last several decades, there is a strong need for a tissue engineering strategy that alleviates pain and restores spine function by directly addressing the underlying biological causes of disc degeneration. Numerous studies that are currently showing potential have been conducted on developing regenerative and reparative strategies for treating this condition. In this study, to numerically describe the anisotropic mechanical damage behaviour of discs, the pseudo-elastic damage model was applied. To experimentally picture the biomechanical response of discs and to study the damage mechanisms as well as the spinal disc herniation, a special bioreactor was evolved. The specimens were obtained from pigs aged six months. A total of eight functional spine units were taken from porcine lumbar spines (L1-L2). Firstly, the experiments were performed by using long-term cyclic uniaxial compression tests. Secondly, the mean value of experimental results with consideration of the different shapes and sizes of the samples was calculated. Afterwards, the experimental results were compared with outcomes of numerical simulations.

Annulus fibrosus functional extrafibrillar and fibrous mechanical behaviour: experimental and computational characterisation

The development of current surgical treatments for intervertebral disc damage could benefit from virtual environment accounting for population variations. For such models to be reliable, a relevant description of the mechanical properties of the different tissues and their role in the functional mechanics of the disc is of major importance. The aims of this work were first to assess the physiological hoop strain in the annulus fibrosus in fresh conditions (n = 5) in order to extract a functional behaviour of the extrafibrillar matrix; then to reverse-engineer the annulus fibrosus fibrillar behaviour (n = 6). This was achieved by performing both direct and global controlled calibration of material parameters, accounting for the whole process of experimental design and in silico model methodology. Direct-controlled models are specimen-specific models representing controlled experimental conditions that can be replicated and directly comparing measurements. Validation was performed on another six specimens and a sensitivity study was performed. Hoop strains were measured as 17 ± 3% after 10 min relaxation and 21 ± 4% after 20-25 min relaxation, with no significant difference between the two measurements. The extrafibrillar matrix functional moduli were measured as 1.5 ± 0.7 MPa. Fibre-related material parameters showed large variability, with a variance above 0.28. Direct-controlled calibration and validation provides confidence that the model development methodology can capture the measurable variation within the population of tested specimens.

Hybrid Rigid-Deformable Model for Prediction of Neighboring Intervertebral Disk Loads During Flexion Movement After Lumbar Interbody Fusion at L3-4 Level

Knowledge of spinal loads in neighboring disks after interbody fusion plays an important role in the clinical decision of this treatment as well as in the elucidation of its effect. However, controversial findings are still noted in the literature. Moreover, there are no existing models for efficient prediction of intervertebral disk stresses within annulus fibrosus (AF) and nucleus pulposus (NP) regions. In this present study, a new hybrid rigid-deformable modeling workflow was established to quantify the mechanical stress behaviors within AF and NP regions of the L1-2, L2-3, and L4-5 disks after interbody fusion at L3-4 level. The changes in spinal loads were compared with results of the intact model without interbody fusion. The fusion outcomes revealed maximal stress changes (10%) in AF region of L1-2 disk and in NP region of L2-3 disk. The minimal stress change (1%) is noted at the NP region of the L1-2 disk. The validation of simulation outcomes of fused and intact lumbar spine models against those of other computational models and in vivo measurements showed good agreements. Thus, this present study may be used as a novel design guideline for a specific implant and surgical scenario of the lumbar spine disorders.

Analysis of Uncertainty and Variability in Finite Element Computational Models for Biomedical Engineering: Characterization and Propagation

Computational modeling has become a powerful tool in biomedical engineering thanks to its potential to simulate coupled systems. However, real parameters are usually not accurately known, and variability is inherent in living organisms. To cope with this, probabilistic tools, statistical analysis and stochastic approaches have been used. This article aims to review the analysis of uncertainty and variability in the context of finite element modeling in biomedical engineering. Characterization techniques and propagation methods are presented, as well as examples of their applications in biomedical finite element simulations. Uncertainty propagation methods, both non-intrusive and intrusive, are described. Finally, pros and cons of the different approaches and their use in the scientific community are presented. This leads us to identify future directions for research and methodological development of uncertainty modeling in biomedical engineering.

Moderately degenerated lumbar motion segments: Are they truly unstable?

The two main load bearing tissues of the intervertebral disc are the nucleus pulposus and the annulus fibrosus. Both tissues are composed of the same basic components, but differ in their organization and relative amounts. With degeneration, the clear distinction between the two tissues disappears. The changes in biochemical content lead to changes in mechanical behaviour of the intervertebral disc. The aim of the current study was to investigate if well-documented moderate degeneration at the biochemical and fibre structure level leads to instability of the lumbar spine. By taking into account biochemical and ultrastructural changes to the extracellular matrix of degenerating discs, a set of constitutive material parameters were determined that described the individual tissue behaviour. These tissue biomechanical models were then used to simulate dynamic behaviour of the degenerated spinal motion segment, which showed instability in axial rotation, while a stabilizing effect in the other two principle bending directions. When a shear load was applied to the degenerated spinal motion segment, no sign of instability was found. This study found that reported changes to the nucleus pulposus and annulus fibrosus matrix during moderate degeneration lead to a more stable spinal motion segment and that such biomechanical considerations should be incorporated into the general pathophysiological understanding of disc degeneration and how its progress could affect low back pain and its treatments thereof.

Musculoskeletal Modeling of the Lumbar Spine to Explore Functional Interactions between Back Muscle Loads and Intervertebral Disk Multiphysics

During daily activities, complex biomechanical interactions influence the biophysical regulation of intervertebral disks (IVDs), and transfers of mechanical loads are largely controlled by the stabilizing action of spine muscles. Muscle and other internal forces cannot be easily measured directly in the lumbar spine. Hence, biomechanical models are important tools for the evaluation of the loads in those tissues involved in low-back disorders. Muscle force estimations in most musculoskeletal models mainly rely, however, on inverse calculations and static optimizations that limit the predictive power of the numerical calculations. In order to contribute to the development of predictive systems, we coupled a predictive muscle model with the passive resistance of the spine tissues, in a L3-S1 musculoskeletal finite element model with osmo-poromechanical IVD descriptions. The model included 46 fascicles of the major back muscles that act on the lower spine. The muscle model interacted with activity-related loads imposed to the osteoligamentous structure, as standing position and night rest were simulated through distributed upper body mass and free IVD swelling, respectively. Calculations led to intradiscal pressure values within ranges of values measured in vivo. Disk swelling led to muscle activation and muscle force distributions that seemed particularly appropriate to counterbalance the anterior body mass effect in standing. Our simulations pointed out a likely existence of a functional balance between stretch-induced muscle activation and IVD multiphysics toward improved mechanical stability of the lumbar spine understanding. This balance suggests that proper night rest contributes to mechanically strengthen the spine during day activity.