Molecular Dynamics Simulation in Solvent of the Estrogen Receptor Protein DNA Binding Domain in Complex with a Non-Consensus Estrogen Response Element DNA Sequence

Abstract

We investigated protein/DNA interactions, using molecular dynamics simulations computed between a 10 Angstom water layer model of the estrogen receptor (ER) protein DNA binding domain (DBD) amino acids and DNA of a non-consensus estrogen response element (ERE) consisting of 29 nucleotide base pairs. This ERE nucleotide sequence occurs naturally upstream of the Xenopus laevis Vitelligenin AI gene. The ER DBD is encoded by three exons. Namely, exons 2 and 3 which encode the two zinc binding motifs and a sequence of exon 4 which encodes a predicted alpha helix. We generated a computer model of the ER DBD using atomic coordinates derived from the average of 30 nuclear magnetic resonance (NMR) spectroscopy coordinate sets. Amino acids on the carboxyl end of the ER DBD were disordered in both X-ray crystallography and NMR determinations and no coordinates were reported. This disordered region includes 10 amino acids of a predicted alpha helix encoded in exon 4 at the exon 3/4 splice junction. These amino acids are known to be important in DNA binding and are also believed to function as a nuclear translocation signal sequence for the ER protein. We generated a computer model of the predicted alpha helix consisting of the 10 amino acids encoded in exon 4 and attached this helix to the carboxyl end of the ER DBD at the exon 3/4 splice junction site. We docked the ER DBD model within the DNA major groove halfsites of the 29 base pair non-consensus ERE and flanking nucleotides. We constructed a solvated model with the ER DBD/ERE complex surrounded by a ten Angstrom water layer and conducted molecular dynamics simulations. Hydrogen bonding interactions were monitored. In addition, van der Waals and electrostatic interaction energies were calculated. Amino acids of the ER DBD DNA recognition helix formed both direct and water mediated hydrogen bonds at cognate codon-anticodon nucleotide base and backbone sites within the ERE DNA right major groove halfsite. Amino acids of the ER DBD exon 4 encoded predicted alpha helix formed direct and water mediated H-bonds with base and backbone sites of their cognate codon-anticodon nucleotides within the minor grooves flanking the ERE DNA major groove halfsites. These interactions together induced bending of the DNA into the protein.     

Conformational dynamics and molecular recognition: backbone dynamics of the estrogen receptor DNA-binding domain.

We examined the internal mobility of the estrogen receptor DNA-binding domain (ER DBD) using NMR15N relaxation measurements and compared it to that of the glucocorticoid receptor DNA-binding domain (GR DBD). The studied protein fragments consist of residues Arg183-His267 of the human ER and residues Lys438-Gln520 of the rat GR. The15N longitudinal (R1) and transverse (R2) relaxation rates and steady state {1H}-15N nuclear Overhauser enhancements (NOEs) were measured at 30 degrees C at1H NMR frequencies of 500 and 600 MHz. The NOE versus sequence profile and calculated order parameters for ER DBD backbone motions indicate enhanced internal dynamics on pico- to nanosecond time-scales in two regions of the core DBD. These are the extended strand which links the DNA recognition helix to the second zinc domain and the larger loop region of the second zinc domain. The mobility of the corresponding regions of the GR DBD, in particular that of the second zinc domain, is more limited. In addition, we find large differences between the ER and GR DBDs in the extent of conformational exchange mobility on micro- to millisecond time-scales. Based on measurements of R2as a function of the15N refocusing (CPMG) delay and quantitative (Lipari-Szabo-type) analysis, we conclude that conformational exchange occurs in the loop of the first zinc domain and throughout most of the second zinc domain of the ER DBD. The conformational exchange dynamics in GR DBD is less extensive and localized to two sites in the second zinc domain. The different dynamical features seen in the two proteins is consistent with previous studies of the free state structures in which the second zinc domain in the ER DBD was concluded to be disordered whereas the corresponding region of the GR DBD adopts a stable fold. Moreover, the regions of the ER DBD that undergo conformational dynamics on the micro- to millisecond time-scales in the free state are involved in intermolecular protein-DNA and protein-protein interactions in the dimeric bound state. Based on the present data and the previously published dynamical and DNA binding properties of a GR DBD triple mutant which recognize an ER binding site on DNA, we argue that the free state dynamical properties of the nuclear receptor DBDs is an important element in molecular recognition upon DNA binding.     

