Analysis of local helix geometry in three B-DNA decamers and eight dodecamers

Local variations in B-DNA helix structure are compared among three decamers and eight dodecamers, which contain examples of all ten base-pair step types. All pairwise combinations of helix parameters are compared by linear regression analysis, in a search for internal relationships as well as correlations with base sequence. The primary conclusions are: (1) Three-center hydrogen bonds between base-pairs occur frequently in the major groove at C-C, C-A, A-A and A-C steps, but are less convincing at C-C and C-T steps in the minor groove. The requirements for large base-pair propeller are (1) that the base-pair should be A.T rather than G.C, and (2) that it be involved in a major groove three-center hydrogen bond with the following base-pair. Either condition alone is insufficient. Hence, a large propeller is expected at the leading base-pair of A-A and A-C steps, but not at A-T, T-A, C-A or C-C steps. (2) A systematic and quantitative linkage exists between helix variables twist, rise, cup and roll, of such strength that the rise between base-pairs can hardly be described as an independent variable at all. Two typical patterns of behavior are observed at steps from one base-pair to the next: high twist profile (HTP), characterized by high twist, low rise, positive cup and negative roll, and low twist profile (LTP), marked by low twist, high rise; negative cup and positive roll. Examples of HTP are steps G-C, G-A and Y-C-A-R, where Y is pyrimidine and R is purine. Examples of LTP steps are C-G, G-G, A-G and C-A steps other than Y-C-A-R. (3) The minor groove is especially narrow across the two base-pairs of the following steps: A-T, T-A, A-A and G-A. (4) In general, base step geometry cannot be correlated solely with the bases that define the step in question; the two flanking steps also must be taken into account. Hence, local helix structure must be studied in the context, not of two base-pairs: A-B, but of four: x-A-B-y. Calladine's rules, although too simple in detail, were correct in defining the length of sequence over which a given perturbation is expressed. Whereas ten different two-base steps are possible, allowing for the identity of complementary sequences, there are 136 different four-base steps. Only 33 of these 136 four-base steps are represented in the decamer and dodecamer structures solved to date, and hence it is premature to try to set up detailed structural algorithms. (5) The sugar-phosphate backbone chains of B-DNA place strong limits on sequence-induced structural variation, damping down most variables within four or five base-pairs, and preventing purine-purine anti-anti mismatches from causing bulges in the double helix. Hence, although short-range sequence-induced deformations (or deformability) are observed, long-range deformations propagated down the helix are not to be expected.     

Structural comparison of the B-DNA dodecamers d(CGCGTTAACGCG) and d(CGCGAATTCGCG) with T2A2 and A2T2 tracts.

The x-ray structure of the deoxy oligonucleotide dodecamer d(CGCGTTAACGCG) recently determined in our laboratory shows that the helical parameters of the central TTAA segment are significantly different compared to the central AATT in d(CGCGAATTCGCG). The roll in the central TA step of the T2A2 dodecamer opens towards the minor groove while the AT step of the A2T2 dodecamer opens towards the major groove. Also, the roll angles at the steps 4 and 8 (GT and AC in T2A2) and (GA and TC in A2T2) are in opposite directions. The high cup and helical twist angles at the central base-pair of T2A2 decreases the base stacking interactions compared to A2T2. Tilt angles within the tetranucleotide segments TTAA and AATT have opposite signs. In spite of the local differences caused by the sequence inversion (TTAA----AATT), the two dodecamers exhibit similar overall bending. The top third is more bent than the bottom third relative to the central segment. This asymmetric bending in the two dodecamers is mainly due to crystal packing interactions.     

DNA bending: the prevalence of kinkiness and the virtues of normality.

