Electrostatic free energy landscapes for nucleic acid helix assembly

Metal ions are crucial for nucleic acid folding. From the free energy landscapes, we investigate the detailed mechanism for ion-induced collapse for a paradigm system: loop-tethered short DNA helices. We find that Na⁺ and Mg²⁺ play distinctive roles in helix–helix assembly. High [Na⁺] (>0.3 M) causes a reduced helix–helix electrostatic repulsion and a subsequent disordered packing of helices. In contrast, Mg²⁺ of concentration >1 mM is predicted to induce helix–helix attraction and results in a more compact and ordered helix–helix packing. Mg²⁺ is much more efficient in causing nucleic acid compaction. In addition, the free energy landscape shows that the tethering loops between the helices also play a significant role. A flexible loop, such as a neutral loop or a polynucleotide loop in high salt concentration, enhances the close approach of the helices in order to gain the loop entropy. On the other hand, a rigid loop, such as a polynucleotide loop in low salt concentration, tends to de-compact the helices. Therefore, a polynucleotide loop significantly enhances the sharpness of the ion-induced compaction transition. Moreover, we find that a larger number of helices in the system or a smaller radius of the divalent ions can cause a more abrupt compaction transition and a more compact state at high ion concentration, and the ion size effect becomes more pronounced as the number of helices is increased.     

Pyrimidine morpholino oligonucleotides form a stable triple helix in the absence of magnesium ions

Oligonucleotides can be used as sequence-specific DNA ligands by forming a local triple helix. In order to form more stable triple-helical structures or prevent their degradation in cells, oligonucleotide analogues that are modified at either the backbone or base level are routinely used. Morpholino oligonucleotides appeared recently as a promising modification for antisense applications. We report here a study that indicates the possibility of a triple helix formation with a morpholino pyrimidine TFO and its comparison with a phosphodiester and a phosphoramidate oligonucleotide. At a neutral pH and in the presence of a high magnesium ion concentration (10 mM), the phosphoramidate oligomer forms the most stable triple helix, whereas in the absence of magnesium ion but at a physiological monovalent cation concentration (0.14 M) only morpholino oligonucleotides form a stable triplex. To our knowledge, this is the first report of a stable triple helix in the pyrimidine motif formed by a noncharged oligonucleotide third strand (the morpholino oligonucleotide) and a DNA duplex. We show here that the structure formed with the morpholino oligomer is a bona fide triple helix and it is destabilized by high concentrations of potassium ions or divalent cations (Mg(2+)).     

On the application of polyelectrolyte limiting laws to the helix-coil transition of DNA. III. TDependence of helix stability on excess univalent salt and on polynucleotide phosphate concentration for variable equivalent ratios of divalent metal ion to phosphate.

Using the free energy difference between double-helix and random-coil forms of DNA as a measure of the stability of the double helix, we calculate the dependence of the stability on excess univalent cation concentration and on polynucleotide phosphate concentration, both as functions of the equivalent ratio r of divalent cation-to-phosphate concentrations. The theoretical tool is merely to compare the free energy of one polyelectrolyte solution, characterized by the polyelectrolyte linear charge density, with the free energy of another, characterized by a different value of the charge density. It is assumed only that the charge density of the double helix is greater than that of the coil form. The calculation represents the only molecular theory given to date (for r ≠ O) for these aspects of helix stability. We find that, as excess univalent cation concentration increases, the helix stability increases if r is small but decreases if r is large (i.e., of the order of unity). Moreover, as the concentration of nucleotide phosphate increases, the helix stability does not change for small values of r but increases for large values. For both effects, a continuous transition as a function of r bridges the low-r and high-r behaviour.     

Raman microspectroscopic study of effects of Na(I) and Mg(II) ions on low pH induced DNA structural changes

