Ultrahigh resolution crystal structures of human carbonic anhydrases I and II complexed with "two-prong" inhibitors reveal the molecular basis of high affinity

The atomic-resolution crystal structures of human carbonic anhydrases I and II complexed with "two-prong" inhibitors are reported. Each inhibitor contains a benzenesulfonamide prong and a cupric iminodiacetate (IDA-Cu(2+)) prong separated by linkers of different lengths and compositions. The ionized NH(-) group of each benzenesulfonamide coordinates to the active site Zn(2+) ion; the IDA-Cu(2+) prong of the tightest-binding inhibitor, BR30, binds to H64 of CAII and H200 of CAI. This work provides the first evidence verifying the structural basis of nanomolar affinity measured for two-prong inhibitors targeting the carbonic anhydrases.     

Two-prong inhibitors for human carbonic anhydrase II

The enzyme inhibitors are usually designed by taking into consideration the overall dimensions of the enzyme's active site pockets. This conventional approach often fails to produce desirable affinities of inhibitors for their cognate enzymes. To circumvent such constraints, we contemplated enhancing the binding affinities of inhibitors by attaching tether groups, which would interact with the surface exposed amino acid residues. This strategy has been tested for the inhibition of human carbonic anhydrase II. Benzenesulfonamide serves as a weak inhibitor for the enzyme, but when it is conjugated to iminodiacetate-Cu2+ (which interacts with the surface-exposed His residues) via a spacer group, its binding affinity is enhanced by about 2 orders of magnitude. This "two-prong" approach is expected to serve as a general strategy for converting weak inhibitors of enzymes into tight-binding inhibitors.     

Models of F.H contacts relevant to the binding of fluoroaromatic inhibitors to carbonic anhydrase II

[formula: see text] Complexes formed between fluorobenzene and N-methylformamide or benzene have been used as models of the interaction of fluoroaromatic drugs with carbonic anhydrase II. These structures have been investigated via ab initio and density functional methods, including HF, B3LYP, and MP2 procedures. The results of the calculations are consistent with the hypothesis, suggested originally by experimental X-ray crystal structures of the drug-receptor complexes, that favorable fluorine-hydrogen interactions affect binding affinity.     

Fluoroaromatic-fluoroaromatic interactions between inhibitors bound in the crystal lattice of human carbonic anhydrase II

Intermolecular interactions of eleven different fluoroaromatic inhibitors are probed within the scaffolding of the crystal lattice of Phe-131-->Val carbonic anhydrase II. The degree and pattern of fluorine substitution on the inhibitor benzyl ring modulate its size, shape, and electronic character. In turn, these properties affect the geometry of intermolecular interactions between the fluoroaromatic rings of two different inhibitor molecules bound in the crystal lattice, as determined by X-ray crystallography. Depending on the degree and pattern of fluorine substitution, we observe a face-to-face (aromatic-aromatic) interaction, an atom-to-face (carbonyl-aromatic) interaction, or no interaction at all. These interaction geometries are analyzed with regard to van der Waals, electrostatic, and possible charge-transfer effects. For the aromatic-aromatic interactions investigated in this study, with aromatic ring quadrupoles specifically "tuned" by the degree and pattern of fluorination, the structural results suggest that London forces and charge-transfer complexation dominate over weakly polar electrostatic interactions in the association of aromatic ring pairs.     

The binding of human carbonic anhydrase II by functionalized folded polypeptide receptors

Several receptors for human carbonic anhydrase II (HCAII) have been prepared by covalently attaching benzenesulfonamide carboxylates via aliphatic aminocarboxylic acid spacers of variable length to the side chain of a lysine residue in a designed 42 residue helix-loop-helix motif. The sulfonamide group binds to the active site zinc ion of human carbonic anhydrase II located in a 15 A deep cleft. The dissociation constants of the receptor-HCAII complexes were found to be in the range from low micromolar to better than 20 nM, with the lowest affinities found for spacers with less than five methylene groups and the highest affinity found for the spacer with seven methylene groups. The results suggest that the binding is a cooperative event in which both the sulfonamide residue and the helix-loop-helix motif contribute to the overall affinity.     

Linear free energy relationships implicate three modes of binding for fluoroaromatic inhibitors to a mutant of carbonic anhydrase II

[figure: see text] Linear free energy relationships between binding affinity and hydrophobicity for a library of fluoroaromatic inhibitors of F131V carbonic anhydrase II (CA) implicate three modes of interaction. X-ray crystal structures suggest that F131 interacts with fluoroaromatic inhibitors, while P202, on the opposite side of the active site cleft, serves as the site of the hydrophobic contact in the case of the F131V mutant. 2-Fluorinated compounds bind more tightly, perhaps due to the field effect of the nearby fluorine on the acidity of the amide proton.     

