DMD Fast Forward. Published on December 23, 2019 as DOI: 10.1124/dmd.119.089680 This article has not been copyedited and formatted. The final version may differ from this version.

DMD # 89680 1

Carboxylesterase 1 and Precision Pharmacotherapy: Pharmacogenetics and Non-

Genetic Regulators

Lucy Her and Hao-Jie Zhu

Department of Clinical Pharmacy (L.H. and H.-J.Z.), University of Michigan, Ann Arbor,

Michigan

Downloaded from dmd.aspetjournals.org at ASPET Journals on October 1, 2021 DMD Fast Forward. Published on December 23, 2019 as DOI: 10.1124/dmd.119.089680 This article has not been copyedited and formatted. The final version may differ from this version.

DMD # 89680 2

Running title: 1 and Precision Pharmacotherapy

Corresponding author: Hao-Jie Zhu Ph.D. Department of Clinical Pharmacy University of Michigan College of Pharmacy 428 Church Street, Room 4565 NUB Ann Arbor, MI 48109-1065 Tel: 734-763-8449, E-mail: [email protected]

Number of text pages: 25 Number of tables: 4 Downloaded from Number of figures: 5 Number of references: 102 Number of words in Abstract: 223 Number of words in Introduction: 166

Number of words in Discussion: 287 dmd.aspetjournals.org

Keywords: , pharmacogenetics, precision pharmacotherapy

List of Abbreviations: at ASPET Journals on October 1, 2021 ADHD attention deficit hyperactivity disorder API active pharmaceutical ingredient AT I angiotensin i AT II angiotensin ii AUC area under the curve BCS biopharmaceutics classification system CES1 carboxylesterase 1 CES2 CI confidence interval Cmax maximum plasma concentration CNVs copy number variations Clinical Pharmacogenetics Implementation CPIC Consortium CYPs cytochrome p450 DL-EPH dl-ethylphenidate DMEs drug-metabolizing FDA Food and Drug Administration GWAS genome-wide association study INH isoniazid LD loading dose LOF loss-of-function DMD Fast Forward. Published on December 23, 2019 as DOI: 10.1124/dmd.119.089680 This article has not been copyedited and formatted. The final version may differ from this version.

DMD # 89680 3

MAF minor allele frequency MD maintenance dose OR odds ratio PCI percutaneous coronary intervention PD pharmacodynamics PK pharmacokinetics PNPA p-nitrophenyl acetate PTM post-transcriptional modification SNPs single-nucleotide polymorphisms T2DM Type 2 Diabetes Mellitus

Downloaded from dmd.aspetjournals.org at ASPET Journals on October 1, 2021 DMD Fast Forward. Published on December 23, 2019 as DOI: 10.1124/dmd.119.089680 This article has not been copyedited and formatted. The final version may differ from this version.

DMD # 89680 4

Abstract

Carboxylesterase 1 (CES1) is the most abundant drug-metabolizing in human livers, comprising approximately 1% of the entire liver proteome. CES1 is responsible for

80-95% of total hydrolytic activity in the liver and plays a crucial role in the metabolism of a wide range of drugs (especially ester-prodrugs), pesticides, environmental pollutants, and endogenous compounds. Expression and activity of CES1 vary markedly among individuals, which is a major contributing factor to interindividual variability in the Downloaded from pharmacokinetics (PK) and pharmacodynamics (PD) of drugs metabolized by CES1.

Both genetic and non-genetic factors contribute to CES1 variability. Here, we discuss dmd.aspetjournals.org genetic polymorphisms, including single-nucleotide polymorphisms (SNPs) and copy number variants, and non-genetic contributors, such as developmental status, genders, and drug-drug interactions, that could influence CES1 functionality and the PK and PD at ASPET Journals on October 1, 2021 of CES1 substrates. Currently, the loss-of-function SNP G143E (rs71647871) is the only clinically significant CES1 variant identified to date, and alcohol is the only potent CES1 inhibitor that could alter the therapeutic outcomes of CES1 medications.

However, G143E and alcohol can only explain a small portion of the interindividual variability in the CES1 function. A better understanding of the regulation of CES1 expression and activity and identification of biomarkers for CES1 function in vivo could lead to the development of a precision pharmacotherapy strategy to improve the efficacy and safety of many CES1 substrate drugs.

DMD Fast Forward. Published on December 23, 2019 as DOI: 10.1124/dmd.119.089680 This article has not been copyedited and formatted. The final version may differ from this version.

DMD # 89680 5

Significance Statement

Carboxylesterase 1 (CES1) is a major phase I drug-metabolizing enzyme responsible for 80-95% of total hydrolytic activity in human livers. CES1 plays a crucial role in metabolizing a wide range of drugs, pesticides, environmental pollutants, and endogenous compounds. The clinical relevance of CES1 has been well demonstrated in various clinical trials with methylphenidate, oseltamivir, and clopidogrel. Here, we discuss genetic polymorphisms and non-genetic contributors that influence CES1 Downloaded from functionality and the PK and PD of drugs metabolized by CES1.

dmd.aspetjournals.org at ASPET Journals on October 1, 2021 DMD Fast Forward. Published on December 23, 2019 as DOI: 10.1124/dmd.119.089680 This article has not been copyedited and formatted. The final version may differ from this version.

DMD # 89680 6

1. Introduction

Carboxylesterase 1 (CES1) is a phase I drug-metabolizing enzyme (DME) responsible for 80-95% of total hydrolytic activity in the liver (Imai et al., 2006); it metabolizes a wide range of drugs, pesticides, environmental pollutants, and endogenous compounds, including lipid esters (Table 1). CES1-mediated metabolism can lead to the biotransformation of a pharmacologically active drug into its inactive metabolite, as exemplified by methylphenidate hydrolysis in the liver. CES1 also plays Downloaded from an important role in activating prodrugs since most ester-containing prodrugs are exclusively dependent on CES1 for their activation. The clinical relevance of CES1 has dmd.aspetjournals.org been well demonstrated in various clinical trials with oseltamivir, methylphenidate, and clopidogrel (Zhu et al., 2008; Tarkiainen et al., 2012; Lewis et al., 2013; Jiang et al.,

2016). Recent studies have also revealed that CES1 acts as a cholesteryl ester at ASPET Journals on October 1, 2021 in lipid metabolism in human macrophages and hepatocytes, and suggest

CES1 as a potential drug target for the treatment of metabolic diseases, such as diabetes and atherosclerosis (Dolinsky et al., 2004; Zhao et al., 2007; Ghosh et al.,

2010; Ross et al., 2010; Lian et al., 2018b).

2. Importance of CES1 in

CES1 plays an important role in metabolizing many clinically significant medications, especially the ester-prodrugs (Table 1). A prodrug refers to an inactive drug molecule that needs to be enzymatically biotransformed in vivo to its active metabolite to produce its intended pharmacological effect (Rautio et al., 2008). Prodrug design offers an attractive method to overcome the issue of low bioavailability for

Biopharmaceutics Classifications System (BCS) class III drug molecules. Drug DMD Fast Forward. Published on December 23, 2019 as DOI: 10.1124/dmd.119.089680 This article has not been copyedited and formatted. The final version may differ from this version.

DMD # 89680 7 molecules can be categorized into four BCS classes based on permeability and solubility, and a BCS class III substance is a hydrophilic compound with low permeability and high solubility (Shah and Amidon, 2014). In particular, hydrophilic compounds with –

OH or –COOH functional groups usually have difficulty being absorbed into the body, and drug developers often mask these functional groups using an ester-prodrug design.

The prodrug market has been growing: 20% of drugs approved in 2015 were prodrugs, compared to ~6% of all currently approved drugs (Rautio et al., 2017). Downloaded from

Two major assumptions behind the ester-prodrug design are that prodrugs are

rapidly activated via unspecific in the body, and that the interindividual dmd.aspetjournals.org variability in activating a prodrug is clinically insignificant. These incorrect assumptions may have stemmed from the fact that many hydrolytic enzymes exist in the body, such

as CES1, CES2, , , paraoxonases, and at ASPET Journals on October 1, 2021 arylesterase. However, these differ in their tissue-specific expression, cellular localization, and most importantly, substrate selectivity (Fukami and Yokoi,

2012). In humans, CES1 is highly abundant in the liver and expressed to a less extent in the lung and brain; CES1 expression is considered negligible in the human intestine, kidney, and plasma. CES1 is substrate-selective towards carboxyl esters with a large ethyl group and a small alcohol group. In comparison, CES2, another major carboxylesterase in humans, is highly expressed in the intestine, kidney, and liver, and is more efficient at metabolizing compounds with a small ethyl group and a large alcohol group (Jewell et al., 2007). Numerous in vivo and in vitro studies have demonstrated the specificity of CES1, and many CES1 substrates cannot be metabolized by other esterases (Table 3). DMD Fast Forward. Published on December 23, 2019 as DOI: 10.1124/dmd.119.089680 This article has not been copyedited and formatted. The final version may differ from this version.

DMD # 89680 8

CES1 expression and activity vary significantly among individuals (Wang et al.,

2016a); this variability could result in treatment failure and unexpected adverse effects of CES1 substrate drugs. A better understanding of the genetic and non-genetic factors contributing to CES1 variability will improve the design and clinical use of many drugs that are metabolized (deactivated/activated) by CES1.

3. Pharmacogenetics of drug-metabolizing enzymes Downloaded from Traditionally, fixed-dose regimens have been used for most medications.

However, different individuals taking the same dose of medication do not necessarily

achieve the same drug exposure and hence drug response. More individualized, patient- dmd.aspetjournals.org centered dosing regimens have been developed based on a patient’s characteristics, such as renal clearance, liver function, body weight, and surface area (DiPiro, 2017). In

addition, genetic polymorphisms of DMEs have been found to play an important role in at ASPET Journals on October 1, 2021 the response to pharmacotherapy, and pharmacogenomics has been increasingly utilized in the clinic to improve the efficacy and safety of drug treatment. DMEs serve to primarily detoxify digested xenobiotics through four general mechanisms: hydrolysis

(e.g., carboxylesterase), reduction (e.g., carbonyl reductase), oxidation (e.g., cytochrome P450), and conjugation (e.g., UDP-glucuronosyltransferase) (Foti and

Dalvie, 2016). The expression and activity of DMEs vary significantly among individuals, and studying pharmacogenomics of DMEs is one means of better understanding interindividual variability in the PK and PD of a drug. For example, the active metabolite of irinotecan, SN-38, is primarily metabolized by the enzyme UDP glucuronosyltransferase family 1 member A1 (UGT1A1 enzyme) (Ando et al., 2000). If a patient carries the common UGT1A1*28 polymorphism, the decrease it causes in DMD Fast Forward. Published on December 23, 2019 as DOI: 10.1124/dmd.119.089680 This article has not been copyedited and formatted. The final version may differ from this version.

DMD # 89680 9

