Structural aspects of ligand binding to and electron transfer in bacterial and fungal P450s

Hierarchical dynamic regulation of Saccharomyces cerevisiae for enhanced lutein biosynthesis

Lutein, as a carotenoid with strong antioxidant capacity and an important component of macular pigment in the retina, has wide applications in pharmaceutical, food, feed, and cosmetics industries. Besides extraction from plant and algae, microbial fermentation using engineered cell factories to produce lutein has emerged as a promising route. However, intra-pathway competition between the lycopene cyclases and the conflict between cell growth and production are two major challenges. In our previous study, de novo synthesis of lutein had been achieved in Saccharomyces cerevisiae by dividing the pathway into two stages (δ-carotene formation and conversion) using temperature as the input signal to realize sequential cyclation of lycopene. However, lutein production was limited to microgram level, which is still too low to meet industrial demand. In this study, a dual-signal hierarchical dynamic regulation system was developed and applied to divide lutein biosynthesis into three stages in response to glucose concentration and culture temperature. By placing the genes involved in δ-carotene formation under the glucose-responsive ADH2 promoter and genes involved in the conversion of δ-carotene to lutein under temperature-responsive GAL promoters, the growth-production conflict and intra-pathway competition were simultaneously resolved. Meanwhile, the rate-limiting lycopene ε-cyclation and carotene hydroxylation reactions were improved by screening for lycopene ε-cyclase with higher activity and fine tuning of the P450 enzymes and their redox partners. Finally, a lutein titer of 19.92 mg/L (4.53 mg/g DCW) was obtained in shake-flask cultures using the engineered yeast strain YLutein-3S-6, which is the highest lutein titer ever reported in heterologous production systems.

Controlling the Substrate Specificity of an Enzyme through Structural Flexibility by Varying the Salt-Bridge Density

Many enzymes, particularly in one single family, with highly conserved structures and folds exhibit rather distinct substrate specificities. The underlying mechanism remains elusive, the resolution of which is of great importance for biochemistry, biophysics, and bioengineering. Here, we performed a neutron scattering experiment and molecular dynamics (MD) simulations on two structurally similar CYP450 proteins; CYP101 primarily catalyzes one type of ligands, then CYP2C9 can catalyze a large range of substrates. We demonstrated that it is the high density of salt bridges in CYP101 that reduces its structural flexibility, which controls the ligand access channel and the fluctuation of the catalytic pocket, thus restricting its selection on substrates. Moreover, we performed MD simulations on 146 different kinds of CYP450 proteins, spanning distinct biological categories including Fungi, Archaea, Bacteria, Protista, Animalia, and Plantae, and found the above mechanism generally valid. We demonstrated that, by fine changes of chemistry (salt-bridge density), the CYP450 superfamily can vary the structural flexibility of its member proteins among different biological categories, and thus differentiate their substrate specificities to meet the specific biological needs. As this mechanism is well-controllable and easy to be implemented, we expect it to be generally applicable in future enzymatic engineering to develop proteins of desired substrate specificities.

Reconciling conformational heterogeneity and substrate recognition in cytochrome P450

Cytochrome P450, the ubiquitous metalloenzyme involved in detoxification of foreign components, has remained one of the most popular systems for substrate-recognition process. However, despite being known for its high substrate specificity, the mechanistic basis of substrate-binding by archetypal system cytochrome P450cam has remained at odds with the contrasting reports of multiple diverse crystallographic structures of its substrate-free form. Here, we address this issue by elucidating the probability of mutual dynamical transition to the other crystallographic pose of cytochrome P450cam and vice versa via unbiased all-atom computer simulation. A robust Markov state model, constructed using adaptively sampled 84-μs-long molecular dynamics simulation trajectories, maps the broad and heterogenous P450cam conformational landscape into five key substates. In particular, the Markov state model identifies an intermediate-assisted dynamic equilibrium between a pair of conformations of P450cam, in which the substrate-recognition sites remain "closed" and "open," respectively. However, the estimate of a significantly higher stationary population of closed conformation, coupled with faster rate of open → closed transition than its reverse process, dictates that the net conformational equilibrium would be swayed in favor of "closed" conformation. Together, the investigation quantitatively infers that although a potential substrate of cytochrome P450cam would, in principle, explore a diverse array of conformations of substrate-free protein, it would mostly encounter a "closed" or solvent-occluded conformation and hence would follow an induced-fit-based recognition process. Overall, the work reconciles multiple precedent crystallographic, spectroscopic investigations and establishes how a statistical elucidation of conformational heterogeneity in protein would provide crucial insights in the mechanism of potential substrate-recognition process.

Functional Analysis of P450 Monooxygenase SrrO in the Biosynthesis of Butenolide-Type Signaling Molecules in Streptomyces rochei

