Email updates

Keep up to date with the latest news and content from BMC Biophysics and BioMed Central.

Open Access Research article

Intrinsic thermodynamics of ethoxzolamide inhibitor binding to human carbonic anhydrase XIII

Lina Baranauskienė and Daumantas Matulis*

Author Affiliations

Department of Biothermodynamics and Drug Design, Vilnius University Institute of Biotechnology, Graičiūno 8, Vilnius LT-02241, Lithuania

For all author emails, please log on.

BMC Biophysics 2012, 5:12  doi:10.1186/2046-1682-5-12

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/2046-1682/5/12


Received:29 November 2011
Accepted:1 May 2012
Published:7 June 2012

© 2012 Baranauskienė and Matulis; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License ( http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Human carbonic anhydrases (CAs) play crucial role in various physiological processes including carbon dioxide and hydrocarbon transport, acid homeostasis, biosynthetic reactions, and various pathological processes, especially tumor progression. Therefore, CAs are interesting targets for pharmaceutical research. The structure-activity relationships (SAR) of designed inhibitors require detailed thermodynamic and structural characterization of the binding reaction. Unfortunately, most publications list only the observed thermodynamic parameters that are significantly different from the intrinsic parameters. However, only intrinsic parameters could be used in the rational design and SAR of the novel compounds.

Results

Intrinsic binding parameters for several inhibitors, including ethoxzolamide, trifluoromethanesulfonamide, and acetazolamide, binding to recombinant human CA XIII isozyme were determined. The parameters were the intrinsic Gibbs free energy, enthalpy, entropy, and the heat capacity. They were determined by titration calorimetry and thermal shift assay in a wide pH and temperature range to dissect all linked protonation reaction contributions.

Conclusions

Precise determination of the inhibitor binding thermodynamics enabled correct intrinsic affinity and enthalpy ranking of the compounds and provided the means for SAR analysis of other rationally designed CA inhibitors.

Background

Carbonic anhydrases (CAs) are zinc-containing metalloenzymes that catalyze reversible hydration of carbon dioxide to the bicarbonate ion. The enzyme has been studied since 1930s and has become a classical model protein in biochemistry [1-5]. Recently, the interest in the enzyme has risen again. It is known now that CAs are present both in prokaryotes and eukaryotes, and are encoded by five gene families: the α-CAs (predominantly in vertebrates), the β-CAs (bacteria, algae, other), γ-CAs (archaea), δ-CAs and ζ-CAs (in some marine diatoms) [6-8]. All human CAs belong to the alpha class. There are 12 catalytically active human CA isoforms: I, II, III, VII, and XIII are cytosolic, IV, IX, XII, and XIV are associated with cell membrane, VA and VB occur in mitochondria, and VI is secreted. There are 3 CA-related proteins VIII, X, and XI that are devoid of catalytic activity due to the absence of zinc. Human CAs were reviewed in detail, including their widely different catalytic and kinetic properties, response to inhibitors, expression patterns, and application for drug design [8].

Human CAs are widely distributed in many tissues and organs. The CAs play a crucial role in CO2 and HCO3- transport, pH and CO2 homeostasis, electrolyte secretion, biosynthetic reactions, and tumor progression. Therefore the CAs became interesting targets for pharmaceutical research [9]. However, the main drawback of available compounds is insufficient selectivity towards particular CA isozyme.

Rational design of compounds with desired binding properties requires detailed investigation of the structure-activity relationships (SAR) of the newly designed compounds. Both the thermodynamics of binding and the structure of protein-ligand complex are required for detailed understanding of the reaction and search and rational design of drug-like molecules. It is important to go beyond the determination of the Ki and determine the full thermodynamic profile, especially the enthalpy of binding [10]. Here we describe the intrinsic thermodynamic analysis of well known CA inhibitors, ethoxzolamide (EZA), trifluoromethanesulfonamide (TFMSA), and acetazolamide (AZM) binding to human CA XIII isozyme.

Human CA XIII isozyme has been characterized [11]. The enzyme expression in human tissues showed that CA XIII is found in several organs including the thymus, kidney, submandibular gland, small intestine, and most notably in reproductive organs suggesting involvement in the fertilization process [7,11]. The CO2 hydration activity showed CA XIII to be a catalyst of medium efficiency [12]. Inhibition profiles demonstrated that CA XIII is similar to CA II [13]. The crystal structure of hCA XIII and its complex with acetazolamide has been solved providing crucial structural data for inhibitor SAR [12].

Manuscripts reporting novel CA inhibitors usually list only the observed thermodynamic parameters of binding or the observed inhibition constants (Ki) [8,14-16]. Such data is useful for the design of novel inhibitors. However, it should not be used for the SAR and correlations with structures. Proper SAR requires correlation of structures only with the intrinsic binding parameters that are dissected from linked reactions occurring simultaneously with the binding reaction. Most often, such linked reactions are protonation-deprotonation reactions occurring upon ligand binding. In the case of CAs, both the inhibitor, and the CA molecule may or may not exhibit a linked protonation reaction upon inhibitor binding. Observed and intrinsic binding parameters coincide only in some rare cases when the inhibitor is deprotonated and the zinc-linked hydroxide anion is protonated. Such situation is possible only with some CA isozymes and only with the inhibitors with the pKas significantly below 7. Most sulfonamide inhibitors have their pKas in the range between 8 and 11. Therefore, only small fraction of the compound is in the form that binds the protein. Despite the fact that the linked protonation reactions are largely worked-out, they are rarely used in the SAR of CA inhibitors [17-19].

By ‘intrinsic’ we mean the parameters that describe actual binding reaction. In the case of sulfonamide inhibitor binding to CAs, the actual binding components are the deprotonated sulfonamide and the CA with protonated Zn-bound water molecule. Only these two forms bind each other. However, in most cases, other forms are the most abundant in solution – sulfonamide is usually protonated and the CA contains unprotonated water molecule. The concentration of the forms that actually bind is small relative to other forms and it takes energy to convert them to active forms. The observed binding energy is therefore diminished and does not resemble true (intrinsic) energetic parameters.

Determination of the intrinsic binding constant (Kb, Gibbs free energy ΔbGintr), intrinsic enthalpy (ΔbHintr), intrinsic entropy (ΔbSintr), and the intrinsic heat capacity of binding (ΔbCp_intr) requires significant effort and a number of measurements at various pHs and temperatures. Here we use isothermal titration calorimetry (ITC) and thermal shift assay (TSA, also called ThermoFluor®, differential scanning fluorimetry) to measure inhibitor binding to CAs. ITC has been routinely used to measure ligand-protein binding thermodynamics [20]. However, it is not appropriate for determining weak (millimolar) or very tight dissociation constants (subnanomolar Kds require displacement ITC [21]) and consumes extensive amounts of protein and time. In contrast, TSA is a rapid screening method used in pharmaceutical industry for the identification of binders and requires lower amounts of protein [22-27]. The method is based on the protein melting temperature (Tm) shift that occurs upon ligand binding. The Tm is observed by following intrinsic or extrinsic fluorescence changes upon heat-induced protein unfolding. The employment of two techniques to determine binding reactions reduces the error of the measurements. Here we apply both techniques to determine the observed binding parameters and then estimate the intrinsic parameters that could be used for CA inhibitor SAR analysis.

Results

Ethoxzolamide and other inhibitor binding to hCA XIII by ITC

In order to draw any thermodynamic energetics-structure correlations, the experimentally determined binding parameters must be intrinsic. In other words, all binding-linked contributing reactions must be accounted for and dissected from the reaction of interest. The most common protein-ligand binding-linked reactions are protonation or deprotonation. Sulfonamide inhibitor binding reaction to carbonic anhydrase (CA) is linked to at least two reactions, the protonation of zinc-bound hydroxide anion, and deprotonation of sulfonamide inhibitor group.

Figure 1 shows the reactions linked to ethoxzolamide (EZA) binding to hCA XIII. First, the sulfonamide group of ethoxzolamide must be deprotonated. Second, the zinc-bound hydroxide must be protonated. Third, the protons needed for CA and inhibitor are absorbed or donated by the buffer of the system. TRIS buffer has a large enthalpy of protonation equal to −47.4 kJ/mol while phosphate exhibits a relatively small enthalpy of approximately −5.1 kJ/mol at 25°C. Therefore, in the presence of binding-linked protonation reactions, the observed enthalpies of binding will vary significantly in these buffers and may even change sign. Obviously, protonation contributions must be dissected from binding reactions in order to obtain the intrinsic binding values.

thumbnailFigure 1. The processes linked to the binding of EZA to hCA XIII. Two central reactions show the deprotonation of the inhibitor sulfonamide group and protonation of the zinc-bound hydroxide to make water molecule and consecutive binding of the two active species. Two top and bottom lines show linked buffer protonation reactions. The numbers give estimates of intrinsic enthalpies for each process multiplied by the number of linked protons (n) yielding the observed enthalpic contribution of each reaction at pH 7.0. The numbers n were estimated by fitting the enthalpy dependence on pH data. The intrinsic enthalpy of binding is shown next to the reaction and is equal to −42.1 kJ/mol. The observed enthalpies are estimated for both buffers in the right upper corner. These values are similar to experimentally observed values listed in Table 1. Zinc atom is shown as grey shaded sphere and the carbonic anhydrase is shown as CA.

A typical isothermal titration calorimetry (ITC) raw data curve of EZA binding to hCA XIII is shown in Figure 2a. The reaction is highly exothermic and the binding is tight as seen from the steepness of the curve. The large negative peaks are due to all binding and protonation reactions while the small peaks after the middle of titration are due to dilution effects. Figure 2b shows integrated ITC curves of EZA binding to hCA XIII at various pHs in sodium phosphate buffer. The observed enthalpies of binding vary in the range of −11 kJ/mol at pH 6.0 to −66 kJ/mol at pH 9.0.

thumbnailFigure 2. The binding of EZA and TFMSA to hCA XIII by isothermal titration calorimetry (ITC). (a) Raw ITC curve of EZA binding to hCA XIII at pH 7.0, TRIS chloride buffer, at 25°C. (b) Integrated ITC curves of EZA binding to hCA XIII in sodium phosphate buffer, 25°C, at several pHs: 6.5 (■), 7.5 (), and 8.5 (▲). (c) Integrated ITC curves of EZA binding to hCA XIII in TRIS chloride (◊) and sodium phosphate (▲) buffer, pH 8.5, at 25°C. (d) Observed integral enthalpies of EZA () and TFMSA (◆) binding to hCA XIII as a function of buffer deprotonation enthalpy, pH 7.0, at 25°C in the following buffers (in the order of increasing enthalpy of deprotonation, shown in the brackets, in kJ/mol): Pi (5.1), PIPES (11.5), MES (15.5), HEPES (21.0), MOPS (22.2), Imidazole (36.6), and TRIS (47.6). (e) Integral observed enthalpies of EZA binding to hCA XIII as a function of temperature. The slope is equal to the observed heat capacity of binding.

Similarly to the effect of pH, the effect of buffer is seen from Figure 2c where the observed enthalpy of binding is approximately 40% more exothermic in phosphate than in TRIS buffer. The observed enthalpies of binding vary linearly with the enthalpy of buffer protonation. Figure 2d shows EZA and trifluoromethanesulfonamide (TFMSA) binding at pH 7.0 to hCA XIII in various buffers: phosphate (5.1), PIPES (11.5), MES (15.5), HEPES (21.0), MOPS (22.2), imidazole (36.6), and TRIS (47.4) (in the order of increasing buffer protonation enthalpy in kJ/mol at 25°C in the brackets). Here, the slope of the dependence shows the net number of uptaken or released protons upon binding. The positive slope indicates that the protons are taken up from buffer upon ligand binding while the negative slope shows that the protons are given from the protein and /or ligand to the buffer. It is clear that EZA binding is linked to protonation reactions and the protons are released from the protein-ligand system to the buffer while TFMSA binding appears to be either not linked to protonation or the net number is approximately equal to zero.

Figure 2e shows the observed EZA binding enthalpy dependence on temperature at pH 7.0. The dependence is linear within the error margin of the enthalpy determination yielding the observed heat capacity of EZA binding equal to −980 J/(mol × K).

Table 1 and Figure 3 show the ITC-observed thermodynamic parameters of EZA binding to hCA XIII at 25°C as a function of pH in sodium phosphate and TRIS chloride buffers. The observed parameters are significantly buffer and pH-dependent. Such parameters are experimentally observed but cannot be used for structure-activity or structure-energetics correlations. The observed enthalpies of binding vary from approximately −10 to −70 kJ/mol for EZA and TFMSA binding to hCA XIII. The values for two buffers cross each other near pH 7. Detailed dissection of linked protonation reactions is necessary to obtain the intrinsic binding thermodynamic parameters.

Table 1. Observed thermodynamic parameters of EZA binding to hCA XIII in phosphate and TRIS buffers are listed as a function of pH as determined by ITC at 25°C

thumbnailFigure 3. Observed enthalpies of EZA (a) and TFMSA (b) binding to hCA XIII in TRIS chloride () and sodium phosphate (◯) buffer, 25°C, are shown as a function of pH. Datapoints represent observed enthalpies obtained by ITC, while the curves are fitted according to Eq. (4).

ITC is a very good method to determine intrinsic binding parameters and dissect the linked protonation reactions. However, it must be used with great care. The main limitation is that the method has rather narrow window of experimentally determinable binding constants. If protein concentration is 10 μM, and if c parameter (the product of binding constant and protein concentration) is to be between 5 and 500, then the accurate observed Kb_obs could be obtained only from 5 × 105 to 5 × 107 M-1. As seen in Table 1, the binding is tighter at many conditions. Therefore, an additional independent method is necessary to accurately determine the intrinsic binding constant Kb_intr, and the derived thermodynamic parameters such as the Gibbs free energy and the entropy of binding.

EZA and TFMSA binding to hCA XIII by TSA and DSC

A useful, simple and widely applicable method to determine ultratight binding constants is the thermal shift assay (TSA). The method measures the shift in protein melting temperature upon addition of a binding ligand. Figure 4a shows raw fluorescence curves as a function of temperature at various added EZA concentrations. The steep increase in fluorescence is due to hCA XIII unfolding and exposure of additional hydrophobic surfaces where the fluorescence of bound ANS increases. The curve midpoints represent the melting temperatures (Tm). These midpoints shift towards higher temperatures upon addition of ligand. Addition of 100 μM EZA shifted the melting temperature by more than 10°C. The melting temperatures can be plotted as a function of temperature in a form of a ligand-dosing curve. Both the fluorescence and ligand-dosing curves can be fit as previously described [28] to accurately determine the thermodynamic parameters of unfolding and the Gibbs free energies of binding. Figure 4b compares EZA and TFMSA dosing curves. The Tm shift is greater for EZA than TFMSA and thus the observed Kb_obs is greater for EZA at the experimental conditions. The dosing curves best fitted with the enthalpy of unfolding 460 ± 100 kJ/mol (at the temperature of the Tm equal to 59.2°C at pH 7.0).

thumbnailFigure 4. Thermal shift assay and DSC data of EZA and TFMSA binding to hCA XIII. (a) Thermal melting fluorescence curves of hCA XIII (concentration 10 μM) in sodium phosphate buffer, pH 7.0, in the presence of various EZA concentrations, 0 μM (◆), 12 μM (□), 26 μM (▲), and 89 μM (◊). Ligand addition shifted the protein melting temperatures upwards. (b) The dose response curves as a function of added concentration of EZA (◊) or TFMSA (■). Datapoints are obtained from fluorescence melting curves such as in panel (a), while the lines are fitted according to [27-29]. The addition of EZA shifted the Tms more than TFMSA at the same concentration. Therefore, the observed Kbs are greater for EZA than for TFMSA. The possible reasons for the bending of the EZA data have been previously discussed [29]. (c) DSC curves of hCA XIII (50 μM) in the absence (solid line) and presence of 500 μM EZA (dashed line) in sodium phosphate buffer at pH 7.0. The addition of EZA shifted the peak to higher temperatures (thus stabilized the protein) and the area under the peak (equal to the enthalpy of unfolding) has increased due to positive heat capacity of unfolding.

Similar Tm shift is observed by DSC (Figure 4c). Addition of EZA shifted the Tm by more than 10°C and the enthalpy of unfolding (the area under the DSC peak) increased due to the positive heat capacity of hCA XIII unfolding. In the absence of ligand, the hCA XIII exhibited the Tm equal to 58.6°C and the enthalpy of unfolding was equal to 573 ± 150 kJ/mol at pH 7.0. In the presence of 0.5 mM EZA the peak shifted to 72.4°C and the enthalpy of unfolding increased to 937 kJ/mol.

Dissection of binding-linked protonation reactions

Here we use the same formalism to dissect protonation reactions as previously described [17-19]. The intrinsic binding constant is related to the observed binding constant and the fractions of active binding species:

<a onClick="popup('http://www.biomedcentral.com/2046-1682/5/12/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/2046-1682/5/12/mathml/M1">View MathML</a>

(1)

Where the fractions deprotonated sulfonamide and protonated CA as a function of pH are:

<a onClick="popup('http://www.biomedcentral.com/2046-1682/5/12/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/2046-1682/5/12/mathml/M2">View MathML</a>

(2)

<a onClick="popup('http://www.biomedcentral.com/2046-1682/5/12/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/2046-1682/5/12/mathml/M3">View MathML</a>

(3)

The intrinsic enthalpy of binding has contributions from the observed binding enthalpy (Δb_obsH), enthalpies of inhibitor (Δb_proton_inhH), CA (Δb_proton_CAH), and the buffer protonation (Δb_proton_bufH):

<a onClick="popup('http://www.biomedcentral.com/2046-1682/5/12/mathml/M4','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/2046-1682/5/12/mathml/M4">View MathML</a>

(4)

Figure 5a shows the EZA and TFMSA -hCA XIII observed Kb_obs dependence on pH. The points are experimentally (by ITC or TSA) determined while the lines are fitted according to Eq. (1). Figure 5b shows the same data for EZA plotted as Gibbs free energy. In addition, the horizontal line shows the level of intrinsic Gibbs free energy of binding and two dashed lines show contributions to the observed Gibbs free energy by the two fractions (Eqs. 2 and 3). Figure 5c compares intrinsic Gibbs free energy of TFMSA binding with the observed values. It is notable, that for TFMSA one can measure experimentally the values that are similar to intrinsic, while EZA observed values are significantly different from the intrinsic ones.

thumbnailFigure 5. The dependence of ligand binding constant on pH. (a) Comparison of TFMSA (filled symbols) and EZA (open symbols) binding pH profiles. Curves show the fits according to Eq. (1) using parameters listed in the Tables  2 and 3. Datapoints show Kbs obtained by ITC in phosphate (diamonds) and TRIS (triangles) buffers and by TSA (squares). (b) The dissection of the Gibbs free energies of EZA binding to hCA XIII as a function of pH. Datapoints are calculated from panel (a). Solid bent line shows the fit according to Eq. (1). Solid straight line shows the position of the intrinsic Gibbs free energy of binding which is independent of pH. Dashed and dotted lines show the contributions of the fractions of deprotonated EZA and protonated CA, respectively. (c) The observed and intrinsic Gibbs free energies of TFMSA binding to hCA XIII as a function of pH. Solid bent line is the fit according to Eq. (1) recalculated to Gibbs free energies. Straight line shows the intrinsic Gibbs free energy of TFMSA binding, independent of pH.

Dissection of contribution from inhibitor protonation requires a set of experiments with pure inhibitor in the absence of protein to determine the pKa and the enthalpy of inhibitor protonation. The deprotonation pKa could be determined by pH titration as previously explained [30]. The enthalpy of protonation is determined by titrating alkaline solution of the inhibitor with acid. Figure 6 shows a typical ITC curve and the enthalpy of protonation dependence on temperature for both EZA and TFMSA. Pure crystalline inhibitor is in protonated electrostatically neutral form. Therefore, 1.5 equivalent of NaOH is added to the inhibitor solution to make it deprotonated. The ITC curve (Figure 6a, b) exhibits two transitions – the first transition is due to the titration of 0.5 equivalent of surplus hydroxide while the second transition is due to inhibitor protonation. This enthalpy of the second transition is the enthalpy of interest. The temperature dependence of EZA and TFMSA protonation is shown in Figure 6c. Table 2 shows the thermodynamic parameters of protonation of both inhibitors.

thumbnailFigure 6. Determination of the EZA protonation enthalpy. (a) ITC raw titration curve of 0.25 mM EZA in water with 1.5 equivalent NaOH titrated with 2.5 mM HNO3 at 37°C. (b) Integrated ITC curves of EZA () and TFMSA (▲) titration with HNO3 at 37°C. The first transition represents the neutralization of 0.5 equivalent of NaOH with HNO3 while the second transition (between the molar ratio of 0.5 and 1.5) represents the titration of the compound. In these experiments the enthalpies of protonation were −21.3 kJ/mol for TFMSA and −27.6 kcal/mol for EZA. (c) The enthalpies of EZA () and TFMSA (▲) protonation as a function of temperature. The slopes yield the heat capacities of protonation equal to 88.7 J/(mol × K) for EZA and 124.7 J/(mol × K) for TFMSA.

Table 2. Thermodynamic parameters of the protonation of inhibitor sulfonamide group and hCA XIII zinc-bound hydroxide anion at 25°C

The pKa of ethoxzolamide protonation has been previously determined to be equal to 8.0 [31]. The enthalpy of ethoxzolamide deprotonation has been determined by ITC (Figure 6). The pKas and enthalpies of other widely used inhibitors, including acetazolamide (AZM), methazolamide, dichlorophenamide, trifluoromethanesulfonamide (TFMSA), and p-amino benzene sulfonamide, have been previously determined [30].

Intrinsic binding parameters

Intrinsic Gibbs free energy of EZA and TFMSA binding was determined by the global fit of data as in Figure 5 while the intrinsic enthalpy was obtained by globally fitting data in Figure 3. The intrinsic entropy of binding was obtained by subtracting Gibbs free energy from the enthalpy using standard thermodynamic equations. The intrinsic heat capacity of binding could be determined by carrying out all or part of ITC and TSA experiments at several temperatures. This was not fully accomplished in this set of experiments. However, because at pH 7.0 most of the hCA XIII zinc-bound water hydroxide is protonated, only a very small error would be obtained by disregarding the heat capacity of protonation of the zinc-bound hydroxide. The resultant intrinsic thermodynamic parameters of binding are listed in Table 3.

Table 3. The intrinsic thermodynamic binding parameters of deprotonated inhibitor anion to the protonated form of zinc-bound water molecule of hCA XIII at 25°C

The intrinsic binding constants were greater than the observed ones and were equal to 1.8 × 109 M-1 for EZA and 4.6 × 107 M-1 for TFMSA. The enthalpies of binding were respectively equal to −42.1 and −33.5 kJ/mol. The entropies of binding of both inhibitors had a relatively small but positive contribution equal to 10.7 and 10.2 kJ/mol respectively. The binding process is therefore primarily enthalpy-driven and slightly supported by favorable entropic contribution for both inhibitors.

The intrinsic heat capacity of EZA binding to hCA XIII was determined to be −750 J/(mol × K). This is significantly different from the observed heat capacity of EZA binding at pH 7.0 (Figure 2e) equal to −980 J/(mol × K). The difference between observed and intrinsic parameters is important even for the heat capacity and again illustrates the importance of determining the intrinsic parameters.

The intrinsic binding thermodynamic parameters of acetazolamide (AZM) to hCA XIII were estimated by applying the same additivity principles from the data at single pH 7.0 and single temperature. The pKa and enthalpy of the protein protonation were determined using EZA and TFMSA as ligands while the pKa and enthalpy of AZM protonation have been previously determined [30] (Table 2). Using these parameters in the additivity scheme yields the intrinsic parameters listed in Table 3. Interestingly, binding of AZM had the largest intrinsic exothermic enthalpy of the three inhibitors, but the Gibbs free energy was least exergonic. Therefore, the entropic contribution to the AZM binding reaction was negative.

Discussion

Detailed analysis with numerous ITC experiments at various pHs is needed to be performed just once with one selected strongly-binding inhibitor. In this study, such inhibitor was selected to be EZA. To confirm the intrinsic binding parameters, additional ligand, TFMSA was studied and showed that the same set of parameters is suitable to describe the binding of this ligand. The binding parameters of the third ligand, AZM, were determined from experiments at single pH 7.0. All other sulfonamide inhibitors to be designed and studied in the future will only need one set of ITC and TSA experiments at a single pH. In addition, the pKa and enthalpy of the sulfonamide protonation will have to be determined.

Most studies, describing novel sulfonamide inhibitors listed only the observed Kds and/or Kis. [8,14,28,32,33]. The observed parameters provide estimates of the potency of each inhibitor at particular experimental conditions. However, these parameters are not suitable for the SAR-type studies and may mislead the rational design of novel inhibitors. Quite often, the observed Kds may differ by 3 orders of magnitude, but the intrinsic Kds may be comparable. In such cases, one may erroneously assign the greater observed potency to the protein-ligand interaction surface, while the true reason is simply the difference in pKa of the ligand sulfonamide group.

A good example is the comparison of intrinsic and observed Kds of TFMSA and methanesulfonamide interaction with CA II. Their observed Kis and Kds at pH 7.0 are 2 × 10-9 M and 3 × 10-4 M, respectively. However, their intrinsic Kds are 1.9 × 10-9 M and 2.4 × 10-8 M, respectively. The difference between the observed and intrinsic parameters for TFMSA is quite negligible, but for methanesulfonamide is very large. Based on the observed parameters, one could erroneously conclude that the three fluorine atoms increase the inhibition constant by 150 000 times as compared to methane group. In reality, the fluorines only increase the intrinsic binding constant by only 12 times. In other words, the fluorines bind to CA II 12 times stronger than methane hydrogens and this value could be used for the SAR of the fluorines. The remaining part is due to the difference in compound pKas and the available fractions of anionic sulfonamide at pH 7.0. Therefore, the studies presenting novel inhibitors should either determine the intrinsic parameters of binding or avoid pursuing observed Kd/Ki-based SAR analysis [26].

Acetazolamide is one the most widely used inhibitor to compare results between various laboratories. Here we determined the observed Kd for AZM binding to hCA XIII at pH 7.0 by both ITC and TSA to be 113 ± 30 nM. The intrinsic Kd is then calculated to be 53 nM. This insignificant difference from the observed value is due to the relatively low pKa of AZM (7.30). Large fraction of AZM is in its active deprotonated form at pH 7.0. Supuran with coworkers lists the observed Ki of AZM binding to murine isozyme mCA XIII to be equal to 17 nM [13]. The difference may be due to the different organism - human and murine isozymes.

The example of AZM shows the case when intrinsic and observed binding parameters are quite similar and may appear to be within the error region. However, most CA inhibitors have their pKas in the range between 9 and 11. Such inhibitors would exhibit observed Kds at pH 7.0 that are 2 to 4 orders of magnitude different from the intrinsic Kds. Large hidden differences in their pKas could mislead the researchers to the erroneous reasons of the binding potency.

Conclusions

Ethoxzolamide binding to hCA XIII is linked to several protonation reactions that must be dissected in order to estimate the intrinsic Gibbs free energies, enthalpies, entropies, and heat capacities of binding. Observed binding parameters for at least one inhibitor should be determined in a wide pH range. Combination of both methods, ITC and TSA is desired when the observed Kds are near the limit of precise determination range of ITC. Intrinsic enthalpy of EZA binding is −42.1 kJ/mol, intrinsic Gibbs free energy −52.8, and the intrinsic entropy of binding is 35.9 J/(mol × K).

Methods

Materials

Ethoxzolamide (EZA) was purchased from Aldrich (Milwaukee, WI, USA), trifluoromethanesulfonamide (TFMSA) was purchased from Alfa Aesar. Human carbonic anhydrase XIII (hCA XIII) was expressed in Escherichia coli and purified as previously described [16]. Ethics approval: no human or animal samples or live animals were used in this research.

Isothermal titration calorimetry (ITC)

ITC experiments were performed using a VP-ITC instrument (Microcal, Inc., GE-Healthcare, Northampton, MA, USA). For ligand-protein binding experiments 5–20 μM (usually 10) protein solution in the cell and 50–200 μM (usually 100) ligand solution in the syringe were used. A typical experiment consisted of 25–30 injections (usually 25), 10 μL each, with 3–5 min (usually 4 min) intervals between injections. Experiments were performed at different temperatures in various buffers containing 50 mM NaCl and 0.5-2% DMSO concentration. All experiments were repeated at least twice. For the determination of ligand binding at various pHs, the solution pH was carefully checked in the presence of protein and ligand, in the cell and syringe, respectively. The buffering capacity is low at some pH values (for example, phosphate at pH 9.5 or TRIS at pH 5.5) and had to be carefully adjusted with acidic or basic components of the buffer.

Thermal Shift Assay (TSA)

TSA experiments were performed with Corbett Rotor-Gene 6000 (QIAGEN Rotor-Gene Q) instrument using the blue channel (excitation 365 ± 20, detection 460 ± 15 nm), applying the heating rate of 1°C/min. Usually, the samples contained 20 μL of 10 μM protein, 0–200 μM ligand, 50 μM 8-anilino-1-naphthalene sulfonate, 50 mM sodium phosphate at pH 7.0, 50 mM NaCl , and 2% DMSO. The pH dependence of the Kb was determined in the universal buffer containing 50 mM sodium phosphate, 50 mM sodium acetate, and 25 mM sodium borate. The fluorescence and ligand dosing curves were fit as previously described [27,28]. TSA experiments were repeated at least three times. The standard deviation was below 10%.

Differential Scanning Calorimetry (DSC)

DSC experiments were performed using MC-2 Scanning Calorimeter (Microcal, Inc., GE-Healthcare, Northampton, MA, USA). The samples of 1.2 mL (cuvette volume) contained 50 μM hCA XIII, 0 or 500 μM ethoxzolamide, 50 mM sodium phosphate (pH 7.0), and 50 mM NaCl. A scan rate of 1°C /min was used.

Abbreviations

AZM: Acetazolamide; CA: Carbonic anhydrase; ΔbCp_intr: Intrinsic heat capacity of ligand binding; Δb_protonG: Gibbs free energy of protonation; ΔbGintr: Intrinsic Gibbs free energy of binding; ΔbGobs: Gibbs free energy of binding; Δb_protonCp: Heat capacity of protonation; Δb_protonH: Enthalpy of protonation; ΔbHintr: Intrinsic enthalpy of binding; ΔbHobs: Observed enthalpy of binding; Δb_protonS: Entropy of protonation; ΔbSintr: Intrinsic entropy of binding; ΔbSobs: Observed entropy of binding; EZA: Ethoxzolamide; ITC: Isothermal titration calorimetry; Ka: Acidity constant (proton dissociation, deprotonation constant); Kb_intr: Intrinsic binding constant; Kb_obs: Observed binding constant; nH: Number of protons in a protonation/deprotonation reaction; Pi: Phosphate buffer; TFMSA: Trifluoromethanesulfonamide; TRIS: TRIS (Trizma) buffer; TSA: Thermal shift assay.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

LB participated in the design of this study, performed all experiments, and participated in the analysis of the data and manuscript writing. DM designed the study and wrote the manuscript. All authors read and approved the final manuscript.

Acknowledgements

This research was funded by a grant (No. LIG-16/2010) from the Research Council of Lithuania.

References

  1. Stadie W, O’Brien H: The catalysis of the hydration of carbon dioxide and dehydration of carbonic acid by an enzyme isolated from red blood cells.

    J Biol Chem 1933, 103:521-529. OpenURL

  2. Khalifah RG: Carbon dioxide hydration activity of carbonic anhydrase: paradoxical consequences of the unusually rapid catalysis.

    Proc Natl Acad Sci U S A 1973, 70(7):1986-1989. OpenURL

  3. Silverman DN, Lindskog S: The catalytic mechanism of carbonic anhydrase: implications of a rate-limiting proteolysis of water.

    Acc Chem Res 1988, 21:30-36. OpenURL

  4. Sly WS, Hu PY: Human carbonic anhydrases and carbonic anhydrase deficiencies.

    Annu Rev Biochem 1995, 64:375-401. OpenURL

  5. Liljas A, Laurberg M: A wheel invented three times. The molecular structures of the three carbonic anhydrases.

    EMBO Rep 2000, 1:16-17. OpenURL

  6. Supuran CT: Carbonic anhydrase inhibition with natural products: novel chemotypes and inhibition mechanisms.

    Mol Divers 2011, 15(2):305-316. OpenURL

  7. Supuran CT, Fiore AD, Alterio V, Monti SM, Simone GD: Recent Advances in Structural Studies of the Carbonic Anhydrase Family: The Crystal Structure of Human CA IX and CA XIII.

    Curr Pharm Des 2010, 16:3346-3254. OpenURL

  8. Supuran CT: Carbonic anhydrases: novel therapeutic applications for inhibitors and activators.

    Nat Rev Drug Discov 2008, 7(2):168-181. OpenURL

  9. Supuran CT, Scozzafava A: Carbonic anhydrases as targets for medicinal chemistry.

    Bioorg Med Chem 2007, 15(13):4336-4350. OpenURL

  10. Ladbury JE, Klebe G, Freire E: Adding calorimetric data to decision making in lead discovery: a hot tip.

    Nat Rev Drug Discov 2010, 9:23-27. OpenURL

  11. Lehtonen J, Shen B, Vihinen M, Casini A, Scozzafava A, Supuran CT, Parkkila AK, Saarnio J, Kivela AJ, Waheed A, Sly WS, Parkkila S: Characterization of CA XIII, a novel member of the carbonic anhydrase isozyme family.

    J Biol Chem 2004, 279(4):2719-2727. OpenURL

  12. Fiore AD, Monti SM, Hilvo M, Parkkila S, Romano V, Scaloni A, Pedone C, Scozzafava A, Supuran CT, Simone GD: Crystal structure of human carbonic anhydrase XIII and its complex with the inhibitor acetazolamide.

    Proteins 2009, 74:164-175. OpenURL

  13. Lehtonen JM, Parkkila S, Vullo D, Casini A, Scozzafava A, Supuran CT: Carbonic anhydrase inhibitors. Inhibition of cytosolic isozyme XIII with aromatic and heterocyclic sulfonamides: a novel target for the drug design.

    Bioorg Med Chem Lett 2004, 14(14):3757-3762. OpenURL

  14. Vernier W, Chong W, Rewolinski D, Greasley S, Pauly T, Shaw M, Dinh D, Ferre RA, Nukui S, Ornelas M, Reyner E: Thioether benzenesulfonamide inhibitors of carbonic anhydrases II and IV: structure-based drug design, synthesis, and biological evaluation.

    Bioorg Med Chem 2010, 18(9):3307-3319. OpenURL

  15. Capkauskaite E, Baranauskiene L, Golovenko D, Manakova E, Grazulis S, Tumkevicius S, Matulis D: Indapamide-like benzenesulfonamides as inhibitors of carbonic anhydrases I, II, VII, and XIII.

    Bioorg Med Chem 2010, 18(21):7357-7364. OpenURL

  16. Sudzius J, Baranauskiene L, Golovenko D, Matuliene J, Michailoviene V, Torresan J, Jachno J, Sukackaite R, Manakova E, Grazulis S, Tumkevicius S, Matulis D: 4-[N-(substituted 4-pyrimidinyl)amino]benzenesulfonamides as inhibitors of carbonic anhydrase isozymes I, II, VII, and XIII.

    Bioorg Med Chem 2010, 18(21):7413-7421. OpenURL

  17. Khalifah RG, Zhang F, Parr JS, Rowe ES: Thermodynamics of binding of the CO2-competitive inhibitor imidazole and related compounds to human carbonic anhydrase I: an isothermal titration calorimetry approach to studying weak binding by displacement with strong inhibitors.

    Biochemistry 1993, 32(12):3058-3066. OpenURL

  18. Baker BM, Murphy KP: Evaluation of linked protonation effects in protein binding reactions using isothermal titration calorimetry.

    Biophys J 1996, 71(4):2049-2055. OpenURL

  19. Krishnamurthy VM, Kaufman GK, Urbach AR, Gitlin I, Gudiksen KL, Weibel DB, Whitesides GM: Carbonic anhydrase as a model for biophysical and physical-organic studies of proteins and protein-ligand binding.

    Chem Rev 2008, 108(3):946-1051. OpenURL

  20. Velazquez Campoy A, Freire E: ITC in the post-genomic era…? Priceless.

    Biophys Chem 2005, 115(2–3):115-124. OpenURL

  21. Velazquez-Campoy A, Freire E: Isothermal titration calorimetry to determine association constants for high-affinity ligands.

    Nat Protoc 2006, 1:186-191. OpenURL

  22. Matulis D, Kranz JK, Salemme FR, Todd MJ: Thermodynamic stability of carbonic anhydrase: measurements of binding affinity and stoichiometry using ThermoFluor.

    Biochemistry 2005, 44(13):5258-5266. OpenURL

  23. Mezzasalma TM, Kranz JK, Chan W, Struble GT, Schalk-Hihi C, Deckman IC, Springer BA, Todd MJ: Enhancing recombinant protein quality and yield by protein stability profiling.

    J Biomol Screen 2007, 12(3):418-428. OpenURL

  24. Zubriene A, Matuliene J, Baranauskiene L, Jachno J, Torresan J, Michailoviene V, Cimmperman P, Matulis D: Measurement of Nanomolar Dissociation Constants by Titration Calorimetry and Thermal Shift Assay - Radicicol Binding to Hsp90 and Ethoxzolamide Binding to CAII.

    Int J Mol Sci 2009, 10(6):2662-2680. OpenURL

  25. Kranz JK, Schalk-Hihi C: Protein thermal shifts to identify low molecular weight fragments.

    Methods Enzymol 2011, 493:277-298. OpenURL

  26. Zubriene A, Kazlauskas E, Baranauskiene L, Petrauskas V, Matulis D: Isothermal titration calorimetry and thermal shift assay in drug design.

    European pharmaceutical review 2011, 16:56-59. OpenURL

  27. Cimmperman P, Matulis D: Protein thermal denaturation measurements via a fluorescent dye. RSC Publishing; 2011:247-274.

    chap. 8

    OpenURL

  28. Baranauskiene L, Hilvo M, Matuliene J, Golovenko D, Manakova E, Dudutiene V, Michailoviene V, Torresan J, Jachno J, Parkkila S, Maresca A, Supuran CT, Grazulis S, Matulis D: Inhibition and binding studies of carbonic anhydrase isozymes I, II and IX with benzimidazo[1,2-c][1,2,3]thiadiazole-7-sulphonamides.

    J Enzyme Inhib Med Chem 2010, 25(6):863-870. OpenURL

  29. Cimmperman P, Baranauskiene L, Jachimoviciute S, Jachno J, Torresan J, Michailoviene V, Matuliene J, Sereikaite J, Bumelis V, Matulis D: A quantitative model of thermal stabilization and destabilization of proteins by ligands.

    Biophys J 2008, 95(7):3222-3231. OpenURL

  30. Matulis D, Todd MJ: Thermodynamics - Structure Correlations of Sulfonamide Inhibitor Binding to Carbonic Anhydrase. Wiley; 2004:107-132. OpenURL

  31. Remko M, von der Lieth CW: Theoretical study of gas-phase acidity, pKa, lipophilicity, and solubility of some biologically active sulfonamides.

    Bioorg Med Chem 2004, 12(20):5395-5403. OpenURL

  32. Innocenti A, Scozzafava A, Supuran CT: Carbonic anhydrase inhibitors. Inhibition of transmembrane isoforms IX, XII, and XIV with less investigated anions including trithiocarbonate and dithiocarbamate.

    Bioorg Med Chem Lett 2010, 20(5):1548-1550. OpenURL

  33. Thiry A, Ledecq M, Cecchi A, Frederick R, Dogne JM, Supuran CT, Wouters J, Masereel B: Ligand-based and structure-based virtual screening to identify carbonic anhydrase IX inhibitors.

    Bioorg Med Chem 2009, 17(2):553-557. OpenURL