Published online 24 April 2021 Nucleic Acids Research, 2021, Vol. 49, No. 9 4919–4933 doi: 10.1093/nar/gkab285 A unified computational view of DNA duplex, triplex, quadruplex and their donor–acceptor interactions Gyuri Park1, Byunghwa Kang1, Soyeon V. Park1, Donghwa Lee1,2,3,* and Seung Soo Oh 1,3,4,*
1Department of Materials Science and Engineering, Pohang University of Science and Technology (POSTECH), Pohang 37673, South Korea, 2Division of Advanced Materials Science, Pohang University of Science and Technology (POSTECH), Pohang 37673, South Korea, 3Institute for Convergence Research and Education in Advanced Technology (I-CREATE), Yonsei University, Incheon 21983, South Korea and 4School of Interdisciplinary Bioscience and Bioengineering, Pohang University of Science and Technology (POSTECH), Pohang 37673, South Korea
Received October 15, 2020; Revised April 07, 2021; Editorial Decision April 07, 2021; Accepted April 14, 2021
ABSTRACT INTRODUCTION DNA can assume various structures as a result Deoxyribonucleic acid (DNA) is a unique material as com- of interactions at atomic and molecular levels posed of nitrogenous bases (adenine (A), thymine (T), (e.g., hydrogen bonds, – stacking interactions, guanine (G) or cytosine (C)), sugar rings, and phosphate and electrostatic potentials), so understanding of groups. DNA is programmable, so it can be rationally the consequences of these interactions could guide designed into molecular structures ranging from simple Watson–Crick base-pairing primers to DNA origami-based development of ways to produce elaborate pro- complex 3D constructs (1,2). Even spatial and temporal grammable DNA for applications in bio- and nan- control of DNA nanostructures is achievable; sophisticated otechnology. We conducted advanced ab initio cal- DNA molecular machines can perform a series of nanome- culations to investigate nucleobase model struc- chanical motions in a controllable manner, so this abil- tures by componentizing their donor-acceptor inter- ity provides unprecedented applications in bio- and nan- actions. By unifying computational conditions, we otechnology (3–6). Such active use of DNA requires funda- compared the independent interactions of DNA du- mental understanding of DNA folding and its stabilization. plexes, triplexes, and quadruplexes, which led us In particular, DNA hybridization programming extensively to evaluate a stability trend among Watson–Crick exploits understanding of various donor-acceptor interac- and Hoogsteen base pairing, stacking, and even tions of DNA, including hydrogen bonds, – stacking in- ion binding. For a realistic solution-like environ- teractions, and electrostatic potentials (7). Hydrogen bonds and – stacking are among the most ment, the influence of water molecules was care- important intra- and inter-molecular interactions in DNA fully considered, and the potassium-ion preference (8–11). In general, nucleobases create specific hydrogen of G-quadruplex was first analyzed at an ab initio bonds between a purine (A or G) and a pyrimidine (C and level by considering both base-base and ion-water T), which yield planar Watson–Crick base pairs (A–T; G– interactions. We devised new structure factors in- C). Between adjacent base pairs, – stacking interactions cluding hydrogen bond length, glycosidic vector an- result in sequential stacking of base pairs. Therefore, hy- gle, and twist angle, which were highly effective for drogen bonds and – stacking drive the formation of a comparison between computationally-predicted and duplex as a basic structure of DNA. However, the nucle- experimentally-determined structures; we clarified obases can be influenced by different interactions due to pH the function of phosphate backbone during nucle- and metal ions, so triplexes and quadruplexes are sometimes obase ordering. The simulated tendency of net inter- produced (12,13). For instance, C is protonated at slightly acidic pH to become a hydrogen-bond acceptor C+, which action energies agreed well with that of real world, binds to the guanine of a G–C pair as a bond donor to and this agreement validates the potential of ab ini- yield a C+•G–C triad in a DNA triplex (14). The C+ can tio study to guide programming of complicated DNA be also paired with a non-protonated form of C, and the re- constructs. sulting C–C+ pairs involve in formation of C-quadruplex, i.e., i-motif (15). Moreover, guanines can interact strongly
*To whom correspondence should be addressed. Tel: +82 54 279 2144; Email: [email protected] Correspondence may also be addressed to Donghwa Lee. Tel: +82 54 279 2160; Email: [email protected]