Nearest-neighbor thermodynamics of deoxyinosine pairs in DNA duplexes

Nearest-neighbor thermodynamic parameters of the ‘universal pairing base’ deoxyinosine were determined for the pairs I·C, I·A, I·T, I·G and I·I adjacent to G·C and A·T pairs. Ultraviolet absorbance melting curves were measured and non-linear regression performed on 84 oligonucleotide duplexes with 9 or 12 bp lengths. These data were combined with data for 13 inosine containing duplexes from the literature. Multiple linear regression was used to solve for the 32 nearest-neighbor unknowns. The parameters predict the T_m for all sequences within 1.2°C on average. The general trend in decreasing stability is I·C > I·A > I·T ≈ I· G > I·I. The stability trend for the base pair 5′ of the I·X pair is G·C > C·G > A·T > T·A. The stability trend for the base pair 3′ of I·X is the same. These trends indicate a complex interplay between H-bonding, nearest-neighbor stacking, and mismatch geometry. A survey of 14 tandem inosine pairs and 8 tandem self-complementary inosine pairs is also provided. These results may be used in the design of degenerate PCR primers and for degenerate microarray probes.     

Influence of flanking sequence context on the conformational flexibility of aminofluorene-modified dG adduct in dA mismatch DNA duplexes

When positioned opposite to a dA in a DNA duplex, the prototype arylamine–DNA adduct [N-(2′-deoxyguanosin-yl)-7-fluoro-2-aminofluorene (FAF)] adopts the so-called ‘wedge’ (W) conformation, in which the carcinogen resides in the minor groove of the duplex. All 16 FAF-modified 12-mer NG*N/NAN dA mismatch duplexes (G* = FAF, N = G, A, C, T) exhibited strongly positive induced circular dichroism in the 290–360 nm range (ICD_{290–360 nm}), which supports the W conformation. The ICD_{290–360 nm} intensities were the greatest for duplexes with a 3′-flanking T. The AG*N duplex series showed little adduct-induced destabilization. An exception was the AG*T duplex, which displayed two well-resolved signals in the ¹⁹F NMR spectra. This was presumably due to a strong lesion-destabilizing effect of the 3′-T. The flanking T effect was substantiated further by findings with the TG*T duplex, which exhibited greater lesion flexibility and nucleotide excision repair recognition. Adduct conformational heterogeneity decreased in order of TG*T > AG*T > CG*T > AG*A > AG*G > AG*C. The dramatic flanking T effect on W-conformeric duplexes is consistent with the strong dependence of the ICD_290-360 on both temperature and salt concentration and could be extended to the arylamine food mutagens that are biologically relevant in humans.     

A unique pattern of intrastrand anomalies in base composition of the DNA in hypotrichs

The 50 non-coding bases immediately internal to the telomeric repeats in the two 5′ ends of macronuclear DNA molecules of a group of hypotrichous ciliates are anomalous in composition, consisting of 61% purines and 39% pyrimidines, A>T (ratio of 44:32), and G>C (ratio of 17:7). These ratio imbalances violate parity rule 2, according to which A should equal T and G should equal C within a DNA strand and therefore pyrimidines should equal purines. The purine-rich and base ratio imbalances are in marked contrast to the rest of the non-coding parts of the molecules, which have the theoretically expected purine content of 50%, with A = T and G = C. The ORFs contain an average of 52% purines as a result of bias in codon usage. The 50 bases that flank the 5′ ends of macronuclear sequences in micronuclear DNA (12 cases) consist of ~50% purines. Thus, the 50 bases in the 5′ ends of macronuclear sequences in micronuclear DNA are islands of purine richness in which A>T and G>C. These islands may serve as signals for the excision of macronuclear molecules during macronuclear development. We have found no published reports of coding or non-coding native DNA with such anomalous base composition.     

Evaluation of complexation of metal-mediated DNA-binding drugs to oligonucleotides via electrospray ionization mass spectrometry

