Thermodynamic and hydration effects for the incorporation of a cationic 3-aminopropyl chain into DNA

The introduction of cationic 5-(ω-aminoalkyl)-2′-deoxypyrimidines into duplex DNA has been shown to induce DNA bending. In order to understand the energetic and hydration contributions for the incorporation of a cationic side chain in DNA a combination of spectroscopy, calorimetry and density techniques were used. Specifically, the temperature unfolding and isothermal formation was studied for a pair of duplexes with sequence d(CGTAGUCG TGC)/d(GCACGACTACG), where U represents 2′-deoxyuridine (‘control’) or 5-(3-aminopropyl)-2′-deoxyuridine (‘modified’). Continuous variation experiments confirmed 1:1 stoichiometries for each duplex and the circular dichroism spectra show that both duplexes adopted the B conformation. UV and differential scanning calorimetry melting experiments reveal that each duplex unfolds in two-state transitions. In low salt buffer, the ‘modified’ duplex is more stable and unfolds with a lower endothermic heat and lower release of counterion and water. This electrostatic stabilization is entropy driven and disappears at higher salt concentrations. Complete thermodynamic profiles at 15°C show that the favorable formation of each duplex results from the compensation of a favorable exothermic heat with an unfavorable entropy contribution. However, the isothermal profiles yielded a differential enthalpy of 8.8 kcal/mol, which is 4.3 kcal/mol higher than the differential enthalpy observed in the unfolding profiles. This indicates that the presence of the aminopropyl chain induces an increase in base stacking interactions in the modified single strand and a decrease in base stacking interactions in the modified duplex. Furthermore, the formation of the ‘control’ duplex releases water while the ‘modified’ duplex takes up water. Relative to the control duplex, formation of the modified duplex at 15°C yielded a marginal differential ΔG° term, positive ΔΔH_ITC–Δ(TΔS) compensation, negative ΔΔV and a net release of counterions. The opposite signs of the differential enthalpy–entropy compensation and differential volume change terms show a net uptake of structural water around polar and non-polar groups. This indicates that incorporation of the aminopropyl chain induces a higher exposure of aromatic bases to the solvent, which may be consistent with a small and local bend in the ‘modified’ duplex.     

DNA hairpins: fuel for autonomous DNA devices

We present a study of the hybridization of complementary DNA hairpin loops, with particular reference to their use as fuel for autonomous DNA devices. The rate of spontaneous hybridization between complementary hairpins can be reduced by increasing the neck length or decreasing the loop length. Hairpins with larger loops rapidly form long-lived kissed complexes. Hairpin loops may be opened by strand displacement using an opening strand that contains the same sequence as half of the neck and a "toehold" complementary to a single-stranded domain adjacent to the neck. We find loop opening via an external toehold to be 10-100 times faster than via an internal toehold. We measure rates of loop opening by opening strands that are at least 1000 times faster than the spontaneous interaction between hairpins. We discuss suitable choices for loop, neck, and toehold length for hairpin loops to be used as fuel for autonomous DNA devices.     

Solvent isotope effect on thermodynamics of hydration

Partial molar heat capacities of five linear alcohols (methanol, ethanol, n-propanol, n-butanol, n-pentanol) and five N-substituted amides (n-propionamide, N-methylformamide, N-methylacetamide, N-methylpropionamide, N-ethylacetamide) in aqueous D(2)O solution have been measured at 25 degrees C. The heat capacities of transfer of these compounds from H(2)O to D(2)O were calculated using previously reported (Makhatadze et al., Biophys. Chem. 64 (1997) 93) values of partial heat capacities of alcohols and amides in aqueous H(2)O solutions. It is shown that the sign and magnitude of the heat capacity change upon transfer from H(2)O to D(2)O depends on the relative amount of polar and non-polar solvent accessible surface areas of solute. Analysis shows that transfer of non-polar surface from H(2)O to D(2)O is accompanied by a positive heat capacity change. In contrast, transfer of polar surface from H(2)O to D(2)O occurs with negative heat capacity change. Estimates show that the solvent isotope effect on the heat capacity changes upon protein unfolding can be predicted using the changes of the polar and non-polar surface area changes upon protein unfolding and the transfer data of model compounds. Analysis of the thermodynamic functions of transfer of non-polar compounds from H(2)O to D(2)O shows puzzling behavior which contradicts current definitions of the hydrophobic effect.     

DNA replication fidelity: kinetics and thermodynamics

