Individualization of treatments with drugs metabolized by CES1: combining genetics and metabolomics

Cannabidiol Inhibits In Vitro Human Liver Microsomal Metabolism of Remdesivir

Introduction: The year 2020 began with the world being flounced with a wave of novel coronavirus (severe acute respiratory syndrome coronavirus 2 [SARS-CoV-2]) disease, named COVID-19. Based on promising pre-clinical and clinical data, remdesivir (RDV) was the first drug to receive FDA approval and so far, it is the most common therapy for treatment of SARS-CoV-2/MERS-CoV. However, following intravenous administration, RDV metabolizes majorly by human liver carboxylesterase 1 (CES1) and marginally by the CYP3A4 enzyme in merely less than an hour. Its resultant active metabolite is a hydrophilic nucleoside with very limited accumulation within lung tissues. Therefore, there is a need to investigate strategies to overcome such premature metabolism issues and improve the antiviral efficacy of RDV at the target site. Objective: Considering the major CES1-mediated metabolism of RDV on systemic administration, we intend to explore the remarkable CES1 plus CYP3A4 inhibitory activity of cannabidiol (CBD) against in vitro microsomal metabolism of RDV to indicate its therapeutic potential as an adjuvant to RDV in the treatment and management of COVID-19. Methods: We investigated the in vitro human liver microsomal metabolism of RDV in the presence of two potential CES1 inhibitors-CBD and nelfinavir, and two standard CYP3A4 inhibitors-ritonavir (RITO) and cyclosporin A. The microsomal metabolism assay was further validated by using a well-characterized CYP3A4-selective substrate, midazolam (MDZ), in the presence of CBD and RITO. Results: Our findings depicted that RDV was rapidly and completely metabolized by human liver microsomes within 60 min. Coincubation with CBD substantially reduced microsomal metabolism of RDV and prolonged its in vitro half-life from 8.93 to 31.07 min. CBD showed significantly higher inhibition of RDV compared with known CES1 and CYP3A4 inhibitors. Inhibition of MDZ metabolism by CBD and RITO further validated the assay. Conclusions: The current study strongly suggests that CBD significantly inhibits human liver microsomal metabolism of RDV and extends its in vitro half-life. Thus, concomitant administration of CBD with RDV intravenous injection could be a promising strategy to prevent premature metabolism in COVID-19 patients.

Novel potent bifunctional carboxylesterase inhibitors based on a polyfluoroalkyl-2-imino-1,3-dione scaffold

An expanded series of alkyl 2-arylhydrazinylidene-3-oxo-3-polyfluoroalkylpropionates (HOPs) 3 was obtained via Cu(OAc)₂-catalyzed azo coupling. All were nanomolar inhibitors of carboxylesterase (CES), while moderate or weak inhibitors of acetylcholinesterase and butyrylcholinesterase. Steady-state kinetics studies showed that HOPs 3 are mixed type inhibitors of the three esterases. Molecular docking studies demonstrated that two functional groups in the structure of HOPs, trifluoromethyl ketone (TFK) and ester groups, bind to the CES active site suggesting subsequent reactions: formation of a tetrahedral adduct, and a slow hydrolysis reaction. The results of molecular modeling allowed us to explain some structure-activity relationships of CES inhibition by HOPs 3: their selectivity toward CES in comparison with cholinesterases and the high selectivity of pentafluoroethyl-substituted HOP 3p to hCES1 compared to hCES2. All compounds were predicted to have good intestinal absorption and blood-brain barrier permeability, low cardiac toxicity, good lipophilicity and aqueous solubility, and reasonable overall drug-likeness. HOPs with a TFK group and electron-donor substituents in the arylhydrazone moiety were potent antioxidants. All compounds possessed low cytotoxicity and low acute toxicity. Overall, a new promising type of bifunctional CES inhibitors has been found that are able to interact with the active site of the enzyme with the participation of two functional groups. The results indicate that HOPs have the potential to be good candidates as human CES inhibitors for biomedicinal applications.

