Structure conservation in cytochromes P450

Multilevel superposition for deciphering the conformational variability of protein ensembles

The dynamics and variability of protein conformations are directly linked to their functions. Many comparative studies of X-ray protein structures have been conducted to elucidate the relevant conformational changes, dynamics and heterogeneity. The rapid increase in the number of experimentally determined structures has made comparison an effective tool for investigating protein structures. For example, it is now possible to compare structural ensembles formed by enzyme species, variants or the type of ligands bound to them. In this study, the author developed a multilevel model for estimating two covariance matrices that represent inter- and intra-ensemble variability in the Cartesian coordinate space. Principal component analysis using the two estimated covariance matrices identified the inter-/intra-enzyme variabilities, which seemed to be important for the enzyme functions, with the illustrative examples of cytochrome P450 family 2 enzymes and class A $\beta$-lactamases. In P450, in which each enzyme has its own active site of a distinct size, an active-site motion shared universally between the enzymes was captured as the first principal mode of the intra-enzyme covariance matrix. In this case, the method was useful for understanding the conformational variability after adjusting for the differences between enzyme sizes. The developed method is advantageous in small ensemble-size problems and hence promising for use in comparative studies on experimentally determined structures where ensemble sizes are smaller than those generated, for example, by molecular dynamics simulations.

Rational Design of Taxadiene Hydroxylase by Ancestral Enzyme Construction and the Elucidation of Key Amino Acids

Cytochrome P450 monooxygenases (CYP450s) play an important role in the biosynthesis of natural products by activating inert C-H bonds and inserting hydroxyl groups. However, the activities of most plant-derived CYP450s are extremely low, limiting the heterologous biosynthesis of natural products. Traditional enzyme engineering methods, either rational or screening-based, are not suitable for CYP450s because of the lack of crystal structures and high-throughput screening methods for this class of enzymes. CYP725A4 is the first hydroxylase involved in the biosynthesis pathway of Taxol. Its low activity, promiscuity, and multispecificity make it a bottleneck in Taxol biosynthesis. Here, we identified key amino acids that affect the in vivo activity of CYP725A4 by constructing the ancestral enzymes of CYP725A4. We obtained positive mutants that showed an improved yield of hydroxylated products based on the key amino acids identified, providing guidance for the modification of other CYP450s involved in the biosynthesis of natural products.

Gene excavation and expression analysis of CYP and UGT related to the post modifying stage of gypenoside biosynthesis in Gynostemma pentaphyllum (Thunb.) Makino by comprehensive analysis of RNA and proteome sequencing

Previous studies have revealed that gypenosides produced from Gynostemma pentaphyllum (Thunb.) Makino are mainly dammarane-type triterpenoid saponins with diverse structures and important biological activities, but the mechanism of diversity for gypenoside biosynthesis is still unclear. In this study, a combination of isobaric tags for relative and absolute quantification (iTRAQ) proteome analysis and RNA sequencing transcriptome analysis was performed to identify the proteins and genes related to gypenoside biosynthesis. A total of 3925 proteins were identified by proteomic sequencing, of which 2537 were quantified. Seventeen cytochrome P450 (CYP) and 11 uridine 5’-diphospho-glucuronosyltransferase (UDP-glucuronosyltransferase, UGT) candidate genes involved in the side chain synthesis and modification of gypenosides were found. Seven putative CYPs (CYP71B19, CYP77A3, CYP86A7, CYP86A8, CYP89A2, CYP90A1, CYP94A1) and five putative UGTs (UGT73B4, UGT76B1, UGT74F2, UGT91C1 and UGT91A1) were selected as candidate structural modifiers of triterpenoid saponins, which were cloned for gene expression analysis. Comprehensive analysis of RNA sequencing and proteome sequencing showed that some CYPs and UGTs were found at both the transcription and translation levels. In this study, an expression analysis of 7 CYPs and 5 UGTs that contributed to gypenoside biosynthesis and distribution in G. pentaphyllum was performed, providing consistent results that will inspire more future research on vital genes/proteins involved in gypenoside biosynthesis.

19F-NMR reveals substrate specificity of CYP121A1 in Mycobacterium tuberculosis

Cytochromes P450 are versatile enzymes that function in endobiotic and xenobiotic metabolism and undergo meaningful structural changes that relate to their function. However, the way in which conformational changes inform the specific recognition of the substrate is often unknown. Here, we demonstrate the utility of fluorine (¹⁹F)-NMR spectroscopy to monitor structural changes in CYP121A1, an essential enzyme from Mycobacterium tuberculosis. CYP121A1 forms functional dimers that catalyze the phenol-coupling reaction of the dipeptide dicyclotyrosine. The thiol-reactive compound 3-bromo-1,1,1-trifluoroacetone was used to label an S171C mutation of the enzyme FG loop, which is located adjacent to the homodimer interface. Substrate titrations and inhibitor-bound ¹⁹F-NMR spectra indicate that ligand binding reduces conformational heterogeneity at the FG loop in both the dimer and in an engineered monomer of CYP121A1. However, only the dimer was found to promote a substrate-bound conformation that was preexisting in the substrate-free spectra, thus confirming a role for the dimer interface in dicyclotyrosine recognition. Moreover, ¹⁹F-NMR spectra in the presence of substrate analogs indicate the hydrogen-bonding feature of the dipeptide aromatic side chain as a dicyclotyrosine specificity criterion. This study demonstrates the utility of ¹⁹F-NMR as applied to a multimeric cytochrome P450, while also revealing mechanistic insights for an essential M. tuberculosis enzyme.

