Nonspherical vibrations of microbubbles in contact with a wall: a pilot study at low mechanical index

An Acoustic Device for Ultra High-Speed Quantification of Cell Strain During Cell-Microbubble Interaction

Microbubbles utilize high-frequency oscillations under ultrasound stimulation to induce a range of therapeutic effects in cells, often through mechanical stimulation and permeabilization of cells. One of the largest challenges remaining in the field is the characterization of interactions between cells and microbubbles at therapeutically relevant frequencies. Technical limitations, such as employing sufficient frame rates and obtaining sufficient image resolution, restrict the quantification of the cell's mechanical response to oscillating microbubbles. Here, a novel methodology was developed to address many of these limitations and improve the image resolution of cell-microbubble interactions at high frame rates. A compact acoustic device was designed to house cells and microbubbles as well as a therapeutically relevant acoustic field while being compatible with a Shimadzu HPV-X camera. Cell viability tests confirmed the successful culture and proliferation of cells, and the attachment of DSPC- and cationic DSEPC-microbubbles to osteosarcoma cells was quantified. Microbubble oscillation was observed within the device at a frame rate of 5 million FPS, confirming suitable acoustic field generation and ultra high-speed image capture. High spatial resolution in these images revealed observable deformation in cells following microbubble oscillation and supported the first use of digital image correlation for strain quantification in a single cell. The novel acoustic device provided a simple, effective method for improving the spatial resolution of cell-microbubble interaction images, presenting the opportunity to develop an understanding of the mechanisms driving the therapeutic effects of oscillating microbubbles upon ultrasound exposure.

Biomechanical Sensing Using Gas Bubbles Oscillations in Liquids and Adjacent Technologies: Theory and Practical Applications

Gas bubbles present in liquids underpin many natural phenomena and human-developed technologies that improve the quality of life. Since all living organisms are predominantly made of water, they may also contain bubbles—introduced both naturally and artificially—that can serve as biomechanical sensors operating in hard-to-reach places inside a living body and emitting signals that can be detected by common equipment used in ultrasound and photoacoustic imaging procedures. This kind of biosensor is the focus of the present article, where we critically review the emergent sensing technologies based on acoustically driven oscillations of bubbles in liquids and bodily fluids. This review is intended for a broad biosensing community and transdisciplinary researchers translating novel ideas from theory to experiment and then to practice. To this end, all discussions in this review are written in a language that is accessible to non-experts in specific fields of acoustics, fluid dynamics and acousto-optics.

Internalization of targeted microbubbles by endothelial cells and drug delivery by pores and tunnels

Ultrasound insonification of microbubbles can locally enhance drug delivery by increasing the cell membrane permeability. To aid development of a safe and effective therapeutic microbubble, more insight into the microbubble-cell interaction is needed. In this in vitro study we aimed to investigate the initial 3D morphology of the endothelial cell membrane adjacent to individual microbubbles (n = 301), determine whether this morphology was affected upon binding and by the type of ligand on the microbubble, and study its influence on microbubble oscillation and the drug delivery outcome. High-resolution 3D confocal microscopy revealed that targeted microbubbles were internalized by endothelial cells, while this was not the case for non-targeted or IgG1-κ control microbubbles. The extent of internalization was ligand-dependent, since α_vβ₃-targeted microbubbles were significantly more internalized than CD31-targeted microbubbles. Ultra-high-speed imaging (~17 Mfps) in combination with high-resolution confocal microscopy (n = 246) showed that microbubble internalization resulted in a damped microbubble oscillation upon ultrasound insonification (2 MHz, 200 kPa peak negative pressure, 10 cycles). Despite damped oscillation, the cell's susceptibility to sonoporation (as indicated by PI uptake) was increased for internalized microbubbles. Monitoring cell membrane integrity (n = 230) showed the formation of either a pore, for intracellular delivery, or a tunnel (i.e. transcellular perforation), for transcellular delivery. Internalized microbubbles caused fewer transcellular perforations and smaller pore areas than non-internalized microbubbles. In conclusion, studying microbubble-mediated drug delivery using a state-of-the-art imaging system revealed receptor-mediated microbubble internalization and its effect on microbubble oscillation and resulting membrane perforation by pores and tunnels.

Effects of extracellular matrix rigidity on sonoporation facilitated by targeted microbubbles: Bubble attachment, bubble dynamics, and cell membrane permeabilization

In this study, we investigated the effects of extracellular matrix rigidity, an important physical property of microenvironments regulating cell morphology and functions, on sonoporation facilitated by targeted microbubbles, highlighting the role of microbubbles. We conducted mechanistic studies at the cellular level on physiologically relevant soft and rigid substrates. By developing a unique imaging strategy, we first resolved details of the 3D attachment configurations between targeted microbubbles and cell membrane. High-speed video microscopy then unveiled bubble dynamics driven by a single ultrasound pulse. Finally, we evaluated the cell membrane permeabilization using a small molecule model drug. Our results demonstrate that: (1) stronger targeted microbubble attachment was formed for cells cultured on the rigid substrate, while six different attachment configurations were revealed in total; (2) more violent bubble oscillation was observed for cells cultured on the rigid substrate, while one third of bubbles attached to cells on the soft substrate exhibited deformation shortly after ultrasound was turned off; (3) higher acoustic pressure was needed to permeabilize the cell membrane for cells on the soft substrate, while under the same ultrasound condition, acoustically-activated microbubbles generated larger pores as compared to cells cultured on the soft substrate. The current findings provide new insights to understand the underlying mechanisms of sonoporation in a physiologically relevant context and may be useful for the clinical translation of sonoporation.

