Visualizing Rare Watson-Crick-Like Tautomeric and Anionic Mismatches in DNA and RNA

The central dogma of molecular biology relies on the correct Watson-Crick (WC) geometry of canonical deoxyribonucleic acid (DNA) dG•dC and dA•dT base pairs to replicate and transcribe genetic information with speed and an astonishing level of fidelity. In addition, the Watson-Crick geometry of canonical ribonucleic acid (RNA) rG•rC and rA•rU base pairs is highly conserved to ensure that proteins are translated with high fidelity. However, numerous other potential nucleobase tautomeric and ionic configurations are possible that can give rise to entirely new pairing modes between the nucleotide bases. Very early on, James Watson and Francis Crick recognized their importance and in 1953 postulated that if bases adopted one of their less energetically disfavored tautomeric forms (and later ionic forms) during replication it could lead to the formation of a mismatch with a Watson-Crick-like geometry and could give rise to “natural mutations.”

Since this time numerous studies have provided evidence in support of this hypothesis and have expanded upon it; computational studies have addressed the energetic feasibilities of different nucleobases’ tautomeric and ionic forms in siico; crystallographic studies have trapped different mismatches with WC-like geometries in polymerase or ribosome active sites. However, no direct evidence has been given for (i) the direct existence of these WC-like mismatches in canonical DNA duplex, RNA duplexes, or non-coding RNAs; (ii) which, if any, tautomeric or ionic form stabilizes the WC-like geometry. This thesis utilizes nuclear magnetic resonance (NMR) spectroscopy and rotating frame relaxation dispersion (R1ρ RD) in combination with density functional theory (DFT), biochemical assays, and targeted chemical perturbations to show that (i) dG•dT mismatches in DNA duplexes, as well as rG•rU mismatches RNA duplexes and non-coding RNAs, transiently adopt a WC-like geometry that is stabilized by (ii) an interconnected network of rapidly interconverting rare tautomers and anionic bases. These results support Watson and Crick’s tautomer hypothesis, but additionally support subsequent hypotheses invoking anionic mismatches and ultimately tie them together. This dissertation shows that a common mismatch can adopt a Watson-Crick-like geometry globally, in both DNA and RNA, and whose geometry is stabilized by a kinetically linked network of rare tautomeric and anionic bases. The studies herein also provide compelling evidence for their involvement in spontaneous replication and translation errors.

Design and Synthesis of Acridine-Based Macrocyclic G-Quadruplex DNA Ligands

DNA sequences that are rich in the guanine nucleic base possess the ability to fold into higher order structures called G-quadruplexes. These higher level structures are formed as a result of two sets of four guanine bases hydrogen-bonding together in a planar arrangement called a guanine quartet. Guanine quartets subsequently stack upon each other to form quadruplexes. G-quadruplexes are mainly localized in telomeres as well as in oncogene promoters. One unique and promising therapeutic approach against cancer involves targeting and stabilizing G-quadruplexes with small molecules, generally in order to suppress oncogene expression and telomerase enzyme activity; the latter has been found to contribute to “out-of control” cell growth in ca. 80-85% of all cancer cells and primary tumours while being absent in normal somatic cells. In this work, we present efforts towards designing and synthesizing acridine-based macrocycles (Mh) and (Mb) with the purpose of providing potential G4 ligands that are suited for selective binding to G4 vs. duplex DNA, and stabilize G-quadruplex structures. Two ligands described in this study include an acridine core which provides an aromatic surface capable of π-π interactions with the surface of G-quadruplexes. The successful synthesis of 4,5-diaminoacridine is described in chapter 2, as an essential fragment of the macrocycles (Mh) and (Mb). In order to investigate the synthetic method for macrocyclization, model compounds composing almost half of the designed macrocycles were explored. As discussed in chapter 3, the synthesis of the model compound for (Mb) turned out to be challenging. However, as a step towards the synthesis of (Mh), the synthesis of the hydrogen-containing model compound, which is almost half of the desired macrocycle (Mh) was achieved in our group and proved to be promising.

DNA damage induced by low energy electrons (LEEs)

