Attenuating midline thalamus bursting to mitigate absence epilepsy.

Advancing the mechanistic understanding of absence epilepsy is crucial for developing new therapeutics, especially for patients unresponsive to current treatments. Utilizing a recently developed mouse model of absence epilepsy carrying the BK gain-of-function channelopathy D434G, here we report that attenuating the burst firing of midline thalamus (MLT) neurons effectively prevents absence seizures. We found that enhanced BK channel activity in the BK-D434G MLT neurons promotes synchronized bursting during the ictal phase of absence seizures. Modulating MLT neurons through pharmacological reagents, optogenetic stimulation, or deep brain stimulation effectively attenuates burst firing, leading to reduced absence seizure frequency and increased vigilance. Additionally, enhancing vigilance by amphetamine, a stimulant medication, or physical perturbation also effectively suppresses MLT bursting and prevents absence seizures. These findings suggest that the MLT is a promising target for clinical interventions. Our diverse approaches offer valuable insights for developing next generation therapeutics to treat absence epilepsy.

Computational models of compound nerve action potentials: Efficient filter-based methods to quantify effects of tissue conductivities, conduction distance, and nerve fiber parameters.

BACKGROUND: Peripheral nerve recordings can enhance the efficacy of neurostimulation therapies by providing a feedback signal to adjust stimulation settings for greater efficacy or reduced side effects. Computational models can accelerate the development of interfaces with high signal-to-noise ratio and selective recording. However, validation and tuning of model outputs against in vivo recordings remains computationally prohibitive due to the large number of fibers in a nerve. METHODS: We designed and implemented highly efficient modeling methods for simulating electrically evoked compound nerve action potential (CNAP) signals. The method simulated a subset of fiber diameters present in the nerve using NEURON, interpolated action potential templates across fiber diameters, and filtered the templates with a weighting function derived from fiber-specific conduction velocity and electromagnetic reciprocity outputs of a volume conductor model. We applied the methods to simulate CNAPs from rat cervical vagus nerve. RESULTS: Brute force simulation of a rat vagal CNAP with all 1,759 myelinated and 13,283 unmyelinated fibers in NEURON required 286 and 15,860 CPU hours, respectively, while filtering interpolated templates required 30 and 38 seconds on a desktop computer while maintaining accuracy. Modeled CNAP amplitude could vary by over two orders of magnitude depending on tissue conductivities and cuff opening within experimentally relevant ranges. Conduction distance and fiber diameter distribution also strongly influenced the modeled CNAP amplitude, shape, and latency. Modeled and in vivo signals had comparable shape, amplitude, and latency for myelinated fibers but not for unmyelinated fibers. CONCLUSIONS: Highly efficient methods of modeling neural recordings quantified the large impact that tissue properties, conduction distance, and nerve fiber parameters have on CNAPs. These methods expand the computational accessibility of neural recording models, enable efficient model tuning for validation, and facilitate the design of novel recording interfaces for neurostimulation feedback and understanding physiological systems.

Reverse-engineered models reveal differential membrane properties of autonomic and cutaneous unmyelinated fibers.

Unmyelinated C-fibers constitute the vast majority of axons in peripheral nerves and play key roles in homeostasis and signaling pain. However, little is known about their ion channel expression, which controls their firing properties. Also, because of their small diameters (~ 1 µm), it has not been possible to characterize their membrane properties using voltage clamp. We developed a novel library of isoform-specific ion channel models to serve as the basis functions of our C-fiber models. We then developed a particle swarm optimization (PSO) framework that used the isoform-specific ion channel models to reverse engineer C-fiber membrane properties from measured autonomic and cutaneous C-fiber conduction responses. Our C-fiber models reproduced experimental conduction velocity, chronaxie, action potential duration, intracellular threshold, and paired pulse recovery cycle. The models also matched experimental activity-dependent slowing, a property not included in model optimization. We found that simple conduction responses, characterizing the action potential, were controlled by similar membrane properties in both the autonomic and cutaneous C-fiber models, but complicated conduction response, characterizing the afterpotenials, were controlled by differential membrane properties. The unmyelinated C-fiber models constitute important tools to study autonomic signaling, assess the mechanisms of pain, and design bioelectronic devices. Additionally, the novel reverse engineering approach can be applied to generate models of other neurons where voltage clamp data are not available.

At home adaptive dual target deep brain stimulation in Parkinson's disease with proportional control.

