An electron transfer competent structural ensemble of membrane-bound cytochrome P450 1A1 and cytochrome P450 oxidoreductase

On the determinants of electron transfer reorganization energy in a cytochrome P450: cytochrome b5 complex. A combined quantum mechanics and molecular dynamics simulation study

The electron transfer steps in the catalytic cycle of cytochrome P450 (CYP) enzymes, ubiquitous proteins with key roles in processes such as drug metabolism and steroidogenesis, are often rate-limiting. To predict ET rates from atomistic molecular dynamics simulations using Marcus theory, values of the reaction free energy ΔG0 and the reorganization free energy λ are required from either experiments or computations. For the reduction of cytochrome P450 17A1 (CYP17A1) by the secondary redox protein cytochrome b5 (CYb5), a critical step in the regulation of steroidogenesis, experimental measurements of λ are not available. We here describe the computation of λ for this system from a combination of molecular mechanics/molecular dynamics simulations and quantum mechanics computations. Our results show that a quantum mechanical treatment of the redox-active cofactors is necessary, even though the surrounding protein and solvent, which are modeled classically, contribute most to the reorganization energy. The values of λ computed for structural ensembles corresponding to two predicted binding modes of the proteins are 1.23 and 1.16 eV. We find that the λ values computed for the individual soluble globular domains of the two proteins sum to approximately the λ values computed for the membrane-bound CYP17A1-CYb5 complex, indicating that additivity can be invoked in a computationally efficient approach to estimating λ values for such protein-protein complexes.

Revealing New Insights into the Dynamics of Human Aromatase Interacting with Cytochrome P450 Reductase in a Realistic Membrane Environment

Molecular dynamics simulations were applied to human aromatase (HA) complexed with cytochrome P450 reductase (CPR) within a realistic endoplasmic reticulum membrane environment to evaluate its structural and dynamical properties. CPR was examined to have a specific point mutation (P281T), where proline was substituted by threonine, which is envisaged to demonstrate a far-reaching influence on its structure, dynamics, and electron transfer behavior. Since CPR plays a key role in the electron transfer to HA, catalyzing steroidogenesis, obtaining detailed information on the mutation effect of CPR on HA was crucial. This compelled us to study the interaction of HA with CPR in its wild-type and mutant forms, enabling the investigation of different properties of CPR and its effect on HA dynamics. Pursuing these objectives, different analytical parameters, notably, root-mean-square deviation, root-mean-square fluctuation, dynamic cross-correlation matrix, and principal component analysis were applied to gain insight into the conformational dynamics of the HA/CPR complex. These analyses demonstrated the CPR shifts in the complexes' dynamics, in the context of the effect of CPR mutation effects on HA behavior. Based on this information, the electron transfer within the HA/CPR complex was also envisaged to be influenced by the contrast dynamics between the two complexes as the mutation was evaluated to significantly alter the dynamics of CPR as well as HA. Furthermore, the electron transfer within the complexes was determined by applying Marcus theory of electron transfer, revealing a contrast between the HA/wild-type and mutant CPR complexes. The latter was found to alter the electron transfer efficiency, demonstrating a direct effect of changes in the protein dynamics observed within the HA/mutant CPR complex. This study therefore provides valuable insights into the conformational dynamics of HA in conjunction with CPR, affecting the electron transfer process and their potential implications for understanding estrogen-related physiological conditions influenced by these proteins.

Modification Strategy of Two-Step Ordered Selective Hydroxylation of DHEA Catalyzed by CYP-cl3

P450 hydroxylation is widely used in steroidal drug synthesis, which can not only activate the oxidation of inert C-H bonds but also increase the pharmacological activity of steroidal compounds after hydroxylation modification. The two-step ordered selective hydroxylation of steroids catalyzed by P450 dihydroxylases often forms the key precursors or intermediates of steroid drugs. However, the poor catalytic performance of P450 dihydroxylase and the imbalance of two-step hydroxylation lead to the accumulation of intermediates, which restricts the efficient synthesis of dihydroxyl products. In this work, the hydroxylation sequence (first at C7 and then at C15) of DHEA to 7α,15α-diOH-DHEA catalyzed by CYP-cl3 from Colletotrichum lini ST-1 was revealed with the aid of computational analysis. To obtain highly C7α- and C15α-selective enzymes, we designed a three-step modification strategy that includes the high-throughput screening of C7-hydroxylation mutants, semirational design of C15-hydroxylation selectivity, and the combination of dominant mutation sites. As a result, we successfully obtained a dominant quadruple mutant A83P/E264I/V281A/T315P, whose the proportion of 7α,15α-diOH-DHEA reached 99.9%. The mechanism of key amino acid residues in improving the catalytic performance of CYP-cl3 was revealed by molecular docking and molecular dynamics simulation analysis. Our study guides the performance improvement of other P450 hydroxylase.

Geometry-encoded molecular dynamics enables deep learning insights into P450 regiospecificity control

