Mechanism of Membrane Poration by Shock Wave Induced Nanobubble Collapse: A Molecular Dynamics Study

Characterization of the Dynamic Behavior of Multinanobubble System under Shock Wave Influence

Shockwave-induced changes in nanobubbles cause cavitation erosion and membrane damage but can also be applied to biocarrier transport. Currently, research focuses on single nanobubbles; however, in reality, nanobubbles usually appear as a multibubble system. Therefore, this study proposes a method based on cutting and replicating to construct a multibubble model. This method can be widely applied to molecular dynamics (MD) models and enhance the customization capabilities of MD models. The dynamic behavior of a multinanobubble system with different numbers and arrangements of nanobubbles is investigated with the MD method under the influence of shock waves in a liquid argon system. The study also explores the range of influence between nanobubbles. The results show that in the case of two nanobubbles, when the distance between the bubbles is constant, the smaller the angle between the direction of the shock wave and the line connecting the bubbles, the greater is the influence between nanobubbles, and the moment of collapse of the nanobubbles farther away from the shock wave is slower. When three nanobubbles are arranged with a right offset, after the first bubble collapses, the effect on the other two bubbles is similar to the changes in bubbles when the angle of arrangement is 30° or 60°. Under a different arrangement, the change of shock wave velocity on the nanobubble size only affects its collapse time and contraction collapse rate. When the shock wave with a radian of about 2.87 or greater than 2.87 touches the bubbles, the collapse of the second nanobubble will not be affected.

The impact of low-velocity shock waves on the dynamic behaviour characteristics of nanobubbles

Low-velocity shock wave-induced contraction and expansion of nanobubbles can be applied to biocarriers and microfluidic systems. Although experiments have been conducted to study the application effects, the dynamic behavior characteristics of nanobubbles remain unexplored. In this work, we utilize molecular dynamics (MD) simulations to investigate the dynamic behavior characteristics of nanobubbles influenced by low-velocity shock waves in a liquid argon system. The DBSCAN (Density-Based Spatial Clustering of Applications with Noise) machine learning method is used to calculate the equivalent radius of nanobubbles. Two statistical methods are then utilized to predict the time series changes in the equivalent radius of nanobubbles without rebound shock waves. The piston velocity is analyzed using the bisection method to obtain the critical impact states of the nanobubble. The results show that at the low velocity shock wave (piston velocity of 0.1 km s-1), the shock wave pressure is small, the non-vacuum nanobubbles contract and expand in a circular shape, and the gas particles inside the bubble are not dispersed. In contrast, the vacuum nanobubbles collapse directly. As the shock wave rebounds upon impact, it triggers periodic contraction and expansion of the nanobubbles. The predictions indicate that the equivalent radius will vary within a small range according to the pre-predicted values in the absence of the rebound shock wave. Nanobubbles are present in four critical impact states: dispersed gaps, multiple smaller bubbles, two split bubbles, and a concave bubble.

Effect of Shock-Wave-Mediated Collapse on Nanobubble Characteristics

To enhance the comprehension of the cavitation mechanism and explore its practical use in industrial production, this study developed models involving oxygen, varying bubble radii, and bubble quantities. This study uses molecular dynamics simulations coupled with the momentum mirror method to examine the collapse characteristics of bubbles during ultrasonic cavitation. The investigation uncovers patterns in the fluctuation of the maximum local density of water molecules, released pressure, and temperature. The findings demonstrate that, when oxygen-containing bubbles collapse at identical radii, the local density is notably higher and diminishes more rapidly. Moreover, the changes in the shape exhibit greater regularity. During the bubble collapse, a depression forms on the bubble's surface, coinciding with a notable surge in local density around the depression. As bubble radii and quantities increase, so does the local density along with a concurrent rise in the maximum pressure. Intriguingly, the model demonstrates the lowest pressure at Z = 35 Å accompanied by the emergence of a small crescent-shaped region with a reduced density. Throughout the pressure ascension phase, the rate of the maximum pressure change escalates with an increase in the number of bubbles. Conversely, during the pressure descent phase, the rate of the maximum pressure change diminishes with a growing number of bubbles. However, it is important to note that the maximum pressure does not exhibit a direct correlation with the number of bubbles. Ultimately, this study provides valuable technical guidance and a theoretical foundation for the integration of ultrasonic cavitation in industrial production processes.

Molecular dynamics of the spontaneous generation mechanism of natural gas hydrates during methane nanobubble rupture

Natural gas hydrates have garnered significant attention as a potential new source of alternative energy, and understanding their formation mechanism is of paramount importance for efficient utilization and pipeline transportation. However, there is no consensus among academics on the formation mechanism of natural gas hydrates. In this paper, we propose a method for promoting the rapid formation of natural gas hydrates based on the rupture of methane nanobubbles, which creates local high temperature and pressure to facilitate the mixing of methane and water. The rapid decrease in system temperature and pressure during the process further enhances the formation of gas hydrates. Using molecular dynamics simulations, we theoretically verify the formation of natural gas hydrates. Our results indicate that the instantaneous rupture of methane nanobubbles induced by shock waves leads to a dramatic increase in the local molecular motion velocity around the bubbles. This results in extreme local high temperature and high pressure, leading to complete mixing of methane and water and rapid formation of gas hydrates during the cooling and pressure drop of the mixture. We confirm our findings by analyzing F3-order parameters, F4-order parameters, and water cage statistics.

Molecular modelling of shockwave-mediated delivery of paclitaxel aggregates across the neuronal plasma membrane

