Ligand selectivity of estrogen receptors by a molecular dynamics study

Estrogen receptors α (ERα) and β (ERβ) have different physiological functions and expression levels in different tissues. ERα and ERβ are highly homologous and have only two residue substitutions in the binding pocket. This high similarity at the active site stimulates the interests for discovering subtype selective ligands. In this study, molecular dynamics (MD) simulations combined with molecular mechanics generalized Born surface area (MM-GBSA) method have been carried out to analyze the basis of selectivity of three ligands (659, 818 and 041). The calculated binding free energies show that all the ligands bind more tightly to ERβ than to ERα. The dominant free energy components of selectivity for 659 are similar to that for 041, but different from that for 818. The decompositions of free energy contributions and structural analysis imply that there are eight residues primarily contributing to the selectivity for 659, five residues for 041, as well as two residues for 818. The structural analysis implies that two residue substitutions in binding packet cause the position of 659 in ERβ-659 complex to shift relative to that in ERα-659 complex and also cause the conformational changes of other residues in the binding pocket. The higher selectivity for 041 is mainly caused by three residues, Ile373 (Met421), His475 (His524) and Leu476 (Leu525).