Crystal structure of cindoxin, the P450cin redox partner

Roles of Ferredoxin-NADP+ Oxidoreductase and Flavodoxin in NAD(P)H-Dependent Electron Transfer Systems

Distinct isoforms of FAD-containing ferredoxin-NADP⁺ oxidoreductase (FNR) and ferredoxin (Fd) are involved in photosynthetic and non-photosynthetic electron transfer systems. The FNR (FAD)-Fd [2Fe-2S] redox pair complex switches between one- and two-electron transfer reactions in steps involving FAD semiquinone intermediates. In cyanobacteria and some algae, one-electron carrier Fd serves as a substitute for low-potential FMN-containing flavodoxin (Fld) during growth under low-iron conditions. This complex evolves into the covalent FNR (FAD)-Fld (FMN) pair, which participates in a wide variety of NAD(P)H-dependent metabolic pathways as an electron donor, including bacterial sulfite reductase, cytochrome P450 BM3, plant or mammalian cytochrome P450 reductase and nitric oxide synthase isoforms. These electron transfer systems share the conserved Ser-Glu/Asp pair in the active site of the FAD module. In addition to physiological electron acceptors, the NAD(P)H-dependent diflavin reductase family catalyzes a one-electron reduction of artificial electron acceptors such as quinone-containing anticancer drugs. Conversely, NAD(P)H: quinone oxidoreductase (NQO1), which shares a Fld-like active site, functions as a typical two-electron transfer antioxidant enzyme, and the NQO1 and UDP-glucuronosyltransfease/sulfotransferase pairs function as an antioxidant detoxification system. In this review, the roles of the plant FNR-Fd and FNR-Fld complex pairs were compared to those of the diflavin reductase (FAD-FMN) family. In the final section, evolutionary aspects of NAD(P)H-dependent multi-domain electron transfer systems are discussed.

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.

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.

Efficient production of oxidized terpenoids via engineering fusion proteins of terpene synthase and cytochrome P450

The functionalization of terpenes using cytochrome P450 enzymes is a versatile route to the production of useful derivatives that can be further converted to value-added products. Many terpenes are hydrophobic and volatile making their availability as a substrate for P450 enzymes significantly limited during microbial production. In this study, we developed a strategy to improve the accessibility of terpene molecules for the P450 reaction by linking terpene synthase and P450 together. As a model system, fusion proteins of 1,8-cineole synthase (CS) and P450_cin were investigated and it showed an improved hydroxylation of the monoterpenoid 1,8-cineole up to 5.4-fold. Structural analysis of the CS-P450_cin fusion proteins by SEC-SAXS indicated a dimer formation with preferred orientations of the active sites of the two domains. We also applied the enzyme fusion strategy to the oxidation of a sesquiterpene epi-isozizaene and the fusion enzymes significantly improved albaflavenol production in engineered E. coli. From the analysis of positive and negative examples of the fusion strategy, we proposed key factors in structure-based prediction and evaluation of fusion enzymes. Developing fusion enzymes for terpene synthase and P450 presents an efficient strategy toward oxidation of hydrophobic terpene compounds. This strategy could be widely applicable to improve the biosynthetic titer of the functionalized products from hydrophobic terpene intermediates.

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.

Evolutionary Relationships Between Low Potential Ferredoxin and Flavodoxin Electron Carriers

Proteins from the ferredoxin (Fd) and flavodoxin (Fld) families function as low potential electrical transfer hubs in cells, at times mediating electron transfer between overlapping sets of oxidoreductases. To better understand protein electron carrier (PEC) use across the domains of life, we evaluated the distribution of genes encoding [4Fe-4S] Fd, [2Fe-2S] Fd, and Fld electron carriers in over 7,000 organisms. Our analysis targeted genes encoding small PEC genes encoding proteins having ≤200 residues. We find that the average number of small PEC genes per Archaea (~13), Bacteria (~8), and Eukarya (~3) genome varies, with some organisms containing as many as 54 total PEC genes. Organisms fall into three groups, including those lacking genes encoding low potential PECs (3%), specialists with a single PEC gene type (20%), and generalists that utilize multiple PEC types (77%). Mapping PEC gene usage onto an evolutionary tree highlights the prevalence of [4Fe-4S] Fds in ancient organisms that are deeply rooted, the expansion of [2Fe-2S] Fds with the advent of photosynthesis and a concomitant decrease in [4Fe-4S] Fds, and the expansion of Flds in organisms that inhabit low-iron host environments. Surprisingly, [4Fe-4S] Fds present a similar abundance in aerobes as [2Fe-2S] Fds. This bioinformatic study highlights understudied PECs whose structure, stability, and partner specificity should be further characterized.

Molecular mechanism of metabolic NAD(P)H-dependent electron-transfer systems: The role of redox cofactors