Conformational analysis by NMR and distance geometry techniques of a peptide mimetic of the third helix of the Antennapedia homeodomain.

The Antennapedia homeodomain structure consists of four helices. The helices II and III are connected by a tripeptide that forms a turn, and constitute the well-known helix-turn-helix motif. The recognition helix penetrates the DNA major groove, gives specific protein-DNA contacts and forms direct, or water-mediated, intermolecular hydrogen bonds. It was suggested that helix III (and perhaps also helix IV) might represent the recognition helix of Antennapedia homeodomain, which makes contact with the surface of the major groove of the DNA. In an attempt to clarify the helix III capabilities of assuming an helical conformation when separated from the rest of the protein, we carried out the structural determination of the recognition helix III in different solvent media. The conformational study of fragments 42-53, where residues W48 and F49, not involved in the protein-DNA interaction, were substituted by two alanines, was conducted in sodium dodecyl sulfate (SDS), trifluoroethanol (TFE) and TFE/water, using circular dichroism, nuclear magnetic resonance (NMR) and distance geometry (DG) techniques. The fragment assumes a well-defined secondary structure in TFE and in TFE/water (90/10, v/v) with an alpha-helix encompassing residues 4-9, while in TFE/water (70/30, v/v) a less regular structure was found. The DG results in the micellar system evidence the presence of a distorted alpha-helical conformation involving residues 4-8. Our results reveal that the isolated Antennapedia recognition helix III tend to preserve in solution the alpha-helical conformation even if separated from the rest of the molecule.     

DNA-binding mechanism of the monomeric orphan nuclear receptor NGFI-B.

The 2.7 A X-ray crystal structure of the DNA-binding domain (DBD) of the orphan nuclear receptor, nerve growth factor-induced-B (NGFI-B), complexed to its high-affinity DNA target, represents the first structure analysis of a nuclear receptor DBD bound as a monomer to DNA. The structure of the core DBD and its interactions with the major groove of the DNA are similar to previously crystallographically solved DBD-DNA complexes in this superfamily; however, residues C-terminal to this core form a separate and unique substructure that interacts extensively and in a sequence-specific way with the minor groove of its DNA target, in particular with the characteristic 3 A-T base-pair identity element that extends 5' to the usual nuclear receptor half-site (AGGTCA).     

Monomeric complex of human orphan estrogen related receptor-2 with DNA: a pseudo-dimer interface mediates extended half-site recognition

While most nuclear receptors bind DNA as homo or heterodimers, the human estrogen related receptors (hERRs) are members of a subfamily of orphan receptors that bind DNA as monomers. We have determined the solution structure of the DNA binding domain (DBD) of hERR2 bound to its cognate DNA. The structure and base interactions of the core DBD are similar to those of other nuclear receptors. However, high-affinity, sequence-specific DNA binding as a monomer necessitates formation of additional base contacts outside the core DBD. This is accomplished using a modified guanosine-binding "AT-hook" within the C-terminal extension (CTE) flanking the DBD, which makes base-specific minor groove interactions. The structure of the CTE is stabilized both by interactions with the DNA and by packing against a region of the core DBD normally reserved for dimerization. This pseudo-dimer interface provides a basis for the expansion of DNA recognition and suggests a mechanism through which dimerization may have evolved from an ancestral monomeric receptor.     

Identification of protein contact sites within the glucocorticoid/progestin response element

