Structural insights into the semiquinone form of human Cytochrome P450 reductase by DEER distance measurements between a native flavin and a spin labelled non-canonical amino acid.

The flavoprotein Cytochrome P450 reductase (CPR) is the unique electron pathway from NADPH to Cytochrome P450 (CYPs). The conformational dynamics of human CPR in solution, which involves transitions from a "locked/closed" to an "unlocked/open" state, is crucial for electron transfer. To date, however, the factors guiding these changes remain unknown. By Site-Directed Spin Labelling coupled to Electron Paramagnetic Resonance spectroscopy, we have incorporated a non-canonical amino acid onto the flavin mononucleotide (FMN) and flavin adenine dinucleotide (FAD) domains of soluble human CPR, and labelled it with a specific nitroxide spin probe. Taking advantage of the endogenous FMN cofactor, we successfully measured for the first time, the distance distribution by DEER between the semiquinone state FMNH• and the nitroxide. The DEER data revealed a salt concentration-dependent distance distribution, evidence of an "open" CPR conformation at high salt concentrations exceeding previous reports. We also conducted molecular dynamics simulations which unveiled a diverse ensemble of conformations for the "open" semiquinone state of the CPR at high salt concentration. This study unravels the conformational landscape of the one electron reduced state of CPR, which had never been studied before.

Inferring assembly-curving trends of bacterial micro-compartment shell hexamers from crystal structure arrangements.

Bacterial microcompartments (BMC) are complex macromolecular assemblies that participate in varied chemical processes in about one fourth of bacterial species. BMC-encapsulated enzymatic activities are segregated from other cell contents by means of semipermeable shells, justifying why BMC are viewed as prototype nano-reactors for biotechnological applications. Herein, we undertook a comparative study of bending propensities of BMC hexamers (BMC-H), the most abundant shell constituents. Published data show that some BMC-H, like ß-carboxysomal CcmK, tend to assemble flat whereas other BMC-H often build curved objects. Inspection of available crystal structures presenting BMC-H in tiled arrangements permitted us to identify two major assembly modes with a striking connection with experimental trends. All-atom molecular dynamics (MD) supported that BMC-H bending is triggered robustly only from the arrangement adopted in crystals by BMC-H that experimentally form curved objects, leading to very similar arrangements to those found in structures of recomposed BMC shells. Simulations on triplets of planar-behaving hexamers, which were previously reconfigured to comply with such organization, confirmed that bending propensity is mostly defined by the precise lateral positioning of hexamers, rather than by BMC-H identity. Finally, an interfacial lysine was pinpointed as the most decisive residue in controlling PduA spontaneous curvature. Globally, results presented herein should contribute to improve our understanding of the variable mechanisms of biogenesis characterized for BMC, and of possible strategies to regulate BMC size and shape.

Interaction Modes of Microsomal Cytochrome P450s with Its Reductase and the Role of Substrate Binding.

The activity of microsomal cytochromes P450 (CYP) is strictly dependent on the supply of electrons provided by NADPH cytochrome P450 oxidoreductase (CPR). The variant nature of the isoform-specific proximal interface of microsomal CYPs implies that the interacting interface between the two proteins is degenerated. Recently, we demonstrated that specific CPR mutations in the FMN-domain (FD) may induce a gain in activity for a specific CYP isoform. In the current report, we confirm the CYP isoform dependence of CPR's degenerated binding by demonstrating that the effect of four of the formerly studied FD mutants are indeed exclusive of a specific CYP isoform, as verified by cytochrome c inhibition studies. Moreover, the nature of CYP's substrate seems to have a modulating role in the CPR:CYP interaction. In silico molecular dynamics simulations of the FD evidence that mutations induces very subtle structural alterations, influencing the characteristics of residues formerly implicated in the CPR:CYP interaction or in positioning of the FMN moiety. CPR seems therefore to be able to form effective interaction complexes with its structural diverse partners via a combination of specific structural features of the FD, which are functional in a CYP isoform dependent manner, and dependent on the substrate bound.

Accurate Determination of Human CPR Conformational Equilibrium by smFRET Using Dual Orthogonal Noncanonical Amino Acid Labeling.