NAD(P)H-dependent electron-transfer (ET) systems require three functional components: a flavin-containing NAD(P)H-dehydrogenase, one-electron carrier and metal-containing redox center. In principle, these ET systems consist of one-, two- and three-components, and the electron flux from pyridine nucleotide cofactors, NADPH or NADH to final electron acceptor follows a linear pathway: NAD(P)H → flavin → one-electron carrier → metal containing redox center. In each step ET is primarily controlled by one- and two-electron midpoint reduction potentials of protein-bound redox cofactors in which the redox-linked conformational changes during the catalytic cycle are required for the domain-domain interactions. These interactions play an effective ET reactions in the multi-component ET systems. The microsomal and mitochondrial cytochrome P450 (cyt P450) ET systems, nitric oxide synthase (NOS) isozymes, cytochrome b₅ (cyt b₅) ET systems and methionine synthase (MS) ET system include a combination of multi-domain, and their organizations display similarities as well as differences in their components. However, these ET systems are sharing of a similar mechanism. More recent structural information obtained by X-ray and cryo-electron microscopy (cryo-EM) analysis provides more detail for the mechanisms associated with multi-domain ET systems. Therefore, this review summarizes the roles of redox cofactors in the metabolic ET systems on the basis of one-electron redox potentials. In final Section, evolutionary aspects of NAD(P)H-dependent multi-domain ET systems will be discussed.

Structural and catalytic properties of the peroxygenase P450 enzyme CYP152K6 from Bacillus methanolicus

The CYP152 family of cytochrome P450 enzymes (P450s or CYPs) are bacterial peroxygenases that use hydrogen peroxide to drive hydroxylation and decarboxylation of fatty acid substrates. We have expressed and purified a novel CYP152 family member - CYP152K6 from the methylotroph Bacillus methanolicus MGA3. CYP152K6 was characterized using spectroscopic, analytical and structural methods. CYP152K6, like its peroxygenase counterpart P450SPα (CYP152B1) from Sphingomonas paucimobilis, does not undergo significant fatty acid-induced perturbation to the heme spectrum, with the exception of a minor Soret shift observed on binding dodecanoic acid. However, CYP152K6 purified from an E. coli expression system was crystallized and its structure was determined to 1.3 Å with tetradecanoic acid bound. No lipids were present in conditions used for crystallogenesis, and thus CYP152K6 must form a complex by incorporating the fatty acid from E. coli cells. Turnover studies with dodecanoic acid revealed several products, with 2-hydroxydodecanoic acid as the major product and much smaller quantities of 3-hydroxydodecanoic acid. Secondary turnover products were undec-1-en-1-ol, 2-hydroxydodec-2-enoic acid and 2,3-dihydroxydodecanoic acid. This is the first report of a 2,3-hydroxylated fatty acid product made by a peroxygenase P450, with the dihydroxylated product formed by CYP152K6-catalyzed 3-hydroxylation of 2-hydroxydodecanoic acid, but not by 2-hydroxylation of 3-hydroxydodecanoic acid.

Structural and Functional Characterization of a Short-Chain Flavodoxin Associated with a Noncanonical 1,2-Propanediol Utilization Bacterial Microcompartment

Bacterial microcompartments (BMCs) are proteinaceous organelles that encapsulate enzymes involved in CO₂ fixation (carboxysomes) or carbon catabolism (metabolosomes). Metabolosomes share a common core of enzymes and a distinct signature enzyme for substrate degradation that defines the function of the BMC (e.g., propanediol or ethanolamine utilization BMCs, or glycyl-radical enzyme microcompartments). Loci encoding metabolosomes also typically contain genes for proteins that support organelle function, such as regulation, transport of substrate, and cofactor (e.g., vitamin B₁₂) synthesis and recycling. Flavoproteins are frequently among these ancillary gene products, suggesting that these redox active proteins play an undetermined function in many metabolosomes. Here, we report the first characterization of a BMC-associated flavodoxin (Fld1C), a small flavoprotein, derived from the noncanonical 1,2-propanediol utilization BMC locus (PDU1C) of Lactobacillus reuteri. The 2.0 Å X-ray structure of Fld1C displays the α/β flavodoxin fold, which noncovalently binds a single flavin mononucleotide molecule. Fld1C is a short-chain flavodoxin with redox potentials of -240 ± 3 mV oxidized/semiquinone and -344 ± 1 mV semiquinone/hydroquinone versus the standard hydrogen electrode at pH 7.5. It can participate in an electron transfer reaction with a photoreductant to form a stable semiquinone species. Collectively, our structural and functional results suggest that PDU1C BMCs encapsulate Fld1C to store and transfer electrons for the reactivation and/or recycling of the B₁₂ cofactor utilized by the signature enzyme.