Continuous deep brain stimulation (cDBS) of the subthalamic nucleus (STN) or globus pallidus is an effective treatment for the motor symptoms of Parkinson's disease. The relative benefit of one region over the other is of great interest but cannot usually be compared in the same patient. Simultaneous DBS of both regions may synergistically increase the therapeutic benefit. Continuous DBS is limited by a lack of responsiveness to dynamic, fluctuating symptoms intrinsic to the disease. Adaptive DBS (aDBS) adjusts stimulation in response to biomarkers to improve efficacy, side effects, and efficiency. We combined bilateral DBS of both STN and globus pallidus (dual target DBS) in a prospective within-participant, clinical trial in six patients with Parkinson's disease (n = 6, 55-65 years, n = 2 females). Dual target cDBS was tested for Parkinson's disease symptom control annually over 2 years, measured by motor rating scales, on time without dyskinesia, and medication reduction. Random amplitude experiments probed system dynamics to estimate parameters for aDBS. We then implemented proportional-plus-integral aDBS using a novel distributed (off-implant) architecture. In the home setting, we collected tremor and dyskinesia scores as well as individualized ß and DBS amplitudes. Dual target cDBS reduced motor symptoms as measured by Unified Parkinson's Disease Rating Scale (UPDRS) to a greater degree than either region alone (P < 0.05, linear mixed model) in the cohort. The amplitude of ß-oscillations in the STN correlated to the speed of hand grasp movements for five of six participants (P < 0.05, Pearson correlation). Random amplitude experiments provided insight into temporal windowing to avoid stimulation artefacts and demonstrated a correlation between STN ß amplitude and DBS amplitude. Proportional plus integral control of aDBS reduced average power, while preserving UPDRS III scores in the clinic (P = 0.28, Wilcoxon signed rank), and tremor and dyskinesia scores during blinded testing at home (n = 3, P > 0.05, Wilcoxon ranked sum). In the home setting, DBS power reductions were slight but significant. Dual target cDBS may offer an improvement in treatment of motor symptoms of Parkinson's disease over DBS of either the STN or globus pallidus alone. When combined with proportional plus integral aDBS, stimulation power may be reduced, while preserving the increased benefit of dual target DBS.

The Therapeutic Frequency Profile of Subthalamic Nucleus Deep Brain Stimulation in Rats Is Shaped by Antidromic Spike Failure.

The therapeutic mechanisms of subthalamic nucleus (STN) deep brain stimulation (DBS) may depend on antidromic activation of cortex via the hyperdirect pathway. However, hyperdirect pathway neurons cannot reliably follow high-stimulation frequencies, and the spike failure rate appears to correlate with symptom relief as a function of stimulation frequency. We hypothesized that antidromic spike failure contributes to the cortical desynchronization caused by DBS. We measured in vivo evoked cortical activity in female Sprague Dawley rats and developed a computational model of cortical activation from STN DBS. We modeled stochastic antidromic spike failure to determine how spike failure affected the desynchronization of pathophysiological oscillatory activity in cortex. We found that high-frequency STN DBS desynchronized pathologic oscillations via the masking of intrinsic spiking through a combination of spike collision, refractoriness, and synaptic depletion. Antidromic spike failure shaped the parabolic relationship between DBS frequency and cortical desynchronization, with maximum desynchronization at ∼130 Hz. These findings reveal that antidromic spike failure plays a critical role in mediating the dependency of symptom relief on stimulation frequency.SIGNIFICANCE STATEMENT Deep brain stimulation (DBS) is a highly effective neuromodulation therapy, yet it remains uncertain why conventionally used stimulation frequencies (e.g., ∼130 Hz) are optimal. In this study, we demonstrate a potential explanation for the stimulation frequency dependency of DBS through a combination of in vivo experimental measurements and computational modeling. We show that high-frequency stimulation can desynchronize pathologic firing patterns in populations of neurons by inducing an informational lesion. However, sporadic spike failure at these high frequencies limits the efficacy of the informational lesion, yielding a parabolic profile with optimal effects at ∼130 Hz. This work provides a potential explanation for the therapeutic mechanism of DBS, and highlights the importance of considering spike failure in mechanistic models of DBS.

Relative contributions of different neural sources to the EEG.

