Tremor and other oscillations in neuromuscular systems

A model has been analyzed which is based on recent experimental evidence concerning the properties of muscles and the sensory feedback pathways from muscles. Damped oscillations can arise in the absence of sensory feedback due to the interaction of a muscle with inertial loads. These mechanical oscillations can have a wide range of frequencies depending on the inertial and elastic loads that are attached to the muscle. Small amounts of sensory feedback will tend to reduce deviations from a steady muscle length, but larger amounts of feedback can produce oscillations. The frequency of these reflex oscillations is determined by the properties of the muscle and feedback pathway, and is rather independent of load. If the strength of the sensory feedback is sufficient, either the mechanical oscillations or the reflex oscillations or both can grow, rather than decay, with time. The growth of these oscillations is limited by saturation non-linearities in the muscle receptors and the muscle itself, so that the oscillations approach a steady amplitude and frequency. Using typical properties of muscles and spinal reflex pathways, the frequency of reflex oscillations will be within the range 8–12 Hz found for physiological tremor. With the longer latency found for supraspinal reflexes, oscillations will occur in the range 4–6 Hz which is characteristic of Parkinson's and cerebellar diseases. The role of longer latency reflexes in the generation of these tremors is discussed.     

Effects of elastic loads on the contractions of cat muscles

The nerves to plantaris and soleus muscles in the cat were stimulated with maximal single shocks and with random stimulus trains which produced partially fused contractions. In order to obtain information on the mechanism of muscular contraction, the effects of allowing the muscles to shorten against various elastic loads were studied in the time domain and in the frequency domain. When springs of increasing stiffness were placed in series with the muscle, the twitch tension increased greatly. The gain of the frequency response curve was also much greater with stiffer springs. The shape of the frequency response curve for plantaris muscle could usually be described by that expected for a second-order system with two real time constants or rate constants. The rate constants changed in qualitatively similar ways in response to increased stiffness of an elastic load, increased muscle length and increased mean rate of nerve stimulation. These results are in agreement with the hypothesis that the linear responses of muscles working against elastic loads are determined by the values of two rate constants. Thus, of the many processes associated with contraction, only two are rate-limiting: one associated with the viscoelastic properties of muscle and the second associated with the reuptake of Ca into the sarcoplasmic reticulum. Non-linear aspects of muscular contraction are also discussed. These are more prominent in soleus muscle than in plantaris muscle.Graduate student of the Medical Research Council of Canada.Formerly a Post-doctoral Fellow of the Muscular Dystrophy Association of Canada.     

Mathematical model that predicts lower leg motion in response to electrical stimulation

Electrical stimulation of skeletal muscles of patients with upper motor neuron lesions can be used to restore functional movements such as standing or walking. Mathematical muscle models can assist in designing stimulation patterns that will enable patients to perform particular tasks more efficiently. In this study we extend our previous model to allow us to predict changes in knee joint angle in response to electrical stimulation of the human quadriceps femoris muscle. The model was tested both with and without inertial loads placed around the ankle joints of healthy subjects. Results showed that the model predicted the knee extensions with a RMS angle error that was generally 相似文献   

Seated whole body vibrations with high-magnitude accelerations--relative roles of inertia and muscle forces

Reliable computation of spinal loads and trunk stability under whole body vibrations with high acceleration contents requires accurate estimation of trunk muscle activities that are often overlooked in existing biodynamic models. A finite element model of the spine that accounts for nonlinear load- and direction-dependent properties of lumbar segments, complex geometry and musculature of the spine, and dynamic characteristics of the trunk was used in our iterative kinematics-driven approach to predict trunk biodynamics in measured vehicle's seat vibrations with shock contents of about 4g (g: gravity acceleration of 9.8m/s(2)) at frequencies of about 4 and 20Hz. Muscle forces, spinal loads and trunk stability were evaluated for two lumbar postures (erect and flexed) with and without coactivity in abdominal muscles. Estimated peak spinal loads were substantially larger under 4Hz excitation frequency as compared to 20Hz with the contribution of muscle forces exceeding that of inertial forces. Flattening of the lumbar lordosis from an erect to a flexed posture and antagonistic coactivity in abdominal muscles, both noticeably increased forces on the spine while substantially improving trunk stability. Our predictions clearly demonstrated the significant role of muscles in trunk biodynamics and associated risk of back injuries. High-magnitude accelerations in seat vibration, especially at near-resonant frequency, expose the vertebral column to large forces and high risk of injury by significantly increasing muscle activities in response to equilibrium and stability demands.     

