Crystal structure of a thermophilic cytochrome P450 from the archaeon Sulfolobus solfataricus

Predicting flexible loop regions that interact with ligands: the challenge of accurate scoring

Flexible loop regions play a critical role in the biological function of many proteins and have been shown to be involved in ligand binding. In the context of structure-based drug design, using or predicting an incorrect loop configuration can be detrimental to the study if the loop is capable of interacting with the ligand. Three protein systems, each with at least one flexible loop region in close proximity to the known binding site, were selected for loop prediction using the CorLps program; a six residue loop region from phosphoribosylglycinamide formyltransferase (GART), two nine residue loop regions from cytochrome P450 (CYP) 119, and an 11 residue loop region from enolase were selected for loop prediction. The results of this study indicate that the statistically based DFIRE scoring function implemented in the CorLps program did not accurately rank native-like predicted loop configurations in any protein system. In an attempt to improve the ranking of the native-like predicted loop configurations, the MM/GBSA and the optimized MM/GBSA-dsr scoring functions were used to re-rank the predicted loops with and without bound ligand. In general, single snapshot MM/GBSA scoring provided the best ranking of native-like loop configurations. Based on the scoring function analyses presented, the optimal ranking of native-like loop configurations is still a difficult challenge and the choice of the "best" scoring function appears to be system dependent.

Potentially increasing the metabolic stability of drug candidates via computational site of metabolism prediction by CYP2C9: The utility of incorporating protein flexibility via an ensemble of structures

Cytochrome P450 enzymes are responsible for metabolizing many endogenous and xenobiotic molecules encountered by the human body. It has been estimated that 75% of all drugs are metabolized by cytochrome P450 enzymes. Thus, predicting a compound's potential sites of metabolism (SOM) is highly advantageous early in the drug development process. We have combined molecular dynamics, AutoDock Vina docking, the neighboring atom type (NAT) reactivity model, and a solvent-accessible surface-area term to form a reactivity-accessibility model capable of predicting SOM for cytochrome P450 2C9 substrates. To investigate the importance of protein flexibility during the ligand-binding process, the results of SOM prediction using a static protein structure for docking were compared to SOM prediction using multiple protein structures in ensemble docking. The results reported here indicate that ensemble docking increases the number of ligands that can be docked in a bioactive conformation (ensemble: 96%, static: 85%) but only leads to a slight improvement (49% vs. 44%) in predicting an experimentally known SOM in the top-1 position for a ligand library of 75 CYP2C9 substrates. Using ensemble docking, the reactivity-accessibility model accurately predicts SOM in the top-1 ranked position for 49% of the ligand library and considering the top-3 predicted sites increases the prediction success rate to approximately 70% of the ligand library. Further classifying the substrate library according to K(m) values leads to an improvement in SOM prediction for substrates with low K(m) values (57% at top-1). While the current predictive power of the reactivity-accessibility model still leaves significant room for improvement, the results illustrate the usefulness of this method to identify key protein-ligand interactions and guide structural modifications of the ligand to increase its metabolic stability.

Efficient incorporation of protein flexibility and dynamics into molecular docking simulations

Flexibility and dynamics are protein characteristics that are essential for the process of molecular recognition. Conformational changes in the protein that are coupled to ligand binding are described by the biophysical models of induced fit and conformational selection. Different concepts that incorporate protein flexibility into protein-ligand docking within the context of these two models are reviewed. Several computational studies that discuss the validity and possible limitations of such approaches will be presented. Finally, different approaches that incorporate protein dynamics, e.g., configurational entropy, and solvation effects into docking will be highlighted.

Crystal structure of H2O2-dependent cytochrome P450SPalpha with its bound fatty acid substrate: insight into the regioselective hydroxylation of fatty acids at the alpha position

Cytochrome P450(SPα) (CYP152B1) isolated from Sphingomonas paucimobilis is the first P450 to be classified as a H(2)O(2)-dependent P450. P450(SPα) hydroxylates fatty acids with high α-regioselectivity. Herein we report the crystal structure of P450(SPα) with palmitic acid as a substrate at a resolution of 1.65 Å. The structure revealed that the C(α) of the bound palmitic acid in one of the alternative conformations is 4.5 Å from the heme iron. This conformation explains the highly selective α-hydroxylation of fatty acid observed in P450(SPα). Mutations at the active site and the F-G loop of P450(SPα) did not impair its regioselectivity. The crystal structures of mutants (L78F and F288G) revealed that the location of the bound palmitic acid was essentially the same as that in the WT, although amino acids at the active site were replaced with the corresponding amino acids of cytochrome P450(BSβ) (CYP152A1), which shows β-regioselectivity. This implies that the high regioselectivity of P450(SPα) is caused by the orientation of the hydrophobic channel, which is more perpendicular to the heme plane than that of P450(BSβ).

Active-site residues move independently from the rest of the protein in a 200 ns molecular dynamics simulation of cytochrome P450 CYP119

