Accounting for the stability of microbubbles to multi-pulse excitation using a lipid-shedding model

Dependence of sonoporation efficiency on microbubble size: An in vitro monodisperse microbubble study

Sonoporation is the process where intracellular drug delivery is facilitated by ultrasound-driven microbubble oscillations. Several mechanisms have been proposed to relate microbubble dynamics to sonoporation including shear and normal stress. The present work aims to gain insight into the role of microbubble size on sonoporation and thereby into the relevant mechanism(s) of sonoporation. To this end, we measured the sonoporation efficiency while varying microbubble size using monodisperse microbubble suspensions. Sonoporation experiments were performed in vitro on cell monolayers using a single ultrasound pulse with a fixed frequency of 1 MHz while the acoustic pressure amplitude and pulse length were varied at 250, 500, and 750 kPa, and 10, 100, and 1000 cycles, respectively. Sonoporation efficiency was quantified using flow cytometry by measuring the FITC-dextran (4 kDa and 2 MDa) fluorescence intensity in 10,000 cells per experiment to average out inherent variations in the bioresponse. Using ultra-high-speed imaging at 10 million frames per second, we demonstrate that the bubble oscillation amplitude is nearly independent of the equilibrium bubble radius at acoustic pressure amplitudes that induce sonoporation (≥ 500 kPa). However, we show that sonoporation efficiency is strongly dependent on the equilibrium bubble size and that under all explored driving conditions most efficiently induced by bubbles with a radius of 4.7 μm. Polydisperse microbubbles with a typical ultrasound contrast agent size distribution perform almost an order of magnitude lower in terms of sonoporation efficiency than the 4.7-μm bubbles. We elucidate that for our system shear stress is highly unlikely the mechanism of action. By contrast, we show that sonoporation efficiency correlates well with an estimate of the bubble-induced normal stress.

Photoacoustic Impulse Response of Lipid-Coated Ultrasound Contrast Agents

The utility of ultrasound imaging and therapy with microbubbles may be greatly enhanced by determining their impulse-response dynamics as a function of size and composition. Prior methods for microbubble characterization utilizing high-speed cameras, acoustic transducers and laser-based techniques typically scan a limited frequency range. Here, we report on the use of a novel photoacoustic technique to measure the impulse response of single microbubbles. Individual microbubbles are driven with a broadband photoacoustic wave generated by a nanosecond-pulse laser illuminating an optical absorber. The resulting microbubble oscillations were detected by following transmission of a second laser as it passes twice through the microbubble. The system could even resolve oscillations resulting from a single-shot. As a proof-of-concept study, the size-dependent, linear impulse response of lipid-coated microbubbles was characterized using this technique. This unique method of microbubble characterization with exceptional spatiotemporal resolution opens new avenues for capturing and analyzing microbubble system dynamics.

Nonlinear dynamics and bifurcation structure of ultrasonically excited lipid coated microbubbles

In many applications, microbubbles (MBs) are encapsulated by a lipid coating to increase their stability. However, the complex behavior of the lipid coating including buckling and rupture sophisticates the dynamics of the MBs and as a result the dynamics of the lipid coated MBs (LCMBs) are not well understood. Here, we investigate the nonlinear behavior of the LCMBs by analyzing their bifurcation structure as a function of acoustic pressure. We show that, the LC can enhance the generation of period 2 (P2), P3, higher order subharmonics (SH), superharmonics and chaos at very low excitation pressures (e.g. 1 kPa). For LCMBs sonicated by their SH resonance frequency and in line with experimental observations with increasing pressure, P2 oscillations exhibit three stages: generation at low acoustic pressures, disappearance and re-generation. Within non-destructive oscillation regimes and by pressure amplitude increase, LCMBs can also exhibit two saddle node (SN) bifurcations resulting in possible abrupt enhancement of the scattered pressure. The first SN resembles the pressure dependent resonance phenomenon in uncoated MBs and the second SN resembles the pressure dependent SH resonance. Depending on the initial surface tension of the LCMBs, the nonlinear behavior may also be suppressed for a wide range of excitation pressures.

Ultrasound Contrast Agent Modeling: A Review

