Extracting thermodynamic data from equilibrium melting curves for oligonucleotide order-disorder transitions

Developing Community Resources for Nucleic Acid Structures

In this review, we describe the creation of the Nucleic Acid Database (NDB) at Rutgers University and how it became a testbed for the current infrastructure of the RCSB Protein Data Bank. We describe some of the special features of the NDB and how it has been used to enable research. Plans for the next phase as the Nucleic Acid Knowledgebase (NAKB) are summarized.

Vicinal abasic site impaired processing of a Tg:G mismatch and 8-oxoguanine lesions in three-component bistranded clustered DNA damage

The occurrence of 7,8-dihydro-8-oxo-2′deoxyguanosine (8-oxodG), thymine glycol:guanine (Tg:G) mismatch and abasic site DNA damage lesions in close proximity induce repair refractive multicomponent clustered DNA damage. Herein, the influence of abasic sites in the processing of 8-oxodG lesion and Tg:G mismatch bistranded cluster is evaluated. Abasic sites are found to impart conformational destabilization that appreciably hinders the repair activity of the other lesions whenever present in a cluster combination. The repair process reduces the formation of double strand breaks (DSBs) and renders this three-lesion combination a non-DSB forming cluster. The stability of the DNA duplex harbouring these three lesions is highly compromised due to altered base helicity and base stacking phenomena leading to impaired repair.

The occurrence of 7,8-dihydro-8-oxo-2′deoxyguanosine (8-oxodG), thymine glycol:guanine (Tg:G) mismatch and abasic site DNA damage lesions in close proximity induce repair refractive non-DSB cluster.

Versatile Method for the Site-Specific Modification of DNA with Boron Clusters: Anti-Epidermal Growth Factor Receptor (EGFR) Antisense Oligonucleotide Case

A general and convenient approach for the incorporation of different types of boron clusters into specific locations of the DNA-oligonucleotide chain based on the automated phosphoramidite method of oligonucleotide synthesis and post-synthetic "click chemistry" modification has been developed. Pronounced effects of boron-cluster modification on the physico- and biochemical properties of the antisense oligonucleotides were observed. The silencing activity of antisense oligonucleotides bearing a single boron cluster modification in the middle of the oligonucleotide chain was substantially higher than that of unmodified oligonucleotides. This finding may be of importance for the design of therapeutic nucleic acids with improved properties. The proposed synthetic methodology broadens the availability of nucleic acid-boron cluster conjugates and opens up new avenues for their potential practical use.

Mode of Action of a Designed Antimicrobial Peptide: High Potency against Cryptococcus neoformans

There is a significant need for developing compounds that kill Cryptococcus neoformans, the fungal pathogen that causes meningoencephalitis in immunocompromised individuals. Here, we report the mode of action of a designed antifungal peptide, VG16KRKP (VARGWKRKCPLFGKGG) against C. neoformans. It is shown that VG16KRKP kills fungal cells mainly through membrane compromise leading to efflux of ions and cell metabolites. Intracellular localization, inhibition of in vitro transcription, and DNA binding suggest a secondary mode of action for the peptide, hinting at possible intracellular targets. Atomistic structure of the peptide determined by NMR experiments on live C. neoformans cells reveals an amphipathic arrangement stabilized by hydrophobic interactions among A2, W5, and F12, a conventional folding pattern also known to play a major role in peptide-mediated Gram-negative bacterial killing, revealing the importance of this motif. These structural details in the context of live cell provide valuable insights into the design of potent peptides for effective treatment of human and plant fungal infections.

Sequence-dependent theory of oligonucleotide hybridization kinetics

A theoretical approach to the prediction of the sequence and temperature-dependent rate constants for oligonucleotide hybridization reactions has been developed based on the theory of relaxation kinetics. One-sided and two-sided melting reaction mechanisms for oligonucleotide hybridization reactions have been considered, analyzed, modified, and compared to select a physically consistent as well as robust model for prediction of the relaxation times of DNA hybridization reactions that agrees with the experimental evidence. The temperature- and sequence-dependent parameters of the proposed model have been estimated using available experimental data. The relaxation time model that we developed has been combined with the nearest neighbor model of hybridization thermodynamics to estimate the temperature- and sequence-dependent rate constants of an oligonucleotide hybridization reaction. The model-predicted rate constants are compared to experimentally determined rate constants for the same oligonucleotide hybridization reactions. Finally, we consider a few important applications of kinetically controlled DNA hybridization reactions.

APE1 incision activity at abasic sites in tandem repeat sequences

