Hoogsteen Research Articles

How stop codon pseudouridylation induces nonsense suppression

Nonsense mutation can be an inherited disorder or caused by genetic mutation under certain conditions such as stress or exposure. Nonsense suppression, can have a broad potential as a therapeutic approach to nonsense mutation. Pseudouridine is proposed to result in nonsense suppression if induced in one of the three nonsense codons, Amber, Ochre, or Opal. Understanding how pseudouridine in stop codon induces nonsense suppression is vital for preventing nonsense mutation. Using in-silico analysis with DFT, we have simulated different near and non-cognate codon-anticodon interactions with pseudouridylated and normal nonsense codons. The structural analysis show that pseudouridylation of the nonsense codons in the first position leads to the mRNA-tRNA helix major groove conformation and formation of unusual Hoogsteen base pair in the purine bases at the second or third position. This leads to stabilized base-pair interaction and π…π stacking for the pseudouridylated nonsense codon compared its normal nonsense codon counterpart.

Hydrogen‐Bonding Interactions of Methylated Adenine Derivatives

AbstractHydrogen bonding between nucleobases is a crucial noncovalent interaction for the life on Earth. Apart from Watson‐Crick binding, Hoogsteen pairing has been found in many structures of nucleic acids. Methylation of nucleic acids (NAs) is a post‐replication or post‐transcription mechanism that can modulate the structure and function of a NA without changing its sequence. Methylation of adenine at the N6 position to form m6A is one of the most important and common epigenetic markers. This paper describes an investigation of intermolecular H‐bonding interactions of methylated derivatives of adenine with its complementary partner, thymine. Adenine derivatives with (di)methylamino groups in positions 2 or 6 have been prepared and their interactions with a thymine derivative have been studied by NMR spectroscopy and DFT calculations. It has been found that Hoogsteen pairing is preferred for adenine derivatives, which offer two hydrogen‐bond sites on both Watson‐Crick and Hoogsteen sides of the molecule. Methylation of the N6 position leads to further stabilization of the Hoogsteen pair.

Probing the Hydrogen-Bonding Environment of Individual Bases in DNA Duplexes with Isotope-Edited Infrared Spectroscopy.

Measuring the strength of the hydrogen bonds between DNA base pairs is of vital importance for understanding how our genetic code is physically accessed and recognized in cells, particularly during replication and transcription. Therefore, it is important to develop probes for these key hydrogen bonds (H-bonds) that dictate events critical to cellular function, such as the localized melting of DNA. The vibrations of carbonyl bonds are well-known probes of their H-bonding environment, and their signals can be observed with infrared (IR) spectroscopy. Yet, pinpointing a single bond of interest in the complex IR spectrum of DNA is challenging due to the large number of carbonyl signals that overlap with each other. Here, we develop a method using isotope editing and infrared (IR) spectroscopy to isolate IR signals from the thymine (T) C2=O carbonyl. We use solvatochromatic studies to show that the TC2=O signal’s position in the IR spectrum is sensitive to the H-bonding capacity of the solvent. Our results indicate that C2=O of a single T base within DNA duplexes experiences weak H-bonding interactions. This finding is consistent with the existence of a third, noncanonical CH···O H-bond between adenine and thymine in both Watson–Crick and Hoogsteen base pairs in DNA.

Probing the Conformational States of a pH-Sensitive DNA Origami Zipper via Label-Free Electrochemical Methods.

DNA origami structures represent an exciting class of materials for use in a wide range of biotechnological applications. This study reports the design, production, and characterization of a DNA origami “zipper” structure, which contains nine pH-responsive DNA locks. Each lock consists of two parts that are attached to the zipper’s opposite arms: a DNA hairpin and a single-stranded DNA that are able to form a DNA triplex through Hoogsteen base pairing. The sequences of the locks were selected in a way that the zipper adopted a closed configuration at pH 6.5 and an open state at pH 8.0 (transition pKa 7.6). By adding thiol groups, it was possible to immobilize the zipper structure onto gold surfaces. The immobilization process was characterized electrochemically to confirm successful adsorption of the zipper. The open and closed states were then probed using differential pulse voltammetry and electrochemical impedance spectroscopy with solution-based redox agents. It was found that after immobilization, the open or closed state of the zipper in different pH regimes could be determined by electrochemical interrogation. These findings pave the way for development of DNA origami-based pH monitoring and other pH-responsive sensing and release strategies for zipper-functionalized gold surfaces.