The conformational dynamics of cytochrome P450 enzymes are critical to their catalytic activity. In this study, the correlated motion between residues in a 200 ns molecular dynamics trajectory of the thermophilic CYP119 was analyzed to parse out conformational relationships. Residues that are structurally related, for example residues within a helix, generally have highly correlated motion. In addition, clusters of non-adjacent residues that show correlated motion ("hot spots") are seen in various regions, including at the base of the F and G helices that make up the most dynamic region of the enzyme. A modified k-means algorithm that clusters residues based on their correlated motion indicates that functionally related residues are in the same cluster (e.g., the catalytic threonines and the heme). Tightly coupled clusters form a solvent-exposed "shell" around the enzyme, whereas less coupling between clusters is seen in regions that are critical to ligand interactions, redox partner interactions, and catalysis. Most notably, we find that residues in the active site move independently from the rest of the enzyme, effectively insulating the catalytic machinery from other regions of the protein.

Zinc binding to the Tyr402 and His402 allotypes of complement factor H: possible implications for age-related macular degeneration

The Tyr402His polymorphism of complement factor H (FH) with 20 short complement regulator (SCR) domains is associated with age-related macular degeneration (AMD). How FH contributes to disease pathology is not clear. Both FH and high concentrations of zinc are found in drusen deposits, the key feature of AMD. Heterozygous FH is inhibited by zinc, which causes FH to aggregate. Here, zinc binding to homozygous FH was studied. By analytical ultracentrifugation, large amounts of oligomers were observed with both the native Tyr402 and the AMD-risk His402 homozygous allotypes of FH and both the recombinant SCR-6/8 allotypes with Tyr/His402. X-ray scattering also showed that both FH and SCR-6/8 allotypes strongly aggregated at >10 μM zinc. The SCR-1/5 and SCR-16/20 fragments were less likely to bind zinc. These observations were supported by bioinformatics predictions. Starting from known zinc binding sites in crystal structures, we predicted 202 putative partial surface zinc binding sites in FH, most of which were in SCR-6. Metal site prediction web servers also suggested that SCR-6 and other domains bind zinc. Predicted SCR-6/8 dimer structures showed that zinc binding sites could be formed at the protein-protein interface that would lead to daisy-chained oligomers. It was concluded that zinc binds weakly to FH at multiple surface locations, most probably within the functionally important SCR-6/8 domains, and this explains why zinc inhibits FH activity. Given the high pathophysiological levels of bioavailable zinc present in subretinal deposits, we discuss how zinc binding to FH may contribute to deposit formation and inflammation associated with AMD.

Protein dynamics in cytochrome P450 molecular recognition and substrate specificity using 2D IR vibrational echo spectroscopy

Cytochrome (cyt) P450s hydroxylate a variety of substrates that can differ widely in their chemical structure. The importance of these enzymes in drug metabolism and other biological processes has motivated the study of the factors that enable their activity on diverse classes of molecules. Protein dynamics have been implicated in cyt P450 substrate specificity. Here, 2D IR vibrational echo spectroscopy is employed to measure the dynamics of cyt P450(cam) from Pseudomonas putida on fast time scales using CO bound at the active site as a vibrational probe. The substrate-free enzyme and the enzyme bound to both its natural substrate, camphor, and a series of related substrates are investigated to explicate the role of dynamics in molecular recognition in cyt P450(cam) and to delineate how the motions may contribute to hydroxylation specificity. In substrate-free cyt P450(cam), three conformational states are populated, and the structural fluctuations within a conformational state are relatively slow. Substrate binding selectively stabilizes one conformational state, and the dynamics become faster. Correlations in the observed dynamics with the specificity of hydroxylation of the substrates, the binding affinity, and the substrates' molecular volume suggest that motions on the hundreds of picosecond time scale contribute to the variation in activity of cyt P450(cam) toward different substrates.

Experimentally restrained molecular dynamics simulations for characterizing the open states of cytochrome P450cam

Residual dipolar couplings (RDCs) were used as restraints in fully solvated molecular dynamics simulations of reduced substrate- and carbonmonoxy-bound cytochrome P450(cam) (CYP101A1), a 414-residue soluble monomeric heme-containing camphor monooxygenase from the soil bacterium Pseudomonas putida. The (1)D(NH) residual dipolar couplings used as restraints were measured in two independent alignment media. A soft annealing protocol was used to heat the starting structures while incorporating the RDC restraints. After production dynamics, structures with the lowest total violation energies for RDC restraints were extracted to identify ensembles of conformers accessible to the enzyme in solution. The simulations result in substrate orientations different from that seen in crystallographic structures and a more open and accessible enzyme active site and largely support previously reported differences between the open and closed states of CYP101A1.

The structure of CYP101D2 unveils a potential path for substrate entry into the active site

The cytochrome P450 CYP101D2 from Novosphingobium aromaticivorans DSM12444 is closely related to CYP101D1 from the same bacterium and to P450cam (CYP101A1) from Pseudomonas putida. All three are capable of oxidizing camphor stereoselectively to 5-exo-hydroxycamphor. The crystal structure of CYP101D2 revealed that the likely ferredoxin-binding site on the proximal face is largely positively charged, similar to that of CYP101D1. However, both the native and camphor-soaked forms of CYP101D2 had open conformations with an access channel. In the active site of the camphor-soaked form, the camphor carbonyl interacted with the haem-iron-bound water. Two other potential camphor-binding sites were also identified from electron densities in the camphor-soaked structure: one located in the access channel, flanked by the B/C and F/G loops and the I helix, and the other in a cavity on the surface of the enzyme near the F helix side of the F/G loop. The observed open structures may be conformers of the CYP101D2 enzyme that enable the substrate to enter the buried active site via a conformational selection mechanism. The second and third binding sites may be intermediate locations of substrate entry and translocation into the active site, and provide insight into a multi-step substrate-binding mechanism.

