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.

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.

Early Regenerative Capacity in the Porcine Heart.

BACKGROUND: The adult mammalian heart has limited ability to repair itself after injury. Zebrafish, newts, and neonatal mice can regenerate cardiac tissue, largely by cardiac myocyte (CM) proliferation. It is unknown whether hearts of young large mammals can regenerate. METHODS: We examined the regenerative capacity of the pig heart in neonatal animals (ages 2, 3, or 14 days postnatal) after myocardial infarction or sham procedure. Myocardial scar and left ventricular function were determined by cardiac magnetic resonance imaging and echocardiography. Bromodeoxyuridine pulse-chase labeling, histology, immunohistochemistry, and Western blotting were performed to study cell proliferation, sarcomere dynamics, and cytokinesis and to quantify myocardial fibrosis. RNA-sequencing was also performed. RESULTS: After myocardial infarction, there was early and sustained recovery of cardiac function and wall thickness in the absence of fibrosis in 2-day-old pigs. In contrast, older animals developed full-thickness myocardial scarring, thinned walls, and did not recover function. Genome-wide analyses of the infarct zone revealed a strong transcriptional signature of fibrosis in 14-day-old animals that was absent in 2-day-old pigs, which instead had enrichment for cytokinesis genes. In regenerating hearts of the younger animals, up to 10% of CMs in the border zone of the myocardial infarction showed evidence of DNA replication that was associated with markers of myocyte division and sarcomere disassembly. CONCLUSIONS: Hearts of large mammals have regenerative capacity, likely driven by cardiac myocyte division, but this potential is lost immediately after birth.

Durable engraftment of genetically modified FVIII-secreting autologous bone marrow stromal cells in the intramedullary microenvironment.

Genetically modified FVIII-expressing autologous bone marrow-derived mesenchymal stromal cells (BMSCs) could cure haemophilia A. However, culture-expanded BMSCs engraft poorly in extramedullary sites. Here, we compared the intramedullary cavity, skeletal muscle, subcutaneous tissue and systemic circulation as tissue microenvironments that could support durable engraftment of FVIII-secreting BMSC in vivo. A zinc finger nuclease integrated human FVIII transgene into PPP1R12C (intron 1) of culture-expanded primary canine BMSCs. FVIII-secretory capacity of implanted BMSCs in each dog was expressed as an individualized therapy index (number of viable BMSCs implanted × FVIII activity secreted/million BMSCs/24 hours). Plasma samples before and after implantation were assayed for transgenic FVIII protein using an anti-human FVIII antibody having negligible cross-reactivity with canine FVIII. Plasma transgenic FVIII persisted for at least 48 weeks after implantation in the intramedullary cavity. Transgenic FVIII protein levels were low after intramuscular implantation and undetectable after both intravenous infusion and subcutaneous implantation. All plasma samples were negative for anti-human FVIII antibodies. Plasma concentrations and durability of transgenic FVIII secretion showed no correlation with the therapy index. Thus, the implantation site microenvironment is crucial. The intramedullary microenvironment, but not extramedullary tissues, supported durable engraftment of genetically modified autologous FVIII-secreting BMSCs.

Color changing bioadhesive barrier for peripherally inserted central catheters.

Bacteria migration at catheter insertion sites presents a serious complication (bacteraemia) with high mortality rates. One strategy to mediate bacteraemia is a physical barrier at the skin-catheter interface. Herein a colorimetric biosensor adhesive (CathoGlu) is designed and evaluated for both colorimetric detection of bacterial infection and application as a bacteria barrier. The design intent combines viscous, hydrophobic bioadhesive with an organic pH indicator (bromothymol blue). Visual observation can then distinguish healthy skin at pH = ∼5 from an infected catheter insertion site at pH = ∼8. The liquid-to-biorubber transition of CathoGlu formulation occurs via a brief exposure to UVA penlight, providing an elastic barrier to the skin flora. Leachates from CathoGlu demonstrate no genotoxic and skin sensitization effect, assessed by OECD-recommended in vitro and in chemico assays. The CathoGlu formulation was found non-inferior against clinically approved 2-octyl-cyanoacrylate (Dermabond™), and adhesive tape (Micropore™) within an in vivo porcine model. CathoGlu skin adhesive provides new opportunities to prevent sepsis in challenging clinical situations.

