Machine learning to optimize spinal cord stimulation

An example of a system may include a processor and a memory device comprising instructions, which when executed by the processor, cause the processor to: access a patient metric of a subject; use the patient metric as an input to a machine learning algorithm, the machine learning algorithm to search a plurality of neuromodulation parameter sets and to identify a candidate neuromodulation parameter set of the plurality of neuromodulation parameter sets, the candidate neuromodulation parameter set designed to produce a non-regular waveform that varies over a time domain and a space domain; and program a neuromodulator using the candidate neuromodulation parameter set to stimulate the subject.

TECHNICAL FIELD

This document relates generally to medical devices, and more particularly, to systems, devices, and methods for delivering neural modulation.

BACKGROUND

Neuromodulation, which includes neurostimulation, has been proposed as a therapy for a number of conditions. Examples of neurostimulation include Spinal Cord Stimulation (SCS), Deep Brain Stimulation (DBS), Peripheral Nerve Stimulation (PNS), and Functional Electrical Stimulation (FES). Implantable neurostimulation systems have been applied to deliver such a therapy. An implantable neurostimulation system may include an implantable neurostimulator, also referred to as an implantable pulse generator (IPG), and one or more implantable leads each including one or more electrodes. The implantable neurostimulator delivers neurostimulation energy through one or more electrodes placed on or near a target site in the nervous system. An external programming device is used to program the implantable neurostimulator with stimulation parameters controlling the delivery of the neurostimulation energy.

The neurostimulation energy may be delivered in the form of electrical neurostimulation pulses. The delivery is controlled using stimulation parameters that specify spatial (where to stimulate), temporal (when to stimulate), and informational (patterns of pulses directing the nervous system to respond as desired) aspects of a pattern of neurostimulation pulses. Many current neurostimulation systems are programmed to deliver periodic pulses with one or a few uniform waveforms continuously or in bursts. However, neural signals may include more sophisticated patterns to communicate various types of information, including sensations of pain, pressure, temperature, etc. The nervous system may interpret an artificial stimulation with a simple pattern of stimuli as an unnatural phenomenon, and respond with an unintended and undesirable sensation and/or movement. For example, some neurostimulation therapies are known to cause paresthesia and/or vibration of non-targeted tissue or organ.

Recent research has shown that the efficacy and efficiency of certain neurostimulation therapies can be improved, and their side-effects can be reduced, by using patterns of neurostimulation pulses that emulate natural patterns of neural signals observed in the human body.

SUMMARY

Example 1 includes subject matter (such as a method, means for performing acts, machine readable medium including instructions that when performed by a machine cause the machine to performs acts, or an apparatus to perform) comprising: accessing, at a computerized system, a patient metric of a subject; using the patient metric as an input to a machine learning algorithm executing on the computerized system, the machine learning algorithm to search a plurality of neuromodulation parameter sets and to identify a candidate neuromodulation parameter set of the plurality of neuromodulation parameter sets, the candidate neuromodulation parameter set designed to produce a non-regular waveform that varies over a time domain and a space domain; and programming a neuromodulator using the candidate neuromodulation parameter set to stimulate the subject.

In Example 2, the subject matter of Example 1 may include, wherein the subject is a patient.

In Example 3, the subject matter of any one of Examples 1 to 2 may include, wherein the subject is an animal from a preclinical trial.

In Example 4, the subject matter of any one of Examples 1 to 3 may include, wherein the machine learning algorithm comprises a genetic algorithm.

In Example 5, the subject matter of any one of Examples 1 to 4 may include, wherein the machine learning algorithm comprises a neural network.

In Example 6, the subject matter of any one of Examples 1 to 5 may include, wherein the computerized system is a cloud-based system, and wherein programming the neuromodulator comprises transmitting the candidate neuromodulation parameter set to a client device of the cloud-based system.

In Example 7, the subject matter of any one of Examples 1 to 6 may include, wherein the patient metric is an objective pain metric.

In Example 8, the subject matter of any one of Examples 1 to 7 may include, wherein the objective pain metric is a physiological indication sensed by a sensor worn by the subject.

In Example 9, the subject matter of any one of Examples 1 to 8 may include, wherein the patient metric is a subjective pain metric.

In Example 10, the subject matter of any one of Examples 1 to 9 may include, wherein the subjective pain metric is obtained from querying the subject.

In Example 11, the subject matter of any one of Examples 1 to 10 may include, wherein the space domain includes at least one parameter related to fractionalization or polarity.

In Example 12, the subject matter of any one of Examples 1 to 11 may include, wherein the candidate neuromodulation parameter set is designed to produce a pulse burst that varies pulse-by-pulse.

In Example 13, the subject matter of any one of Examples 1 to 12 may include, receiving from a user, a selected parameter; and focusing the machine learning algorithm using the selected parameter.

Example 14 includes a machine-readable medium including instructions, which when executed by a machine, cause the machine to perform operations of any of the Examples 1-13.

Example 15 includes an apparatus comprising means for performing any of the Examples 1-13.

Example 16 includes subject matter (such as a device, apparatus, or machine) comprising: a processor; and a memory device comprising instructions, which when executed by the processor, cause the processor to: access a patient metric of a subject; use the patient metric as an input to a machine learning algorithm, the machine learning algorithm to search a plurality of neuromodulation parameter sets and identify a candidate neuromodulation parameter set of the plurality of neuromodulation parameter sets, the candidate neuromodulation parameter set designed to produce a non-regular waveform that varies over a time domain and a space domain; and program a neuromodulator using the candidate neuromodulation parameter set to stimulate the subject.

In Example 17, the subject matter of Example 16 may include, wherein the subject is a patient.

In Example 18, the subject matter of any one of Examples 16 to 17 may include, wherein the subject is an animal from a preclinical trial.

In Example 19, the subject matter of any one of Examples 16 to 18 may include, wherein the machine learning algorithm comprises a genetic algorithm.

In Example 20, the subject matter of any one of Examples 16 to 19 may include, wherein the machine learning algorithm comprises a neural network.

In Example 21, the subject matter of any one of Examples 16 to 20 may include, wherein the system is a cloud-based system, and wherein the instructions to provide the candidate neuromodulation parameter set comprise instructions to transmit the candidate neuromodulation parameter set to a client device of the cloud-based system.

In Example 22, the subject matter of any one of Examples 16 to 21 may include, wherein the patient metric is an objective pain metric.

In Example 23, the subject matter of any one of Examples 16 to 22 may include, wherein the objective pain metric is a physiological indication sensed by a sensor worn by the subject.

In Example 24, the subject matter of any one of Examples 16 to 23 may include, wherein the patient metric is a subjective pain metric.

In Example 25, the subject matter of any one of Examples 16 to 24 may include, wherein the subjective pain metric is obtained from querying the subject.

In Example 26, the subject matter of any one of Examples 16 to 25 may include, wherein the space domain includes at least one parameter related to fractionalization or polarity.

In Example 27, the subject matter of any one of Examples 16 to 26 may include, wherein the candidate neuromodulation parameter set is designed to produce a pulse burst that varies pulse-by-pulse.