Engineering Sequence and Selectivity of Late-Stage C-H Oxidation in the MycG Iterative Cytochrome P450

MycG is a multifunctional P450 monooxygenase that catalyzes sequential hydroxylation and epoxidation or a single epoxidation in mycinamicin biosynthesis. In the mycinamicin-producing strain Micromonospora griseorubida A11725, very low-level accumulation of mycinamicin V generated by the initial C-14 allylic hydroxylation of MycG is observed due to its subsequent epoxidation to generate mycinamicin II, the terminal metabolite in this pathway. Herein, we investigated whether MycG can be engineered for production of the mycinamicin II intermediate as the predominant metabolite. Thus, mycG was subject to random mutagenesis and screening was conducted in Escherichia coli whole-cell assays. This enabled efficient identification of amino acid residues involved in reaction profile alterations, which included MycG R111Q/V358L, W44R, and V135G/E355K with enhanced monohydroxylation to accumulate mycinamicin V. The MycG V135G/E355K mutant generated 40-fold higher levels of mycinamicin V compared to wild-type M. griseorubida A11725. In addition, the E355K mutation showed improved ability to catalyze sequential hydroxylation and epoxidation with minimal mono-epoxidation product mycinamicin I compared to the wild-type enzyme. These approaches demonstrate the ability to selectively coordinate the catalytic activity of multifunctional P450s and efficiently produce the desired compounds.

Single Nucleotide Polymorphism Induces Divergent Dynamic Patterns in CYP3A5: A Microsecond Scale Biomolecular Simulation of Variants Identified in Sub-Saharan African Populations

Pharmacogenomics aims to reveal variants associated with drug response phenotypes. Genes whose roles involve the absorption, distribution, metabolism, and excretion of drugs, are highly polymorphic between populations. High coverage whole genome sequencing showed that a large proportion of the variants for these genes are rare in African populations. This study investigated the impact of such variants on protein structure to assess their functional importance. We used genetic data of CYP3A5 from 458 individuals from sub-Saharan Africa to conduct a structural bioinformatics analysis. Five missense variants were modeled and microsecond scale molecular dynamics simulations were conducted for each, as well as for the CYP3A5 wildtype and the Y53C variant, which has a known deleterious impact on enzyme activity. The binding of ritonavir and artemether to CYP3A5 variant structures was also evaluated. Our results showed different conformational characteristics between all the variants. No significant structural changes were noticed. However, the genetic variability seemed to act on the plasticity of the protein. The impact on drug binding might be drug dependant. We concluded that rare variants hold relevance in determining the pharmacogenomics properties of populations. This could have a significant impact on precision medicine applications in sub-Saharan Africa.

Cytochrome P450 expression and regulation in the brain

The regulation of brain cytochrome P450 enzymes (CYPs) is different compared with respective hepatic enzymes. This may result from anatomical bases and physiological functions of the two organs. The brain is composed of a variety of functional structures built of different interconnected cell types endowed with specific receptors that receive various neuronal signals from other brain regions. Those signals activate transcription factors or alter functioning of enzyme proteins. Moreover, the blood-brain barrier (BBB) does not allow free penetration of all substances from the periphery into the brain. Differences in neurotransmitter signaling, availability to endogenous and exogenous active substances, and levels of transcription factors between neuronal and hepatic cells lead to differentiated expression and susceptibility to the regulation of CYP genes in the brain and liver. Herein, we briefly describe the CYP enzymes of CYP1-3 families, their distribution in the brain, and discuss brain-specific regulation of CYP genes. In parallel, a comparison to liver CYP regulation is presented. CYP enzymes play an essential role in maintaining the levels of bioactive molecules within normal ranges. These enzymes modulate the metabolism of endogenous neurochemicals, such as neurosteroids, dopamine, serotonin, melatonin, anandamide, and exogenous substances, including psychotropics, drugs of abuse, neurotoxins, and carcinogens. The role of these enzymes is not restricted to xenobiotic-induced neurotoxicity, but they are also involved in brain physiology. Therefore, it is crucial to recognize the function and regulation of CYP enzymes in the brain to build a foundation for future medicine and neuroprotection and for personalized treatment of brain diseases.

New insights into diosgenin biosynthesis pathway and its regulation in Trigonella foenum-graecum L

INTRODUCTION

Throughout history, thousands of medicinal and aromatic plants have been widely utilised by people worldwide. Owing to them possessing of valuable compounds with little side effects in comparison with chemical drugs, herbs have been of interest to humans for a number of purposes. Diosgenin, driven from fenugreek, Trigonella foenum-graecum L., has extensively drawn scientist's attention owing to having curable properties and being a precursor of steroid hormones synthesis. Nonetheless, complete knowledge about the biosynthesis pathway of this metabolite is still elusive.

OBJECTIVE

In the present research, we isolated the full-length CDS of 14 genes involving in diosgenin formation and measured their expression rate in various genotypes, which had illustrated different amount of diosgenin.

METHODOLOGY

The genes were successfully isolated, and functional motifs were also assessed using in silico approaches.

RESULTS

Moreover, combining transcript and metabolite analysis revealed that there are many genes playing the role in diosgenin formation, some of which are highly influential. Among them, ∆²⁴ -reductase, which converts cycloartenol to cycloartanol, is the first-committed and rate-limiting enzyme in this pathway. Additionally, no transcripts indicating to the presence or expression of lanosterol synthase were detected, contradicting the previous hypothesis about the biosynthetic pathway of diosgenin in fenugreek.

