Neutron Diffraction of Acetazolamide-Bound Human Carbonic Anhydrase II Reveals Atomic Details of Drug Binding
S. Zoë Fisher,† Mayank Aggarwal,‡ Andrey Y. Kovalevsky,† David N. Silverman,§ and Robert McKenna*,‡
† ‡

Bioscience Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States Department of Biochemistry and Molecular Biology, University of Florida, P.O. Box 100245, Gainesville, Florida 32610, United States § Department of Pharmacology and Therapeutics, University of Florida, P.O. Box 100247, Gainesville, Florida 32610, United States
S * Supporting Information

ABSTRACT: Carbonic anhydrases (CAs) catalyze the hydration of CO2 forming HCO3− and a proton, an important reaction for many physiological processes including respiration, fluid secretion, and pH regulation. As such, CA isoforms are prominent clinical targets for treating various diseases. The clinically used acetazolamide (AZM) is a sulfonamide that binds with high affinity to human CA isoform II (HCA II). There are several X-ray structures available of AZM bound to various CA isoforms, but these complexes do not show the charged state of AZM or the hydrogen atom positions of the protein and solvent. Neutron diffraction is a useful technique for directly observing H atoms and the mapping of H-bonding networks that can greatlycontribute to rational drug design. To this end, the neutron structure of H/D exchanged HCA II crystals in complex with AZM was determined. The structure reveals the molecular details of AZM binding and the charged state of the bound drug. This represents the first determined neutron structure of a clinically used drug bound to its target.

arbonic anhydrase (CA) is a ubiquitous metalloenzyme found in all kingdoms of life, from plants to humans, catalyzing CO2 hydration (and bicarbonate dehydration). Humans have 15 isoforms that are expressed in diverse tissues. These isoforms can be cytosolic, transmembrane, or membrane-bound.1 Human CA II (HCA II) is a monomeric 29 kDa soluble cytosolic isoform that contains zinc in the active site. HCA II is found predominantly in red blood cells and has the fastest catalytic turnover (106 s−1) among all the CAs characterized to date. Most CAs use a zinc-hydroxide (Zn− OH−) mechanism to reversibly convert CO2 to HCO3− and a proton, as illustrated in reactions 1 and 2:2
EZnOH− + CO2 EZnHCO3− EZnH 2O + HCO3−
(1)

C

EZnH 2O + B EZnOH + BH


(2)

CA inhibitors (CAIs) have been used for decades as diuretics and antiglaucoma drugs. CAIs can be broadly classified into two groups based on their binding mechanism: (a) binding to the tetrahedral configuration of the Zn metal center, thus displacing the Zn-boundsolvent (e.g., sulfonamides), or (b) by forming a
© 2012 American Chemical Society

trigonal-bipyramidal species through expansion of the metal coordination geometry without displacing the Zn-bound solvent (e.g., cyanates). The use of CAIs has now been expanded to include treatments for convulsions, obesity, and cancer, as well as being developed for use as diagnostic tools.3,4 In cancer, CA IX is believed to be responsible for the acidification of the extracellular matrix (ECM) when tumors become hypoxic.5 Due to the sequence conservation among CA isoforms in humans, clinically used CAIs target many isoforms. It should be noted that HCA II inhibition dominates as its expression is widespread, and HCA II has one of the highest catalytic rates. Most structure-based CAI designs have focused on HCA II Xray crystallographic studies. There are numerous crystal structures in the Protein Data Bank (www.rcsb.pdb.org) of CAs bound to a variety of CAIs: the sulfonamides, their bioisosteres, small molecule anions, phenols, coumarins, and polyamines.6,7 Sulfonamides are still of significant interest with more than 30 derivatives being used clinically.8 The primary sequence conservation among CA isoforms causes cross-reactivity of CAIs, and this necessitates the development of isoform-specific drugs.4 A detailed understanding of the water patterns and H-bonding in theactive site provided by this neutron diffraction study gives a new avenue for the rational structure-based drug design effort. Acetazolamide (AZM) dissolved in water has three possible protonation states with two associated pKa values (7.2 and 8.7) that are relevant to physiological pH (Figure 1). It is thought that any of these forms can bind to HCA II.9,10 There have been no crystallographic studies to definitively observe the charged state of AZM bound to HCA II, and it is unknown which form dominates in the crystalline state. This is despite the availability of very high-resolution (1.1 Šresolution) X-ray data.11 However, results from past 15N NMR studies have shown that it is the sulfonamido anion (Form 3 in Figure 1) that binds to HCA II in solution.12 There is limited explicit information available about the Hbonding interactions between enzymes and inhibitors/drugs and the role of solvent orientations. Besides the current neutron study of HCA II:AZM, there are only two other reports of a similar nature: Escherichia coli dihydrofolate
Received: July 12, 2012 Published: August 28, 2012
14726
dx.doi.org/10.1021/ja3068098 | J. Am. Chem. Soc. 2012, 134, 14726−14729