Comparison of Bubble Size Distributions Inferred from Acoustic, Optical Visualisation, and Laser Diffraction

Bubble measurement has been widely discussed in the literature and comparison studies have been widely performed to validate the results obtained for various forms of bubble size inferences. This paper explores three methods used to obtain a bubble size distribution—optical detection, laser diffraction and acoustic inferences—for a bubble cloud. Each of these methods has advantages and disadvantages due to their intrinsic inference methodology or design flaws due to lack of specificity in measurement. It is clearly demonstrated that seeing bubbles and hearing them are substantially and quantitatively different. The main hypothesis being tested is that for a bubble cloud, acoustic methods are able to detect smaller bubbles compared to the other techniques, as acoustic measurements depend on an intrinsic bubble property, whereas photonics and optical methods are unable to “see” a smaller bubble that is behind a larger bubble. Acoustic methods provide a real-time size distribution for a bubble cloud, whereas for other techniques, appropriate adjustments or compromises must be made in order to arrive at robust data. Acoustic bubble spectrometry consistently records smaller bubbles that were not detected by the other techniques. The difference is largest for acoustic methods and optical methods, with size differences ranging from 5–79% in average bubble size. Differences in size between laser diffraction and optical methods ranged from 5–68%. The differences between laser diffraction and acoustic methods are less, and range between 0% (i.e., in agreement) up to 49%. There is a wider difference observed between the optical method, laser diffraction and acoustic methods whilst good agreement between laser diffraction and acoustic methods. The significant disagreement between laser diffraction and acoustic method (35% and 49%) demonstrates the hypothesis, as there is a higher proportion of smaller bubbles in these measurements (i.e., the smaller bubbles ‘hide’ during measurement via laser diffraction). This study, which shows that acoustic bubble spectrometry is able to detect smaller bubbles than laser diffraction and optical techniques. This is supported by heat and mass transfer studies that show enhanced performance due to increased interfacial area of microbubbles, compared to fine bubbles.

Sonoprinting liposomes on tumor spheroids by microbubbles and ultrasound

Ultrasound-triggered drug-loaded microbubbles have great potential for drug delivery due to their ability to locally release drugs and simultaneously enhance their delivery into the target tissue. We have recently shown that upon applying ultrasound, nanoparticle-loaded microbubbles can deposit nanoparticles onto cells grown in 2D monolayers, through a process that we termed "sonoprinting". However, the rigid surfaces on which cell monolayers are typically growing might be a source of acoustic reflections and aspherical microbubble oscillations, which can influence microbubble-cell interactions. In the present study, we aim to reveal whether sonoprinting can also occur in more complex and physiologically relevant tissues, by using free-floating 3D tumor spheroids as a tissue model. We show that both monospheroids (consisting of tumor cells alone) and cospheroids (consisting of tumor cells and fibroblasts, which produce an extracellular matrix) can be sonoprinted. Using doxorubicin-liposome-loaded microbubbles, we show that sonoprinting allows to deposit large amounts of doxorubicin-containing liposomes to the outer cell layers of the spheroids, followed by doxorubicin release into the deeper layers of the spheroids, resulting in a significant reduction in cell viability. Sonoprinting may become an attractive approach to deposit drug patches at the surface of tissues, thereby promoting the delivery of drugs into target tissues.

The Role of Ultrasound-Driven Microbubble Dynamics in Drug Delivery: From Microbubble Fundamentals to Clinical Translation

In the last couple of decades, ultrasound-driven microbubbles have proven excellent candidates for local drug delivery applications. Besides being useful drug carriers, microbubbles have demonstrated the ability to enhance cell and tissue permeability and, as a consequence, drug uptake herein. Notwithstanding the large amount of evidence for their therapeutic efficacy, open issues remain. Because of the vast number of ultrasound- and microbubble-related parameters that can be altered and the variability in different models, the translation from basic research to (pre)clinical studies has been hindered. This review aims at connecting the knowledge gained from fundamental microbubble studies to the therapeutic efficacy seen in in vitro and in vivo studies, with an emphasis on a better understanding of the response of a microbubble upon exposure to ultrasound and its interaction with cells and tissues. More specifically, we address the acoustic settings and microbubble-related parameters (i.e., bubble size and physicochemistry of the bubble shell) that play a key role in microbubble-cell interactions and in the associated therapeutic outcome. Additionally, new techniques that may provide additional control over the treatment, such as monodisperse microbubble formulations, tunable ultrasound scanners, and cavitation detection techniques, are discussed. An in-depth understanding of the aspects presented in this work could eventually lead the way to more efficient and tailored microbubble-assisted ultrasound therapy in the future.

SPIO labeling of endothelial cells using ultrasound and targeted microbubbles at diagnostic pressures

In vivo cell tracking of therapeutic, tumor, and endothelial cells is an emerging field and a promising technique for imaging cardiovascular disease and cancer development. Site-specific labeling of endothelial cells with the MRI contrast agent superparamagnetic iron oxide (SPIO) in the absence of toxic agents is challenging. Therefore, the aim of this in vitro study was to find optimal parameters for efficient and safe SPIO-labeling of endothelial cells using ultrasound-activated CD31-targeted microbubbles for future MRI tracking. Ultrasound at a frequency of 1 MHz (10,000 cycles, repetition rate of 20 Hz) was used for varying applied peak negative pressures (10–160 kPa, i.e. low mechanical index (MI) of 0.01–0.16), treatment durations (0–30 s), time of SPIO addition (-5 min– 15 min with respect to the start of the ultrasound), and incubation time after SPIO addition (5 min– 3 h). Iron specific Prussian Blue staining in combination with calcein-AM based cell viability assays were applied to define the most efficient and safe conditions for SPIO-labeling. Optimal SPIO labeling was observed when the ultrasound parameters were 40 kPa peak negative pressure (MI 0.04), applied for 30 s just before SPIO addition (0 min). Compared to the control, this resulted in an approximate 12 times increase of SPIO uptake in endothelial cells in vitro with 85% cell viability. Therefore, ultrasound-activated targeted ultrasound contrast agents show great potential for effective and safe labeling of endothelial cells with SPIO.

