Alternative Paradigm of Selective Vagus Nerve Stimulation Tested on an Isolated Porcine Vagus Nerve

Abstract

Alternative paradigm for spatial and fibre-type selective vagus nerve stimulation (VNS) was developed using realistic structural topography and tested in an isolated segment of a porcine cervical left vagus nerve (LVN). A spiral cuff (cuff) containing a matrix of ninety-nine electrodes was developed for selective VNS. A quasitrapezoidal stimulating pulse (stimulus) was applied to the LVN via an appointed group of three electrodes (triplet). The triplet and stimulus were configured to predominantly stimulate the B-fibres, minimizing stimulation of the A-fibres and by-passing the stimulation of the C-fibres. To assess which fibres made the most probable contribution to the neural response (NR) during selective VNS, the distribution of conduction velocity (CV) within the LVN was considered. Experimental testing of the paradigm showed the existence of certain parameters and waveforms of the stimulus, for which the contribution of the A-fibres to the NR was slightly reduced and that of the B-fibres was slightly enlarged. The cuff provided satisfactory fascicle discrimination in selective VNS as well as satisfactory fascicle discrimination during NR recording. However, in the present stage of development, fibre-type VNS remained rather limited.

Academic Editor:D. J. Moore and Academic Editor:A. K. Moschovakis

ITIS d.o.o. Ljubljana, Center for Implantable Technology and Sensors, Lepi Pot 11, 1000 Ljubljana, Slovenia

Received 26 August 2013; Accepted 2 October 2013; 6 February 2014

This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1. Introduction

In recent decades, considerable scientific and technological efforts have been devoted to the development of neuroprostheses that interface the human autonomic nervous system with electronic implantable devices. Particular attention has been paid to VNS techniques that are to be used to treat, among others, a number of nervous system disorders, neuropsychiatric disorders, eating disorders, sleep disorders, cardiac disorders, endocrine disorders, and pain [1-5].

In the current research and clinical practice, however, the stimulation of LVN is more frequently used than the stimulation of right vagus nerve. In adults, the right vagus nerve innervates the sinoatrial (SA) node, the atrial muscle and, to a much lesser degree, the atrioventricular (AV) node, whereas the LVN innervates the SA node and atrial muscle to a lesser degree than it innervates the AV node [6].

In many heart diseases, for instance, such as hypertension and congestive heart failure, cardiac vagal activity is diminished and unresponsive. Numerous studies have therefore proposed VNS also as a method for treating various heart conditions, including supraventricular arrhythmias, angina pectoris, congestive heart failure, atrial fibrillation, myocardial ischemia, and variant angina (spastic coronary arteries) [2, 7-12].

In addition, brain-stimulation methods are making significant inroads into psychiatric practice. For instance, VNS has also been approved as a therapy for medication-resistant depression [13, 14]. Finally, VNS is used worldwide as a non-pharmacological treatment to control seizures in epilepsy patients [15].

However, the method commonly used is in general a non-selective VNS of the LVN, which in turn causes frequent occurrence of undesirable side effects [16]. In addition, afferent VNS leads to reflex-excitation of vagal efferent activity and inhibition of sympathetic efferent activity [17].

To alleviate such problems, various selective VNS models and electrode systems that selectively stimulate certain features have been developed, that is, intermediate-diameter B-fibres in a nerve while avoiding the stimulation of A-fibres and C-fibres [18].

In the field of functional electrical stimulation, cuffs have been used as stimulation electrodes as well as electrodes for NR recording for more than 35 years [19-21]. Since then, theoretical considerations and different models have accelerated the development of cuffs [22-24].

Prosthetic devices that are controlled by cuffs recording neural activity are a reality today. However, the amount of information that can be extracted from these recordings is still limited [25].

For reliable and safe VNS using cuffs, the response of neural elements to stresses that may occur during the complex interactions taking place between the electrode and stimulated nerve, must be understood [26]. The successful long-term use of cuffs for selective VNS may be restricted by unreliable NR recordings, owing to changes in the complex impedances of the stimulating and recording electrodes [27]. Furthermore, currently there does not exist a type of multi-electrode system that adequately combines spatial and fibre-type VNS with recording.

This paper presents the hypothesis concerning both selective activation of B-fibres and selective recording of NR, wherein both stimulation and recording were performed in a single fascicle with a single-part ninety-nine-electrode cuff fitted around an isolated porcine LVN and maintained at physiological temperature in oxygenated artificial cerebrospinal fluid. One specific aim was to determine which precise parameters of the applied pulses were most effective for fibre-type VNS.

2. Methods

2.1. Model of Porcine LVN Stimulation

In this study, the porcine LVN as an animal model was adopted. The LVN is especially suited for evaluation of fiber-type selective methods because it is composed of the following three distinct groups of myelinated fibres: Aβ , hereafter referred to as A-fibres (diameter of the fibre (D ) = 5-12 μ m), myelinated Aδ -fibres, also designated as B-fibres (D = 1-5 μ m), and unmyelinated C-fibres (D = 0.5-2 μ m) [28]. The vast majority are C-fibres, whereas A- and B-fibres are in minority. A- and B-fibres in vagus nerve, however, are primarily related to innervation of cardiovascular and respiratory systems [29-38].

2.2. Multielectrode Cuff Design

The cuff was designed considering both published results of studies modelling the selective stimulation of peripheral nerves [21-23, 39] and a realistic structural topography of the porcine LVN [18]. It should be noted, however, that models of myelinated nerve fibers, used to investigate electrical nerve stimulation, are highly stylized abstractions of real fibers.

