Regiospecific inhibition of DNA duplication by antisense phosphate-methylated oligodeoxynucleotides.

A new synthesis route for long phosphate-methylated oligodeoxynucleotides is described, which were used as antisense inhibitors of the DNA replication. Phosphate-methylated oligomers hybridize more strongly with natural DNA than their natural analogues, due to the absence of electrostatic interstrand repulsions. Compared with phosphate-ethylated and methyl phosphonate systems, phosphate-methylated systems are preferable as antisense DNA, which was concluded from the high Tm values and sharp melting transitions of duplexes of phosphate-methylated and natural DNA. By using the Sanger dideoxy technique, it was shown that a complementary phosphate-methylated 18-mer can effectively and site-specifically block the DNA replication process at room temperature.     

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.     

The thermodynamic advantage of DNA oligonucleotide 'stacking hybridization' reactions: energetics of a DNA nick.

'Stacking hybridization reactions' wherein two or more short DNA oligomers hybridize in a contiguous tandem orientation onto a longer complementary DNA single strand have been employed to enhance a variety of analytical oligonucleotide hybridization schemes. If the short oligomers anneal in perfect head-to-tail register the resulting duplex contains a nick at every boundary between hybridized oligomers. Alternatively, if the short oligomers do not hybridize precisely in register, i.e. single strand regions on the longer strand are left unbound, gaps are formed between regions where short oligomers bind. The resulting gapped DNA duplexes are considerably less stable than their nicked duplex analogs. Formation of base pair stacking interactions between neighboring oligomers at the nicks that do not occur in gapped duplexes has been proposed as the source of the observed added stability. However, quantitative evidence supporting this hypothesis for DNA has not been reported. Until now, a direct comparison of the thermodynamics of DNA nicks versus DNA gaps has not been performed. In this communication we report such a comparison. Analysis of optical melting experiments in a well defined molecular context enabled quantitative evaluations of the relative thermodynamic difference between nicked and gapped DNA duplexes. Results of the analysis reveal that a nick may be energetically favored over a gap by at least 1.4 kcal/mol and perhaps as much as 2.4 kcal/mol. The presence of a 5'phosphate at a nick or gap fails to significantly affect their stabilities.     

Quantitative analysis of DNA secondary structure from solvent-accessible surfaces: the B- to Z-DNA transition as a model