Streptomyces rochei 7434AN4 produces two structurally unrelated polyketide antibiotics lankacidin and lankamycin, and their biosynthesis is tightly controlled by butenolide-type signaling molecules SRB1 and SRB2. SRBs are synthesized by SRB synthase SrrX, and induce lankacidin and lankamycin production at 40 nM concentration. We here investigated the role of a P450 monooxygenase gene srrO (orf84), which is located adjacent to srrX (orf85), in SRB biosynthesis. An srrO mutant KA54 accumulated lankacidin and lankamycin at a normal level when compared with the parent strain. To elucidate the chemical structures of the signaling molecules accumulated in KA54 (termed as KA54-SRBs), this mutant was cultured (30 L) and the active components were purified. Two active components (KA54-SRB1 and KA54-SRB2) were detected in ESI-MS and chiral HPLC analysis. The molecular formulae for KA54-SRB1 and KA54-SRB2 are C₁₅H₂₆O₄ and C₁₆H₂₈O₄, whose values are one oxygen smaller and two hydrogen larger when compared with those for SRB1 and SRB2, respectively. Based on extensive NMR analysis, the signaling molecules in KA54 were determined to be 6'-deoxo-SRB1 and 6'-deoxo-SRB2. Gel shift analysis indicated that a ligand affinity of 6'-deoxo-SRB1 to the specific receptor SrrA was 100-fold less than that of SRB1. We performed bioconversion of the synthetic 6'-deoxo-SRB1 in the Streptomyces lividans recombinant carrying SrrO-expression plasmid. Substrate 6'-deoxo-SRB1 was converted through 6'-deoxo-6'-hydroxy-SRB1 to SRB1 in a time-dependent manner. Thus, these results clearly indicated that SrrO catalyzes the C-6' oxidation at a final step in SRB biosynthesis.

Structural basis for plant lutein biosynthesis from α-carotene

Two cytochrome P450 enzymes, CYP97A3 and CYP97C1, catalyze hydroxylations of the β- and ε-rings of α-carotene to produce lutein. Chirality is introduced at the C-3 atom of both rings, and the reactions are both pro-3R-stereospecific. We determined the crystal structures of CYP97A3 in substrate-free and complex forms with a nonnatural substrate and the structure of CYP97C1 in a detergent-bound form. The structures of CYP97A3 in different states show the substrate channel and the structure of CYP97C1 bound with octylthioglucoside confirms the binding site for the carotenoid substrate. Biochemical assays confirm that the ferredoxin-NADP⁺ reductase (FNR)-ferredoxin pair is used as the redox partner. Details of the pro-3R stereospecificity are revealed in the retinal-bound CYP97A3 structure. Further analysis indicates that the CYP97B clan bears similarity to the β-ring-specific CYP97A clan. Overall, our research describes the molecular basis for the last steps of lutein biosynthesis.

A Continuing Career in Biocatalysis: Frances H. Arnold

On the occasion of Professor Frances H. Arnold's recent acceptance of the 2018 Nobel Prize in Chemistry, we honor her numerous contributions to the fields of directed evolution and biocatalysis. Arnold pioneered the development of directed evolution methods for engineering enzymes as biocatalysts. Her highly interdisciplinary research has provided a ground not only for understanding the mechanisms of enzyme evolution but also for developing commercially viable enzyme biocatalysts and biocatalytic processes. In this Account, we highlight some of her notable contributions in the past three decades in the development of foundational directed evolution methods and their applications in the design and engineering of enzymes with desired functions for biocatalysis. Her work has created a paradigm shift in the broad catalysis field.

Conformational Changes in Cytochrome P450cam and the Effector Role of Putidaredoxin

The cytochromes P450 form an enormous family of over 20 000 enzyme variants found in all branches of life. They catalyze the O2 dependent monooxygenation of a wide range of substrates in reactions important to drug metabolism, biosynthesis and energy utilization. Understanding how they function is important for biomedical science and requires a full description of their notorious propensity for specificity and promiscuity. The bacterial P450cam is an unusual example, having the most well characterized chemical mechanism of all of the forms. It also undergoes an increasingly well characterized structural change upon substrate binding, which may be similar to to that displayed by some, but not all forms of P450. Finally, P450cam is one of the rare forms that have a strict requirement for a particular electron donor, putidaredoxin (pdx). Pdx provides the required electrons for enzyme turnover, but it also induces specific changes in the enzyme to allow enzyme turnover, long known as its effector role. This review summarizes recent crystallographic and double electron–electron resonance studies that have revealed the effects of substrate and pdx binding on the structure of P450cam. We describe an emerging idea for how pdx exerts its effector function by inducing a conformational change in the enzyme. This change then propagates to the active site to enable cleavage of the ferric–hydroperoxy bond during catalysis, and appears to provide a very elegant approach for P450cam to attain both high efficiency and protection from oxidative damage.

Effector Roles of Putidaredoxin on Cytochrome P450cam Conformational States

In this study, the effector role of Pdx (putidaredoxin) on cytochrome P450cam conformation is refined by attaching two different spin labels, MTSL or BSL (bifunctional spin-label) onto the F or G helices and using DEER (double electron-electron resonance) to measure the distance between labels. Recent EPR and crystallographic studies have observed that oxidized Pdx induces substrate-bound P450cam to change from the closed to the open state. However, this change was not observed by DEER in the reduced Pdx complex with carbon-monoxide-bound P450cam (Fe(2+)CO). In addition, recent NMR studies have failed to observe a change in P450cam conformation upon binding Pdx. Hence, resolving these issues is important for a full understanding the effector role of Pdx. Here we show that oxidized Pdx induces camphor-bound P450cam to shift from the closed to the open conformation when labeled on either the F or G helices with MTSL. BSL at these sites can either narrow the distance distribution widths dramatically or alter the extent of the conformational change. In addition, we report DEER spectra on a mixed oxidation state containing oxidized Pdx and ferrous CO-bound P450cam, showing that P450cam remains closed. This indicates that CO binding to the heme prevents P450cam from opening, overriding the influence exerted by Pdx binding. Finally, we report the open form P450cam crystal structure with substrate bound, which suggests that crystal packing effects may prevent conformational conversion. Using multiple labeling approaches, DEER provides a unique perspective to resolve how the conformation of P450cam depends on Pdx and ligand states.

Mechanistic study of allopurinol oxidation using aldehyde oxidase, xanthine oxidase and cytochrome P450 enzymes

