In-silico calculations of DNA damage induced by α-particles in the 224Ra DaRT decay chain for a better understanding of the radiobiological effectiveness of this treatment.

Diffusing alpha-emitters radiation Therapy (DaRT) is an interstitial brachytherapy technique using 224Ra seeds. For accurate treatment planning a good understanding of the early DNA damage due to α-particles is required. Geant4-DNA was used to calculate the initial DNA damage and radiobiological effectiveness due to α-particles with linear energy transfer (LET) values in the range 57.5-225.9 keV/µm from the 224Ra decay chain. The impact of DNA base pair density on DNA damage has been modelled, as this parameter varies between human cell lines. Results show that the quantity and complexity of DNA damage changes with LET as expected. Indirect damage, due to water radical reactions with the DNA, decreases and becomes less significant at higher LET values as shown in previous studies. As expected, the yield of complex double strand breaks (DSBs), which are harder for a cell to repair, increases approximately linearly with LET. The level of complexity of DSBs and radiobiological effectiveness have been found to increase with LET as expected. The quantity of DNA damage has been shown to increase for increased DNA density in the expected base pair density range of human cells. The change in damage yield as a function of base pair density is largest for higher LET α-particles, an increase of over 50% for individual strand breaks between 62.7 and 127.4 keV/µm. This change in yield shows that the DNA base pair density is an important parameter for modelling DNA damage particularly at higher LET where the DNA damage is greatest and most complex.

Geant4-DNA simulation of human cancer cells irradiation with helium ion beams.

PURPOSE: This study aimed to develop a computational environment for the accurate simulation of human cancer cell irradiation using Geant4-DNA. New cell geometrical models were developed and irradiated by alpha particle beams to induce DNA damage. The proposed approach may help further investigation of the benefits of external alpha irradiation therapy. METHODS: The Geant4-DNA Monte Carlo (MC) toolkit allows the simulation of cancer cell geometries that can be combined with accurate modelling of physical, physicochemical and chemical stages of liquid water irradiation, including radiolytic processes. Geant4-DNA is used to calculate direct and non-direct DNA damage yields, such as single and double strand breaks, produced by the deposition of energy or by the interaction of DNA with free radicals. RESULTS: In this study, the "molecularDNA" example application of Geant4-DNA was used to quantify early DNA damage in human cancer cells upon irradiation with alpha particle beams, as a function of linear energy transfer (LET). The MC simulation results are compared to experimental data, as well as previously published simulation data. The simulation results agree well with the experimental data on DSB yields in the lower LET range, while the experimental data on DSB yields are lower than the results obtained with the "molecularDNA" example in the higher LET range. CONCLUSION: This study explored and demonstrated the possibilities of the Geant4-DNA toolkit together with the "molecularDNA" example to simulate the helium beam irradiation of cancer cell lines, to quantify the early DNA damage, or even the following DNA damage response.

The general-purpose Geant4 Monte Carlo toolkit and its Geant4-DNA extension to investigate mechanisms underlying the FLASH effect in radiotherapy: Current status and challenges.

FLASH radiotherapy is a promising approach to cancer treatment that offers several advantages over conventional radiotherapy. With this novel technique, high doses of radiation are delivered in a short period of time, inducing the so-called FLASH effect - a phenomenon characterized by healthy tissue sparing without alteration of tumor control. The mechanisms behind the FLASH effect remain unknown. One way to approach this problem is to gain insight into the initial parameters that can distinguish FLASH from conventional irradiation by simulating particle transport in aqueous media using the general-purpose Geant4 Monte Carlo toolkit and its Geant4-DNA extension. This review article discusses the current status of Geant4 and Geant4-DNA simulations to investigate mechanisms underlying the FLASH effect, as well as the challenges faced in this research field. One of the primary challenges is to accurately simulate the experimental irradiation parameters. Another challenge is the temporal extension of the simulations. This review also focuses on two hypotheses to explain the FLASH effect - namely the oxygen depletion hypothesis and the inter-track interactions hypothesis - and discusses how the Geant4 toolkit can be used to investigate them. The aim of this review is to provide an overview of Geant4 and Geant4-DNA simulations for FLASH radiotherapy and to highlight the challenges that need to be overcome in order to better study the FLASH effect.

Space radiation quality factor for Galactic Cosmic Rays and typical space mission scenarios using a microdosimetric approach.

Space radiation exposure from omnipresent Galactic Cosmic Rays (GCRs) in interplanetary space poses a serious carcinogenic risk to astronauts due to the-limited or absent-protective effect of the Earth's magnetosphere and, in particular, the terrestrial atmosphere. The radiation risk is directly influenced by the quality of the radiation, i.e., its pattern of energy deposition at the micron/DNA scale. For stochastic biological effects, radiation quality is described by the quality factor, [Formula: see text], which can be defined as a function of Linear Energy Transfer (LET) or the microdosimetric lineal energy ([Formula: see text]). In the present work, the average [Formula: see text] of GCR for different mission scenarios was calculated using a modified version of the microdosimetric Theory of Dual Radiation Action (TDRA). NASA's OLTARIS platform was utilized to generate the radiation environment behind different aluminum shielding (0-30 g/cm2) for a typical mission scenario in low-earth orbit (LEO) and in deep space. The microdosimetric lineal energy spectra of ions ([Formula: see text]) in 1 µm liquid water spheres were calculated by a generalized analytical model which considers energy-loss fluctuations and δ-ray transport inside the irradiated medium. The present TDRA-based [Formula: see text]-values for the LEO and deep space missions were found to differ by up to 10% and 14% from the corresponding ICRP-based [Formula: see text]-values and up to 3% and 6% from NASA's [Formula: see text]-model. In addition, they were found to be in good agreement with the [Formula: see text]-values measured in the International Space Station (ISS) and by the Mars Science Laboratory (MSL) Radiation Assessment Detector (RAD) which represent, respectively, a LEO and deep space orbit.

Modeling of scavenging systems in water radiolysis with Geant4-DNA.

PURPOSE: This paper presents the capabilities of the Geant4-DNA Monte Carlo toolkit to simulate water radiolysis with scavengers using the step-by-step (SBS) or the independent reaction times (IRT) methods. It features two examples of application areas: (1) computing the escape yield of H2O2 following a 60Co γ-irradiation and (2) computing the oxygen depletion in water irradiated with 1 MeV electrons. METHODS: To ease the implementation of the chemical stage in Geant4-DNA, we developed a user interface that helps define the chemical reactions and set the concentration of scavengers. The first application area example required two computational steps to perform water radiolysis using NO2- and NO3- as scavengers and a 60Co irradiation. The oxygen depletion computation technique for the second application area example consisted of simulating track segments of 1 MeV electrons and determining the radio-induced loss and gain of oxygen molecules. RESULTS: The production of H2O2 under variable scavenging levels is consistent with the literature; the mean relative difference between the SBS and IRT methods is 7.2 % ± 0.5 %. For the oxygen depletion 1 µs post-irradiation, the mean relative difference between both methods is equal to 9.8 % ± 0.3 %. The results in the microsecond scale depend on the initial partial pressure of oxygen in water. In addition, the computed oxygen depletions agree well with the literature. CONCLUSIONS: The Geant4-DNA toolkit makes it possible to simulate water radiolysis in the presence of scavengers. This feature offers perspectives in radiobiology, with the possibility of simulating cell-relevant scavenging mechanisms.

Sub-keV corrections to binary encounter cross section models for electron ionization of liquid water with application to the Geant4-DNA Monte Carlo code.

INTRODUCTION: The electron ionization cross section of water is one of the most important input in Monte Carlo studies of cellular radiobiological effects. Analytical cross section models of the binary-encounter type have the potential of reducing simulation time and facilitate application to a variety of biological materials (other than water). The Binary-Encounter-Bethe (BEB) and Binary-Encounter-Dipole (BED) models of NIST are perhaps the most popular of such models giving reliable results for atoms and molecules in the gas-phase over a wide energy range. However, the use of such models to sub-keV electron energies in liquid water raises concerns due to the neglect of condensed phase effects that leads to a significant overestimation when compared to medium-specific dielectric models. PURPOSE: To modify the BEB and BED models towards better agreement with the recommended low-energy dielectric model of Geant4-DNA (Option 4). To implement the new modifications to the existing BEB model of the Option 6 physics constructor of Geant4-DNA and re-evaluate fundamental transport quantities for sub-keV electrons. METHODS: In analogy to a Yukawa potential a simple, yet physically-motivated, modification of the Burgess correction term is proposed to account for the reduction of the Coulomb interaction due to the polarizability of the target. The magnitude of the correction is guided by the dielectric-based ionization cross section implemented in Option 4. RESULTS: Differential, total and stopping ionization cross sections for low-energy electrons in liquid water are presented. When combined with the Vriens correction (which is not included in Option 6), the proposed modification to the BEB and BED models brings the ionization and stopping cross sections in much better agreement against those used in the Option 4 dielectric model of Geant4-DNA, with up to 30% and 10% deviation, respectively. Implementation of the new correction to the Option 6 constructor of Geant4-DNA and re-evaluation of fundamental transport quantities, such as electron penetration ranges and dose-point-kernels, reduced the discrepancies from Option 4 at sub-keV energies from 20 to 100% (or more) to well below 10% in most cases. CONCLUSIONS: A simple modification to the BEB and BED analytic models was found to improve their performance for sub-keV electrons in liquid water medium. Implementation of the new modification to the Option 6 constructor of Geant4-DNA significantly improved the agreement with the recommended low-energy Option 4 constructor for a variety of fundamental quantities related to electron transport.

