Coumarin phytoestrogens, as one of the important classes of phytoestrogens, have been proved to play an important role in various fields of human life. In this study, molecular simulation method including molecular docking and molecular dynamics methods were performed to explore the various effects between four classical coumarin phytoestrogens (coumestrol, 4-methoxycoumestrol, psoralen and isopsoralen), and estrogen receptors (ERα, ERβ), respectively. The calculated results not only proved that the four coumarin phytoestrogens have weaker affinity than 17β-estradiol to both ERα, and ERβ, but also pointed out that the selective affinity for ERβ is greater than ERα. In addition, the binding mode indicated that the formation of hydrogen bond and hydrophobic interaction have an important effect on the stability of the complexes. Further, the calculation and decomposition of binding free energy explored the main contribution interactions to the total free energy.
Keywords: binding free energy; coumarin; estrogen receptor; molecular dynamics simulation; phytoestrogens.