Exploring the role of glutamine 50 in the homeodomain-DNA interface: crystal structure of engrailed (Gln50 --> ala) complex at 2.0 A

In Silico Characterization of African Swine Fever Virus Nucleoprotein p10 Interaction with DNA

African swine fever virus (ASFV) is the etiological agent of a highly contagious, hemorrhagic infectious swine disease, with a tremendous sanitary and economic impact on a global scale. Currently, there are no globally available vaccines or treatments. The p10 protein, a structural nucleoprotein encoded by ASFV, has been previously described as capable of binding double-stranded DNA (dsDNA), which may have implications for viral replication. However, the molecular mechanism that governs this interaction is still unknown, mostly due to the lack of a structural model for this protein. In this work, we have generated an ab initio model of the p10 protein and performed extensive structural characterization, using molecular dynamics simulations to identify the motifs and residues regulating DNA recognition. The helix-turn-helix motif identified at the C-terminal region of the protein was shown to be crucial to the dsDNA-binding efficiency. As with other DNA-binding proteins, two distinct serine and lysine-rich regions found in the two helices were identified as key players in the binding to DNA, whose importance was later validated using experimental binding assays. Altogether, these findings may contribute to a better understanding of the p10 function in ASFV replication.

Lysine Side-Chain Dynamics in the Binding Site of Homeodomain/DNA Complexes As Observed by NMR Relaxation Experiments and Molecular Dynamics Simulations

An important but poorly characterized contribution to the thermodynamics of protein-DNA interactions is the loss of entropy that occurs from restricting the conformational freedom of amino acid side chains. The effect of restricting the flexibility of several side chains at a protein-DNA interface may be comparable in many cases to the other factors that determine the binding thermodynamics and may, therefore, play a key role in dictating the binding affinity and/or specificity. Because the entropic contributions, including the presence and influence of side-chain dynamics, are especially difficult to estimate based on structural information, it is important to pursue experimental and theoretical studies that can provide direct information regarding these issues. We report on studies of a model system, the homeodomain/DNA complex, focusing on the Lys50 class of homeodomains where a key lysine residue in position 50 was shown previously to be critical for binding site specificity. NMR methodology was employed for determining the dynamics of lysine side-chain amino groups via 15N relaxation measurements in the Lys50-class homeodomains from the Drosophila protein Bicoid and the human protein Pitx2. In the case of Pitx2, complexes with both a consensus and a nonconsensus DNA binding site were examined. NMR-derived order parameters indicated moderate to substantial conformational freedom for the lysine NH3+ group in the complexes studied. To complement the experimental NMR measurements, molecular dynamics simulations were performed for the consensus complexes to gain further, detailed insights regarding the dynamics of the Lys50 side chain and other important residues in the protein-DNA interface.

Dynamic Local Polymorphisms in the Gbx1 Homeodomain Induced by DNA Binding

The Gastrulation Brain Homeobox 1 (Gbx1) gene encodes the Gbx1 homeodomain that targets TAATTA motifs in double-stranded DNA (dsDNA). Residues Glu17 and Arg52 in Gbx1 form a salt bridge, which is preserved in crystal structures and molecular dynamics simulations of homologous homeodomain-DNA complexes. In contrast, our nuclear magnetic resonance (NMR) studies show that DNA binding to Gbx1 induces dynamic local polymorphisms, which include breaking of the Glu17-Arg52 salt bridge. To study this interaction, we produced a variant with Glu17Arg and Arg52Glu mutations, which exhibited the same fold as the wild-type protein, but a 2-fold reduction in affinity for dsDNA. Analysis of the NMR structures of the Gbx1 homeodomain in the free form, the Gbx1[E17R,R52E] variant, and a Gbx1 homeodomain-DNA complex showed that stabilizing interactions of the Arg52 side chain with the DNA backbone are facilitated by transient breakage of the Glu17-Arg52 salt bridge in the DNA-bound Gbx1.

