Crystal Structure of Human Cytochrome P450 2D6

Understanding Cytochrome P450 Enzyme Substrate Inhibition and Prospects for Elimination Strategies

Cytochrome P450 (CYP450) enzymes, which are widely distributed and pivotal in various biochemical reactions, catalyze diverse processes such as hydroxylation, epoxidation, dehydrogenation, dealkylation, nitrification, and bond formation. These enzymes have been applied in drug metabolism, antibiotic production, bioremediation, and fine chemical synthesis. Recent research revealed that CYP450 catalytic kinetics deviated from the classic Michaelis-Menten model. A notable substrate inhibition phenomenon that affects the catalytic efficiency of CYP450 at high substrate concentrations was identified. However, the substrate inhibition of various reactions catalyzed by CYP450 enzymes have not been comprehensively reviewed. This review describes CYP450 substrate inhibition examples and atypical Michaelis-Menten kinetic models, and provides insight into mechanisms of these enzymes. We also reviewed 3D structure and dynamics of CYP450 with substrate binding. Outline methods for alleviating substrate inhibition in CYP450 and other enzymes, including traditional fermentation approaches and protein engineering modifications. The comprehensive analysis presented in this study lays the foundation for enhancing the catalytic efficiency of CYP450 by deregulating substrate inhibition.

Roles of Individual Human Cytochrome P450 Enzymes in Drug Metabolism

Our knowledge of the roles of individual cytochrome P450 (P450) enzymes in drug metabolism has developed considerably in the past 30 years, and this base has been of considerable use in avoiding serious issues with drug interactions and issues due to variations. Some newer approaches are being considered for "phenotyping" metabolism reactions with new drug candidates. Endogenous biomarkers are being used for noninvasive estimation of levels of individual P450 enzymes. There is also the matter of some remaining "orphan" P450s, which have yet to be assigned reactions. Practical problems that continue in drug development include predicting drug-drug interactions, predicting the effects of polymorphic and other P450 variations, and evaluating interspecies differences in drug metabolism, particularly in the context of "metabolism in safety testing" regulatory issues ["disproportionate (human) metabolites"]. SIGNIFICANCE STATEMENT: Cytochrome P450 enzymes are the major catalysts involved in drug metabolism. The characterization of their individual roles has major implications in drug development and clinical practice.

Rational Molecular Design of a Fluorescent Probe for Selectively Sensing Human Cytochrome P450 2D6

Cytochrome P450 2D6 (CYP2D6) is a key enzyme that mediates the metabolism of various drugs and endogenous substances in humans. However, its biological role in drug-drug interactions especially mechanism-based inactivation (MBI), and various diseases remains poorly understood, owing to the lack of molecular tools suitable for selectively monitoring CYP2D6 in complex biological systems. Herein, using a tailored molecular strategy, we developed a fluorescent probe BDPM for CYP2D6. BDPM exhibits excellent specificity and imaging capability for CYP2D6, making it suitable for the real-time monitoring of endogenous CYP2D6 activity in living bio-samples. Therefore, our tailored strategy proved useful for constructing the highly selective and enzyme-activated fluorescent probes. BDPM as a molecular tool to explore the critical roles of CYP2D6 in the pathogenesis of diseases, high-throughput screening of inhibitors and intensive investigation of CYP2D6-induced MBI in natural systems.

Aromatase inhibitors for the treatment of breast cancer: An overview (2019-2023)

Aromatase inhibition is considered a legitimate approach for the treatment of ER-positive (ER+) breast cancer as it accounts for more than 70% of breast cancer cases. Aromatase inhibitor therapy has been demonstrated to be highly effective in decreasing tumour size, increasing survival rates, and lowering the chance of cancer recurrence. The present review deliberates the pathophysiology and the role of aromatase in estrogen biosynthesis. Estrogen biosynthesis, various androgens, and their function in the human body have also been discussed. The salient aspects of the aromatase active site, its mode of action, and AIs, along with their intended interactions with presently FDA-approved inhibitors, have been briefly discussed. It has been detailed how different reported AIs were designed, their SAR investigations, in silico analysis, and biological evaluations. Various AIs from multiple origins, such as synthetic and semi-synthetic, have also been discussed.

Effect of Lipid Bilayer Anchoring on the Conformational Properties of the Cytochrome P450 2D6 Binding Site

Human cytochrome P450 (CYP) proteins metabolize 75% of small-molecule pharmaceuticals, which makes structure-based modeling of CYP metabolism and inhibition, bolstered by the current availability of X-ray crystal structures of CYP globular catalytic domains, an attractive prospect. Accounting for this broad metabolic capacity is a combination of the existence of multiple different CYP proteins and the capacity of a single CYP protein to metabolize multiple different small molecules. It is thought that structural plasticity and flexibility contribute to this latter property; therefore, incorporating diverse conformational states of a particular CYP is likely an important consideration in structure-based CYP metabolism and inhibition modeling. While all-atom explicit-solvent molecular dynamics simulations can be used to generate conformational ensembles under biologically relevant conditions, existing CYP crystal structures are of the globular domain only, whereas human CYPs contain N-terminal transmembrane and linker peptides that anchor the globular catalytic domain to the endoplasmic reticulum. To determine whether this can cause significant differences in the sampled binding site conformations, microsecond scale all-atom explicit-solvent molecular dynamics simulations of the CYP2D6 globular domain in an aqueous environment were compared with those of the full-length protein anchored in a POPC lipid bilayer. While bilayer-anchoring damped some structural fluctuations in the globular domain relative to the aqueous simulations, none of the affected residues included binding site pocket residues. Furthermore, clustering of molecular dynamics snapshots based on either pairwise binding site pocket RMSD or volume differences demonstrated a lack of separation of snapshots from the two simulation conditions into different clusters. These results suggest the substantially simpler and computationally cheaper aqueous simulation approach can be used to generate a relevant conformational ensemble of the CYP2D6 binding site for structure-based metabolism and inhibition modeling.

The genomic landscape of CYP2D6 variation in the Indian population

Aim: The CYP2D6 gene is highly polymorphic, causing large interindividual variability in the metabolism of several clinically important drugs. Materials & methods: The authors investigated the diversity and distribution of CYP2D6 alleles in Indians using whole genome sequences (N = 1518). Functional consequences were assessed using pathogenicity scores and molecular dynamics simulations. Results: The analysis revealed population-specific CYP2D6 alleles (*86, *7, *111, *112, *113, *99) and remarkable differences in variant and phenotype frequencies with global populations. The authors observed that one in three Indians could benefit from a dose alteration for psychiatric drugs with accurate CYP2D6 phenotyping. Molecular dynamics simulations revealed large conformational fluctuations, confirming the predicted reduced function of *86 and *113 alleles. Conclusion: The findings emphasize the utility of comprehensive CYP2D6 profiling for aiding precision public health.

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.

Interactions studies of CYP2D6 with quercetin and hyperoside by spectral analysis and molecular dynamics simulations