Cytochrome P450 1A2, as many isoenzymes, can generate multiple metabolites from a single substrate. A loose coupling between substrate binding and oxygen activation makes possible substrate reorientations at the active site prior to catalysis. In the present work, caffeine oxidation to alternative bioactive compounds was used to decipher this pluripotency. A model involving two interacting subsites capable of sequentially accommodating one or two caffeine molecules was considered. Molecular dynamics was used to characterize subsite interactions and feed a dedicated geometric encoding of trajectories that was coupled to dimensional reductions and differential machine learning. The two subsites differentially control caffeine orientations and can exchange substrate through a phenylalanine gated mechanism. This exchange can be locked by the presence of a second bound molecule. Complementary roles of subsites in progressively determining the caffeine orientation during its approach to active oxygen were examined. Interestingly, substrate face flipping becomes impaired upon entry into the rather flat active site. This makes the mechanisms that define the orientation of caffeine relative to active oxygen dependent on the substrate face oriented toward heme. Globally, this evidenced that P450 1A2 regioselectivity results from local determinants combined with subsite interactions and caffeine face preselection at a longer distance.

Protein subinteractomes of human microsomal cytochromes P450

Microsomal cytochromes P450 (micCYPs) are monooxygenases located in the endoplasmic reticulum and other endomembranes of human cells. micCYPs receive electrons from specific redox partners and perform enzymatic transformations of drugs and different endogenous substrates. The large biodiversity of micCYPs leads to the idea that protein-protein interactions (PPIs) involving micCYPs are not limited to classical redox partners. This review aims to perform a systems biology analysis of the complete set of PPIs for all 33 micCYPs studied, as well as to examine the subinteractome of each micCYP. We have retrieved 287 PPIs from interactomic databases, involving 246 unique protein interactors that share a similar profile of subcellular localization with micCYPs. The number of protein interactors per micCYP unevenly varies from one to 47. Interactors of micCYPs are involved in cellular metabolism, signal transduction, cell-cell junctions, cytoskeleton organization, and intracellular or transmembrane transport. Notably, up to one-third of all interactors belong to the latter group, half of which consists of membrane transporters of compounds, metabolites, and ions (e.g., CACNA2D1, ORAI1, SCN3B, SLC7A2, SLC19A3, and SLC11A2). The CYP2C8 subinteractome is enriched with proteins involved in autophagy; CYP2S1- ERBB2 and EPH-Ephrin signaling; CYP3A4- glucuronidation. Proteins UBC, PGRMC1, and FANCG are the most frequent common interactors across various micCYPs. Nine and 12 interactors of micCYPs are involved in phosphorylation and ubiquitination, respectively; 20 interactors are 'moonlighting' proteins that are represented in the CYP3A4 subinteractome. Furthermore, micCYPs such as CYP2C9, 3A5, 2E1, 2A6, 4F2, and 4A11 may be involved in potentially binary interactions with other micCYPs. The functional implication of these CYP-CYP pairs is likely associated with modulation of their activity. Analysis of transcriptomic data revealed that some micCYP/interactor pairs exhibit tissue-, time-, and disease-specific gene expression patterns. Drugs that are metabolized by micCYPs in some cases can influence the expression of corresponding interactors at the gene or protein levels. These findings suggest that micCYPs may play roles in functions beyond their monooxygenase activity, as indicated by the spectrum of PPIs analyzed.

A Novel A105Y Mutant of CYP17A1 Exhibits Almost Perfect Regioselectivity in the Production of 17α-Hydroxyprogesterone

17α-Hydroxyprogesterone (17α-OHP) is a steroid hormone with significant biological activity that can be obtained by catalyzing progesterone (PROG), the main product of sitosterol, through CYP17A1. However, increasing the catalytic specificity of HCYP17A1 for C17 hydroxylation of progesterone (PROG) poses a formidable challenge due to the close proximity of the C16 and C17 positions. In this study, a rational design was utilized to alter the spatial configuration of the substrate channel, leading to the complete abolition of its 16-hydroxylation activity. Subsequent molecular dynamics simulations revealed that the A105Y mutation heightened the rigidity of the G95-I112 region of CYP17A1, consequently regulating the direction of the entry of PROG into the catalytic pocket. Moreover, the establishment of hydrogen bonding between Y105 and R239, coupled with Pi-stacking of A105Y with F114, effectively immobilizes the substrate PROG in a fixed position, explaining the practically perfect regioselectivity observed in A105Y. Finally, a multifaceted enzymatic cascade system, incorporating A105Y, cytochrome P450 reductase (CPR), and glucose-6-phosphate dehydrogenase (ZWF) for NADPH cofactor regeneration, was constructed in Pichia pastoris GS115. The resulting biocatalyst produced 106 ± 3.2 mg L-1 17α-OHP, a 4.6-fold increase compared with A105Y alone. Thus, this study provides valuable insights for improving the regioselectivity and activity of P450 enzymes.

Multiresolution molecular dynamics simulations reveal the interplay between conformational variability and functional interactions in membrane-bound cytochrome P450 2B4

