Transcriptomic profiling reveals differential cellular response to copper oxide nanoparticles and polystyrene nanoplastics in perfused human placenta.

The growing nanoparticulate pollution (e.g. engineered nanoparticles (NPs) or nanoplastics) has been shown to pose potential threats to human health. In particular, sensitive populations such as pregnant women and their unborn children need to be protected from harmful environmental exposures. However, developmental toxicity from prenatal exposure to pollution particles is not yet well studied despite evidence of particle accumulation in human placenta. Our study aimed to investigate how copper oxide NPs (CuO NPs; 10-20 nm) and polystyrene nanoplastics (PS NPs; 70 nm) impact on gene expression in ex vivo perfused human placental tissue. Whole genome microarray analysis revealed changes in global gene expression profile after 6 h of perfusion with sub-cytotoxic concentrations of CuO (10 µg/mL) and PS NPs (25 µg/mL). Pathway and gene ontology enrichment analysis of the differentially expressed genes suggested that CuO and PS NPs trigger distinct cellular response in placental tissue. While CuO NPs induced pathways related to angiogenesis, protein misfolding and heat shock responses, PS NPs affected the expression of genes related to inflammation and iron homeostasis. The observed effects on protein misfolding, cytokine signaling, and hormones were corroborated by western blot (accumulation of polyubiquitinated proteins) or qPCR analysis. Overall, the results of the present study revealed extensive and material-specific interference of CuO and PS NPs with placental gene expression from a single short-term exposure which deserves increasing attention. In addition, the placenta, which is often neglected in developmental toxicity studies, should be a key focus in the future safety assessment of NPs in pregnancy.

The Use of Episodic Memory Tests for Screening in Clinical Trials for Early Alzheimer's Disease: A Comparison of the Free and Cued Selective Reminding Test (FCSRT) and the Repeatable Battery for the Assessment of Neuropsychological Status (RBANS).

BACKGROUND: Screening procedures for early Alzheimer's disease (AD) trials seek to efficiently identify participants who fulfill clinical and biomarker criteria for AD and enrich for those most likely to experience significant clinical progression during the study. Episodic memory performance is often assessed in screening, but the utility of different memory tests for optimizing screening efficiency and/or rates of clinical progression remains uncertain. OBJECTIVES: Cross-study comparisons of the effects of inclusion criteria based on performance on the Free and Cued Selective Reminding Test (FCSRT) or the Repeatable Battery for the Assessment of Neuropsychological Status (RBANS) on screen-failure rates for episodic memory and ß-amyloid (Aß) positivity (by CSF or PET) and on subsequent rates of clinical disease progression in randomized participants across three clinical trials in early (prodromal-to-mild) AD. DESIGN: Secondary analyses of cross-sectional and longitudinal clinical trial data. SETTING: Multi-center international clinical trials. PARTICIPANTS: Individuals with prodromal-to-mild AD screened and/or randomized in clinical trials for crenezumab (CREAD, CREAD2) or semorinemab (Tauriel). Cross-sectional analyses of screening data for episodic memory impairment included participants from CREAD2 (n=2897) and Tauriel (n=887) and for Aß positivity included participants from CREAD (n=1138), CREAD2 (n=1119), and Tauriel (n=483). Longitudinal analyses of rates of clinical progression included participants from CREAD (n=779), CREAD2 (n=773), and Tauriel (n=331). MEASUREMENTS: Cross-sectional analyses examined eligibility rates per cutoffs defined for the FCSRT (CREAD, CREAD2) or RBANS (Tauriel) and per Aß positivity using CSF and/or PET biomarkers. Longitudinal analyses examined rates of clinical progression on the Clinical Dementia Rating-Sum of Boxes (CDR-SB), the 13-item version of the Alzheimer's Disease Assessment Scale-Cognitive Subscale (ADAS-Cog13), and Alzheimer's Disease Cooperative Study-Activities of Daily Living Scale (ADCS-ADL). RESULTS: Lower rates of study eligibility per episodic memory criteria were seen with the FCSRT (CREAD2) relative to the RBANS (Tauriel), but similar rates of eligibility per Aß positivity criteria were seen amongst participants with episodic memory impairment per the cutoffs used on either assessment. Similar rates of clinical decline over 18 months on the CDR-SB, ADAS-Cog13, and ADCS-ADL were observed in study populations enriched using the FCSRT (CREAD, CREAD2) or the RBANS (Tauriel). CONCLUSIONS: Cutoffs for episodic memory impairment on the FCSRT used in the CREAD and CREAD2 studies are more stringent than those on the RBANS used in the Tauriel study, resulting in lower rates of eligibility. However, given that study enrichment with either test yields similar rates of Aß positivity and clinical progression, considerations beyond these factors may drive the decision of which assessment to use for screening in early AD clinical trials.

Auto-commissioning of a Monte Carlo electron beam model with application to photon MLC shaped electron fields.

