GNRA/receptor interacting modules: versatile modular units for natural and artificial RNA architectures.

Junya Ishikawa,Yuki Fujita,Y. Maeda,H. Furuta,Y. Ikawa

DOI: https://doi.org/10.1016/j.ymeth.2010.12.011

IF: 4.647

2011-06-01

Methods

Abstract:

What problem does this paper attempt to address?

Matbed_B&C: A 3-Dimensional Biologically Effective Dose Analytic Approach for the Retrospective Study of Gamma Knife Radiosurgery in a B&C Model

Ke Tang,Nan Zhang,Xiaodong Yuan,Zenghui Qian,Yang Li,Xu Feng

DOI: https://doi.org/10.1016/j.mex.2023.102320

2023-01-01

MethodsX

Abstract:The biological effect of irradiation is not solely determined by the physical dose. Gamma knife radiosurgery may be influenced by dose rate, beam-on-time, numbers of iso-centers, the gap between the individual iso-centers, and the dose‒response of various tissues. The biologically effective dose (BED) for radiosurgery considers these issues. Millions of patients treated with Models B and C provide a vast database to mine BED-related information. This research aims to develop MatBED_B&C, a 3-dimensional (3D) BED analytic approach, to generate a BED for individual voxels in the calculation matrix with related parameters extracted from Gammaplan. This approach calculates the distribution profiles of the BED in radiosurgical targets and organs at risk. A BED calculated on a voxel-by-voxel basis can be used to show the 3D morphology of the iso-BED surface and visualize the BED spatial distribution in the target. A 200 × 200 × 200 matrix can cover a greater range of the organ at risk. The BED calculated by MatBED_B&C can also be used to form BED-volume histograms to generate plan quality metrics, which will be studied in a retrospective study of gamma knife radiosurgery to guide future BED planning.•We develop MatBED_B&C to calculate the 3D BED in radiosurgical targets.•The BED of MatBED_B&C can visualize the BED spatial distribution profiles.•The BED of MatBED_B&C will generate plan quality metrics studied in a retrospective study.
Toolkit for Determination of Dose-response Relations, Validation of Radiobiological Parameters and Treatment Plan Optimization Based on Radiobiological Measures

Panayiotis Mavroidis,Athanasios Tzikas,Nikos Papanikolaou,Bengt K. Lind

DOI: https://doi.org/10.1177/153303461000900511

2010-10-01

Abstract:Accurately determined dose-response relations of the different tumors and normal tissues should be estimated and used in the clinic. The aim of this study is to demonstrate developed tools that are necessary for determining the dose-response parameters of tumors and normal tissues, for clinically verifying already published parameter sets using local patient materials and for making use of all this information in the optimization and comparison of different treatment plans and radiation techniques. One of the software modules (the Parameter Determination Module) is designed to determine the dose-response parameters of tumors and normal tissues. This is accomplished by performing a maximum likelihood fitting to calculate the best estimates and confidence intervals of the parameters used by different radiobiological models. Another module of this software (the Parameter Validation Module) concerns the validation and compatibility of external or reported dose-response parameters describing tumor control and normal tissue complications. This is accomplished by associating the expected response rates, which are calculated using different models and published parameter sets, with the clinical follow-up records of the local patient population. Finally, the last module of the software (the Radiobiological Plan Evaluation Module) is used for estimating and optimizing the effectiveness a treatment plan in terms of complication-free tumor control, P + . The use of the Parameter Determination Module is demonstrated by deriving the dose-response relation of proximal esophagus from head & neck cancer radiotherapy. The application of the Parameter Validation Module is illustrated by verifying the clinical compatibility of those dose-response parameters with the examined treatment methodologies. The Radiobiological Plan Evaluation Module is demonstrated by evaluating and optimizing the effectiveness of head & neck cancer treatment plans. The results of the radiobiological evaluation are compared against dosimetric criteria. The presented toolkit appears to be very convenient and efficient for clinical implementation of radiobiological modeling. It can also be used for the development of a clinical data and health information database for assisting the performance of epidemiological studies and the collaboration between different institutions within research and clinical frameworks.

oncology
SU‐FF‐T‐05: A Monte Carlo Based Dose Calculation and Evaluation Toolkit for Electronic Brachytherapy: Feasibility of IMBT

B. Guo,C. Cheng,B. Director,T. Rusch,C. Esquivel,S. Stathakis,C. Shi,N. Papanikolaou

DOI: https://doi.org/10.1118/1.3181476

2009-01-01

