Substrate Binding Mode Research Articles

X-ray crystallography of mutant GDSL esterase S12A of Photobacterium marinum J15.

GDSL esterase is designated as a member of Family II of lipolytic enzymes known to catalyse the synthesis and hydrolysis of ester bonds. The enzyme possesses a highly conserved motif Ser-Gly-Asn-His in the four conserved blocks I, II, III and V respectively. The enzyme characteristics, such as region-, chemo-, and enantioselectivity, help in resolving the racemic mixture of single-isomer chiral drugs. Recently, crystal structure of GDSL esterase from Photobacterium J15 has been reported (PDB ID: 5XTU) but not in complex with substrate. Therefore, GDSL in complex with substrate could provide insights into the binding mode of substrate towards inactive form of GDSL esterase (S12A) and identify the hot spot residues for the designing of a better binding pocket. Insight into molecular mechanisms is limited due to the lack of crystal structure of GDSL esterase-substrate complex. In this paper, the crystallization of mutant GDSL esterase (S12A) (PDB ID: 8HWO) and its complex with butyric acid (PDB ID: 8HWP) are reported. The optimized structure would be vital in determining hot spot residue for GDSL esterase. This preliminary study provides an understanding of the interactions between enzymes and hydrolysed p-nitro-phenyl butyrate. The information could guide in the rational design of GDSL esterase in overcoming the medical limitations associated with racemic mixture.

Structure and function of the metagenomic plastic-degrading polyester hydrolase PHL7 bound to its product

The recently discovered metagenomic-derived polyester hydrolase PHL7 is able to efficiently degrade amorphous polyethylene terephthalate (PET) in post-consumer plastic waste. We present the cocrystal structure of this hydrolase with its hydrolysis product terephthalic acid and elucidate the influence of 17 single mutations on the PET-hydrolytic activity and thermal stability of PHL7. The substrate-binding mode of terephthalic acid is similar to that of the thermophilic polyester hydrolase LCC and deviates from the mesophilic IsPETase. The subsite I modifications L93F and Q95Y, derived from LCC, increased the thermal stability, while exchange of H185S, derived from IsPETase, reduced the stability of PHL7. The subsite II residue H130 is suggested to represent an adaptation for high thermal stability, whereas L210 emerged as the main contributor to the observed high PET-hydrolytic activity. Variant L210T showed significantly higher activity, achieving a degradation rate of 20 µm h−1 with amorphous PET films.

Structural insights into the substrate binding sites of O-carbamoyltransferase VtdB from Streptomyces sp. NO1W98

The O-carbamoyltransferase VtdB catalyzes the carbamoylation of venturicidin B, which is essential for the biosynthesis of the antibiotic venturicidin A. Here, the crystal structures of VtdB and VtdB in complex with the intermediate carbamoyladenylate (VtdBCAO) were determined at resolutions of 2.99 Å and 2.90 Å, respectively. The structures resemble the conserved YrdC-like and specific Kae1-like domains. A magnesium ion and the intermediate carbamoyladenylate were also observed in the Kae1-like domain of VtdB. The structure of VtdBCAO in complex with the substrate venturicidin B was modeled by a molecular docking method to better understand the substrate binding mode, revealing a novel venturicidin B binding pocket.

Complete bio-degradation of poly(butylene adipate-co-terephthalate) via engineered cutinases

Poly(butylene adipate-co-terephthalate) (PBAT), a polyester made of terephthalic acid (TPA), 1,4-butanediol, and adipic acid, is extensively utilized in plastic production and has accumulated globally as environmental waste. Biodegradation is an attractive strategy to manage PBAT, but an effective PBAT-degrading enzyme is required. Here, we demonstrate that cutinases are highly potent enzymes that can completely decompose PBAT films in 48 h. We further show that the engineered cutinases, by applying a double mutation strategy to render a more flexible substrate-binding pocket exhibit higher decomposition rates. Notably, these variants produce TPA as a major end-product, which is beneficial feature for the future recycling economy. The crystal structures of wild type and double mutation of a cutinase from Thermobifida fusca in complex with a substrate analogue are also solved, elucidating their substrate-binding modes. These structural and biochemical analyses enable us to propose the mechanism of cutinase-mediated PBAT degradation.

Structural and Mechanistic Studies on Substrate and Stereoselectivity of the Indole Monooxygenase VpIndA1: New Avenues for Biocatalytic Epoxidations and Sulfoxidations.

