Characterization of the bacterially expressed Drosophila engrailed homeodomain.

To investigate the mechanism of DNA recognition by the homeodomain, truncated proteins containing the entire homeodomain encoded by the Drosophila engrailed gene were expressed in Escherichia coli. Each protein was accumulated to an amount representing more than 40% of the total bacterial protein and recovered in the soluble fraction. Of the three truncated proteins, the shortest one (71 amino acid residues) was further purified by conventional chromatography. The purified engrailed homeodomain (En-HD) protected a DNA sequence, TTAATT, the core element of consensus sequences recognized by many other homeodomain proteins, from DNase I digestion. UV-CD spectra of the En-HD showed that it mainly consisted of alpha-helix. Based on one-dimensional 1H-NMR spectra, the tertiary structure of the En-HD was shown to be stable against temperature up to 50 degrees C and low pH. The low pH resistance of the protein was also demonstrated by UV-CD measurement. Thus, the current over-production system provides an active and stable homeodomain, which is suitable for structure-function analysis.     

Members of TALE and WUS subfamilies of homeodomain proteins with potentially important functions in development form dimers within each subfamily in rice

Transacting factors often form homo- and heterodimers and regulate various targets, the type of regulation depending on the dimeric combination. The WUS and TALE subfamilies are two atypical homeodomains in plants. A homeodomain mediates sequence-specific binding to its target DNA and usually consists of 60 amino acid residues, whereas atypical homeodomains have extra amino acid residues in the well-conserved region. The genes OsWUS and OsPRS, which encode atypical homeodomain proteins from the WUS subfamily, and OsBEL and OSH15, which encode those from the TALE subfamily, were isolated from rice and tested for their interactions by yeast two-hybrid analysis. OsWUS and OsPRS formed homodimers and formed heterodimers with each other but did not form dimers with the TALE family homeodomain proteins OSH15 or OsBEL. Likewise, OSH15 and OsBEL formed homodimers and heterodimers but did not form dimers with the WUS family homeodomain proteins OsWUS and OsPRS. These findings suggest that the combinations of dimers are well correlated with the classification of these proteins on the basis of sequence similarity. RT-PCR analysis revealed that expression of OsWUS and OsPRS was detected in the same organs, namely floral buds, roots, and suspension cells. Therefore, it is possible that the proteins encoded by both of these genes function as homo- and heterodimers in planta. These results suggest that, during the evolution of these subfamilies, various combinations of dimers within proteins encoded by paralogous genes were formed and generated independent regulatory networks that enabled complex patterns of plant development.     

Structure of the adenovirus E4 Orf6 protein predicted by fold recognition and comparative protein modeling

To facilitate investigation of the molecular and biochemical functions of the adenovirus E4 Orf6 protein, we sought to derive three-dimensional structural information using computational methods, particularly threading and comparative protein modeling. The amino acid sequence of the protein was used for secondary structure and hidden Markov model (HMM) analyses, and for fold recognition by the ProCeryon program. Six alternative models were generated from the top-scoring folds identified by threading. These models were examined by 3D-1D analysis and evaluated in the light of available experimental evidence. The final model of the E4 protein derived from these and additional threading calculations was a chimera, with the tertiary structure of its C-terminal 226 residues derived from a TIM barrel template and a mainly alpha-nonbundle topology for its poorly conserved N-terminal 68 residues. To assess the accuracy of this model, additional threading calculations were performed with E4 Orf6 sequences altered as in previous experimental studies. The proposed structural model is consistent with the reported secondary structure of a functionally important C-terminal sequence and can account for the properties of proteins carrying alterations in functionally important sequences or of those that disrupt an unusual zinc-coordination motif.     

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.     

A test of proposed rules for helix capping: implications for protein design

