Frequency-dependent changes in neuromuscular responses to cyclic lumbar flexion

Repetitive lifting in the workplace has been identified to be a cause of low back disorders. Epidemiologic data further supports an hypothesis that higher repetition rate (i.e. frequency) is an added risk factor. The objective of this study was to provide experimental data testing the above hypothesis. An in vivo feline model was subjected to 20-min of cyclic lumbar loading at frequencies of 0.1 Hz and 0.5 Hz while monitoring the EMG from the L-3/4-L-5/6 multifidus muscles and the creep at the L-4/5 level. Seven hours of rest were allowed after the cyclic flexion/extension was terminated. During this rest period, a single test cycle was performed every hour to assess recovery of EMG and lumbar creep. The results demonstrate that cyclic lumbar flexion elicits a transient neuromuscular disorder consisting of EMG spasms during the cyclic loading and initial and delayed muscular hyperexcitabilities during the rest period. Cyclic loading at 0.5 Hz resulted in significant (p<0.05) increase in the hyperexcitability magnitude and duration during the recovery period. It was concluded that repetitive lumbar loading at fast rates is indeed a risk factor as it induces larger creep in the lumbar viscoelastic tissues which in turn intensify the resulting neuromuscular disorder.     

Use of additive models to represent trends in a barley field trial

Even when a field experiment has been designed with care, subsequent examination of the plot values may reveal additional unforeseen trends. In this paper we examine data from a barley pathology field trial and show that additive models provide a flexible representation of environmental trends, in one or two dimensions. Such models smooth out noise in the observed data, rather than fit an equation specified in advance. This approach tends to increase the precision of treatment comparisons relative to a classical analysis of variance. We recommend the use of residual plots to explore experimental data for underlying trends, and additive models to display these trends and estimate treatment effects.     

Position sense in the lumbar spine with torso flexion and loading

Proprioception plays an important role in appropriate sensation of spine position, movement, and stability. Previous research has demonstrated that position sense error in the lumbar spine is increased in flexed postures. This study investigated the change in position sense as a function of altered trunk flexion and moment loading independently. Reposition sense of lumbar angle in 17 subjects was assessed. Subjects were trained to assume specified lumbar angles using visual feedback. The ability of the subjects to reproduce this curvature without feedback was then assessed. This procedure was repeated for different torso flexion and moment loading conditions. These measurements demonstrated that position sense error increased significantly with the trunk flexion (40%, p < .05) but did not increase with moment load (p = .13). This increased error with flexion suggests a loss in the ability to appropriately sense and therefore control lumbar posture in flexed tasks. This loss in proprioceptive sense could lead to more variable lifting coordination and a loss in dynamic stability that could increase low back injury risk. This research suggests that it is advisable to avoid work in flexed postures.     

Conclusions on motor control depend on the type of model used to represent the periphery

Within the field of motor control, there is no consensus on which kinematic and kinetic aspects of movements are planned or controlled. Perturbing goal-directed movements is a frequently used tool to answer this question. To be able to draw conclusions about motor control from kinematic responses to perturbations, a model of the periphery (i.e., the skeleton, muscle–tendon complexes, and spinal reflex circuitry) is required. The purpose of the present study was to determine to what extent such conclusions depend on the level of simplification with which the dynamical properties of the periphery are modeled. For this purpose, we simulated fast goal-directed single-joint movement with four existing types of models. We tested how three types of perturbations affected movement trajectory if motor commands remained unchanged. We found that the four types of models of the periphery showed different robustness to the perturbations, leading to different predictions on how accurate motor commands need to be, i.e., how accurate the knowledge of external conditions needs to be. This means that when interpreting kinematic responses obtained in perturbation experiments the level of error correction attributed to adaptation of motor commands depends on the type of model used to describe the periphery.     

Validation and discovery from computational biology models

Simulation software is often a fundamental component in systems biology projects and provides a key aspect of the integration of experimental and analytical techniques in the search for greater understanding and prediction of biology at the systems level. It is important that the modelling and analysis software is reliable and that techniques exist for automating the analysis of the vast amounts of data which such simulation environments generate. A rigorous approach to the development of complex modelling software is needed. Such a framework is presented here together with techniques for the automated analysis of such models and a process for the automatic discovery of biological phenomena from large simulation data sets. Illustrations are taken from a major systems biology research project involving the in vitro investigation, modelling and simulation of epithelial tissue.     

Biomechanical implications of lumbar spinal ligament transection

Many lumbar spine surgeries either intentionally or inadvertently damage or transect spinal ligaments. The purpose of this work was to quantify the previously unknown biomechanical consequences of isolated spinal ligament transection on the remaining spinal ligaments (stress transfer), vertebrae (bone remodelling stimulus) and intervertebral discs (disc pressure) of the lumbar spine. A finite element model of the full lumbar spine was developed and validated against experimental data and tested in the primary modes of spinal motion in the intact condition. Once a ligament was removed, stress increased in the remaining spinal ligaments and changes occurred in vertebral strain energy, but disc pressure remained similar. All major biomechanical changes occurred at the same spinal level as the transected ligament, with minor changes at adjacent levels. This work demonstrates that iatrogenic damage to spinal ligaments disturbs the load sharing within the spinal ligament network and may induce significant clinically relevant changes in the spinal motion segment.     