Journal of the American Chemical Society

Communication

Figure 2. Stick representation of AZM bound to HCA II. Zinc is shown as a magenta sphere, D atoms are in cyan (indicated byarrows), and H atoms in white (positions calculated). The nuclear density map is shown in yellow and is contoured at 1.5 σ, and the electron density map is shown in blue and is contoured at 2.0 σ.

Figure 1. Ionization and pKa of AZM in water. These different charged states represent the relevant forms that may bind HCA II under physiological conditions. Figure was drawn with ChemDraw (PerkinElmer) and adapted from Supuran and Winum, 2009.10

reductase (DHFR):methotrexate and HIV protease:KNI272.13,14 While methotrexate is a widely used chemotherapy agent, the reported neutron structure is in complex with a bacterial DHFR homologue. Similarly, the inhibitor KNI-272 bound to HIV protease has no clinical significance as it is not a drug currently in use. To determine the charged state and binding mode of AZM, a single H/D exchanged crystal of HCA II, cocrystallized with AZM, was prepared. The room temperature (RT) time-of-flight neutron diffraction data set was collected on a large crystal (volume 2 mm3) at the neutron Protein Crystallography Station at Los Alamos National Laboratory. Neutron and X-ray diffraction data extended to a resolution of 2.0 and 1.6 Å, respectively. Details of sample preparation, data collection, and joint X-ray and neutron structure refinement are described in the Supporting Information. All data collection and model refinement statistics aregiven in Table S1. The AZM molecule was clearly observed in omit Fo−Fc electron density maps and after placement and refinement was found to bind in the same overall configuration, as reported previously.11 A series of omit Fo−Fc and 2Fo−Fc nuclear density maps reveals the positions of the two exchangeable D atoms of AZM (Figure 2), information not obtainable from the X-ray data alone. Neutron diffraction gives an advantage over conventional Xray diffraction techniques in that even at medium resolution (1.5−2.0 Å), it is the only direct method available for unambiguously visualizing H (or D) atomic positions. Neutrons scattering off atomic nuclei lead to H/D atoms being as visible as heavier N, C, or O atoms. In contrast, X-ray diffraction occurs from the electron cloud making it very challenging to observe light H/D atoms with any confidence. A systematic study by Gardberg et al. carried out on a limited set of X-ray and neutron structures revealed a trend that at subatomic resolution ( 1 Å), only a fraction of H atoms are observed in electron density maps. This is in dramatic contrast to the case for nuclear density maps where almost all H/D atoms are observed at medium resolution.15
14727

The joint X-ray and neutron structure refinement approach exploits the strengths of each technique: the use of X-ray data to refine “heavy” (non-H) atoms in proteins andcomplementing this with neutron data to refine “light” (H/D) atoms (Figure 2).16 Highly complementary joint X-ray and neutron studies give very accurate and elusive details regarding Hbonding, solvent molecule orientation, and the protonation states of residues.17,18 The neutron structure reported here clearly shows the charged state of AZM as well as H-bonding between HCA II and AZM, information not obtainable from numerous highresolution X-ray studies (Figure 2).11,19 Omit Fo−Fc nuclear density maps calculated without the two exchanged D atomsiƒ•one each on the acetoamido and sulfonamido groupsiƒ•clearly revealed that AZM was in the anionic form, with the negatively charged sulfonamido group coordinated to the zinc (form 3 in Figure 1; also Figures 2 and 3).

Figure 3. Stick representation of the neutron structure of the active site of HCA II:AZM (PDB ID 4g0c). Amino acid residues and waters are as labeled; zinc is shown as a magenta sphere; exchanged D atoms are shown in cyan; and unexchanged H atoms in AZM are in white. Hbonds as observed in the nuclear maps are indicated by dashed lines with distances as indicated.

Four H-bonded waters are displaced upon AZM binding (waters S1−4; Figure S1).17 It was noted that these waters serve as a chemical template for where the AZM binds, in that the O of the waters superimpose with the N and S atoms of AZM.dx.doi.org/10.1021/ja3068098 | J. Am. Chem. Soc. 2012, 134, 14726−14729


