Natalia Estrada-Ortiz1, Constantinos G Neochoritis1, Aleksandra Twarda-Clapa2,3, Bogdan Musielak4, Tad A Holak3,4, Alexander Dömling1. 1. Department of Drug Design, University of Groningen, A. Deusinglaan 1, Groningen 9700AV, The Netherlands. 2. Faculty of Biochemistry, Biophysics and Biotechnology, Jagiellonian University, Gronostajowa 7, 30-387 Krakow, Poland. 3. Malopolska Centre of Biotechnology, Jagiellonian University, Gronostajowa 7a, 30-387 Krakow, Poland. 4. Department of Chemistry, Jagiellonian University, Ingardena 3, 30-060 Krakow, Poland.
Abstract
Based on a combination of an Ugi four component reaction and a ring closing metathesis, a library of novel artificial macrocyclic inhibitors of the p53-MDM2 interaction was designed and synthesized. These macrocycles, alternatively to stapled peptides, target for the first time the large hydrophobic surface area formed by Tyr67, Gln72, His73, Val93, and Lys94 yielding derivatives with affinity to MDM2 in the nanomolar range. Their binding affinity with MDM2 was evaluated using fluorescence polarization (FP) assay and 1H-15N two-dimensional HSQC nuclear magnetic resonance experiments.
Based on a combination of an Ugi four component reaction and a ring closing metathesis, a library of novel artificial macrocyclic inhibitors of the p53-MDM2 interaction was designed and synthesized. These macrocycles, alternatively to stapled peptides, target for the first time the large hydrophobic surface area formed by Tyr67, Gln72, His73, Val93, and Lys94 yielding derivatives with affinity to MDM2 in the nanomolar range. Their binding affinity with MDM2 was evaluated using fluorescence polarization (FP) assay and 1H-15N two-dimensional HSQC nuclear magnetic resonance experiments.
Entities:
Keywords:
MDM2; Macrocycles; Ugi reaction; cancer; p53; protein−protein interaction
In nature,
macrocycles are not uncommon and often exhibit interesting biological
activities.[1] Often macrocyclic compounds
have distinct advantages over their open chain analogues including
higher affinity and selectivity,[2] preferable
entropic signature, better membrane permeation, and oral bioavailability
or higher stability.[1,2] Thus, macrocycles can increase
affinity and selectivity for a specific target.[1] In spite of their potential, macrocycles pose considerable
synthesis problems, and also, the accurate prediction of their conformation
makes it difficult to predict activity.[1,2]The tumor
suppressor protein p53 is involved in controlling pathways of cell
cycle, apoptosis, angiogenesis metabolism, senescence, and autophagia.[3,4] The TP53 gene is one of the most frequently mutated genes in a multitude
of humancancer.[5] Additionally, in multiple
cases where TP53 is intact, p53’s function is impaired by its
negative regulators, mouse double minute 2 and 4 homologues (MDM2
and MDMX), due to amplification or enhanced expression of their coding
genes.[3−5] Multiple potent and selective compound classes to
inhibit the p53–MDM2 interaction have been discovered, described,
and evaluated in early clinical trials.[6] However, the pharmacokinetic and pharmacodynamics properties of
the studied scaffolds could still be optimized to minimize the side
effects. Thus, the discovery of new p53–MDM2 inhibitors with
diverse structures to improve their properties is still of importance.[6−8] The three finger pharmacophore model for p53–MDM2 is recognized
as responsible for the binding of small molecules and peptides to
the MDM2.[9,10] We already described several series of potent
p53–MDM2 antagonists, proposing an extended four finger model;
the intrinsically disordered MDM2 N-terminus is ordered
by certain small molecules, which can be obtained by multicomponent
reaction chemistry,[11,12] as shown by cocrystallization.[13,14] Recently, several macrocyclic stapled peptides have been described
with great affinity toward MDM2 and MDMX.[6] ALRN-6924 (Aileron Therapeutics) is currently undergoing phase I
and II clinical trials in patients suffering from solid tumors, lymphoma,
and myeloid leukemias (ClinicalTrials.gov ID: NCT02264613 and NCT02909972).Here, we propose a novel series of nonpeptidic artificial macrocyclic
compounds that inhibit the p53–MDM2 interaction, which might
have a different activity profile from the currently available scaffolds.
Our synthesis strategy is shown in Figure . Based on the Ugi scaffold, we introduced
two terminal ene-functionalities via the carboxylic acid and the isocyanide
component, and we were able to cyclize the compounds by ring closing
metathesis (RCM, Figure ).[15,16]
Figure 1
Macrocyclization strategy to inhibit the p53–MDM2
interaction.
Macrocyclization strategy to inhibit the p53–MDM2
interaction.Our previously introduced
three and four-point pharmacophore models were the basis of the discovery
and development of the current inhibitors. Thus, we used as the starting
point our formerly described α-aminoacylamide (YH300, shown
in cyan sticks, Figure ) with a Ki of 600 nM and its crystal
structure in complex with MDM2 receptor (PDB ID 4MDN).[13] Accordingly, the 6-chloroindole-2-carboxylic acid was used
as an “anchor” in order to mimic the Trp23 amino acid
and constrain the position of other substituents, and three additional
binding sites were defined, Phe19, Leu26, and the induced Leu26 subpocket.
Figure 2
Four-
and three-point pharmacophore modeling of macrocyclic compounds. Left:
YH300 (PDB ID: 4MDN, cyan lines), Tyr67, Gln72, His73, Val93, and Lys94 (blue sticks),
proposed macrocycle to explore the four-pharmacophore point (pink
sticks). Right: SAH-p53–8 stapled-peptide (PDB ID: 3V3B, cyan lines), Tyr67,
Gln72, His73, Val93, and Lys94 (blue sticks), proposed macrocycle
to explore the three-pharmacophore point (pink sticks).
Four-
and three-point pharmacophore modeling of macrocyclic compounds. Left:
YH300 (PDB ID: 4MDN, cyan lines), Tyr67, Gln72, His73, Val93, and Lys94 (blue sticks),
proposed macrocycle to explore the four-pharmacophore point (pink
sticks). Right: SAH-p53–8 stapled-peptide (PDB ID: 3V3B, cyan lines), Tyr67,
Gln72, His73, Val93, and Lys94 (blue sticks), proposed macrocycle
to explore the three-pharmacophore point (pink sticks).Our design of the macrocyclic linker aims to bind
on top of the loop linking α2′ and α1′ helices,
covering a large hydrophobic surface area formed by Tyr67, Gln72,
His73, Val93, and Lys94 (Figure ) and potentially increasing the affinity to the receptor.
Moreover, we reasoned that the addition of a methyl group on the ring
could mimic the tert-butyl group of YH300,[13] likely leading to increased affinity.The retrosynthetic plan of the designed macrocycle 1 foresees (Scheme ) a ring closing metathesis reaction (RCM), followed by a classical
Ugi four-component (U-4CR). The Ugi adduct is formed of the anchoring
6-chloro-3-carboxaldehyde 6, suitably substituted benzylamines 5, and the long-chained aliphaticcarboxylic acids 8 along with the isocyanides 7 incorporating terminal
double bonds.
Scheme 1
Retrosynthetic Plan
Based on a U-4CR and a Ring Closing Metathesis
Aldehyde 6 was synthesized from
the
6-chloro-indole derivative using the Vilsmeier–Haack formylation
reaction.[17,18] For the preliminary structure–activity
relationship (SAR) analysis of the Leu26 and the induced pockets,
the benzylamines 5 were used as commercially available
(5b–h) or obtained via Williamson
ether synthesis (5a). Probing both the Phe19 pocket and
the larger hydrophobic surface area formed by Tyr67, Gln72, His73,
Val93, and Lys94, the isocyanides 7a,b were synthesized
from the corresponding formamides, whereas the carboxylic acids are
commercially available. In particular, the amine 5a with
the oxygen linker that was designed to probe the induced pocket was
synthesized through a Williamson ether synthesis of the protected
4-hydroxybenzylamine with the 3,4-benzyl chloride under basic conditions
(Supporting Information (SI), Scheme ). The isocyanides 7a,b were synthesized from the formamides via the revised
Leuckart–Wallach reaction of the corresponding carbonyl compounds
(SI, Scheme ).[19]
Scheme 2
Synthesis of the
Macrocycle Library Based on the U-4CR/RCM Strategy
Next, we proceeded
in the Ugi four-component
reaction. Equimolar mixture of the substituted benzylamine 5, aldehyde 6, isocyanide 7, and carboxylic
acid 8 in 2,2,2-trifluoroethanol (TFE) was irradiated
at 120 °C for 1 h in a microwave oven yielding compounds 4 (Scheme ). Afterward, we successfully performed the ring closing metathesis
with the second generation of Grubbs catalyst in dichloromethane (DCM)
affording compounds 3 as mixture of isomers (E and Z). Due to the fact that the existence
of the double bond gives rise to two possible isomers and most importantly
reduces the flexibility of the macrocycle, we decided to subject the
mixture to hydrogenation on Pd/C isolating compounds 2. The last step was the ester hydrolysis obtaining the final screening
compounds 1 (Scheme ). Performing a preliminary SAR, we built a small library
of macrocycles of various ring sizes (12, 13, 18, 19, and 24 number
of atoms) targeting the hydrophobic region around Tyr67, Gln72, His73,
Val93, and Lys94. We maintained the anchoring indole group for the
Trp23 and the phenyl group for Leu26 and further explored the induced
pocket with the extended dichlorobenzyloxy moiety.Two complementary
assays based on independent physicochemical principles, fluorescence
polarization (FP) and 1H–15N 2D HSQC
nuclear magnetic resonance (NMR) were used to measure affinity and
to exclude false positive hits. FP assay was employed to determine
the inhibitory affinities (Ki) of the
macrocycles against MDM2 as previously described,[20] and the results are presented in Table .
Table 1
Results of the Evaluation
of Inhibitory Activity of
Macrocyclic Inhibitors toward MDM2 Using FP Assaya
entry
compound
ring size
Ki (μM)
1
1a
12
0.35
2
1b
13
0.32
3
1c
18
0.10
4
1d
18
N/A
5
1e
18
5.3
6
1f
18
N/A
7
1g
18
1.3
8
1h
18
0.08
9
1i
18
1.8
10
1j
19
0.14 (rac-1j)
0.09 {(+)-1ja}
0.70 {(−)-1jb}
11
1k
24
1.9
12
3j
19
0.34
13
3k
24
2.5
14
2a
12
3.0
15
2k
24
N/A
N/A: not active; all measurements performed in duplicates; standard
errors up to 10%.
N/A: not active; all measurements performed in duplicates; standard
errors up to 10%.Examining
both the three- and four-finger
pharmacophore model, it seems that most of the obtained macrocycles
are active toward MDM2, many of them demonstrating an affinity below
1 μM. Although it is a preliminary SAR study, we can conclude
that the ideal ring size for the four-finger model seems to be around
18 (entry 3). The affinity improves while increasing the size from
12 (1a, entry 1) to 13 (1b, entry 2) and
eventually 18 (1c, entry 3) atom ring size. Moreover,
compound 1c has an affinity of 100 nM as a mixture of
diastereomers, whereas the larger 24-membered macrocycle 1k (entry 11) shows decreased activity (1.91 μM).In addition,
we focused on the three-finger-pharmacophore model, which characterizes
the vast majority of the currently available small-molecule MDM2 inhibitors.[6] We synthesized various macrocycles with a different
substitution pattern. The position of halogens on the phenyl group
seems to play a significant role since the para fluoro
(1d, entry 4) or chloro substituted (1f,
entry 6) derivatives are only slightly or not active at all. On the
contrary, the addition of a second chlorine in o-
or m-position influenced the binding mode, the 2,4-dichloro
derivative (1e) showed activity of 5.3 μM (entry
5,) and the corresponding 3,4-dichloro compound 1h (entry
8) displayed an activity of 80 nM. Whereas placing the donor group
−OMe in compound 1g (entry 7) improved the affinity
with 1.3 μM. Nonetheless, m-OMe substitution,
keeping the same ring size of 18 (1i, entry 9), did not
show any significant difference, exhibiting an activity of 1.8 μM.
Due to our previous results,[21] where we
were able to synthesize fluorinated phenyls such as compound 9 with Ki up to 100 nM, we employed
the 3,4,5-trifluorobenzylamine and we synthesized compound 1j (entry 10), which exhibited an interesting affinity as a racemic
mixture of 140 nM. Enantiomeric separation of the racemic mixture
via chiral SFC provided the enantiomers (+)-1ja and (−)-1jb with affinities of 90 and 700 nM, respectively.
As it was expected, the separated enantiomers showed a significant
increase of the activity compared to the racemic mixture.The
existence of the double bond in the macrocycles 3j and 3k, as anticipated, reduced significantly the affinity to
340 nM and 2.5 μM, respectively (entries 12 and 13), compared
with the corresponding hydrogenated compounds 1j (Ki = 140 nM) and 1k (Ki = 1.9 μM). The corresponding esters, 2a and 2k (entries 14 and 15), are orders of magnitude
less reactive or inactive comparing with the acids in accordance with
our previous experience (entries 1 and 11).[11,12,22−25] Interestingly, the acyclic Ugi-adduct 4j was also proven to be practically inactive both in the
ester and acid forms (SI, Table S1, entry
15). Moreover, changing the anchor to a nonsubstituted indole moiety
(compound 10) or to the 3- or 4-phenyl moiety (compound 11), resulted in nearly no activity.The expected ligand-induced
perturbations in 1H–15N 2D HSQC NMR spectra
were indeed observed (Figure ). The 15N-labeled MDM2 was titrated with increasing
concentration of the compound. Since all cross peaks in the MDM2 spectrum
were assigned to particular amino acid residues,[26] it was possible to analyze the interaction within the MDM2/1j complex. Particularly, Val93 is clearly involved in the
interaction, as its cross peak shifted between titration steps for
MDM2/1j molar ratios equal to 2:1 and 1:1. After 1:1
step, the peak remained in the same position. NMR titration also confirmed
the tight binding of 1j, as, e.g., for Arg29, NMR signal
splitting was observed (Figure ), which indicated strong interaction with MDM2 at Kd below 1 μM (and a slow chemical exchange).
Figure 3
Spectrum
of the 15N-labeled MDM2 (blue) superimposed with spectrum
after addition
of 1j in a MDM2/1j molar ratio equal to
2:1 (red) and 1:1 (green). The close-up view shows selected peaks
assigned to Val93 and Arg29. For Arg29, NMR signal splitting indicates
strong interaction at Kd below 1 μM.
Spectrum
of the 15N-labeled MDM2 (blue) superimposed with spectrum
after addition
of 1j in a MDM2/1j molar ratio equal to
2:1 (red) and 1:1 (green). The close-up view shows selected peaks
assigned to Val93 and Arg29. For Arg29, NMR signal splitting indicates
strong interaction at Kd below 1 μM.Three of the macrocyclic compounds
(1c, 1h, and 1j) obtained demonstrated
improved binding affinities (Ki ≈
100 nM) over the lead acyclic molecule, YH300 (Ki = 600 nM). In order to rationalize the tight receptor ligand
interaction, we exploit modeling studies using MOLOC[27] based on the HSQC binding data having as template a known
cocrystal structure (PDB ID: 3TU1)[21] and the small network
analysis using Scorpion software (Figure ).[28] It revealed
the existence of van der Waals interactions of the aliphatic handle
with Tyr67 and His73, the expected alignment of the 6-chloro-indole
moiety of the designed compounds with the p53Trp23 pocket,
whereas the 3,4,5-trifluorophenyl ring occupied the p53Leu26 hydrophobic pocket. Moreover, the π–π interaction
of His96 with the 3,4,5-trifluorophenyl fragment and several van der
Waals interactions with Leu54, Ile61, Phe86, Phe91, Val93, His96,
and Tyr100 are depicted. These findings support our initial hypothesis
of the divergent hydrophobic handle position compared to staple peptides
shown before,[29−34] suggesting a new approach to improve and diversify the extensive
collection of MDM2/X inhibitors.
Figure 4
Small network analysis of 1j (white
sticks) modeled into the MDM2 receptor (PDB ID: 3TU1, cyan sticks): π–π
and van der Waals interactions are shown in orange and yellow dotted
lines, respectively.
Small network analysis of 1j (white
sticks) modeled into the MDM2 receptor (PDB ID: 3TU1, cyan sticks): π–π
and van der Waals interactions are shown in orange and yellow dotted
lines, respectively.To analyze and compare the
physicochemical properties of our newly synthesized macrocycles with
the orally available macrocycle drugs,[35] we plotted molecular weight, clogP, topological polar surface area
(TPSA), number of hydrogen bond donors (HBDs), hydrogen bond acceptors
(HBAs), and number of rotatable bonds (RB) (Figure ). Interestingly, the values of the properties
in most of our macrocycles are set in the appropriate range, demonstrating
the significance of this new strategy to develop potentially oral
bioavailable macrocycles targeting the p53–MDM2 interaction.
In the case of 1c and 1k, clogP goes off
the limits as expected since we were targeting a very lipophilic surface.
However, this will be overcome in the future with a strategy to incorporate
heteroatoms (oxygens) on both the acid and isocyanide linker. After
this initial SAR study, we will in the future synthesize libraries
of novel macrocycles as potent p53–MDM2 inhibitors with higher
diversity and complexity.
Figure 5
Physicochemical properties of the synthesized
macrocycles, compared
with the oral macrocycle marketed drugs on a hexagon radar graph.
The dark gray area contains the low limits of the oral macrocycle
drugs, whereas the light gray shows the highest limits.
Physicochemical properties of the synthesized
macrocycles, compared
with the oral macrocycle marketed drugs on a hexagon radar graph.
The dark gray area contains the low limits of the oral macrocycle
drugs, whereas the light gray shows the highest limits.
Authors: Yijun Huang; Siglinde Wolf; David Koes; Grzegorz M Popowicz; Carlos J Camacho; Tad A Holak; Alexander Dömling Journal: ChemMedChem Date: 2011-09-27 Impact factor: 3.466
Authors: Anna Czarna; Grzegorz M Popowicz; Aleksandra Pecak; Siglinde Wolf; Grzegorz Dubin; Tad A Holak Journal: Cell Cycle Date: 2009-04-16 Impact factor: 4.534
Authors: Anna Czarna; Barbara Beck; Stuti Srivastava; Grzegorz M Popowicz; Siglinde Wolf; Yijun Huang; Michal Bista; Tad A Holak; Alexander Dömling Journal: Angew Chem Int Ed Engl Date: 2010-07-19 Impact factor: 15.336
Authors: Constantinos G Neochoritis; Kan Wang; Natalia Estrada-Ortiz; Eberhardt Herdtweck; Katarzyna Kubica; Aleksandra Twarda; Krzysztof M Zak; Tad A Holak; Alexander Dömling Journal: Bioorg Med Chem Lett Date: 2015-12-15 Impact factor: 2.823
Authors: Constantinos G Neochoritis; Jack Atmaj; Aleksandra Twarda-Clapa; Ewa Surmiak; Lukasz Skalniak; Lisa-Maria Köhler; Damian Muszak; Katarzyna Kurpiewska; Justyna Kalinowska-Tłuścik; Barbara Beck; Tad A Holak; Alexander Dömling Journal: Eur J Med Chem Date: 2019-08-06 Impact factor: 6.514
Authors: Lukasz Skalniak; Aleksandra Twarda-Clapa; Constantinos G Neochoritis; Ewa Surmiak; Monika Machula; Aneta Wisniewska; Beata Labuzek; Ameena M Ali; Sylwia Krzanik; Grzegorz Dubin; Matthew Groves; Alexander Dömling; Tad A Holak Journal: FEBS J Date: 2019-02-16 Impact factor: 5.542
Authors: Eman M M Abdelraheem; Katarzyna Kurpiewska; Justyna Kalinowska-Tłuścik; Alexander Dömling Journal: J Org Chem Date: 2016-09-15 Impact factor: 4.354
Authors: Constantinos G Neochoritis; Maryam Kazemi Miraki; Eman M M Abdelraheem; Ewa Surmiak; Tryfon Zarganes-Tzitzikas; Beata Łabuzek; Tad A Holak; Alexander Dömling Journal: Beilstein J Org Chem Date: 2019-02-20 Impact factor: 2.883
Authors: Atilio Reyes Romero; Angel Jonathan Ruiz-Moreno; Matthew R Groves; Marco Velasco-Velázquez; Alexander Dömling Journal: J Chem Inf Model Date: 2020-12-03 Impact factor: 6.162