The interactions of self-complementary oligonucleotides with a group of metal-mediated DNA-binding drugs, including chromomycin A₃, mithramycin and the novel compound UK-1, were examined via electrospray ionization quadrupole ion trap mass spectrometry. Both chromomycin and mithramycin were shown to bind preferentially to GC-rich oligonucleotide duplexes in a 2:1 drug:metal ratio, while UK-1 was shown to bind in a 1:1 drug:metal stoichiometric ratio without a strong sequence preference. These trends were observed in the presence of Co²⁺, Ni²⁺ and Zn²⁺, with the exception that chromomycin–Zn²⁺ complexes were not readily observed. The binding stoichiometries as well as the sequence specificities are in agreement with literature reports for solution studies. Binding selectivities and stabilities of the complexes were also probed using electrospray ionization mass spectrometry. Both of the GC-rich oligomers 5′-GCGCGC-3′ and 5′-GCGCATGCGC-3′ exhibited a binding preference for chromomycin over mithramycin in the presence of Co²⁺ and Ni²⁺. Energy-variable collisionally activated dissociation of the complexes was employed to determine the stabilities of the complexes. The relative metal-dependent binding energies were Ni²⁺ > Zn²⁺ > Co²⁺ for UK-1–oligomer complexes and Ni²⁺ > Co²⁺ for both mithramycin and chromomycin complexes.     

hnRNP H binding at the 5' splice site correlates with the pathological effect of two intronic mutations in the NF-1 and TSHbeta genes

We have recently reported a disease-causing substitution (+5G > C) at the donor site of NF-1 exon 3 that produces its skipping. We have now studied in detail the splicing mechanism involved in analyzing RNA–protein complexes at several 5′ splice sites. Characteristic protein patterns were observed by pulldown and band-shift/super-shift analysis. Here, we show that hnRNP H binds specifically to the wild-type GGGgu donor sequence of the NF-1 exon 3. Depletion analyses shows that this protein restricts the accessibility of U1 small nuclear ribonucleoprotein (U1snRNA) to the donor site. In this context, the +5G > C mutation abolishes both U1snRNP base pairing and the 5′ splice site (5′ss) function. However, exon recognition in the mutant can be rescued by disrupting the binding of hnRNP H, demonstrating that this protein enhances the effects of the +5G > C substitution. Significantly, a similar situation was found for a second disease-causing +5G > A substitution in the 5′ss of TSHβ exon 2, which harbors a GGgu donor sequence. Thus, the reason why similar nucleotide substitutions can be either neutral or very disruptive of splicing function can be explained by the presence of specific binding signatures depending on local contexts.     

Ligation reaction specificities of an NAD(+)-dependent DNA ligase from the hyperthermophile Aquifex aeolicus

An NAD⁺-dependent DNA ligase from the hyperthermophilic bacterium Aquifex aeolicus was cloned, expressed in Escherichia coli and purified to homogeneity. The enzyme is most active in slightly alkaline pH conditions with either Mg²⁺ or Mn²⁺ as the metal cofactor. Ca²⁺ and Ni²⁺ mainly support formation of DNA–adenylate intermediates. The catalytic cycle is characterized by a low k_cat value of 2 min^–1 with concomitant accumulation of the DNA–adenylate intermediate when Mg²⁺ is used as the metal cofactor. The ligation rates of matched substrates vary by up to 4-fold, but exhibit a general trend of T/A ≤ G/C < C/G < A/T on both the 3′- and 5′-side of the nick. Consistent with previous studies on Thermus ligases, this Aquifex ligase exhibits greater discrimination against a mismatched base pair on the 3′-side of the nick junction. The requirement of 3′ complementarity for a ligation reaction is reaffirmed by results from 1 nt insertions on either the 3′- or 5′-side of the nick. Furthermore, most of the unligatable 3′ mismatched base pairs prohibit formation of the DNA–adenylate intermediate, indicating that the substrate adenylation step is also a control point for ligation fidelity. Unlike previously studied ATP ligases, gapped substrates cannot be ligated and intermediate accumulation is minimal, suggesting that complete elimination of base pair complementarity on one side of the nick affects substrate adenylation on the 5′-side of the nick junction. Relationships among metal cofactors, ligation products and intermediate, and ligation fidelity are discussed.     

Pronounced instability of tandem IU base pairs in RNA

