When two laser-generated microbubbles oscillate in anti-phase next to an adherent cell, their coupled collapse drives a directional microjet that punches localized, subcellular pores (200 nm to 2 µm) in the cell membrane, as demonstrated by Sankin, Yuan and Zhong (Phys. Rev. Lett. 105, 078101, 2010). The source paper left the physics of this poration unsettled: the relative contributions of direct jet impact (water-hammer pressure), microstreaming shear stress, and jet/vortex-induced membrane stretching and cell bending were not quantified, and no thresholds or scaling laws existed to predict pore size, uptake probability, or resealing from controllable bubble parameters. A follow-up series by the same group (PNAS 2015; AIP Conf. Proc. 2017; PNAS 2018) established that poration correlates with the membrane area strain integral produced by the shear of the jetting flow and clarified downstream calcium responses, but did not partition the three candidate mechanisms or deliver a validated general scaling law. The open problem is a mechanistically partitioned, predictive quantitative account of tandem-bubble poration, including pore-size statistics and resealing under physiological extracellular calcium.

Background

Cavitation — the formation and violent collapse of vapor bubbles in a liquid — converts concentrated mechanical energy into shock waves and high-speed liquid jets, and is deliberately harnessed in therapeutic ultrasound (drug and gene delivery), in lithotripsy (breaking kidney stones with focused shock waves), and in histotripsy (mechanically fractioning tissue with cavitation) to transiently permeabilize cell membranes, a process called sonoporation. Prentice et al. (Nature Physics 1, 107, 2005) provided direct observational evidence that a single optically triggered microbubble collapsing near a cell disrupts the membrane through the combined action of several distinct processes; the water-hammer pressure is the impulsive pressure spike produced when a moving liquid column is suddenly arrested at the membrane. Sankin, Yuan and Zhong (Phys. Rev. Lett. 105, 078101, 2010) introduced a more controllable variant: two microbubbles of about 50 µm maximum diameter, generated by nanosecond laser pulses 40 µm apart in a 25-µm-deep microfluidic chamber and driven with a controlled interpulse delay so that they oscillate in anti-phase. Their coupled oscillation alternately forms directional microjets (microstreaming velocities up to about 10 m/s) and counter-rotating vortices; the microstreaming is the bulk flow induced by the oscillating bubbles, mapped by micro-Particle Image Velocimetry (μ\muPIV, in which tracer-particle image pairs yield the velocity field). A cell placed on the bubble axis at stand-off distance SdS_d (the gap between the nearer bubble and the cell) is bent, stretched, and porated in an orientation- and proximity-dependent manner, with maximum dye-uptake probability at Sd=20 μmS_d=20\,\mu\mathrm{m} and pores of 200 nm to 2 µm; a single oscillating bubble of the same size at the same stand-off produces no poration. The authors could only estimate order-of-magnitude stresses: a wall shear stress τ≈μVj/h\tau\approx\mu V_j/h of about 1 kPa (with μ\mu the viscosity, VjV_j the jet speed, hh the half-channel height), bracketed by the roughly 0.1 kPa shear known to detach cells and the roughly 3 kPa water-hammer pressure the Letter cites for jet-impact membrane rupture. Because the treated cells sat in a solution of Trypan blue (a membrane-impermeant dye whose uptake reports poration) without extracellular Ca2+\mathrm{Ca}^{2+}, which the membrane-repair machinery requires, the observed pore sizes were acknowledged as likely upper limits, and resealing under physiological ionic conditions was unmeasured. Independent single-bubble work sharpened the threshold question: the classical quasi-static membrane rupture areal strain of 2–4% does not govern impulsive loading, under which erythrocyte membranes survive areal strains up to roughly 40% when stretched within tens of microseconds (Li, Chan and Ohl, Biophys. J., 2013), so any poration threshold must depend on loading rate and duration. In the tandem-bubble follow-up, Yuan, Yang and Zhong (PNAS 112, E7039, 2015) patterned HeLa cells on fibronectin islands (adhesive, extracellular-matrix-protein-coated patches that fix cell shape and anchorage), tracked membrane deformation through triads of attached beads, and defined the area strain εA\varepsilon_A (relative change of a small membrane surface patch) and its integral ∫εA2 dt\int\varepsilon_A^2\,dt, which scales with strain-energy-density impulse; outcomes (necrosis, repairable poration, nonporation) correlated best with this area strain integral rather than with peak strain, with a nominal poration threshold of about 55 µs for HeLa cells and 50% survival at integral values of 56–123 µs. The same group characterized the tandem-bubble jet and flow field (Yang et al., 2017) and the downstream intracellular calcium waves initiated either at poration sites or through mechanosensitive ion channels (Li et al., PNAS, 2018), yet stated that full characterization of the flow field and of the shear stress actually applied to the cell, and numerical modeling of the tandem-bubble–cell interaction, remained open problems.

Problem Statement

