An NMR study of polymorphous behaviour of the mismatched DNA octamer d(m5C-G-m5C-G-A-G-m5C-G) in solution. The B-duplex and hairpin forms

Effect of benzyl alcohol on recombinant human interleukin-1 receptor antagonist structure and hydrogen-deuterium exchange

Benzyl alcohol, a preservative commonly added to multidose therapeutic protein formulations, can accelerate aggregation of recombinant human interleukin-1 receptor antagonist (rhIL-1ra). To investigate the interactions between benzyl alcohol and rhIL-1ra, we used nuclear magnetic resonance to observe the effect of benzyl alcohol on the chemical shifts of amide resonances of rhIL-1ra and to measure hydrogen-deuterium exchange rates of individual rhIL-1ra residues. Addition of 0.9% benzyl alcohol caused significant chemical shifts of amide resonances for residues 90-97, suggesting that these solvent-exposed residues participate in the binding of benzyl alcohol. In contrast, little perturbation of exchange rates was observed in the presence of either sucrose or benzyl alcohol.

Probing microsecond time scale dynamics in proteins by methyl (1)H Carr-Purcell-Meiboom-Gill relaxation dispersion NMR measurements. Application to activation of the signaling protein NtrC(r)

To study microsecond processes by relaxation dispersion NMR spectroscopy, low power deposition and short pulses are crucial and encourage the development of experiments that employ (1)H Carr-Purcell-Meiboom-Gill (CPMG) pulse trains. Herein, a method is described for the comprehensive study of microsecond to millisecond time scale dynamics of methyl groups in proteins, exploiting their high abundance and favorable relaxation properties. In our approach, protein samples are produced using [(1)H, (13)C]-d-glucose in ∼100% D(2)O, which yields CHD(2) methyl groups for alanine, valine, threonine, isoleucine, leucine, and methionine residues with high abundance, in an otherwise largely deuterated background. Methyl groups in such samples can be sequence-specifically assigned to near completion, using (13)C TOCSY NMR spectroscopy, as was recently demonstrated (Otten, R.; et al. J. Am. Chem. Soc. 2010, 132, 2952-2960). In this Article, NMR pulse schemes are presented to measure (1)H CPMG relaxation dispersion profiles for CHD(2) methyl groups, in a vein similar to that of backbone relaxation experiments. Because of the high deuteration level of methyl-bearing side chains, artifacts arising from proton scalar coupling during the CPMG pulse train are negligible, with the exception of Ile-δ1 and Thr-γ2 methyl groups, and a pulse scheme is described to remove the artifacts for those residues. Strong (13)C scalar coupling effects, observed for several leucine residues, are removed by alternative biochemical and NMR approaches. The methodology is applied to the transcriptional activator NtrC(r), for which an inactive/active state transition was previously measured and the motions in the microsecond time range were estimated through a combination of backbone (15)N CPMG dispersion NMR spectroscopy and a collection of experiments to determine the exchange-free component to the transverse relaxation rate. Exchange contributions to the (1)H line width were detected for 21 methyl groups, and these probes were found to collectively report on a local structural rearrangement around the phosphorylation site, with a rate constant of (15.5 ± 0.5) × 10(3) per second (i.e., τ(ex) = 64.7 ± 1.9 μs). The affected methyl groups indicate that, already before phosphorylation, a substantial, transient rearrangement takes place between helices 3 and 4 and strands 4 and 5. This conformational equilibrium allows the protein to gain access to the active, signaling state in the absence of covalent modification through a shift in a pre-existing dynamic equilibrium. Moreover, the conformational switching maps exactly to the regions that differ between the solution NMR structures of the fully inactive and active states. These results demonstrate that a cost-effective and quantitative study of protein methyl group dynamics by (1)H CPMG relaxation dispersion NMR spectroscopy is possible and can be applied to study functional motions on the microsecond time scale that cannot be accessed by backbone (15)N relaxation dispersion NMR. The use of methyl groups as dynamics probes extends such applications also to larger proteins.

