Theoretical Study of Ethidium Intercalation in Triple-Stranded DNA and at Triplex-Duplex Junctions

Abstract

The contribution of different factors in the interaction of ethidium intercalated into various sequences of a triple helix, or in the region of the junction between the double- and triple-stranded DNA has been studied by energy minimization. It is found that in the total energy of the ethidium - triple helix complexes, a particular electrostatic contribution emerges due to the presence of protonated cytosines in the triple helix. This parameter is determinant in the sequence-specificity of ethidium binding to the triple helix. The preferred intercalation sites of ethidium in the triple helix are proposed. The interaction of ethidium at the triplex-duplex junction, and its effects are also discussed. This study is aimed at searching for new drugs specific for the triple helix, or for the triplex-duplex junctions.     

Sequence-specific recognition of the major groove of DNA by oligodeoxynucleotides via triple helix formation. Footprinting studies.

Homopyrimidine oligodeoxynucleotides recognize the major groove of the DNA double helix at homopurine.homopyrimidine sequences by forming local triple helices. The oligonucleotide is bound parallel to the homopurine strand of the duplex. This binding can be revealed by a footprinting technique using copper-phenanthroline as a cleaving reagent. Oligonucleotide binding in the major groove prevents cleavage by copper-phenanthroline. The cleavage patterns on opposite strands of the duplex at the boundaries of the triple helix are asymmetric. They are shifted to the 3'-side, indicating that the copper-phenanthroline chelate binds in the minor groove of the duplex structure. Binding of the chelate at the junction between the triple and the double helix is not perturbed on the 5'-side of the bound homopyrimidine oligonucleotide. In contrast, a strong enhancement of cleavage is observed on the purine-containing strand at the triplex-duplex junction on the 3'-side of the homopyrimidine oligonucleotide.     

Intercalation of ethidium bromide into a triple-stranded oligonucleotide.

We have examined the ability of a cationic planar chromophore, ethidium bromide, to intercalate into a short, defined triple helix. Using UV absorption, fluorescence spectroscopy and a gel retardation assay we demonstrate that ethidium bromide is able to bind to a triple helix with a lower affinity than to the corresponding duplex. Energy transfer from base triplets to ethidium shows that ethidium is intercalated into the triple helix. The spectroscopic characteristics of ethidium intercalated into a triplex are similar to those observed for intercalation into duplex DNA.     

Binding of ethidium bromide to a DNA triple helix. Evidence for intercalation

The interaction of ethidium, a DNA intercalator, with the poly(dA).poly(dT) duplex and the poly (dA).2poly(dT) triplex has been investigated by a variety of spectrophotometric and hydrodynamic techniques. The fluorescence of ethidium is increased when either the duplex or triplex form is present. Binding constants, determined from absorbance measurements, indicate that binding to the triple helical form is substantially stronger than to the duplex, with a larger binding site size (2.8 base triplets compared to 2.4 base pairs). Furthermore, while binding to poly(dA).poly(dT) shows strong positive cooperativity, binding to the triplex is noncooperative. Thermal denaturation experiments demonstrate that ethidium stabilizes the triple helix. Binding to either form induces a weak circular dichroism band in the visible wavelength region, while in the region around 310 nm, there is a band that is strongly dependent on the degree of saturation of the duplex, and which is positive for the duplex but negative for the triplex. Both fluorescence energy transfer and quenching studies provide evidence of intercalation of ethidium in both duplex and triplex complexes. Binding of ethidium leads to an initial decrease in viscosity for both the duplex and triplex structures, followed by an increase, which is greater for the duplex. Taken together, these results strongly suggest that ethidium binds to the poly (dA).2poly(dT) triple helix via an intercalative mechanism.     

Oligodeoxynucleotide-directed photo-induced cross-linking of HIV proviral DNA via triple-helix formation.

