Role of protein-protein interactions in cytochrome P450-mediated drug metabolism and toxicity

Identification of Protein-Protein Interaction Associated Functions Based on Gene Ontology

Protein-protein interactions (PPIs) involve the physical or functional contact between two or more proteins. Generally, proteins that can interact with each other always have special relationships. Some previous studies have reported that gene ontology (GO) terms are related to the determination of PPIs, suggesting the special patterns on the GO terms of proteins in PPIs. In this study, we explored the special GO term patterns on human PPIs, trying to uncover the underlying functional mechanism of PPIs. The experimental validated human PPIs were retrieved from STRING database, which were termed as positive samples. Additionally, we randomly paired proteins occurring in positive samples, yielding lots of negative samples. A simple calculation was conducted to count the number of positive samples for each GO term pair, where proteins in samples were annotated by GO terms in the pair individually. The similar number for negative samples was also counted and further adjusted due to the great gap between the numbers of positive and negative samples. The difference of the above two numbers and the relative ratio compared with the number on positive samples were calculated. This ratio provided a precise evaluation of the occurrence of GO term pairs for positive samples and negative samples, indicating the latent GO term patterns for PPIs. Our analysis unveiled several nuclear biological processes, including gene transcription, cell proliferation, and nutrient metabolism, as key biological functions. Interactions between major proliferative or metabolic GO terms consistently correspond with significantly reported PPIs in recent literature.

Combined impact of hypoalbuminemia and pharmacogenomic variants on voriconazole trough concentration: data from a real-life clinical setting in the Chinese population

Voriconazole (VRC) displays highly variable pharmacokinetics impacting treatment efficacy and safety. To provide evidence for optimizing VRC therapy regimens, the authors set out to determine the factors impacting VRC steady-state trough concentration (Cmin) in patients with various albumin (Alb) level. A total of 275 blood samples of 120 patients and their clinical characteristics and genotypes of CYP2C19, CYP3A4, CYP3A5, CYP2C9, FMO3, ABCB1, POR, NR1I2 and NR1I3 were included in this study. Results of multivariate linear regression analysis demonstrated that C-reactive protein (CRP) and total bilirubin (T-Bil) were predictors of the VRC Cmin adjusted for dose in patients with hypoalbuminemia (Alb < 35 g/L) (R2 = 0.16, P < 0.001). Additionally, in patients with normal albumin level (Alb ≥ 35 g/L), it resulted in a significant model containing factors of the poor metabolizer (PM) CYP2C19 genotype and CRP level (R2 = 0.26, P < 0.001). Therefore, CRP and T-Bil levels ought to receive greater consideration than genetic factors in patients with hypoalbuminemia.

Pleiotropy of Progesterone Receptor Membrane Component 1 in Modulation of Cytochrome P450 Activity

Progesterone receptor membrane component 1 (PGRMC1) is one of few proteins that have been recently described as direct modulators of the activity of human cytochrome P450 enzymes (CYP)s. These enzymes form a superfamily of membrane-bound hemoproteins that metabolize a wide variety of physiological, dietary, environmental, and pharmacological compounds. Modulation of CYP activity impacts the detoxification of xenobiotics as well as endogenous pathways such as steroid and fatty acid metabolism, thus playing a central role in homeostasis. This review is focused on nine main topics that include the most relevant aspects of past and current PGRMC1 research, focusing on its role in CYP-mediated drug metabolism. Firstly, a general overview of the main aspects of xenobiotic metabolism is presented (I), followed by an overview of the role of the CYP enzymatic complex (IIa), a section on human disorders associated with defects in CYP enzyme complex activity (IIb), and a brief account of cytochrome b5 (cyt b5)'s effect on CYP activity (IIc). Subsequently, we present a background overview of the history of the molecular characterization of PGRMC1 (III), regarding its structure, expression, and intracellular location (IIIa), and its heme-binding capability and dimerization (IIIb). The next section reflects the different effects PGRMC1 may have on CYP activity (IV), presenting a description of studies on the direct effects on CYP activity (IVa), and a summary of pathways in which PGRMC1's involvement may indirectly affect CYP activity (IVb). The last section of the review is focused on the current challenges of research on the effect of PGRMC1 on CYP activity (V), presenting some future perspectives of research in the field (VI).

CYP2C9, CYP3A and CYP2C19 metabolize Δ9-tetrahydrocannabinol to multiple metabolites but metabolism is affected by human liver fatty acid binding protein (FABP1)

Δ9-tetrahydrocannabinol (THC) is the psychoactive constituent of cannabis. It is cleared predominantly via metabolism. Metabolism to 11-OH-THC by cytochrome P450 (CYP) 2C9 has been proposed as the main clearance pathway of THC, with the estimated fraction metabolized (fm) about 70%. The remaining clearance pathways are not well established, and it is unknown how THC is eliminated in individuals with reduced CYP2C9 activity. The goal of this study was to systematically identify the CYP enzymes contributing to THC clearance and characterize the metabolites formed. Further, this study aimed to characterize the impact of liver fatty acid binding protein (FABP1) on THC metabolism by human CYPs. THC was metabolized to at least four different metabolites including 11-OH-THC in human liver microsomes (HLMs) and with recombinant CYPs. 11-OH-THC was formed by recombinant CYP2C9 (Km,u = 0.77 nM, kcat = 12 min-1) and by recombinant CYP2C19 (Km,u = 2.2 nM, kcat = 14 min-1). The other three major metabolites were likely hydroxylations in the cyclohexenyl ring and were formed mainly by recombinant CYP3A4/5 (Km,u > 10 nM). HLM experiments confirmed the contributions of CYP2C9, CYP2C19 and CYP3A to THC metabolism. The presence of FABP1 and THC binding to FABP1 altered THC metabolism by recombinant CYPs and HLMs in an enzyme and metabolite specific manner. This suggests that FABP1 may interact with CYP enzymes and alter the fm by CYPs towards THC metabolism. In conclusion, this study is the first to systematically establish the metabolic profile of THC by human CYPs and characterize how FABP1 binding alters CYP mediated THC metabolism.

Population Pharmacokinetics of Voriconazole and Dose Optimization in Elderly Chinese Patients

Voriconazole is commonly recommended as a first-line therapy for invasive aspergillosis infections. Elderly patients are susceptible to infectious diseases owing to their decreased physical function and immune system. Our study aims to establish a population pharmacokinetics model for elderly patients receiving intravenous voriconazole, and to optimize dosing protocols through a simulated approach. An accurate fit to the concentration-time profile of voriconazole was achieved by employing a 1-compartment model featuring first-order elimination. The typical clearance rate of voriconazole was found to be 3.22 L/h, with a typical volume of distribution of 194 L. The covariate analysis revealed that albumin (ALB), gamma-glutamyl transpeptidase, and direct bilirubin had significant impacts on voriconazole clearance. Additionally, body weight was found to be associated with the volume of distribution. Individualized dosing regimens were recommended for different ALB levels based on population pharmacokinetics model prediction. The proposed dosing regimens could provide a rationale for dosage individualization, improve the clinical outcomes, and minimize drug-related toxicities.

Computational identification and analysis of deleterious non-synonymous single nucleotide polymorphisms (nsSNPs) in the human POR gene: a structural and functional impact

Cytochrome P450 oxidoreductase (POR) protein is essential for steroidogenesis, and POR gene mutations are frequently associated with P450 Oxidoreductase Deficiency (PORD), a disorder of hormone production. To our knowledge, no previous attempt has been made to identify and analyze the deleterious/pathogenic non-synonymous single nucleotide polymorphisms (nsSNPs) in the human POR gene through an extensive computational approach. Computational algorithms and tools were employed to identify, characterize, and validate the pathogenic SNPs associated with certain diseases. To begin with, all the high-confidence SNPs were collected, and their structural and functional impacts on the protein structures were explored. The results of various in silico analyses affirm that the A287P and R457H variants of POR could destabilize the interactions between the amino acids and the hydrogen bond networks, resulting in functional deviations of POR. The literature study further confirms that the pathogenic mutations (A287P and R457H) are associated with the onset of PORD. Molecular dynamics simulations (MDS) and essential dynamics (ED) studies characterized the structural consequences of prioritized deleterious mutations, representing the structural destabilization that might disrupt POR biological function. The identified deleterious mutations at the cofactor's binding domains might interfere with the essential interactions between the protein and cofactors, thus inhibiting POR catalytic activity. The consolidated insights from the computational analyses can be used to predict potential deleterious mutants and understand the disease's pathological basis and the molecular mechanism of drug metabolism for the application of personalized medication. HIGHLIGHTSNADPH cytochrome P450 oxidoreductase (POR) mutations are associated with a broad spectrum of human diseasesIdentified and analyzed the most deleterious nsSNPs of POR through the sequence and structure-based prediction toolsInvestigated the structural and functional impacts of the most significant mutations (A287P and R457H) associated with PORDMolecular dynamics and PCA-based FEL analysis were utilized to probe the mutation-induced structural alterations in PORCommunicated by Ramaswamy H. Sarma.

Production of Gibberellins via a Non-Natural Pathway Using Steviol as a Substrate

Gibberellins (GAs) are plant hormones widely used in agriculture. At present, GAs are produced by fermentation of the fungus Fusarium fujikuroi. However, fungal growth is too slow, resulting in slow fungal fermentation and a low yield. Here, to develop an alternative production source of GAs, an artificial pathway was engineered in Escherichia coli. By selecting and combining enzymes derived from plants and bacteria, a novel 4-enzyme pathway was successfully constructed to produce GAs using steviol, a readily available and less valuable byproduct during enzymatic refining of rebaudioside A, as a feedstock. Whole-cell biotransformation with E. coli strain expressing the novel pathway produced 71.17 ± 2.00 mg/L GA1 and a trace amount of GA3 from steviol in 48 h. This report presents a significant advancement in the fast production of GAs and establishes a method for the metabolism of terpenoids to produce target products in microbial hosts.

Combinatorial optimization and spatial remodeling of CYPs to control product profile

Activating inert substrates is a challenge in nature and synthetic chemistry, but essential for creating functionally active molecules. In this work, we used a combinatorial optimization approach to assemble cytochrome P450 monooxygenases (CYPs) and reductases (CPRs) to achieve a target product profile. By creating 110 CYP-CPR pairs and iteratively screening different pairing libraries, we demonstrated a framework for establishing a CYP network that catalyzes six oxidation reactions at three different positions of a chemical scaffold. Target product titer was improved by remodeling endoplasmic reticulum (ER) size and spatially controlling the CYPs' configuration on the ER. Out of 47 potential products that could be synthesized, 86% of the products synthesized by the optimized network was our target compound quillaic acid (QA), the aglycone backbone of many pharmaceutically important saponins, and fermentation achieved QA titer 2.23 g/L.

Impact of Heterotropic Allosteric Modulation on the Time-Dependent Inhibition of Cytochrome P450 3A4

The current study was designed to investigate the influence of allosteric effectors on the metabolism of the prototypical cytochrome P450 (CYP) 3A4 substrate midazolam (MDZ), and on the determination in vitro time-dependent inhibition (TDI) of CYP3A4 using human liver microsomes (HLM). As the concentration of midazolam increased to 250 µM in HLMs, homotropic cooperativity resulted in a decrease in the 1'-hydroxymidazolam to 4-hydroxymidazolam ratio to a maximum of 1.1. The presence of varying concentrations of testosterone, progesterone (PGS), or carbamazepine (CBZ) in HLMs with MDZ could recapitulate the effect of homotropic cooperativity such that the formation rates of the 1'hydroxymidazolam and 4-hydroxymidazolam were equal even at low concentrations of MDZ. The presence of PGS (10 or 100 µM) and CBZ (100 or 1000 µM) in in vitro TDI determination of four known CYP3A4 time-dependent inactivators (clarithromycin, troleandomycin, mibefradil, raloxifene) simultaneously decreased potency and inactivation rate constant, resulting in fold changes in inactivation efficiency on average of 1.6-fold and 13-fold for the low and high concentrations of allosteric modulator tested, respectively. The formation of a metabolic-intermediate complex (MIC) for clarithromycin and troleandomycin decreased in the presence of the allosteric modulators in a concentration-dependent manner, reaching a new steady state formation that could not be overcome with increased incubation time. Maximum reduction of the MIC formed by clarithromycin was up to ∼91%, while troleandomycin MIC decreased up to ∼31%. These findings suggest that the absence of endogenous allosteric modulators may contribute to the poor translation of HLM-based drug-drug interaction predictions. SIGNIFICANCE STATEMENT: The reported overprediction of in vitro human liver microsome time-dependent inhibition of CYP3A4 and observed drug interactions in vivo remains an issue in drug development. We provide characterization of allosteric modulators on the CYP3A4 metabolism of the prototypical substrate midazolam, demonstrating the ability of the modulators to recapitulate the homotropic cooperativity of midazolam. Furthermore, we demonstrate that allosteric heterotropic cooperativity of CYP3A4 can impact the time-dependent inhibition kinetics of known mechanisms-based inhibitors, providing a potential mechanism to explain the overprediction.

Using machine learning to predict protein-protein interactions between a zombie ant fungus and its carpenter ant host

Parasitic fungi produce proteins that modulate virulence, alter host physiology, and trigger host responses. These proteins, classified as a type of "effector," often act via protein-protein interactions (PPIs). The fungal parasite Ophiocordyceps camponoti-floridani (zombie ant fungus) manipulates Camponotus floridanus (carpenter ant) behavior to promote transmission. The most striking aspect of this behavioral change is a summit disease phenotype where infected hosts ascend and attach to an elevated position. Plausibly, interspecific PPIs drive aspects of Ophiocordyceps infection and host manipulation. Machine learning PPI predictions offer high-throughput methods to produce mechanistic hypotheses on how this behavioral manipulation occurs. Using D-SCRIPT to predict host-parasite PPIs, we found ca. 6000 interactions involving 2083 host proteins and 129 parasite proteins, which are encoded by genes upregulated during manipulated behavior. We identified multiple overrepresentations of functional annotations among these proteins. The strongest signals in the host highlighted neuromodulatory G-protein coupled receptors and oxidation-reduction processes. We also detected Camponotus structural and gene-regulatory proteins. In the parasite, we found enrichment of Ophiocordyceps proteases and frequent involvement of novel small secreted proteins with unknown functions. From these results, we provide new hypotheses on potential parasite effectors and host targets underlying zombie ant behavioral manipulation.

A frequent CYP2D6 variant promotes skipping of exon 3 and reduces CYP2D6 protein expression in human liver samples

CYP2D6 is one of the most polymorphic drug-metabolizing enzymes in the liver. While genetic CYP2D6 variants serve as clinical biomarkers to predict CYP2D6 activity, large inter-person variability in CYP2D6 expression remains unaccounted for. Previous results suggest that there is variable expression of a CYP2D6 splice isoform with an in-frame deletion of exon 3 (CYP2D6ΔE3) encoding a protein lacking numerous active site residues. Here, using fragment analysis and RT-qPCR, we revealed that rs1058164 G (MAF = 27%-43%) is associated with increased formation of CYP2D6∆E3 in human liver samples (1.4-2.5-fold) and transfected cells. Furthermore, western blots showed that rs1058164 G was associated with a 50% decrease in full-length hepatic CYP2D6 protein expression. In addition, by studying a larger liver cohort, we confirmed our previous results that rs16947 (CYP2D6*2) reduces full-length CYP2D6 mRNA by increasing the production of an unstable splice isoform lacking exon 6 (CYP2D6ΔE6) and that the impact of CYP2D6ΔE6 is offset in carriers of the downstream enhancer variant rs5758550. The three frequent SNPs (rs1058164, rs16947, and rs5758550) form various 3-SNP-haplotypes, each with distinct CYP2D6 expression characteristics. Using an expression score (ES) system, we tested the impact of the 3-SNP-haplotype on improving the standard model to predict hepatic CYP2D6 protein expression based on genotype. A model that incorporates the 3-SNP-haplotype provided the best fit for CYP2D6 expression and also accounted for more variability in CYP2D6 protein levels (59%) than a model based on the accepted standard (36%) or one that only adds rs16947 and rs5758550 (42%). Clinical studies are needed to determine whether including the 3-SNP-haplotype alongside current standard CYP2D6 models improves the predictive value of CYP2D6 panels.

