Computer simulation of the motoneuron pool-muscle complex. I. Input system and motoneuron pool.
The aim of the present study was to simulate the input system and the motoneuron (MN) pool of the MN pool-muscle complex (MNPMC). Input fibers, which can originate from command centers in the central nervous system or from sensory organs, activate the MN pool. They generate sequences of action potentials, the frequency of which is proportional to a time-dependent activation factor (which is an input to the model). Different connection patterns between the input fibers and motor units (MUs) are allowed. For simplicity and since no precise experimental data are available, 70 input fibers and 4 boutons per fiber and MN are simulated (this corresponds approximately to the monosynaptic group-Ia input of the cat medial gastrocnemius muscle). Each bouton generates the same conductance change in the postsynaptic membrane. The MNs are modeled with a single compartment and a homogeneous membrane. According to experimental data, the membrane leakage conductance and capacitance are MU dependent. Since the precise relation is unknown: (a) the computed relation between MU contraction force and the MN leakage conductance was taken from a steady-state MNPMC model, and (b) the capacitance was assumed to be proportional to the leakage conductance. The MN membrane includes time- and voltage-dependent ionic channels (fast and slow K(+) and low- and high-threshold Ca(2+) channels). The density and time constant of the slow K(+) channels and the density of the Ca(2+) channels were fitted to approximate afterhyperpolarization characteristics and frequency-injected current relations of type-identified cat MNs. If the membrane reaches a voltage threshold the MNs generate action potentials, which were simulated by voltage pulses. The activation of the MN pool of the human first dorsal interosseus muscle was simulated with injected and synaptic currents in order to illustrate the size principle, synaptic noise, and other features of muscle activation. It is concluded that the present model reproduces the main properties of the input-output relations of different MN types within a muscle. Together with the simulation of the muscle force and the surface EMG, which will be published in subsequent papers, it will be a powerful tool for reproducing experiments on the motor system and investigating functional mechanisms of motor control.
Nussbaumer
R M
RM
Department of Physiology, University of Fribourg, Rue du Musée 5, 1700 Fribourg, Switzerland.
Ruegg
D G
DG
Studer
L M
LM
Gabriel
J-P
JP
D000344:Afferent Pathways / Q000502:physiology; D000465:Algorithms; D002490:Central Nervous System / Q000502:physiology; D003198:Computer Simulation; D019706:Excitatory Postsynaptic Potentials; D006801:Humans; D008959:Models, Neurological; D009046:Motor Neurons / Q000502:physiology*; D009119:Muscle Contraction / Q000502:physiology*; D018482:Muscle, Skeletal / Q000294:innervation*; D009412:Nerve Fibers / Q000502:physiology*