Relation between DNA ionization potentials, single base substitutions and pathogenic variants

The Influence of Spirodi(Iminohydantoin) on Charge Transfer through ds-DNA Containing 8-OXO-dG: A Theoretical Approach

Genetic information stored in a DNA base sequence is continuously exposed to harmful factors. It has been determined that 9 × 10⁴ different DNA damage events occur in a single human cell every 24 h. Of these, 7,8-dihydro-8-oxo-guanosine (^OXOG) is one of the most abundant and can undergo further transformations towards spirodi(iminohydantoin) (Sp). Sp is highly mutagenic in comparison to its precursor if not repaired. In this paper, the influence of both Sp diastereomers 4R and 4S as well as their anti and syn conformers on charge transfer through the double helix was taken into theoretical consideration. In addition, the electronic properties of four modelled double-stranded oligonucleotides (ds-oligos) were also discussed, i.e., d[A₁Sp₂A₃^oxoG₄A₅] * [T₅C₄T₃C₂T₁]. Throughout the study, the M06-2X/6-31++G** level theory was used. Solvent-solute non-equilibrated and equilibrated interactions were also considered. The subsequent results elucidated that the 7,8-dihydro-8-oxo-guanosine:cytidine (^OXOGC) base pair is the settled point of a migrated radical cation in each of the discussed cases, due to its low adiabatic ionization potential, i.e., _~5.55 [eV]. The opposite was noted for excess electron transfer through ds-oligos containing anti (R)-Sp or anti (S)-Sp. The radical anion was found on the ^OXOGC moiety, whereas in the presence of syn (S)-Sp or syn (R)-Sp, an excess electron was found on the distal A₁T₅ or A₅T₁ base pair, respectively. Furthermore, a spatial geometry analysis of the discussed ds-oligos revealed that the presence of syn (R)-Sp in the ds-oligo caused only a slight deformation to the double helix, while syn (S)-Sp formed an almost ideal base pair with a complementary dC. The above results are in strong agreement with the final charge transfer rate constant, as calculated according to Marcus' theory. In conclusion, DNA damage such as spirodi(iminohydantoin), especially when becoming part of clustered DNA damage, can affect the effectiveness of other lesion recognition and repair processes. This can lead to the acceleration of undesired and deleterious processes such as carcinogenesis or aging. However, in terms of anticancer radio-/chemo- or combined therapy, the slowing down of the repair machinery can result in increased effectiveness. With this in mind, the influence of clustered damage on charge transfer and its subsequent effect on single-damage recognition by glycosylases justifies future investigation.

Estimating the Vertical Ionization Potential of Single-Stranded DNA Molecules

The electronic properties of DNA molecules, defined by the sequence-dependent ionization potentials of nucleobases, enable long-range charge transport along the DNA stacks. This has been linked to a range of key physiological processes in the cells and to the triggering of nucleobase substitutions, some of which may cause diseases. To gain molecular-level understanding of the sequence dependence of these phenomena, we estimated the vertical ionization potential (vIP) of all possible nucleobase stacks in B-conformation, containing one to four Gua, Ade, Thy, Cyt, or methylated Cyt. To do this, we used quantum chemistry calculations and more precisely the second-order Møller-Plesset perturbation theory (MP2) and three double-hybrid density functional theory methods, combined with several basis sets for describing atomic orbitals. The calculated vIP of single nucleobases were compared to experimental data and those of nucleobase pairs, triplets, and quadruplets, to observed mutability frequencies in the human genome, reported to be correlated with vIP values. This comparison selected MP2 with the 6-31G* basis set as the best of the tested calculation levels. These results were exploited to set up a recursive model, called vIPer, which estimates the vIP of all possible single-stranded DNA sequences of any length based on the calculated vIPs of overlapping quadruplets. vIPer's vIP values correlate well with oxidation potentials measured by cyclic voltammetry and activities obtained through photoinduced DNA cleavage experiments, further validating our approach. vIPer is freely available on the github.com/3BioCompBio/vIPer repository.

Single-strand mismatch and damage patterns revealed by single-molecule DNA sequencing

Mutations accumulate in the genome of every cell of the body throughout life, causing cancer and other genetic diseases1-4. Almost all of these mosaic mutations begin as nucleotide mismatches or damage in only one of the two strands of the DNA prior to becoming double-strand mutations if unrepaired or misrepaired5. However, current DNA sequencing technologies cannot resolve these initial single-strand events. Here, we developed a single-molecule, long-read sequencing method that achieves single-molecule fidelity for single-base substitutions when present in either one or both strands of the DNA. It also detects single-strand cytosine deamination events, a common type of DNA damage. We profiled 110 samples from diverse tissues, including from individuals with cancer-predisposition syndromes, and define the first single-strand mismatch and damage signatures. We find correspondences between these single-strand signatures and known double-strand mutational signatures, which resolves the identity of the initiating lesions. Tumors deficient in both mismatch repair and replicative polymerase proofreading show distinct single-strand mismatch patterns compared to samples deficient in only polymerase proofreading. In the mitochondrial genome, our findings support a mutagenic mechanism occurring primarily during replication. Since the double-strand DNA mutations interrogated by prior studies are only the endpoint of the mutation process, our approach to detect the initiating single-strand events at single-molecule resolution will enable new studies of how mutations arise in a variety of contexts, especially in cancer and aging.

Large-scale in silico mutagenesis experiments reveal optimization of genetic code and codon usage for protein mutational robustness

BACKGROUND

How, and the extent to which, evolution acts on DNA and protein sequences to ensure mutational robustness and evolvability is a long-standing open question in the field of molecular evolution. We addressed this issue through the first structurome-scale computational investigation, in which we estimated the change in folding free energy upon all possible single-site mutations introduced in more than 20,000 protein structures, as well as through available experimental stability and fitness data.

RESULTS

At the amino acid level, we found the protein surface to be more robust against random mutations than the core, this difference being stronger for small proteins. The destabilizing and neutral mutations are more numerous in the core and on the surface, respectively, whereas the stabilizing mutations are about 4% in both regions. At the genetic code level, we observed smallest destabilization for mutations that are due to substitutions of base III in the codon, followed by base I, bases I+III, base II, and other multiple base substitutions. This ranking highly anticorrelates with the codon-anticodon mispairing frequency in the translation process. This suggests that the standard genetic code is optimized to limit the impact of random mutations, but even more so to limit translation errors. At the codon level, both the codon usage and the usage bias appear to optimize mutational robustness and translation accuracy, especially for surface residues.

CONCLUSION

Our results highlight the non-universality of mutational robustness and its multiscale dependence on protein features, the structure of the genetic code, and the codon usage. Our analyses and approach are strongly supported by available experimental mutagenesis data.