Risk of recurrent Helicobacter pylori infection 1 year after initial eradication therapy in 7 Latin American communities.

IMPORTANCE: The long-term effectiveness of Helicobacter pylori eradication programs for preventing gastric cancer will depend on recurrence risk and individual and community factors. OBJECTIVE: To estimate risk of H. pylori recurrence and assess factors associated with successful eradication 1 year after treatment. DESIGN, SETTING, AND PARTICIPANTS: Cohort analysis of 1463 randomized trial participants aged 21 to 65 years from 7 Latin American communities, who were treated for H. pylori and observed between September 2009 and July 2011. INTERVENTIONS: Randomization to 1 of 3 treatment groups: 14-day lansoprazole, amoxicillin, and clarithromycin (triple therapy); 5-day lansoprazole and amoxicillin followed by 5-day lansoprazole, clarithromycin, and metronidazole (sequential); or 5-day lansoprazole, amoxicillin, clarithromycin, and metronidazole (concomitant). Participants with a positive (13)C-urea breath test (UBT) 6 to 8 weeks posttreatment were offered voluntary re-treatment with 14-day bismuth-based quadruple therapy. MEASUREMENTS: Recurrent infection after a negative posttreatment UBT and factors associated with successful eradication at 1-year follow-up. RESULTS: Among participants with UBT-negative results who had a 1-year follow-up UBT (n=1091), 125 tested UBT positive, a recurrence risk of 11.5% (95% CI, 9.6%-13.5%). Recurrence was significantly associated with study site (P = .03), nonadherence to initial therapy (adjusted odds ratio [AOR], 2.94; 95% CI, 1.31-6.13; P = .01), and children in the household (AOR, 1.17; 95% CI, 1.01-1.35 per child; P = .03). Of the 281 with positive posttreatment UBT results, 138 completed re-treatment, of whom 93 tested UBT negative at 1 year. Among the 1340 who had a 1-year UBT, 80.4% (95% CI, 76.4%-83.9%), 79.8% (95% CI, 75.8%-83.5%), and 77.8% (95% CI, 73.6%-81.6%) had UBT-negative results in the triple, sequential, and concomitant groups, respectively (P = .61), with 79.3% overall effectiveness (95% CI, 77.1%-81.5%). In a single-treatment course analysis that ignored the effects of re-treatment, the percentage of UBT-negative results at 1 year was 72.4% (95% CI, 69.9%-74.8%) and was significantly associated with study site (P < .001), adherence to initial therapy (AOR, 0.26; 95% CI, 0.15-0.42; P < .001), male sex (AOR, 1.63; 95% CI, 1.25-2.13; P < .001), and age (AOR, 1.14; 95% CI, 1.02-1.27 per decade; P = .02). One-year effectiveness among all 1463 enrolled participants, considering all missing UBT results as positive, was 72.7% (95% CI, 70.3%-74.9%). CONCLUSIONS AND RELEVANCE: One year after treatment for H. pylori infection, recurrence occurred in 11.5% of participants who had negative posttreatment UBT results. Recurrence determinants (ie, nonadherence and demographics) may be as important as specific antibiotic regimen in determining the long-term success of H. pylori eradication interventions. Study findings are relevant to the feasibility of programs for the primary prevention of gastric cancer in high-incidence regions of Latin America. TRIAL REGISTRATION: clinicaltrials.gov Identifier: NCT01061437.

Improved Methods for Preparing Transverse Sections and Unrolled Whole Mounts of Maize Leaf Primordia for Fluorescence and Confocal Imaging.

