Combined quantum mechanical and molecular mechanical reaction pathway calculation for aromatic hydroxylation by p-hydroxybenzoate-3-hydroxylase

Phenolic hydroxylases

Many flavin-dependent phenolic hydroxylases (monooxygenases) have been extensively investigated. Their crystal structures and reaction mechanisms are well understood. These enzymes belong to groups A and D of the flavin-dependent monooxygenases and can be classified as single-component and two-component flavin-dependent monooxygenases. The insertion of molecular oxygen into the substrates catalyzed by these enzymes is beneficial for modifying the biological properties of phenolic compounds and their derivatives. This chapter provides an in-depth discussion of the structural features of single-component and two-component flavin-dependent phenolic hydroxylases. The reaction mechanisms of selected enzymes, including 3-hydroxy-benzoate 4-hydroxylase (PHBH) and 3-hydroxy-benzoate 6-hydroxylase as representatives of single-component enzymes and 3-hydroxyphenylacetate 4-hydroxylase (HPAH) as a representative of two-component enzymes, are discussed in detail. This chapter comprises the following four main parts: general reaction, structures, reaction mechanisms, and enzyme engineering for biocatalytic applications. Enzymes belonging to the same group catalyze similar reactions but have different unique structural features to control their reactivity to substrates and the formation and stabilization of C4a-hydroperoxyflavin. Protein engineering has been employed to improve the ability to use these enzymes to synthesize valuable compounds. A thorough understanding of the structural and mechanistic features controlling enzyme reactivity is useful for enzyme redesign and enzyme engineering for future biocatalytic applications.

Monooxygenation of aromatic compounds by flavin-dependent monooxygenases

Many flavoenzymes catalyze hydroxylation of aromatic compounds especially phenolic compounds have been isolated and characterized. These enzymes can be classified as either single-component or two-component flavin-dependent hydroxylases (monooxygenases). The hydroxylation reactions catalyzed by the enzymes in this group are useful for modifying the biological properties of phenolic compounds. This review aims to provide an in-depth discussion of the current mechanistic understanding of representative flavin-dependent monooxygenases including 3-hydroxy-benzoate 4-hydroxylase (PHBH, a single-component hydroxylase), 3-hydroxyphenylacetate 4-hydroxylase (HPAH, a two-component hydroxylase), and other monooxygenases which catalyze reactions in addition to hydroxylation, including 2-methyl-3-hydroxypyridine-5-carboxylate oxygenase (MHPCO, a single-component enzyme that catalyzes aromatic-ring cleavage), and HadA monooxygenase (a two-component enzyme that catalyzes additional group elimination reaction). These enzymes have different unique structural features which dictate their reactivity toward various substrates and influence their ability to stabilize flavin intermediates such as C4a-hydroperoxyflavin. Understanding the key catalytic residues and the active site environments important for governing enzyme reactivity will undoubtedly facilitate future work in enzyme engineering or enzyme redesign for the development of biocatalytic methods for the synthesis of valuable compounds.

Hydroxyl Radical-Coupled Electron-Transfer Mechanism of Flavin-Dependent Hydroxylases

Class A flavin-dependent hydroxylases (FdHs) catalyze the hydroxylation of organic compounds in a site- and stereoselective manner. In stark contrast, conventional synthetic routes require environmentally hazardous reagents and give modest yields. Thus, understanding the detailed mechanism of this class of enzymes is essential to their rational manipulation for applications in green chemistry and pharmaceutical production. Both electrophilic substitution and radical intermediate mechanisms have been proposed as interpretations of FdH hydroxylation rates and optical spectra. While radical mechanistic steps are often difficult to examine directly, modern quantum chemistry calculations combined with statistical mechanical approaches can yield detailed mechanistic models providing insights that can be used to differentiate reaction pathways. In the current work, we report quantum mechanical/molecular mechanical (QM/MM) calculations on the fungal TropB enzyme that shows an alternative reaction pathway in which hydroxylation through a hydroxyl radical-coupled electron-transfer mechanism is significantly favored over electrophilic substitution. Furthermore, QM/MM calculations on several modified flavins provide a more consistent interpretation of the experimental trends in the reaction rates seen experimentally for a related enzyme, para-hydroxybenzoate hydroxylase. These calculations should guide future enzyme and substrate design strategies and broaden the scope of biological spin chemistry.

The reaction mechanism of Zika virus NS2B/NS3 serine protease inhibition by dipeptidyl aldehyde: a QM/MM study

