DNA conformation and configuration in protein-DNA complexes

Conformational changes in DNA-binding proteins: relationships with precomplex features and contributions to specificity and stability

Both Proteins and DNA undergo conformational changes in order to form functional complexes and also to facilitate interactions with other molecules. These changes have direct implications for the stability and specificity of the complex, as well as the cooperativity of interactions between multiple entities. In this work, we have extensively analyzed conformational changes in DNA-binding proteins by superimposing DNA-bound and unbound pairs of protein structures in a curated database of 90 proteins. We manually examined each of these pairs, unified the authors' annotations, and summarized our observations by classifying conformational changes into six structural categories. We explored a relationship between conformational changes and functional classes, binding motifs, target specificity, biophysical features of unbound proteins, and stability of the complex. In addition, we have also investigated the degree to which the intrinsic flexibility can explain conformational changes in a subset of 52 proteins with high quality coordinate data. Our results indicate that conformational changes in DNA-binding proteins contribute significantly to both the stability of the complex and the specificity of targets recognized by them. We also conclude that most conformational changes occur in proteins interacting with specific DNA targets, even though unbound protein structures may have sufficient information to interact with DNA in a nonspecific manner.

Tautomeric equilibrium of uracil and thymine in model protein-nucleic acid contacts. Spectroscopic and quantum chemical approach

This work deals with tautomeric transformations of uracil (Ura) and thymine (Thy) in their model complexes with the deprotonated carboxylic group. Essential changes in the UV spectra of the bases upon their interaction with NaAc, vanishing signals of both imino protons in (1)H NMR spectra, and a perceptible decrease in intensity of both IR bands, related to the stretching vibrations nu(C=O) of the carbonyl groups, imply involvement of enolic tautomers. Results of quantum chemical calculations of the double complexes of the Ura(Thy) tautomers with CH(3)COO(-) at the MP2/6-311++G(2df,pd)//B3LYP/6-311++G(d,p) level of theory proved to be incompatible with the spectral features: despite the fact that the complexes of the enolic tautomers are much closer in energy to the diketo ones as compared to isolated tautomers, the energy gap between them is such that in tautomeric equilibrium dominate diketo forms. Calculations of triple complexes of the type CH(3)COO(-):Ura(Thy) tautomer:Na(+), taking into account the effect of the Na(+) coordination with tautomers, show that three triple complexes formed by enolic tautomers appeared more stable than those formed by diketo ones. This makes the UV and (1)H NMR data understandable, but the high residual intensity of the nu(C=O) bands in the IR spectra remains unclear. At that ion, Na(+) itself was not able to disturb the tautomeric equilibrium in the coordination complexes of the type Ura(Thy) tautomer:Na(+). To evaluate the DMSO effect, the CPCM solvation model was applied to triple complexes of the Ura tautomers. It appeared that in the solution there is coexistence between the diketo and enolic tautomers in a ratio of 53%:47%. This makes possible reconciliation of our experimental data. The biological significance of high-energy tautomers of nucleotide bases is discussed.

Directional shape complementarity at the protein-DNA interface

Nature utilizes various styles of architecture for DNA-binding proteins to recognize diverse DNA sequences, a process facilitated by a complementary surface between protein and DNA. However, the extent and ways this 'shape complementarity' occurs at the protein-DNA interface have yet to be characterized. Here, by analyzing a set of diverse protein-DNA complexes of known three-dimensional structures, we investigated whether the normal vectors of a protein surface at the interface exhibited any relationship with DNA conformation. Generally, the normal vectors of a DNA-contacting protein surface distinctly preferred certain angles, enabling them to align with certain axes characterizing the conformation of DNA. Thus, a new geometric property of DNA-binding protein is demonstrated, i.e. the "shape complementarity" of protein-DNA recognition clearly bears the property of "directionality".

A-form conformational motifs in ligand-bound DNA structures

