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

6S RNA Mimics B-Form DNA to Regulate Escherichia coli RNA Polymerase

Noncoding RNAs (ncRNAs) regulate gene expression in all organisms. Bacterial 6S RNAs globally regulate transcription by binding RNA polymerase (RNAP) holoenzyme and competing with promoter DNA. Escherichia coli (Eco) 6S RNA interacts specifically with the housekeeping σ⁷⁰-holoenzyme (Eσ⁷⁰) and plays a key role in the transcriptional reprogramming upon shifts between exponential and stationary phase. Inhibition is relieved upon 6S RNA-templated RNA synthesis. We report here the 3.8 Å resolution structure of a complex between 6S RNA and Eσ⁷⁰ determined by single-particle cryo-electron microscopy and validation of the structure using footprinting and crosslinking approaches. Duplex RNA segments have A-form C3' endo sugar puckers but widened major groove widths, giving the RNA an overall architecture that mimics B-form promoter DNA. Our results help explain the specificity of Eco 6S RNA for Eσ⁷⁰ and show how an ncRNA can mimic B-form DNA to directly regulate transcription by the DNA-dependent RNAP.

Mg2+ in the major groove modulates B-DNA structure and dynamics

This study investigates the effect of Mg(2+) bound to the DNA major groove on DNA structure and dynamics. The analysis of a comprehensive dataset of B-DNA crystallographic structures shows that divalent cations are preferentially located in the DNA major groove where they interact with successive bases of (A/G)pG and the phosphate group of 5'-CpA or TpG. Based on this knowledge, molecular dynamics simulations were carried out on a DNA oligomer without or with Mg(2+) close to an ApG step. These simulations showed that the hydrated Mg(2+) forms a stable intra-strand cross-link between the two purines in solution. ApG generates an electrostatic potential in the major groove that is particularly attractive for cations; its intrinsic conformation is well-adapted to the formation of water-mediated hydrogen bonds with Mg(2+). The binding of Mg(2+) modulates the behavior of the 5'-neighboring step by increasing the BII (ε-ζ>0°) population of its phosphate group. Additional electrostatic interactions between the 5'-phosphate group and Mg(2+) strengthen both the DNA-cation binding and the BII character of the 5'-step. Cation binding in the major groove may therefore locally influence the DNA conformational landscape, suggesting a possible avenue for better understanding how strong DNA distortions can be stabilized in protein-DNA complexes.

Electron donor-acceptor interactions with flanking purines influence the efficiency of thymine photodimerization

Quantum yields for thymine photodimerization (Φ(TT)) have been determined for a series of short DNA single-strand and base-paired hairpin structures possessing a single thymine-thymine step with flanking purines. Values of Φ(TT) are strongly dependent upon the oxidation potential of the flanking purine, decreasing in the order: inosine > adenine > guanine > deazaguanine. The dependence of Φ(TT) on the ionization potential of the flanking purine is more pronounced when the purine of lower oxidation potential is located at the 5'- versus 3'-position in either a single strand or a hairpin. Molecular dynamics simulations for hairpin structures indicate that the TT step is π-stacked with both the 5' and 3' purine, but that there is little π-stacking with either purine in single-strand structures. The observation of moderately intense long-wavelength UV absorption features for hairpins having 5'-Z or G flanking purines suggests that excitation of ground state donor-acceptor complexes may account for more extensive quenching of dimerization by 5'- versus 3'-purines. The "purine effect" on Φ(TT) is attributed to a combination of ground state conformation, ground state electron donor-acceptor interactions, and excited state exciplex formation.

Local elasticity of strained DNA studied by all-atom simulations

Genomic DNA is constantly subjected to various mechanical stresses arising from its biological functions and cell packaging. If the local mechanical properties of DNA change under torsional and tensional stress, the activity of DNA-modifying proteins and transcription factors can be affected and regulated allosterically. To check this possibility, appropriate steady forces and torques were applied in the course of all-atom molecular dynamics simulations of DNA with AT- and GC-alternating sequences. It is found that the stretching rigidity grows with tension as well as twisting. The torsional rigidity is not affected by stretching, but it varies with twisting very strongly, and differently for the two sequences. Surprisingly, for AT-alternating DNA it passes through a minimum with the average twist close to the experimental value in solution. For this fragment, but not for the GC-alternating sequence, the bending rigidity noticeably changes with both twisting and stretching. The results have important biological implications and shed light on earlier experimental observations.

Sequence-specific B-DNA flexibility modulates Z-DNA formation

Conversion of right-handed B-DNA into left-handed Z-DNA is one of the largest structural transitions in biology that plays fundamental roles in gene expression and regulation. Z-DNA segments must form within genomes surrounded by a sea of B-DNA and require creation of energetically costly B/Z junctions. Here, we show using a combination of natural abundance NMR R(1ρ) carbon relaxation measurements and CD spectroscopy that sequence-specific B-DNA flexibility modulates the thermodynamic propensity to form Z-DNA and the location of B/Z junctions. We observe sequence-specific flexibility in B-DNA spanning fast (ps-ns) and slow (μs-ms) time scales localized at the site of B/Z junction formation. Further, our studies show that CG-repeats play an active role tuning this intrinsic B-DNA flexibility. Taken together, our results suggest that sequence-specific B-DNA flexibility may provide a mechanism for defining the length and location of Z-DNA in genomes.

Extensive protein and DNA backbone sampling improves structure-based specificity prediction for C2H2 zinc fingers

Sequence-specific DNA recognition by gene regulatory proteins is critical for proper cellular functioning. The ability to predict the DNA binding preferences of these regulatory proteins from their amino acid sequence would greatly aid in reconstruction of their regulatory interactions. Structural modeling provides one route to such predictions: by building accurate molecular models of regulatory proteins in complex with candidate binding sites, and estimating their relative binding affinities for these sites using a suitable potential function, it should be possible to construct DNA binding profiles. Here, we present a novel molecular modeling protocol for protein-DNA interfaces that borrows conformational sampling techniques from de novo protein structure prediction to generate a diverse ensemble of structural models from small fragments of related and unrelated protein-DNA complexes. The extensive conformational sampling is coupled with sequence space exploration so that binding preferences for the target protein can be inferred from the resulting optimized DNA sequences. We apply the algorithm to predict binding profiles for a benchmark set of eleven C₂H₂ zinc finger transcription factors, five of known and six of unknown structure. The predicted profiles are in good agreement with experimental binding data; furthermore, examination of the modeled structures gives insight into observed binding preferences.

