Numerical solution of an optimal control model of dendritic cell treatment of a growing tumour

A new optimal control model of the interactions between a growing tumour and the host immune system along with an immunotherapy treatment strategy is presented. The model is based on an ordinary differential equation model of interactions between the growing tu- mour and the natural killer, cytotoxic T lymphocyte and dendritic cells of the host immune system, extended through the addition of a control function representing the application of a dendritic cell treat- ment to the system. The numerical solution of this model, obtained from a multi species Runge–Kutta forward-backward sweep scheme, is described. We investigate the effects of varying the maximum al- lowed amount of dendritic cell vaccine administered to the system and find that control of the tumour cell population is best effected via a high initial vaccine level, followed by reduced treatment and finally cessation of treatment. We also found that increasing the strength of the dendritic cell vaccine causes an increase in the number of natural killer cells and lymphocytes, which in turn reduces the growth of the tumour.

Segmentation of cone-beam CT using a hidden Markov random field with informative priors

Purpose: Flat-detector, cone-beam computed tomography (CBCT) has enormous potential to improve the accuracy of treatment delivery in image-guided radiotherapy (IGRT). To assist radiotherapists in interpreting these images, we use a Bayesian statistical model to label each voxel according to its tissue type. Methods: The rich sources of prior information in IGRT are incorporated into a hidden Markov random field (MRF) model of the 3D image lattice. Tissue densities in the reference CT scan are estimated using inverse regression and then rescaled to approximate the corresponding CBCT intensity values. The treatment planning contours are combined with published studies of physiological variability to produce a spatial prior distribution for changes in the size, shape and position of the tumour volume and organs at risk (OAR). The voxel labels are estimated using the iterated conditional modes (ICM) algorithm. Results: The accuracy of the method has been evaluated using 27 CBCT scans of an electron density phantom (CIRS, Inc. model 062). The mean voxel-wise misclassification rate was 6.2%, with Dice similarity coefficient of 0.73 for liver, muscle, breast and adipose tissue. Conclusions: By incorporating prior information, we are able to successfully segment CBCT images. This could be a viable approach for automated, online image analysis in radiotherapy.

Developing an adaptive radiotherapy technique for virally mediated head and neck cancer

Purpose Virally mediated head and neck cancers (VMHNC) often present with nodal involvement and are highly radioresponsive, meaning that treatment plan adaptation during radiotherapy (RT) in a subset of patients is required. We sought to determine potential risk profiles and a corresponding adaptive treatment strategy for these patients. Methodology 121 patients with virally mediated, node positive nasopharyngeal (Epstein Barr Virus positive) or oropharyngeal (Human Papillomavirus positive) cancers, receiving curative intent RT were reviewed. The type, frequency and timing of adaptive interventions, including source-to-skin distance (SSD) corrections, re-scanning and re-planning, were evaluated. Patients were reviewed based on the maximum size of the dominant node to assess the need for plan adaptation. Results Forty-six patients (38%) required plan adaptation during treatment. The median fraction at which the adaptive intervention occurred was 26 for SSD corrections and 22 for re-planning CTs. A trend toward 3 risk profile groupings was discovered: 1) Low risk with minimal need (< 10%) for adaptive intervention (dominant pre-treatment nodal size of ≤ 35 mm), 2) Intermediate risk with possible need (< 20%) for adaptive intervention (dominant pre-treatment nodal size of 36 mm – 45 mm) and 3) High-risk with increased likelihood (> 50%) for adaptive intervention (dominant pre-treatment nodal size of ≥ 46 mm). Conclusion In this study, patients with VMHNC and a maximum dominant nodal size of > 46 mm were identified at a higher risk of requiring re-planning during a course of definitive RT. Findings will be tested in a future prospective adaptive RT study.

Evaluating the dosimetric effect of treatment-induced changed in virally mediated head and neck cancer patients