Solvent structure and its interactions have been suggested to play a critical role in defining the conformation of polynucleotides and other macromolecules. In this work, we attempt to quantitate solvent effects on the well-studied conformational transition between right-handed B- and left-handed Z-DNA. The solvent-accessible surfaces of the hexamer sequences d(m5CG)3, d(CG)3, d(CA)3, and d(TA)3 were calculated in their B- and Z-DNA conformations. The difference in hydration free energies between the Z and the B conformations (delta delta GH(Z-B] was determined from these surfaces to be -0.494 kcal/mol for C-5 methylated d(CG), 0.228 kcal/mol for unmethylated d(CG), 0.756 kcal/mol for d(CA)-d(TG), and 0.896 kcal/mol for d(TA) dinucleotides. These delta delta GH(Z-B) values were compared to the experimental B- to Z-DNA transition energies of -0.56 kcal/mol that we measured for C-5 methylated d(CG), 0.69-1.30 kcal/mol reported for unmethylated d(CG), 1.32-1.48 kcal/mol reported for d(CA)-d(TG), and 2.3-2.4 kcal/mol for d(TA) dinucleotides. From this comparison, we found that the calculated delta delta GH(Z-B) of these dinucleotides could account for the previous observation that the dinucleotides were ordered as d(m5CG) greater than d(CG) greater than d(CA)-d(TG) greater than d(TA) in stability as Z-DNA. Furthermore, we predicted that one of the primary reasons for the inability of d(TA) sequences to form Z-DNA results from a decrease in exposed hydrophilic surfaces of adjacent base pairs due to the C-5 methyl group of thymine; thus, d(UA) dinucleotides should be more stable as Z-DNA than the analogous d(TA) dinucleotides.(ABSTRACT TRUNCATED AT 250 WORDS)     

Effect of thiazole orange doubly labeled thymidine on DNA duplex formation

Nucleic acid oligonucleotides are widely used in hybridization experiments for specific detection of complementary nucleic acid sequences. For design and application of oligonucleotides, an understanding of their thermodynamic properties is essential. Recently, exciton-controlled hybridization-sensitive fluorescent oligonucleotides (ECHOs) were developed as uniquely labeled DNA oligomers containing commonly one thymidine having two covalently linked thiazole orange dye moieties. The fluorescent signal of an ECHO is strictly hybridization-controlled, where the dye moieties have to intercalate into double-stranded DNA for signal generation. Here we analyzed the hybridization thermodynamics of ECHO/DNA duplexes, and thermodynamic parameters were obtained from melting curves of 64 ECHO/DNA duplexes measured by ultraviolet absorbance and fluorescence. Both methods demonstrated a substantial increase in duplex stability (ΔΔG°(37) ～ -2.6 ± 0.7 kcal mol(-1)) compared to that of DNA/DNA duplexes of the same sequence. With the exception of T·G mismatches, this increased stability was mostly unaffected by other mismatches in the position opposite the labeled nucleotide. A nearest neighbor model was constructed for predicting thermodynamic parameters for duplex stability. Evaluation of the nearest neighbor parameters by cross validation tests showed higher predictive reliability for the fluorescence-based than the absorbance-based parameters. Using our experimental data, a tool for predicting the thermodynamics of formation of ECHO/DNA duplexes was developed that is freely available at http://genome.gsc.riken.jp/echo/thermodynamics/ . It provides reliable thermodynamic data for using the unique features of ECHOs in fluorescence-based experiments.     

The thermal stability of DNA fragments with tandem mismatches at a d(CXYG).d(CY'X'G) site.

Temperature-Gradient Gel Electrophoresis (TGGE) was employed to determine the thermal stabilities of 28 DNA fragments, 373 bp long, with two adjacent mismatched base pairs, and eight DNAs with Watson-Crick base pairs at the same positions. Heteroduplex DNAs containing two adjacent mismatches were formed by melting and reannealing pairs of homologous 373 bp DNA fragments differing by two adjacent base pairs. Product DNAs were separated based on their thermal stability by parallel and perpendicular TGGE. The polyacrylamide gel contained 3.36 M urea and 19.2 % formamide to lower the DNA melting temperatures. The order of stability was determined in the sequence context d(CXYG).d(CY'X'G) where X.X' and Y.Y" represent the mismatched or Watson-Crick base pairs. The identity of the mismatched bases and their stacking interactions influence DNA stability. Mobility transition melting temperatures (T u) of the DNAs with adjacent mismatches were 1.0-3.6 degrees C (+/-0.2 degree C) lower than the homoduplex DNA with the d(CCAG).d(CTGG) sequence. Two adjacent G.A pairs, d(CGAG).d(CGAG), created a more stable DNA than DNAs with Watson-Crick A.T pairs at the same sites. The d(GA).d(GA) sequence is estimated to be 0.4 (+/-30%) kcal/mol more stable in free energy than d(AA).d(TT) base pairs. This result confirms the unusual stability of the d(GA).d(GA) sequence previously observed in DNA oligomers. All other DNAs with adjacent mismatched base pairs were less stable than Watson-Crick homoduplex DNAs. Their relative stabilities followed an order expected from previous results on single mismatches. Two homoduplex DNAs with identical nearest neighbor sequences but different next-nearest neighbor sequences had a small but reproducible difference in T u value. This result indicates that sequence dependent next neighbor stacking interactions influence DNA stability.     

Spectroscopic and calorimetric investigation on the DNA triplex formed by d(CTCTTCTTTCTTTTCTTTCTTCTC) and d(GAGAAGAAAGA) at acidic pH.

The equimolar mixture of d(CTCTTCTTTCTTTTCTTTCTTCTC) (dY24) and d(GAGAAGAAAGA) (dR11) [designated (dY24).(dR11)], forms at pH = 5 a DNA triplex, which mimicks the H-DNA structure. The DNA triplex was identified by the following criteria: (i) dY24 and dR11 co-migrate in a poly-acrylamide gel, with a mobility and a retardation coefficient comparable to those observed for an 11-triad DNA triplex, previously characterized in our laboratories (1); (ii) the intercalator ethidium bromide shows a poor affinity for (dR11).(dY24) at pH = 5, and a high affinity at pH = 8; (iii) the (dR11).(dY24) mixture is not a substrate for DNase I at pH = 5; (iv) the CD spectrum of (dR11).(dY24), at pH = 5, is consistent with those previously reported for triple-stranded DNA. The (dR11).(dY24) mixture exhibits a thermally induced co-operative transition, which appears to be monophasic, reversible and concentration dependent. This transition is attributed to the disruption of the DNA triplex into single strands. The enthalpy change of the triplex-coil transition was measured by DSC (delta Hcal = 129 +/- 6 kcal/mol) and, assuming a two-state model, by analysis of UV-denaturation curves (average of two methods delta HUV = 137 +/- 13 kcal/mol). Subtracting from delta Hcal of triplex formation the contributions due to the Watson-Crick helix and to the protonation of the C-residues, we found that each pyrimidine binding into the major groove of the duplex, through a Hoogsteen base pair, is accompanied by an average delta H = -5.8 +/- 0.6 kcal/mol. The effect on the stability of the (dR11).(dY24) triplex due to the substitution of a T:A:T triad with a T:T:T one was also investigated.     

Coupling of sequential transitions in a DNA double hairpin: energetics, ion binding, and hydration

In an effort to evaluate the relative contributions of sequence, ion binding, and hydration to the thermodynamic stability of nucleic acids, we have investigated the melting behavior of a double hairpin and that of its component single hairpins. Temperature-dependent UV absorption and differential scanning calorimetry techniques have been used to characterize the helix-coil transitions of three deoxyoligonucleotides: d(GTACT5GTAC), d(GCGCT5GCGC), and d(GCGCT5GCGCGTACT5GTAC). The first two oligomers melt with transition temperatures equal to 28 and 69 degrees C, respectively, in 10 mM dibasic sodium phosphate at pH 7.0. The Tm's are independent of strand concentration, strongly indicating the presence of single-stranded hairpin structures at low temperatures. The third oligomer, with a sequence corresponding to the joined sequences of the first two oligomers, melts with two apparently independent monomolecular transitions with Tm's of 41 and 69 degrees C. These transitions correspond to the melting of a double hairpin. In the salt range of 10-100 mM in NaCl, we obtain average enthalpies of 24 and 38 kcal/mol for the transitions in the single-hairpin molecules. Each transition in the double hairpin has an enthalpy of 32 kcal/mol. In addition, dtm/d log [Na+] for the transitions are 4.1 and 4.7 degrees C for the single hairpins and 12.6 and 11.2 degrees C for each transition in the double hairpin. The differential ion binding parameter between the double hairpin and that of the sum of single hairpins is roughly equal to 1.1 mol of Na+ ions/mol of double hairpin and is consistent with an increase in the electrostatic behavior of the stem phosphates of this molecule.     

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.     

Thermodynamic and base-pairing studies of matched and mismatched DNA dodecamer duplexes containing cis-syn, (6-4) and Dewar photoproducts of TT.

Cis-syn dimers, (6-4) products and their Dewar valence isomers are the major photoproducts of DNA and have different mutagenic properties and rates of repair. To begin to understand the physical basis for these differences, the thermal stability and base pairing properties of the corresponding photoproducts of the TT site in d(GAGTATTATGAG) were investigated. The (6-4) and Dewar products destabilize the duplex form by approximately 6 kcal/mol of free energy at 37 degreesC relative to the parent, whereas a cis-syn dimer only destabilizes the duplex form by 1.5 kcal/mol. Duplexes with G opposite the 3'-T of the (6-4) and Dewar products are more stable than those with A by approximately 0.4 kcal/mol, whereas the cis-syn dimer prefers A over G by 0.7 kcal/mol. Proton NMR suggests that wobble base pairing takes place between the 3'-T of the cis-syn dimer and an opposed G, whereas there is no evidence of significant H-bonding between these two bases in the (6-4) product. The thermodynamic and H-bonding data for the (6-4) product are consistent with a 4 nt interior loop structure which may facilitate flipping of the photoproduct in and out of the helix.     

Comparison of thermodynamic stabilities between PNA (peptide nucleic acid)/DNA hybrid duplexes and DNA/DNA duplexes.

We have examined quantitatively stabilities of PNA/DNA hybrid duplexes with identical nearest-neighbor base pairs and compared stabilities between PNA/DNA and DNA/DNA. The average difference of stabilization energy of the short PNA/DNA was 0.9 kcal mol(-1), which suggests that the stability of the hybrids with identical nearest-neighbor base pairs can be predicted with the nearest-neighbor model as well as those of nucleic acid duplexes.     

The conformational study of two carbocyclic nucleosides: why carbocyclic nucleic acids (CarNAs) form more stable duplexes with RNA than DNA does

Replacement of the furanose moiety of DNA with a cyclopentane ring produces a modified sugar analogue: Carbocyclic nucleic acid (CarNA). UV melting-temperature experiments demonstrate that the incorporation of 2'-deoxycarbaguanosine ((c)G) and 2'-deoxyaristeromycin ((c)A) of carbocyclic nucleosides into a DNA strand increases the stability of the CarNA/RNA hybrid. Circular Dichroism (CD) study indicates that the CarNA/RNA hybrid adopts an A-like conformation. To elucidate the molecular basis of the increased stability of the CarNA/RNA, the conformation of (c)G and (c)A were examined by (1)H NMR conformational analysis of (3)J(HH) coupling constants and ab initio molecular orbital (MO) calculations. These results show that the populations of N-type of (c)G and (c)A are higher than those of dG and dA, respectively, at different temperatures [For example, 37% (N%) of (c)G vs. 28%of dG, 36% (N%) of (c)A vs. 25% of dA at 278 K], which suggest that the cyclopentane rings of (c)G and (c)A prefer the N-type conformation in two-state N-S pseudorotional equilibrium in comparison with the furanose rings of dG and dA. The DeltaH degrees of (c)G (DeltaH degrees = - 0.43 kcal mol(-1)) and (c)A (DeltaH degrees = - 0.41kcal mol(-1)) are lower than that of dG (dG = - 1.8 kcal mol(-1)) and dA (dA = - 1.0 kcal mol(-1)), respectively, which suggest that the gauche effect in the (c)A and (c)G driving N-S pseudorotional equilibrium to S-type is reduced by replacement of the 4'-oxygen by a CH(2) group. These results suggest that the preferred N-type of the (c)G and (c)A leads to the A-like conformation, which contributes to the stability of CarNA/RNA hybrid.     

Thermodynamic parameters for DNA sequences with dangling ends

The thermodynamic contributions to duplex formation of all 32 possible single-nucleotide dangling ends on a Watson-Crick pair are reported. In most instances, dangling ends are stabilizing with free energy contributions ranging from +0.48 (GT(A)) to-0.96 kcal/mol (). In comparison, Watson-Crick nearest-neighbor increments range from -0. 58 (TA/AT) to -2.24 (GC/CG) kcal/mol. Hence, in some cases, a dangling end contributes as much to duplex stability as a Watson-Crick A-T base pair. The implications of these results for DNA probe design are discussed. Analysis of the sequence dependence of dangling-end stabilities show that the nature of the closing base pair largely determines the stabilization. For a given closing base pair, however, adenine dangling ends are always more or equally as stable as the other dangling nucleotides. Moreover, 5' dangling ends are more or equally as stabilizing as their 3' counterparts. Comparison of DNA with RNA dangling-end motifs shows that DNA motifs with 5' dangling ends contribute to stability equally or more than their RNA counterparts. Conversely, RNA 3' dangling ends contribute to stability equally or more than their DNA counterparts. This data set has been incorporated into a DNA secondary structure prediction algorithm (DNA MFOLD) (http://mfold2.wustl.edu/mfold/dna/for m1.cgi) as well as a DNA hybridization prediction algorithm (HYTHERtrade mark) (http://jsl1.chem.wayne.edu/Hyther/hythermenu .html).     

Spectroscopic analysis of DNA base-pair opening by Escherichia coli RNA polymerase. Temperature and ionic strength effects

The interaction of Escherichia coli RNA polymerase with poly[d(A-T)] and poly[d-(I-C)] was studied by difference absorption spectroscopy at temperatures, from 5 to 45 degrees C in the absence and presence of Mg2+. The effect of KCl concentration, at a fixed temperature, was studied from 12.5 to 400 mM. Difference absorption experiments permitted calculation of the extent of DNA opening induced by RNA polymerase and estimation of the equilibrium constant associated with the isomerization from a closed to an open RNA polymerase-DNA complex. delta H0 and delta S0 for the closed-to-open transition with poly[d(A-T)] or poly[d(I-C)] complexed with RNA polymerase are significantly lower than the values associated with the helix-to-coil transition for the free polynucleotides. For the RNA polymerase complexes with poly[d(A-T)] and poly[d(I-C)] in 50 mM KCl, delta H0 approximately 15-16 kcal/mol (63-67 kJ/mol) and delta S0 approximately 50-57 cal/K per mol (209-239 J/K per mol). The presence of Mg2+ does not change these parameters appreciably for the RNA polymerase-poly[d(A-T)] complex, but for the RNA polymerase-poly[d(I-C)] complex in the presence of Mg2+, the delta H0 and delta S0 values are larger and temperature-dependent, with delta H0 approximately 22 kcal/mol (92 kJ/mol) and delta S0 approximately 72 cal/K per mol (approx. 300 J/K per mol) at 25 degrees C, and delta Cp0 approximately 2 kcal/K per mol (approx. 8.3 kJ/K per mol). The circular dichroism (CD) changes observed for helix opening induced by RNA polymerase are qualitatively consistent with the thermally induced changes observed for the free polynucleotides, supporting the difference absorption method. The salt-dependent studies indicate that two monovalent cations are released upon helix opening. For poly[d(A-T)], the temperature-dependence of enzyme activity correlates well with the helix opening, implying this step to be the rate-determining step. In the case of poly[d(I-C)], the same is not true, and so the rate-determining step must be a process subsequent to helix opening.     

Enthalpy and heat capacity changes for formation of an oligomeric DNA duplex: interpretation in terms of coupled processes of formation and association of single-stranded helices.

The thermodynamics of self-assembly of a 14 base pair DNA double helix from complementary strands have been investigated by titration (ITC) and differential scanning (DSC) calorimetry, in conjunction with van't Hoff analysis of UV thermal scans of individual strands. These studies demonstrate that thermodynamic characterization of the temperature-dependent contributions of coupled conformational equilibria in the individual "denatured" strands and in the duplex is essential to understand the origins of duplex stability and to derive stability prediction schemes of general applicability. ITC studies of strand association at 293 K and 120 mM Na+ yield an enthalpy change of -73 +/- 2 kcal (mol of duplex)-1. ITC studies between 282 and 312 K at 20, 50, and 120 mM Na+ show that the enthalpy of duplex formation is only weakly salt concentration-dependent but is very strongly temperature-dependent, decreasing approximately linearly with increasing temperature with a heat capacity change (282-312 K) of -1.3 +/- 0.1 kcal K-1 (mol of duplex)-1. From DSC denaturation studies in 120 mM Na+, we obtain an enthalpy of duplex formation of -120 +/- 5 kcal (mol of duplex)-1 and an estimate of the corresponding heat capacity change of -0.8 +/- 0.4 kcal K-1 (mol of duplex)-1 at the Tm of 339 K. van't Hoff analysis of UV thermal scans on the individual strands indicates that single helix formation is noncooperative with a temperature-independent enthalpy change of -5.5 +/- 0.5 kcal at 120 mM Na+. From these observed enthalpy and heat capacity changes, we obtain the corresponding thermodynamic quantities for two fundamental processes: (i) formation of single helices from disordered strands, involving only intrastrand (vertical) interactions between neighboring bases; and (ii) formation of double helices by association (docking) of single helical strands, involving interstrand (horizontal and vertical) interactions. At 293 K and 120 mM Na+, we calculate that the enthalpy change for association of single helical strands is approximately -64 kcal (mol of duplex)-1 as compared to -210 kcal (mol of duplex)-1 calculated for duplex formation from completely unstructured single strands and to the experimental ITC value of -73 kcal (mol of duplex)-1. The intrinsic heat capacity change for association of single helical strands to form the duplex is found to be small and positive [ approximately 0.1 kcal K-1 (mol of duplex)-1], in agreement with the result of a surface area analysis, which also predicts an undetectably small heat capacity change for single helix formation.     

The mechanism of recognition of templates by DNA polymerases from pro- and eukaryotes as revealed by affinity modification data.

Pt(2+)-containing derivatives of oligodeoxyribonucleotides were used to evaluate the ligand affinity to the template sites of Klenow fragment of DNA polymerase I from E. coli and DNA polymerase alpha from human placenta. The values of Kd and Gibb's energy (delta G degree) for the complexes of oligodeoxyribonucleotides and their derivatives with the template sites of these enzymes were determined from the effects protecting the enzyme from inactivation by Pt(2+)-containing oligonucleotides. Kd and delta G degree values of the complexes made by DNA polymerases and orthophosphate, triethylphosphate, d(pC)n, d(pT)n, d(pG)n, d(pA)n (where n = 1-25), heterooligonucleotides of various length and structure, and oligothymidylates with partially and completely ethylated internucleotide phosphates were evaluated. The obtained data enabled us to suggest 19-20 mononucleotide units of the template to interact with the protein. Only one template internucleotide phosphate forms a Me(2+)-dependent electrostatic contact (delta G = -1.1...-1.7 kcal/mol) and a hydrogen bond (delta G = -4.4...-4.9 kcal/mol) with the enzyme. It is likely that the mononucleoside units of the template form hydrophobic contacts with the enzymes. The efficiency of such interaction changes with the hydrophobicity of the bases: C less than T less than G approximately A. For both homo- and heterooligonucleotides the contributions of nucleoside units to the affinity of the templates to the enzymes is due to the complementary interactions with the primers. A hypothetical model for the template-primer interaction with DNA polymerases is suggested.     

Binding free energy calculation for duocarmycin/DNA complex based on the QPLD-derived partial charge model

The 3ns unrestrained MD simulations were carried out on the DNA/duocarmycin complex based on (1) the classic RESP charge model, and (2) the QM-polarized ligand docking (QPLD)-based charge model. The RMSDs of the trajectories and the DeltaG(bind) of the QPLD model perform much better than the RESP model, with the DeltaG(bind) estimation for QPLD model (-16.11 kcal/mol) versus DeltaG(bind) estimation for RESP model (-10.05 kcal/mol).     

Melting of a self-complementary DNA minicircle. Comparison of optical melting theory with exchange broadening of the nuclear magnetic resonance spectrum

Melting curves are calculated for the 16-base-pair duplex DNA sequence 5' GTATCCGTACGGATAC 3' linked on the ends by TTTT single-strand loops. The equilibrium statistical thermodynamic theory of DNA melting is modified to include effects of end-loops on the melting transition. An excellent fit of the experimental melting curve in 0.2 M-NaCl is obtained using two adjustable parameters, one for end-loop formation and the other for formation of the complete 40-base single-strand loop. The best-fit calculated melting curve permits evaluation of these parameters. The free energy to close a TTTT end-loop is 2.12 kcal/mol (1 cal = 4.184 J). A TTTT end-loop or hairpin loop is significantly more stable than an internal loop of comparable size sandwiched between two helical regions, even after allowing for the different stacking contributions. Reasons for this increased stability are presented. The loop free energy of the 40-base single-strand open minicircle is evaluated to be +1.27 kcal/mol, thus favoring the melting of two end-loops into the large open minicircle. The present results are compared with those of others for d(T-A) oligomers. The sequence TTTT forms a more stable end-loop, or hairpin, than TATA by about 2.0 kcal/mol. Theoretical rate constants for the proton-transfer step in the standard hydrogen-exchange model are calculated by extending the theory of diffusion-controlled reactions to take account of the electrostatic potential of the DNA. The predicted ratios of rate constants for different pairs of catalysts exchanging an A.T proton agree satisfactorily with the available experimental data for a 14-base-pair linear duplex, which confirms the diffusion-control of the proton-transfer step. Data presented here for the 16 base-pair duplex of the minicircle are consistent with catalysis-limited exchange in which the proton-transfer step is likewise diffusion-controlled. Under catalysis-limited conditions, the imino proton exchange rates are predicted from the catalytic rate constants, prevailing buffer catalyst concentrations, and the equilibrium constants to form the unstacked open state of optical melting theory. The observed exchange rates of the A.T base-pairs show no sign of the strong predicted end-melting trend, and exceed the predicted values by factors of 10 to 400. Moreover, the succession of "melting" in the nuclear magnetic resonance line-broadening deviates from that predicted by optical melting theory.(ABSTRACT TRUNCATED AT 400 WORDS)