Pharmacometabolomics Informs About Pharmacokinetic Profile of Methylphenidate

Carboxylesterase 1 (CES1) metabolizes methylphenidate and other drugs. CES1 gene variation only partially explains pharmacokinetic (PK) variability. Biomarkers predicting the PKs of drugs metabolized by CES1 are needed. We identified lipids in plasma from 44 healthy subjects that correlated with CES1 activity as determined by PK parameters of methylphenidate including a ceramide (q value = 0.001) and a phosphatidylcholine (q value = 0.005). Carriers of the CES1 143E allele had decreased methylphenidate metabolism and altered concentration of this phosphatidylcholine (q value = 0.040) and several high polyunsaturated fatty acid lipids (PUFAs). The half‐maximal inhibitory concentration (IC₅₀) values of chenodeoxycholate and taurocholate were 13.55 and 19.51 μM, respectively, consistent with a physiological significance. In silico analysis suggested that bile acid inhibition of CES1 involved both binding to the active and superficial sites of the enzyme. We initiated identification of metabolites predicting PKs of drugs metabolized by CES1 and suggest lipids to regulate or be regulated by this enzyme.

Genetic determinants of dabigatran safety (CES1 gene rs2244613 polymorphism) in the Russian population: multi-ethnic analysis

This study was aimed to investigate the prevalence of the CES1 gene (c.1168-33A > C, rs2244613) polymorphism among 12 different ethnic groups living in Russia to provide a basis for future clinical studies concerning genetic determinants of dabigatran safety. The study involved 1630 apparently healthy, unrelated, and chronic medication-free volunteers of both genders from 12 different ethnic groups in Russia: 136 Russians, 90 Avars, 50 Dargins, 46 Laks, 120 Kabardians, 112 Balkars, 244 Ossetians, 206 Mari, 204 Mordvinians, 238 Chuvashes, 114 Buryats and 70 Nanays. Genotyping was performed by using real-time polymerase chain reaction-based methods. The allelic prevalence of the ethnic groups was compared with Caucasus population participating in the RE-LY study. Statistically significant differences for the following gene polymorphism were found between all ethnic groups and RE-LY participants. Based on obtained results, it can be assumed that patients of all ethnic groups living in Russia taking dabigatran have a lower risk of bleeding.

Functional Study of Carboxylesterase 1 Protein Isoforms

Carboxylesterase 1 (CES1) is a primary human hepatic hydrolase involved in hydrolytic biotransformation of numerous medications. Considerable interindividual variability in CES1 expression and activity has been consistently reported. Four isoforms of the CES1 protein are produced by alternative splicing (AS). In the current study, the activity and expression of each CES1 isoform are examined using transfected cell lines, and CES1 isoform composition and its impact on CES1 activity in human livers are determined. In transfected cells, isoforms 3 and 4 show mRNA and protein expressions comparable to isoforms 1 and 2, but have significantly impaired activity when hydrolyzing enalapril and clopidogrel. In individual human liver samples, isoforms 1 and 2 are the major forms, contributing 73-90% of the total CES1 protein expression. In addition, the protein expression ratios of isoforms 1 and 2 to isoforms 3 and 4 are positively associated with CES1 activity in the liver, suggesting that CES1 isoform composition is a factor contributing to the variability in hepatic CES1 function. Further investigations of the regulation of CES1 AS would improve the understanding of CES1 variability and help develop a strategy to optimize the pharmacotherapy of many CES1 substrate medications.

Genetics of the patent ductus arteriosus (PDA) and pharmacogenetics of PDA treatment