Some ingredients from herbal medicine can significantly affect the activity of CYP2D6, thus leading to serious interactions between herbs and drugs. Quercetin and hyperoside are active ingredients widely found in vegetables, fruits, and herbal medicines. Quercetin and hyperoside have many biological activities. In this work, the characteristic bindings of CYP2D6 with quercetin/hyperoside are revealed by multi-spectroscopy analysis, molecular docking, and molecular dynamics simulations. The fluorescence of CYP2D6 is statically quenched by quercetin and hyperoside. The binding constant (Ka ) values of CYP2D6-quercetin/hyperoside range from 104  L mol-1 , which indicates that these two flavonoids bind moderately to CYP2D6. Meanwhile, quercetin has a stronger quenching ability to CYP2D6 than that of hyperoside. The secondary structure of CYP2D6 is obviously changed by binding with quercetin/hyperoside. The docking results reveal that the quercetin/hyperoside enters the active site of CYP2D6 near heme and binds to CYP2D6 by hydrogen bonds and van der Waals forces. The molecular dynamics simulation results indicate that the binding of quercetin/hyperoside can stabilize the two complexes, enhance the flexibility of CYP2D6 backbone atoms, and make a more unfolded and looser structure of CYP2D6.

Structural and vibrational investigations and molecular docking studies of a vinca alkoloid, vinorelbine

Vinorelbine, a vinca alkaloid, is an antimitotic drug that inhibits polymerisation process of tubulins to microtubules, and is widely used in cancer chemotherapy. Due to the importance of the structure-activity relationship, in this work the conformational preferences of the vinorelbine molecule were surched by PM3 method. The obtained lowest energy conformer was then optimized at DFT/B3LYP/6-31G(d,p) level of theory and the structural characteristics were determined. Frontier orbital (HOMO, LUMO) and molecular electrostatic potential (MEP) analyses were performed for the optimized structure. The experimental FT-IR, Raman and UV-VIS spectral data of vinorelbine along with the theoretical DFT/B3LYP/6-31G(d,p) calculations were investigated in detail. The vibrational wavenumbers were assigned based on the calculated potential energy distribution (PED) of the vibrational modes. To shed light into the anticancer property of vinorelbine as microtubule destabilizer, the most favourable binding mode and the interaction details between vinorelbine and tubulin were revealed by molecular docking studies of vinorelbine into the α,β-tubulin (PDB IDs: 4O2B; 1SA0; 7CNN) and binding free energies were calculated by the combination of Molecular Mechanics/Generalized Born Surface Area (MMGBSA) and Molecular Mechanics/Poisson-Boltzmann Surface Area (MM-PBSA) methods {MM/PB(GB)SA}. The calculated vinorelbine-7CNN binding free energy, using by MM/PB(GB)SA approach, was found to be the best (-50.39 kcal/mol), and followed by vinorelbine-4O2B (-28.5 kcal/mol) and vinorelbine-1SA0 (-17.59 kcal/mol) systems. Moreover, the interaction of vinorelbine with the cytochrome P450 enzymes (CYP), which are known to help in the metabolism of many drugs in the body, was investigated by docking studies against CYP2D6 and CYP3A4 targets.Communicated by Ramaswamy H. Sarma.

Catalytic Site Constraints in the P450s Mediating Loganic Acid (7DLH) and Secologanic Acid Synthesis (SLAS) in Camptotheca

Terpene indole alkaloids (TIAs) are plant-derived natural products synthesized in low levels in medicinal plants such as Catharanthus roseus and Camptotheca acuminata. TIA pathways species utilize several CYP72A subfamily members to form loganic acid from 7-deoxyloganic acid (a simple hydroxylation) as well as secologanin and secologanic acid from loganin and loganic acid (a C-C bond scission). Divergences in the specificities of these P450s have allowed Camptotheca secologanic acid synthases (SLASs) to become bifunctional enzymes capable of performing both reactions. In contrast, Catharanthus 7-deoxyloganic acid hydroxylase (7DLH) and secologanin synthase (SLS) have remained monofunctional enzymes capable either of monooxygenation or C-C bond scission. Our in vitro reconstitutions have now demonstrated that Camptotheca also contains a monofunctional 7DLH capable only of hydroxylating 7-deoxyloganic acid. Mutageneses aimed at evaluating residues important for the tight specificity of Camptotheca 7DLH (CYP72A729) and the broad specificity of SLAS (CYP72A564) have identified several residues where reciprocal switches substantially affect their activities: Lys128His in 7DLH increases hydroxylation of 7-deoxyloganic acid, and His132Lys in SLAS decreases this hydroxylation and C-C bond scissions of loganic acid and loganin; Gly321Ser in 7DLH does not affect hydroxylation of 7-deoxyloganic acid, whereas Ser324Gly in SLAS significantly increases C-C bond scission of loganic acid; Asp332Glu in the acid-alcohol pair of 7DLH increases hydroxylation of 7-deoxyloganic acid, whereas Glu335Asp in SLAS completely eliminates both of its activities. These mutations that enhance or eliminate these respective activities have significant potential to aid engineering efforts aimed at increasing TIA production in cell cultures, microbial systems, and/or other plants.

Selective steroidogenic cytochrome P450 haem iron ligation by steroid-derived isonitriles

Alkyl isonitriles, R-NC, have previously been shown to ligate the heme (haem) iron of cytochromes P450 in both accessible oxidation states (ferrous, Fe2+, and ferric, Fe3+). Herein, the preparation of four steroid-derived isonitriles and their interactions with several P450s, including the steroidogenic CYP17A1 and CYP106A2, as well as the more promiscuous drug metabolizers CYP3A4 and CYP2D6, is described. It was found that successful ligation of the heme iron by the isonitrile functionality for a given P450 depends on both the position and stereochemistry of the isonitrile on the steroid skeleton. Spectral studies indicate that isonitrile ligation of the ferric heme is stable upon reduction to the ferrous form, with reoxidation resulting in the original complex. A crystallographic structure of CYP17A1 with an isonitrile derived from pregnanalone further confirmed the interaction and identified the absolute stereochemistry of the bound species.

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.

Recombinant Technologies Facilitate Drug Metabolism, Pharmacokinetics, and General Biomedical Research

The development of safe and effective medications requires a profound understanding of their pharmacokinetic (PK) and pharmacodynamic properties. PK studies have been built through investigation of enzymes and transporters that drive drug absorption, distribution, metabolism, and excretion (ADME). Like many other disciplines, the study of ADME gene products and their functions has been revolutionized through the invention and widespread adoption of recombinant DNA technologies. Recombinant DNA technologies use expression vectors such as plasmids to achieve heterologous expression of a desired transgene in a specified host organism. This has enabled the purification of recombinant ADME gene products for functional and structural characterization, allowing investigators to elucidate their roles in drug metabolism and disposition. This strategy has also been used to offer recombinant or bioengineered RNA (BioRNA) agents to investigate the posttranscriptional regulation of ADME genes. Conventional research with small noncoding RNAs such as microRNAs (miRNAs) and small interfering RNAs has been dependent on synthetic RNA analogs that are known to carry a range of chemical modifications expected to improve stability and PK properties. Indeed, a novel transfer RNA fused pre-miRNA carrier-based bioengineering platform technology has been established to offer consistent and high-yield production of unparalleled BioRNA molecules from Escherichia coli fermentation. These BioRNAs are produced and processed inside living cells to better recapitulate the properties of natural RNAs, representing superior research tools to investigate regulatory mechanisms behind ADME. SIGNIFICANCE STATEMENT: This review article summarizes recombinant DNA technologies that have been an incredible boon in the study of drug metabolism and PK, providing investigators with powerful tools to express nearly any ADME gene products for functional and structural studies. It further overviews novel recombinant RNA technologies and discusses the utilities of bioengineered RNA agents for the investigation of ADME gene regulation and general biomedical research.

Development of an efficient insecticide substrate and inhibitor screening system of insect P450s using fission yeast

