Interaction of λ Cro repressor and its V55C mutant S-S dimer with symmetrical and asymmetrical DNA

Interaction of phage λ repressor Cro and its mutant form CroV55C, an S-S-bonded dimer that arises spontaneously as a result of substitution of Val55 with Cys, with synthetic oligonucleotide duplexes was studied using competition with the antibiotic distamycin A. The equilibrium binding constants of Cro with phage λ operator O_R3 and various oduplexes containing right and left O_R3 half-operators with single-basepair deletions and insertions in the center of duplexes were calculated. The highest binding constant was obtained when Cro was adsorbed on the symmetrical 17-bp complex O_S, composed of two left half-operators. Deletion of the central GC pair or insertion of an additional GC in O_S decreased the binding constant by a factor of 4–5. The Cro binding constant for the symmetrical 17-bp duplex O_D consisting of two right half-operators and its analog with an additional central GC was five times lower than that for O_S. The thermodynamic binding characteristics of some oligonucleotide duplexes with CroV55C were determined. The binding constant for CroV55C with O_R3 was lower by an order of magnitude than that of Cro, whereas with O_D the CroV55C and Cro binding constants were much the same. For CroV55C and a symmetrical analog of O_S with the central GC deleted, the binding constant increased by an order of magnitude over that for CroV55C and O_R3.     

Thermodynamics of RNA duplexes with tandem mismatches containing a uracil-uracil pair flanked by C.G/G.C or G.C/A.U closing base pairs

The thermodynamics governing the denaturation of RNA duplexes containing 8 bp and a central tandem mismatch or 10 bp were evaluated using UV absorbance melting curves. Each of the eight tandem mismatches that were examined had one U-U pair adjacent to another noncanonical base pair. They were examined in two different RNA duplex environments, one with the tandem mismatch closed by G.C base pairs and the other with G.C and A.U closing base pairs. The free energy increments (Delta Gdegrees(loop)) of the 2 x 2 loops were positive, and showed relatively small differences between the two closing base pair environments. Assuming temperature-independent enthalpy changes for the transitions, (Delta Gdegrees(loop)) for the 2 x 2 loops varied from 0.9 to 1.9 kcal/mol in 1 M Na(+) at 37 degrees C. Most values were within 0.8 kcal/mol of previously estimated values; however, a few sequences differed by 1.2-2.0 kcal/mol. Single strands employed to form the RNA duplexes exhibited small noncooperative absorbance increases with temperature or transitions indicative of partial self-complementary duplexes. One strand formed a partial self-complementary duplex that was more stable than the tandem mismatch duplexes it formed. Transitions of the RNA duplexes were analyzed using equations that included the coupled equilibrium of self-complementary duplex and non-self-complementary duplex denaturation. The average heat capacity change (DeltaC(p)) associated with the transitions of two RNA duplexes was estimated by plotting DeltaH degrees and DeltaS degrees evaluated at different strand concentrations as a function of T(m) and ln T(m), respectively. The average DeltaC(p) was 70 +/- 5 cal K(-)(1) (mol of base pairs)(-)(1). Consideration of this heat capacity change reduced the free energy of formation at 37 degrees C of the 10 bp control RNA duplexes by 0.3-0.6 kcal/mol, which may increase Delta Gdegrees(loop) values by similar amounts.     

Parallel-stranded duplex DNA containing blocks of trans purine-purine and purine-pyrimidine base pairs.