A duplex DNA model with regular inter-base-pair hydrogen bonds

It is well known that base-pair stacking is the main factor in stabilizing DNA duplex and plays an important role in determining DNA sequence-dependence. What is the dominant force in base-pair stacking? This fundamental biological question remains a challenging problem. Here, based on recent studies about the non-planarity of amino groups on DNA bases, we propose a new duplex DNA model, in which all base amino groups are non-planar and participate in forming regular inter-base-pair hydrogen bonds (IBP H-bonds). This model implies that IBP H-bonds are the dominant force stabilizing base-pair stacking and play a crucial role in determining the geometry and physical properties of sequence-dependent twisted stacking between adjacent base pairs. The model presents a new insight into the link, through regular IBP H-bonds, between base-sequence, fine structure and physical properties at dinucleotide step level, and provides an attractively concise, uniform and quantitative interpretation for various experimentally observed DNA sequence-dependent properties in terms of regular IBP H-bonds. It would provide a new approach to understanding the dynamics and underlying mechanisms of DNA sequence-dependent biological processes, sequence-structure-property relationships, DNA strand separation during replication and transcriptions, etc.

Conformations of the sugar-phosphate backbone in helical DNA crystal structures

The DNA backbone geometry was analyzed for 96 crystal structures of oligodeoxynucleotides. The ranges and mean values of the torsion angles for the observed subclasses of the A-, B-, and Z-DNA conformational types were determined by analyzing distributions of the torsion angles and scattergrams relating pairs of angles.

Acceleration of 5-methylcytosine deamination in cyclobutane dimers by G and its implications for UV-induced C-to-T mutation hotspots

Sunlight-induced C-->T mutation hotspots occur most frequently at methylated CpG sites in tumor suppressor genes and are thought to arise from translesion synthesis past deaminated cyclobutane pyrimidine dimers (CPDs). While it is known that methylation enhances CPD formation in sunlight, little is known about the effect of methylation and sequence context on the deamination of 5-methylcytosine ((m)C) and its contribution to mutagenesis at these hotspots. Using an enzymatic method, we have determined the yields and deamination rates of C and (m)C in CPDs and find that the frequency of UVB-induced CPDs correlates with the oxidation potential of the flanking bases. We also found that the deamination of T(m)C and (m)CT CPDs is about 25-fold faster when flanked by G's than by A's, C's or T's in duplex DNA and appears to involve catalysis by the O6 group of guanine. In contrast, the first deamination of either C or (m)C in AC(m)CG with a flanking G was much slower (t(1/2) >250 h) and rate limiting, while the second deamination was much faster. The observation that C(m)CG dimers deaminate very slowly but at the same time correlate with C-->T mutation hotspots suggests that their repair must be slow enough to allow sufficient time for deamination. There are, however, a greater number of single C-->T mutations than CC-->TT mutations at C(m)CG sites even though the second deamination is very fast, which could reflect faster repair of doubly deaminated dimers.

The structure of LNA:DNA hybrids from molecular dynamics simulations: the effect of locked nucleotides

Locked nucleic acids (LNAs) exhibit a modified sugar fragment that is restrained to the C3'-endo conformation. LNA-containing duplexes are rather stable and have a more rigid structure than DNA duplexes, with a propensity for A-conformation of the double helix. To gain detailed insight into the local structure of LNA-modified DNA oligomers (as a foundation for subsequent exploration of the electron-transfer capabilities of such modified duplexes), we carried out molecular dynamics simulations on a set of LNA:DNA 9-mer duplexes and analyzed the resulting structures in terms of base step parameters and the conformations of the sugar residues. The perturbation introduced by a single locked nucleotide was found to be fairly localized, extending mostly to the first neighboring base pairs; such duplexes featured a B-type helix. With increasing degree of LNA modification the structure gradually changed; the duplex with one complete LNA strand assumed a typical A-DNA structure. The relative populations of the sugar conformations agreed very well with NMR data, lending credibility to the validity of the computational protocol.

Effect of cesium on the volume of the helix-coil transition of dA.dT polymers and their ligand complexes

The pressure dependence of the helix-coil transition of poly(dA).poly(dT) and poly[d(A-T)].poly[d(A-T)] in aqueous solutions of NaCl and CsCl at concentrations between 10 and 200 mM is reported and used to calculate the accompanying volume change. We also investigated the binding parameters and volume change of ethidium bromide binding with poly(dA).poly(dT) and poly[d(A-T)].poly[d(A-T)] in aqueous solutions of these two salts. The volume change of helix-coil transition of poly(dA).poly(dT) in Cs(+)-containing solutions differs by less than 1 cm(3) mol(-1) from the value measured when Na(+) is the counter-ion. We propose that this insensitivity towards salt type arises if the counter-ions are essentially fully hydrated around DNA and the DNA conformation is not significantly altered by salt types. Circular dichroism spectroscopy showed that the previously observed large volumetric disparity for the helix-coil transition of poly[d(A-T)].poly[d(A-T)] in solutions containing Na(+) and Cs(+) is likely result of a Cs(+)-induced conformation change that is specific for poly[d(A-T)].poly[d(A-T)]. This cation-specific conformation difference is mostly absent for poly(dA).poly(dT) and EB bound poly[d(A-T)].poly[d(A-T)].

Temperature dependence of the volumetric parameters of drug binding to poly[d(A-T)].Poly[d(A-T)] and Poly(dA).Poly(dT)

