Using Constrained Square-Root Cubature Kalman Filter for Quantifying the Severity of Epileptic Activities in Mice.

(1) Background: Quantification of severity of epileptic activities, especially during electrical stimulation, is an unmet need for seizure control and evaluation of therapeutic efficacy. In this study, a parameter ratio derived from constrained square-root cubature Kalman filter (CSCKF) was formulated to quantify the excitability of local neural network and compared with three commonly used indicators, namely, band power, Teager energy operator, and sample entropy, to objectively determine their effectiveness in quantifying the severity of epileptiform discharges in mice. (2) Methods: A set of one normal and four types of epileptic EEGs was generated by a mathematical model. EEG data of epileptiform discharges during two types of electrical stimulation were recorded in 20 mice. Then, EEG segments of 5 s in length before, during and after the real and sham stimulation were collected. Both simulated and experimental data were used to compare the consistency and differences among the performance indicators. (3) Results: For the experimental data, the results of the four indicators were inconsistent during both types of electrical stimulation, although there was a trend that seizure severity changed with the indicators. For the simulated data, when the simulated EEG segments were used, the results of all four indicators were similar; however, this trend did not match the trend of excitability of the model network. In the model output which retained the DC component, except for the CSCKF parameter ratio, the results of the other three indicators were almost identical to those using the simulated EEG. For CSCKF, the parameter ratio faithfully reflected the excitability of the neural network. (4) Conclusion: For common EEG, CSCKF did not outperform other commonly used performance indicators. However, for EEG with a preserved DC component, CSCKF had the potential to quantify the excitability of the neural network and the associated severity of epileptiform discharges.

Morphological and hemodynamic changes of sciatic nerves and their vasa nervorum during circular compression and relaxation.

The main aim of this study was to evaluate the biomechanical and hemodynamic responses of vasa nervorum under transverse circular compression. In situ compress-and-hold experiments were performed on the sciatic nerves of healthy and diabetic rats, and the blood flow within the vasa nervorum was observed using Doppler-optical coherence tomography. A new technique was developed to obtain the time-course of the cross sectional area and the morphology of the vasa nervorum from the tomographic images. A quasi-linear viscoelastic model was used to investigate the overall biomechanical properties of the nerves, and a two-dimensional three-layered finite element model was constructed to analyze the distribution of stress and the morphological changes during the compression-relaxation process. The results showed that the lumenal area of vasa nervorum was reduced in the compression stage, especially for the diabetic nerves. The reduction was greater than 70% when the reduction of the nerve diameter was only 10%. The quasi-linear viscoelastic model showed that normal nerves were more elastic but less viscous than the diabetic nerves. The finite element analyses demonstrated that perineurium could sustain more stress than other layers, while epineurium served as a cushion to protect vasa nervora. In addition, there were regions within epineurium with less stress, so that vasa nervora in these saddle regions were less deformed. The vasa nervorum in diabetic rats was more prone to compression and reduction of blood flow than that of the normal rats. The histological studies supported the simulation results.

Cytotoxic and biomechanical effects of clinical dosing schemes of paclitaxel on neurons and cancer cells.

PURPOSE: Chemotherapy-induced peripheral neuropathy often results in a reduction in drug dose. However, the serum level of anticancer drugs varies with time after intravenous infusion, and this factor has seldom been considered in previous in vitro studies. The goals of this study were to build an automatic dosage control system and to evaluate the influence of drug infusion rate on the cells. METHODS: Neurons and melanoma cells were used as the samples. The 3-h (average and peak concentration: 0.024 and 0.287 µM) and 24-h infusion (average and peak concentration: 0.020 and 0.042 µM) schemes were investigated. For evaluations, cell indentation tests by an atomic force microscope, serial immunofluorescent images, and cell viability analysis was performed. RESULTS: For the neurons, Young's modulus first increased and then remained unchanged in the 3-h scheme, but was stationary throughout the observation period in the 24-h scheme. For the cancer cells, Young's modulus increased in both infusion schemes, and the increase was larger in the 3-h scheme. Morphologically, axons swelled and shortened, and the number of their branches decreased in the 3-h scheme. In contrast, there was only slowed growth of axons without obvious morphological changes in the 24-h scheme. Viability analysis of the cancer cells revealed that the 3-h scheme had a better anticancer effect. CONCLUSION: A dosage-control system simulating the pharmacodynamic changes of drugs was successfully constructed for in vitro cell cultures. The 3-h scheme of paclitaxel showed better anticancer effects but more adverse effects on neuronal growth and morphology.