Zika virus (ZIKV) infection has become a global public health problem, associated with microcephaly in newborns and Guillain-Barré syndrome in adults. Currently, there are no commercially available anti-ZIKV drugs. The viral protease NS2B/NS3, which is involved in viral replication and maturation, is a potential drug target. Peptidomimetic aldehyde inhibitors bind covalently to the catalytic S135 of the NS3 protease. Here, we apply hybrid quantum mechanics/molecular mechanics (QM/MM) free-energy simulations at the PDDG-PM3/ff14SB level to investigate the inhibition mechanism of the ZIKV protease by a dipeptidyl aldehyde inhibitor (acyl-KR-aldehyde). The results show that proton transfer from the catalytic S135 to H51 occurs in concert with nucleophilic addition on the aldehyde warhead by S135. The anionic covalent complex between the dipeptidyl aldehyde and the ZIKV protease is analogous to the tetrahedral intermediate for substrate hydrolysis. Spontaneous protonation by H51 forms the hemiacetal. In addition, we use correlated ab initio QM/MM potential energy path calculations at levels up to LCCSD(T)/(aug)-cc-pVTZ to obtain accurate potential energy profiles of the reaction, which also support a concerted mechanism. These results provide detailed insight into the mechanism of ZIKV protease inhibition by a peptidyl aldehyde inhibitor, which will guide in the design of inhibitors.

Mechanism of Kynurenine 3-Monooxygenase-Catalyzed Hydroxylation Reaction: A Quantum Cluster Approach

The mechanism of the hydroxylation reaction between l-Kyn and model flavin adenine dinucleotide (FAD)-hydroperoxide was investigated via density functional theory (DFT) calculations in the absence and in the presence of the kynurenine 3-monooxygenase (KMO) enzyme by considering possible pathways that can lead to the product 3-hydroxykynurenine (3-HK). Crystal structure (pdb code: 5NAK )-based calculations involved a quantum cluster model in which the active site of the enzyme with the substrate l-Kyn was represented with 348 atoms. According to the deduced mechanism, KMO-catalyzed hydroxylation reaction takes place with four transformations. In the initial transition state, FAD delivers its peroxy hydroxyl to the l-Kyn ring, creating an sp³-hybridized carbon center. Then, the hydrogen on the hydroxyl moiety is immediately transferred back to the proximal oxygen that remained on FAD. These consequent transformations are in line with the somersault rearrangement previously described for similar enzymatic systems. The second step corresponds to a hydride shift from the sp³-hybridized carbon of the substrate ring to its adjacent carbon, producing the keto form of 3-HK. Then, keto-3-HK is transformed into its enol form (3-HK) with a water-assisted tautomerization. Lastly, FAD is oxidized with a water-assisted dehydration, which also involves 3-HK as a catalyst. In the proposed pathway, Asn54, Pro318, and a crystal water molecule were seen to play significant roles in the proton relays. The energies obtained via the cluster approach were calculated at the B3LYP/6-311+G(2d,2p)//B3LYP/6-31G(d,p) level with solvation (polarizable continuum model) and dispersion (DFT-D3(BJ)) corrections.

Oxidative dehalogenation and denitration by a flavin-dependent monooxygenase is controlled by substrate deprotonation

Enzymes that are capable of detoxifying halogenated phenols (HPs) and nitrophenols (NPs) are valuable for bioremediation and waste biorefining. HadA monooxygenase was found to perform dual functions of oxidative dehalogenation (hydroxylation plus halide elimination) and denitration (hydroxylation plus nitro elimination). Rate constants associated with individual steps of HadA reactions with phenol, halogenated phenols and nitrophenols were measured using combined transient kinetic approaches of stopped-flow absorbance/fluorescence and rapid-quench flow techniques. Density functional theory was used to calculate the thermodynamic and electronic parameters associated with hydroxylation and group elimination steps. These parameters were correlated with the rate constants of hydroxylation, group elimination, and overall product formation to identify factors controlling individual steps. The results indicated that the hydroxylation rate constant is higher when the pK a of the phenolic group is lower, i.e. it is more easily deprotonated, but not higher when the energy gap between the E LUMO of the C4a-hydroperoxy-FAD intermediate and the E HOMO of the phenolate substrate is lower. These data suggest that the substrate deprotonation has a higher energy barrier than the -OH transfer, and thus controls the hydroxylation step. For the group elimination, the process is controlled by the ability of the C-X bond to break. For the overall product formation (hydroxylation and group elimination combined), this analysis showed that the rate constant of product formation is dependent on the pK a value of the substrate, indicating that the overall reaction is controlled by substrate deprotonation. This step also likely has the highest energy barrier and thus controls the overall process of oxidative dehalogenation and denitration by HadA. This report is the first to identify a key mechanistic factor controlling the enzymatic processes of oxidative dehalogenation and denitration.