Shock-assisted paclitaxel (PTX) transport across the blood-brain barrier offers a promising treatment strategy for brain tumors. Here, based on a realistically complex human brain plasma membrane (PM) model, we investigated the dynamic transmembrane behavior of a PTX cluster by shock induced bubble collapse, focusing on the effect of impulse (I), bubble diameter (D) and arrays. The results show that all three factors can control the transport depth (ΔDPM) of PTX. For a fixed D, the ΔDPM grows exponentially with I, ΔDPM ∼ exp (I), and eventually reaches a critical depth. But the depth, ΔDPM, can be adjusted linearly in a wider range of D. This mainly depends on the size of jets from bubble collapse. For bubble arrays, the bubbles in series can transport PTX deeper than a single bubble, while the parallel does the opposite. In addition, only PTX clusters in the range of jet action can be successfully transported. Finally, the absorption of PTX clusters was examined via recovery simulation. Not all PTX clusters across the membrane can be effectively absorbed by cells. The shallow PTX clusters are quickly attracted by the membrane and embedded into it. The critical depth at which PTX clusters can be effectively absorbed is about 20 nm. These molecular-level mechanisms and dynamic processes of PTX clusters crossing the PM membrane may be helpful in optimizing the application of shock-induced bubble collapse for the delivery of PTX to tumor cells.

A Deep Potential model for liquid-vapor equilibrium and cavitation rates of water

Computational studies of liquid water and its phase transition into vapor have traditionally been performed using classical water models. Here, we utilize the Deep Potential methodology-a machine learning approach-to study this ubiquitous phase transition, starting from the phase diagram in the liquid-vapor coexistence regime. The machine learning model is trained on ab initio energies and forces based on the SCAN density functional, which has been previously shown to reproduce solid phases and other properties of water. Here, we compute the surface tension, saturation pressure, and enthalpy of vaporization for a range of temperatures spanning from 300 to 600 K and evaluate the Deep Potential model performance against experimental results and the semiempirical TIP4P/2005 classical model. Moreover, by employing the seeding technique, we evaluate the free energy barrier and nucleation rate at negative pressures for the isotherm of 296.4 K. We find that the nucleation rates obtained from the Deep Potential model deviate from those computed for the TIP4P/2005 water model due to an underestimation in the surface tension from the Deep Potential model. From analysis of the seeding simulations, we also evaluate the Tolman length for the Deep Potential water model, which is (0.091 ± 0.008) nm at 296.4 K. Finally, we identify that water molecules display a preferential orientation in the liquid-vapor interface, in which H atoms tend to point toward the vapor phase to maximize the enthalpic gain of interfacial molecules. We find that this behavior is more pronounced for planar interfaces than for the curved interfaces in bubbles. This work represents the first application of Deep Potential models to the study of liquid-vapor coexistence and water cavitation.

Cell-cycle dependent nuclear gene delivery enhances the effects of E-cadherin against tumor invasion and metastasis

Gene delivery is the process by which foreign DNA is transferred to host cells, released from intracellular vesicles, and transported to the nuclei for transcription. This process is frequently inefficient and difficult to control spatiotemporally. We developed a gene delivery strategy that uses ultrasound to directly deliver plasmid DNA into nuclei via gas vesicles (GVs)-based intracellular cavitation. pDNA-binding GVs can be taken up by cells and cause intracellular cavitation when exposed to acoustic irradiation and delivering their pDNA payloads into nuclei. Importantly, GVs can remain stable in the cytoplasm in the absence of acoustic irradiation, allowing for temporally controlled nuclear gene delivery. We were able to achieve spatiotemporal control of E-cadherin nuclear gene delivery in this manner, demonstrating its efficacy in tumor invasion and metastasis inhibition. Interestingly, we discovered that nuclear gene delivery of E-cadherin during the G2/M phase of the cell cycle in C6 tumor cells inhibited tumor invasion and metastasis more effectively than during the G1 and S phases. The gene delivery of E-cadherin at the G2/M phase resulted in significantly lower expression of Fam50a, which reduced Fam50a/Runx2 interaction and led to reduced transactivation of MMP13, an important factor for epithelial-mesenchymal transition, as observed in a molecular mechanism assay. Thus, using remote acoustic control of intracellular cavitation of pDNA-GVs, we developed a high spatiotemporally controllable gene delivery strategy and achieved stronger tumor invasion and metastasis inhibition effects by delivering the E-cadherin gene at the G2/M phase.

Anomalous Vapor and Ice Nucleation in Water at Negative Pressures: A Classical Density Functional Theory Study

In contrast to the abundance of work on the anomalous behavior of water, the relationship between the water's thermodynamic anomalies and kinetics of phase transition from metastable water is relatively unexplored. In this work, we have employed classical density functional theory to provide a unified and coherent picture of nucleation (both vapor and ice) from metastable water at negative pressure conditions. Our results suggest a peculiar nonmonotonic temperature dependence of vapor-liquid surface tension at temperatures where vapor-liquid coexistence is metastable with respect to the ice phase. The vapor nucleation barrier on isochoric cooling also shows a nonmonotonic temperature dependence. We further report that, for low density isochores, the temperature of the minimum vapor nucleation barrier (TΔΩv/min*) does not coincide with the temperature of maximum density (TMD) where metastability is maximum. The difference between the TΔΩv/min* and the TMD, however, decreases with increasing the density of the isochore. The vapor nucleation barrier along isobars shows an interesting crossover behavior in the vicinity of the Widom line on lowering the temperature. Our results on the ice nucleation suggest an anomalous retracing behavior of the nucleation barrier along isotherms at negative pressures and theoretically validate the recent findings that the reentrant ice(Ih)-liquid coexistence line can induce a drastic change in the kinetics of ice nucleation. Thus, this study establishes a direct connection between the metastable water's thermodynamic anomalies and the (vapor and ice) nucleation kinetics. In addition, this study provides deeper insights into the origin of the isothermal compressibility maximum on isochoric cooling.

