: This article presents an overview of recent computational studies dedicated to the analysis of binding between galectins and small-molecule ligands. We first present a summary of the most popular simulation techniques adopted for calculating binding poses and binding energies and then discuss relevant examples reported in the literature for the three main classes of galectins (dimeric, tandem, and chimera). We show that simulation of galectin-ligand interactions is a mature field that has proven invaluable for completing and unraveling experimental observations. Future perspectives to further improve the accuracy and cost-effectiveness of existing computational approaches will involve the development of new schemes to account for solvation and entropy effects, which represent the main current limitations to the accuracy of computational results.