Inertial artifact in viscoelastic measurements of striated muscle: Modeling and experimental results

Viscoelastic properties of striated muscle are often measured using length perturbation analysis and quantified as a complex modulus, whose elastic and viscous components reflect the energy-storage and energy-absorbing properties of the tissue, respectively. The energy stored as inertia is commonly ignored due to the small size of samples examined, typically <1 mm. Considering recent advances in tissue engineering to generate muscle tissues of larger sizes, we questioned whether ignoring the inertial artifact was still reasonable in these samples. To answer this question, we derived and solved the one-dimensional wave equation that describes the propagation of strain along the length of a sample. The inertial artifact was predicted to contaminate the elastic modulus with (2πf)²L₀²ρ/6, where f is perturbation frequency, L₀ is muscle length, and ρ is muscle density. We then measured viscoelastic properties up to 500 Hz in mouse skeletal muscle fibers at long (4.8 mm) and short (<1 mm) lengths and up to 100 Hz in rat cardiac slices at long (10–12 mm) and short (<2 mm) lengths. We found the elastic modulus of long preparations was elevated as frequency increased and was about half the magnitude of that predicted by the model. While the prediction tended to overestimate the measured inertial artifact, these results provided some validity to the model. We used the predicted artifact as an overly conservative estimate of error that might arise in a mechanics assay of mammalian striated muscle, whose nominal resting stiffness is on the order 100 kN m^?2. We found that muscle lengths of <1 mm resulted in negligible inertial artifact (<0.5% error) for perturbation frequencies under 250 Hz. Muscle samples longer than 5 mm, on the other hand, would result in >5% error at frequencies of 200 Hz and higher.     

Bone conduction and seismic sensitivity in golden moles (Chrysochloridae)

Some genera of golden moles are known to possess enormously hypertrophied auditory ossicles. These structures have been implicated as potentially mediating a form of inertial bone conduction, used by the golden mole to detect seismic vibrations. A simple model of ossicular inertial bone conduction, based on an existing model of the human middle ear from the literature, was used in an attempt to examine vibrational sensitivity in these animals. Those golden moles with hypertrophied ossicles are predicted to possess a sensitive inertial bone conduction response at frequencies below a few hundred hertz, whereas species lacking these middle ear adaptations are predicted to have a far less sensitive response in this ecologically important frequency range. An alternative mode of inertial bone conduction in golden moles, potentially conferring sensitivity to vertically-polarized seismic vibrations such as Rayleigh waves, is proposed. Certain behaviours of golden moles described in the literature are interpreted as augmenting seismic sensitivity.     

Directional specificity and encoding of muscle forces and loads by stick insect tibial campaniform sensilla,including receptors with round cuticular caps

