Structure Based Discovery of Inhibitors of CYP125 and CYP142 from Mycobacterium tuberculosis

Mycobacterium tuberculosis (Mtb) was responsible for approximately 1.6 million deaths in 2021. With the emergence of extensive drug resistance, novel therapeutic agents are urgently needed, and continued drug discovery efforts required. Host-derived lipids such as cholesterol not only support Mtb growth, but are also suspected to function in immunomodulation, with links to persistence and immune evasion. Mtb cytochrome P450 (CYP) enzymes facilitate key steps in lipid catabolism and thus present potential targets for inhibition. Here we present a series of compounds based on an ethyl 5-(pyridin-4-yl)-1H-indole-2-carboxylate pharmacophore which bind strongly to both Mtb cholesterol oxidases CYP125 and CYP142. Using a structure-guided approach, combined with biophysical characterization, compounds with micromolar range in-cell activity against clinically relevant drug-resistant isolates were obtained. These will incite further development of much-needed additional treatment options and provide routes to probe the role of CYP125 and CYP142 in Mtb pathogenesis.

Computation-Aided Engineering of Cytochrome P450 for the Production of Pravastatin

CYP105AS1 is a cytochrome P450 from Amycolatopsis orientalis that catalyzes monooxygenation of compactin to 6-epi-pravastatin. For fermentative production of the cholesterol-lowering drug pravastatin, the stereoselectivity of the enzyme needs to be inverted, which has been partially achieved by error-prone PCR mutagenesis and screening. In the current study, we report further optimization of the stereoselectivity by a computationally aided approach. Using the CoupledMoves protocol of Rosetta, a virtual library of mutants was designed to bind compactin in a pro-pravastatin orientation. By examining the frequency of occurrence of beneficial substitutions and rational inspection of their interactions, a small set of eight mutants was predicted to show the desired selectivity and these variants were tested experimentally. The best CYP105AS1 variant gave >99% stereoselective hydroxylation of compactin to pravastatin, with complete elimination of the unwanted 6-epi-pravastatin diastereomer. The enzyme-substrate complexes were also examined by ultrashort molecular dynamics simulations of 50 × 100 ps and 5 × 22 ns, which revealed that the frequency of occurrence of near-attack conformations agreed with the experimentally observed stereoselectivity. These results show that a combination of computational methods and rational inspection could improve CYP105AS1 stereoselectivity beyond what was obtained by directed evolution. Moreover, the work lays out a general in silico framework for specificity engineering of enzymes of known structure.

A new strategy for hit generation: Novel in cellulo active inhibitors of CYP121A1 from Mycobacterium tuberculosis via a combined X-ray crystallographic and phenotypic screening approach (XP screen)

There is a pressing need for new drugs against tuberculosis (TB) to combat the growing resistance to current antituberculars. Herein a novel strategy is described for hit generation against promising TB targets involving X-ray crystallographic screening in combination with phenotypic screening. This combined approach (XP Screen) affords both a validation of target engagement as well as determination of in cellulo activity. The utility of this method is illustrated by way of an XP Screen against CYP121A1, a cytochrome P450 enzyme from Mycobacterium tuberculosis (Mtb) championed as a validated drug discovery target. A focused screening set was synthesized and tested by such means, with several members of the set showing promising activity against Mtb strain H37Rv. One compound was observed as an X-ray hit against CYP121A1 and showed improved activity against Mtb strain H37Rv under multiple assay conditions (pan-assay activity). Data obtained during X-ray crystallographic screening were utilized in a structure-based campaign to design a limited number of analogues (less than twenty), many of which also showed pan-assay activity against Mtb strain H37Rv. These included the benzo[b][1,4]oxazine derivative (MIC₉₀ 6.25 μM), a novel hit compound suitable as a starting point for a more involved hit to lead candidate medicinal chemistry campaign.

•

CYP121 from M.tuberculosis has been previously shown to be a crucial target for the survival of the mycobacteria.

•

Strategies previously employed have identified high affinity inhibitors however these have lacked activity on M.tuberculosis.

•

The strategy reported here uses a combination of X-ray crystallography and phenotypic screening (XP Screen) to identify compounds.

•

The XP screen approach identified a number of compounds which show good affinity (up to 3.2 μM) and MIC against M.tuberculosis (up to 6.25 μM).

Catalytic Mechanism of Aromatic Nitration by Cytochrome P450 TxtE: Involvement of a Ferric-Peroxynitrite Intermediate

The cytochromes P450 are heme-dependent enzymes that catalyze many vital reaction processes in the human body related to biodegradation and biosynthesis. They typically act as mono-oxygenases; however, the recently discovered P450 subfamily TxtE utilizes O₂ and NO to nitrate aromatic substrates such as L-tryptophan. A direct and selective aromatic nitration reaction may be useful in biotechnology for the synthesis of drugs or small molecules. Details of the catalytic mechanism are unknown, and it has been suggested that the reaction should proceed through either an iron(III)-superoxo or an iron(II)-nitrosyl intermediate. To resolve this controversy, we used stopped-flow kinetics to provide evidence for a catalytic cycle where dioxygen binds prior to NO to generate an active iron(III)-peroxynitrite species that is able to nitrate l-Trp efficiently. We show that the rate of binding of O₂ is faster than that of NO and also leads to l-Trp nitration, while little evidence of product formation is observed from the iron(II)-nitrosyl complex. To support the experimental studies, we performed density functional theory studies on large active site cluster models. The studies suggest a mechanism involving an iron(III)-peroxynitrite that splits homolytically to form an iron(IV)-oxo heme (Compound II) and a free NO₂ radical via a small free energy of activation. The latter activates the substrate on the aromatic ring, while compound II picks up the ipso-hydrogen to form the product. The calculations give small reaction barriers for most steps in the catalytic cycle and, therefore, predict fast product formation from the iron(III)-peroxynitrite complex. These findings provide the first detailed insight into the mechanism of nitration by a member of the TxtE subfamily and highlight how the enzyme facilitates this novel reaction chemistry.

