Proton transfer function of carbonic anhydrase: Insights from QM/MM simulations

Recent QM/MM analyses of proton transfer function of human carbonic anhydrase II (CAII) are briefly reviewed. The topics include a preliminary analysis of nuclear quadrupole coupling constant calculations for the zinc ion and more detailed analyses of microscopic pK_a of the zinc-bound water and free energy profile for the proton transfer. From a methodological perspective, our results emphasize that performing sufficient sampling is essential to the calculation of all these quantities, which reflects the well solvated nature of CAII active site. From a mechanistic perspective, our analyses highlight the importance of electrostatics in shaping the energetics and kinetics of proton transfer in CAII for its function. We argue that once the pK_a for the zinc-bound water is modulated to be in the proper range (~ 7.0), proton transfer through a relatively well solvated cavity towards/from the protein surface (His64) does not require any major acceleration. Therefore, although structural details like the length of the water wire between the donor and acceptor groups still may make a non-negligible contribution, our computational results and the framework of analysis suggest that the significance of such “fine-tuning” is likely secondary to the modulation of pK_a of the zinc-bound water. We encourage further experimental analysis with mutation of (charged) residues not in the immediate neighborhood of the zinc ion to quantitatively test this electrostatics based framework; in particular, Φ analysis based on these mutations may shed further light into the relative importance of the classical Grotthus mechanism and the “proton hole” pathway that we have proposed recently for CAII.     

Proton transfer from exogenous donors in catalysis by human carbonic anhydrase II

In the site-specific mutant of human carbonic anhydrase in which the proton shuttle His64 is replaced with alanine, H64A HCA II, catalysis can be activated in a saturable manner by the proton donor 4-methylimidazole (4-MI). From 1H NMR relaxivities, we found 4-MI bound as a second-shell ligand of the tetrahedrally coordinated cobalt in Co(II)-substituted H64A HCA II, with 4-MI located about 4.5 A from the metal. Binding constants of 4-MI to H64A HCA II were estimated from: (1) NMR relaxation of the protons of 4-MI by Co(II)-H64A HCA II, (2) the visible absorption spectrum of Co(II)-H64A HCA II in the presence of 4-MI, (3) the inhibition by 4-MI of the catalytic hydration of CO2, and (4) from the catalyzed exchange of 18O between CO2 and water. These experiments along with previously reported crystallographic and catalytic data help identify a range of distances at which proton transfer is efficient in carbonic anhydrase II.     

Role of protein motions on proton transfer pathways in human carbonic anhydrase II

We report here a theoretical study on the formation of long-range proton transfer pathways in proteins due to side chain conformational fluctuations of amino acid residues and reorganization of interior hydration positions. The proton transfer pathways in such systems may be modeled as fluctuating hydrogen-bonded networks with both short- and long-lived connections between the networked nodes, the latter being formed by polar protein atoms and water molecules. It is known that these fluctuations may extend over several decades of time ranging from a few femtoseconds to a few milliseconds. We have shown in this article how the use of a variety of theoretical methods may be utilized to detect a generic set of pathways and assess the feasibility of forming one or more transient connections. We demonstrate the application of these methods to the enzyme human carbonic anhydrase II and its mutants. Our results reveal several alternative pathways in addition to the one mediated by His-64. We also probe at length the mechanism of key conformational fluctuations contributing to the formation of the detected pathways.     

Structure and catalysis by carbonic anhydrase II: role of active-site tryptophan 5

The tryptophan residue Trp5, highly conserved in the α class of carbonic anhydrases including human carbonic anhydrase II (HCA II), is positioned at the entrance of the active site cavity and forms a π-stacking interaction with the imidazole ring of the proton shuttle His64 in its outward orientation. We have observed that replacement of Trp5 in HCA II caused significant structural changes, as determined by X-ray diffraction, in the conformation of 11 residues at the N-terminus and in the orientation of the proton shuttle residue His64. Most significantly, two variants W5H and W5E HCA II had His64 predominantly outward in orientation, while W5F and wild type showed the superposition of both outward and inward orientations in crystal structures. Although Trp5 influences the orientation of the proton shuttle His64, this orientation had no significant effect on the rate constant for proton transfer near 1 μs⁻¹, determined by exchange of ¹⁸O between CO₂ and water measured by mass spectrometry. The apparent values of the pK_a of the zinc-bound water and the proton shuttle residue suggest that different active-site conformations influence the two stages of catalysis, the proton transfer stage and the interconversion of CO₂ and bicarbonate.     