Patent ductus arteriosus (PDA) is a frequent, complex, and difficult to treat clinical syndrome among preterm infants in the neonatal intensive care unit. In addition to known clinical risk factors, there are emerging data about genetic predisposition to PDA in both animal and human models. Clinical response and toxicity from drugs used to treat PDA are highly variable. Developmental and genetic aspects of pharmacokinetics and pharmacodynamics influence exposure and response to pharmacologic therapies. Given the variable efficacy and toxicity of known drug therapies, novel therapeutic targets for PDA treatment offer the promise of precision medicine. This review addresses the known genetic contributions to prolonged ductal patency, variability in response to drug therapy for PDA, and potential novel drug targets for future PDA treatment discovery.

Carboxylesterase 1A2 encoding gene with increased transcription and potential rapid drug metabolism in Asian populations

The carboxylesterase 1 gene (CES1) encodes a hydrolase implicated in the metabolism of commonly used drugs. CES1A2, a hybrid of CES1 and a CES1-like pseudogene, has a promoter that is weak in most individuals. However, some individuals harbor a promoter haplotype of this gene with two overlapping Sp1 sites that confer significantly increased transcription potentially leading to rapid drug metabolism. This CES1A2 haplotype has previously been reported to be common among Asians. Using polymerase chain reaction followed by sequencing, the present study examined variation in the promoter and 5' untranslated region of CES1A2 in 120 Han Chinese and 120 Japanese people enrolled in the 1000 Genomes Project. We identified 11 single nucleotide variations, two of which were novel, in 145 of the individuals who were found to carry CES1A2. Alignment analysis indicated that the CES1A2 haplotype with the overlapping Sp1 sites has been generated by incorporation of a segment of CES1. All minor allele frequencies were equal to or below 0.022 and the frequencies of the minor haplotypes were up to 40-fold lower than previously reported, including that of the haplotype with the overlapping Sp1 sites. This information is novel and suggests that the pharmacogenetic relevance of CES1A2 is limited in Asians.

A Comprehensive Functional Assessment of Carboxylesterase 1 Nonsynonymous Polymorphisms

Carboxylesterase 1 (CES1) is the predominant human hepatic hydrolase responsible for the metabolism of many clinically important medications. CES1 expression and activity vary markedly among individuals; and genetic variation is a major contributing factor to CES1 interindividual variability. In this study, we comprehensively examined the functions of CES1 nonsynonymous single nucleotide polymorphisms (nsSNPs) and haplotypes using transfected cell lines and individual human liver tissues. The 20 candidate variants include CES1 nsSNPs with a minor allele frequency >0.5% in a given population or located in close proximity to the CES1 active site. Five nsSNPs, including L40Ter (rs151291296), G142E (rs121912777), G147C (rs146456965), Y170D (rs148947808), and R171C (rs201065375), were loss-of-function variants for metabolizing the CES1 substrates clopidogrel, enalapril, and sacubitril. In addition, A158V (rs202121317), R199H (rs2307243), E220G (rs200707504), and T290M (rs202001817) decreased CES1 activity to a lesser extent in a substrate-dependent manner. Several nsSNPs, includingL40Ter (rs151291296), G147C (rs146456965), Y170D (rs148947808), and R171C (rs201065375), significantly reduced CES1 protein and/or mRNA expression levels in the transfected cells. Functions of the common nonsynonymous haplotypes D203E-A269S and S75N-D203E-A269S were evaluated using cells stably expressing the haplotypes and a large set of the human liver. Neither CES1 expression nor activity was affected by the two haplotypes. In summary, this study revealed several functional nsSNPs with impaired activity on the metabolism of CES1 substrate drugs. Clinical investigations are warranted to determine whether these nsSNPs can serve as biomarkers for the prediction of therapeutic outcomes of drugs metabolized by CES1.

Reappraisal of the genetic diversity and pharmacogenetic assessment of CES1

