Competition for hydrogen-bond formation in the helix-coil transition and protein folding.

The problem of the helix-coil transition of biopolymers in explicit solvents, such as water, with the ability for hydrogen bonding with a solvent is addressed analytically using a suitably modified version of the Generalized Model of Polypeptide Chains. Besides the regular helix-coil transition, an additional coil-helix or reentrant transition is also found at lower temperatures. The reentrant transition arises due to competition between polymer-polymer and polymer-water hydrogen bonds. The balance between the two types of hydrogen bonding can be shifted to either direction through changes not only in temperature, but also by pressure, mechanical force, osmotic stress, or other external influences. Both polypeptides and polynucleotides are considered within a unified formalism. Our approach provides an explanation of the experimental difficulty of observing the reentrant transition with pressure and underscores the advantage of pulling experiments for studies of DNA. Results are discussed and compared with those reported in a number of recent publications with which a significant level of agreement is obtained.

Microscopic formulation of the Zimm-Bragg model for the helix-coil transition.

A microscopic spin model is proposed for the phenomenological Zimm-Bragg model for the helix-coil transition in biopolymers. This model is shown to provide the same thermophysical properties of the original Zimm-Bragg model and it allows a very convenient framework to compute statistical quantities. Physical origins of this spin model are made transparent by an exact mapping into a one-dimensional Ising model with an external field. However, the dependence on temperature of the reduced external field turns out to differ from the standard one-dimensional Ising model and hence it gives rise to different thermophysical properties, despite the exact mapping connecting them. We discuss how this point has been frequently overlooked in the recent literature.

Intersegment interactions and helix-coil transition within the generalized model of polypeptide chains approach.

The generalized model of polypeptide chains is extended to describe the helix-coil transition in a system comprised of two chains interacting side-by-side. The Hamiltonian of the model takes into account four possible types of interactions between repeated units of the two chains, i.e., helix-helix, helix-coil, coil-helix, and coil-coil. Analysis reveals when the energy I(hh)+I(cc) of (h-h, c-c) interactions overwhelms the energy I(hc)+I(ch) of mixed (h-c, c-h) interactions, the correlation length rises substantially, resulting in narrowing of the transition interval. In the opposite case, when I(hh)+I(cc)

Two scale generalized model of polypeptide chains.

The generalized model of polypeptide chains (GMPC) is expanded to simultaneously consider two types of interactions occurring over different scales. This new two scale GMPC is applied in several specific cases to examine: The combined influence of stacking or antistacking and hydrogen bonding, or spatial restrictions on the length of helical segments, on the cooperativity and temperature interval of the helix-coil transition of duplex DNA. For the cases of stacking or antistacking in combination with hydrogen bonding the model reduces to the basic uniscale model with a redefined scaling parameter Delta. Antistacking increases the cooperativity, while stacking decreases it. In each case, explanations are given in terms of different lengths of helical segments. Restrictions on the length of helical regions result in the appearance of antiferromagnetic-type correlations where there is no apparent link between cooperativity and transition interval.

Stacking heterogeneity: a model for the sequence dependent melting cooperativity of duplex DNA.

A microscopic Potts-like one-dimensional model with many particle interactions [referred as the generalized model of polypeptide chains (GMPCs)] is developed to investigate cooperativity of DNA sequence dependent melting. For modeling sequence, regular homogeneous sequences were arranged in heterogeneous blocks of various lengths. Within the framework of the GMPC the authors show that the inclusion of stacking interaction heterogeneity relative to homogeneous hydrogen bond interactions leads to an unexpected and quite remarkable increase in melting cooperativity for small blocks. In some cases this tendency persists for long blocks having sharp sequence heterogeneity.

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.

The helix-coil transition in heterogeneous double stranded DNA: microcanonical method.

A microscopic Potts-like one-dimensional model with many-particle interactions is developed to construct a statistical mechanical description of the melting of heterogeneous sequence duplex DNA. For this model, referred as the generalized model of polypeptide chains (GMPC), a closed-form expression for the free energy is derived. The characteristic equation of the model enables estimates on the melting temperature and transition interval, consistent with results obtained from more classical approaches. From the characteristic equation of the model, the temperature-dependent statistical weight parameter for helical states is evaluated. This parameter is shown to change throughout the transition from a harmonic form in early regions of the transition to an arithmetic form in later stages. The GMPC is extended to consider the influence of sequence heterogeneity in the melting of duplex DNA.

A semiflexible polymer model applied to loop formation in DNA hairpins.