To define the relationship between the structural topography and the physical model, the distributions of fibre diameters and physical dimensions in three porcine LVNs were measured. For this purpose, approximately 8 cm long and 3 mm thick LVN was taken out from the midcervical neck of a pork (about 80 kg in weight) and fixated for histological examination. When fixated, three microtome sections were cut from each LVN at the following sites: at the middle of the LVN, at a distance of 14 mm left from the middle of the LVN, and at a distance of 14 mm right from the middle of the LVN.

A section cut from the middle of the LVN was considered as the site of stimulation origin, while the other two sections represented the sites of NR recording. For the presented paradigm, only microtome sections cut at the middle of the LVNs were analyzed. For the expression of myelin in microtome sections an immunohistochemical marker of neural tissue, namely, anti-S-100 protein antibody (PAP method), was used. Accordingly, great majority of the nerve fibers in the LVN were not expressed, since they were unmyelinated (C-fibres). Out of counted 500 myelinated fibres, 200 of them were of the diameter larger than 5 μ m (A-fibres), while 300 of them were of the diameter smaller than 5 μ m (B-fibres). Corresponding results of histological examinations and physical measurements in the LVNs, adopted in the paradigm, are shown in Table 1.

Table 1: Structural topography of the porcine LVN considered in the paradigm.

Diameter of the LVN	[approximate] 3 mm
Circumference of the LVN	[approximate] 9 mm
Type of nerve fibres	A, B, and C
Number of B-fibres (%)	300 (60%)
Average area of a single B-fibre	16 μ m2
Average diameter of a single B-fibre	[approximate] 4.6 μ m
Number of A-fibres (%)	200 (40%)
Number of C-fibres	Majority (not counted)

The cuff was produced by bonding two 0.05 mm thick silicone sheets together (BioPlexus Corporation, Ventura, CA, U.S.A.). One sheet, stretched and fixed in position, was covered with an adhesive layer (Applied Silicone Corporation, Santa Paula, CA, U.S.A.). The second, unstretched sheet was placed on top of the adhesive, and the composite was compressed to a thickness of 0.15 mm. When released, the composite curled into a spiral tube as the stretched sheet contracted to its natural length [20, 21].

Ninety-nine rectangular electrodes made of 45 μ m thick platinum ribbon (99.99% purity) were mounted on a third silicone sheet. These were arranged in a matrix of nine parallel groups, each containing eleven electrodes, and divided into the following sections (Figure 1(a)): stimulating section, containing eleven groups of three electrodes (triplet 1-11), positioned at the middle of the matrix; two blocking sections, having eleven electrodes each, positioned bilaterally next to the stimulating section; and two recording sections, having eleven spiral couples (couple 1-11), positioned next to the blocking sections.

Figure 1: (a) Matrix of ninety-nine electrodes, (b) fabricated cuff, (c) cross-section through the cuff, and (d) distances between spiral groups of electrodes and times used in calculations of the CV in the A- and B-fibres.

[figure omitted; refer to PDF]

Afterwards, the silicone sheet containing the matrix of electrodes was adhered to the inner side of the spiral tube. Figure 1 shows a fabricated cuff 44 mm in length (b) and schematic cross-section through the cuff (c).

2.3. Selective VNS

In myelinated nerve fibres, action potential (AP) is regenerated only at the nodes of Ranvier and propagate by jumping from one node to a subsequent node at rapid CV [40, 41]. One assumption made was that of an anatomical feature previously explained by Rushton [42], that is, in myelinated fibres, the internode length and the thickness of the myelin sheath are both proportional to D, and that the spike CV shows an approximately linear relationship to axon diameter throughout the entire range of myelinated fibres. Data considering the CV of different nerve fibres was adopted from the literature [43-46].

The differing properties of each individual fibre result in different AP thresholds, refractory periods, and the duration of the AP [47]. Consequently, any change in the amplitude and waveform of the NR amplitude is due to a change in the number of fibres that are firing. With increasing intensities of a stimulus, however, the number of axons firing is equivalent to the sum of all those, whose thresholds are met by a given input. The latency between the application of the stimulus and the onset of the NR is a function of the events during the depolarization and the distance between the recording site and the site of the stimulation.

A difficulty observed in previous studies of nerve stimulation is that the difference in stimulating (cathodic) intensity between the threshold excitation and the maximum recruitment of myelinated fibres is not large [48, 49]. A more selective recruitment of myelinated fibres can therefore be obtained if it is possible to exploit the difference in the threshold between the different fibre diameters. However, the conventional stimulation of nerves, using biphasic rectangular stimuli, unavoidably excites the larger nerve fibres before the smaller ones [50]. One possibility to overcome this difficulty is to introduce an exponentially decaying anodal block of nerve conduction, as described in the model of Accornero et al. [51]. Later on, it was also demonstrated by Fang and Mortimer [39] and others [52, 53] that the A-fibres could be activated and blocked at lower currents than the B-fibres. This is because the B-fibres have more nodes per unit length than the A-fibres, so they require a higher excitation voltage and thus a higher injected current to activate the AP as well as to elicit a block. With this regard it was proposed that the stimulating intensity, required for the activation AP in the A-fibres, would be slightly lower than the threshold for the B-fibre activation. Such activation would enable also the use of an intensity that is sufficiently strong to activate the B-fibres. Taking into account the above-mentioned differences in the excitation threshold and the CV of the A- and B-fibres it was assumed that the latest AP of the A-fibres arrives at an inner edge of the triplet anode at a time when the hyperpolarizing effect at the anode is still strong and should be blocked by hyperpolarization. This, in effect, produces an anodal block that prevents the propagation of the AP of the A-fibres towards the recording couple. However, the hyperpolarization should be decayed exponentially to prevent any anodal break excitation, which had previously been a problem in the stimulation with conventional rectangular pulses [39]. In addition, the intensity of the anodic phase of the stimulus had to be lower than the threshold of the anodic activation of the A-fibres in order to prevent the excitation due to the reversed polarity on the cathode during the anodic phase. According to the paradigm, an exponentially decaying part of cathodic phase of the stimulus should last until the latest AP of the A-fibres passes the outer edge of the blocking electrode. The AP in the B-fibres, however, should reach an inner edge of the triplet anode after the cathodic intensity has exponentially decayed to a point that is inefficient at blocking their conduction. The selectivity of the technique for fibres of different CV could therefore only be achieved by changing the time constant of the exponential decay and by adjusting the stimulating cathodic intensity.