Hexamer oligonucleotide topology and assembly under solution phase NMR and theoretical modeling scrutiny

The entire family of noncomplementary hexamer oligodeoxyribonucleotides d(GCXYGC) (X and Y = A, G, C, or T) were assessed for topological indicators and equilibrium thermodynamics using a priori molecular modeling and solution phase NMR spectroscopy. Feasible modeled hairpin structures formed a basis from which solution structure and equilibria for each oligonucleotide were considered. ¹H and ³¹P variable temperature-dependent (VT) and concentration-dependent NMR data, NMR signal assignments, and diffusion parameters led to d(GCGAGC) and d(GCGGGC) being understood as exceptions within the family in terms of self-association and topological character. A mean diffusion coefficient D(298 K) = (2.0 ± 0.07) × 10⁻¹⁰ m² s⁻¹ was evaluated across all hexamers except for d(GCGAGC) (D(298 K) = 1.7 × 10⁻¹⁰ m² s⁻¹) and d(GCGGGC) (D(298 K) = 1.2 × 10⁻¹⁰ m² s⁻¹). Melting under VT analysis (T(m) = 323 K) combined with supporting NMR evidence confirmed d(GCGAGC) as the shortest tandem sheared GA mismatched duplex. Diffusion measurements were used to conclude that d(GCGGGC) preferentially exists as the shortest stable quadruplex structure. Thermodynamic analysis of all data led to the assertion that, with the exception of XY = GA and GG, the remaining noncomplementary oligonucleotides adopt equilibria between monomer and duplex, contributed largely by monomer random-coil forms. Contrastingly, d(GCGAGC) showed preference for tandem sheared GA mismatch duplex formation with an association constant K = 3.9 × 10⁵M⁻¹. No direct evidence was acquired for hairpin formation in any instance although its potential existence is considered possible for d(GCGAGC) on the basis of molecular modeling studies.

Chemical-shift referencing and resolution stability in gradient LC-NMR (acetonitrile:water)

Wide ranges of solvent conditions are generated during solvent-gradient LC-NMR. This complicates the referencing of the chemical-shift scale of the resulting NMR data. The problems that arise when performing LC-NMR in acetonitrile:water - particularly when the mobile-phase composition can range anywhere from 0% to 100% - are examined here, and the reliability of the secondary reference signals are evaluated. It is shown that under these conditions the use of the acetonitrile signal is superior to the use of the water signal in any form (either the (1)H or the (2)H signal) as a secondary reference, a lock signal, and a signal for shimming. The limitations of the referencing methods and other experimental parameters, and the limitations on the solvent-gradient ramp parameters, primarily as they affect lineshapes, are all shown. These results are compared to the way some other publications have referenced the (1)H chemical-shift axis (when using acetonitrile:water mixtures to perform reversed-phase chromatography LC-NMR).

Structural Evaluation of a Novel Pro-apoptotic Peptide Coupled to CNGRC Tumor Homing Sequence by NMR

Hunter-killer peptides (HKPs) are synthetic peptides that target specific cell types for apoptosis. These studies report functional and structural characteristics of HKP9, an hunter-killer peptide that specifically targets tumor vasculature with a new apoptotic sequence. Vesicle leakage experiments were performed as a model for membrane perturbing activity. Placement of the homing sequence reduces both cell toxicity and vesicle leakage activity. NMR studies elucidate the conformation and orientation of HKP9 in micelles. The positively charged end of the HKP9 killing sequence is solvent exposed; however, the central portion of the peptide is helical and buried in dodecylphosphorylcholine micelles. The homing sequence is less solvent exposed than in a previously reported tumor-homing peptide. The results suggest that solvent accessibility of the homing sequence should be considered in design of future peptides.

Calcium-modulated S100 protein-phospholipid interactions. An NMR study of calbindin D9k and DPC

