Free energy calculations give insight into the stereoselective hydroxylation of α-ionones by engineered cytochrome P450 BM3 mutants

Mechanisms of Electron Transfer Rate Modulations in Cytochrome P450 BM3

Bacterial cytochromes P450 BM3 (CYP450 BM3) catalyze reactions of industrial importance. Despite many successful biotransformations, robust (re)design for novel applications remains challenging. Rational design and evolutionary approaches are not always successful, highlighting a lack of complete understanding of the mechanisms of electron transfer (ET) modulations. Thus, the full potential of CYP450 reactions remains under-exploited. In this work, we report the first molecular dynamics (MD)-based explicit prediction of BM3 ET parameters (reorganization energies; λ and ET free energies; ΔG°), and log ET rates (log k_ET) using the Marcus theory. Overall, the calculated ET rates for the BM3 wild-type (WT), mutants (F393 and L86), ligand-bound state, and ion concentrations agree well with experimental data. In ligand-free (LF) BM3, mutations modulate k_ET via ET ΔG°. Simulations show that the experimental ET rate enhancement is due to increased driving force (more negative ΔG°) upon ligation. This increase is related to the protein reorganization required to accommodate the ligand in the binding pocket rather than binding interactions with the ligand. Our methodology (CYPWare 1.0) automates all the stages of the MD simulation step-up, energy calculations, and estimation of ET parameters. CYPWare 1.0 and this work thus represent an important advancement in the CYP450 ET rate predictions, which has the potential to guide the redesign of ET enzymes. This program and a Web tool are available on GitHub for academic research.

Insights into regioselective metabolism of mefenamic acid by cytochrome P450 BM3 mutants through crystallography, docking, molecular dynamics, and free energy calculations

Cytochrome P450 BM3 (CYP102A1) mutant M11 is able to metabolize a wide range of drugs and drug-like compounds. Among these, M11 was recently found to be able to catalyze formation of human metabolites of mefenamic acid and other nonsteroidal anti-inflammatory drugs (NSAIDs). Interestingly, single active-site mutations such as V87I were reported to invert regioselectivity in NSAID hydroxylation. In this work, we combine crystallography and molecular simulation to study the effect of single mutations on binding and regioselective metabolism of mefenamic acid by M11 mutants. The heme domain of the protein mutant M11 was expressed, purified, and crystallized, and its X-ray structure was used as template for modeling. A multistep approach was used that combines molecular docking, molecular dynamics (MD) simulation, and binding free-energy calculations to address protein flexibility. In this way, preferred binding modes that are consistent with oxidation at the experimentally observed sites of metabolism (SOMs) were identified. Whereas docking could not be used to retrospectively predict experimental trends in regioselectivity, we were able to rank binding modes in line with the preferred SOMs of mefenamic acid by M11 and its mutants by including protein flexibility and dynamics in free-energy computation. In addition, we could obtain structural insights into the change in regioselectivity of mefenamic acid hydroxylation due to single active-site mutations. Our findings confirm that use of MD and binding free-energy calculation is useful for studying biocatalysis in those cases in which enzyme binding is a critical event in determining the selective metabolism of a substrate.

Selective oxidation of carotenoid-derived aroma compounds by CYP260B1 and CYP267B1 from Sorangium cellulosum So ce56

Due to their bioactive properties as well as their application as precursors in chemical synthesis, hydroxylated isoprenoids and norisoprenoids are very valuable compounds. The efficient hydroxylation of such compounds remains a challenge in organic chemistry caused by the formation of a variety of side products and lack of overall regio- and stereoselectivity. In contrast, cytochromes P450 are known for their selective oxidation under mild conditions. Here, we demonstrate for the first time the ability of myxobacterial CYP260B1 and CYP267B1 from Sorangium cellulosum So ce56 to oxidize such carotenoid-derived aroma compounds. A focused library of 14 substrates such as ionones, damascones, as well as some of their isomers and derivatives was screened in vitro. Both P450s were capable of an efficient oxidation of all tested compounds. CYP260B1-dependent conversions mainly formed multiple products, whereas conversions by CYP267B1 resulted predominantly in a single product. To identify the main products by NMR spectroscopy, an Escherichia coli-based whole-cell system was used. CYP267B1 showed a hydroxylase activity towards the formation of allylic alcohols. Likewise, CYP260B1 performed the allylic hydroxylation of β-damascone [(E)-1-(2,6,6-trimethylcyclohex-1-enyl)but-2-en-1-one] and δ-damascone [(E)-1-(2,6,6-trimethylcyclohex-3-enyl)but-2-en-1-one]. Moreover, CYP260B1 showed an epoxidase activity towards β-ionone [(E)-4-(2,6,6-trimethylcyclohex-1-enyl)but-3-en-2-one] as well as the methyl-substituted α-ionone derivatives raldeine [(E)-1-(2,6,6-trimethylcyclohex-2-enyl)pent-1-en-3-one] and isoraldeine [(E)-3-methyl-4-(2,6,6-trimethylcyclohex-2-enyl)but-3-en-2-one]. In addition, to known products, also novel products such as 2-OH-δ-damascone [(E)-1-(5-hydroxy-2,6,6-trimethylcyclohex-3-enyl)but-2-en-1-one], 3-OH-allyl-α-ionone [(E)-1-(4-hydroxy-2,6,6-trimethylcyclohex-2-enyl)hepta-1,6-dien-3-one], and 4-OH-allyl-β-ionone [(E)-1-(3-hydroxy-2,6,6-trimethylcyclohex-1-enyl)hepta-1,6-dien-3-one] were identified during our studies.

