Structure of a cytochrome P450-redox partner electron-transfer complex

Biosynthesis of eriodictyol in citrus waster by endowing P450BM3 activity of naringenin hydroxylation

The flavonoid naringenin is abundantly present in pomelo peels, and the unprocessed naringenin in wastes is not friendly for the environment once discarded directly. Fortunately, the hydroxylated product of eriodictyol from naringenin exhibits remarkable antioxidant and anticancer properties. The P450s was suggested promising for the bioconversion of the flavonoids, but less naturally existed P450s show hydroxylation activity to C3' of the naringenin. By well analyzing the catalytic mechanism and the conformations of the naringenin in P450, we proposed that the intermediate Cmpd I ((porphyrin)Fe = O) is more reasonable as key conformation for the hydrolyzation, and the distance between C3'/C5' of naringenin to the O atom of CmpdI determines the hydroxylating activity for the naringenin. Thus, the "flying kite model" that gradually drags the C-H bond of the substrate to the O atom of CmpdI was put forward for rational design. With ab initio design, we successfully endowed the self-sufficient P450-BM3 hydroxylic activity to naringenin and obtained mutant M5-5, with kcat, Km, and kcat/Km values of 230.45 min-1, 310.48 µM, and 0.742 min-1 µM-1, respectively. Furthermore, the mutant M4186 was screened with kcat/Km of 4.28-fold highly improved than the reported M13. The M4186 also exhibited 62.57% yield of eriodictyol, more suitable for the industrial application. This study provided a theoretical guide for the rational design of P450s to the nonnative compounds. KEY POINTS: •The compound I is proposed as the starting point for the rational design of the P450BM3 •"Flying kite model" is proposed based on the distance between O of Cmpd I and C3'/C5' of naringenin •Mutant M15-5 with 1.6-fold of activity than M13 was obtained by ab initio modification.

Understanding the stability of a plastic-degrading Rieske iron oxidoreductase system

Rieske oxygenases (ROs) are a diverse metalloenzyme class with growing potential in bioconversion and synthetic applications. We postulated that ROs are nonetheless underutilized because they are unstable. Terephthalate dioxygenase (TPADO PDB ID 7Q05) is a structurally characterized heterohexameric α3β3 RO that, with its cognate reductase (TPARED), catalyzes the first intracellular step of bacterial polyethylene terephthalate plastic bioconversion. Here, we showed that the heterologously expressed TPADO/TPARED system exhibits only ~300 total turnovers at its optimal pH and temperature. We investigated the thermal stability of the system and the unfolding pathway of TPADO through a combination of biochemical and biophysical approaches. The system's activity is thermally limited by a melting temperature (Tm) of 39.9°C for the monomeric TPARED, while the independent Tm of TPADO is 50.8°C. Differential scanning calorimetry revealed a two-step thermal decomposition pathway for TPADO with Tm values of 47.6 and 58.0°C (ΔH = 210 and 509 kcal mol-1, respectively) for each step. Temperature-dependent small-angle x-ray scattering and dynamic light scattering both detected heat-induced dissociation of TPADO subunits at 53.8°C, followed by higher-temperature loss of tertiary structure that coincided with protein aggregation. The computed enthalpies of dissociation for the monomer interfaces were most congruent with a decomposition pathway initiated by β-β interface dissociation, a pattern predicted to be widespread in ROs. As a strategy for enhancing TPADO stability, we propose prioritizing the re-engineering of the β subunit interfaces, with subsequent targeted improvements of the subunits.

Calcitriol Treatment Is Safe and Increases Frataxin Levels in Friedreich Ataxia Patients

BACKGROUND

Calcitriol, the active form of vitamin D (also known as 1,25-dihydroxycholecalciferol), improves the phenotype and increases frataxin levels in cell models of Friedreich ataxia (FRDA).

OBJECTIVES

Based on these results, we aimed measuring the effects of a calcitriol dose of 0.25 mcg/24h in the neurological function and frataxin levels when administered to FRDA patients for a year.

METHODS

20 FRDA patients where recluted and 15 patients completed the treatment for a year. Evaluations of neurological function changes (SARA scale, 9-HPT, 8-MWT, PATA test) and quality of life (Barthel Scale and Short Form (36) Health Survey [SF-36] quality of life questionnaire) were performed. Frataxin amounts were measured in isolated platelets obtained from these FRDA patients, from heterozygous FRDA carriers (relatives of the FA patients) and from non-heterozygous sex and age matched controls.

RESULTS

Although the patients did not experience any observable neurological improvement, there was a statistically significant increase in frataxin levels from initial values, 5.5 to 7.0 pg/μg after 12 months. Differences in frataxin levels referred to total protein levels were observed among sex- and age-matched controls (18.1 pg/μg), relative controls (10.1 pg/μg), and FRDA patients (5.7 pg/μg). The treatment was well tolerated by most patients, and only some of them experienced minor adverse effects at the beginning of the trial.

CONCLUSIONS

Calcitriol dosage used (0.25 mcg/24 h) is safe for FRDA patients, and it increases frataxin levels. We cannot rule out that higher doses administered longer could yield neurological benefits. © 2024 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

On the engineering of reductase-based-monooxygenase activity in CYP450 peroxygenases

Recent bioengineering of CYP450OleT shows that peroxide-based CYP450OleT can be converted to a reductase-based self-sufficient enzyme, which is capable of showing efficient hydroxylation and decarboxylation activity for a wide range of substrates. The so-generated enzyme creates several mechanistic puzzles: (A) as CYP450 peroxygenases lack the conventional acid-alcohol pair, what is the source of two protons that are required to create the ultimate oxidant Cpd I? (B) Why is it only CYP450OleT that shows the reductase-based activity but no other CYP members? The present study provides a mechanistic solution to these puzzles using comprehensive MD simulations and hybrid QM/MM calculations. We show that the fusion of the reductase domain to the heme-binding domain triggers significant conformational rearrangement, which is gated by the propionate side chain, which constitutes a new water aqueduct via the carboxylate end of the substrate that ultimately participates in Cpd I formation. Importantly, such well-synchronized choreographies are controlled by remotely located Tyr359, which senses the fusion of reductase and communicates to the heme domain via non-covalent interactions. These findings provide crucial insights and a broader perspective which enables us to make a verifiable prediction: thus, the catalytic activity is not only limited to the first or second catalytic shell of an enzyme. Furthermore, it is predicted that reinstatement of tyrosine at a similar position in other members of CYP450 peroxygenases can convert these enzymes to reductase-based monooxygenases.

Structure, dynamics, and redox reactivity of an all-purpose flavodoxin