Dogma dictates that the EEG signal is generated by postsynaptic currents (PSCs) because there are an enormous number of synapses in the brain, and PSCs have relatively long durations. However, PSCs are not the only potential source of electric fields in the brain. Action potentials, afterpolarizations, and presynaptic activity can also generate electric fields. Experimentally it is exceedingly difficult to delineate the contributions of different sources because they are casually linked. However, using computational modeling, we can interrogate the relative contributions of different neural elements to the EEG. We used a library of neuron models with morphologically realistic axonal arbors to quantify the relative contributions of PSCs, action potentials, and presynaptic activity to the EEG signal. Consistent with prior assertions, PSCs were the largest contributor to the EEG, but action potentials and afterpolarizations can also make appreciable contributions. For a population of neurons generating simultaneous PSCs and action potentials, we found that the action potentials accounted for up to 20% of the source strength while PSCs accounted for the other 80% and presynaptic activity negligibly contributed. Additionally, L5 PCs generated the largest PSC and action potential signals indicating that they the dominant EEG signal generator. Further, action potentials and afterpolarizations were sufficient to generate physiological oscillations, indicating that they are valid source contributors to the EEG. The EEG emerges from a combination of multiple different source, and, while PSCs are the largest contributor, other sources are non-negligible and should be included in modeling, analysis and interpretation of the EEG.

Rapid estimation of cortical neuron activation thresholds by transcranial magnetic stimulation using convolutional neural networks.

BACKGROUND: Transcranial magnetic stimulation (TMS) can modulate neural activity by evoking action potentials in cortical neurons. TMS neural activation can be predicted by coupling subject-specific head models of the TMS-induced electric field (E-field) to populations of biophysically realistic neuron models; however, the significant computational cost associated with these models limits their utility and eventual translation to clinically relevant applications. OBJECTIVE: To develop computationally efficient estimators of the activation thresholds of multi-compartmental cortical neuron models in response to TMS-induced E-field distributions. METHODS: Multi-scale models combining anatomically accurate finite element method (FEM) simulations of the TMS E-field with layer-specific representations of cortical neurons were used to generate a large dataset of activation thresholds. 3D convolutional neural networks (CNNs) were trained on these data to predict thresholds of model neurons given their local E-field distribution. The CNN estimator was compared to an approach using the uniform E-field approximation to estimate thresholds in the non-uniform TMS-induced E-field. RESULTS: The 3D CNNs estimated thresholds with mean absolute percent error (MAPE) on the test dataset below 2.5% and strong correlation between the CNN predicted and actual thresholds for all cell types (R2 > 0.96). The CNNs estimated thresholds with a 2-4 orders of magnitude reduction in the computational cost of the multi-compartmental neuron models. The CNNs were also trained to predict the median threshold of populations of neurons, speeding up computation further. CONCLUSION: 3D CNNs can estimate rapidly and accurately the TMS activation thresholds of biophysically realistic neuron models using sparse samples of the local E-field, enabling simulating responses of large neuron populations or parameter space exploration on a personal computer.

Pathways and parameters of sacral neuromodulation in rats.

The stimulation paradigm for sacral neuromodulation has remained largely unchanged since its inception. We sought to determine, in rats, whether stimulation-induced increases in bladder capacity correlated with the proportion of sensory pudendal (PudS) neurons at each stimulated location (L6, S1). If supported, this finding could guide the choice of stimulation side (left/right) and level (S2, S3, S4) in humans. Unexpectedly, we observed that acute stimulation at clinically relevant (low) amplitudes [1-1.5 × motor threshold (Tm)], did not increase bladder capacity, regardless of stimulus location (L6 or S1). More importantly for the ability to test our hypothesis, there was little anatomic variation, and S1 infrequently contributed nerve fibers to the PudS nerve. During mapping studies we noticed that large increases in PudS nerve activation occurred at amplitudes exceeding 2Tm. Thus, additional cystometric studies were conducted, this time with stimulation of the L6-S1 trunk, to examine further the relationship between stimulation amplitude and cystometric parameters. Stimulation at 1Tm to 6Tm evoked increases in bladder capacity and decreases in voiding efficiency that mirrored those produced by PudS nerve stimulation. Many animal studies involving electrical stimulation of nerves of the lower urinary tract use stimulation amplitudes that exceed those used clinically (∼1Tm). Our results confirm that high amplitudes generate immediate changes in cystometric parameters; however, the relationship to low-amplitude chronic stimulation in humans remains unclear. Additional studies are needed to understand changes that occur with chronic stimulation, how these changes relate to therapeutic outcomes, and the contribution of specific nerve fibers to these changes.NEW & NOTEWORTHY Acute low-amplitude electrical stimulation of sacral nerve (sacral neuromodulation) did not increase bladder capacity in anesthetized CD, obese-prone, or obese-resistant rats. Increasing stimulation amplitude correlated with increases in bladder capacity and pudendal sensory nerve recruitment. It is unclear how the high-amplitude acute stimulation that is commonly used in animal experiments to generate immediate effects compares mechanistically to the chronic low-amplitude stimulation used clinically.

