Molecular Dynamics and QM/MM Calculations Predict the Substrate-Induced Gating of Cytochrome P450 BM3 and the Regio- and Stereoselectivity of Fatty Acid Hydroxylation

A promiscuous cytochrome P450 aromatic O-demethylase for lignin bioconversion

Microbial aromatic catabolism offers a promising approach to convert lignin, a vast source of renewable carbon, into useful products. Aryl-O-demethylation is an essential biochemical reaction to ultimately catabolize coniferyl and sinapyl lignin-derived aromatic compounds, and is often a key bottleneck for both native and engineered bioconversion pathways. Here, we report the comprehensive characterization of a promiscuous P450 aryl-O-demethylase, consisting of a cytochrome P450 protein from the family CYP255A (GcoA) and a three-domain reductase (GcoB) that together represent a new two-component P450 class. Though originally described as converting guaiacol to catechol, we show that this system efficiently demethylates both guaiacol and an unexpectedly wide variety of lignin-relevant monomers. Structural, biochemical, and computational studies of this novel two-component system elucidate the mechanism of its broad substrate specificity, presenting it as a new tool for a critical step in biological lignin conversion.

Simulation-Guided Design of Cytochrome P450 for Chemo- and Regioselective Macrocyclic Oxidation

Engineering high chemo-, regio-, and stereoselectivity is a prerequisite for enzyme usage in organic synthesis. Cytochromes P450 can oxidize a broad range of substrates, including macrocycles, which are becoming popular scaffolds for therapeutic agents. However, a large conformational space explored by macrocycles not only reduces the selectivity of oxidation but also impairs computational enzyme design strategies based on docking and molecular dynamics (MD) simulations. We present a novel design workflow that uses enhanced-sampling Hamiltonian replica exchange (HREX) MD and focuses on quantifying the substrate binding for suggesting the mutations to be made. This computational approach is applied to P450 BM3 with the aim to shift regioselectively toward one of the numerous possible positions during β-cembrenediol oxidation. The predictions are experimentally tested and the resulting product distributions validate our design strategy, as single mutations led up to 5-fold regioselectivity increases. We thus conclude that the HREX-MD-based workflow is a promising tool for the identification of positions for mutagenesis aiming at P450 enzymes with improved regioselectivity.

Exploring the structure characteristics and major channels of cytochrome P450 2A6, 2A13, and 2E1 with pilocarpine

The majority of cytochromes P450 play a critical role in metabolism of endogenous and exogenous substrates, some of its products are carcinogens. Therefore, inhibition of P450 enzymes activity can promote the detoxification and elimination of chemical carcinogens. In this study, molecular dynamics (MD) simulations and adaptive steered molecular dynamics (ASMD) simulations were performed to explore the structure features and channel dynamics of three P450 isoforms 2A6, 2A13, and 2E1 bound with the common inhibitor pilocarpine. The binding free energy results combined with the PMF calculations give a reasonable ranking of binding affinity, which are consistent with the experimental data. Our results uncover how a sequence divergence of different CYP2 enzymes causes individual variations in major channel selections. On the basis of channel bottleneck and energy decomposition analysis, we propose a gating mechanism of their respective major channels in three enzymes, which may be attributed to a reversal of Phe209 in CYP2A6/2A13, as well as the rotation of Phe116 and Phe298 in CYP2E1. The hydrophobic residues not only make strong hydrophobic interactions with inhibitor, but also act as gatekeeper to regulate the opening of channel. The present study provides important insights into the structure-function relationships of three cytochrome P450s and the molecular basis for development of potent inhibitors.

Density Functional Theory Study on the Demethylation Reaction between Methylamine, Dimethylamine, Trimethylamine, and Tamoxifen Catalyzed by a Fe(IV)-Oxo Porphyrin Complex

In this work, we studied computationally the N-demethylation reaction of methylamine, dimethylamine, and trimethylamine as archetypal examples of primary, secondary, and tertiary amines catalyzed by high-field low-spin Fe-containing enzymes such as cytochromes P450. Using DFT calculations, we found that the expected C-H hydroxylation process was achieved for trimethylamine. When dimethylamine and methylamine were studied, two different reaction mechanisms (C-H hydroxylation and a double hydrogen atom transfer) were computed to be energetically accessible and both are equally preferred. Both processes led to the formation of formaldehyde and the N-demethylated substrate. Finally, as an illustrative example, the relative contribution of the three primary oxidation routes of tamoxifen was rationalized through energetic barriers obtained from density functional calculations and docking experiments involving CYP3A4 and CYP2D6 isoforms. We found that the N-demethylation process was the intrinsically favored one, whereas other oxidation reactions required most likely preorganization imposed by the residues close to the active sites.