Flavoprotein monooxygenases are a versatile group of enzymes for biocatalytic transformations. Among these, group E monooxygenases (GEMs) catalyze enantioselective epoxidation and sulfoxidation reactions. Here, we describe the crystal structure of an indole monooxygenase from the bacterium Variovorax paradoxus EPS, a GEM designated as VpIndA1. Complex structures with substrates reveal productive binding modes that, in conjunction with force-field calculations and rapid mixing kinetics, reveal the structural basis of substrate and stereoselectivity. Structure-based redesign of the substrate cavity yielded variants with new substrate selectivity (for sulfoxidation of benzyl phenyl sulfide) or with greatly enhanced stereoselectivity (from 35.1 % to 99.8 % ee for production of (1S,2R)-indene oxide). This first determination of the substrate binding mode of GEMs combined with structure-function relationships opens the door for structure-based design of these powerful biocatalysts.

To the Understanding of Catalysis by D-Amino Acid Transaminases: A Case Study of the Enzyme from Aminobacterium colombiense

Pyridoxal-5'-phosphate (PLP)-dependent transaminases are highly efficient biocatalysts for stereoselective amination. D-amino acid transaminases can catalyze stereoselective transamination producing optically pure D-amino acids. The knowledge of substrate binding mode and substrate differentiation mechanism in D-amino acid transaminases comes down to the analysis of the transaminase from Bacillus subtilis. However, at least two groups of D-amino acid transaminases differing in the active site organization are known today. Here, we present a detailed study of D-amino acid transaminase from the gram-negative bacterium Aminobacterium colombiense with a substrate binding mode different from that for the transaminase from B. subtilis. We study the enzyme using kinetic analysis, molecular modeling, and structural analysis of holoenzyme and its complex with D-glutamate. We compare the multipoint binding of D-glutamate with the binding of other substrates, D-aspartate and D-ornithine. QM/MM MD simulation reveals that the substrate can act as a base and its proton can be transferred from the amino group to the α-carboxylate group. This process occurs simultaneously with the nucleophilic attack of the PLP carbon atom by the nitrogen atom of the substrate forming gem-diamine at the transimination step. This explains the absence of the catalytic activity toward (R)-amines that lack an α-carboxylate group. The obtained results clarify another substrate binding mode in D-amino acid transaminases and underpinned the substrate activation mechanism.

Structural insight into the substrate-binding mode and catalytic mechanism for MlrC enzyme of Sphingomonas sp. ACM-3962 in linearized microcystin biodegradation.

Removing microcystins (MCs) safely and effectively has become an urgent global problem because of their extremely hazardous to the environment and public health. Microcystinases derived from indigenous microorganisms have received widespread attention due to their specific MC biodegradation function. However, linearized MCs are also very toxic and need to be removed from the water environment. How MlrC binds to linearized MCs and how it catalyzes the degradation process based on the actual three-dimensional structure have not been determined. In this study, the binding mode of MlrC with linearized MCs was explored using a combination of molecular docking and site-directed mutagenesis methods. A series of key substrate binding residues, including E70, W59, F67, F96, S392 and so on, were identified. Sodium dodecane sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) was used to analyze samples of these variants. The activity of MlrC variants were measured using high performance liquid chromatography (HPLC). We used fluorescence spectroscopy experiments to research the relationship between MlrC enzyme (E), zinc ion (M), and substrate (S). The results showed that MlrC enzyme, zinc ion and substrate formed E-M-S intermediates during the catalytic process. The substrate-binding cavity was made up of N and C-terminal domains and the substrate-binding site mainly included N41, E70, D341, S392, Q468, S485, R492, W59, F67, and F96. The E70 residue involved in both substrate catalysis and substrate binding. In conclusion, a possible catalytic mechanism of the MlrC enzyme was further proposed based on the experimental results and a literature survey. These findings provided new insights into the molecular mechanisms of the MlrC enzyme to degrade linearized MCs, and laid a theoretical foundation for further biodegradation studies of MCs.

Elucidating the Roles of Distinct Chemical Factors in the Hydrolytic Activities of Hetero- and Homonuclear Synthetic Analogues of Binuclear Metalloenzymes

