Overview: Nucleation of clathrate hydrates

Molecular level knowledge of nucleation and growth of clathrate hydrates is of importance for advancing fundamental understanding on the nature of water and hydrophobic hydrate formers, and their interactions that result in the formation of ice-like solids at temperatures higher than the ice-point. The stochastic nature and the inability to probe the small length and time scales associated with the nucleation process make it very difficult to experimentally determine the molecular level changes that lead to the nucleation event. Conversely, for this reason, there have been increasing efforts to obtain this information using molecular simulations. Accurate knowledge of how and when hydrate structures nucleate will be tremendously beneficial for the development of sustainable hydrate management strategies in oil and gas flowlines, as well as for their application in energy storage and recovery, gas separation, carbon sequestration, seawater desalination, and refrigeration. This article reviews various aspects of hydrate nucleation. First, properties of supercooled water and ice nucleation are reviewed briefly due to their apparent similarity to hydrates. Hydrate nucleation is then reviewed starting from macroscopic observations as obtained from experiments in laboratories and operations in industries, followed by various hydrate nucleation hypotheses and hydrate nucleation driving force calculations based on the classical nucleation theory. Finally, molecular simulations on hydrate nucleation are discussed in detail followed by potential future research directions.

Correction: Unravelling the impact of hydrocarbon structure on the fumarate addition mechanism – a gas-phase ab initio study

Correction for ‘Unravelling the impact of hydrocarbon structure on the fumarate addition mechanism – a gas-phase ab initio study’ by Vivek S. Bharadwaj et al., Phys. Chem. Chem. Phys., 2015, 17, 4054–4066.

Quantum-Mechanical Study of the Reaction Mechanism for 2π-2π Cycloaddition of Fluorinated Methylene Groups

Perfluorocyclobutyl polymers are thermally and chemically stable, may be produced without a catalyst via thermal 2π-2π cycloaddition, and can form block structures, making them suitable for commercialization of specialty polymers. Thermal 2π-2π cycloaddition is a rare reaction that begins in the singlet state and proceeds through a triplet intermediate to form an energetically stable four-membered ring in the singlet state. This reaction involves two changes in spin state and, thus, two spin-crossover transitions. Presented here are density functional theory calculations that evaluate the energetics and reaction mechanisms for the dimerizations of two different polyfluorinated precursors, 1,1,2-trifluoro-2-(trifluoromethoxy)ethane and hexafluoropropylene. The spin-crossover transition states are thoroughly investigated, revealing important kinetics steps and an activation energy for the gas-phase cycloaddition of two hexafluoropropene molecules of 36.9 kcal/mol, which is in good agreement with the experimentally determined value of 34.3 kcal/mol. It is found that the first carbon-carbon bond formation is the rate-limiting step, followed by a rotation about the newly formed bond in the triplet state that results in the formation of the second carbon-carbon bond. Targeting the rotation of the C-C bond, a set of parameters were obtained that best produce high molecular weight polymers using this chemistry.

A Novel Optical Diagnostic for In Situ Measurements of Lithium Polysulfides in Battery Electrolytes

An optical diagnostic technique to determine the order and concentration of lithium polysulfides in lithium-sulfur (Li-S) battery electrolytes has been developed. One of the major challenges of lithium-sulfur batteries is the problem of polysulfide shuttling between the electrodes, which leads to self-discharge and loss of active material. Here we present an optical diagnostic for quantitative in situ measurements of lithium polysulfides using attenuated total reflection Fourier transform infrared (FT-IR) spectroscopy. Simulated infrared spectra of lithium polysulfide molecules were generated using computational quantum chemistry routines implemented in Gaussian 09. The theoretical spectra served as a starting point for experimental characterization of lithium polysulfide solutions synthesized by the direct reaction of lithium sulfide and sulfur. Attenuated total reflection FT-IR spectroscopy was used to measure absorption spectra. The lower limit of detection with this technique is 0.05 M. Measured spectra revealed trends with respect to polysulfide order and concentration, consistent with theoretical predictions, which were used to develop a set of equations relating the order and concentration of lithium polysulfides in a sample to the position and area of a characteristic infrared absorption band. The diagnostic routine can measure the order and concentration to within 5% and 0.1 M, respectively.

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.

In silico insights into the solvation characteristics of the ionic liquid 1-methyltriethoxy-3-ethylimidazolium acetate for cellulosic biomass