Single Mutations in Cytochrome P450 Oxidoreductase Can Alter the Specificity of Human Cytochrome P450 1A2-Mediated Caffeine Metabolism

A unique cytochrome P450 (CYP) oxidoreductase (CPR) sustains activities of human microsomal CYPs. Its function requires toggling between a closed conformation enabling electron transfers from NADPH to FAD and then FMN cofactors and open conformations forming complexes and transferring electrons to CYPs. We previously demonstrated that distinct features of the hinge region linking the FAD and FMN domain (FD) modulate conformer poses and their interactions with CYPs. Specific FD residues contribute in a CYP isoform-dependent manner to the recognition and electron transfer mechanisms that are additionally modulated by the structure of CYP-bound substrate. To obtain insights into the underlying mechanisms, we analyzed how hinge region and FD mutations influence CYP1A2-mediated caffeine metabolism. Activities, metabolite profiles, regiospecificity and coupling efficiencies were evaluated in regard to the structural features and molecular dynamics of complexes bearing alternate substrate poses at the CYP active site. Studies reveal that FD variants not only modulate CYP activities but surprisingly the regiospecificity of reactions. Computational approaches evidenced that the considered mutations are generally in close contact with residues at the FD-CYP interface, exhibiting induced fits during complexation and modified dynamics depending on caffeine presence and orientation. It was concluded that dynamic coupling between FD mutations, the complex interface and CYP active site exist consistently with the observed regiospecific alterations.

Tissue-preferential recruitment of electron transfer chains for cytochrome P450-catalyzed phenolic biosynthesis

Cytochrome P450 system consists of P450 monooxygenase and redox pattern(s). While the importance of monooxygenases in plant metabolism is well documented, the metabolic roles of the related redox components have been largely overlooked. Here, we show that distinct electron transfer chains are recruited in phenylpropanoid-monolignol P450 systems to support the synthesis and distribution of different classes of phenolics in different plant tissues. While Arabidopsis cinnamate 4-hydroxylase adopts conventional NADPH-cytochrome P450 oxidoreductase (CPR) electron transfer chain for its para-hydroxylation reaction, ferulate 5-hydroxylase uses both NADPH-CPR-cytochrome b₅ (CB5) and NADH-cytochrome b₅ reductase-CB5 chains to support benzene ring 5-hydroxylation, in which the former route is primarily recruited in the stem for syringyl lignin synthesis, while the latter dominates in the syntheses of 5-hydroxylated phenolics in seeds and seed coat suberin. Our study unveils an additional layer of complexity and versatility of P450 system that the plants evolved for diversifying phenolic repertoires.

Cytochrome b 5: A versatile electron carrier and regulator for plant metabolism

Cytochrome b ₅ (CB5) is a small heme-binding protein, known as an electron donor delivering reducing power to the terminal enzymes involved in oxidative reactions. In plants, the CB5 protein family is substantially expanded both in its isoform numbers and cellular functions, compared to its yeast and mammalian counterparts. As an electron carrier, plant CB5 proteins function not only in fatty acid desaturation, hydroxylation and elongation, but also in the formation of specialized metabolites such as flavonoids, phenolic esters, and heteropolymer lignin. Furthermore, plant CB5s are found to interact with different non-catalytic proteins such as ethylene signaling regulator, cell death inhibitor, and sugar transporters, implicating their versatile regulatory roles in coordinating different metabolic and cellular processes, presumably in respect to the cellular redox status and/or carbon availability. Compared to the plentiful studies on biochemistry and cellular functions of mammalian CB5 proteins, the cellular and metabolic roles of plant CB5 proteins have received far less attention. This article summarizes the fragmentary information pertaining to the discovery of plant CB5 proteins, and discusses the conventional and peculiar functions that plant CB5s might play in different metabolic and cellular processes. Gaining comprehensive insight into the biological functions of CB5 proteins could offer effective biotechnological solutions to tailor plant chemodiversity and cellular responses to environment stimuli.

Biological Role of Pazopanib and Sunitinib Aldehyde Derivatives in Drug-Induced Liver Injury

Tyrosine kinase inhibitors pazopanib and sunitinib are both used to treat advanced renal cell carcinoma but expose patients to an increased risk of hepatotoxicity. We have previously identified two aldehyde derivatives for pazopanib and sunitinib (P-CHO and S-CHO, respectively) in liver microsomes. In this study, we aimed to decipher their role in hepatotoxicity by treating HepG2 and HepaRG hepatic cell lines with these derivatives and evaluating cell viability, mitochondrial dysfunction, and oxidative stress accumulation. Additionally, plasma concentrations of P-CHO were assessed in a cohort of patients treated with pazopanib. Results showed that S-CHO slightly decreased the viability of HepG2, but to a lesser extent than sunitinib, and affected the maximal respiratory capacity of the mitochondrial chain. P-CHO decreased viability and ATP production in HepG2. Traces of P-CHO were detected in the plasma of patients treated with pazopanib. Overall, these results showed that P-CHO and S-CHO affect hepatocyte integrity and could be involved in the pazopanib and sunitinib hepatotoxicity.

In vitro-in silico-based prediction of inter-individual and inter-ethnic variations in the dose-dependent cardiotoxicity of R- and S-methadone in humans

New approach methodologies predicting human cardiotoxicity are of interest to support or even replace in vivo-based drug safety testing. The present study presents an in vitro–in silico approach to predict the effect of inter-individual and inter-ethnic kinetic variations in the cardiotoxicity of R- and S-methadone in the Caucasian and the Chinese population. In vitro cardiotoxicity data, and metabolic data obtained from two approaches, using either individual human liver microsomes or recombinant cytochrome P450 enzymes (rCYPs), were integrated with physiologically based kinetic (PBK) models and Monte Carlo simulations to predict inter-individual and inter-ethnic variations in methadone-induced cardiotoxicity. Chemical specific adjustment factors were defined and used to derive dose–response curves for the sensitive individuals. Our simulations indicated that Chinese are more sensitive towards methadone-induced cardiotoxicity with Margin of Safety values being generally two-fold lower than those for Caucasians for both methadone enantiomers. Individual PBK models using microsomes and PBK models using rCYPs combined with Monte Carlo simulations predicted similar inter-individual and inter-ethnic variations in methadone-induced cardiotoxicity. The present study illustrates how inter-individual and inter-ethnic variations in cardiotoxicity can be predicted by combining in vitro toxicity and metabolic data, PBK modelling and Monte Carlo simulations. The novel methodology can be used to enhance cardiac safety evaluations and risk assessment of chemicals.

Exploring the Interactome of Cytochrome P450 2E1 in Human Liver Microsomes with Chemical Crosslinking Mass Spectrometry

Aiming to elucidate the system-wide effects of the alcohol-induced increase in the content of cytochrome P450 2E1 (CYP2E1) on drug metabolism, we explored the array of its protein-protein interactions (interactome) in human liver microsomes (HLM) with chemical crosslinking mass spectrometry (CXMS). Our strategy employs membrane incorporation of purified CYP2E1 modified with photoreactive crosslinkers benzophenone-4-maleimide and 4-(N-succinimidylcarboxy)benzophenone. Exposure of bait-incorporated HLM samples to light was followed by isolating the His-tagged bait protein and its crosslinked aggregates on Ni-NTA agarose. Analyzing the individual bands of SDS-PAGE slabs of thereby isolated protein with the toolset of untargeted proteomics, we detected the crosslinked dimeric and trimeric complexes of CYP2E1 with other drug-metabolizing enzymes. Among the most extensively crosslinked partners of CYP2E1 are the cytochromes P450 2A6, 2C8, 3A4, 4A11, and 4F2, UDP-glucuronosyltransferases (UGTs) 1A and 2B, fatty aldehyde dehydrogenase (ALDH3A2), epoxide hydrolase 1 (EPHX1), disulfide oxidase 1α (ERO1L), and ribophorin II (RPN2). These results demonstrate the exploratory power of the proposed CXMS strategy and corroborate the concept of tight functional integration in the human drug-metabolizing ensemble through protein-protein interactions of the constituting enzymes.

A Cytochrome B5-Like Heme/Steroid Binding Domain Protein, PlCB5L1, Regulates Mycelial Growth, Pathogenicity and Oxidative Stress Tolerance in Peronophythora litchii

