Structural and kinetic study of the extended active site for proton transfer in human carbonic anhydrase II

pH-Dependent Structure and Dynamics of the Catalytic Domains of Human Carbonic Anhydrase II and IX

Extensive computer simulation studies have been carried out to probe the pH-dependent structure and dynamics of the two most efficient isoenzymes II and IX of human carbonic anhydrase (HCA) that control the pH in the human body. The equilibrium structure and hydration of their catalytic domains are found to be largely unaffected by the variation of pH in the range studied, in close agreement with the known experimental results. In contrast, a significant effect of the change in pH is observed for the first time on the local electrostatic potential of the active site walls and the dynamics of active site water molecules. We also report for the first time the free energy and kinetics of coupled fluctuations of orientation and protonation states of the well-known His-mediated proton shuttle (His-64) in both isozymes at pH 7 and 8. The transitions between different tautomers of in or out conformations of His-64 side chain range between 109 and 106 s-1 depending on pH. Possible implications of these results on conformation-dependent pKa of His-64 side chain and its role in driving the catalysis toward hydration of CO2 or dehydration of HCO3- with varying pH are discussed.

Efficient proton shuttle makes SazCA an excellent CO2 hydration enzyme

The fastest member of the carbonic anhydrase family catalysing the reversible hydration of carbon dioxide to bicarbonate ions has been recently reported to be SazCA. While thermostable, this enzyme shows exceptional activity at 353 K for the reaction. This study explores the molecular basis for the exceptional activity of SazCA, in contrast to SspCA, probed using molecular dynamics simulations. Our simulations, carried out at different temperatures, indicate the presence of efficient proton shuttle between the active zinc centre and His64 residue in the two enzymes. The proton accepting His64 residue was identified to have in and out conformations with the in conformations being supportive to proton acceptance. Our simulations show a large population of in conformations in SazCA making the enzyme exhibit an exceptional activity. The RMSF and H-bonds analysis confirmed the role of His2 and His207 in supporting the attainment of in conformations in SazCA resulting in exceptional activity.Communicated by Ramaswamy H. Sarma.

Molecular modelling and dynamics of CA2 missense mutations causative to carbonic anhydrase 2 deficiency syndrome

Carbonic anhydrase 2 (CA2) enzyme deficiency caused by CA2 gene mutations is an inherited disorder characterized by symptoms like osteopetrosis, renal tubular acidosis, and cerebral calcification. This study has collected the CA2 deficiency causal missense mutations and assessed their pathogenicity using diverse computational programs. The 3D protein models for all missense mutations were built, and analyzed for structural divergence, protein stability, and molecular dynamics properties. We found M-CAP as the most sensitive prediction method to measure the deleterious potential of CA2 missense mutations. Free energy dynamics of tertiary structure models of CA2 mutants with DUET, mCSM, and SDM based consensus methods predicted only 50% of the variants as destabilizing. Superimposition of native and mutant CA2 models revealed the minor structural fluctuations at the amino acid residue level but not at the whole protein structure level. Near native molecular dynamic simulation analysis indicated that CA2 causative missense variants result in residue level fluctuation pattern in the protein structure. This study expands the understanding of genotype-protein phenotype correlations underlying CA2 variant pathogenicity and presents a potential avenue for modifying the CA2 deficiency by targeting biophysical structural features of CA2 protein. Communicated by Ramaswamy H. Sarma.

Exploration of the residues modulating the catalytic features of human carbonic anhydrase XIII by a site-specific mutagenesis approach

Carbonic anhydrases (CAs) are ubiquitous metallo-enzymes that catalyse the reversible hydration of carbon dioxide to bicarbonate and proton. In humans there are 15 isoforms among which only 12 are catalytically active. Since active human (h) CAs show different efficiency, the understanding of the molecular determinants affecting it is a matter of debate. Here we investigated, by a site-specific mutagenesis approach, residues modulating the catalytic features of one of the least investigated cytosolic isoform, i.e. hCA XIII. Results showed that residues assisting the formation of an ordered solvent network within the catalytic site as well as those forming a histidine cluster on the protein surface are important to guarantee an efficient proton transfer.

Thermodynamic, kinetic, and structural parameterization of human carbonic anhydrase interactions toward enhanced inhibitor design