Metabolic resistance is one of the most frequent mechanisms of insecticide resistance, characterized by an increased expression of several important enzymes and transporters, especially cytochrome P450s (CYPs). Due to the large number of P450s in pests, determining the precise relationship between these enzymes and the insecticide substrates is a challenge. Herein, we developed a luminescence-based screening system for efficient identification of insecticide substrates and insect P450 inhibitors. We recombinantly expressed Bemisia tabaci CYP6CM1vQ (Bt CYP6CM1vQ) in the fission yeast Schizosaccharomyces pombe and subsequently permeabilized the yeast cells to convert them into "enzyme bags". We exploited these enzyme bags to screen the activity of twelve luciferin substrates and identified Luciferin-FEE as the optimal competing probe that was further used to characterize the metabolism of eight candidate commercial insecticides. Among them, Bt CYP6CM1vQ exhibited notable activity against pymetrozine and imidacloprid. Their binding modes were predicted by homology modeling and molecular docking, revealing the mechanisms of the metabolism. We also tested the inhibitory effect of eight known P450 inhibitors using our system and identified letrozole and 1-benzylimidazole as showing significant activity against Bt CYP6CM1vQ, with IC₅₀ values of 23.74 μM and 1.30 μM, respectively. Their potential to be developed as an insecticide synergist was further proven by an in vitro toxicity assay using imidacloprid-resistant Bemisia tabaci. Overall, our luciferin-based enzyme bag method is capable of providing a robust and efficient screening of insect P450 substrates and, more importantly, inhibitors to overcome the resistance.

Key regulators in the architecture of substrate access/egress channels in mammalian cytochromes P450 governing flexibility in substrate oxyfunctionalization

Cytochrome P450s (CYP) represent a superfamily of b-type hemoproteins catalyzing oxifunctionalization of a vast array of endogenous and exogenous compounds. The present review focuses on assessment of the topology of prospective determinants in substrate entry and product release channels of mammalian P450s, steering the conformational dynamics of substrate accessibility and productive ligand orientation toward the iron-oxene core. Based on a generalized, CYP3A4-related construct, the sum of critical elements from diverse target enzymes was found to cluster within the known substrate recognition sites. The majority of prevalent substrate access/egress tunnels revealed to be of fairly balanced functional importance. The hydrophobicity profile of the candidates revealed to be the most salient feature in functional interaction throughout the conduits, while bulkiness of the residues imposes steric restrictions on substrate traveling. Thus, small amino acids such as prolines and glycines serve as hinges, driving conformational flexibility in ligand passage. Similarly, bottlenecks in the tunnel architecture, being narrowest encounter points within the CYP3A4 model, have a vital function in substrate selectivity along with clusters of aromatic amino acids acting as gatekeepers. In addition, peripheral patches in conduits may house determinants modulating allosteric cooperativity between remote and central domains in the P450 structure. Remarkably, the bulk critical residues lining tunnels in the various isozymes reside in helices B'/C and F/G inclusive of their interhelical turns as well as in helix I. This suggests these regions to represent hotspots for targeted genetic engineering to tailor more sophisticated mammalian P450s exploitable in industrial, biotechnological and medicinal areas.

Combined and independent effects of OCT1 and CYP2D6 on the cellular disposition of drugs

The organic cation transporter 1 (OCT1) mediates the cell uptake and cytochrome P450 2D6 (CYP2D6) the metabolism of many cationic substrates. Activities of OCT1 and CYP2D6 are affected by enormous genetic variation and frequent drug-drug interactions. Single or combined deficiency of OCT1 and CYP2D6 might result in dramatic differences in systemic exposure, adverse drug reactions, and efficacy. Thus, one should know what drugs are affected to what extent by OCT1, CYP2D6 or both. Here, we compiled all data on CYP2D6 and OCT1 drug substrates. Among 246 CYP2D6 substrates and 132 OCT1 substrates, we identified 31 shared substrates. In OCT1 and CYP2D6 single and double-transfected cells, we studied which, OCT1 or CYP2D6, is more critical for a given drug and whether there are additive, antagonistic or synergistic effects. In general, OCT1 substrates were more hydrophilic than CYP2D6 substrates and smaller in size. Inhibition studies showed unexpectedly pronounced inhibition of substrate depletion by shared OCT1/CYP2D6 inhibitors. In conclusion, there is a distinct overlap in the OCT1/CYP2D6 substrate and inhibitor spectra, so in vivo pharmacokinetics and -dynamics of shared substrates may be significantly affected by frequent OCT1- and CYP2D6-polymorphisms and by comedication with shared inhibitors.

Acetylcholinesterase inhibition of Alzheimer's disease: identification of potential phytochemicals and designing more effective derivatives to manage disease condition

Alzheimer's disease (AD) is a brain disease characterized by gradual memory loss and cognitive impairments. Acetylcholinesterase (AChE) inhibitors-such as donepezil, memantine, and tacrine-are FDA-approved medications for AD treatment. Due to the lack of their efficacy and higher side effects, many researchers have been searching for effective and safer alternatives. In this study, experimentally proved phytochemicals against brain diseases were screened based on their binding energies to the target site of AChE, pharmacokinetic properties, and drug-likeness. Although some phytochemicals showed higher binding affinities than the control drug (donepezil), they did not show permeability across the blood-brain barrier (BBB). However, berberine, anthocyanin, and diterpene alkaloid can cross the BBB and showed good binding affinities of -10.3, -10.1, and -10.2 kcal/mol, respectively. MD simulation and PCA of the simulation data of the protein and protein-ligand complexes proved that the complexes are stable in the biological environment. A total of 16 derivatives of berberine and 3 derivatives of anthocyanin also showed higher binding energies compared to the binding affinity (-11.5 kcal/mol) of the donepezil. The derivatives were designed by substituting -F, -CF3, -CN, and -NH2, and provided higher docking scores due to increasing of nonbonding interactions. MM/GBSA calculations show that the binding free energies of the best predicted derivatives of diterpene alkaloid, anthocyanin, and berberine (DA22, AC11, and BB40) are -100.4 ± 8.4, -79.3 ± 8.7, and -78.3 ± 10.7 kcal/mol, respectively, with the protein. Overall, this study was successful in finding new, highly effective, and possibly safer inhibitors of AChE.Communicated by Ramaswamy H. Sarma.

Evaluation of tea (Camellia sinensis L.) phytochemicals as multi-disease modulators, a multidimensional in silico strategy with the combinations of network pharmacology, pharmacophore analysis, statistics and molecular docking

Tea (Camellia sinensis L.) is considered as to be one of the most consumed beverages globally and a reservoir of phytochemicals with immense health benefits. Despite numerous advantages, tea compounds lack a robust multi-disease target study. In this work, we presented a unique in silico approach consisting of molecular docking, multivariate statistics, pharmacophore analysis, and network pharmacology approaches. Eight tea phytochemicals were identified through literature mining, namely gallic acid, catechin, epigallocatechin gallate, epicatechin, epicatechin gallate (ECG), quercetin, kaempferol, and ellagic acid, based on their richness in tea leaves. Further, exploration of databases revealed 30 target proteins related to the pharmacological properties of tea compounds and multiple associated diseases. Molecular docking experiment with eight tea compounds and all 30 proteins revealed that except gallic acid all other seven phytochemicals had potential inhibitory activities against these targets. The docking experiment was validated by comparing the binding affinities (Kcal mol^-1) of the compounds with known drug molecules for the respective proteins. Further, with the aid of the application of statistical tools (principal component analysis and clustering), we identified two major clusters of phytochemicals based on their chemical properties and docking scores (Kcal mol^-1). Pharmacophore analysis of these clusters revealed the functional descriptors of phytochemicals, related to the ligand-protein docking interactions. Tripartite network was constructed based on the docking scores, and it consisted of seven tea phytochemicals (gallic acid was excluded) targeting five proteins and ten associated diseases. Epicatechin gallate (ECG)-hepatocyte growth factor receptor (PDB id 1FYR) complex was found to be highest in docking performance (10 kcal mol^-1). Finally, molecular dynamic simulation showed that ECG-1FYR could make a stable complex in the near-native physiological condition.

