Interaction of molecules with nucleic acids. II. Two pairs of families of intercalation sites, unwinding angles, and the neighbor-exclusion principle

Intercalation-site geometries are generated for a tetramer duplex extracted from B-DNA. Glycosidic angles and puckers of the deoxyribose sugar groups bonded to base pairs BP₁ and BP₄, namely, those at either end of the tetramer duplex, are assumed to be those of B-DNA to insure continuity. All possible geometrical conformations for combinations of C(2′)-endo, C(3′)-endo, C(2′)-exo, and C(3′)-exo sugar puckers are determined for the tetranucleotide backbone. Those with minimum energy are selected as candidates for intercalation sites. Calculations reveal two pairs of physically meaningful families of intercalation sites which occur in two distinct regions, I and II, of helical angles which orient BP₂ relative to BP₃ and with the helical axis disjointed between these base pairs. For each site I and II within BP₂ and BP₃, there are two distinct backbone conformations, A and B, connecting BP₃ to BP₄ or BP₁ to BP₂ which do not disrupt backbone conformations connecting BP₂ to BP₃. Hence two pairs, IA and IB, and IIA and IIB, of intercalation sites exist in which the sugar puckers along the backbone of the tetramer alternate from C(2′)-endo to C(3′)-endo on the backbone (5′p3′) connecting BP₂ to BP₃. The glycosidic angles of the C(3′)-endo sugar χ₃^γ are, coincidentally, 80° ± 2° for both conformations γ = A and B connecting BP₃ to BP₄ along the phosphate backbone (5′p3′). Consistent with the theoretical results, the experimental unwinding angles can be grouped into two categories with absolute values of 18° and 26°. The theoretical unwinding angles for sites IA and IB of 16° and for sites IIA and IIB of 20° occur for a displacement of -0.8 Å in the helical axes of BP₂ and BP₃ and for a 100% G·C composition, with a decrease depending on the amount of A·T base pairs present. Ratios of theoretical unwinding angles of sites I and II, which range from 0.75 to 0.84 for the two principal sites, compare well with the experimental value of 0.71. The theoretical results, in agreement with experimental observation, provide a new interpretation of the nature and conformation of the possible binding sites. Conformations obtained from these studies of intercalation sites in a tetramer duplex are used to rationalize the well-known neighbor-exclusion principle. The possibility of violation of this principle is demonstrated by the existence of two families of physically meaningful conformations. Conformations of unconstrained dimer duplexes are also obtained, one of which corresponds to the experimental crystal structure of ethidium–dinucleoside complexes, but these cannot be joined to the B-DNA structure. Backbone conformations of the tetramer duplex can be constructed until the base-pair separation reaches 8.25 Å, which may limit the molecules that can intercalate.     

Interactions of molecules with nucleic acids. I. An algorithm to generate nucleic acid structures with an application to the B-DNA structure and a counterclockwise helix.

An algorithm is developed that enables the routine determination of backbone conformations of nucleic acids. All atomic positions including hydrogen are specified in accord with experimental bond lengths and angles but with theoretically determined conformational angles. For two Watson-Crick base pairs at a separation of 3.38 Å, and perpendicular to a common helical axis, minimum energy configurations are found for all 10 combinations at helical angles of α ～ 36°–38°, corresponding to the B-DNA structure with C(2′)-endo sugar puckers. Backbone configurations exist only within the range 35.5° ? α ? 42°, which suggests the origin of the 10-fold helix. Calculated stacking energies for the B-DNA structure increases for each of the clustered groups of base pairs: G·C with G·C, G·C with A·T, and A·T with A·T, and they are in approximate agreement with experimental observations. The counter-clockwise helix is examined, and physically meaningful structures are found only when the helical axes of successive base pairs are disjointed.     

Theoretical studies of nucleic acid interactions. I. Estimates of conformational mobility in intercalated chains