Dynamic mechanism of proton transfer in mannitol 2-dehydrogenase from Pseudomonas fluorescens: mobile GLU292 controls proton relay through a water channel that connects the active site with bulk solvent

The active site of mannitol 2-dehydrogenase from Pseudomonas fluorescens (PfM2DH) is connected with bulk solvent through a narrow protein channel that shows structural resemblance to proton channels utilized by redox-driven proton pumps. A key element of the PfM2DH channel is the "mobile" Glu(292), which was seen crystallographically to adopt distinct positions up and down the channel. It was suggested that the "down → up" conformational change of Glu(292) could play a proton relay function in enzymatic catalysis, through direct proton shuttling by the Glu or because the channel is opened for water molecules forming a chain along which the protons flow. We report evidence from site-directed mutagenesis (Glu(292) → Ala) substantiated by data from molecular dynamics simulations that support a role for Glu(292) as a gate in a water chain (von Grotthuss-type) mechanism of proton translocation. Occupancy of the up and down position of Glu(292) is influenced by the bonding and charge state of the catalytic acid base Lys(295), suggesting that channel opening/closing motions of the Glu are synchronized to the reaction progress. Removal of gatekeeper control in the E292A mutant resulted in a selective, up to 120-fold slowing down of microscopic steps immediately preceding catalytic oxidation of mannitol, consistent with the notion that formation of the productive enzyme-NAD(+)-mannitol complex is promoted by a corresponding position change of Glu(292), which at physiological pH is associated with obligatory deprotonation of Lys(295) to solvent. These results underscore the important role of conformational dynamics in the proton transfer steps of alcohol dehydrogenase catalysis.     

Dual role of CO2/HCO3(-) buffer in the regulation of intracellular pH of three-dimensional tumor growths

Intracellular pH (pH(i)), a major modulator of cell function, is regulated by acid/base transport across membranes. Excess intracellular H(+) ions (e.g. produced by respiration) are extruded by transporters such as Na(+)/H(+) exchange, or neutralized by HCO(3)(-) taken up by carriers such as Na(+)-HCO(3)(-) cotransport. Using fluorescence pH(i) imaging, we show that cancer-derived cell lines (colorectal HCT116 and HT29, breast MDA-MB-468, pancreatic MiaPaca2, and cervical HeLa) extrude acid by H(+) efflux and HCO(3)(-) influx, largely sensitive to dimethylamiloride and 4,4'-diisothiocyanatostilbene-2,2'-disulfonate (DIDS), respectively. The magnitude of HCO(3)(-) influx was comparable among the cell lines and may represent a constitutive element of tumor pH(i) regulation. In contrast, H(+) efflux varied considerably (MDA-MB-468 > HCT116 > HT29 > MiaPaca2 > HeLa). When HCO(3)(-) flux was pharmacologically inhibited, acid extrusion in multicellular HT29 and HCT116 spheroids (～10,000 cells) was highly non-uniform and produced low pH(i) at the core. With depth, acid extrusion became relatively more DIDS-sensitive because the low extracellular pH at the spheroid core inhibits H(+) flux more than HCO(3)(-) flux. HCO(3)(-) flux inhibition also decelerated HCT116 spheroid growth. In the absence of CO(2)/HCO(3)(-), acid extrusion by H(+) flux in HCT116 and MDA-MB-468 spheroids became highly non-uniform and inadequate at the core. This is because H(+) transporters require extracellular mobile pH buffers, such as CO(2)/HCO(3)(-), to overcome low H(+) ion mobility and chaperone H(+) ions away from cells. CO(2)/HCO(3)(-) exerts a dual effect: as substrate for membrane-bound HCO(3)(-) transporters and as a mobile buffer for facilitating extracellular diffusion of H(+) ions extruded from cells. These processes can be augmented by carbonic anhydrase activity. We conclude that CO(2)/HCO(3)(-) is important for maintaining uniformly alkaline pH(i) in small, non-vascularized tumor growths and may be important for cancer disease progression.     

