Evoked resonant neural activity long-term dynamics can be reproduced by a computational model with vesicle depletion.

Subthalamic deep brain stimulation (DBS) robustly generates high-frequency oscillations known as evoked resonant neural activity (ERNA). Recently the importance of ERNA has been demonstrated through its ability to predict the optimal DBS contact in the subthalamic nucleus in patients with Parkinson's disease. However, the underlying mechanisms of ERNA are not well understood, and previous modelling efforts have not managed to reproduce the wealth of published data describing the dynamics of ERNA. Here, we aim to present a minimal model capable of reproducing the characteristics of the slow ERNA dynamics published to date. We make biophysically-motivated modifications to the Kuramoto model and fit its parameters to the slow dynamics of ERNA obtained from data. Our results demonstrate that it is possible to reproduce the slow dynamics of ERNA (over hundreds of seconds) with a single neuronal population, and, crucially, with vesicle depletion as one of the key mechanisms behind the ERNA frequency decay in our model. We further validate the proposed model against experimental data from Parkinson's disease patients, where it captures the variations in ERNA frequency and amplitude in response to variable stimulation frequency, amplitude, and to stimulation pulse bursting. We provide a series of predictions from the model that could be the subject of future studies for further validation.

Evoked resonant neural activity in subthalamic local field potentials reflects basal ganglia network dynamics.

Evoked resonant neural activity (ERNA) is induced by subthalamic deep brain stimulation (DBS) and was recently suggested as a marker of lead placement and contact selection in Parkinson's disease. Yet, its underlying mechanisms and how it is modulated by stimulation parameters are unclear. Here, we recorded local field potentials from 27 Parkinson's disease patients, while leads were externalised to scrutinise the ERNA. First, we show that ERNA in the time series waveform and spectrogram likely represent the same activity, which was contested before. Second, our results show that the ERNA has fast and slow dynamics during stimulation, consistent with the synaptic failure hypothesis. Third, we show that ERNA parameters are modulated by different DBS frequencies, intensities, medication states and stimulation modes (continuous DBS vs. adaptive DBS). These results suggest the ERNA might prove useful as a predictor of the best DBS frequency and lowest effective intensity in addition to contact selection. Changes with levodopa and DBS mode suggest that the ERNA may indicate the state of the cortico-basal ganglia circuit making it a putative biomarker to track clinical state in adaptive DBS.

Spectral Topography of the Subthalamic Nucleus to Inform Next-Generation Deep Brain Stimulation.

BACKGROUND: The landscape of neurophysiological symptoms and behavioral biomarkers in basal ganglia signals for movement disorders is expanding. The clinical translation of sensing-based deep brain stimulation (DBS) also requires a thorough understanding of the anatomical organization of spectral biomarkers within the subthalamic nucleus (STN). OBJECTIVES: The aims were to systematically investigate the spectral topography, including a wide range of sub-bands in STN local field potentials (LFP) of Parkinson's disease (PD) patients, and to evaluate its predictive performance for clinical response to DBS. METHODS: STN-LFPs were recorded from 70 PD patients (130 hemispheres) awake and at rest using multicontact DBS electrodes. A comprehensive spatial characterization, including hot spot localization and focality estimation, was performed for multiple sub-bands (delta, theta, alpha, low-beta, high-beta, low-gamma, high-gamma, and fast-gamma (FG) as well as low- and fast high-frequency oscillations [HFO]) and compared to the clinical hot spot for rigidity response to DBS. A spectral biomarker map was established and used to predict the clinical response to DBS. RESULTS: The STN shows a heterogeneous topographic distribution of different spectral biomarkers, with the strongest segregation in the inferior-superior axis. Relative to the superiorly localized beta hot spot, HFOs (FG, slow HFO) were localized up to 2 mm more inferiorly. Beta oscillations are spatially more spread compared to other sub-bands. Both the spatial proximity of contacts to the beta hot spot and the distance to higher-frequency hot spots were predictive for the best rigidity response to DBS. CONCLUSIONS: The spatial segregation and properties of spectral biomarkers within the DBS target structure can additionally be informative for the implementation of next-generation sensing-based DBS. © 2023 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

Mean-Field Approximations With Adaptive Coupling for Networks With Spike-Timing-Dependent Plasticity.

