Optimality of position commands to horizontal eye muscles: A test of the minimum-norm rule

Three-dimensional ocular kinematics underlying binocular single vision

We have analyzed the binocular coordination of the eyes during far-to-near refixation saccades based on the evaluation of distance ratios and angular directions of the projected target images relative to the eyes' rotation centers. By defining the geometric point of binocular single vision, called Helmholtz point, we found that disparities during fixations of targets at near distances were limited in the subject's three-dimensional visual field to the vertical and forward directions. These disparities collapsed to simple vertical disparities in the projective binocular image plane. Subjects were able to perfectly fuse the vertically disparate target images with respect to the projected Helmholtz point of single binocular vision, independent of the particular location relative to the horizontal plane of regard. Target image fusion was achieved by binocular torsion combined with corrective modulations of the differential half-vergence angles of the eyes in the horizontal plane. Our findings support the notion that oculomotor control combines vergence in the horizontal plane of regard with active torsion in the frontal plane to achieve fusion of the dichoptic binocular target images.

Cerebellar-inspired algorithm for adaptive control of nonlinear dielectric elastomer-based artificial muscle

Electroactive polymer actuators are important for soft robotics, but can be difficult to control because of compliance, creep and nonlinearities. Because biological control mechanisms have evolved to deal with such problems, we investigated whether a control scheme based on the cerebellum would be useful for controlling a nonlinear dielectric elastomer actuator, a class of artificial muscle. The cerebellum was represented by the adaptive filter model, and acted in parallel with a brainstem, an approximate inverse plant model. The recurrent connections between the two allowed for direct use of sensory error to adjust motor commands. Accurate tracking of a displacement command in the actuator's nonlinear range was achieved by either semi-linear basis functions in the cerebellar model or semi-linear functions in the brainstem corresponding to recruitment in biological muscle. In addition, allowing transfer of training between cerebellum and brainstem as has been observed in the vestibulo-ocular reflex prevented the steady increase in cerebellar output otherwise required to deal with creep. The extensibility and relative simplicity of the cerebellar-based adaptive-inverse control scheme suggests that it is a plausible candidate for controlling this type of actuator. Moreover, its performance highlights important features of biological control, particularly nonlinear basis functions, recruitment and transfer of training.

Do muscle synergies reduce the dimensionality of behavior?

The muscle synergy hypothesis is an archetype of the notion of Dimensionality Reduction (DR) occurring in the central nervous system due to modular organization. Toward validating this hypothesis, it is important to understand if muscle synergies can reduce the state-space dimensionality while maintaining task control. In this paper we present a scheme for investigating this reduction utilizing the temporal muscle synergy formulation. Our approach is based on the observation that constraining the control input to a weighted combination of temporal muscle synergies also constrains the dynamic behavior of a system in a trajectory-specific manner. We compute this constrained reformulation of system dynamics and then use the method of system balancing for quantifying the DR; we term this approach as Trajectory Specific Dimensionality Analysis (TSDA). We then investigate the consequence of minimization of the dimensionality for a given task. These methods are tested in simulations on a linear (tethered mass) and a non-linear (compliant kinematic chain) system. Dimensionality of various reaching trajectories is compared when using idealized temporal synergies. We show that as a consequence of this Minimum Dimensional Control (MDC) model, smooth straight-line Cartesian trajectories with bell-shaped velocity profiles emerged as the optima for the reaching task. We also investigated the effect on dimensionality due to adding via-points to a trajectory. The results indicate that a trajectory and synergy basis specific DR of behavior results from muscle synergy control. The implications of these results for the synergy hypothesis, optimal motor control, motor development, and robotics are discussed.

Adaptive filters and internal models: multilevel description of cerebellar function

