CAP binding sites reveal pyrimidine-purine pattern characteristic of DNA bending

To investigate the intrinsic bending of DNA at sites where proteins bind, we analyzed catabolite gene activator protein (CAP) binding sites and various operators from the viewpoint of DNA bending flexibility. Theoretical conformational analysis. DNase I digestion and x-ray crystallography data indicate that bending of B-DNA is highly anisotropic and sequence-dependent. Certain dimers prefer to bend into the major groove ("major-philic") and others prefer to bend into the minor groove ("minor-philic" dimers). From these data we considered TA, CG, CA:TG and GG:CC as major-philic dimers and AT,AA:TT and GT:AC as minor-philic ones. Analysis of 31 CAP binding sites has identified strong major-philic tendencies 5-7 base pairs (bp) away from the center. In addition, we found minor-philic poly-A tracts extending 4-5 bp away from the proposed major-philic bends. Finally, to analyze the central regions we followed the lead of Shumilov and classified the DNA sites by their spacer lengths [V.Y. Shumilov, Mol. Biol. (Mosk) 21, 168-187 (1987)]. In this way, we identified two subsets of CAP binding sites: one with 6 bp between the TGTGA:TCACA consensus boxes (N6-set) and one with 8 central bp (N8-set). We discovered that the dimer at the center of an N6-set site was usually major-philic, whereas at the center of an N8-set site more often minor-philic. Analysis of phages 434, P22 lambda and trp operators revealed similar results. In conclusion, our data show that CAP binding sites have major-philic and minor-philic dimers at specific positions; the location of these dimers may facilitate wrapping of DNA around CAP. A similar pattern is seen in nucleosomes.     

Atomic-resolution crystal structures of B-DNA reveal specific influences of divalent metal ions on conformation and packing.

Crystal structures of B-form DNA have provided insights into the global and local conformational properties of the double helix, the solvent environment, drug binding and DNA packing. For example, structures of the duplex with sequence CGCGAATTCGCG, the Dickerson-Drew dodecamer (DDD), established a unique geometry of the central A-tract and a hydration spine in the minor groove. However, our knowledge of the various interaction modes between metal ions and DNA is very limited and almost no information exists concerning the origins of the different effects on DNA conformation and packing exerted by individual metal ions.Crystallization of the DDD duplex in the presence of Mg(2+)and Ca(2+)yields different crystal forms. The structures of the new Ca(2+)-form and isomorphous structures of oligonucleotides with sequences GGCGAATTCGCG and GCGAATTCGCG were determined at a maximum resolution of 1.3 A. These and the 1.1 A structure of the DDD Mg(2+)-form have revealed the most detailed picture yet of the ionic environment of B-DNA. In the Mg(2+)and Ca(2+)-forms, duplexes in the crystal lattice are surrounded by 13 magnesium and 11 calcium ions, respectively.Mg(2+)and Ca(2+)generate different DNA crystal lattices and stabilize different end-to-end overlaps and lateral contacts between duplexes, thus using different strategies for reducing the effective repeat length of the helix to ten base-pairs. Mg(2+)crystals allow the two outermost base-pairs at either end to interact laterally via minor groove H-bonds, turning the 12-mer into an effective 10-mer. Ca(2+)crystals, in contrast, unpair the outermost base-pair at each end, converting the helix into a 10-mer that can stack along its axis. This reduction of a 12-mer into a functional 10-mer is followed no matter what the detailed nature of the 5'-end of the chain: C-G-C-G-A-ellipsis, G-G-C-G-A-ellipsis, or a truncated G-C-G-A-ellipsis Rather than merely mediating close contacts between phosphate groups, ions are at the origin of many well-known features of the DDD duplex structure. A Mg(2+)coordinates in the major groove, contributing to kinking of the duplex at one end. While Ca(2+)resides in the minor groove, coordinating to bases via its hydration shell, two magnesium ions are located at the periphery of the minor groove, bridging phosphate groups from opposite strands and contracting the groove at one border of the A-tract.     

Efficient pH-independent sequence-specific DNA binding by pseudoisocytosine-containing bis-PNA.