Characterization of cytochrome P450 monooxygenase CYP154H1 from the thermophilic soil bacterium Thermobifida fusca

Cytochrome P450 monooxygenases are valuable biocatalysts due to their ability to hydroxylate unactivated carbon atoms using molecular oxygen. We have cloned the gene for a new cytochrome P450 monooxygenase, named CYP154H1, from the moderately thermophilic soil bacterium Thermobifida fusca. The enzyme was overexpressed in Escherichia coli at up to 14% of total soluble protein and purified to homogeneity in three steps. CYP154H1 activity was reconstituted using putidaredoxin reductase and putidaredoxin from Pseudomonas putida DSM 50198 as surrogate electron transfer partners. In biocatalytic reactions with different aliphatic and aromatic substrates of varying size, the enzyme converted small aromatic and arylaliphatic compounds like ethylbenzene, styrene, and indole. Furthermore, CYP154H1 also accepted different arylaliphatic sulfides as substrates chemoselectively forming the corresponding sulfoxides and sulfones. The enzyme is moderately thermostable with an apparent melting temperature of 67°C and exhibited still 90% of initial activity after incubation at 50°C.

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.

Using molecular dynamics to probe the structural basis for enhanced stability in thermal stable cytochromes P450

High-temperature molecular dynamics (MD) has been used to assess if MD can be employed as a useful tool for probing the structural basis for enhanced stability in thermal stable cytochromes P450. CYP119, the most thermal stable P450 known, unfolds more slowly during 500 K MD simulations than P450s that melt at lower temperatures, P450cam and P450cin. A comparison of the 500 K MD trajectories shows that the Cys ligand loop, a critically important structural feature just under the heme, in both P450cin and P450cam completely unfolds while this region is quite stable in CYP119. In CYP119, this region is stabilized by tight nonpolar interactions involving Tyr26 and Leu308. The corresponding residues in P450cam are Gly and Thr, respectively. The in silico generated Y26A/L308A CYP119 double mutant is substantially less stable than wild-type CYP119, and the Cys ligand loop unfolds in a manner similar to that of P450cam. The MD thus has identified a potential "hot spot" important for stability. As an experimental test of the MD results, the Y26A/L308A double mutant was prepared, and thermal melting curves show that the double mutant exhibits a melting temperature (T(m)) 16 degrees C lower than that of wild-type CYP119. Control mutations that were predicted by MD not to destabilize the protein were also generated, and the experimental melting temperature was not significantly different from that of the wild-type enzyme. Therefore, high-temperature MD is a useful tool in predicting the structural underpinnings of thermal stability in P450s.

New cytochrome P450 mechanisms: implications for understanding molecular basis for drug toxicity at the level of the cytochrome

IMPORTANCE OF THE FIELD

Cytochrome (CYP) P450 is a collective name for a very large group of heme enzymes, which catalyze largely oxidative reactions, including those of pharmacological and toxicological importance. Their efficient operation requires coupling of specific electron donor and O(2) consumption and substrate hydroxylation. Many drug oxidation reactions are partially uncoupled, leading to the formation of highly toxic reactive oxygen species, which can cause unpredictable toxic effects on the cell. Rational approaches to avoid uncoupling require knowledge of the underlying mechanisms.

AREAS COVERED IN THIS REVIEW

In this communication, attempts have been made to bring together past as well as present information indicating that i) the P450 active site has two differently accessible allosterically interacting subsites geared for entirely different types of functionally relevant interactions; and ii) substrate binding to the specific protein residues (Site I) forming the reducible high-spin complex and product binding at L(6) (Site II) of the heme iron forming inhibited low-spin complex can regulate the functional state of the enzyme during catalysis.

WHAT THE READER WILL GAIN

Since P450 enzymes catalyze a wide variety of reactions, understanding the molecular basis for their efficient operation is of interest to many fields, including rational approaches to design safer drugs, tailoring P450 for a given task (e.g., bioremediation).

TAKE HOME MESSAGE

It is important to take into account that the two sub-sites function as interacting sites rather than parts of a site functioning as single site for rational approaches to P450 mechanisms. This is important especially in regard to interpretation of the observed effects of drugs, products and inhibitors on these enzymes.

The Mycobacterium tuberculosis cytochrome P450 system

Tuberculosis remains a leading cause of human mortality. The emergence of strains of Mycobacterium tuberculosis, the causative agent, that are resistant to the major frontline antitubercular drugs increases the urgency for the development of new therapeutic agents. Sequencing of the M. tuberculosis genome revealed the existence of 20 cytochrome P450 enzymes, some of which are potential candidates for drug targeting. The recent burst of studies reporting microarray-based gene essentiality and transcriptome analyses under in vitro, ex vivo and in vivo conditions highlight the importance of selected P450 isoforms for M. tuberculosis viability and pathogenicity. Current knowledge of the structural and biochemical properties of the M. tuberculosis P450 enzymes and their putative redox partners is reviewed, with an emphasis on findings related to their physiological function(s) as well as their potential as drug targets.

Cooperative properties of cytochromes P450

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

Crystal structures of cytochrome P450 2B4 in complex with the inhibitor 1-biphenyl-4-methyl-1H-imidazole: ligand-induced structural response through alpha-helical repositioning

