A Case Study on the Use of Binding Free Energies to Screen Inhibitors of Human Carbonic Anhydrase II.

We present in this article a case study on the thermodynamics of binding to human carbonic anhydrase II (HCA II) by three well-known inhibitors, viz. (a) acetazolamide (AZM) that directly binds to the catalytic Zn(II) ion at the active site, (b) non-zinc binding 6-hydroxy-2-thioxocoumarin (FC5) (c) 2-[(S)-benzylsulfinyl]benzoic acid (3G1). In each case, the crystal structure or its analogue of inhibitor-bound HCA II has been used to perform classical molecular dynamics (MD) simulation in water till 1 µ s ${1\hskip0.33em\mu s}$ . AZM and FC5 are found to undergo repeated binding and unbinding with markedly different dynamics from the partially buried, substrate-binding hydrophobic pocket near the active site. 3G1, on the other hand, is found to remain mostly at its crystallographic binding site occluded from the active site of HCA II. The associated binding free energies ( Δ G b i n d , s o l v ${{\rm \Delta }{G}_{bind,solv}}$ ) have been computed using the known MM/GBSA method and compared to the available experimental data. Our results show that Δ G b i n d , s o l v ${{\rm \Delta }{G}_{bind,solv}}$ encounters several issues including limited sampling of multiple binding sites and incorrect prediction of the affinity of the chosen ligands. Possible use of the simulation results in further construction of Markov state models is also discussed.

In the footsteps of an inhibitor unbinding from the active site of human carbonic anhydrase II.

The crystal structure of human carbonic anhydrase (HCA) II bound to an inhibitor molecule, 6-hydroxy-2-thioxocoumarin (FC5), shows FC5 to be located in a hydrophobic pocket at the active site. The present work employs classical molecular dynamics (MD) simulation to follow the FC5 molecule for 1 µs as it unbinds from its binding location, adopts the path of substrate/product diffusion (path 1) to leave the active site at around 75 ns. It is then found to undergo repeated binding and unbinding at different locations on the surface of the enzyme in water. Several transient excursions through different regions of the enzyme are also observed prior to its exit from the active site. These transient paths are combined with functionally relevant cavities/channels to enlist five additional pathways (path 2-6). Pathways 1-6 are subsequently explored using steered MD and umbrella sampling simulations. A free energy barrier of 0.969 kcal mol-1 is encountered along path 1, while barriers in the range of 0.57-2.84 kcal mol-1 are obtained along paths 2, 4 and 5. We also analyze in detail the interaction between FC5 and the enzyme along each path as the former leaves the active site of HCA II. Our results indicate path 1 to be the major exit pathway for FC5, although competing contributions may also come from the paths 2, 4 and 5.Communicated by Ramaswamy H. Sarma.

Structure and Dynamics of the Isozymes II and IX of Human Carbonic Anhydrase.

Human carbonic anhydrases (HCAs) are responsible for the pH control and sensing in our body and constitute key components in the central pH paradigm connected to cancer therapeutics. However, little or no molecular level studies are available on the pH-dependent stability and functional dynamics of the known isozymes of HCA. The main objective of this Article is to report the first bench-marking study on the structure and dynamics of the two most efficient isozymes, HCA II and IX, at neutral pH using classical molecular dynamics (MD) and constant pH MD (CpHMD) simulations combined with umbrella sampling, transition path sampling, and Markov state models. Starting from the known crystal structures of HCA II and the monomeric catalytic domain of HCA IX (labeled as HCA IX-c), we have generated classical MD and CpHMD trajectories (of length 1 µs each). In all cases, the overall stability, RMSD, and secondary structure segments of the two isozymes are found to be quite similar. Functionally important dynamics of these two enzymes have been probed in terms of active site hydration, coordination of the Zn(II) ion to a transient excess water, and the formation of putative proton transfer paths. The most important difference between the two isozymes is observed for the side-chain fluctuations of His-64 that is expected to shuttle an excess proton out of the active site as a part of the rate-determining intramolecular proton transfer reaction. The relative stability of the stable inward and outward conformations of the His-64 side-chain and the underlying free energy surfaces are found to depend strongly on the isozyme. In each case, a lower free energy barrier is detected between predominantly inward conformations from predominantly outward ones when simulated under constant pH conditions. The kinetic rate constants of interconversion between different free energy basins are found to span 107-108 s-1 with faster conformational transitions predicted at constant pH condition. The estimated rate constants and free energies are expected to validate if the fluctuation of the His-64 side-chain in HCA IX may have a significance similar to that known in the multistep catalytic cycle of HCA II.