Force sensor in simulated skin and neural model mimic tactile SAI afferent spiking response to ramp and hold stimuli
Journal of NeuroEngineering and Rehabilitation volume 9, Article number: 45 (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.
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
The first component of the spiking-sensor model was the physical force sensor-elastic substrate. A single piezo-resistive sensor (Flexiforce A201; Tekscan Inc., South Boston, MA) was embedded within a silicone-elastomer substrate (Figure1). This sensor is commercially available, robust to damage when embedded within a silicone substrate, and integrated relatively easily with peripheral hardware and software. It responds to normal force within a range of 0 to 4.44 N over its thin (height = 0.20 mm) and circular (diameter = 9.53 mm) pressure sensitive area.
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
The second component of the spiking-sensor model was the mathematical transduction sub-model. Force detected at the sensor in the substrate (Figure2a) was transformed into current (Figure2d), similar to how stress and/or strain applied at an SAI afferent’s end organ is transformed into current across its membrane.
Unlike previous work by Lesniak and Gerling, which transformed strain energy density into transmembrane current using a sigmoidal function, the functions (1) and (2) developed here linearly convert sensor-detected force and change in detected force, f’(t) in N/ms, into current, I(t) in mA. Its three coefficient terms are the intercept constant ß in mA, the static gain k s in mA/N, and the dynamic gain k d in mA·s/N. The ß term is intended to account for the varying baseline between sensors. The change in detected force is calculated with a step-size resolution h of 10 ms, given 100 Hz sampling rate.
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].
This leaky integrate-and-fire behavior was mathematically defined by (3) and (4), where the membrane’s time constant τ, in units of ms, is the product of the membrane’s resistance and capacitance.
These two equations were combined and rewritten as a single differential equation
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
A custom-built mechanical indenter and contactor tip were used to deliver precise ramp-and-hold stimuli in the z–direction to the force sensor-elastic substrate (Figure3).
A motion controller (ESP300; Newport Corporation, Irvine, CA) commanded a high-precision, mechanical linear z-stage (ILS100CC; maximum velocity = 100 mm/s; Newport Corporation) and reported the stage's position (resolution = 0.0001 mm), as shown in the top panel of Figure4.
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
A mechanical indenter and contactor tip was used to deliver vertical ramp-and-hold stimuli to the receptive field of an SAI afferent in an ex vivo skin-nerve preparation from a single adult female mouse (Figure5). Thirty-eight stimulations were delivered over several hours. The first three preconditioning stimuli were omitted from analysis. Details of the experimental procedure are described in the Methods: Common Experimental Procedure section. All animal use was conducted according to the National Institutes of Health Guide for the Care and Use of Laboratory Animals and was approved by the Institutional Animal Care and Use Committee of Baylor College of Medicine and the Department of Defense.
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 I, II, and III thereby shared the same commanded magnitude of force in the sustained hold portion via a constant 1.2 mm displacement, but delivered increasing levels of ramp-up rate of change – as shown qualitatively in Figure6 through increases in achieved ramp-up slope and quantitatively through increases in achieved maximum rate of force change, dForce. Although the stimuli were displacement-controlled, similar force levels were maintained between 1000 and 5000 ms of indentation (Figure7).
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.
or simply the time difference between two consecutive spikes, is calculated by dividing the sum of all ISI s within a given time window by the total number of ISI s contained within the same window. Thus, the dynamic (7) or
is calculated by dividing the sum of all ISI s from the time of stimulus onset to the time of peak indentation by the total number of ISI s, m, within the same window. Stimulus onset, for the spiking-sensor model, was defined as when the magnitude of applied force from the load cell exceeded 1 N. Stimulus onset, for the mouse SAI, was defined as when dForce from the load cell exceeded 0.5 mN/ms. Time of peak indentation for both set-ups was defined as when magnitude of applied force reached maximum. The static (8) or
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 parameters were fit using response surface methodology (RSM). The general RSM process is typically driven by the goal of optimizing a response, which is influenced by one or many factors. This work strives to maximize the fractional sum of squares (FSS), previously used by Phillips and Johnson (1981), between predicted and observed spike times, which is influenced by the six model parameters: ß, k s , k d , τ, C, and. Specifically, FSS is a complement of normalized sum of squared deviations and is contributed equally by dynamic (ω dyn = 0.5) and static (ω stat = 0.5) comparisons of spiking response (9).
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
Table1 focuses on the second session and provides the values of the six model parameters (ß, k s , k d , τ, C, and) and their respective calculated FSS for each of the iteration steps. We demonstrated that RSM is sensitive to starting conditions by the different paths taken by the two fitting sessions. The first session ended with a parameter set that produced a worse fit (FSS = 0.829 with parameter values 4.33E-08 mA, 5.74E-07 mA/N, 1.01E-03 mA· s/N, 71.592 ms, 1.01E-06 mF, 50.723 mV) than that produced by the second session’s set of parameters (FSS = 0.936 with parameter values 2.72E-08 mA, 6.20E-07 mA/N, 2.71E-04 mA· s/N, 71.409 ms, 9.70E-07 mF, 47.300 mV).
All following analyses focus on using the parameter values obtained during the second parameter fitting session. Using the final model parameter values obtained in the last iteration step of the second parameter fitting session, Figure8 portrays single translations of sensor-detected force into trains of spike times across indentation types. These qualitative results demonstrate increased spike density in the static phase as the indentation’s magnitude increases. In addition, increased spike density and earlier first spike latency can be seen in the dynamic phase as the indentation’s ramp-up rate of change increases.
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 ramp-up rates, Figure9a compares the mean and (±SD) between predicted and observed responses. The Tukey method’s multiple pair-wise comparisons of and means across the indentation types I, II, and III (increasing ramp-up rate from low, medium, to high, respectively) suggest that at a significance level of 0.05, only the I vs. III type pairs of means significantly differed (difference > 5.01 ms, given q0.05, 5, 20 = 4.24) for the spiking-sensor. None of the pairs of means observed in the SAI afferent were significantly different (all differences < 11.53 ms, given q0.05, 5, 24 = 4.17). For both the spiking-sensor and SAI afferent, all type pairs I vs. II, I vs. III, and II vs. III of means were not significantly different (all differences < 27.19 and < 38.50 ms, respectively).
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
The onset response of the Flexiforce sensor, apart from the silicone-elastomer substrate was also analyzed. As is observable in Figure10, the response of the load cell’s trace lags the displacement trace slightly, as does that of the Flexiforce sensor. However, the onset timing of the responses of Flexiforce sensor and load cell are identical. This indicates that the temporal lag in sensor response (Figure4, lower trace) is predominantly due to its encasement in an elastic substrate, analogous to the skin’s effects on mechanoreceptor response in the biology.
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.
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.005
Cipriani C, Controzzi M, Carrozza MC: The SmartHand transradial prosthesis. J Neuroeng Rehabil 2011, 8: 29. 10.1186/1743-0003-8-29
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-53
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-17
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-6
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.
Johnson KO: The roles and functions of cutaneous mechanoreceptors. Curr Opin Neurobiol 2001, 11: 455-461. 10.1016/S0959-4388(00)00234-8
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-8
Krishna GM, Rajanna K: Tactile sensor based on piezoelectric resonance. IEEE Sens J 2004, 4: 691-697. 10.1109/JSEN.2004.833505
Lee MH: Tactile sensing: new directions, new challenges. Int J Robot Res 2000, 19: 636-643.
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.
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.
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.1613673
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.2005
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.94
Gerling GJ, Thomas GW: Fingerprint lines may not directly affect SA-I mechanoreceptor response. Somatosens Mot Res 2008, 25: 61-76. 10.1080/08990220701838996
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.
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.004
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.2795945
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.007
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.11730
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.2009
Gerstner W, Kistler WM: Spiking Neuron Models: Single Neurons, Populations, and Plasticity. Cambridge University Press, Cambridge, UK; 2002.
Izhikevich EM: Simple model of spiking neurons. IEEE Trans Neural Netw 2003, 14: 1569-1572. 10.1109/TNN.2003.820440
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.
Bensmaia S: A transduction model of the meissner corpuscle. Math Biosci 2002, 176: 203-217. 10.1016/S0025-5564(02)00089-5
Kandel ER, Schwartz JH, Jessell TM: Principles of Neural Science. 4th edition. McGraw-Hill, New York; 2000.
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.
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.
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.
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.1166467
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.
Pawluk DTV: A Viscoelastic Model of the Human Fingerpad and a Holistic Model of Human Touch. 1997.
Bensmaia SJ, Leung YY, Hsiao SS, Johnson KO: Vibratory adaptation of cutaneous mechanoreceptive afferents. J Neurophysiol 2005, 94: 3023-3036. 10.1152/jn.00002.2005
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/nn1177
Horch KW, Whitehorn D, Burgess PR: Impulse generation in type I cutaneous mechanoreceptors. J Neurophysiol 1974, 37: 267-281.
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.
Pubols BH: Factors affecting cutaneous mechanoreceptor response. I. constant-force versus constant-displacement stimulation. J Neurophysiol 1982, 47: 530-543.
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.856072
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.
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.
All authors declare that they have no competing interests.
All authors have made substantial contributions in conception, design, data acquisition, analysis, and interpretation of data. All authors have been involved in drafting the manuscript and revising it critically. All authors have read and approved of the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Kim, E.K., Wellnitz, S.A., Bourdon, S.M. et al. Force sensor in simulated skin and neural model mimic tactile SAI afferent spiking response to ramp and hold stimuli. J NeuroEngineering Rehabil 9, 45 (2012). https://doi.org/10.1186/1743-0003-9-45