Flexor and extensor muscle tone evaluated using the quantitative pendulum test in stroke and parkinsonian patients.

The aim of this study was to evaluate the flexor and extensor muscle tone of the upper limbs in patients with spasticity or rigidity and to investigate the difference in hypertonia between spasticity and rigidity. The two experimental groups consisted of stroke patients and parkinsonian patients. The control group consisted of age and sex-matched normal subjects. Quantitative upper limb pendulum tests starting from both flexed and extended joint positions were conducted. System identification with a simple linear model was performed and model parameters were derived. The differences between the three groups and two starting positions were investigated by these model parameters and tested by two-way analysis of variance. In total, 57 subjects were recruited, including 22 controls, 14 stroke patients and 21 parkinsonian patients. While stiffness coefficient showed no difference among groups, the number of swings, relaxation index and damping coefficient showed changes suggesting significant hypertonia in the two patient groups. There was no difference between these two patient groups. The test starting from the extended position constantly manifested higher muscle tone in all three groups. In conclusion, the hypertonia of parkinsonian and stroke patients could not be differentiated by the modified pendulum test; the elbow extensors showed a higher muscle tone in both control and patient groups; and hypertonia of both parkinsonian and stroke patients is velocity dependent.

Finite element modeling of hyper-viscoelasticity of peripheral nerve ultrastructures.

The mechanical characteristics of ultrastructures of rat sciatic nerves were investigated through animal experiments and finite element analyses. A custom-designed dynamic testing apparatus was used to conduct in vitro transverse compression experiments on the nerves. The optical coherence tomography (OCT) was utilized to record the cross-sectional images of nerve during the dynamic testing. Two-dimensional finite element models of the nerves were built based on their OCT images. A hyper-viscoelastic model was employed to describe the elastic and stress relaxation response of each ultrastructure of the nerve, namely the endoneurium, the perineurium and the epineurium. The first-order Ogden model was employed to describe the elasticity of each ultrastructure and a generalized Maxwell model for the relaxation. The inverse finite element analysis was used to estimate the material parameters of the ultrastructures. The results show the instantaneous shear modulus of the ultrastructures in decreasing order is perineurium, endoneurium, and epineurium. The FE model combined with the first-order Ogden model and the second-order Prony series is good enough for describing the compress-and-hold response of the nerve ultrastructures. The integration of OCT and the nonlinear finite element modeling may be applicable to study the viscoelasticity of peripheral nerve down to the ultrastructural level.

Discontinuous Galerkin finite element method for solving population density functions of cortical pyramidal and thalamic neuronal populations.

Compared with the Monte Carlo method, the population density method is efficient for modeling collective dynamics of neuronal populations in human brain. In this method, a population density function describes the probabilistic distribution of states of all neurons in the population and it is governed by a hyperbolic partial differential equation. In the past, the problem was mainly solved by using the finite difference method. In a previous study, a continuous Galerkin finite element method was found better than the finite difference method for solving the hyperbolic partial differential equation; however, the population density function often has discontinuity and both methods suffer from a numerical stability problem. The goal of this study is to improve the numerical stability of the solution using discontinuous Galerkin finite element method. To test the performance of the new approach, interaction of a population of cortical pyramidal neurons and a population of thalamic neurons was simulated. The numerical results showed good agreement between results of discontinuous Galerkin finite element and Monte Carlo methods. The convergence and accuracy of the solutions are excellent. The numerical stability problem could be resolved using the discontinuous Galerkin finite element method which has total-variation-diminishing property. The efficient approach will be employed to simulate the electroencephalogram or dynamics of thalamocortical network which involves three populations, namely, thalamic reticular neurons, thalamocortical neurons and cortical pyramidal neurons.

Suppression of acute seizures by theta burst electrical stimulation of the hippocampal commissure using a closed-loop system.