CONCLUSION

Considering all these, therefore, we propose the most possible pathway of diosgenin. This knowledge will then pave the way toward cloning the genes as well as engineering the diosgenin biosynthesis pathway.

Post-Translational Regulation of CYP450s Metabolism As Revealed by All-Atoms Simulations of the Aromatase Enzyme

Phosphorylation by kinases enzymes is a widespread regulatory mechanism able of rapidly altering the function of target proteins. Among these are cytochrome P450s (CYP450), a superfamily of enzymes performing the oxidation of endogenous and exogenous substrates thanks to the electron supply of a redox partner. In spite of its pivotal role, the molecular mechanism by which phosphorylation modulates CYP450s metabolism remains elusive. Here by performing microsecond-long all-atom molecular dynamics simulations, we disclose how phosphorylation regulates estrogen biosynthesis, catalyzed by the Human Aromatase (HA) enzyme. Namely, we unprecedentedly propose that HA phosphorylation at Y361 markedly stabilizes its adduct with the flavin mononucleotide domain of CYP450s reductase (CPR), the redox partner of microsomal CYP450s, and a variety of other proteins. With CPR present at physiological conditions in a limiting ratio with respect to its multiple oxidative partners, the enhanced stability of the CPR/HA adduct may favor HA in the competition with the other proteins requiring CPR's electron supply, ultimately facilitating the electron transfer and estrogen biosynthesis. As a result, our work elucidates at atomic-level the post-translational regulation of CYP450s catalysis. Given the potential for rational clinical management of diseases associated with steroid metabolism disorders, unraveling this mechanism is of utmost importance, and raises the intriguing perspective of exploiting this knowledge to devise novel therapies.

Characterization of the cytochrome P450 monooxygenase genes (P450ome) from the carotenogenic yeast Xanthophyllomyces dendrorhous

BACKGROUND

The cytochromes P450 (P450s) are a large superfamily of heme-containing monooxygenases involved in the oxidative metabolism of an enormous diversity of substrates. These enzymes require electrons for their activity, and the electrons are supplied by NAD(P)H through a P450 electron donor system, which is generally a cytochrome P450 reductase (CPR). The yeast Xanthophyllomyces dendrorhous has evolved an exclusive P450-CPR system that specializes in the synthesis of astaxanthin, a carotenoid with commercial potential. For this reason, the aim of this work was to identify and characterize other potential P450 genes in the genome of this yeast using a bioinformatic approach.

RESULTS

Thirteen potential P450-encoding genes were identified, and the analysis of their deduced proteins allowed them to be classified in ten different families: CYP51, CYP61, CYP5139 (with three members), CYP549A, CYP5491, CYP5492 (with two members), CYP5493, CYP53, CYP5494 and CYP5495. Structural analyses of the X. dendrorhous P450 proteins showed that all of them have a predicted transmembrane region at their N-terminus and have the conserved domains characteristic of the P450s, including the heme-binding region (FxxGxRxCxG); the PER domain, with the characteristic signature for fungi (PxRW); the ExxR motif in the K-helix region and the oxygen-binding domain (OBD) (AGxDTT); also, the characteristic secondary structure elements of all the P450 proteins were identified. The possible functions of these P450s include primary, secondary and xenobiotic metabolism reactions such as sterol biosynthesis, carotenoid synthesis and aromatic compound degradation.

CONCLUSIONS

The carotenogenic yeast X. dendrorhous has thirteen P450-encoding genes having potential functions in primary, secondary and xenobiotic metabolism reactions, including some genes of great interest for fatty acid hydroxylation and aromatic compound degradation. These findings established a basis for future studies about the role of P450s in the carotenogenic yeast X. dendrorhous and their potential biotechnological applications.

Reactivity Patterns of (Protonated) Compound II and Compound I of Cytochrome P450: Which is the Better Oxidant?

The cytochromes P450 are versatile enzymes in human physiology that perform substrate hydroxylation reactions extremely efficiently. In this work, we present results of a computational study on the reactivity patterns of Compound I, Compound II, and protonated Compound II with model substrates, and we address the question of which of these compounds is the most effective oxidant? All calculations, regardless of the substrate, implicated that Compound I is the superior oxidant of the three. However, Compound II and protonated Compound II were found to react with free energies of activation that are only a few kcal mol^-1 higher in energy than those obtained with Compound I. Therefore, Compound II and protonated Compound II should be able to react with aliphatic groups with moderate C-H bond strengths. We have analysed all results in detail and have given electronic, thermochemical, valence bond, and molecular orbital rationalizations on the reactivity differences and explained experimental product distributions. Overall, the findings implied that alternative oxidants could operate alongside Compound I in complex reaction mechanisms of enzymatic and synthetic iron porphyrinoid complexes.

Recent progress on structural bioinformatics research of cytochrome P450 and its impact on drug discovery

Cytochrome P450 is predominantly responsible for human drug metabolism, which is of critical importance for drug discovery and development. Structural bioinformatics focuses on analysis and prediction of three-dimentional structure of biological macromolecules and elucidation of structure-function relationship as well as identification of important binding interactions. Rapid advancement of structural bioinformatics has been made over the last decade. With more information available for CYP structures, the methods of structural bioinformatics may be used in the CYP field. In this review, we demonstrate three previous studies on CYP using the methods of structural bioinformatics, including the investigation of reasons for decrease of enzymatic activity of CYP1A2 caused by a peripheral mutation, the construction of a pharmacophore model specific to active site of CYP1A2 and the prediction of the functional consequences of single residue mutation in CYP. By illustrating these studies we attempt to show the potential role of structural bioinformatics in CYP research and help better understanding the importance of structural bioinformatics in drug designing.