The HIV proviral genome contains two copies of a 16 bp homopurine.homopyrimidine sequence which overlaps the recognition and cleavage site of the Dra I restriction enzyme. Psoralen was attached to the 16-mer homopyrimidine oligonucleotide, d5'(TTTTCT-TTTCCCCCCT)3', which forms a triple helix with this HIV proviral sequence. Two plasmids, containing part of the HIV proviral DNA, with either one (pLTR) or two (pBT1) copies of the 16-bp homopurine.homopyrimidine sequence and either 4 or 14 Dra I cleavage sites, respectively, were used as substrates for the psoralen-oligonucleotide conjugate. Following UV irradiation the two strands of the DNA targeted sequence were cross-linked at the triplex-duplex junction. The psoralen-oligonucleotide conjugate selectively inhibited Dra I enzymatic cleavage at sites overlapping the two triple helix-forming sequences. A secondary triplex-forming site of 8 contiguous base pairs was observed on the pBT1 plasmid when binding of the 16 base-long oligonucleotide was allowed to take place at high oligonucleotide concentrations. Replacement of a stretch of six cytosines in the 16-mer oligomer by a stretch of six guanines increased binding to the primary sites and abolished binding to the secondary site under physiological conditions. These results demonstrate that oligonucleotides can be designed to selectively recognize and modify specific sequences in HIV proviral DNA.     

Site-specific intercalation at the triplex-duplex junction induces a conformational change which is detectable by hypersensitivity to diethylpyrocarbonate.

Using site-specific intercalation directed by intermolecular triplex formation, the conformation of an intercalation site in DNA was examined by footprinting with the purine-specific (A much greater than G) reagent diethylpyrocarbonate. Site specific intercalation was achieved by covalently linking an intercalator to the 5' end of a homopyrimidine oligodeoxynucleotide, which bound to a homopurinehomopyrimidine stretch in a recombinant plasmid via intermolecular triplex formation. This directs intercalation to a single site in 3kb of DNA at the 5' triplex-duplex junction. Footprinting with diethylpyrocarbonate and dimethylsulphate revealed strong protection from modification of adenine residues within the triple-helix in concordance with their Hoogsteen pairing with the third strand, and a strong hypersensitivity to diethylpyrocarbonate at the first adenine of the duplex. This result indicates that intercalation at this site induces a conformational change at the 5' triplex-duplex junction. Furthermore, the same diethlypyrocarbonate hypersensitivity was observed with an unmodified triple-strand forming oligonucleotide and a range of intercalating molecules present in solution. Thus the 5' triplex-duplex junction is a strong binding site for some intercalating molecules and the junction undergoes a conformational change which is sensitive to diethylpyrocarbonate upon insertion of the planar aromatic chromophore. This conformational change can be used to direct a single-strand cut in duplex DNA to a defined site.     

Spectroscopic studies on the formation and thermal stability of DNA triplexes with a benzoannulated delta-carboline-oligonucleotide conjugate

A benzoannulated delta-carboline with a phenyl substituent has been covalently tethered to the 3'-end of a triplex-forming oligonucleotide and its ability to bind and stabilize DNA triple helices has been examined by various spectroscopic methods. UV thermal melting experiments were conducted with different hairpin duplexes and with a complementary single-stranded oligonucleotide as targets for the conjugate. The delta-carboline ligand preferentially binds triplexes over duplexes and leads to a temperature increase of the triplex-to-duplex transition by up to 23 degrees C. The results obtained from UV, CD and fluorescence measurements suggest that the delta-carboline ligand exhibits specific interactions with a triplex and favors binding by intercalation at the triplex-duplex junction.     

Structural hierarchy controls deformation behavior of collagen

