Computational modelling of the therapeutic outputs of photodynamic therapy on spheroid-on-chip models

Abstract

Photodynamic therapy (PDT) is a medical radio chemotherapeutic method that uses light, photosensitizing agents, and oxygen to produce cytotoxic compounds, which eliminate malignant cells. Recently, Microfluidic systems have been used to analyse photosensitizers (PSs) due to their potential to replicate in vivo environments. While prior studies have established a strong correlation between reacted singlet oxygen concentration and PDT-induced cellular death, the effects that the ambient fluid flow might have on the concentration of oxygen and PS have been disregarded in many, which limits the reliability of the results. Herein, we coupled the transport of oxygen and PS throughout the ambient medium and within the spheroidal multicellular aggregate to initially study the profiles of oxygen and PS concentration alongside PDT-induced cellular death throughout the spheroid before and after radiation. The attained results indicate that the PDT-induced cellular death initiates on the surface of the spheroids and subsequently spreads to the neighbouring regions, which is in great accordance with experimental results. Afterward, the effects that drug-light interval (DLI), fluence rate, PS composition, microchannel height, and inlet flow rate have on the therapeutic outcomes are studied. The findings show that adequate DLI is critical to ensure uniform distribution of PS throughout the medium, and a value of 5 h was found to be sufficient. The composition of PS is critical, as ALA-PpIX induces earlier cell death but accelerates oxygen consumption, especially in the outer layers, depriving the inner layers of oxygen necessary for PDT, which in turn disrupts and prolongs the exposure time compared to mTHPC and Photofrin. Despite the fluence rate directly influencing the singlet oxygen generation rate, increasing the fluence rate by 189 mW/cm² would not significantly benefit us. Microwell height and inlet flow rate involve competing phenomena—increasing height or decreasing flow reduces oxygen supply and increases PS “washout” and its concentration.