Introduction Patients with virally mediated head and neck cancer (VMHNC) often present with advanced nodal disease that is highly radioresponsive as demonstrated by tumour and nodal regression during treatment. The resultant changes may impact on the planned dose distribution and so adversely affect the therapeutic ratio. The aim of this study was to evaluate the dosimetric effect of treatment-induced anatomical changes in VMHNC patients who had undergone a re-plan. Methods Thirteen patients with virally mediated oropharyngeal or nasopharyngeal cancer who presented for definitive radiotherapy between 2005 and 2010 and who had a re-plan generated were investigated. The dosimetric effect of anatomical changes, was quantified by comparing dose volume histograms (DVH) of primary and nodal gross target volumes and organs at risk (OAR), including spinal cord and parotid glands, from the original plan and a comparison plan. Results Eleven 3DCRT and 2 IMRT plans were evaluated. Dose to the spinal cord and brainstem increased by 4.1% and 2.6%, respectively. Mean dose to the parotid glands also increased by 3.5%. In contrast, the dose received by 98% of the primary and nodal gross tumour volumes decreased by 0.15% and 0.3%, respectively when comparing the initial treatment plan to the comparison plan. Conclusion In this study, treatment-induced anatomical changes had the greatest impact on OAR dose with negligible effect on the dose to nodal gross tumour volumes. In the era of intensity modulated radiotherapy (IMRT), accounting for treatment-induced anatomical changes is important as focus is placed on minimising the acute and long-term side effects of treatment.

The effect of organ motion of radiotherapy plans

Organ motion as a result of respiration is an important field of research for medical physics. Knowledge of magnitude and direction of this motion is necessary to allow for more accurate radiotherapy treatment planning. This will result in higher doses to the tumour whilst sparing healthy tissue. This project involved human trials, where the radiation therapy patient's kidneys were CT scanned under three different conditions; whilst free breathing (FB), breath-hold at normal tidal inspiration (BHIN), and breath-hold at normal tidal expiration (BHEX). The magnitude of motion was measured by recording the outline of the kidney from a Beam's Eye View (BEV). The centre of mass of this 2D shape was calculated for each set using "ImageJ" software and the magnitude of movement determined from the change in the centroid's coordinates between the BHIN and BHEX scans. The movement ranged from, for the left and right kidneys, 4-46mm and 2-44mm in the superior/inferior (axial) plane, 1-21mm and 2- 16mm in the anterior/posterior (coronal) plane, and 0-6mm and 0-8mm in the lateral/medial (sagittal) plane. From exhale to inhale, the kidneys tended to move inferiorly, anteriorly and laterally. A standard radiotherapy plan, designed to treat the para-aortics with opposed lateral fields was performed on the free breathing (planning) CT set. The field size and arrangement was set up using the same parameters for each subject. The prescription was to deliver 45 Gray in 25 fractions. This field arrangement and prescription was then copied over to the breath hold CT sets, and the dosimetric differences were compared using Dose Volume Histograms (DVH). The point of comparison for the three sets was recorded as the percentage volume of kidney receiving less than or equal to 10 Gray. The QUASAR respiratory motion phantom was used with the range of motion determined from the human study. The phantom was imaged, planned and treated with a linear accelerator with dose determined by film. The effect of the motion was measured by the change in the penumbra of the film and compared to the penumbra from the treatment planning system.

A feasibility study of multislice x-ray CT imaging of gel dosimeters using the “zero scan” method

This study extends the ‘zero scan’ method for CT imaging of polymer gel dosimeters to include multi-slice acquisitions. Multi slice CT images consisting of 24 slices of 1.2 mm thickness were acquired of an irradiated polymer gel dosimeter, and processed with the zero scan technique. The results demonstrate that zero scan based gel readout can be successfully applied to generate a three dimensional image of the irradiated gel field. Compared to the raw CT images the processed figures and cross gel profiles demonstrated reduced noise and clear visibility of the penumbral region. Moreover these improved results further highlight the suitability of this method in volumetric reconstruction with reduced CT data acquisition per slice. This work shows that 3D volumes of irradiated polymer gel dosimeters can be acquired and processed with x-ray CT.

Constructing a clinical decision-making framework for image-guided radiotherapy using a Bayesian Network

