The von Hippel-Lindau tumor suppressor protein is the substrate binding subunit of the VHL E3 ubiquitin ligase, which targets hydroxylated α subunit of hypoxia inducible factors (HIFs) for ubiquitination and subsequent proteasomal degradation. VHL is a potential target for treating anemia and ischemic diseases, motivating the development of inhibitors of the VHL:HIF-α protein-protein interaction. Additionally, bifunctional proteolysis targeting chimeras (PROTACs) containing a VHL ligand can hijack the E3 ligase activity to induce degradation of target proteins. We report the structure-guided design and group-based optimization of a series of VHL inhibitors with low nanomolar potencies and improved cellular permeability. Structure-activity relationships led to the discovery of potent inhibitors 10 and chemical probe VH298, with dissociation constants <100 nM, which induced marked HIF-1α intracellular stabilization. Our study provides new chemical tools to probe the VHL-HIF pathways and new VHL ligands for next-generation PROTACs.
The von Hippel-Lindau tumor suppressor protein is the substrate binding subunit of the VHL E3 ubiquitin ligase, which targets hydroxylated α subunit of hypoxia inducible factors (HIFs) for ubiquitination and subsequent proteasomal degradation. VHL is a potential target for treating anemia and ischemic diseases, motivating the development of inhibitors of the VHL:HIF-α protein-protein interaction. Additionally, bifunctional proteolysis targeting chimeras (PROTACs) containing a VHL ligand can hijack the E3 ligase activity to induce degradation of target proteins. We report the structure-guided design and group-based optimization of a series of VHL inhibitors with low nanomolar potencies and improved cellular permeability. Structure-activity relationships led to the discovery of potent inhibitors 10 and chemical probe VH298, with dissociation constants <100 nM, which induced marked HIF-1α intracellular stabilization. Our study provides new chemical tools to probe the VHL-HIF pathways and new VHL ligands for next-generation PROTACs.
Proteins are constantly
synthesized and targeted for degradation
during cellular homeostasis. The major pathway to protein degradation
is the ubiquitin–proteasome system (UPS), first reported in
the early 1980s.[1] In a highly regulated
enzymatic signaling cascade, ubiquitin is covalently attached to a
substrate protein as a monomer or as a polyubiquitin chain, promoting
proteasome-catalyzed target degradation.[2,3] The complexity
and biological importance of the UPS, together with implication of
this pathway in a wide range of diseases, highlight its importance
as a therapeutic target.[4,5] The market release of
the pioneering proteasome inhibitors bortezomib and carfilzomib, for
the treatment of hematopoietic and lymphoid malignancies, validated
the UPS as a focus for therapeutic intervention and opened the door
to a new series of proteasome inhibitors.[5−7] Despite their
success, proteasome inhibitors lack substrate specificity, which is
a significant limitation to their use as chemical probes of specific
biological pathways and leads to unwanted side effects in their use
as drugs.[7,8]Specific binding of a target substrate
to an E3 ubiquitin ligase
is required for substrate ubiquitination and is therefore a key step
in the ubiquitination process.[3,9] E3 ubiquitin ligases
determine specificity of substrate ubiquitination and thus could be
regarded as more attractive targets to center drug discovery efforts
over the proteasome. However, small-molecule modulators of E3 ligases
need to target protein–protein interactions (PPIs), either
directly or allosterically.[10] The typically
flat and featureless nature of many PPI interfaces can be a limiting
factor in the development of potent and selective inhibitors.[11,12] To cover considerable protein surface area that often lacks well-defined
pockets, most PPI inhibitors tend to be larger in size than classical
enzyme inhibitors or receptor antagonists, forcing medicinal chemists
to break conventional drug-like quality metrics.[13−15] These observations
consequently led to E3 ligases being perceived as untractable target
space.[16] On the other hand, the identification
of small-sized hotspots for certain PPIs, particularly those characterized
by the recognition of primary-type linear epitopes, has led to the
successful observation of small-molecule-sized patches on protein
surfaces that could be adequately targeted by drug-like molecules.[17−19] Today, several small molecules that modulate PPIs have been developed
as chemical probes, and many have entered clinical trials.[19,20] These realizations have reinvigorated drug discovery efforts for
a wide range of PPI-based targets, including E3 ubiquitin ligases.[21]The von Hippel–Lindau protein (VHL)
is a member of the Cullin-RING
ligase family of E3 ubiquitin ligases.[21] The major substrate of VHL is the hypoxia inducible factor 1α
(HIF-1α), a transcription factor that drives the transcriptional
program of many human genes,[22] mainly involved
in hypoxia adaptation.[23,24] When tissue oxygen levels are
normal, HIF-1α expression levels are tightly controlled. Iron-
and oxygen-dependent hydroxylation occurs at two specific proline
residues within the HIF-1α oxygen-dependent degradation domain
(ODD) by prolyl hydroxylase domain (PHD) enzymes. This modification
leads to HIF-1α specific recognition and ubiquitination by VHL
and subsequent degradation via the UPS (Supporting Information, Figure 1).[25−28]The importance of this pathway has been uncovered
in a wide range
of diseases, including conditions characterized by anemia, ischemia,
inflammation, chronic neurodegeneration, and more recently mitochondrial
dysfunction.[29−32] Small-molecule inhibition of this pathway could activate HIF-1α
expression, upregulating genes involved in the hypoxic response, consequently
providing a potential therapeutic strategy. Indeed, small-molecule
PHD inhibitors have been developed that have shown potential in a
number of pathologies.[33−35] Among these, the candidate FG-4592 (N-[(4-hydroxy-1-methyl-7-phenoxy-3-isoquinolinyl)carbonyl]glycine,
also known as roxadustat or ASP1517) has progressed to a phase III
clinical trial to evaluate efficacy and safety in hemodialysis chronic
kidney patients with anemia.[36,37] Despite their clinical
potential, no PHD inhibitor shows selectivity between the three different
PHD paralogues (PHD1, PHD2, and PHD3).[38] In addition, several non-HIF substrates of PHD enzymes have been
identified.[39,40] Lack of selectivity could promote
unwanted side effects due to the different substrate specificities
as well as cellular and tissue localization of PHD enzymes.[38]In contrast to PHD inhibitors, inhibitors
of the VHL:HIF-1α
PPI would allow blocking the pathway downstream of HIF hydroxylation
by PHD enzymes.[41] This approach could avoid
HIF-independent off-target effects and deliver new lead compounds
as hypoxia mimetics.[42] To validate the
chemical intervention on VHL and establish the biological consequences
of modulating its activity, small molecules must effectively penetrate
cells and engage with VHL with enough potency to productively block
HIF-1α binding. In a different approach, VHL ligands can be
conjugated to target ligands of interest to design VHL-recruiting
bivalent chimeric ligands (also known as PROTACs).[43] PROTACs can efficiently hijack the VHL ligase to induce
intracellular degradation of the desired target proteins. Previous
work from our laboratory and others have demonstrated potent and selective
activities of VHL-recruiting PROTACs against a wide range of target
proteins, including Brd4,[44−48] RIPK2,[49] and c-Abl kinase among others.[50]Co-crystal structures of VHL bound with
a hydroxyproline peptide
from human HIF-1α carboxy-terminal oxygen-dependent degradation
(CODD) motif identified a groove on the pVHL surface, which accommodated
the substrate peptide.[28,51] Recognizing that this primary
type PPI could be suitably targeted by small-molecule disruptors,
in pioneering work, our laboratory in collaboration with the Crews
laboratory developed a first-generation of VHL ligands with single-digit
micromolar binding affinities based around the central hydroxyproline
core fragment.[17,52−54] However, these
initial compounds proved inactive as PPI inhibitors in cells. Subsequent
structure-guided optimization led to more potent inhibitors with improved
nanomolar binding affinities.[55] The best
inhibitor of this series, compound 1 (VH032, Figure A) was able to disrupt
the VHL:HIF-1α PPI in cells, illustrated by its capacity to
stabilize hydroxylated HIF-1α.[41] Nevertheless,
this inhibitor showed low passive cell permeability and exhibited
a large mismatch between in vitro potency in solution and in cells.[41]
Figure 1
Initial optimization strategy. (A) X-ray crystal structure
of VHL,
elongin B, and elongin C (VBC) protein complex with inhibitor 1 (purple carbons, PDB 4W9H).[55] VHL is
shown as a pale-green surface and the VHL residues forming the binding
pocket as orange stick representations. (B) Designed group modifications
on inhibitor 1 to improve activity.
Initial optimization strategy. (A) X-ray crystal structure
of VHL,
elongin B, and elongin C (VBC) protein complex with inhibitor 1 (purple carbons, PDB 4W9H).[55] VHL is
shown as a pale-green surface and the VHL residues forming the binding
pocket as orange stick representations. (B) Designed group modifications
on inhibitor 1 to improve activity.Motivated by a desire to drive further inhibitor optimization,
analysis of the X-ray crystal structure of the complex composed of
VHL, elongin B, and elongin C (VBC) with 1 suggested
an attractive pocket at the left-hand side (LHS) of the VHL:HIF interface
that was only partially filled by the terminal acetamide group of 1 and could therefore provide space to further grow and optimize
the molecule (Figure A). Herein, we describe a systematic group-based optimization of
VHL inhibitors starting from 1, improving binding affinity,
cell membrane passive permeability, and ultimately cellular activity,
a campaign that led to the discovery of inhibitor 15 (VH298),
a potent and selective probe of the hypoxia HIF-1α pathway.[41]
Results and Discussion
VHL Inhibitor SAR
We aimed to optimize the physicochemical
properties and binding affinity of VHL inhibitors in order to minimize
the observed mismatches between biophysical and cellular potency.[41] The information gathered from the co-crystal
structure of inhibitor 1 bound to VBC protein complex
offered a starting point for the design of a new series of inhibitors
by adding small modifications on the LHS of the inhibitor (Figure A).In the
first new series of inhibitors, we decided to keep the carbonyl on
the LHS in order to maintain the hydrogen bond between this group
and the structural water in this pocket. Our design strategy was to
replace the three hydrogens of the acetamide methyl group, one at
a time, using alkyl groups to better fill the LHS pocket or adding
electron-withdrawing groups at the α position of the LHS amide
in order to lock its conformation.[56,57] It was also
decided to evaluate the effect of removing the LHS acetamide hydrogen-bond
donor group (NH) via cyclization, which we hypothesized could improve
cellular permeability (Figure B).[58] The synthesized compounds
were evaluated biophysically in a direct binding assay using isothermal
titration calorimetry (ITC) and for their ability to displace a high-affinity
HIF peptide using fluorescence polarization (FP). In parallel, the
compounds cellular activity was evaluated in HeLa cells by monitoring
protein levels of HIF-1α using Western blot (Table ).
Table 1
Chemical
Structure, FP Kds Back-Calculated from
IC50s,[17] ITC Measured Kds, and Group
Efficiencies (GE) [Based on FP Results and Calculated as GE = ΔpKd/ΔHA (Heavy Atoms)] of Inhibitors 1 and 2–12
HeLa cells were treated with 50
μM of the respective inhibitor and HIF-1α stabilization
levels were measured by Western blot after 2 h treatment (protein
levels normalized to HIF-1α stabilization level observed with 1).
HeLa cells were treated with 50
μM of the respective inhibitor and HIF-1α stabilization
levels were measured by Western blot after 2 h treatment (protein
levels normalized to HIF-1α stabilization level observed with 1).Substitution
of a single hydrogen on the acetamide group for a
hydroxyl (3) maintained binding affinity but led to a
substantial loss in cellular potency. On the other hand, expansion
of the alkyl chain from methyl to ethyl (4) allowed the
retention of inhibitor cellular potency despite a 2-fold loss in inhibitor
binding affinity. These two observations suggested that tuning hydrogen
bond donors and lipophilicity could play as important a role as binding
affinity in determining cellular activity. Replacing all three hydrogens
with fluorine (7) or chlorine groups (8)
did not impact binding affinity and cellular potency. In contrast,
introduction of a bulky tert-butyl group (5) resulted in a loss in binding affinity and in cellular potency,
a similar effect observed when replacing all hydrogens for two methyl
groups together with one fluorine group (9). However,
substitution of two hydrogens with a restrained cyclopropyl group
(6) promoted a small improvement in both binding affinity
and cellular activity. Further replacement of the tertiary hydrogen
of 6 with a fluorine (10) led to a marked
increase in both binding affinity and cellular potency compared to 1 and 6 (Table ). Figure illustrates representative ITC and FP binding curve for 10, the most potent inhibitor of the series, and the parent
inactive epimer cis-10, which due to
a change in the proline OH group orientation no longer binds to VHL.
Figure 2
Biophysical
characterization of inhibitor 10 binding
to VBC. (A) Competitive FP binding assay curve, monitoring the displacement
of a 20-mer FAM-labeled HIF-1α peptide bound from VBC by inhibitor 10 (Kd = 3 nM). (B) ITC titrations
of 300 μM inhibitor 10 (blue) or 300 μM of
its cis inactive epimer (green) into 30 μM
VBC protein.
Biophysical
characterization of inhibitor 10 binding
to VBC. (A) Competitive FP binding assay curve, monitoring the displacement
of a 20-mer FAM-labeled HIF-1α peptide bound from VBC by inhibitor 10 (Kd = 3 nM). (B) ITC titrations
of 300 μM inhibitor 10 (blue) or 300 μM of
its cis inactive epimer (green) into 30 μM
VBC protein.Attempts to modify other
regions of the chemical structure of 1 did not yield
improvements. For example, modification of
the tert-leucine side chain for a less bulky alanine
side chain (11), or removal of the acetamide hydrogen
bond donor group through cyclization into a pyrrolidinone (12) led to significant loss of inhibitor binding affinity and cellular
potency. Last, acetamide methyl group substitution for a Boc protecting
group (2) was not tolerated, resulting in the greatest
loss in binding and cellular potency.To better understand the
SAR of this new inhibitor series, the
X-ray crystal structures of inhibitors 3, 6, 10, and 11 bound to VBC were obtained.
All solved structures revealed a conserved inhibitor binding mode
at the VHL:HIF interface (Figure ), consistent with the binding of 1.[55] The crystal structure of 3 shows
two different conformations upon binding to VBC, each with 50% of
occurrence. In the first inhibitor conformation, the added OH group
points upward, forming an extra hydrogen bond with a structural water
that promotes the formation of an interaction network between Tyr112
side chain and the amide carbonyl on the LHS of inhibitor’s
hydroxyproline (Hyp) (Figure A). In the second conformation, the hydroxyl group formed
a bidentate interaction with the structural water present on the far
left side of the binding pocket, together with LHS amide carbonyl
(Figure B). These
favorable interactions contribute to a small improvement in inhibitor
binding affinity (Table ).
Figure 3
Co-crystal structures of first-series inhibitors. Crystal structures
of VBC in complex with 3 (A,B), 6 (C), 10 (D), and 11 (E). VHL inhibitors are shown
as sticks with purple carbons. VHL is shown as a pale-green surface
and the VHL residues forming the binding pocket as orange stick representations.
Inhibitor-binding waters are shown as red spheres. Hydrogen bond interactions
between inhibitors, bound waters, and VHL pocket residues are shown
as dashed red lines.
Co-crystal structures of first-series inhibitors. Crystal structures
of VBC in complex with 3 (A,B), 6 (C), 10 (D), and 11 (E). VHL inhibitors are shown
as sticks with purple carbons. VHL is shown as a pale-green surface
and the VHL residues forming the binding pocket as orange stick representations.
Inhibitor-binding waters are shown as red spheres. Hydrogen bond interactions
between inhibitors, bound waters, and VHL pocket residues are shown
as dashed red lines.Crystal structures of 6 (Figure C) and 10 (Figure D) showed the cyclopropyl group
occupying the far LHS of the VHL–HIF interface, maintaining
the hydrogen bond interactions with the structural water. The addition
of cyclopropyl induces a conformational modification in the residues
on this side of the pocket to accommodate this moiety. The most notable
conformational change is observed on the side chain of Arg69, which
adopts a bent conformation to accommodate the new group on the LHS.
From the crystal structure of inhibitor 10, the fluorine
atom at the α position is found anti to the
amide carbonyl. This matches a known minimum-energy conformation for
α-fluoroamides.[56,57] This effect allows the inhibitor
to be preorganized in its bound conformation prior to binding, thus
minimizing the entropic penalty to binding and as a result increasing
binding affinity.Substitution of the tert-leucine
group for a less
bulky alanine in 11 was meant to more closely mimic the
sequence of the HIF-1α peptide (Supporting Information, Figure 2). This modification however led to a
less rigid conformation on the inhibitor LHS. In the co-crystal structure
of 11 (Figure E), the ligand bends upward, increasing the distance between
the acetamide carbonyl and the structural water from 2.8 ± 0.1
to 3.8 ± 0.5 Å (Supporting Information, Figure 3), consequently weakening this hydrogen bond and decreasing
the binding affinity.The first inhibitor series led us to the
discovery of 10, a VHL inhibitor with double-digit nanomolar
binding affinity, good
cell membrane permeability, and high activity in cells (Table and Figure A), which are some of the key parameters
needed for a good chemical probe.[59] Despite
meeting these benchmarks, further investigation revealed a degree
of cytotoxicity for 10, which notably decreased cell
proliferation relative to DMSO vehicle (Figure B,C). The inactive epimer cis-10 still presented a degree of cytotoxicity (Figure B,C). Together, these
findings suggested an off-target toxicity for 10 (Figure ).
Figure 4
Cellular activity and
toxicity of 10. (A) Activity
of inhibitor 10 in HeLa cells expressing a hypoxia responsive
element (HRE)-luciferase reporter and treated with indicated concentrations
for 24 h. Luminescence reports luciferase expression as proxy of HIF-α
activity.[41] Colony formation (B) and cell-proliferation
assays (C) showing toxicity in HeLa cells after treatment with 150
μM of inhibitor 10 and its cis epimer.
Cellular activity and
toxicity of 10. (A) Activity
of inhibitor 10 in HeLa cells expressing a hypoxia responsive
element (HRE)-luciferase reporter and treated with indicated concentrations
for 24 h. Luminescence reports luciferase expression as proxy of HIF-α
activity.[41] Colony formation (B) and cell-proliferation
assays (C) showing toxicity in HeLa cells after treatment with 150
μM of inhibitor 10 and its cis epimer.The observed cytotoxicity of 10 motivated us to identify
a more suitable chemical probe. On the basis of the information furnished
by the co-crystal structure, we designed a second series of VHL inhibitors
(Figure ).
Figure 5
Second series
optimization strategy. (A) X-ray crystal structure
of VBC protein complex with inhibitor 10. VHL is shown
as a pale-green surface and the VHL residues forming the binding pocket
as orange stick representations. (B) Designed group modifications
on inhibitor structure to improve inhibitor activity and optimize
inhibitors toward a chemical probe.
Second series
optimization strategy. (A) X-ray crystal structure
of VBC protein complex with inhibitor 10. VHL is shown
as a pale-green surface and the VHL residues forming the binding pocket
as orange stick representations. (B) Designed group modifications
on inhibitor structure to improve inhibitor activity and optimize
inhibitors toward a chemical probe.In this second series, we initially studied the impact of
changing
the fluorine atom for other electron withdrawing groups or promoters
of intramolecular hydrogen bonds. Second, to exploit the flexibility
of Arg69 in the LHS pocket, we expanded the cyclopropyl ring by adding
larger and more lipophilic groups. In addition, we considered moieties
that could potentially form a hydrogen bond with the Arg69 side chain
or, alternatively, together with the LHS carbonyl amide, could form
bidentate hydrogen bonds with the structural water present in the
LHS pocket. Finally, we decided to test the impact of methylating
the solvent-exposed NH group of the LHS acetamide and of replacing
the tert-leucine side chain for threonine to explore
potential hydrogen bonds at the top LHS pocket, as observed previously
(Figure B).[53]Removing the hydrogen bond donor group
through methylation (13) resulted in a major loss in
binding affinity and cellular
potency when compared with 1, presumably due to destabilization
of the required trans-amide bond conformation. Substitution
of tert-leucine on 10 for a threonine
(14) side chain group also led to a loss in binding affinity
and cellular activity. Conversion of the fluorine group into an electron-withdrawing
group such a cyano (15) revealed a 2-fold increase in
binding affinity and in cellular potency compared to 1. Changing the fluorine group for a ketone moiety (16) also improved both binding affinity and cellular activity. Modification
of the fluorine atom with an acetamide (17) retained
a comparable affinity when compared with 1 but yielded
a considerable loss in cellular potency, presumably as a result of
poorer cell membrane passive permeability. Expansion of the cyclopropyl
ring into a cyclobutyl (18) led to a small loss of inhibitor
binding affinity but a greater cellular activity, presumably due to
the increased lipophilicity and permeability. Conversion of the cyclopropyl
group into an oxetane (19) retained binding affinity
but led to a 4-fold loss in cellular activity. Further conversion
of the cyclopropyl group into a cyclobutanone (20) or
into an acetylated azetidine (21) led to a 2- and 20-fold
loss in binding affinity and a 5- and 11-fold decrease in cellular
potency, respectively.In summary, while we were not able to
improve binding affinity
beyond that of inhibitor 10, most of the designed inhibitors
retained in vitro binding potencies around 100 nM and also retained
high cellular activity, in many cases greater than starting inhibitor 1. In particular, 15 and 10 were
the only compounds that showed Kd <
100 nM and positive group efficiencies (Table ). Inhibitor 15 was shown to
be highly cell permeable (vide infra, Table ) and exhibited cellular activity comparable
with 10 and significantly greater than 1. Further investigation with 15 revealed no toxicity
in cells (Figure ),
supporting further evaluation of it as a chemical probe.
Table 2
Chemical Structure, FP Back-Calculated Kds, ITC Measured Kds, and Group Efficiencies
(GE) [Based on FP Results and Calculated
as GE = ΔpKd/ΔHA (Heavy Atoms)]
of Inhibitors 1, 10, and 13–21a
HeLa cells were treated with 50
μM inhibitors and HIF-1α stabilization levels were measured
by Western blot after 2 h treatment (protein levels normalized to
HIF-1α stabilization level observed with 1).
Table 3
Structures, Kds Measured by ITC, Calculated Ligand Efficiency
(LE) [LE =
−ΔG/NHA = −RT ln Kd/NHA],[65] and Ligand Lipophilic
Efficiency (LLE) [LLE = pKd–Log D][65] On the Basis of ITC Results,
PAMPA Measured Passive Permeability, Chromatographic Hydrophobicity
Index log D (CHILogD7.4), Number of Hydrogen
Bond Donors, Calculated Topological Polar Surface Area (TPSA), and
Number of Rotatable Bonds of Inhibitors 1, 3, 6, and 15–19
Permeability of
inhibitors were
determined at room temperature by PAMPA, using propanolol (medium
permeability −62.7 nm/s) as a control.
Values calculated using StarDrop
software.
Figure 6
Inhibitor 15 is not cytotoxic.
Colony formation (A)
and cell-proliferation assays (B) showing no toxicity in HeLa cells
after treatment with 150 μM of inhibitor 15 and
its cis epimer.
HeLa cells were treated with 50
μM inhibitors and HIF-1α stabilization levels were measured
by Western blot after 2 h treatment (protein levels normalized to
HIF-1α stabilization level observed with 1).Inhibitor 15 is not cytotoxic.
Colony formation (A)
and cell-proliferation assays (B) showing no toxicity in HeLa cells
after treatment with 150 μM of inhibitor 15 and
its cis epimer.To support the SAR results of the second inhibitor series,
the
X-ray crystal structures of 15 and 16–19 bound to VBC were obtained. These crystal structures revealed
that most of the key interactions and structural features observed
for 10 were maintained. Described herein are the key
changes in interactions observed at the LHS from the different groups.
The cyano group on inhibitor 15 (Figure A) forms a hydrogen bond with a water molecule
that allows the formation of a water network similar to the one observed
in the crystal structure with inhibitor 3 bound (Figure A). In the crystal
structure of inhibitor 16 (Figure B), the ketone adopts an optimal distance
(2.5 ± 0.1 Å) and acceptable orientation (angles formed
between the C=O bond and amide nitrogen of 90.7 ± 2.0°
and between the C–N bond and ketone oxygen of 91.5 ± 0.8°)
for the formation of an intramolecular hydrogen bond with the NH of
inhibitor LHS amide (Supporting Information, Figure 5).[60] The observed intramolecular
interaction could stabilize the conformation of the LHS of the inhibitor,
improving affinity and cell permeability, thus contributing to the
observed high cellular activity.[61,62]
Figure 7
Co-crystal
structures of second-series inhibitors. Crystal structure
of VBC in complex with 15 (A, PDB 5LLI), 16 (B), 17 (C), 18 (D), and 19 (E) (purple carbons). VHL is shown as a pale-green surface and the
VHL residues forming the binding pocket as orange stick representations.
Inhibitor-binding waters are shown as red spheres. Hydrogen bond interactions
between inhibitors, bound waters, and VHL pocket residues are shown
as dashed red lines.
Co-crystal
structures of second-series inhibitors. Crystal structure
of VBC in complex with 15 (A, PDB 5LLI), 16 (B), 17 (C), 18 (D), and 19 (E) (purple carbons). VHL is shown as a pale-green surface and the
VHL residues forming the binding pocket as orange stick representations.
Inhibitor-binding waters are shown as red spheres. Hydrogen bond interactions
between inhibitors, bound waters, and VHL pocket residues are shown
as dashed red lines.The crystal structure of 17 (Figure C) did not show any major changes
in the
ligand binding mode and in the conformation of pocket residues when
compared with the crystal structure of the fluoro-analogue 10 (Figure C). The
terminal acetamide of 17 points upward toward the solvent,
not forming any new interactions with the protein nor an intramolecular
hydrogen bond with the NH of the LHS inhibitor amide. These features,
together with the extra amide H-bond donor group, likely account for
the poor cellular activity of 17 (Table ).The structure of inhibitor 18 bound to VBC did not
reveal any critical changes in either the ligand or pocket residue
conformations (Figure D). The cyclobutyl moiety fills the far LHS of the pocket, inducing
a conformational change to the side chain of Arg69, as previously
described.[41] Replacement of the cyclobutyl
with an oxetane group in inhibitor 19 showed that the
ether oxygen was able to replace the LHS amide in the ability to form
a hydrogen bond with the LHS structural water (Figure E). A hydrogen bond interaction from the
oxetane oxygen (2.9 ± 0.1 Å distance) leads to an increase
in the distance between the amide carbonyl oxygen and the structural
water (3.7 ± 0.2 Å) (Supporting Information, Figure 6). This is incompatible with the bidentate
interaction mode previously observed with inhibitor 3, a result of the increased distance between the two oxygen atoms
in 19 relative to 3 and the steric constraints
of the oxetane ring.In summary we describe a group-based optimization
of 1 that led us to the discovery of a more potent inhibitor 10. Further optimization led us to the discovery of inhibitor 15 which showed equivalent cellular potency when compared
with 10 and was not cytotoxic in cells.
Chemistry
Intermediate 22, 24, and 27 were obtained following the protocol previously
described by us.[41,55] Intermediate 22 was
used for the synthesis of final inhibitors 4 and 6–8 via acylation or inhibitors 3, 5, 9, 10, and 16–20 by amide coupling reaction (Scheme ). Further treatment
of intermediate 22 with 4-chlorobutanoyl chloride, followed
by intramolecular cyclization, led to the synthesis of inhibitor 12 (Scheme ). Reductive amination of 22 with formaldehyde followed
by acylation yielded inhibitor 13 (Scheme ). An amide coupling reaction between intermediate 22 and 1-(tert-butoxycarbonyl)azetidine-3-carboxylic
acid led to the intermediate 23, which after deprotection,
followed by acylation, yielded inhibitor 21 (Scheme ).
Reagents and conditions:
(i)
TFA:DCM, rt, 30 min; (ii) Boc-l-Ala, HATU, DIPEA, DMF, rt,
1 h; (iii) TFA:DCM, rt, 30 min; (iv) N(Et)3, acetic anhydride,
DCM, rt, 2 h; (v) Boc-l-Thr, HATU, DIPEA, DMF, rt, 1 h; (vi)
1-fluorocyclopropanecarboxylic acid, HATU, DIPEA, DMF, rt, 1 h.The cis epimer of inhibitor 10 was
obtained from intermediate 27 following the same synthetic
pathway used for its trans conformer, using Boc-Hyp cis epimer to obtain the desired final inhibitor (Scheme ).
Scheme 3
Reagents and conditions: (i)
TFA:DCM, rt, 30 min; (ii) 1-fluorocyclopropanecarboxylic acid, HATU,
DIPEA, DMF, rt, 1 h.
Reagents and conditions: (i)
TFA:DCM, rt, 30 min; (ii) 1-fluorocyclopropanecarboxylic acid, HATU,
DIPEA, DMF, rt, 1 h.
Evaluation of VHL Inhibitor
Permeability and Binding Kinetics
Cell membrane passive permeability
is known to have an important
role in inhibitor optimization, as it can lead to an increase in intracellular
concentration of inhibitor available to interact with the desired
target, which would be expected to increase cellular activity. Optimization
of lipophilicity at physiological pH, topological polar surface area
(TPSA), and number of rotatable bonds are known to directly impact
cell membrane passive permeability.[63,64] From the SAR
results, it was possible to observe that the increase in cellular
potency was not only related with the increase of inhibitor binding
affinity to VBC (Figure A). Therefore, to evaluate how lipophilicity and cell membrane passive
permeability contributed to VHL inhibitors activity in cells, these
parameters were experimentally measured for a selection of the most
potent inhibitors (Table ).
Figure 8
Correlations between cellular activity and physicochemical
parameters
in vitro. Reported values are from Tables and 3 for inhibitors 1, 3, 6, 10, 15, and 16–19. (A) ITC pKd vs activity in cells. (B) CHILogD7.4 vs PAMPA permeability. (C) relative cellular activity vs permeability.
HIF-1α protein levels measured using Western blot, as proxy
of cellular activity, and normalized to protein levels with inhibitor 10.
Correlations between cellular activity and physicochemical
parameters
in vitro. Reported values are from Tables and 3 for inhibitors 1, 3, 6, 10, 15, and 16–19. (A) ITC pKd vs activity in cells. (B) CHILogD7.4 vs PAMPA permeability. (C) relative cellular activity vs permeability.
HIF-1α protein levels measured using Western blot, as proxy
of cellular activity, and normalized to protein levels with inhibitor 10.Permeability of
inhibitors were
determined at room temperature by PAMPA, using propanolol (medium
permeability −62.7 nm/s) as a control.Values calculated using StarDrop
software.Inhibitors 10 (LE = 0.27, LLE = 5.73), 15 (LE = 0.26, LLE
= 5.65), and 16 (LE = 0.25, LLE = 5.63)
with permeabilities of 27.1, 19.4, and 15.0 nm/s, respectively, were
found to be the most permeable inhibitors. At the other extreme were
inhibitors 3 (LE = 0.27, LLE = 6.28), 17 (LE = 0.24, LLE = 6.28), and 19 (LE = 0.25, LLE = 5.92)
with the lowest permeability values. A direct relationship was observed
between inhibitor permeability, lipophilicity, and activity in cells
(Figure B,C). Inhibitors 6, 10, 15, 16, and 18 were found to be more lipophilic than starting inhibitor 1, contributing to their higher cellular permeability, which
was thus reflected in their reasonable intracellular free compound
concentration (Supporting Information, Table 1) that leads to a higher cellular potency.[66] In contrast, analysis of TPSA, number of rotatable bonds, and number
of hydrogen bond donors (HBDs) revealed no correlation with cell permeability
(Table ). Because
of limitation of our PPI-targeting pharmacophore, these parameters
in VHL inhibitors closely approach the maximal limits for desired
membrane permeability (TPSA < 140 Å2, HBDs <
5).[64,67] Clearly, for the most active inhibitors,
these high values are compensated for by a higher lipophilicity, resulting
in an overall increase of inhibitors permeability and activity in
cells. The data, however, suggests that careful monitoring of these
parameters is warranted. For example, adding an extra HBD group (from
3 to 4) in inhibitors 3 and 17 led to marked
decrease in permeability and cellular potency compared to 1, despite the comparable CHILogD7.4 values (Tables and 3). Similarly, breaking beyond TPSA of 140 Å2 in inhibitor 17 resulted in the lowest cell permeability and cellular activity
amonst the tested inhibitors (Tables and 3). These results argue
for close monitoring of all these physicochemical properties during
optimization of a PPI inhibitor.In addition to permeability
and binding affinity, binding kinetics
are also key parameters that can influence cellular activity. Longer
drug residence at the target can promote the duration of pharmacological
on-target effect. To evaluate how the binding kinetics could influence
VHL inhibitor activity in cells, binding kinetics parameters were
determined by surface plasmon resonance for 1, 10, 15, 16, and 18,
which all showed high passive membrane permeability and high potency
in cells (Table ). The binding experiments were performed
at 10 and 20 °C in order to compare the kinetics and affinity
of the inhibitors at different temperatures.
Table 4
Structures,
Association (kon), and Dissociation (koff) Rate Constants Determined by Surface Plasmon
Resonance (SPR), Kds Determined by SPR
(Kd = koff/kon), Calculated Dissociative Half-Times
(t1/2), and Microsomal and Plasma Stability
of Inhibitors 1, 10, 15, 16, and 18
Dissociative half-life of inhibitor–VBC
binary complex calculated based on inhibitor dissociation rate constants
(koff), t1/2 = ln2/koff.[68]
Dissociative half-life of inhibitor–VBC
binary complex calculated based on inhibitor dissociation rate constants
(koff), t1/2 = ln2/koff.[68]A close relationship was
observed between dissociation rate constants
and dissociation constants from VBC. Inhibitors 10, 15, and 16, with higher binding affinity than 1, also showed slower dissociation rates relative to 1 (Table ).
After comparison of the dissociative half-life (from VBC complex)
of these inhibitors, it was possible to observe higher dissociation
times correlate with higher cellular potencies. Inhibitor 10 and 15, the most active inhibitors in cells, had the
highest dissociation half-lives among the compounds evaluated (Figure ). It is worth pointing
out that even the slowest off rates measured for our VHL inhibitors
are still relatively fast overall. The resulting residence times are
consequently all on the order of seconds, rather than minutes, suggesting
that there is much scope to optimize this parameter in the future.
Nevertheless, our data suggests that even these relatively short residence
times suffice for a pharmacological effect on VHL, as shown with other
drug–target interactions.[69]
Figure 9
Activity in
cells vs half-life dissociation from VBC protein complex.
Values shown are for inhibitors 1, 3, 6, 10, 15, and 16–19 at 10 °C (squares) and at 20 °C (triangles).
HIF-1α protein levels measured using Western blot, as proxy
of cellular activity, and normalized to protein levels with inhibitor 10.
Activity in
cells vs half-life dissociation from VBC protein complex.
Values shown are for inhibitors 1, 3, 6, 10, 15, and 16–19 at 10 °C (squares) and at 20 °C (triangles).
HIF-1α protein levels measured using Western blot, as proxy
of cellular activity, and normalized to protein levels with inhibitor 10.To begin to assess potential
metabolic liabilities and inform potential
use on animal models, we evaluated the microsomal and plasma stability
of the most promising VHL inhibitors (Table ). The data revealed slow microsomal clearance
(0.7–3.1 mL × min–1 × g–1) and high plasma metabolic stability (inhibitors half-life >180
min) with all inhibitors, supporting their use as chemical probes
in vivo in animal models.Overall, the data support 15 as optimal inhibitor,
with high binding affinity and cellular potency, good cell membrane
permeability, slow dissociation from its target, low toxicity in cells,
and excellent stability. On the basis of these results, 15 was elected as a novel chemical probe selectively blocking VHL in
the hypoxia-signaling pathway downstream of HIF hydroxylation.[41] The compound is available on the “Chemical
Probes Portal” (http://www.chemicalprobes.org).[59]
Conclusions
We
report the structure-guided optimization of a new series of
ligands targeting the VHL E3 ubiquitin ligase. The best inhibitors
have shown affinities in the nanomolar range and have on-target cellular
activity. An increase in inhibitor lipophilicity yielded a higher
cell membrane passive permeability, which together with slower dissociation
rate constants (koff) and increased inhibitor/VBC
dissociation half-times (t1/2), clearly
led to increased cellular potency. Additionally, the most promising
inhibitors have also shown excellent microsomal and plasma stability,
supporting their use in vivo with animal models. Our best inhibitors 10 and 15, to our knowledge, are the first inhibitors
to reach double-digit nanomolar affinities for binding to VHL, as
well as the most potent inhibitors of the VHL:HIF-1α PPI inside
cells reported to date. Inhibitor 15 was elected as the
final chemical probe (http://www.chemicalprobes.org/vh298). The compound, and its
inactive cis epimer, are now available from various commercial vendors.Our work revealed new VHL inhibitors that can act as probes of
the hypoxia signaling pathway, an approach recently validated by us.[41] Future optimization of in vivo activity is warranted,
which could provide lead compounds with an alternative mode of action
as potential therapeutics against diseases where mimicking a hypoxic
response has proven to be beneficial. Additionally, the reported novel
VHL ligands can be explored for the design of improved VHL-recruiting
bifunctional chimeric molecules (PROTACs) to efficiently hijack the
VHL E3 ubiquitin ligase to induce intracellular degradation of target
proteins.[44,46,49] Importantly,
we provide an exemplary study for optimizing PPI-targeting E3 ligase
ligands via similar strategies, which could inspire current and future
efforts to develop inhibitors or PROTACs against other E3 ligases.
Experimental Section
Materials
and Methods
Commercially available starting
reagents for each reaction were purchased from Sigma-Aldrich, Fluorochem,
Apollo Scientific, or Manchester Organics and used without further
purification. All reactions were carried out using anhydrous solvents.
Analytical thin-layer chromatography (TLC) was performed on precoated
TLC plates (layer 0.20 mm silica gel 60 with fluorescent indicator
(UV 254: Merck)). The TLC plates were air-dried and revealed under
UV lamp (254/365 nm). Flash-column chromatography was performed using
prepacked silica gel cartridges (230–400 mesh, 40–63
mm; SiliCycle) using a Teledyne ISCO Combiflash Companion or Combiflash
Retrieve using the solvent mixtures stated for each synthesis as mobile
phase.Liquid chromatography–mass spectrometry (LC-MS)
analyses were performed with either an Agilent HPLC 1100 series connected
to a Bruker Daltonics MicroTOF or an Agilent Technologies 1200 series
HPLC connected to an Agilent Technologies 6130 quadrupole spectrometer
or a Waters 2795 connected to a Waters ZQ Micromass spectrometer,
where all instruments were connected to a diode array detector. All
the final compounds used in all the experiments were evaluated after
preparative LC-MS separations with a Waters X-bridge C18 column (50
mm × 2.1 mm × 3.5 mm particle size); flow rate, 0.5 mL/min
with a mobile phase of water/MeCN + 0.1% CHOOH or water/MeCN + 0.1%
NH3; 95/5 water/MeCN was initially held for 0.5 min followed
by a linear gradient from 95/5 to 5/95 water/MeCN over 3.5 min which
was then held for 2 min. The purity of all the compounds was evaluated
using the analytical LC-MS system described before and yield a purity
>95%.High-resolution electrospray measurements were performed
on a Bruker
Daltonics MicroTOF mass spectrometer. 1H NMR and 13C NMR spectra were recorded on a Bruker Avance II 500 spectrometer
(1H at 500.1 MHz, 13C at 125.8 MHz) or on a
Bruker DPX-400 cryospectrometer (1H at 400.1 MHz, 13C at 101 MHz). Chemical shifts (δ) are expressed in
ppm reported using residual solvent as the internal reference in all
cases. Signal splitting patterns are described as singlet (s), doublet
(d), triplet (t), multiplet (m), or a combination thereof. Coupling
constants (J) are quoted to the nearest 0.1 Hz.Intermediates 22, 24, and 27, and final inhibitors 1 and 15 were synthesized
as described elsewhere.[41,55]
General Methodology
for the Synthesis of VHL Inhibitors 3–11, 13–21, cis-10, and Intermediates 23, 25, and 26
General Method A (Synthesis by Acylation)
A solution
of compound 2 (100 mg, 0.19 mmol) in 1:1 TFA:DCM (6 mL)
was stirred at room temperature for 30 min. The solvents were evaporated
under reduced pressure to give the corresponding deprotected intermediate
(TFA salt–22) as a brown oil that was used in
the following reactions without further purification (102 mg, 0.19
mmol). To a solution of the deprotected intermediate 22 (102 mg, 0.19 mmol, 1 equiv) in DCM was added triethylamine (57
mg, 79 μL, 0.57 mmol, 3 equiv). After stirring the mixture for
10 min at room temperature, acetic anhydride derivative (1.5 equiv)
was added and the resulting mixture was then stirred 2 h at room temperature.
The solvents were evaporated under reduced pressure to afford the
corresponding crude compound that was purified by flash column chromatography
using a gradient of 10% to 70% acetone in heptane to yield the final
compounds as solids.
General Method B (Synthesis by HATU Assisted
Amide Coupling)
A solution of compound 2 (100
mg, 0.19 mmol) in 1:1
TFA:DCM (6 mL) was stirred at room temperature for 30 min. The mixture
was evaporated under reduced pressure to give the corresponding deprotected
intermediate (TFA salt–22) as a brown oil without
further purification (102 mg, 0.19 mmol). To a solution of the deprotected
intermediate 22 (102 mg, 0.19 mmol, 1 equiv) in DMF was
added the carboxylic acid derivative (1 equiv). DIPEA (97 mg, 129
μL, 0.75 mmol, 4 equiv) was added dropwise, and the mixture
was stirred for 5 min at room temperature. HATU (78 mg, 0.21 mmol,
1.1 equiv) was added, and the mixture was stirred at room temperature
for 1 h. Water was added, and the mixture was extracted with ethyl
acetate (3×). The combined organic phases were washed with brine
(2×), dried over MgSO4, and evaporated to afford the
corresponding crude compound that was purified by flash column chromatography
using a gradient of 10% to 70% acetone in heptane to yield the final
compounds as solids.
To a solution of intermediate 22 (205 mg, 0.38 mmol) and 4-chlorobutanoyl chloride (0.069 mg, 56
μL, 0.49 mmol, 1.3 equiv) in DCM (15 mL), a cold 1N solution
of NaOH (6 mL) was added and the resulting mixture was stirred vigorously
for 20 min at room temperature. The organic layer was collected and
dried over MgSO4 and evaporated to afford the intermediate
crude compound. To a solution of the intermediate crude compound in
THF cooled to 0 °C was added potassium tert-butoxide
(128 mg, 1.14 mmol, 3 equiv), and the mixture was allowed to warm
to room temperature and stirred overnight. Water was added, and the
mixture was extracted with ethyl acetate (3×). The combined organic
phases were dried over MgSO4 and evaporated to afford the
corresponding crude compound that was purified by flash column chromatography
using gradient elution of 10–70% acetone in heptane to yield
the final compound 12 as a white powder (120 mg, 0.24
mmol, 63%). 1H NMR (CDCl3, 500 MHz): δ
8.68 (s, 1H), 7.37–7.32 (m, 5H), 4.74 (s, 1H), 4.70 (t, 1H, J = 10.0 Hz), 4.56–4.51 (m, 2H), 4.38 (dd, 1H, J = 15.0, 5.0 Hz), 4.01 (d, 1H, J = 10
Hz), 3.75–3.64 (m, 2H), 3.58 (dd, 1H, J =
10.0, 5.0 Hz), 2.58–2.51 (m, 4H), 2.41–2.35 (m, 1H),
2.31–2.25 (m, 1H), 2.11–2.0 (m, 1H), 1.99–193
(m, 2H), 1.00 (s, 9H). 13C NMR (CDCl3, 125 MHz):
δ 177.0, 171.9, 170.2, 150.5, 148.6, 138.2, 131.7, 131.1, 129.7,
128.2, 70.0, 59.2, 58.3, 56.6, 47.2, 43.4, 35.9, 35.8, 30.7, 27.7,
19.0, 16.2. HRMS (ESI) m/z: [M+ + 1] calculated for C26H35N4O4S, 499.2379; observed, 499.2382. HPLC: tR = 2.9 min, k′ = 8.7.
Ethanol (2 mL)
was added to a round-bottom flask and cooled to 0 °C, then thionyl
chloride (706 mg, 430 μL, 5.94 mmol, 2 equiv) was added dropwise
and the mixture was stirred for 15 min. 1-Aminocyclopropanecarboxylic
(300 mg, 2.97 mmol, 1 equiv) acid was added, and the mixture was refluxed
for 2 h. Solvent was removed, and the final compound 5 was obtained as a white powder (363 mg, 2.82 mmol, 95%). 1H NMR (DMSO-d6, 400 MHz): δ 9.11
(s, 2H), 4.16 (q, 2H, J = 12.0, 4.0 Hz), 1.51–1.47
(m, 2H), 1.39–1.36 (m, 2H), 1.21 (t, 3H, J = 8.0 Hz). 13C NMR (DMSO-d6, 101 MHz): δ 170.0, 62.3, 34.0, 14.4, 13.6.
Ethyl 1-Acetamidocyclopropanecarboxylate
Following
the general protocol for the synthesis of final inhibitors by acylation,
from intermediate 5 (100 mg, 0.77 mmol, 1 equiv) and
acetic anhydride (118 mg, 110 μL, 1.15 mmol, 1.5 equiv), compound 6 was obtained as a pale-lime powder (80 mg, 0.47 mmol, 62%). 1H NMR (DMSO-d6, 400 MHz): δ
8.44 (s, 1H), 4.02 (q, 2H, J = 12.0, 4.0 Hz), 1.79
(s, 3H), 1.33–1.30 (m, 2H), 1.14 (t, 3H, J = 8.0 Hz), 0.98–0.95 (m, 2H). 13C NMR (DMSO-d6, 101 MHz): δ 172.1, 170.2, 60.5, 32.8,
22.4, 16.5, 14.1.
1-Acetamidocyclopropanecarboxylic Acid
To a solution
of compound 6 (70 mg, 0.41 mmol, 1 equiv) in methanol
was added an aqueous solution of NaOH (1N), and the resulting mixture
was then heated at 100 °C for 4 h. Methanol was evaporated, and
the resulting solution was acidified and extracted with DCM (5×).
The combined organic phases were dried over MgSO4 and evaporated
to afford the corresponding final compound 7 as a white
powder (53 mg, 0.37 mmol, 91%). 1H NMR (DMSO-d6, 400 MHz): δ 12.25 (s, 1H), 8.35 (s, 1H), 1.77
(s, 3H), 1.30–1.28 (m, 2H), 0.94–0.91 (m, 2H). 13C NMR (DMSO-d6, 101 MHz): δ
174.0, 170.1, 32.5, 22.3, 16.3.
FP competitive binding
experiments were performed on a PHERAstar FS (BMG LABTECH) in 384-well
plates (Corning 3575), with an excitation wavelength (λ) at
485 nm and emission λ at 520 nm. Each well solution (15 μL)
contained 15 nM of VBC protein, 10 nM of FAM-labeled HIF-1α
peptide (FAM-DEALAHypYIPMDDDFQLRSF, Kd = 3 nM as measured by a direct FP titration), and decreasing concentrations
of compound (14-point serial 2-fold dilutions starting from 50 μM)
in 100 mM Bis-tris, 100 mM NaCl, 1 mM DTT, pH 7. Control wells contained
VBC and peptide in the absence of compound (maximum signal) and peptide
in the absence of protein (background signal). Data were obtained
in triplicate, and the percentage of displacement was determined and
graphed against log[VHL inhibitors]. Average IC50 values
and the standard error of the mean (SEM) were determined for each
titration using Prism 6. Dissociation constants Kd were back-calculated from the measured IC50 values using a displacement binding model, as described previously.[17]
Isothermal Titration Calorimetry (ITC)
ITC experiments
were carried in an ITC200 microcalorimeter (GE Healthcare). The compounds
were diluted from DMSO stock solution to 300 μM in a buffer
of 20 mM Bis-Tris propane, 150 mM NaCl, 1 mM DTT, pH 7. The compounds
were titrated against 30 μM VBC complex, equilibrated in the
same buffer. The final concentration of DMSO in each experiment was
3% (v/v). The titrations consisted of 20 injections of 2 μL
of ligand solution at a rate of 0.5 μL/s at 120 s time intervals.
An initial injection of ligand (0.4 μL) was made and discarded
during data analysis. All experiments were performed at 25°C while stirring at 750 rpm. The data were fitted to a single-binding-site
model using the Microcal LLC ITC200 Origin software provided by the
manufacturer to obtain the stoichiometry n, the dissociation
constant Kd, and the enthalpy of binding
ΔH.
Cell Culture and Cell Treatment
Human cervical carcinoma
cell HeLa was obtained from ATCC and propagated in DMEM supplemented
with 10% fetal bovine serum (FBS), l-glutamine, and 100 μg/mL
of penicillin/streptomycin at 37 °C. HeLa cells were routinely
tested for mycoplasma contamination using MycoAlert kit from Lonza.
Cells were treated with VHL inhibitors at 50 μM for 2 h in fresh
medium, and 1% DMSO was used as vehicle control.
Immunoblotting
Cells were lysed in RIPA buffer (50
mM Tris pH 8, 150 mM NaCl, 1% NP-40, 0.5% sodium deoxycholate, 0.1%
SDS, 250 M Na3VO4, 10 mM NaF) and a protease
inhibitor cocktail (Roche) per 10 mL of buffer. Proteins were resolved
using sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS–PAGE),
transferred onto nitrocellulose membranes, and detected using primary
antibodies, with β-actin as loading control. Primary antibodies
were used at following dilutions for mammalian cells: anti-HIF-1α
(BD Biosciences, 610958, 1:1000), antihydroxy-HIF-1α (Hyp564)
(Cell Signaling Technology, no. 3434, 1:1000), and anti-β-actin
(Cell Signaling Technology, no. 3700s, 1:10000). Following primary
antibody incubation, a goat antimouse (LI-COR, 926-32210, 1:10,000)
or donkey antirabbit (LI-COR, 926–32213, 1:10000) conjugated
to IRDye 800CW secondary antibody were used for detection on the LI-COR
Odyssey (LI-COR, Bad Homburgh, Germany).
Proliferation Assay
First, 1.5 × 105 cells were seeded in six-well plates
1 day prior to treatments with
VHL compounds. At treatment times of 0, 24, 48, and 72 h, cells were
then trypsinized using 400 μL of trypsin, followed by 600 μL
of media. Cells were then counted using hemocytometer.
Colony Formation
Assay
First, 500–2000 cells
were seeded in a six-well plate 1 day prior to treatments with VHL
compounds. After 24 h of treatment, media was then changed and cells
were further incubated for 1 week. Cells were washed with PBS twice
and stained with 0.005% (w/v) crystal violet in 25% (v/v) methanol
for 10 min at room temperature. The staining was washed off with water,
and plates were left to dry overnight. Numbers of colonies formed
were counted with M-tools suite from OMERO on the scanned plates.
X-ray Crystallography
The VBC ternary complex was purified
and crystallized as described previously.[17,52] Equal volume solutions of VBC (∼5 mg/mL) and liquor solution
were mixed in the hanging-drop vapor diffusion method at 18 °C.
The liquor solution contained 0.1 mM sodium cacodylate, pH 6.2–6.6,
16–18% polyethylene glycol 3350, 0.2 M magnesium acetate, and
10 mM DTT. The drop was streaked with seeds of disrupted VBC crystals,
and a 2–3 mm layer of Al’s Oil (Hampton Research) was
applied on top of the liquor solution to slow the vapor diffusion
rate. To obtain the structures of VHL inhibitors bound to VBC, crystals
were soaked overnight in a 1 mM solution of inhibitor in 1% DMSO,
4% 2-propanol, and 95% liquor solution. Crystals were screened using
an in-house Rigaku M007HF X-ray generator and Saturn 944HG+ CCD detector.
X-ray data were collected at 100 K at Diamond Light Source beamline
I04–1. Indexing and integration of reflections was performed
using XDS with the XDSGUI interface[70] and
scaling and merging with AIMLESS in CCP4i.[71,72] The isomorphous data sets were refined using REFMAC5 (refs (73,74)) and COOT (ref (75)) using a template structure derived from the
Protein Data Bank (PDB) entry 1vcb (ref (76)). Ligand structures and restraints were generated
using the PRODRG server.[77] The MOLPROBITY
server was used to validate the geometry and steric clashes in the
structures.[78] The structures have been
deposited in the PDB with accession codes 5NVV, 5NVW, 5NVX, 5NVY, 5NVZ, 5NW0, 5NW1, and 5NW2, and data collection and refinement statistics
are presented in Supporting Information, Table 1.
PAMPA was performed using a 96-well precoated
BD Gentest PAMPA plate
(BD Biosciences, UK). Each well was divided into two chambers: donor
and acceptor, separated by a lipid–oil–lipid trilayer
constructed in a porous filter. The effective permeability, Pe, of the compound was measured at pH 7.4. Stock
solutions (5 mM) of the compound were prepared in DMSO. The compound
was then further diluted to 10 μM in PBS, pH 7.4. The final
DMSO concentration did not exceed 5% (v/v). The compound dissolved
in PBS was then added to the donor side of the membrane and PBS without
compound was added to the acceptor side. The PAMPA plate was left
at room temperature for 5 h, after which time, an aliquot (100 μL)
was removed from both acceptor and donor compartments and mixed with
acetonitrile (80 μL) containing an internal standard. The samples
were centrifuged (10 min, 5 °C, 3270 g) to sediment precipitated
protein and then sealed prior to UPLC-MS/MS analysis using a Quattro
Premier XE (Waters Corp, USA). Pe was
calculated as shown in the equation below:where: CA(t) = peak area of compound
present in acceptor well at time t = 18000 s, Cequiv = [CD(t) × VD + CA(t)
× VA]/(VD + VA), CD(t) = peak area of compound present in donor well
at time t = 18000 s, A = filter
area, VD = donor-well volume, VA = acceptor-well volume, t = incubation time. Recovery of compound from donor and acceptor
wells was calculated, and data was only accepted when recovery exceeded
70%.
CHILogD7.4 Measurement
The CHIlogD (chromatographic
hydrophobicity index log D) at pH 7.4 was determined
using retention time measurements on a HPLC Dionex system (Thermo
Fisher) with a Luna C18 column (Phenomenex). Test samples in DMSO
(10 mM) were diluted to a concentration of 0.25 mM using 50:50 acetonitrile:water.
Mobile phase A was 10 mM ammonium acetate solution (pH 7.4), and mobile
phase B was acetonitrile. HPLC method was as follows: 1 mL/min flow,
temperature 20 °C, injection volume 10 μL, gradient 0–10.5
min 100% A, 10.5–14 min 100% B, 14–15 min 100% A. A
calibration line was generated using a test mix of compounds (paracetamol,
theophylline, caffeine, benzimidiazole, colchicine, carbamazepine,
indole, propiophenone, butyrophenone, valerophenone, and heptanophenone).
The CHIlogD was calculated as previously described.[79,80]
Intracellular Drug Concentration.[66]
Briefly, HeLa cells (1 × 106 cells/2 mL
per well, 6-well plate) were incubated with known concentration of
compounds (50 μM) for 2 h at 37 °C before being washed
and lysed using lysis buffer (20 mM Tris, 150 mM, NaCl, 1% Triton
X100, and 1 pill per 10 mL of proteases inhibitors (Roche)). The lysate
was resuspended in 2 mL of ice-cold PBS. Resuspended lysate (100 μL)
was then subjected to solvent crash in a 1:2 ratio of lysate to acetonitrile
containing internal standard 5 ng/mL of donepezil. The concentration
was determined with the aid of an appropriate calibration curve and
UPLC-MS/MS. The unbound intracellular compound concentration (free
fraction) was determined by dialyzing 150 μL of resuspended
cell lysate against isotonic phosphate buffer in an equilibrium dialysis
equipment. This was also subjected to UPLC-MS/MS.
Surface Plasmon
Resonance (SPR)
VHL inhibitors were
dissolved in DMSO (1 mM) and then diluted 20-fold in DMSO to achieve
a 50 μM final stock concentration. Ligand stock solution was
serially diluted 2-fold (five times) in DMSO, and the obtained solutions
were then diluted individually in SPR buffer (20 mM HEPES, 150 mM
NaCl, 1 mM DDT, 0.005% Tween P20, pH 7.0) to obtain the final 2% (v/v)
DMSO concentration series from 1 μM to 31.25 nM (2-fold dilutions)
and transferred to a 96-well plate. The experiment was conducted in
a Biacore T100 (GE Healthcare, Biacore, Uppsala, Sweden) at 10 and
20 °C, and solutions were injected individually using 60 and
160 s association and dissociation times, respectively. Data were
treated using Biacore T100 evaluation software provided by the manufacturer.
Reference flow-cell response was subtracted from the sample response
with immobilized VBC protein to correct for systematic noise and baseline
drift. Data were solvent corrected, and the response from the blank
injections was used to double reference the binding data. The data
were normalized by molecular weight, and rate constants kon and koff values were obtained
using a 1:1 binding model fit.
Intrinsic Clearance (Cli)
Experiments
Test compound
(0.5 μM) was incubated with female CD1 mouse liver microsomes
(Xenotech LLC; 0.5 mg/mL 50 mM potassium phosphate buffer, pH 7.4)
and the reaction started with addition of excess NADPH (8 mg/mL 50
mM potassium phosphate buffer, pH 7.4). Immediately, at time zero,
then at 3, 6, 9, 15, and 30 min, an aliquot (50 μL) of the incubation
mixture was removed and mixed with acetonitrile (100 μL) to
stop the reaction. Internal standard was added to all samples, the
samples were centrifuged to sediment precipitated protein, and the
plates then sealed prior to UPLCMSMS analysis using a Quattro Premier
XE (Waters Corporation, USA).XLfit (IDBS, UK) was used to calculate
the exponential decay and consequently the rate constant (k) from the ratio of peak area of test compound to internal
standard at each time point. The rate of intrinsic clearance (CLi) of each test compound was then calculated using the following
calculation:where V (mL/mg
protein) is
the incubation volume/mg protein added and microsomal protein yield
is taken as 52.5 mg protein/g liver. Verapamil (0.5 μM) was
used as a positive control to confirm acceptable assay performance.
Plasma Stability (Plas Stab) Experiments
Test compound
(50 μM) was incubated in prewarmed plasma at 37 °C (that
is buffered to pH 7.4 in ratio of 70:30 plasma to buffer). Immediately,
at time zero, then at 30, 60, 120, and 180 min, a 50 μL aliquot
of the incubation mixture was removed and mixed with 200 μL
of acetonitrile containing Donepezil as the internal standard (50
ng/mL) to stop the reaction. The samples were centrifuged to sediment
the precipitated protein and the plates then sealed prior to UPLC-MS/MS
analysis using a Quattro Premier XE (Waters Corporation, USA).XLfit (IDBS, UK) was used to calculate the exponential decay and
consequently the rate constant (k) from the ratio
of peak area of test compound to internal standard at each time point.
The half-life was calculated for each test compound from the rate
by using the following calculation:
Authors: Garib N Murshudov; Pavol Skubák; Andrey A Lebedev; Navraj S Pannu; Roberto A Steiner; Robert A Nicholls; Martyn D Winn; Fei Long; Alexei A Vagin Journal: Acta Crystallogr D Biol Crystallogr Date: 2011-03-18
Authors: Martyn D Winn; Charles C Ballard; Kevin D Cowtan; Eleanor J Dodson; Paul Emsley; Phil R Evans; Ronan M Keegan; Eugene B Krissinel; Andrew G W Leslie; Airlie McCoy; Stuart J McNicholas; Garib N Murshudov; Navraj S Pannu; Elizabeth A Potterton; Harold R Powell; Randy J Read; Alexei Vagin; Keith S Wilson Journal: Acta Crystallogr D Biol Crystallogr Date: 2011-03-18
Authors: Matthias R Bauer; Paolo Di Fruscia; Simon C C Lucas; Iacovos N Michaelides; Jennifer E Nelson; R Ian Storer; Benjamin C Whitehurst Journal: RSC Med Chem Date: 2021-01-07
Authors: Haibin Zhou; Longchuan Bai; Renqi Xu; Yujun Zhao; Jianyong Chen; Donna McEachern; Krishnapriya Chinnaswamy; Bo Wen; Lipeng Dai; Praveen Kumar; Chao-Yie Yang; Zhaomin Liu; Mi Wang; Liu Liu; Jennifer L Meagher; Han Yi; Duxin Sun; Jeanne A Stuckey; Shaomeng Wang Journal: J Med Chem Date: 2019-12-10 Impact factor: 7.446