Dlx5 Homeodomain:DNA Complex: Structure, Binding and Effect of Mutations Related to Split Hand and Foot Malformation Syndrome

The Dlx5 homeodomain is a transcription factor related to the Drosophila distal-less gene that is associated with breast and lung cancer, lymphoma, Rett syndrome and osteoporosis in humans. Mutations in the DLX5 gene have been linked to deficiencies in craniofacial and limb development in higher eukaryotes, including split hand and foot malformation 1 in humans. Our characterization of a Dlx5 homeodomain:(CGACTAATTAGTCG)2 complex by NMR spectroscopy paved the way for determination of its crystal structure at 1.85Å resolution that enabled rationalization of the effects of disease-related mutations on the protein function. A Q186H mutation linked to split hand and foot malformation 1 likely affects affinity of DNA binding by disrupting water-mediated interactions with the DNA major groove. A more subtle effect is implicated for the Q178P mutation, which is not in direct contact with the DNA. Our data indicate that these mutations diminish the ability of the Dlx5 homeodomain to recognize and bind target DNAs, and they likely destabilize the formation of functional complexes.

Structure solution of DNA-binding proteins and complexes with ARCIMBOLDO libraries

Protein-DNA interactions play a major role in all aspects of genetic activity within an organism, such as transcription, packaging, rearrangement, replication and repair. The molecular detail of protein-DNA interactions can be best visualized through crystallography, and structures emphasizing insight into the principles of binding and base-sequence recognition are essential to understanding the subtleties of the underlying mechanisms. An increasing number of high-quality DNA-binding protein structure determinations have been witnessed despite the fact that the crystallographic particularities of nucleic acids tend to pose specific challenges to methods primarily developed for proteins. Crystallographic structure solution of protein-DNA complexes therefore remains a challenging area that is in need of optimized experimental and computational methods. The potential of the structure-solution program ARCIMBOLDO for the solution of protein-DNA complexes has therefore been assessed. The method is based on the combination of locating small, very accurate fragments using the program Phaser and density modification with the program SHELXE. Whereas for typical proteins main-chain α-helices provide the ideal, almost ubiquitous, small fragments to start searches, in the case of DNA complexes the binding motifs and DNA double helix constitute suitable search fragments. The aim of this work is to provide an effective library of search fragments as well as to determine the optimal ARCIMBOLDO strategy for the solution of this class of structures.

Human HOX gene disorders

The Hox genes are an evolutionarily conserved family of genes, which encode a class of important transcription factors that function in numerous developmental processes. Following their initial discovery, a substantial amount of information has been gained regarding the roles Hox genes play in various physiologic and pathologic processes. These processes range from a central role in anterior-posterior patterning of the developing embryo to roles in oncogenesis that are yet to be fully elucidated. In vertebrates there are a total of 39 Hox genes divided into 4 separate clusters. Of these, mutations in 10 Hox genes have been found to cause human disorders with significant variation in their inheritance patterns, penetrance, expressivity and mechanism of pathogenesis. This review aims to describe the various phenotypes caused by germline mutation in these 10 Hox genes that cause a human phenotype, with specific emphasis paid to the genotypic and phenotypic differences between allelic disorders. As clinical whole exome and genome sequencing is increasingly utilized in the future, we predict that additional Hox gene mutations will likely be identified to cause distinct human phenotypes. As the known human phenotypes closely resemble gene-specific murine models, we also review the homozygous loss-of-function mouse phenotypes for the 29 Hox genes without a known human disease. This review will aid clinicians in identifying and caring for patients affected with a known Hox gene disorder and help recognize the potential for novel mutations in patients with phenotypes informed by mouse knockout studies.

Binding polymorphism in the DNA bound state of the Pdx1 homeodomain