In many systems, loads are detected as the resistance to muscle contractions. We studied responses to loads and muscle forces in stick insect tibial campaniform sensilla, including a subgroup of receptors (Group 6B) with unusual round cuticular caps in oval-shaped collars. Loads were applied in different directions and muscle contractions were emulated by applying forces to the tibial flexor muscle tendon (apodeme). All sensilla 1) were maximally sensitive to loads applied in the plane of joint movement and 2) encoded muscle forces but did not discharge to unresisted movements. Identification of 6B sensilla by stimulation of cuticular caps demonstrated that receptor responses were correlated with their morphology. Sensilla with small cuticular collars produced small extracellular potentials, had low thresholds and strong tonic sensitivities that saturated at moderate levels. These receptors could effectively signal sustained loads. The largest spikes, derived from sensilla with large cuticular collars, had strong dynamic sensitivities and signaled a wide range of muscle forces and loads. Tibial sensilla are apparently tuned to produce no responses to inertial forces, as occur in the swing phase of walking. This conclusion is supported by tests in which animals 'stepped' on a compliant surface and sensory discharges only occurred in stance.     

Body segment inertial parameters and low back load in individuals with central adiposity

There is a paucity of information regarding the impact of central adiposity on the inertial characteristics of body segments. Deriving low back loads during lifting requires accurate estimate of inertial parameters. The purpose was to determine the body segment inertial parameters of people with central adiposity using a photogrammetric technique, and then to evaluate the impact on lumbar spine loading. Five participants with central adiposity (waist:hip ratio>0.9, waist circumference>102 cm) were compared to a normal BMI group. A 3D wireframe model of the surface topography was constructed, partitioned into 8 body segments and then body segment inertial parameters were calculated using volumetric integration assuming uniform segment densities for the segments. Central adiposity dependent increases in body segment parameters ranged from 12 to 400%, varying across segments (greatest for trunk) and parameters. The increase in mass distribution to the trunk was accompanied by an anterior and inferior shift of the centre of mass. A proximal shift in centre of mass was detected for the extremities, along with a reduction in mass distribution to the lower extremity. L5/S1 torques (392 vs 263 Nm) and compressive forces (5918 vs 3986 N) were substantially elevated in comparison to the normal BMI group, as well as in comparison to torques and forces predicted using published BSIP equations. Central adiposity resulted in substantial but non-uniform increases in inertial parameters resulting in task specific increases in torque and compressive loads arising from different inertial and physical components.     

Predicting optimal electrical stimulation for repetitive human muscle activation.

Functional electrical stimulation is the use of electrical currents to activate paralyzed muscles to produce functional movements. Muscle force output must meet or exceed the external load to maintain a posture or produce movements. A mathematical force-fatigue modeling system that predicts muscle force responses during repetitive electrical stimulation has been developed in our laboratory to help identify stimulation patterns that optimize force output for individual subjects. This study tests how well this model predicts the number of contractions that can be maintained above a required force level (successful contractions) during repetitive activation of a muscle. Healthy human quadriceps muscles were tested isometrically on 12 subjects. Data were first collected and used to parameterize the model. Next, the model was used to predict the number of successful contractions that were produced by trains with frequencies ranging from 5 to 100 Hz while the pulse durations and amplitudes were held constant. Finally, three clinically relevant stimulation frequencies were selected and tested to verify the model's predictions. Under these conditions, the model accurately predicted the number of successful contractions for clinically relevant stimulation frequencies. Furthermore, the model appears to have the potential to identify the stimulation frequency that maximizes muscle force output and minimizes fatigue for each subject.     

The mechanical response of the active human triceps brachii muscle to very rapid stretch and shortening

Very rapid, small amplitude, ramp-and-hold rotations were imposed on the braced forearms of three normal adult male subjects who were isometrically contracting their elbow extensors. By carefully accounting for inertial and viscoelastic coupling effects in the experimental system it was possible to compute the time course of the muscle-moment evoked by these mechanical perturbations. The muscle-moment responses, and their dependence on rotation amplitude and direction, as well as tonic contraction level, are described. These responses are also compared to the predictions of a simple muscle model which we have proposed previously on the basis of frequency-response tests. The results indicate that: at a given tonic contraction level, triceps may be stiffer in an isometric state than in an oscillatory steady state, and high frequency fluctuations in the myoelectric activity are very ineffective in generating corresponding muscle-force fluctuations.     

