Induction of duplex to G-quadruplex transition in the c-myc promoter region by a small molecule

A major control element of the human c-myc oncogene is the nuclease-hypersensitive purine/pyrimidine-rich sequence. This double-stranded DNA fragment, corresponding to the 27-base pair segment in the nuclease-hypersensitive element of the c-myc promoter region, forms a stable Watson-Crick double helix under physiological conditions. However, this duplex DNA can be effectively converted to G-quadruplex DNA by a small molecular weight ligand. Both intermolecular and intramolecular G-quadruplex forms can be induced by this ligand. Similar transitional changes are also observed with the duplex telomeric sequence from the Oxytricha species. These results provide additional support to the idea that G-quadruplex structures may play structural roles in vivo and also provide insight into novel methodologies for rational drug design. These structurally altered DNA elements might serve as regulatory signals in gene expression or in telomere dynamics and hence are promising targets for drug action.     

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

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

Folding and stability of a primitive protein

We have previously attempted to simulate domain creation in early protein evolution by recombining polypeptide segments from non-homologous proteins, and we have described the structure of one such de novo protein, 1b11, a segment-swapped tetramer with novel architecture. Here, we have analyzed the thermodynamic stability and folding kinetics of the 1b11 tetramer and its monomeric and dimeric intermediates, and of 1b11 mutants with changes at the domain interface. Denatured 1b11 polypeptides fold into transient, folded monomers with marginal stability (DeltaG<1kcalmol(-1)) which convert rapidly ( approximately 6x10(4)M(-1)s(-1)) into dimers (DeltaG=9.8kcal/mol) and then more slowly ( approximately 3M(-1)s(-1)) into tetramers (DeltaG=28kcalmol(-1)). Segment swapping takes place during dimerization, as suggested by mass spectroscopic analysis of covalently linked peptides derived from proteolysis of a disulfide-linked dimer. Our results confirm that segment swapping and associated oligomerization are both powerful ways of stabilizing proteins, and we suggest that this may have been a feature of early protein evolution.     

Tetraplex DNA transitions within the human c-myc promoter detected by multivariate curve resolution of fluorescence resonance energy transfer

The nuclease hypersensitive element (NHE) III(I) of the c-myc promoter regulates the expression of oncogene c-myc and hence is an important anti-cancer target. Paranemic secondary structure formation within the promoter has been implicated in mechanistic regulation models. Here, it is shown that two monomeric tetraplexes form within the c-myc promoter, which coexist in solution. The development and application of a new experimental approach for detection of conformation transitions in nucleic acids [which exploits the sensitivity of fluorescence resonance energy transfer (FRET) for theoretical spectral resolution by multivariate curve resolution-alternating least-squares (MCR-ALS) method] has been used for this study. The pK(a) for tetraplex transitions are centered around 5.9 +/- 0.2 (between two intercalation topologies) and 6.8 +/- 0.1 (tetraplex to random coil). The presence of two tetraplexes has been further confirmed by S1 nuclease digestion. Finally, it is established that MCR-ALS analysis of FRET at different temperatures, pH, and salt concentrations allows resolution of pure species. Results are discussed in the light of recent observations implicating paranemic DNA motifs within the c-myc NHE in regulation of the oncogene. This method has several advantages over other methods vis-à-vis, high sensitivity and linear detection over a wide concentration range and, particularly, potential applications in intracellular probing.     

Hoechst 33258 binds to G-quadruplex in the promoter region of human c-myc

In vitro binding of Hoechst 33258 to the promoter region of human c-myc, d(GG GGAGGG TGG GGA GGG TGG GGA AGG TGG GG) which forms G-quadruplex, both in vitro and in vivo in the presence of metal ions, was investigated by equilibrium absorption, fluorescence, and kinetic surface plasmon resonance methods. Hypochromic effect in UV absorption spectra and blue shift in fluorescence emission maxima of Hoechst in the presence of quadruplex revealed that Hoechst binds to the quadruplex. Analysis of UV and fluorescence titration data revealed that Hoechst binds to quadruplex with binding affinity of the order of 10(6). Anisotropy measurements and higher lifetime obtained from time-resolved decay experiments revealed that quadruplex-bound Hoechst is rotationally restricted in a less polar environment than the bulk buffer medium. From surface plasmon resonance studies, we obtained kinetic association (k(a)) and dissociation (k(d)) of 1.23+/-0.04 x 10(5)M(-1)s(-1) and 0.686+/-0.009 s(-1), respectively. As Hoechst is known to bind A-T-rich region of duplex DNA, here we propose the likelihood of Hoechst interacting with the AAGGT loop of the quadruplex.     

