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

Cytochrome b5 (b5) plays an important role in the enhancement of the activity of specific cytochrome P450s, especially CYP3A4. Contradictory interpretations were proposed that alternatively defend conformational or redox-based mechanisms. A semi-synthetic system involving CYP3A4, its reductase (CPR), and native or modified b5 bound to the ER membrane was used to clarify this question. This took advantage of the use of b5 analogs featuring cofactors with modified redox properties, in combination with rapid kinetic experiments, activity determinations, and structural predictions. The results supported a redox-based mechanism of activation in which the b5 molecules remain partially confined with specific CYP3A4 partners throughout the catalytic cycle, acting both as an electron acceptor and donor. Structural models of corresponding complexes were built using AlphaFold predictions supplemented by holoenzyme reconstruction and molecular dynamics. Results supported the formation of mutually exclusive CPR-P450 and b5-P450 binary complexes, each with electron transfer capacity. The potential formation of an unusual ternary complex allowing a more delocalized electron exchange between P450 heme, FAD, FMN of CPR and b5 heme was also considered. Modeled formation of these alternate complexes also provided potential structural interpretations for non-redox complementary mechanisms that could contribute to the global activation mechanism of CYP3A4.

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.

BACKGROUND:Physarum polycephalum is an unusual macroscopic myxomycete expressing a large range of glycosyl hydrolases. Among them, enzymes from the GH18 family can hydrolyze chitin, an important structural component of the cell walls in fungi and in the exoskeleton of insects and crustaceans.METHODS:Low stringency sequence signature search in transcriptomes was used to identify GH18 sequences related to chitinases. Identified sequences were expressed in E. coli and corresponding structures modelled. Synthetic substrates and in some cases colloidal chitin were used to characterize activities.RESULTS:Catalytically functional hits were sorted and their predicted structures compared. All share the TIM barrel structure of the GH18 chitinase catalytic domain, optionally fused to binding motifs, such as CBM50, CBM18, and CBM14, involved in sugar recognition. Assessment of the enzymatic activities following deletion of the C-terminal CBM14 domain of the most active clone evidenced a significant contribution of this extension to the chitinase activity. A classification based on module organization, functional and structural criteria of characterized enzymes was proposed.CONCLUSIONS:Physarum polycephalum sequences encompassing a chitinase like GH18 signature share a modular structure involving a structurally conserved catalytic TIM barrels decorated or not by a chitin insertion domain and optionally surrounded by additional sugar binding domains. One of them plays a clear role in enhancing activities toward natural chitin.GENERAL SIGNIFICANCE:Myxomycete enzymes are currently poorly characterized and constitute a potential source for new catalysts. Among them glycosyl hydrolases have a strong potential for valorization of industrial waste as well as in therapeutic field.

Large range structural dynamics plays a critical role for the function of electron transfer proteins. This information is generally not available from crystallographic structures, while cryo-electron microscopy (cryo-EM) can provide some elements but frequently with a degraded spatial resolution. Recently, AlphaFold-based structural modelling was extended to the prediction of protein complexes. In this work, bacterial CYP102A1 from Priestia megaterium was used as a test case to evaluate the capability of AlphaFold2 to predict alternative structures critical for catalysis. CYP102A1 monooxygenase, a NADPH-supported fatty acid hydroxylase, works as a soluble homodimer, each monomer harboring two flavins (FAD and FMN) and one heme cofactors. Large conformational changes are required during catalytic cycle to allow successive electron transfers from FAD to FMN and finally heme iron. We used the recently released AlphaFold2_advanced notebook (AF2A), to predict the possible alternate conformations supporting electron transfers in CYP102A1 homodimer. Challenging AF2A-derived models with previously reported experimental data revealed an unforeseen domain connectivity of the diflavin reductase part of the enzyme. Intermolecular crossed complex constitutes a novel type of structural organization never previously described. The predicted formation within the dimer of a stable complex between the heme containing domains was challenged and found consistent with uninterpreted features of reported crystallographic structures and cryo-EM imaging. The particularly efficient CYP102A1 catalytic mechanism was revisited to the light of the new evidenced connectivity in which the FMN-binding domain of each monomer oscillates on themselves to alternatively receive and transfer electrons without needing large structural change in the dimer. Such model was found explanatory for previously contradictory reported biochemical data. Possibility to mimic CYP102A1 structural organization into bicomponent eukaryotic P450 systems was evaluated by designing and modeling in silico synthetic reductase domains built from composite sequence segments from P. megaterium and human origins. More generally, this work illustrates how the ability of AF2A to predict alternate complex structures can enlighten and explain conformational changes critical for bio-assemblies.

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.