Kinetic and crystallographic studies of the role of tyrosine 7 in the active site of human carbonic anhydrase II

The rate limiting step in catalysis of bicarbonate dehydration by human carbonic anhydrase II (HCA II) is an intramolecular proton transfer from His64 to the zinc-bound hydroxide. We have examined the role of Tyr7 using site-specific mutagenesis and measuring catalysis by the ¹⁸O exchange method using membrane inlet mass spectrometry. The side chain of Tyr7 in HCA II extends into the active-site cavity about 7 Å from the catalytic zinc atom. Replacement of Tyr7 with eight other amino acids had no effect on the interconversion of bicarbonate and CO₂, but in some cases caused enhancements in the rate constant of proton transfer by nearly 10-fold. The variant Y7I HCA II enhanced intramolecular proton transfer approximately twofold; its structure was determined by X-ray crystallography at 1.5 Å resolution. No changes were observed in the ordered solvent structure in the active-site cavity or in the conformation of the side chain of the proton shuttle His64. However, the first 11 residues of the amino-terminal chain in Y7I HCA II assumed an alternate conformation compared with the wild type. Differential scanning calorimetry showed variants at position 7 had a melting temperature approximately 8 °C lower than that of the wild type.     

Drug protein interaction at the molecular level: A study of sulphonamide carbonic anhydrase complexes

The molecular mechanism of drug action has been studied by X-ray diffraction analysis of human carbonic anhydrase I complexed with two different sulphonamides. The acetazolamide and amino benzene sulphonamide are found to bind to the catalytically essential zinc ion thereby inhibiting the function of the enzyme. The inhibitor molecules are stabilized in the active site of the protein by van der Waals interaction with a number of protein side chain groups.     

Water-assisted proton transfer in ferredoxin I

The role of water molecules in assisting proton transfer (PT) is investigated for the proton-pumping protein ferredoxin I (FdI) from Azotobacter vinelandii. It was shown previously that individual water molecules can stabilize between Asp(15) and the buried [3Fe-4S](0) cluster and thus can potentially act as a proton relay in transferring H(+) from the protein to the μ(2) sulfur atom. Here, we generalize molecular mechanics with proton transfer to studying proton transfer reactions in the condensed phase. Both umbrella sampling simulations and electronic structure calculations suggest that the PT Asp(15)-COOH + H(2)O + [3Fe-4S](0) → Asp(15)-COO(-) + H(2)O + [3Fe-4S](0) H(+) is concerted, and no stable intermediate hydronium ion (H(3)O(+)) is expected. The free energy difference of 11.7 kcal/mol for the forward reaction is in good agreement with the experimental value (13.3 kcal/mol). For the reverse reaction (Asp(15)-COO(-) + H(2)O + [3Fe-4S](0)H(+) → Asp(15)-COOH + H(2)O + [3Fe-4S](0)), a larger barrier than for the forward reaction is correctly predicted, but it is quantitatively overestimated (23.1 kcal/mol from simulations versus 14.1 from experiment). Possible reasons for this discrepancy are discussed. Compared with the water-assisted process (ΔE ≈ 10 kcal/mol), water-unassisted proton transfer yields a considerably higher barrier of ΔE ≈ 35 kcal/mol.     

Squeezing at Entrance of Proton Transport Pathway in Proton-translocating Pyrophosphatase upon Substrate Binding

