Skip to main content
Advertisement

Main menu

  • Home
  • Articles
    • Current Issue
    • Fast Forward
    • Latest Articles
    • Archive
  • Information
    • Instructions to Authors
    • Submit a Manuscript
    • FAQs
    • For Subscribers
    • Terms & Conditions of Use
    • Permissions
  • Editorial Board
  • Alerts
    • Alerts
    • RSS Feeds
  • Virtual Issues
  • Feedback
  • Other Publications
    • Drug Metabolism and Disposition
    • Journal of Pharmacology and Experimental Therapeutics
    • Molecular Pharmacology
    • Pharmacological Reviews
    • Pharmacology Research & Perspectives
    • ASPET

User menu

  • My alerts
  • Log in
  • My Cart

Search

  • Advanced search
Journal of Pharmacology and Experimental Therapeutics
  • Other Publications
    • Drug Metabolism and Disposition
    • Journal of Pharmacology and Experimental Therapeutics
    • Molecular Pharmacology
    • Pharmacological Reviews
    • Pharmacology Research & Perspectives
    • ASPET
  • My alerts
  • Log in
  • My Cart
Journal of Pharmacology and Experimental Therapeutics

Advanced Search

  • Home
  • Articles
    • Current Issue
    • Fast Forward
    • Latest Articles
    • Archive
  • Information
    • Instructions to Authors
    • Submit a Manuscript
    • FAQs
    • For Subscribers
    • Terms & Conditions of Use
    • Permissions
  • Editorial Board
  • Alerts
    • Alerts
    • RSS Feeds
  • Virtual Issues
  • Feedback
  • Visit jpet on Facebook
  • Follow jpet on Twitter
  • Follow jpet on LinkedIn
Research ArticleMetabolism, Transport, and Pharmacogenetics
Open Access

Targeting Enteropeptidase with Reversible Covalent Inhibitors To Achieve Metabolic Benefits

Weimei Sun, Xuqing Zhang, Maxwell D. Cummings, Kamal Albarazanji, Jiejun Wu, Mina Wang, Richard Alexander, Bin Zhu, YueMei Zhang, James Leonard, James Lanter and James Lenhard
Journal of Pharmacology and Experimental Therapeutics December 2020, 375 (3) 510-521; DOI: https://doi.org/10.1124/jpet.120.000219
Weimei Sun
DPDS Discovery Technology and Molecular Pharmacology, Spring House, Pennsylvania (W.S., M.W., R.A.); DPDS Analytical Sciences, La Jolla, California (J.W.); Discovery Chemistry, Spring House, Pennsylvania (X.Z., M.D.C., B.Z., Y.Z., J.La.); CVM Discovery, Spring House, Pennsylvania (K.A., J.Leo., J.Len.); and Janssen Research & Development, LLC, Spring House, Pennsylvania
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Xuqing Zhang
DPDS Discovery Technology and Molecular Pharmacology, Spring House, Pennsylvania (W.S., M.W., R.A.); DPDS Analytical Sciences, La Jolla, California (J.W.); Discovery Chemistry, Spring House, Pennsylvania (X.Z., M.D.C., B.Z., Y.Z., J.La.); CVM Discovery, Spring House, Pennsylvania (K.A., J.Leo., J.Len.); and Janssen Research & Development, LLC, Spring House, Pennsylvania
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Maxwell D. Cummings
DPDS Discovery Technology and Molecular Pharmacology, Spring House, Pennsylvania (W.S., M.W., R.A.); DPDS Analytical Sciences, La Jolla, California (J.W.); Discovery Chemistry, Spring House, Pennsylvania (X.Z., M.D.C., B.Z., Y.Z., J.La.); CVM Discovery, Spring House, Pennsylvania (K.A., J.Leo., J.Len.); and Janssen Research & Development, LLC, Spring House, Pennsylvania
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Kamal Albarazanji
DPDS Discovery Technology and Molecular Pharmacology, Spring House, Pennsylvania (W.S., M.W., R.A.); DPDS Analytical Sciences, La Jolla, California (J.W.); Discovery Chemistry, Spring House, Pennsylvania (X.Z., M.D.C., B.Z., Y.Z., J.La.); CVM Discovery, Spring House, Pennsylvania (K.A., J.Leo., J.Len.); and Janssen Research & Development, LLC, Spring House, Pennsylvania
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Jiejun Wu
DPDS Discovery Technology and Molecular Pharmacology, Spring House, Pennsylvania (W.S., M.W., R.A.); DPDS Analytical Sciences, La Jolla, California (J.W.); Discovery Chemistry, Spring House, Pennsylvania (X.Z., M.D.C., B.Z., Y.Z., J.La.); CVM Discovery, Spring House, Pennsylvania (K.A., J.Leo., J.Len.); and Janssen Research & Development, LLC, Spring House, Pennsylvania
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Mina Wang
DPDS Discovery Technology and Molecular Pharmacology, Spring House, Pennsylvania (W.S., M.W., R.A.); DPDS Analytical Sciences, La Jolla, California (J.W.); Discovery Chemistry, Spring House, Pennsylvania (X.Z., M.D.C., B.Z., Y.Z., J.La.); CVM Discovery, Spring House, Pennsylvania (K.A., J.Leo., J.Len.); and Janssen Research & Development, LLC, Spring House, Pennsylvania
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Richard Alexander
DPDS Discovery Technology and Molecular Pharmacology, Spring House, Pennsylvania (W.S., M.W., R.A.); DPDS Analytical Sciences, La Jolla, California (J.W.); Discovery Chemistry, Spring House, Pennsylvania (X.Z., M.D.C., B.Z., Y.Z., J.La.); CVM Discovery, Spring House, Pennsylvania (K.A., J.Leo., J.Len.); and Janssen Research & Development, LLC, Spring House, Pennsylvania
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Bin Zhu
DPDS Discovery Technology and Molecular Pharmacology, Spring House, Pennsylvania (W.S., M.W., R.A.); DPDS Analytical Sciences, La Jolla, California (J.W.); Discovery Chemistry, Spring House, Pennsylvania (X.Z., M.D.C., B.Z., Y.Z., J.La.); CVM Discovery, Spring House, Pennsylvania (K.A., J.Leo., J.Len.); and Janssen Research & Development, LLC, Spring House, Pennsylvania
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
YueMei Zhang
DPDS Discovery Technology and Molecular Pharmacology, Spring House, Pennsylvania (W.S., M.W., R.A.); DPDS Analytical Sciences, La Jolla, California (J.W.); Discovery Chemistry, Spring House, Pennsylvania (X.Z., M.D.C., B.Z., Y.Z., J.La.); CVM Discovery, Spring House, Pennsylvania (K.A., J.Leo., J.Len.); and Janssen Research & Development, LLC, Spring House, Pennsylvania
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
James Leonard
DPDS Discovery Technology and Molecular Pharmacology, Spring House, Pennsylvania (W.S., M.W., R.A.); DPDS Analytical Sciences, La Jolla, California (J.W.); Discovery Chemistry, Spring House, Pennsylvania (X.Z., M.D.C., B.Z., Y.Z., J.La.); CVM Discovery, Spring House, Pennsylvania (K.A., J.Leo., J.Len.); and Janssen Research & Development, LLC, Spring House, Pennsylvania
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
James Lanter
DPDS Discovery Technology and Molecular Pharmacology, Spring House, Pennsylvania (W.S., M.W., R.A.); DPDS Analytical Sciences, La Jolla, California (J.W.); Discovery Chemistry, Spring House, Pennsylvania (X.Z., M.D.C., B.Z., Y.Z., J.La.); CVM Discovery, Spring House, Pennsylvania (K.A., J.Leo., J.Len.); and Janssen Research & Development, LLC, Spring House, Pennsylvania
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
James Lenhard
DPDS Discovery Technology and Molecular Pharmacology, Spring House, Pennsylvania (W.S., M.W., R.A.); DPDS Analytical Sciences, La Jolla, California (J.W.); Discovery Chemistry, Spring House, Pennsylvania (X.Z., M.D.C., B.Z., Y.Z., J.La.); CVM Discovery, Spring House, Pennsylvania (K.A., J.Leo., J.Len.); and Janssen Research & Development, LLC, Spring House, Pennsylvania
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • Article
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF
Loading

Abstract

Inhibition of the serine protease enteropeptidase (EP) opens a new avenue to the discovery of chemotherapeutics for the treatment of metabolic diseases. Camostat has been used clinically for treating chronic pancreatitis in Japan; however, the mechanistic basis of the observed clinical efficacy has not been fully elucidated. We demonstrate that camostat is a potent reversible covalent inhibitor of EP, with an inhibition potency (kinact/KI) of 1.5 × 104 M−1s−1. High-resolution liquid chromatography-mass spectrometry (LC-MS) showed addition of 161.6 Da to EP after the reaction with camostat, consistent with insertion of the carboxyphenylguanidine moiety of camostat. Covalent inhibition of EP by camostat is reversible, with an enzyme reactivation half-life of 14.3 hours. Formation of a covalent adduct was further supported by a crystal structure resolved to 2.19 Å, showing modification of the catalytic serine of EP by a close analog of camostat, leading to formation of the carboxyphenylguanidine acyl enzyme identical to that expected for the reaction with camostat. Of particular note, minor structural modifications of camostat led to changes in the mechanism of inhibition. We observed from other studies that sustained inhibition of EP is required to effect a reduction in cumulative food intake and body weight, with concomitant improved blood glucose levels in obese and diabetic leptin-deficient mice. Thus, the structure-activity relationship needs to be driven by not only the inhibition potency but also the mechanistic and kinetic characterization. Our findings support EP as a target for the treatment of metabolic diseases and demonstrate that reversible covalent EP inhibitors show clinically relevant efficacy.

