The role of phenylalanine 483 in cytochrome P450 2D6 is strongly substrate dependent

Removal of BPA by Pseudomonas asiatica P1: Synergistic response mechanism of toxicity resistance and biodegradation

Bisphenol A (BPA) is a globally concerning toxic pollutant, and microbial degradation is considered an effective method to treat BPA contamination. However, the inherent microbial toxicity of BPA is often overlooked, particularly the microbial mechanisms of resistance and detoxification against BPA. This study found that under the toxic stress of BPA, cbb3-type cytochrome c oxidase (cbb3-Cox) in the cells of Pseudomonas asiatica P1 (P. asiatica P1) was the first to resist the toxicity. Genes such as ccoNOQPG showed significant upregulation with an average log2FC value of 3.56. Subsequently, genes that are related to metal ion binding, transport, and DNA repair were upregulated in the middle to later phase, which enhanced the metabolic functions of the strains and induced strain mutations to assist P. asiatica P1 in resisting the BPA toxicity. Meanwhile, three potential BPA degradation genes were identified, among which sdrP1 was crucial to the BPA degradation and detoxification. After genetic recombination, sdrP1 achieved a degradation rate of 92.52 % for BPA. Furthermore, through various methods such as alkyl interactions, sdrP1 exhibited oxidation and demethylation to form lower toxic intermediate products and complete the biological detoxification of BPA. This study provides a systematic analysis of the toxicity resistance, biodegradation, and detoxification processes in bacterial BPA removal, refines the mechanism of BPA biodegradation and contributes to a more comprehensive and systematic understanding of the overall process of microbial removal of toxic pollutants.

Molecular Mechanisms Underlying Single Nucleotide Polymorphism-Induced Reactivity Decrease in CYP2D6

Cytochrome P450 2D6 (CYP2D6) is one of the most important enzymes involved in drug metabolism. Genetic polymorphism can influence drug metabolism by CYP2D6 such that a therapy is seriously affected by under- or overdosing of drugs. However, a general explanation at the atomistic level for poor activity is missing so far. Here we show for the 20 most common single nucleotide polymorphisms (SNPs) of CYP2D6 that poor metabolism is driven by four mechanisms. We found in extensive all-atom molecular dynamics simulations that the rigidity of the I-helix (central helix), distance between central phenylalanines (stabilizing bound substrate), availability of basic residues on the surface of CYP2D6 (binding of cytochrome P450 reductase), and position of arginine 132 (electron transfer to heme) are essential for an extensive function of the enzyme. These results were applied to SNPs with unknown effects, and potential SNPs that may lead to poor drug metabolism were identified. The revealed molecular mechanisms might be important for other drug-metabolizing cytochrome P450 enzymes.

Differential Interactions of Selected Phytocannabinoids with Human CYP2D6 Polymorphisms

Cytochrome P450 2D6 (CYP2D6) is primarily expressed in the liver and in the central nervous system. It is known to be highly polymorphic in nature. It metabolizes several endogenous substrates such as anandamide (AEA). Concomitantly, it is involved in phase 1 metabolism of several antidepressants, antipsychotics, and other drugs. Research in the field of phytocannabinoids (pCBs) has recently accelerated owing to their legalization and increasing medicinal use for pain and inflammation. The primary component of cannabis is THC, which is well-known for its psychotropic effects. Since CYP2D6 is an important brain and liver P450 and is known to be inhibited by CBD, we investigated the interactions of four important highly prevalent CYP2D6 polymorphisms with selected phytocannabinoids (CBD, THC, CBDV, THCV, CBN, CBG, CBC, β-carophyllene) that are rapidly gaining popularity. We show that there is differential binding of CYP2D6*17 to pCBs as compared to WT CYP2D6. We also perform a more detailed comparison of WT and *17 CYP2D6, which reveals the possible regulation of AEA metabolism by CBD. Furthermore, we use molecular dynamics to delineate the mechanism of this binding, inhibition, and regulation. Taken together, we have found that the interactions of CYP2D6 with pCBs vary by polymorphism and by specific pCB class.

Tryptophan-75 Is a Low-Energy Channel-Gating Residue that Facilitates Substrate Egress/Access in Cytochrome P450 2D6

