Computational study of evolutionary selection pressure on rainbow trout estrogen receptors

Molecular dynamics simulations were used to determine the binding affinities between the hormone 17-estradiol (E2) and different estrogen receptor (ER) isoforms in the rainbow trout, Oncorhynchus mykiss. Previous phylogenetic analysis indicates that a whole genome duplication prior to the divergence of ray-finned fish led to two distinct ER isoforms, ER and ER, and the recent whole genome duplication in the ancestral salmonid created two ER isoforms, ER and ER. The objective of our computational studies is to provide insight into the underlying evolutionary pressures on these isoforms. For the ER subtype our results show that E2 binds preferentially to ER over ER. Tests of lineage specific N/S ratios indicate that the ligand binding domain of the ER gene is evolving under relaxed selection relative to all other ER genes. Comparison with the highly conserved DNA binding domain suggests that ER may be undergoing neofunctionalization possibly by binding to another ligand. By contrast, both ER and ER bind similarly to E2 and the best fitting model of selection indicates that the ligand binding domain of all ER genes are evolving under the same level of purifying selection, comparable to ER.