Structural and Spectroscopic Characterization of P450 BM3 Mutants with Unprecedented P450 Heme Iron Ligand Sets

Applications of the Newly Developed Force-Field Parameters Uncover a Dynamic Nature of Ω-Loop C in the Lys-Ligated Alkaline Form of Cytochrome c

Lys-ligated cytochromes make up an emerging family of heme proteins. Density functional theory calculations on the amine/imidazole-ligated c-type ferric heme were employed to develop force-field parameters for molecular dynamics (MD) simulations of structural and dynamic features of these proteins. The new force-field parameters were applied to the alkaline form of yeast iso-1 cytochrome c to rationalize discrepancies resulting from distinct experimental conditions in prior structural studies and to provide insights into the mechanisms of the alkaline transition. Our simulations have revealed the dynamic nature of Ω-loop C in the Lys-ligated protein and its unfolding in the Lys-ligated conformer having this loop in the same position as in the native Met-ligated protein. The proximity of Tyr67 or Tyr74 to the Lys ligand of ferric heme iron suggests a possible mechanism of the backward alkaline transition where a proton donor Tyr assists in Lys dissociation. The developed force-field parameters will be useful in structural and dynamic characterization of other native or engineered Lys-ligated heme proteins.

Directed evolution of enzymatic silicon-carbon bond cleavage in siloxanes

Volatile methylsiloxanes (VMS) are man-made, nonbiodegradable chemicals produced at a megaton-per-year scale, which leads to concern over their potential for environmental persistence, long-range transport, and bioaccumulation. We used directed evolution to engineer a variant of bacterial cytochrome P450BM3 to break silicon-carbon bonds in linear and cyclic VMS. To accomplish silicon-carbon bond cleavage, the enzyme catalyzes two tandem oxidations of a siloxane methyl group, which is followed by putative [1,2]-Brook rearrangement and hydrolysis. Discovery of this so-called siloxane oxidase opens possibilities for the eventual biodegradation of VMS.

Tryptophan extends the life of cytochrome P450

Powerfully oxidizing enzymes need protective mechanisms to prevent self-destruction. The flavocytochrome P450 BM3 from Priestia megaterium (P450BM3) is a self-sufficient monooxygenase that hydroxylates fatty acid substrates using O2 and NADPH as co-substrates. Hydroxylation of long-chain fatty acids (≥C14) is well coupled to O2 and NADPH consumption, but shorter chains (≤C12) are more poorly coupled. Hydroxylation of p-nitrophenoxydodecanoic acid by P450BM3 produces a spectrophotometrically detectable product wherein the coupling of NADPH consumption to product formation is just 10%. Moreover, the rate of NADPH consumption is 1.8 times that of O2 consumption, indicating that an oxidase uncoupling pathway is operative. Measurements of the total number of enzyme turnovers before inactivation (TTN) indicate that higher NADPH concentrations increase TTN. At lower NADPH levels, added ascorbate increases TTN, while a W96H mutation leads to a decrease. The W96 residue is about 7 Å from the P450BM3 heme and serves as a gateway residue in a tryptophan/tyrosine (W/Y) hole transport chain from the heme to a surface tyrosine residue. The data indicate that two oxidase pathways protect the enzyme from damage by intercepting the powerfully oxidizing enzyme intermediate (Compound I) and returning it to its resting state. At high NADPH concentrations, reducing equivalents from the flavoprotein are delivered to Compound I by the usual reductase pathway. When NADPH is not abundant, however, oxidizing equivalents from Compound I can traverse a W/Y chain, arriving at the enzyme surface where they are scavenged by reductants. Ubiquitous tryptophan/tyrosine chains in highly oxidizing enzymes likely perform similar protective functions.

Characterisation of the heme aqua-ligand coordination environment in an engineered peroxygenase cytochrome P450 variant

The cytochrome P450 enzymes (CYPs) are heme-thiolate monooxygenases that catalyse the insertion of an oxygen atom into the C-H bonds of organic molecules. In most CYPs, the activation of dioxygen by the heme is aided by an acid-alcohol pair of residues located in the I-helix of the enzyme. Mutation of the threonine residue of this acid-alcohol pair of CYP199A4, from the bacterium Rhodospeudomonas palustris HaA2, to a glutamate residue induces peroxygenase activity. In the X-ray crystal structures of this variant an interaction of the glutamate side chain and the distal aqua ligand of the heme was observed and this results in this ligand not being readily displaced in the peroxygenase mutant on the addition of substrate. Here we use a range of bulky hydrophobic and nitrogen donor containing ligands in an attempt to displace the distal aqua ligand of the T252E mutant of CYP199A4. Ligand binding was assessed by UV-visible absorbance spectroscopy, native mass spectrometry, electron paramagnetic resonance and X-ray crystallography. None of the ligands tested, even the nitrogen donor ligands which bind directly to the iron in the wild-type enzyme, resulted in displacement of the aqua ligand. Therefore, modification of the I-helix threonine residue to a glutamate residue results in a significant strengthening of the ferric distal aqua ligand. This ligand was not displaced using any of the ligands during this study and this provides a rationale as to why this mutant can shutdown the monooxygenase pathway of this enzyme and switch to peroxygenase activity.

Biosynthesis of High-Active Hemoproteins by the Efficient Heme-Supply Pichia Pastoris Chassis