A 30 base pair parallel-stranded (ps) duplex ps-L1.L2 composed of two adjoined purine-purine and purine-pyrimidine sequence blocks has been characterized thermodynamically and spectroscopically. The 5'-terminal 15 residues in both strands ('left-half') consisted of the alternating d(GA)7G sequence that forms a ps homoduplex secondary structure stabilized by d(G.G) and d(A.A) base pairs. The 3'-terminal 15 positions of the sequence ('right-half') were combinations of A and T with complementary reverse Watson-Crick d(A.T) base pairing between the two strands. The characteristics of the full length duplex were compared to those of the constituent left and right halves in order to determine the compatibility of the two ps helical forms. The thermal denaturation curves and hyperchromicity profiles of all three duplexes determined by UV absorption spectroscopy were characteristic of ps-DNA, in accordance with previous studies. The thermodynamic properties of the 30 bp duplex corresponded within experimental error to the linear combination of the two 15-mers. Thus, the Tm and delta HvH of ps-L1.L2 in 10 mM MgCl2, derived from analyses according to a statistical mechanical formulation for the helix-coil transition, were 43 degrees C and 569 kJ mol-1, compared to 21 degrees C, 315 kJ mol-1 (ps-F5.F6) and 22 degrees C, 236 kJ mol-1 (ps-GA15). The UV absorption and CD spectra of ps-L1.L2 and the individual 15-mer ps motifs were also compared quantitatively. The sums of the two constituent native spectra (left+right halves) accurately matched that of the 30 bp duplex, with only small deviations in the 195-215 nm (CD) and 220-240 nm (absorption) regions. Based on analysis by native gel electrophoresis, the sequences studied formed duplex structures exclusively; there were no indications of higher order species. Chemical modification with diethyl pyrocarbonate showed no hyperreactivity of the junctional bases, indicating a smooth transition between the two parallel-stranded conformations. We conclude that under given salt conditions, oligonucleotides with normal primary chemical structures can readily form a parallel-stranded double helix based on blocks of very disparate non-canonical purine-purine and purine-pyrimidine base pairs and without perceptible destabilization at the junction. There are biological implications of these findings in relation to genetic structure and expression.     

Crystal structure of the highly distorted chimeric decamer r(C)d(CGGCGCCG)r(G).spermine complex--spermine binding to phosphate only and minor groove tertiary base-pairing.

The crystal structure of the self-complementary chimeric decamer duplex r(C)d(CGGCGCCG)r(G), with RNA base pairs at both termini, has been solved at 1.9 A resolution by the molecular replacement method and refined to an R value of 0.145 for 2,314 reflections. The C3'-endo sugar puckers of the terminal riboses apparently drive the entire chimeric duplex into an A-DNA conformation, in contrast to the B-DNA conformation adopted by the all-deoxy decamer of the same sequence. Five symmetry related duplexes encapsulate a spermine molecule which interacts with ten phosphate groups, both directly and through water molecules to form multiple ionic and hydrogen bonding interactions. The spermine interaction severely bends the duplexes by 31 degrees into the major groove at the fourth base pair G(4).C(17), jolts it and slides the 'base plate' into the minor groove. This base pair, together with the adjacent base pair in the top half and the corresponding pseudo two-fold related base pairs in the bottom half, form four minor groove base-paired multiples with the terminal base pairs of two neighboring duplexes.     

Thermodynamic and structural properties of the specific binding between Ag ion and C:C mismatched base pair in duplex DNA to form C–Ag–C metal-mediated base pair

Metal ion-nucleic acid interactions have attracted considerable interest for their involvement in structure formation and catalytic activity of nucleic acids. Although interactions between metal ion and mismatched base pair duplex are important to understand mechanism of gene mutations related to heavy metal ions, they have not been well-characterized. We recently found that the Ag⁺ ion stabilized a C:C mismatched base pair duplex DNA. A C–Ag–C metal-mediated base pair was supposed to be formed by the binding between the Ag⁺ ion and the C:C mismatched base pair to stabilize the duplex. Here, we examined specificity, thermodynamics and structure of possible C–Ag–C metal-mediated base pair. UV melting indicated that only the duplex with the C:C mismatched base pair, and not of the duplexes with the perfectly matched and other mismatched base pairs, was specifically stabilized on adding the Ag⁺ ion. Isothermal titration calorimetry demonstrated that the Ag⁺ ion specifically bound with the C:C base pair at 1:1 molar ratio with a binding constant of 10⁶ M⁻¹, which was significantly larger than those for nonspecific metal ion-DNA interactions. Electrospray ionization mass spectrometry also supported the specific 1:1 binding between the Ag⁺ ion and the C:C base pair. Circular dichroism spectroscopy and NMR revealed that the Ag⁺ ion may bind with the N3 positions of the C:C base pair without distorting the higher-order structure of the duplex. We conclude that the specific formation of C–Ag–C base pair with large binding affinity would provide a binding mode of metal ion-DNA interactions, similar to that of the previously reported T-Hg-T base pair. The C–Ag–C base pair may be useful not only for understanding of molecular mechanism of gene mutations related to heavy metal ions but also for wide variety of potential applications of metal-mediated base pairs in various fields, such as material, life and environmental sciences.     

HU binding to bent DNA: a fluorescence resonance energy transfer and anisotropy study