The subtle effects of DNA-protein recognition are illustrated in the homeodomain fold. This is one of several small DNA binding motifs that, in spite of limited DNA binding specificity, adopts crucial, specific roles when incorporated in a transcription factor. The homeodomain is composed of a 3-helix domain and a mobile N-terminal arm. Helix 3 (the recognition helix) interacts with the DNA bases through the major groove, while the N-terminal arm becomes ordered upon binding a specific sequence through the minor groove. Although many structural studies have characterized the DNA binding properties of homeodomains, the factors behind the binding specificity are still difficult to elucidate. A crystal structure of the Pdx1 homeodomain bound to DNA (PDB 2H1K) obtained previously in our lab shows two complexes with differences in the conformation of the N-terminal arm, major groove contacts, and backbone contacts, raising new questions about the DNA recognition process by homeodomains. Here, we carry out fully atomistic Molecular Dynamics simulations both in crystal and aqueous environments in order to elucidate the nature of the difference in binding contacts. The crystal simulations reproduce the X-ray experimental structures well. In the absence of crystal packing constraints, the differences between the two complexes increase during the solution simulations. Thus, the conformational differences are not an artifact of crystal packing. In solution, the homeodomain with a disordered N-terminal arm repositions to a partially specific orientation. Both the crystal and aqueous simulations support the existence of different stable binding conformers identified in the original crystallographic data with different degrees of specificity. We propose that protein-protein and protein-DNA interactions favor a subset of the possible conformations. This flexibility in DNA binding may facilitate multiple functions for the same transcription factor.

Discovery, optimization and validation of an optimal DNA-binding sequence for the Six1 homeodomain transcription factor

The Six1 transcription factor is a homeodomain protein involved in controlling gene expression during embryonic development. Six1 establishes gene expression profiles that enable skeletal myogenesis and nephrogenesis, among others. While several homeodomain factors have been extensively characterized with regards to their DNA-binding properties, relatively little is known of the properties of Six1. We have used the genomic binding profile of Six1 during the myogenic differentiation of myoblasts to obtain a better understanding of its preferences for recognizing certain DNA sequences. DNA sequence analyses on our genomic binding dataset, combined with biochemical characterization using binding assays, reveal that Six1 has a much broader DNA-binding sequence spectrum than had been previously determined. Moreover, using a position weight matrix optimization algorithm, we generated a highly sensitive and specific matrix that can be used to predict novel Six1-binding sites with highest accuracy. Furthermore, our results support the idea of a mode of DNA recognition by this factor where Six1 itself is sufficient for sequence discrimination, and where Six1 domains outside of its homeodomain contribute to binding site selection. Together, our results provide new light on the properties of this important transcription factor, and will enable more accurate modeling of Six1 function in bioinformatic studies.

Exploring the DNA-recognition potential of homeodomains

The recognition potential of most families of DNA-binding domains (DBDs) remains relatively unexplored. Homeodomains (HDs), like many other families of DBDs, display limited diversity in their preferred recognition sequences. To explore the recognition potential of HDs, we utilized a bacterial selection system to isolate HD variants, from a randomized library, that are compatible with each of the 64 possible 3' triplet sites (i.e., TAANNN). The majority of these selections yielded sets of HDs with overrepresented residues at specific recognition positions, implying the selection of specific binders. The DNA-binding specificity of 151 representative HD variants was subsequently characterized, identifying HDs that preferentially recognize 44 of these target sites. Many of these variants contain novel combinations of specificity determinants that are uncommon or absent in extant HDs. These novel determinants, when grafted into different HD backbones, produce a corresponding alteration in specificity. This information was used to create more explicit HD recognition models, which can inform the prediction of transcriptional regulatory networks for extant HDs or the engineering of HDs with novel DNA-recognition potential. The diversity of recovered HD recognition sequences raises important questions about the fitness barrier that restricts the evolution of alternate recognition modalities in natural systems.

Origins of specificity in protein-DNA recognition