The CES1 gene encodes a hydrolase that metabolizes important drugs. Variants generated by exchange of segments with CES1P1 complicate genotyping of CES1. Using a highly specific procedure we examined DNA samples from 200 Caucasians and identified 46 single nucleotide variants (SNVs) in CES1 and 21 SNVs in CES1A2, a hybrid composed of CES1 and CES1P1. Several of these SNVs were novel. The frequencies of SNVs with a potential functional impact were below 0.02 suggesting limited pharmacogenetic potential for CES1 genotyping. In silico PCR revealed that the majority of the primer pairs for amplification of CES1 or CES1A2 in three previous studies lacked specificity, which partially explains a limited overlap with our findings.

Pharmacodynamic Impact of Carboxylesterase 1 Gene Variants in Patients with Congestive Heart Failure Treated with Angiotensin-Converting Enzyme Inhibitors

BACKGROUND

Variation in the carboxylesterase 1 gene (CES1) may contribute to the efficacy of ACEIs. Accordingly, we examined the impact of CES1 variants on plasma angiotensin II (ATII)/angiotensin I (ATI) ratio in patients with congestive heart failure (CHF) that underwent ACEI dose titrations. Five of these variants have previously been associated with drug response or increased CES1 expression, i.e., CES1 copy number variation, the variant of the duplicated CES1 gene with high transcriptional activity, rs71647871, rs2244613, and rs3815583. Additionally, nine variants, representatives of CES1Var, and three other CES1 variants were examined.

METHODS

Patients with CHF, and clinical indication for ACEIs were categorized according to their CES1 genotype. Differences in mean plasma ATII/ATI ratios between genotype groups after ACEI dose titration, expressed as the least square mean (LSM) with 95% confidence intervals (CIs), were assessed by analysis of variance.

RESULTS

A total of 200 patients were recruited and 127 patients (63.5%) completed the study. The mean duration of the CHF drug dose titration was 6.2 (SD 3.6) months. After ACEI dose titration, there was no difference in mean plasma ATII/ATI ratios between subjects with the investigated CES1 variants, and only one previously unexplored variation (rs2302722) qualified for further assessment. In the fully adjusted analysis of effects of rs2302722 on plasma ATII/ATI ratios, the difference in mean ATII/ATI ratio between the GG genotype and the minor allele carriers (GT and TT) was not significant, with a relative difference in LSMs of 0.67 (95% CI 0.43-1.07; P = 0.10). Results of analyses that only included enalapril-treated patients remained non-significant after Bonferroni correction for multiple parallel comparisons (difference in LSM 0.60 [95% CI 0.37-0.98], P = 0.045).

CONCLUSION

These findings indicate that the included single variants of CES1 do not significantly influence plasma ATII/ATI ratios in CHF patients treated with ACEIs and are unlikely to be primary determinants of ACEI efficacy.

Association of Oseltamivir Activation with Gender and Carboxylesterase 1 Genetic Polymorphisms

Oseltamivir, an inactive anti-influenza virus prodrug, is activated (hydrolysed) in vivo by carboxylesterase 1 (CES1) to its active metabolite oseltamivir carboxylate. CES1 functions are significantly associated with certain CES1 genetic variants and some non-genetic factors. The purpose of this study was to investigate the effect of gender and several CES1 genetic polymorphisms on oseltamivir activation using a large set of individual human liver samples. CES1-mediated oseltamivir hydrolysis and CES1 genotypes, including the G143E (rs71647871), rs2244613, rs8192935, the -816A>C (rs3785161) and the CES1P1/CES1P1VAR, were determined in 104 individual human livers. The results showed that hepatic CES1 protein expression in females was 17.3% higher than that in males (p = 0.039), while oseltamivir activation rate in the livers from female donors was 27.8% higher than that from males (p = 0.076). As for CES1 genetic polymorphisms, neither CES1 protein expression nor CES1 activity on oseltamivir activation was significantly associated with the rs2244613, rs8192935, -816A>C or CES1P1/CES1P1VAR genotypes. However, oseltamivir hydrolysis in the livers with the genotype 143G/E was approximately 40% of that with the 143G/G genotype (0.7 ± 0.2 versus 1.8 ± 1.1 nmole/mg protein/min, p = 0.005). In summary, the results suggest that hepatic oseltamivir activation appears to be more efficient in females than that in males, and the activation can be impaired by functional CES1 variants, such as the G143E. However, clinical implication of CES1 gender differences and pharmacogenetics in oseltamivir pharmacotherapy warrants further investigations.