Focal opening of the neuronal plasma membrane by shock-induced bubble collapse for drug delivery: a coarse-grained molecular dynamics simulation

Cell permeabilization using shock-induced bubble collapse provides an attractive choice for drug delivery systems. In this work, based on a realistically human brain plasma membrane (PM) model, we investigated the focal opening of this complex model by the jet from cavitation, focusing on the effect of characteristic membrane components, particle velocity (u_p) and bubble diameters (D). Both high levels of cholesterol and specific cerebrosides in the PM model limit the pore opening of cavitation jets. Sphingomyelin is the opposite, but has little effect due to its low content. Two adjustable parameters of u_p and D can be coupled to control the opening size. The relationship between them and the maximum pore area was provided for the first time. The maximum pore area increases with the u_p (or the impulse that is positively related to u_p) in the low-speed range, which agrees with the experimentally observed impulse determinism. However, the maximum area drops in the high-speed range. Combined with D, we proposed that the jet size determines the pore size, not the impulse. Larger bubbles that can create a larger pore in the membrane have a larger jet size, but their impulse is relatively small. Finally, the recovery simulation shows that the membrane with a small pore can be quickly recovered within 300 ps, while that with a larger pore did not recover until 2 μs. These rules from this work may be helpful to optimize the choice of shock waves for the delivery of different drugs across membranes.

Molecular dynamics simulation of cancer cell membrane perforated by shockwave induced bubble collapse

It has been widely accepted that cancer cells are softer than their normal counterparts. This motivates us to propose, as a proof-of-concept, a method for the efficient delivery of therapeutic agents into cancer cells, while normal cells are less affected. The basic idea of this method is to use a water jet generated by the collapse of the bubble under shockwaves to perforate pores in the cell membrane. Given a combination of shockwave and bubble parameters, the cancer membrane is more susceptible to bending, stretching, and perforating than the normal membrane because the bending modulus of the cancer cell membrane is smaller than that of the normal cell membrane. Therefore, the therapeutic agent delivery into cancer cells is easier than in normal cells. Adopting two well-studied models of the normal and cancer membranes, we perform shockwave induced bubble collapse molecular dynamics simulations to investigate the difference in the response of two membranes over a range of shockwave impulse 15-30 mPa s and bubble diameter 4-10 nm. The simulation shows that the presence of bubbles is essential for generating a water jet, which is required for perforation; otherwise, pores are not formed. Given a set of shockwave impulse and bubble parameters, the pore area in the cancer membrane is always larger than that in the normal membrane. However, a too strong shockwave and/or too large bubble results in too fast disruption of membranes, and pore areas are similar between two membrane types. The pore closure time in the cancer membrane is slower than that in the normal membrane. The implications of our results for applications in real cells are discussed in some details. Our simulation may be useful for encouraging future experimental work on novel approaches for cancer treatment.

Light triggered nanoscale biolistics for efficient intracellular delivery of functional macromolecules in mammalian cells

Biolistic intracellular delivery of functional macromolecules makes use of dense microparticles which are ballistically fired onto cells with a pressurized gun. While it has been used to transfect plant cells, its application to mammalian cells has met with limited success mainly due to high toxicity. Here we present a more refined nanotechnological approach to biolistic delivery with light-triggered self-assembled nanobombs (NBs) that consist of a photothermal core particle surrounded by smaller nanoprojectiles. Upon irradiation with pulsed laser light, fast heating of the core particle results in vapor bubble formation, which propels the nanoprojectiles through the cell membrane of nearby cells. We show successful transfection of both adherent and non-adherent cells with mRNA and pDNA, outperforming electroporation as the most used physical transfection technology by a factor of 5.5–7.6 in transfection yield. With a throughput of 10⁴-10⁵ cells per second, biolistic delivery with NBs offers scalable and highly efficient transfections of mammalian cells.

Ballistic delivery with micro/nano-particles has been successfully used to transfect plant cells, however, has failed in mammalian cells due to toxic effects. Here, the authors report on a self-assembled nano-ballistic delivery system for the delivery of functional macromolecules and demonstrate efficient transfection of mammalian cells.

A Novel Gating Mechanism of Aquaporin-4 Water Channel Mediated by Blast Shockwaves for Brain Edema

As the principal water channel in the brain, aquaporin-4 (AQP4) plays a vital role in brain edema, but its role in blast brain edema is unclear. On the basis of molecular simulations, we reveal the atomically detailed picture of AQP4 in response to blast shockwaves. The results show that the shockwave alone closes the AQP4 channel; however, shock-induced bubble collapse opens it. The jet from bubble collapse forcefully increases the distance between helices and the tilt angles of six helices relative to the membrane vertical direction in a very short time. The average channel size increases about 2.6 times, and the water flux rate is nearly 20 times higher than for normal states. It is responsible for abnormal water transport and a potential cause of acute blast brain edema. Additionally, the open AQP4 channel quickly returns to its normal state, which is in turn helpful for edema absorption. Thus, a novel gating mechanism for AQP4 related to the secondary structure change has been provided, which is different from the previous residue-mediated gating mechanism.

Molecular Dynamics Simulations of Shockwave Affected STMV Virus to Measure the Frequencies of the Oscillatory Response