The flavodoxin of Rhodopseudomonas palustris CGA009 (Rp9Fld) supplies highly reducing equivalents to crucial enzymes such as hydrogenase, especially when the organism is iron-restricted. By acquiring those electrons from photodriven electron flow via the bifurcating electron transfer flavoprotein, Rp9Fld provides solar power to vital metabolic processes. To understand Rp9Fld's ability to work with diverse partners, we solved its crystal structure. We observed the canonical flavodoxin (Fld) fold and features common to other long-chain Flds but not all the surface loops thought to recognize partner proteins. Moreover, some of the loops display alternative structures and dynamics. To advance studies of protein-protein associations and conformational consequences, we assigned the 19F NMR signals of all five tyrosines (Tyrs). Our electrochemical measurements show that incorporation of 3-19F-Tyr in place of Tyr has only a modest effect on Rp9Fld's redox properties even though Tyrs flank the flavin on both sides. Meanwhile, the 19F probes demonstrate the expected paramagnetic effect, with signals from nearby Tyrs becoming broadened beyond detection when the flavin semiquinone is formed. However, the temperature dependencies of chemical shifts and linewidths reveal dynamics affecting loops close to the flavin and regions that bind to partners in a variety of systems. These coincide with patterns of amino acid type conservation but not retention of specific residues, arguing against detailed specificity with respect to partners. We propose that the loops surrounding the flavin adopt altered conformations upon binding to partners and may even participate actively in electron transfer.

Structural and spectroscopic characterization of RufO indicates a new biological role in rufomycin biosynthesis

Rufomycins constitute a class of cyclic heptapeptides isolated from actinomycetes. They are secondary metabolites that show promising treatment against Mycobacterium tuberculosis infections by inhibiting a novel drug target. Several nonproteinogenic amino acids are integrated into rufomycins, including a conserved 3-nitro-tyrosine. RufO, a cytochrome P450 (CYP)-like enzyme, was proposed to catalyze the formation of 3-nitro-tyrosine in the presence of O2 and NO. To define its biological function, the interaction between RufO and the proposed substrate tyrosine is investigated using various spectroscopic methods that are sensitive to the structural change of a heme center. However, a low- to high-spin state transition and a dramatic increase in the redox potential that are commonly found in CYPs upon ligand binding have not been observed. Furthermore, a 1.89-Å crystal structure of RufO shows that the enzyme has flexible surface regions, a wide-open substrate access tunnel, and the heme center is largely exposed to solvent. Comparison with a closely related nitrating CYP reveals a spacious and hydrophobic distal pocket in RufO, which is incapable of stabilizing a free amino acid. Molecular docking validates the experimental data and proposes a possible substrate. Collectively, our results disfavor tyrosine as the substrate of RufO and point to the possibility that the nitration occurs during or after the assembly of the peptides. This study indicates a new function of the unique nitrating enzyme and provides insights into the biosynthesis of nonribosomal peptides.

The Versatile Biocatalyst of Cytochrome P450 CYP102A1: Structure, Function, and Engineering

Wild-type cytochrome P450 CYP102A1 from Bacillus megaterium is a highly efficient monooxygenase for the oxidation of long-chain fatty acids. The unique features of CYP102A1, such as high catalytic activity, expression yield, regio- and stereoselectivity, and self-sufficiency in electron transfer as a fusion protein, afford the requirements for an ideal biocatalyst. In the past three decades, remarkable progress has been made in engineering CYP102A1 for applications in drug discovery, biosynthesis, and biotechnology. The repertoire of engineered CYP102A1 variants has grown tremendously, whereas the substrate repertoire is avalanched to encompass alkanes, alkenes, aromatics, organic solvents, pharmaceuticals, drugs, and many more. In this article, we highlight the major advances in the past five years in our understanding of the structure and function of CYP102A1 and the methodologies used to engineer CYP102A1 for novel applications. The objective is to provide a succinct review of the latest developments with reference to the body of CYP102A1-related literature.

Structures and Functions of Human Placental Aromatase and Steroid Sulfatase, Two Key Enzymes in Estrogen Biosynthesis

Cytochrome P450 aromatase (AROM) and steroid sulfatase (STS) are the two key enzymes for the biosynthesis of estrogens in human, and maintenance of the critical balance between androgens and estrogens. Human AROM, an integral membrane protein of the endoplasmic reticulum, is a member of the cytochrome P450 superfamily. It is the only enzyme to catalyze the conversion of androgens with non-aromatic A-rings to estrogens characterized by the aromatic A-ring. Human STS, also an integral membrane protein of the endoplasmic reticulum, is a Ca²⁺-dependent enzyme that catalyzes the hydrolysis of sulfate esters of estrone and dehydroepiandrosterone to the unconjugated steroids, the precursors of the most potent forms of estrogens and androgens, namely, 17β-estradiol, 16α,17β-estriol, testosterone and dihydrotestosterone. Expression of these steroidogenic enzymes locally within organs and tissues of the endocrine, reproductive, and central nervous systems is the key for maintaining high levels of the reproductive steroids. The enzymes have been drug targets for the prevention and treatment of diseases associated with steroid hormone excesses, especially in breast, endometrial and prostate malignancies. Both enzymes have been the subjects of vigorous research for the past six decades. In this article, we review the important findings on their structure-function relationships, specifically, the work that began with unravelling of the closely guarded secrets, namely, the 3-D structures, active sites, mechanisms of action, origins of substrate specificity and the basis of membrane integration. Remarkably, these studies were conducted on the enzymes purified in their pristine forms from human placenta, the discarded and their most abundant source. The purification, assay, crystallization, and structure determination methodologies are described. Also reviewed are their functional quaternary organizations, post-translational modifications and the advancements made in the structure-guided inhibitor design efforts. Outstanding questions that still remain open are summarized in closing.

Improved Bioproduction of the Nylon 12 Monomer by Combining the Directed Evolution of P450 and Enhancing Heme Synthesis

The nylon 12 (PA12) monomer ω-aminododecanoic acid (ω-AmDDA) could be synthesized from lauric acid (DDA) through multi-enzyme cascade transformation using engineered E. coli, with the P450 catalyzing terminal hydroxylation of DDA as a rate-limiting enzyme. Its activity is jointly determined by the heme domain and the reductase domain. To obtain a P450 mutant with higher activity, directed evolution was conducted using a colorimetric high-throughput screening (HTS) system with DDA as the real substrate. After two rounds of directed evolution, a positive double-site mutant (R14R/D629G) with 90.3% higher activity was obtained. Molecular docking analysis, kinetic parameter determination and protein electrophoresis suggested the improved soluble expression of P450 resulting from the synonymous mutation near the N-terminus and the shortened distance of the electron transfer between FMN and FAD caused by D629G mutation as the major reasons for activity improvement. The significantly increased k_cat and unchanged K_m provided further evidence for the increase in electron transfer efficiency. Considering the important role of heme in P450, its supply was strengthened by the metabolic engineering of the heme synthesis pathway. By combining P450-directed evolution and enhancing heme synthesis, 2.02 ± 0.03 g/L of ω-AmDDA was produced from 10 mM DDA, with a yield of 93.6%.

Axial N Ligand-Modulated Ultrahigh Activity and Selectivity Hyperoxide Activation over Single-Atoms Nanozymes