Objective.Presently electron beam treatments are delivered using dedicated applicators. An alternative is the usage of the already installed photon multileaf collimator (pMLC) enabling efficient electron treatments. Currently, the commissioning of beam models is a manual and time-consuming process. In this work an auto-commissioning procedure for the Monte Carlo (MC) beam model part representing the beam above the pMLC is developed for TrueBeam systems with electron energies from 6 to 22 MeV.Approach.The analytical part of the electron beam model includes a main source representing the primary beam and a jaw source representing the head scatter contribution each consisting of an electron and a photon component, while MC radiation transport is performed for the pMLC. The auto-commissioning of this analytical part relies on information pre-determined from MC simulations, in-air dose profiles and absolute dose measurements in water for different field sizes and source to surface distances (SSDs). For validation calculated and measured dose distributions in water were compared for different field sizes, SSDs and beam energies for eight TrueBeam systems. Furthermore, a sternum case in an anthropomorphic phantom was considered and calculated and measured dose distributions were compared at different SSDs.Main results.Instead of the manual commissioning taking up to several days of calculation time and several hours of user time, the auto-commissioning is carried out in a few minutes. Measured and calculated dose distributions agree generally within 3% of maximum dose or 2 mm. The gamma passing rates for the sternum case ranged from 96% to 99% (3% (global)/2 mm criteria, 10% threshold).Significance.The auto-commissioning procedure was successfully implemented and applied to eight TrueBeam systems. The newly developed user-friendly auto-commissioning procedure allows an efficient commissioning of an MC electron beam model and eases the usage of advanced electron radiotherapy utilizing the pMLC for beam shaping.

Efficiency enhancements of a Monte Carlo beamlet based treatment planning process: implementation and parameter study.

Objective.The computational effort to perform beamlet calculation, plan optimization and final dose calculation of a treatment planning process (TPP) generating intensity modulated treatment plans is enormous, especially if Monte Carlo (MC) simulations are used for dose calculation. The goal of this work is to improve the computational efficiency of a fully MC based TPP for static and dynamic photon, electron and mixed photon-electron treatment techniques by implementing multiple methods and studying the influence of their parameters.Approach.A framework is implemented calculating MC beamlets efficiently in parallel on each available CPU core. The user can specify the desired statistical uncertainty of the beamlets, a fractional sparse dose threshold to save beamlets in a sparse format and minimal distances to the PTV surface from which 2 × 2 × 2 = 8 (medium) or even 4 × 4 × 4 = 64 (large) voxels are merged. The compromise between final plan quality and computational efficiency of beamlet calculation and optimization is studied for several parameter values to find a reasonable trade-off. For this purpose, four clinical and one academic case are considered with different treatment techniques.Main results.Setting the statistical uncertainty to 5% (photon beamlets) and 15% (electron beamlets), the fractional sparse dose threshold relative to the maximal beamlet dose to 0.1% and minimal distances for medium and large voxels to the PTV to 1 cm and 2 cm, respectively, does not lead to substantial degradation in final plan quality compared to using 2.5% (photon beamlets) and 5% (electron beamlets) statistical uncertainty and no sparse format nor voxel merging. Only OAR sparing is slightly degraded. Furthermore, computation times are reduced by about 58% (photon beamlets), 88% (electron beamlets) and 96% (optimization).Significance.Several methods are implemented improving computational efficiency of beamlet calculation and plan optimization of a fully MC based TPP without substantial degradation in final plan quality.

Enabling non-isocentric dynamic trajectory radiotherapy by integration of dynamic table translations.

Objective.The purpose of this study is to develop a treatment planning process (TPP) for non-isocentric dynamic trajectory radiotherapy (DTRT) using dynamic gantry rotation, collimator rotation, table rotation, longitudinal, vertical and lateral table translations and intensity modulation and to validate the dosimetric accuracy.Approach.The TPP consists of two steps. First, a path describing the dynamic gantry rotation, collimator rotation and dynamic table rotation and translations is determined. Second, an optimization of the intensity modulation along the path is performed. We demonstrate the TPP for three use cases. First, a non-isocentric DTRT plan for a brain case is compared to an isocentric DTRT plan in terms of dosimetric plan quality and delivery time. Second, a non-isocentric DTRT plan for a craniospinal irradiation (CSI) case is compared to a multi-isocentric intensity modulated radiotherapy (IMRT) plan. Third, a non-isocentric DTRT plan for a bilateral breast case is compared to a multi-isocentric volumetric modulated arc therapy (VMAT) plan. The non-isocentric DTRT plans are delivered on a TrueBeam in developer mode and their dosimetric accuracy is validated using radiochromic films.Main results.The non-isocentric DTRT plan for the brain case is similar in dosimetric plan quality and delivery time to the isocentric DTRT plan but is expected to reduce the risk of collisions. The DTRT plan for the CSI case shows similar dosimetric plan quality while reducing the delivery time by 45% in comparison with the IMRT plan. The DTRT plan for the breast case showed better treatment plan quality in comparison with the VMAT plan. The gamma passing rates between the measured and calculated dose distributions are higher than 95% for all three plans.Significance.The versatile benefits of non-isocentric DTRT are demonstrated with three use cases, namely reduction of collision risk, reduced setup and delivery time and improved dosimetric plan quality.