In this study, hydrolytic activities of hetero- and homobinuclear metallovariants of an asymmetric (I) or symmetric (II) ligand with the FeIII-ZnII, ZnII-ZnII, and CuII-CuII cores (i.e., IFZ, IZZ, and ICC or IIFZ, IIZZ, and IICC, respectively) are investigated using DFT calculations through four distinct mechanisms: dissociative (DA), substrate-assisted (SA), water-assisted (WA), and associative (AS). Additionally, the effects of different nucleophiles (μ-OH, terminal-OH, and -O2H3), coordination numbers, para substituents (-CH3, -Cl, and -NO2) of the linker, and an external electric field on the energetics of these reactions are computed. The geometries, spin ground states, and substrate binding modes of the three metal centers for both asymmetric (I) and symmetric (II) ligands differ from each other. There is no Lewis acid activation of the bis(2,4-dinitrophenyl) phosphate (BDNPP) substrate, and hydrolysis is predominantly controlled by the nucleophilicity of the metal-bound hydroxyl ion. The electronic nature of the metal ions determines the activities of their complexes, and the homobinuclear IZZ is found to be the most active complex. However, complexes formed with ligand I are not more active than their ligand-II-containing counterparts for all metal ions. The DFT calculations suggest that the DA pathway is the energetically most feasible among the four pathways. The terminal hydroxyl group is the strongest nucleophile, and the electron-donating -CH3 group is the most suitable para substituent in the linker. Whereas the introduction of an external electric field along the reaction axis lowers the barrier for IFZ, it leaves that unchanged for IZZ and increases that for ICC. These combined results in this study highlight the influence of distinct critical chemical factors such as the electronic nature of the metal ions, the ligand environment, as well as the linker and nucleophile on phosphoester hydrolysis. Insights gained will guide the design of the next generation of versatile metal complexes for a wide range of reactions and applications.

Conformational heterogeneity suggests multiple substrate binding modes in CYP106A2

CYP106A2 (cytochrome P450meg) is a bacterial enzyme originally isolated from B. megaterium, and has been shown to hydroxylate a wide variety of substrates, including steroids. The regio- and stereochemistry of CYP106A2 hydroxylation has been shown to be dependent on a variety of factors, and hydroxylation often occurs at more than one site and/or with lack of stereospecificity for some substrates. Comprehensive backbone 15N, 1H and 13C resonance assignments based on multidimensional nuclear magnetic resonance (NMR) experiments performed with uniform and selective isotopically labeled CYP106A2 samples are reported herein, and broadening and splitting of resonances assigned to regions of the enzyme shown to be affected by substrate binding in other P450 enzymes indicate that substrate binding does not reduce structural heterogeneity as has been observed previously in P450 enzymes CYP101A1 and MycG. Paramagnetic relaxation enhancement (PRE) due to proximity between substrate protons and the heme iron were measured for three different substrates, and the relatively uniform nature of the PREs support the proposal that multiple substrate binding modes are occupied at saturating substrate concentrations.

New insights into the substrate inhibition of human 17β-hydroxysteroid dehydrogenase type 1

Human type 1 17β-hydroxysteroid dehydrogenase (17β-HSD1),a member of the short-chain dehydrogenase/reductase family, catalyzes the last step in the bioactivation of the most potent estrogen estradiol with high specificity and is thus involved in estrogen-dependent diseases. As an oxidoreductase, 17β-HSD1 can utilize both triphosphate and diphosphate cofactors in reaction at the molecular level, but more specific with triphosphate cofactor. The NADPH is much higher than NADP+ in living cells leading to preliminary reduction action. The enzyme also showed substrate-induced inhibition unprecedented in other members of 17β-HSDs. Our previous study elucidated the structural mechanism of substrate inhibition is due to the reversely bound estrone (E1) in the substrate-binding pocket of the enzyme resulting in a dead-end complex. However, the effect of the cofactor preference on the substrate inhibition of the enzyme is not yet clear. In the present study, we solved the ternary crystal structures of 17β-HSD1 in complex with E1 and cofactor analog NAD+ . Combined with molecular dynamics simulation using the enzyme with NADH/NADPH and different oriented E1 (normally oriented, E1N; reversely oriented, E1R), such ternary structure provides a complete picture of enzyme-substrate-cofactor interactions. The results reveal that different cofactors and substrate binding mode affect the allosteric effect between the two subunits of the enzyme. And the results from MD simulations confirmed that His221 plays a key role in the formation of dead-end complex in NADPH complex, and the absence of stable interaction between His221 and E1R in the NADH complex should be the main reason for its lack of substrate inhibition.