A combined geometric and potential-energy analysis has been carried out to identify the torsional arrangements of the nucleic acid chain that can accommodate the intercalation of small planar moieties. In contrast to previous theoretical efforts, which detail local conformations after adjacent bases are positioned in space, the likely geometries are found here on the basis of the base orientations that result from all feasible combinations of the nine torsional variables of the basic dinucleotide intercalation unit. The relatively mobile nature of the sugar-phosphate backbone, together with the fairly long stretches of chemical bonds between adjacent units, is apparently responsible for the large number of feasible binding geometries. Some previously overlooked conformations with unusual sugar-puckering combinations and various phosphodiester arrangements are found in the survey. A large proportion of the energetically favored intercalation states are closely related to the backbone conformations of familiar double-helical models such as A-, B-, and Z-DNA, as well as the Watson-Crick model. Moreover, the intercalated forms are found to interconvert smoothly along a continuous conformational pathway. The intercalation structures derived from x-ray crystallographic analyses of drug-oligonucleotide complexes, in contrast, are stiff three-dimensional forms essentially frozen in a single domain of conformation space. Specific ligand-nucleic acid interactions that may be responsible for the experimental observations are not included in this study. The classical intramolecular potential energies reported here are highly approximate, providing only rough gauges of the relative importance of the many competing conformations.     

Conformational characteristics and correlations in crystal structures of nucleic acid oligonucleotides: evidence for sub-states

Sugar phosphate backbone conformations are a structural element inextricably involved in a complete understanding of specific recognition nucleic acid ligand interactions, from early stage discrimination of the correct target to complexation per se, including any structural adaptation on binding. The collective results of high resolution DNA, RNA and protein/DNA crystal structures provide an opportunity for an improved and enhanced statistical analysis of standard and unusual sugar-phosphate backbone conformations together with corresponding dinucleotide sequence effects as a basis for further exploration of conformational effects on binding. In this study, we have analyzed the conformations of all relevant crystal structures in the nucleic acids data base, determined the frequency distribution of all possible epsilon, zeta, alpha, beta and gamma backbone angle arrangements within four nucleic acid categories (A-RNA and A-DNA, free and bound B-DNA) and explored the relationships between backbone angles, sugar puckers and selected helical parameters. The trends in the correlations are found to be similar regardless of the nucleic acid category. It is interesting that specific structural effects exhibited by the different unusual backbone sub-states are in some cases contravariant. Certain alpha/gamma changes are accompanied by C3' endo (north) sugars, small twist angles and positive values of base pair roll, and favor a displacement of nucleotide bases towards the minor groove compared to that of canonical B form structures. Unusual epsilon/zeta combinations occur with C2' (south) sugars, high twist angles, negative values of base pair roll, and base displacements towards the major groove. Furthermore, any unusual backbone correlates with a reduced dispersion of equilibrium structural parameters of the whole double helix, as evidenced by the reduced standard deviations of almost all conformational parameters. Finally, a strong sequence effect is displayed in the free oligomers, but reduced somewhat in the ligand bound forms. The most variable steps are GpA and CpA, and, to a lesser extent, their partners TpC and TpG. The results provide a basis for considering if the variable and non-variable steps within a biological active sequence precisely determine morphological structural features as the curvature direction, the groove depth, and the accessibility of base pair for non covalent associations.     

Interactions of molecules with nucleic acids. XII. Theoretical model for the interaction of a fragment of bleomycin with DNA

An intercalation model of a complex between DNA and a bleomycin fragment (BLMF), consisting of the bithiazole core and an amide and a protonated amino substituent, is presented. The model, which shows a preference for BLMF with the protonated amine in the minor groove and the acetyl terminal inserted into either the minor and major grooves, respectively, agrees with recently obtained nmr data. The selection of sites I and II, which have the smallest unwinding of the three theoretical intercalation sites, is consistent with the experimental unwinding angle of 12°. The bithiazole moiety stacks between two base pairs of the double helix, while the protonated substituent interacts ionically with the negatively charged regions of the backbone in the minor groove of the DNA. The protonated amine also forms an intramolecular hydrogen bond with the carbonyl oxygen of the amide group on the same substituent. Analysis of drug complexes with different base-pair sequences reveal four energetically defined groups. The relative energy of the dimer duplex complexes of BLMF correlates with bleomycin's observed base-sequence specificity upon cleavage. The most stable intercalation complexes form adjacent to the bases cleaved most readily. This correlation suggests a primary connection between intercalation and cleavage. A model cleavage site based on these preliminary theoretical calculations and the experimental observations is proposed. It consists of an intercalation site in a trimer duplex. Pyrimidine(p)purine sequences are the predominant sites for intercalation, and the base adjacent to the site at the (3′) end is cleaved.     

