Force sensor in simulated skin and neural model mimic tactile SAI afferent spiking response to ramp and hold stimuli
© Kim et al.; licensee BioMed Central Ltd. 2012
Received: 3 October 2011
Accepted: 5 July 2012
Published: 23 July 2012
The next generation of prosthetic limbs will restore sensory feedback to the nervous system by mimicking how skin mechanoreceptors, innervated by afferents, produce trains of action potentials in response to compressive stimuli. Prior work has addressed building sensors within skin substitutes for robotics, modeling skin mechanics and neural dynamics of mechanotransduction, and predicting response timing of action potentials for vibration. The effort here is unique because it accounts for skin elasticity by measuring force within simulated skin, utilizes few free model parameters for parsimony, and separates parameter fitting and model validation. Additionally, the ramp-and-hold, sustained stimuli used in this work capture the essential features of the everyday task of contacting and holding an object.
This systems integration effort computationally replicates the neural firing behavior for a slowly adapting type I (SAI) afferent in its temporally varying response to both intensity and rate of indentation force by combining a physical force sensor, housed in a skin-like substrate, with a mathematical model of neuronal spiking, the leaky integrate-and-fire. Comparison experiments were then conducted using ramp-and-hold stimuli on both the spiking-sensor model and mouse SAI afferents. The model parameters were iteratively fit against recorded SAI interspike intervals (ISI) before validating the model to assess its performance.
Model-predicted spike firing compares favorably with that observed for single SAI afferents. As indentation magnitude increases (1.2, 1.3, to 1.4 mm), mean ISI decreases from 98.81 ± 24.73, 54.52 ± 6.94, to 41.11 ± 6.11 ms. Moreover, as rate of ramp-up increases, ISI during ramp-up decreases from 21.85 ± 5.33, 19.98 ± 3.10, to 15.42 ± 2.41 ms. Considering first spikes, the predicted latencies exhibited a decreasing trend as stimulus rate increased, as is observed in afferent recordings. Finally, the SAI afferent’s characteristic response of producing irregular ISIs is shown to be controllable via manipulating the output filtering from the sensor or adding stochastic noise.
This integrated engineering approach extends prior works focused upon neural dynamics and vibration. Future efforts will perfect measures of performance, such as first spike latency and irregular ISIs, and link the generation of characteristic features within trains of action potentials with current pulse waveforms that stimulate single action potentials at the peripheral afferent.
KeywordsTactile Force sensor Leaky integrate-and-fire Mechanoreceptor SAI Neural Prosthetic Biomechanics Electrophysiology Skin Elasticity
Our sense of touch helps us perform activities of daily living, such as grasping a glass, discerning the structure of a coin, and buttoning a shirt. Completing these tasks proves difficult for the 541,000 U.S. citizens living with upper limb loss. In the near future, however, advanced prosthetics may help reconstitute their motor and sensory function[2–4]. It is likely that sensory connectivity to the nervous system, via either peripheral afferent or cortex, will originate from mechanical sensors within an elastic skin. One major hurdle lies in modulating the delivery of the appropriate signals to the afferent, specifically both the shape of current pulses to elicit action potentials and the features within trains of action potentials. Work herein is focused upon the latter, in particular, those essential features captured by the slowly adapting type I (SAI) afferent in its response to the everyday task of contacting and holding an object. The intent is to complement and extend prior work with vibratory stimuli.
Among mechanoreceptive afferent types, SAI afferents respond to sustained compression as well as movement. SAI afferents exhibit firing rates about 10 times greater during stimulus movement than sustained hold, whereas their firing rates during held stimuli increases linearly with indentation depth from 0.1 to 2.0 mm. SAI afferents exhibit comparatively high variability in the time interval between spikes.
Although both skin mechanoreceptors and artificial pressure sensors respond to compression, their outputs are quite different. Mechanoreceptive afferents in the skin produce trains of discrete action potentials, or “spikes,” whose rates represent salient stimulus features such as magnitude or velocity. In contrast, artificial force sensors output a continuous analog signal, typically with changes in resistance or capacitance that represent intensity[8, 9]. The rate of intensity change is decoded from signal slope, rather than spike frequency. Other sensor technologies, including optical fibers and tracking, output similar analog signals. Another difference between biological and artificial sensors is that mechanoreceptors respond to stress and/or strain at the location of natural tactile end organs, which lie embedded within the skin. Mechanics models clearly indicate that the skin influences the propagation of forces from its surface to the end organs, shaping the resultant neural response[12–16].
Work herein seeks to address these gaps and contribute to the knowledge base by computationally replicating how the slowly adapting type I (SAI) mechanoreceptor in its skin environment converts ramp-and-hold, sustained stimuli into spike trains that capture the features of contacting and holding an object. These features include the SAI’s increased neural firing to changes in stimulus magnitude and rate of change, the SAI afferent’s continuous firing response to sustained stimuli and greater response to moving stimuli. To address these gaps, this systems integration effort builds upon research that includes models of membrane transduction and neural dynamics by integrating a physical force sensor-elastic substrate with a mathematical model of neuronal spiking to transform sensor output to spike trains that capture the SAI response’s characteristic features. Model predictions are compared to electrophysiological recordings from mouse SAI afferents for similar ramp-and-hold stimuli. The rigorous engineering technique of response surface methodology was used to optimize the model’s six free parameters by following a gradient method of steepest descent to match predicted spike response with those observed in the SAI afferent, before separately validating the model. In addition, the work begins to consider the SAI afferent’s variable interspike interval and first spike latency.
Spiking-sensor model: force sensor-elastic substrate
When applied load is increased over its sensitive area, the sensor's resistance decreases, allowing its conductance to increase linearly with respect to applied force. When integrated with custom-built circuitry, the voltage passed through the sensor is amplified and filtered (INA114; Burr-Brown Corporation, Tucson, AZ) before being fed to an analog-to-digital converter (DAQCard-6036E; National Instruments Corporation, Austin, TX), which collects data at a 100 Hz sampling rate. Supporting software (LabVIEW 8.5 Professional; National Instruments Corporation) was developed to record, calibrate, and translate analog voltage, v(t) in V, into force detected at the sensor's location, f(t) in N.
The voltage was linearly translated to force, per Tekscan's specifications and our own sensor calibration experiments. Noise from the sensor force output was filtered using a low-pass Gaussian filter to remove frequencies >15 Hz. For more information on the onset response of the Flexiforce sensor, see section Discussion: Responsiveness of Sensor and Load Cell.
The sensor was embedded within an elastic substrate to emulate the skin environment that typically surrounds mechanoreceptors. During the pouring and curing process of a cylindrical (diameter = 30 mm; height = 10 mm) silicone-elastomer substrate (TC5005 A/B-C; BJB Enterprises Inc., Tustin, CA), the sensor was embedded at a location in the center of the x–y plane and at a depth 1.0 mm below the surface. The specific depth of 1.0 mm, roughly emulating the depth of 0.5 to 1.0 mm of the SAI afferent’s Merkel cells in human skin, ensured that the embedded sensor responded over a maximal surface area (sensor diameter = 9.53 mm; surface receptive field diameter = 20 mm), while also minimizing its depth from the surface. Furthermore, the stiffness of the silicone-elastomer, controlled by varying the percentage of cross-linker (0.98%), was manipulated to match the modulus value (Young’s modulus = 136 kPa) reported for epidermal fingertip skin of human cadavers. Although cadaver skin is stiffer than in vivo or ex vivo tissue, minimal empirical data on such stiffness exists. This modulus value is supported by finite element analysis[15, 16, 18, 19].
Spiking-sensor model: transduction sub-model
Although a sigmoid may more realistically reflect the behavior of sensory cells, the linear function minimizes the complexity of both the model and computational engine. Furthermore, the derivative represents the dynamic ramp-up phase of spike firing using only half the number of parameters required with a sigmoid.
The dynamic term, k d f’(t) as shown in Figure1c, responds to a first-order change in sensor-detected force, and therefore dominates I(t) (Figure2d) during both the ramp-up (<500 ms) and retraction portions of indentation, while the static term, k s f(t) as shown in Figure1b, responds to the magnitude of force and contributes mainly during the sustained hold. Thus, the transduction sub-model accounts for stimulus adaptation. While the full-wave model suggests the retraction of the stimulus contributes to the vigorous elicitation of action potentials in this phase, a phenomenon exhibited in neural recordings[22, 23], we did not perform an in-depth analysis here. The sub-model was implemented in C# and the values for parameters ß, k s , and k d were determined through parameter fitting as described in section Methods: Data analysis: parameter fitting.
Spiking-sensor model: neural dynamics sub-model
The third component of the spiking-sensor model was the mathematical neural dynamics sub-model. The transmembrane current was transformed into spike times (Figure2e) by abstracting leaky integrate-and-fire behavior as a resistive-capacitive (RC) circuit[24, 25]. First, current I(t) passes through a membrane with resistance R in Ohm and capacitance C in mF, to set the membrane potential u(t) in mV.
Once the membrane potential exceeds a threshold in mV, an action potential is elicited. Upon firing, the time of the spike is recorded, the membrane potential is reset to rest, and the absolute refractory period is entered, during which no spike may be elicited. This entire process iterates until stimulus offset. While the biologically equivalent values of these parameters are known (i.e., typical resting membrane potential −65 mV), this sub-model simplifies resting membrane potential to 0 mV, so that relative to this baseline, the threshold is a positive value. The absolute refractory period was set to 1 ms[27, 28].
where, g(t,u), the change in membrane potential is a function of time-dependent membrane potential, u(t), and current, I(t).
With the aid of a numerical, fourth-order Runge Kutta method for solving ordinary differential equations, the membrane potential in the next timestep u ti+1 – a linear combination of the differential equation evaluated at various points in time and membrane potential – is calculated iteratively until the future membrane potential u tN reaches or surpasses the threshold. The resolution of simulated time increments was set to a constant step size of 0.01 ms. The initial membrane potential, u t0 which represents the baseline resting membrane potential, was set to zero.
The neural dynamics sub-model implemented in C# does not account for intrinsic neuronal adaptation mechanisms. Parameters τ, C, and were determined through parameter fitting as described in section Methods: Data analysis: parameter fitting.
Experimental set-up of spiking-sensor model
Applied force was also monitored using a low-profile load cell (Sensotec 11 subminiature; maximum load = 44 N; Honeywell Inc., Columbus, OH). Although the load cell typically outputs voltage, it was calibrated by applying small indentations to a mass scale to derive normal force (resolution = 10-3 N), as shown in the middle panel of Figure4. The bottom panel of Figure4 illustrates an example of the force sensor's response to the applied indentation as a filtered continuously-detected trace with a sampling rate of 100 Hz.
The cylindrical Delrin contactor tip, threaded to the bottom of the load cell, was oriented so that its flat end provided a rigid contact against the surface of the elastic substrate (Figure2). The cylindrical tip’s dimensions were 14 mm height by 20 mm diameter. The use of this particular tip size and shape ensured that an equal state of normal stress was applied over the entire sensing area surface of the embedded sensor and minimized the out-of-plane shear stress. These conditions thereby matched those for the mouse SAI afferent.
Directly below the indenter tip, a manually-controlled and low-profile x–y stage (443 dual-mounted; Newport Corporation) supported a square aluminum tray, on which the force sensor-elastic substrate was placed and positioned in the local horizontal plane during the experiments (Figure3).
Experimental set-up of mouse SAI afferent
Single SAI afferents were isolated and identified as described in Wellnitz, Lesniak, Gerling, and Lumpkin, summarized here. Both a portion of nerve and the surrounding hairy skin of the mouse's hind paw were dissected and pinned to a 5 mm thick silicone-elastomer substrate within a custom recording and perfusion chamber. Synthetic interstitial fluid was perfused beneath the skin while individual nerve fibers were separated and draped over an electrode for differential recording. A calibrated force fiber was then used to locate the receptive fields of single afferents and estimate their mechanical thresholds. Responses were classified as SAI afferents based on conduction velocity (>9 m/s), mechanical threshold (<1 mN), punctate receptive field (diameter <0.5 mm), slow adaptation, irregular firing pattern, low spontaneous firing rate, and directional insensitivity to stretch. SAI receptor identity was confirmed by fluorescent detection of Merkel cells in the receptive field.
Over the entire receptive field of the identified afferent, a linear actuator (D-A.25AB-HT17-2-BR/4; Ultra Motion, Cutchogue, NY) commanded a flat-cylindrical MACOR contactor tip (diameter = 3.42 mm; edge round radius = 0.32 mm) to provide calibrated axial displacements. A load cell (Sensotec 31 subminiature; Honeywell Inc.) was mounted to the contactor tip to monitor reaction force in Newtons. Due to noise present in the load cell’s force readings during movements with high acceleration, a low-band pass Gaussian filter, a combination of a fast Fourier transform algorithm and its inverse, eliminated all frequencies greater than 15 Hz to reconstruct cleaned force readings.
During displacements, extracellular action potentials were recorded through a differential amplifier (1800; A-M Systems Inc., Sequim, WA) and an A/D card (DT304; Data Translation Inc., Marlboro, MA) with supporting software (SciWorks Experimenter 6.0; DataWave Technologies, Loveland, CO). Spikes were collected from single afferents isolated using Sciworks.
Common experimental procedure
Given both set-ups, a procedure with a randomized complete block design was performed on the respective experimental platforms, where treatments of ramp-and-hold indentations (coded as five unique indentation types) were blocked by replication. Each indentation type was defined by a unique combination of the rate of change of indentation during the ramp-up phase (3 levels: types I, II, & III) and the magnitude of indentation during the hold phase (3 levels: types III, IV, & V).
Indentation types IV and V shared the same ramp-up rate of change as III, but delivered increasing magnitudes of force in the sustained hold portion via 1.3 and 1.4 mm displacements – shown qualitatively in Figure6 through increases in achieved plateau height, but with similar ramp-up slope. Thirty ramp-and-hold indentations were presented to the force sensor-elastic substrate (five indentation types × six replication blocks), while 35 were presented to the mouse SAI afferent (five types × seven blocks).
Following this experimental design, the indenter tip was centered and positioned above each respective experimental platform's receptive field before conducting the replication blocks. Within each replication block, all five types were performed in random order. All indentations shared the same starting position of 0.6 mm above the experimental platform’s surface, as described in Wellnitz, et al.. Additionally, all indentations consisted of a ramp-up phase followed by a sustained ~5 s hold phase and were presented at 60 s intervals. Periods of no stimulation reduced stimulus adaptation and allowed sufficient time for the experimenter to control the indenter.
Prior to running each experiment, platforms were preconditioned. During a preconditioning session, the indenter tip – at a velocity of 6 to 8 mm/s – applied three indentations to a depth of 1.2 mm into the experimental platform from a starting gap of 0.6 mm above the platform’s surface. Preconditioning indentations were not included in the analysis. In the preconditioning process, it was evident that after three indentations, the force recorded at the probe tip was consistent between subsequent indentations of equal depth, as had been considered elsewhere. Over the first three indentations, the force recorded while the probe was held at its commanded depth decreased significantly between trials, due to initial relaxation of the elastic substrate and skin specimens.
Data analysis: independent and dependent metrics
We experimentally controlled for rate of ramp-up and magnitude of hold to elicit continuously detected force over time from the force sensor and spikes over time from the SAI afferent. The independent metrics are coded treatments of indentation types I, II, III, IV, and V and the dependent metric is average interspike interval, or, during separate phases of stimulus ramp-up () and hold (), in Eq. 6 – 8.
is calculated by dividing the sum of all ISI s from 2 to 5 seconds by the number of those ISI s, p - n. This selected window incorporates the later portion of typical double exponential adaptation of an SAI afferent. Thus, t n is the time of spike upon entering static hold phase, while t p is the spike time before ending static hold phase. Data were analyzed using ISI of the time domain, since this requires less transformation of the data through our model than a measure such as firing rate in the frequency domain. In addition, we preliminarily examined first spike latency in ms.
Figure7 illustrates the dependent metrics used as input into parameter fitting and model validation procedures. From the experiment with the spiking-sensor model, the 30 traces of sensor-detected force were randomly split amongst a training set (two-thirds of the total data set: 20 traces) and test set (remaining one-third: 10 traces). The training set was used for parameter fitting while the test set was reserved for model validation. In the experiment with the mouse SAI afferent, the 35 calculations each of and were made (five indentation types with seven replications), from which we selected five and five to represent a typical SAI afferent’s response respectively to the dynamic ramp-up and static hold, given each unique application of an indentation type (Figure7, bottom panel).
Data analysis: parameter fitting
The spiking-sensor model’s neural firing predictions for dynamic ramp-up and static hold phases were fit to those observed in the mouse SAI afferent. The parameters of the mathematical neuronal spiking model were fit using a training set, before validating via test set the fitted model's spike time predictions against those observed in the SAI afferent. Two parameter fitting sessions were performed with a different set of parameter values.
The terms and are the and pairs predicted by the spiking-sensor during replication j (amongst the 20 traces in the training set) in response to indentation type i (indentation type conditions I–V) using a specified set of model parameter values.
Data analysis: model validation
The final FSS achieved in parameter fitting was compared to an FSS calculated using the two stimulations reserved for the test set, using the fitting session’s final set of parameters. Although a high FSS was desired, the test set FSS was hypothesized to be lower than that calculated using the training set, since the alternative would indicate overfitting the model.
Next, the effect of the rate during dynamic ramp-up and the magnitude during static hold on predicted were analyzed. Since both experiments presented to the spiking-sensor model and SAI afferent were of randomized complete block design, the general treatment effect of indentation type on both and was analyzed before further pair-wise comparisons were made. An analysis of variance (ANOVA) was conducted to test the null hypothesis of no treatment effect on and. ANOVA was conducted on both the observed and predicted spike firing responses. Although the five indentation types were blocked by the order of experimental blocks, only the treatment effect was tested. The blocking effect on spiking response addresses the adaptation of the spiking-sensor model and SAI afferent to repeated stimuli. Although an intriguing topic, it is not within the scope of this work. As a follow-up post hoc analysis, a Tukey’s test was conducted for multiple pair-wise comparisons to test the individual null hypotheses of equal treatment means. All statistical analyses were conducted at a 0.05 significance level.
Parameter fitting results
Parameter fitting session
k s (mA/N)
k d (mA· s/N)
Model validation results
FSS calculated from both the training and test set were compared. Comparing the predicted spiking response to those observed upon reaching the end of parameter fitting resulted in a FSS of 0.936, which consisted of dynamic ramp-up and static hold components, 0.904 and 0.968 respectively. In comparison, the differences between predicted and observed spiking response when using the test set resulted in a lower FSS of 0.888, where the dynamic ramp-up and static hold components were 0.895 and 0.881 respectively. This decrease in performance does not suggest a poor fit since the test FSS was still greater than 0.8. Instead, it suggests the model was not overfit.
Next, the general treatment effects of indentation type on and were tested on both predicted and observed spiking response. At a significance level of 0.05, the null hypothesis, that indentation type has no effect on predicted and observed, is rejected (FNullPredicted = 16.45 > F0.05,4,20 = 2.87; FNullObserved = 3.86 > F0.05,4,24 = 2.78); thus, the means differ when both predicted by the spiking-sensor model (p-value < 0.01) and observed in the SAI afferent (p-value < 0.025). At a significance level of 0.05, the null hypothesis (HNull = indentation type has no effect on predicted and observed) is rejected (FNullPredicted = 16.52 > 2.87; FNullObserved = 11.49 > 2.78); thus, the means differ when both predicted by the spiking-sensor (p-value < 0.01) and observed in the SAI (p-value < 0.01).
Across increasing stimulus hold magnitudes, Figure9b compares the and (± SD) between predicted and observed responses. The Tukey method’s multiple pair-wise comparisons of and means across the indentation types III, IV, and V (increasing hold magnitude from 1.2, 1.3, to 1.4 mm, respectively) suggest that at a significance level of 0.05, none of the pairs ofmeans significantly differed for those predicted by the spiking-sensor (all differences < 5.01 ms, given q0.05, 5, 20 = 4.24) and those observed in the SAI afferent (all differences < 11.53 ms, given q0.05, 5, 24 = 4.17). For both spiking-sensor and SAI afferent, the pairs of means from III vs. IV and III vs. V were significantly different (differences > 27.19 and > 38.50 ms, respectively).
First spike latency results
As rate of stimulus ramp-up increased, mean first spike latencies (± SD) decreased for those predicted by the spiking-sensor (244.90 ± 8.81, 192.28 ± 6.19, 154.96 ± 3.17 ms). As magnitude of stimulus hold increased, for indentation types III, IV, and V (1.2, 1.3, 1.4 mm respectively), mean first spike latencies (± SD) were not different for those predicted by the spiking-sensor (154.96 ± 3.17, 155.525 ± 4.55, 159.055 ± 5.86 ms) although they exhibited an increasing trend. These predictions are an order of magnitude greater than first spike latencies recorded from the mouse SAI afferent, which exhibited a range from 0.81 to 96.92 ms across all indentation types, a point considered in the Discussion section.
Interspike interval irregularity results
In addition, we qualitatively compared the effect of low-pass filtering out frequencies greater than 15 Hz on the production of irregular ISI s. For instance, when spike firing was predicted from raw force sensor data, a coefficient of variation (CV) for static ISI of 0.16 was reported. When the same raw force sensor data trace was low-pass filtered (frequencies > 15 Hz), the predicted CV for static ISI decreased to 0.07. Under similar stimulation, the mouse SAI afferent reports mean CV value of 0.78 ± 0.09 for static ISI. This suggests that unfiltered data from the force sensor aids in the generation of irregular ISI s, though the resultant CV value is still not as great as that observed in the SAI afferent case.
This effort takes an integrated engineering approach to artificially replicate the neural firing behavior of an SAI afferent in its response to both magnitude and rate of indentation force by integrating a physical force sensor, embedded in an elastic substrate, with mathematical models of membrane transduction and neural dynamics. We sought to mimic particular SAI neural firing features which include increased response to changes in stimulus magnitude and rate of change, as well as the SAI afferent’s steady response to sustained stimuli and differentially greater response to moving stimuli. These responses to ramp-and-hold, sustained stimuli tie to tasks of contacting and holding an object. Importantly, the input which drives the models was the force recorded within the elastic substrate, of distinct form and delay between force or displacement at the tip of the probe and the surface of the sensor.
The effort is unique because it begins to account for skin elasticity by measuring force from within simulated skin, utilizes only six free parameters, and separates parameter fitting and model validation using response surface methodology. It also focuses upon the paradigm of the ramp-and-hold stimulus rather than vibration. By addressing these gaps, the effort extends prior works. First, prior efforts to predict the characteristic features within trains of action potentials have not taken into account the mechanics of the skin in the propagation of forces through elastic skin substrates toward the locations of the natural end organs, but rather use as input the position of the stimulus (as well as velocity, acceleration and jerk) relative to the surface contact or reaction force at the stimulus tip. However, it is clear from many past skin mechanics studies that the skin’s structure significantly influences and reshapes this propagation of forces. The forces at the tip of the indenting object are quite different from the forces at the site of end organ mechanotransduction. Moreover, particular quantities of stress and/or strain near the locations of end organs (e.g., max compressive strain and strain energy density) have been correlated to spike firing rate during the stimulus hold phase. Though the exact quantity remains an elusive question, no practical sensor can sense complex stress and strain properties, nor be robust enough to be embedded within an elastic skin-like substrate. Therefore, while work here is a start, further efforts are needed to understand and account for the many complexities of the skin (e.g., exact stress/strain quantity, rheological and tribiological characteristics, undulating ridge structure[16, 32], etc.). Second, prior efforts focused upon vibratory stimuli, whereas this work sought to mimic the elicitation of trains of action potentials, in particular, those essential features captured by the slowly adapting type I (SAI) afferent in response to everyday tasks of contacting and holding an object. Cues used here in the differentiation of various stimulus pressures, are of great importance for prosthetics users who want to hold a child’s hand or grasp a glass (i.e., “on-off” or “how much”) and match with the SAI afferent’s characteristic response to held indentation with sustained firing, and to increased hold magnitude and rate of indentation with linear increases in firing rate. The ramp-and-hold paradigm is also central to important and long-standing questions regarding SAI adaptation, first spike timing, and discrimination of displacement depth and force. Finally, this work utilizes and fits six free parameters in comparison to the 100–150 used in other works.
The model predictions fit well with observations of the characteristic features within trains of SAI action potentials elicited in the mouse. In specific, the computational model replicated basic findings regarding the mouse SAI afferent’s discrimination of levels of stimulus magnitude in the static hold phase. As the magnitude of indentation increased from 1.2 mm (indentation type III), 1.3 mm (IV), to 1.4 mm (V) the spiking-sensor model predicted a trend of increased spike firing, evident from decreasing average ISI s (98.81, 54.51, 41.11 ms respectively) within the static hold phase (Figure8b). These predictions fit favorably (mean stat_FSS = 0.925) to those observed average static ISI s (93.31, 36.28, 29.89 ms respectively) and agree with works by Mountcastle et al. and Ge and Khalsa. Mountcastle et al. demonstrated a roughly positive and linear relationship between stimulus displacement magnitude and firing rate reported by SAI afferents in the macaque monkey. An increase in magnitude from 0.2 to 0.9 mm related to an increase from an average of 0 to 70 impulses/s, or a decrease in average ISI s from 1000 to 4.29 ms. Ge and Khalsa’s work demonstrated that an increase in indentation magnitude from 0.05 to 0.2 mm related to a linear increase (R2 = 0.72) in firing rate response from 0 to 20 impulses/s, or a decrease in average ISI s from 1000 to 50 ms. While the findings from both works did not share the exact magnitude of stimuli, the trend of response to increasing magnitude is consistent. Inconsistencies in response range may be due to, but not limited to, different tactile sensitivities across animal species.
Furthermore, this work replicated the mouse SAI afferent’s discrimination of rates of stimulus ramp-up. Regarding the fit of predicted and observed ISI s in the dynamic ramp-up phase, as the rate of indentation increased during the ramp-up phase from 1.43 N/ms (indentation type I), 1.78 N/ms (II), to 2.19 N/ms (III), the spiking-sensor model predicted a trend of increased spike firing, inferred from decreasing average ISI s (21.85, 19.98, 15.42 ms respectively) in the dynamic ramp-up phase (Figure8a). These predictions fit favorably (mean dyn_FSS = 0.900) to those observed average dynamic ramp-up ISI s (24.58, 18.32, 17.86 ms respectively). The closest connection to these findings is work by Pawluk in which ramp-up stimuli were used. However, that work did not provide quantitative descriptions, with which one could directly compare. Indirectly related to the rates of stimulus ramp-up used here are vibration studies, if one considers that increased frequency corresponds indirectly. In those studies employing repeated sinusoidal vibration, increases in frequency from 10 to 60 Hz lead evoke increases the SAI afferent response from ~23 to 42 impulses/s (or a decrease ins from 43.48 to 23.81 ms) in a macaque monkey.
First spike latency
In addition to ISI, others have sought to tie the timing of the first few neural spikes to perception. We have preliminarily analyzed the time elapsed between stimulus onset and the first spike elicited, or first spike latency, as stimulus magnitude and ramp-up rate are changed. Our predicted first spike latencies exhibit decreasing trends as both ramp-up rate and hold magnitude was increased. In addition, the predictions are an order of magnitude greater than both those recorded from the mouse SAI afferent and reported by Johannson for human SAI afferents (range of ~20 to 90 ms).
We believe the discrepancy to be due to our model fitting procedure’s exclusion of this metric, filtering artifacts, responsiveness of the sensor in silicone (see section Responsiveness of sensor and load cell), and sensor location and size. Furthermore, current force sensors may not respond on the timescales of biological mechanoreceptors, which exceeds even our analog load cell. Future efforts will seek to better match first detection of stimuli to biological scales through improvements in modeling and sensor miniaturization.
Interspike interval irregularity
In addition to the presented findings, future efforts might recreate the irregularity of neuronal spike firing. Although our model predicts regular ISIs, it may be useful to match to the highly-variable ISI s that are characteristic of SAI responses[22, 37]. Noise in the nervous system may originate from innate signal transduction mechanisms, or synaptic transmission, might be a consequence of SAI neurite arbor integration, or may stem from other unknown sources. It may be simulated using a combination of or one of two methods: adding noise or manipulating sensor-output filtering. The first method could simulate the introduction of noise directly into the computation model, as demonstrated by Russell et al., where Gaussian noise is added after the filtering processes. The second method could reconsider our filtering procedure of the sensor analog voltage output and reintroduce a small band of frequencies and thereby a small level of variance in computations of both current and membrane potential.
Responsiveness of sensor and load cell
Direct nerve interfacing
Finally, toward the goal of direct nerve interfacing, we note that work has been done to deliver trains of current pulses by directly inserting electrodes into sensory fibers of the median nerve. Current pulses depolarized nerves to elicit action potentials. Some of these researchers validated action potential elicitation against expected signal trains while others recorded psychophysical graded perceptions. At present, the majority of their work has focused upon the transformation of signal from force sensor to current pulse not informed by native transduction mechanisms. While our work herein has sought to generate the characteristic features within trains of action potentials, the next steps are to link this model with current pulse waveforms that stimulate single action potentials at the peripheral afferent.
Elmer K. Kim is a research engineer in the University of Virginia's Department of Systems and Information Engineering. His research interests focus on human factors and human-computer interaction. Currently, he is developing technologies for sensory neural prosthetics. He earned his B.S and M.S. degrees in Systems Engineering from University of Virginia.
Gregory J. Gerling received his Ph.D. degree from the Department of Mechanical and Industrial Engineering, the University of Iowa, Iowa City. He is an associate professor in the Department of Systems and Information Engineering, University of Virginia, Charlottesville. His research interests include haptics, computational neuroscience, human factors/ergonomics, and human–computer interaction.
Sarah M. Bourdon is a graduate student in the Department of Neuroscience at Baylor College of Medicine. She worked as a research technician also at Baylor College of Medicine. During this time, she studied how the sense of touch is encoded. As a graduate student, she now focuses on neurodegeneration, in hopes of understanding diseases such as Alzheimer’s.
Scott A. Wellnitz and Kim, Wellnitz, Bourdon, Lumpkin, Gerling received his Ph.D. in 2010 from the Department of Neuroscience at Baylor College of Medicine. He is currently a postdoctoral associate at the J. David Gladstone Institute for Neurological Disease. His current work focuses on molecular mechanisms of Alzheimer's disease pathology.
Ellen A. Lumpkin performed her Ph.D. research in sensory neuroscience with Dr. A. James Hudspeth at UT Southwestern Medical Center and The Rockefeller University. Her research program aims to define molecular and cellular mechanisms of cutaneous touch sensation in mammals. She is an Associate Professor of Dermatology and of Physiology and Cellular Biophysics at Columbia University College of Physicians and Surgeons.
This work was supported by grants from the Defense Advanced Research Projects Agency (Grant Number HR0011-08-1-0072 and N11AP20002 to GJG), National Library of Medicine (Grant Number T15LM009462 to GJG) and National Institutes of Health (NINDS R01NS073119 to EAL and GJG, NIAMS R01 AR051219 to EAL). Portions of this study were conducted in the Depts. of Neuroscience and Molecular Physiology and Biophysics at Baylor College of Medicine. The content is solely the responsibility of the authors and does not necessarily represent the official views of DARPA, the National Library of Medicine, or the National Institutes of Health.
- Ziegler-Graham K, MacKenzie EJ, Ephraim PL, Travis TG, Brookmeyer R: Estimating the prevalence of limb loss in the united states: 2005 to 2050. Arch Phys Med Rehabil 2008, 89: 422-429. 10.1016/j.apmr.2007.11.005View ArticlePubMedGoogle Scholar
- Cipriani C, Controzzi M, Carrozza MC: The SmartHand transradial prosthesis. J Neuroeng Rehabil 2011, 8: 29. 10.1186/1743-0003-8-29PubMed CentralView ArticlePubMedGoogle Scholar
- Micera S, Rossini P, Rigosa J, Citi L, Carpaneto J, Raspopovic S, Tombini M, Cipriani C, Assenza G, Carrozza M, Hoffmann K-P, Yoshida K, Navarro X, Dario P: Decoding of grasping information from neural signals recorded using peripheral intrafascicular interfaces. J Neuroeng Rehabil 2011, 8: 53. 10.1186/1743-0003-8-53PubMed CentralView ArticlePubMedGoogle Scholar
- Raspopovic S, Carpaneto J, Udina E, Navarro X, Micera S: On the identification of sensory information from mixed nerves by using single-channel cuff electrodes. J Neuroeng Rehabil 2010, 7: 17. 10.1186/1743-0003-7-17PubMed CentralView ArticlePubMedGoogle Scholar
- Johnson KO, Popovic D, Riso RR, Koris M, Vandoren C, Kantor C: Perspectives on the role of afferent signals in control of motor neuroprostheses. Med Eng Phys 1995, 17: 481-496. 10.1016/1350-4533(95)00003-6View ArticlePubMedGoogle Scholar
- Kim SS, Sripati AP, Vogelstein RJ, Armiger RS, Russell AF, Bensmaia SJ: Conveying tactile feedback in sensorized hand neuroprostheses using a biofidelic model of mechanotransduction. IEEE Trans Biomed Circuits Syst 2009, 3: 398-404.PubMed CentralView ArticleGoogle Scholar
- Johnson KO: The roles and functions of cutaneous mechanoreceptors. Curr Opin Neurobiol 2001, 11: 455-461. 10.1016/S0959-4388(00)00234-8View ArticlePubMedGoogle Scholar
- Ounaies Z, Park C, Wise KE, Siochi EJ, Harrison JS: Electrical properties of single wall carbon nanotube reinforced polyimide composites. Comp Sci Tech 2003, 63: 1637-1646. 10.1016/S0266-3538(03)00067-8View ArticleGoogle Scholar
- Krishna GM, Rajanna K: Tactile sensor based on piezoelectric resonance. IEEE Sens J 2004, 4: 691-697. 10.1109/JSEN.2004.833505View ArticleGoogle Scholar
- Lee MH: Tactile sensing: new directions, new challenges. Int J Robot Res 2000, 19: 636-643.Google Scholar
- Sato K, Kamiyama K, Kawakami N, Tachi S: Finger-shaped gelforce: sensor for measuring surface traction fields for robotic hand. IEEE Trans Haptics 2010, 3: 37-47.View ArticleGoogle Scholar
- Phillips JR, Johnson KO: Tactile spatial resolution. III. A continuum mechanics model of skin predicting mechanoreceptor responses to bars, edges, and gratings. J Neurophysiol 1981, 46: 1204-1225.PubMedGoogle Scholar
- Dandekar K, Raju BI, Srinivasan MA: 3-D finite-element models of human and monkey fingertips to investigate the mechanics of tactile sense. J Biomech Eng Trans ASME 2003, 125: 682-691. 10.1115/1.1613673View ArticleGoogle Scholar
- Sripati AP, Bensmaia SJ, Johnson KO: A continuum mechanical model of mechanoreceptive afferent responses to indented spatial patterns. J Neurophysiol 2006, 95: 3852-3864. 10.1152/jn.01240.2005PubMed CentralView ArticlePubMedGoogle Scholar
- Maeno T, Kobayashi K, Yamazaki N: Relationship between the structure of human finger tissue and the location of tactile receptors. JSME Inter J 1998, 41: 94-100. 10.1299/jsmec.41.94View ArticleGoogle Scholar
- Gerling GJ, Thomas GW: Fingerprint lines may not directly affect SA-I mechanoreceptor response. Somatosens Mot Res 2008, 25: 61-76. 10.1080/08990220701838996View ArticlePubMedGoogle Scholar
- Quilliam TA: The structure of finger print skin. In Active Touch: The Mechanism of Recognition of Objects by Manipulation: A Multidisciplinary Approach. Pergamon, Oxford; 1978:1-18.Google Scholar
- Wu JZ, Dong RG, Rakheja S, Schopper AW, Smutz WP: A structural fingertip model for simulating of the biomechanics of tactile sensation. Med Eng Phys 2004, 26: 165-175. 10.1016/j.medengphy.2003.09.004View ArticlePubMedGoogle Scholar
- Srinivasan MA, Dandekar K: An investigation of the mechanics of tactile sense using two-dimensional models of the primate fingertip. J Biomech Eng Trans ASME 1996, 118: 48-55. 10.1115/1.2795945View ArticleGoogle Scholar
- Lesniak DR, Gerling GJ: Predicting SA-I mechanoreceptor spike times with a skin-neuron model. Math Biosci 2009, 220: 15-23. 10.1016/j.mbs.2009.03.007PubMed CentralView ArticlePubMedGoogle Scholar
- Holt JR, Corey DP: Two mechanisms for transducer adaptation in vertebrate hair cells. Proc Natl Acad Sci 2000, 97: 11730-11735. 10.1073/pnas.97.22.11730PubMed CentralView ArticlePubMedGoogle Scholar
- Wellnitz SA, Lesniak DR, Gerling GJ, Lumpkin EA: The regularity of sustained firing reveals two populations of slowly adapting touch receptors in mouse hairy skin. J Neurophysiol 2010, 103: 3378-3388. 10.1152/jn.00810.2009PubMed CentralView ArticlePubMedGoogle Scholar
- Gerstner W, Kistler WM: Spiking Neuron Models: Single Neurons, Populations, and Plasticity. Cambridge University Press, Cambridge, UK; 2002.View ArticleGoogle Scholar
- Izhikevich EM: Simple model of spiking neurons. IEEE Trans Neural Netw 2003, 14: 1569-1572. 10.1109/TNN.2003.820440View ArticlePubMedGoogle Scholar
- Nicholls JG, Martin AR, Wallace BG, Fuchs: From Neuron to Brain: A Cellular and Molecular Approach to the Function of the Nervous System. 4th edition. Sinauer Associates, Sunderland Mass; 2001.Google Scholar
- Bensmaia S: A transduction model of the meissner corpuscle. Math Biosci 2002, 176: 203-217. 10.1016/S0025-5564(02)00089-5View ArticlePubMedGoogle Scholar
- Kandel ER, Schwartz JH, Jessell TM: Principles of Neural Science. 4th edition. McGraw-Hill, New York; 2000.Google Scholar
- Ge WQ, Khalsa PS: Encoding of compressive stress during indentation by slowly adapting type i mechanoreceptors in rat hairy skin. J Neurophysiol 2002, 87: 1686-1693.PubMedGoogle Scholar
- Williams AL, Gerling GJ, Wellnitz SA, Lumpkin EA: Skin relaxation predicts neural firing rate adaptation in sai touch receptors. Proceedings of the 32nd Annual International Conference of the IEEE EMBS. Buenos Aires, Argentina 2010, 6678-6681.Google Scholar
- Iggo A, Muir AR: The structure and function of a slowly adapting touch corpuscle in hairy skin. J Physiol 1969, 1969/02/01: 763-796.View ArticleGoogle Scholar
- Scheibert J, Leurent S, Prevost A, Debregeas G: The role of fingerprints in the coding of tactile information probed with a biomimetic sensor. Science 2009, 323: 1503-1506. 10.1126/science.1166467View ArticlePubMedGoogle Scholar
- Mountcastle VB, Talbot WH, Kornhuber HH: The neural transformation of mechanical stimuli delivered to the monkey’s hand. In Touch, Heat and Pain. Edited by: Reuck AV, Knight J. Ciba Foundation, Churchill, London; 1966:325-345.Google Scholar
- Pawluk DTV: A Viscoelastic Model of the Human Fingerpad and a Holistic Model of Human Touch. 1997.Google Scholar
- Bensmaia SJ, Leung YY, Hsiao SS, Johnson KO: Vibratory adaptation of cutaneous mechanoreceptive afferents. J Neurophysiol 2005, 94: 3023-3036. 10.1152/jn.00002.2005PubMed CentralView ArticlePubMedGoogle Scholar
- Johansson RS, Birznieks I: First spikes in ensembles of human tactile afferents code complex spatial fingertip events. Nat Neurosci 2004, 7: 170-177. 10.1038/nn1177View ArticlePubMedGoogle Scholar
- Horch KW, Whitehorn D, Burgess PR: Impulse generation in type I cutaneous mechanoreceptors. J Neurophysiol 1974, 37: 267-281.PubMedGoogle Scholar
- Russell AF, Armiger RS, Vogelstein RJ, Bensmaia SJ, Etienne-Cummings R: Real-time implementation of biofidelic sa1 model for tactile feedback. Proceedings of the 31st Annual International Conference of the IEEE EMBS. Minneapolis, MN, USA 2009, 185-188.Google Scholar
- Pubols BH: Factors affecting cutaneous mechanoreceptor response. I. constant-force versus constant-displacement stimulation. J Neurophysiol 1982, 47: 530-543.PubMedGoogle Scholar
- Dhillon GS, Horch KW: Direct neural sensory feedback and control of a prosthetic arm. IEEE Trans Neural Syst Rehabil Eng 2005, 13: 468-472. 10.1109/TNSRE.2005.856072View ArticlePubMedGoogle Scholar
- Shimojo M, Suzuki T, Namiki A, Saito T, Kunimoto M, Makino R, Ogawa H, Mabuchi K: Development of a system for experiencing tactile sensation from a robot hand by electrically stimulating sensory nerve fiber. In Proceedings of the IEEE International Conference on Robotics & Automation. Volume 1; 2003:1264-1270.Google Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.