The aim of rational drug design is to develop small molecules using a quantitative approach to optimize affinity. This should enhance the development of chemical compounds that would specifically, selectively, reversibly, and with high affinity interact with a target protein. It is not yet possible to develop such compounds using computational (i.e., in silico) approach and instead the lead molecules are discovered in high-throughput screening searches of large compound libraries. The main reason why in silico methods are not capable to deliver is our poor understanding of the compound structure-thermodynamics and structure-kinetics correlations. There is a need for databases of intrinsic binding parameters (e.g., the change upon binding in standard Gibbs energy (ΔGint), enthalpy (ΔHint), entropy (ΔSint), volume (ΔVintr), heat capacity (ΔCp,int), association rate (ka,int), and dissociation rate (kd,int)) between a series of closely related proteins and a chemically diverse, but pharmacophoric group-guided library of compounds together with the co-crystal structures that could help explain the structure-energetics correlations and rationally design novel compounds. Assembly of these data will facilitate attempts to provide correlations and train data for modeling of compound binding. Here, we report large datasets of the intrinsic thermodynamic and kinetic data including over 400 primary sulfonamide compound binding to a family of 12 catalytically active human carbonic anhydrases (CA). Thermodynamic parameters have been determined by the fluorescent thermal shift assay, isothermal titration calorimetry, and by the stopped-flow assay of the inhibition of enzymatic activity. Kinetic measurements were performed using surface plasmon resonance. Intrinsic thermodynamic and kinetic parameters of binding were determined by dissecting the binding-linked protonation reactions of the protein and sulfonamide. The compound structure-thermodynamics and kinetics correlations reported here helped to discover compounds that exhibited picomolar affinities, hour-long residence times, and million-fold selectivities over non-target CA isoforms. Drug-lead compounds are suggested for anticancer target CA IX and CA XII, antiglaucoma CA IV, antiobesity CA VA and CA VB, and other isoforms. Together with 85 X-ray crystallographic structures of 60 compounds bound to six CA isoforms, the database should be of help to continue developing the principles of rational target-based drug design.

The Crystal Structure of a hCA VII Variant Provides Insights into the Molecular Determinants Responsible for Its Catalytic Behavior

Although important progress has been achieved in understanding the catalytic mechanism of Carbonic Anhydrases, a detailed picture of all factors influencing the catalytic efficiency of the various human isoforms is still missing. In this paper we report a detailed structural study and theoretical pKa calculations on a hCA VII variant. The obtained data were compared with those already known for another thoroughly investigated cytosolic isoform, hCA II. Our structural studies show that in hCA VII the network of ordered water molecules, which connects the zinc bound solvent molecule to the proton shuttle His64, is altered compared to hCA II, causing a reduction of the catalytic efficiency. Theoretical calculations suggest that changes in solvent network are related to the difference in pKa of the proton shuttle in the two enzymes. The residue that plays a major role in determining the diverse pKa values of the proton shuttle is the one in position four, namely His for hCA II and Gly for hCA VII. This residue is located on the protein surface, outside of the active site cavity. These findings are in agreement with our previous studies that highlighted the importance of histidines on the protein surface of hCA II (among which His4) as crucial residues for the high catalytic efficiency of this isoform.

Crystal Structure of Carbonic Anhydrase II in Complex with an Activating Ligand: Implications in Neuronal Function

Carbonic anhydrase (CA) plays a key role in neuronal signaling, providing bicarbonate and proton ions for GABAergic and glutamatergic neuronal function. Activation of CA isoforms expressed in neurons have been shown to have implications in the prognosis of Alzheimer's disease and dementia, while inhibitors of CAs are clinically used in the treatment of epilepsy, emphasizing the importance of this family of enzymes in both disease and normal neuronal function. Previously, compounds have been reported to enhance activity of CAs in an aging rat model, but their mechanism of action was not known. We report the 1.6 Å resolution structure of an imidazole-based CA activator in complex with the ubiquitously-expressed human CA II. Based on the structure, a proposed mechanism of CA activation by the compound and its potential applications in the neurobiology of aging are discussed.

Carbonic anhydrase under pressure

Investigations of the rapid enzyme carbonic anhydrase have now been extended by crystallographic analysis at high CO2 pressures to examine the movements of water molecules in different steps of the catalysis. The rate of catalysis seems well explained by the assembled observations.

In Silico Studies of Small Molecule Interactions with Enzymes Reveal Aspects of Catalytic Function