Cerebellar function is increasingly discussed in terms of engineering schemes for motor control and signal processing that involve internal models. To address the relation between the cerebellum and internal models, we adopt the chip metaphor that has been used to represent the combination of a homogeneous cerebellar cortical microcircuit with individual microzones having unique external connections. This metaphor indicates that identifying the function of a particular cerebellar chip requires knowledge of both the general microcircuit algorithm and the chip's individual connections. Here we use a popular candidate algorithm as embodied in the adaptive filter, which learns to decorrelate its inputs from a reference ('teaching', 'error') signal. This algorithm is computationally powerful enough to be used in a very wide variety of engineering applications. However, the crucial issue is whether the external connectivity required by such applications can be implemented biologically. We argue that some applications appear to be in principle biologically implausible: these include the Smith predictor and Kalman filter (for state estimation), and the feedback-error-learning scheme for adaptive inverse control. However, even for plausible schemes, such as forward models for noise cancellation and novelty-detection, and the recurrent architecture for adaptive inverse control, there is unlikely to be a simple mapping between microzone function and internal model structure. This initial analysis suggests that cerebellar involvement in particular behaviours is therefore unlikely to have a neat classification into categories such as 'forward model'. It is more likely that cerebellar microzones learn a task-specific adaptive-filter operation which combines a number of signal-processing roles.

Dynamics of primate oculomotor plant revealed by effects of abducens microstimulation

Despite their importance for deciphering oculomotor commands, the mechanics of the extraocular muscles and orbital tissues (oculomotor plant) are poorly understood. In particular, the significance of plant nonlinearities is uncertain. Here primate plant dynamics were investigated by measuring the eye movements produced by stimulating the abducens nucleus with brief pulse trains of varying frequency. Statistical analysis of these movements indicated that the effects of stimulation lasted about 40 ms after the final pulse, after which the eye returned passively toward its position before stimulation. Behavior during the passive phase could be approximated by a linear plant model, corresponding to Voigt elements in series, with properties independent of initial eye position. In contrast, behavior during the stimulation phase revealed a sigmoidal relation between stimulation frequency and estimated steady-state tetanic tension, together with a frequency-dependent rate of tension increase, that appeared very similar to the nonlinearities previously found for isometric-force production in primate lateral rectus muscle. These results suggest that the dynamics of the oculomotor plant have an approximately linear component related to steady-state viscoelasticity and a nonlinear component related to changes in muscle activation. The latter may in part account for the nonlinear relations observed between eye-movement parameters and single-unit firing patterns in the abducens nucleus. These findings point to the importance of recruitment as a simplifying factor for motor control with nonlinear plants.

Cerebellar-inspired adaptive control of a robot eye actuated by pneumatic artificial muscles

In this paper, a model of cerebellar function is implemented and evaluated in the control of a robot eye actuated by pneumatic artificial muscles. The investigated control problem is stabilization of the visual image in response to disturbances. This is analogous to the vestibuloocular reflex (VOR) in humans. The cerebellar model is structurally based on the adaptive filter, and the learning rule is computationally analogous to least-mean squares, where parameter adaptation at the parallel fiber/Purkinje cell synapse is driven by the correlation of the sensory error signal (carried by the climbing fiber) and the motor command signal. Convergence of the algorithm is first analyzed in simulation on a model of the robot and then tested online in both one and two degrees of freedom. The results show that this model of neural function successfully works on a real-world problem, providing empirical evidence for validating: 1) the generic cerebellar learning algorithm; 2) the function of the cerebellum in the VOR; and 3) the signal transmission between functional neural components of the VOR.

Effects of initial eye position on saccade-related behavior of abducens nucleus neurons in the primate

Previous work suggests that when the eye starts at different orbital initial positions (IPs), the saccade control system is faced with significant nonlinearities. Here we studied the effects of IP on saccade-related firing of monkey abducens neurons by either isolating saccade variables behaviorally or applying a multiple linear regression analysis. Over a 50 degrees range of IPs, we could select 10 degrees horizontal saccades with identical velocity profiles, which would require identical control signals in a linear system. The bursts accompanying ipsiversive saccades for IPs above the threshold for steady firing were quite similar. The excess burst rate above steady firing was either constant or decreased with ipsiversive IP, and both the number of excess spikes in the burst and burst duration were nearly constant. However, for ipsiversive saccades from IPs below threshold, both peak burst rate (6.82 +/- 1.38 spikes.s(-1).deg(-1)) and burst duration (0.67 +/- 0.28 ms/deg) increased substantially with ipsiversive IPs. Moreover, the pause associated with contraversive saccades shortened considerably with ipsiversive IPs (mean 1.2 ms/deg). This pattern of results for pauses and for bursts below threshold suggests the presence of a significant nonlinearity. Abducting saccades are produced by the net force of agonist lateral rectus (LR) and antagonist medial rectus (MR) muscles. We suggest that the decreasing force in the MR muscle with IPs in the abducting direction requires a more vigorous burst in LR motoneurons, which appears to be generated by a combination of saturating and nonsaturating burst commands and the recruitment of additional abducens neurons.

