Optimal design and experimental investigation of surfactant encapsulated microbubbles

The diagnostic capabilities of ultrasound imaging can be improved with contrast-specific nonlinear imaging modalities such as harmonic and subharmonic imaging. The nonlinear response of an encapsulated microbubble in an acoustic field is strongly influenced by the shell viscoelastic properties that are determined by the shell composition and thickness. In this paper, the subharmonic performance of a surfactant encapsulated microbubble was optimized by choosing the appropriate composition of shell material with the aid of theoretical model. To study the effects of viscoelastic properties of microbubble shell materials on the nonlinear scattered response of microbubbles, a theoretical model-modified Herring equation for the oscillation of encapsulated microbubbles in the ultrasound field was employed. Based on this model, a computer aided design system was developed to optimize and analyze the acoustic properties, particularly subharmonic responses, of microbubbles under different shell parameters. Furthermore, surfactant encapsulated microbubbles with different viscoelastic properties were prepared by changing the shell composition. Their shell viscoelastic behavior was measured indirectly as dilational modulus of monolayer film formed with surfactant molecular. Moreover, in vitro quantitative acoustic properties measurements of these microbubbles were carried out to evaluate their subharmonic performance. Both of the theoretical simulation and acoustic measurement showed that the surfactant encapsulated microbubbles with good subharmonic properties could be designed and prepared by adjusting the shell material composition with the guide of the computer aided design system.     

Frequency relationships for ultrasonic activation of free microbubbles, encapsulated microbubbles, and gas-filled micropores.

The ultrasonic activation of free microbubbles, encapsulated microbubbles, and gas-filled micropores was explored using available linear theory. Encapsulated microbubbles, used in contrast agents for diagnostic ultrasound, have relatively high resonance frequencies and damping. At 2 MHz the resonance radii are 1.75 microns for free microbubbles, 4.0 microns for encapsulated microbubbles, and 1.84 microns for gas-filled micropores. Higher-pressure amplitudes are needed to elicit equivalent subharmonic, fundamental, or second-harmonic responses from the encapsulated microbubbles, and this behavior increases for higher frequencies. If an encapsulated microbubble becomes destabilized during exposure,the resulting liberated microbubble would be about twice the linear resonance size, which would be likely to produce subharmonic signals. Scattered signals used for medical imaging purposes may be indicative of bioeffects potential: The second harmonic signal is proportional to local shear stress for a microbubble on a boundary, and a strong subharmonic signal may imply destabilization and nucleation of free-microbubble cavitation activity. The potential for bioeffects from contrast agent gas bodies decreases rapidly with increasing frequency. This information should be valuable for understanding of the etiology of bioeffects related to contrast agents and for developing exposure indices and risk management strategies for their use in diagnostic ultrasound.     

Effects of medium viscoelasticity on bubble collapse strength of interacting polydisperse bubbles

Due to its physical and/or chemical effects, acoustic cavitation plays a crucial role in various emerging applications ranging from advanced materials to biomedicine. The cavitation bubbles usually undergo oscillatory dynamics and violent collapse within a viscoelastic medium, which are closely related to the cavitation-associated effects. However, the role of medium viscoelasticity on the cavitation dynamics has received little attention, especially for the bubble collapse strength during multi-bubble cavitation with the complex interactions between size polydisperse bubbles. In this study, modified Gilmore equations accounting for inter-bubble interactions were coupled with the Zener viscoelastic model to simulate the dynamics of multi-bubble cavitation in viscoelastic media. Results showed that the cavitation dynamics (e.g., acoustic resonant response, nonlinear oscillation behavior and bubble collapse strength) of differently-sized bubbles depend differently on the medium viscoelasticity and each bubble is affected by its neighboring bubbles to a different degree. More specifically, increasing medium viscosity drastically dampens the bubble dynamics and weakens the bubble collapse strength, while medium elasticity mainly affects the bubble resonance at which the bubble collapse strength is maximum. Differently-sized bubbles can achieve resonances and even subharmonic resonances at high driving acoustic pressures as the elasticity changes to certain values, and the resonance frequency of each bubble increases with the elasticity increasing. For the interactions between the size polydisperse bubbles, it indicated that the largest bubble generally has a dominant effect on the dynamics of smaller ones while in turn it is almost unaffected, exhibiting a pattern of destructive and constructive interactions. This study provides a valuable insight into the acoustic cavitation dynamics of multiple interacting polydisperse bubbles in viscoelastic media, which may offer a potential of controlling the medium viscoelasticity to appropriately manipulate the dynamics of multi-bubble cavitation for achieving proper cavitation effects according to the desired application.     