In maize (Zea mays) and other grasses (Poaceae), the leaf primordia are deeply ensheathed and rolled within the leaf whorl, making it difficult to study early leaf development. Here, we describe methods for preparing transverse sections and unrolled whole mounts of maize leaf primordia for fluorescence and confocal imaging. The first method uses a wire stripper to remove the upper portions of older leaves, exposing the tip of the leaf primordium and allowing its measurement for more accurate transverse section sampling. The second method uses clear, double-sided nano tape to unroll and mount whole-leaf primordia for imaging. We show the utility of the two methods in visualizing and analyzing fluorescent protein reporters in maize. These methods provide a solution to the challenges presented by the distinctive morphology of maize leaf primordia and will be useful for visualizing and quantifying leaf anatomical and developmental traits in maize and other grass species.

14-day triple, 5-day concomitant, and 10-day sequential therapies for Helicobacter pylori infection in seven Latin American sites: a randomised trial.

BACKGROUND: Evidence from Europe, Asia, and North America suggests that standard three-drug regimens of a proton-pump inhibitor plus amoxicillin and clarithromycin are significantly less effective for eradication of Helicobacter pylori infection than are 5-day concomitant and 10-day sequential four-drug regimens that include a nitroimidazole. These four-drug regimens also entail fewer antibiotic doses than do three-drug regimens and thus could be suitable for eradication programmes in low-resource settings. Few studies in Latin America have been done, where the burden of H pylori-associated diseases is high. We therefore did a randomised trial in Latin America comparing the effectiveness of four-drug regimens given concomitantly or sequentially with that of a standard 14-day regimen of triple therapy. METHODS: Between September, 2009, and June, 2010, we did a randomised trial of empiric 14-day triple, 5-day concomitant, and 10-day sequential therapies for H pylori in seven Latin American sites: Chile, Colombia, Costa Rica, Honduras, Nicaragua, and Mexico (two sites). Participants aged 21-65 years who tested positive for H pylori by a urea breath test were randomly assigned by a central computer using a dynamic balancing procedure to: 14 days of lansoprazole, amoxicillin, and clarithromycin (standard therapy); 5 days of lansoprazole, amoxicillin, clarithromycin, and metronidazole (concomitant therapy); or 5 days of lansoprazole and amoxicillin followed by 5 days of lansoprazole, clarithromycin, and metronidazole (sequential therapy). Eradication was assessed by urea breath test 6-8 weeks after randomisation. The trial was not masked. Our primary outcome was probablity of H pylori eradication. Our analysis was by intention to treat. This trial is registered with ClinicalTrials.gov, registration number NCT01061437. FINDINGS: 1463 participants aged 21-65 years were randomly allocated a treatment: 488 were treated with 14-day standard therapy, 489 with 5-day concomitant therapy, and 486 with 10-day sequential therapy. The probability of eradication with standard therapy was 82·2% (401 of 488), which was 8·6% higher (95% adjusted CI 2·6-14·5) than with concomitant therapy (73·6% [360 of 489]) and 5·6% higher (-0·04% to 11·6) than with sequential therapy (76·5% [372 of 486]). Neither four-drug regimen was significantly better than standard triple therapy in any of the seven sites. INTERPRETATION: Standard 14-day triple-drug therapy is preferable to 5-day concomitant or 10-day sequential four-drug regimens as empiric therapy for H pylori infection in diverse Latin American populations. FUNDING: Bill & Melinda Gates Foundation, US National Institutes of Health.

[Uterine myoma in remnant cervix]. / Miomatosis uterina en cérvix residual.

Parasitic myomas are rare and their ethiopathogenesis is uncertain. They may develop from a detached fibroid adhering to an extrauterine surface in order to obtain new blood supply. It has been stated that they form from uterine or myoma fragments left behind after morcellation in the abdominopelvic cavity and thus are called "iatrogenic". Surgeons must be aware of this recently reported complication related to the increasing number of laparoscopic procedures. Thorough inspection and washing of the abdominal cavity are recommended. A case of a patient with iatrogenic parasitic myomas, which appeared six years after a laparoscopic supracervical hysterectomy involving a morcellator, is reported.

