Drug delivery through nanoparticles in solid tumors: a mechanistic understanding

Agent-based modeling in cancer biomedicine: applications and tools for calibration and validation

Computational models are not just appealing because they can simulate and predict the development of biological phenomena across multiple spatial and temporal scales, but also because they can integrate information from well-established in vitro and in vivo models and test new hypotheses in cancer biomedicine. Agent-based models and simulations are especially interesting candidates among computational modeling procedures in cancer research due to the capability to, for instance, recapitulate the dynamics of neoplasia and tumor - host interactions. Yet, the absence of methods to validate the consistency of the results across scales can hinder adoption by turning fine-tuned models into black boxes. This review compiles relevant literature that explores strategies to leverage high-fidelity simulations of multi-scale, or multi-level, cancer models with a focus on verification approached as simulation calibration. We consolidate our review with an outline of modern approaches for agent-based models' validation and provide an ambitious outlook toward rigorous and reliable calibration.

Radiopharmaceutical transport in solid tumors via a 3-dimensional image-based spatiotemporal model

Lutetium-177 prostate-specific membrane antigen (177Lu-PSMA)-targeted radiopharmaceutical therapy is a clinically approved treatment for patients with metastatic castration-resistant prostate cancer (mCRPC). Even though common practice reluctantly follows "one size fits all" approach, medical community believes there is significant room for deeper understanding and personalization of radiopharmaceutical therapies. To pursue this aim, we present a 3-dimensional spatiotemporal radiopharmaceutical delivery model based on clinical imaging data to simulate pharmacokinetic of 177Lu-PSMA within the prostate tumors. The model includes interstitial flow, radiopharmaceutical transport in tissues, receptor cycles, association/dissociation with ligands, synthesis of PSMA receptors, receptor recycling, internalization of radiopharmaceuticals, and degradation of receptors and drugs. The model was studied for a range of values for injection amount (100-1000 nmol), receptor density (10-500 nmol•l-1), and recycling rate of receptors (10-4 to 10-1 min-1). Furthermore, injection type, different convection-diffusion-reaction mechanisms, characteristic time scales, and length scales are discussed. The study found that increasing receptor density, ligand amount, and labeled ligands improved radiopharmaceutical uptake in the tumor. A high receptor recycling rate (0.1 min-1) increased radiopharmaceutical concentration by promoting repeated binding to tumor cell receptors. Continuous infusion results in higher radiopharmaceutical concentrations within tumors compared to bolus administration. These insights are crucial for advancing targeted therapy for prostate cancer by understanding the mechanism of radiopharmaceutical distribution in tumors. Furthermore, measures of characteristic length and advection time scale were computed. The presented spatiotemporal tumor transport model can analyze different physiological parameters affecting 177Lu-PSMA delivery.

Spatiotemporal modeling of nano-delivered chemotherapeutics for synergistic microwave ablation cancer therapy

BACKGROUND AND OBJECTIVE

The effectiveness of current microwave ablation (MWA) therapies is limited. Administration of thermosensitive liposomes (TSLs) which release drugs in response to heat has presented a significant potential for enhancing the efficacy of thermal ablation treatment, and the benefits of targeted drug delivery. However, a complete knowledge of the mechanobiological processes underlying the drug release process, especially the intravascular drug release mechanism and its distribution in response to MWA needs to be improved. Multiscale computational-based modeling frameworks, integrating different biophysical phenomena, have recently emerged as promising tools to decipher the mechanobiological events in combo therapies. The present study aims to develop a novel multiscale computational model of TSLs delivery following MWA implantation.

METHODS

Due to the complex interplay between the heating procedure and the drug concentration maps, a computational model is developed to determine the intravascular release of doxorubicin from TSL, its transvascular transport into the interstitium, transport in the interstitium, and cell uptake. Computational models can estimate the interplays among liposome and drug properties, tumor perfusion, and heating regimen to examine the impact of essential parameters and to optimize a targeted drug delivery platform.

RESULTS