Recognition and biochemical processing of DNA requires that proteins and other ligands are able to distinguish their DNA binding sites from other parts of the molecule. In addition to the direct recognition elements embedded in the linear sequence of bases (i.e. hydrogen bonding sites), these molecular agents seemingly sense and/or induce an "indirect" conformational response in the DNA base-pairs that facilitates close intermolecular fitting. As part of an effort to decipher this sequence-dependent structural code, we have analyzed the extent of B-->A conformational conversion at individual base-pair steps in protein and drug-bound DNA crystal complexes. We take advantage of a novel structural parameter, the position of the phosphorus atom in the dimer reference frame, as well as other documented measures of local helical structure, e.g. torsion angles, base-pair step parameters. Our analysis pinpoints ligand-induced conformational changes that are difficult to detect from the global perspective used in other studies of DNA structure. The collective data provide new structural details on the conformational pathway connecting A and B-form DNA and illustrate how both proteins and drugs take advantage of the intrinsic conformational mechanics of the double helix. Significantly, the base-pair steps which exhibit pure A-DNA conformations in the crystal complexes follow the scale of A-forming tendencies exhibited by synthetic oligonucleotides in solution and the known polymorphism of synthetic DNA fibers. Moreover, most crystallographic examples of complete B-to-A deformations occur in complexes of DNA with enzymes that perform cutting or sealing operations at the (O3'-P) phosphodiester linkage. The B-->A transformation selectively exposes sugar-phosphate atoms, such as the 3'-oxygen atom, ordinarily buried within the chain backbone for enzymatic attack. The forced remodeling of DNA to the A-form also provides a mechanism for smoothly bending the double helix, for controlling the widths of the major and minor grooves, and for accessing the minor groove edges of individual base-pairs.

Modeling helix-turn-helix protein-induced DNA bending with knowledge-based distance restraints

A crucial element of many gene functions is protein-induced DNA bending. Computer-generated models of such bending have generally been derived by using a presumed bending angle for DNA. Here we describe a knowledge-based docking strategy for modeling the structure of bent DNA recognized by a major groove-inserting alpha-helix of proteins with a helix-turn-helix (HTH) motif. The method encompasses a series of molecular mechanics and dynamics simulations and incorporates two experimentally derived distance restraints: one between the recognition helix and DNA, the other between respective sites of protein and DNA involved in chemical modification-enabled nuclease scissions. During simulation, a DNA initially placed at a distance was "steered" by these restraints to dock with the binding protein and bends. Three prototype systems of dimerized HTH DNA binding were examined: the catabolite gene activator protein (CAP), the phage 434 repressor (Rep), and the factor for inversion stimulation (Fis). For CAP-DNA and Rep-DNA, the root mean square differences between model and x-ray structures in nonhydrogen atoms of the DNA core domain were 2.5 A and 1.6 A, respectively. An experimental structure of Fis-DNA is not yet available, but the predicted asymmetrical bending and the bending angle agree with results from a recent biochemical analysis.

Protein-DNA interactions: A structural analysis

A detailed analysis of the DNA-binding sites of 26 proteins is presented using data from the Nucleic Acid Database (NDB) and the Protein Data Bank (PDB). Chemical and physical properties of the protein-DNA interface, such as polarity, size, shape, and packing, were analysed. The DNA-binding sites shared common features, comprising many discontinuous sequence segments forming hydrophilic surfaces capable of direct and water-mediated hydrogen bonds. These interface sites were compared to those of protein-protein binding sites, revealing them to be more polar, with many more intermolecular hydrogen bonds and buried water molecules than the protein-protein interface sites. By looking at the number and positioning of protein residue-DNA base interactions in a series of interaction footprints, three modes of DNA binding were identified (single-headed, double-headed and enveloping). Six of the eight enzymes in the data set bound in the enveloping mode, with the protein presenting a large interface area effectively wrapped around the DNA.A comparison of structural parameters of the DNA revealed that some values for the bound DNA (including twist, slide and roll) were intermediate of those observed for the unbound B-DNA and A-DNA. The distortion of bound DNA was evaluated by calculating a root-mean-square deviation on fitting to a canonical B-DNA structure. Major distortions were commonly caused by specific kinks in the DNA sequence, some resulting in the overall bending of the helix. The helix bending affected the dimensions of the grooves in the DNA, allowing the binding of protein elements that would otherwise be unable to make contact. From this structural analysis a preliminary set of rules that govern the bending of the DNA in protein-DNA complexes, are proposed.

RecA binding to a single double-stranded DNA molecule: a possible role of DNA conformational fluctuations

Most genetic regulatory mechanisms involve protein-DNA interactions. In these processes, the classical Watson-Crick DNA structure sometimes is distorted severely, which in turn enables the precise recognition of the specific sites by the protein. Despite its key importance, very little is known about such deformation processes. To address this general question, we have studied a model system, namely, RecA binding to double-stranded DNA. Results from micromanipulation experiments indicate that RecA binds strongly to stretched DNA; based on this observation, we propose that spontaneous thermal stretching fluctuations may play a role in the binding of RecA to DNA. This has fundamental implications for the protein-DNA binding mechanism, which must therefore rely in part on a combination of flexibility and thermal fluctuations of the DNA structure. We also show that this mechanism is sequence sensitive. Theoretical simulations support this interpretation of our experimental results, and it is argued that this is of broad relevance to DNA-protein interactions.

Thermodynamics of specific and non-specific DNA binding by the c-Myb DNA-binding domain

The thermodynamics of the c-Myb DNA-binding domain (R2R3) interaction with its target DNA have been analyzed using isothermal titration calorimetry and amino acid mutagenesis. The enthalpy of association between the standard R2R3, the Cys130 mutant substituted with Ile, and the cognate DNA is -12.5 (+/- 0.1) kcal mol-1 at pH 7.5 and at 20 degrees C, and this interaction is enthalpically driven throughout the physiological temperature range. In order to understand the DNA recognition mechanism, several pairs of interactions were investigated using single and multiple-base alterations with single and multiple-amino acid substituted mutants. The interactions between the standard R2R3 and many non-cognate DNAs were accompanied by binding enthalpy changes and heat capacity changes, although their affinities were reduced. The roles of the electrostatic interactions in binding to the cognate and the non-cognate DNAs were also analyzed from the dependency of the thermodynamic parameters on the salt concentration. The heat capacity change was found to be significantly dependent upon the salt concentration. Several mutant proteins bound to the multiple-base altered DNA with very small enthalpy changes, although they bound to the cognate and the single-base altered DNAs with detectable enthalpy and heat capacity changes. From the thermodynamic cycles derived from the DNA binding of the amino acid substituted R2R3 to the base substituted DNA duplexes, the individual thermodynamic mechanisms of the specific DNA recognition of R2R3 were dissected. The local folding mechanism was highlighted by the substitution of Pro with either Gly or Ala at the linker between R2 and R3. The characteristic thermodynamic features of specific and non-specific DNA binding are discussed.

DNA bending by GCN4 mutants bearing cationic residues

Transcription activation is thought to require DNA bending to promote the interaction of upstream activators and the basal transcription machinery. Previous experiments have shown that some members of the bZIP family of DNA binding proteins bend DNA, while others do not. We are exploring the possibility that electrostatic effects play a role in these differences. The yeast bZIP transcription factor GCN4 does not induce DNA bending in vitro. Substitution of basic residues for three neutral amino acids of GCN4 confers the ability to bend DNA. This result is consistent with a model of induced DNA bending wherein excess positive charge in proximity to one face of the double helix neutralizes local phosphate diester anions resulting in a laterally asymmetric charge distribution along the DNA. Previous data suggest that such an unbalanced charge distribution results in collapse of the DNA toward the neutralized surface. Interpretations of the present data are discussed. Our result supports the hypothesis that electrostatic interactions can play a key role in DNA bending by bZIP proteins.

Investigation of the pyrimidine preference by the c-Myb DNA-binding domain at the initial base of the consensus sequence

The principal determinant of the pyrimidine preference by the c-Myb DNA-binding domain at the initial base of the consensus sequence was investigated by mutation of both the protein and the DNA base pairs, with analysis by a filter binding assay. Amino acid residue 187 was revealed to interact with the pyrimidine base position, as estimated from our previous complex structure. Unexpectedly, since the pyrimidine preference is retained even in the Gly187 mutant, the principal origin of the base specificity should not occur via the direct-readout mechanism, but by an indirect-readout mechanism, namely in the intrinsic "bendability" of the pyrimidine-purine step of the DNA duplex. A significant but rather small positive base pair roll is detectable in the conformation of DNA in complex with the c-Myb DNA-binding domain. Following the conventional chemical rules of the direct-readout mechanism, amino acid mutagenesis at position 187 yielded several new base preferences for the protein.

Physical basis of a protein-DNA recognition code

Can a stereochemical recognition code explain sequence-specific protein-nucleic acid interactions? Whereas a code that is generally applicable to DNA-binding proteins of all known structural families is unattainable, the indications are that a code can describe at least some of the interactions of classical zinc fingers with DNA. The crystal structures of related zinc finger-DNA complexes reveal a remarkable mode of interaction that sets the framework for this code, and recent biochemical studies have elucidated the intermolecular contacts (contingent on this framework) that result in specificity.

Common DNA structural features exhibited by eukaryotic ribosomal gene promoters

Nucleotide sequences of DNA regions containing eukaryotic ribosomal promoters were analysed using strategies designed to reveal sequence-directed structural features. DNA curvature, duplex stability and pattern of twist angle variation were studied by computer modelling. Although ribosomal promoters are known to lack sequence homology (unless very closely related species are considered), investigation of these structural characteristics uncovered striking homologies in all the taxonomic groups examined so far. This wide conservation of DNA structures, while DNA sequence is not conserved, suggests that the determined structures are fundamental for ribosomal promoter function. Moreover, this result agrees well with the recent observations showing that RNA polymerase I transcription factors have not evolved as intensively as previously suspected.

Oct-1 POU domain-DNA interactions: cooperative binding of isolated subdomains and effects of covalent linkage

Structural and biochemical studies of Oct-1 POU domain-DNA interactions have raised important questions about cooperativity and the role of the linker connecting the POU-specific domain and the POU homeo domain. To analyze these interactions, we have studied binding of the isolated domains. Surprisingly, we find that two unlinked polypeptides corresponding to the POU-specific domain and the POU homeo domain bind cooperatively to the octamer site and have a coupling energy of 1.6 kcal/mole. We suggest that overlapping DNA contacts near the center of the octamer site may be the source of this cooperativity, as there are no protein-protein contacts between the domains in the crystal structure of the Oct-1 POU domain-DNA complex. These studies also have allowed us to describe the thermodynamic contribution of the linker (present in the intact POU domain) in terms of an effective concentration (3.6 mM). The broader implications for understanding cooperativity in protein-DNA recognition and gene regulation are discussed.

Crystal structure of LacI member, PurR, bound to DNA: minor groove binding by alpha helices

The three-dimensional structure of a ternary complex of the purine repressor, PurR, bound to both its corepressor, hypoxanthine, and the 16-base pair purF operator site has been solved at 2.7 A resolution by x-ray crystallography. The bipartite structure of PurR consists of an amino-terminal DNA-binding domain and a larger carboxyl-terminal corepressor binding and dimerization domain that is similar to that of the bacterial periplasmic binding proteins. The DNA-binding domain contains a helix-turn-helix motif that makes base-specific contacts in the major groove of the DNA. Base contacts are also made by residues of symmetry-related alpha helices, the "hinge" helices, which bind deeply in the minor groove. Critical to hinge helix-minor groove binding is the intercalation of the side chains of Leu54 and its symmetry-related mate, Leu54', into the central CpG-base pair step. These residues thereby act as "leucine levers" to pry open the minor groove and kink the purF operator by 45 degrees.

DNA structure, hydration and dynamics

Although the double helical model of DNA structure is now 40 years old, there is still considerable effort being made to elucidate the range of conformations that can be adopted by this flexible molecule. We review the current state of our knowledge of DNA structure which is available from both experimental and computational approaches.

Determination of common structural features in Escherichia coli promoters by computer analysis

Escherichia coli promoters show a large degree of sequence variation. However, they are all recognized specifically by RNA polymerase as the sites for transcription initiation, suggesting that they share common basic structural features distinguishing them from the rest of the sequence. Our hypothesis is that the promoter is determined not only by the two consensus sequences at -10 and -35, but also by the surrounding nucleotides, and that it is not only the identity of the nucleotides that is important for promoter function but the presence of specific physical-chemical and structural characteristics that are sequence dependent. This approach is supported by accumulating evidence indicating the role that the DNA conformation may play in modulating protein-DNA interaction. In this study, four intrinsic sequence-dependent characteristics are examined in E. coli promoter regions: helix stability, helix flexibility, and two conformational parameters represented by the DNA tendencies for B-->Z and B-->A transition. The promoter is defined by the consensus sequences and their vicinity and the examined properties are compared between promoter and random sequences. It is demonstrated that both the consensus and flanking regions are less stable, more flexible and show a higher tendency for the B conformation in comparison to random sequences. Discriminant analysis is used to evaluate the relative contributions of the various characteristics.

Modelling DNA conformational mechanics

Using a modelling technique specifically adapted to DNA helices, we have systematically studied the impact of base sequence on the geometry of the double helix. The results obtained show that each repetitive base sequence leads to several stable conformations belonging to the B-DNA family. These conformational sub-states generally have similar stabilities, but often differ considerably in terms of their helical and backbone parameters. Each sub-state can be characterised by the puckering of its sugar rings. Surface energy mapping and combinatorial search techniques are used to further understand the DNA conformational hypersurface and to extend our study from dinucleotide repeats to tetranucleotide sequences. The resulting structural database should be useful for predicting the properties of longer and more irregular base sequences and thus should contribute to understanding how DNA target sites are recognised.

Spatial translational motions of base pairs in DNA molecules: application of the extended matrix generator method

We have used the elementary generator matrices outlined in the preceding paper to examine the conformational plasticity of the nucleic acid double helix. Here we investigate kinked DNA structures made up of alternating B- and A-type helices and intrinsically curved duplexes perturbed by the intercalation of ligands. We model the B-to-A transition by the lateral translation of adjacent base pairs, and the intercalation of ligands by the vertical displacement of neighboring residues. We report a complete set of average configuration-dependent parameters, ranging from scalars (i.e., persistence lengths) to first- and second-order tensor parameters (i.e., average second moments of inertia), as well as approximations of the associated spatial distributions of the DNA and their angular correlations. The average structures of short chains (of lengths less than 100 base pairs) with local kinks or intrinsically curved sequences are essentially rigid rods. At the smallest chain lengths (10 base pairs), the kinked and curved chains exhibit similar average properties, although they are structurally perturbed compared to the standard B-DNA duplex. In contrast, at lengths of 200 base pairs, the curved and kinked chains are more compact on average and are located in a different space from the standard B- or A-DNA helix. While A-DNA is shorter and thicker than B-DNA in x-ray models, the long flexible A-DNA helix is thinner and more extended on average than its B-DNA counterpart because of more limited fluctuations in local structure. Curved polymers of 50 base pairs or longer also show significantly greater asymmetry than other DNAs (in terms of the distribution of base pairs with respect to the center of gravity of the chain). The intercalation of drugs in the curved DNA straightens and extends the smoothly deformed template. The dimensions of the average ellipsoidal boundaries defining the configurations of the intercalated polymers are roughly double those of the intrinsically curved chain. The altered proportions and orientations of these density functions reflect the changing shape and flexibility of the double helix. The calculations shed new light on the possible structural role of short A-DNA fragments in long B-type duplexes and also offer a model for understanding how GC-specific intercalative ligands can straighten naturally curved DNA. The mechanism is not immediately obvious from current models of DNA curvature, which attribute the bending of the chain to a perturbed structure in repeating tracts of A.T base pairs.

Crystal structure of a five-finger GLI-DNA complex: new perspectives on zinc fingers

Zinc finger proteins, of the type first discovered in transcription factor IIIA (TFIIIA), are one of the largest and most important families of DNA-binding proteins. The crystal structure of a complex containing the five Zn fingers from the human GLI oncogene and a high-affinity DNA binding site has been determined at 2.6 A resolution. Finger one does not contact the DNA. Fingers two through five bind in the major groove and wrap around the DNA, but lack the simple, strictly periodic arrangement observed in the Zif268 complex. Fingers four and five of GLI make extensive base contacts in a conserved nine base-pair region, and this section of the DNA has a conformation intermediate between B-DNA and A-DNA. Analyzing the GLI complex and comparing it with Zif268 offers new perspectives on Zn finger-DNA recognition.

Protein interactions with a gender-specific gene of Schistosoma mansoni: characterization by DNase I footprinting, band shift and UV cross-linking

A schistosome gender specific gene (F-10), was used as a probe to characterize DNA binding proteins from adult male and female Schistoma mansoni. Using the band-shift and DNase I footprinting methods, it was found that proteins from male and female worms bound to the intact F-10 gene and to restriction fragments corresponding to different domains of the gene, generating relatively long protected sites. Clear differences between male and female proteins were only observed when nuclear proteins were tested. Thus, gender-specific binding was detected in fragments corresponding to the 5' and 3' ends. UV-induced cross-linking between schistosome proteins and a synthetic oligonucleotide bearing a steroid response element present in the 3' untranslated end of the F-10 gene, revealed a major DNA binding protein with a molecular mass of 30 kDa, in both male and female worms. These results suggested that the activation of transcription of the F-10 gene may depend essentially on nuclear proteins.

Protein-induced bending as a transcriptional switch

The question of whether protein-induced DNA bending can act as a switch factor when placed upstream of an array of promoters located in tandem was investigated in vivo. The catabolite activating protein binding site of the fur operon was replaced by the binding site of the RepA repressor protein, which is able to bend DNA immediately after binding. Appropriately phased induced bending could act as a transcriptional switch factor in vivo.