Prediction of DNA rejoining kinetics and cell survival after proton irradiation for V79 cells using Geant4-DNA.

PURPOSE: Track structure Monte Carlo (MC) codes have achieved successful outcomes in the quantitative investigation of radiation-induced initial DNA damage. The aim of the present study is to extend a Geant4-DNA radiobiological application by incorporating a feature allowing for the prediction of DNA rejoining kinetics and corresponding cell surviving fraction along time after irradiation, for a Chinese hamster V79 cell line, which is one of the most popular and widely investigated cell lines in radiobiology. METHODS: We implemented the Two-Lesion Kinetics (TLK) model, originally proposed by Stewart, which allows for simulations to calculate residual DNA damage and surviving fraction along time via the number of initial DNA damage and its complexity as inputs. RESULTS: By optimizing the model parameters of the TLK model in accordance to the experimental data on V79, we were able to predict both DNA rejoining kinetics at low linear energy transfers (LET) and cell surviving fraction. CONCLUSION: This is the first study to demonstrate the implementation of both the cell surviving fraction and the DNA rejoining kinetics with the estimated initial DNA damage, in a realistic cell geometrical model simulated by full track structure MC simulations at DNA level and for various LET. These simulation and model make the link between mechanistic physical/chemical damage processes and these two specific biological endpoints.

Electron tracks simulation in water: Performance comparison between GPU CPU and the EUMED grid installation.

PURPOSE: We explored different technologies to minimize simulation time of the Monte-Carlo method for track generation following the Geant4-DNA processes for electrons in water. METHODS: A GPU software tool is developed for electron track simulations. A similar CPU version is also developed using the same collision models. CPU simulations were carried out on a single user desktop computer and on the computing grid France Grilles using 10 and 100 computing nodes. Computing time results for CPU, GPU, and grid simulations are compared with those using Geant4-DNA processes. RESULTS: The CPU simulations better performs when the number of electrons is less than 104 with 100 eV initial energy, this number decreases as the energy increases. The GPU simulations gives better results when the number of electrons is more than 104 with initial energy of 100 eV, this number decreases to 103 for electrons with 10KeV and increases back with higher energy. The use of the grid introduces an additional queuing time which slows down the overall simulation performance. Thus, the Grid gives better performance when the number of electrons is over 105 with initial energy of 10KeV, and this number decreases as the energy increases. CONCLUSIONS: The CPU is best suited for small numbers of primary incident electrons. The GPU is best suited when the number of primary incident particles occupies sufficient resources on GPU card in order to get an important computing power. The grid is best suited for simulations with high number of primary incident electrons with high initial energy.

Intercomparison of nanodosimetric distributions in nitrogen simulated with Geant4 and PTra track structure codes.

To facilitate the use of Geant4-DNA for radiation transport simulations in micro- and nanodosimeters, which are physically operated with tissue-equivalent gases such as nitrogen (and propane), this work aims to extend the cross section data available in Geant4-DNA to include those of nitrogen for electron energies ranging from 1 MeV down to the ionisation threshold. To achieve this, interaction cross section data for nitrogen that have been used with the in-house PTB PTra track structure code have been implemented in the current state-of-the-art Geant4-DNA simulation toolkit. An intercomparison has been performed between the two codes to validate this implementation. To quantify the agreement between the cross section models for nitrogen adopted in PTra and those implemented in Geant4-DNA, the simulation results of both codes were analysed using three physical parameters describing the ionisation cluster size distribution (ICSD): mean ionisation cluster size, variance of the cluster size and the probability to obtain a single ionisation within the target. Statistical analysis of the results indicates that the interaction cross section models for nitrogen used in PTra (elastic scattering, impact ionisations and electronic excitations) have been successfully implemented in Geant4-DNA. In addition, simulated ICSDs were compared to those measured with the Jet Counter nanodosimeter for energies between 100 and 2000 eV. For greater energies, the ICRP data for LET and particle range were used as a reference. The modified Geant4-DNA code and data successfully passed all these benchmarks fulfilling the requirement for their public release in the next version of the Geant4 toolkit.