alpha-helices within proteins are often terminated (capped) by distinctive configurations of the polypeptide chain. Two common arrangements are the Schellman motif and the alternative alpha(L) motif. Rose and coworkers developed stereochemical rules to identify the locations of such motifs in proteins of unknown structure based only on their amino acid sequences. To check the effectiveness of these rules, they made specific predictions regarding the structural and thermodynamic consequences of certain mutations in T4 lysozyme. We have constructed these mutants and show here that they have neither the structure nor the stability that was predicted. The results show the complexity of the protein-folding problem. Comparison of known protein structures may show that a characteristic sequence of amino acids (a sequence motif) corresponds to a conserved structural motif. In any particular protein, however, changes in other parts of the sequence may result in a different conformation. The structure is determined by sequence as a whole, not by parts considered in isolation.     

Cloning and sequencing of a leech homolog to the Drosophila engrailed gene.

We have cloned and sequenced a homolog (ht-en) to the Drosophila engrailed (en) gene from the glossiphoniid leech, Helobdella triserialis. Amino acid comparisons of the ht-en homeodomain and C-terminal residues with the corresponding residues encoded by en-class genes of other species reveal 75-79% sequence identity. In addition, the ht-en sequence appears to have a serine-rich region 16 residues C-terminal from the homeodomain, which by analogy to Drosophila may be a target site for phosphorylation. The leech gene encodes some amino acid substitutions for residues that are highly conserved in other species. These are found within the second and third of the three putative helices of the homeodomain, and in both of the intervening turn regions.     

Conformational changes in the PBX homeodomain and C-terminal extension upon binding DNA and HOX-derived YPWM peptides

PBX is a member of the three amino acid loop extension (TALE) class of homeodomains. PBX binds DNA cooperatively with HOX homeodomain proteins that contain a conserved YPWM motif. The amino acids immediately C-terminal to the PBX homeodomain increase the affinity of the homeodomain for its DNA site and HOX proteins. We have determined the structure of the free PBX homeodomain using NMR spectroscopy. Both the PBX homeodomain and the extended PBX homeodomain make identical contacts with a 5'-TGAT-3' DNA site and a YPWM peptide. A fourth alpha-helix, which forms upon binding to DNA, stabilizes the extended PBX structure. Variations in DNA sequence selectivity of heterodimeric PBX-HOX complexes depend on the HOX partner; however, a comparison of five different HOX-derived YPWM peptides showed that each bound to PBX in the same way, differing only in the strength of the association.     

Use of a structural alphabet to find compatible folds for amino acid sequences

The structural annotation of proteins with no detectable homologs of known 3D structure identified using sequence‐search methods is a major challenge today. We propose an original method that computes the conditional probabilities for the amino‐acid sequence of a protein to fit to known protein 3D structures using a structural alphabet, known as “Protein Blocks” (PBs). PBs constitute a library of 16 local structural prototypes that approximate every part of protein backbone structures. It is used to encode 3D protein structures into 1D PB sequences and to capture sequence to structure relationships. Our method relies on amino acid occurrence matrices, one for each PB, to score global and local threading of query amino acid sequences to protein folds encoded into PB sequences. It does not use any information from residue contacts or sequence‐search methods or explicit incorporation of hydrophobic effect. The performance of the method was assessed with independent test datasets derived from SCOP 1.75A. With a Z‐score cutoff that achieved 95% specificity (i.e., less than 5% false positives), global and local threading showed sensitivity of 64.1% and 34.2%, respectively. We further tested its performance on 57 difficult CASP10 targets that had no known homologs in PDB: 38 compatible templates were identified by our approach and 66% of these hits yielded correctly predicted structures. This method scales‐up well and offers promising perspectives for structural annotations at genomic level. It has been implemented in the form of a web‐server that is freely available at http://www.bo‐protscience.fr/forsa .     

The DxDxDG motif for calcium binding: multiple structural contexts and implications for evolution