The cellular functions of several S100 proteins involve specific interactions with phospholipids and the cell membrane. The interactions between calbindin D(9k) (S100D) and the detergent dodecyl phosphocholine (DPC) were studied using NMR spectroscopy. In the absence of Ca(2+), the protein associates with DPC micelles. The micelle-associated state has intact helical secondary structures but no apparent tertiary fold. At neutral pH, Ca(2+)-loaded calbindin D(9k) does not associate with DPC micelles. However, a specific interaction is observed with individual DPC molecules at a site close to the linker between the two EF-hands. Binding to this site occurs only when Ca(2+) is bound to the protein. A reduction in pH in the absence of Ca(2+) increases the stability of the micelle-associated state. This along with the corresponding reduction in Ca(2+) affinity causes a transition to the micelle-associated state also in the presence of Ca(2+) when the pH is lowered. Site-specific analysis of the data indicates that calbindin D(9k) has a core of three tightly packed helices (A, B, and D), with a dynamic fourth helix (C) more loosely associated. Evidence is presented that the Ca(2+)-binding characteristics of the two EF-hands are distinctly different in a micelle environment. The role of calbindin D(9k) in the cell is discussed, along with the broader implications for the function of the S100 protein family.

Orientation and helical conformation of a tissue-specific hunter-killer peptide in micelles

Hunter-killer peptides are chimeric synthetic peptides that selectively target specific cell types for an apoptotic death. These peptides, which are models for potential therapeutics, contain a homing sequence for receptor-mediated interactions and a pro-apoptotic sequence. Homing domains have been designed to target angiogenic tumor cells, prostate cells, arthritic tissue and, most recently, adipose tissue. After a receptor-mediated internalization, the apoptotic sequence, which contains D-enantiomer amino acids, initiates apoptosis through mitochondrial membrane disruption. We have begun structure and functional studies on a peptide (HKP1) that specifically targets angiogenic tumor cells for apoptosis. As a model for mitochondrial membrane disruption, we have examined peptide-induced leakage of a calcein fluorophore from large unilamellar vesicles. These experiments demonstrate more potent leakage activity by HKP1 than the peptide lacking the homing domain. Circular dichroism and 2D homonuclear NMR experiments demonstrate that this tumor-specific HKP adopts a left-handed amphipathic helix in association with dodecylphosphorylcholine micelles in a parallel orientation to the lipid-water interface with the homing domain remaining exposed to solvent. The amphipathic helix of the apoptotic domain orients with nonpolar leucine and alanine residues inserting most deeply into the lipid environment.

The depletion of protein signals in metabonomics analysis with the WET-CPMG pulse sequence

Nuclear magnetic resonance (NMR) spectroscopy is a powerful analytical tool capable of providing a comprehensive metabolic profile of biofluids such as urine, plasma, and serum. Unfortunately, when measuring serum and plasma, the high protein concentration can obscure the signals originating from low molecular weight metabolites. We evaluated the use of different parameters within the Carr-Purcell-Meiboom-Gill (CPMG) pulse train of fast spin-echoes to remove the macromolecular signal contribution in one-dimensional proton (1H) NMR spectra. Experimental parameters such as the refocusing delay in the CPMG pulse train, pulse miscalibration, and recycle time were examined to assess the ability to remove the protein signals from the spectrum without causing a deleterious effect on the signals originating from free, low molecular weight metabolites. The 1H-NMR spectra of a variety of serum samples spiked with 2'-deoxyadenosine were acquired using various acquisition parameters. Our results show that the delay used in the CPMG spin-echo and the combination of the acquisition pulse flip angle and recycle time are the two major factors affecting the observed metabolite signal amplitudes in the resulting 1H-NMR spectrum.

Association between the first two immunoglobulin-like domains of the neural cell adhesion molecule N-CAM

The extracellular domain of N-CAM contains five immunoglobulin-like (Ig) and two fibronectin type III-like domains and facilitates cell-cell binding through multiple, weak interdomain interactions. NMR spectroscopy indicated that the two N-terminal Ig-like domains from chicken N-CAM (Ig I and Ig II) interact with millimolar affinity. Physico-chemical studies show that this interaction is significantly amplified when the domains are covalently linked, consistent with an antiparallel domain arrangement. The binding of the two individual domains and the dimerization of the concatenated protein were essentially independent of salt, up to a concentration of 200 mM. The residues in Ig I involved in the interaction map to the BED strands of the beta sandwich, and delineate a largely hydrophobic patch.

