Acting during phase II metabolism, sulfotransferases (SULTs) serve detoxification by transforming a broad spectrum of compounds from pharmaceutical, nutritional, or environmental sources into more easily excretable metabolites. However, SULT activity has also been shown to promote formation of reactive metabolites that may have genotoxic effects. SULT subtype 1E1 (SULT1E1) was identified as a key player in estrogen homeostasis, which is involved in many physiological processes and the pathogenesis of breast and endometrial cancer. The development of an in silico prediction model for SULT1E1 ligands would therefore support the development of metabolically inert drugs and help to assess health risks related to hormonal imbalances. Here, we report on a novel approach to develop a model that enables prediction of substrates and inhibitors of SULT1E1. Molecular dynamics simulations were performed to investigate enzyme flexibility and sample protein conformations. Pharmacophores were developed that served as a cornerstone of the model, and machine learning techniques were applied for prediction refinement. The prediction model was used to screen the DrugBank (a database of experimental and approved drugs): 28% of the predicted hits were reported in literature as ligands of SULT1E1. From the remaining hits, a selection of nine molecules was subjected to biochemical assay validation and experimental results were in accordance with the in silico prediction of SULT1E1 inhibitors and substrates, thus affirming our prediction hypotheses.Metabolism plays an important role in drug discovery and appropriate metabolic profiles of new drug candidates should be taken into consideration during the process of drug development. Acting in phase I metabolism, the enzyme family of cytochrome P450 is estimated to transform about 75% of marketed drugs (1) and numerous in silico approaches for the prediction of cytochrome P450-mediated metabolism have emerged to date (2, 3). Although the majority of metabolism prediction studies focuses on phase I, the significance of phase II metabolism is generally underestimated (4) and to this day, computer-based models for the prediction of phase II metabolism remain scarce (5).Among the predominant phase II enzyme families are the soluble sulfotransferases that form a gene superfamily termed SULT. These enzymes regulate the sulfonation of smaller molecules such as endogenous hormones, neurotransmitters, and xenobiotic substances from pharmaceutical, nutritional, or environmental sources. Based on sequence similarity, functional human SULTs are divided into two main families (SULT1 and SULT2) and further into subfamilies that exhibit individual, but somewhat overlapping substrate specificities (6). Influencing the level of female sex hormones (estrogens), SULT subtype 1E1 (SULT1E1) shows a specific substrate preference for physiological estrogenic compounds (K m ϭ 5 nM for estradiol (7)) and has been extensively investigated in experimental studies.In general, sulfonation reactions in which a sulfona...