Diffusing and colliding: the atomic level folding/unfolding pathway of a small helical protein

Proteins with ultra-fast folding/unfolding kinetics are excellent candidates for study by molecular dynamics. Here, we describe such simulations of a three helix bundle protein, the engrailed homeodomain (En-HD), which folds via the diffusion-collision model. The unfolding pathway of En-HD was characterized by seven simulations of the protein and 12 simulations of its helical fragments yielding over 1.1 micros of simulation time in water. Various conformational states along the unfolding pathway were identified. There is the compact native-like transition state, a U-shaped helical intermediate and an unfolded state with dynamic helical segments. Each of these states is in good agreement with experimental data. Examining these states as well as the transitions between them, we find the role of long-range tertiary contacts, specifically salt-bridges, important in the folding/unfolding pathway. In the folding direction, charged residues form long-range tertiary contacts before the hydrophobic core is formed. The formation of HII is assisted by a specific salt-bridge and by non-specific (fluctuating) tertiary contacts, which we call contact-assisted helix formation. Salt-bridges persist as the protein approaches the transition state, stabilizing HII until the hydrophobic core is formed. To complement this information, simulations of fragments of En-HD illustrate the helical propensities of the individual segments. By thermal denaturation, HII proved to be the least stable helix, unfolding in less than 450 ps at high temperature. We observed the low helical propensity of C-terminal residues from HIII in fragment simulations which, when compared to En-HD unfolding simulations, link the unraveling of HIII to the initial event that drives the unfolding of En-HD.     

Drosophila engrailed-1,10-phenanthroline chimeras as probes of homeodomain-DNA complexes.

We have converted the Drosophila engrailed homeodomain into a sequence-specific nuclease by linking the protein to the chemical nuclease 1,10-phenanthroline-copper (OP-Cu). Unique cysteines were introduced at six positions into the homeodomain by site-directed mutagenesis for the covalent attachment of OP-Cu. The varied DNA-binding affinity and specificity of these mutants and the DNA cleavage pattern of their OP-Cu derivatives allowed us to assess the crystal structure of the engrailed homeodomain-DNA complex. We have also achieved site-specific double-stranded DNA scission with one of the homeodomain mutants, E28C, which has the potential of being used to identify engrailed binding sites in the genome. Because the homeodomain is so well conserved among members of the homeodomain-containing protein family, other homeodomain proteins can be converted into nucleases by attaching OP-Cu at position 28 of their homeodomains.     

Molecular dynamics analysis of the engrailed homeodomain-DNA recognition

Molecular dynamics (MD) simulations were performed for investigating the role of Gln50 in the engrailed homeodomain-DNA recognition. Employing the crystal structure of free engrailed homeodomain and homeodomain-DNA complex as a starting structure, we carried out MD simulations of: (i) the complex between engrailed homeodomain and a 20 base-pair DNA containing TAATTA core sequence; (ii) the free engrailed homeodomain. The simulations show that homeodomain flexibility does not depend on its ligation state. The engrailed homeodomain shows similar flexibility, and the recognition helix-3 shows very similar characteristic of high rigidity and limited conformational space in two complexation states. At the same time, DNA structure has also no obvious conformational fluctuations. These results preclude the possibility of the side chain of Gln50 forming direct hydrogen bonds to the core DNA bases. MD simulations confirm a few well-conserved sites for water-mediated hydrogen bonds from protein to DNA are occupied by water molecules, and Gln50 interacts with corresponding core DNA bases through water-mediated hydrogen bonds. So Gln50 plays a relatively modest role in determining the affinity and specificity of the engrailed homeodomain. In addition, the electrostatic interaction between homeodomain and phosphate backbone of the DNA is a main factor for N- and C-terminal arm becoming ordered upon DNA binding.     

Mutational analysis of the engrailed homeodomain recognition helix by phage display.

The homeodomain (HD) is a ubiquitous protein fold that confers DNA binding function on a superfamily of eukaryotic gene regulatory proteins. Here, the DNA binding of recognition helix variants of the HD from the engrailed gene of Drosophila melanogaster was investigated by phage display. Nineteen different combinations of pairwise mutations at positions 50 and 54 were screened against a panel of four DNA sequences consisting of the engrailed consensus, a non-specific DNA control based on the lambda repressor operator OR1 and two model sequence targets con-taining imperfect versions of the 5'-TAAT-3' consensus. The resulting mutant proteins could be divided into four groups that varied with respect to their affinity for DNA and specificity for the engrailed consensus. The altered specificity phenotypes of several mutant proteins were confirmed by DNA mobility shift analysis. Lys50/Ala54 was the only mutant protein that exhibited preferential binding to a sequence other than the engrailed consensus. Arginine was also demonstrated to be a functional replacement for Ala54. The functional combinations at 50 and 54 identified by these experiments recapitulate the distribution of naturally occurring HD sequences and illustrate how the engrailed HD can be used as a framework to explore covariation among DNA binding residues.     

Crystallization and preliminary X-ray diffraction studies of the engrailed homeodomain and of an engrailed homeodomain/DNA complex

The homeodomain from the engrailed protein of Drosophila has been crystallized from ammonium phosphate at pH 6.8. The crystals form in space group P6(1)22 (or P6(5)22), with cell dimensions a = b = 44.8 A and c = 118.2 A. These crystals diffract to 1.8 A resolution. A complex containing the engrailed homeodomain and a duplex DNA site also has been crystallized. The cocrystals form in space group C2 with a = 131.2 A, b = 45.5 A, c = 72.9 A and beta = 119.0 degrees. These crystals diffract to 2.6 A resolution.     

The sequence specificity of homeodomain-DNA interaction

