Orientation of tendons in vivo with active and passive knee muscles

Tendon orientations in knee models are often taken from cadaver studies. The aim of this study was to investigate the effect of muscle activation on tendon orientation in vivo. Magnetic resonance imaging (MRI) images of the knee were made during relaxation and isometric knee extensions and flexions with 0 degrees , 15 degrees and 30 degrees of knee joint flexion. For six tendons, the orientation angles in sagittal and frontal plane were calculated. In the sagittal plane, muscle activation pulled the patellar tendon to a more vertical orientation and the semitendinosus and sartorius tendons to a more posterior orientation. In the frontal plane, the semitendinosus had a less lateral orientation, the biceps femoris a more medial orientation and the patellar tendon less medial orientation in loaded compared to unloaded conditions. The knee joint angle also influenced the tendon orientations. In the sagittal plane, the patellar tendon had a more anterior orientation near full extension and the biceps femoris had an anterior orientation with 0 degrees and 15 degrees flexions and neutral with 30 degrees flexions. Within 0 degrees to 30 degrees of flexion, the biceps femoris cannot produce a posterior shear force and the anterior angle of the patellar tendon is always larger than the hamstring tendons. Therefore, co-contraction of the hamstring and quadriceps is unlikely to reduce anterior shear forces in knee angles up to 30 degrees . Finally, inter-individual variation in tendon angles was large. This suggests that the amount of shear force produced and the potential to counteract shear forces by co-contraction is subject-specific.     

Dynamic measurement of pennation angle of gastrocnemius muscles during contractions based on ultrasound imaging

ABSTRACT: BACKGROUND: Muscle fascicle pennation angle (PA) is an important parameter related to musculoskeletal functions, and ultrasound imaging has been widely used for measuring PA, but manually and frame by frame in most cases. We have earlier reported an automatic method to estimate aponeurosis orientation based on Gabor transform and Revoting Hough Transform (RVHT). METHODS: In this paper, we proposed a method to estimate the overall orientation of muscle fascicles in a region of interest, in order to complete computing the orientation of the other side of the pennation angle, but the side found by RVHT. The measurements for orientations of both fascicles and aponeurosis were conducted in each frame of ultrasound images, and then the dynamic change of pennation angle during muscle contraction was obtained automatically. The method for fascicle orientation estimation was evaluated using synthetic images with different noise levels and later on 500 ultrasound images of human gastrocnemius muscles during isometric plantarflexion. RESULTS: The muscle fascicle orientations were also estimated manually by two operators. From the results it's found that the proposed automatic method demonstrated a comparable performance to the manual method. CONCLUSIONS: With the proposed methods, ultrasound measurement for muscle pennation angles can be more widely used for functional assessment of muscles.     

Evaluation of regression-based 3-D shoulder rhythms

The movements of the humerus, the clavicle, and the scapula are not completely independent. The coupled pattern of movement of these bones is called the shoulder rhythm. To date, multiple studies have focused on providing regression-based 3-D shoulder rhythms, in which the orientations of the clavicle and the scapula are estimated by the orientation of the humerus. In this study, six existing regression-based shoulder rhythms were evaluated by an independent dataset in terms of their predictability. The datasets include the measured orientations of the humerus, the clavicle, and the scapula of 14 participants over 118 different upper arm postures. The predicted orientations of the clavicle and the scapula were derived from applying those regression-based shoulder rhythms to the humerus orientation. The results indicated that none of those regression-based shoulder rhythms provides consistently more accurate results than the others. For all the joint angles and all the shoulder rhythms, the RMSE are all greater than 5°. Among those shoulder rhythms, the scapula lateral/medial rotation has the strongest correlation between the predicted and the measured angles, while the other thoracoclavicular and thoracoscapular bone orientation angles only showed a weak to moderate correlation. Since the regression-based shoulder rhythm has been adopted for shoulder biomechanical models to estimate shoulder muscle activities and structure loads, there needs to be further investigation on how the predicted error from the shoulder rhythm affects the output of the biomechanical model.     

Estimation of instantaneous moment arms of lower-leg muscles