Understanding the effect of spike-timing-dependent plasticity (STDP) is key to elucidating how neural networks change over long timescales and to design interventions aimed at modulating such networks in neurological disorders. However, progress is restricted by the significant computational cost associated with simulating neural network models with STDP and by the lack of low-dimensional description that could provide analytical insights. Phase-difference-dependent plasticity (PDDP) rules approximate STDP in phase oscillator networks, which prescribe synaptic changes based on phase differences of neuron pairs rather than differences in spike timing. Here we construct mean-field approximations for phase oscillator networks with STDP to describe part of the phase space for this very high-dimensional system. We first show that single-harmonic PDDP rules can approximate a simple form of symmetric STDP, while multiharmonic rules are required to accurately approximate causal STDP. We then derive exact expressions for the evolution of the average PDDP coupling weight in terms of network synchrony. For adaptive networks of Kuramoto oscillators that form clusters, we formulate a family of low-dimensional descriptions based on the mean-field dynamics of each cluster and average coupling weights between and within clusters. Finally, we show that such a two-cluster mean-field model can be fitted to synthetic data to provide a low-dimensional approximation of a full adaptive network with symmetric STDP. Our framework represents a step toward a low-dimensional description of adaptive networks with STDP, and could for example inform the development of new therapies aimed at maximizing the long-lasting effects of brain stimulation.

Optimal closed-loop deep brain stimulation using multiple independently controlled contacts.

Deep brain stimulation (DBS) is a well-established treatment option for a variety of neurological disorders, including Parkinson's disease and essential tremor. The symptoms of these disorders are known to be associated with pathological synchronous neural activity in the basal ganglia and thalamus. It is hypothesised that DBS acts to desynchronise this activity, leading to an overall reduction in symptoms. Electrodes with multiple independently controllable contacts are a recent development in DBS technology which have the potential to target one or more pathological regions with greater precision, reducing side effects and potentially increasing both the efficacy and efficiency of the treatment. The increased complexity of these systems, however, motivates the need to understand the effects of DBS when applied to multiple regions or neural populations within the brain. On the basis of a theoretical model, our paper addresses the question of how to best apply DBS to multiple neural populations to maximally desynchronise brain activity. Central to this are analytical expressions, which we derive, that predict how the symptom severity should change when stimulation is applied. Using these expressions, we construct a closed-loop DBS strategy describing how stimulation should be delivered to individual contacts using the phases and amplitudes of feedback signals. We simulate our method and compare it against two others found in the literature: coordinated reset and phase-locked stimulation. We also investigate the conditions for which our strategy is expected to yield the most benefit.

Average beta burst duration profiles provide a signature of dynamical changes between the ON and OFF medication states in Parkinson's disease.

Parkinson's disease motor symptoms are associated with an increase in subthalamic nucleus beta band oscillatory power. However, these oscillations are phasic, and there is a growing body of evidence suggesting that beta burst duration may be of critical importance to motor symptoms. This makes insights into the dynamics of beta bursting generation valuable, in particular to refine closed-loop deep brain stimulation in Parkinson's disease. In this study, we ask the question "Can average burst duration reveal how dynamics change between the ON and OFF medication states?". Our analysis of local field potentials from the subthalamic nucleus demonstrates using linear surrogates that the system generating beta oscillations is more likely to act in a non-linear regime OFF medication and that the change in a non-linearity measure is correlated with motor impairment. In addition, we pinpoint the simplest dynamical changes that could be responsible for changes in the temporal patterning of beta oscillations between medication states by fitting to data biologically inspired models, and simpler beta envelope models. Finally, we show that the non-linearity can be directly extracted from average burst duration profiles under the assumption of constant noise in envelope models. This reveals that average burst duration profiles provide a window into burst dynamics, which may underlie the success of burst duration as a biomarker. In summary, we demonstrate a relationship between average burst duration profiles, dynamics of the system generating beta oscillations, and motor impairment, which puts us in a better position to understand the pathology and improve therapies such as deep brain stimulation.

Predicting the effects of deep brain stimulation using a reduced coupled oscillator model.

Deep brain stimulation (DBS) is known to be an effective treatment for a variety of neurological disorders, including Parkinson's disease and essential tremor (ET). At present, it involves administering a train of pulses with constant frequency via electrodes implanted into the brain. New 'closed-loop' approaches involve delivering stimulation according to the ongoing symptoms or brain activity and have the potential to provide improvements in terms of efficiency, efficacy and reduction of side effects. The success of closed-loop DBS depends on being able to devise a stimulation strategy that minimizes oscillations in neural activity associated with symptoms of motor disorders. A useful stepping stone towards this is to construct a mathematical model, which can describe how the brain oscillations should change when stimulation is applied at a particular state of the system. Our work focuses on the use of coupled oscillators to represent neurons in areas generating pathological oscillations. Using a reduced form of the Kuramoto model, we analyse how a patient should respond to stimulation when neural oscillations have a given phase and amplitude, provided a number of conditions are satisfied. For such patients, we predict that the best stimulation strategy should be phase specific but also that stimulation should have a greater effect if applied when the amplitude of brain oscillations is lower. We compare this surprising prediction with data obtained from ET patients. In light of our predictions, we also propose a new hybrid strategy which effectively combines two of the closed-loop strategies found in the literature, namely phase-locked and adaptive DBS.