DNA-induced conformational changes are the basis for cooperative dimerization by the DNA binding domain of the retinoid X receptor

Dimerization of the DNA-binding domains of nuclear hormone receptors occurs in a manner that is highly cooperative with DNA binding. We have investigated the molecular basis for this cooperativity through an NMR study of the interaction between the monomeric DNA-binding domain (DBD) of the retinoid-X-receptor (RXR) and a single DNA half-site. Major changes were observed in the chemical shifts of the backbone resonances and in the pattern of medium-range nuclear Overhauser enhancement connectivities of the RXR upon binding to DNA, indicating that the DNA induces conformational changes in the monomer. Binding to DNA induces and stabilizes the structure in a region of the second zinc binding domain that forms the dimerization interface when RXR binds as a dimer to a direct repeat recognition element. These studies provide direct experimental evidence that DNA-induced protein conformational changes constitute the molecular basis for cooperative enhancement of dimer formation and DNA binding by the nuclear hormone receptor DBDs. In contrast to the localized folding induced in the dimerization interface, DNA binding leads to unfolding of the C-terminal helix found in the free RXR DBD. Unwinding of this helix may facilitate homodimer formation by maximizing interactions between the two DNA-bound RXR domains.

High-resolution solution structure of the retinoid X receptor DNA-binding domain

The retinoid X receptor (RXR) is a member of the nuclear hormone receptor superfamily of transcriptional regulators and plays a central role in the retinoid and, through its ability to heterodimerize with other nuclear hormone receptors, non-steroid signaling pathways. The DNA-binding and recognition functions of RXR are located in a conserved 83 amino acid residue domain that recognizes the consensus sequence AGGTCA. In order to provide a detailed picture of its structure, we have calculated a high-resolution solution structure of the C195A RXRalpha DNA-binding domain. Structures were calculated using 1131 distance and dihedral angle constraints derived from 1H, 13C and 15N NMR spectra. The structures reveal a perpendicularly packed, "loop-helix" fold similar to other nuclear hormone receptor DNA-binding domains and confirm the existence of the C-terminal helix, which was first observed in the low-resolution NMR structure. The C-terminal helix is well formed and is stabilized by packing interactions with residues in the hydrophobic core. The solution structure of RXR is very similar to that determined by X-ray crystallographic studies of the RXR-TR heterodimer complex with DNA, except that in the latter case no electron density was observed for residues corresponding to the C-terminal helix. Other differences between the X-ray and NMR structures occur in the second zinc-binding loop, which is disordered in solution. Heteronuclear 15N NOE measurements suggest that this loop has enhanced flexibility in the free protein.

Kinetics of DNA hydration

The hydration of the d(CGCGAATTCGCG) B-DNA duplex in solution was studied by nuclear magnetic relaxation dispersion (NMRD) of the water nuclei 1H, 2H, and 17O, and by nuclear Overhauser effects (NOEs) in high-resolution two-dimensional 1H NMR spectra. By comparing results from the free duplex with those from its complex with netropsin, water molecules in the "spine of hydration" in the AATT region of the minor groove could be distinguished from hydration water elsewhere in the duplex. The 2H and 17O relaxation dispersions yield a model-independent residence time of 0.9(+/-0.1) ns at 4 degrees C for five highly ordered water molecules in the spine. When corrected for frequency offset effects, the NOE data yield the same residence time as the NMRD data, giving credence to both methods. At 27 degrees C, the residence time is estimated to 0.2 ns, a factor of 40 shorter than the tumbling time of the duplex. The NMRD data show that all water molecules associated with the duplex, except the five molecules in the spine, have residence times significantly shorter than 1 ns at 4 degrees C. There is thus no long-lived hydration structure associated with the phosphate backbone. In contrast to 2H and 17O, the 1H relaxation dispersion is dominated by labile DNA protons and therefore provides little information about DNA hydration.