Two different ligand occupancy structures of cytochrome P450 2B4 (CYP2B4) in complex with 1-biphenyl-4-methyl-1H-imidazole (1-PBI) have been determined by X-ray crystallography. 1-PBI belongs to a series of tight binding, imidazole-based CYP2B4 inhibitors. 1-PBI binding to CYP2B4 yields a type II spectrum with a K(s) value of 0.23 microM and inhibits enzyme activity with an IC(50) value of 0.035 microM. Previous CYP2B4 structures have shown a large degree of structural movement in response to ligand size. With two phenyl rings, 1-PBI is larger than 1-(4-chlorophenyl)imidazole (1-CPI) and 4-(4-chlorophenyl)imidazole (4-CPI) but smaller than bifonazole, which is branched and contains three phenyl rings. The CYP2B4-1-PBI complex is a structural intermediate to the closed CPI and the open bifonazole structures. The B/C-loop reorganizes itself to include two short partial helices while closing one side of the active site. The F-G-helix cassette pivots over the I-helix in direct response to the size of the ligand in the active site. A cluster of Phe residues at the fulcrum of this pivot point allows for dramatic repositioning of the cassette with only a relatively small amount of secondary structure rearrangement. Comparisons of ligand-bound CYP2B4 structures reveal trends in plastic region mobility that could allow for predictions of their position in future structures based on ligand shape and size.

Identification of PpoA from Aspergillus nidulans as a fusion protein of a fatty acid heme dioxygenase/peroxidase and a cytochrome P450

The homothallic ascomycete Aspergillus nidulans serves as model organism for filamentous fungi because of its ability to propagate with both asexual and sexual life cycles, and fatty acid-derived substances regulate the balance between both cycles. These so-called psi (precocious sexual inducer) factors are produced by psi factor-producing oxygenases (Ppo enzymes). Bioinformatic analysis predicted the presence of two different heme domains in Ppo proteins: in the N-terminal region, a fatty acid heme dioxygenase/peroxidase domain is predicted, whereas in the C-terminal region, a P450 heme thiolate domain is predicted. To analyze the reaction catalyzed by Ppo enzymes, PpoA was expressed in Escherichia coli as an active enzyme. The protein was purified by 62-fold and identified as a homotetrameric ferric heme protein that metabolizes mono- as well as polyunsaturated C(16) and C(18) fatty acids at pH approximately 7.25. The presence of thiolate-ligated heme was confirmed on the basis of sequence alignments and the appearance of a characteristic 450 nm CO-binding spectrum. Studies on its reaction mechanism revealed that PpoA uses different heme domains to catalyze two separate reactions. Within the heme peroxidase domain, linoleic acid is oxidized to (8R)-hydroperoxyoctadecadienoic acid by abstracting a H-atom from C-8 of the fatty acid, yielding a carbon-centered radical that reacts with molecular dioxygen. In the second reaction step, 8-hydroperoxyoctadecadienoic acid is isomerized within the P450 heme thiolate domain to 5,8-dihydroxyoctadecadienoic acid. We identify PpoA as a bifunctional P450 fusion protein that uses a previously unknown reaction mechanism for forming psi factors.

Crystal structures of cytochrome P450 105P1 from Streptomyces avermitilis: conformational flexibility and histidine ligation state

The polyene macrolide antibiotic filipin is widely used as a probe for cholesterol in biological membranes. The filipin biosynthetic pathway of Streptomyces avermitilis contains two position-specific hydroxylases, C26-specific CYP105P1 and C1'-specific CYP105D6. In this study, we describe the three X-ray crystal structures of CYP105P1: the ligand-free wild-type (WT-free), 4-phenylimidazole-bound wild-type (WT-4PI), and ligand-free H72A mutant (H72A-free) forms. The BC loop region in the WT-free structure has a unique feature; the side chain of His72 within this region is ligated to the heme iron. On the other hand, this region is highly disordered and widely open in WT-4PI and H72A-free structures, respectively. Histidine ligation of wild-type CYP105P1 was not detectable in solution, and a type II spectral change was clearly observed when 4-phenylimidazole was titrated. The H72A mutant showed spectroscopic characteristics that were almost identical to those of the wild-type protein. In the H72A-free structure, there is a large pocket that is of the same size as the filipin molecule. The highly flexible feature of the BC loop region of CYP105P1 may be required to accept a large hydrophobic substrate.

Ligand-induced conformational heterogeneity of cytochrome P450 CYP119 identified by 2D NMR spectroscopy with the unnatural amino acid (13)C-p-methoxyphenylalanine

Conformational dynamics are thought to play an important role in ligand binding and catalysis by cytochrome P450 enzymes, but few techniques exist to examine them in molecular detail. Using a unique isotopic labeling strategy, we have site specifically inserted a (13)C-labeled unnatural amino acid residue, (13)C-p-methoxyphenylalanine (MeOF), into two different locations in the substrate binding region of the thermophilic cytochrome P450 enzyme CYP119. Surprisingly, in both cases the resonance signal from the ligand-free protein is represented by a doublet in the (1)H,(13)C-HSQC spectrum. Upon binding of 4-phenylimidazole, the signals from the initial resonances are reduced in favor of a single new resonance, in the case of the F162MeOF mutant, or two new resonances, in the case of the F153MeOF mutant. This represents the first direct physical evidence for the ligand-dependent existence of multiple P450 conformers simultaneously in solution. This general approach may be used to further illuminate the role that conformational dynamics plays in the complex enzymatic phenomena exhibited by P450 enzymes.