Network model of nociceptive processing in the superficial spinal dorsal horn reveals mechanisms of hyperalgesia, allodynia, and spinal cord stimulation.

The spinal dorsal horn (DH) processes sensory information and plays a key role in transmitting nociception to supraspinal centers. Loss of DH inhibition during neuropathic pain unmasks a pathway from nonnociceptive Aß-afferent inputs to superficial dorsal horn (SDH) nociceptive-specific (NS) projection neurons, and this change may contribute to hyperalgesia and allodynia. We developed and validated a computational model of SDH neuronal circuitry that links nonnociceptive Aß-afferent inputs in lamina II/III to a NS projection neuron in lamina I via a network of excitatory interneurons. The excitatory pathway and the NS projection neuron were in turn gated by inhibitory interneurons with connections based on prior patch-clamp recordings. Changing synaptic weights in the computational model to replicate neuropathic pain states unmasked a low-threshold excitatory pathway to NS neurons similar to experimental recordings. Spinal cord stimulation (SCS) is an effective therapy for neuropathic pain, and accumulating experimental evidence indicates that NS neurons in the SDH also respond to SCS. Accounting for these responses may inform therapeutic improvements, and we quantified responses to SCS in the SDH network model and examined the role of different modes of inhibitory control in modulating NS neuron responses to SCS. We combined the SDH network model with a previously published model of the deep dorsal horn (DDH) and identified optimal stimulation frequencies across different neuropathic pain conditions. Finally, we found that SCS-generated inhibition did not completely suppress model NS activity during simulated pinch inputs, providing an explanation of why SCS does not eliminate acute pain.NEW & NOTEWORTHY Chronic pain is a severe public health problem that reduces the quality of life for those affected and exacts an enormous socio-economic burden worldwide. Spinal cord stimulation (SCS) is an effective treatment for chronic pain, but SCS efficacy has not significantly improved over time, in part because the mechanisms of action remain unclear. Most preclinical studies investigating pain and SCS mechanisms have focused on the responses of deep dorsal horn (DDH) neurons, but neural networks in the superficial dorsal horn (SDH) are also important for processing nociceptive information. This work synthesizes heterogeneous experimental recordings from the SDH into a computational model that replicates experimental responses and that can be used to quantify neuronal responses to SCS under neuropathic pain conditions.

Spatiotemporal parameters for energy efficient kilohertz-frequency nerve block with low onset response.

BACKGROUND: Electrical nerve conduction block has great potential for treatment of disease through reversible and local inactivation of somatic and autonomic nerves. However, the relatively high energy requirements and the presence of undesired excitation at the onset of the kilohertz-frequency (KHF) signals used for block pose obstacles to effective translation. Frequency, electrode geometry, and waveform shape are known to influence block threshold and onset response, but available data provide a limited understanding of how to select these parameters to optimize nerve block. METHODS: We evaluated KHF nerve block in rat tibial nerve across frequencies (5-60 kHz), electrode geometries (monopolar, bipolar, and tripolar), and waveform shapes. We present a novel Fourier-based method for constructing composite signals that systematically sample the KHF waveform design space. RESULTS: The lowest frequencies capable of blocking (5-16 kHz) were not the most energy-efficient among the tested frequencies. Further, bipolar cuffs required the largest current and power to block, monopolar cuffs required the lowest current, and both tripolar and monopolar cuffs required the lowest power. Tripolar cuffs produced the smallest onset response across frequencies. Composite signals comprised of a first harmonic sinusoid at fundamental frequency (f0) superposed on a second harmonic sinusoid at 2f0 could block at lower threshold and lower onset response compared to the constituent sinusoids alone. This effect was strongly dependent on the phase of the second harmonic and on the relative amplitudes of the first and second harmonics. This effect was also dependent on electrode geometry: monopolar and tripolar cuffs showed clear composite signal effects in most experiments; bipolar cuffs showed no clear effects in most experiments. CONCLUSIONS: Our data provide novel information about block threshold and onset response at the boundary of frequencies that can block. Our results also show an interaction between spatial (cuff geometry) and temporal (frequency and waveform shape) parameters. Finally, while previous studies suggested that temporal parameters could reduce onset response only in exchange for increased block threshold (or vice versa), our results show that waveform shape influences KHF response in ways that can be exploited to reduce both energy and onset responses.

