Energetic basis of molecular recognition in a DNA aptamer

The thermal stability and ligand binding properties of the L-argininamide-binding DNA aptamer (5'-GATCGAAACGTAGCGCCTTCGATC-3') were studied by spectroscopic and calorimetric methods. Differential calorimetric studies showed that the uncomplexed aptamer melted in a two-state reaction with a melting temperature T(m)=50.2+/-0.2 degrees C and a folding enthalpy DeltaH(0)(fold)=-49.0+/-2.1 kcal mol(-1). These values agree with values of T(m)=49.6 degrees C and DeltaH(0)(fold)=-51.2 kcal mol(-1) predicted for a simple hairpin structure. Melting of the uncomplexed aptamer was dependent upon salt concentration, but independent of strand concentration. The T(m) of aptamer melting was found to increase as L-argininamide concentrations increased. Analysis of circular dichroism titration data using a single-site binding model resulted in the determination of a binding free energy DeltaG(0)(bind)=-5.1 kcal mol(-1). Isothermal titration calorimetry studies revealed an exothermic binding reaction with DeltaH(0)(bind)=-8.7 kcal mol(-1). Combination of enthalpy and free energy produce an unfavorable entropy of -TDeltaS(0)=+3.6 kcal mol(-1). A molar heat capacity change of -116 cal mol(-1) K(-1) was determined from calorimetric measurements at four temperatures over the range of 15-40 degrees C. Molecular dynamics simulations were used to explore the structures of the unligated and ligated aptamer structures. From the calculated changes in solvent accessible surface areas of these structures a molar heat capacity change of -125 cal mol(-1) K(-1) was calculated, a value in excellent agreement with the experimental value. The thermodynamic signature, along with the coupled CD spectral changes, suggest that the binding of L-argininamide to its DNA aptamer is an induced-fit process in which the binding of the ligand is thermodynamically coupled to a conformational ordering of the nucleic acid.     

Calorimetric unfolding of the bimolecular and i-motif complexes of the human telomere complementary strand, d(C(3)TA(2))(4)

A combination of spectroscopic and calorimetric techniques is used to determine the unfolding thermodynamics of the complexes formed by the complementary sequence of the human telomere, d(C(3)TA(2))(4), in the pH range of 4.2 to 6. Calorimetric melting curves show biphasic transitions; both transitions are shifted to higher temperatures as the pH is decreased, indicative of cytosine protonation, which favors the formation of C*C(+) base pairs. Furthermore, the transition temperature, T(M), of the lower transition depends on strand concentration, while the T(M) of the higher transition is independent of strand concentration, indicating the following sequential melting: bimolecular complex(s)-->intramolecular complex-->random coil. The thermodynamic profiles for the formation of each complex, bimolecular and i-motif reveals small favorable free energy terms resulting from favorable enthalpy-unfavorable entropy compensations, uptake of protons, marginal uptake of counterions (i-motif) and marginal release of water molecules (i-motif). Furthermore, an enthalpy of 3.2 kcal/mol (bimolecular complex) and 5.0 kcal/mol (i-motif) is estimated for a single C*C(+)/C*C(+) base-pair stack.     

Effects of sodium ions on DNA duplex oligomers: improved predictions of melting temperatures

Melting temperatures, T(m), were systematically studied for a set of 92 DNA duplex oligomers in a variety of sodium ion concentrations ranging from 69 mM to 1.02 M. The relationship between T(m) and ln [Na(+)] was nonlinear over this range of sodium ion concentrations, and the observed melting temperatures were poorly predicted by existing algorithms. A new empirical relationship was derived from UV melting data that employs a quadratic function, which better models the melting temperatures of DNA duplex oligomers as sodium ion concentration is varied. Statistical analysis shows that this improved salt correction is significantly more accurate than previously suggested algorithms and predicts salt-corrected melting temperatures with an average error of only 1.6 degrees C when tested against an independent validation set of T(m) measurements obtained from the literature. Differential scanning calorimetry studies demonstrate that this T(m) salt correction is insensitive to DNA concentration. The T(m) salt correction function was found to be sequence-dependent and varied with the fraction of G.C base pairs, in agreement with previous studies of genomic and polymeric DNAs. The salt correction function is independent of oligomer length, suggesting that end-fraying and other end effects have little influence on the amount of sodium counterions released during duplex melting. The results are discussed in the context of counterion condensation theory.     