The structure of collagen, the most abundant protein in mammals, consists of a triple helix composed of three helical polypeptide chains. The deformation behavior of collagen is governed by molecular mechanisms that involve the interaction between different helical hierarchies found in collagen. Here, we report results of Steered Molecular Dynamics study of the full-length collagen molecule (～290 nm). The collagen molecule is extended at various pulling rates ranging from 0.00003/ps to 0.012/ps. These simulations reveal a new level of hierarchy exhibited by collagen: helicity of the triple chain. This level of hierarchy is apparent at the 290 nm length and cannot be observed in the 7-9 nm models often described to evaluate collagen mechanics. The deformation mechanisms in collagen are governed by all three levels of hierarchy, helicity of single chain (level-1), helical triple helix (level-2), and hereby described helicity of the triple chain (level-3). The mechanics resulting from the three levels is described by an interlocking gear analogy. In addition, remarkably, the full-length collagen does not show much unwinding of triple helix unlike that exhibited by short collagen models. Further, the full-length collagen does not show significant unwinding of the triple helix, unlike that exhibited by short collagen. Also reported is that the interchain hydrogen bond energy in the full-length collagen is significantly smaller than the overall interchain nonbonded interaction energies, suggesting that the nonbonded interactions have far more important role than hydrogen bonds in the mechanics of collagen. However, hydrogen bonding is essential for the triple helical conformation of the collagen. Hence, although mechanics of collagen is controlled by nonbonded interchain interaction energies, the confirmation of collagen is attributed to the interchain hydrogen bonding.     

Free energy of hydration of collagen models and the enthalpy of the transition between the triple-helical coiled-coil and single-stranded conformations

Interactions with water make an important contribution to the free energy of stabilization of the collagen triple helix, but they do not alter the structure of the triple helix, i.e., the packing geometry of the three strands. Conformational energy computations have been carried out on poly(tripeptide) analogues of collagen, with the introduction of a newly developed form of a hydration shell model to compute the free energy of hydration. The most stable triple helix formed by poly(Gly-Pro-Pro), obtained earlier from conformational energy computations [M. H. Miller & H. A. Scheraga (1976) J. Polym. Sci. Polym. Symp. 54, 171], with a structure that is very closely similar to the observed structure, is strongly favored over other three-strand complexes, both in the absence and the presence of hydration. The hydration shell model does not provide an explanation for the increased stability of the poly(Gly-Pro-Hyp) triple helix as compared to poly(Gly-Pro-Pro). It appears that the difference should be attributed to specific binding of water, and effect that is not yet included in the present version of the hydration shell model. On the other hand, this model accounts for the observed enthalpy of unfolding of a poly(Gly-Pro-Pro) triple helix to isolated single chains in solution in terms of intramolecular noncovalent interactions and the free energy of hydration.     

Solution structural studies on human erythrocyte alpha-spectrin tetramerization site

We have determined the solution NMR structure of a recombinant peptide that consists of the first 156 residues of erythroid alpha-spectrin. The first 20 residues preceding the first helix (helix C') are in a disordered conformation. The subsequent three helices (helices A1, B1, and C1) form a triple helical bundle structural domain that is similar, but not identical, to previously published structures for spectrin from Drosophila and chicken brain. Paramagnetic spin label-induced NMR resonance broadening shows that helix C', the partial domain involved in alpha- and beta-spectrin association, exhibits little interaction with the structural domain. Surprisingly, helix C' is connected to helix A1 of the structural domain by a segment of 7 residues (the junction region) that exhibits a flexible disordered conformation, in contrast to the predicted rigid helical structure. We suggest that the flexibility of this particular junction region may play an important role in modulating the association affinity of alpha- and beta-spectrin at the tetramerization site of different isoforms, such as erythroid spectrin and brain spectrin. These findings may provide insight for explaining various physiological and pathological conditions that are a consequence of varying alpha- and beta-subunit self-association affinities in their formation of the various spectrin tetramers.     

Stabilization of two-stranded ribohomopolymer helices and destabilization of a three-stranded helix by ethidium bromide (Short Communication)

Thermal-denaturation profiles of helical polynucleotides have been measured in the presence of increasing concentrations of ethidium bromide. The poly(A).poly(U) helix is strongly stabilized by binding of ethidium, to much the same extent as is DNA, but the stabilizing effect on poly(I).poly(C) is much smaller. In the poly(A).2poly(U) system the drug selectively destabilizes and eventually destroys the triple helix, leaving only the double-helix-to-coil transition.     

The fluorescence anisotropy decay due to energy transfers occuring in the ethidium bromide-DNA complex. Determination of the deformation angle of the DNA helix

It has been shown in a preceding work that the fluorescence anisotropy decay of ethidium bromide-DNA complex is accelerated by energy migration between dyes bound to the same DNA molecule. In the present work, this result is confirmed. A quantitative analysis has been performed in the following way. The spectroscopic term of the transfer rate constant has been accurately reevaluated by quantum yield and spectral measurements. One assumes that the dye intercalates between two adjacent base pairs and that its distribution is random along the DNA molecule. One introduces the deformation angle δ of the DNA helix induced by the ethidium bromide intercalation. For several values of δ, the energy migration contribution to the anisotropy decay is computed by a Monte Carlo method. In multiplying these computed functions by the measured brownian anisotropy, one obtains the anisotropy decay curve. Comparison with the experimental data leads to the conclusion that the ethidium bromide molecule unwinds the DNA helix by an angle δ = ?16°. This result is m agreement with the work of other authors. We think that the method used here may provide accurate information on the spatial distribution of an array of chromophores bound to a rigid structure.     

A calorimetric study of the triple helix formation: interaction of poly(C) - guanosine - protonated poly(C).]