Acoustic shockwaves are of interest as a possible means of the selective inactivation of viruses. It has been proposed that such inactivation may be enhanced by driving the virus particles at frequencies matching the characteristic frequency corresponding to acoustic modes of the viral structures, setting up a resonant response. Characteristic frequencies of viruses have been previously studied through opto-mechanical techniques. In contrast to optical excitation, shockwaves may be able to probe acoustic modes without the limitation of optical selection rules. This work explores molecular dynamics simulations of shockwaves interacting with a single STMV virus structure, in full atomistic detail, in order to measure the frequency of the response of the overall structure. Shockwaves of varying energy were set up in a water box containing the STMV structure by assigning water molecules at the edge of the box with an elevated velocity inward—in the direction of the virus. It was found that the structure compressed and stretched in a periodic oscillation of frequency 65 ± 6.5 GHz. This measured frequency did not show strong dependency on the energy of the shockwave perturbing the structure, suggesting the frequency is a characteristic of the structure. The measured frequency is also consistent with values predicted from elastic theory. Additionally, it was found that subjecting the virus to repeated shockwaves led to further deformation of the structure and the magnitude of the overall deformation could be altered by varying the time delay between repeated shockwave pulses.

Generation and characterization of nanobubbles in ionic liquid for a green extraction of polyphenols from Carya cathayensis Sarg

Nanobubbles (NBs) generated-nanojets membrane poration have gained enormous attention. In this study, NBs were fabricated as a novel green approach to assist ionic liquid (IL) [C₄C₁im][BF₄] extraction of polyphenols from Carya cathayensis Sarg. husk. NBs were successfully generated with mean size of 85.47 ± 5 nm, zeta potential of +39 ± 2.24 mV, and concentration of 21.15 ± 0.75 × 10⁸ particles/mL (stable for over 48 h in IL solution). Compared to common solutions extract, IL-NBs extract showed significantly higher (p < 0.05) antioxidant activity and polyphenols yields with a total polyphenol, total flavonoid, and total tannins contents of 85.67 ± 2.05 mg GAE/g DW, 42.44 ± 1.17 mg CE/g DW, and 8.2 ± 0.05 mg TAE/g DW, respectively. The SEM results confirmed that NBs' nanojets caused morphological destruction of the husk powder. Overall, IL-NBs solution showed better extraction efficiency of polyphenols than other solutions, giving insight into a new "green" nanotechnology-based extraction method.

Shock-Induced Damage Mechanism of Perineuronal Nets

The perineuronal net (PNN) region of the brain’s extracellular matrix (ECM) surrounds the neural networks within the brain tissue. The PNN is a protective net-like structure regulating neuronal activity such as neurotransmission, charge balance, and action potential generation. Shock-induced damage of this essential component may lead to neuronal cell death and neurodegenerations. The shock generated during a vehicle accident, fall, or improvised device explosion may produce sufficient energy to damage the structure of the PNN. The goal is to investigate the mechanics of the PNN in reaction to shock loading and to understand the mechanical properties of different PNN components such as glycan, GAG, and protein. In this study, we evaluated the mechanical strength of PNN molecules and the interfacial strength between the PNN components. Afterward, we assessed the PNN molecules’ damage efficiency under various conditions such as shock speed, preexisting bubble, and boundary conditions. The secondary structure altercation of the protein molecules of the PNN was analyzed to evaluate damage intensity under varying shock speeds. At a higher shock speed, damage intensity is more elevated, and hyaluronan (glycan molecule) is most likely to break at the rigid junction. The primary structure of the protein molecules is least likely to fail. Instead, the molecules’ secondary bonds will be altered. Our study suggests that the number of hydrogen bonds during the shock wave propagation is reduced, which leads to the change in protein conformations and damage within the PNN structure. As such, we found a direct connection between shock wave intensity and PNN damage.

Impact of Lipid Peroxidation on the Response of Cell Membranes to High-Speed Equibiaxial Stretching: A Computational Study

The difference between diseased and healthy cellular membranes in response to mechanical stresses is crucial for biology, as well as in the development of medical devices. However, the biomolecular mechanisms by which mechanical stresses interact with diseased cellular components remain largely unknown. In this work, we focus on the response of diseased cellular membranes with lipid peroxidation to high-speed tensile loadings. We find that the critical areal strain (ξ_c, when the pore forms) is highly sensitive to lipid peroxidation. For example, ξ_c of a fully oxidized bilayer is only 64 and 69% of the nonoxidized one at the stretching speed of 0.1 and 0.6 m/s, respectively. ξ_c decreases with the increase in the oxidized lipid ratio, regardless of the speeds. Also, the critical rupture tension of membranes exhibits a similar change. It is obvious that the oxidized membranes are more easily damaged than normal ones by high-speed stretching, which coincides with experimental findings. The reason is that peroxidation introduces a polar group to the tail of lipids, increases the hydrophilicity of tails, and warps the tails to the membrane-water interface, which causes loose accumulation and disorder of lipid tails. This can be deduced from the variation in the area per lipid and order parameter. In addition, the lowering stretching modulus and line tension of membranes (i.e., softening) after lipid peroxidation is also a significant factor. We reveal the difference between the peroxidized (diseased) and normal membrane in response to high-speed stretching, give the ξ_c value in the pore formation of membranes and analyze the influence of the stretching speed, peroxidation ratio, and molecular structure of phospholipids. We hope that the molecular-level information will be useful for the development of biological and medical devices in the future.

Stimuli-responsive nanobubbles for biomedical applications

Stimuli-responsive nanobubbles have received increased attention for their application in spatial and temporal resolution of diagnostic techniques and therapies, particularly in multiple imaging methods, and they thus have significant potential for applications in the field of biomedicine. This review presents an overview of the recent advances in the development of stimuli-responsive nanobubbles and their novel applications. Properties of both internal- and external-stimuli responsive nanobubbles are highlighted and discussed considering the potential features required for biomedical applications. Furthermore, the methods used for synthesis and characterization of nanobubbles are outlined. Finally, novel biomedical applications are proposed alongside the advantages and shortcomings inherent to stimuli-responsive nanobubbles.