Optical melting was used to determine the stabilities of three series of RNA oligomers containing tandem XU base pairs, GGCXUGCC (5′XU3′), GGCUXGCC (5′UX3′) and GGCXXGGC/CCGUUCCG (5′XX3′), where X is either A, G or I (inosine). The helices containing tandem AU base pairs were the most stable in the first two series (5′XU3′ and 5′UX3′), with an average melting temperature ~11°C higher than the helices with tandem 5′GU3′ base pairs and 25°C higher than the helices with tandem 5′IU3′ base pairs. For the third series (5′XX3′), the helix containing tandem GG is the most stable, with an average melting temperature ~2°C higher than the helix with tandem AA base pairs and ~24°C higher than the helix with tandem II base pairs. The thermodynamic stability of the oligomers with tandem IU base pairs was also investigated as a function of magnesium ion concentration. As with normal A–U or G–U tandem duplexes, the data could best be interpreted as non-specific binding of magnesium ions to the inosine-containing RNA oligonucleotides.     

Gel-based oligonucleotide microarray approach to analyze protein–ssDNA binding specificity

Gel-based oligonucleotide microarray approach was developed for quantitative profiling of binding affinity of a protein to single-stranded DNA (ssDNA). To demonstrate additional capabilities of this method, we analyzed the binding specificity of ribonuclease (RNase) binase from Bacillus intermedius (EC 3.1.27.3) to ssDNA using generic hexamer oligodeoxyribonucleotide microchip. Single-stranded octamer oligonucleotides were immobilized within 3D hemispherical gel pads. The octanucleotides in individual pads 5′-{N}N₁N₂N₃N₄N₅N₆{N}-3′ consisted of a fixed hexamer motif N₁N₂N₃N₄N₅N₆ in the middle and variable parts {N} at the ends, where {N} represent A, C, G and T in equal proportions. The chip has 4096 pads with a complete set of hexamer sequences. The affinity was determined by measuring dissociation of the RNase–ssDNA complexes with the temperature increasing from 0°C to 50°C in quasi-equilibrium conditions. RNase binase showed the highest sequence-specificity of binding to motifs 5′-NNG(A/T/C)GNN-3′ with the order of preference: GAG > GTG > GCG. High specificity towards G(A/T/C)G triplets was also confirmed by measuring fluorescent anisotropy of complexes of binase with selected oligodeoxyribonucleotides in solution. The affinity of RNase binase to other 3-nt sequences was also ranked. These results demonstrate the applicability of the method and provide the ground for further investigations of nonenzymatic functions of RNases.     

Impact of the six nucleotides downstream of the stop codon on translation termination

The efficiency of translation termination is influenced by local contexts surrounding stop codons. In Saccharomyces cerevisiae, upstream and downstream sequences act synergistically to influence the translation termination efficiency. By analysing derivatives of a leaky stop codon context, we initially demonstrated that at least six nucleotides after the stop codon are a key determinant of readthrough efficiency in S. cerevisiae. We then developed a combinatorial-based strategy to identify poor 3′ termination contexts. By screening a degenerate oligonucleotide library, we identified a consensus sequence –CA(A/G)N(U/C/G)A–, which promotes >5% readthrough efficiency when located downstream of a UAG stop codon. Potential base pairing between this stimulatory motif and regions close to helix 18 and 44 of the 18S rRNA provides a model for the effect of the 3′ stop codon context on translation termination.     

Synthesis and physical and physiological properties of 4'-thioRNA: application to post-modification of RNA aptamer toward NF-kappaB