Choreography of the Reductase and P450BM3 Domains Toward Electron Transfer Is Instigated by the Substrate

The driving force for the electron transfer (ET) step in the catalytic cycle of cytochrome P450_BM3 is investigated using three sets of 1 μs molecular dynamic simulations for the resting state of P450 in complex with the flavin (FMN) in the reductase domain. These sets involve the following: (i) substrate-free (SF), (ii) substrate (N-palmitoyl glycine, i.e., NPG)-bound (SB), and (iii) SB with the semiquinone radical anion (SQ^-) of FMN. Starting from the X-ray structure of the SF heme domain, we observe that the α1-helix (of the reductase) and the C-helix (of the heme) undergo reorientation to a parallel orientation, which is the thermochemically stable form. The reorientation of the helices pushes away the FMN to a distance of 18.4 Å from the heme's center. When the substrate binds it causes the I-helix of the heme domain to kink and push the C-helix toward the α1-helix, thereby locking the latter two into a stabilized perpendicular conformation, wherein the FMN-heme distance is 12 Å. The distance drops further in the SQ^- form, and upon QM/MM geometry optimization the two moieties approach 8.8 Å, which enhances the ET rate (by 10⁴-10⁶ fold) to the heme's Fe³⁺ ion. These motions are driven by hydrogen bond strengthening between the C- and the α1-helices. Finally, substrate binding leads to formation of an organized water chain connecting the FMN and heme moieties. The water channel assists the ET and couples it to the proton transfer steps that should activate O₂ and create the oxo-iron active species.

Combined Quantum Mechanics and Molecular Mechanics Studies of Enzymatic Reaction Mechanisms

The combined quantum mechanics/molecular mechanics (QM/MM) methods have become a valuable tool in computational biochemistry and received versatile applications for studying the reaction mechanisms of enzymes. The approach combines the calculations of the electronic structure of the active site by QM, with modeling of the protein environment using MM force field, which allows the long-range electrostatics and steric effects on the enzyme reactivity to be accounted for. In this review, we review some key theoretical and computational aspects of the method and we also present some applications to particular enzymatic reactions such as tryptophan-7-halogenase, cyclooxygenase-1, and the epidermal growth factor receptor.

Mechanistic insight into the C7-selective C–H functionalization of N-acyl indole catalyzed by a rhodium complex: a theoretical study

The role of the additive AgNTf₂ and the origins of the reaction are clarified through our calculations.

Oxidation reactivity of 1,2-bis(2,4,6-tribromophenoxy)ethane (BTBPE) by Compound I model of cytochrome P450s

Alternative brominated flame retardants (BFRs) have become prevalent as a consequence of restrictions on the use of polybrominated diphenyl ethers (PBDEs). For risk assessment of these alternatives, knowledge of their metabolism via cytochrome P450 enzymes is needed. We have previously proved that density functional theory (DFT) is able to predict the metabolism of PBDEs by revealing the molecular mechanisms. In the current study, the reactivity of 1,2-bis(2,4,6-tribromophenoxy)ethane and structurally similar chemicals with the Compound I model representing the active site of P450 enzymes was investigated. The DFT calculations delineated reaction pathways which lead to reasonable explanations for products that were detected by wet experiments, meanwhile intermediates which cannot be determined were also proposed. Results showed that alkyl hydrogen abstraction will lead to bis(2,4,6-tribromophenoxy)ethanol, which may undergo hydrolysis yielding 2,4,6-tribromophenol, a neurotoxic compound. In addition, a general pattern of oxidation reactivity regarding the 2,4,6-tribromophenyl moiety was observed among several model compounds. Our study has provided insights for convenient evaluation of the metabolism of other structurally similar BFRs.

Conformational Heterogeneity and the Affinity of Substrate Molecular Recognition by Cytochrome P450cam

The broad and variable substrate specificity of cytochrome P450 enzymes makes them a model system for studying the determinants of protein molecular recognition. The archetypal cytochrome P450cam (P450cam) is a relatively specific P450, a feature once attributed to the high rigidity of its active site. However, increasingly studies have provided evidence of the importance of conformational changes to P450cam activity. Here we used infrared (IR) spectroscopy to investigate the molecular recognition of P450cam. Toward this goal, and to assess the influence of a hydrogen bond (H-bond) between active site residue Y96 and substrates, two variants in which Y96 is replaced by a cyanophenyl (Y96CNF) or phenyl (Y96F) group were characterized in complexes with the substrates camphor, isoborneol, and camphane. These combinations allow for a comparison of complexes in which the moieties on both the protein and substrate can serve as a H-bond donor, acceptor, or neither. The IR spectra of heme-bound CO and the site-specifically incorporated CN of Y96CNF were analyzed to characterize the number and nature of environments in each protein, both in the free and bound states. Although the IR spectra do not support the idea that protein-substrate H-bonding is central to P450cam recognition, the data altogether suggest that the differing conformational heterogeneity in the active site of the P450cam variants and changes in heterogeneity upon binding of different substrates likely contribute to their variable affinities via a conformational selection mechanism. This study further extends our understanding of the molecular recognition of archetypal P450cam and demonstrates the application of IR spectroscopy combined with selective protein modification to delineate protein-ligand interactions.