DNA polymerase λ promotes error-free replication through Watson–Crick impairing N1-methyl-deoxyadenosine adduct in conjunction with DNA polymerase ζ

In a previous study, we showed that replication through the N1-methyl-deoxyadenosine (1-MeA) adduct in human cells is mediated via three different Polι/Polθ, Polη, and Polζ-dependent pathways. Based on biochemical studies with these Pols, in the Polι/Polθ pathway, we inferred a role for Polι in the insertion of a nucleotide (nt) opposite 1-MeA and of Polθ in extension of synthesis from the inserted nt; in the Polη pathway, we inferred that this Pol alone would replicate through 1-MeA; in the Polζ pathway, however, the Pol required for inserting an nt opposite 1-MeA had remained unidentified. In this study, we provide biochemical and genetic evidence for a role for Polλ in inserting the correct nt T opposite 1-MeA, from which Polζ would extend synthesis. The high proficiency of purified Polλ for inserting a T opposite 1-MeA implicates a role for Polλ—which normally uses W-C base pairing for DNA synthesis—in accommodating 1-MeA in a syn confirmation and forming a Hoogsteen base pair with T. The potential of Polλ to replicate through DNA lesions by Hoogsteen base pairing adds another novel aspect to Polλ’s role in translesion synthesis in addition to its role as a scaffolding component of Polζ. We discuss how the action mechanisms of Polλ and Polζ could be restrained to inserting a T opposite 1-MeA and extending synthesis thereafter, respectively.

A unified computational view of DNA duplex, triplex, quadruplex and their donor-acceptor interactions.

DNA can assume various structures as a result of interactions at atomic and molecular levels (e.g., hydrogen bonds, π–π stacking interactions, and electrostatic potentials), so understanding of the consequences of these interactions could guide development of ways to produce elaborate programmable DNA for applications in bio- and nanotechnology. We conducted advanced ab initio calculations to investigate nucleobase model structures by componentizing their donor-acceptor interactions. By unifying computational conditions, we compared the independent interactions of DNA duplexes, triplexes, and quadruplexes, which led us to evaluate a stability trend among Watson–Crick and Hoogsteen base pairing, stacking, and even ion binding. For a realistic solution-like environment, the influence of water molecules was carefully considered, and the potassium-ion preference of G-quadruplex was first analyzed at an ab initio level by considering both base-base and ion-water interactions. We devised new structure factors including hydrogen bond length, glycosidic vector angle, and twist angle, which were highly effective for comparison between computationally-predicted and experimentally-determined structures; we clarified the function of phosphate backbone during nucleobase ordering. The simulated tendency of net interaction energies agreed well with that of real world, and this agreement validates the potential of ab initio study to guide programming of complicated DNA constructs.

Molecular Dynamical Investigation of the 7,8-dihydro-8-oxoguanine Mutation in dsDNA

Background: The oxidization of a Guanine (G) base pair to 7,8-dihydro-8-oxoguanine (OG) is one of the most common DNA mutations. OG mutations can undergo a regular Watson-Crick base-pairing, or a reverse Hoogsteen (HG) base-pairing, especially in OG:A mismatches. While the causes of these mutations are well-understood, the kinetic and energetic characteristics of this new pseudo-base have never been fully investigated, especially at temperatures around biological function (17-37°C). Methods: We created a simulation to derive the Free Energy Surface (FES) of OG:C and OG:A Hoogsteen to Watson-Crick base-pair (bp) transitions under multiple temperatures, relative to 2 collective geometric variables: the dihedral Chi and the pseudo-dihedral CPD angle. To make the simulation, we used the relatively recent Metadynamics algorithms in conjunction with GROMACS 2020.2. Results: The lowest free energy increased linearly with increasing temperatures (17-37°C). Major Chi and CPD rotations at these minima varied heavily for 27°C and 32°C (the largest was seen in the former), but stayed relatively similar for other temperatures, indicating a highly sensitive relationship to temperature, likely due to DNA flexibility, quantum mechanical (QM) effects, and hydrogen bonding. Free energies had a weak negative linear relationship, and free energy hypersurfaces were given for studied temperatures of 17-37°C. Human body temperature (37°C) results were also included and explained. The simulations showed why OG:A Hoogsteen bps often occur in organisms and are energetically preferable to standard Watson-Crick. OG:C HG base pairings are determined to likely be not as common as OG:A HG. Limitations: Future investigations must focus on discovering rate constants of these base-pairs, as time constraints did not permit them to be done here, as well as more QM-focused simulations.