Subharmonic behavior of phospholipid-coated ultrasound contrast agent microbubbles

Coated microbubbles, unlike tissue are able to scatter sound subharmonically. Therefore, the subharmonic behavior of coated microbubbles can be used to enhance the contrast in ultrasound contrast imaging. Theoretically, a threshold amplitude of the driving pressure can be calculated above which subharmonic oscillations of microbubbles are initiated. Interestingly, earlier experimental studies on coated microbubbles demonstrated that the threshold for these bubbles is much lower than predicted by the traditional linear viscoelastic shell models. This paper presents an experimental study on the subharmonic behavior of differently sized individual phospholipid coated microbubbles. The radial subharmonic response of the microbubbles was recorded with the Brandaris ultra high-speed camera as a function of both the amplitude and the frequency of the driving pulse. Threshold pressures for subharmonic generation as low as 5 kPa were found near a driving frequency equal to twice the resonance frequency of the bubble. An explanation for this low threshold pressure is provided by the shell buckling model proposed by Marmottant et al. [J. Acoust. Soc. Am. 118, 3499-3505 (2005)]. It is shown that the change in the elasticity of the bubble shell as a function of bubble radius as proposed in this model, enhances the subharmonic behavior of the microbubbles.     

Resonance frequencies of lipid-shelled microbubbles in the regime of nonlinear oscillations

Knowledge of resonant frequencies of contrast microbubbles is important for the optimization of ultrasound contrast imaging and therapeutic techniques. To date, however, there are estimates of resonance frequencies of contrast microbubbles only for the regime of linear oscillation. The present paper proposes an approach for evaluating resonance frequencies of contrast agent microbubbles in the regime of nonlinear oscillation. The approach is based on the calculation of the time-averaged oscillation power of the radial bubble oscillation. The proposed procedure was verified for free bubbles in the frequency range 1-4 MHz and then applied to lipid-shelled microbubbles insonified with a single 20-cycle acoustic pulse at two values of the acoustic pressure amplitude, 100 kPa and 200 kPa, and at four frequencies: 1.5, 2.0, 2.5, and 3.0 MHz. It is shown that, as the acoustic pressure amplitude is increased, the resonance frequency of a lipid-shelled microbubble tends to decrease in comparison with its linear resonance frequency. Analysis of existing shell models reveals that models that treat the lipid shell as a linear viscoelastic solid appear may be challenged to provide the observed tendency in the behavior of the resonance frequency at increasing acoustic pressure. The conclusion is drawn that the further development of shell models could be improved by the consideration of nonlinear rheological laws.     

Ultrasonic bubbles in medicine: influence of the shell

Ultrasound contrast agents consist of microscopically small bubbles encapsulated by an elastic shell. These microbubbles oscillate upon ultrasound insonification, and demonstrate highly nonlinear behavior, ameliorating their detectability. (Potential) medical applications involving the ultrasonic disruption of contrast agent microbubble shells include release-burst imaging, localized drug delivery, and noninvasive blood pressure measurement. To develop and enhance these techniques, predicting the cracking behavior of ultrasound-insonified encapsulated microbubbles has been of importance. In this paper, we explore microbubble behavior in an ultrasound field, with special attention to the influence of the bubble shell. A bubble in a sound field can be considered a forced damped harmonic oscillator. For encapsulated microbubbles, the presence of a shell has to be taken into account. In models, an extra damping parameter and a shell stiffness parameter have been included, assuming that Hooke's Law holds for the bubble shell. At high acoustic amplitudes, disruptive phenomena have been observed, such as microbubble fragmentation and ultrasonic cracking. We analyzed the occurrence of ultrasound contrast agent fragmentation, by simulating the oscillating behavior of encapsulated microbubbles with various sizes in a harmonic acoustic field. Fragmentation occurs exclusively during the collapse phase and occurs if the kinetic energy of the collapsing microbubble is greater than the instantaneous bubble surface energy, provided that surface instabilities have grown big enough to allow for break-up. From our simulations it follows that the Blake critical radius is not a good approximation for a fragmentation threshold. We demonstrated how the phase angle differences between a damped radially oscillating bubble and an incident sound field depend on shell parameters.     