In Example 28, the subject matter of any one of Examples 16 to 27 may include, wherein the memory device further comprises instructions, which when executed by the processor, cause the processor to: receive from a user, a selected parameter; and focus the machine learning algorithm using the selected parameter.

Example 29 includes subject matter (such as a method, means for performing acts, machine readable medium including instructions that when performed by a machine cause the machine to performs acts, or an apparatus to perform) comprising: accessing, at a computerized system, a patient metric of a subject; using the patient metric as an input to a machine learning algorithm executing on the computerized system, the machine learning algorithm to search a plurality of neuromodulation parameter sets and to identify a candidate neuromodulation parameter set of the plurality of neuromodulation parameter sets, the candidate neuromodulation parameter set designed to produce a non-regular waveform that varies over a time domain and a space domain; programming a neuromodulator using the candidate neuromodulation parameter set to stimulate the subject; and producing the non-regular waveform.

In Example 30, the subject matter of Example 29 may include, wherein the subject is a patient.

In Example 31, the subject matter of any one of Examples 29 to 30 may include, wherein the subject is an animal from a preclinical trial.

In Example 32, the subject matter of any one of Examples 29 to 31 may include, wherein the machine learning algorithm comprises a genetic algorithm.

In Example 33, the subject matter of any one of Examples 29 to 32 may include, wherein the candidate neuromodulation parameter set is designed to produce a pulse burst that varies pulse-by-pulse.

In Example 34, the subject matter of any one of Examples 29 to 33 may include, receiving from a user, a selected parameter; and focusing the machine learning algorithm using the selected parameter.

Example 35 includes subject matter (such as a machine-readable medium, computer-readable medium, or other stored instructions) comprising instructions, which when executed, cause a device to: access a patient metric of a subject; use the patient metric as an input to a machine learning algorithm, the machine learning algorithm to search a plurality of neuromodulation parameter sets and to identify a candidate neuromodulation parameter set of the plurality of neuromodulation parameter sets, the candidate neuromodulation parameter set designed to produce a non-regular waveform that varies over a time domain and a space domain; and program a neuromodulator using the candidate neuromodulation parameter set to stimulate the subject.

DETAILED DESCRIPTION

Various embodiments described herein involve spinal cord modulation. A brief description of the physiology of the spinal cord and related apparatus is provided herein to assist the reader.FIG. 1illustrates, by way of example, a portion of a spinal cord100including white matter101and gray matter102of the spinal cord. The gray matter102includes cell bodies, synapse, dendrites, and axon terminals. Thus, synapses are located in the gray matter. White matter101includes myelinated axons that connect gray matter areas. A typical transverse section of the spinal cord includes a central “butterfly” shaped central area of gray matter102substantially surrounded by an ellipse-shaped outer area of white matter101. The white matter of the dorsal column (DC)103includes mostly large myelinated axons that form afferent fibers that run in an axial direction. The dorsal portions of the “butterfly” shaped central area of gray matter are referred to as dorsal horns (DH)104. In contrast to the DC fibers that run in an axial direction, DH fibers can be oriented in many directions, including perpendicular to the longitudinal axis of the spinal cord. Examples of spinal nerves105are also illustrated, including a dorsal root (DR)105, dorsal root ganglion107and ventral root108. The dorsal root105mostly carries sensory signals into the spinal cord, and the ventral root functions as an efferent motor root. The dorsal and ventral roots join to form mixed spinal nerves105.

SCS has been used to alleviate pain. A therapeutic goal for conventional SCS programming has been to maximize stimulation (i.e., recruitment) of the DC fibers that run in the white matter along the longitudinal axis of the spinal cord and minimal stimulation of other fibers that run perpendicular to the longitudinal axis of the spinal cord (dorsal root fibers, predominantly), as illustrated inFIG. 1. The white matter of the DC includes mostly large myelinated axons that form afferent fibers. While the full mechanisms of pain relief are not well understood, it is believed that the perception of pain signals is inhibited via the gate control theory of pain, which suggests that enhanced activity of innocuous touch or pressure afferents via electrical stimulation creates interneuronal activity within the DH of the spinal cord that releases inhibitory neurotransmitters (Gamma-Aminobutyric Acid (GABA), glycine), which in turn, reduces the hypersensitivity of wide dynamic range (WDR) sensory neurons to noxious afferent input of pain signals traveling from the dorsal root (DR) neural fibers that innervate the pain region of the patient, as well as treating general WDR ectopy. Consequently, the large sensory afferents of the DC nerve fibers have been targeted for stimulation at an amplitude that provides pain relief. Current implantable neuromodulation systems typically include electrodes implanted adjacent, i.e., resting near, or upon the dura, to the dorsal column of the spinal cord of the patient and along a longitudinal axis of the spinal cord of the patient.

Activation of large sensory DC nerve fibers also typically creates the paresthesia sensation that often accompanies conventional SCS therapy. Some embodiments deliver therapy where the delivery of energy is perceptible due to sensations such as paresthesia. Although alternative or artifactual sensations, such as paresthesia, are usually tolerated relative to the sensation of pain, patients sometimes report these sensations to be uncomfortable, and therefore, they can be considered an adverse side-effect to neuromodulation therapy in some cases. Some embodiments deliver sub-perception therapy that is therapeutically effective to treat pain, for example, but the patient does not sense the delivery of the modulation field (e.g. paresthesia). Sub-perception therapy may include higher frequency modulation (e.g. about 1500 Hz or above) of the spinal cord that effectively blocks the transmission of pain signals in the afferent fibers in the DC. Some embodiments herein selectively modulate DH tissue or DR tissue over DC tissue to provide sub-perception therapy. Such selective modulation is not delivered at these higher frequencies. For example, the selective modulation may be delivered at frequencies less than 1,200 Hz. The selective modulation may be delivered at frequencies less than 1,000 Hz in some embodiments. In some embodiments, the selective modulation may be delivered at frequencies less than 500 Hz. In some embodiments, the selective modulation may be delivered at frequencies less than 350 Hz. In some embodiments, the selective modulation may be delivered at frequencies less than 130 Hz. The selective modulation may be delivered at low frequencies (e.g. as low as 2 Hz). The selective modulation may be delivered even without pulses (e.g. 0 Hz) to modulate some neural tissue. By way of example and not limitation, the selective modulation may be delivered within a frequency range selected from the following frequency ranges: 2 Hz to 1,200 Hz; 2 Hz to 1,000 Hz, 2 Hz to 500 Hz; 2 Hz to 350 Hz; or 2 Hz to 130 Hz. Systems may be developed to raise the lower end of any these ranges from 2 Hz to other frequencies such as, by way of example and not limitation, 10 Hz, 20 Hz, 50 Hz or 100 Hz. By way of example and not limitation, it is further noted that the selective modulation may be delivered with a duty cycle, in which stimulation (e.g. a train of pulses) is delivered during a Stimulation ON portion of the duty cycle, and is not delivered during a Stimulation OFF portion of the duty cycle. By way of example and not limitation, the duty cycle may be about 10%±5%, 20%±5%, 30%±5%, 40%±5%, 50%±5% or 60%±5%. For example, a burst of pulses for 10 ms during a Stimulation ON portion followed by 15 ms without pulses corresponds to a 40% duty cycle.

