Numerous biological processes and mechanisms depend on details of base pairing and hydrogen bonding in DNA. Hydrogen bonds are challenging to quantify by X-ray crystallography and cryo-EM due to difficulty of visualizing hydrogen atom locations but can be probed with site specificity by NMR spectroscopy in solution and the solid state with the latter particularly suited to large, slowly tumbling DNA complexes. Recently, we showed that low-temperature dynamic nuclear polarization (DNP) enhanced solid-state NMR is a valuable tool for distinguishing Hoogsteen base pairs (bps) from canonical Watson-Crick bps in various DNA systems under native-like conditions. Here, using a model 12-mer DNA duplex containing two central adenine-thymine (A-T) bps in either Watson-Crick or Hoogsteen confirmation, we demonstrate DNP solid-state NMR measurements of thymine N3-H3 bond lengths, which are sensitive to details of N-H···N hydrogen bonding and permit hydrogen bonds for the two bp conformers to be systematically compared within the same DNA sequence context. For this DNA duplex, effectively identical TN3-H3 bond lengths of 1.055 ± 0.011 Å and 1.060 ± 0.011 Å were found for Watson-Crick A-T and Hoogsteen A (syn)-T base pairs, respectively, relative to a reference amide bond length of 1.015 ± 0.010 Å determined for N-acetyl-valine under comparable experimental conditions. Considering that prior quantum chemical calculations which account for zero-point motions predict a somewhat longer effective peptide N-H bond length of 1.041 Å, in agreement with solution and solid-state NMR studies of peptides and proteins at ambient temperature, to facilitate direct comparisons with these earlier studies TN3-H3 bond lengths for the DNA samples can be readily scaled appropriately to yield 1.083 Å and 1.087 Å for Watson-Crick A-T and Hoogsteen A (syn)-T bps, respectively, relative to the 1.041 Å reference peptide N-H bond length. Remarkably, in the context of the model DNA duplex, these results indicate that there are no significant differences in N-H···N A-T hydrogen bonds between Watson-Crick and Hoogsteen bp conformers. More generally, high precision measurements of N-H bond lengths by low-temperature DNP solid-state NMR based methods are expected to facilitate detailed comparative analysis of hydrogen bonding for a range of DNA complexes and base pairing environments.
- Award ID(s):
- 1715174
- NSF-PAR ID:
- 10374086
- Date Published:
- Journal Name:
- Proceedings of the National Academy of Sciences
- Volume:
- 119
- Issue:
- 30
- ISSN:
- 0027-8424
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
-
Genetic information is encoded in the DNA double helix, which, in its physiological milieu, is characterized by the iconical Watson-Crick nucleo-base pairing. Recent NMR relaxation experiments revealed the transient presence of an alternative, Hoogsteen (HG) base pairing pattern in naked DNA duplexes, and estimated its relative stability and lifetime. In contrast with DNA, such structures were not observed in RNA duplexes. Understanding HG base pairing is important because the underlying "breathing" motion between the two conformations can significantly modulate protein binding. However, a detailed mechanistic insight into the transition pathways and kinetics is still missing. We performed enhanced sampling simulation (with combined metadynamics and adaptive force-bias method) and Markov state modeling to obtain accurate free energy, kinetics, and the intermediates in the transition pathway between Watson-Crick and HG base pairs for both naked B-DNA and A-RNA duplexes. The Markov state model constructed from our unbiased MD simulation data revealed previously unknown complex extrahelical intermediates in the seemingly simple process of base flipping in B-DNA. Extending our calculation to A-RNA, for which HG base pairing is not observed experimentally, resulted in relatively unstable, single-hydrogen-bonded, distorted Hoogsteen-like bases. Unlike B-DNA, the transition pathway primarily involved base paired and intrahelical intermediates with transition timescales much longer than that of B-DNA. The seemingly obvious flip-over reaction coordinate (i.e., the glycosidic torsion angle) is unable to resolve the intermediates. Instead, a multidimensional picture involving backbone dihedral angles and distance between hydrogen bond donor and acceptor atoms is required to gain insight into the molecular mechanism.more » « less
-
more » « less
Abstract Artificially Expanded Genetic Information Systems (AEGIS) add independently replicable unnatural nucleotide pairs to the natural G:C and A:T/U pairs found in native DNA, joining the unnatural pairs through alternative modes of hydrogen bonding. Whether and how AEGIS pairs are recognized and processed by multi-subunit cellular RNA polymerases (RNAPs) remains unknown. Here, we show that
E. coli RNAP selectively recognizes unnatural nucleobases in a six-letter expanded genetic system. High-resolution cryo-EM structures of three RNAP elongation complexes containing template-substrate UBPs reveal the shared principles behind the recognition of AEGIS and natural base pairs. In these structures, RNAPs are captured in an active state, poised to perform the chemistry step. At this point, the unnatural base pair adopts a Watson-Crick geometry, and the trigger loop is folded into an active conformation, indicating that the mechanistic principles underlying recognition and incorporation of natural base pairs also apply to AEGIS unnatural base pairs. These data validate the design philosophy of AEGIS unnatural basepairs. Further, we provide structural evidence supporting a long-standing hypothesis that pair mismatch during transcription occurs via tautomerization. Together, our work highlights the importance of Watson-Crick complementarity underlying the design principles of AEGIS base pair recognition. -
In this work hydrogen bonding in a diverse set of 36 unnatural and the three natural Watson Crick base pairs adenine (A)–thymine (T), adenine (A)–uracil (U) and guanine (G)–cytosine (C) was assessed utilizing local vibrational force constants derived from the local mode analysis, originally introduced by Konkoli and Cremer as a unique bond strength measure based on vibrational spectroscopy. The local mode analysis was complemented by the topological analysis of the electronic density and the natural bond orbital analysis. The most interesting findings of our study are that (i) hydrogen bonding in Watson Crick base pairs is not exceptionally strong and (ii) the N–H⋯N is the most favorable hydrogen bond in both unnatural and natural base pairs while O–H⋯N/O bonds are the less favorable in unnatural base pairs and not found at all in natural base pairs. In addition, the important role of non-classical C–H⋯N/O bonds for the stabilization of base pairs was revealed, especially the role of C–H⋯O bonds in Watson Crick base pairs. Hydrogen bonding in Watson Crick base pairs modeled in the DNA via a QM/MM approach showed that the DNA environment increases the strength of the central N–H⋯N bond and the C–H⋯O bonds, and at the same time decreases the strength of the N–H⋯O bond. However, the general trends observed in the gas phase calculations remain unchanged. The new methodology presented and tested in this work provides the bioengineering community with an efficient design tool to assess and predict the type and strength of hydrogen bonding in artificial base pairs.more » « less
-
more » « less
Abstract The rational design of nanoscopic DNA tiles has yielded highly ordered crystalline matter in 2D and 3D. The most well‐studied 3D tile is the DNA tensegrity triangle, which is known to self‐assemble into macroscopic crystals. However, contemporary rational design parameters for 3D DNA crystals nearly universally invoke integer numbers of DNA helical turns and Watson–Crick (WC) base pairs. In this study, 24‐bp edges are substituted into a previously 21‐bp (two helical turns of DNA) tensegrity triangle motif to explore whether such unconventional motif can self‐assemble into 3D crystals. The use of noncanonical base pairs in the sticky ends results in a cubic arrangement of tensegrity triangles with exceedingly high symmetry, assembling a lattice from winding helical axes and diamond‐like tessellation patterns. Reverting this motif to sticky ends with Watson–Crick pairs results in a trigonal hexagonal arrangement, replicating this diamond arrangement in a hexagonal context. These results showcase that the authors can generate unexpected, highly complex, pathways for materials design by testing modifications to 3D tiles without prior knowledge of the ensuing symmetry. This study expands the rational design toolbox for DNA nanotechnology; and it further illustrates the existence of yet‐unexplored arrangements of crystalline soft matter.
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1073/pnas.2200681119
