Cytochrome P450 structure–function: insights from molecular dynamics simulations

Small angle X-ray scattering analysis of thermophilic cytochrome P450 CYP119 and the effects of the N-terminal histidine tag

Combining size exclusion chromatography-small angle X-ray scattering (SEC-SAXS) and molecular dynamics (MD) analysis is a promising approach to investigate protein behavior in solution, particularly for understanding conformational changes due to substrate binding in cytochrome P450s (CYPs). This study investigates conformational changes in CYP119, a thermophilic CYP from Sulfolobus acidocaldarius that exhibits structural flexibility similar to mammalian CYPs. Although the crystal structure of ligand-free (open state) and ligand-bound (closed state) forms of CYP119 is known, the overall structure of the enzyme in solution has not been explored until now. It was found that theoretical scattering profiles from the crystal structures of CYP119 did not align with the SAXS data, but conformers from MD simulations, particularly starting from the open state (46 % of all frames), agreed well. Interestingly, a small percentage of closed-state conformers also fit the data (9 %), suggesting ligand-free CYP119 samples ligand-bound conformations. Ab initio SAXS models for N-His tagged CYP119 revealed a tail-like unfolded structure impacting protein flexibility, which was confirmed by in silico modeling. SEC-SAXS analysis of N-His CYP119 indicated pentameric structures in addition to monomers in solution, affecting the stability and activity of the enzyme. This study adds insights into the conformational dynamics of CYP119 in solution.

Assessing the role of residue Phe108 of cytochrome P450 3A4 in allosteric effects of midazolam metabolism

Cytochrome P450 3A4 (CYP3A4) is involved in the metabolism of more drugs in clinical use than any other xenobiotic-metabolizing enzyme. CYP3A4-mediated drug metabolism is usually allosterically modulated by substrate concentration (homotropic allostery) and other drugs (heterotropic allostery), exhibiting unusual kinetic profiles and regiospecific metabolism. Recent studies suggest that residue Phe108 (F108) of CYP3A4 may have an important role in drug metabolism. In this work, residue mutations were coupled with well-tempered metadynamics simulations to assess the importance of F108 in the allosteric effects of midazolam metabolism. Comparing the simulation results of the wild-type and mutation systems, we identify that the π-π interaction and steric effect between the F108 side chain and midazolam is favorable for the stable binding of substrate in the active site. F108 also plays an important role in the transition of substrate binding mode, which mainly induces the transition of substrate binding mode by forming π-π interactions with multiple aromatic rings of the substrate. Moreover, the side chain of F108 is closely related to the radius and depth of the 2a and 2f channels, and F108 may further regulate drug metabolism by affecting the pathway, orientation, or time of substrate entry into the CYP3A4 active site or product egress from the active site. Altogether, we suggest that F108 affects drug metabolism and regulatory mechanisms by affecting substrate binding stability, binding mode transition, and channel characteristics of CYP3A4. Our findings could promote the understanding of complicated allosteric mechanisms in CYP3A4-mediated drug metabolism, and the knowledge could be used for drug development and disease treatment.

Water Exchange from the Buried Binding Sites of Cytochrome P450 Enzymes 1A2, 2D6, and 3A4 Correlates with Conformational Fluctuations

Human cytochrome P450 enzymes (CYPs) are critical for the metabolism of small-molecule pharmaceuticals (drugs). As such, the prediction of drug metabolism by and drug inhibition of CYP activity is an important component of the drug discovery and design process. Relative to the availability of a wide range of experimental atomic-resolution CYP structures, the development of structure-based CYP activity models has been limited. To better characterize the role of CYP conformational fluctuations in CYP activity, we perform multiple microsecond-scale all-atom explicit-solvent molecular dynamics (MD) simulations on three CYP isoforms, 1A2, 2D6, and 3A4, which together account for the majority of CYP-mediated drug metabolism. The MD simulations employ a variety of positional restraints, ranging from keeping all CYP atoms close to their experimentally determined coordinates to allowing full flexibility. We find that, with full flexibility, large fluctuations in the CYP binding sites correlate with efficient water exchange from these buried binding sites. This is especially true for 1A2, which, when restrained to its crystallographic conformation, is unable to exchange water between the binding site and bulk solvent. These findings imply that, in addition to crystal structures, a representative ensemble of conformational states ought to be included when developing structure-based CYP activity models.

Computational and artificial intelligence-based approaches for drug metabolism and transport prediction

Drug metabolism and transport, orchestrated by drug-metabolizing enzymes (DMEs) and drug transporters (DTs), are implicated in drug-drug interactions (DDIs) and adverse drug reactions (ADRs). Reliable and precise predictions of DDIs and ADRs are critical in the early stages of drug development to reduce the rate of drug candidate failure. A variety of experimental and computational technologies have been developed to predict DDIs and ADRs. Recent artificial intelligence (AI) approaches offer new opportunities for better predicting and understanding the complex processes related to drug metabolism and transport. We summarize the role of major DMEs and DTs, and provide an overview of current progress in computational approaches for the prediction of drug metabolism, transport, and DDIs, with an emphasis on AI including machine learning (ML) and deep learning (DL) modeling.

The molecular basis of dapsone activation of CYP2C9-catalyzed nonsteroidal anti-inflammatory drug oxidation

Positive heterotropic cooperativity, or "activation," results in an instantaneous increase in enzyme activity in the absence of an increase in protein expression. Thus, cytochrome P450 (CYP) enzyme activation presents as a potential drug-drug interaction mechanism. It has been demonstrated previously that dapsone activates the CYP2C9-catalyzed oxidation of a number of nonsteroidal anti-inflammatory drugs in vitro. Here, we conducted molecular dynamics simulations (MDS) together with enzyme kinetic investigations and site-directed mutagenesis to elucidate the molecular basis of the activation of CYP2C9-catalyzed S-flurbiprofen 4'-hydroxylation and S-naproxen O-demethylation by dapsone. Supplementation of incubations of recombinant CYP2C9 with dapsone increased the catalytic efficiency of flurbiprofen and naproxen oxidation by 2.3- and 16.5-fold, respectively. MDS demonstrated that activation arises predominantly from aromatic interactions between the substrate, dapsone, and the phenyl rings of Phe114 and Phe476 within a common binding domain of the CYP2C9 active site, rather than involvement of a distinct effector site. Mutagenesis of Phe114 and Phe476 abrogated flurbiprofen and naproxen oxidation, and MDS and kinetic studies with the CYP2C9 mutants further identified a pivotal role of Phe476 in dapsone activation. MDS additionally showed that aromatic stacking interactions between two molecules of naproxen are necessary for binding in a catalytically favorable orientation. In contrast to flurbiprofen and naproxen, dapsone did not activate the 4'-hydroxylation of diclofenac, suggesting that the CYP2C9 active site favors cooperative binding of nonsteroidal anti-inflammatory drugs with a planar or near-planar geometry. More generally, the work confirms the utility of MDS for investigating ligand binding in CYP enzymes.

Pb exposure causes non-linear accumulation of Pb in D. melanogaster controlled by metallothionein B and exerts ecological effects

