Mechanically simulated muscle forces strongly stabilize intact and injured upper cervical spine specimens

Although muscles are assumed to be capable of stabilizing the spinal column in vivo, they have only rarely been simulated in vitro. Their effect might be of particular importance in unstable segments. The present study therefore tests the hypothesis that mechanically simulated muscle forces stabilize intact and injured cervical spine specimens. In the first step, six human occipito-cervical spine specimens were loaded intact in a spine tester with pure moments in lateral bending (+/- 1.5 N m), flexion-extension (+/- 1.5 N m) and axial rotation (+/- 0.5 N m). In the second step, identical flexibility tests were carried out during constant traction of three mechanically simulated muscle pairs: splenius capitits (5 N), semispinalis capitis (5 N) and longus colli (15 N). Both steps were repeated after unilateral and bilateral transection of the alar ligaments. The muscle forces strongly stabilized C0-C2 in all loading and injury states. This was most obvious in axial rotation, where a reduction of range of motion (ROM) and neutral zone to <50% (without muscles=100%) was observed. With increasing injury the normalized ROM (intact condition=100%) increased with and without muscles approximately to the same extend. With bilateral injury this increase was 125-132% in lateral bending, 112%-119% in flexion-extension and 103-116% in axial rotation. Mechanically simulated cervical spine muscles strongly stabilized intact and injured cervical spine specimens. Nevertheless, it could be shown that in vitro flexibility tests without muscle force simulation do not necessarily lead to an overestimation of spinal instability if the results are normalized to the intact state.     

Sensitivity of muscle and intervertebral disc force computations to variations in muscle attachment sites

Abstract

The current paper aims at assessing the sensitivity of muscle and intervertebral disc force computations against potential errors in modeling muscle attachment sites. We perturbed each attachment location in a complete and coherent musculoskeletal model of the human spine and quantified the changes in muscle and disc forces during standing upright, flexion, lateral bending, and axial rotation of the trunk. Although the majority of the muscles caused minor changes (less than 5%) in the disc forces, certain muscle groups, for example, quadratus lumborum, altered the shear and compressive forces as high as 353% and 17%, respectively. Furthermore, percent changes were higher in the shear forces than in the compressive forces. Our analyses identified certain muscles in the rib cage (intercostales interni and intercostales externi) and lumbar spine (quadratus lumborum and longissimus thoracis) as being more influential for computing muscle and disc forces. Furthermore, the disc forces at the L4/L5 joint were the most sensitive against muscle attachment sites, followed by T6/T7 and T12/L1 joints. Presented findings suggest that modeling muscle attachment sites based on solely anatomical illustrations might lead to erroneous evaluation of internal forces and promote using anatomical datasets where these locations were accurately measured. When developing a personalized model of the spine, certain care should also be paid especially for the muscles indicated in this work.     

Optimised loads for the simulation of axial rotation in the lumbar spine

Simplified loading modes (pure moment, compressive force) are usually applied in the in vitro studies to simulate flexion-extension, lateral bending and axial rotation of the spine. The load magnitudes for axial rotation vary strongly in the literature. Therefore, the results of current investigations, e.g. intervertebral rotations, are hardly comparable and may involve unrealistic values. Thus, the question 'which in vitro applicable loading mode is the most realistic' remains open. A validated finite element model of the lumbar spine was employed in two sensitivity studies to estimate the ranges of results due to published load assumptions and to determine the input parameters (e.g. torsional moment), which mostly affect the spinal load and kinematics during axial rotation. In a subsequent optimisation study, the in vitro applicable loading mode was determined, which delivers results that fit best with available in vivo measurements. The calculated results varied widely for loads used in the literature with potential high deviations from in vivo measured values. The intradiscal pressure is mainly affected by the magnitude of the compressive force, while the torsional moment influences mainly the intervertebral rotations and facet joint forces. The best agreement with results measured in vivo were found for a compressive follower force of 720N and a pure moment of 5.5Nm applied to the unconstrained vertebra L1. The results reveal that in many studies the assumed loads do not realistically simulate axial rotation. The in vitro applicable simplified loads cannot perfectly mimic the in vivo situation. However, the optimised values lead to the best agreement with in vivo measured values. Their consequent application would lead to a better comparability of different investigations.     