The oxidation reaction of allopurinol to its active metabolite (oxypurinol) is investigated using the AO and P450 enzymes. To the contrary of AO (and XO), the P450 enzyme can metabolize the allopurinol with a not self-inhibitory mechanism.

Coupling Oxygen Consumption with Hydrocarbon Oxidation in Bacterial Multicomponent Monooxygenases

A fundamental goal in catalysis is the coupling of multiple reactions to yield a desired product. Enzymes have evolved elegant approaches to address this grand challenge. A salient example is the biological conversion of methane to methanol catalyzed by soluble methane monooxygenase (sMMO), a member of the bacterial multicomponent monooxygenase (BMM) superfamily. sMMO is a dynamic protein complex of three components: a hydroxylase, a reductase, and a regulatory protein. The active site, a carboxylate-rich non-heme diiron center, is buried inside the 251 kDa hydroxylase component. The enzyme processes four substrates: O2, protons, electrons, and methane. To couple O2 activation to methane oxidation, timely control of substrate access to the active site is critical. Recent studies of sMMO, as well as its homologues in the BMM superfamily, have begun to unravel the mechanism. The emerging and unifying picture reveals that each substrate gains access to the active site along a specific pathway through the hydroxylase. Electrons and protons are delivered via a three-amino-acid pore located adjacent to the diiron center; O2 migrates via a series of hydrophobic cavities; and hydrocarbon substrates reach the active site through a channel or linked set of cavities. The gating of these pathways mediates entry of each substrate to the diiron active site in a timed sequence and is coordinated by dynamic interactions with the other component proteins. The result is coupling of dioxygen consumption with hydrocarbon oxidation, avoiding unproductive oxidation of the reductant rather than the desired hydrocarbon. To initiate catalysis, the reductase delivers two electrons to the diiron(III) center by binding over the pore of the hydroxylase. The regulatory component then displaces the reductase, docking onto the same surface of the hydroxylase. Formation of the hydroxylase-regulatory component complex (i) induces conformational changes of pore residues that may bring protons to the active site; (ii) connects hydrophobic cavities in the hydroxylase leading from the exterior to the diiron active site, providing a pathway for O2 and methane, in the case of sMMO, to the reduced diiron center for O2 activation and substrate hydroxylation; (iii) closes the pore, as well as a channel in the case of four-component BMM enzymes, restricting proton access to the diiron center during formation of "Fe2O2" intermediates required for hydrocarbon oxidation; and (iv) inhibits undesired electron transfer to the Fe2O2 intermediates by blocking reductase binding during O2 activation. This mechanism is quite different from that adopted by cytochromes P450, a large class of heme-containing monooxygenases that catalyze reactions very similar to those catalyzed by the BMM enzymes. Understanding the timed enzyme control of substrate access has implications for designing artificial catalysts. To achieve multiple turnovers and tight coupling, synthetic models must also control substrate access, a major challenge considering that nature requires large, multimeric, dynamic protein complexes to accomplish this feat.

Structure, dynamics, and function of the monooxygenase P450 BM-3: insights from computer simulations studies

The monooxygenase P450 BM-3 is a NADPH-dependent fatty acid hydroxylase enzyme isolated from soil bacterium Bacillus megaterium. As a pivotal member of cytochrome P450 superfamily, it has been intensely studied for the comprehension of structure-dynamics-function relationships in this class of enzymes. In addition, due to its peculiar properties, it is also a promising enzyme for biochemical and biomedical applications. However, despite the efforts, the full understanding of the enzyme structure and dynamics is not yet achieved. Computational studies, particularly molecular dynamics (MD) simulations, have importantly contributed to this endeavor by providing new insights at an atomic level regarding the correlations between structure, dynamics, and function of the protein. This topical review summarizes computational studies based on MD simulations of the cytochrome P450 BM-3 and gives an outlook on future directions.

Evaluation of Luminescent P450 Analysis for Directed Evolution of Human CYP4A11

Cytochrome P450 4A11 (CYP4A11) is a fatty acid hydroxylase enzyme expressed in human liver. It catalyzes not only the hydroxylation of saturated and unsaturated fatty acids, but the conversion of arachidonic acid to 20-hydroxyeicosatetraenoic acid (20-HETE), a regulator of blood pressure. In this study, we performed a directed evolution analysis of CYP4A11 using the luminogenic assay system. A random mutant library of CYP4A11, in which mutations were made throughout the entire coding region, was screened with luciferase activity to detect the demethylation of luciferin-4A (2-[6-methoxyquinolin-2-yl]-4,5-dihydrothiazole-4-carboxylic acid) of CYP4A11 mutants in Escherichia coli. Consecutive rounds of random mutagenesis and screening yielded three improved CYP4A11 mutants, CP2600 (A24T/T263A), CP2601 (T263A), and CP2616 (A24T/T263A/V430E) with ~3-fold increase in whole cells and >10-fold increase in purified proteins on the luminescence assay. However, the steady state kinetic analysis for lauric acid hydroxylation showed the significant reductions in enzymatic activities in all three mutants. A mutant, CP2600, showed a 51% decrease in catalytic efficiency (k_cat/K_m) for lauric acid hydroxylation mainly due to an increase in K_m. CP2601 and CP2616 showed much greater reductions (>75%) in the catalytic efficiency due to both a decrease in k_cat and an increase in K_m. These decreased catalytic activities of CP2601 and CP2616 can be partially attributed to the changes in substrate affinities. These results suggest that the enzymatic activities of CYP4A11 mutants selected from directed evolution using a luminogenic P450 substrate may not demonstrate a direct correlation with the hydroxylation activities of lauric acid.

Controlled oxidation of remote sp3 C-H bonds in artemisinin via P450 catalysts with fine-tuned regio- and stereoselectivity