Muscle moment arms at the human knee and ankle were estimated from muscle length changes measured as a function of joint flexion angle in cadaver specimens. Nearly all lower-leg muscles were studied: extensor digitorum longus, extensor hallucis longus, flexor digitorum longus, flexor hallucis longus, gastrocnemius lateralis, gastrocnemius medialis, peroneus brevis, peroneus longus, peroneus tertius, plantaris, soleus, tibialis anterior, and tibialis posterior. Noise in measured muscle length was filtered by means of quintic splines. Moment arms of the mm. gastrocnemii appear to be much more dependent on joint flexion angles than was generally assumed by other investigators. Some consequences for earlier analyses are mentioned.     

Reproducibility of gastrocnemius medialis muscle architecture during treadmill running

The purpose of this study was to assess the reproducibility of fascicle length (FL) and pennation angle (PA) of gastrocnemius medialis (GM) muscle during running in vivo. Twelve male recreational long distance runners (mean ± SD; age: 24 ± 3 years, mass: 76 ± 7 kg) ran on a treadmill at a speed of 3.0 m/s, wearing their own running shoes, for two different 10 min sessions that were at least 2 days apart. For each test day 10 acceptable trials were recorded. Ankle and knee joint angle data were recorded by a Vicon 624 system with three cameras operating at 120 Hz. B-mode ultrasonography was used to examine fascicle length and pennation angle of gastrocnemius medialis muscle. The ultrasound probe was firmly secured on the muscle belly using a lightweight foam fixation. The results indicated that fascicle length and pennation angle demonstrated high reproducibility values during treadmill running both for within and between test days. The root mean square scores between the repeated waveforms of pennation angle and fascicle length were small (∼2° and ∼3.5 mm, respectively). However, ∼14 trials for pennation angle and ∼9 trials for fascicle length may be required in order to record accurate data from muscle architecture parameters. In conclusion, ultrasound measurements may be highly reproducible during dynamic movements such as treadmill running, provided that a proper fixation is used in order to assure the constant location and orientation of the ultrasound probe throughout the movement.     

A predictive model of moment-angle characteristics in human skeletal muscle: application and validation in muscles across the ankle joint

In the present work, a generic model for the prediction of moment-angle characteristics in individual human skeletal muscles is presented. The model's prediction is based on the equation M = V x Lo(-1)sigma c cos phi x d, where M, V, and Lo are the moment-generating potential of the muscle, the muscle volume and the optimal muscle fibre length, respectively, and sigma, phi and d are the stress-generating potential of the muscle fibres, their pennation angle and the tendon moment arm length, respectively, at any given joint angle. The input parameters V, Lo, sigma, phi and d can be measured or derived mechanistically. This eliminates the common problem of the necessity to estimate one or more of the input parameters in the model by fitting its outcome to experimental results often inappropriate for the function modelled. The model's output was validated by comparisons with the moment-angle characteristics of the gastrocnemius (GS) and tibialis anterior (TA) muscles in six men, determined experimentally using voluntary contractions at several combinations of ankle and knee joint angles for the GS muscle and electrical stimulation for the TA muscle. Although the model predicted realistically the pattern of moment-angle relationship in both muscles, it consistently overestimated the GS muscle M and consistently underestimated the TA muscle M, with the difference gradually increasing from dorsiflexion to plantarflexion in both cases. The average difference between predicted and measured M was 14% for the GS muscle and 10% for the TA muscle. Approximating the muscle fibres as a single sarcomere in both muscles and failing to achieve complete TA muscle activation by electrical stimulation may largely explain the differences between theory and experiment.     

A noninvasive technique for in vivo measurement of joint torques of biarticular muscles.

The objective of this work was to develop a noninvasive method to measure the joint torques produced by biarticular muscles at two joints simultaneously. During intramuscular stimulation of the cat medial gastrocnemius (MG) muscle, torques at the ankle and knee joints were calculated from forces measured in two dimensions at the end point of the cat paw under isometric conditions. The method was verified by the known anatomical properties of cat MG muscle and the tibialis anterior (TA) muscle. The MG muscle was shown to produce a significant flexion torque at the knee, besides an extension torque at the ankle. This was in agreement with its anatomical arrangement. The TA muscle produced primarily an ankle flexion torque. The small knee torque, due to measurement errors, yielded an estimate of measurement accuracy of 3.0 +/- 2.1% (n = 52). The coupling ratio of the MG muscle, defined as T(ankle)/T(knee), varied significantly with both knee and ankle angles. The profile of MG mechanical coupling agreed qualitatively with changes in limb configuration. The method can be used to measure recruitment properties of electrically stimulated biarticular muscles, and may potentially be used to study the biomechanics of biarticular coupling.     

