The binding to the T cell receptor of wild-type and variant HTLV-1 Tax peptide complexed to the major histocompatibility complex has been investigated by means of molecular dynamics simulations. The binding free energy difference is calculated using the molecular mechanics Poisson-Boltzmann surface area and linear interaction energy methods. These methods extract useful information on the binding energetics from simulations of the physical states of the ligands, which are more computationally expedient than the commonly used thermodynamic integration method. The successful reproduction of the relative binding free energies shows that these methods can be useful for free energy calculations and the rational design of drugs and vaccines.