On the human CYP2C9*13 variant activity reduction: a molecular dynamics simulation and docking study

The relationship between CYP2C9 gene polymorphisms and azilsartan metabolism in vitro

BACKGROUND

The gene polymorphisms of the CYP2C9, as well as the substrate specificity of the enzyme, result in different clearances for different substrates by CYP2C9 variants.

RESEARCH DESIGNAND METHODS

The CYP2C9 wild type and 38 CYP2C9 variants, expressed in insectmicrosomes, were incubated with azilsartan. The resulting metabolite,O-desethyl azilsartan, was determined by HPLC-MS/MS. The enzyme kineticparameters of the 38 variants were calculated and compared with the wild type.Subsequently, we selected CYP2C9*1, *2, and  *3 as target proteins for molecular docking with azilsartan to elucidate the mechanisms underlying changes in enzyme function.

RESULTS

Compared with CYP2C9*1, three variants (CYP2C9*29, *39, and *49) exhibited markedlyincreased CLint values (from 170%-275%, *p < 0.05), whereas 28 variants exhibited significantly decreased CLint values (from 3-63%,*p < 0.05). The molecular docking results showed that the binding energy of CYP2C9*2 and *3 was lower than that of the wild type.

CONCLUSION

Thisassessment revealed the effect of CYP2C9 gene polymorphisms on azilsartan metabolism, establishing a theoretical basis for further in-vivo studies and clinical applications. This study will help expand the database of CYP2C9 gene-drug pairs and identify appropriate treatment strategies for azilsartan, contributing to the field of precision medicine.

APF2: an improved ensemble method for pharmacogenomic variant effect prediction

Lack of efficacy or adverse drug response are common phenomena in pharmacological therapy causing considerable morbidity and mortality. It is estimated that 20-30% of this variability in drug response stems from variations in genes encoding drug targets or factors involved in drug disposition. Leveraging such pharmacogenomic information for the preemptive identification of patients who would benefit from dose adjustments or alternative medications thus constitutes an important frontier of precision medicine. Computational methods can be used to predict the functional effects of variant of unknown significance. However, their performance on pharmacogenomic variant data has been lackluster. To overcome this limitation, we previously developed an ensemble classifier, termed APF, specifically designed for pharmacogenomic variant prediction. Here, we aimed to further improve predictions by leveraging recent key advances in the prediction of protein folding based on deep neural networks. Benchmarking of 28 variant effect predictors on 530 pharmacogenetic missense variants revealed that structural predictions using AlphaMissense were most specific, whereas APF exhibited the most balanced performance. We then developed a new tool, APF2, by optimizing algorithm parametrization of the top performing algorithms for pharmacogenomic variations and aggregating their predictions into a unified ensemble score. Importantly, APF2 provides quantitative variant effect estimates that correlate well with experimental results (R2 = 0.91, p = 0.003) and predicts the functional impact of pharmacogenomic variants with higher accuracy than previous methods, particularly for clinically relevant variations with actionable pharmacogenomic guidelines. We furthermore demonstrate better performance (92% accuracy) on an independent test set of 146 variants across 61 pharmacogenes not used for model training or validation. Application of APF2 to population-scale sequencing data from over 800,000 individuals revealed drastic ethnogeographic differences with important implications for pharmacotherapy. We thus think that APF2 holds the potential to improve the translation of genetic information into pharmacogenetic recommendations, thereby facilitating the use of Next-Generation Sequencing data for stratified medicine.

Enzymatic activity of 38 CYP2C9 genotypes on ibuprofen

BACKGROUND AND OBJECTIVE

Ibuprofen, a common non-steroidal anti-inflammatory drug, is used clinically for pain relief and antipyretic treatment worldwide. However, regular or long-term use of ibuprofen may lead to a series of adverse reactions, including gastrointestinal bleeding, hypertension and kidney injury. Previous studies have shown that CYP2C9 gene polymorphism plays an important role in the elimination of various drugs, which leads to the variation in drug efficacy. This study aimed to evaluate the effect of 38 CYP2C9 genotypes on ibuprofen metabolism.

METHODS