Shock-induced collapse of surface nanobubbles

The collapse of cavitation bubbles often releases high-speed liquid jets capable of surface damage, with applications in drug delivery, cancer treatment, and surface cleaning. Spherical cap-shaped surface nanobubbles have previously been found to exist on immersed substrates. Despite being known nucleation sites for cavitation, their collapsing dynamics are currently unexplored. Here, we use molecular dynamics simulations to model the shock-induced collapse of different surface nanobubble sizes and contact angles. Comparisons are made with additional collapsing spherical nanobubble simulations near a substrate, to investigate the differences in their jet formation and resulting substrate pitting damage. Our main finding is that the pitting damage in the surface nanobubble simulations is greatly reduced, when compared to the spherical nanobubbles, which is primarily caused by the weaker jets formed during their collapse. Furthermore, the pit depths for surface nanobubble collapse do not depend on bubble size, unlike in the spherical nanobubble cases, but instead depend only on their contact angle. We also find a linear scaling relationship for all bubble cases between the final substrate damage and the peak pressure impulse at the impact centre, which can now be exploited to assess the relative damage in other computational studies of collapsing bubbles. We anticipate the more controlled surface-damage features produced by surface nanobubble cavitation jets will open up new applications in advanced manufacturing, medicine, and precision cleaning.

Molecular Mechanism of Ultrasound-Induced Structural Defects in Liposomes: A Nonequilibrium Molecular Dynamics Simulation Study

The use of ultrasound in combination with liposomes is a promising approach to improve drug delivery. To achieve an optimal drug release rate, it is important to understand how ultrasound induces pathways on the liposome surface where drugs can be released from the liposome. To this end, we carry out large-scale ultrasound-induced molecular dynamics simulations for three single lipid component liposomes formed from the commonly used phospholipids: 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC), 1,2-dipalmitoylphosphatidylcholine (DPPC), or phosphatidylcholine (POPC). The results show that ultrasound induces the detachment of two leaflets of the DOPC surface, suggesting that the drug release pathway may be through the low lipid packing areas on the stretched surface. In contrast, ultrasound induces pore formation on the surface of DPPC and DOPC, where drugs could escape from the liposomes. While the leaflet detachment and transient pore formation are the mechanisms of DOPC and DPPC, respectively, in both liquid-ordered and liquid-disordered phases, the leaflet detachment mechanism is switched to the transient pore formation mechanism on going from the liquid-ordered phase to the liquid-disordered phase in the POPC liposome. By adding 30% mol cholesterol, the leaflet detachment mechanism is observed in all liposomes. We found that the molecular origin that determines a mechanism is the competition between the intraleaflet and interleaflet interacting energy of lipids. The connection to experimental and theoretical modeling is discussed in some detail.

Impact of Shock-Induced Cavitation Bubble Collapse on the Damage of Cell Membranes with Different Lipid Peroxidation Levels

Although the interaction mechanism between shock waves and cells is critical for advancing the medical applications of shock waves, we still have little understanding about it. This work aims to study the response of diseased cells subjected to lipid peroxidation to the nanojet from shock wave-induced bubble collapse by using the coarse-grained molecular dynamics simulation. Factors considered in the simulations include the shock velocity (up), movement time of piston (τp), bubble size (R), and peroxidation level of membranes. Here, we mainly focus on the role of peroxidation levels, that is, the degree (%) and the distribution of oxidized lipids in membranes. The results indicate that the shock damage threshold (up at which the pore in membranes is formed) of peroxidation membranes is less than that of normal membranes and decreases with the peroxidation degree. Importantly, the distribution of oxidized lipids has more effect on the damage threshold than the peroxidation degree. The threshold of membrane with 33% localized oxidized lipids is lower than that of membrane with 50% average oxidized lipids. The results can be explained by the stretching modulus (κs) and bending modulus (κb) of cell membranes. For example, the κb value (4.3 × 10-20 J) of 100% peroxidation membrane is about half of that (8.4 × 10-20 J) of a membrane without peroxidation. A lower modulus means high deformation under the same impact. Further analysis shows that peroxidation introduces a polar hydrophobic group to the tail of phospholipids that increases the hydrophilicity of tails and warps the tail of phospholipids toward the membrane-water interface, resulting in looser accumulation. This can be confirmed by the increased average phospholipid area with peroxidation levels. Indeed, most of the pores formed during the shock can heal. However, the permeation of water molecules across the healing membrane still increased. All these membrane-level information obtained from this study will be useful for improving the biomedical applications of shock waves.

Molecular mechanism of ultrasound interaction with a blood brain barrier model

The brain is strictly protected by the blood brain barrier preventing the crossing of therapeutics to treat brain diseases. The high and low intensity focused ultrasound methods have been used to temporarily open the blood brain barrier, facilitating the transport of drugs. The methods are very promising because the opening is transient, localized, and noninvasive. However, the molecular mechanism of the opening is unknown, and this limits the development and application of these methods. With this in mind, we carry out a molecular dynamics simulation study to understand the interaction of ultrasound with the cell membrane and the tight junction. Our minimal blood brain barrier model is composed of two lipid bilayers, mimicking two portions of neighboring cells, connected together by a tight junction formed by a pair of two cis-dimers of the claudin-5 protein. Using an experimental ultrasound frequency of 50 MHz, simulations show that at low intensities, ultrasound does not impact the structure of the cell membranes and tight junction, implying that the direct interaction of ultrasound with the blood brain barrier is not responsible for the experimentally observed opening. At high intensities, the ultrasound pulls the monolayers of individual cell membrane lipid bilayers apart, creating air compartments inside the bilayers. This reduces the free energy barrier for the translocation of drugs across the lipid bilayer and enhances drug permeability. At very high intensities, the two monolayers are largely separated, resulting in cell damage and implying that the blood brain barrier is primarily opened at the experimentally observed damaged areas.

