Ranking of Molecular Biomarker Interaction with Targeted DNA Nucleobases via Full Atomistic Molecular Dynamics

DNA-based sensors can detect disease biomarkers, including acetone and ethanol for diabetes and H2S for cardiovascular diseases. Before experimenting on thousands of potential DNA segments, we conduct full atomistic steered molecular dynamics (SMD) simulations to screen the interactions between different DNA sequences with targeted molecules to rank the nucleobase sensing performance. We study and rank the strength of interaction between four single DNA nucleotides (Adenine (A), Guanine (G), Cytosine (C), and Thymine (T)) on single-stranded DNA (ssDNA) and double-stranded DNA (dsDNA) with acetone, ethanol, H2S and HCl. By sampling forward and reverse interaction paths, we compute the free-energy profiles of eight systems for the four targeted molecules. We find that dsDNA react differently than ssDNA to the targeted molecules, requiring more energy to move the molecule close to DNA as indicated by the potential of mean force (PMF). Comparing the PMF values of different systems, we obtain a relative ranking of DNA base for the detection of each molecule. Via the same procedure, we could generate a library of DNA sequences for the detection of a wide range of chemicals. A DNA sensor array built with selected sequences differentiating many disease biomarkers can be used in disease diagnosis and monitoring.