Strategies for Assessing Acceptable Intakes for Novel N-Nitrosamines Derived from Active Pharmaceutical Ingredients

The detection of N-nitrosamines, derived from solvents and reagents and, on occasion, the active pharmaceutical ingredient (API) at higher than acceptable levels in drug products, has led regulators to request a detailed review for their presence in all medicinal products. In the absence of rodent carcinogenicity data for novel N-nitrosamines derived from amine-containing APIs, a conservative class limit of 18 ng/day (based on the most carcinogenic N-nitrosamines) or the derivation of acceptable intakes (AIs) using structurally related surrogates with robust rodent carcinogenicity data is recommended. The guidance has implications for the pharmaceutical industry given the vast number of marketed amine-containing drugs. In this perspective, the rate-limiting step in N-nitrosamine carcinogenicity, involving cytochrome P450-mediated α-carbon hydroxylation to yield DNA-reactive diazonium or carbonium ion intermediates, is discussed with reference to the selection of read-across analogs to derive AIs. Risk-mitigation strategies for managing putative N-nitrosamines in the preclinical discovery setting are also presented.

Generation of new inhibitors of selected cytochrome P450 subtypes- In silico study

Physicochemical and pharmacokinetic compound profile has crucial impact on compound potency to become a future drug. Ligands with desired activity profile cannot be used for treatment if they are characterized by unfavourable physicochemical or ADMET properties. In the study, we consider metabolic stability and focus on selected subtypes of cytochrome P450 - proteins, which take part in the first phase of compound transformations in the organism. We develop a protocol for generation of new potential inhibitors of selected cytochrome isoforms. Its subsequent stages are composed of generation and assessment of new derivatives of known cytochrome inhibitors, docking and evaluation of the compound possible inhibition on the basis of the obtained ligand-protein complexes. Besides the library of new potential agents inhibiting particular cytochrome subtypes, we also prepare a graph neural network that predicts the change in activity for all modifications of the starting molecule. In addition, we perform a systematic statistical study on the influence of particular substitutions on the potential inhibition properties of generated compounds (both mono- and di-substitutions are considered), provide explanations of the inhibitory predictions and prepare an on-line visualization platform enabling manual inspection of the results. The developed methodology can greatly support the design of new cytochrome P450 inhibitors with the overarching goal of generation of new metabolically stable compounds. It enables instant evaluation of possible compound-cytochrome interactions and selection of ligands with the highest potential of possessing desired biological activity.

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.

Single mutations toggle the substrate selectivity of multifunctional Camptotheca secologanic acid synthases (CYP72As)

Terpene indole alkaloids (TIAs) are plant-derived specialized metabolites with widespread use in medicine. Species-specific pathways derive various TIAs from common intermediates, strictosidine or strictosidinic acid, produced by coupling tryptamine with secologanin or secologanic acid. The penultimate reaction in this pathway is catalyzed by either secologanin synthase (SLS) or secologanic acid synthase (SLAS) according to whether plants produce secologanin from loganin or secologanic acid from loganic acid. Previous work has identified SLSs and SLASs from different species, but the determinants of selectivity remain unclear. Here, combining molecular modeling, ancestral sequence reconstruction, and biochemical methodologies, we identified key residues that toggle SLS and SLAS selectivity in two CYP72A (cytochrome P450) subfamily enzymes from Camptotheca acuminata. We found that the positions of foremost importance are in substrate recognition sequence 1 (SRS1), where mutations to either of two adjacent histidine residues switched selectivity; His131Phe selects for and increases secologanin production whereas His132Asp selects for secologanic acid production. Furthermore, a change in SRS3 in the predicted substrate entry channel (Arg/Lys270Thr) and another in SRS4 at the start of the I-helix (Ser324Glu) decreased enzyme activity toward either substrate. We propose that the Camptotheca SLASs have maintained the broadened activities found in a common asterid ancestor, even as the Camptotheca lineage lost its ability to produce loganin while the campanulid and lamiid lineages specialized to produce secologanin by acquiring mutations in SRS1. The identification here of the residues essential for the broad substrate scope of SLASs presents opportunities for more tailored heterologous production of TIAs.

Molecular physiology of pumiliotoxin sequestration in a poison frog

Poison frogs bioaccumulate alkaloids for chemical defense from their arthropod diet. Although many alkaloids are accumulated without modification, some poison frog species can metabolize pumiliotoxin (PTX 251D) into the more potent allopumiliotoxin (aPTX 267A). Despite extensive research characterizing the chemical arsenal of poison frogs, the physiological mechanisms involved in the sequestration and metabolism of individual alkaloids remain unclear. We first performed a feeding experiment with the Dyeing poison frog (Dendrobates tinctorius) to ask if this species can metabolize PTX 251D into aPTX 267A and what gene expression changes are associated with PTX 251D exposure in the intestines, liver, and skin. We found that D. tinctorius can metabolize PTX 251D into aPTX 267A, and that PTX 251D exposure changed the expression level of genes involved in immune system function and small molecule metabolism and transport. To better understand the functional significance of these changes in gene expression, we then conducted a series of high-throughput screens to determine the molecular targets of PTX 251D and identify potential proteins responsible for metabolism of PTX 251D into aPTX 267A. Although screens of PTX 251D binding human voltage-gated ion channels and G-protein coupled receptors were inconclusive, we identified human CYP2D6 as a rapid metabolizer of PTX 251D in a cytochrome P450 screen. Furthermore, a CYP2D6-like gene had increased expression in the intestines of animals fed PTX, suggesting this protein may be involved in PTX metabolism. These results show that individual alkaloids can modify gene expression across tissues, including genes involved in alkaloid metabolism. More broadly, this work suggests that specific alkaloid classes in wild diets may induce physiological changes for targeted accumulation and metabolism.

The impact of legacy and novel perfluoroalkyl substances on human cytochrome P450: An in vitro study on the inhibitory potential and underlying mechanisms