Effect of ankle joint position and electrode placement on the estimation of the antagonistic moment during maximal plantarflexion.

During maximal efforts, antagonistic activity can significantly influence the joint moment. During maximal voluntary "isometric" contractions, certain joint rotation can not be avoided. This can influence the estimation of the antagonistic moment from the EMG activity. Our study aimed to quantify the influence on the calculated agonistic moment produced during maximal voluntary isometric plantarflexions (a) when estimating antagonistic moments at different ankle angles and (b) when placing the EMG electrodes at different portions over the m. tibialis anterior. Ten subjects performed maximal voluntary isometric plantarflexions at 90 degrees ankle angle. In order to estimate the antagonistic moment, submaximal isometric dorsiflexions were performed at various ankle angles. Moment and EMG signals from mm. triceps surae and tibialis anterior were measured. The RMS differences between plantarflexors moment calculated considering the antagonistic cocontraction estimated at the same ankle angle at which the maximal plantarflexion moment was achieved and at different ankle angles ranged from 0.10 to 2.94 Nm. The location of the electrodes led to greater RMS differences (2.35-5.18 Nm). In conclusion, an angle 10 degrees greater than the initial plantarflexion angle is enough to minimize the effect of the change in length of the m. tibialis anterior during the plantarflexion on the estimation of the plantarflexors moment. The localisation of the electrodes over the m. tibialis anterior can influence the estimation of its cocontraction during maximal plantarflexion efforts.     

A planar model of the knee joint to characterize the knee extensor mechanism

A simple planar static model of the knee joint was developed to calculate effective moment arms for the quadriceps muscle. A pathway for the instantaneous center of rotation was chosen that gives realistic orientations of the femur relative to the tibia. Using the model, nonlinear force and moment equilibrium equations were solved at one degree increments for knee flexion angles from 0 (full extension) to 90 degrees, yielding patellar orientation, patellofemoral contact force and patellar ligament force and direction with respect to both the tibial insertion point and the tibiofemoral contact point. The computer-derived results from this two-dimensional model agree with results from more complex models developed previously from experimentally obtained data. Due to our model's simplicity, however, the operation of the patellar mechanism as a lever as well as a spacer is clearly illustrated. Specifically, the thickness of the patella was found to increase the effective moment arm significantly only at flexions below 35 degrees even though the actual moment arm exhibited an increase throughout the flexion range. Lengthening either the patella or the patellar ligament altered the force transmitted from the quadriceps to the patellar ligament, significantly increasing the effective moment arm at flexions greater than 25 degrees. We conclude that the levering action of the patella is an essential mechanism of knee joint operation at moderate to high flexion angles.     

Quantitative diffusion tensor MRI-based fiber tracking of human skeletal muscle.

Diffusion-tensor magnetic resonance imaging (DT-MRI) offers great potential for understanding structure-function relationships in human skeletal muscles. The purposes of this study were to demonstrate the feasibility of using in vivo human DT-MRI fiber tracking data for making pennation angle measurements and to test the hypothesis that heterogeneity in the orientation of the tibialis anterior (TA) muscle's aponeurosis would lead to heterogeneity in pennation angle. Eight healthy subjects (5 male) were studied. T(1)-weighted anatomical MRI and DT-MRI data were acquired of the TA muscle. Fibers were tracked from the TA's aponeurosis by following the principal eigenvector. The orientations of the aponeurosis and muscle fiber tracts in the laboratory frame of reference and the orientation of the fiber tracts with respect to the aponeurosis [i.e., the pennation angle (theta)] were determined. The muscle fiber orientations, when expressed relative to the laboratory frame of reference, did not change as functions of superior-to-inferior position. The sagittal and coronal orientations of the aponeurosis did not change in practically significant manners either, but the aponeurosis' axial orientation changed by approximately 40 degrees . As a result, the mean value for theta decreased from 16.3 (SD 6.9) to 11.4 degrees (SD 5.0) along the muscle's superior-to-inferior direction. The mean value of theta was greater in the deep than in the superficial compartment. We conclude that pennation angle measurements of human muscle made using DT-MRI muscle fiber tracking are feasible and reveal that in the foot-head direction, there is heterogeneity in the pennation properties of the human TA muscle.     