Abstract:Purpose: To develop a Monte Carlo based dose calculation and evaluation toolkit for electronic brachytherapy sources, capable of simulating Intensity Modulated Brachytherapy (IMBT). Material and methods: This toolkit used Monte Carlo code EGS4 to calculate the dose distribution for a treatment in a realistic virtual human phantom converted from the patient's CTimages and contours. An in‐house Matlab program was developed to analyze the dose distribution and generate DVHs and isodoses. Results: The system was benchmarked by comparing the calculated radial dose functions (in water) with experimental data published. Difference was within 3%. A typical intracavitary accelerated partial breast irradiation (APBI) treatment plan using Xoft Axxent electronic brachytherapy source was simulated using this toolkit. DVHs and isodoses revealed that the dose to the ribs were high in this plan due to proximity of balloon to the chest wall and the high absorption coefficient of bone to low energy X‐rays. A simple IMBT plan using partial block could conform the isodose distribution to the target, reduce the dose to ribs and chest wall without compromising the dose homogeneity to the target or increasing the dose to other critical structures. Conclusion: A Monte Carlo based dose calculation and evaluation system was developed for electronic brachytherapy sources. Benchmarked through published data, this system is capable of producing reliable and detailed dose distributions for both isotropic and intensity modulated sources. The feasibility of intensity modulation in improving the plan quality was proved.
Dose-response modeling and treatment plan assessment with a python software toolkit

Athanasios Tzikas,Eleftherios Lavdas,Dimitrios Kechagias,Panayiotis Mavroidis

DOI: https://doi.org/10.1016/j.meddos.2024.04.005

IF: 1.531

2024-05-24

Medical Dosimetry

Abstract:This software assistant aims at calculating the dose-response relations of tumors and normal tissues, or clinically assessing already determined values by other researchers. It can also indicate the optimal dose prescription by optimizing the expected treatment outcome. The software is developed solely in python programming language, and it employs PSFL license for its Graphical User Interface (GUI), NUMPY, MATPLOTLIB, and SCIPY libraries. It comprises of two components. The first is the Dose-response relations derivation component, which takes as input the dose volume histograms (DVHs) of patients and their recorded responses regarding a given clinical endpoint to determine the parameters of different tumor control probability (TCP) or normal tissue complication probability (NTCP) models. The second is the Treatment Plan Assessment component, which uses the DVHs of a plan and the dose-response parameters values of the involved tumors and organs at risk (OARs) to calculate their expected responses. Additionally, the overall probabilities of benefit ( P B ), injury ( P I ) and complication-free tumor control ( P + ) are calculated. The software calculates rapidly the corresponding generalized equivalent uniform doses (gEUD) and biologically effective uniform doses ( D ̅ ̅ ) for the Lyman-Kutcher-Burman (LKB), parallel volume (PV) and relative seriality (RS) models respectively, determining the model parameters. In the Dose-Response Relations Derivation component, the software plots the dose-response curves of the irradiated organ with the relevant confidence internals along with the data of the patients with and without toxicity. It also calculates the odds ratio (OR) and the area under the curve (AUC) of different dose metrics or model parameter values against the individual patient outcomes to determine their discrimination capacity. It also performs a goodness-of-fit evaluation of any model parameter set. The user has the option of viewing plots like Scatter, 3D surfaces, and Bootstrap plots. In the Treatment Plan Assessment part, the software calculates the TCP and NTCP values of the involved tumors and OARs, respectively. Furthermore, it plots the dose-response curves of the TCPs, NTCPs, P B , P I , and P + for a range of prescription doses for different treatment plans. The presented software is ideal for efficiently conducting studies of radiobiological modeling. Furthermore, it is ideal for performing treatment plan assessment, comparison, and optimization studies.

oncology,radiology, nuclear medicine & medical imaging
Clinical treatment planning for kilovoltage radiotherapy using EGSnrc and Python

Mihails Nikandrovs,Brendan McClean,Laura Shields,Patrick McCavana,Luis León Vintró

DOI: https://doi.org/10.1002/acm2.13832

2022-11-30

Journal of applied clinical medical physics [electronic resource] / American College of Medical Physics

Abstract:Kilovoltage radiotherapy dose calculations are generally performed with manual point dose calculations based on water dosimetry. Tissue heterogeneities, irregular surfaces, and introduction of lead cutouts for treatment are either not taken into account or crudely approximated in manual calculations. Full Monte Carlo (MC) simulations can account for these limitations but require a validated treatment unit model, accurately segmented patient tissues and a treatment planning interface (TPI) to facilitate the simulation setup and result analysis. EGSnrc was used in this work to create a model of Xstrahl kilovoltage unit extending the range of energies, applicators, and validation parameters previously published. The novel functionality of the Python‐based framework developed in this work allowed beam modification using custom lead cutouts and shields, commonly present in kilovoltage treatments, as well as absolute dose normalization using the output of the unit. 3D user‐friendly planning interface of the developed framework facilitated non‐co‐planar beam setups for CT phantom MC simulations in DOSXYZnrc. The MC models of 49 clinical beams showed good agreement with measured and reference data, to within 2% for percentage depth dose curves, 4% for beam profiles at various depths, 2% for backscatter factors, 0.5 mm of absorber material for half‐value layers, and 3% for output factors. End‐to‐end testing of the framework using custom lead cutouts resulted in good agreement to within 3% of absolute dose distribution between simulations and EBT3 GafChromic film measurements. Gamma analysis demonstrated poor agreement at the field edges which was attributed to the limitations of simulating smooth cutout shapes. Dose simulated in a heterogeneous phantom agreed to within 7% with measured values converted using the ratio of mass energy absorption coefficients of appropriate tissues and air.
A Monte-Carlo-based and GPU-accelerated 4D-dose calculator for a pencil-beam scanning proton therapy system