Effect of muscle wrapping on model estimates of neck muscle strength

Of the computational models of the cervical spine reported in the literature, not one takes into account the changes in muscle paths due to the underlying vertebrae. Instead, all model the individual muscle paths as straight-line segments. The major aim of this study was to quantify the changes in muscle moment arm, muscle force and joint moment due to muscle wrapping in the cervical spine. Five muscles in a straight-line model of the cervical spine were wrapped around underlying vertebrae, and the results obtained from this model were compared against the original. The two models were then validated against experimental and computational data. Results show that muscle wrapping has a significant effect on muscle moment arms and therefore joint moments and should not be neglected.     

Effect of muscle wrapping on model estimates of neck muscle strength

Of the computational models of the cervical spine reported in the literature, not one takes into account the changes in muscle paths due to the underlying vertebrae. Instead, all model the individual muscle paths as straight-line segments. The major aim of this study was to quantify the changes in muscle moment arm, muscle force and joint moment due to muscle wrapping in the cervical spine. Five muscles in a straight-line model of the cervical spine were wrapped around underlying vertebrae, and the results obtained from this model were compared against the original. The two models were then validated against experimental and computational data. Results show that muscle wrapping has a significant effect on muscle moment arms and therefore joint moments and should not be neglected.     

Finite element analysis of moment-rotation relationships for human cervical spine

A comprehensive, geometrically accurate, nonlinear C0-C7 FE model of head and cervical spine based on the actual geometry of a human cadaver specimen was developed. The motions of each cervical vertebral level under pure moment loading of 1.0 Nm applied incrementally on the skull to simulate the movements of the head and cervical spine under flexion, tension, axial rotation and lateral bending with the inferior surface of the C7 vertebral body fully constrained were analysed. The predicted range of motion (ROM) for each motion segment were computed and compared with published experimental data. The model predicted the nonlinear moment-rotation relationship of human cervical spine. Under the same loading magnitude, the model predicted the largest rotation in extension, followed by flexion and axial rotation, and least ROM in lateral bending. The upper cervical spines are more flexible than the lower cervical levels. The motions of the two uppermost motion segments account for half (or even higher) of the whole cervical spine motion under rotational loadings. The differences in the ROMs among the lower cervical spines (C3-C7) were relatively small. The FE predicted segmental motions effectively reflect the behavior of human cervical spine and were in agreement with the experimental data. The C0-C7 FE model offers potentials for biomedical and injury studies.     

Validation of lumbar spine loading from a musculoskeletal model including the lower limbs and lumbar spine

Low back mechanics are important to quantify to study injury, pain and disability. As in vivo forces are difficult to measure directly, modeling approaches are commonly used to estimate these forces. Validation of model estimates is critical to gain confidence in modeling results across populations of interest, such as people with lower-limb amputation. Motion capture, ground reaction force and electromyographic data were collected from ten participants without an amputation (five male/five female) and five participants with a unilateral transtibial amputation (four male/one female) during trunk-pelvis range of motion trials in flexion/extension, lateral bending and axial rotation. A musculoskeletal model with a detailed lumbar spine and the legs including 294 muscles was used to predict L4-L5 loading and muscle activations using static optimization. Model estimates of L4-L5 intervertebral joint loading were compared to measured intradiscal pressures from the literature and muscle activations were compared to electromyographic signals. Model loading estimates were only significantly different from experimental measurements during trunk extension for males without an amputation and for people with an amputation, which may suggest a greater portion of L4-L5 axial load transfer through the facet joints, as facet loads are not captured by intradiscal pressure transducers. Pressure estimates between the model and previous work were not significantly different for flexion, lateral bending or axial rotation. Timing of model-estimated muscle activations compared well with electromyographic activity of the lumbar paraspinals and upper erector spinae. Validated estimates of low back loading can increase the applicability of musculoskeletal models to clinical diagnosis and treatment.     