HU, an architectural DNA-binding protein, either stabilizes DNA in a bent conformation or induces a bend upon binding to give other proteins access to the DNA. In this study, HU binding affinity for a bent DNA sequence relative to a linear sequence was investigated using fluorescence anisotropy measurements. A static bend was achieved by the introduction of two phased A4T4 tracts in a 20 bp duplex. Binding affinity for 20 bp duplexes containing two phased A-tracts in either a 5'-3' or 3'-5' orientation was found to be almost 10-fold higher than HU binding to a random sequence 20 bp duplex (6.1 vs 0.68 microM(-1)). The fluorescence technique of resonance energy transfer was used to quantitatively determine the static bend of the DNA duplexes and the HU-induced bend. DNA molecules were 5'-end labeled with fluorescein as the donor or rhodamine as the acceptor. From the efficiency of energy transfer, the end-to-end distance of the DNA duplexes was calculated. The end-to-end distance relative to DNA contour length (R/R(C)) yields a bend angle for the A-tract duplex of 45 +/- 7 degrees in the absence of HU and 70 +/- 3 degrees in the presence of HU. The bend angle calculated for the T4A4 tract duplex was 62 +/- 4 degrees after binding two HU dimers. Fluorescence anisotropy measurements reveal that HU binds in a 1:1 stoichiometry to the A4T4 tract duplex but a 2:1 stoichiometry to the T4A4 tract and random sequence duplex. These findings suggest that HU binding and recognition of DNA may be governed by a structural mechanism.     

The absence of non-local conformational changes in OR3 operator DNA on complexing with the cro repressor

The Interaction of the cro protein of lambda phage with a synthetic OR3 operator having 17 base pairs in length and with its 9 bp fragment has been studied using the circular dichroism (CD) method. In both cases, a considerable change in the CD of the samples was found in the region 260-300 nm upon the addition of the cro protein. The stoichiometry obtained by the CD titration was identical for OR3 and its 9 bp fragment: one duplex per dimeric cro. NaCl addition makes the complexes dissociate so that the 9 bp fragment becomes free at [NaCl] greater than 0.2 M while the whole OR3 becomes free at [NaCl] greater than 0.5 M. The CD spectra of both the free duplexes show a typical B-form conservative pattern with a positive CD band (270 nm) and a negative one (250 nm). The specific complexing of both the duplexes results in a substantial CD depression in the positive band. The most pronounced effect occurs at 280 nm. This spectral change is quite distinct from those in the B to A transition and in the non-cooperative winding of the DNA within the B-family of forms. The interaction of the cro protein with the non-operator DNAs, calf thymus DNA and a synthetic 10 bp duplex, reveals no visible CD changes at all. An inference is drawn that the CD change in the specific complexes is mainly due to the induced CD in tyr-26 upon its interaction with a specific base pair in the operator or its fragment, the operator DNA conformation being conserved in a B-like form as a whole.(ABSTRACT TRUNCATED AT 250 WORDS)     

HU binding to DNA: evidence for multiple complex formation and DNA bending

HU, a nonspecific histone-like DNA binding protein, participates in a number of genomic events as an accessory protein and forms multiple complexes with DNA. The HU-DNA binding interaction was characterized by fluorescence, generated with the guanosine analogue 3-methyl-8-(2-deoxy-beta-D-ribofuranosyl)isoxanthopterin (3-MI) directly incorporated into DNA duplexes. The stoichiometry and equilibrium binding constants of complexes formed between HU and 13 and 34 bp DNA duplexes were determined using fluorescence anisotropy and analytical ultracentrifugation. These measurements reveal that three HU molecules bind to the 34 bp duplexes, while two HU molecules bind to the 13 bp duplex. The data are well described by an independent binding site model, and the association constants for the first binding event for both duplexes are similar (approximately 1 x 10(6) M(-1)), indicating that HU binding affinity is independent of duplex length. Further analysis of the binding curves in terms of a nonspecific binding model is indicative that HU binding to DNA exhibits little to no cooperativity. The fluorescence intensity also increases upon HU binding, consistent with decreased base stacking and increased solvent exposure of the 3-MI fluorescence probe. These results are suggestive of a local bending or unwinding of the DNA. On the basis of these results we propose a model in which bending of DNA accompanies HU binding. Up to five complex bands are observed in gel mobility shift assays of HU binding to the 34 bp duplexes. We suggest that protein-induced bending of the DNA leads to the observation of complexes in the gel, which have the same molecular weight but different relative mobilities.     