Calcium ions regulate many cellular processes and have important structural roles in living organisms. Despite the great variety of calcium-binding proteins (CaBPs), many of them contain the same Ca(2+)-binding helix-loop-helix structure, referred to as the EF-hand. In the canonical EF-hand, the loop contains three calcium-binding aspartic acid residues, which form the DxDxDG sequence motif, and is flanked by two alpha-helices. Recently, other CaBPs containing the same motif, but lacking one or both helices, have been described. Here, structural motif searches were used to analyse the full diversity of structural context in the known set of DxDxDG-containing CaBPs, including those where the structural resemblance of a given DxDxDG motif to that of EF-hands had not been noted. The results obtained indicate that the EF-hand represents but one, among many, structural context for the DxDxDG-like Ca(2+)-binding loops. While the structural similarity of the binuclear calcium-binding sites in anthrax protective antigen and human thrombospondin suggests that they are homologous, evolutionary relationships for mononuclear sites are harder to discern. The possible scenarios for the evolution of DxDxDG motif-containing calcium-binding loops in a variety of non-homologous proteins suggested loop transplant as a mechanism perhaps responsible for much of the diversity in structural contexts of present day DxDxDG-type CaBPs. Additionally, while it can be shown that existence of a DxDxDG sequence is not enough to confer a conformation suitable for calcium binding, local convergent evolution may still have a role. The analysis presented here has consequences for the prediction of calcium binding from sequence alone.     

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.     

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.     

Structure-based prediction of DNA target sites by regulatory proteins

Regulatory proteins play a critical role in controlling complex spatial and temporal patterns of gene expression in higher organism, by recognizing multiple DNA sequences and regulating multiple target genes. Increasing amounts of structural data on the protein-DNA complex provides clues for the mechanism of target recognition by regulatory proteins. The analyses of the propensities of base-amino acid interactions observed in those structural data show that there is no one-to-one correspondence in the interaction, but clear preferences exist. On the other hand, the analysis of spatial distribution of amino acids around bases shows that even those amino acids with strong base preference such as Arg with G are distributed in a wide space around bases. Thus, amino acids with many different geometries can form a similar type of interaction with bases. The redundancy and structural flexibility in the interaction suggest that there are no simple rules in the sequence recognition, and its prediction is not straightforward. However, the spatial distributions of amino acids around bases indicate a possibility that the structural data can be used to derive empirical interaction potentials between amino acids and bases. Such information extracted from structural databases has been successfully used to predict amino acid sequences that fold into particular protein structures. We surmised that the structures of protein-DNA complexes could be used to predict DNA target sites for regulatory proteins, because determining DNA sequences that bind to a particular protein structure should be similar to finding amino acid sequences that fold into a particular structure. Here we demonstrate that the structural data can be used to predict DNA target sequences for regulatory proteins. Pairwise potentials that determine the interaction between bases and amino acids were empirically derived from the structural data. These potentials were then used to examine the compatibility between DNA sequences and the protein-DNA complex structure in a combinatorial "threading" procedure. We applied this strategy to the structures of protein-DNA complexes to predict DNA binding sites recognized by regulatory proteins. To test the applicability of this method in target-site prediction, we examined the effects of cognate and noncognate binding, cooperative binding, and DNA deformation on the binding specificity, and predicted binding sites in real promoters and compared with experimental data. These results show that target binding sites for several regulatory proteins are successfully predicted, and our data suggest that this method can serve as a powerful tool for predicting multiple target sites and target genes for regulatory proteins.     

A phage display selection of engrailed homeodomain mutants and the importance of residue Q50

Mutants of engrailed homeodomain (HD) that retain DNA-binding activity were isolated using a phage display selection. This selection was used to enrich for active DNA-binding clones from a complex library consisting of over a billion members. A more focused library of mutant homeodomains consisting of all possible amino acid combinations at two DNA-contacting residues (I47 and Q50) was constructed and screened for members capable of binding tightly and specifically to the engrailed consensus sequence, TAATTA. The isolated mutants largely recapitulated the distribution of amino acids found at these positions in natural homeodomains thus validating the in vitro selection conditions. In particular, the unequivocal advantage enjoyed by glutamine at residue 50 is surprising in light of reports that minimize the importance of this residue. Here, the subtle contributions of residue Q50 are demonstrated to play a functionally important role in specific recognition of DNA. These results highlight the complex subtlety of protein–DNA interactions, underscoring the value of the first reported in vitro selection of a homeodomain.