Mark D. Pepin,Erik Tryggestad,H. Wan Chan Tseung,Jedediah E. Johnson,Michael G. Herman,Chris Beltran

DOI: https://doi.org/10.1002/mp.13182

2018-09-13

Abstract:Purpose: The presence of respiratory motion during radiation treatment leads to degradation of the expected dose distribution, both for target coverage and healthy-tissue sparing, particularly for techniques like pencil-beam scanning proton therapy which have dynamic delivery systems. While tools exist to estimate this degraded four-dimensional (4D) dose, they typically have one or more deficiencies such as ... Methods: To quickly compute the 4D-dose, the three main tasks of the calculator were run on graphics processing units (GPUs). These tasks were: simulating the delivery of the plan using measured delivery parameters to distribute the plan amongst 4DCT phases characterizing the patient breathing, using an in-house Monte Carlo simulation (MC) dose calculator to determine the dose delivered to each breathing phase, and accumulating the doses from the various breathing phases onto a single phase for evaluation. The accumulation was performed by individually transferring the energy and mass of dose-grid subvoxels, a technique models the transfer of dose in a more physically realistic manner. The calculator was run ... Results: 4D doses were successfully computed for the three test cases with computation times ranging from 4-6 min on a server with eight NVIDIA Titan X graphics cards; the most time-consuming component was the MC dose engine. The subvoxel-based dose-accumulation technique produced stable 4D-dose distributions at subvoxel scales of 0.5-1.0 mm without impairing the total computation time. The uncertainties in the beam-delivery simulation ... Conclusions: A MC-based and GPU-accelerated 4D-dose calculator was developed to estimate the effects of respiratory motion on pencil-beam scanning proton therapy treatments. The calculator can currently be used ...

Medical Physics
RapidBrachyMCTPS 2.0: A Comprehensive and Flexible Monte Carlo-Based Treatment Planning System for Brachytherapy Applications

Harry Glickman,Majd Antaki,Christopher Deufel,Shirin A. Enger

DOI: https://doi.org/10.48550/arXiv.2007.02902

2020-07-07

Abstract:We have previously described RapidBrachyMCTPS, a brachytherapy treatment planning toolkit consisting of a graphical user interface (GUI) and a Geant4-based Monte Carlo (MC) dose calculation engine. This work describes the tools that have recently been added to RapidBrachyMCTPS, such that it now serves as the first stand-alone application for MC-based brachytherapy treatment planning. Notable changes include updated applicator import and positioning, three-plane contouring tools, and updated dose optimization algorithms that, in addition to optimizing dwell position and dwell time, also optimize the rotating shield angles in intensity modulated brachytherapy. The main modules of RapidBrachyMCTPS were validated including DICOM import, applicator import and positioning, contouring, material assignment, source specification, catheter reconstruction, EGSphant generation, interface with the MC code, and dose optimization and analysis tools. Two patient cases were simulated to demonstrate these principles, illustrating the control and flexibility offered by RapidBrachyMCTPS for all steps of the treatment planning pathway. RapidBrachyMCTPS is now a stand-alone application for brachytherapy treatment planning, and offers a user-friendly interface to access powerful MC calculations. It can be used to validate dose distributions from clinical treatment planning systems or model-based dose calculation algorithms, and is also well suited to testing novel combinations of radiation sources and applicators, especially those shielded with high-Z materials.

Medical Physics
Calculation of Organ Dose Distribution (in-field and Out-of-field) in Breast Cancer Radiotherapy on RANDO Phantom Using GEANT4 Application for Tomographic Emission (Gate) Monte Carlo Simulation

Marziyeh Behmadi,Mohammad Taghi Bahreyni Toossi,Shahrokh Nasseri,Mohammad Ehsan Ravari,Mahdi Momennezhad,Hamid Gholamhosseinian,Mohammad Mohammadi,Sibusiso Mdletshe

DOI: https://doi.org/10.4103/jmss.jmss_25_23

2024-07-10