CYP2D6 is a major drug metabolizing enzyme with a buried active site. Channels leading to the active site from various enzyme surfaces are believed to facilitate ligand egress and access to the active site. The present study used molecular dynamics (MD) and in vitro studies with CYP2D6*1 and a Trp75-to-Ala mutant to examine channel gating in CYP2D6 by Trp75. MD simulations measured energy landscapes of Trp75 conformations and simulated substrate passage within channel 2b using bufuralol as a model substrate. Trp75 alternated between multiple stable states that supported substrate transport along channel 2b with low-energy barriers between states (∼ -1 kcal/mol). Trp75 conformations were stabilized primarily by hydrogen bonding between Trp75 and Glu222, Asn226, Ala225, or Gln72. Energy barriers were low between Trp75 conformations, allowing Trp75 to easily move between various conformations over time and to function in both binding to and moving substrates in the 2b channel of CYP2D6. Michaelis-Menten kinetic studies completed with purified enzyme in a reconstituted system showed overall reduced enzyme efficiency for metabolism of bufuralol and dextromethorphan by the Trp75Ala mutant compared with CYP2D6*1. In stopped-flow measurements, k off for dextromethorphan was decreased in the absence of Trp75. Our results support a role for Trp75 in substrate shuttling to the active site of CYP2D6. SIGNIFICANCE STATEMENT: Using combined molecular dynamics and in vitro assays, this study shows for the first time a role for Trp75 as a channel entrance gating residue in the mechanism of substrate binding/unbinding in CYP2D6. Energy landscapes derived from molecular dynamics were used to quantitate the strength of gating, and kinetics assays showed the impact on enzyme efficiency and k off of a Trp75Ala mutation.

Sympathomimetic amine compounds and hepatotoxicity: Not all are alike-Key distinctions noted in a short review

Sympathomimetic amine compounds are often pooled together and incorrectly assumed to be interchangeable with respect to potential adverse effects. A brief and specific review of sympathomimetic compounds and one instance (i.e., hepatotoxicity) where these compounds have been improperly grouped together is covered. A review of the proposed mechanisms through which known hepatotoxic sympathomimetic agents (e.g., 3,4-methylenedioxymethamphetamine or MDMA, methamphetamine and amphetamine) cause liver injury, along with a corresponding review of in vitro data, interventional data, animal model studies and observational data allow for a comparison/contrast of different agents and reveals a lack of potential toxicity for some agents (e.g., pseudoephedrine, phenylephrine, ephedrine, 1,3-dimethylamylamine, phentermine) in this broad category. Data show that compounds within the broad group of sympathomimetics display divergent pharmacological and toxicological profiles and can be clearly distinguished with respect to liver injury. These data serve as a reminder to clinicians and others, that even small structural differences between molecules can lead to drastically different pharmacological/toxicological profiles and that one should not assume that all sympathomimetic agents are hepatotoxic. Such assumptions could lead to diagnostic errors and incorrect or insufficient treatment.

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.

Double site saturation mutagenesis of the human cytochrome P450 2D6 results in regioselective steroid hydroxylation

The human cytochrome P450 2D6 (CYP2D6) is one of the major human drug metabolizing enzymes and acts preferably on substrates containing a basic nitrogen atom. Testosterone - just as other steroids - is an atypical substrate and only poorly metabolized by CYP2D6. The present study intended to investigate the influence of the two active site residues 216 and 483 on the capability of CYP2D6 to hydroxylate steroids such as for example testosterone. All 400 possible combinatorial mutations at these two positions have been generated and expressed individually in Pichia pastoris. Employing whole-cell biotransformations coupled with HPLC-MS analysis the testosterone hydroxylase activity and regioselectivity of every single CYP2D6 variant was determined. Covering the whole sequence space, CYP2D6 variants with improved activity and so far unknown regio-preference in testosterone hydroxylation were identified. Most intriguingly and in contrast to previous literature reports about mutein F483I, the mutation F483G led to preferred hydroxylation at the 2β-position, while the slow formation of 6β-hydroxytestosterone, the main product of wild-type CYP2D6, was further reduced. Two point mutations have already been sufficient to convert CYP2D6 into a steroid hydroxylase with the highest ever reported testosterone hydroxylation rate for this enzyme, which is of the same order of magnitude as for the conversion of the standard substrate bufuralol by wild-type CYP2D6. Furthermore, this study is also an example for efficient human CYP engineering in P. pastoris for biocatalytic applications and to study so far unknown pharmacokinetic effects of individual and combined mutations in these key enzymes of the human drug metabolism.

Rationalization of stereospecific binding of propranolol to cytochrome P450 2D6 by free energy calculations

Abstract

Cytochrome P450 2D6 is a major drug-metabolising enzyme with a wide substrate range. A single-point mutation introduced in this enzyme induces stereoselective binding of R and S-propranolol whereas the wild type has no preference. The system has previously been studied both experimentally and computationally (de Graaf et al. in Eur Biophys J 36:589–599, 2007a). The in silico study reported hysteresis and significant deviations from closure of thermodynamic cycles, probably because of lack of sampling. Here, we focus on the effect of prolonged simulation time and enhanced sampling methods, such as Hamiltonian replica exchange, to reduce these problems and to improve the precision of free energy calculations. Finally we rationalize the results at a molecular level and compare data with experimental findings and previously estimated free energies.