The selective oxyfunctionalization of isolated sp(3) C-H bonds in complex molecules represents a formidable challenge in organic chemistry. Here, we describe a rational, systematic strategy to expedite the development of P450 oxidation catalysts with refined regio- and stereoselectivity for the hydroxylation of remote, unactivated C-H sites in a complex scaffold. Using artemisinin as model substrate, we demonstrate how a three-tier strategy involving first-sphere active site mutagenesis, high-throughput P450 fingerprinting, and fingerprint-driven P450 reactivity predictions enabled the rapid evolution of three efficient biocatalysts for the selective hydroxylation of a primary and a secondary C-H site (with both S and R stereoselectivity) in a relevant yet previously inaccessible region of this complex natural product. The evolved P450 variants could be applied to provide direct access to the desired hydroxylated derivatives at preparative scales (0.4 g) and in high isolated yields (>90%), thereby enabling further elaboration of this molecule. As an example, enantiopure C7-fluorinated derivatives of the clinical antimalarial drugs artesunate and artemether, in which a major metabolically sensitive site is protected by means of a C-H to C-F substitution, were afforded via P450-mediated chemoenzymatic synthesis.

Double electron-electron resonance shows cytochrome P450cam undergoes a conformational change in solution upon binding substrate

Although cytochrome P450cam from Pseudomonas putida, the archetype for all heme monooxygenases, has long been known to have a closed active site, recent reports show that the enzyme can also be crystallized in at least two clusters of open conformations. This suggests that the enzyme may undergo significant conformational changes during substrate binding and catalytic turnover. However, these conformations were observed in the crystalline state, and information is needed about the conformations that are populated in solution. In this study, double electron-electron resonance experiments were performed to observe substrate-induced changes in distance as measured by the dipolar coupling between spin labels introduced onto the surface of the enzyme on opposite sides of the substrate access channel. The double electron-electron resonance data show a decrease of 0.8 nm in the distance between spin labels placed at S48C and S190C upon binding the substrate camphor. A rotamer distribution model based on the crystal structures adequately describes the observed distance distributions. These results demonstrate conclusively that, in the physiologically relevant solution state, the substrate-free enzyme exists in the open P450cam-O conformation and that camphor binding results in conversion to the closed P450cam-C form. This approach should be useful for investigating many other P450s, including mammalian forms, in which the role of conformational change is of central importance but not well understood.

A novel type of allosteric regulation: functional cooperativity in monomeric proteins

Cooperative functional properties and allosteric regulation in cytochromes P450 play an important role in xenobiotic metabolism and define one of the main mechanisms of drug-drug interactions. Recent experimental results suggest that ability to bind simultaneously two or more small organic molecules can be the essential feature of cytochrome P450 fold, and often results in rich and complex pattern of allosteric behavior. Manifestations of non-Michaelis kinetics include homotropic and heterotropic activation and inhibition effects depending on the stoichiometric ratios of substrate and effector, changes in the regio- and stereospecificity of catalytic transformations, and often give rise to the clinically important drug-drug interactions. In addition, functional response of P450 systems is modulated by the presence of specific and non-specific effector molecules, metal ions, membrane incorporation, formation of homo- and hetero-oligomers, and interactions with the protein redox partners. In this article we briefly overview the main factors contributing to the allosteric effects in cytochromes P450 with the main focus on the sources of cooperative behavior in xenobiotic metabolizing monomeric heme enzymes with their conformational flexibility and extremely broad substrate specificity. The novel mechanism of functional cooperativity in P450 enzymes does not require substantial binding cooperativity, rather it implies the presence of one or more binding sites with higher affinity than the single catalytically active site in the vicinity of the heme iron.

Structure and function of CYP108D1 from Novosphingobium aromaticivorans DSM12444: an aromatic hydrocarbon-binding P450 enzyme

CYP108D1 from Novosphingobium aromaticivorans DSM12444 binds a range of aromatic hydrocarbons such as phenanthrene, biphenyl and phenylcyclohexane. Its structure, which is reported here at 2.2 Å resolution, is closely related to that of CYP108A1 (P450terp), an α-terpineol-oxidizing enzyme. The compositions and structures of the active sites of these two enzymes are very similar; the most significant changes are the replacement of Glu77 and Thr103 in CYP108A1 by Thr79 and Val105 in CYP108D1. Other residue differences lead to a larger and more hydrophobic access channel in CYP108D1. These structural features are likely to account for the weaker α-terpineol binding by CYP108D1 and, when combined with the presence of three hydrophobic phenylalanine residues in the active site, promote the binding of aromatic hydrocarbons. The haem-proximal surface of CYP108D1 shows a different charge distribution and topology to those of CYP101D1, CYP101A1 and CYP108A1, including a pronounced kink in the proximal loop of CYP108D1, which may result in poor complementarity with the [2Fe-2S] ferredoxins Arx, putidaredoxin and terpredoxin that are the respective redox partners of these three P450 enzymes. The unexpectedly low reduction potential of phenylcyclohexane-bound CYP108D1 (-401 mV) may also contribute to the low activity observed with these ferredoxins. CYP108D1 appears to function as an aromatic hydrocarbon hydroxylase that requires a different electron-transfer cofactor protein.

Miniaturized, microarray-based assays for chemical proteomic studies of protein function

Systematic analysis of protein and enzyme function typically requires scale-up of protein expression and purification prior to assay development; this can often be limiting. Miniaturization of assays provides an alternative approach, but simple, generic methods are in short supply. Here we show how custom microarrays can be adapted to this purpose. We discuss the different routes to array fabrication and describe in detail one facile approach in which the purification and immobilization procedures are combined into a single step, significantly simplifying the array fabrication process. We illustrate this approach by reference to the creation of arrays of human protein kinases and of human cytochrome P450s. We discuss methods for both ligand-binding and turnover-based assays, as well as data analysis on such arrays.