Microbubble Axial Localization Errors in Ultrasound Super-Resolution Imaging

Acoustic super-resolution imaging has allowed the visualization of microvascular structure and flow beyond the diffraction limit using standard clinical ultrasound systems through the localization of many spatially isolated microbubble signals. The determination of each microbubble position is typically performed by calculating the centroid, finding a local maximum, or finding the peak of a 2-D Gaussian function fit to the signal. However, the backscattered signal from a microbubble depends not only on diffraction characteristics of the waveform, but also on the microbubble behavior in the acoustic field. Here, we propose a new axial localization method by identifying the onset of the backscattered signal. We compare the accuracy of localization methods using in vitro experiments performed at 7-cm depth and 2.3-MHz center frequency. We corroborate these findings with simulation results based on the Marmottant model. We show experimentally and in simulations that detecting the onset of the returning signal provides considerably increased accuracy for super-resolution. Resulting experimental cross-sectional profiles in super-resolution images demonstrate at least 5.8 times improvement in contrast ratio and more than 1.8 times reduction in spatial spread (provided by 90% of the localizations) for the onset method over centroiding, peak detection, and 2-D Gaussian fitting methods. Simulations estimate that these latter methods could create errors in relative bubble positions as high as at these experimental settings, while the onset method reduced the interquartile range of these errors by a factor of over 2.2. Detecting the signal onset is, therefore, expected to considerably improve the accuracy of super-resolution.

Focal areas of increased lipid concentration on the coating of microbubbles during short tone-burst ultrasound insonification

Acoustic behavior of lipid-coated microbubbles has been widely studied, which has led to several numerical microbubble dynamics models that incorporate lipid coating behavior, such as buckling and rupture. In this study we investigated the relationship between microbubble acoustic and lipid coating behavior on a nanosecond scale by using fluorescently labeled lipids. It is hypothesized that a local increased concentration of lipids, appearing as a focal area of increased fluorescence intensity (hot spot) in the fluorescence image, is related to buckling and folding of the lipid layer thereby highly influencing the microbubble acoustic behavior. To test this hypothesis, the lipid microbubble coating was fluorescently labeled. The vibration of the microbubble (n = 177; 2.3–10.3 μm in diameter) upon insonification at an ultrasound frequency of 0.5 or 1 MHz at 25 or 50 kPa acoustic pressure was recorded with the UPMC Cam, an ultra-high-speed fluorescence camera, operated at ~4–5 million frames per second. During short tone-burst excitation, hot spots on the microbubble coating occurred at relative vibration amplitudes > 0.3 irrespective of frequency and acoustic pressure. Around resonance, the majority of the microbubbles formed hot spots. When the microbubble also deflated acoustically, hot spot formation was likely irreversible. Although compression-only behavior (defined as substantially more microbubble compression than expansion) and subharmonic responses were observed in those microbubbles that formed hot spots, both phenomena were also found in microbubbles that did not form hot spots during insonification. In conclusion, this study reveals hot spot formation of the lipid monolayer in the microbubble’s compression phase. However, our experimental results show that there is no direct relationship between hot spot formation of the lipid coating and microbubble acoustic behaviors such as compression-only and the generation of a subharmonic response. Hence, our hypothesis that hot spots are related to acoustic buckling could not be verified.

Loss of gas from echogenic liposomes exposed to pulsed ultrasound

The destruction of echogenic liposomes (ELIP) in response to pulsed ultrasound excitations has been studied acoustically previously. However, the mechanism underlying the loss of echogenicity due to cavitation nucleated by ELIP has not been fully clarified. In this study, an ultra-high speed imaging approach was employed to observe the destruction phenomena of single ELIP exposed to ultrasound bursts at a center frequency of 6 MHz. We observed a rapid size reduction during the ultrasound excitation in 139 out of 397 (35%) ultra- high-speed recordings. The shell dilation rate, which is defined as the microbubble wall velocity divided by the instantaneous radius, [Formula: see text] /R, was extracted from the radius versus time response of each ELIP, and was found to be correlated with the deflation. Fragmentation and surface mode vibrations were also observed and are shown to depend on the applied acoustic pressure and initial radius. Results from this study can be utilized to optimize the theranostic application of ELIP, e.g. by tuning the size distribution or the excitation frequency.

In Vitro Investigation of the Individual Contributions of Ultrasound-Induced Stable and Inertial Cavitation in Targeted Drug Delivery

Ultrasound-mediated targeted drug delivery is a therapeutic modality under development with the potential to treat cancer. Its ability to produce local hyperthermia and cell poration through cavitation non-invasively makes it a candidate to trigger drug delivery. Hyperthermia offers greater potential for control, particularly with magnetic resonance imaging temperature measurement. However, cavitation may offer reduced treatment times, with real-time measurement of ultrasonic spectra indicating drug dose and treatment success. Here, a clinical magnetic resonance imaging-guided focused ultrasound surgery system was used to study ultrasound-mediated targeted drug delivery in vitro. Drug uptake into breast cancer cells in the vicinity of ultrasound contrast agent was correlated with occurrence and quantity of stable and inertial cavitation, classified according to subharmonic spectra. During stable cavitation, intracellular drug uptake increased by a factor up to 3.2 compared with the control. Reported here are the value of cavitation monitoring with a clinical system and its subsequent employment for dose optimization.