Mechanisms that control the fidelity of DNA replication are discussed. Data are reviewed for 3 steps in a fidelity pathway: nucleotide insertion, exonucleolytic proofreading, and extension from matched and mismatched 3′-primer termini. Fidelity mechanisms that involve predominately K_m discrimination, V_max discrimination, or a combination of the two are analyzed in the context of a simple model for fidelity. Each fidelity step is divided into 2 components, thermodynamics and kinetic. The thermodynamic component, which relates to free-energy differences between right and wrong base pair, is associated with a K_m discrimination mechanism for polymerase. The kinetic component, which represents the enzyme's ability to select bases for insertion and excision to achieve fidelity greater than that availablek from base pairing free-energy differences, is associated with a V_max discrimination mechanism for polymerase. Currently available fidelity data for nucleotide insertion and primer extension in the absence of proofreading appears to have relatively large K_m and small V_max components. An important complication can arise when analyzing data from polymerases containing an associated 3′-exonuclease activity. In the presence of proofreading, a V_max discrimination mechanisms is likely to occur, but this may be the result of two K_m discrimination mechanisms acting serially, one for nucleotide insertion and other for excision. Possible relationships between base pairing free energy differences measured in aqueous solution and those defined within the polymerase active cleft are considered in the context of the enzyme's ability to exclude water, at least partially, from the vicinity of its active site.     

Equilenin-derived DNA adducts to cytosine in DNA duplexes: structures and thermodynamics

The drug Premarin is the most widely used formula for hormone replacement therapy. However, long-term exposure to estrogens from the Premarin drug increases the risk of breast cancer. Equilin and equilenin, major components of Premarin, are predominantly metabolized to 4-hydroxyequilenin (4-OHEN). The quinoids produced by 4-OHEN oxidation react with dG, dA, and dC to form unusual stable cyclic bulky adducts, with four stereoisomers identified for each base adduct. The 4-OHEN-dC adducts are most predominant. They are mutagenic in vitro and have been found in human tumor tissue. We have carried out molecular modeling and molecular dynamics simulations to investigate structures and thermodynamics of the four 4-OHEN-dC stereoisomeric adducts in DNA duplexes. Our results show that the structure of each stereoisomer adduct in duplex DNA is specifically governed by its unique stereochemistry. The bulky adducts, with an obstructed Watson-Crick edge and an equilenin ring system near perpendicular to the damaged cytosine, are located in the B-DNA major or minor groove, with the modified cytosine in the syn or anti conformation, respectively. The DNA duplex structures are distorted, in terms of Watson-Crick pairing at and near the lesion, stacking interactions, and groove dimensions. Stereochemistry determines the orientation of the equilenin rings with respect to the 5'- to 3'-direction of the modified strand, as well as the positioning of the equilenin moiety's methyl and hydroxyl groups for each stereoisomer. The unusual structures and the stereochemical effects underlie their biological processing as miscoding DNA lesions whose mutagenic properties may contribute to breast cancer.     

Stable sheared A.C pair in DNA hairpins

Single-residue d(Pu1NPu2) (Pu1.Pu2=G.A, G.G or A.A) hairpin loops can be stably closed by sheared purine.purine pairs. These special motifs have been found in several important biological systems. We now extend these loop-closing base-pairs to a sheared purine. pyrimidine (A.C) pair at a neutral pH condition. High-resolution NMR spectroscopy, distance geometry, and molecular dynamics methods were used to study d(GTACANCGTAC) oligomers. Numerous idiosyncratic nuclear Overhauser enhancements, especially those across the A.C base-pair between C4NH2left and right arrow AH1', C4NH2left and right arrow AH2, and CH5left and right arrow AH2 proton pairs, clearly define the novel sheared nature of the closing A.C base-pair. This novel base-pair is possibly present in several biological systems and in two single-stranded DNA aptamers selected from oligonucleotide libraries.     

Immobilized DNA hairpins for assay of sequential breaking and joining of DNA backbones

Immobilized DNA hairpins are exploited in a novel approach to assay DNA ligases and nucleases. A fundamental characteristic of the assay is that a fluorophore at the remote terminus of the hairpin reports on the integrity of the DNA backbone. The functionality of the protocol is confirmed using ATP- and NAD+-dependent DNA ligases and the nicking enzyme N.BbvCIA. The assay format is amenable to high-throughput analysis and quantitation of enzyme activity, and it is shown to be in excellent agreement with the more laborious electrophoretic approaches that are widely used for such analyses. Significantly, the assay is used to demonstrate sequential breaking and rejoining of a specific nucleic acid. Thus, a simple platform for biochemically innovative studies of pathways in cellular nucleic acid metabolism is demonstrated.     