Effect of coupled oscillations on microbubble behavior

Ultrasound contrast agents are encapsulated microbubbles whose nonlinear acoustic scattering signatures have been the foundation of their use in diagnostic imaging. The coupled oscillations of microbubbles along their lines of center are investigated theoretically using radial equations in the monopole approximation and an energy balance, which is obtained for the system. Coupled microbubble pairs of different initial radii are investigated numerically relative to the normal modes for the linearized system. For microbubble pairs of different size bubbles driven below the mode of the smaller bubble and above the mode of the larger bubble, it is shown that oscillations of the smaller agent are affected substantially more by the coupling than those of the larger one. For separation distances of 10 and 500 microns, a difference of approximately 10 dB occurs in the second harmonic output of a 1.0-micron radius agent coupled with a 2.2-micron radius agent forced at 2.0 MHz and 0.3 MPa. The subharmonic spectral peak is shown to decrease approximately 19 dB for the coupling of 1.5- and 2.2-micron radius agents at 10- and 500-micron distances under the same acoustic forcing conditions. These coupling effects on the radiated pressure and its spectral power are highlighted for contrast agent imaging applications.     

Harmonic responses and cavitation activity of encapsulated microbubbles coupled with magnetic nanoparticles

Encapsulated microbubbles coupled with magnetic nanoparticles, one kind of hybrid agents that can integrate both ultrasound and magnetic resonance imaging/therapy functions, have attracted increasing interests in both research and clinic communities. However, there is a lack of comprehensive understanding of their dynamic behaviors generated in diagnostic and therapeutic applications. In the present work, a hybrid agent was synthesized by integrating superparamagnetic iron oxide nanoparticles (SPIOs) into albumin-shelled microbubbles (named as SPIO-albumin microbubbles). Then, both the stable and inertial cavitation thresholds of this hybrid agent were measured at varied SPIO concentrations and ultrasound parameters (e.g., frequency, pressure amplitude, and pulse length). The results show that, at a fixed acoustic driving frequency, both the stable and inertial cavitation thresholds of SPIO-albumin microbubble should decrease with the increasing SPIO concentration and acoustic driving pulse length. The inertial cavitation threshold of SPIO-albumin microbubbles also decreases with the raised driving frequency, while the minimum sub- and ultra-harmonic thresholds appear at twice and two thirds resonance frequency, respectively. It is also noticed that both the stable and inertial cavitation thresholds of SonoVue microbubbles are similar to those measured for hybrid microbubbles with a SPIO concentration of 114.7 μg/ml. The current work could provide better understanding on the impact of the integrated SPIOs on the dynamic responses (especially the cavitation activities) of hybrid microbubbles, and suggest the shell composition of hybrid agents should be appropriately designed to improve their clinical diagnostic and therapeutic performances of hybrid microbubble agents.     

Numerical modeling of microbubble backscatter to optimize ultrasound particle image velocimetry imaging: initial studies

We have developed a promising non-invasive ultrasound-based method for performing particle image velocimetry (PIV) in vivo. This method, termed echo PIV, provides multi-component blood velocity data with good (2 ms) temporal resolution. The method takes advantage of the non-linear ultrasound backscatter characteristics of small gas-filled microbubbles (ultrasound contrast) that are seeded into the blood stream. In this study, we use a numerical model to explore potential areas to focus future work in echo PIV.