Effects of Adenine Methylation on the Structure and Thermodynamic Stability of a DNA Minidumbbell.

DNA methylation is a prevalent regulatory modification in prokaryotes and eukaryotes. N1-methyladenine (m1A) and N6-methyladenine (m6A) have been found to be capable of altering DNA structures via disturbing Watson–Crick base pairing. However, little has been known about their influences on non-B DNA structures, which are associated with genetic instabilities. In this work, we investigated the effects of m1A and m6A on both the structure and thermodynamic stability of a newly reported DNA minidumbbell formed by two TTTA tetranucleotide repeats. As revealed by the results of nuclear magnetic resonance spectroscopic studies, both m1A and m6A favored the formation of a T·m1A and T·m6A Hoogsteen base pair, respectively. More intriguingly, the m1A and m6A modifications brought about stabilization and destabilization effects on the DNA minidumbbell, respectively. This work provides new biophysical insights into the effects of adenine methylation on the structure and thermodynamic stability of DNA.

The biofilm matrix scaffold of Pseudomonas aeruginosa contains G-quadruplex extracellular DNA structures

Extracellular DNA, or eDNA, is recognised as a critical biofilm component; however, it is not understood how it forms networked matrix structures. Here, we isolate eDNA from static-culture Pseudomonas aeruginosa biofilms using ionic liquids to preserve its biophysical signatures of fluid viscoelasticity and the temperature dependency of DNA transitions. We describe a loss of eDNA network structure as resulting from a change in nucleic acid conformation, and propose that its ability to form viscoelastic structures is key to its role in building biofilm matrices. Solid-state analysis of isolated eDNA, as a proxy for eDNA structure in biofilms, reveals non-canonical Hoogsteen base pairs, triads or tetrads involving thymine or uracil, and guanine, suggesting that the eDNA forms G-quadruplex structures. These are less abundant in chromosomal DNA and disappear when eDNA undergoes conformation transition. We verify the occurrence of G-quadruplex structures in the extracellular matrix of intact static and flow-cell biofilms of P. aeruginosa, as displayed by the matrix to G-quadruplex-specific antibody binding, and validate the loss of G-quadruplex structures in vivo to occur coincident with the disappearance of eDNA fibres. Given their stability, understanding how extracellular G-quadruplex structures form will elucidate how P. aeruginosa eDNA builds viscoelastic networks, which are a foundational biofilm property.

Free-Energy Landscape of a pH-Modulated G·C Base Pair Transition from Watson-Crick to Hoogsteen State in Duplex DNA.

In double-helical DNAs, the most stable Watson-Crick (WC) base pair (bp) can be in thermal equilibrium with much less abundant Hoogsteen (HG) bp by the spontaneous rotation of the glycosidic angle in purine bases. Previous experimental studies showed that in the case of a G·C bp, the population of the transient HG is enhanced as a protonated form (HG+) through the protonation of the cytosine base under weakly acidic conditions. Hence, pH is a key factor that can modulate this transition event from the WC to HG+ bp. In this study, to computationally probe the overall free-energy landscapes of this pH-modulated G·C HG breathing, a comprehensive classical molecular dynamics (MD) simulation protocol is proposed using an enhanced sampling MD in conjunction with the standard thermodynamic integration method. From this MD protocol proposed, the free-energy surfaces of the G·C bp transition from the WC to HG bp were constructed successfully at any pH range, producing pH-dependent free-energy quantities in close agreement with previously reported experimental results. The simulation protocol is expected to provide valuable atomistic insight into the DNA bp transition events coupled with protonation or tautomeric shift in a target bp.

Discovering Collective Variables of Molecular Transitions via Genetic Algorithms and Neural Networks.

With the continual improvement of computing hardware and algorithms, simulations have become a powerful tool for understanding all sorts of (bio)molecular processes. To handle the large simulation data sets and to accelerate slow, activated transitions, a condensed set of descriptors, or collective variables (CVs), is needed to discern the relevant dynamics that describes the molecular process of interest. However, proposing an adequate set of CVs that can capture the intrinsic reaction coordinate of the molecular transition is often extremely difficult. Here, we present a framework to find an optimal set of CVs from a pool of candidates using a combination of artificial neural networks and genetic algorithms. The approach effectively replaces the encoder of an autoencoder network with genes to represent the latent space, i.e., the CVs. Given a selection of CVs as input, the network is trained to recover the atom coordinates underlying the CV values at points along the transition. The network performance is used as an estimator of the fitness of the input CVs. Two genetic algorithms optimize the CV selection and the neural network architecture. The successful retrieval of optimal CVs by this framework is illustrated at the hand of two case studies: the well-known conformational change in the alanine dipeptide molecule and the more intricate transition of a base pair in B-DNA from the classic Watson–Crick pairing to the alternative Hoogsteen pairing. Key advantages of our framework include the following: optimal interpretable CVs, avoiding costly calculation of committor or time-correlation functions, and automatic hyperparameter optimization. In addition, we show that applying a time-delay between the network input and output allows for enhanced selection of slow variables. Moreover, the network can also be used to generate molecular configurations of unexplored microstates, for example, for augmentation of the simulation data.