Model-generated decomposition of unfused tetani of motor units evoked by random stimulation

Unfused tetani of motor units (MUs) evoked by stimulation at variable interpulse intervals at mean frequencies of 20, 25, 33, 40 and 50 Hz were studied using ten functionally isolated fast-type MUs from the medial gastrocnemius muscle of adult Wistar rats. A previously proposed algorithm and computer program for mathematical decomposition of unfused tetani into a series of twitches, representing responses to individual pulses, were used. Analysis of the parameters of the decomposed twitches showed considerable variability in force of successive contractions. These twitches were extremely variable with up to 2-fold higher forces and longer contraction times than a single twitch evoked by one stimulus. However, when the stimulation frequency was decreased, the decomposed twitches became similar to the single twitch with respect to amplitude and contraction time. It was found that the basic contractile parameters of decomposed twitches could be predicted with high accuracy on the basis of the tetanus force level at which the next contraction begins. This analysis of the parameters of decomposed twitches demonstrated that the contractile responses of the muscle fibers to successive action potentials generated by motoneurons are highly variable and depend on the previous MU state.     

Optimal Workloop Energetics of Muscle-Actuated Systems: An Impedance Matching View

Integrative approaches to studying the coupled dynamics of skeletal muscles with their loads while under neural control have focused largely on questions pertaining to the postural and dynamical stability of animals and humans. Prior studies have focused on how the central nervous system actively modulates muscle mechanical impedance to generate and stabilize motion and posture. However, the question of whether muscle impedance properties can be neurally modulated to create favorable mechanical energetics, particularly in the context of periodic tasks, remains open. Through muscle stiffness tuning, we hypothesize that a pair of antagonist muscles acting against a common load may produce significantly more power synergistically than individually when impedance matching conditions are met between muscle and load. Since neurally modulated muscle stiffness contributes to the coupled muscle-load stiffness, we further anticipate that power-optimal oscillation frequencies will occur at frequencies greater than the natural frequency of the load. These hypotheses were evaluated computationally by applying optimal control methods to a bilinear muscle model, and also evaluated through in vitro measurements on frog Plantaris longus muscles acting individually and in pairs upon a mass-spring-damper load. We find a 7-fold increase in mechanical power when antagonist muscles act synergistically compared to individually at a frequency higher than the load natural frequency. These observed behaviors are interpreted in the context of resonance tuning and the engineering notion of impedance matching. These findings suggest that the central nervous system can adopt strategies to harness inherent muscle impedance in relation to external loads to attain favorable mechanical energetics.     

Modified Bilston nonlinear viscoelastic model for finite element head injury studies

This paper proposes a modified nonlinear viscoelastic Bilston model (Bilston et al., 2001, Biorheol., 38, pp. 335-345). for the modeling of brain tissue constitutive properties. The modified model can be readily implemented in a commercial explicit finite element (FE) code, PamCrash. Critical parameters of the model have been determined through a series of rheological tests on porcine brain tissue samples and the time-temperature superposition (TTS) principle has been used to extend the frequency to a high region. Simulations by using PamCrash are compared with the test results. Through the use of the TTS principle, the mechanical and rheological behavior at high frequencies up to 10(4) rads may be obtained. This is important because the properties of the brain tissue at high frequencies and impact rates are especially relevant to studies of traumatic head injury. The averaged dynamic modulus ranges from 130 Pa to 1500 Pa and loss modulus ranges from 35 Pa to 800 Pa in the frequency regime studied (0.01 rads to 3700 rads). The errors between theoretical predictions and averaged relaxation test results are within 20% for strains up to 20%. The FEM simulation results are in good agreement with experimental results. The proposed model will be especially useful for application to FE analysis of the head under impact loads. More realistic analysis of head injury can be carried out by incorporating the nonlinear viscoelastic constitutive law for brain tissue into a commercial FE code.     