Polymer principles of protein calorimetric two-state cooperativity

The experimental calorimetric two-state criterion requires the van't Hoff enthalpy DeltaH(vH) around the folding/unfolding transition midpoint to be equal or very close to the calorimetric enthalpy DeltaH(cal) of the entire transition. We use an analytical model with experimental parameters from chymotrypsin inhibitor 2 to elucidate the relationship among several different van't Hoff enthalpies used in calorimetric analyses. Under reasonable assumptions, the implications of these DeltaH(vH)'s being approximately equal to DeltaH(cal) are equivalent: Enthalpic variations among denatured conformations in real proteins are much narrower than some previous lattice-model estimates, suggesting that the energy landscape theory "folding to glass transition temperature ratio" T(f) /T(g) may exceed 6.0 for real calorimetrically two-state proteins. Several popular three-dimensional lattice protein models, with different numbers of residue types in their alphabets, are found to fall short of the high experimental standard for being calorimetrically two-state. Some models postulate a multiple-conformation native state with substantial pre-denaturational energetic fluctuations well below the unfolding transition temperature, or predict a significant post-denaturational continuous conformational expansion of the denatured ensemble at temperatures well above the transition point, or both. These scenarios either disagree with experiments on protein size and dynamics, or are inconsistent with conventional interpretation of calorimetric data. However, when empirical linear baseline subtractions are employed, the resulting DeltaH(vH)/DeltaH(cal)'s for some models can be increased to values closer to unity, and baseline subtractions are found to correspond roughly to an operational definition of native-state conformational diversity. These results necessitate a re-assessment of theoretical models and experimental interpretations.     

Preliminary dielectric measurement and analysis protocol for determining the melting temperature and binding energy of short sequences of DNA in solution

Measurement of the real dielectric constant of bulk buffer solutions containing short sequences of DNA as a function of temperature through the DNA melting or denaturiztion transition can be used to determine melting temperatures, T(m), and to estimate the binding energy of the complimentary strands. We describe a preliminary dielectric measurement and analysis protocol to determine these parameters and its application to two known short sequences. The relative real dielectric constant for the bulk solutions was determined over the frequency range of 50 Hz-20 kHz and temperature range of <40-65 degrees C. The measurements were performed on dilute solutions and utilized low electric field strengths. Based on fits to the data by modified sigmoid functions, the melting temperatures, width of transition, and binding energy for the two sequences in solution were estimated. It was observed that the order of the transition appeared to be second order. The results were then compared against predictions of a number of models from the literature that provide theoretical estimates for the melting temperatures of known short sequences of DNA.     

Melting behavior of a covalently closed, single-stranded, circular DNA

We synthesized the 26-residue deoxynucleotide sequence d(TTCCT5GGAATTCCT5GGAA) which folds intramolecularly to form a dumbbell-shaped, double-hairpin structure with a gap between the 3' and the 5' ends. We used T4 polynucleotide kinase to phosphorylate the 5' end followed by T4 DNA ligase to close the 3' and 5' ends. Melting of the dumbbell structure formed by this ligated sequence produces a covalently closed, single-stranded, circular final state. We employed calorimetric and spectroscopic techniques to characterize thermodynamically the melting behavior of the ligated molecule and compared it with the corresponding melting behavior of its unligated precursor. This comparison allowed us to characterize uniquely the influence of single-stranded ring closure on intramolecular duplex melting. The data reveal that ring closure produces a thermally more stable structure which exhibits significantly altered melting thermodynamics. We rationalize these thermodynamic differences in terms of differential solvation and differential counterion association between the ligated and unligated molecules. We also note the importance of such constrained dumbbell structures as models for hairpins, cruciforms, and locally melted domains within naturally occurring DNA polymers.     

Study of the influence of metal ions on tRNA(Phe) thermal unfolding equilibria by UV spectroscopy and multivariate curve resolution