Sounds of marine seeps: a study of bubble activity near a rigid boundary

A passive acoustic method for detecting environmentally dangerous gas leaks from pipelines and methane naturally leaking from the seabed has been investigated. Gas escape involves the formation and release of bubbles of different sizes. Each bubble emits a sound at a specific frequency. Determination of the bubble radius from the frequency of its signature passive acoustic emission by use of so-called Minnaert formula has a restricted area of applicability near the seabed. The point is that the inertial mass and the damping constant of the birthing bubble are markedly different from those of a free bubble. The theoretical model for the bubble volume oscillations near the seabed has been proposed and an analytical solution has been derived. It was shown that the bispherical coordinates provide separation of variables and are more suitable for analysis of the volume oscillations of these constrained bubbles. Explicit formulas have been derived, which describe the dependence of the bubble emission near a rigid wall on its size and the separation distance between the bubble and the boundary.

Lipid shedding from single oscillating microbubbles

Lipid-coated microbubbles are used clinically as contrast agents for ultrasound imaging and are being developed for a variety of therapeutic applications. The lipid encapsulation and shedding of the lipids by acoustic driving of the microbubble has a crucial role in microbubble stability and in ultrasound-triggered drug delivery; however, little is known about the dynamics of lipid shedding under ultrasound excitation. Here we describe a study that optically characterized the lipid shedding behavior of individual microbubbles on a time scale of nanoseconds to microseconds. A single ultrasound burst of 20 to 1000 cycles, with a frequency of 1 MHz and an acoustic pressure varying from 50 to 425 kPa, was applied. In the first step, high-speed fluorescence imaging was performed at 150,000 frames per second to capture the instantaneous dynamics of lipid shedding. Lipid detachment was observed within the first few cycles of ultrasound. Subsequently, the detached lipids were transported by the surrounding flow field, either parallel to the focal plane (in-plane shedding) or in a trajectory perpendicular to the focal plane (out-of-plane shedding). In the second step, the onset of lipid shedding was studied as a function of the acoustic driving parameters, for example, pressure, number of cycles, bubble size and oscillation amplitude. The latter was recorded with an ultrafast framing camera running at 10 million frames per second. A threshold for lipid shedding under ultrasound excitation was found for a relative bubble oscillation amplitude >30%. Lipid shedding was found to be reproducible, indicating that the shedding event can be controlled.

Acoustic behavior of microbubbles and implications for drug delivery

Ultrasound contrast agents are valuable in diagnostic ultrasound imaging, and they increasingly show potential for drug delivery. This review focuses on the acoustic behavior of flexible-coated microbubbles and rigid-coated microcapsules and their contribution to enhanced drug delivery. Phenomena relevant to drug delivery, such as non-spherical oscillations, shear stress, microstreaming, and jetting will be reviewed from both a theoretical and experimental perspective. Further, the two systems for drug delivery, co-administration and the microbubble as drug carrier system, are reviewed in relation to the microbubble behavior. Finally, future prospects are discussed that need to be addressed for ultrasound contrast agents to move from a pre-clinical tool into a clinical setting.

Circulatory bubble dynamics: from physical to biological aspects

Bubbles can form in the body during or after decompression from pressure exposures such as those undergone by scuba divers, astronauts, caisson and tunnel workers. Bubble growth and detachment physics then becomes significant in predicting and controlling the probability of these bubbles causing mechanical problems by blocking vessels, displacing tissues, or inducing an inflammatory cascade if they persist for too long in the body before being dissolved. By contrast to decompression induced bubbles whose site of initial formation and exact composition are debated, there are other instances of bubbles in the bloodstream which are well-defined. Gas emboli unwillingly introduced during surgical procedures and ultrasound microbubbles injected for use as contrast or drug delivery agents are therefore also discussed. After presenting the different ways that bubbles can end up in the human bloodstream, the general mathematical formalism related to the physics of bubble growth and detachment from decompression is reviewed. Bubble behavior in the bloodstream is then discussed, including bubble dissolution in blood, bubble rheology and biological interactions for the different cases of bubble and blood composition considered.

Encapsulated microbubbles and echogenic liposomes for contrast ultrasound imaging and targeted drug delivery

Micron- to nanometer-sized ultrasound agents, like encapsulated microbubbles and echogenic liposomes, are being developed for diagnostic imaging and ultrasound mediated drug/gene delivery. This review provides an overview of the current state of the art of the mathematical models of the acoustic behavior of ultrasound contrast microbubbles. We also present a review of the in vitro experimental characterization of the acoustic properties of microbubble based contrast agents undertaken in our laboratory. The hierarchical two-pronged approach of modeling contrast agents we developed is demonstrated for a lipid coated (Sonazoid™) and a polymer shelled (poly D-L-lactic acid) contrast microbubbles. The acoustic and drug release properties of the newly developed echogenic liposomes are discussed for their use as simultaneous imaging and drug/gene delivery agents. Although echogenicity is conclusively demonstrated in experiments, its physical mechanisms remain uncertain. Addressing questions raised here will accelerate further development and eventual clinical approval of these novel technologies.

Ultrafast dynamics of the acoustic vaporization of phase-change microdroplets