Interactions of molecules with nucleic acids. III. Steric and electrostatic energy contours for the principal intercalation sites,prerequisites for binding,and the exclusion of essential metabolites from intercalation

Based on steric and electrostatic considerations, the prerequisites for binding to DNA via the intercalation mechanism are proposed. Steric contour energy curves are presented to demonstrate the region inaccessible to an intercalant. They are calculated with a 6-n (n = 14) potential. This method is a soft potential analog of an excluded-volume approach. Electrostatic contours on the steric surface illustrate the relatively positive and negative regions of the binding site. The principal intercalation sites, predicted to fit into B-DNA via a tetramer-duplex unit, and the unconstrained dimer-duplex units, obtained in crystal structures, are examined. These contours illustrate the requirements of size, conformation, and net atomic charges necessary for intercalation and optimum binding. Based on the limited space available for intercalation by the presence of the backbone and the maximum base-pair separation of 8.25 Å, an Essential Metabolite Exclusion Hypothesis is presented.     

A Monte Carlo method for generating structures of short single-stranded DNA sequences

A Monte Carlo method has been developed for generating the conformations of short single-stranded DNAs from arbitrary starting states. The chain conformers are constructed from energetically favorable arrangements of the constituent mononucleotides. Minimum energy states of individual dinucleotide monophosphate molecules are identified using a torsion angle minimizer. The glycosyl and acyclic backbone torsions of the dimers are allowed to vary, while the sugar rings are held fixed in one of the two preferred puckered forms. A total of 108 conformationally distinct states per dimer are considered in this first stage of minimization. The torsion angles within 5 kcal/mole of the global minimum in the resulting optimized states are then allowed to vary by ±10° in an effort to estimate the breadth of the different local minima. The energies of a total of 2187 (3⁷) angle combinations are examined per local conformational minimum. Finally, the energies of all dinucleotide conformers are scaled so that the populations of differently puckered sugar rings in the theoretical sample match those found in nmr solution studies. This last step is necessitated by limitations in the theoretical methods to predict DNA sugar puckering accurately. The conformer populations of the individual acyclic torsion angles in the composite dimer ensembles are found to be in good agreement with the distributions of backbone conformations deduced from nmr coupling constants and the frequencies of glycosyl conformations in x-ray crystal structures, suggesting that the low energy states are reasonable. The low energy dimer forms (consisting of 150–325 conformational states per dimer step) are next used as variables in a Monte Carlo algorithm, which generates the conformations of single-stranded d(CX_nG) chains, where X = A, T and n = 3, 4, 5. The oligonucleotides are built sequentially from the 5′ end of the chain using random numbers to select the conformations of overlapping dimer units. The simulations are very fast, involving a total of 10⁶ conformations per chain sequence. The potential errors in the buildup procedure are minimized by taking advantage of known rotational interdependences in the sugar–phosphate backbone. The distributions of oligonucleotide conformations are examined in terms of the magnitudes, positions, and orientations of the end-to-end vectors of the chains. The differences in overall flexibility and extension of the oligomers are discussed in terms of the conformations of the constituent dinucleotide steps, while the general methodology is discussed and compared with other nucleic acid model building techniques. © 1993 John Wiley & Sons, Inc.     

Interactions of molecules with nucleic acids. IV. Binding energies and conformations of acridine and phenanthridine compounds in the two principal and in several unconstrained dimer-duplex intercalation sites

