Risk assessment for molecular toxicity endpoints of environmental matrices may be a pressing issue. Here, we combined chemical analysis with species sensitivity distributions (SSD) and in silico docking for multi-species estrogen receptor mediated-risk assessment in water from Dongjiang River, China. The water contains high levels of phenolic endocrine-disrupting chemicals (PEDCs) and phthalic acid esters (PAEs). The concentration of ∑4PEDCs and ∑6PAEs ranged from 2202 to 3404 ng/L and 834-4368 ng/L, with an average of 3241 and 2215 ng/L, respectively. The SSD approach showed that 4-NP, BPA, E2 of PEDCs, and DBP, DOP, and DEHP could severely threaten the aquatic ecosystems, while most other target compounds posed low-to-medium risks. Moreover, binding affinities from molecular docking among PEDCs, PAEs, and estrogen receptors (ERα, Erβ, and GPER) were applied as toxic equivalency factors. Estrogen receptor-mediated risk suggested that PEDCs were the main contributors, containing 53.37-69.79% of total risk. They potentially pose more severe estrogen-receptor toxicity to zebrafish, turtles, and frogs. ERβ was the major contributor, followed by ERα and GPER. This study is the first attempt to assess the estrogen receptor-mediated risk of river water in multiple aquatic organisms. The in silico simulation approach could complement toxic effect evaluations in molecular endpoints.