The synthesis and DNA binding properties of bis-PNA (peptide nucleic acid) are reported. Two PNA segments each of seven nucleobases in length were connected in a continuous synthesis via a flexible linker composed of three 8-amino-3,6-dioxaoctanoic acid units. The sequence of the first strand was TCTCTTT (C- to N-terminal), while the second strand was TTTCTCT or TTTJTJT, where J is pseudoisocytosine. These bis-PNAs form triple-stranded complexes of somewhat higher thermal stability than monomeric PNA with complementary oligonucleotides and the thermal melting transition shows very little hysteresis. When the J base is placed in the strand parallel to the DNA complement ('Hoogsteen strand'), the DNA binding was pH independent. The bis-PNAs were also superior to monomeric PNAs for targeting double-stranded DNA by strand invasion.     

High resolution crystal structures of piscine transthyretin reveal different binding modes for triiodothyronine and thyroxine

Transthyretin (TTR) is an extracellular transport protein involved in the distribution of thyroid hormones and vitamin A. So far, TTR has only been found in vertebrates, of which piscine TTR displays the lowest sequence identity with human TTR (47%). Human and piscine TTR bind both thyroid hormones 3,5,3'-triiodo-l-thyronine (T(3)) and 3,5,3',5'-tetraiodo-l-thyronine (thyroxine, T(4)). Human TTR has higher affinity for T(4) than T(3), whereas the reverse holds for piscine TTR. X-ray structures of Sparus aurata (sea bream) TTR have been determined as the apo-protein at 1.75 A resolution and bound to ligands T(3) and T(4), both at 1.9 A resolution. The apo structure is similar to human TTR with structural changes only at beta-strand D. This strand forms an extended loop conformation similar to the one in chicken TTR. The piscine TTR.T(4) complex shows the T(4)-binding site to be similar but not identical to human TTR, whereas the TTR.T(3) complex shows the I3' halogen situated at the site normally occupied by the hydroxyl group of T(4). The significantly wider entrance of the hormone-binding channel in sea bream TTR, in combination with its narrower cavity, provides a structural explanation for the different binding affinities of human and piscine TTR to T(3) and T(4).     

Thermodynamics of sequence-specific binding of PNA to DNA

For further characterization of the hybridization properties of peptide nucleic acids (PNAs), the thermodynamics of hybridization of mixed sequence PNA-DNA duplexes have been studied. We have characterized the binding of PNA to DNA in terms of binding affinity (perfectly matched duplexes) and sequence specificity of binding (singly mismatched duplexes) using mainly absorption hypochromicity melting curves and isothermal titration calorimetry. For perfectly sequence-matched duplexes of varying lengths (6-20 bp), the average free energy of binding (DeltaG degrees ) was determined to be -6.5+/-0.3 kJ mol(-1) bp(-1), corresponding to a microscopic binding constant of about 14 M(-1) bp(-1). A variety of single mismatches were introduced in 9- and 12-mer PNA-DNA duplexes. Melting temperatures (T(m)) of 9- and 12-mer PNA-DNA duplexes with a single mismatch dropped typically 15-20 degrees C relative to that of the perfectly matched sequence with a corresponding free energy penalty of about 15 kJ mol(-1) bp(-1). The average cost of a single mismatch is therefore estimated to be on the order of or larger than the gain of two matched base pairs, resulting in an apparent binding constant of only 0.02 M(-1) per mismatch. The impact of a mismatch was found to be dependent on the neighboring base pairs. To a first approximation, increasing the stability of the surrounding region, i.e., the distribution of A.T and G.C base pairs, decreases the effect of the introduced mismatch.     

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.     

Structural features of B-DNA dodecamer crystal structures: influence of crystal packing versus base sequence