In-vivo determination of 3D muscle architecture of human muscle using free hand ultrasound

Muscle architecture is an important parameter affecting the muscle function. Most of the previous studies on in-vivo muscle architecture have used in 2D ultrasound. The importance of the third dimension has not been much explored due to lack of appropriate methods. DT-MRI has been used to study muscle architecture in 3D, however, due to long scan times of about 15 min DT-MRI has not been suitable to study active muscle contractions. The purpose of this study was to develop and validate methods to determine in-vivo muscle fascicle orientations in 3D using ultrasound. We have used 2D ultrasound and a 3D position tracker system to find the 3D fascicle orientation in 3D space. 2D orientations were obtained by using automated methods developed in our previous studies and we have extended these in the current study to obtain the 3D muscle fascicle orientation in 3D space. The methods were validated using the physical phantom and we found that the mean error in the measurement was less than 0.5° in each of the three co-ordinate planes. These methods can be achieved with short scan times (less than 2 min for the gastrocnemii) and will thus enable future studies to quantify 3D muscle architecture during sub-maximal voluntary contractions.     

Arabidopsis kinetochore fiber-associated MAP65-4 cross-links microtubules and promotes microtubule bundle elongation

The acentrosomal plant mitotic spindle is uniquely structured in that it lacks opposing centrosomes at its poles and is equipped with a connective preprophase band that regulates the spatial framework for spindle orientation and mobility. These features are supported by specialized microtubule-associated proteins and motors. Here, we show that Arabidopsis thaliana MAP65-4, a non-motor microtubule associated protein (MAP) that belongs to the evolutionarily conserved MAP65 family, specifically associates with the forming mitotic spindle during prophase and with the kinetochore fibers from prometaphase to the end of anaphase. In vitro, MAP65-4 induces microtubule (MT) bundling through the formation of cross-bridges between adjacent MTs both in polar and antipolar orientations. The association of MAP65-4 with an MT bundle is concomitant with its elongation. Furthermore, MAP65-4 modulates the MT dynamic instability parameters of individual MTs within a bundle, mainly by decreasing the frequency of catastrophes and increasing the frequency of rescue events, and thereby supports the progressive lengthening of MT bundles over time. These properties are in line with its role of initiating kinetochore fibers during prospindle formation.     

Age does not influence muscle fiber length adaptation to increased excursion.

Muscle fiber length adaptation to static stretch or shortening depends on age, with sarcomere addition in young muscle being dependent on mobility. Series sarcomere number can also increase in young animals in response to increased muscle excursion, but it is not clear whether adult muscles respond similarly. The ankle flexor retinaculum was transected in neonatal and adult rats to increase tibialis anterior muscle excursion. Sarcomere number in tibialis anterior was determined after 8 wk of adaptation. Muscle moment arm and excursion were increased 30% (P < 0.01) in both age groups. Muscle cross-sectional area was reduced by 12% (P < 0.01) in response to the increased mechanical advantage, and this reduction was unaffected by age. Fiber length change was also unaffected by age, with both groups showing a trend (P < 0.10) for slightly (6%) increased fiber length. Retinaculum transection results in shorter muscle length in all joint configurations, so this trend opposes the fiber length decrease predicted by an adaptation to muscle length and indicates that fiber length is influenced by dynamic mechanical signals in addition to static length.     

Automated tracking of muscle fascicle orientation in B-mode ultrasound images

