PETRA: A pencil beam trimming algorithm for analytical proton therapy dose calculations with the dynamic collimation system

BACKGROUND

The Dynamic Collimation System (DCS) has been shown to produce superior treatment plans to uncollimated pencil beam scanning (PBS) proton therapy using an in-house treatment planning system (TPS) designed for research. Clinical implementation of the DCS requires the development and benchmarking of a rigorous dose calculation algorithm that accounts for pencil beam trimming, performs monitor unit calculations to produce deliverable plans at all beam energies, and is ideally implemented with a commercially available TPS.

PURPOSE

To present an analytical Pencil bEam TRimming Algorithm (PETRA) for the DCS, with and without its range shifter, implemented in the Astroid TPS (.decimal, Sanford, Florida, USA).

MATERIALS

PETRA was derived by generalizing an existing pencil beam dose calculation model to account for the DCS-specific effects of lateral penumbra blurring due to the nickel trimmers in two different planes, integral depth dose variation due to the trimming process, and the presence and absence of the range shifter. Tuning parameters were introduced to enable agreement between PETRA and a measurement-validated Dynamic Collimation Monte Carlo (DCMC) model of the Miami Cancer Institute's IBA Proteus Plus system equipped with the DCS. Trimmer position, spot position, beam energy, and the presence or absence of a range shifter were all used as variables for the characterization of the model. The model was calibrated for pencil beam monitor unit calculations using procedures specified by International Atomic Energy Agency Technical Report Series 398 (IAEA TRS-398).

RESULTS

The integral depth dose curves (IDDs) for energies between 70 MeV and 160 MeV among all simulated trimmer combinations, with and without the ranger shifter, agreed between PETRA and DCMC at the 1%/1 mm 1-D gamma criteria for 99.99% of points. For lateral dose profiles, the median 2-D gamma pass rate for all profiles at 1.5%/1.5 mm was 99.99% at the water phantom surface, plateau, and Bragg peak depths without the range shifter and at the surface and Bragg peak depths with the range shifter. The minimum 1.5%/1.5 mm gamma pass rates for the 2-D profiles at the water phantom surface without and with the range shifter were 98.02% and 97.91%, respectively, and, at the Bragg peak, the minimum pass rates were 97.80% and 97.5%, respectively.

CONCLUSION

The PETRA model for DCS dose calculations was successfully defined and benchmarked for use in a commercially available TPS.

Integration and dosimetric validation of a dynamic collimation system for pencil beam scanning proton therapy

Objective.To integrate a Dynamic Collimation System (DCS) into a pencil beam scanning (PBS) proton therapy system and validate its dosimetric impact.Approach.Uncollimated and collimated treatment fields were developed for clinically relevant targets using an in-house treatment plan optimizer and an experimentally validated Monte Carlo model of the DCS and IBA dedicated nozzle (DN) system. The dose reduction induced by the DCS was quantified by calculating the mean dose in 10- and 30-mm two-dimensional rinds surrounding the target. A select number of plans were then used to experimentally validate the mechanical integration of the DCS and beam scanning controller system through measurements with the MatriXX-PT ionization chamber array and EBT3 film. Absolute doses were verified at the central axis at various depths using the IBA MatriXX-PT and PPC05 ionization chamber.Main results.Simulations demonstrated a maximum mean dose reduction of 12% for the 10 mm rind region and 45% for the 30 mm rind region when utilizing the DCS. Excellent agreement was observed between Monte Carlo simulations, EBT3 film, and MatriXX-PT measurements, with gamma pass rates exceeding 94.9% for all tested plans at the 3%/2 mm criterion. Absolute central axis doses showed an average verification difference of 1.4% between Monte Carlo and MatriXX-PT/PPC05 measurements.Significance.We have successfully dosimetrically validated the delivery of dynamically collimated proton therapy for clinically relevant delivery patterns and dose distributions with the DCS. Monte Carlo simulations were employed to assess dose reductions and treatment planning considerations associated with the DCS.

Design, testing and characterization of a proton central axis alignment device for the Dynamic Collimation System

OBJECTIVE