Hidden in plain sight: subtle effects of the 8-oxoguanine lesion on the structure, dynamics, and thermodynamics of a 15-base pair oligodeoxynucleotide duplex

The base lesion 8-oxoguanine is formed readily by oxidation of DNA, potentially leading to G → T transversion mutations. Despite the apparent similarity of 8-oxoguanine-cytosine base pairs to normal guanine-cytosine base pairs, cellular base excision repair systems effectively recognize the lesion base. Here we apply several techniques to examine a single 8-oxoguanine lesion at the center of a nonpalindromic 15-mer duplex oligonucleotide in an effort to determine what, if anything, distinguishes an 8-oxoguanine-cytosine (8oxoG-C) base pair from a normal base pair. The lesion duplex is globally almost indistinguishable from the unmodified parent duplex using circular dichroism spectroscopy and ultraviolet melting thermodynamics. The DNA mismatch-detecting photocleavage agent Rh(bpy)(2)chrysi(3+) cleaves only weakly and nonspecifically, revealing that the 8oxoG-C pair is locally stable at the level of the individual base pairs. Nuclear magnetic resonance spectra are also consistent with a well-conserved B-form duplex structure. In the two-dimensional nuclear Overhauser effect spectra, base-sugar and imino-imino cross-peaks are strikingly similar between parent and lesion duplexes. Changes in chemical shift due to the 8oxoG lesion are localized to its complementary cytosine and to the 2-3 bp immediately flanking the lesion on the lesion strand. Residues further removed from the lesion are shown to be unperturbed by its presence. Notably, imino exchange experiments indicate that the 8-oxoguanine-cytosine pair is strong and stable, with an apparent equilibrium constant for opening equal to that of other internal guanine-cytosine base pairs, on the order of 10(-6). This collection of experiments shows that the 8-oxoguanine-cytosine base pair is incredibly stable and similar to the native pair.     

Correcting for heat capacity and 5'-TA type terminal nearest neighbors improves prediction of DNA melting temperatures using nearest-neighbor thermodynamic models

Nearest-neighbor thermodynamic (NNT) models currently provide some of the most accurate predictions of melting thermodynamics, including melting temperature (T(m)) values, for short DNA duplexes. Inherent to all existing NNT models is the assumption that ΔH° and ΔS° for the helix-to-coil transition are temperature invariant. Here we investigate the impact that this zero-ΔC(p) assumption has on the accuracy of T(m) predictions for 128 DNA duplexes. Previous and new melting thermodynamic data are analyzed to establish an estimate of ΔC(p)(bp), the heat capacity change per base pair, of 42 ± 16 cal mol(-1) K(-1) bp(-1), as well as an optimal thermodynamic reference temperature (T(ref)) of 53 ± 5 °C. These results were used to modify the unified NNT model to properly account for the temperature dependence of ΔH° and ΔS° and thereby extend the range over which T(m) is accurately predicted. This new approach is shown to be especially useful for duplexes that melt at a T(m) greater than 70 °C. Thermodynamic data collected by differential scanning calorimetry (DSC) for 16 duplexes designed to melt over a broad temperature range were used to verify the values of ΔC(p)(bp) and T(ref) and to show that ΔC(p)(bp) is essentially constant above 37 °C. Additional DSC analysis of 12 duplex sequences containing all 10 nearest neighbors allowed for errors associated with different terminal nearest neighbors to be examined and showed that duplexes containing one or more terminal 5'-TA groups are significantly more stable than predicted by the unified NNT model. A correction to improve prediction of the hybridization thermodynamics of duplexes with terminal 5'-TA groups is provided.     

Thermodynamics of sequence-specific binding of PNA to DNA

For further characterization of the hybridization properties of peptide nucleic acids (PNAs), the thermodynamics of hybridization of mixed sequence PNA-DNA duplexes have been studied. We have characterized the binding of PNA to DNA in terms of binding affinity (perfectly matched duplexes) and sequence specificity of binding (singly mismatched duplexes) using mainly absorption hypochromicity melting curves and isothermal titration calorimetry. For perfectly sequence-matched duplexes of varying lengths (6-20 bp), the average free energy of binding (DeltaG degrees ) was determined to be -6.5+/-0.3 kJ mol(-1) bp(-1), corresponding to a microscopic binding constant of about 14 M(-1) bp(-1). A variety of single mismatches were introduced in 9- and 12-mer PNA-DNA duplexes. Melting temperatures (T(m)) of 9- and 12-mer PNA-DNA duplexes with a single mismatch dropped typically 15-20 degrees C relative to that of the perfectly matched sequence with a corresponding free energy penalty of about 15 kJ mol(-1) bp(-1). The average cost of a single mismatch is therefore estimated to be on the order of or larger than the gain of two matched base pairs, resulting in an apparent binding constant of only 0.02 M(-1) per mismatch. The impact of a mismatch was found to be dependent on the neighboring base pairs. To a first approximation, increasing the stability of the surrounding region, i.e., the distribution of A.T and G.C base pairs, decreases the effect of the introduced mismatch.     