A decision-making framework for image-guided radiotherapy (IGRT) is being developed using a Bayesian Network (BN) to graphically describe, and probabilistically quantify, the many interacting factors that are involved in this complex clinical process. Outputs of the BN will provide decision-support for radiation therapists to assist them to make correct inferences relating to the likelihood of treatment delivery accuracy for a given image-guided set-up correction. The framework is being developed as a dynamic object-oriented BN, allowing for complex modelling with specific sub-regions, as well as representation of the sequential decision-making and belief updating associated with IGRT. A prototype graphic structure for the BN was developed by analysing IGRT practices at a local radiotherapy department and incorporating results obtained from a literature review. Clinical stakeholders reviewed the BN to validate its structure. The BN consists of a sub-network for evaluating the accuracy of IGRT practices and technology. The directed acyclic graph (DAG) contains nodes and directional arcs representing the causal relationship between the many interacting factors such as tumour site and its associated critical organs, technology and technique, and inter-user variability. The BN was extended to support on-line and off-line decision-making with respect to treatment plan compliance. Following conceptualisation of the framework, the BN will be quantified. It is anticipated that the finalised decision-making framework will provide a foundation to develop better decision-support strategies and automated correction algorithms for IGRT.

Retrospective evaluation of dosimetric quality for prostate carcinomas treated with 3D conformal, intensity-modulated and volumetric-modulated arc radiotherapy

Introduction This study examines and compares the dosimetric quality of radiotherapy treatment plans for prostate carcinoma across a cohort of 163 patients treated across 5 centres: 83 treated with three-dimensional conformal radiotherapy (3DCRT), 33 treated with intensity-modulated radiotherapy (IMRT) and 47 treated with volumetric-modulated arc therapy (VMAT). Methods Treatment plan quality was evaluated in terms of target dose homogeneity and organ-at-risk sparing, through the use of a set of dose metrics. These included the mean, maximum and minimum doses; the homogeneity and conformity indices for the target volumes; and a selection of dose coverage values that were relevant to each organ-at-risk. Statistical significance was evaluated using two-tailed Welch’s T-tests. The Monte Carlo DICOM ToolKit software was adapted to permit the evaluation of dose metrics from DICOM data exported from a commercial radiotherapy treatment planning system. Results The 3DCRT treatment plans offered greater planning target volume dose homogeneity than the other two treatment modalities. The IMRT and VMAT plans offered greater dose reduction in the organs-at-risk: with increased compliance with recommended organ-at-risk dose constraints, compared to conventional 3DCRT treatments. When compared to each other, IMRT and VMAT did not provide significantly different treatment plan quality for like-sized tumour volumes. Conclusions This study indicates that IMRT and VMAT have provided similar dosimetric quality, which is superior to the dosimetric quality achieved with 3DCRT.

A simple method for EPID-based in vivo dosimetry for radiotherapy treatments of head-and-neck cancers

This study used a homogeneous water-equivalent model of an electronic portal imaging device (EPID), contoured as a structure in a radiotherapy treatment plan, to produce reference dose images for comparison with in vivo EPID dosimetry images. Head and neck treatments were chosen as the focus of this study, due to the heterogeneous anatomies involved and the consequent difficulty of rapidly obtaining reliable reference dose images by other means. A phantom approximating the size and heterogeneity of a typical neck, with a maximum radiological thickness of 8.5 cm, was constructed for use in this study. This phantom was CT scanned and a simple treatment including five square test fields and one off-axis IMRT field was planned. In order to allow the treatment planning system to calculate dose in a model EPID positioned a distance downstream from the phantom to achieve a source-to-detector distance (SDD) of 150 cm, the CT images were padded with air and the phantom’s “body” contour was extended to encompass the EPID contour. Comparison of dose images obtained from treatment planning calculations and experimental irradiations showed good agreement, with more than 90% of points in all fields passing a gamma evaluation, at γ (3%, 3mm )Similar agreement was achieved when the phantom was over-written with air in the treatment plan and removed from the experimental beam, suggesting that water EPID model at 150 cm SDD is capable of providing accurate reference images for comparison with clinical IMRT treatment images, for patient anatomies with radiological thicknesses ranging from 0 up to approximately 9 cm. This methodology therefore has the potential to be used for in vivo dosimetry during treatments to tissues in the neck as well as the oral and nasal cavities, in the head-and-neck region.

Treatment plan complexity metrics for predicting IMRT pre-treatment quality assurance results

The planning of IMRT treatments requires a compromise between dose conformity (complexity) and deliverability. This study investigates established and novel treatment complexity metrics for 122 IMRT beams from prostate treatment plans. The Treatment and Dose Assessor software was used to extract the necessary data from exported treatment plan files and calculate the metrics. For most of the metrics, there was strong overlap between the calculated values for plans that passed and failed their quality assurance (QA) tests. However, statistically significant variation between plans that passed and failed QA measurements was found for the established modulation index and for a novel metric describing the proportion of small apertures in each beam. The ‘small aperture score’ provided threshold values which successfully distinguished deliverable treatment plans from plans that did not pass QA, with a low false negative rate.