Pb pollution can harm human health and the ecosystem. Therefore, it is worthwhile to study the metabolic processes of heavy metals in individual bodies and their influence on ecological systems. In this work, we analyzed the genetic responses and physiological changes of D. melanogaster which took diets exposed to different doses of Pb using transcriptomic analysis, ICP-MS, and various other physiological methods. We found that the Pb accumulated in D. melanogaster in a nonlinear pattern with the increase of Pb content in food. Metallothioneins (Mtns), especially the MtnB directly affects the accumulation and excretion of metal Pb in D. melanogaster, and causes the nonlinear accumulation. Metal regulatory transcription factor-1 (MTF-1) is involved in the regulation of Pb-induced high expressions of Mtns. Furthermore, an interaction between the metal metabolism pathway and xenobiotic response pathway leads to the cross-tolerances of Pb-exposed D. melanogaster to insecticides and other toxins. The oxidative stress induced by Pb toxicity may be the bridge between them. Our findings provide a physiological and molecular genetic basis for further study of the accumulation and metabolism of Pb in D. melanogaster.

Differences of Atomic-Level Interactions between Midazolam and Two CYP Isoforms 3A4 and 3A5

CYP 3A4 and CYP 3A5 are two important members of the human cytochrome P450 family. Although their overall structures are similar, the local structures of the active site are different, which directly leads to obvious individual differences in drug metabolic efficacy and toxicity. In this work, midazolam (MDZ) was selected as the probe substrate, and its interaction with two proteins, CYP 3A4 and CYP 3A5, was studied by molecular dynamics simulation (MD) along with the calculation of the binding free energy. The results show that two protein-substrate complexes have some similarities in enzyme-substrate binding; that is, in both complexes, Ser119 forms a high occupancy hydrogen bond with MDZ, which plays a key role in the stability of the interaction between MDZ and the enzymes. However, the complex formed by CYP 3A4 and MDZ is more stable, which may be attributed to the sandwich structure formed by the fluorophenyl group of the substrate with Leu216 and Leu482. Our study interprets the binding differences between two isoform-substrate complexes and reveals a structure-function relationship from the atomic perspective, which is expected to provide a theoretical basis for accurately measuring the effectiveness and toxicity of drugs for individuals in the era of precision medicine.

Hydration and Structural Adaptations of the Human CYP1A1, CYP1A2, and CYP1B1 Active Sites by Molecular Dynamics Simulations