Specific interactions between proteins and DNA are fundamental to many biological processes. In this review, we provide a revised view of protein-DNA interactions that emphasizes the importance of the three-dimensional structures of both macromolecules. We divide protein-DNA interactions into two categories: those when the protein recognizes the unique chemical signatures of the DNA bases (base readout) and those when the protein recognizes a sequence-dependent DNA shape (shape readout). We further divide base readout into those interactions that occur in the major groove from those that occur in the minor groove. Analogously, the readout of the DNA shape is subdivided into global shape recognition (for example, when the DNA helix exhibits an overall bend) and local shape recognition (for example, when a base pair step is kinked or a region of the minor groove is narrow). Based on the >1500 structures of protein-DNA complexes now available in the Protein Data Bank, we argue that individual DNA-binding proteins combine multiple readout mechanisms to achieve DNA-binding specificity. Specificity that distinguishes between families frequently involves base readout in the major groove, whereas shape readout is often exploited for higher resolution specificity, to distinguish between members within the same DNA-binding protein family.

Novel structural features in two ZHX homeodomains derived from a systematic study of single and multiple domains

BACKGROUND

Zhx1 to 3 (zinc-fingers and homeoboxes) form a set of paralogous genes encoding multi-domain proteins. ZHX proteins consist of two zinc fingers followed by five homeodomains. ZHXs have biological roles in cell cycle control by acting as co-repressors of the transcriptional regulator Nuclear Factor Y. As part of a structural genomics project we have expressed single and multi-domain fragments of the different human ZHX genes for use in structure determination.

RESULTS

A total of 30 single and multiple domain ZHX1-3 constructs selected from bioinformatics protocols were screened for soluble expression in E. coli using high throughput methodologies. Two homeodomains were crystallized leading to structures for ZHX1 HD4 and ZHX2 HD2. ZHX1 HD4, although closest matched to homeodomains from 'homez' and 'engrailed', showed structural differences, notably an additional C-terminal helix (helix V) which wrapped over helix I thereby making extensive contacts. Although ZHX2 HD2-3 was successfully expressed and purified, proteolysis occurred during crystallization yielding crystals of just HD2. The structure of ZHX2 HD2 showed an unusual open conformation with helix I undergoing 'domain-swapping' to form a homodimer.

CONCLUSIONS

Although multiple-domain constructs of ZHX1 selected by bioinformatics studies could be expressed solubly, only single homeodomains yielded crystals. The crystal structure of ZHX1 HD4 showed additional hydrophobic interactions relative to many known homeodomains via extensive contacts formed by the novel C-terminal helix V with, in particular, helix I. Additionally, the replacement of some charged covariant residues (which are commonly observed to form salt bridges in non-homeotherms such as the Drosophila 'engrailed' homeodomain), by apolar residues further increases hydrophobic contacts within ZHX1 HD4, and potentially stability, relative to engrailed homeodomain. ZHX1 HD4 helix V points away from the normally observed DNA major groove binding site on homeodomains and thus would not obstruct the putative binding of nucleic acid. In contrast, for ZHX2 HD2 the observed altered conformation involving rearrangement of helix I, relative to the canonical homeodomain fold, disrupts the normal DNA binding site, although protein-protein binding is possible as observed in homodimer formation.

Comparison of multiple crystal structures with NMR data for engrailed homeodomain

Two methods are currently available to solve high resolution protein structures-X-ray crystallography and nuclear magnetic resonance (NMR). Both methods usually produce highly similar structures, but small differences between both solutions are always observed. Here the raw NMR data as well as the solved NMR structure were compared to the multiple crystal structures solved for the WT 60 residue three helix bundle engrailed homeodomain (EnHD) and single point mutants. There was excellent agreement between TALOS-predicted and crystal structure-observed dihedral angles and a good agreement for the (3) J(H ( N ) H ( alpha )) couplings for the multiple crystal structures. Around 1% of NOEs were violated for any crystal structure, but no NOE was inconsistent with all of the crystal structures. Violations usually occurred for surface residues or for residues for which multiple discreet conformations were observed between the crystal structures. Comparison of the disorder shown in the multiple crystal structures shows little correlation with dynamics under native conditions for this protein.

A one-dimensional reaction coordinate for identification of transition states from explicit solvent P(fold)-like calculations