FIG. 2illustrates an embodiment of a neuromodulation system. The illustrated system210includes electrodes211, a modulation device212, and a programming device213. The electrodes211are configured to be placed on or near one or more neural targets in a patient. The modulation device212is configured to be electrically connected to electrodes211and deliver neuromodulation energy, such as in the form of electrical pulses, to the one or more neural targets though electrodes211. The delivery of the neuromodulation is controlled by using a plurality of modulation parameters, such as modulation parameters specifying the electrical pulses and a selection of electrodes through which each of the electrical pulses is delivered. In various embodiments, at least some parameters of the plurality of modulation parameters are programmable by a user, such as a physician or other caregiver. The programming device213provides the user with accessibility to the user-programmable parameters. In various embodiments, the programming device213is configured to be communicatively coupled to modulation device via a wired or wireless link. In various embodiments, the programming device213includes a graphical user interface (GUI)214that allows the user to set and/or adjust values of the user-programmable modulation parameters.

FIG. 3illustrates an embodiment of a modulation device312, such as may be implemented in the neuromodulation system210ofFIG. 2. The illustrated embodiment of the modulation device312includes a modulation output circuit315and a modulation control circuit316. Those of ordinary skill in the art will understand that the neuromodulation system210may include additional components such as sensing circuitry for patient monitoring and/or feedback control of the therapy, telemetry circuitry and power. The modulation output circuit315produces and delivers neuromodulation pulses. The modulation control circuit316controls the delivery of the neuromodulation pulses using the plurality of modulation parameters. The lead system317includes one or more leads each configured to be electrically connected to modulation device312and a plurality of electrodes311-1to311-N distributed in an electrode arrangement using the one or more leads. Each lead may have an electrode array consisting of two or more electrodes, which also may be referred to as contacts. Multiple leads may provide multiple electrode arrays to provide the electrode arrangement. Each electrode is a single electrically conductive contact providing for an electrical interface between modulation output circuit315and tissue of the patient, where N≥2. The neuromodulation pulses are each delivered from the modulation output circuit315through a set of electrodes selected from the electrodes311-1to311-N. The number of leads and the number of electrodes on each lead may depend on, for example, the distribution of target(s) of the neuromodulation and the need for controlling the distribution of electric field at each target. In one embodiment, by way of example and not limitation, the lead system includes two leads each having eight electrodes.

The neuromodulation system may be configured to modulate spinal target tissue or other neural tissue. The configuration of electrodes used to deliver electrical pulses to the targeted tissue constitutes an electrode configuration, with the electrodes capable of being selectively programmed to act as anodes (positive), cathodes (negative), or left off (zero). In other words, an electrode configuration represents the polarity being positive, negative, or zero. Other parameters that may be controlled or varied include the amplitude, pulse width, and rate (or frequency) of the electrical pulses. Each electrode configuration, along with the electrical pulse parameters, can be referred to as a “modulation parameter set.” Each set of modulation parameters, including fractionalized current distribution to the electrodes (as percentage cathodic current, percentage anodic current, or off), may be stored and combined into a modulation program that can then be used to modulate multiple regions within the patient.

The number of electrodes available combined with the ability to generate a variety of complex electrical pulses, presents a huge selection of available modulation parameter sets to the clinician or patient. For example, if the neuromodulation system to be programmed has sixteen electrodes, millions of modulation parameter sets may be available for programming into the neuromodulation system. Furthermore, for example SCS systems may have thirty-two electrodes which exponentially increases the number of modulation parameters sets available for programming. To facilitate such selection, the clinician generally programs the modulation parameters sets through a computerized programming system to allow the optimum modulation parameters to be determined based on patient feedback or other means and to subsequently program the desired modulation parameter sets. A closed-loop mechanism may be used to identify and test modulation parameter sets, receive patient or clinician feedback, and further revise the modulation parameter sets to attempt to optimize stimulation paradigms for pain relief. The patient or clinician feedback may be objective and/or subjective metrics reflecting pain, paresthesia coverage, or other aspects of patient satisfaction with the stimulation.

FIG. 4illustrates an embodiment of a programming device413, such as may be implemented as the programming device213in the neuromodulation system ofFIG. 2. The programming device413includes a storage device418, a programming control circuit419, and a GUI414. The programming control circuit419generates the plurality of modulation parameters that controls the delivery of the neuromodulation pulses according to the pattern of the neuromodulation pulses. In various embodiments, the GUI414includes any type of presentation device, such as interactive or non-interactive screens, and any type of user input devices that allow the user to program the modulation parameters, such as touchscreen, keyboard, keypad, touchpad, trackball, joystick, and mouse. The storage device418may store, among other things, modulation parameters to be programmed into the modulation device. The programming device413may transmit the plurality of modulation parameters to the modulation device. In some embodiments, the programming device413may transmit power to the modulation device (e.g., modulation device312ofFIG. 3). The programming control circuit419may generate the plurality of modulation parameters. In various embodiments, the programming control circuit419may check values of the plurality of modulation parameters against safety rules to limit these values within constraints of the safety rules.

In various embodiments, circuits of neuromodulation, including its various embodiments discussed in this document, may be implemented using a combination of hardware, software, and firmware. For example, the circuit of GUI414, modulation control circuit316, and programming control circuit419, including their various embodiments discussed in this document, may be implemented using an application-specific circuit constructed to perform one or more particular functions or a general-purpose circuit programmed to perform such function(s). Such a general-purpose circuit includes, but is not limited to, a microprocessor or a portion thereof, a microcontroller or a portion thereof, and a programmable logic circuit or a portion thereof.

FIG. 5illustrates, by way of example, an implantable neuromodulation system and portions of an environment in which system may be used. The system is illustrated for implantation near the spinal cord. However, neuromodulation system may be configured to modulate other neural targets. The system520includes an implantable system521, an external system522, and a telemetry link523providing for wireless communication between implantable system521and external system522. The implantable system521is illustrated as being implanted in the patient's body. The implantable system521includes an implantable modulation device (also referred to as an implantable pulse generator, or IPG)512, a lead system517, and electrodes511. The lead system517includes one or more leads each configured to be electrically connected to the modulation device512and a plurality of electrodes511distributed in the one or more leads. In various embodiments, the external system402includes one or more external (non-implantable) devices each allowing a user (e.g. a clinician or other caregiver and/or the patient) to communicate with the implantable system521. In some embodiments, the external system522includes a programming device intended for a clinician or other caregiver to initialize and adjust settings for the implantable system521and a remote control device intended for use by the patient. For example, the remote control device may allow the patient to turn a therapy on and off and/or adjust certain patient-programmable parameters of the plurality of modulation parameters. The remote control device may also provide a mechanism for the patient to provide feedback on the operation of the implantable neuromodulation system. Feedback may be metrics reflecting perceived pain, effectiveness of therapies, or other aspects of patient comfort or condition.