Structural differences between soluble and membrane bound cytochrome P450s

The superfamily of cytochrome P450s forms a large class of heme monooxygenases with more than 13,000 enzymes represented in organisms from all biological kingdoms. Despite impressive variability in sizes, sequences, location, and function, all cytochrome P450s from various organisms have very similar tertiary structures within the same fold. Here we show that systematic comparison of all available X-ray structures of cytochrome P450s reveals the presence of two distinct structural classes of cytochrome P450s. For all membrane bound enzymes, except the CYP51 family, the beta-domain and the A-propionate heme side chain are shifted towards the proximal side of the heme plane, which may result in an increase of the volume of the substrate binding pocket and an opening of a potential channel for the substrate access and/or product escape directly into the membrane. This structural feature is also observed in several soluble cytochrome P450s, such as CYP108, CYP151, and CYP158A2, which catalyze transformations of bulky substrates. Alternatively, both beta-domains and the A-propionate side chains in the soluble isozymes extend towards the distal site of the heme. This difference between the structures of soluble and membrane bound cytochrome P450s can be rationalized through the presence of several amino acid inserts in the latter class which are involved in direct interactions with the membrane, namely the F'- and G'-helices. Molecular dynamics using the most abundant human cytochrome P450, CYP3A4, incorporated into a model POPC bilayer reveals the facile conservation of a substrate access channel, directed into the membrane between the B-C loop and the beta domain, and the closure of the peripheral substrate access channel directed through the B-C loop. This is in contrast to the case when the same simulation is run in buffer, where no such channel closing occurs. Taken together, these results reveal a key structural difference between membrane bound and soluble cytochrome P450s with important functional implications induced by the lipid bilayer.

P450(BM3) (CYP102A1): connecting the dots

P450(BM3) (CYP102A1), a fatty acid hydroxylase from Bacillus megaterium, has been extensively studied over a period of almost forty years. The enzyme has been redesigned to catalyse the oxidation of non-natural substrates as diverse as pharmaceuticals, terpenes and gaseous alkanes using a variety of engineering strategies. Crystal structures have provided a basis for several of the catalytic effects brought about by mutagenesis, while changes to reduction potentials, inter-domain electron transfer rates and catalytic parameters have yielded functional insights. Areas of active research interest include drug metabolite production, the development of process-scale techniques, unravelling general mechanistic aspects of P450 chemistry, methane oxidation, and improving selectivity control to allow the synthesis of fine chemicals. This review draws together the disparate research themes and places them in a historical context with the aim of creating a resource that can be used as a gateway to the field.

Three clusters of conformational states in p450cam reveal a multistep pathway for closing of the substrate access channel

Conformational changes in the substrate access channel have been observed for several forms of cytochrome P450, but the extent of conformational plasticity exhibited by a given isozyme has not been completely characterized. Here we present crystal structures of P450cam bound to a library of 12 active site probes containing a substrate analogue tethered to a variable linker. The structures provide a unique view of the range of protein conformations accessible during substrate binding. Principal component analysis of a total of 30 structures reveals three discrete clusters of conformations: closed (P450cam-C), intermediate (P450cam-I), and fully open (P450cam-O). Relative to P450cam-C, the P450cam-I state results predominantly from a retraction of helix F, while both helices F and G move in concert to reach the fully open P450cam-O state. Both P450cam-C and P450cam-I are well-defined states, while P450cam-O shows evidence of a somewhat broader distribution of conformations and includes the open form recently seen in the absence of substrate. The observed clustering of protein conformations over a wide range of ligand variants suggests a multistep closure of the enzyme around the substrate that begins by conformational selection from an ensemble of open conformations and proceeds through a well-defined intermediate, P450cam-I, before full closure to the P450cam-C state in the presence of small substrates. This multistep pathway may have significant implications for a full understanding of substrate specificity, kinetics, and coupling of substrate binding to P450 function.

P450cam visits an open conformation in the absence of substrate

P450cam from Pseudomonas putida is the best characterized member of the vast family of cytochrome P450s, and it has long been believed to have a more rigid and closed active site relative to other P450s. Here we report X-ray structures of P450cam crystallized in the absence of substrate and at high and low [K(+)]. The camphor-free structures are observed in a distinct open conformation characterized by a water-filled channel created by the retraction of the F and G helices, disorder of the B' helix, and loss of the K(+) binding site. Crystallization in the presence of K(+) alone does not alter the open conformation, while crystallization with camphor alone is sufficient for closure of the channel. Soaking crystals of the open conformation in excess camphor does not promote camphor binding or closure, suggesting resistance to conformational change by the crystal lattice. This open conformation is remarkably similar to that seen upon binding large tethered substrates, showing that it is not the result of a perturbation by the ligand. Redissolved crystals of the open conformation are observed as a mixture of P420 and P450 forms, which is converted to the P450 form upon addition of camphor and K(+). These data reveal that P450cam can dynamically visit an open conformation that allows access to the deeply buried active site without being induced by substrate or ligand.

Investigating the structural plasticity of a cytochrome P450: three-dimensional structures of P450 EryK and binding to its physiological substrate