Nanodosimetric Calculations of Radiation-Induced DNA Damage in a New Nucleus Geometrical Model Based on the Isochore Theory.

Double-strand breaks (DSBs) in nuclear DNA represents radiation-induced damage that has been identified as particularly deleterious. Calculating this damage using Monte Carlo track structure modeling could be a suitable indicator to better assess and anticipate the side-effects of radiation therapy. However, as already demonstrated in previous work, the geometrical description of the nucleus and the DNA content used in the simulation significantly influence damage calculations. Therefore, in order to obtain accurate results, this geometry must be as realistic as possible. In this study, a new geometrical model of an endothelial cell nucleus and DNA distribution according to the isochore theory are presented and used in a Monte Carlo simulation chain based on the Geant4-DNA toolkit. In this theory, heterochromatin and euchromatin compaction are distributed along the genome according to five different families (L1, L2, H1, H2, and H3). Each of these families is associated with a different hetero/euchromatin rate related to its compaction level. In order to compare the results with those obtained using a previous nuclear geometry, simulations were performed for protons with linear energy transfers (LETs) of 4.29 keV/µm, 19.51 keV/µm, and 43.25 keV/µm. The organization of the chromatin fibers at different compaction levels linked to isochore families increased the DSB yield by 6-10%, and it allowed the most affected part of the genome to be identified. These new results indicate that the genome core is more radiosensitive than the genome desert, with a 3-8% increase in damage depending on the LET. This work highlights the importance of using realistic distributions of chromatin compaction levels to calculate radio-induced damage using Monte Carlo simulation methods.

Implementation of the EPICS2017 database for photons in Geant4.

This paper describes in detail the implementation of Geant4 Livermore electromagnetic physics models based on the EPICS2017 database for the low energy transport of photons. These models describe four photon processes: gamma conversion, Compton scattering, photoelectric effect and Rayleigh scattering. New parameterizations based on EPICS2017 were performed for scattering functions of Compton effect, subshell cross-sections of the photoelectric effect and form factors of Rayleigh scattering, in order to improve the precision of fitted values compared to tabulated values. Comparisons between new and old parameterizations were also carried out to evaluate the precision of the new parameterizations. The models were tested through a comparative study, in which the mass attenuation coefficient was calculated for both total photon interaction and each process using Geant4 simulations based on EPICS2017 and EPDL97 respectively. The results obtained from the simulations were found in good agreement with the XCOM reference data.

A Geant4 simulation of X-ray emission for three-dimensional proton imaging of microscopic samples.

PURPOSE: Proton computed microtomography is a technique that reveals the inner content of microscopic samples. The density distribution of the material (in g·cm-3) is obtained from proton transmission tomography (STIM: Scanning Transmission Ion Microscopy) and the element content from X-ray emission tomography (PIXE: Particle Induced X-ray Emission). A precise quantification of chemical elements is difficult for thick samples, because of the variations of X-ray production cross-sections and of X-ray absorption. Both phenomena are at the origin of an attenuation of the measured X-ray spectra, which leads to an underestimation of the element content. Our aim is to quantify the accuracy of a specific correction method that we designed for thick samples. METHODS: In this study, we describe how the 3D variations in the mass density were taken into account in the reconstruction code, in order to quantify the correction according to the position of the proton beam and the position and aperture angle of the X-ray detector. Moreover, we assess the accuracy of the reconstructed densities using Geant4 simulations on numerical phantoms, used as references. RESULTS: The correction process was successfully applied and led, for the largest regions of interest (little affected by partial volume effects), to an accuracy ≤ 4% for phosphorus (compared to about 40% discrepancy without correction). CONCLUSION: This study demonstrates the accuracy of the correction method implemented in the tomographic reconstruction code for thick samples. It also points out some advantages offered by Geant4 simulations: i) they produce projection data that are totally independent of the inversion method used for the image reconstruction; ii) one or more physical processes (X-ray absorption, proton energy loss) can be artificially turned off, in order to precisely quantify the effect of the different phenomena involved in the attenuation of X-ray spectra.