A properly identified transition state ensemble (TSE) in a molecular dynamics (MD) simulation can reveal a tremendous amount about how a protein folds and offer a point of comparison to experimentally derived Phi(F) values, which reflect the degree of structure in these transient states. In one such method of TSE identification, dubbed P(fold), MD simulations of individual protein structures taken from an unfolding trajectory are used to directly assess an input structure's probability of folding before unfolding, and P(fold) is, by definition, 0.5 for the TSE. Other, less computationally intensive methods, such as multidimensional scaling (MDS) of the pairwise root mean-squared deviation (RMSD) matrix of the conformations sampled in a thermal unfolding trajectory, have also been used to identify the TSE. Identification of the TSE is made from the original MD simulation without the need to run further simulations. Here we present a P(fold)-like study and describe methods for identification of the TSE through the derivation of a high fidelity, bounded, one-dimensional reaction coordinate for protein folding. These methods are applied to the engrailed homeodomain. The TSE identified by this approach is essentially identical to the TSE identified previously by MDS of the pairwise RMSD matrix. However, the cost of performing P(fold), or even our reduced P(fold)-like calculations, is at least 36,000 times greater than the MDS method.

Energetics of the protein-DNA-water interaction

Background

To understand the energetics of the interaction between protein and DNA we analyzed 39 crystallographically characterized complexes with the HINT (Hydropathic INTeractions) computational model. HINT is an empirical free energy force field based on solvent partitioning of small molecules between water and 1-octanol. Our previous studies on protein-ligand complexes demonstrated that free energy predictions were significantly improved by taking into account the energetic contribution of water molecules that form at least one hydrogen bond with each interacting species.

Results

An initial correlation between the calculated HINT scores and the experimentally determined binding free energies in the protein-DNA system exhibited a relatively poor r² of 0.21 and standard error of ± 1.71 kcal mol^-1. However, the inclusion of 261 waters that bridge protein and DNA improved the HINT score-free energy correlation to an r² of 0.56 and standard error of ± 1.28 kcal mol^-1. Analysis of the water role and energy contributions indicate that 46% of the bridging waters act as linkers between amino acids and nucleotide bases at the protein-DNA interface, while the remaining 54% are largely involved in screening unfavorable electrostatic contacts.

Conclusion

This study quantifies the key energetic role of bridging waters in protein-DNA associations. In addition, the relevant role of hydrophobic interactions and entropy in driving protein-DNA association is indicated by analyses of interaction character showing that, together, the favorable polar and unfavorable polar/hydrophobic-polar interactions (i.e., desolvation) mostly cancel.

Structure and properties of a re-engineered homeodomain protein-DNA interface

The homeodomain (HD)-DNA interface has been conserved over 500 million years of evolution. Despite this conservation, we have successfully re-engineered the engrailed HD to specifically recognize an unnatural nucleotide using a phage display selection. Here we report the synthesis of novel nucleosides and the selection of mutant HDs that bind these nucleotides using phage display. The high-resolution crystal structure of one mutant in complex with modified and unmodified DNA demonstrates that, even with the substantial perturbation to the interface, this selected mutant retains a canonical HD structure. Dissection of the contributions due to each of the selected mutations reveals that the majority of the modification-specific binding is accomplished by a single mutation (I47G) but that the remaining mutations retune the stability of the HD. These results afford a detailed look at a re-engineered protein-DNA interaction and provide insight into the opportunities for re-engineering highly conserved interfaces.

The solution structure of the native K50 Bicoid homeodomain bound to the consensus TAATCC DNA-binding site

