Flanking Base Pairs Research Articles

JASPAR 2024: 20th anniversary of the open-access database of transcription factor binding profiles.

JASPAR (https://jaspar.elixir.no/) is a widely-used open-access database presenting manually curated high-quality and non-redundant DNA-binding profiles for transcription factors (TFs) across taxa. In this 10th release and 20th-anniversary update, the CORE collection has expanded with 329 new profiles. We updated three existing profiles and provided orthogonal support for 72 profiles from the previous release's UNVALIDATED collection. Altogether, the JASPAR 2024 update provides a 20% increase in CORE profiles from the previous release. A trimming algorithm enhanced profiles by removing low information content flanking base pairs, which were likely uninformative (within the capacity of the PFM models) for TFBS predictions and modelling TF-DNA interactions. This release includes enhanced metadata, featuring a refined classification for plant TFs' structural DNA-binding domains. The new JASPAR collections prompt updates to the genomic tracks of predicted TFbinding sites (TFBSs) in 8 organisms, with human and mouse tracks available as native tracks in the UCSC Genome browser. All data are available through the JASPAR web interface and programmatically through its API and the updated Bioconductor and pyJASPAR packages. Finally, a new TFBS extraction tool enables users to retrieve predicted JASPAR TFBSs intersecting their genomic regions of interest.

Interaction of an anticancer benzopyrane derivative with DNA: Biophysical, biochemical, and molecular modeling studies

BackgroundSIMR1281 is a potent anticancer lead candidate with multi- target activity against several proteins; however, its mechanism of action at the molecular level is not fully understood. Revealing the mechanism and the origin of multitarget activity is important for the rational identification and optimization of multitarget drugs. MethodsWe have used a variety of biophysical (circular dichroism, isothermal titration calorimetry, viscosity, and UV DNA melting), biochemical (topoisomerase I & II assays) and computational (molecular docking and MD simulations) methods to study the interaction of SIMR1281 with duplex DNA structures. ResultsThe biophysical results revealed that SIMR1281 binds to dsDNA via an intercalation-binding mode with an average binding constant of 3.1 × 106 M−1. This binding mode was confirmed by the topoisomerases' inhibition assays and molecular modeling simulations, which showed the intercalation of the benzopyrane moiety between DNA base pairs, while the remaining moieties (thiazole and phenyl rings) sit in the minor groove and interact with the flanking base pairs adjacent to the intercalation site. ConclusionsThe DNA binding characteristics of SIMR1281, which can disrupt/inhibit DNA function as confirmed by the topoisomerases' inhibition assays, indicate that the observed multi-target activity might originate from ligand intervention at nucleic acids level rather than due to direct interactions with multiple biological targets at the protein level. General significanceThe findings of this study could be helpful to guide future optimization of benzopyrane-based ligands for therapeutic purposes.

Revisiting DNA Sequence-Dependent Deformability in High-Resolution Structures: Effects of Flanking Base Pairs on Dinucleotide Morphology and Global Chain Configuration.

DNA carries more than the list of biochemical ingredients that drive the basic functions of living systems. The sequence of base pairs includes a multitude of structural and energetic signals, which determine the degree to which the long, threadlike molecule moves and how it responds to proteins and other molecules that control its processing and govern its packaging. The chemical composition of base pairs directs the spatial disposition and fluctuations of successive residues. The observed arrangements of these moieties in high-resolution protein–DNA crystal structures provide one of the best available estimates of the natural, sequence-dependent structure and deformability of the double-helical molecule. Here, we update the set of knowledge-based elastic potentials designed to describe the observed equilibrium structures and configurational fluctuations of the ten unique base-pair steps. The large number of currently available structures makes it possible to characterize the configurational preferences of the DNA base-pair steps within the context of their immediate neighbors, i.e., tetrameric context. Use of these knowledge-based potentials shows promise in accounting for known effects of sequence in long chain molecules, e.g., the degree of curvature reported in classic gel mobility studies and the recently reported sequence-dependent responses of supercoiled minicircles to nuclease cleavage.

2-Amino-1,3-benzothiazole-6-carboxamide Preferentially Binds the Tandem Mismatch Motif r(UY:GA).

RNA helices are often punctuated with non-Watson-Crick features that may be targeted by chemical compounds, but progress toward identifying such compounds has been slow. We embedded a tandem UU:GA mismatch motif (5'-UG-3':5'-AU-3') within an RNA hairpin stem to identify compounds that bind the motif specifically. The three-dimensional structure of the RNA hairpin and its interaction with a small molecule identified through virtual screening are presented. The G-A mismatch forms a sheared pair upon which the U-U base pair stacks. The hydrogen bond configuration of the U-U pair involves O2 of the U adjacent to the G and O4 of the U adjacent to the A. The G-A and U-U pairs are flanked by A-U and G-C base pairs, respectively, and the stability of the mismatch is greater than when the motif is within the context of other flanking base pairs or when the 5'-3' orientation of the G-A and U-U pairs is swapped. Residual dipolar coupling constants were used to generate an ensemble of structures against which a virtual screen of 64480 small molecules was performed. The tandem mismatch was found to be specific for one compound, 2-amino-1,3-benzothiazole-6-carboxamide, which binds with moderate affinity but extends the motif to include the flanking A-U and G-C base pairs. The finding that the affinity for the UU:GA mismatch is dependent on flanking sequence emphasizes the importance of the motif context and potentially increases the number of small noncanonical features within RNA that can be specifically targeted by small molecules.

Melting temperature measurement and mesoscopic evaluation of single, double and triple DNA mismatches.

Unlike the canonical base pairs AT and GC, the molecular properties of mismatches such as hydrogen bonding and stacking interactions are strongly dependent on the identity of the neighbouring base pairs. As a result, due to the sheer number of possible combinations of mismatches and flanking base pairs, only a fraction of these have been studied in varying experiments or theoretical models. Here, we report on the melting temperature measurement and mesoscopic analysis of contiguous DNA mismatches in nearest-neighbours and next-nearest neighbour contexts. A total of 4032 different mismatch combinations, including single, double and triple mismatches were covered. These were compared with 64 sequences containing all combinations of canonical base pairs in the same location under the same conditions. For a substantial number of single mismatch configurations, 15%, the measured melting temperatures were higher than the least stable AT base pair. The mesoscopic calculation, using the Peyrard–Bishop model, was performed on the set of 4096 sequences, and resulted in estimates of on-site and nearest-neighbour interactions that can be correlated to hydrogen bonding and base stacking. Our results confirm many of the known properties of mismatches, including the peculiar sheared stacking of tandem GA mismatches. More intriguingly, it also reveals that a number of mismatches present strong hydrogen bonding when flanked on both sites by other mismatches. To highlight the applicability of our results, we discuss a number of practical situations such as enzyme binding affinities, thymine DNA glycosylase repair activity, and trinucleotide repeat expansions.

Molecular Dynamics Study of One-Base Deletion Duplexes Containing the Major DNA Adduct Formed by Ochratoxin A: Effects of Sequence Context and Adduct Ionization State on Lesion Site Structure and Mutagenicity.

Ochratoxin A (OTA) is a ubiquitous food toxin associated with chronic nephropathy in humans and renal carcinogenicity in rodents. The mutational spectra of cells exposed to OTA reveal that one-base deletions comprise the largest percentage (73%) of the total mutations that occur upon OTA exposure. To contribute toward understanding the prevalence of OTA-induced one-base deletion mutations, the present work uses molecular dynamics (MD) simulations to analyze the conformational preferences of one-base deletion duplexes containing OT-G, the major OTA adduct (addition product) at the C8-site of guanine. Specifically, the influence of OT-G in four possible ionization states and three sequence contexts (G1, G2 and G3 in the NarI (5'-G1G2CG3CC-3'), a prokaryotic mutational hotspot sequence) on the structure of the adducted DNA is investigated. Our data reveal that the damaged helices are stable in two (B-type (B) and stacked (S)) conformations that are structurally similar to those adopted by common N-linked C8-guanine lesions. However, the adduct ionization state and sequence context affect the degree of helical distortion and the B/S conformational heterogeneity, which will impact the lesion repair and replication outcomes. This finding correlates with the experimentally reported tissue-specific mutagenicity of OTA exposure. Furthermore, regardless of the adduct conformation, ionization state, or sequence context, more stable lesion-site interactions and lack of disruption of the flanking base pairs in the one-base deletion duplexes compared to the corresponding two-base deletion helices rationalize the greater abundance of OTA induced one-base deletions. Overall, our work provides valuable structural insights that help explain the experimentally observed mutagenicity associated with OTA.

Characterization of mutational processes in ER+ metastatic breast cancer.

1019 Background: Diverse mutational processes are active at different stages of cancer development and cause the mutational burden in cancer cells. These processes leave distinctive signatures, characterized by varying frequencies of nucleotide substitutions and their flanking base-pairs. Mutational signatures found in early stage breast cancer have been described, but less is known about the mutational processes in metastatic breast cancer (MBC). Methods: We performed whole exome sequencing (WES) on 209 metastatic tumor biopsies from patients with ER+ MBC. In a subset of patients, we obtained and performed WES on the matched primary biopsy (n = 62). Decomposition of mutational signatures was done by determining the linear combination of 30 previously reported mutational signatures (Alexandrov et al., 2013) that most accurately reconstructs the mutational profile of each tumor sample. For comparison, signature decomposition was also performed on 695 treatment-naïve primary ER+ breast cancers form The Cancer Genome Atlas (TCGA). Results: Cohort-level comparison of signature activity between our cohort of MBC biopsies versus the primary ER+ breast cancers from TCGA demonstrated enrichment of Signature 13 in MBC samples (q = 0.002), attributed to activity of the APOBEC family of cytidine deaminases and implicated in resistance to endocrine therapies. Activity of Signature 22, associated to activity of the nucleotide-excision repair mechanism, was depleted in MBC tumors relative to TCGA ER+ breast tumors (q = 0.02). Decomposition of mutational signatures at the clone level in the subset of patients with both metastatic and matched primary samples allowed us to infer signatures acquired in the metastatic setting. Among mutations in acquired clones, we found enrichment of APOBEC signatures (Sig. 2, q = 0.018; Sig. 13, q = 0.003). Conversely, activity of Signature 1, attributed to the normal aging process, was enriched in mutations shared between primary and MBC samples (q = 0.004). This suggests that mutations developed early in tumor evolution may be driven primarily by natural processes, while later mutations may reflect selective pressures imposed by therapies and cellular defects acquired during disease progression. Conclusions: The profile of somatic mutations in MBC differs from early stage breast cancer. APOBEC-related mutations arising late during tumor evolution are one possible mechanism for this difference. Further studying the associations between clinical features and therapies to mutational processes could inform design of clinical trials and management of MBC.

Signals of Variation in Human Mutation Rate at Multiple Levels of Sequence Context

Our understanding of the human mutation rate helps us build evolutionary models and interpret patterns of genetic variation observed in human populations. Recent work indicates that the frequencies of specific polymorphism types have been elevated in Europe, and that many more, subtler signatures of global polymorphism variation may yet remain unidentified. Here, we present an analysis of the 1000 Genomes Project supported by analysis in the Simons Genome Diversity Panel, suggesting additional putative signatures of mutation rate variation across populations and the extent to which they are shaped by local sequence context. First, we compiled a list of the most significantly variable polymorphism types in a cross-continental statistical test. Clustering polymorphisms together, we observe three sets that showed distinct shared patterns of relative enrichment among ancestral populations, and we characterize each one of these putative “signatures” of polymorphism variation. For three of these signatures, we found that a single flanking base pair of sequence context was sufficient to determine the majority of enrichment or depletion of a polymorphism type. However, local genetic context up to 2–3 bp away contributes additional variability and may help to interpret a previously noted enrichment of certain polymorphism types in some East Asian groups. Moreover, considering broader local genetic context highlights patterns of polymorphism variation, which were not captured by previous approaches. Building our understanding of mutation rate in this way can help us to construct more accurate evolutionary models and better understand the mechanisms that underlie genetic change.

Structural insights into the mechanism of double strand break formation by Hermes, a hAT family eukaryotic DNA transposase.

Some DNA transposons relocate from one genomic location to another using a mechanism that involves generating double-strand breaks at their transposon ends by forming hairpins on flanking DNA. The same double-strand break mode is employed by the V(D)J recombinase at signal-end/coding-end junctions during the generation of antibody diversity. How flanking hairpins are formed during DNA transposition has remained elusive. Here, we describe several co-crystal structures of the Hermes transposase bound to DNA that mimics the reaction step immediately prior to hairpin formation. Our results reveal a large DNA conformational change between the initial cleavage step and subsequent hairpin formation that changes which strand is acted upon by a single active site. We observed that two factors affect the conformational change: the complement of divalent metal ions bound by the catalytically essential DDE residues, and the identity of the –2 flanking base pair. Our data also provides a mechanistic link between the efficiency of hairpin formation (an A:T basepair is favored at the –2 position) and Hermes' strong target site preference. Furthermore, we have established that the histidine residue within a conserved C/DxxH motif present in many transposase families interacts directly with the scissile phosphate, suggesting a crucial role in catalysis.

BpRNA: large-scale automated annotation and analysis of RNA secondary structure.

While RNA secondary structure prediction from sequence data has made remarkable progress, there is a need for improved strategies for annotating the features of RNA secondary structures. Here, we present bpRNA, a novel annotation tool capable of parsing RNA structures, including complex pseudoknot-containing RNAs, to yield an objective, precise, compact, unambiguous, easily-interpretable description of all loops, stems, and pseudoknots, along with the positions, sequence, and flanking base pairs of each such structural feature. We also introduce several new informative representations of RNA structure types to improve structure visualization and interpretation. We have further used bpRNA to generate a web-accessible meta-database, ‘bpRNA-1m’, of over 100 000 single-molecule, known secondary structures; this is both more fully and accurately annotated and over 20-times larger than existing databases. We use a subset of the database with highly similar (≥90% identical) sequences filtered out to report on statistical trends in sequence, flanking base pairs, and length. Both the bpRNA method and the bpRNA-1m database will be valuable resources both for specific analysis of individual RNA molecules and large-scale analyses such as are useful for updating RNA energy parameters for computational thermodynamic predictions, improving machine learning models for structure prediction, and for benchmarking structure-prediction algorithms.

Sensitive and specific detection of microRNAs based on two-stage amplification reaction using molecular beacons as turn-on probes

In this study, a rapid, sensitive, and specific assay for detecting miRNAs was developed based on a two-stage amplification reaction (TSAR) using molecular beacons (MBs) as turn-on probes. In the TSAR, different miRNAs can be converted to the same reporter oligonucleotides (Y), which can hybridize with the same MB. Therefore, in combination with specific templates, this method can be applied to multiplex miRNA detection by simply using the same MB. The loop region of the MB was screened by computer simulation methods. In particular, to improve the specificity of the MB in real sample analysis, the maximum similarity of the MB loop region to the human genome and human transcriptome is less than 70%. Two MBs were designed in this study. MB I, with nine flanking base pairs in its stem region, was used for real-time monitoring of the production of Y during the TSAR. MB II, with five flanking base pairs in its stem region, was used to detect the production of Y at the end of the TSAR. This assay exhibited high sensitivity with a limit of detection of 2.0 × 10−16M and 6.7 × 10−16M using MB I and MB II as turn-on probes, respectively. In addition, this assay can clearly discriminate single base differences in miRNA sequences, and the TSAR can be completed under isothermal conditions. Accordingly, the isothermal reaction conditions and simple fluorescence measurement can greatly contribute to the development of a fast point-of-care detection system.

Structure and Dynamics of the Tetra-A Loop and (A-A)-U Sequence Motif within the Coliphage GA Replicase RNA Operator.

The three-dimensional structure of a RNA hairpin containing the RNA operator binding site for bacteriophage GA coat protein is presented. The phage GA operator contains the asymmetric (A-A)-U sequence motif and is capped by a four-adenine (tetra-A) loop. The uridine of the (A-A)-U motif preferentially pairs with the 5'-proximal cross-strand adenine, and the 3'-proximal adenine stacks into the helix. The tetra-A loop is well-ordered with adenine residues 2-4 forming a 3' stack. This loop conformation stands in contrast to the structure of the 5'-AUUA loop of the related phage MS2 operator in which residues 1 and 2 form a 5' stack. The context dependence of the (A-A)-U sequence motif conformation was examined using structures of 76 unique occurrences from the Protein Data Bank. The motif almost always has one adenine bulged and the other adenine adopting an A-U base pair. In the case in which the (A-A)-U motif is flanked by only one Watson-Crick base pair, the adenine adjacent to the flanking base pair tends to bulge; 80% of motifs with a 3' flanking pair have a 3' bulged adenine, and 84% of motifs with a 5' flanking pair have a 5' bulged adenine. The frequencies of 3'- and 5'-proximal adenines bulging are 33 and 67%, respectively, when the (A-A)-U motif is flanked by base pairs on both sides. Although a 3' flanking cytidine correlates (88%) with bulging of the 5'-proximal adenine, no strict dependence on flanking nucleotide identity was identified for the 5' side.

DNA Sequence Determinants Controlling Affinity, Stability and Shape of DNA Complexes Bound by the Nucleoid Protein Fis.

The abundant Fis nucleoid protein selectively binds poorly related DNA sequences with high affinities to regulate diverse DNA reactions. Fis binds DNA primarily through DNA backbone contacts and selects target sites by reading conformational properties of DNA sequences, most prominently intrinsic minor groove widths. High-affinity binding requires Fis-stabilized DNA conformational changes that vary depending on DNA sequence. In order to better understand the molecular basis for high affinity site recognition, we analyzed the effects of DNA sequence within and flanking the core Fis binding site on binding affinity and DNA structure. X-ray crystal structures of Fis-DNA complexes containing variable sequences in the noncontacted center of the binding site or variations within the major groove interfaces show that the DNA can adapt to the Fis dimer surface asymmetrically. We show that the presence and position of pyrimidine-purine base steps within the major groove interfaces affect both local DNA bending and minor groove compression to modulate affinities and lifetimes of Fis-DNA complexes. Sequences flanking the core binding site also modulate complex affinities, lifetimes, and the degree of local and global Fis-induced DNA bending. In particular, a G immediately upstream of the 15 bp core sequence inhibits binding and bending, and A-tracts within the flanking base pairs increase both complex lifetimes and global DNA curvatures. Taken together, our observations support a revised DNA motif specifying high-affinity Fis binding and highlight the range of conformations that Fis-bound DNA can adopt. The affinities and DNA conformations of individual Fis-DNA complexes are likely to be tailored to their context-specific biological functions.

DNA Distortion Caused by Uracil-Containing Intrastrand Cross-Links.

Four uracil-containing intrastrand cross-links have been detected in human cells upon UV irradiation of 5-bromouracil-containing DNA, namely 5'-G[8-5]U-3', 5'-U[5-8]G-3', 5'-A[8-5]U-3', and 5'-A[2-5]U-3'. These lesions feature unique composition and connectivity compared with other intrastrand cross-links reported in the literature. For the first time, structural information obtained using molecular dynamics (MD) simulations reveal that all four lesions distort the DNA helix, which can involve an extrahelical location of the cross-link, changes in the helical interactions of the complementary nucleotides, or disruption of hydrogen bonding in the flanking base pairs up to two positions from the cross-linked site; however, the degree of distortion varies between the cross-links, being affected by the sequence, nucleobase-nucleobase connectivity, and the purine involved. Most importantly, the relative distortion of the damaged DNA provides the first structural explanation for the observed abundances of the four uracil-containing cross-links. Furthermore, the highly distorted conformations suggest that these lesions will likely have severe implications for DNA replication and repair processes in cells.

Conformational Insights into the Mechanism of Acetylaminofluorene-dG-Induced Frameshift Mutations in the NarI Mutational Hotspot.

Frameshift mutagenesis encompasses the gain or loss of DNA base pairs, resulting in altered genetic outcomes. The NarI restriction site sequence 5'-G1G2CG3CX-3' in Escherichia coli is a well-known mutational hotspot, in which lesioning of acetylaminofluorene (AAF) at G3* induces a greater -2 deletion frequency than that at other guanine sites. Its mutational efficiency is modulated by the nature of the nucleotide in the X position (C ∼ A > G ≫ T). Here, we conducted a series of polymerase-free solution experiments that examine the conformational and thermodynamic basis underlying the propensity of adducted G3 to form a slipped mutagenic intermediate (SMI) and its sequence dependence during translesion synthesis (TLS). Instability of the AAF-dG3:dC pair at the replication fork promoted slippage to form a G*C bulge-out SMI structure, consisting of S- ("lesion stacked") and B-SMI ("lesion exposed") conformations, with conformational rigidity increasing as a function of primer elongation. We found greater stability of the S- compared to the B-SMI conformer throughout TLS. The dependence of their population ratios was determined by the 3'-next flanking base X at fully elongated bulge structures, with 59% B/41% S and 86% B/14% S for the dC and dT series, respectively. These results indicate the importance of direct interactions of the hydrophobic AAF lesion with the 3'-next flanking base pair and its stacking fit within the -2 bulge structure. A detailed conformational understanding of the SMI structures and their sequence dependence may provide a useful model for DNA polymerase complexes.

Probing the run-on oligomer of activated SgrAI bound to DNA.

SgrAI is a type II restriction endonuclease with an unusual mechanism of activation involving run-on oligomerization. The run-on oligomer is formed from complexes of SgrAI bound to DNA containing its 8 bp primary recognition sequence (uncleaved or cleaved), and also binds (and thereby activates for DNA cleavage) complexes of SgrAI bound to secondary site DNA sequences which contain a single base substitution in either the 1st/8th or the 2nd/7th position of the primary recognition sequence. This modulation of enzyme activity via run-on oligomerization is a newly appreciated phenomenon that has been shown for a small but increasing number of enzymes. One outstanding question regarding the mechanistic model for SgrAI is whether or not the activating primary site DNA must be cleaved by SgrAI prior to inducing activation. Herein we show that an uncleavable primary site DNA containing a 3’-S-phosphorothiolate is in fact able to induce activation. In addition, we now show that cleavage of secondary site DNA can be activated to nearly the same degree as primary, provided a sufficient number of flanking base pairs are present. We also show differences in activation and cleavage of the two types of secondary site, and that effects of selected single site substitutions in SgrAI, as well as measured collisional cross-sections from previous work, are consistent with the cryo-electron microscopy model for the run-on activated oligomer of SgrAI bound to DNA.

The effect of sequence context on the activity of cytosine DNA glycosylases.

We have prepared single (N204D) and double (N204D:L272A) mutants of human uracil DNA glycosylase (hUDG), generating two cytosine DNA glycosylases (hCDG and hCYDG). Both these enzymes are able to excise cytosine (but not 5-methylcytosine), when this base is part of a mismatched base pair. hCDG is more active than the equivalent E. coli enzyme (eCYDG) and also has some activity when the cytosine is paired with guanine, unlike eCYDG. hCDG also has some activity against single stranded DNA, while having poor activity towards an unnatural base pair that forces the cytosine into an extrahelical conformation (in contrast to eCYDG for which a bulky base enhances the enzyme's activity). We also examined how sequence context affects the activity of these enzymes, determining the effect of flanking base pairs on cleavage efficiency. An abasic site or a hexaethylene glycol linker placed opposite the target cytosine, also causes an increase in activity compared with an AC mismatch. Flanking an AC mismatch with GC base pairs resulted in a 100-fold decrease in excision activity relative to flanking AT base pairs and the 5'-flanking base pair had a greater effect on the rate of cleavage. However, this effect is not simply due to the stability of the flanking base pairs as adjacent GT mismatches also produce low cleavage efficiency.

An isocytidine derivative with a 2-amino-6-methylpyridine unit for selective recognition of the CG interrupting site in an antiparallel triplex DNA.

Sequence-specific recognition of duplex DNA mediated by triple helix formation offers a potential basis for oligonucleotide therapy and biotechnology. However, triplex formation is limited mostly to homopurine strands, due to poor stabilization at CG or TA base pairs in the target duplex DNA sequences. Several non-natural nucleosides have been designed for the recognition of CG or TA base pairs within an antiparallel triplex DNA. Nevertheless, problems including low selectivity and high dependence on the neighboring bases remain unsolved. We thus synthesized N(2)-arylmethyl isodC derivatives and incorporated them into triplex-forming oligonucleotides (TFOs) for the selective recognition of the CG base pair within antiparallel triplex DNA. It was shown that an isodC derivative bearing a 2-amino-6-methylpyridine moiety (AP-isodC) recognizes the CG base pair with high selectivity in antiparallel triplex DNA irrespective of the flanking base pairs.

PKa Shifting in Double-Stranded RNA Is Highly Dependent upon Nearest Neighbors and Bulge Positioning

Shifting of pKa's in RNA is important for many biological processes; however, the driving forces responsible for shifting are not well understood. Herein, we determine how structural environments surrounding protonated bases affect pKa shifting in double-stranded RNA (dsRNA). Using (31)P NMR, we determined the pKa of the adenine in an A(+)·C base pair in various sequence and structural environments. We found a significant dependence of pKa on the base pairing strength of nearest neighbors and the location of a nearby bulge. Increasing nearest neighbor base pairing strength shifted the pKa of the adenine in an A(+)·C base pair higher by an additional 1.6 pKa units, from 6.5 to 8.1, which is well above neutrality. The addition of a bulge two base pairs away from a protonated A(+)·C base pair shifted the pKa by only ~0.5 units less than a perfectly base paired hairpin; however, positioning the bulge just one base pair away from the A(+)·C base pair prohibited formation of the protonated base pair as well as several flanking base pairs. Comparison of data collected at 25 °C and 100 mM KCl to biological temperature and Mg(2+) concentration revealed only slight pKa changes, suggesting that similar sequence contexts in biological systems have the potential to be protonated at biological pH. We present a general model to aid in the determination of the roles protonated bases may play in various dsRNA-mediated processes including ADAR editing, miRNA processing, programmed ribosomal frameshifting, and general acid-base catalysis in ribozymes.

Thermodynamic Profiles and Nuclear Magnetic Resonance Studies of Oligonucleotide Duplexes Containing Single Diastereomeric Spiroiminodihydantoin Lesions

The spiroiminodihydantoins (Sp) are highly mutagenic oxidation products of guanine and 8-oxo-7,8-dihydroguanine in DNA. The Sp lesions have recently been detected in the liver and colon of mice infected with Helicobacter hepaticus that induces inflammation and the development of liver and colon cancers in murine model systems [Mangerich, A., et al. (2012) Proc. Natl. Acad. Sci. U.S.A. 109, E1820-E1829]. The impact of Sp lesions on the thermodynamic characteristics and the effects of the diastereomeric Sp-R and Sp-S lesions on the conformational features of double-stranded 11-mer oligonucleotide duplexes have been studied by a combination of microcalorimetric methods, analysis of DNA melting curves, and two-dimensional nuclear magnetic resonance methods. The nonplanar, propeller-like shapes of the Sp residues strongly diminish the extent of local base stacking interactions that destabilize the DNA duplexes characterized by unfavorable enthalpy contributions. Relative to that of an unmodified duplex, the thermally induced unfolding of the duplexes with centrally positioned Sp-R and Sp-S lesions into single strands is accompanied by a smaller release of cationic counterions (Δn(Na⁺) = 0.6 mol of Na⁺/mol of duplex) and water molecules (Δn(w) = 17 mol of H₂O/mol of duplex). The unfolding parameters are similar for the Sp-R and Sp-S lesions, although their orientations in the duplexes are different. The structural disturbances radiate one base pair beyond the flanking C:G pair, although Watson-Crick hydrogen bonding is maintained at all flanking base pairs. The observed relatively strong destabilization of B-form DNA by the physically small Sp lesions is expected to have a significant impact on the processing of these lesions in biological environments.