Cytochrome P450s are heme-containing proteins that catalyze the oxidative metabolism of many physiological endogenous compounds. Because of their unique oxygen chemistry and their key role in drug and xenobiotic metabolism, particular attention has been devoted in elucidating their mechanism of substrate recognition. In this work, we analyzed the three-dimensional structures of a monomeric cytochrome P450 from Saccharopolyspora erythraea, commonly called EryK, and the binding kinetics to its physiological ligand, erythromycin D. Three different structures of EryK were obtained: two ligand-free forms and one in complex with its substrate. Analysis of the substrate-bound structure revealed the key structural determinants involved in substrate recognition and selectivity. Interestingly, the ligand-free structures of EryK suggested that the protein may explore an open and a closed conformation in the absence of substrate. In an effort to validate this hypothesis and to investigate the energetics between such alternative conformations, we performed stopped-flow absorbance experiments. Data demonstrated that EryK binds erythromycin D via a mechanism involving at least two steps. Contrary to previously characterized cytochrome P450s, analysis of double jump mixing experiments confirmed that this complex scenario arises from a pre-existing equilibrium between the open and closed subpopulations of EryK, rather than from an induced-fit type mechanism.

Probing intermediates in the activation cycle of [NiFe] hydrogenase by infrared spectroscopy: the Ni-SIr state and its light sensitivity

The [NiFe] hydrogenase from the sulphate-reducing bacterium Desulfovibrio vulgaris Miyazaki F is reversibly inhibited in the presence of molecular oxygen. A key intermediate in the reactivation process, Ni-SI_r, provides the link between fully oxidized (Ni-A, Ni-B) and active (Ni-SI_a, Ni-C and Ni-R) forms of hydrogenase. In this work Ni-SI_r was found to be light-sensitive (T ≤ 110 K), similar to the active Ni-C and the CO-inhibited states. Transition to the final photoproduct state (Ni-SL) was shown to involve an additional transient light-induced state (Ni-SI₁₉₆₁). Rapid scan kinetic infrared measurements provided activation energies for the transition from Ni-SL to Ni-SI_r in protonated as well as in deuterated samples. The inhibitor CO was found not to react with the active site of the Ni-SL state. The wavelength dependence of the Ni-SI_r photoconversion was examined in the range between 410 and 680 nm. Light-induced effects were associated with a nickel-centred electronic transition, possibly involving a change in the spin state of nickel (Ni²⁺). In addition, at T ≤ 40 K the CN⁻ stretching vibrations of Ni-SL were found to be dependent on the colour of the monochromatic light used to irradiate the species, suggesting a change in the interaction of the hydrogen-bonding network of the surrounding amino acids. A possible mechanism for the photochemical process, involving displacement of the oxygen-based ligand, is discussed.

Cooperative properties of cytochromes P450

Cytochromes P450 form a large and important class of heme monooxygenases with a broad spectrum of substrates and corresponding functions, from steroid hormone biosynthesis to the metabolism of xenobiotics. Despite decades of study, the molecular mechanisms responsible for the complex non-Michaelis behavior observed with many members of this superfamily during metabolism, often termed 'cooperativity', remain to be fully elucidated. Although there is evidence that oligomerization may play an important role in defining the observed cooperativity, some monomeric cytochromes P450, particularly those involved in xenobiotic metabolism, also display this behavior due to their ability to simultaneously bind several substrate molecules. As a result, formation of distinct enzyme-substrate complexes with different stoichiometry and functional properties can give rise to homotropic and heterotropic cooperative behavior. This review aims to summarize the current understanding of cooperativity in cytochromes P450, with a focus on the nature of cooperative effects in monomeric enzymes.

Structural characterization of CalO2: a putative orsellinic acid P450 oxidase in the calicheamicin biosynthetic pathway

Although bacterial iterative Type I polyketide synthases are now known to participate in the biosynthesis of a small set of diverse natural products, the subsequent downstream modification of the resulting polyketide products remains poorly understood. Toward this goal, we report the X-ray structure determination at 2.5 A resolution and preliminary characterization of the putative orsellenic acid P450 oxidase (CalO2) involved in calicheamicin biosynthesis. These studies represent the first crystal structure for a P450 involved in modifying a bacterial iterative Type I polyketide product and suggest the CalO2-catalyzed step may occur after CalO3-catalyzed iodination and may also require a coenzyme A- (CoA) or acyl carrier protein- (ACP) bound substrate. Docking studies also reveal a putative docking site within CalO2 for the CLM orsellinic acid synthase (CalO5) ACP domain which involves a well-ordered helix along the CalO2 active site cavity that is unique compared with other P450 structures.

Evolutionary history of a specialized p450 propane monooxygenase

The evolutionary pressures that shaped the specificity and catalytic efficiency of enzymes can only be speculated. While directed evolution experiments show that new functions can be acquired under positive selection with few mutations, the role of negative selection in eliminating undesired activities and achieving high specificity remains unclear. Here we examine intermediates along the 'lineage' from a naturally occurring C12-C20 fatty acid hydroxylase (P450BM3) to a laboratory-evolved P450 propane monooxygenase (P450PMO) having 20 heme domain substitutions compared to P450BM3. Biochemical, crystallographic, and computational analyses show that a minimal perturbation of the P450BM3 fold and substrate-binding pocket accompanies a significant broadening of enzyme substrate range and the emergence of propane activity. In contrast, refinement of the enzyme catalytic efficiency for propane oxidation (approximately 9000-fold increase in kcat/Km) involves profound reshaping and partitioning of the substrate access pathway. Remodeling of the substrate-recognition mechanisms ultimately results in remarkable narrowing of the substrate profile around propane and enables the acquisition of a basal iodomethane dehalogenase activity as yet unknown in natural alkane monooxygenases. A highly destabilizing L188P substitution in a region of the enzyme that undergoes a large conformational change during catalysis plays an important role in adaptation to the gaseous alkane. This work demonstrates that positive selection alone is sufficient to completely respecialize the cytochrome P450 for function on a nonnative substrate.