The glucocorticoid receptor (GR) and the progestin receptor (PR) bind specifically to a variety of DNA sequences, glucocorticoid/progestin response elements (GRE/PRE), located in the proximity of responsive gene promoters. Using the isolated recombinant GR DNA-binding domain (DBD), it has recently been shown that GR interacts with the GRE/PRE, a 15-basepair partially palindromic consensus sequence, as a dimer. In this study an investigation into the GR-GRE/PRE and PR-GRE/PRE interaction has been performed using missing base contact analysis with the tyrosine aminotransferase GREII (TATII) and recombinant GR DBD as well as a fusion protein consisting of the PR DBD fused to Staph. aureus protein-A. GR and PR had identical base contact points, localized within two consecutive major grooves, binding to the same face of the DNA. Ethylation interference was also performed on the GR DBD-TATII interaction. The contact points with the backbone phosphate groups flank the contacts within the major groove for each of the two half-sites. Knowledge of the contact points within the DNA sequence together with the three-dimensional structure of the protein enables modelling of the protein-DNA interaction.     

Interchange of DNA-binding modes in the deformed and ultrabithorax homeodomains: a structural role for the N-terminal arm

The deformed (Dfd) and ultrabithorax (Ubx) homeoproteins regulate developmental gene expression in Drosophila melanogaster by binding to specific DNA sequences within its genome. DNA binding is largely accomplished via a highly conserved helix-turn-helix DNA-binding domain that is known as a homeodomain (HD). Despite nearly identical DNA recognition helices and similar target DNA sequence preferences, the in vivo functions of the two proteins are quite different. We have previously revealed differences between the two HDs in their interactions with DNA. In an effort to define the individual roles of the HD N-terminal arm and recognition helix in sequence-specific binding, we have characterized the structural details of two Dfd/Ubx chimeric HDs in complex with both the Dfd and Ubx-optimal-binding site sequences. We utilized hydroxyl radical cleavage of DNA to assess the positioning of the proteins on the binding sites. The effects of missing nucleosides and purine methylation on HD binding were also analyzed. Our results show that both the Dfd and Ubx HDs have similar DNA-binding modes when in complex with the Ubx-optimal site. There are subtle but reproducible differences in these modes that are completely interchanged when the Dfd N-terminal arm is replaced with the corresponding region of the Ubx HD. In contrast, we showed previously that the Dfd-optimal site sequence elicits a very different binding mode for the Ubx HD, while the Dfd HD maintains a mode similar to that elicited by the Ubx-optimal site. Our current methylation interference studies suggest that this alternate binding mode involves interaction of the Ubx N-terminal arm with the minor groove on the opposite face of DNA relative to the major groove that is occupied by the recognition helix. As judged by hydroxyl radical footprinting and the missing nucleoside experiment, it appears that interaction of the Ubx recognition helix with the DNA major groove is reduced. Replacing the Dfd N-terminal arm with that of Ubx does not elicit a complete interchange of the DNA-binding mode. Although the position of the chimera relative to DNA, as judged by hydroxyl radical footprinting, is similar to that of the Dfd HD, the missing nucleoside and methylation interference patterns resemble those of the Ubx HD. Repositioning of amino acid side-chains without wholesale structural alteration in the polypeptide appears to occur as a function of N-terminal arm identity and DNA-binding site sequence. Complete interchange of binding modes was achieved only by replacement of the Dfd N-terminal arm and the recognition helix plus 13 carboxyl-terminal residues with the corresponding residues of Ubx. The position of the N-terminal arm in the DNA minor groove appears to differ in a manner that depends on the two base-pair differences between the Dfd and Ubx-optimal-binding sites. Thus, N-terminal arm position dictates the binding mode and the interaction of the recognition helix with nucleosides in the major groove.     

Sequence-specific deoxyribonucleic acid (DNA) recognition by steroidogenic factor 1: a helix at the carboxy terminus of the DNA binding domain is necessary for complex stability