Learning and studying the structure-activity relationship in the bio-enzymes is conducive to the design of nanozymes for energy and environmental application. Herein, Fe single-atom nanozymes (Fe-SANs) with Fe-N5 site, inspired by the structure of cytochromes P450 (CYPs), are developed and characterized. Similar to the CYPs, the hyperoxide can activate the Fe(III) center of Fe-SANs to generate Fe(IV)O intermediately, which can transfer oxygen to the substrate with ultrafast speed. Particularly, using the peroxymonosulfate (PMS)-activated Fe-SANs to oxidize sulfamethoxazole, a typical antibiotic contaminant, as the model hyperoxides activation reaction, the excellent activity within 284 min-1 g-1 (catalyst) mmol-1 (PMS) oxidation rate and 91.6% selectivity to the Fe(IV)O intermediate oxidation are demonstrated. More importantly, instead of promoting PMS adsorption, the axial N ligand modulates the electron structure of FeN5 SANs for the lower reaction energy barrier and promotes electron transfer to PMS to produce Fe(IV)O intermediate with high selectivity. The highlight of the axial N coordination in the nanozymes in this work provides deep insight to guide the design and development of nanozymes nearly to the bio-enzyme with excellent activity and selectivity.

Mechanisms of Electron Transfer Rate Modulations in Cytochrome P450 BM3

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

Unravelling the role of transient redox partner complexes in P450 electron transfer mechanics

The molecular evolution of cytochromes P450 and associated redox-driven oxidative catalysis remains a mystery in biology. It is widely believed that sterol 14α-demethylase (CYP51), an essential enzyme of sterol biosynthesis, is the ancestor of the whole P450 superfamily given its conservation across species in different biological kingdoms. Herein we have utilized X-ray crystallography, molecular dynamics simulations, phylogenetics and electron transfer measurements to interrogate the nature of P450-redox partner binding using the naturally occurring fusion protein, CYP51-ferredoxin found in the sterol-producing bacterium Methylococcus capsulatus. Our data advocates that the electron transfer mechanics in the M. capsulatus CYP51-ferredoxin fusion protein involves an ensemble of ferredoxin molecules in various orientations and the interactions are transient. Close proximity of ferredoxin, however, is required to complete the substrate-induced large-scale structural switch in the P450 domain that enables proton-coupled electron transfer and subsequent oxygen scission and catalysis. These results have fundamental implications regarding the early evolution of electron transfer proteins and for the redox reactions in the early steps of sterol biosynthesis. They also shed new light on redox protein mechanics and the subsequent diversification of the P450 electron transfer machinery in nature.

Confrontation of AlphaFold models with experimental structures enlightens conformational dynamics supporting CYP102A1 functions

Conformational dynamics plays a critical role for the function of multidomain electron transfer complexes. While crystallographic or NMR approaches allow detailed insight into structures, lower resolution methods like cryo-electron microscopy can provide more information on dynamics. In silico structure modelling using AlphaFold was recently successfully extended to the prediction of protein complexes but its capability to address large conformational changes involved in catalysis remained obscure. We used bacterial CYP102A1 monooxygenase homodimer as a test case to design a competitive modelling approach (CMA) for assessing alternate conformations of multi-domain complexes. Predictions were confronted with published crystallographic and cryo-EM data, evidencing consistencies but also permitting some reinterpretation of experimental data. Structural determinants stabilising the new type of domain connectivity evidenced in this bacterial self-sufficient monooxygenase were analysed by CMA and used for in silico retro-engineering applied to its eukaryotic bi-component counterparts.

Tartary buckwheat FtF3'H1 as a metabolic branch switch to increase anthocyanin content in transgenic plant

Tartary buckwheat (TB) is a pseudocereal rich in flavonoids, mainly including flavonols and anthocyanins. The flavonoid 3'-hydroxylase (F3'H) is a key enzyme in flavonoid biosynthesis and is encoded by two copies in TB genome. However, its biological function and effects on flavonol and anthocyanin synthesis in TB have not been well validated yet. In this study, we cloned the full-length FtF3'H1 gene highly expressed in all tissues (compared with FtF3'H2) according to TB flowering transcriptome data. The corresponding FtF3'H1 protein contains 534 amino acids with the molecular properties of the typical plant F3'H and belongs to the CYP75B family. During the flowering stage, the FtF3'H1 expression was highest in flowers, and its expression pattern showed a significant and positive correlation with the total flavonoids (R ² > 0.95). The overexpression of FtF3'H1 in Arabidopsis thaliana, Nicotiana tabacum and TB hairy roots resulted in a significant increase in anthocyanin contents (p < 0.05) but a decrease in rutin (p < 0.05). The average anthocyanin contents were 2.94 mg/g (fresh weight, FW) in A. thaliana (about 135% increase), 1.18 mg/g (FW) in tobacco (about 17% increase), and 1.56 mg/g (FW) TB hairy roots (about 44% increase), and the rutin contents were dropped to about 53.85, 14.99, 46.31%, respectively. However, the expression of genes involved in anthocyanin (DFRs and ANSs) and flavonol (FLSs) synthesis pathways were significantly upregulated (p < 0.05). In particular, the expression level of DFR, a key enzyme that enters the anthocyanin branch, was upregulated thousand-fold in A. thaliana and in N. tabacum. These results might be attributed to FtF3'H1 protein with a higher substrate preference for anthocyanin synthesis substrates. Altogether, we identified the basic biochemical activity of FtF3'H1 in vivo and investigated its involvement in anthocyanin and flavonol metabolism in plant.

A disease-associated missense mutation in CYP4F3 affects the metabolism of leukotriene B4 via disruption of electron transfer

BACKGROUND

Cytochrome P450 4F3 (CYP4F3) is an ω-hydroxylase that oxidizes leukotriene B4 (LTB4), prostaglandins, and fatty acid epoxides. LTB4 is synthesized by leukocytes and acts as a chemoattractant for neutrophils, making it an essential component of the innate immune system. Recently, involvement of the LTB4 pathway was reported in various immunological disorders such as asthma, arthritis, and inflammatory bowel disease. We report a 26-year-old female with a complex immune phenotype, mainly marked by exhaustion, muscle weakness, and inflammation-related conditions. The molecular cause is unknown, and symptoms have been aggravating over the years.

METHODS

Whole exome sequencing was performed and validated; flow cytometry and enzyme-linked immunosorbent assay were used to describe patient's phenotype. Function and impact of the mutation were investigated using molecular analysis: co-immunoprecipitation, western blot, and enzyme-linked immunosorbent assay. Capillary electrophoresis with ultraviolet detection was used to detect LTB4 and its metabolite and in silico modelling provided structural information.

RESULTS