Microbial synthesis of valuable hemoproteins has become a popular research topic, and Pichia pastoris is a versatile platform for the industrial production of recombinant proteins. However, the inadequate supply of heme limits the synthesis of high-active hemoproteins. Here a strategy for enhancing intracellular heme biosynthesis to improve the titers and functional activities of hemoproteins is reported. After selecting a suitable expressional strategy for globins, the efficient heme-supply P. pastoris chassis is established by removing the spatial segregation during heme biosynthesis, optimizing precursor synthesis, assembling rate-limiting enzymes using protein scaffolds, and inhibiting heme degradation. This robust chassis produces several highly active hemoproteins, including porcine myoglobin, soy hemoglobin, Vitreoscilla hemoglobin, and P450-BM3, which can be used in the development of artificial meat, high-cell-density fermentation, and whole-cell catalytic synthesis of high-value-added compounds. Furthermore, the engineered chassis strain has great potential for producing and applying other hemoproteins with high activities in various fields.

Tryptophan-96 in cytochrome P450 BM3 plays a key role in enzyme survival

Flavocytochrome P450 from Bacillus megaterium (P450BM3 ) is a natural fusion protein containing reductase and heme domains. In the presence of NADPH and dioxygen the enzyme catalyses the hydroxylation of long-chain fatty acids. Analysis of the P450BM3 structure reveals chains of closely spaced tryptophan and tyrosine residues that might serve as pathways for high-potential oxidizing equivalents to escape from the heme active site when substrate oxidation is not possible. Our investigations of the total number of enzyme turnovers before deactivation have revealed that replacement of selected tryptophan and tyrosine residues with redox inactive groups leads to a twofold reduction in enzyme survival time. Tryptophan-96 is critical for prolonging enzyme activity, suggesting a key protective role for this residue.

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

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

Genetic fusion of P450 BM3 and formate dehydrogenase towards self-sufficient biocatalysts with enhanced activity

Fusion of multiple enzymes to multifunctional constructs has been recognized as a viable strategy to improve enzymatic properties at various levels such as stability, activity and handling. In this study, the genes coding for cytochrome P450 BM3 from B. megaterium and formate dehydrogenase from Pseudomonas sp. were fused to enable both substrate oxidation catalyzed by P450 BM3 and continuous cofactor regeneration by formate dehydrogenase within one construct. The order of the genes in the fusion as well as the linkers that bridge the enzymes were varied. The resulting constructs were compared to individual enzymes regarding substrate conversion, stability and kinetic parameters to examine whether fusion led to any substantial improvements of enzymatic properties. Most noticeably, an activity increase of up to threefold was observed for the fusion constructs with various substrates which were partly attributed to the increased diflavin reductase activity of the P450 BM3. We suggest that P450 BM3 undergoes conformational changes upon fusion which resulted in altered properties, however, no NADPH channeling was detected for the fusion constructs.

Functional and protective hole hopping in metalloenzymes

Electrons can tunnel through proteins in microseconds with a modest release of free energy over distances in the 15 to 20 Å range. To span greater distances, or to move faster, multiple charge transfers (hops) are required. When one of the reactants is a strong oxidant, it is convenient to consider the movement of a positively charged "hole" in a direction opposite to that of the electron. Hole hopping along chains of tryptophan (Trp) and tyrosine (Tyr) residues is a critical function in several metalloenzymes that generate high-potential intermediates by reactions with O2 or H2O2, or by activation with visible light. Examination of the protein structural database revealed that Tyr/Trp chains are common protein structural elements, particularly among enzymes that react with O2 and H2O2. In many cases these chains may serve a protective role in metalloenzymes by deactivating high-potential reactive intermediates formed in uncoupled catalytic turnover.

Advances in enzymatic oxyfunctionalization of aliphatic compounds

Selective oxyfunctionalizations of aliphatic compounds are difficult chemical reactions, where enzymes can play an important role due to their stereo- and regio-selectivity and operation under mild reaction conditions. P450 monooxygenases are well-known biocatalysts that mediate oxyfunctionalization reactions in different living organisms (from bacteria to humans). Unspecific peroxygenases (UPOs), discovered in fungi, have arisen as "dream biocatalysts" of great biotechnological interest because they catalyze the oxyfunctionalization of aliphatic and aromatic compounds, avoiding the necessity of expensive cofactors and regeneration systems, and only depending on H₂O₂ for their catalysis. Here, we summarize recent advances in aliphatic oxyfunctionalization reactions by UPOs, as well as the molecular determinants of the enzyme structures responsible for their activities, emphasizing the differences found between well-known P450s and the novel fungal peroxygenases.

Assessing the Covalent Attachment and Energy Transfer Capabilities of Upconverting Phosphors With Cofactor Containing Bioactive Enzymes

Upconverting phosphors (UCPs) convert multiple low energy photons into higher energy emission via the process of photon upconversion and offer an attractive alternative to organic fluorophores for use as luminescent probes. Examples of biosensors utilizing the apparent energy transfer of UCPs and nanophosphors (UCNPs) with biomolecules have started to appear in the literature but very few exploit the covalent anchoring of the biomolecule to the surface of the UCP to improve the sensitivity of the systems. Here, we demonstrate a robust and versatile method for the covalent attachment of biomolecules to the surface of a variety of UCPs and UCNPs in which the UCPs were capped with functionalized silica in order to provide a surface to covalently conjugate biomolecules with surface-accessible cysteines. Variants of BM3Heme, cytochrome C, glucose oxidase, and glutathione reductase were then attached via maleimide-thiol coupling. BM3Heme, glucose oxidase, and glutathione reductase were shown to retain their activity when coupled to the UCPs potentially opening up opportunities for biosensing applications.

Crystals in Minutes: Instant On-Site Microcrystallisation of Various Flavours of the CYP102A1 (P450BM3) Haem Domain