Abstract : The major objective of our study is to investigate DNA damage induced by soft X-rays (1.5 keV) and low-energy electrons (˂ 30 eV) using a novel irradiation system created by Prof. Sanche’s group. Thin films of double-stranded DNA are deposited on either glass and tantalum substrates and irradiated under standard temperature and pressure surrounded by a N[subscript 2] environment. Base release (cytosine, thymine, adenine and guanine) and base modifications (8-oxo-7,8-dihydro -2’-deoxyguanosine, 5-hydroxymethyl-2’-deoxyuridine, 5-formyl-2’-deoxyuridine, 5,6-dihydrothymidine and 5,6-dihydro-2’-deoxy uridine) are analyzed and quantified by LC-MS/MS. Our results reveal larger damage yields in the sample deposited on tantalum than those on glass. This can be explained by an enhancement of damage due to low-energy electrons, which are emitted from the metal substrate. From a comparison of the yield of products, base release is the major type of damage especially for purine bases, which are 3-fold greater than base modifications. A proposed pathway leading to base release involves the formation of a transient negative ion (TNI) followed by dissociative electron attachment (DEA) at the N-g lycosidic bond. On the other hand, base modification products consist of two major types of chemical modifications, which include thymine methyl oxidation products that likely arises from DEA from the methyl group of thymine, and 5,6-dihydropyrimidine that can involve the initial addition of electrons, H atoms, or hydride ions to the 5,6-pyrimidine double bond.

Development of a simple and fast “DNA extraction kit” for sea food identification and marine species

Seafood products fraud, the misrepresentation of them, have been discovered all around the world in different forms as false labeling, species substitution, short-weighting or over glazing in order to hide the correct identity, origin or weight of the seafood products. Due to the value of seafood products such as canned tuna, swordfish or grouper, these species are the subject of the commercial fraud is mainly there placement of valuable species with other little or no value species. A similar situation occurs with the shelled shrimp or shellfish that are reduced into pieces for the commercialization. Food fraud by species substitution is an emerging risk given the increasingly global food supply chain and the potential food safety issues. Economic food fraud is committed when food is deliberately placed on the market, for financial gain deceiving consumers (Woolfe, M. & Primrose, S. 2004). As a result of the increased demand and the globalization of the seafood supply, more fish species are encountered in the market. In this scenary, it becomes essential to unequivocally identify the species. The traditional taxonomy, based primarily on identification keys of species, has shown a number of limitations in the use of the distinctive features in many animal taxa, amplified when fish, crustacean or shellfish are commercially transformed. Many fish species show a similar texture, thus the certification of fish products is particularly important when fishes have undergone procedures which affect the overall anatomical structure, such as heading, slicing or filleting (Marko et al., 2004). The absence of morphological traits, a main characteristic usually used to identify animal species, represents a challenge and molecular identification methods are required. Among them, DNA-based methods are more frequently employed for food authentication (Lockley & Bardsley, 2000). In addition to food authentication and traceability, studies of taxonomy, population and conservation genetics as well as analysis of dietary habits and prey selection, also rely on genetic analyses including the DNA barcoding technology (Arroyave & Stiassny, 2014; Galimberti et al., 2013; Mafra, Ferreira, & Oliveira, 2008; Nicolé et al., 2012; Rasmussen & Morrissey, 2008), consisting in PCR amplification and sequencing of a COI mitochondrial gene specific region. The system proposed by P. Hebert et al. (2003) locates inside the mitochondrial COI gene (cytochrome oxidase subunit I) the bioidentification system useful in taxonomic identification of species (Lo Brutto et al., 2007). The COI region, used for genetic identification - DNA barcode - is short enough to allow, with the current technology, to decode sequence (the pairs of nucleotide bases) in a single step. Despite, this region only represents a tiny fraction of the mitochondrial DNA content in each cell, the COI region has sufficient variability to distinguish the majority of species among them (Biondo et al. 2016). This technique has been already employed to address the demand of assessing the actual identity and/or provenance of marketed products, as well as to unmask mislabelling and fraudulent substitutions, difficult to detect especially in manufactured seafood (Barbuto et al., 2010; Galimberti et al., 2013; Filonzi, Chiesa, Vaghi, & Nonnis Marzano, 2010). Nowadays,the research concerns the use of genetic markers to identify not only the species and/or varieties of fish, but also to identify molecular characters able to trace the origin and to provide an effective control tool forproducers and consumers as a supply chain in agreementwith local regulations.

Analysis of radiation induced DNA damage by LC-MS/MS in isolated and cellular DNA