A hybrid column generation and simulated annealing algorithm for direct aperture optimization.

The purpose of this work was to develop a hybrid column generation (CG) and simulated annealing (SA) algorithm for direct aperture optimization (H-DAO) and to show its effectiveness in generating high quality treatment plans for intensity modulated radiation therapy (IMRT) and mixed photon-electron beam radiotherapy (MBRT). The H-DAO overcomes limitations of the CG-DAO with two features improving aperture selection (branch-feature) and enabling aperture shape changes during optimization (SA-feature). The H-DAO algorithm iteratively adds apertures to the plan. At each iteration, a branch is created for each field provided. First, each branch determines the most promising aperture of its assigned field and adds it to a copy of the current apertures. Afterwards, the apertures of each branch undergo an MU-weight optimization followed by an SA-based simultaneous shape and MU-weight optimization and a second MU-weight optimization. The next H-DAO iteration continues the branch with the lowest objective function value. IMRT and MBRT treatment plans for an academic, a brain and a head and neck case generated using the CG-DAO and H-DAO were compared. For every investigated case and both IMRT and MBRT, the H-DAO leads to a faster convergence of the objective function value with number of apertures compared to the CG-DAO. In particular, the H-DAO needs about half the apertures to reach the same objective function value as the CG-DAO. The average aperture areas are 27% smaller for H-DAO than for CG-DAO leading to a slightly larger discrepancy between optimized and final dose. However, a dosimetric benefit remains. The H-DAO was successfully developed and applied to IMRT and MBRT. The faster convergence with number of apertures of the H-DAO compared to the CG-DAO allows to select a better compromise between plan quality and number of apertures.

TriB-RT: Simultaneous optimization of photon, electron and proton beams.

PURPOSE: To develop a novel treatment planning process (TPP) with simultaneous optimization of modulated photon, electron and proton beams for improved treatment plan quality in radiotherapy. METHODS: A framework for fluence map optimization of Monte Carlo (MC) calculated beamlet dose distributions is developed to generate treatment plans consisting of photon, electron and spot scanning proton fields. Initially, in-house intensity modulated proton therapy (IMPT) plans are compared to proton plans created by a commercial treatment planning system (TPS). A triple beam radiotherapy (TriB-RT) plan is generated for an exemplary academic case and the dose contributions of the three particle types are investigated. To investigate the dosimetric potential, a TriB-RT plan is compared to an in-house IMPT plan for two clinically motivated cases. Benefits of TriB-RT for a fixed proton beam line with a single proton field are investigated. RESULTS: In-house optimized IMPT are of at least equal or better quality than TPS-generated proton plans, and MC-based optimization shows dosimetric advantages for inhomogeneous situations. Concerning TriB-RT, for the academic case, the resulting plan shows substantial contribution of all particle types. For the clinically motivated case, improved sparing of organs at risk close to the target volume is achieved compared to IMPT (e.g. myelon and brainstem [Formula: see text] -37%) at cost of an increased low dose bath (healthy tissue V 10% +22%). In the scenario of a fixed proton beam line, TriB-RT plans are able to compensate the loss in degrees of freedom to substantially improve plan quality compared to a single field proton plan. CONCLUSION: A novel TPP which simultaneously optimizes photon, electron and proton beams was successfully developed. TriB-RT shows the potential for improved treatment plan quality and is especially promising for cost-effective single-room proton solutions with a fixed beamline in combination with a conventional linac delivering photon and electron fields.

Development of an extended Macro Monte Carlo method for efficient and accurate dose calculation in magnetic fields.