Ultrasound is extensively used in medical imaging, being safe and inexpensive and operating in real time. Its scope of applications has been widely broadened by the use of ultrasound contrast agents (UCAs) in the form of microscopic bubbles coated by a biocompatible shell. Their increased use has motivated a large amount of research to understand and characterize their physical properties as well as their interaction with the ultrasound field and their surrounding environment. Here we review the theoretical models that have been proposed to study and predict the behavior of UCAs. We begin with a brief introduction on the development of UCAs. We then present the basics of free-gas-bubble dynamics upon which UCA modeling is based. We review extensively the linear and non-linear models for shell elasticity and viscosity and present models for non-spherical and asymmetric bubble oscillations, especially in the presence of surrounding walls or tissue. Then, higher-order effects such as microstreaming, shedding and acoustic radiation forces are considered. We conclude this review with promising directions for the modeling and development of novel agents.

Lipid Shell Retention and Selective Binding Capability Following Repeated Transient Acoustic Microdroplet Vaporization

Targeted therapy and molecular imaging using ultrasound have been widely explored using microbubble contrast agents, and more recently, activatable droplet contrast agents that vaporize when exposed to focused ultrasound have been explored. These droplets are coated with a stabilizing, functionalizable shell, typically comprised of fully saturated phospholipids. While the shedding of the lipid shell under ultrasound exposure is a well-studied phenomenon in microbubbles, it has not been fully explored in droplet-based contrast agents, particularly in those that undergo a reversible phase change and recondense following vaporization. Here, we investigate the retention of the lipid shell following repeated transient vaporization events. Two separate fluorescent markers were used to track individual lipid subpopulations: PEGylated lipids, to which targeting ligands are typically bound, and non-PEGylated lipids, which primarily contribute to droplet stability. Following confirmation of the homogeneous surface distribution of each subpopulation of shell lipids using confocal microscopy, high-speed optical imaging provided visual evidence of the ability to repeatedly induce vaporization and recondensation in micron-scale droplets using 5.208 MHz, 3.17 MPa focused ultrasound pulses transmitted from an imaging transducer. Flow cytometry analysis indicated that while PEGylated lipids were fully retained following repeated transient phase change events, 20% of the bulk lipids were shed. While this likely contributed to an observed significant reduction in the average droplet diameter, the selective binding capabilities of droplets functionalized with an RGD peptide, targeted to the integrin αvβ3, were not affected. These results indicate that repeated droplet activation may promote shifts in the droplet size distribution but will not influence the accuracy of targeting for therapy or molecular imaging.

Microbubble Agents: New Directions

Microbubble ultrasound contrast agents have now been in use for several decades and their safety and efficacy in a wide range of diagnostic applications have been well established. Recent progress in imaging technology is facilitating exciting developments in techniques such as molecular, 3-D and super resolution imaging and new agents are now being developed to meet their specific requirements. In parallel, there have been significant advances in the therapeutic applications of microbubbles, with recent clinical trials demonstrating drug delivery across the blood-brain barrier and into solid tumours. New agents are similarly being tailored toward these applications, including nanoscale microbubble precursors offering superior circulation times and tissue penetration. The development of novel agents does, however, present several challenges, particularly regarding the regulatory framework. This article reviews the developments in agents for diagnostic, therapeutic and "theranostic" applications; novel manufacturing techniques; and the opportunities and challenges for their commercial and clinical translation.

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.

High-precision acoustic measurements of the nonlinear dilatational elasticity of phospholipid coated monodisperse microbubbles

The acoustic response of phospholipid-coated microbubble ultrasound contrast agents (UCA) is dramatically affected by their stabilizing shell. The interfacial shell elasticity increases the resonance frequency, the shell viscosity increases damping, and the nonlinear shell viscoelasticity increases the generation of harmonic echoes that are routinely used in contrast-enhanced ultrasound imaging. To date, the surface area-dependent interfacial properties of the phospholipid coating have never been measured due to the extremely short time scales of the MHz frequencies at which the microscopic bubbles are driven. Here, we present high-precision acoustic measurements of the dilatational nonlinear viscoelastic shell properties of phospholipid-coated microbubbles. These highly accurate measurements are now accessible for the first time by tuning the surface dilatation, that is, the lipid packing density, of well-controlled monodisperse bubble suspensions through the ambient pressure. Upon compression, the shell elasticity of bubbles coated with DPPC and DPPE-PEG5000 was found to increase up to an elasticity of 0.6 N m-1 after which the monolayer collapses and the elasticity vanishes. During bubble expansion, the elasticity drops monotonically in two stages, first to an elasticity of 0.35 N m-1, and then more rapidly to zero. Integration of the elasticity vs. surface area curves showed that, indeed, a phospholipid-coated microbubble is in a tensionless state upon compression, and that it reaches the interfacial tension of the surrounding medium upon expansion. The measurements presented in this work reveal the detailed features of the nonlinear dilatational shell behavior of micron-sized lipid-coated bubbles.