The neuromodulation lead(s) of the lead system517may be placed adjacent, e.g., resting near, or upon the dura, adjacent to the spinal cord area to be stimulated. For example, the neuromodulation lead(s) may be implanted along a longitudinal axis of the spinal cord of the patient. Due to the lack of space near the location where the neuromodulation lead(s) exit the spinal column, the implantable modulation device512may be implanted in a surgically-made pocket either in the abdomen or above the buttocks, or may be implanted in other locations of the patient's body. The lead extension(s) may be used to facilitate the implantation of the implantable modulation device512away from the exit point of the neuromodulation lead(s).

Application of Machine Learning to Optimize Stimulation Parameters

Identifying optimized stimulation patterns to relive pain in neuromodulation (e.g., SCS) is challenging because the parameter space is so large, objective metrics to assess patient pain are not used, and response latency may be long (e.g., up to days). A machine learning system may be used in a closed-loop hybrid clinical or pre-clinical process to develop optimized stimulation patterns for pain management. The stimulation patterns may be modulated in both the time and space domains.

An initial set of stimulation patterns may be generated from a domain of all available stimulation patterns. The initial set may be based on preclinical pain trials. The initial set may be obtained using one or more machine learning or optimization algorithms to search for and identify effective patterns.

In the clinical system, a patient may be provided one or more stimulation patterns, which may be tested by the patient with or without clinician supervision. Objective pain metrics, subjective pain metrics, or both objective and subjective pain metrics may be received from the patient, which are used in the machine learning or optimization algorithms to develop further sets of patterns. Objective pain metrics include those that are physiologically expressed, such as EEG activity, heart rate, heart rate variability, galvanic skin response, or the like. Subjective pain metrics may be provided by the patient and be expressed as “strong pain,” “lower pain,” or numerically in a range, for example. The pain metrics may be communicated using various communication mechanisms, such as wireless networks, tethered communication, short-range telemetry, or combinations of such mechanisms. The patient may manually input some information (e.g., subjective pain scores).

A non-exhaustive list of pain metrics is provided herein. One example of a pain metric is EEG activity (e.g., Theta activity in the somatosensory cortex and alpha and gamma activity in the prefrontal cortex have been shown to correlate with pain). Another example pain metric is fMRI (activity in the anterior cingulate cortex and insula have been shown to correlate with changes in chronic pain). Another example pain metric is fMRI (activity in the pain matrix, which consists of the thalamus, primary somatosensory cortex, anterior cingulate cortex, prefrontal cortex, and cerebellum and is activated in pain conditions). Another example pain metric is heart rate variability, galvanic skin response, cortisol level, and other measures of autonomic system functioning (autonomic system health has been shown to correlate with pain). Another example pain metric is physical activity (amount of physical activity has been shown to correlate with pain). Another example pain metric is pain scores (may be inputted through an interface where the patient selects a point on a visual analog scale, or clicks a number on a numerical rating scale). Another example pain metric is quantitative sensory testing [e.g., spatial discrimination (two-point, location, diameter), temporal discrimination, detection threshold (mechanical, thermal, electrical), pain threshold (mechanical, thermal, electrical), temporal summation, thermal grill] (QST measures have been shown to correlate with pain). Another example pain metric is somatosensory evoked potentials, contact heat evoked potentials (these have been shown to be correlated with pain). Another example pain metric is H-reflex, nociceptive flexion reflex (these have been shown to be reduced by SCS). Another example pain metric is conditioned place preference (e.g., in one chamber, stimulate with one paradigm1, in other chamber, stimulate with paradigm2. The chamber where the animal spends the most time wins and continues to the next round.). Another example pain metric is local field potential recordings in the pain matrix (recordings of neural activity in these areas are possible with invasive electrodes in a preclinical model).

Some pain metrics are primarily preclinical in nature (e.g., conditioned place preference and local field potential recordings), while others are primarily clinical in nature (e.g., pain scores and quantitative sensory testing). However, it is understood that the pain metrics may be obtained in either preclinical or clinical settings.

Pain metrics may be continuously or repeatedly collected from patients and fed into the machine learning or optimization algorithms to refine or alter the stimulation patterns. For example, the patients may interact with a programmer, remote control, bedside monitor, or other patient device to record physical condition, pain, medication dosages, etc. The patient device may be wired or wirelessly connected to the system with the machine learning system. This closed-loop mechanism provides an advantage of reducing the search domain during repeated iterations of the machine learning or optimization algorithm. By reducing the search domain, a clinician is able to more quickly identify efficacious patterns and a patient may be subjected to shorter programming sessions, which produce less discomfort.

The physical system may take on many different forms. Data collected from the patient or pre-clinical model may be measured using wearable sensors (e.g., heart rate monitor, accelerometer, EEG headset). The pain metrics requiring manual input may be entered via remote control or other external device used by the patient (e.g. cellular phone).

FIG. 6illustrates, by way of example, an embodiment of an SCS system600. The SCS system600generally comprises a plurality of neurostimulation leads12(in this case, two percutaneous leads12aand12b), an implantable pulse generator (IPG)14, an external remote control (RC)16, a Clinician's Programmer (CP)18, an External Trial Stimulator (ETS)20, and an external charger22.

The IPG14is physically connected via two lead extensions24to the neurostimulation leads12, which carry a plurality of electrodes26arranged in an array. In the illustrated embodiment, the neurostimulation leads12are percutaneous leads, and to this end, the electrodes26are arranged in-line along the neurostimulation leads12. The number of neurostimulation leads12illustrated is two, although any suitable number of neurostimulation leads12can be provided, including only one. Alternatively, a surgical paddle lead can be used in place of one or more of the percutaneous leads. As will also be described in further detail below, the IPG14includes pulse generation circuitry that delivers electrical stimulation energy in the form of a pulsed electrical waveform (i.e., a temporal series of electrical pulses) to the electrode array26in accordance with a set of stimulation parameters. The IPG14and neurostimulation leads12can be provided as an implantable neurostimulation kit, along with, e.g., a hollow needle, a stylet, a tunneling tool, and a tunneling straw.

The ETS20may also be physically connected via percutaneous lead extensions28or external cable30to the neurostimulation lead12. The ETS20, which has similar pulse generation circuitry as the IPG14, also delivers electrical stimulation energy in the form of a pulsed electrical waveform to the electrode array26in accordance with a set of stimulation parameters. The major difference between the ETS20and the IPG14is that the ETS20is a non-implantable device that is used on a trial basis after the neurostimulation lead12has been implanted and prior to implantation of the IPG14, to test the responsiveness of the stimulation that is to be provided. Thus, any functions described herein with respect to the IPG14can likewise be performed with respect to the ETS20.

The RC16may be used to telemetrically control the ETS20via a bi-directional RF communications link32. Once the IPG14and stimulation leads12are implanted, the RC16may be used to telemetrically control the IPG14via a bi-directional RF communications link34. Such control allows the IPG14to be turned on or off and to be programmed with different stimulation programs after implantation. Once the IPG14has been programmed, and its power source has been charged or otherwise replenished, the IPG14may function as programmed without the RC16being present.