Lignocellulosic biomass is a domestically grown, sustainable, and potentially carbon-neutral feedstock for the production of liquid fuels and other value added chemicals. This underutilized renewable feedstock has the potential to alleviate some of the current socio-economic dependence on foreign petroleum supplies while stimulating rural economies. Unfortunately, the potential of biomass has largely been underdeveloped due to the recalcitrant nature of lignocellulosic materials. Task-specific ionic liquids (ILs) have shown considerable promise as an alternative non-aqueous solvent for solvation and deconstruction of lignocellulose in the presence of metal chloride catalyst or enzymes. Recently it has been hypothesized that adding oxygen atoms to the tail of an imidazolium cation would alleviate some of the negative characteristics of the ILs by increasing mass transport properties, and decreasing IL deactivation of enzymes, while at the same time retaining favorable solvation characteristics for lignocellulose. Reported here are fully atomistic molecular dynamic simulations of 1-methyltriethoxy-3-ethylimidazolium acetate ([Me-(OEt)3-Et-IM(+)] [OAc(-)]) that elucidate promising molecular-level details pertaining to the solvation characteristics of model compounds of cellulose, and IL-induced side-chain and ring puckering conformations. It is found that the anion interactions with the saccharide induce alternate ring puckering conformations from those seen in aqueous environments (i.e.(1)C4), while the cation interactions are found to influence the conformation of the ω dihedral. These perturbations in saccharide structures are discussed in the context of their contribution to the disruption of hydrogen bonding in cellulosic architecture and their role in solvation.

Dynamic free energy surfaces for sodium diffusion in type II silicon clathrates

Earth abundant semiconducting type II Si clathrates have attracted attention as photovoltaic materials due to their wide band gaps. To realize the semiconducting properties of these materials, guest species that arise during the synthesis process must be completely evacuated from the host cage structure post synthesis. A common guest species utilized in the synthesis of Si clathrates is Na (metal), which templates the clathrate cage formation. Previous experimental investigations have identified that it is possible to evacuate Na from type II clathrates to an occupancy of less than 1 Na per unit cell. This work investigates the energetics, kinetics, and resulting mechanism of Na diffusion through type II Si clathrates by means of biased molecular dynamics and kinetic Monte Carlo simulations. Well-tempered metadynamics has been used to determine the potential of mean force for Na moving between clathrate cages, from which the thermodynamic preferences and transition barrier heights have been obtained. Kinetic Monte Carlo simulations based on the metadynamics results have identified the mechanism of Na diffusion in type II Si clathrates. The overall mechanism consists of a coupled diffusive process linked via electrostatic guest-guest interactions. The large occupied hexakaidechedral cages initially empty their Na guests to adjacent empty large cages, thereby changing the local electrostatic environment around the occupied small pentagonal dodecahedral cages and increasing the probability of Na guests to leave the small cages. This coupled process continues through the cross-over point that is identified as the point where large and small cages are equally occupied by Na guests. Further Na removal results in the majority of guests residing in the large cages as opposed to the small cages, in agreement with experiments, and ultimately a Na free structure.

Unravelling the impact of hydrocarbon structure on the fumarate addition mechanism--a gas-phase ab initio study

The fumarate addition reaction mechanism is central to the anaerobic biodegradation pathway of various hydrocarbons, both aromatic (e.g., toluene, ethyl benzene) and aliphatic (e.g., n-hexane, dodecane). Succinate synthase enzymes, which belong to the glycyl radical enzyme family, are the main facilitators of these biochemical reactions. The overall catalytic mechanism that converts hydrocarbons to a succinate molecule involves three steps: (1) initial H-abstraction from the hydrocarbon by the radical enzyme, (2) addition of the resulting hydrocarbon radical to fumarate, and (3) hydrogen abstraction by the addition product to regenerate the radical enzyme. Since the biodegradation of hydrocarbon fuels via the fumarate addition mechanism is linked to bio-corrosion, an improved understanding of this reaction is imperative to our efforts of predicting the susceptibility of proposed alternative fuels to biodegradation. An improved understanding of the fuel biodegradation process also has the potential to benefit bioremediation. In this study, we consider model aromatic (toluene) and aliphatic (butane) compounds to evaluate the impact of hydrocarbon structure on the energetics and kinetics of the fumarate addition mechanism by means of high level ab initio gas-phase calculations. We predict that the rate of toluene degradation is ∼100 times faster than butane at 298 K, and that the first abstraction step is kinetically significant for both hydrocarbons, which is consistent with deuterium isotope effect studies on toluene degradation. The detailed computations also show that the predicted stereo-chemical preference of the succinate products for both toluene and butane are due to the differences in the radical addition rate constants for the various isomers. The computational and kinetic modeling work presented here demonstrates the importance of considering pre-reaction and product complexes in order to accurately treat gas phase systems that involve intra and inter-molecular non-covalent interactions.