A single G-to-C change causes human centromere TGGAA repeats to fold back into hairpins

Recently, we established that satellite III (TGGAA)n tandem repeats, which occur at the centromeres of human chromosomes, pair with themselves to form an unusual "self-complementary" antiparallel duplex containing (GGA)2 motifs in which two unpaired guanines from opposite strands intercalate between sheared G.A base pairs. In separate studies, we have also established that the GCA triplet does not form bimolecular (GCA)2 motifs but instead promotes the formation of hairpins containing a GCA-turn motif in which the loop contains a single cytidine closed by a sheared G.A pair. Since TGCAA is the most frequent variant of TGGAA found in satellite III repeats, we reasoned that the potential of this variant to form GCA-turn miniloop fold-back structures might be an important factor in modulating the local structure in natural (TGGAA)n repeats. We report here the NMR-derived solution structure of the heptadecadeoxynucleotide (G)TGGAATGCAATGGAA(C) in which a central TGCAA pentamer is flanked by two TGGAA pentamers. This 17-mer forms a rather unusual and very stable hairpin structure containing eight base pairs in the stem, only four of which are Watson-Crick pairs, and a loop consisting of a single cytidine residue. The stem contains a (GGA)2 motif with intercalative 14G/4G stacking between two sheared G.A base pairs; the loop end of the stem consists of a sheared 8G.10A closing pair with the cytosine base of the 9C loop stacked on 8G. The remarkable stability of this unusual hairpin structure (Tm = 63 degrees C) suggests that it probably plays an important role in modulating the folding of satellite III (TGGAA)n repeats at the centromere.

Slow conformational exchange in DNA minihairpin loops: a conformational study of the circular dumbbell d<pCGC-TT-GCG-TT>

In recent years various examples of highly stable two-residue hairpin loops (miniloops) in DNA have been encountered. As the detailed structure and stability of miniloops appear to be determined not only by the nature and sequence of the two bases in the loop, but also by the closing base pair, it is desirable to carry out in-depth studies of especially designed small model DNA compounds. Therefore, a circular DNA dumbbell-like molecule is tailored to consist of a stem of three Watson-Crick base pairs, flanked on each side by a minihairpin loop. The resulting circular DNA decamer 5'-d<pCGC-TT-GCG-TT>-3' (I) is studied in solution by means of nmr spectroscopy. At a temperature of 269 K the molecule occurs in a 50/50 mixture of two dumbbell structures (denoted L2L2 and L2L4). L2L2 contains three Watson-Crick C-G base pairs and two two-residue loops (H2-family type) in opposite parts of the molecule. On raising the temperature from 269 to 314 K, the L2L4 conformer becomes increasingly dominant (95% at 314 K). This conformer has a partially disrupted closing G-C base pair in the 5'-GTTC-3' loop with only one remaining solvent-accessible hydrogen bond between NH alpha of the cytosine C(1) and O6 of the guanine G(8), whereas the opposite 5'-CTTG-3' loop remains stable. The disruption of the C(1)-G(8) base pair in the L2L4 form is correlated with the presence of a syn orientation for the C(1) base at the 5'-3' loop-stem junction in the 5'-GTTC-3' loop. The two conformers, L2L2 and L2L4, occur in slow equilibrium (2-20 s-1). Moderate line broadening of specific 1H, 13C, and 31P resonances of residues C(1), G(8), T(9), and T(10) at low temperatures, due to chemical exchange between L2L2 and L2L4, show that the interconversion from an anti to syn conformer in residue C(1) has a small local effect on the structure of the dumbbell. T1 relaxation measurements, chemical-shift considerations, and complete band-shape calculations of the exchange process of the G(8) imino proton reveal a possibility for the existence of multiconformational states in the anti-syn equilibrium.

Hairpin opening by single-strand-specific nucleases