Modes of heme binding and substrate access for cytochrome P450 CYP74A revealed by crystal structures of allene oxide synthase

Cytochrome P450s exist ubiquitously in all organisms and are involved in many biological processes. Allene oxide synthase (AOS) is a P450 enzyme that plays a key role in the biosynthesis of oxylipin jasmonates, which are involved in signal and defense reactions in higher plants. The crystal structures of guayule (Parthenium argentatum) AOS (CYP74A2) and its complex with the substrate analog 13(S)-hydroxyoctadeca-9Z,11E-dienoic acid have been determined. The structures exhibit a classic P450 fold but possess a heme-binding mode with an unusually long heme binding loop and a unique I-helix. The structures also reveal two channels through which substrate and product may access and leave the active site. The entrances are defined by a loop between beta3-2 and beta3-3. Asn-276 in the substrate binding site may interact with the substrate's hydroperoxy group and play an important role in catalysis, and Lys-282 at the entrance may control substrate access and binding. These studies provide both structural insights into AOS and related P450s and a structural basis to understand the distinct reaction mechanism.

Characterization of the peroxidase activity of CYP119, a thermostable P450 from Sulfolobus acidocaldarius

We report the cloning, expression, and purification of CYP119, a thermostable enzyme previously thought to derive from Sulfolobus solfataricus. Sequence analysis suggested that, in contrast to the conclusions of earlier studies, the enzyme stems from the closely related Sulfolobus acidocaldarius, and we were indeed able to clone the gene from the genomic DNA of this organism. For the first time, we report here on the peroxidase activity of this enzyme and the optimization of the associated reaction parameters. The optimized reaction conditions were then applied to the biocatalytic epoxidation of styrene. The values obtained for k(cat) (78.2+/-20.6 min(-1)) and K(M) (9.2+/-4.3 mM) indicated an approximately 100-fold increased catalytic activity over previously reported results.

Modulation of redox potential and alteration in reactivity via the peroxide shunt pathway by mutation of cytochrome P450 around the proximal heme ligand

In the thermophilic cytochrome P450 from the thermoacidophilic crenarchaeon Sulfolobus tokodaii strain 7 (P450st), a phenylalanine residue at position 310 and an alanine residue at position 320 are located close to the heme thiolate ligand, Cys317. Single site-directed mutants F310A and A320Q and double mutant F310A/A320Q have been constructed. All mutant enzymes as well as wild-type (WT) P450st were expressed at high levels. The substitution of F310 with Ala and of A320 with Gln induced shifts in redox potential and blue shifts in Soret absorption of ferrous-CO forms, while spectral characterization showed that in the resting state, the mutants almost retained the structural integrity of the active site. The redox potential of the heme varied as follows: -481 mV (WT), -477 mV (A320Q), -453 mV (F310A), and -450 mV (F310A/A320Q). The trend in the Soret band of the ferrous-CO form was as follows: 450 nm (WT) < 449 nm (A320Q) < 446 nm (F310A) < 444 nm (F310A/A320Q). These results established that the reduction potential and electron density on the heme iron are modulated by the Phe310 and Ala320 residues in P450st. The electron density on the heme decreases in the following order: WT > A320Q > F310A > F310A/A320Q. The electron density on the heme iron infers an essential role in P450 activity. The decrease in electron density interferes with the formation of a high-valent oxo-ferryl species called Compound I. However, steady-state turnover rates of styrene epoxidation with H2O2 show the following trend: WT approximately equal to A320Q < F310A approximately equal to F310A/A320Q. The shunt pathway which can provide the two electrons and oxygen required for a P450 reaction instead of NAD(P)H and dioxygen can rule out the first and second heme reduction in the catalytic process. Because the electron density on the heme iron might be deeply involved in the k cat values in this system, the intermediate Compound 0 which is the precursor species of Compound I mainly appears to participate dominantly in epoxidation with H2O2.

Crystal structure of inhibitor-bound P450BM-3 reveals open conformation of substrate access channel

P450BM-3 is an extensively studied P450 cytochrome that is naturally fused to a cytochrome P450 reductase domain. Crystal structures of the heme domain of this enzyme have previously generated many insights into features of P450 structure, substrate binding specificity, and conformational changes that occur on substrate binding. Although many P450s are inhibited by imidazole, this compound does not effectively inhibit P450BM-3. Omega-imidazolyl fatty acids have previously been found to be weak inhibitors of the enzyme and show some unusual cooperativity with the substrate lauric acid. We set out to improve the properties of these inhibitors by attaching the omega-imidazolyl fatty acid to the nitrogen of an amino acid group, a tactic that we used previously to increase the potency of substrates. The resulting inhibitors were significantly more potent than their parent compounds lacking the amino acid group. A crystal structure of one of the new inhibitors bound to the heme domain of P450BM-3 reveals that the mode of interaction of the amino acid group with the enzyme is different from that previously observed for acyl amino acid substrates. Further, required movements of residues in the active site to accommodate the imidazole group provide an explanation for the low affinity of imidazole itself. Finally, the previously observed cooperativity with lauric acid is explained by a surprisingly open substrate-access channel lined with hydrophobic residues that could potentially accommodate lauric acid in addition to the inhibitor itself.