Proton therapy conformity has improved over the years by evolving from passive scattering to spot scanning delivery technologies with smaller proton beam spot sizes. Ancillary collimation devices, such the Dynamic Collimation System (DCS), further improves high dose conformity by sharpening the lateral penumbra. However, as spot sizes are reduced, collimator positional errors play a significant impact on the dose distributions and hence accurate collimator to radiation field alignment is critical. 
Approach. The purpose of this work was to develop a system to align and verify coincidence between the center of the DCS and the proton beam central axis. The Central Axis Alignment Device (CAAD) is composed of a camera and scintillating screen-based beam characterization system. Within a light-tight box, a 12.3-megapixel camera monitors a P43/Gadox scintillating screen via a 45⁰ first-surface mirror. When a collimator trimmer of the DCS is placed in the uncalibrated center of the field, the proton radiation beam continuously scans a 7x7 cm² square field across the scintillator and collimator trimmer while a 7 second exposure is acquired. From the relative positioning of the trimmer to the radiation field, the true center of the radiation field can be calculated.
Main results. The CAAD can calculate the offset between the proton beam radiation central axis and the DCS central axis within 0.054 mm accuracy and 0.075 mm reproducibility.
Significance. Using the CAAD, the DCS is now able to be aligned accurately to the proton radiation beam central axis and no longer relies on an x-ray source in the gantry head which is only validated to within 1.0 mm of the proton beam.&#xD.

Dosimetric delivery validation of dynamically collimated pencil beam scanning proton therapy

Objective. Pencil beam scanning (PBS) proton therapy target dose conformity can be improved with energy layer-specific collimation. One such collimator is the dynamic collimation system (DCS), which consists of four nickel trimmer blades that intercept the scanning beam as it approaches the lateral extent of the target. While the dosimetric benefits of the DCS have been demonstrated through computational treatment planning studies, there has yet to be experimental verification of these benefits for composite multi-energy layer fields. The objective of this work is to dosimetrically characterize and experimentally validate the delivery of dynamically collimated proton therapy with the DCS equipped to a clinical PBS system.Approach. Optimized single field, uniform dose treatment plans for 3 × 3 × 3 cm³target volumes were generated using Monte Carlo dose calculations with depths ranging from 5 to 15 cm, trimmer-to-surface distances ranging from 5 to 18.15 cm, with and without a 4 cm thick polyethylene range shifter. Treatment plans were then delivered to a water phantom using a prototype DCS and an IBA dedicated nozzle system and measured with a Zebra multilayer ionization chamber, a MatriXX PT ionization chamber array, and Gafchromic™ EBT3 film.Main results. For measurements made within the SOBPs, average 2D gamma pass rates exceeded 98.5% for the MatriXX PT and 96.5% for film at the 2%/2 mm criterion across all measured uncollimated and collimated plans, respectively. For verification of the penumbra width reduction with collimation, film agreed with Monte Carlo with differences within 0.3 mm on average compared to 0.9 mm for the MatriXX PT.Significance. We have experimentally verified the delivery of DCS-collimated fields using a clinical PBS system and commonly available dosimeters and have also identified potential weaknesses for dosimeters subject to steep dose gradients.

On the perturbation effect and LET dependence of beam quality correction factors in carbon ion beams

BACKGROUND

In a recent study, we reported beam quality correction factors, f_Q , in carbon ion beams using Monte Carlo (MC) methods for a cylindrical and a parallel-plate ionization chamber (IC). A non-negligible perturbation effect was observed; however, the magnitude of the perturbation correction due to the specific IC subcomponents was not included. Furthermore, the stopping power data presented in the International Commission on Radiation Units and Measurements (ICRU) report 73 were used, whereas the latest stopping power data have been reported in the ICRU report 90.

PURPOSE

The aim of this study was to extend our previous work by computing f_Q correction factors using the ICRU 90 stopping power data and by reporting IC-specific perturbation correction factors. Possible energy or linear energy transfer (LET) dependence of the f_Q correction factor was investigated by simulating both pristine beams and spread-out Bragg peaks (SOBPs).

METHODS

The TOol for PArticle Simulation (TOPAS)/GEANT4 MC code was used in this study. A 30 × 30 × 50 cm³ water phantom was simulated with a uniform 10 × 10 cm² parallel beam incident on the surface. A Farmer-type cylindrical IC (Exradin A12) and two parallel-plate ICs (Exradin P11 and A11) were simulated in TOPAS using the manufacturer-provided geometrical drawings. The f_Q correction factor was calculated in pristine carbon ion beams in the 150-450 MeV/u energy range at 2 cm depth and in the middle of the flat region of four SOBPs. The k_Q correction factor was calculated by simulating the f_Qo correction factor in a ⁶⁰ Co beam at 5 cm depth. The perturbation correction factors due to the presence of the individual IC subcomponents, such as the displacement effect in the air cavity, collecting electrode, chamber wall, and chamber stem, were calculated at 2 cm depth for monoenergetic beams only. Additionally, the mean dose-averaged and track-averaged LET was calculated at the depths at which the f_Q was calculated.