Structural and functional characterization of a cytochrome P450 2B4 F429H mutant with an axial thiolate-histidine hydrogen bond

The structural basis of the regulation of microsomal cytochrome P450 (P450) activity was investigated by mutating the highly conserved heme binding motif residue, Phe429, on the proximal side of cytochrome P450 2B4 to a histidine. Spectroscopic, pre-steady-state and steady-state kinetic, thermodynamic, theoretical, and structural studies of the mutant demonstrate that formation of an H-bond between His429 and the unbonded electron pair of the Cys436 axial thiolate significantly alters the properties of the enzyme. The mutant lost >90% of its activity; its redox potential was increased by 87 mV, and the half-life of the oxyferrous mutant was increased ∼37-fold. Single-crystal electronic absorption and resonance Raman spectroscopy demonstrated that the mutant was reduced by a small dose of X-ray photons. The structure revealed that the δN atom of His429 forms an H-bond with the axial Cys436 thiolate whereas the εN atom forms an H-bond with the solvent and the side chain of Gln357. The amide of Gly438 forms the only other H-bond to the tetrahedral thiolate. Theoretical quantification of the histidine-thiolate interaction demonstrates a significant electron withdrawing effect on the heme iron. Comparisons of structures of class I-IV P450s demonstrate that either a phenylalanine or tryptophan is often found at the location corresponding to Phe429. Depending on the structure of the distal pocket heme, the residue at this location may or may not regulate the thermodynamic properties of the P450. Regardless, this residue appears to protect the thiolate from solvent, oxidation, protonations, and other deleterious reactions.

MutaCYP: Classification of missense mutations in human cytochromes P450

BACKGROUND

Cytochrome P450 monooxygenases (CYPs) represent a large and diverse family of enzymes involved in various biological processes in humans. Individual genome sequencing has revealed multiple mutations in human CYPs, and many missense mutations have been associated with variety of diseases. Since 3D structures are not resolved for most human CYPs, there is a need for a reliable sequence-based prediction that discriminates benign and disease causing mutations.

METHODS

A new prediction method (MutaCYP) has been developed for scoring de novo missense mutations to have a deleterious effect. The method utilizes only five features, all of which are sequence-based: predicted relative solvent accessibility (RSA), variance of predicted RSA among the residues in close sequence proximity, Z-score of Shannon entropy for a given position, difference in similarity scores and weighted difference in size between wild type and new amino acids. The method is based on a single neural network.

RESULTS

MutaCYP achieves MCC = 0.70, Q2 = 88.52%, Recall = 93.40% with Precision = 91.09%, and AUC = 0.909. Comparative evaluation with other existing methods indicates that MutaCYP outperforms SIFT and PolyPhen-2. Predictions by MutaCYP appear to be orthogonal to predictions by the evaluated methods. Potential issues on reliability of annotations of mutations in the existing databases are discussed.

CONCLUSIONS

A new accurate method, MutaCYP, for classification of missense mutations in human CYPs is presented. The prediction model consists of only five sequence-based features, including a real-valued predicted relative solvent accessibility. The method is publicly available at http://research.cchmc.org/MutaSense/.

Prediction of amino acid residues participated in substrate recognition by cytochrome P450 subfamilies with broad substrate specificity

Cytochromes P450 comprise a large superfamily and several of their isoforms play a crucial role in metabolism of xenobiotics, including drugs. Although these enzymes demonstrate broad and cross-substrate specificity, different cytochrome P450 subfamilies exhibit certain selectivity for some types of substrates. Analysis of amino acid residues of the active sites of six cytochrome subfamilies (CYP1А, CYP2А, CYP2С, CYP2D, CYP2E and CYP3А) enables to define subfamily-specific patterns that consist of four residues. These residues are located on the periphery of the active sites of these cytochromes. We suggest that they can form a primary binding site at the entrance to the active site, defining cytochrome substrate recognition.

CYP1B1, MYOC, and LTBP2 mutations in primary congenital glaucoma patients in the United States

PURPOSE

To screen primary congenital glaucoma patients in the United States for sequence variants within the CYP1B1, LTBP2, and MYOC genes using Sanger and whole exome sequencing.

DESIGN

Retrospective case-control study.

METHODS

Fifty-seven primary congenital glaucoma patients (47 families), 71 unaffected family members of the primary congenital glaucoma probands, and 101 healthy unrelated individuals were recruited from a single institution. Sanger sequencing of the primary congenital glaucoma gene, CYP1B1, was performed on 47 proband deoxyribonucleic acid samples. Simultaneously, whole exome sequencing was conducted on 3 families, each including more than 1 affected individual. Concurrently, 33 of 47 primary congenital glaucoma probands with extended family deoxyribonucleic acid samples were screened for LTBP2 and MYOC gene mutations. Exome-sequenced variations were validated by additional Sanger sequencing to confirm segregation of filtered disease-causing single nucleotide variations.

RESULTS

Seven primary congenital glaucoma families (14.9%) manifested disease phenotypes attributable to CYP1B1 mutations. One primary congenital glaucoma family possessed homozygous mutant alleles, whereas 6 families carried compound heterozygous mutations. Five novel combinations of compound heterozygous mutations were identified, of which 2 combinations were found with whole exome sequencing. No disease-causing mutations within the LTBP2 and MYOC genes were discovered.

CONCLUSIONS