B-mode ultrasound can be used to non-invasively image muscle fascicles during both static and dynamic contractions. Digitizing these muscle fascicles can be a timely and subjective process, and usually studies have used the images to determine the linear fascicle lengths. However, fascicle orientations can vary along each fascicle (curvature) and between fascicles. The purpose of this study was to develop and test two methods for automatically tracking fascicle orientation. Images were initially filtered using a multiscale vessel enhancement (a technique used to enhance tube-like structures), and then fascicle orientations quantified using either the Radon transform or wavelet analysis. Tests on synthetic images showed that these methods could identify fascicular orientation with errors of less than 0.06°. Manual digitization of muscle fascicles during a dynamic contraction resulted in a standard deviation of angle estimates of 1.41° across ten researchers. The Radon transform predicted fascicle orientations that were not significantly different from the manually digitized values, whilst the wavelet analysis resulted in angles that were 1.35° less, and reasons for these differences are discussed. The Radon transform can be used to identify the dominant fascicular orientation within an image, and thus used to estimate muscle fascicle lengths. The wavelet analysis additionally provides information on the local fascicle orientations and can be used to quantify fascicle curvatures and regional differences with fascicle orientation across an image.     

A maximum entropy analysis of protein orientations using fluorescence polarization data from multiple probes

Techniques have recently become available to label protein subunits with fluorescent probes at predetermined orientation relative to the protein coordinates. The known local orientation enables quantitative interpretation of fluorescence polarization experiments in terms of orientation and motions of the protein within a larger macromolecular assembly. Combining data obtained from probes placed at several distinct orientations relative to the protein structure reveals functionally relevant information about the axial and azimuthal orientation of the labeled protein segment relative to its surroundings. Here we present an analytical method to determine the protein orientational distribution from such data. The method produces the broadest distribution compatible with the data by maximizing its informational entropy. The key advantages of this approach are that no a priori assumptions are required about the shape of the distribution and that a unique, exact fit to the data is obtained. The relative orientations of the probes used for the experiments have great influence on information content of the maximum entropy distribution. Therefore, the choice of probe orientations is crucial. In particular, the probes must access independent aspects of the protein orientation, and two-fold rotational symmetries must be avoided. For a set of probes, a "figure of merit" is proposed, based on the independence among the probe orientations. With simulated fluorescence polarization data, we tested the capacity of maximum entropy analysis to recover specific protein orientational distributions and found that it is capable of recovering orientational distributions with one and two peaks. The similarity between the maximum entropy distribution and the test distribution improves gradually as the number of independent probe orientations increases. As a practical example, ME distributions were determined with experimental data from muscle fibers labeled with bifunctional rhodamine at known orientations with respect to the myosin regulatory light chain (RLC). These distributions show a complex relationship between the axial orientation of the RLC relative to the fiber axis and the azimuthal orientation of the RLC about its own axis. Maximum entropy analysis reveals limitations in available experimental data and supports the design of further probe angles to resolve details of the orientational distribution.     

Changes in muscle architecture and performance during a competitive season in female softball players

ABSTRACT: Nimphius, S, McGuigan, MR, and Newton, RU. Changes in muscle architecture and performance during a competitive season in female softball players. J Strength Cond Res 26(10): 2655-2666, 2012-The purpose of this research was (a) to examine the performance changes that occur in elite female softball players during 20 weeks of softball training (that included 14 weeks of periodized resistance training [RT]) and (b) to examine the relationship between percent change (%change) in muscle architecture variables and %change in strength, speed, and change of direction performance. Ten female softball players (age = 18.1 ± 1.6 years, height = 166.5 ± 8.9 cm, weight = 72.4 ± 10.8 kg) from a state Australian Institute of Sport softball team were tested for maximal lower-body strength using a 3 repetition maximum for a predicted 1 repetition maximum (1RM) and peak force, peak velocity (PV), and peak power (PP) were measured during jump squats (JS) unloaded and loaded. In addition, the first base (1B) and the second base (2B) sprint performance, change of direction (505) on dominant (D) and nondominant (ND) sides, aerobic capacity, and muscle architecture characteristics of vastus lateralis (VL) including muscle thickness (MT), fascicle length (FL), and pennation angle (θp) were examined. The testing sessions occurred pre, mid, and post training (total 20 week pre- and in-season training period). Changes over time were analyzed by repeated-measures analysis of variance. The relationship between %change in muscle architecture variables and strength, speed, and change of direction variables from pre to post were assessed by Pearson product-moment correlation coefficient. Significant improvements in PV and PP occurred at all JS loads pre- to mid-testing and pre- to post-testing. Significant increases occurred pre-post in absolute 1RM, relative 1RM, 505 ND, and 2B sprint. The strongest relationships were found between %change in VL MT and 1B sprint (r = -0.80, p = 0.06), %change in VL FL and 2B sprint (r = -0.835, p = 0.02), and %change in relative 1RM and 505 D (r = -0.70, p = 0.04). In conclusion, gains in strength, power, and performance can occur during the season in elite softball players who are also engaged in a periodized RT program. Furthermore, changes in performance measures are associated with changes in muscle architecture.     