We present the first report of a patient with a heterozygous de novo missense mutation c.C1123 > G;p.L375V in CYP4F3 that severely impairs its activity by 50% (P < 0.0001), leading to reduced metabolization of the pro-inflammatory LTB4. Systemic LTB4 levels (1034.0 ± 75.9 pg/mL) are significantly increased compared with healthy subjects (305.6 ± 57.0 pg/mL, P < 0.001), and immune phenotyping shows increased total CD19+ CD27- naive B cells (25%) and decreased total CD19+ CD27+ IgD- switched memory B cells (19%). The mutant CYP4F3 protein is stable and binding with its electron donors POR and Cytb5 is unaffected (P > 0.9 for both co-immunoprecipitation with POR and Cytb5). In silico modelling of CYP4F3 in complex with POR and Cytb5 suggests that the loss of catalytic activity of the mutant CYP4F3 is explained by a disruption of an α-helix that is crucial for the electron shuffling between the electron carriers and CYP4F3. Interestingly, zileuton still inhibits ex vivo LTB4 production in patient's whole blood to 2% of control (P < 0.0001), while montelukast and fluticasone do not (99% and 114% of control, respectively).

CONCLUSIONS

A point mutation in the catalytic domain of CYP4F3 is associated with high leukotriene B4 plasma levels and features of a more naive adaptive immune response. Our data provide evidence for the pathogenicity of the CYP4F3 variant as a cause for the observed clinical features in the patient. Inhibitors of the LTB4 pathway such as zileuton show promising effects in blocking LTB4 production and might be used as a future treatment strategy.

Designing Artificial Metalloenzymes by Tuning of the Environment beyond the Primary Coordination Sphere

Metalloenzymes catalyze a variety of reactions using a limited number of natural amino acids and metallocofactors. Therefore, the environment beyond the primary coordination sphere must play an important role in both conferring and tuning their phenomenal catalytic properties, enabling active sites with otherwise similar primary coordination environments to perform a diverse array of biological functions. However, since the interactions beyond the primary coordination sphere are numerous and weak, it has been difficult to pinpoint structural features responsible for the tuning of activities of native enzymes. Designing artificial metalloenzymes (ArMs) offers an excellent basis to elucidate the roles of these interactions and to further develop practical biological catalysts. In this review, we highlight how the secondary coordination spheres of ArMs influence metal binding and catalysis, with particular focus on the use of native protein scaffolds as templates for the design of ArMs by either rational design aided by computational modeling, directed evolution, or a combination of both approaches. In describing successes in designing heme, nonheme Fe, and Cu metalloenzymes, heteronuclear metalloenzymes containing heme, and those ArMs containing other metal centers (including those with non-native metal ions and metallocofactors), we have summarized insights gained on how careful controls of the interactions in the secondary coordination sphere, including hydrophobic and hydrogen bonding interactions, allow the generation and tuning of these respective systems to approach, rival, and, in a few cases, exceed those of native enzymes. We have also provided an outlook on the remaining challenges in the field and future directions that will allow for a deeper understanding of the secondary coordination sphere a deeper understanding of the secondary coordintion sphere to be gained, and in turn to guide the design of a broader and more efficient variety of ArMs.

Modulating the Coupling Efficiency of P450 BM3 by Controlling Water Diffusion through Access Tunnel Engineering

Cytochromes P450 have gained much interest for their broad substrate scope in the catalysis of oxidation reactions for pharmaceuticals, plastics, and hormones. However, achieving high coupling efficiency by the engineering of P450s is still a big challenge. The presence of extra water around the active site is deemed to be related to uncoupling. In this study, the access tunnels of P450 BM3 from Bacillus megaterium are engineered to control water access from bulk solvent to the active site. Nine residues located in tunnels are investigated by site-saturation mutagenesis to reduce water diffusion, thereby improving the coupling efficiency. The recombined variant N319L/T411V/T436A shows improved coupling efficiency (from 31.2 % to 52.6 %). Tunnel polarity analysis and molecular dynamics simulation further indicate that reduced water molecules around the active site lead to higher coupling efficiency. Overall, this study provides valuable insight on improving coupling efficiency by controlling water diffusion through tunnel engineering.

Valorization of Small Alkanes by Biocatalytic Oxyfunctionalization

The oxidation of alkanes into valuable chemical products is a vital reaction in organic synthesis. This reaction, however, is challenging, owing to the inertness of C-H bonds. Transition metal catalysts for C-H functionalization are frequently explored. Despite chemical alternatives, nature has also evolved powerful oxidative enzymes (e. g., methane monooxygenases, cytochrome P450 oxygenases, peroxygenases) that are capable of transforming C-H bonds under very mild conditions, with only the use of molecular oxygen or hydrogen peroxide as electron acceptors. Although progress in alkane oxidation has been reviewed extensively, little attention has been paid to small alkane oxidation. The latter holds great potential for the manufacture of chemicals. This Minireview provides a concise overview of the most relevant enzyme classes capable of small alkanes (C_<6 ) oxyfunctionalization, describes the essentials of the catalytic mechanisms, and critically outlines the current state-of-the-art in preparative applications.

Radiothermometric Study of the Effect of Amino Acid Mutation on the Characteristics of the Enzymatic System

The radiothermometry (RTM) study of a cytochrome-containing system (CYP102 A1) has been conducted in order to demonstrate the applicability of RTM for monitoring changes in the functional activity of an enzyme in case of its point mutation. The study has been performed with the example of the wild-type cytochrome (WT) and its mutant type A264K. CYP102 A1 is a nanoscale protein-enzymatic system of about 10 nm in size. RTM uses a radio detector and can record the corresponding brightness temperature (T_br) of the nanoscale enzyme solution within the 3.4–4.2 GHz frequency range during enzyme functioning. It was found that the enzymatic reaction during the lauric acid hydroxylation at the wild-type CYP102 A1 (WT) concentration of ~10⁻⁹ M is accompanied by T_br fluctuations of ~0.5–1 °C. At the same time, no T_br fluctuations are observed for the mutated forms of the enzyme CYP102 A1 (A264K), where one amino acid was replaced. We know that the activity of CYP102 A1 (WT) is ~4 orders of magnitude higher than that of CYP102 A1 (A264K). We therefore concluded that the disappearance of the fluctuation of T_br CYP102 A1 (A264K) is associated with a decrease in the activity of the enzyme. This effect can be used to develop new methods for testing the activity of the enzyme that do not require additional labels and expensive equipment, in comparison with calorimetry and spectral methods. The RTM is beginning to find application in the diagnosis of oncological diseases and for the analysis of biochemical processes.

Molecular Lego of Human Cytochrome P450: The Key Role of Heme Domain Flexibility for the Activity of the Chimeric Proteins