Conjugation of fluorescent dyes to proteins-a prerequisite for the study of conformational dynamics by single-molecule (sm) FRET-can lead to substantial changes in a dye's photophysical properties, ultimately biasing the determination of inter-dye distances. In particular, cyanine dyes and their derivatives, the most commonly used dyes in smFRET experiments, exhibit such behavior. To overcome this, we developed a general strategy to equip proteins site-specifically with FRET pairs through chemoselective reactions with two distinct noncanonical amino acids simultaneously incorporated through genetic code expansion in Escherichia coli. Application of this technique to human NADPH-cytochrome P450 reductase (CPR) demonstrated the importance of homogenously labeled samples for accurate determination of FRET efficiencies and unveiled the effect of NADP+ on the ionic-strength-dependent modulation of the conformational equilibrium of CPR. Thanks to its generality and accuracy, the presented methodology establishes a new benchmark for deciphering of complex molecular dynamics in single molecules.

Functional analysis of isoprenoid precursors biosynthesis by quantitative metabolomics and isotopologue profiling.

INTRODUCTION: Isoprenoids are amongst the most abundant and diverse biological molecules and are involved in a broad range of biological functions. Functional understanding of their biosynthesis is thus key in many fundamental and applicative fields, including systems biology, medicine and biotechnology. However, available methods do not yet allow accurate quantification and tracing of stable isotopes incorporation for all the isoprenoids precursors. OBJECTIVES: We developed and validated a complete methodology for quantitative metabolomics and isotopologue profiling of isoprenoid precursors in the yeast Saccharomyces cerevisiae. METHODS: This workflow covers all the experimental and computational steps from sample collection and preparation to data acquisition and processing. It also includes a novel quantification method based on liquid chromatography coupled to high-resolution mass spectrometry. Method validation followed the Metabolomics Standards Initiative guidelines. RESULTS: This workflow ensures accurate absolute quantification (RSD < 20%) of all mevalonate and prenyl pyrophosphates intermediates with a high sensitivity over a large linear range (from 0.1 to 50 pmol). In addition, we demonstrate that this workflow brings crucial information to design more efficient phytoene producers. Results indicate stable turnover rates of prenyl pyrophosphate intermediates in the constructed strains and provide quantitative information on the change of the biosynthetic flux of phytoene precursors. CONCLUSION: This methodology fills one of the last technical gaps for functional studies of isoprenoids biosynthesis and should be applicable to other eukaryotic and prokaryotic (micro)organisms after adaptation of some organism-dependent steps. This methodology also opens the way to 13C-metabolic flux analysis of isoprenoid biosynthesis.

Ligand Access Channels in Cytochrome P450 Enzymes: A Review.

Quantitative structure-activity relationships may bring invaluable information on structural elements of both enzymes and substrates that, together, govern substrate specificity. Buried active sites in cytochrome P450 enzymes are connected to the solvent by a network of channels exiting at the distal surface of the protein. This review presents different in silico tools that were developed to uncover such channels in P450 crystal structures. It also lists some of the experimental evidence that actually suggest that these predicted channels might indeed play a critical role in modulating P450 functions. Amino acid residues at the entrance of the channels may participate to a first global ligand recognition of ligands by P450 enzymes before they reach the buried active site. Moreover, different P450 enzymes show different networks of predicted channels. The plasticity of P450 structures is also important to take into account when looking at how channels might play their role.

Probing the Role of the Hinge Segment of Cytochrome P450 Oxidoreductase in the Interaction with Cytochrome P450.

NADPH-cytochrome P450 reductase (CPR) is the unique redox partner of microsomal cytochrome P450s (CYPs). CPR exists in a conformational equilibrium between open and closed conformations throughout its electron transfer (ET) function. Previously, we have shown that electrostatic and flexibility properties of the hinge segment of CPR are critical for ET. Three mutants of human CPR were studied (S243P, I245P and R246A) and combined with representative human drug-metabolizing CYPs (isoforms 1A2, 2A6 and 3A4). To probe the effect of these hinge mutations different experimental approaches were employed: CYP bioactivation capacity of pre-carcinogens, enzyme kinetic analysis, and effect of the ionic strength and cytochrome b5 (CYB5) on CYP activity. The hinge mutations influenced the bioactivation of pre-carcinogens, which seemed CYP isoform and substrate dependent. The deviations of Michaelis-Menten kinetic parameters uncovered tend to confirm this discrepancy, which was confirmed by CYP and hinge mutant specific salt/activity profiles. CPR/CYB5 competition experiments indicated a less important role of affinity in CPR/CYP interaction. Overall, our data suggest that the highly flexible hinge of CPR is responsible for the existence of a conformational aggregate of different open CPR conformers enabling ET-interaction with structural varied redox partners.