Repetitive DNA sequences, such as those present in microsatellites and minisatellites, telomeres, and trinucleotide repeats (linked to fragile X syndrome, Huntington disease, etc.), account for nearly 30% of the human genome. These domains exhibit enhanced susceptibility to oxidative attack to yield base modifications, strand breaks, and abasic sites; have a propensity to adopt non-canonical DNA forms modulated by the positions of the lesions; and, when not properly processed, can contribute to genome instability that underlies aging and disease development. Knowledge on the repair efficiencies of DNA damage within such repetitive sequences is therefore crucial for understanding the impact of such domains on genomic integrity. In the present study, using strategically designed oligonucleotide substrates, we determined the ability of human apurinic/apyrimidinic endonuclease 1 (APE1) to cleave at apurinic/apyrimidinic (AP) sites in a collection of tandem DNA repeat landscapes involving telomeric and CAG/CTG repeat sequences. Our studies reveal the differential influence of domain sequence, conformation, and AP site location/relative positioning on the efficiency of APE1 binding and strand incision. Intriguingly, our data demonstrate that APE1 endonuclease efficiency correlates with the thermodynamic stability of the DNA substrate. We discuss how these results have both predictive and mechanistic consequences for understanding the success and failure of repair protein activity associated with such oxidatively sensitive, conformationally plastic/dynamic repetitive DNA domains.

Impact of bulge loop size on DNA triplet repeat domains: Implications for DNA repair and expansion

Repetitive DNA sequences exhibit complex structural and energy landscapes, populated by metastable, noncanonical states, that favor expansion and deletion events correlated with disease phenotypes. To probe the origins of such genotype-phenotype linkages, we report the impact of sequence and repeat number on properties of (CNG) repeat bulge loops. We find the stability of duplexes with a repeat bulge loop is controlled by two opposing effects; a loop junction-dependent destabilization of the underlying double helix, and a self-structure dependent stabilization of the repeat bulge loop. For small bulge loops, destabilization of the underlying double helix overwhelms any favorable contribution from loop self-structure. As bulge loop size increases, the stabilizing loop structure contribution dominates. The role of sequence on repeat loop stability can be understood in terms of its impact on the opposing influences of junction formation and loop structure. The nature of the bulge loop affects the thermodynamics of these two contributions differently, resulting in unique differences in repeat size-dependent minima in the overall enthalpy, entropy, and free energy changes. Our results define factors that control repeat bulge loop formation; knowledge required to understand how this helix imperfection is linked to DNA expansion, deletion, and disease phenotypes.

Molecular crowding inhibits U-insertion/deletion RNA editing in vitro: consequences for the in vivo reaction

Mitochondrial pre-mRNAs in African trypanosomes are edited to generate functional transcripts. The reaction is typified by the insertion and deletion of U nucleotides and is catalyzed by a macromolecular complex, the editosome. Editosomes bind pre-edited mRNA/gRNA pairs and the reaction can be recapitulated in vitro by using pre-mRNA- and gRNA-mimicking oligoribonucleotides together with enriched editosome preparations. Although the in vitro assay has been instrumental in unraveling the basic steps of the editing cycle it is performed at dilute solvent conditions. This ignores the fact that editing takes place inside the highly crowded mitochondria. Here we investigate the effects of molecular crowding on RNA editing. By using neutral, macromolecular cosolutes we generate defined dilute, semidilute and crowded solvent properties and we demonstrate different thermodynamic stabilities of the pre-mRNA/gRNA hybrid RNAs at these conditions. Crowded conditions stabilize the RNAs by -30 kJ/mol. Furthermore, we show that the rate constants for the association and dissociation (kass/kdiss) of substrate RNAs to editosomes decrease, ultimately inhibiting the in vitro reaction. The data demonstrate that the current RNA editing in vitro system is sensitive to molecular crowding, which suggests that the in vivo reaction cannot rely on a diffusion-controlled, collision-based mechanism. Possible non-diffusional reaction pathways are discussed.

Thermodynamic consequences of the hyperoxidized guanine lesion guanidinohydantoin in duplex DNA

Guanidinohydantoin (Gh) is a hyperoxidized DNA lesion produced by oxidation of 8-oxo-7,8-dihydroguanine (8-oxoG). Previous work has shown that Gh is potently mutagenic in both in vitro and in vivo coding for G → T and G → C transversion mutations. In this work, analysis by circular dichroism shows that the Gh lesion does not significantly alter the global structure of a 15-mer duplex and that the DNA remains in the B-form. However, we find that Gh causes a large decrease in the thermal stability, decreasing the duplex melting temperature by ~17 °C relative to an unmodified duplex control. Using optical melting analysis and differential scanning calorimetry, the thermodynamic parameters describing duplex melting were also determined. We find that the Gh lesion causes a dramatic decrease in the enthalpic stability of the duplex. This enthalpic destabilization is somewhat tempered by entropic stabilization; yet, Gh results in an overall decrease in thermodynamic stability of the duplex relative to a control that lacks DNA damage, with a ΔΔG° of -7 kcal/mol. These results contribute to our understanding of the consequences of hyperoxidation of G and provide insight into how the thermal and thermodynamic destabilization caused by Gh may influence replication and/or repair of the lesion.

Effects of locked nucleic acid substitutions on the stability of oligonucleotide hairpins

