Estrogen receptors exist as two subtypes ER alpha and ER beta, which are characterized by various distributions in human tissues and diverse transcription regulation. Ligands capable of selective ER beta activation show positive effects in treatment of such diseases as certain cancers, endometriosis, inflammatory diseases, and assist in maintaining cardiovascular and nervous system health. Thus far, there are no pharmaceutical drugs available acting on this target. In order to provide new treatment for such diseases, a new generation of selective estrogen receptor modulators is required. This remains an unsolved task due to several difficulties. It is known that minor modifications of ER agonists can influence the selectivity of their binding. The majority of designed ligands acting on ER possess chiral centers thus exist as stereoisomers. Unfortunately, not every spatial isomer is individually considered in experimental research. The molecular docking was applied to investigate the structural basis of diverse selectivity and binding affinity of selected estrogen receptor beta agonists. Docking simulations revealed that terminal aromatic rings positioned in the A- and D-ring regions are a factor that determines binding affinity of ER beta agonists. This positioning can be ascribed to the presence of two terminal hydroxyl groups, a rigid linker, and the introduction of aliphatic substituents. The side substituents of underlined molecular scaffold should adopt inside characterized cavities I and II in order to provide selectivity. The bulkiness, attachment to linker and stereochemistry of the substituents affect ER beta selectivity. These molecular features should be considered during search and design of new improved ER beta agonists.