“…The effect of the rest of the protein environment was considered with the CPCM implicit solvation model using a dielectric constant of ε = 4.0, as suggested by Himo and co-workers [28], and a dielectric constant of ε = 78.4 for the aqueous solution, all in line with our previous reports [13], where we also demonstrated that a potential increase in the former dielectric constant to ε = 20.0 lowers the accuracy of the obtained results and even predicts wrong trends among ligands. Additionally, in our experience, such a truncated cluster-continuum model of the entire protein turned very useful in rationalizing various aspects of the catalytic activity [29], selectivity [30] and inhibition [31] of the monoamine oxidase family of enzymes, and is broadly used by different groups to decipher various biological phenomena [32][33][34][35][36], which justifies its use here. We note in passing that, despite its practical usefulness, the proposed value of ε = 4 is, in some cases, apparently too small to prevent the proton transfer from protonated to anionic residues, which then occurred spontaneously during the geometry optimization, provided the involved pair is in close vicinity.…”