Characterizing Watson-Crick versus Hoogsteen Base Pairing in a DNA-Protein Complex Using Nuclear Magnetic Resonance and Site-Specifically 13C- and 15N-Labeled DNA

Probing Watson-Crick and Hoogsteen base pairing in duplex DNA using dynamic nuclear polarization solid-state NMR spectroscopy

The majority of base pairs in double-stranded DNA exist in the canonical Watson-Crick geometry. However, they can also adopt alternate Hoogsteen conformations in various complexes of DNA with proteins and small molecules, which are key for biological function and mechanism. While detection of Hoogsteen base pairs in large DNA complexes and assemblies poses considerable challenges for traditional structural biology techniques, we show here that multidimensional dynamic nuclear polarization-enhanced solid-state NMR can serve as a unique spectroscopic tool for observing and distinguishing Watson-Crick and Hoogsteen base pairs in a broad range of DNA systems based on characteristic NMR chemical shifts and internuclear dipolar couplings. We illustrate this approach using a model 12-mer DNA duplex, free and in complex with the antibiotic echinomycin, which features two central adenine-thymine base pairs with Watson-Crick and Hoogsteen geometry, respectively, and subsequently extend it to the ∼200 kDa Widom 601 DNA nucleosome core particle.

1H NMR Chemical Exchange Techniques Reveal Local and Global Effects of Oxidized Cytosine Derivatives

5-Carboxycytosine (5caC) is a rare epigenetic modification found in nucleic acids of all domains of life. Despite its sparse genomic abundance, 5caC is presumed to play essential regulatory roles in transcription, maintenance and base-excision processes in DNA. In this work, we utilize nuclear magnetic resonance (NMR) spectroscopy to address the effects of 5caC incorporation into canonical DNA strands at multiple pH and temperature conditions. Our results demonstrate that 5caC has a pH-dependent global destabilizing and a base-pair mobility enhancing local impact on dsDNA, albeit without any detectable influence on the ground-state B-DNA structure. Measurement of hybridization thermodynamics and kinetics of 5caC-bearing DNA duplexes highlighted how acidic environment (pH 5.8 and 4.7) destabilizes the double-stranded structure by ∼10–20 kJ mol^–1 at 37 °C when compared to the same sample at neutral pH. Protonation of 5caC results in a lower activation energy for the dissociation process and a higher barrier for annealing. Studies on conformational exchange on the microsecond time scale regime revealed a sharply localized base-pair motion involving exclusively the modified site and its immediate surroundings. By direct comparison with canonical and 5-formylcytosine (5fC)-edited strands, we were able to address the impact of the two most oxidized naturally occurring cytosine derivatives in the genome. These insights on 5caC’s subtle sensitivity to acidic pH contribute to the long-standing questions of its capacity as a substrate in base excision repair processes and its purpose as an independent, stable epigenetic mark.

Revealing A-T and G-C Hoogsteen base pairs in stressed protein-bound duplex DNA

Watson–Crick base pairs (bps) are the fundamental unit of genetic information and the building blocks of the DNA double helix. However, A-T and G-C can also form alternative ‘Hoogsteen’ bps, expanding the functional complexity of DNA. We developed ‘Hoog-finder’, which uses structural fingerprints to rapidly screen Hoogsteen bps, which may have been mismodeled as Watson–Crick in crystal structures of protein–DNA complexes. We uncovered 17 Hoogsteen bps, 7 of which were in complex with 6 proteins never before shown to bind Hoogsteen bps. The Hoogsteen bps occur near mismatches, nicks and lesions and some appear to participate in recognition and damage repair. Our results suggest a potentially broad role for Hoogsteen bps in stressed regions of the genome and call for a community-wide effort to identify these bps in current and future crystal structures of DNA and its complexes.

Z-form extracellular DNA is a structural component of the bacterial biofilm matrix

Biofilms are community architectures adopted by bacteria inclusive of a self-formed extracellular matrix that protects resident bacteria from diverse environmental stresses and, in many species, incorporates extracellular DNA (eDNA) and DNABII proteins for structural integrity throughout biofilm development. Here, we present evidence that this eDNA-based architecture relies on the rare Z-form. Z-form DNA accumulates as biofilms mature and, through stabilization by the DNABII proteins, confers structural integrity to the biofilm matrix. Indeed, substances known to drive B-DNA into Z-DNA promoted biofilm formation whereas those that drive Z-DNA into B-DNA disrupted extant biofilms. Importantly, we demonstrated that the universal bacterial DNABII family of proteins stabilizes both bacterial- and host-eDNA in the Z-form in situ. A model is proposed that incorporates the role of Z-DNA in biofilm pathogenesis, innate immune response, and immune evasion.

