Reduced catalytic activity of human CYP2C9 natural alleles for gliclazide: Molecular dynamics simulation and docking studies

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.

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.

The determination of the effect(s) of solute carrier family 22-member 2 (SLC22A2) haplotype variants on drug binding via molecular dynamic simulation systems

Single nucleotide polymorphisms detected in the solute carrier member family-22 has been shown to result in a variable response in the treatment of type 2 diabetes mellitus with Metformin. This study predicted a three-dimensional protein structure for the SLC22A2 protein sequence using AlphaFold 2 and modelled five haplotypes within SLC22A2 protein structure observed in the Xhosa population of South Africa. The protein models were used to determine the effect(s) of haplotype variations on the transport function of Metformin and 10 other drugs by the SLC22A2 protein. Molecular dynamic simulation studies, molecular docking and interaction analysis of the five SLC22A2 haplotypes were performed in complex with the ligand 5RE in a POPC lipid bilayer to understand the mechanism of drug binding. Weakest binding free energy was found between 5RE and haplotype 1. Molecular docking studies indicated the top binding ligands as well as Metformin to bind inside the transport channel in all haplotypes increasing the probability of Metformin inhibition during co-administration of drugs. Metformin showed reduced binding affinity and number of interactions compared to the top four binding molecules. Molecular dynamic simulation analysis indicated that haplotypes 1, 3 and 4 were less stable than 2 and 5. The findings suggest haplotypes 4 and 5 having stronger preference for large inhibitor molecule binding in the active site and this could result in haplotypes 4 and 5 demonstrating reduced Metformin clearance via the SLC22A2 transporter during co-administration of drugs. The current study is the first to investigate the potential effect(s) of haplotype variation on the protein structure of SLC22A2 to assess its ability to transport Metformin in an indigenous South African population.

PAMPA model of gliclazide permeability: The impact of probiotic bacteria and bile acids

Gut microbiota and bile acids possess the ability to modify absorption and pharmacokinetic profile of numerous drugs. Since the variability of gliclazide response in patients cannot be explained only by genetic factors, the influence of gut microbiota and bile acids should be considered. The aim of this study was to determine the effects of probiotic bacteria and bile acids on the gliclazide permeability. The permeability of gliclazide with and without probiotic bacteria and bile acids (cholic acid, CA and deoxycholic acid, DCA) was tested using in vitro PAMPA model, at three different pH values (5.8, 6.5 and 7.4). Concentrations of gliclazide were determined by HPLC analysis. The interactions of gliclazide and bile acids were also investigated by molecular mechanics calculations (MM2). Probiotic bacteria significantly increased the permeability of gliclazide across the PAMPA membrane at all observed pH values while the total amount of gliclazide during incubation with bacteria was significantly reduced at pH 7.4, which could be a consequence of partial metabolism of the drug by enzymes of probiotic bacteria. Bile acids decreased the permeability of gliclazide through PAMPA membrane, with more pronounced effects of DCA, by forming more stable complexes with gliclazide. Given that probiotic bacteria and bile acids are naturally present in the gut and that each individual has a specific bacterial fingerprint, future research should extend the explanation of their effect on the gliclazide bioavailability and therapy individualization in in vivo conditions.

UV-Vis, FTIR, 1H, 13C NMR spectra and thermal studies of charge transfer complexes formed in the reaction of Gliclazide with π- and σ-electron acceptors

Charge transfer interactions (CT) between a gliclazide (GLC) donor and a picric acid (PA) π acceptor or iodine σ acceptor, were studied in a chloroform solution and in the solid state. UV-vis spectroscopy elucidated the formation of the complexes, and allowed determination of the stoichiometry, stability constants (K), and thermodynamic quantities (ΔG°, ΔH°, and ΔS°), and spectroscopic properties such as the molar extinction coefficient (ε_CT), oscillator strength (f), transition dipole moment (μ_EN), and ionization potential (I_p). Beer's law was obeyed over the 2-8 and 4-12 μg mL^-1 concentration ranges for GLC with PA (method A) and I₂ (method B), respectively, with correlation coefficients of 0.9986 and 0.9989. The limits of detection (LOD) and limits of quantification (LOQ) have also been reported. The 1:1 stoichiometric CT complexes were synthesized and characterized by FTIR, ¹H, and ¹³C NMR spectroscopy. The results indicated a favorable proton migration from PA to the donor molecule, and an interaction between the NH of GLC and iodine. Thermogravimetric analysis techniques (TGA/DTA) and differential scanning calorimetry (DSC) were used to determine the thermal stability of the synthesized CT complex. The kinetic parameters (ΔG*, ΔH*, and ΔS*) were calculated from thermal decomposition data using the Coats-Redfern method.

Altered metabolism of synthetic cannabinoid JWH-018 by human cytochrome P450 2C9 and variants