Characterization of the first beta-class carbonic anhydrase from an arthropod (Drosophila melanogaster) and phylogenetic analysis of beta-class carbonic anhydrases in invertebrates

Background

The β-carbonic anhydrase (CA, EC 4.2.1.1) enzymes have been reported in a variety of organisms, but their existence in animals has been unclear. The purpose of the present study was to perform extensive sequence analysis to show that the β-CAs are present in invertebrates and to clone and characterize a member of this enzyme family from a representative model organism of the animal kingdom, e.g., Drosophila melanogaster.

Results

The novel β-CA gene, here named DmBCA, was identified from FlyBase, and its orthologs were searched and reconstructed from sequence databases, confirming the presence of β-CA sequences in 55 metazoan species. The corresponding recombinant enzyme was produced in Sf9 insect cells, purified, kinetically characterized, and its inhibition was investigated with a series of simple, inorganic anions. Holoenzyme molecular mass was defined by dynamic light scattering analysis and gel filtration, and the results suggested that the holoenzyme is a dimer. Double immunostaining confirmed predictions based on sequence analysis and localized DmBCA protein to mitochondria. The enzyme showed high CO₂ hydratase activity, with a k_cat of 9.5 × 10⁵ s^-1 and a k_cat/K_M of 1.1 × 10⁸ M^{- 1}s^{- 1}. DmBCA was appreciably inhibited by the clinically-used sulfonamide acetazolamide, with an inhibition constant of 49 nM. It was moderately inhibited by halides, pseudohalides, hydrogen sulfide, bisulfite and sulfate (K_I values of 0.67 - 1.36 mM) and more potently by sulfamide (K_I of 0.15 mM). Bicarbonate, nitrate, nitrite and phenylarsonic/boronic acids were much weaker inhibitors (K_Is of 26.9 - 43.7 mM).

Conclusions

The Drosophila β-CA represents a highly active mitochondrial enzyme that is a potential model enzyme for anti-parasitic drug development.     

Residual ligand entropy in the binding of p-substituted benzenesulfonamide ligands to bovine carbonic anhydrase II

In studies on the thermodynamics of ligand-protein interactions, it is often assumed that the configurational and conformational entropy of the ligand is zero in the bound state (i.e., the ligand is rigidly fixed in the binding pocket). However, there is little direct experimental evidence for this assumption, and in the case of binding of p-substituted benzenesulfonamide inhibitors to bovine carbonic anhydrase II (BCA II), the observed thermodynamic binding signature derived from isothermal titration calorimetry experiments leads indirectly to the conclusion that a considerable degree of residual entropy remains in the bound ligand. Specifically, the entropy of binding increases with glycine chain length n, and strong evidence exists that this thermodynamic signature is not driven by solvent reorganization. By use of heteronuclear (15)N NMR relaxation measurements in a series (n = 1-6) of (15)N-glycine-enriched ligands, we find that the observed thermodynamic binding signature cannot be explained by residual ligand dynamics in the bound state, but rather results from the indirect influence of ligand chain length on protein dynamics.     

Determination of acidity constants,ionic mobilities,and hydrodynamic radii of carborane-based inhibitors of carbonic anhydrases by capillary electrophoresis