This study investigated the effects of electrical stimulation with theta burst stimulation (eTBS) on seizure suppression. Optimal parameters of eTBS were determined through open-loop stimulation experiments and then implemented in a close-loop seizure control system. For the experiments, 4-aminopyridine (4-AP) was injected into the right hippocampus of Sprague-Dawley rats to induce an acute seizure. eTBS was applied on the ventral hippocampal commissure and the effects of eTBS with different combinations of burst frequency and number of pulses per burst were analyzed in terms of seizure suppression. A closed-loop seizure control system was then implemented based on optimal eTBS parameters. The efficiency of the closed-loop eTBS was evaluated and compared to that of high frequency stimulation. The results show that eTBS induced global suppression in the hippocampus and this was sustained even after the application of eTBS. The optimal parameter of eTBS in the open-loop stimulation experiments was a burst frequency at 100Hz with nine pulses in a burst. The eTBS integrated with the on-off control law yielded less actions and cumulative delivered charge, but induced longer after-effects of seizure suppression compared to continuous high frequency stimulation (cHFS). To conclude, eTBS has suppressive effects on 4-AP induced seizure. A closed-loop eTBS system provides a more effective way of suppressing seizure and requires less effort compared to cHFS. eTBS may be a novel stimulation protocol for effective seizure control.

Using a discrete preisach model for hysteresis in ankle joint passive moment.

The steady-state passive joint moment was considered as a nonlinear elasticity in the past. However, we found that it was path dependent and the estimation error could be large if the commonly used path-independent functions were adopted. The aim of this study was to develop a model to describe the movement history-dependent passive moment in the steady state. The steady-state passive ankle moments of the rabbit were measured by a series of ramp-and-hold angle changes (stairway angle trajectory). A customized discrete Preisach model was constructed and a commonly adopted double-exponential function was also implemented. Two sets of data with different angle paths (major loop and inward loop trajectories) were acquired for model validation. The performance of the two models was compared. The results showed that the proposed model could accurately estimate the steady-state passive moment for both sets of validation data. The estimated error of the proposed model was approximately 50% smaller than that of the double-exponential function approach. It is expected that this new approach, by reducing the error of estimating passive joint moment, may contribute to the active control of joint moments.

On-off control of burst high frequency electrical stimulation to suppress 4-AP induced seizures.

OBJECTIVE: The goal of this study was to investigate, using model simulations and animal experiments, the efficiency and the side effects of burst high frequency stimulation combined with on-off control in seizure suppression. APPROACH: A modified mathematical hippocampal seizure model was created to provide evidence of the eligibility of this approach. In the experimental setup, two recording electrodes were inserted into bilateral septal CA1 of the hippocampus, and a stimulation electrode was placed on the ventral hippocampal commissure of a rat. After seizures had been induced by 4-aminopyridine treatment, on-off control stimulation was used to suppress the seizures at 20 s intervals. The stimulation time, cumulative charge and post-stimulation suppression were used to assess the effects of burst duration. MAIN RESULTS: The results showed that burst stimulation could suppress the seizures during the control period and burst stimulation of a shorter duration could keep the seizure suppressed with less effort. By decreasing the burst duration, the cumulative stimulation time became shorter, the delivered cumulative charge became lower, and the cumulative time of post-stimulation suppression became longer. SIGNIFICANCE: The on-off control stimulation not only prolonged the duration of suppression but also avoided the side effects of the conversion of seizure patterns. In particular, decreasing the specified burst duration increased the efficiency of the burst stimulation.

Effects of levetiracetam on µ rhythm in persons with epilepsy.

