Structural and dynamic differences of the estrogen receptor DNA-binding domain, binding as a dimer and as a monomer to DNA: molecular dynamics simulation studies.

Molecular dynamics (MD) simulations of the estrogen receptor DNA-binding domain (ERDBD) as a dimer in complex with its DNA response element (ERE) show a significant difference in both structure and dynamics, compared to a MD simulation of monomeric ERDBD bound to its half-site response element (EREH). The C-terminal zinc binding domain (ZnII), including a region (helix II) which is in a helical conformation in ERE-(ERDBD)2, is considerably more flexible in EREH-ERDBD than in the dimeric complex. In EREH-ERDBD, all helical hydrogen bonds in helix II are broken and the entire ZnII region is detached from a hydrogen bonding network that in ERE-(ERDBD)2 connects to other parts of the protein as well as to the DNA. The regions that become flexible in EREH-ERDBD are identical to the regions where the NMR solution structure of free ERDBD is poorly ordered. This strongly suggests that dimerisation of ERDBD is required for ordering of the ZnII region and that monomeric binding to DNA is not sufficient for the ordering. This contrasts to the glucocorticoid receptor DNA-binding domain (GRDBD) which has essentially the same mobility (uniform and limited), regardless of whether it is free as a monomer in solution, bound as a monomer to its half-site response element or in a dimeric complex with the full response element. The hydrogen bonding network that connects ZnII with other parts of the protein and to DNA is almost identical in ERDBD and GRDBD. However, in GRDBD there is also a serine (in the N-terminal zinc coordinating region) with a central role in this network, connecting to the ZnII region. This serine is replaced by a glycine in ERDBD and we suggest that this substitution is sufficient for destabilisation of the network, thus leading to a more flexible ZnII region, which becomes ordered first upon forming a complex with another ERDBD and DNA.