The resulting stimulus shown in Figure 2 was a current, biphasic, charge-balanced, and asymmetric pulse, composed of a precisely determined quasitrapezoidal cathodic phase with a square leading edge of intensity _ic , a plateau of the cathodic phase with the width of _tc , followed by an exponentially decaying phase with the width of _texp and the time constant _τexp , and ended by a wide, rectangular, anodic phase with the width _ta and intensity of _ia .

Figure 2: Parameters and waveform of the stimulus determined in the paradigm.

[figure omitted; refer to PDF]

As a result of the cuff design, the following dimensions of the cuff were adopted:

: nominal nerve diameter: 2.5 mm,

: cuff length: 44 mm,

: electrode length: 2 mm,

: electrode width: 0.5 mm,

: circumferential separation between electrodes: 0.5 mm,

: longitudinal separation between electrodes: 2 mm.

It was presumed that the A- and B-fibres could be activated at any site on the cathode and the corresponding AP would propagate simultaneously in both directions. The population of closely spaced nerve fibres above an appointed triplet with roughly synchronous firing pattern was considered as a single pathway.

It was also assumed in the paradigm that an additional blocking electrode situated next to the triplet anode would provide an additional blocking effect of the A-fibres. In the two blocking sections of eleven electrodes, each one was galvanically connected to the neighboring triplet anode.

The efficacy of the selective VNS and recording of the NR are strongly dependent on the physical proximity of the deployed electrodes and the distance between the electrodes and the nerve fibres. To stimulate a certain group of fibres in a fascicle, while also avoiding causing injury associated with high charge density, a well-defined electrical charge should be applied to preselected locations [54]. To keep the electrode-electrolyte interface within capacitive mechanisms, the cathode was dimensioned so as not to exceed the limits for a reversible charge injection [55].

In the proposed paradigm, the view of Sunderland and Bedbrook [56], wherein fibres in fascicles extending peripherally along a nerve extensively interweave so that they converge and diverge into new and different fascicular assemblies, was adopted. For this reason, the distance between the stimulating cathode and the couple within the row of nine electrodes was fixed to be well below that at which the confluence of fibres significantly changes (Figures 1(a) and 1(b)).

Only myelinated A-fibres (CV = 30-70 m/s) and myelinated B-fibres (CV = 3-15 m/s) were considered in the paradigm. The unmyelinated C-fibres (CV = 0.5-2 m/s) however, could not be activated by the proposed stimuli, so they were not considered.

Taking into account the differences in the CV of the A- and B-fibres that could potentially contribute to the NR during the selective VNS, Table 2, describing the propagation of AP in the A- and B-fibres from the triplet cathodes to different groups of electrodes within the cuff, was compiled (Figure 1(d)). For this purpose, a previously reported calculation method, using length and latency differences between stimulated and recording sites on a LVN, was employed [57, 58].

Table 2: Timing of AP propagation in A- and B-fibres between longitudinal electrodes.

Fibre type	t 1 (μ s)	t 2 (μ s)	t 3 (μ s)	t 4 (μ s)
Max	Min	Max	Min	Max	Min	Max	Min
A	133.3	57.1	266.7	114.3	400.0	171.4	533.3	228.6
B	1333.3	266.7	2666.7	533.3	4000.0	800.0	5333.3	1066.7
Δ t (B mean-A mean)	Δ t 1 = 704.8	Δ t 2 = 1992.9	Δ t 3 = 2114.3	Δ t 4 = 3200.0

_{t 1} : time to reach an inner edge of a triplet anode;

_{t 2} : time to reach an inner edge of a blocking electrode;

_{t 3} : time to reach an inner edge of a positive recording electrode;

_{t 4} : time to reach an inner edge of a negative recording electrode.

2.4. Selective NR Recording

The fact that different contacts in a cuff preferentially record different pathways in the nerve, demonstrated by numerous authors, was considered [23, 59-62].

It was obvious that fibre distribution, electrode separation in the couples, and the distance between the active fibres and the recording electrode would all have major roles in determining the peak-peak NR amplitude [63].

According to the paradigm, positive NR deflection was expected to occur at the moment when the eventually nonblocked A-fibre and particularly fastest B-fibre AP reached the inner edge of a positive recording electrode, lasting until the last of the B-fibre AP passed the outer edge of a positive recording electrode. However, negative NR deflection was expected to occur at the moment when the eventually nonblocked A-fibre and fastest B-fibre AP reached the inner edge of a negative recording electrode, and lasting until the last of the B-fibre AP passed the outer edge of a negative recording electrode.

2.5. Preparation of Nerves

Experiments were performed in four Slovenian male Landrace pigs weighing about 150 kg each. The animals were killed according to established protocols for mature animals, sows, or boars (Penetrating Captive Bolt and immediate discharge of blood). The obtained neural tissue was treated in accordance with the approval provided by the ethics committee at the Veterinary Administration, Ministry of Agriculture, Forestry and Food, Republic of Slovenia (number: 34401-27/2010/3).