EMG responses resulting from transient and steady-state dynamic isometric loading of the human biceps can be distinguished.

In order to determine if differences in human muscle EMG response to steady-state (STD) dynamic load and transient (TRN) load could be detected, two distinct sinusoidal loads having identical amplitude and frequency were applied to the forearm, which was held in a fixed posture. The first condition used application of a constant amplitude, constant frequency, STD sinusoid, and the second condition used application of a single sinusoid of identical amplitude and frequency as the STD load. Time-domain parameters of Half-Mass, Dispersion Moments, Asymmetry Moments, Peak Ratio (ratio of cycle peak RMS EMG to cycle mean RMS EMG) and Time-to-Peak Ratio (ratio of cycle time to reach peak RMS EMG to cycle time duration) were used to quantify responses in the time domain. A single frequency-domain parameter, Dispersion Moments, was used to quantify responses in the frequency domain. Individual ANOVA's of the time domain parameter results revealed that p = 0.689 for Dispersion Moments, p = 0.111 for Half-Mass, p = 0.007 for Asymmetry Moments, p = 0.001 for Peak Ratio, and p = 0.001 for Time-to-Peak Ratio. The time-domain MANOVA and frequency-domain ANOVA revealed that EMG responses resulting from STD and TRN applied loads were statistically different in the time domain (p = 0.001), but not in the frequency domain (p = 0.810). Thus, use of Asymmetry Moments, Peak Ratio, and Time-to-Peak Ratio parameters can distinguish EMG responses to STD and TRN loads.     

Frequency responses of cat rapidly adapting mechanoreceptive fibers

The frequency responses of 11 rapidly adapting (RA) fibers in cat were studied by representing the average firing rate as a function of sinusoidal stimulus amplitude and stimulus frequency. Specifically, rate-intensity functions at different stimulation frequencies were fitted by four-parameter (a0, a1, a2, a3), piece-wise linear functions using nonlinear regression (n = 59; R2 > 0.877). Rate-intensity functions at intermediate frequencies were found by linear interpolation. The result of this analysis is rate-amplitude-frequency functions plotted as two-dimensional surfaces. The surfaces consist of five regions separated and sufficiently defined by four space curves. At 14 different frequencies, the statistical distribution of each rate-intensity-function parameter could be approximated by a particular lognormal distribution (n = 56; R2 > 0.796). The Kolmogorov-Smirnov test fails to reject this hypothesis for each combination of frequency and parameter (56 tests; p > 0.39). Therefore, at a given frequency, the variation of the parameters can be represented by lognormal distributions with specific means and standard deviations. Responses of six RA fibers, which are different from the data-set used for modeling, were compared with the stochastic model at different frequencies. The parameters of those fibers were tested against the null hypotheses that they were sampled from the particular parameter distributions dictated by the model. The Kolmogorov-Smirnov test fails to reject all the hypotheses at the alpha = 0.05 level (44 tests). At the alpha = 0.10 level, only a few test parameters were found to be departing from the model (a0 and a1 at 5 Hz; a2 at 20 Hz; a2 and a3 at 50 Hz). The remaining test parameters could be accurately described by the model. Having confirmed the validity of the model, the logarithmic means and the logarithmic standard deviations of the lognormally distributed rate-intensity-function parameters were estimated in the frequency range of 4-200 Hz. The rate-amplitude-frequency surfaces sampled from the established stochastic model completely characterize the rate responses of RA fibers to sinusoidal stimuli and are superior to tuning curves which require selecting criterion responses. The current rate-response model is promising for future computational work, especially on population modeling.     

Optimising muscle parameters in musculoskeletal models using Monte Carlo simulation