This study analyzed CYP1B1, LTBP2, and MYOC mutations in a cohort of primary congenital glaucoma patients from the United States, applying whole exome sequencing as a complementary tool to Sanger sequencing. Whole exome sequencing, coupled with Sanger sequencing, may identify novel genes in primary congenital glaucoma patients who have no mutations in known primary congenital glaucoma genes.

Origins of P450 diversity

The P450 enzymes maintain a conserved P450 fold despite a considerable variation in sequence. The P450 family even includes proteins that lack the single conserved cysteine and are therefore no longer haem-thiolate proteins. The mechanisms of successive gene duplications leading to large families in plants and animals are well established. Comparisons of P450 CYP gene clusters in related species illustrate the rapid changes in CYPome sizes. Examples of CYP copy number variation with effects on fitness are emerging, and these provide an opportunity to study the proximal causes of duplication or pseudogenization. Birth and death models can explain the proliferation of CYP genes that is amply illustrated by the sequence of every new genome. Thus, the distribution of P450 diversity within the CYPome of plants and animals, a few families with many genes (P450 blooms) and many families with few genes, follows similar power laws in both groups. A closer look at some families with few genes shows that these, often single member families, are not stable during evolution. The enzymatic prowess of P450 may predispose them to switch back and forth between metabolism of critical structural or signal molecules and metabolism dedicated to environmental response.

Mass spectrometric analyses of microsomal cytochrome P450 isozymes isolated from β-naphthoflavone-treated Atlantic cod (Gadus morhua) liver reveal insights into the cod CYPome

Four cytochrome P450 (CYP) isozymes were purified earlier from liver microsomes of Atlantic cod (Gadus morhua) exposed to β-naphthoflavone. These isozymes were named P-450a-d and described with regard to optical properties, enzymatic activity, molecular weight and immunological reactivity. Subsequent analyses of the individual CYP fractions have shown that P-450c corresponds to a CYP1A isozyme, and immunochemical analyses have indicated that P-450b belongs to the CYP3A subfamily. However, no sequence data has been obtained to confirm these results, and the identities of the P-450a and P-450d have also remained unknown. The sequencing effort of the cod genome and expanding cod EST-databases (www.codgenome.no) have substantially increased the number of protein sequences available from cod. Mass spectrometric techniques were therefore applied to further characterize the proteins in historically archived samples of P-450a-d fractions. These analyses revealed large heterogeneities of CYPs within the purified samples. The most prominent CYP isozymes present include members of the CYP1A, CYP1C, CYP3A and CYP4V families. In total, 29 unique CYPs belonging to 9 CYP gene families and 15 subfamilies were identified with mass spectrometry. This analysis was also accompanied with genomic mining as the first step to unveil the full suite of cod CYP genes. In total, 55 CYP genes were predicted from the cod genome, distributed among 16 CYP gene families that are also present in other fish as well as mammals. Importantly, the majority of the CYPs revealed in this study have not previously been reported from cod, and represents the first proteomic survey to uncover the expressed complement of CYPs in any non-mammalian species.

Molecular characterizations of two cytochrome P450 genes encoding CYP6A41 and CYP6EK1 from the oriental fruit fly, Bactrocera dorsalis (Diptera: Tephritidae)

Two P450 genes encoding CYP6A41 and CYP6EK1 were cloned from the oriental fruit fly using polymerase chain reaction (PCR) and rapid amplification of cDNA ends (RACE) techniques. CYP6A41 and CYP6EK1 contained open reading frames of 1,530 and 1,524 nucleotides that encode 510 and 508 amino acid residues, respectively. The putative proteins shared 44% identity with each other. Phylogenetic analysis showed that CYP6A41 and CYP6EK1 were most closely related to Ceratitis capitata CYP6A10 and CYP6A subfamily. Expression patterns of the two genes in different geographical populations (Yunnan, Hainan, Dongguang, and Guangzhou), developmental stages (eggs, larvae, pupae, and adults), and tissues (midguts, fat bodies, and Malpighian tubules) were analyzed by real-time quantitative PCR (RT-qPCR) methods. The results showed that the expression levels of CYP6EK1 were significantly different among the four populations, but were not different for CYP6A41. Both the expressions of CYP6A41 and CYP6EK1 were development specific and had significantly higher levels in the larval stage. The expression of CYP6A41 did not vary among the midgut, fat body, or Malpighian tubules; however, CYP6EK1 expression was higher in the Malpighian tubules. The results suggest that CYP6A41 and CYP6EK1 might be involved in detoxification of xenobiotic compounds that were harmful to larval flies or development. Moreover, high expression of CYP6EK1 in the Malpighian tubules also implied participation in detoxification.

A fungal P450 (CYP5136A3) capable of oxidizing polycyclic aromatic hydrocarbons and endocrine disrupting alkylphenols: role of Trp(129) and Leu(324)