A myoelectrically based dynamic three-dimensional model to predict loads on lumbar spine tissues during lateral bending.

This work describes a dynamic model of the low back that incorporates extensive anatomical detail of the musculo-ligamentous-skeletal system to predict the load time histories of individual tissues. The dynamic reaction moment about L4/L5 was determined during lateral bending from a linked-segment model. This reaction moment was partitioned into restorative components provided by the disc, ligament strain, and active-muscle contraction using a second model of the spine that incorporated a detailed representation of the anatomy. Muscle contraction forces were estimated using both information from surface electromyograms, collected from 12 sites, and consideration of the modulating effects of muscle length, cross-sectional area and passive elasticity. This modelling technique is sensitive to the different ways in which individuals recruit their musculature to satisfy moment constraints. Time histories of muscle forces are provided. High muscle loads are consistent with the common clinical observation of muscle strain often produced by load handling. Furthermore, the coactivation measured in muscles on both sides of the trunk suggests that muscles are recruited to satisfy the lateral bending reaction torque in addition to performing other mechanical roles such as spine stabilization. If an estimate of the intervertebral joint compression is desired for assessment of lateral bends in industry, then a single equivalent lateral muscle with a moment arm of approximately 3.0-4.0 cm would conservatively capture the effects of muscle co-contraction quantified in this study.     

Muscle moment arm and normalized moment contributions as reference data for musculoskeletal elbow and wrist joint models

A geometric musculoskeletal model of the elbow and wrist joints was developed to calculate muscle moment arms throughout elbow flexion/extension, forearm pronation/supination, wrist flexion/extension and radial/ulnar deviation. Model moment arms were verified with data from cadaver specimen studies and geometric models available in the literature. Coefficients of polynomial equations were calculated for all moment arms as functions of joint angle, with special consideration to coupled muscles as a function of two joint angles. Additionally, a “normalized potential moment (NPM)” contribution index for each muscle across the elbow and wrist joints in four degrees-of-freedom was determined using each muscle's normalized physiological cross-sectional area (PCSA) and peak moment arm (MA). We hypothesize that (a) a geometric model of the elbow and wrist joints can represent the major attributes of MA versus joint angle from many literature sources of cadaver and model data and (b) an index can represent each muscle's normalized moment contribution to each degree-of-freedom at the elbow and wrist. We believe these data serve as a simple, yet comprehensive, reference for how the primary 16 muscles across the elbow and wrist contribute to joint moment and overall joint performance.     

A three-dimensional model of the human cervical spine for impact simulation

A three-dimensional analytical model of the cervical spine is described. The cervical vertebrae and the head are modeled as rigid bodies which are interconnected by deformable elements representing the intervertebral disks, facet joints, ligaments and muscles. A special pentahedral continuum element for representing the articular facets is described which effectively maintains stability of the cervical spine in both lateral and frontal plane accelerations, which is very difficult with multi-spring models of the facets. A simplified representation is used for the spine and body below the level of T1. The neck musculature is modeled by over 100 muscle elements representing 22 major muscle groups in the neck. The model has been validated for frontal and sideways impact accelerations by simulating published experimental data. Results are also presented to show the effects of the stretch reflex response on the dynamics of the head and neck under moderate acceleration.     

Abdominal muscles dominate contributions to vertebral joint stiffness during the push-up