An understanding of the stability of nucleic acid folding is critical for applications involving RNA viruses, small molecule-RNA binding, and therapeutics, for example. To explore factors that affect this stability, hairpins made from oligonucleotides containing both a GAAA tetraloop and three to five complements in the stem have been used as models where locked nucleic acids (LNAs) have been substituted into the sequence. UV spectroscopy was used to obtain melting curves in 20% by volume formamide, and the enthalpies and entropies of melting were determined. Although LNA substitutions typically increase the stability of a hybrid, we have found a decrease in stability for DNA and RNA GAAA hairpins when LNA is substituted into the loop. Tetraloops synthesized from natural bases show higher enthalpies and entropies of melting compared to the LNA substituted sequences indicating that LNA substitutions can destabilize a hairpin but stabilize the corresponding double stranded structure.

Energy landscapes of dynamic ensembles of rolling triplet repeat bulge loops: implications for DNA expansion associated with disease states

DNA repeat domains can form ensembles of canonical and noncanonical states, including stable and metastable DNA secondary structures. Such sequence-induced structural diversity creates complex conformational landscapes for DNA processing pathways, including those triplet expansion events that accompany replication, recombination, and/or repair. Here we demonstrate further levels of conformational complexity within repeat domains. Specifically, we show that bulge loop structures within an extended repeat domain can form dynamic ensembles containing a distribution of loop positions, thereby yielding families of positional loop isomers, which we designate as "rollamers". Our fluorescence, absorbance, and calorimetric data are consistent with loop migration/translocation between sites within the repeat domain ("rollamerization"). We demonstrate that such "rollameric" migration of bulge loops within repeat sequences can invade and disrupt previously formed base-paired domains via an isoenthalpic, entropy-driven process. We further demonstrate that destabilizing abasic lesions alter the loop distributions so as to favor "rollamers" with the lesion positioned at the duplex/loop junction, sites where the flexibility of the abasic "universal hinge" relaxes unfavorable interactions and/or facilitates topological accommodation. Another strategic siting of an abasic site induces directed loop migration toward denaturing domains, a phenomenon that merges destabilizing domains. In the aggregate, our data reveal that dynamic ensembles within repeat domains profoundly impact the overall energetics of such DNA constructs as well as the distribution of states by which they denature/renature. These static and dynamic influences within triplet repeat domains expand the conformational space available for selection and targeting by the DNA processing machinery. We propose that such dynamic ensembles and their associated impact on DNA properties influence pathways that lead to DNA expansion.

Structural properties of DNA oligomers containing (GACX)(n) and (GAXC)(n) tandem repeats

The antisense DNA sequence of mature mouse micro RNA, miR341, includes three repeats of the tetranucleotide (GACC). The -GAC- repeat is known to form a parallel duplex, in acidic environments. The thermal melting profile of miR341 DNA, at pH 4, 5, and 6 indicates the formation of a very stable structure, which loses its stability when pH is increased. Thus, the addition of a cytosine at the 3' end of the (GAC) motif preserves the molecule's potential to fold into an unusual structure at low pH. The effect of modifying the nucleotide composition of the GACC sequence on the secondary structures formed by oligomers containing seven tandem repeats of the altered motifs was examined here. UV melting profiles were determined, as a function of pH, for 28-mers of the two series (GAXC)(7) and (GACX)(7) (X= A/C/T/G)(.) The sequence (GACC)(7) was found to be extremely sensitive to pH variations, with a stable structure formed at pH 5 (T(m) ≥ 60°C). NMR spectroscopy established that the low pH structure is not B-DNA. (GACA)(7) and (GACT)(7) also formed stable structures at low pH but the addition of guanine at the 3'end, as seen in the (GACG) series resulted in the loss of this property. Introducing a break in the 5'-GAC-3' motif, explored in the (GAXC) series, also inhibits formation of stable structures under acidic conditions.

Interstrand dipole-dipole interactions can stabilize the collagen triple helix

The amino acid sequence of collagen is composed of GlyXaaYaa repeats. A prevailing paradigm maintains that stable collagen triple helices form when (2S)-proline (Pro) or Pro derivatives that prefer the C(γ)-endo ring pucker are in the Xaa position and Pro derivatives that prefer the C(γ)-exo ring pucker are in the Yaa position. Anomalously, an amino acid sequence in an invertebrate collagen has (2S,4R)-4-hydroxyproline (Hyp), a C(γ)-exo-puckered Pro derivative, in the Xaa position. In certain contexts, triple helices with Hyp in the Xaa position are now known to be hyperstable. Most intriguingly, the sequence (GlyHypHyp)(n) forms a more stable triple helix than does the sequence (GlyProHyp)(n). Competing theories exist for the physicochemical basis of the hyperstability of (GlyHypHyp)(n) triple helices. By synthesizing and analyzing triple helices with different C(γ)-exo-puckered proline derivatives in the Xaa and Yaa positions, we conclude that interstrand dipole-dipole interactions are the primary determinant of their additional stability. These findings provide a new framework for understanding collagen stability.

Stability of intramolecular quadruplexes: sequence effects in the central loop