In this work a confocal Raman microspectrometer is used to investigate the influence of Na(+) and Mg(2+) ions on the DNA structural changes induced by low pH. Measurements are carried out on calf thymus DNA at neutral pH (7) and pH 3 in the presence of low and high concentrations of Na(+) and Mg(2+) ions, respectively. It is found that low concentrations of Na(+) ions do not protect DNA against binding of H(+). High concentrations of monovalent ions can prevent protonation of the DNA double helix. Our Raman spectra show that low concentrations of Mg(2+) ions partly protect DNA against protonation of cytosine (line at 1262 cm(-1)) but do not protect adenine and guanine N(7) against binding of H(+) (characteristic lines at 1304 and 1488 cm(-1), respectively). High concentrations of Mg(2+) can prevent protonation of cytosine and protonation of adenine (disruption of AT pairs). By analyzing the line at 1488 cm(-1), which obtains most of its intensity from a guanine vibration, high magnesium salt protect the N(7) of guanine against protonation. A high salt concentration can prevent protonation of guanine, cytosine, and adenine in DNA. Higher salt concentrations cause less DNA protonation than lower salt concentrations. Magnesium ions are found to be more effective in protecting DNA against binding of H(+) as compared with calcium ions presented in a previous study. Divalent metal cations (Mg(2+), Ca(2+)) are more effective in protecting DNA against protonation than monovalent ions (Na(+)).     

P-Postlabeling high-performance liquid chromatographic improvements to characterize DNA adduct stereoisomers from benzo[a]pyrene and benzo[c]phenanthrene, and to separate DNA adducts from 7,12-dimethylbenz[a]anthracene

Single compounds can generate complex DNA adduct patterns by reactions through different pathways, with different target nucleotides and through different configurations of the products. DNA adduct analysis by ³²P-HPLC was improved by adding an isocratic plateau in an otherwise linear gradient, thereby enhancing resolution of predictable retention time intervals. This enhanced ³²P-HPLC technique was used to analyze and at least partly resolve 14 out of 16 available benzo[c]phenanthrene deoxyadenosine and deoxyguanosine adduct standards, 8 out of 8 available benzo[a]pyrene deoxyadenosine and deoxyguanosine adduct standards, and 51 peaks from 7,12-dimethylbenz[a]anthracene-calf thymus DNA reaction products. The same type of gradient modifications could be used to enhance resolution in analyses of other complex DNA adduct mixtures, e.g., in vivo in humans.     

Nucleic acid helix stability: effects of salt concentration, cation valence and size, and chain length

Metal ions play crucial roles in thermal stability and folding kinetics of nucleic acids. For ions (especially multivalent ions) in the close vicinity of nucleic acid surface, interion correlations and ion-binding mode fluctuations may be important. Poisson-Boltzmann theory ignores these effects whereas the recently developed tightly bound ion (TBI) theory explicitly accounts for these effects. Extensive experimental data demonstrate that the TBI theory gives improved predictions for multivalent ions (e.g., Mg2+) than the Poisson-Boltzmann theory. In this study, we use the TBI theory to investigate how the metal ions affect the folding stability of B-DNA helices. We quantitatively evaluate the effects of ion concentration, ion size and valence, and helix length on the helix stability. Moreover, we derive practically useful analytical formulas for the thermodynamic parameters as functions of finite helix length, ion type, and ion concentration. We find that the helix stability is additive for high ion concentration and long helix and nonadditive for low ion concentration and short helix. All these results are tested against and supported by extensive experimental data.     

Abolition of Intrinsically Bent DNA Structure Components in AT Clusters by Netropsin Interaction; Titration Viscometric Investigations

Abstract

It is argued that the enhancement of the apparent DNA contour length by the specifically binding non-intercalating drug netropsin (Nt) (Reinert et al., NAR 9,2335, 1981) at very low Nt/DNA-phosphate ratios essentially is the result of an abolition of periodically arranged intrinsic helix bends in A · T rich tracts of base pairs.

In the preceding paper the existence of pronounced DNA tertiary structure components has been postulated for (two species of) natural eukaryotic DNA. The resulting model suggests local apparent solenoid-related DNA tertiary structure components at high sodium ion concentration c_s, partly/totally molten out at 45/60 C. With decreasing c_s the tertiary structure components have been found to be gradually reduced, at least below c_s = 0.010 M, as titration viscometrically revealed by a gradual rise of the apparent DNA contour length (Reinert et al., JBSD 9, 537, 1991).