Abstract:Introduction: Organ dose distribution calculation in radiotherapy and knowledge about its side effects in cancer etiology is the most concern for medical physicists. Calculation of organ dose distribution for breast cancer treatment plans with Monte Carlo (MC) simulation is the main goal of this study. Materials and methods: Elekta Precise linear accelerator (LINAC) photon mode was simulated and verified using the GEANT4 application for tomographic emission. Eight different radiotherapy treatment plans on RANDO's phantom left breast were produced with the ISOgray treatment planning system (TPS). The simulated plans verified photon dose distribution in clinical tumor volume (CTV) with TPS dose volume histogram (DVH) and gamma index tools. To verify photon dose distribution in out-of-field organs, the point dose measurement results were compared with the same point doses in the MC simulation. Eventually, the DVHs for out-of-field organs that were extracted from the TPS and MC simulation were compared. Results: Based on the implementation of gamma index tools with 2%/2 mm criteria, the simulated LINAC output demonstrated high agreement with the experimental measurements. Plan simulation for in-field and out-of-field organs had an acceptable agreement with TPS and experimental measurement, respectively. There was a difference between DVHs extracted from the TPS and MC simulation for out-of-field organs in low-dose parts. This difference is due to the inability of the TPS to calculate dose distribution in out-of-field organs. Conclusion and discussion: Based on the results, it was concluded that the treatment plans with the MC simulation have a high accuracy for the calculation of out-of-field dose distribution and could play a significant role in evaluating the important role of dose distribution for second primary cancer estimation.
SU-E-T-513: Clinical Implementation of a GPU Accelerated Pencil Beam Dose Calculation Algorithm.

H. Chen,Z. Xiao,D. Du,D. Yang,S. Mutic,H. Li

DOI: https://doi.org/10.1118/1.4735602

2012-01-01

Abstract:PURPOSE:This work reports a clinical implementation of a GPU accelerated pencil beam dose calculation algorithm (GPU-PB).METHODS:Model parameters were determined using in-house scripts written in MATLAB. Dose distributions in a water phantom were calculated using a Pinnacle TPS for various open field sizes. Lateral profiles at 2-mm incremental depths were used to calculate PB kernel parameters. Weighted sum of squares were used for least-squares parameter fitting utilizing a Levenberg-marquardt algorithm. Weightings were adjusted based on goodness of fitting and in accordance with field size to suppress deviations due to horn effects. The scale factor was fitted iteratively. The calculated doses for two patient cases were analyzed with a 3D gamma method (3%/3mm).RESULTS:Excellent agreement between Pinnacle calculation and GPU-PB calculation was achieved regarding PDD, profiles and output factor. For a head-neck 10-field step-and-shoot IMRT case, gamma passing rate was over 99% and maximum absolute dose value is 75.6 Gy vs. 73.9 Gy, differing by 2.3%. Similar results were obtained for a SBRT lung case. Gamma passing rate is a function of PB kernel cut-off distance and beamlet resolution. It appears that if the highest accuracy is desired, a resolution of 10 mm or better in the direction parallel to MLC travel and a cutoff distance of 10 cm or better should be used. The calculation time increases with both cut-off distance and beamlet resolution. For example, GPU calculation time is 1.36 seconds for 5 cm cut-off distance, and increases to 4.84 s for 15 cm cut-off distance.CONCLUSIONS:A GPU accelerated PB dose calculation algorithm has been implemented using clinical measurement data. Excellent agreement with Pinnacle TPS has been achieved. Beamlet size and PB cut-off distance should be chosen according to desired dose calculation accuracy and speed.
SU-GG-T-38: Dosimetry and Inverse Treatment Planning for 3D Intensity Modulated Brachytherapy

B. Guo,C. Shi,C. Cheng,C. Esquivel,T. Eng,N. Papanikolaou

DOI: https://doi.org/10.1118/1.3468424

IF: 4.506

2010-01-01

Medical Physics

Abstract:Purpose: To propose a dosimetry algorithm for three-dimensional (3D) intensity modulated brachytherapy (IMBT) and to develop an inverse treatment planning method applying the dosimetry algorithm. Method and Materials: A Matlab based prototype 3D IMBT treatment planning system was developed. The system consisted of three main components: (1) a comprehensive source commissioning method for intensity modulated sources based on Monte Carlo (EGSnrc) simulations; (2) a “modified TG43” (mTG43) dose calculation algorithm for IMBT dosimetry; (3) an inverse IMBT treatment planning method based on Dose Volume Histogram (DVH) or Dose Surface Histogram (DSH) constraints and simulated annealing optimization algorithm. The system was applied for planning of an intracavitary accelerated partial breast irradiation (APBI) case treated with Xoft Axxent electronic brachytherapy. Plan quality, planning and delivery time of the IMBT plan were compared with the original plan used for the patient's treatment. Results: For the patient studied, IMBT plan showed better plan quality compared with the original plan. With similar coverage to the target, high dose region V200 was decrease by 16.1%. Maximum doses to skin and ribs were reduced by 56 cGy and 104 cGy in one fraction respectively. Mean dose to ipsilateral and contralateral breasts and lungs were also slightly reduced by IMBT. Conclusion: Application of three-dimensional intensity modulation in brachytherapy treatment planning is both feasible and promising. 3D IMBT improves the quality of APBI brachytherapy treatment plan, increasing dose uniformity in target and reducing the dose to critical structures.
Interactive Treatment Planning in High Dose-Rate Brachytherapy for Gynecological Cancer

Huan Liu,Chang M Ma,Xun Jia,Chenyang Shen,Peter Klages,Kevin Albuquerque

DOI: https://doi.org/10.48550/arXiv.2109.05081

2021-09-11

