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 [...].


Introduction
The effectiveness of tumor treatment heavily relies on the successful delivery of anticancer drugs.Solid tumors present physiological barriers, such as a dense extracellular matrix (ECM), enhanced solid stress and vessel compression, and irregularly shaped and leaky microvascular networks, reducing drug delivery efficacy [1][2][3].These barriers and an absent or impaired lymphatic function contribute to elevated interstitial fluid pressure (IFP), hindering systematic drug delivery due to outward convection counteracting inward drug diffusion [1][2][3][4][5].Each tumor possesses unique structural and physiochemical characteristics, such as the surface area of the vessel wall per unit volume for transport into the tissue, capillary network patterns, and hydraulic conductivity, making drug delivery and treatment highly complex and case-specific.The quantitative prediction of drug delivery to solid tumors using patient-specific data is crucial for clinical decision-making in therapeutic planning.
Mathematical modeling and simulation aim to enhance our understanding of tumor behaviors, ultimately improving treatment outcomes.Two major approaches for simulating drug delivery to solid tumors exist: the macroscopic and microscopic approaches.The macroscopic approach emphasizes the distribution of variables such as interstitial fluid velocity (IFV), IFP, and drug concentration over the tumor radius length scale.In the microscopic approach, characteristics like the microvascular network structure, blood flow within microvessels, and the interaction between blood flow and the microvascular wall are incorporated into the model.Given several transport processes involved in the delivery of therapeutic/diagnostic cancer agents and the complexity of the tumor microenvironment, sophisticated mathematical/computational modeling can be used to study the limitations of these approaches in treating/diagnosing cancer.Recently, multi-physics and multi-scale models have been applied to aid in the development of therapeutic/diagnostic agent delivery approaches.A wide range of mathematical models have been developed to mimic in vitro/in vivo/human biological and physiological environments and simulate these agents' behaviors [6,7].Various drug delivery systems (e.g., conventional chemotherapy [8][9][10][11][12][13][14], implantable systems [15,16], intraperitoneal chemotherapy [17][18][19][20], intratumoral injection [21], and anti-angiogenesis/vascular normalization combined with another therapy [22][23][24]) and diagnostic agents (e.g., the FDG PET radiotracer [25], the FMISO PET radiotracer [26], drug-loaded microbubbles [27], and radiopharmaceuticals [28]) have been explored using mathematical approaches.
On the other hand, medical imaging and image processing methods can be used to estimate, in principle, patient-specific variables that can be used to increase the accuracy of computational simulations.Consequently, there is considerable interest in predicting therapeutic/diagnostic agent delivery techniques through image-based computational modeling, facilitating personalized predictions.Different outputs of imaging modalities, particularly magnetic resonance imaging (MRI) [29][30][31], computed tomography (CT) [32], and diffusion tensor imaging (DTI) [33,34], can be used as inputs for these computational personalized models, enabling the timely and cost-effective optimization of drug delivery plans.
Nanoparticles have shown promise in preclinical studies for drug delivery in cancer treatment.Their drug-release kinetics and drug-loading capacity are crucial properties that determine treatment efficacy.In addition, stimulus-responsive nanoparticles, using external (magnetic, electric, light, laser, ultrasound, etc.) and/or internal (pH, redox, hypoxia, etc.) stimuli, significantly reduce the side effects of chemotherapeutics by controlling their release kinetics.Upon reaching the target site and responding to an internal or external stimulus, nanoparticles initiate drug release, either gradually (sustainable release) or rapidly (burst release), or as a combination of both [35,36].However, nanoparticles encounter various barriers, such as rapid clearance and low transvascular rates.Some nanoparticles rely on the enhanced permeability and retention (EPR) effect for passive targeting, leading to their accumulation in tissues.Taking into account parameters associated with therapeutic agents (drugs and nanoparticles) and tumor-related factors, as well as the multiple scales involved (tissue, cell, and nano scales) and the physical nature of the transport process, computational modeling can be used to design efficient nano-sized drug delivery systems for solid tumors [6,[37][38][39][40][41][42][43].
Computational models combine concepts and ideas from various disciplines, including pharmacokinetics, pharmacodynamics, fluid mechanics/dynamics, tissue mechanics, mass transport, heat transfer and biochemical processes.Multi-physics modeling allows for integrating information from various disciplines to inform the development of therapeutic/diagnostic agents [44].Delivery processes can be evaluated separately or by combining all stages, aiding in identifying opportunities to maximize delivery outcomes and treatment effectiveness.This approach can reveal limitations in the transport and delivery mechanisms of the agents and can help guide therapeutic/diagnostic agents' development.On the other hand, to improve the prediction capability of drug delivery, recent efforts in multiscale modeling to develop advanced models that can simulate transport processes across multiple length and time scales have emerged, integrating discrete and continuum modeling strategies [6,45].These approaches promise to establish bottom-up in silico pipelines for drug design and delivery.Nevertheless, challenges persist in model parametrization and validation, particularly in the presence of variability introduced by various levels of heterogeneity in diseased tissues.When coupled with patient-specific data, these multi-physics and multi-scale models possess considerable potential for validation and, in some cases, personalization (such as creating digital twins) to facilitate their translation into routine clinical applications.

Highlights of the Special Issue
The current Special Issue aims to fill knowledge gaps, introduce a variety of mathematical and computational frameworks, and highlight the state-of-the-art research on multi-scale and multi-physics models in therapeutic/diagnostic agent development for cancer, aiming to demonstrate their potential clinical impact.In this Special Issue, original research articles and reviews are invited for submission.The topics covered include, but are not limited to, mathematical/computational modeling of delivery systems for theranostics, targeted delivery and nanomedicine more broadly.
Eight manuscripts were submitted for consideration in this Special Issue, and each underwent a rigorous review process by Cancers.Ultimately, five papers were accepted for publication and inclusion in this Special Issue.Concise descriptions of the published articles are listed below: This study evaluates a novel treatment approach utilizing antiangiogenic and cytotoxic drugs, supported by mathematical modeling, revealing the potential benefits of this combined therapy in achieving a more uniform intratumoral environment and advancing the development of more effective brain tumor treatments.
Contribution 1 introduces a sophisticated numerical model to investigate drug delivery with an anti-angiogenesis effect, surpassing previous research complexity.The model simulates intravascular and interstitial fluid flow, considering non-Newtonian blood behavior and a dynamically adapting microvascular network.The model evaluates solute exposure and uniformity in tumors of various sizes using a convection-diffusion-reaction model, incorporating anti-angiogenesis simulations.The findings demonstrate that anti-angiogenic therapy reduces drug wash-out in the tumor periphery.Notably, the microvascular structure and modifications induced by vascular normalization significantly impact drug delivery quality, showing a 39% improvement in uniformity with specific modifications in tumors with a radius of 0.2 cm.
Intraperitoneal (IP) chemotherapy holds promise for treating peritoneal carcinomatosis, yet limited drug penetration into tumors remains a challenge.Contribution 2 employs a numerical model incorporating an actual tumor image with heterogeneous vasculature to investigate IP chemotherapy.The tumor's geometry is reconstructed using image processing techniques, and the model accounts for drug binding and uptake by cancer cells.After 60 min of IP treatment with doxorubicin (DOX), the area under the curve (AUC) for the average free drug concentration reaches 295.18 mol•m −3 •s −1 , and the half-width parameter W 1/2 , reflecting drug penetration, ranges from 0.11 to 0.14 mm.The treatment resulted in a fraction of killed cells reaching 20.4% by the end of the procedure.The findings emphasize the importance of considering specific vascular networks in modeling IP chemotherapy, showcasing the potential for patient-specific applications.
In contribution 3, a novel multiphysics simulation is developed to model LIPUSinduced doxorubicin (DOX) release from gold nanoparticle (GNP) drug carriers in an ex vivo tissue model.Transmission electronic microscopy imaging is performed before and after LIPUS exposure, and significant aggregation of the GNPs is observed upon DOX release.A numerical model is 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 total interaction potential is found before and after LIPUS exposure, and an energy barrier for aggregation is determined.It is concluded that the interaction potential energy threshold for GNP aggregation (and, as a result, DOX release) is equal to 0.24 k B T. This simulation model can predict whether DOX release will occur for a series of GNP and LIPUS parameters.An extensive parametric study is also performed with the simulation model to study how GNP and LIPUS parameters affect DOXloaded GNP colloidal stability.This model could optimize future ultrasound-triggered GNP-based drug delivery systems before clinical implementation.
Radiosensitizers enhance radiotherapy outcomes, and their effective distribution within the tumor is crucial due to their limited range.Contribution 4 introduces a computational model for nanoparticle transport within tumors, considering fluid velocity, particle deposition, and the convection-diffusion equation.Evaluating the impact of particle surface charge and injection locations reveals that negatively charged particles achieve a more uniform distribution, covering 100% of the tumor volume in the fluid and 44.5% of the deposited volume.Shifting the injection location towards the tumor's center enhances the particle distribution.The findings suggest that negatively charged radiosensitizing particles and strategic injection sites can optimize their spread and penetration within the tumor.
Contribution 5 presents an image-based computational study of interstitial fluid flow and transport of liposomal cytotoxic and antiangiogenic drugs using an integrated approach that couples the key drug delivery processes in convection-enhanced delivery (CED) treatments.The novel contribution is the development of an integrated drug transport model to predict the dynamic changes of interstitial fluid flow and drug transport upon anti-angiogenesis.More importantly, the tumor's 3D geometry and its heterogeneous characteristics, including its microvascular density, porosity, and cell density, are extracted from patient-specific dynamic MR images and applied as a model input to accommodate the realistic tumor properties.This study demonstrates that the co-infusion of antiangiogenic drugs can effectively improve cytotoxic drug delivery outcomes, facilitating enhanced drug penetration and accumulation.This can be attributed to the reduced microvascular density that inhibits the fluid exchange between the blood and tissue, mitigating cytotoxic drug dilution and drug loss due to blood drainage.

Conclusions and Future Research Need
Although submissions for this Special Issue are now closed, there is an ongoing need for sustained research and development in the realm of multi-scale and multi-physics models concerning the release, transport, and delivery of anticancer drugs (as well as diagnostic agents) to cancer.Considering the distinctive features of computational models, which include their widespread accessibility, cost-effectiveness, and the ability to integrate artificial intelligence, machine learning, and data inputs from imaging modalities, it is expected that multi-scale and multi-physics models for anticancer drug delivery will play a more prominent role in future cancer treatment strategies.
The primary hurdles in these models involve achieving a high resolution and accuracy in quantitative predictions while maintaining computational speed.Establishing a modeling framework and workflow that can concurrently address both objectives is imperative.The accuracy of these models is primarily dependent on parameterization related to the tumor, drug, nanoparticle, external source, and the physical models used to simulate the entire process.In practical terms, the intricacies of the drug delivery process require many model parameters, posing a distinct challenge in measuring all the parameter values.Validating computational models involving diverse physics introduces additional complexity.Simplifications and assumptions within the models are other factors that must be carefully verified before making final decisions.

Table 1 .
Hornsby, T.K.; Kashkooli, F.M.; Jakhmola, A.; Kolios, M.C.; Tavakkoli, J. Multiphysics Modeling of Low-Intensity Pulsed Ultrasound Induced Chemotherapeutic Drug Release from the Surface of Gold Nanoparticles.Cancers 2023, 15, 523.https://doi.The contents of this Special Issue are outlined in Table1, showcasing published papers spanning four distinct research domains: drug delivery, nano-sized drug delivery, drug release, and anti-angiogenesis.Table1provides in-depth insights into the diverse aspects of the physics, scales, and numerical approaches employed across these studies.Analysis of the published contributions in the Special Issue.