RESULTS

The ICRU 90 f_Q correction factors were reported. The p_dis correction factor was found to be significant for the cylindrical IC with magnitudes up to 1.70%. The individual perturbation corrections for the parallel-plate ICs were <1.0% except for the A11 p_cel correction at the lowest energy. The f_Q correction for the P11 IC exhibited an energy dependence of >1.00% and displayed differences up to 0.87% between pristine beams and SOBPs. Conversely, the f_Q for A11 and A12 displayed a minimal energy dependence of <0.50%. The energy dependence was found to manifest in the LET dependence for the P11 IC. A statistically significant LET dependence was found only for the P11 IC in pristine beams only with a magnitude of <1.10%.

CONCLUSIONS

The perturbation and k_Q correction factor should be calculated for the specific IC to be used in carbon ion beam reference dosimetry as a function of beam quality.

Investigating aperture-based approximations to model a focused dynamic collimation system for pencil beam scanning proton therapy

Purpose. The Dynamic Collimation System (DCS) is an energy layer-specific collimation device designed to reduce the lateral penumbra in pencil beam scanning proton therapy. The DCS consists of two pairs of nickel trimmers that rapidly and independently move and rotate to intercept the scanning proton beam and an integrated range shifter to treat targets less than 4 cm deep. This work examines the validity of a single aperture approximation to model the DCS, a commonly used approximation in commercial treatment planning systems, as well as higher-order aperture-based approximations for modeling DCS-collimated dose distributions.Methods. An experimentally validated TOPAS/Geant4-based Monte Carlo model of the DCS integrated with a beam model of the IBA pencil beam scanning dedicated nozzle was used to simulate DCS- and aperture-collimated 100 MeV beamlets and composite treatment plans. The DCS was represented by three different aperture approximations: a single aperture placed halfway between the upper and lower trimmer planes, two apertures located at the upper and lower trimmer planes, and four apertures, located at both the upstream and downstream faces of each pair of trimmers. Line profiles and three-dimensional regions of interest were used to evaluate the validity and limitations of the aperture approximations investigated.Results. For pencil beams without a range shifter, minimal differences were observed between the DCS and single aperture approximation. For range shifted beamlets, the single aperture approximation yielded wider penumbra widths (up to 18%) in the X-direction and sharper widths (up to 9.4%) in the Y-direction. For the example treatment plan, the root-mean-square errors (RMSEs) in an overall three-dimensional region of interest were 1.7%, 1.3%, and 1.7% for the single aperture, two aperture, and four aperture models, respectively. If the region of interest only encompasses the lateral edges outside of the target, the resulting RMSEs were 1.7%, 1.1%, and 0.5% single aperture, two aperture, and four aperture models, respectively.Conclusions. Monte Carlo simulations of the DCS demonstrated that a single aperture approximation is sufficient for modeling pristine fields at the Bragg depth while range shifted fields require a higher-order aperture approximation. For the treatment plan considered, the double aperture model performed the best overall, however, the four-aperture model most accurately modeled the lateral field edges at the expense of increased dose differences proximal to and within the target.

Development and validation of the Dynamic Collimation Monte Carlo simulation package for pencil beam scanning proton therapy

PURPOSE

The aim of this work was to develop and experimentally validate a Dynamic Collimation Monte Carlo (DCMC) simulation package specifically designed for the simulation of collimators in pencil beam scanning proton therapy (PBS-PT). The DCMC package was developed using the TOPAS Monte Carlo platform and consists of a generalized PBS source model and collimator component extensions.

METHODS

A divergent point-source model of the IBA dedicated nozzle (DN) at the Miami Cancer Institute (MCI) was created and validated against on-axis commissioning measurements taken at MCI. The beamline optics were mathematically incorporated into the source to model beamlet deflections in the X and Y directions at the respective magnet planes. Off-axis measurements taken at multiple planes in air were used to validate both the off-axis spot size and divergence of the source model. The DCS trimmers were modeled and incorporated as TOPAS geometry extensions that linearly translate and rotate about the bending magnets. To validate the collimator model, a series of integral depth dose (IDD) and lateral profile measurements were acquired at MCI and used to benchmark the DCMC performance for modeling both pristine and range shifted beamlets. The water equivalent thickness (WET) of the range shifter was determined by quantifying the shift in the depth of the 80% dose point distal to the Bragg peak between the range shifted and pristine uncollimated beams.

RESULTS