Results indicated that the synergy of TSLs with MWA allows more localized drug delivery with lower side effects. The drug release rate and tumor permeability play crucial roles in the efficacy of TSLs during MWA treatment. The computational model predicted an unencapsulated drug lime around the ablated zone, which can destroy more cancer cells compared to MWA alone by 40%. Administration of TSLs with a high release rate capacity can improve the percentage of killed cancer cells by 24%. Since the heating duration in MWA is less than 15 min, the presented combination therapy showed better performance for highly permeable tumors.

CONCLUSION

This study highlights the potential of the proposed computational framework to address complex and realistic scenarios in cancer treatment, which can serve as the future research foundation, including advancements in nanomedicine and optimizing the pair of TSL and MWA for both preclinical and clinical studies. The present model could be as a valuable tool for patient-specific calibration of essential parameters.

Effect of tumor heterogeneity on enhancing drug delivery to vascularized tumors using thermo-sensitive liposomes triggered by hyperthermia: A multi-scale and multi-physics computational model

In this study, a novel multi-scale and multi-physics image-based computational model is introduced to assess the delivery of doxorubicin (Dox) loaded temperature-sensitive liposomes (TSLs) in the presence of hyperthermia. Unlike previous methodologies, this approach incorporates capillary network geometry extracted from images, resulting in a more realistic physiological tumor model. This model holds significant promise in advancing personalized medicine by integrating patient-specific tumor properties. The finite element method is employed to solve the equations governing intravascular and interstitial fluid flows, as well as the transport of therapeutic agents within the tissue. Realistic biological conditions and intricate processes like intravascular pressure, drug binding to cells, and cellular uptake are also considered to enhance the model's accuracy. The results underscore the significant impact of vascular architecture on treatment outcomes. Variation in vascular network pattern yielded changes of up to 38 % in the fraction of killed cells (FKCs) parameter under identical conditions. Pressure control of the parent vessels can also improve FKCs by approximately 17 %. Tailoring the treatment plan based on tumor-specific parameters emerged as a critical factor influencing treatment efficacy. For instance, changing the time interval between the administration of Dox-loaded TSLs and hyperthermia can result in a 48 % improvement in treatment outcomes. Additionally, devising a customized heating schedule led to a 20 % increase in treatment efficacy. Our proposed model highlights the significant effect of tumor characteristics and vascular network structure on the final treatment outcomes of the presented combination treatment.

Spatiotemporal modeling of radiopharmaceutical transport in solid tumors: Application to 177Lu-PSMA therapy of prostate cancer

BACKGROUND AND OBJECTIVE

177Lu-labeled prostate-specific membrane antigen (PSMA) radiopharmaceutical therapy (RPT) represents a pivotal advancement in addressing prostate cancer. However, existing therapies, while promising, remain incompletely understood and optimized. Computational models offer potential insights into RPTs, aiding in clinical drug delivery enhancement. In this study, we investigate the impact of various physiological parameters on the delivery of 177Lu-PSMA-617 RPT using the convection-diffusion-reaction (CDR) model.

METHODS

Our investigation encompasses tumor geometry and surrounding tissue, characterized by well-defined boundaries and initial conditions. Utilizing the finite element method, we solve governing equations across a range of parameters: dissociation constant KD (1, 0.1, 0.01 [nM]), internalization rate (0.01-0.0001 [min-1]), diverse tumor shapes, and variable necrotic zone sizes. This model can provide an accurate analysis of radiopharmaceutical delivery from the injection site to the tumor cell, including drug transport in the vascular, interstitial, and intracellular spaces, and considering important parameters (e.g., drug extravasation from microvessels or to lymphatic vessels, the extracellular matrix, receptors, and intracellular space).

RESULTS