The model white rot fungus Phanerochaete chrysosporium, which is known for its versatile pollutant-biodegradation ability, possesses an extraordinarily large repertoire of P450 monooxygenases in its genome. However, the majority of these P450s have hitherto unknown function. Our initial studies using a genome-wide gene induction strategy revealed multiple P450s responsive to individual classes of xenobiotics. Here we report functional characterization of a cytochrome P450 monooxygenase, CYP5136A3 that showed common responsiveness and catalytic versatility towards endocrine-disrupting alkylphenols (APs) and mutagenic/carcinogenic polycyclic aromatic hydrocarbons (PAHs). Using recombinant CYP5136A3, we demonstrated its oxidation activity towards APs with varying alkyl side-chain length (C3-C9), in addition to PAHs (3–4 ring size). AP oxidation involves hydroxylation at the terminal carbon of the alkyl side-chain (ω-oxidation). Structure-activity analysis based on a 3D model indicated a potential role of Trp¹²⁹ and Leu³²⁴ in the oxidation mechanism of CYP5136A3. Replacing Trp¹²⁹ with Leu (W129L) and Phe (W129F) significantly diminished oxidation of both PAHs and APs. The W129L mutation caused greater reduction in phenanthrene oxidation (80%) as compared to W129F which caused greater reduction in pyrene oxidation (88%). Almost complete loss of oxidation of C3-C8 APs (83–90%) was observed for the W129L mutation as compared to W129F (28–41%). However, the two mutations showed a comparable loss (60–67%) in C9-AP oxidation. Replacement of Leu³²⁴ with Gly (L324G) caused 42% and 54% decrease in oxidation activity towards phenanthrene and pyrene, respectively. This mutation also caused loss of activity towards C3-C8 APs (20–58%), and complete loss of activity toward nonylphenol (C9-AP). Collectively, the results suggest that Trp¹²⁹ and Leu³²⁴ are critical in substrate recognition and/or regio-selective oxidation of PAHs and APs. To our knowledge, this is the first report on an AP-oxidizing P450 from fungi and on structure-activity relationship of a eukaryotic P450 for fused-ring PAHs (phenanthrene and pyrene) and AP substrates.

Prediction and analysis of the modular structure of cytochrome P450 monooxygenases

Background

Cytochrome P450 monooxygenases (CYPs) form a vast and diverse family of highly variable sequences. They catalyze a wide variety of oxidative reactions and are therefore of great relevance in drug development and biotechnological applications. Despite their differences in sequence and substrate specificity, the structures of CYPs are highly similar. Although being in research focus for years, factors mediating selectivity and activity remain vague.

Description

This systematic comparison of CYPs based on the Cytochrome P450 Engineering Database (CYPED) involved sequence and structure analysis of more than 8000 sequences. 31 structures have been applied to generate a reliable structure-based HMM profile in order to predict structurally conserved regions. Therefore, it was possible to automatically transfer these modules on CYP sequences without any secondary structure information, to analyze substrate interacting residues and to compare interaction sites with redox partners.

Conclusions

Functionally relevant structural sites of CYPs were predicted. Regions involved in substrate binding were analyzed in all sequences among the CYPED. For all CYPs that require a reductase, two reductase interaction sites were identified and classified according to their length. The newly gained insights promise an improvement of engineered enzyme properties for potential biotechnological application. The annotated sequences are accessible on the current version of the CYPED. The prediction tool can be applied to any CYP sequence via the web interface at http://www.cyped.uni-stuttgart.de/cgi-bin/strpred/dosecpred.pl.

Plasticity of cytochrome P450 2B4 as investigated by hydrogen-deuterium exchange mass spectrometry and X-ray crystallography

Crystal structures of the xenobiotic metabolizing cytochrome P450 2B4 have demonstrated markedly different conformations in the presence of imidazole inhibitors or in the absence of ligand. However, knowledge of the plasticity of the enzyme in solution has remained scant. Thus, hydrogen-deuterium exchange mass spectrometry (DXMS) was utilized to probe the conformations of ligand-free P450 2B4 and the complex with 4-(4-chlorophenyl)imidazole (4-CPI) or 1-biphenyl-4-methyl-1H-imidazole (1-PBI). The results of DXMS indicate that the binding of 4-CPI slowed the hydrogen-deuterium exchange rate over the B'- and C-helices and portions of the F-G-helix cassette compared with P450 2B4 in the absence of ligands. In contrast, there was little difference between the ligand-free and 1-PBI-bound exchange sets. In addition, DXMS suggests that the ligand-free P450 2B4 is predominantly open in solution. Interestingly, a new high resolution structure of ligand-free P450 2B4 was obtained in a closed conformation very similar to the 4-CPI complex. Molecular dynamics simulations performed with the closed ligand-free structure as the starting point were used to probe the energetically accessible conformations of P450 2B4. The simulations were found to equilibrate to a conformation resembling the 1-PBI-bound P450 2B4 crystal structure. The results indicate that conformational changes observed in available crystal structures of the promiscuous xenobiotic metabolizing cytochrome P450 2B4 are consistent with its solution structural behavior.

Five novel mutations in CYP11B2 gene detected in patients with aldosterone synthase deficiency type I: Functional characterization and structural analyses

CONTEXT

Aldosterone synthase deficiency (ASD) is an important differential diagnosis of diseases associated with salt wasting in early infancy.

OBJECTIVE

The objective of this study was to investigate the molecular basis for the disorder by (1) molecular genetic analysis in the CYP11B2 from patients suffering from ASD type I. (2) Functional characterization of the missense mutant gene products. (3) Structural simulation of the missense mutations.

RESULTS

Patient 1 was a homozygous carrier of a novel mutation located in exon 4 causing a premature stop codon (p.W260X). Patient 2 was analyzed to be compound heterozygous for two novel mutations: The first was an insertion mutation (p.G206WfsX51), and the second was a deletion mutation (p.L496SfsX169). Two siblings (patients 3 and 4) were compound heterozygous carriers of two novel missense mutations (p.S315R, p.R374W). The expression studies of the mutant proteins in COS-1 cells showed a complete absence of CYP11B2 activity of p.S315R and p.R374W mutants for the conversion of 11-deoxycorticosterone to aldosterone. A 3-D model of CYP11B2 p.S315R and p.R374W indicated a change of the hydrogen bond network which might explain the cause of the dysfunction.