Active-site structure, binding and redox activity of the heme-thiolate enzyme CYP2D6 immobilized on coated Ag electrodes: a surface-enhanced resonance Raman scattering study

Surface-enhance resonance Raman scattering spectra of the heme-thiolate enzyme cytochrome P450 2D6 (CYP2D6) adsorbed on Ag electrodes coated with 11-mercaptoundecanoic acid (MUA) were obtained in various experimental conditions. An analysis of these spectra, and a comparison between them and the RR spectra of CYP2D6 in solution, indicated that the enzyme's active site retained its nature of six-coordinated low-spin heme upon immobilization. Moreover, the spectral changes detected in the presence of dextromethorphan (a CYP2D6 substrate) and imidazole (an exogenous heme axial ligand) indicated that the immobilized enzyme also preserved its ability to reversibly bind a substrate and form a heme-imidazole complex. The reversibility of these processes could be easily verified by flowing alternately solutions of the various compounds and the buffer through a home-built spectroelectrochemical flow cell which contained a sample of immobilized protein, without the need to disassemble the cell between consecutive spectral data acquisitions. Despite immobilized CYP2D6 being effectively reduced by a sodium dithionite solution, electrochemical reduction via the Ag electrode was not able to completely reduce the enzyme, and led to its extensive inactivation. This behavior indicated that although the enzyme's ability to exchange electrons is not altered by immobilization per se, MUA-coated electrodes are not suited to perform direct electrochemistry of CYP2D6.

Filling a Hole in Cytochrome P450 BM3 Improves Substrate Binding and Catalytic Efficiency

Cytochrome P450BM3 (CYP102A1) from Bacillus megaterium, a fatty acid hydroxylase, is a member of a very large superfamily of monooxygenase enzymes. The available crystal structures of the enzyme show non-productive binding of substrates with their omega-end distant from the iron in a hydrophobic pocket at one side of the active site. We have constructed and characterised mutants in which this pocket is filled by large hydrophobic side-chains replacing alanine at position 82. The mutants having phenylalanine or tryptophan at this position have very much (approximately 800-fold) greater affinity for substrate, with a greater conversion of the haem iron to the high-spin state, and similarly increased catalytic efficiency. The enzyme as isolated contains bound palmitate, reflecting this much higher affinity. We have determined the crystal structure of the haem domain of the Ala82Phe mutant with bound palmitate; this shows that the substrate is binding differently from the wild-type enzyme but still distant from the haem iron. Detailed analysis of the structure indicates that the tighter binding in the mutant reflects a shift in the conformational equilibrium of the substrate-free enzyme towards the conformation seen in the substrate complex rather than differences in the enzyme-substrate interactions. On this basis, we outline a sequence of events for the initial stages of the catalytic cycle. The Ala82Phe and Ala82Trp mutants are also very much more effective catalysts of indole hydroxylation than the wild-type enzyme, suggesting that they will be valuable starting points for the design of mutants to catalyse synthetically useful hydroxylation reactions.

Conformational dynamics in the F/G segment of CYP51 from Mycobacterium tuberculosis monitored by FRET

A cysteine was introduced into the FG-loop (P187C) of CYP51 from Mycobacterium tuberculosis (MT) for selective labeling with BODIPY and fluorescence energy transfer (FRET) analysis. Förster radius for the BODIPY-heme pair was calculated assuming that the distance between the heme and Cys187 in solution corresponds to that in the crystal structure of ligand free MTCYP51. Interaction of MTCYP51 with azole inhibitors ketoconazole and fluconazole or the substrate analog estriol did not influence the fluorescence, but titration with the substrate lanosterol quenched BODIPY emission, the effect being proportional to the portion of substrate bound MTCYP51. The detected changes correspond to approximately 10A decrease in the calculated distance between BODIPY-Cys187 and the heme. The results confirm (1) functional importance of conformational motions in the MTCYP51 F/G segment and (2) applicability of FRET to monitor them in solution.

Neutral genetic drift can alter promiscuous protein functions, potentially aiding functional evolution

Background

Many of the mutations accumulated by naturally evolving proteins are neutral in the sense that they do not significantly alter a protein's ability to perform its primary biological function. However, new protein functions evolve when selection begins to favor other, "promiscuous" functions that are incidental to a protein's original biological role. If mutations that are neutral with respect to a protein's primary biological function cause substantial changes in promiscuous functions, these mutations could enable future functional evolution.

Results

Here we investigate this possibility experimentally by examining how cytochrome P450 enzymes that have evolved neutrally with respect to activity on a single substrate have changed in their abilities to catalyze reactions on five other substrates. We find that the enzymes have sometimes changed as much as four-fold in the promiscuous activities. The changes in promiscuous activities tend to increase with the number of mutations, and can be largely rationalized in terms of the chemical structures of the substrates. The activities on chemically similar substrates tend to change in a coordinated fashion, potentially providing a route for systematically predicting the change in one activity based on the measurement of several others.

Conclusion

Our work suggests that initially neutral genetic drift can lead to substantial changes in protein functions that are not currently under selection, in effect poising the proteins to more readily undergo functional evolution should selection favor new functions in the future.

Reviewers

This article was reviewed by Martijn Huynen, Fyodor Kondrashov, and Dan Tawfik (nominated by Christoph Adami).

Conformational dynamics of the cytochrome P450 BM3/N-palmitoylglycine complex: the proposed "proximal-distal" transition probed by temperature-jump spectroscopy

