Production of the PET isotope 11C in hadron therapy via neutron knockout

PurposePositron emitting isotopes such as ¹¹C and ¹⁰C can be used for vital dose verification in hadron therapy. These isotopes are produced when the high energy ¹²C primary beam particles undergo nuclear reactions within the patient.MethodsWe discuss a model for calculating cross sections for the production ¹¹C in ¹²C + ¹²C collisions, applicable at hadron therapy energies.ResultsGood agreement with the available cross section measurements is found for ¹²C(−1n), though more detailed, systematic measurements would be very valuable.ConclusionsNuclear structure plays a crucial role in the reactions of light nuclei, particularly when those reactions are peripheral and involve only a few nucleons. For such reactions, nuclear structure has a strong influence on the energy and angular distribution of the cross section, and is an important consideration for reliable dose verification using ¹¹C in hadron therapy.     

The determination of a dose deposited in reference medium due to (p,n) reaction occurring during proton therapy

AimThe aim of the investigation was to determine the undesirable dose coming from neutrons produced in reactions (p,n) in irradiated tissues represented by water.BackgroundProduction of neutrons in the system of beam collimators and in irradiated tissues is the undesirable phenomenon related to the application of protons in radiotherapy. It makes that proton beams are contaminated by neutrons and patients receive the undesirable neutron dose.Materials and methodsThe investigation was based on the Monte Carlo simulations (GEANT4 code). The calculations were performed for five energies of protons: 50 MeV, 55 MeV, 60 MeV, 65 MeV and 75 MeV. The neutron doses were calculated on the basis of the neutron fluence and neutron energy spectra derived from simulations and by means of the neutron fluence–dose conversion coefficients taken from the ICRP dosimetry protocol no. 74 for the antero-posterior irradiation geometry.ResultsThe obtained neutron doses are much less than the proton ones. They do not exceed 0.1%, 0.4%, 0.5%, 0.6% and 0.7% of the total dose at a given depth for the primary protons with energy of 50 MeV, 55 MeV, 60 MeV, 65 MeV and 70 MeV, respectively.ConclusionsThe neutron production takes place mainly along the central axis of the beam. The maximum neutron dose appears at about a half of the depth of the maximum proton dose (Bragg peak), i.e. in the volume of a healthy tissue. The doses of neutrons produced in the irradiated medium (water) are about two orders of magnitude less than the proton doses for the considered range of energy of protons.     

Monte Carlo calculations of radiotherapy dose in “homogeneous” anatomy

Given the substantial literature on the use of Monte Carlo (MC) simulations to verify treatment planning system (TPS) calculations of radiotherapy dose in heterogeneous regions, such as head and neck and lung, this study investigated the potential value of running MC simulations of radiotherapy treatments of nominally homogeneous pelvic anatomy. A pre-existing in-house MC job submission and analysis system, built around BEAMnrc and DOSXYZnrc, was used to evaluate the dosimetric accuracy of a sample of 12 pelvic volumetric arc therapy (VMAT) treatments, planned using the Varian Eclipse TPS, where dose was calculated with both the Analytical Anisotropic Algorithm (AAA) and the Acuros (AXB) algorithm. In-house TADA (Treatment And Dose Assessor) software was used to evaluate treatment plan complexity, in terms of the small aperture score (SAS), modulation index (MI) and a novel exposed leaf score (ELS/ELA). Results showed that the TPS generally achieved closer agreement with the MC dose distribution when treatments were planned for smaller (single-organ) targets rather than larger targets that included nodes or metastases. Analysis of these MC results with reference to the complexity metrics indicated that while AXB was useful for reducing dosimetric uncertainties associated with density heterogeneity, the residual TPS dose calculation uncertainties resulted from treatment plan complexity and TPS model simplicity. The results of this study demonstrate the value of using MC methods to recalculate and check the dose calculations provided by commercial radiotherapy TPSs, even when the treated anatomy is assumed to be comparatively homogeneous, such as in the pelvic region.     

Proton microbeam radiotherapy with scanned pencil-beams – Monte Carlo simulations

