No Arabic abstract
Dose painting of hypoxic tumour sub-volumes using positron-emission tomography (PET) has been shown to improve tumour control in silico in several sites. Pancreatic cancer presents a more stringent challenge, given its proximity to critical organs-at-risk (OARs) and anatomic motion. A radiobiological model was developed to estimate clonogen survival fraction (SF), using 18F-fluoroazomycin arabinoside PET (FAZA PET) images from ten patients with pancreatic cancer to quantify oxygen enhancement effects. For each patient, four simulated five-fraction stereotactic body radiotherapy (SBRT) plans were generated: 1) a standard SBRT plan aiming to cover the planning target volume with 40 Gy, 2) dose painting plans delivering escalated doses to FAZA-avid hypoxic sub-volumes, 3) dose painting plans with simulated spacer separating the duodenum and pancreatic head, and 4), plans with integrated boosts to geometric contractions of the tumour (GTV). All plans saturated at least one OAR dose limit. SF was calculated for each plan and sensitivity of SF to simulated hypoxia quantification errors was evaluated. Dose painting resulted in a 55% reduction in SF as compared to standard SBRT; 78% with spacer. Integrated boosts to hypoxia-blind geometric contractions resulted in a 41% reduction in SF. The reduction in SF for dose-painting plans persisted for all hypoxia quantification parameters studied, including registration and rigid motion errors that resulted in shifts and rotations of the GTV and hypoxic sub-volumes by as much as 1 cm and 10 degrees. Although proximity to OARs ultimately limited dose escalation, with estimated SFs (~10^-5) well above levels required to completely ablate a ~10 cm^3 tumour, dose painting robustly reduced clonogen survival when accounting for expected treatment and imaging uncertainties and thus, may improve local response and associated morbidity.
A new variant of the pencil-beam (PB) algorithm for dose distribution calculation for radiotherapy with protons and heavier ions, the grid-dose spreading (GDS) algorithm, is proposed. The GDS algorithm is intrinsically faster than conventional PB algorithms due to approximations in convolution integral, where physical calculations are decoupled from simple grid-to-grid energy transfer. It was effortlessly implemented to a carbon-ion radiotherapy treatment planning system to enable realistic beam blurring in the field, which was absent with the broad-beam (BB) algorithm. For a typical prostate treatment, the slowing factor of the GDS algorithm relative to the BB algorithm was 1.4, which is a great improvement over the conventional PB algorithms with a typical slowing factor of several tens. The GDS algorithm is mathematically equivalent to the PB algorithm for horizontal and vertical coplanar beams commonly used in carbon-ion radiotherapy while dose deformation within the size of the pristine spread occurs for angled beams, which was within 3 mm for a single proton pencil beam of $30^circ$ incidence, and needs to be assessed against the clinical requirements and tolerances in practical situations.
Tumor motion plays a key role in the safe delivery of Stereotactic Body Radiotherapy (SBRT) for pancreatic cancer. The purpose of this study was to use tumor motion data measured in patients to establish limits on motion magnitude for safe delivery of pancreatic SBRT. Using 91 sets of pancreatic tumor motion data measured in patients, we calculated motion-convolved dose for 25 pancreatic cancer patients, and established the maximum amount of motion allowable while satisfying error thresholds on key dose metrics. In our patient cohort, the mean [min-max] allowable motion for 33/40/50 Gy to the PTV was 11.9 [6.3-22.4], 10.4 [5.2-19.1] and 9.0 [4.2-16.0] mm, respectively. Maximum allowable motion decreased as dose was escalated, and was smaller in patients with larger tumors. The effects of motion on pancreatic SBRT are highly variable between patients and there is potential to allow more motion in certain patients, even in dose-escalated scenarios. In our dataset, a conservative limit of 6.3 mm would ensure safe treatment of all patients treated to 33 Gy in 5 fractions.
Stereotactic body radiation therapy (SBRT) for pancreatic cancer requires a skillful approach to deliver ablative doses to the tumor while limiting dose to the highly sensitive duodenum, stomach, and small bowel. Here, we develop knowledge-based artificial neural network dose models (ANN-DMs) to predict dose distributions that would be approved by experienced physicians. Using dose distributions calculated by a commercial treatment planning system (TPS), physician-approved treatment plans were used to train ANN-DMs that could predict physician-approved dose distributions based on a set of geometric parameters (vary from voxel to voxel) and plan parameters (constant across all voxels for a given patient). Differences between TPS and ANN-DM dose distributions were used to evaluate model performance. ANN-DM design, including neural network structure and parameter choices, were evaluated to optimize dose model performance. Mean dose errors were less than 5% at all distances from the PTV, and mean absolute dose errors were on the order of 5%, but no more than 10%. Dose-volume histogram errors demonstrated good model performance above 25 Gy, but much larger errors were seen at lower doses. ANN-DM dose distributions showed excellent overall agreement with TPS dose distributions, and accuracy was substantially improved when each physicians treatment approach was taken into account by training their own dedicated models. In this manner, one could feasibly train ANN-DMs that could predict the dose distribution desired by a given physician for a given treatment site.
In carbon-ion radiotherapy, single-beam delivery each day in alternate directions has been commonly practiced for operational efficiency, taking advantage of the Bragg peak and the relative biological effectiveness (RBE) for uniform dose conformation to a tumor. The treatment plans are usually evaluated with total RBE-weighted dose, which is however deficient in relevance to the biological effect in the linear-quadratic model due to its quadratic-dose term, or the dose-fractionation effect. In this study, we reformulate the extrapolated response dose (ERD), or synonymously BED, which normalizes the dose-fractionation and cell-repopulation effects as well as the RBE of treating radiation, based on inactivation of a single model cell system and a typical treating radiation in carbon-ion RT. The ERD distribution virtually represents the biological effect of the treatment regardless of radiation modality or fractionation scheme. We applied the ERD formulation to simplistic model treatments and to a preclinical survey for hypofractionation based on an actual prostate-cancer treatment of carbon-ion radiotherapy. The proposed formulation was demonstrated to be practical and to offer theoretical implications. In the prostate-cancer case, the ERD distribution was very similar to the RBE-weighted-dose distribution of the actual treatment in 12 fractions. With hypofractionation, while the RBE-weighted-dose distribution varied significantly, the ERD distribution was nearly invariant, implying that the carbon-ion radiotherapy would be insensitive to fractionation. However, treatment evaluation with simplistic biological dose is intrinsically limited and must be complemented in practice somehow by clinical experiences and biology experiments.
Pancreas stereotactic body radiotherapy treatment planning requires planners to make sequential, time consuming interactions with the treatment planning system (TPS) to reach the optimal dose distribution. We seek to develop a reinforcement learning (RL)-based planning bot to systematically address complex tradeoffs and achieve high plan quality consistently and efficiently. The focus of pancreas SBRT planning is finding a balance between organs-at-risk sparing and planning target volume (PTV) coverage. Planners evaluate dose distributions and make planning adjustments to optimize PTV coverage while adhering to OAR dose constraints. We have formulated such interactions between the planner and the TPS into a finite-horizon RL model. First, planning status features are evaluated based on human planner experience and defined as planning states. Second, planning actions are defined to represent steps that planners would commonly implement to address different planning needs. Finally, we have derived a reward system based on an objective function guided by physician-assigned constraints. The planning bot trained itself with 48 plans augmented from 16 previously treated patients and generated plans for 24 cases in a separate validation set. All 24 bot-generated plans achieve similar PTV coverages compared to clinical plans while satisfying all clinical planning constraints. Moreover, the knowledge learned by the bot can be visualized and interpreted as consistent with human planning knowledge, and the knowledge maps learned in separate training sessions are consistent, indicating reproducibility of the learning process.