Abstract:High dose-rate brachytherapy (HDRBT) is widely used for gynecological cancer treatment. Although commercial treatment planning systems (TPSs) have inverse optimization modules, it takes several iterations to adjust planning objectives to achieve a satisfactory plan. Interactive plan-modification modules enable modifying the plan and visualizing results in real time, but they update plans based on simple geometrical or heuristic algorithms, which cannot ensure resulting plan optimality. This project develops an interactive plan optimization module for HDRBT of gynecological cancer. By efficiently solving an optimization problem in real time, it allows a user to visualize a plan and interactively modify it to improve quality. We formulated an optimization problem with an objective function containing a weighted sum of doses to normal organs subject to user-specified target coverage. A user interface was developed that allows a user to adjust organ weights using scroll bars. With a simple mouse click, the optimization problem is solved in seconds with a highly efficient alternating-direction method of multipliers and a warm start optimization strategy. Resulting clinically relevant D2cc of organs are displayed immediately. This allows a user to intuitively adjust plans with satisfactory quality. We tested the effectiveness of our development in cervix cancer cases treated with a tandem-and-ovoid applicator. It took a maximum of 3 seconds to solve the optimization problem in each instance. With interactive optimization capability, a satisfactory plan can be obtained in <1 min. In our clinic, although the time for plan adjustment was typically <10min with simple interactive plan modification tools in TPS, the resulting plans do not ensure optimality. Our plans achieved on average 5% lower D2cc than clinical plans, while maintaining target coverage.

Medical Physics
An interface tool to parametrize treatment plans for the TrueBeam radiotherapy system into TOPAS parameter control files for Monte Carlo simulation

Ramon Ortiz,Daren Sawkey,Bruce Faddegon,Naoki D-Kondo,José Ramos-Méndez

DOI: https://doi.org/10.1016/j.ejmp.2024.104485

Abstract:Purpose: The Monte Carlo (MC) method, the gold standard method for radiotherapy dose calculations, is underused in clinical research applications mainly due to computational speed limitations. Another reason is the time-consuming and error prone conversion of treatment plan specifications into MC parameters. To address this issue, we developed an interface tool that creates a set of TOPAS parameter control files (PCF) from information exported from a clinical treatment planning system (TPS) for plans delivered by the TrueBeam radiotherapy system. Methods: The interface allows the user to input DICOM-RT files, exported from a TPS and containing the plan parameters, and choose different multileaf-collimator models, variance reduction technique parameters, scoring quantities and simulation output formats. Radiation sources are precomputed phase space files obtained from Varian. Based on this information, ready-to-run TOPAS PCF that incorporate the position and angular rotation of the TrueBeam dynamic collimation devices, gantry, couch, and patient according to treatment plan specifications are created. Results: Dose distributions computed using these PCF were compared against predictions from commercial TPS for different clinical treatment plans and techniques (3D-CRT, IMRT step-and-shoot and VMAT) to evaluate the performance of the interface. The agreement between dose distributions from TOPAS and TPS (>98 % pass ratio in the gamma test) confirmed the correct parametrization of treatment plan specifications into MC PCF. Conclusions: This interface tool is expected to widen the use of MC methods in the clinical medical physics field by facilitating the straightforward transfer of treatment plan parameters from commercial TPS into MC PCF.
Conservation of Pyramidal Tract in Radiosurgery for Brain Metastases of Lung Adenocarcinoma: Three-dimensional Analysis of Biologically Effective Dose

Ke Tang,Nan Zhang,Xiaodong Yuan,Zenghui Qian,Yang Li,Xu Feng

DOI: https://doi.org/10.1016/j.radonc.2022.109451

IF: 6.901

2022-12-30

Radiotherapy and Oncology

Abstract:Background Gamma knife surgery (GKS) for brain metastases (BMs) adjacent to the pyramidal tract (PT) is still a challenge to conduct. PT visualization and biologically effective dose (BED) calculation on a voxel-by-voxel basis may provide data to establish clinically safe values. We aimed to assess the relationship of parameters extracted from the BED-volume histogram with outcomes of PT after GKS-treating target (adjacent BM of lung adenocarcinoma). Methods We formed BED-volume histograms for 672 BMs in a retrospective cohort, using 3-dimensional (3D) coordinate values of PT, target, and each iso-centre to calculate the 3D BED distribution in a 200×200×200 matrix. PT conservation failure (PTCF) was judged clinically and radiologically and classified as lesion progression and radionecrosis. Cox proportional hazards models were used to analyse 3D BED parameters. Internal validation of models was performed by bootstrapping. Results There were 116 (17.3%) subjects with PTCF in the cohort, of which 74 (11.0%) and 42 (6.3%) were caused by lesion progression and radionecrosis, respectively. Multivariate analysis showed that D Lesion_min BED and D Lesion_90% BED significantly predicted lesion progression (P< .001). D PT_Max BED and V PT_ BED40 significantly predicted radionecrosis (P< .001). The model predicting PTCF showed fair discrimination and calibration of D Lesion_min BED + D Lesion_90% BED and D PT_Max BED + V PT_ BED40 . Conclusions The conservation of PT in GKS for BMs of lung adenocarcinoma depends on the combination of PT-tolerated BED and target effective control BED. Therefore, a BED-volume histogram with a 3D BED algorithm is proposed to assess plan quality.