Ultrasound backscatter from encapsulated microbubbles was modeled using a modified Rayleigh–Plesset equation (Church model, 1995), taking into account the protein/lipid shell layer as a thick, mass-conserving incompressible fluid surrounded by incompressible blood-like fluid. The equation of motion was solved numerically to characterize the fundamental and second harmonic components of the backscattered pressure. Results show a significant advantage in using the second harmonic component for echo PIV, especially for small bubble sizes less than 3 μm in diameter at 2.2 MHz frequency. The effect of the shell thickness ranging from 10 to 500 nm on the vibration amplitude of the bubble was examined and it is shown that the presence of the shell requires mechanical index (MI)>0.2 of incident pressure amplitude to improve bubble detectability. Analysis of the effect of pulse length shows a tradeoff between axial resolution (short pulse length) and bubble detectability (longer pulse length) will most likely be required. The effect of varying MI between 0.1 and 0.6 was also studied at a center frequency of 2.2 MHz and the results indicate that the resonance of the second harmonic is maximized for bubbles with diameter of approximately 2.75 μm. Bubble non-linearities at MI>0.2 induced a resonant frequency shift away from the integer multiple of the incident frequency in the second harmonic backscatter. For a given bubble size, there is a combination of optimal incident frequency and mechanical index range that maximizes the ratio of the second harmonic compared to the fundamental. This resonant frequency decreases with increasing bubble radius. Further, a narrow bandwidth pulse is shown to increase signal strength. Both of these effects may cause conflict with factors governing spatial resolution. Optimization of the incident frequency, microbubble size and mechanical index to enhance bubble detectability will depend on the particular clinical application. These theoretical predictions provide further understanding of the physics behind our echo PIV technique, and should be useful for guiding the design of echo PIV systems.     

The viscoelasticity of lipid shell and the hysteresis of subharmonic in liquid containing microbubbles

The viscoelasticity and subharmonic generation of a kind of lipid ultrasound contrast agent are investigated. Based on the measurement of the sound attenuation spectrum, the viscoelasticity of the lipid shell is estimated by use of an optimization method. Shear modulus G_S=10MPa and shear viscosity \mu _S=1.49N\cdotS/m² are obtained. The nonlinear oscillation of the encapsulated microbubble is studied with Church's model theoretically and experimentally. Especially, the dependence of subharmonic on the incident acoustic pressure is studied. The results reveal that the development of the subharmonic undergoes three stages, i.e. occurrence, growth and saturation, and that hysteresis appears in descending ramp insonation.     

Modeling subharmonic response from contrast microbubbles as a function of ambient static pressure

Variation of subharmonic response from contrast microbubbles with ambient pressure is numerically investigated for non-invasive monitoring of organ-level blood pressure. Previously, several contrast microbubbles both in vitro and in vivo registered approximately linear (5-15 dB) subharmonic response reduction with 188 mm Hg change in ambient pressure. In contrast, simulated subharmonic response from a single microbubble is seen here to either increase or decrease with ambient pressure. This is shown using the code BUBBLESIM for encapsulated microbubbles, and then the underlying dynamics is investigated using a free bubble model. The ratio of the excitation frequency to the natural frequency of the bubble is the determining parameter--increasing ambient pressure increases natural frequency thereby changing this ratio. For frequency ratio below a lower critical value, increasing ambient pressure monotonically decreases subharmonic response. Above an upper critical value of the same ratio, increasing ambient pressure increases subharmonic response; in between, the subharmonic variation is non-monotonic. The precise values of frequency ratio for these three different trends depend on bubble radius and excitation amplitude. The modeled increase or decrease of subharmonic with ambient pressure, when one happens, is approximately linear only for certain range of excitation levels. Possible reasons for discrepancies between model and previous experiments are discussed.     

超声造影剂微气泡的包膜黏弹特性的定量表征研究

包膜黏弹特性显著影响微气泡超声造影剂的诊断及治疗应用效果. 本文结合原子力显微镜技术及声衰减特性测量提出了一种对微气泡造影剂包膜黏弹特性定量表征的新方法. 首先采用原子力显微镜技术进行机械特性分析得到包膜微气泡的有效硬度及体弹性模量; 然后测量声衰减特性, 基于微气泡动力学理论, 计算包膜微气泡的体黏度系数. 为验证方法的有效性, 实验制备了直径为1-5 μm的白蛋白包膜微气泡造影剂, 原子力显微镜测量的有效硬度和体弹性模量分别为0.149±0.012 N/m和8.31±0.667 MPa, 并与粒径无关. 声衰减特性测量和动力学理论拟合的包膜微气泡的体黏度系数为0.374±0.003 Pa·s. 该方法可推广至其他种类包膜微气泡的黏弹特性表征, 对超声造影剂的制备及其诊断和治疗应用有积极意义.     