We report the temperature and salt dependence of the volume change (DeltaVb) associated with the binding of ethidium bromide and netropsin with poly(dA).poly(dT) and poly[d(A-T)].poly[d(A-T)]. The DeltaV(b) of binding of ethidium with poly(dA).poly(dT) was much more negative at temperatures approximately 70 degrees C than at 25 degrees C, whereas the difference is much smaller in the case of binding with poly[d(A-T)].poly[d(A-T)]. We also determined the volume change of DNA-drug interaction by comparing the volume change of melting of DNA duplex and DNA-drug complex. The DNA-drug complexes display helix-coil transition temperatures (Tm several degrees above those of the unbound polymers, e.g., the Tm of the netropsin complex with poly(dA)poly(dT) is 106 degrees C. The results for the binding of ethidium with poly[d(A-T)].poly[d(A-T)] were accurately described by scaled particle theory. However, this analysis did not yield results consistent with our data for ethidium binding with poly(dA).poly(dT). We hypothesize that heat-induced changes in conformation and hydration of this polymer are responsible for this behavior. The volumetric properties of poly(dA).poly(dT) become similar to those of poly[d(A-T)].poly[d(A-T)] at higher temperatures.

TBP flanking sequences: asymmetry of binding, long-range effects and consensus sequences

We carried out in vitro selection experiments to systematically probe the effects of TATA-box flanking sequences on its interaction with the TATA-box binding protein (TBP). This study validates our previous hypothesis that the effect of the flanking sequences on TBP/TATA-box interactions is much more significant when the TATA box has a context-dependent DNA structure. Several interesting observations, with implications for protein-DNA interactions in general, came out of this study. (i) Selected sequences are selection-method specific and TATA-box dependent. (ii) The variability in binding stability as a function of the flanking sequences for (T-A)4 boxes is as large as the variability in binding stability as a function of the core TATA box itself. Thus, for (T-A)4 boxes the flanking sequences completely dominate and determine the binding interaction. (iii) Binding stabilities of all but one of the individual selected sequences of the (T-A)4 form is significantly higher than that of their mononucleotide-based consensus sequence. (iv) Even though the (T-A)4 sequence is symmetric the flanking sequence pattern is asymmetric. We propose that the plasticity of (T-A)n sequences increases the number of conformationally distinct TATA boxes without the need to extent the TBP contact region beyond the eight-base-pair long TATA box.

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.

Interfacial water as a "hydration fingerprint" in the noncognate complex of BamHI

The molecular code of specific DNA recognition by proteins as a paradigm in molecular biology remains an unsolved puzzle primarily because of the subtle interplay between direct protein-DNA interaction and the indirect contribution from water and ions. Transformation of the nonspecific, low affinity complex to a specific, high affinity complex is accompanied by the release of interfacial water molecules. To provide insight into the conversion from the loose to the tight form, we characterized the structure and energetics of water at the protein-DNA interface of the BamHI complex with a noncognate sequence and in the specific complex. The fully hydrated models were produced with Grand Canonical Monte Carlo simulations. Proximity analysis shows that water distributions exhibit sequence dependent variations in both complexes and, in particular, in the noncognate complex they discriminate between the correct and the star site. Variations in water distributions control the number of water molecules released from a given sequence upon transformation from the loose to the tight complex as well as the local entropy contribution to the binding free energy. We propose that interfacial waters can serve as a "hydration fingerprint" of a given DNA sequence.

Stability of DNA Duplexes with Watson−Crick Base Pairs: A Predicted Model

The conformational stability (difference between the free energies of the folded and unfolded states, DeltaG degrees ) of a DNA duplex is considered as a function of component energy terms, hydrophobic, base stacking, hydrogen bonding, van der Waals, and electrostatic, and a trinucleotide-level helix stiffness parameter measured in terms of its Young's modulus. Hydrophobic and base stacking energy components were determined with the use of the crystal structure data of 30 DNA duplexes judicially selected within a resolution of 1.5 A, and hydrogen bonding, van der Waals and electrostatic terms were determined through an extensive review of experimental and theoretical studies. The stiffness indices for the trinucleotides were the ones realized by M. M. Gromiha [(2000) J. Biol. Phys. 26, 43-50] using the crystal structure data of 70 DNA duplexes. The unfolded state was treated in the classical way to determine its stability. Thermodynamically determined DeltaG degrees values for 111 DNA duplexes, with the number of base pairs ranging from 4 to 16, were selected in two sets, and the regression equation formed with one set was used to predict the stabilities of the other set, taking the energy components and the stiffness parameter to be independent variables. The computed energy terms indicate that the base stacking and hydrogen bonding forces are the dominant and the hydrophobic and electrostatic forces the weak partners in imparting stability to the duplexes. This model predicts DeltaG degrees values for DNA duplexes examined with a level of accuracy similar to that used for predictions made by the widely used nearest-neighbor models. The uniqueness of this model is that it combines the crystal and thermodynamic data for interpretation of conformational stability.

The role of water in protein-DNA recognition

Is it by design or by default that water molecules are observed at the interfaces of some protein-DNA complexes? Both experimental and theoretical studies on the thermodynamics of protein-DNA binding overwhelmingly support the extended hydrophobic view that water release from interfaces favors binding. Structural and energy analyses indicate that the waters that remain at the interfaces of protein-DNA complexes ensure liquid-state packing densities, screen the electrostatic repulsions between like charges (which seems to be by design), and in a few cases act as linkers between complementary charges on the biomolecules (which may well be by default). This review presents a survey of the current literature on water in protein-DNA complexes and a critique of various interpretations of the data in the context of the role of water in protein-DNA binding and principles of protein-DNA recognition in general.

Cooperative dimerization of a heterocyclic diamidine determines sequence-specific DNA recognition

In the course of a program aimed at discovering novel DNA-targeted antiparasitic drugs, the phenylfuran-benzimidazole unfused aromatic dication DB293 was identified as the first diamidine capable of forming stacked dimers in the DNA minor groove of GC-containing sequences. Its preferred binding sequence encompasses the tetranucleotide 5'-ATGA.5'-TCAT to which DB293 binds tightly with a strong positive cooperativity. Here we have investigated the influence of the DNA sequence on drug binding using two complementary technical approaches: surface plasmon resonance and DNase I footprinting. The central dinucleotide of the primary ATGA motif was systematically varied to represent all of the eight possible combinations (AXGA and ATYA, where X or Y = A, T, G, or C). Binding affinities for each site were precisely measured by SPR, and the extent of cooperative drug binding was also determined. The sequence recognition process was found to be extremely dependent on the nature of the central dinucleotide pair. Modification of the central TG step decreases binding affinity by a factor varying from 2 to over 500 depending on the base substitution. However, the diminished binding affinity does not affect the unique binding mode. In nearly all cases, the SPR titrations revealed a positive cooperativity in complex formation which reflects the ease of the dication to form stacked dimeric motifs in the DNA minor groove. DNase I footprinting served to identify additional binding sites for DB293 in the context of long DNA sequences offering a large variety of randomly distributed or specifically designed sites. The ATGA motif provided the best receptor for the drug, but lower affinity sequences were also identified. The design of two DNA fragments composed of various targeted tetranucleotide binding sites separated by an "insulator" (nonbinding) sequence allowed us to delineate further the influence of DNA sequence on drug binding and to identify a novel high-affinity site: 5'-ACAA.5'-TTGT. Collectively, the SPR and footprinting results show that the consensus sequence 5'-(A/T)-TG-(A/T) represents the optimal site for cooperative dimerization of the heterocyclic diamidine DB293.

Sequence-selective targeting of long stretches of the DNA minor groove by a novel dimeric bis-benzimidazole

A dimeric bis-benzimidazole molecule has been designed by computer modeling to bind to a DNA sequence via the DNA minor groove that covers a complete turn of B-DNA. A series of bis-benzimidazole dimers incorporating a -O-(CH(2))(n)()-X-CH(2))(n)()-O- linker, with n = 2 or 3 and X = O or N(+)H(Me), were screened for their capacity to fit the DNA minor groove. The modeling studies enabled an optimal linker to be devised (n = 3, X = N(+)H(Me)), and the synthesis of the predicted "best" molecule, N-methyl-N,N-bis-3,3-[4'-[5' '-(2' "-p-methoxyphenyl)-5' "-1H-benzimidazolyl]-2' '-1H-benzimidazolyl]phenoxypropylamine (5), is reported. The optimized linker permits the two symmetric bis-benzimidazole motifs to maintain hydrogen-bonded contacts with the floor of the DNA minor groove. DNase I footprinting studies have shown that this ligand binds with high affinity to sequences representing approximately a complete turn of B-DNA, represented by the [A.T](4)-[G.C]-[A.T](4) motif, and only poorly to sequences of half this site size, in accord with the computer modeling studies. Compound 5 does not show acute cellular cytotoxicity, in contrast with its monomeric bis-benzimidazole precursors, yet is rapidly taken up into cells.