The influence of metal ions (Na(+), Mg(2+) and Cd(2+)) on the thermal unfolding of phenylalanine transfer ribonucleic acid (tRNA(Phe)) was studied by UV spectroscopy-monitored melting experiments. Absorbance data were obtained during the unfolding process in the range 220-340 nm and later analyzed by a multivariate curve resolution approach (MCR-ALS) based on factor analysis. This procedure determines the number of spectroscopically distinct conformations present during the unfolding process and reveals their concentration profiles and pure spectra, without any initial assumption having to be made about the number of steps in the unfolding pathway. From the concentration profiles and pure spectra, information such as T(m) values can be recovered. The results were compared with those obtained previously in spectroscopic and calorimetric unfolding experiments, showing that the multivariate approach recovers information that complements that obtained in traditional spectroscopic melting experiments.     

A FRET-based screening assay for nucleic acid ligands

Some of the most serious diseases are characterized by the presence of a specific secondary structure within DNA or RNA, often in the promoter or the coding region of the responsible gene, that enhances or disrupts expression of the protein. Structural elements that impact cellular function may also be formed in other genomic regions such as telomeres. Compounds that interact with such structural elements may be useful in diagnosis or treatment of patients. In this report, we present a FRET melting assay that allows testing of libraries of compounds against four different nucleic acid structures. Compounds are tested to determine whether they stabilize preformed secondary structures (i.e., whether they cause an increase in melting temperature (T(m))). This property is described by the ΔT(m) parameter, which is the difference between the T(m) of the compound-stabilized structure and the T(m) of the unbound structure. Model oligonucleotides are labeled with FAM as a fluorescent donor and TAMRA as an acceptor. The intensity of FAM fluorescence is recorded as a function of temperature. Melting temperatures are determined by the FRET method in 96-well plates; this assay could easily be converted into 384-well format.     

Conformational changes in single-strand DNA as a function of temperature by SANS

Small-angle neutron scattering (SANS) measurements were performed on a solution of single-strand DNA, 5'-ATGCTGATGC-3', in sodium phosphate buffer solution at 10 degrees C temperature increments from 25 degrees C to 80 degrees C. Cylindrical, helical, and random coil shape models were fitted to the SANS measurements at each temperature. All the shapes exhibited an expansion in the diameter direction causing a slightly shortened pitch from 25 degrees C to 43 degrees C, an expansion in the pitch direction with a slight decrease in the diameter from 43 degrees C to 53 degrees C, and finally a dramatic increase in the pitch and diameter from 53 degrees C to 80 degrees C. Differential scanning calorimeter scans of the sequence in solution exhibited a reversible two-state transition profile with a transition temperature of 47.5 +/- 0.5 degrees C, the midpoint of the conformational changes observed in the SANS measurements, and a calorimetric transition enthalpy of 60 +/- 3 kJ mol(-1) that indicates a broad transition as is observed in the SANS measurements. A transition temperature of 47 +/- 1 degrees C was also obtained from ultraviolet optical density measurements of strand melting scans of the single-strand DNA. This transition corresponds to unstacking of the bases of the sequence and is responsible for the thermodynamic discrepancy between its binding stability to its complementary sequence determined directly at ambient temperatures and determined from extrapolated values of the melting of the duplex at high temperature.     

Melting, glass transition, and apparent heat capacity of α-D-glucose by thermal analysis

The thermal behaviors of α-D-glucose in the melting and glass transition regions were examined utilizing the calorimetric methods of standard differential scanning calorimetry (DSC), standard temperature-modulated differential scanning calorimetry (TMDSC), quasi-isothermal temperature-modulated differential scanning calorimetry (quasi-TMDSC), and thermogravimetric analysis (TGA). The quantitative thermal analyses of experimental data of crystalline and amorphous α-D-glucose were performed based on heat capacities. The total, apparent and reversingheat capacities, and phase transitions were evaluated on heating and cooling. The melting temperature (T(m)) of a crystalline carbohydrate such as α-D-glucose, shows a heating rate dependence, with the melting peak shifted to lower temperature for a lower heating rate, and with superheating of around 25K. The superheating of crystalline α-D-glucose is observed as shifting the melting peak for higher heating rates, above the equilibrium melting temperature due to of the slow melting process. The equilibrium melting temperature and heat of fusion of crystalline α-D-glucose were estimated. Changes of reversing heat capacity evaluated by TMDSC at glass transition (T(g)) of amorphous and melting process at T(m) of fully crystalline α-D-glucose are similar. In both, the amorphous and crystalline phases, the same origin of heat capacity changes, in the T(g) and T(m) area, are attributable to molecular rotational motion. Degradation occurs simultaneously with the melting process of the crystalline phase. The stability of crystalline α-D-glucose was examined by TGA and TMDSC in the melting region, with the degradation shown to be resulting from changes of mass with temperature and time. The experimental heat capacities of fully crystalline and amorphous α-D-glucose were analyzed in reference to the solid, vibrational, and liquid heat capacities, which were approximated based on the ATHAS scheme and Data Bank.     