The LVN was removed from the approximately 10 cm long porcine midcervical neck, carefully freed from excessive fat tissue and shortened to about 8 cm (Figure 3(a)). Afterwards, an experimenter spread the cuff, placed the LVN into the cuff, and subsequently wetted the LVN with cotton wool dipped into an artificial cerebrospinal fluid (comprising (in mM): MgCl₂ 2, CaCl₂ 2, KCl 2.5, NaCl 126, glucose 10, NaH₂ PO₄ ·H₂ O 1.25, and NaHCO₃ 26), oxygenated at room temperature with a mixture of oxygen and carbon dioxide (95% O₂ /5% CO₂ ). The cuff was snugly fitted to the wet LVN, covering the entire nerve perimeter, and closed. Finally, the resulting composition was placed into the experimental chamber (Figure 3(b)), which was heated to 37°C using a precision water circulator (Perfectherm PFV, Boehringer, Labor Manheim GmbH für Labortechnik, Germany).

Figure 3: (a) Isolated LVN and (b) LVN within the cuff mounted into an experimental chamber.

[figure omitted; refer to PDF]

2.6. Selective VNS Procedure

LVN was stimulated with a custom-designed stimulator, delivering single stimuli (1 Hz) to the appointed triplet 5. The experiments carried out on each LVN consisted of four chronological subexperiments, referred to as Phases 1-4, where the stimulus preset at the stimulator was recorded simultaneously with the elicited NR. At the beginning of Phases 1-4 the following parameters were fixed: _ic =4 mA, _tc =155 μ s, _texp =100 μ s, _τexp =30 μ s, _ia =0.45 mA, and _ta =490 μ s.

In Phase 1, variable parameter was _ic (range between 0.8 and 4.2 mA); in Phase 2, variable parameter was _tc (range between 60 and 300 μ s); and in Phase 3, variable parameter was _τexp (range between 20 and to 60 μ s); while other parameters remained fixed. In Phase 4, however, variable parameter was the site of NR recording (couple 1-11) (Figure 1(c)) relative to the appointed triplet 5, while stimulating parameters were not of importance because only spatial selectivity was tested in this phase.

Table 3 shows numerical values of the parameters and waveforms that resulted in the most indicative alteration of NR recorded simultaneously in each phase. Regarding the preset stimuli delivered in all four phases to the triplet 5, Table 3 shows the values of the charge _Qc injected in the cathodic phase, as well as values of the charge _Qa injected by a triplet anode in the anodic phase of the stimuli. The triplet anodes and neighboring blocking electrodes were galvanically connected so that the charge density at the anode would be approximately four times lower than that at the cathode. To test the influences of different stimuli on offsets in NR that might be elicited during the selective VNS, the corresponding degrees of imbalance between _Qc and _Qa were calculated and compared.

Table 3: Parameters of the stimuli and values of recorded NR.

Variable	Example (a)	Example (b)	Example (c)	Example (d)
i c (mA)	3.10	4.00	4.00	4.00
t c (μ s)	155	155	155	265
t exp (μ s)	100	100	100	100
t au (μ s)	30	30	45	30
i a (mA)	-0.4	-0.45	-0.45	-0.45
t a (μ s)	490	490	490	490
Q c (nAs)	539	770	800	1100
Q a (nAs)	-163	-208	-208	-200
Integral 1 (nVs)	512.20	665.22	536.22	n/a
NRpeak (mV)	2.76	3.17	2.5	n/a
NR w (μ s)	190	210	222	n/a

For the assessment of the population of fibres that most probably contributed to the NR during selective VNS and, consequently, to extract the relevant stimulation and block parameters, patterns recorded in above-mentioned four phases were analyzed.

2.7. NR Measurement and Analysis

Unidirectional NR measured with couples 5 or 9 was analysed. The measured NR and voltage drops on the precision serial resistor at the stimulator output were amplified at a differential amplifier (A = 100 for the neural response and at A = 10 for _ic ). Afterwards, the data for both signals were gathered at 200 kHz using a high-performance data-acquisition system (DEWE-43, DEWESOFT d. o. o., Republic of Slovenia) and proprietary acquisition software (DEWESoft 7.0.2); the data was then stored on a Lenovo T61 portable computer (Lenovo, Singapore). Offline signal analysis was performed using MATLAB R2007a software (The Mathworks Inc., USA).

To test the paradigm, the components in the NR that had potentially originated from the A- and B-fibres needed to be identified. For this purpose, the integral of the NR observed under the cathodic phase of the stimulus was calculated. This integral actually represented the cumulative contribution of the AP originating from both the A- and B-fibres as well as an ensemble artefact.

To test both, spatial and fibre-type selectivity, the LVN was stimulated with triplet 5, while in Phase 4 the resulting NR were measured from couple 9, actually situated at a site opposite to couple 5 (Figure 1(c)).

In addition, the times considered needed for the AP in the A- and B-fibres to propagate between the different groups of electrodes (Figure 1(d)) are indicated in Table 2. For this purpose, the time separations between the AP of A and B-fibres to reach an inner edge of an anode (Δ_t1 ), to reach an inner edge of a blocking electrode (Δ_t2 ), to reach an inner edge of a positive couple electrode (Δ_t3 ), and to reach a negative couple electrode (Δ_t4 ) were calculated.

In this stage of the study, the data obtained from the different animals were variable and so could not be effectively combined. Therefore, data from only one animal are presented.

3. Results and Discussion

The most indicative alterations obtained in each Phase 1-4, were extracted and presented as the examples of interest (a)-(d). Accordingly, Figure 4 shows examples (a)-(d) comprising preset stimulus waveforms and the proprietary NR.

Examples of preset stimulus waveforms and the proprietary NR. Solid line represents stimulating intensity and dashed line represents neural response NR. NR_peak is obtained at the top of the bell-shaped NR in the cathodic phase of the stimulus; N_Rw is made at 50% of NR_peak .