A Porcine Model Using Adipose Stem Cell-Loaded Scaffolds for Alveolar Ridge Augmentation.

Tooth loss greatly affects a person's quality of life and many turn to dental implants to replace lost teeth. The success of a dental implant depends on the amount of alveolar bone supporting the implant, and thus, bone augmentation is often necessary to preserve or build up bone volume in the alveolar ridge. Bone can be augmented with autogenous bone, allografts, or xenografts, but the limitations of such natural bone grafts prompt researchers to develop synthetic scaffolds supplemented with cells and/or bioactive agents as alternative bone grafts. The translation of these combination scaffolds from the laboratory to the clinic requires reliable experimental models that can simulate the clinical conditions in human patients. In this article, we describe the use of a porcine alveolar defect model as a platform to evaluate the efficacy of a novel combination of a three-dimensional-printed polycaprolactone-tricalcium phosphate (PCL-TCP) scaffold and adipose-derived mesenchymal stem cells (AD-MSCs) in lateral alveolar augmentation. The surgical protocol for the defect creation and regenerative surgery, as well as analytical methods to determine the extent of tissue regeneration, are described and discussed.

Establishment of a Captive Cave Nectar Bat (Eonycteris spelaea) Breeding Colony in Singapore.

Bats are known natural reservoirs of several highly pathogenic zoonotic viruses, including Hendra virus, Nipah virus, rabies virus, SARS-like coronaviruses, and suspected ancestral reservoirs of SARS-CoV-2 responsible for the ongoing COVID-19 pandemic. The capacity to survive infections of highly pathogenic agents without severe disease, together with many other unique features, makes bats an ideal animal model for studying the regulation of infection, cancer, and longevity, which is likely to translate into human health outcomes. A key factor that limits bat research is lack of breeding bat colonies. To address this need, a captive bat colony was established in Singapore from 19 wild-caught local cave nectar bats. The bats were screened for specific pathogens before the start of captive breeding. Custom-made cages and an optimized diet inclusive of Wombaroo dietary formula, liquid diet, and supplement of fruits enabled the bats to breed prolifically in our facility. Cages are washed daily and disinfected once every fortnight. Bats are observed daily to detect any sick bat or abnormal behavior. In addition, bats undergo a thorough health check once every 3 to 4 mo to check on their overall wellbeing, perform sampling, and document any potential pregnancy. The current colony houses over 80 bats that are successfully breeding, providing a valuable resource for research in Singapore and overseas.

Non-monotonic kilohertz frequency neural block thresholds arise from amplitude- and frequency-dependent charge imbalance.

Reversible block of nerve conduction using kilohertz frequency electrical signals has substantial potential for treatment of disease. However, the ability to block nerve fibers selectively is limited by poor understanding of the relationship between waveform parameters and the nerve fibers that are blocked. Previous in vivo studies reported non-monotonic relationships between block signal frequency and block threshold, suggesting the potential for fiber-selective block. However, the mechanisms of non-monotonic block thresholds were unclear, and these findings were not replicated in a subsequent in vivo study. We used high-fidelity computational models and in vivo experiments in anesthetized rats to show that non-monotonic threshold-frequency relationships do occur, that they result from amplitude- and frequency-dependent charge imbalances that cause a shift between kilohertz frequency and direct current block regimes, and that these relationships can differ across fiber diameters such that smaller fibers can be blocked at lower thresholds than larger fibers. These results reconcile previous contradictory studies, clarify the mechanisms of interaction between kilohertz frequency and direct current block, and demonstrate the potential for selective block of small fiber diameters.

Individual Magnetoencephalography Response Profiles to Short-Duration L-Dopa in Parkinson's Disease.

