Molecular dynamics simulations (MD) have been performed on variant crystal and NMR-derived structures of the glucocorticoid receptor DNA-binding domain (GR DBD). A loop region five residues long, the so-called D-box, exhibits significant flexibility, and transient perturbations of the tetrahedral geometry of two structurally important Cys4 zinc finger are seen, coupled to conformational changes in the D-box. In some cases, one of the Cys ligands to zinc exchanges with water, although no global distortion of the protein structure is observed. Thus, from MD simulation, dynamics of the D-box could partly be explained by solvent effects in conjunction with structural reformation of the zinc finger.