Despite CYP102A1 (P450BM3) representing one of the most extensively researched metalloenzymes, crystallisation of its haem domain upon modification can be a challenge. Crystal structures are indispensable for the efficient structure-based design of P450BM3 as a biocatalyst. The abietane diterpenoid derivative N-abietoyl-l-tryptophan (AbiATrp) is an outstanding crystallisation accelerator for the wild-type P450BM3 haem domain, with visible crystals forming within 2 hours and diffracting to a near-atomic resolution of 1.22 Å. Using these crystals as seeds in a cross-microseeding approach, an assortment of P450BM3 haem domain crystal structures, containing previously uncrystallisable decoy molecules and diverse artificial metalloporphyrins binding various ligand molecules, as well as heavily tagged haem-domain variants, could be determined. Some of the structures reported herein could be used as models of different stages of the P450BM3 catalytic cycle.

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

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

Living with Oxygen

Work on the electronic structures of metal–oxo complexes began in Copenhagen over 50 years ago. This work led to the prediction that tetragonal multiply bonded transition metal–oxos would not be stable beyond the iron–ruthenium–osmium oxo wall in the periodic table and that triply bonded metal–oxos could not be protonated, even in the strongest Brønsted acids. In this theory, only double bonded metal–oxos could attract protons, with basicities being a function of the electron donating ability of ancillary ligands. Such correlations of electronic structure with reactivity have gained importance in recent years, most notably owing to the widespread recognition that high-valent iron–oxos are intermediates in biological reactions critical to life on Earth.

In this Account, we focus attention on the oxygenations of inert organic substrates by cytochromes P450, as these reactions involve multiply bonded iron–oxos. We emphasize that P450 iron–oxos are strong oxidants, so strong that they would destroy nearby amino acids if substrates are not oxygenated rapidly; it is our view that these high-valent iron–oxos are such dangerous reactive oxygen species that Nature surely found ways to disable them. Looking more deeply into this matter, mainly by examining many thousands of structures in the Protein Data Bank, we have found that P450s and other enzymes that require oxygen for function have chains of tyrosines and tryptophans that extend from active-site regions to protein surfaces. Tyrosines are near the heme active sites in bacterial P450s, whereas tryptophan is closest in most human enzymes. High-valent iron–oxo survival times taken from hole hopping maps range from a few nanoseconds to milliseconds, depending on the distance of the closest Trp or Tyr residue to the heme. In our proposed mechanism, multistep hole tunneling (hopping) through Tyr/Trp chains guides the damaging oxidizing hole to the protein surface, where it can be quenched by soluble protein or small molecule reductants. As the Earth’s oxygenic atmosphere is believed to have developed about 2.5 billion years ago, the increase in occurrence frequency of tyrosine and tryptophan since the last universal evolutionary ancestor may be in part a consequence of enzyme protective functions that developed to cope with the environmental toxin, O₂.

Homology Modeling and Molecular Dynamics Simulation Combined with X-ray Solution Scattering Defining Protein Structures of Thromboxane and Prostacyclin Synthases

A combination of molecular dynamics (MD) simulations and X-ray scattering (SAXS) has emerged as the approach of choice for studying protein structures and dynamics in solution. This approach has potential applications for membrane proteins that neither are soluble nor form crystals easily. We explore the water-coupled dynamic structures of thromboxane synthase (TXAS) and prostacyclin synthase (PGIS) from scanning HPLC-SAXS measurements combined with MD ensemble analyses. Both proteins are heme-containing enzymes in the cytochrome P450 family, known as prostaglandin H₂ (PGH₂) isomerase, with counter-functions in regulation of platelet aggregation. Currently, the X-ray crystallographic structures of PGIS are available, but those for TXAS are not. The use of homology modeling of the TXAS structure with ns-μs explicit water solvation MD simulations allows much more accurate estimation of the configuration space with loop motion and origin of the protein behaviors in solution. In contrast to the stability of the conserved PGIS structure in solution, the pronounced TXAS flexibility has been revealed to have unstructured loop regions in connection with the characteristic P450 structural elements. The MD-derived and experimental-solution SAXS results are in excellent agreement. The significant protein internal motions, whole-molecule structures, and potential problems with protein folding, crystallization, and functionality are examined.

Determinants of thermostability in the cytochrome P450 fold

Cytochromes P450 are found throughout the biosphere in a wide range of environments, serving a multitude of physiological functions. The ubiquity of the P450 fold suggests that it has been co-opted by evolution many times, and likely presents a useful compromise between structural stability and conformational flexibility. The diversity of substrates metabolized and reactions catalyzed by P450s makes them attractive starting materials for use as biocatalysts of commercially useful reactions. However, process conditions impose different requirements on enzymes to those in which they have evolved naturally. Most natural environments are relatively mild, and therefore most P450s have not been selected in Nature for the ability to withstand temperatures above ~40°C, yet industrial processes frequently require extended incubations at much higher temperatures. Thus, there has been considerable interest and effort invested in finding or engineering thermostable P450 systems. Numerous P450s have now been identified in thermophilic organisms and analysis of their structures provides information as to mechanisms by which the P450 fold can be stabilized. In addition, protein engineering, particularly by directed or artificial evolution, has revealed mutations that serve to stabilize particular mesophilic enzymes of interest. Here we review the current understanding of thermostability as it applies to the P450 fold, gleaned from the analysis of P450s characterized from thermophilic organisms and the parallel engineering of mesophilic forms for greater thermostability. We then present a perspective on how this information might be used to design stable P450 enzymes for industrial application. This article is part of a Special Issue entitled: Cytochrome P450 biodiversity and biotechnology, edited by Erika Plettner, Gianfranco Gilardi, Luet Wong, Vlada Urlacher, Jared Goldstone.

In vitro characterization of CYP102G4 from Streptomyces cattleya: A self-sufficient P450 naturally producing indigo