The binding positions and relative minimum binding energies are calculated for complexes of 9-aminoacridine, proflavine, N-methylphenanthridinium, and ethidium in theoretically determined intercalation sites in B-DNA (sites I and II) and in unconstrained dimer-duplex sites. The selection of site I in B-DNA by these compounds agrees with the theoretical interpretation of studies of unwinding angles in closed circular DNA in all cases but ethidium, which is predicted to select site II. The most stable binding positions of the acridines and ethidium in unconstrained dimer-duplex units agree with experimental results of intercalation complexes of dinucleoside monophosphate units. Base-pair specificity for Watson-Crick pairing is examined. The energy of an intercalation complex is partitioned into ΔE₂₃, the energy required to open base pairs BP₂ and BP₃ in B-DNA to a site, and ΔE_In, the energy change when a free molecular intercalates. ΔE₂₃ depends strongly on the base-pair sequence, whereas ΔE_In for the four molecules studied does not. The three most stable sequences contain (pyrimidine)p(purine) units, and this provides a rationale for the exclusive formation of crystals of intercalation complexes with these units. In spite of this selectivity, the distribution of G?C and A?T base pairs is equal for these three units and persists as the more unstable sequences are included. Therefore, specificity arises from the interaction between the base pairs and the 2′-deoxyribose 5′-monophosphate backbone for the opening of B-DNA to an intercalation site and not from the interaction between the chromophore and the DNA.     

Intercalation conformations in single- and double-stranded nucleic acids

Two regions in the crystal structure of yeast phenylalanine tRNA, where single-stranded loops interact by intercalation, have been examined in detail. There are four examples of a nucleotide base from one loop intercalating between two sequential bases of another loop in these two regions. These four dinucleoside phosphate conformations serve as models for intercalation in single-stranded nucleic acids. Double-stranded DNA and RNA polymers were constructed by computer model building methods, which incorporated the dinucleoside phosphate conformations found in these single-stranded, intercalation regions in otherwise standard double-helices. The results suggest that it is unlikely that there is a unique intercalation geometry for either single- or double-stranded nucleic acids, but that nucleic acids may assume one of a variety of intercalation geometries which will best accommodate a particular intercalating agent for a particular base sequence.     

The influence of helix morphology on co-operative polyamide backbone conformational flexibility in peptide nucleic acid complexes.

A systematic analysis of peptide nucleic acid (PNA) complexes deposited in the Protein Data Bank has been carried out using a set of contiguous atom torsion angle definitions. The analysis is complemented by molecular mechanics adiabatic potential energy calculations on hybrid PNA-nucleic acid model systems. Hitherto unobserved correlations in the values of the (alpha and epsilon) dihedral angles flanking the backbone secondary amide bond are found. This dihedral coupling forms the basis of a PNA backbone conformation classification scheme. Six conformations are thus characterised in experimental structures. Helix morphology is found to exert a significant influence on backbone conformation and flexibility: Watson-Crick PNA strands in complexes with DNA and RNA, that possess A-like base-pair stacking, adopt backbone conformations distinct from those in PNA.DNA-PNA triplex and PNA-PNA duplex P-helix forms. Solvation effects on Watson-Crick PNA backbone conformation in heterotriplexes are discussed and the possible involvement of inter-conformational transitions and dihedral angle uncoupling in asymmetric heteroduplex base-pair breathing is suggested.     

Mutagen–nucleic acid complexes at the polynucleotide duplex level in solution: Intercalation of proflavine into poly(dA-dT) and the melting transition of the complex

The nmr chemical shifts and line widths of the nucleic acid base and sugar proton resonances and the proflavine ring protons can be monitored through the melting transition of the proflavine + poly(dA-dT) complex, phosphate/dye (P/D) ratio = 24 and 8 in 1M salt solution. The nucleic acid and mutagen protons in the complex are in fast exchange between duplex and strand states with the midpoint of the melting transition monitored at the nucleic acid resonances increasing from 72.6°C for poly(dA-dT) to 78.1°C for the P/D = 24 complex and 83.4°C for the P/D = 8 complex in 1M salt solution. The melting transition monitored by the proflavine resonances were 80.0°C for the P/D = 24 complex and 84.3°C for the P/D = 8 complex in 1M salt solution. Since the nucleic acid is in excess at high P/D ratios, the nucleic acid transitions are an average for the opening of mutagen-free and mutagen-bound base-pair regions, while the proflavine transitions monitor the melting of mutagen-bound base-pair regions. The observed 0.75 to 0.95 ppm unfield shift at all four proflavine protons on formation of the complex with poly(dA-dT) provides direct evidence for intercalation of the mutagen between base pairs of the nucleic acid duplex. We have deduced the approximate overlap geometry between the proflavine ring and nearest-neighbor base pairs at the intercalation site from a comparison between experimental proflavine complexation shifts and those calculated for various stacking orientations. The experimental chemical shift of the poly(dA-dT) adenine H-2 resonance in the duplex state in the absence and presence of proflavine suggests that intercalation occurs preferentially at dT-dA sites. The selective chemical shift changes at the sugar H-2′,2″ and H-3′ resonances of the poly(dA-dT) duplex on complex formation demonstrates changes in the sugar pucker and/or torsion angles of the sugar phosphate backbone at the intercalation site.     