SIGNIFICANCE STATEMENT Interest in targeted covalent drugs has expanded in recent years, particularly so for kinase targets, but also more broadly. This study demonstrates that reversible covalent inhibition of the serine protease enteropeptidase is a therapeutically viable approach to the treatment of metabolic diseases and that mechanistic details of inhibition are relevant to clinical efficacy. Our mechanistic and kinetic studies outline a framework for detailed inhibitor characterization that is proving essential in guiding discovery efforts in this area.

Introduction

Enteropeptidase (EP, EC. 3.4.21.9) is a serine protease exclusively localized to the apical membrane of the proximal small intestine (Yuan et al., 1998; Imamura and Kitamoto, 2003). The initiation of intestinal nutrient digestion occurs in the lumen of the upper duodenum through the proteolytic processing of trypsinogen to catalytically active trypsin by EP (Light and Janska, 1989; Kitamoto et al., 1995; Zheng et al., 2009). This initiates a cascade of events leading to zymogen activation and engagement in digestive processes to regulate nutrient signaling, absorption, and metabolism. A brush-border enzyme, duodenase, has been characterized (Zamolodchikova et al., 1997) and shown to activate recombinant bovine EP (Zamolodchikova et al., 2000); the same research group reported an analogous protein from human Brunner’s glands (Zamolodchikova et al., 2013). However, a different group speculated that trypsin, present in pancreatic fluid, is the true physiologic activator of proenteropeptidase (Zheng et al., 2009). This latter contention was supported by data showing that the rate of activation of purified bovine recombinant proenteropeptidase by purified duodenase was approximately 100-fold slower than that by trypsin (Zamolodchikova et al., 2000). The process by which EP is activated has not been thoroughly characterized and remains unclear.

Full-length human EP consists of a proposed transmembrane domain, followed by heavy and light chains linked by a disulfide bond. The light chain comprises a trypsin-like serine protease domain with a typical Asp-His-Ser catalytic triad (Light and Janska, 1989; Zheng et al., 2009). The structural basis for the observed classic trypsin-like substrate specificity at P1 (Schechter and Berger, 1967) as well as the unusual P2-P4 Asp-Asp-Asp recognition sequence has been thoroughly discussed (Lu et al., 1999). The protease domain of EP shows high structural similarity with and >40% sequence identity to related proteases, including matriptase, plasma kallikrein, hepsin, FXIa, epitheliasin, and prostasin. In addition to the N-terminal transmembrane domain that serves to anchor EP in the brush border, the heavy chain includes several additional structural domains of varying function (Kitamoto et al., 1994; Zheng et al., 2009).

Camostat (Fig. 1) was reported to inhibit trypsin-like serine proteases. Diabetic rats treated with camostat showed significant reduction in blood glucose and body weight as well as improved insulin resistance (Ito et al., 2007). Camostat is clinically approved for treating pancreatitis and reflux esophagitis in Japan. Patients treated with camostat for pancreatitis for 8 years showed significantly lower incidence of diabetes (Kitagawa and Hayakawa, 2007; Bronze-da-Rocha and Santos-Silva, 2018; Ramsey et al., 2019). However, the pharmacological mechanism underlying the observed clinical efficacy has not been fully elucidated.

Fig. 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 1.

Molecular structure of camostat.

We observed that camostat is a potent inhibitor of EP. Our study establishes EP as a viable therapeutic target for obesity and diabetes based on pharmacologic results observed with leptin-deficient mice (ob/ob mice) treated with camostat and camostat-like inhibitors. Additional evidence has recently emerged illustrating that inhibition of EP can lead to therapeutically useful metabolic effects in animal models (Sasaki et al., 2019; Yashiro et al., 2019).

To leverage the established safety profile of this clinically approved molecule, an extensive chemical design effort around the camostat scaffold was undertaken to identify new EP inhibitors with increased inhibition potency and enhanced pharmacological properties. Our approach is to target protein digestion in the duodenum by inhibiting EP locally. This gastrointestinal tract–restricted modality offers both the challenge of durably inhibiting EP and the opportunity to mitigate potential off-target effects by limiting drug exposure to the alimentary organs.

The mechanism of EP inhibition by camostat was thoroughly characterized to establish a benchmark and outline a framework for kinetic and mechanistic characterization aimed at reliable and pharmacologically relevant ranking of new inhibitors. The results of our structure-activity relationship (SAR) and pharmacological studies will be reported in separate publications. In the present work, we focus on kinetic and mechanistic characterization of EP inhibition by camostat as well as new inhibitors we are evaluating. Our explorations have yielded a variety of EP inhibitors with differentiated mechanisms. It is clearly demonstrated that detailed in vitro biochemical characterization of inhibition kinetics and mechanism is critical to SAR generation and comprehension, which in turn is essential in the selection of lead compounds most likely to show the desired pharmacological effects.

Materials and Methods

Recombinant full-length human EP, expressed and purified from CHO cells, was purchased from PeproTech (catalog number 450-48C). Recombinant light chain of bovine EP (residues 801–1035), expressed and purified from Escherichia coli, was purchased from R&D Systems (catalog number 4139-SE). This form of the light chain of EP contains the C-terminal fragment of the heavy chain, Cys788–Lys800 with an N-terminal Ala, and Ile801-His1035 of the light chain. Peptide Gly-Asp-Asp-Asp-Asp-Lys-β-naphthylamide, purity greater than 95% by high performance liquid chromatography (HPLC) analysis, was purchased from Sigma (catalog number G5261); the powder was dissolved in DMSO and stored at −20°C. β-Naphthylamine (NA, catalog number 31618), camostat mesylate (catalog number SML0057), Triton X-100 (catalog number 93443), and urea (catalog number 1084880250) were all purchased from Sigma-Aldrich. NaCl solution (5 M; catalog number S0250) and 1 M Tris-HCl at pH 7.4 (catalog number T1074) were purchased from Teknova. LC-MS–grade B&J Brand water (catalog number LC365-4), water with 0.1% formic acid (catalog number LC452-1), and acetonitrile with 0.1% formic acid (catalog number LC441-1) were purchased from Honeywell International Inc. (Muskegon, MI). Anhydrous DMSO was purchased from EMD Millipore Corporation (catalog number MX1457-7). The 384-well nonbinding surface black plates, used for biochemical assays, were purchased from Perkin Elmer (catalog number 6007279). A Labcyte Echo 550 was used for nanoliter compound transfers, an Infinite M1000 plate reader (Tecan Group) was used for kinetic reading of fluorescence intensity, and an Eppendorf centrifuge 5810 was used to spin assay plates to ensure thorough mixing of solutions and to remove bubbles generated from pipetting. Micro Bio-SpinTM P-6 Gel columns (catalog number 732-6200) were purchased from Bio-Rad Laboratories.

In Vitro Biochemical Assay.

Recombinant full-length human EP was used in the biochemical assay, with the quenched fluorescently labeled peptide Gly-Asp-Asp-Asp-Asp-Lys-β-naphthylamide as substrate in assay buffer [50 mM Tris-HCl (pH 7.4), 200 mM NaCl, 0.01% Triton X-100]. Assay buffer was prepared from purchased stock solutions and was diluted with LC-MS–grade water. Progress curves of protease activity were measured by monitoring the fluorescence intensity of the cleavage product NA at excitation and emission wavelengths of 340 and 410 nm, respectively.

Inhibition of EP.

Concentration-dependent inhibition studies involved mixing equal volumes of 600 μM substrate and 2.5 nM EP solutions. Compounds were dissolved to 10 mM in DMSO immediately prior to use in inhibition studies, and 2-fold serial dilutions from 10 mM were made with DMSO, for a total of 11 inhibitor concentrations; a DMSO control (no inhibitor) was also included. Inhibitor solutions were transferred to assay plates with an Echo acoustic liquid dispenser. Reaction was initiated by addition of the enzyme solution to the mixture of substrate and compound. Reaction progress curve data were collected continuously at 1-minute intervals for 60 minutes. The observed inactivation rate constant (kobs) was derived by fitting the reaction progress curve to eq. 1.Embedded Image(1)YT is the detection signal of the reaction at time t; vi and vs are the initial and the steady-state velocities of the reaction, respectively; and t is the reaction time. Y0 is the background signal at t = 0.

The derived inactivation rate constant (kobs) is dependent on inhibitor concentration. A linear dependence indicates a one-step mechanism (Scheme A), with the slope of the linear fit equivalent to the second-order rate constant of compound binding (k1), which is used to represent the compound inhibition potency (kinact/KI). Alternatively, a hyperbolic fit of kobs versus [I] indicates a two-step binding mechanism (Scheme B), wherein the plateau of the curve is the forward rate constant of the isomerization step (k3), and KI is the inhibitor concentration that gives half-maximal inactivation. Inhibition potency (kinact/KI) is reported as the ratio k3/KI, with substrate competition corrected for by the term (1 + [S]/KM).