Hundreds of thousands of putative quadruplex sequences have been found in the human genome. It is important to understand the rules that govern the stability of these intramolecular structures. In this report, we analysed sequence effects in a 3-base-long central loop, keeping the rest of the quadruplex unchanged. A first series of 36 different sequences were compared; they correspond to the general formula GGGTTTGGGHNHGGGTTTGGG. One clear rule emerged from the comparison of all sequence motifs: the presence of an adenine at the first position of the loop was significantly detrimental to stability. In contrast, adenines have no detrimental effect when present at the second or third position of the loop. Cytosines may either have a stabilizing or destabilizing effect depending on their position. In general, the correlation between the T(m) or DeltaG degrees in sodium and potassium was weak. To determine if these sequence effects could be generalized to different quadruplexes, specific loops were tested in different sequence contexts. Analysis of 26 extra sequences confirmed the general destabilizing effect of adenine as the first base of the loop(s). Finally, analysis of some of the sequences by microcalorimetry (DSC) confirmed the differences found between the sequence motifs.

Divalent metal ions promote the formation of the 5'-splice site recognition complex in a self-splicing group II intron

Group II introns are ribozymes occurring in genes of plants, fungi, lower eukaryotes, and bacteria. These large RNA molecular machines, ranging in length from 400 to 2500 nucleotides, are able to catalyze their own excision from pre-mRNA, as well as to reinsert themselves into RNA or sometimes even DNA. The intronic domain 1 contains two sequences (exon binding sites 1 and 2, EBS1 and EBS2) that pair with their complementary regions at the 3'-end of the 5'-exon (intron binding sites 1 and 2, IBS1 and IBS2) such defining the 5'-splice site. The correct recognition of the 5'-splice site stands at the beginning of the two steps of splicing and is thus crucial for catalysis. It is known that metal ions play an important role in folding and catalysis of ribozymes in general. Here, we characterize the specific metal ion requirements for the formation of the 5'-splice site recognition complex from the mitochondrial yeast group II intron Sc.ai5gamma. Circular dichroism studies reveal that the formation of the EBS1.IBS1 duplex does not necessarily require divalent metal ions, as large amounts of monovalent metal ions also promote the duplex, albeit at a 5000 times higher concentration. Nevertheless, micromolar amounts of divalent metal ions, e.g. Mg2+ or Cd2+, strongly promote the formation of the 5'-splice site. These observations illustrate that a high charge density independent of the nature of the ion is needed for binding EBS1 to IBS1, but divalent metal ions are presumably the better players.

Electrical detection of the temperature induced melting transition of a DNA hairpin covalently attached to gold interdigitated microelectrodes

The temperature induced melting transition of a self-complementary DNA strand covalently attached at the 5' end to the surface of a gold interdigitated microelectrode (GIME) was monitored in a novel, label-free, manner. The structural state of the hairpin was assessed by measuring four different electronic properties of the GIME (capacitance, impedance, dissipation factor and phase angle) as a function of temperature from 25 degrees C to 80 degrees C. Consistent changes in all four electronic properties of the GIME were observed over this temperature range, and attributed to the transition of the attached single-stranded DNA (ssDNA) from an intramolecular, folded hairpin structure to a melted ssDNA. The melting curve of the self-complementary single strand was also measured in solution using differential scanning calorimetry (DSC) and UV absorbance spectroscopy. Temperature dependent electronic measurements on the surface and absorbance versus temperature values measured in solution experiments were analyzed assuming a two-state process. The model analysis provided estimates of the thermodynamic transition parameters of the hairpin on the surface. Two-state analyses of optical melting data and DSC measurements provided evaluations of the thermodynamic transition parameters of the hairpin in solution. Comparison of surface and solution measurements provided quantitative evaluation of the effect of the surface on the thermodynamics of the melting transition of the DNA hairpin.

Effects of the minor groove pyrimidine nucleobase functional groups on the stability of duplex DNA: The impact of uncompensated minor groove amino groups

DNA sequences containing four types of analog nucleosides are described. All four are pyridine derivatives constructed as C-nucleosides so that they mimic the pyrimidine derivatives 2'-deoxyuridine, thymidine or 2'-deoxycytidine, but in all cases the analogs lack the corresponding O2-carbonyls that in duplex DNA are located in the minor groove. In place of the O2-carbonyl is a hydrogen atom, a polar fluorine atom, or a nonpolar methyl group. The described C-nucleosides have native-like bidentate Watson-Crick hydrogen-bonding faces and can form essentially normal W-C base pairs of varying stability with A or G. In each modified base pair, two inter-residue hydrogen bonds should be present. In spite of a common number of interstrand hydrogen bonds, the thermodynamic stabilities of the prepared duplexes, each containing two analog base pairs, vary dramatically. Most notably, base pairs containing uncompensated purine amino groups (those lacking a hydrogen-bonding partner) in the minor groove exhibit the most dramatic reductions in thermodynamic stability. Removal of such uncompensated amino groups results in increased duplex stability. Base pairs containing fluorine in the minor groove positioned adjacent to an amino group seem to enhance duplex stability marginally (relative to --H or --CH(3)), but there is little evidence to suggest that fluorine is an effective hydrogen-bonding partner in these systems. The presence of minor groove methyl groups results in the least stable duplexes in each series of sequences.