Characterization of the structure and interactions of P450 BM3 using hybrid mass spectrometry approaches

The cytochrome P450 monooxygenase P450 BM3 (BM3) is a biotechnologically important and versatile enzyme capable of producing important compounds such as the medical drugs pravastatin and artemether, and the steroid hormone testosterone. BM3 is a natural fusion enzyme comprising two major domains: a cytochrome P450 (heme-binding) catalytic domain and a NADPH-cytochrome P450 reductase (CPR) domain containing FAD and FMN cofactors in distinct domains of the CPR. A crystal structure of full-length BM3 enzyme is not available in its monomeric or catalytically active dimeric state. In this study, we provide detailed insights into the protein-protein interactions that occur between domains in the BM3 enzyme and characterize molecular interactions within the BM3 dimer by using several hybrid mass spectrometry (MS) techniques, namely native ion mobility MS (IM-MS), collision-induced unfolding (CIU), and hydrogen-deuterium exchange MS (HDX-MS). These methods enable us to probe the structure, stoichiometry, and domain interactions in the ∼240 kDa BM3 dimeric complex. We obtained high-sequence coverage (88–99%) in the HDX-MS experiments for full-length BM3 and its component domains in both the ligand-free and ligand-bound states. We identified important protein interaction sites, in addition to sites corresponding to heme-CPR domain interactions at the dimeric interface. These findings bring us closer to understanding the structure and catalytic mechanism of P450 BM3.

Clobetasol Propionate Is a Heme-Mediated Selective Inhibitor of Human Cytochrome P450 3A5

The human cytochrome P450 (CYP) enzymes CYP3A4 and CYP3A5 metabolize most drugs and have high similarities in their structure and substrate preference. Whereas CYP3A4 is predominantly expressed in the liver, CYP3A5 is upregulated in cancer, contributing to drug resistance. Selective inhibitors of CYP3A5 are, therefore, critical to validating it as a therapeutic target. Here we report clobetasol propionate (clobetasol) as a potent and selective CYP3A5 inhibitor identified by high-throughput screening using enzymatic and cell-based assays. Molecular dynamics simulations suggest a close proximity of clobetasol to the heme in CYP3A5 but not in CYP3A4. UV-visible spectroscopy and electron paramagnetic resonance analyses confirmed the formation of an inhibitory type I heme-clobetasol complex in CYP3A5 but not in CYP3A4, thus explaining the CYP3A5 selectivity of clobetasol. Our results provide a structural basis for selective CYP3A5 inhibition, along with mechanistic insights, and highlight clobetasol as an important chemical tool for target validation.

MhuD from Mycobacterium tuberculosis: Probing a Dual Role in Heme Storage and Degradation

The Mycobacterium tuberculosis (Mtb) heme oxygenase MhuD liberates free iron by degrading heme to the linear tetrapyrrole mycobilin. The MhuD dimer binds up to two hemes within the active site of each monomer. Binding the first solvent-exposed heme allows heme degradation and releases free iron. Binding a second heme renders MhuD inactive, allowing heme storage. Native-mass spectrometry revealed little difference in binding affinity between solvent-exposed and solvent-protected hemes. Hence, diheme-MhuD is formed even when a large proportion of the MhuD population is in the apo form. Apomyoglobin heme transfer assays showed MhuD-diheme dissociation is far slower than monoheme dissociation at ∼0.12 min^-1 and ∼0.25 s^-1, respectively, indicating that MhuD has a strong affinity for diheme. MhuD has not evolved to preferentially occupy the monoheme form and, through formation of a diheme complex, it functions as part of a larger network to tightly regulate both heme and iron levels in Mtb.

Structure-Activity Relationships of cyclo(l-Tyrosyl-l-tyrosine) Derivatives Binding to Mycobacterium tuberculosis CYP121: Iodinated Analogues Promote Shift to High-Spin Adduct

A series of analogues of cyclo(l-tyrosyl-l-tyrosine), the substrate of the Mycobacterium tuberculosis enzyme CYP121, have been synthesized and analyzed by UV-vis and electron paramagnetic resonance spectroscopy and by X-ray crystallography. The introduction of iodine substituents onto cyclo(l-tyrosyl-l-tyrosine) results in sub-μM binding affinity for the CYP121 enzyme and a complete shift to the high-spin state of the heme Fe^III. The introduction of halogens that are able to interact with heme groups is thus a feasible approach to the development of next-generation, tight binding inhibitors of the CYP121 enzyme, in the search for novel antitubercular compounds.

Design and Synthesis of Imidazole and Triazole Pyrazoles as Mycobacterium Tuberculosis CYP121A1 Inhibitors

