Histone deacetylases (HDACs) are responsible for the removal of acetyl groups from histones, resulting in gene silencing. Overexpression of HDACs is associated with cancer, and their inhibitors are of particular interest as chemotherapeutics. However, HDACs remain a target of mechanistic debate. HDAC class 8 is the most studied HDAC, and of particular importance due to its human oncological relevance. HDAC8 has traditionally been considered to be a Zn-dependent enzyme. However, recent experimental assays have challenged this assumption and shown that HDAC8 is catalytically active with a variety of different metals, and that it may be a Fe-dependent enzyme in vivo. We studied two opposing mechanisms utilizing a series of divalent metal ions in physiological abundance (Zn2+, Fe2+, Co2+, Mn2+, Ni2+, and Mg2+). Extensive sampling of the entire protein with different bound metals was done with the mixed quantum-classical QM/DMD method. Density functional theory (DFT) on an unusually large cluster model was used to describe the active site and reaction mechanism. We have found that the reaction profile of HDAC8 is similar among all metals tested, and follows one of the previously published mechanisms, but the rate-determining step is different from the one previously claimed. We further provide a scheme for estimating the metal binding affinities to the protein. We use the quantum theory of atoms in molecules (QTAIM) to understand the different binding affinities for each metal in HDAC8 as well as the ability of each metal to bind and properly orient the substrate for deacetylation. The combination of this data with the catalytic rate constants is required to reproduce the experimentally observed trend in metal-depending performance. We predict Co2+ and Zn2+ to be the most active metals in HDAC8, followed by Fe2+, and Mn2+ and Mg2+ to be the least active.