Structure, thermodynamics and cooperativity of the glucocorticoid receptor DNA-binding domain in complex with...

Molecular dynamics simulations and free energy perturbation calculations have been performed on the glucocorticoid receptor DNA-binding domain (GRDBD) in complex with three different half-site response elements (RE). These are the glucocorticoid receptor RE (GREH), the estrogen receptor RE (EREH) and an "intermediate" RE (GRE2H), where the base-pair 6A.T in GREH has been replaced by 6G.C. The molecular dynamics/free energy perturbation calculations reveal the same order of stability of the three complexes, i.e. GREH-GRDBD > GRE2H-GRDBD > EREH-GRDBD, as found experimentally, with reasonable quantitative agreement of the relative stabilities with the experiments. The more favourable entropy in associating GRDBD to GREH than to GRE2H could be explained by the additional non-polar surface that is removed from the solvent in the former complex formation. Part of this additional surface originates from the 6T methyl group in GREH. Sequence-dependent hydration patterns of the DNA major groove, which is interrupted by the GRDBD binding, could also contribute to the relatively more favourable entropy in forming GREH-GRDBD than GRE2H-GRDBD. The tenfold reduced cooperativity in the dimeric GRDBD binding to a partially palindromic ERE (palERE) compared to palGRE/GRE2 could be attributed to conformational changes in the residues of the "D-box" (Cys476 to Cys482) observed in the molecular dynamics simulations of EREH-GRDBD, but not in the other two complexes. Moreover, the conformations of the base-pairs at positions 5 and 6 are distorted (compared to free B-DNA) in GREH and GRE2H, due to the binding of GRDBD. Since the corresponding conformational distortions were not found in EREH-GRDBD, it is plausible that the distortions mediate cooperative binding through a propagation of the distortions to the other half site, thus facilitating binding of the second GRDBD.