The solution structure of the homeodomain of the Drosophila morphogenic protein Bicoid (Bcd) complexed with a TAATCC DNA site is described. Bicoid is the only known protein that uses a homeodomain to regulate translation, as well as transcription, by binding to both RNA and DNA during early Drosophila development; in addition, the Bcd homeodomain can recognize an array of different DNA sites. The dual functionality and broad recognition capabilities signify that the Bcd homeodomain may possess unique structural/dynamic properties. Bicoid is the founding member of the K50 class of homeodomain proteins, containing a lysine residue at the critical 50th position (K50) of the homeodomain sequence, a residue required for DNA and RNA recognition; Bcd also has an arginine residue at the 54th position (R54), which is essential for RNA recognition. Bcd is the only known homeodomain with the K50/R54 combination of residues. The Bcd structure indicates that this homeodomain conforms to the conserved topology of the homeodomain motif, but exhibits a significant variation from other homeodomain structures at the end of helix 1. A key result is the observation that the side-chains of the DNA-contacting residues K50, N51 and R54 all show strong signs of flexibility in the protein-DNA interface. This finding is supportive of the adaptive-recognition theory of protein-DNA interactions.

Structural basis of Prospero-DNA interaction: implications for transcription regulation in developing cells

The crystal structure of a complex between the novel homeodomain of the neural transcription factor Prospero and DNA shows that the invariant residues Lys1290, Asn1294, and Asp1297 make specific contacts with the noncanonical DNA binding site. The overall structure includes the homeodomain and the adjacent Prospero domain and confirms that they act as a single structural unit, a Homeo-Prospero domain. The Prospero domain facilitates the proper alignment of the protein on the DNA. Knowledge of the structure reconciles two different DNA sequences that have been proposed as transcriptional targets for Prospero. As in the apo structure, the C terminus of the Prospero domain shields a short helix within the homeodomain that includes a nuclear export signal (NES). The structural results suggest that exposure of the NES is not coupled directly to DNA binding. We propose a DNA recognition mechanism specific to Prospero-type homeodomains in developing cells.

Homeodomain revisited: a lesson from disease-causing mutations

The homeodomain is a highly conserved DNA-binding motif that is found in numerous transcription factors throughout a large variety of species from yeast to humans. These gene-specific transcription factors play critical roles in development and adult homeostasis, and therefore, any germline mutations associated with these proteins can lead to a number of congenital abnormalities. Although much has been revealed concerning the molecular architecture and the mechanism of homeodomain-DNA interactions, the study of disease-causing mutations can further provide us with instructive information as to the role of particular residues in a conserved mode of action. In this paper, I have compiled the homeodomain missense mutations found in various human diseases and re-examined the functional role of the mutational "hot spot" residues in light of the structures obtained from crystallography. These findings should be useful in understanding the essential components of the homeodomain and in attempts to design agonist or antagonists to modulate their activity and to reverse the effects caused by the mutations.

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.

Specific DNA recognition by theAntp homeodomain: MD simulations of specific and nonspecific complexes

Four molecular dynamics simulation trajectories of complexes between the wild-type or a mutant Antennapedia homeodomain and 2 DNA sequences were generated in order to probe the mechanisms governing the specificity of DNA recognition. The starting point was published affinity measurements showing that a single protein mutation combined with a replacement of 2 base pairs yields a new high-affinity complex, whereas the other combinations, with changes on only 1 macromolecule, exhibited lower affinity. The simulations of the 4 complexes yielded fluctuating networks of interaction. On average, these networks differ significantly, explaining the switch of affinity caused by the alterations in the macromolecules. The network of mostly hydrogen-bonding interactions involving several water molecules, which was suggested both by X-ray and NMR structures of the wild-type homeodomain and its DNA operator sequence, could be reproduced in the trajectory. More interestingly, the high-affinity complex with alterations in both the protein and the DNA yielded again a dynamic but very tight network of intermolecular interactions, however, attributing a significantly stronger role to direct hydrophobic interactions at the expense of water bridges. The other 2 homeodomain-DNA complexes, with only 1 molecule altered, show on average over the trajectories a clearly reduced number of protein-DNA interactions. The observations from these simulations suggest specific experiments and thus close the circle formed by biochemical, structural, and computational studies. The shift from a water-dominated to a more "dry" interface may prove important in the design of proteins binding DNA in a specific manner.

Crystal structures of engrailed homeodomain mutants: implications for stability and dynamics