The goal of this study was to quantify the relative contributions of each muscle group surrounding the spine to vertebral joint rotational stiffness (VJRS) during the push-up exercise. Upper-body kinematics, three-dimensional hand forces and lumbar spine postures, and 14 channels (bilaterally from rectus abdominis, external oblique, internal oblique, latissimus dorsi, thoracic erector spinae, lumbar erector spinae, and multifidus) of trunk electromyographic (EMG) activity were collected from 11 males and used as inputs to a biomechanical model that determined the individual contributions of 10 muscle groups surrounding the lumbar spine to VJRS at five lumbar vertebral joints (L1-L2 to L5-S1). On average, the abdominal muscles contributed 64.32 +/- 8.50%, 86.55 +/- 1.13%, and 83.84 +/- 1.95% to VJRS about the flexion/extension, lateral bend, and axial twist axes, respectively. Rectus abdominis contributed 43.16 +/- 3.44% to VJRS about the flexion/extension axis at each lumbar joint, and external oblique and internal oblique, respectively contributed 52.61 +/- 7.73% and 62.13 +/- 8.71% to VJRS about the lateral bend and axial twist axes, respectively, at all lumbar joints with the exception of L5-S1. Owing to changes in moment arm length, the external oblique and internal oblique, respectively contributed 55.89% and 50.01% to VJRS about the axial twist and lateral bend axes at L5-S1. Transversus abdominis, multifidus, and the spine extensors contributed minimally to VJRS during the push-up exercise. The push-up challenges the abdominal musculature to maintain VJRS. The orientation of the abdominal muscles suggests that each muscle primarily controls the rotational stiffness about a single axis.     

Static and dynamic bending responses of the human cervical spine.

The quasi-static and dynamic bending responses of the human mid-lower cervical spine were determined using cadaver intervertebral joints fixed at the base to a six-axis load cell. Flexion bending moment was applied to the superior end of the specimen using an electrohydraulic piston. Each specimen was tested under three cycles of quasi-static load-unload and one high-speed dynamic load. A total of five specimens were included in this study. The maximum intervertebral rotation ranged from 11.0 to 15.4 deg for quasi-static tests and from 22.9 to 34.4 deg for dynamic tests. The resulting peak moments at the center of the intervertebral joint ranged from 3.8 to 6.9 Nm for quasi-static tests and from 14.0 to 31.8 Nm for dynamic tests. The quasi-static stiffness ranged from 0.80 to 1.35 Nm/deg with a mean of 1.03 Nm/deg (+/- 0.11 Nm/deg). The dynamic stiffness ranged from 1.08 to 2.00 Nm/deg with a mean of 1.50 Nm/deg (+/- 0.17 Nm/deg). The differences between the two stiffnesses were statistically significant (p < 0.01). Exponential functions were derived to describe the quasi-static and dynamic moment-rotation responses. These results provide input data for lumped-parameter models and validation data for finite element models to better investigate the biomechanics of the human cervical spine.     

Moment-rotation relationships of the ligamentous occipito-atlanto-axial complex

The relationships between applied pure moments at the occiput (C0) and the resulting rotations at the atlanto-occipital (C0-C1) and atlanto-axial (C1-C2) joints are quantified. In axial twist, with a moment of 0.3 Nm, a mean rotation of about 2.5 degrees and 23.3 degrees was observed at C0-C1 and C1-C2 units respectively. Both the atlas and axis contributed to produce lateral bending motion. The ratio between extension and flexion rotations at C0-C1 was 2.5:1. Lateral bending and axial rotations were strongly coupled to each other. The occipito-atlanto-axial complex exhibited a large 'neutral zone' compared to lower cervical spine segments. The likely clinical significance of these findings are discussed.     

The effect of follower load on the intersegmental coupled motion characteristics of the human thoracic spine: An in vitro study using entire rib cage specimens