Clinical responses to dopamine replacement therapy for individuals with Parkinson's disease (PD) are often difficult to predict. We characterized changes in MDS-UPDRS motor factor scores resulting from a short-duration L-Dopa response (SDR), and investigated how the inter-subject clinical differences could be predicted from motor cortical magnetoencephalography (MEG). MDS-UPDRS motor factor scores and resting-state MEG recordings were collected during SDR from twenty individuals with a PD diagnosis. We used a novel subject-specific strategy based on linear support vector machines to quantify motor cortical oscillatory frequency profiles that best predicted medication state. Motor cortical profiles differed substantially across individuals and showed consistency across multiple data folds. There was a linear relationship between classification accuracy and SDR of lower limb bradykinesia, although this relationship did not persist after multiple comparison correction, suggesting that combinations of spectral power features alone are insufficient to predict clinical state. Factor score analysis of therapeutic response and novel subject-specific machine learning approaches based on subject-specific neuroimaging provide tools to predict outcomes of therapies for PD.

Quantitative comparisons of block thresholds and onset responses for charge-balanced kilohertz frequency waveforms.

OBJECTIVE: There is growing interest in delivering kilohertz frequency (KHF) electrical signals to block conduction in peripheral nerves for treatment of various diseases. Previous studies used different KHF waveforms to achieve block, and it remains unclear how waveform affects nerve block parameters. APPROACH: We quantified the effects of waveform on KHF block of the rat tibial nerve in vivo and in computational models. We compared block thresholds and onset responses across current-controlled sinusoids and charge-balanced rectangular waveforms with different asymmetries and duty cycles. MAIN RESULTS: Sine waves had higher block thresholds than square waves, but used less power at block threshold. Block threshold had an inverse relationship with duty cycle of rectangular waveforms irrespective of waveform asymmetry. Computational model results were consistent with relationships measured in vivo, although the models underestimated the effect of duty cycle on increasing thresholds. The axonal membrane substantially filtered waveforms, the filter transfer function was strikingly similar across waveforms, and filtering resulted in post-filtered rms block thresholds that were approximately constant across waveforms in silico and in vivo. Onset response was not consistently affected by waveform shape, but onset response was smaller at amplitudes well above block threshold. Therefore, waveforms with lower block thresholds (e.g. sine waves or square waves) could be more readily increased to higher amplitudes relative to block threshold to reduce onset response. We also observed a reduction in onset responses across consecutive trials after initial application of supra-block threshold amplitudes. SIGNIFICANCE: Waveform had substantial effects on block thresholds, and the amplitude relative to block threshold had substantial effects on onset response. These data inform choice of waveform in subsequent studies and clinical applications, enhance effective use of block in therapeutic applications, and facilitate the design of parameters that achieve block with minimal onset responses.

Multi-objective particle swarm optimization for postoperative deep brain stimulation targeting of subthalamic nucleus pathways.

OBJECTIVE: The effectiveness of deep brain stimulation (DBS) therapy strongly depends on precise surgical targeting of intracranial leads and on clinical optimization of stimulation settings. Recent advances in surgical targeting, multi-electrode designs, and multi-channel independent current-controlled stimulation are poised to enable finer control in modulating pathways within the brain. However, the large stimulation parameter space enabled by these technologies also poses significant challenges for efficiently identifying the most therapeutic DBS setting for a given patient. Here, we present a computational approach for programming directional DBS leads that is based on a non-convex optimization framework for neural pathway targeting. APPROACH: The algorithm integrates patient-specific pre-operative 7 T MR imaging, post-operative CT scans, and multi-objective particle swarm optimization (MOPSO) methods using dominance based-criteria and incorporating multiple neural pathways simultaneously. The algorithm was evaluated on eight patient-specific models of subthalamic nucleus (STN) DBS to identify electrode configurations and stimulation amplitudes to optimally activate or avoid six clinically relevant pathways: motor territory of STN, non-motor territory of STN, internal capsule, superior cerebellar peduncle, thalamic fasciculus, and hyperdirect pathway. MAIN RESULTS: Across the patient-specific models, single-electrode stimulation showed significant correlations across modeled pathways, particularly for motor and non-motor STN efferents. The MOPSO approach was able to identify multi-electrode configurations that achieved improved targeting of motor STN efferents and hyperdirect pathway afferents than that achieved by any single-electrode monopolar setting at equivalent power levels. SIGNIFICANCE: These results suggest that pathway targeting with patient-specific model-based optimization algorithms can efficiently identify non-trivial electrode configurations for enhancing activation of clinically relevant pathways. However, the results also indicate that inter-pathway correlations can limit selectivity for certain pathways even with directional DBS leads.