Synthetic cannabinoids (SCBs), synonymous with 'K2', 'Spice' or 'synthetic marijuana', are psychoactive drugs of abuse that frequently result in clinical effects and toxicity more severe than those classically associated with Δ⁹-tetrahydrocannabinol such as extreme agitation, hallucinations, supraventricular tachycardia, syncope, and seizures. JWH-018 is one of the earliest compounds identified in various SCB products, and our laboratory previously demonstrated that JWH-018 undergoes extensive metabolism by cytochromes P450 (P450), binds to, and activates cannabinoid receptors (CBRs). The major enzyme involved in the metabolism of JWH-018 is CYP2C9, a highly polymorphic enzyme found largely in the intestines and liver, with *1 being designated as the wild type, and *2 and *3 as the two most common variants. Three different major products have been identified in human urine and plasma: JWH-018 (ω)-OH, JWH-018 (ω-1)-OH(R), and JWH-018 (ω-1)-OH(S). The (ω-1)-OH metabolite of JWH-018 is a chiral molecule, and is thus designated as either (ω-1)-OH(R) or (ω-1)-OH(S). Here, in vitro enzyme kinetic assays performed with human recombinant CYP2C9 variants (*1, *2, and *3) revealed that oxidative metabolism by CYP2C9*3 resulted in significantly less formation of (ω)-OH and (ω-1)-OH metabolites. Surprisingly, CYP2C9*2 was roughly 3.6-fold more efficient as the CYP2C9*1 enzyme based on V_max/K_m, increasing the rate of JWH-018 metabolism and allowed for a much more rapid elimination. These results suggest that genetic polymorphisms of P450 enzymes result in the production of varying levels of biologically active JWH-018 metabolites in some individuals, offering a mechanistic explanation for the diverse clinical toxicity often observed following JWH-018 abuse.

Pharmacogenomics of CYP2C9: Functional and Clinical Considerations

CYP2C9 is the most abundant CYP2C subfamily enzyme in human liver and the most important contributor from this subfamily to drug metabolism. Polymorphisms resulting in decreased enzyme activity are common in the CYP2C9 gene and this, combined with narrow therapeutic indices for several key drug substrates, results in some important issues relating to drug safety and efficacy. CYP2C9 substrate selectivity is detailed and, based on crystal structures for the enzyme, we describe how CYP2C9 catalyzes these reactions. Factors relevant to clinical response to CYP2C9 substrates including inhibition, induction and genetic polymorphism are discussed in detail. In particular, we consider the issue of ethnic variation in pattern and frequency of genetic polymorphisms and clinical implications. Warfarin is the most well studied CYP2C9 substrate; recent work on use of dosing algorithms that include CYP2C9 genotype to improve patient safety during initiation of warfarin dosing are reviewed and prospects for their clinical implementation considered. Finally, we discuss a novel approach to cataloging the functional capabilities of rare 'variants of uncertain significance', which are increasingly detected as more exome and genome sequencing of diverse populations is conducted.

Gliclazide

Gliclazide is a second-generation oral hypoglycemic drug used for the treatment of noninsulin-dependent diabetes mellitus. It belongs to the sulfonylurea class that stimulates insulin secretion from pancreatic β-cells by inhibiting ATP-dependent potassium channels. Gliclazide also possesses unique antioxidant properties and other beneficial hemobiological effects. This profile represents a comprehensive description of the physical properties, chemical synthesis, spectroscopic characterization (FTIR, ¹H NMR, ¹³C NMR, UV, and single-crystal X-ray), methods of analysis, pharmacological actions, and pharmacokinetic and pharmacodynamic properties of the title drug.

Mechanistic insights into the effect of CYP2C9*2 and CYP2C9*3 variants on the 7-hydroxylation of warfarin

AIM

To evaluate the impact of CYP2C9*2 and CYP2C9*3 variants on binding and hydroxylation of warfarin.

MATERIALS & METHODS

Multiple linear regression model of warfarin pharmacokinetics was developed from the dataset of patients (n = 199). Pymol based in silico models were developed for the genetic variants.

RESULTS

CYP2C9*2 and CYP2C9*3 variants exhibited high warfarin/7-hydroxywarfarin (multiple linear regression model), dose-dependent disruption of hydrogen bonds with warfarin, dose-dependent increase in the distance between C7 of S-warfarin and Fe-O of CYP2C9, dose-dependent decrease in the glide scores (in silico).

CONCLUSION

CYP2C9*2 and CYP2C9*3 variants result in disruption of hydrogen bonding interactions with warfarin and longer distance between C7 and Fe-O thus impairing warfarin 7-hydroxylation due to lower binding affinity of warfarin. Original submitted 7 May 2014; Revision submitted 30 October 2014.

Molecular functionality of CYP2C9 polymorphisms and their influence on drug therapy

metabolizes approximately 20% of clinically used drugs, including the narrow therapeutic window drugs warfarin and phenytoin. More than 16,000 variants have been reported in the National Center for Biotechnology Information

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.