“…First, the optimal orientation of metal–ligand complexes within a protein is not well reproduced by docking programs because of the complexity of the coordination geometries that transition metals can exploit. , Recent efforts have been devoted to the development of ad hoc programs for zinc-dependent proteins, as AutoDock4Zn and MpSDockZn, but their usage is limited to this class of proteins. Second, commonly used force field-based, empirical or knowledge-based scoring functions do not account for electronic effects, and consequently, they are not able to capture the nature of organometallic bonds, preventing a correct ranking of ligand poses. − More accurate binding energies from quantum mechanics (QM)-based scoring functions ,− seem to correlate with experimental data, but the application of these methodologies to high throughput screening is restricted by the computational cost. These limitations in the modeling of metal–ligand binding explain the general preference for experimental techniques compared to in silico studies for drug discovery of jumonji demethylases. − …”