Kinetic Data Analysis.

Kinetic data were fit to the appropriate rate equations by GraphPad Prism 7.0.

Confirmation of Covalent Modification of EP by LC-MS.

The recombinant light chain of bovine EP (residues 801–1035) was studied by LC-MS to explore covalent modification of EP. Samples for LC-MS study were mixtures of 1 µM EP and 10 µM inhibitor in assay buffer, with a control sample that contained EP alone.

Intact protein LC-MS data were acquired on an Agilent 1290 UHPLC/6550 qTOF system. Two Phenomenex Aeris widepore XB-C18 columns (50 × 2.1 mm, 3.6 μm; catalog number 00B-4482-AN) were used: one was maintained at 40°C and used with the normal sample injection, whereas the other column was maintained at 70°C and the sample injection was preceded with 20 μl of 7 M urea aqueous solution injected to create a denaturing condition. The mobile phases used were (A) water with 0.1% formic acid and (B) acetonitrile with 0.1% formic acid; the flow rate was 0.4 ml/min. The LC linear gradient segments were as follows: 0 minutes, 5% B; 4 minutes, 99% B; 5 minutes, 99% B; 5.5 minutes, 5% B; 6 minutes, 5% B. The qTOF was set at MS acquisition electrospray positive ion mode and calibrated to within 5 ppm mass error. The MS acquisition range was from 300 to 2500 m/z with a 2 Hz scanning rate. Intact protein molecular weight was deconvoluted from the raw data by using the Agilent MassHunter (version 4.0) Maximal Entropy deconvolution method. A ligand or ligand fragment covalently attached to EP shows an increase in molecular weight, and the m/z measured with both columns is expected to be equivalent if the attachment is covalent and stable.

Inhibition Half-Life.

Hydrolysis of EP-compound adduct was determined by measuring recovery of EP activity as a function of time. Samples were made by mixing EP (0.51 μM final) with an excess amount of compound (23.3 μM final) in the assay buffer and incubating at room temperature for 50 minutes. As a control, EP was incubated with DMSO. After incubation, 75 μl of the sample was transferred to a pre-equilibrated Bio-Rad microspin column to remove free compound. Eluted EP-compound adduct was immediately diluted 100-fold in reaction buffer by mixing with 300 and 600 μM (final concentrations) of substrate. EP activity was observed to recover slowly with time. The time traces were fit to eq. 1 to derive the first-order rate constant for reactivation. The final steady-state velocity was calculated from the control sample and set as a fixed parameter in the curve, fitting for the derivation of the dissociation rate constant of compound, kobs. The half-life (t1/2) of EP recovery was derived from the value of ln (2)/kobs, i.e., 0.693/kobs.

Protein Crystallography.

Supercharged light-chain human EP (residues 785–1019) containing amino acid substitutions N6D, G21D, G22D, C112S, N142D, and K210E was expressed in E. coli, refolded from inclusion bodies, and purified as described previously (Simeonov et al., 2012). The protein at 10 mg/ml was crystallized from 22% (w/v) PEG20,000, 0.1 M HEPES/NaOH (pH 7.25), by mixing 0.5 µl of the protein with 0.5 µl of the reservoir solution using hanging drop vapor diffusion at 18°C. For complex formation, apo crystals were soaked at 10 mM from a 100 mM DMSO stock solution of compound 6 for 3 hours in reservoir solution at pH 6.5 at 20°C. Diffraction data were collected at the Diamond Light Source at beamline I02 (Didcot, Oxfordshire, UK) at a wavelength of 0.9159 Å. Data were processed with XDS and XSCALE (Kabsch, 2010) and CCP4 (Winn et al., 2011). The complex structure was solved by molecular replacement using the program MOLREP (Vagin and Teplyakov, 2010) and refined with REFMAC5 (Murshudov et al., 2011). Model building was performed with COOT (Emsley et al., 2010).

Effect of Camostat on Food Intake, Body Weight, and Blood Glucose in Obese and Diabetic ob/ob Mice.

