-
A Unified Dynamic Model for the Decomposition of Skin Conductance and the Inference of Sudomotor Nerve Activities
Hui S. WangStacy MarsellaMisha Pavel
Keywords:Mathematical modelsSkinBiological system modelingDynamical systemsPhysiologyEstimationPsychologySystem identificationHeuristic algorithmsAdaptation modelsDynamic ModelSudomotorSystem DynamicsSystem IdentificationLow ArousalAspects Of StressArousal IndexSparse OptimizationLinear Time-invariant ModelTime ConstantLatent VariablesModel In OrderFixed PointExpectation MaximizationPhase ResponseCubic SplineSystem OutputLong-term ResponseSparse RepresentationLevels In The AbsenceTonic LevelsFree ResponseForced ResponseComplex PolesBi-exponential FunctionShort-term ResponseSparse SignalNeutral ImagesLinear Time-invariant SystemsSparsity LevelBiomedical signal processingdynamic systemsoptimizationphysiologysparse representationsympathetic nervous systemsystem identification
Abstracts:Electrodermal activity (EDA), commonly measured as skin conductance (SC), is a widely used physiological signal in psychological research and behavioral health applications. EDA is considered an indicator of arousal, a key aspect of emotion and stress. This work proposes a data-driven dynamic system model that characterizes the temporal dynamics of skin conductance and infers the latent arousal signal, utilizing techniques from system identification and sparse optimization. It introduces a fourth-order, linear time-invariant model for the overall skin conductance signal, including both the tonic and phasic components. The model was applied to a large dataset of over 200 participants to evaluate model fit. Furthermore, a three-component decomposition of skin conductance is introduced, based on mathematical definitions derived from the model, which provides key insights into the temporal dynamics of skin conductance. Comparative evaluation shows that the estimated latent neural signal effectively differentiates between high and low arousal states, while maintaining expected physiological properties. This work lays the foundation for numerous behavioral health applications and paves the road for designing physiology-based interventions aimed at regulating arousal.
-
Cortical Activation Patterns Determine Effectiveness of rTMS-Induced Motor Imagery Decoding Enhancement in Stroke Patients
Tianyu JiaLinhong MoCiarán McGeadyJingyao SunAixian LiuLinhong JiJianing XiChong Li
Keywords:MotorsStroke (medical condition)Medical treatmentLesionsElectroencephalographyDecodingNeuromodulationBiomedical imagingBiomedical engineeringTranscranial magnetic stimulationActivity PatternsStroke PatientsCortical ActivityMotor ImageryCortical Activation PatternsEffects Of Motor ImageryMotor Imagery DecodingPatients In GroupIndividual DifferencesRecovery RateTranscranial Magnetic StimulationNeural PatternsRepetitive Transcranial Magnetic StimulationEffect Of Combination TherapyCortical LesionsDecoding PerformanceFugl-Meyer AssessmentNeural ReorganizationHz rTMSNeuromodulation TherapyEvent-related DesynchronizationEffects Of rTMSContralesional HemisphereNon-invasive NeuromodulationSubcortical LesionsNeural ActivitySensorimotor AreasSimilar PhenomenonBrain ActivityMotor RecoveryMotor imageryneurorehabilitationstroketranscranial magnetic stimulation
Abstracts:Combination therapy with motor imagery (MI)-based brain-computer interface (BCI) and repetitive transcranial magnetic stimulation (rTMS) is a promising therapy for poststroke neurorehabilitation. However, with patients’ individual differences, the clinical effects vary greatly. This study aims to explore the hypothesis that stroke patients show individualized cortical response to rTMS treatments, which determine the effectiveness of rTMS-induced MI decoding enhancement. We applied four kinds of rTMS treatments respectively to four groups of subacute stroke patients, twenty-six patients in total, and observed their EEG dynamics, MI decoding performance, and Fugl-Meyer assessment changes following 2-week neuromodulation. Four treatments consisted of ipsilesional 10 Hz rTMS, contralesional 1 Hz rTMS, ipsilesional 1 Hz rTMS, and sham stimulation. Results showed stroke patients with different neural reorganization patterns responded individually to rTMS therapy. Patients with cortical lesions mostly showed contralesional recruitment and patients without cortical lesions mostly presented ipsilesional focusing. Significant activation increases in the ipsilesional hemisphere (pre: −15.7% ∓ 8.2%, post: −17.3% ∓ 8.1%, p = 0.037) and MI decoding accuracy enhancement (pre: 76.3 ± 13.8%, post: 86.6 ± 8.2%, p = 0.037) were concurrently found in no-cortical-lesion patients with ipsilesional activation treatment. In the group of patients without cortical lesions, recovery rate in those receiving ipsilesional activation therapy (23.5 ± 10.4%) was higher than those receiving ipsilesional suppression therapy (9.9 ± 9.3%) (p = 0.041). This study reveals that tailoring neuromodulation therapy by recognizing cortical activation patterns is promising for improving effectiveness of the combination therapy with BCI and rTMS.
-
A Growing Bubble Speller Paradigm for Brain–Computer Interface Based on Event-Related Potentials
Jing JinXueqing ZhaoIan DalyShurui LiXingyu WangAndrzej CichockiTzyy-Ping Jung
Keywords:Image color analysisElectroencephalographyVisualizationFeature extractionDecodingFacesElectric potentialBrain-computer interfacesSupport vector machinesConvolutional neural networksParietal CortexClassification PerformanceGlobal ModeCentral RegionConvolutional Neural NetworkSupport Vector MachineFeature MapsLinear Discriminant AnalysisNetwork ParametersAdam OptimizerOccipital LobeAccuracy ScoresEarly StoppingRed SquareVisual Evoked PotentialsAverage Pooling LayerEvent-related Potentials AnalysisStimulus Onset AsynchronyBinomial CoefficientStimulus PatternBrain-computer Interface SystemExponential Linear UnitBrain-computer interface (BCI)event-related potential (ERP)speller paradigmgrowing bubblemultiple windows
Abstracts:Objective: Event-related potentials (ERPs) reflect electropotential changes within specific cortical regions in response to specific events or stimuli during cognitive processes. The P300 speller is an important application of ERP-based brain-computer interfaces (BCIs), offering potential assistance to individuals with severe motor disabilities by decoding their electroencephalography (EEG) to communicate. Methods: This study introduced a novel speller paradigm using a dynamically growing bubble (GB) visualization as the stimulus, departing from the conventional flash stimulus (TF). Additionally, we proposed a “Lock a Target by Two Flashes” (LT2F) method to offer more versatile stimulus flash rules, complementing the row and column (RC) and single character (SC) modes. We applied the “Sub and Global” multi-window mode to EEGNet (mwEEGNet) to enhance classification and explored the performance of eight other representative algorithms. Results: Twenty healthy volunteers participated in the experiments. Our analysis revealed that our proposed pattern elicited more pronounced negative peaks in the parietal and occipital brain regions between 200 ms and 230 ms post-stimulus onset compared with the TF pattern. Compared to the TF pattern, the GB pattern yielded a 2.00% increase in online character accuracy (ACC) and a 5.39 bits/min improvement in information transfer rate (ITR) when using mwEEGNet. Furthermore, results demonstrated that mwEEGNet outperformed other methods in classification performance. Conclusion and Significance: These results underscore the significance of our work in advancing ERP-based BCIs.
-
Gait Symmetric Adaptation and Aftereffect Through Concurrent Split-Belt Treadmill Walking and Explicit Visual Feedback Distortion
Seung-Jae KimOmik SaveEmily TannerArianna MarquezHyunglae Lee
Keywords:Legged locomotionVisualizationDistortionTrainingBarsBeltsMotorsProtocolsPerturbation methodsMotion captureVisual FeedbackVisual HallucinationsTreadmill WalkingExplicit FeedbackSplit-belt TreadmillStep LengthMotor LearningAdaptation PeriodGait TrainingTreadmill BeltBelt SpeedGait SymmetryWalking SpeedMotion CaptureLeft LegTraining ModalitiesEffective RehabilitationGait PatternAsymmetric PatternImplicit LearningPreferred SpeedExplicit StrategyBar HeightChange In SymmetryVisuomotor AdaptationHeel StrikeImplicit ProcessesCombination TrialsDecay Behavior10-minute PeriodGait adaptationgait symmetrymotor retentionsplit-belt trainingvisual feedback distortion
Abstracts:Objective: Gait asymmetry can be improved with various gait training methods. Combining split-belt treadmill walking (SB) with visual feedback distortion (VD) could enhance motor learning, improving gait symmetry adaptation and retention. This study compared step length symmetry adaptation and aftereffects between SB-only and combined explicit VD with SB, as well as between explicit VD-only and combined explicit VD with SB. Method: The 28-minute trials included three phases: a 3-minute baseline, a 10-minute adaptation, and a 15-minute post-adaptation. In the VD trial, two bars representing step lengths were displayed, with the right bar gradually decreasing by 3% to prompt participants to consciously correct their steps to match the heights of the two bars. In the SB trial, the right treadmill belt speed was incrementally increased by 5%. The VD+SB trial combined both perturbations. After the removal of these perturbations, the aftereffect of the adapted asymmetric step length was evaluated in the post-adaptation period. Results: During the adaptation period, the step length symmetry ratio shifted negatively in the SB trial, while it increased positively in the VD trial, indicating longer right steps than left. In the VD+SB trial, subjects extended their right step more than their left. Notably, the VD+SB trial demonstrated a longer aftereffect compared to both the SB-only and VD-only trials. Conclusion: The visual distortion paradigm can be explicitly applied and integrated with split-belt treadmill walking to enhance the efficacy of symmetric gait adaptation, resulting in more sustained effects on the retention of newly learned motor patterns.
-
Label-Free Non-Contact Vascular Imaging Using Photon Absorption Remote Sensing
James A. Tummon SimmonsSarah J. WerezakBenjamin R. EcclestoneJames E. D. TweelHager GaoudaParsin Haji Reza
Keywords:ImagingMicroscopyOptical imagingAbsorptionScatteringIn vivoOptical fibersOptical scatteringMirrorsLaser excitationRemote SensingVascular ImagingNon-contact ImagingRed Blood CellsDisease PreventionImaging TechniquesFunctional ImagingChicken EmbryosAlignment MethodCapillary NetworkSignal ExtractionMouse EarBlood OxygenLaser ExcitationLarge VesselGraphical User InterfaceCarbon FiberExcitatory EffectsOptical ScanningAvalanche PhotodiodeRadiative RelaxationNon-radiative RelaxationAbsorption ContrastFrames Per Second2D CorrelationNon-radiative PathwaysCapillary StructureChorioallantoic MembraneCondenser LensFirst-order CorrectionLabel free microscopyvascular imagingfunctional imagingphoton absorption remote sensing (PARS)
Abstracts:Objective: Functional vascular imaging is a critical method for early detection and prevention of disease. Established non-contact vascular imaging techniques capture predominantly structural information. In this study, a novel non-contact label-free in vivo Photon Absorption Remote Sensing (PARS) microscope is developed for structural and functional vascular imaging. Methods: The presented in vivo PARS microscope captures the endogenous absorption of green (532nm) light to form a complete picture of vasculature and surrounding tissues. Imaging system repeatability is enhanced through robust transient absorption signal extraction, and state-of-the-art real-time alignment methods. Results: Detailed imaging of vascular structure is demonstrated through in vivo microscopy of two established animal models: mouse ear and chicken embryo. Preliminary functional contrast is realized through video rate imaging of red blood cell dynamics in the capillary networks of chicken embryos. Conclusion: The presented in vivo PARS microscope successfully captures detailed structural and functional vascular contrast. Significance: This innovative non-contact label-free imaging technique holds promise as a tool for preventative medical care, as functional change often precedes structural change.
-
Ventricular Arrhythmia Classification Using Similarity Maps and Hierarchical Multi-Stream Deep Learning
Qing LinDino OglićMichael J. CurtisHak-Keung LamZoran Cvetković
Keywords:ElectrocardiographyArrhythmiaHeart beatFeature extractionDatabasesAccuracyFiltersSensitivityFibrillationConvolutionDeep LearningVentricular ArrhythmiasSimilarity MapArrhythmia ClassificationClassification AccuracyCardiac ArrestSudden DeathSudden Cardiac DeathSecondary LossVentricular TachycardiaTranslational MedicineVentricular FibrillationAverage SensitivityCorrect DetectionTraining SetTraining DataDifferences In DistributionWindow SizeInput FeaturesCharacteristic SignalsECG SignalsQRS ComplexMaximum Mean DiscrepancyCyclic ShiftT-distributed Stochastic Neighbor EmbeddingTask AccuracyWaveform FeaturesSecond DerivativeLow-dimensional FeatureFeature ChannelsCardiac arrhythmiashierarchical classificationresidual convolutional neural networksventricular fibrillationventricular tachycardia
Abstracts:Objective: Ventricular arrhythmias are the primary arrhythmias that cause sudden cardiac death. We address the problem of classification between ventricular tachycardia (VT), ventricular fibrillation (VF) and non-ventricular rhythms (NVR). Methods: To address the challenging problem of the discrimination between VT and VF, we develop similarity maps – a novel set of features designed to capture regularity within an ECG trace. These similarity maps are combined with features extracted through learnable Parzen band-pass filters and derivative features to discriminate between VT, VF, and NVR. To combine the benefits of these different features, we propose a hierarchical multi-stream ResNet34 architecture. Results: Our empirical results demonstrate that the similarity maps significantly improve the accuracy of distinguishing between VT and VF. Overall, the proposed approach achieves an average class sensitivity of 89.68%, and individual class sensitivities of 81.46% for VT, 89.29% for VF, and 98.28% for NVR. Conclusion: The proposed method achieves a high accuracy of ventricular arrhythmia detection and classification. Significance: Correct detection and classification of ventricular fibrillation and ventricular tachycardia are essential for effective intervention and for the development of new therapies and translational medicine.
-
LMI-Enabled Absolutely Stabilizing PID Control of Pharmacological Systems for Closed-Loop Automated Intravenous Drug Administration
Weidi YinDrew X. HohenhausAli TivayRajesh RajamaniJin-Oh Hahn
Keywords:DrugsPD controlControl designPI controlAnesthesiaStability criteriaLinear matrix inequalitiesUncertaintyLyapunov methodsBiomedical engineeringDrug AdministrationIntravenous AdministrationInjection DrugDose-dependentStability Of SystemControl DesignParameter UncertaintyClosed-loop SystemAsymptotically StableMatrix InequalitiesExtensive EvaluationClosed-loop ControlLinear Matrix InequalitiesPlant DynamicsClosed-loop Control SystemDynamic UncertaintiesNon-linear Dose ResponseStability Of Control SystemAbsolute StabilityNon-linear Dose-response RelationshipController Design MethodLinear DynamicsDesign MethodConcentration Of PropofolInter-individual VariabilityVirtual PatientsPositive Definite MatrixSystem In FigDC GainSymmetric Positive Definite MatrixLinear matrix inequalityabsolute stabilitypharmacological systemsintravenous druganesthesia
Abstracts:Objective: We developed a linear matrix inequality-enabled absolutely stabilizing proportional-integral-derivative control design approach for pharmacological systems applicable to intravenous drug administration. Methods: We developed a proportional-integral- derivative control design approach that does not require detailed knowledge of the dose-response relationship other than its sector bound. It repetitively solves a set of linear matrix inequalities, which encapsulate the Lyapunov stability conditions against unknown dose-response relationship, over a broad proportional-integral-derivative gain space. The linear matrix inequality-feasible proportional-integral-derivative gains guarantee the absolute stability of the closed-loop control system against unknown yet sector-bounded dose-response relationship. The proof-of-concept of the approach was shown in silico using intravenous propofol anesthesia as a practical case scenario. Results: The in silico evaluation results demonstrated the robustness and performance of the proportional-integral- derivative controllers designed with the proposed control design approach against unknown sector-bounded nonlinear dose-response relationship and parametric uncertainty in the plant dynamics. Conclusion: Pending follow-up development and extensive evaluation in various complex intravenous drug administration problems, the proposed approach may find applications in various closed-loop automated intravenous drug administration problems with complex and highly nonlinear dose-response relationships. Significance: The proposed control design approach provides a systematic way to absolutely stabilize pharmacological systems against unknown, nonlinear, and time- varying dose-response relationship, perhaps for the first time.
-
A Computational Study on the Activation of Neural Transmission in Deep Brain Stimulation
Arash GolmohammadiJan Philipp PayonkUrsula van RienenRevathi Appali
Keywords:AxonsComputational modelingMathematical modelsBrain modelingCapacitanceOptical fiber cablesMembrane potentialsDeep brain stimulationTrajectoryTemperature measurementElectrodeStimulation ProtocolNotion Of ActionPrediction ModelAction PotentialCross-sectional AreaBasal GangliaPulse WidthFirst-principlesMyelin SheathLateral SurfaceSource TermAction Potential GenerationTransmembrane PotentialInternal CapsuleNeuronal DynamicsVolume ConductionGreat CircleNeuronal CommunicationAntidromicAxoplasmAxonal ResponseRefractory StatusFiber BundlesSpecific CapacityActivation functionAxonal communicationcomputational modelingcable modeldeep brain stimulation
Abstracts:Deep brain stimulation (DBS) is an established treatment for neurodegenerative movement disorders such as Parkinson's disease that mitigates symptoms by overwriting pathological signals from the central nervous system to the motor system. Nearly all computational models of DBS, directly or indirectly, associate clinical improvements with the extent of fiber activation in the vicinity of the stimulating electrode. However, it is not clear how such activation modulates information transmission. Here, we use the exact cable equation for straight or curved axons and show that DBS segregates the signaling pathways into one of the three communicational modes: complete information blockage, uni-, and bi-directional transmission. Furthermore, all these modes respond to the stimulating pulse in an asynchronous but frequency-locked fashion. Asynchrony depends on the geometry of the axon, its placement and orientation, and the stimulation protocol. At the same time, the electrophysiology of the nerve determines frequency-locking. Such a trimodal response challenges the notion of activation as a binary state and studies that correlate it with the DBS outcome. Importantly, our work suggests that a mechanistic understanding of DBS action relies on distinguishing between these three modes of information transmission.
-
Cluster Neuronal Firing Induced by Uniform Pulses of High-Frequency Stimulation on Axons in Rat Hippocampus
Yue YuanZhouyan FengZhaoxiang Wang
Keywords:NeuronsFiringRecordingElectrodesHafniumAxonsComputational modelingBiomedical engineeringRatsAction potentialsNeuronal FiringHigh-frequency StimulationAction PotentialComputational ModelPyramidal NeuronsNonlinear DynamicsClustering PatternsIndividual NeuronsCA1 RegionFiring PatternsElectrode ArraySpiking Neural NetworksNeuronal SomaHippocampal CA1 RegionStimulation SiteRecording ElectrodesRecording SitesNeuromodulation TechniquesUnit SpikeCell BodiesBurst FiringAction Potential AmplitudeAction Potential InitiationSodium ChannelSpike AmplitudeSingle PulseWeak IntensityBiphasic PulsesFiring RateChanges In AmplitudeHigh-frequency stimulationburstclusterfiring patternunit spikeCA1 pyramidal neuronaxons
Abstracts:Objective: High-frequency stimulation (HFS) of electrical pulse sequences has been used in various neuromodulation techniques to treat certain disorders. Here, we test the hypothesis that HFS sequences with purely periodic pulses could directly generate non-uniform firing in directly stimulated neurons. Methods: In vivo experiments were conducted in the rat hippocampal CA1 region. A stimulation electrode was placed on the alveus fibers, and a recording electrode array was inserted into the CA1 region upstream of the stimulation site. Antidromic-HFS (A-HFS) of 100 Hz pulses was applied to the alveus to antidromically activate the soma of pyramidal neurons around the recording site. By minimizing the interferences of population spikes, the evoked unit spikes of individual pyramidal neurons were obtained during A-HFS. Additionally, a computational model of pyramidal neuron was used to simulate the neuronal responses to A-HFS, revealing possible mechanisms underlying the different firing patterns. Results: Of the total 54 pyramidal neurons recorded during 2-min 100 Hz A-HFS, 38 (70%) neurons fired in a cluster pattern with alternating periods of intensive spikes and silence. The remaining 16 (30%) neurons fired in a non-cluster pattern with regular spikes. Modeling simulations showed that under the situation of HFS-induced intermittent block, conduction failure and generation failure of action potentials along the axons resulted in the cluster and non-cluster firing. Conclusion: Sustained axonal A-HFS with periodic pulses can induce non-uniform firing in directly stimulated neurons. Significance: This finding provides new evidence for the nonlinear dynamics of neuronal firing, even under uniform stimulation.
-
A Bipedal Walking Model Considering Trunk Pitch Angle for Estimating the Influence of Suspension Load on Human Biomechanics
Qinhao ZhangWenbin ChenHanwen ZhangSiyuan LinCaihua Xiong
Keywords:Load modelingLegged locomotionBiological system modelingHipCostsForceBiomechanicsTorqueSuspensions (mechanical systems)DynamicsPitch AngleSuspended LoadBipedal WalkingHuman BiomechanicsTrunk PitchRoot Mean Square ErrorHumoral ResponseLoading ConditionsEffect Of LoadEnergetic CostBiomechanical ResponseExtension TorqueFlexion TorqueTrunk ExtensionLeg StiffnessSimulation ResultsOptimization ProblemOptimal ControlVertical DirectionPhysical WorkTrunk AngleOptimal Control TheoryVertical AccelerationHuman GaitRelative MovementMetabolic CostMechanical PowerLoad ForceMoment Of InertiaGait CycleBiped walking modelsuspended backpacktrunk pitch angletrunk torque
Abstracts:Objective: Suspended loads have been shown to improve loaded-walking economy. Establishing a biped walking model with dynamic trunk pitch angles can provide more comprehensive estimates of the human biomechanical response under suspended loads. Methods: We developed the trunk-load- hip dynamics, modified the spring-loaded-inverted-pendulum (SLIP) model, and optimized the loaded-walking pattern for minimal energetic cost. 9 subjects participated in experiments using a powered backpack to validate the model's performance, conducting two trials: Load-Suspended (LS) and Load-Locked (LL). Results: The averaged correlation coefficient of simulated and experimental hip trajectory, vertical and horizontal GRFs, and individual leg mechanical (ILM) powers are 0.96, 0.97, 0.93, and 0.81, respectively. The RMS error between simulated and experimental peaks of vertical GRFs, braking peaks of horizontal GRFs, and energetic costs was under 10%. Simulation also provides observation on the effect of suspended load on dynamic trunk pitch angles and torques, and leg stiffness. Both the simulation and experiment demonstrated the advantages of LS in reducing GRFs and energetic cost. Additionally, the simulation shows the peaks of trunk flexion and extension torque are reduced by 34.77% (p < 0.05) and 37.88% (p < 0.05) in LS. Conclusion: The model effectively estimates hip trajectory, vertical and horizontal GRFs, ILM powers, and energetic cost, and provides observations on trunk behavior under different load conditions. The model also supports the advantages of suspension load. Significance: Appropriate models could comprehensively reveal the mechanism between the mechanical systems and human biomechanics responses, guide the design of carrying load devices, and provide rapid evaluation of its effects.