Tuning lambda6-85 towards downhill folding at its melting temperature

The five-helix bundle lambda6-85* is a fast two-state folder. Several stabilized mutants have been reported to fold kinetically near-downhill or downhill. These mutants undergo a transition to two-state folding kinetics when heated. It has been suggested that this transition is caused by increased hydrophobicity at higher temperature. Here we investigate two histidine-containing mutants of lambda6-85* to see if a weaker hydrophobic core can extend the temperature range of downhill folding. The very stable lambdaHA is the fastest-folding lambda repressor to date (k(f)(-1) approximately k(obs)(-1)=2.3 micros at 44 degrees C). It folds downhill at low temperature, but transits back to two-state folding at its unfolding midpoint. lambdaHG has a weakened hydrophobic core. It is less stable than some slower folding mutants of lambda6-85*, and it has more exposed hydrophobic surface area in the folded state. This mutant nonetheless folds very rapidly, and has the non-exponential folding kinetics of an incipient downhill folder even at the unfolding midpoint (k(m)(-1) approximately 2 micros, k(a)(-1)=15 micros at 56 degrees C). We also compare the thermodynamic melting transition of lambdaHG with the nominal two-state folding mutant lambdaQG, which has a similar melting temperature. Unlike lambdaQG, lambdaHG yields fluorescence wavelength-dependent cooperativities and probe-dependent melting temperatures. This result combined with previous work shows that the energy landscapes of lambda repressor mutants support all standard folding mechanisms.     

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

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

Melting of individual lipid components in binary lipid mixtures studied by FTIR spectroscopy, DSC and Monte Carlo simulations

Monte Carlo (MC) simulations, Differential Scanning Calorimetry (DSC) and Fourier Transform InfraRed (FTIR) spectroscopy were used to study the melting behavior of individual lipid components in two-component membranes made of DMPC and DSPC. We employed Monte Carlo simulations based on parameters obtained from DSC profiles to simulate the melting of the different lipids as a function of temperature. The simulations show good agreement with the FTIR data recorded for deuterated and non-deuterated lipids, which demonstrates that the information on the differential melting of the individual components is already contained in the calorimetric profiles. In mixtures, both lipids melt over a wide temperature range. As expected, the lipid melting events of the lipid with the lower melting temperature occur on average at lower temperatures. The simulations also yield information on the lateral distribution of the lipids that is neither directly contained in the DSC nor in the FTIR data. In the phase coexistence region, liquid disordered domains are typically richer in the lower-melting-temperature lipid species.     

Estimation of the diversity between DNA calorimetric profiles,differential melting curves and corresponding melting temperatures

The Poland–Fixman–Freire formalism was adapted for modeling of calorimetric DNA melting profiles, and applied to plasmid pBR 322 and long random sequences. We studied the influence of the difference (H_GC?H_AT) between the helix‐coil transition enthalpies of AT and GC base pairs on the calorimetric melting profile and on normalized calorimetric melting profile. A strong alteration of DNA calorimetrical profile with H_GC?H_AT was demonstrated. In contrast, there is a relatively slight change in the normalized profiles and in corresponding ordinary (optical) normalized differential melting curves (DMCs). For fixed H_GC?H_AT, the average relative deviation (S) between DMC and normalized calorimetric profile, and the difference between their melting temperatures (T_cal?T_m) are weakly dependent on peculiarities of the multipeak fine structure of DMCs. At the same time, both the deviation S and difference (T_cal?T_m) enlarge with the temperature melting range of the helix‐coil transition. It is shown that the local deviation between DMC and normalized calorimetric profile increases in regions of narrow peaks distant from the melting temperature.     

Temperature and guanidine hydrochloride dependence of the structural stability of ribonuclease T1.

