Effect of loop sequence and size on DNA aptamer stability

The thrombin aptamer is a 15-mer oligodeoxyribonucleotide that folds into a unimolecular quadruplex consisting of a stack of two guanine quartets connected by two external loops and one central loop and possesses a high affinity for thrombin. We have undertaken a systematic examination, in KCl, of the thermodynamic stability of thrombin aptamer analogues containing sequence modifications in one or more of the loops, as well as in the number of quartets. UV melting studies have been carried out to obtain the relevant thermodynamic parameters for these aptamers. van't Hoff analysis of these data, with a two-state model for unimolecular denaturation, gave excellent fits to the experimental observations. Thermodynamic analysis indicates that the central loop sequence in the parent aptamer is optimal for stability. Modifications in this or other loops can effect either DeltaH degrees, DeltaS degrees, or both. Addition of a single G at the 5'-end decreases stability while addition of a G at the 3'-end increases stability. Differential scanning calorimetry experiments on the thrombin aptamer reveal that a heat capacity change, not detected by UV measurements, accompanies the unfolding of the aptamer.     

Influence of loop size on the stability of intramolecular DNA quadruplexes

We have determined the stability of intramolecular DNA quadruplexes in which the four G₃-tracts are connected by non-nucleosidic linkers containing propanediol, octanediol or hexaethylene glycol, replacing the TTA loops in the human telomeric repeat sequence. We find that these sequences all fold to form intramolecular complexes, which are stabilized by lithium < sodium < potassium. Quadruplex stability increases in the order propanediol < hexaethylene glycol < octanediol. The shallower shape of the melting profile with propanediol linkers and its lower dependency on potassium concentration suggests that this complex contains fewer stacks of G-quartets. The sequence with octanediol linkers displays a biphasic melting profile, suggesting that it can adopt more than one stable structure. All these complexes display melting temperatures above 310 K in the presence of 10 mM lithium, without added potassium, in contrast to the telomeric repeat sequence. These complexes also fold much faster than the telomeric repeat and there is little or no hysteresis between their melting and annealing profiles. In contrast, the human telomeric repeat sequence and a complex containing two hexaethylene glycol groups in each loop, are less stable and fold more slowly. The melting and annealing profiles for the latter sequence show significant differences, even when heated at 0.2°C min^–1. CD spectra for the oligonucleotides containing non-nucleosidic linkers show positive maxima at 264 nm, with negative minima ~244 nm, which are characteristic of parallel quadruplex structures. These results show that the structure and stability of intramolecular quadruplexes is profoundly influenced by the length and composition of the loops.     

Protein:DNA interactions at chromosomal loop attachment sites

We have recently identified an evolutionarily conserved class of sequences that organize chromosomal loops in the interphase nucleus, which we have termed "matrix association regions" (MARs). MARs are about 200 bp long, AT-rich, contain topoisomerase II consensus sequences and other AT-rich sequence motifs, often reside near cis-acting regulatory sequences, and their binding sites are abundant (greater than 10,000 per mammalian nucleus). Here we demonstrate that the interactions between the mouse kappa immunoglobulin gene MAR and topoisomerase II or the "nuclear matrix" occur between multiple and sometimes overlapping binding sites. Interestingly, the sites most susceptible to topoisomerase II cleavage are localized near the breakpoints of a previously described illegitimate recombination event. The presence of multiple binding sites within single MARs may allow DNA and RNA polymerase passage without disrupting primary loop organization.     

Thermodynamics of DNA binding and distortion by the hyperthermophile chromatin protein Sac7d

