Computational Study of the Sulfonylated Amino Acid Hydroxamates Binding to the Zinc Ion within the Active Site of Carbonic Anhydrase

Abstract

The B3LYP/6–311+G(d,p) method and three ONIOM extrapolation methods ONI-OM (B3LYP/6–311+G(d,p): AM1); ONIOM(B3LYP/6–311+G(d,p): MNDO); ONIOM (B3LYP/6–311+G(d,p): HF/3-21G(d)) were used to characterize the complexes of Zn²⁺ cation with anionic sulfonylated amino acid hydroxamates (RSO₂NH-AA-CON(-)OH), possessing an unsubstituted RSO₂NH—amino acyl moiety. According to the R moiety we distinguish between pentafluorophenyl and 4-methoxyphenyl derivates. The amino acid hydroxamates included in the study were the Gly, Ala, and Leu derivates. Of the inhibitors investigated, the weakest zinc affinity exhibits the pentafluorophenyl derivate with Gly amino acid and the strongest affinity the 4-methoxyphenyl derivate with Leu amino acid. The inhibitors form bidentate coordination bonds with the zinc cation by means of the sulfonyl oxygen and the ionized hydroxamate nitrogen atoms, respectively. The zinc affinities computed using the B3LYP/6–311 +G(d,p)//HF/6–31 +G(d,p) method are in very good agreement with the full density functional theory (DFT) B3LYP/6–311+G(d,p)//B3LYP/6- 311+G(d,p) method and this method can be adopted to model larger complexes of inhibitors with the active site of carbonic anhydrase.