Atomic Dislocations and Bond Rupture Govern Dissolution Enhancement under Acoustic Stimulation

By focusing the power of sound, acoustic stimulation (i.e., often referred to as sonication) enables numerous "green chemistry" pathways to enhance chemical reaction rates, for instance, of mineral dissolution in aqueous environments. However, a clear understanding of the atomistic mechanism(s) by which acoustic stimulation promotes mineral dissolution remains unclear. Herein, by combining nanoscale observations of dissolving surface topographies using vertical scanning interferometry, quantifications of mineral dissolution rates via analysis of solution compositions using inductively coupled plasma optical emission spectrometry, and classical molecular dynamics simulations, we reveal how acoustic stimulation induces dissolution enhancement. Across a wide range of minerals (Mohs hardness ranging from 3 to 7, surface energy ranging from 0.3 to 7.3 J/m², and stacking fault energy ranging from 0.8 to 10.0 J/m²), we show that acoustic fields enhance mineral dissolution rates (reactivity) by inducing atomic dislocations and/or atomic bond rupture. The relative contributions of these mechanisms depend on the mineral's underlying mechanical properties. Based on this new understanding, we create a unifying model that comprehensively describes how cavitation and acoustic stimulation processes affect mineral dissolution rates.

Damage to Polystyrene Polymer Film by Shock Wave Induced Bubble Collapse

Metallic surfaces that are in contact with solutions are commonly used in numerous applications where these surfaces can be damaged by shock wave induced bubble collapse. Use of polymer films that coat such surfaces to prevent them from damage requires a better understanding of how much harm collapsing bubbles produce in the films. In this study, we report the results from coarse-grained molecular dynamics simulations to study the damage to polystyrene (PS) films coating a hard surface. The damage was caused by a collapsing nanobubble located in the proximity of the film and interacting with an impinging shock wave. This collapse produces a high-speed water jet that impacts the PS film with a greater force than the shock front and creates cavities/pits in the PS film. We observed that polymer molecules located in the jet vicinity undergo conformational extension in the direction perpendicular to the jet motion, while chain molecules in the rest of the film undergo compression. We also observed that damage to the film is sensitive to the strength of the shock wave.

Echographic and physical characterization of albumin-stabilized nanobubbles

There has been increasing interest in using nanobubbles (NBs) for ultrasound mediated drug delivery as well as for ultrasound imaging. Albumin NBs are especially attractive for its potential of becoming a versatile platform for drug carriers and molecular targeted therapy agents. However, physical characterization of NBs is generally considered to be difficult due to various technical issues, such as concentration limitations, nanoparticle contamination, etc. In the present study, we measured the size distribution, concentration and weight density of albumin stabilized NBs by means of multiple nanoscale measurement modalities. Laser nanoparticle tracking analysis, multicolor flow cytometry, resonance mass evaluation showed consistent measurement results of the NBs with low mass weight density and diameter size ranging from 100 nm to 400 nm. Furthermore, the NB solution showed excellent images by high frequency ultrasound (30–50 MHz) in flow model acoustic phantoms. The NBs also induced acute cell disruption by low intensity ultrasound (0.8 W/cm²) irradiation. We successfully fabricated and characterized albumin stabilized NBs which could serve as an effective platform for future theranositic agents.

Interaction mechanism between the focused ultrasound and lipid membrane at the molecular level

Focused ultrasound (FUS) has a wide range of medical applications. Nowadays, the diagnostic and therapeutic ultrasound procedures are routinely used; effects of ultrasound on biological systems at the molecular level are, however, not fully understood. Experimental results on the interaction of the cell membrane, a simplest but important system component, with ultrasound are controversial. Molecular dynamics (MD) simulations could provide valuable insights, but there is no single study on the mechanism of the FUS induced structural changes in cell membranes. With this in mind, we develop a simple method to include FUS into a standard MD simulation. Adopting the 1,2-dioleoyl-sn-glycero-3-phosphocholine lipid membrane as a representative model described by the MARTINI coarse-grained force field, and using experimental values of the ultrasound frequency and intensity, we show that the heat and bubble cavitation are not the primary direct mechanisms that cause structural changes in the membrane. The spatial pressure gradients between the focused and free regions and between the parallel and perpendicular directions to the membrane are the origin of the mechanism. These gradients force lipids to move out of the focused region, forming a lipid flow along the membrane diagonal. Lipids in the free region move in the opposite direction due to the conservation of the total momentum. These opposite motions create wrinkles along the membrane diagonal at low FUS intensities and tear up the membrane at high FUS intensities. Once the membrane is torn up, it is not easy to reform. The implication of our findings in the FUS-induced drug delivery is discussed in some detail.

Damage and Failure of Axonal Microtubule under Extreme High Strain Rate: An In-Silico Molecular Dynamics Study