The Drosophila developmental gene, engrailed, encodes a sequence-specific DNA binding activity. Using deletion constructs expressed as fusion proteins in E. coli, we localized this activity to the conserved homeodomain (HD). The binding site consensus, TCAATTAAAT, is found in clusters in the engrailed regulatory region. Weak binding of the En HD to one copy of a synthetic consensus is enhanced by adjacent copies. The distantly related HD encoded by fushi tarazu binds to the same sites as the En HD, but differs in its preference for related sites. Both HDs bind a second type of sequence, a repeat of TAA. The similarity in sequence specificity of En and Ftz HDs suggests that, within families of DNA binding proteins, close relatives will exhibit similar specificities. Competition among related regulatory proteins might govern which protein occupies a given binding site and consequently determine the ultimate effect of cis-acting regulatory sites.     

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.     

Structural studies of the engrailed homeodomain.

The structure of the Drosophila engrailed homeodomain has been solved by molecular replacement and refined to an R-factor of 19.7% at a resolution of 2.1 A. This structure offers a high-resolution view of an important family of DNA-binding proteins and allows comparison to the structure of the same protein bound to DNA. The most significant difference between the current structure and that of the 2.8-A engrailed-DNA complex is the close packing of an extended strand against the rest of the protein in the unbound protein. Structural features of the protein not previously noted include a "herringbone" packing of 4 aromatic residues in the core of the protein and an extensive network of salt bridges that covers much of the helix 1-helix 2 surface. Other features that may play a role in stabilizing the native state include the interaction of buried carbonyl oxygen atoms with the edge of Phe 49 and a bias toward statistically preferred side-chain dihedral angles. There is substantial disorder at both ends of the 61 amino acid protein. A 51-amino acid variant of engrailed (residues 6-56) was synthesized and shown by CD and thermal denaturation studies to be structurally and thermodynamically similar to the full-length domain.     

Fourier transform infrared/vibrational circular dichroism spectroscopy as an informative tool for the investigation of large supramolecular complexes of biological macromolecules

A combination of ultraviolet (UV) and infrared (IR) absorption and circular dichroism (CD) spectroscopy was applied to investigate the structure and formation of large supramolecular DNA-protein complexes. This combination of techniques was used to overcome limitations of UV-CD (electronic, or ECD) spectroscopy due to considerable light scattering in such solutions. Based on the analysis of FTIR and UV-CD spectra, the interaction of DNA with nonhistone chromatin protein HMGB1 and linker histone H1 was studied. The data obtained showed that under the conditions of the experiment (15 mM NaCl, protein/DNA ratio r < 1 w/w) the proteins did not reveal any AT or GC specificity in binding to DNA. In the presence of both proteins, mainly interactions in the DNA minor groove were observed, which were attributed to HMGB1 binding. Histone H1 facilitated binding of HMGB1 to DNA by interacting with the negatively charged groups of the sugar-phosphate backbone and binding of aspartic and glutamic amino acid residues of HMGB1. Acting together, HMGB1 and H1 stimulated the assemblage of supramolecular DNA-protein structures. The structural organization of the ternary complexes depended not only on the properties of the protein-DNA interactions but also on the interactions between HMGB1 and H1 molecules.     

Full-sequence computational design and solution structure of a thermostable protein variant

Computational protein design procedures were applied to the redesign of the entire sequence of a 51 amino acid residue protein, Drosophila melanogaster engrailed homeodomain. Various sequence optimization algorithms were compared and two resulting designed sequences were experimentally evaluated. The two sequences differ by 11 mutations and share 22% and 24% sequence identity with the wild-type protein. Both computationally designed proteins were considerably more stable than the naturally occurring protein, with midpoints of thermal denaturation greater than 99 degrees C. The solution structure was determined for one of the two sequences using multidimensional heteronuclear NMR spectroscopy, and the structure was found to closely match the original design template scaffold.     

Expression of engrailed proteins in arthropods, annelids, and chordates

engrailed is a homeobox gene that has an important role in Drosophila segmentation. Genes homologous to engrailed have been identified in several other organisms. Here we describe a monoclonal antibody that recognizes a conserved epitope in the homeodomain of engrailed proteins of a number of different arthropods, annelids, and chordates; we use this antibody to isolate the grasshopper engrailed gene. In Drosophila embryos, the antibody reveals engrailed protein in the posterior portion of each segment during segmentation, and in a segmentally reiterated subset of neuronal cells during neurogenesis. Other arthropods, including grasshopper and two crustaceans, have similar patterns of engrailed expression. However, these patterns of expression are not shared by the annelids or chordates we examined. Our results provide the most comprehensive view that has been obtained of how expression patterns of a regulatory gene vary during evolution. On the basis of these patterns, we suggest that engrailed is a gene whose ancestral function was in neurogenesis and whose function was co-opted during the evolution of segmentation in the arthropods, but not in the annelids and chordates.     

Structure of a HoxB1-Pbx1 heterodimer bound to DNA: role of the hexapeptide and a fourth homeodomain helix in complex formation

Hox homeodomain proteins are developmental regulators that determine body plan in a variety of organisms. A majority of the vertebrate Hox proteins bind DNA as heterodimers with the Pbx1 homeodomain protein. We report here the 2.35 A structure of a ternary complex containing a human HoxB1-Pbx1 heterodimer bound to DNA. Heterodimer contacts are mediated by the hexapeptide of HoxB1, which binds in a pocket in the Pbx1 protein formed in part by a three-amino acid insertion in the Pbx1 homeodomain. The Pbx1 DNA-binding domain is larger than the canonical homeodomain, containing an additional alpha helix that appears to contribute to binding of the HoxB1 hexapeptide and to stable binding of Pbx1 to DNA. The structure suggests a model for modulation of Hox DNA binding activity by Pbx1 and related proteins.     

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.