Assessment of skin dose in breast cancer radiotherapy: on-phantom measurement and Monte Carlo simulation

AimThe main purpose of the present study is assessment of skin dose in breast cancer radiotherapy.BackgroundAccurate assessment of skin dose in radiotherapy can provide useful information for clinical considerations.Materials and methodsA RANDO phantom was irradiated using a 6 MV Siemens Primus linac with medial and tangential radiotherapy fields for simulating breast cancer treatment. Dosimetry was also performed on various positions across the fields using an EBT3 radiochromic film. Similar conditions of measurement on the RANDO phantom including field size, irradiation angle, number of fields, etc. were subsequently simulated via the Monte Carlo N-Particle Transport code (MCNP). Ultimately, dose values for corresponding points from both methods were compared.ResultsConsidering dosimetry using radiochromic films on the RANDO phantom, there were points having underdose and overdose based on the prescribed dose and skin tolerance levels. In this respect, 81.25% and 18.75% of the points had underdose and overdose, respectively. In some cases, several differences were observed between the measurement and the MCNP simulation results associated with skin dose.ConclusionBased on the results of the points which had underdose, it was suggested that a bolus should be used for the given points. With regard to overdose points, it was advocated to consider skin tolerance dose in treatment planning. Differences between the measurement and the MCNP simulation results might be due to voxel size of tally cells in simulations, effect of beam’s angle of incidence, validation time of linac’s head, lack of electronic equilibrium in the build-up region, as well as MCNP tally type.     

Impact of lung density on isolated lung tumor dose in VMAT using inline MR-Linac

PurposeThis study investigated the impact of lung density on the isolated lung tumor dose for volumetric modulated arc therapy (VMAT) in an inline magnetic resonance linear accelerator (MR-Linac) using the Monte Carlo (MC) simulation.MethodsCT images of the thorax phantoms with lung tumors of 1, 2, and 3 cm diameters were converted into voxel-base phantoms with lung densities of 0.1, 0.2, and 0.3 g/cm³, respectively. The dose distributions were calculated for partial-arc VMAT. The dose distributions were compared using dose differences, dose volume histograms, and dose volume indices.ResultsIn all cases, the inline magnetic field significantly enhanced the lung tumor dose compared to that at 0 T. For the 1 cm lung tumor, the inline magnetic field of 1 T increased the minimum dose of 95% of the Planning target volume (PTV D₉₅) by 14.0% in 0.1 g/cm³ lung density as compared to that in 0.3 g/cm³ at 0 T. In contrast, at 0 and 0.5 T, the PTV D₉₅ in 0.3 g/cm³ lung density was larger than that in lung density of 0.1 g/cm³. For the 2 cm lung tumor, a similar tendency to 1 cm was observed, whereas the dose impact of lung density was smaller than that for 1 cm. For the 3 cm lung tumor, the lung tumor dose was independent of lung density at 0.5 T and 1.0 T.ConclusionThe inline MR-Linac with the magnetic field over 1 T can enhance the PTV D₉₅ for VMAT regardless of the lung density.     

An on-lattice agent-based Monte Carlo model simulating the growth kinetics of multicellular tumor spheroids

PurposeTo develop an on-lattice agent-based model describing the growth of multicellular tumor spheroids using simple Monte Carlo tools.MethodsCells are situated on the vertices of a cubic grid. Different cell states (proliferative, hypoxic or dead) and cell evolution rules, driven by 10 parameters, and the effects of the culture medium are included. About twenty spheroids of MCF-7 human breast cancer were cultivated and the experimental data were used for tuning the model parameters.ResultsSimulated spheroids showed adequate sizes of the necrotic nuclei and of the hypoxic and proliferative cell phases as a function of the growth time, mimicking the overall characteristics of the experimental spheroids. The relation between the radii of the necrotic nucleus and the whole spheroid obtained in the simulations was similar to the experimental one and the number of cells, as a function of the spheroid volume, was well reproduced. The statistical variability of the Monte Carlo model described the whole volume range observed for the experimental spheroids. Assuming that the model parameters vary within Gaussian distributions it was obtained a sample of spheroids that reproduced much better the experimental findings.ConclusionsThe model developed allows describing the growth of in vitro multicellular spheroids and the experimental variability can be well reproduced. Its flexibility permits to vary both the agents involved and the rules that govern the spheroid growth. More general situations, such as, e. g., tumor vascularization, radiotherapy effects on solid tumors, or the validity of the tumor growth mathematical models can be studied.     

One size does not fit all: Customizing MCMC methods for hierarchical models using NIMBLE