Per- and polyfluoroalkyl substances (PFASs) are a group of synthetic compounds with a wide range of industrial applications. PFOA and PFOS have been the most extensively studied and have been associated with hepatotoxicity. Recently, the interaction with cytochrome P450 (CYP) has been proposed as a potential key molecular event leading to PFAS-induced hepatotoxicity. In the present study, we aimed to determine a structure-activity relationship between thirteen PFASs and their inhibitory potential on the activities of four CYPs (CYP2E1, CYP2D6, CYP3A4 and CYP2C19). The influence of PFASs (5- 3200 µM) on CYP enzyme activities was measured using the Vivid® P450 metabolism assays. Using the same assays, Michaelis-Menten saturation curves were determined to explore the type of PFAS-induced CYP inhibition. Most PFASs were capable of inhibiting activity of the tested CYPs, as shown by their IC50 values. CYP2E1 is particularly inhibited by 3:1 FTOH, PFOA, and PFOS, whereas CYP2D6 is inhibited by PFHxS, PFHpA, PFOA, PFOS, PFNA, and PFDA. Additionally, CYP3A4 is most strongly inhibited by PFHxS, PFOA, PFOS, PFNA, and PFDA. Finally, CYP2C19 is inhibited by PFBS, PFHxS, PFHpA, PFOA, PFOS, PFNA, and PFDA. Interestingly, PFHxA and PFHxS induced an increase in CYP2E1 activity, whereas 4:2 FTOH strongly induced CYP2D6 activity. The mechanism of inhibition of CYPs by PFASs differed per CYP isoenzyme. CYP3A4 was competitively inhibited by PFBS, PFHxS, PFOS, PFNA and PFDA and non-competitively by PFOA. Additionally, CYP2C19 was competitively inhibited by PFHxA, PFOS and PFNA, whereas PFBS and PFHxS induced a mixed inhibition. Inhibition of CYP2C19 by PFHpA was atypical with an increased Vmax and a decreased Km. Finally, PFHxS competitively inhibited CYP2D6, whereas PFBS, PFOA, PFOS, PFDA and PFNA induced an atypical inhibition. Our results show that CYP inhibition by PFASs appears to be structure-dependent as well as CYP dependent. Inhibition of CYP2D6, CYP2C19 and CYP3A4 increased with increasing chain-lengths between six and nine carbons. The PFTOHs were only able to inhibit CYP2E1 and did not affect any of the other CYPS. Some PFASs remarkably induced the enzyme activity of CYPs. These results indicate that in addition to PFOA and PFOS, multiple novel PFASs may alter drug metabolism by the interference with CYPs.

Interaction study of engeletin toward cytochrome P450 3A4 and 2D6 by multi-spectroscopy and molecular docking

The inhibitory effects of engeletin on the activities of human cytochrome P450 3A4 and 2D6 (CYP3A4 and CYP2D6) were investigated by enzyme kinetics, multi-spectroscopy and molecular docking. Engeletin was found to strongly inhibit CYP3A4 and CYP2D6, with the IC₅₀ of 1.32 μM and 2.87 μM, respectively. The inhibition modes of engeletin against CYP3A4 and CYP2D6 were a competitive type and a mixed type, respectively. The fluorescence of the two CYPs was quenched statically by engeletin, which was bound to CYP3A4 stronger than to CYP2D6 at the same temperature. Circular dichroism spectroscopy, three-dimensional fluorescence, ultraviolet-visible spectroscopy and synchronous fluorescence confirmed that the conformation and micro-environment of the two CYPs protein were changed after binding with engeletin. Molecular docking, ultraviolet-visible spectroscopy and the fluorescence data revealed that engeletin had strong binding affinity to the two CYPs through hydrogen and van der Waals forces. The findings here suggested that engeletin may cause the herb-drug interactions for its inhibition of CYP3A4 and CYP2D6 activities.

Role of P34S, G169R, R296C, and S486T Substitutions in Ligand Access and Catalysis for Cytochrome P450 2D6 Allelic Variants CYP2D6*14A and CYP2D6*14B

BACKGROUND

Genetic polymorphism of cytochrome P450 (CYP) contributes to variability in drug metabolism, clearance, and response. This study aimed to investigate the functional and molecular basis for altered ligand binding and catalysis in CYP2D6*14A and CYP2D6*14B, two unique alleles common in the Asian population.

METHODS

CYP proteins expressed in Escherichia coli were studied using the substrate 3-cyano-7- ethoxycoumarin (CEC) and inhibitor probes (quinidine, fluoxetine, paroxetine, terbinafine) in the enzyme assay. Computer modelling was additionally used to create three-dimensional structures of the CYP2D6*14 variants.

RESULTS

Kinetics data indicated significantly reduced intrinsic clearance in CYP2D6*14 variants, suggesting that P34S, G169R, R296C, and S486T substitutions worked cooperatively to alter the conformation of the active site that negatively impacted the deethylase activity of CYP2D6. For the inhibition studies, IC50 values decreased in quinidine, paroxetine, and terbinafine but increased in fluoxetine, suggesting a varied ligand-specific susceptibility to inhibition. Molecular docking further demonstrated the role of P34S and R296C in altering access channel dimensions, thereby affecting ligand access and binding and subsequently resulting in varied inhibition potencies.

CONCLUSION

In summary, the differential selectivity of CYP2D6*14 variants for the ligands (substrate and inhibitor) was governed by the alteration of the active site and access channel architecture induced by the natural mutations found in the alleles.

Unraveling the Structural Basis of Selective Inhibition of Human Cytochrome P450 3A5

The human cytochrome P450 (CYP) CYP3A4 and CYP3A5 enzymes metabolize more than one-half of marketed drugs. They share high structural and substrate similarity and are often studied together as CYP3A4/5. However, CYP3A5 preferentially metabolizes several clinically prescribed drugs, such as tacrolimus. Genetic polymorphism in CYP3A5 makes race-based dosing adjustment of tacrolimus necessary to minimize acute rejection after organ transplantation. Moreover, the differential tissue distribution and expression levels of CYP3A4 and CYP3A5 can aggravate toxicity during treatment. Therefore, selective inhibitors of CYP3A5 are needed to distinguish the role of CYP3A5 from that of CYP3A4 and serve as starting points for potential therapeutic development. To this end, we report the crystal structure of CYP3A5 in complex with a previously reported selective inhibitor, clobetasol propionate (CBZ). This is the first CYP3A5 structure with a type I inhibitor, which along with the previously reported substrate-free and type II inhibitor-bound structures, constitute the main CYP3A5 structural modalities. Supported by structure-guided mutagenesis analyses, the CYP3A5-CBZ structure showed that a unique conformation of the F-F' loop in CYP3A5 enables selective binding of CBZ to CYP3A5. Several polar interactions, including hydrogen bonds, stabilize the position of CBZ to interact with this unique F-F' loop conformation. In addition, functional and biophysical assays using CBZ analogs highlight the importance of heme-adjacent moieties for selective CYP3A5 inhibition. Our findings can be used to guide further development of more potent and selective CYP3A5 inhibitors.

The Central Role of Cytochrome P450 in Xenobiotic Metabolism-A Brief Review on a Fascinating Enzyme Family

Human Cytochrome P450 (CYP) enzymes constitute a superfamily of membrane-bound hemoproteins that are responsible for the metabolism of a wide variety of clinically, physiologically, and toxicologically important compounds. These heme-thiolate monooxygenases play a pivotal role in the detoxification of xenobiotics, participating in the metabolism of many structurally diverge compounds. This short-review is intended to provide a summary on the major roles of CYPs in Phase I xenobiotic metabolism. The manuscript is focused on eight main topics that include the most relevant aspects of past and current CYP research. Initially, (I) a general overview of the main aspects of absorption, distribution, metabolism, and excretion (ADME) of xenobiotics are presented. This is followed by (II) a background overview on major achievements in the past of the CYP research field. (III) Classification and nomenclature of CYPs is briefly reviewed, followed by (IV) a summary description on CYP's location and function in mammals. Subsequently, (V) the physiological relevance of CYP as the cornerstone of Phase I xenobiotic metabolism is highlighted, followed by (VI) reviewing both genetic determinants and (VI) nongenetic factors in CYP function and activity. The last topic of the review (VIII) is focused on the current challenges of the CYP research field.

Uncovering of cytochrome P450 anatomy by SecStrAnnotator

