2.1 Cells, lipid, and lipid bubbles
In this study, we used bovine-derived carotid endothelial cells (Ito, et al., 2022; Ogawa, et al., 2024; Noguchi, et al., 2024) (referred to as “cells”, hereinafter), which were cultured at 37°C with a CO2 concentration of 5%, using Dulbecco's modified Eagle medium (DMEM, Thermo Fisher Scientific) with 10% fetal calf serum. We also prepared a phospholipid, which consisted of 1,2-distearoyl-sn-glycero-3-phosphatidyl-choline (DSPC) and 1,2-distearoyl-sn-glycero-3-phosphatidyl-ethanolamine-polyethylene glycol 2000 (DSPE-PEG2k) (Suzuki, et al., 2008; Wada, et al., 2016; Negishi, et al., 2013). To compose lipid bubbles (referred to as “LBs”, hereinafter) containing perfluoropropane (PFP, C3F8) gas, DSPC:DSPE-PEG2k was dissolved (molar ratio of 94:6) in mixed organic solvents (each containing 4 ml of chloroform). Then, phosphate-buffered saline (PBS) was added to the lipid solution and sonicated before the organic solvent evaporated. The obtained LBs had the size distribution (Ogawa, et al., 2024), which occurs naturally through the manufacturing process, with an average diameter of 1 µm and were encapsulated in a liposome using the phosphate buffer solution. Finally, we modified the LBs by conjugating a ligand of cyclic-RGD (cRGD) peptides, which covalently adhere to vascular endothelial cells via the DSPE-PEG on the LB surfaces. The conjugation of DSPE-PEG3.4k-Mal and cRGD ligand was performed during a 1-hour incubation using a rotary mixer. The molar ratio of DSPC:DSPE-PEG2k:DSPE-PEG3.4k-cRGD was 94:5:1. To prepare the LB suspension, 2 ml of liposome suspension (lipid concentration: 1 mg/ml) was diluted to yield the desired PBS concentration. The concentration [mg/ml] was derived from the weight of the lipid divided by the volume of the suspension.
In our previous studies (Chikaarashi, et al., 2022; Ito, et al., 2022; Noguchi, et al., 2024), the controllability of the cells, which was evaluated as the amount of cell retention on the wall surface of a flow path, increased in proportion to the LB concentration. However, the presence of surrounding LBs may affect not only the cells damage via cavitation but also the subsequent cell engraftment and culture process. In the following experiments, the cell concentration in the suspension was set to 1.0 × 105 /ml, and the LBs were diluted at a concentration within the range of 0–0.5 mg/ml. Here, the cell concentration was expressed with the cell number per ml, whereas the LBs concentration was with lipid weight per ml (Kajita, et al., 2023; Ogawa, et al., 2024).
2.2 Experimental setup
Figure 2 (a) shows the experimental setup for observing the behavior of the cells in BSCs passing through the flow path during ultrasound exposure (Chikaarashi, et al., 2022; Kajita, et al., 2023; Ito, et al., 2022; Ogawa, et al., 2024). This setup comprised a fluorescence microscope (Olympus, BXFM with DP74) and the 2D array transducer, which has a central frequency of 3 MHz and 128 elements (Takano, et al., 2020; Ushimizu, et al., 2018) (12 × 12 elements forming a square, and the four elements at each corner were lacked). The driving equipment (ES1144-1, Microsonic Co., Ltd.) (Chikaarashi, et al., 2022; Suzuki, et al., 2017; Hosaka, et al., 2013; Koda, et al., 2014) generated a burst wave with a minimum pulse repetition time and delay pitch of 10 ms and 5 ns, respectively, to electrically oscillate the elements individually. The maximum sound pressure of 400 kPa-pp with a duty ratio of 60% is the output limitation of this equipment.
The flow path was placed at the surface of the water. The transducer targeted the observation area in the flow path from a distance of l = 60 mm at the elevation angle of θ = 60°. Figure 2 (b) shows the schematic of the flow path made of polydimethylsiloxane (PDMS), which has the acoustic impedance of 1.03 MRayl at 3 MHz and 37 ˚C (Guillermic, et al., 2019). PDMS easily deforms in the behavior of elastic tissues. The inner cavity of the flow path has a length of 30 mm and a rectangular cross-section with a width wp = 2.1 mm and a height hp = 1.5 mm, respectively, as shown in Fig. 2 (c). From the viewpoint of fabricating a tubular structure, the cross-sectional shape of the flow path should be circular. However, in order to accurately measure the distribution of cells on the inner wall using microscopic images, a rectangular cross-section was produced so that the effect of depth of field could be ignored. The thickness of the upper and lower wall of the path was dup = 2.6 mm and dup = 1.4 mm, respectively. To coat the inner surface of the path with a basement membrane, the cavity of the flow path was filled with a mixture of Cellmatrix Type I-C (Nitta Gelatin Co., Ltd.) (Kojima, et al., 2020) and pure water, which was allowed to stand for 24 h to enable the collagen fibers to precipitate. The thickness of the membrane was estimated to be less than 10 µm.
A BSC suspension was injected during the artificial flow of the medium through the flow path maintained at a velocity of 10 mm/s. The BSCs in the suspension passing through the flow path were exposed to an acoustic field for up to 60 s. The behavior of the BSCs, which were propelled to adhere to the upper surface of the cavity, was optically observed using a fluorescence microscope. To evaluate the cell retention performance in the observation area, the cells were stained using Calcein-AM solution (Dojindo) (Ogawa, et al., 2024; Uggeri, et al., 2004; Tawakoli, et al., 2013) in advance. The cells emitted fluorescence with a wavelength of 515 nm for 490 nm excitation light to distinguish them. The fluorescent images were binarized by applying a threshold based on the average background brightness in 256-scale (8-bit) to obtain a retained area of the cells as Sret. The images were obtained with a 1.25× objective lens (NA 0.04, Olympus).
Following the retention process, the flow path was extracted from the water tank to be maintained in a CO2 incubator for culturing 24 h without turning the flow path upside down. To evaluate the cultured area of cells, the cells were not stained in the above-mentioned retention process. After the suspension in the path was drained to remove floating cells, the flow path was turned upside down to stain the living cells engrafted on the basement membrane. The living cells in the fluorescent images were distinguished from the background and dead cells to obtain a cultured area of the cells as Scul, and the images were obtained with a 5× objective lens (NA 0.15, Olympus).
2.3 Formation of acoustic fields
We designed several acoustic fields using the 2D array transducer to form an arbitrary shape of the acoustic field (Xu, et al., 2020; Riis, et al., 2024). First, a single focal (SF) acoustic field formed a single focal point in the center of the flow path. To extend the irradiation area along the path, a multifocal (MF) acoustic field was produced using tempo-spatial division emission (Takano, et al., 2020; Ushimizu, et al., 2018; Suzuki, et al., 2017), as shown in Fig. 3 (a), where the applied acoustic energy was distributed temporally and spatially. A focal point translates linearly as Fi (2 ≤ i ≤ n) along the y-direction with the spatial interval, r. Figure 3 (b) shows the time chart of the position transition of the focal point, where the repetition transition of the focal point is from F1 to Fn, each of which has emission duration, τ [s]. Number n indicates the number of preservable focal points, where one period of the tempo-spatial division emission is nτ [s]. In the following experiments, the MF acoustic field is denoted by MFn (e.g., MF3 includes three focal points).
For comparison with the MF acoustic field, a bar-shaped (BS) acoustic field was produced using the time-reversal method (Fink, 2008; Jing, et al., 2012), which refocuses an incident acoustic field back to the position of the source, irrespective of the complexity of the propagation medium, using the k-Wave mathematical platform (Treeby, et al., 2012). A BS virtual sound source was established in the space to generate the acoustic field as shown in Fig. 4. Thereafter, the ultrasound propagation was calculated to record the amplitudes and phases at the position of the 2D array. Finally, the recorded amplitudes and phases in all the elements were reproduced. In the following experiments, the BS acoustic field is denoted by BSn, where the bar length in the y-direction corresponds to (n − 1)r, to be compared with the length of MFn acoustic field.
Here, we defined the applied acoustic intensity of the aforementioned acoustic fields by measuring the distributions of sound pressure (Chikaarashi, et al., 2022; Noguchi, et al., 2024). The measured distribution (x-y plane) was divided into small grid areas with a width of w. Assuming a mean sound pressure of Pij [Pa-pp] in a small area (Sij), an applied acoustic power E [W] can be estimated as follows:

In the above equation, Z denotes the acoustic impedance of the medium. Coefficient D indicates the duty ratio in a small area. The applied acoustic intensity of spatial average temporal average (SATA) (Izadifar, et al., 2017; Acevedo, et al., 2002) was defined as ISATA = E/S [mW/cm2]. S is the area of the acoustic field defined by the area, where the relative amplitude exceeds − 20 dB of the maximum point. In the following experiments, an identical coordinate of w = 0.5 mm was used. Additionally, the duty ratio was set to D = 0.6, which was fixed due to the device limitation as mentioned above.
2.4 Measurement of acoustic fields
Figure 5 shows the normalized distributions of (a) SF, (b) MF3, (c) BS2, (d) BS3, and (e) BS4 acoustic fields. They were measured using an acoustic intensity measurement system (AIMS III, Onda Co., Ltd.) by translating the hydrophone (HNR-1000, Onda Co., Ltd.) (Seki, et al., 2019; Chikaarashi, et al., 2022; Kajita, et al., 2023) in degassed water without the flow path installed. The coordinate in Fig. 5 corresponds to that of the x-y plane (z = 0), where the flow path was installed in Fig. 2.
In the MF acoustic fields, the emission duration, τ, was set to 0.01 ms, which is faster enough than cell motion. In the MF and BS acoustic fields, the spatial interval, r, was fixed at 3 mm. The beam width, defined as the half-width of the distribution in SF and MF3, was measured at 2.5 mm. In Fig. 5 (b), the three focal points appeared to occur simultaneously since the scanning speed of the hydrophone was approximately 10 mm/s. The BS2 acoustic field has one peak in the distribution shown in Fig. 5 (c), whereas there were two peaks appeared in the BS3 and BS4 acoustic fields as shown in Figs. 5 (d) and (e).
Figure 6 shows the measured profiles of acoustic fields of SF, BS2, BS3, and BS4 along the (a) y-axis (x = z = 0) and (b) x-axis (y = z = 0). When the maximum sound pressure of SF was 400 kPa-pp, the acoustic intensity of ISATA was 143 mW/cm2. The maximum sound pressures in BS acoustic fields were adjusted to establish similar acoustic intensities, which of BS2, BS3 and BS4 were 138, 139, and 141 mW/cm2, respectively. Figure 7 shows the measured profiles of the acoustic fields of MF2, MF3, and MF4, where the SF acoustic field translated linearly in the y-direction, preserving the distribution of the SF field.