Clinical deep brain stimulation strategies for orientation-selective pathway activation.

OBJECTIVE: This study investigated stimulation strategies to increase the selectivity of activating axonal pathways within the brain based on their orientations relative to clinical deep brain stimulation (DBS) lead implants. APPROACH: Previous work has shown how varying electrode shape and controlling the primary electric field direction through preclinical electrode arrays can produce orientation-selective axonal stimulation. Here, we significantly extend those results using computational models to evaluate the degree to which clinical DBS leads can direct stimulus-induced electric fields and generate orientation-selective activation of fiber pathways in the brain. Orientation-selective pulse paradigms were evaluated in conceptual models and in patient-specific models of subthalamic nucleus (STN)-DBS for treating Parkinson's disease. MAIN RESULTS: Single-contact monopolar or two-contact bipolar stimulation through clinical DBS leads with cylindrical electrodes primarily activated axons orientated parallel to the lead. Conversely, multi-contact monopolar stimulation with a cathode-leading pulse waveform selectively activated axons perpendicular to the DBS lead. Clinical DBS leads with segmented rows of electrodes and a single current source provided additional angular resolution for activating axons oriented 0°, ±22.5°, ±45°, ±67.5°, or 90° relative to the lead shaft. Employing multiple independent current sources to deliver unequal amounts of current through these leads further increased the angular resolution of activation relative to the lead shaft. The patient-specific models indicated that multi-contact cathode configurations, which are rarely used in clinical practice, could increase activation of the hyperdirect pathway collaterals projecting into STN (a putative therapeutic target), while minimizing direct activation of the corticospinal tract of internal capsule, which can elicit sensorimotor side-effects when stimulated. SIGNIFICANCE: When combined with patient-specific tissue anisotropy and patient-specific anatomical morphologies of neural pathways responsible for therapy and side effects, orientation-selective DBS approaches show potential to significantly improve clinical outcomes of DBS therapy for a range of existing and investigational clinical indications.

Particle swarm optimization for programming deep brain stimulation arrays.

OBJECTIVE: Deep brain stimulation (DBS) therapy relies on both precise neurosurgical targeting and systematic optimization of stimulation settings to achieve beneficial clinical outcomes. One recent advance to improve targeting is the development of DBS arrays (DBSAs) with electrodes segmented both along and around the DBS lead. However, increasing the number of independent electrodes creates the logistical challenge of optimizing stimulation parameters efficiently. APPROACH: Solving such complex problems with multiple solutions and objectives is well known to occur in biology, in which complex collective behaviors emerge out of swarms of individual organisms engaged in learning through social interactions. Here, we developed a particle swarm optimization (PSO) algorithm to program DBSAs using a swarm of individual particles representing electrode configurations and stimulation amplitudes. Using a finite element model of motor thalamic DBS, we demonstrate how the PSO algorithm can efficiently optimize a multi-objective function that maximizes predictions of axonal activation in regions of interest (ROI, cerebellar-receiving area of motor thalamus), minimizes predictions of axonal activation in regions of avoidance (ROA, somatosensory thalamus), and minimizes power consumption. MAIN RESULTS: The algorithm solved the multi-objective problem by producing a Pareto front. ROI and ROA activation predictions were consistent across swarms (<1% median discrepancy in axon activation). The algorithm was able to accommodate for (1) lead displacement (1 mm) with relatively small ROI (⩽9.2%) and ROA (⩽1%) activation changes, irrespective of shift direction; (2) reduction in maximum per-electrode current (by 50% and 80%) with ROI activation decreasing by 5.6% and 16%, respectively; and (3) disabling electrodes (n = 3 and 12) with ROI activation reduction by 1.8% and 14%, respectively. Additionally, comparison between PSO predictions and multi-compartment axon model simulations showed discrepancies of <1% between approaches. SIGNIFICANCE: The PSO algorithm provides a computationally efficient way to program DBS systems especially those with higher electrode counts.