Cytochrome P450 2B4 (CYP 2B4) is one of the best-characterized CYPs and serves as a key model system for understanding the mechanisms of microsomal class II CYPs, which metabolize most known drugs. The highly flexible nature of CYP 2B4 is apparent from crystal structures that show the active site with either a wide open or a closed heme binding cavity. Here, we investigated the conformational ensemble of the full-length CYP 2B4 in a phospholipid bilayer, using multiresolution molecular dynamics (MD) simulations. Coarse-grained MD simulations revealed two predominant orientations of CYP 2B4's globular domain with respect to the bilayer. Their refinement by atomistic resolution MD showed adaptation of the enzyme's interaction with the lipid bilayer, leading to open configurations that facilitate ligand access to the heme binding cavity. CAVER analysis of enzyme tunnels, AquaDuct analysis of water routes, and Random Acceleration Molecular Dynamics simulations of ligand dissociation support the conformation-dependent passage of molecules between the active site and the protein surroundings. Furthermore, simulation of the re-entry of the inhibitor bifonazole into the open conformation of CYP 2B4 resulted in binding at a transient hydrophobic pocket within the active site cavity that may play a role in substrate binding or allosteric regulation. Together, these results show how the open conformation of CYP 2B4 facilitates the binding of substrates from and release of products to the membrane, whereas the closed conformation prolongs the residence time of substrates or inhibitors and selectively allows the passage of smaller reactants via the solvent and water channels.

The Dissociation Process of NADP+ from NADPH-Cytochrome P450 Reductase Studied by Molecular Dynamics Simulation

NADPH-cytochrome P450 reductase (CPR) plays a vital role as a redox partner for mammalian cytochrome P450 enzymes (P450s), facilitating the transfer of two electrons from NADPH to the P450 heme center in a sequential manner. Previous experimental studies revealed substantial domain movements of CPR, transitioning between closed and open states during the electron transfer (ET) cycle. These transitions are essential and are influenced by the binding of NADPH or the release of NADP+. However, the intricate molecular mechanisms governing the CPR-mediated ET cycle have largely remained elusive. This study employed molecular dynamics (MD) simulation techniques to explore the dissociation of NADP+ from CPR, a crucial step preceding the initial ET from CPR to a P450. Alongside the binding structure of NADP+ observed in a crystal structure (pose I), our MD simulations identified an alternative binding structure (pose II). Although pose II exhibits slightly lower stability than pose I, it can be formed through an approximate 210° counterclockwise rotation of the adenine group, with a free energy barrier of only 2.76 kcal/mol. The simulation results further suggest that NADP+ dissociation involves a tentative formation of pose II from pose I before complete dissociation, and that the binding of NADP+ to CPR is primarily governed by nonbonded interactions within the adenosine binding pocket.

Effect of Lipid Bilayer Anchoring on the Conformational Properties of the Cytochrome P450 2D6 Binding Site

Human cytochrome P450 (CYP) proteins metabolize 75% of small-molecule pharmaceuticals, which makes structure-based modeling of CYP metabolism and inhibition, bolstered by the current availability of X-ray crystal structures of CYP globular catalytic domains, an attractive prospect. Accounting for this broad metabolic capacity is a combination of the existence of multiple different CYP proteins and the capacity of a single CYP protein to metabolize multiple different small molecules. It is thought that structural plasticity and flexibility contribute to this latter property; therefore, incorporating diverse conformational states of a particular CYP is likely an important consideration in structure-based CYP metabolism and inhibition modeling. While all-atom explicit-solvent molecular dynamics simulations can be used to generate conformational ensembles under biologically relevant conditions, existing CYP crystal structures are of the globular domain only, whereas human CYPs contain N-terminal transmembrane and linker peptides that anchor the globular catalytic domain to the endoplasmic reticulum. To determine whether this can cause significant differences in the sampled binding site conformations, microsecond scale all-atom explicit-solvent molecular dynamics simulations of the CYP2D6 globular domain in an aqueous environment were compared with those of the full-length protein anchored in a POPC lipid bilayer. While bilayer-anchoring damped some structural fluctuations in the globular domain relative to the aqueous simulations, none of the affected residues included binding site pocket residues. Furthermore, clustering of molecular dynamics snapshots based on either pairwise binding site pocket RMSD or volume differences demonstrated a lack of separation of snapshots from the two simulation conditions into different clusters. These results suggest the substantially simpler and computationally cheaper aqueous simulation approach can be used to generate a relevant conformational ensemble of the CYP2D6 binding site for structure-based metabolism and inhibition modeling.

Sterol 14-alpha demethylase (CYP51) activity in Leishmania donovani is likely dependent upon cytochrome P450 reductase 1