PNA-DNA hybridization study using labeled streptavidin by voltammetry and surface plasmon fluorescence spectroscopy

Using ferrocene-streptavidin conjugates as amplifiers, we recently have demonstrated the simultaneous detection of DNA hybridization to peptide nucleic acid (PNA)-modified gold surfaces at the femtomole level by electrochemical and surface plasmon resonance techniques (Liu, J.; Tian, S.; Tiefenauer, L.; Nielsen, P. E.; Knoll, W. Anal. Chem. 2005, 77, 2756-2761). In this paper, a detailed study of the binding behavior of PNA-DNA is presented by square wave voltammetry and surface plasmon field-enhanced fluorescence spectroscopy (SPFS). The different binding constants for fully matched and single-mismatched DNA were obtained. The effect of the buffer concentration on the PNA-DNA hybrids was investigated using labeled streptavidin by cyclic voltammetry (CV) and SPFS. At high ionic strength, both the CV and SPFS signals were restrained dramatically, which is most probably due to a conformational change of the short-strand PNA-DNA helices on the surface. We conclude that the combination of electrochemical techniques with SPFS is very useful for the study of short DNA structure transformation.

Replication of 2-hydroxyadenine-containing DNA and recognition by human MutSalpha

2-Hydroxyadenine (2-OH-A), a product of DNA oxidation, is a potential source of mutations. We investigated how representative DNA polymerases from the A, B and Y families dealt with 2-OH-A in primer extension experiments. A template 2-OH-A reduced the rate of incorporation by DNA polymerase alpha (Pol alpha) and Klenow fragment (Kf(exo-)). Two Y family DNA polymerases, human polymerase eta (Pol eta) and the archeal Dpo4 polymerase were affected differently. Bypass by Pol eta was very inefficient whereas Dpo4 efficiently replicated 2-OH-A. Replication of a template 2-OH-A by both enzymes was mutagenic and caused base substitutions. Dpo4 additionally introduced single base deletions. Thermodynamic analysis showed that 2-OH-A forms stable base pairs with T, C and G, and to a lesser extent with A. Oligonucleotides containing 2-OH-A base pairs, including the preferred 2-OH-A:T, were recognized by the human MutSalpha mismatch repair (MMR). MutSalpha also recognized 2-OH-A located in a repeat sequence that mimics a frameshift intermediate.

Statistical thermodynamics and kinetics of DNA multiplex hybridization reactions

A general analytical description of the equilibrium and reaction kinetics of DNA multiplex hybridization has been developed. In this approach, multiplex hybridization is considered to be a competitive multichannel reaction process: a system wherein many species can react both specifically and nonspecifically with one another. General equations are presented that can consider equilibrium and kinetic models of multiplex hybridization systems comprised, in principle, of any number of targets and probes. Numerical solutions to these systems for both equilibrium and kinetic behaviors are provided. Practical examples demonstrate clear differences between results obtained from more common simplex methods, in which individual hybridization reactions are considered to occur in isolation; and multiplex hybridization, where desired and competitive cross-hybrid reactions between all possible pairs of strands are considered. In addition, sensitivities of the hybridization process of the perfect match duplex, to temperature, target concentration, and existence of sequence homology with other strands, are examined. This general approach also considers explicit sequence-dependent interactions between targets and probes involved in the reactions. Sequence-dependent stabilities of all perfect match and mismatch duplex complexes are explicitly considered and effects of relative stability of cross-hybrid complexes are also explored. Results reveal several interdependent factors that strongly influence DNA multiplex hybridization behavior. These include: relative concentrations of all probes and targets; relative thermodynamic stability of all perfect match and mismatch complexes; sensitivity to temperature, particularly for mismatches; and amount of sequence homology shared by the probe and target strands in the multiplex mix.

Analysis of thermal melting curves

T(m) is defined as Temperature of melting or, more accurately, as temperature of midtransition. This term is often used for nucleic acids (DNA and RNA, oligonucleotides and polynucleotides). A thermal denaturation experiment determines the stability of the secondary structure of a DNA or RNA and aids in the choice of the sequences for antisense oligomers or PCR primers. Beyond a simple numerical value (the T(m)), a thermal denaturation experiment, in which the folded fraction of a structure is plotted vs. temperature, yields important thermodynamic information. We present the classic problems encountered during these experiments and try to demonstrate that a number of useful pieces of information can be extracted from these experimental curves.

Fluorescence anisotropy in the frequency domain by an optical microscope