As a major cytoskeleton element of the axon, the breaking of microtubules (MTs) has been considered as a major cause of the axon degeneration. High strain rate loading is considered as one of the key factors in microtubule breaking. Due to the small size of microtubule, the real-time behavior of microtubule breaking is hard to capture. This study employs fully-atomistic molecular dynamics (MD) simulation to determine the failure modes of microtubule under different loadings conditions such as, unidirectional stretching, bending and hydrostatic expansion. For each loading conditions, MT is subjected to extreme high strain rate (108-109 s-1) loading. We argue that such level of high strain rate may be realized during cavitation bubble implosion. For each loading type, we have determined the critical energy for MT rupture. The associated rupture mechanisms are also discussed. We observed that the stretching has the lowest energy barrier to break the MT at the nanosecond time scale. Moreover, the breakage between the dimers starts at ~16% of total strain when stretched, which is much smaller compared to the reported strain-at-failure (50%) for lower strain rate loading. It suggests that MT fails at a significantly smaller strain states when loaded at higher strain rates.

Cavitation-threshold Determination and Rheological-parameters Estimation of Albumin-stabilized Nanobubbles

Nanobubbles (NBs) are of high interest for ultrasound (US) imaging as contrast agents and therapy as cavitation nuclei. Because of their instability (Laplace pressure bubble catastrophe) and low sensitivity to US, reducing the size of commonly used microbubbles to submicron-size is not trivial. We introduce stabilized NBs in the 100-250-nm size range, manufactured by agitating human serum albumin and perfluoro-propane. These NBs were exposed to 3.34- and 5.39-MHz US, and their sensitivity to US was proven by detecting inertial cavitation. The cavitation-threshold information was used to run a numerical parametric study based on a modified Rayleigh-Plesset equation (with a Newtonian rheology model). The determined values of surface tension ranged from 0 N/m to 0.06 N/m. The corresponding values of dilatational viscosity ranged from 5.10-10 Ns/m to 1.10-9 Ns/m. These parameters were reported to be 0.6 N/m and 1.10-8 Ns/m for the reference microbubble contrast agent. This result suggests the possibility of using albumin as a stabilizer for the nanobubbles that could be maintained in circulation and presenting satisfying US sensitivity, even in the 3-5-MHz range.

Shock wave-induced permeabilization of mammalian cells

Controlled permeabilization of mammalian cell membranes is fundamental to develop gene and cell therapies based on macromolecular cargo delivery, a process that emerged against an increasing number of health afflictions, including genetic disorders, cancer and infections. Viral vectors have been successfully used for macromolecular delivery; however, they may have unpredictable side effects and have been limited to life-threatening cases. Thus, several chemical and physical methods have been explored to introduce drugs, vaccines, and nucleic acids into cells. One of the most appealing physical methods to deliver genes into cells is shock wave-induced poration. High-speed microjets of fluid, emitted due to the collapse of microbubbles after shock wave passage, represent the most significant mechanism that contributes to cell membrane poration by this technique. Herein, progress in shock wave-induced permeabilization of mammalian cells is presented. After covering the main concepts related to molecular strategies whose applications depend on safer drug delivery methods, the physics behind shock wave phenomena is described. Insights into the use of shock waves for cell membrane permeation are discussed, along with an overview of the two major biomedical applications thereof-i.e., genetic modification and anti-cancer shock wave-assisted chemotherapy. The aim of this review is to summarize 30 years of data showing underwater shock waves as a safe, noninvasive method for macromolecular delivery into mammalian cells, encouraging the development of further research, which is still required before the introduction of this promising tool into clinical practice.

Targeting Early Dementia: Using Lipid Cubic Phase Nanocarriers to Cross the Blood⁻Brain Barrier

Over the past decades, a frequent co-morbidity of cerebrovascular pathology and Alzheimer's disease has been observed. Numerous published studies indicate that the preservation of a healthy cerebrovascular endothelium can be an important therapeutic target. By incorporating the appropriate drug(s) into biomimetic (lipid cubic phase) nanocarriers, one obtains a multitasking combination therapeutic, which targets certain cell surface scavenger receptors, mainly class B type I (i.e., SR-BI), and crosses the blood⁻brain barrier. This targeting allows for various cell types related to Alzheimer's to be simultaneously searched out for localized drug treatment in vivo.

Nanotherapy for Alzheimer's disease and vascular dementia: Targeting senile endothelium

Due to the complexity of Alzheimer's disease, multiple cellular types need to be targeted simultaneously in order for a given therapy to demonstrate any major effectiveness. Ultrasound-sensitive coated microbubbles (in a targeted lipid nanoemulsion) are available. Versatile small molecule drug(s) targeting multiple pathways of Alzheimer's disease pathogenesis are known. By incorporating such drug(s) into the targeted "lipid-coated microbubble" [LCM]/"nanoparticle-derived" [ND] (or LCM/ND) nanoemulsion type, one obtains a multitasking combination therapeutic for translational medicine. This multitasking therapeutic targets cell-surface scavenger receptors (mainly class B type I), or SR-BI, making possible for various Alzheimer's-related cell types to be simultaneously searched out for localized drug treatment in vivo. Besides targeting cell-surface SR-BI, the proposed LCM/ND-nanoemulsion combination therapeutic(s) include a characteristic lipid-coated microbubble [LCM] subpopulation (i.e., a stable LCM suspension); such film-stabilized microbubbles are well known to substantially reduce the acoustic power levels needed for accomplishing temporary noninvasive (transcranial) ultrasound treatment, or sonoporation, if additionally desired for the Alzheimer's patient.

Alzheimer’s Disease, Brain Injury, and C.N.S. Nanotherapy in Humans: Sonoporation Augmenting Drug Targeting