oncology,radiology, nuclear medicine & medical imaging
Fast full patient and radioisotope Monte Carlo simulations of targeted radionuclide therapy: introducing egs_mird

Martin P. Martinov,Chidera Opera,Rowan M. Thomson,Ting-Yim Lee

DOI: https://doi.org/10.48550/arXiv.2204.01838

IF: 4.506

2022-04-04

Medical Physics

Abstract:BACKGROUND: Targeted Radionuclide Therapy (TRT) is a fast-growing field garnering much interest, however, dosimetry calculation techniques remain relatively simple. PURPOSE: To introduce egs_mird, a new Monte Carlo application built in EGSnrc which allows users to model full patient tissue and density (using clinical CT images) and radionuclide distribution (using clinical PET images) for fast and detailed TRT dose calculation. METHODS: The novel application egs_mird is introduced, including its structure and variation reduction approaches. A new egs++ source class egs_internal_source and a modified version of egs_radionuclide_source are described. The new code is compared to other MC calculations of S-value kernels, along with self-validation using the electron Fano test. Full patient prostate 177Lu TRT cancer treatment simulations are performed using a single set of patient DICOM CT and [18F]-DCFPyL PET data. RESULTS: Good agreement is found between S-value kernels calculated using egs_mird and those found in the literature. The Fano test is satisfied to 0.1%. Patient prostate, rectum, bone marrow, and bladder dose volume histogram results did not vary significantly when using the track-length estimator and not modelling electron transport, modelling bone marrow explicitly (rather than using generic tissue compositions), and reducing activity to voxels containing partial or full calcifications to half or none, respectively. Simulations using the track-length estimator can be completed in under 15 minutes. CONCLUSION: This work shows egs mird to be a reliable MC code for computing TRT doses as realistically as the patient CT and PET data allow, supporting the use of egs mird for dose calculations in TRT.
Regulation of tyrosine hydroxylase activity in pheochromocytoma PC-12 cells by bradykinin.

H. Houchi,J. Masserano,J. Bowyer,N. Weiner

Molecular Pharmacology

Abstract:Tyrosine hydroxylase is activated and phosphorylated following treatment of PC-12 cells with bradykinin. In order to determine the mechanisms by which this occurs, we have evaluated the second messenger systems that may be responsible for this activation and phosphorylation. Inositol phosphates appear to play an important role in the activation and phosphorylation of tyrosine hydroxylase because bradykinin treatment significantly increased the formation of [3H]inositol phosphates and the concentration of intracellular free calcium ([Ca2+]i) in PC-12 cells. The uptake of extracellular 45Ca2+ into PC-12 cells at 1 min was significantly increased (107%) by bradykinin treatment and this increase was blocked by La3+, an inorganic calcium channel inhibitor, but not by nifedipine, an inhibitor of voltage-dependent calcium channels. The activation of tyrosine hydroxylase in PC-12 cells following bradykinin treatment was partially inhibited by La3+. Additivity experiments were performed to evaluate whether the activation and phosphorylation of tyrosine hydroxylase in PC-12 cells following treatment with bradykinin (10 microM) was similar to the activation and phosphorylation of tyrosine hydroxylase in PC-12 cells following treatment with dibutyryl cAMP (2 mM), 4 beta-phorbol-12 beta-myristate-13 alpha-acetate (PMA) (2 microM), and high K+ (56 mM). The combination of bradykinin and PMA produced additive effects, indicating that the activation of tyrosine hydroxylase by treatment with these two compounds was through different mechanisms. Furthermore, exposure of PC-12 cells to bradykinin did not increase intracellular cAMP levels. The combination of bradykinin and PMA treatments produced only partial additivity in tyrosine hydroxylase activity and phosphorylation. No additivity was produced with bradykinin and high K-treatment. Phosphopeptide analysis was performed on tyrosine hydroxylase obtained from PC-12 cells treated with bradykinin. Bradykinin treatment produced a significant incorporation of [32P]-phosphate into two phosphopeptides of tryptically digested tyrosine hydroxylase. One of these peptides corresponds to a peptide obtained by trypsinization of purified tyrosine hydroxylase that is phosphorylated by purified calcium/calmodulin-dependent protein kinase. The other 32P-tyrosine hydroxylase-peptide obtained from PC-12 cells treated with bradykinin corresponds to the phosphorylation site obtained during PMA stimulation of PC-12 cells. These results indicate that bradykinin treatment increases intracellular inositol phosphates, calcium, and possibly diacylglycerol levels in PC-12 cells. These effects could then increase calcium/calmodulin-dependent protein kinase activity and possibly calcium/phospholipid-dependent protein (protein kinase C) activity, resulting in increased phosphorylation and activity of tyrosine hydroxylase.
Efficiently processing (p,ε)-approximate join aggregation on massive data