Abstract: It is well established that ionizing radiation induces a variety of damage in DNA by direct effects that are mediated by one-electron oxidation and indirect effects that are mediated by the reaction of water radiolysis products, e.g., hydroxyl radicals (•OH). In cellular DNA, direct and indirect effects appear to have about an equal effect toward DNA damage. We have shown that ϒ-(gamma) ray irradiation of aqueous solutions of DNA, during which •OH is the major damaging ROS can lead to the formation several lesions. On the other hand, the methylation and oxidative demethylation of cytosine in CpG dinucleotides plays a critical role in the gene regulation. The C5 position of cytosine in CG dinucleotides is frequently methylated by DNA methyl transferees (DNMTs) and constitutes 4-5% of the total cytosine. Here, my PhD research work focuses on the analysis of oxidative base modifications of model compounds of methylated and non methylated oligonucleotides, isolated DNA (calf-thymus DNA) and F98 cultured cell by gamma radiation. In addition, we identified a series of modifications of the 2-deoxyribose moiety of DNA arising from the exposure of isolated and cellular DNA to ionizing radiation. We also studied one electron oxidation of cellular DNA in cultured human HeLa cells initiated by intense nanosecond 266 nm laser pulse irradiation, which produces cross-links between guanine and thymine bases (G*-T*). To achieve these goals, we developed several methods based on mass spectrometry to analyze base modifications in isolated DNA and cellular DNA.

Análise do efeito de polimorfismos não-sinônimos em genes de reparo de DNA da via BER na resposta inflamatória da meningite

In vitro and in animal models, APE1, OGG1, and PARP-1 have been proposed as being involved with inflammatory response. In this work, we have investigated if the SNPs APE1 Asn148Glu, OGG1 Ser326Cys, and PARP-1 Val762Ala are associated to meningitis and also developed a system to enable the functional analysis of polymorphic proteins. Patients with bacterial meningitis (BM), aseptic meningitis (AM) and controls (non-infected) genotypes were investigated by PIRA-PCR or PCR-RFLP. DNA damages were detected in genomic DNA by Fpg treatment. IgG and IgA were measured from plasma and the cytokines and chemokines were measured from cerebrospinal fluid samples using Bio-Plex assays. The levels of NF-κB and c-Jun were measured in CSF by dot blot assays. A significant (P<0.05) increase in the frequency of APE1 148Glu allele in BM and AM patients was observed. A significant increase in the genotypes Asn/Asn in control group and Asn/Glu in BM group was also found. For the SNP OGG1 Ser326Cys, the genotype Cys/Cys was more frequent (P<0.05) in BM group. The frequency of PARP-1 Val/Val genotype was higher in control group (P<0.05). The occurrence of combined SNPs increased significantly in BM patients, indicating that these SNPs may be associated to the disease. Increasing in sensitive sites to Fpg was observed in carriers of APE1 148Glu allele or OGG1 326Cys allele, suggesting that SNPs affect DNA repair activity. Alterations in IgG production were observed in the presence of SNPs APE1Asn148Glu, OGG1Ser326Cys or PARP-1Val762Ala. Reductions in the levels ofIL-6, IL-1Ra, MCP-1/CCL2and IL-8/CXCL8 were observed in the presence of APE1148Glu allele in BM patients, however no differences were observed in the levels of NF-κB and c-Jun considering genotypes and analyzed groups. Using APE1 as model, a system to enable the analysis of cellular effects and functional characterization of polymorphic proteins was developed using strategies of cloning APE1 cDNA in pIRES2-EGFP vector, cellular transfection of the construction obtained, siRNA for endogenous APE1 and cellular cultures genotyping. In conclusion, we obtained evidences of an effect of SNPs in DNA repair genes on the regulation of immune response. This is a pioneering work in the field that shows association of BER variant enzymes with an infectious disease in human patients, suggesting that the SNPs analyzed may affect immune response and damage by oxidative stress level during brain infection. Considering these data, new approaches of functional characterization must be developed to better analysis and interactions of polymorphic proteins in response to this context

Dna Extraction From Filter-paper Spots Of Vaginal Samples Collected After Sexual Violence.