A source model of the IBA DN system was successfully commissioned against on- and off-axis IDD and lateral profile measurements performed at MCI. The divergence of the source model was matched through an optimization of the source-to-axis distance and comparison against in-air spot profiles. The DCS model was then benchmarked against collimated IDD and in-air and in-phantom lateral profile measurements. Gamma analysis was used to evaluate the agreement between measured and simulated lateral profiles and IDDs with 1%/1 mm criteria and a 1% dose threshold. For the pristine collimated beams, the average 1%/1 mm gamma pass rates across all collimator configurations investigated were 99.8% for IDDs and 97.6% and 95.2% for in-air and in-phantom lateral profiles. All range shifted collimated IDDs passed at 100% while in-air and in-phantom lateral profiles had average pass rates of 99.1% and 99.8%, respectively. The measured and simulated WET of the polyethylene range shifter was determined to be 40.9 and 41.0 mm, respectively.

CONCLUSIONS

We have developed a TOPAS-based Monte Carlo package for modeling collimators in PBS-PT. This package was then commissioned to model the IBA DN system and DCS located at MCI using both uncollimated and collimated measurements. Validation results demonstrate that the DCMC package can be used to accurately model other aspects of a DCS implementation via simulation.

Immunogenicity of the hepatitis A vaccine 20 years after infant immunization

To determine the duration of immunity provided by the Hepatitis A vaccination (HepA), we evaluated a cohort of participants in Alaska 20 years after being immunized as infants. At recruitment, participants received two doses of inactivated HepA vaccine on one of three schedules. We conducted hepatitis A antibody (anti-HAV) testing for participants at the 20-year time-point. Seventy-five of the original 183 participants (41%) were available for follow-up. The overall anti-HAV geometric mean concentration was 29.9 mIU/mL (95% CI 22.4 mIU/mL, 39.7 mIU/mL) and 50 participants (68%) remained seropositive (titer ≥ 20 mIU/mL). Using a fractional polynomial model, the predicted percent seropositive at 25 years was 55.3%, 49.8% at 30 years and 45.7% at 35 years, suggesting that the percent sero-positive could drop below 50% earlier than previously expected. Further research is necessary to understand if protection continues after seropositivity diminishes or if a HepA booster dose may become necessary.

Design of a focused collimator for proton therapy spot scanning using Monte Carlo methods

PURPOSE

When designing a collimation system for pencil beam spot scanning proton therapy, a decision must be made whether or not to rotate, or focus, the collimator to match beamlet deflection as a function of off-axis distance. If the collimator is not focused, the beamlet shape and fluence will vary as a function of off-axis distance due to partial transmission through the collimator. In this work, we quantify the magnitude of these effects and propose a focused dynamic collimation system (DCS) for use in proton therapy spot scanning.

METHODS

This study was done in silico using a model of the Miami Cancer Institute's (MCI) IBA Proteus Plus system created in Geant4-based TOPAS. The DCS utilizes rectangular nickel trimmers mounted on rotating sliders that move in synchrony with the pencil beam to provide focused collimation at the edge of the target. Using a simplified setup of the DCS, simulations were performed at various off-axis locations corresponding to beam deflection angles ranging from 0° to 2.5°. At each off-axis location, focused (trimmer rotated) and unfocused (trimmer not rotated) simulations were performed. In all simulations, a 4 cm water equivalent thickness range shifter was placed upstream of the collimator, and a voxelized water phantom that scored dose was placed downstream, each with 4 cm airgaps.

RESULTS

Increasing the beam deflection angle for an unfocused trimmer caused the collimated edge of the beamlet profile to shift 0.08-0.61 mm from the baseline 0° simulation. There was also an increase in low-dose regions on the collimated edge ranging from 14.6% to 192.4%. Lastly, the maximum dose, D max , was 0-5% higher for the unfocused simulations. With a focused trimmer design, the profile shift and dose increases were all eliminated.

CONCLUSIONS

We have shown that focusing a collimator in spot scanning proton therapy reduces dose at the collimated edge compared to conventional, unfocused collimation devices and presented a simple, mechanical design for achieving focusing for a range of source-to-collimator distances.

An overview of internet biosurveillance

Internet biosurveillance utilizes unstructured data from diverse web-based sources to provide early warning and situational awareness of public health threats. The scope of source coverage ranges from local media in the vernacular to international media in widely read languages. Internet biosurveillance is a timely modality that is available to government and public health officials, healthcare workers, and the public and private sector, serving as a real-time complementary approach to traditional indicator-based public health disease surveillance methods. Internet biosurveillance also supports the broader activity of epidemic intelligence. This overview covers the current state of the field of Internet biosurveillance, and provides a perspective on the future of the field.