Response linearity determined by recruitment strategy in detailed model of nictitating membrane control

Many models of eyeblink conditioning assume that there is a simple linear relationship between the firing patterns of neurons in the interpositus nucleus and the time course of the conditioned response (CR). However, the complexities of muscle behaviour and plant dynamics call this assumption into question. We investigated the issue by implementing the most detailed model available of the rabbit nictitating membrane response (Bartha and Thompson in Biol Cybern 68:135-143, 1992a and in Biol Cybern 68:145-154, 1992b), in which each motor unit of the retractor bulbi muscle is represented by a Hill-type model, driven by a non-linear activation mechanism designed to reproduce the isometric force measurements of Lennerstrand (J Physiol 236:43-55, 1974). Globe retraction and NM extension are modelled as linked second order systems. We derived versions of the model that used a consistent set of SI units, were based on a physically realisable version of calcium kinetics, and used simulated muscle cross-bridges to produce force. All versions showed similar non-linear responses to two basic control strategies. (1) Rate-coding with no recruitment gave a sigmoidal relation between control signal and amplitude of CR, reflecting the measured relation between isometric muscle force and stimulation frequency. (2) Recruitment of similar strength motor units with no rate coding gave a sublinear relation between control signal and amplitude of CR, reflecting the increase in muscle stiffness produced by recruitment. However, the system response could be linearised by either a suitable combination of rate-coding and recruitment, or by simple recruitment of motor units in order of (exponentially) increasing strength. These plausible control strategies, either alone or in combination, would in effect present the cerebellum with the simplified virtual plant that is assumed in many models of eyeblink conditioning. Future work is therefore needed to determine the extent to which motor neuron firing is in fact linearly related to the nictitating membrane response.

Consistency of Listing's law and reciprocal innervation with pseudo-inverse control of eye position in 3-D

Pseudo-inverse kinematics, under which small movements are produced by the least possible sum square changes in motor command, has been proposed as a unifying principle for the elimination of redundancy in general biological motor control systems (Pellionsz 1984) and in particular in the oculomotor system (Daunicht 1988, 1991). We have noted elsewhere (Dean et al, 1999) that this principle is incomplete without first specifying a parameterisation of motor command space and we proposed that the relevant motor-command parameter is summed motor unit firing rate. Under this assumption we were able to show that pseudo-inverse control of the horizontal extraocular muscles is consistent with available motor pool firing rate data. In this paper we extend this result to three dimensions and six extraocular muscles, showing that pseudo-inverse control is consistent with published firing rate date for a realistic model of oculomotor kinematics. We suggest that pseudo-inverse control may represent a common currency for modular control of many degree of freedom systems while its implementation may be a consequences of the minimisation of a more ecologically relevant parameter such as post-saccadic retinal slip.

A simple control law generates Listing's positions in a detailed model of the extraocular muscle system

The neural commands for maintaining static Listing's positions were identified using a detailed model of extraocular muscle based on Miller and Shamaeva (Orbit 1.5 gaze mechanics simulation 1995). The commands were approximately separable, suggesting a simple control law whereby independent horizontal and vertical commands are combined to generate tertiary positions. Tests showed that this control law (i) generated Listing' s positions to reasonable accuracy over+/-30 deg, provided pulleys were represented in the model; (ii) if driven by retinal coordinates, produced errors close to the theoretical minimum for a commutative system. The proposed commands appear consistent with electrophysiological evidence.