DNA polymorphism: a comparison of force fields for nucleic acids

The improvements of the force fields and the more accurate treatment of long-range interactions are providing more reliable molecular dynamics simulations of nucleic acids. The abilities of certain nucleic acid force fields to represent the structural and conformational properties of nucleic acids in solution are compared. The force fields are AMBER 4.1, BMS, CHARMM22, and CHARMM27; the comparison of the latter two is the primary focus of this paper. The performance of each force field is evaluated first on its ability to reproduce the B-DNA decamer d(CGATTAATCG)(2) in solution with simulations in which the long-range electrostatics were treated by the particle mesh Ewald method; the crystal structure determined by Quintana et al. (1992) is used as the starting point for all simulations. A detailed analysis of the structural and solvation properties shows how well the different force fields can reproduce sequence-specific features. The results are compared with data from experimental and previous theoretical studies.

Characterization of divalent cation localization in the minor groove of the A(n)T(n) and T(n)A(n) DNA sequence elements by (1)H NMR spectroscopy and manganese(II)

The localization of Mn(2+) in A-tract DNA has been studied by (1)H NMR spectroscopy using a series of self-complementary dodecamer oligonucleotides that contain the sequence motifs A(n)(n) and T(n)A(n), where n = 2, 3, or 4. Mn(2+) localization in the minor groove is observed for all the sequences that have been studied, with the position and degree of localization being highly sequence-dependent. The site most favored for Mn(2+) localization in the minor groove is near the 5'-most ApA step for both the T(n)A(n) and the A(n)T(n) series. For the T(n)A(n) series, this results in two closely spaced symmetry-related Mn(2+) localization sites near the center of each duplex, while for the A(n)T(n) series, the two symmetry-related sites are separated by as much as one half-helical turn. The degree of Mn(2+) localization in the minor groove of the T(n)A(n) series decreases substantially as the AT sequence element is shortened from T(4)A(4) to T(2)A(2). The A(n)T(n) series also exhibits length-dependent Mn(2+) localization; however, the degree of minor groove occupancy by Mn(2+) is significantly less than that observed for the T(n)A(n) series. For both A(n)T(n) and T(n)A(n) sequences, the 3'-most AH2 resonance is the least broadened of the AH2 resonances. This is consistent with the observation that the minor groove of A-tract DNA narrows in the 5' to 3' direction, apparently becoming too narrow after two base pairs for the entry of a fully hydrated divalent cation. The results that are reported illustrate the delicate interplay that exists between DNA nucleotide sequence, minor groove width, and divalent cation localization. The proposed role of cation localization in helical axis bending by A-tracts is also discussed.

Solvation and hydration of proteins and nucleic acids: a theoretical view of simulation and experiment

Many theoretical, computational, and experimental techniques recently have been successfully used for description of the solvent distribution around macromolecules. In this Account, we consider recent developments in the areas of protein and nucleic acid solvation and hydration as seen by experiment, theory, and simulations. We find that in most cases not only the general phenomena of solvation but even local hydration patterns are more accurately discussed in the context of water distributions rather than individual molecules of water. While a few localized or high-residency waters are often associated with macromolecules in solution (or crystals from aqueous liquors), these are readily and accurately included in this more general description. The goal of this Account is to review the theoretical models used for the description of the interfacial solvent structure on the border near DNA and protein molecules. In particular, we hope to highlight the progress in this field over the past five years with a focus on comparison of simulation and experimental results.

Thymine-methyl/pi interaction implicated in the sequence-dependent deformability of DNA

The crystal structures of deoxy-oligonucleotides were retrieved from the Nucleic Acid Database and analyzed with the use of our program CHPI. The structure of 5'-ApTpApT-3' has been shown to be stabilized by the 5-methyl group in the thymidine moiety that favorably interacts with the adenine pi-ring preceding it. H2' of the deoxyribose in adenine also interacts with the thymine ring next to it. Since a 5'-ApT-3' sequence is accompanied by another 5'-ApT-3' in the complementary strand, the interaction is duplicated, thus forming a 'twin A/T-Me interaction'. Coordinates of oligonucleotides with A-T rich sequences were retrieved and analyzed. In almost every case, the thymidine 5-methyl group favorably interacts with an adenine ring in the same strand. The structure of duplexes incorporating A-tracts was also analyzed. The 5-methyl group in the thymidine moiety has been found to interact favorably with the base pi-ring before it. Since an A-tract is lined with an oligo-T sequence in the complementary strand, a successive N/T-Me stacking may contribute in making the A-tracts robust and straight. The possible involvement of the N/T-Me and the twin A/T-Me motif in the deformability of DNA has been suggested. The role of methyl groups in modified DNA has been discussed on a similar basis.

Influence of compound structure on affinity, sequence selectivity, and mode of binding to DNA for unfused aromatic dications related to furamidine