Determining Electrode Placement for Transcranial Direct Current Stimulation: A Comparison of EEG- Versus TMS-Guided Methods.

Transcranial direct current stimulation (tDCS) is increasingly researched as an adjuvant to motor rehabilitation for children with hemiparesis. The optimal method for the primary motor cortex (M1) somatotopic localization for tDCS electrode placement has not been established. The objective, therefore, was to determine the location of the M1 derived using the 10/20 electroencephalography (EEG) system and transcranial magnetic stimulation (TMS) in children with hemiparesis (CWH) and a comparison group of typically developing children (TDC). We hypothesized a difference in location for CWH but not for TDC. The 2 locations were evaluated in 47 children (21 CWH, 26 TDC). Distances between the locations were measured pending presence of a motor evoked potential. Distances between the EEG and TMS locations that exceeded the 2.5 cm × 2.5 cm rubber electrode area are reported in percentages [95% confidence interval] in CWH-nonlesioned hemisphere was 68.8% [41.3-89.0], lesioned: 85.7% [57.2-98.2]; TDC-dominant hemisphere 73.9% [51.6-89.8], nondominant: 82.6% [61.2-95.0]. Distances that exceeded the 3 × 5 cm electrode sponge area in CWH-nonlesioned was 25.0% [7.3-52.4], lesioned was 28.6% [8.4-58.1]; TDC-dominant was 52.2% [30.6-73.2], nondominant was 43.5 [23.2-65.5]). Distances that exceeded the 5 × 7 cm electrode sponge area in CWH-nonlesioned was 18.8% [4.0-45.6] and lesioned was 21.4% [4.7-50.8]; TDC-dominant was 21.7% [7.5-43.7] and nondominant was 26.1% [10.2-48.4]. Individual variability in brain somatotopic organization may influence surface scalp localization of underlying M1 in children regardless of neurologic impairment. Findings suggest further investigation of optimal tDCS electrode placement. EEG and TMS methods reveal variability in localizing M1 in children regardless of stroke diagnosis. This study was registered on clinicaltrials.gov NCT02015338.

Theoretical Optimization of Stimulation Strategies for a Directionally Segmented Deep Brain Stimulation Electrode Array.

Programming deep brain stimulation (DBS) systems currently involves a clinician manually sweeping through a range of stimulus parameter settings to identify the setting that delivers the most robust therapy for a patient. With the advent of DBS arrays with a higher number and density of electrodes, this trial and error process becomes unmanageable in a clinical setting. This study developed a computationally efficient, model-based algorithm to estimate an electrode configuration that will most strongly activate tissue within a volume of interest. The cerebellar-receiving area of motor thalamus, the target for treating essential tremor with DBS, was rendered from imaging data and discretized into grid points aligned in approximate afferent and efferent axonal pathway orientations. A finite-element model (FEM) was constructed to simulate the volumetric tissue voltage during DBS. We leveraged the principle of voltage superposition to formulate a convex optimization-based approach to maximize activating function (AF) values at each grid point (via three different criteria), hence increasing the overall probability of action potential initiation and neuronal entrainment within the target volume. For both efferent and afferent pathways, this approach achieved global optima within several seconds. The optimal electrode configuration and resulting AF values differed across each optimization criteria and between axonal orientations. This approach only required a set of FEM simulations equal to the number of DBS array electrodes, and could readily accommodate anisotropic-inhomogeneous tissue conductances or other axonal orientations. The algorithm provides an efficient, flexible determination of optimal electrode configurations for programming DBS arrays.