Small molecules, such as solvent, substrate, and cofactor molecules, are key players in enzyme catalysis. Computational methods are powerful tools for exploring the dynamics and thermodynamics of these small molecules as they participate in or contribute to enzymatic processes. In-depth knowledge of how small molecule interactions and dynamics influence protein conformational dynamics and function is critical for progress in the field of enzyme catalysis. Although numerous computational studies have focused on enzyme-substrate complexes to gain insight into catalytic mechanisms, transition states and reaction rates, the dynamics of solvents, substrates, and cofactors are generally less well studied. Also, solvent dynamics within the biomolecular solvation layer play an important part in enzyme catalysis, but a full understanding of its role is hampered by its complexity. Moreover, passive substrate transport has been identified in certain enzymes, and the underlying principles of molecular recognition are an area of active investigation. Enzymes are highly dynamic entities that undergo different conformational changes, which range from side chain rearrangement of a residue to larger-scale conformational dynamics involving domains. These events may happen nearby or far away from the catalytic site, and may occur on different time scales, yet many are related to biological and catalytic function. Computational studies, primarily molecular dynamics (MD) simulations, provide atomistic-level insight and site-specific information on small molecule interactions, and their role in conformational pre-reorganization and dynamics in enzyme catalysis. The review is focused on MD simulation studies of small molecule interactions and dynamics to characterize and comprehend protein dynamics and function in catalyzed reactions. Experimental and theoretical methods available to complement and expand insight from MD simulations are discussed briefly.

Structure-Activity Relationships of Benzenesulfonamide-Based Inhibitors towards Carbonic Anhydrase Isoform Specificity

Carbonic anhydrases (CAs) are implicated in a wide range of diseases, including the upregulation of isoforms CA IX and XII in many aggressive cancers. However, effective inhibition of disease-implicated CAs should minimally affect the ubiquitously expressed isoforms, including CA I and II, to improve directed distribution of the inhibitors to the cancer-associated isoforms and reduce side effects. Four benzenesulfonamide-based inhibitors were synthesized by using the tail approach and displayed nanomolar affinities for several CA isoforms. The crystal structures of the inhibitors bound to a CA IX mimic and CA II are presented. Further in silico modeling was performed with the inhibitors docked into CA I and XII to identify residues that contributed to or hindered their binding interactions. These structural studies demonstrated that active-site residues lining the hydrophobic pocket, especially positions 92 and 131, dictate the positional binding and affinity of inhibitors, whereas the tail groups modulate CA isoform specificity. Geometry optimizations were performed on each ligand in the crystal structures and showed that the energetic penalties of the inhibitor conformations were negligible compared to the gains from active-site interactions. These studies further our understanding of obtaining isoform specificity when designing small molecule CA inhibitors.

The Structure of Carbonic Anhydrase IX Is Adapted for Low-pH Catalysis

Human carbonic anhydrase IX (hCA IX) expression in many cancers is associated with hypoxic tumors and poor patient outcome. Inhibitors of hCA IX have been used as anticancer agents with some entering Phase I clinical trials. hCA IX is transmembrane protein whose catalytic domain faces the extracellular tumor milieu, which is typically associated with an acidic microenvironment. Here, we show that the catalytic domain of hCA IX (hCA IX-c) exhibits the necessary biochemical and biophysical properties that allow for low pH stability and activity. Furthermore, the unfolding process of hCA IX-c appears to be reversible, and its catalytic efficiency is thought to be correlated directly with its stability between pH 3.0 and 8.0 but not above pH 8.0. To rationalize this, we determined the X-ray crystal structure of hCA IX-c to 1.6 Å resolution. Insights from this study suggest an understanding of hCA IX-c stability and activity in low-pH tumor microenvironments and may be applicable to determining pH-related effects on enzymes.

Inhibitory effects of nitrite on the reactions of bovine carbonic anhydrase II with CO2 and bicarbonate consistent with zinc-bound nitrite

