Sarcomere length organization as a design for cooperative function amongst all lumbar spine muscles

The functional design of spine muscles in part dictates their role in moving, loading, and stabilizing the lumbar spine. There have been numerous studies that have examined the isolated properties of these individual muscles. Understanding how these muscles interact and work together, necessary for the prediction of muscle function, spine loading, and stability, is lacking. The objective of this study was to measure sarcomere lengths of lumbar muscles in a neutral cadaveric position and predict the sarcomere operating ranges of these muscles throughout full ranges of spine movements. Sarcomere lengths of seven lumbar muscles in each of seven cadaveric donors were measured using laser diffraction. Using published anatomical coordinate data, superior muscle attachment sites were rotated about each intervertebral joint and the total change in muscle length was used to predict sarcomere length operating ranges. The extensor muscles had short sarcomere lengths in a neutral spine posture and there were no statistically significant differences between extensor muscles. The quadratus lumborum was the only muscle with sarcomere lengths that were optimal for force production in a neutral spine position, and the psoas muscles had the longest lengths in this position. During modeled flexion the extensor, quadratus lumborum, and intertransversarii muscles lengthened so that all muscles operated in the approximate same location on the descending limb of the force-length relationship. The intrinsic properties of lumbar muscles are designed to complement each other. The extensor muscles are all designed to produce maximum force in a mid-flexed posture, and all muscles are designed to operate at similar locations of the force-length relationship at full spine flexion.     

A 3D model of the Achilles tendon to determine the mechanisms underlying nonuniform tendon displacements

The Achilles is the thickest tendon in the body and is the primary elastic energy-storing component during running. The form and function of the human Achilles is complex: twisted structure, intratendinous interactions, and differential motor control from the triceps surae muscles make Achilles behavior difficult to intuit. Recent in vivo imaging of the Achilles has revealed nonuniform displacement patterns that are not fully understood and may result from complex architecture and musculotendon interactions. In order to understand which features of the Achilles tendon give rise to the nonuniform deformations observed in vivo, we used computational modeling to predict the mechanical contributions from different features of the tendon. The aims of this study are to: (i) build a novel computational model of the Achilles tendon based on ultrashort echo time MRI, (ii) compare simulated displacements with published in vivo ultrasound measures of displacement, and (iii) use the model to elucidate the effects of tendon twisting, intratendon sliding, retrocalcaneal insertion, and differential muscle forces on tendon deformation. Intratendon sliding and differential muscle forces were found to be the largest factors contributing to displacement nonuniformity between tendon regions. Elimination of intratendon sliding or muscle forces reduced displacement nonuniformity by 96% and 85%, respectively, while elimination of tendon twist and the retrocalcaneal insertion reduced displacement nonuniformity by only 35% and 3%. These results suggest that changes in the complex internal structure of the tendon alter the interaction between muscle forces and tendon behavior and therefore may have important implications on muscle function during movement.     

A multisegmental cross-bridge kinetics model of the myofibril

Striated muscle is a mechanical system that develops force and generates power in serving vital activities in the body. Striated muscle is a complex biological system; a single mammalian muscle fibre contains up to hundred or even more myofibrils in parallel connected via an inter-myofibril filament network. In one single myofibril thousands of sarcomeres are lined up as a series of linear motors. We recently demonstrated that half-sarcomeres (hS) in a single myofibril operate non-uniformly. We outline a mathematical framework based on cross-bridge kinetics for the simulation of the force response and length change of individual hS in a myofibril. The model describes the muscle myofibril in contraction experiments under various conditions. The myofibril is modeled as a multisegmental mechanical system of hS models, which have active and viscoelastic properties. In the first approach, a two-state cross-bridge formalism relates the hS force to the chemical kinetics of ATP hydrolysis, as first described by Huxley [1957. Muscle structure and theories of contraction. Prog. Biophys. Mol. Biol. 7, 255-318]. Two possible types of biological variability are introduced and modeled. Numerical simulations of a myofibril composed of four to eight hS show a non-uniform hS length distribution and complex internal dynamics upon activation. We demonstrate that the steady-state approximation holds only in restricted time zones during activation. Simulations of myofibril contraction experiments that reproduce the classic steady-state force-length and force-velocity relationships, strictly constrained or “clamped” in either end-held isometric or isotonic contraction conditions, reveal a small but conspicuous effect of hS dynamics on force.     