MOTIVATION: Progress in the field of magnetic resonance (MR)-guided radiotherapy has triggered the need for fast and accurate dose calculation in presence of magnetic fields. The aim of this work is to satisfy this need by extending the macro Monte Carlo (MMC) method to enable dose calculation for photon, electron, and proton beams in a magnetic field. METHODS: The MMC method is based on the transport of particles in macroscopic steps through an absorber by sampling the relevant physical quantities from a precalculated database containing probability distribution functions. To enable MMC particle transport in a magnetic field, a transformation accounting for the Lorentz force is applied for each macro step by rotating the sampled position and direction around the magnetic field vector. The transformed position and direction distributions on local geometries are validated against full MC for electron and proton pencil beams. To enable photon dose calculation, an in-house MC algorithm is used for photon transport and interaction. Emerging secondary charged particles are passed to MMC for transport and energy deposition. The extended MMC dose calculation accuracy and efficiency is assessed by comparison with EGSnrc (photon and electron beams) and Geant4 (proton beam) calculated dose distributions of different energies and homogeneous magnetic fields for broad beams impinging on water phantoms with bone and lung inhomogeneities. RESULTS: The geometric transformation on the local geometries is able to reproduce the results of full MC for all investigated settings (difference in mean value and standard deviation <1%). Macro Monte Carlo calculated dose distributions in a homogeneous magnetic field are in agreement with EGSnrc and Geant4, respectively, with gamma passing rates >99.6% (global 2%, 2 mm and 10% threshold criteria) for all situations. MMC achieves a substantial efficiency gain of up to a factor of 21 (photon beam), 66 (electron beam), and 356 (proton beam) compared to EGSnrc or Geant4. CONCLUSION: Efficient and accurate dose calculation in magnetic fields was successfully enabled by utilizing the developed extended MMC transport method for photon, electron, and proton beams.

Towards MR-guided electron therapy: Measurement and simulation of clinical electron beams in magnetic fields.

PURPOSE: In the current era of MRI-linac radiotherapy, dose optimization with arbitrary dose distributions is a reality. For the first time, we present new and targeted experiments and modeling to aid in evaluating the potential dose improvements offered with an electron beam mode during MRI-linac radiotherapy. METHODS: Small collimated (1 cm diameter and 1.5 × 1.5 cm2 square) electron beams (6, 12 and 20 MeV) from a clinical linear accelerator (Varian Clinac 2100C) are incident perpendicular and parallel to the strong and localized magnetic fields (0-0.7 T) generated by a permanent magnet device. Gafchromic EBT3 film is placed inside a slab phantom to measure two-dimensional dose distributions. A benchmarked and comprehensive Monte Carlo model (Geant4) is established to directly compare with experiments. RESULTS: With perpendicular fields a 5% narrowing of the beam FWHM and a 10 mm reduction in the 15% isodose penetration is seen for the 20 MeV beam. In the inline setup the penumbral width is reduced by up to 20%, and a local central dose enhancement of 100% is observed. Monte Carlo simulations are in agreement with the measured dose distributions (2% or 2 mm). CONCLUSION: A new range of experiments have been performed to offer insight into how an electron beam mode could offer additional choices in MRI-linac radiotherapy. The work extends on historic studies to bring a successful unified experimental and Monte Carlo modeling approach for studying small field electron beam dosimetry inside magnetic fields. The results suggest further work, particularly on the inline magnetic field scenario.

Benchmarking Monte-Carlo dose calculation for MLC CyberKnife treatments.

BACKGROUND: Vendor-independent Monte Carlo (MC) dose calculation (IDC) for patient-specific quality assurance of multi-leaf collimator (MLC) based CyberKnife treatments is used to benchmark and validate the commercial MC dose calculation engine for MLC based treatments built into the CyberKnife treatment planning system (Precision MC). METHODS: The benchmark included dose profiles in water in 15 mm depth and depth dose curves of rectangular MLC shaped fields ranging from 7.6 mm × 7.7 mm to 115.0 mm × 100.1 mm, which were compared between IDC, Precision MC and measurements in terms of dose difference and distance to agreement. Dose distributions of three phantom cases and seven clinical lung cases were calculated using both IDC and Precision MC. The lung PTVs ranged from 14 cm3 to 93 cm3. Quantitative comparison of these dose distributions was performed using dose-volume parameters and 3D gamma analysis with 2% global dose difference and 1 mm distance criteria and a global 10% dose threshold. Time to calculate dose distributions was recorded and efficiency was assessed. RESULTS: Absolute dose profiles in 15 mm depth in water showed agreement between Precision MC and IDC within 3.1% or 1 mm. Depth dose curves agreed within 2.3% / 1 mm. For the phantom and clinical lung cases, mean PTV doses differed from - 1.0 to + 2.3% between IDC and Precision MC and gamma passing rates were > =98.1% for all multiple beam treatment plans. For the lung cases, lung V20 agreed within ±1.5%. Calculation times ranged from 2.2 min (for 39 cm3 PTV at 1.0 × 1.0 × 2.5 mm3 native CT resolution) to 8.1 min (93 cm3 at 1.1 × 1.1 × 1.0 mm3), at 2% uncertainty for Precision MC for the 7 examined lung cases and 4-6 h for IDC, which, however, is not optimized for efficiency but used as a gold standard for accuracy. CONCLUSIONS: Both accuracy and efficiency of Precision MC in the context of MLC based planning for the CyberKnife M6 system were benchmarked against MC based IDC framework. Precision MC is used in clinical practice at our institute.

Investigation on the resolution of a micro cone beam CT scanner scintillating detector using Monte Carlo methods.