Acoustically sensitive emulsions are a promising tool for medical applications such as localized drug delivery. The physical mechanisms underlying the ultrasound-triggered nucleation and subsequent vaporization of these phase-change emulsions are largely unexplored. Here, the acoustic vaporization of individual micron-sized perfluoropentane (PFP) droplets is studied at a nanoseconds timescale. Highly diluted emulsions of PFP-in-water and oil-in-PFP-in-water droplets, ranging from 3.5 to 11 μm in radius, were prepared and the nucleation and growth of the vapor bubbles was imaged at frame rates of up to 20 Mfps. The droplet vaporization dynamics was observed to have three distinct regimes: (1) prior to nucleation, a regime of droplet deformation and oscillatory translations within the surrounding fluid along the propagation direction of the applied ultrasound; (2) a regime characterized by the rapid growth of a vapor bubble enhanced by ultrasound-driven rectified heat transfer; and (3) a final phase characterized by a relatively slow expansion, after ultrasound stops, that is fully dominated by heat transfer. A method to measure the moment of inception of the nucleation event with respect to the phase of the ultrasound wave is proposed. A simple physical model captures quantitatively all of the features of the subsequent vapor bubble growth.

Application of nonlinear sliding mode control to ultrasound contrast agent microbubbles

A sliding mode control system is developed and applied to a spherical model of a contrast agent microbubble that simulates its radial response to ultrasound. The model uses a compressible form of the Rayleigh-Plesset equation combined with a thin-shell model. A nonlinear control law for the second-order model is derived and used to design and simulate the controller. The effect of the controller on the contrast agent response is investigated for various control scenarios. This work demonstrates the feasibility of using a nonlinear control system to modulate the dynamic response of ultrasound contrast agents, but highlights the need for improved feedback mechanisms and control input methods. Possible applications of the nonlinear control system to contrast agents illustrated in this work include radius stabilization in the presence of an acoustic wave, radial growth and subsequent collapse, and generation of periodic radial oscillations while a contrast agent is within an acoustic forcing regime known to cause a chaotic response.

Dynamical analysis of the nonlinear response of ultrasound contrast agent microbubbles

The nonlinear response of spherical ultrasound contrast agent microbubbles is investigated to understand the effects of common shells on the dynamics. A compressible form of the Rayleigh-Plesset equation is combined with a thin-shell model developed by Lars Hoff to simulate the radial response of contrast agents subject to ultrasound. The responses of Albunex, Sonazoid, and polymer shells are analyzed through the application of techniques from dynamical systems theory such as Poincaré sections, phase portraits, and bifurcation diagrams to illustrate the qualitative dynamics and transition to chaos that occurs under certain changes in system parameters. Corresponding calculations of Lyapunov exponents provide quantitative data on the system dynamics. The results indicate that Albunex and polymer shells sufficiently stabilize the response to prevent transition to the chaotic regime throughout typical clinical ranges of ultrasound pressure and frequency. By contrast, Sonazoid shells delay the onset of chaos relative to an unshelled bubble but do not prevent it. A contour plot identifying regions of periodic and chaotic behavior over clinical ranges of ultrasound pressure and frequency is provided for Sonazoid. This work characterizes the nonlinear response of various ultrasound contrast agents, and shows that shell properties have a profound influence on the dynamics.

Single bubble acoustic characterization and stability measurement of adherent microbubbles

This article examines how the acoustic and stability characteristics of single lipid-shelled microbubbles (MBs) change as a result of adherence to a target surface. For individual adherent and non-adherent MBs, the backscattered echo from a narrowband 2-MHz, 90-kPa peak negative pressure interrogation pulse was obtained. These measurements were made in conjunction with an increasing amplitude broadband disruption pulse. It was found that, for the given driving frequency, adherence had little effect on the fundamental response of an MB. Examination of the second harmonic response indicated an increase of the resonance frequency for an adherent MB: resonance radius increasing of 0.3 ± 0.1 μm for an adherent MB. MB stability was seen to be closely related to MB resonance and gave further evidence of a change in the resonance frequency due to adherence.

20 years of ultrasound contrast agent modeling

The merits of ultrasound contrast agents (UCAs) were already known in the 1960s. It was, however, not until the 1990s that UCAs were clinically approved and marketed. In these years, it was realized that the UCAs are not just efficient ultrasound scatterers, but that their main constituent, the coated gas microbubble, acts as a nonlinear resonator and, as such, is capable of generating harmonic energy. Subharmonic, ultraharmonic, and higher harmonic frequencies of the transmitted ultrasound frequency have been reported. This opened up new prospects for their use and several detection strategies have been developed to exploit this harmonic energy to discriminate the contrast bubbles from surrounding tissue. This insight created a need for tools to study coated bubble behavior in an ultrasound field and the first models were developed. Since then, 20 years have elapsed, in which a broad range of UCAs and UCA models have been developed. Although the models have helped in understanding the responses of coated bubbles, the influence of the coating has not been fully elucidated to date and UCA models are still being improved. The aim of this review paper is to offer an overview in these developments and indicate future directions for research.

Effect of ultrasound on adherent microbubble contrast agents

An investigation into the effect of clinical ultrasound exposure on adherent microbubbles is described. A flow phantom was constructed in which targeted microbubbles were attached using biotin-streptavidin linkages. Microbubbles were insonated by broadband imaging pulses (centred at 2.25 MHz) over a range of pressures (peak negative pressure (PNP) = 60-375 kPa). Individual adherent bubbles were observed optically and classified as either being isolated or with a single neighbouring bubble. It is found that bubble detachment and deflation are two significant effects, even during low amplitude ultrasound exposure. Specifically, while at very low acoustic pressure (PNP < 75 kPa) 95% of the bubbles were not affected, at medium pressure (151 kPa < P < 225 kPa) 53% of the bubbles detached and at higher pressures (301 kPa < P < 375 kPa) 96% of the bubbles detached. In addition, more than 50% of the bubbles underwent deflation at pressures between 301 and 375 kPa. At pressures between 226 and 300 kPa, more adherent bubbles detached when there was a neighbouring bubble, suggesting the role of multiple scattering and secondary Bjerknes force on bubble detachment. The flow shear, primary and secondary Bjerknes forces exerted on each bubble were calculated and compared to the estimated forces acting on the bubble due to oscillations. The oscillation force is shown to be much higher than other forces. The mechanisms of bubble detachment are discussed.