Thermodynamic, spectroscopic, and equilibrium binding studies of DNA sequence context effects in four 40 base pair deoxyoligonucleotides

Effects of different end sequences on melting, circular dichroism spectra (CD), and enzyme binding properties were investigated for four 40 base pair, non-self-complementary duplex DNA oligomers. The center sequences of these oligoduplexes have either of two 22 base pair modules flanked on both sides by sequences differing in AT content. Temperature-induced melting transitions monitored by differential scanning calorimetry (DSC) and ultraviolet absorbance were measured for the six duplexes in buffered 115 mM Na(+) solutions. Values of the melting transition enthalpy, DeltaH(cal), and entropy, DeltaS(cal), were obtained directly from DSC experiments. Melting transition parameters, DeltaH(vH) and DeltaS(vH), were also estimated from a van't Hoff analysis of optical melting curves collected as a function of DNA concentration, assuming that the melting transition is two-state. Melting free energies (20 degrees C) evaluated from DSC melting experiments on the four duplex DNAs ranged from -52.2 to -77.5 kcal/mol. Free energies based on the van't Hoff analysis were -37.9 to -58.8 kcal/mol. Although the values are different, trends in the melting free energies of the four duplex DNAs as a function of sequence were identical in both DSC and optical analyses. Subject to several assumptions, values for the initiation free energy were estimated for each duplex, defined as DeltaG(int) = DeltaG(cal) - DeltaG(pred), where DeltaG(cal) is the experimental free energy at 20 degrees C determined from the experimentially measured values of the transition enthalpy, DeltaH(cal), and entropy, DeltaS(cal). The predicted free energy of the sequence, DeltaG(pred)(20 degrees C), is obtained using published nearest-neighbor sequence stability values. For three of the four duplexes, values of DeltaG(int) are essentially nil. In contrast, the duplex with 81.8% GC has a considerably higher estimate of DeltaG(int) = 7.1 kcal/mol. The CD spectra for the six duplexes collected over the wavelength range from 200 to 320 nm are also sequence-dependent. Factor analysis of the CD spectra by singular value decomposition revealed that the experimental CD spectra could be reconstructed from linear combinations of two minor and one major subspectra. Changes in the coefficients of the major subspectrum for different sequences reflect incremental sequence-dependent variations of the CD spectra. Equilibrium binding by BamHI restriction endonuclease to the 40 base pair DNAs whose central eight base pairs contain the recognition sequence for BamHI restriction enzyme bounded by A.T base pairs, 5'-A-GGATCC-A-3' was investigated. Binding assays were performed by titering BamHI against a constant concentration of each of the duplex DNA substrates, in the absence of Mg(2+), followed by analysis by gel retardation. Under the conditions employed, the enzyme binds but does not cleave the DNAs. Results of the assays revealed two binding modes with retarded gel mobilities. Binding isotherms for the fraction of bound DNA species versus enzyme concentration for each binding mode were constructed and analyzed with a simple two-step equilibrium binding model. This analysis provided semiquantitative estimates on the equilibrium binding constants for BamHI to the four DNAs. Values obtained for the binding constants varied only 7-fold and ranged from 6 x 10(-)(8) to 42 x 10(-)(8) M, with binding free energies from -8.6 to -9.7 (+/- 0.2) kcal/mol depending on the sequence that flanks the enzyme binding site. Unlike what was found earlier in binding studies of the 22 base pair duplexes that constitute the core modules of the present 40-mers [Riccelli, P. V., Vallone, P. M., Kashin, I., Faldasz, B. D., Lane, M. J., and Benight, A. S. (1999) Biochemistry 38, 11197-11208], no obvious relationship between binding and stability was found for these longer DNAs. Apparently, effects of flanking sequence stability on restriction enzyme binding may only be measurable in very short duplex deoxyoligonucl     

