Martin G Strebl1,2, Arthur J Campbell3, Wen-Ning Zhao4,5, Frederick A Schroeder1, Misha M Riley1, Peter S Chindavong4,5, Thomas M Morin1,6, Stephen J Haggarty4,5, Florence F Wagner3, Tobias Ritter2,7,8, Jacob M Hooker1,7. 1. Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital and Harvard Medical School, Charlestown, Massachusetts 02129, United States. 2. Department of Chemistry and Chemical Biology, Harvard University, 12 Oxford Street, Cambridge, Massachusetts 02138, United States. 3. Stanley Center for Psychiatric Research, Broad Institute of MIT and Harvard, 75 Ames Street, Cambridge, Massachusetts 02142, United States. 4. Chemical Neurobiology Laboratory, Center for Genomic Medicine, Department of Neurology, Massachusetts General Hospital, Harvard Medical School, Boston, Massachusetts 02114, United States. 5. Department of Psychiatry, Massachusetts General Hospital, Harvard Medical School, Boston, Massachusetts 02114, United States. 6. Tufts University, 419 Boston Avenue, Medford, Massachusetts 02155, United States. 7. Division of Nuclear Medicine and Molecular Imaging, Department of Radiology, Massachusetts General Hospital, Boston, Massachusetts 02144, United States. 8. Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany.
Abstract
Histone deacetylase 6 (HDAC6) function and dysregulation have been implicated in the etiology of certain cancers and more recently in central nervous system (CNS) disorders including Rett syndrome, Alzheimer's and Parkinson's diseases, and major depressive disorder. HDAC6-selective inhibitors have therapeutic potential, but in the CNS drug space the development of highly brain penetrant HDAC inhibitors has been a persistent challenge. Moreover, no tool exists to directly characterize HDAC6 and its related biology in the living human brain. Here, we report a highly brain penetrant HDAC6 inhibitor, Bavarostat, that exhibits excellent HDAC6 selectivity (>80-fold over all other Zn-containing HDAC paralogues), modulates tubulin acetylation selectively over histone acetylation, and has excellent brain penetrance. We further demonstrate that Bavarostat can be radiolabeled with 18F by deoxyfluorination through in situ formation of a ruthenium π-complex of the corresponding phenol precursor: the only method currently suitable for synthesis of [18F]Bavarostat. Finally, by using [18F]Bavarostat in a series of rodent and nonhuman primate imaging experiments, we demonstrate its utility for mapping HDAC6 in the living brain, which sets the stage for first-in-human neurochemical imaging of this important target.
Histone deacetylase 6 (HDAC6) function and dysregulation have been implicated in the etiology of certain cancers and more recently in central nervous system (CNS) disorders including Rett syndrome, Alzheimer's and Parkinson's diseases, and major depressive disorder. HDAC6-selective inhibitors have therapeutic potential, but in the CNS drug space the development of highly brain penetrant HDAC inhibitors has been a persistent challenge. Moreover, no tool exists to directly characterize HDAC6 and its related biology in the living human brain. Here, we report a highly brain penetrant HDAC6 inhibitor, Bavarostat, that exhibits excellent HDAC6 selectivity (>80-fold over all other Zn-containing HDAC paralogues), modulates tubulin acetylation selectively over histone acetylation, and has excellent brain penetrance. We further demonstrate that Bavarostat can be radiolabeled with 18F by deoxyfluorination through in situ formation of a ruthenium π-complex of the corresponding phenol precursor: the only method currently suitable for synthesis of [18F]Bavarostat. Finally, by using [18F]Bavarostat in a series of rodent and nonhuman primate imaging experiments, we demonstrate its utility for mapping HDAC6 in the living brain, which sets the stage for first-in-human neurochemical imaging of this important target.
Histone deacetylases
have emerged as a pharmaceutical target with
a broad range of potential indications. Several pan-HDAC inhibitors
(not paralogue selective) of different chemotypes, which target several
of the 11 paralogues of Zn-dependent HDACs and span different chemotypes,
are approved by the FDA or are currently in clinical trials.[1] However, these nonselective agents typically
lead to undesired side effects, and considerable efforts are still
spent toward the development of more selective chemical probes and
drug candidates.[2−4]HDAC6 has received recent interest because
its structure is unique
among the paralogues. Its cytosolic location reduces pharmaceutical
interference with genetic material, and HDAC6-selective inhibitors
have shown promise preclinically in terms of efficacy and safety.[5] Paralogue selectivity has been generally difficult
to engineer, but HDAC6 (along with HDAC3 and HDAC8 to a lesser extent[2,6]) has enough structural distinction from other paralogues to offer
a starting point for rational design of selective inhibitors.[7,8] Such paralogue-selective chemical probes are needed to understand
the contribution of HDAC6 to disease: aberrant HDAC6 expression levels
have been heavily implicated in the pathophysiology of glioblastoma
multiforme,[9−11] Rett syndrome,[12,13] and Alzheimer’s[14,15] and Parkinson’s diseases.[16−18] Nonetheless, the design
of brain-penetrant HDAC6-selective agents has proven challenging,
and high doses are needed to achieve any functional effects of HDAC6
inhibition in the CNS.[19]As a consequence,
our understanding of the function and dysfunction
of HDAC6 and the therapeutic potential of its inhibition in the living
human brain remains limited. A positron emission tomography (PET)
probe to study HDAC6 has great potential not only to gain insights
into the molecular underpinnings of brain function and disease but
also in the validation of therapeutic targets and therapeutic small
molecules. Herein, we describe the de novo development
of a brain-penetrant, selective HDAC6 inhibitor and its application
in rodent and nonhuman primate PET imaging, enabled by a recently
described deoxyfluorination method that leverages the reactivity of
a phenol–ruthenium π-complex intermediate.
Results and Discussion
Inhibitor
Design
Several strategies exist to achieve
selectivity for HDAC6 over other HDAC paralogues with hydroxamate
inhibitors. These include the use of (1) large polycyclic capping
groups,[3,20,21] (2) a substituted
sp2 hybridized α-carbon relative to the hydroxamate,[22] and (3) sterically demanding substitution close
to the chelating group.[7,8] On the basis of these strategies,
combined with our experience in PET radiotracer design for HDACs,[23−26] we developed a small library of potential CNS-penetrant HDAC6 inhibitors.
Of these, Bavarostat (Figure A) emerged as unique. In Figure B,C, the structure of Bavarostat is shown
docked into the catalytic domain 2 (CD2) of humanHDAC6, providing
a model to the contribution of structural elements to its selectivity
for HDAC6 over all other HDAC paralogues (Figure D) and exemplified by comparison to HDAC2.
Figure 1
A: Structure
of Bavarostat (1). B: Alignment of hHDAC6’s
catalytic domain 2 (CD2) (white and cyan protein, PDB 5EDU) and the 15 zHDAC6
CD2 crystal structures available in the PDB database (wheat colored
proteins), showing similar L1 loop and L1 helix conformations. HDAC6-selective
inhibitor Bavarostat was docked into the hHDAC6 complex described
in the Experimental Section. Bavarostat is
predicted to form a hydrogen bond (yellow dashed line) with Ser568
on the hHDAC6 L1 loop segment represented in cyan. C: Alignment of
hHDAC6 CD2 (white and cyan) and the five hHDAC2 catalytic site crystal
structures (pink proteins) in PDB. Large movement of both the L1 loop
segment and the L1 helix is represented by arrows (arrows depict movement
of analogous amino acids from hHDAC6 to hHDAC2: left arrow depicts
hHDAC6’s Ile571 going to corresponding hHDAC2’s Asp100′,
and right arrow depicts hHDAC6’s Arg557 going to corresponding
hHDAC2’s Ser86′). Prime indicates hHDAC2. The hypothesized
steric and electrostatic clash between Bavarostat and the hHDAC2 L1
loop is represented in red. D: IC50 values for Bavarostat
inhibition of HDAC1–11 were determined in a fluorescence assay
(BPS Bioscience, San Diego, CA) measuring acetylation of a synthetic
substrate. E: Bavarostat docked into a hHDAC6 CD2 complex incorporating
a conserved water. Key amino acids chelating the catalytic zinc or
lining the 10 Å channel are labeled. F/G: Significant topography
differences between the surface of the catalytic domains of hHDAC6
(F) and hHDAC2 (G) are highlighted by the pink contouring of the L1
loop segment of hHDAC2.
A: Structure
of Bavarostat (1). B: Alignment of hHDAC6’s
catalytic domain 2 (CD2) (white and cyan protein, PDB 5EDU) and the 15 zHDAC6
CD2 crystal structures available in the PDB database (wheat colored
proteins), showing similar L1 loop and L1 helix conformations. HDAC6-selective
inhibitor Bavarostat was docked into the hHDAC6 complex described
in the Experimental Section. Bavarostat is
predicted to form a hydrogen bond (yellow dashed line) with Ser568
on the hHDAC6 L1 loop segment represented in cyan. C: Alignment of
hHDAC6 CD2 (white and cyan) and the five hHDAC2 catalytic site crystal
structures (pink proteins) in PDB. Large movement of both the L1 loop
segment and the L1 helix is represented by arrows (arrows depict movement
of analogous amino acids from hHDAC6 to hHDAC2: left arrow depicts
hHDAC6’s Ile571 going to corresponding hHDAC2’s Asp100′,
and right arrow depicts hHDAC6’s Arg557 going to corresponding
hHDAC2’s Ser86′). Prime indicates hHDAC2. The hypothesized
steric and electrostatic clash between Bavarostat and the hHDAC2 L1
loop is represented in red. D: IC50 values for Bavarostat
inhibition of HDAC1–11 were determined in a fluorescence assay
(BPS Bioscience, San Diego, CA) measuring acetylation of a synthetic
substrate. E: Bavarostat docked into a hHDAC6 CD2 complex incorporating
a conserved water. Key amino acids chelating the catalytic zinc or
lining the 10 Å channel are labeled. F/G: Significant topography
differences between the surface of the catalytic domains of hHDAC6
(F) and hHDAC2 (G) are highlighted by the pink contouring of the L1
loop segment of hHDAC2.Bavarostat was docked into an hHDAC6 complex to rationalize
its
selectivity (see Experimental Section for
further details). Briefly, the hydroxamic acid of Bavarostat was modeled
to chelate the catalytic zinc and form a hydrogen bond to an ordered
water in a similar fashion as seen in the zHDAC6 structure (PDB 5EF7) recently published
by Hai et al.[27] Indeed, the rigidity introduced
by the linker phenyl (sp2-hybridized hydroxamate α-carbon
atom) in Bavarostat prevents an effective κ2-binding
mode of the hydroxamate warhead. Since a high-resolution crystal structure
of zHDAC6 bound to a phenyl-linked HDAC6-selective inhibitor, HPOB
in PDB 5EF7,
showed a unique mode of binding to the catalytic zinc ion via a conserved
water molecule,[27] we applied it to Bavarostat.
The presence of this water molecule greatly improved the fit of Bavarostat
to the protein, suggesting that the rigid phenyl linker moiety is
partially conveying HDAC6 selectivity. Additional discriminating factors,
such as sterically demanding moieties, were employed to convey paralogue
selectivity and potency for HDAC6. First, the fluorophenyl of Bavarostat
occupying the 10 Å channel leading to the protein surface forms
a favorable noncovalent interaction with Phe620 via a parallel-displaced π-stacking
configuration (Figure E). The fluorine substituent, installed with imaging applications
in mind in a position not expected to interfere with HDAC binding,[26] contributes to favorable van der Waals interactions
by vectoring in a small divot in the 10 Å channel (Figure SI3).Second, the proton of the
tertiary amine off of the linker phenyl
ring directly hydrogen bonds to Ser568 on the L1 loop (Figures B and 1E). Finally, the adamantyl group in Bavarostat is predicted to fill
a hydrophobic surface groove adjacent to the L1 loop in hHDAC6 complex’s CD2 (Figure F). To understand if the conformation of the L1 loop
and L1 helix found in the hHDAC6 CD2 structure (PDB 5EDU) was an anomaly
due to the crystallization technique used or the result of the bound
ligand that was crystallized with it, all zHDAC6 CD2 structures available
were aligned to hHDAC6 CD2[27] (Figure B). The L1 loop and
L1 helix in CD2 of hHDAC6 and CD2 of zHDAC6 structures consistently
folded in the same conformation, suggesting that the L1 loop and L1
helix of hHDAC6 are in a stable low energy conformation. The selectivity
profile of Bavarostat (Figure D) shows that it is highly selective for hHDAC6 over hHDAC2.
To rationalize how Bavarostat may have gained such selectivity, the
catalytic site of hHDAC2 was used as a comparison and a representative
for class I HDACs. We aligned the CD2 structure of hHDAC6 (PDB 5EDU) to all available
CD1 hHDAC2 structures to see if the L1 loop and L1 helix were in a
similar conformation since the adamantyl of Bavarostat sits on a hydrophobic
groove next to the L1 loop. The L1 loop and L1 helix of hHDAC2 were
found to be in a consistent conformation across all available structures,
suggesting that this folded state is also stable. This conformation,
however, is significantly different relative to CD2 hHDAC6 as exemplified
by the large movement of both the L1 loop segment and the L1 helix depicted by arrows in Figures B and 1C). The L1 loop/L1 helix conformation of CD1 of hHDAC2 would
sterically clash with the bound Bavarostat compound (Figure C), preventing its positioning
into the analogous hydrophobic groove (Figure G).Our modeling suggests that bulky
hydrophobic adamantyl group is
best suited to mediate nonpolar surface interactions while affording
HDAC6 selectivity and favorable physicochemical properties increasing
Bavarostat’s likelihood to penetrate the blood–brain
barrier (BBB). Further, the high protein density of HDAC6 in human
brain (5–10 pmol per mg tissue)[28] makes the IC50 of 60 nM promising for imaging applications.
Functional Selectivity for HDAC6 in Human iPSC-Derived Neuronal
Cells
HDAC6 is a class IIb histone deacetylase, structurally
unique among other HDACs, predominantly located in the cytosol where
it displays distinct substrate selectivity.[29−31] Most prominently
α-tubulin and HSP90 are deacetylated by HDAC6, while nuclear
proteins such as histones 3 and 4 (H3, H4) remain unaltered.[5,31] Therefore, the inhibition of HDAC6 should increase the acetylation
of α-tubulin (α-tubulin-ac) but not histone acetylation
such as the acetylation of histone H3 lysine 9 (H3K9ac) and histone
H4 lysine 12 (H4K12ac). To assess the selective inhibition of HDAC6
by Bavarostat in intact cells, we evaluated the cellular acetylation
levels by Western blot analysis of a known substrate (α-tubulin-ac)
along with two nonsubstrates of HDAC6 (H3K9ac and H4K12ac) in human
induced pluripotent stem cell (iPSC) derived neural progenitor cells
(NPCs).[28] We treated NPCs with Bavarostat
(10 μM) as well as several HDACi chemical probes including CI-994[32] (a class I selective HDACi, 10 μM), ACY1215[5,33] (a mixed selectivity inhibitor, 5 μM), and Tubastatin A[7] (a selective HDAC6 inhibitor, 10 μM) for
6 h (Figure ).
Figure 2
In
vitro functional selectivity of Bavarostat
for HDAC6. Histone 3 (lysine 9), histone 4 (lysine 12), and α-tubulin
protein acetylation levels in IPS neuroprogenitor cells after treatment
with Bavarostat (10 μM), a class I HDAC inhibitor (CI-994, 10
μM), an HDAC6-selective inhibitor (Tubastatin A, 10 μM),
and a mixed selectivity inhibitor (ACY1215, 5 μM), in comparison
to DMSO-treated cells, indicating functional HDAC6 selectivity for
Bavarostat.
In
vitro functional selectivity of Bavarostat
for HDAC6. Histone 3 (lysine 9), histone 4 (lysine 12), and α-tubulin
protein acetylation levels in IPS neuroprogenitor cells after treatment
with Bavarostat (10 μM), a class I HDAC inhibitor (CI-994, 10
μM), an HDAC6-selective inhibitor (Tubastatin A, 10 μM),
and a mixed selectivity inhibitor (ACY1215, 5 μM), in comparison
to DMSO-treated cells, indicating functional HDAC6 selectivity for
Bavarostat.As expected, CI-994,
the class I HDACi, showed a significant increase
in H3K9ac and H4K12ac but not α-tubulin-ac. ACY-1215 also increased
H3K9ac and H4K12ac, and additionally, it effectively increased α-tubulin-ac,
consistent with its mixed selective nature for class I HDACs and HDAC6.
Tubastatin A, the highly HDAC6 selective inhibitor, significantly
increased α-tubulin-ac but not H3K9-ac or H4K12-ac. Bavarostat,
the novel HDAC6-selective inhibitor, exhibited specific increased
α-tubulin-ac only. Our data indicate that Bavarostat is able
to access the cytoplasmic compartment of human neuronal cells, engage
with endogenous HDAC6, and inhibit its deacetylase activity. In conclusion,
the changes in cellular protein acetylation in response to HDACi treatment
confirm the functional selectivity of Bavarostat for HDAC6 versus
class I HDACs.
Radiochemical Synthesis
Radiolabeling
of hydroxamic
acids, particularly scaffolds like Bavarostat, has been difficult
to achieve. Routes were lengthy, inefficient, low-yielding, and notoriously
difficult to automate.[26] Enabled by innovative
ruthenium-mediated radiofluorination[34] (Figure A), a highly efficient
radiosynthesis afforded [18F]Bavarostat in high yield and
high specific radioactivity. Using air-stable ruthenium complex 3, the labeling precursor 2 was converted into
an η6-coordinated ruthenium–phenol complex
solution and served without further purification to elute [18F]fluoride off an anion exchange cartridge with 83% elution efficiency.
Within 30 min at 130 °C, the labeling proceeds with high conversion
(>70% radiochemical yield estimated by TLC). Subsequent transacylation
in the same pot afforded [18F]Bavarostat in 8.1% (n = 2) overall nondecay corrected radiochemical yield and
high specific activity (4 mCi·nmol–1 at end
of synthesis) after HPLC purification and reformulation.
Figure 3
Rodent PET
imaging with [18F]Bavarostat. A: Radiochemical
synthesis of [18F]Bavarostat via ruthenium-mediated deoxyfluorination,
ClIM = N,N-bis(2,6-diisopropylphenyl)-1-chloroimidazolium
chloride. B: Representative sagittal slices summed from 30 to 90 min.
C: Averaged (n = 3) time–activity curves of
a whole-brain ROI of Sprague–Dawley rats injected with [18F]Bavarostat. In the blocked animals, 1.0 mg·kg–1 unlabeled Bavarostat was injected immediately prior
to radiotracer administration, baseline animals treated with vehicle.
Rodent PET
imaging with [18F]Bavarostat. A: Radiochemical
synthesis of [18F]Bavarostat via ruthenium-mediated deoxyfluorination,
ClIM = N,N-bis(2,6-diisopropylphenyl)-1-chloroimidazolium
chloride. B: Representative sagittal slices summed from 30 to 90 min.
C: Averaged (n = 3) time–activity curves of
a whole-brain ROI of Sprague–Dawley rats injected with [18F]Bavarostat. In the blocked animals, 1.0 mg·kg–1 unlabeled Bavarostat was injected immediately prior
to radiotracer administration, baseline animals treated with vehicle.
Rodent Imaging
To evaluate the potential of [18F]Bavarostat to serve
as an HDAC6 radiotracer, PET imaging in rodents
was performed. Sprague–Dawley rats were either pretreated with
vehicle (1 mL·kg–1, saline:Tween 80:DMSO 8:1:1)
or unlabeled Bavarostat (1 mg·kg–1, 1 mg·mL–1) immediately before radiotracer injection, and the
dynamic PET data was acquired. As shown in the time–activity
curves in Figure C,
the compound exhibits excellent brain uptake and retention. Treatment
of the animals with nonradioactive Bavarostat at 1 mg·kg–1 led to solid blocking of brain uptake, indicative
of specific binding (Figure B). Normalized time–activity curves were highly reproducible,
and signal magnitude was therefore set to 100% at t = 2 min to compare the outcome of several replicates. More work
to determine a suitable input function will be needed to appropriately
model the tracer kinetics for more in-depth analysis.
Autoradiography
Heterologous in vivo blocking studies were inconclusive
due to the poor brain penetrance
of known HDAC6-selective inhibitors. To circumvent the issue of blood–brain
barrier penetrance, selectivity of the probe in brain tissue was demonstrated
using autoradiography (Figure ). Sagittal slices of rat brain tissue were exposed to [18F]Bavarostat in the presence of [19F]Bavarostat
or Tubastatin A. At 10 μM [19F]Bavarostat, binding
was reduced by almost 40% in the cerebellum, relative to corpus callosum,
a region of high nonspecific binding. Tubastatin A, an HDAC6-selective
compound, reduced the amount of bound radioactivity to a comparable
level as Bavarostat. These data indicate that the binding to brain
tissue occurs with selectivity for HDAC6.
Figure 4
Autoradiography with
[18F]Bavarostat. A: Signal intensity
of cerebellum relative to corpus callosum at baseline and in sagittal
slices of Sprague–Dawley rat brain pretreated with 10 μM
Bavarostat or Tubastatin A respectively. Blocking index = 1 –
[(PTCB·PTCC–1)/(BLCB·BLCC–1)],
PT = pretreated, BL = vehicle treatment, CB = cerebellum, CC = corpus
callosum. ****P < 0.0001. B: Representative sagittal
slices of Sprague–Dawley rat brains pretreated with vehicle,
10 μM Bavarostat, or Tubastatin A respectively.
Autoradiography with
[18F]Bavarostat. A: Signal intensity
of cerebellum relative to corpus callosum at baseline and in sagittal
slices of Sprague–Dawley rat brain pretreated with 10 μM
Bavarostat or Tubastatin A respectively. Blocking index = 1 –
[(PTCB·PTCC–1)/(BLCB·BLCC–1)],
PT = pretreated, BL = vehicle treatment, CB = cerebellum, CC = corpus
callosum. ****P < 0.0001. B: Representative sagittal
slices of Sprague–Dawley rat brains pretreated with vehicle,
10 μM Bavarostat, or Tubastatin A respectively.
Nonhuman Primate (NHP) Imaging
To
assess the translational
potential of [18F]Bavarostat for human imaging, the tracer
properties were assessed in baboon (Papio anubis)
through a series of preliminary experiments (Figure ). The molecule exhibits excellent brain
uptake, with a whole-brain standard uptake value of 3 around 30 min
after tracer injection and high retention of signal over time. The
uptake of [18F]Bavarostat is particularly high in subcortical
areas and low in white matter, distinct from other HDAC imaging agents
assessed in NHPs[23,25,26,35−38] (see the Supporting Information). Particularly noteworthy is the low
amount of nonspecific binding observed after pretreatment with 1 mg·kg–1 nonlabeled Bavarostat, and the localization of the
residual signal mostly to white matter, which is also observed in
autoradiography (Figure ). An investigation of in vivo metabolism and kinetic
modeling of [18F]Bavarostat binding is currently ongoing.
In a preliminary analysis, a good correlation was obtained between
standardized uptake value (SUV) and distribution volume (VT) derived by Logan plot from an arterial plasma input
function (Supporting Information). In summary,
the NHP imaging data are promising and warrant translation of [18F]Bavarostat for humanHDAC6-imaging.
Figure 5
Nonhuman primate imaging
with [18F]Bavarostat. A: Sagittal
and coronal slices of PET images acquired with [18F]Bavarostat
in the baboon brain, baseline and pretreated with 1.0 mg·kg–1 unlabeled Bavarostat (averaged 60–120 min).
B: Whole-brain SUV time–activity curves of [18F]Bavarostat
in baboon brain, baseline and pretreated with 1.0 mg·kg–1 unlabeled Bavarostat.
Nonhuman primate imaging
with [18F]Bavarostat. A: Sagittal
and coronal slices of PET images acquired with [18F]Bavarostat
in the baboon brain, baseline and pretreated with 1.0 mg·kg–1 unlabeled Bavarostat (averaged 60–120 min).
B: Whole-brain SUV time–activity curves of [18F]Bavarostat
in baboon brain, baseline and pretreated with 1.0 mg·kg–1 unlabeled Bavarostat.Conceptually, the animal imaging experiments reveal two important
pieces of information: (i) that the labeled small molecule [18F]Bavarostat at a microdose can easily access the brain compartment
in living animals with high specific binding and (ii) that nonradioactive
Bavarostat at a putative pharmacological dose of 1 mg·kg–1 is well tolerated by two species at least by acute
iv treatment. The development of Bavarostat, a brain-penetrant and
highly selective HDAC6 inhibitor, offers a powerful chemical probe
for preclinical testing to reveal the impact of HDAC6 activity on
brain chemistry and behavior.
Conclusion
We
show that Bavarostat selectively inhibits
HDAC6 in recombinant and cellular assays. Radiolabeling of [18F]Bavarostat was possible with innovative radiofluorination, and in vivo imaging in both rodent and nonhuman primate demonstrated
the high brain uptake of the molecule. Its engagement of HDAC6 was
confirmed by autoradiography. Bavarostat shows great promise as a
radiotracer in rodent and nonhuman primate models and is expected
to enable the study of HDAC6 in the living human brain.
Experimental
Section
Chemical Synthesis
The synthesis of Bavarostat, the
labeling precursor 2, and ruthenium complex 3 is described in the Supporting Information.
IC50
IC50 measurements were conducted
by BPS Biosciences with an established fluorescence assay, using proprietary
class-specific substrates and recombinant enzymes.
Molecular Modeling
Protein
Preparation
The three HDAC structures studied
in the modeling section of this work, humanHDAC2 (hHDAC2, PDB 3MAX),[39] zebrafishHDAC6 (zHDAC6, PDB 5EF7),[27] and humanHDAC6 (hHDAC6, PDB 5EDU),[27] were aligned with the Protein Structure
Alignment function of Maestro, Schrödinger Release 2017-1 (Maestro,
Schrödinger, LLC, New York, NY, 2017) using the default settings.The HDAC6 system used for the docking studies described in this
work (referred to as the hHDAC6 complex from here on) was created
from two structures, zHDAC6 (PDB 5EF7) and hHDAC6 (PDB 5EDU), to exploit key
features from both. The zebrafishHDAC6 (zHDAC6, PDB 5EF7) system was selected
since it is a holo structure bound to small molecule HPOB, which has
the same warhead and linker moiety, a phenyl hydroxamic acid, as Bavarostat.
The hydroxamic acid of HPOB in the zHDAC6 PDB 5EF7 structure makes
two specific nonbonded contacts in the active site, one chelating
the zinc and the other forming a hydrogen bond to an ordered water.
In this regard, from the zHDAC6 structure (PDB 5EF7) the zinc, ligand
(HPOB), and ordered water adjacent to the zinc were used to build
the hHDAC6 complex. The protein and sequence from the unique hHDAC6
(PDB 5EDU) was
also used to create the hHDAC6 complex. This new hHDAC6 complex system
was then run through the Protein Preparation Wizard[40] in Maestro, Schrödinger Release 2017-1 (Protein
Preparation Wizard, Schrödinger, LLC, New York, NY, 2017),
to protonate it and add any missing side chain residues using the
OPLS3 force field.[41,42] This preparation was then followed
by an all atom minimization using the default settings to create a
final chimeric hHDAC6 complex system.
Receptor Preparation for
Docking
The final minimized hHDAC6 complex described above was prepared for docking using the
Receptor Grid Generation program in Maestro. The HPOB ligand was removed
from the hHDAC6 complex for docking. To enhance binding mode prediction
a select number of binding pocket hydroxyls and thiols in the side
chains were allowed to rotate during the docking simulation: these
residues were Ser568, Cys618, and Tyr782. To address modeling into
the tight receptor of the HDAC6 binding pocket channel, a GRID scaling
factor of 0.6 with a partial charge cutoff of 0.25 was selected.
Ligand Preparation and Docking
The Bavarostat ligand
studied in this work was prepared for docking using LigPrep in Maestro,
Schrödinger Release 2017-1 (LigPrep, Schrödinger, LLC,
New York, NY, 2017) and the OPLS3 force field.[41,42] To enable ligand docking into the tight channel of the hHDAC6 complex
a van der Waals radius of 0.7 and a partial charge cutoff of 0.15
were selected for the ligand. A targeted pH of 7.0 ±1.0 was selected.
Glide, Schrödinger Release 2017-1 (Glide, Schrödinger,
LLC, New York, NY, 2017), was used for docking of Bavarostat into
the hHDAC6 complex described above.[43] To
gain accuracy in the docking simulation into the hHDAC6 complex system,
the 3D coordinates of the four heavy atoms of the hydroxamic acid
of Bavarostat were constrained to their equivalent atoms in HPOB in
zHDAC6 (PDB 5EF7). The top scoring docking pose was selected for analysis. All modeling
figures in this work were made with the PyMOL Molecular Graphics System,
Version 1.8 Schrödinger, LLC.
Cellular Protein Acetylation
Levels
Human iPSC derived
neural progenitor cells from a healthy control subject fibroblast
cell line GM08330 (Coriell Institute for Medical Research) were generated
as described previously[28] and treated with
DMSO or a solution of HDAC inhibitor (ACY1215, final concentration
5 μM; Bavarostat, CI-994, Tubastatin A, final concentration
10 μM) for 6 h at 37 °C. Lysis of cell pellets (n = 3 per condition) was performed in radioimmunoprecipitation
assay (RIPA) buffer (Boston BioProducts #BP-115) with EDTA-free protease
inhibitors (Sigma #4693159001). The lysates were centrifuged at 18,000
rpm at 4 °C for 15 min, and the supernatants were collected.
The protein concentration was determined by BCA assay (Thermo Scientific
#23227). Western blot analysis was conducted on samples adjusted to
6 μg of total protein/replicate.
Radiosynthesis of [18F]Bavarostat
Aqueous
[18F]fluoride obtained from a cyclotron was passed through
a SPE Chromafix 30-PS-HCO3 cartridge that had been previously
conditioned with 5.0 mg·mL–1 aqueous potassium
carbonate and then washed with 18 mL of Millipore Milli-Q water. The
captured [18F]fluoride was washed by passing 1 mL of ethanol
through the cartridge. At the beginning of the synthesis, 657 mCi
was measured on the cartridge.5 mg of methyl 4-((((adamantan-1-yl)methyl)(methyl)amino)methyl)-3-hydroxybenzoate,
10 mg of CpRu(cod)Cl, and 30 mg of N,N-bis(2,6-diisopropyl)phenyl-2-chloroimidazolium chloride were heated
in 250 μL of ethanol at 85 °C for 30 min. The resulting
solution was passed through the anion exchange cartridge and collected
into a dram vial. The cartridge was flushed with 400 μL of acetonitrile
and 400 μL of DMSO and collected into the same vial, which was
subsequently sealed with a Teflon lined cap and heated to 130 °C
for 30 min. Then, 1 mL of THF/MeOH (1:1), 0.4 mL of 50% aqueous NH2OH, and 0.1 mL of 5 M NaOH were added at room temperature
and the reaction was stirred for 5 min.The solution was diluted
with water to 10 mL, loaded onto an OASIS
MAX SPE cartridge (60 mg), washed with 5 mL of water, eluted with
2 mL of ethanol/0.1 M AcOH (1:1), and purified by semipreparative
HPLC (Agilent Eclipse C-18, 9.4 × 250 mm, 5 μm; flow ramp
0.5 mL·min–1 to 5 mL·min–1 from 0 to 4 min, then 5 mL·min–1; 5% ACN
in 0.01 N NaOH from 0 to 4 min, then ramp to 70% ACN in 0.01 NaOH
at 45 min). The isolated fraction was reformulated on an OASIS MAX
SPE cartridge (60 mg), washed with 5 mL of water, eluted with 2 mL
of ethanol/0.1 M AcOH (1:1), diluted with 8 mL of 0.9% saline, and
neutralized with 0.1 N NaOH to pH 5. Overall, 53.3 mCi was isolated
(8.1% activity yield, 4 mCi·nmol–1 specific
activity) within 94 min. For more analytical data, refer to the Supporting Information.Sagittal slices of Sprague–Dawley
rat brains were sectioned (10 μm) with a cryostat (Thermo Scientific
HM550), mounted onto ColorFrost Plus microscope slides (Fisher Scientific
12-550-18), and stored at −20 °C. The slides were submerged
in baths of 50 mL of buffer (100 mM Tris, 50 mM NaCl, pH 7.5, 1% BSA)
at room temperature and equilibrated for 15 min. The slides were then
transferred into baths (50 mL of buffer) containing either 50 μL
of DMSO or a solution of blocking compound in DMSO (Bavarostat at
final concentrations of 100 nM, 1 μM, and 10 μM and Tubastatin
A at a final concentration of 100 μM). After 15 min, 50 μCi
of [18F]Bavarostat was added to each bath. After 15 min, all slides were washed by dipping 10× into buffer
and subsequently submerged in 50 mL baths of buffer at 4 °C for
3 min. The slides were washed by dipping 10× into wash buffer
(100 mM Tris, 50 mM NaCl, pH 7.5) and were dried at 35 °C in
a vacuum chamber. A phosphorus screen (PerkinElmer 7001723) was exposed
with the slides for 45 min and subsequently imaged with a Cyclone
Plus Storage Phosphor (PerkinElmer) detector. ImageJ was used to apply
a Gaussian blur (3.0 radius) smoothing, and a lookup table (Royal)
with equivalent thresholds for brightness was applied. Raw intensity
values from gray and white matter were quantified with the ImageJ
measurement tool.
Animal Preparation
Rodent
Six male
Sprague–Dawley rats (Charles
River Laboratories) were used for PET imaging. According to institutional
policy, animals were pair-housed until a weight of 500 g, and a 12
h/12 h light/dark cycle was applied. All treatment and imaging experiments
were performed in accordance with procedures approved by the Institutional
Animal Care and Use Committee at the Massachusetts General Hospital.
Anesthesia was achieved with isoflurane in medical oxygen carrier
(3% for induction, 2% for maintenance). For intravenous administration,
a catheter with an extension line was placed in a lateral tail vein.
Each animal received a bolus injection of either vehicle (1 mL·kg–1, 1:1:8 DMSO/Tween80/saline) or blocking agent in
solution (1 mg·kg–1, 1 mg·mL–1 in 1:1:8 DMSO/Tween80/saline) immediately prior to injection of
the radiotracer.
NHP
Two PET/MR studies (baseline/pretreatment)
were
carried out with one baboon (female, 16 kg, Papio Anubis) as approved by the Institutional Animal Care and Use Committee
at the Massachusetts General Hospital. Nil per os was instructed 12 h prior to the study. Anesthesia was induced with
intramuscular (im) ketamine (1.5 mL, 100 mg·mL-1) + xylazine (0.35 mL, 20 mg·mL-1), and refreshed 45 min later with
ketamine (0.2 mL, 100 mg/mL) + xylazine (0.1 mL, 20 mg/mL). Anesthesia
was continued during the study with 1–1.5% isoflurane in oxygen
carrier. Pretreatment (1 mg·kg–1 Bavarostat,
5 mg·mL–1 in 10% DMSO, 10% Tween 80, 80% saline)
was administered in one of the scans 5 min prior to radiotracer injection
through a catheter in the saphenous vein. Radiotracer injection was
performed through the same catheter. Vital signs (end-tidal CO2, oxygen saturation, heart rate, and respiration rate) were
under continuous surveillance to maintain a normal physiological range,
documented every 15 min.
PET/CT Image Acquisition
After injection of a radiotracer bolus (ca.
700 μCi of [18F]Bavarostat), a 90 min dynamic PET
scan was acquired in pairs for all animals. PET scans were performed
on a GammaMedica Triumph PET/CT/SPECT scanner, corrected for attenuation
with a μ-map derived from the corresponding CT image, which
was acquired immediately following the PET scan. The dynamic PET data
was binned into 38 time frames (8 × 15 s, 8 × 1 min, 10
× 2 min, 12 ×
5 min) and reconstructed individually via an iterative MLEM (maximum
likelihood expectation maximization) algorithm in 16 iterations.PET/MRI acquisition was performed on a 3T Siemens
TIM-Trio with a BrainPET insert (Siemens, Erlangen, Germany). A PET/MRI
compatible eight-channel array coil customized for nonhuman primate
brain imaging to increase image signal and quality was employed. After
administration of the radiotracer (5.2 mCi baseline, 2.9 mCi block),
dynamic PET image acquisition was initiated. Dynamic PET data were
collected and stored in list mode for 120 min. Image reconstruction
was performed using the 3D ordinary Poisson expectation maximization
algorithm with detector efficiency, decay, dead time, attenuation,
and scatter corrections. PET data were binned in 29 frames (6 ×
10 s, 6 × 20 s, 2 × 30 s, 1 × 1 min, 5 × 5 min,
9 × 10 min). Image volumes were eventually reconstructed into
76 slices with 128 × 128 pixels and a 2.5 mm isotropic voxel
size. Thirty minutes after scanner start, a high-resolution anatomical
scan using multiecho MPRAGE sequence (TR = 2530 ms, TE1/TE2/TE3/TE4
= 1.64/3.5/5.36/7.22 ms, TI = 1200 ms, flip angle = 7°, and 1
mm isotropic) was acquired.
Image Analysis
PET images were coregistered to the CT image
acquired from the same animal using AMIDE. Data sets were cropped,
and all further image analysis was conducted using PMOD 3.3 (PMOD
Technologies, Zurich, Switzerland). For maximum consistency, the data
were manually coregistered to the Schiffer Px Rat rat brain template[44] and data was derived from a whole-brain VOI
(volume of interest) for time–activity curves.Using the Martinos jip analysis toolkit (www.nitrc.org/projects/jip), MPRAGE data was coregistered to the black baboon atlas[45] and the transformation matrix applied to the
dynamic PET data. Using PMOD 3.3, time–activity curves for
VOIs were extracted, and time-averaged SUV values were obtained. Kinetic
modeling was performed with a metabolite corrected arterial plasma
input function in PMOD.
Authors: Richard A Friesner; Robert B Murphy; Matthew P Repasky; Leah L Frye; Jeremy R Greenwood; Thomas A Halgren; Paul C Sanschagrin; Daniel T Mainz Journal: J Med Chem Date: 2006-10-19 Impact factor: 7.446
Authors: Hsin-Hsien Yeh; Mei Tian; Rainer Hinz; Daniel Young; Alexander Shavrin; Uday Mukhapadhyay; Leo G Flores; Julius Balatoni; Suren Soghomonyan; Hwan J Jeong; Ashutosh Pal; Rajesh Uthamanthil; James N Jackson; Ryuichi Nishii; Hiroshi Mizuma; Hirotaka Onoe; Shinya Kagawa; Tatsuya Higashi; Nobuyoshi Fukumitsu; Mian Alauddin; William Tong; Karl Herholz; Juri G Gelovani Journal: Neuroimage Date: 2012-09-17 Impact factor: 6.556
Authors: Mar Cuadrado-Tejedor; Carolina Garcia-Barroso; Juan A Sánchez-Arias; Obdulia Rabal; Marta Pérez-González; Sara Mederos; Ana Ugarte; Rafael Franco; Victor Segura; Gertrudis Perea; Julen Oyarzabal; Ana Garcia-Osta Journal: Neuropsychopharmacology Date: 2016-08-23 Impact factor: 7.853
Authors: Florence F Wagner; David E Olson; Jennifer P Gale; Taner Kaya; Michel Weïwer; Nadia Aidoud; Méryl Thomas; Emeline L Davoine; Bérénice C Lemercier; Yan-Ling Zhang; Edward B Holson Journal: J Med Chem Date: 2013-02-18 Impact factor: 7.446
Authors: Alan P Kozikowski; Sida Shen; Marta Pardo; Maurício T Tavares; Dora Szarics; Veronick Benoy; Chad A Zimprich; Zsófia Kutil; Guiping Zhang; Cyril Bařinka; Matthew B Robers; Ludo Van Den Bosch; James H Eubanks; Richard S Jope Journal: ACS Chem Neurosci Date: 2018-12-14 Impact factor: 4.418
Authors: Nicholas J Porter; Jeremy D Osko; Daniela Diedrich; Thomas Kurz; Jacob M Hooker; Finn K Hansen; David W Christianson Journal: J Med Chem Date: 2018-08-17 Impact factor: 7.446
Authors: Paris R Watson; Ping Bai; Changning Wang; Abigail D Cragin; Jacob M Hooker; David W Christianson Journal: Biochemistry Date: 2022-09-08 Impact factor: 3.321
Authors: Gengyang Yuan; Maeva Dhaynaut; Yu Lan; Nicolas J Guehl; Dalena Huynh; Suhasini M Iyengar; Sepideh Afshar; Manish Kumar Jain; Julie E Pickett; Hye Jin Kang; Hao Wang; Sung-Hyun Moon; Mary Jo Ondrechen; Changning Wang; Timothy M Shoup; Georges El Fakhri; Marc D Normandin; Anna-Liisa Brownell Journal: J Med Chem Date: 2022-01-28 Impact factor: 8.039
Authors: Sofie Celen; Johanna Rokka; Tonya M Gilbert; Michel Koole; Isabeau Vermeulen; Kim Serdons; Frederick A Schroeder; Florence F Wagner; Tom Bleeser; Baileigh G Hightower; Jiyun Hu; Dania Rahal; M Hassan Beyzavi; Wim Vanduffel; Koen Van Laere; Janice E Kranz; Jacob M Hooker; Guy Bormans; Christopher J Cawthorne Journal: ACS Chem Neurosci Date: 2020-03-19 Impact factor: 4.418
Authors: Jay S Wright; Liam S Sharninghausen; Sean Preshlock; Allen F Brooks; Melanie S Sanford; Peter J H Scott Journal: J Am Chem Soc Date: 2021-04-29 Impact factor: 15.419