Our findings reveal significant enhancements in tumor-absorbed doses as KD decreases. This outcome can be attributed to the higher affinity of radiopharmaceuticals for PSMA receptors as KD diminishes, facilitating a more efficient binding and retention of the therapeutic agent within the tumor microenvironment. Additionally, tumor-absorbed doses for KD ∼ 1 [nM] show an upward trend with higher internalization rates. This observation can be rationalized by considering that a greater internalization rate would result in a higher proportion of radiopharmaceuticals being taken up by tumor cells after binding to receptors on the cell surface. Notably, tumor shape and necrotic zone size exhibit limited influence on tumor absorbed dose.

CONCLUSIONS

The present study employs the CDR model to explore the role of physiological parameters in shaping 177Lu-PSMA-617 RPT delivery. These findings provide insights for improving prostate cancer therapy by understanding radiopharmaceutical transport dynamics. This computational approach contributes to advancing our understanding of radiopharmaceutical delivery mechanisms and has implications for enhancing treatment efficacy.

Ultrasound-mediated nano-sized drug delivery systems for cancer treatment: Multi-scale and multi-physics computational modeling

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.

Multi-Scale and Multi-Physics Models of the Transport of Therapeutic/Diagnostic Cancer Agents

The effectiveness of tumor treatment heavily relies on the successful delivery of anticancer drugs [...].

Stimuli-sensitive nano-drug delivery with programmable size changes to enhance accumulation of therapeutic agents in tumors

Nano-based drug delivery systems hold significant promise for cancer therapies. Presently, the poor accumulation of drug-carrying nanoparticles in tumors has limited their success. In this study, based on a combination of the paradigms of intravascular and extravascular drug release, an efficient nanosized drug delivery system with programmable size changes is introduced. Drug-loaded smaller nanoparticles (secondary nanoparticles), which are loaded inside larger nanoparticles (primary nanoparticles), are released within the microvascular network due to temperature field resulting from focused ultrasound. This leads to the scale of the drug delivery system decreasing by 7.5 to 150 times. Subsequently, smaller nanoparticles enter the tissue at high transvascular rates and achieve higher accumulation, leading to higher penetration depths. In response to the acidic pH of tumor microenvironment (according to the distribution of oxygen), they begin to release the drug doxorubicin at very slow rates (i.e., sustained release). To predict the performance and distribution of therapeutic agents, a semi-realistic microvascular network is first generated based on a sprouting angiogenesis model and the transport of therapeutic agents is then investigated based on a developed multi-compartment model. The results show that reducing the size of the primary and secondary nanoparticles can lead to higher cell death rate. In addition, tumor growth can be inhibited for a longer time by enhancing the bioavailability of the drug in the extracellular space. The proposed drug delivery system can be very promising in clinical applications. Furthermore, the proposed mathematical model is applicable to broader applications to predict the performance of drug delivery systems.

A comparative study between conventional chemotherapy and photothermal activated nano-sized targeted drug delivery to solid tumor

Delivery of chemotherapeutic medicines to solid tumors is critical for optimal therapeutic success and minimal adverse effects. We mathematically developed a delivery method using thermosensitive nanocarriers activated by light irradiation. To assess its efficacy and identify critical events and parameters affecting therapeutic response, we compared this method to bolus and continuous infusions of doxorubicin for both single and multiple administrations. A hybrid sprouting angiogenesis approach generates a semi-realistic microvascular network to evaluate therapeutic drug distribution and microvascular heterogeneity. A pharmacodynamics model evaluates treatment success based on tumor survival cell percentage. The study found that whereas bolus injection boosted extracellular drug concentration levels by 90%, continuous infusion improved therapeutic response due to improved bioavailability. Cancer cell death increases by 6% with several injections compared to single injections due to prolonged chemotherapeutic medication exposure. However, responsive nanocarriers supply more than 2.1 times more drug than traditional chemotherapy in extracellular space, suppressing tumor development longer. Also, controlled drug release decreases systemic side effects substantial through diminishing the concentration of free drug in the circulation. The primary finding of this work highlights the significance of high bioavailability in treatment response. The results indicate that responsive nanocarriers contribute to increased bioavailability, leading to improved therapeutic benefits. By including drug delivery features in a semi-realistic model, this numerical study sought to improve drug-bio interaction comprehension. The model provides a good framework for understanding preclinical and clinical targeted oncology study outcomes.