To detect the presence of male DNA in vaginal samples collected from survivors of sexual violence and stored on filter paper. A pilot study was conducted to evaluate 10 vaginal samples spotted on sterile filter paper: 6 collected at random in April 2009 and 4 in October 2010. Time between sexual assault and sample collection was 4-48hours. After drying at room temperature, the samples were placed in a sterile envelope and stored for 2-3years until processing. DNA extraction was confirmed by polymerase chain reaction for human β-globin, and the presence of prostate-specific antigen (PSA) was quantified. The presence of the Y chromosome was detected using primers for sequences in the TSPY (Y7/Y8 and DYS14) and SRY genes. β-Globin was detected in all 10 samples, while 2 samples were positive for PSA. Half of the samples amplified the Y7/Y8 and DYS14 sequences of the TSPY gene and 30% amplified the SRY gene sequence of the Y chromosome. Four male samples and 1 female sample served as controls. Filter-paper spots stored for periods of up to 3years proved adequate for preserving genetic material from vaginal samples collected following sexual violence.

Analysis Of The Dna Fourier Transform-infrared Microspectroscopic Signature Using An All-reflecting Objective.

The Fourier transform-infrared (FT-IR) signature of dry samples of DNA and DNA-polypeptide complexes, as studied by IR microspectroscopy using a diamond attenuated total reflection (ATR) objective, has revealed important discriminatory characteristics relative to the PO2(-) vibrational stretchings. However, DNA IR marks that provide information on the sample's richness in hydrogen bonds have not been resolved in the spectral profiles obtained with this objective. Here we investigated the performance of an all reflecting objective (ARO) for analysis of the FT-IR signal of hydrogen bonds in DNA samples differing in base richness types (salmon testis vs calf thymus). The results obtained using the ARO indicate prominent band peaks at the spectral region representative of the vibration of nitrogenous base hydrogen bonds and of NH and NH2 groups. The band areas at this spectral region differ in agreement with the DNA base richness type when using the ARO. A peak assigned to adenine was more evident in the AT-rich salmon DNA using either the ARO or the ATR objective. It is concluded that, for the discrimination of DNA IR hydrogen bond vibrations associated with varying base type proportions, the use of an ARO is recommended.

Prognostic Value Of Dna And Mrna E6/e7 Of Human Papillomavirus In The Evolution Of Cervical Intraepithelial Neoplasia Grade 2.

This study aimed at evaluating whether human papillomavirus (HPV) groups and E6/E7 mRNA of HPV 16, 18, 31, 33, and 45 are prognostic of cervical intraepithelial neoplasia (CIN) 2 outcome in women with a cervical smear showing a low-grade squamous intraepithelial lesion (LSIL). This cohort study included women with biopsy-confirmed CIN 2 who were followed up for 12 months, with cervical smear and colposcopy performed every three months. Women with a negative or low-risk HPV status showed 100% CIN 2 regression. The CIN 2 regression rates at the 12-month follow-up were 69.4% for women with alpha-9 HPV versus 91.7% for other HPV species or HPV-negative status (P < 0.05). For women with HPV 16, the CIN 2 regression rate at the 12-month follow-up was 61.4% versus 89.5% for other HPV types or HPV-negative status (P < 0.05). The CIN 2 regression rate was 68.3% for women who tested positive for HPV E6/E7 mRNA versus 82.0% for the negative results, but this difference was not statistically significant. The expectant management for women with biopsy-confirmed CIN 2 and previous cytological tests showing LSIL exhibited a very high rate of spontaneous regression. HPV 16 is associated with a higher CIN 2 progression rate than other HPV infections. HPV E6/E7 mRNA is not a prognostic marker of the CIN 2 clinical outcome, although this analysis cannot be considered conclusive. Given the small sample size, this study could be considered a pilot for future larger studies on the role of predictive markers of CIN 2 evolution.

Phylogenetics, Ancestral State Reconstruction, And A New Infrafamilial Classification Of The Pantropical Ochnaceae (medusagynaceae, Ochnaceae S.str., Quiinaceae) Based On Five Dna Regions.