DNA bending in 86 complexes with sequence-specific proteins has been examined using normal vector plots, matrices of normal vector angles between all base pairs in the helix, and one-digit roll/slide/twist tables. FREEHELIX, a new program especially designed to analyze severely bent and kinked duplexes, generates the foregoing quantities plus local roll, tilt, twist, slide, shift and rise parameters that are completely free of any assumptions about an overall helix axis. In nearly every case, bending results from positive roll at pyrimidine-purine base pair steps: C-A (= T-G), T-A, or less frequently C-G, in a direction that compresses the major groove. Normal vector plots reveal three well-defined types of bending among the 86 examples: (i) localized kinks produced by positive roll at one or two discrete base pairs steps, (ii) three-dimensional writhe resulting from positive roll at a series of adjacent base pairs steps, or (iii) continuous curvature produced by alternations of positive and negative roll every 5 bp, with side-to-side zig-zag roll at intermediate position. In no case is tilt a significant component of the bending process. In sequences with two localized kinks, such as CAP and IHF, the dihedral angle formed by the three helix segments is a linear function of the number of base pair steps between kinks: dihedral angle = 36 degrees x kink separation. Twenty-eight of the 86 examples can be described as major bends, and significant elements in the recognition of a given base sequence by protein. But even the minor bends play a role in fine-tuning protein/DNA interactions. Sequence-dependent helix deformability is an important component of protein/DNA recognition, alongside the more generally recognized patterns of hydrogen bonding. The combination of FREEHELIX, normal vector plots, full vector angle matrices, and one-digit roll/slide/twist tables affords a rapid and convenient method for assessing bending in DNA.     

Analysis of local helix bending in crystal structures of DNA oligonucleotides and DNA-protein complexes.

Sequence-dependent bending of the helical axes in 112 oligonucleotide duplex crystal structures resident in the Nucleic Acid Database have been analyzed and compared with the use of bending dials, a computer graphics tool. Our analysis includes structures of both A and B forms of DNA and considers both uncomplexed forms of the double helix as well as those bound to drugs and proteins. The patterns in bending preferences in the crystal structures are analyzed by base pair steps, and emerging trends are noted. Analysis of the 66 B-form structures in the Nucleic Acid Database indicates that uniform trends within all pyrimidine-purine and purine-pyrimidine steps are not necessarily observed but are found particularly at CG and GC steps of dodecamers. The results support the idea that AA steps are relatively straight and that larger roll bends occur at or near the junctions of these A-tracts with their flanking sequences. The data on 16 available crystal structures of protein-DNA complexes indicate that the majority of the DNA bends induced via protein binding are sharp localized kinks. The analysis of the 30 available A-form DNA structures indicates that these structures are also bent and show a definitive preference for bending into the deep major groove over the shallow minor groove.     

The structure of B-helical C-G-A-T-C-G-A-T-C-G and comparison with C-C-A-A-C-G-T-T-G-G. The effect of base pair reversals

The crystal structure of the DNA decamer C-G-A-T-C-G-A-T-C-G has been solved to a resolution of 1.5 A, with a final R-factor of 16.1% for 5,107 two-sigma reflections. Crystals are orthorhombic space group P2(1)2(1)2(1), with cell dimensions a = 38.93 A, b = 39.63 A, c = 33.30 A, and 10 base pairs/asymmetric unit. The final structure contains 404 DNA atoms, 142 water molecules treated as oxygen atoms, and two Mg(H2O)6(2+) complexes. Decamers stack atop one another to simulate continuous helical columns through the crystal, as with three previously solved monoclinic decamers, but the lateral contacts between columns are quite different in the orthorhombic and monoclinic cells. Narrow and wide regions of the minor groove exhibit a single spine or two ribbons of hydration, respectively, and the minor groove is widest when BII phosphate conformations are opposed diagonally across the groove. Phosphate conformation, in turn, appears to have a base sequence dependence. Twist, rise, cup, and roll are linked as has been observed in the three monoclinic decamers and can be characterized by high or low twist profiles. In all five known decamer crystal structures and eight representative dodecamers, a high twist profile is observed with G-C and G-A steps whereas all other R-R steps are low twist profiles (R = purine). A-T and A-C steps are intermediate in character whereas C-A and C-G exhibit behavior that is strongly influenced by the profiles of the preceding and following steps. When sufficient data are in hand, sequence/structure relationships for all helix parameters probably should be considered in a 4-base pair context. At this stage of limited information the problem is compounded because there are 136 unique 4-base steps x-A-B-y in a double helix as compared with only 10 2-base steps A-B.     

Structure of a B-DNA decamer with a central T-A step: C-G-A-T-T-A-A-T-C-G.