Hence, we performed titration viscometric analyses about Nt interaction with calf thymus DNA (ctDNA) at c_s = 0.075 M, 0.010 M and 0.004 M Na⁺. The concomitant DNA conformational changes are quantitatively described in terms of the relative changes of both DNA persistence length and hydrodynamically operative apparent DNA contour length for the three first resolved interaction modes below a Nt/DNA-P ratio of 0.03.

These experiments, together with previous respective analyses at c_s = 0.20 M Na⁺ and different temperatures (I.e.), suggest that those DNA sites binding Nt most strongly predominantly are responsible for the formation of solenoid-related DNA tertiary structure components. Most probably these are A tract-containing sequences. As the essential factor for their apparent elongation effect at low Na⁺ concentrations, a gradual alteration of the number of base pairs per helix turn seems to occur below c_s = 0.010 M Na⁺ and, concomitantly, a change in phasing between intrinsic helix bends and helix screw.     

Microscopic Mechanics of Hairpin DNA Translocation through Synthetic Nanopores

Nanoscale pores have proved useful as a means to assay DNA and are actively being developed as the basis of genome sequencing methods. Hairpin DNA (hpDNA), having both double-helical and overhanging coil portions, can be trapped in a nanopore, giving ample time to execute a sequence measurement. In this article, we provide a detailed account of hpDNA interaction with a synthetic nanopore obtained through extensive all-atom molecular dynamics simulations. For synthetic pores with minimum diameters from 1.3 to 2.2 nm, we find that hpDNA can translocate by three modes: unzipping of the double helix and—in two distinct orientations—stretching/distortion of the double helix. Furthermore, each of these modes can be selected by an appropriate choice of the pore size and voltage applied transverse to the membrane. We demonstrate that the presence of hpDNA can dramatically alter the distribution of ions within the pore, substantially affecting the ionic current through it. In experiments and simulations, the ionic current relative to that in the absence of DNA can drop below 10% and rise beyond 200%. Simulations associate the former with the double helix occupying the constriction and the latter with accumulation of DNA that has passed through the constriction.     

X-Ray conformational study of the DNA duplex in solution

Earlier x-ray studies on dissolved linear DNA molecules were interpreted on the assumption that the molecules scattered as rigid rods. Improvement in equipment and advances in theory of the scattering from randomly oriented helices prompted us into a reinvestigation of this problem. Careful measurements were made on the scattering from both linear calf thymus DNA and from circular plasmid C0P608 superhelical DNA. Contrary to the earlier work, we find that the scattering patterns show a helical character, with maxima corresponding to those of a helix with pitch angle of 62°, close to that of the C-W helix. The patterns for both types of DNA, although similar, show a 5% displacement of the maximum in the superhelical form, just that expected when the C-W helix is superimposed on a superhelix axis. Introduction of intercalators (PtTS) causes a progressive extension of the helix, as shown by a shift to larger angles and a fading out of the maximum. In the concentration range of 40 mg/mL, interfernce peaks develop, the result of an apparent stacking of the chains, with an interchain distance of ～35 Å.     

Transition Metal Ligands as Novel DNA-Base Substitutes

Abstract

The base modified nucleoside dBP, carrying a non-hydrogen-bonding non-shape complementary base was incorporated into oligonucleotides (Brotschi, C.; Häberli, A.; Leumann C.J. Angew. Chem. Int. Ed. 2001, 40, 3012–3014). This base was designed to coordinate transition metal ions into well defined positions within a DNA double helix. Melting experiments revealed that the stability of a dBP: dBP base couple in a DNA duplex is similar to a dG: dC base pair even in the absence of transition metal ions. In the presence of transition metal ions, melting experiments revealed a decrease in duplex stability which is on a similar order for all metal ions (Mn²⁺, Cu²⁺, Zn²⁺, Ni²⁺) tested.     

Nanoscale time-lapse AFM imaging in solution for DNA aggregation