(a) [figure omitted; refer to PDF]

(b) [figure omitted; refer to PDF]

(c) [figure omitted; refer to PDF]

(d) [figure omitted; refer to PDF]

Example (d), wherein the NR was measured using couple 9, is shown as a reference used for the validity of the paradigm. It was presumed that if the paradigm was valid, then the corresponding NR should not contain the AP of fibres activated with triplet 5. In other words, the integral of the NR manifested under the cathodic phase of the stimulus, representing the cumulative contribution of the AP to the NR, could be interpreted exclusively as an ensemble artefact.

In examples (a)-(c) two components of interest, the amplitude of cathodic peak of the NR (NR_peak ) and the width of the cathodic NR (N_Rw ), were identified. NR_peak was obtained at the top of the bell-shaped NR in the cathodic phase of the stimulus, whereas _NRw was made at 50% of NR_peak . With regard to the NR, Table 3 also shows the calculated values of integral 1, manifested under a part of NR belonging to the cathodic phase of the stimulus. These values quantitatively show how the NR changed with different pre-set parameters of the stimuli during selective VNS of the LVN.

From an electrochemical point of view, electrochemical reactions that occured at the cathode owing to the charge _Qc injected via the cathodic phase were reversed only in part by the charge _Qa , injected via the anodic phase. Since the triplet anodes and neighboring blocking electrodes were galvanically connected, the charge densities _Qa presented in Table 3 correspond to the charge applied to a single anode within the tested triplet. In examples (a)-(d) (see Figure 4), however, the measured _ia corresponds to the sum of currents delivered to each of four galvanically connected electrodes.

By comparison the examples in Table 3, the largest change in the bell-shaped NR was identified. It was observed in example (a) that _ic influenced both the NR_peak as well as _NRw . Precisely, lower _ic elicited slightly lower NR_peak and slightly smaller _NRw . Therefore, it is presumed that change of _ic influenced both A- and B-fibres. However, the axons had differing CV and additional numbers of axons firing contributed to the width of the bell-shaped NR, and not only to the peak amplitude. Therefore, _NRw was indicative of the additional contribution of the B-fibres to the bell-shaped NR when specific parameters and waveforms of the stimulus were used.

It was presumed that the B-fibres contributed the AP that fell towards the region between the apex and the tail section of the bell-shaped NR and that the A-fibres, if not completely blocked, contributed the AP that fell towards the beginning.

It was noticed that parameter _tc did not influence the bell-shaped NR significantly.

One parameter obviously the most important for the confirmation of the paradigm was _τexp . It could be seen in example (c) when compared to example (b) that larger _τexp (45 μ s) exhibited significant influence on _NRw . It could be assumed that the AP of the A-fibres was more effectively blocked, whereas the AP of the B-fibres passed through. A certain lack of AP from A-fibres however could also be confirmed by the lower NR_peak .

It could be noticed that integrals 1 belonging to examples (a) and (c) were significantly smaller compared to integral 1 in example (b). With this regard, it could be speculated that in example (a) less A-fibres were activated, while in example (c) more of them were blocked. It could also be noticed that integral 1 in example (c) is larger when compared to example (a) and is probably a result of easier passage of AP in B-fibres using preset _τexp (shown as larger _NRw ), while the contribution of A fibres was smaller (shown as lower NR_peak ).

Regarding the offset of the NR, it could be seen that in all four examples (a)-(d) the bell-shaped NR was not greatly influenced by imbalances between _Qc and _Qa .

From initial analysis it appeared that the hypothesis stated in the paradigm was shown to be loosely true. Namely, there was a large stimulation artefact superimposed in the NR that greatly obscured the components of the NR, and various components did overlap. These were considered as an ensemble artefact, where the strongest components were thought to arise from the stimulating pulse via the transient response characteristics of an electrode/neural tissue interface and in part from the inherent capacitance of the LVN. An actual ratio between the contribution of the NR and the contribution of the ensemble artefact, relating the waveform of the pulse to the recorded signal, was not defined within this study. In addition, the recordings did not show separate peaks corresponding to the A- and B-fibre types. However, if the distance between the stimulating and recording spiral section was larger, the accuracy of NR measurements would certainly improve [63].

Similar results regarding the selectivity and activation thresholds of the A- and B-fibres were obtained in recent "in vivo " studies [64, 65]. Differences in the selective activating and blocking efficacies obtained could be attributed to the different shapes and dimensions of the electrodes used in these studies. If successfully developed further, this paradigm would enable significant improvement of current neuroprostheses. Namely, for neuroprosthetic applications, the usefulness of the paradigm as a framework for determining the activation of combinations of several pathways could certainly be contingent on future developments. Furthermore, a preference towards nerve-fibre type stimulation is advantageous in applications where organ-specific stimulation is required and the side-effect profile related to the propagation of the AP of the A-fibres towards the CNS needs to be minimized. However, more detailed investigations of neural control systems, considering the realistic structural topography of the nerve and the presence of a spatiotemporal constraint based on the electrophysiology of myelinated nerve fibres, should be carried out based on this work. Ultimately, selective NR recording from LVN fibres could be effectively used for closed-loop control of implantable stimulators selectively activating different neural pathways.

4. Conclusions

The major finding of the tests carried out on the proposed paradigm, wherein stimulating pulses were applied to preselected locations along an insulated porcine LVN via an appointed stimulating triplet, was the activation of AP in the A- and B-fibres within the corresponding pathways and the slight inhibition of AP in the A-fibres. Activation and inhibition were noticed from the widths and amplitudes of the measured NR to the VNS (bell-shaped NR). Namely, it can be seen in Figure 4 and Table 3 that the bell-shaped NR slightly changed as different parameters and waveforms of the stimulus were chosen. Despite the above-mentioned drawbacks of the paradigm approach used in this study, we achieved an encouraging correlation between the expected and measured NR. In this regard, the directions our further work could take are the following:

(i) subtraction of the artifacts modifying the components in the NR;

(ii) further development of the cuff;

(iii): testing the ability of the cuff to steer _ic in a desired direction within the LVN;

(iv) testing the ability of the cuff to perform selective VNS in an orthodromic or an antidromic direction in the LVN;

(v) an independent measurement of compound action potential (CAP) using external hook electrodes;

(vi) comparison of NR measured by the cuff to CAP measured by external hook electrodes;

(vii): performing more experiments to show group statistics.

Acknowledgments

This work was funded by the Slovenian Research Agency, research Programme P3-0171. The authors thank Mr. Marjan Kavcic s.p., Drenov Gric 9, Vrhnika, Slovenia, for providing the vagus nerve examples. The research was performed also as partial fulfillment of the requirements for the Doctor of Philosophy degree (Polona Peclin, MD).

Conflict of Interests

The authors state that there is no conflicts of interests related to the results presented within the paper. None of the authors received any financial or other compensation as a result of the study or during its execution. No personal relationships with other people or organizations inappropriately influenced the work.

References

[1] A. Carlsten, B. Folkow, C. A. Hamberger, "Cardiovascular effects of direct vagal stimulation in man," Acta physiologica Scandinavica , vol. 41, no. 1, pp. 68-76, 1957.

[2] A. M. Bilgutay, I. M. Bilgutay, F. K. Merkel, C. W. Lillehei, "Vagal tuning. A new concept in the treatment of supraventricular arrhythmias, angina pectoris, and heart failure," Journal of Thoracic and Cardiovascular Surgery , vol. 56, no. 1, pp. 71-82, 1968.

[3] N. P. Xenopoulos, R. J. Applegate, "The effect of vagal stimulation on left ventricular systolic and diastolic performance," American Journal of Physiology--Heart and Circulatory Physiology , vol. 266, no. 6, part 2, pp. H2167-H2173, 1994.

[4] H.-R. Berthoud, W. L. Neuhuber, "Functional and chemical anatomy of the afferent vagal system," Autonomic Neuroscience , vol. 85, no. 1-3, pp. 1-17, 2000.

[5] M. N. Levy, P. Martin, W. C. Randall, "Parasympathetic control of the heart," Nervous Control of Cardiovascular Function , pp. 68-94, Oxford University Press, New York, NY, USA, 1984.

[6] R. E. Klabunde Cardiovascular Physiology Concepts , Wolters Kluwer Health, Lippincott Williams & Wilkins, Philadelphia, PA, USA, 2012.

[7] P. Parker, B. G. Celler, E. K. Potter, D. I. McCloskey, "Vagal stimulation and cardiac slowing," Journal of the Autonomic Nervous System , vol. 11, no. 2, pp. 226-231, 1984.

[8] J. A. Armour, J. A. Armour, J. L. Ardell, "Peripheral autonomic neuronal interactions in cardiac regulation," Neurocardiology , pp. 219-249, Oxford University Press, New York, NY, USA, 1994.

[9] A. V. Zamotrinsky, B. Kondratiev, J. W. de Jong, "Vagal neurostimulation in patients with coronary artery disease," Autonomic Neuroscience , vol. 88, no. 1-2, pp. 109-116, 2001.

[10] Y. Zhang, H. Yamada, S. Bibevski, S. Zhuang, K. A. Mowrey, D. W. Wallick, S. Oh, T. N. Mazgalev, "Chronic atrioventricular nodal vagal stimulation: first evidence for long-term ventricular rate control in canine atrial fibrillation model," Circulation , vol. 112, no. 19, pp. 2904-2911, 2005.

[11] T. A. Anholt, S. Ayal, J. A. Goldberg, "Recruitment and blocking properties of the CardioFit stimulation lead," Journal of Neural Engineering , vol. 8, no. 3, 2011.

[12] V. Pavlov, G. Spassov, T. Mazgalev, S. Mileva, "Effects of vagal stimuli of different duration on heart rate and atrioventricular conduction," Acta Physiologica et Pharmacologica Bulgarica , vol. 2, no. 4, pp. 21-26, 1976.

[13] C. Frick, D. Verstappen, T. E. Schlaepfer, "Time course of effects of vagus nerve stimulation on Hamilton depression rating scale symptom levels," Biological Psychiatry , vol. 61, no. 8, pp. 233S-233S, 2007.

[14] D. M. Labiner, G. L. Ahern, "Vagus nerve stimulation therapy in depression and epilepsy: therapeutic parameter settings," Acta Neurologica Scandinavica , vol. 115, no. 1, pp. 23-33, 2007.

[15] C. Heck, S. L. Helmers, C. M. DeGiorgio, "Vagus nerve stimulation therapy, epilepsy, and device parameters: scientific basis and recommendations for use," Neurology , vol. 59, no. 6, supplement 4, pp. S31-S37, 2002.

[16] I. I. Ali, N. A. Pirzada, Y. Kanjwal, B. Wannamaker, A. Medhkour, M. T. Koltz, B. V. Vaughn, "Complete heart block with ventricular asystole during left vagus nerve stimulation for epilepsy," Epilepsy and Behavior , vol. 5, no. 5, pp. 768-771, 2004.

[17] P. J. Schwartz, "Vagal stimulation for heart diseases: from animals to men--an example of translational cardiology," Circulation Journal , vol. 75, no. 1, pp. 20-27, 2011.