Self-sufficient CYP102As possess outstanding hydroxylating activity to fatty acids such as myristic acid. Other CYP102 subfamily members share substrate specificity of CYP102As, but, occasionally, unusual characteristics of its own subfamily have been found. In this study, only one self-sufficient cytochrome P450 from Streptomyces cattleya was renamed from CYP102A_scat to CYP102G4, purified and characterized. UV-Vis spectrometry pattern, FAD/FMN analysis, and protein sequence comparison among CYP102s have shown that CYP102 from Streptomyces cattleya belongs to CYP102G subfamily. It showed hydroxylation activity toward fatty acids generating ω-1, ω-2, and ω-3-hydroxyfatty acids, which is similar to the general substrate specificity of CYP102 family. Unexpectedly, however, expression of CYP102G4 showed indigo production in LB medium batch flask culture, and high catalytic activity (k_cat/K_m) for indole was measured as 6.14±0.10min^-1mM^-1. Besides indole, CYP102G4 was able to hydroxylate aromatic compounds such as flavone, benzophenone, and chloroindoles. Homology model has shown such ability to accept aromatic compounds is due to its bigger active site cavity. Unlike other CYP102s, CYP102G4 did not have biased cofactor dependency, which was possibly determined by difference in NAD(P)H binding residues (Ala984, Val990, and Tyr1064) compared to CYP102A1 (Arg966, Lys972 and Trp1046). Overall, a self-sufficient CYP within CYP102G subfamily was characterized using purified enzymes, which appears to possess unique properties such as an only prokaryotic CYP naturally producing indigo.

Crystallographic insights into a cobalt (III) sepulchrate based alternative cofactor system of P450 BM3 monooxygenase

P450 BM3 is a multi-domain heme-containing soluble bacterial monooxygenase. P450 BM3 and variants are known to oxidize structurally diverse substrates. Crystal structures of individual domains of P450 BM3 are available. However, the spatial organization of the full-length protein is unknown. In this study, crystal structures of the P450 BM3 M7 heme domain variant with and without cobalt (III) sepulchrate are reported. Cobalt (III) sepulchrate acts as an electron shuttle in an alternative cofactor system employing zinc dust as the electron source. The crystal structure shows a binding site for the mediator cobalt (III) sepulchrate at the entrance of the substrate access channel. The mediator occupies an unusual position which is far from the active site and distinct from the binding of the natural redox partner (FAD/NADPH binding domain).

Hole Hopping through Tryptophan in Cytochrome P450

Electron-transfer kinetics have been measured in four conjugates of cytochrome P450 with surface-bound Ru-photosensitizers. The conjugates are constructed with enzymes from Bacillus megaterium (CYP102A1) and Sulfolobus acidocaldarius (CYP119). A W96 residue lies in the path between Ru and the heme in CYP102A1, whereas H76 is present at the analogous location in CYP119. Two additional conjugates have been prepared with (CYP102A1)W96H and (CYP119)H76W mutant enzymes. Heme oxidation by photochemically generated Ru³⁺ leads to P450 compound II formation when a tryptophan residue is in the path between Ru and the heme; no heme oxidation is observed when histidine occupies this position. The data indicate that heme oxidation proceeds via two-step tunneling through a tryptophan radical intermediate. In contrast, heme reduction by photochemically generated Ru⁺ proceeds in a single electron tunneling step with closely similar rate constants for all four conjugates.

Influence of haem environment on the catalytic properties of the tetrathionate reductase TsdA from Campylobacter jejuni

In the present study, we provide a detailed analysis of the catalytic properties of the bifunctional thiosulfate dehydrogenases/tetrathionate reductases (TsdA) of the human food-borne pathogen Campylobacter jejuni. Structural differences in the immediate environment of Haem 2 were shown to influence the reaction directionality.

