Neurostimulation success rate of repetitive‐pulse focused ultrasound in an in vivo giant axon model: An acoustic parametric study
Jérémy Vion‐Bailly, Ivan M. Suarez‐Castellanos, Jean‐Yves Chapelon, Alexandre Carpentier, W. Apoutou N'Djin
Medical Physics 2022, 49, 682-701 · 10.1002/mp.15358
Abstract
Purpose Focused ultrasound (FUS) is a promising tool to develop new modalities of therapeutic neurostimulation. The ability of FUS to stimulate the nervous system, in a noninvasive and spatiotemporally precise manner, has been demonstrated in animals and human subjects, but the underlying biomechanisms are not fully understood yet. The objective of the present study was to investigate the bioeffects involved in the generation of trains of action potentials (APs) by repetitive-pulse FUS stimuli in a simple invertebrate neural model. Methods The respective influences of different acoustic parameters on the neurostimulation success rate (NSR), defined as the rate of FUS stimuli capable of evoking at least one AP, were explored using the system of afferent nerves and giant fibers of Lumbricus terrestris as neural model. Each parameter was studied independently by administering random FUS sequences while keeping all but one FUS parameter constant. The NSR was evaluated as a function of (i) the spatial-average pulse-average intensity (I sapa ); (ii) the pulse duration (PD); (iii) the pulse repetition frequency (PRF); iv) the number of cycles per pulse (N cycles ); (v) two ultrasound frequencies, f 0 = 1.1 MHz and f 3 = 3.3 MHz, corresponding to the fundamental and third-harmonic resonant frequencies of the FUS transducer, respectively (spherical, radius of curvature: 50 mm); and (vi) levels of emerging stable cavitation and inertial cavitation. Results The NSR associated to 1.1 MHz repetitive-pulse FUS stimuli was found to increase as a function of increasing I sapa , PD, PRF, and N cycles . When evaluating each parameter at f = 1.1 MHz, it was observed that NSRs close to 100% were achieved when sufficiently elevating their respective values. When computing the NSR as a function of the spatial-average, temporal-average intensity (I sata ), defined as the product of PRF, PD, and I sapa , a significant elevation of the NSR from 0% to close to 100% was measured by increasing I sata from values approximate to 4 W/cm 2 to values higher than 12 W/cm 2 . No clear and consistent trend was observed in trials aimed at exploring the effects of different levels of stable and inertial acoustic cavitation on the NSR. Finally, the feasibility of inducing neural responses with 3.3 MHz repetitive-pulse FUS stimuli was also demonstrated with NSRs reaching up to 60%, in the range of FUS parameters studied. Conclusion The time-averaged value of the radiation force per unit volume of tissue is proportional to the acoustic intensity. As a result, the observations from this study suggest that the neural structure responding to the stimulus is sensitive to the mean radiation force carried by the FUS sequence, regardless of the combination of FUS parameters giving rise to such force. The results from this study further revealed the existence of a minimal activation threshold with regard to I sapa .
Abstract via europepmc.
Exposures
Exposure 1: Posterior end of earthworm (giant fibers), fundamental frequency 1.1 MHz, repetitive-pulse FUS with varying intensity/PD/PRF/pulse count
Target: earthworm neurons — “posterior end of ventral nerve cord / giant fibers (MGF/LGF) of Lumbricus terrestris”
Device: other named manufacturer · Ferroperm · PZ28 ✓
| Waveform | pulsed | |
|---|---|---|
| Fundamental frequency (kHz) | 1,100 | ✓✓✓ |
| Pulse duration (ms) | 0.02, 0.04, 0.08, 0.16, 0.32, 0.48, 0.64, 0.8, 0.96swept | ✓✓✓ |
| Pulse repetition frequency (Hz) | 25, 125swept | ✓✓✓ |
| Duty cycle (%) | 0.25, 12swept | ✓✓✓ |
| Sonication duration (s) | not reported | ⚑ |
| Free-field pressure (kPa) | 1,800, 5,200swept | ✓✓✓⚑ |
|---|---|---|
| Free-field Isppa (W/cm²) | 70, 260, 560, 930, 1,260, 1,520, 1,750, 1,860swept | ✓✓✓⚑ |
| Free-field Ispta (W/cm²) | 2.5, 33swept | ✓✓✓ |
| In-situ estimate | not reported | |
| In-situ pressure (kPa) | not reported | |
| In-situ Isppa (W/cm²) | not reported | |
| In-situ Ispta (W/cm²) | not reported |
A single ultrasound burst was composed of several pulses repeated at a fixed pulse repetition frequency (PRF); the burst of reference was f=f0 (1.1 MHz), PD=160 microseconds, Isapa=0.93 kW/cm2, PRF=125 Hz, Npulses=20. Pulse duration (PD) was swept over 20-960 microseconds at fixed Isapa values in n=4 animals; PRF was compared at 125 Hz vs 25 Hz (n=6 animals) with PD=160 microseconds, Isapa=0.56 kW/cm2, Npulses=20 fixed; the number of pulses per burst (Npulses) was varied among 1, 2, 3, 5, 20 and 200 in a single animal (n=1) with PD=480 microseconds, Isapa=0.26 kW/cm2, PRF=125 Hz fixed.
Exposure 2: Posterior end of earthworm (giant fibers), third-harmonic frequency 3.3 MHz, repetitive-pulse FUS with varying intensity/PD
Target: earthworm neurons — “posterior end of ventral nerve cord / giant fibers (MGF/LGF) of Lumbricus terrestris”
Device: other named manufacturer · Ferroperm · PZ28 ✓
| Waveform | pulsed | |
|---|---|---|
| Fundamental frequency (kHz) | 3,300 | ✓✓✓ |
| Pulse duration (ms) | 0.48, 0.64, 0.8, 0.96swept | ✓✓✓ |
| Pulse repetition frequency (Hz) | 125swept | ✓✓✓ |
| Duty cycle (%) | 6, 12swept | ✓✓✓⚑ |
| Sonication duration (s) | not reported | ⚑ |
| Free-field pressure (kPa) | 1,800, 5,200swept | ✓? |
|---|---|---|
| Free-field Isppa (W/cm²) | 140, 430, 770, 1,040, 1,190swept | ✓✓✓⚑ |
| Free-field Ispta (W/cm²) | 8.4, 83.2swept | ✓✓✓ |
| In-situ estimate | not reported | |
| In-situ pressure (kPa) | not reported | |
| In-situ Isppa (W/cm²) | not reported | |
| In-situ Ispta (W/cm²) | not reported |
Using the same transducer driven at its third harmonic (f3 ~ 3*f0 = 3.3 MHz), bursts were built from different combinations of Isapa and PD (PD=480-960 microseconds; Isapa=0.14-1.19 kW/cm2) with PRF=125 Hz and Npulses=20 fixed.
Flags from extraction
n_subjects— Different sub-experiments used separate cohorts of earthworms (n=4 for PD/intensity, n=6 for PRF, n=1 for pulse-count, n=2 for cavitation); the 3.3 MHz frequency study does not state an n. Recorded as a list of the stated group sizes; total unique animals unclear due to possible overlap between sub-studies.randomised— Randomisation refers to the random ordering of FUS burst types within each stimulation sequence (to avoid order/fatigue bias), not allocation to independent treatment arms; there is no sham/control group.exposures[0].free_field.isppa_w_cm2— Paper reports Isapa (spatial-average pulse-average intensity), not Isppa (spatial-peak pulse-average); recorded here as the closest available pulse-average intensity metric.exposures[1].free_field.isppa_w_cm2— Paper reports Isapa (spatial-average pulse-average intensity), not Isppa (spatial-peak pulse-average); recorded here as the closest available pulse-average intensity metric.exposures[0].free_field.pressure_kpa— The 1.8-5.2 MPa RMS pressure range is stated once for the overall FUS exposures used in the study; it is not explicitly broken down by the 1.1 MHz vs 3.3 MHz conditions, so the same range is assigned to both exposures.exposures[1].timing.duty_cycle_pct— Paper states duty cycles 'ranging from 0.6 to 0.12' for the 3.3 MHz PD sweep (480-960 microseconds at PRF=125 Hz); this ordering (max before min) and the value 0.6 (60%) appear inconsistent with PD*PRF arithmetic (which gives ~6-12%), suggesting a possible typo (0.06 instead of 0.6) in the source text. Recorded as stated (60, 12).exposures[0].timing.sonication_duration_s— The paper gives the number of pulses per burst (Npulses) and PRF but does not explicitly state a total burst/train duration in seconds; not computed here per protocol.exposures[1].timing.sonication_duration_s— As for exposure 0, total burst duration is not explicitly stated (only Npulses=20 and PRF=125 Hz); not computed here.