Universal Equations for the Coalescence Probability and Long-Term Size Stability of Phospholipid-Coated Monodisperse Microbubbles Formed by Flow Focusing

Resonantly driven monodisperse phospholipid-coated microbubbles are expected to substantially increase the sensitivity and efficiency in contrast-enhanced ultrasound imaging and therapy. They can be produced in a microfluidic flow-focusing device, but questions remain as to the role of the device geometry, the liquid and gas flow, and the phospholipid formulation on bubble stability. Here, we develop a model based on simple continuum mechanics equations that reveals the scaling of the coalescence probability with the key physical parameters. It is used to characterize short-term coalescence behavior and long-term size stability as a function of flow-focusing geometry, bulk viscosity, lipid cosolvent mass fraction, lipid concentration, lipopolymer molecular weight, and lipopolymer molar fraction. All collected data collapse on two master curves given by universal equations for the coalescence probability and the long-term size stability. This work is therefore a route to a more fundamental understanding of the physicochemical monolayer properties of microfluidically formed bubbles and their coalescence behavior in a flow-focusing device.

Ultrasound propagation through dilute polydisperse microbubble suspensions

In a fully nonlinear model of wave propagation through bubbly media, computational complexity arises when the medium contains a polydisperse bubble population. This is because a nonlinear ordinary differential equation governing the bubble response must be solved for the current radius of each bubble size present at every spatial location and at every time step. In biomedical ultrasound imaging, commercial contrast agents typically possess a wide range of bubble sizes that exhibit a variety of differing behaviours at ultrasound frequencies of clinical interest. Despite the advent of supercomputing resources, the simulation of ultrasound propagation through microbubble populations still represents a formidable numerical task. Consequently, efficient computational algorithms that have the potential to be implemented in real time on clinical scanners remain highly desirable. In this work, a numerical approach is investigated that computes only a single ordinary differential equation at each spatial location which can potentially reduce significantly the computational effort. It is demonstrated that, under certain parameter regimes, the approach replicates the fully nonlinear model of an incident ultrasound pulse propagating through a polydisperse population of bubbles with a high degree of accuracy.

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.

Characterization of Contrast Agent Microbubbles for Ultrasound Imaging and Therapy Research

The high efficiency with which gas microbubbles can scatter ultrasound compared with the surrounding blood pool or tissues has led to their widespread employment as contrast agents in ultrasound imaging. In recent years, their applications have been extended to include super-resolution imaging and the stimulation of localized bio-effects for therapy. The growing exploitation of contrast agents in ultrasound and in particular these recent developments have amplified the need to characterize and fully understand microbubble behavior. The aim in doing so is to more fully exploit their utility for both diagnostic imaging and potential future therapeutic applications. This paper presents the key characteristics of microbubbles that determine their efficacy in diagnostic and therapeutic applications and the corresponding techniques for their measurement. In each case, we have presented information regarding the methods available and their respective strengths and limitations, with the aim of presenting information relevant to the selection of appropriate characterization methods. First, we examine methods for determining the physical properties of microbubble suspensions and then techniques for acoustic characterization of both suspensions and single microbubbles. The next section covers characterization of microbubbles as therapeutic agents, including as drug carriers for which detailed understanding of their surface characteristics and drug loading capacity is required. Finally, we discuss the attempts that have been made to allow comparison across the methods employed by various groups to characterize and describe their microbubble suspensions and promote wider discussion and comparison of microbubble behavior.

Analysis of the Uncertainty in Microbubble Characterization