Collective bubble dynamics near a surface in a weak acoustic standing wave field

The transport of bubbles to a neighboring surface is very important in surface chemistry, bioengineering, and ultrasonic cleaning, etc. This paper proposes a multi-bubble transport method by using an acoustic standing wave field and establishes a model that explains the multi-bubble translation by expressing the balance between Bjerknes forces and hydrodynamic forces on a bubble in a liquid medium. Results indicated that the influence of primary Bjerknes force, secondary Bjerknes force, and buoyancy force on the bubble translation depends on the position of the target bubble in the acoustic field. Moreover, it was found that increasing the size of a bubble or pressure amplitude can accelerate the bubble motion and enhance the bubble-bubble interaction. The secondary Bjerknes force between two bubbles can switch from an attractive one when they oscillate in phase to a repulsive one when the bubble oscillations are out of phase. These findings provide an insight into the multi-bubble translation near a surface and can be applied to future bubble motion control studies, especially in drug delivery, sonoporation, and ultrasonic cleaning.

Model for bubble pulsation in liquid between parallel viscoelastic layers

A model is presented for a pulsating spherical bubble positioned at a fixed location in a viscous, compressible liquid between parallel viscoelastic layers of finite thickness. The Green's function for particle displacement is found and utilized to derive an expression for the radiation load imposed on the bubble by the layers. Although the radiation load is derived for linear harmonic motion it may be incorporated into an equation for the nonlinear radial dynamics of the bubble. This expression is valid if the strain magnitudes in the viscoelastic layer remain small. Dependence of bubble pulsation on the viscoelastic and geometric parameters of the layers is demonstrated through numerical simulations.

Nonlinear dynamic behavior of microscopic bubbles near a rigid wall

The nonlinear dynamic behavior of microscopic bubbles near a rigid wall is investigated. Oscillations are driven by the ultrasonic pressure field that arises in various biomedical applications such as ultrasound imaging or targeted drug delivery. It is known that, when bubbles approach a blood-vessel wall, their linear dynamic response is modified. This modification may be very useful for real-time detection of bubbles that have found targets; in future therapeutic technologies, it may be useful for controlled release of medical agents encapsulating microbubbles. In this paper, the nonlinear response of microbubbles near a wall is studied. The Keller-Miksis-Parlitz equation is adopted, but modified to account for the presence of a rigid wall. This base model describes the time evolution of the bubble surface, which is assumed to remain spherical, and accounts for the effect of acoustic radiation losses owing to liquid compressibility in the momentum conservation. Two situations are considered: the base case of an isolated bubble in an unbounded medium, and a bubble near a rigid wall. In the latter case, the wall influence is modeled by including a symmetrically oscillating image bubble. The bubble dynamics is traced using a numerical solution of the model equation. Subsequently, Floquet theory is used to accurately detect the bifurcation point where bubble oscillations stop following the driving ultrasound frequency and undergo period-changing bifurcations. Of particular interest is the detection of the subcritical period-tripling and -quadrupling transition. The parametric bifurcation maps are obtained as functions of nondimensional parameters representing the bubble radius, the frequency and pressure amplitude of the driving ultrasound field, and the distance from the wall. It is shown that the presence of the wall generally stabilises the bubble dynamics, so that much larger values of the pressure amplitude are needed to generate nonlinear responses. Thus, a clinical protocol in which selected nonlinear harmonics are examined under varying insonation amplitudes may be useful in detecting microbubble proximity to walls.

Theoretical and experimental characterisation of magnetic microbubbles

In addition to improving image contrast, microbubbles have shown great potential in molecular imaging and drug/gene delivery. Previous work by the authors showed that considerable improvements in gene transfection efficiency were obtained using microbubbles loaded with magnetic nanoparticles under simultaneous exposure to ultrasound and magnetic fields. The aim of this study was to characterise the effect of nanoparticles on the dynamic and acoustic response of the microbubbles. High-speed video microscopy indicated that the amplitude of oscillation was very similar for magnetic and nonmagnetic microbubbles of the same size for the same ultrasound exposure (0.5 MHz, 100 kPa, 12-cycle pulse) and that this was minimally affected by an imposed magnetic field. The linear scattering to attenuation ratio (STAR) was also similar for suspensions of both bubble types although the nonlinear STAR was ~50% lower for the magnetic microbubbles. Both the video and acoustic data were supported by the results from theoretical modelling.

Combined optical and acoustical detection of single microbubble dynamics

A detailed understanding of the response of single microbubbles subjected to ultrasound is fundamental to a full understanding of the contrast-enhancing abilities of microbubbles in medical ultrasound imaging, in targeted molecular imaging with ultrasound, and in ultrasound-mediated drug delivery with microbubbles. Here, single microbubbles are isolated and their ultrasound-induced radial dynamics recorded with an ultra-high-speed camera at up to 25 million frames per second. The sound emission is recorded simultaneously with a calibrated single element transducer. It is shown that the sound emission can be predicted directly from the optically recorded radial dynamics, and vice versa, that the nanometer-scale radial dynamics can be predicted from the acoustic response recorded in the far field.

Forced vibrations of a bubble in a liquid-filled elastic vessel

