Abstract
The active site architecture of heme enzymes has evolved to control the formation of highly reactive intermediates for oxidative catalysis. Finely tuned proton delivery to heme is essential, yet the mechanisms of proton delivery and the sources of protons are poorly understood. Here, we identify routes and drivers of proton delivery in a heme peroxidase enzyme (ascorbate peroxidase) employing a range of computational molecular modeling approaches, including molecular dynamics (MD) simulations, density functional theory (DFT) and combined quantum mechanics/molecular mechanics (QM/MM) calculations, enhanced sampling techniques, and local electric field (LEF) analyses. Our results show that networks of active site water molecules facilitate proton ex-change with Arg38, which may serve as a transient proton carrier at the γ-edge of the heme, where the substrate (ascorbate) binds. The distal His42 residue aids proton transfer into the active site via solvent from the δ-heme edge. MD simulations of three peroxidases identify hydrated channels in the structure – pointing towards both the γ- and δ-edges of the heme – through which protons from solvent can access the active site. Comparison with structures of 8 other heme peroxidases shows that these channels are conserved features. LEF calculations show a continuous electrostatic funnel, running from positive to negative, that pulls protons from the protein surface along the channels towards the heme iron from both the γ- and δ-heme edges. This electric field is also found to be a conserved feature in several peroxidases; its direction and shape is conserved but the gradient varies between differ-ent enzymes, suggesting that nature preorganises the electric field as a determinant of relative activity. Rather than a simplistic model for proton delivery (in which a single substrate provides a single proton (and an electron) to the heme) the data instead show conserved solvent channels and electrostatic funnels, which provide multiple well defined routes for proton delivery in peroxidase catalysis.
Supplementary materials
Title
Supporting Information
Description
Details of DFT, QM(DFT)/MM calculations, MM MD simula-tions, QM(DFTB2)/MM MD simulations and LEF calcula-tions and visualization; Tables S1-S4; Figure S1-S6. Excel file containing Mulliken charges and spin densities for each atom in the DFT/MM-optimized structures, calculated with and without the MM point charges. This dataset highlights the effect of protein polarization on the reaction centre.
Complementary videos for Figure 2 (movie S1-S2) and for the local electric field lines visualization in VR (movie S3-S5). Coordinates of the DFT optimized largest cluster models and QM/MM-optimized structures for different tautomers and spin multiplicities. Custom analysis scripts for identifying consecutive water chains in MD.
Actions
Title
Excel file containing Mulliken charges and spin densities for each atom in the DFT/MM-optimized structures, calculated with and without the MM external point charges
Description
Excel file containing Mulliken charges and spin densities for each atom in the DFT/MM-optimized structures, calculated with and without the MM external point charges. This dataset high-lights the effect of distant protein polarization on the reaction centre.
Actions
Title
Complementary videos for figure 2
Description
Complementary videos for figure 2A
Actions
Title
Complementary videos for figure 2
Description
Complementary videos for figure 2B
Actions
Title
Coordinates of the DFT optimized largest cluster models
Description
Coordinates of the DFT optimized largest cluster models
Actions
Title
QM/MM-optimized structures
Description
QM/MM-optimized structures for different tautomers and spin multiplicities
Actions
Title
Complementary videos for the local electric field lines
Description
Complementary videos for the local electric field lines visualization in VR
Actions
Title
Complementary videos for the local electric field lines
Description
Complementary videos for the local electric field lines visualization in VR
Actions
Title
Complementary videos for the local electric field lines
Description
Complementary videos for the local electric field lines visualization in VR
Actions
Title
Tcl scripts for MD analysis
Description
Tcl scripts for MD analysis
Actions