We report herein full details of the preparation of 4′-thiouridine, -cytidine, -adenosine and -guanosine phosphoramidites based on our synthetic protocol via the Pummerer reaction. Fully modified 4′-thioRNAs containing four kinds of 4′-thioribonucleoside units were prepared according to the standard RNA synthesis. The T_m values and thermodynamic parameters of a series of duplexes were determined by UV melting and differential scanning calorimetry (DSC) measurements. The resulting overall order of thermal stabilities for the duplexes was 4′-thioRNA:4′-thioRNA >> 4′-thioRNA:RNA > RNA:RNA > RNA:DNA > 4′-thioRNA:DNA. In addition, it was shown that the dominant factor in the stability of the duplexes consisting of 4′-thioRNA was enthalpic in character. The CD spectra of not only 4′-thioRNA:RNA and 4′-thioRNA:4′-thioRNA but also 4′-thioRNA:DNA were all similar to those of duplexes in the A-conformation. The stability of 4′-thioRNA in human serum was 600 times greater than that of natural RNA. Neither the RNA:RNA nor the 4′-thioRNA:4′-thioRNA duplexes were digested under the same conditions. The first example of a post-modification of an RNA aptamer by 4′-thioribonucleoside units was demonstrated. Full modification of the aptamer thioRNA3 resulted in complete loss of binding activity. In contrast, modifications at positions other than the binding site were tolerated without loss of binding activity. The post-modified RNA aptamer thioRNA5 was thermally stabilized and resistant toward nuclease digestion. The results presented in this paper will, it is hoped, contribute to the development of 4′-thioRNA as a new generation of artificial RNA.     

Theoretical Investigation on the Binding Specificity of Sialyldisaccharides with Hemagglutinins of Influenza A Virus by Molecular Dynamics Simulations

Recognition of cell-surface sialyldisaccharides by influenza A hemagglutinin (HA) triggers the infection process of influenza. The changes in glycosidic torsional linkage and the receptor conformations may alter the binding specificity of HAs to the sialylglycans. In this study, 10-ns molecular dynamics simulations were carried out to examine the structural and dynamic behavior of the HAs bound with sialyldisaccharides Neu5Acα(2–3)Gal (N23G) and Neu5Acα(2–6)Gal (N26G). The analysis of the glycosidic torsional angles and the pair interaction energy between the receptor and the interacting residues of the binding site reveal that N23G has two binding modes for H1 and H5 and a single binding mode for H3 and H9. For N26G, H1 and H3 has two binding modes, and H5 and H9 has a single binding mode. The direct and water-mediated hydrogen bonding interactions between the receptors and HAs play dominant roles in the structural stabilization of the complexes. It is concluded from pair interaction energy and Molecular Mechanic-Poisson-Boltzmann Surface Area calculations that N26G is a better receptor for H1 when compared with N23G. N23G is a better receptor for H5 when compared with N26G. However, H3 and H9 can recognize N23G and N26G in equal binding specificity due to the marginal energy difference (≈2.5 kcal/mol). The order of binding specificity of N23G is H3 > H5 > H9 > H1 and N26G is H1 > H3 > H5 > H9, respectively. The proposed conformational models will be helpful in designing inhibitors for influenza virus.     

GAD2 on chromosome 10p12 is a candidate gene for human obesity

The gene GAD2 encoding the glutamic acid decarboxylase enzyme (GAD65) is a positional candidate gene for obesity on Chromosome 10p11–12, a susceptibility locus for morbid obesity in four independent ethnic populations. GAD65 catalyzes the formation of γ-aminobutyric acid (GABA), which interacts with neuropeptide Y in the paraventricular nucleus to contribute to stimulate food intake. A case-control study (575 morbidly obese and 646 control subjects) analyzing GAD2 variants identified both a protective haplotype, including the most frequent alleles of single nucleotide polymorphisms (SNPs) +61450 C>A and +83897 T>A (OR = 0.81, 95% CI [0.681–0.972], p = 0.0049) and an at-risk SNP (−243 A>G) for morbid obesity (OR = 1.3, 95% CI [1.053–1.585], p = 0.014). Furthermore, familial-based analyses confirmed the association with the obesity of SNP +61450 C>A and +83897 T>A haplotype (χ² = 7.637, p = 0.02). In the murine insulinoma cell line βTC3, the G at-risk allele of SNP −243 A>G increased six times GAD2 promoter activity (p < 0.0001) and induced a 6-fold higher affinity for nuclear extracts. The −243 A>G SNP was associated with higher hunger scores (p = 0.007) and disinhibition scores (p = 0.028), as assessed by the Stunkard Three-Factor Eating Questionnaire. As GAD2 is highly expressed in pancreatic β cells, we analyzed GAD65 antibody level as a marker of β-cell activity and of insulin secretion. In the control group, −243 A>G, +61450 C>A, and +83897 T>A SNPs were associated with lower GAD65 autoantibody levels (p values of 0.003, 0.047, and 0.006, respectively). SNP +83897 T>A was associated with lower fasting insulin and insulin secretion, as assessed by the HOMA-B% homeostasis model of β-cell function (p = 0.009 and 0.01, respectively). These data support the hypothesis of the orexigenic effect of GABA in humans and of a contribution of genes involved in GABA metabolism in the modulation of food intake and in the development of morbid obesity.     

Thermodynamics of RNA duplexes modified with unlocked nucleic acid nucleotides

Thermodynamics provides insights into the influence of modified nucleotide residues on stability of nucleic acids and is crucial for designing duplexes with given properties. In this article, we introduce detailed thermodynamic analysis of RNA duplexes modified with unlocked nucleic acid (UNA) nucleotide residues. We investigate UNA single substitutions as well as model mismatch and dangling end effects. UNA residues placed in a central position makes RNA duplex structure less favourable by 4.0–6.6 kcal/mol. Slight destabilization, by ∼0.5–1.5 kcal/mol, is observed for 5′- or 3′-terminal UNA residues. Furthermore, thermodynamic effects caused by UNA residues are extremely additive with ΔG°₃₇ conformity up to 98%. Direct mismatches involving UNA residues decrease the thermodynamic stability less than unmodified mismatches in RNA duplexes. Additionally, the presence of UNA residues adjacent to unpaired RNA residues reduces mismatch discrimination. Thermodynamic analysis of UNA 5′- and 3′-dangling ends revealed that stacking interactions of UNA residues are always less favourable than that of RNA residues. Finally, circular dichroism spectra imply no changes in overall A-form structure of UNA–RNA/RNA duplexes relative to the unmodified RNA duplexes.     

Misincorporation of dAMP opposite 2-hydroxyadenine, an oxidative form of adenine.

Nucleotide incorporation opposite an oxidative form of adenine, 2-hydroxyadenine (2-OH-Ade) was investigated. When a primed template with 2-OH-Ade was treated with an exonuclease-deficient Klenow fragment of Escherichia coli DNA polymerase I (KFexo-), recombinant rat DNA polymerase beta (pol beta) or calf thymus DNA polymerase alpha (pol alpha), incorporation of dTMP and dAMP was observed. In addition, KFexo- inserted dGMP as well. A steady-state kinetic study indicated that the insertion of dAMP and dTMP opposite the DNA lesion occurred with similar frequency with KFexo- and pol beta. Insertion of dTMP opposite 2-OH-Ade was favored to that of dAMP by pol alpha. Chain extension from the A.2-OH-Ade pair is less favored than that from the T.2-OH-Ade pair by all three DNA polymerase. Analysis of full-length products of in vitro DNA synthesis showed that dTMP and dAMP were incorporated by DNA polymerases and that exonuclease-proficient and -deficient Klenow fragments also inserted dGMP opposite 2-OH-Ade. These results suggest that formation of 2-OH-Ade from A in DNA will induce A-->T and A-->C transversions in cells.     

Conservation and Diversity of Influenza A H1N1 HLA-Restricted T Cell Epitope Candidates for Epitope-Based Vaccines

Background

The immune-related evolution of influenza viruses is exceedingly complex and current vaccines against influenza must be reformulated for each influenza season because of the high degree of antigenic drift among circulating influenza strains. Delay in vaccine production is a serious problem in responding to a pandemic situation, such as that of the current H1N1 strain. Immune escape is generally attributed to reduced antibody recognition of the viral hemagglutinin and neuraminidase proteins whose rate of mutation is much greater than that of the internal non-structural proteins. As a possible alternative, vaccines directed at T cell epitope domains of internal influenza proteins, that are less susceptible to antigenic variation, have been investigated.

Methodology/Principal Findings

HLA transgenic mouse strains expressing HLA class I A*0201, A*2402, and B*0702, and class II DRB1*1501, DRB1*0301 and DRB1*0401 were immunized with 196 influenza H1N1 peptides that contained residues of highly conserved proteome sequences of the human H1N1, H3N2, H1N2, H5N1, and avian influenza A strains. Fifty-four (54) peptides that elicited 63 HLA-restricted peptide-specific T cell epitope responses were identified by IFN-γ ELISpot assay. The 54 peptides were compared to the 2007–2009 human H1N1 sequences for selection of sequences in the design of a new candidate H1N1 vaccine, specifically targeted to highly-conserved HLA-restricted T cell epitopes.