The X-ray crystal structure analysis of the decamer C-G-A-T-T-A-A-T-C-G has been carried out to a resolution of 1.5 A. The crystals are space group P2(1)2(1)2(1), cell dimensions a = 38.60 A, b = 39.10 A, c = 33.07 A. The structure was solved by molecular replacement and refined with X-PLOR and NUCLSQ. The final R factor for a model with 404 DNA atoms, 108 water molecules and one magnesium hexahydrate cation is 15.7%. The double helix is essentially isostructural with C-G-A-T-C-G-A-T-C-G, with closely similar local helix parameters. The structure of the T-T-A-A center differs from that found in C-G-C-G-T-T-A-A-C-G-C-G in that the minor groove in our decamer is wide at the central T-A step rather than narrow, and the twist angle of the T-A step is small (31.1 degrees) rather than large. Whereas the tetrad model provides a convenient framework for discussing local DNA helix structure, it cannot be the entire story. The articulated helix model of DNA structure proposes that certain sequence regions of DNA show preferential twisting or bending properties, whereas other regions are less capable of deformation, in a manner that may be useful in sequence recognition by drugs and protein. Further crystal structure analyses should help to delineate the precise nature of sequence-dependent articulation in the DNA double helix.     

A-Tract bending: insights into experimental structures by computational models

While solution structures of adenine tract (A-tract) oligomers have indicated a unique bend direction equivalent to negative global roll (commonly termed "minor-groove bending"), crystallographic data have not unambiguously characterized the bend direction; nevertheless, many features are shared by all A-tract crystal and solution structures (e.g. propeller twisting, narrow minor grooves, and localized water spines). To examine the origin of bending and to relate findings to the crystallographic and solution data, we analyze molecular dynamics trajectories of two solvated A-tract dodecamers: 1D89, d(CGCGA(6)CG), and 1D98, d(CGCA(6)GCG), using a new general global bending framework for analyzing bent DNA and DNA/protein complexes. It is significant that the crystallographically-based initial structures are converted from dissimilar to similar bend directions equivalent to negative global roll, with the average helical-axis bend ranging from 10.5 degrees to 14.1 degrees. The largest bend occurs as positive roll of 12 degrees on the 5' side of the A-tracts (supporting a junction model) and is reinforced by gradual curvature at each A-tract base-pair (bp) step (supporting a wedge model). The precise magnitude of the bend is subtly sequence dependent (consistent with a curved general sequence model). The conversion to negative global roll only requires small local changes at each bp, accumulated over flexible moieties both outside and inside the A-tract. In contrast, the control sequence 1BNA, d(CGCGA(2)TTCGCG), bends marginally (only 6.9 degrees ) with no preferred direction. The molecular features that stabilize the bend direction in the A-tract dodecamers include propeller twisting of AT base-pairs, puckering differences between A and T deoxyriboses, a narrow minor groove, and a stable water spine (that extends slightly beyond the A-tract, with lifetimes approaching 0.2 ns). The sugar conformations, in particular, are proposed as important factors that support bent DNA. It is significant that all these curvature-stabilizing features are also observed in the crystallographic structures, but yield overall different bending paths, largely due to the effects of sequences outside the A-tract. These results merge structural details reported for A-tract structures by experiment and theory and lead to structural and dynamic insights into sequence-dependent DNA flexibility, as highlighted by the effect of an A-tract variant of a TATA-box element on bending and flexibility required for TBP binding.     

A Stochastic Model for Helix Bending in B-DNA

Abstract

Bending in double-helical B-DNA apparently occurs only by rolling adjacent base pairs over one another along their long axes. The lifting apart of ends that would be required by tilt or wedge angle contributions is too costly in free energy and does not occur. Roll angles at base steps can be positive (compression of major groove) or negative (compression of minor groove); with the former somewhat easier.

Individual steps may advance or oppose the overall direction of bend, or make lateral excursions, but the result of this series of “random roll” steps is the production of a net bending in the helix axis. Because the natural roll points for bending in a given plane occur every 5 base pairs, one would expect that double-helical DNA wrapped around a nucleosome core would exhibit bends with the same periodicity. Alternate bends might be particularly acute where the major groove faced the nucleosome core and was compressed against it.