A statistical mechanical "zipper" model is applied to describe the equilibrium melting of short DNA hairpins with poly(dT) loops ranging from 4 to 12 bases in the loop. The free energy of loop formation is expressed in terms of the persistence length of the chain. This method provides a new measurement of the persistence length of single-stranded DNA, which is found to be approximately 1.4 nm for poly(dT) strands in 100 mM NaCl. The free energy of the hairpin relative to the random coil state is found to scale with the loop size with an apparent exponent of > or = 7, much larger than the exponent of approximately 1.5-1.8 expected from considerations of loop entropy alone. This result indicates a strong dependence of the excess stability of the hairpins, from stacking interactions of the bases within the loop, on the size of the loop. We interpret this excess stability as arising from favorable hydrophobic interactions among the bases in tight loops and which diminish as the loops get larger. Free energy profiles along a generalized reaction coordinate are calculated from the equilibrium zipper model. The transition state for hairpin formation is identified as an ensemble of looped conformations with one basepair closing the loop, and with a lower enthalpy than the random coil state. The equilibrium model predicts apparent activation energy of approximately -11 kcal/mol for the hairpin closing step, in remarkable agreement with the value obtained from kinetics measurements.

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.

Novel signal amplification technology with applications in DNA and protein detection systems.

A non-enzymatic approach to signal amplification has practical advantages over conventional target amplification methods. We have designed a simple, cost-efficient signal amplification system that can be used to enhance the detection of nucleic acids or protein. The signal amplification process requires initial capture of analyte by a specific probe, which, depending on the analyte, can be an oligomer or an antibody. Once the analyte is captured, amplification moieties are applied to significantly enhance the sensitivity of analyte detection. Nucleic acid amplification is typically greater than 1000-fold, increasing the sensitivity of target detection to less than 1 amol/100 microL. This amplification strategy presents a very flexible system with components that are easily altered to accommodate diverse assay requirements.

Interactions of meso-tetra-(4-N-oxyethylpyridyl) porphyrin, its 3-N analog and their metallocomplexes with duplex DNA.

Interactions of meso-tetra-(4-N-oxyethylpyridyl) porphyrin (TOEPyP(4)), its 3-N analog (TOEPyP(3)) and their Co, Cu, Ni, Zn metallocomplexes with duplex DNA have been investigated by uv/visible absorbance and circular dichrosim spectroscopies. Results reveal the interactions of these complexes with duplex DNA are of two types. (1) External binding of duplex DNA by metalloporphyrins containing Zn and Co, and (2) Binding of duplex DNA both externally and internally (by intercalation) by porphyrins not containing metals, and metalloporphyrins containing Cu and Ni. Results indicate that (4N-oxyethylpyridyl) porphyrins intercalate more preferably in the structure of duplex DNA and have weaker external binding than 3N-porphyrins.

Thermal unfolding of ribonuclease A in phosphate at neutral pH: deviations from the two-state model.

The thermal denaturation of ribonuclease A (RNase A) in the presence of phosphate at neutral pH was studied by differential scanning calorimetry (DSC) and a combination of optical spectroscopic techniques to probe the existence of intermediate states. Fourier transform infrared (FTIR) spectra of the amide I' band and far-uv circular dichroism (CD) spectra were used to monitor changes in the secondary structure. Changes in the tertiary structure were monitored by near-uv CD. Spectral bandshape changes with change in temperature were analyzed using factor analysis. The global unfolding curves obtained from DSC confirmed that structural changes occur in the molecule before the main thermal denaturation transition. The analysis of the far-uv CD and FTIR spectra showed that these lower temperature-induced modifications occur in the secondary structure. No pretransition changes in the tertiary structure (near-uv CD) were observed. The initial changes observed in far-uv CD were attributed to the fraying of the helical segments, which would explain the loss of spectral intensity with almost no modification of spectral bandshape. Separate analyses of different regions of the FTIR amide I' band indicate that, in addition to alpha-helix, part of the pretransitional change also occurs in the beta-strands.

Hybridization of single-stranded DNA targets to immobilized complementary DNA probes: comparison of hairpin versus linear capture probes.