The triple helix formation of poly(C) - guanosine - poly(C+) was investigated by the help of an LKB scanning micro-calorimeter. The existence of the triple helix could also be shown by recording the melting curves. The ultraviolet absorption at different wave lengths namely 275 nm, 260 nm, and 245 nm was plotted as a function of the temperature. Furthermore formation of the triple helix was shown by plotting the ultraviolet absorption at 245 nm during the increasing addition of guanosine solution to a fixed amount of poly(C) in the solution. Finally the formation of the triple helix was demonstrated by plotting the ultraviolet absorption at 245 nm of a certain mixture of the components while the pH value of the solution was continuously lowered. All these methods show that the monomer interacts with the polymer double helix to form a triple helix. The calorimetric measurements show that the reaction enthalpy is concentration dependent. Above a threshold concentration a rapid increase of the reaction enthalpy is observed. This increase occurs in a very narrow concentration interval. Above this interval a final value of the reaction enthalpy is reached. The amount of the reaction enthalpy for the interaction of guanosine with poly(C) - poly(C+) double helix is 5.5 Kcal (mol base triplet)-1.     

A computational and experimental study of the bending induced at a double-triple helix junction.

We have studied the conformation of a 17 base-pair homopyrimidine.homopurine triple helix formed on a fragment of duplex DNA derived from Simian Virus SV40. Gel retardation assays indicate that an 80 base-pair fragment has an altered conformation when the triple helix is formed, which is most likely to result from an induced bend in the DNA. Investigation of the detailed conformation of the double helix-triple helix junctions has been performed by means of molecular modelling. Bending on the 5' and 3' sides of the third strand oligonucleotide are not located at equivalent positions with respect to the junctions, which is explained in terms of base stacking. The junction effects on DNA structure, induced by the requirement for cytosine protonation in the Hoogsteen-bonded strand to form CGC+ base triplets, are also discussed.     

Deoxymyoglobin studied by the conformational normal mode analysis. II. The conformational change upon oxygenation

The conformational change taking place in myoglobin concomitantly with the observed geometrical change at the heme-His(F8) linkage upon oxygenation is studied by normal mode analysis, which is based on the quadratic approximation of the conformational energy function. The heme-globin interaction energy increases for this change by 8.114 kcal/mol when both the heme group and the globin molecule are held rigid. When they are permitted flexibility, the interaction energy relaxes by 7.038 kcal/mol, and the difference (1.076 kcal/mol) is distributed as strain energy within the molecule. This increase is the work necessary for the heme group to move against the force exerted by the globin. If the force is assumed to be invariable during this move, the work is small, 0.276 kcal/mol, meaning that the force is strongly variable. Furthermore, this means that the heme group is located near the equilibrium point of the potential energy of the heme-globin interaction. The change in the localized strain energy stored in the force field at the linkage between the heme and the imidazole of HisF8 is estimated to be of the same order of magnitude as the distributed energy. The largest atomic displacements are observed at the region from the F helix to the beginning of the G helix, and secondary large displacements occur at several regions, i.e, the A helix, from the C helix to the CD corner, the E helix, and the C-terminal side of the H helix. All of these regions have strong dynamic interactions with the heme group, either directly or indirectly. Their secondary structures show complex deformations. In other parts, relatively rigid segments undergo rotational and/or bending changes in a way consistent with the large changes described above and close atomic packing within the molecule. The calculated conformational change is decomposed to vibrational normal modes of deoxymyoglobin. The magnitude of the conformational change, measured by the mass-weighted mean-square atomic displacement, is accounted for up to 92.0% by the 151 normal modes with frequencies lower than 40 cm-1. In descending order of contribution, the first six modes, each of which has a frequency lower than 12 cm-1, account for up to 57.4%. This means that the functionally important conformational change can well be expressed in terms of a relatively small number of collective low frequency normal modes.     