Liposomal amphotericin B is an important frontline drug for the treatment of visceral leishmaniasis, a neglected disease of poverty. The mechanism of action of amphotericin B (AmB) is thought to involve interaction with ergosterol and other ergostane sterols, resulting in disruption of the integrity and key functions of the plasma membrane. Emergence of clinically refractory isolates of Leishmania donovani and L. infantum is an ongoing issue and knowledge of potential resistance mechanisms can help to alleviate this problem. Here we report the characterisation of four independently selected L. donovani clones that are resistant to AmB. Whole genome sequencing revealed that in three of the moderately resistant clones, resistance was due solely to the deletion of a gene encoding C24-sterol methyltransferase (SMT1). The fourth, hyper-resistant resistant clone (>60-fold) was found to have a 24 bp deletion in both alleles of a gene encoding a putative cytochrome P450 reductase (P450R1). Metabolic profiling indicated these parasites were virtually devoid of ergosterol (0.2% versus 18% of total sterols in wild-type) and had a marked accumulation of 14-methylfecosterol (75% versus 0.1% of total sterols in wild-type) and other 14-alpha methylcholestanes. These are substrates for sterol 14-alpha demethylase (CYP51) suggesting that this enzyme may be a bona fide P450R specifically involved in electron transfer from NADPH to CYP51 during catalysis. Deletion of P450R1 in wild-type cells phenocopied the metabolic changes observed in our AmB hyper-resistant clone as well as in CYP51 nulls. Likewise, addition of a wild type P450R1 gene restored sterol profiles to wild type. Our studies indicate that P450R1 is essential for L. donovani amastigote viability, thus loss of this gene is unlikely to be a driver of clinical resistance. Nevertheless, investigating the mechanisms underpinning AmB resistance in these cells provided insights that refine our understanding of the L. donovani sterol biosynthetic pathway.

Charged Amino Acids in the Transmembrane Helix Strongly Affect the Enzyme Activity of Aromatase

Estrogens play critical roles in embryonic development, gonadal sex differentiation, behavior, and reproduction in vertebrates and in several human cancers. Estrogens are synthesized from testosterone and androstenedione by the endoplasmic reticulum membrane-bound P450 aromatase/cytochrome P450 oxidoreductase complex (CYP19/CPR). Here, we report the characterization of novel mammalian CYP19 isoforms encoded by CYP19 gene copies. These CYP19 isoforms are all defined by a combination of mutations in the N-terminal transmembrane helix (E42K, D43N) and in helix C of the catalytic domain (P146T, F147Y). The mutant CYP19 isoforms show increased androgen conversion due to the KN transmembrane helix. In addition, the TY substitutions in helix C result in a substrate preference for androstenedione. Our structural models suggest that CYP19 mutants may interact differently with the membrane (affecting substrate uptake) and with CPR (affecting electron transfer), providing structural clues for the catalytic differences.

Mutation to a cytochrome P450 -like gene alters the leaf color by affecting the heme and chlorophyll biosynthesis pathways in Brassica napus

The regulated biosynthesis of chlorophyll is important because of its effects on plant photosynthesis and dry biomass production. In this study, a map-based cloning approach was used to isolate the cytochrome P450 -like gene BnaC08g34840D (BnCDE1) from a chlorophyll-deficient mutant (cde1) of Brassica napus obtained by ethyl methanesulfonate (EMS) mutagenization. Sequence analyses revealed that BnaC08g34840D in the cde1 mutant (BnCDE1I320T ) encodes a substitution at amino acid 320 (Ile320Thr) in the conserved region. The over-expression of BnCDE1I320T in ZS11 (i.e., gene-mapping parent with green leaves) recapitulated a yellow-green leaf phenotype. The CRISPR/Cas9 genome-editing system was used to design two single-guide RNAs (sgRNAs) targeting BnCDE1I320T in the cde1 mutant. The knockout of BnCDE1I320T in the cde1 mutant via a gene-editing method restored normal leaf coloration (i.e., green leaves). These results indicate that the substitution in BnaC08g34840D alters the leaf color. Physiological analyses showed that the over-expression of BnCDE1I320T leads to decreases in the number of chloroplasts per mesophyll cell and in the contents of the intermediates of the chlorophyll biosynthesis pathway in leaves, while it increases heme biosynthesis, thereby lowering the photosynthetic efficiency of the cde1 mutant. The Ile320Thr mutation in the highly conserved region of BnaC08g34840D inhibited chlorophyll biosynthesis and disrupted the balance between heme and chlorophyll biosynthesis. Our findings may further reveal how the proper balance between the chlorophyll and heme biosynthesis pathways is maintained.

Influence of cholesterol on kinetic parameters for human aromatase (P450 19A1) in phospholipid nanodiscs

Cholesterol, a significant constituent of the endoplasmic reticulum membrane, exerts a substantial effect on the membrane's biophysical and mechanical properties. Cholesterol, however, is often neglected in model systems used to study membrane-bound proteins. For example, the influence of cholesterol on the enzymatic functions of type 2 cytochromes P450, which require a phospholipid bilayer and the redox partner P450-oxidoreductase (POR) for activity, are rarely investigated. Human aromatase (P450 19A1) catalyzes three sequential oxygenations of 19‑carbon steroids to estrogens and is widely expressed across various tissues, which are characterized by varying cholesterol compositions. Our study examined the impact of cholesterol on the functionality of the P450 19A1 complex with POR. Nanodiscs containing P450 19A1 with 20% cholesterol/80% phospholipid had similar rates and affinity of androstenedione binding as phospholipid-only P450 19A1 nanodiscs, and rates of product formation were indistinguishable among these conditions. In contrast, the rate of the first electron transfer from POR to P450 19A1 was 3-fold faster in cholesterol-containing nanodiscs than in phospholipid-only nanodiscs. These results suggest that cholesterol influences some aspects of POR interaction with P450 19A1 and might serve as an additional regulatory mechanism in this catalytic system.