Characterisation and substrate binding modes of exopolygalacturonase PGQ1 from Saccharobesus litoralis

Among the enzymes required for the efficient utilisation of pectin is polygalacturonase. Saccharobesus litoralis harbours two polygalacturonases belonging to glycoside hydrolase family 28 (GH28). One of them, PGQ1, cleaved polygalacturonate exolytically at the non-reducing end into monomeric units. It was most active at 60 °C and pH 8, with K m and k cat values of 2.3 mg/ml and 6.4 s−1 respectively. Its homology model of a right-handed parallel β-helix core consisted of Asp297 as the general acid and either Asp276 or Asp298 as the general base. By inferring the substrate binding modes at the −1 and +1 subsites from known crystal structures, a hexagalacturonate could be docked into the highly electropositive binding cleft. Interestingly, while no residues were present in the vicinity to make up the +2 and +4 subsites, Arg361 and Arg430 could readily bind to the carboxyl groups of the galacturonates at the +3 and +5 subsites respectively. Structural comparison suggested that this binding pattern with missing subsites might be unique to closely related exopolygalacturonases. As S. litoralis grew much more slowly on extracellular galacturonate due to the lack of a transporter for the monosaccharide, PGQ1 probably functioned in the periplasm to help degrade oligopectates completely. Communicated by Ramaswamy H. Sarma

Rationally Controlling Selective Steroid Hydroxylation via Scaffold Sampling of a P450 Family

Many steroids are important pharmaceutically active compounds, while cytochrome P450 monooxygenases (CYPs) are attractive enzymes for applications in steroidal drug synthesis. However, the catalytic efficiency of existing P450s is not routinely high enough, as well as the molecular basis for selectivity control is unclear, which severely restrict their real applications. Here, a 16β steroid-hydroxylase CYP109B4 from Bacillus sonorensis is identified with excellent selectivity and activity. The crystallization and structural analysis of CYP109B4 reveal potential three “hotspot” residues (V84, V292, and S387) responsible for selectivity control. Then, guided by the sequence–function relationships revealed from the mutability landscape construction on the three residues, focused rational iterative site-specific mutagenesis (FRISM) and limited iterative saturation mutagenesis were performed, which provide variant B4-M7 (L240V/S387F/V84L/V292S/I291T/M290F/F294I) with completely switched regioselectivity from 16β to 15β. The subsequent computational analysis uncovers insights into the substrate binding modes in CYP109B4 and its variants, which further confirms the critical role of the “hotspot” residues for selectivity control. Finally, the generality of conserved-“hotspots”-mediated selectivity control is demonstrated by performing scaffold sampling between a panel of CYP109B members. Overall, in addition to the present chemical results, our study provides guidance in rationally designing more excellent P450 biocatalysts for potential practical (industrial) applications.

Structural and Functional Characterization of the Ureidoacrylate Amidohydrolase RutB from Escherichia coli

We present a detailed structure-function analysis of the ureidoacrylate amidohydrolase RutB from Eschericha coli, which is an essential enzyme of the Rut pathway for pyrimidine utilization. Crystals of selenomethionine-labeled RutB were produced, which allowed us to determine the first structure of the enzyme at a resolution of 1.9 Å and to identify it as a new member of the isochorismatase-like hydrolase family. RutB was co-crystallized with the substrate analogue ureidopropionate, revealing the mode of substrate binding. Mutation of residues constituting the catalytic triad (D24A, D24N, K133A, C166A, C166S, C166T, C166Y) resulted in complete inactivation of RutB, whereas mutation of other residues close to the active site (Y29F, Y35F, N72A, W74A, W74F, E80A, E80D, S92A, S92T, S92Y, Q105A, Y136A, Y136F) leads to distinct changes of the turnover number (kcat) and/or the Michaelis constant (KM). The results of our structural and mutational studies allowed us to assign specific functions to individual residues and to formulate a plausible reaction mechanism for RutB.

New small-molecule alcohol synthesis by breaking the space limitation of the "aromatic cage" in Pseudomonas sp. AK1 BBOX.

