Existing crystal structure data has indicated that 1 alpha,25-dihydroxyvitamin D-3 (1 alpha,25(OH)(2)D-3) and its analogues bind the ligand-binding pocket (LBP) of the human vitamin D receptor in a very similar fashion. Because docking of a ligand into the LBP is a more flexible process than crystallography can monitor, we analyzed 1 alpha,25(OH)(2)D-3, its 20-epi derivative MC1288, the two side-chain analogues Gemini and Ro43-83582 ( a hexafluoro-derivative) by molecular dynamics simulations in a complex with the vitamin D receptor ligand-binding domain and a co-activator peptide. Superimposition of the structures showed that the side chain of MC1288, the first side chain of the conformation II of Gemini, the second side chain of Ro43-83582 in conformation I and the first side chain of Ro43-83582 in conformation II take the same agonistic position as the side chain of 1 alpha,25(OH)(2)D-3. Compared with the LBP of the natural hormone MC1288 reduced the volume by 17%, and Gemini expanded it by 19%. The shrinking of the LBP of MC1288 and its expansion to accommodate the second side chain of Gemini or Ro43-83582 is the combined result of minor movements of more than 30 residues and major movements of a few critical amino acids. The agonist-selective recognition of anchoring OH groups by the conformational flexible residues Ala-303, Leu-309, and His-397 was confirmed by in vitro assays. In summary, variations in the volume of agonists lead to adaptations in the volume of the LBP and alternative contacts of anchoring OH-groups.