DNA molecules with covalently sealed (hairpin) ends are probable intermediates in V(D)J recombination. According to current models hairpin ends are opened to produce short single-stranded extensions that are thought to be precursors of a particular type of extra nucleotides, termed P nucleotides, which are frequently present at recombination junctions. Nothing is known about the activities responsible for hairpin opening. We have used two single-strand-specific nucleases to explore the effects of loop sequence on the hairpin opening reaction. Here we show that a variety of hairpin ends are opened by P1 nuclease and mung bean nuclease (MBN) to leave short, 1-2 nt single-stranded extensions. Analysis of 22 different hairpin sequences demonstrates that the terminal 4 nt of the hairpin loop strongly influence the sites of cleavage. Correlation of the nuclease digestion patterns with structural (NMR) data for some of the hairpin loops studied here provides new insights into the structural features recognized by these enzymes.

Oligonucleotides containing G.A pairs: effect of flanking sequences on structure and stability

Sixteen oligodeoxyribonucleotides, 5'd(GXGAYC)3', X and Y = G, A, C, or T, have been synthesized and studied by UV melting and 1H and 31P NMR methods. By varying X and Y, the sixteen resulting oligonucleotides can theoretically form 10 duplexes with all possible Watson-Crick base pairs flanking the center two G.A base pairs. Two-dimensional 1H NMR data on 5'd(GCGAGC)3' revealed that the center bases G and A pair through G amino hydrogen bonding and that the two consecutive G.A pairs form excellent purine-purine stacks. The concurrent appearance of one or more upfield-shifted imino proton peaks (approximately 10.5 ppm) and both upfield- and downfield- shifted 31P signals (approximately -2 and approximately -5.1 ppm) was a unique characteristic in imino 1H and 31P NMR spectra and was used as a conformational probe for this type of G.A pairs. Using this probe, seven out of 10 duplexes of 5'GXGAYC3' were found to adopt the G.A base pairing with G amino proton bonding and G to G and A to A base stacking. Three were in the group comprising 5'pyrimidine-GA-purine3', and four were in the group comprising 5'purine-GA-purine3'/5'pyrimidine-GA-pyrimidine3'. The G.A pairs in 5'purine-GA-pyrimidine3' adopted a totally different conformation. Thermodynamic analysis indicated that duplexes in the 5'pyrimidine-GA-purine3' group were more stable than the duplexes in the 5'purine-GA-pyrimidine3' group. Overall, G.C base pairs were preferred as neighbors to this type of G.A pairs.

Sequence dependent effects of CpG cytosine methylation. A joint 1H-NMR and 31P-NMR study

The impact of cytosine methylation in the central CpG step of two closely related octanucleotide duplexes d(CATCGATG)2 and d(CTTCGAAG)2 was examined by 1H-NMR and 31P-NMR experiments, and a quantitative structural analysis was performed using the NOE-derived distances, the sugar puckers and the epsilon torsion angles. The two starting oligonucleotides displayed a B-DNA conformation with, however, significant local structure differences. Although the methylated oligonucleotides retained their B-DNA conformation, different structural and thermal stability effects were observed. The magnitude of the methylation effects was to depend on the initial conformation of the CpG site, which is governed by the nature of the dinucleotide AT or TT located on the CpG flanks. As an example of sequence dependence, the methylation of CpG entailed larger conformational variation in d(CATCGATG)2 than in d(CTTCGAAG)2. In this study, the 1H and 31P chemical-shift parameters averred as extremely sensitive probes for detecting subtle conformational changes. Finally, our comparative results may aid our understanding of the structural and related biological effects produced by cytosine methylation in DNA.

Structural characterization of d(CAACCCGTTG) and d(CAACGGGTTG) mini-hairpin loops by heteronuclear NMR: the effects of purines versus pyrimidines in DNA hairpins