The cortical evoked potential corresponds with deep brain stimulation efficacy in rats.

Deep brain stimulation (DBS) of the subthalamic nucleus (STN) antidromically activates the motor cortex (M1), and this cortical activation appears to play a role in the treatment of hypokinetic motor behaviors (Gradinaru V, Mogri M, Thompson KR, Henderson JM, Deisseroth K. Science 324: 354-359, 2009; Yu C, Cassar IR, Sambangi J, Grill WM. J Neurosci 40: 4323-4334, 2020). The synchronous antidromic activation takes the form of a short-latency cortical evoked potential (cEP) in electrocorticography (ECoG) recordings of M1. We assessed the utility of the cEP as a biomarker for STN DBS in unilateral 6-hydroxydopamine-lesioned female Sprague Dawley rats, with stimulating electrodes implanted in the STN and the ECoG recorded above M1. We quantified the correlations of the cEP magnitude and latency with changes in motor behavior from DBS and compared them to the correlation between motor behaviors and several commonly used spectral-based biomarkers. The cEP features correlated strongly with motor behaviors and were highly consistent across animals, whereas the spectral biomarkers correlated weakly with motor behaviors and were highly variable across animals. The cEP may thus be a useful biomarker for assessing the therapeutic efficacy of DBS parameters, as its features strongly correlate with motor behavior, it is consistent across time and subjects, it can be recorded under anesthesia, and it is simple to quantify with a large signal-to-noise ratio, enabling rapid, real-time evaluation. Additionally, our work provides further evidence that antidromic cortical activation mediates changes in motor behavior from STN DBS and that the dependence of DBS efficacy on stimulation frequency may be related to antidromic spike failure.NEW & NOTEWORTHY We characterize a new potential biomarker for deep brain stimulation (DBS), the cortical evoked potential (cEP), and demonstrate that it exhibits a robust correlation with motor behaviors as a function of stimulation frequency. The cEP may thus be a useful clinical biomarker for changes in motor behavior. This work also provides insight into the cortical mechanisms of DBS, suggesting that motor behaviors are strongly affected by the rate of antidromic spike failure during DBS.

ASCENT (Automated Simulations to Characterize Electrical Nerve Thresholds): A pipeline for sample-specific computational modeling of electrical stimulation of peripheral nerves.

Electrical stimulation and block of peripheral nerves hold great promise for treatment of a range of disease and disorders, but promising results from preclinical studies often fail to translate to successful clinical therapies. Differences in neural anatomy across species require different electrodes and stimulation parameters to achieve equivalent nerve responses, and accounting for the consequences of these factors is difficult. We describe the implementation, validation, and application of a standardized, modular, and scalable computational modeling pipeline for biophysical simulations of electrical activation and block of nerve fibers within peripheral nerves. The ASCENT (Automated Simulations to Characterize Electrical Nerve Thresholds) pipeline provides a suite of built-in capabilities for user control over the entire workflow, including libraries for parts to assemble electrodes, electrical properties of biological materials, previously published fiber models, and common stimulation waveforms. We validated the accuracy of ASCENT calculations, verified usability in beta release, and provide several compelling examples of ASCENT-implemented models. ASCENT will enable the reproducibility of simulation data, and it will be used as a component of integrated simulations with other models (e.g., organ system models), to interpret experimental results, and to design experimental and clinical interventions for the advancement of peripheral nerve stimulation therapies.

Activation of ENS Circuits in Mouse Colon: Coordination in the Mouse Colonic Motor Complex as a Robust, Distributed Control System.

The characteristic motor patterns of the colon are coordinated by the enteric nervous system (ENS) and involve enterochromaffin (EC) cells, enteric glia, smooth muscle fibers, and interstitial cells. While the fundamental control mechanisms of colonic motor patterns are understood, greater complexity in the circuitry underlying motor patterns has been revealed by recent advances in the field. We review these recent advances and new findings from our laboratories that provide insights into how the ENS coordinates motor patterns in the isolated mouse colon. We contextualize these observations by describing the neuromuscular system underling the colonic motor complex (CMC) as a robust, distributed control system. Framing the colonic motor complex as a control system reveals a new perspective on the coordinated motor patterns in the colon. We test the control system by applying electrical stimulation in the isolated mouse colon to disrupt the coordination and propagation of the colonic motor complex.