Improved efficiency of Markov chain Monte Carlo facilitates all aspects of statistical analysis with Bayesian hierarchical models. Identifying strategies to improve MCMC performance is becoming increasingly crucial as the complexity of models, and the run times to fit them, increases. We evaluate different strategies for improving MCMC efficiency using the open‐source software NIMBLE (R package nimble) using common ecological models of species occurrence and abundance as examples. We ask how MCMC efficiency depends on model formulation, model size, data, and sampling strategy. For multiseason and/or multispecies occupancy models and for N‐mixture models, we compare the efficiency of sampling discrete latent states vs. integrating over them, including more vs. fewer hierarchical model components, and univariate vs. block‐sampling methods. We include the common MCMC tool JAGS in comparisons. For simple models, there is little practical difference between computational approaches. As model complexity increases, there are strong interactions between model formulation and sampling strategy on MCMC efficiency. There is no one‐size‐fits‐all best strategy, but rather problem‐specific best strategies related to model structure and type. In all but the simplest cases, NIMBLE's default or customized performance achieves much higher efficiency than JAGS. In the two most complex examples, NIMBLE was 10–12 times more efficient than JAGS. We find NIMBLE is a valuable tool for many ecologists utilizing Bayesian inference, particularly for complex models where JAGS is prohibitively slow. Our results highlight the need for more guidelines and customizable approaches to fit hierarchical models to ensure practitioners can make the most of occupancy and other hierarchical models. By implementing model‐generic MCMC procedures in open‐source software, including the NIMBLE extensions for integrating over latent states (implemented in the R package nimbleEcology), we have made progress toward this aim.     

A detailed model and Monte Carlo simulation for predicting DIP genome length distribution in baculovirus infection of insect cells

Baculoviruses have enormous potential for use as biopesticides to control insect pest populations without the adverse environmental effects posed by the widespread use of chemical pesticides. However, continuous baculovirus production is susceptible to DNA mutation and the subsequent production of defective interfering particles (DIPs). The amount of DIPs produced and their genome length distribution are of great interest not only for baculoviruses but for many other DNA and RNA viruses. In this study, we elucidate this aspect of virus replication using baculovirus as an example system and both experimental and modeling studies. The existing mathematical models for the virus replication process consider DIPs as a lumped quantity and do not consider the genome length distribution of the DIPs. In this study, a detailed population balance model for the cell‐virus culture is presented, which predicts the genome length distribution of the DIP population along with their relative proportion. The model is simulated using the kinetic Monte Carlo algorithm, and the results agree well with the experimental results. Using this model, a practical strategy to maintain the DIP fraction to near to its maximum and minimum limits has been demonstrated.     

Comparison of the RayStation photon Monte Carlo dose calculation algorithm against measured data under homogeneous and heterogeneous irradiation geometries

PurposeThis work compares Monte Carlo dose calculations performed using the RayStation treatment planning system against data measured on a Varian Truebeam linear accelerator with 6 MV and 10 MV FFF photon beams.MethodsThe dosimetric performance of the RayStation Monte Carlo calculations was evaluated in a variety of irradiation geometries employing homogeneous and heterogeneous phantoms. Profile and depth dose comparisons against measurement were carried out in relative mode using the gamma index as a quantitative measure of similarity within the central high dose regions.ResultsThe results demonstrate that the treatment planning system dose calculation engine agrees with measurement to within 2%/1 mm for more than 95% of the data points in the high dose regions for all test cases. A systematic underestimation was observed at the tail of the profile penumbra and out of field, with mean differences generally <0.5 mm or 1% of curve dose maximum respectively. Out of field agreement varied between evaluated beam models.ConclusionsThe RayStation implementation of photon Monte Carlo dose calculations show good agreement with measured data for the range of scenarios considered in this work and is deemed sufficiently accurate for introduction into clinical use.     

CPOP: An open source C++ cell POPulation modeler for radiation biology applications