Single Mutations in Cytochrome P450 Oxidoreductase Can Alter the Specificity of Human Cytochrome P450 1A2-Mediated Caffeine Metabolism

A unique cytochrome P450 (CYP) oxidoreductase (CPR) sustains activities of human microsomal CYPs. Its function requires toggling between a closed conformation enabling electron transfers from NADPH to FAD and then FMN cofactors and open conformations forming complexes and transferring electrons to CYPs. We previously demonstrated that distinct features of the hinge region linking the FAD and FMN domain (FD) modulate conformer poses and their interactions with CYPs. Specific FD residues contribute in a CYP isoform-dependent manner to the recognition and electron transfer mechanisms that are additionally modulated by the structure of CYP-bound substrate. To obtain insights into the underlying mechanisms, we analyzed how hinge region and FD mutations influence CYP1A2-mediated caffeine metabolism. Activities, metabolite profiles, regiospecificity and coupling efficiencies were evaluated in regard to the structural features and molecular dynamics of complexes bearing alternate substrate poses at the CYP active site. Studies reveal that FD variants not only modulate CYP activities but surprisingly the regiospecificity of reactions. Computational approaches evidenced that the considered mutations are generally in close contact with residues at the FD-CYP interface, exhibiting induced fits during complexation and modified dynamics depending on caffeine presence and orientation. It was concluded that dynamic coupling between FD mutations, the complex interface and CYP active site exist consistently with the observed regiospecific alterations.

β-Cryptoxanthin Production in Escherichia coli by Optimization of the Cytochrome P450 CYP97H1 Activity

Cytochromes P450, forming a superfamily of monooxygenases containing heme as a cofactor, show great versatility in substrate specificity. Metabolic engineering can take advantage of this feature to unlock novel metabolic pathways. However, the cytochromes P450 often show difficulty being expressed in a heterologous chassis. As a case study in the prokaryotic host Escherichia coli, the heterologous synthesis of β-cryptoxanthin was addressed. This carotenoid intermediate is difficult to produce, as its synthesis requires a monoterminal hydroxylation of β-carotene whereas most of the classic carotene hydroxylases are dihydroxylases. This study was focused on the optimization of the in vivo activity of CYP97H1, an original P450 β-carotene monohydroxylase. Engineering the N-terminal part of CYP97H1, identifying the matching redox partners, defining the optimal cellular background and adjusting the culture and induction conditions improved the production by 400 times compared to that of the initial strain, representing 2.7 mg/L β-cryptoxanthin and 20% of the total carotenoids produced.

Metabolic Engineering of Yarrowia lipolytica for Terpenoid Production: Tools and Strategies

Terpenoids are a diverse group of compounds with isoprene units as basic building blocks. They are widely used in the food, feed, pharmaceutical, and cosmetic industries due to their diverse biological functions such as antioxidant, anticancer, and immune enhancement. With an increase in understanding the biosynthetic pathways of terpenoids and advances in synthetic biology techniques, microbial cell factories have been built for the heterologous production of terpenoids, with the oleaginous yeast Yarrowia lipolytica emerging as an outstanding chassis. In this paper, recent progress in the development of Y. lipolytica cell factories for terpenoid production with a focus on the advances in novel synbio tools and metabolic engineering strategies toward enhanced terpenoid biosynthesis is reviewed.

Cytochrome P450 Surface Domains Prevent the β-Carotene Monohydroxylase CYP97H1 of Euglena gracilis from Acting as a Dihydroxylase

Molecular biodiversity results from branched metabolic pathways driven by enzymatic regioselectivities. An additional complexity occurs in metabolites with an internal structural symmetry, offering identical extremities to the enzymes. For example, in the terpene family, β-carotene presents two identical terminal closed-ring structures. Theses cycles can be hydroxylated by cytochrome P450s from the CYP97 family. Two sequential hydroxylations lead first to the formation of monohydroxylated β-cryptoxanthin and subsequently to that of dihydroxylated zeaxanthin. Among the CYP97 dihydroxylases, CYP97H1 from Euglena gracilis has been described as the only monohydroxylase. This study aims to determine which enzymatic domains are involved in this regioselectivity, conferring unique monohydroxylase activity on a substrate offering two identical sites for hydroxylation. We explored the effect of truncations, substitutions and domain swapping with other CYP97 members and found that CYP97H1 harbours a unique N-terminal globular domain. This CYP97H1 N-terminal domain harbours a hydrophobic patch at the entrance of the substrate channel, which is involved in the monohydroxylase activity of CYP97H1. This domain, at the surface of the enzyme, highlights the role of distal and non-catalytic domains in regulating enzyme specificity.