Development of a two-dimension manifold to represent high dimension mathematical models of the intracellular Mammalian circadian clock

A new focus for mathematical models of the circadian pacemaker involves the encapsulation within the models of detailed biological processes responsible for generating those circadian rhythms. Representing greater biological detail requires more mathematical equations, which pose a greater challenge for the analysis of such systems. Development of a method that retains the predominant dynamics while still providing biologically detailed information is advantageous. Two high-dimension mathematical models of intracellular mammalian circadian pacemakers, Leloup-Goldbeter and Forger-Peskin, with 19 and 73 differential equations, respectively, have been published. The authors projected each of these high-dimension models onto their respective manifold using proper orthogonal functions (POFs) obtained from the empirical decomposition of the model's phase space to obtain a 2-dimension model. The resulting 2-dimension model, represented by 2 differential equations, predicts most of the salient characteristics of a biological clock including approximately 24-h oscillations, entrainment to an LD cycle, phase response curves, and the amplitude recovery dynamics that emerge following amplitude suppression. The manifold representation simplifies the mathematical analysis, since only 2 variables need to be observed and analyzed to understand the behavior of the biological clock. This reduced model derived from a model based on biological variables can be used for the development and analysis of mathematical models of the coupled mammalian oscillators to understand the dynamics of the integrated circadian pacemaker.     

Structural behavior of human lumbar spinal motion segments

The objectives of this study were to obtain linearized stiffness matrices, and assess the linearity and hysteresis of the motion segments of the human lumbar spine under physiological conditions of axial preload and fluid environment. Also, the stiffness matrices were expressed in the form of an 'equivalent' structure that would give insights into the structural behavior of the spine. Mechanical properties of human cadaveric lumbar L2-3 and L4-5 spinal motion segments were measured in six degrees of freedom by recording forces when each of six principal displacements was applied. Each specimen was tested with axial compressive preloads of 0, 250 and 500 N. The displacements were four slow cycles of +/-0.5mm in anterior-posterior and lateral displacements, +/-0.35 mm axial displacement, +/-1.5 degrees lateral rotation and +/-1 degrees flexion-extension and torsional rotations. There were significant increases with magnitude of preload in the stiffness, hysteresis area (but not loss coefficient) and the linearity of the load-displacement relationship. The mean values of the diagonal and primary off-diagonal stiffness terms for intact motion segments increased significantly relative to values with no preload by an average factor of 1.71 and 2.11 with 250 and 500 N preload, respectively (all eight tests p<0.01). Half of the stiffness terms were greater at L4-5 than L2-3 at higher preloads. The linearized stiffness matrices at each preload magnitude were expressed as an equivalent structure consisting of a truss and a beam with a rigid posterior offset, whose geometrical properties varied with preload. These stiffness properties can be used in structural analyses of the lumbar spine.     

Descending neuronal projections to the lumbar segment of the spinal cord of cats

In the experiments, performed on cats by means of retrograde axonal horseradish peroxidase transport, localization of the sources of descending supra- and propriospinal projections to the area of the lumbar parts of the spinal cord, where the generator of locomotor movements is located. After local administration of horseradish peroxidase into the spinal cord, the greatest amount of the labelled neurons is observed in the ipsilateral reticular formation and in the Koelliker-Fuse nucleus. Comparison of the number of neurons, forming the descending projections from the I cervical up to the IV lumbar segments demonstrates an essential predominance of short propriospinal connections over the long, as well as over the supraspinal ones.     

Comparison of erector spinae and hamstring muscle activities and lumbar motion during standing knee flexion in subjects with and without lumbar extension rotation syndrome

The aim of this study was to compare the activity of the erector spinae (ES) and hamstring muscles and the amount and onset of lumbar motion during standing knee flexion between individuals with and without lumbar extension rotation syndrome. Sixteen subjects with lumbar extension rotation syndrome (10 males, 6 females) and 14 healthy subjects (8 males, 6 females) participated in this study. During the standing knee flexion, surface electromyography (EMG) was used to measure muscle activity, and surface EMG electrodes were attached to both the ES and hamstring (medial and lateral) muscles. A three-dimensional motion analysis system was used to measure kinematic data of the lumbar spine. An independent-t test was conducted for the statistical analysis. The group suffering from lumbar extension rotation syndrome exhibited asymmetric muscle activation of the ES and decreased hamstring activity. Additionally, the group with lumbar extension rotation syndrome showed greater and earlier lumbar extension and rotation during standing knee flexion compared to the control group. These data suggest that asymmetric ES muscle activation and a greater amount of and earlier lumbar motion in the sagittal and transverse plane during standing knee flexion may be an important factor contributing to low back pain.     