Prediction of cytochrome P450 mediated metabolism

Cytochrome P450 enzymes (CYPs) form one of the most important enzyme families involved in the metabolism of xenobiotics. CYPs comprise many isoforms, which catalyze a wide variety of reactions, and potentially, a large number of different metabolites can be formed. However, it is often hard to rationalize what metabolites these enzymes generate. In recent years, many different in silico approaches have been developed to predict binding or regioselective product formation for the different CYP isoforms. These comprise ligand-based methods that are trained on experimental CYP data and structure-based methods that consider how the substrate is oriented in the active site or/and how reactive the part of the substrate that is accessible to the heme group is. We will review key aspects for various approaches that are available to predict binding and site of metabolism (SOM), what outcome can be expected from the predictions, and how they could potentially be improved.

Predicting drug metabolism: experiment and/or computation?

Drug metabolism can produce metabolites with physicochemical and pharmacological properties that differ substantially from those of the parent drug, and consequently has important implications for both drug safety and efficacy. To reduce the risk of costly clinical-stage attrition due to the metabolic characteristics of drug candidates, there is a need for efficient and reliable ways to predict drug metabolism in vitro, in silico and in vivo. In this Perspective, we provide an overview of the state of the art of experimental and computational approaches for investigating drug metabolism. We highlight the scope and limitations of these methods, and indicate strategies to harvest the synergies that result from combining measurement and prediction of drug metabolism.

Pyranose dehydrogenase ligand promiscuity: a generalized approach to simulate monosaccharide solvation, binding, and product formation

The flavoenzyme pyranose dehydrogenase (PDH) from the litter decomposing fungus Agaricus meleagris oxidizes many different carbohydrates occurring during lignin degradation. This promiscuous substrate specificity makes PDH a promising catalyst for bioelectrochemical applications. A generalized approach to simulate all 32 possible aldohexopyranoses in the course of one or a few molecular dynamics (MD) simulations is reported. Free energy calculations according to the one-step perturbation (OSP) method revealed the solvation free energies (ΔG_solv) of all 32 aldohexopyranoses in water, which have not yet been reported in the literature. The free energy difference between β- and α-anomers (ΔG_β-α) of all d-stereoisomers in water were compared to experimental values with a good agreement. Moreover, the free-energy differences (ΔG) of the 32 stereoisomers bound to PDH in two different poses were calculated from MD simulations. The relative binding free energies (ΔΔG_bind) were calculated and, where available, compared to experimental values, approximated from K_m values. The agreement was very good for one of the poses, in which the sugars are positioned in the active site for oxidation at C1 or C2. Distance analysis between hydrogens of the monosaccharide and the reactive N5-atom of the flavin adenine dinucleotide (FAD) revealed that oxidation is possible at HC1 or HC2 for pose A, and at HC3 or HC4 for pose B. Experimentally detected oxidation products could be rationalized for the majority of monosaccharides by combining ΔΔG_bind and a reweighted distance analysis. Furthermore, several oxidation products were predicted for sugars that have not yet been tested experimentally, directing further analyses. This study rationalizes the relationship between binding free energies and substrate promiscuity in PDH, providing novel insights for its applicability in bioelectrochemistry. The results suggest that a similar approach could be applied to study promiscuity of other enzymes.

Generally, enzymes are perceived as being specific for both their substrates and the reaction they catalyze. This standard paradigm started to shift and currently enzyme promiscuity towards various substrates is perceived rather as the rule than the exception. Enzyme promiscuity seems to be vital for proteins to acquire new functions, and therefore for evolution itself. The driving forces for promiscuity are manifold and consequently challenging to study. Binding free energies, which can be calculated from computer simulations, represent a convenient measure for them. Here, we investigate the binding free energies between the enzyme pyranose dehydrogenase (PDH) and a sugar-substrate by computational means. PDH has an extraordinarily promiscuous substrate-specificity, making it interesting for e.g. bioelectrochemical applications. By introducing modifications to the sugar-structure used for the molecular dynamics simulations, we could simultaneously study all 32 possible aldohexopyranoses from a single simulation. This saves costly computational resources and time for setting up and analyzing the simulations. We could nicely reproduce experimental results and predict so far undetected sugar-oxidation products, directing further experiments. This study gives novel insights into PDH's substrate promiscuity and the enzyme's applicability. A similar approach could be applied to study the promiscuity of other enzymes.