There is increasing interest in the use of microbubble contrast agents for quantitative imaging applications such as perfusion and blood pressure measurement. The response of a microbubble to ultrasound excitation is, however, extremely sensitive to its size, the properties of its coating and the characteristics of the sound field and surrounding environment. Hence the results of microbubble characterization experiments can be significantly affected by experimental uncertainties, and this can limit their utility in predictive modelling. The aim of this study was to attempt to quantify these uncertainties and their influence upon measured microbubble characteristics. Estimates for the parameters characterizing the microbubble coating were obtained by fitting model data to numerical simulations of microbubble dynamics. The effect of uncertainty in different experimental parameters was gauged by modifying the relevant input values to the fitting process. The results indicate that even the minimum expected uncertainty in, for example, measurements of microbubble radius using conventional optical microscopy, leads to variations in the estimated coating parameters of ∼20%. This should be taken into account in designing microbubble characterization experiments and in the use of data obtained from them.

Stability of Monodisperse Phospholipid-Coated Microbubbles Formed by Flow-Focusing at High Production Rates

Monodisperse microbubble ultrasound contrast agents may dramatically increase the sensitivity and efficiency in ultrasound imaging and therapy. They can be produced directly in a microfluidic flow-focusing device, but questions remain as to the interfacial chemistry, such as the formation and development of the phospholipid monolayer coating over time. Here, we demonstrate the synthesis of monodisperse bubbles with radii of 2-10 μm at production rates ranging from 10(4) to 10(6) bubbles/s. All bubbles were found to dissolve to a stable final radius 2.55 times smaller than their initial radius, independent of the nozzle size and shear rate, indicating that the monolayer self-assembles prior to leaving the nozzle. The corresponding decrease in surface area by a factor 6.6 reveals that lipid molecules are adsorbed to the gas-liquid interface in the disordered expanded state, and they become mechanically compressed by Laplace pressure-driven bubble dissolution to a more ordered condensed state with near zero surface tension. Acoustic characterization of the stabilized microbubbles revealed that their shell stiffness gradually increased from 0.8 to 2.5 N/m with increasing number of insonations through the selective loss of the more soluble lipopolymer molecules. This work therefore demonstrates high-throughput production of clinically relevant monodisperse contrast microbubbles with excellent control over phospholipid monolayer elasticity and microbubble resonance.

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.

Measuring absolute blood pressure using microbubbles

Gas microbubbles are highly compressible, which makes them very efficient sound scatterers. As another consequence of their high compressibility, the radii of the microbubbles are affected by the pressure of the fluid around them, which changes their resonance frequency. Although the pressures present within the human body cause only minor variations in the radii of uncoated microbubbles (∼0.2% per 10 mmHg) and, therefore, very small variations in the resonance frequency (∼1 kHz per 10 mmHg), it was found in the work described here, through both simulations and in vitro measurements, that large changes in resonance frequency can occur in phospholipid-coated microbubbles for small blood pressure variations because of the exotic buckling dynamics of phospholipid monolayers (up to 240 kHz per 10 mmHg). This method should allow non-invasive measurement of the gauge blood pressure in deep blood vessels as long as the microbubble physical properties are well controlled.

Surfactant shedding and gas diffusion during pulsed ultrasound through a microbubble contrast agent suspension

Interest in coated microbubbles as agents for therapeutic and quantitative imaging applications in biomedical ultrasound has increased the need for their accurate theoretical characterization. Effects such as gas diffusion, variation in the properties of the coating and the resulting changes in bubble behavior under repeated exposure to ultrasound pulses are, however, still not well understood. In this study, a revised equation for microbubble motion is proposed that includes the effects of gas diffusion, as well as adsorption, desorption and shedding of a surfactant from the bubble surface. This is incorporated into a nonlinear wave propagation model to account for these additional time dependent effects in the response of microbubble populations. The results from the model indicate there can be significant changes in both bubble behavior and the propagated pulse over time. This is in agreement with existing experimental data but is not predicted by existing propagation models. The analysis indicates that changes in bubble dynamics are dominated by surfactant shedding on the timescale of a diagnostic ultrasound pulse and gas diffusion over the timescale of the pulse repetition frequency. The implications of these results for the development of more accurate algorithms for quantitative imaging and for therapeutic applications are discussed.