Bifunctional dihaem cytochrome c thiosulfate dehydrogenases/tetrathionate reductases (TsdA) exhibit different catalytic properties depending on the source organism. In the human food-borne intestinal pathogen Campylobacter jejuni, TsdA functions as a tetrathionate reductase enabling respiration with tetrathionate as an alternative electron acceptor. In the present study, evidence is provided that Cys¹³⁸ and Met²⁵⁵ serve as the sixth ligands of Haem 1 and Haem 2 respectively, in the oxidized CjTsdA wt protein. Replacement of Cys¹³⁸ resulted in a virtually inactive enzyme, confirming Haem 1 as the active site haem. Significantly, TsdA variants carrying amino acid exchanges in the vicinity of the electron-transferring Haem 2 (Met²⁵⁵, Asn²⁵⁴ and Lys²⁵²) exhibited markedly altered catalytic properties of the enzyme, showing these residues play a key role in the physiological function of TsdA. The growth phenotypes and tetrathionate reductase activities of a series of ΔtsdA/*tsdA complementation strains constructed in the original host C. jejuni 81116, showed that in vivo, the TsdA variants exhibited the same catalytic properties as the pure, recombinantly produced enzymes. However, variants that catalysed tetrathionate reduction more effectively than the wild-type enzyme did not allow better growth.

An Evolved Orthogonal Enzyme/Cofactor Pair

We introduce a strategy that expands the functionality of hemoproteins through orthogonal enzyme/heme pairs. By exploiting the ability of a natural heme transport protein, ChuA, to promiscuously import heme derivatives, we have evolved a cytochrome P450 (P450BM3) that selectively incorporates a nonproteinogenic cofactor, iron deuteroporphyrin IX (Fe-DPIX), even in the presence of endogenous heme. Crystal structures show that selectivity gains are due to mutations that introduce steric clash with the heme vinyl groups while providing a complementary binding surface for the smaller Fe-DPIX cofactor. Furthermore, the evolved orthogonal enzyme/cofactor pair is active in non-natural carbenoid-mediated olefin cyclopropanation. This methodology for the generation of orthogonal enzyme/cofactor pairs promises to expand cofactor diversity in artificial metalloenzymes.

Electron flow through biological molecules: does hole hopping protect proteins from oxidative damage?

Biological electron transfers often occur between metal-containing cofactors that are separated by very large molecular distances. Employing photosensitizer-modified iron and copper proteins, we have shown that single-step electron tunneling can occur on nanosecond to microsecond timescales at distances between 15 and 20 Å. We also have shown that charge transport can occur over even longer distances by hole hopping (multistep tunneling) through intervening tyrosines and tryptophans. In this perspective, we advance the hypothesis that such hole hopping through Tyr/Trp chains could protect oxygenase, dioxygenase, and peroxidase enzymes from oxidative damage. In support of this view, by examining the structures of P450 (CYP102A) and 2OG-Fe (TauD) enzymes, we have identified candidate Tyr/Trp chains that could transfer holes from uncoupled high-potential intermediates to reductants in contact with protein surface sites.

Mutating a Highly Conserved Residue in Diverse Cytochrome P450s Facilitates Diastereoselective Olefin Cyclopropanation

Cytochrome P450s and other heme-containing proteins have recently been shown to have promiscuous activity for the cyclopropanation of olefins using diazoacetate reagents. Despite the progress made thus far, engineering selective catalysts for all possible stereoisomers for the cyclopropanation reaction remains a considerable challenge. Previous investigations of a model P450 (P450BM3 ) revealed that mutation of a conserved active site threonine (Thr268) to alanine transformed the enzyme into a highly active and selective cyclopropanation catalyst. By incorporating this mutation into a diverse panel of P450 scaffolds, we were able to quickly identify enantioselective catalysts for all possible diastereomers in the model reaction of styrene with ethyl diazoacetate. Some alanine variants exhibited selectivities that were markedly different from the wild-type enzyme, with a few possessing moderate to high diastereoselectivity and enantioselectivities up to 97 % for synthetically challenging cis-cyclopropane diastereomers.

Insights into regioselective metabolism of mefenamic acid by cytochrome P450 BM3 mutants through crystallography, docking, molecular dynamics, and free energy calculations

Cytochrome P450 BM3 (CYP102A1) mutant M11 is able to metabolize a wide range of drugs and drug-like compounds. Among these, M11 was recently found to be able to catalyze formation of human metabolites of mefenamic acid and other nonsteroidal anti-inflammatory drugs (NSAIDs). Interestingly, single active-site mutations such as V87I were reported to invert regioselectivity in NSAID hydroxylation. In this work, we combine crystallography and molecular simulation to study the effect of single mutations on binding and regioselective metabolism of mefenamic acid by M11 mutants. The heme domain of the protein mutant M11 was expressed, purified, and crystallized, and its X-ray structure was used as template for modeling. A multistep approach was used that combines molecular docking, molecular dynamics (MD) simulation, and binding free-energy calculations to address protein flexibility. In this way, preferred binding modes that are consistent with oxidation at the experimentally observed sites of metabolism (SOMs) were identified. Whereas docking could not be used to retrospectively predict experimental trends in regioselectivity, we were able to rank binding modes in line with the preferred SOMs of mefenamic acid by M11 and its mutants by including protein flexibility and dynamics in free-energy computation. In addition, we could obtain structural insights into the change in regioselectivity of mefenamic acid hydroxylation due to single active-site mutations. Our findings confirm that use of MD and binding free-energy calculation is useful for studying biocatalysis in those cases in which enzyme binding is a critical event in determining the selective metabolism of a substrate.

Could tyrosine and tryptophan serve multiple roles in biological redox processes?

Single-step electron tunnelling reactions can transport charges over distances of 15-20 Åin proteins. Longer-range transfer requires multi-step tunnelling processes along redox chains, often referred to as hopping. Long-range hopping via oxidized radicals of tryptophan and tyrosine, which has been identified in several natural enzymes, has been demonstrated in artificial constructs of the blue copper protein azurin. Tryptophan and tyrosine serve as hopping way stations in high-potential charge transport processes. It may be no coincidence that these two residues occur with greater-than-average frequency in O(2)- and H(2)O(2)-reactive enzymes. We suggest that appropriately placed tyrosine and/or tryptophan residues prevent damage from high-potential reactive intermediates by reduction followed by transfer of the oxidizing equivalent to less harmful sites or out of the protein altogether.

Thiosulfate dehydrogenase (TsdA) from Allochromatium vinosum: structural and functional insights into thiosulfate oxidation

Although the oxidative condensation of two thiosulfate anions to tetrathionate constitutes a well documented and significant part of the natural sulfur cycle, little is known about the enzymes catalyzing this reaction. In the purple sulfur bacterium Allochromatium vinosum, the reaction is catalyzed by the periplasmic diheme c-type cytochrome thiosulfate dehydrogenase (TsdA). Here, we report the crystal structure of the "as isolated" form of A. vinosum TsdA to 1.98 Å resolution and those of several redox states of the enzyme to different resolutions. The protein contains two typical class I c-type cytochrome domains wrapped around two hemes axially coordinated by His(53)/Cys(96) and His(164)/Lys(208). These domains are very similar, suggesting a gene duplication event during evolution. A ligand switch from Lys(208) to Met(209) is observed upon reduction of the enzyme. Cys(96) is an essential residue for catalysis, with the specific activity of the enzyme being completely abolished in several TsdA-Cys(96) variants. TsdA-K208N, K208G, and M209G variants were catalytically active in thiosulfate oxidation as well as in tetrathionate reduction, pointing to heme 2 as the electron exit point. In this study, we provide spectroscopic and structural evidence that the TsdA reaction cycle involves the transient presence of heme 1 in the high-spin state caused by movement of the Sγ atom of Cys(96) out of the iron coordination sphere. Based on the presented data, we draw important conclusions about the enzyme and propose a possible reaction mechanism for TsdA.

Fluorescence quenching of (dimethylamino)naphthalene dyes Badan and Prodan by tryptophan in cytochromes P450 and micelles

Fluorescence of 2-(N,N-dimethylamino)-6-propionylnaphthalene dyes Badan and Prodan is quenched by tryptophan in Brij 58 micelles as well as in two cytochrome P450 proteins (CYP102, CYP119) with Badan covalently attached to a cysteine residue. Formation of nonemissive complexes between a dye molecule and tryptophan accounts for about 76% of the fluorescence intensity quenching in micelles, the rest is due to diffusive encounters. In the absence of tryptophan, fluorescence of Badan-labeled cytochromes decays with triexponential kinetics characterized by lifetimes of about 100 ps, 700–800 ps, and 3 ns. Site mutation of a histidine residue in the vicinity of the Badan label by tryptophan results in shortening of all three decay lifetimes. The relative amplitude of the fastest component increases at the expense of the two slower ones. The average quenching rate constants are 4.5 × 10⁸ s^–1 (CYP102) and 3.7 × 10⁸ s^–1 (CYP119), at 288 K. Cyclic voltammetry of Prodan in MeCN shows a reversible reduction peak at −1.85 V vs NHE that becomes chemically irreversible and shifts positively upon addition of water. A quasireversible reduction at −0.88 V was observed in an aqueous buffer (pH 7.3). The excited-state reduction potential of Prodan (and Badan) is estimated to vary from about +0.6 V (vs NHE) in polar aprotic media (MeCN) to approximately +1.6 V in water. Tryptophan quenching of Badan/Prodan fluorescence in CYPs and Brij 58 micelles is exergonic by ≤0.5 V and involves tryptophan oxidation by excited Badan/Prodan, coupled with a fast reaction between the reduced dye and water. Photoreduction is a new quenching mechanism for 2-(N,N-dimethylamino)-6-propionylnaphthalene dyes that are often used as solvatochromic polarity probes, FRET donors and acceptors, as well as reporters of solvation dynamics.

A facile and versatile methodology for cysteine specific labeling of proteins with octahedral polypyridyl d⁶ metal complexes

We have synthesized and characterized four octahedral polypyridyl d(6) metal complexes bearing the 5,6-epoxy-5,6-dihydro-[1,10]phenanthroline ligand (L1) as cysteine specific labeling reagents. The proposed synthetic pathways allow the preparation of the metal complexes containing Re(I), Ru(II), Os(II) and Ir(III) while preserving the epoxide functionality. The complexes were characterized by (1)H and (13)C NMR, mass spectrometry, UV-visible and luminescence spectroscopies as well as cyclic voltammetry. As proof of concept, a set of non-native single cysteine P450 BM3 heme domain mutants previously developed in our laboratory was used to study the labeling reaction. We demonstrate that the proposed labels can selectively react, often in high yield, with cysteine residues of the protein via the nucleophilic thiol ring opening of the epoxide moiety. In addition, under basic conditions, subsequent loss of a water molecule led to the aromatization of the phenanthroline ring on the protein-bound label compounds, as observed by mass spectrometry and luminescence measurements.

Human P450-like oxidation of diverse proton pump inhibitor drugs by ‘gatekeeper’ mutants of flavocytochrome P450 BM3

Structurally destabilizing and conformation altering ‘gatekeeper’ mutants of the biotechnologically important flavocytochrome P450 BM3 mono-oxygenase were shown to catalyse oxidation of several proton pump inhibitor drugs, and to generate products identical to those formed by major human drug metabolizing P450s.

Strength of axial water ligation in substrate-free cytochrome P450s is isoform dependent

The heme-containing cytochrome P450s exhibit isoform-dependent ferric spin equilibria in the resting state and differential substrate-dependent spin equilibria. The basis for these differences is not well understood. Here, magnetic circular dichroism (MCD) reveals significant differences in the resting low spin ligand field of CYPs 3A4, 2E1, 2C9, 125A1, and 51B1, which indicates differences in the strength of axial water ligation to the heme. The near-infrared bands that specifically correspond to charge-transfer porphyrin-to-metal transitions span a range of energies of nearly 2 kcal/mol. In addition, the experimentally determined MCD bands are not entirely in agreement with the expected MCD energies calculated from electron paramagnetic resonance parameters, thus emphasizing the need for the experimental data. MCD marker bands of the high spin heme between 500 and 680 nm were also measured and suggest only a narrow range of energies for this ensemble of high spin Cys(S(-)) → Fe(3+) transitions among these isoforms. The differences in axial ligand energies between CYP isoforms of the low spin states likely contribute to the energetics of substrate-dependent spin state perturbation. However, these ligand field energies do not correlate with the fraction of high spin vs low spin in the resting state enzyme, suggestive of differences in water access to the heme or isoform-dependent differences in the substrate-free high spin states as well.

Dynamics and flexibility of human aromatase probed by FTIR and time resolved fluorescence spectroscopy

Human aromatase (CYP19A1) is a steroidogenic cytochrome P450 converting androgens into estrogens. No ligand-free crystal structure of the enzyme is available to date. The crystal structure in complex with the substrate androstenedione and the steroidal inhibitor exemestane shows a very compact conformation of the enzyme, leaving unanswered questions on the conformational changes that must occur to allow access of the ligand to the active site. As H/D exchange kinetics followed by FTIR spectroscopy can provide information on the conformational changes in proteins where solvent accessibility is affected, here the amide I region was used to measure the exchange rates of the different elements of the secondary structure for aromatase in the ligand-free form and in the presence of the substrate androstenedione and the inhibitor anastrozole. Biphasic exponential functions were found to fit the H/D exchange data collected as a function of time. Two exchange rates were assigned to two populations of protons present in different flexible regions of the protein. The addition of the substrate androstenedione and the inhibitor anastrozole lowers the H/D exchange rates of the α-helices of the enzyme when compared to the ligand-free form. Furthermore, the presence of the inhibitor anastrozole lowers exchange rate constant (k₁) for β-sheets from 0.22±0.06 min⁻¹ for the inhibitor-bound enzyme to 0.12±0.02 min⁻¹ for the free protein. Dynamics effects localised in helix F were studied by time resolved fluorescence. The data demonstrate that the fluorescence lifetime component associated to Trp224 emission undergoes a shift toward longer lifetimes (from ≈5.0 to ≈5.5 ns) when the substrate or the inhibitor are present, suggesting slower dynamics in the presence of ligands. Together the results are consistent with different degrees of flexibility of the access channel and therefore different conformations adopted by the enzyme in the free, substrate- and inhibitor-bound forms.

Structural evidence: a single charged residue affects substrate binding in cytochrome P450 BM-3

Cytochrome P450 BM-3 is a bacterial enzyme with sequence similarity to mammalian P450s that catalyzes the hydroxylation of fatty acids with high efficiency. Enzyme-substrate binding and dynamics has been an important topic of study for cytochromes P450 because most of the crystal structures of substrate-bound structures show the complex in an inactive state. We have determined a new crystal structure for cytochrome P450 BM-3 in complex with N-palmitoylglycine (NPG), which unexpectedly showed a direct bidentate ion pair between NPG and arginine 47 (R47). We further explored the role of R47, the only charged residue in the binding pocket in cytochrome P450 BM-3, through mutagenesis and crystallographic studies. The mutations of R47 to glutamine (R47Q), glutamic acid (R47E), and lysine (R47K) were designed to investigate the role of its charge in binding and catalysis. The oppositely charged R47E mutation had the greatest effect on activity and binding. The crystal structure of R47E BMP shows that the glutamic acid side chain is blocking the entrance to the binding pocket, accounting for NPG's low binding affinity and charge repulsion. For R47Q and R47K BM-3, the mutations caused only a slight change in kcat and a large change in Km and Kd, which suggests that R47 mostly is involved in binding and that our crystal structure, 4KPA , represents an initial binding step in the P450 cycle.

Key mutations alter the cytochrome P450 BM3 conformational landscape and remove inherent substrate bias

Cytochrome P450 monooxygenases (P450s) have enormous potential in the production of oxychemicals, due to their unparalleled regio- and stereoselectivity. The Bacillus megaterium P450 BM3 enzyme is a key model system, with several mutants (many distant from the active site) reported to alter substrate selectivity. It has the highest reported monooxygenase activity of the P450 enzymes, and this catalytic efficiency has inspired protein engineering to enable its exploitation for biotechnologically relevant oxidations with structurally diverse substrates. However, a structural rationale is lacking to explain how these mutations have such effects in the absence of direct change to the active site architecture. Here, we provide the first crystal structures of BM3 mutants in complex with a human drug substrate, the proton pump inhibitor omeprazole. Supported by solution data, these structures reveal how mutation alters the conformational landscape and decreases the free energy barrier for transition to the substrate-bound state. Our data point to the importance of such "gatekeeper" mutations in enabling major changes in substrate recognition. We further demonstrate that these mutants catalyze the same 5-hydroxylation reaction as performed by human CYP2C19, the major human omeprazole-metabolizing P450 enzyme.

Structure-guided directed evolution of highly selective p450-based magnetic resonance imaging sensors for dopamine and serotonin

New tools that allow dynamic visualization of molecular neural events are important for studying the basis of brain activity and disease. Sensors that permit ligand-sensitive magnetic resonance imaging (MRI) are useful reagents due to the noninvasive nature and good temporal and spatial resolution of MR methods. Paramagnetic metalloproteins can be effective MRI sensors due to the selectivity imparted by the protein active site and the ability to tune protein properties using techniques such as directed evolution. Here, we show that structure-guided directed evolution of the active site of the cytochrome P450-BM3 heme domain produces highly selective MRI probes with submicromolar affinities for small molecules. We report a new, high-affinity dopamine sensor as well as the first MRI reporter for serotonin, with which we demonstrate quantification of neurotransmitter release in vitro. We also present a detailed structural analysis of evolved cytochrome P450-BM3 heme domain lineages to systematically dissect the molecular basis of neurotransmitter binding affinity, selectivity, and enhanced MRI contrast activity in these engineered proteins.

Involvement of a natural fusion of a cytochrome P450 and a hydrolase in mycophenolic acid biosynthesis

Mycophenolic acid (MPA) is a fungal secondary metabolite and the active component in several immunosuppressive pharmaceuticals. The gene cluster coding for the MPA biosynthetic pathway has recently been discovered in Penicillium brevicompactum, demonstrating that the first step is catalyzed by MpaC, a polyketide synthase producing 5-methylorsellinic acid (5-MOA). However, the biochemical role of the enzymes encoded by the remaining genes in the MPA gene cluster is still unknown. Based on bioinformatic analysis of the MPA gene cluster, we hypothesized that the step following 5-MOA production in the pathway is carried out by a natural fusion enzyme MpaDE, consisting of a cytochrome P450 (MpaD) in the N-terminal region and a hydrolase (MpaE) in the C-terminal region. We verified that the fusion gene is indeed expressed in P. brevicompactum by obtaining full-length sequence of the mpaDE cDNA prepared from the extracted RNA. Heterologous coexpression of mpaC and the fusion gene mpaDE in the MPA-nonproducer Aspergillus nidulans resulted in the production of 5,7-dihydroxy-4-methylphthalide (DHMP), the second intermediate in MPA biosynthesis. Analysis of the strain coexpressing mpaC and the mpaD part of mpaDE shows that the P450 catalyzes hydroxylation of 5-MOA to 4,6-dihydroxy-2-(hydroxymethyl)-3-methylbenzoic acid (DHMB). DHMB is then converted to DHMP, and our results suggest that the hydrolase domain aids this second step by acting as a lactone synthase that catalyzes the ring closure. Overall, the chimeric enzyme MpaDE provides insight into the genetic organization of the MPA biosynthesis pathway.

The crystal structure of the FAD/NADPH-binding domain of flavocytochrome P450 BM3

We report the crystal structure of the FAD/NADPH-binding domain (FAD domain) of the biotechnologically important Bacillus megaterium flavocytochrome P450 BM3, the last domain of the enzyme to be structurally resolved. The structure was solved in both the absence and presence of the ligand NADP(+), identifying important protein interactions with the NADPH 2'-phosphate that helps to dictate specificity for NADPH over NADH, and involving residues Tyr974, Arg966, Lys972 and Ser965. The Trp1046 side chain shields the FAD isoalloxazine ring from NADPH, and motion of this residue is required to enable NADPH-dependent FAD reduction. Multiple binding interactions stabilize the FAD cofactor, including aromatic stacking with the adenine group from the side chains of Tyr860 and Trp854, and several interactions with FAD pyrophosphate oxygens, including bonding to tyrosines 828, 829 and 860. Mutagenesis of C773 and C999 to alanine was required for successful crystallization, with C773A predicted to disfavour intramolecular and intermolecular disulfide bonding. Multiangle laser light scattering analysis showed wild-type FAD domain to be near-exclusively dimeric, with dimer disruption achieved on treatment with the reducing agent dithiothreitol. By contrast, light scattering showed that the C773A/C999A FAD domain was monomeric. The C773A/C999A FAD domain structure confirms that Ala773 is surface exposed and in close proximity to Cys810, with this region of the enzyme's connecting domain (that links the FAD domain to the FMN-binding domain in P450 BM3) located at a crystal contact interface between FAD domains. The FAD domain crystal structure enables molecular modelling of its interactions with its cognate FMN (flavodoxin-like) domain within the BM3 reductase module.

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.

Identification of mutant Asp251Gly/Gln307His of cytochrome P450 BM3 for the generation of metabolites of diclofenac, ibuprofen and tolbutamide

The soluble, catalytically self-sufficient cytochrome P450 BM3 from Bacillus megaterium is a good candidate as biocatalyst for the synthesis of drug metabolites. To this end, error-prone polymerase chain reaction (PCR) was used to generate a library of P450 BM3 mutants with novel activities toward drugs. The double mutant Asp251Gly/Gln307His (A2) with activities towards diclofenac, ibuprofen and tolbutamide was identified by screening with the alkali method. This is based on the detection of NADPH oxidation during enzymatic turnover on whole Escherichia coli cells heterologously expressing the P450 BM3 mutants in the presence of the target substrates. The three drugs screened are marker substrates of human liver cytochromes P450 belonging to the 2C subfamily. Interestingly the mutations Asp251Gly/Gln307His are located on the protein surface and they are not directly involved in substrate binding and turnover. Dissociation constants and K(M) values of mutant A2 for diclofenac, ibuprofen and tolbutamide are in the micromolar range. Catalysis leads to hydroxylations in specific positions, producing 4'-hydroxydiclofenac, 2-hydroxyibuprofen and 4-hydroxytolbutamide, respectively.

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.

Finding homes for orphan cytochrome P450s: CYP4V2 and CYP4F22 in disease states

The cytochrome P450 (CYP) 4 family of enzymes contains several recently identified membersthat are referred to as “orphan P450s” because their endogenous substrates are unknown.Human CYP4V2 and CYP4F22 are two such orphan P450s that are strongly linked to ocular andskin disease, respectively. Genetic analyses have identified a wide spectrum of mutations in the CYP4V2gene from patients suffering from Bietti’s crystalline corneoretinal dystrophy, and mutations in theCYP4F22 gene have been linked to lamellar ichthyosis. The strong gene–disease associations provideunique opportunities for elucidating the substrate specificity of these orphan P450s and unraveling thebiochemical pathways that may be impacted in patients with CYP4V2 and CYP4F22 functional deficits.

Conformational plasticity and structure/function relationships in cytochromes P450

The cytochrome P450s are a superfamily of enzymes that are found in all kingdoms of living organisms, and typically catalyze the oxidative addition of atomic oxygen to an unactivated C-C or C-H bond. Over 8000 nonredundant sequences of putative and confirmed P450 enzymes have been identified, but three-dimensional structures have been determined for only a small fraction of these. While all P450 enzymes for which structures have been determined share a common global fold, the flexibility and modularity of structure around the active site account for the ability of P450 enzymes to accommodate a vast number of structurally dissimilar substrates and support a wide range of selective oxidations. In this review, known P450 structures are compared, and some structural criteria for prediction of substrate selectivity and reaction type are suggested. The importance of dynamic processes such as redox-dependent and effector-induced conformational changes in determining catalytic competence and regio- and stereoselectivity is discussed, and noncrystallographic methods for characterizing P450 structures and dynamics, in particular, mass spectrometry and nuclear magnetic resonance spectroscopy are reviewed.