Homodimeric proton-translocating pyrophosphatase (H⁺-PPase; EC 3.6.1.1) is indispensable for many organisms in maintaining organellar pH homeostasis. This unique proton pump couples the hydrolysis of PP_i to proton translocation across the membrane. H⁺-PPase consists of 14–16 relatively hydrophobic transmembrane domains presumably for proton translocation and hydrophilic loops primarily embedding a catalytic site. Several highly conserved polar residues located at or near the entrance of the transport pathway in H⁺-PPase are essential for proton pumping activity. In this investigation single molecule FRET was employed to dissect the action at the pathway entrance in homodimeric Clostridium tetani H⁺-PPase upon ligand binding. The presence of the substrate analog, imidodiphosphate mediated two sites at the pathway entrance moving toward each other. Moreover, single molecule FRET analyses after the mutation at the first proton-carrying residue (Arg-169) demonstrated that conformational changes at the entrance are conceivably essential for the initial step of H⁺-PPase proton translocation. A working model is accordingly proposed to illustrate the squeeze at the entrance of the transport pathway in H⁺-PPase upon substrate binding.     

Molecular dynamics simulation and bioinformatics study on yeast aquaporin Aqy1 from Pichia pastoris

In the present study, an equilibrated system for the Aqy1 tetramer was developed, and molecular biophysics modeling showed that the Aqy1 channel was blocked by Tyr-31 in the N-terminus, which was also supported by the free energy profiles. However, bioinformatics analysis of the amino acid sequence of Aqy1 indicated this Tyr-31 is not conserved across all fungi. Analysis of the equilibrated structure showed that the central pore along the four-fold axis of the tetramers is formed with hydrophobic amino acid residues. In particular, Phe-90, Trp-198, and Phe-202 form the narrowest part of the pore. Therefore, water molecules are not expected to translocate through the central pore, a hypothesis that we confirmed by molecular dynamics simulations. Each monomer of the Aqy1 tetramers forms a channel whose walls consist mostly of hydrophilic residues, transporting through the selectivity filter containing Arg-227, His-212, Phe-92, and Ala-221, and the two conserved Asn-Pro-Ala (NPA) motifs containing asparagines 224 and 112. In summary, not all fungal aquaporins share the same gating mechanism by a tyrosine residue in the N-terminus, and the structural analysis in the present study should aid our understanding of aquaporin structure and its functional implications.     

Comparative physiology and molecular analysis of carbonic anhydrase from the red blood cells of teleost fish

This study investigates the early evolution of vertebrate red blood cell (rbc) carbonic anhydrase (CA) by examining the physiological and molecular properties of rbc CA in teleost fish. When representatives of four different families of teleosts were compared, it was found that differences in overall rbc CA activity were due to different concentrations of CA, rather than differences in the enzymes kinetic properties. Additional molecular analysis of CA from the rbcs of rainbow trout provided further evidence that critical elements of the enzyme, such as the active site, have been highly conserved during vertebrate evolution. The active site of the trout CA differed from that of gar rbc CA at only two amino acid positions. The rainbow trout rbc CA sequence also showed high sequence homology with CA sequences from other fish tissues, and fits into an emerging group of fish CAs that are basal to mammalian CA I, II and III. Northern blot analysis of the tissue expression of the sequenced CA indicated that it is primarily found in the rbcs, but high amounts of cytosolic CA activity were also found in the gill, suggesting the presence of other cytosolic CA isozymes in this species.Abbreviations Az acetazolamide - CA carbonic anhydrase - MP maximum parsimony - NJ neighbour joining - RACE rapid amplification of cDNA ends - rbc red blood cellCommunicated by L.C.-H. Wang     

Mitochondrial carbonic anhydrase in the nervous system: expression in neuronal and glial cells