Thirty-eight recombinant human CYP2C9 microsomal enzymes were obtained using a frugiperda 21 insect expression system according to a previously described method. Assessment of the catalytic function of these variants was completed via a mature incubation system: 5 pmol CYP2C9*1 and 38 CYP2C9 variants recombinant human microsomes, 5 μL cytochrome B5, ibuprofen (5-1000 μM), and Tris-HCl buffer (pH 7.4). The ibuprofen metabolite contents were determined using HPLC analysis. HPLC analysis included a UV detector, Plus-C18 column, and mobile phase [50% acetonitrile and 50% water (containing 0.05% trifluoroacetic acid)]. The kinetic parameters of the CYP2C9 genotypes were obtained by Michaelis-Menten curve fitting.

RESULTS

The intrinsic clearance (CL_int) of eight variants was not significantly different from CYP2C9*1; four CYP2C9 variants (CYP2C9*38, *44, *53 and *59) showed significantly higher CL_int (increase by 35%-230%) than that of the wild-type; the remaining twenty-six variants exhibited significantly reduced CL_int (reduced by 30%-99%) compared to that of the wild-type.

CONCLUSION

This is the first systematic evaluation of the catalytic characteristics of 38 CYP2C9 genotypes involved ibuprofen metabolism. Our results provide a corresponding supplement to studies on CYP2C9 gene polymorphisms and kinetic characteristics of different variants. We need to focus on poor metabolizers (PMs) with severely abnormal metabolic functions, because they are more susceptible to drug exposure.

Enzymatic activity on valsartan of 38 CYP2C9 variants from the Chinese population

BACKGROUND AND OBJECTIVE

Valsartan is widely used for the treatment of moderate hypertension. However, previous studies have found that efficacy of the valsartan depends on the dose and intake. Cytochrome P450 (CYP) 2C9 metabolizes ∼15% of the clinical drugs. Genetic polymorphisms of CYP2C9 markedly affect the safety and effectiveness of many drugs, which might lead to adverse reactions and therapeutic failure. Twenty-four novel CYP2C9 variants (*36-*60) had been previously discovered via gene sequencing in the Han population. Our study aims to evaluate the impact of 38 CYP2C9 variants from the Chinese population on valsartan metabolism compared with CYP2C9*1 in vitro.

METHODS

Wild-type CYP2C9*1 and other CYP2C9 variants were expressed in Spodoptera frugiperda 21 insect cells. Incubations were performed at 37 °C with 20-2000 μM substrate for 30 min. The metabolite 4-OH valsartan was determined via UPLC-MS/MS.

RESULTS

Among the 38 CYP2C9 variants, the enzymatic activities of most variants were significantly altered compared with the wild-type. Three variants (CYP2C9*27, *40 and *49) exhibited increased intrinsic clearance values (134-153% relative clearance). However, 12 variants (CYP *8, *13, *16, *19, *33, *36, *42, *43, *45, *52, *54, *58) caused >90% decreases in the relative clearance of valsartan compared to CYP2C9*1.

CONCLUSIONS

Our research provides systematic data for evaluating the effects of CYP2C9 variants on valsartan metabolism in the Chinese population. These results will expand our understanding of the impact of CYP2C9 genetic polymorphisms on valsartan metabolism and will contribute to precision medicine.

miR-26a desensitizes non-small cell lung cancer cells to tyrosine kinase inhibitors by targeting PTPN13

Epidermal growth factor receptor (EGFR)-targeted tyrosine kinase inhibitors (TKIs) have emerged as first-line drugs for non-small cell lung cancers (NSCLCs). However, the resistance to TKIs represents the key limitation for their therapeutic efficacy. We found that miR-26a was upregulated in gefitinib-refractory NSCLCs; miR-26a is downstream of EGFR signaling and directly targets and silences protein tyrosine phosphatase non-receptor type 13 (PTPN13) to maintain the activation of Src, a dephosphorylation substrate of PTPN13, thus reinforcing EGFR pathway in a regulatory circuit. miR-26a inhibition significantly improved NSCLC responses to gefitinib. These data revealed a novel mechanism of NSCLC resistance to TKI treatment.

TOPK promotes lung cancer resistance to EGFR tyrosine kinase inhibitors by phosphorylating and activating c-Jun