Conclusions/Significance

Seventeen (17) T cell epitopes in PB1, PB2, and M1 were selected as vaccine targets based on sequence conservation over the past 30 years, high functional avidity, non-identity to human peptides, clustered localization, and promiscuity to multiple HLA alleles. These candidate vaccine antigen sequences may be applicable to any avian or human influenza A virus.     

The crystal structure of an alternating RNA heptamer r(GUAUACA) forming a six base-paired duplex with 3'-end adenine overhangs

The crystal structure of an alternating RNA heptamer r(GUAUACA) has been determined to 2.0 Å resolution and refined to an R_work of 17.1% and R_free of 18.5% using 2797 reflections. The heptamer crystallized in the space group C222 with a unit cell of a = 25.74, b = 106.58, c = 30.26 Å and two independent strands in the asymmetric unit. Each heptamer forms a duplex with its symmetry-related strand and each duplex contains six Watson–Crick base pairs and 3′-end adenosine overhangs. Therefore, two kinds of duplex (duplex 1 and duplex 2) are formed. Duplexes 1 stack on each other forming a pseudo-continuous column, which is typical of the RNA packing mode, while duplex 2 is typical of A-DNA packing with its termini in abutting interactions. Overhang adenine residues stack within the duplexes with C3′-endo sugar pucker and C2′-endo sugar pucker in duplexes 1 and 2, respectively. A Na⁺ ion in the crystal lattice is water bridged to two N1 atoms of symmetry-related A7 bases.     

Delineation of the mechanisms of aberrant splicing caused by two unusual intronic mutations in the RSK2 gene involved in Coffin-Lowry syndrome

Coffin–Lowry syndrome (CLS) is caused by mutations in the RSK2 gene encoding a protein kinase of the Ras signalling pathway. We have studied two point mutations which cause aberrant splicing but do not concern the invariant GT or AG nucleotides of splice sites. The first, an A→G transition at position +3 of the 5′ splice site of exon 6, results in vivo and in vitro in exon skipping and premature translation termination. The natural 5′ splice site, although intrinsically weak, is not transactivated under normal conditions. Consequently, replacement of an A/U by a G/U base pairing with U1 snRNA reduces its strength below a critical threshold. The second mutation, an A→G transition 11 nt upstream of exon 5, creates a new AG near the natural 3′ splice site. In vitro this synthetic 3′ AG is used exclusively by the splicing machinery. In vivo this splicing event is also observed, but is underestimated because the resulting RSK2 mRNA contains premature stop codons which trigger the nonsense-mediated decay process. We show that a particular mechanism is involved in the aberrant splicing of exon 5, implying involvement of the natural 3′ AG during the first catalytic step and the new 3′ AG during the second step. Thus, our results explain how these mutations cause severe forms of CLS.     

New insights into Hoogsteen base pairs in DNA duplexes from a structure-based survey

Hoogsteen (HG) base pairs (bps) provide an alternative pairing geometry to Watson–Crick (WC) bps and can play unique functional roles in duplex DNA. Here, we use structural features unique to HG bps (syn purine base, HG hydrogen bonds and constricted C1′–C1′ distance across the bp) to search for HG bps in X-ray structures of DNA duplexes in the Protein Data Bank. The survey identifies 106 A•T and 34 G•C HG bps in DNA duplexes, many of which are undocumented in the literature. It also uncovers HG-like bps with syn purines lacking HG hydrogen bonds or constricted C1′–C1′ distances that are analogous to conformations that have been proposed to populate the WC-to-HG transition pathway. The survey reveals HG preferences similar to those observed for transient HG bps in solution by nuclear magnetic resonance, including stronger preferences for A•T versus G•C bps, TA versus GG steps, and also suggests enrichment at terminal ends with a preference for 5′-purine. HG bps induce small local perturbations in neighboring bps and, surprisingly, a small but significant degree of DNA bending (∼14°) directed toward the major groove. The survey provides insights into the preferences and structural consequences of HG bps in duplex DNA.