A microtiter-based assay system is described in which DNA hairpin probes with dangling ends and single-stranded, linear DNA probes were immobilized and compared based on their ability to capture single-strand target DNA. Hairpin probes consisted of a 16 bp duplex stem, linked by a T(2)-biotin.dT-T(2) loop. The third base was a biotinylated uracil (U(B)) necessary for coupling to avidin coated microtiter wells. The capture region of the hairpin was a 3' dangling end composed of either 16 or 32 bases. Fundamental parameters of the system, such as probe density and avidin adsorption capacity of the plates were characterized. The target DNA consisted of 65 bases whose 3' end was complementary to the dangling end of the hairpin or to the linear probe sequence. The assay system was employed to measure the time dependence and thermodynamic stability of target hybridization with hairpin and linear probes. Target molecules were labeled with either a 5'-FITC, or radiolabeled with [gamma-(33)P]ATP and captured by either linear or hairpin probes affixed to the solid support. Over the range of target concentrations from 10 to 640 pmol hybridization rates increased with increasing target concentration, but varied for the different probes examined. Hairpin probes displayed higher rates of hybridization and larger equilibrium amounts of captured targets than linear probes. At 25 and 45 degrees C, rates of hybridization were better than twice as great for the hairpin compared with the linear capture probes. Hairpin-target complexes were also more thermodynamically stable. Binding free energies were evaluated from the observed equilibrium constants for complex formation. Results showed the order of stability of the probes to be: hairpins with 32 base dangling ends > hairpin probes with l6 base dangling ends > 16 base linear probes > 32 base linear probes. The physical characteristics of hairpins could offer substantial advantages as nucleic acid capture moieties in solid support based hybridization systems.

Thermodynamic, spectroscopic, and equilibrium binding studies of DNA sequence context effects in four 40 base pair deoxyoligonucleotides.

Effects of different end sequences on melting, circular dichroism spectra (CD), and enzyme binding properties were investigated for four 40 base pair, non-self-complementary duplex DNA oligomers. The center sequences of these oligoduplexes have either of two 22 base pair modules flanked on both sides by sequences differing in AT content. Temperature-induced melting transitions monitored by differential scanning calorimetry (DSC) and ultraviolet absorbance were measured for the six duplexes in buffered 115 mM Na(+) solutions. Values of the melting transition enthalpy, DeltaH(cal), and entropy, DeltaS(cal), were obtained directly from DSC experiments. Melting transition parameters, DeltaH(vH) and DeltaS(vH), were also estimated from a van't Hoff analysis of optical melting curves collected as a function of DNA concentration, assuming that the melting transition is two-state. Melting free energies (20 degrees C) evaluated from DSC melting experiments on the four duplex DNAs ranged from -52.2 to -77.5 kcal/mol. Free energies based on the van't Hoff analysis were -37.9 to -58.8 kcal/mol. Although the values are different, trends in the melting free energies of the four duplex DNAs as a function of sequence were identical in both DSC and optical analyses. Subject to several assumptions, values for the initiation free energy were estimated for each duplex, defined as DeltaG(int) = DeltaG(cal) - DeltaG(pred), where DeltaG(cal) is the experimental free energy at 20 degrees C determined from the experimentially measured values of the transition enthalpy, DeltaH(cal), and entropy, DeltaS(cal). The predicted free energy of the sequence, DeltaG(pred)(20 degrees C), is obtained using published nearest-neighbor sequence stability values. For three of the four duplexes, values of DeltaG(int) are essentially nil. In contrast, the duplex with 81.8% GC has a considerably higher estimate of DeltaG(int) = 7.1 kcal/mol. The CD spectra for the six duplexes collected over the wavelength range from 200 to 320 nm are also sequence-dependent. Factor analysis of the CD spectra by singular value decomposition revealed that the experimental CD spectra could be reconstructed from linear combinations of two minor and one major subspectra. Changes in the coefficients of the major subspectrum for different sequences reflect incremental sequence-dependent variations of the CD spectra. Equilibrium binding by BamHI restriction endonuclease to the 40 base pair DNAs whose central eight base pairs contain the recognition sequence for BamHI restriction enzyme bounded by A.T base pairs, 5'-A-GGATCC-A-3' was investigated. Binding assays were performed by titering BamHI against a constant concentration of each of the duplex DNA substrates, in the absence of Mg(2+), followed by analysis by gel retardation. Under the conditions employed, the enzyme binds but does not cleave the DNAs. Results of the assays revealed two binding modes with retarded gel mobilities. Binding isotherms for the fraction of bound DNA species versus enzyme concentration for each binding mode were constructed and analyzed with a simple two-step equilibrium binding model. This analysis provided semiquantitative estimates on the equilibrium binding constants for BamHI to the four DNAs. Values obtained for the binding constants varied only 7-fold and ranged from 6 x 10(-)(8) to 42 x 10(-)(8) M, with binding free energies from -8.6 to -9.7 (+/- 0.2) kcal/mol depending on the sequence that flanks the enzyme binding site. Unlike what was found earlier in binding studies of the 22 base pair duplexes that constitute the core modules of the present 40-mers [Riccelli, P. V., Vallone, P. M., Kashin, I., Faldasz, B. D., Lane, M. J., and Benight, A. S. (1999) Biochemistry 38, 11197-11208], no obvious relationship between binding and stability was found for these longer DNAs. Apparently, effects of flanking sequence stability on restriction enzyme binding may only be measurable in very short duplex deoxyoligonucl