Xixian Han,Jianzhong Li,Hong Gao

DOI: https://doi.org/10.1016/j.ins.2014.03.091

IF: 8.1

2014-09-10

Information Sciences

Abstract:
Patient-specific dosimetric endpoints based treatment plan quality control in radiotherapy

Ting Song,David Staub,Mingli Chen,Weiguo Lu,Zhen Tian,Xun Jia,Yongbao Li,Linghong Zhou,Steve B Jiang,Xuejun Gu

DOI: https://doi.org/10.1088/0031-9155/60/21/8213

IF: 3.5

2015-10-08

Physics in Medicine and Biology

Abstract:In intensity modulated radiotherapy (IMRT), the optimal plan for each patient is specific due to unique patient anatomy. To achieve such a plan, patient-specific dosimetric goals reflecting each patient's unique anatomy should be defined and adopted in the treatment planning procedure for plan quality control. This study is to develop such a personalized treatment plan quality control tool by predicting patient-specific dosimetric endpoints (DEs). The incorporation of patient specific DEs is realized by a multi-OAR geometry-dosimetry model, capable of predicting optimal DEs based on the individual patient's geometry. The overall quality of a treatment plan is then judged with a numerical treatment plan quality indicator and characterized as optimal or suboptimal. Taking advantage of clinically available prostate volumetric modulated arc therapy (VMAT) treatment plans, we built and evaluated our proposed plan quality control tool. Using our developed tool, six of twenty evaluated plans were identified as sub-optimal plans. After plan re-optimization, these suboptimal plans achieved better OAR dose sparing without sacrificing the PTV coverage, and the dosimetric endpoints of the re-optimized plans agreed well with the model predicted values, which validate the predictability of the proposed tool. In conclusion, the developed tool is able to accurately predict optimally achievable DEs of multiple OARs, identify suboptimal plans, and guide plan optimization. It is a useful tool for achieving patient-specific treatment plan quality control.

engineering, biomedical,radiology, nuclear medicine & medical imaging
Fred: a GPU-accelerated fast-Monte Carlo code for rapid treatment plan recalculation in ion beam therapy

A Schiavi,M Senzacqua,S Pioli,A Mairani,G Magro,S Molinelli,M Ciocca,G Battistoni,V Patera

DOI: https://doi.org/10.1088/1361-6560/aa8134

IF: 3.5

2017-09-05

Physics in Medicine and Biology

Abstract:Ion beam therapy is a rapidly growing technique for tumor radiation therapy. Ions allow for a high dose deposition in the tumor region, while sparing the surrounding healthy tissue. For this reason, the highest possible accuracy in the calculation of dose and its spatial distribution is required in treatment planning. On one hand, commonly used treatment planning software solutions adopt a simplified beam-body interaction model by remapping pre-calculated dose distributions into a 3D water-equivalent representation of the patient morphology. On the other hand, Monte Carlo (MC) simulations, which explicitly take into account all the details in the interaction of particles with human tissues, are considered to be the most reliable tool to address the complexity of mixed field irradiation in a heterogeneous environment. However, full MC calculations are not routinely used in clinical practice because they typically demand substantial computational resources. Therefore MC simulations are usually only used to check treatment plans for a restricted number of difficult cases. The advent of general-purpose programming GPU cards prompted the development of trimmed-down MC-based dose engines which can significantly reduce the time needed to recalculate a treatment plan with respect to standard MC codes in CPU hardware. In this work, we report on the development of fred, a new MC simulation platform for treatment planning in ion beam therapy. The code can transport particles through a 3D voxel grid using a class II MC algorithm. Both primary and secondary particles are tracked and their energy deposition is scored along the trajectory. Effective models for particle-medium interaction have been implemented, balancing accuracy in dose deposition with computational cost. Currently, the most refined module is the transport of proton beams in water: single pencil beam dose-depth distributions obtained with fred agree with those produced by standard MC codes within 1-2% of the Bragg peak in the therapeutic energy range. A comparison with measurements taken at the CNAO treatment center shows that the lateral dose tails are reproduced within 2% in the field size factor test up to 20 cm. The tracing kernel can run on GPU hardware, achieving 10 million primary [Formula: see text] on a single card. This performance allows one to recalculate a proton treatment plan at 1% of the total particles in just a few minutes.

engineering, biomedical,radiology, nuclear medicine & medical imaging
A GPU-based Finite-Size Pencil Beam Algorithm with 3D-Density Correction for Radiotherapy Dose Calculation.

X. Gu,U. Jelen,J. Li,X. Jia,S. Jiang

DOI: https://doi.org/10.1088/0031-9155/56/11/010

2011-01-01