Cytochromes CYP1A1, CYP1A2, and CYP1B1, the members of the cytochrome P450 family 1, catalyze the metabolism of endogenous compounds, drugs, and non-drug xenobiotics which include substances involved in the process of carcinogenesis, cancer chemoprevention, and therapy. In the present study, the interactions of three selected polymethoxy-trans-stilbenes, analogs of a bioactive polyphenol trans-resveratrol (3,5,4'-trihydroxy-trans-stilbene) with the binding sites of CYP1 isozymes were investigated with molecular dynamics (MD) simulations. The most pronounced structural changes in the CYP1 binding sites were observed in two substrate recognition sites (SRS): SRS2 (helix F) and SRS3 (helix G). MD simulations show that the number and position of water molecules occurring in CYP1 APO and in the structures complexed with ligands are diverse. The presence of water in binding sites results in the formation of water-protein, water-ligand, and bridging ligand-water-protein hydrogen bonds. Analysis of the solvent and substrate channels opening during the MD simulation showed significant differences between cytochromes in relation to the solvent channel and the substrate channels 2c, 2ac, and 2f. The results of this investigation lead to a deeper understanding of the molecular processes that occur in the CYP1 binding sites and may be useful for further molecular studies of CYP1 functions.

Therapeutic treatments for diabetes mellitus-induced liver injury by regulating oxidative stress and inflammation

Diabetes mellitus (DM) is a metabolic disease that affects all systems in the body, including the liver. Numerous studies have reported that chronic DM etiology and pathogenesis complications implicate oxidative stress, generating reactive oxygen species, such as superoxide anions and free radicals. In addition, pro-inflammatory reactions are also underlying functions closely related to oxidative stress that further exacerbate pathological DM states. The liver is especially susceptible to hyperglycemia-induced oxidative stress and the related inflammation. Thus, anti-oxidation and anti-inflammation therapies are promising strategies for treating liver damage. This review summarizes therapeutic treatments attenuating the generation of oxidative stress and pro-inflammation, which also cause DM-induced liver injury. Although the treatments have several impediments to be solved, these remedies may have clinically important implications under the absence of effective drugs for the damaged liver in DM patients.

Prediction of Sites of Metabolism of CYP3A4 Substrates Utilizing Docking-Derived Geometric Features

Cytochrome P450 3A4 (CYP3A4) is one of the major drug-metabolizing enzymes in the human body and is responsible for the metabolism of ∼50% of clinically used drugs. Therefore, the identification of the compound's sites of metabolism (SOMs) mediated by CYP3A4 is of utmost importance in the early stage of drug discovery and development. Herein, docking-based approaches incorporating geometric features were used for SOMs prediction of CYP3A4 substrates. The cross-docking poses of a relatively large data set containing 474 substrates were analyzed in depth, and a widely observed geometric pattern called the close proximity of SOMs was derived from the poses. On the basis of the close proximity, several structure-based models have been constructed, which demonstrated better performance than those structure-based models using the criterion of Fe-SOM distance. For further improving the prediction performance, the structure-based models were also combined with the well-known ligand-based model SMARTCyp. One combined model exhibited good performance on the SOMs prediction of an external substrate set containing kinase inhibitors, PROTACs, approved drugs, and some lead compounds.

Bioefficacy of Nga-Mon (Perilla frutescens) Fresh and Dry Leaf: Assessment of Antioxidant, Antimutagenicity, and Anti-Inflammatory Potential

Perilla leaves are known to be a rich source of polyphenols, which have been shown to exhibit various biological effects. This study aimed to compare the bioefficacies and bioactivities of fresh (PLE_f) and dry (PLE_d) Thai perilla (Nga-mon) leaf extracts. Phytochemical analysis indicated that both PLE_f and PLE_d were abundant in rosmarinic acid and bioactive phenolic compounds. PLE_d, which had higher levels of rosmarinic acid but lower levels of ferulic acid and luteolin than PLE_f, exhibited greater effectiveness in a free radical scavenging assay. Furthermore, both extracts were found to suppress intracellular ROS generation and exhibit antimutagenic activity against food-borne carcinogens in S. typhimurium. They also attenuated lipopolysaccharide-induced inflammation in RAW 264.7 cells by inhibiting the expression of nitric oxide, iNOS, COX-2, TNF-α, IL-1β, and IL-6 through the suppression of NF-κB activation and translocation. However, PLE_f exhibited a higher ability to suppress cellular ROS production and higher antimutagenic and anti-inflammatory activities than PLE_d, which can be attributed to its combination of phytochemical components. Overall, PLE_f and PLE_d have the potential to serve as natural bioactive antioxidant, antimutagenic, and anti-inflammatory agents to achieve potential health benefits.

Cytochrome P450 superfamily in spotted sea bass: Genome-wide identification and expression profiles under trichlorfon and environmental stresses

Cytochrome P450s (CYPs), as one of the most diverse enzyme superfamilies in nature, play critical functions in antioxidant reactions against endogenous and exogenous compounds. In this study, we performed genome-wide characterization of CYP superfamily members and analyzed their expression patterns under several abiotic stresses in spotted sea bass, which is known as an economically important fish species in the Chinese aquaculture industry. A total of 55 CYP genes were identified and divided into 17 families within 10 clans. The analysis of phylogeny, gene structure, and syntenic relationships provided evidence for the evolution of CYP genes and confirmed their annotation and orthology. The expression of CYP genes was examined in the liver during trichlorfon stress using quantitative real-time PCR. The results showed that 20 tested CYP genes displayed significant mRNA expression changes, indicating that they may play crucial roles in the metabolism of trichlorfon and can be potential biomarkers for trichlorfon pollution. Moreover, by screening transcriptomic databases, 10, 3 and 19 CYP genes exhibited differential expression patterns in response to hypoxia, alkalinity and heat stress, respectively. Taken together, this study provided insights into the regulation of CYP genes by toxicological and environmental stresses, laid basis for extensive functional studies of the CYP superfamily in spotted sea bass and other teleost species.

Comprehensive in vitro and in silico assessments of metabolic capabilities of 24 genomic variants of CYP2C19 using two different substrates

Introduction: Most hepatically cleared drugs are metabolized by cytochromes P450 (CYPs), and Clinical Pharmacogenetics Implementation Consortium (CPIC) guidelines provide curated clinical references for CYPs to apply individual genome data for optimized drug therapy. However, incorporating novel pharmacogenetic variants into guidelines takes considerable time. Methods: We comprehensively assessed the drug metabolizing capabilities of CYP2C19 variants discovered through population sequencing of two substrates, S-mephenytoin and omeprazole. Results: Based on established functional assays, 75% (18/24) of the variants not yet described in Pharmacogene Variation (PharmVar) had significantly altered drug metabolizing capabilities. Of them, seven variants with inappreciable protein expression were evaluated as protein damaging by all three in silico prediction algorithms, Sorting intolerant from tolerant (SIFT), Polymorphism Phenotyping v2 (PolyPhen-2), and Combined annotation dependent depletion (CADD). The five variants with decreased metabolic capability (<50%) of wild type for either substrates were evaluated as protein damaging by all three in silico prediction algorithms, except CADD exact score of NM_000769.4:c.593T>C that was 19.68 (<20.0). In the crystal structure of the five polymorphic proteins, each altered residue of all those proteins was observed to affect the key structures of drug binding specificity. We also identified polymorphic proteins indicating different tendencies of metabolic capability between the two substrates (5/24). Discussion: Therefore, we propose a methodology that combines in silico prediction algorithms and functional assays on polymorphic CYPs with multiple substrates to evaluate the changes in the metabolism of all possible genomic variants in CYP genes. The approach would reinforce existing guidelines and provide information for prescribing appropriate medicines for individual patients.

Effects of transgenic cotton lines expressing dsAgCYP6CY3-P1 on the growth and detoxification ability of Aphis gossypii glover

BACKGROUND

The pest Aphis gossypii Glover globally causes considerable economic losses on various crops by its feeding damage and disease transmission. Transgenic plants that produce double-stranded RNA (dsRNA) targeted to insect genes are being developed as a pest control strategy. In this study, we evaluated the effects of transgenic cotton-mediated RNA interference (RNAi) on the growth and detoxification ability of A. gossypii after the transgenic cotton lines expressing dsAgCYP6CY3-P1 (the TG cotton lines) were obtained on the basis of exploring the functions of CYP6CY3 in our previous research.

RESULTS

The developmental time of third- and fourth-instar nymphs which fed on the TG cotton lines were significantly prolonged. Life table parameters showed that the fitness of cotton aphids from the TG cotton lines decreased. Additionally, the relative expression level of CYP6CY3 in cotton aphids which fed on the TG cotton lines was significantly reduced by 47.3 % at 48 h compared with that from the nontransgenic cotton (the NT cotton). Bioassay showed that silencing of CYP6CY3 increased mortality of the nymphs to imidacloprid by 28.49 % (at 24 h) and to acetamiprid by 73.77 % (at 48 h), respectively.

CONCLUSION

These results indicated that the TG cotton lines delayed the growth and development of A. gossypii, but also decreased population density and increased its sensitivity to imidacloprid and acetamiprid, respectively. The results provide further support for the development and application of plant-mediated RNAi. © 2022 Society of Chemical Industry.

Identification and in vitro functional assessment of 10 CYP2C9 variants found in Chinese Han subjects

Cytochrome P450 2C9 (CYP2C9) participates in about 15% of clinical drug metabolism, and its polymorphism is associated with individual drug metabolism differences, which may lead to the adverse drug reactions (ADRs). In this study, 1163 Chinese Han individuals were recruited to investigate their distribution pattern of CYP2C9 gene and find out the variants that may affect their drug metabolic activities. We successfully developed a multiplex PCR amplicon sequencing method and used it for the genetic screening of CYP2C9 in a large scale. Besides the wild type CYP2C9*1, totally 26 allelic variants of CYP2C9 were detected, which included 16 previously reported alleles and 10 new non-synonymous variants that had not been listed on the PharmVar website. The characteristics of these newly detected CYP2C9 variants were then evaluated after co-expressing them with CYPOR in S. cerevisiae microsomes. Immunoblot analysis revealed that except for Pro163Ser, Glu326Lys, Gly431Arg and Ile488Phe, most of newly detected variants showed comparable protein expression levels to wild type in yeast cells. Two typical CYP2C9 probe drugs, losartan and glimepiride, were then used for the evaluation of metabolic activities of variants. As a result, 3 variants Thr301Met, Glu326Lys, and Gly431Arg almost lost their catalytic activities and most of other variants exhibited significantly elevated activities for drug metabolism. Our data not only enriches the knowledge of naturally occurring CYP2C9 variants in the Chinese Han population, but also provides the fundamental evidence for its potential clinical usage for personalized medicine in the clinic.

Predicting Regioselectivity of AO, CYP, FMO, and UGT Metabolism Using Quantum Mechanical Simulations and Machine Learning

Unexpected metabolism in modification and conjugation phases can lead to the failure of many late-stage drug candidates or even withdrawal of approved drugs. Thus, it is critical to predict the sites of metabolism (SoM) for enzymes, which interact with drug-like molecules, in the early stages of the research. This study presents methods for predicting the isoform-specific metabolism for human AOs, FMOs, and UGTs and general CYP metabolism for preclinical species. The models use semi-empirical quantum mechanical simulations, validated using experimentally obtained data and DFT calculations, to estimate the reactivity of each SoM in the context of the whole molecule. Ligand-based models, trained and tested using high-quality regioselectivity data, combine the reactivity of the potential SoM with the orientation and steric effects of the binding pockets of the different enzyme isoforms. The resulting models achieve κ values of up to 0.94 and AUC of up to 0.92.

A gene cluster in Ginkgo biloba encodes unique multifunctional cytochrome P450s that initiate ginkgolide biosynthesis

The ginkgo tree (Ginkgo biloba) is considered a living fossil due to its 200 million year's history under morphological stasis. Its resilience is partly attributed to its unique set of specialized metabolites, in particular, ginkgolides and bilobalide, which are chemically complex terpene trilactones. Here, we use a gene cluster-guided mining approach in combination with co-expression analysis to reveal the primary steps in ginkgolide biosynthesis. We show that five multifunctional cytochrome P450s with atypical catalytic activities generate the tert-butyl group and one of the lactone rings, characteristic of all G. biloba trilactone terpenoids. The reactions include scarless C-C bond cleavage as well as carbon skeleton rearrangement (NIH shift) occurring on a previously unsuspected intermediate. The cytochrome P450s belong to CYP families that diversifies in pre-seed plants and gymnosperms, but are not preserved in angiosperms. Our work uncovers the early ginkgolide pathway and offers a glance into the biosynthesis of terpenoids of the Mesozoic Era.

Computer-Aided (In Silico) Modeling of Cytochrome P450-Mediated Food–Drug Interactions (FDI)

Modifications of the activity of Cytochrome 450 (CYP) enzymes by compounds in food might impair medical treatments. These CYP-mediated food–drug interactions (FDI) play a major role in drug clearance in the intestine and liver. Inter-individual variation in both CYP expression and structure is an important determinant of FDI. Traditional targeted approaches have highlighted a limited number of dietary inhibitors and single-nucleotide variations (SNVs), each determining personal CYP activity and inhibition. These approaches are costly in time, money and labor. Here, we review computational tools and databases that are already available and are relevant to predicting CYP-mediated FDIs. Computer-aided approaches such as protein–ligand interaction modeling and the virtual screening of big data narrow down hundreds of thousands of items in databanks to a few putative targets, to which the research resources could be further directed. Structure-based methods are used to explore the structural nature of the interaction between compounds and CYP enzymes. However, while collections of chemical, biochemical and genetic data are available today and call for the implementation of big-data approaches, ligand-based machine-learning approaches for virtual screening are still scarcely used for FDI studies. This review of CYP-mediated FDIs promises to attract scientists and the general public.

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

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

Molecular Basis of the Recognition of Cholesterol by Cytochrome P450 46A1 along the Major Access Tunnel

CYP46A1 is an important potential target for the treatment of Alzheimer's disease (AD), which is the most common neurodegenerative disease among older individuals. However, the binding mechanism between CYP46A1 and substrate cholesterol (CH) has not been clarified and will not be conducive to the research of relevant drug molecules. In this study, we integrated molecular docking, molecular dynamics (MD) simulations, and adaptive steered MD simulations to explore the recognition and binding mechanism of CYP46A1 with CH. Two key factors affecting the interaction between CH and CYP46A1 are determined: one is a hydrophobic cavity formed by seven hydrophobic residues (F80, Y109, L112, I222, W368, F371, and T475), which provides nonpolar interactions to stabilize CH, and the other is a hydrogen bond formed by H81 and CH, which ensures the binding direction of CH. In addition, the tunnel analysis results show that tunnel 2a is identified as the primary pathway of CH. The entry of CH induces tunnel 2e to close and tunnel w to open. Our results may provide effective clues for the design of drugs based on the substrate for AD and improve our understanding of the structure-function of CYP46A1.

Identification of Novel Molecular Targets of Four Microcystin Variants by High-Throughput Virtual Screening

Highly toxic microcystins (MCs) perform complex interactions with many proteins that induce cellular dysregulation, leading to the development of several diseases including cancer. There is significant diversity and chemical complexity among MC congeners, which makes it difficult to identify structure-dependent toxicity outcomes and their long-term effects. The aim of this study was to exploratory identify likely molecular targets of the main MC variants (MC-LA, MC-LR, MC-RR, and MC-LY) by conducting a computational binding affinity analysis using AutoDock Vina to evaluate the interaction of the toxins with 1000 proteins related to different biological functions. All four variants showed strong in silico interactions with proteins that regulate metabolism/immune system, CD38 (top scoring hit, −11.5 kcal/mol); inflammation, TLR4 (−11.4 kcal/mol) and TLR8 (−11.5 kcal/mol); neuronal conduction, BChE; renin–angiotensin signaling, (ACE); thyroid hormone homeostasis (TTR); and cancer-promoting processes, among other biochemical activities. The results show MCs have the potential to bind onto distinct molecular targets which could generate biochemical alterations through a number of signal transduction pathways. In short, this study broadens our knowledge about the mechanisms of action of different variants of microcystins and provides information for future direct experimentation.

In vitro evaluation of the inhibition potential of echinacoside on human cytochrome P450 isozymes

Background

Echinacoside (ECH) possesses a wide range of biological activity. This present study analyzes the effect of ECH on cytochrome P450 isozymes (CYPs) activities of human liver microsomes.

Methods

The effect of ECH on CYPs enzyme activities were studied using the enzyme-selective substrates phenacetin (1A2), chlorzoxazone (2E1), S-mephenytoin (2C19), testosterone (3A4), coumarin (2A6), diclofenac (2C9), paclitaxel (2C8), and dextromethorphan (2D6). The IC50 values for CYP1A2, CYP2E1, CYP2C19, and CYP3A4 isoforms were examined to express the strength of inhibition. Further, the inhibition of CYPs was checked for time-dependent or not, and then fitted with competitive or non-competitive inhibition models. The corresponding parameters were also obtained.

Results

ECH caused inhibitions on CYP1A2, CYP2E1, CYP2C19 and CYP3A4 enzyme activities in HLMs with IC50 of 21.23, 19.15, 8.70 and 55.42 μM, respectively. The obtained results showed that the inhibition of ECH on CYP3A4 was time-dependent with the KI/K_inact value of 6.63/0.066 min^− 1·μM^− 1. Moreover, ECH inhibited the activity of CYP1A2 and CYP2E1 via non-competitive manners (K_i = 10.90 μM and K_i = 14.40 μM, respectively), while ECH attenuated the CYP2C19 activity via a competitive manner (K_i = 4.41 μM).

Conclusions

The results of this study indicate that ECH inhibits CYP1A2, CYP2E1, CYP2C19 and CYP3A4 activities in vitro. In vivo and clinical studies are warranted to verify the relevance of these interactions.

Perspective of structural flexibility on selective inhibition towards CYP1B1 over CYP1A1 by α-naphthoflavone analogs

Research on action selectivity between CYP1A1 and CYP1B1 is particularly valuable for cancer chemoprevention and chemotherapy. However, they share a very close similarity in their ligand-binding pockets that α-naphthoflavone (ANF) is the co-crystal ligand for both isoforms, which poses a major challenge in revealing their selectivity mechanism. Therefore, three selective CYP1B1 inhibitors derived from ANF were selected to illustrate the structural basis for the selectivity between the two isoforms via a comprehensive computational strategy. It was found that the sustainability of the π-π stacking interactions with the phenylalanine residues of the two isoforms, namely, Phe123, Phe224, and Phe258 for CYP1A1, and Phe134, Phe231, and Phe268 for CYP1B1, played a crucial role in determining the selectivity of ligands with a classic aromatic conjugation system like ANF and its derivatives for CYP1B1 versus CYP1A1. Of note, the structural flexibility of the corresponding protein domains mainly orchestrated the sustainability of the corresponding π-π stacking interactions, thereby determining the binding selectivity. Therefore, the structure modification of naphthoflavone lead compounds into preferable binding configurations to satisfy the π-π stacking interactions of the key phenylalanine residues within CYP1B1 would be an inspiring strategy devised to improve the inhibitory selectivity towards CYP1B1. Collectively, this study revealed valuable insight into understanding the selective mechanism between CYP1A1 and CYP1B1 from the perspective of structural flexibility, which sheds light on the future rational design of CYP1B1 selective inhibitors.

Massively parallel characterization of CYP2C9 variant enzyme activity and abundance

CYP2C9 encodes a cytochrome P450 enzyme responsible for metabolizing up to 15% of small molecule drugs, and CYP2C9 variants can alter the safety and efficacy of these therapeutics. In particular, the anti-coagulant warfarin is prescribed to over 15 million people annually and polymorphisms in CYP2C9 can affect individual drug response and lead to an increased risk of hemorrhage. We developed click-seq, a pooled yeast-based activity assay, to test thousands of variants. Using click-seq, we measured the activity of 6,142 missense variants in yeast. We also measured the steady-state cellular abundance of 6,370 missense variants in a human cell line by using variant abundance by massively parallel sequencing (VAMP-seq). These data revealed that almost two-thirds of CYP2C9 variants showed decreased activity and that protein abundance accounted for half of the variation in CYP2C9 function. We also measured activity scores for 319 previously unannotated human variants, many of which may have clinical relevance.

Tryptophan-75 Is a Low-Energy Channel-Gating Residue that Facilitates Substrate Egress/Access in Cytochrome P450 2D6

CYP2D6 is a major drug metabolizing enzyme with a buried active site. Channels leading to the active site from various enzyme surfaces are believed to facilitate ligand egress and access to the active site. The present study used molecular dynamics (MD) and in vitro studies with CYP2D6*1 and a Trp75-to-Ala mutant to examine channel gating in CYP2D6 by Trp75. MD simulations measured energy landscapes of Trp75 conformations and simulated substrate passage within channel 2b using bufuralol as a model substrate. Trp75 alternated between multiple stable states that supported substrate transport along channel 2b with low-energy barriers between states (∼ -1 kcal/mol). Trp75 conformations were stabilized primarily by hydrogen bonding between Trp75 and Glu222, Asn226, Ala225, or Gln72. Energy barriers were low between Trp75 conformations, allowing Trp75 to easily move between various conformations over time and to function in both binding to and moving substrates in the 2b channel of CYP2D6. Michaelis-Menten kinetic studies completed with purified enzyme in a reconstituted system showed overall reduced enzyme efficiency for metabolism of bufuralol and dextromethorphan by the Trp75Ala mutant compared with CYP2D6*1. In stopped-flow measurements, k off for dextromethorphan was decreased in the absence of Trp75. Our results support a role for Trp75 in substrate shuttling to the active site of CYP2D6. SIGNIFICANCE STATEMENT: Using combined molecular dynamics and in vitro assays, this study shows for the first time a role for Trp75 as a channel entrance gating residue in the mechanism of substrate binding/unbinding in CYP2D6. Energy landscapes derived from molecular dynamics were used to quantitate the strength of gating, and kinetics assays showed the impact on enzyme efficiency and k off of a Trp75Ala mutation.

Catalytic approach to in vivo metabolism of atractylenolide III using biomimetic iron–porphyrin complexes

In vivo oxygenation of atractylenolide III during the metabolism and in situ oxidation mechanism by an iron–porphyrin iron(iv)-oxo porphyrin π-cation–radical complex.

Deciphering Reaction Determinants of Altered-Activity CYP2D6 Variants by Well-Tempered Metadynamics Simulation and QM/MM Calculations

The xenobiotic metabolizing enzyme CYP2D6 is the P450 cytochrome family member with the highest rate of polymorphism. This causes changes in the enzyme activity and specificity, which can ultimately lead to adverse reactions during drug treatment. To avoid or lower CYP-related toxicity risks, prediction of the most likely positions within a molecule where a metabolic reaction might occur is paramount. In order to obtain accurate predictions, it is crucial to understand all phenomena within the active site of the enzyme that contribute to an efficient substrate recognition and the subsequent catalytic reaction together with their relative weight within the overall thermodynamic context. This study aims to define the weight of the driving forces upon the C-H bond activation within CYP2D6 wild-type and a clinically relevant allelic variant with increased activity (CYP2D6*53) featuring two amino acid mutations in close vicinity of the heme. First, we investigated the steric and electrostatic complementarity of the substrate bufuralol using well-tempered metadynamics simulations with the aim to obtain the free energy profiles for each site of metabolism (SoM) within the different active sites. Second, the stereoelectronic complementarity was determined for each SoM within the two different active-site environments. Relying on the well-tempered metadynamics simulation energy profiles of each SoM, we identified the binding mode that was closest to the preferred transition-state geometry for efficient C-H bond activation. The binding modes were then used as starting structures for the quantum mechanics/molecular mechanics calculations performed to quantify the corresponding activation barriers. Our results show the relevance of the steric component in orienting the SoM in an energetically accessible position toward the heme. However, the corresponding intrinsic reactivity and electronic complementarity within the active site must be accurately evaluated in order to obtain a meaningful reaction prediction, from which the predominant SoM can be determined. The F120I mutation lowered the activation barrier for the major site and one of the minor SoMs. However, it had an impact neither on the CYP2D6 enantioselectivity preference of the oxidation reaction nor on the stereoselectivity from the substrate point of view.

Toxicity of Dithiothreitol (DTT) to Drosophila melanogaster

The thiol-containing compound Dithiothreitol (DTT) has been shown to be toxic to cultured cells by inducing the generation of reactive oxygen species that ultimately cause cell death. However, its effects on multicellular organisms and the environment have not been investigated yet in detail. In this work, we tested the toxicity of DTT to the model insect Drosophila melanogaster. We show that DTT is lethal to D. melanogaster by topical application but not through feeding. DTT treatment triggers the transcription of the canonical apoptosis regulators grim, hid and rpr at low amounts. The amplitude of this induction declines with elevating DTT amounts. By live microscopy, we observe apoptotic cells especially in the gut of DTT treated flies. In parallel, low DTT amounts also activate the expression of the cuticle barrier component gene snsl. This indicates that a physical defence response is launched upon DTT contact. This combined measure is seemingly successful in preventing fly death. The expression of a number of known detoxification genes including cyp6a2, cyp6a8, cyp12d1 and GstD2 is also enhanced through DTT contact. The degree of upregulation of these genes is proportional to the applied DTT amounts. Despite this effort, flies exposed to high amounts of DTT eventually die. Together, D. melanogaster is able to sense DTT toxicity and adjust the defence response successfully at least at low concentrations. This is the first time to analyse the molecular consequences of DTT exposure in a multicellular organism. Our work provides a new model to discuss the physiological response of animals against thiol toxins and to resurvey the effect of redox agents on the environment.

Roles of the variable P450 substrate recognition sites SRS1 and SRS6 in esfenvalerate metabolism by CYP6AE subfamily enzymes in Helicoverpa armigera

The cotton bollworm P450s of the clustered CYP6AE subfamily share high sequence identities but differ dramatically in their capacity to metabolize xenobiotics, especially esfenvalerate. Among them, CYP6AE17 has the highest sequence identity with CYP6AE18 but shows ~7-fold higher metabolic efficiency. CYP6AE11 is most active towards esfenvalerate but CYP6AE20 is inactive even though the enzymes share 54.8% sequence identity. Sequence analysis revealed the SRS1 (Substrate Recognition Site) and SRS6 between CYP6AE17 and CYP6AE18, and SRS1 between CYP6AE11 and CYP6AE20 are the most variable among all six SRSs. In order to identify the key factors that underlie the observed catalytic difference, we exchanged these SRS sequences between two pairs of P450s and studied the activity of the resulting hybrid mutants or chimeras. In vitro metabolism showed that the CYP6AE17/18 chimeras had 2- and 14-fold decreased activities and the CYP6AE18/17 chimeras had 6- and 10-fold increased activities to esfenvalerate. Meanwhile, after exchanging SRS1 with each other, the CYP6AE11/20 chimera folded incorrectly but the CYP6AE20/11 chimera gained moderate activity to esfenvalerate. Molecular modelling showed that amino acids variants within SRS1 or SRS6 change the shape and chemical environment of the active sites, which may affect the ligand-binding interactions. These results indicate that the protein structure variation resulting from the sequence diversity of SRSs promotes the evolution of insect chemical defense and contributes to the development of insect resistance to pesticides.

Arginine-259 of UGT2B7 Confers UDP-Sugar Selectivity

Enzymes of the human UDP-glycosyltransferase (UGT) superfamily typically catalyze the covalent addition of the sugar moiety from a UDP-sugar cofactor to relatively low-molecular weight lipophilic compounds. Although UDP-glucuronic acid (UDP-GlcUA) is most commonly employed as the cofactor by UGT1 and UGT2 family enzymes, UGT2B7 and several other enzymes can use both UDP-GlcUA and UDP-glucose (UDP-Glc), leading to the formation of glucuronide and glucoside conjugates. An investigation of UGT2B7-catalyzed morphine glycosidation indicated that glucuronidation is the principal route of metabolism because the binding affinity of UDP-GlcUA is higher than that of UDP-Glc. Currently, it is unclear which residues in the UGT2B7 cofactor binding domain are responsible for the preferential binding of UDP-GlcUA. Here, molecular dynamics (MD) simulations were performed together with site-directed mutagenesis and enzyme kinetic studies to identify residues within the UGT2B7 binding site responsible for the selective cofactor binding. MD simulations demonstrated that Arg259, which is located within the N-terminal domain, specifically interacts with UDP-GlcUA, whereby the side chain of Arg259 H-bonds and forms a salt bridge with the carboxylate group of glucuronic acid. Consistent with the MD simulations, substitution of Arg259 with Leu resulted in the loss of morphine, 4-methylumbelliferone, and zidovudine glucuronidation activity, but morphine glucosidation was preserved. SIGNIFICANCE STATEMENT: Despite the importance of uridine diphosphate glycosyltransferase (UGT) enzymes in drug and chemical metabolism, cofactor binding interactions are incompletely understood, as is the molecular basis for preferential glucuronidation by UGT1 and UGT2 family enzymes. The study demonstrated that long timescale molecular dynamics (MD) simulations with a UGT2B7 homology model can be used to identify critical binding interactions of a UGT protein with UDP-sugar cofactors. Further, the data provide a basis for the application of MD simulations to the elucidation of UGT-aglycone interactions.

Molecular Dynamics Investigations of Binding Mechanism for Triazoles Inhibitors to CYP51

The sterol 14α demethylase enzyme (CYP51) is an important target of fungal infections. However, the molecular mechanism between triazoles inhibitors and CYP51 remains obscure. In this study, we have investigated the binding mechanism and tunnel characteristic upon four triazoles inhibitors with CYP51 based on the molecular docking and molecular dynamics simulations. The results indicate the four inhibitors stabilize in the binding cavity of CYP51 in a similar binding mode. We discover a hydrophobic cavity (F58, Y64, Y118, L121, Y132, L376, S378, S506, S507, and M508) and the hydrophobic interaction is the main driving force for inhibitors binding to CYP51. The long-tailed inhibitors (posaconazole and itraconazole) have stronger binding affinities than short-tailed inhibitors (fluconazole and voriconazole) because long-tailed inhibitors can form more hydrophobic interactions with CYP51. The tunnel 2f is the predominant pathway for inhibitors ingress/egress protein, which is similar to the other works of CYP51. This study could provide the theoretical basis for the development of efficient azoles inhibitors and may lead a better insight into structure-function relationships of CYP51.

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.

Slow-, Tight-Binding Inhibition of CYP17A1 by Abiraterone Redefines Its Kinetic Selectivity and Dosing Regimen

Substantial evidence underscores the clinical efficacy of inhibiting CYP17A1-mediated androgen biosynthesis by abiraterone for treatment of prostate oncology. Previous structural analysis and in vitro assays revealed inconsistencies surrounding the nature and potency of CYP17A1 inhibition by abiraterone. Here, we establish that abiraterone is a slow-, tight-binding inhibitor of CYP17A1, with initial weak binding preceding the subsequent slow isomerization to a high-affinity CYP17A1-abiraterone complex. The in vitro inhibition constant of the final high-affinity CYP17A1-abiraterone complex ( ( K _i * = 0.39 nM )yielded a binding free energy of -12.8 kcal/mol that was quantitatively consistent with the in silico prediction of -14.5 kcal/mol. Prolonged suppression of dehydroepiandrosterone (DHEA) concentrations observed in VCaP cells after abiraterone washout corroborated its protracted CYP17A1 engagement. Molecular dynamics simulations illuminated potential structural determinants underlying the rapid reversible binding characterizing the two-step induced-fit model. Given the extended residence time (42 hours) of abiraterone within the CYP17A1 active site, in silico simulations demonstrated sustained target engagement even when most abiraterone has been eliminated systemically. Subsequent pharmacokinetic-pharmacodynamic (PK-PD) modeling linking time-dependent CYP17A1 occupancy to in vitro steroidogenic dynamics predicted comparable suppression of downstream DHEA-sulfate at both 1000- and 500-mg doses of abiraterone acetate. This enabled mechanistic rationalization of a clinically reported PK-PD disconnect, in which equipotent reduction of downstream plasma DHEA-sulfate levels was achieved despite a lower systemic exposure of abiraterone. Our novel findings provide the impetus for re-evaluating the current dosing paradigm of abiraterone with the aim of preserving PD efficacy while mitigating its dose-dependent adverse effects and financial burden. SIGNIFICANCE STATEMENT: With the advent of novel molecularly targeted anticancer modalities, it is becoming increasingly evident that optimal dose selection must necessarily be predicated on mechanistic characterization of the relationships between target exposure, drug-target interactions, and pharmacodynamic endpoints. Nevertheless, efficacy has always been perceived as being exclusively synonymous with affinity-based measurements of drug-target binding. This work demonstrates how elucidating the slow-, tight-binding inhibition of CYP17A1 by abiraterone via in vitro and in silico analyses was pivotal in establishing the role of kinetic selectivity in mediating time-dependent CYP17A1 engagement and eventually downstream efficacy outcomes.

Molecular Basis for Metabolic Regioselectivity and Mechanism of Cytochrome P450s toward Carcinogenic 4-(Methylnitrosamino)-(3-pyridyl)-1-butanone

As an abundantly present tobacco component, carcinogen 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone (NNK) has also been detected in atmospheric particulate matter, suggesting the ineluctable exposure risk of this contaminant. NNK metabolic activation by cytochrome P450 enzymes (CYPs) is a prerequisite to exerting its genotoxicity, but the metabolic regioselectivity and mechanism are still unknown. Here the binding feature and regioselectivity of CYPs 1A1, 1A2, 2A6, 2A13, 2B6, and 3A4 toward NNK are unraveled through molecular docking and molecular dynamics (MD) simulations. Binding mode analyses reveal that 1A2 and 2B6 have definite preferences for NNK α-methyl hydroxylation, while the other four CYPs preferentially catalyze α-methylene hydroxylation. The binding affinities between NNK and CYPs evaluated by the binding free energies follow the order 2A13 > 2B6 > 1A2 > 2A6 > 1A1 > 3A4. Density functional theory (DFT) calculations are further performed to characterize the mechanism of NNK biotransformation. Results show that the α-hydroxyNNK generated from α-hydroxylation may undergo nonenzymatic decomposition to form genotoxic diazohydroxide and aldehyde, and further oxidation by P450 to yield nitrosamide, which mainly contributes to NNK toxification capacity. Meanwhile the pyridine N-oxidation and denitrosation of C_α-radical intermediate play an important role in detoxifying NNK. Overall, the present study provides the molecular basis for CYP-catalyzed regioselectivity and mechanism of NNK biotransformation, which can enable the identification of metabolites for assessing the health risk of individual NNK exposure.

How do mutations affect the structural characteristics and substrate binding of CYP21A2? An investigation by molecular dynamics simulations

CYP21A2 mutations affect the activity of the protein leading to CAH disease.

Influence of Selected Carbon Nanostructures on the CYP2C9 Enzyme of the P450 Cytochrome

Carbon nanostructures have recently gained significant interest from scientists due to their unique physicochemical properties and low toxicity. They can accumulate in the liver, which is the main expression site of cytochrome P450 (CYP450) enzymes. These enzymes play an important role in the metabolism of exogenous compounds, such as drugs and xenobiotics. Altered activity or expression of CYP450 enzymes may lead to adverse drug effects and toxicity. The objective of this study was to evaluate the influence of three carbon nanostructures on the activity and expression at the mRNA and protein levels of CYP2C9 isoenzyme from the CYP2C subfamily: Diamond nanoparticles, graphite nanoparticles, and graphene oxide platelets. The experiments were conducted using two in vitro models. A microsome model was used to assess the influence of the three-carbon nanostructures on the activity of the CYP2C9 isoenzyme. The CYP2C9 gene expression at the mRNA and protein levels was determined using a hepatoma-derived cell line HepG2. The experiments have shown that all examined nanostructures inhibit the enzymatic activity of the studied isoenzymes. Moreover, a decrease in the expression at the mRNA and protein levels was also observed. This indicates that despite low toxicity, the nanostructures can alter the enzymatic function of CYP450 enzymes, and the molecular pathways involved in their expression.

Evaluation of the time-dependent antiproliferative activity and liver microsome stability of 3 phenyl 4-(2-oxo-3-alkylimidazolidin-1-yl)benzenesulfonates as promising CYP1A1-dependent antimicrotubule prodrugs

OBJECTIVES

In this study, the antiproliferative activity of 3 phenyl 4-(2-oxo-3-alkylimidazolidin-1-yl)benzenesulfonates (PAIB-SOs) was assessed in a time-dependent manner together with their hepatic stability and metabolism using human, mouse and rat liver microsomes.

METHODS

CEU-818, -820 and -913 were selected as promising hit compounds. Their antiproliferative activity on human breast carcinoma MCF-7 cells was evaluated using escalating concentrations of drugs at 24, 36 and 48 h and the sulforhodamine B assay. Their hepatic stability was evaluated by HPLC-UV of extracts obtained from human, mouse and rat liver microsomes.

KEY FINDINGS

The antiproliferative activity of PAIB-SOs is concentration and time-dependent and requires between 24 and 36 h of contact with MCF-7 cells to detect a significant antiproliferative activity. PAIB-SOs stability in microsomes usually decreases following this order: human ≈ (rat > mouse). The CEU-913 exhibits the longest half-life in rat and human liver microsomes while the CEU-820 exhibits the longest half-life in mouse liver microsomes.

CONCLUSIONS

Our in vitro results suggest that PAIB-SOs should have a minimum contact time of 24 h with the tumour to trigger significant antitumoural activity. The activity of mouse liver microsomes towards PAIB-SOs is higher than rat microsomes and tends to be higher than human liver microsomes.

Transcriptome-Based Identification and Molecular Evolution of the Cytochrome P450 Genes and Expression Profiling under Dimethoate Treatment in Amur Stickleback (Pungitius sinensis)

Cytochrome P450s (CYPs) are a family of membrane-bound mono-oxygenase proteins, which are involved in cell metabolism and detoxification of various xenobiotic substances. In this study, we identified 58 putative CYP genes in Amur stickleback (Pungitius sinensis) based on the transcriptome sequencing. Conserved motif distribution suggested their functional relevance within each group. Some present recombination events have accelerated the evolution of this gene family. Moreover, a few positive selection sites were identified, which may have accelerated the functional divergence of this family of proteins. Expression patterns of these CYP genes were investigated and indicated that most were affected by dimethoate treatment, suggesting that CYPs were involved in the detoxication of dimethoate. This study will provide a foundation for the further functional investigation of CYP genes in fishes.

Inhibition of Hepatic CYP2D6 by the Active N-Oxide Metabolite of Sorafenib

The multikinase inhibitor sorafenib (SOR) is used to treat patients with hepatocellular and renal carcinomas. SOR undergoes CYP-mediated biotransformation to a pharmacologically active N-oxide metabolite (SNO) that has been shown to accumulate to varying extents in individuals. Kinase inhibitors like SOR are frequently coadministered with a range of other drugs to improve the efficacy of anticancer drug therapy and to treat comorbidities. Recent evidence has suggested that SNO is more effective than SOR as an inhibitor of CYP3A4-mediated midazolam 1'-hydroxylation. CYP2D6 is also reportedly inhibited by SOR. The present study assessed the possibility that SNO might contribute to CYP2D6 inhibition. The inhibition kinetics of CYP2D6-mediated dextromethorphan O-demethylation were analyzed in human hepatic microsomes, with SNO found to be ~ 19-fold more active than SOR (K_is 1.8 ± 0.3 μM and 34 ± 11 μM, respectively). Molecular docking studies of SOR and SNO were undertaken using multiple crystal structures of CYP2D6. Both molecules mediated interactions with key amino acid residues in putative substrate recognition sites of CYP2D6. However, a larger number of H-bonding interactions was noted between the N-oxide moiety of SNO and active site residues that account for its greater inhibition potency. These findings suggest that SNO has the potential to contribute to pharmacokinetic interactions involving SOR, perhaps in those individuals in whom SNO accumulates.

P450 Monooxygenases Enable Rapid Late-Stage Diversification of Natural Products via C-H Bond Activation

The biological potency of natural products has been exploited for decades. Their inherent structural complexity and natural diversity might hold the key to efficiently address the urgent need for the development of novel pharmaceuticals. At the same time, it is that very complexity, which impedes necessary chemical modifications such as structural diversification, to improve the effectiveness of the drug. For this purpose, Cytochrome P450 enzymes, which possess unique abilities to activate inert sp3-hybridised C-H bonds in a late-stage fashion, offer an attractive synthetic tool. In this review the potential of cytochrome P450 enzymes in chemoenzymatic lead diversification is illustrated discussing studies reporting late-stage functionalisations of natural products and other high-value compounds. These enzymes were proven to extend the synthetic toolbox significantly by adding to the flexibility and efficacy of synthetic strategies of natural product chemists, and scientists of other related disciplines.

Can multiscale simulations unravel the function of metallo-enzymes to improve knowledge-based drug discovery?

Metallo-enzymes are a large class of biomolecules promoting specialized chemical reactions. Quantum-classical quantum mechanics/molecular mechanics molecular dynamics, describing the metal site at quantum mechanics level, while accounting for the rest of system at molecular mechanics level, has an accessible time-scale limited by its computational cost. Hence, it must be integrated with classical molecular dynamics and enhanced sampling simulations to disentangle the functions of metallo-enzymes. In this review, we provide an overview of these computational methods and their capabilities. In particular, we will focus on some systems such as CYP19A1 a Fe-dependent enzyme involved in estrogen biosynthesis, and on Mg2+-dependent DNA/RNA processing enzymes/ribozymes and the spliceosome, a protein-directed ribozyme. This information may guide the discovery of drug-like molecules and genetic manipulation tools.

Spin Interconversion of Heme-Peroxo-Copper Complexes Facilitated by Intramolecular Hydrogen-Bonding Interactions

Synthetic peroxo-bridged high-spin (HS) heme-(μ-η2:η1-O22-)-Cu(L) complexes incorporating (as part of the copper ligand) intramolecular hydrogen-bond (H-bond) capabilities and/or steric effects are herein demonstrated to affect the complex's electronic and geometric structure, notably impacting the spin state. An H-bonding interaction with the peroxo core favors a low-spin (LS) heme-(μ-η1:η1-O22-)-Cu(L) structure, resulting in a reversible temperature-dependent interconversion of spin state (5 coordinate HS to 6 coordinate LS). The LS state dominates at low temperatures, even in the absence of a strong trans-axial heme ligand. Lewis base addition inhibits the H-bond facilitated spin interconversion by competition for the H-bond donor, illustrating the precise H-bonding interaction required to induce spin-crossover (SCO). Resonance Raman spectroscopy (rR) shows that the H-bonding pendant interacts with the bridging peroxide ligand to stabilize the LS but not the HS state. The H-bond (to the Cu-bound O atom) acts to weaken the O-O bond and strengthen the Fe-O bond, exhibiting ν(M-O) and ν(O-O) values comparable to analogous known LS complexes with a strong donating trans-axial ligand, 1,5-dicyclohexylimidazole, (DCHIm)heme-(μ-η1:η1-O22-)-Cu(L). Variable-temperature (-90 to -130 °C) UV-vis and 2H NMR spectroscopies confirm the SCO process and implicate the involvement of solvent binding. Examining a case of solvent binding without SCO, thermodynamic parameters were obtained from a van't Hoff analysis, accounting for its contribution in SCO. Taken together, these data provide evidence for the H-bond group facilitating a core geometry change and allowing solvent to bind, stabilizing a LS state. The rR data, complemented by DFT analysis, reveal a stronger H-bonding interaction with the peroxo core in the LS compared to the HS complexes, which enthalpically favors the LS state. These insights enhance our fundamental understanding of secondary coordination sphere influences in metalloenzymes.

Dynamic Effects of CYP2D6 Genetic Variants in a Set of Poor Metaboliser Patients with Infiltrating Ductal Cancer Under Treatment with Tamoxifen

Breast cancer is a group of multigenic diseases. It is the most common cancer diagnosed among women worldwide and is often treated with tamoxifen. Tamoxifen is catalysed by cytochrome P450 2D6 (CYP2D6), and inter-individual variations in the enzyme due to single nucleotide polymorphisms (SNPs) could alter enzyme activity. We evaluated SNPs in patients from Colombia in South America who were receiving tamoxifen treatment for breast cancer. Allelic diversity in the CYP2D6 gene was found in the studied population, with two patients displaying the poor-metaboliser phenotype. Molecular dynamics and trajectory analyses were performed for CYP2D6 from these two patients, comparing it with the common allelic form (CYP2D6*1). Although we found no significant structural change in the protein, its dynamics differ significantly from those of CYP2D6*1, the effect of such differential dynamics resulting in an inefficient enzyme with serious implications for tamoxifen-treated patients, increasing the risk of disease relapse and ineffective treatment.

From Endogenous Compounds as Biomarkers to Plasma-Derived Nanovesicles as Liquid Biopsy; Has the Golden Age of Translational Pharmacokinetics-Absorption, Distribution, Metabolism, Excretion-Drug-Drug Interaction Science Finally Arrived?

It is now established that a drug's pharmacokinetics (PK) absorption, distribution, metabolism, excretion (ADME) and drug-drug interaction (DDI) profile can be modulated by age, disease, and genotype. In order to facilitate subject phenotyping and clinical DDI assessment, therefore, various endogenous compounds (in plasma and urine) have been pursued as drug-metabolizing enzyme and transporter biomarkers. Compared with biomarkers, however, the topic of circulating extracellular vesicles as "liquid biopsy" has received little attention within the ADME community; most organs secrete nanovesicles (e.g., exosomes) into the blood that contain luminal "cargo" derived from the originating organ (proteins, messenger RNA, and microRNA). As such, ADME profiling of plasma exosomes could be leveraged to better define genotype-phenotype relationships and the study of ontogeny, disease, and complex DDIs. If methods to support the isolation of tissue-derived plasma exosomes are successfully developed and validated, it is envisioned that they will be used jointly with genotyping, biomarkers, and modeling tools to greatly progress translational PK-ADME-DDI science.