Rapid assessment of Watson-Crick to Hoogsteen exchange in unlabeled DNA duplexes using high-power SELOPE imino 1H CEST

In duplex DNA, Watson-Crick A-T and G-C base pairs (bp's) exist in dynamic equilibrium with an alternative Hoogsteen conformation, which is low in abundance and short-lived. Measuring how the Hoogsteen dynamics varies across different DNA sequences, structural contexts and physiological conditions is key for identifying potential Hoogsteen hot spots and for understanding the potential roles of Hoogsteen base pairs in DNA recognition and repair. However, such studies are hampered by the need to prepare 13 C or 15 N isotopically enriched DNA samples for NMR relaxation dispersion (RD) experiments. Here, using SELective Optimized Proton Experiments (SELOPE) 1 H CEST experiments employing high-power radiofrequency fields (B 1  >  250 Hz) targeting imino protons, we demonstrate accurate and robust characterization of Watson-Crick to Hoogsteen exchange, without the need for isotopic enrichment of the DNA. For 13 residues in three DNA duplexes under different temperature and pH conditions, the exchange parameters deduced from high-power imino 1 H CEST were in very good agreement with counterparts measured using off-resonance 13 C /  15 N spin relaxation in the rotating frame (R 1 ρ ). It is shown that 1 H-1 H NOE effects which typically introduce artifacts in 1 H-based measurements of chemical exchange can be effectively suppressed by selective excitation, provided that the relaxation delay is short (≤  100 ms). The 1 H CEST experiment can be performed with ∼  10× higher throughput and ∼  100× lower cost relative to 13 C /  15 N R 1 ρ and enabled Hoogsteen chemical exchange measurements undetectable by R 1 ρ . The results reveal an increased propensity to form Hoogsteen bp's near terminal ends and a diminished propensity within A-tract motifs. The 1 H CEST experiment provides a basis for rapidly screening Hoogsteen breathing in duplex DNA, enabling identification of unusual motifs for more in-depth characterization.

Contingency and chance erase necessity in the experimental evolution of ancestral proteins

The roles of chance, contingency, and necessity in evolution are unresolved because they have never been assessed in a single system or on timescales relevant to historical evolution. We combined ancestral protein reconstruction and a new continuous evolution technology to mutate and select proteins in the B-cell lymphoma-2 (BCL-2) family to acquire protein–protein interaction specificities that occurred during animal evolution. By replicating evolutionary trajectories from multiple ancestral proteins, we found that contingency generated over long historical timescales steadily erased necessity and overwhelmed chance as the primary cause of acquired sequence variation; trajectories launched from phylogenetically distant proteins yielded virtually no common mutations, even under strong and identical selection pressures. Chance arose because many sets of mutations could alter specificity at any timepoint; contingency arose because historical substitutions changed these sets. Our results suggest that patterns of variation in BCL-2 sequences – and likely other proteins, too – are idiosyncratic products of a particular and unpredictable course of historical events.

One of the most fundamental and unresolved questions in evolutionary biology is whether the outcomes of evolution are predictable. Is the diversity of life we see today the expected result of organisms adapting to their environment throughout history (also known as natural selection) or the product of random chance? Or did chance events early in history shape the paths that evolution could take next, determining the biological forms that emerged under natural selection much later?

These questions are hard to study because evolution happened only once, long ago. To overcome this barrier, Xie, Pu, Metzger et al. developed an experimental approach that can evolve reconstructed ancestral proteins that existed deep in the past. Using this method, it is possible to replay evolution multiple times, from various historical starting points, under conditions similar to those that existed long ago. The end products of the evolutionary trajectories can then be compared to determine how predictable evolution actually is.

Xie, Pu, Metzger et al. studied proteins belonging to the BCL-2 family, which originated some 800 million years ago. These proteins have diversified greatly over time in both their genetic sequences and their ability to bind to specific partner proteins called co-regulators. Xie, Pu, Metzger et al. synthesized BCL-2 proteins that existed at various times in the past. Each ancestral protein was then allowed to evolve repeatedly under natural selection to acquire the same co-regulator binding functions that evolved during history.

At the end of each evolutionary trajectory, the genetic sequence of the resulting BCL-2 proteins was recorded. This revealed that the outcomes of evolution were almost completely unpredictable: trajectories initiated from the same ancestral protein produced proteins with very different sequences, and proteins launched from different ancestral starting points were even more dissimilar.

Further experiments identified the mutations in each trajectory that caused changes in coregulator binding. When these mutations were introduced into other ancestral proteins, they did not yield the same change in function. This suggests that early chance events influenced each protein’s evolution in an unpredictable way by opening and closing the paths available to it in the future.

This research expands our understanding of evolution on a molecular level whilst providing a new experimental approach for studying evolutionary drivers in more detail. The results suggest that BCL-2 proteins, in all their various forms, are unique products of a particular, unpredictable course of history set in motion by ancient chance events.

Using quantum chemistry to estimate chemical shifts in biomolecules

An automated fragmentation quantum mechanics/molecular mechanics approach (AFNMR) has shown promising results in chemical shift calculations for biomolecules. Sample results for ubiquitin, and an RNA hairpin and helix are presented, and used to recent directions in quantum calculations. Trends in chemical shift are stable with regards to change in density functional or basis sets, and the use of the small "pcSseg-0" basis, which was optimized for chemical shift prediction [1], opens the way to more extensive conformational averaging, which can often be necessary, even for fairly well-defined structures.

Hoogsteen base pairs increase the susceptibility of double-stranded DNA to cytotoxic damage

As the Watson-Crick faces of nucleobases are protected in dsDNA, it is commonly assumed that deleterious alkylation damage to the Watson-Crick faces of nucleobases predominantly occurs when DNA becomes single-stranded during replication and transcription. However, damage to the Watson-Crick faces of nucleobases has been reported in dsDNA in vitro through mechanisms that are not understood. In addition, the extent of protection from methylation damage conferred by dsDNA relative to ssDNA has not been quantified. Watson-Crick base pairs in dsDNA exist in dynamic equilibrium with Hoogsteen base pairs that expose the Watson-Crick faces of purine nucleobases to solvent. Whether this can influence the damage susceptibility of dsDNA remains unknown. Using dot-blot and primer extension assays, we measured the susceptibility of adenine-N1 to methylation by dimethyl sulfate (DMS) when in an A-T Watson-Crick versus Hoogsteen conformation. Relative to unpaired adenines in a bulge, Watson-Crick A-T base pairs in dsDNA only conferred ∼130-fold protection against adenine-N1 methylation, and this protection was reduced to ∼40-fold for A(syn)-T Hoogsteen base pairs embedded in a DNA-drug complex. Our results indicate that Watson-Crick faces of nucleobases are accessible to alkylating agents in canonical dsDNA and that Hoogsteen base pairs increase this accessibility. Given the higher abundance of dsDNA relative to ssDNA, these results suggest that dsDNA could be a substantial source of cytotoxic damage. The work establishes DMS probing as a method for characterizing A(syn)-T Hoogsteen base pairs in vitro and also lays the foundation for a sequencing approach to map A(syn)-T Hoogsteen and unpaired adenines genome-wide in vivo.

Infrared Spectroscopic Observation of a G-C+ Hoogsteen Base Pair in the DNA:TATA-Box Binding Protein Complex Under Solution Conditions

Hoogsteen DNA base pairs (bps) are an alternative base pairing to canonical Watson-Crick bps and are thought to play important biochemical roles. Hoogsteen bps have been reported in a handful of X-ray structures of protein-DNA complexes. However, there are several examples of Hoogsteen bps in crystal structures that form Watson-Crick bps when examined under solution conditions. Furthermore, Hoogsteen bps can sometimes be difficult to resolve in DNA:protein complexes by X-ray crystallography due to ambiguous electron density and by solution-state NMR spectroscopy due to size limitations. Here, using infrared spectroscopy, we report the first direct solution-state observation of a Hoogsteen (G-C⁺ ) bp in a DNA:protein complex under solution conditions with specific application to DNA-bound TATA-box binding protein. These results support a previous assignment of a G-C⁺ Hoogsteen bp in the complex, and indicate that Hoogsteen bps do indeed exist under solution conditions in DNA:protein complexes.