Ordered chimerogenesis applied to CYP2B P450 enzymes.

BACKGROUND: Structural studies on CYP2B enzymes identified some of the features that are related to their high plasticity. The aim of this work was to understand further the possible relationships between combinations of structural elements and functions by linking shift in substrate specificity with sequence element swaps between CYP2B6 and CYP2B11. METHODS: A series of 15 chimeras in which a small CYP2B6 sequence segment was swapped with its equivalent in CYP2B11 were constructed. All chimeras produced were thus mostly of CYP2B11 sequence. Time course studies were carried out with two typical CYP2B substrates, cyclophosphamide and 7-ethoxy-4-trifluoromethylcoumarin. Steady-state kinetic parameters were determined for all chimeras expressed in yeast. RESULTS: Most of the chimeras exhibit a high affinity for cyclophosphamide, as CYP2B11 does. A few exhibit an affinity similar to that of CYP2B6 without altered behavior toward the other substrate assayed. The swapped elements that control this specificity shift are discussed in terms of F'/G' cassette role and substrate access channels. CONCLUSIONS: Some sequence segments control precisely the shift in affinity for cyclophosphamide between CYP2B6, which has a typical low affinity, and CYP2B11 which has a typical high affinity. GENERAL SIGNIFICANCE: The result provides a new basis for determining the structural elements that control functions in complex enzymes.

Access channels to the buried active site control substrate specificity in CYP1A P450 enzymes.

BACKGROUND: A cytochrome P450 active site is buried within the protein molecule and several channels connect the catalytic cavity to the protein surface. Their role in P450 catalysis is still matter of debate. The aim of this study was to understand the possible relations existing between channels and substrate specificity. METHODS: Time course studies were carried out with a collection of polycyclic substrates of increasing sizes assayed with a library of wild-type and chimeric CYP1A enzymes. This resulted in a matrix of activities sufficiently large to allow statistical analysis. Multivariate statistical tools were used to decipher the correlation between observed activity shifts and sequence segment swaps. RESULTS: The global kinetic behavior of CYP1A enzymes toward polycyclic substrates is significantly different depending on the size of the substrate. Mutations which are close or lining the P450 channels significantly affect this discrimination, whereas mutations distant from the P450 channels do not. CONCLUSIONS: Size discrimination is taking place for polycyclic substrates at the entrance of the different P450 access channels. It is thus hypothesized that channels differentiate small from large substrates in CYP1A enzymes, implying that residues located at the surface of the protein may be implied in this differential recognition. GENERAL SIGNIFICANCE: Catalysis thus occurs after a two-step recognition process, one at the surface of the protein and the second within the catalytic cavity in enzymes with a buried active site.

A well-balanced preexisting equilibrium governs electron flux efficiency of a multidomain diflavin reductase.

Diflavin reductases are bidomain electron transfer proteins in which structural reorientation is necessary to account for the various intramolecular and intermolecular electron transfer steps. Using small-angle x-ray scattering and nuclear magnetic resonance data, we describe the conformational free-energy landscape of the NADPH-cytochrome P450 reductase (CPR), a typical bidomain redox enzyme composed of two covalently-bound flavin domains, under various experimental conditions. The CPR enzyme exists in a salt- and pH-dependent rapid equilibrium between a previously described rigid, locked state and a newly characterized, highly flexible, unlocked state. We further establish that maximal electron flux through CPR is conditioned by adjustable stability of the locked-state domain interface under resting conditions. This is rationalized by a kinetic scheme coupling rapid conformational sampling and slow chemical reaction rates. Regulated domain interface stability associated with fast stochastic domain contacts during the catalytic cycle thus provides, to our knowledge, a new paradigm for improving our understanding of multidomain enzyme function.

Dynamics of α-Hb chain binding to its chaperone AHSP depends on heme coordination and redox state.