The emergence of untreatable drug-resistant strains of Mycobacterium tuberculosis is a major public health problem worldwide, and the identification of new efficient treatments is urgently needed. Mycobacterium tuberculosis cytochrome P450 CYP121A1 is a promising drug target for the treatment of tuberculosis owing to its essential role in mycobacterial growth. Using a rational approach, which includes molecular modelling studies, three series of azole pyrazole derivatives were designed through two synthetic pathways. The synthesized compounds were biologically evaluated for their inhibitory activity towards M. tuberculosis and their protein binding affinity (K D). Series 3 biarylpyrazole imidazole derivatives were the most effective with the isobutyl (10 f) and tert-butyl (10 g) compounds displaying optimal activity (MIC 1.562 μg/mL, K D 0.22 μM (10 f) and 4.81 μM (10 g)). The spectroscopic data showed that all the synthesised compounds produced a type II red shift of the heme Soret band indicating either direct binding to heme iron or (where less extensive Soret shifts are observed) putative indirect binding via an interstitial water molecule. Evaluation of biological and physicochemical properties identified the following as requirements for activity: LogP >4, H-bond acceptors/H-bond donors 4/0, number of rotatable bonds 5-6, molecular volume >340 Å3, topological polar surface area <40 Å2.

Synthesis and biological evaluation of novel cYY analogues targeting Mycobacterium tuberculosis CYP121A1

The rise in multidrug resistant (MDR) cases of tuberculosis (TB) has led to the need for the development of TB drugs with different mechanisms of action. The genome sequence of Mycobacterium tuberculosis (Mtb) revealed twenty different genes coding for cytochrome P450s. CYP121A1 catalyzes a CC crosslinking reaction of dicyclotyrosine (cYY) producing mycocyclosin and current research suggests that either mycocyclosin is essential or the overproduction of cYY is toxic to Mtb. A series of 1,4-dibenzyl-2-imidazol-1-yl-methylpiperazine derivatives were designed and synthesised as cYY mimics. The derivatives substituted in the 4-position of the phenyl rings with halides or alkyl group showed promising antimycobacterial activity (MIC 6.25 μg/mL), with the more lipophilic branched alkyl derivatives displaying optimal binding affinity with CYP121A1 (ⁱPr K_D = 1.6 μM; ^tBu K_D = 1.2 μM). Computational studies revealed two possible binding modes within the CYP121A1 active site both of which would effectively block cYY from binding.

Novel insights into P450 BM3 interactions with FDA-approved antifungal azole drugs

Flavocytochrome P450 BM3 is a natural fusion protein constructed of cytochrome P450 and NADPH-cytochrome P450 reductase domains. P450 BM3 binds and oxidizes several mid- to long-chain fatty acids, typically hydroxylating these lipids at the ω-1, ω-2 and ω-3 positions. However, protein engineering has led to variants of this enzyme that are able to bind and oxidize diverse compounds, including steroids, terpenes and various human drugs. The wild-type P450 BM3 enzyme binds inefficiently to many azole antifungal drugs. However, we show that the BM3 A82F/F87V double mutant (DM) variant binds substantially tighter to numerous azole drugs than does the wild-type BM3, and that their binding occurs with more extensive heme spectral shifts indicative of complete binding of several azoles to the BM3 DM heme iron. We report here the first crystal structures of P450 BM3 bound to azole antifungal drugs - with the BM3 DM heme domain bound to the imidazole drugs clotrimazole and tioconazole, and to the triazole drugs fluconazole and voriconazole. This is the first report of any protein structure bound to the azole drug tioconazole, as well as the first example of voriconazole heme iron ligation through a pyrimidine nitrogen from its 5-fluoropyrimidine ring.

Cytochrome P450 1A1 opens up to new substrates

The cytochromes P450 (CYPs) oxidatively transform a huge number of substrates in both prokaryotic and eukaryotic organisms, but the mechanisms by which they accommodate these diverse molecules remain unclear. A new study by Bart and Scott reports two co-crystal structures of CYP1A1 that reveal structural rearrangements and flexible interaction networks that explain how the active site cavity shapes itself around new ligands. These data open the door to an increased understanding of fundamental enzyme behavior and improved searches for anti-cancer compounds.

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.

Characterization of Cytochrome P450 Enzymes and Their Applications in Synthetic Biology

The cytochrome P450 monooxygenase enzymes (P450s) catalyze a diverse array of chemical transformations, most originating from the insertion of an oxygen atom into a substrate that binds close to the P450 heme. The oxygen is delivered by a highly reactive heme iron-oxo species (compound I) and, according to the chemical nature of the substrate and its position in the active site, the P450 can catalyze a wide range of reactions including, e.g., hydroxylation, reduction, decarboxylation, sulfoxidation, N- and O-demethylation, epoxidation, deamination, CC bond formation and breakage, nitration, and dehalogenation. In this chapter, we describe the structural, biochemical, and catalytic properties of the P450s, along with spectroscopic and analytical methods used to characterize P450 enzymes and their redox partners. Important uses of P450 enzymes are highlighted, including how various P450s have been exploited for applications in synthetic biology.

Novel Aryl Substituted Pyrazoles as Small Molecule Inhibitors of Cytochrome P450 CYP121A1: Synthesis and Antimycobacterial Evaluation

Three series of biarylpyrazole imidazole and triazoles are described, which vary in the linker between the biaryl pyrazole and imidazole/triazole group. The imidazole and triazole series with the short -CH₂- linker displayed promising antimycobacterial activity, with the imidazole-CH₂- series (7) showing low MIC values (6.25-25 μg/mL), which was also influenced by lipophilicity. Extending the linker to -C(O)NH(CH₂)₂- resulted in a loss of antimycobacterial activity. The binding affinity of the compounds with CYP121A1 was determined by UV-visible optical titrations with K_D values of 2.63, 35.6, and 290 μM, respectively, for the tightest binding compounds 7e, 8b, and 13d from their respective series. Both binding affinity assays and docking studies of the CYP121A1 inhibitors suggest type II indirect binding through interstitial water molecules, with key binding residues Thr77, Val78, Val82, Val83, Met86, Ser237, Gln385, and Arg386, comparable with the binding interactions observed with fluconazole and the natural substrate dicyclotyrosine.