The ferric spin state equilibrium of the heme iron was analyzed in wild-type cytochrome P450 BM3 and its F87G mutant by using temperature (T)-jump relaxation spectroscopy in combination with static equilibrium experiments. No relaxation process was measurable in the substrate-free enzyme indicating a relaxation process with a rate constant>10,000 s(-1). In contrast, a slow spin state transition process was observed in the N-palmitoylglycine (NPG)-bound enzyme species. This transition occurred with an observed rate constant (298 K) of approximately 800 s(-1) in the wild-type, and approximately 2500 s(-1) in the F87G mutant, suggesting a significant contribution of the phenylalanine side chain to a reaction step rate limiting the actual spin state transition. These findings are discussed in terms of an equilibrium between different binding modes of the substrate, including a position 7.5 A away from the heme iron ("distal") and the catalytically relevant "proximal" binding site, and are in accordance with results from X-ray crystallography, NMR studies, and molecular dynamics simulations.

Oxidative Phenol Coupling Reactions Catalyzed by OxyB: A Cytochrome P450 from the Vancomycin Producing Organism. Implications for Vancomycin Biosynthesis

OxyB is a cytochrome P450 enzyme that catalyzes the first phenol coupling reaction during the biosynthesis of vancomycin-like glycopeptide antibiotics. The phenol coupling reaction occurs on a linear peptide intermediate linked as a C-terminal thioester to a peptide carrier protein (PCP) domain within the multidomain glycopeptide nonribosomal peptide synthetase (NRPS). Using model peptides with the sequence (R)(NMe)Leu-(R)Tyr-(S)Asn-(R)Hpg-(R)Hpg-(S)Tyr-S-PCP and (R)(NMe)Leu-(R)Tyr-(S)Asn-(R)Hpg-(R)Hpg-(S)Tyr-(S)Dpg-S-PCP (where Hpg = 4-hydroxyphenylglycine, and Dpg = 3,5-dihydroxyphenylglycine), or containing (R)Leu instead of (R)(NMe)Leu, attached to recombinant PCPs derived from modules-6 and -7 in the vancomycin NRPS, we show that cross-linking of Hpg4 and Tyr6 by OxyB can occur in both hexapeptide- and heptapeptide-PCP conjugates. Thus, whereas OxyB may act preferentially on a hexapeptide still linked to the PCP-6 in NRPS subunit-2, it is possible that a linear heptapeptide intermediate linked to PCP-7 in NRPS subunit-3 may also be transformed into monocyclic product. For turnover, OxyB requires electrons, which in vitro can be supplied by spinach ferredoxin and E. coli flavodoxin reductase. Turnover is also dependent upon the presence of molecular oxygen. The model substrate (R)(NMe)Leu-(R)Tyr-(S)Asn-(R)Hpg-(R)Hpg-(S)Tyr-S-PCP is transformed into cross-linked product by OxyB with a kcat of 0.1 s-1 and Km in the range 4-13 muM. Equilibrium binding of this substrate to OxyB, monitored by UV-vis, is accompanied by a typical low-to-high spin state change in the heme, characterized with a Kd of 17 +/- 5 muM.

CO migration pathways in cytochrome P450cam studied by molecular dynamics simulations

Previous laser flash photolysis investigations between 100 and 300 K have shown that the kinetics of CO rebinding with cytochrome P450(cam)(camphor) consist of up to four different processes revealing a complex internal dynamics after ligand dissociation. In the present work, molecular dynamics simulations were undertaken on the ternary complex P450(cam)(cam)(CO) to explore the CO migration pathways, monitor the internal cavities of the protein, and localize the CO docking sites. One trajectory of 1 nsec with the protein in a water box and 36 trajectories of 1 nsec in the vacuum were calculated. In each trajectory, the protein contained only one CO ligand on which no constraints were applied. The simulations were performed at 200, 300, and 320 K. The results indicate the presence of seven CO docking sites, mainly hydrophobic, located in the same moiety of the protein. Two of them coincide with xenon binding sites identified by crystallography. The protein matrix exhibits eight persistent internal cavities, four of which corresponding to the ligand docking sites. In addition, it was observed that water molecules entering the protein were mainly attracted into the polar pockets, far away from the CO docking sites. Finally, the identified CO migration pathways provide a consistent interpretation of the experimental rebinding kinetics.

Combining substrate dynamics, binding statistics, and energy barriers to rationalize regioselective hydroxylation of octane and lauric acid by CYP102A1 and mutants

Hydroxylations of octane and lauric acid by Cytochrome P450-BM3 (CYP102A1) wild-type and three active site mutants--F87A, L188Q/A74G, and F87V/L188Q/A74G--were rationalized using a combination of substrate orientation from docking, substrate binding statistics from molecular dynamics simulations, and barrier energies for hydrogen atom abstraction from quantum mechanical calculations. Wild-type BM3 typically hydroxylates medium- to long-chain fatty acids on subterminal (omega-1, omega-2, omega-3) but not the terminal (omega) positions. The known carboxylic anchoring site Y51/R47 for lauric acid, and hydrophobic interactions and steric exclusion, mainly by F87, for octane as well as lauric acid, play a role in the binding modes of the substrates. Electrostatic interactions between the protein and the substrate strongly modulate the substrate's regiodependent activation barriers. A combination of the binding statistics and the activation barriers of hydrogen-atom abstraction in the substrates is proposed to determine the product formation. Trends observed in experimental product formation for octane and lauric acid by wild-type BM3 and the three active site mutants were qualitatively explained. It is concluded that the combination of substrate binding statistics and hydrogen-atom abstraction barrier energies is a valuable tool to rationalize substrate binding and product formation and constitutes an important step toward prediction of product ratios.