Glucocerebrosidase (GCase, acid β-Glucosidase) hydrolyzes the sphingolipid glucosylceramide into glucose and ceramide. Mutations in this enzyme lead to a lipid metabolism disorder known as Gaucher disease. The design of competitive inhibitors of GCase is a promising field of research for the design of pharmacological chaperones as new therapeutic agents. Using a series of recently reported molecules with experimental binding affinities for GCase in the nanomolar to micromolar range, we here report an extensive theoretical analysis of their binding mode. On the basis of molecular docking, molecular dynamics, and binding free energy calculations using the linear interaction energy method (LIE), we provide details on the molecular interactions supporting ligand binding in the different families of compounds. The applicability of other computational approaches, such as the COMBINE methodology, is also investigated. The results show the robustness of the standard parametrization of the LIE method, which reproduces the experimental affinities with a mean unsigned error of 0.7 kcal/mol. Several structure-activity relationships are established using the computational models here provided, including the identification of hot spot residues in the binding site. The models derived are envisaged as important tools in ligand-design programs for GCase inhibitors.