Melting of cross-linked DNA IV. Methods for computer modeling of total influence on DNA melting of monofunctional adducts, intrastrand and interstrand cross-links formed by molecules of an antitumor drug.

A theoretical method is developed for calculation of melting curves of covalent complexes of DNA with antitumor drugs. The method takes into account all the types of chemical modifications of the double helix caused by platinum compounds and DNA alkylating agents: 1) monofunctional adducts bound to one nucleotide; 2) intrastrand cross-links which appear due to bidentate binding of a drug molecule to two nucleotides that are included into the same DNA strand; 3) interstrand cross-links caused by bidentate binding of a molecule to two nucleotides of different strands. The developed calculation method takes into account the following double helix alterations at sites of chemical modifications: 1) a change in stability of chemically modified base pairs and neighboring ones, that is caused by all the types of chemical modifications; 2) a change in the energy of boundaries between helical and melted regions at sites of chemical modification (local alteration of the factor of cooperativity of DNA melting), that is caused by all the types of chemical modifications, too; 3) a change in the loop entropy factor of melted regions that include interstrand cross-links; 4) the prohibition of divergence of DNA strands in completely melted DNA molecules, which is caused by interstrand cross-links only. General equations are derived, and three calculation methods are proposed to calculate DNA melting curves and the parameters that characterize the helix-coil transition.

The role of the loop in binding of an actinomycin D analog to hairpins formed by single-stranded DNA.

Our recent work has indicated that the potent antibiotic and antitumor agent actinomycin D has the ability to selectively bind and stabilize single-stranded DNA that is capable of adopting a hairpin conformation. This mechanism of DNA binding has been implicated in the drug's ability to inhibit transcription by HIV reverse transcriptase from single-stranded DNA templates. In this report, we studied the importance of the hairpin loop on the ability of the 7-amino analog of actinomycin D to selectively bind DNA hairpins. Binding dissociation constant (Kd) values were determined to be 0.22 +/- 0.11 microM for the hairpin formed from the single-stranded DNA 5'-AAAAAAATAGTTTTAAATATTTTTTT-3' (dubbed HP1). The hairpin stem without the loop resulted in binding with Kd = 2.6 +/- 0.9 microM. The drug showed low affinity for the HP1 strand fully duplexed to its complementary sequence (estimated to be at least Kd > 21 microM). Evaluation of 7-aminoactinomycin D binding to a library of thermodynamically characterized DNA hairpins revealed an affinity for the hairpin-forming sequence 5'-GGATACCCCCGTATCC-3' (dubbed ACC4) of Kd = 6.8 +/- 2.2 microM. Replacement of the terminal guanines of this sequence to generate 5'-ATATACCCCCGTATAT-3' resulted in a 10-fold increase in affinity for this hairpin compared to ACC4, to Kd = 0.74 +/- 0.06 microM. A molecular model of the ACC4actinomycin D complex reveals that significant interactions between the hairpin loop and the pentapeptide rings of the drug must occur during drug binding. Taken together, our data indicate that the composition of the stem-loop interface is critical for the selectivity of actinomycin D and its 7-amino analog for DNA hairpins and suggests that novel drugs may be designed based on selection for the desired hairpin composition.

Thermodynamic, spectroscopic, and equilibrium binding studies of DNA sequence context effects in six 22-base pair deoxyoligonucleotides.