[18] C. R. Butson, I. O. Miller, R. A. Normann, G. A. Clark, "Selective neural activation in a histologically derived model of peripheral nerve," Journal of Neural Engineering , vol. 8, no. 3, 2011.

[19] R. B. Stein, D. Charles, L. Davis, "Principles underlying new methods for chronic neural recording," Canadian Journal of Neurological Sciences , vol. 2, no. 3, pp. 235-244, 1975.

[20] G. G. Naples, J. D. Sweeny, J. T. Mortimer, "Implantable cuff, Method and manufacture, and Method of installation," Inventors, U. S. Patent #4, 602, 624, 1986

[21] J. D. Sweeney, D. A. Ksienski, J. T. Mortimer, "A nerve cuff technique for selective excitation of peripheral nerve trunk regions," IEEE Transactions on Biomedical Engineering , vol. 37, no. 7, pp. 706-715, 1990.

[22] E. V. Goodall, L. M. Kosterman, J. Holsheimer, J. J. Struijk, "Modeling study of activation and propagation delays during stimulation of peripheral nerve fibers with a tripolar cuff electrode," IEEE Transactions on Rehabilitation Engineering , vol. 3, no. 3, pp. 272-282, 1995.

[23] M. Sahin, D. M. Durand, "Selective recording with a multi-contact nerve cuff electrode," in Proceedings of the 18th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBS '96), pp. 369-370, October-November 1996.

[24] A. Q. Choi, J. K. Cavanaugh, D. M. Durand, "Selectivity of multiple-contact nerve cuff electrodes: a simulation analysis," IEEE Transactions on Biomedical Engineering , vol. 48, no. 2, pp. 165-172, 2001.

[25] R. Rieger, J. Taylor, E. Comi, N. Donaldson, M. Russold, C. M. O. Mahony, J. A. McLaughlin, E. McAdams, A. Demosthenous, J. C. Jarvis, "Experimental determination of compound action potential direction and propagation velocity from multi-electrode nerve cuffs," Medical Engineering and Physics , vol. 26, no. 6, pp. 531-534, 2004.

[26] H. H. Hoffman, A. Kuntz, "Vagus nerve components," The Anatomical Record , vol. 127, no. 3, pp. 551-567, 1957.

[27] J. C. Williams, J. A. Hippensteel, J. Dilgen, W. Shain, D. R. Kipke, "Complex impedance spectroscopy for monitoring tissue responses to inserted neural implants," Journal of Neural Engineering , vol. 4, no. 4, pp. 410-423, 2007.

[28] H. S. Gasser, "The classification of nerve fibers," Ohio Journal of Science , vol. 41, no. 3, pp. 145-159, 1941.

[29] E. Agostoni, J. E. Chinnock, M. de Burgh Daly, J. G. MURRAY, "Functional and histological studies of the vagus nerve and its branches to the heart, lungs and abdominal viscera in the cat," The Journal of Physiology , vol. 135, no. 1, pp. 182-205, 1957.

[30] H. N. Schnitzlein, L. C. Rowe, H. H. Hoffman, "The myelinated component of the vagus nerve in man," Anatomical Record , vol. 131, pp. 649-666, 1958.

[31] R. M. McAllen, K. M. Spyer, "Two types of vagal preganglionic motoneurones projecting to the heart and lungs," Journal of Physiology , vol. 282, pp. 353-364, 1978.

[32] L. Rossi, M. N. Levy, "Histology of cardiac vagal innervation in man," Vagal Control of the Heart: Experimental Basis and Clinical Implications , pp. 3-20, Futura Publishing, Armonk, NY, USA, 1994.

[33] K. M. Spyer, P. A. Brooks, P. N. Izzo, M. N. Levy, "Vagal preganglionic neurons supplying the heart," Vagal Control of the Heart: Experimental Basis and Clinical Implications , pp. 45-64, Futura Publishing, Armonk, NY, USA, 1994.

[34] S. Nosaka, K. Yasunaga, S. Tamai, "Vagal cardiac preganglionic neurons: distribution, cell types, and reflex discharges," American Journal of Physiology--Regulatory Integrative and Comparative Physiology , vol. 12, no. 1, pp. R92-R98, 1982.

[35] J. F. X. Jones, Y. Wang, D. Jordan, "Dorsal medullary neurones activated by stimulation of the cardiac vagal branch of the anaesthetized rat, and their behaviour during the pulmonary chemoreflex," Journal of Physiology , vol. 483, pp. 89-90, 1995.

[36] J. F. X. Jones, Y. Wang, D. Jordan, "Activity of C fibre cardiac vagal efferents in anaesthetized cats and rats," Journal of Physiology , vol. 507, no. 3, pp. 869-880, 1998.

[37] Z. Cheng, T. L. Powley, "Nucleus ambiguus projections to cardiac ganglia of rat atria: an anterograde tracing study," Journal of Comparative Neurology , vol. 424, no. 4, pp. 588-606, 2000.

[38] H. H. Hoffman, H. N. Schnitzlein, "The numbers of nerve fibers in the vagus nerve of man," The Anatomical Record , vol. 139, pp. 429-435, 1961.

[39] Z.-P. Fang, J. T. Mortimer, "Selective activation of small motor axons by quasitrapezoidal current pulses," IEEE Transactions on Biomedical Engineering , vol. 38, no. 2, pp. 168-174, 1991.

[40] R. FitzHugh, "Computation of impulse initiation and saltatory conduction in a myelinated nerve fiber," Biophysical Journal , vol. 2, no. 1, pp. 11-21, 1962.

[41] R. Stampfli, "Saltatory conduction in nerve," Physiological Reviews , vol. 34, no. 1, pp. 101-112, 1954.

[42] W. A. H. Rushton, "A theory of the effects of fibre size in medullated nerve," Journal of Physiology , vol. 115, no. 1, pp. 101-122, 1951.