Sacubitril Is Selectively Activated by Carboxylesterase 1 (CES1) in the Liver and the Activation Is Affected by CES1 Genetic Variation

Sacubitril was recently approved by the Food and Drug Administration for use in combination with valsartan for the treatment of patients with heart failure with reduced ejection fraction. As a prodrug, sacubitril must be metabolized (hydrolyzed) to its active metabolite sacubitrilat (LBQ657) to exert its intended therapeutic effects. Thus, understanding the determinants of sacubitril activation will lead to the improvement of sacubitril pharmacotherapy. The objective of this study was to identify the enzyme(s) responsible for the activation of sacubitril, and determine the impact of genetic variation on sacubitril activation. First, an incubation study of sacubitril with human plasma and the S9 fractions of human liver, intestine, and kidney was conducted. Sacubitril was found to be activated by human liver S9 fractions only. Moreover, sacubitril activation was significantly inhibited by the carboxylesterase 1 (CES1) inhibitor bis-(p-nitrophenyl) phosphate in human liver S9. Further incubation studies with recombinant human CES1 and carboxylesterase 2 confirmed that sacubitril is a selective CES1 substrate. The in vitro study of cell lines transfected with wild-type CES1 and the CES1 variant G143E (rs71647871) demonstrated that G143E is a loss-of-function variant for sacubitril activation. Importantly, sacubitril activation was significantly impaired in human livers carrying the G143E variant. In conclusion, sacubitril is selectively activated by CES1 in human liver. The CES1 genetic variant G143E can significantly impair sacubitril activation. Therefore, CES1 genetic variants appear to be an important contributing factor to interindividual variability in sacubitril activation, and have the potential to serve as biomarkers to optimize sacubitril pharmacotherapy.

Prognostic impact of carboxylesterase 1 gene variants in patients with congestive heart failure treated with angiotensin-converting enzyme inhibitors

OBJECTIVE

Most angiotensin-converting enzyme inhibitors (ACEIs) are prodrugs activated by carboxylesterase 1 (CES1). We investigated the prognostic importance of CES1 gene (CES1) copy number variation and the rs3815583 single-nucleotide polymorphism in CES1 among ACEI-treated patients with congestive heart failure (CHF).

METHODS

Danish patients with chronic CHF enrolled in the previously reported Echocardiography and Heart Outcome Study were categorized according to their CES1 variants and followed up for up to 10 years. Risk for cardiovascular death and all-cause death was modeled by Cox proportional hazard analyses.

RESULTS

A total of 491 ACEI-treated patients were included in the analyses. After a mean follow-up of 5.5 years, we found no difference in the risk for cardiovascular death and all-cause death between patients having three [hazard ratios (HRs) 1.06 (95% confidence interval (CI) 0.77-1.45) and 1.16 (95% CI 0.88-1.52)] or four [HRs 0.88 (95% CI 0.39-2.01) and 1.37 (95% CI 0.74-2.54)] CES1 copies and those with two copies, respectively. Similarly, no difference in the risk for cardiovascular and all-cause death was found for patients heterozygous [HRs 0.91 (95% CI 0.70-1.19) and 0.88 (95% CI 0.69-1.12)] or homozygous [HRs 0.58 (95% CI 0.30-1.15) and 0.82 (95% CI 0.48-1.39)] for the rs3815583 minor allele versus patients homozygous for the major allele. The active promoter of CES1A2 and the rs71647871 single-nucleotide polymorphism minor allele were detected at very low frequencies.

CONCLUSION

This study did not support the use of CES1 copy number variation or rs3815583 as a predictor of fatal outcomes in ACEI-treated patients with CHF.