How to design optimal brain stimulation to modulate phase-amplitude coupling?

Objective.Phase-amplitude coupling (PAC), the coupling of the amplitude of a faster brain rhythm to the phase of a slower brain rhythm, plays a significant role in brain activity and has been implicated in various neurological disorders. For example, in Parkinson's disease, PAC between the beta (13-30 Hz) and gamma (30-100 Hz) rhythms in the motor cortex is exaggerated, while in Alzheimer's disease, PAC between the theta (4-8 Hz) and gamma rhythms is diminished. Modulating PAC (i.e. reducing or enhancing PAC) using brain stimulation could therefore open new therapeutic avenues. However, while it has been previously reported that phase-locked stimulation can increase PAC, it is unclear what the optimal stimulation strategy to modulate PAC might be. Here, we provide a theoretical framework to narrow down the experimental optimisation of stimulation aimed at modulating PAC, which would otherwise rely on trial and error.Approach.We make analytical predictions using a Stuart-Landau model, and confirm these predictions in a more realistic model of coupled neural populations.Main results.Our framework specifies the critical Fourier coefficients of the stimulation waveform which should be tuned to optimally modulate PAC. Depending on the characteristics of the amplitude response curve of the fast population, these components may include the slow frequency, the fast frequency, combinations of these, as well as their harmonics. We also show that the optimal balance of energy between these Fourier components depends on the relative strength of the endogenous slow and fast rhythms, and that the alignment of fast components with the fast rhythm should change throughout the slow cycle. Furthermore, we identify the conditions requiring to phase-lock stimulation to the fast and/or slow rhythms.Significance.Together, our theoretical framework lays the foundation for guiding the development of innovative and more effective brain stimulation aimed at modulating PAC for therapeutic benefit.

How to entrain a selected neuronal rhythm but not others: open-loop dithered brain stimulation for selective entrainment.

Objective.While brain stimulation therapies such as deep brain stimulation for Parkinson's disease (PD) can be effective, they have yet to reach their full potential across neurological disorders. Entraining neuronal rhythms using rhythmic brain stimulation has been suggested as a new therapeutic mechanism to restore neurotypical behaviour in conditions such as chronic pain, depression, and Alzheimer's disease. However, theoretical and experimental evidence indicate that brain stimulation can also entrain neuronal rhythms at sub- and super-harmonics, far from the stimulation frequency. Crucially, these counterintuitive effects could be harmful to patients, for example by triggering debilitating involuntary movements in PD. We therefore seek a principled approach to selectively promote rhythms close to the stimulation frequency, while avoiding potential harmful effects by preventing entrainment at sub- and super-harmonics.Approach.Our open-loop approach to selective entrainment, dithered stimulation, consists in adding white noise to the stimulation period.Main results.We theoretically establish the ability of dithered stimulation to selectively entrain a given brain rhythm, and verify its efficacy in simulations of uncoupled neural oscillators, and networks of coupled neural oscillators. Furthermore, we show that dithered stimulation can be implemented in neurostimulators with limited capabilities by toggling within a finite set of stimulation frequencies.Significance.Likely implementable across a variety of existing brain stimulation devices, dithering-based selective entrainment has potential to enable new brain stimulation therapies, as well as new neuroscientific research exploiting its ability to modulate higher-order entrainment.

When do bursts matter in the primary motor cortex? Investigating changes in the intermittencies of beta rhythms associated with movement states.

Brain activity exhibits significant temporal structure that is not well captured in the power spectrum. Recently, attention has shifted to characterising the properties of intermittencies in rhythmic neural activity (i.e. bursts), yet the mechanisms that regulate them are unknown. Here, we present evidence from electrocorticography recordings made over the motor cortex to show that the statistics of bursts, such as duration or amplitude, in the beta frequency (14-30 Hz) band, significantly aid the classification of motor states such as rest, movement preparation, execution, and imagery. These features reflect nonlinearities not detectable in the power spectrum, with states increasing in nonlinearity from movement execution to preparation to rest. Further, we show using a computational model of the cortical microcircuit, constrained to account for burst features, that modulations of laminar specific inhibitory interneurons are responsible for the temporal organisation of activity. Finally, we show that the temporal characteristics of spontaneous activity can be used to infer the balance of cortical integration between incoming sensory information and endogenous activity. Critically, we contribute to the understanding of how transient brain rhythms may underwrite cortical processing, which in turn, could inform novel approaches for brain state classification, and modulation with novel brain-computer interfaces.