Pathways for Arene Oxidation in Non-Heme Diiron Enzymes: Lessons from Computational Studies on Benzoyl Coenzyme A Epoxidase

Oxygenation of aromatic rings using O₂ is catalyzed by several non-heme carboxylate-bridged diiron enzymes. In order to provide a general mechanistic description for these reactions, computational studies were carried out at the ONIOM(B3LYP/BP86/Amber) level on the non-heme diiron enzyme benzoyl coenzyme A epoxidase, BoxB. The calculations revealed four possible pathways for attacking the aromatic ring: (a) electrophilic (2e^-) attack by a bis(μ-oxo)-diiron(IV) species (Q pathway); (b) electrophilic (2e^-) attack via the σ* orbital of a μ-η²:η²-peroxo-diiron(III) intermediate (Pσ* pathway); (c) radical (1e^-) attack via the π*-orbital of a superoxo-diiron(II,III) species (Pπ* pathway); (d) radical (1e^-) attack of a partially quenched bis(μ-oxo)-diiron(IV) intermediate (Q' pathway). The results allowed earlier work of de Visser on olefin epoxidation by diiron complexes and QM-cluster studies of Liao and Siegbahn on BoxB to be put into a broader perspective. Parallels with epoxidation using organic peracids were also examined. Specifically for the BoxB enzyme, the Q pathway was found to be the most preferred, but the corresponding bis(μ-oxo)-diiron(IV) species is significantly destabilized and not expected to be directly observable. Epoxidation via the Pσ* pathway represents an energetically somewhat higher lying alternative; possible strategies for experimental discrimination are discussed. The selectivity toward epoxidation is shown to stem from a combination of inherent electronic properties of the thioacyl substituent and enzymatic constraints. Possible implications of the results for toluene monooxygenases are considered as well.

High-level QM/MM calculations support the concerted mechanism for Michael addition and covalent complex formation in thymidylate synthase

Thymidylate synthase (TS) is a promising cancer target, due to its crucial function in thymine synthesis. It performs the reductive methylation of 2'-deoxyuridine-5'-phosphate (dUMP) to thymidine-5'-phosphate (dTMP), using N-5,10-methylene-5,6,7,8-tetrahydrofolate (mTHF) as a cofactor. After the formation of the dUMP/mTHF/TS noncovalent complex, and subsequent conformational activation, this complex has been proposed to react via nucleophilic attack (Michael addition) by Cys146, followed by methylene-bridge formation to generate the ternary covalent intermediate. Herein, QM/MM (B3LYP-D/6-31+G(d)-CHARMM27) methods are used to model the formation of the ternary covalent intermediate. A two-dimensional potential energy surface reveals that the methylene-bridged intermediate is formed via a concerted mechanism, as indicated by a single transition state on the minimum energy pathway and the absence of a stable enolate intermediate. A range of different QM methods (B3LYP, MP2 and SCS-MP2, and different basis sets) are tested for the calculation of the activation energy barrier for the formation of the methylene-bridged intermediate. We test convergence of the QM/MM results with respect to size of the QM region. Inclusion of Arg166, which interacts with the nucleophilic thiolate, in the QM region is important for reliable results; the MM model apparently does not reproduce energies for distortion of the guanidinium side chain correctly. The spin component scaled-Møller-Plessett perturbation theory (SCS-MP2) approach was shown to be in best agreement (within 1.1 kcal/mol) while the results obtained with MP2 and B3LYP also yielded acceptable values (deviating by less than 3 kcal/mol) compared with the barrier derived from experiment. Our results indicate that using a dispersion-corrected DFT method, or a QM method with an accurate treatment of electron correlation, increases the agreement between the calculated and experimental activation energy barriers, compared with the semiempirical AM1 method. These calculations provide important insight into the reaction mechanism of TS and may be useful in the design of new TS inhibitors.

Form follows function: structural and catalytic variation in the class a flavoprotein monooxygenases

Flavoprotein monooxygenases (FPMOs) exhibit an array of mechanistic solutions to a common chemical objective; the monooxygenation of a target substrate. Each FPMO efficiently couples reduction of a flavin cofactor by NAD(P)H to oxygenation of the target substrate via a (hydro)peroxyflavin intermediate. This purpose of this review is to describe in detail the Class A flavoprotein hydroxylases (FPMO) in the context of the other FPMO classes (B–F). Both one and two component FPMOs are found in nature. Two-component enzymes require, in addition to the monooxygenase, the involvement of a reductase that first catalyzes the reduction of the flavin by NAD(P)H. The Class A and B FPMOs are single-component and manage to orchestrate the same net reaction within a single peptide. The Class A enzymes have, by some considerable margin, the most complete research record. These enzymes use choreographed movements of the flavin ring that facilitate access of the organic substrates to the active site, provide a means for interaction of NADPH with the flavin, offer a mechanism to sequester the dioxygen reduction chemistry from solvent and a means to release the product. The majority of the discrete catalytic events of the catalytic cycle can be observed directly in exquisite detail using spectrophotometric kinetic methods and many of the key mechanistic conclusions are further supported by structural data. This review attempts to compile each of the key observations made for both paradigm and newly discovered examples of Class A FPMOs into a complete catalytic description of one enzymatic turnover.

QM/MM methods for biomolecular systems

Combined quantum-mechanics/molecular-mechanics (QM/MM) approaches have become the method of choice for modeling reactions in biomolecular systems. Quantum-mechanical (QM) methods are required for describing chemical reactions and other electronic processes, such as charge transfer or electronic excitation. However, QM methods are restricted to systems of up to a few hundred atoms. However, the size and conformational complexity of biopolymers calls for methods capable of treating up to several 100,000 atoms and allowing for simulations over time scales of tens of nanoseconds. This is achieved by highly efficient, force-field-based molecular mechanics (MM) methods. Thus to model large biomolecules the logical approach is to combine the two techniques and to use a QM method for the chemically active region (e.g., substrates and co-factors in an enzymatic reaction) and an MM treatment for the surroundings (e.g., protein and solvent). The resulting schemes are commonly referred to as combined or hybrid QM/MM methods. They enable the modeling of reactive biomolecular systems at a reasonable computational effort while providing the necessary accuracy.

Analysis of Classical and Quantum Paths for Deprotonation of Methylamine by Methylamine Dehydrogenase

The hydrogen-transfer reaction catalysed by methylamine dehydrogenase (MADH) with methylamine (MA) as substrate is a good model system for studies of proton tunnelling in enzyme reactions--an area of great current interest--for which atomistic simulations will be vital. Here, we present a detailed analysis of the key deprotonation step of the MADH/MA reaction and compare the results with experimental observations. Moreover, we compare this reaction with the related aromatic amine dehydrogenase (AADH) reaction with tryptamine, recently studied by us, and identify possible causes for the differences observed in the measured kinetic isotope effects (KIEs) of the two systems. We have used combined quantum mechanics/molecular mechanics (QM/MM) techniques in molecular dynamics simulations and variational transition state theory with multidimensional tunnelling calculations averaged over an ensemble of paths. The results reveal important mechanistic complexity. We calculate activation barriers and KIEs for the two possible proton transfers identified-to either of the carboxylate oxygen atoms of the catalytic base (Asp428beta)-and analyse the contributions of quantum effects. The activation barriers and tunnelling contributions for the two possible proton transfers are similar and lead to a phenomenological activation free energy of 16.5+/-0.9 kcal mol(-1) for transfer to either oxygen (PM3-CHARMM calculations applying PM3-SRP specific reaction parameters), in good agreement with the experimental value of 14.4 kcal mol(-1). In contrast, for the AADH system, transfer to the equivalent OD1 was found to be preferred. The structures of the enzyme complexes during reaction are analysed in detail. The hydrogen bond of Thr474beta(MADH)/Thr172beta(AADH) to the catalytic carboxylate group and the nonconserved active site residue Tyr471beta(MADH)/Phe169beta(AADH) are identified as important factors in determining the preferred oxygen acceptor. The protein environment has a significant effect on the reaction energetics and hence on tunnelling contributions and KIEs. These environmental effects, and the related clearly different preferences for the two carboxylate oxygen atoms (with different KIEs) in MADH/MA and AADH/tryptamine, are possible causes of the differences observed in the KIEs between these two important enzyme reactions.

Tunneling and Classical Paths for Proton Transfer in an Enzyme Reaction Dominated by Tunneling: Oxidation of Tryptamine by Aromatic Amine Dehydrogenase

Proton tunneling dominates the oxidative deamination of tryptamine catalyzed by the enzyme aromatic amine dehydrogenase. For reaction with the fast substrate tryptamine, a H/D kinetic isotope effect (KIE) of 55 +/- 6 has been reported-one of the largest observed in an enzyme reaction. We present here a computational analysis of this proton-transfer reaction, applying combined quantum mechanics/molecular mechanics (QM/MM) methods (PM3-SRP//PM3/CHARMM22). In particular, we extend our previous computational study (Masgrau et al. Science 2006, 312, 237) by using improved energy corrections, high-level QM/MM methods, and an ensemble of paths to estimate the tunneling contributions. We have carried out QM/MM molecular dynamics simulations and variational transition state theory calculations with small-curvature tunneling corrections. The results provide detailed insight into the processes involved in the reaction. Transfer to the O2 oxygen of the catalytic base, Asp128beta, is found to be the favored reaction both thermodynamically and kinetically, even though O1 is closer in the reactant complex. Comparison of quantum and classical models of proton transfer allows estimation of the contribution of hydrogen tunneling in lowering the barrier to reaction in the enzyme. A reduction of the activation free energy due to tunneling of 3.1 kcal mol-1 is found, which represents a rate enhancement due to tunneling by 2 orders of magnitude. The calculated KIE of 30 is significantly elevated over the semiclassical limit, in agreement with the experimental observations; a semiclassical value of 6 is obtained when tunneling is omitted. A polarization of the C-H bond to be broken is observed due to the close proximity of the catalytic aspartate and the (formally) positively charged imine nitrogen. A comparison is also made with the related quinoprotein methylamine dehydrogenase (MADH)-the much lower KIE of 11 that we obtain for the MADH/methylamine system is found to arise from a more endothermic potential energy surface for the MADH reaction.

Crystal structure of 3-hydroxybenzoate hydroxylase from Comamonas testosteroni has a large tunnel for substrate and oxygen access to the active site

The 3-hydroxybenzoate hydroxylase (MHBH) from Comamonas testosteroni KH122-3s is a single-component flavoprotein monooxygenase, a member of the glutathione reductase (GR) family. It catalyzes the conversion of 3-hydroxybenzoate to 3,4-dihydroxybenzoate with concomitant requirements for equimolar amounts of NADPH and molecular oxygen. The production of dihydroxy-benzenoid derivative by hydroxylation is the first step in the aerobic degradation of various phenolic compounds in soil microorganisms. To establish the structural basis for substrate recognition, the crystal structure of MHBH in complex with its substrate was determined at 1.8 A resolution. The enzyme is shown to form a physiologically active homodimer with crystallographic 2-fold symmetry, in which each subunit consists of the first two domains comprising an active site and the C-terminal domain involved in oligomerization. The protein fold of the catalytic domains and the active-site architecture, including the FAD and substrate-binding sites, are similar to those of 4-hydroxybenzoate hydroxylase (PHBH) and phenol hydroxylase (PHHY), which are members of the GR family, providing evidence that the flavoprotein aromatic hydroxylases share similar catalytic actions for hydroxylation of the respective substrates. Structural comparison of MHBH with the homologous enzymes suggested that a large tunnel connecting the substrate-binding pocket to the protein surface serves for substrate transport in this enzyme. The internal space of the large tunnel is distinctly divided into hydrophilic and hydrophobic regions. The characteristically stratified environment in the tunnel interior and the size of the entrance would allow the enzyme to select its substrate by amphiphilic nature and molecular size. In addition, the structure of the Xe-derivative at 2.5 A resolution led to the identification of a putative oxygen-binding site adjacent to the substrate-binding pocket. The hydrophobic nature of the xenon-binding site extends to the solvent through the tunnel, suggesting that the tunnel could be involved in oxygen transport.

Computer-aided modeling of pentachlorophenol 4-monooxygenase and site-directed mutagenesis of its active site

Homology modeling was used to construct a model of the three-dimensional structure of pentachlorophenol 4-monooxygenase (PcpB). A PSI-BLAST homology search was initially performed to identify the 3D structure of proteins homologous with PcpB. The feasibility of modeled structures of PcpB was evaluated by Verify3D, which calculated structural compatibility scores based on 3D-1D profiles. The predicted structure of PcpB had an acceptable 3D-1D self-compatibility score, beyond the incorrect fold score threshold. A PcpB-pentachlorophenol (PCP) complex was then constructed utilizing the modeled PcpB structure. After energy minimization of the complex, and successive minimizations of the system that consisted of the complex and the water layer surrounding the complex, the molecular dynamics of the system were simulated. The active-site residues of PcpB were identified on the basis of the modeled structure, and PcpB mutants were then designed to change the active site residues, expressed, and purified by affinity chromatography. The mutant activity was compared with that of the wild-type to investigate the validity of the modeled structure. The experimental results suggested that Phe85, Tyr216, and Arg235 were relevant to enzyme activity, and that Tyr397 and Phe87 were important for stabilization of the structure of PcpB.

Studies of the mechanism of phenol hydroxylase: effect of mutation of proline 364 to serine

An active site residue in phenol hydroxylase (PHHY), Pro364, was mutated to serine to investigate its role in enzymatic catalysis. In the presence of phenol, the reaction between the reduced flavin of P364S and oxygen is very fast, but only 13% of the flavin is utilized to hydroxylate the substrate, compared to nearly 100% for the wild-type enzyme. The oxidative half-reaction of PHHY using m-cresol as a substrate is similarly affected by the mutation. Pro364 was suggested to be important in stabilizing the transition state of the oxygen transfer step by forming a hydrogen bond between its carbonyl oxygen and the C4a-hydroperoxyflavin [Ridder, L., Mullholland, A. J., Rietjens, I. M. C. M., and Vervoort, J. (2000) J. Am. Chem. Soc. 122, 8728-8738]. The P364S mutation may weaken this interaction by increasing the flexibility of the peptide chain; hence, the transition state would be destabilized to result in a decreased level of hydroxylation of phenol. However, when the oxidative half-reaction was studied using resorcinol as a substrate, the P364S mutant form was not significantly different from the wild-type enzyme. The rate constants for all the reaction steps as well as the hydroxylation efficiency (coupling between NADPH oxidation and resorcinol consumption) are comparable to those of the wild-type enzyme. It is suggested that the function of Pro364 in catalysis, stabilization of the transition state, is not as important in the reaction with resorcinol, possibly because the position of hydroxylation is different with resorcinol than with phenol and m-cresol.

Quantum mechanical methods for enzyme kinetics

This review discusses methods for the incorporation of quantum mechanical effects into enzyme kinetics simulations in which the enzyme is an explicit part of the model. We emphasize three aspects: (a) use of quantum mechanical electronic structure methods such as molecular orbital theory and density functional theory, usually in conjunction with molecular mechanics; (b) treating vibrational motions quantum mechanically, either in an instantaneous harmonic approximation, or by path integrals, or by a three-dimensional wave function coupled to classical nuclear motion; (c) incorporation of multidimensional tunneling approximations into reaction rate calculations.

Computational enzymology

Recent advances in computational methods and the availability of fast, affordable computers have made the modeling of enzymatic reactions practical. The remaining challenges include achieving the accuracy level at which thermodynamic parameters and kinetic constants for different substrates, mutant enzymes, or in the presence of allosteric effectors can be predicted quantitatively.

Modelling flavin and substrate substituent effects on the activation barrier and rate of oxygen transfer by p-hydroxybenzoate hydroxylase

The simulation of enzymatic reactions, using computer models, is becoming a powerful tool in the most fundamental challenge in biochemistry: to relate the catalytic activity of enzymes to their structure. In the present study, various computed parameters were correlated with the natural logarithm of experimental rate constants for the hydroxylation of various substrate derivatives catalysed by wild-type para-hydroxybenzoate hydroxylase (PHBH) as well as for the hydroxylation of the native substrate (p-hydroxybenzoate) by PHBH reconstituted with a series of 8-substituted flavins. The following relative parameters have been calculated and tested: (a) energy barriers from combined quantum mechanical/molecular mechanical (QM/MM) (AM1/CHARMM) reaction pathway calculations, (b) gas-phase reaction enthalpies (AM1) and (c) differences between the HOMO and LUMO energies of the isolated substrate and cofactor molecules (AM1 and B3LYP/6-31+G(d)). The gas-phase approaches yielded good correlations, as long as similarly charged species are involved. The QM/MM approach resulted in a good correlation, even including differently charged species. This indicates that the QM/MM model accounts quite well for the solvation effects of the active site surroundings, which vary for differently charged species. The correlations obtained demonstrate quantitative structure activity relationships for an enzyme-catalysed reaction including, for the first time, substitutions on both substrate and cofactor.