The thermal unfolding of ribonuclease T1 has been studied by high-sensitivity differential scanning calorimetry as a function of temperature, [GuHCl], and scanning rate. The destabilizing effect of GuHCl has revealed that the kinetics of the unfolding transition become extremely slow as the transition temperature decreases. At pH 5.3 and zero GuHCl, the unfolding transition is centered at 59.1 degrees C; upon increasing the GuHCl concentration, the transition occurs at lower temperatures and exhibits progressively slower kinetics; so, for example, at 3 M GuHCl, the transition temperature is 40.6 degrees C and is characterized by a time constant close to 10 min. Under all conditions studied (pH 5.3, pH 7.0, [GuHCl] < 3 M), the transition is thermodynamically reversible. The slow kinetics of the transition induce significant distortions in the shape of the transition profiles that can be mistakenly interpreted as deviations from a two-state mechanism. Determination of the thermodynamic parameters from the calorimetric data has required the development of an analytical formalism that explicitly includes the thermodynamics as well as the kinetics of the transition. Using this formalism, it is shown that a two-state slow-kinetics model is capable of accurately describing the structural stability of ribonuclease T1 as a function of temperature, GuHCl concentration, and scanning rate. Multidimensional analysis of the calorimetric data has been used to estimate the intrinsic thermodynamic parameters for protein stability, the interaction parameters with GuHCl, and the time constant for the unfolding transition and its temperature dependence.     

Melting of local ordered structures in yeast 5S ribosomal RNA in aqueous salts.

Equilibrium and kinetics of thermal melting of yeast 5S ribosomal RNA in aqueous NaCl with or without Mg2+ were investigated by differential thermal melting and temperature jump methods. Two peaks (1 and 2) and a shoulder were observed in each of melting curves at ionic strength I=0.002-0.5 and linearity between each of melting temperatures T1m and T2m and log I was found at I=0.01-0.5 in the Mg2+-free solution. The local structures were found to be stabilized considerably by Mg2+. The temperature jump measurements gave the kinetic melting curve of the structure 1 at I=0.03 without Mg2+ or with 0.5 mM Mg2+. The kinetic Tm coincided well with the corresponding static Tm. For the structure 1, various parameters were calculated from the kinetic data, which indicated a double helical character of the structure 1. In terms of the values of Tm, G-C content, and enthalpy change of the transition of the structure 1 or 2, appropriateness of each of the secondary structure models of eukaryotic 5S RNA proposed previously was discussed.     

Entropy and heat capacity of DNA melting from temperature dependence of single molecule stretching

When a single molecule of double-stranded DNA is stretched beyond its B-form contour length, the measured force shows a highly cooperative overstretching transition. We have measured the force at which this transition occurs as a function of temperature. To do this, single molecules of DNA were captured between two polystyrene beads in an optical tweezers apparatus. As the temperature of the solution surrounding a captured molecule was increased from 11 degrees C to 52 degrees C in 500 mM NaCl, the overstretching transition force decreased from 69 pN to 50 pN. This reduction is attributed to a decrease in the stability of the DNA double helix with increasing temperature. These results quantitatively agree with a model that asserts that DNA melting occurs during the overstretching transition. With this model, the data may be analyzed to obtain the change in the melting entropy DeltaS of DNA with temperature. The observed nonlinear temperature dependence of DeltaS is a result of the positive change in heat capacity of DNA upon melting, which we determine from our stretching measurements to be DeltaC(p) = 60 +/- 10 cal/mol K bp, in agreement with calorimetric measurements.     

Thermodynamics of (dG--dC)3 double-helix formation in water and deuterium oxide

The thermodynamics of double-helix formation by (dG-dC)3 in H2O and D2O were measured spectrophotometrically and calorimetrically. Spectrophotometric values based on plots of inverse melting temperature vs. log of concentration are in good agreement with calorimetry. However, spectrophotometric values based on the shape of absorbance vs. temperature curves are in good agreement with calorimetry only if linear pre- and posttransition base lines are subtracted before analysis. The averages of the enthalpies measured by the three most reliable spectroscopic methods are -56.9 and -62.7 kcal/mol of helix in H2O and D2O, respectively. The corresponding calorimetric values are -59.6 and -65.8 kcal/mol of helix.