The DNA decamers, d(CAACCCGTTG) and d(CAACGGGTTG) were studied in solution by proton and heteronuclear NMR. Under appropriate conditions of pH, temperature, salt concentration and DNA concentration, both decamers form hairpin conformations with similar stabilities [Avizonis and Kearns (1995) Biopolymers, 35, 187-200]. Both decamers adopt mini-hairpin loops, where the first and last four nucleotides are involved in Watson-Crick hydrogen bonding and the central two nucleotides, CC or GG respectively, form the loop. Through the use of proton-proton, proton-phosphorus and natural abundance proton-carbon NMR experiments, backbone torsion angles (beta, gamma and epsilon), sugar puckers and interproton distances were measured. The nucleotides forming the loops of these decamers were found to stack upon one another in an L1 type of loop conformation. Both show gamma tr and unusual beta torsion angles in the loop-closing nucleotide G7, as expected for mini-hairpin loop formation. Our results indicate that the beta and epsilon torsion angles of the fifth and sixth nucleotides that form the loop and the loop-closing nucleotide G7 are not in the standard trans conformation as found in B-DNA. Although the loop structures calculated from NMR-derived constraints are not well defined, the stacking of the bases in the two different hairpins is different. This difference in the base stacking of the loop may provide an explanation as to why the cytosine-containing hairpin is thermodynamically more stable than the guanine-containing hairpin.

Conformational studies of hairpin sequences from the ColE1 cruciform

Inverted repeat sequences derived from the ColE1 cruciform were investigated by nuclear magnetic resonance (NMR) and UV spectroscopy. It was shown that 15 different sequences exist as stable hairpin structures over a range of buffer conditions and DNA concentrations. Experiments with six oligomers (1-6) containing the native stem sequence and five base loops, found that the two hairpins with the wild-type loops (1-2) served as upper and lower bounds for the thermodynamic stability of all the other sequences. NMR experiments, including rotational correlation time measurements and NOESY spectra, were then performed on 1, the most stable hairpin sequence to begin to uncover a structural basis of its stability.

An NMR study of the conformation and thermodynamics of the circular dumbbell d [formula: see text] Slow exchange between two- and four-membered hairpin loops

The circular DNA decamer 5'-d [formula: see text] 3' is studied in solution by means of NMR spectroscopy. At low temperature the molecule adopts a dumbbell structure with three Watson-Crick C-G base pairs and two two-residue loops in opposite parts of the molecule. On raising the temperature another conformer appears, in which the closing C-G base pair in the 5'-GTTC-3' loop is disrupted, whereas the opposite 5'-CTTG-3' loop remains stable. The two conformers are in slow equilibrium over a limited temperature range.

Conformational consequences of the incorporation of arabinofuranosylcytidine in DNA. An NMR study of the DNA fragments d(CGCTAGCG) and d(CGaCTAGCG) in solution

The self-complementary octamers d(CGCTAGCG) and d(CGaCTAGCG) (aC, arabinofuranosylcytidine) were studied by means of NMR spectroscopy. It is shown that d(CGaCTAGCG), under suitable conditions of oligonucleotide concentration, ionic strength and temperature, exclusively adopts a hairpin structure. However, under the same experimental conditions (5 mM DNA, no added salt, 295 K) d(CGCTAGCG) mainly adopts a B-DNA-type duplex. At lower temperatures (less than or equal to 290 K) the hairpin form of d(CGaCTAGCG) occurs in slow exchange with an intact B-DNA-type duplex. When the DNA concentration of d(CGCTAGCG) is dramatically reduced (less than or equal to 0.5 mM) the hairpin form becomes highly favoured at the expense of the dimer. Moreover, proton-chemical-shift considerations indicate that the structural features of the hairpin structure of d(CGCTAGCG) mimic, in part, those of the modified octamer d(CGaCTAGCG), i.e. a loop comprising only the two central residues with the thymine located into the minor groove (Pieters, J. M. L., de Vroom, E., van der Marel, G. A., van Boom, J. H., Koning, T. M. G., Kaptein, R. and Altona, C. unpublished results). Thermodynamic analysis of d(CGCTAGCG) yields an average Tmd value of 342 K (1 M DNA) and a delta Hod value of -266 kJ/mol for the dimer/coil transition and an average Tmh value of 321 K and delta Hoh - 102 kJ/mol for the hairpin/coil equilibrium. For the duplex/coil equilibrium of d(CGaCTAGCG) an average Tmd value of 336 K (1 M DNA) and delta Hod value of -253 kJ/mol are deduced. The hairpin/coil transition of d(CGaCTAGCG) is characterized by a delta Hoh value of -104 kJ/mol and an average Tmh value of 331 K. It is concluded that incorporation of an arabinofuranosylcytidine in the octamer d(CGaCTAGCG) results in stabilization of the hairpin form, whereas the dimer is destablized by two aC.dG base pairs.

Thermodynamics of the various forms of the dodecamer d(ATTACCGGTAAT) and of its constituent hexamers from proton nmr chemical shifts and UV melting curves: three-state and four-state thermodynamic models

Chemical shifts of base and H1' protons of the single-stranded hexamers d(ATTACC) and d(GGTAAT), of the 1:1 mixtures of these complementary hexamers, and of the self-complementary dodecamer d(ATTACCGGTAAT) were measured at various temperatures in aqueous solution. Four different sample concentrations were used in the case of the dodecamer and of the mixture of the complementary hexamers; the individual hexamers were measured at two different DNA concentrations. Absorbance temperature profiles at five different NaCl concentrations were measured for the dodecamer in order to quantify the effect of the ionic strength on the duplex formation. Under suitable conditions of nucleotide concentration, temperature, and ionic strength, the dodecamer adopts either a B-DNA duplex or a hairpin-loop structure. Chemical shift vs temperature profiles, constructed for all samples, were used to obtain thermodynamic parameters either for the various stacking interactions in the single strands or for the duplex or the hairpin-loop formation. In the analysis of the duplex formation of the hexamers, a two-state approach appeared too simple, because systematic deviations were revealed. Therefore, a new three-state model (DUPSTAK) was developed. In order to investigate the magnitude of error arising from the use of the two-state approach in cases where the DUPSTAK model appears more appropriate, a series of test calculations was made. The magnitude of error in the enthalpy and in the entropy of duplex melting is found to depend linearly upon the actual melting temperature and not upon the individual delta Hd degrees and delta Sd degrees values. Thermodynamic analysis of the chemical shift vs temperature profiles in D2O solution (no added salt) yields an average Tmd value of 341 K (1M DNA) and delta Hd degrees of - 121 kJ.mol-1 for the dimer/random-coil transition of the hexamer duplex d(ATTACC).d(GGTAAT). For the duplex in equilibrium random-coil transition of the 12-mer d(ATTACCGGTAAT) an average Tmd value of 336 K (1M DNA) and delta Hd degrees of -372 kJ.mol-1 are found. The hairpin/random-coil transition of d(ATTACCGGTAAT) is characterized by a rather large delta Hh degrees value, -130 kJ.mol-1, and an average Tmh value of 304 K.

Unusual duplex formation in purine rich oligodeoxyribonucleotides

The purine rich oligodeoxyribonucleotides 1C, d(ATGACGGAATA) and 2C, d(ATGAGCGAATA) alone exhibit highly cooperative melting transitions. Analysis of the concentration dependence of melting, and electrophoretic studies indicate that these oligomers can form an unusual purine rich offset double helix. The unusual duplex is predicted to contain four A.T, two G.C, and four G.A mismatch base pairs as well as a single A base stacked on the 3' end of each chain of the helix. Other possible models for the duplex are unlikely because they are predicted to contain many base pairs of low stability. Changing the central sequence to CGG or GGG should destabilize the duplex and this is observed. The unusual duplex of 2C is more stable than the duplex of 1C indicating that the stability of G.A base pairs is quite sensitive to the surrounding sequence. Addition of 1C and 2C to their complementary pyrimidine strands results in normal duplexes of similar stability. We feel that the unusual duplexes are significantly stabilized by the intrinsic stacking tendency of purine bases.