High-yield expression and purification of isotopically labeled cytochrome P450 monooxygenases for solid-state NMR spectroscopy

Cytochrome P450 monooxygenases (P450s), which represent the major group of drug metabolizing enzymes in humans, also catalyze important synthetic and detoxicative reactions in insects, plants and many microbes. Flexibilities in their catalytic sites and membrane associations are thought to play central roles in substrate binding and catalytic specificity. To date, Escherichia coli expression strategies for structural analysis of eukaryotic membrane-bound P450s by X-ray crystallography have necessitated full or partial removal of their N-terminal signal anchor domain and, often, replacement of residues more peripherally associated with the membrane (such as the F-G loop region). Even with these modifications, investigations of P450 structural flexibility remain challenging with multiple single crystal conditions needed to identify spatial variations between substrate-free and different substrate-bound forms. To overcome these limitations, we have developed methods for the efficient expression of 13C- and 15N-labeled P450s and analysis of their structures by magic-angle spinning solid-state NMR (SSNMR) spectroscopy. In the presence of co-expressed GroEL and GroES chaperones, full-length (53 kDa) Arabidopsis 13C,15N-labeled His4CYP98A3 is expressed at yields of 2-4 mg per liter of minimal media without the necessity of generating side chain modifications or N-terminal deletions. Precipitated His4CYP98A3 generates high quality SSNMR spectra consistent with a homogeneous, folded protein. These data highlight the potential of these methodologies to contribute to the structural analysis of membrane-bound proteins.

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

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

Crystal structures and catalytic mechanism of cytochrome P450 StaP that produces the indolocarbazole skeleton

Staurosporine isolated from Streptomyces sp. TP-A0274 is a member of the family of indolocarbazole alkaloids that exhibit strong antitumor activity. A key step in staurosporine biosynthesis is the formation of the indolocarbazole core by intramolecular C-C bond formation and oxidative decarboxylation of chromopyrrolic acid (CPA) catalyzed by cytochrome P450 StaP (StaP, CYP245A1). In this study, we report x-ray crystal structures of CPA-bound and -free forms of StaP. Upon substrate binding, StaP adopts a more ordered conformation, and conformational rearrangements of residues in the active site are also observed. Hydrogen-bonding interactions of two carboxyl groups and T-shaped pi-pi interactions with indole rings hold the substrate in the substrate-binding cavity with a conformation perpendicular to the heme plane. Based on the crystal structure of StaP-CPA complex, we propose that C-C bond formation occurs through an indole cation radical intermediate that is equivalent to cytochrome c peroxidase compound I [Sivaraja M, Goodin DB, Smith M, Hoffman BM (1989) Science 245:738-740]. The subsequent oxidative decarboxylation reaction is also discussed based on the crystal structure. Our crystallographic study shows the first crystal structures of enzymes involved in formation of the indolocarbazole core and provides valuable insights into the process of staurosporine biosynthesis, combinatorial biosynthesis of indolocarbazoles, and the diversity of cytochrome P450 chemistry.

Cytochrome P450 systems—biological variations of electron transport chains

Cytochromes P450 (P450) are hemoproteins encoded by a superfamily of genes nearly ubiquitously distributed in different organisms from all biological kingdoms. The reactions carried out by P450s are extremely diverse and contribute to the biotransformation of drugs, the bioconversion of xenobiotics, the bioactivation of chemical carcinogens, the biosynthesis of physiologically important compounds such as steroids, fatty acids, eicosanoids, fat-soluble vitamins and bile acids, the conversion of alkanes, terpenes and aromatic compounds as well as the degradation of herbicides and insecticides. Cytochromes P450 belong to the group of external monooxygenases and thus receive the necessary electrons for oxygen cleavage and substrate hydroxylation from different redox partners. The classical as well as the recently discovered P450 redox systems are compiled in this paper and classified according to their composition.

Protein dynamics and imidazole binding in cytochrome P450 enzymes

P450 (cytochrome P450) enzymes have major roles in the biosynthesis of endogenous factors such as steroids and eicosanoids, in the termination of the action of endogenous factors such as retinoic acid, in the metabolism of most drugs and xenobiotics and in the generation of toxic and carcinogenic products. Understanding the determinants of the substrate and inhibitor specificities of these enzymes is important for drug design. The crystallographic analysis of the deformability of two bacterial P450 active sites associated with the binding of azole (a class of inhibitors with an imidazole or triazole ring that co-ordinates to the haem iron) inhibitors described in the present study illustrates the importance of protein conformational malleability in the binding of imidazole derivatives.

The structural basis for substrate anchoring, active site selectivity, and product formation by P450 PikC from Streptomyces venezuelae