The cytochrome P450 superfamily are heme-thiolate enzymes able to carry out monooxygenase reactions. Several studies have demonstrated the feasibility of using a soluble bacterial reductase from Bacillus megaterium, BMR, as an artificial electron transfer partner fused to the human P450 domain in a single polypeptide chain in an approach known as ‘molecular Lego’. The 3A4-BMR chimera has been deeply characterized biochemically for its activity, coupling efficiency, and flexibility by many different biophysical techniques leading to the conclusion that an extension of five glycines in the loop that connects the two domains improves all the catalytic parameters due to improved flexibility of the system. In this work, we extend the characterization of 3A4-BMR chimeras using differential scanning calorimetry to evaluate stabilizing role of BMR. We apply the ‘molecular Lego’ approach also to CYP19A1 (aromatase) and the data show that the activity of the chimeras is very low (<0.003 min−1) for all the constructs tested with a different linker loop length: ARO-BMR, ARO-BMR-3GLY, and ARO-BMR-5GLY. Nevertheless, the fusion to BMR shows a remarkable effect on thermal stability studied by differential scanning calorimetry as indicated by the increase in Tonset by 10 °C and the presence of a cooperative unfolding process driven by the BMR protein domain. Previously characterized 3A4-BMR constructs show the same behavior of ARO-BMR constructs in terms of thermal stabilization but a higher activity as a function of the loop length. A comparison of the ARO-BMR system to 3A4-BMR indicates that the design of each P450-BMR chimera should be carefully evaluated not only in terms of electron transfer, but also for the biophysical constraints that cannot always be overcome by chimerization.

Functional and structural insight into the flexibility of cytochrome P450 reductases from Sorghum bicolor and its implications for lignin composition

Plant NADPH-dependent cytochrome P450 reductase (CPR) is a multidomain enzyme that donates electrons for hydroxylation reactions catalyzed by class II cytochrome P450 monooxygenases involved in the synthesis of many primary and secondary metabolites. These P450 enzymes include trans-cinnamate-4-hydroxylase, p-coumarate-3′-hydroxylase, and ferulate-5-hydroxylase involved in monolignol biosynthesis. Because of its role in monolignol biosynthesis, alterations in CPR activity could change the composition and overall output of lignin. Therefore, to understand the structure and function of three CPR subunits from sorghum, recombinant subunits SbCPR2a, SbCPR2b, and SbCPR2c were subjected to X-ray crystallography and kinetic assays. Steady-state kinetic analyses demonstrated that all three CPR subunits supported the oxidation reactions catalyzed by SbC4H1 (CYP73A33) and SbC3′H (CYP98A1). Furthermore, comparing the SbCPR2b structure with the well-investigated CPRs from mammals enabled us to identify critical residues of functional importance and suggested that the plant flavin mononucleotide–binding domain might be more flexible than mammalian homologs. In addition, the elucidated structure of SbCPR2b included the first observation of NADP⁺ in a native CPR. Overall, we conclude that the connecting domain of SbCPR2, especially its hinge region, could serve as a target to alter biomass composition in bioenergy and forage sorghums through protein engineering.

Prediction of Monomeric and Dimeric Structures of CYP102A1 Using AlphaFold2 and AlphaFold Multimer and Assessment of Point Mutation Effect on the Efficiency of Intra- and Interprotein Electron Transfer

The three-dimensional structure of monomers and homodimers of CYP102A1/WT (wild-type) proteins and their A83F and A83I mutant forms was predicted using the AlphaFold2 (AF2) and AlphaFold Multimer (AFMultimer) programs, which were compared with the rate constants of hydroxylation reactions of these enzyme forms to determine the efficiency of intra- and interprotein electron transport in the CYP102A1 hydroxylase system. The electron transfer rate constants (k_et), which determine the rate of indole hydroxylation by the CYP102A1 system, were calculated based on the distances (R) between donor-acceptor prosthetic groups (PG) FAD→FMN→HEME of these proteins using factor β, which describes an exponential decay from R the speed of electron transport (ET) according to the tunnelling mechanism. It was shown that the structure of monomers in the homodimer, calculated using the AlpfaFold Multimer program, is in good agreement with the experimental structures of globular domains (HEME-, FMN-, and FAD-domains) in CYP102A1/WT obtained by X-ray structural analysis, and the structure of isolated monomers predicted in AF2 does not coincide with the structure of monomers in the homodimer, although a high level of similarity in individual domains remains. The structures of monomers and homodimers of A83F and A83I mutants were also calculated, and their structures were compared with the wild-type protein. Significant differences in the structure of all isolated monomers with respect to the structures of monomers in homodimers were also found for them, and at the same time, insignificant differences were revealed for all homodimers. Comparative analysis for CYP102A1/WT between the calculated intra- and interprotein distances FAD→FMN→HEME and the rate constants of hydroxylation in these proteins showed that the distance between prosthetic groups both in the monomer and in the dimer allows the implementation of electron transfer between PGs, which is consistent with experimental literature data about k_cat. For the mutant form of monomer A83I, an increase in the distance between PGs was obtained, which can restrict electron transportation compared to WT; however, for the dimer of this protein, a decrease in the distance between PGs was observed compared to the WT form, which can lead to an increase in the electron transfer rate constant and, accordingly, k_cat. For the monomer and homodimer of the A83F mutant, the calculations showed an increase in the distance between the PGs compared to the WT form, which should have led to a decrease in the electron transfer rate, but at the same time, for the homodimer, the approach of the aromatic group F262 with heme can speed up transportation for this form and, accordingly, the rate of hydroxylation.

Updating the Paradigm: Redox Partner Binding and Conformational Dynamics in Cytochromes P450

This Account summarizes recent findings centered on the role that redox partner binding, allostery, and conformational dynamics plays in cytochrome P450 proton coupled electron transfer. P450s are one of Nature's largest enzyme families and it is not uncommon to find a P450 wherever substrate oxidation is required in the formation of essential molecules critical to the life of the organism or in xenobiotic detoxification. P450s can operate on a remarkably large range of substrates from the very small to the very large, yet the overall P450 three-dimensional structure is conserved. Given this conservation of structure, it is generally assumed that the basic catalytic mechanism is conserved. In nearly all P450s, the O₂ O-O bond must be cleaved heterolytically enabling one oxygen atom, the distal oxygen, to depart as water and leave behind a heme iron-linked O atom as the powerful oxidant that is used to activate the nearby substrate. For this process to proceed efficiently, externally supplied electrons and protons are required. Two protons must be added to the departing O atom while an electron is transferred from a redox partner that typically contains either a Fe₂S₂ or FMN redox center. The paradigm P450 used to unravel the details of these mechanisms has been the bacterial CYP101A1 or P450cam. P450cam is specific for its own Fe₂S₂ redox partner, putidaredoxin or Pdx, and it has long been postulated that Pdx plays an effector/allosteric role by possibly switching P450cam to an active conformation. Crystal structures, spectroscopic data, and direct binding experiments of the P450cam-Pdx complex provide some answers. Pdx shifts the conformation of P450cam to a more open state, a transition that is postulated to trigger the proton relay network required for O₂ activation. An essential part of this proton relay network is a highly conserved Asp (sometimes Glu) that is known to be critical for activity in a number of P450s. How this Asp and proton delivery networks are connected to redox partner binding is quite simple. In the closed state, this Asp is tied down by salt bridges, but these salt bridges are ruptured when Pdx binds, leaving the Asp free to serve its role in proton transfer. An alternative hypothesis suggests that a specific proton relay network is not really necessary. In this scenario, the Asp plays a structural role in the open/close transition and merely opening the active site access channel is sufficient to enable solvent protons in for O₂ protonation. Experiments designed to test these various hypotheses have revealed some surprises in both P450cam and other bacterial P450s. Molecular dynamics and crystallography show that P450cam can undergo rather significant conformational gymnastics that result in a large restructuring of the active site requiring multiple cis/trans proline isomerizations. It also has been found that X-ray driven substrate hydroxylation is a useful tool for better understanding the role that the essential Asp and surrounding residues play in catalysis. Here we summarize these recent results which provide a much more dynamic picture of P450 catalysis.