Protein structural families are groups of homologous proteins defined by the organization of secondary structure elements (SSEs). Nowadays, many families contain vast numbers of structures, and the SSEs can help to orient within them. Communities around specific protein families have even developed specialized SSE annotations, always assigning the same name to the equivalent SSEs in homologous proteins. A detailed analysis of the groups of equivalent SSEs provides an overview of the studied family and enriches the analysis of any particular protein at hand. We developed a workflow for the analysis of the secondary structure anatomy of a protein family. We applied this analysis to the model family of cytochromes P450 (CYPs)-a family of important biotransformation enzymes with a community-wide used SSE annotation. We report the occurrence, typical length and amino acid sequence for the equivalent SSE groups, the conservation/variability of these properties and relationship to the substrate recognition sites. We also suggest a generic residue numbering scheme for the CYP family. Comparing the bacterial and eukaryotic part of the family highlights the significant differences and reveals a well-known anomalous group of bacterial CYPs with some typically eukaryotic features. Our workflow for SSE annotation for CYP and other families can be freely used at address https://sestra.ncbr.muni.cz .

Descriptors of Cytochrome Inhibitors and Useful Machine Learning Based Methods for the Design of Safer Drugs

Roughly 2.8% of annual hospitalizations are a result of adverse drug interactions in the United States, representing more than 245,000 hospitalizations. Drug-drug interactions commonly arise from major cytochrome P450 (CYP) inhibition. Various approaches are routinely employed in order to reduce the incidence of adverse interactions, such as altering drug dosing schemes and/or minimizing the number of drugs prescribed; however, often, a reduction in the number of medications cannot be achieved without impacting therapeutic outcomes. Nearly 80% of drugs fail in development due to pharmacokinetic issues, outlining the importance of examining cytochrome interactions during preclinical drug design. In this review, we examined the physiochemical and structural properties of small molecule inhibitors of CYPs 3A4, 2D6, 2C19, 2C9, and 1A2. Although CYP inhibitors tend to have distinct physiochemical properties and structural features, these descriptors alone are insufficient to predict major cytochrome inhibition probability and affinity. Machine learning based in silico approaches may be employed as a more robust and accurate way of predicting CYP inhibition. These various approaches are highlighted in the review.

In Silico Structural, Functional, and Phylogenetic Analysis of Cytochrome (CYPD) Protein Family

Cytochrome (CYP) enzymes catalyze the metabolic reactions of endogenous and exogenous compounds. The superfamily of enzymes is found across many organisms, regardless of type, except for plants. Information was gathered about CYP2D enzymes through protein sequences of humans and other organisms. The secondary structure was predicted using the SOPMA. The structural and functional study of human CYP2D was conducted using ProtParam, SOPMA, Predotar 1.03, SignalP, TMHMM 2.0, and ExPASy. Most animals shared five central motifs according to motif analysis results. The tertiary structure of human CYP2D, as well as other animal species, was predicted by Phyre2. Human CYP2D proteins are heavily conserved across organisms, according to the findings. This indicates that they are descended from a single ancestor. They calculate the ratio of alpha-helices to extended strands to beta sheets to random coils. Most of the enzymes are alpha-helix, but small amounts of the random coil were also found. The data were obtained to provide us with a better understanding of mammalian proteins' functions and evolutionary relationships.

In vitro inhibitory effects of glucosamine, chondroitin and diacerein on human hepatic CYP2D6

OBJECTIVES

Glucosamine, chondroitin and diacerein are natural compounds commonly used in treating osteoarthritis. Their concomitant intake may trigger drug-natural product interactions. Cytochrome P450 (CYP) has been implicated in such interactions. Cytochrome P450 2D6 (CYP2D6) is a major hepatic CYP involved in metabolism of 25% of the clinical drugs. This study aimed to investigate the inhibitory effect of these antiarthritic compounds on CYP2D6.

METHODS

CYP2D6 was heterologously expressed in Escherichia coli. CYP2D6-antiarthritic compound interactions were studied using in vitro enzyme kinetics assay and molecular docking.

RESULTS

The high-performance liquid chromatography (HPLC)-based dextromethorphan O-demethylase assay was established as CYP2D6 marker. All glucosamines and chondroitins weakly inhibited CYP2D6 (IC₅₀ values >300 µM). Diacerein exhibited moderate inhibition with IC₅₀ and K _i values of 34.99 and 38.27 µM, respectively. Its major metabolite, rhein displayed stronger inhibition potencies (IC₅₀=26.22 μM and K _i =32.27 μM). Both compounds exhibited mixed-mode of inhibition. In silico molecular dockings further supported data from the in vitro study. From in vitro-in vivo extrapolation, rhein presented an area under the plasma concentration-time curve (AUC) ratio of 1.5, indicating low potential to cause in vivo inhibition.

CONCLUSIONS

Glucosamine, chondroitin and diacerein unlikely cause clinical interaction with the drug substrates of CYP2D6. Rhein, exhibits only low potential to cause in vivo inhibition.

In vitro inhibitory effects of glucosamine, chondroitin and diacerein on human hepatic CYP2D6

OBJECTIVES

Glucosamine, chondroitin and diacerein are natural compounds commonly used in treating osteoarthritis. Their concomitant intake may trigger drug-natural product interactions. Cytochrome P450 (CYP) has been implicated in such interactions. Cytochrome P450 2D6 (CYP2D6) is a major hepatic CYP involved in metabolism of 25% of the clinical drugs. This study aimed to investigate the inhibitory effect of these antiarthritic compounds on CYP2D6.

METHODS

CYP2D6 was heterologously expressed in Escherichia coli. CYP2D6-antiarthritic compound interactions were studied using in vitro enzyme kinetics assay and molecular docking.

RESULTS

The high-performance liquid chromatography (HPLC)-based dextromethorphan O-demethylase assay was established as CYP2D6 marker. All glucosamines and chondroitins weakly inhibited CYP2D6 (IC50 values >300 µM). Diacerein exhibited moderate inhibition with IC50 and K i values of 34.99 and 38.27 µM, respectively. Its major metabolite, rhein displayed stronger inhibition potencies (IC50=26.22 μM and K i =32.27 μM). Both compounds exhibited mixed-mode of inhibition. In silico molecular dockings further supported data from the in vitro study. From in vitro-in vivo extrapolation, rhein presented an area under the plasma concentration-time curve (AUC) ratio of 1.5, indicating low potential to cause in vivo inhibition.

CONCLUSIONS

Glucosamine, chondroitin and diacerein unlikely cause clinical interaction with the drug substrates of CYP2D6. Rhein, exhibits only low potential to cause in vivo inhibition.

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.

P450 variations bifurcate the early terpene indole alkaloid pathway in Catharanthus roseus and Camptotheca acuminata

Divergent terpene indole alkaloid (TIA) pathways in Catharanthus roseus and Camptotheca acuminata generate vinblastine and vincristine, and camptothecin, respectively. In contrast to Catharanthus which feeds secologanin (from methylated loganin) into its species-specific late pathway, Camptotheca feeds secologanic acid (from unmethylated loganic acid) into its late pathway. Having identified putative Camptotheca secologanic acid synthases (SLASs) and cytochrome P450 reductases (CPRs) in transcriptome databases, we have demonstrated that two P450s, CYP72A564 and CYP72A565, are capable of utilizing both loganic acid and loganin to generate secologanic acid and secologanin. We have extended the previous report of these activities by CYP72A565 and CYP72A610 (Yang et al., 2019) by demonstrating that both Arabidopsis CPRs (ATR1, ATR2) couple with these CYP72A proteins in yeast microsomal assays and that purified Camptotheca CPR1 couples with them in in vitro reconstitution assays. Kinetic analyses of purified full-length Camptotheca SLASs have indicated that both process loganic acid with nearly identical catalytic rates and efficiencies as measured by their k_cat and k_cat/K_M. In contrast, CYP72A564 processes loganin with two-fold greater efficiency than CYP72A565 correlating with the former's 3-fold greater affinity for loganin. The closely-related CYP72A730 does not bind or process either compound. Molecular modeling of these three proteins and comparisons with Catharanthus secologanin synthase (SLS) have identified key differences that likely determine their SLAS versus SLS selectivities. Our ability to reconstitute these SLAS/SLS activities provides valuable tools for further examinations of the residues involved in substrate recognition and determinations of their unusual mechanism of C-C bond scission.