Twente spine model: A complete and coherent dataset for musculo-skeletal modeling of the lumbar region of the human spine

Musculo-skeletal modeling can greatly help in understanding normal and pathological functioning of the spine. For such models to produce reliable muscle and joint force estimations, an adequate set of musculo-skeletal data is necessary. In this study, we present a complete and coherent dataset for the lumbar spine, based on medical images and dissection measurements from one embalmed human cadaver. We divided muscles into muscle-tendon elements, digitized their attachments at the bones and measured morphological parameters. In total, we measured 11 muscles from one body side, using 96 elements. For every muscle element, we measured three-dimensional coordinates of its attachments, fiber length, tendon length, sarcomere length, optimal fiber length, pennation angle, mass, and physiological cross-sectional area together with the geometry of the lumbar spine. Results were consistent with other anatomical studies and included new data for the serratus posterior inferior muscle. The dataset presented in this paper enables a complete and coherent musculo-skeletal model for the lumbar spine and will improve the current state-of-the art in predicting spinal loading.     

Diastolic ventricular aspiration: a mechanism supporting the rapid filling phase of the human ventricles

During the rapid filling phase of the heart cycle, the internal volumes of the two ventricular cavities approximately double, while the intraventricular pressures rise typically only by an amount of less than 1 kPa. Such a small pressure increase cannot be the sole driving mechanism for the large inflow of blood associated with ventricular expansion during this period. Instead, the rapid filling phase is to be interpreted as being mediated primarily by the heart recoiling elastically from its contracted state, causing blood to be aspirated rapidly into the ventricles. In order to study the role of this mechanism, elastic finite element (FE) simulations of ventricular expansion were performed, taking into account the large deformations occurring during this period and the effective compressibility of the myocardium due to intramural fluid flow. Thereby, a realistic three-dimensional geometry derived from magnetic resonance imaging (MRI) measurements of both human ventricles was used. To validate our FE analyses, the results were compared with published measurements relating to the rapid filling phase of the human left ventricle. Our study shows that, under normal physiological conditions, ventricular aspiration plays a key role in the ventricular filling process.     

Development of a statistical shape-function model of the implanted knee for real-time prediction of joint mechanics

Outcomes of total knee arthroplasty (TKA) are dependent on surgical technique, patient variability, and implant design. Non-optimal design or alignment choices may result in undesirable contact mechanics and joint kinematics, including poor joint alignment, instability, and reduced range of motion. Implant design and surgical alignment are modifiable factors with potential to improve patient outcomes, and there is a need for robust implant designs that can accommodate patient variability. Our objective was to develop a statistical shape-function model (SFM) of a posterior stabilized implanted knee to instantaneously predict joint mechanics in an efficient manner. Finite element methods were combined with Latin hypercube sampling and regression analyses to produce modeling equations relating nine implant design and six surgical alignment parameters to tibiofemoral (TF) joint mechanics outcomes during a deep knee bend. A SFM was developed and TF contact mechanics, kinematics, and soft tissue loads were instantaneously predicted from the model. Average normalized root-mean-square error predictions were between 2.79% and 9.42%, depending on the number of parameters included in the model. The statistical shape-function model generated instantaneous joint mechanics predictions using a maximum of 130 training simulations, making it ideally suited for integration into a patient-specific design and alignment optimization pipeline. Such a tool may be used to optimize kinematic function to achieve more natural motion or minimize implant wear, and may aid the engineering and clinical communities in improving patient satisfaction and surgical outcomes.     

Passive diastolic modelling of human ventricles: Effects of base movement and geometrical heterogeneity

Left-ventricular (LV) remodelling, associated with diastolic heart failure, is driven by an increase in myocardial stress. Therefore, normalisation of LV wall stress is the cornerstone of many therapeutic treatments. However, information regarding such regional stress–strain for human LV is still limited. Thus, the objectives of our study were to determine local diastolic stress–strain field in healthy LVs, and consequently, to identify the regional variations amongst them due to geometric heterogeneity. Effects of LV base movement on diastolic model predictions, which were ignored in the literature, were further explored. Personalised finite-element modelling of five normal human bi-ventricles was carried out using subject-specific myocardium properties. Model prediction was validated individually through comparison with end-diastolic volume and a new shape-volume based measurement of LV cavity, extracted from magnetic resonance imaging. Results indicated that incorporation of LV base movement improved the model predictions (shape-volume relevancy of LV cavity), and therefore, it should be considered in future studies. The LV endocardium always experienced higher fibre stress compared to the epicardium for all five subjects. The LV wall near base experienced higher stress compared to equatorial and apical locations. The lateral LV wall underwent greater stress distribution (fibre and sheet stress) compared to other three regions. In addition, normal ranges of different stress–strain components in different regions of LV wall were reported for five healthy ventricles. This information could be used as targets for future computational studies to optimise diastolic heart failure treatments or design new therapeutic interventions/devices.     

