An Integrated Monte Carlo Model for Heterogeneous Glioblastoma Treated with Boron Neutron Capture Therapy

Glioblastomas (GBMs) are notorious for their aggressive features, e.g., intrinsic radioresistance, extensive heterogeneity, hypoxia, and highly infiltrative behaviours. The prognosis has remained poor despite recent advances in systemic and modern X-ray radiotherapy. Boron neutron capture therapy (BNCT) represents an alternative radiotherapy technique for GBM. Previously, a Geant4 BNCT modelling framework was developed for a simplified model of GBM.

The current work expands on the previous model by applying a more realistic in silico GBM model with heterogeneous radiosensitivity and anisotropic microscopic extensions (ME).

Each cell within the GBM model was assigned an α/β value associated with different GBM cell lines and a 10B concentration. Dosimetry matrices corresponding to various MEs were calculated and combined to evaluate cell survival fractions (SF) using clinical target volume (CTV) margins of 2.0 & 2.5 cm. SFs for the BNCT simulation were compared with external X-ray radiotherapy (EBRT) SFs.

The SFs within the beam region decreased by more than two times compared to EBRT. It was demonstrated that BNCT results in markedly reduced SFs for both CTV margins compared to EBRT. However, the SF reduction as a result of the CTV margin extension using BNCT was significantly lower than using X-ray EBRT for one MEP distribution, while it remained similar for the other two MEP models.

Although the efficiency of BNCT in terms of cell kill is superior to EBRT, the extension of the CTV margin by 0.5 cm may not increase the BNCT treatment outcome significantly.