Journal of the American Chemical Society The lone pair of sulfonamide N is involved in a coordinating bond with the zinc of 2.4 Å distance (Figures 3 and 4). This leads to a H-bonding scheme where the single D atom on the sulfonamide group acts as a H-bond donor to Thr199, which in turn acts as a H-bond donor to Glu106.

Communication

Figure 5. Ball-and-stick representation of the AZM weak hydrophobic interactions within the active site of HCA II. D (cyan) and H (white; positions calculated) atoms. Zinc is shown as a magenta sphere.

Figure 4. Complementarity of electron and nuclear density maps. Stereo stick representation of the active site of HCA II with AZM bound (PDB ID 4g0c). D atoms on AZM are indicated by black arrows. (a) Nuclear density map is shown in yellow and is contoured at 1.5 σ. (b) Electron density map is shown in blue and is contoured at 2.0 σ. Residue labeling is as shown in Figure 3.

The −OD of Thr200 participates as a bifurcated H-bond donor to both W1120 and the carbonyl of Pro201. W1120 forms an H-bond bridge between AZM and HCA II, acting as H-bond acceptor from the protonated acetoamido group and as a donor to the backbone carbonyl of Pro201 (Figures 3 and 4). Figure 4 illustrates the complementarity between electron and nuclear density maps. The neutronscattering length for S is quite small at only 2.7 fm. This is evident in the nuclear density maps where there is no nuclear density for the two S atoms of AZM at a 1.5 σ contour level (Figure 4a). In contrast to this and as expected, there is very strong density for S in the electron density maps (Figure 4b). The neutron scattering magnitude for D is comparatively large at 6.7 fm, and this is evident by the excellent density observed for both D atoms. As expected, the terminal −CH3 group was not visible in the neutron densities as there is significant signal cancellation at 2.0 Å resolution due to the presence of nonexchangeable H atoms (negative 3.7 fm scattering length) in such groups. However, it is not visible in the electron density map either. Although, electron density for the terminal −CH3 is seen in the highest resolution structure (PDB ID: 3hs4)11 of HCA II bound to AZM. In this case, the higher resolution X-ray data was not more informative for a heavy atom C compared to the neutron data. Combining the previous X-ray11 and current neutron structures, a full picture emerges of AZM binding to HCA II (Figures 3 and 5). Besides coordination to the zinc and Hbonding to HCA II, AZM also interacts through two very weak hydrophobic interactions. The first is a type of −CH···π interaction ( 3.5 Å distance) between Leu198 and the thiadiazole ring of AZM. Theother is a somewhat distorted −CH···π interaction between the terminal −CH3 of AZM and Phe131. As weak as these interactions are, they contribute to the overall binding of AZM to HCA II. The KI for AZM
14728

inhibition of HCA II and HCA IX is 12 and 25 nM, respectively.20 A direct comparison of the HCA II and HCA IX active sites shows that the only differences are at residues Phe131 and Val135, with Val and Leu at these positions in HCA IX, respectively (Figure S2). Phe131 (in HCA II) is involved in very weak hydrophobic interactions with AZM and its absence in HCA IX could, in part, explain the small 2-fold difference in binding constants for AZM between HCA II and IX.4 In conclusion, we have solved the first neutron structure of HCA II in complex with a sulfonamide inhibitor. This is the first example of a clinically used inhibitor bound to a human enzyme target to be studied with neutron diffraction. This study adds explicit details about the charged form of the drug that binds to HCA II, H-bonding to the target, and water displacement when AZM binds. A comparison of the current neutron structure with AZM bound to another neutron structure of unbound HCA II (PDB ID: 3tmj)17 reveals that four water molecules in the active site (S1−4, Figure S1) that were bound to either the zinc or H-bonded to other water molecules in the active site are displaced uponbinding of AZM (Figure 3). This structure provides insightful observations regarding H-bonds and hydrophobic interactions that play a key role in binding of AZM to HCA II. Neutron diffraction is the only technique that can directly reveal these details and is expected to contribute greatly to structure-based drug design in the future.

a–  a–  a– 

S * Supporting Information

ASSOCIATED CONTENT

Experimental data and coordinates were deposited with PDB ID 4g0c. Additional figures and experimental details of crystallization, data collection, and structure refinement. Table S1 contains crystallographic data collection and refinement statistics. This material is available free of charge via the Internet at https://pubs.acs.org.

AUTHOR INFORMATION

Corresponding Author

rmckenna@ufl.edu
Notes

The authors declare no competing financial interest.

ACKNOWLEDGMENTS The Protein Crystallography Station is funded by Department of Energy’s Office of Science Biological and Environmental Research. S.Z.F. is partially funded by LDRD Early Career
dx.doi.org/10.1021/ja3068098 | J. Am. Chem. Soc. 2012, 134, 14726−14729