A redox-mediated Kemp eliminase

The acid/base-catalysed Kemp elimination of 5-nitro-benzisoxazole forming 2-cyano-4-nitrophenol has long served as a design platform of enzymes with non-natural reactions, providing new mechanistic insights in protein science. Here we describe an alternative concept based on redox catalysis by P450-BM3, leading to the same Kemp product via a fundamentally different mechanism. QM/MM computations show that it involves coordination of the substrate's N-atom to haem-Fe(II) with electron transfer and concomitant N-O heterolysis liberating an intermediate having a nitrogen radical moiety Fe(III)-N· and a phenoxyl anion. Product formation occurs by bond rotation and H-transfer. Two rationally chosen point mutations cause a notable increase in activity. The results shed light on the prevailing mechanistic uncertainties in human P450-catalysed metabolism of the immunomodulatory drug leflunomide, which likewise undergoes redox-mediated Kemp elimination by P450-BM3. Other isoxazole-based pharmaceuticals are probably also metabolized by a redox mechanism. Our work provides a basis for designing future artificial enzymes.

Structural-functional adaptations of porcine CYP1A1 to metabolize polychlorinated dibenzo-p-dioxins

Polychlorinated dibenzo-p-dioxins (PCDDs) are widespread by-products of human industrial activity. They accumulate in tissues of animals and humans, exerting numerous adverse effects on different systems. In living organisms, dioxins are metabolized by enzymes of the cytochrome P450 family, including CYP1A1. Particular dioxin congeners differ in their toxicity level and ability to undergo biodegradation. Since the molecular mechanisms underlying dioxin susceptibility or resistance to biodegradation are unknown, in the present study the molecular interactions between five selected dioxins and porcine CYP1A1 protein were investigated. It was found that the ability of a dioxin to undergo CYP1A1-mediated degradation is associated mainly with the number and position of chlorine atoms in the dioxin molecule. Among all examined congeners, 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) demonstrated the highest affinity to CYP1A1 and, at the same time, the greatest distance to the active site of the enzyme. Interestingly, in contrast to other dioxins, the binding of the TCDD molecule to the porcine CYP1A1 active site resulted in a rapid and continuous closure of substrate channels. All the information may help to explain the extended half-life of TCDD in living organisms as well as its high toxicity.

Control of stereoselectivity of benzylic hydroxylation catalysed by wild-type cytochrome P450BM3 using decoy molecules

The benzylic hydroxylation of non-native substrates was catalysed by cytochrome P450BM3, wherein “decoy molecules” controlled the stereoselectivity of the reactions.

Cleavage mechanism of the aliphatic C–C bond catalyzed by 2,4′-dihydroxyacetophenone dioxygenase from Alcaligenes sp. 4HAP: a QM/MM study

Calculations suggest that the reactant complex may firstly undergo a triplet–quintet crossing to initiate the reaction and then the subsequent chemistry occurs on the multiple-states surfaces. The key C–C bond cleavage is accompanied by an insertion reaction of oxygen radical.

Predicting Regioselectivity and Lability of Cytochrome P450 Metabolism Using Quantum Mechanical Simulations

We describe methods for predicting cytochrome P450 (CYP) metabolism incorporating both pathway-specific reactivity and isoform-specific accessibility considerations. Semiempirical quantum mechanical (QM) simulations, parametrized using experimental data and ab initio calculations, estimate the reactivity of each potential site of metabolism (SOM) in the context of the whole molecule. Ligand-based models, trained using high-quality regioselectivity data, correct for orientation and steric effects of the different CYP isoform binding pockets. The resulting models identify a SOM in the top 2 predictions for between 82% and 91% of compounds in independent test sets across seven CYP isoforms. In addition to predicting the relative proportion of metabolite formation at each site, these methods estimate the activation energy at each site, from which additional information can be derived regarding their lability in absolute terms. We illustrate how this can guide the design of compounds to overcome issues with rapid CYP metabolism.

Mechanism for the enhanced reactivity of 4-mercaptoprolyl thioesters in native chemical ligation

Ring-strain-precluded strategy benefiting from entropy effects and n → π* orbital interaction, enhances the reactivity of C-terminal prolyl thioesters in NCL.