As an electron transport component, cytochrome b₅ is an essential component of the Class II cytochrome P450 monooxygenation system and widely present in animals, plants, and fungi. However, the roles of Cyt-b₅ domain proteins in pathogenic oomycetes remain unknown. Peronophythora litchii is an oomycete pathogen that causes litchi downy blight, the most destructive disease of litchi. In this study, we identified a gene, designated PlCB5L1, that encodes a Cyt-b₅ domain protein in P. litchii, and characterized its function. PlCB5L1 is highly expressed in the zoospores, cysts, germinated cysts, and during early stages of infection. PlCB5L1 knockout mutants showed reduced growth rate and β-sitosterol utilization. Importantly, we also found that PlCB5L1 is required for the full pathogenicity of P. litchii. Compared with the wild-type strain, the PlCB5L1 mutants exhibited significantly higher tolerance to SDS and sorbitol, but impaired tolerance to cell wall stress, osmotic stress, and oxidative stress. Further, the expression of genes involved in oxidative stress tolerance, including peroxidase, cytochrome P450, and laccase genes, were down-regulated in PlCB5L1 mutants under oxidative stress. This is the first report that a Cyt-b₅ domain protein contributes to the development, stress response, and pathogenicity in plant pathogenic oomycetes.

Redox metabolism for improving whole-cell P450-catalysed terpenoid biosynthesis

The growing preference for producing cytochrome P450-mediated natural products in microbial systems stems from the challenging nature of the organic chemistry approaches. The P450 enzymes are redox-dependent proteins, through which they source electrons from reducing cofactors to drive their activities. Widely researched in biochemistry, most of the previous studies have extensively utilised expensive cell-free assays to reveal mechanistic insights into P450 functionalities in presence of commercial redox partners. However, in the context of microbial bioproduction, the synergic activity of P450- reductase proteins in microbial systems have not been largely investigated. This is mainly due to limited knowledge about their mutual interactions in the context of complex systems. Hence, manipulating the redox potential for natural product synthesis in microbial chassis has been limited. As the potential of redox state as crucial regulator of P450 biocatalysis has been greatly underestimated by the scientific community, in this review, we re-emphasize their pivotal role in modulating the in vivo P450 activity through affecting the product profile and yield. Particularly, we discuss the applications of widely used in vivo redox engineering methodologies for natural product synthesis to provide further suggestions for patterning on P450-based terpenoids production in microbial platforms.

The detoxification effect of cytochrome P450 3A4 on gelsemine-induced toxicity

Gelsemine (GA), the principal alkaloid in Gelsemium elegans Benth, exhibits potent and specific antinociception in chronic pain without the induction of apparent tolerance. However, GA also exerts neurotoxicity and hepatotoxicity when overdosed, and potential detoxification pathways are urgently needed. Cytochrome P450 enzymes (CYPs) are important phase I enzymes involved in the detoxification of xenobiotic compounds. The study aimed to investigate the role of CYPs-mediated metabolism in GA-induced toxicity. Microsomes, chemical special inhibitors and human recombinant CYPs indicated that GA was mainly metabolized by CYP3A4/5. The major metabolite of GA was isolated and identified as 4-N-demethyl-GA by high-resolution mass spectrometry and nuclear magnetic resonance technology. The CYP3A4 inhibitor ketoconazole significantly inhibited the metabolism of GA. This drastically increased GA toxicity which is caused by increasing the level of malondialdehyde and decreasing the level of the superoxide dismutase in mice. In contrast, the CYP3A4 inducer dexamethasone significantly increased GA metabolism and markedly decreased GA toxicity in mice. Notably, in CYP3A4-humanized mice, the toxicity of GA was significantly reduced compared to normal mice. These findings demonstrated that CYP3A4-mediated metabolism is a robust detoxification pathway for GA-induced toxicity.

Assembling the P450 puzzle: on the sources of nonadditivity in drug metabolism

There is an increasing number of indications of an oversimplification in the premise that the cumulative properties of the human drug-metabolizing ensemble represent a simple aggregate of the properties of the constituting enzymes. Recent studies of the functional effects of hetero-association of multiple cytochrome P450 species and their interactions with metabolically related enzymes revealed a tight integration in the drug-metabolizing ensemble. In our opinion, the sources of interindividual variability in drug metabolism can be elucidated only when considering this ensemble as a multienzyme system, the functional parameters of which are determined by interactions between its constituents. In this article, we present a conceptual model providing a mechanistic explanation for the functional effects of the interactions between multiple P450 species and propose a clue to understanding the nonadditive behavior of the drug-metabolizing ensemble.

Numerical Methods for Modeling Enzyme Kinetics

Differential equations are used to describe time-dependent changes in enzyme kinetics and pharmacokinetics. Analytical and numerical methods can be used to solve differential equations. This chapter describes the use of numerical methods in solving differential equations and its applications in characterizing the complexities observed in enzyme kinetics. A discussion is included on the use of numerical methods to overcome limitations of explicit equations in the analysis of metabolism kinetics, reversible inhibition kinetics, and inactivation kinetics. The chapter describes the advantages of using numerical methods when Michaelis-Menten assumptions do not hold.

Structural Dynamics of Cytochrome P450 3A4 in the Presence of Substrates and Cytochrome P450 Reductase