Performance Evaluation for Repair of HSGc-C5 Carcinoma Cell Using Geant4-DNA.

Track-structure Monte Carlo simulations are useful tools to evaluate initial DNA damage induced by irradiation. In the previous study, we have developed a Gean4-DNA-based application to estimate the cell surviving fraction of V79 cells after irradiation, bridging the gap between the initial DNA damage and the DNA rejoining kinetics by means of the two-lesion kinetics (TLK) model. However, since the DNA repair performance depends on cell line, the same model parameters cannot be used for different cell lines. Thus, we extended the Geant4-DNA application with a TLK model for the evaluation of DNA damage repair performance in HSGc-C5 carcinoma cells which are typically used for evaluating proton/carbon radiation treatment effects. For this evaluation, we also performed experimental measurements for cell surviving fractions and DNA rejoining kinetics of the HSGc-C5 cells irradiated by 70 MeV protons at the cyclotron facility at the National Institutes for Quantum and Radiological Science and Technology (QST). Concerning fast- and slow-DNA rejoining, the TLK model parameters were adequately optimized with the simulated initial DNA damage. The optimized DNA rejoining speeds were reasonably agreed with the experimental DNA rejoining speeds. Using the optimized TLK model, the Geant4-DNA simulation is now able to predict cell survival and DNA-rejoining kinetics for HSGc-C5 cells.

A Geant4-DNA Evaluation of Radiation-Induced DNA Damage on a Human Fibroblast.

Accurately modeling the radiobiological mechanisms responsible for the induction of DNA damage remains a major scientific challenge, particularly for understanding the effects of low doses of ionizing radiation on living beings, such as the induction of carcinogenesis. A computational approach based on the Monte Carlo technique to simulate track structures in a biological medium is currently the most reliable method for calculating the early effects induced by ionizing radiation on DNA, the primary cellular target of such effects. The Geant4-DNA Monte Carlo toolkit can simulate not only the physical, but also the physico-chemical and chemical stages of water radiolysis. These stages can be combined with simplified geometric models of biological targets, such as DNA, to assess direct and indirect early DNA damage. In this study, DNA damage induced in a human fibroblast cell was evaluated using Geant4-DNA as a function of incident particle type (gammas, protons, and alphas) and energy. The resulting double-strand break yields as a function of linear energy transfer closely reproduced recent experimental data. Other quantities, such as fragment length distribution, scavengeable damage fraction, and time evolution of damage within an analytical repair model also supported the plausibility of predicting DNA damage using Geant4-DNA.The complete simulation chain application "molecularDNA", an example for users of Geant4-DNA, will soon be distributed through Geant4.

Geant4-DNA simulation of the pre-chemical stage of water radiolysis and its impact on initial radiochemical yields.

This paper demonstrates the impact of the pre-chemical stage, especially the dissociation scheme and the associated probabilities, on water radiolysis simulation using the Geant4-DNA Monte Carlo track structure simulation toolkit. The models and parameters provided by TRACs have been collected and implemented into Geant4-DNA. In order to evaluate their influence on water radiolysis simulation, the radiochemical yields (G-values) are evaluated as a function of time and LET using the "chem6" Geant4-DNA example, and they are compared with published experimental and calculated data. The new pre-chemical models lead to a better agreement with literature data than the default pre-chemical models of Geant4-DNA, especially for OH radicals and H2O2. The revised chemistry constructor "G4EmDNAChemistry_option3" is available in Geant4-DNA version 10.7.

Geant4-DNA Modeling of Water Radiolysis beyond the Microsecond: An On-Lattice Stochastic Approach.

In this work, we use the next sub-volume method (NSM) to investigate the possibility of using the compartment-based ("on-lattice") model to simulate water radiolysis. We, first, start with a brief description of the reaction-diffusion master equation (RDME) in a spatially discretized simulation volume ("mesh"), which is divided into sub-volumes (or "voxels"). We then discuss the choice of voxel size and merging technique of a given mesh, along with the evolution of the system using the hierarchical algorithm for the RDME ("hRDME"). Since the compartment-based model cannot describe high concentration species of early radiation-induced spurs, we propose a combination of the particle-based step-by-step ("SBS") Brownian dynamics model and the compartment-based model ("SBS-RDME model") for the simulation. We, finally, use the particle-based SBS Brownian dynamics model of Geant4-DNA as a reference to test the model implementation through several benchmarks. We find that the compartment-based model can efficiently simulate the system with a large number of species and for longer timescales, beyond the microsecond, with a reasonable computing time. Our aim in developing this model is to study the production and evolution of reactive oxygen species generated under irradiation with different dose rate conditions, such as in FLASH and conventional radiotherapy.