The impact of several physical quantities on the spatial resolution of an X-ray scintillating pixel detector for a micro cone beam CT (µCBCT) is investigated and discussed. The XtremeCT from SCANCO Medical AG was simulated using the EGSnrc/EGS++ Monte Carlo (MC) framework and extensively benchmarked in a previous work. The resolution of the detector was determined by simulating a titanium knife-edge to obtain the edge spread function (ESF) and the modulation transfer function (MTF). Propagation of the scintillation light through the scintillator and its coupling into the fiber optics system was taken into account. The contribution of particles scattered in the main scanner components to the detector signal is very low and does not affect the spatial resolution of the detector. The resolution obtained from the energy deposition in the scintillator without any blurring due to the propagation of the scintillation light into the fiber optics array was 31 µm. By assuming isotropic light propagation in the scintillator, the resolution degraded to 360 µm. A simple light propagation model taking into account the impact of the scintillator's columnar microstructures was developed and compared with the MANTIS Monte Carlo simulation package. By reducing the width of the model's light propagation kernel by a factor of 2 compared to the isotropic case, the detector resolution can be improved to 83 µm, which corresponds well to the measured resolution of 86 µm. The resolution of the detector is limited mainly by the propagation of the scintillation light through the scintillator layer. It offers the greatest potential to improve the resolution of the µCBCT imaging system.

Accuracy assessment of a potential clinical use of navigation-guided intra-operative liver metastasis brachytherapy-a planning study.

For patients with inoperable liver metastases, intra-operative liver high dose-rate brachytherapy (HDR-BT) is a promising technology enabling delivery of a high radiation dose to the tumor, while sparing healthy tissue. Liver brachytherapy has been described in the literature as safe and effective for the treatment of primary or secondary hepatic malignancies. It is preferred over other ablative techniques for lesions that are either larger than 4 cm or located in close proximity to large vessels or the common bile duct. In contrast to external beam radiation techniques, organ movements do not affect the size of the irradiated volume in intra-operative HDR-BT and new technical solutions exist to support image guidance for intra-operative HDR-BT. We have retrospectively analyzed anonymized CT datasets of 5 patients who underwent open liver surgery (resection and/or ablation) in order to test whether the accuracy of a new image-guidance method specifically adapted for intra-operative HDR-BT is high enough to use it in similar situations and whether patients could potentially benefit from navigation-guided intra-operative needle placement for liver HDR-BT.

Part 2: Dynamic mixed beam radiotherapy (DYMBER): Photon dynamic trajectories combined with modulated electron beams.

PURPOSE: The purpose of this study was to develop a treatment technique for dynamic mixed beam radiotherapy (DYMBER) utilizing increased degrees of freedom (DoF) of a conventional treatment unit including different particle types (photons and electrons), intensity and energy modulation and dynamic gantry, table, and collimator rotations. METHODS: A treatment planning process has been developed to create DYMBER plans combining photon dynamic trajectories (DTs) and step and shoot electron apertures collimated with the photon multileaf collimator (pMLC). A gantry-table path is determined for the photon DTs with minimized overlap of the organs at risk (OARs) with the target. In addition, an associated dynamic collimator rotation is established with minimized area between the pMLC leaves and the target contour. pMLC sequences of photon DTs and electron pMLC apertures are then simultaneously optimized using direct aperture optimization (DAO). Subsequently, the final dose distribution of the electron pMLC apertures is calculated using the Swiss Monte Carlo Plan (SMCP). The pMLC sequences of the photon DTs are then re-optimized with a finer control point resolution and with the final electron dose distribution taken into account. Afterwards, the final photon dose distribution is calculated also using the SMCP and summed together with the one of the electrons. This process is applied for a brain and two head and neck cases. The resulting DYMBER dose distributions are compared to those of dynamic trajectory radiotherapy (DTRT) plans consisting only of photon DTs and clinically applied VMAT plans. Furthermore, the deliverability of the DYMBER plans is verified in terms of dosimetric accuracy, delivery time and collision avoidance. For this purpose, The DYMBER plans are delivered to Gafchromic EBT3 films placed in an anthropomorphic head phantom on a Varian TrueBeam linear accelerator. RESULTS: For each case, the dose homogeneity in the target is similar or better for DYMBER compared to DTRT and VMAT. Averaged over all three cases, the mean dose to the parallel OARs is 16% and 28% lower, D2% to the serial OARs is 17% and 37% lower and V10% to normal tissue is 12% and 4% lower for the DYMBER plans compared to the DTRT and VMAT plans, respectively. The DYMBER plans are delivered without collision and with a 4-5 min longer delivery time than the VMAT plans. The absolute dose measurements are compared to calculation by gamma analysis using 2% (global)/2 mm criteria with passing rates of at least 99%. CONCLUSIONS: A treatment technique for DYMBER has been successfully developed and verified for its deliverability. The dosimetric superiority of DYMBER over DTRT and VMAT indicates utilizing increased DoF to be the key to improve brain and head and neck radiation treatments in future.