Carbonic anhydrase (CA) V is a mitochondrial enzyme that has been reported in several tissues of the gastrointestinal tract. In liver, it participates in ureagenesis and gluconeogenesis by providing bicarbonate ions for two other mitochondrial enzymes: carbamyl phosphate synthetase I and pyruvate carboxylase. This study presents evidence of immunohistochemical localization of CA V in the rodent nervous tissue. Polyclonal rabbit antisera against a polypeptide of 17 C-terminal amino acids of rat CA V and against purified recombinant mouse isozyme were used in western blotting and immunoperoxidase stainings. Immunohistochemistry showed that CA V is expressed in astrocytes and neurons but not in oligodendrocytes, which are rich in CA II, or capillary endothelial cells, which express CA IV on their plasma face. The specificity of the immunohistochemical results was confirmed by western blotting, which identified a major 30-kDa polypeptide band of CA V in mouse cerebral cortex, hippocampus, cerebellum, spinal cord, and sciatic nerve. The expression of CA V in astrocytes and neurons suggests that this isozyme has a cell-specific, physiological role in the nervous system. In astrocytes, CA V may play an important role in gluconeogenesis by providing bicarbonate ions for the pyruvate carboxylase. The neuronal CA V could be involved in the regulation of the intramitochondrial calcium level, thus contributing to the stability of the intracellular calcium concentration. CA V may also participate in bicarbonate ion-induced GABA responses by regulating the bicarbonate homeostasis in neurons, and its inhibition could be the basis of some neurotropic effects of carbonic anhydrase inhibitors.     

Esterase D and carbonic anhydrase 2 in a natural population ofCebus apella from Brazil

The electrophoretic patterns of esterase D (ESD; E.C.3.1.1.1) and carbonic anhydrase 2 (CA 2, E.C.4.2.1.1) were studied in 147 specimens ofCebus apella. Three phenotypes were detected at the esterase D system,ESD ¹ allele showing a frequency of 44%, markedly different from those observed in Old World monkeys. CA2 also proved to be polymorphic, with three alleles being detected at the following frequencies:CA2 ¹, 98%;CA2 ² andCA2 ³, both 1%. The CA2 activity was absent in newborn animals and in fetuses.     

Proton transport behavior through the influenza A M2 channel: insights from molecular simulation

The structural properties of the influenza A virus M2 transmembrane channel in dimyristoylphosphatidylcholine bilayer for each of the four protonation states of the proton-gating His-37 tetrad and their effects on proton transport for this low-pH activated, highly proton-selective channel are studied by classical molecular dynamics with the multistate empirical valence-bond (MS-EVB) methodology. The excess proton permeation free energy profile and maximum ion conductance calculated from the MS-EVB simulation data combined with the Poisson-Nernst-Planck theory indicates that the triply protonated His-37 state is the most likely open state via a significant side-chain conformational change of the His-37 tetrad. This proposed open state of M2 has a calculated proton permeation free energy barrier of 7 kcal/mol and a maximum conductance of 53 pS compared to the experimental value of 6 pS. By contrast, the maximum conductance for Na(+) is calculated to be four orders of magnitude lower, in reasonable agreement with the experimentally observed proton selectivity. The pH value to activate the channel opening is estimated to be 5.5 from dielectric continuum theory, which is also consistent with experimental results. This study further reveals that the Ala-29 residue region is the primary binding site for the antiflu drug amantadine (AMT), probably because that domain is relatively spacious and hydrophobic. The presence of AMT is calculated to reduce the proton conductance by 99.8% due to a significant dehydration penalty of the excess proton in the vicinity of the channel-bound AMT.     

Identification of potential tumour-associated carbonic anhydrase isozyme IX inhibitors: atom-based 3D-QSAR modelling,pharmacophore-based virtual screening and molecular docking studies

Abstract

Tumour hypoxia results in dramatic changes in the gene expression, proliferation and survival of tumour cells. The tumour cells shift towards anaerobic glycolysis which results in change of pH in their microenvironment. In response to this stress, over expression of carbonic anhydrase IX (CA IX) genes is observed in many solid tumours. So, selective inhibition of CA IX can be a promising target for anti-cancer drugs. In this work in silico tools like atom-based 3D-QSAR modelling, pharmacophore-based virtual screening and molecular docking were used to identify potential CA IX inhibitors. Based on the training set used in the QSAR model, twenty pharmacophore models were generated. Out of these, HHHR_1, AHHR_1, DHHHR_1, AHHHR_1 model was used to screen a database of 1,50,000 compounds retrieved from ZINC 15 database. R² and Q² was 0.9864 and 0.8799, respectively, for the developed QSAR model. 163 compounds showed a phase screen score above 2.4 in which ZINC02260669 was the highest ranked (screen score, 2.852058) compound in all the four models. Built QSAR model was used to predict the activity of all these 163 compounds and ZINC72370966 showed the highest predicted activity with pKi value of 7.649. These compounds were docked against CA IX (human) protein (PDB ID 5FL6) and molecular docking results showed favourable binding interactions for the best ten identified hits. This work gives design insights and some potential scaffolds which can be developed as CA IX inhibitors.

