- Research article
- Open Access
A role of proton transfer in peroxidase-catalyzed process elucidated by substrates docking calculations
© Kulys and Ziemys; licensee BioMed Central Ltd. 2001
Received: 5 July 2001
Accepted: 28 August 2001
Published: 28 August 2001
Previous kinetic investigations of fungal-peroxidase catalyzed oxidation of N-aryl hydroxamic acids (AHAs) and N-aryl-N- hydroxy urethanes (AHUs) revealed that the rate of reaction was independent of the formal redox potential of substrates. Moreover, the oxidation rate was 3–5 orders of magnitude less than for oxidation of physiological phenol substrates, though the redox potential was similar.
To explain the unexpectedly low reactivity of AHAs and AHUs we made ab initio calculations of the molecular structure of the substrates following in silico docking in the active center of the enzyme.
AHAs and AHUs were docked at the distal side of heme in the sites formed by hydrophobic amino acid residues that retarded a proton transfer and finally the oxidation rate. The analogous phenol substrates were docked at different sites permitting fast proton transfer in the relay of distal His and water that helped fast substrate oxidation.
Heme peroxidases are widespread in bacteria, fungi, plants and mammalians . Therefore, peroxidase-catalyzed oxidation of organic compounds is a paramount electron transfer process in molecular biology. Phenols and anilines are generally recognized as substrates of the heme peroxidases (donor: H2O2 oxidoreductases EC 1.11.17). The peroxidases catalyze oxidation of the substrates by hydrogen peroxide or alkyl peroxides, usually but not always, via free-radical intermediates [1, 2]. Non-phenolic compounds, such as indole-3-acetic acid, phenylenediamines, ferrocenes, phenothiazines, phenoxazines, have also been investigated as peroxidase substrates [3–5]. Steady-state kinetics of peroxidase action has been described as a ping-pong scheme with compound I and compound II formation .
The principal question in enzyme-catalyzed processes is the dependence of the reaction rate on the substrate and the enzyme active center structure. In peroxidase-catalyzed reactions it has been demonstrated that the reaction rate correlated with the substitution parameters of phenols . The reactivity of the horseradish peroxidase toward phenols and non-phenolic substrates, i.e. indole-3-acetic acids, was compared on the basis of the thermodynamic driving force of the reaction . These observations as well as kinetics of phenylenediamines, phenothiazines and phenoxazines oxidation revealed that reactivity of peroxidases (at least within the same substrate type) were depend on redox properties of the compounds [3–6].
The peroxidase substrates investigated in [2, 5] contained phenol or aniline structure in which -OH or -NH groups were conjugated with aromatic system. The oxidation rate of these substrates approached diffusion limit. However, during investigations of N-arylhydroxamic acids and N-aryl-N-hydroxyurethanes, containing redox core Ar-N(OH)-COR and Ar-N(OH)-CO-OR, unexpectedly low reaction rate has been established [7, 8]. In order to explain the sharp difference of reactivity of phenols and non phenolic N-OH substrates in this study peroxidase kinetic results were compared with molecular parameters of substrates and their binding in the enzyme active centre.
pH dependence of redox potential has not been associated with acid-base dissociation of AHA 1a since absorption spectra didn't change at these pH values. Only nitro compounds, i.e.AHA 1b andAHA 2c changed absorption spectrum at different pH with an apparent pKa of transitions 8.50 ± 0.1 and 8.0 ± 0.4, respectively . The absorbance change was associated with heterolytic NO-H bond dissociation at alkaline pH.
The dependence of redox potential on pH could be explained by different oxidation mechanism. At pH larger than 7.5 the electrochemical conversion was associated with single electron transfer followed by fast proton-release (equations 1,2). At these pH values the process was not limited by proton transfer.
CH3-CO-N(OH)-Ph → CH3-CO-N(O•+H)-Ph + e (1)
CH3-CO-N(O•+H)-Ph → CH3-CO-NO•-Ph + H+ (2)
In pH interval between pH 3.7–6.3 the oxidation proceeded by electron/proton transfer. At acidic pH (pH<3.7) AHA 1a oxidation was linked to electron transfer following the radical cation formation. The high redox potential of this conversion has been limited by slow proton transfer. The formal potentials of analyzed conversions were 368 ± 2, 481 ± 9 and 671 ± 170 mV (vs SCE), respectively. The role for proton transfer was confirmed during AHA 1a oxidation in methylene chloride. The oxidation of AHA 1a in this aprotic solvent proceeded at 1.0 –1.2 V.
370 ± 2
(7.1 ± 0.2) · 103
375 ± 6
(9.9 ± 0.3) · 104
411 ± 2
(2.5 ± 0.1) · 105
313 ± 2
(1.5 ± 0.1) · 105
406 ± 1
(8.5 ± 0.6) · 103
375 ± 4
(7.1 ± 0.4) · 104
393 ± 2
(9.2 ± 2.1) · 104
391 ± 3
(5.0 ± 0.3) · 104
369 ± 1
(1.5 ± 0.2) · 107
364 ± 2
(8.3 ± 0.6) · 103
403 ± 1
(9.2 ± 0.6) · 103
307 ± 2
(8.7 ± 0.4) · 105
345 ± 1
(2.3 ± 0.2) · 104
337 ± 2
(7.5 ± 0.5) · 104
328 ± 2
(3.2 ± 0.2) · 104
351 ± 4
(4.2 ± 0.2) · 104
396 ± 5
(2.5 ± 0.5) · 104
The calculated redox potentials of the substrates AHU 1-4 changed in the range between 328 and 351 mV (Table 1). Formal potential of AHU 1 was also dependent on pH (Figure 3). pKa of transitions were 5.3 ± 0.2 and 7.9 ± 0.2. The pH dependence of transitions might be explained in the same manner as AHA 1a conversion (equations 1,2). The calculated values of formal potentials of each reaction were 595 ± 10, 430 ± 12 and 289 ± 8 mV.
The steady state kinetics of the rCiP-catalyzed oxidation of substrates revealed that bimolecular constants of AHAs and AHUs oxidation, with exception of AHA 1e covered range from 8.3 · 103 to 8.7 · 105 M-1s-1 (Table 1). This rate is associated with compound II reduction [7, 8]. Oxidation rate of AHA1e was 2100 times faster in comparison with AHA 1a (Table 1).
A typical approach to evaluate structure activity relationship (SAR) in enzymatic redox reactions is the comparison of reaction rate with redox potential of substrates [2–6]. If electron transfer is predominantly chemically controlled, the reactivity of substrates should be related to the free energy of reaction (to the thermodynamic driving force) as it is the case of other electron transfer reactions . However, only a tendency of the correlation of the log of bimolecular constants with the respective redox potentials could be observed for AHAs benzamides as predicted by electron transfer theory . For acetamides an opposite tendency was found: the constant increased with increasing substrate potential. The reactivity of AHUs didn't correlate with redox potential too. In the attempt to find explanation of unusually low reactivity and independent reactivity to redox potential the docking of AHAs and AHUs in the active center of rCiP was calculated. The docking of the substrates in the active center of rCiP was calculated by taking from the bank the crystal structure of Arthromyces ramosus peroxidase (ARP) since these two enzymes vary only by one terminal amino acid .
In the second cluster, AHA 1a was rotated 180° around N-phenyl axis and a little bit twisted while the water molecule interacted with the C=O group (Figure 4).
The calculations of the docking of AHU 1 revealed that about 95% of the substrate combines in a small pocket at Pro 91 and Ile 153. The mean docking energy of these clusters changes from -21.87 to -26.77 kJ mol-1 or from -21.83 to -25.39 kJ mol-1. The small amount of the substrate was also combined at the residues of Asp 98, Thr 99, and Lys 49. In this complex the N-aryl fragment was located at a distance of 14 Å from the heme and water molecule. The docking of AHU4 produced 3 clusters near Pro 91 and Ile153 at the distal side of heme. A mean docking energy of these clusters was -29.54, -28.9 kJ mol-1 and -27.70 kJ mol-1. The rest cluster was located in the heme pocket. The docking energy of this cluster was -26.90 kJ mol-1. The N-aryl fragment was located at the distance of 4.72 Å from water molecule and O-aryl fragment was 3.73 Å from heme.
The calculations showed that other AHAs, with exception of AHA 1e, and AHUs were docked in the same sites as AHA 1a, AHA 2a, AHU 1 and AHU 4. These sites were formed by hydrophobic amino acid residues. Dockings showed that hydrogen from N-OH group of AHA 1a and AHU 1 does not form strong hydrogen bonds with polypeptide chain. In the case of AHA 2a and AHU 4 N-OH group may form a weak hydrogen bong with Arg 52.
AHA 1e substrate was combined in such manner that hydroxyl group at para position of phenyl ring formed hydrogen bond with proximal His. Different AHA 1e docking in the active center and different reactivity of phenolic and N-OH substrates permits to hypothesize that peroxidase-catalyzed reactivity depends on proton transfer rate in the active center. The electrochemical investigations demonstrated that proton transfer rate determined formal redox potential of the substrates in water solution. In acidic solutions or in aprotic non-polar solvent like methylene chloride the proton transfer was retarded and oxidation of substrate proceeded at high potential. The complexation of substrates in the hydrophobic site of the enzyme also retards proton transfer during oxidation. If proton transfer limits the reaction the local surrounding of N-OH group can be critical to the reaction rate. The combining of AHAs and AHUs in hydrophobic surrounding does not permit fast proton transfer. The larger reactivity of AHA 2a or benzamides may be explained by hydrogen bond formation with Arg 52 and facilitated proton transfer. The combining of AHA 1e in the site where proton transfer is fast due to hydrogen bonds relay to distal His and proximal water makes fast substrate oxidation.
The role of proton transfer in enzymatic oxidoreductase catalysis was described in many publications . During C-H bond oxidation of specific substrates the proton is tunneling along the reaction coordinate and the rate of overall process depends on this process. Recent investigations of peroxidases catalysis permitted to make conclusions that proton transfer with participation of distal water molecule also plays an important role in peroxidases-catalyzed oxidation of physiological substrates. Recently Henriksen and co-workers proved that the oxidation of ferulic acid by plant peroxidase is accompanied by proton transfer with distal water molecule being involved into the process . The conclusion was made that in the case of compound I reduction the final destination of proton was the distal histidine whereas ferryl oxygen was the proton acceptor in the case of compound II reduction. The detail analysis of roles of water made by Jones  permitted to hypothesize that presence or absence of "resident" water molecule at distal site of heme may switch between peroxidase and catalase activities.
The hypothesis formulated in this work states that docking of non-specific substrates may cause troubles of proton transfer relay formation in the active centre of peroxidases and may control the rate of the substrates oxidation if proton transfer is associated with the oxidation reaction. This hypothesis permits to explain unusually low activity of some type of substrates and opens possibilities for new substrates design. Although the experimental results are consistent with the docking calculations more directly target studies such as determination of ternary structure of peroxidase and substrate complex with cyanide are highly desirable.
Materials and Methods
Ab initio calculations of electronic structure and the energies of substrates were performed using a Gaussian 94W package . Geometry optimization of the substrates was accomplished using HF (Hartree-Fock) theory and a 3–21G basis set. The optimized geometry of molecules was used for energies and charges calculations with a 6–31G basis set using RHF and B3PW91 (Density Functional Theory). The atom charges of molecules were calculated with Mulliken algorithm.
The simulations of substrates docking in the active center of Arthromyces ramosus peroxidase (ARP) were performed with AutoDock 3.0 . The crystal data of the complex of ARP with benzhydroxamic acid (HSR)  was downloaded from the Protein Data Bank. All water molecules in the active center of ARP were removed with the exception of a water molecule in the distal side of the heme. Atomic interaction energy grid maps were calculated with 0.25 Å grid spacing and 100 grid points forming 25 Å cubic box centered on the active site of peroxidases on heme side exposed to water. The space of the cubic box covered the active site of peroxidases and the space beyond. The electrostatic interaction energy grid used a distance-dependent dielectric function of Mehler-Solmajer. The docking was accomplished using Lamarckian genetic algorithm. The number of individuals in populations was set to 50. The maximum number of energy evaluations that the genetic algorithm should make was 250 000. Maximum number of generation was 27000. The number of the top individuals that are guaranteed to survive into the next generation was 1. Crossover rate and mutation rate were 0.02 and 0.80.
The authors express sincere thanks to H.-J. Deussen, K. Krikstopaitis and R. Lolck who have made synthesis and have determined the parameters of the compounds. We also thank P. Schneider for the help and support throughout the project and A.J. Olson for the supply of AutoDock 3.0 program.
- Dunford HB: Horseradish peroxidase: structure and kinetic properties. In: Peroxidases in Chemistry and Biology Volume 2 (Edited by Everse J, Everse K), CRC Press, Boca Raton, Ann Arbor, Boston, 1991, 1–25.Google Scholar
- Candeias LP, Folkes LK, Wardman P: Factors controlling the substrate specificity of peroxidases: kinetics and thermodynamics of the reaction of horseradish peroxidase compound I with phenols and indole-3-acetic acids. Biochemistry 1997, 36: 7081–7085. 10.1021/bi970384eView ArticlePubMedGoogle Scholar
- Krikstopaitis K, Kulys J, Pedersen AH, Schneider P: N-substituted p-phenylenediamines as peroxidase and laccase substrates. Acta Chem.Scand 1998, 52: 469–474.View ArticleGoogle Scholar
- Ryabov AD, Goral VN, Ivanova EV, Reshetova MD, Hradsky A, Bildstein B: Linear free-energy relationships and inverted Marcus region in the horseradish peroxidase-catalyzed oxidation of ferrocenes by hydrogen peroxide. J. Organ.Chem 1999, 589: 85–91. 10.1016/S0022-328X(99)00321-6View ArticleGoogle Scholar
- Kulys J, Krikstopaitis K, Ziemys A: Kinetics and thermodynamics of peroxidase-and laccase-catalyzed oxidation of N-substituted phenothiazines and phenoxazines. J. Biol.Inorg.Chem. 2000, 5: 333–340.View ArticlePubMedGoogle Scholar
- Sakurada J, Sekiguchi R, Sato K, Hosoya T: Kinetic and molecular orbital studies on the rate of oxidation of monosubstituted phenols and anilines by horseradish peroxidase compound II. Biochemistry 1990, 29: 4093–4098.View ArticlePubMedGoogle Scholar
- Kulys J, Deussen H-J, Krikstopaitis K, Schneider P, Ziemys A: N-aryl-N-hydroxy urethanes as peroxidase substrates. Mon.Chem 2001, 132: 295–304. 10.1007/s007060170138View ArticleGoogle Scholar
- Kulys J, Deussen H-J, Krikstopaitis K, Lolck R, Schneider P, Ziemys A: N-aryl hydroxamic acids as novel oxidoreductases substrates. Eur. J. Org. Chem. 2001, 3475–3484. Publisher Full Text 10.1002/1099-0690(200109)2001:18<3475::AID-EJOC3475>3.0.CO;2-KGoogle Scholar
- Bixon M, Jortner J: Electron transfer-from isolated molecules to biomolecules. In:Advances in Chemical Physics Volume 106, Part one (Edited by Jortner J, Bixon M) 1999, 35–202.Google Scholar
- Kunishima N, Fukuyama K, Matsubara H, Hatanaka H, Shibano Y, Amachi T: Crystal structure of the fungal peroxidase from Arthromyces ramosus at 1.9 angstroms resolution: structural comparison with the lignin and cytchrome C peroxidases. J.Mol.Biol. 1994, 235: 331–334.View ArticlePubMedGoogle Scholar
- Itakura H, Oda Y, Fukuyama K: Binding made of benzohydroxamic acid to Arthromyces ramosus peroxidase shown by X-ray crystallographic analysis of the complex at 1.6 Å resolution. FEBS Lett 1997, 412: 107–110. 10.1016/S0014-5793(97)00751-5View ArticlePubMedGoogle Scholar
- Kohen A, Klinman JP: Enzyme catalysis: Beyond classical paradigms. Acc. Chem. Res. 1998, 31: 397–404. 10.1021/ar9701225View ArticleGoogle Scholar
- Henriksen A, Smith AT, Gajhede M: Structures of the horseradish peroxidase c-ferulic acid complex and the ternary complex with cyanide suggest how peroxidases oxidize small phenolic substrates. J.Biol. Chem. 1999, 274: 35005–35011. 10.1074/jbc.274.49.35005View ArticlePubMedGoogle Scholar
- Jones P: Roles of water in heme peroxidase and catalase mechanisms. J.Biol.Chem. 2001, 276: 13791–13796.PubMedGoogle Scholar
- Frisch MJ, Trucks GW, Schlegel HB, Gill PMW, Johnson BG, Robb MA, Cheeseman JR, Keith T, Petersson GA, Montgomery JA, et al.: Gaussian 94, Revision E.1,. Gaussian, Inc., Pittsburgh PA 1995.Google Scholar
- Morris GM, Goodsell DS, Huey R, Olson AJ: Distributed automated docking of flexible ligands to proteins: Parallel applications of AudoDock 2.4. J. Comp.-Aided. Mol.Des. 1996, 10: 293–304.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article: verbatim copying and redistribution of this article are permitted in all media for any purpose, provided this notice is preserved along with the article's original URL.