Phylogenetic Analysis and Protein Modelling of Isoflavonoid Synthase Highlights Key Catalytic Sites towards Realising New Bioengineering Endeavours

Isoflavonoid synthase (IFS) is a critical enzyme for the biosynthesis of over 2400 isoflavonoids. Isoflavonoids are an important class of plant secondary metabolites that have a range of pharmaceutical and nutraceutical properties. With growing interest in isoflavonoids from both research and industrial perspectives, efforts are being forwarded to enhance isoflavonoid production in-planta and ex-planta; therefore, in-silico analysis and characterisation of available IFS protein sequences are needed. The present study is the first-ever attempt toward phylogenetic analysis and protein modelling of available IFS protein sequences. Phylogenetic analysis has shown that IFS amino acid sequences have 86.4% pairwise identity and 26.5% identical sites, and the sequences were grouped into six different clades. The presence of a β-hairpin and extra loop at catalytic sites of Trifolium pratense, Beta vulgaris and Medicago truncatula, respectively, compared with Glycyrrhiza echinata are critical structural differences that may affect catalytic function. Protein docking highlighted the preference of selected IFS for liquiritigenin compared with naringenin and has listed T. pratense as the most efficient candidate for heterologous biosynthesis of isoflavonoids. The in-silico characterisation of IFS represented in this study is vital in realising the new bioengineering endeavours and will help in the characterisation and selection of IFS candidate enzymes for heterologous biosynthesis of isoflavonoids.

Non-Ionic Inulin-Based Polymer Nanodiscs Enable Functional Reconstitution of a Redox Complex Composed of Oppositely Charged CYP450 and CPR in a Lipid Bilayer Membrane

Although polymer-based lipid nanodiscs are increasingly used in the structural studies of membrane proteins, the charge of the belt-forming polymer is a major limitation for functional reconstitution of membrane proteins possessing an opposite net charge to that of the polymer. This limitation also rules out the reconstitution of a protein-protein complex composed of oppositely charged membrane proteins. In this study, we report the first successful functional reconstitution of a membrane-bound redox complex constituting a cationic cytochrome P450 (CYP450) and an anionic cytochrome P450 reductase (CPR) in non-ionic inulin-based lipid nanodiscs. The gel-to-liquid-crystalline phase-transition temperature (Tm) of DMPC:DMPG (7:3 w/w) lipids in polymer nanodiscs was determined by differential scanning calorimetry (DSC) and 31P NMR experiments. The CYP450-CPR redox complex reconstitution in polymer nanodiscs was characterized by size-exclusion chromatography (SEC), and the electron transfer kinetics was carried out using the stopped-flow technique under anaerobic conditions. The Tm of DMPC:DMPG (7:3 w/w) in polymer nanodiscs measured from 31P NMR agrees with that obtained from DSC and was found to be higher than that for liposomes due to the decreased cooperativity of lipids present in the nanodiscs. The stopped-flow measurements revealed the CYP450-CPR redox complex reconstituted in nanodiscs to be functional, and the electron transfer kinetics was found to be temperature-dependent. Based on the successful demonstration of the use of non-ionic inulin-based polymer nanodiscs reported in this study, we expect them to be useful in studying the function and structures of a variety of membrane proteins/complexes irrespective of the charge of the molecular components. Since the polymer nanodiscs were shown to align in an externally applied magnetic field, they can also be used to measure residual dipolar couplings (RDCs) and residual quadrupolar couplings (RQCs) for various molecules ranging from small molecules to soluble proteins and nucleic acids.

Use of engineered cytochromes P450 for accelerating drug discovery and development

Numerous steps in drug development, including the generation of authentic metabolites and late-stage functionalization of candidates, necessitate the modification of often complex molecules, such as natural products. While it can be challenging to make the required regio- and stereoselective alterations to a molecule using purely chemical catalysis, enzymes can introduce changes to complex molecules with a high degree of stereo- and regioselectivity. Cytochrome P450 enzymes are biocatalysts of unequalled versatility, capable of regio- and stereoselective functionalization of unactivated CH bonds by monooxygenation. Collectively they catalyze over 60 different biotransformations on structurally and functionally diverse organic molecules, including natural products, drugs, steroids, organic acids and other lipophilic molecules. This catalytic versatility and substrate range makes them likely candidates for application as potential biocatalysts for industrial chemistry. However, several aspects of the P450 catalytic cycle and other characteristics have limited their implementation to date in industry, including: their lability at elevated temperature, in the presence of solvents, and over lengthy incubation times; the typically low efficiency with which they metabolize non-natural substrates; and their lack of specificity for a single metabolic pathway. Protein engineering by rational design or directed evolution provides a way to engineer P450s for industrial use. Here we review the progress made to date toward engineering the properties of P450s, especially eukaryotic forms, for industrial application, and including the recent expansion of their catalytic repertoire to include non-natural reactions.

Molecular Basis of the Recognition of Cholesterol by Cytochrome P450 46A1 along the Major Access Tunnel