The pikromycin (Pik)/methymycin biosynthetic pathway of Streptomyces venezuelae represents a valuable system for dissecting the fundamental mechanisms of modular polyketide biosynthesis, aminodeoxysugar assembly, glycosyltransfer, and hydroxylation leading to the production of a series of macrolide antibiotics, including the natural ketolides narbomycin and pikromycin. In this study, we describe four x-ray crystal structures and allied functional studies for PikC, the remarkable P450 monooxygenase responsible for production of a number of related macrolide products from the Pik pathway. The results provide important new insights into the structural basis for the C10/C12 and C12/C14 hydroxylation patterns for the 12-(YC-17) and 14-membered ring (narbomycin) macrolides, respectively. This includes two different ligand-free structures in an asymmetric unit (resolution 2.1 A) and two co-crystal structures with bound endogenous substrates YC-17 (resolution 2.35 A)or narbomycin (resolution 1.7 A). A central feature of the enzyme-substrate interaction involves anchoring of the desosamine residue in two alternative binding pockets based on a series of distinct amino acid residues that form a salt bridge and a hydrogen-bonding network with the deoxysugar C3' dimethylamino group. Functional significance of the salt bridge was corroborated by site-directed mutagenesis that revealed a key role for Glu-94 in YC-17 binding and Glu-85 for narbomycin binding. Taken together, the x-ray structure analysis, site-directed mutagenesis, and corresponding product distribution studies reveal that PikC substrate tolerance and product diversity result from a combination of alternative anchoring modes rather than an induced fit mechanism.

Designing better drugs: predicting cytochrome P450 metabolism

Many 3D ligand-based and structure-based computational approaches have been used to predict, and thus help explain, the metabolism catalyzed by the enzymes of the cytochrome P450 superfamily (P450s). P450s are responsible for >90% of the metabolism of all drugs, so the computational prediction of metabolism can help to design out drug-drug interactions in the early phases of the drug discovery process. Computational methodologies have focused on a few P450s that are directly involved in drug metabolism. The recently derived crystal structures for human P450s enable better 3D modelling of these important metabolizing enzymes. Models derived for P450s have evolved from simple comparisons of known substrates to more-elaborate experiments that require considerable computer power involving 3D overlaps and docking experiments. These models help to explain and, more importantly, predict the involvement of P450s in the metabolism of specific compounds and guide the drug-design process.

Functional interaction of nitrogenous organic bases with cytochrome P450: A critical assessment and update of substrate features and predicted key active-site elements steering the access, binding, and orientation of amines

The widespread use of nitrogenous organic bases as environmental chemicals, food additives, and clinically important drugs necessitates precise knowledge about the molecular principles governing biotransformation of this category of substrates. In this regard, analysis of the topological background of complex formation between amines and P450s, acting as major catalysts in C- and N-oxidative attack, is of paramount importance. Thus, progress in collaborative investigations, combining physico-chemical techniques with chemical-modification as well as genetic engineering experiments, enables substantiation of hypothetical work resulting from the design of pharmacophores or homology modelling of P450s. Based on a general, CYP2D6-related construct, the majority of prospective amine-docking residues was found to cluster near the distal heme face in the six known SRSs, made up by the highly variant helices B', F and G as well as the N-terminal portion of helix C and certain beta-structures. Most of the contact sites examined show a frequency of conservation < 20%, hinting at the requirement of some degree of conformational versatility, while a limited number of amino acids exhibiting a higher level of conservation reside close to the heme core. Some key determinants may have a dual role in amine binding and/or maintenance of protein integrity. Importantly, a series of non-SRS elements are likely to be operative via long-range effects. While hydrophobic mechanisms appear to dominate orientation of the nitrogenous compounds toward the iron-oxene species, polar residues seem to foster binding events through H-bonding or salt-bridge formation. Careful uncovering of structure-function relationships in amine-enzyme association together with recently developed unsupervised machine learning approaches will be helpful in both tailoring of novel amine-type drugs and early elimination of potentially toxic or mutagenic candidates. Also, chimeragenesis might serve in the construction of more efficient P450s for activation of amine drugs and/or bioremediation.

Conformational flexibility of mammalian cytochrome P450 2B4 in binding imidazole inhibitors with different ring chemistry and side chains. Solution thermodynamics and molecular modeling

Recent x-ray structures of cytochrome P450 2B4 (CYP2B4) reveal an open form that undergoes a large-scale structural transition to a closed form upon binding to 4-(4-chlorophenyl)imidazole (4-CPI). Here, we report for the first time a complete solution thermodynamic study using isothermal titration calorimetry supported by spectroscopic studies to elucidate the conformational flexibility of CYP2B4 in binding imidazole inhibitors with different ring chemistry and side chains: 4-CPI, 1-benzylimidazole (1-BI), 1-CPI, 4-phenylimidazole (4-PI), 1-(2-(benzyloxy)ethyl)imidazole (BEI), and 1-PI. Each of the inhibitors induced type II spectral changes, and IC50 values for enzyme inhibition ranged from 0.1 to 2.4 microM, following the order 1-BI < 4-CPI < 1-CPI < 4-PI < BEI < 1-PI. Calorimetric titrations using monomeric enzyme yielded a 1:1 binding stoichiometry, with the associated KD values ranging from 0.3 to 4.8 microM and following the same rank order as the IC50 values. Changes in enthalpy at 25 degrees C ranged from -6.5 to -8.8 kcal mol(-1). The largest difference in binding entropy (+5.9 versus -4.1 cal mol(-1) K(-1)) was observed between 4-CPI and BEI, respectively, with a 2-fold difference in heat capacity changes (-604 versus -331 cal mol(-1) K(-1)), which is inferred to result from the reduction of apolar surface area of the enzyme ensuing from a conformational change upon 4-CPI binding. Accessibility to acrylamide of the only tryptophan (Trp121), which is located in helix C, was greatly decreased only in protein bound to 4-CPI. Steric restrictions hindered the perfect docking of only BEI to the closed conformation of the enzyme. The thermodynamic signature obtained for structurally similar inhibitors suggests remarkable plasticity of CYP2B4.