CONCLUSION

We have identified the first CYP11B2 gene defects in two Polish families associated with phenotypes of ASD type I. Analysis of the enzymatic function as a complementary procedure to genotyping revealed data for understanding the clinical phenotype of ASD. Molecular modeling of the mutated enzyme provided a rational basis for understanding the changed activities of the mutant proteins.

Theoretical characterization of substrate access/exit channels in the human cytochrome P450 3A4 enzyme: involvement of phenylalanine residues in the gating mechanism

The human cytochrome P450 3A4 mono-oxygenates ∼50% of all drugs. Its substrates/products enter/leave the active site by access/exit channels. Here, we perform steered molecular dynamics simulations, pulling the products temazepam and testosterone-6βOH out of the P450 3A4 enzyme in order to identify the preferred substrate/product pathways and their gating mechanism. We locate six different egress pathways of products from the active site with different exit preferences for the two products and find that there is more than just one access/exit channel in CYP3A4. The so-called solvent channel manifests the largest opening for both tested products, thereby identifying this channel as a putative substrate channel. Most channels consist of one or two π-stacked phenylalanine residues that serve as gate keepers. The oxidized drug breaks the hydrophobic interactions of the gating residues and forms mainly hydrophobic contacts with the gate. We argue that product exit preferences in P450s are regulated by protein−substrate specificity.

Two-dimensional NMR and all-atom molecular dynamics of cytochrome P450 CYP119 reveal hidden conformational substates

Cytochrome P450 enzymes are versatile catalysts involved in a wide variety of biological processes from hormonal regulation and antibiotic synthesis to drug metabolism. A hallmark of their versatility is their promiscuous nature, allowing them to recognize a wide variety of chemically diverse substrates. However, the molecular details of this promiscuity have remained elusive. Here, we have utilized two-dimensional heteronuclear single quantum coherence NMR spectroscopy to examine a series of mutants site-specific labeled with the unnatural amino acid, [(13)C]p-methoxyphenylalanine, in conjunction with all-atom molecular dynamics simulations to examine substrate and inhibitor binding to CYP119, a P450 from Sulfolobus acidocaldarius. The results suggest that tight binding hydrophobic ligands tend to lock the enzyme into a single conformational substate, whereas weak binding low affinity ligands bind loosely in the active site, resulting in a distribution of localized conformers. Furthermore, the molecular dynamics simulations suggest that the ligand-free enzyme samples ligand-bound conformations of the enzyme and, therefore, that ligand binding may proceed largely through a process of conformational selection rather than induced fit.

New insights into the structural characteristics and functional relevance of the human cytochrome P450 2D6 enzyme

To date, the crystal structures of at least 12 human CYPs (1A2, 2A6, 2A13, 2C8, 2C9, 2D6, 2E1, 2R1, 3A4, 7A1, 8A1, and 46A1) have been determined. CYP2D6 accounts for only a small percentage of all hepatic CYPs (< 2%), but it metabolizes approximately 25% of clinically used drugs with significant polymorphisms. CYP2D6 also metabolizes procarcinogens and neurotoxins, such as 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine, 1,2,3,4-tetrahydroquinoline, and indolealkylamines. Moreover, the enzyme utilizes hydroxytryptamines and neurosteroids as endogenous substrates. Typical CYP2D6 substrates are usually lipophilic bases with an aromatic ring and a nitrogen atom, which can be protonated at physiological pH. Substrate binding is generally followed by oxidation (5-7 A) from the proposed nitrogen-Asp301 interaction. A number of homology models have been constructed to explore the structural features of CYP2D6, while antibody studies also provide useful structural information. Site-directed mutagenesis studies have demonstrated that Glu216, Asp301, Phe120, Phe481, and Phe483 play important roles in determining the binding of ligands to CYP2D6. The structure of human CYP2D6 has been recently determined and shows the characteristic CYP fold observed for other members of the CYP superfamily. The lengths and orientations of the individual secondary structural elements in the CYP2D6 structure are similar to those seen in other human CYP2 members, such as CYP2C9 and 2C8. The 2D6 structure has a well-defined active-site cavity located above the heme group with a volume of approximately 540 A(3), which is larger than equivalent cavities in CYP2A6 (260 A(3)), 1A2 (375 A(3)), and 2E1 (190 A(3)), but smaller than those in CYP3A4 (1385 A(3)) and 2C8 (1438 A(3)). Further studies are required to delineate the molecular mechanisms involved in CYP2D6 ligand interactions and their implications for drug development and clinical practice.

The cytochrome P450 engineering database: Integration of biochemical properties

Background

Cytochrome P450 monooxygenases (CYPs) form a vast and diverse enzyme class of particular interest in drug development and a high biotechnological potential. Although very diverse in sequence, they share a common structural fold. For the comprehensive and systematic comparison of protein sequences and structures the Cytochrome P450 Engineering Database (CYPED) was established. It was built up based on an extensible data model that enables its functions readily enhanced.

Description

The new version of the CYPED contains information on sequences and structures of 8613 and 47 proteins, respectively, which strictly follow Nelson's classification rules for homologous families and superfamilies. To gain biochemical information on substrates and inhibitors, the CYPED was linked to the Cytochrome P450 Knowledgebase (CPK). To overcome differences in the data model and inconsistencies in the content of CYPED and CPK, a metric was established based on sequence similarity to link protein sequences as primary keys. In addition, the annotation of structurally and functionally relevant residues was extended by a reliable prediction of conserved secondary structure elements and by information on the effect of single nucleotide polymorphisms.