The “annealed kinking” model proposed by Fratini et al. (J. Biol. Chem. 257, 14686 (1982) was suggested from the observation that a major bend at a natural roll point is flanked by decreasing roll angles at the steps to either side, as though local strain was being minimized by somewhat blurring the bend out rather than keeping it localized. The random walk model suggested in this paper would describe this as a decreased roll angle as the helix step rotates toward a direction perpendicular to the overall bend. Bending of DNA is seen to be a more stochastic process than had been suspected. Detailed analysis of every helix step reveals both side excursions and backward or retrograde motion, as in any random walk situation. Yet these isolated steps counteract one another, to leave behind a residuum of overall bending in a specific direction.     

Solution structure of an A-tract DNA bend

The solution structure of a DNA dodecamer d(GGCAAAAAACGG)/d(CCGTTTTTTGCC) containing an A-tract has been determined by NMR spectroscopy with residual dipolar couplings. The structure shows an overall helix axis bend of 19 degrees in a geometry consistent with solution and gel electrophoresis experiments. Fourteen degrees of the bending occurs in the GC regions flanking the A-tract. The remaining 5 degrees is spread evenly over its six AT base-pairs. The A-tract is characterized by decreasing minor groove width from the 5' to the 3' direction along the A strand. This is a result of propeller twist in the AT pairs and the increasing negative inclination of the adenine bases at the 3' side of the run of adenine bases. The four central thymine bases all have negative inclination throughout the A-tract with an average value of -6.1 degrees. Although this negative inclination makes the geometry of the A-tract different from all X-ray structures, the proton on N6 of adenine and the O4 of thymine one step down the helix are within distance to form bifurcated hydrogen bonds. The 5' bend of 4 degrees occurs at the junction between the GC flank and the A-tract through a combination of tilt and roll. The larger 3' bend, 10 degrees, occurs in two base steps: the first composed of tilt, -4.1 degrees, and the second a combination of tilt, -4.2 degrees, and roll, 6.0 degrees. This second step is a direct consequence of the change in inclination between an adjacent cytosine base, which has an inclination of -12 degrees, and the next base, a guanine, which has 3 degrees inclination. This bend is a combination of tilt and roll. The large change in inclination allows the formation of a hydrogen bond between the protons of N4 of the 3' cytosine and the O6 of the next 3' base, a guanine, stabilizing the roll component in the bend. These structural features differ from existing models for A-tract bends.For comparison, we also determined the structure of the control sequence, d(GGCAAGAAACGG)/d(CCGTTTCTTGCC), with an AT to GC transition in the center of the A-tract. This structure has no negative inclination in most of the bases within the A-tract, resulting in a bend of only 9 degrees. When ligated in phase, the control sequence has nearly normal mobility in gel electrophoresis experiments.     

A unified model for the origin of DNA sequence-directed curvature

The fine structure of the DNA double helix and a number of its physical properties depend upon nucleotide sequence. This includes minor groove width, the propensity to undergo the B-form to A-form transition, sequence-directed curvature, and cation localization. Despite the multitude of studies conducted on DNA, it is still difficult to appreciate how these fundamental properties are linked to each other at the level of nucleotide sequence. We demonstrate that several sequence-dependent properties of DNA can be attributed, at least in part, to the sequence-specific localization of cations in the major and minor grooves. We also show that effects of cation localization on DNA structure are easier to understand if we divide all DNA sequences into three principal groups: A-tracts, G-tracts, and generic DNA. The A-tract group of sequences has a peculiar helical structure (i.e., B*-form) with an unusually narrow minor groove and high base-pair propeller twist. Both experimental and theoretical studies have provided evidence that the B*-form helical structure of A-tracts requires cations to be localized in the minor groove. G-tracts, on the other hand, have a propensity to undergo the B-form to A-form transition with increasing ionic strength. This property of G-tracts is directly connected to the observation that cations are preferentially localized in the major groove of G-tract sequences. Generic DNA, which represents the vast majority of DNA sequences, has a more balanced occupation of the major and minor grooves by cations than A-tracts or G-tracts and is thereby stabilized in the canonical B-form helix. Thus, DNA secondary structure can be viewed as a tug of war between the major and minor grooves for cations, with A-tracts and G-tracts each having one groove that dominates the other for cation localization. Finally, the sequence-directed curvature caused by A-tracts and G-tracts can, in both cases, be explained by the cation-dependent mismatch of A-tract and G-tract helical structures with the canonical B-form helix of generic DNA (i.e., a cation-dependent junction model).     

Novel Hoogsteen-like bases for configurational recognition of the T-A base pair by DNA triplex formation

Effective sequence-specific recognition of duplex DNA is possible by triplex formation with natural oligonucleotides via Hoogsteen H-bonding. However, triplex formation is in practice limited to pyrimidine oligonucleotides binding duplex A-T or G-C base-pair DNA sequences specifically at homopurine sites in the major groove as T·A-T and C⁺·G-C triplets. Here we report the successful modeling of novel unnatural nucleosides that recognize the T-A DNA base pair by Hoogsteen interaction. Since the DNA triplex can be considered to assume an A-type or B-type conformation, these novel Hoogsteen nucleotides are tested within model A-type and B-type conformation triplex structures. A triplet consisting of the T-A base pair and one of the novel Hoogsteen nucleotides replaces the central T·A-T triplet in the triplex using the same deoxyribose-phosphodiester and base-deoxyribose dihedral angle configuration. The entire triplex is energy minimized and the presence of any structural or energetic perturbations due to the central triplet is assessed with respect to the unmodified energy-minimized (T·A-T)₁₁ proposed starting structures. Incorporation of these novel triplets into both A-type and B-type natural triplex structures provokes minimal change in the configuration of the central and adjacent triplets. The plan is to produce a series of Hoogsteen-like bases that preferentially bind the T-A major groove in either an A-type or B-type conformation. Selective recognition of the T-A major groove with respect to the G-C major groove, which presents similar keto and amine placement, is also assessed with configurational preference. Evaluation of the triplex solution structure by using these unnatural bases as binding conformational probes is a prerequisite to the further design of triplet forming bases. © 1996 John Wiley & Sons, Inc.     

DNA A-tracts bending: polarization effects on electrostatic interactions across their minor groove

Bending by the DNA A-tracts constitutes a contentious issue, suggesting deficiencies in the physics employed so far. Here, we inquire as to the importance in this bending of many-body polarization effects on the electrostatic interactions across their narrow minor groove. We have done this on the basis of the findings of Jarque and Buckingham who developed a procedure based on a Monte Carlo simulation for two charges of the same sign embedded in a polarizable medium. Remarkably, the present analysis reveals that for compact DNA conformations, which result from dynamic effects, an overall attractive interaction operates between the phosphate charges; this interaction is especially strong for the narrow minor groove of the A-tracts, suggesting a tendency for DNA to bend toward this groove. This tendency is in agreement with the conclusions of electrophoretic and NMR solution studies. The present analysis is also consistent with the experimental observations that the minor groove is much more easily compressible than the major groove and the bending propensity of the A-tracts is greatly reduced at “premelting” temperatures. By contrast, the dielectric screening model predicts a repulsion between the phosphate charges and is not consistent with the aforementioned bending tendency or experimental observations.     

Uranyl photoprobing of nonbent A/T- and bent A-tracts. A difference of flexibility?

In this study, we have systematically compared the uranyl photocleavage of a range of bent A-tracts and nonbent TA-tracts as well as interrupted A-tracts. We demonstrate that uranyl photocleavage of A-tracts and TA-tracts is almost identical, indicating a very similar minor groove conformation. Furthermore, a 10 base pair A-tract is divided into two independent tracts by an intervening TA or GC step. Uranyl probing also clearly distinguishes the bent A4T4 and the nonbent T4A4 sequences as adopting different structures, and our interpretation of the data is consistent with a structure for the bent A4T4 sequence that resembles a continuous A-tract, whereas the nonbent T4A4 sequences are closer to two independent and opposite A-tracts that cancel each other in terms of macroscopic bending. Finally, we also note that even single TA and TAT steps are highly sensitive to uranyl photocleavage and propose that in addition to average minor groove width, uranyl also senses DNA helix flexibility/deformability. Thus, the structural difference of TA-tracts and A-tracts may to a large extent reflect a difference in flexibility, and DNA curvature may consequently require a rigid narrow minor groove conformation that creates distinct A-tract-B-DNA junctions as the predominant cause of the bending.     

Sequence-dependent variation in DNA minor groove width dictates orientational preference of Hoechst 33258 in A-tract recognition: solution NMR structure of the 2:1 complex with d(CTTTTGCAAAAG)(2)

The solution structure of the dodecamer duplex d(CTTTTGCAAAAG)₂ and its 2:1 complex with the bis-benzimidazole Hoechst 33258 has been investigated by NMR and NOE-restrained molecular dynamics (rMD) simulations. Drug molecules are bound in each of the two A-tracts with the bulky N-methylpiperazine ring of each drug located close to the central TG (CA) step, binding essentially to the narrow minor groove of each A-tract. MD simulations over 1 ns, using an explicit solvation model, reveal time-averaged sequence-dependent narrowing of the minor groove from the 3′-end towards the 5′-end of each TTTT sequence. Distinct junctions at the TpG (CpA) steps, characterised by large positive roll, low helical and propeller twists and rapid AT base pair opening rates, add to the widening of the groove at these sites and appear to account for the bound orientation of the two drug molecules with the N-methylpiperazine ring binding in the wider part of the groove close to the junctions. Comparisons between the free DNA structure and the 2:1 complex (heavy atom RMSD 1.55 Å) reveal that these sequence-dependent features persist in both structures. NMR studies of the sequence d(GAAAAGCTTTTC)₂, in which the A-tracts have been inverted with the elimination of the TpG junctions, results in loss of orientational specificity of Hoechst 33258 and formation of multiple bound species in solution, consistent with the drug binding in a number of different orientations.     

Structure of a B-DNA dodecamer. II. Influence of base sequence on helix structure

Detailed examination of the structure of the B-DNA dodecamer C-G-C-G-A-A-T-T-C-G-C-G, obtained by single-crystal X-ray analysis (Drew et al., 1981), reveals that the local helix parameters, twist, tilt and roll, are much more strongly influenced by base sequence than by crystal packing or any other external forces. The central EcoRI restriction endonuclease recognition site, G-A-A-T-T-C, is a B helix with an average of 9.8 base-pairs per turn. It is flanked on either side by single-base-pair steps having aspects of an A-like helix character. The dodecamer structure suggests several general principles, whose validity must be tested by other B-DNA analyses. (1) When an external bending moment is applied to a B-DNA double helix, it bends smoothly, without kinks or breaks, and with relatively little effect on local helix parameters. (2) Purine-3′,5′-pyrimidine steps open their base planes towards the major groove, pyrimidine-purine steps open toward the minor groove, and homopolymer (Pur-Pur, Pyr-Pyr) steps resist rolling in either direction. This behavior is related to the preference of pyrimidines for more negative glycosyl torsion angles. (3) CpG steps have smaller helical twist angles than do GpC, as though in compensation for their smaller intrinsic base overlap. Data on A-T steps are insufficient for generalization. (4) G.C base-pairs have smaller propellor twist than A · T, and this arises mainly from interstrand base overlap rather than the presence of the third hydrogen bond. (5) DNAase I cuts preferentially at positions of high helical twist, perhaps because of increased exposure of the backbone to attack. The correlation of the digestion patterns in solution and helical twist in the crystal argues for the essential identity of the helix structure in the two environments. (6) In the two places where the sequence TpCpG occurs, the C slips from under T in order to stack more efficiently over G. At the paired bases of this CpG step, the G and C are tilted so the angle between base planes is splayed out to the outside of the helix. This TpC is the most favored cutting site for DNAase I by a factor of 4.5 (Lomonossoff et al., 1981). (7) The EcoRI restriction endonuclease and methylase both appear to prefer a cutting site of the type purine-purine-A-T-T-pyrimidine, involving two adjacent homopolymer triplets, and this may be a consequence of the relative stiffness of homopolymer base-stacking observed in the dodecamer.     

Induced fit DNA recognition by a minor groove binding analogue of Hoechst 33258: fluctuations in DNA A tract structure investigated by NMR and molecular dynamics simulations

NMR analysis and molecular dynamics simulations of d(GGTAATTACC)₂ and its complex with a tetrahydropyrimidinium analogue of Hoechst 33258 suggest that DNA minor groove recognition in solution involves a combination of conformational selection and induced fit, rather than binding to a preorganised site. Analysis of structural fluctuations in the bound and unbound states suggests that the degree of induced fit observed is primarily a consequence of optimising van der Waals contacts with the walls of the minor groove resulting in groove narrowing through: (i) changes in base step parameters, including increased helical twist and propeller twist; (ii) changes to the sugar–phosphate backbone conformation to engulf the bound ligand; (iii) suppression of bending modes at the TpA steps. In contrast, the geometrical arrangement of hydrogen bond acceptors on the groove floor appears to be relatively insensitive to DNA conformation (helical twist and propeller twist). We suggest that effective recognition of DNA sequences (in this case an A tract structure) appears to depend to a significant extent on the sequence being flexible enough to be able to adopt the geometrically optimal conformation compatible with the various binding interactions, rather than involving ‘lock and key’ recognition.     

(+)-CC-1065 as a probe for intrinsic and protein-induced bending of DNA

(+)-CC -1065 is biologically potent DNA-reactive antitumor antibiotic produced by Streptomyces zelensis. This antibiotic covalently modifies DNA by alkylation of N-3 of a adenine in the minor groove. As a Structural consequence of covalent modification of DNA, the helix axis id bent into the minor groove. The drug-induced bending of DNA has similarities to intrinsic. A-tract bending and the 3′ adenine of A-tracts shows a unique reactivity to alkylation by (+) -CC-1065. Upon covalent modification of A-tracts, the magnitude of bending is increased and helix is stiffened. Using high-field NMR, hydroxyl-radical footprinting and gel electrophoresis, the molecular basis for the high reactivity of the bonding sequence 5′ - AGTTA* (an asterisk indicates the covalent modification site) to (+)-CC-1065 has been shown to involve the inherent conformational flexibility of this sequence. Furthermore, these studies also demonstrate that after alkylation the drug-induced bending is focused over the TT region. By analogy with the junction bend model for A-tracts, a ‘truncated junction bend model’ is proposed for this structure. Last, the application of (+)-CC-1065 entrapped/induced bending of DNA as a probe for the Sp1-induced bending of the 21-base-pair repeat an Mu transpose bending of the att L3 sequence is described.     

Sequence dependence of DNA conformational flexibility

By using conformational free energy calculations, we have studied the sequence dependence of flexibility and its anisotropy along various conformational variables of DNA base pairs. The results show the AT base step to be very flexible along the twist coordinate. On the other hand, homonucleotide steps, GG(CC) and AA(TT), are among the most rigid sequences. For the roll motion that would correspond to a bend, the TA step is most flexible, while the GG(CC) step is least flexible. The flexibility of roll is quite anisotropic; the ratio of fluctuations toward the major and minor grooves is the largest for the GC step and the smallest for the AA(TT) and CG steps. Propeller twisting of base pairs is quite flexible, especially of A.T base pairs; propeller twist can reach 19 degrees by thermal fluctuation. We discuss the effect of electrostatic parameters, comparison with available experimental results, and biological relevance of these results.     

Indirect readout of DNA sequence at the primary-kink site in the CAP-DNA complex: recognition of pyrimidine-purine and purine-purine steps

The catabolite activator protein (CAP) bends DNA in the CAP-DNA complex, typically introducing a sharp DNA kink, with a roll angle of approximately 40 degrees and a twist angle of approximately 20 degrees, between positions 6 and 7 of the DNA half-site, 5'-A1A2A3T4G5T6G7A8T9C10T11 -3' ("primary kink"). In previous work, we showed that CAP recognizes the nucleotide immediately 5' to the primary-kink site, T6, through an "indirect-readout" mechanism involving sequence effects on energetics of primary-kink formation. Here, to understand further this example of indirect readout, we have determined crystal structures of CAP-DNA complexes containing each possible nucleotide at position 6. The structures show that CAP can introduce a DNA kink at the primary-kink site with any nucleotide at position 6. The DNA kink is sharp with the consensus pyrimidine-purine step T6G7 and the non-consensus pyrimidine-purine step C6G7 (roll angles of approximately 42 degrees, twist angles of approximately 16 degrees ), but is much less sharp with the non-consensus purine-purine steps A6G7 and G6G7 (roll angles of approximately 20 degrees, twist angles of approximately 17 degrees). We infer that CAP discriminates between consensus and non-consensus pyrimidine-purine steps at positions 6-7 solely based on differences in the energetics of DNA deformation, but that CAP discriminates between the consensus pyrimidine-purine step and non-consensus purine-purine steps at positions 6-7 both based on differences in the energetics of DNA deformation and based on qualitative differences in DNA deformation. The structures further show that CAP can achieve a similar, approximately 46 degrees per DNA half-site, overall DNA bend through a sharp DNA kink, a less sharp DNA kink, or a smooth DNA bend. Analysis of these and other crystal structures of CAP-DNA complexes indicates that there is a large, approximately 28 degrees per DNA half-site, out-of-plane component of CAP-induced DNA bending in structures not constrained by end-to-end DNA lattice interactions and that lattice contacts involving CAP tend to involve residues in or near biologically functional surfaces.