Abstract:Targeting at the development of an accurate and efficient dose calculation engine for online adaptive radiotherapy, we have implemented a finite-size pencil beam (FSPB) algorithm with a 3D-density correction method on graphics processing unit (GPU). This new GPU-based dose engine is built on our previously published ultrafast FSPB computational framework (Gu et al 2009 Phys. Med. Biol. 54 6287–97). Dosimetric evaluations against Monte Carlo dose calculations are conducted on ten IMRT treatment plans (five head-and-neck cases and five lung cases). For all cases, there is improvement with the 3D-density correction over the conventional FSPB algorithm and for most cases the improvement is significant. Regarding the efficiency, because of the appropriate arrangement of memory access and the usage of GPU intrinsic functions, the dose calculation for an IMRT plan can be accomplished well within 1 s (except for one case) with this new GPU-based FSPB algorithm. Compared to the previous GPU-based FSPB algorithm without 3D-density correction, this new algorithm, though slightly sacrificing the computational efficiency (∼5–15% lower), has significantly improved the dose calculation accuracy, making it more suitable for online IMRT replanning.
Further development of spinal cord retreatment dose estimation: including radiotherapy with protons and light ions

Joshua W. Moore,Thomas E. Woolley,John W. Hopewell,Bleddyn Jones

DOI: https://doi.org/10.1101/2020.12.20.20248585

2020-12-22

Abstract:Abstract A new graphical user interface (GUI) was developed to aid in the assessment of changes in the radiation tolerance of spinal cord/similar central nervous system tissues with time between two treatment courses. The GUI allows any combination of photons, protons (or ions) to be used as the initial, or retreatment, course. Allowances for clinical circumstances, of reduced tolerance, can also be made. The radiobiological model was published previously and has been incorporated with additional checks and safety features, to be as conservative as possible. The proton option includes use of a fixed RBE of 1.1 (set as the default), or a variable RBE, the latter depending on the proton linear energy transfer (LET) for organs at risk. This second LET-based approach can also be used for ions, by changing the LET parameters. GUI screenshots are used to show the input and output parameters for clinical situations used in worked examples from previous publications, where the proton and ion treatments required additional ‘longhand’ calculations. The results from the GUI are in agreement with previously published calculations, but the results are now rapidly available without tedious and error-prone manual computations. The software outputs provide a maximum dose limit boundary, which should not be exceeded. Clinicians may also choose a lower number of treatment fractions, whilst using the same dose per fraction (or conversely a lower dose per fraction but with the same number of fractions) in order to achieve the intended clinical benefit. The new GUI will allow rational estimations of time related radiation tolerance changes in the spinal cord and similar central nervous tissues (optic chiasm, brainstem), which can be used to guide the choice of retreatment dose fractionation schedules.

GNRA/receptor interacting modules: versatile modular units for natural and artificial RNA architectures.

Matbed_B&C: A 3-Dimensional Biologically Effective Dose Analytic Approach for the Retrospective Study of Gamma Knife Radiosurgery in a B&C Model

Toolkit for Determination of Dose-response Relations, Validation of Radiobiological Parameters and Treatment Plan Optimization Based on Radiobiological Measures

SU‐FF‐T‐05: A Monte Carlo Based Dose Calculation and Evaluation Toolkit for Electronic Brachytherapy: Feasibility of IMBT

Dose-response modeling and treatment plan assessment with a python software toolkit

Clinical treatment planning for kilovoltage radiotherapy using EGSnrc and Python

A Monte-Carlo-based and GPU-accelerated 4D-dose calculator for a pencil-beam scanning proton therapy system

RapidBrachyMCTPS 2.0: A Comprehensive and Flexible Monte Carlo-Based Treatment Planning System for Brachytherapy Applications

Calculation of Organ Dose Distribution (in-field and Out-of-field) in Breast Cancer Radiotherapy on RANDO Phantom Using GEANT4 Application for Tomographic Emission (Gate) Monte Carlo Simulation

SU-E-T-513: Clinical Implementation of a GPU Accelerated Pencil Beam Dose Calculation Algorithm.

SU-GG-T-38: Dosimetry and Inverse Treatment Planning for 3D Intensity Modulated Brachytherapy

Interactive Treatment Planning in High Dose-Rate Brachytherapy for Gynecological Cancer

An interface tool to parametrize treatment plans for the TrueBeam radiotherapy system into TOPAS parameter control files for Monte Carlo simulation

Conservation of Pyramidal Tract in Radiosurgery for Brain Metastases of Lung Adenocarcinoma: Three-dimensional Analysis of Biologically Effective Dose

Fast full patient and radioisotope Monte Carlo simulations of targeted radionuclide therapy: introducing egs_mird

Regulation of tyrosine hydroxylase activity in pheochromocytoma PC-12 cells by bradykinin.

Efficiently processing (p,ε)-approximate join aggregation on massive data

Patient-specific dosimetric endpoints based treatment plan quality control in radiotherapy

Fred: a GPU-accelerated fast-Monte Carlo code for rapid treatment plan recalculation in ion beam therapy

A GPU-based Finite-Size Pencil Beam Algorithm with 3D-Density Correction for Radiotherapy Dose Calculation.

Further development of spinal cord retreatment dose estimation: including radiotherapy with protons and light ions