The CP18provides the user detailed stimulation parameters for programming the IPG14and ETS20in the operating room and in follow-up sessions. The CP18may perform this function by indirectly communicating with the IPG14or ETS20, through the RC16, via an IR communications link36. Alternatively, the CP18may directly communicate with the IPG14or ETS20via an RF communications link (not shown).

The external charger22is a portable device used to transcutaneously charge the IPG14via an inductive link38. Once the IPG14has been programmed, and its power source has been charged by the external charger22or otherwise replenished, the IPG14may function as programmed without the RC16or CP18being present.

For the purposes of this specification, the terms “neurostimulator,” “stimulator,” “neurostimulation,” and “stimulation” generally refer to the delivery of electrical energy that affects the neuronal activity of neural tissue, which may be excitatory or inhibitory; for example by initiating an action potential, inhibiting or blocking the propagation of action potentials, affecting changes in neurotransmitter/neuromodulator release or uptake, and inducing changes in neuro-plasticity or neurogenesis of tissue. For purposes of brevity, the details of the RC16, ETS20, and external charger22will not be described herein.

The algorithm may reside on the CP, the IPG, the ETS, the RC or other external device used by the patient, or in the cloud or remote servers connected to patient external via Wi-Fi, Bluetooth, cellular data, or other wired/wireless scheme. There may be a GUI on the CP, remote control, or other external device, that enables selection of algorithm as well as manual input. Training of the algorithm may take place in the clinic or in daily life, and may be set to be execute continually or only at certain times. Optimization data may be stored in the cloud so that optimized patterns and history can be transferred when the patient moves from trial to permanent implant and also if the IPG is replaced.

FIG. 7illustrates, by way of example, an embodiment of data and control flow in a system that utilizes machine learning to optimize neurostimulation patterns. A preclinical pain model may be used or alternatively, the feedback from a previous stimulation may be used to initialize the system (block700). The patient is stimulated and one or more patient metrics are obtained (block702). The patient metrics may be obtained via passive or active participation with the patient. For example, the patient metrics may be derived from sensing and correlating physiological changes or states in the patient. In this case, the patient provides the patient metrics passively—without active participation. Alternatively, the patient may be prompted (e.g., by a clinician or with an electronic graphical user interface) to provide a patient metric. Patient metrics gathered during passive participation may be referred to as objective pain measurements because they are typically physiological responses that are mostly uncontrollable by the patient, such as EEG activity or heart rate variability. Some patient metrics gathered via active participation with the patient may be referred to as subjective patient metrics, where the patient is asked to describe the pain. The patient metrics may include various aspects of pain, such as the severity as measured with a numerical value, the location(s) of pain, the sensation of pain (e.g., numbness, shape acute pain, throbbing, etc.), the duration of pain, or other aspects of pain. The patient metrics may also include results of questionnaires, responses to queries about a general state of wellness, results of memory tests (e.g., working memory tasks), rating scales, and the like.

The patient metric(s) are used as input to the machine learning or optimization algorithm (block704). Various algorithms may be used, such as genetic algorithms, neural networks, or reinforcement learning strategies (e.g., Q-learning, Temporal Difference learning, or Markov decision processes). The machine learning or optimization algorithm may modify one or more variables, such as time delays, amplitudes of one or more pulses in a stimulation train, pulse shape, fractionalization, or pulse-by-pulse changes in spatial location of stimulation.

In a preclinical setting, objective patient metrics may be gathered from a study (e.g., an animal study) and used to identify a subset of parameter settings that may be used as an initial setting in a clinical environment. One benefit of the preclinical-machine learning closed-loop system is that much more time can be spent on optimization. Another benefit is that more invasive measures of pain can be used (e.g. electrodes implanted in the brain). As such, some embodiments use the preclinical system first to identify stimulation paradigms that are better than others, and then use the clinical system to further hone these identified stimulation patterns in patients. In addition to genetic algorithms, machine learning algorithms that require large data sets such as neural networks may be employed in the preclinical model. The preclinical model will also enable exploration in a larger parameter space due to the increased amount of testing time. This information can be used to characterize stimulation parameters and design the optimization algorithm in the clinical system.

At block704, the machine learning or optimization algorithm may be a genetic algorithm. Stimulation parameters that describe arbitrary waveforms in the time and space domain are modulated during the optimization search. An objective function uses patient metrics to evaluate the stimulation waveform tested. The patient metrics may be objective or subjective feedback. Based on the objective function values, new stimulation parameters are selected for subsequent testing. The objective function may include several feedback components such as the pain score and a side effects score. The objective function may then be customized to fit the needs of individual patients by adjusting the objective function component weights to emphasize one aspect more than the other (e.g., pain or side effects). The genetic algorithm seeks to minimize the objective function value due to changing the stimulation parameters during the search.

In an embodiment, the genetic algorithm resides and executes on the remote control. The patient is able to decide a convenient time to use the “optimize” mode to modulate their stimulation. A stimulation waveform is selected and the patient can increase the amplitude to their comfort level. For each iteration of testing, the amount of time the stimulation is administered may be patient adjustable, for example 5 minutes. After 5 minutes, the amplitude may be decreased to 0, a new stimulation is then queued based on the patient's feedback, the amplitude is increased to the patient's comfort level, and the patient's feedback is recorded. This feedback mechanism may continue allowing the patient to optimize their stimulation. After the patient opts to finish or termination criteria is met (e.g., the objective function is satisfied), the best stimulation measured by their feedback metric may be saved and available for use. The patient may later return and optimize further using their remote control at their convenience.

Output of algorithm-derived stimulation patterns are produced at block706. These stimulation patterns are capable of having pulse-by-pulse changes in timing such as pulse start time, pulse end time, pulse duration, and pulse interval, pulse-by-pulse changes of pulse amplitudes and pulse shape, and of having pulse-by-pulse changes in spatial location of stimulation including selection of active electrodes for each pulse, polarity and fractionalization of the modulation energy using the active electrodes. The patterns may be regular, repeating patterns, or may be non-regular in at least some of the modulation parameters (amplitude, pulse width, pulse interval, pulse polarity, and the like) in pulses groups and in groups of pulses groups.

Although some of the present disclosure discusses SCS, it is understood that other types of neuromodulation may be controlled using the systems and methods described herein. Follow is a non-exhaustive list of potential use cases.

In a neuropsychiatric use case, a number of neuropsychiatric disorders exist including major depression, obsessive compulsive disorder, addiction, anorexia, bipolar disorder, Tourette's syndrome, and the like. Brain stimulation has been investigated as a method to reduce or eliminate symptoms of such disorders.

Neuropsychiatric symptoms may be difficult to quantitatively measure, but correlations have been discovered between depression and EEG features, and between depression and autonomic features such as heart rate variability. Other autonomic measures that may be candidates include properties of respiration, galvanic skin response, or other measures related to autonomic balance. Because EEG is a candidate, it is expected that local field potential signatures also likely contain information that may be used as a feedback signal when properly processed.