The energetics of i-DNA tetraplex structures formed intermolecularly by d(TC5) and intramolecularly by d[(C5T3)3C5

Cytosine-rich DNA at low pH adopts an antiparallel tetraplex structure via the intercalation of two partially protonated, parallel stranded duplexes. This intriguing structural motif has been named i-DNA. We have used a combination of spectroscopic and calorimetric techniques to characterize the properties of an intermolecular i-DNA formed by d(TC(5)) and an intramolecular i-DNA formed by d[(C(5)T(3))(3)C(5)]. Our measurements reveal that both i-DNA complexes are enthalpically stabilized by 6.5-7.0 kcal/mol(base) and entropically destabilized by 20 cal/mol(base)/K. These values are about 50% larger than the corresponding enthalpy and entropy values per base for Watson and Crick duplexes and for Hoogsteen triplexes, while being similar to per base enthalpy and entropy values reported for G-quadruplexes. Our data also reveal a positive heat capacity change between 20 and 30 cal/mol(base)/K, values similar to that reported for polymeric Watson & Crick DNA duplexes. Solution-dependent studies reveal the overall thermal and thermodynamic stability of i-DNA complexes to be dictated by an interplay between pH and ionic strength. Based on the thermodynamic data measured, we discuss the feasibility of i-DNA formation in the context of conventional DNA sequences, while commenting on potential roles for this structural motif in biological regulatory mechanisms.     

Volume and enthalpy profiles of CO rebinding to horse heart myoglobin

Carbon monoxide binding to myoglobin was characterized using the photothermal beam deflection method. The volume and enthalpy changes coupled to CO dissociation were found to be 9.3+/-0.8 mL x mol(-1) and 7.4+/-2.8 kcal x mol(-1), respectively. The corresponding values observed for CO rebinding have the same magnitude but opposite sign: Delta V=-8.6+/-0.9 mL x mol(-1) and Delta H=-5.8+/-2.9 kcal x mol(-1). Ligand rebinding occurs as a single conformational step with a rate constant of 5 x 10(5) M(-1) s(-1) and with activation enthalpy of 7.1+/-0.8 kcal x mol(-1) and activation entropy of -22.4+/-2.8 cal x mol(-1) K(-1). Activation parameters for the ligand binding correspond to the activation parameters previously obtained using the transient absorption methods. Hence, at room temperature the CO binding to Mb can be described as a two-state model and the observed volume contraction occurs during CO-Fe bond formation. Comparing these results with CO dissociation reactions, for which two discrete intermediates were characterized, indicates differences in mechanism by which the protein modulates ligand association and dissociation.     

Structure-based thermodynamic analysis of a coupled metal binding-protein folding reaction involving a zinc finger peptide

The thermodynamics of metal binding by the prototypical Cys(2)His(2) zinc finger peptide CP-1 has been examined through the use of isothermal titration calorimetry. In cholamine buffer at pH 7.0, the binding of zinc(II) to CP-1 shows an enthalpy change of DeltaH degrees (obs) = -33.7 +/- 0.8 kcal/mol. Between one and two protons appear to be released accompanying the metal binding process. The heat of protonation of the cholamine buffer used is quite large (-11.5 kcal/mol), indicating that a portion of the observed metal binding enthalpy is due to buffer protonation. Structure-based thermodynamic analysis including the effect of water release from zinc(II) appears to account for the entropy associated with the coupled metal binding-protein folding process semiquantitatively. The strongest driving force for the reaction is the enthalpy associated with the four bonds from zinc(II) to cysteinate and histidine residues, compared with the bonds from zinc(II) to water. The binding of cobalt(II) to CP-1 is less enthalpically driven than the binding of zinc(II) by -7.6 kcal/mol. This value is approximately equal to, but slightly larger than, the expectation based on considerations of ligand field stabilization energy.     

Molecular crowding regulates the structural switch of the DNA G-quadruplex

Almost all biochemical reactions in vitro have been investigated through numerous experiments conducted in dilute solutions containing low concentrations of solutes. However, biomacromolecules such as nucleic acids, proteins, and polysaccharides are designed to function and/or form their native structures in a living cell containing high concentrations of biomacromolecules, substrates, cofactors, salts, and so on. In the present study, we have demonstrated quantitatively the effect of molecular crowding on structures and stabilities of the G-quadruplex of d(G(4)T(4)G(4)). Molecular crowding with poly(ethylene glycol) (PEG) induced a structural transition from the antiparallel to the parallel G-quadruplex of d(G(4)T(4)G(4)), while molecular crowding with polycations did not alter the structure of the antiparallel G-quadruplex. The binding constants of putrescine, one of the polycations, for d(G(4)T(4)G(4)) in the absence and presence of Na(+) are calculated to be 277 and 2.5 M(-)(1), respectively. This indicates that the polycations coordinate to d(G(4)T(4)G(4)) with electrostatic interactions. The thermodynamic parameters of the antiparallel G-quadruplex formation under the crowding and noncrowding conditions induced by putrescine were also estimated. The stability of the antiparallel G-quadruplex decreased (-DeltaG degrees (25) decreased from 28 to 22 kcal mol(-)(1)) with molecular crowding by putrescine. Also, enthalpy and entropy changes in the structural formation under crowding and noncrowding conditions clearly showed that destabilization was entropy-driven. These quantitative parameters indicated that both the volume excluded by PEG and chemical interactions such as electrostatic interaction with solute polycations are critical for determining how molecular crowding affects the structure and stability of highly ordered DNA structures.     

Rapid cooperative two-state folding of a miniature alpha-beta protein and design of a thermostable variant

There is a great deal of interest in developing small stably folded miniature proteins. A limited number of these molecules have been described, however they typically have not been characterized in depth. In particular, almost no detailed studies of the thermodynamics and folding kinetics of these proteins have been reported. Here we describe detailed studies of the thermodynamics and kinetics of folding of a 39 residue mixed alpha-beta protein (NTL9(1-39)) derived from the N-terminal domain of the ribosomal protein L9. The protein folds cooperatively and rapidly in a two-state fashion to a native state typical of those found for normal globular proteins. At pH 5.4 in 20mM sodium acetate, 100mM NaCl the temperature of maximum stability is 6 degrees C, the t(m) is 65.3 degrees C, deltaH degrees (t(m)) is between 24.6 kcalmol(-1) and 26.3 kcalmol(-1), and deltaC(p) degrees is 0.38 kcalmol(-1)deg(-1). The thermodynamic parameters are in the range expected on the basis of per residue values determined from databases of globular proteins. H/2H exchange measurements reveal a set of amides that exchange via global unfolding, exactly as expected for a normal cooperatively folded globular protein. Kinetic measurements show that folding is two-state folding. The folding rate is 640 s(-1) and the value of deltaG degrees calculated from the folding and unfolding rates is in excellent agreement with the equilibrium value. A designed thermostable variant, generated by mutating K12 to M, was characterized and found to have a t(m) of 82 degrees C. Equilibrium and kinetic measurements demonstrate that its folding is cooperative and two-state.     

Thermodynamics and folding pathway of tetraloop receptor-mediated RNA helical packing

Little is known about the thermodynamic forces that drive the folding pathways of higher-order RNA structure. In this study, we employ calorimetric [isothermal titration calorimetry (ITC) and differential scanning calorimetry (DSC)] and spectroscopic (NMR and UV) methods to characterize the thermodynamics of the GAAA tetraloop-receptor interaction, utilizing a previously described bivalent construct. ITC studies indicate that the bivalent interaction is enthalpy driven and highly stable, with a binding constant (K(obs)) of 5.5x10(6) M(-1) and enthalpy (DeltaH(obs)(o)) of -33.8 kcal/mol at 45 degrees C in 20 mM KCl and 2 mM MgCl(2). Thus, we derive the DeltaH(obs)(o) for a single tetraloop-receptor interaction to be -16.9 kcal/mol at these conditions. UV absorbance data indicate that an increase in base stacking quality contributes to the enthalpy of complex formation. These highly favorable thermodynamics are consistent with the known critical role for the tetraloop-receptor motif in the folding of large RNAs. Additionally, a significant heat capacity change (DeltaC(p,obs)(o)) of -0.24 kcal mol(-1) K(-1) was determined by ITC. DSC and UV-monitored thermal denaturation experiments indicate that the bivalent tetraloop-receptor construct follows a minimally five-state unfolding pathway and suggest the observed DeltaC(p,obs)(o) for the interaction results from a temperature-dependent unbound receptor RNA structure.     

Duplex-tetraplex equilibrium between a hairpin and two interacting hairpins of d(A-G)10 at neutral pH.

d(A-G)10 forms two helical structures at neutrality, at low ionic strength a single-hairpin duplex, and at higher ionic strength a double-hairpin tetraplex. An ionic strength-dependent equilibrium between these forms is indicated by native PAGE, which also reveals additional single-stranded species below 0.3 M Na+, probably corresponding to partially denatured states. The equilibrium also depends upon oligomer concentration: at very low concentrations, d(A-G)10 migrates faster than the random coil d(C-T)10, probably because it is a more compact single hairpin; at high concentrations, it co-migrates with the linear duplex d(A-G)10 x d(C-T)10, probably because it is a two-hairpin tetraplex. Molecular weights measured by equilibrium sedimentation in 0.1 M Na+, pH 7, reveal a mixture of monomer and dimer species at 1 degree C, but only a monomer at 40 degrees C; in 0.6 M Na+, pH 7, only a dimer species is observed at 4 degrees C. That the single- and double-stranded species are hairpin helices, is indicated by preferential S1 nuclease cleavage at the center of the oligomer(s), i.e., the loop of the hairpin(s). The UV melting transition below 0.3 M Na+ or K+, exhibits a dTm/dlog[Na+/K+] of 33 or 36 degrees C, respectively, consistent with conversion of a two-hairpin tetraplex to a single-hairpin duplex with extrahelical residues. When [Na+/K+] > or = 0.3 M, dTm/dlog [Na+/K+] is 19 or 17 degrees C, respectively, consistent with conversion of a two-hairpin tetraplex directly to single strands. A two-hairpin structure stabilized by G-tetrads is indicated by differential scanning calorimetry in 0.15 M Na+/5 mM Mg2+, with deltaH of formation per mole of the two-hairpin tetraplex of -116.9 kcal or -29.2 kcal/mol of G-tetrad.     

Equilibrium structure and folding of a helix-forming peptide: circular dichroism measurements and replica-exchange molecular dynamics simulations

We have performed experimental measurements and computer simulations of the equilibrium structure and folding of a 21-residue alpha-helical heteropeptide. Far ultraviolet circular dichroism spectroscopy is used to identify the presence of helical structure and to measure the thermal unfolding curve. The observed melting temperature is 296 K, with a folding enthalpy of -11.6 kcal/mol and entropy of -39.6 cal/(mol K). Our simulations involve 45 ns of replica-exchange molecular dynamics of the peptide, using eight replicas at temperatures between 280 and 450 K, and the program CHARMM with a continuum solvent model. In a 30-ns simulation started from a helical structure, conformational equilibrium at all temperatures was reached after 15 ns. This simulation was used to calculate the peptide melting curve, predicting a folding transition with a melting temperature in the 330-350 K range, enthalpy change of -10 kcal/mol, and entropy change of -30 cal/(mol K). The simulation results were also used to analyze the peptide structural fluctuations and the free-energy surface of helix unfolding. In a separate 15-ns replica-exchange molecular dynamics simulation started from the extended structure, the helical conformation was first attained after approximately 2.8 ns, and equilibrium was reached after 10 ns of simulation. These results showed a sequential folding process with a systematic increase in the number of hydrogen bonds until the helical state is reached, and confirmed that the alpha-helical state is the global free-energy minimum for the peptide at low temperatures.     

Increase in the stability and helical content of estrogen receptor alpha in the presence of the estrogen response element: analysis by circular dichroism spectroscopy

Ligand-dependent stabilization of the estrogen receptor (ER) is often postulated, with limited support from experimental data. We studied the thermal unfolding of recombinant ERalpha by circular dichroism (CD) spectroscopy. The T(M) of unfolding of ERalpha was 38 +/- 2.4 degrees C, and the van't Hoff enthalpy of unfolding was 31.7 +/- 3.4 kcal/mol in the absence of ligands. Addition of estradiol (E(2)) increased the T(M) to 43.6 +/- 2.3 degrees C, while addition of E(2) and an oligonucleotide harboring the estrogen response element (ERE) increased the T(M) to 47.9 +/- 1.6 degrees C. Addition of the antiestrogen 4-hydroxytamoxifen (HT) alone did not increase the T(M); however, a combination of HT and the ERE increased the T(M) to 48.9 +/- 1.0 degrees C. The ERE alone increased the T(M) to 46.1 +/- 0.9 degrees C. Addition of E(2) alone had no effect on the apparent enthalpy of unfolding; however, the ERE alone increased the apparent enthalpy from 31.7 to 36.1 kcal/mol. ERalpha samples containing the ERE also exhibited an increase in the negative ellipticity at 208 and 222 nm, relative to that of ligand-free ERalpha, suggesting a stabilization of the alpha-helix. CD data analysis further showed that the presence of the ERE caused a large increase in alpha-helical content of ERalpha in both the presence and absence of the ligands. This increase in alpha-helical content of ERalpha was not observed in the presence of a nonspecific oligonucleotide. These results show that the ERE can increase the thermal stability of ERalpha, enhance its alpha-helical content, and facilitate the cooperativity of the folding transition.