Utilizing Biological Experimental Data and Molecular Dynamics for the Classification of Mutational Hotspots Through Machine Learning

Overview

Journal Bioinform Adv

Publisher Oxford University Press

Specialty Biology

Date 2024 Sep 6

PMID 39239360

Authors

James G Davies

Georgina E Menzies

Affiliations

Soon will be listed here.

Abstract

Motivation: Benzo[]pyrene, a notorious DNA-damaging carcinogen, belongs to the family of polycyclic aromatic hydrocarbons commonly found in tobacco smoke. Surprisingly, nucleotide excision repair (NER) machinery exhibits inefficiency in recognizing specific bulky DNA adducts including Benzo[]pyrene Diol-Epoxide (BPDE), a Benzo[]pyrene metabolite. While sequence context is emerging as the leading factor linking the inadequate NER response to BPDE adducts, the precise structural attributes governing these disparities remain inadequately understood. We therefore combined the domains of molecular dynamics and machine learning to conduct a comprehensive assessment of helical distortion caused by BPDE-Guanine adducts in multiple gene contexts. Specifically, we implemented a dual approach involving a random forest classification-based analysis and subsequent feature selection to identify precise topological features that may distinguish adduct sites of variable repair capacity. Our models were trained using helical data extracted from duplexes representing both BPDE hotspot and nonhotspot sites within the gene, then applied to sites within , , and genes.

Results: We show our optimized model consistently achieved exceptional performance, with accuracy, precision, and f1 scores exceeding 91%. Our feature selection approach uncovered that discernible variance in regional base pair rotation played a pivotal role in informing the decisions of our model. Notably, these disparities were highly conserved among and duplexes and appeared to be influenced by the regional GC content. As such, our findings suggest that there are indeed conserved topological features distinguishing hotspots and nonhotpot sites, highlighting regional GC content as a potential biomarker for mutation.

Availability And Implementation: Code for comparing machine learning classifiers and evaluating their performance is available at https://github.com/jdavies24/ML-Classifier-Comparison, and code for analysing DNA structure with Curves+ and Canal using Random Forest is available at https://github.com/jdavies24/ML-classification-of-DNA-trajectories.

References

Feng Z, Hu W, Chen J, Pao A, Li H, Rom W . Preferential DNA damage and poor repair determine ras gene mutational hotspot in human cancer. J Natl Cancer Inst. 2002; 94(20):1527-36. DOI: 10.1093/jnci/94.20.1527. View

Yu D, Berlin J, Penning T, Field J . Reactive oxygen species generated by PAH o-quinones cause change-in-function mutations in p53. Chem Res Toxicol. 2002; 15(6):832-42. DOI: 10.1021/tx010177m. View

Yella V, Bhimsaria D, Ghoshdastidar D, Rodriguez-Martinez J, Ansari A, Bansal M . Flexibility and structure of flanking DNA impact transcription factor affinity for its core motif. Nucleic Acids Res. 2018; 46(22):11883-11897. PMC: 6294565. DOI: 10.1093/nar/gky1057. View

Hodgkinson A, Eyre-Walker A . Variation in the mutation rate across mammalian genomes. Nat Rev Genet. 2011; 12(11):756-66. DOI: 10.1038/nrg3098. View

Beal M, Gagne R, Williams A, Marchetti F, Yauk C . Characterizing Benzo[a]pyrene-induced lacZ mutation spectrum in transgenic mice using next-generation sequencing. BMC Genomics. 2015; 16:812. PMC: 4617527. DOI: 10.1186/s12864-015-2004-4. View

Li C, Golebiowski F, Onishi Y, Samara N, Sugasawa K, Yang W . Tripartite DNA Lesion Recognition and Verification by XPC, TFIIH, and XPA in Nucleotide Excision Repair. Mol Cell. 2015; 59(6):1025-34. PMC: 4617536. DOI: 10.1016/j.molcel.2015.08.012. View