In the course of a program aimed at developing sequence-specific gene-regulatory small organic molecules, we have investigated the DNA interactions of a new series of nine diphenylfuran dications related to the antiparasitic drug furamidine (DB75). Two types of structural modifications were tested: the terminal amidine groups of DB75 were shifted from the para to the meta position, and the amidines were replaced with imidazoline or dimethyl-imidazoline groups, to test the importance of both the position and nature of positively charged groups on DNA interactions. The interactions of these compounds with DNA and oligonucleotides were studied by a combination of biochemical and biophysical techniques. Absorption and CD measurements suggested that the drugs bind differently to AT and GC sequences in DNA. The para-para dications, like DB75, bind into the minor groove of poly(dAT)(2) and intercalate between the base pairs of poly(dGC)(2), as revealed by electric linear dichroism experiments. In contrast, the meta-meta compounds exhibit a high tendency to intercalate into DNA whatever the target sequence. The lack of sequence selectivity of the meta-meta compounds containing amidines or dimethyl-imidazoline groups was also evident from DNase I footprinting and surface plasmon resonance (SPR) experiments. Accurate binding measurements using the BIAcore SPR method revealed that all nine compounds bind with similar affinity to an immobilized GC sequence DNA hairpin but exhibit very distinct affinities for the corresponding AT hairpin oligonucleotide. The minor groove-binding para-para compounds have a high specificity for AT sequences. The biophysical data clearly indicate that shifting the cationic substituents from the para to the meta position results in a loss of specificity and change in binding mode. The strong AT selectivity of the para-para compounds was independently confirmed by DNase I footprinting experiments performed with a range of DNA restrictions fragments. In terms of AT selectivity, the compounds rank in the order para-para > para-meta > meta-meta. The para dications bind preferentially to sequences containing four contiguous AT base pairs. Additional footprinting experiments with substrates containing the 16 possible [A.T](4) blocks indicated that the presence of a TpA step within an [A.T] (4) block generally reduces the extent of binding. The diverse methods, from footprinting to SPR to dichroism, provide a consistent model for the interactions of the diphenylfuran dications with DNA of different sequences. Altogether, the results attest unequivocally that the binding mode for unfused aromatic cations can change completely depending on substituent position and DNA sequence. These data provide a rationale to explain the relationships between sequence selectivity and mode of binding to DNA for unfused aromatic dications related to furamidine.

Base-sequence specificity of Hoechst 33258 and DAPI binding to five (A/T)4 DNA sites with kinetic evidence for more than one high-affinity Hoechst 33258-AATT complex

The binding of Hoechst 33258 and DAPI to five different (A/T)4 sequences in a stable DNA hairpin was studied exploiting the substantial increase in dye fluorescence upon binding. The two dyes have comparable affinities for the AATT site (e.g. association constant K(a)=5.5 x 10(8) M(-1) for DAPI), and their affinities decrease in the series AATT >> TAAT approximately equal to ATAT > TATA approximately equal to TTAA. The extreme values of K(a) differ by a factor of 200 for Hoechst 33258 but only 30 for DAPI. The binding kinetics of Hoechst 33258 were measured by stopped-flow under pseudo-first order conditions with an (A/T)4 site in excess. The lower-resolution experiments can be well represented by single exponential processes, corresponding to a single-step binding mechanism. The calculated association-rate parameters for the five (A/T)4 sites are similar (2.46 x 10(8) M(-1) s(-1) to 0.86 x 10(8) M(-1) s(-1)) and nearly diffusion-controlled, while the dissociation-rate parameters vary from 0.42 s(-1) to 96 s(-1). Thus the association constants are kinetically controlled and are close to their equilibrium-determined values. However, when obtained with increased signal-to-noise ratio, the kinetic traces for Hoechst 33258 binding at the AATT site reveal two components. The concentration dependencies of the two time constants and amplitudes are consistent with two different kinetically equivalent two-step models. In the first model, fast bimolecular binding is followed by an isomerization of the initial complex. In the second model, two single-step associations form two complexes that mutually exclude each other. For both models the four reaction-rate parameters are calculated. Finally, specific dissociation kinetics, using poly[d(A-5BrU)], show that the kinetics are even more complex than either two-step model. We correlate our results with the different binding orientations and locations of Hoechst 33258 in the DNA minor groove found in several structural studies in the literature.

Do water molecules mediate protein-DNA recognition?

A comprehensive analysis of interfacial water molecules in the structures of 109 unique protein-DNA complexes is presented together with a new view on their role in protein-DNA recognition. Location of interfacial water molecules as reported in the crystal structures and as emerging from a series of molecular dynamics studies on protein-DNA complexes with explicit solvent and counterions, was analyzed based on their acceptor, donor hydrogen bond relationships with the atoms and residues of the macromolecules, electrostatic field calculations and packing density considerations. Water molecules for the purpose of this study have been categorized into four classes: viz. (I) those that contact both the protein and the DNA simultaneously and thus mediate recognition directly; (II) those that contact either the protein or the DNA exclusively via hydrogen bonds solvating each solute separately; (III) those that contact the hydrophobic groups in either the protein or the DNA; and, lastly (IV) those that contact another water molecule. Of the 17,963 crystallographic water molecules under examination, about 6% belong to class I and 76% belong to class II. About three-fourths of class I and class II water molecules are exclusively associated with hydrogen bond acceptor atoms of both protein and DNA. Noting that DNA is polyanionic, it is significant that a majority of the crystallographically observed water molecules as well as those from molecular dynamics simulations should be involved in facilitating binding by screening unfavorable electrostatics. Less than 2% of the reported water molecules occur between hydrogen bond donor atoms of protein and acceptor atoms of DNA. These represent cases where protein atoms cannot reach out to DNA to make favorable hydrogen bond interactions due to packing/structural restrictions and interfacial water molecules provide an extension to side-chains to accomplish hydrogen bonding.

Intrinsic bending and deformability at the T-A step of CCTTTAAAGG: a comparative analysis of T-A and A-T steps within A-tracts