Three-dimensional structures of bulge-containing DNA fragments.

The three-dimensional structure of a DNA tridecamer d(CGCAGAATTCGCG)2 containing bulged adenine bases was determined by single crystal X-ray diffraction methods, at 120 K, to 2.6 A resolution. The structure is a B-DNA type double helix with a single duplex in the asymmetric unit. One of the bulged adenine bases loops out from the double helix, while the other stacks in to it. This is in contrast to our preliminary finding, which indicated that both adenine bases were looped out. This revised model was confirmed by the use of a covalently bound heavy-atom derivative. The conformation of the looped-out bulge hardly disrupts base stacking interactions of the bases flanking it. This is achieved by the backbone making a "loop-the-loop" curve with the extra adenine flipping over with respect to the other nucleotides in the strand. The looped-out base intercalates into the stacked-in bulge site of a symmetrically related duplex. The looped-out and stacked-in bases form an A.A reversed Hoogsteen base-pair that stacks between the surrounding base-pairs, thus stabilizing both bulges. The double helix is frayed at one end with the two "melted" bases participating in intermolecular interactions. A related structure, of the same tridecamer, after soaking the crystals with proflavin, was determined to 3.2 A resolution. The main features of this B-DNA duplex are basically similar to the native tridecamer but differ in detail especially in the conformation of the bulged-out base. Accommodation of a large perturbation such as that described here with minimal disruption of the double helix shows both the flexibility and resiliency of the DNA molecule.     

Interactions of Molecules with Nucleic Acids. X. Covalent Intercalat e Binding of the Carcinogenic BPDE I(+) to Kinked DNA

Abstract

A theoretical model is proposed for the covalent binding of (+) 7 β,8α-dihydroxy-9α, 10α- epoxy-7,8,9,10-tetrahydrobenzo[a]pyrene denoted by BPDE I(+), to N2 on guanine. The DNA must kink a minimum of 39° to allow proper hybrid configurations about the C10 and N2 atoms involved in bond formation and to allow stacking of the pyrene moiety with the non-bonded adjacent base pair. Conservative (same sugar puckers and glycosidic angles as in B-DNA) and non-conservative (alternating sugar puckers as in intercalation sites) conformations are found and they are proposed structures in pathways connecting B-DNA, an intercalation site, and a kink site in the formation of a covalently intercalative bound adduct of BPDE I(+) to N2 on guanine. Stereographic projections are presented for (3′) and (5′) binding in the DNA. Experimental data for bending of DNA by BPDE, orientation of BPDE in DNA and unwinding of superhelical DNA is explained. The structure of a covalent intercalative complex is predicted to result from the reaction. Also, an anti ? syn transition of guanine results in a structure which allows the DNA to resume its overall B-form. The only change is that guanine has been rotated by 200° about its glycosidic bond so that the BPDE I(+) is bound in the major groove. The latter step may allow the DNA to be stored with an adduct which may produce an error in the genetic code.     

DNA fine structure and dynamics in crystals and in solution: the impact of BI/BII backbone conformations

Sugar-phosphate backbone conformations are an important structural element for a complete understanding of specific recognition in nucleic acid-protein interactions. They can be involved both in early stages of target discrimination and in structural adaptation upon binding. In the first part of this study, we have analyzed high-resolution structures of double-stranded B-DNA either isolated or bound to proteins, and explored the impact of both the standard BI and the unusual BII phosphate backbone conformations on neighboring sugar puckers and on selected helical parameters. Correlations are found to be similar for free and bound DNA, and in both categories, the possible facing backbone conformations (BI.BI, BI.BII, and BII.BII) define well-characterized substates in the B-DNA conformational space. Notably, BII.BII steps are characterized by specific, and sequence-independent, structural effects involving reduced standard deviations for almost all conformational parameters. In the second part of this work, we analyze four 10 ns molecular dynamics simulations in explicit solvent on the DNA targets of NF-kappaB and bovine papillomavirus E2 proteins, highlighting the multiplicity of backbone dynamical behavior. These results show sequence effects on the percentages of BI and BII conformers, the preferential state of facing backbones, the occurrence of coupled transitions. The backbone states can consequently be seen as a mechanism for transmitting information from the bases to the phosphate groups and thus for modulating the overall structural properties of the target DNA.     