Electron beam collimation with a photon MLC for standard electron treatments.

Standard electron treatments are currently still performed using standard or molded patient-specific cut-outs placed in the electron applicator. Replacing cut-outs and electron applicators with a photon multileaf collimator (pMLC) for electron beam collimation would make standard electron treatments more efficient and would facilitate advanced treatment techniques like modulated electron radiotherapy (MERT) and mixed beam radiotherapy (MBRT). In this work, a multiple source Monte Carlo beam model for pMLC shaped electron beams commissioned at a source-to-surface distance (SSD) of 70 cm is extended for SSDs of up to 100 cm and validated for several Varian treatment units with field sizes typically used for standard electron treatments. Measurements and dose calculations agree generally within 3% of the maximal dose or 2 mm distance to agreement. To evaluate the dosimetric consequences of using pMLC collimated electron beams for standard electron treatments, pMLC-based and cut-out-based treatment plans are created for a left and a right breast boost, a sternum, a testis and a parotid gland case. The treatment plans consist of a single electron field, either alone (1E) or in combination with two 3D conformal tangential photon fields (1E2X). For each case, a pMLC plan with similar treatment plan quality in terms of dose homogeneity to the target and absolute mean dose values to the organs at risk (OARs) compared to a cut-out plan is found. The absolute mean dose to an OAR is slightly increased for pMLC-based compared to cut-out-based 1E plans if the OAR is located laterally close to the target with respect to beam direction, or if a 6 MeV electron beam is used at an extended SSD. In conclusion, treatment plans using cut-out collimation can be replaced by plans of similar treatment plan quality using pMLC collimation with accurately calculated dose distributions.

Independent Monte-Carlo dose calculation for MLC based CyberKnife radiotherapy.

This work aims to develop, implement and validate a Monte Carlo (MC)-based independent dose calculation (IDC) framework to perform patient-specific quality assurance (QA) for multi-leaf collimator (MLC)-based CyberKnife® (Accuray Inc., Sunnyvale, CA) treatment plans. The IDC framework uses an XML-format treatment plan as exported from the treatment planning system (TPS) and DICOM format patient CT data, an MC beam model using phase spaces, CyberKnife MLC beam modifier transport using the EGS++ class library, a beam sampling and coordinate transformation engine and dose scoring using DOSXYZnrc. The framework is validated against dose profiles and depth dose curves of single beams with varying field sizes in a water tank in units of cGy/Monitor Unit and against a 2D dose distribution of a full prostate treatment plan measured with Gafchromic EBT3 (Ashland Advanced Materials, Bridgewater, NJ) film in a homogeneous water-equivalent slab phantom. The film measurement is compared to IDC results by gamma analysis using 2% (global)/2 mm criteria. Further, the dose distribution of the clinical treatment plan in the patient CT is compared to TPS calculation by gamma analysis using the same criteria. Dose profiles from IDC calculation in a homogeneous water phantom agree within 2.3% of the global max dose or 1 mm distance to agreement to measurements for all except the smallest field size. Comparing the film measurement to calculated dose, 99.9% of all voxels pass gamma analysis, comparing dose calculated by the IDC framework to TPS calculated dose for the clinical prostate plan shows 99.0% passing rate. IDC calculated dose is found to be up to 5.6% lower than dose calculated by the TPS in this case near metal fiducial markers. An MC-based modular IDC framework was successfully developed, implemented and validated against measurements and is now available to perform patient-specific QA by IDC.

Simultaneous optimization of photons and electrons for mixed beam radiotherapy.

The aim of this work is to develop and investigate an inverse treatment planning process (TPP) for mixed beam radiotherapy (MBRT) capable of performing simultaneous optimization of photon and electron apertures. A simulated annealing based direct aperture optimization (DAO) is implemented to perform simultaneous optimization of photon and electron apertures, both shaped with the photon multileaf collimator (pMLC). Validated beam models are used as input for Monte Carlo dose calculations. Consideration of photon pMLC transmission during DAO and a weight re-optimization of the apertures after deliverable dose calculation are utilized to efficiently reduce the differences between optimized and deliverable dose distributions. The TPP for MBRT is evaluated for an academic situation with a superficial and an enlarged PTV in the depth, a left chest wall case including the internal mammary chain and a squamous cell carcinoma case. Deliverable dose distributions of MBRT plans are compared to those of modulated electron radiotherapy (MERT), photon IMRT and if available to those of clinical VMAT plans. The generated MBRT plans dosimetrically outperform the MERT, photon IMRT and VMAT plans for all investigated situations. For the clinical cases of the left chest wall and the squamous cell carcinoma, the MBRT plans cover the PTV similarly or more homogeneously than the VMAT plans, while OARs are spared considerably better with average reductions of the mean dose to parallel OARs and D 2% to serial OARs by 54% and 26%, respectively. Moreover, the low dose bath expressed as V 10% to normal tissue is substantially reduced by up to 45% compared to the VMAT plans. A TPP for MBRT including simultaneous optimization is successfully implemented and the dosimetric superiority of MBRT plans over MERT, photon IMRT and VMAT plans is demonstrated for academic and clinical situations including superficial targets with and without deep-seated part.