Paroxetine-Overview of the Molecular Mechanisms of Action

In the 21st century and especially during a pandemic, the diagnosis and treatment of depression is an essential part of the daily practice of many family doctors. It mainly affects patients in the age category 15–44 years, regardless of gender. Anxiety disorders are often diagnosed in children and adolescents. Social phobias can account for up to 13% of these diagnoses. Social anxiety manifests itself in fear of negative social assessment and humiliation, which disrupts the quality of social functioning. Treatment of the above-mentioned disorders is based on psychotherapy and pharmacotherapy. Serious side effects or mortality from antidepressant drug overdose are currently rare. Recent studies indicate that paroxetine (ATC code: N06AB), belonging to the selective serotonin reuptake inhibitors, has promising therapeutic effects and is used off-label in children and adolescents. The purpose of this review is to describe the interaction of paroxetine with several molecular targets in various points of view including the basic chemical and pharmaceutical properties. The central point of the review is focused on the pharmacodynamic analysis based on the molecular mechanism of binding paroxetine to various therapeutic targets.

Molecular level insights into the direct health impacts of some organic aerosol components

Quantum chemistry and biomodeling indicate that the studied organic aerosol components cannot directly cause oxidative stress or mutagenicity/carcinogenicity.

Crystal Structures of Drug-Metabolizing CYPs

The complex enzyme kinetics displayed by drug-metabolizing cytochrome P450 enzymes (CYPs) (see Chapter 9 ) can, in part, be explained by an examination of their crystallographic protein structures. Fortunately, despite low sequence similarity between different families of drug-metabolizing CYPs, there exists a high degree of structural homology within the superfamily. This similarity in the protein fold allows for a direct comparison of the structural features of CYPs that contribute toward differences in substrate binding, heterotropic and homotropic cooperativity, and genetic variability in drug metabolism. In this chapter, we first provide an overview of the nomenclature and the role of structural features that are common in all CYPs. We then apply these definitions to understand the different substrate specificities and functions in the CYP3A, CYP2C, and CYP2D families of enzymes.

Molecular determinant of substrate binding and specificity of cytochrome P450 2J2

Cytochrome P450 2J2 (CYP2J2) is responsible for the epoxidation of endogenous arachidonic acid, and is involved in the metabolism of exogenous drugs. To date, no crystal structure of CYP2J2 is available, and the proposed structural basis for the substrate recognition and specificity in CYP2J2 varies with the structural models developed using different computational protocols. In this study, we developed a new structural model of CYP2J2, and explored its sensitivity to substrate binding by molecular dynamics simulations of the interactions with chemically similar fluorescent probes. Our results showed that the induced-fit binding of these probes led to the preferred active poses ready for the catalysis by CYP2J2. Divergent conformational dynamics of CYP2J2 due to the binding of each probe were observed. However, a stable hydrophobic clamp composed of residues I127, F310, A311, V380, and I487 was identified to restrict any substrate access to the active site of CYP2J2. Molecular docking of a series of compounds including amiodarone, astemizole, danazol, ebastine, ketoconazole, terfenadine, terfenadone, and arachidonic acid to CYP2J2 confirmed the role of those residues in determining substrate binding and specificity of CYP2J2. In addition to the flexibility of CYP2J2, the present work also identified other factors such as electrostatic potential in the vicinity of the active site, and substrate strain energy and property that have implications for the interpretation of CYP2J2 metabolism.

Galeon: A Biologically Active Molecule with In Silico Metabolite Prediction, In Vitro Metabolic Profiling in Rat Liver Microsomes, and In Silico Binding Mechanisms with CYP450 Isoforms

Galeon, a natural cyclic-diarylheptanoid (CDH), which was first isolated from Myrica gale L., is known to have potent cytotoxicity against A549 cell lines, anti-tubercular activity against Mycobacterium tuberculosis H37Rv, chemo-preventive potential, and moderate topoisomerase inhibitory activity. Here, in silico metabolism and toxicity prediction of galeon by CYP450, in vitro metabolic profiling study in rat liver microsomes (RLMs), and molecular interactions of galeon-CYP450 isoforms were performed. An in silico metabolic prediction study showed demethyl and mono-hydroxy galeon were the metabolites with the highest predictability. Among the predicted metabolites, mono-hydroxy galeon was found to have plausible toxicities such as skin sensitization, thyroid toxicity, chromosome damage, and carcinogenicity. An in vitro metabolism study of galeon, incubated in RLMs, revealed eighteen Phase-I metabolites, nine methoxylamine, and three glutathione conjugates. Identification of possible metabolites and confirmation of their structures were carried out using ion-trap tandem mass spectrometry. In silico docking analysis of galeon demonstrated significant interactions with active site residues of almost all CYP450 isoforms.

Inhibitory effects of astilbin, neoastilbin and isoastilbin on human cytochrome CYP3A4 and 2D6 activities

Astilbin, neoastilbin and isoastilbin are three flavonoid isomers from Smilacis glabrae Roxb. (S. glabrae). Several studies have shown that consumption of flavonoids can increase the risk of food/drug-drug interaction by affecting the activities of human cytochrome CYP3A4 and 2D6. In the present study, an ultrahigh-performance liquid chromatography and triple quadrupole mass spectrometry method was developed for the determination of the interaction between three flavonoid isomers and two CYPs. Under the optimized reaction conditions, the K_m values were 18.9 and 36.4 μM and the V_max values were 0.02 and 0.20 μM/min for CYP3A4 and 2D6 in vitro, respectively. Astilbin showed the strongest inhibition on CYP3A4, followed by isoastilbin and neoastilbin with IC₅₀ values of 2.63, 3.03 and 6.51 μM. Neoastilbin showed the strongest inhibition on CYP2D6, followed by isoastilbin and astilbin, with IC₅₀ values of 1.48, 11.87 and 14.16 μM, respectively. The three isomers showed reversible inhibition on both enzymes. Neoastilbin and astilbin were noncompetitive type for CYP3A4 and 2D6, isoastilbin was a mixture and noncompetitive type for CYP3A4 and 2D6, respectively. Our study suggests that the three isomers may increase the risk of food/drug-drug interactions by affecting the activities of CYP3A4 and 2D6.

CYP2J2 Molecular Recognition: A New Axis for Therapeutic Design

Cytochrome P450 (CYP) epoxygenases are a special subset of heme-containing CYP enzymes capable of performing the epoxidation of polyunsaturated fatty acids (PUFA) and the metabolism of xenobiotics. This dual functionality positions epoxygenases along a metabolic crossroad. Therefore, structure-function studies are critical for understanding their role in bioactive oxy-lipid synthesis, drug-PUFA interactions, and for designing therapeutics that directly target the epoxygenases. To better exploit CYP epoxygenases as therapeutic targets, there is a need for improved understanding of epoxygenase structure-function. Of the characterized epoxygenases, human CYP2J2 stands out as a potential target because of its role in cardiovascular physiology. In this review, the early research on the discovery and activity of epoxygenases is contextualized to more recent advances in CYP epoxygenase enzymology with respect to PUFA and drug metabolism. Additionally, this review employs CYP2J2 epoxygenase as a model system to highlight both the seminal works and recent advances in epoxygenase enzymology. Herein we cover CYP2J2's interactions with PUFAs and xenobiotics, its tissue-specific physiological roles in diseased states, and its structural features that enable epoxygenase function. Additionally, the enumeration of research on CYP2J2 identifies the future needs for the molecular characterization of CYP2J2 to enable a new axis of therapeutic design.