CYP46A1 is an important potential target for the treatment of Alzheimer's disease (AD), which is the most common neurodegenerative disease among older individuals. However, the binding mechanism between CYP46A1 and substrate cholesterol (CH) has not been clarified and will not be conducive to the research of relevant drug molecules. In this study, we integrated molecular docking, molecular dynamics (MD) simulations, and adaptive steered MD simulations to explore the recognition and binding mechanism of CYP46A1 with CH. Two key factors affecting the interaction between CH and CYP46A1 are determined: one is a hydrophobic cavity formed by seven hydrophobic residues (F80, Y109, L112, I222, W368, F371, and T475), which provides nonpolar interactions to stabilize CH, and the other is a hydrogen bond formed by H81 and CH, which ensures the binding direction of CH. In addition, the tunnel analysis results show that tunnel 2a is identified as the primary pathway of CH. The entry of CH induces tunnel 2e to close and tunnel w to open. Our results may provide effective clues for the design of drugs based on the substrate for AD and improve our understanding of the structure-function of CYP46A1.

Functional expression and regulation of eukaryotic cytochrome P450 enzymes in surrogate microbial cell factories

Cytochrome P450 (CYP) enzymes play crucial roles during the evolution and diversification of ancestral monocellular eukaryotes into multicellular eukaryotic organisms due to their essential functionalities including catalysis of housekeeping biochemical reactions, synthesis of diverse metabolites, detoxification of xenobiotics, and contribution to environmental adaptation. Eukaryotic CYPs with versatile functionalities are undeniably regarded as promising biocatalysts with great potential for biotechnological, pharmaceutical and chemical industry applications. Nevertheless, the modes of action and the challenges associated with these membrane-bound proteins have hampered the effective utilization of eukaryotic CYPs in a broader range. This review is focused on comprehensive and consolidated approaches to address the core challenges in heterologous expression of membrane-bound eukaryotic CYPs in different surrogate microbial cell factories, aiming to provide key insights for better studies and applications of diverse eukaryotic CYPs in the future. We also highlight the functional significance of the previously underrated cytochrome P450 reductases (CPRs) and provide a rational justification on the progression of CPR from auxiliary redox partner to function modulator in CYP catalysis.

Role of oxidative stress in the pathogenesis of COPD

Chronic inhalation of cigarette smoke is a prominent cause of chronic obstructive pulmonary disease (COPD) and provides an important source of exogenous oxidants. In addition, several inflammatory and structural cells are a source of endogenous oxidants in the lower airways of COPD patients, even in former smokers. This suggests that oxidants play a key role in the pathogenesis of COPD. This oxidative stress is counterbalanced by the protective effects of the various endogenous antioxidant defenses of the lower airways. A large amount of data from animal models and patients with COPD have shown that both the stable phase of the disease, and during exacerbations, have increased oxidative stress in the lower airways compared with age-matched smokers with normal lung function. Thus, counteracting the increased oxidative stress may produce clinical benefits in COPD patients. Smoking cessation is currently the most effective treatment of COPD patients and reduces oxidative stress in the lower airways. In addition, many drugs used to treat COPD have some antioxidant effects, however, it is still unclear if their clinical efficacy is related to pharmacological modulation of the oxidant/antioxidant balance. Several new antioxidant compounds are in development for the treatment of COPD.

Forty Years of Directed Evolution and its Continuously Evolving Technology Toolbox - A Review of the Patent Landscape

Generating functional protein variants with novel or improved characteristics has been a goal of the biotechnology industry and life sciences, for decades. Rational design and directed evolution are two major pathways to achieve the desired ends. Whilst rational protein design approach has made substantial progress, the idea of using a method based on cycles of mutagenesis and natural selection to develop novel binding proteins, enzymes and structures has attracted great attention. Laboratory evolution of proteins/enzymes requires new tools and analytical approaches to create genetic diversity and identifying variants with desired traits. In this pursuit, construction of sufficiently large libraries of target molecules to search for improved variants and the need for new protocols to alter the properties of target molecules has been a continuing challenge in the directed evolution experiments. This review will discuss the in vivo and in vitro gene diversification tools, library screening or selection approaches, and artificial intelligence/machine-learning-based strategies to mutagenesis developed in the last forty years to accelerate the natural process of evolution in creating new functional protein variants, optimization of microbial strains and transformation of enzymes into industrial machines. Analyzing patent position over these techniques and mechanisms also constitutes an integral and distinctive part of this review. The aim is to provide an up-to-date resource/technology toolbox for research-based and pharmaceutical companies to discover the boundaries of competitor's intellectual property (IP) portfolio, their freedom-to-operate in the relevant IP landscape, and the need for patent due diligence analysis to rule out whether use of a particular patented mutagenesis method, library screening/selection technique falls outside the safe harbor of experimental use exemption. While so doing, we have referred to some recent cases that emphasize the significance of selecting a suitable gene diversification strategy in directed evolution experiments. This article is protected by copyright. All rights reserved.

Redox metabolism for improving whole-cell P450-catalysed terpenoid biosynthesis

