Despite the availability of high resolution crystal structure of the catalytic domain of human carbonic anhydrase IX (denoted as HCA IX-c), its solution structure and functional dynamics are mostly unknown. In this article, we report classical molecular dynamics (MD) simulation studies on the structure and non-reactive dynamics of dimeric HCA IX-c (d-HCA IX-c, with chains A and B) in water and compare them to each chain simulated in water as a monomeric solute (mA/mB). While the average solution structures appear to be largely similar in all systems, d-HCA IX-c shows significant deviations from mA and mB in terms of non-reactive reorganization of the active site needed to initiate the rate determining proton transfer step between a zinc-bound water and one of the sidechain N-atoms of a distal histidine residue (His-64). First, 2−6 water molecules are found to form fluctuating hydrogen bonded networks that would serve as proton paths across the active site. A minimally frustrated core of each monomer in the dimeric protein is found to be associated with a significantly higher likelihood of forming complete proton transfer pathways mediated by 2 water molecules synchronously at the two active sites hosted by the monomeric chains. The relative stability and kinetics of transition between two major sidechain rotamers of His-64 are also affected by dimerization. The sidechain rotation of His-64 is found to involve at least one transition as slow as 105s−1 in chain A of d-HCA IX-c thereby rendering it slower than the actual rate determining proton transfer step. This observation indicates an important deviation from the widely accepted mechanism of catalysis by HCAs.
Read full abstract