Sac7d is a hyperthermophile chromatin protein which binds non-specifically to the minor groove of duplex DNA and induces a sharp kink of 66 degrees with intercalation of valine and methionine side-chains. We have utilized the thermal stability of Sac7d and the lack of sequence specificity to define the thermodynamics of DNA binding over a wide temperature range. The binding affinity for poly(dGdC) was moderate at 25 degrees C (Ka = 3.5(+/-1.6) x 10(6) M(-1)) and increased by nearly an order of magnitude from 10 degrees C to 80 degrees C. The enthalpy of binding was unfavorable at 25 degrees C, and decreased linearly from 5 degrees C to 60 degrees C. A positive binding heat at 25 degrees C is attributed in part to the energy of distorting DNA, and ensures that the temperature of maximal binding affinity (75.1+/-5.6 degrees C) is near the growth temperature of Sulfolobus acidocaldarius. Truncation of the two intercalating residues to alanine led to a decreased ability to bend and unwind DNA at 25 degrees C with a small decrease in binding affinity. The energy gained from intercalation is slightly greater than the free energy penalty of bending duplex DNA. Surprisingly, reduced distortion from the double alanine substitution did not lead to a significant decrease in the heat of binding at 25 degrees C. In addition, an anomalous positive DeltaCp of binding was observed for the double alanine mutant protein which could not be explained by the change in polar and apolar accessible surface areas. Both the larger than expected binding enthalpy and the positive heat capacity can be explained by a temperature dependent structural transition in the protein-DNA complex with a Tm of 15-20 degrees C and a DeltaH of 15 kcal/mol. Data are discussed which indicate that the endothermic transition in the complex is consistent with DNA distortion.     

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. © 2013 Wiley Periodicals, Inc. Biopolymers 101: 1–12, 2014.     

The solution structure of the Sac7d/DNA complex: a small-angle X-ray scattering study.

Small-angle X-ray scattering has been used to study the structure of the multimeric complexes that form between double-stranded DNA and the archaeal chromatin protein Sac7d from Sulfolobus acidocaldarius. Scattering data from complexes of Sac7d with a defined 32-mer oligonucleotide, with poly[d(GC)], and with E. coli DNA indicate that the protein binds along the surface of an extended DNA structure. Molecular models of fully saturated Sac7d/DNA complexes were constructed using constraints from crystal structure and solution binding data. Conformational space was searched systematically by varying the parameters of the models within the constrained set to find the best fits between the X-ray scattering data and simulated scattering curves. The best fits were obtained for models composed of repeating segments of B-DNA with sharp kinks at contiguous protein binding sites. The results are consistent with extrapolation of the X-ray crystal structure of a 1:1 Sac7d/octanucleotide complex [Robinson, H., et al. (1998) Nature 392, 202-205] to polymeric DNA. The DNA conformation in our multimeric Sac7d/DNA model has the base pairs tilted by about 35 degrees and displaced 3 A from the helix axis. There is a large roll between two base pairs at the protein-induced kink site, resulting in an overall bending angle of about 70 degrees for Sac7d binding. Regularly repeating bends in the fully saturated complex result in a zigzag structure with negligible compaction of DNA. The Sac7d molecules in the model form a unique structure with two left-handed helical ribbons winding around the outside of the right-handed duplex DNA.     

Probing the DNA kink structure induced by the hyperthermophilic chromosomal protein Sac7d

Sac7d, a small, abundant, sequence-general DNA-binding protein from the hyperthermophilic archaeon Sulfolobus acidocaldarius, causes a single-step sharp kink in DNA (~60°) via the intercalation of both Val26 and Met29. These two amino acids were systematically changed in size to probe their effects on DNA kinking. Eight crystal structures of five Sac7d mutant–DNA complexes have been analyzed. The DNA-binding pattern of the V26A and M29A single mutants is similar to that of the wild-type, whereas the V26A/M29A protein binds DNA without side chain intercalation, resulting in a smaller overall bending (~50°). The M29F mutant inserts the Phe29 side chain orthogonally to the C2pG3 step without stacking with base pairs, inducing a sharp kink (~80°). In the V26F/M29F-GCGATCGC complex, Phe26 intercalates deeply into DNA bases by stacking with the G3 base, whereas Phe29 is stacked on the G15 deoxyribose, in a way similar to those used by the TATA box-binding proteins. All mutants have reduced DNA-stabilizing ability, as indicated by their lower T_m values. The DNA kink patterns caused by different combinations of hydrophobic side chains may be relevant in understanding the manner by which other minor groove-binding proteins interact with DNA.     