The Birth-and-Death Evolution of Cytochrome P450 Genes in Bees

The birth-and-death model of multigene family evolution describes how gene families evolve and diversify through duplication and deletion. The cytochrome P450s are one of the most diverse and well-studied multigene families, involved in both physiological and xenobiotic functions. Extensive studies of insect P450 genes have demonstrated their role in insecticide resistance. Bees are thought to experience toxin exposure through their diet of nectar and pollen, as well as the resin-collecting behavior exhibited by some species. Here, we describe the repertoire of P450 genes in the orchid bee Euglossa dilemma. Male orchid bees form perfume bouquets used in courtship displays by collecting volatile compounds, resulting in exposure to compounds known to be toxic. In addition, we conducted phylogenetic and selection analyses across ten bee species encompassing three bee families. We find that social behavior and resin collection are not correlated with the repertoire of P450 present in a bee species. However, our analyses revealed that P450 clades can be classified as stable and unstable, and that genes involved in xenobiotic metabolism are more likely to belong to unstable clades. Furthermore, we find that unstable clades are under more dynamic evolutionary pressures and exhibit signals of adaptive evolution. This work highlights the complexity of multigene family evolution, revealing that multiple factors contribute to the diversification, stability, and dynamics of this gene family. Furthermore, we provide a resource for future detailed studies investigating the function of different P450s in economically important bee species.

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.

A Modeling-Based Design to Engineering Protein Hydrogels with Random Copolymers

Protein enzymes have shown great potential in numerous technological applications. However, the design of supporting materials is needed to preserve protein functionality outside their native environment. Direct enzyme-polymer self-assembly offers a promising alternative to immobilize proteins in an aqueous solution, achieving higher control of their stability and enzymatic activity in industrial applications. Herein, we propose a modeling-based design to engineering hydrogels of cytochrome P450 and of PETase with styrene/2-vinylpyridine (2VP) random copolymers. By tuning the copolymer fraction of polar groups and of charged groups via quaternization of 2VP for coassembly with cytochrome P450 and via sulfonation of styrene for coassembly with PETase, we provide quantitative guidelines to select either a protein-polymer hydrogel structure or a single-protein encapsulation. The results highlight that, regardless of the protein surface domains, the presence of polar interactions and hydration effects promote the formation of a more elongated enzyme-polymer complex, suggesting a membrane-like coassembly. On the other hand, the effectiveness of a single-protein encapsulation is reached by decreasing the fraction of polar groups and by increasing the charge fraction up to 15%. Our computational analysis demonstrates that the enzyme-polymer assemblies are first promoted by the hydrophobic interactions which lead the protein nonpolar residues to achieve the maximum coverage and to play the role of the most robust contact points. The mechanisms of coassembly are unveiled in the light of both protein and polymer physical-chemistry, providing bioconjugate phase diagrams for the optimal material design.

Nitrene transfers mediated by natural and artificial iron enzymes

Amines are ubiquitous in biology and pharmacy. As a consequence, introducing N functionalities in organic molecules is attracting strong continuous interest. The past decade has witnessed the emergence of very efficient and selective catalytic systems achieving this goal thanks to engineered hemoproteins. In this review, we examine how these enzymes have been engineered focusing rather on the rationale behind it than the methodology employed. These studies are put in perspective with respect to in vitro and in vivo nitrene transfer processes performed by cytochromes P450. An emphasis is put on mechanistic aspects which are confronted to current molecular knowledge of these reactions. Forthcoming developments are delineated.

Enzyme modification using mutation site prediction method for enhancing the regioselectivity of substrate reaction sites

Enzymes with low regioselectivity of substrate reaction sites may produce multiple products from a single substrate. When a target product is produced industrially using these enzymes, the production of non-target products (byproducts) causes adverse effects such as increased processing costs for purification and the amount of raw material. Thus it is required the development of modified enzymes to reduce the amount of byproducts’ production. In this paper, we report a method called mutation site prediction for enhancing the regioselectivity of substrate reaction sites (MSPER). MSPER takes conformational data for docking poses of an enzyme and a substrate as input and automatically generates a ranked list of mutation sites to destabilize docking poses for byproducts while maintaining those for target products in silico. We applied MSPER to the enzyme cytochrome P450 CYP102A1 (BM3) and the two substrates to enhance the regioselectivity for four target products with different reaction sites. The 13 of the total 14 top-ranked mutation sites predicted by MSPER for the four target products succeeded in selectively enhancing the regioselectivity up to 6.4-fold. The results indicate that MSPER can distinguish differences of substrate structures and the reaction sites, and can accurately predict mutation sites to enhance regioselectivity without selection by directed evolution screening.

Structure-based virtual screening of CYP1A1 inhibitors: towards rapid tier-one assessment of potential developmental toxicants

Cytochrome P450 1A1 (CYP1A1) metabolizes estrogens, melatonin, and other key endogenous signaling molecules critical for embryonic/fetal development. The enzyme has increasing expression during pregnancy, and its inhibition or knockout increases embryonic/fetal lethality and/or developmental problems. Here, we present a virtual screening model for CYP1A1 inhibitors based on the orthosteric and predicted allosteric sites of the enzyme. Using 1001 reference compounds with CYP1A1 activity data, we optimized the decision thresholds of our model and classified the training compounds with 68.3% balanced accuracy (91.0% sensitivity and 45.7% specificity). We applied our final model to 11 known CYP1A1 orthosteric binders and related compounds, and found that our ranking of the known orthosteric binders generally agrees with the relative activity of CYP1A1 in metabolizing these compounds. We also applied the model to 22 new test compounds with unknown/unclear CYP1A1 inhibitory activity, and predicted 16 of them are CYP1A1 inhibitors. The CYP1A1 potency and modes of inhibition of these 22 compounds were experimentally determined. We confirmed that most predicted inhibitors, including drugs contraindicated during pregnancy (amiodarone, bicalutamide, cyproterone acetate, ketoconazole, and tamoxifen) and environmental agents suspected to be endocrine disruptors (bisphenol A, diethyl and dibutyl phthalates, and zearalenone), are indeed potent inhibitors of CYP1A1. Our results suggest that virtual screening may be used as a rapid tier-one method to screen for potential CYP1A1 inhibitors, and flag them out for further experimental evaluations.

Insights on intermolecular FMN-heme domain interaction and the role of linker length in cytochrome P450cin fusion proteins