DNA double-strand breaks in cancer cells as a function of proton linear energy transfer and its variation in time.

PURPOSE: The complex relationship between linear energy transfer (LET) and cellular response to radiation is not yet fully elucidated. To better characterize DNA damage after irradiations with therapeutic protons, we monitored formation and disappearance of DNA double-strand breaks (DNA DSB) as a function of LET and time. Comparisons with conventional γ-rays and high LET carbon ions were also performed. MATERIALS AND METHODS: In the present work, we performed immunofluorescence-based assay to determine the amount of DNA DSB induced by different LET values along the 62 MeV therapeutic proton Spread out Bragg peak (SOBP) in three cancer cell lines, i.e. HTB140 melanoma, MCF-7 breast adenocarcinoma and HTB177 non-small lung cancer cells. Time dependence of foci formation was followed as well. To determine irradiation positions, corresponding to the desired LET values, numerical simulations were carried out using Geant4 toolkit. We compared γ-H2AX foci persistence after irradiations with protons to that of γ-rays and carbon ions. RESULTS: With the rise of LET values along the therapeutic proton SOBP, the increase of γ-H2AX foci number is detected in the three cell lines up to the distal end of the SOBP, while there is a decrease on its distal fall-off part. With the prolonged incubation time, the number of foci gradually drops tending to attain the residual level. For the maximum number of DNA DSB, irradiation with protons attain higher level than that of γ-rays. Carbon ions produce more DNA DSB than protons but not substantially. The number of residual foci produced by γ-rays is significantly lower than that of protons and particularly carbon ions. Carbon ions do not produce considerably higher number of foci than protons, as it could be expected due to their physical properties. CONCLUSIONS: In situ visualization of γ-H2AX foci reveal creation of more lesions in the three cell lines by clinically relevant proton SOBP than γ-rays. The lack of significant differences in the number of γ-H2AX foci between the proton and carbon ion-irradiated samples suggests an increased complexity of DNA lesions and slower repair kinetics after carbon ions compared to protons. For all three irradiation types, there is no major difference between the three cell lines shortly after irradiations, while later on, the formation of residual foci starts to express the inherent nature of tested cells, therefore increasing discrepancy between them.

A benchmarking study of Geant4 for Auger electrons emitted by medical radioisotopes.

Auger emitting radioisotopes are of great interest in targeted radiotherapy because, once internalised in the tumour cells, they can deliver dose locally to the radiation sensitive targets, while not affecting surrounding cells. Geant4 is a Monte Carlo code widely used to characterise the physics mechanism at the basis of targeted radiotherapy. In this work, we benchmarked the modelling of the emission of Auger electrons in Geant4 deriving from the decay of 123I, 124I, 125I radionuclides against existing theoretical approaches. We also compared Geant4 against reference data in the case of 131Cs, which is of interest for brachytherapy. In the case of 125I and 131Cs, the simulation results are compared to experimental measurements as well. Good agreement was found between Geant4 and the reference data. As far as we know, this is the first study aimed to benchmark against experimental measurements the emission of Auger electrons in Geant4 for radiotherapy applications.

Monte Carlo Simulation of SARS-CoV-2 Radiation-Induced Inactivation for Vaccine Development.

Immunization with an inactivated virus is one of the strategies currently being tested towards developing a SARS-CoV-2 vaccine. One of the methods used to inactivate viruses is exposure to high doses of ionizing radiation to damage their nucleic acids. While gamma (γ) rays effectively induce lesions in the RNA, envelope proteins are also highly damaged in the process. This in turn may alter their antigenic properties, affecting their capacity to induce an adaptive immune response able to confer effective protection. Here, we modeled the effect of sparsely and densely ionizing radiation on SARS-CoV-2 using the Monte Carlo toolkit Geant4-DNA. With a realistic 3D target virus model, we calculated the expected number of lesions in the spike and membrane proteins, as well as in the viral RNA. Our findings showed that γ rays produced significant spike protein damage, but densely ionizing charged particles induced less membrane damage for the same level of RNA lesions, because a single ion traversal through the nuclear envelope was sufficient to inactivate the virus. We propose that accelerated charged particles produce inactivated viruses with little structural damage to envelope proteins, thereby representing a new and effective tool for developing vaccines against SARS-CoV-2 and other enveloped viruses.