Watson-Crick versus Hoogsteen Base Pairs: Chemical Strategy to Encode and Express Genetic Information in Life.

Nucleic acids typically form a double helix structure through Watson-Crick base-pairing. In contrast, non-Watson-Crick base pairs can form other three-dimensional structures. Although it is well-known that Watson-Crick base pairs may be more unstable than non-Watson-Crick base pairs under some conditions, the importance of non-Watson-Crick base pairs has not been widely examined. Hoogsteen base pairs, the non-Watson-Crick base pairs, contain important hydrogen-bond patterns that form the helices of nucleic acids, such as in Watson-Crick base pairs, and can form non-double helix structures such as triplexes and quadruplexes. In recent years, non-double helix structures have been discovered in cells and were reported to considerably influence gene expression. The complex behavior of these nucleic acids in cells is gradually being revealed, but the underlying mechanisms remain almost unknown.Quantitatively analyzing the structural stability of nucleic acids is important for understanding their behavior. A nucleic acid is an anionic biopolymer composed of a sugar, base, and phosphoric acid. The physicochemical factors that determine the stability of nucleic acid structures include those derived from the interactions of nucleic acid structures and those derived from the environments surrounding nucleic acids. The Gibbs free energy change (ΔG) of structure formation is the most commonly used physicochemical parameter for analyzing quantitative stability. Quantitatively understanding the intracellular behavior of nucleic acids involves describing the formation of nucleic acid structures and related reactions as ΔG. Based on this concept, we quantitatively analyzed the stability of double helix and non-double helix structures and found that decreased water activity, an important factor in crowded cellular conditions, significantly destabilize the formation of Watson-Crick base pairs but stabilizes Hoogsteen base pairs.Here, we describe a physicochemical approach to understand the regulation of gene expressions based on the stability of nucleic acid structures. We developed new methods for predicting the stability of double and non-double helices in various molecular environments by mimicking intracellular environments. Furthermore, the physicochemical approach used for analyzing gene expression regulated by non-double helix structures is useful for not only determining how gene expression is controlled by cellular environments but also for developing new technologies to chemically regulate gene expression by targeting non-double helix structures. We discuss the roles of Watson-Crick and Hoogsteen base pairs in cells based on our results and why both types of base pairing are required for life. Finally, a new concept in nucleic acid science beyond that of Watson and Crick base pairing is introduced.

Hoogsteen Base Pairs Increase the Susceptibility of Double-Stranded DNA to Cytotoxic Damage

Watson-Crick base pairing is believed to protect the Watson-Crick face of nucleobases in double-stranded DNA (dsDNA) from cytotoxic damage. Lesions on the Watson-Crick face are believed to occur when the DNA is transiently single stranded during replication and transcription. However, damage to the Watson-Crick face in dsDNA has been reported in vitro, but the mechanisms rendering dsDNA susceptible to damage remain poorly understood. in this work, we explored the role of Hoogsteen base pairs (bps) in mediating cytotoxic damage to dsDNA.

RNA architecture influences plant biology

The majority of the genome is transcribed to RNA in living organisms. RNA transcripts can form astonishing arrays of secondary and tertiary structures via Watson-Crick, Hoogsteen, or wobble base pairing. In vivo, RNA folding is not a simple thermodynamic event of minimizing free energy. Instead, the process is constrained by transcription, RNA-binding proteins, steric factors, and the microenvironment. RNA secondary structure (RSS) plays myriad roles in numerous biological processes, such as RNA processing, stability, transportation, and translation in prokaryotes and eukaryotes. Emerging evidence has also implicated RSS in RNA trafficking, liquid-liquid phase separation, and plant responses to environmental variations such as temperature and salinity. At molecular level, RSS is correlated with splicing, polyadenylation, protein synthesis, and miRNA biogenesis and functions. In this review, we summarize newly reported methods for probing RSS in vivo and functions and mechanisms of RSS in plant physiology.