Introduction of a T-A or pyrimidine-purine step into a straight and rigid A-tract can cause a positive roll deformation that kinks the DNA helix at that step. In CCTTTAAAGG, the central T-A step has an 8.6 degrees bend toward the major groove. We report the structural analysis of CCTTTAAAGG and a comparison with 25 other representative crystal structures from the NDB containing at least four consecutive A or T bases. On average, more local bending occurs at the disruptive T-A step (8.21 degrees ) than at an A-T step (5.71 degrees ). In addition, A-tracts containing an A-T step are more bent than are pure A-tracts, and hence A-A and A-T steps are not equivalent. All T-A steps examined exhibit positive roll, bending towards the major groove, while A-T steps display negative roll and bend slightly towards the minor groove. This illustrates how inherent negative and positive roll are, respectively, at A-T and T-A steps within A-tracts. T-A steps are more deformable, showing larger and more variable deformations of minor groove width, rise, cup, twist, and buckle. Standard deviations of twist, rise, and cup for T-A steps are 6.66 degrees, 0.55 A, and 15.90 degrees, versus 2.28 degrees, 0.21 A, and 2.99 degrees for A-T steps. Packing constraints determine which local values of these helical parameters an individual T-A step will adopt. For instance, with CCTTTAAAGG and three isomorphous structures, CGATTAATCG, CGATATATCG, and CGATCGATCG, crystal packing forces lead to a series of correlated changes: widened minor groove, large slide, low twist, and large rise. The difference in helical parameters between A-T steps lying within A-tracts, versus A-T steps within alternating AT sequences, demonstrates the importance of neighboring steps on the conformation of a given dinucleotide step.

Binding of bis-linked netropsin derivatives in the parallel-stranded hairpin form to DNA

Cis-diammine Pt(II)- bridged bis-netropsin and oligomethylene-bridged bis-netropsin in which two monomers are linked in a tail-to-tail manner bind to the DNA oligomer with the sequence 5'-CCTATATCC-3' in a parallel-stranded hairpin form with a stoichiometry 1:1. The difference circular dichroism (CD) spectra characteristic of binding of these ligands in the hairpin form are similar. They differ from CD patterns obtained for binding to the same duplex of another bis-netropsin in which two netropsin moieties were linked in a head-to-tail manner. This reflects the fact that tail-to-tail and head-to-tail bis-netropsins use parallel and antiparallel side-by-side motifs, respectively, for binding to DNA in the hairpin forms. The binding affinity of cis-diammine Pt(II)-bridged bis-netropsin in the hairpin form to DNA oligomers with nucleotide sequences 5'-CCTATATCC-3' (I), 5'-CCTTAATCC-3' (II), 5'-CCTTATTCC-3' (III), 5'-CCTTTTTCC-3' (IV) and 5'-CCAATTTCC-3' (V) decreases in the order I = II > III > IV > V . The binding of oligomethylene-bridged bis-netropsin in the hairpin form follows a similar hierarchy. An opposite order of sequence preferences is observed for partially bonded monodentate binding mode of the synthetic ligand.

Direct and indirect readout in mutant Met repressor-operator complexes

BACKGROUND

The methionine repressor, MetJ, represses the transcription of genes involved in methionine biosynthesis by binding to arrays of two to five adjacent copies of an eight base-pair 'metbox' sequence. Naturally occurring operators differ from the consensus sequence to a greater extent as the number of metboxes increases. MetJ, while accommodating this sequence variation in natural operators, is very sensitive to particular base changes, even where bases are not directly contacted in the crystal structure of a complex formed between the repressor and consensus operator.

RESULTS

Here we report the high-resolution structure of a MetJ mutant, Q44K, bound to the consensus operator sequence (Q44Kwt19) and two related sequences containing mutations at sites believed to be important for indirect readout at non-contacted bases. The overall structure of the Q44Kwt19 complex is very similar to the wild-type complex, but there are small variations in sugar-phosphate backbone conformation and direct contacts to the DNA bases. The mutant complexes show a mixture of direct and indirect readout of sequence variations, with differences in direct contacts and DNA conformation.

CONCLUSIONS

Comparison of the wild-type and mutant repressor-operator complexes shows that the repressor makes sufficiently strong interactions with the sugar-phosphate backbone to accommodate some variation in operator sequence with minor changes in direct bases contacts. The reduction in repressor affinity for the two mutant repressor complexes can be partially attributed to a loss in direct contacts to the DNA. In one case, however, the replacement of a flexible TA base-step leads to an unfavourable DNA conformation that reduces the stability of the repressor-operator complex.

Sequence-dependent conformational perturbation in DNA duplexes containing an epsilonA.T mismatch using molecular dynamics simulation

Previous experiments from this laboratory showed that 1, N:(6)-ethenoadenine (epsilonA) in 15mer DNA oligonucleotide duplexes with GGepsilonAGG and CCepsilonACC central sequences is repaired 3-5-fold more efficiently than in duplexes containing AAepsilonAAA and TTepsilonATT central sequences. This sequence dependence in repair rates appeared to correlate with the observed thermodynamic stability of these duplexes [Hang et al. (1998) J. Biol. Chem., 273, 33406-33413]. In the present work, unrestrained molecular dynamics was used to evaluate the sequence-dependent structural features of these duplexes. Explicit solvent and the particle mesh Ewald method were applied for the accurate representation of the electrostatic interactions. The differences observed in the axis- and intra-base pair parameters were primarily localized at the epsilonA*T mismatch in all sequences and indicate conformational diversity between the structures. However, all four structures remained in the B-conformational family. In the tip, tilt and propeller twist parameters for the five central base pairs, larger perturbations were found for the two duplexes with epsilonA flanked by A or T bases than for duplexes with epsilonA flanked by G or C bases. As a result of these perturbations, the average global curvature of the AAepsilonAAA and TTvarepsilonATT DNA duplexes was larger by approximately 12 degrees than that of the duplexes with the GGepsilonAGG and CCepsilonACC central sequences. The observed conformational differences between the duplexes containing A or T and G or C neighbors of epsilonA may contribute to the observed differential enzymatic repair of the same sequences.

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.

Crystal studies of B-DNA: the answers and the questions

The history of DNA crystallography is reviewed and is followed by discussion of the methods used for structure determination. The features of B-DNA molecular and crystal structures are described.

Structure of Poly (U).poly (A).poly (U)

The molecular structure of poly (U).poly (A).poly (U) has been determined and refined using the continuous x-ray intensity data on layer lines in the diffraction pattern obtained from an oriented fiber of the RNA. The final R-value for the preferred structure is 0.24, far lower than that for the plausible alternatives. The polymer forms an 11-fold right-handed triple-helix of pitch 33.5A and each base triplet is stabilized by Crick-Watson-Hoogsteen hydrogen bonds. The ribose rings in the three strands have C3'-endo, C2'-endo and C2'-endo conformations, respectively. The helix derives additional stability through systematic interchain hydrogen bonds involving ribose hydroxyls and uracil bases. The relatively grooveless cylindrical shape of the triple-helix is consistent with the lack of lateral organization.