The use of musculoskeletal simulation software has become a useful tool for modelling joint and muscle forces during human activity, including in reduced gravity because direct experimentation is difficult. Knowledge of muscle and joint loads can better inform the design of exercise protocols and exercise countermeasure equipment. In this study, the LifeModeler? (San Clemente, CA, USA) biomechanics simulation software was used to model a squat exercise. The initial model using default parameters yielded physiologically reasonable hip-joint forces but no activation was predicted in some large muscles such as rectus femoris, which have been shown to be active in 1-g performance of the activity. Parametric testing was conducted using Monte Carlo methods and combinatorial reduction to find a muscle parameter set that more closely matched physiologically observed activation patterns during the squat exercise. The rectus femoris was predicted to peak at 60.1% activation in the same test case compared to 19.2% activation using default parameters. These results indicate the critical role that muscle parameters play in joint force estimation and the need for exploration of the solution space to achieve physiologically realistic muscle activation.     

Resonant frequencies of arms and legs identify different walking patterns

The present study is aimed at investigating changes in the coordination of arm and leg movements in young healthy subjects. It was hypothesized that with changes in walking velocity there is a change in frequency and phase coupling between the arms and the legs. In addition, it was hypothesized that the preferred frequencies of the different coordination patterns can be predicted on the basis of the resonant frequencies of arms and legs with a simple pendulum model. The kinematics of arms and legs during treadmill walking in seven healthy subjects were recorded with accelerometers in the sagittal plane at a wide range of different velocities (i.e., 0.3-1. 3m/s). Power spectral analyses revealed a statistically significant change in the frequency relation between arms and legs, i.e., within the velocity range 0.3-0.7m/s arm movement frequencies were dominantly synchronized with the step frequency, whereas from 0.8m/s onwards arm frequencies were locked onto stride frequency. Significant effects of walking speed on mean relative phase between leg and arm movements were found. All limb pairs showed a significantly more stable coordination pattern from 0.8 to 1.0m/s onwards. Results from the pendulum modelling demonstrated that for most subjects at low-velocity preferred movement frequencies of the arms are predicted by the resonant frequencies of individual arms (about 0.98Hz), whereas at higher velocities these are predicted on the basis of the resonant frequencies of the individual legs (about 0.85Hz). The results support the above-mentioned hypotheses, and suggest that different patterns of coordination, as shown by changes in frequency coupling and phase relations, can exist within the human walking mode.     

A viscoelastic model of mechanically induced and spontaneous contractions of vascular smooth muscle

A new viscoelastic model was developed for the mathematical characterisation of mechanically induced and intrinsic contractional responses of the vascular smooth muscle. To this end, the elastic and viscous analogue elements were supplemented with a new active element generating stress proportional to its momentary elongation. The four-element model consisting of an active element, a parallel viscous element and both series and parallel elastic elements predicted biphasic or damped oscillatory stress relaxation and creep responses which were similar to that found experimentally earlier. Above a certain exciting frequency the model exhibited dissipative and below energy producing behaviour, as indicated by the sign change of the hysteresis loop area. In the case of sinusoidal modulation of the stress generation parameter the model showed parametric resonance, which was regarded as a simulation of intrinsic oscillation of the smooth muscle.     

Modeling the length dependence of isometric force in human quadriceps muscles

Functional electrical stimulation is used to restore movement and function of paralyzed muscles by activating skeletal muscle artificially. An accurate and predictive mathematical model can facilitate the design of stimulation patterns that produce the desired force. The present study is a first step in developing a mathematical model for non-isometric muscle contractions. The goals of this study were to: (1) identify how our isometric force model's parameters vary with changes in knee joint angle, (2) identify the best knee flexion angle to parameterize this model, and (3) validate the model by comparing experimental data to predictions in response to a wide range of stimulation frequencies and muscle lengths. Results showed that by parabolically varying one of the free parameters with knee joint angle and fixing the other parameters at the values identified at 40 degrees of knee flexion, the model could predict the force responses to a wide range of stimulation frequencies and patterns at different muscle lengths. This work showed that the current isometric force model is capable of predicting the changes in skeletal muscle force at different muscle lengths.