Carbonic anhydrase (CA) is a zinc enzyme that catalyzes hydration of carbon dioxide (CO2) and dehydration of bicarbonate in red blood cells, thus facilitating CO2 transport and excretion. Bovine CA II may also react with nitrite to generate nitric oxide, although nitrite is a known inhibitor of the CO2 hydration reaction. To address the potential in vivo interference of these reactions and the nature of nitrite binding to the enzyme, we here investigate the inhibitory effect of 10-30 mM nitrite on Michaelis-Menten kinetics of CO2 hydration and bicarbonate dehydration by stopped-flow spectroscopy. Our data show that nitrite significantly affects the apparent dissociation constant KM for CO2 (11 mM) and bicarbonate (221 mM), and the turnover number kcat for the CO2 hydration (1.467 × 10(6) s(-1)) but not for the bicarbonate dehydration (7.927 × 10(5) s(-1)). These effects demonstrate mixed and competitive inhibition for the reaction with CO2 and bicarbonate, respectively, and are consistent with nitrite binding to the active site zinc. The high apparent dissociation constant found here for CO2, bicarbonate and nitrite (16-120 mM) are all overall consistent with published data and reveal a large capacity of free enzyme available for binding each of the three substrates at their in vivo levels, with little or no significant interference among reactions. The low affinity of the enzyme for nitrite suggests that the in vivo interaction between red blood cell CA II and nitrite requires compartmentalization at the anion exchanger protein of the red cell membrane to be physiologically relevant.

Crystal structure of the most catalytically effective carbonic anhydrase enzyme known, SazCA from the thermophilic bacterium Sulfurihydrogenibium azorense

Two thermostable α-carbonic anhydrases (α-CAs) isolated from thermophilic Sulfurihydrogenibium spp., namely SspCA (from S. yellowstonensis) and SazCA (from S. azorense), were shown in a previous work to possess interesting complementary properties. SspCA was shown to have an exceptional thermal stability, whereas SazCA demonstrated to be the most active α-CA known to date for the CO2 hydration reaction. Here we report the crystallographic structure of SazCA and the identification of the structural features responsible for its high catalytic activity, by comparing it with SspCA structure. These data are of relevance for the design of engineered proteins showing higher stability and catalytic activity than other α-CAs known to date.

Catalytic mechanism of α-class carbonic anhydrases: CO2 hydration and proton transfer

The carbonic anhydrases (CAs; EC 4.2.1.1) are a family of metalloenzymes that catalyze the reversible hydration of carbon dioxide (CO2) and dehydration of bicarbonate (HCO3 (-)) in a two-step ping-pong mechanism: [Formula: see text] CAs are ubiquitous enzymes and are categorized into five distinct classes (α, β, γ, δ and ζ). The α-class is found primarily in vertebrates (and the only class of CA in mammals), β is observed in higher plants and some prokaryotes, γ is present only in archaebacteria whereas the δ and ζ classes have only been observed in diatoms.The focus of this chapter is on α-CAs as the structure-function relationship is best understood for this class, in particular for humans. The reader is referred to other reviews for an overview of the structure and catalytic mechanism of the other CA classes. The overall catalytic site structure and geometry of α-CAs are described in the first section of this chapter followed by the kinetic studies, binding of CO2, and the proton shuttle network.

Water networks in fast proton transfer during catalysis by human carbonic anhydrase II

Variants of human carbonic anhydrase II (HCA II) with amino acid replacements at residues in contact with water molecules in the active-site cavity have provided insights into the proton transfer rates in this protein environment. X-ray crystallography and (18)O exchange measured by membrane inlet mass spectrometry have been used to investigate structural and catalytic properties of variants of HCA II containing replacements of Tyr7 with Phe (Y7F) and Asn67 with Gln (N67Q). The rate constants for transfer of a proton from His64 to the zinc-bound hydroxide during catalysis were 4 and 9 μs(-1) for Y7F and Y7F/N67Q, respectively, compared with a value of 0.8 μs(-1) for wild-type HCA II. These higher values observed for Y7F and Y7F/N67Q HCA II could not be explained by differences in the values of the pK(a) of the proton donor (His64) and acceptor (zinc-bound hydroxide) or by the orientation of the side chain of the proton shuttle residue His64. They appeared to be associated with a reduced level of branching in the networks of hydrogen-bonded water molecules between proton shuttle residue His64 and the zinc-bound solvent molecule as observed in crystal structures at 1.5-1.6 Å resolution. Moreover, Y7F/N67Q HCA II is unique among the variants studied in having a direct, hydrogen-bonded chain of water molecules between the zinc-bound solvent and N(ε) of His64. This study provides the clearest example to date of the relevance of ordered water structure to rate constants for proton transfer in catalysis by carbonic anhydrase.