Steroidogenic factor 1 (SF1) is a member of the NR5A subfamily of nuclear hormone receptors and is considered a master regulator of reproduction because it regulates a number of genes encoding reproductive hormones and enzymes involved in steroid hormone biosynthesis. Like other NR5A members, SF1 harbors a highly conserved approximately 30-residue segment called the FTZ-F1 box C-terminal to the core DNA binding domain (DBD) common to all nuclear receptors and binds to 9-bp DNA sequences as a monomer. Here we describe the solution structure of the SF1 DBD in complex with an atypical sequence in the proximal promoter region of the inhibin-alpha gene that encodes a subunit of a reproductive hormone. SF1 forms a specific complex with the DNA through a bipartite motif binding to the major and minor grooves through the core DBD and the N-terminal segment of the FTZ-F1 box, respectively, in a manner previously described for two other monomeric receptors, nerve growth factor-induced-B and estrogen-related receptor 2. However, unlike these receptors, SF1 harbors a helix in the C-terminal segment of the FTZ-F1 box that interacts with both the core DBD and DNA and serves as an important determinant of stability of the complex. We propose that the FTZ-F1 helix along with the core DBD serves as a platform for interactions with coactivators and other DNA-bound factors in the vicinity.     

Crystal structure of an engrailed homeodomain-DNA complex at 2.8 A resolution: a framework for understanding homeodomain-DNA interactions

The crystal structure of a complex containing the engrailed homeodomain and a duplex DNA site has been determined at 2.8 A resolution and refined to a crystallographic R factor of 24.4%. In this complex, two separate regions of the 61 amino acid polypeptide contact a TAAT subsite. An N-terminal arm fits into the minor groove, and the side chains of Arg-3 and Arg-5 make contacts near the 5' end of this "core consensus" binding site. An alpha helix fits into the major groove, and the side chains of IIe-47 and Asn-51 contact base pairs near the 3' end of the TAAT site. This "recognition helix" is part of a structurally conserved helix-turn-helix unit, but these helices are longer than the corresponding helices in the lambda repressor, and the relationship between the helix-turn-helix unit and the DNA is significantly different.     

1H NMR studies of DNA recognition by the glucocorticoid receptor: complex of the DNA binding domain with a half-site response element.

The complex of the rat glucocorticoid receptor (GR) DNA binding domain (DBD) and half-site sequence of the consensus glucocorticoid response element (GRE) has been studied by two-dimensional 1H NMR spectroscopy. The DNA fragment is a 10 base-pair oligonucleotide, 5'd(GCTGTTCTGC)3'.5'd-(GCAGAACAGC)3', containing the stronger binding GRE half-site hexamer, with GC base pairs at each end. The 93-residue GR-DBD contains an 86-residue segment corresponding to residues 440-525 of the rat GR. Eleven NOE cross peaks between the protein and DNA have been identified, and changes in the chemical shift of the DNA protons upon complex formation have been analyzed. Using these protein-DNA contact points, it can be concluded that (i) the "recognition helix" formed by residues C460-E469 lies in the major groove of the DNA; (ii) the GR-DBD is oriented on the GRE half-site such that residues A477-D481, forming the so-called D-loop, are available for protein-protein interaction in the GR-DBD dimer on the intact consensus GRE; and (iii) the 5-methyl of the second thymine in the half-site and valine 462 interact, confirming indirect evidence [Truss et al. (1990) Proc. Natl. Acad. Sci. U.S.A. 87, 7180-7184; Mader et al. (1989) Nature 338, 271-274] that both play an important role in GR-DBD DNA binding. These findings are consistent with the model proposed by H?rd et al. [(1990) Science 249, 157-160] and the X-ray crystallographic complex structure determined by Luisi et al. [(1991) Nature 352, 497-505].     

Structure of Bacterial Transcription Factor SpoIIID and Evidence for a Novel Mode of DNA Binding

SpoIIID is evolutionarily conserved in endospore-forming bacteria, and it activates or represses many genes during sporulation of Bacillus subtilis. An SpoIIID monomer binds DNA with high affinity and moderate sequence specificity. In addition to a predicted helix-turn-helix motif, SpoIIID has a C-terminal basic region that contributes to DNA binding. The nuclear magnetic resonance (NMR) solution structure of SpoIIID in complex with DNA revealed that SpoIIID does indeed have a helix-turn-helix domain and that it has a novel C-terminal helical extension. Residues in both of these regions interact with DNA, based on the NMR data and on the effects on DNA binding in vitro of SpoIIID with single-alanine substitutions. These data, as well as sequence conservation in SpoIIID binding sites, were used for information-driven docking to model the SpoIIID-DNA complex. The modeling resulted in a single cluster of models in which the recognition helix of the helix-turn-helix domain interacts with the major groove of DNA, as expected. Interestingly, the C-terminal extension, which includes two helices connected by a kink, interacts with the adjacent minor groove of DNA in the models. This predicted novel mode of binding is proposed to explain how a monomer of SpoIIID achieves high-affinity DNA binding. Since SpoIIID is conserved only in endospore-forming bacteria, which include important pathogenic Bacilli and Clostridia, whose ability to sporulate contributes to their environmental persistence, the interaction of the C-terminal extension of SpoIIID with DNA is a potential target for development of sporulation inhibitors.     