The five-to-six-coordination transition of ferric human serum heme-albumin is allosterically-modulated by ibuprofen and warfarin: a combined XAS and MD study

Human serum albumin (HSA) is involved physiologically in heme scavenging; in turn, heme-albumin (HSA-heme-Fe) displays globin-like properties. Here, the allosteric effect of ibuprofen and warfarin on the local atomic structure around the ferric heme-Fe (heme-Fe(III)) atom of HSA-heme-Fe (HSA-heme-Fe(III)) has been probed by Fe-K edge X-ray absorption spectroscopy (XAS). The quantitative analysis of the Fe-K edge extended X-ray absorption fine structure (EXAFS) signals and modeling of the near edge (XANES) spectral features demonstrated that warfarin and ibuprofen binding modify the local structure of the heme-Fe(III). Combined XAS data analysis and targeted molecular dynamics (MD) simulations provided atomic resolution insights of protein structural rearrangements required to accommodate the heme-Fe(III) upon ibuprofen and warfarin binding. In the absence of drugs, the heme-Fe(III) atom is penta-coordinated having distorted 4+1 configuration made by the nitrogen atoms of the porphyrin ring and the oxygen phenoxy atom of the Tyr161 residue. MD simulations show that ibuprofen and warfarin association to the secondary fatty acid (FA) binding site 2 (FA2) induces a reorientation of domain I of HSA-heme-Fe(III), this leads to the redirection of the His146 residue providing an additional bond to the heme-Fe(III) atom, providing the 5+1 configuration. The comparison of Fe-K edge XANES spectra calculated using MD structures with those obtained experimentally confirms the reliability of the proposed structural model. As a whole, combining XAS and MD simulations it has been possible to provide a reliable model of the heme-Fe(III) atom coordination state and to understand the complex allosteric transition occurring in HSA-heme-Fe(III) upon ibuprofen and warfarin binding.

Fast, metadynamics-based method for prediction of the stereochemistry-dependent relative free energies of ligand-receptor interactions

The computational approach applicable for the molecular dynamics (MD)-based techniques is proposed to predict the ligand-protein binding affinities dependent on the ligand stereochemistry. All possible stereoconfigurations are expressed in terms of one set of force-field parameters [stereoconfiguration-independent potential (SIP)], which allows for calculating all relative free energies by only single simulation. SIP can be used for studying diverse, stereoconfiguration-dependent phenomena by means of various computational techniques of enhanced sampling. The method has been successfully tested on the β2-adrenergic receptor (β2-AR) binding the four fenoterol stereoisomers by both metadynamics simulations and replica-exchange MD. Both the methods gave very similar results, fully confirming the presence of stereoselective effects in the fenoterol-β2-AR interactions. However, the metadynamics-based approach offered much better efficiency of sampling which allows for significant reduction of the unphysical region in SIP.

Entropic and enthalpic contributions to stereospecific ligand binding from enhanced sampling methods

The stereoselective binding of R- and S-propranolol to the metabolic enzyme cytochrome P450 2D6 and its mutant F483A was studied using various computational approaches. Previously reported free-energy differences from Hamiltonian replica exchange simulations, combined with thermodynamic integration, are compared to the one-step perturbation approach, combined with local-elevation enhanced sampling, and an excellent agreement between methods was obtained. Further, the free-energy differences are interpreted in terms of enthalpic and entropic contributions where it is shown that exactly compensating contributions obscure a molecular interpretation of differences in the affinity while various reduced terms allow a more detailed analysis, which agree with heuristic observations on the interactions.

Influence of variation of a side chain on the folding equilibrium of a β-peptide: limitations of one-step perturbation

In a recent study (Lin et al., Helv. Chim. Acta 2011, 94, 597), the one-step perturbation method was applied to tackle a challenging computational problem, that is, the calculation of the folding free enthalpies ΔGF,U of six hepta-β-peptides with different, Ala, Val, Leu, Ile, Ser, or Thr, side chains in the fifth residue. The ΔGF,U values obtained using one-step perturbation based on a single molecular dynamics simulation of a judiciously chosen reference state with soft-core atoms in the side chain of the fifth residue showed an overall accuracy of about kB T for the four peptides with nonpolar side chains, but twice as large deviations were observed for the peptides with polar side chains. Here, alternative reference-state Hamiltonians that better cover the conformational space relevant to these peptides are investigated, and post simulation rotational sampling of the χ1 and χ2 torsional angles of the fifth residue is carried out to sample different orientations of the side chain. A reference state with rather soft atoms yields accurate ΔGF,U values for the peptides with the Ser and Thr side chains, but it failed to correctly predict the folding free enthalpy for one peptide with a nonpolar side chain, that is, Leu. Based on the results and those of earlier studies, possible ways to improve the accuracy of the efficient one-step perturbation technique to compute free enthalpies of folding are discussed.