We report the crystal structures and biophysical characterization of two stabilized mutants of the Drosophila Engrailed homeodomain that have been engineered to minimize electrostatic repulsion. Four independent copies of each mutant occupy the crystal lattice, and comparison of these structures illustrates variation that can be partly ascribed to networks of correlated conformational adjustments. Central to one network is leucine 26 (Leu26), which occupies alternatively two side chain rotameric conformations (-gauche and trans) and different positions within the hydrophobic core. Similar sets of conformational substates are observed in other Engrailed structures and in another homeodomain. The pattern of structural adjustments can account for NMR relaxation data and sequence co-variation networks in the wider homeodomain family. It may also explain the dysfunction associated with a P26L mutation in the human ARX homeodomain protein. Finally, we observe a novel dipolar interaction between a conserved tryptophan and a water molecule positioned along the normal to the indole ring. This interaction may explain the distinctive fluorescent properties of the homeodomain family.

Molecular characterization of NF-HEV, a nuclear factor preferentially expressed in human high endothelial venules

Lymphocyte homing to secondary lymphoid tissue and lesions of chronic inflammation is directed by multi-step interactions between the circulating cells and the specialized endothelium of high endothelial venules (HEVs). In this study, we used the PCR-based method of suppression subtractive hybridization (SSH) to identify novel HEV genes by comparing freshly purified HEV endothelial cells (HEVECs) with nasal polyp-derived microvascular endothelial cells (PMECs). By this approach, we cloned the first nuclear factor preferentially expressed in HEVECs, designated nuclear factor from HEVs (NF-HEV). Virtual Northern and Western blot analyses showed strong expression of NF-HEV in HEVECs, compared to human umbilical vein endothelial cells (HUVECs) and PMECs. In situ hybridization and immunohistochemistry revealed that NF-HEV mRNA and protein are expressed at high levels and rather selectively by HEVECs in human tonsils, Peyers's patches, and lymph nodes. The NF-HEV protein was found to contain a bipartite nuclear localization signal, and was targeted to the nucleus when ectopically expressed in HUVECs and HeLa cells. Furthermore, endogenous NF-HEV was found in situ to be confined to the nucleus of tonsillar HEVECs. Finally, threading and molecular modeling studies suggested that the amino-terminal part of NF-HEV (aa 1-60) corresponds to a novel homeodomain-like Helix-Turn-Helix (HTH) DNA-binding domain. Similarly to the atypical homeodomain transcription factor Prox-1, which plays a critical role in the induction of the lymphatic endothelium phenotype, NF-HEV may be one of the key nuclear factors that controls the specialized HEV phenotype.

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.

Structural and thermodynamic characterization of the DNA binding properties of a triple alanine mutant of MATalpha2

Triply mutated MATalpha2 protein, alpha2-3A, in which all three major groove-contacting residues are mutated to alanine, is defective in binding DNA alone or in complex with Mcm1 yet binds with MATa1 with near wild-type affinity and specificity. To gain insight into this unexpected behavior, we determined the crystal structure of the a1/alpha2-3A/DNA complex. The structure shows that the triple mutation causes a collapse of the alpha2-3A/DNA interface that results in a reorganized set of alpha2-3A/DNA contacts, thereby enabling the mutant protein to recognize the wild-type DNA sequence. Isothermal titration calorimetry measurements reveal that a much more favorable entropic component stabilizes the a1/alpha2-3A/DNA complex than the alpha2-3A/DNA complex. The combined structural and thermodynamic studies provide an explanation of how partner proteins influence the sequence specificity of a DNA binding protein.

Testing water-mediated DNA recognition by the Hin recombinase

