Use of binding energy in comparative molecular field analysis of isoform selective estrogen receptor ligands
A diverse set of 30 estrogen receptor ligands whose relative binding affinities (RBA) with respect to 17?-estradiol were available in both isoforms of the nuclear estrogen receptor (ER?, ER?) were studied with a combination of comparative molecular field analysis (CoMFA) and binding energy calculations. The ligands were docked inside the ligand-binding domain (LBD) of both ER? and ER? utilizing the docking program Gold. The binding energy (?E) and corresponding non-bonded interactions (NB) of the subsequent proteinligand complexes were calculated in both the gas-phase and implicit aqueous solution using the generalized born surface area (GB/SA) model. A partial least-squares analysis of the calculated energies indicated that the NB(g) were sufficiently predictive in ER?, but performed poorly in ER?. Further analysis of the calculated energies by dissecting the ligands into two distinct classes, estrogen-like and heterocyclic, yielded more predictive models. In particular the ?E calculated in solution proved particularly predictive for the estrogen-like ligands in ER?. Finally the estrogen subtype selective nature RBA (ER?/ER?) of a test-set consisting of six of the original ligands was predicted. The combined CoMFA and non-bonded interaction energy model ranked correctly the ligands in order of increasing RBA (ER?/ER?), illustrating the utility of this method as a prescreening tool in the development of novel estrogen receptor subtype selective ligands.
binding energy estrogen comfa subtype selective