A great number of sesquiterpenes are reported in the available literature as good antileishmanial leads. However, their mode of action at the molecular level has not been elucidated. The lack of molecular studies could be considered an impediment for studies seeking to improve sesquiterpene-based drug design. The present in silico study allows us to make important observations about the molecular details of the binding modes of a set of antileishmanial sesquiterpenes against four drug-enzyme targets [pteridine reductase-1 (PTR1), N-myristoyl transferase (NMT), cysteine synthase (CS), trypanothione synthetase (TryS)]. Through molecular docking it was found that two sesquiterpene coumarins are promising leads for the PTR1 and TryS inhibition purposes, and some xanthanolides also exhibited better affinity towards PTR1 and CS binding. In addition, the affinity values were clustered by Principal Component Analysis and drug-like properties were analyzed for the strongest-docking sesquiterpenes. The results are an excellent starting point for future studies of structural optimization of this kind of compounds.