Irradiation, delivered by a synchrotron facility, using a set of highly collimated, narrow and parallel photon beams spaced by 1 mm or less, has been termed Microbeam Radiation Therapy (MRT). The tolerance of healthy tissue after MRT was found to be better than after standard broad X-ray beams, together with a more pronounced response of malignant tissue. The microbeam spacing and transverse peak-to-valley dose ratio (PVDR) are considered to be relevant biological MRT parameters. We investigated the MRT concept for proton microbeams, where we expected different depth-dose profiles and PVDR dependences, resulting in skin sparing and homogeneous dose distributions at larger beam depths, due to differences between interactions of proton and photon beams in tissue. Using the FLUKA Monte Carlo code we simulated PVDR distributions for differently spaced 0.1 mm (sigma) pencil-beams of entrance energies 60, 80, 100 and 120 MeV irradiating a cylindrical water phantom with and without a bone layer, representing human head. We calculated PVDR distributions and evaluated uniformity of target irradiation at distal beam ranges of 60–120 MeV microbeams. We also calculated PVDR distributions for a 60 MeV spread-out Bragg peak microbeam configuration. Application of optimised proton MRT in terms of spot size, pencil-beam distribution, entrance beam energy, multiport irradiation, combined with relevant radiobiological investigations, could pave the way for hypofractionation scenarios where tissue sparing at the entrance, better malignant tissue response and better dose conformity of target volume irradiation could be achieved, compared with present proton beam radiotherapy configurations.     

Configuration and validation of an analytical model predicting secondary neutron radiation in proton therapy using Monte Carlo simulations and experimental measurements

PurposeThis study focuses on the configuration and validation of an analytical model predicting leakage neutron doses in proton therapy.MethodsUsing Monte Carlo (MC) calculations, a facility-specific analytical model was built to reproduce out-of-field neutron doses while separately accounting for the contribution of intra-nuclear cascade, evaporation, epithermal and thermal neutrons. This model was first trained to reproduce in-water neutron absorbed doses and in-air neutron ambient dose equivalents, H*(10), calculated using MCNPX. Its capacity in predicting out-of-field doses at any position not involved in the training phase was also checked. The model was next expanded to enable a full 3D mapping of H*(10) inside the treatment room, tested in a clinically relevant configuration and finally consolidated with experimental measurements.ResultsFollowing the literature approach, the work first proved that it is possible to build a facility-specific analytical model that efficiently reproduces in-water neutron doses and in-air H*(10) values with a maximum difference less than 25%. In addition, the analytical model succeeded in predicting out-of-field neutron doses in the lateral and vertical direction. Testing the analytical model in clinical configurations proved the need to separate the contribution of internal and external neutrons. The impact of modulation width on stray neutrons was found to be easily adjustable while beam collimation remains a challenging issue. Finally, the model performance agreed with experimental measurements with satisfactory results considering measurement and simulation uncertainties.ConclusionAnalytical models represent a promising solution that substitutes for time-consuming MC calculations when assessing doses to healthy organs.     

All-Atom Monte Carlo Approach to Protein-Peptide Binding

We develop a procedure for exploring the free energy landscape of protein-peptide binding at atomic detail and apply it to PDZ domain-peptide interactions. The procedure involves soft constraints on receptor proteins providing limited chain flexibility, including backbone motions. Peptide chains are left fully flexible and kept in spatial proximity of the protein through periodic boundary conditions. By extensive Monte Carlo simulations, full representative conformational ensembles at temperatures where bound and unbound states coexist are obtained. To make this approach computationally feasible, we develop an effective all-atom energy function centering on hydrophobicity, hydrogen bonding, and electrostatic interactions. Our initial focus is a set of 11 PDZ domain-peptide pairs with experimentally determined complex structures. Minimum-energy conformations are found to be highly similar to the respective native structures in eight of the cases (all-atom peptide RMSDs < 6 Å). Having achieved that, we turn to a more complete characterization of the bound peptide state through a clustering scheme applied on the full ensembles of peptide structures. We find a significant diversity among bound peptide conformations for several PDZ domains, in particular involving the N terminal side of the peptide chains. Our computational model is then tested further on a set of nine PDZ domain-peptide pairs where the peptides are not originally present in the experimentally determined structures. We find a similar success rate in terms of the nativeness of minimum-energy conformations. Finally, we investigate the ability of our approach to capture variations in binding affinities for different peptide sequences. This is done in particular for a set of related sequences binding to the third PDZ domain of PSD-95 with encouraging results.