Atomic-scale analysis of the solvation thermodynamics of hydrophobic hydration.

Molecular dynamics simulations are used to model the transfer thermodynamics of krypton from the gas phase into water. Extra long, nanosecond simulations are required to reduce the statistical uncertainty of the calculated "solvation" enthalpy to an acceptable level. Thermodynamic integration is used to calculate the "solvation" free energy, which together with the enthalpy is used to calculate the "solvation" entropy. A comparison series of simulations are conducted using a single Lennard-Jones sphere model of water to identify the contribution of hydrogen bonding to the thermodynamic quantities. In contrast to the classical "iceberg" model of hydrophobic hydration, the favorable enthalpy change for the transfer process at room temperature is found to be due primarily to the strong van der Waals interaction between the solute and solvent. Although some stabilization of hydrogen bonding does occur in the solvation shell, this is overshadowed by a destabilization due to packing constraints. Similarly, whereas some of the unfavorable change in entropy is attributed to the reduced rotational motion of the solvation shell waters, the major component is due to a decrease in the number of positional arrangements associated with the translational motions.     

Interaction of cationic surfactants with DNA: a single-molecule study

The interaction of cationic surfactants with single dsDNA molecules has been studied using force-measuring optical tweezers. For hydrophobic chains of length 12 and greater, pulling experiments show characteristic features (e.g. hysteresis between the pulling and relaxation curves, force-plateau along the force curves), typical of a condensed phase (compaction of a long DNA into a micron-sized particle). Depending on the length of the hydrophobic chain of the surfactant, we observe different mechanical behaviours of the complex (DNA-surfactants), which provide evidence for different binding modes. Taken together, our measurements suggest that short-chain surfactants, which do not induce any condensation, could lie down on the DNA surface and directly interact with the DNA grooves through hydrophobic–hydrophobic interactions. In contrast, long-chain surfactants could have their aliphatic tails pointing away from the DNA surface, which could promote inter-molecular interactions between hydrophobic chains and subsequently favour DNA condensation.     

DNA phosphate crowding correlates with protein cationic side chain density and helical curvature in protein/DNA crystal structures

Sequence-specific binding of proteins to their DNA targets involves a complex spectrum of processes that often induce DNA conformational variation in the bound complex. The forces imposed by protein binding that cause the helical deformations are intimately interrelated and difficult to parse or rank in importance. To investigate the role of electrostatics in helical deformation, we quantified the relationship between protein cationic residue density (Cpc) and DNA phosphate crowding (Cpp). The correlation between Cpc and Cpp was then calculated for a subset of 58 high resolution protein–DNA crystal structures. Those structures containing strong Cpc/Cpp correlation (>±0.25) were likely to contain DNA helical curvature. Further, the correlation factor sign predicted the direction of helical curvature with positive (16 structures) and negative (seven structures) correlation containing concave (DNA curved toward protein) and convex (DNA curved away from protein) curvature, respectively. Protein–DNA complexes without significant Cpc/Cpp (36 structures) correlation (-0.25<0<0.25) tended to contain DNA without significant curvature. Interestingly, concave and convex complexes also include more arginine and lysine phosphate contacts, respectively, whereas linear complexes included essentially equivalent numbers of Lys/Arg phosphate contacts. Together, these findings suggest an important role for electrostatic interactions in protein–DNA complexes involving helical curvature.     

Incorporation of oxyethylene units between hydrocarbon chain and pseudoglyceryl backbone in cationic lipid potentiates gene transfection efficiency in the presence of serum.

Four novel cationic lipids with different numbers of oxyethylene units at the linkage region between the pseudoglyceryl backbone and the hydrocarbon chains have been synthesized and used as mixtures with 1,2-dioleoyl-L-alpha-glycero-3-phosphatidyl ethanolamine (DOPE) for liposome-mediated gene transfection. Incorporation of different numbers of oxyethylene (-CH(2)CH(2)O-) units between long hydrocarbon chain at the C-1 and C-2 positions of the pseudoglyceryl skeleton improved the transfection efficiency considerably compared to the one in which the chains were connected via simple ether links. A pronounced improvement in the gene transfer efficiency was observed with the unsymmetrical cationic lipid 3 in which the long hydrocarbon at the C-1 position of the pseudoglyceryl segment is connected via two (-CH(2)CH(2)O-) units. Notably, the transfection ability of lipid 3 with DOPE in the presence of serum was significantly greater than LIPOFECTAMINE. This suggests that introduction of oxyethylene units between long hydrocarbon chains at the C-1 and C-2 positions of the pseudoglyceryl skeleton provides a novel strategy to achieve efficient gene transfer, especially in conditions where the presence of serum is critical.     