Establish a quantitative, mechanistically partitioned, and predictive account of single-cell membrane poration by tandem-bubble microjets, for the geometry of Sankin, Yuan and Zhong (2010): two anti-phase, laser-generated microbubbles of controllable size and interpulse phase near an adherent cell. The account must (i) quantify the relative contributions of direct jet impact (water-hammer pressure), microstreaming shear stress, and jet/vortex-induced membrane stretching and cell bending to pore nucleation and growth, and (ii) provide validated thresholds and scaling laws, expressed through controllable bubble parameters (bubble size, interpulse delay, inter-bubble distance, stand-off distance, orientation) and derived loading quantities (jet speed, shear stress, membrane area strain and its time integral), that predict the pore-size distribution, the probability of molecular uptake, and the resealing outcome — including under physiological extracellular calcium. The two parts serve one overarching objective: a causal, predictive chain from bubble dynamics through flow and membrane deformation to poration outcome, so that pore formation can be designed rather than tuned empirically. The 2015 area-strain-integral correlation for a single cell type, bubble size, and chip geometry counts as partial progress toward this objective, not as its completion, since it neither separates the three candidate mechanisms nor demonstrates predictive validity beyond its calibration configuration.

References

  1. G. N. Sankin, F. Yuan, P. Zhong, Pulsating tandem microbubble for localized and directional single-cell membrane poration, Physical Review Letters 105, 078101 (2010). DOI: 10.1103/physrevlett.105.078101, https://doi.org/10.1103/physrevlett.105.078101
  2. Fang Yuan, Chen Yang, Pei Zhong, Cell membrane deformation and bioeffects produced by tandem bubble-induced jetting flow, Proceedings of the National Academy of Sciences USA 112 (51), E7039–E7047 (2015). DOI: 10.1073/pnas.1518679112, https://doi.org/10.1073/pnas.1518679112
  3. C. Yang, A. S. Koff, F. Yuan, P. Zhong, C.-T. Hsiao, G. L. Chahine, Characterization of jet formation and flow field produced by tandem bubbles, AIP Conference Proceedings 1821, 090003 (2017). DOI: 10.1063/1.4977637, https://doi.org/10.1063/1.4977637
  4. F. Li, C. Yang, F. Yuan, D. Liao, T. Li, F. Guilak, P. Zhong, Dynamics and mechanisms of intracellular calcium waves elicited by tandem bubble-induced jetting flow, Proceedings of the National Academy of Sciences USA 115 (3), E353–E362 (2018). DOI: 10.1073/pnas.1713905115, https://doi.org/10.1073/pnas.1713905115
  5. Z. Fan, H. Liu, M. Mayer, C. X. Deng, Spatiotemporally controlled single cell sonoporation, Proceedings of the National Academy of Sciences USA 109 (41), 16486–16491 (2012). DOI: 10.1073/pnas.1208198109, https://doi.org/10.1073/pnas.1208198109
  6. F. Li, C. U. Chan, C.-D. Ohl, Yield strength of human erythrocyte membranes to impulsive stretching, Biophysical Journal 105 (4), 872–879 (2013). DOI: 10.1016/j.bpj.2013.06.045, https://doi.org/10.1016/j.bpj.2013.06.045
  7. Z. G. Li, A. Q. Liu, E. Klaseboer, J. B. Zhang, C.-D. Ohl, Single cell membrane poration by bubble-induced microjets in a microfluidic chip, Lab on a Chip 13 (6), 1144–1150 (2013). DOI: 10.1039/c3lc41252k, https://doi.org/10.1039/c3lc41252k
  8. Junjie Zhao, Jingzhu Wang, Shunxiang Cao, Dynamics of tandem bubble interaction near tissue, Physics of Fluids 36 (5), 053324 (2024). DOI: 10.1063/5.0201811, https://doi.org/10.1063/5.0201811
  9. P. Prentice, A. Cuschieri, K. Dholakia, M. Prausnitz, P. Campbell, Membrane disruption by optically controlled microbubble cavitation, Nature Physics 1 (2), 107–110 (2005). DOI: 10.1038/nphys148, https://doi.org/10.1038/nphys148
  10. M. P. Stewart, R. Langer, K. F. Jensen, Intracellular delivery by membrane disruption: mechanisms, strategies, and concepts, Chemical Reviews 118 (16), 7409–7531 (2018). DOI: 10.1021/acs.chemrev.7b00678, https://doi.org/10.1021/acs.chemrev.7b00678
  11. S. M. Chowdhury, L. Abou-Elkacem, T. Lee, J. Dahl, A. M. Lutz, Ultrasound and microbubble mediated therapeutic delivery: underlying mechanisms and future outlook, Journal of Controlled Release 326, 75–90 (2020). DOI: 10.1016/j.jconrel.2020.06.008, https://doi.org/10.1016/j.jconrel.2020.06.008
  12. Q. Zeng, H. An, C.-D. Ohl, Wall shear stress from jetting cavitation bubbles: influence of the stand-off distance and liquid viscosity, Journal of Fluid Mechanics 932, A14 (2021). DOI: 10.1017/jfm.2021.997, https://doi.org/10.1017/jfm.2021.997