1 A crystal structures of B-DNA reveal sequence-specific binding and groove-specific bending of DNA by magnesium and calcium

The 1 A resolution X-ray crystal structures of Mg(2+) and Ca(2+) salts of the B-DNA decamers CCAACGTTGG and CCAGCGCTGG reveal sequence-specific binding of Mg(2+) and Ca(2+) to the major and minor grooves of DNA, as well as non-specific binding to backbone phosphate oxygen atoms. Minor groove binding involves H-bond interactions between cross-strand DNA base atoms of adjacent base-pairs and the cations' water ligands. In the major groove the cations' water ligands can interact through H-bonds with O and N atoms from either one base or adjacent bases, and in addition the softer Ca(2+) can form polar covalent bonds bridging adjacent N7 and O6 atoms at GG bases. For reasons outlined earlier, localized monovalent cations are neither expected nor found.Ultra-high atomic resolution gives an unprecedented view of hydration in both grooves of DNA, permits an analysis of individual anisotropic displacement parameters, and reveals up to 22 divalent cations per DNA duplex. Each DNA helix is quite anisotropic, and alternate conformations, with motion in the direction of opening and closing the minor groove, are observed for the sugar-phosphate backbone. Taking into consideration the variability of experimental parameters and crystal packing environments among these four helices, and 24 other Mg(2+) and Ca(2+) bound B-DNA structures, we conclude that sequence-specific and strand-specific binding of Mg(2+) and Ca(2+) to the major groove causes DNA bending by base-roll compression towards the major groove, while sequence-specific binding of Mg(2+) and Ca(2+) in the minor groove has a negligible effect on helix curvature. The minor groove opens and closes to accommodate Mg(2+) and Ca(2+) without the necessity for significant bending of the overall helix.The program Shelxdna was written to facilitate refinement and analysis of X-ray crystal structures by Shelxl-97 and to plot and analyze one or more Curves and Freehelix output files.     

HU protein employs similar mechanisms of minor-groove recognition in binding to different B-DNA sites: demonstration by Raman spectroscopy

The sequence isomers d(CGCAAATTTGCG) and d(TCAAGGCCTTGA) form self-complementary duplexes that present distinct targets for binding of the homodimeric architectural protein HU of Bacillus stearothermophilus (HUBst). Raman spectroscopy shows that although each duplex structure is of the B-DNA type, there are subtle conformational dissimilarities between them, involving torsion angles of the phosphodiester backbone and the arrangements of stacked bases. Each DNA duplex forms a stable stoichiometric (1:1) complex with HUBst, in which the structure of the HUBst dimer is largely conserved. However, the Raman signature of each DNA duplex is perturbed significantly and similarly with HUBst binding, as reflected in marker bands assigned to localized vibrations of the phosphodiester moieties and base residues. The spectral perturbations identify a reorganization of the DNA backbone and partial unstacking of bases with HUBst binding, which is consistent with non-sequence-specific minor-groove recognition. Prominent among the HUBst-induced perturbations of B-DNA are a conversion of approximately one-third of the alpha/beta/gamma torsions from the canonical g(-)/t/g(+) conformation to an alternative conformation, an equivalent conversion of deoxyadenosyl moieties from the C2'-endo/anti to the C3'-endo/anti conformation, and appreciable unstacking of purines. The results imply that each solution complex is characterized by structural perturbations extending throughout the 12-bp sequence. Comparison with previously studied protein/DNA complexes suggests that binding of HUBst bends DNA by approximately 70 degrees.     

Characterization of the 6-methyl isoxanthopterin (6-MI) base analog dimer, a spectroscopic probe for monitoring guanine base conformations at specific sites in nucleic acids

