The field of neuromodulation encompasses a wide spectrum of interventional technologies that modify pathological activity within the nervous system to achieve a therapeutic effect. Therapies including deep brain stimulation (DBS), intracranial cortical stimulation (ICS), transcranial direct current stimulation (tDCS), and transcranial magnetic stimulation (TMS) have all shown promising results across a range of neurological and neuropsychiatric disorders. While the mechanisms of therapeutic action are invariably different amongst these approaches, there are several fundamental neuroengineering challenges that are commonly applicable to improving neuromodulation efficacy. This article reviews the state-of-the-art of neuromodulation for brain disorders and discusses the challenges and opportunities available for clinicians and researchers interested in advancing neuromodulation therapies.

Beta frequency oscillations (15 to 35 Hz) in cortical and basal ganglia circuits become abnormally synchronized in Parkinson’s disease (PD). How excessive beta oscillations emerge in these circuits is unclear. We addressed this issue by defining the firing properties of basal ganglia neurons around the emergence of cortical beta bursts (β bursts), transient (50 to 350 ms) increases in the beta amplitude of cortical signals. In PD patients, the phase locking of background spiking activity in the subthalamic nucleus (STN) to frontal electroencephalograms preceded the onset and followed the temporal profile of cortical β bursts, with conditions of synchronization consistent within and across bursts. Neuronal ensemble recordings in multiple basal ganglia structures of parkinsonian rats revealed that these dynamics were recapitulated in STN, but also in external globus pallidus and striatum. The onset of consistent phase-locking conditions was preceded by abrupt phase slips between cortical and basal ganglia ensemble signals. Single-unit recordings demonstrated that ensemble-level properties of synchronization were not underlain by changes in firing rate but, rather, by the timing of action potentials in relation to cortical oscillation phase. Notably, the preferred angle of phase-locked action potential firing in each basal ganglia structure was shifted during burst initiation, then maintained stable phase relations during the burst. Subthalamic, pallidal, and striatal neurons engaged and disengaged with cortical β bursts to different extents and timings. The temporal evolution of cortical and basal ganglia synchronization is cell type-selective, which could be key for the generation/ maintenance of excessive beta oscillations in parkinsonism.

Synchronized oscillations within and between brain areas facilitate normal processing, but are often amplified in disease. A prominent example is the abnormally sustained beta-frequency (∼20 Hz) oscillations recorded from the cortex and subthalamic nucleus of Parkinson's disease patients. Computational modeling suggests that the amplitude of such oscillations could be modulated by applying stimulation at a specific phase. Such a strategy would allow selective targeting of the oscillation, with relatively little effect on other activity parameters. Here, activity was recorded from 10 awake, parkinsonian patients (6 male, 4 female human subjects) undergoing functional neurosurgery. We demonstrate that stimulation arriving on a particular patient-specific phase of the beta oscillation over consecutive cycles could suppress the amplitude of this pathophysiological activity by up to 40%, while amplification effects were relatively weak. Suppressive effects were accompanied by a reduction in the rhythmic output of subthalamic nucleus (STN) neurons and synchronization with the mesial cortex. While stimulation could alter the spiking pattern of STN neurons, there was no net effect on firing rate, suggesting that reduced beta synchrony was a result of alterations to the relative timing of spiking activity, rather than an overall change in excitability. Together, these results identify a novel intrinsic property of cortico-basal ganglia synchrony that suggests the phase of ongoing neural oscillations could be a viable and effective control signal for the treatment of Parkinson's disease. This work has potential implications for other brain diseases with exaggerated neuronal synchronization and for probing the function of rhythmic activity in the healthy brain.SIGNIFICANCE STATEMENT In Parkinson's disease (PD), movement impairment is correlated with exaggerated beta frequency oscillations in the cerebral cortex and subthalamic nucleus (STN). Using a novel method of stimulation in PD patients undergoing neurosurgery, we demonstrate that STN beta oscillations can be suppressed when consecutive electrical pulses arrive at a specific phase of the oscillation. This effect is likely because of interrupting the timing of neuronal activity rather than excitability, as stimulation altered the firing pattern of STN spiking without changing overall rate. These findings show the potential of oscillation phase as an input for “closed-loop” stimulation, which could provide a valuable neuromodulation strategy for the treatment of brain disorders and for elucidating the role of neuronal oscillations in the healthy brain.

