Abstract
The recently discovered metal-free carbonic anhydrase (CA) enzyme may have a significant impact on the global carbon dioxide (CO2) cycle, as it can irreversibly perform the CO2 hydration reaction. In this study, we investigate several key aspects of metal-free CA, including the identification of the catalytic site, determination of the CO2 binding site, and the mechanism of catalysis. This is achieved through the use of classical molecular dynamics (MD) simulations, quantum chemical density functional theory (DFT), and hybrid quantum mechanics/molecular mechanics (QM/MM) calculations. Our study indicates that the experimental structure based on X-ray crystallography, which shows the ‘bicarbonate (HCO3−) product’ trapped in the hydrophilic region of the metal-free CA, might not accurately depict the actual enzyme-substrate picture. Instead, the simulation reveals that the CO2 prefers the hydrophobic zone, which serves as the primary catalytic site. It also highlights the strategic role of a gatekeeper residue (Phe504), which assists in regulating the transportation of CO2 by tilting its aromatic plane. Additionally, the hybrid QM/MM calculations establish that CO2 hydration is catalyzed within the hydrophobic zone by a deprotonated tyrosine with the help of an organized water chain.
Supplementary materials
Title
Supplementary information contains optimised coordinates and other relevant information related to work.
Description
We have provided following information:
Adaptive Molecular Dynamics Simulation Methodology
Representative MD snapshots and corresponding RMSD plots
Snapshots from the simulation of CO 2 in the hydrophobic zone
Water Assisted CO 2 Hydration
Scheme for water assisted CO 2 hydration and corresponding reaction profile
Figure Proposed schematic pathways for CO 2 hydration
Occupancy plots of potential proton acceptor amino acids
Radial distribution plot for amino acid around 5 Å from CO 2
Search of well-known proton acceptor residues
QM/MM optimized snapshots and distance plots
QM/MM explored CO 2 hydration with His584 QM/MM explored CO 2 hydration with Glu498 Radial distribution function between Tyr476 and water
Water population plot
QM/MM optimized MD snapshots and corresponding energy profile diagram
Calculation
QM coordinates of QM/MM Optimized geometries
Actions