Revealing the behavior of biofunctional molecules (i.e., nucleic acids, nucleic acid binding reagents, enzymatic proteins, etc.) by monitoring them in solution is important for understanding the nanoscale dynamism of their interactions. Atomic force microscope (AFM) imaging with a dynamic force mode (DFM, i.e., tapping mode) in aqueous solution, has many advantages for the imaging of DNA morphological change at a single molecule scale. Hoechst 33258 (H33258) induces DNA condensation in the presence of its excess concentration. To have a better understanding of the condensation process of DNA with excess H33258, we tried to find the optimum conditions for carrying out time-lapse AFM imaging in aqueous solution. To immobilize DNA on the substrate surface, the mica was modified with the various concentrations of 3-aminopropyltriethoxysilane (APTES) solution. We observed that DNA was minimally immobilized on 0.002% APTES-modified mica surface. Then, we determined that the movement of DNA on the mica surface could be observed in the presence of 500 mM NaCl in 10 mM PBS (pH 7.0). Moreover, after the injection of 5 μM H33258, the partial condensation of DNA was observed.     

Premelting effects in DNA under saltfree conditions

We have studied the electrical conductivity of NaDNA solutions under “saltfree” conditions at temperatures well below the melting point of DNA, using radio-frequency dielectric and noise measurements. A conductivity discontinuity is observed at a temperature well below that at which the usual denaturation processes and trans conformation may commence. The radio-frequency permittivity also exhibits a discontinuity at the same temperature. For the premelting phase, the conductivity versus temperature curves consist of two linear regions with a change in slope occurring at 23°C. This effect is related to the behavior of the ionic sheath covering the DNA macromolecule. The activation energy of the alternative current conductivity as well as that the equivalent noise conductivity results as 3.11 kcal/mole below and 4.08 kcal/mole.     

Different families of double-stranded conformations of DNA as revealed by computer calculations