We propose a novel, closed-loop approach to tuning deep brain stimulation (DBS) for Parkinson’s disease (PD). The approach, termed Phasic Burst Stimulation (PhaBS), applies a burst of stimulus pulses over a range of phases predicted to disrupt pathological oscillations seen in PD. Stimulation parameters are optimized based on phase response curves (PRCs), which would be measured from each patient. This approach is tested in a computational model of PD with an emergent population oscillation. We show that the stimulus phase can be optimized using the PRC, and that PhaBS is more effective at suppressing the pathological oscillation than a single phasic stimulus pulse. PhaBS provides a closed-loop approach to DBS that can be optimized for each patient.

53BP1 controls a specialized non-homologous end joining (NHEJ) pathway that is essential for adaptive immunity, yet oncogenic in BRCA1 mutant cancers. Intra-chromosomal DNA double-strand break (DSB) joining events during immunoglobulin class switch recombination (CSR) require 53BP1. However, in BRCA1 mutant cells, 53BP1 blocks homologous recombination (HR) and promotes toxic NHEJ, resulting in genomic instability. Here, we identify the protein dimerization hub—DYNLL1—as an organizer of multimeric 53BP1 complexes. DYNLL1 binding stimulates 53BP1 oligomerization, and promotes 53BP1’s recruitment to, and interaction with, DSB-associated chromatin. Consequently, DYNLL1 regulates 53BP1-dependent NHEJ: CSR is compromised upon deletion of Dynll1 or its transcriptional regulator Asciz, or by mutation of DYNLL1 binding motifs in 53BP1; furthermore, Brca1 mutant cells and tumours are rendered resistant to poly-ADP ribose polymerase (PARP) inhibitor treatments upon deletion of Dynll1 or Asciz. Thus, our results reveal a mechanism that regulates 53BP1-dependent NHEJ and the therapeutic response of BRCA1-deficient cancers.

Efficacy of deep brain stimulation (DBS) for motor signs of Parkinson’s disease (PD) depends in part on post-operative programming of stimulus parameters. There is a need for a systematic approach to tuning parameters based on patient physiology. We used a physiologically realistic computational model of the basal ganglia network to investigate the emergence of a 34 Hz oscillation in the PD state and its optimal suppression with DBS. Discrete time transfer functions were fit to post-stimulus time histograms (PSTHs) collected in open-loop, by simulating the pharmacological block of synaptic connections, to describe the behavior of the basal ganglia nuclei. These functions were then connected to create a mean-field model of the closed-loop system, which was analyzed to determine the origin of the emergent 34 Hz pathological oscillation. This analysis determined that the oscillation could emerge from the coupling between the globus pallidus external (GPe) and subthalamic nucleus (STN). When coupled, the two resonate with each other in the PD state but not in the healthy state. By characterizing how this oscillation is affected by subthreshold DBS pulses, we hypothesize that it is possible to predict stimulus frequencies capable of suppressing this oscillation. To characterize the response to the stimulus, we developed a new method for estimating phase response curves (PRCs) from population data. Using the population PRC we were able to predict frequencies that enhance and suppress the 34 Hz pathological oscillation. This provides a systematic approach to tuning DBS frequencies and could enable closed-loop tuning of stimulation parameters.

Computational modeling can be a powerful tool for an experimentalist, providing a rigorous mathematical model of the system you are studying. This can be valuable in testing your hypotheses and developing experimental protocols prior to experimenting. This paper reviews models of seizures and epilepsy at different scales, including cellular, network, cortical region, and brain scales by looking at how they have been used in conjunction with experimental data. At each scale, models with different levels of abstraction, the extraction of physiological detail, are presented. Varying levels of detail are necessary in different situations. Physiologically realistic models are valuable surrogates for experimental systems because, unlike in an experiment, every parameter can be changed and every variable can be observed. Abstract models are useful in determining essential parameters of a system, allowing the experimentalist to extract principles that explain the relationship between mechanisms and the behavior of the system. Modeling is becoming easier with the emergence of platforms dedicated to neuronal modeling and databases of models that can be downloaded. Modeling will never be a replacement for animal and clinical experiments, but it should be a starting point in designing experiments and understanding their results.

We develop a methodology to design a stimulus optimized to entrain nonlinear, noisy limit cycle oscillators with uncertain properties. Conditions are derived which guarantee that the stimulus will entrain the oscillators despite these uncertainties. Using these conditions, we develop an energy optimal control strategy to design an efficient entraining stimulus and apply it to numerical models of noisy phase oscillators and to in vitro hippocampal neurons. In both instances, the optimal stimuli outperform other similar but suboptimal entraining stimuli. Because this control strategy explicitly accounts for both noise and inherent uncertainty of model parameters, it could have experimental relevance to neural circuits where robust spike timing plays an important role.