Cytochrome P450s are an important group of enzymes catalyzing hydroxylation, and epoxidations reactions. In this work we describe the characterization of the CinA-CinC fusion enzyme system of a previously reported P450 using genetically fused heme (CinA) and FMN (CinC) enzyme domains from Citrobacter braaki. We observed that mixing individually inactivated heme (-) with FMN (-) domain in the CinA-10aa linker - CinC fusion constructs results in recovered activity and the formation of (2S)-2β-hydroxy,1,8-cineole (174 µM), a similar amount when compared to the fully functional fusion protein (176 µM). We also studied the effect of the fusion linker length in the activity complementation assay. Our results suggests an intermolecular interaction between heme and FMN parts from different CinA-CinC fusion protein similar to proposed mechanisms for P450 BM3 on the other hand, linker length plays a crucial influence on the activity of the fusion constructs. However, complementation assays show that inactive constructs with shorter linker lengths have functional subunits, and that the lack of activity might be due to incorrect interaction between fused enzymes.

Uncovering of cytochrome P450 anatomy by SecStrAnnotator

Protein structural families are groups of homologous proteins defined by the organization of secondary structure elements (SSEs). Nowadays, many families contain vast numbers of structures, and the SSEs can help to orient within them. Communities around specific protein families have even developed specialized SSE annotations, always assigning the same name to the equivalent SSEs in homologous proteins. A detailed analysis of the groups of equivalent SSEs provides an overview of the studied family and enriches the analysis of any particular protein at hand. We developed a workflow for the analysis of the secondary structure anatomy of a protein family. We applied this analysis to the model family of cytochromes P450 (CYPs)-a family of important biotransformation enzymes with a community-wide used SSE annotation. We report the occurrence, typical length and amino acid sequence for the equivalent SSE groups, the conservation/variability of these properties and relationship to the substrate recognition sites. We also suggest a generic residue numbering scheme for the CYP family. Comparing the bacterial and eukaryotic part of the family highlights the significant differences and reveals a well-known anomalous group of bacterial CYPs with some typically eukaryotic features. Our workflow for SSE annotation for CYP and other families can be freely used at address https://sestra.ncbr.muni.cz .

Structural Basis for the Diminished Ligand Binding and Catalytic Ability of Human Fetal-Specific CYP3A7

Cytochrome P450 3A7 (CYP3A7) is a fetal/neonatal liver enzyme that participates in estriol synthesis, clearance of all-trans retinoic acid, and xenobiotic metabolism. Compared to the closely related major drug-metabolizing enzyme in adult liver, CYP3A4, the ligand binding and catalytic capacity of CYP3A7 are substantially reduced. To better understand the structural basis for these functional differences, the 2.15 Å crystal structure of CYP3A7 has been solved. Comparative analysis of CYP3A enzymes shows that decreased structural plasticity rather than the active site microenvironment defines the ligand binding ability of CYP3A7. In particular, a rotameric switch in the gatekeeping amino acid F304 triggers local and long-range rearrangements that transmit to the F-G fragment and alter its interactions with the I-E-D-helical core, resulting in a more rigid structure. Elongation of the β₃-β₄ strands, H-bond linkage in the substrate channel, and steric constraints in the C-terminal loop further increase the active site rigidity and limit conformational ensemble. Collectively, these structural distinctions lower protein plasticity and change the heme environment, which, in turn, could impede the spin-state transition essential for optimal reactivity and oxidation of substrates.

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.

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.

Self-Sufficient Class VII Cytochromes P450: From Full-Length Structure to Synthetic Biology Applications

Members of class VII cytochromes P450 are catalytically self-sufficient enzymes containing a phthalate dioxygenase reductase-like domain fused to the P450 catalytic domain. Among these, CYP116B46 is the first enzyme for which the 3D structure of the whole polypeptide chain has been solved, shedding light on the interaction between its domains, which is crucial for catalysis. Most of these enzymes have been isolated from extremophiles or detoxifying bacteria that can carry out regio- and enantioselective oxidation of compounds of biotechnological interest. Protein engineering has generated mutants that can perform challenging organic reactions such as the anti-Markovnikov alkene oxidation. This potential, combined with the detailed 3D structure, forms the basis for further directed evolution studies aimed at widening their biotechnological exploitation.

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.

Calcitriol increases frataxin levels and restores mitochondrial function in cell models of Friedreich Ataxia

Friedreich ataxia (FA) is a neurodegenerative disease caused by the deficiency of frataxin, a mitochondrial protein. In primary cultures of dorsal root ganglia neurons, we showed that frataxin depletion resulted in decreased levels of the mitochondrial calcium exchanger NCLX, neurite degeneration and apoptotic cell death. Here, we describe that frataxin-deficient dorsal root ganglia neurons display low levels of ferredoxin 1 (FDX1), a mitochondrial Fe/S cluster-containing protein that interacts with frataxin and, interestingly, is essential for the synthesis of calcitriol, the active form of vitamin D. We provide data that calcitriol supplementation, used at nanomolar concentrations, is able to reverse the molecular and cellular markers altered in DRG neurons. Calcitriol is able to recover both FDX1 and NCLX levels and restores mitochondrial membrane potential indicating an overall mitochondrial function improvement. Accordingly, reduction in apoptotic markers and neurite degeneration was observed and, as a result, cell survival was also recovered. All these beneficial effects would be explained by the finding that calcitriol is able to increase the mature frataxin levels in both, frataxin-deficient DRG neurons and cardiomyocytes; remarkably, this increase also occurs in lymphoblastoid cell lines derived from FA patients. In conclusion, these results provide molecular bases to consider calcitriol for an easy and affordable therapeutic approach for FA patients.

AFM and FTIR Investigation of the Effect of Water Flow on Horseradish Peroxidase

Atomic force microscopy (AFM)-based fishing is a promising method for the detection of low-abundant proteins. This method is based on the capturing of the target proteins from the analyzed solution onto a solid substrate, with subsequent counting of the captured protein molecules on the substrate surface by AFM. Protein adsorption onto the substrate surface represents one of the key factors determining the capturing efficiency. Accordingly, studying the factors influencing the protein adsorbability onto the substrate surface represents an actual direction in biomedical research. Herein, the influence of water motion in a flow-based system on the protein adsorbability and on its enzymatic activity has been studied with an example of horseradish peroxidase (HRP) enzyme by AFM, attenuated total reflection Fourier-transform infrared spectroscopy (ATR-FTIR) and conventional spectrophotometry. In the experiments, HRP solution was incubated in a setup modeling the flow section of a biosensor communication. The measuring cell with the protein solution was placed near a coiled silicone pipe, through which water was pumped. The adsorbability of the protein onto the surface of the mica substrate has been studied by AFM. It has been demonstrated that incubation of the HRP solution near the coiled silicone pipe with flowing water leads to an increase in its adsorbability onto mica. This is accompanied by a change in the enzyme's secondary structure, as has been revealed by ATR-FTIR. At the same time, its enzymatic activity remains unchanged. The results reported herein can be useful in the development of models describing the influence of liquid flow on the properties of enzymes and other proteins. The latter is particularly important for the development of biosensors for biomedical applications-particularly for serological analysis, which is intended for the early diagnosis of various types of cancer and infectious diseases. Our results should also be taken into account in studies of the effects of protein aggregation on hemodynamics, which plays a key role in human body functioning.

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

