Focused Ultrasound (FUS)-triggered nano-sized drug delivery, as a smart stimuli-responsive system for treating solid tumors, is computationally investigated to enhance localized delivery of drug and treatment efficacy. Integration of thermosensitive liposome (TSL), as a doxorubicin (DOX)-loaded nanocarrier, and FUS, provides a promising drug delivery system. A fully coupled partial differential system of equations, including the Helmholtz equation for FUS propagation, bio-heat transfer, interstitial fluid flow, drug transport in tissue and cellular spaces, and a pharmacodynamic model is first presented for this treatment approach. Equations are then solved by finite element methods to calculate intracellular drug concentration and treatment efficacy. The main objective of this study is to present a multi-physics and multi-scale model to simulate drug release, transport, and delivery to solid tumors, followed by an analysis of how FUS exposure time and drug release rate affect these processes. Our findings not only show the capability of model to replicate this therapeutic approach, but also confirm the benefits of this treatment with an improvement of drug aggregation in tumor and reduction of drug delivery in healthy tissue. For instance, the survival fraction of tumor cells after this treatment dropped to 62.4%, because of a large amount of delivered drugs to cancer cells. Next, a combination of three release rates (ultrafast, fast, and slow) and FUS exposure times (10, 30, and 60 min) was examined. Area under curve (AUC) results show that the combination of 30 min FUS exposure and rapid drug release leads to a practical and effective therapeutic response.

Although transcranial photoacoustic imaging has been previously investigated by several groups, there are many unknowns about the distorting effects of the skull due to the impedance mismatch between the skull and underlying layers. The current computational methods based on finite-element modeling are slow, especially in the cases where fine grids are defined for a large 3-D volume. We develop a very fast modeling/simulation framework based on deterministic ray-tracing. The framework considers a multilayer model of the medium, taking into account the frequency-dependent attenuation and dispersion effects that occur in wave reflection, refraction, and mode conversion at the skull surface. The speed of the proposed framework is evaluated. We validate the accuracy of the framework using numerical phantoms and compare its results to k-Wave simulation results. Analytical validation is also performed based on the longitudinal and shear wave transmission coefficients. We then simulated, using our method, the major skull-distorting effects including amplitude attenuation, time-domain signal broadening, and time shift, and confirmed the findings by comparing them to several ex vivo experimental results. It is expected that the proposed method speeds up modeling and quantification of skull tissue and allows the development of transcranial photoacoustic brain imaging.

Computational modeling enables researchers to study and understand various complex biological phenomena in anticancer drug delivery systems (DDSs), especially nano‐sized DDSs (NSDDSs). The combination of NSDDSs and therapeutic ultrasound (TUS), that is, focused ultrasound and low‐intensity pulsed ultrasound, has made significant progress in recent years, opening many opportunities for cancer treatment. Multiple parameters require tuning and optimization to develop effective DDSs, such as NSDDSs, in which mathematical modeling can prove advantageous. In silico computational modeling of ultrasound‐responsive DDS typically involves a complex framework of acoustic interactions, heat transfer, drug release from nanoparticles, fluid flow, mass transport, and pharmacodynamic governing equations. Owing to the rapid development of computational tools, modeling the different phenomena in multi‐scale complex problems involved in drug delivery to tumors has become possible. In the present study, we present an in‐depth review of recent advances in the mathematical modeling of TUS‐mediated DDSs for cancer treatment. A detailed discussion is also provided on applying these computational models to improve the clinical translation for applications in cancer treatment.This article is categorized under: Nanotechnology Approaches to Biology > Nanoscale Systems in Biology Therapeutic Approaches and Drug Discovery > Nanomedicine for Oncologic Disease

Currently, no numerical model for low-intensity pulsed ultrasound (LIPUS)-triggered anticancer drug release from gold nanoparticle (GNP) drug carriers exists in the literature. In this work, LIPUS-induced doxorubicin (DOX) release from GNPs was achieved in an ex vivo tissue model. Transmission electronic microscopy (TEM) imaging was performed before and after LIPUS exposure, and significant aggregation of the GNPs was observed upon DOX release. Subsequently, GNP surface potential was determined before and after LIPUS-induced DOX release, using a Zetasizer. A numerical model was then created to predict GNP aggregation, and the subsequent DOX release, via combining a thermal field simulation by solving the bioheat transfer equation (in COMSOL) and the Derjaguin, Landau, Verwey, and Overbeek (DLVO) total interaction potential (in MATLAB). The DLVO model was applied to the colloidal DOX-loaded GNPs by summing the attractive van der Waals and electrostatic repulsion interaction potentials for any given GNP pair. DLVO total interaction potential was found before and after LIPUS exposure, and an energy barrier for aggregation was determined. The DLVO interaction potential peak amplitude was found to drop from 1.36 k_BT to 0.24 k_BT after LIPUS exposure, translating to an 82.4% decrease in peak amplitude value. It was concluded that the interaction potential energy threshold for GNP aggregation (and, as a result, DOX release) was equal to 0.24 k_BT.

This study investigated the dynamic changes of tissue attenuation coefficients before, during, and after high-intensity focused ultrasound (HIFU) treatment at different total acoustic powers (TAP) in ex vivo porcine muscle tissue. It further assessed the reliability of employing changes in tissue attenuation coefficient parameters as potential indicators of tissue thermal damage.