Tyrosine kinase inhibitors (TKIs) targeting the epidermal growth factor receptor (EGFR) have shown promising clinical efficacy in non-squamous non-small cell lung cancer (NSCLC); however, resistance is frequently observed in malignant cells, operating through a mechanism that remains largely unknown. The present study shows that T-lymphokine-activated killer cell-originated protein kinase (TOPK) is upregulated in NSCLC and excessively activated in TKI-refractory cells. TOPK dictates the responsiveness of lung cancers to the EGFR-targeted TKI gefitinib through the transcription factor AP-1 component c-Jun. TOPK binds directly to and phosphorylates c-Jun, which consequently activates the transcription of AP-1 target genes, including CCND1 and CDC2. TOPK silencing sensitizes EGFR-TKI-resistant lung cancer cells to gefitinib and increases gefitinib efficacy in preclinical lung adenocarcinoma xenograft models. These findings represent a novel mechanism of lung cancer resistance to TKIs and suggest that TOPK may have value both as a predictive biomarker and as a therapeutic target: TOPK-targeted therapy may synergize with EGFR-targeted therapy in lung cancers.

Effects of CYP2C9 genetic polymorphisms on the pharmacokinetics of celecoxib and its carboxylic acid metabolite

Celecoxib, a selective cyclooxygenase (COX)-2 inhibitor, is used for the treatment of rheumatoid arthritis and osteoarthritis. The predominant hepatic metabolism of celecoxib to celecoxib carboxylic acid (CCA) is mediated mainly by CYP2C9. We investigated the effects of the major CYP2C9 genetic variants in Asian populations, CYP2C9*3 and CYP2C9*13, on the pharmacokinetics of celecoxib and its carboxylic acid metabolite in healthy Korean subjects. A single 200-mg oral dose of celecoxib was given to 52 Korean subjects with different CYP2C9 genotypes: CYP2C9EM (n = 26; CYP2C9*1/*1), CYP2C9IM (n = 24; CYP2C9*1/*3 and *1/*13), and CYP2C9PM (n = 2; CYP2C9*3/*3). Celecoxib and CCA concentrations in plasma samples collected up to 48 or 96 h after drug intake were determined by HPLC-MS/MS. The mean area under the plasma concentration-time curve (AUC_0-∞) of celecoxib was increased 1.63-fold (P < 0.001), and the apparent oral clearance (CL/F) of celecoxib was decreased by 39.6% in the CYP2C9IM genotype group compared with that of CYP2C9EM (P < 0.001). The overall pharmacokinetic parameters for celecoxib in CYP2C9*1/*13 subjects were similar to those in CYP2C9*1/*3 subjects. Two subjects with CYP2C9PM genotype both showed markedly higher AUC_0-∞, prolonged half-life, and lower CL/F for celecoxib than did subjects with CYP2C9EM and IM genotypes. CYP2C9*3 and CYP2C9*13 variant alleles significantly affected the plasma concentration of celecoxib.

Effects of CYP2C9 genetic polymorphisms on the pharmacokinetics of zafirlukast

Zafirlukast, a cysteinyl leukotriene receptor antagonist, is indicated for the treatment of patients with mild to moderate asthma. Zafirlukast is metabolized mainly by CYP3A4 and CYP2C9. We investigated the effects of the major CYP2C9 variant alleles in Asian populations, CYP2C9*3 and CYP2C9*13, on the pharmacokinetics of zafirlukast in healthy Korean subjects. A single 20-mg oral dose of zafirlukast was given to 23 Korean male subjects divided into two genotype groups according to CYP2C9 genotypes, CYP2C9EM (n = 11; CYP2C9*1/*1) and CYP2C9IM (n = 12; 9 and 3 carriers of CYP2C9*1/*3 and *1/*13, respectively). Zafirlukast concentrations were determined using a validated HPLC-MS/MS analytical method in plasma samples collected after the drug intake. Compared with the CYP2C9EM group, the Cmax and AUCinf of zafirlukast in the CYP2C9IM group were 1.44- and 1.70-fold higher, respectively (p < 0.01 and p < 0.0001). The CL/F of zafirlukast was 42.8 % lower in the CYP2C9IM group compared with the CYP2C9EM group (p < 0.001). Slightly higher Cmax and AUC, and lower CL/F of zafirlukast were observed in subjects with the CYP2C9*1/*13 genotype compared with the CYP2C9*1/*3 genotype subjects. CYP2C9*3 and CYP2C9*13 alleles significantly affected the plasma concentrations of zafirlukast.

Linobiflavonoid inhibits human lung adenocarcinoma A549 cells: effect on tubulin protein

The antitumor bioactivities of linobiflavonoid were studied through evaluating its in vitro cytotoxicity against several cell lines (A549, H1975, SMMC-7721, HEP-2 and Vero cells), with the aid of 3-(4,5)-dimethylthiazoly1)-3,5-diphenytetrazolium bromide (MTT) assay. It was found that linobiflavonoid shows more notable inhibiting activity against A549 cells, with IC50 value of 4.67 μM. Furthermore, western blot analysis revealed that linobiflavonoid is able to increase the expression of β-tubulin, whereas not α-tubulin. In virtuale simulations indicated that linobiflavonoid specifically interacts with the binding pocket which is located at the top of β-tubulin, due to the presence of strong hydrophobic effects between the core templates and the hydrophobic surface of the tubulin protein (TB) binding site. The binding energy (E inter ) was calculated to be -140.47 kcal/mol. Results above suggest that linobiflavonoid possesses anti-A549 properties relating to β-tubulin depolymerization inhibition.

Distal effect of amino acid substitutions in CYP2C9 polymorphic variants causes differences in interatomic interactions against (S)-warfarin

Cytochrome P450 2C9 (CYP2C9) is crucial in excretion of commonly prescribed drugs. However, changes in metabolic activity caused by CYP2C9 polymorphisms inevitably result in adverse drug effects. CYP2C9*2 and *3 are prevalent in Caucasian populations whereas CYP2C9*13 is remarkable in Asian populations. Single amino acid substitutions caused by these mutations are located outside catalytic cavity but affect kinetic activities of mutants compared to wild-type enzyme. To relate distal effects of these mutations and defective drug metabolisms, simulations of CYP2C9 binding to anti-coagulant (S)-warfarin were performed as a system model. Representative (S)-warfarin-bound forms of wild-type and mutants were sorted and assessed through knowledge-based scoring function. Interatomic interactions towards (S)-warfarin were predicted to be less favorable in mutant structures in correlation with larger distance between hydroxylation site of (S)-warfarin and reactive oxyferryl heme than wild-type structure. Using computational approach could delineate complication of CYP polymorphism in management of drug therapy.

EXPLORING THE MOLECULAR BASIS OF H5N1 HEMAGGLUTININ BINDING WITH CATECHINS IN GREEN TEA: A FLEXIBLE DOCKING AND MOLECULAR DYNAMICS STUDY

The influenza A (H5N1) virus attracts a worldwide attention and calls for the urgent development of novel antiviral drugs. In this study, explicitly solvated flexible docking and molecular dynamics (MD) simulations were used to study the interactions between the H5N1 sub-type hemagglutinin (HA) and various catechin compounds, including EC ([–]-epicatechin), EGC ([–]-epigallocatechin), ECG ([–]-epicatechin gallate) and EGCG ([–]-epigallocatechin gallate). The four compounds have respective binding specificities and their interaction energies with HA decrease in the order of EGCG (-133.52) > ECG (-111.11) > EGC (-97.94) > EC (-83.39). Units in kcal mol-1. Residues IleA267, LysA269, ArgB68 and GluB78 play important roles during all the binding processes. EGCG has the best bioactivity and shows potential as a lead compound. Besides, the importance was clarified for the functional groups it was revealed that the C5′ hydroxyl and trihydroxybenzoic acid groups are crucial for the catechin inhibitory activities, especially the latter. Combined with the structural and property analyses, this work also proposed the way to effectively modify the functional groups of EGCG. The experimental efforts are expected in order to actualize the catechin derivatives as novel anti-influenza agents in the near future.

MOLECULAR DOCKING STUDY OF THE AFFINITY OF CYP2C9 AND CYP2D6 FOR IMRECOXIB

With the aid of the automatic molecular docking, the affinity of CYP2C9 and CYP2D6 for imrecoxib was studied by InsightII/Affinity program. The results indicate that CYP2C9–imrecoxib complex has higher stability and stronger affinity because CYP2C9 has more favorable interaction energy (-62.72 kcal/mol) and higher Ludi score (610) with imrecoxib than CYP2D6 (-50.22 kcal/mol and 551) and this is consistent with the results of the kinetic experiments by Li et al. By analyzing the theoretical results combined with the experimental ones, we suggest that the affinity difference is caused by the difference of the structure between CYP2C9 and CYP2D6, and the most important residues for enzyme–substrate complexes, such as Phe476, Asn204, Phe100, Leu366 and Arg108 of CYP2C9 and Phe120, Glu216, and Phe483 of CYP2D6 were also identified.

Effects of CYP2C9*1/*13 on the pharmacokinetics and pharmacodynamics of meloxicam

AIMS

To determine the effects of the CYP2C9*1/*13 genotype on the pharmacokinetics and pharmacodynamics of meloxicam in Korean subjects.

METHODS

Meloxicam (15 mg) was orally administered to 21 healthy Korean volunteers with either the CYP2C9*1/*1 or the CYP2C9*1/*13 genotype. Plasma meloxicam concentrations were analysed by HPLC-UV for 72 h after drug administration. The pharmacodynamic effects of meloxicam were determined by measuring TXB(2) generated in blood.

RESULTS

The AUC(0,∞) and C(max) of meloxicam were 2.43- and 1.46-fold higher in the CYP2C9*1/*13 group than in the CYP2C9*1/*1 group, respectively. The oral clearance of meloxicam was significantly lower in the CYP2C9*1/*13 group (37.9% of wild type) than in the CYP2C9*1/*1 group. The t(1/2) of meloxicam was 1.84-fold longer in the CYP2C9*1/*13 group than in the CYP2C9*1/*1 group. The rate of TXB(2) production was significantly lower in the CYP2C9*1/*13 group than in the CYP2C9*1/*1 group.

CONCLUSIONS

The CYP2C9*1/*13 genotype is associated with decreased metabolism and increased pharmacodynamic effects of meloxicam.

Pathways and mechanisms for product release in the engineered haloalkane dehalogenases explored using classical and random acceleration molecular dynamics simulations

Eight mutants of the DhaA haloalkane dehalogenase carrying mutations at the residues lining two tunnels, previously observed by protein X-ray crystallography, were constructed and biochemically characterized. The mutants showed distinct catalytic efficiencies with the halogenated substrate 1,2,3-trichloropropane. Release pathways for the two dehalogenation products, 2,3-dichloropropane-1-ol and the chloride ion, and exchange pathways for water molecules, were studied using classical and random acceleration molecular dynamics simulations. Five different pathways, denoted p1, p2a, p2b, p2c, and p3, were identified. The individual pathways showed differing selectivity for the products: the chloride ion releases solely through p1, whereas the alcohol releases through all five pathways. Water molecules play a crucial role for release of both products by breakage of their hydrogen-bonding interactions with the active-site residues and shielding the charged chloride ion during its passage through a hydrophobic tunnel. Exchange of the chloride ions, the alcohol product, and the waters between the buried active site and the bulk solvent can be realized by three different mechanisms: (i) passage through a permanent tunnel, (ii) passage through a transient tunnel, and (iii) migration through a protein matrix. We demonstrate that the accessibility of the pathways and the mechanisms of ligand exchange were modified by mutations. Insertion of bulky aromatic residues in the tunnel corresponding to pathway p1 leads to reduced accessibility to the ligands and a change in mechanism of opening from permanent to transient. We propose that engineering the accessibility of tunnels and the mechanisms of ligand exchange is a powerful strategy for modification of the functional properties of enzymes with buried active sites.

Mechanism of CYP2C9 inhibition by flavones and flavonols

This article describes an in vitro investigation of the inhibition of cytochrome P450 (P450) 2C9 by a series of flavonoids made up of flavones (flavone, 6-hydroxyflavone, 7-hydroxyflavone, chrysin, baicalein, apigenin, luteolin, scutellarein, and wogonin) and flavonols (galangin, fisetin, kaempferol, morin, and quercetin). With the exception of flavone, all flavonoids were shown to inhibit CYP2C9-mediated diclofenac 4'-hydroxylation in the CYP2C9 RECO system, with K(i) value <or= 2.2 microM. In terms of the mechanism of inhibition, 6-hydroxyflavone was found to be a noncompetitive inhibitor of CYP2C9, whereas the other flavonoids were competitive inhibitors. Computer docking simulation and constructed mutants substituted at residue 100 of CYP2C9.1 indicate that the noncompetitive binding site of 6-hydroxyflavone lies beside Phe100, similar to the reported allosteric binding site of warfarin. The other flavonoids exert competitive inhibition through interaction with the substrate binding site of CYP2C9 accessed by flurbiprofen. These results suggest flavonoids can participate in interactions with drugs that act as substrates for CYP2C9 and provide a possible molecular basis for understanding cooperativity in human P450-mediated drug-drug interactions.