A finite element study relating to the rapid filling phase of the human ventricles

During the rapid diastolic filling phase at rest, the ventricles of the human heart double approximately in volume. In order to investigate whether the ventricular filling pressures measured under physiological conditions can give rise to such an extensive augmentation in ventricular volumes, a finite element model of the human right and left ventricles has been developed, taking into account the nonlinear mechanical behavior and effective compressibility of the myocardial tissue. The results were compared with the filling phase of the human left ventricle as extrapolated from measurements documented in the literature. We arrived at the conclusion that the ventricular pressures measured during the rapid filling phase cannot be the sole cause of the rise of the observed ventricular volumes. We rather advocate the assumption that further dilating mechanisms might be part of ventricular activity thus heralding a multiple function of the ventricular muscle body. A further result indicates that under normal conditions the influence of the viscoelasticity of the tissue should not be disregarded in ventricular mechanics.     

Influence of viscosity on myocardium mechanical activity: a mathematical model

We have previously proposed and validated a mathematical model of myocardium contraction-relaxation cycle based on current knowledge of regulatory role of Ca2+ and cross-bridge kinetics in cardiac cell. That model did not include viscous elements. Here we propose a modification of the model, in which two viscous elements are added, one in parallel to the contractile element, and one more in parallel to the series elastic element. The modified model allowed us to simulate and explain some subtle experimental data on relaxation velocity in isotonic twitches and on a mismatch between the time course of sarcomere shortening/lengthening and the time course of active force generation in isometric twitches. Model results were compared with experimental data obtained from 28 rat LV papillary muscles contracting and relaxing against various loads. Additional model analysis suggested contribution of viscosity to main inotropic and lusitropic characteristics of myocardium performance.     

Twente spine model: A complete and coherent dataset for musculo-skeletal modeling of the thoracic and cervical regions of the human spine

Musculo-skeletal modeling could play a key role in advancing our understanding of the healthy and pathological spine, but the credibility of such models are strictly dependent on the accuracy of the anatomical data incorporated. In this study, we present a complete and coherent musculo-skeletal dataset for the thoracic and cervical regions of the human spine, obtained through detailed dissection of an embalmed male cadaver. We divided the muscles into a number of muscle-tendon elements, digitized their attachments at the bones, and measured morphological muscle parameters. In total, 225 muscle elements were measured over 39 muscles. For every muscle element, we provide the coordinates of its attachments, fiber length, tendon length, sarcomere length, optimal fiber length, pennation angle, mass, and physiological cross-sectional area together with the skeletal geometry of the cadaver. Results were consistent with similar anatomical studies. Furthermore, we report new data for several muscles such as rotatores, multifidus, levatores costarum, spinalis, semispinalis, subcostales, transversus thoracis, and intercostales muscles. This dataset complements our previous study where we presented a consistent dataset for the lumbar region of the spine (Bayoglu et al., 2017). Therefore, when used together, these datasets enable a complete and coherent dataset for the entire spine. The complete dataset will be used to develop a musculo-skeletal model for the entire human spine to study clinical and ergonomic applications.     

A new cell-based FE model for the mechanics of embryonic epithelia

In order to overcome a significant stiffening artefact associated with current finite element (FE) models for the mechanics of embryonic epithelia, two new FE formulations were developed. Cell–cell interfacial tensions γ are represented by constant-force rod elements as in previous models. However, the viscosity of the cytoplasm with its embedded organelles and filament networks is modeled using viscous triangular elements, it is modeled using either radial and circumferential dashpots or an orthogonal dashpot system rather than the viscous triangular elements typical of previous two-dimensional FE models. The models are tested against tissue (epithelium) stretching because it gives rise to significant changes in cell shape and against cell sorting because it involves high rates of cell rearrangement. The orthogonal dashpot system is found to capture cell size and shape effects well, give the model cells characteristics that are consistent with those of real cells, provide high computational efficiency and hold promise for future three-dimensional analyses.     