Indirect effects of heterospecific interactions on progeny size through maternal stress

Maternal effects are increasingly being recognized as an important pre-natal source of life history variation in the next generation. The present study uses a field experiment to explore the influence of heterospecific interactions on the reproductive output and offspring characteristics of a common Indo-Pacific damselfish, Pomacentrus amboinensis . On the Great Barrier Reef pairs of breeding P. amboinensis were placed on isolated patch reefs and to half of the pairs resource competitors (other planktivorous damselfishes), and predators of eggs and juveniles were added. Females inhabiting patches with heterospecifics had more aggressive interactions and higher levels of the stress hormone, cortisol. Neither the number of clutches nor number of eggs produced differed among treatments. The size of larvae at hatching was found to be reduced as a result of the stress associated with increased interactions with heterospecific and the transfer of cortisol to offspring. This stress-associated mechanism appears to be an important and directional source of life history variability, but the individual nature of the maternal response is likely to result in a conclusion of a diversified bet hedging reproductive strategy when viewed at the local population level. These findings highlight the complex determinants of individual success and the important role of parental well-being in the population dynamics of the next generation.     

Crystal structures of the chromosomal proteins Sso7d/Sac7d bound to DNA containing T-G mismatched base-pairs

Sso7d and Sac7d are two small chromatin proteins from the hyperthermophilic archaeabacterium Sulfolobus solfataricus and Sulfolobus acidocaldarius, respectively. The crystal structures of Sso7d-GTGATCGC, Sac7d-GTGATCGC and Sac7d-GTGATCAC have been determined and refined at 1.45 A, 2.2 A and 2.2 A, respectively, to investigate the DNA binding property of Sso7d/Sac7d in the presence of a T-G mismatch base-pair. Detailed structural analysis revealed that the intercalation site includes the T-G mismatch base-pair and Sso7d/Sac7d bind to that mismatch base-pair in a manner similar to regular DNA. In the Sso7d-GTGATCGC complex, a new inter-strand hydrogen bond between T2O4 and C14N4 is formed and well-order bridging water molecules are found. The results suggest that the less stable DNA stacking site involving a T-G mismatch may be a preferred site for protein side-chain intercalation.     

Chromosome length and DNA loop size during early embryonic development of Xenopus laevis

The looped organization of the eukaryotic genome mediated by a skeletal framework of non-histone proteins is conserved throughout the cell cycle. The radial loop/scaffold model envisages that the higher order architecture of metaphase chromosomes relies on an axial structure around which looped DNA domains are radially arranged through stable attachment sites. In this light we investigated the relationship between the looped organization and overall morphology of chromosomes. In developing Xenopus laevis embryos at gastrulation, the bulk of the loops associated with histone-depleted nuclei exhibit a significant size increase, as visualized by fluorescence microscopy of the fully extended DNA halo surrounding high salt treated, ethidium bromide stained nuclei. This implies a reduction in the number of looped domains anchored to the supporting nucleoskeletal structure. The cytological analysis of metaphase plates from acetic acid fixed whole embryos, carried out in the absence of drugs inducing chromosome condensation, reveals a progressive thickening and shortening of metaphase chromosomes during development. We interpret these findings as a strong indication that the size and number of DNA loops influence the thickness and length of the chromosomes, respectively. The quantitative analysis of chromosome length distributions at different developmental stages suggests that the shortening is timed differently in different embryonic cells.     

The acid-induced folded state of Sac7d is the native state