Assessing dose rate distributions in VMAT plans.

Dose rate is an essential factor in radiobiology. As modern radiotherapy delivery techniques such as volumetric modulated arc therapy (VMAT) introduce dynamic modulation of the dose rate, it is important to assess the changes in dose rate. Both the rate of monitor units per minute (MU rate) and collimation are varied over the course of a fraction, leading to different dose rates in every voxel of the calculation volume at any point in time during dose delivery. Given the radiotherapy plan and machine specific limitations, a VMAT treatment plan can be split into arc sectors between Digital Imaging and Communications in Medicine control points (CPs) of constant and known MU rate. By calculating dose distributions in each of these arc sectors independently and multiplying them with the MU rate, the dose rate in every single voxel at every time point during the fraction can be calculated. Independently calculated and then summed dose distributions per arc sector were compared to the whole arc dose calculation for validation. Dose measurements and video analysis were performed to validate the calculated datasets. A clinical head and neck, cranial and liver case were analyzed using the tool developed. Measurement validation of synthetic test cases showed linac agreement to precalculated arc sector times within ±0.4 s and doses ±0.1 MU (one standard deviation). Two methods for the visualization of dose rate datasets were developed: the first method plots a two-dimensional (2D) histogram of the number of voxels receiving a given dose rate over the course of the arc treatment delivery. In similarity to treatment planning system display of dose, the second method displays the dose rate as color wash on top of the corresponding computed tomography image, allowing the user to scroll through the variation over time. Examining clinical cases showed dose rates spread over a continuous spectrum, with mean dose rates hardly exceeding 100 cGy min(-1) for conventional fractionation. A tool to analyze dose rate distributions in VMAT plans with sub-second accuracy was successfully developed and validated. Dose rates encountered in clinical VMAT test cases show a continuous spectrum with a mean less than or near 100 cGy min(-1) for conventional fractionation.

Beamlet based direct aperture optimization for MERT using a photon MLC.

PURPOSE: A beamlet based direct aperture optimization (DAO) for modulated electron radiotherapy (MERT) using photon multileaf collimator (pMLC) shaped electron fields is developed and investigated. METHODS: The Swiss Monte Carlo Plan (SMCP) allows the calculation of dose distributions for pMLC shaped electron beams. SMCP is interfaced with the Eclipse TPS (Varian Medical Systems, Palo Alto, CA) which can thus be included into the inverse treatment planning process for MERT. This process starts with the import of a CT-scan into Eclipse, the contouring of the target and the organs at risk (OARs), and the choice of the initial electron beam directions. For each electron beam, the number of apertures, their energy, and initial shape are defined. Furthermore, the DAO requires dose-volume constraints for the structures contoured. In order to carry out the DAO efficiently, the initial electron beams are divided into a grid of beamlets. For each of those, the dose distribution is precalculated using a modified electron beam model, resulting in a dose list for each beamlet and energy. Then the DAO is carried out, leading to a set of optimal apertures and corresponding weights. These optimal apertures are now converted into pMLC shaped segments and the dose calculation for each segment is performed. For these dose distributions, a weight optimization process is launched in order to minimize the differences between the dose distribution using the optimal apertures and the pMLC segments. Finally, a deliverable dose distribution for the MERT plan is obtained and loaded back into Eclipse for evaluation. For an idealized water phantom geometry, a MERT treatment plan is created and compared to the plan obtained using a previously developed forward planning strategy. Further, MERT treatment plans for three clinical situations (breast, chest wall, and parotid metastasis of a squamous cell skin carcinoma) are created using the developed inverse planning strategy. The MERT plans are compared to clinical standard treatment plans using photon beams and the differences between the optimal and the deliverable dose distributions are determined. RESULTS: For the idealized water phantom geometry, the inversely optimized MERT plan is able to obtain the same PTV coverage, but with an improved OAR sparing compared to the forwardly optimized plan. Regarding the right-sided breast case, the MERT plan is able to reduce the lung volume receiving more than 30% of the prescribed dose and the mean lung dose compared to the standard plan. However, the standard plan leads to a better homogeneity within the CTV. The results for the left-sided thorax wall are similar but also the dose to the heart is reduced comparing MERT to the standard treatment plan. For the parotid case, MERT leads to lower doses for almost all OARs but to a less homogeneous dose distribution for the PTV when compared to a standard plan. For all cases, the weight optimization successfully minimized the differences between the optimal and the deliverable dose distribution. CONCLUSIONS: A beamlet based DAO using multiple beam angles is implemented and successfully tested for an idealized water phantom geometry and clinical situations.

