Person:

Michmizos, Konstantinos P.

Little is known about whether our knowledge of how the central nervous system controls the upper extremities (UE), can generalize, and to what extent to the lower limbs. Our continuous efforts to design the ideal adaptive robotic therapy for the lower limbs of stroke patients and children with cerebral palsy highlighted the importance of analyzing and modeling the kinematics of the lower limbs, in general, and those of the ankle joints, in particular. We recruited 15 young healthy adults that performed in total 1,386 visually evoked, visually guided, and target-directed discrete pointing movements with their ankle in dorsal–plantar and inversion–eversion directions. Using a non-linear, least-squares error-minimization procedure, we estimated the parameters for 19 models, which were initially designed to capture the dynamics of upper limb movements of various complexity. We validated our models based on their ability to reconstruct the experimental data. Our results suggest a remarkable similarity between the top-performing models that described the speed profiles of ankle pointing movements and the ones previously found for the UE both during arm reaching and wrist pointing movements. Among the top performers were the support-bounded lognormal and the beta models that have a neurophysiological basis and have been successfully used in upper extremity studies with normal subjects and patients. Our findings suggest that the same model can be applied to different “human” hardware, perhaps revealing a key invariant in human motor control. These findings have a great potential to enhance our rehabilitation efforts in any population with lower extremity deficits by, for example, assessing the level of motor impairment and improvement as well as informing the design of control algorithms for therapeutic ankle robots.

Distributed cortical solutions of magnetoencephalography (MEG) and electroencephalography (EEG) exhibit complex spatial and temporal dynamics. The extraction of patterns of interest and dynamic features from these cortical signals has so far relied on the expertise of investigators. There is a definite need in both clinical and neuroscience research for a method that will extract critical features from high-dimensional neuroimaging data in an automatic fashion. We have previously demonstrated the use of optical flow techniques for evaluating the kinematic properties of motion field projected on non-flat manifolds like in a cortical surface. We have further extended this framework to automatically detect features in the optical flow vector field by using the modified and extended 2-Riemannian Helmholtz–Hodge decomposition (HHD). Here, we applied these mathematical models on simulation and MEG data recorded from a healthy individual during a somatosensory experiment and an epilepsy pediatric patient during sleep. We tested whether our technique can automatically extract salient dynamical features of cortical activity. Simulation results indicated that we can precisely reproduce the simulated cortical dynamics with HHD; encode them in sparse features and represent the propagation of brain activity between distinct cortical areas. Using HHD, we decoded the somatosensory N20 component into two HHD features and represented the dynamics of brain activity as a traveling source between two primary somatosensory regions. In the epilepsy patient, we displayed the propagation of the epileptic activity around the margins of a brain lesion. Our findings indicate that HHD measures computed from cortical dynamics can: (i) quantitatively access the cortical dynamics in both healthy and disease brain in terms of sparse features and dynamic brain activity propagation between distinct cortical areas, and (ii) facilitate a reproducible, automated analysis of experimental and clinical MEG/EEG source imaging data.