BACKGROUND: AHSP is an erythroid molecular chaperone of the α-hemoglobin chains (α-Hb). Upon AHSP binding, native ferric α-Hb undergoes an unprecedented structural rearrangement at the heme site giving rise to a 6th coordination bond with His(E7). METHODS: Recombinant AHSP, WT α-Hb:AHSP and α-Hb(HE7Q):AHSP complexes were expressed in Escherichia coli. Thermal denaturation curves were measured by circular dichroism for the isolated α-Hb and bound to AHSP. Kinetics of ligand binding and redox reactions of α-Hb bound to AHSP as well as α-Hb release from the α-Hb:AHSP complex were measured by time-resolved absorption spectroscopy. RESULTS: AHSP binding to α-Hb is kinetically controlled to prevail over direct binding with ß-chains and is also thermodynamically controlled by the α-Hb redox state and not the liganded state of the ferrous α-Hb. The dramatic instability of isolated ferric α-Hb is greatly decreased upon AHSP binding. Removing the bis-histidyl hexacoordination in α-HbH58(E7)Q:AHSP complex reduces the stabilizing effect of AHSP binding. Once the ferric α-Hb is bound to AHSP, the globin can be more easily reduced by several chemical and enzymatic systems compared to α-Hb within the Hb-tetramer. CONCLUSION: α-Hb reduction could trigger its release from AHSP toward its final Hb ß-chain partner producing functional ferrous Hb-tetramers. This work indicates a preferred kinetic pathway for Hb-synthesis. GENERAL SIGNIFICANCE: The cellular redox balance in Hb-synthesis should be considered as important as the relative proportional synthesis of both Hb-subunits and their heme cofactor. The in vivo role of AHSP is discussed in the context of the molecular disorders observed in thalassemia.

Exploring linker's sequence diversity to fuse carotene cyclase and hydroxylase for zeaxanthin biosynthesis.

Fusion of catalytic domains can accelerate cascade reactions by bringing enzymes in close proximity. However, the design of a protein fusion and the choice of a linker are often challenging and lack of guidance. To determine the impact of linker parameters on fusion proteins, a library of linkers featuring various lengths, secondary structures, extensions and hydrophobicities was designed. Linkers were used to fuse the lycopene cyclase (crtY) and ß-carotene hydroxylase (crtZ) from Pantoea ananatis to create fusion proteins to produce zeaxanthin. The fusion efficiency was assessed by comparing the carotenoids content in a carotenoid-producer Escherichia coli strain. It was shown that in addition to the orientation of the enzymes and the size of the linker, the first amino acid of the linker is also a key factor in determining the efficiency of a protein fusion. The wide range of sequence diversity in our linker library enables the fine tuning of protein fusion and this approach can be easily transferred to other enzyme couples.

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

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

Optogenetic control of beta-carotene bioproduction in yeast across multiple lab-scales.

Optogenetics arises as a valuable tool to precisely control genetic circuits in microbial cell factories. Light control holds the promise of optimizing bioproduction methods and maximizing yields, but its implementation at different steps of the strain development process and at different culture scales remains challenging. In this study, we aim to control beta-carotene bioproduction using optogenetics in Saccharomyces cerevisiae and investigate how its performance translates across culture scales. We built four lab-scale illumination devices, each handling different culture volumes, and each having specific illumination characteristics and cultivating conditions. We evaluated optogenetic activation and beta-carotene production across devices and optimized them both independently. Then, we combined optogenetic induction and beta-carotene production to make a light-inducible beta-carotene producer strain. This was achieved by placing the transcription of the bifunctional lycopene cyclase/phytoene synthase CrtYB under the control of the pC120 optogenetic promoter regulated by the EL222-VP16 light-activated transcription factor, while other carotenogenic enzymes (CrtI, CrtE, tHMG) were expressed constitutively. We show that illumination, culture volume and shaking impact differently optogenetic activation and beta-carotene production across devices. This enabled us to determine the best culture conditions to maximize light-induced beta-carotene production in each of the devices. Our study exemplifies the stakes of scaling up optogenetics in devices of different lab scales and sheds light on the interplays and potential conflicts between optogenetic control and metabolic pathway efficiency. As a general principle, we propose that it is important to first optimize both components of the system independently, before combining them into optogenetic producing strains to avoid extensive troubleshooting. We anticipate that our results can help designing both strains and devices that could eventually lead to larger scale systems in an effort to bring optogenetics to the industrial scale.

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

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

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

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