Determination of metal ion binding sites within the hairpin ribozyme domains by NMR

Cations play an important role in RNA folding and stabilization. The hairpin ribozyme is a small catalytic RNA consisting of two domains, A and B, which interact in the transition state in an ion-dependent fashion. Here we describe the interaction of mono-, di-, and trivalent cations with the domains of the ribozyme, as studied by homo- and heteronuclear NMR spectroscopy. Paramagnetic line broadening, chemical shift mapping, and intermolecular NOEs indicate that the B domain contains four to five metal binding sites, which bind Mn(2+), Mg(2+), and Co(NH(3))(6)(3+). There is no significant structural change in the B domain upon the addition of Co(NH(3))(6)(3+) or Mg(2+). No specific monovalent ion binding sites exist on the B domain, as determined by (15)NH(4)(+) binding studies. In contrast to the B domain, there are no observable metal ion interactions within the internal loop of the A domain. Model structure calculations of Mn(2+) interactions at two sites within the B domain indicate that the binding sites comprise major groove pockets lined with functional groups oriented so that multiple hydrogen bonds can be formed between the RNA and Mn(H(2)O)(6)(2+) or Co(NH(3))(6)(3+). Site 1 is very similar in geometry to a site within the P4-P6 domain of the Tetrahymena group I intron, while site 2 is unique among known ion binding sites. The site 2 ion interacts with a catalytically essential nucleotide and bridges two phosphates. Due to its location and geometry, this ion may play an important role in the docking of the A and B domains.

Hydration of [d(CGC)r(aaa)d(TTTGCG)](2)

We have studied the hydration and dynamics of RNA C2'-OH in a DNA. RNA hybrid chimeric duplex [d(CGC)r(aaa)d(TTTGCG)](2). Long-lived water molecules with correlation time tau(c) larger than 0.3 ns were found close to the RNA adenine H2 and H1' protons in the hybrid segment. A possible long-lived water molecule was also detected close to the methyl group of 7T in the RNA-DNA junction but not to the other two thymine bases (8T and 9T). This result correlates with the structural studies that only DNA residue 7T in the RNA-DNA junction adopts an O4'-endo sugar conformation (intermediate between B-form and A-form), while the other DNA residues including 3C in the DNA-RNA junction, adopt C1'-exo or C2'-endo conformations (in the B-form domain). Based on the NOE cross-peak patterns, we have found that RNA C2'-OH tends to orient toward the O3' direction, forming a possible hydrogen bond with the 3'-phosphate group. The exchange rates for RNA C2'-OH were found to be around 5-20 s(-1), compared to 26.7(+/-13.8) s(-1) reported previously for the other DNA.RNA hybrid duplex. This slow exchange rate may be due to the narrow minor groove width of [d(CGC)r(aaa)d(TTTGCG)](2), which may trap the water molecules and restrict the dynamic motion of hydroxyl protons. The distinct hydration patterns of the RNA adenine H2 and H1' protons and the DNA 7T methyl group in the hybrid segment, as well as the orientation and dynamics of the RNA C2'-OH protons, may provide a molecular basis for further understanding the structure and recognition of DNA.RNA hybrid and chimeric duplexes.

Water and ions in a high resolution structure of B-DNA

A detailed picture of hydration and counterion location in the B-DNA duplex d(GCGAATTCG) is presented. Detailed data have been obtained by single crystal x-ray diffraction at atomic resolution (0.89 A) in the presence of Mg(2+). The latter is the highest resolution ever obtained for a B-DNA oligonucleotide. Minor groove hydration is compared with that found in the Na(+) and Ca(2+) crystal forms of the related dodecamer d(CGCGAATTCGCG). High resolution data (1.45 A) of the Ca(2+) form obtained in our laboratory are used for that purpose. The central GAATTC has a very stable hydration spine identical in all cases, independent of duplex length and crystallization conditions (counterions, space group). However, the organization of the water molecules (tertiary and quaternary layers) associated with the central spine vary in each case.

Modeling high-resolution hydration patterns in correlation with DNA sequence and conformation

Hydration around the DNA fragment d(C5T5).(A5G5) is presented from two molecular dynamics simulations of 10 and 12 ns total simulation time. The DNA has been simulated as a flexible molecule with both the CHARMM and AMBER force fields in explicit solvent including counterions and 0.8 M additional NaCl salt. From the previous analysis of the DNA structure B-DNA conformations were found with the AMBER force-field and A-DNA conformations with CHARMM parameters. High-resolution hydration patterns are compared between the two conformations and between C.G and T.A base-pairs from the homopolymeric parts of the simulated sequence. Crystallographic results from a statistical analysis of hydration sites around DNA crystal structures compare very well with the simulation results. Differences between the crystal sites and our data are explained by variations in conformation, sequence, and limitations in the resolution of water sites by crystal diffraction. Hydration layers are defined from radial distribution functions and compared with experimental results. Excellent agreement is found when the measured experimental quantities are compared with the equivalent distribution of water molecules in the first hydration shell. The number of water molecules bound to DNA was found smaller around T.A base-pairs and around A-DNA as compared to B-DNA. This is partially offset by a larger number of water molecules in hydrophobic contact with DNA around T.A base-pairs and around A-DNA. The numbers of water molecules in minor and major grooves have been correlated with helical roll, twist, and inclination angles. The data more fully explain the observed B-->A transition at low humidity.

Targeting of human retrotransposon integration is directed by the specificity of the L1 endonuclease for regions of unusual DNA structure

L1 elements are polyA retrotransposons which inhabit the human genome. Recent work has defined an endonuclease (L1 EN) encoded by the L1 element required for retrotransposition. We report the sequence specificity of this nicking endonuclease and the physical basis of its DNA recognition. L1 endonuclease is specific for the unusual DNA structural features found at the TpA junction of 5'(dTn-dAn) x 5'(dTn-dAn) tracts. Within the context of this sequence, substitutions which generate a pyrimidine-purine junction are tolerated, whereas purine-pyrimidine junctions greatly reduce or eliminate nicking activity. The A-tract conformation of the DNA substrate 5' of the nicked site is required for L1 EN nicking. Chemical or physical unwinding of the DNA helix enhances L1 endonuclease activity, while disruption of the adenine mobility associated with TpA junctions reduces it. Akin to the protein-DNA interactions of DNase I, L1 endonuclease DNA recognition is likely mediated by minor groove interactions. Unlike several of its homologues, however, L1 EN exhibits no AP endonuclease activity. Finally, we speculate on the implications of the specificity of the L1 endonuclease for the parasitic relationship between retroelements and the human genome.