There is increasing demand for accurate characterization of the in vivo behavior of microbubble agents used for ultrasound imaging and therapy. This study examines bubble-vessel interaction, in particular the propagation of disturbances along the vessel wall. Finite element simulations of a 3 μm radius microbubble suspended in a viscous liquid and enclosed in a 4 μm radius elastic vessel were performed, and the results compared with existing analytical results for wave propagation in elastic liquid-filled tubes. The vessel wall was shown to have a significant effect upon the amplitude of bubble oscillation and hence acoustic radiation from it, as well as distension of the vessel wall. It was found that the most important factor was the ratio of the excitation frequency to the natural "ring" frequency of the vessel which in turn depends upon its dimensions and mechanical properties. As this ratio increases, the motion of the vessel wall becomes increasingly localized to the site of the bubble. It was also shown that the validity of the results obtained using the applied model of vessel elasticity is limited to frequencies below the ring frequency, and this should be taken into account in the development of protocols for ultrasound safety and/or therapeutic procedures.

Nonspherical shape oscillations of coated microbubbles in contact with a wall

In this experimental study, the nonspherical and translational behavior of individual coated microbubbles of different sizes, in contact with a 20-μm thickness cellulose wall, are observed and categorized systematically. Images from two orthogonally positioned microscopes are merged and then recorded with an ultra-fast framing camera. Large nonspherical deformations were found with 2.25 MHz frequency ultrasound pulses having driving pressures from 80 to 325 kPa. A parametric model combining potential flow theory with a viscous boundary layer at the wall is developed and used to calculate stresses from the optically recorded microbubble oscillations. Peak shear stress of up to 300 kPa and normal stresses of up to 1 MPa are estimated when microbubbles are insonifed with a 2.25 MHz pulse at 325 kPa. The clinical relevance of these results is discussed.

A new transducer receive transfer function calibration method: application to microbubble backscattering cross-section measurements at high frequency

When comparing acoustic scattering experiments with theory, the relationship between the pressure generated by a scatterer at the surface of a transducer and the induced voltage must be known. Methods have been previously proposed to measure the receive transfer function that rely on several assumptions. A new, experimental method for measuring the acoustic response of a spherically-focused transducer, using a hydrophone at twice the focal distance, is proposed that requires a minimum number of assumptions and calculations. The receive transfer function of a spherically-focused, high-frequency transducer was calculated, and found to be within 10% of the receive transfer function calculated assuming reciprocity. Further, using the receive transfer function, the effective backscattering cross-section of bound microbubbles interrogated at 30 MHz was measured to be, on average, 65% of the geometric backscattering cross-section, with significant size-independent variability. These results give insight into selecting the optimal microbubble size distribution for linear microbubble imaging at high frequencies.

Sonoporation of endothelial cells by vibrating targeted microbubbles

Molecular imaging using ultrasound makes use of targeted microbubbles. In this study we investigated whether these microbubbles could also be used to induce sonoporation in endothelial cells. Lipid-coated microbubbles were targeted to CD31 and insonified at 1 MHz at low peak negative acoustic pressures at six sequences of 10 cycle sine-wave bursts. Vibration of the targeted microbubbles was recorded with the Brandaris-128 high-speed camera (~13 million frames per second). In total, 31 cells were studied that all had one microbubble (1.2-4.2 micron in diameter) attached per cell. After insonification at 80 kPa, 30% of the cells (n=6) had taken up propidium iodide, while this was 20% (n=1) at 120 kPa and 83% (n=5) at 200 kPa. Irrespective of the peak negative acoustic pressure, uptake of propidium iodide was observed when the relative vibration amplitude of targeted microbubbles was greater than 0.5. No relationship was found between the position of the microbubble on the cell and induction of sonoporation. This study shows that targeted microbubbles can also be used to induce sonoporation, thus making it possible to combine molecular imaging and drug delivery.

Surface oscillation and jetting from surface attached acoustic driven bubbles

We report on an experimental study of the onset of surface oscillation and jetting of bubbles attached to a rigid surface. The driving frequency is 16.27 kHz and the radius of the spherical capped bubble is 160 ± 5 μm. The acoustic amplitude is increased from 0 to 0.085 bar while the oscillation is recorded with a high-speed camera at 180,000 frames/s over 8100 periods of oscillations. The radial and surface modes are analyzed from a Fourier decomposition. With increasing pressure amplitude we find three regimes: pure radial oscillation, development of surface oscillations, and a chaotic surface oscillation regime. These regimes appear abrupt and are repeatable. In the chaotic regime, fast liquid jetting towards the rigid surface is observed.

Nonlinear emission from individual bound microbubbles at high frequencies

Targeted microbubbles detected with high-frequency ultrasound can establish the molecular expression of blood vessels with submillimeter resolution. To improve microbubble-specific imaging at high frequencies, the subharmonic and second harmonic signal from individual microbubbles were measured as a function of size and pressure. Single phospholipid-shell microbubbles (1.1 to 5.0 microm in diameter) bound to gelatin, co-aligned with an optical microscope and transducer, were insonated with 30 MHz Gaussian-enveloped pulses at pressures from 20 kPa to 1 MPa with -6 dB one-way bandwidths of 11%, 20% and 45%. A subharmonic signal (15 MHz) was detected above a pressure threshold of 110 kPa--independent of bandwidth. The signal peaked for microbubbles 1.60 microm in diameter subject to 20% and 11% bandwidth pulses, and 1.80 microm for 45% bandwidth pulses, for pressures up to 400 kPa, agreeing with the notion that microbubbles insonated at twice their resonant frequency preferentially emit a subharmonic component. For pressures between 400 kPa and 1 MPa, a broader range of microbubbles emitted a subharmonic signal, and microbubbles below 1.70 mum in diameter were disrupted. The second harmonic signal measured, within the limited experimental conditions, was consistent with nonlinear propagation. Further, the results shed light on the effect of the shell on the phase of the subharmonic signal with respect to the fundamental signal.