Role of the interface between the FMN and FAD domains in the control of redox potential and electronic transfer of NADPH-cytochrome P450 reductase.

CPR (NADPH-cytochrome P450 reductase) is a multidomain protein containing two flavin-containing domains joined by a connecting domain thought to control the necessary movements of the catalytic domains during electronic cycles. We present a detailed biochemical analysis of two chimaeric CPRs composed of the association of human or yeast FMN with the alternative connecting/FAD domains. Despite the assembly of domains having a relatively large evolutionary distance between them, our data support the idea that the integrity of the catalytic cycle is conserved in our chimaeric enzymes, whereas the recognition, interactions and positioning of both catalytic domains are probably modified. The main consequences of the chimaerogenesis are a decrease in the internal electron-transfer rate between both flavins correlated with changes in the geometry of chimaeric CPRs in solution. Results of the present study highlight the role of the linker and connecting domain in the recognition at the interfaces between the catalytic domains and the impact of interdomain interactions on the redox potentials of the flavins, the internal electron-transfer efficiency and the global conformation and dynamic equilibrium of the CPRs.

Dynamic control of electron transfers in diflavin reductases.

Diflavin reductases are essential proteins capable of splitting the two-electron flux from reduced pyridine nucleotides to a variety of one electron acceptors. The primary sequence of diflavin reductases shows a conserved domain organization harboring two catalytic domains bound to the FAD and FMN flavins sandwiched by one or several non-catalytic domains. The catalytic domains are analogous to existing globular proteins: the FMN domain is analogous to flavodoxins while the FAD domain resembles ferredoxin reductases. The first structural determination of one member of the diflavin reductases family raised some questions about the architecture of the enzyme during catalysis: both FMN and FAD were in perfect position for interflavin transfers but the steric hindrance of the FAD domain rapidly prompted more complex hypotheses on the possible mechanisms for the electron transfer from FMN to external acceptors. Hypotheses of domain reorganization during catalysis in the context of the different members of this family were given by many groups during the past twenty years. This review will address the recent advances in various structural approaches that have highlighted specific dynamic features of diflavin reductases.

Structure of the open conformation of a functional chimeric NADPH cytochrome P450 reductase.

Two catalytic domains, bearing FMN and FAD cofactors, joined by a connecting domain, compose the core of the NADPH cytochrome P450 reductase (CPR). The FMN domain of CPR mediates electron shuttling from the FAD domain to cytochromes P450. Together, both enzymes form the main mixed-function oxidase system that participates in the metabolism of endo- and xenobiotic compounds in mammals. Available CPR structures show a closed conformation, with the two cofactors in tight proximity, which is consistent with FAD-to-FMN, but not FMN-to-P450, electron transfer. Here, we report the 2.5 A resolution crystal structure of a functionally competent yeast-human chimeric CPR in an open conformation, compatible with FMN-to-P450 electron transfer. Comparison with closed structures shows a major conformational change separating the FMN and FAD cofactors from 86 A.

A tripartite carbohydrate-binding module to functionalize cellulose nanocrystals.

The development of protein and microorganism engineering have led to rising expectations of biotechnology in the design of emerging biomaterials, putatively of high interest to reduce our dependence on fossil carbon resources. In this way, cellulose, a renewable carbon based polysaccharide and derived products, displays unique properties used in many industrial applications. Although the functionalization of cellulose is common, it is however limited in terms of number and type of functions. In this work, a Carbohydrate-Binding Module (CBM) was used as a central core to provide a versatile strategy to bring a large diversity of functions to cellulose surfaces. CBM3a from Clostridium thermocellum, which has a high affinity for crystalline cellulose, was flanked through linkers with a streptavidin domain and an azide group introduced through a non-canonical amino acid. Each of these two extra domains was effectively produced and functionalized with a variety of biological and chemical molecules. Structural properties of the resulting tripartite chimeric protein were investigated using molecular modelling approaches, and its potential for the multi-functionalization of cellulose was confirmed experimentally. As a proof of concept, we show that cellulose can be labelled with a fluorescent version of the tripartite protein grafted to magnetic beads and captured using a magnet.