Changes in fascicle lengths and pennation angles do not contribute to residual force enhancement/depression in voluntary contractions

Force enhancement following muscle stretching and force depression following muscle shortening are well-accepted properties of skeletal muscle contraction. However, the factors contributing to force enhancement/depression remain a matter of debate. In addition to factors on the fiber or sarcomere level, fiber length and angle of pennation affect the force during voluntary isometric contractions in whole muscles. Therefore, we hypothesized that differences in fiber lengths and angles of pennation between force-enhanced/depressed and reference states may contribute to force enhancement/depression during voluntary contractions. The purpose of this study was to test this hypothesis. Twelve subjects participated in this study, and force enhancement/depression was measured in human tibialis anterior. Fiber lengths and angles of pennation were quantified using ultrasound imaging. Neither fiber lengths nor angles of pennation were found to differ between the isometric reference contractions and any of the force-enhanced or force-depressed conditions. Therefore, we rejected our hypothesis and concluded that differences in fiber lengths or angles of pennation do not contribute to the observed force enhancement/depression in human tibialis anterior, and speculate that this result is likely true for other muscles too.     

Reliability and robustness of muscle architecture measurements obtained using diffusion tensor imaging with anatomically constrained tractography

For detailed analyses of muscle adaptation mechanisms during growth, ageing or disease, reliable measurements of muscle architecture are required. Diffusion tensor imaging (DTI) and DTI tractography have been used to reconstruct the architecture of human muscles in vivo. However, muscle architecture measurements reconstructed with conventional DTI techniques are often anatomically implausible because the reconstructed fascicles do not terminate on aponeuroses, as real muscle fascicles are known to do. In this study, we tested the reliability of an anatomically constrained DTI-based method for measuring three-dimensional muscle architecture. Anatomical magnetic resonance images and diffusion tensor images were obtained from the left legs of eight healthy participants on two occasions one week apart. Muscle volumes, fascicle lengths, pennation angles and fascicle curvatures were measured in the medial and lateral gastrocnemius, soleus and the tibialis anterior muscles. Averaged across muscles, the intraclass correlation coefficient was 0.99 for muscle volume, 0.81 for fascicle length, 0.73 for pennation angle and 0.76 for fascicle curvature. Measurements of muscle architecture obtained using conventional DTI tractography were highly sensitive to variations in the stopping criteria for DTI tractography. The application of anatomical constraints reduced this sensitivity significantly. This study demonstrates that anatomically constrained DTI tractography can provide reliable and robust three-dimensional measurements of whole-muscle architecture. The algorithms used to constrain tractography have been made publicly available.     

Predictability of in vivo changes in pennation angle of human tibialis anterior muscle from rest to maximum isometric dorsiflexion

The aim of this study was to assess the predictability of in vivo, ultrasound-based changes in human tibialis anterior (TA) pennation angle from rest to maximum isometric dorsiflexion (MVC) using a planimetric model assuming constant thickness between aponeuroses and straight muscle fibres. Sagittal sonographs of TA were taken in six males at ankle angles of -15 degrees (dorsiflexion direction), 0 degrees (neutral position), + 15 (plantarflexion direction) and + 30 degrees both at rest and during dorsiflexor MVC trials performed on an isokinetic dynamometer. At all four ankle angles scans were taken from the TA proximal, central and distal regions. TA architecture did not differ (P > 0.05) neither between its two unipennate parts nor along the scanned regions over its length at a given ankle angle and state of contraction. Comparing MVC with rest at any given ankle angle, pennation angle was larger (62-71%, P < 0.01), fibre length smaller (37-40%, P < 0.01) and muscle thickness unchanged (P > 0.05). The model used estimated accurately (P > 0.05) changes in TA pennation angle occurring in the transition from rest to MVC and therefore its use is encouraged for estimating the isometric TA ankle moment and force generating capacity using musculoskeletal modelling.