Ultrasound-mediated nano drug delivery for treating cancer: Fundamental physics to future directions

The application of biocompatible nanocarriers in medicine has provided several benefits over conventional treatment methods. However, achieving high treatment efficacy and deep penetration of nanocarriers in tumor tissue is still challenging. To address this, stimuli-responsive nano-sized drug delivery systems (DDSs) are an active area of investigation in delivering anticancer drugs. While ultrasound is mainly used for diagnostic purposes, it can also be applied to affect cellular function and the delivery/release of anticancer drugs. Therapeutic ultrasound (TUS) has shown potential as both a stand-alone anticancer treatment and a method to induce targeted drug release from nanocarrier systems. TUS approaches have been used to overcome various physiological obstacles, including endothelial barriers, the tumor microenvironment (TME), and immunological hurdles. Combining nanomedicine and ultrasound as a smart DDS can increase in situ drug delivery and improve access to impermeable tissues. Furthermore, smart DDSs can perform targeted drug release in response to distinctive TMEs, external triggers, or dual/multi-stimulus. This results in enhanced treatment efficacy and reduced damage to surrounding healthy tissue or organs at risk. Integrating DDSs and ultrasound is still in its early stages. More research and clinical trials are required to fully understand ultrasound's underlying physical mechanisms and interactions with various types of nanocarriers and different types of cells and tissues. In the present review, ultrasound-mediated nano-sized DDS, specifically focused on cancer treatment, is presented and discussed. Ultrasound interaction with nanoparticles (NPs), drug release mechanisms, and various types of ultrasound-sensitive NPs are examined. Additionally, in vitro, in vivo, and clinical applications of TUS are reviewed in light of the critical challenges that need to be considered to advance TUS toward an efficient, secure, straightforward, and accessible cancer treatment. This study also presents effective TUS parameters and safety considerations for this treatment modality and gives recommendations about system design and operation. Finally, future perspectives are considered, and different TUS approaches are examined and discussed in detail. This review investigates drug release and delivery through ultrasound-mediated nano-sized cancer treatment, both pre-clinically and clinically.

Mathematical analysis to the variation of tumor cells density according to recovery of five organs: Kidney, liver, heart, spleen, and lung

BACKGROUND AND OBJECTIVE

Harmonious interactions of five representative organs: kidney, liver, heart, spleen, and lung, improve metastasis and cell divisions, and abnormal cell division causes cancer cell development. The research is processed through a mathematical approach based on win-win principle of five organs to generate medicine in blood vessel. The variations of solute medicine amount in blood vessel with respect to the flow rates of injected drugs are interpreted. The alterations of tumor cells density and tumor angiogenesis factor concentration are described according to the recovery of five organs' functions.

METHODS

A compartmental analysis is applied to obtain medicine concentration in blood vessel by the functional recovery of five organs considering time level t_i, the reaction rate coefficient R_j, and the medicine flow rate α. Random motility and chemotaxis in response to tumor angiogenesis factor gradients are comprised to derive mathematical governing equations for tumor cells motion and a finite volume method with time-changing is adopted to obtain numerical solutions due to the complexity of the governing equations.

RESULTS

Drug concentration in blood vessel grows as heart reaction rate increases, and the medicine made through the functional enhancements of five representative organs is highly influential to restrain the activity of tumor angiogenesis factor. With the growth of medicine concentration in blood vessel according to the decline of reaction rate and medicine flow rate, tumor cells reacts hypersensitively at the moment of medicines injection and the density of tumor cells approached to zero.

CONCLUSIONS

Consequently, reaction rate, time level, and medicine flow rate are crucial factors in the determination of medicine amount in blood vessel and to control tumor angiogenesis factor concentration, and harmonious balanced functions among five organs based on win-win principle contribute to control the activity of tumor cells.