The emerging structural complexity of G-quadruplex RNAs.

G-quadruplexes (G4s) are four-stranded nucleic acid structures that arise from the stacking of G-quartets, cyclic arrangements of four guanines engaged in Hoogsteen base-pairing. Until recently, most RNA G4 structures were thought to conform to a sequence pattern in which guanines stacking within the G4 would also be contiguous in sequence (e.g., four successive guanine trinucleotide tracts separated by loop nucleotides). Such a sequence restriction, and the stereochemical constraints inherent to RNA (arising, in particular, from the presence of the 2′-OH), dictate relatively simple RNA G4 structures. Recent crystallographic and solution NMR structure determinations of a number of in vitro selected RNA aptamers have revealed RNA G4 structures of unprecedented complexity. Structures of the Sc1 aptamer that binds an RGG peptide from the Fragile-X mental retardation protein, various fluorescence turn-on aptamers (Corn, Mango, and Spinach), and the spiegelmer that binds the complement protein C5a, in particular, reveal complexity hitherto unsuspected in RNA G4s, including nucleotides in syn conformation, locally inverted strand polarity, and nucleotide quartets that are not all-G. Common to these new structures, the sequences folding into G4s do not conform to the requirement that guanine stacks arise from consecutive (contiguous in sequence) nucleotides. This review highlights how emancipation from this constraint drastically expands the structural possibilities of RNA G-quadruplexes.

Experimental and theoretical assessment of protonated Hoogsteen 9-methylguanine-1-methylcytosine base-pair dissociation: kinetics within a statistical reaction framework.

We investigated the collision-induced dissociation (CID) reactions of a protonated Hoogsteen 9-methylguanine-1-methylcytosine base pair (HG-[9MG·1MC + H]+), which aims to address the mystery of the literature reported "anomaly" in product ion distributions and compare the kinetics of a Hoogsteen base pair with its Watson-Crick isomer WC-[9MG·1MC + H]+ (reported recently by Sun et al.; Phys. Chem. Chem. Phys., 2020, 22, 24986). Product ion cross sections and branching ratios were measured as a function of center-of-mass collision energy using guided-ion beam tandem mass spectrometry, from which base-pair dissociation energies were determined. Product structures and energetics were assessed using various theories, of which the composite DLPNO-CCSD(T)/aug-cc-pVTZ//ωB97XD/6-311++G(d,p) was adopted as the best-performing method for constructing a reaction potential energy surface. The statistical Rice-Ramsperger-Kassel-Marcus theory was found to provide a useful framework for rationalizing the dominating abundance of [1MC + H]+ over [9MG + H]+ in the fragment ions of HG-[9MG·1MC + H]+. The kinetics analysis proved the necessity for incorporating into kinetics modeling not only the static properties of reaction minima and transition states but more importantly, the kinetics of individual base-pair conformers that have formed in collisional activation. The analysis also pinpointed the origin of the statistical kinetics of HG-[9MG·1MC + H]+vs. the non-statistical behavior of WC-[9MG·1MC + H]+ in terms of their distinctively different intra-base-pair hydrogen-bonds and consequently the absence of proton transfer between the N1 position of 9MG and the N3' of 1MC in the Hoogsteen base pair. Finally, the Hoogsteen base pair was examined in the presence of a water ligand, i.e., HG-[9MG·1MC + H]+·H2O. Besides the same type of base-pair dissociation as detected in dry HG-[9MG·1MC + H]+, secondary methanol elimination was observed via the SN2 reaction of water with nucleobase methyl groups.

Novel mechanisms of the conformational transformations of the biologically important G·C nucleobase pairs in Watson-Crick, Hoogsteen and wobble configurations via the mutual rotations of the bases around the intermolecular H-bonds: a QM/QTAIM study.