Fluorescence anisotropy decay spectroscopy is a suitable tool for investigating the size and the shape of biological molecules. We coupled this technique to an optical microscope in order to reduce the excitation volume and to allow its application to spatially inhomogeneous samples. Phase modulated measurements of the fluorescence anisotropy decay were performed by feeding an intensity modulated linearly polarized laser beam to the epifluorescence port of a microscope. Here we report the test of the dynamic response of the microscope by comparing the lifetime and fluorescence polarization anisotropy decays obtained in cuvettes in a standard phase modulation fluorometer and on tiny drops on the microscope stage. We show that once a correction factor for the objective depolarization is introduced in the best-fit functions for the data analysis of the decays, the results obtained on the two setups are comparable. Some applications are reported here on long DNA tracts as well on short DNA fragments containing structural anomalies.

EPR spectroscopic analysis of TAR RNA-metal ion interactions

Metal ion-induced changes in HIV-1 TAR RNA internal dynamics were determined by the changes in EPR spectral width for TAR RNAs containing spin-labeled nucleotides (U23, U25, U38, and U40). This gave a dynamic signature for each of 10 metal ions studied, which fell into one of three distinct groups. While Li(+) and K(+) had little effect on TAR RNA internal dynamics, Na(+) unexpectedly had a dynamic signature that was similar to Ca(2+) and Sr(2+), with a decrease in mobility at U23 and U38, little or no change at U25, and an increase in mobility at U40. Mg(2+), Co(2+), Ni(2+), Zn(2+), and Ba(2+) had similar effects on U23, U38, and U40, but the mobility of U25 was markedly increased. Our results show that RNA dynamics change upon metal binding to the TAR RNA bulge, indicating that RNA structure adapts to accommodate metal ions of different size and coordination properties.

Influence of DNA torsional rigidity on excision of 7,8-dihydro-8-oxo-2'-deoxyguanosine in the presence of opposing abasic sites by human OGG1 protein

The human protein OGG1 (hOGG1) targets the highly mutagenic base 7,8-dihydro-8-oxo-2'-deoxyguanosine (8-oxodG) and shows a high specificity for the opposite DNA base. Abasic sites can arise in DNA in close opposition to 8-oxodG either during repair of mismatched bases (i.e. 8-oxodG/A mismatches) or, more frequently, as a consequence of ionizing radiation exposure. Bistranded DNA lesions may remain unrepaired and lead to cell death via double-strand break formation. In order to explore the role of damaged-DNA dynamics in recognition/excision by the hOGG1 repair protein, specific oligonucleotides containing an 8-oxodG opposite an abasic site, at different relative distances on the complementary strand, were synthesized. Rotational dynamics were studied by means of fluorescence polarization anisotropy decay experiments and the torsional elastic constant as well as the hydrodynamic radius of the DNA fragments were evaluated. Efficiency of excision of 8-oxodG was tested using purified human glycosylase. A close relation between the twisting flexibility of the DNA fragment and the excision efficiency of the oxidative damage by hOGG1 protein within a cluster was found.

Conformational energetics of stable and metastable states formed by DNA triplet repeat oligonucleotides: implications for triplet expansion diseases

We have embedded the hexameric triplet repeats (CAG)(6) and (CTG)(6) between two (GC)(3) domains to produce two 30-mer hairpins with the sequences d[(GC)(3)(CAG)(6)(GC)(3)] and d[(GC)(3)(CTG)(6)(GC)(3)]. This construct reduces the conformational space available to these repetitive DNA sequences. We find that the (CAG)(6) and (CTG)(6) repeats form stable, ordered, single-stranded structures. These structures are stabilized at 62 degrees C by an average enthalpy per base of 1.38 kcal.mol(-1) for the CAG triplet and 2.87 kcal.mol(-1) for the CTG triplet, while being entropically destabilized by 3.50 cal.K(-1).mol(-1) for the CAG triplet and 7.6 cal.K(-1).mol(-1) for the CTG triplet. Remarkably, these values correspond, respectively, to 1/3 (for CAG) and 2/3 (for CTG) of the enthalpy and entropy per base values associated with Watson-Crick base pairs. We show that the presence of the loop structure kinetically inhibits duplex formation from the two complementary 30-mer hairpins, even though the duplex is the thermodynamically more stable state. Duplex formation, however, does occur at elevated temperatures. We propose that this thermally induced formation of a more stable duplex results from thermal disruption of the single-stranded order, thereby allowing the complementary domains to associate (perhaps via "kissing hairpins"). Our melting profiles show that, once duplex formation has occurred, the hairpin intermediate state cannot be reformed, consistent with our interpretation of kinetically trapped hairpin structures. The duplex formed by the two complementary oligonucleotides does not have any unusual optical or thermodynamic properties. By contrast, the very stable structures formed by the individual single-stranded triplet repeat sequences are thermally and thermodynamically unusual. We discuss this stable, triplet repeat, single-stranded structure and its interconversion with duplex in terms of triplet expansion diseases.

Importance of minor groove functional groups for the stability of DNA duplexes