As such, in embodiments, a patient with a neuropsychiatric disorder, such as depression, may receive a stimulator of a neural tissue or an endocrine tissue (e.g., brain stimulator, nerve stimulator, or spinal cord stimulator). A machine learning method may be used to optimize a pattern of stimulation, where optimization includes minimizing or maximizing a cost function of one or more objective metrics. In a related embodiment, subjective metrics, such as scores from an e-diary that a patient or caregiver records electronically into the system, may be used. Subjective measurement scores may be from standard questionnaires (e.g., the Hamilton depression rating scale) or responses to queries about general well-being or state of wellness. In yet another related embodiment, both subjective scores and objective scores may be combined and used. In yet another related embodiment, the subjective responses constitute a “ground truth” and the objective metrics may be used in conjunction with that ground truth to adapt the machine learned model. Such hybrid use of subjective and objective metrics are useful because biological systems are often not static and change over time. Such an embodiment allows the system to respond to phenomena like accommodation and habituation to the benefit of the patient.

In a cardiac function use case, it is known that stimulation of a peripheral nerve (e.g., vagus), spinal cord (e.g., SCS for angina), or DBS may affect cardiac function. Some direct and indirect measures of cardiac function are able to be measured objectively, including heart rate, heart rate variability, blood oxygen perfusion, blood pressure, patient activity, EKG properties. In embodiments, electrodes in the heart for patients with a pacemaker or defibrillator may enable sensing of atrial-ventricular coordination, cardiac EMGs, and other measures. Such measurements may serve as a candidate output target of the machine learning algorithm in addition to the neurostimulator. In some cases, these metrics may be obtained from wearable devices.

As it relates to the present disclosure, in an embodiment, a patient with a cardiac condition, such as hypertension, may receive a stimulator of a neural tissue or an endocrine tissue, and a machine learning method may be used to optimize a pattern of stimulation, where optimization includes minimizing or maximizing a cost function comprised of one or more objective metrics such as blood pressure. In a related embodiment, a preclinical model of a cardiac condition such as hypertension is configured with a stimulator of neural tissue and an optimization algorithm (such as machine learning) that has access to the quantitative metric (e.g., blood pressure) and control of a stimulation pattern and is used to optimize the stimulation pattern to achieve an ideal blood pressure. This stimulation pattern or a similar pattern may then be subsequently used as an initial setting in a patient being stimulated for the same or a related condition using a same or related stimulation target (e.g., neural tissue).

In an epilepsy use case, it is understood that epilepsy is a central nervous system disorder in which neuronal activity in the brain becomes abnormal without apparent reason, and manifests in recurring seizures or periods of unusual behavior and sensations and loss of consciousness. Seizures may be of variable duration and may vary from being nearly undetectable to significant shaking. More than 30% of patients suffer persistent seizures despite maximum antiepileptic drug therapy. There is a pressing need for alternative treatments.

Deep brain stimulation (DBS), spinal cord stimulation (SCS), vagus nerve stimulation (VNS), and peripheral nerve stimulation (PNS) are all subjects of active investigation for potential therapeutic benefit in epilepsy. In addition to their potential for treatment of patients who are refractory to treatment with drugs, these neurostimulation treatments may be used to use physiological signals of seizure to deliver therapy at opportune times, thereby reducing unnecessary stimulation and potentially improving therapy.

Patient metrics include physiological signals that are of interest in this context, such as EEG, near-infrared spectroscopy (NIRS) or other non-invasive measures of brain activity, EMG or alternate measures of muscle activity, EKG or other measures of heart activity, and clinician, patient or caregiver reports. These measurements either singly or together may be used to derive a severity metric for the pathology. For instance, spectral analysis of the scalp EEG can be used to determine seizure onset and offset, and the duration and frequency of seizures, and the latter can be used to derive an epilepsy severity metric. The metric resulting from this analysis may be used to quantify the patient's response to stimulation settings and this data is then used as input to a machine learning algorithm. The output of this algorithm includes a candidate stimulation setting that reduces disease severity.

In an overactive bladder (OAB) use case, certain nerve stimulation may be used as a treatment. OAB is a cluster of symptoms related to urinary function, the chief among which is a sudden urge to urinate. Some 30% of American men and 40% of American women suffer from overactive bladder. Behavioral interventions, medications, bladder injections, and nerve stimulation are all used as treatments of OAB. Evidence also indicates that both deep brain and spinal cord stimulation may have efficacy in treating some OAB symptoms.

OAB may be quantified by objective and subjective metrics, including void frequency, void volume, frequency of painful voids, subjective patient report of well-being and urinary urgency, etc. Effect of nerve stimulation on OAB severity may be assessed using these criteria and resulting data is fed as input to a machine learning algorithm. The output of this algorithm includes a candidate stimulation setting that may reduce symptom severity.

Other types of disorders may be treated by neuromodulation, such as movement disorders or cognitive disorders. Movement disorders include ailments like Parkinson's disease. Parkinson's disease is characterized by the cardinal symptoms of tremor, bradykinesia, and rigidity. Deep brain stimulation of the thalamus, STN (subthalamic nucleus), or GPi (Globus pallidus) is often used to improve these symptoms. However, adjustment of DBS by a neurologist is traditionally done through a serial process where the neurologist makes a program adjustment, observes a certain symptom (e.g., tremor, arm rigidity), task (e.g., finger-tapping, rapidly alternating movement), or side effect (e.g., dysarthria, muscle twitches), and then makes further adjustments. This is time-consuming and may fail to optimize the stimulation settings across all symptoms in all areas of the body.