Sac7d unfolds at low pH in the absence of salt, with the greatest extent of unfolding obtained at pH 2. We have previously shown that the acid unfolded protein is induced to refold by decreasing the pH to 0 or by addition of salt (McCrary BS, Bedell J. Edmondson SP, Shriver JW, 1998, J Mol Biol 276:203-224). Both near-ultraviolet circular dichroism spectra and ANS fluorescence enhancements indicate that the acid- and salt-induced folded states have a native fold and are not molten globular. 1H,15N heteronuclear single quantum coherence NMR spectra confirm that the native, acid-, and salt-induced folded states are essentially identical. The most significant differences in amide 1H and 15N chemical shifts are attributed to hydrogen bonding to titrating carboxyl side chains and through-bond inductive effects. The 1H NMR chemical shifts of protons affected by ring currents in the hydrophobic core of the acid- and salt-induced folded states are identical to those observed in the native. The radius of gyration of the acid-induced folded state at pH 0 is shown to be identical to that of the native state at pH 7 by small angle X-ray scattering. We conclude that acid-induced collapse of Sac7d does not lead to a molten globule but proceeds directly to the native state. The folding of Sac7d as a function of pH and anion concentration is summarized with a phase diagram that is similar to those observed for other proteins that undergo acid-induced folding except that the A-state is encompassed by the native state. These results demonstrate that formation of a molten globule is not a general property of proteins that are refolded by acid.     

Habitat edges and predator-prey interactions: effects on critical patch size.

We use partial differential equation models to examine the effects of cross-edge incursions by a predator on the persistence or extinction of a patch-resident prey species. For each of two predator-incursion profiles (namely, a constant incursion distance and a constant loss rate for predators during incursions), we examine the conditions under which the predator can (and cannot) influence the critical patch size of a prey species.     

Effect of mutation of the Sac7d intercalating residues on the temperature dependence of DNA distortion and binding thermodynamics

Sac7d is a small chromatin protein from the hyperthermophile Sulfolobus acidocaldarius which kinks duplex DNA by approximately 66 degrees at a single base pair step with intercalation of V26 and M29 side chains. Site-directed mutagenesis coupled with calorimetric and spectroscopic data has been used to characterize the influence of the intercalating side chains on the structure and thermodynamics of the DNA complex from 5 to 85 degrees C. Two single-alanine substitutions (V26A and M29A) and five double-glycine, -alanine, -leucine, -phenylalanine, and -tryptophan substitutions of the surface residues have been created. NMR and fluorescence titrations indicated that the substitutions had little effect on the structure of the protein or DNA binding site size. Each of the mutant proteins demonstrated a temperature-dependent binding enthalpy which was correlated with a similar temperature dependence in the structure of the complex reflected by changes in fluorescence and circular dichroism. A positive heat capacity change (DeltaC(p)) for DNA binding was observed for only those mutants which also demonstrated a thermotropic structural transition in the complex, and the temperature range for the positive DeltaC(p) coincided with that observed for the structural transition. The thermodynamic data are interpreted using a model in which binding is linked to an endothermic distortion of the DNA in the complex. The results support the proposal that the unfavorable enthalpy of binding of Sac7d at 25 degrees C is due in part to the distortion of DNA.     

Partial B-to-A DNA transition upon minor groove binding of protein Sac7d monitored by Raman spectroscopy

Members of the Sso7d/Sac7d protein family and other related proteins are believed to play an important role in DNA packaging and maintenance in archeons. Sso7d/Sac7d are small, abundant, basic, and nonspecific DNA-binding proteins of the hyperthermophilic archeon Sulfolobus. Structures of several complexes of Sso7d/Sac7d with DNA octamers are known. These structures are characterized by sequence unspecific minor groove binding of the proteins and sharp kinking of the double helix. Corresponding Raman vibrational signatures have been identified in this study. A Raman spectroscopic analysis of Sac7d binding to the oligonucleotide decamer d(GAGGCGCCTC)(2) reveals large conformational perturbations in the DNA structure upon complex formation. Perturbed Raman bands are associated with the vibrational modes of the sugar phosphate backbone and frequency shifts of bands assigned to nucleoside vibrations. Large changes in the DNA backbone and partial B- to A-form DNA transitions are indicated that are closely associated with C2'-endo/anti to C3'-endo/anti conversion of the deoxyadenosyl moiety upon Sac7d binding. The major spectral feature of Sac7d binding is kinking of the DNA. Raman markers of minor groove binding do not largely contribute to spectral differences; however, clear indications for minor groove binding come from G-N2 and G-N3 signals that are supported by Trp24 features. Trp24 is the only tryptophan present in Sac7d and binds to guanine N3, as has been demonstrated clearly in X-ray structures of Sac7d-DNA complexes. No changes of the Sac7d secondary structure have been detected upon DNA binding.     