Crystal structure of the CENP-B protein-DNA complex: the DNA-binding domains of CENP-B induce kinks in the CENP-B box DNA.

The human centromere protein B (CENP-B), one of the centromere components, specifically binds a 17 bp sequence (the CENP-B box), which appears in every other alpha-satellite repeat. In the present study, the crystal structure of the complex of the DNA-binding region (129 residues) of CENP-B and the CENP-B box DNA has been determined at 2.5 A resolution. The DNA-binding region forms two helix-turn-helix domains, which are bound to adjacent major grooves of the DNA. The DNA is kinked at the two recognition helix contact sites, and the DNA region between the kinks is straight. Among the major groove protein-bound DNAs, this 'kink-straight-kink' bend contrasts with ordinary 'round bends' (gradual bending between two protein contact sites). The larger kink (43 degrees ) is induced by a novel mechanism, 'phosphate bridging by an arginine-rich helix': the recognition helix with an arginine cluster is inserted perpendicularly into the major groove and bridges the groove through direct interactions with the phosphate groups. The overall bending angle is 59 degrees, which may be important for the centromere-specific chromatin structure.     

A helix-turn-helix structure unit in human centromere protein B (CENP-B).

CENP-B has been suggested to organize arrays of centromere satellite DNA into a higher order structure which then directs centromere formation and kinetochore assembly in mammalian chromosomes. The N-terminal portion of CENP-B is a 15 kDa DNA binding domain (DBD) consisting of two repeating units, RP1 and RP2. The DBD specifically binds to the CENP-B box sequence (17 bp) in centromere DNA. We determined the solution structure of human CENP-B DBD RP1 by multi-dimensional 1H, 13C and 15N NMR methods. The CENP-B DBD RP1 structure consists of four helices and has a helix-turn-helix structure. The overall folding is similar to those of some other eukaryotic DBDs, although significant sequence homology with these proteins was not found. The DBD of yeast RAP1, a telomere binding protein, is most similar to CENP-B DBD RP1. We studied the interaction between CENP-B DBD RP1 and the CENP-B box by the use of NMR chemical shift perturbation. The results suggest that CENP-B DBD RP1 interacts with one of the essential regions of the CENP-B box DNA, mainly at the N-terminal basic region, the N-terminal portion of helix 2 and helix 3.     

A model for the LexA repressor DNA complex

A structural model for the interaction of the LexA repressor DNA binding domain (DBD) with operator DNA is derived by means of Monte Carlo docking. Protein–DNA complexes were generated by docking the LexA repressor DBD NMR solution structure onto both rigid and bent B-DNA structures while giving energy bonuses for contacts in agreement with experimental data. In the resulting complexes, helix III of the LexA repressor DBD is located in the major groove of the DNA and residues Asn-41, Glu-44, and Glu-45 form specific hydrogen bonds with bases of the CTGT DNA sequence. Ser-39, Ala-42, and Asn-41 are involved in a hydrophobic interaction with the methyl group of the first thymine base. Residues in the loop region connecting the two β-sheet strands are involved in nonspecific contacts near the dyad axis of the operator. The contacts observed in the docked complexes cover the entire consensus CTGT half-site DNA operator, thus explaining the specificity of the LexA repressor for such sequences. In addition, a large number of nonspecific interactions between protein and DNA is observed. The agreement between the derived model for the LexA repressor DBD/DNA complex and experimental biochemical results is discussed. © 1995 Wiley-Liss, Inc.