A study of the viscoelastic effect in a bone remodeling model

The present paper addresses the following question can a simple regulatory bone remodeling model predict effects of viscosity on the trabecular morphology? For that, we propose an extension of a previous bone remodeling model by taking into account the viscosity properties of the tissue. Zener’s law is used to describe the mechanical behavior of the bone and a specific law of the apparent bone density rate is proposed. Based on stability analysis, numerical simulations are then performed to investigate the viscosity role on simulations of the bone remodeling process. We show that the viscous contribution affects the evolution of the apparent bone density, by slowing down the adaptation process, which seems to be confirmed by simulations with real data obtained from rat tibia.     

Strains at the myotendinous junction predicted by a micromechanical model

The goal of this work was to create a finite element micromechanical model of the myotendinous junction (MTJ) to examine how the structure and mechanics of the MTJ affect the local micro-scale strains experienced by muscle fibers. We validated the model through comparisons with histological longitudinal sections of muscles fixed in slack and stretched positions. The model predicted deformations of the A-bands within the fiber near the MTJ that were similar to those measured from the histological sections. We then used the model to predict the dependence of local fiber strains on activation and the mechanical properties of the endomysium. The model predicted that peak micro-scale strains increase with activation and as the compliance of the endomysium decreases. Analysis of the models revealed that, in passive stretch, local fiber strains are governed by the difference of the mechanical properties between the fibers and the endomysium. In active stretch, strain distributions are governed by the difference in cross-sectional area along the length of the tapered region of the fiber near the MTJ. The endomysium provides passive resistance that balances the active forces and prevents the tapered region of the fiber from undergoing excessive strain. These model predictions lead to the following hypotheses: (i) the increased likelihood of injury during active lengthening of muscle fibers may be due to the increase in peak strain with activation and (ii) endomysium may play a role in protecting fibers from injury by reducing the strains within the fiber at the MTJ.     

3D finite element model of the diabetic neuropathic foot: A gait analysis driven approach

Diabetic foot is an invalidating complication of diabetes that can lead to foot ulcers. Three-dimensional (3D) finite element analysis (FEA) allows characterizing the loads developed in the different anatomical structures of the foot in dynamic conditions. The aim of this study was to develop a subject specific 3D foot FE model (FEM) of a diabetic neuropathic (DNS) and a healthy (HS) subject, whose subject specificity can be found in term of foot geometry and boundary conditions. Kinematics, kinetics and plantar pressure (PP) data were extracted from the gait analysis trials of the two subjects with this purpose. The FEM were developed segmenting bones, cartilage and skin from MRI and drawing a horizontal plate as ground support. Materials properties were adopted from previous literature. FE simulations were run with the kinematics and kinetics data of four different phases of the stance phase of gait (heel strike, loading response, midstance and push off). FEMs were then driven by group gait data of 10 neuropathic and 10 healthy subjects. Model validation focused on agreement between FEM-simulated and experimental PP.     

Implementation of reflex loops in a biomechanical finite element model

In the field of biomechanics, the offer of models which are more and more realistic requires to integrate a physiological response, in particular, the controlled muscle bracing and the reflexes. The following work aims to suggest a unique methodology which couples together a sensory and motor loop with a finite element model. Our method is applied to the study of the oscillation of the elbow in the case of a biceps brachial stretch reflex. The results obtained are promising in the purpose of the development of reactive human body models.     

3D non-linear analysis of the acetabular construct following impaction grafting

The study investigates the short-term behaviour of the acetabular construct following revision hip arthroplasty, carried out using the Slooff–Ling impaction grafting technique; using 3D finite element analyses. An elasto-plastic material model is used to describe the constitutive behaviour of morsellised cortico-cancellous bone (MCB) graft, since it has been shown that MCB undergoes significant plastic deformation under normal physiological loads. Based on previous experimental studies carried out by the authors and others, MCB is modelled using non-linear elasticity and Drucker Prager Cap (DPC) plasticity. Loading associated with walking, sitting down, and standing up is applied to the acetabular cup through a femoral head using smooth sliding surfaces. The analyses yield distinctive patterns of migration and rotation due to different activities. These are found to be similar to those observed in the clinical setting.