Cytochrome P450 3A4 (CYP3A4) is the most important drug-metabolizing enzyme in humans and has been associated with harmful drug interactions. The activity of CYP3A4 is known to be modulated by several compounds and by the electron transfer partner, cytochrome P450 reductase (CPR). The underlying mechanism of these effects, however, is poorly understood. We have used hydrogen-deuterium exchange mass spectrometry to investigate the impact of binding of CPR and of three different substrates (7-benzyloxy-4-trifluoromethyl-coumarin, testosterone, and progesterone) on the conformational dynamics of CYP3A4. Here, we report that interaction of CYP3A4 with substrates or with the oxidized or reduced forms of CPR leads to a global rigidification of the CYP3A4 structure. This was evident from the suppression of deuterium exchange in several regions of CYP3A4, including regions known to be involved in protein-protein interactions (helix C) and substrate binding and specificity (helices B' and E, and loop K/β1). Furthermore, the bimodal isotopic distributions observed for some CYP3A4-derived peptides were drastically impacted upon binding to CPR and/or substrates, suggesting the existence of stable CYP3A4 conformational populations that are perturbed by ligand/CPR binding. The results have implications for understanding the mechanisms of ligand binding, allostery, and catalysis in CYP enzymes.

Structure-based virtual screening of CYP1A1 inhibitors: towards rapid tier-one assessment of potential developmental toxicants

Cytochrome P450 1A1 (CYP1A1) metabolizes estrogens, melatonin, and other key endogenous signaling molecules critical for embryonic/fetal development. The enzyme has increasing expression during pregnancy, and its inhibition or knockout increases embryonic/fetal lethality and/or developmental problems. Here, we present a virtual screening model for CYP1A1 inhibitors based on the orthosteric and predicted allosteric sites of the enzyme. Using 1001 reference compounds with CYP1A1 activity data, we optimized the decision thresholds of our model and classified the training compounds with 68.3% balanced accuracy (91.0% sensitivity and 45.7% specificity). We applied our final model to 11 known CYP1A1 orthosteric binders and related compounds, and found that our ranking of the known orthosteric binders generally agrees with the relative activity of CYP1A1 in metabolizing these compounds. We also applied the model to 22 new test compounds with unknown/unclear CYP1A1 inhibitory activity, and predicted 16 of them are CYP1A1 inhibitors. The CYP1A1 potency and modes of inhibition of these 22 compounds were experimentally determined. We confirmed that most predicted inhibitors, including drugs contraindicated during pregnancy (amiodarone, bicalutamide, cyproterone acetate, ketoconazole, and tamoxifen) and environmental agents suspected to be endocrine disruptors (bisphenol A, diethyl and dibutyl phthalates, and zearalenone), are indeed potent inhibitors of CYP1A1. Our results suggest that virtual screening may be used as a rapid tier-one method to screen for potential CYP1A1 inhibitors, and flag them out for further experimental evaluations.

Mechanism of electron transfers mediated by cytochromes c and b5 in mitochondria and endoplasmic reticulum: classical and murburn perspectives

We explore the mechanism of electron transfers mediated by cytochrome c, a soluble protein involved in mitochondrial oxidative phosphorylation and cytochrome b₅, a microsomal membrane protein acting as a redox aide in xenobiotic metabolism. We found minimal conservation in the sequence and surface amino acid residues of cytochrome c/b₅ proteins among divergent species. Therefore, we question the evolutionary logic for electron transfer (ET) occurring through affinity binding via recognition of specific surface residues/topography. Also, analysis of putative protein-protein interactions in the crystal structures of these proteins and their redox partners did not point to any specific interaction logic. A comparison of the kinetic and thermodynamic constants of wildtype vs. mutants did not provide strong evidence to support the binding-based ET paradigm, but indicated support for diffusible reactive species (DRS)-mediated process. Topographically divergent cytochromes from one species have been substituted for reaction with proteins from other species, implying the involvement of non-specific interactions. We provide a viable alternative (murburn concept) to classical protein-protein binding-based long range ET mechanism. To account for the promiscuity of interactions and solvent-accessible hemes, we propose that the two proteins act as non- specific redox capacitors, mediating one-electron redox equilibriums involving DRS and unbound ions.Communicated by Ramaswamy H. Sarma.

Conformational Landscape of Cytochrome P450 Reductase Interactions

Oxidative reactions catalyzed by Cytochrome P450 enzymes (CYPs), which constitute the most relevant group of drug-metabolizing enzymes, are enabled by their redox partner Cytochrome P450 reductase (CPR). Both proteins are anchored to the membrane of the endoplasmic reticulum and the CPR undergoes a conformational change in order to interact with the respective CYP and transfer electrons. Here, we conducted over 22 microseconds of molecular dynamics (MD) simulations in combination with protein-protein docking to investigate the conformational changes necessary for the formation of the CPR-CYP complex. While some structural features of the CPR and the CPR-CYP2D6 complex that we highlighted confirmed previous observations, our simulations revealed additional mechanisms for the conformational transition of the CPR. Unbiased simulations exposed a movement of the whole protein relative to the  membrane, potentially to facilitate interactions with its diverse set of redox partners. Further, we present a structural mechanism for the susceptibility of the CPR to different redox states based on the flip of a glycine residue disrupting the local interaction network that maintains inter-domain proximity. Simulations of the CPR-CYP2D6 complex pointed toward an additional interaction surface of the FAD domain and the proximal side of CYP2D6. Altogether, this study provides novel structural insight into the mechanism of CPR-CYP interactions and underlying conformational changes, improving our understanding of this complex machinery Cytochrome P450 reductase; CPR; conformational; dynamicsrelevant for drug metabolism.

Bioinformatic analysis of potential hub genes in gastric adenocarcinoma

Gastric adenocarcinoma is the most common histologic type of gastric cancer; however, the pathogenic mechanisms remain unclear. To improve mechanistic understanding and identify new treatment targets or diagnostic biomarkers, we used bioinformatic tools to predict the hub genes related to the process of gastric adenocarcinoma development from public datasets, and explored their prognostic significance. We screened differentially expressed genes between gastric adenocarcinoma and normal gastric tissues in Gene Expression Omnibus datasets (GSE79973, GSE118916, and GSE29998) using the GEO2R tool, and their functions were annotated with Gene Ontology and Kyoto Encyclopedia of Genes and Genomes signaling pathway enrichment analyses in the DAVID database. Hub genes were identified based on the protein-protein network constructed in the STRING database with Cytoscape software. A total of 10 hub genes were selected for further analysis, and their expression patterns in gastric adenocarcinoma patients were investigated using the Oncomine GEPIA database. The expression levels of ATP4A, CA9, FGA, ALDH1A1, and GHRL were reduced, whereas those of TIMP1, SPP1, CXCL8, THY1, and COL1A1 were increased in gastric adenocarcinoma. The Kaplan-Meier online plotter tool showed associations of all hub genes except for CA9 with prognosis in gastric adenocarcinoma patients; CXCL8 and ALDH1A1 were positively correlated with survival, and the other genes were negatively correlated with survival. These 10 hub genes may be involved in important processes in gastric adenocarcinoma development, providing new directions for research to clarify the role of these genes and offer insight for improved treatment.

Population Pharmacokinetics of Voriconazole in Patients With Invasive Aspergillosis: Serum Albumin Level as a Novel Marker for Clearance and Dosage Optimization

BACKGROUND

Voriconazole (VRCZ) is an antifungal triazole recommended as an effective first-line agent for treating invasive aspergillosis.

OBJECTIVES

To develop a population pharmacokinetic model of VRCZ and trough concentration-based dosing simulation for dynamic patient conditions.

METHODS

The authors combined plasma VRCZ data from intensive sampling, and retrospective trough concentration monitoring for analysis. Nonlinear mixed-effects modeling with subsequent model validation was performed. The recommended dosage regimens were simulated based on the developed model.

RESULTS

The study participants included 106 patients taking oral VRCZ. A linear one-compartment model with first-order elimination and absorption best described the observed data. The CYP2C19 phenotypes did not influence the pharmacokinetic parameters. Serum albumin (SA) levels and gamma-glutamyl transferase significantly correlated with the VRCZ clearance rate, whereas the actual body weight influenced the volume. A visual predictive check showed good consistency with the observed data, whereas SA levels across the treatment course correlated with linear clearance, irrespective of the CYP2C19 phenotype. Patients with SA levels ≤30 g/L had lower linear clearance than that in patients with SA levels >30 g/L. Dosing simulation based on the developed model indicated that patients with SA levels of ≤30 g/L required a lower daily maintenance dose to attain the therapeutic trough level.

CONCLUSIONS

SA level was identified as a novel marker associated with VRCZ clearance. This marker may be a practical choice for physicians to perform therapeutic drug monitoring and optimize VRCZ dosage.

Interaction Modes of Microsomal Cytochrome P450s with Its Reductase and the Role of Substrate Binding

The activity of microsomal cytochromes P450 (CYP) is strictly dependent on the supply of electrons provided by NADPH cytochrome P450 oxidoreductase (CPR). The variant nature of the isoform-specific proximal interface of microsomal CYPs implies that the interacting interface between the two proteins is degenerated. Recently, we demonstrated that specific CPR mutations in the FMN-domain (FD) may induce a gain in activity for a specific CYP isoform. In the current report, we confirm the CYP isoform dependence of CPR's degenerated binding by demonstrating that the effect of four of the formerly studied FD mutants are indeed exclusive of a specific CYP isoform, as verified by cytochrome c inhibition studies. Moreover, the nature of CYP's substrate seems to have a modulating role in the CPR:CYP interaction. In silico molecular dynamics simulations of the FD evidence that mutations induces very subtle structural alterations, influencing the characteristics of residues formerly implicated in the CPR:CYP interaction or in positioning of the FMN moiety. CPR seems therefore to be able to form effective interaction complexes with its structural diverse partners via a combination of specific structural features of the FD, which are functional in a CYP isoform dependent manner, and dependent on the substrate bound.

The Role of the FMN-Domain of Human Cytochrome P450 Oxidoreductase in Its Promiscuous Interactions With Structurally Diverse Redox Partners

NADPH cytochrome P450 oxidoreductase (CPR) is the obligatory electron supplier that sustains the activity of microsomal cytochrome P450 (CYP) enzymes. The variant nature of the isoform-specific proximal interface of microsomal CYPs indicates that CPR is capable of multiple degenerated interactions with CYPs for electron transfer, through different binding mechanisms, and which are still not well-understood. Recently, we showed that CPR dynamics allows formation of open conformations that can be sampled by its structurally diverse redox partners in a CYP-isoform dependent manner. To further investigate the role of the CPR FMN-domain in effective binding of CPR to its diverse acceptors and to clarify the underlying molecular mechanisms, five different CPR-FMN-domain random mutant libraries were created. These libraries were screened for mutants with increased activity when combined with specific CYP-isoforms. Seven CPR-FMN-domain mutants were identified, supporting a gain in activity for CYP1A2 (P117H, G144C, A229T), 2A6 (P117L/L125V, G175D, H183Y), or 3A4 (N151D). Effects were evaluated using extended enzyme kinetic analysis, cytochrome b₅ competition, ionic strength effect on CYP activity, and structural analysis. Mutated residues were located either in or adjacent to several acidic amino acid stretches – formerly indicated to be involved in CPR:CYP interactions – or close to two tyrosine residues suggested to be involved in FMN binding. Several of the identified positions co-localize with mutations found in naturally occurring CPR variants that were previously shown to cause CYP-isoform-dependent effects. The mutations do not seem to significantly alter the geometry of the FMN-domain but are likely to cause very subtle alterations leading to improved interaction with a specific CYP. Overall, these data suggest that CYPs interact with CPR using an isoform specific combination of several binding motifs of the FMN-domain.

The Biosynthesis of Enzymatically Oxidized Lipids

Enzymatically oxidized lipids are a specific group of biomolecules that function as key signaling mediators and hormones, regulating various cellular and physiological processes from metabolism and cell death to inflammation and the immune response. They are broadly categorized as either polyunsaturated fatty acid (PUFA) containing (free acid oxygenated PUFA "oxylipins", endocannabinoids, oxidized phospholipids) or cholesterol derivatives (oxysterols, steroid hormones, and bile acids). Their biosynthesis is accomplished by families of enzymes that include lipoxygenases (LOX), cyclooxygenases (COX), cytochrome P450s (CYP), and aldo-keto reductases (AKR). In contrast, non-enzymatically oxidized lipids are produced by uncontrolled oxidation and are broadly considered to be harmful. Here, we provide an overview of the biochemistry and enzymology of LOXs, COXs, CYPs, and AKRs in humans. Next, we present biosynthetic pathways for oxylipins, oxidized phospholipids, oxysterols, bile acids and steroid hormones. Last, we address gaps in knowledge and suggest directions for future work.

Genome-Wide Identification and Analysis of the Cytochrome B5 Protein Family in Chinese Cabbage (Brassica rapa L. ssp. Pekinensis)

Cytochrome B5 (CB5) family proteins play an important role in various oxidation/reduction reactions in cells as the electron donor and are involved in a variety of biotic and abiotic stress processes. However, the function of the CB5s in Brassica rapa is still unclear. In this study, we carried out genome-wide identification, characterization, and expression analysis of BrCB5s in different tissues under adversities and stresses. It was identified that fifteen BrCB5s were distributed on different chromosomes, which were classified into seven groups (A-G) according to its phylogenetic relationship. Phylogenetic analysis of the CB5 protein sequences from six species showed that the BrCB5s conduct a close evolutionary process with the CB5s of Arabidopsis thaliana and far from those of Oryza sativa. Protein interaction analysis showed that 40 interaction patterns were predicted including two Sucrose Transporter 4 subfamily proteins (SUT 4) and Fatty Acid Hydroxylase 2 protein (FAH 2) can interact with most members of BrCB5s. The expression profile analysis indicated that BrCB5s were differentially expressed in different tissues, and the transcript abundances were significantly different under various abiotic stresses and plant hormone treatments. Our study provides a basis for a better understanding of the characteristics and biological functions of the CB5 family genes in Chinese cabbage during plant development, especially in plant responses to multiple stresses.

Transcriptomic Analysis Reveals the Roles of Detoxification Systems in Response to Mercury in Chromera velia

Heavy metal pollution is an increasing global concern. Among heavy metals, mercury (Hg) is especially dangerous because of its massive release into the environment and high toxicity, especially for aquatic organisms. The molecular response mechanisms of algae to Hg exposure are mostly unknown. Here, we combine physiological, biochemical, and transcriptomic analysis to provide, for the first time, a comprehensive view on the pathways activated in Chromera velia in response to toxic levels of Hg. Production of hydrogen peroxide and superoxide anion, two reactive oxygen species (ROS), showed opposite patterns in response to Hg²⁺ while reactive nitrogen species (RNS) levels did not change. A deep RNA sequencing analysis generated a total of 307,738,790 high-quality reads assembled in 122,874 transcripts, representing 89,853 unigenes successfully annotated in databases. Detailed analysis of the differently expressed genes corroborates the biochemical results observed in ROS production and suggests novel putative molecular mechanisms in the algal response to Hg²⁺. Moreover, we indicated that important transcription factor (TF) families associated with stress responses differentially expressed in C. velia cultures under Hg stress. Our study presents the first in-depth transcriptomic analysis of C. velia, focusing on the expression of genes involved in different detoxification defense systems in response to heavy metal stress.

Cytochrome b 5 Is an Obligate Electron Shuttle Protein for Syringyl Lignin Biosynthesis in Arabidopsis

Angiosperms have evolved the metabolic capacity to synthesize p-hydroxyphenyl, guaiacyl (G), and syringyl (S) lignin subunits in their cell walls to better adapt to the harsh terrestrial environment. The structural characteristics of lignin subunits are essentially determined by three cytochrome P450-catalzyed reactions. NADPH-dependent cytochrome P450 oxidoreductase (CPR) is commonly regarded as the electron carrier for P450-catalyzed reactions during monolignol biosynthesis. Here, we show that cytochrome b 5 isoform D (CB5D) is an indispensable electron shuttle protein specific for S-lignin biosynthesis. Arabidopsis (Arabidopsis thaliana) CB5D localizes to the endoplasmic reticulum membrane and physically associates with monolignol P450 enzymes. Disrupting CB5D in Arabidopsis resulted in a >60% reduction in S-lignin subunit levels but no impairment in G-lignin formation compared with the wild type, which sharply contrasts with the impaired G- and S-lignin synthesis observed after disrupting ATR2, encoding Arabidopsis CPR. The defective S-lignin synthesis in cb5d mutants was rescued by the expression of the gene encoding CB5D but not with mutant CB5D devoid of its electron shuttle properties. Disrupting ATR2 suppressed the catalytic activity of both cinnamic acid 4-hydroxylase and ferulate 5-hydroxylase (F5H), but eliminating CB5D specifically depleted the latter's activity. Therefore, CB5D functions as an obligate electron shuttle intermediate that specifically augments F5H-catalyzed reactions, thereby controlling S-lignin biosynthesis.

Lipid composition and macromolecular crowding effects on CYP2J2-mediated drug metabolism in nanodiscs

Lipid composition and macromolecular crowding are key external effectors of protein activity and stability whose role varies between different proteins. Therefore, it is imperative to study their effects on individual protein function. CYP2J2 is a membrane-bound cytochrome P450 in the heart involved in the metabolism of fatty acids and xenobiotics. In order to facilitate this metabolism, cytochrome P450 reductase (CPR), transfers electrons to CYP2J2 from NADPH. Herein, we use nanodiscs to show that lipid composition of the membrane bilayer affects substrate metabolism of the CYP2J2-CPR nanodisc (ND) system. Differential effects on both NADPH oxidation and substrate metabolism by CYP2J2-CPR are dependent on the lipid composition. For instance, sphingomyelin containing nanodiscs produced more secondary substrate metabolites than discs of other lipid compositions, implying a possible conformational change leading to processive metabolism. Furthermore, we demonstrate that macromolecular crowding plays a role in the lipid-solubilized CYP2J2-CPR system by increasing the K_m and decreasing the V_max , and effect that is size-dependent. Crowding also affects the CYP2J2-CPR-ND system by decreasing both the K_m and V_max for Dextran-based macromolecular crowding agents, implying an increase in substrate affinity but a lack of metabolism. Finally, protein denaturation studies show that crowding agents destabilize CYP2J2, while the multidomain protein CPR is stabilized. Overall, these studies are the first report on the role of the surrounding lipid environment and macromolecular crowding in modulating enzymatic function of CYP2J2-CPR membrane protein system.

SPR-Based study of affinity of cytochrome P450s / redox partners interactions modulated by steroidal substrates

The goal of this work was to test the hypothesis that the affinity of protein-protein interactions in the cytochrome P450-dependent monooxygenase system is modulated by the low-molecular-weight compounds (substrates or inhibitors). The surface plasmon resonance (SPR) based study was carried out using the recombinant protein preparations of three microsomal cytochromes P450 (CYP17A1, CYP21A2, and CYP2C19) and their redox partners: cytochrome b5 (CYB5A), NADPH - cytochrome P450 reductase (CPR), and also iron-sulfur protein adrenodoxin (Adx). As a result, we have revealed some specificity of the influence of the steroid substrates on the binding affinity of CYPs with their redox partners, namely: the lack of effect on CPR/CYPs and Adx/CYP complex formation, and a significant effect on interactions between CYB5A and steroidogenic CYPs. The equilibrium dissociation constant (Kd) value of the CYB5A/CYP17A1 complex decreased by 5 times in the presence of progesterone (P4), which was due to a 10 times increase in the association rate constant (k_on). In this case, a twofold increase in the dissociation rate constant (k_off) value of CYB5A/CYP17A1 complex formation was observed. It was also demonstrated that the affinity of CYB5A/CYP17A1 interaction increased in the presence of two other steroidal substrates 17α-hydroxyprogesterone and pregnenolone and that effect was comparable with P4. In contrast, only the twofold decrease in the affinity of CYB5A/CYP21A2 interaction in the presence of P4 was caused by a slight increase in the k_off value (the k_on value of the complex did not change). This indicates a different format of the steroidal substrates effects expressed in a change in the stability of the CYB5A/CYPs complexes. Thus, it was found that P4 modulated the both kinetic and equilibrium constants of CYB5A/CYP17A1 and CYB5/CYP21A2 complex formation and complexes, while not affecting the CYB5A/CYP2C19 interaction (2C19 is the cytochrome P450 isoenzyme possessing broad substrate specificity), thereby indicating a specific influence of steroidal substrates on interactions involving steroidogenic CYPs. Our results are consistent with current understanding of the role of CYB5A as a regulator of cytochrome P450 activity in P450-dependent monooxygenase system.

Probing the Role of the Hinge Segment of Cytochrome P450 Oxidoreductase in the Interaction with Cytochrome P450

NADPH-cytochrome P450 reductase (CPR) is the unique redox partner of microsomal cytochrome P450s (CYPs). CPR exists in a conformational equilibrium between open and closed conformations throughout its electron transfer (ET) function. Previously, we have shown that electrostatic and flexibility properties of the hinge segment of CPR are critical for ET. Three mutants of human CPR were studied (S243P, I245P and R246A) and combined with representative human drug-metabolizing CYPs (isoforms 1A2, 2A6 and 3A4). To probe the effect of these hinge mutations different experimental approaches were employed: CYP bioactivation capacity of pre-carcinogens, enzyme kinetic analysis, and effect of the ionic strength and cytochrome b₅ (CYB5) on CYP activity. The hinge mutations influenced the bioactivation of pre-carcinogens, which seemed CYP isoform and substrate dependent. The deviations of Michaelis-Menten kinetic parameters uncovered tend to confirm this discrepancy, which was confirmed by CYP and hinge mutant specific salt/activity profiles. CPR/CYB5 competition experiments indicated a less important role of affinity in CPR/CYP interaction. Overall, our data suggest that the highly flexible hinge of CPR is responsible for the existence of a conformational aggregate of different open CPR conformers enabling ET-interaction with structural varied redox partners.

Bioinformatics tools for the identification of gene clusters that biosynthesize specialized metabolites

Specialized metabolites (also called natural products or secondary metabolites) derived from bacteria, fungi, marine organisms and plants constitute an important source of antibiotics, anti-cancer agents, insecticides, immunosuppressants and herbicides. Many specialized metabolites in bacteria and fungi are biosynthesized via metabolic pathways whose enzymes are encoded by clustered genes on a chromosome. Metabolic gene clusters comprise a group of physically co-localized genes that together encode enzymes for the biosynthesis of a specific metabolite. Although metabolic gene clusters are generally not known to occur outside of microbes, several plant metabolic gene clusters have been discovered in recent years. The discovery of novel metabolic pathways is being enabled by the increasing availability of high-quality genome sequencing coupled with the development of powerful computational toolkits to identify metabolic gene clusters. To provide a comprehensive overview of various bioinformatics methods for detecting gene clusters, we compare and contrast key aspects of algorithmic logic behind several computational tools, including 'NP.searcher', 'ClustScan', 'CLUSEAN', 'antiSMASH', 'SMURF', 'MIDDAS-M', 'ClusterFinder', 'CASSIS/SMIPS' and 'C-Hunter' among others. We also review additional tools such as 'NRPSpredictor' and 'SBSPKS' that can infer substrate specificity for previously identified gene clusters. The continual development of bioinformatics methods to predict gene clusters will help shed light on how organisms assemble multi-step metabolic pathways for adaptation to various ecological niches.

Human hepatocytes and cytochrome P450-selective inhibitors predict variability in human drug exposure more accurately than human recombinant P450s

BACKGROUND AND PURPOSE

Drugs metabolically eliminated by several enzymes are less vulnerable to variable compound exposure in patients due to drug-drug interactions (DDI) or if a polymorphic enzyme is involved in their elimination. Therefore, it is vital in drug discovery to accurately and efficiently estimate and optimize the metabolic elimination profile.

EXPERIMENTAL APPROACH

CYP3A and/or CYP2D6 substrates with well described variability in vivo in humans due to CYP3A DDI and CYP2D6 polymorphism were selected for assessment of fraction metabolized by each enzyme (fmCYP ) in two in vitro systems: (i) human recombinant P450s (hrP450s) and (ii) human hepatocytes combined with selective P450 inhibitors. Increases in compound exposure in poor versus extensive CYP2D6 metabolizers and by the strong CYP3A inhibitor ketoconazole were mathematically modelled and predicted changes in exposure were compared with in vivo data.

KEY RESULTS

Predicted changes in exposure were within twofold of reported in vivo values using fmCYP estimated in human hepatocytes and there was a strong linear correlation between predicted and observed changes in exposure (r2  = 0.83 for CYP3A, r2  = 0.82 for CYP2D6). Predictions using fmCYP in hrP450s were not as accurate (r2  = 0.55 for CYP3A, r2  = 0.20 for CYP2D6).

CONCLUSIONS AND IMPLICATIONS

The results suggest that variability in human drug exposure due to DDI and enzyme polymorphism can be accurately predicted using fmCYP from human hepatocytes and CYP-selective inhibitors. This approach can be efficiently applied in drug discovery to aid optimization of candidate drugs with a favourable metabolic elimination profile and limited variability in patients.

Use of Phenoxyaniline Analogues To Generate Biochemical Insights into the Interactio n of Polybrominated Diphenyl Ether with CYP2B Enzymes

Human hepatic cytochromes P450 (CYP) are integral to xenobiotic metabolism. CYP2B6 is a major catalyst of biotransformation of environmental toxicants, including polybrominated diphenyl ethers (PBDEs). CYP2B substrates tend to contain halogen atoms, but the biochemical basis for this selectivity and for species specific determinants of metabolism has not been identified. Spectral binding titrations and inhibition studies were performed to investigate interactions of rat CYP2B1, rabbit CYP2B4, and CYP2B6 with a series of phenoxyaniline (POA) congeners that are analogues of PBDEs. For most congeners, there was a <3-fold difference between the spectral binding constants (K_S) and IC₅₀ values. In contrast, large discrepancies between these values were observed for POA and 3-chloro-4-phenoxyaniline. CYP2B1 was the enzyme most sensitive to POA congeners, so the Val-363 residue from that enzyme was introduced into CYP2B4 or CYP2B6. This substitution partially altered the protein-ligand interaction profiles to make them more similar to that of CYP2B1. Addition of cytochrome P450 oxidoreductase (POR) to titrations of CYP2B6 with POA or 2'4'5'TCPOA decreased the affinity of both ligands for the enzyme. Addition of cytochrome b₅ to a recombinant enzyme system containing POR and CYP2B6 increased the POA IC₅₀ value and decreased the 2'4'5'TCPOA IC₅₀ value. Overall, the inconsistency between K_S and IC₅₀ values for POA versus 2'4'5'TCPOA is largely due to the effects of redox partner binding. These results provide insight into the biochemical basis of binding of diphenyl ethers to human CYP2B6 and changes in CYP2B6-mediated metabolism that are dependent on POA congener and redox partner identity.

Evidence for a complex formation between CYP2J5 and mEH in living cells by FRET analysis of membrane protein interaction in the endoplasmic reticulum (FAMPIR)

The potential complex formation between microsomal epoxide hydrolase (mEH) and cytochrome P450-dependent monooxygenase (CYP) has been a subject of research for many decades. Such an association would enable efficient substrate channeling between CYP and mEH and as such represent an attractive strategy to prevent deleterious accumulation of harmful metabolic by-products such as CYP-generated epoxide intermediates. However, such complex formation is experimentally difficult to prove, because CYP and mEH are membrane-bound proteins that are prone to unspecific aggregation after solubilization. Here, we report the development of a FRET-based procedure to analyze the mEH–CYP interaction in living cells by fluorescence-activated cell sorting. With this non-invasive procedure, we demonstrate that CYP2J5 and mEH associate in the endoplasmic reticulum of recombinant HEK293 cells to the same extent as do CYP2J5 and its indispensible redox partner cytochrome P450 reductase. This presents final proof for a very close proximity of CYP and mEH in the endoplasmic reticulum, compatible with and indicative of their physical interaction. In addition, we provide with FAMPIR a robust and easy-to-implement general method for analyzing the interaction of ER membrane-resident proteins that share a type I topology.

Toward a systems approach to the human cytochrome P450 ensemble: interactions between CYP2D6 and CYP2E1 and their functional consequences

Functional cross-talk among human drug-metabolizing cytochrome P450 through their association is a topic of emerging importance. Here, we studied the interactions of human CYP2D6, a major metabolizer of psychoactive drugs, with one of the most prevalent human P450 enzymes, ethanol-inducible CYP2E1. Detection of P450-P450 interactions was accomplished through luminescence resonance energy transfer between labeled proteins incorporated into human liver microsomes and the microsomes of insect cells containing NADPH-cytochrome P450 reductase. The potential of CYP2D6 to form oligomers in the microsomal membrane is among the highest observed with human cytochrome P450 studied up to date. We also observed the formation of heteromeric complexes of CYP2D6 with CYP2E1 and CYP3A4, and found a significant modulation of these interactions by 3,4-methylenedioxymethylamphetamine, a widespread drug of abuse metabolized by CYP2D6. Our results demonstrate an ample alteration of the catalytic properties of CYP2D6 and CYP2E1 caused by their association. In particular, we demonstrated that preincubation of microsomes containing co-incorporated CYP2D6 and CYP2E1 with CYP2D6-specific substrates resulted in considerable time-dependent activation of CYP2D6, which presumably occurs via a slow substrate-induced reorganization of CYP2E1-CYP2D6 hetero-oligomers. Furthermore, we demonstrated that the formation of heteromeric complexes between CYP2E1 and CYP2D6 affects the stoichiometry of futile cycling and substrate oxidation by CYP2D6 by means of decreasing the electron leakage through the peroxide-generating pathways. Our results further emphasize the role of P450-P450 interactions in regulatory cross-talk in human drug-metabolizing ensemble and suggest a role of interactions of CYP2E1 with CYP2D6 in pharmacologically important instances of alcohol-drug interactions.

The functional effects of physical interactions involving cytochromes P450: putative mechanisms of action and the extent of these effects in biological membranes

Cytochromes P450 represent a family of enzymes, which are responsible for the oxidative metabolism of a wide variety of xenobiotics. Although the mammalian P450s require interactions with their redox partners in order to function, more recently, P450 system proteins have been shown to exist as multi-protein complexes that include the formation of P450•P450 complexes. Evidence has shown that the metabolism of some substrates by a given P450 can be influenced by the specific interaction of the enzyme with other forms of P450. Detailed kinetic analysis of these reactions in vitro has shown that the P450-P450 interactions can alter metabolism by changing the ability of a P450 to bind to its cognate redox partner, NADPH-cytochrome P450 reductase; by altering substrate binding to the affected P450; and/or by changing the rate of a catalytic step of the reaction cycle. This review summarizes the known examples of P450-P450 interactions that have been shown in vitro to influence metabolism and categorizes them according to the mechanism(s) causing the effects. P450-P450 interactions have the potential to cause major changes in the metabolism and elimination of drugs in vivo. This review summarizes the evidence that the P450-P450 interactions influence metabolism in biological membranes and discusses the studies, which will provide further insight into the extent of these effects in the future.

Advances in the Understanding of Protein-Protein Interactions in Drug Metabolizing Enzymes through the Use of Biophysical Techniques

In recent years, a growing appreciation has developed for the importance of protein-protein interactions to modulate the function of drug metabolizing enzymes. Accompanied with this appreciation, new methods and technologies have been designed for analyzing protein-protein interactions both in vitro and in vivo. These technologies have been applied to several classes of drug metabolizing enzymes, including: cytochrome P450's (CYPs), monoamine oxidases (MAOs), UDP-glucuronosyltransferases (UGTs), glutathione S-transferases (GSTs), and sulfotransferases (SULTs). In this review, we offer a brief description and assessment of the impact of many of these technologies to the study of protein-protein interactions in drug disposition. The still expanding list of these techniques and assays has the potential to revolutionize our understanding of how these enzymes carry out their important functions in vivo.

Organic anion transporting polypeptide 1B3 can form homo- and hetero-oligomers

OATP1B3 is a 12 transmembrane domain protein expressed at the basolateral membrane of human hepatocytes where it mediates the uptake of numerous drugs and endogenous compounds. Previous western blot results suggest the formation of OATP1B3 multimers. In order to better understand the function of OATP1B3 under normal physiological conditions, we investigated its oligomerization status. We transiently transfected OATP1B3 with a C-terminal His-, FLAG- or HA-tag in HEK293 cells and used co-immunoprecipitation and a Proximity Ligation Assay to detect interactions between the different constructs. All three constructs retained similar transport rates as wild-type OATP1B3. Immunofluorescence experiments indicated that in contrast to wild-type, His- and FLAG-tagged OATP1B3, where the C-terminal end is on the cytoplasmic side of the membrane, the C-terminal end of HA-tagged OATP1B3 is extracellular. After cross-linking, anti-FLAG antibodies were able to pull down FLAG-tagged OATP1B3 (positive control) and co-transfected His- or HA-tagged OATP1B3, demonstrating the formation of homo-oligomers and suggesting that the C-terminal part is not involved in oligomer formation. We confirmed co-localization of His- and FLAG-tagged OATP1B3 in transfected HEK293 cells with the Proximity Ligation Assay. Transport studies with a non-functional OATP1B3 mutant suggest that the individual subunits and not the whole oligomer are the functional units in the homo-oligomers. In addition, we also detected OATP1B3-FLAG co-localization with OATP1B1-His or NTCP-His, suggesting that OATP1B3 also hetero-oligomerizes with other transport proteins. Using the Proximity Ligation Assay with transporter specific antibodies, we demonstrate close association of OATP1B3 with NTCP in frozen human liver tissue. These findings demonstrate that OATP1B3 can form homo- and hetero-oligomers and suggest a potential co-regulation of the involved transporters.

Membrane Associated Progesterone Receptors: Promiscuous Proteins with Pleiotropic Functions - Focus on Interactions with Cytochromes P450

Membrane-associated progesterone receptors (MAPR) are a group of four rather small, partially homologous proteins, which share a similar non-covalent heme-binding domain that is related to cytochrome b5, a well-known functional interaction partner of microsomal cytochrome P450 (CYP) monooxygenase systems. Apart from their structural similarities the four proteins progesterone membrane component 1 (PGRMC1, also referred to as IZA, sigma-2 receptor, Dap1), PGRMC2, neudesin (NENF) and neuferricin (CYB5D2) display surprisingly divergent and multifunctional physiological properties related to cholesterol/steroid biosynthesis, drug metabolism and response, iron homeostasis, heme trafficking, energy metabolism, autophagy, apoptosis, cell cycle regulation, cell migration, neural functions, and tumorigenesis and cancer progression. The purpose of this mini-review is to briefly summarize the structural and functional properties of MAPRs with particular focus on their interactions with the CYP system. For PGRMC1, originally identified as a non-canonical progesterone-binding protein that mediates some immediate non-genomic actions of progesterone, available evidence indicates mainly activating interactions with steroidogenic CYPs including CYP11A1, CYP21A2, CYP17, CYP19, CYP51A1, and CYP61A1, while interactions with drug metabolizing CYPs including CYP2C2, CYP2C8, CYP2C9, CYP2E1, and CYP3A4 were either ineffective or slightly inhibitory. For the other MAPRs the evidence is so far less conclusive. We also point out that experimental limitations question some of the previous conclusions. Use of appropriate model systems should help to further clarify the true impact of these proteins on CYP-mediated metabolic pathways.