Elucidation of the proton transport mechanism in human carbonic anhydrase II

Insight into the activation mechanism of carbonic anhydrase(II) through 2-(2-aminoethyl)-pyridine: a promising pathway for enhanced enzymatic activity

Activation of human carbonic anhydrase II (hCA II) holds great promise for treating memory loss symptoms associated with Alzheimer's disease. Despite its importance, the activation mechanism of hCA II has been largely overlooked in favor of the well-studied inhibition mechanism. To address this unexplored realm, we use first-principles calculations to tease out the activation mechanism of hCA II using 2-(2-aminoethyl)-pyridine (2-2AEPy), a promising in vitro activator. We explored both stepwise and concerted mechanisms via both available nitrogen sites of 2-2AEPy: (i) aminoethyl group (Nα) and (ii) pyridine ring (Nβ). Our results show that a concerted mechanism via Nα holds the key to hCA II activation. The activation process of the concerted mechanism exhibits the characteristics of an exergonic reaction, wherein the transition state resembles the reactant with a notably low imaginary frequency of 452.4i cm-1 and barrier height of 5.2 kcal mol-1. Such meager transition barriers propel the activation of hCA II at in vivo temperatures. These findings initiate future research into hCA II activation mechanisms and the development of efficient activators, which may lead to promising therapeutic interventions for Alzheimer's disease.

XFEL structure of carbonic anhydrase II: a comparative study of XFEL, NMR, X-ray and neutron structures