Hamidi E, Hajeb P, Selamat J, Razis A . Polycyclic Aromatic Hydrocarbons (PAHs) and their Bioaccessibility in Meat: a Tool for Assessing Human Cancer Risk. Asian Pac J Cancer Prev. 2016; 17(1):15-23. DOI: 10.7314/apjcp.2016.17.1.15. View

Yoon J, Smith L, Feng Z, Tang M, Lee C, Pfeifer G . Methylated CpG dinucleotides are the preferential targets for G-to-T transversion mutations induced by benzo[a]pyrene diol epoxide in mammalian cells: similarities with the p53 mutation spectrum in smoking-associated lung cancers. Cancer Res. 2001; 61(19):7110-7. View

Toyooka S, Tsuda T, Gazdar A . The TP53 gene, tobacco exposure, and lung cancer. Hum Mutat. 2003; 21(3):229-39. DOI: 10.1002/humu.10177. View

10.

Derreumaux S, Chaoui M, Tevanian G, Fermandjian S . Impact of CpG methylation on structure, dynamics and solvation of cAMP DNA responsive element. Nucleic Acids Res. 2001; 29(11):2314-26. PMC: 55717. DOI: 10.1093/nar/29.11.2314. View

11.

DiGiovanna J, Kraemer K . Shining a light on xeroderma pigmentosum. J Invest Dermatol. 2012; 132(3 Pt 2):785-96. PMC: 3279615. DOI: 10.1038/jid.2011.426. View

12.

Bae T, Tomasini L, Mariani J, Zhou B, Roychowdhury T, Franjic D . Different mutational rates and mechanisms in human cells at pregastrulation and neurogenesis. Science. 2017; 359(6375):550-555. PMC: 6311130. DOI: 10.1126/science.aan8690. View

13.

Ellegren H, Smith N, Webster M . Mutation rate variation in the mammalian genome. Curr Opin Genet Dev. 2003; 13(6):562-8. DOI: 10.1016/j.gde.2003.10.008. View

14.

Muheim R, Buterin T, Colgate K, Kolbanovsij A, Geacintov N, Naegeli H . Modulation of human nucleotide excision repair by 5-methylcytosines. Biochemistry. 2003; 42(11):3247-54. DOI: 10.1021/bi0268504. View

15.

Volker M, Mone M, Karmakar P, van Hoffen A, Schul W, Vermeulen W . Sequential assembly of the nucleotide excision repair factors in vivo. Mol Cell. 2001; 8(1):213-24. DOI: 10.1016/s1097-2765(01)00281-7. View

16.

Gunz D, Hess M, Naegeli H . Recognition of DNA adducts by human nucleotide excision repair. Evidence for a thermodynamic probing mechanism. J Biol Chem. 1996; 271(41):25089-98. DOI: 10.1074/jbc.271.41.25089. View

17.

Donny-Clark K, Broyde S . Influence of local sequence context on damaged base conformation in human DNA polymerase iota: molecular dynamics studies of nucleotide incorporation opposite a benzo[a]pyrene-derived adenine lesion. Nucleic Acids Res. 2009; 37(21):7095-109. PMC: 2790882. DOI: 10.1093/nar/gkp745. View

18.

Forbes S, Beare D, Gunasekaran P, Leung K, Bindal N, Boutselakis H . COSMIC: exploring the world's knowledge of somatic mutations in human cancer. Nucleic Acids Res. 2014; 43(Database issue):D805-11. PMC: 4383913. DOI: 10.1093/nar/gku1075. View

19.

Hoogstraten D, Bergink S, Ng J, Verbiest V, Luijsterburg M, Geverts B . Versatile DNA damage detection by the global genome nucleotide excision repair protein XPC. J Cell Sci. 2008; 121(Pt 17):2850-9. DOI: 10.1242/jcs.031708. View

20.

Denissenko M, Pao A, Pfeifer G, Tang M . Slow repair of bulky DNA adducts along the nontranscribed strand of the human p53 gene may explain the strand bias of transversion mutations in cancers. Oncogene. 1998; 16(10):1241-7. DOI: 10.1038/sj.onc.1201647. View