The farnesoid X receptor (FXR) is a crucial therapeutic target for treating non-alcoholic steatohepatitis (NASH). Although obeticholic acid (OCA) as a FXR agonist presents good efficacy, the safety data such as severe pruritus should be carefully considered. To discover new medications, we screen and choose the optimal compounds from ZINC15 database that may agonistically interact with FXR. We utilized the DS19 software to assist us in conducting the computer-aided structure based virtual screening to discover potential FXR agonists. After LibDock scores were determined by screening, their absorption, distribution, metabolism, excretion and toxicity predictions were examined. To determine the binding affinity between the chosen drugs and FXR, molecule docking was utilized. Molecular dynamics simulation was utilized to evaluate the stabilization of the ligand-FXR complex in its native environment. Higher binding affinity and stability with FXR were observed for ZINC000013374322 and ZINC000006036327, as two novel natural compounds, with lower rodent carcinogenicity, Ames mutagenicity, no hepatotoxicity and non-inhibitors of CYP2D6. They could stably exist in the environment, possess favorable potential energy and exert pharmacological effects at lower doses. Furthermore, ZINC000006036327 had lower skin irritancy and sensitization potential compared to OCA, also suggest the possibility of improved skin itching occurrence. ZINC000013374322 and ZINC000006036327 were found to be the best leading compounds to be FXR agonists. They are chosen as safe candidates for FXR target medicine, which play comparable pharmacological effects at lower doses.