超声造影剂微泡非线性动力学响应的机理及相关应用

随着生命科学及现代医学的发展, 一体化无创精准诊疗已经日益成为人们关注的焦点问题, 而关于超声造影剂微泡的非线性效应的相关机理、动力学建模及其在超声医学领域中的应用研究也得到了极大的推动. 本文对下列课题进行了总结和讨论, 包括: 1)基于Mie散射技术和流式细胞仪对造影剂微泡参数进行定征的一体化解决方案; 2)通过对微泡包膜的黏弹特性进行非线性修正, 构建新的包膜微泡动力学模型; 3)探索造影剂惯性空化阈值与其包膜参数之间的相关性; 以及4)研究超声联合造影剂微泡促进基因/药物转染效率并有效降低其生物毒性的相关机理.     

多层膜磁性微泡的非线性声振动特性

超顺磁性氧化铁纳米粒子与造影剂微泡结合形成磁性微泡,用于产生多模态造影剂,以增强医学超声和磁共振成像.将装载有纳米磁性颗粒的微泡包膜层看作由磁流体膜与磷脂膜组合而成的双层膜结构,同时考虑磁性纳米颗粒体积分数a对膜密度及黏度的影响,从气泡动力学基本理论出发,构建多层膜结构磁性微泡非线性动力学方程.数值分析了驱动声压和频率等声场参数、颗粒体积分数、膜层厚度以及表面张力等膜壳参数对微泡声动力学行为的影响.结果表明,当磁性颗粒体积分数较小且a≤0.1时,磁性微泡声响应特性与普通包膜微泡相似,微泡的声频响应与其初始尺寸和驱动压有关;当驱动声场频率f为磁性微泡共振频率f0的2倍(f=2f0)时,微泡振动失稳临界声压最低;磁性颗粒的存在抑制了泡的膨胀和收缩但抑制效果非常有限;磁性微泡外膜层材料的表面张力参数K及膜层厚度d也会影响微泡的振动,当表面张力参数及膜厚取值分别为0.2—0.4 N/m及50—150 nm时,可观察到气泡存在不稳定振动响应区.     

Difference-frequency ultrasound generation from microbubbles under dual-frequency excitation

The difference-frequency (DF) ultrasound generated by using parametric effect promises to improve detection depth owing to its low attenuation, which is beneficial for deep tissue imaging. With ultrasound contrast agents infusion, the harmonic components scattered from the microbubbles, including DF, can be generated due to the nonlinear vibration. A theoretical study on the DF generation from microbubbles under the dual-frequency excitation is proposed in formula based on the solution of the RPNNP equation. The optimisation of the DF generation is discussed associated with the applied acoustic pressure, frequency, and the microbubble size. Experiments are performed to validate the theoretical predictions by using a dual-frequency signal to excite microbubbles. Both the numerical and experimental results demonstrate that the optimised DF ultrasound can be achieved as the difference frequency is close to the resonance frequency of the microbubble and improve the contrast-to-tissue ratio in imaging.     

Resonance frequency of microbubbles: effect of viscosity

The transmitted frequency at which a gas bubble of millimeter or submillimeter size oscillates resonantly in a low-viscosity liquid is approximately equal to the undamped natural frequency (referred to as the Minnaert frequency if surface tension effects are disregarded). Based on a theoretical analysis of bubble oscillation, this paper shows that such an approximation cannot be validated for microbubbles used in contrast-enhanced ultrasound imaging. The contrast-agent microbubbles represent either encapsulated bubbles of size less than 10 microm or free (nonencapsulated) bubbles of submicron size. The resonance frequency of the microbubbles deviates significantly from the undamped natural frequency over the whole range of microbubble sizes due to the increased viscous damping coefficient. The difference between these two frequencies is shown to have a tremendous impact on the resonant backscatter by the microbubbles. In particular, the first and second harmonics of the backscattered signal from the microbubbles are characterized by their own resonance frequencies, equal to neither the microbubble resonance frequency nor the undamped natural frequency.     

Using light scattering to measure the response of individual ultrasound contrast microbubbles subjected to pulsed ultrasound in vitro

