Development of an integrated Monte Carlo model for glioblastoma multiforme treated with boron neutron capture therapy

Glioblastomas (GBM) are notorious for their high fatality rate. Boron Neutron Capture Therapy (BNCT) being a biochemically targeted type of radiotherapy is a potent modality for GBM. In the current work, a BNCT treatment modelling framework for GBM was developed. Optimal Clinical Target Volume (CTV) margins for GBM-BNCT and the BNCT efficacy have been investigated. The model integrated a cell-based dosimetry model, an in-house-developed epithermal neutron beam model and previously-developed Microscopic Extension Probability (MEP) model. The system was defined as a cubic ICRP-brain phantom divided into 20 μm side voxels. The corresponding ¹⁰B concentrations in GBM and normal brain cells were applied. The in-silico model was irradiated with the epithermal neutron beam using 2 and 2.5 cm CTV margins. Results from the cell-based dosimetry and the MEP models were combined to calculate GBM cell survival fractions (SF) post BNCT and compared to x-ray radiotherapy (XRT) SFs. Compared to XRT, the SF within the beam decreased by five orders of magnitudes and the total SF was reduced three times following BNCT. CTV extension by 0.5 cm reduced the SF by additional (53.8 ± 0.3)%. In conclusion, BNCT results in a more efficient cell kill. The extension of the CTV margin, however, may not increase the treatment outcome significantly.