Farnesoid X receptor (FXR, NR1H4) is a member of a nuclear receptor superfamily, which plays important roles in bile acid homeostasis, lipoprotein and glucose metabolism, and hepatic regeneration. GW4064 is a potent and selective FXR agonist and has become a tool compound to probe the physiological functions of FXR. Until now, the mechanism of GW4064 entering and leaving the FXR pocket is still poorly understood. Here, we report a computational study of GW4064 unbinding pathways from FXR by using several molecular dynamics (MD) simulation techniques. Based on the crystal structure of FXR in complex with GW4064, conventional MD was first used to refine the binding and check the stability of GW4064 in the FXR pocket. Random acceleration MD simulations were then performed to explore the possible unbinding pathways of GW4064 from FXR. Four main pathway clusters were found, among which three subpathways, namely Paths 2A, 2B, and 1B, were observed most frequently. Multiple steered MD simulations were further employed to estimate the maximum rupture force and the sum of the forces and to characterize the intermediate states of the ligand unbinding process. By comparing the average force profiles and structural changes, Paths 2A and 2B were identified to be the most favorable unbinding pathways. The former is located between the H1-H2 loop and the H5-H6 loop, and the latter is located in the cleft formed by the H5-H6 loop, H6, and H7. Moreover, the residues lining the pathways were analyzed for their roles in ligand unbinding. Based on our results, the possible structural modification strategies on GW4064 were also proposed.