[43] H. S. Gasser, J. Erlanger, "The role played by the sizes of the constituent fibers of a nerve-trunk in determining the form of its action potencial wave," American Journal of Physiology , vol. 80, pp. 522-547, 1927.

[44] D. F. Stegeman, J. P. C. de Weerd, "Modelling compound action potentials of peripheral nerves in situ. I. Model description: evidence for a non-linear relation between fibre diameter and velocity," Electroencephalography and Clinical Neurophysiology , vol. 54, no. 4, pp. 436-448, 1982.

[45] E. R. Kandel, J. F. Schwartz, T. M. Jessell Principles of Neural Science , Elsevier Science, New York, NY, USA, 1991., 3rd.

[46] G. M. Manzano, L. M. P. Giuliano, J. A. M. Nóbrega, "A brief historical note on the classification of nerve fibers," Arquivos de Neuro-Psiquiatria , vol. 66, no. 1, pp. 117-119, 2008.

[47] J. G. Whitwam, "Classification of peripheral nerve fibres. An historical perspective," Anaesthesia , vol. 31, no. 4, pp. 494-503, 1976.

[48] D. R. McNeal, D. A. Teicher, F. T. Hambrecht, J. B. Reswick, "Effect of electrode placement on threshold and initial site of excitation of a myelinated nerve fibre," Functional Electrical Stimulation , pp. 405-412, Marcel Dekker, New York, NY, USA, 1977.

[49] C. van den Honert, J. T. Mortimer, "The response of the myelinated nerve fiber to short duration biphasic stimulating currents," Annals of Biomedical Engineering , vol. 7, no. 2, pp. 117-125, 1979.

[50] P. H. Gorman, J. T. Mortimer, "The effect of stimulus parameters on the recruitment characteristics of direct nerve stimulation," IEEE Transactions on Biomedical Engineering , vol. 30, no. 7, pp. 407-414, 1983.

[51] N. Accornero, G. Bini, G. L. Lenzi, M. Manfredi, "Selective activation of peripheral nerve fiber groups of different diameter by triangular shaped stimulus pulses," Journal of Physiology , vol. 273, no. 3, pp. 539-560, 1977.

[52] K. Fukushima, O. Yahara, M. Kato, "Differential blocking of motor fibers by direct current," Pflügers Archiv European Journal of Physiology , vol. 358, no. 3, pp. 235-242, 1975.

[53] D. R. McNeal, "Analysis of a model for excitation of myelinated nerve," IEEE Transactions on Biomedical Engineering , vol. 23, no. 4, pp. 329-337, 1976.

[54] W. F. Agnew, D. B. McCreery, "Considerations for safety with chronically implanted nerve electrodes," Epilepsia , vol. 31, no. 2, pp. S27-S32, 1990.

[55] B. Onaral, H. H. Sun, H. P. Schwan, "ELECTRICAL PROPERTIES OF BIOELECTRODES," IEEE Transactions on Biomedical Engineering , vol. 31, no. 12, pp. 827-832, 1984.

[56] S. Sunderland, G. M. Bedbrook, "The cross-sectional area of peripheral nerve trunks occupied by the fibres representing individual muscular and cutaneous branches," Brain , vol. 72, no. 4, pp. 613-624, 1949.

[57] D. Gu, R. E. Gander, E. C. Crichlow, "Determination of nerve conduction velocity distribution from sampled compound action potential signals," IEEE Transactions on Biomedical Engineering , vol. 43, no. 8, pp. 829-838, 1996.

[58] W. H. Olsen, S. L. BeMent, L. J. Dorfman, K. L. Cummins, L. J. Leifer, "Compound action potential reconstructions and predicted fiber diameter distributions," In Conduction Velocity Distributions , pp. 57-83, Alan R. Liss, New York, NY, USA, 1981.

[59] J. J. Struijk, M. K. Haugland, M. Thomsen, "Fascicle selective recording with a nerve cuff electrode," in Proceedings of the 18th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBS '96), pp. 361-362, October-November 1996.

[60] P. R. Christensen, Y. Chen, K. D. Strange, J. A. Hoffer, "Multichannel recordings from peripheral nerves: 3. Evaluation of selectivity using mechanical stimulation of individual digits," in Proceedings 2nd Annual Conference of the International Functional Electrical Stimulation Society (IFESS '97), pp. 110, August 1997.

[61] J. Zariffa, M. K. Nagai, M. Schuettler, T. Stieglitz, Z. J. Daskalakis, M. R. Popovic, "Use of an experimentally derived leadfield in the peripheral nerve pathway discrimination problem," IEEE Transactions on Neural Systems and Rehabilitation Engineering , vol. 19, no. 2, pp. 147-156, 2011.

[62] B. Wodlinger, D. M. Durand, "Selective recovery of fascicular activity in peripheral nerves," Journal of Neural Engineering , vol. 8, no. 5, 2011.

[63] I. F. Triantis, A. Demosthenous, N. Donaldson, "On cuff imbalance and tripolar ENG amplifier configurations," IEEE Transactions on Biomedical Engineering , vol. 52, no. 2, pp. 314-320, 2005.

[64] A. Vuckovic, M. Tosato, J. J. Struijk, "A comparative study of three techniques for diameter selective fiber activation in the vagal nerve: anodal block, depolarizing prepulses and slowly rising pulses," Journal of Neural Engineering , vol. 5, no. 3, pp. 275-286, 2008.

[65] S. C. M. A. Ordelman, L. Kornet, R. Cornelussen, H. P. J. Buschman, P. H. Veltink, "An indirect component in the evoked compound action potential of the vagal nerve," Journal of Neural Engineering , vol. 7, no. 6, 2010.