Coordination Dynamics of Zinc Triggers the Rate Determining Proton Transfer in Human Carbonic Anhydrase II

Controllable Solid Electrolyte Interphase by Ionic Environment Regulation for Stable Zn-Ion Battery

Artificial solid electrolyte interphase in organic solutions is effective and facile for long-cycling aqueous zinc ion batteries. However, the specific effects on different ionic environments have not been thoroughly investigated. Herein, pyromellitic acid (PA) are employed as organic ligand to coordinate with Zn2+ under various ionic environments. The connection between the ionic environment and reaction spontaneity is analyzed to provide insights into the reasons behind the effectiveness of the SEI layer and to characterize its protective impact on the zinc anode. Notably, the PA solution (pH4) lacking OH- contributes to the formation of a dense and ultrathin SEI with Zn-PA coordination, preventing direct contact between the anode and electrolyte. Moreover, the presence of organic functional groups facilitates a uniform flux of Zn2+. These advantages enable stable cycling of the PA4-Zn symmetric cell at a current density of 3 mA cm-2 for over 3500 h. The PA4-Zn//MVO full cell demonstrates excellent electrochemical reversibility. Investigating the influence of the ionic environment on SEI generation informs the development of novel SEI strategies.

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.

Perspective: Path Sampling Methods Applied to Enzymatic Catalysis

This Perspective reviews the use of Transition Path Sampling methods to study enzymatically catalyzed chemical reactions. First applied by our group to an enzymatic reaction over 15 years ago, the method has uncovered basic principles in enzymatic catalysis such as the protein promoting vibration, and it has also helped harmonize such ideas as electrostatic preorganization with dynamic views of enzyme function. It is now being used to help uncover principles of protein design necessary to artificial enzyme creation.

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.

In the footsteps of an inhibitor unbinding from the active site of human carbonic anhydrase II

The crystal structure of human carbonic anhydrase (HCA) II bound to an inhibitor molecule, 6-hydroxy-2-thioxocoumarin (FC5), shows FC5 to be located in a hydrophobic pocket at the active site. The present work employs classical molecular dynamics (MD) simulation to follow the FC5 molecule for 1 μs as it unbinds from its binding location, adopts the path of substrate/product diffusion (path 1) to leave the active site at around 75 ns. It is then found to undergo repeated binding and unbinding at different locations on the surface of the enzyme in water. Several transient excursions through different regions of the enzyme are also observed prior to its exit from the active site. These transient paths are combined with functionally relevant cavities/channels to enlist five additional pathways (path 2-6). Pathways 1-6 are subsequently explored using steered MD and umbrella sampling simulations. A free energy barrier of 0.969 kcal mol^-1 is encountered along path 1, while barriers in the range of 0.57-2.84 kcal mol^-1 are obtained along paths 2, 4 and 5. We also analyze in detail the interaction between FC5 and the enzyme along each path as the former leaves the active site of HCA II. Our results indicate path 1 to be the major exit pathway for FC5, although competing contributions may also come from the paths 2, 4 and 5.Communicated by Ramaswamy H. Sarma.

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.

Free Energy Landscape and Proton Transfer Pathways of the Transimination Reaction at the Active site of the Serine Hydroxymethyltransferase Enzyme in Aqueous Medium

Serine hydroxymethyltransferase (SHMT) is a ubiquitous enzyme belonging to the fold type I or aspartate aminotransferase (AspAT) family of the pyridoxal 5'-phosphate (PLP)-dependent enzymes. Like other PLP-dependent enzymes, SHMT also undergoes the so-called transimination reaction before exhibiting its enzymatic activity. The transimination process constitutes an important pre-step for all PLP-dependent enzymes, where an internal aldimine of the PLP-enzyme complex gets converted to an external aldimine of the substrate-PLP complex at the active site of the enzyme. In case of the transimination reaction involving SHMT, the PLP molecule bound to the active site lysine residue of SHMT (internal aldimine) gets detached from the enzyme by a serine substrate to produce an external aldimine complex, where the PLP is now bound to the serine substrate. In the current study, the free energy surfaces and reaction pathways of different steps of the transimination reaction at the active site of SHMT are investigated by employing hybrid quantum mechanical/molecular mechanical (QM/MM) simulations combined with metadynamics methods of rare event sampling. It is found that the process of transimination involving serine and PLP at the active site of the SHMT enzyme takes place through different elementary steps such as the formation of the first geminal diamine intermediate (GDI1), transfer of a proton from the substrate serine to the phenolic oxygen of PLP, followed by another proton transfer from PLP to the amine nitrogen of lysine with the formation of the second geminal diamine intermediate (GDI2), and finally, detachment of the active site lysine residue from PLP to produce the external aldimine.

Unexpected Deprotonation from a Chemically Inert OH Group Promoted by Metal Ions in Lanthanide-Erythritol Complexes

Single-crystal structures of five lanthanide-erythritol complexes are reported. The analysis of the chemical compositions and scrutinization of structural features in the single-crystal data of the complexes led us to find that unexpected deprotonation occurs on the OH group of erythritol of three complexes. Considering these complexes were prepared in acidic environments, where spontaneous ionization on an OH group is suppressed, we suggest metal ions play an important role in promoting the proton transfer. To find out why the chemically inert OH is activated, the single-crystal structures of 63 rare-earth complexes containing organic ligands with multiple hydroxyl groups (OLMHs) were surveyed. The formation of μ₂-bridges turns out to be directly relevant to the occurrence of deprotonation. When an OH group from an OLMH molecule participates in the formation of a μ₂-bridge, the polarization ability of the metal ions becomes strong enough to promote the deprotonation on the OH group. The above structural characteristics may be useful in the rational design of catalysts that can activate the chemically inert OH group and promote the relevant chemical conversions.