Fe(II)/2OG-dependent oxygenase γ-butyrobetaine hydroxylase (BBOX) stereoselectively hydroxylates inactive C-H bonds and produces L-carnitine. It has potential applications in the biosynthesis of L-carnitine and the synthesis of other small molecule alcohols. In this paper, we systematically explore the substrate range of Pseudomonas sp. AK1 BBOX (psBBOX), with emphasis on the quaternary ammonium portion of γ-butyrobetaine (γ-BB). The space limitation of the "aromatic cage" in psBBOX in the hydroxylation of large quaternary ammonium analogues was studied, and the role of four aromatic amino acid residues in the substrate binding mode was analyzed. Consequently, the F188A mutant was developed with the ability to hydroxylate cyclic quaternary ammonium analogues and generate new alcohol compounds by breaking the limitation of the "aromatic cage".

Structure prediction-based insights into the patatin family of Mycobacterium tuberculosis.

Despite its genome sequencing more than two decades ago, the majority of the genes of Mycobacterium tuberculosis remain functionally uncharacterized. Patatins are one such class of proteins that, despite undergoing an expansion in this pathogenic species compared to their non-pathogenic cousins, remain largely unstudied. Recent advances in protein structure prediction using machine learning tools such as AlphaFold2 have provided high-confidence predicted structures for all M. tuberculosis proteins. Here we present detailed analyses of the patatin family of M. tuberculosis using AlphaFold-predicted structures, providing insights into likely modes of regulation, membrane interaction and substrate binding. Regulatory domains within this family of proteins include cyclic nucleotide binding, lid-like domains and other helical domains. Using structural homologues, we identified the likely membrane localization mechanisms and substrate-binding sites. These analyses reveal diversity in their regulatory capacity, mechanisms of membrane binding and likely length of fatty acid substrates. Together, this analysis suggests unique roles for the eight predicted patatins of M. tuberculosis.

Crystal Structure of the RNA Lariat Debranching Enzyme Dbr1 with Hydrolyzed Phosphorothioate RNA Product.

The RNA lariat debranching enzyme is the sole enzyme responsible for hydrolyzing the 2'-5' phosphodiester bond in RNA lariats produced by the spliceosome. Here, we test the ability of Dbr1 to hydrolyze branched RNAs (bRNAs) that contain a 2'-5'-phosphorothioate linkage, a modification commonly used to resist degradation. We attempted to cocrystallize a phosphorothioate-branched RNA (PS-bRNA) with wild-type Entamoeba histolytica Dbr1 (EhDbr1) but observed in-crystal hydrolysis of the phosphorothioate bond. The crystal structure revealed EhDbr1 in a product-bound state, with the hydrolyzed 2'-5' fragment of the PS-bRNA mimicking the binding mode of the native bRNA substrate. These findings suggest that product inhibition may contribute to the kinetic mechanism of Dbr1. We show that Dbr1 enzymes cleave phosphorothioate linkages at rates ∼10,000-fold more slowly than native phosphate linkages. This new product-bound crystal structure offers atomic details, which can aid inhibitor design. Dbr1 inhibitors could be therapeutic or investigative compounds for human diseases such as human immunodeficiency virus (HIV), amyotrophic lateral sclerosis (ALS), cancer, and viral encephalitis.

A Combined Computational–Experimental Study on the Substrate Binding and Reaction Mechanism of Salicylic Acid Decarboxylase

Salicylic acid decarboxylase (SDC) from the amidohydrolase superfamily (AHS) catalyzes the reversible decarboxylation of salicylic acid to form phenol. In this study, the substrate binding mode and reaction mechanism of SDC were investigated using computational and crystallographic methods. Quantum chemical calculations show that the enzyme follows the general mechanism of AHS decarboxylases. Namely, the reaction begins with proton transfer from a metal-coordinated aspartic acid residue (Asp298 in SDC) to the C1 of salicylic acid, which is followed by the C–C bond cleavage, to generate the phenol product and release CO2. Interestingly, the calculations show that SDC is a Mg-dependent enzyme rather than the previously proposed Zn-dependent, and the substrate is shown to be bidentately coordinated to the metal center in the catalysis, which is also different from the previous proposal. These predictions are corroborated by the crystal structure of SDC solved in complex with the substrate analogue 2-nitrophenol. The mechanistic insights into SDC in the present study provide important information for the rational design of the enzyme.

Insights into the Structure of the Highly Glycosylated Ffase from Rhodotorula dairenensis Enhance Its Biotechnological Potential