Owing to the complexity of neurodegenerative diseases, multiple cellular types need to be targeted simultaneously in order for a given therapy to demonstrate any major effectiveness. Ultrasound-sensitive coated microbubbles (in a targeted nanoemulsion) are available. Versatile small-molecule drug(s) targeting multiple pathways of Alzheimer’s disease pathogenesis are known. By incorporating such drug(s) into the targeted lipid-coated microbubble/nanoparticle-derived (LCM/ND) lipid nanoemulsion type, one obtains a multitasking combination therapeutic for translational medicine. This multitasking therapeutic targets cell-surface scavenger receptors (mainly scavenger receptor class B type I (SR-BI)), making it possible for various Alzheimer’s-related cell types to be simultaneously sought for localized drug treatment in vivo. Besides targeting cell-surface SR-BI, the proposed LCM/ND-nanoemulsion combination therapeutic(s) include a characteristic lipid-coated microbubble (LCM) subpopulation (i.e., a stable LCM suspension); such LCM substantially reduce the acoustic power levels needed for accomplishing temporary noninvasive (transcranial) ultrasound treatment, or sonoporation, if additionally desired for the Alzheimer’s patient.

Effect of Shock-Induced Cavitation Bubble Collapse on the damage in the Simulated Perineuronal Net of the Brain

The purpose of this study is to conduct modeling and simulation to understand the effect of shock-induced mechanical loading, in the form of cavitation bubble collapse, on damage to the brain's perineuronal nets (PNNs). It is known that high-energy implosion due to cavitation collapse is responsible for corrosion or surface damage in many mechanical devices. In this case, cavitation refers to the bubble created by pressure drop. The presence of a similar damage mechanism in biophysical systems has long being suspected but not well-explored. In this paper, we use reactive molecular dynamics (MD) to simulate the scenario of a shock wave induced cavitation collapse within the perineuronal net (PNN), which is the near-neuron domain of a brain's extracellular matrix (ECM). Our model is focused on the damage in hyaluronan (HA), which is the main structural component of PNN. We have investigated the roles of cavitation bubble location, shockwave intensity and the size of a cavitation bubble on the structural evolution of PNN. Simulation results show that the localized supersonic water hammer created by an asymmetrical bubble collapse may break the hyaluronan. As such, the current study advances current knowledge and understanding of the connection between PNN damage and neurodegenerative disorders.

Molecular mechanism for cavitation in water under tension

Despite its relevance in biology and engineering, the molecular mechanism driving cavitation in water remains unknown. Using computer simulations, we investigate the structure and dynamics of vapor bubbles emerging from metastable water at negative pressures. We find that in the early stages of cavitation, bubbles are irregularly shaped and become more spherical as they grow. Nevertheless, the free energy of bubble formation can be perfectly reproduced in the framework of classical nucleation theory (CNT) if the curvature dependence of the surface tension is taken into account. Comparison of the observed bubble dynamics to the predictions of the macroscopic Rayleigh-Plesset (RP) equation, augmented with thermal fluctuations, demonstrates that the growth of nanoscale bubbles is governed by viscous forces. Combining the dynamical prefactor determined from the RP equation with CNT based on the Kramers formalism yields an analytical expression for the cavitation rate that reproduces the simulation results very well over a wide range of pressures. Furthermore, our theoretical predictions are in excellent agreement with cavitation rates obtained from inclusion experiments. This suggests that homogeneous nucleation is observed in inclusions, whereas only heterogeneous nucleation on impurities or defects occurs in other experiments.

Biophysical insight into mechanisms of sonoporation

This study presents a unique approach to understanding the biophysical mechanisms of ultrasound-triggered cell membrane disruption (i.e., sonoporation). We report direct correlations between ultrasound-stimulated encapsulated microbubble oscillation physics and the resulting cellular membrane permeability by simultaneous microscopy of these two processes over their intrinsic physical timescales (microseconds for microbubble dynamics and seconds to minutes for local macromolecule uptake and cell membrane reorganization). We show that there exists a microbubble oscillation-induced shear-stress threshold, on the order of kilopascals, beyond which endothelial cellular membrane permeability increases. The shear-stress threshold exhibits an inverse square-root relation to the number of oscillation cycles and an approximately linear dependence on ultrasound frequency from 0.5 to 2 MHz. Further, via real-time 3D confocal microscopy measurements, our data provide evidence that a sonoporation event directly results in the immediate generation of membrane pores through both apical and basal cell membrane layers that reseal along their lateral area (resealing time of ∼<2 min). Finally, we demonstrate the potential for sonoporation to indirectly initiate prolonged, intercellular gaps between adjacent, confluent cells (∼>30-60 min). This real-time microscopic approach has provided insight into both the physical, cavitation-based mechanisms of sonoporation and the biophysical, cell-membrane-based mechanisms by which microbubble acoustic behaviors cause acute and sustained enhancement of cellular and vascular permeability.

Nanobubbles, cavitation, shock waves and traumatic brain injury

Shock wave induced cavitation denaturates blood–brain barrier tight junction proteins; this may result in various neurological complications.

Opening of the blood-brain barrier tight junction due to shock wave induced bubble collapse: a molecular dynamics simulation study

Passage of a shock wave across living organisms may produce bubbles in the blood vessels and capillaries. It was suggested that collapse of these bubbles imposed by an impinging shock wave can be responsible for the damage or even destruction of the blood-brain barrier. To check this possibility, we performed molecular dynamics computer simulations on systems that contained a model of tight junction from the blood-brain barrier. In our model, we represent the tight junction by two pairs of interacting proteins, claudin-15. Some of the simulations were done in the absence of a nanobubble, some in its presence. Our simulations show that when no bubble is present in the system, no damage to tight junction is observed when the shock wave propagates across it. In the presence of a nanobubble, even when the impulse of the shock wave is relatively low, the implosion of the bubble causes serious damage to our model tight junction.