Golgi α-mannosidase II (GM) is a pharmaceutical target for the design of inhibitors with anticancer activity. The known potent GM inhibitors undergo complex interactions with Zn2+ ions and the active-site amino acids, many of which contain ionisable functional groups. Herein, the physical insight into the ligand⋯receptor interactions has been provided based on energy decomposition techniques: SAPT (symmetry adapted perturbation theory) and FMO-PIEDA (fragment molecular orbital-pair interaction energy decomposition analysis) for a large GM active-site cluster. Protonation-dependent molecular recognition in Golgi α-mannosidase was demonstrated for five inhibitors, mannose, and its transition state. Zn2+ ion and Asp472 induce the key interactions with the deprotonated inhibitors (bearing an amino group in the neutral state), followed by Asp92 and Asp341. This interaction pattern is consistent for all the studied inhibitors and is similar to the interaction pattern of the enzyme native substrate – mannose. The interactions with Zn2+ ion become repulsive for the protonated states of the inhibitors (bearing an amino group with +1 charge) and the mannosyl transition state. The importance of Asp92 and Asp204 considerably increases, while the interactions with Asp472 and Asp341 are slightly modified. The interaction pattern for the protonated ligands seems to have an oxocarbenium transition state-like character, rather than a Michaelis complex of GM. The electrostatic interactions with amino acids coordinating zinc ion are of key importance for both the neutral and protonated states of the inhibitors. The ligand's diol group has a dual role as an electron donor, coordinating zinc ion, and as an electron acceptor, interacting with Asp92 and Asp472 via strong hydrogen bonds. This interaction pattern is an essential structural feature of the potent GM inhibitors, which is consistent with the experimental findings. Based on the calculations, either the protonated or deprotonated state of the ligand may be the active form of the GM inhibitor, exhibiting different interacting patterns.
Quantum mechanical (QM) calculations at the level of density-functional tight-binding are applied to a protein–DNA complex (PDB: 2o8b) consisting of 3763 atoms, averaging 100 snapshots from molecular dynamics simulations. A detailed comparison of QM and force field (Amber) results is presented. It is shown that, when solvent screening is taken into account, the contributions of the backbones are small, and the binding of nucleotides in the double helix is governed by the base–base interactions. On the other hand, the backbones can make a substantial contribution to the binding of amino acid residues to nucleotides and other residues. The effect of charge transfer on the interactions is also analyzed, revealing that the actual charge of nucleotides and amino acid residues can differ by as much as 6 and 8% from the formal integer charge, respectively. The effect of interactions on topological models (protein -residue networks) is elucidated.