Eight oligonucleotide duplexes have been prepared with four pairs of selected complementary pairs of native/analogue heterocyclic bases incorporated at a selected test site. The base pairs vary in the nature of their functionality in the minor groove. Each pair has a minor groove purine amino group present or absent, and correspondingly has a minor grove pyrimidine carbonyl present or absent. Loss of duplex stability is most notable when the minor groove pyrimidine carbonyl is absent although in other respects normal Watson-Crick hydrogen bonding is maintained in these sequences. These differences in stability are discussed in terms of possible variations in minor groove hydration.

Phosphate backbone neutralization increases duplex DNA flexibility: a model for protein binding

An important component of protein-DNA recognition is the charge neutralization of DNA backbone phosphates and subsequent protein-induced DNA bending. Replacement of phosphates by neutral methylphosphonates has previously been shown to be a model for protein-induced bending. In addition to bending, the neutralization process may change the inherent flexibility of the DNA--a feature never before tested. We have developed a method to measure the differential flexibility of duplex DNA when methylphosphonate substitutions are made and find that the local flexibility is increased up to 40%. These results imply that backbone-neutralization-dependent DNA flexibility augments DNA-binding motifs in protein-DNA recognition processes.

Calculating sequence-dependent melting stability of duplex DNA oligomers and multiplex sequence analysis by graphs

The analytical methods for characterizing DNA sequence-dependent thermodynamic stability have been reviewed. A set of n-n sequence stability parameters is presented. Examples in which these values are used to calculate the thermodynamic stability of short duplex DNA oligomers are presented. The problem of determining sets of isothermal sequences is addressed by representing DNA sequences as graphs. Representing DNA sequences by a graph descriptor with special mathematical properties minimizes the computational difficulty of determining the number of DNA sequences with identical predicted thermodynamic stability. This is achieved by replacement of a whole set of sequences by a single representative. Applications of this concept were demonstrated for sequences assembled from individual bases and sequences assembled from oligomeric blocks.

DNA, RNA and hybrid RNA-DNA oligomers of identical sequence: structural and dynamic differences

A 27-mer sequence was synthesised as DNA duplex (DD), RNA duplex (RR), and RNA-DNA (RD) hybrid in order to characterise their structural and dynamic features. The hydrodynamic radius (Rh) and the rise (b) values of the three samples were consistent with the conformations predicted by CD analysis. The value of the torsional constant (alpha) of the samples containing RNA was approximately twice that of the DD sample and followed the order: DD < RD < RR. The same order was observed in the thermodynamic stability and in the reduction of the electrophoretic mobility. gamma-Ray footprinting analysis was carried out to resolve the individual strand conformation in the hybrid. The RNA strand preserved its conformation, while the DNA strand showed local deformations mainly at TA and TG steps.

Predicting sequence-dependent melting stability of short duplex DNA oligomers

Many important applications of DNA sequence-dependent hybridization reactions have recently emerged. This has sparked a renewed interest in analytical calculations of sequence-dependent melting stability of duplex DNA. In particular, for many applications it is often desirable to accurately predict the transition temperature, or tm of short duplex DNA oligomers (approximately 20 base pairs or less) from their sequence and concentration. The thermodynamic analytical method underlying these predictive calculations is based on the nearest-neighbor model. At least 11 sets of nearest-neighbor sequence-dependent thermodynamic parameters for DNA have been published. These sets are compared. Use of the nearest-neighbor sets in predicting tm from the DNA sequence is demonstrated, and the ability of the nearest-neighbor parameters to provide accurate predictions of experimental tm's of short duplex DNA oligomers is assessed.

The effects of 5R-5,6-dihydro-5-hydroxythymidine on duplex DNA stability and structure

An improved method for the chemical synthesis of oligodeoxynucleotides containing 5 R -5,6-dihydro-5-hydroxythymidine (1) at a defined site is reported. UV melting studies carried out on duplexes containing1synthesized in this manner correlate with previously reported enzyme inhibition experiments, as well as computational studies. The melting experiments suggest that1destabilizes duplex DNA, but that the lesion preferentially base pairs to deoxyadenosine. These experiments also suggest that the presence of1in a duplex disrupts base pairing at the 5'-adjacent nucleotide and results in the thermally preferred misincorporation of purines opposite the 5'-deoxyadenosine stacked above 1 at this position. Despite the disruptions in base stacking, the UV melting experiments and enzymatic ligation/electrophoretic migration assays are consistent with the predicted macroscopic duplex structure containing intrahelical nucleotides.

Structural and dynamical properties of two DNA oligomers with the same base composition and different sequence

We compared the structural and dynamical properties of two DNA fragments, 27 bp long, having the same base composition but a different sequence. This work aims to understand how the base sequence on a purine rich strand in a double helix, which is important for many biological functions, is related to structural features and to measurable physical quantities. Structural characterization of the two samples was performed both by conventional spectroscopic methods (circular dichroism and UV denaturation experiments) and by means of a gamma-ray footprinting technique which gives information on fine conformational differences. Dynamical features of the samples were studied by fluorescence polarization anisotropy (FPA) measurements which allow the evaluation of some hydrodynamic parameters, such as the hydrodynamic radius and the elastic torsion constant of DNA. Using a gamma-ray footprinting technique, we observed that the interruption of the long homopurine-homopyrimidine run in the control sample, due to the 'scrambling' operation, alters the DNA three-dimensional structure, also at nucleotide level. Besides, an increase in thermal stability and in the torsional rigidity of the 'scrambled' sample was observed. A possible association between base-stacking interaction and torsional rigidity was inferred from the comparison of the two samples.