Microbubble characterization through acoustically induced deflation

Ultrasound contrast agents (UCA) populations are typically polydisperse and contain microbubbles with radii over a given range. Although the behavior of microbubbles of certain sizes might be masked by the behavior of others, the acoustic characterization of UCA is typically made on full populations. In this paper, we have combined acoustic and optical methods to investigate the response of isolated lipid-shelled microbubbles to low-pressure (49 and 62 kPa peak negative pressure) ultrasound tone bursts. These bursts induced slow deflation of the microbubbles. The experimental setup included a microscope connected to a fast camera acquiring one frame per pulse transmitted by a single-element transducer. The behavior of each bubble was measured at multiple frequencies, by cyclically changing the transmission frequency over the range of 2 to 4 MHz during subsequent pulse repetition intervals. The bubble echoes were captured by a second transducer and coherently recorded. More than 50 individual microbubbles were observed. Microbubbles with radii larger than 3 mum did not experience any size reduction. Smaller bubbles slowly deflated, generally until the radius reached a value around 1.4 microm, independent of the initial microbubble size. The detected pressure amplitude backscattered at 2.5 cm distance was very low, decreasing from about 5 Pa down to 1 Pa at 2 MHz as the bubbles deflated. The resonant radius was evaluated from the echo amplitude normalized with respect to the geometrical cross section. At 2-MHz excitation, deflating microbubbles showed highest normalized echo when the radius was 2.2 microm while at higher excitation frequencies, the resonant radius was lower. The relative phase shift of the echo during the deflation process was also measured. It was found to exceed pi/2 in all cases. A heuristic procedure based on the analysis of multiple bubbles of a same population was used to estimate the undamped natural frequency. It was found that a microbubble of 1.7 microm has an undamped natural frequency of 2 MHz. The difference between this size and the resonant radius is discussed as indicative of significant damping.

Ultrasonic characterization of ultrasound contrast agents

The main constituent of an ultrasound contrast agent (UCA) is gas-filled microbubbles. An average UCA contains billions per ml. These microbubbles are excellent ultrasound scatterers due to their high compressibility. In an ultrasound field they act as resonant systems, resulting in harmonic energy in the backscattered ultrasound signal, such as energy at the subharmonic, ultraharmonic and higher harmonic frequencies. This harmonic energy is exploited for contrast enhanced imaging to discriminate the contrast agent from surrounding tissue. The amount of harmonic energy that the contrast agent bubbles generate depends on the bubble characteristics in combination with the ultrasound field applied. This paper summarizes different strategies to characterize the UCAs. These strategies can be divided into acoustic and optical methods, which focus on the linear or nonlinear responses of the contrast agent bubbles. In addition, the characteristics of individual bubbles can be determined or the bubbles can be examined when they are part of a population. Recently, especially optical methods have proven their value to study individual bubbles. This paper concludes by showing some examples of optically observed typical behavior of contrast bubbles in ultrasound fields.

Ultrasound triggered image-guided drug delivery

The integration of therapeutic interventions with diagnostic imaging has been recognized as one of the next technological developments that will have a major impact on medical treatments. Important advances in this field are based on a combination of progress in guiding and monitoring ultrasound energy, novel drug classes becoming available, the development of smart delivery vehicles, and more in depth understanding of the mechanisms of the cellular and molecular basis of diseases. Recent research demonstrates that both pressure sensitive and temperature sensitive delivery systems hold promise for local treatment. The use of ultrasound for the delivery of drugs has been demonstrated in particular the field of cardiology and oncology for a variety of therapeutics ranging from small drug molecules to biologics and nucleic acids.

Acoustic characterization of single ultrasound contrast agent microbubbles

Individual ultrasound contrast agent microbubbles (BR14) were characterized acoustically. The bubbles were excited at a frequency of 2 MHz and at peak-negative pressure amplitudes of 60 and 100 kPa. By measuring the transmit and receive transfer functions of both the transmit and receive transducers, echoes of individual bubbles were recorded quantitatively and compared to simulated data. At 100 kPa driving pressure, a second harmonic response was observed for bubbles with a size close to their resonance size. Power spectra were derived from the echo waveforms of bubbles of different sizes. These spectra were in good agreement with those calculated from a Rayleigh-Plesset-type model, incorporating the viscoelastic properties of the phospholipid shell. Small bubbles excited below their resonance frequency have a response dominated by the characteristics of their phospholipid shell, whereas larger bubbles, excited above resonance, have a response identical to those of uncoated bubbles of similar size.

Nonspherical oscillations of ultrasound contrast agent microbubbles

The occurrence of nonspherical oscillations (or surface modes) of coated microbubbles, used as ultrasound contrast agents in medical imaging, is investigated using ultra-high-speed optical imaging. Optical tweezers designed to micromanipulate single bubbles in 3-D are used to trap the bubbles far from any boundary, enabling a controlled study of the nonspherical oscillations of free-floating bubbles. Nonspherical oscillations appear as a parametric instability and display subharmonic behavior: they oscillate at half the forcing frequency, which was fixed at 1.7 MHz in this study. Surface modes are shown to preferentially develop for a bubble radius near the resonance of radial oscillations. In the studied range of acoustic pressures, the growth of surface modes saturates at a level far below bubble breakage. With the definition of a single, dimensionless deformation parameter, the amplitude of nonspherical deformation is quantified as a function of the bubble radius (between 1.5 and 5 microm) and of the acoustic pressure (up to 200 kPa).