Frequency-Specific Optogenetic Deep Brain Stimulation of Subthalamic Nucleus Improves Parkinsonian Motor Behaviors.

Deep brain stimulation (DBS) of the subthalamic nucleus (STN) is an effective therapy for the motor symptoms of Parkinson's disease (PD). However, the neural elements mediating symptom relief are unclear. A previous study concluded that direct optogenetic activation of STN neurons was neither necessary nor sufficient for relief of parkinsonian symptoms. However, the kinetics of the channelrhodopsin-2 (ChR2) used for cell-specific activation are too slow to follow the high rates required for effective DBS, and thus the contribution of activation of STN neurons to the therapeutic effects of DBS remains unclear. We quantified the behavioral and neuronal effects of optogenetic STN DBS in female rats following unilateral 6-hydroxydopamine (6-OHDA) lesion using an ultrafast opsin (Chronos). Optogenetic STN DBS at 130 pulses per second (pps) reduced pathologic circling and ameliorated deficits in forelimb stepping similarly to electrical DBS, while optogenetic STN DBS with ChR2 did not produce behavioral effects. As with electrical DBS, optogenetic STN DBS exhibited a strong dependence on stimulation rate; high rates produced symptom relief while low rates were ineffective. High-rate optogenetic DBS generated both increases and decreases in firing rates of single neurons in STN, globus pallidus externa (GPe), and substantia nigra pars reticular (SNr), and disrupted ß band oscillatory activity in STN and SNr. High-rate optogenetic STN DBS can indeed ameliorate parkinsonian motor symptoms through reduction of abnormal oscillatory activity in the STN-associated neural circuit, and these results highlight that the kinetic properties of opsins have a strong influence on the effects of optogenetic stimulation.SIGNIFICANCE STATEMENT Whether STN local cells contribute to the therapeutic effects of subthalamic nucleus (STN) deep brain stimulation (DBS) in Parkinson's disease (PD) remains unclear. We re-examined the role of STN local cells in mediating the symptom-relieving effects of STN DBS using cell type-specific optogenetic stimulation with a much faster opsin, Chronos. Direct optogenetic stimulation of STN neurons was effective in treating the symptoms of parkinsonism in the 6-hydroxydopamine (6-OHDA) lesion rat. These results highlight that the kinetic properties of opsins can have a strong influence on the effects of optogenetic activation/inhibition and must be considered when employing optogenetic to study high-rate neural stimulation.

Voiding behavior in awake unrestrained untethered spontaneously hypertensive and Wistar control rats.

The spontaneously hypertensive rat (SHR), a genetic model of high blood pressure, has also been studied as a potential model of overactive bladder. In vivo studies have confirmed the presence of surrogate markers of overactive bladder, including detrusor overactivity, increased urinary frequency, decreased bladder capacity and voided volume (VV), and afferent hypersensitivity to bladder irritation. However, these observations were during awake cystometry using implanted bladder catheters tethered to an infusion pump and artificially filled. We conducted experiments in awake unrestrained untethered age-matched female SHRs and Wistar rats to quantify naïve consumption and voiding behavior and the effect of capsaicin desensitization on consumption and voiding behavior. Food and water consumption, body weight, voiding frequency, and VV were recorded. Rats were placed in metabolism cages for 24 h, up to twice a week, from 17 to 37 wk of age. Compared with Wistar rats, SHRs exhibited decrease in VV and did not exhibit diurnal variation in VV between light and dark periods, suggesting that SHRs may have bladder hypersensitivity. Furthermore, SHRs may also have smaller bladder capacities, as they consumed less water, voided less volume (regardless of light cycle), and had equal urinary frequencies compared with age-matched Wistar rats. We detected no change in SHR voiding behavior following capsaicin desensitization, which was in contrast to a prior awake in vivo cystometry study describing increased VV and micturition interval in SHRs and suggests that C-fiber activity may not contribute to bladder hypersensitivity in SHRs.NEW & NOTEWORTHY We characterized the long-term (20 wk) voiding, defecation, and consumption behavior of age-matched spontaneously hypertensive and Wistar rats without the influence of anesthesia or catheters. Spontaneously hypertensive rats exhibited bladder hypersensitiviy that persisted for the 20-wk duration and was unaffected by capsacin desensitization.

Effects of intravesical prostaglandin E2 on bladder function are preserved in capsaicin-desensitized rats.