Rat models of traumatic spinal cord injury to assess motor recovery

Devastating motor, sensory, and autonomic dysfunctions render long-term personal hardships to the survivors of traumatic spinal cord injury (SCI). The suffering also extends to the survivors' families and friends, who endure emotional, physical, and financial burdens in providing for necessary surgeries, care, and rehabilitation. After the primary mechanical SCI, there is a complex secondary injury cascade that leads to the progressive death of otherwise potentially viable axons and cells and that impairs endogenous recovery processes. Investigations of possible cures and of ways to alleviate the hardships of traumatic SCI include those of interventions that attenuate or overcome the secondary injury cascade, enhance the endogenous repair mechanisms, regenerate axons, replace lost cells, and rehabilitate. These investigations have led to the creation of laboratory animal models of the different types of traumatic human SCI and components of the secondary injury cascade. However, no particular model completely addresses all aspects of traumatic SCI. In this article, we describe adult rat SCI models and the motor, and in some cases sensory and autonomic, deficits that each produces. Importantly, as researchers in this area move toward clinical trials to alleviate the hardships of traumatic SCI, there is a need for standardized small and large animal SCI models as well as quantitative behavioral and electrophysiological assessments of their outcomes so that investigators testing various interventions can directly compare their results and correlate them with the molecular, biochemical, and histological alterations.     

Markerless Roentgen Stereophotogrammetric Analysis for in vivo implant migration measurement using three dimensional surface models to represent bone

Recent studies have shown that model-based RSA using implant surface models to detect in vivo migration is as accurate as the classical marker-based RSA method. Use of bone surface models would be a further advancement of the model-based method by decreasing complications arising from marker insertion. The aim of this pilot investigation was to assess the feasibility of a "completely markerless" model-based RSA in detecting migration of an implant using bone surface models instead of bone markers. A total knee arthroplasty (TKA) was performed on a human cadaver knee, which was subsequently investigated by repeated RSA measurements performed by one observer. The cadaver knee was CT scanned prior to implantation of the TKA. Tibia-fibular surface models were created using two different commercially available software packages to investigate the effect of segmentation software on the accuracy of repeated migration measures of zero displacement by one observer. Reverse engineered surface models of the TKA tibial component were created. The analysis of the RSA images was repeated 10 times by one individual observer. For the markerless method, the greatest apparent migration observed about the three anatomical axes investigated was between -2.08 and 1.35 mm (SD ≤ 0.88) for z-axis translation, and -4.57° to 7.86° (SD ≤ 3.17) for R(y)-axis rotation, which were well beyond out of the range of what is typically considered adequate for clinically relevant RSA measurements. Use of tibia-fibular surface models of the bone instead of markers could provide practical advantages in evaluating implant migration. However, we found the accuracy and precision of the markerless approach to be lower than that of marker-based RSA, to a degree which precludes the use of this method for measuring implant migration in its present form.     

On computer simulation methods used to study models of two-component lipid bilayers

We show that the use of a computer simulation method introduced to calculate the equilibrium thermodynamic properties of a model of a two-component lipid bilayer membrane [Freire, E., & Snyder, B. (1980) Biochemistry 19, 88-94] is incorrect. This is done by comparing the method to that of Metropolis, which has been proven to generate equilibrium distribution of that model, and by showing that back-processes have been omitted in the implicit master equation of Freire and Snyder. We have illustrated this explicitly by first generating distributions according to the method of Freire and Snyder and then allowing the system to relax via the Kawasaki method, which uses the technique of Metropolis. We show that relaxation to a different distribution occurs. We also remark that the cluster distributions generated by the Freire-Snyder method are substantially different from those occurring in equilibrium distributions. Thus, conclusions about equilibrium thermodynamic properties such as specific heats and transition enthalpies or about transport properties or cluster properties at equilibrium cannot be drawn from the results obtained by using this method. Finally, we point out that the method of Freire and Snyder is appropriate to so-called aggregation models, which have been used to study irreversible growth; and we suggest biological systems that might be simulated by their method.     

Validation of protein models by a neural network approach

Background  

The development and improvement of reliable computational methods designed to evaluate the quality of protein models is relevant in the context of protein structure refinement, which has been recently identified as one of the bottlenecks limiting the quality and usefulness of protein structure prediction.     

FEBio finite element models of the human lumbar spine

Finite element analysis has proven to be a viable method for assessing many structure-function relationships in the human lumbar spine. Several validated models of the spine have been published, but they typically rely on commercial packages and are difficult to share between labs. The goal of this study is to present the development of the first open-access models of the human lumbar spine in FEBio. This modeling framework currently targets three deficient areas in the field of lumbar spine modeling: 1) open-access models, 2) accessibility for multiple meshing schemes, and 3) options to include advanced hyperelastic and biphasic constitutive models.