Investigating the sensitivity of OSLD responses for varying field sizes across matched linear accelerators

When radiation therapy centres are equipped with two or more linear accelerators from the same vendor, they are usually beam-matched. This work tested the sensitivity of optically stimulated luminescence dosimeters (OSLDs) across matched linear accelerators. The responses were compared with an unshielded diode detector for varying field sizes. Clinical studies are currently done with thermoluminescent dosimeters (TLD), which absorb radiation then emit some levels of light determined by the radiation absorption when heated.

Anti-proliferative and cytotoxic activity of pentadactylin isolated from Leptodactylus labyrinthicus on melanoma cells

Nowadays, the emergence of resistance to the current available chemotherapeutic drugs by cancer cells makes the development of new agents imperative. The skin secretion of amphibians is a natural rich source of antimicrobial peptides (AMP), and researchers have shown that some of these wide spectrum molecules are also toxic to cancer cells. The aim of this study was to verify a putative anticancer activity of the AMP pentadactylin isolated for the first time from the skin secretion of the frog Leptodactylus labyrinthicus and also to study its cytotoxic mechanism to the murine melanoma cell line B16F10. The results have shown that pentadactylin reduces the cell viability of B16F10 cells in a dose-dependent manner. It was also cytotoxic to normal human fibroblast cells; nevertheless, pentadactylin was more potent in the first case. The studies of action mechanism revealed that pentadactylin causes cell morphology alterations (e.g., round shape and shrinkage morphology), membrane disruption, DNA fragmentation, cell cycle arrest at the S phase, and alteration of mitochondrial membrane potential, suggesting that B16F10 cells die by apoptosis. The exact mechanism that causes reduction of cell viability and cytotoxicity after treatment with pentadactylin is still unknown. In conclusion, as cancer cells become resilient to death, it is worthwhile the discovery of new drugs such as pentadactylin that induces apoptosis.

The cost-effectiveness of total laparoscopic hysterectomy compared to total abdominal hysterectomy for the treatment of early stage endometrial cancer

Objective To summarise how costs and health benefits will change with the adoption of total laparoscopic hysterectomy compared to total abdominal hysterectomy for the treatment of early stage endometrial cancer. Design Cost-effectiveness modelling using the information from a randomised controlled trial. Participants Two hypothetical modelled cohorts of 1000 individuals undergoing total laparoscopic hysterectomy and total abdominal hysterectomy. Outcome measures Surgery costs; hospital bed days used; total healthcare costs; quality-adjusted life years; and net monetary benefits. Results For 1000 individuals receiving total laparoscopic hysterectomy surgery, the costs were $509 575 higher, 3548 hospital fewer bed days were used and total health services costs were reduced by $3 746 221. There were 39.13 more quality-adjusted life years for a 5 year period following surgery. Conclusions The adoption of total laparoscopic hysterectomy is almost certainly a good decision for health services policy makers. There is 100% probability that it will be cost saving to health services, a 86.8% probability that it will increase health benefits and a 99.5% chance that it returns net monetary benefits greater than zero.

Transforming the future of treatment for ovarian cancer

As for many other cancers, metastasis is the leading cause of death of patients with ovarian cancer. Vigorous basic and clinical research is being performed to initiate more efficacious treatment strategies to improve the poor outcome of women with this cancer. Current treatment for ovarian cancer includes advanced cyto-reductive surgery and traditional platinum and taxane combined chemotherapy. Clinical trials using novel cytotoxic reagents and tyrosine kinase inhibitors have also been progressing. In parallel, the application of robust unbiased high throughput research platforms using transcriptomic and proteomic approaches has identified that not only individual cell signalling pathways, but a network of molecular pathways, play an important role in the biology of ovarian cancer. Furthermore, intensive genomic and epigenetic analyses have also revealed single nucleotide polymorphisms associated with risk and/or aetiology of this cancer including patient response to treatment. Taken together, these approaches, that are advancing our understanding, will have an impact on the generation of new therapeutic approaches and strategies for improving the outcome and quality of life of patients with ovarian cancer in the near future.