Intramolecular DNA triplexes in supercoiled plasmids. I. Effect of loop size on formation and stability

The capacity of four oligopurine.oligopyrimidine (pur.pyr) sequences with different lengths of interruptions in the center [GAA)4(N)n(GAA)4G) (n = 3, 5, 7, and 9) to adopt intramolecular DNA triplexes was evaluated in recombinant plasmids. The hyperreactive patterns of the pur.pyr inserts to specific chemical probes (OsO4, diethyl pyrocarbonate, and dimethyl sulfate) at the base pair level demonstrate that intramolecular triplexes with identical 12-base triads in the stem but with different loop sizes (4, 6, 8, and 10 bases) can form in supercoiled plasmids. Furthermore, the extent of OsO4 modification was measured as a function of temperature and of average negative supercoil density. In addition, the transition free energy of B-DNA to triplexes at pH 4.5 was determined by two-dimensional electrophoresis. These comparative studies show that longer loops require more supercoil energy for triplex formation and are less thermostable than triplexes with shorter loops. Also, it may be that not only the loop size but the base composition of the loop region affects the structural transition and triplex stability. Thus, these results significantly broaden the range of natural pur.pyr sequences that may adopt triplexes.     

Bead size effects on protein-mediated DNA looping in tethered-particle motion experiments

Tethered particle motion (TPM) has become an important tool for single-molecule studies of biomolecules; however, concerns remain that the method may alter the dynamics of the biophysical process under study. We investigate the effect of the attached microsphere on an illustrative biological example: the formation and breakdown of protein-mediated DNA loops in the lac repressor system. By comparing data from a conventional TPM experiment with 800 nm polystyrene beads and dark-field TPM using 50 nm Au nanoparticles, we found that the lifetimes of the looped and unlooped states are only weakly modified, less than two-fold, by the presence of the large bead. This is consistent with our expectation of weak excluded-volume effects and hydrodynamic surface interactions from the cover glass and microsphere.     

The role of loop 7 in mediating calcineurin regulation

Calcineurin (CN) is a heterodimer protein consisting of a 61 kDa catalytic subunit A and a 19 kDa regulatory subunit B. It plays a critical role in T-cell activation and is involved in many cellular processes. Regulation of CN is rather complex, including a number of factors such as divalent metal ions (primarily Ca(2+) and Mn(2+)), calmodulin (CaM) and autoinhibition (AI) segment. Previously, we reported that a loop 7 deletion mutant (V314) in subunit A exhibited high phosphatase activity, although the mechanism for the surprising activity enhancement and whether the activity change applies to other loop 7 residues were not known. In order to probe the role of loop 7, we have carried out extensive mutagenesis experiments, followed by systematic activity assays under a number of regulatory conditions. All mutants, including single deletion mutants Y315, N316 and double deletion mutant V314Y315, showed increased phosphatase activity. Significantly, activities of the mutants containing the V314 deletion, namely V314 and V314Y315, were no longer regulated by regulatory subunit B. These results, along with the structure analysis, suggest that loop 7 as a whole plays an important role in mediating CN's regulation through bridging the regulatory subunit and catalytic core and interaction with the AI segment of CN.     

Characterization of DNA-Hv1 histone interactions; discrimination of DNA size and shape

We have studied the formation of histone Hv1-DNA complexes using an acoustic biosensor and AFM imaging. Our results show that DNA and histone molecules aggregate into amorphous accumulations which form a compact rigid layer on the sensor’s surface. By measuring changes in the acoustic wave amplitude, it was possible to titrate surface bound DNA with Hv1 and discriminate between DNA molecules of different size and shape. From the kinetic analysis of real time data, K_eq was found equal to 3 × 10⁵ M⁻¹.