The growing preference for producing cytochrome P450-mediated natural products in microbial systems stems from the challenging nature of the organic chemistry approaches. The P450 enzymes are redox-dependent proteins, through which they source electrons from reducing cofactors to drive their activities. Widely researched in biochemistry, most of the previous studies have extensively utilised expensive cell-free assays to reveal mechanistic insights into P450 functionalities in presence of commercial redox partners. However, in the context of microbial bioproduction, the synergic activity of P450- reductase proteins in microbial systems have not been largely investigated. This is mainly due to limited knowledge about their mutual interactions in the context of complex systems. Hence, manipulating the redox potential for natural product synthesis in microbial chassis has been limited. As the potential of redox state as crucial regulator of P450 biocatalysis has been greatly underestimated by the scientific community, in this review, we re-emphasize their pivotal role in modulating the in vivo P450 activity through affecting the product profile and yield. Particularly, we discuss the applications of widely used in vivo redox engineering methodologies for natural product synthesis to provide further suggestions for patterning on P450-based terpenoids production in microbial platforms.

Mapping protein-protein interactions in homodimeric CYP102A1 by crosslinking and mass spectrometry

Covalent crosslinking and mass spectrometry techniques hold great potential in the study of multiprotein complexes, but a major challenge is the inability to differentiate intra- and inter- protein crosslinks in homomeric complexes. In the current study we use CYP102A1, a well-characterized homodimeric P450, to examine a subtractive method that utilizes limited crosslinking with disuccinimidyl dibutyric urea (DSBU) and isolation of the monomer, in addition to the crosslinked dimer, to identify inter-monomer crosslinks. The utility of this approach was examined with the use of MS-cleavable crosslinker DSBU and recently published cryo-EM based structures of the CYP102A1 homodimer. Of the 31 unique crosslinks found, 26 could be fit to the reported structures whereas 5 exceeded the spatial constraints. Not only did these crosslinks validate the cryo-EM structure, they point to new conformations of CYP102A1 that bring the flavins in closer proximity to the heme.

Amphipol-facilitated elucidation of the functional tetrameric complex of full-length cytochrome P450 CYP2B4 and NADPH-cytochrome P450 oxidoreductase

Interactions of membrane-bound mammalian cytochromes P450 (CYPs) with NADPH-cytochrome P450 oxidoreductase (POR), which are required for metabolism of xenobiotics, are facilitated by membrane lipids. A variety of membrane mimetics, such as phospholipid liposomes and nanodiscs, have been used to simulate the membrane to form catalytically active CYP:POR complexes. However, the exact mechanism(s) of these interactions are unclear because of the absence of structural information of full-length mammalian CYP:POR complexes in membranes. Herein, we report the use of amphipols (APols) to form a fully functional, soluble, homogeneous preparation of full-length CYP:POR complexes amenable to biochemical and structural study. Incorporation of CYP2B4 and POR into APols resulted in a CYP2B4:POR complex with a stoichiometry of 1:1, which was fully functional in demethylating benzphetamine at a turnover rate of 37.7 ± 2.2 min^-1, with a coupling efficiency of 40%. Interestingly, the stable complex had a molecular weight (M_w) of 338 ± 22 kDa determined by multiangle light scattering, suggestive of a tetrameric complex of 2CYP2B4:2POR embedded in one APol nanoparticle. Moreover, negative stain electron microscopy (EM) validated the homogeneity of the complex and allowed us to generate a three-dimensional EM map and model consistent with the tetramer observed in solution. This first report of the full-length mammalian CYP:POR complex by transmission EM not only reveals the architecture that facilitates electron transfer but also highlights a potential use of APols in biochemical and structural studies of functional CYP complexes with redox partners.

Conformational Landscape of Cytochrome P450 Reductase Interactions

Oxidative reactions catalyzed by Cytochrome P450 enzymes (CYPs), which constitute the most relevant group of drug-metabolizing enzymes, are enabled by their redox partner Cytochrome P450 reductase (CPR). Both proteins are anchored to the membrane of the endoplasmic reticulum and the CPR undergoes a conformational change in order to interact with the respective CYP and transfer electrons. Here, we conducted over 22 microseconds of molecular dynamics (MD) simulations in combination with protein-protein docking to investigate the conformational changes necessary for the formation of the CPR-CYP complex. While some structural features of the CPR and the CPR-CYP2D6 complex that we highlighted confirmed previous observations, our simulations revealed additional mechanisms for the conformational transition of the CPR. Unbiased simulations exposed a movement of the whole protein relative to the  membrane, potentially to facilitate interactions with its diverse set of redox partners. Further, we present a structural mechanism for the susceptibility of the CPR to different redox states based on the flip of a glycine residue disrupting the local interaction network that maintains inter-domain proximity. Simulations of the CPR-CYP2D6 complex pointed toward an additional interaction surface of the FAD domain and the proximal side of CYP2D6. Altogether, this study provides novel structural insight into the mechanism of CPR-CYP interactions and underlying conformational changes, improving our understanding of this complex machinery Cytochrome P450 reductase; CPR; conformational; dynamicsrelevant for drug metabolism.