Diabetes is affecting the life of millions of people. A large proportion of diabetic patients suffer from severe complications such as neuropathic pain, and current treatments for these complications have deleterious side effects. Thus, alternate therapeutic strategies are needed. Recently, the elevation of epoxy-fatty acids through inhibition of soluble epoxide hydrolase (sEH) was shown to reduce diabetic neuropathic pain in rodents. In this report, we describe a series of newly synthesized sEH inhibitors with at least 5-fold higher potency and doubled residence time inside both the human and rodent sEH enzyme than previously reported inhibitors. These inhibitors also have better physical properties and optimized pharmacokinetic profiles. The optimized inhibitor selected from this new series displayed improved efficacy of almost 10-fold in relieving pain perception in diabetic neuropathic rats as compared to the approved drug, gabapentin, and previously published sEH inhibitors. Therefore, these new sEH inhibitors could be an attractive alternative to treat diabetic neuropathy in humans.
Diabetes is affecting the life of millions of people. A large proportion of diabeticpatients suffer from severe complications such as neuropathic pain, and current treatments for these complications have deleterious side effects. Thus, alternate therapeutic strategies are needed. Recently, the elevation of epoxy-fatty acids through inhibition of soluble epoxide hydrolase (sEH) was shown to reduce diabetic neuropathic pain in rodents. In this report, we describe a series of newly synthesized sEH inhibitors with at least 5-fold higher potency and doubled residence time inside both the human and rodent sEH enzyme than previously reported inhibitors. These inhibitors also have better physical properties and optimized pharmacokinetic profiles. The optimized inhibitor selected from this new series displayed improved efficacy of almost 10-fold in relieving pain perception in diabetic neuropathicrats as compared to the approved drug, gabapentin, and previously published sEH inhibitors. Therefore, these new sEH inhibitors could be an attractive alternative to treat diabetic neuropathy in humans.
A recent survey from
the Centers for Disease Control and Prevention
indicates that diabetes affects 25.8 million people in the United
States which is 8.3% of the U.S. population.[1] Most diabeticpatients will ultimately develop kidney failure, hypertension,
and/or suffer stroke. In addition, about two-thirds of diabeticpatients
will develop peripheral neuropathy.[2,3] People suffering
from diabetic neuropathic pain experience spontaneous pain (pain sensation
in the absence of stimulation), hyperalgesia (increased pain sensation
to painful stimuli), and allodynia (pain sensation to innocuous stimuli),
which greatly affect their quality of life. Hyperglycemia has been
suggested to be the initiating cause of peripheral nerve fiber degeneration,
which results in pain.[4] However, aggressive
glycemic control can only control the progression of neuronal degeneration
but not reverse the neuropathy.[4] Current
treatments of diabetic neuropathy include tricyclic antidepressants,
anticonvulsants, selective serotonin reuptake inhibitors, and opioids,
however they often have side effects that limit their use.[5] Therefore, an alternative therapy with no or
greatly reduced side effects is still imperative for these patients
often suffering multiple comorbid conditions.Epoxy fatty acids
(EpFAs), formed by cytochrome P450 (CYP450) from
polyunsaturated fatty acids, are important lipid mediators.[6] Epoxy-eicosatrienoic acids (EETs), epoxy-eicosatetraenoic
acids (EpETEs), and epoxy-docosapentaenoic acids (EpDPEs), from arachidonic
acid, eicosapentaenoic acid, and docosahexaenoic acid, respectively,
have analgesic properties in inflammatory pain models.[7,8] Although these EpFAs are very potent lipid mediators, they are rapidly
metabolized by soluble epoxide hydrolase (sEH EC 3.3.2.10) to their
corresponding 1,2-diols and to a lesser extent by other enzymes in
vivo.[9] The in vivo biological activities
of these natural chemical mediators appear limited by their rapid
degradation. Stabilization of EpFAs through inhibition of sEH has
shown anti-inflammatory, antihypertensive, and analgesic effects.
Recent studies also indicate that sEH inhibition is analgesic in chronic
diabetic neuropathic pain in animal models. In fact, it is more efficacious
than gabapentin, a clinically approved drug for this condition.[10,11] In nonmodel species, sEH inhibitors have reduced the inflammatory
and devastating neuropathic pain in laminitis horses,[12] reduced blood pressure in forearm blood flow studies in
man,[13] and reduced neuropathic pain in
humandiabetics (www.sphaerapharma.com). Thus, sEH is a
potentially important pharmaceutical target.[6,8,9,12,14−20]Over the years, several groups have reported the synthesis
and
evaluation of sEH inhibitors with different central pharmacophores
with potency varying from micromolar to nanomolar ranges.[21−27] The 1,3-disubstituted urea is one of the more potent central pharmacophores
being used to inhibit sEH because the urea forms tight hydrogen bonds
with the active residue Asp335 and the chemistry is easily accessible.[21,23,28−30] The physical
properties of many of the most potent compounds are generally poor.
Efforts to improve physical properties including water solubility,
hydrophilicity, decreased clogP, and lowered melting point of sEH
inhibitors have generally resulted in a decrease in potency and less
desirable pharmacokinetics. These physical properties can also result
in poor absorption and inferior pharmacokinetic properties and can
demand heroic formulation.[26,30−32] Therefore, it is necessary to further optimize the structures of
the inhibitors and improve the oral bioavailability of the sEH inhibitors
carrying a 1,3-disubstituted urea as a central pharmacophores.Recent reports in drug discovery suggest that the residence time
of a drug in its target is an important parameter to predict in vivo
drug efficacy.[33] Residence time is defined
as the duration of time which the target, either enzyme or receptor,
is occupied by the ligand.[33] The traditional
IC50 and Ki or Kd is determined in a closed in vitro system in which the
concentrations of the ligand and the target are constant. However,
an in vivo system is an open system, thus the target is exposed to
a varying concentration of ligand after dosing because of circulation,
metabolism, and excretion. Therefore, drug efficacy is no longer correlated
with the in vitro potency (IC50 or Ki) that is determined in a closed system but rather depends
on the duration the target is occupied by the ligand. This residence
time can be calculated from the reciprocal of the dissociation rate
constant (koff) of the target–ligand
complex.In this report, inhibitors with improved potency have
been designed
and synthesized based on the holo-structure of the
recombinant humansEH with published sEH inhibitor: TPPU (UC1770)[23] (Figure 1B). Their residence times have been determined and their in vivo
efficacies have been tested in a diabetic neuropathic pain model.
Figure 1
(A) The
general scaffold of sEH inhibitors used in this report.
(B) The structure of inhibitor 18 (TPPU/UC1770).
(C) The name structures of previously published inhibitors.
(A) The
general scaffold of sEH inhibitors used in this report.
(B) The structure of inhibitor 18 (TPPU/UC1770).
(C) The name structures of previously published inhibitors.
Results and Discussion
Inhibitor Design and Synthesis
Piperidyl-urea inhibitors
of sEH were first described in 2006 and have the general scaffold
shown in Figure 1A.[23,32] It has been demonstrated that large amide substituents at R2 improve the potency of the sEH inhibitors (Figure 1A).[23] However, sEH inhibitors
with large amide substituents at R2 (Figure 1A) exhibit poor pharmacokinetic profiles in dogs.[34] To further improve both the potency and pharmacokinetic
profiles of the sEH inhibitors, we started with a published inhibitor
having a small substituent at R2: TPPU (inhibitor 18) (Figure 1B), which resulted in
good potency and a good pharmacokinetic profile. Because of the favorable
properties of inhibitor 18, an X-ray structure of humansEH with inhibitor 18 was obtained and used to predict
structural modifications of inhibitors leading to improved potency
and properties. The X-ray structure indicates that there is a small
secondary binding site (valley) next to the α-carbon of the
amide of inhibitor 18 (Figure 2A,B), which can provide additional binding possibilities. We therefore
hypothesized that the potency of the inhibitors could be significantly
improved by incorporating a small hydrophobic substituent at the α-carbon
of the amide (Figure 2B). In addition, the
4-trifluoromethoxyphenyl– group on the urea of inhibitor 18 fits closely to the “right side” of the binding
pocket with only a little room for additional binding (Figure 2C,D). Thus, only substituents which are similar
to the size of the 4-trifluoromethoxyphenyl– group were used
for inhibitor optimization. On the basis of these findings, a series
of 30 piperidyl-ureas was made using four different methods described
previously.[23] These methods with slight
modifications are summarized in Scheme 1 (detailed
synthetic procedures are provided in the Supporting
Information).
Figure 2
(A) Holo-crystal structure of human sEH (green)
with inhibitor 18 (TPPU) (cyan) (PDB code: 4OD0). (B) The left side
of the tunnel of human sEH with inhibitor 18 (cyan). The arrow indicated the valley of the left side of the tunnel
for potential additional binding for new inhibitors. (C,D) The right
binding pocket of human sEH with UC1770 from the
view of the front and back (cyan). The graphics were prepared by the
PyMOL Molecular Graphics System, version 1.3 edu, Schrodinger, LCC.
Scheme 1
Synthetic Schemes for sEH Inhibitors Synthesis
(A) Holo-crystal structure of humansEH (green)
with inhibitor 18 (TPPU) (cyan) (PDB code: 4OD0). (B) The left side
of the tunnel of humansEH with inhibitor 18 (cyan). The arrow indicated the valley of the left side of the tunnel
for potential additional binding for new inhibitors. (C,D) The right
binding pocket of humansEH with UC1770 from the
view of the front and back (cyan). The graphics were prepared by the
PyMOL Molecular Graphics System, version 1.3 edu, Schrodinger, LCC.
Optimization of the Potency
(Ki)
of sEH Inhibitors
A new series of sEH inhibitors was synthesized
with various substituents at R2 on the amide (Table 1) while maintaining R1 as a 4-trifluoromethylphenyl–
group. In general, as the size of the substituent increases, the potency
of the inhibitors against humansEH increases. As we hypothesized,
addition of a methyl– group on the α-carbon of the amide
(Figure 1A) greatly improves the potency by
almost 5-fold (inhibitors 2 vs 4 and inhibitors 3 vs 6). Interestingly, replacing the isopropyl–
(4) with cyclopropyl– group (5) also
results in a slight increase in potency, an increase in water solubility,
and a decrease in melting point. The elongation of the aliphatic chain
by one carbon (inhibitors 2 vs 3) slightly
enhances (30%) the potency (Table 1). To investigate
whether the enhanced potency of inhibitor 6 is stereospecific,
we tested the corresponding S-isomer (inhibitor 7), which is 2-fold more potent than its racemic mixture (Table 1), suggesting that the R-isomer
is less active. The related inhibitor 8 is less potent
than inhibitors 6 and 7, which indicates
that the addition of methyl– group not only is stereospecific
but is also regiospecific. An X-ray structure of humansEH with inhibitor 4 shows that the methyl– group on the α-carbon
of the amide of inhibitor 4 effectively fits into the
additional binding site on the right side of the binding pocket (Figure 3). This additional binding site was predicted from
the crystal structure of humansEH with inhibitor 18 (Figure 2B).
Table 1
Physical Properties
and Potency of
sEH Inhibitors against Human sEH (Modification of R2)e
Solubility was
measured with sodium
phosphate buffer (0.1 M, pH 7.4) according to the method described
by Tsai et al. and described in detail in Supporting
Information.[34]
N.D. means “not determined”.
elogP was determined by HPLC method
calibrated with elogP of six selected inhibitors determined by the
shake-flask method (Supporting Information, Figure
S3).
Ki was
determined by FRET-based displacement assay described by Lee et al.[38] The results are the average of duplicates with
± SEM.
Abbreviation:
elogP stands for experimental
Log P.
Figure 3
Overlay structures of human sEH with
inhibitor 18 (cyan) and inhibitor 4 (orange).
This figure
suggests that the design principle is valid and the methyl–
group at α-position of the amide provides extra binding toward
the valley of the left binding pocket. The graphics were prepared
by the PyMOL Molecular Graphics System, version 1.3 edu, Schrodinger,
LCC.
Overlay structures of humansEH with
inhibitor 18 (cyan) and inhibitor 4 (orange).
This figure
suggests that the design principle is valid and the methyl–
group at α-position of the amide provides extra binding toward
the valley of the left binding pocket. The graphics were prepared
by the PyMOL Molecular Graphics System, version 1.3 edu, Schrodinger,
LCC.Solubility was
measured with sodium
phosphate buffer (0.1 M, pH 7.4) according to the method described
by Tsai et al. and described in detail in Supporting
Information.[34]N.D. means “not determined”.elogP was determined by HPLC method
calibrated with elogP of six selected inhibitors determined by the
shake-flask method (Supporting Information, Figure
S3).Ki was
determined by FRET-based displacement assay described by Lee et al.[38] The results are the average of duplicates with
± SEM.Abbreviation:
elogP stands for experimental
Log P.Although the addition of an alkyl group enhances the potency of
this series of sEH inhibitors, their metabolic stability will likely
decrease because the added alkyl group at the α-position is
expected to be metabolized faster. This also increases the lipophilicity
of the inhibitors, which may increase their affinity to CYP450 enzymes
that are responsible for drug oxidation.[35] Because fluoride replacement is known to decrease metabolism,[36] we synthesized inhibitor 9,[36] and it shows potency similar to inhibitor 8. The size of the CF3– group is similar
to the size of the isopropyl–, group and the potency of inhibitor 9 is maintained. This could provide an approach to increase
the stability of the inhibitors without hampering the potency of the
inhibitors.We then further optimized the potency of sEH inhibitors
by modifying
the substituent at R1 while maintaining the 2-methyl butanoyl
group at R2 (Table 2). The binding
pocket is very hydrophobic (Supporting Information,
Figure S2) with limited space (Figure 2C,D). Therefore, the size of R1 was modified slightly
in order to test if increasing the hydrophobic surface of the inhibitors
could enhance their potency. Consequently, inhibitors 11 and 12 were synthesized (Table 2). The result indicates that there is steric hindrance within the
right binding pocket (Figure 2C,D). Thus, the
potency of the inhibitor dropped when the size of the substituent
at R1 (Figure 1A) was increased
from a 4-iso-propylphenyl group to a 4-tert-butylphenyl- group (inhibitors 11 and 12) (Table 2). It is interesting that the inhibitor 6, which has CF3– group in place of an iso-propyl– and tert-butyl–
group,[37] is more potent as an inhibitor
when compared to inhibitors 11 and 12. To
investigate whether the enhanced potency of inhibitor 6 is due to fluorine induced interactions or better occupancy of the
binding pocket, the 4-t-butylphenyl– group
of inhibitor 12 was replaced by its isostere: a 4-heptafluoro-iso-propyl-phenyl group at R1 of inhibitor 13.[37] This compound is 20-fold
more potent than inhibitor 6, which indicates that the
fluorine induced interactions in the left binding pocket are very
strong (Table 2). These data suggest that the
left side of the binding pocket is likely fluorophilic. Interestingly,
replacement of the 4-trifluoromethylphenyl– at R1 of inhibitor 6 by a 4-trifluoromethoxyphenyl–
group (inhibitor 14) did not alter potency, but the replacement
improved solubility by 10-fold (Table 1). Furthermore,
inhibitors with a cycloalkyl– group at R1 were synthesized,
and the addition of carbon atoms (inhibitor 15 to 17) enhances the binding toward humansEH by more than 4 times
(Table 2) and placement of cycloalkyl–
group at R1 greatly enhances their solubility.
Table 2
Physical Properties and Potency of
sEH Inhibitors against Human sEH (Modification of R1)e
Solubility was
measured with sodium
phosphate buffer (0.1 M, pH 7.4) according to the method described
by Tsai et al. and described in detail in Supporting
Information.[34]
N.D. means “not determined”.
elogP was determined by HPLC method
calibrated with elogP of six selected inhibitors determined by the
shake-flask method (Supporting Information, Figure
S3).
Ki was
determined by FRET-based displacement assay described by Lee et al.[38] The results are the average of duplicates with
± SEM.
Abbreviation:
elogP stands for experimental
Log P.
Solubility was
measured with sodium
phosphate buffer (0.1 M, pH 7.4) according to the method described
by Tsai et al. and described in detail in Supporting
Information.[34]N.D. means “not determined”.elogP was determined by HPLC method
calibrated with elogP of six selected inhibitors determined by the
shake-flask method (Supporting Information, Figure
S3).Ki was
determined by FRET-based displacement assay described by Lee et al.[38] The results are the average of duplicates with
± SEM.Abbreviation:
elogP stands for experimental
Log P.It has been reported that inhibitors with 4-trifluoromethoxyphenyl–
at R1 are more potent than inhibitors with 4-trifluoromethylphenyl–
at R1 and have better physical properties.[23] To investigate whether such substitution could enhance
the potency of the inhibitors and improve physical properties in general,[38] a series of inhibitors with 4-trifluoromethylphenyl
group at R1 was replaced with a 4-trifluoromethoxyphenyl–
group (Figure 1A). The results suggest that
substitution of 4-trifluoromethylphenyl– with 4-trifluoromethoxyphenyl–
enhances the potency of certain inhibitors (2 vs 18 and 4 vs 19 in Tables 1 and 3), but there is no
substantial improvement with inhibitors with the 2-methyl butanoyl
group at R2 (6 vs 14 and 7 vs 21 in Table 1–3) or the cyclopropyl group at R2 (5 Vs 20). This may be due to the possibility
of different binding orientations than the one shown in Figure 2A and in previous crystallography.[39] Therefore, two inhibitors which are structurally similar
to 21 but with urea substituted by the amide (22 and 23) were synthesized. In general, structure–activity
relationships developed with the urea central pharmacophore are predictive
of enzyme inhibition with amide and carbamate. Both inhibitors 22 and 23 are less potent than their corresponding
urea (21) (Table 3). Compound 22, which has amidenitrogen (N″H) on the piperidine
side (Figure 1A), is far less potent than the
amide 23, which has an amidenitrogen (N′H) on
the R1 side (Figure 1A) (Table 2). It indicates that the urea N′H–Asp(O)
interaction is stronger than the urea N″H–Asp(O) interaction
by 5.7 kJ mol–1, calculated based on the Ki of both inhibitors. This result is consistent
with the crystal structure where the distance between urea N′H
and Asp(O) is shorter than the distance between urea N″H and
Asp(O) hydrogen bonding (Figure 2A). In some
cases, the placement of 4-trifluoromethoxyphenyl– at R1 (Figure 1A) enhances the potency of
sEH inhibitors, but in all cases, this substitution decreases the
melting point and improves the solubility of the inhibitors (Tables 1–3). These improvements
of physical properties ease the drug formulation process.
Table 3
Physical Properties and Potency of
sEH Inhibitors against Human sEH (Modification of R1)e
Solubility was
measured with sodium
phosphate buffer (0.1 M, pH 7.4) according to the method described
by Tsai et al. and described in detail in Supporting
Information.[34]
N.D. means “not determined”.
elogP was determined by HPLC method
calibrated with elogP of six selected inhibitors determined by the
shake-flask method (Supporting Information, Figure
S3).
Ki was
determined by FRET-based displacement assay described by Lee et al.[38] The results are the average of duplicates with
± SEM.
Abbreviation:
elogP stands for experimental
Log P.
Solubility was
measured with sodium
phosphate buffer (0.1 M, pH 7.4) according to the method described
by Tsai et al. and described in detail in Supporting
Information.[34]N.D. means “not determined”.elogP was determined by HPLC method
calibrated with elogP of six selected inhibitors determined by the
shake-flask method (Supporting Information, Figure
S3).Ki was
determined by FRET-based displacement assay described by Lee et al.[38] The results are the average of duplicates with
± SEM.Abbreviation:
elogP stands for experimental
Log P.It was reported that replacement of amide with sulfonamide at R2 was shown to enhance the potency of the inhibitors. However,
such observations were based on very few comparisons.[23,34] Therefore, a new series of inhibitors with a different sulfonamide
at R2 was synthesized for a more detailed study (Table 4). In general, the potency of the inhibitors increases
with the size of R2 (7-fold better potency from a methyl 25 to a butyl 31). However, unlike the inhibitor
with an amide at R2 (Table 1), the
placement of an iso-propylsulfonamide at R2 (29 vs 26) does not significantly enhance
the potency against humansEH. This is probably because the sulfonamide
exists as a tetrahedron while amide is trigonal planar and the S(O)2–C bond is at least 0.28 Å longer than the C(O)–C
bond.[40] These data further indicate that
the enhanced potency of the methyl substitution at the α-position
of the amide is structurally specific. Overall, unlike previously
reported, the inhibitors with a sulfonamide at R2 are less
potent than the inhibitors with amides at R2. In addition,
the physical properties (solubility and melting point) of inhibitors
with a sulfonamide at R2 are generally poor as compared
to inhibitors with amides at R2.
Table 4
Physical
Properties and Potency of
sEH Inhibitors against Human sEH (Modification of R1 and
R2 with sulfonamide at R2)e
Solubility was measured with sodium
phosphate buffer (0.1 M, pH 7.4) according to the method described
by Tsai et al. and described in detail in Supporting
Information.[34]
N.D. means “not determined”.
elogP was determined by HPLC method
calibrated with elogP of six selected inhibitors determined by the
shake-flask method (Supporting Information, Figure
S3).
Ki was
determined by FRET-based displacement assay described by Lee et al.[38] The results are the average of duplicates with
± SEM.
Abbreviation:
Sol. stands for Solubility;
elogP stands for experimental Log P.
Solubility was measured with sodium
phosphate buffer (0.1 M, pH 7.4) according to the method described
by Tsai et al. and described in detail in Supporting
Information.[34]N.D. means “not determined”.elogP was determined by HPLC method
calibrated with elogP of six selected inhibitors determined by the
shake-flask method (Supporting Information, Figure
S3).Ki was
determined by FRET-based displacement assay described by Lee et al.[38] The results are the average of duplicates with
± SEM.Abbreviation:
Sol. stands for Solubility;
elogP stands for experimental Log P.Here, we have identified several
structural changes that can significant
enhance the potency and improve the physical properties of the inhibitors.
We have also demonstrated that the sulfonamide at R2 is
less attractive than the corresponding amide.
New sEH Inhibitors Have
Improved Dissociation Rate Constants
(koff) against Human sEH
Recent
studies have suggested that koff, a kinetic
parameter on enzyme inhibition, from the enzyme is a better indicator
for in vivo potency than Ki.[33] This is because inhibitors are only effective
in blocking catalysis when the target proteins are occupied by the
inhibitors. The koff can provide more
detailed information, about the duration of time the inhibitors are
bound to the target enzyme (target occupancy) than the Ki, an equilibrium parameter on enzyme inhibition, and
ultimately this translates into in vivo efficacy.[40] Therefore, a small set of new potent inhibitors together
with several potent published inhibitors (Figure 1C) were selected to determine their koff against sEH using a recently developed FRET-based assay.[38]The result indicates that the koff of inhibitors decreases with the size of
R2 increased (Table 5 and Figure 1C, inhibitor 4 < 7; TPAU < 18 (TPPU) < 19 < 21; and TUPS < 32). However, the koff does not substantially
change when the R1 was varied among a 4-trifluoromethylphenyl–
group (inhibitor 6), a 4-trifluoromethoxyphenyl–
group (inhibitor 14), and a 4-isopropylphenyl–
group (inhibitor 11) (Table 5).
The 4-heptafluoroisopropylphenyl– group at R1 (inhibitor 13 and inhibitor 24) significantly decreased
the koff of the inhibitors. This further
supports our hypothesis that the pocket is fluorophilic and the added
fluorines can induce several interactions with the nearby residues
within the binding pocket. The koffs of
new inhibitors are slower overall, indicating a longer residence time
in the target and have at least a 2-fold slower off rate than any
of the previously published inhibitors (Table 5 and Figure 1C, APAU, TPAU, 18 (TPPU), TUPS, t-TUCB).
Table 5
Study in Vitro Target
Occupancy (koff) of Selected sEH Inhibitors
against Human
sEH
entry
inhibitor
koff (× 10–4 s–1)a
t1/2 (min)b
entry
inhibitor
koff (× 10–4 s–1)a
t1/2 (min)b
1
APAU
19.2 ± 0.70
6.03 ± 0.2
11
17
3.51 ± 0.20
33.0 ± 1.9
2
18 (TPPU)
10.5 ± 0.20
11.0 ± 0.2
12
19
6.14 ± 0.18
18.8 ± 0.6
3
4
6.57 ± 0.30
17.6 ± 0.8
13
21
5.05 ± 0.02
22.9 ± 0.1
4
5
7.91 ± 0.31
14.6 ± 0.6
14
24
4.39 ± 0.43
26.5 ± 2.6
5
6
5.76 ± 0.26
20.1 ± 0.9
15
25
23.1 ± 1.1
5.02 ± 0.3
6
7
5.19 ± 0.09
22.3 ± 0.4
16
31
10.3 ± 0.1
11.0 ± 0.1
7
9
4.75 ± 0.11
24.3 ± 0.6
17
32
8.90 ± 0.35
13.0 ± 0.5
8
11
5.79 ± 0.43
20.0 ± 1.5
18
33 (TPAU)
26.4 ± 2.30
4.40 ± 0.4
9
13
3.13 ± 0.06
37.0 ± 0.7
19
TUPS
20.0 ± 0.64
5.79 ± 0.2
10
14
5.39 ± 0.39
21.5 ± 1.6
20
t-TUCB
7.19 ± 0.36
16.1 ± 0.8
koff was determined by FRET-based displacement assay described
by Lee
et al.[38] and described in detail in Supporting Information. Brieftly, a preincubated
human sEH–inhibitor complex (8 μM) was diluted by ×40
times by fluorescent reporter–APCU (2 μM,
0.1 M sodium phosphate, pH 7.4). The fluorescent enhancement (λexcit = 280 nm, λemit = 450) was measured
over time (5100 s). The results are the average of triplicates with
± SD
t1/2 = ln(2)/koff, which describes
the half-life
of enzyme–inhibitor complex.
koff was determined by FRET-based displacement assay described
by Lee
et al.[38] and described in detail in Supporting Information. Brieftly, a preincubated
humansEH–inhibitor complex (8 μM) was diluted by ×40
times by fluorescent reporter–APCU (2 μM,
0.1 M sodium phosphate, pH 7.4). The fluorescent enhancement (λexcit = 280 nm, λemit = 450) was measured
over time (5100 s). The results are the average of triplicates with
± SDt1/2 = ln(2)/koff, which describes
the half-life
of enzyme–inhibitor complex.To investigate relationship of koff with Ki and other physical
properties
of the inhibitors, correlation graphs of koff with these parameters were plotted (Supporting
Information, Figure S1A). The results show no correlation (R2 = 0.21) between the number of non-hydrogen
atoms and koff. However, there is a trend
showing (R2 = 0.52) that an increase in
elogP results in a decrease of koff (Supporting Information, Figure S1B). This may
due to the fact that the binding pocket of the humansEH is hydrophobic
(Supporting Information, Figure S2). An
increase of hydrophobicity of the ligand increases the lipophilic
interactions with the binding pocket. Therefore, it requires higher
activation energy to break up the interactions between protein and
the bound inhibitors. In addition, a correlation between Ki and koff was plotted (Supporting Information, Figures S1C, S1D). The
plot indicates that there is a good correlation between Ki and koff (R2 = 0.88) over a wide range of potencies (from 0 to 20
nM). However, when we focused on a narrower range of potencies (from
0 to 1.4 nM), the correlation is moderate (R2 = 0.44). Because Ki is a ratio
of koff over kon, Ki should be inversely proportional
to koff. Therefore, the poor correlation
over a moderate range Ki values suggests
that the differences are due to the kon. These results indicate that it is possible to specifically modulate koff without greatly affecting Ki.
An Improved Pharmacokinetic Profile of New
Series of sEH Inhibitors
To investigate the oral bioavailability
of sEH inhibitors, we determined
the pharmacokinetic profiles post oral administration. Cassette dosing
was used as a screening tool to select the compounds for more detailed
study to allocate limited resources to the most promising compounds.
There are cautions with cassette dosing including changes in pharmacokinetic
behavior due to competition for xenobiotic metabolism. The high potency
and thus low doses used in this study make this artifact less likely.
A comparison of pharmacokinetic profiles of sEH inhibitors between
the cassette dose and single dose carried out in nonhuman primate.
The results show that there was no statistically significant difference
in the pharmacokinetic behaviors between cassette and individual dosing
of sEH inhibitors; therefore, cassette dosing was predictive for the
compounds reported here.[41]The pharmacokinetic
profiles of this series of selected potent compounds also determined
the effect of the position or addition of aliphatic carbons on the
inhibitors. In general, the results indicate that the sEH inhibitors
with a 4-trifluoromethylphenyl– group or a 4-trifluoromethoxyphenyl–
group at R1 (Figure 1A) have good
drug exposure levels based on the area under the curve of the pharmacokinetic
kinetic profile (PK-AUC) after oral administration (Table 6). As we replaced the 4-trifluoromethylphenyl–
group or 4-trifluoromethoxyphenyl– group at R1 with
a 4-isopropylphenyl– group (inhibitors 6 or 14 vs 11), the potency of inhibitor 14 decreases (Table 2). As we anticipated, the
pharmacokinetic T1/2 (PK-T1/2) and oral drug exposure level estimated by PK-AUC
also decreases (Table 6). Further replacement
of the phenyl group at R1 with a cycloalkyl– group
(inhibitors 16 and 17) greatly decreases
the oral drug exposure level (Table 6, Supporting Information, Figures S6, S8, and S9). It is likely that each addition of an aliphatic carbon renders
the molecule more susceptible to metabolism by CYP450 enzymes.[42] Such a phenomenon was also observed as the alkyl
chain length varied at R2. As we hypothesized, the PK-T1/2 of the inhibitors also decrease when the
alkyl chain length at R2 of the amide and sulfonamide series
increases (Table 6 for amide inhibitors, 18, 19, and 21; for sulfonamide
inhibitors, 25 and TUPS (Figure 1C) vs 31 and 32). The
result also show that the PK-AUC estimated drug exposure levels of
the sulfonamide inhibitors 31 and 32 are,
in general, worse than the amide inhibitors (Table 6).
Table 6
Pharmacokinetic Parameters of Selected
sEH Inhibitors after Oral Dosing on Micea
inhibitor
PK-AUCa (nM·h)b
Cmax (nM)c
T1/2b (h)
Tmax(h)c
APAU(10)
183
30
3
1
TUPS(18)
5300
370
15
3.3
18 (TPPU)
10650
495
12.1
8
4
19650
1160
6
4
5
24630
960
17
6
6
3295
320
6
0.5
7
4940
435
5
1
9
8845
330
15.4
8
11
200
90
2
0.5
13
1600
135
4.5
2.5
14
3300
235
6
2
16
75
3
8.4
0.5
17
37
5
12
0.5
19
19500
900
9
5
21
2530
195
7.7
4
24
9600
365
18.9
8
25
4675
130
22
8
30
4270
100
25
8
31
245
30
7.4
1
32
220
65
2
1.2
The mice (n = 4)
were treated by oral dosing with a cassette of 3 to 5 compounds (0.3
mg/kg per each compounds dissolved in 20% PEG400 in oleic acid rich
triglycerides).
The pharmacokinetic
profiles of
the inhibitors were calculated by Winonlin based on the model of one
compartmental analysis.
The pharmacokinetic profiles of
the inhibitors were calculated by Winonlin based on the model of noncompartmental
analysis.
The mice (n = 4)
were treated by oral dosing with a cassette of 3 to 5 compounds (0.3
mg/kg per each compounds dissolved in 20% PEG400 in oleic acid rich
triglycerides).The pharmacokinetic
profiles of
the inhibitors were calculated by Winonlin based on the model of one
compartmental analysis.The pharmacokinetic profiles of
the inhibitors were calculated by Winonlin based on the model of noncompartmental
analysis.We then synthesized
inhibitors 9 and 24 in order to study whether
the addition of fluorine could enhance
the stability of our inhibitors as has been previously suggested.[36] The replacement of a terminal methyl group with
a CF3– group at the R2 of the amide inhibitors
(inhibitors 9 vs TPPU and 24 vs 13) not only improves the PK-T1/2 of the inhibitors but also increases their oral drug exposure
levels (Table 6). A similar result was obtained
with the sulfonamide inhibitor. Inhibitor 30 shows at
least a 20-fold better bioavailability based on oral PK-AUC with a
substantially longer PK-T1/2 as compared
to inhibitors 31 and 32 (Table 6). These data support the hypothesis that the CF3– substitution effectively blocks metabolism at R2. However, addition of fluorine to the fluorinated substituents
could not further enhance the stability of the inhibitors. When the
inhibitors with 4-trifluoromethyl– groups at R1 were
substituted by a heptafluoro-isopropyl– group (inhibitor 6 vs 13), both the PK-T1/2 and bioavailability decreased (Table 6). This is probably due to the fact that the increase of elogP enhances
the affinity toward xenobiotic metabolizing enzymes and cell membranes.
These results suggest that in order to enhance the stability of the
inhibitors, addition of fluorine should be carefully positioned.[35,43]In general, inhibitors with 4-trifluoromethylphenyl–,
4-trifluoromethoxyphenyl–,
and 4-(heptafluoro-iso-propyl)phenyl– groups
at R1 with amide substituents at R2 showed good
pharmacokinetic profiles. Inhibitors 16, 17, 31, and 32, although very potent, are
not optimal candidates for testing in animal chronic disease models
because of poor oral drug exposure levels. Our data also suggest that
these inhibitors (16, 17, 31, and 32) may not be top candidates to be optimized
for human drugs. The inhibitors with a piperidyl moiety attached at
N2 of the urea show greatly improved oral pharmacokinetic
profiles compared to inhibitors with a cyclohexyl moiety attached
at N2 of the urea also (Supporting
Information, Figure S6).[30] Overall,
our data suggested that several newly synthesized sEH inhibitors have
an improved and optimized pharmacokinetic profile for diseases requiring
chronic treatment.
New Series of sEH Inhibitors Show No Significant
Nonspecific
Binding toward Other Pharmacologically Important Proteins
Several of the newly synthesized sEH inhibitors selected for potency
and good oral drug exposure levels were tested against a set of proteins
to evaluate potential off-target side effects and nonspecific binding.
Inhibitor 24 was excluded because of its relatively high
elogP (Table 3), which was suspected to lead
to high nonspecific binding. High lipophilicity is avoided when possible
in medicinal chemistry because it is associated with poor pharmacokinetic
behavior and nontargeted related side effects. Although not conclusive,
the pharmacokinetic studies indicated that the 4-(heptafluoro-iso-propyl-)phenyl group likely enhances the affinity of
inhibitors toward xenobiotic metabolizing enzymes and cell membranes,
potentially leading to the poor oral drug exposure of inhibitor 13 (Table 6).It was reported
that very high plasma protein binding can not only hamper the efficacy
of drugs or signaling molecules but also can alter their pharmacokinetic
profile.[44,45] However, modulating plasma protein binding
can help drug solubilization and drug distribution.[46] Therefore, the plasma protein binding of the selected inhibitors
were measured with an assay carried out with a rapid equilibrium dialysis
device based on manufacturer’s protocol (Table 7). These new inhibitors demonstrated a moderate level of plasma
protein binding ranging from 85 to 96%. This level of plasma protein
binding is unlikely to alter their bioavailability. Given the inhibitors
had Kis approaching the subnanomolar range,
the plasma proteins could be considered as carrier proteins that facilitate
the their distribution.
Table 7
Testing for Nonspecific
Protein Binding
with Selected sEH Inhibitors
inhibitor
hERG channel
inhibition at 50 μM (%)a
human plasma
protein binding at 1 μM (%)a
human CYP
2J remaining activity at 10 μM (%)b,d
human CYP2C
remaining activity at 10 μM (%)b,d
APAU
4.5 ± 0.5
50.0 ± 2.0
140.0 ± 5.0
104.0 ± 1.0
18 (TPPU)
26.0 ± 1.0
79.0 ± 1.0
91.9 ± 2.2
118.0 ± 2.4
4
16.5 ± 1.5
85.0 ± 0.5
81.6 ± 3.0
110.0 ± 2.0
5
41.5 ± 1.5
95.1 ± 0.5
N.D.c
N.D.c
7
N.D.c
96.3 ± 0.2
82.6 ± 3.1
107.0 ± 3.0
19
32.5 ± 1.5
89.0 ± 2.5
73.4 ± 0.3
103.0 ± 1.0
21
N.D.c
92.4 ± 0.1
60.6 ± 8.2
108.0 ± 1.0
The result is the
average of duplicates
with standard error shown.
The result is the average of triplicates
with standard derivation shown.
N.D. means “not determined”.
The experiments were conducted according
to the Graves et al. procedure.[50]
The result is the
average of duplicates
with standard error shown.The result is the average of triplicates
with standard derivation shown.N.D. means “not determined”.The experiments were conducted according
to the Graves et al. procedure.[50]Previously, sEH inhibitors with
similar structures failed to show
significant inhibition of the CYP450s which are highly involved in
xenobiotic metabolism.[47] In this study,
we examined the inhibition of CYP2C and CYP2J2 because CYP2C is an
important drug metabolizing enzyme but principally because both enzymes
are implicated in the synthesis of epoxy fatty acids.[48−50] The results indicated that only minor inhibition of both CYP450
enzymes is observed at 10 μM of the sEH inhibitors (Table 7). These data suggested that the sEH inhibitors
do not affect the biosynthesis of EpFAs by inhibition of CYP2C and
CYP2J2. This in turn suggests that the in vivo efficacy of the inhibitors
is not due to inhibition upstream in the CYP450 pathway.Inhibition
of the hERG channel is a very important toxicology screen
due to demonstrated association with cardiotoxicity. Most of the tested
inhibitors show very minor inhibition on hERG at 50 μM except
for inhibitors 5 and 19 (Table 7). However, even inhibitors 5 and 19, which have 0.49 and 0.31 nM Ki against recombinant humansEH, have more than an 10000-fold selectivity
over the hERG channel. Therefore, these selected inhibitors are considered
not to present a risk with hERG inhibition.In summary, we have
demonstrated that our sEH inhibitors having Ki value in the lower subnanomolar are unlikely
to induce unanticipated side effects due to nonspecific binding to
other pharmacologically important proteins. However, this does not
completely rule out idiopathic off-target effects.
sEH Inhibitors
Are Effective against an Animal Model of Type
1 Diabetic Neuropathic Pain
Several new potent sEH inhibitors
with good pharmacokinetic profiles and no significant nonspecific
binding on human proteins and enzymes were selected for further in
vitro testing prior to in vivo studies. Because the in vivo studies
were to be conducted in rat, we determined the potency (Ki and koff) against the recombinant
ratsEH (Table 8). The data suggest that the
inhibitors have a slightly different structure–activity relationship
(SAR) on the recombinant ratsEH compared to humansEH. The inhibitors
with a 4-trifluoromethoxyphenyl– group at R1 are
less potent than the inhibitors with a 4-trifluoromethylphenyl–
group at R1 by at least 10-fold (Table 8). These results are opposite to the SAR obtained from the
human recombinant sEH (Table 1 and 3). In addition, there is not a clear SAR based on
the koff of the inhibitors (Table 8). However, the results indicate that the new inhibitors
are more potent, both in IC50 and koff, than the previously published inhibitors 18 (TPPU) and TPAU (Figure 1C and Table 8). Additionally, the pharmacokinetic
profiles in rat were obtained for this new set of inhibitors (Supporting Information, Figure S11). All of the
selected inhibitors show good oral bioavailability (PK-AUC) except
inhibitor 21, which demonstrated at least 2× lower
blood concentration than the other inhibitors (Table 9). Elimination is an important parameter of overall pharmacokinetics,
and it has been previously reported that inhibitors with short elimination
PK-T1/2 are not well suited to use for
chronic treatment.[51,52] Several of the new inhibitors
demonstrate good pharmacokinetic profiles with moderate elimination
PK-T1/2 (≥5 h) in both mice and
rats. Therefore, these inhibitors were selected for further in vivo
efficacy test on rat.
Table 8
In Vitro Protency
of Selected sEH
Inhibitors against Rat sEH
inhibitor
IC50 (nM)a,b
koff (× 10–4 s–1)a,c
t1/2 (min)a,d
18 (TPPU)
29.1 ± 4.5
8.52 ± 0.47
13.6 ± 0.8
TPAU
79 ± 2.5
9.70 ± 0.14
11.9 ± 0.18
4
1.8 ± 0.1
3.26 ± 0.72
36.6 ± 8.2
7
<1.25
5.37 ± 0.45
21.6 ± 1.8
19
13.3 ± 1.5
6.87 ± 0.11
16.8 ± 0.3
21
20 ± 0.1
4.89 ± 0.28
23.7 ± 1.4
The result is the
average of triplicates
with standard derivation (SD) shown.
IC50 is determined by
radiometric assay using [3H]-t-DPPO (50 μM) as a
substrate and rat sEH (2.5 nM) incubated at 30 °C for 10 min
koff was determined by FRET-based displacement assay described by Lee
et al.[38] Briefly, a preincubated rat sEH–inhibitor
complex (8 μM) was diluted by ×40 times by fluorescent
reporter–APCU (2 μM, 0.1 M sodium phosphate,
pH 7.4). The fluorescent enhancement (λexcit = 280
nm, λemit = 450) was measured over time (5100 s).
The results are the average of triplicates with ± SD.
t1/2 = ln(2)/koff, which describes the half-life
of enzyme–inhibitor complex.
Table 9
Pharmacokinetic Parameters of Selected
sEH Inhibitors Followed by Oral Dosing on Rata
structure
PK-AUCa (nM·h)b
Cmax (nM)b
T1/2b (h)
Tmax(h)b
4
7304
510
5.4
4
7
5716
790
4.9
1
19
12000
890
3.5
1
21
2590
297
4.2
4
The rat (n = 3
or 4) were treated by oral dosing with a cassette of four compounds
(0.3 mg/kg per each compound dissolved in 20% PEG400 in oleic acid
rich triglycerides).
The
pharmacokinetic profiles of
the inhibitors were calculated by Winonlin based on the best fit model
of one compartmental analysis.
The result is the
average of triplicates
with standard derivation (SD) shown.IC50 is determined by
radiometric assay using [3H]-t-DPPO (50 μM) as a
substrate and ratsEH (2.5 nM) incubated at 30 °C for 10 minkoff was determined by FRET-based displacement assay described by Lee
et al.[38] Briefly, a preincubated ratsEH–inhibitor
complex (8 μM) was diluted by ×40 times by fluorescent
reporter–APCU (2 μM, 0.1 M sodium phosphate,
pH 7.4). The fluorescent enhancement (λexcit = 280
nm, λemit = 450) was measured over time (5100 s).
The results are the average of triplicates with ± SD.t1/2 = ln(2)/koff, which describes the half-life
of enzyme–inhibitor complex.The rat (n = 3
or 4) were treated by oral dosing with a cassette of four compounds
(0.3 mg/kg per each compound dissolved in 20% PEG400 in oleic acid
rich triglycerides).The
pharmacokinetic profiles of
the inhibitors were calculated by Winonlin based on the best fit model
of one compartmental analysis.Because inhibitor 7 has the best in vitro potency
(IC50 and koff) against ratsEH with a good pharmacokinetic profile (Table 8 and 9), it was tested in a model of type
1 diabetic neuropathic pain. At as low as 0.1 mg/kg, inhibitor 7 effectively increases mechanical withdrawal thresholds (reduces
pain) in neuropathicrats (Figure 4A). Importantly,
this reduction in pain behavior increases dose dependently from 0.1
to 0.3 mg/kg (Figure 4A,B, p < 0.05). In addition, there is a good correlation between the
blood concentration of inhibitor 7 and increased withdrawal
thresholds (Figure 4C,D).
Figure 4
(A) Inhibitor 7 (IC50 rat sEH = <1.25
nM) improves mechanical withdrawal thresholds in a model of diabetic
neuropathy. Oral dosing of 0.1 and 0.3 mg/kg dose dependently increased
MWT, indicating pain relief (n = 5, mean ± SEM).
The neuropathic baseline is normalized to 100% to show the response
to a single dose of sEH inhibitor over the time course. The response
to treatments depended on the time, but when the treatments were compared,
there was a statistically significant increase from the 0.1 to the
0.3 mg/kg dose (Mann–Whitney rank sum test, U = 454.5, n1 = n2 = 45, p = <0.001). (B) Inhibitor 7 shows dose dependent pain relief in diabetic rats. AUC describes
an area under the curve of the MWT post-oral dosing of the inhibitor
vs time. The bar chart depicts the AUC of MWTs after oral dosing of
inhibitor 7. An increased in AUC of MWT is interpreted
as an increase in pain relief. When the AUC for these doses were compared,
this relationship maintained statistical significance (Mann–Whitney
rank sum test, U = 2.00, n1 = n2 = 5, p = 0.032).
(C) Efficacy in nociceptive assays relates to blood concentration.
A plot of the nociceptive responses (n = 5, mean
± SEM) vs blood concentration (n = 4, mean ±
SEM) reveals increasing efficacy with increasing blood concentration.
(D) The efficacy of inhibitor 7 is dependent on blood
concentration. When compared, the nociceptive responses (n = 5, mean ± SEM) and the pharmacokinetic profile after oral
dosing (n = 4, mean ± SEM) of 7 follow the same trend while revealing a slight delay to in the behavioral
assay. The graphics and statistics were prepared by KaleidaGraph version
4.1 (Synergy Software, Reading, PA).
(A) Inhibitor 7 (IC50 ratsEH = <1.25
nM) improves mechanical withdrawal thresholds in a model of diabeticneuropathy. Oral dosing of 0.1 and 0.3 mg/kg dose dependently increased
MWT, indicating pain relief (n = 5, mean ± SEM).
The neuropathic baseline is normalized to 100% to show the response
to a single dose of sEH inhibitor over the time course. The response
to treatments depended on the time, but when the treatments were compared,
there was a statistically significant increase from the 0.1 to the
0.3 mg/kg dose (Mann–Whitney rank sum test, U = 454.5, n1 = n2 = 45, p = <0.001). (B) Inhibitor 7 shows dose dependent pain relief in diabeticrats. AUC describes
an area under the curve of the MWT post-oral dosing of the inhibitor
vs time. The bar chart depicts the AUC of MWTs after oral dosing of
inhibitor 7. An increased in AUC of MWT is interpreted
as an increase in pain relief. When the AUC for these doses were compared,
this relationship maintained statistical significance (Mann–Whitney
rank sum test, U = 2.00, n1 = n2 = 5, p = 0.032).
(C) Efficacy in nociceptive assays relates to blood concentration.
A plot of the nociceptive responses (n = 5, mean
± SEM) vs blood concentration (n = 4, mean ±
SEM) reveals increasing efficacy with increasing blood concentration.
(D) The efficacy of inhibitor 7 is dependent on blood
concentration. When compared, the nociceptive responses (n = 5, mean ± SEM) and the pharmacokinetic profile after oral
dosing (n = 4, mean ± SEM) of 7 follow the same trend while revealing a slight delay to in the behavioral
assay. The graphics and statistics were prepared by KaleidaGraph version
4.1 (Synergy Software, Reading, PA).Our data indicate that inhibitor 7 shows better
in
vitro efficacy in terms of both IC50 and koff than the previously published inhibitor TPAU (Figure 1C) as well as having good pharmacokinetic
profiles in both mouse and rat (Table 6 and 9). Therefore, we compared the in vivo efficacy of
inhibitor 7 as previously reported results for TPAU in the same model. Our results indicate that using a
10-fold lower dose of inhibitor 7 is as efficacious as TPAU in the nociceptive bioassay (Figure 5).[11] It is also reassuring that
the relative in vivo potency of these two sEH inhibitors in the nociceptive
assays correlates with their relative in vitro potency on the target
enzyme. These pain assays in the rat appear valuable for broad ranking
of the analgesic activity of sEH inhibitors. However, the differences
in pharmacokinetic and target site affinities caution against extrapolating
these data to predict the fine ranking of the potency of different
sEH inhibitors in man.
Figure 5
Newly optimized inhibitor 7 shows better
in vivo efficacy
in disease model. AUC describes an area under the curve of the withdrawal
threshold post-oral dosing of the diabetic rat with sEH inhibitors
vs time. A comparison of the newly synthesized sEH inhibitor 7 (IC50 rat sEH = <1.25 nM, t1/2 = 21.6 min) to previously published analogue TPAU (IC50 rat sEH = 79 nM, t1/2 = 11.9 min) showed a significantly higher response of inhibitor 7 at 0.3 mg/kg compared to the same 0.3 mg/kg dose of TPAU in a model of diabetic neuropathy (t test, t = 2.31 with 9 degrees of freedom, p = 0.046, * shows significant difference to 0.3 mpk of
inhibitor 7). However, there was no significant difference
between the 10-fold lower dose of inhibitor 7 at 0.3
mg/kg and TPAU at 3 mg/kg (Mann–Whitney rank sum
test, U = 14.00, n1 =
5 n2 = 6, ns p = 0.931).
The graphic and statistics were prepared by SigmaPlot (SysTat Software,
San Jose, CA).
Newly optimized inhibitor 7 shows better
in vivo efficacy
in disease model. AUC describes an area under the curve of the withdrawal
threshold post-oral dosing of the diabeticrat with sEH inhibitors
vs time. A comparison of the newly synthesized sEH inhibitor 7 (IC50 ratsEH = <1.25 nM, t1/2 = 21.6 min) to previously published analogue TPAU (IC50 ratsEH = 79 nM, t1/2 = 11.9 min) showed a significantly higher response of inhibitor 7 at 0.3 mg/kg compared to the same 0.3 mg/kg dose of TPAU in a model of diabetic neuropathy (t test, t = 2.31 with 9 degrees of freedom, p = 0.046, * shows significant difference to 0.3 mpk of
inhibitor 7). However, there was no significant difference
between the 10-fold lower dose of inhibitor 7 at 0.3
mg/kg and TPAU at 3 mg/kg (Mann–Whitney rank sum
test, U = 14.00, n1 =
5 n2 = 6, ns p = 0.931).
The graphic and statistics were prepared by SigmaPlot (SysTat Software,
San Jose, CA).Overall, our data show
that the new inhibitor 7 exhibits
a strong correlation between in vivo efficacy with dose and drug level
in blood, is more efficacious than the previously reported TPAU, and has promise as a therapy for use in treating chronic pain conditions.
Conclusion
Here, a new series of 1,3-disubstituentd ureasEH inhibitors was
synthesized with the design based on the recently obtained holo-structure of humansEH with inhibitor 18 (TPPU). The SAR of this new series indicates that the
right side binding pocket of the sEH enzyme has limited space for
optimization and is fluorophilic (Figures 1A, 2A,B), with an additional binding site
identified in the left side binding pocket (Figures 1A, 2D). The newly synthesized inhibitors
are at least 10 times more potent against humansEH than a previously
published inhibitor (TPAU, Figure 1C, Ki(humansEH) = 9.4 ± 0.3
nM) with a minimum 2-fold longer residence time on the humansEH.
These new inhibitors have been demonstrated selective for sEH with
no significant specific binding toward several other human proteins
(plasma protein, CYP450 enzymes, and hERG channel protein). The poor
physical properties limit the oral bioavailability of inhibitors with
a 1,3-disubstituted urea.[18,30,31,34] This new series of inhibitors
has improved physical properties translating into good pharmacokinetic
profiles in the rat and mouse (Supporting Information,
Figure S5–S11) increases their druglikeness compared
to previous sEH inhibitors.Inhibiting sEH has been demonstrated
to be more efficacious than
gabapentin and celecoxib in alleviating modeled diabetic neuropathicpain in rat with no apparent side effect such as impaired mobility,
cognition, or motor skill.[10,11] However, the sEH inhibitors
reported in those studies had either poor physical properties (poor
water solubility and high melting point) or poor in vivo stability.
In this report, the new more potent sEH inhibitors are close to 10
times more efficacious than TPAU in a diabetic neuropathicpain model. In addition, the selected inhibitor 7 demonstrated
strong correlation between drug level in blood and dose with in vivo
efficacy. This new series of inhibitors has demonstrated enhanced
potency with slow koff, improved pharmacokinetic
profiles (moderate to long elimination T1/2 and high AUC), and more importantly, improved efficacy against diabeticneuropathic pain in a rat model. With increased potency and bioavailability,
there are decreases in the required effective dose and greatly simplified
formulation. The new sEH inhibitors are good candidates for chronic
treatment of diabetic neuropathic pain.
Experimental
Section
General
All reagents and solvent were purchased from
commercial suppliers and were used directly without further purifications.
All syntheses were carried out in a dry nitrogen atmosphere unless
otherwise specified. Reactions were monitored by thin-layer chromatography
(TLC) on Merck F254 silica gel 60 aluminum sheets, and
spots were either visible under light or UV light (254 mm) or stained
with an oxidizing solution (KMnO4 stain). The same TLC
system was used to test purity, and all final products showed a single
spot on TLC. Column chromatography was performed with silica gel.1H NMR spectra were recorded on a Varian QE-300 spectrometer
with deuterated chloroform (CDCl3; δ = 7.24 ppm)
or deuterated dimethyl sulfoxide (DMSO-d6) containing TMS an internal standard. 13C NMR spectra
were recorded on a Varian QE-300 spectrometer at 75 MHz.The
purity of the inhibitors reported in this manuscript was determined
either by (1) HPLC-UV using Agilent 1200 series HPLC series equipped
with Phenomenex Luna2 C18 reverse phase column (C18, 4.6 mm ×
150 mm, 5 μm) coupled with Agilent G1314 UV–vis detector
(detection at 230 nm) with isocratic flow at methanol:water (2:1 by
volume) for 90 min, or by (2) H NMR. The lowest obtained purity was
reported. The inhibitor was dissolved in EtOH at 100 μM and
10 μL was injected on HPLC. Purity was based on the percent
of total peak area at 230 nm using HPLC-UV. This purity estimate was
compared with that from the H NMR. The presence of anilines in the
final product was estimated from H NMR. The lowest obtained purity
was reported. The purity was also further supported as described in
the Supporting Information by LC/MS with
monitoring of total ion current, TLC in several systems, a sharp melting
point, and occasional other technique. The elemental analysis was
conducted by MIDWESTMICRO lab, LCC.The synthesis of tert-butyl 4-(3-(4-(trifluoromethyl)phenyl)ureido)piperidine-1-carboxylate,
1-(piperidin-4-yl)-3-(4-(trifluoromethyl)phenyl)urea, 1-(piperidin-4-yl)-3-(4-(trifluoromethoxy)phenyl)urea
and tert-butyl 4-(3-(4-(trifluoromethoxy)phenyl)ureido)piperidine-1-carboxylate
have been reported elsewhere.[23,30,53] The synthesis of inhibitors 1,[38]2,[23]3,[23]6,[38]APAU,[32]TPAU,[23]18 (TPPU),[23]20 (TPCU),[23] and t-TUCB[34] were reported elsewhere. Experimental
of each individual inhibitor is described in Supporting
Information in detail.The experimental procedures of
enzyme preparation, IC50 determination for sEH inhibitors,
and Ki determination for sEH inhibitors
followed the published procedures
and are described in detail in Supporting Information.[38,54−58]
Synthetic Method 1
Step 1
Corresponding isocyanate (1 equiv)
and 4-amino-1-Boc-piperidine (1.1 equiv) were dissolved in CH2Cl2 (50 mM, corresponding to isocyanate) and stirred
at rt for 12 h. The reaction was quenched by addition of water. The
organic layer was isolated, and the aqueous layer was extracted with
EtOAc (EtOAc:aqueous layer/1:1) four times. The combined organic layer
was dried over anhydrous magnesium sulfate and was concentrated under
vacuo and was further purified by flash chromatography, yielding corresponding
Boc-protected urea.
Step 2
The BOC protected
urea from the
step 1 was dissolved in HCl solution (2M, MeOH) to make reaction mixture
(186 mM, BOC protected urea). The resulting solution was refluxed
for 2 h. The solvent was removed in vacuo, and the crude reaction
product was adjusted to pH 12 with NaOH solution (6N). The precipitates
were filtered and dried under high vacuum. The final product unprotected
urea was served as a scaffold for the next step of synthesis.
Step 3
Unless specified, the unprotected
urea (1 equiv) from step 2, EDCI (1.5 equiv), DMAP (1.5 equiv), and
corresponding carboxylic acid (1.5 equiv) were dissolved in CH2Cl2 (8.3 mM, unprotected urea) and stirred overnight
(12 h) at rt. The reaction was quenched by addition of HCl solution
(1M). The organic layer was collected, and the aqueous layer was extracted
with EtOAc (EtOAc:aqueous layer/1:1) four times. The combined organic
layer was dried over anhydrous magnesium sulfate and was concentrated
in vacuo and further purified by flash chromatography.
Synthetic
Method 2
The corresponding isocyanate
(1 equiv) was added to a suspension of targeted piperidine (1.1 equiv)
in CH2Cl2 (20 mM, corresponding isocyanate).
The reaction was stirred overnight (12 h) at rt. The reaction was
quenched with the addition of HCl solution (2M). The organic layer
was collected, and the aqueous layer was further extracted with EtOAc
(EtOAc:aqueous layer/1:1) three times. The combined organic layer
was washed with satd NaCl solution. The organic layer was dried over
anhydrous magnesium sulfate and was concentrated in vacuo. The product
was purified by flash chromatography.
Synthetic Method 3
Corresponding amine
(1 equiv) and triethylamine (1.2 equiv) was dissolved in CH2Cl2 (54 mM corresponding to amine) and stirred at −78
°C. Triphosgene (0.37 equiv) dissolved in CH2Cl2 (20 mM, corresponding triphosgene) was added dropwise at
−78 °C. The reaction was then warm to rt and was stirred
for 30 min. The reaction was cooled to 0 °C. Corresponding piperidine
(1.1 equiv) dissolved in CH2Cl2 (54 mM, corresponding
piperidine) was added slowly, and the reaction was further stirred
at rt for 12 h. The reaction was quenched with the addition of HCl
solution (2M). The organic layer was collected, and the aqueous layer
was further extracted with EtOAc (EtOAc:aqueous layer/1:1) three times.
The combined organic layer was washed with satd NaCl solution. The
organic layer was dried over anhydrous magnesium sulfate and was concentrated
in vacuo. The product was purified by flash chromatography.
Synthetic
Method 4
The first two steps
are the same as synthetic method 1, steps 1 and 2, unless specified.The unprotected urea (1 equiv)
and triethylamine (1.2 equiv) was dissolved in CH2Cl2 (8.3 mM, corresponding unprotected urea), and corresponding
sulfonyl chloride was added dropwise at 0 °C and the reaction
was stirred overnight (12 h) at rt. The reaction was quenched by addition
of HCl solution (1M). The organic layer was collected, and the aqueous
layer was extracted with EtOAc (EtOAc:aqueous layer/1:1) four times.
The combined organic layer was dried over anhydrous magnesium sulfate
and was concentrated in vacuo and further purified by flash chromatography.
Protein Crystallization
Crystals of the enzyme were
obtained using the hanging drop vapor-diffusion method by mixing equal
volumes of protein (8–12 mg/mL concentration in 100 mM sodium
phosphate pH 7.4, 3 mM DTT) and the reservoir solution (30% PEG 3350,
0–10% sucrose) at 4 °C. The crystals grew in approximately
1 week and belonged to the hexagonal space group P6522.Complexes of sEH with inhibitors 8/UC1770 or 4 have been obtained by soaking sEH crystals
grown as described above in modified mother liquor (35% PEG 3350,
50 mM sodium phosphate pH 7.4) supplemented with 1 mM solution of
inhibitor for 1–7 days
Data Collection
Prior to the data
collection, a suitable
crystal was dipped for 30 s in a modified mother liquor solution (35%
PEG 3350) with the addition of 10% glycerol as a cryoprotectant. Diffraction
data were collected at 100 K at the XP station at the Center for Advance
Microstructures and Devices at Louisiana State University with a MAR
charge-coupled device camera (structure TPPU/UC1770)
or the NE-CAT beamline 24-ID-C at the Advanced Photon Source equipped
with the Pilatus 6 M detector (structure 4). The images
were processed and scaled using the HKL2000 (structure TPPU/UC1770)[59] or Xia2 program suit[60] (structure 4). Data collection
and data processing statistics are given in Table 1.
Crystal Structure Determination
The molecular replacement
procedure was applied to locate a solution using the program MOLREP.[61] A monomer of humansEH (PDB accession code 1S8O) was used as a search
model. The positioned MR model was refined using the maximum likelihood
refinement in REFMAC[61] with the TLS parameters
generated by the TLSMD server.[62] Program
Coot was used for model building throughout the refinement.[63]
HERG Patch Clamp Measurement
Cell Culture
for Channel Analysis
HEK-293 cells stably
expressing hKv11.1 (hERG) under G418 selection were a generous gift
from Craig January (University of Wisconsin, Madison). Cells were
cultured in DMEM containing 10% fetal bovine serum, 2 mM glutamine,
1 mM Na+ pyruvate, 100 U/mL penicillin and 100 μg/mL
streptomycin, and 500 mg/mL G418. Before electrophysiological experiments,
cells were grown to 60% confluency and 10 mM astemizole was added
to the culture for 24 h to increase Kv11.1 surface expression.
Electrophysiology
All experiments were conducted with
an EPC-10 amplifier (HEKA, Lambrecht/Pfalz, Germany) in the whole-cell
configuration of the patch-clamp technique with a holding potential
of −80 mV. Pipette resistances averaged 2.0 MΩ. Compound
solutions in Na+-Ringer were freshly prepared during the experiments
from 10 mM stock solutions in DMSO. The final DMSO concentration never
exceeded 1%. The external solution contained in mM: 160 NaCl, 4.5
KCl, 2 CaCl2, 1 MgCl2, and 10 HEPES, pH 7.4,
osmolality 300 mmol/kg. The internal pipet solution contained in mM:
120 KCl, 10 HEPES, 4 Na2ATP, 10 EGTA, 5.374 CaCl2, 1,75 MgCl2, pH 7.2, osmolality 295 mmol/kg (free Ca2+ concentration 150 nM calculated with MaxChelator assuming
a temperature of 25 °C, a pH of 7.2, and an ionic strength of
160 mM). HERG (Kv11.1) currents were elicited with a two-step pulse
(applied every 10 s) from −80 mV first to 20 mV for 1 s and
then to −50 mV for 2 s, and the percent reduction of both peak
and tail current by the drug were determined.
Pharmacokinetic
Study of Inhibitors Using Oral Dosing in Mice
All the animal
experiments were performed according to the protocols
approved by the Animal Use and Care Committee of University of California—Davis.
Male Swiss Webster mice (8 week old, 24–30 g) purchased from
Charles River Laboratories were used for PK studies. Inhibitors for
oral administration were dissolved in oleic acid-rich triglyceride
containing 20% PEG400 (v/v) to give a clear solution. Blood (10 μL)
was collected from the tail vein using a pipet tip rinsed with 7.5%
EDTA(K3) at 0, 0.5, 1, 1.5, 2, 4, 6, 8, 24, and 48 h after
administration of the inhibitor in a cassette of 3–5 compounds
(Supporting Information, Table S1) (0.3
mg/kg per compounds, 100–110 uL). Each group contained 3–4
animals. Each blood sample was immediately transferred to a tube containing
50 μL of water containing 0.1% EDTA. After being mixed strongly
on a Vortex for 1 min, all samples were stored at −80 °C
until analysis. The blood samples were prepared for the measurement
of sEH inhibitors according to the previously reported method by Liu
et al.[18] The details LC/MS–MS methods
are described in detail in the Supporting Information.
Pharmacokinetic Study of Inhibitors Using Oral Dosing in Rat
All the animal experiments were performed according to the protocols
approved by the Animal Use and Care Committee of University of California—Davis.
Male Sprague–Dawley rats (n = 4, 8 week old,
250–300 g) were used for pharmacokinetic study for sEH inhibitors.
A cassette of four inhibitors (inhibitors 4, 7, 19, and 21, 0.3 mg/kg per inhibitors,
0.9–1.2 mL) was given by oral administration. Inhibitor was
dissolved in oleic oil containing 5% polyethylene glycol 400 to form
a clear solution. Blood (10 μL) was collected from the tail
vein by using a pipet tip rinsed with 7.5% EDTA(K3) at
0, 0.5, 1, 1.5, 2, 4, 6, 8, and 24 h after oral dosing with the inhibitor.
Each blood sample was immediately transferred to a tube containing
50 μL of water and mixed by vortex for 1 min, and all samples
were stored at −80 °C until analysis. The blood samples
were prepared for the measurement of sEH inhibitors according to the
previously reported method by Liu et al.[18] The details LC/MS–MS methods are described in detail in the Supporting Information.
Diabetic Neuropathic Pain
Model
A diabetic neuropathicpain modeled was generated using streptozocin which targets and kills
the pancreatic beta islet cells, rendering the animals with type I
diabetes. The rats were acclimated for 1 h and tested for baseline
thresholds before inducing diabetes. The baseline mechanical withdrawal
thresholds were established using the von Frey mechanical nociceptive
test with an electronic anesthesiometer (IITC, Woodland Hills, CA).
Subsequently, streptozocin (55 mg/kg) in saline was injected via tail
vein per previously reported methods.[64] After 5 days, the allodynia of diabeticrats was confirmed with
the von Frey nociceptive assay. Rats were placed in clear acrylic
chambers on a steel mesh floor. The hind paw of the rat was probed
through the mesh with a rigid tip probe connected to the electronic
readout pressure meter set to the maximum hold setting. The withdrawal
thresholds per rat were measured 3–5 times at 1 min intervals
for each time point.The baseline diabetic allodynia was quantified
again at the beginning of all test days. The rats were then administered
vehicle or sEH inhibitor via oral gavage and tested at 30 min, 1,
2, 3, 4, 5, 6, and 8 h for mechanical withdrawal thresholds. The reported
scores are the grams of force required to elicit a hind paw withdrawal
averaged with standard error of the mean (SEM) per a group of rats
(n = 5). The baseline diabetic neuropathic scores
were normalized to 100% to reflect the response to treatments which
are reported as % of post diabetic neuropathic baseline.
Authors: Paul D Jones; Nicola M Wolf; Christophe Morisseau; Paul Whetstone; Bertold Hock; Bruce D Hammock Journal: Anal Biochem Date: 2005-08-01 Impact factor: 3.365
Authors: Christophe Morisseau; Bora Inceoglu; Kara Schmelzer; Hsing-Ju Tsai; Steven L Jinks; Christine M Hegedus; Bruce D Hammock Journal: J Lipid Res Date: 2010-07-27 Impact factor: 5.922
Authors: Arzu Ulu; Todd R Harris; Christophe Morisseau; Christina Miyabe; Hiromi Inoue; Gertrud Schuster; Hua Dong; Ana-Maria Iosif; Jun-Yan Liu; Robert H Weiss; Nipavan Chiamvimonvat; John D Imig; Bruce D Hammock Journal: J Cardiovasc Pharmacol Date: 2013-09 Impact factor: 3.105
Authors: Hong C Shen; Fa-Xiang Ding; Siyi Wang; Qiaolin Deng; Xiaoping Zhang; Yuli Chen; Gaochao Zhou; Suoyu Xu; Hsuan-Shen Chen; Xinchun Tong; Vincent Tong; Kaushik Mitra; Sanjeev Kumar; Christine Tsai; Andra S Stevenson; Lee-Yuh Pai; Magdalena Alonso-Galicia; Xiaoli Chen; Stephen M Soisson; Sophie Roy; Bei Zhang; James R Tata; Joel P Berger; Steven L Colletti Journal: J Med Chem Date: 2009-08-27 Impact factor: 7.446
Authors: Kerstin Hiesinger; Jan S Kramer; Janosch Achenbach; Daniel Moser; Julia Weber; Sandra K Wittmann; Christophe Morisseau; Carlo Angioni; Gerd Geisslinger; Astrid S Kahnt; Astrid Kaiser; Anna Proschak; Dieter Steinhilber; Denys Pogoryelov; Karen Wagner; Bruce D Hammock; Ewgenij Proschak Journal: ACS Med Chem Lett Date: 2019-05-29 Impact factor: 4.345
Authors: W Swardfager; M Hennebelle; D Yu; B D Hammock; A J Levitt; K Hashimoto; A Y Taha Journal: Neurosci Biobehav Rev Date: 2018-02-02 Impact factor: 8.989
Authors: Kin Sing Stephen Lee; Niel M Henriksen; Connie J Ng; Jun Yang; Weitao Jia; Christophe Morisseau; Armann Andaya; Michael K Gilson; Bruce D Hammock Journal: Arch Biochem Biophys Date: 2016-10-29 Impact factor: 4.013
Authors: Yue Zhang; Gina Hong; Kin Sing Stephen Lee; Bruce D Hammock; Debebe Gebremedhin; David R Harder; Raymond C Koehler; Adam Sapirstein Journal: J Neurochem Date: 2017-01-23 Impact factor: 5.372
Authors: Anamitra Ghosh; Michele M Comerota; Debin Wan; Fading Chen; Nicholas E Propson; Sung Hee Hwang; Bruce D Hammock; Hui Zheng Journal: Sci Transl Med Date: 2020-12-09 Impact factor: 17.956
Authors: Bora Inceoglu; Ahmed Bettaieb; Carlos A Trindade da Silva; Kin Sing Stephen Lee; Fawaz G Haj; Bruce D Hammock Journal: Proc Natl Acad Sci U S A Date: 2015-07-06 Impact factor: 11.205
Authors: Seiya Kitamura; Kelli L Hvorecny; Jun Niu; Bruce D Hammock; Dean R Madden; Christophe Morisseau Journal: J Med Chem Date: 2016-05-04 Impact factor: 7.446
Authors: Kin Sing Stephen Lee; Jen C Ng; Jun Yang; Sung-Hee Hwang; Christophe Morisseau; Karen Wagner; Bruce D Hammock Journal: Bioorg Med Chem Date: 2020-08-31 Impact factor: 3.641