Temperature-dependent biphasic shrinkage of lipid-coated bubbles in ultrasound

Lipid-coated microbubbles and emulsions are of interest as possible ultrasound-mediated drug delivery vehicles and for their interesting behaviors and fundamental properties. We and others have noted that bubbles coated with the long chain saturated phospholipid distearoylphosphatidylcholine (DSPC) rapidly shrink to a quasistable size when repeatedly insonated with short ultrasound pulses; such stability may adversely affect the bubble's subsequent ability to deliver its pharmacological cargo. Bubbles coated with the unsaturated lipid dioleoylphosphatidylcholine (DOPC) did not show stability but did undergo an abrupt change from rapid initial shrinkage to a slow persistent shrinkage, leading ultimately to dissolution or dispersion. As DOPC and DSPC differ not only in chain saturation but also phase behavior, we performed additional studies using dimyristoyl PC (DMPC) as a coat lipid and controlled the solution temperature to study bubble behavior on exposure to repeated ultrasound pulses for the same coat, in both fluid and gel phases. We find, first, that essentially all bubbles show an initially rapid shrinkage, in which gas loss exceeds the limit imposed by gas diffusion into the surrounding medium; this rapid shrinkage may be evidence of nanoscopic bubble fragmentation. Second, upon reaching a fraction of their initial size, bubbles begin a slower shrinkage with a shrinkage rate that depends on the resting phase state of the coat lipid: fluid DMPC monolayers give a more rapid shrinkage than gel phase. DOPC-coated bubbles showed no temperature-dependent responses in the same temperature range. The results are especially interesting in that bubble compression during the pulse is likely to adiabatically heat the bubble and fluidize the coat, regardless of its initial phase state; thus, some structural feature of the resting coat, such as defect lines in the gel phase, may be important in the subsequent response to the ~3 μs ultrasound pulse.

Rapid shrinkage of lipid-coated bubbles in pulsed ultrasound

Many researchers have observed rapid shrinkage of lipid-coated microbubbles subjected to brief, MHz ultrasound pulses. The shrinkage is sometimes, but not always, accompanied by the shedding of visible fragments of the coat. It has been suggested that the shedding of the lipid coat alone is sufficient to explain the rapid shrinkage, as that loss increases bubble surface tension and, thus, internal pressure, increasing gas loss even between pulses. We have determined, however, that the shedding of the coat lipid must also entrain some of the gas content of the bubble, to account for the observed shrinkage rates. The evidence for this is that insonated bubbles typically shrink much faster than the Epstein-Plesset (diffusion) limit for gas dissolution and diffusion, whereas uncoated quiescent bubbles shrink more slowly. We have also modeled the diffusion of gas in the moving liquid surrounding the bubble and find no advective enhancement of diffusive loss of gas from the bubble. Thus, bubble gas loss through diffusion alone is insufficient to account for rapid shrinkage.

Lipid monolayer collapse and microbubble stability

Microbubbles are micrometer-size gaseous particles suspended in water, and they are often stabilized by a lipid monolayer shell. Natural microbubbles are found in freshwater and saltwater systems, and engineered microbubbles have a variety of applications in food sciences, biotechnology and medicine. Lipid-coated microbubbles are found to have remarkable stability and mechanical behavior owing to the resistance of the lipid monolayer encapsulation to collapse. The purpose of this review is to tie in recent observations of lipid-coated microbubble dissolution and gas exchange with current literature on the physics of lipid monolayer collapse in the context of lung surfactant. Based on this analysis, we conclude that microbubble shells collapse through the nucleation of microscopic folds, which then catalyze the formation and aggregation of new folds, leading to macroscopic folding events. This process results in a cyclic behavior of crumple-to-smooth transitions, which can be modulated through lipid composition. Eventually, the microbubbles stabilize at 1-2 μm diameter, regardless of initial size or lipid composition, and various mechanisms for this stabilization are postulated. Our ultimate goal is to inspire the reader to consider lipid monolayer collapse as the main long-term stabilizing mechanism for lipid-coated microbubbles, and to stimulate the use of microbubbles as a platform for studying monolayer collapse phenomena.