Sub-harmonic entrainment of cortical gamma oscillations to deep brain stimulation in Parkinson's disease: Model based predictions and validation in three human subjects.

OBJECTIVES: The exact mechanisms of deep brain stimulation (DBS) are still an active area of investigation, in spite of its clinical successes. This is due in part to the lack of understanding of the effects of stimulation on neuronal rhythms. Entrainment of brain oscillations has been hypothesised as a potential mechanism of neuromodulation. A better understanding of entrainment might further inform existing methods of continuous DBS, and help refine algorithms for adaptive methods. The purpose of this study is to develop and test a theoretical framework to predict entrainment of cortical rhythms to DBS across a wide range of stimulation parameters. MATERIALS AND METHODS: We fit a model of interacting neural populations to selected features characterising PD patients' off-stimulation finely-tuned gamma rhythm recorded through electrocorticography. Using the fitted models, we predict basal ganglia DBS parameters that would result in 1:2 entrainment, a special case of sub-harmonic entrainment observed in patients and predicted by theory. RESULTS: We show that the neural circuit models fitted to patient data exhibit 1:2 entrainment when stimulation is provided across a range of stimulation parameters. Furthermore, we verify key features of the region of 1:2 entrainment in the stimulation frequency/amplitude space with follow-up recordings from the same patients, such as the loss of 1:2 entrainment above certain stimulation amplitudes. CONCLUSION: Our results reveal that continuous, constant frequency DBS in patients may lead to nonlinear patterns of neuronal entrainment across stimulation parameters, and that these responses can be predicted by modelling. Should entrainment prove to be an important mechanism of therapeutic stimulation, our modelling framework may reduce the parameter space that clinicians must consider when programming devices for optimal benefit.

Optimizing deep brain stimulation based on isostable amplitude in essential tremor patient models.

OBJECTIVE: Deep brain stimulation is a treatment for medically refractory essential tremor. To improve the therapy, closed-loop approaches are designed to deliver stimulation according to the system's state, which is constantly monitored by recording a pathological signal associated with symptoms (e.g. brain signal or limb tremor). Since the space of possible closed-loop stimulation strategies is vast and cannot be fully explored experimentally, how to stimulate according to the state should be informed by modeling. A typical modeling goal is to design a stimulation strategy that aims to maximally reduce the Hilbert amplitude of the pathological signal in order to minimize symptoms. Isostables provide a notion of amplitude related to convergence time to the attractor, which can be beneficial in model-based control problems. However, how isostable and Hilbert amplitudes compare when optimizing the amplitude response to stimulation in models constrained by data is unknown. APPROACH: We formulate a simple closed-loop stimulation strategy based on models previously fitted to phase-locked deep brain stimulation data from essential tremor patients. We compare the performance of this strategy in suppressing oscillatory power when based on Hilbert amplitude and when based on isostable amplitude. We also compare performance to phase-locked stimulation and open-loop high-frequency stimulation. MAIN RESULTS: For our closed-loop phase space stimulation strategy, stimulation based on isostable amplitude is significantly more effective than stimulation based on Hilbert amplitude when amplitude field computation time is limited to minutes. Performance is similar when there are no constraints, however constraints on computation time are expected in clinical applications. Even when computation time is limited to minutes, closed-loop phase space stimulation based on isostable amplitude is advantageous compared to phase-locked stimulation, and is more efficient than high-frequency stimulation. SIGNIFICANCE: Our results suggest a potential benefit to using isostable amplitude more broadly for model-based optimization of stimulation in neurological disorders.

Case Report: Embedding "Digital Chronotherapy" Into Medical Devices-A Canine Validation for Controlling Status Epilepticus Through Multi-Scale Rhythmic Brain Stimulation.