Graphical Abstract

Propranolol (PPD) binds stereospecifically to cytochrome P450 2D6, if the mutation F483A is introduced, by reducing binding affinity to R-PPD. This observation has previously been studied by molecular dynamics, by use of several different thermodynamic cycles. Previous computational results suffered from internal inconsistencies, probably because of insufficient sampling. Here we use prolonged simulation time and Hamiltonian replica exchange to overcome the inconsistencies, and give a possible rationalization of the phenomenon at a molecular level.

New insights into the structural characteristics and functional relevance of the human cytochrome P450 2D6 enzyme

To date, the crystal structures of at least 12 human CYPs (1A2, 2A6, 2A13, 2C8, 2C9, 2D6, 2E1, 2R1, 3A4, 7A1, 8A1, and 46A1) have been determined. CYP2D6 accounts for only a small percentage of all hepatic CYPs (< 2%), but it metabolizes approximately 25% of clinically used drugs with significant polymorphisms. CYP2D6 also metabolizes procarcinogens and neurotoxins, such as 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine, 1,2,3,4-tetrahydroquinoline, and indolealkylamines. Moreover, the enzyme utilizes hydroxytryptamines and neurosteroids as endogenous substrates. Typical CYP2D6 substrates are usually lipophilic bases with an aromatic ring and a nitrogen atom, which can be protonated at physiological pH. Substrate binding is generally followed by oxidation (5-7 A) from the proposed nitrogen-Asp301 interaction. A number of homology models have been constructed to explore the structural features of CYP2D6, while antibody studies also provide useful structural information. Site-directed mutagenesis studies have demonstrated that Glu216, Asp301, Phe120, Phe481, and Phe483 play important roles in determining the binding of ligands to CYP2D6. The structure of human CYP2D6 has been recently determined and shows the characteristic CYP fold observed for other members of the CYP superfamily. The lengths and orientations of the individual secondary structural elements in the CYP2D6 structure are similar to those seen in other human CYP2 members, such as CYP2C9 and 2C8. The 2D6 structure has a well-defined active-site cavity located above the heme group with a volume of approximately 540 A(3), which is larger than equivalent cavities in CYP2A6 (260 A(3)), 1A2 (375 A(3)), and 2E1 (190 A(3)), but smaller than those in CYP3A4 (1385 A(3)) and 2C8 (1438 A(3)). Further studies are required to delineate the molecular mechanisms involved in CYP2D6 ligand interactions and their implications for drug development and clinical practice.

Metabolism ofN-substituted 7-methoxy-4-(aminomethyl) -coumarins by cytochrome P450 2D6 mutants and the indication of additional substrate interaction points

Previous studies have shown the critical roles residues F120 and F483 play in the oxidative metabolism of 7-methoxy-4-(aminomethyl)-coumarin (MAMC) by cytochrome P450 2D6 (CYP2D6). In the present study, a series of N-alkyl-7-methoxy-4-(aminomethyl)-coumarins (MAMC analogues) were used as substrates for the F120A and F483A mutants in order to probe the CYP2D6 active site. The F120A and F483A mutants of CYP2D6 displayed significant activity towards the MAMC analogues. Automated docking studies of the MAMC analogues in a CYP2D6 homology model suggested a distal hydrophobic active site binding cleft for the substrate N-alkyl chains, consisting of the residues L213 and V308.

Impact of plasticity and flexibility on docking results for cytochrome P450 2D6: a combined approach of molecular dynamics and ligand docking

Cytochrome P450s (CYPs) exhibit a large plasticity and flexibility in the active site allowing for the binding of a large variety of substrates. The impact of plasticity and flexibility on ligand binding is investigated by docking 65 known CYP2D6 substrates to an ensemble of 2500 protein structures. The ensemble was generated by molecular dynamics simulations of CYP2D6 in complex with five representative substrates. The effect of induced fit, the conformation of Phe483, and thermal motion on the accuracy of site of metabolism (SOM) predictions is analyzed. For future predictions, the three most essential CYP2D6 structures were selected which are suitable for different kinds of ligands. We have developed a binary decision tree to decide which protein structure to dock the ligand into, such that each ligand needs to be docked only once, leading to successful SOM prediction in 80% of the substrates.

AmineDB: Large scale docking of amines with CYP2D6 and scoring for druglike properties—towards defining the scope of the chemical defense against foreign amines in humans

Organic amines are prevalent in nature and in drugs, especially the psychotherapeutic agents, and a major defense against potentially toxic amines is metabolism by CYP2D6. In order to understand better the constraints on the broad specificity of CYP2D6, 4207 amines were docked into the binding site of this enzyme. Docking poses were found predominantly with the positively charged amino groups closest to Asp301, with aromatic rings close to Phe120 and sometimes extending as far as Phe483. Organic amines that bind best to CYP2D6 tend to have larger molecular weights and logP values. Organic amines that score highly as being druglike, based on a Bayesian model constructed using a 5223-drug training set, are least likely to bind to CYP2D6. This correlation suggests that the set of known drugs, which have been largely designed or selected to avoid high affinity CYP binding, partially encodes the binding site preferences (or rather anti-preferences) of CYP2D6. Finally, in order to benchmark our docking and druglike scoring procedures, an analysis of psychotherapeutic agents is presented. All of these data, including the 4207 AM1-optimized ligand structures in proper ionization states, docking poses and scores, Druglike Scores and Lipinski properties, can be viewed from an online database, the AmineDB.

Analysis of CYP2D6 substrate interactions by computational methods

Cytochrome P450 CYP2D6 is involved in the oxidation of well over 150 drugs and, in general, those which contain a basic nitrogen atom in the molecule. To clarify how the residues of CYP2D6 are utilized for orientating a wide range of its specific substrates and distinguishing them from a variety of other organic compounds, docking studies by AutoDock and molecular dynamics (MD) simulations were conducted. Specific ligands were docked to both the homology model and crystal structures optimally to estimate the site of reaction on the ligand molecule and the binding energy for the complex, which were generally in good agreement with the experimental data. MD simulation for the CYP2D6-propranolol complex was then carried out to reveal the amino acid residues interacting with the substrate at the active site. Phe-120, Glu-216, Asp-301, and Phe-483 are identified as the substrate-binding residues in agreement with previously reported site-directed mutagenesis data and the crystal structure reported recently (PDB code: 2F9Q). As well as these residues, our theoretical prediction suggests that Phe-219 and Glu-222 are also important residues for mediating oxidation of substrates, especially propranolol.

The bacterial P450 BM3: a prototype for a biocatalyst with human P450 activities

The use of cytochrome P450 (P450 or CYP) enzymes as biocatalysts for the production of fine chemicals, including pharmaceuticals, has been of increasing interest, primarily owing to their catalytic diversity and broad substrate range. CYP102A1 (P450 BM3) from Bacillus megaterium integrates an entire monooxygenase system into one polypeptide and represents an appropriate prokaryotic model for industrial applications of mammalian P450 activities. CYP102A1 not only exhibits the highest catalytic activity ever detected in a P450 monooxygenase but also provides a potentially versatile biocatalyst for the production of human P450 metabolites. CYP102A1 can be further engineered to be a drug-metabolizing enzyme, making it a promising candidate to use as a biocatalyst in drug discovery and synthesis.

Free energies of binding of R- and S-propranolol to wild-type and F483A mutant cytochrome P450 2D6 from molecular dynamics simulations

Detailed molecular dynamics (MD) simulations have been performed to reproduce and rationalize the experimental finding that the F483A mutant of CYP2D6 has lower affinity for R-propranolol than for S-propranolol. Wild-type (WT) CYP2D6 does not show this stereospecificity. Four different approaches to calculate the free energy differences have been investigated and were compared to the experimental binding data. From the differences between calculations based on forward and backward processes and the closure of thermodynamic cycles, it was clear that not all simulations converged sufficiently. The approach that calculates the free energies of exchanging R-propranolol with S-propranolol in the F483A mutant relative to the exchange free energy in WT CYP2D6 accurately reproduced the experimental binding data. Careful inspection of the end-points of the MD simulations involved in this approach, allowed for a molecular interpretation of the observed differences.

Synergistic Use of Compound Properties and Docking Scores in Neural Network Modeling of CYP2D6 Binding: Predicting Affinity and Conformational Sampling

Cytochrome P450 2D6 (CYP2D6) is used to develop an approach for predicting affinity and relevant binding conformation(s) for highly flexible binding sites. The approach combines the use of docking scores and compound properties as attributes in building a neural network (NN) model. It begins by identifying segments of CYP2D6 that are important for binding specificity, based on structural variability among diverse CYP enzymes. A family of distinct, low-energy conformations of CYP2D6 are generated using simulated annealing (SA) and a collection of 82 compounds with known CYP2D6 affinities are docked. Interestingly, docking poses are observed on the backside of the heme as well as in the known active site. Docking scores for the active site binders, along with compound-specific attributes, are used to train a neural network model to properly bin compounds as strong binders, moderate binders, or nonbinders. Attribute selection is used to preselect the most important scores and compound-specific attributes for the model. A prediction accuracy of 85+/-6% is achieved. Dominant attributes include docking scores for three of the 20 conformations in the ensemble as well as the compound's formal charge, number of aromatic rings, and AlogP. Although compound properties were highly predictive attributes (12% improvement over baseline) in the NN-based prediction of CYP2D6 binders, their combined use with docking score attributes is synergistic (net increase of 23% above baseline). Beyond prediction of affinity, attribute selection provides a way to identify the most relevant protein conformation(s), in terms of binding competence. In the case of CYP2D6, three out of the ensemble of 20 SA-generated structures are found to be the most predictive for binding.

Catalytic site prediction and virtual screening of cytochrome P450 2D6 substrates by consideration of water and rescoring in automated docking

Automated docking strategies successfully applied to binding mode predictions of ligands in Cyt P450 crystal structures in an earlier study (de Graaf et al. J. Med. Chem. 2005, 7, 2308-2318) were used for the catalytic site prediction (CSP) of 65 substrates in a CYP2D6 homology model. The consideration of water molecules at predicted positions in the active site and the rescoring of pooled docking poses from four different docking programs (AutoDock, FlexX, GOLD-Goldscore, and GOLD-Chemscore) with the SCORE scoring function enabled the successful prediction of experimentally reported sites of catalysis of more than 80% of the substrates. Three docking algorithms (FlexX, GOLD-Goldscore, and GOLD-Chemscore) were subsequently used in combination with six scoring functions (Chemscore, DOCK, FlexX, GOLD, PMF, and SCORE) to assess the ability of docking-based virtual screening methods to prioritize known CYP2D6 substrates seeded into a drug-like chemical database (in the absence and presence of active-site water molecules). Finally, the optimal docking strategy in terms of virtual screening accuracy, GOLD-Chemscore with the consideration of active-site water (60% of known substrates recovered in the top 5% of the ranked drug-like database), was verified experimentally; it was successfully used to identify high-affinity CYP2D6 ligands among a larger proprietary database.

Crystal Structure of Human Cytochrome P450 2D6

Cytochrome P450 2D6 is a heme-containing enzyme that is responsible for the metabolism of at least 20% of known drugs. Substrates of 2D6 typically contain a basic nitrogen and a planar aromatic ring. The crystal structure of human 2D6 has been solved and refined to 3.0A resolution. The structure shows the characteristic P450 fold as seen in other members of the family, with the lengths and orientations of the individual secondary structural elements being very similar to those seen in 2C9. There are, however, several important differences, the most notable involving the F helix, the F-G loop, the B'helix, beta sheet 4, and part of beta sheet 1, all of which are situated on the distal face of the protein. The 2D6 structure has a well defined active site cavity above the heme group, containing many important residues that have been implicated in substrate recognition and binding, including Asp-301, Glu-216, Phe-483, and Phe-120. The crystal structure helps to explain how Asp-301, Glu-216, and Phe-483 can act as substrate binding residues and suggests that the role of Phe-120 is to control the orientation of the aromatic ring found in most substrates with respect to the heme. The structure has been compared with published homology models and has been used to explain much of the reported site-directed mutagenesis data and help understand the metabolism of several compounds.

Role of the conserved threonine 309 in mechanism of oxidation by cytochrome P450 2D6

Based on sequence alignments and homology modeling, threonine 309 in cytochrome P450 2D6 (CYP2D6) is proposed to be the conserved I-helix threonine, which is supposed to be involved in dioxygen activation by CYPs. The T309V mutant of CYP2D6 displayed a strong shift from O-dealkylation to N-dealkylation reactions in oxidation of dextromethorphan and 3,4-methylenedioxymethylamphetamine. This may be explained by an elevated ratio of hydroperoxo-iron to oxenoid-iron of the oxygenating species. In consistence, using cumene hydroperoxide, which directly forms the oxenoid-iron, the T309V mutant again selectively catalyzed the O-dealkylation reactions. The changed ratio of oxygenating species can also explain the decreased activity and changed regioselectivity that were observed in 7-methoxy-4-(aminomethyl)-coumarin and bufuralol oxidation, respectively, by the T309V mutant. Interestingly, the T309V mutant always showed a significantly increased, up to 75-fold, higher activity compared to that of the wild-type when using cumene hydroperoxide. These results indicate that T309 in CYP2D6 is involved in maintaining the balance of multiple oxygenating species and thus influences substrate and regioselectivity.