Forward treatment planning for modulated electron radiotherapy (MERT) employing Monte Carlo methods.

PURPOSE: This paper describes the development of a forward planning process for modulated electron radiotherapy (MERT). The approach is based on a previously developed electron beam model used to calculate dose distributions of electron beams shaped by a photon multi leaf collimator (pMLC). METHODS: As the electron beam model has already been implemented into the Swiss Monte Carlo Plan environment, the Eclipse treatment planning system (Varian Medical Systems, Palo Alto, CA) can be included in the planning process for MERT. In a first step, CT data are imported into Eclipse and a pMLC shaped electron beam is set up. This initial electron beam is then divided into segments, with the electron energy in each segment chosen according to the distal depth of the planning target volume (PTV) in beam direction. In order to improve the homogeneity of the dose distribution in the PTV, a feathering process (Gaussian edge feathering) is launched, which results in a number of feathered segments. For each of these segments a dose calculation is performed employing the in-house developed electron beam model along with the macro Monte Carlo dose calculation algorithm. Finally, an automated weight optimization of all segments is carried out and the total dose distribution is read back into Eclipse for display and evaluation. One academic and two clinical situations are investigated for possible benefits of MERT treatment compared to standard treatments performed in our clinics and treatment with a bolus electron conformal (BolusECT) method. RESULTS: The MERT treatment plan of the academic case was superior to the standard single segment electron treatment plan in terms of organs at risk (OAR) sparing. Further, a comparison between an unfeathered and a feathered MERT plan showed better PTV coverage and homogeneity for the feathered plan, with V95% increased from 90% to 96% and V107% decreased from 8% to nearly 0%. For a clinical breast boost irradiation, the MERT plan led to a similar homogeneity in the PTV compared to the standard treatment plan while the mean body dose was lower for the MERT plan. Regarding the second clinical case, a whole breast treatment, MERT resulted in a reduction of the lung volume receiving more than 45% of the prescribed dose when compared to the standard plan. On the other hand, the MERT plan leads to a larger low-dose lung volume and a degraded dose homogeneity in the PTV. For the clinical cases evaluated in this work, treatment plans using the BolusECT technique resulted in a more homogenous PTV and CTV coverage but higher doses to the OARs than the MERT plans. CONCLUSIONS: MERT treatments were successfully planned for phantom and clinical cases, applying a newly developed intuitive and efficient forward planning strategy that employs a MC based electron beam model for pMLC shaped electron beams. It is shown that MERT can lead to a dose reduction in OARs compared to other methods. The process of feathering MERT segments results in an improvement of the dose homogeneity in the PTV.

Monte Carlo based beam model using a photon MLC for modulated electron radiotherapy.

PURPOSE: Modulated electron radiotherapy (MERT) promises sparing of organs at risk for certain tumor sites. Any implementation of MERT treatment planning requires an accurate beam model. The aim of this work is the development of a beam model which reconstructs electron fields shaped using the Millennium photon multileaf collimator (MLC) (Varian Medical Systems, Inc., Palo Alto, CA) for a Varian linear accelerator (linac). METHODS: This beam model is divided into an analytical part (two photon and two electron sources) and a Monte Carlo (MC) transport through the MLC. For dose calculation purposes the beam model has been coupled with a macro MC dose calculation algorithm. The commissioning process requires a set of measurements and precalculated MC input. The beam model has been commissioned at a source to surface distance of 70 cm for a Clinac 23EX (Varian Medical Systems, Inc., Palo Alto, CA) and a TrueBeam linac (Varian Medical Systems, Inc., Palo Alto, CA). For validation purposes, measured and calculated depth dose curves and dose profiles are compared for four different MLC shaped electron fields and all available energies. Furthermore, a measured two-dimensional dose distribution for patched segments consisting of three 18 MeV segments, three 12 MeV segments, and a 9 MeV segment is compared with corresponding dose calculations. Finally, measured and calculated two-dimensional dose distributions are compared for a circular segment encompassed with a C-shaped segment. RESULTS: For 15 × 34, 5 × 5, and 2 × 2 cm(2) fields differences between water phantom measurements and calculations using the beam model coupled with the macro MC dose calculation algorithm are generally within 2% of the maximal dose value or 2 mm distance to agreement (DTA) for all electron beam energies. For a more complex MLC pattern, differences between measurements and calculations are generally within 3% of the maximal dose value or 3 mm DTA for all electron beam energies. For the two-dimensional dose comparisons, the differences between calculations and measurements are generally within 2% of the maximal dose value or 2 mm DTA. CONCLUSIONS: The results of the dose comparisons suggest that the developed beam model is suitable to accurately reconstruct photon MLC shaped electron beams for a Clinac 23EX and a TrueBeam linac. Hence, in future work the beam model will be utilized to investigate the possibilities of MERT using the photon MLC to shape electron beams.