Cytochrome P450 (CYP) heme monooxygenases require two electrons for their catalytic cycle. For mammalian microsomal CYPs, key enzymes for xenobiotic metabolism and steroidogenesis and important drug targets and biocatalysts, the electrons are transferred by NADPH-cytochrome P450 oxidoreductase (CPR). No structure of a mammalian CYP-CPR complex has been solved experimentally, hindering understanding of the determinants of electron transfer (ET), which is often rate-limiting for CYP reactions. Here, we investigated the interactions between membrane-bound CYP 1A1, an antitumor drug target, and CPR by a multiresolution computational approach. We find that upon binding to CPR, the CYP 1A1 catalytic domain becomes less embedded in the membrane and reorients, indicating that CPR may affect ligand passage to the CYP active site. Despite the constraints imposed by membrane binding, we identify several arrangements of CPR around CYP 1A1 that are compatible with ET. In the complexes, the interactions of the CPR FMN domain with the proximal side of CYP 1A1 are supplemented by more transient interactions of the CPR NADP domain with the distal side of CYP 1A1. Computed ET rates and pathways agree well with available experimental data and suggest why the CYP-CPR ET rates are low compared to those of soluble bacterial CYPs.

Advanced Understanding of the Electron Transfer Pathway of Cytochrome P450s

Cytochrome P450s are heme-thiolate enzymes that participate in carbon source assimilation, natural compound biosynthesis and xenobiotic metabolism in all kingdoms of life. P450s can catalyze various reactions by using a wide range of organic compounds, thus exhibiting great potential in biotechnological applications. The catalytic reactions of P450s are driven by electron equivalents that are sourced from pyridine nucleotides and delivered by cognate or matching redox partners (RPs). The electron transfer (ET) route from RPs to P450s involves one or more redox center-containing domains. As the rate of ET is one of the main determinants of P450 efficacy, an in-depth understanding of the P450 ET pathway should increase our knowledge of these important enzymes and benefit their further applications. Here, the various P450 RP systems along with current understanding of their ET routes will be reviewed. Notably, state-of-the-art structural studies of the two main types of self-sufficient P450 will also be summarized.

Structure-based mechanisms: On the way to apply alcohol dehydrogenases/reductases to organic-aqueous systems

Alcohol dehydrogenases/reductases catalyze enantioselective syntheses of versatile chiral compounds relying on direct hydride transfer from cofactor to substrates, or to an intermediate and then to substrates. Since most of the substrates catalyzed by alcohol dehydrogenases/reductases are insoluble in aqueous solutions, increasing interest has been turning to organic-aqueous systems. However, alcohol dehydrogenases/reductases are normally instable in organic solvents, leading to the unsatisfied enantioselective synthesis efficiency. The behaviors of these enzymes in organic solvents at an atomic level are unclear, thus it is of great importance to understand its structure-based mechanisms in organic-aqueous systems to improve their relative stability. Here, we summarized the accessible structures of alcohol dehydrogenases/reductases in Protein Data Bank crystallized in organic-aqueous systems, and compared the structures of alcohol dehydrogenases/reductases which have different tolerance towards organic solvents. By understanding the catalytic behaviors and mechanisms of these enzymes in organic-aqueous systems, the efficient enantioselective syntheses mediated by alcohol dehydrogenases/reductases and further challenges are also discussed through solvent engineering and enzyme-immobilization in the last decade.

Ru(II)-diimine complexes and cytochrome P450 working hand-in-hand

With a growing interest in utilizing visible light to drive biocatalytic processes, several light-harvesting units and approaches have been employed to harness the synthetic potential of heme monooxygenases and carry out selective oxyfunctionalization of a wide range of substrates. While the fields of cytochrome P450 and Ru(II) photochemistry have separately been prolific, it is not until the turn of the 21st century that they converged. Non-covalent and subsequently covalently attached Ru(II) complexes were used to promote rapid intramolecular electron transfer in bacterial P450 enzymes. Photocatalytic activity with Ru(II)-modified P450 enzymes was achieved under reductive conditions with a judicious choice of a sacrificial electron donor. The initial concept of Ru(II)-modified P450 enzymes was further improved using protein engineering, photosensitizer functionalization and was successfully applied to other P450 enzymes. In this review, we wish to present the recent contributions from our group and others in utilizing Ru(II) complexes coupled with P450 enzymes in the broad context of photobiocatalysis, protein assemblies and chemoenzymatic reactions. The merging of chemical catalysts with the synthetic potential of P450 enzymes has led to the development of several chemoenzymatic approaches. Moreover, strained Ru(II) compounds have been shown to selectively inhibit P450 enzymes by releasing aromatic heterocycle containing molecules upon visible light excitation taking advantage of the rapid ligand loss feature in those complexes.

Protein Surface Printer for Exploring Protein Domains

The surface of proteins is vital in determining protein functions. Herein, a program, Protein Surface Printer (PSP), is built that performs multiple functions in quantifying protein surface domains. Two proteins, PETase and cytochrome P450, are used to validate that the program supports atomistic simulations with different combinations of programs and force fields. A case study is conducted on the structural analysis of the spike proteins of SARS-CoV-2 and SARS-CoV and the human cell receptor ACE2. Although the surface domains of both spike proteins are highly similar, their receptor-binding domains (RBDs) and the O-linked glycan domains are structurally different. The O-linked glycan domain of SARS-CoV-2 is highly positively charged, which may promote binding to negatively charged human cells.

Nature's Machinery, Repurposed: Expanding the Repertoire of Iron-Dependent Oxygenases

Iron is an especially important redox-active cofactor in biology because of its ability to mediate reactions with atmospheric O2. Iron-dependent oxygenases exploit this earth-abundant transition metal for the insertion of oxygen atoms into organic compounds. Throughout the astounding diversity of transformations catalyzed by these enzymes, the protein framework directs reactive intermediates toward the precise formation of products, which, in many cases, necessitates the cleavage of strong C-H bonds. In recent years, members of several iron-dependent oxygenase families have been engineered for new-to-nature transformations that offer advantages over conventional synthetic methods. In this Perspective, we first explore what is known about the reactivity of heme-dependent cytochrome P450 oxygenases and nonheme iron-dependent oxygenases bearing the 2-His-1-carboxylate facial triad by reviewing mechanistic studies with an emphasis on how the protein scaffold maximizes the catalytic potential of the iron-heme and iron cofactors. We then review how these cofactors have been repurposed for abiological transformations by engineering the protein frameworks of these enzymes. Finally, we discuss contemporary challenges associated with engineering these platforms and comment on their roles in biocatalysis moving forward.