The impact of active site protonation on substrate ring conformation in Melanocarpus albomyces cellobiohydrolase Cel7B

Understanding how the protonation state of active site residues impacts the enzyme's structure and substrate conformation is important for improving the efficiency and economic viability of the degradation of cellulosic materials as feedstock for liquid fuel and value-added chemicals.

Elucidating the conformational energetics of glucose and cellobiose in ionic liquids

The energetics driving the conformational preference of the ω dihedral of glucose and the φ–ψ dihedrals of cellobiose solvated in imidazolium acetate ionic liquids and water are elucidated and compared.

Homology modeling and molecular docking for the science curriculum

DockoMatic 2.0 is a powerful open source software program (downloadable from sourceforge.net) that allows users to utilize a readily accessible computational tool to explore biomolecules and their interactions. This manuscript describes a practical tutorial for use in the undergraduate curriculum that introduces students to macromolecular structure creation, ligand binding calculations, and visualization of docking results. A student procedure is provided that illustrates the use of DockoMatic to create a homology model for the amino propeptide region (223 amino acids with two disulfide bonds) of collagen α1 (XI), followed by molecular docking of the commercial drug Arixtra(®) to the homology model of α1 (XI), and finally, analysis of the results of the docking experiment. The activities and Supporting Information described are intended to educate students in the use of computational tools to create and investigate homology models for other systems of interest and to train students to perform and analyze molecular docking studies. The tutorial also serves as a foundation for investigators seeking to explore the viability of using computational biochemistry to study their receptor-ligand binding motifs. © 2013 by The International Union of Biochemistry and Molecular Biology, 42(2):179-182, 2014.

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.

Acetylcholine promotes binding of α-conotoxin MII at α3 β2 nicotinic acetylcholine receptors

α-Conotoxin MII (α-CTxMII) is a 16-residue peptide with the sequence GCCSNPVCHLEHSNLC, containing Cys2-Cys8 and Cys3-Cys16 disulfide bonds. This peptide, isolated from the venom of the marine cone snail Conus magus, is a potent and selective antagonist of neuronal nicotinic acetylcholine receptors (nAChRs). To evaluate the impact of channel-ligand interactions on ligand-binding affinity, homology models of the heteropentameric α3β2-nAChR were constructed. The models were created in MODELLER with the aid of experimentally characterized structures of the Torpedo marmorata-nAChR (Tm-nAChR, PDB ID: 2BG9) and the Aplysia californica-acetylcholine binding protein (Ac-AChBP, PDB ID: 2BR8) as templates for the α3- and β2-subunit isoforms derived from rat neuronal nAChR primary amino acid sequences. Molecular docking calculations were performed with AutoDock to evaluate interactions of the heteropentameric nAChR homology models with the ligands acetylcholine (ACh) and α-CTxMII. The nAChR homology models described here bind ACh with binding energies commensurate with those of previously reported systems, and identify critical interactions that facilitate both ACh and α-CTxMII ligand binding. The docking calculations revealed an increased binding affinity of the α3β2-nAChR for α-CTxMII with ACh bound to the receptor, and this was confirmed through two-electrode voltage clamp experiments on oocytes from Xenopus laevis. These findings provide insights into the inhibition and mechanism of electrostatically driven antagonist properties of the α-CTxMIIs on nAChRs.

pKa determination of histidine residues in α-conotoxin MII peptides by 1H NMR and constant pH molecular dynamics simulation