The Hin recombinase specifically recognizes its DNA-binding site by means of both major and minor groove interactions. A previous X-ray structure, together with new structures of the Hin DNA-binding domain bound to a recombination half-site that were solved as part of the present study, have revealed that two ordered water molecules are present within the major groove interface. In this report, we test the importance of these waters directly by X-ray crystal structure analysis of complexes with four mutant DNA sequences. These structures, combined with their Hin-binding properties, provide strong support for the critical importance of one of the intermediate waters. A lesser but demonstrable role is ascribed to the second water molecule. The mutant structures also illustrate the prominent roles of thymine methyls both in stabilizing intermediate waters and in interfering with water or amino acid side chain interactions with DNA.

Altering the DNA-binding specificity of the yeast Matalpha 2 homeodomain protein

Homeodomain proteins are a highly conserved class of DNA-binding proteins that are found in virtually every eukaryotic organism. The conserved mechanism that these proteins use to bind DNA suggests that there may be at least a partial DNA recognition code for this class of proteins. To test this idea, we have investigated the sequence-specific requirements for DNA binding and repression by the yeast alpha2 homeodomain protein in association with its cofactors, Mcm1 and Mata1. We have determined the contribution for each residue in the alpha2 homeodomain that contacts the DNA in the co-crystal structures of the protein. We have also engineered mutants in the alpha2 homeodomain to alter the DNA-binding specificity of the protein. Although we were unable to change the specificity of alpha2 by making substitutions at residues 47, 54, and 55, we were able to alter the DNA-binding specificity by making substitutions at residue 50 in the homeodomain. Since other homeodomain proteins show similar changes in specificity with substitutions at residue 50, this suggests that there is at least a partial DNA recognition code at this position.

Spectrum of mutations and genotype-phenotype analysis in Currarino syndrome

The triad of a presacral tumour, sacral agenesis and anorectal malformation constitutes the Currarino syndrome which is caused by dorsal-ventral patterning defects during embryonic development. The syndrome occurs in the majority of patients as an autosomal dominant trait associated with mutations in the homeobox gene HLXB9 which encodes the nuclear protein HB9. However, genotype-phenotype analyses have been performed only in a few families and there are no reports about the specific impact of HLXB9 mutations on HB9 function. We performed a mutational analysis in 72 individuals from nine families with Currarino syndrome. We identified a total of five HLXB9 mutations, four novel and one known mutation, in four out of four families and one out of five sporadic cases. Highly variable phenotypes and a low penetrance with half of all carriers being clinically asymptomatic were found in three families, whereas affected members of one family showed almost identical phenotypes. However, an obvious genotype-phenotype correlation was not found. While HLXB9 mutations were diagnosed in 23 patients, no mutation or microdeletion was detected in four sporadic patients with Currarino syndrome. The distribution pattern of here and previously reported HLXB9 mutations indicates mutational predilection sites within exon 1 and the homeobox. Furthermore, sequence homology to Drosophila homeobox genes suggest that some of these mutations located within the homeobox may alter the DNA-binding specificity of HB9 while those in sequences homologous to a recently identified NLS motif of the human homeobox gene PDX-1 may impair nuclear translocation of the mutated protein.

Regulation of an endogenous locus using a panel of designed zinc finger proteins targeted to accessible chromatin regions. Activation of vascular endothelial growth factor A

We have mapped conserved regions of enhanced DNase I accessibility within the endogenous chromosomal locus of vascular endothelial growth factor A (VEGF-A). Synthetic zinc finger protein (ZFP) transcription factors were designed to target DNA sequences contained within the DNase I-hypersensitive regions. These ZFPs, when fused to either VP16 or p65 transcriptional activation domains, were able to activate expression of the VEGF-A gene as assayed by mRNA accumulation and VEGF-A protein secretion through a range exceeding that induced by hypoxic stress. Importantly, multiple splice variants of VEGF-A mRNA with defined physiological functions were induced by a single engineered ZFP transcription factor. We present evidence for an enhanced activation of VEGF-A gene transcription by ZFP transcription factors fused to VP16 and p65 targeted to two distinct chromosomal sites >500 base pairs upstream or downstream of the transcription start site. Our strategy provides a novel approach for dissecting the requirements for gene regulation at a distance without altering the DNA sequence of the endogenous target locus.