We here characterize local conformations of site-specifically placed pairs of guanine (G) residues in RNA and DNA, using 6-methyl isoxanthopterin (6-MI) as a conformational probe. 6-MI is a base analog of G and spectroscopic signals obtained from pairs of adjacent 6-MI residues reflect base-base interactions that are sensitive to the sequence context, local DNA conformation and solvent environment of the probe bases. CD signals show strong exciton coupling between stacked 6-MI bases in double-stranded (ds) DNA; this coupling is reduced in single-stranded (ss) DNA sequences. Solvent interactions reduce the fluorescence of the dimer probe more efficiently in ssDNA than dsDNA, while self-quenching between 6-MI bases is enhanced in dsDNA. 6-MI dimer probes closely resemble adjacent GG residues, in that these probes have minimal effects on the stability of dsDNA and on interactions with solvent additive betaine. They also serve as effective template bases, although further polymerase-dependent extension of DNA primers past 6-MI template bases is significantly inhibited. These probes are also used to monitor DNA 'breathing' at model replication forks. The 6-MI dimer probe can serve in many contexts as a useful tool to investigate GG conformations at specific sites within the nucleic acid frameworks of functioning macromolecular machines in solution.     

The role of molecular structure of sugar‐phosphate backbone and nucleic acid bases in the formation of single‐stranded and double‐stranded DNA structures

Our previous DFT computations of deoxydinucleoside monophosphate complexes with Na⁺‐ions (dDMPs) have demonstrated that the main characteristics of Watson‐Crick (WC) right‐handed duplex families are predefined in the local energy minima of dDMPs. In this work, we study the mechanisms of contribution of chemically monotonous sugar‐phosphate backbone and the bases into the double helix irregularity. Geometry optimization of sugar‐phosphate backbone produces energy minima matching the WC DNA conformations. Studying the conformational variability of dDMPs in response to sequence permutation, we found that simple replacement of bases in the previously fully optimized dDMPs, e.g. by constructing Pyr‐Pur from Pur‐Pyr, and Pur‐Pyr from Pyr‐Pur sequences, while retaining the backbone geometry, automatically produces the mutual base position characteristic of the target sequence. Based on that, we infer that the directionality and the preferable regions of the sugar‐phosphate torsions, combined with the difference of purines from pyrimidines in ring shape, determines the sequence dependence of the structure of WC DNA. No such sequence dependence exists in dDMPs corresponding to other DNA conformations (e.g., Z‐family and Hoogsteen duplexes). Unlike other duplexes, WC helix is unique by its ability to match the local energy minima of the free single strand to the preferable conformations of the duplex. © 2013 Wiley Periodicals, Inc. Biopolymers 101: 640–650, 2014.     

Interactions of molecules with nucleic acids. VII. Intercalation and T.A specificity of daunomycin in DNA

Intercalation complexes of daunomicin(+1) with tetramer duplexes in DNA are studied with the theoretically determined intercalation sites (I, ?0.4), (II, ?0.4), and (III, ?1.4). These sites occur with base pairs separated by 6.76 Å for helical angles of 26°, 22°, and 8° about the intercalation site. Site I is preferred, and this is in agreement with experimental unwinding angles. Optimum binding positions and conformations are established, and these are in agreement with experimental results from crystal structures. A systematic procedure is devised to study base-pair and base-sequence specificity, which results in the demonstration that the most stable sequences are mainly ↑BP₁, T·A, DAUN, A·T, BP₄↓ and ↑BP₁, T·A, DAUN, G·C, BP₄↓, i.e., with the TpA and CpG (pyrimidine)p(purine) sequences about the intercalation site. These 32 possible sequences are found among the 40 most stable complexes. These theoretical calculations of intercalation complexes with daunomicin(+1) provide the first example in which a drug specifically selects the base pair T·A and prefers it in a particular sequence about the intercalation site. This specificity is in agreement with some experimental results. Problems associated with the interpretation of specificity are discussed in terms of the base, base-pair, and base-sequence resulting from the DNA site and the DNA–drug interactions. T·A specificity is rationalized by noting that the 2′deoxyribo-5′-monophosphate backbone attached to A is slightly more negative than that on the other nucleotides. Hence, a preference exists for binding to the protonated daunosamine (+1) groups. Stereographic projections of daunomycinone and daunomycin(+1) in a bond model and in a space-filling model with steric contours illustrate the results.