Expression, Purification, and Biochemical Characterization of the Flavocytochrome P450 CYP505A30 from Myceliophthora thermophila

The cytochrome P450/P450 reductase fusion enzyme CYP505A30 from the thermophilic fungus Myceliophthora thermophila and its heme (P450) domain were expressed in Escherichia coli and purified using affinity, ion exchange, and size exclusion chromatography. CYP505A30 binds straight chain fatty acids (from ∼C10 to C20), with highest affinity for tridecanoic acid (K_D = 2.7 μM). Reduced nicotinamide adenine dinucleotide phosphate is the preferred reductant for CYP505A30 (K_M = 3.1 μM compared to 330 μM for reduced nicotinamide adenine dinucleotide in cytochrome c reduction). Electron paramagnetic resonance confirmed cysteine thiolate coordination of heme iron in CYP505A30 and its heme domain. Redox potentiometry revealed an unusually positive midpoint potential for reduction of the flavin adenine dinucleotide and flavin mononucleotide cofactors (E⁰' ∼ -118 mV), and a large increase in the CYP505A30 heme domain Fe^III/Fe^II redox couple (ca. 230 mV) on binding arachidonic acid substrate. This switch brings the ferric heme iron potential into the same range as that of the reductase flavins. Multiangle laser light scattering analysis revealed CYP505A30's ability to dimerize, whereas the heme domain is monomeric. These data suggest CYP505A30 may function catalytically as a dimer (as described for Bacillus megaterium P450 BM3), and that binding interactions between CYP505A30 heme domains are not required for dimer formation. CYP505A30 catalyzed hydroxylation of straight chain fatty acids at the ω-1 to ω-3 positions, with a strong preference for ω-1 over ω-3 hydroxylation in the oxidation of dodecanoic and tetradecanoic acids (88 vs 2% products and 63 vs 9% products, respectively). CYP505A30 has important structural and catalytic similarities to P450 BM3 but distinct regioselectivity of lipid substrate oxidation with potential biotechnological applications.

Fragment Profiling Approach to Inhibitors of the Orphan M. tuberculosis P450 CYP144A1

Similarity between the ligand binding profiles of enzymes may aid functional characterization and be of greater relevance to inhibitor development than sequence similarity or structural homology. Fragment screening is an efficient approach for characterization of the ligand binding profile of an enzyme and has been applied here to study the family of cytochrome P450 enzymes (P450s) expressed by Mycobacterium tuberculosis (Mtb). The Mtb P450s have important roles in bacterial virulence, survival, and pathogenicity. Comparing the fragment profiles of seven of these enzymes revealed that P450s which share a similar biological function have significantly similar fragment profiles, whereas functionally unrelated or orphan P450s exhibit distinct ligand binding properties, despite overall high structural homology. Chemical structures that exhibit promiscuous binding between enzymes have been identified, as have selective fragments that could provide leads for inhibitor development. The similarity between the fragment binding profiles of the orphan enzyme CYP144A1 and CYP121A1, a characterized enzyme that is important for Mtb viability, provides a case study illustrating the subsequent identification of novel CYP144A1 ligands. The different binding modes of these compounds to CYP144A1 provide insight into structural and dynamic aspects of the enzyme, possible biological function, and provide the opportunity to develop inhibitors. Expanding this fragment profiling approach to include a greater number of functionally characterized and orphan proteins may provide a valuable resource for understanding enzyme-ligand interactions.

Catalytic Determinants of Alkene Production by the Cytochrome P450 Peroxygenase OleTJE

The Jeotgalicoccus sp. peroxygenase cytochrome P450 OleT_JE (CYP152L1) is a hydrogen peroxide-driven oxidase that catalyzes oxidative decarboxylation of fatty acids, producing terminal alkenes with applications as fine chemicals and biofuels. Understanding mechanisms that favor decarboxylation over fatty acid hydroxylation in OleT_JE could enable protein engineering to improve catalysis or to introduce decarboxylation activity into P450s with different substrate preferences. In this manuscript, we have focused on OleT_JE active site residues Phe⁷⁹, His⁸⁵, and Arg²⁴⁵ to interrogate their roles in substrate binding and catalytic activity. His⁸⁵ is a potential proton donor to reactive iron-oxo species during substrate decarboxylation. The H85Q mutant substitutes a glutamine found in several peroxygenases that favor fatty acid hydroxylation. H85Q OleT_JE still favors alkene production, suggesting alternative protonation mechanisms. However, the mutant undergoes only minor substrate binding-induced heme iron spin state shift toward high spin by comparison with WT OleT_JE, indicating the key role of His⁸⁵ in this process. Phe⁷⁹ interacts with His⁸⁵, and Phe⁷⁹ mutants showed diminished affinity for shorter chain (C10–C16) fatty acids and weak substrate-induced high spin conversion. F79A OleT_JE is least affected in substrate oxidation, whereas the F79W/Y mutants exhibit lower stability and cysteine thiolate protonation on reduction. Finally, Arg²⁴⁵ is crucial for binding the substrate carboxylate, and R245E/L mutations severely compromise activity and heme content, although alkene products are formed from some substrates, including stearic acid (C18:0). The results identify crucial roles for the active site amino acid trio in determining OleT_JE catalytic efficiency in alkene production and in regulating protein stability, heme iron coordination, and spin state.

Structural Characterization and Ligand/Inhibitor Identification Provide Functional Insights into the Mycobacterium tuberculosis Cytochrome P450 CYP126A1

The Mycobacterium tuberculosis H37Rv genome encodes 20 cytochromes P450, including P450s crucial to infection and bacterial viability. Many M. tuberculosis P450s remain uncharacterized, suggesting that their further analysis may provide new insights into M. tuberculosis metabolic processes and new targets for drug discovery. CYP126A1 is representative of a P450 family widely distributed in mycobacteria and other bacteria. Here we explore the biochemical and structural properties of CYP126A1, including its interactions with new chemical ligands. A survey of azole antifungal drugs showed that CYP126A1 is inhibited strongly by azoles containing an imidazole ring but not by those tested containing a triazole ring. To further explore the molecular preferences of CYP126A1 and search for probes of enzyme function, we conducted a high throughput screen. Compounds containing three or more ring structures dominated the screening hits, including nitroaromatic compounds that induce substrate-like shifts in the heme spectrum of CYP126A1. Spectroelectrochemical measurements revealed a 155-mV increase in heme iron potential when bound to one of the newly identified nitroaromatic drugs. CYP126A1 dimers were observed in crystal structures of ligand-free CYP126A1 and for CYP126A1 bound to compounds discovered in the screen. However, ketoconazole binds in an orientation that disrupts the BC-loop regions at the P450 dimer interface and results in a CYP126A1 monomeric crystal form. Structural data also reveal that nitroaromatic ligands "moonlight" as substrates by displacing the CYP126A1 distal water but inhibit enzyme activity. The relatively polar active site of CYP126A1 distinguishes it from its most closely related sterol-binding P450s in M. tuberculosis, suggesting that further investigations will reveal its diverse substrate selectivity.

An oxidative N-demethylase reveals PAS transition from ubiquitous sensor to enzyme

The universal Per-ARNT-Sim (PAS) domain functions as a signal transduction module involved in sensing diverse stimuli such as small molecules, light, redox state and gases. The highly evolvable PAS scaffold can bind a broad range of ligands, including haem, flavins and metal ions. However, although these ligands can support catalytic activity, to our knowledge no enzymatic PAS domain has been found. Here we report characterization of the first PAS enzyme: a haem-dependent oxidative N-demethylase. Unrelated to other amine oxidases, this enzyme contains haem, flavin mononucleotide, 2Fe-2S and tetrahydrofolic acid cofactors, and specifically catalyses the NADPH-dependent oxidation of dimethylamine. The structure of the α subunit reveals that it is a haem-binding PAS domain, similar in structure to PAS gas sensors. The dimethylamine substrate forms part of a highly polarized oxygen-binding site, and directly assists oxygen activation by acting as both an electron and proton donor. Our data reveal that the ubiquitous PAS domain can make the transition from sensor to enzyme, suggesting that the PAS scaffold can support the development of artificial enzymes.

Drug targeting of heme proteins in Mycobacterium tuberculosis

TB, caused by the human pathogen Mycobacterium tuberculosis (Mtb), causes more deaths than any other infectious disease. Iron is crucial for Mtb to infect the host and to sustain infection, with Mtb encoding large numbers of iron-binding proteins. Many of these are hemoproteins with key roles, including defense against oxidative stress, cellular signaling and regulation, host cholesterol metabolism, and respiratory processes. Various heme enzymes in Mtb are validated drug targets and/or products of genes essential for bacterial viability or survival in the host. Here, we review the structure, function, and druggability of key Mtb heme enzymes and strategies used for their inhibition.

Substrate Fragmentation for the Design of M. tuberculosis CYP121 Inhibitors

The cyclo-dipeptide substrates of the essential M. tuberculosis (Mtb) enzyme CYP121 were deconstructed into their component fragments and screened against the enzyme. A number of hits were identified, one of which exhibited an unexpected inhibitor-like binding mode. The inhibitory pharmacophore was elucidated, and fragment binding affinity was rapidly improved by synthetic elaboration guided by the structures of CYP121 substrates. The resulting inhibitors have low micromolar affinity, good predicted physicochemical properties and selectivity for CYP121 over other Mtb P450s. Spectroscopic characterisation of the inhibitors' binding mode provides insight into the effect of weak nitrogen-donor ligands on the P450 heme, an improved understanding of factors governing CYP121-ligand recognition and speculation into the biological role of the enzyme for Mtb.

Structural characterization of CYP144A1 - a cytochrome P450 enzyme expressed from alternative transcripts in Mycobacterium tuberculosis

Mycobacterium tuberculosis (Mtb) causes the disease tuberculosis (TB). The virulent Mtb H37Rv strain encodes 20 cytochrome P450 (CYP) enzymes, many of which are implicated in Mtb survival and pathogenicity in the human host. Bioinformatics analysis revealed that CYP144A1 is retained exclusively within the Mycobacterium genus, particularly in species causing human and animal disease. Transcriptomic annotation revealed two possible CYP144A1 start codons, leading to expression of (i) a “full-length” 434 amino acid version (CYP144A1-FLV) and (ii) a “truncated” 404 amino acid version (CYP144A1-TRV). Computational analysis predicted that the extended N-terminal region of CYP144A1-FLV is largely unstructured. CYP144A1 FLV and TRV forms were purified in heme-bound states. Mass spectrometry confirmed production of intact, His₆-tagged forms of CYP144A1-FLV and -TRV, with EPR demonstrating cysteine thiolate coordination of heme iron in both cases. Hydrodynamic analysis indicated that both CYP144A1 forms are monomeric. CYP144A1-TRV was crystallized and the first structure of a CYP144 family P450 protein determined. CYP144A1-TRV has an open structure primed for substrate binding, with a large active site cavity. Our data provide the first evidence that Mtb produces two different forms of CYP144A1 from alternative transcripts, with CYP144A1-TRV generated from a leaderless transcript lacking a 5′-untranslated region and Shine-Dalgarno ribosome binding site.

Fragment-Based Approaches to the Development of Mycobacterium tuberculosis CYP121 Inhibitors

The essential enzyme CYP121 is a target for drug development against antibiotic resistant strains of Mycobacterium tuberculosis. A triazol-1-yl phenol fragment 1 was identified to bind to CYP121 using a cascade of biophysical assays. Synthetic merging and optimization of 1 produced a 100-fold improvement in binding affinity, yielding lead compound 2 (KD = 15 μM). Deconstruction of 2 into its component retrofragments allowed the group efficiency of structural motifs to be assessed, the identification of more LE scaffolds for optimization and highlighted binding affinity hotspots. Structure-guided addition of a metal-binding pharmacophore onto LE retrofragment scaffolds produced low nanomolar (KD = 15 nM) CYP121 ligands. Elaboration of these compounds to target binding hotspots in the distal active site afforded compounds with excellent selectivity against human drug-metabolizing P450s. Analysis of the factors governing ligand potency and selectivity using X-ray crystallography, UV-vis spectroscopy, and native mass spectrometry provides insight for subsequent drug development.

Applications of microbial cytochrome P450 enzymes in biotechnology and synthetic biology

Cytochrome P450 enzymes (P450s) are a superfamily of monooxygenase enzymes with enormous potential for synthetic biology applications. Across Nature, their substrate range is vast and exceeds that of other enzymes. The range of different chemical transformations performed by P450s is also substantial, and continues to expand through interrogation of the properties of novel P450s and by protein engineering studies. The ability of P450s to introduce oxygen atoms at specific positions on complex molecules makes these enzymes particularly valuable for applications in synthetic biology. This review focuses on the enzymatic properties and reaction mechanisms of P450 enzymes, and on recent studies that highlight their broad applications in the production of oxychemicals. For selected soluble bacterial P450s (notably the high-activity P450-cytochrome P450 reductase enzyme P450 BM3), variants with a multitude of diverse substrate selectivities have been generated both rationally and by random mutagenesis/directed evolution approaches. This highlights the robustness and malleability of the P450 fold, and the capacity of these biocatalysts to oxidise a wide range of chemical scaffolds. This article reviews recent research on the application of wild-type and engineered P450s in the production of important chemicals, including pharmaceuticals and drug metabolites, steroids and antibiotics. In addition, the properties of unusual members of the P450 superfamily that do not follow the canonical P450 catalytic pathway are described.

The structure, function and properties of sirohaem decarboxylase--an enzyme with structural homology to a transcription factor family that is part of the alternative haem biosynthesis pathway

Some bacteria and archaea synthesize haem by an alternative pathway, which involves the sequestration of sirohaem as a metabolic intermediate rather than as a prosthetic group. Along this pathway the two acetic acid side-chains attached to C12 and C18 are decarboxylated by sirohaem decarboxylase, a heterodimeric enzyme composed of AhbA and AhbB, to give didecarboxysirohaem. Further modifications catalysed by two related radical SAM enzymes, AhbC and AhbD, transform didecarboxysirohaem into Fe-coproporphyrin III and haem respectively. The characterization of sirohaem decarboxylase is reported in molecular detail. Recombinant versions of Desulfovibrio desulfuricans, Desulfovibrio vulgaris and Methanosarcina barkeri AhbA/B have been produced and their physical properties compared. The D. vulgaris and M. barkeri enzyme complexes both copurify with haem, whose redox state influences the activity of the latter. The kinetic parameters of the D. desulfuricans enzyme have been determined, the enzyme crystallized and its structure has been elucidated. The topology of the enzyme reveals that it shares a structural similarity to the AsnC/Lrp family of transcription factors. The active site is formed in the cavity between the two subunits and a AhbA/B-product complex with didecarboxysirohaem has been obtained. A mechanism for the decarboxylation of the kinetically stable carboxyl groups is proposed.

Biofragments: an approach towards predicting protein function using biologically related fragments and its application to Mycobacterium tuberculosis CYP126

We present a novel fragment-based approach that tackles some of the challenges for chemical biology of predicting protein function. The general approach, which we have termed biofragments, comprises two key stages. First, a biologically relevant fragment library (biofragment library) can be designed and constructed from known sets of substrate-like ligands for a protein class of interest. Second, the library can be screened for binding to a novel putative ligand-binding protein from the same or similar class, and the characterization of hits provides insight into the basis of ligand recognition, selectivity, and function at the substrate level. As a proof-of-concept, we applied the biofragments approach to the functionally uncharacterized Mycobacterium tuberculosis (Mtb) cytochrome P450 isoform, CYP126. This led to the development of a tailored CYP biofragment library with notable 3D characteristics and a significantly higher screening hit rate (14%) than standard drug-like fragment libraries screened previously against Mtb CYP121 and 125 (4% and 1%, respectively). Biofragment hits were identified that make both substrate-like type-I and inhibitor-like type-II interactions with CYP126. A chemical-fingerprint-based substrate model was built from the hits and used to search a virtual TB metabolome, which led to the discovery that CYP126 has a strong preference for the recognition of aromatics and substrate-like type-I binding of chlorophenol moieties within the active site near the heme. Future catalytic analyses will be focused on assessing CYP126 for potential substrate oxidative dehalogenation.

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.

Structure and biochemical properties of the alkene producing cytochrome P450 OleTJE (CYP152L1) from the Jeotgalicoccus sp. 8456 bacterium

The production of hydrocarbons in nature has been documented for only a limited set of organisms, with many of the molecular components underpinning these processes only recently identified. There is an obvious scope for application of these catalysts and engineered variants thereof in the future production of biofuels. Here we present biochemical characterization and crystal structures of a cytochrome P450 fatty acid peroxygenase: the terminal alkene forming OleTJE (CYP152L1) from Jeotgalicoccus sp. 8456. OleTJE is stabilized at high ionic strength, but aggregation and precipitation of OleTJE in low salt buffer can be turned to advantage for purification, because resolubilized OleTJE is fully active and extensively dissociated from lipids. OleTJE binds avidly to a range of long chain fatty acids, and structures of both ligand-free and arachidic acid-bound OleTJE reveal that the P450 active site is preformed for fatty acid binding. OleTJE heme iron has an unusually positive redox potential (-103 mV versus normal hydrogen electrode), which is not significantly affected by substrate binding, despite extensive conversion of the heme iron to a high spin ferric state. Terminal alkenes are produced from a range of saturated fatty acids (C12-C20), and stopped-flow spectroscopy indicates a rapid reaction between peroxide and fatty acid-bound OleTJE (167 s(-1) at 200 μm H2O2). Surprisingly, the active site is highly similar in structure to the related P450BSβ, which catalyzes hydroxylation of fatty acids as opposed to decarboxylation. Our data provide new insights into structural and mechanistic properties of a robust P450 with potential industrial applications.

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.

Heme sensor proteins

Heme is a prosthetic group best known for roles in oxygen transport, oxidative catalysis, and respiratory electron transport. Recent years have seen the roles of heme extended to sensors of gases such as O2 and NO and cell redox state, and as mediators of cellular responses to changes in intracellular levels of these gases. The importance of heme is further evident from identification of proteins that bind heme reversibly, using it as a signal, e.g. to regulate gene expression in circadian rhythm pathways and control heme synthesis itself. In this minireview, we explore the current knowledge of the diverse roles of heme sensor proteins.

Unusual cytochrome p450 enzymes and reactions

Cytochrome P450 enzymes primarily catalyze mixed-function oxidation reactions, plus some reductions and rearrangements of oxygenated species, e.g. prostaglandins. Most of these reactions can be rationalized in a paradigm involving Compound I, a high-valent iron-oxygen complex (FeO(3+)), to explain seemingly unusual reactions, including ring couplings, ring expansion and contraction, and fusion of substrates. Most P450s interact with flavoenzymes or iron-sulfur proteins to receive electrons from NAD(P)H. In some cases, P450s are fused to protein partners. Other P450s catalyze non-redox isomerization reactions. A number of permutations on the P450 theme reveal the diversity of cytochrome P450 form and function.

Bacillus megaterium has both a functional BluB protein required for DMB synthesis and a related flavoprotein that forms a stable radical species

Despite the extensive study of the biosynthesis of the complex molecule B12 (cobalamin), the mechanism by which the lower ligand 5,6-dimethylbenzimidazole (DMB) is formed has remained something of a mystery. However, recent work has identified and characterized a DMB-synthase (BluB) responsible for the oxygen-dependent, single enzyme conversion of FMN to DMB. In this work, we have identified BluB homologs from the aerobic purple, nonsulfur, photosynthetic bacterium Rhodobacter capsulatus and the aerobic soil bacterium Bacillus megaterium and have demonstrated DMB synthesis by the use of a novel complementation assay in which a B12 deficient strain, substituted with the precursor cobinamide is recovered either by the addition of DMB or by the recombinant expression of a bluB gene. The DMB-synthetic activity of the purified recombinant BluB enzymes was further confirmed in vitro by providing the enzyme with FMNH2 and oxygen and observing the formation of DMB by HPLC. The formation of a 4a-peroxyflavin intermediate, the first step in the oxygen dependent mechanism of DMB biosynthesis, is reported here and is the first intermediate in the enzyme catalysed reaction to be demonstrated experimentally to date. The identification and characterization of an FMN-binding protein found on the cobI operon of B. megaterium, CbiY, is also detailed, revealing an FMN-containing enzyme which is able to stabilize a blue flavin semiquinone upon reduction with a 1-electron donor.

Unusual spectroscopic and ligand binding properties of the cytochrome P450-flavodoxin fusion enzyme XplA

The Rhodococcus rhodochrous strain 11Y XplA enzyme is an unusual cytochrome P450-flavodoxin fusion enzyme that catalyzes reductive denitration of the explosive hexahydro-1,3,5-trinitro-1,3,5-triazene (RDX). We show by light scattering that XplA is a monomeric enzyme. XplA has high affinity for imidazole (K(d) = 1.6 μM), explaining previous reports of a red-shifted XplA Soret band in pure enzyme. The true Soret maximum of XplA is at 417 nm. Similarly, unusually weak XplA flavodoxin FMN binding (K(d) = 1.09 μM) necessitates its purification in the presence of the cofactor to produce hallmark flavin contributions absent in previously reported spectra. Structural and ligand-binding data reveal a constricted active site able to accommodate RDX and small inhibitory ligands (e.g. 4-phenylimidazole and morpholine) while discriminating against larger azole drugs. The crystal structure also identifies a high affinity imidazole binding site, consistent with its low K(d), and shows active site penetration by PEG, perhaps indicative of an evolutionary lipid-metabolizing function for XplA. EPR studies indicate heterogeneity in binding mode for RDX and other ligands. The substrate analog trinitrobenzene does not induce a substrate-like type I optical shift but creates a unique low spin EPR spectrum due to influence on structure around the distal water heme ligand. The substrate-free heme iron potential (-268 mV versus NHE) is positive for a low spin P450, and the elevated potential of the FMN semiquinone/hydroquinone couple (-172 mV) is also an adaptation that may reflect (along with the absence of a key Thr/Ser residue conserved in oxygen-activating P450s) the evolution of XplA as a specialized RDX reductase catalyst.

Special issue: Cytochrome P450 structure and function: introduction

The 17th International Conference on Cytochrome P450 Biochemistry, Biophysics and Structure was held in Manchester, UK from 26-30 June 2011. This issue of FEBS J. contains review and primary research articles reflecting the breadth of science covered at this conference, and reflecting the impact of P450-related research in fields as diverse as steroid metabolism, plant biochemistry, structural biology and biotechnology.

Characterization of Cupriavidus metallidurans CYP116B1--a thiocarbamate herbicide oxygenating P450-phthalate dioxygenase reductase fusion protein

The novel cytochrome P450/redox partner fusion enzyme CYP116B1 from Cupriavidus metallidurans was expressed in and purified from Escherichia coli. Isolated CYP116B1 exhibited a characteristic Fe(II)CO complex with Soret maximum at 449 nm. EPR and resonance Raman analyses indicated low-spin, cysteinate-coordinated ferric haem iron at both 10 K and ambient temperature, respectively, for oxidized CYP116B1. The EPR of reduced CYP116B1 demonstrated stoichiometric binding of a 2Fe-2S cluster in the reductase domain. FMN binding in the reductase domain was confirmed by flavin fluorescence studies. Steady-state reduction of cytochrome c and ferricyanide were supported by both NADPH/NADH, with NADPH used more efficiently (K(m[NADPH]) = 0.9 ± 0.5 μM and K(m[NADH]) = 399.1 ± 52.1 μM). Stopped-flow studies of NAD(P)H-dependent electron transfer to the reductase confirmed the preference for NADPH. The reduction potential of the P450 haem iron was -301 ± 7 mV, with retention of haem thiolate ligation in the ferrous enzyme. Redox potentials for the 2Fe-2S and FMN cofactors were more positive than that of the haem iron. Multi-angle laser light scattering demonstrated CYP116B1 to be monomeric. Type I (substrate-like) binding of selected unsaturated fatty acids (myristoleic, palmitoleic and arachidonic acids) was shown, but these substrates were not oxidized by CYP116B1. However, CYP116B1 catalysed hydroxylation (on propyl chains) of the herbicides S-ethyl dipropylthiocarbamate (EPTC) and S-propyl dipropylthiocarbamate (vernolate), and the subsequent N-dealkylation of vernolate. CYP116B1 thus has similar thiocarbamate-oxidizing catalytic properties to Rhodoccocus erythropolis CYP116A1, a P450 involved in the oxidative degradation of EPTC.

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.

Analysis of the oxidation of short chain alkynes by flavocytochrome P450 BM3

Bacillus megaterium flavocytochrome P450 BM3 (BM3) is a high activity fatty acid hydroxylase, formed by the fusion of soluble cytochrome P450 and cytochrome P450 reductase modules. Short chain (C6, C8) alkynes were shown to be substrates for BM3, with productive outcomes (i.e. alkyne hydroxylation) dependent on position of the carbon-carbon triple bond in the molecule. Wild-type P450 BM3 catalyses ω-3 hydroxylation of both 1-hexyne and 1-octyne, but is suicidally inactivated in NADPH-dependent turnover with non-terminal alkynes. A F87G mutant of P450 BM3 also undergoes turnover-dependent heme destruction with the terminal alkynes, pointing to a key role for Phe87 in controlling regioselectivity of alkyne oxidation. The terminal alkynes access the BM3 heme active site led by the acetylene functional group, since hydroxylated products are not observed near the opposite end of the molecules. For both 1-hexyne and 1-octyne, the predominant enantiomeric product formed (up to ∼90%) is the (S)-(-)-1-alkyn-3-ol form. Wild-type P450 BM3 is shown to be an effective oxidase catalyst of terminal alkynes, with strict regioselectivity of oxidation and potential biotechnological applications. The absence of measurable octanoic or hexanoic acid products from oxidation of the relevant 1-alkynes is also consistent with previous studies suggesting that removal of the phenyl group in the F87G mutant does not lead to significant levels of ω-oxidation of alkyl chain substrates.

A novel intermediate in the reaction of seleno CYP119 with m-chloroperbenzoic acid

Cytochrome P450-mediated monooxygenation generally proceeds via a reactive ferryl intermediate coupled to a ligand radical [Fe(IV)═O]+• termed Compound I (Cpd I). The proximal cysteine thiolate ligand is a critical determinant of the spectral and catalytic properties of P450 enzymes. To explore the effect of an increased level of donation of electrons by the proximal ligand in the P450 catalytic cycle, we recently reported successful incorporation of SeCys into the active site of CYP119, a thermophilic cytochrome P450. Here we report relevant physical properties of SeCYP119 and a detailed analysis of the reaction of SeCYP119 with m-chloroperbenzoic acid. Our results indicate that the selenolate anion reduces rather than stabilizes Cpd I and also protects the heme from oxidative destruction, leading to the generation of a new stable species with an absorbance maximum at 406 nm. This stable intermediate can be returned to the normal ferric state by reducing agents and thiols, in agreement with oxidative modification of the selenolate ligand itself. Thus, in the seleno protein, the oxidative damage shifts from the heme to the proximal ligand, presumably because (a) an increased level of donation of electrons more efficiently quenches reactive species such as Cpd I and (b) the protection of the thiolate ligand provided by the protein active site structure is insufficient to shield the more oxidizable selenolate ligand.