Mu rhythm can be suppressed by movements, the so called event-related desynchronization (ERD). Levetiracetam (LEV) is a newer type of antiepileptic drug. A previous study reported that LEV might enhance mu rhythm and caused mu status in one subject. The main purpose of this study was to investigate the effects of LEV on EEG frequency contents and ERD. Seventeen patients with epileptic foci outside the Rolandic area were recruited. The following studies were performed before and after chronically taking LEV. An electroencephalogram (EEG) with 10 minutes of resting state and 5 minutes covering 10 right thumb movements were recorded. Reaction time was evaluated with a simple visual reaction time test. EEG data were analyzed by S-transformation and relative band powers were calculated. The results showed that the relative powers of theta, alpha and beta band in frontal (F3 and F4) and occipital (O1 and O2) leads and mu band in the centro-parietal (C3, C4, P3 and P4) leads were not changed by chronically taking LEV. No mu status was observed in any subject. However, the mean group ERD was enhanced at C3, Cz and P4 leads. Reaction time was similar before and after taking LEV. In conclusion, chronically taking LEV did not change the frequency contents of EEG and did not cause drowsiness, but enhanced ERD. The results suggest that chronically taking medication, such as LEV, is a plausible method to broaden the applicability of ERD-based brain-computer interfaces.

High frequency stimulation can suppress globally seizures induced by 4-AP in the rat hippocampus: an acute in vivo study.