Role of Trp19 and Tyr200 in catalysis by the γ-class carbonic anhydrase from Methanosarcina thermophila

Although widely distributed in Nature, only two γ class carbonic anhydrases are reported besides the founding member (Cam). Although roles for active-site residues important for catalysis have been identified in Cam, second shell residues have not been investigated. Two residues (Trp19 and Tyr200), positioned distant from the catalytic metal, were investigated by structural and kinetic analyses of replacement variants. Steady-state k(cat)/K(m) and k(cat) values decreased 3- to 10-fold for the Trp19 variants whereas the Y200 variants showed up to a 5-fold increase in k(cat). Rate constants for proton transfer decreased up to 10-fold for the Trp19 variants, and an increase of ~2-fold for Y200F. The pK(a) values for the proton donor decreased 1-2 pH units for Trp19 and Y200 variants. The variant structures revealed a loop composed of residues 62-64 that occupies a different conformation than previously reported. The results show that, although Trp19 and Y200 are non-essential, they contribute to an extended active-site structure distant from the catalytic metal that fine tunes catalysis. Trp19 is important for both CO(2)/bicarbonate interconversion, and the proton transfer step of catalysis.

Subunit III-depleted cytochrome c oxidase provides insight into the process of proton uptake by proteins

We review studies of subunit III-depleted cytochrome c oxidase (CcO III (-)) that elucidate the structural basis of steady-state proton uptake from solvent into an internal proton transfer pathway. The removal of subunit III from R. sphaeroides CcO makes proton uptake into the D pathway a rate-determining step, such that measurements of the pH dependence of steady-state O(2) consumption can be used to compare the rate and functional pK(a) of proton uptake by D pathways containing different initial proton acceptors. The removal of subunit III also promotes spontaneous suicide inactivation by CcO, greatly shortening its catalytic lifespan. Because the probability of suicide inactivation is controlled by the rate at which the D pathway delivers protons to the active site, measurements of catalytic lifespan provide a second method to compare the relative efficacy of proton uptake by engineered CcO III (-) forms. These simple experimental systems have been used to explore general questions of proton uptake by proteins, such as the functional value of an initial proton acceptor, whether an initial acceptor must be surface-exposed, which side chains will function as initial proton acceptors and whether multiple acceptors can speed proton uptake.

Proton delivery to ferryl heme in a heme peroxidase: enzymatic use of the Grotthuss mechanism

We test the hypothesized pathway by which protons are passed from the substrate, ascorbate, to the ferryl oxygen in the heme enzyme ascorbate peroxidase (APX). The role of amino acid side chains and bound solvent is demonstrated. We investigated solvent kinetic isotope effects (SKIE) for the wild-type enzyme and several site-directed replacements of the key residues which form the proposed proton path. Kinetic constants for H(2)O(2)-dependent enzyme oxidation to Compound I, k(1), and subsequent reduction of Compound II, k(3), were determined in steady-state assays by variation of both H(2)O(2) and ascorbate concentrations. A high value of the SKIE for wild type APX ((D)k(3) = 4.9) as well as a clear nonlinear dependence on the deuterium composition of the solvent in proton inventory experiments suggest the simultaneous participation of several protons in the transition state for proton transfer. The full SKIE and the proton inventory data were modeled by applying Gross-Butler-Swain-Kresge theory to a proton path inferred from the known structure of APX. The model has been tested by constructing and determining the X-ray structures of the R38K and R38A variants and accounts for their observed SKIEs. This work confirms APX uses two arginine residues in the proton path. Thus, Arg38 and Arg172 have dual roles, both in the formation of the ferryl species and binding of ascorbate respectively and to facilitate proton transfer between the two.

Identification of ligand binding sites of proteins using the Gaussian Network Model

The nonlocal nature of the protein-ligand binding problem is investigated via the Gaussian Network Model with which the residues lying along interaction pathways in a protein and the residues at the binding site are predicted. The predictions of the binding site residues are verified by using several benchmark systems where the topology of the unbound protein and the bound protein-ligand complex are known. Predictions are made on the unbound protein. Agreement of results with the bound complexes indicates that the information for binding resides in the unbound protein. Cliques that consist of three or more residues that are far apart along the primary structure but are in contact in the folded structure are shown to be important determinants of the binding problem. Comparison with known structures shows that the predictive capability of the method is significant.