α-Conotoxin MII (α-CTxMII) is a potent and selective peptide antagonist of neuronal nicotinic acetylcholine receptors (nAChR's). Studies have shown that His9 and His12 are significant determinants of toxin binding affinity for nAChR, while Glu11 may dictate differential toxin affinity between nAChR isoforms. The protonation state of these histidine residues and therefore the charge on the α-CTx may contribute to the observed differences in binding affinity and selectivity. In this study, we assess the pH dependence of the protonation state of His9 and His12 by (1)H NMR spectroscopy and constant pH molecular dynamics (CpHMD) in α-CTxMII, α-CTxMII[E11A], and the triple mutant, α-CTxMII[N5R:E11A:H12K]. The E11A mutation does not significantly perturb the pKa of His9 or His12, while N5R:E11A:H12K results in a significant decrease in the pKa value of His9. The pKa values predicted by CpHMD simulations are in good agreement with (1)H NMR spectroscopy, with a mean absolute deviation from experiment of 0.3 pKa units. These results support the use of CpHMD as an efficient and inexpensive predictive tool to determine pKa values and structural features of small peptides critical to their function.

Chemical rescue of enzymes: proton transfer in mutants of human carbonic anhydrase II

In human carbonic anhydrase II (HCA II), the mutation of position 64 from histidine to alanine (H64A) disrupts the rate limiting proton transfer (PT) event, resulting in a reduction of the catalytic activity of the enzyme as compared to the wild-type. Potential of mean force (PMF) calculations utilizing the multistate empirical valence bond (MS-EVB) methodology for H64A HCA II yields a PT free energy barrier significantly higher than that found in the wild-type enzyme. This high barrier, determined in the absence of exogenous buffer and assuming no additional ionizable residues in the PT pathway, indicates the likelihood of alternate enzyme pathways that utilize either ionizable enzyme residues (self-rescue) and/or exogenous buffers (chemical rescue). It has been shown experimentally that the catalytic activity of H64A HCA II can be chemically rescued to near wild-type levels by the addition of the exogenous buffer 4-methylimidazole (4MI). Crystallographic studies have identified two 4MI binding sites, yet site-specific mutations intended to disrupt 4MI binding have demonstrated these sites to be nonproductive. In the present work, MS-EVB simulations show that binding of 4MI near Thr199 in the H64A HCA II mutant, a binding site determined by NMR spectroscopy, results in a viable chemical rescue pathway. Additional viable rescue pathways are also identified where 4MI acts as a proton transport intermediary from the active site to ionizable residues on the rim of the active site, revealing a probable mode of action for the chemical rescue pathway.

The self-consistent charge density functional tight binding method applied to liquid water and the hydrated excess proton: benchmark simulations

The self-consistent charge density functional tight binding (SCC-DFTB) method is a relatively new approximate electronic structure method that is increasingly used to study biologically relevant systems in aqueous environments. There have been several gas phase cluster calculations that indicate, in some instances, an ability to predict geometries, energies, and vibrational frequencies in reasonable agreement with high level ab initio calculations. However, to date, there has been little validation of the method for bulk water properties, and no validation for the properties of the hydrated excess proton in water. Presented here is a detailed SCC-DFTB analysis of the latter two systems. This work focuses on the ability of the original SCC-DFTB method, and a modified version that includes a hydrogen bonding damping function (HBD-SCC-DFTB), to describe the structural, energetic, and dynamical nature of these aqueous systems. The SCC-DFTB and HBD-SCC-DFTB results are compared to experimental data and Car-Parrinello molecular dynamics (CPMD) simulations using the HCTH/120 gradient-corrected exchange-correlation energy functional. All simulations for these systems contained 128 water molecules, plus one additional proton in the case of the excess proton system, and were carried out in a periodic simulation box with Ewald long-range electrostatics. The liquid water structure for the original SCC-DFTB is shown to poorly reproduce bulk water properties, while the HBD-SCC-DFTB somewhat more closely represents bulk water due to an improved ability to describe hydrogen bonding energies. Both SCC-DFTB methods are found to underestimate the water dimer interaction energy, resulting in a low heat of vaporization and a significantly elevated water oxygen diffusion coefficient as compared to experiment. The addition of an excess hydrated proton to the bulk water resulted in the Zundel cation (H(5)O(2)(+)) stabilized species being the stable form of the charge defect, which diffuses at a rate similar to the underlying water diffusion. These SCC-DFTB results differ significantly from other existing computational descriptions of the hydrated excess proton in water, as well as from the available experimental data.

Effect of active-site mutation at Asn67 on the proton transfer mechanism of human carbonic anhydrase II

The rate-limiting proton transfer (PT) event in the site-specific mutant N67L of human carbonic anhydrase II (HCA II) has been examined by kinetic, X-ray, and simulation approaches. The X-ray crystallography studies, which were previously reported, and molecular dynamics (MD) simulations indicate that the proton shuttling residue, His64, predominantly resides in the outward orientation with a significant disruption of the ordered water in the active site for the dehydration pathway. While disorder is seen in the active-site water, water cluster analysis indicates that the N67L mutant may form water clusters similar to those seen in the wild-type (WT). For the hydration pathway of the enzyme, the active site water cluster analysis reveals an inability of the N67L mutant to stabilize water clusters when His64 is in the inward orientation, thereby favoring PT when His64 is in the outward orientation. The preference of the N67L mutant to carry out the PT when His64 is in the outward orientation for both the hydration and dehydration pathway is reasoned to be the main cause of the observed reduction in the overall rate. To probe the mechanism of PT, solvent H/D kinetic isotope effects (KIEs) were experimentally studied with catalysis measured by the exchange of (18)O between CO(2) and water. The values obtained from the KIEs were determined as a function of the deuterium content of solvent, using the proton inventory method. No differences were detected in the overarching mechanism of PT between WT and N67L HCA II, despite changes in the active-site water structure and/or the orientation of His64.

Elucidation of the proton transport mechanism in human carbonic anhydrase II

Human carbonic anhydrase II (HCA II) is one of the fastest known enzymes, which utilizes a rate-limiting proton transport (PT) step in its enzymatic reaction. To evaluate the PT event at an atomistic level, the multistate empirical valence bond (MS-EVB) method has been utilized in this work. It is observed that the PT event in HCA II exploits a transient active site water cluster to transport the excess proton between the catalytic zinc-bound water/hydroxide and the proton shuttling residue, His64. This PT event is found to be dependent on the enzyme's ability to form and stabilize the active site water cluster in addition to its ability to orient His64 in a favorable conformation. Evaluation of the PT free energy barrier for different orientations of His64 reveals this residue's vital role as a proton transporter and elucidates its direct effect on the barrier to PT through the active site water. It is suggested that the rate-limiting step oscillates between the active site water PT event to His64 and the de/protonation of His64 depending on the exogenous buffer concentration and the orientation of His64. In the absence of a PT acceptor/donor at position 64, it is found that the excess proton will utilize one of three distinct paths to enter/leave the active site. This latter result not only allows for an increased understanding of how enzymes capitalize on the protein/solvent interface to guide excess protons to/from areas of interest, it also provides valuable insight into the chemical rescue experiments on HCA II mutants.

Origins of enhanced proton transport in the Y7F mutant of human carbonic anhydrase II

Human carbonic anhydrase II (HCA II), among the fastest enzymes known, catalyzes the reversible hydration of CO 2 to HCO 3 (-). The rate-limiting step of this reaction is believed to be the formation of an intramolecular water wire and transfer of a proton across the active site cavity from a zinc-bound solvent to a proton shuttling residue (His64). X-ray crystallographic studies have shown this intramolecular water wire to be directly stabilized through hydrogen bonds via a small well-defined set of amino acids, namely, Tyr7, Asn62, Asn67, Thr199, and Thr200. Furthermore, X-ray crystallographic and kinetic studies have shown that the mutation of tyrosine 7 to phenylalanine, Y7F HCA II, has the effect of increasing the proton transfer rate by 7-fold in the dehydration direction of the enzyme reaction compared to wild-type (WT). This increase in the proton transfer rate is postulated to be linked to the formation of a more directional, less branched, water wire. To evaluate this proposal, molecular dynamics simulations have been employed to study water wire formation in both the WT and Y7F HCA II mutant. These studies reveal that the Y7F mutant enhances the probability of forming small water wires and significantly extends the water wire lifetime, which may account for the elevated proton transfer seen in the Y7F mutant. Correlation analysis of the enzyme and intramolecular water wire indicates that the Y7F mutant significantly alters the interaction of the active site waters with the enzyme while occupancy data of the water oxygens reveals that the Y7F mutant stabilizes the intramolecular water wire in a manner that maximizes smaller water wire formation. This increase in the number of smaller water wires is likely to elevate the catalytic turnover of an already very efficient enzyme.

A multistate empirical valence bond description of protonatable amino acids

The multistate empirical valence bond (MS-EVB) model, which was developed for molecular dynamics simulations of proton transport in water and biomolecular systems, is extended for the modeling of protonatable amino acid residues in aqueous environments, specifically histidine and glutamic acid. The parameters of the MS-EVB force field are first determined to reproduce the geometries and energetics of the gas phase amino acid-water clusters. These parameters are then optimized to reproduce experimental pK(a) values. The free energy profiles for acid ionization and the corresponding pK(a) values are calculated by MS-EVB molecular dynamics simulations utilizing the umbrella sampling technique, with the center of excess charge coordinate chosen as the dissociation reaction coordinate. A general procedure for fitting the MS-EVB parameters is formulated, which allows for the parametrization of other amino acid residues with protonatable groups and the subsequent use of the MS-EVB approach for molecular dynamics simulations of proton transfer processes in proteins involving protonation/deprotonation of the protonatable amino acid groups.

Proton solvation and transport in aqueous and biomolecular systems: insights from computer simulations

The excess proton in aqueous media plays a pivotal role in many fundamental chemical (e.g., acid-base chemistry) and biological (e.g., bioenergetics and enzyme catalysis) processes. Understanding the hydrated proton is, therefore, crucial for chemistry, biology, and materials sciences. Although well studied for over 200 years, excess proton solvation and transport remains to this day mysterious, surprising, and perhaps even misunderstood. In this feature article, various efforts to address this problem through computer modeling and simulation will be described. Applications of computer simulations to a number of important and interesting systems will be presented, highlighting the roles of charge delocalization and Grotthuss shuttling, a phenomenon unique in many ways to the excess proton in water.

Preferred orientations of His64 in human carbonic anhydrase II

Histidine at position 64 (His64) in human carbonic anhydrase II (HCA II) is believed to be the proton acceptor in the hydration direction and the proton donor in the dehydration direction for the rate-limiting proton transfer (PT) event. Although the biochemical effect of histidine at position 64 has been thoroughly investigated, the role of its orientation in the PT event is a topic of considerable debate. X-ray data of HCA II suggests that His64 can adopt either an "in" or "out" orientation. The "in" orientation is believed to be favored for the hydration direction PT event because the Ndelta of His64 is closer to the catalytic zinc. This orientation allows for smaller water bridges, which are postulated to be more conducive to PT. In the present work, classical molecular dynamics simulations have been conducted to elucidate the role that the His64 orientation may play in its ability to act as a proton donor/acceptor in HCA II. The free energy profile for the orientation of His64 suggests that the histidine will adopt an "in" orientation in the hydration direction, which brings Ndelta in close proximity to the catalytic zinc. When the histidine becomes protonated, it then rotates to an "out" orientation, creating a more favorable solvation environment for the protonated His64. In this "out" orientation, the imidazole ring releases the delta nitrogen's excess proton into the bulk environment. After the second PT event and when the zinc-bound water is regenerated, the His64 is again favored to reorient to the "in" orientation, completing the catalytic cycle.

Atomic Crystal and Molecular Dynamics Simulation Structures of Human Carbonic Anhydrase II: Insights into the Proton Transfer Mechanism†,‡

Human carbonic anhydrase II (HCA II) is a zinc-metalloenzyme that catalyzes the reversible interconversion of CO2 and HCO3-. The rate-limiting step of this catalysis is the transfer of a proton between the Zn-bound solvent molecule and residue His64. In order to fully characterize the active site structural features implicated in the proton transfer mechanism, the refined X-ray crystal structure of uncomplexed wild type HCA II to 1.05 A resolution with an Rcryst value of 12.0% and an Rfree value of 15.1% has been elucidated. This structure provides strong clues as to the pathway of the intramolecular proton transfer between the Zn-bound solvent and His64. The structure emphasizes the role of the solvent network, the unique positioning of solvent molecule W2, and the significance of the dual conformation of His64 in the active site. The structure is compared with molecular dynamics (MD) simulation calculations of the Zn-bound hydroxyl/His64+ (charged) and the Zn-bound water/His64 (uncharged) HCA II states. A comparison of the crystallographic anisotropic atomic thermal parameters and MD simulation root-mean-square fluctuation values show excellent agreement in the atomic motion observed between the two methods. It is also interesting that the observed active site solvent positions in the crystal structure are also the most probable positions of the solvent during the MD simulations. On the basis of the comparative study of the MD simulation results, the HCA II crystal structure observed is most likely in the Zn-bound water/His64 state. This conclusion is based on the following observations: His64 is mainly (80%) orientated in an inward conformation; electron density omit maps infer that His64 is not charged in an either inward or outward conformation; and the Zn-bound solvent is most likely a water molecule.