At the MP2/6-311++G(2df,pd)//B3LYP/6-311++G(d,p) level of quantum-mechanical theory, we provide for the first time a comprehensive investigation of the physico-chemical mechanisms of the 55 conformational transformations of the biologically-important G·C nucleobase pairs – Watson–Crick (WC), reverse Watson–Crick (rWC), Hoogsteen (H), reverse Hoogsteen (rH), wobble (w) and reverse wobble (rw) base pairs by the participation of the G and C bases in the canonical and rare tautomeric forms (“r” – means reverse configuration of the base pair). It was established that all these G·C nucleobase pairs can conformationally transform into each other without the changing of the tautomeric status of the G and C bases. These transitions occur through significantly non-planar transition states via the mutual rotation of the G and C bases relative to each other within the G·C nucleobase pair around the upper, middle or lower intermolecular H-bonds: WC ↔ rWC, WC ↔ rwWC, rWC ↔ WC, rWC ↔ wWC, wWC ↔ rwWC, H ↔ rH, H ↔ rwH, rH ↔ H, rH ↔ wH, wH ↔ rwH. Gibbs free energies ΔG of activation for these conformational transformations are ΔG = 2.96–19.04/3.58–13.36 kcal mol−1 (in vacuum under normal conditions (T = 298.15 K)), which means that these reactions proceed quite fast. Obtained conformational transformations are accompanied by the disruption and further formation of the intermolecular specific contacts in the G·C nucleobase pairs (H-bonds and attractive van der Waals contacts). As a result, 76 conformers of the G·C nucleobase pairs were established – 48 base pairs in WC, rWC, wWC and rwWC configurations and 28 base pairs in H, rH, wH and rwH configurations with relative Gibbs free ΔG/electronic ΔE energies in the range ΔG/ΔE = 0.00–44.73/0.00–46.99 and ΔG/ΔE = 0.00–37.52/0.00–38.54 kcal mol−1, respectively (in vacuum under normal conditions). Experimental investigation and verification of the novel G·C nucleobase pairs are promising tasks for the future research. Based on the obtained data, biologically important conclusions were made about the importance of the conformational mobility of the G·C nucleobase pairs for the understanding of the functioning of the DNA and RNA molecules and their transition from the parallel into the anti-parallel duplexes and vice versa.

Slow motions in A\xb7T rich DNA sequence

In free B-DNA, slow (microsecond-to-millisecond) motions that involve equilibrium between Watson–Crick (WC) and Hoogsteen (HG) base-pairing expand the DNA dynamic repertoire that could mediate DNA–protein assemblies. R1ρ relaxation dispersion NMR methods are powerful tools to capture such slow conformational exchanges in solution using 13C/15 N labelled DNA. Here, these approaches were applied to a dodecamer containing a TTAAA element that was assumed to facilitate nucleosome formation. NMR data and inferred exchange parameters assign HG base pairs as the minor, transient conformers specifically observed in three successive A·T base pairs forming the TAA·TTA segment. The abundance of these HG A·T base pairs can be up to 1.2% which is high compared to what has previously been observed. Data analyses support a scenario in which the three adenines undergo non-simultaneous motions despite their spatial proximity, thus optimising the probability of having one HG base pair in the TAA·TTA segment. Finally, revisiting previous NMR data on H2 resonance linewidths on the basis of our results promotes the idea of there being a special propensity of A·T base pairs in TAA·TTA tracts to adopt HG pairing. In summary, this study provides an example of a DNA functional element submitted to slow conformational exchange. More generally, it strengthens the importance of the role of the DNA sequence in modulating its dynamics, over a nano- to milli-second time scale.

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.

P.026 Role of orexin-1 receptor within medial prefrontal cortex in acquisition and expression of morphine- and natural induced conditioned place preference in rats

The ends of eukaryotic chromosomes typically contain a 3′ ssDNA G-rich protrusion (G-overhang). This overhang must be protected against detrimental activities of nucleases and of the DNA damage response machinery and participates in the regulation of telomerase, a ribonucleoprotein complex that maintains telomere integrity. These functions are mediated by DNA-binding proteins, such as Cdc13 in Saccharomyces cerevisiae, and the propensity of G-rich sequences to form various non-B DNA structures. Using CD and NMR spectroscopies, we show here that G-overhangs of S. cerevisiae form distinct Hoogsteen pairing–based secondary structures, depending on their length. Whereas short telomeric oligonucleotides form a G-hairpin, their longer counterparts form parallel and/or antiparallel G-quadruplexes (G4s). Regardless of their topologies, non-B DNA structures exhibited impaired binding to Cdc13 in vitro as demonstrated by electrophoretic mobility shift assays. Importantly, whereas G4 structures formed relatively quickly, G-hairpins folded extremely slowly, indicating that short G-overhangs, which are typical for most of the cell cycle, are present predominantly as single-stranded oligonucleotides and are suitable substrates for Cdc13. Using ChIP, we show that the occurrence of G4 structures peaks at the late S phase, thus correlating with the accumulation of long G-overhangs. We present a model of how time- and length-dependent formation of non-B DNA structures at chromosomal termini participates in telomere maintenance.