In Silico Pharmacogenetics CYP2D6 Study Focused on the Pharmacovigilance of Herbal Antidepressants

The annual increase in depression worldwide together with an upward trend in the use of alternative medicine as treatment asks for developing reliable safety profiles of herbal based medicine. A considerable risk on adverse reactions exists when herbal remedies are combined with prescription medication. Around 25% of the drugs, including many antidepressants, depend on the activity of CYP2D6 for their metabolism and corresponding efficacy. Therefore, probing CYP2D6 inhibition by the active substances in herbal based medicine within the wild-type enzyme and clinically relevant allelic variants is crucial to avoid toxicity issues. In this in silico study several compounds with herbal origin suggested to have antidepressant activity were analyzed on their CYP2D6 wild-type and CYP2D6*53 inhibition potential using molecular docking. In addition, several pharmacokinetic properties were evaluated to assess their probability to cross the blood brain barrier and subsequently reach sufficient brain bioavailability for the modulation of central nervous system targets as well as characteristics which may hint toward potential safety issues.

Conversion of amino acids to aryl/heteroaryl ethanol metabolites using human CYP2D6-expressing live baker's yeast

Different natural aromatic/heterocyclic l-amino acids were biotransformed into aryl/heteroaryl ethanol metabolites via oxidative deamination, decarboxylation and reduction cascades using live baker's yeast cells producing intracellular human CYP2D6 enzyme. Among the three yeast strains expressing 3 different CYP2D6 variants, CYP2D6(2) (i.e. CYP2D6 wild-type) provided the best result under neutral pH conditions at RT. We have successfully converted six natural amino acids into their corresponding alcohols, having one carbon atom less, with moderate yields. Some of the downstream products like tryptophol and tyrosol induced the pTrKB (Tropomyosin receptor kinase B) activation pattern similar to that of BDNF (brain-derived neurotrophic factor), thereby depicting potential antidepressant activity. Control experiments and molecular modelling studies revealed that this tandem bio-transformation probably happens via a pyruvate intermediate. This study establishes that CYP2D6-expressing live yeast cells can be a powerful tool for the enzymatic C-N, C-C bond cleavage of amino-acids.

Computational prediction of cytochrome P450 inhibition and induction

Cytochrome P450 (CYP) enzymes play an important role in the phase I metabolism of many xenobiotics. Most drug-drug interactions (DDIs) associated with CYP are caused by either CYP inhibition or induction. The early detection of potential DDIs is highly desirable in the pharmaceutical industry because DDIs can cause serious adverse events, which can lead to poor patient health and drug development failures. Recently, many computational studies predicting CYP inhibition and induction have been reported. The current computational modeling approaches for CYP metabolism are classified as ligand- and structure-based; various techniques, such as quantitative structure-activity relationships, machine learning, docking, and molecular dynamic simulation, are involved in both the approaches. Recently, combining these two approaches have resulted in improvements in the prediction accuracy of DDIs. In this review, we present important, recent developments in the computational prediction of the inhibition of four clinically crucial CYP isoforms (CYP1A2, 2C9, 2D6, and 3A4) and three nuclear receptors (aryl hydrocarbon receptor, constitutive androstane receptor, and pregnane X receptor) involved in the induction of CYP1A2, 2B6, and 3A4, respectively.

The marmoset cytochrome P450 superfamily: Sequence/phylogenetic analyses, genomic structure, and catalytic function

The common marmoset (Callithrix jacchus) is a New World monkey that has attracted much attention as a potentially useful primate model for preclinical testing. A total of 36 marmoset cytochrome P450 (P450) isoforms in the P450 1-51 subfamilies have been identified and characterized by the application of genome analysis and molecular functional characterization. In this mini-review, we provide an overview of the genomic structures, sequence identities, and substrate selectivities of marmoset P450s compared with those of human P450s. Based on the sequence identity, phylogeny, and genomic organization of marmoset P450s, orthologous relationships were established between human and marmoset P450s. Twenty-four members of the marmoset P450 1A, 2A, 2B, 2C, 2D, 2E, 3A, 4A, and 4F subfamilies shared high degrees of homology in terms of cDNA (>89%) and amino acid sequences (>85%) with the corresponding human P450s; P450 2C76 was among the exceptions. Phylogenetic analysis using amino acid sequences revealed that marmoset P450s in the P450 1-51 families were located in the same clades as their human and macaque P450 homologs. This finding underlines the evolutionary closeness of marmoset P450s to their human and macaque homologs. Most marmoset P450 1-4 enzymes catalyzed the typical drug-metabolizing reactions of the corresponding human P450 homologs, except for some differences of P450 2A6 and 2B6. Consequently, it appears that the substrate specificities of enzymes in the P450 1-4 families are generally similar in marmosets and humans. The information presented here supports a better understanding of the functional characteristics of marmoset P450s and their similarities and differences with human P450s. It is hoped that this mini-review will facilitate the successful use of marmosets as primate models in drug metabolism and pharmacokinetic studies.

Effects of 24 CYP2D6 variants found in Chinese population on the metabolism of clonidine in vitro

BACKGROUND AND OBJECTIVES

Clonidine has been clinically used to treat Tourette's syndrome for decades. There was research finding that clonidine possessed the best risk-benefit ratio, especially for patients associated with attention deficit hyperactivity disorder. CYP2D6 is a significant member of Cytochrome P450 enzymes. The genetic polymorphisms of CYP2D6 greatly affect the clinical effects of drugs even lead to side effects and medical malpractice. Our goal is to research the effect of CYP2D6 genetic polymorphism on the metabolism of clonidine and evaluate the functions of 22 CYP2D6 allelic variants in vitro, which were discovered in Chinese Han population recently.

METHODS

This study was carried out through a mature incubation system. The wild-type CYP2D6*1 and 24 variants (CYP2D6*2, CYP2D6*10 and 22 novel CYP2D6 variants) were expressed in insect cells, and the catalytic activity of all the variants were assessed by substrate clonidine. Metabolite 4-OH clonidine was accurately detected via ultra-performance liquid-chromatography tandem mass spectrometry to evaluate the effect of CYP2D6 genetic polymorphism on the clonidine.

RESULT

Among the 22 novel CYP2D6 variants, the intrinsic clearance (Vmax/Km) of 21 variants were significantly decreased (from 1.53% to 83.25%) compared to the wild-type. In particular, the following seven variants (CYP2D6* 2, CYP2D6* 10, CYP2D6* 93, CYP2D6* 95, E215K, V327 M and R497C) attract more attention, of which the intrinsic clearance decreased more than 70% compared to the wild-type. Because the variants with significantly reduced intrinsic clearance are more likely to cause adverse reactions than the variants with increased or little changed intrinsic clearance. In addition, the related pharmacokinetic parameters of CYP2D6*92 and CYP2D6*96 could not be acquired for the defect of CYP2D6 nucleotide.

CONCLUSION

We comprehensively evaluated the effect of 22 novel CYP2D6 variants on the metabolism of clonidine for the first time and hoped corresponding data provide a reference for metabolism of clonidine for further studies in vivo, and extend our understanding of the clinical drug toxicity or ineffectiveness by CYP2D6 genetic polymorphism.