The combination of X-ray free-electron lasers (XFELs) with serial femtosecond crystallography represents cutting-edge technology in structural biology, allowing the study of enzyme reactions and dynamics in real time through the generation of `molecular movies'. This technology combines short and precise high-energy X-ray exposure to a stream of protein microcrystals. Here, the XFEL structure of carbonic anhydrase II, a ubiquitous enzyme responsible for the interconversion of CO2 and bicarbonate, is reported, and is compared with previously reported NMR and synchrotron X-ray and neutron single-crystal structures.

Water, Protons, and the Gating of Voltage-Gated Potassium Channels

Ion channels are ubiquitous throughout all forms of life. Potassium channels are even found in viruses. Every cell must communicate with its surroundings, so all cells have them, and excitable cells, in particular, especially nerve cells, depend on the behavior of these channels. Every channel must be open at the appropriate time, and only then, so that each channel opens in response to the stimulus that tells that channel to open. One set of channels, including those in nerve cells, responds to voltage. There is a standard model for the gating of these channels that has a section of the protein moving in response to the voltage. However, there is evidence that protons are moving, rather than protein. Water is critical as part of the gating process, although it is hard to see how this works in the standard model. Here, we review the extensive evidence of the importance of the role of water and protons in gating these channels. Our principal example, but by no means the only example, will be the Kv1.2 channel. Evidence comes from the effects of D2O, from mutations in the voltage sensing domain, as well as in the linker between that domain and the gate, and at the gate itself. There is additional evidence from computations, especially quantum calculations. Structural evidence comes from X-ray studies. The hydration of ions is critical in the transfer of ions in constricted spaces, such as the gate region and the pore of a channel; we will see how the structure of the hydrated ion fits with the structure of the channel. In addition, there is macroscopic evidence from osmotic experiments and streaming current measurements. The combined evidence is discussed in the context of a model that emphasizes the role of protons and water in gating these channels.

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.

Structure and Dynamics of the Isozymes II and IX of Human Carbonic Anhydrase

Human carbonic anhydrases (HCAs) are responsible for the pH control and sensing in our body and constitute key components in the central pH paradigm connected to cancer therapeutics. However, little or no molecular level studies are available on the pH-dependent stability and functional dynamics of the known isozymes of HCA. The main objective of this Article is to report the first bench-marking study on the structure and dynamics of the two most efficient isozymes, HCA II and IX, at neutral pH using classical molecular dynamics (MD) and constant pH MD (CpHMD) simulations combined with umbrella sampling, transition path sampling, and Markov state models. Starting from the known crystal structures of HCA II and the monomeric catalytic domain of HCA IX (labeled as HCA IX-c), we have generated classical MD and CpHMD trajectories (of length 1 μs each). In all cases, the overall stability, RMSD, and secondary structure segments of the two isozymes are found to be quite similar. Functionally important dynamics of these two enzymes have been probed in terms of active site hydration, coordination of the Zn(II) ion to a transient excess water, and the formation of putative proton transfer paths. The most important difference between the two isozymes is observed for the side-chain fluctuations of His-64 that is expected to shuttle an excess proton out of the active site as a part of the rate-determining intramolecular proton transfer reaction. The relative stability of the stable inward and outward conformations of the His-64 side-chain and the underlying free energy surfaces are found to depend strongly on the isozyme. In each case, a lower free energy barrier is detected between predominantly inward conformations from predominantly outward ones when simulated under constant pH conditions. The kinetic rate constants of interconversion between different free energy basins are found to span 10⁷-10⁸ s^-1 with faster conformational transitions predicted at constant pH condition. The estimated rate constants and free energies are expected to validate if the fluctuation of the His-64 side-chain in HCA IX may have a significance similar to that known in the multistep catalytic cycle of HCA II.

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.

Towards solvent regulated self-activation of N-terminal disulfide bond oxidoreductase-D

N-terminal disulfide bond oxidoreductase-D (nDsbD), an essential redox enzyme in Gram-negative bacteria, consists of a single disulfide bond (Cys₁₀₃-Cys₁₀₉) in its active site. The enzymatic functions are believed to be regulated by an electron transfer mediated redox switching of the disulfide bond, which is vital in controlling bacterial virulence factors. In light of the disulfide bond's inclination towards nucleophilic cleavage, it is also plausible that an internal nucleophile could second the existing electron transfer mechanism in nDsbD. Using QM/MM MD metadynamics simulations, we explore different possibilities of generating an internal nucleophile near the nDsbD active site, which could serve as a fail-over mechanism in cleaving the disulfide bond. The simulations show the formation of the internal nucleophile Tyr₄₂O^- (F ≈ 9 kcal mol^-1) and its stabilization through the solvent medium. The static gas-phase calculations show that Tyr₄₂O^- could be a potential nucleophile for cleaving the S-S bond. Most strikingly, it is also seen that Tyr₄₂O^- and Asp₆₈OH communicate with each other through a proton-hole like water wire (F ≈ 12 kcal mol^-1), thus modulating the nucleophile formation. Accordingly, we propose the role of a solvent in regulating the internal nucleophilic reactions and the subsequent self-activation of nDsbD. We believe that this could be deterministic while designing enzyme-targeted inhibitor compounds.

Nonlinear Reaction Coordinate of an Enzyme Catalyzed Proton Transfer Reaction

We present an in-depth study on the theoretical calculation of an optimum reaction coordinate as a linear or nonlinear combination of important collective variables (CVs) sampled from an ensemble of reactive transition paths for an intramolecular proton transfer reaction catalyzed by the enzyme human carbonic anhydrase (HCA) II. The linear models are optimized by likelihood maximization for a given number of CVs. The nonlinear models are based on an artificial neural network with the same number of CVs and optimized by minimizing the root-mean-square error in comparison to a training set of committor estimators generated for the given transition. The nonlinear reaction coordinate thus obtained yields the free energy of activation and rate constant as 9.46 kcal mol^-1 and 1.25 × 10⁶ s^-1, respectively. These estimates are found to be in quantitative agreement with the known experimental results. We have also used an extended autoencoder model to show that a similar analysis can be carried out using a single CV only. The resultant free energies and kinetics of the reaction slightly overestimate the experimental data. The implications of these results are discussed using a detailed microkinetic scheme of the proton transfer reaction catalyzed by HCA II.

Hydrogen-deuterium exchange mass spectrometry reveals three unique binding responses of mAbs directed to the catalytic domain of hCAIX

Human carbonic anhydrase (hCAIX), an extracellular enzyme that catalyzes the reversible hydration of CO₂, is often overexpressed in solid tumors. This enzyme is instrumental in maintaining the survival of cancer cells in a hypoxic and acidic tumor microenvironment. Absent in most normal tissues, hCAIX is a promising therapeutic target for detection and treatment of solid tumors. Screening of a library of anti-hCAIX monoclonal antibodies (mAbs) previously identified three therapeutic candidates (mAb c2C7, m4A2 and m9B6) with distinct biophysical and functional characteristics. Selective binding to the catalytic domain was confirmed by yeast surface display and isothermal calorimetry, and deeper insight into the dynamic binding profiles of these mAbs upon binding were highlighted by bottom-up hydrogen-deuterium exchange mass spectrometry (HDX-MS). Here, a conformational and allosterically silent epitope was identified for the antibody-drug conjugate candidate c2C7. Unique binding profiles are described for both inhibitory antibodies, m4A2 and m9B6. M4A2 reduces the ability of the enzyme to hydrate CO₂ by steric gating at the entrance of the catalytic cavity. Conversely, m9B6 disrupts the secondary structure that is necessary for substrate binding and hydration. The synergy of these two inhibitory mechanisms is demonstrated in in vitro activity assays and HDX-MS. Finally, the ability of m4A2 to modulate extracellular pH and intracellular metabolism is reported. By highlighting three unique modes by which hCAIX can be targeted, this study demonstrates both the utility of HDX-MS as an important tool in the characterization of anti-cancer biotherapeutics, and the underlying value of CAIX as a therapeutic target.

Post-translational modifications in tumor-associated carbonic anhydrases

Human carbonic anhydrases IX (hCA IX) and XII (hCA XII) are two proteins associated with tumor formation and development. These enzymes have been largely investigated both from a biochemical and a functional point of view. However, limited data are currently available on the characterization of their post-translational modifications (PTMs) and the functional implication of these structural changes in the tumor environment. In this review, we summarize existing literature data on PTMs of hCA IX and hCA XII, such as disulphide bond formation, phosphorylation, O-/N-linked glycosylation, acetylation and ubiquitination, highlighting, when possible, their specific role in cancer pathological processes.

Conformational Change of H64 and Substrate Transportation: Insight Into a Full Picture of Enzymatic Hydration of CO2 by Carbonic Anhydrase

The enzymatic hydration of CO₂ into HCO₃ ^- by carbonic anhydrase (CA) is highly efficient and environment-friendly measure for CO₂ sequestration. Here extensive MM MD and QM/MM MD simulations were used to explore the whole enzymatic process, and a full picture of the enzymatic hydration of CO₂ by CA was achieved. Prior to CO₂ hydration, the proton transfer from the water molecule (WT1) to H64 is the rate-limiting step with the free energy barrier of 10.4 kcal/mol, which leads to the ready state with the Zn-bound OH^-. The nucleophilic attack of OH^- on CO₂ produces HCO₃ ^- with the free energy barrier of 4.4 kcal/mol and the free energy release of about 8.0 kcal/mol. Q92 as the key residue manipulates both CO₂ transportation to the active site and release of HCO₃ ^-. The unprotonated H64 in CA prefers in an inward orientation, while the outward conformation is favorable energetically for its protonated counterpart. The conformational transition of H64 between inward and outward correlates with its protonation state, which is mediated by the proton transfer and the product release. The whole enzymatic cycle has the free energy span of 10.4 kcal/mol for the initial proton transfer step and the free energy change of -6.5 kcal/mol. The mechanistic details provide a comprehensive understanding of the entire reversible conversion of CO₂ into bicarbonate and roles of key residues in chemical and nonchemical steps for the enzymatic hydration of CO₂.

Cooperative protein-solvent tuning of proton transfer energetics: carbonic anhydrase as a case study

We investigate the coupling between the proton transfer (PT) energetics and the protein-solvent dynamics using the intra-molecular PT in wild type (wt) human carbonic anhydrase II and its ten-fold faster mutant Y7F/N67Q as a test case. We calculate the energy variation upon PT, and from that we also calculate the PT reaction free energy, making use of a hybrid quantum mechanics/molecular dynamics approach. In agreement with the experimental data, we obtain that the reaction free energy is basically the same in the two systems. Yet, we show that the instantaneous PT energy is on average lower in the mutant possibly contributing to the faster PT rate. Analysis of the contribution to the PT energetics of the solvent and of each protein residue, also not in the vicinity of the active site, provides evidence for electrostatic tuning of the PT energy arising from the combined effect of the solvent and the protein environment. These findings open up a way to the more general task of the rational design of mutants with either enhanced or reduced PT rate.

Structural and biochemical characterization of novel carbonic anhydrases from Phaeodactylum tricornutum

Carbonic anhydrases (CAs) are a well characterized family of metalloenzymes that are highly efficient in facilitating the interconversion between carbon dioxide and bicarbonate. Recently, CA activity has been associated with the LCIB (limiting CO₂-inducible protein B) protein family, which has been an interesting target in aquatic photosynthetic microorganisms. To gain further insight into the catalytic mechanism of this new group of CAs, the X-ray structure of a highly active LCIB homolog (PtLCIB3) from the diatom Phaeodactylum tricornutum was determined. The CA activities of PtLCIB3, its paralog PtLCIB4 and a variety of their mutants were also measured. It was discovered that PtLCIB3 has a classic β-CA fold and its overall structure is highly similar to that of its homolog PtLCIB4. Subtle structural alterations between PtLCIB3 and PtLCIB4 indicate that an alternative proton-shuttle cavity could perhaps be one reason for their remarkable difference in CA activity. A potential alternative proton-shuttle route in the LCIB protein family is suggested based on these results.

Structure and mechanism of copper-carbonic anhydrase II: a nitrite reductase

Nitric oxide (NO) promotes vasodilation through the activation of guanylate cyclase, resulting in the relaxation of the smooth muscle vasculature and a subsequent decrease in blood pressure. Therefore, its regulation is of interest for the treatment and prevention of heart disease. An example is pulmonary hypertension which is treated by targeting this NO/vasodilation pathway. In bacteria, plants and fungi, nitrite (NO2 -) is utilized as a source of NO through enzymes known as nitrite reductases. These enzymes reduce NO2 - to NO through a catalytic metal ion, often copper. Recently, several studies have shown nitrite reductase activity of mammalian carbonic anhydrase II (CAII), yet the molecular basis for this activity is unknown. Here we report the crystal structure of copper-bound human CAII (Cu-CAII) in complex with NO2 - at 1.2 Å resolution. The structure exhibits Type 1 (T-1) and 2 (T-2) copper centers, analogous to bacterial nitrite reductases, both required for catalysis. The copper-substituted CAII active site is penta-coordinated with a 'side-on' bound NO2 -, resembling a T-2 center. At the N terminus, several residues that are normally disordered form a porphyrin ring-like configuration surrounding a second copper, acting as a T-1 center. A structural comparison with both apo- (without metal) and zinc-bound CAII (Zn-CAII) provides a mechanistic picture of how, in the presence of copper, CAII, with minimal conformational changes, can function as a nitrite reductase.

Mechanism of Action of Non-Synonymous Single Nucleotide Variations Associated with α-Carbonic Anhydrase II Deficiency

Human carbonic anhydrase II (CA-II) is a Zinc (Zn 2 + ) metalloenzyme responsible for maintenance of acid-base balance within the body through the reversible hydration of CO 2 to produce protons (H + ) and bicarbonate (BCT). Due to its importance, alterations to the amino acid sequence of the protein as a result of single nucleotide variations (nsSNVs) have detrimental effects on homeostasis. Six pathogenic CA-II nsSNVs, K18E, K18Q, H107Y, P236H, P236R and N252D were identified, and variant protein models calculated using homology modeling. The effect of each nsSNV was analyzed using motif analysis, molecular dynamics (MD) simulations, principal component (PCA) and dynamic residue network (DRN) analysis. Motif analysis identified 11 functionally important motifs in CA-II. RMSD data indicated subtle SNV effects, while PCA analysis revealed that the presence of BCT results in greater conformational sampling and free energy in proteins. DRN analysis showed variant allosteric effects, and the average betweenness centrality (BC) calculations identified Glu117 as the most important residue for communication in CA-II. The presence of BCT was associated with a reduction to Glu117 usage in all variants, suggesting implications for Zn 2 + dissociation from the CA-II active site. In addition, reductions to Glu117 usage are associated with increases in the usage of the primary and secondary Zn 2 + ligands; His94, His96, His119 and Asn243 highlighting potential compensatory mechanisms to maintain Zn 2 + within the active site. Compared to traditional MD simulation investigation, DRN analysis provided greater insights into SNV mechanism of action, indicating its importance for the study of missense mutation effects in proteins and, in broader terms, precision medicine related research.

Kinetic Study of CO2 Hydration by Small-Molecule Catalysts with A Second Coordination Sphere that Mimic the Effect of the Thr-199 Residue of Carbonic Anhydrase

Zinc complexes were synthesized as catalysts that mimic the ability of carbonic anhydrase (CA) for the CO2 hydration reaction (H2O + CO2 → H+ + HCO3-). For these complexes, a tris(2-pyridylmethyl)amine (TPA) ligand mimicking only the active site, and a 6-((bis(pyridin-2-ylmethyl)amino)methyl)pyridin-2-ol (TPA-OH) ligand mimicking the hydrogen-bonding network of the secondary coordination sphere of CA were used. Potentiometric pH titration was used to determine the deprotonation ability of the Zn complexes, and their pKa values were found to be 8.0 and 6.8, respectively. Stopped-flow spectrophotometry was used to confirm the CO2 hydration rate. The rate constants were measured to be 648.4 and 730.6 M-1s-1, respectively. The low pKa value was attributed to the hydrogen-bonding network of the secondary coordination sphere of the catalyst that mimics the behavior of CA, and this was found to increase the CO2 hydration rate of the catalyst.

Kinetic study of catalytic CO2 hydration by metal-substituted biomimetic carbonic anhydrase model complexes

The rapid rise of the CO₂ level in the atmosphere has spurred the development of CO₂ capture methods such as the use of biomimetic complexes that mimic carbonic anhydrase. In this study, model complexes with tris(2-pyridylmethyl)amine (TPA) were synthesized using various transition metals (Zn²⁺, Cu²⁺ and Ni²⁺) to control the intrinsic proton-donating ability. The pK_a of the water coordinated to the metal, which indicates its proton-donating ability, was determined by potentiometric pH titration and found to increase in the order [(TPA)Cu(OH₂)]²⁺ < [(TPA)Ni(OH₂)]²⁺ < [(TPA)Zn(OH₂)]²⁺. The effect of pK_a on the CO₂ hydration rate was investigated by stopped-flow spectrophotometry. Because the water ligand in [(TPA)Zn(OH₂)]²⁺ had the highest pK_a, it would be more difficult to deprotonate it than those coordinated to Cu²⁺ and Ni²⁺. It was, therefore, expected that the complex would have the slowest rate for the reaction of the deprotonated water with CO₂ to form bicarbonate. However, it was confirmed that [(TPA)Zn(OH₂)]²⁺ had the fastest CO₂ hydration rate because the substitution of bicarbonate with water (bicarbonate release) occurred easily.

Orthogonal order parameters to model the reaction coordinate of an enzyme catalyzed reaction

The choice of suitable collective variables in formulating an optimal reaction coordinate is a challenging task for activated transitions between a pair of stable states especially when dealing with biochemical changes such as enzyme catalyzed reactions. A detailed benchmarking study is carried out on the choice of collective variables that can distinguish between the stable states unambiguously. We specifically address the issue if these variables may be directly used to model the optimal reaction coordinate, or if it would be better to use their orthogonalized counterparts. The proposed computational scheme is applied to the rate determining intramolecular proton transfer step in the enzyme human carbonic anhydrase II. The optimum reaction coordinate is determined with and without orthogonalization of the collective variables pertinent to a key conformational fluctuation and the actual proton transfer step at the active site of the enzyme. Suitability of the predicted reaction coordinates in different processes is examined in terms of the free energy profile projected along the reaction coordinate, the rate constant of transition and the underlying molecular mechanism of barrier crossing. Our results indicate that a better agreement with earlier simulation and experimental data is obtained when the orthogonalized collective variables are used to model the reaction coordinate.

Crystallography and Its Impact on Carbonic Anhydrase Research

X-ray and neutron crystallography are powerful techniques utilized to study the structures of biomolecules. Visualization of enzymes in complex with substrate/product and the capture of intermediate states can be related to activity to facilitate understanding of the catalytic mechanism. Subsequent analysis of small molecule binding within the enzyme active site provides insight into mechanisms of inhibition, supporting the design of novel inhibitors using a structure-guided approach. The first X-ray crystal structures were determined for small, ubiquitous enzymes such as carbonic anhydrase (CA). CAs are a family of zinc metalloenzymes that catalyze the hydration of CO2, producing HCO3 - and a proton. The CA structure and ping-pong mechanism have been extensively studied and are well understood. Though the function of CA plays an important role in a variety of physiological functions, CA has also been associated with diseases such as glaucoma, edema, epilepsy, obesity, and cancer and is therefore recognized as a drug target. In this review, a brief history of crystallography and its impact on CA research is discussed.

Analysis of Phosphoryl-Transfer Enzymes with QM/MM Free Energy Simulations

We discuss the application of quantum mechanics/molecular mechanics (QM/MM) free energy simulations to the analysis of phosphoryl transfers catalyzed by two enzymes: alkaline phosphatase and myosin. We focus on the nature of the transition state and the issue of mechanochemical coupling, respectively, in the two enzymes. The results illustrate unique insights that emerged from the QM/MM simulations, especially concerning the interpretation of experimental data regarding the nature of enzymatic transition states and coupling between global structural transition and catalysis in the active site. We also highlight a number of technical issues worthy of attention when applying QM/MM free energy simulations, and comment on a number of technical and mechanistic issues that require further studies.

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.

Inhibition of carbonic anhydrase IX targets primary tumors, metastases, and cancer stem cells: Three for the price of one

Human carbonic anhydrase (CA) IX is a tumor-associated protein, since it is scarcely present in normal tissues, but highly overexpressed in a large number of solid tumors, where it actively contributes to survival and metastatic spread of tumor cells. Due to these features, the characterization of its biochemical, structural, and functional features for drug design purposes has been extensively carried out, with consequent development of several highly selective small molecule inhibitors and monoclonal antibodies to be used for different purposes. Aim of this review is to provide a comprehensive state-of-the-art of studies performed on this enzyme, regarding structural, functional, and biomedical aspects, as well as the development of molecules with diagnostic and therapeutic applications for cancer treatment. A brief description of additional pharmacologic applications for CA IX inhibition in other diseases, such as arthritis and ischemia, is also provided.

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.

Metal Ion Modeling Using Classical Mechanics

Metal ions play significant roles in numerous fields including chemistry, geochemistry, biochemistry, and materials science. With computational tools increasingly becoming important in chemical research, methods have emerged to effectively face the challenge of modeling metal ions in the gas, aqueous, and solid phases. Herein, we review both quantum and classical modeling strategies for metal ion-containing systems that have been developed over the past few decades. This Review focuses on classical metal ion modeling based on unpolarized models (including the nonbonded, bonded, cationic dummy atom, and combined models), polarizable models (e.g., the fluctuating charge, Drude oscillator, and the induced dipole models), the angular overlap model, and valence bond-based models. Quantum mechanical studies of metal ion-containing systems at the semiempirical, ab initio, and density functional levels of theory are reviewed as well with a particular focus on how these methods inform classical modeling efforts. Finally, conclusions and future prospects and directions are offered that will further enhance the classical modeling of metal ion-containing systems.

Coupling Protein Dynamics with Proton Transport in Human Carbonic Anhydrase II

The role of protein dynamics in enzyme catalysis is one of the most highly debated topics in enzymology. The main controversy centers around what may be defined as functionally significant conformational fluctuations and how, if at all, these fluctuations couple to enzyme catalyzed events. To shed light on this debate, the conformational dynamics along the transition path surmounting the highest free energy barrier have been herein investigated for the rate limiting proton transport event in human carbonic anhydrase (HCA) II. Special attention has been placed on whether the motion of an excess proton is correlated with fluctuations in the surrounding protein and solvent matrix, which may be rare on the picosecond and subpicosecond time scales of molecular motions. It is found that several active site residues, which do not directly participate in the proton transport event, have a significant impact on the dynamics of the excess proton. These secondary participants are shown to strongly influence the active site environment, resulting in the creation of water clusters that are conducive to fast, moderately slow, or slow proton transport events. The identification and characterization of these secondary participants illuminates the role of protein dynamics in the catalytic efficiency of HCA II.

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.

Computationally Efficient Multiscale Reactive Molecular Dynamics to Describe Amino Acid Deprotonation in Proteins

An important challenge in the simulation of biomolecular systems is a quantitative description of the protonation and deprotonation process of amino acid residues. Despite the seeming simplicity of adding or removing a positively charged hydrogen nucleus, simulating the actual protonation/deprotonation process is inherently difficult. It requires both the explicit treatment of the excess proton, including its charge defect delocalization and Grotthuss shuttling through inhomogeneous moieties (water and amino residues), and extensive sampling of coupled condensed phase motions. In a recent paper (J. Chem. Theory Comput. 2014, 10, 2729-2737), a multiscale approach was developed to map high-level quantum mechanics/molecular mechanics (QM/MM) data into a multiscale reactive molecular dynamics (MS-RMD) model in order to describe amino acid deprotonation in bulk water. In this article, we extend the fitting approach (called FitRMD) to create MS-RMD models for ionizable amino acids within proteins. The resulting models are shown to faithfully reproduce the free energy profiles of the reference QM/MM Hamiltonian for PT inside an example protein, the ClC-ec1 H(+)/Cl(-) antiporter. Moreover, we show that the resulting MS-RMD models are computationally efficient enough to then characterize more complex 2-dimensional free energy surfaces due to slow degrees of freedom such as water hydration of internal protein cavities that can be inherently coupled to the excess proton charge translocation. The FitRMD method is thus shown to be an effective way to map ab initio level accuracy into a much more computationally efficient reactive MD method in order to explicitly simulate and quantitatively describe amino acid protonation/deprotonation in proteins.

NMR studies of active-site properties of human carbonic anhydrase II by using (15) N-labeled 4-methylimidazole as a local probe and histidine hydrogen-bond correlations

By using a combination of liquid and solid-state NMR spectroscopy, (15) N-labeled 4-methylimidazole (4-MI) as a local probe of the environment has been studied: 1) in the polar, wet Freon CDF3 /CDF2 Cl down to 130 K, 2) in water at pH 12, and 3) in solid samples of the mutant H64A of human carbonic anhydrase II (HCA II). In the latter, the active-site His64 residue is replaced by alanine; the catalytic activity is, however, rescued by the presence of 4-MI. For the Freon solution, it is demonstrated that addition of water molecules not only catalyzes proton tautomerism but also lifts its quasidegeneracy. The possible hydrogen-bond clusters formed and the mechanism of the tautomerism are discussed. Information about the imidazole hydrogen-bond geometries is obtained by establishing a correlation between published (1) H and (15) N chemical shifts of the imidazole rings of histidines in proteins. This correlation is useful to distinguish histidines embedded in the interior of proteins and those at the surface, embedded in water. Moreover, evidence is obtained that the hydrogen-bond geometries of His64 in the active site of HCA II and of 4-MI in H64A HCA II are similar. Finally, the degeneracy of the rapid tautomerism of the neutral imidazole ring His64 reported by Shimahara et al. (J. Biol. Chem.- 2007, 282, 9646) can be explained with a wet, polar, nonaqueous active-site conformation in the inward conformation, similar to the properties of 4-MI in the Freon solution. The biological implications for the enzyme mechanism are discussed.

Multiscale simulation reveals a multifaceted mechanism of proton permeation through the influenza A M2 proton channel

The influenza A virus M2 channel (AM2) is crucial in the viral life cycle. Despite many previous experimental and computational studies, the mechanism of the activating process in which proton permeation acidifies the virion to release the viral RNA and core proteins is not well understood. Herein the AM2 proton permeation process has been systematically characterized using multiscale computer simulations, including quantum, classical, and reactive molecular dynamics methods. We report, to our knowledge, the first complete free-energy profiles for proton transport through the entire AM2 transmembrane domain at various pH values, including explicit treatment of excess proton charge delocalization and shuttling through the His37 tetrad. The free-energy profiles reveal that the excess proton must overcome a large free-energy barrier to diffuse to the His37 tetrad, where it is stabilized in a deep minimum reflecting the delocalization of the excess charge among the histidines and the cost of shuttling the proton past them. At lower pH values the His37 tetrad has a larger total charge that increases the channel width, hydration, and solvent dynamics, in agreement with recent 2D-IR spectroscopic studies. The proton transport barrier becomes smaller, despite the increased charge repulsion, due to backbone expansion and the more dynamic pore water molecules. The calculated conductances are in quantitative agreement with recent experimental measurements. In addition, the free-energy profiles and conductances for proton transport in several mutants provide insights for explaining our findings and those of previous experimental mutagenesis studies.

Multiscale Reactive Molecular Dynamics for Absolute pKa Predictions and Amino Acid Deprotonation

Accurately calculating a weak acid's pK_a from simulations remains a challenging task. We report a multiscale theoretical approach to calculate the free energy profile for acid ionization, resulting in accurate absolute pK_a values in addition to insights into the underlying mechanism. Importantly, our approach minimizes empiricism by mapping electronic structure data (QM/MM forces) into a reactive molecular dynamics model capable of extensive sampling. Consequently, the bulk property of interest (the absolute pK_a) is the natural consequence of the model, not a parameter used to fit it. This approach is applied to create reactive models of aspartic and glutamic acids. We show that these models predict the correct pK_a values and provide ample statistics to probe the molecular mechanism of dissociation. This analysis shows changes in the solvation structure and Zundel-dominated transitions between the protonated acid, contact ion pair, and bulk solvated excess proton.

Structural insight into activity enhancement and inhibition of H64A carbonic anhydrase II by imidazoles

Human carbonic anhydrases (CAs) are zinc metalloenzymes that catalyze the hydration and dehydration of CO2 and HCO3 (-), respectively. The reaction follows a ping-pong mechanism, in which the rate-limiting step is the transfer of a proton from the zinc-bound solvent (OH(-)/H2O) in/out of the active site via His64, which is widely believed to be the proton-shuttling residue. The decreased catalytic activity (∼20-fold lower with respect to the wild type) of a variant of CA II in which His64 is replaced with Ala (H64A CA II) can be enhanced by exogenous proton donors/acceptors, usually derivatives of imidazoles and pyridines, to almost the wild-type level. X-ray crystal structures of H64A CA II in complex with four imidazole derivatives (imidazole, 1--methylimidazole, 2--methylimidazole and 4-methylimidazole) have been determined and reveal multiple binding sites. Two of these imidazole binding sites have been identified that mimic the positions of the 'in' and 'out' rotamers of His64 in wild-type CA II, while another directly inhibits catalysis by displacing the zinc-bound solvent. The data presented here not only corroborate the importance of the imidazole side chain of His64 in proton transfer during CA catalysis, but also provide a complete structural understanding of the mechanism by which imidazoles enhance (and inhibit when used at higher concentrations) the activity of H64A CA II.

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.

Carbonic anhydrase binding site parameterization in OPLS-AA force field

The parameterization of carbonic anhydrase binding site in OPLS-AA force field was performed using quantum chemistry calculations. Both OH2 and OH(-) forms of the binding site were considered, showing important differences in terms of atomic partial charges. Three different parameterization protocols were used, and the results obtained highlighted the importance of including an extended binding site in the charge calculation. The force field parameters were subsequently validated using standard molecular dynamics simulations. The results presented in this work should greatly facilitate the use of molecular dynamics simulations for studying the carbonic anhydrase, and more generally, the metallo-enzymes.

Modeling the structure and proton transfer pathways of the mutant His-107-Tyr of human carbonic anhydrase II

We present molecular modeling of the structure and possible proton transfer pathways from the surface of the protein to the zinc-bound water molecule in the active site of the mutant His-107-Tyr of human carbonic anhydrase II (HCAII). No high-resolution structure or crystal structure is available till now for this particular mutant due to its lack of stability at physiological temperature. Our analysis utilizes as starting point a series of structures derived from high-resolution crystal structure of the wild type protein. While many of the structures investigated do not reveal a complete path between the zinc bound water and His-64, several others do indicate the presence of a transient connection even when His-64 is present in its outward conformation. Mutation at the residue 107 also reveals the formation of a new path into the active site. Competing contributions from His-64 sidechain rotation from its outward conformation are also evaluated in terms of optimal path analysis. No indication of a lower catalytic efficiency of the mutant is evident from our results under the condition of thermal stability of the mutant.

Structure, function and applications of carbonic anhydrase isozymes

The carbonic anhydrases enzymes (CAs, EC 4.2.1.1) are zinc containing metalloproteins, which efficiently catalyse the reversible conversion of carbon dioxide to bicarbonate and release proton. These enzymes are essentially important for biological system and play several important physiological and patho-physiological functions. There are 16 different alpha-carbonic anhydrase isoforms studied, differing widely in their cellular localization and biophysical properties. The catalytic domains of all CAs possess a conserved tertiary structure fold, with predominately β-strands. We performed an extensive analysis of all 16 mammalian CAs for its structure and function in order to establish a structure-function relationship. CAs have been a potential therapeutic target for many diseases. Sulfonamides are considered as a strong and specific inhibitor of CA, and are being used as diuretics, anti-glaucoma, anti-epileptic, anti-ulcer agents. Currently CA inhibitors are widely used as a drug for the treatment of neurological disorders, anti-glaucoma drugs, anti-cancer, or anti-obesity agents. Here we tried to emphasize how CAs can be used for drug discovery, design and screening. Furthermore, we discussed the role of CA in carbon capture, carbon sensor and metabolon. We hope this review provide many useful information on structure, function, mechanism, and applications of CAs in various discipline.

Expanding the view of proton pumping in cytochrome c oxidase through computer simulation

In cytochrome c oxidase (CcO), a redox-driven proton pump, protons are transported by the Grotthuss shuttling via hydrogen-bonded water molecules and protonatable residues. Proton transport through the D-pathway is a complicated process that is highly sensitive to alterations in the amino acids or the solvation structure in the channel, both of which can inhibit proton pumping and enzymatic activity. Simulations of proton transport in the hydrophobic cavity showed a clear redox state dependence. To study the mechanism of proton pumping in CcO, multi-state empirical valence bond (MS-EVB) simulations have been conducted, focusing on the proton transport through the D-pathway and the hydrophobic cavity next to the binuclear center. The hydration structures, transport pathways, effects of residues, and free energy surfaces of proton transport were revealed in these MS-EVB simulations. The mechanistic insight gained from them is herein reviewed and placed in context for future studies.