Homology-modeled ligand-binding domains of medaka estrogen receptors and androgen receptors: A model system for the study of reproduction

Estrogen and androgen and their receptors play critical roles in physiological processes such as sexual differentiation and development. Using the available structural models for the human estrogen receptors alpha and beta and androgen receptor as templates, we designed in silico agonist and antagonist models of medaka estrogen receptor (meER) alpha, beta-1, and beta-2, and androgen receptor (meAR) alpha and beta. Using these models, we studied (1) the structural relationship between the ligand-binding domains (LBDs) of ERs and ARs of human and medaka, and (2) whether medaka ER and AR can be potential models for studying the ligand-binding activities of various agonists and antagonists of these receptors by docking analysis. A high level of conservation was observed between the sequences of the ligand-binding domains of meERα and huERα, meERβ1 and huERβ, meERβ2, and huERβ with 62.8%, 66.4%, and 65.1% identity, respectively. The sequence conservation between meARα and huAR, meARβ, and huAR was found with 70.1% and 61.0% of identity, respectively. Thirty-three selected endocrine disrupting chemicals (EDCs), including both agonists and antagonists, were docked into the LBD of ER and AR, and the corresponding docking score for medaka models and human templates were calculated. In order to confirm the conservation of the overall geometry and the binding pocket, the backbone root mean square deviation (RMSD) for Cα atoms was derived from the structure superposition of all 10 medaka homology models to the six human templates. Our results suggested conformational conservation between the ERs and ARs of medaka and human, Thus, medaka could be highly useful as a model system for studies involving estrogen and androgen interaction with their receptors.