Nanomedicine: New Frontiers in Fighting Microbial Infections

Microbes have dominated life on Earth for the past two billion years, despite facing a variety of obstacles. In the 20th century, antibiotics and immunizations brought about these changes. Since then, microorganisms have acquired resistance, and various infectious diseases have been able to avoid being treated with traditionally developed vaccines. Antibiotic resistance and pathogenicity have surpassed antibiotic discovery in terms of importance over the course of the past few decades. These shifts have resulted in tremendous economic and health repercussions across the board for all socioeconomic levels; thus, we require ground-breaking innovations to effectively manage microbial infections and to provide long-term solutions. The pharmaceutical and biotechnology sectors have been radically altered as a result of nanomedicine, and this trend is now spreading to the antibacterial research community. Here, we examine the role that nanomedicine plays in the prevention of microbial infections, including topics such as diagnosis, antimicrobial therapy, pharmaceutical administration, and immunizations, as well as the opportunities and challenges that lie ahead.

Multiphysics Modeling of Low-Intensity Pulsed Ultrasound Induced Chemotherapeutic Drug Release from the Surface of Gold Nanoparticles

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.

HYDRHA: Hydrogels of hyaluronic acid. New biomedical approaches in cancer, neurodegenerative diseases, and tissue engineering

In the last decade, hyaluronic acid (HA) has attracted an ever-growing interest in the biomedical engineering field as a biocompatible, biodegradable, and chemically versatile molecule. In fact, HA is a major component of the extracellular matrix (ECM) and is essential for the maintenance of cellular homeostasis and crosstalk. Innovative experimental strategies in vitro and in vivo using three-dimensional (3D) HA systems have been increasingly reported in studies of diseases, replacement of tissue and organ damage, repairing wounds, and encapsulating stem cells for tissue regeneration. The present work aims to give an overview and comparison of recent work carried out on HA systems showing advantages, limitations, and their complementarity, for a comprehensive characterization of their use. A special attention is paid to the use of HA in three important areas: cancer, diseases of the central nervous system (CNS), and tissue regeneration, discussing the most innovative experimental strategies. Finally, perspectives within and beyond these research fields are discussed.

Application of a multicomponent model of convectional reaction-diffusion to description of glucose gradients in a neurovascular unit

A supply of glucose to a nervous tissue is fulfilled by a cerebrovascular network, and further diffusion is known to occur at both an arteriolar and a microvascular level. Despite a direct relation, a blood flow dynamic and reaction-diffusion of metabolites are usually considered separately in the mathematical models. In the present study they are coupled in a multiphysical approach which allows to evaluate the effects of capillary blood flow changes on near-vessels nutrient concentration gradients evidently. Cerebral blood flow (CBF) was described by the non-steady-state Navier-Stokes equations for a non-Newtonian fluid whose constitutive law is given by the Carreau model. A three-level organization of blood–brain barrier (BBB) is modelled by the flux dysconnectivity functions including densities and kinetic properties of glucose transporters. The velocity of a fluid flow in brain extracellular space (ECS) was estimated using Darcy’s law. The equations of reaction-diffusion with convection based on a generated flow field for continues and porous media were used to describe spatial-time gradients of glucose in the capillary lumen and brain parenchyma of a neurovascular unit (NVU), respectively. Changes in CBF were directly simulated using smoothing step-like functions altering the difference of intracapillary pressure in time. The changes of CBF cover both the decrease (on 70%) and the increase (on 50%) in a capillary flow velocity. Analyzing the dynamics of glucose gradients, it was shown that a rapid decrease of a capillary blood flow yields an enhanced level of glucose in a near-capillary nervous tissue if the contacts between astrocytes end-feet are not tight. Under the increased CBF velocities the amplitude of glucose concentration gradients is always enhanced. The introduced approach can be used for estimation of blood flow changes influence not only on glucose but also on other nutrients concentration gradients and for the modelling of distributions of their concentrations near blood vessels in other tissues as well.