Incorporation in lipid bilayers of a large conductance cationic channel from mitochondrial membranes.

Membranes from subcellular fractions of adrenal medulla were incorporated in phospholipid bilayers formed at the tip of microelectrodes. Current fluctuations recorded in the presence of a transmembrane potential revealed the existence of a voltage-dependent channel of large conductance. This channel is characterized by fast kinetics and four conductance levels separated by jumps of 100, 220 and 220 pS in 150 mM NaCl. It is permeant to Na+,K+, tetraethylammonium, Cl- and acetate and has some cation selectivity. Exposure to trypsin or pronase abolished the voltage-dependence. Upon subcellular fractionation, the activity was found to be associated with mitochondria. A similar activity was observed in mitochondrial fractions from other organs. By its kinetics, its selectivity and its potential-dependence, this channel differs from the voltage-dependent anion channel of outer mitochondrial membranes.     

Extracting kinetics from single-molecule force spectroscopy: nanopore unzipping of DNA hairpins

Single-molecule force experiments provide powerful new tools to explore biomolecular interactions. Here, we describe a systematic procedure for extracting kinetic information from force-spectroscopy experiments, and apply it to nanopore unzipping of individual DNA hairpins. Two types of measurements are considered: unzipping at constant voltage, and unzipping at constant voltage-ramp speeds. We perform a global maximum-likelihood analysis of the experimental data at low-to-intermediate ramp speeds. To validate the theoretical models, we compare their predictions with two independent sets of data, collected at high ramp speeds and at constant voltage, by using a quantitative relation between the two types of measurements. Microscopic approaches based on Kramers theory of diffusive barrier crossing allow us to estimate not only intrinsic rates and transition state locations, as in the widely used phenomenological approach based on Bell's formula, but also free energies of activation. The problem of extracting unique and accurate kinetic parameters of a molecular transition is discussed in light of the apparent success of the microscopic theories in reproducing the experimental data.     

Structure and hydration of BamHI DNA recognition site: a molecular dynamics investigation

The results of a 3-ns molecular dynamics simulation of the dodecamer duplex d(TATGGATCCATA)(2) recognized by the BamHI endonuclease are presented here. The DNA has been simulated as a flexible molecule using an AMBER force field and the Ewald summation method, which eliminates the undesired effects of truncation and permits evaluation of the full effects of electrostatic forces. The starting B conformation evolves toward a configuration quite close to that observed through x-ray diffraction in its complex with BamHI. This configuration is fairly stable and the Watson-Crick hydrogen bonds are well maintained over the simulation trajectory. Hydration analysis indicates a preferential hydration for the phosphate rather than for the ester oxygens. Hydration shells in both the major and minor groove were observed. In both grooves the C-G pairs were found to be more hydrated than A-T pairs. The "spine of hydration" in the minor groove was clear. Water residence times are longer in the minor groove than in the major groove, although relatively short in both cases. No special long values are observed for sites where water molecules were observed by x-ray diffraction, indicating that water molecules having a high probability of being located in a specific site are also fast-exchanging.     

Drug recognition of a DNA single strand break: nogalamycin intercalation between coaxially stacked hairpins.