A machine learning algorithm may be used to improve DBS programming by taking a plurality of assessments (e.g., local field potential measurements from the implanted leads, anatomical placement of the leads based on MRI and CT images, motor diary information, unified Parkinson's disease rating scale (UPDRS) scores, quantitative assessments using wearable accelerometers, speech recordings, timed motor tests), at a plurality of neurostimulation settings, and then recommending a setting based on the entirety of these data, rather than making serial adjustments after one or two observations.

Similar approaches may be used for other movement disorders such as dystonia, which is further complicated by the slow onset of DBS response and resulting difficulty of adjusting using serial observations; or essential tremor, which may show optimal tremor control for different parts of the body at slightly different stimulation settings.

Cognitive disorders include ailments such as Alzheimer's disease or Parkinson's-related dementia. DBS is also used to affect structures, such as the fornix, nucleus basalis Meynert, or entorhinal cortex. Cognitive performance is complex and may be assessed through a wide variety of methods, including working memory tasks (e.g. N-back tests, mini-cog), questionnaires and rating scales (mini-mental state examination (MMSE), Mattis Dementia Rating Scale, Alzheimer's Disease Assessment Scale—Cognition (ADAS-cog), etc.), brain imaging, mood assessments, and dual motor-cognitive tasks. Furthermore, like movement disorder evaluations, cognitive performance may be time-consuming to assess and does not lend itself to programming through serial observations. As such, a machine learning algorithm may be used to improve DBS programming by taking a plurality of assessments (e.g. local field potential measurements from the implanted leads, anatomical placement of the leads based on MRI and CT images, working memory tasks, questionnaires, rating scales, etc.) at a plurality of neurostimulation settings, and then recommending a setting based on the entirety of these data, rather than making serial adjustments after one or two observations.

FIG. 8illustrates, by way of example, another embodiment of data and control flow in a system that utilizes machine learning to optimize neurostimulation patterns.

At stage800, initial input parameters are accessed to construct a waveform (stage802). Input parameters may include active contact fractionalization and waveform characteristics, such as amplitude, time delay, pulse width, or shape of each pulse in a stimulation pulse train. The constructed waveform is used to stimulate a patient via an implanted device (stage804).

At stage806, the patient's feedback is received. Depending on the disease or disorder being treated by the neuromodulation, the patient feedback may include various objective or subjective metrics. Objective metrics, such as the patient's EEG, medication therapy, heart rate, heart rate variability, blood oxygen perfusion, blood pressure, patient activity scores, seizure frequency or duration, seizure duration, near-infrared spectroscopy (NIRS) of brain activity, EMG or other measures of muscle activity, assessments using wearable accelerometers, timed motor tests, memory test, mood assessments, brain imaging, dual motor cognitive tasks, and the like may be used. Subjective patient feedback may include responses to questionnaires, summaries of diaries, clinician reports, or rating scales (e.g., MMSE, Mattis Dementia Rating Scale, ADAS-cog (Alzheimer's Disease Assessment Scale.—Cognition), etc.). Thus, in related embodiments, the patient feedback may be a patient metric, such as a pain score or biomarkers.

At stage808, the patient's feedback is analyzed to determine whether additional modification to the stimulation parameters is needed. A genetic algorithm is used at stage810to identify one or more stimulation parameters. A genetic algorithm is a general purpose search algorithm based on the principle of evolution. The population of candidate solutions (called individuals, creatures, or phenotypes) to an optimization problem is evolved toward better solutions. Each candidate solution has a set of properties (its chromosomes or genotype) that can be mutated or altered. A genetic algorithm combines the operations of selection (stage810-1), crossover (stage810-2), and mutation (stage810-3) with the goal of finding the best solution. The selection-crossover-mutation process favors better solutions from generation to generation. When the genetic algorithm reaches a termination state (e.g., a maximum number of generations) or has convergence to a solution (e.g., minimizes or maximizes an objective function) (stage810-4), then the stimulation parameters for optimized therapy are considered to be reached (stage812). The optimized stimulation parameters may be used until the patient's condition changes or until the parameters are no longer effective. When the genetic algorithm has not reached termination criteria or convergence state, then a next set of waveform parameters is identified and tested (stage814). The parameters are used to construct a waveform (stage802) and the cycle continues.

In the selection stage (810-1), chromosomes of higher fitness are selected to become the co-founders of the next generation of chromosomes. The probability of a chromosome being selected is based on a fitness function. The fitness function is used to choose which chromosomes will be elected to survive to the next generation. The fitness function is a form of an objective function used in machine learning applications. One example design of a fitness function may be to maximize pain reduction while minimizing the area affected by paresthesia. Examples of fitness selection functions include fitness proportionate selection, Bolzmann selection, tournament selection, rank selection, steady state selection, truncation selection, and local selection.

The crossover stage (810-2) is used to select two chromosomes (e.g., parents) and combine them to produce a new chromosome (e.g., offspring). Various crossover techniques may be used, such as one-point crossover, two-point crossover, or a uniform crossover. A one-point crossover function uses a single point in the parent chromosomes, and uses one parent for the first portion of the offspring chromosome and the other parent for the remainder of the offspring chromosome (after the crossover point). A two-point crossover uses two points on the chromosome to create three sections in the offspring chromosome(s), using one parent for the first and third portion of the offspring chromosome and the other parent for the second portion of the offspring chromosome. Crossover points may be chosen randomly or uniformly among chromosomes. Other functions may be used to combine chromosomes. One example function may be to equally weight each gene in the parent chromosomes and randomly select a gene from one parent chromosome for the offspring chromosome.

The mutation stage (810-3) alters one or more gene values in an offspring chromosome. The gene value may be selected randomly from a chromosome in the general gene pool. The mutated gene value may be constrained to a range (e.g., ±25% of the original gene value). Mutation prevents a selected chromosome population from becoming stagnant.

The termination of the genetic algorithm (810-4) may be achieved by executing a certain number of iterations (e.g., producing a certain number of generations), or by producing a chromosome that meets a fitness requirement. Other ccommon terminating conditions are reaching an allocated budget (computation time/money), reaching a highest ranking solution's fitness or reaching a plateau of chromosomes such that successive iterations no longer produce better results, or combinations of these conditions.

For the purposes of this discussion, a stimulation protocol may be considered as a construction of building blocks beginning with a pulse. A pulse is single waveform and typically has a timescale in the millisecond range. A burst is a sequence of pulses and may have a timescale on the millisecond to second range. A train is a sequence of bursts and may have a timescale of millisecond, seconds, or even minutes depending on the programming used. A programming sequence is a combination of pulses, bursts, and trains. The programming sequence may also include pauses; periods with no electrical stimulation. A programming sequence may be cyclical over short durations or be non-cyclical over a short duration, but repeat over some longer “macropulse” duration.

In a pulse burst or a pulse train, the intervals between pulses may be regular or irregular. In general, the time domain includes stimulation parameters that control the timing, size, or shape of pulses. Time domain parameters include, but are not limited to, the pulse rate, pulse amplitude, pulse shape, pulse width, and interpulse delay (e.g., between bursts or trains).

In addition to the characteristics of the pulses, the location and direction of stimulation may be controlled using stimulation parameters in the space domain. Various spatial domain parameters include, but are not limited to, lead activation (e.g., which lead(s) are active/inactive), electrode activation (e.g., which electrode(s) in a lead are active/inactive) and active contact fractionalization (e.g., of the active electrodes, how much current is supplied to each active electrode in a lead).

FIG. 9illustrates, by way of example, an embodiment of constructing stimulation waveforms in space and time domains. A pulse burst may be represented as a space-time vector900. The space-time vector900includes parameters for a spatial elements902and temporal elements904. In this example, the space-time vector900parameterizes a pulse burst of three pulses. In the space domain, the spatial elements902include a fractionalization for each contact index is indicated. In the time domain, the temporal elements904include the pulse width, time delay, scaled amplitude, and shape index for each of Pulse #1, Pulse #2, and Pulse #3. The temporal element of the shape index is indicated using a “Shape Index,” where for this example, a shape index of 1 indicates a square wave shape, 2 indicates a Gaussian wave shape, and 3 indicates an increasing triangle wave shape. It is understood that more or fewer wave shapes may be used.

Each element in the space-time vector900may be considered a gene and the space-time vector900may be optimized using a genetic algorithm, e.g., by selection, mutation, and crossover of the genes.

Based on the elements (parameters) in the space-time vector900, the contact fractionalization on the second and fourth electrodes (from the distal end)906,908are fractionalized with a “−1” and a “1,” respectively. In the time domain, the pulses910-1,910-2,910-3are produced according the time domain parameters.

FIG. 10illustrates, by way of example, an embodiment of a system1000that utilizes machine learning to optimize neurostimulation patterns. The system1000may take on one of many forms. The system1000may be a remote control or other external device used by a patient or clinician. Alternatively, the system1000may be a server or cloud-based device, a network appliance, or other networked device connected via a network (or combination of networks) to a user device. The networks may include local, short-range, or long-range networks, such as Bluetooth, cellular, Wi-Fi, or other wired or wireless networks.

The system1000includes a processor1002and a memory1004. The processor may be any single processor or group of processors that act cooperatively. The memory1004may be any type of memory, including volatile or non-volatile memory. The memory1004may include instructions, which when executed by the processor1002, cause the processor1002to access a patient metric of a subject. The patient metric may be stored in memory1004or remote from the system1000. For example, the patient metric may be accessed from a remote storage system using various access protocols. In an embodiment, the subject is a patient. In such an embodiment, the patient may be using the system1000in a clinical context. In another embodiment, the subject is an animal from a preclinical trial. The patient metric may be used to quantify a level of pain or discomfort felt by the patient and obtained via passive or active participation with the patient.

The processor1002may further use the patient metric as an input to a machine learning algorithm, the machine learning algorithm to search a plurality of neuromodulation parameter sets and to identify a candidate neuromodulation parameter set of the plurality of neuromodulation parameter sets, the candidate neuromodulation parameter set designed to produce a non-regular waveform that varies over a time domain and a space domain. In an embodiment, the machine learning algorithm comprises a genetic algorithm. In such an embodiment, the neuromodulation parameter set may be represented as genes in a chromosome. By using selection, crossover, and mutation, the genetic algorithm may derive an optimal neuromodulation parameter set from a pool of parameter sets. In another embodiment, the machine learning algorithm comprises a neural network.

The processor1002may further program a neuromodulator using the candidate neuromodulation parameter set to stimulate the subject.

In an embodiment, the system1000is a cloud-based system, and in such an embodiment, programming the neuromodulator includes transmitting the candidate neuromodulation parameter set to a client device of the cloud-based system. The client device may be the CP, the IPG, the ETS, or the RC in an SCS system.

In an embodiment, the patient metric is an objective pain metric. In a further embodiment, the objective pain metric is a physiological indication sensed by a sensor worn by the subject. Examples of objective pain metrics include, but are not limited to EEG activity, heart rate, heart rate variability, galvanic skin response, or the like. It is understood that other types of objective metrics may be used, such as patient activity, EEG, EKG, or EMG measurements, etc., and that the system's use of objective patient data is not limited to objective pain metrics.

In an embodiment, the patient metric is a subjective pain metric. In a further embodiment, the subjective pain metric is obtained from querying the subject. The subject may be a patient, and querying the patient may be performed via a patient device. The patient device may include a graphical user interface within which one or more prompts may be displayed to interrogate the patient about the function, performance, or efficacy of the IPG. The patient may be asked to describe the severity or location of pain, for example, using textual input, a body map, or other user interface elements to indicate subjective pain metrics. It is understood that other types of subjective metrics may be used, such as patient responses to clinicians' queries, patient diaries, questionnaires and rating scales (MMSE, Mattis Dementia Ratings Scale, ADAS-cog), etc., and that the system's use of subjective patient data is not limited to subjective pain metrics.

In an embodiment, the space domain includes at least one parameter related to fractionalization or polarity. Other aspects of the spatial domain may be included, such as those described with respect toFIG. 8.

In an embodiment, the candidate neuromodulation parameter set is designed to produce a pulse burst that varies pulse-by-pulse. Altering the amplitude, pulse shape, pulse width, or other aspects from pulse to pulse may provide additional pain management that are not available with a repeating pulse.

In an embodiment, the processor may further receive from a user, a selected parameter and focus the machine learning algorithm using the selected parameter. The selected parameter may be one or more aspects of a pulse in the time or space domain, such as amplitude.

FIG. 11illustrates, by way of example, an embodiment of a method1100that utilizes machine learning to optimize neurostimulation patterns. At1102, a patient metric of a subject is accessed at a computerized system. In an embodiment, the subject is a patient. In an embodiment, the subject is an animal from a preclinical trial.

At1104, the patient metric is used as an input to a machine learning algorithm executing on the computerized system, the machine learning algorithm to search a plurality of neuromodulation parameter sets and to identify a candidate neuromodulation parameter set of the plurality of neuromodulation parameter sets, the candidate neuromodulation parameter set designed to produce a non-regular waveform that varies over a time domain and a space domain.

In an embodiment, the machine learning algorithm comprises a genetic algorithm. In an embodiment, the machine learning algorithm comprises a neural network.

In an embodiment, the patient metric is an objective pain metric. In a further embodiment, the objective pain metric is a physiological indication sensed by a sensor worn by the subject. In an embodiment, the patient metric is a subjective pain metric. In a further embodiment, the subjective pain metric is obtained from querying the subject.

In an embodiment, the space domain includes at least one parameter related to fractionalization or polarity.

In an embodiment, the candidate neuromodulation parameter set is designed to produce a pulse burst that varies pulse-by-pulse.

At1106, a neuromodulator is programmed using the candidate neuromodulation parameter set to stimulate the subject. In an embodiment, the computerized system is a cloud-based system, and programming the neuromodulator comprises transmitting the candidate neuromodulation parameter set to a client device of the cloud-based system. The neuromodulator may then produce the non-regular waveform.

In an embodiment, the method1100further comprises receiving from a user, a selected parameter and focusing the machine learning algorithm using the selected parameter.

Example computer system1200includes at least one processor1202(e.g., a central processing unit (CPU), a graphics processing unit (GPU) or both, processor cores, compute nodes, etc.), a main memory1204and a static memory1206, which communicate with each other via a link1208(e.g., bus). The computer system1200may further include a video display unit1210, an alphanumeric input device1212(e.g., a keyboard), and a user interface (UI) navigation device1214(e.g., a mouse). In one embodiment, the video display unit1210, input device1212and UI navigation device1214are incorporated into a touch screen display. The computer system1200may additionally include a storage device1216(e.g., a drive unit), a signal generation device1218(e.g., a speaker), a network interface device1220, and one or more sensors (not shown), such as a global positioning system (GPS) sensor, compass, accelerometer, or other sensor.

The storage device1216includes a machine-readable medium1222on which is stored one or more sets of data structures and instructions1224(e.g., software) embodying or utilized by any one or more of the methodologies or functions described herein. The instructions1224may also reside, completely or at least partially, within the main memory1204, static memory1206, and/or within the processor1202during execution thereof by the computer system1200, with the main memory1204, static memory1206, and the processor1202also constituting machine-readable media.