Unusual structural features of hydantoin lesions translate into efficient recognition by Escherichia coli Fpg

Oxidation of guanine (G) and 8-oxoguanine (OG) with a wide variety of oxidants yields the hydantoin lesions, guanidinohydantoin (Gh) and spiroiminodihydantoin (Sp). These two lesions have garnered much recent attention due to their unusual structures and high mutagenic potential. We have previously shown that duplexes containing Gh and Sp are substrates for the base excision repair glycosylase Escherichia coli Fpg (EcFpg). To evaluate the recognition features of these unusual lesions, binding and footprinting experiments were performed using a glycosylase inactive variant, E3Q EcFpg, and 30 bp duplexes containing the embedded lesions. Surprisingly, E3Q EcFpg was found to bind significantly more tightly ( approximately 1000-fold) to duplexes containing Gh or Sp over the corresponding duplexes containing OG. This may be a consequence of the helix-destabilizing nature of the hydantoin lesions that facilitates their recognition within duplex DNA. Though DNA binding affinities of E3Q EcFpg with Gh- and Sp-containing duplexes were found to be similar to each other, hydroxyl radical footprinting using methidium-propyl-EDTA (MPE)-Fe(II) revealed subtle differences between binding of E3Q EcFpg to the two lesions. Most notably, in the presence of E3Q EcFpg, the Sp nucleotide (nt) is hyperreactive toward cleavage by MPE-Fe(II)-generated hydroxyl radicals, suggestive of the formation of an intercalation site for the MPE-Fe(II) reagent at the Sp nt. Interestingly, increasing the duplex length from 18 to 30 bp enhanced the excision efficiency of Gh and Sp paired with C, G, or T by EcFpg such that these substrates are processed as efficiently as the signature substrate lesion, OG. Moreover, the base removal activity with these two lesions was more efficient than removal of OG when in a base pairing context opposite A. The high affinity and efficient activity of EcFpg toward the hydantoin lesions suggest that EcFpg mediates repair of the lesions in vivo. Notably, the facile activity of EcFpg toward Gh and Sp in base pairing contexts with G and A, which are likely to be present after DNA replication, would be detrimental and enhance mutagenesis.     

T4-DNA ligase: substrate properties of synthetic DNA-duplexes with structural anomalies

Structural defects, affecting T4 DNA ligase function, were revealed with the help of synthetic DNA duplexes, containing modifications at single nick. Changes of configuration at C2' and C3' atoms of furanose in the acceptor terminus lead to total blocking of the nick sealing activity of T4 DNA ligase. On the contrary, substitution of 3'-terminal deoxyribonucleotide for ribonucleotide doesn't affect the enzyme's action. The duplex looses all of it's substrate activity if the next from the nick G.C pair is substituted for the noncomplementary G.C pair. In DNA duplexes containing an unpaired base in the nick, elimination of the extrahelical nucleotide proceeds the ligation step. In these cases the duplex substrate activity decreases depending on the extent of extrahelical base stacking into the double stranded DNA.     

Differential hydration of dA.dT base pairing and dA and dT bulges in deoxyoligonucleotides

The role of water in the formation of stable duplexes of nucleic acids is being studied by determining the concurrent volume change, heats, and counterion uptake that accompany the duplexation process. The variability of the volume contraction that we have observed in the formation of a variety of homoduplexes suggests that sequence and conformation acutely affect the degree of hydration. We have used a combination of densimetric and calorimetric techniques to measure the change in volume and enthalpy resulting from the mixing of two complementary strands to form (a) fully paired duplexes with 10 or 11 base pairs and (b) bulged decameric duplexes with an extra dA or dT unmatched residue. We also monitored absorbance vs temperature profiles as a function of strand and salt concentration for all four duplexes. Relative to the decamer duplex, insertion of an extra dA.dT base pair to form an undecamer duplex results in a favorable enthalpy of -5.6 kcal/mol that is nearly compensated by an unfavorable entropy term of -5.1 kcal/mol. This enthalpy difference correlates with a differential uptake of water molecules, corresponding to an additional hydration of 16 mol of water molecules/mol of base pair. Relative to the fully paired duplexes, both bulged duplexes are 12-16 degrees C less stable and exhibit marginally larger counterion uptake on forming the duplex. The enthalpy change is slightly lower for the T-bulge duplex and less still for the A-bulge duplex. The volume change results indicate that an unmatched residue increases the amount of coulombic and/or structural hydration. The combined results strongly suggest that the destabilizing forces in bulged duplexes are partially compensated by an increase in hydration levels.     