We have analyzed the set of inter and intra base pair parameters for each dinucleotide step in single crystal structures of dodecamers, solved at high and medium resolution and all crystallized in P2(1)2(1)2(1) space group. The objective was to identify whether all the structures which have either the Drew-Dickerson (DD) sequence d[CGCGAATTCGCG] with some base modification or related sequence (non-DD), would display the same sequence dependent structural variability about its palindromic sequence, despite the molecule being bent at one end because of similar crystal lattice packing effect. Most of the local doublet parameters for base pairs steps G2-C3 and G10-C11 positions, symmetrically situated about the lateral two-fold, were significantly correlated between themselves. In non-DD sequences, significant correlations between these positional parameters were absent. The different range of local step parameter values at each sequence position contributed to the gross feature of smooth helix axis bending in all structures. The base pair parameters in some of the positions, for medium resolution DD sequence, were quite unlike the high-resolution set and encompassed a higher range of values. Twist and slide are the two main parameters that show wider conformational range for the middle region of non-DD sequence structures in comparison to DD sequence structures. On the contrary, the minor and major groove features bear good resemblance between DD and non-DD sequence crystal structure datasets. The sugar-phosphate backbone torsion angles are similar in all structures, in sharp contrast to base pair parameter variation for high and low resolution DD and non-DD sequence structures, consisting of unusual (epsilon = g-, xi = t) BII conformation at the 10th position of the dodecamer sequence. Thus examining DD and non-DD sequence structures packed in the same crystal lattice arrangement, we infer that inter and intra base pair parameters are as symmetrically equivalent in its value as the symmetry related step for the palindromic DD sequence about lateral two-fold axis. This feature would lead us to agree with the conclusion that DNA conformation is not substantially affected by end-to-end or lateral inter-molecular interaction due to crystal lattice packing effect. Non-DD sequence structures acquire step parameter values which reflect the altered sequence at each of the dodecamer sequence position in the orthorhombic lattice while showing similar gross features of DD sequence structures.     

Zfy1 encodes a nuclear sequence-specific DNA binding protein

Zfy1 is a mouse Y chromosomal gene encoding a zincfinger protein which is thought to have some function during spermatogenesis. Here we show that, when introduced into tissue culture cells, Zfy1 is targeted to the nucleus. Two independent signals are present within the protein for nuclear localization. This nuclear Zfy1 protein is able to bind strongly to DNA-cellulose and, using site-selection assays, we have identified specific Zfy1 DNA binding sites. Taken together these results suggest that Zfy1 is a nuclear-located sequence-specific DNA binding protein which functions during spermatogenesis.     

DNA recognition by F factor TraI36: highly sequence-specific binding of single-stranded DNA

The TraI protein has two essential roles in transfer of conjugative plasmid F Factor. As part of a complex of DNA-binding proteins, TraI introduces a site- and strand-specific nick at the plasmid origin of transfer (oriT), cutting the DNA strand that is transferred to the recipient cell. TraI also acts as a helicase, presumably unwinding the plasmid strands prior to transfer. As an essential feature of its nicking activity, TraI is capable of binding and cleaving single-stranded DNA oligonucleotides containing an oriT sequence. The specificity of TraI DNA recognition was examined by measuring the binding of oriT oligonucleotide variants to TraI36, a 36-kD amino-terminal domain of TraI that retains the sequence-specific nucleolytic activity. TraI36 recognition is highly sequence-specific for an 11-base region of oriT, with single base changes reducing affinity by as much as 8000-fold. The binding data correlate with plasmid mobilization efficiencies: plasmids containing sequences bound with lower affinities by TraI36 are transferred between cells at reduced frequencies. In addition to the requirement for high affinity binding to oriT, efficient in vitro nicking and in vivo plasmid mobilization requires a pyrimidine immediately 5' of the nick site. The high sequence specificity of TraI single-stranded DNA recognition suggests that despite its recognition of single-stranded DNA, TraI is capable of playing a major regulatory role in initiation and/or termination of plasmid transfer.     

A geminivirus replication protein is a sequence-specific DNA binding protein.