The mechanical coupling behaviour of the thoracic spine is still not fully understood. For the validation of numerical models of the thoracic spine, however, the coupled motions within the single spinal segments are of importance to achieve high model accuracy. In the present study, eight fresh frozen human thoracic spinal specimens (C7-L1, mean age 54 ± 6 years) including the intact rib cage were loaded with pure bending moments of 5 Nm in flexion/extension (FE), lateral bending (LB), and axial rotation (AR) with and without a follower load of 400 N. During loading, the relative motions of each vertebra were monitored. Follower load decreased the overall ROM (T1-T12) significantly (p < 0.01) in all primary motion directions (extension: −46%, left LB: −72%, right LB: −72%, left AR: −26%, right AR: −26%) except flexion (−36%). Substantial coupled motion was found in lateral bending with ipsilateral axial rotation, which increased after a follower load was applied, leading to a dominant axial rotation during primary lateral bending, while all other coupled motions in the different motion directions were reduced under follower load. On the monosegmental level, the follower load especially reduced the ROM of the upper thoracic spine from T1-T2 to T4-T5 in all motion directions and the ROM of the lower thoracic spine from T9-T10 to T11-T12 in primary lateral bending. The facet joints, intervertebral disc morphologies, and the sagittal curvature presumably affect the thoracic spinal coupled motions depending on axial compressive preloading. Using these results, the validation of numerical models can be performed more accurately.     

Predicted pattern of human muscle activity during clenching derived from a computer assisted model: symmetric vertical bite forces

A computer assisted three-dimensional model of the jaw, based on linear programming, is presented. The upper and lower attachments of the muscles of mastication have been measured on a single human skull and divided into thirteen independent units on each side--a total of 26 muscle elements. The direction (in three dimensions) and maximum forces that could be developed by each muscle element, the bite reaction and two joint reactions are included in the model. It is shown for symmetrical biting that a model which minimizes the sum of the muscle forces used to produce a given bite force activates muscles in a way which corresponds well with previous observations on human subjects. A model which minimizes the joint reactions behaves differently and is rejected. An analysis of the way the chosen model operates suggests that there are two types of jaw muscles, power muscles and control muscles. Power muscles (superficial masseter, medial pterygoid and some of temporalis) produce the bite force but tend to displace the condyle up or down the articular eminence. This displacement is prevented by control muscles (oblique temporalis and lateral pterygoid) which have very poor moment arms for generating usual bite forces, but are efficient for preventing condylar slide. The model incorporates the concept that muscles consist of elements which can contract independently. It predicts that those muscle elements with longer moment arms relative to the joint are the first to be activated and, as the bite force increases, a ripple of activity spreads into elements with shorter moment arms. In general, the model can be used to study the three-dimensional activity in any system of joints and muscles.     

In vivo determination of the direction of rotation and moment-angle relationship of individual elbow muscles.

The direction of rotation (DOR) of individual elbow muscles, defined as the direction in which a muscle rotates the forearm relative to the upper arm in three-dimensional space, was studied in vivo as a function of elbow flexion and forearm rotation. Electrical stimulation was used to activate an individual muscle selectively, and the resultant flexion-extension, supination-pronation, and varus-valgus moments were used to determine the DOR. Furthermore, multi-axis moment-angle relationships of individual muscles were determined by stimulating the muscle at a constant submaximal level across different joint positions, which was assumed to result in a constant level of muscle activation. The muscles generate significant moments about axes other than flexion-extension, which is potentially important for actively controlling joint movement and maintaining stability about all axes. Both the muscle DOR and the multi axis moments vary with the joint position systematically. Variations of the DOR and moment-angle relationship across muscle twitches of different amplitudes in a subject were small, while there were considerable variations between subjects.     

Artificial Cervical Vertebra and Intervertebral Complex Replacement through the Anterior Approach in Animal Model: A Biomechanical and In Vivo Evaluation of a Successful Goat Model