Structural biology of heme monooxygenases

Over the past few years the number of crystal structures available for heme monooxygenases has substantially increased. Those most closely related to one another are cytochrome P450, nitric oxide synthase, and heme oxygenase. The present mini-review provides a summary of some recently published work on how crystallography and solution studies have provided new insights on function and especially the oxygen activation process. It now appears that in all three monooxygenases highly ordered solvent in the active site serves as direct proton donors to the iron-linked dioxygen; a requirement for splitting the O-O bond. This is in sharp contrast to the related peroxidase family of enzymes where strategically positioned amino acid side chains serve the function of shuttling protons. The P450cam-oxy-complex as well as various mutants in a complex with either oxygen or carbon monoxide have enabled a fairly detailed picture to be developed on the role of specific amino acids and conformational changes in both electron transfer and oxygen activation.

Thermophilic cytochrome P450 enzymes

Thermophilic cytochrome P450 enzymes are of potential interest from structural, mechanistic, and biotechnological points of view. The structures and properties of two such enzymes, CYP119 and CYP175A1, have been investigated and provide the foundation for future work on thermophilic P450 enzymes.

Cytochrome P450 biosensors—a review

Cytochrome P450 (CYP) is a large family of enzymes containing heme as the active site. Since their discovery and the elucidation of their structure, they have attracted the interest of scientist for many years, particularly due to their catalytic abilities. Since the late 1970s attempts have concentrated on the construction and development of electrochemical sensors. Although sensors based on mediated electron transfer have also been constructed, the direct electron transfer approach has attracted most of the interest. This has enabled the investigation of the electrochemical properties of the various isoforms of CYP. Furthermore, CYP utilized to construct biosensors for the determination of substrates important in environmental monitoring, pharmaceutical industry and clinical practice.

Structure and direct electrochemistry of cytochrome P450 from the thermoacidophilic crenarchaeon, Sulfolobus tokodaii strain 7

Cytochrome P450 from thermoacidophilic crenarchaeon, Sulfolobus tokodaii strain 7 (P450st) has been expressed in Escherichia coli and purified at high homogeneity. P450st was crystallized in an orthorhombic system with the space group P2(1)2(1)2(1) and cell dimensions of a=53.6 A, b=55.1 A, and c=130.9 A, and the structure was determined at a 3.0 A resolution. The final R-factor was 0.194 (Rfree=0.235). Structural comparison with cytochrome P450 from S. solfataricus (CYP119) suggests that the region composed of the F to G helices and the Cl- binding site is responsible for the affinity for a ligand coordinating heme iron. Direct electrochemistry of P450st in a didodecyldimethylammonium bromide (DDAB) film on a plastic formed carbon (PFC) electrode has also been demonstrated. A quasi-reversible redox response has been observed even at elevated temperatures of up to 80 degrees C.

The Sulfolobus solfataricus electron donor partners of thermophilic CYP119: an unusual non-NAD(P)H-dependent cytochrome P450 system

CYP119 from Sulfolobus solfataricus is the first well-characterized thermophilic cytochrome P450 enzyme. The endogenous substrate for this enzyme is not known but it hydroxylates lauric acid in a reaction supported by surrogate mesophilic electron donors. However, reconstitution of a high-temperature catalytic system requires identification of the normal thermophilic electron donor partners of CYP119. Here, we describe cloning, expression in Escherichia coli, and characterization of the requisite electron donor partners from S. solfataricus. One is a thermostable ferredoxin and the second a 2-oxoacid-ferredoxin oxidoreductase that utilizes pyruvic acid rather than NAD(P)H as the source of reducing equivalents. CYP119 is the only cytochrome P450 to date known to obtain electrons from a non-NAD(P)H-dependent protein. The two thermophilic partners have been used to reconstitute a catalytic system that hydroxylates lauric acid at 70 degrees C, and the optimal conditions for this system have been defined. This first high-temperature in vitro catalytic system represents an important step in the development of industrially relevant catalysts.

Structural and functional diversity in heme monooxygenases

Recent advances in understanding structure-function relationships in cytochrome P450 (P450), nitric-oxide synthase (NOS), and heme oxygenase are summarized. Of particular importance is the role that dynamics plays in P450 function, where the active site undergoes large open/close motions to enable substrates to bind and products to leave. In sharp contrast, the heme-containing active site of NOS is rigid and remains relatively exposed compared with P450s. This difference in dynamics and active site exposure requires that the O(2) activation machinery operate somewhat differently in P450 and NOS. Owing to the open NOS active site, the NOS-oxy complex could be subject to nonspecific protonation that short-circuits the normal reaction path. One working hypothesis holds that NOS recruited the cofactor, tetrahydrobiopterin, to bind near the heme for very rapid coupled electron/proton transfer to the oxy complex, which avoids indiscriminate reaction with bulk solvent. Despite these differences, P450, NOS, and also heme oxygenase use a very similar network of H-bonded water molecules in the active site that are required for oxygen activation. Both P450 and NOS are important drug targets. With NOS, the structural basis for isoform-selective inhibition by a class of dipeptide inhibitors has been worked out, thus providing the basis for structure-based drug design.