BACKGROUND: High frequency stimulation (HFS) on the hippocampus can locally suppress epileptiform activity in-vitro and decrease seizure frequency in vivo. In-vitro HFS on the ventral commissural tract, a novel target, was shown to block the axonal conduction and suppress activity in the CA1 and CA3 neuron. OBJECTIVE: To study the spatial extent of seizure suppression by HFS applied on the tract and focus site in an in vivo experiment. METHODS: Five adult Sprague-Dawley rats were used for the study. Six electrodes were placed on the septal, middle, and temporal hippocampus bilaterally to simultaneously record seizure activity in the entire hippocampus. Seizure activity was induced by injecting 4-aminopyridine (4-AP) into the right middle part of the hippocampus. Following induction, HFS (100 Hz) was applied to the tract and the focus site at 100, 300 and 500 µA. RESULTS: The induced seizure activity was dominated by two patterns, high frequency spiking and pseudo-periodic spikes. Either tract or focus site stimulation could generate suppression of only the pseudo-periodic spikes. The suppression rates were dependent on stimulation amplitude (P < 0.005, chi square test). However, HFS also caused conversion of the seizure pattern. The conversion rates increased with higher stimulation amplitudes and were higher with focus site stimulation (P < 0.01, Fisher's exact test). CONCLUSIONS: The results of this study have two practical implications [1], both tract and focus site stimulation can produce global suppression of hippocampus and [2] the choice of stimulation parameters is critical in order to produce suppression, but not conversion, of seizure pattern.

Myelinated nerve bundles developed on the plano-concave fibers containing nerve conduit.

Morphologically and chemically modified plano-concave fibers (PCFs) are designed as a unit of guided channels for supporting Schwann cells to facilitate mass transport and promote nerve regeneration. The surface-modified PCFs are imprinted with linearly patterned grooves (LPGs) to guide adherent Schwann cell elongation and axon extension. After being cocultured with PC12 neuron-like cells, Schwann cells differentiate into the myelinated type and interact with PC12 axons. The myelinated axons aggregate as a linear bundle and extend along the direction of LPGs on a PCF. The cross section of a myelin structure is examined using a transmission electron microscope. The PCFs can potentially bridge gaps in injured nerves, improving the therapeutic efficacy of nerve regeneration.

Segmentation of nerve fibers using multi-level gradient watershed and fuzzy systems.

OBJECTIVE: This paper presents an algorithm based on multi-level watershed segmentation combined with three fuzzy systems to segment a large number of myelinated nerve fibers in microscope images. The method can estimate various geometrical parameters of myelinated nerve fibers in peripheral nerves. It is expected to be a promising tool for the quantitative assessment of myelinated nerve fibers in related research. MATERIALS AND METHODS: A novel multi-level watershed scheme iteratively detects pre-candidate nerve fibers. At each immersion level, watershed segmentation extracts the initial axon locations and obtains meaningful myelinated nerve fiber features. Thereafter, according to a priori characteristics of the myelinated nerve fibers, fuzzy rules reject unlikely pre-candidates and collect a set of candidates. Initial candidate boundaries are then refined by a fuzzy active contour model, which flexibly deforms contours according to the observed features of each nerve fiber. A final scan with a different set of fuzzy rules based on the a priori properties of the myelinated nerve fibers removes false detections. A particle swarm optimization method is employed to efficiently train the large number of parameters in the proposed fuzzy systems. RESULTS: The proposed method can automatically segment the transverse cross-sections of nerve fibers obtained from optical microscope images. Although the microscope image is usually noisy with weak or variable levels of contrast, the proposed system can handle images with a large number of myelinated nerve fibers and achieve a high fiber detection ratio. As compared to manual segmentation by experts, the proposed system achieved an average accuracy of 91% across different data sets. CONCLUSION: We developed an image segmentation system that automatically handles myelinated nerve fibers in microscope images. Experimental results showed the efficacy of this system and its superiority to other nerve fiber segmentation approaches. Moreover, the proposed method can be extended to other applications of automatic segmentation of microscopic images.

A robust algorithm for removing artifacts in EEG recorded during FMRI/EEG study.

The main purpose of this study was to propose a robust algorithm for removing artifacts from the electroencephalographic (EEG) data collected during magnetic resonance imaging (MRI). The core idea of the proposed method was to remove the main gradient artifacts by the maximum cross-correlation method and to remove the residual artifacts by the rolling-ball algorithm and lowpass filtering. The results showed that the proposed algorithm had a better performance and was robust in the sense that its performance was maintained when the sampling rate of EEG data was decreased from 10KHz to 200Hz.

Enhanced Schwann cell adhesion and elongation on a topographically and chemically modified poly(L-lactic acid) film surface.

A dense poly-L-lactic acid (PLLA) film was employed as the primary material and hot-embossed with the formation of microgrooves (g-PLLA). A thin layer of Au was then deposited on the film to obtain a morphologically modified substrate (Au/g-PLLA). The Au/g-PLLA film surface was then chemically modified by imprinting octadecanethiolate (ODT) self-assembled monolayers on the upper surface (ODT/Au/g-PLLA), followed by Arg-Gly-Asp (RGD) peptide sequences on the microgrooves (RGD_ODT/Au/g-PLLA). The surface chemistry of the as-prepared RGD_ODT/Au/g-PLLA samples was examined. The bioactivity and spreading function of Schwann cells cultured on the morphologically and chemically modified surfaces were assessed. The results demonstrate that Schwann cells adhered to the RGD/Au/g-PLLA surface and proliferated along the microgrooves without crossing over the ODT/Au/PLLA surface. The proposed film surface can be used for manipulating the outgrowth of axons by modifying and arranging a selected region to induce cell growth and to prevent cells from spreading out nondirectionally.

Description and computational modeling of the whole course of status epilepticus induced by low dose lithium-pilocarpine in rats.

The main purposes of this study were (1) to describe the whole course of status epilepticus induced by a low dose lithium-pilocarpine model in rats, including depth-EEG from the hippocampus, ECoG from cortex and gross behaviors, and (2) to investigate the possible changes of the intrinsic neural network in the hippocampus during the status epilepticus by model simulation. The course of the induced status epilepticus was divided into baseline, pre-ictal, episodic, onset, continuous, and convalescence stages. At the pre-ictal stage, the main component of the spectrum of the depth-EEG shifted before seizure activity first appeared in the hippocampus at the episodic stage and propagated to the cortex at the onset stage. Model simulation indicated that the changes of depth-EEG of the whole course could be simulated by changes of the loop gains and the inputs of the neural network. The excitatory and inhibitory loops were first enhanced in the pre-ictal stage. The isolated seizures in the episodic stage were caused by variations of inputs from other pyramidal cells. The gain of slow inhibitory loop gradually decreased, leading to the onset of status epilepticus. Different types of waveforms in the status epilepticus corresponded to the different ratio of excitatory and inhibitory influence, caused by different levels of inhibition. Finally, in the convalescence stage, the gains of the excitatory loop and the inhibitory loop both changed toward the base line, recovering the balance. The changes of the parameters in the model were compatible with the experimental results in the literature.

Depth-sensing nano-indentation on a myelinated axon at various stages.

A nano-mechanical characterization of a multi-layered myelin sheath structure, which enfolds an axon and plays a critical role in the transmission of nerve impulses, is conducted. Schwann cells co-cultured in vitro with PC12 cells for various co-culture times are differentiated to form a myelinated axon, which is then observed using a transmission electron microscope. Three major myelination stages, with distinct structural characteristics and thicknesses around the axon, can be produced by varying the co-culture time. A dynamic contact module and continuous depth-sensing nano-indentation are used on the myelinated structure to obtain the load-on-sample versus measured displacement curve of a multi-layered myelin sheath, which is used to determine the work required for the nano-indentation tip to penetrate the myelin sheath. By analyzing the harmonic contact stiffness versus the measured displacement profile, the results can be used to estimate the three stages of the multi-layered structure on a myelinated axon. The method can also be used to evaluate the development stages of myelination or demyelination during nerve regeneration.

In vivo impedance evaluation of Au/PI microelectrode with surface modulated by alkanethiolate self-assembled monolayers.

The goal of this study was to verify that a fully implanted microelectrode with modulated surface may have a reduced rising rate of total impedance and a longer life time. In the previous work, alkanethiolate self-assembled monolayers (SAMs) surface as protein-resistant spacer or cell-repulsive dense-packed spacer has been verified from in vitro experiments. In this study, microelectrodes with the same surface modulation were implanted into the subcutaneous layers of Wistar rats. Nine rats were implanted with the microelectrodes and the total impedance data were measured every 24 h for 2 weeks after implantation. An equivalent electrical circuit model of the electrode-tissue interface was established and parameters were estimated by using an optimization algorithm. Four out of nine rats had manifested acute inflammation reaction and the rests revealed only slight tissue response. Histological examination for the inflammatory group showed fibroblasts, macrophages, and polymorphonuclear leukocytes in adjacent to the electrode contact surface. In the inflammatory group, no significantly difference in total impedance was found in both types of electrodes. However, the trend of total impedance of SAMs-treated electrodes could maintain a steady state value after 1 week. For the non-inflammatory group, both types of electrodes could reduce the impedance value within implanted days. The tissue resistance might be related to the thickness of cells adhered upon the electrode contacts.

Estimation of ankle joint angle from peroneal and tibial electroneurograms based on muscle spindle model.

The main goal of this study was to develop a new method of estimating the angle of the passively stretched ankle joint, based on structural muscle spindle models of the tibial and peroneal electroneurograms (ENG). Passive ramp-and-hold and alternating stretches of the ankle joint were performed in a rabbit. Simultaneously, two cuff electrodes were used to record the ENGs of peroneal and tibial nerves. Based on the two ENGs and the joint angle trajectory, two muscle spindle models were constructed and their inverse models were integrated to compute angle estimates. The model parameters were optimized. The performance of our approach was compared with those of the adaptive neuro-fuzzy inference system and artificial neural network model. The results revealed that our model had a better performance of estimating the ankle joint angle in large-range movements and smaller tracking errors. This study provides a new estimation algorithm to extract the joint angle from the information conveyed in a nerve.

Neuro-rehabilitation robot-assisted assessments of synergy patterns of forearm, elbow and shoulder joints in chronic stroke patients.

BACKGROUND: Abnormal synergy is one of the major motor deficits in stroke patients. Abnormal muscle synergies, in conjunction with weakness and spasticity, interfere with voluntary movements and restrict the range of motion. This study aimed to quantify abnormal synergies in the affected upper limbs of chronic stroke patients by using a neuro-rehabilitation robot. METHODS: Twelve chronic stroke patients and eight age-matched control subjects were recruited to perform rectilinear tracking movements in four horizontal directions (back-forth, two oblique directions at 45 degrees , and right-left). Kinematic, kinetic and electromyogram data were recorded and used to develop two biomechanical indices and one electromyogram assessment index based on principal component analysis. FINDINGS: Significant differences between upper limbs of control subjects and the affected side of stroke patients were observed in all three assessment indices. Higher correlation between the elbow joint angle and the forearm pronation/supination torque, higher variation of the forearm torque, and abnormal co-contraction of the elbow and shoulder muscles were observed in the affected limbs of stroke patients. The difference was more prominent in the right-left direction and the oblique direction contra-proximal to ipsi-distal. INTERPRETATION: The proposed assessment indices could be employed to quantify the abnormal synergies in stroke patients. Rectilinear tracking along the right-left direction and the oblique direction of contra-proximal to ipsi-distal is more suitable for assessing abnormal synergies.