Capillary electrophoresis (CE) has been applied for determination of the thermodynamic acidity constants (pK_a) of the sulfamidoalkyl and sulfonamidoalkyl groups, the actual and limiting ionic mobilities and hydrodynamic radii of important compounds, eight carborane-based inhibitors of carbonic anhydrases, which are potential new anticancer drugs. Two types of carboranes were investigated, (i) icosahedral cobalt bis(dicarbollide)(1-) ion with sulfamidoalkyl moieties, and (ii) 7,8-nido-dicarbaundecaborate with sulfonamidoalkyl side chains. First, the mixed acidity constants, pK_a^mix, of the sulfamidoalkyl and sulfonamidoalkyl groups of the above carboranes and their actual ionic mobilities were determined by nonlinear regression analysis of the pH dependences of their effective electrophoretic mobility measured by capillary electrophoresis in the pH range 8.00−12.25, at constant ionic strength (25 mM), and constant temperature (25°C). Second, the pK_a^mix were recalculated to the thermodynamic pK_as using the Debye–Hückel theory. The sulfamidoalkyl and sulfonamidoalkyl groups were found to be very weakly acidic with the pK_as in the range 10.78−11.45 depending on the type of carborane cluster and on the position and length of the alkyl chain on the carborane scaffold. These pK_as were in a good agreement with the pK_as (10.67−11.27) obtained by new program AnglerFish (freeware at https://echmet.natur.cuni.cz ), which provides thermodynamic pK_as and limiting ionic mobilities directly from the raw CE data. The absolute values of the limiting ionic mobilities of univalent and divalent carborane anions were in the range 18.3−27.8 TU (Tiselius unit, 1 × 10⁻⁹ m²/Vs), and 36.4−45.9 TU, respectively. The Stokes hydrodynamic radii of univalent and divalent carborane anions varied in the range 0.34−0.52 and 0.42−0.52 nm, respectively.     

Bicarbonate as a proton donor in catalysis by Zn(II)- and Co(II)-containing carbonic anhydrases.

Catalysis of (18)O exchange between CO(2) and water catalyzed by a Co(II)-substituted mutant of human carbonic anhydrase II is analyzed to show the rate of release of H(2)(18)O from the active site. This rate, measured by mass spectrometry, is dependent on proton transfer to the metal-bound (18)O-labeled hydroxide, and was observed in a site-specific mutant of carbonic anhydrase II in which a prominent proton shuttle residue His64 was replaced by alanine, which does not support proton transport. Upon increasing the concentration of bicarbonate, the rate of release of H(2)(18)O increased in a saturable manner to a maximum of 4 x 10(5) s(-)(1), consistent with proton transfer from bicarbonate to the Co(II)-bound hydroxide. The same mutant of carbonic anhydrase containing Zn(II) had the rate of release of H(2)(18)O smaller by 10-fold, but rate of interconversion of CO(2) and HCO(3)(-) about the same as the Co(II)-containing enzyme. These data as well as solvent hydrogen isotope effects suggest that the bicarbonate transferring the proton is bound to the cobalt in the enzyme. The enhancement of (18)O exchange caused by increasing bicarbonate concentration during catalysis by the Zn(II)-containing carbonic anhydrase from the archaeon Methanosarcina thermophila suggests that a very similar mechanism for proton donation by bicarbonate occurs with this wild-type enzyme.     

Carbonic anhydrase inhibitors: X-ray crystallographic studies for the binding of N-substituted benzenesulfonamides to human isoform II

N-substituted benzenesulfonamides, incorporating the N-amino-, N-hydroxy- and N-methoxy-moieties at the sulfonamide zinc binding group, have been investigated as CAIs by means of inhibition and structural studies, unraveling interesting aspects related to their inhibition mechanism.     

Identification of proton-transfer pathways in human carbonic anhydrase II

We investigate the probable proton-transfer pathways from the surface of human carbonic anhydrase II into the active site cavity through His-64 that has been widely implicated as a key residue along the proton-transfer path. A recursive analysis of hydrogen-bonded clusters in the static crystallographic structure shows that there is no complete path through His-64 in either of its experimentally detected conformations. Side chain conformational fluctuation of His-64 from its outward conformation toward the active site is found to provide a crucial dynamic connectivity needed to complete the path coupled to local reorganization of the protein structure and hydration. The energy and free energy barriers along the detected pathway have been estimated to derive the mechanism of His-64 rotation toward the active site. We also investigate a dynamical connectivity map that highlights networks of disordered water molecules that may promote a direct (and probably transient) access of the solvent to the active site. Our studies reveal how such solvent access channels may be related to the putative proton shuttle mediated by His-64. The paths thus identified can be potentially used as reaction coordinates for further studies on the molecular mechanism of enzyme action.     

Screening and docking studies of natural phenolic inhibitors of carbonic anhydrase II

Carbonic anhydrase II (CA II) is an important enzyme complex with Zn²⁺, which is involved in many physiological and pathological processes, such as calcification, glaucoma and tumorigenicity. In order to search for novel inhibitors of CA II, inhibition assay of carbonic anhydrase II was performed, by which seven natural phenolic compounds, including four phenolics (grifolin, 4-O-methyl-grifolic acid, grifolic acid, and isovanillic acid) and three flavones (eriodictyol, quercetin and puerin A), showed inhibitory activities against CA II with IC₅₀s in the range of 6.37–71.73 μmol/L. Grifolic acid is the most active one with IC₅₀ of 6.37 _μmol/L. These seven phenolic compounds were proved to be novel natural carbonic anhydrase II inhibitors, which were obtained in flexible docking study with GOLD 3.0 software. Results indicated that the aliphatic chain and polar groups of hydroxyl and carboxyl are important to their inhibitory activities, providing a new insight into study on CA II potent inhibitors. Authors with the equal contribution Supported by the National Natural Science Foundation of China (Grant No. 30725048) and the Foundation of Chinese Academy of Sciences (West Light Program).     

Calculation of binding free energies of inhibitors to plasmepsin II

An understanding at the atomic level of the driving forces of inhibitor binding to the protein plasmepsin (PM) II would be of interest to the development of drugs against malaria. To this end, three state of the art computational techniques to compute relative free energies-thermodynamic integration (TI), Hamiltonian replica-exchange (H-RE) TI, and comparison of bound versus unbound ligand energy and entropy-were applied to a protein-ligand system of PM II and several exo-3-amino-7-azabicyclo[2.2.1]heptanes and the resulting relative free energies were compared with values derived from experimental IC(50) values. For this large and flexible protein-ligand system, the simulations could not properly sample the relevant parts of the conformational space of the bound ligand, resulting in failure to reproduce the experimental data. Yet, the use of Hamiltonian replica exchange in conjunction with thermodynamic integration resulted in enhanced convergence and computational efficiency compared to standard thermodynamic integration calculations. The more approximate method of calculating only energetic and entropic contributions of the ligand in its bound and unbound states from conventional molecular dynamics (MD) simulations reproduced the major trends in the experimental binding free energies, which could be rationalized in terms of energetic and entropic characteristics of the different structural and physico-chemical properties of the protein and ligands.     

Rapid ion-exchange chromatography for preparative separation of proteins. Application to porcine and bovine carbonic anhydrases

A method of rapid ion-exchange chromatography of DEAE-cellulose for preparative purposes is described. Basically, the flow-rate is increased by applying an air pressure on the column. By this technique it is possible to purify gram quantities of protein in 2-4 h with acceptable resolution. In preparations of bovine and porcine carbonic anhydrases the elution times were reduced by a factor of about ten compared to those of conventional methods. The enzymes purified in this way showed a high degree of homogeneity. The method should be generally applicable in protein purification, and especially advantageous in purification of unstable proteins where time-consuming separations often give rise to low yields of active material.     

Using covalent dimers of human carbonic anhydrase II to model bivalency in immunoglobulins

This paper describes the development of a new bivalent system comprising synthetic dimers of carbonic anhydrase linked chemically through thiol groups of cysteine residues introduced by site-directed mutagenesis. These compounds serve as models with which to study the interaction of bivalent proteins with ligands presented at the surface of mixed self-assembled monolayers (SAMs). Monovalent carbonic anhydrase (CA) binds to benzenesulfonamide ligands presented on the surface of the SAM with K(d)(surf) = 89 nM. The synthetic bivalent proteins--inspired by the structure of immunoglobulins--bind bivalently to the sulfonamide-functionalized SAMs with low nanomolar avidities (K(d)(avidity,surf) = 1-3 nM); this difference represents a ~50-fold enhancement of bivalent over monovalent association. The paper describes dimers of CA having (i) different lengths of the covalent linker that joined the two proteins and (ii) different points of attachment of the linker to the protein (either near the active site (C133) or distal to the active site (C185)). Comparison of the thermodynamics of their interactions with SAMs presenting arylsulfonamide groups demonstrated that varying the length of the linker between the molecules of CA had virtually no effect on the rate of association, or on the avidity of these dimers with ligand-presenting surfaces. Varying the point of attachment of the linker between monomeric CA's also had almost no effect on the avidity of the dimers, although changing the point of attachment affected the rates of binding and unbinding. These observations indicate that the avidities of these bivalent proteins, and by inference the avidities of structurally similar bivalent proteins such as IgG, are unexpectedly insensitive to the structure of the linker connecting them.     

Molecular charge distribution and chemical binding in five-membered heterocycles. I

Electron density maps for the heterocycles thiophene, furan, and pyrrole are determined from ab initio 4-31G wavefunctions. The charge distributions in these molecules are analyzed in terms of the total molecular density and difference density maps and their profiles. The atomiclike core, especially the L core of sulfur, is found to play an important role, via its polarization and interaction, in determining the extent and direction of valence density transfer from the carbon to the heteroatom. The changes in the charge distributions that occur in the immediate vicinity of the heteroatoms and the relation of density quantities to binding and antibinding characteristics are discussed. The quantum topological features of the molecular charge distributions of the three heterocycles are analyzed and discussed and the different bonding situations, e.g., ring strain, ionic and covalent binding, etc., are compared in a model-independent way.