PurposeMulticellular tumor spheroids are realistic in-vitro systems used in radiation biology research to study the effect of anticancer drugs or to evaluate the resistance of cancer cells under specific conditions. When combining the modeling of spheroids together with the simulation of radiation using Monte Carlo methods, one could estimate cell and DNA damage to be compared with experimental data. We developed a Cell Population (CPOP) modeler combined to Geant4 simulations in order to tackle how energy depositions are allocated to cells, especially when enhancing radiation outcomes using high-Z nanoparticles. CPOP manages to model large three-dimensional cell populations with independent deformable cells described with their nucleus, cytoplasm and membranes together with force law systems to manage cell–cell interactions.MethodsCPOP is an opensource platform written in C++. It is divided into two main libraries: a “Modeler” library, for cell geometry modeling using meshes, and a Multi Agent System (MAS) library, simulating all agent (cell) interactions among the population. CPOP is fully interfaced with the Geant4 Monte Carlo toolkit and is able to directly launch Geant4 simulations after compilation.We modeled a full and realistic 3D cell population from SK-MEL28 melanoma cell population cultured experimentally. The spheroid diameter of 550 ± 40 µm corresponds to a population of approximately 1000 cells having a diameter of 17.2 ± 2.5 µm and a nucleus diameter of 11.2 ± 2.0 µm. We decided to reproduce cell irradiations performed with a X-RAD 320 Biological Irradiator (Precision XRay Inc., North Branford, CT).ResultsWe simulated the energy spectrum of secondary particles generated in the vicinity of the spheroid and plotted the different energy spectra recovered internally to the spheroid. We evaluated also the impact of AGuIX (Gadolinium) nanoparticles modeled into the spheroid with their corresponding secondary energy spectra.ConclusionsWe succeeded into modeling cell populations and combined them with Geant4 simulations. The next step will be to integrate DNA geometrical models into cell nuclei and to use the Geant4-DNA physics and radiolysis modeling capabilities in order to evaluate early strand breaks induced on DNA.     

Ricerche Morfologiche e Anatomiche Su Mimosa Pudica L.

Summary

Firstly, the author describes the anatomical and histological structures of Mimosa pudica L. during the stages of development of the young plant prior to the appearing of cambium. She then discusses, in the light of previous literature on the subject, the three principal features of Mimosa pudica: herbaceous bearing, stable tetrarchy with intercotyledonal bundles along the plumule, and small sized germination.

Finally, the author concludes that Mimosa pudica appears to be an individual with tree-like rather than herbaceous symmetry (stable tetrarchy), in which the size of germination has diminished parallelly to reduction in size of the adult individual.     

Effect of anode angle on photon beam spectra and depth dose characteristics for X-RAD320 orthovoltage unit

Background

In radiation therapy with orthovoltage units, the tube design has a crucial effect on its dosimetric features.

Aim

In this study, the effect of anode angle on photon beam spectra, depth dose and photon fluence per initial electron was studied for a commercial orthovoltage unit of X-RAD320 biological irradiator.

Materials and methods

The MCNPX MC code was used for modeling in the current study. We used the Monte Carlo method to model the X-RAD320 X-ray unit based on the manufacturer provided information. The MC model was validated by comparing the MC calculated photon beam spectra with the results of SpekCalc software. The photon beam spectra were calculated for anode angles from 15 to 35 degrees. We also calculated the percentage depth doses for some angles to verify the impact of anode angle on depth dose. Additionally, the heel effect and its relation with anode angle were studied for X-RAD320 irradiator.

Results

Our results showed that the photon beam spectra and their mean energy are changed significantly with anode angle and the optimum anode angle of 30 degrees was selected based on less heel effect and appropriate depth dose and photon fluence per initial electron.

Conclusion

It can be concluded that the anode angle of 30 degrees for X-RAD320 unit used by manufacturer has been selected properly considering the heel effect and dosimetric properties.     

Hemocyanin gene family evolution in spiders (Araneae), with implications for phylogenetic relationships and divergence times in the infraorder Mygalomorphae

Hemocyanins are multimeric copper-containing hemolymph proteins involved in oxygen binding and transport in all major arthropod lineages. Most arachnids have seven primary subunits (encoded by paralogous genes a–g), which combine to form a 24-mer (4 × 6) quaternary structure. Within some spider lineages, however, hemocyanin evolution has been a dynamic process with extensive paralog duplication and loss. We have obtained hemocyanin gene sequences from numerous representatives of the spider infraorders Mygalomorphae and Araneomorphae in order to infer the evolution of the hemocyanin gene family and estimate spider relationships using these conserved loci. Our hemocyanin gene tree is largely consistent with the previous hypotheses of paralog relationships based on immunological studies, but reveals some discrepancies in which paralog types have been lost or duplicated in specific spider lineages. Analyses of concatenated hemocyanin sequences resolved deep nodes in the spider phylogeny and recovered a number of clades that are supported by other molecular studies, particularly for mygalomorph taxa. The concatenated data set is also used to estimate dates of higher-level spider divergences and suggests that the diversification of extant mygalomorphs preceded that of extant araneomorphs. Spiders are diverse in behavior and respiratory morphology, and our results are beneficial for comparative analyses of spider respiration. Lastly, the conserved hemocyanin sequences allow for the inference of spider relationships and ancient divergence dates.