This was an in vitro and in vivo study to develop a novel artificial cervical vertebra and intervertebral complex (ACVC) joint in a goat model to provide a new method for treating degenerative disc disease in the cervical spine. The objectives of this study were to test the safety, validity, and effectiveness of ACVC by goat model and to provide preclinical data for a clinical trial in humans in future. We designed the ACVC based on the radiological and anatomical data on goat and human cervical spines, established an animal model by implanting the ACVC into goat cervical spines in vitro prior to in vivo implantation through the anterior approach, and evaluated clinical, radiological, biomechanical parameters after implantation. The X-ray radiological data revealed similarities between goat and human intervertebral angles at the levels of C2-3, C3-4, and C4-5, and between goat and human lordosis angles at the levels of C3-4 and C4-5. In the in vivo implantation, the goats successfully endured the entire experimental procedure and recovered well after the surgery. The radiological results showed that there was no dislocation of the ACVC and that the ACVC successfully restored the intervertebral disc height after the surgery. The biomechanical data showed that there was no significant difference in range of motion (ROM) or neural zone (NZ) between the control group and the ACVC group in flexion-extension and lateral bending before or after the fatigue test. The ROM and NZ of the ACVC group were greater than those of the control group for rotation. In conclusion, the goat provides an excellent animal model for the biomechanical study of the cervical spine. The ACVC is able to provide instant stability after surgery and to preserve normal motion in the cervical spine.     

Contact pressure in the facet joint during sagittal bending of the cadaveric cervical spine

The facet joint contributes to the normal biomechanical function of the spine by transmitting loads and limiting motions via articular contact. However, little is known about the contact pressure response for this joint. Such information can provide a quantitative measure of the facet joint's local environment. The objective of this study was to measure facet pressure during physiologic bending in the cervical spine, using a joint capsule-sparing technique. Flexion and extension bending moments were applied to six human cadaveric cervical spines. Global motions (C2-T1) were defined using infra-red cameras to track markers on each vertebra. Contact pressure in the C5-C6 facet was also measured using a tip-mounted pressure transducer inserted into the joint space through a hole in the postero-inferior region of the C5 lateral mass. Facet contact pressure increased by 67.6 ± 26.9 kPa under a 2.4 Nm extension moment and decreased by 10.3 ± 9.7 kPa under a 2.7 Nm flexion moment. The mean rotation of the overall cervical specimen motion segments was 9.6 ± 0.8° and was 1.6 ± 0.7° for the C5-C6 joint, respectively, for extension. The change in pressure during extension was linearly related to both the change in moment (51.4 ± 42.6 kPa/Nm) and the change in C5-C6 angle (18.0 ± 108.9 kPa/deg). Contact pressure in the inferior region of the cervical facet joint increases during extension as the articular surfaces come in contact, and decreases in flexion as the joint opens, similar to reports in the lumbar spine despite the difference in facet orientation in those spinal regions. Joint contact pressure is linearly related to both sagittal moment and spinal rotation. Cartilage degeneration and the presence of meniscoids may account for the variation in the pressure profiles measured during physiologic sagittal bending. This study shows that cervical facet contact pressure can be directly measured with minimal disruption to the joint and is the first to provide local pressure values for the cervical joint in a cadaveric model.     

Morphological aspects and biomechanical properties of the vertebroaxial joint (C2-C3).

The C2-C3 intervertebral joint must be regarded as a transitional area situated between the upper cervical spine where most rotation of the neck and little flexion and extension occur and the lower cervical spaces where chiefly motion in the sagittal plane and also somewhat rotation take place. Under normal circumstances the range of flexion-extension reaches 11 degrees, slighter than below (19.5 degrees at C5-C6); on the opposite, the range of rotation attains 7 degrees; less than above but much more than below (0 degrees at C5-C6). The motion in the sagittal and coronal planes is relatively poor because of the location of Penning's motor-axis of C2 which runs far from the vertebral body and the lowness of the intervertebral disc. However, the rotation of C2 with respect to C3 is fair by the peculiar inclination of the articular facets which slope sagittally but also coronally and trace a sphere whereupon C2 may move around its motor-centre in any plane. When C2-C3 is surgically fixed by bone graft, the lack of motion is completed by a "compensatory movement" in the upper cervical spaces and especially at the atlantooccipital joint for flexion-extension. In the same way, C2-C3 may improve its mobility especially in the sagittal plane when the inferior partner is blocked by surgical or arthritic fusion.