Rhodotorula dairenensis β-fructofuranosidase is a highly glycosylated enzyme with broad substrate specificity that catalyzes the synthesis of 6-kestose and a mixture of the three series of fructooligosaccharides (FOS), fructosylating a variety of carbohydrates and other molecules as alditols. We report here its three-dimensional structure, showing the expected bimodular arrangement and also a unique long elongation at its N-terminus containing extensive O-glycosylation sites that form a peculiar arrangement with a protruding loop within the dimer. This region is not required for activity but could provide a molecular tool to target the dimeric protein to its receptor cellular compartment in the yeast. A truncated inactivated form was used to obtain complexes with fructose, sucrose and raffinose, and a Bis-Tris molecule was trapped, mimicking a putative acceptor substrate. The crystal structure of the complexes reveals the major traits of the active site, with Asn387 controlling the substrate binding mode. Relevant residues were selected for mutagenesis, the variants being biochemically characterized through their hydrolytic and transfructosylating activity. All changes decrease the hydrolytic efficiency against sucrose, proving their key role in the activity. Moreover, some of the generated variants exhibit redesigned transfructosylating specificity, which may be used for biotechnological purposes to produce novel fructosyl-derivatives.

The protein conformational basis of isoflavone biosynthesis

Isoflavonoids play important roles in plant defense and also exhibit a range of mammalian health-promoting activities. Their biosynthesis is initiated by two enzymes with unusual catalytic activities; 2-hydroxyisoflavanone synthase (2-HIS), a membrane-bound cytochrome P450 catalyzing a coupled aryl-ring migration and hydroxylation, and 2-hydroxyisoflavanone dehydratase (2-HID), a member of a large carboxylesterase family that paradoxically catalyzes dehydration of 2-hydroxyisoflavanones to isoflavone. Here we report the crystal structures of 2-HIS from Medicago truncatula and 2-HID from Pueraria lobata. The 2-HIS structure reveals a unique cytochrome P450 conformation and heme and substrate binding mode that facilitate the coupled aryl-ring migration and hydroxylation reactions. The 2-HID structure reveals the active site architecture and putative catalytic residues for the dual dehydratase and carboxylesterase activities. Mutagenesis studies revealed key residues involved in substrate binding and specificity. Understanding the structural basis of isoflavone biosynthesis will facilitate the engineering of new bioactive isoflavonoids.

Enabling Aromatic Hydroxylation in a Cytochrome P450 Monooxygenase Enzyme through Protein Engineering.

The cytochrome P450 (CYP) family of heme monooxygenases catalyse the selective oxidation of C-H bonds under ambient conditions. The CYP199A4 enzyme from Rhodopseudomonas palustris catalyses aliphatic oxidation of 4-cyclohexylbenzoic acid but not the aromatic oxidation of 4-phenylbenzoic acid, due to the distinct mechanisms of aliphatic and aromatic oxidation. The aromatic substrates 4-benzyl-, 4-phenoxy- and 4-benzoyl-benzoic acid and methoxy-substituted phenylbenzoic acids were assessed to see if they could achieve an orientation more amenable to aromatic oxidation. CYP199A4 could catalyse the efficient benzylic oxidation of 4-benzylbenzoic acid. The methoxy-substituted phenylbenzoic acids were oxidatively demethylated with low activity. However, no aromatic oxidation was observed with any of these substrates. Crystal structures of CYP199A4 with 4-(3'-methoxyphenyl)benzoic acid demonstrated that the substrate binding mode was like that of 4-phenylbenzoic acid. 4-Phenoxy- and 4-benzoyl-benzoic acid bound with the ether or ketone oxygen atom hydrogen-bonded to the heme aqua ligand. We also investigated whether the substitution of phenylalanine residues in the active site could permit aromatic hydroxylation. Mutagenesis of the F298 residue to a valine did not significantly alter the substrate binding position or enable the aromatic oxidation of 4-phenylbenzoic acid; however the F182L mutant was able to catalyse 4-phenylbenzoic acid oxidation generating 2'-hydroxy-, 3'-hydroxy- and 4'-hydroxy metabolites in a 83 : 9 : 8 ratio, respectively. Molecular dynamics simulations, in which the distance and angle of attack were considered, demonstrated that in the F182L variant, in contrast to the wild-type enzyme, the phenyl ring of 4-phenylbenzoic acid attained a productive geometry for aromatic oxidation to occur.