Ochnaceae s.str. (Malpighiales) are a pantropical family of about 500 species and 27 genera of almost exclusively woody plants. Infrafamilial classification and relationships have been controversial partially due to the lack of a robust phylogenetic framework. Including all genera except Indosinia and Perissocarpa and DNA sequence data for five DNA regions (ITS, matK, ndhF, rbcL, trnL-F), we provide for the first time a nearly complete molecular phylogenetic analysis of Ochnaceae s.l. resolving most of the phylogenetic backbone of the family. Based on this, we present a new classification of Ochnaceae s.l., with Medusagynoideae and Quiinoideae included as subfamilies and the former subfamilies Ochnoideae and Sauvagesioideae recognized at the rank of tribe. Our data support a monophyletic Ochneae, but Sauvagesieae in the traditional circumscription is paraphyletic because Testulea emerges as sister to the rest of Ochnoideae, and the next clade shows Luxemburgia+Philacra as sister group to the remaining Ochnoideae. To avoid paraphyly, we classify Luxemburgieae and Testuleeae as new tribes. The African genus Lophira, which has switched between subfamilies (here tribes) in past classifications, emerges as sister to all other Ochneae. Thus, endosperm-free seeds and ovules with partly to completely united integuments (resulting in an apparently single integument) are characters that unite all members of that tribe. The relationships within its largest clade, Ochnineae (former Ochneae), are poorly resolved, but former Ochninae (Brackenridgea, Ochna) are polyphyletic. Within Sauvagesieae, the genus Sauvagesia in its broad circumscription is polyphyletic as Sauvagesia serrata is sister to a clade of Adenarake, Sauvagesia spp., and three other genera. Within Quiinoideae, in contrast to former phylogenetic hypotheses, Lacunaria and Touroulia form a clade that is sister to Quiina. Bayesian ancestral state reconstructions showed that zygomorphic flowers with adaptations to buzz-pollination (poricidal anthers), a syncarpous gynoecium (a near-apocarpous gynoecium evolved independently in Quiinoideae and Ochninae), numerous ovules, septicidal capsules, and winged seeds with endosperm are the ancestral condition in Ochnoideae. Although in some lineages poricidal anthers were lost secondarily, the evolution of poricidal superstructures secured the maintenance of buzz-pollination in some of these genera, indicating a strong selective pressure on keeping that specialized pollination system.

Diet Carotenoid Lutein Modulates The Expression Of Genes Related To Oxygen Transporters And Decreases Dna Damage And Oxidative Stress In Mice.

Lutein (LT) is a carotenoid obtained by diet and despite its antioxidant activity had been biochemically reported, few studies are available concerning its influence on the expression of antioxidant genes. The expression of 84 genes implicated in antioxidant defense was quantified using quantitative reverse transcription polymerase chain reaction array. DNA damage was measured by comet assay and glutathione (GSH) and thiobarbituric acid reactive substances (TBARS) were quantified as biochemical parameters of oxidative stress in mouse kidney and liver. cDDP treatment reduced concentration of GSH and increased TBARS, parameters that were ameliorated in treatment associated with LT. cDDP altered the expression of 32 genes, increasing the expression of GPx2, APC, Nqo1 and CCs. LT changed the expression of 37 genes with an induction of 13 mainly oxygen transporters. In treatments associating cDDP and LT, 30 genes had their expression changed with a increase of the same genes of the cDDP treatment alone. These results suggest that LT might act scavenging reactive species and also inducing the expression of genes related to a better antioxidant response, highlighting the improvement of oxygen transport. This improved redox state of the cell through LT treatment could be related to the antigenotoxic and antioxidant effects observed.

Changes In Liver Cell Dna Methylation Status In Diabetic Mice Affect Its Ft-ir Characteristics.

Lower levels of cytosine methylation have been found in the liver cell DNA from non-obese diabetic (NOD) mice under hyperglycemic conditions. Because the Fourier transform-infrared (FT-IR) profiles of dry DNA samples are differently affected by DNA base composition, single-stranded form and histone binding, it is expected that the methylation status in the DNA could also affect its FT-IR profile. The DNA FT-IR signatures obtained from the liver cell nuclei of hyperglycemic and normoglycemic NOD mice of the same age were compared. Dried DNA samples were examined in an IR microspectroscope equipped with an all-reflecting objective (ARO) and adequate software. Changes in DNA cytosine methylation levels induced by hyperglycemia in mouse liver cells produced changes in the respective DNA FT-IR profiles, revealing modifications to the vibrational intensities and frequencies of several chemical markers, including νas -CH3 stretching vibrations in the 5-methylcytosine methyl group. A smaller band area reflecting lower energy absorbed in the DNA was found in the hyperglycemic mice and assumed to be related to the lower levels of -CH3 groups. Other spectral differences were found at 1700-1500 cm(-1) and in the fingerprint region, and a slight change in the DNA conformation at the lower DNA methylation levels was suggested for the hyperglycemic mice. The changes that affect cytosine methylation levels certainly affect the DNA-protein interactions and, consequently, gene expression in liver cells from the hyperglycemic NOD mice.