The genome of the geminivirus tomato golden mosaic virus (TGMV) consists of two circular DNA molecules designated as components A and B. The A component encodes the only viral protein, AL1, that is required for viral replication. We showed that AL1 interacts specifically with TGMV A and B DNA by using an immunoprecipitation assay for AL1:DNA complex formation. In this assay, a monoclonal antibody against AL1 precipitated AL1:TGMV DNA complexes, whereas an unrelated antibody failed to precipitate the complexes. Competition assays with homologous and heterologous DNAs established the specificity of AL1:DNA binding. AL1 produced by transgenic tobacco plants and by baculovirus-infected insect cells exhibited similar DNA binding activity. The AL1 binding site maps to 52 bp on the left side of the common region, a 235-bp region that is highly conserved between the two TGMV genome components. The AL1:DNA binding site does not include the putative hairpin structure that is conserved in the common regions or the equivalent 5' intergenic regions of all geminiviruses. These studies demonstrate that a geminivirus replication protein is a sequence-specific DNA binding protein, and the studies have important implications for the role of this protein in virus replication.     

High resolution crystal structures and molecular dynamics studies reveal substrate binding in the porin Omp32

The porin Omp32 is the major outer membrane protein of the bacterium Delftia acidovorans. The crystal structures of the strongly anion-selective porin alone and in complex with the substrate malate were solved at 1.5 and 1.45 A resolution, respectively, and revealed a malate-binding motif adjacent to the channel constriction zone. Binding is mediated by interaction with a cluster of two arginine residues and two threonines. This binding site is specific for Omp32 and reflects the physiological adaptation of the organism to organic acids. Structural studies are combined with a 7-ns unbiased molecular dynamics simulation of the trimeric channel in a model membrane. Molecular dynamics trajectories show how malate ions are efficiently captured from the surrounding bulk solution by the electrostatic potential of the channel, translocated to the binding site region, and immobilized in the constriction zone. In accordance with these results, conductance measurements with Omp32 inserted in planar lipid membranes revealed binding of malate. The anion-selective channel Omp32 is the first reported example of a porin with a 16-stranded beta-barrel and proven substrate specificity. This finding suggests a new view on the correlation of porin structure with substrate binding in specific channels.     

Isolation and sequence-specific DNA binding of the Antennapedia homeodomain.

The homeodomain encoded by the Antennapedia (Antp) gene of Drosophila was overproduced in a T7 expression vector in Escherichia coli. The corresponding polypeptide of 68 amino acids was purified to homogeneity. The homeodomain was analysed by ultracentrifugation and assayed for DNA binding. The secondary structure of the isolated homeodomain was determined by nuclear magnetic resonance spectroscopy. DNA-binding studies indicate that the isolated homeodomain binds to DNA in vitro. It selectively binds to the same sites as a longer Antp polypeptide and a full-length fushi tarazu (ftz) protein. Therefore, the homeodomain represents the DNA-binding domain of the homeotic proteins.     

Molecular dynamics simulations of A. T-rich oligomers: sequence-specific binding of Na+ in the minor groove of B-DNA

Molecular dynamics simulations have been employed to probe the sequence-specific binding of sodium ions to the minor groove of B-DNA of three A. T-rich oligomers having identical compositions but different orders of the base pairs: C(AT)(4)G, CA(4)T(4)G, and CT(4)A(4)G. Recent experimental investigations, either in crystals or in solution, have shown that monovalent cations bind to DNA in a sequence-specific mode, preferentially in the narrow minor groove regions of uninterrupted sequences of four or more adenines (A-tracts), replacing a water molecule of the ordered hydration structure, the hydration spine. Following this evidence, it has been hypothesized that in A-tracts these events may be responsible for structural peculiarities such as a narrow minor groove and a curvature of the helix axis. The present simulations confirm a sequence specificity of the binding of sodium ions: Na(+) intrusions in the first layer of hydration of the minor groove, with long residence times, up to approximately 3 ns, are observed only in the minor groove of A-tracts but not in the alternating sequence. The effects of these intrusions on the structure of DNA depend on the ion coordination: when the ion replaces a water molecule of the spine, the minor groove becomes narrower. Ion intrusions may also disrupt the hydration spine modifying the oligomer structure to a large extent. However, in no case intrusions were observed to locally bend the axis toward the minor groove. The simulations also show that ions may reside for long time periods in the second layer of hydration, particularly in the wider regions of the groove, often leading to an opening of the groove.