Two DNA hairpin motifs (5'-GCGAAGC-3' and 5'-ACGA AGT-3'), both stabilized by a 5'-GAA loop, have been used to design novel intramolecular double hairpin structures (5'-GCGAAGCACGAAGT-3' and 5'-ACGAAGTGCG AAGC-3') in which coaxial stacking of the two hairpin components generates a double-stranded stem region effectively with a single-strand break in the middle of the sequence at either the TG or CA step between unconnected 3' and 5' terminal bases. We have investigated by NMR the conformation and dynamics of the DNA at the strand break site. We show that mutual stacking significantly enhances the stability of each hairpin. Further, the anthracycline antibiotic nogalamycin binds cleanly to the 5'-TG (5'-CA) site formed by the mutually stacked hairpins despite the break in the sugar-phosphate backbone on one strand. The complex resembles the structure of nogalamycin-DNA complexes with the drug bound at 5'-TG sites in intact duplex sequences, with pi-stacking interactions probably the single dominant stabilizing interaction.     

A thermodynamic study of unusually stable RNA and DNA hairpins.

About 70% of the RNA tetra-loop sequences identified in ribosomal RNAs from different organisms fall into either (UNCG) or (GNRA) families (where N = A, C, G, or U; and R = A or G). RNA hairpins with these loop sequences form unusually stable tetra-loop structures. We have studied the RNA hairpin GGAC(UUCG)GUCC and several sequence variants to determine the effect of changing the loop sequence and the loop-closing base pair on the thermodynamic stability of (UNCG) tetra-loops. The hairpin GGAG(CUUG)CUCC with the conserved loop G(CUUG)C was also unusually stable. We have determined melting temperatures (Tm), and obtained thermodynamic parameters for DNA hairpins with sequences analogous to stable RNA hairpins with (UNCG), C(GNRA)G, C(GAUA)G, and G(CUUG)C loops. DNA hairpins with (TTCG), (dUdUCG), and related sequences in the loop, unlike their RNA counterparts, did not form unusually stable hairpins. However, DNA hairpins with the consensus loop sequence C(GNRA)G were very stable compared to hairpins with C(TTTT)G or C(AAAA)G loops. The C(GATA)G and G(CTTG)C loops were also extra stable. The relative stabilities of the unusually stable DNA hairpins are similar to those observed for their RNA analogs.     

Predicting RNA folding thermodynamics with a reduced chain representation model

Based on the virtual bond representation for the nucleotide backbone, we develop a reduced conformational model for RNA. We use the experimentally measured atomic coordinates to model the helices and use the self-avoiding walks in a diamond lattice to model the loop conformations. The atomic coordinates of the helices and the lattice representation for the loops are matched at the loop-helix junction, where steric viability is accounted for. Unlike the previous simplified lattice-based models, the present virtual bond model can account for the atomic details of realistic three-dimensional RNA structures. Based on the model, we develop a statistical mechanical theory for RNA folding energy landscapes and folding thermodynamics. Tests against experiments show that the theory can give much more improved predictions for the native structures, the thermal denaturation curves, and the equilibrium folding/unfolding pathways than the previous models. The application of the model to the P5abc region of Tetrahymena group I ribozyme reveals the misfolded intermediates as well as the native-like intermediates in the equilibrium folding process. Moreover, based on the free energy landscape analysis for each and every loop mutation, the model predicts five lethal mutations that can completely alter the free energy landscape and the folding stability of the molecule.     

On the hydration of DNA. I. Preferential hydration and stability of DNA in concentrated trifluoroacetate solution

The hydration of DNA is an important factor in the stability of its secondary structure. Methods for measuring the hydration of DNA in solution and the results of various techniques are compared and discussed critically. The buoyant density of native and denatured T-7 bacteriophage DNA in potassium trifluoroacetate (KTFA) solution has been measured as a function of temperature between 5 and 50°C. The buoyant density of native DNA increased linearly with temperature, with a dependence of (2.3 ± 0.5) × 10^?4 g/cc-°C. DNA which has been heat denatured and quenched at 0°C in the salt solution shows a similar dependence of buoyant density on temperature at temperatures far below the T_m, and above the T_m. However, there is an inflection region in the buoyant density versus T curve over a wide range of temperatures below the T_m. Optical density versus temperature studies showed that this is due to the. inhibition by KTFA of recovery of secondary structure on quenching. If the partial specific volume is assumed to be the same for native and denatured DNA, the loss of water of hydration on denaturation is calculated to be about 20% in KTFA at a water activity of 0.7 at 25°C. By treating the denaturation of DNA as a phase transition, an equation has immmi derived relating the destabilizing effect of trifluoroacetate to the loss of hydration on denaturation. The hydration of native DNA is abnormally high in the presence of this anion, and the loss of hydration on denaturation is greater than in CsCl. In addition, trifluoroacetate appears to decrease the ΔHof denaturation.     

Sequence dependent hydration of DNA

The transitions between the different helical conformations of DNA depend on the base sequence and the ambient conditions such as humidity and counter-ion concentration. In this study energy minimization techniques have been used to locate water molecule sites around nucleotides especially those which form hydrogen bonds between two or more nucleotide atoms and thus form solvent mediated bridges. We have studied several sequences and find that those which are known not to exist in the low hydration ‘A’ form have very similar number of bridging sites in both ‘A’ and ‘B’ conformations. Those sequences which are found in the ‘A’ conformation have considerably more bridging sites in this low hydration form than in the ‘B’ conformation. Sequence related solvent effects for a given conformation have also been analysed.