The native and the heat-induced denatured states of alpha-chymotrypsinogen A: thermodynamic and spectroscopic studies

We report the first protein phase-diagram characterized by a combination of volumetric, calorimetric, and spectroscopic techniques. More specifically, we use ultrasonic velocimetry, densimetry, and differential scanning calorimetry, in conjunction with UV absorbance and CD spectroscopy to detect and to characterize the conformational transitions of alpha-chymotrypsinogen A as a function of both pH and temperature. As judged by the CD spectra, we find that, at room temperature, the protein remains in the native state over the entire pH range investigated (pH 1 to 10). The melting profiles of the native state reveal three distinct pH domains in which protein denaturation produces different final states. Below pH 3.1, we find the heat-induced denatured state of the protein to be molten globule (MG), lacking the native-like tertiary structure, while exhibiting significant secondary structural elements. At neutral and alkaline pH, we find the heat-induced denatured state to be unfolded (U), lacking both tertiary and secondary structures, while being structurally similar to the urea-unfolded state. At intermediate pH values (between pH 3.1 and 7), we find the heat-induced denatured state to exhibit properties characteristic of both the MG and U states. Although at room temperature the protein remains native within the whole pH range studied (pH 1 to 10), our volumetric data reveal that the native state slightly "softens" at low pH, probably, due to pH-induced alterations in electrostatic forces causing the packing of the protein interior at low pH and room temperature to become less "tight". This softening of the protein at low pH is reflected in an 8% increase in the intrinsic compressibility, kM, of the protein "native" state. Our volumetric data also allow us to conclude that the heat-induced MG state retains a liquid-like, water-inaccessible core, with a volume that corresponds to about 40% of the solvent-inaccessible core of the native state. By contrast, our volumetric data are consistent with the U state of the protein being essentially unfolded, with the majority of its constituent atomic groups being solvent exposed and, therefore, strongly hydrated.

Binding of a hairpin polyamide in the minor groove of DNA: sequence-specific enthalpic discrimination

Hairpin polyamides are synthetic ligands for sequence-specific recognition in the minor groove of double-helical DNA. A thermodynamic characterization of the DNA-binding properties exhibited by a six-ring hairpin polyamide, ImPyPy-gamma-PyPyPy-beta-Dp (where Im = imidazole, Py = pyrrole, gamma = gamma-aminobutyric acid, beta = beta-alanine, and Dp = dimethylaminopropylamide), reveals an approximately 1-2 kcal/mol greater affinity for the designated match site, 5'-TGTTA-3', relative to the single base pair mismatch sites, 5'-TGGTA-3' and 5'-TATTA-3'. The enthalpy and entropy data at 20 degrees C reveal this sequence specificity to be entirely enthalpic in origin. Correlations between the thermodynamic driving forces underlying the sequence specificity exhibited by ImPyPy-gamma-PyPyPy-beta-Dp and the structural properties of the heterodimeric complex of PyPyPy and ImPyPy bound to the minor groove of DNA provide insight into the molecular forces that govern the affinity and specificity of pyrrole-imidazole polyamides.

Stereodifferentiation--the effect of P chirality of oligo(nucleoside phosphorothioates) on the activity of bacterial RNase H

P stereoregular phosphorothioate analogs of pentadecamer 5'-d(AGATGTTTGAGCTCT)-3' were synthesized by the oxathiaphospholane method. Their diastereomeric purity was assigned by means of enzymatic degradation with nuclease P1 and, independently, with snake venom phosphodiesterase. DNA-RNA hybrids formed by phosphorothioate oligonucleotides (PS-oligos) with the corresponding complementary pentadecaribonucleotide were treated with bacterial RNase H. The DNA-RNA complex containing the PS-oligo of [all-RP] configuration was found to be more susceptible to RNase H-dependent degradation of the pentadecaribonucleotide compared with hybrids containing either the [all-SP] counterpart or the so called 'random mixture of diastereomers' of the pentadeca(nucleoside phosphorothioate). This stereodependence of RNase H action was also observed for a polyribonucleotide (475 nt) hybridized with these phosphorothioate oligonucleotides. The results of melting studies of PS-oligo-RNA hybrids allowed a rationalization of the observed stereodifferentiation in terms of the higher stability of heterodimers formed between oligoribonucleotides and [all-RP]-oligo(nucleoside phosphorothioates), compared with the less stable heterodimers formed with [all-SP]-oligo(nucleoside phosphorothioates) or the random mixture of diastereomers.