Application of droplet-based microfluidics for the screening of microbial libraries is one of the important ongoing developments in functional genomics/metagenomics. In this article, we propose a new method that can be employed for high-throughput profiling of cell growth. It consists of light-driven labelling droplets that contain growing cells directly in a microfluidics observation chamber, followed by recovery of the labelled cells. This method is based on intracellular expression of green-to-red switchable fluorescent proteins. The proof of concept is established here for two commonly used biological models, E. coli and S. cerevisiae. Growth of cells in droplets was monitored under a microscope and, depending on the targeted phenotype, the fluorescence of selected droplets was switched from a "green" to a "red" state. Red fluorescent cells from labelled droplets were then successfully detected, sorted with the Fluorescence Activated Cell Sorting machine and recovered. Finally, the application of this method for different kind of screenings, in particular of metagenomic libraries, is discussed and this idea is validated by the analysis of a model mini-library.

Adjustment of plasmid copy number resulting from the balance between positive and negative impacts of borne synthetic genes, plays a critical role in the global efficiency of multistep metabolic engineering. Differential expression of co-expressed engineered genes is frequently observed depending on growth phases, metabolic status and triggered adjustments of plasmid copy numbers, constituting a dynamic process contributing to minimize global engineering burden. A yeast model involving plasmid based expression of phosphoribulokinase (PRKp), a key enzyme for the reconstruction of synthetic Calvin cycle, was designed to gain further insights into such a mechanism. A conditional PRK expression cassette was cloned either onto a low (ARS-CEN based) or a high (2-micron origin based) copy number plasmid using complementation of a trp1 genomic mutation as constant positive selection. Evolution of plasmid copy numbers, PRKp expressions, and cell growth rates were dynamically monitored following gene de-repression through external doxycycline concentration shifts. In the absence of RubisCO encoding gene permitting metabolic recycling, PRKp expression that led to depletion of ribulose phosphate, a critical metabolite for aromatic amino-acids biosynthesis, and accumulation of the dead-end diphosphate product contribute to toxicity. Triggered copy number adjustment was found to be a dynamic process depending both on plasmid types and levels of PRK induction. With the ARS-CEN plasmid, cell growth was abruptly affected only when level PRKp expression exceeded a threshold value. In contrast, a proportional relationship was observed with the 2-micron plasmid consistent with large copy number adjustments. Micro-compartment partitioning of bulk cultures by embedding individual cells into inverse culture medium/oil droplets, revealed the presence of slow and fast growing subpopulations that differ in relative proportions for low and high copy number plasmids.

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.

CcmK proteins are major constituents of icosahedral shells of β-carboxysomes, a bacterial microcompartment that plays a key role for CO2 fixation in nature. Supported by the characterization of bidimensional (2D) layers of packed CcmK hexamers in crystal and electron microscopy structures, CcmK are assumed to be the major components of icosahedral flat facets. Here, we reassessed the validity of this model by studying CcmK isoforms from Synechocystis sp. PCC6803. Native mass spectrometry studies confirmed that CcmK are hexamers in solution. Interestingly, potential pre-assembled intermediates were also detected with CcmK2. Atomic-force microscopy (AFM) imaging under quasi-physiological conditions confirmed the formation of canonical flat sheets with CcmK4. Conversely, CcmK2 formed both canonical and striped-patterned patches, while CcmK1 assembled into remarkable supra-hexameric curved honeycomb-like mosaics. Mutational studies ascribed the propensity of CcmK1 to form round assemblies to a combination of two features shared by at least one CcmK isoform in most β-cyanobacteria: a displacement of an α helical portion towards the hexamer edge, where a potential phosphate binding funnel forms between packed hexamers, and the presence of a short C-terminal extension in CcmK1. All-atom molecular dynamics supported a contribution of phosphate molecules sandwiched between hexamers to bend CcmK1 assemblies. Formation of supra-hexameric curved structures could be reproduced in coarse-grained simulations, provided that adhesion forces to the support were weak. Apart from uncovering unprecedented CcmK self-assembly features, our data suggest the possibility that transitions between curved and flat assemblies, following cargo maturation, could be important for the biogenesis of β-carboxysomes, possibly also of other BMC.

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-trifiuoromethylcoumarin. 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. (C) 2016 Elsevier B.V. All rights reserved.

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.

Fine tuning of individual enzyme expression level is necessary to alleviate metabolic imbalances in synthetic heterologous pathways. A known approach consists of choosing a suitable combination of promoters, based on their characterized strengths in model conditions. We questioned whether each step of a multiple-gene synthetic pathway could be independently tunable at the transcription level. Three open reading frames, coding for enzymes involved in a synthetic pathway, were combinatorially associated to different promoters on an episomal plasmid in Saccharomyces cerevisiae. We quantified the mRNA levels of the three genes in each strain of our generated combinatorial metabolic library. Our results evidenced that the ORF nature, position, and orientation induce strong discrepancies between the previously reported promoters' strengths and the observed ones. We conclude that, in the context of metabolic reconstruction, the strength of usual promoters can be dramatically affected by many factors. Among them, transcriptional interference and ORF nature seem to be predominant.

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. (C) 2014 Elsevier B.V. All rights reserved.