Prostaglandin E2 (PGE2) instilled into the bladder generates symptoms of urinary urgency in healthy women and reduces bladder capacity and urethral pressure in both humans and female rats. Systemic capsaicin desensitization, which causes degeneration of C-fibers, prevented PGE2-mediated reductions in bladder capacity, suggesting that PGE2 acts as an irritant (Maggi CA, Giuliani S, Conte B, Furio M, Santicioli P, Meli P, Gragnani L, Meli A. Eur J Pharmacol 145: 105-112, 1988). In the present study, we instilled PGE2 in female rats after capsaicin desensitization but without the hypogastric nerve transection that was conducted in the Maggi et al. study. One week after capsaicin injection (125 mg/kg sc), rats underwent cystometric and urethral perfusion testing under urethane anesthesia with saline and 100 µM PGE2. Similar to naïve rats, capsaicin-desensitized rats exhibited a reduction in bladder capacity from 1.23 ± 0.08 mL to 0.70 ± 0.10 mL (P = 0.002, n = 9), a reduction in urethral perfusion pressure from 19.3 ± 2.1 cmH2O to 10.9 ± 1.2 cmH2O (P = 0.004, n = 9), and a reduction in bladder compliance from 0.13 ± 0.020 mL/cmH2O to 0.090 ± 0.014 mL/cmH2O (P = 0.011, n = 9). Thus, changes in bladder function following the instillation of PGE2 were not dependent on capsaicin-sensitive pathways. Further, these results suggest that urethral relaxation/weakness and/or increased detrusor pressure as a result of decreased compliance may contribute to urinary urgency and highlight potential targets for new therapies for overactive bladder.

Excitation properties of computational models of unmyelinated peripheral axons.

Biophysically based computational models of nerve fibers are important tools for designing electrical stimulation therapies, investigating drugs that affect ion channels, and studying diseases that affect neurons. Although peripheral nerves are primarily composed of unmyelinated axons (i.e., C-fibers), most modeling efforts focused on myelinated axons. We implemented the single-compartment model of vagal afferents from Schild et al. (1994) (Schild JH, Clark JW, Hay M, Mendelowitz D, Andresen MC, Kunze DL. J Neurophysiol 71: 2338-2358, 1994) and extended the model into a multicompartment axon, presenting the first cable model of a C-fiber vagal afferent. We also implemented the updated parameters from the Schild and Kunze (1997) model (Schild JH, Kunze DL. J Neurophysiol 78: 3198-3209, 1997). We compared the responses of these novel models with those of three published models of unmyelinated axons (Rattay F, Aberham M. IEEE Trans Biomed Eng 40: 1201-1209, 1993; Sundt D, Gamper N, Jaffe DB. J Neurophysiol 114: 3140-3153, 2015; Tigerholm J, Petersson ME, Obreja O, Lampert A, Carr R, Schmelz M, Fransén E. J Neurophysiol 111: 1721-1735, 2014) and with experimental data from single-fiber recordings. Comparing the two models by Schild et al. (1994, 1997) revealed that differences in rest potential and action potential shape were driven by changes in maximum conductances rather than changes in sodium channel dynamics. Comparing the five model axons, the conduction speeds and strength-duration responses were largely within expected ranges, but none of the models captured the experimental threshold recovery cycle-including a complete absence of late subnormality in the models-and their action potential shapes varied dramatically. The Tigerholm et al. (2014) model best reproduced the experimental data, but these modeling efforts make clear that additional data are needed to parameterize and validate future models of autonomic C-fibers.NEW & NOTEWORTHY Peripheral nerves are primarily composed of unmyelinated axons, and there is growing interest in electrical stimulation of the autonomic nervous system to treat various diseases. We present the first cable model of an unmyelinated vagal nerve fiber and compare its ion channel isoforms and conduction responses with other published models of unmyelinated axons, establishing important tools for advancing modeling of autonomic nerves.

Control of colonic motility using electrical stimulation to modulate enteric neural activity.