Event-based biosurveillance of respiratory disease in Mexico, 2007–2009: connection to the 2009 influenza A(H1N1) pandemic?

Binary file ES_Abstracts_Final_ECDC.txt matches

Event-based biosurveillance of respiratory disease in Mexico, 2007-2009: connection to the 2009 influenza A(H1N1) pandemic?

The emergence of the 2009 pandemic influenza A(H1N1) virus in North America and its subsequent global spread highlights the public health need for early warning of infectious disease outbreaks. Event-based biosurveillance, based on local- and regional-level Internet media reports, is one approach to early warning as well as to situational awareness. This study analyses media reports in Mexico collected by the Argus biosurveillance system between 1 October 2007 and 31 May 2009. Results from Mexico are compared with the United States and Canadian media reports obtained from the HealthMap system. A significant increase in reporting frequency of respiratory disease in Mexico during the 2008-9 influenza season relative to that of 2007-8 was observed (p<0.0001). The timing of events, based on media reports, suggests that respiratory disease was prevalent in parts of Mexico, and was reported as unusual, much earlier than the microbiological identification of the pandemic virus. Such observations suggest that abnormal respiratory disease frequency and severity was occurring in Mexico throughout the winter of 2008-2009, though its connection to the emergence of the 2009 pandemic influenza A(H1N1) virus remains unclear.

Landscape of international event-based biosurveillance

Event-based biosurveillance is a scientific discipline in which diverse sources of data, many of which are available from the Internet, are characterized prospectively to provide information on infectious disease events. Biosurveillance complements traditional public health surveillance to provide both early warning of infectious disease events and situational awareness. The Global Health Security Action Group of the Global Health Security Initiative is developing a biosurveillance capability that integrates and leverages component systems from member nations. This work discusses these biosurveillance systems and identifies needed future studies.

Localization of the factor IX propeptide binding site on recombinant vitamin K dependent carboxylase using benzoylphenylalanine photoaffinity peptide inactivators

The propeptide binding/activation site on the vitamin K dependent carboxylase has been localized to a region of carboxylase between residues Arg +50 and Glu +225 by photoinactivation studies using [125I]tyrosyl-labeled benzoylphenylalanine (Bpa)-containing analogs of proFIX19, a peptide containing residues -18 to +1 of factor IX. Four proFIX19 analogs with Bpa substituents at -16, -13, -7, and -6 were synthesized. These peptides were specific photoinactivators of carboxylase and were used to label a His6-carboxylase construct produced in baculovirus-infected insect cells. Fragments of the labeled carboxylase produced by V8 protease digestion were analyzed by peptide-specific antibodies and by autoradiography. The propeptide recognition site was localized to the N-terminal one-third of the 94 kDa carboxylase. This is consistent with previous studies using a carboxylase substrate affinity label, N-(bromoacetyl)-FLEELY [Kuliopulos, A., Nelson, N.P., Yamada, M., Walsh, C.T., Furie, B., Furie, B.C., & Roth, D.A. (1994) J. Biol. Chem. 269, 21364-21370], indicating that the propeptide binding site and the FLEEL binding site are both located within the N-terminal one-third of the vitamin K dependent carboxylase.

Localization of the affinity peptide-substrate inactivator site on recombinant vitamin K-dependent carboxylase

A recombinant His6-tagged vitamin K-dependent gamma-glutamyl carboxylase has been produced in baculovirus-infected insect cells. The His6-carboxylase shares nearly identical kinetic properties with the wild-type enzyme from bovine liver microsomes. The His6-carboxylase was irreversibly inactivated by the N-bromoacetyl-FLEEL-125I-Y peptide substrate/affinity label under pseudo-first order conditions. This inactivation could be abolished by coincubation with a high affinity peptide substrate consistent with an active site-directed inactivation. The inactivated His6-carboxylase-Ac-FLEEL-125I-Y, purified under denaturing conditions by Ni-chelation chromatography followed by preparative polyacrylamide gel electrophoresis, was subjected to proteolytic digestions with either Glu-C or Lys-C endoproteinases. The resulting polypeptide fragments were probed with three regiospecific antibodies which recognized epitopes present at the extreme N terminus (residues -23 to -13), at the hydrophobic N-terminal region (residues 86-99), and at the hydrophilic C-terminal region (residues 661-673). The site of attachment to the 125I-affinity label is located within the first 218 amino acid residues of the 758-residue carboxylase. This is the first evidence for the involvement of either the putative membrane-anchoring hydrophobic region (residues 50-314) or possibly the N-terminal hydrophilic region (residues 1-50) in gamma-carboxylation of glutamate-peptide substrates.