Identification, characterization and molecular adaptation of class I redox systems for the production of hydroxylated diterpenoids

Background

De novo production of multi-hydroxylated diterpenoids is challenging due to the lack of efficient redox systems.

Results

In this study a new reductase/ferredoxin system from Streptomyces afghaniensis (AfR·Afx) was identified, which allowed the Escherichia coli-based production of the trihydroxylated diterpene cyclooctatin, a potent inhibitor of human lysophospholipase. This production system provides a 43-fold increase in cyclooctatin yield (15 mg/L) compared to the native producer. AfR·Afx is superior in activating the cylcooctatin-specific class I P450s CotB3/CotB4 compared to the conventional Pseudomonas putida derived PdR·Pdx model. To enhance the activity of the PdR·Pdx system, the molecular basis for these activity differences, was examined by molecular engineering.

Conclusion

We demonstrate that redox system engineering can boost and harmonize the catalytic efficiency of class I hydroxylase enzyme cascades. Enhancing CotB3/CotB4 activities also provided for identification of CotB3 substrate promiscuity and sinularcasbane D production, a functionalized diterpenoid originally isolated from the soft coral Sinularia sp.

Electronic supplementary material

The online version of this article (doi:10.1186/s12934-016-0487-6) contains supplementary material, which is available to authorized users.

"Bind and Crawl" Association Mechanism of Leishmania major Peroxidase and Cytochrome c Revealed by Brownian and Molecular Dynamics Simulations

Leishmania major, the parasitic causative agent of leishmaniasis, produces a heme peroxidase (LmP), which catalyzes the peroxidation of mitochondrial cytochrome c (LmCytc) for protection from reactive oxygen species produced by the host. The association of LmP and LmCytc, which is known from kinetics measurements to be very fast (∼10(8) M(-1) s(-1)), does not involve major conformational changes and has been suggested to be dominated by electrostatic interactions. We used Brownian dynamics simulations to investigate the mechanism of formation of the LmP-LmCytc complex. Our simulations confirm the importance of electrostatic interactions involving the negatively charged D211 residue at the LmP active site, and reveal a previously unrecognized role in complex formation for negatively charged residues in helix A of LmP. The crystal structure of the D211N mutant of LmP reported herein is essentially identical to that of wild-type LmP, reinforcing the notion that it is the loss of charge at the active site, and not a change in structure, that reduces the association rate of the D211N variant of LmP. The Brownian dynamics simulations further show that complex formation occurs via a "bind and crawl" mechanism, in which LmCytc first docks to a location on helix A that is far from the active site, forming an initial encounter complex, and then moves along helix A to the active site. An atomistic molecular dynamics simulation confirms the helix A binding site, and steady state activity assays and stopped-flow kinetics measurements confirm the role of helix A charges in the association mechanism.

Cytochrome P450cin (CYP176A1)

Cytochrome P450cin (P450cin) (CYP176A1) is a bacterial P450 enzyme that catalyses the enantiospecific hydroxylation of 1,8-cineole to (1R)-6β-hydroxycineole when reconstituted with its natural reduction-oxidation (redox) partner cindoxin, E. coli flavodoxin reductase, and NADPH as a source of electrons. This catalytic system has become a useful tool in the study of P450s as not only can large quantities of P450cin be prepared and rates of oxidation up to 1,500 min(-1) achieved, but it also displays a number of unusual characteristics. These include an asparagine residue in P450cin that has been found in place of the usual conserved threonine residue observed in most P450s. In general, this conserved threonine controls oxygen activation to create the potent ferryl (Fe(IV=O) porphyrin cation radical required for substrate oxidation. Another atypical characteristic of P450cin is that it utilises an FMN-containing redoxin (cindoxin) rather than a ferridoxin as is usually observed with other bacterial P450s (e.g. P450cam). This chapter will review what is currently known about P450cin and how this enzyme has provided a greater understanding of P450s in general.

webSDA: a web server to simulate macromolecular diffusional association

Macromolecular interactions play a crucial role in biological systems. Simulation of diffusional association (SDA) is a software for carrying out Brownian dynamics simulations that can be used to study the interactions between two or more biological macromolecules. webSDA allows users to run Brownian dynamics simulations with SDA to study bimolecular association and encounter complex formation, to compute association rate constants, and to investigate macromolecular crowding using atomically detailed macromolecular structures. webSDA facilitates and automates the use of the SDA software, and offers user-friendly visualization of results. webSDA currently has three modules: 'SDA docking' to generate structures of the diffusional encounter complexes of two macromolecules, 'SDA association' to calculate bimolecular diffusional association rate constants, and 'SDA multiple molecules' to simulate the diffusive motion of hundreds of macromolecules. webSDA is freely available to all users and there is no login requirement. webSDA is available at http://mcm.h-its.org/webSDA/.