The vnd/NK-2 homeodomain: thermodynamics of reversible unfolding and DNA binding for wild-type and with residue replacements H52R and H52R/T56W in helix III

The conformational stabilities of the vnd (ventral nervous system defective)/NK-2 homeodomain [HD(wt); residues 1-80 that encompass the 60-residue homeodomain] and those harboring mutations in helix III of the DNA recognition site [HD(H52R) and HD(H52R/T56W)] have been investigated by differential scanning calorimetry (DSC) and ellipticity changes at 222 nm. Thermal unfolding reactions at pH 7.4 are reversible and repeatable in the presence of 50-500 mM NaCl with DeltaC(p) = 0.52 +/- 0.04 kcal K(-1) mol(-1). A substantial stabilization of HD(wt) is produced by 50 mM phosphate or by the addition of 100-500 mM NaCl to 50 mM Hepes, pH 7.4, buffer (from T(m) = 35.5 degrees C to T(m) 43-51 degrees C; DeltaH(vH) congruent with 47 +/- 5 kcal mol(-1)). The order of stability is HD(H52R/T56W) > HD(H52R) > HD(wt), irrespective of the anions present. Progress curves for ellipticity changes at 222 nm as a function of increasing temperature are fitted well by a two-state unfolding model, and the cooperativity of secondary structure changes is greater for mutant homeodomains than for HD(wt) and also is increased by adding 100 mM NaCl to Hepes buffer. A 33% quench of the intrinsic tryptophanyl residue fluorescence of HD(wt) by phosphate binding (K(D)' = 2.6 +/- 0.3 mM phosphate) is reversed approximately 60% by DNA binding. Thermodynamic parameters for vnd/NK-2 homeodomain proteins binding sequence-specific 18 bp DNA have been determined by isothermal titration calorimetry (10-30 degrees C). Values of DeltaC(p) are +0.25, -0.17, and -0.10 +/- 0.04 kcal K(-1) mol(-1) for HD(wt), HD(H52R), and HD(H52R/T56W) binding duplex DNA, respectively. Interactions of homeodomains with DNA are enthalpically controlled at 298 K and pH 7.4 with corresponding DeltaH values of -6.6 +/- 0.5, -10.8 +/- 0.1, and -9.0 +/- 0.6 kcal mol(-1) and DeltaG' values of -11.0 +/- 0.1, -11.0 +/- 0.1, and -11.3 +/- 0.3 kcal mol(-1) with a binding stoichiometry of 1.0 +/- 0.1. Thermodynamic parameters for DNA binding are not predicted from homeodomain structural changes that occur upon complexing to DNA and must reflect also solvent and possibly DNA rearrangements.     

The three-dimensional structure of the vnd/NK-2 homeodomain-DNA complex by NMR spectroscopy.

The three-dimensional solution structure obtained by NMR of the complex formed between the uniformly singly15N and doubly13C/15N-labeled vnd/NK-2 homeodomain and its consensus 16 base-pair DNA binding sequence was determined. This work was carried out using the accepted repertoire of experiments augmented with a novel implementation of the water flipback technique to enhance signals from exchangeable amide protons. The results using this new technique confirm the existence of hydrogen bonding between the invariant Asn51 and the second adenine of the DNA binding sequence, as seen in crystal structures of other homeodomain-DNA complexes, but never before detected by NMR. Hydrogen bonding by Arg5 and Lys3 in the minor groove of the DNA appears to be responsible for two unusually upfield-shifted ribose H1' resonances. The DNA duplex is nearly straight and its structure is primarily that of B -DNA. A detailed comparison is presented for all available homeodomain-DNA structures including the vnd/NK-2 DNA complex, which demonstrates that homology is maintained in the protein structure, whereas for the orientation of the homeodomain relative to DNA, small but significant variations are observed. Interactions are described involving certain residues in specific positions of the homeodomain, namely Leu7, Thr41, and Gln50 of vnd/NK-2, where single amino acid residue mutations lead to dramatic developmental alterations. The availability of our previously determined three- dimensional structure of the vnd/NK-2 homeodomain in the absence of DNA allows us to assess structural changes in the homeodomain induced by DNA binding.     

Cardiac-specific Nkx2.5 homeodomain: conformational stability and specific DNA binding of Nkx2.5(C56S)

The cardiac-specific Nkx2.5 homeodomain has been expressed as a 79-residue protein with the oxidizable Cys(56) replaced with Ser. The Nkx2.5 or Nkx2.5(C56S) homeodomain is 73% identical in sequence to and has the same NMR structure as the vnd (ventral nervous system defective)/NK-2 homeodomain of Drosophila when bound to the same specific DNA. The thermal unfolding of Nkx2.5(C56S) at pH 6.0 or 7.4 is a reversible, two-state process with unit cooperativity, as measured by differential scanning calorimetry (DSC) and far-UV circular dichroism. Adding 100 mM NaCl to Nkx2.5(C56S) at pH 7.4 increases T(m) from 44 to 54 +/- 0.2 degrees C and DeltaH from 34 to 45 +/- 2 kcal/mol (giving a DeltaC(p) of approximately 1.2 kcal K(-)(1) mol(-)(1) for homeodomain unfolding). DSC profiles of Nkx2.5 indicate fluctuating nativelike structures at <37 degrees C. Titrations of specific 18 bp DNA with Nkx2.5(C56S) in buffer at pH 7.4 with 100 mM NaCl yield binding constants of 2-6 x 10(8) M(-)(1) from 10 to 37 degrees C and a stoichiometry of 1:1 for homeodomain binding DNA, using isothermal titration calorimetry. The DNA binding reaction of Nkx2.5 is enthalpically controlled, and the temperature dependence of DeltaH gives a DeltaC(p) of -0.18 +/- 0.01 kcal K(-)(1) mol(-)(1). This corresponds to 648 +/- 36 A(2) of buried apolar surface upon Nkx2.5(C56S) binding duplex B-DNA. Thermodynamic parameters differ for Nkx2.5 and vnd/NK-2 homeodomains binding specific DNA. Unbound NK-2 is more flexible than Nkx2.5.     

Quantitative measurement of water diffusion lifetimes at a protein/DNA interface by NMR

Hydration site lifetimes of slowly diffusing water molecules at the protein/DNA interface of the vnd/NK-2 homeodomain DNA complex were determined using novel three-dimensional NMR techniques. The lifetimes were calculated using the ratios of ROE and NOE cross-relaxation rates between the water and the protein backbone and side chain amides. This calculation of the lifetimes is based on a model of the spectral density function of the water-protein interaction consisting of three timescales of motion: fast vibrational/rotational motion, diffusion into/out of the hydration site, and overall macromolecular tumbling. The lifetimes measured ranged from approximately 400 ps to more than 5 ns, and nearly all the slowly diffusing water molecules detected lie at the protein/DNA interface. A quantitative analysis of relayed water cross-relaxation indicated that even at very short mixing times, 5 ms for ROESY and 12 ms for NOESY, relay of magnetization can make a small but detectable contribution to the measured rates. The temperature dependences of the NOE rates were measured to help discriminate direct dipolar cross-relaxation from chemical exchange. Comparison with several X-ray structures of homeodomain/DNA complexes reveals a strong correspondence between water molecules in conserved locations and the slowly diffusing water molecules detected by NMR. A homology model based on the X-ray structures was created to visualize the conserved water molecules detected at the vnd/NK-2 homeodomain DNA interface. Two chains of water molecules are seen at the right and left sides of the major groove, adjacent to the third helix of the homeodomain. Two water-mediated hydrogen bond bridges spanning the protein/DNA interface are present in the model, one between the backbone of Phe8 and a DNA phosphate, and one between the side chain of Asn51 and a DNA phosphate. The hydrogen bond bridge between Asn51 and the DNA might be especially important since the DNA contact made by the invariant Asn51 residue, seen in all known homeodomain/DNA structures, is critical for binding affinity and specificity.     

Interaction of the PHD-finger homeodomain protein HAT3.1 from Arabidopsis thaliana with DNA. Specific DNA binding by a homeodomain with histidine at position 51

HAT3.1 is a member of the PHD-finger homeodomain protein family. The HAT3.1 homeodomain is highly divergent in sequence even at positions that are almost invariable among homeodomains. In this work, we have applied the random oligonucleotide selection technique to investigate if the HAT3.1 homeodomain is able to recognize specific DNA sequences. Analysis of the selected molecules followed by hydroxyl radical footprinting experiments and yeast one-hybrid assays indicated that HAT3.1 shows a preference for the sequence T(A/G)(A/C)ACCA, different from those bound by other homeodomains. Binding was dependent on homeodomain residues located at positions 47, 50, 51, and 54, the same positions that usually participate in DNA binding in most homeodomains. The study of the interaction of mutants at these positions with DNA carrying nucleotide changes at specific sites suggested that H51 and K50 most likely interact with nucleotides 2 to 4 and 5 to 6, respectively, while W54 would establish contacts with position 4. The presence of H51 and W54 represents an innovation among homeodomain structures. The fact that the HAT3.1 homeodomain is able to interact with specific DNA sequences is evidence of the inherent plasticity of the homeodomain as a DNA binding unit.     

The interaction with DNA of wild-type and mutant fushi tarazu homeodomains.

The in vitro DNA binding properties of wild-type and mutant fushi tarazu homeodomains (ftz HD) have been analysed. The DNA binding properties of the ftz HD are very similar to those of the Antp HD. In interference experiments with mutant ftz HDs, close approaches between specific portions of the ftz HD peptide and specific regions of the binding site DNA were mapped. A methylation interference, G7 on the beta strand of BS2, is absent from the interference pattern with a mutant ftz HD [ftz (R43A) HD] in which the Arg43 at the second position of helix III (the recognition helix) is replaced by an Ala. This indicated that Arg43 of the ftz HD is in close proximity to the N7 of G7 of the beta strand of BS2 in the major groove. The methylation and ethylation interference patterns with the ftz (NTD) HD, in which the first six amino acids of the homeodomain were deleted, were extensively altered relative to the ftz HD patterns. Methylation of A11 and G12 of the alpha strand and ethylation of the phosphate of nucleotide A12 of the alpha strand no longer interfere with binding. This indicated that the first six amino acids of the homeodomain of ftz interact with A11 of the alpha strand in the minor groove, the phosphate of the nucleotide A13 on the alpha strand and G12 of the alpha strand in the adjacent major groove of BS2. In a binding study using a change of specificity mutation [ftz (Q50K) HD], in which the Gln50 at the ninth position of the third helix is exchanged for a Lys (as in the bicoid HD), and variant binding sites, we concluded that position 50 of the ftz HD and the ftz (Q50K) HD peptides interacts with base pairs at positions 6 and 7 of BS2. These three points of contact allowed us to propose a crude orientation of the ftz HD within the protein-DNA complex. We find that the ftz HD and the Antp HD peptides contact DNA in a similar way.     

Crystal structure of a MAT alpha 2 homeodomain-operator complex suggests a general model for homeodomain-DNA interactions.

The MAT alpha 2 homeodomain regulates the expression of cell type-specific genes in yeast. We have determined the 2.7 A resolution crystal structure of the alpha 2 homeodomain bound to a biologically relevant DNA sequence. The DNA in this complex is contacted primarily by the third of three alpha-helices, with additional contacts coming from an N-terminal arm. Comparison of the yeast alpha 2 and the Drosophila engrailed homeodomain-DNA complexes shows that the protein fold is highly conserved, despite a 3-residue insertion in alpha 2 and only 27% sequence identity between the two homeodomains. Moreover, the orientation of the recognition helix on the DNA is also conserved. This docking arrangement is maintained by side chain contacts with the DNA--primarily the sugar-phosphate backbone--that are identical in alpha 2 and engrailed. Since these residues are conserved among all homeodomains, we propose that the contacts with the DNA are also conserved and suggest a general model for homeodomain-DNA interactions.     

Insights into nonspecific binding of homeodomains from a structure of MATalpha2 bound to DNA

The 2.1-A resolution crystal structure of the MATalpha2 homeodomain bound to DNA reveals the unexpected presence of two nonspecifically bound alpha2 homeodomains, in addition to the two alpha2 homeodomains bound to canonical alpha2 binding sites. One of the extra homeodomains makes few base-specific contacts, while the other extra homeodomain binds to DNA in a previously unobserved manner. This unusually bound homeodomain is rotated on the DNA, making possible major groove contacts by side-chains that normally do not contact the DNA. This alternate docking may represent one way in which homeodomains sample nonspecific DNA sequences.     

A single amino acid can determine the DNA binding specificity of homeodomain proteins

Many Drosophila developmental genes contain a DNA binding domain encoded by the homeobox. This homeodomain contains a region distantly homologous to the helix-turn-helix motif present in several prokaryotic DNA binding proteins. We investigated the nature of homeodomain-DNA interactions by making a series of mutations in the helix-turn-helix motif of the Drosophila homeodomain protein Paired (Prd). This protein does not recognize sequences bound by the homeodomain proteins Fushi tarazu (Ftz) or Bicoid (Bcd). We show that changing a single amino acid at the C-terminus of the recognition helix is both necessary and sufficient to confer the DNA binding specificity of either Ftz or Bcd on Prd. This simple rule indicates that the amino acids that determine the specificity of homeodomains are different from those mediating protein-DNA contacts in prokaryotic proteins. We further show that Prd contains two DNA binding activities. The Prd homeodomain is responsible for one of them while the other is not dependent on the recognition helix.     

Role of salt bridges in homeodomains investigated by structural analyses and molecular dynamics simulations

Homeodomains are a class of helix-turn-helix DNA-binding protein motifs that play an important role in the control of cellular development in eukaryotes. They fold in a three alpha-helix structural module, where the third helix is the recognition helix that fits into the major groove of DNA. Structural analysis of the members of the homeodomain family led to the identification of interactions likely to stabilize the protein domains. Linking the helices pairwise, three salt bridges were found to be well preserved within the family. Also well conserved were two cation-pi interactions between aromatic and positively charged side chains. To analyze the structural role of the salt bridges, molecular dynamics simulations (MD) were carried out on the wild-type homeodomain from the Drosophila paired protein (1fjl) and on three mutants, which lack one or two salt bridges and mimic natural mutations in other homeodomains. Analysis of the trajectories revealed only small structural rearrangements of the three helices in all MD simulations, thereby suggesting that the salt bridges have no essential stabilizing role at room temperature, but rather might be important for improving thermostability. The latter hypothesis is supported by a good correlation between the melting midpoint temperatures of several homeodomains and the number of salt bridges and cation-pi interactions that connect secondary structures.