Light scattering was used to measure the radial pulsations of individual ultrasound contrast microbubbles subjected to pulsed ultrasound. Highly diluted Optison or Sonazoid microbubbles were injected into either a water bath or an aqueous solution containing small quantities of xanthan gum. Individual microbubbles were insonified by ultrasound pulses from either a commercial diagnostic ultrasound machine or a single element transducer. The instantaneous response curves of the microbubbles were measured. Linear and nonlinear microbubble oscillations were observed. Good agreement was obtained by fitting a bubble dynamics model to the data. The pulse-to-pulse evolution of individual microbubbles was investigated, the results of which suggest that the shell can be semipermeable, and possibly weaken with subsequent pulses. There is a high potential that light scattering can be used to optimize diagnostic ultrasound techniques, understand microbubble evolution, and obtain specific information about shell parameters.     

Microstreaming velocity field and shear stress created by an oscillating encapsulated microbubble near a cell membrane

Sonoporation mediated by microbubbles is being extensively studied as a promising technology to facilitate gene/drug delivery to cells. However, the theoretical study regarding the mechanisms involved in sonoporation is still in its infancy. Microstreaming generated by pulsating microbubble near the cell membrane is regarded as one of the most important mechanisms in the sonoporation process. Here, based on an encapsulated microbubble dynamic model with considering nonlinear rheological effects of both shell elasticity and viscosity, the microstreaming velocity field and shear stress generated by an oscillating microbubble near the cell membrane are theoretically simulated. Some factors that might affect the behaviors of microstreaming are thoroughly investigated, including the distance between the bubble center and cell membrane (d), shell elasticity (χ), and shell viscosity (κ). The results show that (i) the presence of cell membrane can result in asymmetric microstreaming velocity field, while the constrained effect of the membrane wall decays with increasing the bubble-cell distance; (ii) the bubble resonance frequency increases with the increase in d and χ, and the decrease in κ, although it is more dominated by the variation of shell elasticity; and (iii) the maximal microstreaming shear stress on the cell membrane increases rapidly with reducing the d, χ, and κ. The results suggest that microbubbles with softer and less viscous shell materials might be preferred to achieve more efficient sonoporation outcomes, and it is better to have bubbles located in the immediate vicinity of the cell membrane.     

超声微泡低频动力学特性及其生物效应

理论上利用有耗散函数的Lagrange方程,建立了有壳微泡的R(t)运动方程,开展了自由空间中有壳微泡动力学特性的研究,表明微泡内外半径增量随声压的增大、超声频率的降低、初始内径的增大及壳厚的减薄而迅速增大。实验上,利用Mie散射技术在80°散射角和前向Mie散射检测新技术实验测量了微泡R(t)曲线;利用体视显微镜,实时观察了超声微泡对动物活体微血管损伤,开展了超声微泡生物效应的动物和细胞试验研究。结果表明:(1)超声作用下,微泡引起肿瘤中微血管壁周期性膨胀收缩而发生管壁破裂,形成血栓和微血管栓塞,抑制了肿瘤生长;(2)超声联合微泡可以破坏微血管内皮生长因子(VEGF)和肝癌细胞,可以减少肿瘤血管和癌细胞再生,因此,低频超声联合微泡技术是一种值得探索抑制肿瘤生长的新技术。     

Nonlinear response to ultrasound of encapsulated microbubbles

The acoustic backscatter of encapsulated gas-filled microbubbles immersed in a weak compressible liquid and irradiated by ultrasound fields of moderate to high pressure amplitudes is investigated theoretically. The problem is formulated by considering, for the viscoelastic shell of finite thickness, an isotropic hyperelastic neo-Hookean model for the elastic contribution in addition to a Newtonian viscous component. First and second harmonic scattering cross-sections have been evaluated and the quantitative influence of the driving pressure amplitude on the harmonic resonance frequencies for different initial equilibrium bubble sizes and for different encapsulating physical properties has been determined. Conditions for optimal second harmonic imaging have been also investigated and some regions in the parameters space where the second harmonic intensity is dominant over the fundamental have been identified. Results have been obtained for albumin, lipid and polymer encapsulating shells, respectively.