Journal of the American Chemical Society award no. 20110535ER. A.Y.K. is partially funded by LDRD ER no. 20120256ER. R.M. and D.N.S. are partially funded by National Institutes of Health grant GM25154. All the figures were made using PyMOL (DelanoScientific).
(1) Gilmour, K. Comp. Biochem. Physiol., Part A: Mol. Integr. Physiol. 2010, 157, 193−197. (2) Silverman, D. N.; McKenna, R. Acc. Chem. Res. 2007, 40, 669− 675. (3) Hynninen, P.; Parkkila, S.; Huhtala, H.; Pastorekova, S.; Pastorek, J.; Waheed, A.; Sly, W.; Tomas, E. APMIS 2012, 120, 117−129. (4) Carbonic Anhydraseiƒ•Its Inhibitors and Activators; Supuran, C. T., Scozzafava, A., Conway, J., Eds.; CRC Press: Boca Raton, FL, 2004; Vol. 1. (5) Ditte, P.; Dequiedt, F.; Svastova, E.; Hulikova, A.; OhradanovaRepic, A.; Zatovicova, M.; Csaderova, L.; Kopacek, J.; Supuran, C. T.; Pastorekova, S.; Pastorek, J. Cancer Res. 2011, 71, 7558−7567. (6) Di Fiore, A.; Maresca, A.; Alterio, V.; Supuran, C. T.; De Simone, G. Chem. Commun. (Camb.) (7) Alterio, V.; Monti, S. M.; Truppo, E.; Pedone, C.; Supuran, C. T.; De Simone, G. Org. Biomol. Chem. 2010, 8, 3528−3533. (8) Supuran, C. T. Nat. Rev. Drug Discovery 2008, 7, 168−181. (9) Coleman, J. E. Ann. Rev. Pharmacol. (10) Supuran, C. T., Winum, J.-Y. Drug Design of Zinc-Enzyme Inhibitors: Functional, Structural, and Disease Applications; John Wiley & Sons: Hoboken, NJ, 2009; pp 439−446. (11) Sippel, K.; Robbins, A. H.; Domsic, J.; Genis, C.; AgbandjeMcKenna, M.; McKenna, R. Acta Crystallogr. 2009, F65, 992−995. (12) Kanamori, K.; Roberts, J. D. Biochemistry 1983, 22, 2658−2664. (13)Bennett, B.; Langan, P.; Coates, L.; Mustyakimov, M.; Schoenborn, B.; Howell, E. E.; Dealwis, C. Proc. Natl. Acad. Sci. U.S.A. 2006, 103, 18493−18498. (14) Adachi, M.; Ohhara, T.; Kurihara, K.; Tamada, T.; Honjo, E.; Okazaki, N.; Arai, S.; Shoyama, Y.; Kimura, K.; Matsumura, H.; Sugyiama, S.; Adachi, H.; Takano, K.; Mori, Y.; Hidaka, K.; Kimura, T.; Hayashi, Y.; Kiso, Y.; Kuroki, R. Proc. Natl. Acad. Sci. U.S.A. 2009, 106, 4641−4646. (15) Gardberg, A. S.; Del Castillo, A. R.; Weiss, K. L.; Meilleur, F.; Blakeley, M. P.; Myles, D. A. Acta Crystallogr. 2010, D66, 558−567. (16) Adams, P. D.; Mustyakimov, M.; Afonine, P. V.; Langan, P. Acta Crystallogr. 2009, D65, 567−573. (17) Fisher, S. Z.; Kovalevsky, A. Y.; Mustyakimov, M.; Silverman, D. N.; McKenna, R.; Langan, P. Biochem. (18) Kovalevsky, A. Y.; Hanson, B. L.; Mason, S. A.; Yoshida, T.; Fisher, S. Z.; Mustyakimov, M.; Forsyth, V. T.; Blakeley, M. P.; Keen, D. A.; Langan, P. Angew. Chem., Int. Ed. 2011, 50, 7520−7523. (19) Lesburg, C. A.; Huang, C.; Christianson, D. W.; Fierke, C. A. Biochemistry 1997, 36, 15780−15791. (20) Vullo, D.; Franchi, M.; Gallori, E.; Pastorek, J.; Scozzafava, A.; Pastorekova, S.; Supuran, C. T. Bioorg. Med. Chem. Lett.

Communication

a– 

REFERENCES

14729

dx.doi.org/10.1021/ja3068098 | J. Am. Chem. Soc. 2012, 134, 14726−14729