Electrical stimulation of the enteric nervous system (ENS) is an attractive approach to modify gastrointestinal transit. Colonic motor complexes (CMCs) occur with a periodic rhythm, but the ability to elicit a premature CMC depends, at least in part, upon the intrinsic refractory properties of the ENS, which are presently unknown. The objectives of this study were to record myoelectric complexes (MCs, the electrical correlates of CMCs) in the smooth muscle and 1) determine the refractory periods of MCs, 2) inform and evaluate closed-loop stimulation to repetitively evoke MCs, and 3) identify stimulation methods to suppress MC propagation. We dissected the colon from male and female C57BL/6 mice, preserving the integrity of intrinsic circuitry while removing the extrinsic nerves, and measured properties of spontaneous and evoked MCs in vitro. Hexamethonium abolished spontaneous and evoked MCs, confirming the necessary involvement of the ENS for electrically evoked MCs. Electrical stimulation reduced the mean interval between evoked and spontaneous CMCs (24.6 ± 3.5 vs. 70.6 ± 15.7 s, P = 0.0002, n = 7). The absolute refractory period was 4.3 s (95% confidence interval (CI) = 2.8-5.7 s, R2 = 0.7315, n = 8). Electrical stimulation applied during fluid distention-evoked MCs led to an arrest of MC propagation, and following stimulation, MC propagation resumed at an increased velocity (n = 9). The timing parameters of electrical stimulation increased the rate of evoked MCs and the duration of entrainment of MCs, and the refractory period provides insight into timing considerations for designing neuromodulation strategies to treat colonic dysmotility.NEW & NOTEWORTHY Maintained physiological distension of the isolated mouse colon induces rhythmic cyclic myoelectric complexes (MCs). MCs evoked repeatedly by closed-loop electrical stimulation entrain MCs more frequently than spontaneously occurring MCs. Electrical stimulation delivered at the onset of a contraction temporarily suppresses the propagation of MC contractions. Controlled electrical stimulation can either evoke MCs or temporarily delay MCs in the isolated mouse colon, depending on timing relative to ongoing activity.

Kilohertz waveforms optimized to produce closed-state Na+ channel inactivation eliminate onset response in nerve conduction block.

The delivery of kilohertz frequency alternating current (KHFAC) generates rapid, controlled, and reversible conduction block in motor, sensory, and autonomic nerves, but causes transient activation of action potentials at the onset of the blocking current. We implemented a novel engineering optimization approach to design blocking waveforms that eliminated the onset response by moving voltage-gated Na+ channels (VGSCs) to closed-state inactivation (CSI) without first opening. We used computational models and particle swarm optimization (PSO) to design a charge-balanced 10 kHz biphasic current waveform that produced conduction block without onset firing in peripheral axons at specific locations and with specific diameters. The results indicate that it is possible to achieve onset-free KHFAC nerve block by causing CSI of VGSCs. Our novel approach for designing blocking waveforms and the resulting waveform may have utility in clinical applications of conduction block of peripheral nerve hyperactivity, for example in pain and spasticity.

Dysfunctional voiding behavior and impaired muscle contractility in a rat model of detrusor underactivity.

AIMS: Detrusor underactivity (DU) is an understudied health concern with inadequate clinical management. The pathophysiology of DU is unclear, and current therapies fail to improve symptoms. The current studies characterized voiding function and contractility of bladder and urethral tissues in a novel rat model of DU. METHODS: Female obese prone (OP) and obese resistant (OR) rats were fed a 60 kcal% fat diet at 8 weeks old. A subset of rats (n = 4/strain) underwent uroflowmetry biweekly for 18 weeks in metabolic cages. At 40-56 weeks old, rats (n = 9-10/strain) underwent instrumented cystometry under urethane anesthesia. Following cystometry, bladder and urethral tissues (n = 8-9/strain) were harvested for in vitro assessments of contractility in response to carbachol, electric field stimulation, atropine, alpha, beta-methylene ATP, and caffeine. RESULTS: OP rats exhibited increased urinary frequency (p = 0.0031), decreased voided volume (p = 0.0093), and urine flow rate (p = 0.0064) compared to OR rats during uroflowmetry. Bethanechol (10 mg/kg) did not alter uroflowmetry parameters. During cystometry, OP rats exhibited decreased bladder emptying efficiency (p < 0.0001), decreased pressure to generate a void (p < 0.0001), and increased EUS activity during filling (p = 0.0011). Bladder contractility was decreased in OP rats when exposed to carbachol (p < 0.0003) and ATP (p = 0.0004), whereas middle urethral contractility was increased when exposed to carbachol (p = 0.0014), EFS (p = 0.0289), and caffeine (p = 0.0031). CONCLUSION: Impaired cholinergic and purinergic signaling in the bladder may contribute to poor voiding function in OP rats. In addition, increased urethral activity may engage a guarding reflex to augment continence and exacerbate incomplete emptying.