UGT1A1 enzymatic activity would impede the metabolism of SN-38, leading towards the accrual of toxic concentrations. Accordingly, the Food and Drug Administration (FDA) recommended that patients with UGT1A1*28/*28 start irinotecan at a lower dose

(Innocenti et al., 2004). However, given that both genetic and environmental factors contribute to DME function, we should also pay close attention to non-genetic contributors when studying the variability of DMEs.

4. CES1 pharmacogenetics Downloaded from

Although CES1 plays a critical role in the metabolism of many clinically important

medications, CES1 pharmacogenetics is understudied relative to other major DMEs dmd.aspetjournals.org

(e.g., CYP450s). CES1 is encoded by the CES1 , consisting of 14 exons located on 16q13-q22.1. CES1 VAR is a variation of the CES1 gene that differs in

exon 1 DNA sequences and has an average minor allele frequency (MAF) of 17%. at ASPET Journals on October 1, 2021

Although one study claimed that CES1 VAR mRNA was undetectable (Fukami et al.,

2008), an in vitro human liver study showed that the protein expressions of CES1 and

CES1 VAR were not statistically different (Wang et al., 2016a). CES1P1 is a pseudogene due to a premature stop codon in exon 4 and lies tail-to-tail with CES1

(Figure 1) (Wang et al., 2016a). Interestingly, a CES1P1 variant named CES1P1 VAR is a functional coding gene with a DNA sequence identical to CES1 VAR. However, the transcription efficiency of CES1P1 VAR is only 2% of that of CES1, due to the transcription factor Sp1, and the enhancer-binding protein C/EBPα preferring to bind to the CES1 promoter over the CES1P1 VAR promoter (Hosokawa et al., 2008; Yoshimura et al., 2008). Due to the existence of the CES1 VAR and CES1P1 VAR variants, four

CES1/CES1P1 haplotypes can be formed (Figure 1). In addition to these structural DMD Fast Forward. Published on December 23, 2019 as DOI: 10.1124/dmd.119.089680 This article has not been copyedited and formatted. The final version may differ from this version.

DMD # 89680 10 variations, there are over 7000 CES1 SNPs registered in the NCBI SNP database, and approximately 300 of them have MAFs over 1%. These common CES1 variants (MAF >

1%) are distributed in various regions of the gene, including 13 in 5’-UTR and 3’-UTRs,

14 in exons, and 308 in introns. Of the exonic SNPs, 12 are non-synonymous SNPs and two are synonymous SNPs. In the following section, we discuss the clinical findings and mechanistic bases of functional CES1 variants identified to date. Downloaded from 4.1 Pharmacogenetics of the first loss-of-function CES1 variant G143E

(rs71647871) dmd.aspetjournals.org In SNP notation, G143E indicates an change from glycine to glutamic acid at amino acid position 143. G143E is also termed 428G>A, indicating that the nucleotide guanine is changed to adenine at position 428 of the CES1 mRNA at ASPET Journals on October 1, 2021 (DiPiro, 2017). The MAF of G143E is 3.7%, 4.3%, and 2%, in White, Hispanic, and

African American populations, respectively, while the SNP is extremely rare in Asian populations (Zhu et al., 2008; Suzaki et al., 2013a).

G143E is a non-conservative amino acid substitution located near the active-site triad residues of CES1 (serine 221, glutamic acid 354, and histidine 468). Serine hydrolases shares similar catalytic mechanism involving (1) nucleophilic attack from oxygen in the serine residue on a substrate ester bond, (2) formation of a tetrahedral intermediate where the deprotonated oxygen is stabilized via an , (3) formation of an acyl-enzyme intermediate, and (4) water-catalyzed hydrolysis (Satoh and Hosokawa, 2006). For CES1 to maintain its enzymatic function, the and oxyanion hole need to be conserved (Zhu et al., 2008; Arena de Souza et al., 2015).

The change from glycine (hydrophobic residue) to glutamic acid (electrostatic residue) at DMD Fast Forward. Published on December 23, 2019 as DOI: 10.1124/dmd.119.089680 This article has not been copyedited and formatted. The final version may differ from this version.

DMD # 89680 11 codon 143 disrupts the hydrophobicity needed for the oxyanion hole (Gly 141-131), resulting in a complete loss of function of CES1. The G143E is only CES1 SNP that has been subjected to in vitro kinetics studies, in which the variant exhibited null catalytic activity on all tested CES1 substrates except for oseltamivir (Table 2). The Vmax of

G143E on oseltamivir hydrolysis was 37 nmol/min/mg with catalytic efficiency of 17.2

μl/min/mg protein — this was approximately 16% of wild type CES1 catalytic efficiency

(Zhu and Markowitz, 2009). Downloaded from

4.1.1 Discovery of G143E and its impacts on methylphenidate PK and PD

G143E is the first loss-of-function (LOF) variant known for CES1 and was dmd.aspetjournals.org originally discovered in a methylphenidate (Ritalin®) PK study in healthy volunteers.

Methylphenidate is a central nervous system stimulant, the most commonly prescribed

medication for ADHD treatment and has high abuse potential when used with alcohol at ASPET Journals on October 1, 2021

(Neos Therapeutics, 2017). Its drug comes as a racemic mixture of d- and l- methylphenidate hydrochloride; d-methylphenidate is approximately 10 times more pharmacologically potent than l-methylphenidate (Heal and Pierce, 2006).

Methylphenidate is metabolized by de-esterification via CES1 to ritalinic acid, an inactive metabolite which accounts for approximately 80% of the recovered dose in human urine (Figure 2) (Laizure et al., 2013; Neos Therapeutics, 2017). In 2007, a prospective single-dose (0.3 mg/kg) PK study was conducted in twenty healthy volunteers to examine the DDI between methylphenidate and alcohol (Patrick et al.,

2007). During this study, the researchers unexpectedly found a participant that showed significantly elevated pharmacokinetic parameters (e.g., AUC, Cmax) of methylphenidate.

Specifically, dl-methylphenidate Cmax was seven times higher and l-methylphenidate DMD Fast Forward. Published on December 23, 2019 as DOI: 10.1124/dmd.119.089680 This article has not been copyedited and formatted. The final version may differ from this version.

DMD # 89680 12

Cmax was 100-fold higher in this poor metabolizer compared to the rest participants.

Later analysis found that this poor metabolizer carried the G143E polymorphism in

CES1 and the D260fs polymorphism in CES1P1 (Zhu et al., 2008). This study also concluded that while CES1 metabolism is substantially stereoselective towards l- methylphenidate, d-methylphenidate metabolism is also significantly impacted by CES1 dysfunction. Downloaded from Following the discovery of the G143E variant, a retrospective study was conducted to examine the methylphenidate response in Hungarian ADHD patients with

G143E (n=7) and non-carrier patients (n=115). Even though the CES1 genotype could dmd.aspetjournals.org not explain the entire interindividual variability between responders (n=90) and non- responders (n=32), the study demonstrated an association between G143E

polymorphism and methylphenidate dose reduction: five responders who had the G143E at ASPET Journals on October 1, 2021 polymorphism required lower doses of methylphenidate for symptom reduction (0.410 vs

0.572 mg/kg, P=0.022) (Nemoda et al., 2009). In 2017, a healthy volunteer study confirmed the significance of G143E in the PK of methylphenidate. In this open-label, prospective clinical trial (n=22), study participants carrying the G143E SNP (n=6) had approximately 152.4% higher median AUC of d-methylphenidate (53.3ng×ml-1×h-1) than the non-carrier group (21.4 ng×ml-1×h-1) (P<0.0001) (Stage et al., 2017a).

The above studies suggest that G143E carriers may be at high risk of being exposed to a toxic methylphenidate concentration. This result is clinically impactful because methylphenidate is considered as the first-line pharmacotherapy for ADHD, with approximately 40 million prescriptions dispensed every year (Schubert et al., 2010).

This result could potentially explain why many patients have an unsatisfactory response DMD Fast Forward. Published on December 23, 2019 as DOI: 10.1124/dmd.119.089680 This article has not been copyedited and formatted. The final version may differ from this version.

DMD # 89680 13 to the treatment. Further clinical studies in ADHD patients with larger sample sizes are needed to fully understand the effect of CES1 variants on the efficacy and toxicity of methylphenidate, and how methylphenidate doses should be adjusted based on a patient’s CES1 genotypes.

4.1.2 G143E and clopidogrel (Plavix®)

Clopidogrel is a P2Y12 inhibitor and has several clinical indications, including myocardial infarction prophylaxis, cerebrovascular accident prophylaxis, and peripheral Downloaded from arterial occlusive disease prophylaxis. Clopidogrel is usually considered as the first-line antiplatelet agent due to its proven efficacy and cost-effectiveness (Wiviott et al., 2007; dmd.aspetjournals.org Wallentin et al., 2009; Roe et al., 2012). Clopidogrel is a non-ester-prodrug that needs to be activated by two oxidation reactions via several CYPs (Figure 3). CYP2C19 pharmacogenetics and its impact on clopidogrel activation have been extensively at ASPET Journals on October 1, 2021 studied. The Clinical Pharmacogenetics Implementation Consortium (CPIC) guideline and the FDA both recommend intermediate and poor metabolizers of CYP2C19 to use an alternative antiplatelet agent, such as ticagrelor or prasugrel (Scott et al., 2013).

Clopidogrel, its intermediate, and active metabolites are all CES1 substrates and metabolized by CES1 to inactive hydrolytic metabolites (Figure 3). Approximately, 85% of clopidogrel is hydrolyzed by CES1, and only 15% clopidogrel enters the CYPs- mediated activation pathway (Zhu et al., 2013). Thus, patients with CES1 dysfunction would have a higher concentration of clopidogrel active metabolite compared to normal

CES1 metabolizers when taking the same dose. However, the impact of CES1 on the

PK and PD of clopidogrel is less studied than the impacts of CYPs.

Two clinical trials support that CES1 G143E carriers have significantly higher plasma concentrations of clopidogrel active metabolite compared to non-carriers. A DMD Fast Forward. Published on December 23, 2019 as DOI: 10.1124/dmd.119.089680 This article has not been copyedited and formatted. The final version may differ from this version.

DMD # 89680 14 retrospective sub-analysis was performed on participants of the Pharmacogenomics of

Antiplatelet Intervention (PAPI) Study (n=506) and on clopidogrel-treated patients at

Sinai Hospital (n=350) to examine the effect of CES1 G143E on clopidogrel metabolism.

Study participants received a 300mg loading dose of clopidogrel followed by a 75mg maintenance dose for six days, and platelet aggregation was measured as a PD marker.

A 50% higher active metabolite concentration was observed in G143E carriers (n=7,

30.3 ng/ml) compared to non-carriers (n=499, 19.0 ng/mL) (P=0.001). In addition, the Downloaded from inhibition of adenosine diphosphate (ADP)-induced platelet aggregation was 24% higher in G143E carriers (reduced to 71% from baseline) relative to non-carriers (reduced to dmd.aspetjournals.org 57% from baseline) (P=0.003) (Lewis et al., 2013; Bozzi et al., 2016; Jiang et al., 2016).

Another prospective, single-dose, healthy volunteer (n=22) clinical study was conducted by Tarkiainen et al. to determine the effect of CES1 G143E on clopidogrel metabolism. at ASPET Journals on October 1, 2021

The authors found that the AUC0–∞ ratio of clopidogrel carboxylic acid (inactive metabolite (1) in Figure 3) to clopidogrel was 53% less in G143E carriers (n=10) than non-carriers (n=12) (P=0.009). The G143E carriers also exhibited significantly higher plasma concentrations of the parent compound clopidogrel (P = 0.004) and its active metabolite (P = 0.009) compared to non-carriers. In agreement with the PK findings, the average inhibition of P2Y12-mediated platelet aggregation in the carriers was 19 percentage points higher than in non-carriers (P = 0.036) (Zhu et al., 2013; Tarkiainen et al., 2015a). The findings of the above two studies are especially important for patients on triple antithrombotic therapy with a high bleeding risk (Mehta et al., 2001; Steinhubl et al., 2002; Shmyr et al., 2017). Clopidogrel dose adjustment may be necessary to prevent potential toxicity (i.e., bleeding) in patients with CES1 dysfunction. DMD Fast Forward. Published on December 23, 2019 as DOI: 10.1124/dmd.119.089680 This article has not been copyedited and formatted. The final version may differ from this version.

DMD # 89680 15

4.1.3 G143E and angiotensin-converting enzyme inhibitors (ACEIs)

Angiotensin-converting enzyme inhibitors (ACEIs) are generally considered to be the first-line therapy for heart failure and hypertension, and approximately 150 million

ACEI prescriptions are filled in the US annually (Mahmoudpour et al., 2015). Currently, eight out of ten FDA approved ACEIs are ester-containing prodrugs, and all ACEI prodrugs need to be activated by CES1 in order to exert their intended therapeutic effects (Chaturvedi, 2004; Yancy et al., 2017). The activation is essential for the Downloaded from pharmacological effects as the active metabolites are 10-1,000 times more potent than their prodrug forms (Foye et al., 2013). Therefore, patients with CES1 dysfunction would dmd.aspetjournals.org have a lower concentration of the ACEI active metabolite relative to normal CES1 metabolizers (Figure 4).

A prospective, single-dose pharmacokinetic clinical study was conducted in at ASPET Journals on October 1, 2021 healthy volunteers to examine the effect of the G143E variant on the activation of the

ACEI prodrugs enalapril and quinapril. The AUC0–∞ of the enalapril active metabolite enalaprilat was found to be 20% lower in the G143E carriers (n=10) than in non-carriers

(n=12) (P=0.049) (Tarkiainen et al., 2015b). This finding is consistent with an in vitro study that showed that enalapril activation was impaired in liver samples carrying the

G143E variant (Wang et al., 2016a). However, the AUCs0-∞ of the quinapril and its active metabolite (quinaprilat) were not significantly different between carriers and non-carriers

(P=0.114). Further investigations are warranted to fully understand the effect of CES1 variants on the PK and PD of ACEI prodrugs.

4.1.4 G143E and oseltamivir (Tamiflu®) DMD Fast Forward. Published on December 23, 2019 as DOI: 10.1124/dmd.119.089680 This article has not been copyedited and formatted. The final version may differ from this version.

DMD # 89680 16

Oseltamivir is an antiviral drug that has an FDA indication for influenza types A and B infections. Even though oseltamivir is rarely effective due to its specific administration requirement (i.e., this medication should be taken within 2 days of onset of symptoms in order to reduce flu duration by approximately one day), oseltamivir remains one of the most prescribed drug products due to flu epidemics (Singh et al.,

2003; Dahlgren et al., 2018). As an ester-prodrug, oseltamivir needs to be activated by

CES1 into its active metabolite, oseltamivir carboxylate (Shi et al., 2006). An in vitro Downloaded from study based on cell lines stably transfected with CES1 variants suggested the G143E

SNP markedly impaired CES1 activity in oseltamivir activation (Zhu and Markowitz, dmd.aspetjournals.org 2009).

To examine the effect of G143E on oseltamivir PK and activation, a prospective, single-dose pharmacokinetic clinical study was conducted in healthy volunteers at ASPET Journals on October 1, 2021 consisting of nine G143E heterozygotes, one G143E homozygote, and 12 non-carriers.

The AUC0–∞ ratio of oseltamivir carboxylate (active metabolite) to oseltamivir (parent molecule) was 23% lower in G143E heterozygotes compared to non-carriers (P =

0.006). The one G143E homozygous individual had an AUC0–∞ of oseltamivir that was approximately 360% greater than the non-carriers, indicating that loss of CES1 activity could profoundly impair oseltamivir activation (Tarkiainen et al., 2012).

4.1.5 G143E and dabigatran and sacubitril

Dabigatran and sacubitril are both prodrugs that need to be activated by CES1 in the liver (Figure 5). In vitro studies showed that the formation rates of the active metabolites of dabigatran and sacubitril were significantly lower in human livers carrying the G143E variant than in non-carrier samples (Shi et al., 2016c; Shi et al., 2016d). DMD Fast Forward. Published on December 23, 2019 as DOI: 10.1124/dmd.119.089680 This article has not been copyedited and formatted. The final version may differ from this version.

DMD # 89680 17

However, it remains undetermined whether the variant can affect the activation and therapeutic response of these two drugs in patients.

4.2 Pharmacogenetics of other CES1 genetic variants

In addition to G143E, many other CES1 variants have been studied for their effects on the PK and PD of CES1 substrate drugs. However, the results were generally inconclusive, and further studies are needed to determine the clinical significance of Downloaded from these variants.

4.2.1 E220G (rs200707504) dmd.aspetjournals.org

A nonsynonymous variant E220G, commonly referred to as c.662A>G, was suggested to decrease CES1 enzymatic activity in an in silico analysis (Oh et al., 2017).

In agreement with that prediction, an in vitro study on transfected cell lines found E220G at ASPET Journals on October 1, 2021 markedly decreased CES1 activity and the metabolisms of several CES1 substrates including enalapril, clopidogrel, and sacubitril (Wang et al., 2017a). Notably, E220G has a MAF of 0.55% in East Asians but is rare in other populations. To determine the clinical impact of E220G on the PK of a CES1 substrate, a single-dose oseltamivir (75mg) PK study was conducted in 20 healthy Korean volunteers. In this study, the variant was observed to have a marginal effect on the PK of oseltamivir and its active metabolite

(oseltamivir carboxylate); however, the differences were statistically insignificant. In the

E220G carriers (n=8), the AUC0-48h of oseltamivir was increased by 10% (P = 0.334), and the AUC0-48h of oseltamivir carboxylate was decreased by 5% (P = 0.513) relative to the non-carriers (n=12) (Oh et al., 2017).

4.2.2 S75N (rs2307240) DMD Fast Forward. Published on December 23, 2019 as DOI: 10.1124/dmd.119.089680 This article has not been copyedited and formatted. The final version may differ from this version.

DMD # 89680 18

S75N is one of the most common CES1 nonsynonymous SNP with MAFs ranging from 2% to 7% in different populations. A retrospective pharmacodynamics analysis was conducted to examine the effect of CES1 S75N on the outcome of clopidogrel therapy in patients with the coronary syndrome (n=851). The result showed that CES1 S75N carriers (n=372) had higher incidence of cerebrovascular events (P < 0.001), acute myocardial infarction (P < 0.001), and unstable angina (P < 0.001) compared to non- carriers. The study also found that the S75N polymorphism was more frequent in acute Downloaded from coronary syndrome patients (MAF 22%) than the general population (MAF 5%). The authors concluded that there was a significant association between the S75N dmd.aspetjournals.org polymorphism and the outcome of clopidogrel therapy (Xiao et al., 2017). However, this result conflicts with another study that found the S75N variant to be not associated with the outcomes of ADHD patients treated with methylphenidate (Johnson et al., 2013). at ASPET Journals on October 1, 2021 Furthermore, an in vitro study showed the S75N variant did not significantly alter the expression and activity of CES1 in transfected cells and human livers (Wang et al.,

2017a).

4.2.3 -816A>C (rs3785161)

The -816A>C polymorphism is located in the promoter region of CES1P1 VAR and has been suggested as a potential up-regulator of CES1P1 VAR expression

(Yoshimura et al., 2008). A prospective clinical study was conducted to examine the impact of -816A>C on the outcome of the ACEI prodrug (imidapril) therapy in hypertensive patients (n=105). The study found that after eight weeks of imidapril therapy, -816A>C homozygotes and heterozygotes (n=47) had greater systolic blood pressure reduction (24.1 mmHg) compared to non-carriers (17.6 mmHg) (P = 0.0184), DMD Fast Forward. Published on December 23, 2019 as DOI: 10.1124/dmd.119.089680 This article has not been copyedited and formatted. The final version may differ from this version.

DMD # 89680 19 indicating increased CES1 functionality in the carriers. The follow up in vitro study claimed that the -816A>C SNP may have enhanced transcription of the CES1P1 VAR gene (Geshi et al., 2005). The -816A>C SNP was also evaluated for its impact on the outcomes of dual antiplatelet therapy (i.e., aspirin and clopidogrel) in patients with coronary heart diseases (n=162). The -816A>C carriers (n=75) had decreased vasodilator stimulated phosphoprotein-platelet reactivity index (VASP-PRI) (P=0.014), indicating increased CES1 function in the carriers (Xie et al., 2014). Downloaded from

However, conflicting findings were reported by other studies. A study involving the outcome of clopidogrel treatment in patients undergoing percutaneous coronary dmd.aspetjournals.org intervention, -816A>C carriers showed a lower ADP-induced maximum platelet aggregation (21.5%, n=125) compared to non-carriers (31.7%, n=124) (P=0.001), indicating decreased CES1 function (Zou et al., 2014). Zhu et al. also performed a at ASPET Journals on October 1, 2021 retrospective pharmacogenetic analysis of the INternational VErapamil SR Trandolapril

(INVEST) study (n=486) and did not find an association between -816A>C and the blood pressure-lowering effect of trandolapril. The follow up in vitro study also showed -

816A>C genotype was not significantly associated with CES1 protein expression and trandolapril activation in human liver samples (n=100) (Zhu et al., 2016). Other researchers also noted that the CES1P1 VAR gene, which contains -816A>C, is considered functionally insignificant due to its low transcription efficiency (Tanimoto et al., 2007; Hosokawa et al., 2008).

4.2.4 -75G>T (rs3815583)

The -75G>T SNP is located in the promoter region of CES1 and was suspected to alter CES1 expression in the liver; however, the findings are conflicted. A study was DMD Fast Forward. Published on December 23, 2019 as DOI: 10.1124/dmd.119.089680 This article has not been copyedited and formatted. The final version may differ from this version.

DMD # 89680 20 performed to determine the association between the variant and appetite reduction (a side effect of methylphenidate) in children with ADHD (n=213). Appetite reduction was measured by the Barkley Stimulant Side Effect Rating Scale, and methylphenidate dose was titrated up for three months as tolerable. The carrier group (n=129) had worse appetite reduction compared to non-carriers (n=76) (41% vs 77%, P=0.01), indicating that the variant was associated with decreased CES1 function (Bruxel et al., 2013). A study in patients treated with irinotecan, however, showed a contrary finding, suggesting Downloaded from that the -75G>T variant confers greater CES1 function (Sai et al., 2010). CES1 is involved in the conversion of the prodrug irinotecan to its active metabolite, SN-38, and dmd.aspetjournals.org then further metabolized by UGT1As to inactive SN-38G. Following irinotecan treatment, patients who carried the T allele of this variant had higher plasma (SN-38 + SN-

38G)/irinotecan AUC ratios relative to non-carriers (P=0.027) following irinotecan at ASPET Journals on October 1, 2021 treatment (Sai et al., 2010).

Other CES1 substrates, isoniazid, and ACEI prodrugs were also studied in the context of -75G>T, however, no significant relationships were found between the variant and the medication responses. In one such study, the variant was evaluated for its effect on the outcomes of ACEI prodrugs in congestive heart failure patients (n=200) that underwent ACEI prodrug dose titrations. The study reported -75G>T did not significantly impact plasma ATII/ATI ratios, furthermore, the -75G>T variant was not significantly associated with fatal outcomes (i.e. cardiovascular death and all-cause death) (Nelveg-

Kristensen et al., 2016). The study with isoniazid had similar results showing no significant association between the variant and isoniazid-induced hepatotoxicity (n=170)

(Yamada et al., 2010).

4.2.5 1168-33C>A (rs2244613) DMD Fast Forward. Published on December 23, 2019 as DOI: 10.1124/dmd.119.089680 This article has not been copyedited and formatted. The final version may differ from this version.

DMD # 89680 21

Dabigatran (Pradaxa®) is a prodrug that needs to be activated by both CES1 and

CES2 to exert its anticoagulant effect (Figure 5). Paré and associates conducted a genome-wide association study of dabigatran in participants (n=2944) of the

Randomized Evaluation of Long-term Anticoagulation Therapy (RE-LY) clinical trial. The researchers concluded the CES1 intronic variant 1168-33C>A (rs2244613) is associated with lower trough concentrations of the active metabolite (15% decrease per allele; 95%

CI 10-19%) and a lower risk of any bleeding (odds ratio, 0.67; 95% CI 0.55-0.82) Downloaded from compared to non-carriers (Pare et al., 2013). However, an in vitro study did not find the variant to be associated with CES1 protein expression and dabigatran metabolism in dmd.aspetjournals.org human livers (Shi et al., 2016a). A prospective study also examined the impact of 1168-

33C>A in ADHD patients treated with methylphenidate, and found the variant to be associated with the occurrence of sadness, a side effect of short-acting at ASPET Journals on October 1, 2021 methylphenidate. However, researchers concluded this might be due to linkage disequilibrium with two SNPs of the noradrenaline transporter gene (Johnson et al.,

2013).

4.2.6 Copy number variation (i.e. CES1P1/CES1P1 VAR)

Many researchers have studied the impact of CNVs on CES1 functionality; however, the results are conflicted. Stage et al. found that participants with four functional copies of CES1 (n=5) had an increased AUC of d‐methylphenidate relative to the control group with two functional copies of CES1 (n=17) (61% increase, P = 0.011); participants with three copies of CES1 (n=2) had 45% increased AUC compared to the control group (P = 0.028) (Stage et al., 2017a). Stage et al. conducted a similar study with enalapril (n=43), however, they could not find a statistically significant correlation DMD Fast Forward. Published on December 23, 2019 as DOI: 10.1124/dmd.119.089680 This article has not been copyedited and formatted. The final version may differ from this version.

DMD # 89680 22 between CNV and enalapril PK (Stage et al., 2017b). When Sai et al. examined the effect of CNV on the irinotecan exposure, they found patients with multiple CES1 copies

(i.e., three or four) to have 1.24 fold higher irinotecan AUC relative to patients with two copies of CES1 (P = 0.0134) (Sai et al., 2010). Many researchers, however, did not find the relationship between CNVs and CES1 function. Suzaki et al. evaluated the relationship between CNVs of CES1 and oseltamivir PK parameters but did not find any correlation (Suzaki et al., 2013b). Nelveg-Kristensen et al. studied the relationship Downloaded from between CNV and ACEI prodrugs, and again, no association was found (Nelveg-

Kristensen et al., 2016). Moreover, an in vitro study showed CES1 protein expression dmd.aspetjournals.org levels to be comparable among human livers with different copy numbers of functional

CES1 gene (Wang et al., 2016a).

4.2.7 Other CES1 SNPs at ASPET Journals on October 1, 2021

In addition to the polymorphisms discussed above, sporadic reports have stated several CES1 SNPs to be associated with the outcomes of CES1 substrate medications.

For example, the SNP 1315 + 2025A>C (rs8192950) was associated with a decreased risk of ischemic events in patients (n=64) having symptomatic extracranial or intracranial stenosis and receiving dual antiplatelet therapy with clopidogrel for a minimum of five days (Zhao et al., 2016). Another retrospective sub-analysis of a capecitabine clinical study identified associations of 1168-41C>T (rs2244614), 690 + 129del (rs3217164),

95346T>C (rs7187684), -1232A>G (rs1186118) with severe early-onset of capecitabine- induced toxicity (Hamzic et al., 2017). None of these findings has been validated independently. DMD Fast Forward. Published on December 23, 2019 as DOI: 10.1124/dmd.119.089680 This article has not been copyedited and formatted. The final version may differ from this version.

DMD # 89680 23

A rare LOF variant, D260fs (c.70DelT), was reported in a clinical study (Zhu et al.,

2008). D260fs causes a deletion in exon 6, resulting in a frameshift and premature truncation. Moreover, an in vitro study with CES1 variants transfected cell lines examined the SNPs proximate to the CES1 , and identified four LOF nonsynonymous SNPs: G142E, G147C, Y170D, and R171C. However, these variants appear to be clinically insignificant due to their low MAFs (< 0.4%) (Wang et al., 2017b).

The above-mentioned CES1 SNPs and their impacts on the PK and PD of CES1 Downloaded from substrate medications are summarized in Table 3.

5. Non-genetic factors affecting CES1 expression and activity dmd.aspetjournals.org

5.1 Developmental expression of CES1

The developmental expression patterns of CES1 in human and mouse livers at ASPET Journals on October 1, 2021 were similar, and many in vitro studies have suggested that hepatic CES1 protein expression increases with age (Zhu et al., 2009a; Hines et al., 2016; Boberg et al.,

2017). An in vitro study with human liver samples (n=104) demonstrated the adult group

(≥ 18 years of age) to have had higher CES1 expression than children (0 days–10 years); meanwhile child group had higher CES1 expression than fetuses (82–224 gestation days). A follow-up study with liver microsomes showed that, in parallel with expression level, CES1 activity on hydrolyzing its substrate oseltamivir was also positively correlated with age (Yang et al., 2009). The same group did a similar in vitro human liver study with a slightly different age bracket, in which the liver samples were divided into five age groups: 1–31 days old (group 1), 35–70 days old (group 2), 89–119 days old (group 3), 123–198 days old (group 4), and over 18 years old (group 5).

Neonates (group 1) had 10% of the CES1 expression and hydrolysis levels compared to DMD Fast Forward. Published on December 23, 2019 as DOI: 10.1124/dmd.119.089680 This article has not been copyedited and formatted. The final version may differ from this version.

DMD # 89680 24 the adult group (group 5); pediatric groups (Group 2-4) had approximately 50% of the

CES1 expression and hydrolysis levels compared to an adult (Shi et al., 2011). Lastly, a similar in vitro study quantified CES1 protein levels in human liver samples of various ages (n=165). CES1 expression levels were 4.76 pmol/mg from birth to three weeks

(n=36); 15.8 pmol/mg for those aged three weeks to six years (n=90); and 16.6 pmol/mg for ages six years to 18 years (n=36). The study team concluded that the median CES1 expression level is directly correlated with age (P < 0.001) (Hines et al., 2016). Overall, Downloaded from

CES1 expression and activity levels are lower in neonates and pediatric patients; further studies are warranted to investigate the potential effect of CES1 maturation on the dmd.aspetjournals.org treatment outcome of CES1 substrate medications in patients in the early stages of development.

5.2 Sex difference of CES1 expression at ASPET Journals on October 1, 2021

Both in vitro and clinical studies have suggested that CES1 expression is higher in females than in males (Patrick et al., 2007; Zhu et al., 2009a; Shi et al., 2016d). A PK study on healthy volunteers revealed that males had significantly higher exposure to d- methylphenidate than females (Patrick et al., 2007). Nonetheless, females experienced a more pronounced stimulant effect, despite their lower exposure. Shi et al. observed significantly higher CES1 activity in female human liver samples (n = 56) compared to male samples (n = 46). A follow-up in vitro study with dabigatran suggested CES1 activity was higher in females than males (Shi et al., 2016d). However, such difference was not observed in another in vitro study using human liver samples (n=32) and mouse liver samples (n=9) (Zhu et al., 2009a). Further study is needed to examine the impact of sex on the CES1 expression level and the PK and PD of CES1 substrates. DMD Fast Forward. Published on December 23, 2019 as DOI: 10.1124/dmd.119.089680 This article has not been copyedited and formatted. The final version may differ from this version.

DMD # 89680 25

5.3 Drug-drug interactions

5.3.1 CES1 inhibitor - alcohol

To date, ethanol is the only known CES1 inhibitor that has been confirmed in multiple in vivo and in vitro studies. The impact of ethanol on the metabolism of the CES1 substrate, methylphenidate, was tested in healthy volunteers (n=14) (Zhu et al., 2017). D- methylphenidate comes as a single active ingredient (Focalin®) or in combination with l- Downloaded from methylphenidate (racemic mixture, Ritalin®). D-methylphenidate is approximately ten times more pharmacologically potent than l-methylphenidate, while l-methylphenidate is

a more efficient CES1 substrate (Figure 2). This clinical study used a pulsatile dosing dmd.aspetjournals.org regimen with methylphenidate (dl-methylphenidate 40 mg or d-methylphenidate 20 mg) and ethanol (0.6 g/kg, four hours after methylphenidate dose) to eliminate any potential

confounding effect of ethanol on methylphenidate absorption, as the methylphenidate at ASPET Journals on October 1, 2021 drug products (i.e., Ritalin® and Focalin®) might undergo faster gastric dissolution in the stomach if administered with alcohol. When alcohol and d-methylphenidate (Focalin®) were co-administered, the Cmax of d-methylphenidate was elevated by 27% (P = 0.001) and the AUC4→8h was elevated by 20% (P < 0.01); when alcohol and dl-methylphenidate

(Ritalin®) were co-administered, the Cmax of d-methylphenidate was elevated by 35% (P

< 0.01) and the AUC4→8h was elevated by 25% (P < 0.05) (Zhu et al., 2017). These results are consistent with the previous clinical trial by Patrick et al. In that study, when alcohol and d-methylphenidate (Focalin®) were co-administered, the d-methylphenidate

AUC was increased by 14%; when alcohol and dl-methylphenidate (Ritalin®) were co- administered, the d-methylphenidate AUC was increased by 21% (Patrick et al., 2013) .

Patrick and colleagues also showed that the co-administration of alcohol 30 min before DMD Fast Forward. Published on December 23, 2019 as DOI: 10.1124/dmd.119.089680 This article has not been copyedited and formatted. The final version may differ from this version.

DMD # 89680 26 or 30 min after methylphenidate had a similar impact on methylphenidate exposure

(Patrick et al., 2007). Both authors concluded that alcohol is a strong inhibitor of CES1, and the impact of CES1 inhibition is greater for dl-methylphenidate (Ritalin®) than for d- methylphenidate (Focalin®). Additionally, the DDI between methylphenidate and ethanol produced the transesterification metabolites d-ethylphenidate and l-ethylphenidate, and the plasma concentrations of l-ethylphenidate were much higher than d-ethylphenidate due to l-ethylphehidate being a more efficient CES1 substrate (Zhu et al., 2011; Zhu et Downloaded from al., 2017). Other in vivo studies with mice demonstrated similar results (Griffin et al.,

2010; Bell et al., 2011b; Griffin et al., 2013). dmd.aspetjournals.org The impact of alcohol on the CES1 function was also examined in the context of a different CES1 substrate, oseltamivir. A prospective health volunteer PK study (n=18) examined the interaction between oseltamivir 150mg (a recommended daily dose for the at ASPET Journals on October 1, 2021 treatment of influenza) and alcohol. Alcohol increased the oseltamivir AUC0-6h by 27%

(P=0.011) and decreased the AUC0-6h ratio of the active metabolite oseltamivir carboxylate to the parent compound oseltamivir by 34% (p<0.001) (Parker et al., 2015). However, co- administration of alcohol did not significantly affect the AUC0-24h of oseltamivir carboxylate.

These results are consistent with in silico analysis of the DDI between alcohol and oseltamivir

(Hu et al., 2014).

5.3.2 Other CES1 inhibitors: cannabis, protease inhibitors, aripiprazole, isradipine, tacrolimus, valproate

Besides alcohol, many drug products on the market have been suggested to be potent inhibitors of CES1 mainly by in vitro investigations (Table 4). A further clinical DMD Fast Forward. Published on December 23, 2019 as DOI: 10.1124/dmd.119.089680 This article has not been copyedited and formatted. The final version may differ from this version.

DMD # 89680 27 study with a validated CES1 substrate is needed to determine the clinical significance of these CES1 inhibitors.

An in vitro study with CES1 transfected cells suggested that cannabis (i.e., THC,

CBD, and CBN) can act as a potential CES1 inhibitor. The inhibition constant (Ki) values for THC, CBD and CBN were 0.541, 0.974, and 0.263 µM (0.170, 0.306, and 0.0817

µg/ml), respectively (Qian et al., 2019). This result could be clinically impactful as the use of cannabis is expected to increase in the next few years (Hasin, 2018). Downloaded from

Several protease inhibitors (i.e. nelfinavir, amprenavir, atazanavir, ritonavir, and

saquinavir) were identified as CES1 inhibitors by an in silico analysis and later confirmed dmd.aspetjournals.org by an in vitro incubation study. Among those, nelfinavir had a significantly higher inhibitory effect than the other agents. The relative CES1 activity towards PNPA (a

CES1 substrate) was 5.2%, 74.2%, 51.7%, 76.9%, and 67.8% of the control after at ASPET Journals on October 1, 2021 incubation with nelfinavir, ritonavir, amprenavir, saquinavir, and atazanavir, respectively

(Rhoades et al., 2012).

An in vitro study suggested aripiprazole, perphenazine, thioridazine, and fluoxetine to be potent inhibitors of CES1, and a complementary animal study (n=10) with FVB mice demonstrated that co-administration of aripiprazole and methylphenidate

(CES1 substrate) significantly increased the plasma concentrations of dl- methylphenidate (P < 0.01) (Zhu et al., 2010).

Moreover, a total of 27 cardiovascular, antiplatelet, anticoagulant and immunosuppressant drugs have been tested for CES1 inhibition using human liver microsomes and recombinant CES1. The results suggested isradipine (a dihydropyridine calcium antagonist, DHP) and tacrolimus (an immunosuppressive agent) DMD Fast Forward. Published on December 23, 2019 as DOI: 10.1124/dmd.119.089680 This article has not been copyedited and formatted. The final version may differ from this version.

DMD # 89680 28 to be potent CES1 inhibitors. CES1 activity towards PNPA was decreased to 17.6% with isradipine, and 28.4% with tacrolimus (Thomsen et al., 2014).

An In vitro study suggested valproate could inhibit CES1 function and affect rufinamide metabolism in both microsomes and cytosol. This result could be clinically significant as the two antiepileptic medications are often prescribed together when monotherapy is ineffective (Williams et al., 2011).

A combined ensemble docking and machine learning approach was utilized to Downloaded from identify potential CES1 inhibitors from 1114 FDA-approved drugs. Among the identified inhibitor candidates, four drugs including diltiazem, benztropine, iloprost, and treprostinil dmd.aspetjournals.org were found to inhibit CES1 activity in vitro with IC50 values ranging from 13.9 µM to

391.6 µM (Briand et al., 2019).

Lastly, an in vitro study suggested that some naturally occurring oxysterols and at ASPET Journals on October 1, 2021 fatty acids might significantly inhibit CES1 activity with IC50 values within the micromolar range (Crow et al., 2010). These compounds could potentially affect CES1-mediated detoxification and drug metabolism in vivo.

5.3.3 CES1 inducers

Overall, CES1 inducers are understudied relative to its inhibitors. Evidence suggests that various nuclear receptors might be involved in the regulation of CES1 expression (Staudinger et al., 2010). For example, several agonists of peroxisome proliferator-activated receptors (PPARs) induced the mRNA expressions of several

CES1 isoforms in mouse livers (Jones et al., 2013). A moderate increase of CES1 expression was observed in human hepatocytes treated with rifampicin, a prototypical human PXR-activating agent (Shi et al., 2008). An in vivo study with mice suggested DMD Fast Forward. Published on December 23, 2019 as DOI: 10.1124/dmd.119.089680 This article has not been copyedited and formatted. The final version may differ from this version.

DMD # 89680 29 that glucose could induce hepatic CES1 expression by stimulating CES1 promoter activity and increasing acetylation of histone 3 and histone 4 in the CES1 chromatin, indicating a potential role of CES1 in glucose homeostasis (Xu et al., 2014). Moreover, phenobarbital induced CES1 expression in mouse livers, and the inducibility was more prominent in neonatal mice relative to adult mice (Xiao et al., 2012). Again, a further clinical investigation is needed to determine the impacts of CES1 inducers on the PK and PD of CES1 substrate medications. Downloaded from

5.3.4 Drug-drug interactions between CES1 substrates

In addition to CES1 inhibitors and inducers, concomitant use of multiple CES1 substrate dmd.aspetjournals.org drugs can theoretically impact the substrate metabolism by competitively inhibiting the

CES1. This hypothesis has been tested in several studies. An in vitro study suggested

trandolapril and enalapril might increase clopidogrel activation (Kristensen et al., 2014). at ASPET Journals on October 1, 2021

Consistent with the in vitro study, a follow-up retrospective clinical study reported the concomitant use of ACEI prodrugs and clopidogrel increases the risk of clinically important bleeding in patients with myocardial infarction (n=70,934) (P=0.002). The clinical significance of this finding is, however, debatable as the hazard ratio of clinically significant bleeding for patients on concomitant therapy was 1.10 (95% CI 0.97-1.25)

(Kristensen et al., 2014). Another clinical study with the similar design did not report a significant association between the composite cardiovascular outcome and the concomitant use of ACEI prodrugs and clopidogrel in patients with myocardial infarction

(n=45,918). The adjusted odds ratio (aOR) for the perindopril was 0.94 (95% CI 0.76-

1.16), and for ramipril was 0.97 (95% CI 0.80-1.18), relative to lisinopril, an ACEI not metabolized by CES1 (Cressman et al., 2015). DMD Fast Forward. Published on December 23, 2019 as DOI: 10.1124/dmd.119.089680 This article has not been copyedited and formatted. The final version may differ from this version.

DMD # 89680 30

6. Disease states related to CES1

A prospective clinical study was conducted in monozygotic and dizygotic twin subjects (62–83 years) with (n = 48) or without (n = 247) type 2 diabetes mellitus (T2DM) to examine the association of CES1 with adiposity and metabolic function. CES1 mRNA expression level in adipose tissue was positively associated with body-mass index

(P<0.001), fasting glucose level (P = 0.002), insulin (P = 0.006), and triglycerides

(P = 0.003) (Friedrichsen et al., 2013). Recent studies have also found that CES1 Downloaded from function was positively correlated with increased liver lipid storage and plasma lipid concentrations, indicating that CES1 might be heavily involved in lipid metabolism and is dmd.aspetjournals.org a potential drug target for the treatment of human metabolic disorders (Kaddurah-Daouk et al., 2018; Lian et al., 2018a; Lian et al., 2018b).

at ASPET Journals on October 1, 2021 7. Conclusion and future directions

In sum, G143E (rs71647871) is the only clinically significant LOF CES1 variant identified to date, and alcohol is the only potent CES1 inhibitor that significantly affect

CES1-mediated drug metabolism both in vivo and in vitro. However, G143E (MAF 2–

4%, carrier frequency 4-8%) and alcohol-induced DDI can only be able to explain a small portion of the interindividual variability in the CES1 function. Previous in vitro studies have demonstrated marked variability of CES1 activity and expression in human liver samples not carrying G143E (Shi et al., 2016b; Wang et al., 2016a). In fact, analysis of the correlation between CES1 expression and activity revealed that the majority of interindividual variability in the CES1 function is due to variation in CES1 protein expression (Wang et al., 2016a). DMD Fast Forward. Published on December 23, 2019 as DOI: 10.1124/dmd.119.089680 This article has not been copyedited and formatted. The final version may differ from this version.

DMD # 89680 31

Unfortunately, the mechanism by which CES1 protein expression is regulated remains largely unexplored. Notably, most of the existing gene expression regulation studies were based upon the measurement of mRNA expression levels. However, increasing evidence suggests that mRNA expression correlates poorly with protein expression for many including CES1 and most DMEs, which could result in false identification of gene expression regulators (Ohtsuki et al., 2012). Recent advances in

LC-MS/MS-based proteomics have allowed for accurate CES1 protein quantification. Downloaded from

The application of CES1 proteomics in a large set of clinical samples (e.g. human livers) is expected to uncover important factors influencing CES1 expression, such as genetic dmd.aspetjournals.org polymorphisms, disease conditions, inducers, and post-transcriptional modification

(Wang et al., 2016b; He et al., 2019); the findings from such research will lead to the development of an individualized pharmacotherapy approach for improving the efficacy at ASPET Journals on October 1, 2021 and safety of many medications metabolized by CES1.

DMD Fast Forward. Published on December 23, 2019 as DOI: 10.1124/dmd.119.089680 This article has not been copyedited and formatted. The final version may differ from this version.

DMD # 89680 32

Acknowledgements

This work was supported in part by the National Heart, Lung, and Blood Institute [Grant

R01HL126969].

Authorship Contributions

Wrote or contributed to the writing of the manuscript: Lucy Her, Hao-Jie Zhu

Downloaded from

Conflict of interest

The authors state no conflict of interests dmd.aspetjournals.org

at ASPET Journals on October 1, 2021 DMD Fast Forward. Published on December 23, 2019 as DOI: 10.1124/dmd.119.089680 This article has not been copyedited and formatted. The final version may differ from this version.

DMD # 89680 33

Reference

Ando Y, Saka H, Ando M, Sawa T, Muro K, Ueoka H, Yokoyama A, Saitoh S, Shimokata K, and Hasegawa Y (2000) Polymorphisms of UDP-glucuronosyltransferase gene and irinotecan toxicity: a pharmacogenetic analysis. Cancer Res 60:6921-6926. Arena de Souza V, Scott DJ, Nettleship JE, Rahman N, Charlton MH, Walsh MA, and Owens RJ (2015) Comparison of the Structure and Activity of Glycosylated and Aglycosylated Human Carboxylesterase 1. PLoS One 10:e0143919. Bell GH, Griffin WC, 3rd, and Patrick KS (2011a) Oral and transdermal DL-methylphenidate-ethanol interactions in C57BL/6J mice: potentiation of locomotor activity with oral delivery. Pharmacol Biochem Behav 100:264-270. Bell GH, Novak AJ, Griffin WC, and Patrick KS (2011b) Transdermal and Oral dl-Methylphenidate-Ethanol Interactions in C57BL/6J Mice: Transesterification to Ethylphenidate and Elevation of d- Methylphenidate Concentrations. J Pharm Sci-Us 100:2966-2978. Downloaded from Boberg M, Vrana M, Mehrotra A, Pearce RE, Gaedigk A, Bhatt DK, Leeder JS, and Prasad B (2017) Age- Dependent Absolute Abundance of Hepatic (CES1 and CES2) by LC-MS/MS Proteomics: Application to PBPK Modeling of Oseltamivir In Vivo Pharmacokinetics in Infants. Drug Metab Dispos 45:216-223. Bozzi LM, Mitchell BD, Lewis JP, Ryan KA, Herzog WR, O'Connell JR, Horenstein RB, Shuldiner AR, and dmd.aspetjournals.org Yerges-Armstrong LM (2016) The Pharmacogenomics of Anti-Platelet Intervention (PAPI) Study: Variation in Platelet Response to Clopidogrel and Aspirin. Curr Vasc Pharmacol 14:116-124. Briand E, Thomsen R, Linnet K, Rasmussen HB, Brunak S, and Taboureau O (2019) Combined Ensemble Docking and Machine Learning in Identification of Therapeutic Agents with Potential Inhibitory Effect on Human CES1. Molecules 24. Bruxel EM, Salatino-Oliveira A, Genro JP, Zeni CP, Polanczyk GV, Chazan R, Rohde LA, and Hutz MH (2013) at ASPET Journals on October 1, 2021 Association of a carboxylesterase 1 polymorphism with appetite reduction in children and adolescents with attention-deficit/hyperactivity disorder treated with methylphenidate. Pharmacogenomics J 13:476-480. Chaturvedi S (2004) The Seventh Report of the Joint National Committee on Prevention, Detection, Evaluation, and Treatment of High Blood Pressure (JNC 7): is it really practical? Natl Med J India 17:227. Cressman AM, Macdonald EM, Fernandes KA, Gomes T, Paterson JM, Mamdani MM, Juurlink DN, and Effectiveness CDS (2015) A population-based study of the drug interaction between clopidogrel and angiotensin converting enzyme inhibitors. Brit J Clin Pharmaco 80:662-669. Crow JA, Herring KL, Xie S, Borazjani A, Potter PM, and Ross MK (2010) Inhibition of carboxylesterase activity of THP1 monocytes/macrophages and recombinant human carboxylesterase 1 by oxysterols and fatty acids. Biochim Biophys Acta 1801:31-41. Dahlgren FS, Shay DK, Izurieta HS, Forshee RA, Wernecke M, Chillarige Y, Lu Y, Kelman JA, and Reed C (2018) Evaluating oseltamivir prescriptions in Centers for Medicare and Medicaid Services medical claims records as an indicator of seasonal influenza in the United States. Influenza Other Resp 12:465-474. DiPiro JT (2017) Pharmacotherapy : a pathophysiologic approach. McGraw-Hill Education, New York. Dolinsky VW, Gilham D, Alam M, Vance DE, and Lehner R (2004) Triacylglycerol hydrolase: role in intracellular lipid metabolism. Cell Mol Life Sci 61:1633-1651. Foti RS and Dalvie DK (2016) Cytochrome P450 and Non-Cytochrome P450 Oxidative Metabolism: Contributions to the Pharmacokinetics, Safety, and Efficacy of Xenobiotics. Drug Metab Dispos 44:1229-1245. DMD Fast Forward. Published on December 23, 2019 as DOI: 10.1124/dmd.119.089680 This article has not been copyedited and formatted. The final version may differ from this version.

DMD # 89680 34

Foye WO, Lemke TL, and Williams DA (2013) Foye's principles of medicinal chemistry. Wolters Kluwer Health/Lippincott Williams & Wilkins, Philadelphia. Friedrichsen M, Poulsen P, Wojtaszewski J, Hansen PR, Vaag A, and Rasmussen HB (2013) Carboxylesterase 1 gene duplication and mRNA expression in adipose tissue are linked to obesity and metabolic function. PLoS One 8:e56861. Fukami T, Nakajima M, Maruichi T, Takahashi S, Takamiya M, Aoki Y, McLeod HL, and Yokoi T (2008) Structure and characterization of human carboxylesterase 1A1, 1A2, and 1A3 genes. Pharmacogenet Genomics 18:911-920. Fukami T and Yokoi T (2012) The emerging role of human esterases. Drug Metab Pharmacokinet 27:466- 477. Geshi E, Kimura T, Yoshimura M, Suzuki H, Koba S, Sakai T, Saito T, Koga A, Muramatsu M, and Katagiri T (2005) A single nucleotide polymorphism in the carboxylesterase gene is associated with the responsiveness to imidapril medication and the promoter activity. Hypertens Res 28:719-725. Ghosh S, Zhao B, Bie J, and Song J (2010) Macrophage cholesteryl ester mobilization and atherosclerosis. Downloaded from Vascul Pharmacol 52:1-10. Griffin WC, 3rd, McGovern RW, Bell GH, Randall PK, Middaugh LD, and Patrick KS (2013) Interactive effects of methylphenidate and alcohol on discrimination, conditioned place preference and motor coordination in C57BL/6J mice. Psychopharmacology (Berl) 225:613-625. Griffin WC, Novak AJ, Middaugh LD, and Patrick KS (2010) The interactive effects of methylphenidate and dmd.aspetjournals.org ethanol on ethanol consumption and locomotor activity in mice. Pharmacol Biochem Be 95:267- 272. Hamzic S, Kummer D, Milesi S, Mueller D, Joerger M, Aebi S, Amstutz U, and Largiader CR (2017) Novel Genetic Variants in Carboxylesterase 1 Predict Severe Early-Onset Capecitabine-Related Toxicity. Clin Pharmacol Ther 102:796-804. Hasin DS (2018) US Epidemiology of Cannabis Use and Associated Problems. Neuropsychopharmacol at ASPET Journals on October 1, 2021 43:195-212. He B, Shi J, Wang XW, Jiang H, and Zhu HJ (2019) Label-free absolute protein quantification with data- independent acquisition. J Proteomics 200:51-59. Heal DJ and Pierce DM (2006) Methylphenidate and its isomers - Their role in the treatment of attention- deficit hyperactivity disorder using a transdermal delivery system. Cns Drugs 20:713-738. Hines RN, Simpson PM, and McCarver DG (2016) Age-Dependent Human Hepatic Carboxylesterase 1 (CES1) and Carboxylesterase 2 (CES2) Postnatal Ontogeny. Drug Metab Dispos 44:959-966. Hosokawa M, Furihata T, Yaginuma Y, Yamamoto N, Watanabe N, Tsukada E, Ohhata Y, Kobayashi K, Satoh T, and Chiba K (2008) Structural organization and characterization of the regulatory element of the human carboxylesterase (CES1A1 and CES1A2) genes. Drug Metab Pharmacok 23:73-84. Hu ZY, Edginton AN, Laizure SC, and Parker RB (2014) Physiologically Based Pharmacokinetic Modeling of Impaired Carboxylesterase-1 Activity: Effects on Oseltamivir Disposition (vol 53, pg 825, 2014). Clinical Pharmacokinetics 53:959-959. Imai T, Taketani M, Shii M, Hosokawa M, and Chiba K (2006) Substrate specificity of carboxylesterase isozymes and their contribution to hydrolase activity in human liver and small intestine. Drug Metab Dispos 34:1734-1741. Innocenti F, Undevia SD, Iyer L, Chen PX, Das S, Kocherginsky M, Karrison T, Janisch L, Ramirez J, Rudin CM, Vokes EE, and Ratain MJ (2004) Genetic variants in the UDP-glucuronosyltransferase 1A1 gene predict the risk of severe neutropenia of irinotecan. J Clin Oncol 22:1382-1388. Jewell C, Bennett P, Mutch E, Ackermann C, and Williams FM (2007) Inter-individual variability in esterases in human liver. Biochem Pharmacol 74:932-939. DMD Fast Forward. Published on December 23, 2019 as DOI: 10.1124/dmd.119.089680 This article has not been copyedited and formatted. The final version may differ from this version.

DMD # 89680 35

Jiang XL, Samant S, Lewis JP, Horenstein RB, Shuldiner AR, Yerges-Armstrong LM, Peletier LA, Lesko LJ, and Schmidt S (2016) Development of a physiology-directed population pharmacokinetic and pharmacodynamic model for characterizing the impact of genetic and demographic factors on clopidogrel response in healthy adults. Eur J Pharm Sci 82:64-78. Johnson KA, Barry E, Lambert D, Fitzgerald M, McNicholas F, Kirley A, Gill M, Bellgrove MA, and Hawi Z (2013) Methylphenidate Side Effect Profile Is Influenced by Genetic Variation in the Attention- Deficit/Hyperactivity Disorder-Associated CES1 Gene. J Child Adol Psychop 23:655-664. Jones RD, Taylor AM, Tong EY, and Repa JJ (2013) Carboxylesterases Are Uniquely Expressed among Tissues and Regulated by Nuclear Hormone Receptors in the Mouse. Drug Metab Dispos 41:40- 49. Kaddurah-Daouk R, Hankemeier T, Scholl EH, Baillie R, Harms A, Stage C, Dalhoff KP, Jurgens G, Taboureau O, Nzabonimpa GS, Motsinger-Reif AA, Thomsen R, Linnet K, Rasmussen HB, Consortium I, and Network PR (2018) Pharmacometabolomics Informs About Pharmacokinetic Profile of Methylphenidate. Cpt-Pharmacomet Syst 7:525-533. Downloaded from Kristensen KE, Zhu HJ, Wang X, Gislason GH, Torp-Pedersen C, Rasmussen HB, Markowitz JS, and Hansen PR (2014) Clopidogrel Bioactivation and Risk of Bleeding in Patients Cotreated With Angiotensin- Converting Enzyme Inhibitors After Myocardial Infarction: A Proof-of-Concept Study. Clin Pharmacol Ther 96:713-722. Laizure SC, Herring V, Hu ZY, Witbrodt K, and Parker RB (2013) The Role of Human Carboxylesterases in dmd.aspetjournals.org Drug Metabolism: Have We Overlooked Their Importance? Pharmacotherapy 33:210-222. Laizure SC, Parker RB, Herring VL, and Hu ZY (2014) Identification of carboxylesterase-dependent dabigatran etexilate hydrolysis. Drug Metab Dispos 42:201-206. Lewis JP, Horenstein RB, Ryan K, O'Connell JR, Gibson Q, Mitchell BD, Tanner K, Chai S, Bliden KP, Tantry US, Peer CJ, Figg WD, Spencer SD, Pacanowski MA, Gurbel PA, and Shuldiner AR (2013) The functional G143E variant of carboxylesterase 1 is associated with increased clopidogrel active at ASPET Journals on October 1, 2021 metabolite levels and greater clopidogrel response. Pharmacogenet Genom 23:1-8. Lian JH, Bahitham W, Panigrahi R, Nelson R, Li LN, Watts R, Thiesen A, Lemieux MJ, and Lehner R (2018a) Genetic variation in human carboxylesterase CES1 confers resistance to hepatic steatosis. Bba- Mol Cell Biol L 1863:688-699. Lian JH, Nelson R, and Lehner R (2018b) Carboxylesterases in lipid metabolism: from mouse to human. Protein & Cell 9:178-195. Mahmoudpour SH, Asselbergs FW, de Keyser CE, Souverein PC, Hofman A, Stricker BH, de Boer A, and Maitland-van der Zee AH (2015) Change in prescription pattern as a potential marker for adverse drug reactions of angiotensin converting enzyme inhibitors. Int J Clin Pharm-Net 37:1095-1103. Mehta SR, Yusuf S, Peters RJ, Bertrand ME, Lewis BS, Natarajan MK, Malmberg K, Rupprecht H, Zhao F, Chrolavicius S, Copland I, Fox KA, and Clopidogrel in Unstable angina to prevent Recurrent Events trial I (2001) Effects of pretreatment with clopidogrel and aspirin followed by long-term therapy in patients undergoing percutaneous coronary intervention: the PCI-CURE study. Lancet 358:527- 533. Nelveg-Kristensen KE, Madsen MB, Torp-Pedersen C, Kober L, Egfjord M, Hansen T, Pedersen O, Rasmussen HB, and Hansen PR (2016) Prognostic impact of carboxylesterase 1 gene variants in patients with congestive heart failure treated with angiotensin-converting enzyme inhibitors. Pharmacogenet Genomics 26:169-177. Nemoda Z, Angyal N, Tarnok Z, Gadoros J, and Sasvari-Szekely M (2009) Carboxylesterase 1 gene polymorphism and methylphenidate response in ADHD. Neuropharmacology 57:731-733. Neos Therapeutics IGP, TX (2017) Product Information: COTEMPLA XR-ODT(TM) oral extended-release disintegrating tablets, methylphenidate oral extended-release disintegrating tablets. DMD Fast Forward. Published on December 23, 2019 as DOI: 10.1124/dmd.119.089680 This article has not been copyedited and formatted. The final version may differ from this version.

DMD # 89680 36

Oh J, Lee S, Lee H, Cho JY, Yoon SH, Jang IJ, Yu KS, and Lim KS (2017) The novel carboxylesterase 1 variant c.662A>G may decrease the bioactivation of oseltamivir in humans. Plos One 12. Ohtsuki S, Schaefer O, Kawakami H, Inoue T, Liehner S, Saito A, Ishiguro N, Kishimoto W, Ludwig- Schwellinger E, Ebner T, and Terasaki T (2012) Simultaneous Absolute Protein Quantification of Transporters, Cytochromes P450, and UDP-Glucuronosyltransferases as a Novel Approach for the Characterization of Individual Human Liver: Comparison with mRNA Levels and Activities. Drug Metab Dispos 40:83-92. Pare G, Eriksson N, Lehr T, Connolly S, Eikelboom J, Ezekowitz MD, Axelsson T, Haertter S, Oldgren J, Reilly P, Siegbahn A, Syvanen AC, Wadelius C, Wadelius M, Zimdahl-Gelling H, Yusuf S, and Wallentin L (2013) Genetic Determinants of Dabigatran Plasma Levels and Their Relation to Bleeding. Circulation 127:1404-+. Parker RB, Hu ZY, Meibohm B, and Laizure SC (2015) Effects of alcohol on human carboxylesterase drug metabolism. Clin Pharmacokinet 54:627-638. Patrick KS, Straughn AB, Minhinnett RR, Yeatts SD, Herrin AE, DeVane CL, Malcolm R, Janis GC, and Downloaded from Markowitz JS (2007) Influence of ethanol and gender on methylphenidate pharmacokinetics and pharmacodynamics. Clin Pharmacol Ther 81:346-353. Patrick KS, Straughn AB, Reeves OT, Bernstein H, Bell GH, Anderson ER, and Malcolm RJ (2013) Differential Influences of Ethanol on Early Exposure to Racemic Methylphenidate Compared with Dexmethylphenidate in Humans. Drug Metab Dispos 41:197-205. dmd.aspetjournals.org Qian YL, Wang XW, and Markowitz JS (2019) In Vitro Inhibition of Carboxylesterase 1 by Major Cannabinoids and Selected Metabolites. Drug Metab Dispos 47:465-472. Rautio J, Karkkainen J, and Sloan KB (2017) Prodrugs - Recent approvals and a glimpse of the pipeline. Eur J Pharm Sci 109:146-161. Rautio J, Kumpulainen H, Heimbach T, Oliyai R, Oh D, Jarvinen T, and Savolainen J (2008) Prodrugs: design and clinical applications. Nat Rev Drug Discov 7:255-270. at ASPET Journals on October 1, 2021 Rhoades JA, Peterson YK, Zhu HJ, Appel DI, Peloquin CA, and Markowitz JS (2012) Prediction and In Vitro Evaluation of Selected Protease Inhibitor Antiviral Drugs as Inhibitors of Carboxylesterase 1: A Potential Source of Drug-Drug Interactions. Pharm Res-Dordr 29:972-982. Roe MT, Armstrong PW, Fox KAA, White HD, Prabhakaran D, Goodman SG, Cornel JH, Bhatt DL, Clemmensen P, Martinez F, Ardissino D, Nicolau JC, Boden WE, Gurbel PA, Ruzyllo W, Dalby AJ, McGuire DK, Leiva-Pons JL, Parkhomenko A, Gottlieb S, Topacio GO, Hamm C, Pavlides G, Goudev AR, Oto A, Tseng CD, Merkely B, Gasparovic V, Corbalan R, Cinteza M, McLendon RC, Winters KJ, Brown EB, Lokhnygina Y, Aylward PE, Huber K, Hochman JS, Ohman EM, and Investigators TA (2012) Prasugrel versus Clopidogrel for Acute Coronary Syndromes without Revascularization. New Engl J Med 367:1297-1309. Ross MK, Streit TM, and Herring KL (2010) Carboxylesterases: Dual roles in lipid and pesticide metabolism. J Pestic Sci 35:257-264. Sai K, Saito Y, Tatewaki N, Hosokawa M, Kaniwa N, Nishimaki-Mogami T, Naito M, Sawada J, Shirao K, Hamaguchi T, Yamamoto N, Kunitoh H, Tamura T, Yamada Y, Ohe Y, Yoshida T, Minami H, Ohtsu A, Matsumura Y, Saijo N, and Okuda H (2010) Association of carboxylesterase 1A genotypes with irinotecan pharmacokinetics in Japanese cancer patients. Br J Clin Pharmacol 70:222-233. Satoh T and Hosokawa M (2006) Structure, function and regulation of carboxylesterases. Chem Biol Interact 162:195-211. Schubert I, Koster I, and Lehmkuhl G (2010) The Changing Prevalence of Attention-Deficit/Hyperactivity Disorder and Methylphenidate Prescriptions A Study of Data From a Random Sample of Insurees of the AOK Health Insurance Company in the German State of Hesse, 2000-2007. Dtsch Arztebl Int 107:615-U622. DMD Fast Forward. Published on December 23, 2019 as DOI: 10.1124/dmd.119.089680 This article has not been copyedited and formatted. The final version may differ from this version.

DMD # 89680 37

Scott SA, Sangkuhl K, Stein CM, Hulot JS, Mega JL, Roden DM, Klein TE, Sabatine MS, Johnson JA, and Shuldiner AR (2013) Clinical Pharmacogenetics Implementation Consortium Guidelines for CYP2C19 Genotype and Clopidogrel Therapy: 2013 Update. Clin Pharmacol Ther 94:317-323. Shah VP and Amidon GL (2014) G.L. Amidon, H. Lennernas, V.P. Shah, and J.R. Crison. A theoretical basis for a biopharmaceutic drug classification: the correlation of in vitro drug product dissolution and in vivo bioavailability, Pharm Res 12, 413-420, 1995--backstory of BCS. Aaps J 16:894-898. Shi D, Yang D, Prinssen EP, Davies BE, and Yan B (2011) Surge in expression of carboxylesterase 1 during the post-neonatal stage enables a rapid gain of the capacity to activate the anti-influenza prodrug oseltamivir. J Infect Dis 203:937-942. Shi D, Yang J, Yang DF, and Yan B (2008) Dexamethasone suppresses the expression of multiple rat carboxylesterases through transcriptional repression: Evidence for an involvement of the glucocorticoid receptor. Toxicology 254:97-105. Shi DS, Yang J, Yang DF, LeCluyse EL, Black C, You L, Akhlaghi F, and Yan BF (2006) Anti-influenza prodrug oseltamivir is activated by carboxylesterase human carboxylesterase 1, and the activation is Downloaded from inhibited by antiplatelet agent clopidogrel. J Pharmacol Exp Ther 319:1477-1484. Shi J, Wang X, Nguyen JH, Bleske BE, Liang Y, Liu L, and Zhu HJ (2016a) Dabigatran etexilate activation is affected by the CES1 genetic polymorphism G143E (rs71647871) and gender. Biochem Pharmacol 119:76-84. Shi J, Wang XW, Eyler RF, Liang Y, Liu L, Mueller BA, and Zhu HJ (2016b) Association of Oseltamivir dmd.aspetjournals.org Activation with Gender and Carboxylesterase 1 Genetic Polymorphisms. Basic Clin Pharmacol 119:555-561. Shi J, Wang XW, Nguyen J, Wu A, Bleske B, and Zhu HJ (2016c) Sacubitril is selectively activated by carboxylesterase 1 (CES1) in the liver and the activation is affected by CES1 genetic variation. Faseb J 30. Shi J, Wang XW, Nguyen JH, Bleske BE, Liang Y, Liu L, and Zhu HJ (2016d) Dabigatran etexilate activation at ASPET Journals on October 1, 2021 is affected by the CES1 genetic polymorphism G143E (rs71647871) and gender. Biochemical Pharmacology 119:76-84. Shmyr D, Van der Merwe V, Yakiwchuk E, Barry A, and Kosar L (2017) Triple antithrombotic therapy for atrial fibrillation and coronary stents. Can Fam Physician 63:375-381. Singh S, Barghoorn J, Bagdonas A, Adler J, Treanor J, Kinnersley N, and Ward P (2003) Clinical benefits with oseltamivir in treating influenza in adult populations : results of a pooled and subgroup analysis. Clin Drug Investig 23:561-569. Stage C, Jurgens G, Guski LS, Thomsen R, Bjerre D, Ferrero-Miliani L, Lyauk YK, Rasmussen HB, Dalhoff K, and Consortium I (2017a) The impact of CES1 genotypes on the pharmacokinetics of methylphenidate in healthy Danish subjects. Br J Clin Pharmacol 83:1506-1514. Stage C, Jurgens G, Guski LS, Thomsen R, Bjerre D, Ferrero-Miliani L, Lyauk YK, Rasmussen HB, Dalhoff K, and Consortium I (2017b) The Pharmacokinetics of Enalapril in Relation to CES1 Genotype in Healthy Danish Volunteers. Basic Clin Pharmacol 121:487-492. Staudinger JL, Xu C, Cui YJ, and Klaassen CD (2010) Nuclear receptor-mediated regulation of carboxylesterase expression and activity. Expert Opin Drug Metab Toxicol 6:261-271. Steinhubl SR, Berger PB, Mann JT, 3rd, Fry ET, DeLago A, Wilmer C, Topol EJ, and Observation CICftRoED (2002) Early and sustained dual oral antiplatelet therapy following percutaneous coronary intervention: a randomized controlled trial. JAMA 288:2411-2420. Suzaki Y, Uemura N, Hosokawa M, and Ohashi K (2013a) Gly143Glu polymorphism of the human carboxylesterase1 gene in an Asian population. Eur J Clin Pharmacol 69:735-736. Suzaki Y, Uemura N, Takada M, Ohyama T, Itohda A, Morimoto T, Imai H, Hamasaki H, Inano A, Hosokawa M, Tateishi M, and Ohashi K (2013b) The effect of carboxylesterase 1 (CES1) DMD Fast Forward. Published on December 23, 2019 as DOI: 10.1124/dmd.119.089680 This article has not been copyedited and formatted. The final version may differ from this version.

DMD # 89680 38

polymorphisms on the pharmacokinetics of oseltamivir in humans. Eur J Clin Pharmacol 69:21- 30. Tanimoto K, Kaneyasu M, Shimokuni T, Hiyama K, and Nishiyama M (2007) Human carboxylesterase 1A2 expressed from carboxylesterase 1A1 and 1A2 genes is a potent predictor of CPT-11 cytotoxicity in vitro. Pharmacogenet Genom 17:1-10. Tarkiainen EK, Backman JT, Neuvonen M, Neuvonen PJ, Schwab M, and Niemi M (2012) Carboxylesterase 1 polymorphism impairs oseltamivir bioactivation in humans. Clin Pharmacol Ther 92:68-71. Tarkiainen EK, Holmberg MT, Tornio A, Neuvonen M, Neuvonen PJ, Backman JT, and Niemi M (2015a) Carboxylesterase 1 c.428G>A single nucleotide variation increases the antiplatelet effects of clopidogrel by reducing its hydrolysis in humans. Clin Pharmacol Ther 97:650-658. Tarkiainen EK, Tornio A, Holmberg MT, Launiainen T, Neuvonen PJ, Backman JT, and Niemi M (2015b) Effect of carboxylesterase 1 c.428G > A single nucleotide variation on the pharmacokinetics of quinapril and enalapril. Br J Clin Pharmacol 80:1131-1138. Thomsen R, Rasmussen HB, and Linnet K (2014) In Vitro Drug Metabolism by Human Carboxylesterase 1 Downloaded from with Focus on Angiotensin-Converting Enzyme Inhibitors. Drug Metab Rev 45:192-193. Wallentin L, Becker RC, Budaj A, Cannon CP, Emanuelsson H, Held C, Horrow J, Husted S, James S, Katus H, Mahaffey KW, Scirica BM, Skene A, Steg PG, Storey RF, Harrington RA, Investigators P, Freij A, and Thorsen M (2009) Ticagrelor versus clopidogrel in patients with acute coronary syndromes. N Engl J Med 361:1045-1057. dmd.aspetjournals.org Wang X, Rida N, Shi J, Wu A, Bleske B, and Zhu HJ (2017a) A comprehensive functional assessment of carboxylesterase 1 nonsynonymous polymorphisms. Drug Metab Dispos. Wang X, Rida N, Shi J, Wu AH, Bleske BE, and Zhu HJ (2017b) A Comprehensive Functional Assessment of Carboxylesterase 1 Nonsynonymous Polymorphisms. Drug Metab Dispos 45:1149-1155. Wang X, Wang G, Shi J, Aa J, Comas R, Liang Y, and Zhu HJ (2016a) CES1 genetic variation affects the activation of angiotensin-converting enzyme inhibitors. Pharmacogenomics J 16:220-230. at ASPET Journals on October 1, 2021 Wang XW, Liang Y, Liu L, Shi J, and Zhu HJ (2016b) Targeted absolute quantitative proteomics with SILAC internal standards and unlabeled full-length protein calibrators (TAQSI). Rapid Commun Mass Sp 30:553-561. Williams ET, Carlson JE, Lai WG, Wong YN, Yoshimura T, Critchley DJ, and Narurkar M (2011) Investigation of the metabolism of rufinamide and its interaction with valproate. Drug Metab Lett 5:280-289. Wiviott SD, Braunwald E, McCabe CH, Montalescot G, Ruzyllo W, Gottlieb S, Neumann FJ, Ardissino D, De Servi S, Murphy SA, Riesmeyer J, Weerakkody G, Gibson CM, Antman EM, and Investigators T-T (2007) Prasugrel versus clopidogrel in patients with acute coronary syndromes. N Engl J Med 357:2001-2015. Xiao D, Chen YT, Yang D, and Yan B (2012) Age-related inducibility of carboxylesterases by the antiepileptic agent phenobarbital and implications in drug metabolism and lipid accumulation. Biochem Pharmacol 84:232-239. Xiao FY, Luo JQ, Liu M, Chen BL, Cao S, Liu ZQ, Zhou HH, Zhou G, and Zhang W (2017) Effect of carboxylesterase 1 S75N on clopidogrel therapy among acute coronary syndrome patients. Sci Rep-Uk 7. Xie C, Ding X, Gao J, Wang H, Hang Y, Zhang H, Zhang J, Jiang B, and Miao L (2014) The effects of CES1A2 A(-816)C and CYP2C19 loss-of-function polymorphisms on clopidogrel response variability among Chinese patients with coronary heart disease. Pharmacogenet Genomics 24:204-210. Xu JS, Yin LY, Xu Y, Li YY, Zalzala M, Cheng G, and Zhang YQ (2014) Hepatic Carboxylesterase 1 Is Induced by Glucose and Regulates Postprandial Glucose Levels. Plos One 9. DMD Fast Forward. Published on December 23, 2019 as DOI: 10.1124/dmd.119.089680 This article has not been copyedited and formatted. The final version may differ from this version.

DMD # 89680 39

Yamada S, Richardson K, Tang M, Halaschek-Wiener J, Cook VJ, Fitzgerald JM, Elwood K, Marra F, and Brooks-Wilson A (2010) Genetic variation in carboxylesterase genes and susceptibility to isoniazid-induced hepatotoxicity. Pharmacogenomics J 10:524-536. Yancy CW, Jessup M, Bozkurt B, Butler J, Casey DE, Colvin MM, Drazner MH, Filippatos GS, Fonarow GC, Givertz MM, Hollenberg SM, Lindenfeld J, Masoudi FA, McBride PE, Peterson PN, Stevenson LW, and Westlake C (2017) 2017 ACC/AHA/HFSA Focused Update of the 2013 ACCF/AHA Guideline for the Management of Heart Failure A Report of the American College of Cardiology/American Heart Association Task Force on Clinical Practice Guidelines and the Heart Failure Society of America. Circulation 136:E137-+. Yang DF, Pearce RE, Wang XL, Gaedigk R, Wan YJY, and Yan BF (2009) Human carboxylesterases HCE1 and HCE2: Ontogenic expression, inter-individual variability and differential hydrolysis of oseltamivir, aspirin, deltamethrin and permethrin. Biochemical Pharmacology 77:238-247. Yoshimura M, Kimura T, Ishii M, Ishii K, Matsuura T, Geshi E, Hosokawa M, and Muramatsu M (2008) Functional polymorphisms in carboxylesterase1A2 (CES1A2) gene involves specific protein 1 Downloaded from (Sp1) binding sites. Biochem Biophys Res Commun 369:939-942. Zhao B, Song J, St Clair RW, and Ghosh S (2007) Stable overexpression of human macrophage cholesteryl ester hydrolase results in enhanced free cholesterol efflux from human THP1 macrophages. Am J Physiol Cell Physiol 292:C405-412. Zhao Z, Li X, Sun S, Mei S, Ma N, Miao Z, Zhao M, and Peng S (2016) Impact of genetic polymorphisms dmd.aspetjournals.org related to clopidogrel or acetylsalicylic acid pharmacology on clinical outcome in Chinese patients with symptomatic extracranial or intracranial stenosis. Eur J Clin Pharmacol 72:1195- 1204. Zhu HJ, Appel DI, Jiang Y, and Markowitz JS (2009a) Age- and Sex-Related Expression and Activity of Carboxylesterase 1 and 2 in Mouse and Human Liver. Drug Metab Dispos 37:1819-1825. Zhu HJ, Appel DI, Johnson JA, Chavin KD, and Markowitz JS (2009b) Role of carboxylesterase 1 and at ASPET Journals on October 1, 2021 impact of natural genetic variants on the hydrolysis of trandolapril. Biochem Pharmacol 77:1266- 1272. Zhu HJ, Appel DI, Peterson YK, Wang ZC, and Markowitz JS (2010) Identification of selected therapeutic agents as inhibitors of carboxylesterase 1: Potential sources of metabolic drug interactions. Toxicology 270:59-65. Zhu HJ, Langaee TY, Gong Y, Wang X, Pepine CJ, Cooper-DeHoff RM, Johnson JA, and Markowitz JS (2016) CES1P1 variant-816A > C is not associated with hepatic carboxylesterase 1 expression and activity or antihypertensive effect of trandolapril. Eur J Clin Pharmacol 72:681-687. Zhu HJ and Markowitz JS (2009) Activation of the Antiviral Prodrug Oseltamivir Is Impaired by Two Newly Identified Carboxylesterase 1 Variants. Drug Metab Dispos 37:264-267. Zhu HJ, Patrick KS, and Markowitz JS (2011) Enantiospecific determination of DL-methylphenidate and DL-ethylphenidate in plasma by liquid chromatography-tandem mass spectrometry: Application to human ethanol interactions. J Chromatogr B 879:783-788. Zhu HJ, Patrick KS, Straughn AB, Reeves OT, 3rd, Bernstein H, Shi J, Johnson HJ, Knight JM, Smith AT, Malcolm RJ, and Markowitz JS (2017) Ethanol Interactions With Dexmethylphenidate and dl- Methylphenidate Spheroidal Oral Drug Absorption Systems in Healthy Volunteers. J Clin Psychopharmacol 37:419-428. Zhu HJ, Patrick KS, Yuan HJ, Wang JS, Donovan JL, DeVane CL, Malcolm R, Johnson JA, Youngblood GL, Sweet DH, Langaee TY, and Markowitz JS (2008) Two CES1 gene mutations lead to dysfunctional carboxylesterase 1 activity in man: Clinical significance and molecular basis. Am J Hum Genet 82:1241-1248. Zhu HJ, Wang X, Gawronski BE, Brinda BJ, Angiolillo DJ, and Markowitz JS (2013) Carboxylesterase 1 as a determinant of clopidogrel metabolism and activation. J Pharmacol Exp Ther 344:665-672. DMD Fast Forward. Published on December 23, 2019 as DOI: 10.1124/dmd.119.089680 This article has not been copyedited and formatted. The final version may differ from this version.

DMD # 89680 40

Zou J-J, Chen S-L, Fan H-W, Tan J, He B-S, and Xie H-G (2014) The CES1A -816C as a genetic marker to predict greater platelet clopidogrel response in patients with percutaneous coronary intervention. Journal of Cardiovascular Pharmacology 63:178-183.

Downloaded from dmd.aspetjournals.org at ASPET Journals on October 1, 2021 DMD Fast Forward. Published on December 23, 2019 as DOI: 10.1124/dmd.119.089680 This article has not been copyedited and formatted. The final version may differ from this version.

DMD # 89680 41

LEGENDS FOR FIGURES

Figure 1. CES1 gene structure and haplotypes. CES1 gene consists of 14 exons located on chromosome 16q13-q22.1, and CES1P1 is a pseudogene, lying tail-to-tail with CES1. CES1, CES1P1 and their variants CES1 VAR and CES1P1 VAR form four major haplotypes. Red represents where stop codon is located at. Transcription efficiency of CES1P1 VAR is approximately 2% of CES1.

Figure 2. D-methylphenidate comes as a single active ingredient (Focalin®) or in Downloaded from combination with l-methylphenidate (racemic mixture) (Ritalin®). D-methylphenidate is approximately 10 times more pharmacologically potent than l-methylphenidate, while l- dmd.aspetjournals.org methylphenidate is a better CES1 substrate. Ethylphenidate can be formed via transesterification with ethanol.

Figure 3. Clopidogrel metabolic pathway. Clopidogrel is a non-ester-prodrug that needs at ASPET Journals on October 1, 2021 to be activated by two oxidation reactions via CYPs. Clopidogrel and its intermediate and active metabolites are all metabolized (deactivated) by CES1.

Figure 4. ACE inhibitors (enalapril and Trandolapril) metabolism. Enalapril and trandolapril are ester-prodrugs that need to be activated by CES1.

Figure 5. Dabigatran metabolic pathway. Dabigatran is a prodrug that activated by both

CES1 and CES2.

DMD Fast Forward. Published on December 23, 2019 as DOI: 10.1124/dmd.119.089680 This article has not been copyedited and formatted. The final version may differ from this version.

DMD # 89680 42

TABLES

Table 1. List of CES1 Substrates. *Prodrugs that need CES1 activation

ACE Inhibitors CNS Agents Antihyperlidpidemia Agents Enalapril* Methylphenidate Clofibrate Imidapril* Cocaine Fenofibrate Benzapril* Heroin Quinapril* Mepridine Adrenal Glucocorticoid Ramipril* Flumazenil Ciclesonide* Trandolapril* Rufinamide Antiviral Agents Anticancer Agnet Chemical Warfare Agents Downloaded from Oseltamivir* Capecitabine* Sarin Sofosbuvir* Irinotecan* Tenofovir alafenamide* Telotristat etiprate* Tabun Endogenous Compounds Antiplatlets/Anticoagulants Pesticides Cholesterol Clopidogrel trans-permethrin dmd.aspetjournals.org Fatty acid ethyl esters Dabigatran* Para-nitrophenyl valerate ARNi Immunosuppressive Agents Others Sacubitril* Mycophenolate mofetil* Dimethyl fumarate* Oxybutynin at ASPET Journals on October 1, 2021

DMD Fast Forward. Published on December 23, 2019 as DOI: 10.1124/dmd.119.089680 This article has not been copyedited and formatted. The final version may differ from this version.

DMD # 89680 43

Table 2. In vitro kinetics of wild type CES1 in human liver S9 fractions (HLS9). HLS9 Vmax Catalytic Efficiency CES1 Substrates Reference (pmol/min/mg Km (μM) (Vmax/Km, μl/min/mg protein) protein) Clopidogrel 3558.6 62.7 56.8

2-oxoclopidogrel (Zhu et al., (clopidogrel 158.1 2.4 65.9 2013) intermediate) Enalapril 67.5 60.1 1.1

Ramipril 18100 690.4 26.2 Downloaded from (Wang et Perindopril 18100 1767 23.3 al., 2016a) Moexipril 4400 1457 12.7 Fosinopril 1400 471.3 3.0 l-methylphenidate 1701 775.7 2.2

(Zhu et al., dmd.aspetjournals.org d-methylphenidate 177.2 663.5 0.3 2008) (Zhu and Oseltamivir 145000 1380 105.1 Markowitz, 2009) (Zhu et al., Trandolapril 103600 639.9 161.9 at ASPET Journals on October 1, 2021 2009b) (Laizure et Dabigatran 1174 33.5 35.0 al., 2014)

DMD Fast Forward. Published on December 23, 2019 as DOI: 10.1124/dmd.119.089680 This article has not been copyedited and formatted. The final version may differ from this version.

DMD # 89680 44 Table 3. CES1 SNPs and their impacts on the PK and PD of CES1 substrate medications Downloaded from AA/Nucleotide Citation Treatment Population Design/Outcome Result Conclusion Change (db SNP ID) G143E Patrick et Methylphenidate n=20 (with 1 Prospective study The study unexpectedly The later analysis

(rs71647871) al., 2007; carrier) found one volunteer withdmd.aspetjournals.org found this volunteer Zhu et al., Single dose of Aim: To examine the elevated PK parameters had G143E and 2008b 0.3mg/kg Healthy interaction between of methylphenidate; Cmax D260fs SNPs, which volunteers methylphenidate and of l-methylphenidate was resulted in elevated alcohol 100 fold higher (62 ng/mL) plasma concentration compared to the rest of of methylphenidate

participants at ASPET Journals on October 1, 2021 Nemoda Methylphenidate n=122 (with 7 Retrospective study G143E carriers needed G143E impaired et al., carriers) lower doses of methylphenidate 2009 Dose adjusted Outcome: methylphenidate for metabolism in vivo based on Hungarian methylphenidate dose symptom reduction symptom ADHD reduction compared to non-carriers reduction x 1 patients (0.410 vs 0.572 mg/kg , month p=0.022) Stage et Methylphenidate n=22 (with 6 Open labeled, G143E carriers showed G143E carriers had al., 2017a carriers) prospective, PK study 152.4% higher AUC higher exposure to Single dose (53.3ng×ml-1×h-1) methylphenidate 10mg Healthy compared to the non- compared to non- Danish carrier group (21.4 ng×ml- carriers 1 -1 Volunteers ×h ) (P<0.0001) Lewis et Clopidogrel n=506 (with 7 Retrospective (1) A 50% higher active G143E carriers had al., 2013 carriers) subanalysis of two metabolite concentration higher plasma (1) PAPI clinical studies: was observed in G143E concentrations of patients (1) Pharmacogenomics carriers (n=7, 30.3 ng/ml) clopidogrel active received 300mg of Antiplatelet compared to non-carriers metabolites, and LD with 75mg Intervention (PAPI) (n=499, 19.0 ng/mL) consequently had a MD Study (P=0.001) higher antiplatelet (2)Clopidogrel-treated effect (2) Patients n=350 (with 6 patients at Sinai from Sinai carriers) Hospital Hospital received either Patient going Outcome: (2) The inhibition of ADP- 300 or 600 mg through PCI (1) PK parameter: induced platelet LD (n= 204) Clopidogrel and its aggregation effect was with MD, or just active metabolite 24% higher in G143E DMD Fast Forward. Published on December 23, 2019 as DOI: 10.1124/dmd.119.089680 This article has not been copyedited and formatted. The final version may differ from this version.

DMD # 89680 45 received 75mg concentration carriers (reduced to 71%

MD (2)PD parameter: ADP- from baseline) comparedDownloaded from stimulated platelet to non-carriers (reduced aggregation to 57% from baseline) (P=0.003) Tarkiainen Clopidogrel n=22 (with 10 Prospective, PK/PD (1) AUC0–∞ ratios of the G143E carriers had

et al., carriers) study clopidogrel carboxylicdmd.aspetjournals.org acid higher exposure to 2015a Single-dose 600 to clopidogrel was 53% clopidogrel active mg Healthy PD outcome: Inhibition less in G143E carriers metabolite, and volunteers of P2Y12-mediated (P=0.009) consequently had a platelet aggregation (2) Average inhibition of higher antiplatelet P2Y12-mediated platelet aggregation effect aggregation in the carriers was 19 percentage points at ASPET Journals on October 1, 2021 higher in non-carriers (P = 0.036) Tarkiainen Enalapril, n=22 (with 10 Prospective PK study (1) AUC0–∞ of the G143E carriers had a et al., Quinapril carriers) enalapril active metabolite lower enalaprilat 2015b enalaprilat was 20% lower exposure compared to Single Dose 10 Healthy in the G143E carriers non-carriers mg Enalapril or volunteers (n=10) compared to non- Quinapril carriers (n=12) (P=0.049) (2) AUCs0-∞ of the quinapril and its active metabolite (i.e., quinaprilat) were not significantly different between the non-carriers and carriers (P=0.114) Tarkiainen Oseltamivir n=22 (with 9 Prospective PK Study (1) The AUC0–∞ ratio of G143E carriers had et al., G143E oseltamivir carboxylate less exposure to 2012 Single Dose heterozygotes, (active metabolite) to oseltamivir active 75mg 1 G143E oseltamivir (parent metabolite compared homozygote ) molecule) was 23% lower to non-carriers in G143E heterozygotes Healthy compared to non-carriers volunteers (P = 0.006) (2) The one G143E homozygous individual had an AUC0–∞ of oseltamivir that was DMD Fast Forward. Published on December 23, 2019 as DOI: 10.1124/dmd.119.089680 This article has not been copyedited and formatted. The final version may differ from this version.

DMD # 89680 46 approximately 360%

greater than the non-Downloaded from carriers

Shi et al., Dabigatran n=102 human In vitro study with The activation rates of G143E carriers had a 2016c liver samples human liver samples DABE, M1 and M2 in lower dabigatran G143E carriers were dmd.aspetjournals.org 53% activation rate, (P = 0.018), 43% (P = potentially resulting in 0.004), and 37% (P = a lower dabigatran 0.001) of normal carriers active metabolite (after normalized by CES1 concentration in the expression) carriers

Shi et al., Sacubitril n= 53 (with 5 In vitro study with The activation rates of at ASPET Journals on October 1, 2021 G143E carriers had a 2016d carrier) human human liver samples sacubitril were lower in lower sacubitril liver samples the carriers compared to activation rate, the non-carriers (4.2 vs. potentially resulting in 7.2 nmol/mg protein/min, a lower sacubitril P = 0.025) active metabolite plasma concentration in the carriers E220G Oh et al., Oseltamivir n=20 (with 8 Prospective, PK study AUC0-48h of oseltamivir E220G appears to (rs200707504) 2017 carriers) was increased by 10% (P have no significant 75mg single = 0.334) and AUC0-48h of impact on oseltamivir dose oseltamivir carboxylate activation in humans was decreased by 5% (P = 0.513) in carriers S75N Xiao et al., Clopidogrel n=851 (with Retrospective PD CES1 S75N carriers S75N appears to (rs2307240) 2017 372 carriers) analysis (n=372) had more increase the function 75mg x 1 year cerebrovascular events of CES1, resulting in Outcome: (P<0.001), acute the decreased efficacy Cerebrovascular myocardial infarction of clopidogrel events, acute (P<0.001), and unstable myocardial infarction, angina (P<0.001) and unstable angina compared to non-carriers Johnson Methylphenidate n=44 (with 2 Naturalistic, prospective No significant differences S75N does not appear et al., carriers) study in methylphenidate side to affect the function of 2013 Weight based effect were found between CES1 dosing x 6 Children with Outcome: Side effect carriers and non-carriers weeks ADHD reported via behavioral (p = 1) questionnaires DMD Fast Forward. Published on December 23, 2019 as DOI: 10.1124/dmd.119.089680 This article has not been copyedited and formatted. The final version may differ from this version.

DMD # 89680 47 Wang et Enalapril n= 36 (with 3 In vitro study with No statistical difference in S75N does not appear

al., 2017 carriers) human liver samples enalapril activation rateDownloaded from or to affect the function of CES1 protein expression CES1 level between carriers and no carriers -816A>C Geshi et imidapril n=105 (with Prospective clinical Greater systolic blood -816A>C appears to

(rs3785161) al., 2005 47 carriers) study pressure reduction (24.1dmd.aspetjournals.org up-regulate the 5-10 mg x 8 mmHg) was observed CES1P1 VAR weeks Patients with compared to non-carriers expression hypertension (17.6 mmHg) after 8 weeks of imidapril therapy (P = 0.0184) Xie et al., Clopidogrel n=162 (with Retrospective PD The carriers had -816A>C appears to 2014 75 carriers) analysis decreased Vasodilator at ASPET Journals on October 1, 2021 up-regulate the 300 or 600 mg Stimulated CES1P1 VAR (LD) or 75 mg Patient on Outcome: Vasodilator Phosphoprotein-Platelet expression (MD) for dual Stimulated Reactivity Index (VASP- minimum 5 days antiplatelet Phosphoprotein-Platelet PRI) (45.93 vs 53.18%) therapy (i.e., Reactivity Index (VASP- (P=0.014) aspirin and PRI) to measure clopidogrel) platelet reactivity with coronary heart diseases

Zou et al., Clopidogrel n=249 (with Retrospective PD A lower ADP-induced -816A>C appears to 2014 108 analysis maximum platelet down-regulate the 300mg LD + heterozygous aggregation (21.5%, CES1P1 VAR 75mg MD x 3 carrier, 17 Outcome: maximum n=125) was observed expression days homozygous platelet aggregation compared to non-carriers carrier) (MPA) (31.7%, n=124) (P=0.001)

Patient going throught PCI Zhu et al., Trandolapril (1) n=486 (1) Retrospective (1) No association -816A>C does not 2016 (with 109 analysis of the between the -816A>C and appears to be 2-4mg x 104 heterozygous INternational VErapamil the blood pressure associated with overall weeks carriers, 10 SR Trandolapril Study lowering effect of CES1 function homozygous trandolapril carriers) (2) In vitro study with (2) Not associated with (2) n=100 (in human liver samples CES1 protein expression vitro study) and trandolapril activation (26 in human liver samples DMD Fast Forward. Published on December 23, 2019 as DOI: 10.1124/dmd.119.089680 This article has not been copyedited and formatted. The final version may differ from this version.

DMD # 89680 48 heterozygous Outcome: Blood

carriers, 3 Pressure Downloaded from homozygous carriers) -75G>T Bruxel et Methylphenidate n=205 (with Retrospective PD The carriers had worse -75G>T appears to be (rs3815583) al., 2013 129 carriers) analysis appetite reduction associated with

Dose titrated up compared to non-carriersdmd.aspetjournals.org decreased CES1 x 3 months as Outcome: appetite (41% vs 77%, P=0.01) function tolerable reduction - the side effect of methylphenidate Sai et al., Irinotecan n=177 Retrospective PK The carriers had higher -75G>T appears to be 2010 analysis plasma (SN-38 + SN- associated with higher 100 mg m−2 177 Japanese 38G)/irinotecan AUC at ASPET Journals on October 1, 2021 CES1 function weekly or 150 cancer ratios relative to non- mg m−2 patients carriers (P=0.027) biweekly Nelveg- ACEI n=200 Retrospective PD The -75G>T genotypes -75G>T was not Kristensen analysis did not significantly impact associated with CES1 et al., Congestive the plasma ATII/ATI ratios function 2016 heart failure in the study subjects or patients fatal outcomes (i.e. cardiovascular death and all-cause death) 1168-33C>A Pare et Dabigatran n=2944 (with Retrospective GWAS of The carriers had lower 1168-33C>A appears (rs2244613) al., 2013 587 carriers) Randomized Evaluation trough concentrations of to be associated with 110 or 150 mg of Long-term the active metabolite lower CES1 function twice daily Patients with Anticoagulation (15% decrease per allele; atrial Therapy (RE-LY) P=1.2×10-8) and a lower fibrillation clinical trial risk of any bleeding (odds (within 6 ratio, 0.67; P=7×10-5) months) and Outcome: Trough compared to non-carriers additional risk concentrations of factors for dabigatran, bleeding stroke risk Shi et al., Dabigatran n=102 (with In vitro study with No association between 1168-33C>A appears 2016a 29 human liver samples 1168-33C>A and to be not associated heterozygous dabigatran activation with CES1 function carriers and 5 homozygous carriers) DMD Fast Forward. Published on December 23, 2019 as DOI: 10.1124/dmd.119.089680 This article has not been copyedited and formatted. The final version may differ from this version.

DMD # 89680 49

Table 4. Drug-Drug Interaction Summary

CES1 inhibitors CES1 Substrates Interaction Summary Alcohol Methylphenidate Many in vitro and in vivo studies confirmed alcohol inhibits CES1 and mediates biotransformation of methylphenidate to ethylphenidate; methylphenidate plasma concentrations were increased when patients took methylphenidate with alcohol (Griffin et al., 2010; Bell et al., 2011a; Bell et al., 2011b; Zhu et al., 2011; Griffin et al.,

2013; Patrick et al., 2013; Parker et al., Downloaded from 2015; Zhu et al., 2017). Alcohol Oseltamivir When alcohol was administered with oseltamivir in humans, the AUC of oseltamivir increased by 37% (Hu et al.,

2014). dmd.aspetjournals.org Cannabis Oseltamivir In vitro study with CES1-transfected cells suggested THC, CBD and CBN to be the potent CES1 inhibitors. The inhibition constant (Ki) values for THC, CBD and CBN were 0.541, 0.974, and 0.263 µM

(0.170, 0.306, and 0.0817 µg/ml) (Qian et at ASPET Journals on October 1, 2021 al., 2019). Protease Methylphenidate. p- In vitro study showed that protease Inhibitors nitrophenyl acetate inhibitors (i.e., nelfinavir, amprenavir, (PNPA) and p- atazanavir, ritonavir, and saquinavir) nitrophenol (PNP) inhibited the catalytic activity of CES1 (p <0.01). Among protease inhibitors, nelfinavir had a significantly higher inhibitory effect compared to other agents (Rhoades et al., 2012). Aripiprazole Methylphenidate, p- In vitro study suggested aripiprazole, nitrophenyl acetate perphenazine, thioridazine, and (PNPA) fluoxetine to be potent inhibitors of CES1. Among the medications tested, aripiprazole was the most potent inhibitor of CES1, and an in vivo study with FVB mouse confirmed this result (Zhu et al., 2010). Isradipine PNPA, Trandolapril In vitro study with human liver /Tacrolimus microsomes suggested isradipine (dihydropyridine calcium antagonist, DHP) and tacrolimus DMD Fast Forward. Published on December 23, 2019 as DOI: 10.1124/dmd.119.089680 This article has not been copyedited and formatted. The final version may differ from this version.

DMD # 89680 50

(immunosuppressive agent) to be potent CES1 inhibitors (Thomsen et al., 2014). Valproate Rufinamide In vitro study suggested valproate could inhibit CES1 function and affect rufinamide metabolism (Williams et al., 2011). ACEI Clopidogrel ACEIs and clopidogrel are often administered together as both of them are cardiovascular medications; both ACEIs and clopidogrel are suggested to be inhibitors of CES1. A clinical study with myocardial infarction patients

(n=70,934) demonstrated concomitant Downloaded from use of ACEIs increased the rate of clinically significant bleeding compared to the clopidogrel monotherapy (P=0.002) (Kristensen et al., 2014). Another clinical

study with myocardial infarction patients dmd.aspetjournals.org (n=45,918) with clopidogrel showed that concomitant use of clopidogrel and ACEI (perindopril and ramipril) was not associated with the re-infarction, heart failure or death (Cressman et al., 2015).

at ASPET Journals on October 1, 2021

DMD Fast Forward. Published on December 23, 2019 as DOI: 10.1124/dmd.119.089680 This article has not been copyedited and formatted. The final version may differ from this version.

DMD # 89680

Figure 1

Downloaded from dmd.aspetjournals.org

at ASPET Journals on October 1, 2021

DMD Fast Forward. Published on December 23, 2019 as DOI: 10.1124/dmd.119.089680 This article has not been copyedited and formatted. The final version may differ from this version.

DMD # 89680

Figure 2

Downloaded from dmd.aspetjournals.org

at ASPET Journals on October 1, 2021 DMD Fast Forward. Published on December 23, 2019 as DOI: 10.1124/dmd.119.089680 This article has not been copyedited and formatted. The final version may differ from this version.

DMD # 89680

Figure 3

Downloaded from dmd.aspetjournals.org at ASPET Journals on October 1, 2021

DMD Fast Forward. Published on December 23, 2019 as DOI: 10.1124/dmd.119.089680 This article has not been copyedited and formatted. The final version may differ from this version.

DMD # 89680

Figure 4

Downloaded from dmd.aspetjournals.org at ASPET Journals on October 1, 2021

DMD Fast Forward. Published on December 23, 2019 as DOI: 10.1124/dmd.119.089680 This article has not been copyedited and formatted. The final version may differ from this version.

DMD # 89680

Figure 5

Downloaded from dmd.aspetjournals.org at ASPET Journals on October 1, 2021