A directional nucleation-zipping mechanism for triple helix formation

A detailed kinetic study of triple helix formation was performed by surface plasmon resonance. Three systems were investigated involving 15mer pyrimidine oligonucleotides as third strands. Rate constants and activation energies were validated by comparison with thermodynamic values calculated from UV-melting analysis. Replacement of a T·A base pair by a C·G pair at either the 5′ or the 3′ end of the target sequence allowed us to assess mismatch effects and to delineate the mechanism of triple helix formation. Our data show that the association rate constant is governed by the sequence of base triplets on the 5′ side of the triplex (referred to as the 5′ side of the target oligopurine strand) and provides evidence that the reaction pathway for triple helix formation in the pyrimidine motif proceeds from the 5′ end to the 3′ end of the triplex according to the nucleation-zipping model. It seems that this is a general feature for all triple helices formation, probably due to the right-handedness of the DNA double helix that provides a stronger base stacking at the 5′ than at the 3′ duplex–triplex junction. Understanding the mechanism of triple helix formation is not only of fundamental interest, but may also help in designing better triple helix-forming oligonucleotides for gene targeting and control of gene expression.     

Contributions of cation-π interactions to the collagen triple helix stability

Cation–π interactions are found to be an important noncovalent force in proteins. Collagen is a right-handed triple helix composed of three left-handed PPII helices, in which (X–Y-Gly) repeats dominate in the sequence. Molecular modeling indicates that cation–π interactions could be formed between the X and Y positions in adjacent collagen strands. Here, we used a host–guest peptide system: (Pro-Hyp-Gly)₃-(Pro-Y-Gly-X-Hyp-Gly)-(Pro-Hyp-Gly)₃, where X is an aromatic residue and Y is a cationic residue, to study the cation–π interaction in the collagen triple helix. Circular dichroism (CD) measurements and T_m data analysis show that the cation–π interactions involving Arg have a larger contribution to the conformational stability than do those involving Lys, and Trp forms a weaker cation–π interaction with cationic residues than expected as a result of steric effects. The results also show that the formation of cation–π interactions between Arg and Phe depends on their relative positions in the strand. Moreover, the fluorinated and methylated Phe substitutions show that an electron-withdrawing or electron-donating substituent on the aromatic ring can modulate its π–electron density and the cation–π interaction in collagen. Our data demonstrate that the cation–π interaction could play an important role in stabilizing the collagen triple helix.     

Two different heparin-binding domains in the triple-helical domain of ColQ,the collagen tail subunit of synaptic acetylcholinesterase

ColQ, the collagen tail subunit of asymmetric acetylcholinesterase, is responsible for anchoring the enzyme at the vertebrate synaptic basal lamina by interacting with heparan sulfate proteoglycans. To get insights about this function, the interaction of ColQ with heparin was analyzed. For this, heparin affinity chromatography of the complete oligomeric enzyme carrying different mutations in ColQ was performed. Results demonstrate that only the two predicted heparin-binding domains present in the collagen domain of ColQ are responsible for heparin interaction. Despite their similarity in basic charge distribution, each heparin-binding domain had different affinity for heparin. This difference is not solely determined by the number or nature of the basic residues conforming each site, but rather depends critically on local structural features of the triple helix, which can be influenced even by distant regions within ColQ. Thus, ColQ possesses two heparin-binding domains with different properties that may have non-redundant functions. We hypothesize that these binding sites coordinate acetylcholinesterase positioning within the organized architecture of the neuromuscular junction basal lamina.