Quantum chemical study of silanediols as metal binding groups for metalloprotease inhibitors

DFT (B3LYP and M06L) as well as ab initio (MP2) methods with Dunning cc-pVnZ (n=2,3) basis sets are employed for the study of the binding ability of the new class of protease inhibitors, i.e., silanediols, in comparison to the well-known and well-studied class of inhibitors with hydroxamic functionality (HAM). Active sites of metalloproteases are modeled by [R₃M-OH₂]²⁺ complexes, where R stands for ammonia or imidazole molecules and M is a divalent cation, namely zinc, iron or nickel (in their different spin states). The inhibiting activity is estimated by calculating Gibbs free energies of the water displacement by metal binding groups (MBGs) according to: [R₃M-OH₂]²⁺ + MBG → [R₃M-MBG]²⁺ + H₂O. The binding energy of silanediol is only a few kcal mol(-1) inferior to that of HAM for zinc and iron complexes and is even slightly higher for the triplet state of the (NH₃)₃Ni²⁺ complex. For both MBGs studied in the ammonia model the binding ability is nearly the same, i.e., Fe²⁺(t) > Ni²⁺(t) > Fe²⁺(q) > Ni²⁺(s) > Zn²⁺. However, for the imidazole model the order is slightly different, i.e., Ni²⁺(t) > Fe²⁺(t) > Fe²⁺(q) > Ni²⁺(s) ≥ Zn²⁺. Equilibrium structures of the R₃Zn ²⁺ complexes with both HAM and silanediol are characterized by the monodentate binding, but the bidentate character of binding increases on going to iron and nickel complexes. Two types of intermediates of the water displacement reactions for [(NH₃)₃M-OH₂]²⁺ complexes were found which differ by the direction of the attack of the MBG. Hexacoordinated complexes exhibit bidentate bonding of MBGs and are lower in energy for M=Ni and Fe. For Zn penta- and hexacoordinated complexes have nearly the same energy. Intermediate complexes with imidazole ligands have only octahedral structures with bidentate bonding of both HAM and dimethylsilanediol molecules.