Firing rate trajectories of human occipitofrontalis motor units in response to triangular voluntary contraction intensity

A geometric approach to quantifying the neuromodulatory effects of persistent inward currents on individual motor unit discharge patterns

Objective.All motor commands flow through motoneurons, which entrain control of their innervated muscle fibers, forming a motor unit (MU). Owing to the high fidelity of action potentials within MUs, their discharge profiles detail the organization of ionotropic excitatory/inhibitory as well as metabotropic neuromodulatory commands to motoneurons. Neuromodulatory inputs (e.g. norepinephrine, serotonin) enhance motoneuron excitability and facilitate persistent inward currents (PICs). PICs introduce quantifiable properties in MU discharge profiles by augmenting depolarizing currents upon activation (i.e. PIC amplification) and facilitating discharge at lower levels of excitatory input than required for recruitment (i.e. PIC prolongation).Approach. Here, we introduce a novel geometric approach to estimate neuromodulatory and inhibitory contributions to MU discharge by exploiting discharge non-linearities introduced by PIC amplification during time-varying linear tasks. In specific, we quantify the deviation from linear discharge ('brace height') and the rate of change in discharge (i.e. acceleration slope, attenuation slope, angle). We further characterize these metrics on a simulated motoneuron pool with known excitatory, inhibitory, and neuromodulatory inputs and on human MUs (number of MUs; Tibialis Anterior: 1448, Medial Gastrocnemius: 2100, Soleus: 1062, First Dorsal Interosseus: 2296).Main results. In the simulated motor pool, we found brace height and attenuation slope to consistently indicate changes in neuromodulation and the pattern of inhibition (excitation-inhibition coupling), respectively, whereas the paired MU analysis (ΔF) was dependent on both neuromodulation and inhibition pattern. Furthermore, we provide estimates of these metrics in human MUs and show comparable variability in ΔFand brace height measures for MUs matched across multiple trials.Significance. Spanning both datasets, we found brace height quantification to provide an intuitive method for achieving graded estimates of neuromodulatory and inhibitory drive to individual MUs. This complements common techniques and provides an avenue for decoupling changes in the level of neuromodulatory and pattern of inhibitory motor commands.

Firing rate trajectories of human motor units during activity-dependent muscle potentiation

During activity-dependent potentiation (ADP) motor unit firing rates (MUFRs) are lower, however, the mechanism for this response is not known. During increasing torque isometric contractions at low contraction intensities, MUFR trajectories initially accelerate and saturate demonstrating a non-linear response due to the activation of persistent inward currents (PICs) at the motoneuron. The purpose was to assess whether PICs are a factor in the reduction of MUFRs during ADP. To assess this, MUFR trajectories were fit with competing functions of linear regression and a rising exponential (i.e., acceleration and saturation). Using fine-wire electrodes, discrete MU potential trains were recorded in the tibialis anterior during slowly increasing dorsiflexion contractions to 10% of maximal voluntary contraction following both voluntary (post-activation potentiation; PAP) and evoked (post-tetanic potentiation; PTP) contractions. In 8 participants, 25 MUs were recorded across both ADP conditions and compared to the control with no ADP effect. During PAP and PTP, the average MUFRs were 16.4% and 9.2% lower (both P≤ 0.001), respectively. More MUFR trajectories were better fit to the rising exponential during control (16/25) compared to PAP (4/25, P<0.001) and PTP (8/25, P=0.03). The MU samples that had a rising exponential MUFR trajectory during PAP and PTP displayed an ~11% lower initial acceleration compared to control (P<0.05). Thus, synaptic amplification and MUFR saturation due to PIC properties are attenuated during ADP regardless of the type of conditioning contraction. This response may contribute to lower MUFRs and likely occurred because synaptic input is reduced when contractile function is enhanced.