Background: Zebrafish express five cytochrome P450 1 genes: CYP1A, CYP1B1, CYP1C1, CYP1C2, inducible by aryl hydrocarbon receptor agonists, and CYP1D1, a constitutively expressed CYP1A-like gene. We examined substrate selectivity of CYP1s expressed in yeast.Methods: CYP1s were expressed in W(R) yeast, engineered to over-express P450 reductase, via pYES/DEST52 and via pYeDP60. Microsomal fractions from transformed yeast were examined for activity with fluorogenic substrates, benzo[a]pyrene and testosterone. Modeling and docking approaches were used to further evaluate sites of oxidation on benzo[a]pyrene and testosterone.Results: CYP1s expressed in yeast dealkylated ethoxy-, methoxy-, pentoxy- and benzoxy-resorufin (EROD, MROD, PROD, BROD). CYP1A and CYP1C2 had the highest rates of EROD activity, while PROD and BROD activities were low for all five CYP1s. The relative rates of resorufin dealkylation by CYP1C1, CYP1C2 and CYP1D1 expressed via pYeDP60 were highly similar to relative rates obtained with pYES/DEST52-expressed enzymes. CYP1C1 and CYP1C2 dealkylated substituted coumarins and ethoxy-fluorescein-ethylester, while CYP1D1 did not. The CYP1Cs and CYP1D1 co-expressed with epoxide hydrolase oxidized BaP with different rates and product profiles, and all three produced BaP-7,8,9,10-tetrol. The CYP1Cs but not CYP1D1 metabolized testosterone to 6 beta-OH-Testosterone. However, CYP1D1 formed an unidentified testosterone metabolite better than the CYP1Cs. Testosterone and BaP docked to CYP homology models with poses consistent with differing product profiles.Conclusions: Yeast-expressed zebrafish CYP1s will be useful in determining further functionality with endogenous and xenobiotic compounds. General significance: Determining the roles of zebrafish CYP1s in physiology and toxicology depends on knowing the substrate selectivity of these enzymes. (C) 2015 Elsevier B.V. All rights reserved.

The promiscuity of a collection of enzymes consisting of 31 wild-type and synthetic variants of CYP1A enzymes was evaluated using a series of 14 steroids and 2 steroid-like chemicals, namely, nootkatone, a terpenoid, and mifepristone, a drug. For each enzyme-substrate couple, the initial steady-state velocity of metabolite formation was determined at a substrate saturating concentration. For that, a high-throughput approach was designed involving automatized incubations in 96-well microplate with sixteen 6-point kinetics per microplate and data acquisition using LC/MS system accepting 96-well microplate for injections. The resulting dataset was used for multivariate statistics aimed at sorting out the correlations existing between tested enzyme variants and ability to metabolize steroid substrates. Functional classifications of both CYP1A enzyme variants and steroid substrate structures were obtained allowing the delineation of global structural features for both substrate recognition and regioselectivity of oxidation.

Ethoxyresorufin is a classical substrate for vertebrate CYP1A enzymes. In Pterygoplichthys sp. (Loricariidae) this enzyme possesses 48 amino acids substitutions compared to CYP1A sequences from other vertebrate species. These substitutions or a certain subset substitution are responsible for the non-detection of the EROD reaction in this species liver microsomes. In the present study, we investigated the catalytic activity of Pterygoplichthys sp. CYP1A toward 15 potential substrates in order to understand the substrate preferences of this modified CYP1A. The fish gene was expressed in yeast and the accumulation of the protein was confirmed by both the characteristic P450-CO absorbance spectra and by detection with monoclonal antibodies. Catalytic activities were assayed with yeast microsomes and four resorufin ethers, six coumarin derivates, three flavones, resveratrol and ethoxyfluoresceinethylester. Results demonstrated that the initial velocity pattern of this enzyme for the resorufin derivatives is different from the one described for most vertebrate CYP1As. The initial velocity for the activity with the coumarin derivatives is several orders of magnitude higher than with the resorufins, i.e. the turnover number (kcat) for ECOD is 400× higher than for EROD. Nonetheless, the specificity constant (kcat/km) for EROD is only slightly higher than for ECOD. EFEE is degraded at a rate comparable to the resorufins. Pterygoplichthys sp. CYP1A also degrades 7-methoxyflavone and β-naphthoflavone but not resveratrol and chrysin. These results indicate a divergent substrate preference for Pterygoplichthys sp. CYP1A, which may be involved in the adaptation of Loricariidae fish to their particular environment and feeding habits.