Conclusion

The online accessible version of the CYPED at http://www.cyped.uni-stuttgart.de provides a valuable tool for the analysis of sequences, structures and their relationships to biochemical properties.

Identification of selectivity-determining residues in cytochrome P450 monooxygenases: a systematic analysis of the substrate recognition site 5

The large and diverse family of cytochrome P450 monooxygenases (CYPs) was systematically analyzed to identify selectivity- and specificity-determining residues in the substrate recognition site 5, which is located in close vicinity to the heme center. A positively charged heme-interacting residue was identified in the structures of 29 monooxygenases and in 97.7% of the 6379 CYP sequences investigated here. This heme-interacting residue restricts the conformation of the substrate recognition site 5 and is preferentially located at position 10 or 11 after the conserved ExxR motif (in 94.4% of the sequences), in 3.3% of the sequences at position 9 or 12. As a result, a classification by the position of the heme-interacting residue allows to predict residues that are closest to the heme center and restrict its accessibility. In 98.4% of all CYP sequences a preferentially hydrophobic residue is located at position 5 after the ExxR motif that is predicted to point close to the heme center. Replacing this residue by hydrophobic residues of different size has been shown to change substrate specificity and regioselectivity for CYPs of different superfamilies. Twenty-seven percent of all CYPs are predicted to contain a second selectivity-determining residue at position 9 after the ExxR motif that can be identified by the pattern EXXR-X(7)-{P}-x-P-[HKR].

Advances in the interpretation and prediction of CYP2E1 metabolism from a biochemical perspective

BACKGROUND

Cytochrome P450 2E1 (CYP2E1) plays a central role in the metabolism and metabolic activation of a large number of small, mostly xenobiotic compounds. These qualities distinguish CYP2E1 from traditional enzymes and pose significant challenges to understanding the role and consequences of CYP2E1-mediated metabolism.

OBJECTIVE

This review discusses recent advances in kinetic profiling, quantitative structure-activity relationships and structural studies that have furthered the development of tools to interpret and predict CYP2E1 metabolism.

METHODS

Analysis of kinetic profiles by specific mechanisms produces important parameters describing specificity, stoichiometry and metabolism of molecules. Quantitative structure-activity relationships reveal a more specific basis for molecular recognition by CYP2E1. The corresponding protein structures imparting these interactions are the focus of chemical modifications, site-directed mutagenesis and homology modeling studies.

RESULTS

Compilation of kinetic profiling for CYP2E1 substrates established the selectivity for small substrates, whose characteristics could be generalized in parameters for hydrophobicity and steric properties as determined by quantitative structure-activity relationships. The possibility of an effector site for monocyclic compounds added an interesting variable to these modeling efforts. Various structural studies identified important residues contributing to binding and catalysis as well as the volume and location of the active site relative to the heme moiety. Pressure and carbon monoxide-binding experiments also demonstrated the inherent conformational flexibility of CYP2E1 that may contribute to rate-limiting steps during catalytic turnover.

CONCLUSION

Although combinations of these approaches have reinforced important observations, more work is needed to verify findings and seek broader impacts for various interpretative and predictive tools.

Alpha7 nicotinic acetylcholine receptor agonists: prediction of their binding affinity through a molecular mechanics Poisson-Boltzmann surface area approach

A group of agonists for the alpha7 neuronal nicotinic acetylcholine receptors (nAChRs) was investigated, and their free energies of binding DeltaG(bind) were calculated by applying the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) approach. This method, based on molecular dynamics simulations of fully solvated protein-ligand complexes, allowed us to estimate the contribution of both polar and nonpolar terms as well as the entropy to the overall free energy of binding. The calculated results were in a good agreement with the experimentally determined DeltaG(bind) values, thereby pointing to the MM-PBSA protocol as a valuable computational tool for the rational design of specific agents targeting the neuronal alpha7 nAChR subtypes.

Computational prediction of drug binding and rationalisation of selectivity towards cytochromes P450

BACKGROUND

Early in-vitro consideration of metabolism and inhibition of cytochrome P450 has proven its merits over the last 15 years. Simultaneously, many computational drug-design methods have been developed, and are being applied to study the interactions between drug candidates and cytochrome P450 enzymes (P450s).

OBJECTIVE

This review discusses the recent advances of these methods and the implications that are specific for P450s.

METHODS

Mainly focusing on the prediction of binding affinity and ligand selectivity, we outline the applicability of the different methods to answer specific questions. Special emphasis is put on the different levels of theory that are being used in recent computational descriptions of ligand-P450 interactions.

CONCLUSION

P450s offer an additional challenge for computational methods, considering the ambiguities of the catalytic cycle and the significant flexibility of the active site. Different computational methods display different limitations, which is crucial to take into account when choosing the method appropriate to each application.

Human cytochromes P450 in the metabolism of drugs: new molecular models of enzyme-substrate interactions

The overall predictive ability of molecular modelling, as applied to the cytochrome P450 (CYP) system, is analysed in the light of current developments in a variety of techniques, including X-ray crystallography, molecular biology, enzyme kinetics, molecular mechanics and dynamics, in relation to its relevance to drug metabolism in humans. This review demonstrates that it is possible to generate realistic models for the major human CYPs, which metabolise xenobiotics that compare favourably with crystal structures, and thus may be used to derive substrate binding energies that agree closely with experimental K(m) values obtained from enzyme kinetics.