Effects of internal nonbonded bases and a G.U base pair on the stability of a short ribonucleic acid helix.

Proton nuclear magnetic resonance has been used to examine the effect of both noncomplementary and G.U oppositions in the duplexes formed by the synthetic pentaribonucleotides CpApApUpG, CpApUpUpG, CpApGpUpG, and CpApCpUpG. The lack of any sigmoidal behavior in the chemical shift vs. temperature plots of the base protons in the individual pentaribonucleotides indicates that duplexes with noncomplementary base oppositions of the type: formula: (see text), (where X = A, U, G, or C) do not form. Variable temperature spectra of the mixture of CpApGpUpG and CpApUpUpG were recorded over the range of 70--10 degrees C. The chemical shift vs. temperature plot of the purine aromatic protons displayed sigmoidal curves. This demonstrated both duplex formation and the presence of a G.U. base pair. The average Tm of the duplex was found to be 23.4 +/- 2.0 degrees C. This is similar to that of the duplex formed by CpApUpG (24.0 +/- 1.0 degrees C) but less than the Tm of the following duplexes: CpApApUpG:CpApUpUpG (Tm = 28.5 +/- 2.1 degrees C), CpApGpUpG:CpApCpUpG (Tm = 38.4 +/- 0.6 degrees C) and CpApUpApUpG (Tm = 41.5 +/- 1.1 degrees C). The G.U base pair has a Tm (20.0 degrees C) significantly lower than the rest of the duplex (24 +/- 1 degree C) and is a region of local instability within the double helix. This 1H NMR study is the first to investigate both the formation and relative stability of an internal G.U. base pair neighboring regular Watson--Crick base pairs.     

Thermodynamics and kinetics for base-pair opening in the P1 duplex of the Tetrahymena group I ribozyme

The thermodynamics and kinetics for base-pair opening of the P1 duplex of the Tetrahymena group I ribozyme were studied by NMR hydrogen exchange experiments. The apparent equilibrium constants for base pair opening were measured for most of the imino protons in the P1 duplex using the base catalysts NH3, HPO4(2-) or TRIS. These equilibrium constants were also measured for several modified P1 duplexes, and the C-2.G23 base pair was the most stable base pair in all the duplexes. The conserved U-1*G22 base pair is required for activity of the ribozyme and the data here show that this wobble base pair destabilizes neighboring base pairs on only one side of the wobble. A 2'-OMe modification on the U-3 residue stabilized its own base pair but had little effect on the neighboring base pairs. Three base pairs, U-1*G22, C-2*G23 and A2*U21 showed unusual equilibrium constants for opening and possible implications of the opening thermodynamics of these base pairs on the undocking rates of the P1 helix with catalytic core are discussed.     

Conformational features of benzo‐homologated yDNA duplexes by molecular dynamics simulation

yDNA is a base‐modified nucleic acid duplex containing size‐expanded nucleobases. Base‐modified nucleic acids could expand the genetic alphabet and thereby enhance the functional potential of DNA. Unrestrained 100 ns MD simulations were performed in explicit solvent on the yDNA NMR sequence [5′(yA T yA yA T yA T T yA T)₂] and two modeled yDNA duplexes, [5′(yC yC G yC yC G G yC G G)₂] and [(yT5′ G yT A yC yG C yA yG T3′)?(yA5′ C T C yG C G yT A yC A3′)]. The force field parameters for the yDNA bases were derived in consistent with the well‐established AMBER force field. Our results show that DNA backbone can withstand the stretched size of the bases retaining the Watson‐Crick base pairing in the duplexes. The duplexes retained their double helical structure throughout the simulations accommodating the strain due to expanded bases in the backbone torsion angles, sugar pucker and helical parameters. The effect of the benzo‐expansion is clearly reflected in the extended C1′‐C1′ distances and enlarged groove widths. The size expanded base modification leads to reduction in base pair twist resulting in larger overlapping area between the stacked bases, enhancing inter and intra strand stacking interactions in yDNA in comparison with BDNA. This geometry could favour enhanced interactions with the groove binders and DNA binding proteins., 2016. © 2015 Wiley Periodicals, Inc. Biopolymers 105: 55–64, 2016