Seven-week-old male leptin-deficient ob/ob mice (Jax Laboratories, Bar Harbor, ME) were used. Animals were single-housed in a temperature-controlled room with a 12-hour light/dark cycle and allowed ad libitum access to water and powdered chow (modified Purina #5008 irradiated chow with 2% maltose dextrin provided by Dyets, Inc.). Mice were randomized into five groups (n = 8 mice per group) based on baseline average body weights (48.1 ± 0.4 g) and 5-hour fasted blood glucose levels [238.1 ± 9.7 mg/dl, measured with an AlphaTrak2 Glucometer (Abbott, Abbott Park, IL)]. After 4 days of recording baseline food intake and body weights, mice were fed powdered chow containing camostat (w/w) at 0.0, 0.08, 0.25, 0.8 mg/g food (formulated by Dyets, Inc.) for 7 days. In the same study, a group of mice was pair-fed to the high-dose (0.8 mg/g food) group, with these animals receiving the same amount of daily food as the camostat-treated group but fed a camostat-free diet.

Chemistry.

1H NMR spectra were determined with a Bruker Biospin International AG-300 or AG-400 spectrometers at 300 or 400 MHz. Chemical shifts (δ) are reported in parts per million relative to residual chloroform (7.26 ppm), tetramethylsilane (TMS) (0 ppm), or CD3OD (4.87 ppm) as an internal reference with coupling constants (J) reported in hertz. The peak shapes are denoted as follows: singlet (s); doublet (d); triplet; quartet; multiplet (m); broad (br). Electrospray mass spectra were recorded in positive or negative mode on a Micromass Platform spectrometer. The purity of all compounds described here was determined by analytical LC-MS using a Shimadzu LCMS-2110EV with analytical LC using Shiseido Capcell Pak C18 MGS3 column (3.0 × 50 mm, 3.5 μm), 0.05% trifluoroacetic acid (TFA) in water as mobile phase A and 0.05% TFA in acetonitrile as mobile phase B at 1 ml/min flow, gradient from 10% to 100% B in 1.7 minutes, 100% B for 1.5 minutes, 100%–10% B in 0.13 minutes, monitored by UV absorption at 215 nm using photodiode array (PDA) and evaporative light scattering detector (ELSD). The purity of all compounds was found to be >95%. Thin-layer chromatography was performed on Merck PLC prescored plates 60F254. Prep-HPLC conditions: (1#waters 2767-5) column, SunFire Prep C18, 19*150 mm H Prep C-001(T) 18600256819513816414 04; mobile phase, phase A: water with 0.05% TFA, phase B: CH3CN (20% CH3CN up to 50% in 8 minutes, up to 100% in 0.1 minutes, hold 100% in 1.9 minutes, down to 20% in 0.1 minutes, hold 20% in 1.9 minutes); detector, UV 220 and 254 nm. Unless otherwise noted, reagents were obtained from commercial sources and were used without further purification.

Compound 11f (methyl 4-​carbamimidamidobenzo​ate) was purchased from Enamine with >95% purity. Compound 6 was prepared according to PCT International Application (1997), WO 9737969 A1 19971016.

Preparation of 2-(3-Fluoro-4-((4-Guanidino Benzoyl) Oxy) Phenyl) Acetic Acid (11a).

Into a 100-ml round-bottom flask was placed 2-(3-fluoro-4-hydroxyphenyl) acetic acid (3.0 g, 17.6 mmol, 1.0 Eq) in dimethylformamide (DMF) (30 ml); NaH (60% in mineral oil, 776 mg, 19.4 mmol, 1.1 Eq) was added in portions at 0°C. The mixture was stirred for 20 minutes and followed by the addition of benzyl bromide (4.51 g, 26.4 mmol, 1.5 Eq) dropwise. The mixed solution was stirred for 4 hours at 25°C. The reaction was then quenched by the addition of water. The resulting solution was extracted with ethyl acetate three times, and the organic layers were combined. The resulting solution was concentrated under vacuum. The residue was purified by CombiFlash system with silica gel using petroleum ether/ethyl acetate = 4:1 solution. This afforded 1.9 g of benzyl 2-(3-fluoro-4-hydroxyphenyl) acetate as white solid. MS (ESI): mass calculated for C15H13FO3; m/z, 260.3; found, 259.0 [M − H]. Into a 250-ml round-bottom flask was placed 4-guanidinobenzoic acid (1.32 g, 6.15 mmol, 2.0 Eq) in anhydrous pyridine (26 ml), and the mixture solution was stirred for 5 minutes at room temperature. It was followed by the addition of benzyl 2-(3-fluoro-4-hydroxy phenyl) acetate (800 mg, 3.07 mmol, 1.0 Eq) and 4-methylbenzenesulfonic acid (106 mg, 0.62 mmol, 0.2 Eq). After 15 minutes, a solution of dicyclohexylmethanediimine (1.90 g, 9.22 mmol, 3.0 Eq) in dry DMF (13 ml) was added at once, and the mixture was stirred at room temperature for 15 hours. The precipitated solid was filtered and washed with water (20 ml). The resulting solution was extracted with ethyl acetate, and the organic layers were combined. The resulting solution was concentrated under vacuum. The residue was purified by a CombiFlash system with silica gel using dichloromethane (DCM)/MeOH = 5:1 solution. This afforded 800 mg of 4-(2-(benzyloxy)-2-oxoethyl)-2-fluorophenyl 4-guanidino benzoate as white solid. MS (ESI): mass calculated for C23H20FN3O4; m/z, 421.4; found, 422.0 [M + H]+.

Into a 250-ml round-bottom flask charged and maintained with an inert atmosphere of nitrogen was placed 4-(2-(benzyloxy)-2-oxoethyl)-2-fluorophenyl 4-guanidino benzoate (200 mg, 0.48 mmol, 1.0 Eq), MeOH (10 ml), 5% palladium on carbon (Pd/C) (60 mg). H2 was introduced into the mixture at 1 atmosphere. The mixture solution was stirred for 2 hours at 25°C. The resulting mixture was filtered through a pad of Celite. The filtrate was concentrated under vacuum and recrystallized from the mixture of methanol and dichloromethane, which afforded 87.3 mg (95% yield) of 2-(3-fluoro-4-(4-guanidinobenzoyloxy) phenyl) acetic acid as white solid. MS (ESI): mass calculated for C16H14FN3O4; m/z, 331.3; found, 331.9 [M + H]+. 1H NMR (300 MHz, CD3OD) δ = 8.29 (d, J = 2.1 Hz, 2H), 7.50 (d, J = 2.1 Hz, 2H), 7.27–7.30 (m, 2H), 7.20–7.25 (m, 1H), 3.69 (s, 2H).

2-(4-((2-Chloro-4-Guanidinobenzoyl) Oxy) Phenyl) Acetic Acid (7).

Analogous synthesis to 11a. Mass spectrum (ESI): mass calculated for C16H14ClN3O4, 347.7, found 348.0 [M + H]+. 1H NMR (300 MHz, DMSO-d6) δ = 8.10 (d, J = 8.4 Hz, 1H), 7.64–8.03 (m, 2H), 7.35–7.51 (m, 3H), 7.22–7.30 (m, 1H), 7.10–7.21 (m, 2 H), 3.56 (s, 2H).

2-(4-((5-Guanidinopicolinoyl) Oxy) Phenyl) Acetic Acid (8).

Analogous synthesis to 11a. MS (ESI): mass calculated for C15H14N4O4; m/z, 314.3; found, 315.3 [M + H]+. 1H NMR (400 MHz, CD3OD) δ = 7.05 (s, 1H), 6.72 (d, J = 6.5 Hz, 1H), 6.29 (d, J = 7.0 Hz, 1H), 5.72 (d, J = 8.1 Hz, 2H), 5.55 (d, J = 7.9 Hz, 2H), 1.98 (s, 2H).

2-(4-((6-Guanidinonicotinoyl) Oxy) Phenyl) Acetic Acid (9).

Analogous synthesis to 11a. MS (ESI): mass calculated for C15H14N4O4; m/z, 314.3; found, 315.3 [M + H]+. 1H NMR (400 MHz, CD3OD) δ = 7.58 (s, 1H), 6.98 (m, 1H), 5.85 (d, J = 8.0 Hz, 2H), 5.68 (d, J = 8.0 Hz, 2H), 5.60 (m, 1H), 2.10 (s, 2H).

2-(3-Chloro-4-((4-Guanidinobenzoyl) Oxy) Phenyl) Acetic Acid (11b).

Analogous synthesis to 11a. MS (ESI): mass calculated for C16H14ClN3O4; m/z, 347.8; found, 348.9 [M + H]+. 1H NMR (300 MHz, D2O) δ = 8.10–8.13 (m, 2H), 7.32–7.47 (m, 3H), 7.12–7.19 (m, 2H), 3.54 (s, 2H).

2-(3-Trifluoromethyl-4-((4-Guanidino Benzoyl) Oxy) Phenyl) Acetic Acid (11c).

Analogous synthesis to 11a. MS (ESI): mass calculated for C17H14F3N3O4; m/z, 381.3; found, 382.1 [M + H]+. 1H NMR (400 MHz, CD3OD) δ = 8.24–8.26 (m, 2H), 7.72 (s, 1H), 7.66 (d, J = 8.8 Hz, 1H), 7.46–7.49 (m, 2H), 7.42 (d, J = 8.4 Hz, 1H), 3.76 (s, 2H).

4-Chlorophenyl 4-Guanidinobenzoate (11d).

Analogous synthesis to 11a. MS (ESI): mass calculated for C14H12ClN3O2; m/z, 289.7; found, 290.1 [M + H]+. 1H NMR (400 MHz, CD3OD) δ = 7.75 (d, J = 9.5 Hz, 2H), 7.55 (d, J = 10.5 Hz, 2H), 7.48 (d, J = 8.5 Hz, 2H), 6.88 (d, J = 9.5 Hz, 2H).

2-((4-Guanidinobenzoyl) Oxy) Benzoic Acid (11e).

Analogous synthesis to 11a. MS (ESI): mass calculated for C15H13N3O4; m/z, 299.3; found, 300.0 [M + H]+. 1H NMR (400 MHz, CD3OD) δ = 8.31 (d, J = 7.8 Hz, 1H), 8.06 (m, 1H), 7.99 (d, J = 8.1 Hz, 1H), 7.85 (m, 1H), 7.80 (d, J = 9.0 Hz, 2H), 6.94 (d, J = 9.0 Hz, 2H).

1,1,1,3,3,3-Hexafluoropropan-2-yl 4-Guanidinobenzoate (11g).

Analogous synthesis to 11a. MS (ESI): mass calculated for C11H9F6N3O2; m/z, 329.2; found, 330.3 [M + H]+. 1H NMR (400 MHz, DMSO-d6) δ = 10.56 (br, s, 1H), 8.11 (d, J = 7.5 Hz, 2H), 7.92 (s, br, 3H), 7.48 (d, J = 7.5 Hz, 2H), 7.10 (m, 1H). 19F NMR (377 MHz, DMSO-d6) δ = −72.5.

Preparation of 2-(4-((4-Carbamimidoyl Phenoxy) Carbonyl) Phenyl) Acetic Acid (5).

Into a 500-ml flask was added methyl 4-iodobenzoate (10 g, 38.2 mmol, 1 Eq), CuI (727 mg, 3.82 mmol, 0.1 Eq), Cs2CO3 (37.3 g, 114.5 mmol, 3 Eq), and 2-picolinic acid (931.9 mg, 7.63 mmol, 0.2 Eq). The flask was exchanged with N2, then dioxane (100 ml) and tert-butyl 3-oxobutanoate (12.3 ml, 76.3 mmol, 2 Eq) were added. The resulting mixture was stirred for 72 hours at 70°C. The solvent was evaporated under reduced pressure, and then the mixture was diluted with ethyl acetate (EtOAc) (150 ml). The mixture was washed with saturated NH4Cl (150 ml). The organic layer was dried over anhydrous Na2SO4 and concentrated under vacuum. The residue was purified with silica gel column chromatography (petroleum ether/ethyl acetate from 100:0 to 70:30) to afford methyl 4-(2-tert-butoxy-2-oxoethyl) benzoate as yellow oil (4.7 g, 49% yield). MS (ESI): mass calculated for C14H18O4; m/z, 250.3; found, 195.1 [M-t-Bu + H]+.

To a stirred solution of methyl 4-(2-tert-butoxy-2-oxoethyl) benzoate (1.6 g, 6.4 mmol, 1 Eq) in THF/H2O (v/v = 2/1, 30 ml) was added LiOH⋅H2O (402 mg, 9.6 mmol, 1.5 Eq) at 0°C. The mixture was stirred overnight while warming up to room temperature. The mixture was diluted with H2O and then extracted with DCM (100 ml) to remove unreacted starting material. The aqueous phase was acidified to pH 1–2 with 1 N HCl solution, extracted with EtOAc (3 × 100 ml). The combined organic phase was dried over Na2SO4, filtered, and concentrated under vacuum. The residue was purified with silica gel column chromatography (DCM/MeOH from 100:0 to 90:10) to afford 4-(2-tert-butoxy-2-oxoethyl) benzoic acid as white solid (650 mg, 43% yield). MS (ESI): mass calculated for C13H16O4; m/z, 236.3; found, 235.0 [M − H]−.

To a stirred solution of 4-(2-tert-butoxy-2-oxoethyl) benzoic acid (300 mg, 1.27 mmol, 1 Eq) in pyridine was added 4-hydroxybenzimidamide (345 mg, 2.54 mmol, 2 Eq), followed by dicyclohexylmethanediimine (393 mg, 1.91 mmol, 1.5 Eq). The reaction was stirred at room temperature overnight. The solvent was evaporated under vacuum, and the residue was purified by silica gel column chromatography (70:30 EtOAc/petroleum ether, then flush from 0:100 to 15:85 MeOH/DCM mixed solution) to afford 4-carbamimidoylphenyl 4-(2-tert-butoxy-2-oxoethyl)benzoate as white solid (350 mg, 77.7% yield). MS (ESI): mass calculated for C20H22N2O4; m/z, 354.4; found, 355.2 [M + H]+.

To a 100-ml flask was added 4-carbamimidoylphenyl 4-(2-tert-butoxy-2-oxoethyl) benzoate (350 mg, 0.99 mmol) and HCl (4 M in dioxane, 15 ml), and the mixture was stirred overnight at room temperature. The solvent was evaporated under vacuum, and the solid was recrystallized from cold ether to afford 2-(4-((4-carbamimidoylphenoxy) carbonyl) phenyl) acetic acid as white solid (125.1 mg, 36.8% yield). MS (ESI): mass calculated for C16H14N2O4; m/z, 298.3; found, 299.1 [M + H]+. 1H NMR (300 MHz, CD3OD) δ = 8.16 (d, J = 8.4 Hz, 2H), 7.92 (d, J = 8.7 Hz, 2H), 7.45–7.55 (m, 4H), 3.77 (s, 2H).

4-Carbamimidoyl-2-Fluorophenyl Benzoate (4).

Analogous synthesis to 5. Mass spectrum (ESI, m/z): mass calculated for C14H11FN2O2, m/z, 258.3; found, 259.1 [M + H]+. 1H NMR (300 MHz, CD3OD) δ = 8.21 (d, J = 7.5 Hz, 2H), 7.74–7.88 (m, 3H), 7.46–7.68 (m, 3H).

Preparation of 2-(4-((4-Carbamimidoyl Piperazine-1-Carbonyl) Oxy) Phenyl) Acetic Acid (10).

To a solution of benzyl 2-(4-hydroxy phenyl) acetate (1 g, 4.13 mmol, 1.0 Eq) and 4-nitrophenyl carbonochloridate (874 mg, 4.33 mmol, 1.05 Eq) in DCM (50 ml) was added triethylamine (1.25 g, 12.4 mmol, 3.0 Eq) at 0°C. The mixture was warmed to room temperature and stirred for 1 hour. Then, tert-butyl piperidin-4-ylcarbamate (1.24 g, 6.2 mmol, 1.5 Eq) was added into the mixture. The solution was stirred for 2 hours at room temperature. The reaction was then quenched by the addition of water. The resulting solution was extracted with DCM, and the organic layers were combined. The resulting mixture was concentrated under vacuum. The residue was purified with silica gel column chromatography (petroleum ether/EtOAc 2:1) to afford 1.2 g of 1-(4-(2-(benzyloxy)-2-oxoethyl) phenyl) 4-tert-butyl piperazine-1,4-dicarboxylate as white solid. Mass spectrum (ESI, m/z): mass calculated for C25H30N2O6, m/z, 454.5; found, 477.3 [M + Na]+.

Into a 100-ml flask was added 1-(4-(2-(benzyloxy)-2-oxoethyl) phenyl) 4-tert-butyl piperazine-1,4-dicarboxylate (1.2 g, 2.64 mmol, 1.0 Eq) and 20 ml of 10% TFA solution in DCM. The mixture was stirred overnight at room temperature. The mixture was concentrated in vacuum to afford 4-(2-(benzyloxy)-2-oxoethyl) phenyl piperazine-1-carboxylate (862 mg, 69.7% yield) as brown oil. Mass spectrum (ESI, m/z): mass calculated for C20H22N2O4, m/z, 354.4; found, 355.1 [M + H]+.

Into a 250-ml flask was added 4-(2-(benzyloxy)-2-oxoethyl)phenyl piperazine-1-carboxylate (620 mg, 1.75 mmol, 1 Eq), acetonitrile (MeCN) (20 ml), and N,N-diisopropylethylamine (DIPEA) (1.53 ml, 8.75 mmol, 5 Eq), followed by 1H-pyrazole-1-carboximidamide hydrochloride (308 mg, 2.1 mmol, 1.2 Eq). The mixture was stirred at room temperature overnight. The solvent was evaporated in vacuum, and the residue was dissolved in EtOAc (100 ml) and then washed with 2 N HCl (100 ml). The organic layer was dried over Na2SO4, filtered, and concentrated in vacuum. The residue was purified with silica gel column chromatography (MeOH/DCM from 0:100 to 15:85) to afford 4-(2-(benzyloxy)-2-oxoethyl) phenyl 4-carbamimidoylpiperazine-1-carboxylate (501 mg, 72.2% yield) as white solid. Mass spectrum (ESI, m/z): mass calculated for C21H24N4O4, m/z, 396.5; found, 397.2 [M + H]+.

To a 100-ml flask was added 4-(2-(benzyloxy)-2-oxoethyl) phenyl 4-carbam imidoylpiperazine-1-carboxylate (494 mg, 1.25 mmol), 5% Pd/C (200 mg) in MeOH (30 ml). The air in the flask was exchanged with 1 atmosphere of H2, and the mixture was stirred at room temperature under H2 overnight. Pd/C was filtered through a pad of Celite, and the filtrate was evaporated under vacuum. The residue was recrystallized from DMF/HCl/acetonitrile (MeCN) mixed solution to afford 2-(4-(4-carbamimidoylpiperazine-1-carbonyl oxy) phenyl) acetic acid as white solid (124.3 mg, 27.7% yield). Mass spectrum (ESI, m/z): mass calculated for C14H18N4O4, m/z, 306.3; found, 307.1 [M + H]+. 1H NMR (300 MHz, D2O) δ = 7.26 (d, J = 8.1 Hz, 2H), 7.03 (d, J = 8.4 Hz, 2H), 3.76 (br, 2H), 3.66 (s, 2H), 3.44–3.53 (m, 6H).

Results

EP cleaves peptide and protein substrates N-terminal to P1 (Schechter and Berger, 1967) Lys or Arg residues and prefers Asp or Glu at the P2 position. Both the full-length and light-chain–only forms of EP recognize the activation peptide of trypsinogen, Gly-Asp-Asp-Asp-Asp-Lys (GD4K) (Davie and Neurath, 1953; Davie and Neurath, 1955). An enzyme assay with the fluorogenic peptide D4K-NA as substrate was used to determine the protease activity of EP by following the fluorescence intensity of the cleavage product NA at excitation and emission wavelengths of 340 and 410 nm, respectively (Grant and Hermon-Taylor, 1979). Characterization of the protease activity of full-length human EP was studied under steady-state conditions; the Michaelis-Menten constant KM and the apparent turnover number kcat were determined to be 412 ± 22 μM and 11.8 ± 0.4 s−1, respectively (Fig. 2), with the catalytic efficiency of the enzyme therefore 2.86 × 104 M−1s−1. Kinetic characterization of the light-chain–only form of EP showed kinetic properties and catalytic efficiencies identical to those of the full-length enzyme (data not shown).

Fig. 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 2.

Kinetic characterization of EP. (A) EP reaction progress curves were followed in the presence of various concentrations of substrate GDDDDK-NA and 1.25 nM enzyme in buffer of 50 mM Tris-HCl (pH 7.4), 200 mM NaCl, and 0.01% Triton X-100. The steady-state reaction velocities were determined by taking the tangents of the linear curves. (B) EP activity was plotted as a function of substrate concentration. Data were fit to the Michaelis-Menten equation. An average of 412 ± 22 μM and 11.8 ± 0.4 s−1 were obtained for the Michaelis-Menten constant KM and the apparent turnover number kcat, respectively. Data are the averages of two independent experiments.

Inhibition of EP by Camostat.

Camostat inhibits the enzymatic activity of EP in a time-dependent manner (Fig. 3A), with the observed inhibition following a one-step mechanism (Fig. 3B). The apparent inhibition potency (kinact/KIapp) was determined to be 8.7 × 103 M−1s−1, and kinact/KI was therefore 1.5 × 104 M−1s−1 following modification by the competing factor of (1+[S]/KM), which is 1.728 for our assay conditions. The fact that the plot in Fig. 3B bisects the origin indicates that either the inhibition is irreversible or that it is reversible with a long residence time.

Fig. 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 3.

Inhibition of EP by camostat. (A) Time-dependent inhibition of EP by camostat. Inhibition was determined in an assay using 300 µM substrate and 1.25 nM EP in the presence of various concentrations of camostat. First-order inactivation rate constants (kobs) were derived from curve fitting to eq. 1 and then plotted against camostat concentration. (B) Camostat inhibition of EP was shown to follow a one-step mechanism, and the apparent compound potency (kinact/KIapp) was calculated to be 8.7 × 103 M−1s−1. RFU stands for the relative fluorescence units.

The results of an LC-MS study indicated covalent modification of EP by camostat. Preincubation of truncated EP (light chain of EP) with camostat showed an addition of 161.6 Da to the protein, consistent with complex formation between the enzyme and the carboxyphenylguanidine moiety of camostat (Fig. 4; see also Fig. 5).

Fig. 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 4.

High-resolution mass spectra of EP protein with and without camostat. (A) EP protein only; (B) EP protein after incubation with camostat, LC-MS run at 40°C; (C) EP protein plus camostat, LC-MS run at 70°C after one injection of 7 M urea. The addition of 161.6 Da to EP after incubation with camostat was seen, and this addition corresponds to the camostat fragment in the dotted box (B), and the attachment is covalent. MW is the molecular weight.

Fig. 5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 5.

Crystal structure of the EP-compound 6 reaction product adduct, resolved to 2.19 Å (deposited at the PDB with ID code 6ZOV). The carboxymethylphenylguanidine is covalently bound to the catalytic nucleophile Ser187, with the guanidine moiety projecting down into S1 and forming a salt bridge with Asp181; the ester carbonyl does not project directly toward the oxyanion hole region but rather is directed away from the protein and toward bulk solvent. Residue numbering as in the published EP crystal structure (PDB ID 4DGJ). EP is shown as a gray ribbon, with selected side chains shown as purple color-by-atom sticks and the inhibitor moiety shown as green color-by-atom sticks. Selected atom-atom contacts are highlighted with black dashed lines. PDB refers to the Protein Data Bank.

Interestingly, the enzymatic activity of EP could be recovered from fully inhibited EP after an overnight dialysis, consistent with slow release of active enzyme from the covalent adduct (data not shown). The rate of hydrolysis of acylated EP was determined by monitoring the recovery of enzyme activity over time (Fig. 6). The t1/2 of the recovery of EP activity from fully inhibited enzyme was therefore calculated to be 14.3 ± 3.8 hours and was independent of the concentration of substrate used in the experiment.

Fig. 6.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 6.

Recovery of EP from camostat inhibition. EP activity was measured as a function of time in the presence of substrate at 300 and 600 µM, respectively, by a rapid dilution of covalently modified EP. The rate constant of acylated-EP hydrolysis was derived from data fitting analysis, and half-life of the EP recovery was calculated to be 14.3 ± 3.8 hours (average of n = 4, duplicate at each concentration of substrate). FI stands for the fluorescence intensity.

Protein Crystallography.

Initially, we used a publicly available crystal structure of a prostasin-camostat complex (Spraggon et al., 2009) as a reference in advance of our own EP structural results. Prostasin, a trypsin-like serine protease, shows 39% sequence identity across the protease domain and >70% sequence identity in the active site region (results not shown). Based on the extensive structural information available for phenylguanidine and related moieties occupying S1 of trypsin-like serine proteases, we were confident in predicting a similar binding mode for camostat with EP, including covalent labeling of the catalytic serine. This prediction was fully consistent with the results of our mechanistic studies, which are described in the preceding sections, and was subsequently supported by our own structural results (Tables 1 and 2).

View this table:
  • View inline
  • View popup
TABLE 1

Data collection and processing statistics for the EP-compound 6 complex

View this table:
  • View inline
  • View popup
TABLE 2

Refinement statistics for EP-compound 6 complexa

Compound 6 is a close chemical analog of camostat but lacks the terminal dimethylaminoethyl ester moiety of the leaving group of the parent molecule. Furthermore, kinetic mechanistic studies indicate that compound 6 follows the same inhibitory mechanism as camostat (see Fig. 3 and Table 4; inhibition residence time of compound 6 is not shown here). Thus, we expect the structures of the reaction products of these two inhibitors with EP to be identical. The complex structure of the reaction product of the 3-hour soak of EP with compound 6 was determined to a resolution of 2.19 Å (Fig. 5). This structure is fully consistent with the earlier structure of the prostasin-camostat complex mentioned above (Spraggon et al., 2009). In the EP complex, the carboxymethylphenylguanidine moiety occupies the S1 subsite, forming the expected salt bridge with the side chain of Asp181 at the base of this pocket. The side chain of the catalytic residue Ser187 has formed a covalent bond in the observed acyl-enzyme intermediate, unambiguously establishing the covalent nature of inhibition of EP by this class of inhibitors. The carbonyl oxygen of the adduct, positioned near to the oxyanion hole formed by Gly185 and Ser187, is directed away from the protein surface toward bulk solvent.

Mechanism of Inhibition.

The scheme in Fig. 7 summarizes the mechanism of inhibition of EP by camostat. Camostat inhibits EP by forming a covalent adduct (1) between the nucleophile of EP active site serine residue and the acyl moiety of camostat. The acylated EP (1) is fully inactive and slowly hydrolyzed to release active EP with a t1/2 of 14.3 hours.

Fig. 7.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 7.

Scheme of camostat inhibition of EP. Inhibition undergoes a covalent yet slowly reversible mechanism.

SAR Generation of Camostat-Based Chemical Design.

To target EP, exclusively expressed in the small intestine, a series of camostat analogs as biopharmaceutical classification system (BCS) class III structures with high solubility and poor permeability for targeting gut restriction were synthesized and characterized kinetically and mechanistically (Table 4). Removal of the dimethylaminoethyl group in camostat afforded zwitterionic structures, which were poorly permeable to minimize systemic exposure and maximize gastrointestinal tract exposure. Inhibition of EP by these compounds was consistently observed to be time-dependent. Compound potencies were ranked by kinact/KI values, as well as Ki values (derived from the IC50 based on initial velocities). Table 3 provides a telling example, with two analogs having equal potencies (kinact/KI) showing a 5-fold difference in the inhibition Ki. These data demonstrate that accurate SAR for time-dependent inhibitors can only be represented by the value of kinact/KI derived from measuring the inactivation kinetics of the enzyme. Traditional Ki (or IC50) values alone could confuse or mislead the interpretation of SAR studies.

View this table:
  • View inline
  • View popup
TABLE 3

Inhibition of time-dependent inhibitors can only be accurately presented by kinact/KI

SAR of Medicinal Chemistry.

It was observed that structural changes at the compound acyl moiety lead to a change in the inhibition mechanism. Compounds 6 and 7 are two close analogs, differing only at the ortho- position of the guanidine benzoyl ring. Both compounds inhibited EP in a time-dependent fashion with similar potencies (kinact/KI) (Fig. 8). However, compound 6 inhibits EP in a one-step mechanism (Fig. 8A), whereas compound 7 shows a two-step mechanism (Fig. 8B). This is consistent with inhibition of EP by compound 6 being driven largely or completely by reactivity toward EP, whereas both the recognition (or binding) and reactivity components are significant contributors to inhibition by compound 7. Covalent inhibitors that follow a two-step mechanism that includes a significant or major recognition component are highly desired, as they are expected to show superior on-target selectivity and concomitant reduced in vivo off-target liabilities. Conversely, inhibitors driven largely or completely by reactivity are more likely to show a broader range of off-target activities (Cornish-Bowden, 2012; Copeland, 2013).

Fig. 8.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 8.

Analogs of camostat have similar inhibition potencies but follow distinct inhibition mechanisms. (A) Linear dependence of the inactivation rate constant on compound concentration indicated a one-step inhibition mechanism for compound 6; (B) a hyperbolic curve of the plot indicated a two-step inhibition mechanism for compound 7.

Compounds 8 and 9 are regioisomers bearing guanidino-picolinoyl-oxy moieties (Fig. 9), and this relatively minor structural variation in the acyl moiety resulted in a significant difference in the mechanism of inhibition. Linear reaction progress curves in the presence of compound 9 indicated a steady-state fast equilibrium reversible inhibition (Fig. 9B), whereas the nonlinear reaction progress curves reflected time-dependent inhibition by compound 8 (Fig. 9A).

Fig. 9.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 9.

A minor structural change of the inhibitor acyl moiety causes a change of the mechanism of inhibition from (A) thetime-dependent characteristics of compound 8 to (B) the reversible steady-state kinetics of compound 9 as displayed by the inhibition progress curves. FI stands for the fluorescence intensity.

The acyl moiety plays a major role in modulating the compound potency, in addition to its role in regulating the inhibition mechanism. Compound 10 (Fig. 10) incorporates the same leaving group as compounds 8 and 9 (two potent inhibitors in Fig. 9), but replacement of the aromatic acyl moiety with a piperazine carbamate made the compound completely inactive toward EP (<5% inhibition was observed at 50 µM compound).

Fig. 10.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 10.

Molecular structure of compound 10.

The effect of the leaving group of the active esters was also extensively studied (Table 4). The data establish a clear dependence of inhibition potency on the leaving group, while the same one-step inhibition mechanism is conserved. Addition of a 3-F group on the phenyl acetic acid core of compound 6 slightly improved EP potency (11a vs. 6, kinact/KI, 1.68 × 104 vs. 1.02 × 104 M−1s−1). Moreover, analogs 11b and 11c bearing the 3-Cl and 3-CF3 substituents on the phenyl acetic acid core resulted in approximately 3-fold increases in EP potency (kinact/KI, 2.62 × 104 for 11b; 3.32 × 104 M−1 s−1for 11c), indicating electronic withdrawing groups may enhance the enzymatic reactivity of the active esters. Replacement of the acetic acid group with a Cl group on the phenyl ring was well tolerated (11d). In contrast, the para- substitutions between the acid and ester group were essential to EP potency. Compared with 6 bearing para-acetic acid group on the phenyl ring, compound 11e with ortho-carboxylic acid substitution on the phenyl ring significantly attenuated EP potency to kinact/KI, 4.46 × 102 M−1s−1. Removal of the phenyl ring abolished EP potency as shown in the methyl ester 11f. Incorporation of a nonaromatic electronic withdrawing group such as hexafluoropropyl in 11g recovered weak EP activity (kinact/KI, 4.89 × 102 M−1s−1) from the inactive 11f.

View this table:
  • View inline
  • View popup
TABLE 4

SAR of the leaving group of the active esters

Effect of Camostat on Food Intake and Body Weight in ob/ob Mice.

Leptin-deficient ob/ob mice provide a model of obesity and hyperglycemia. The average doses of camostat (0.08, 0.25 and 0.8 mg/g food) intake in feed over the 7 days of treatment were calculated as 9.24, 25.5, and 68.59 mg/kg per day, respectively. Relative to the placebo control group, camostat reduced cumulative food intake by −3.2%, −16.1%, and −34.4% at doses of 9.24, 25.5, and 68.59 mg/kg per day, respectively (see Fig. 11). Camostat reduced the percent change in body weight by −2.8% and −8.3% at doses of 25.5 and 68.59 mg/kg per day, respectively. There was no difference in weight loss between the camostat at high dose and pair-fed group, suggesting that reduced food intake is the primary mechanism for weight loss during the 7 days of treatment (Fig. 12). Baseline 5-hour fasted average blood glucose was 238.1 ± 9.7 mg/dl (n = 40). Treatment with camostat for 7 days reduced 5-hour fasted blood glucose by 15.6% and 36.9% at doses of 25.5 and 68.59 mg/kg per day, respectively. The 5-hour fasted blood glucose was not reduced in the pair-fed group, despite the pair-fed group consuming similar amounts of food and showing weight loss similar to that seen for the camostat-treated group, suggesting that the reduction in body weight or food intake may be not the reason for the improvement in glucose levels in camostat-treated animals (see Fig. 13).

Fig. 11.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 11.

Effects of 7 days of treatment with camostat given in feed on cumulative food intake in ob/ob mice. Data are expressed as means ± S.E.M. (n = 8 per group). *P < 0.05 vs. placebo by one-way ANOVA followed by Dunnett’s multiple comparisons test using Prism Graph Pad 7 software.

Fig. 12.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 12.

Effects of 7 days of treatment with camostat given in feed on percent change in body weight from baseline in ob/ob mice. Data are expressed as means ± S.E.M. (n = 8 per group). *P < 0.05 vs. placebo by one-way ANOVA followed by Dunnett’s multiple comparisons test using Prism Graph Pad 7 software.

Fig. 13.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 13.

Effects of 7 days of treatment with camostat given in feed on Δ 5-hour fasted blood glucose in ob/ob mice. Data are expressed as means ± S.E.M. (n = 8 per group). *P < 0.05 vs. placebo and #P < 0.05 vs. pair-fed group by one-way ANOVA followed by Dunnett’s or Tukey’s multiple comparisons test, respectively, using Prism Graph Pad 7 software.

Discussion

We discovered that camostat is a potent inhibitor of EP, with inhibition following a reversible covalent mechanism and showing a long inhibition half-life of 14.3 hours. A detailed inhibition kinetics study established that camostat inhibition of EP proceeded through a one-step mechanism, consistent with high reactivity toward EP in forming the covalent adduct. The initial binding complex between camostat and EP is weak and kinetically unstable. The mechanism of camostat inhibition of EP explains its poor selectivity among other serine proteases; to achieve high selectivity for EP and minimize in vivo off-target liabilities, an inhibitor with increased binding affinity for EP, therefore following a two-step inhibition mechanism, is preferable. Distinct from an irreversible covalent inhibitor, one advantage of developing a reversible covalent inhibitor drug could be in limiting the provocation of an immune response by avoiding the formation and accumulation of irreversibly modified protein. We have described the surprising discovery of EP inhibitors similar in structure to camostat but showing inhibition mechanisms distinct from that of camostat. We established that sustained inhibition of EP is required to attain a metabolic benefit in animal models (data not shown). In this light, knowledge of the detailed inhibition kinetics and mechanism is proving as important as inhibition potency in guiding our understanding of inhibitor SAR. We have outlined an in vitro biochemical and kinetic characterization framework to guide the discovery of a reversible covalent inhibitor of EP that could lead to new therapies for patients with obesity and diabetes.

Acknowledgments

We thank Janssen Research & Development, LLC, for financial support. The crystal structure of the EP-compound 6 complex (PDB ID 6ZOV) was determined at Proteros biostructures GmbH (Martinsried, Germany), and we appreciate the efforts and support of Georg Kempf, Klaus Maskos, and Stefan Steinbacher at Proteros.

Authorship Contributions

Participated in research design: Sun, X. Zhang, Zhu, Y. Zhang, Leonard, Lanter, Lenhard.

Conducted experiments: Sun, X. Zhang, Cummings, Albarazanji, Wu, Wang.

Performed data analysis: Sun, X. Zhang, Cummings, Albarazanji, Wu, Wang.

Wrote or contributed to the writing of the manuscript: Sun, X. Zhang, Cummings, Albarazanji, Wu.

Footnotes

    • Received July 14, 2020.
    • Accepted September 8, 2020.
  • This paper received no external funding.

  • This work was previously presented as an abstract at the following meeting: Zhang X, Sun W, Lanter J, Zhu B, Cummings MD, Wu J, Albarazanji K, Wang M, Lenhard J, and Leonard J. Targeting enteropeptidase by reversible covalent inhibitors. Abstracts of Papers, 259th ACS National Meeting & Exposition; 2020 March 22–26; Philadelphia, PA. MEDI-0148.

  • https://doi.org/10.1124/jpet.120.000219.

Abbreviations

br
broad
d
doublet
EP
enteropeptidase
J
coupling constant
m
multiplet
NA
β-naphthylamine
ob/ob mice
leptin-deficient mice
s
singlet
SAR
structure-activity relationship
t1/2
half-life
  • Copyright © 2020 The Author(s)

This is an open access article distributed under the CC BY Attribution 4.0 International license.

References

  1. ↵
    1. Bronze-da-Rocha E and
    2. Santos-Silva A
    (2018) Neutrophil elastase inhibitors and chronic kidney disease. Int J Biol Sci 14:1343–1360.
    OpenUrl
  2. ↵
    1. Copeland RA
    (2013) Evaluation of Enzyme Inhibitors in Drug Discovery, 2nd ed, John Wiley, New York.
  3. ↵
    1. Cornish-Bowden A
    (2012) Fundamentals of Enzyme Kinetics, 4th ed, Wiley-VCH Verlag, Berlin.
  4. ↵
    1. Davie EW and
    2. Neurath H
    (1955) Identification of a peptide released during autocatalytic activation of trypsinogen. J Biol Chem 212:515–529.
    OpenUrlFREE Full Text
  5. ↵
    1. Davie EW and
    2. Neurath H
    (1953) Identification of the peptide split from trypsinogen during autocatalytic activation. Biochim Biophys Acta 11:442.
    OpenUrlPubMed
  6. ↵
    1. Emsley P,
    2. Lohkamp B,
    3. Scott WG, and
    4. Cowtan K
    (2010) Features and development of Coot. Acta Crystallogr D Biol Crystallogr 66:486–501.
    OpenUrlCrossRefPubMed
  7. ↵
    1. Grant DA and
    2. Hermon-Taylor J
    (1979) Hydrolysis of artificial substrates by enterokinase and trypsin and the development of a sensitive specific assay for enterokinase in serum. Biochim Biophys Acta 567:207–215.
    OpenUrlPubMed
  8. ↵
    1. Imamura T and
    2. Kitamoto Y
    (2003) Expression of enteropeptidase in differentiated enterocytes, goblet cells, and the tumor cells in human duodenum. Am J Physiol Gastrointest Liver Physiol 285:G1235–G1241.
    OpenUrlCrossRefPubMed
  9. ↵
    1. Ito T,
    2. Otsuki M,
    3. Itoi T,
    4. Shimosegawa T,
    5. Funakoshi A,
    6. Shiratori K,
    7. Naruse S,
    8. Kuroda Y, and Research Committee of Intractable Diseases of the Pancreas
    (2007) Pancreatic diabetes in a follow-up survey of chronic pancreatitis in Japan. J Gastroenterol 42:291–297.
    OpenUrlCrossRefPubMed
  10. ↵
    1. Kabsch W
    (2010) Integration, scaling, space-group assignment and post-refinement. Acta Crystallogr D Biol Crystallogr 66:133–144.
    OpenUrlCrossRefPubMed
  11. ↵
    1. Kitagawa M and
    2. Hayakawa T
    (2007) Antiproteases and the pancreas: basic and clinical update. Introduction. JOP 8(4 Suppl):476–478.
    OpenUrlPubMed
  12. ↵
    1. Kitamoto Y,
    2. Veile RA,
    3. Donis-Keller H, and
    4. Sadler JE
    (1995) cDNA sequence and chromosomal localization of human enterokinase, the proteolytic activator of trypsinogen. Biochemistry 34:4562–4568.
    OpenUrlCrossRefPubMed
  13. ↵
    1. Kitamoto Y,
    2. Yuan X,
    3. Wu Q,
    4. McCourt DW, and
    5. Sadler JE
    (1994) Enterokinase, the initiator of intestinal digestion, is a mosaic protease composed of a distinctive assortment of domains. Proc Natl Acad Sci USA 91:7588–7592.
    OpenUrlAbstract/FREE Full Text
  14. ↵
    1. Light A and
    2. Janska H
    (1989) Enterokinase (enteropeptidase): comparative aspects. Trends Biochem Sci 14:110–112.
    OpenUrlCrossRefPubMed
  15. ↵
    1. Lu D,
    2. Fütterer K,
    3. Korolev S,
    4. Zheng X,
    5. Tan K,
    6. Waksman G, and
    7. Sadler JE
    (1999) Crystal structure of enteropeptidase light chain complexed with an analog of the trypsinogen activation peptide. J Mol Biol 292:361–373.
    OpenUrlCrossRefPubMed
  16. ↵
    1. Murshudov GN,
    2. Skubák P,
    3. Lebedev AA,
    4. Pannu NS,
    5. Steiner RA,
    6. Nicholls RA,
    7. Winn MD,
    8. Long F, and
    9. Vagin AA
    (2011) REFMAC5 for the refinement of macromolecular crystal structures. Acta Crystallogr D Biol Crystallogr 67:355–367.
    OpenUrlCrossRefPubMed
  17. ↵
    1. Ramsey ML,
    2. Nuttall J,
    3. Hart PA, and TACTIC Investigative Team
    (2019) A phase 1/2 trial to evaluate the pharmacokinetics, safety, and efficacy of NI-03 in patients with chronic pancreatitis: study protocol for a randomized controlled trial on the assessment of camostat treatment in chronic pancreatitis (TACTIC). Trials 20:501.
    OpenUrl
  18. ↵
    1. Sasaki M,
    2. Miyahisa I,
    3. Itono S,
    4. Yashiro H,
    5. Hiyoshi H,
    6. Tsuchimori K,
    7. Hamagami KI,
    8. Moritoh Y,
    9. Watanabe M,
    10. Tohyama K, et al.
    (2019) Discovery and characterization of a small-molecule enteropeptidase inhibitor, SCO-792. Pharmacol Res Perspect 7:e00517.
    OpenUrl
  19. ↵
    1. Schechter I and
    2. Berger A
    (1967) On the size of the active site in proteases. I. Papain. Biochem Biophys Res Commun 27:157–162.
    OpenUrlCrossRefPubMed
  20. ↵
    1. Simeonov P,
    2. Zahn M,
    3. Sträter N, and
    4. Zuchner T
    (2012) Crystal structure of a supercharged variant of the human enteropeptidase light chain. Proteins 80:1907–1910.
    OpenUrlPubMed
  21. ↵
    1. Spraggon G,
    2. Hornsby M,
    3. Shipway A,
    4. Tully DC,
    5. Bursulaya B,
    6. Danahay H,
    7. Harris JL, and
    8. Lesley SA
    (2009) Active site conformational changes of prostasin provide a new mechanism of protease regulation by divalent cations. Protein Sci 18:1081–1094.
    OpenUrlCrossRefPubMed
  22. ↵
    1. Vagin A and
    2. Teplyakov A
    (2010) Molecular replacement with MOLREP. Acta Crystallogr D Biol Crystallogr 66:22–25.
    OpenUrlCrossRefPubMed
  23. ↵
    1. Winn MD,
    2. Ballard CC,
    3. Cowtan KD,
    4. Dodson EJ,
    5. Emsley P,
    6. Evans PR,
    7. Keegan RM,
    8. Krissinel EB,
    9. Leslie AG,
    10. McCoy A, et al.
    (2011) Overview of the CCP4 suite and current developments. Acta Crystallogr D Biol Crystallogr 67:235–242.
    OpenUrlCrossRefPubMed
  24. ↵
    1. Yashiro H,
    2. Hamagami K,
    3. Hiyoshi H,
    4. Sugama J,
    5. Tsuchimori K,
    6. Yamaguchi F,
    7. Moritoh Y,
    8. Sasaki M,
    9. Maekawa T,
    10. Yamada Y, et al.
    (2019) SCO-792, an enteropeptidase inhibitor, improves disease status of diabetes and obesity in mice. Diabetes Obes Metab 21:2228–2239.
    OpenUrl
  25. ↵
    1. Yuan X,
    2. Zheng X,
    3. Lu D,
    4. Rubin DC,
    5. Pung CY, and
    6. Sadler JE
    (1998) Structure of murine enterokinase (enteropeptidase) and expression in small intestine during development. Am J Physiol 274:G342–G349.
    OpenUrlCrossRefPubMed
  26. ↵
    1. Zamolodchikova TS,
    2. Scherbakov IT,
    3. Khrennikov BN, and
    4. Svirshchevskaya EV
    (2013) Expression of duodenase-like protein in epitheliocytes of Brunner’s glands in human duodenal mucosa. Biochemistry (Mosc) 78:954–957.
    OpenUrl
  27. ↵
    1. Zamolodchikova TS,
    2. Sokolova EA,
    3. Alexandrov SL,
    4. Mikhaleva II,
    5. Prudchenko IA,
    6. Morozov IA,
    7. Kononenko NV,
    8. Mirgorodskaya OA,
    9. Da U,
    10. Larionova NI, et al.
    (1997) Subcellular localization, substrate specificity and crystallization of duodenase, a potential activator of enteropeptidase. Eur J Biochem 249:612–621.
    OpenUrlPubMed
  28. ↵
    1. Zamolodchikova TS,
    2. Sokolova EA,
    3. Lu D, and
    4. Sadler JE
    (2000) Activation of recombinant proenteropeptidase by duodenase. FEBS Lett 466:295–299.
    OpenUrlCrossRefPubMed
  29. ↵
    1. Zheng XL,
    2. Kitamoto Y, and
    3. Sadler JE
    (2009) Enteropeptidase, a type II transmembrane serine protease. Front Biosci (Elite Ed) 1:242–249.
    OpenUrl
View Abstract
PreviousNext
Back to top

In this issue

Journal of Pharmacology and Experimental Therapeutics: 375 (3)
Journal of Pharmacology and Experimental Therapeutics
Vol. 375, Issue 3
1 Dec 2020
  • Table of Contents
  • Table of Contents (PDF)
  • About the Cover
  • Index by author
  • Editorial Board (PDF)
  • Front Matter (PDF)
Download PDF
Article Alerts
Sign In to Email Alerts with your Email Address
Email Article

Thank you for sharing this Journal of Pharmacology and Experimental Therapeutics article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Targeting Enteropeptidase with Reversible Covalent Inhibitors To Achieve Metabolic Benefits
(Your Name) has forwarded a page to you from Journal of Pharmacology and Experimental Therapeutics
(Your Name) thought you would be interested in this article in Journal of Pharmacology and Experimental Therapeutics.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Citation Tools
Research ArticleMetabolism, Transport, and Pharmacogenetics

Sustained Inhibition of EP Leads to a Potential Therapeutic

Weimei Sun, Xuqing Zhang, Maxwell D. Cummings, Kamal Albarazanji, Jiejun Wu, Mina Wang, Richard Alexander, Bin Zhu, YueMei Zhang, James Leonard, James Lanter and James Lenhard
Journal of Pharmacology and Experimental Therapeutics December 1, 2020, 375 (3) 510-521; DOI: https://doi.org/10.1124/jpet.120.000219

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Share
Research ArticleMetabolism, Transport, and Pharmacogenetics

Sustained Inhibition of EP Leads to a Potential Therapeutic

Weimei Sun, Xuqing Zhang, Maxwell D. Cummings, Kamal Albarazanji, Jiejun Wu, Mina Wang, Richard Alexander, Bin Zhu, YueMei Zhang, James Leonard, James Lanter and James Lenhard
Journal of Pharmacology and Experimental Therapeutics December 1, 2020, 375 (3) 510-521; DOI: https://doi.org/10.1124/jpet.120.000219
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Jump to section

  • Article
    • Abstract
    • Introduction
    • Materials and Methods
    • Results
    • Discussion
    • Acknowledgments
    • Authorship Contributions
    • Footnotes
    • Abbreviations
    • References
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF

Related Articles

Cited By...

More in this TOC Section

  • OATP1A2 Transport of Statins in Human Endothelial Cells
  • Refined Prediction of Potential Kratom-Drug Interactions
  • Amyloid beta 40/42 uptake kinetics at the BBB
Show more Metabolism, Transport, and Pharmacogenetics

Similar Articles

  • Home
  • Alerts
Facebook   Twitter   LinkedIn   RSS

Navigate

  • Current Issue
  • Fast Forward by date
  • Fast Forward by section
  • Latest Articles
  • Archive
  • Search for Articles
  • Feedback
  • ASPET

More Information

  • About JPET
  • Editorial Board
  • Instructions to Authors
  • Submit a Manuscript
  • Customized Alerts
  • RSS Feeds
  • Subscriptions
  • Permissions
  • Terms & Conditions of Use

ASPET's Other Journals

  • Drug Metabolism and Disposition
  • Molecular Pharmacology
  • Pharmacological Reviews
  • Pharmacology Research & Perspectives
ISSN 1521-0103 (Online)

Copyright © 2021 by the American Society for Pharmacology and Experimental Therapeutics