There are numerous quantum chemical methods available to study metal-ligand interaction. Such as Density functional theory (DFT), wave function theory (WFT). DFT is very popular in this category. Free energy perturbation (FEP) method is often used to study the energetics of metal-drug binding with very good accuracy. If the system size is too big, you can also opt for molecular mechanics (MM) or opt for semi empirical methods. The optimized choice depends on your needs and computational expensiveness of the calculations.
For DFT and WFT, you can use ORCA (free licensed software), TURBOMOLE, GAUSSIAN, ADF, MOLCAS, MOLPRO package. For more details, please follow the following link: