Subtle structural changes in the Asp251Gly/Gln307His P450 BM3 mutant responsible for new activity toward diclofenac, tolbutamide and ibuprofen

Sequence diversity in the monooxygenases involved in oxime production in plant defense and signaling: a conservative revision in the nomenclature of the highly complex CYP79 family

Cytochrome P450 monooxygenases of the CYP79 family catalyze conversion of specific amino acids into oximes feeding into a variety of metabolic plant pathways. Here we present an extensive phylogenetic tree of the CYP79 family built on carefully curated sequences collected across the entire plant kingdom. Based on a monophyletic origin of the P450s, a set of evolutionarily distinct branches was identified. Founded on the functionally characterized CYP79 sequences, sequence features of the individual substrate recognition sites (SRSs) were analyzed. Co-evolving amino acid residues were identified using co-evolutionary sequence analysis. SRS4 possesses a specific sequence pattern when tyrosine is a substrate. Except for the CYP79Cs and CYP79Fs, substrate preferences toward specific amino acids could not be assigned to specific subfamilies. The highly diversified CYP79 tree, reflecting recurrent independent evolution of CYP79s, may relate to the different roles of oximes in different plant species. The sequence differences across individual CYP79 subfamilies may facilitate the in vivo orchestration of channeled metabolic pathways based on altered surface charge domains of the CYP79 protein. Alternatively, they may serve to optimize dynamic interactions with oxime metabolizing enzymes to enable optimal ecological interactions. The outlined detailed curation of the CYP79 sequences used for building the phylogenetic tree made it appropriate to make a conservative phylogenetic tree-based revision of the naming of the sequences within this highly complex cytochrome P450 family. The same approach may be used in other complex P450 subfamilies. The detailed phylogeny of the CYP79 family will enable further exploration of the evolution of function in these enzymes.

Choose Your Own Adventure: A Comprehensive Database of Reactions Catalyzed by Cytochrome P450 BM3 Variants

Cytochrome P450 BM3 monooxygenase is the topic of extensive research as many researchers have evolved this enzyme to generate a variety of products. However, the abundance of information on increasingly diversified variants of P450 BM3 that catalyze a broad array of chemistry is not in a format that enables easy extraction and interpretation. We present a database that categorizes variants by their catalyzed reactions and includes details about substrates to provide reaction context. This database of >1500 P450 BM3 variants is downloadable and machine-readable and includes instructions to maximize ease of gathering information. The database allows rapid identification of commonly reported substitutions, aiding researchers who are unfamiliar with the enzyme in identifying starting points for enzyme engineering. For those actively engaged in engineering P450 BM3, the database, along with this review, provides a powerful and user-friendly platform to understand, predict, and identify the attributes of P450 BM3 variants, encouraging the further engineering of this enzyme.

Opportunities for Accelerating Drug Discovery and Development by Using Engineered Drug-Metabolizing Enzymes

The study of drug metabolism is fundamental to drug discovery and development (DDD) since by mediating the clearance of most drugs, metabolic enzymes influence their bioavailability and duration of action. Biotransformation can also produce pharmacologically active or toxic products, which complicates the evaluation of the therapeutic benefit versus liability of potential drugs but also provides opportunities to explore the chemical space around a lead. The structures and relative abundance of metabolites are determined by the substrate and reaction specificity of biotransformation enzymes and their catalytic efficiency. Preclinical drug biotransformation studies are done to quantify in vitro intrinsic clearance to estimate likely in vivo pharmacokinetic parameters, to predict an appropriate dose, and to anticipate interindividual variability in response, including from drug-drug interactions. Such studies need to be done rapidly and cheaply, but native enzymes, especially in microsomes or hepatocytes, do not always produce the full complement of metabolites seen in extrahepatic tissues or preclinical test species. Furthermore, yields of metabolites are usually limiting. Engineered recombinant enzymes can make DDD more comprehensive and systematic. Additionally, as renewable, sustainable, and scalable resources, they can also be used for elegant chemoenzymatic, synthetic approaches to optimize or synthesize candidates as well as metabolites. Here, we will explore how these new tools can be used to enhance the speed and efficiency of DDD pipelines and provide a perspective on what will be possible in the future. The focus will be on cytochrome P450 enzymes to illustrate paradigms that can be extended in due course to other drug-metabolizing enzymes. SIGNIFICANCE STATEMENT: Protein engineering can generate enhanced versions of drug-metabolizing enzymes that are more stable, better suited to industrial conditions, and have altered catalytic activities, including catalyzing non-natural reactions on structurally complex lead candidates. When applied to drugs in development, libraries of engineered cytochrome P450 enzymes can accelerate the identification of active or toxic metabolites, help elucidate structure activity relationships, and, when combined with other synthetic approaches, provide access to novel structures by regio- and stereoselective functionalization of lead compounds.

Enhanced and specific epoxidation activity of P450 BM3 mutants for the production of high value terpene derivatives

Terpenes are natural molecules of valuable interest for different industrial applications. Cytochromes P450 enzymes can functionalize terpenoids to form high value oxidized derivatives in a green and sustainable manner, representing a valid alternative to chemical catalysis. In this work, an enhanced and specific epoxidation activity of cytochrome P450 BM3 mutants was found for the terpenes geraniol and linalool. This is the first report showing the epoxidation of linalool by P450 BM3 and its mutant A2 (Asp251Gly/Gln307His) with the formation of valuable oxide derivatives, highlighting the relevance of this enzymes for industrial applications.

A Promiscuous Bacterial P450: The Unparalleled Diversity of BM3 in Pharmaceutical Metabolism

CYP102A1 (BM3) is a catalytically self-sufficient flavocytochrome fusion protein isolated from Bacillus megaterium, which displays similar metabolic capabilities to many drug-metabolizing human P450 isoforms. BM3's high catalytic efficiency, ease of production and malleable active site makes the enzyme a desirable tool in the production of small molecule metabolites, especially for compounds that exhibit drug-like chemical properties. The engineering of select key residues within the BM3 active site vastly expands the catalytic repertoire, generating variants which can perform a range of modifications. This provides an attractive alternative route to the production of valuable compounds that are often laborious to synthesize via traditional organic means. Extensive studies have been conducted with the aim of engineering BM3 to expand metabolite production towards a comprehensive range of drug-like compounds, with many key examples found both in the literature and in the wider industrial bioproduction setting of desirable oxy-metabolite production by both wild-type BM3 and related variants. This review covers the past and current research on the engineering of BM3 to produce drug metabolites and highlights its crucial role in the future of biosynthetic pharmaceutical production.

Polymorphism on human aromatase affects protein dynamics and substrate binding: spectroscopic evidence

Human aromatase is a member of the cytochrome P450 superfamily, involved in steroid hormones biosynthesis. In particular, it converts androgen into estrogens being therefore responsible for the correct sex steroids balance. Due to its capacity in producing estrogens it has also been considered as a promising target for breast cancer therapy. Two single-nucleotide polymorphisms (R264C and R264H) have been shown to alter aromatase activity and they have been associated to an increased or decreased risk for estrogen-dependent pathologies. Here, the effect of these mutations on the protein dynamics is investigated by UV/FTIR and time resolved fluorescence spectroscopy. H/D exchange rates were measured by FTIR for the three proteins in the ligand-free, substrate- and inhibitor-bound forms and the data indicate that the wild-type enzyme undergoes a conformational change leading to a more compact tertiary structure upon substrate or inhibitor binding. Indeed, the H/D exchange rates are decreased when a ligand is present. In the variants, the exchange rates in the ligand-free and -bound forms are similar, indicating that a structural change is lacking, despite the single amino acid substitution is located in the peripheral shell of the protein molecule. Moreover, the fluorescence lifetimes data show that the quenching effect on tryptophan-224 observed upon ligand binding in the wild-type, is absent in both variants. Since this residue is located in the catalytic pocket, these findings suggest that substrate entrance and/or retention in the active site is partially compromised in both mutants. A contact network analysis demonstrates that the protein structure is organized in two main clusters, whose connectivity is altered by ligand binding, especially in correspondence of helix-G, where the amino acid substitutions occur. Our findings demonstrate that SNPs resulting in mutations on aromatase surface modify the protein flexibility that is required for substrate binding and catalysis. The cluster analysis provides a rationale for such effect, suggesting helix G as a possible target for aromatase inhibition.

Optimisation of Cytochrome P450 BM3 Assisted by Consensus-Guided Evolution

Cytochrome P450 enzymes have attracted much interest over the years given their ability to insert oxygen into saturated carbon-hydrogen bonds, a difficult feat to accomplish by traditional chemistry. Much of the activity in this field has centered on the bacterial enzyme CYP102A1, or BM3, from Bacillus megaterium, as it has shown itself capable of hydroxylating/acting upon a wide range of substrates, thereby producing industrially relevant pharmaceuticals, fine chemicals, and hormones. In addition, unlike most cytochromes, BM3 is both soluble and fused to its natural redox partner, thus facilitating its use. The industrial use of BM3 is however stifled by its instability and its requirement for the expensive NADPH cofactor. In this work, we added several mutations to the BM3 mutant R966D/W1046S that enhanced the turnover number achievable with the inexpensive cofactors NADH and NBAH. These new mutations, A769S, S847G, S850R, E852P, and V978L, are localized on the reductase domain of BM3 thus leaving the oxidase domain intact. For NBAH-driven reactions by new mutant NTD5, this led to a 5.24-fold increase in total product output when compared to the BM3 mutant R966D/W1046S. For reactions driven by NADH by new mutant NTD6, this enhanced total product output by as much as 2.3-fold when compared to the BM3 mutant R966D/W1046S. We also demonstrated that reactions driven by NADH with the NTD6 mutant not only surpassed total product output achievable by wild-type BM3 with NADPH but also retained the ability to use this latter cofactor with greater total product output as well.

Production of metabolites of the anti-cancer drug noscapine using a P450BM3 mutant library

•

Mutants of P450_BM3 can metabolise noscapine.

•

Noscapine is N-demethylated with high selectivity.

•

The metabolites produced are of interest for drug development.

•

The profile of metabolites generated resembles that of mammalian CYP3A4.

Cytochrome P450 enzymes are a promising tool for the late-stage diversification of lead drug candidates and can provide an alternative route to structural modifications that are difficult to achieve with synthetic chemistry. In this study, a library of P450_BM3 mutants was produced using site-directed mutagenesis and the enzymes screened for metabolism of the opium poppy alkaloid noscapine, a drug with anticancer activity. Of the 18 enzyme mutants screened, 12 showed an ability to metabolise noscapine that was not present in the wild-type enzyme. Five noscapine metabolites were detected by LC-MS/MS, with the major metabolite for all mutants being N-demethylated noscapine. The highest observed regioselectivity for N-demethylation was 88%. Two hydroxylated metabolites, a catechol and two C-C cleavage products were also detected. P450-mediated production of hydroxylated and N-demethylated noscapine structures may be useful for the development of noscapine analogues with improved biological activity. The variation in substrate turnover, coupling efficiency and product distribution between the active mutants was considered alongside in silico docking experiments to gain insight into structural and functional effects of the introduced mutations. Selected mutants were identified as targets for further mutagenesis to improve activity and when coupled with an optimised process may provide a route for the preparative-scale production of noscapine metabolites.

Metabolism of liver CYP450 and ultrastructural changes after long-term administration of aspirin and ibuprofen

Worldwide, aspirin and ibuprofen are the most commonly used non-steroidal anti-inflammatory drugs (NSAIDs). Some adverse reactions, including gastrointestinal reactions, have been concerned extensively. Nevertheless, the mechanism of liver injury remains unclear. In the present study, we focused on the metabolism of liver cytochrome P450 (CYP450) and ultrastructural morphology of liver cells. A total of thirty rats were divided into three groups of 10. Rats in the aspirin and ibuprofen groups were given enteric-coated aspirin (15 mg/kg) and ibuprofen (15 mg/kg), respectively by gavage for four weeks. The body weights were recorded every two days. Liver function and metabolic capacity of CYP450 were studied on days 14 and 28. We then conducted ultrastructural examinations. Body weights in the Ibuprofen group were lower than those of the Control group, and ALT and AST levels were significantly higher (P < 0.05). There were no significant differences in terms of body weight, ALT or AST between the Aspirin and Control groups. The metabolic capacity of CYP450 was evaluated using five probe drugs, phenacetin, tolbutamide, metoprolol, midazolam, and bupropion. We found that ibuprofen and aspirin induced metabolism of the probe drugs. Moreover, according to the pharmacokinetic data, the Control, Aspirin and Ibuprofen groups could be discriminated accurately. Ultrastructural examination showed that the number of mitochondria was increased in both the Ibuprofen and Aspirin groups. Long-term administration of enteric-coated aspirin and ibuprofen induced the metabolic activity of the CYP450 enzyme. Aspirin had better tolerability than did ibuprofen, as reflected by pharmacokinetic data of probe drug metabolism.

Molecular Determinants of Substrate Affinity and Enzyme Activity of a Cytochrome P450BM3 Variant

Cytochrome P450_BM3 catalyzes the hydroxylation and/or epoxidation of fatty acids, fatty amides, and alcohols. Protein engineering has produced P450_BM3 variants capable of accepting drug molecules normally metabolized by human P450 enzymes. The enhanced substrate promiscuity has been attributed to the greater flexibility of the lid of the substrate channel. However, it is not well understood how structurally different and highly polar drug molecules can stably bind in the active site nor how the activity and coupling efficiency of the enzyme may be affected by the lack of enzyme-substrate complementarity. To address these important aspects of non-native small molecule binding, this study investigated the binding of drug molecules with different size, charge, polar surface area, and human P450 affinity on the promiscuous R47L/F87V/L188Q/E267V/F81I pentuple mutant of P450_BM3. Binding free energy data and energy decomposition analysis showed that pentuple mutant P450_BM3 stably binds (i.e., negative ΔG_b°) a broad range of substrate and inhibitor types because dispersion interactions with active site residues overcome unfavorable repulsive and electrostatic effects. Molecular dynamics simulations revealed that 1) acidic substrates tend to disrupt the heme propionate A-K69 salt bridge, which may reduce heme oxidizing ability, and 2) the lack of complementarity leads to high substrate mobility and water density in the active site, which may lead to uncoupling. These factors must be considered in future developments of P450_BM3 as a biocatalyst in the large-scale production of drug metabolites.

Inactivation mechanism of N61S mutant of human FMO3 towards trimethylamine

Human flavin-containing monooxygenase 3 (hFMO3) catalyses the oxygenation of a wide variety of compounds including drugs as well as dietary compounds. It is the major hepatic enzyme involved in the production of the N-oxide of trimethylamine (TMAO) and clinical studies have uncovered a striking correlation between plasma TMAO concentration and cardiovascular disease. Certain mutations within the hFMO3 gene cause defective trimethylamine (TMA) N-oxygenation leading to trimethylaminuria (TMAU) also known as fish-odour syndrome. In this paper, the inactivation mechanism of a TMAU-causing polymorphic variant, N61S, is investigated. Transient kinetic experiments show that this variant has a > 170-fold lower NADPH binding affinity than the wild type. Thermodynamic and spectroscopic experiments reveal that the poor NADP⁺ binding affinity accelerates the C4a-hydroperoxyFAD intermediate decay, responsible for an unfavourable oxygen transfer to the substrate. Steady-state kinetic experiments show significantly decreased N61S catalytic activity towards other substrates; methimazole, benzydamine and tamoxifen. The in vitro data are corroborated by in silico data where compared to the wild type enzyme, a hydrogen bond required for the stabilisation of the flavin intermediate is lacking. Taken together, the data presented reveal the molecular basis for the loss of function observed in N61S mutant.

Crystallographic insights into a cobalt (III) sepulchrate based alternative cofactor system of P450 BM3 monooxygenase

P450 BM3 is a multi-domain heme-containing soluble bacterial monooxygenase. P450 BM3 and variants are known to oxidize structurally diverse substrates. Crystal structures of individual domains of P450 BM3 are available. However, the spatial organization of the full-length protein is unknown. In this study, crystal structures of the P450 BM3 M7 heme domain variant with and without cobalt (III) sepulchrate are reported. Cobalt (III) sepulchrate acts as an electron shuttle in an alternative cofactor system employing zinc dust as the electron source. The crystal structure shows a binding site for the mediator cobalt (III) sepulchrate at the entrance of the substrate access channel. The mediator occupies an unusual position which is far from the active site and distinct from the binding of the natural redox partner (FAD/NADPH binding domain).

Human Cytochrome P450 3A4 as a Biocatalyst: Effects of the Engineered Linker in Modulation of Coupling Efficiency in 3A4-BMR Chimeras

Human liver cytochrome P450 3A4 is the main enzyme involved in drug metabolism. This makes it an attractive target for biocatalytic applications, such as the synthesis of pharmaceuticals and drug metabolites. However, its poor solubility, stability and low coupling have limited its application in the biotechnological context. We previously demonstrated that the solubility of P450 3A4 can be increased by creating fusion proteins between the reductase from Bacillus megaterium BM3 (BMR) and the N-terminally modified P450 3A4 (3A4-BMR). In this work, we aim at increasing stability and coupling efficiency by varying the length of the loop connecting the two domains to allow higher inter-domain flexibility, optimizing the interaction between the domains. Starting from the construct 3A4-BMR containing the short linker Pro-Ser-Arg, two constructs were generated by introducing a 3 and 5 glycine hinge (3A4-3GLY-BMR and 3A4-5GLY-BMR). The three fusion proteins show the typical absorbance at 450 nm of the reduced heme-CO adduct as well as the correct incorporation of the FAD and FMN cofactors. Each of the three chimeric proteins were more stable than P450 3A4 alone. Moreover, the 3A4-BMR-3-GLY enzyme showed the highest NADPH oxidation rate in line with the most positive reduction potential. On the other hand, the 3A4-BMR-5-GLY fusion protein showed a Vmax increased by 2-fold as well as a higher coupling efficiency when compared to 3A4-BMR in the hydroxylation of the marker substrate testosterone. This protein also showed the highest rate value of cytochrome c reduction when this external electron acceptor is used to intercept electrons from BMR to P450. The data suggest that the flexibility and the interaction between domains in the chimeric proteins is a key parameter to improve turnover and coupling efficiency. These findings provide important guidelines in engineering catalytically self-sufficient human P450 for applications in biocatalysis.

Impact of R264C and R264H polymorphisms in human aromatase function

The cytochrome P450 aromatase is involved in the last step of sex hormones biosynthesis by converting androgens into estrogens. The human enzyme is highly polymorphic and literature data correlate aromatase single nucleotide polymorphisms to the onset of pathologies such as breast cancer and neurodegenerative diseases. The aims of this study were i) to study the influence of the mutations R264C and R264H on the structure-function of the enzyme also upon phosphorylation by selected kinases and ii) to compare the activity of the variants to that of aromatase wild type in two different cell lines. Far-UV circular dichroism spectroscopy, thermal denaturation experiments and CO-binding assay showed that the two polymorphic variants are correctly folded. Steady-state kinetics experiments showed that rArom R264C and R264H exhibit a 1.5 and 3.4 folds lower catalytic efficiency, respectively, when compared to the wild type protein. Since R264 is part of the consensus motif of PKA and PKG1, phosphorylation experiments were performed to study the effect on aromatase function. Phosphorylation by PKA caused a decrease in activity by 36.2%, 49.3% and 27.9% in the wild type, R264C and R264H proteins respectively. Phosphorylation by PKG1 was also found to decrease the activity by 30.3%, 30.5% and 15.4% in the wild type, R264C and R264H proteins respectively. Experiments performed on the three full-length proteins expressed in human MCF-7 breast cancer cells and rat ST14A neuronal cells showed that, depending on the cell line used, the activity of the proteins is different, implicating different cellular mechanisms modulating aromatase activity. This work demonstrate that R264 polymorphism causes an intrinsic alteration of aromatase activity together with a different consensus for phosphorylation by different kinases, indicating that estrogen production can be different when such mutations are present. These findings are significant in understanding the onset and treatment of pathologies in which aromatase has been shown to be involved.