To account for particle interactions in the extracellular (physical) environment, information from the cell-level Monte Carlo damage simulation (MCDS) for DNA double strand break (DSB) induction has been integrated into the general purpose Monte Carlo N-particle (MCNP) radiation transport code system. The effort to integrate these models is motivated by the need for a computationally efficient model to accurately predict particle relative biological effectiveness (RBE) in cell cultures and in vivo. To illustrate the approach and highlight the impact of the larger scale physical environment (e.g. establishing charged particle equilibrium), we examined the RBE for DSB induction (RBEDSB) of x-rays, (137)Cs γ-rays, neutrons and light ions relative to γ-rays from (60)Co in monolayer cell cultures at various depths in water. Under normoxic conditions, we found that (137)Cs γ-rays are about 1.7% more effective at creating DSB than γ-rays from (60)Co (RBEDSB  =  1.017) whereas 60-250 kV x-rays are 1.1 to 1.25 times more efficient at creating DSB than (60)Co. Under anoxic conditions, kV x-rays may have an RBEDSB up to 1.51 times as large as (60)Co γ-rays. Fission neutrons passing through monolayer cell cultures have an RBEDSB that ranges from 2.6 to 3.0 in normoxic cells, but may be as large as 9.93 for anoxic cells. For proton pencil beams, Monte Carlo simulations suggest an RBEDSB of about 1.2 at the tip of the Bragg peak and up to 1.6 a few mm beyond the Bragg peak. Bragg peak RBEDSB increases with decreasing oxygen concentration, which may create opportunities to apply proton dose painting to help address tumor hypoxia. Modeling of the particle RBE for DSB induction across multiple physical and biological scales has the potential to aid in the interpretation of laboratory experiments and provide useful information to advance the safety and effectiveness of hadron therapy in the treatment of cancer.

A multi-scale Monte Carlo model is proposed to assess the dosimetric and biological impact of iodine-based contrast agents commonly used in computed tomography. As presented, the model integrates the general purpose MCNP6 code system for larger-scale radiation transport and dose assessment with the Monte Carlo damage simulation to determine the sub-cellular characteristics and spatial distribution of initial DNA damage. The repair-misrepair-fixation model is then used to relate DNA double strand break (DSB) induction to reproductive cell death. Comparisons of measured and modeled changes in reproductive cell survival for ultrasoft characteristic k-shell x-rays (0.25-4.55 keV) up to orthovoltage (200-500 kVp) x-rays indicate that the relative biological effectiveness (RBE) for DSB induction is within a few percent of the RBE for cell survival. Because of the very short range of secondary electrons produced by low energy x-ray interactions with contrast agents, the concentration and subcellular distribution of iodine within and near cellular targets have a significant impact on the estimated absorbed dose and number of DSB produced in the cell nucleus. For some plausible models of the cell-level distribution of contrast agent, the model predicts an increase in RBE-weighted dose (RWD) for the endpoint of DSB induction of 1.22-1.40 for a 5-10 mg ml iodine concentration in blood compared to an RWD increase of 1.07  ±  0.19 from a recent clinical trial. The modeled RWD of 2.58  ±  0.03 is also in good agreement with the measured RWD of 2.3  ±  0.5 for an iodine concentration of 50 mg ml relative to no iodine. The good agreement between modeled and measured DSB and cell survival estimates provides some confidence that the presented model can be used to accurately assess biological dose for other concentrations of the same or different contrast agents.

A gradient of negative surface charge based on 1-D spatial variation from surface sulfhydryl to mixed sulfhydryl-sulfonate moities was prepared by controlled UV oxidation of 3-mercaptopropylsilane monolayer on fused silica. Adsorption of three human plasma proteins, albumin (HSA), immunoglobulin G (IgG), and fibrinogen (Fgn) onto such surface gradient was studied using spatially-resolved total internal reflection fluorescence (TIRF) and autoradiography. Adsorption was measured from dilute solutions equivalent to 1/100 (TIRF, autoradiography), and 1/500 and 1/1000 (autoradiography) of protein's physiological concentrations in plasma. All three proteins adsorbed more to the non-oxidized sulfhydryl region than to the oxidized, mixed sulfhydryl-sulfonate region of the gradient. In the case of HSA the adsorption contrast along the gradient was largest when the adsorption took place from more dilute protein solutions. Increasing the concentration to 1/100 of protein plasma concentration eliminated the effect of the gradient on HSA adsorption and to the lesser extent on IgG adsorption. In the case of Fgn the greatest adsorption contrast was observed at the highest concentration used. Based on adsorption kinetics, the estimated binding affinity of HSA for the sulfhydryl region what twice the affinity for the mixed sulfhydryl-sulfonate region of the gradient. For IgG and Fgn the initial adsorption was transport-limited and the initial adsorption rates approached the computed flux of the protein to the surface.

Accurate dosimetry and determination of the biological effectiveness of boron neutron capture therapy (BNCT) is challenging because of the mix of different types and energies of radiation at the cellular and subcellular levels. In this paper, we present a computational, multiscale system of models to better assess the relative biological effectiveness (RBE) and compound biological effectiveness (CBE) of several neutron sources as applied to BNCT using boronophenylalanine (BPA) and a potential monoclonal antibody (mAb) that targets HER-2-positive cells with Trastuzumab. The multiscale model is tested against published in vitro and in vivo measurements of cell survival with and without boron. The combined dosimetric and radiobiological model includes an analytical formulation that accounts for the type of neutron source, the tissue- or cancer-specific dose–response characteristics, and the microdistribution of boron. Tests of the model against results from published experiments with and without boron show good agreement between modeled and experimentally determined cell survival for neutrons alone and in combination with boron. The system of models developed in this work is potentially useful as an aid for the optimization and individualization of BNCT for HER-2-positive cancers, as well as other cancers, that can be targeted with mAb or a conventional BPA compound.

No abstract