A Modified Acidic Approach For Dna Extraction From Plant Species Containing High Levels Of Secondary Metabolites.

Purified genomic DNA can be difficult to obtain from some plant species because of the presence of impurities such as polysaccharides, which are often co-extracted with DNA. In this study, we developed a fast, simple, and low-cost protocol for extracting DNA from plants containing high levels of secondary metabolites. This protocol does not require the use of volatile toxic reagents such as mercaptoethanol, chloroform, or phenol and allows the extraction of high-quality DNA from wild and cultivated tropical species.

Human Herpesvirus-6 And Cytomegalovirus Dna In Liver Donor Biopsies And Their Correlation With Hla Matches And Acute Cellular Rejection.

Herpesvirus reactivation is common after liver transplantation. Analyze the presence of cytomegalovirus (HCMV) and human herpesvirus-6 (HHV-6) DNA in liver donor biopsies, seeking to better understand issues involving human donor leukocyte antigens (HLA)-A, B and DR, as well as correlations with acute cellular rejection. Fifty-nine liver transplantation patients were investigated for the presence of HCMV and HHV-6 DNA in liver donor biopsies, using the Nested-PCR technique. The clinical donor information and HLA matches were obtained from the São Paulo State Transplant System. The recipients' records regarding acute cellular rejection were studied. Seven (11.8%) biopsies were positive for HCMV DNA and 29 (49%) were positive for HHV-6 DNA. In 14 donors with HLA-DR 15 nine had HHV-6 DNA positive liver biopsy with a tendency for significant association (p=0.09), 22 recipients developed acute cellular rejection and 9/22 were positive for HLA-DR 15 (p=0.03; χ(2)=4.51), which was statistically significant in univariate analysis and showed a tendency after multivariate analysis (p=0.08). HHV-6 DNA was prevalent in liver donors studied as well as HLA-DR 15. These findings suggest that patients with HLA-DR 15 in liver donor biopsies develop more rejection after liver transplantation.

On The Consistency Of Monte Carlo Track Structure Dna Damage Simulations.

Monte Carlo track structures (MCTS) simulations have been recognized as useful tools for radiobiological modeling. However, the authors noticed several issues regarding the consistency of reported data. Therefore, in this work, they analyze the impact of various user defined parameters on simulated direct DNA damage yields. In addition, they draw attention to discrepancies in published literature in DNA strand break (SB) yields and selected methodologies. The MCTS code Geant4-DNA was used to compare radial dose profiles in a nanometer-scale region of interest (ROI) for photon sources of varying sizes and energies. Then, electron tracks of 0.28 keV-220 keV were superimposed on a geometric DNA model composed of 2.7 × 10(6) nucleosomes, and SBs were simulated according to four definitions based on energy deposits or energy transfers in DNA strand targets compared to a threshold energy ETH. The SB frequencies and complexities in nucleosomes as a function of incident electron energies were obtained. SBs were classified into higher order clusters such as single and double strand breaks (SSBs and DSBs) based on inter-SB distances and on the number of affected strands. Comparisons of different nonuniform dose distributions lacking charged particle equilibrium may lead to erroneous conclusions regarding the effect of energy on relative biological effectiveness. The energy transfer-based SB definitions give similar SB yields as the one based on energy deposit when ETH ≈ 10.79 eV, but deviate significantly for higher ETH values. Between 30 and 40 nucleosomes/Gy show at least one SB in the ROI. The number of nucleosomes that present a complex damage pattern of more than 2 SBs and the degree of complexity of the damage in these nucleosomes diminish as the incident electron energy increases. DNA damage classification into SSB and DSB is highly dependent on the definitions of these higher order structures and their implementations. The authors' show that, for the four studied models, different yields are expected by up to 54% for SSBs and by up to 32% for DSBs, as a function of the incident electrons energy and of the models being compared. MCTS simulations allow to compare direct DNA damage types and complexities induced by ionizing radiation. However, simulation results depend to a large degree on user-defined parameters, definitions, and algorithms such as: DNA model, dose distribution, SB definition, and the DNA damage clustering algorithm. These interdependencies should be well controlled during the simulations and explicitly reported when comparing results to experiments or calculations.