Introduction: Photothermal hyperthermia is a promising minimally invasive approach to treat surface skin cancers; however, for an effective therapeutic outcome, it is crucial to ensure adequate tumor heating without causing heat damage to the adjacent normal tissues. This paper presents a symmetry-based computational model to analyze microneedle-assisted photothermal cancer treatment using nanoparticles.
Methods: A Pennes-based bioheat equation that incorporates multiple layers in the biological tissues, tumor-dependent blood flow and metabolism, volumetric heat source due to lasers, localized heat generation due to nanoparticles, convective cooling on the surface, interface conditions, and the Arrhenius thermal damage function has been developed. The discretized time-dependent governing equation has been converted into an equivalent Helmholtz equation that is solved by a dual-reciprocity boundary element method. The localized deposition of nanoparticles due to dissolvable microneedles has been represented using Gaussian distribution functions.
Results: Numerical results have proven that factors such as the laser intensity, the concentration of the nanoparticles, the depth of the tumor, and the deposition profile play significant roles in affecting tumor-center temperature, heat dose accumulation, and the spatial selectivity. In contrast to traditional heating using the laser alone, the localized deposition of nanoparticles helps to destroy the thermal symmetry in the temperature distribution and achieve more selective tumor heating.
Conclusions: The newly introduced symmetry-based DRBEM approach serves as a very convenient and accurate instrument for prediction of temperature changes, assessment of thermal damage and parameter optimization for microneedle-based nanoparticle photothermal hyperthermia therapy. Symmetry breaking is found to play a crucial role in making tumors more selective to such therapy.