An algorithm has been developed that permits one to find all possible conformations of the sugar-phosphate backbone for any given disposition of DNA base pairs. For each of the conformations thus obtained, the energy of the helix was calculated by the method of atom-atom potentials. Several isolated regions in the space of the bases′ parameters (Arnott's parameters) have been found for energetically favorable helical structures. Two parameters, the distance of a base pair from the helix axis, D, and the windling angle, τ, allow one to subdivide possible conformations into the families of closely related forms. Two regions (ravines) on the (D, τ) map correspond to the know A and B families. In the B family a continuous transition has been obtained in which the double helix undergoes increasing winding, while the base pairs are moving toward the major (nonglycosidic) groove. Interrelationships between the variables, characterizing the spatial structure of the double helix, D, τ, TL and χ, when going along the bottom of the B ravine, were also obtained. Besides the Known A and B families, several new ones were found to be energetically possible. Among these the strongly underwound helices with the negative D values, as well as the forms with the C4-C5 angle in a trans position, should be mentioned. Biological roles of the different double-stranded conformations, in particular, in protein-nuclei acid interaction are discussed.     

Influence of magnesium ions on helix-coil transition of DNA determined by modified differential scanning calorimeter

The thermal effect of magnesium ions on the helix–coil transition of DNA was studied calorimetrically by a modified differential scanning calorimeter (DSC). It was found that the transition temperature of DNA depends on both the DNA and magnesium ion concentrations. The dependence of the helix–coil transition of DNA on the mole ratio of magnesium ions to DNA(P) can be classified into two groups. When this mole ratio is less than 1, magnesium ions tend to stabilize the double-helix DNA, so that the transition temperature increases linearly and the heat of transition increases significantly with increasing mole ratio. When the mole ratio is more than 1, magnesium ions tend to destabilize the double-helix DNA, so that DNA precipitates when the temperature is raised above the transition temperature. In this case, both the transition temperature and the heat of transition decrease with increasing mole ratio.     

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

Abstract

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.     

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.     

Circular dichroism of polynucleotides: dimers as a function of conformation

Working within the restrictions of a model, we have calculated the circular dichroism of the dinucleoside phosphates ApA, CpC, and CpA for various conformations. Comparing the calculated curves with those measured in aqueous solution we find agreement for (1) ApA as a right-handed helix with both bases either as in B-form DNA, or else rotated 180° around the glycosidic bond, (2) CpC as the right-handed conformation with both bases as in DNA, (3) ApC as either the right-handed conformation with both bases as in DNA, or else as a left-handed helix with both bases rotated 180°, and (4) CpA as either a left-handed helix with both bases in a left-handed DNA, or else in the right-handed conformation with both bases rotated 180°. In addition, we have investigated circular dichroism as a measure of unstacking. We find that opening the bases to a 90° total angle (base planes perpendicular) reduces the intensity of the calculated bands to 20% of their original value. Further, we find that allowing the sliding of one base past the other does not lead to a temperature dependence consistent with experiment.     

Melting of Cross-Linked DNA V. Cross-Linking Effect Caused by Local Stabilization of the Double Helix

Abstract

DNA interstrand cross-links are usually formed due to bidentate covalent or coordination binding of a cross-linking agent to nucleotides of different strands. However interstrand linkages can be also caused by any type of chemical modification that gives rise to a strong local stabilization of the double helix. These stabilized sites conserve their helical structure and prevent local and total strand separation at temperatures above the melting of ordinary AT and GC base pairs. This local stabilization makes DNA melting fully reversible and independent of strand concentration like ordinary covalent interstrand cross-links. The stabilization can be caused by all the types of chemical modifications (interstrand cross-links, intrastrand cross-links or monofunctional adducts) if they give rise to a strong enough local stabilization of the double helix. Our calculation demonstrates that an increase in stability by 25 to 30 kcal in the free energy of a single base pair of the double helix is sufficient for this “cross-linking effect” (i.e. conserving the helicity of this base pair and preventing strand separation after melting of ordinary base pairs). For the situation where there is more then one stabilized site in a DNA duplex (e.g., 1 stabilized site per 1000 bp), a lower stabilization per site is sufficient for the “cross-linking effect” (18–20 kcal). A substantial increase in DNA stability was found in various experimental studies for some metal-based anti-tumor compounds. These compounds may give rise to the effect described above. If ligand induced stabilization is distributed among several neighboring base pairs, a much lower minimum increase per stabilized base pair is sufficient to produce the cross-linking effect (1 bp- 24.4 kcal; 5 bp- 5.3 kcal; 10 bp- 2.9 kcal, 25 bp- 1.4 kcal; 50 bp- 1.0 kcal). The relatively weak non-covalent binding of histones or protamines that cover long regions of DNA (20–40 bp) can also cause this effect if the salt concentration of the solution is sufficiently low to cause strong local stabilization of the double helix. Stretches of GC pairs more than 25 bp in length inserted into poly(AT) DNA also exhibit properties of stabilizing interstrand cross-links.     

Genetic analysis of the LexA repressor: Isolation and characterization of LexA(Def) mutant proteins

Summary We report the isolation of LexA mutant proteins with impaired repressor function. These mutant proteins were obtained by transforming a LexA-deficient recA-lacZ indicator strain with a randomly mutagenized plasmid harbouring the lexA gene and subsequent selection on MacConkey-lactose indicator plates. A total of 24 different lexA(Def) missense mutations were identified. All except three mutant proteins are produced in near-normal amounts suggesting that they are fairly resistant to intracellular proteases. All lexA(Def) missense mutations are situated within the first 67 amino acids of the amino-terminal DNA binding domain. The properties of an intragenic deletion mutant suggest that the part of the amino-terminal domain important for DNA recognition or domain folding should extent at least to amino acids 69 or 70. A recent 2D-NMR study (Lamerichs et al. 1989) has identified three a helices in the DNA binding domain of LexA. The relative orientation of two of them (helices 2 and 3) is reminiscent of, but not identical to, the canonical helix-turn-helix motif suggesting nevertheless that helix 3 might be involved in DNA recognition. The distribution of the lexA(Def) missense mutations along the first 67 amino-terminal amino acids indeed shows some clustering within helix 3, since 8 out of the 24 different missense mutations are found in this helix. However one mutation in front of helix 1 and five mutations between amino acids 61 and 67 suggest that elements other than helices 2 and 3 may be important for DNA binding.