The residence time of the bound water in the hydrophobic minor groove of the carbocyclic-nucleoside analogs of Dickerson-Drew dodecamers

The residence time of the bound water molecules in the antisense oligodeoxyribonucleotides containing 7'-alpha-methyl (TMe) carbocyclic thymidines in duplex (I), d5'(1C2G3C4G5A6A7TMc8TMe9C10G11C12G)23', and 6'-alpha-hydroxy (TOH) carbocyclic thymidines in duplex (II), d5'(1C2G3C4G5AOH6AOH7TOH8 TOH9C10G11C12G)2(3), have been investigated using a combination of NOESY and ROESY experiments. Because of the presence of 7'-alpha-methyl groups of TMe in the centre of the minor groove in duplex (I), the residence time of the bound water molecule is shorter than 0.3 ns. The dramatic reduction of the residence time of the water molecule in the minor groove in duplex (I) compared with the natural counterpart has been attributed to the replacement of second shell of hydration and disruption of hydrogen-bonding with 04' in the minor groove by hydrophobic alpha-methyl groups, as originally observed in the X-ray study. This effect could not be attributed to the change of the width of the minor groove because a comparative NMR study of the duplex (I) and its natural counterpart showed that the widths of their minor grooves are more or less unchanged (r.m.s.d change in the core part is <0.63A). For duplex (II) with polar 6'-alpha-hydroxyl groups pointed to the minor groove, the correlation time is much longer than 0.36 ns as a result of the stabilising hydrogen-bonding interaction with N3 or 02 of the neighbouring nucleotides.

An NMR study of d(CTACTGCTTTAG).d(CTAAAGCAGTAG) showing hydration water molecules in the minor groove of a TpA step

The hydration properties of the non-palindromic duplex d(CTACTGCTTTAG). d(CTAAAGCAGTAG) were investigated by NMR spectroscopy. The oligonucleotide possesses a heterogeneous B-DNA structure. The H2(n)-H1'(m+1) distances reflect a minor groove narrowing within the TTT/AAA segment (approximately 3.9A) and a sudden widening at the T10:A15 base-pair (approximately 5.3A), the standard B-DNA distance being approximately 5A. The facing T10pA11 and T14pA15 steps at the end of the TTTA/AAAT segment have completely different behaviors. Only A15 ending the AAA run displays NMR features comparable to those shown by adenines of TpA steps occupying the central position of TnAn (n> or =2) segments. These involve particular chemical shifts and line broadening of the H2 and H8 protons. Positive NOESY cross-peaks were measured between the water protons and the H2 protons of A15, A16 and A17 reflecting the occurrence of hydration water molecules with residence times longer than 500 picoseconds along the minor groove of the TTT/AAA segment. In contrast no water molecules with long residence times were observed neither for A3, A20 and A23 nor for A11 ending the 5'TTTA run. We confirm thus that the binding of water molecules with long residence time to adenine residues correlates with the minor groove narrowing. In contrast, the widening of the minor groove at the A11:T14 base-pair ending the TTTA/TAAA segment, likely associated to a high negative propeller twist value at this base-pair, prevents the binding of a water molecule with long residence time to A11 but not to A15 of the preceding T10:A15 base-pair. Thus, in our non-palindromic oligonucleotide the water molecules bind differently to A11 and A15 although both adenines are part of a TpA step. The slower motions occurring at A15 compared to A11 are also well explained by the present results.

DNA sequence recognition by bis-linked netropsin and distamycin derivatives

We studied the interaction of cis-diammine Pt(II)-bridged bis-netropsin, cis-diammine Pt(II)-bridged bis-distamycin and oligomethylene-bridged bis-netropsin with synthetic DNA fragments containing pseudosymmetrical AT-rich nucleotide sequences and compared it with the interaction of the parent compounds netropsin and distamycin A. For fragments containing multiple blocks of (AIT)4 and (T/A)4 separated by zero, one, two and three GC-base pairs, DNase I footprinting and CD spectroscopy studies reveal that 5'-TTTTAAAA-3' is the strongest affinity binding site for cis-diammine Pt(II)-bridged bis-netropsin and bis-distamycin. They both bind less strongly to a DNA region containing the sequence 5'-AAAATTTT-3'. Netropsin, distamycin A and oligomethylene-bridged bis-netropsin exhibit far less sequence discrimination.

Hydration of the phosphate group in double-helical DNA

Water distributions around phosphate groups in 59 B-, A-, and Z-DNA crystal structures were analyzed. It is shown that the waters are concentrated in six hydration sites per phosphate and that the positions and occupancies of these sites are dependent on the conformation and type of nucleotide. The patterns of hydration that are characteristic of the backbone of the three DNA helical types can be attributed in part to the interactions of these hydration sites.

The N-terminal methionine is a major determinant of the DNA binding specificity of MEF-2C

Members of the MEF-2 family of transcriptional regulators positively modulate the activity of basic helix-loop-helix proteins in both myogenic and neurogenic cell lineages. Previous work had shown that MEF-2C(2-117), a protein fragment comprising the dimerization and DNA-binding domains of MEF-2C but lacking the N-terminal methionine, bound to AT-rich DNA sequences with high affinity. MEF-2C(2-117) did not discriminate between different AT-rich sequences. We now report the in vitro DNA binding properties of a MEF-2C fragment containing the N-terminal methionine. Measurements of the apparent dissociation constants of the complexes of GG-MEF-2C(1-117) revealed that different AT-rich sequences are bound with different affinities; in particular MEF site containing DNA (CTATAAATAG) is bound preferentially to DNA containing a SRF site (CATAAATG). Strikingly, when the shorter AT run consisted of six alternating thymines and adenines, almost wild-type affinity was observed. Irrespective of the particular DNA sequence, all circular dichroism spectra of the DNA complexes of GG-MEF-2C(1-117) were superimposable and characterized by an identical maximal ellipticity at 269.5 nm, suggesting similar DNA conformations. Bending analysis by circular permutation assay revealed that on complex formation MEF-2C(2-117) induced cognate DNA to bend by 49 degrees, while heterologous DNA remained unbent. In the presence of the N-terminal methionine, however, all DNA sequences were bent by 70 degrees. The above results suggest an important function for the N-terminal methionine in properly orientating MEF-2C on the DNA.