Circadian and other physiological rhythms play a key role in both normal homeostasis and disease processes. Such is the case of circadian and infradian seizure patterns observed in epilepsy. However, these rhythms are not fully exploited in the design of active implantable medical devices. In this paper we explore a new implantable stimulator that implements chronotherapy as a feedforward input to supplement both open-loop and closed-loop methods. This integrated algorithm allows for stimulation to be adjusted to the ultradian, circadian and infradian patterns observed in patients through slowly-varying temporal adjustments of stimulation and algorithm sub-components, while also enabling adaption of stimulation based on immediate physiological needs such as a breakthrough seizure or change of posture. Embedded physiological sensors in the stimulator can be used to refine the baseline stimulation circadian pattern as a "digital zeitgeber," i.e., a source of stimulus that entrains or synchronizes the subject's natural rhythms. This algorithmic approach is tested on a canine with severe drug-resistant idiopathic generalized epilepsy exhibiting a characteristic diurnal pattern correlated with sleep-wake cycles. Prior to implantation, the canine's cluster seizures evolved to status epilepticus (SE) and required emergency pharmacological intervention. The cranially-mounted system was fully-implanted bilaterally into the centromedian nucleus of the thalamus. Using combinations of time-based modulation, thalamocortical rhythm-specific tuning of frequency parameters as well as fast-adaptive modes based on activity, the canine experienced no further SE events post-implant as of the time of writing (7 months). Importantly, no significant cluster seizures have been observed either, allowing the reduction of rescue medication. The use of digitally-enabled chronotherapy as a feedforward signal to augment adaptive neurostimulators could prove a useful algorithmic method in conditions where sensitivity to temporal patterns are characteristics of the disease state, providing a novel mechanism for tailoring a more patient-specific therapy approach.

Phase-dependence of response curves to deep brain stimulation and their relationship: from essential tremor patient data to a Wilson-Cowan model.

Essential tremor manifests predominantly as a tremor of the upper limbs. One therapy option is high-frequency deep brain stimulation, which continuously delivers electrical stimulation to the ventral intermediate nucleus of the thalamus at about 130 Hz. Constant stimulation can lead to side effects, it is therefore desirable to find ways to stimulate less while maintaining clinical efficacy. One strategy, phase-locked deep brain stimulation, consists of stimulating according to the phase of the tremor. To advance methods to optimise deep brain stimulation while providing insights into tremor circuits, we ask the question: can the effects of phase-locked stimulation be accounted for by a canonical Wilson-Cowan model? We first analyse patient data, and identify in half of the datasets significant dependence of the effects of stimulation on the phase at which stimulation is provided. The full nonlinear Wilson-Cowan model is fitted to datasets identified as statistically significant, and we show that in each case the model can fit to the dynamics of patient tremor as well as to the phase response curve. The vast majority of top fits are stable foci. The model provides satisfactory prediction of how patient tremor will react to phase-locked stimulation by predicting patient amplitude response curves although they were not explicitly fitted. We also approximate response curves of the significant datasets by providing analytical results for the linearisation of a stable focus model, a simplification of the Wilson-Cowan model in the stable focus regime. We report that the nonlinear Wilson-Cowan model is able to describe response to stimulation more precisely than the linearisation.

Controlled Sub-Nanometer Epitope Spacing in a Three-Dimensional Self-Assembled Peptide Hydrogel.

Cells in the body use a variety of mechanisms to ensure the specificity and efficacy of signal transduction. One way that this is achieved is through tight spatial control over the position of different proteins, signaling sequences, and biomolecules within and around cells. For instance, the extracellular matrix protein fibronectin presents RGDS and PHSRN sequences that synergistically bind the α5ß1 integrin when separated by 3.2 nm but are unable to bind when this distance is >5.5 nm.1 Building biomaterials to controllably space different epitopes with subnanometer accuracy in a three-dimensional (3D) hydrogel is challenging. Here, we synthesized peptides that self-assemble into nanofiber hydrogels utilizing the ß-sheet motif, which has a known regular spacing along the peptide backbone. By modifying specific locations along the peptide, we are able to controllably space different epitopes with subnanometer accuracy at distances from 0.7 nm to over 6 nm, which is within the size range of many protein clusters. Endothelial cells encapsulated within hydrogels displaying RGDS and PHSRN in the native 3.2 nm spacing showed a significant upregulation in the expression of the alpha 5 integrin subunit compared to those in hydrogels with a 6.2 nm spacing, demonstrating the physiological relevance of the spacing. Furthermore, after 24 h the cells in hydrogels with the 3.2 nm spacing appeared to be more spread with increased staining for the α5ß1 integrin. This self-assembling peptide system can controllably space multiple epitopes with subnanometer accuracy, demonstrating an exciting platform to study the effects of ligand density and location on cells within a synthetic 3D environment.