Effects of different end sequences on stability, circular dichroism spectra (CD), and enzyme binding properties were investigated for six 22-base pair, non-self-complementary duplex DNA oligomers. The center sequences of these deoxyoligonucleotides have 8-14 base pairs in common and are flanked on both sides by sequences differing in context and A-T content. Temperature-induced melting transitions monitored by differential scanning calorimetry (DSC) and ultraviolet absorbance were measured for the six duplexes in buffered 115 mM Na(+) solutions. Values of the melting transition enthalpy, DeltaH(cal), and entropy, DeltaS(cal), were obtained directly from DSC experiments. Melting transition parameters, DeltaH(vH) and DeltaS(vH), were also estimated from van't Hoff analysis of optical melting curves collected as a function of DNA concentration, assuming a two-state melting transition. Melting free energies (20 degrees C) of the six DNAs evaluated from DSC experiments ranged from -18.7 to -32.7 kcal/mol. van't Hoff estimates of the free energies ranged from -18.5 to -48.0 kcal/mol. With either method, the trends in free energy as a function of sequence were identical. Equilibrium binding by BamHI restriction endonuclease to the 22-base pair DNAs was also investigated. The central eight base pairs of all six molecules, 5'-A-GGATCC-A-3', contained a BamHI recognition sequence bounded by A-T base pairs. Magnesium free binding assays were performed by titering BamHI against a constant concentration of each of the deoxyoligonucleotide substrates and analyzing reaction products by gel retardation. Binding isotherms of the total amount of bound DNA versus protein concentration were constructed which provided semiquantitative estimates of the equilibrium dissociation constants for dissociation of BamHI from the six DNA oligomers. Dissociation constants ranged from 0.5 x 10(-)(9) to 12.0 x 10(-)(9) M with corresponding binding free energies of -12.5 to -10.6 (+/-0. 1) kcal/mol. An inverse relationship is found when binding and stability are compared.

Melting studies of short DNA hairpins containing the universal base 5-nitroindole.

Effects of the universal base 5-nitroindole on the thermodynamic stability of DNA hairpins having a 6 bp stem and four base loops were investigated by optical absorbance and differential scanning calorimetry techniques. Melting studies were conducted in buffer containing 115 mM Na(+). Five different modified versions of DNA hairpins containing a 5-nitroindole base or bases substituted at different positions in the stem and loop regions were examined. Thermo-dynamic parameters of the melting transitions estimated from a two-state analysis of optical melting curves and measured directly by calorimetry revealed that the presence of 5-nitroindole bases in the duplex stem or loop regions of short DNA hairpins significantly affects both their enthalpic and entropic melting components in a compensating manner, while the transition free energy varies linearly with the transition temperature. The calorimetrically determined enthalpy and entropy values of the modified hairpins were considerably smaller (43-53%) than the two-state optical parameters, suggesting that solvent effects may be significant in the melting processes of these hairpins. Results of circular dichroism measurements also revealed slight differences between the modified hairpins and the control in both the duplex and melted states, suggesting subtle structural differences between the control and DNA hairpins containing a 5-nitroindole base or bases.

Melting studies of short DNA hairpins: influence of loop sequence and adjoining base pair identity on hairpin thermodynamic stability.

Spectroscopic and calorimetric melting studies of 28 DNA hairpins were performed. These hairpins form by intramolecular folding of 16 base self-complementary DNA oligomer sequences. Sequence design dictated that the hairpin structures have a six base pair duplex linked by a four base loop and that the first five base pairs in the stem are the same in every molecule. Only loop sequence and identity of the duplex base pair closing the loop vary for the set of hairpins. For these DNA samples, melting studies were carried out to investigate effects of the variables on hairpin stability. Stability of the 28 oligomers was ascertained from their temperature-induced melting transitions in buffered 115 mM Na(+) solvent, monitored by ultraviolet absorbance and differential scanning calorimetry (DSC). Experiments revealed the melting temperatures of these molecules range from 32.4 to 60.5 degrees C and are concentration independent over strand concentrations of 0.5 to 260 microM; thus, as expected for hairpins, the melting transitions are apparently unimolecular. Model independent thermodynamic transition parameters, DeltaH(cal), DeltaS(cal), and DeltaG(cal), were determined from DSC measurements. Model dependent transition parameters, DeltaH(vH), DeltaS(vH), and DeltaG(vH) were estimated from a van't Hoff (two-state) analysis of optical melting transitions. Results of these studies reveal a significant sequence dependence to DNA hairpin stability. Thermodynamic parameters evaluated by either procedure reveal the transition enthalpy, DeltaH(cal) (DeltaH(vH)) can differ by as much as 20 kcal/mol depending on sequence. Similarly, values of the transition entropy DeltaS(cal) (DeltaS(vH)) can differ by as much as 60 cal/Kmol (eu) for different molecules. Differences in free energies DeltaG(cal) (DeltaG(vH)) are as large as 4 kcal/mol for hairpins with different sequences. Comparisons between the model independent calorimetric values and the thermodynamic parameters evaluated assuming a two-state model reveal that 10 of the 28 hairpins display non-two-state melting behavior. The database of sequence-dependent melting free energies obtained for the hairpins was employed to extract a set of n-n (nearest-neighbor) sequence dependent loop parameters that were able to reproduce the input data within error (with only two exceptions). Surprisingly, this suggests that the thermodynamic stability of the DNA hairpins can in large part be reasonably represented in terms of sums of appropriate nearest-neighbor loop sequence parameters.