Communicated by Ramaswamy H. Sarma     

Multiscale simulation in polymer science

The paper gives a short overview on recent approaches to link several time and length scales in soft matter simulations. Special attention is given to the fact that in contrast to low molecular weight compounds, intramolecular entropy is as important as the energetic contribution. First applications to industrially relevant problems are mentioned.     

Multiple sources of carbonic anhydrase activity in pea thylakoids: soluble and membrane-bound forms

Carbonic anhydrase (CA) activity of pea thylakoids, thylakoid membranes enriched with photosystem I (PSI-membranes), or photosystem II (PSII-membranes) as well as both supernatant and pellet after precipitation of thylakoids treated with detergent Triton X-100 were studied. CA activity of thylakoids in the presence of varying concentrations of Triton X-100 had two maxima, at Triton/chlorophyll (triton/Chl) ratios of 0.3 and 1.0. CA activities of PSI-membranes and PSII-membranes had only one maximum each, at Triton/Chl ratio 0.3 or 1.0, respectively. Two CAs with characteristics of the membrane-bound proteins and one CA with characteristics of the soluble proteins were found in the medium after thylakoids were incubated with Triton. One of the first two CAs had mobility in PAAG after native electrophoresis the same as that of CA residing in PSI-membranes, and the other CA had mobility the same as the mobility of CA residing in PSII-membranes, but the latter was different from CA situated in PSII core-complex (Ignatova et al. 2006 Biochemistry (Moscow) 71:525–532). The properties of the “soluble” CA removed from thylakoids were different from the properties of the known soluble CAs of plant cell: apparent molecular mass was about 262 kD and it was three orders more sensitive to the specific CA inhibitor, ethoxyzolamide, than soluble stromal CA. The data are discussed as indicating the presence of, at least, four CAs in pea thylakoids.     

Classical Molecular Dynamics Simulation of Kappa Squared Factor in Resonance Energy Transfer for Linear Dipole Models

Molecular dynamics simulations of linear models interacting through a dipolar Kihara intermolecular potential are presented. Molecular orientation correlations are used to calculate the orientational factor kappa squared in the resonance energy transfer (RET) as a function of the intermolecular separation. The distance, R 0 (2/3), at which the simulated systems show an isotropic behavior is calculated and an analysis of the dependence of R 0 (2/3) on microscopic properties (molecular aspect ratio and dipole moment) as well on thermodynamics (temperature and density) is presented. An explanation of the use of metallic cations as probes in RET is given and some relations of our models with biological molecules are pointed out.     

Proton and sucrose transport in isolated tonoplast vesicles from sugarcane stalk tissue

Tonoplast vesicles prepared from immature sugarcane ( Saccharum spp., hybrid cv. H65–7052) tissue and purified on a discontinuous dextran gradient take up sucrose. Uptake was stimulated by MgATP. Evidence that the mechanism is linked to proton transport is derived from "pH jump'data and from inhibition of ATP-stimulated sucrose transport by the protonophore carbonyl cyanide m -chlorophenylhydrazone (CCCP) and by the proton-channel blocker of proton-linked ATPases. N. N '-dicyclo-hexylcarbodiimide (DCCD). A saturable phase of sucrose uptake was found at low substrate concentrations, and a linear phase characterized uptake at higher concentrations. Uptake was specific for sucrose, as demonstrated by competition experiments with various sugars. Sucrose uptake by the vesicle fraction was inhibited by KNO₃, protonophores and protein modifying reagents, whereas sodium orthovanadate had no effect. Overall, the evidence suggests an ATP-hydrolysis-dependent tonoplasl antiport for sucrose transport, although a more direct influence of ATP on conformational changes in relevant tonoplast proteins cannot be ruled out.