Lijia
Chen
a,
Jay
Chauhan
a,
Jeremy L.
Yap
a,
Christopher C.
Goodis
a,
Paul T.
Wilder
bc and
Steven
Fletcher
*ac
aDepartment of Pharmaceutical Sciences, University of Maryland School of Pharmacy, 20 N. Pine St., Baltimore, MD 21201, USA. E-mail: steven.fletcher@rx.umaryland.edu
bUniversity of Maryland School of Medicine, 20 S. Greene St., Baltimore, MD 21201, USA
cUniversity of Maryland Greenebaum Cancer Center, 20 S. Greene St., Baltimore, MD 21201, USA
First published on 27th October 2022
The anti-apoptotic protein MCL-1, which is overexpressed in multiple cancers, is presently a focus for the development of targeted drugs in oncology. We previously discovered inhibitors of MCL-1 based on 1-sulfonylated 1,2,3,4-tetrahydroquinoline-6-carboxylic acids (“1,6-THQs”). However, with the nitrogen atom constrained in the bicyclic ring, we were unable to modify the alkyl portion of the tertiary sulfonamide functionality. Moreover, the introduction of additional functional groups onto the benzene ring portion of the THQ bicycle would not be trivial. Therefore, we elected to deconstruct the piperidine-type ring of the 6-carboxy-THQ lead to create a new 4-aminobenzoic acid scaffold. Given its simplicity, this permitted us to introduce diversity at the sulfonamide nitrogen, as well as vary the positions and substituents of the benzene ring. One of our most potent MCL-1 inhibitors, 6e-OH, exhibited a Ki of 0.778 μM. Heteronuclear single quantum coherence experiments suggested 6e-OH bound in the canonical BH3-binding groove, with significant perturbations of R263, which forms a salt bridge with MCL-1's pro-apoptotic binding partners, as well as residues in the p2 pocket. Selectivity studies indicated that our compounds are dual inhibitors of MCL-1 and BCL-xL, with 17cd the most potent dual inhibitor: Ki = 0.629 μM (MCL-1), 1.67 μM (BCL-xL). Whilst selective inhibitors may be more desirable in certain instances, polypharmacological agents whose additional target(s) address other pathways associated with the disease state, or serve to counter resistance mechanisms to the primary target, may prove particularly effective therapeutics. Since selective MCL-1 inhibition may be thwarted by overexpression of sister anti-apoptotic proteins, including BCL-xL and BCL-2, we believe our work lays a solid foundation towards the development of multi-targeting anti-cancer drugs.
Initial BCL-2 inhibitors were pan-inhibitors, exhibiting similar binding affinities to multiple anti-apoptotic BCL-2 members, such as the pyrogallol derivative TW37.9 Intense drug discovery programs have guided the design of selective inhibitors of these proteins,5,6 including the FDA-approved BCL-2 inhibitor venetoclax,10 the BCL-xL inhibitor A-1331852 (ref. 11) and, more recently MCL-1 inhibitors, such as the indole-2-carboxylic acids12,13 and -2-acylsulfonamides14 (Fig. 1), typified by 1a and 1b, a family of sulfonylated aminonaphthalenes,15e.g. UMI-77 (2), and the macrocyclic compound AZD5991 (3),16 which is currently in clinical trials for haematological malignancies, although one trial remains suspended owing to an adverse cardiac event.17 Our laboratory has also discovered MCL-1 inhibitors based on, among others,18–21,22 a racemic 1,2,3,4-tetrahydroquinoline-3-carboxylic acid (3-THQ) scaffold, such as (±)-4.23 A significant majority of reported MCL-1 inhibitors present a carboxylic acid to bind R263, as well as hydrophobic groups to bind the p2 pocket (and sometimes additional pockets).4–6,8 With selective inhibitors in hand, it becomes possible to interrogate the biological implications of the inhibitions of specific targets. However, ultimately a therapeutic regimen may comprise a drug “cocktail” in which multiple drugs bind multiple targets that together result in a more efficacious treatment of the disease.24 The major determinants of resistance to the BCL-2 drug venetoclax are overexpression of MCL-1 and/or BCL-xL.25 Accordingly, combination therapies of venetoclax with MCL-1 or BCL-xL inhibitors are being investigated.25 An alternative strategy to combination therapy, or polypharmacy, is the development of single agents that are purposely designed to bind multiple targets, i.e. polypharmacology. Polypharmacy and polypharmacology have been reviewed extensively elsewhere, and each carries its own advantages and disadvantages.26 The most significant concern in polypharmacology is the inability to alter the dosing of one pharmacophore over the other, yet that might be outweighed by enhanced additive or synergistic effects. Herein, we describe our continued efforts into the development of inhibitors of the BCL-2 anti-apoptotic proteins, and report the discovery of dual MCL-1/BCL-xL inhibitors based on sulfonylated aminosalicylic acids.
Previously, we reported the 3-THQ derivative (±)-4 as an inhibitor of MCL-1 (Ki = 0.12 μM).23 This compound presented a carboxylic acid to engage R263 through a salt bridge, and a 4-chloro-3,5-dimethylphenoxyphenyl group projected from the nitrogen atom of a sulfonamide linker, which was predicted to probe deep into the p2 pocket. The 3-THQ scaffold was then isomerized into a 1,2,3,4-tetrahydroquinoline-6-carboxylic acid (6-THQ) scaffold to remove the chiral center, as exemplified by compound 5 (Fig. 2), although this led to a significant decrease in MCL-1 binding affinity (Ki = 2.73 μM).27 In the structure of the 6-THQ scaffold, the nitrogen atom is fused in the THQ ring, restricting the ability to alter the N-alkyl component, and modification of the phenyl ring is not trivial. To overcome these barriers to further optimization of the 6-THQ series, we decided to re-engineer the 6-THQ scaffold by deconstructing the piperidine-type ring, thereby transforming the core into a 4-aminobenzoic acid scaffold (e.g.6e) as shown in Fig. 2. This particular scaffold offers an acyclic nitrogen atom at the same para position compared to the 6-THQ scaffold, but now it may be readily doubly-functionalized. Although this modification may be deleterious owing to a loss of rigidity, benefits towards structure–activity relationship studies include (i) a wide range of commercially available 4-aminobenzoic acids, (ii) the symmetry renders the introduction of additional functionality more controlled, and (iii) the location of the nitrogen atom may be easily shifted.
Fig. 2 Dismantling the 6-THQ-based MCL-1 inhibitor 5 into a more flexible and more optimizable ligand. |
A GOLD docking experiment was performed to investigate the possibility of compound 6e binding to the MCL-1 protein. The small-molecule was first MM2 energy minimized, whilst the MCL-1 protein crystal structure (4HW2) was processed to remove the ligand and water molecules. The binding site was defined as a 10 Å radius sphere around residue M231, which is located at the entrance to the p2 pocket. No further constraints were used. One high-scoring result is shown in Fig. 3. As illustrated, the 4-chloro-3,5-dimethylphenyl moiety of compound 6e successfully docked into the p2 binding pocket of MCL-1, there are potential electrostatic interactions and hydrogen bonds between the carboxylic acid and sulfonamide of 6e with R263, and there appears to be a hydrophobic interaction between the isobutyl group of 6e and a shallow hydrophobic pocket formed by A227, F228, M231 and F270, as seen elsewhere.12
Fig. 3 Predicted binding mode of 4-sulfonamidobenzoic acid derivative 6e with MCL-1 (PDB ID: 4HW2) using GOLD software: (a) the MCL-1 protein is shown in grey helical structure with key residues in the p2 pocket shown in sticks (carbon: gray; nitrogen: blue; oxygen: red; sulfur: yellow), 6e is shown in stick structure (carbon: green; nitrogen: blue; oxygen: red; sulfur: yellow); hydrogen bonds are shown in black dashed lines; (b) as in (a) but with the MCL-1 protein shown with a solvent-accessible surface (hydrophobic: grey; basic: blue; acidic: red; sulfur: yellow), and the locations of the pockets p1–p4 are indicated. |
Encouraged by the docking result, a focused library of compounds was synthesized according to Scheme 1. Briefly, 4-nitrobenzoic acid (7) was esterified with methanol and sulfuric acid to yield methyl ester 8. The nitro group of 8 was reduced with tin(II) chloride to deliver aniline 9, which underwent reductive amination with isobutyraldehyde to afford compound 10. Subsequently, sulfonylations of 10 generated the library of compounds 11a–11d, which were then saponified with lithium hydroxide to yield the 4-sulfonamidobenzoic acid final molecules 6a–6d. Alternatively, compound 11d was further elaborated in an SNAr reaction with 4-chloro-3,5-dimethylphenol, followed by saponification to give final compound 6e.
The binding affinities of the final molecules were evaluated using a well-established fluorescence polarization competition assay (FPCA) where the compounds were titrated into a mixture of MCL-1172-327 and FITC-BAK BH3 peptide to monitor their abilities to disrupt this interaction, as described previously.18–20,22,23,28,29 Although the unsubstituted phenyl analogue 6a and the bioisostere 6d demonstrated no inhibitory activity, building off the para-position as in biphenyl analogue 6c resulted in moderate inhibition of MCL-1 with a Ki of 16.7 μM. In contrast to 6c, the naphthyl analogue 6b was inactive (Ki > 500 μM). The data indicate that the geometrical shape of rigid hydrophobic moieties may have significant influence on binding affinities. Analogous to our discovery in the 3-THQ and 6-THQ SAR library studies,23 further extension of the R group with the 4-chloro-3,5-dimethylphenoxy moiety yielded compound 6e as the most potent analogue in this batch, whose binding affinity (Ki = 2.06 μM) was ≥10-fold stronger than 6c. Taken together, these findings suggest that although the binding domain for the sulfonamide R group – possibly the p2 pocket – is not wide, it is deep, which is consistent with previous reports.12
With further investigation of the docking result in Fig. 3, we suspected that installing extra polar functional groups ortho to the carboxylic acid may yield additional interactions with either R263 and/or T266. Accordingly, we transformed 4-aminobenzoic acid derivative 6e into its corresponding 4-aminosalicylic acid derivative “6e-OH” and performed another GOLD docking experiment; one high scoring result is depicted in Fig. 4. In this particular predicted conformation, the carboxylic acid forms a salt bridge with R263, as well as a hydrogen bond with T266. Additionally, the hydroxyl successfully recognized and hydrogen bonded to R263. Therefore, we hypothesized that 4-sulfonamidosalicylic acids, such as 6e, would be stronger binders compared to the corresponding 4-sulfonamidobenzoic acids (Table 1).
Fig. 4 Predicted binding mode of 4-sulfonamidosalicylic acid 6e-OH with MCL-1 (PDB ID: 4HW2) using GOLD: the MCL-1 protein is shown with a solvent-accessible surface (hydrophobic: grey; basic: blue; acidic: red; sulfur: yellow), and the locations of the pockets p1–p4 are indicated; 6e-OH is shown in stick structure (carbon: green; nitrogen: blue; oxygen: red; sulfur: yellow). |
To test this hypothesis, we designed and synthesized a library of 4-sulfonamidosalicylic acid analogues according to Scheme 2. Esterification of 4-aminobenzoic acid (12) yielded methyl ester 13, whose phenolic hydroxyl was protected to afford compound 14. Reductive aminations of aniline 14 with isobutyraldehyde, benzaldehyde or cyclopentanone installed the R1 group, furnishing compounds 15a–15c. Sulfonylations (introduction of the R2 group) of the secondary amines 15 delivered sulfonamides 16, which were then saponified followed by TFA-mediated debenzylation to generate the target molecules 17.30 Alternatively, compounds 16ad, 16bc and 16cc were further elaborated with different substituted phenols R3OH through SNAr chemistry to afford biphenyl ethers 16ah–16ak, 16bd and 16cd. Subsequently, these compounds were subjected to the same deprotection protocol as before to yield the target molecules 17ah–17aj, 6e-OH, 17bd and 17cd. Additionally, compound 14 was directly sulfonylated by 4-phenoxybenzenesulfonyl chloride to yield, after saponification, compound 17da′, which was subsequently debenzylated to furnish target molecule 17da.
Using the same FPCA experiment as described earlier, the MCL-1 binding affinities of the 4-sulfonamidosalicylic acid derivatives were determined (Table 2). Excitingly, all of these compounds exhibited stronger MCL-1 binding affinities compared to their corresponding 4-sulfonamidobenzoic acid counterparts in Table 1, demonstrating the significant contribution from the addition of the OH group, as predicted by the molecular modeling (Fig. 4). Furthermore, increasing the R2 size and hydrophobicity yielded stronger inhibitors, illustrated by the decreasing Kis of almost 50-fold from the phenyl 17aa (Ki = 192 μM), 4-flurophenyl 17ad, 4-methylphenyl 17af, 4-bromophenyl 17ae, 2-naphthyl 17ab to 4-biphenyl 17ac (Ki = 4.11 μM). Further expansion of the R2 moiety from the para position of the benzene sulfonamide to deliver biphenyl ethers produced the group of tighter binders 17ag–17aj and 6e-OH, indicating a larger and kinked R2 group is desirable, affording greater hydrophobic interactions with the BH3-bindng groove, likely enabling deeper access to the p2 pocket. Again, the most potent analogue 6e-OH was equipped with the 4-chloro-3,5-dimethylphenoxy group, inhibiting MCL-1 with a Ki of 778 nM.
In addition to varying the R2 sulfonyl moiety, we modified the R1 group to determine its impact on binding to MCL-1. In each case, consistent increases in binding affinities from naphthyl to 4-biphenyl to 4-chloro-3,5-dimethylphenoxyphenyl were observed within a family of R1 compounds, with R1 = cyclopentyl the worst, and R1 = benzyl the best. Lastly, it is noteworthy that completely eliminating the hydrophobic R1 group confirmed its significance: 17da (R1 = H), Ki = 13.7 μM; 17ag (R1 = isobutyl), Ki = 5.77 μM.
Upon dismantling the 6-THQ ring, as well as relinquishing the constraint on the identity of the nitrogen alkyl appendage, there was also now no restriction on the relative position of the nitrogen atom. Inspection of previous docking results suggested that the optimal relative locations of the carboxylic acid and the sulfonamido groups were para and meta. Thus, we designed 5-sulfonamidosalicylic acids, such as 6e-OH-5, to ascertain which substitution pattern was optimal. As illustrated in Fig. 5, 6e-OH-5 favorably fitted into the MCL-1 binding pocket and developed similar polar interactions with R263 and T266 compared to its predecessor 4-sulfonamidosalicylic acid 6e-OH (Fig. 4). Additionally, the hydrophobic 4-chloro-3,5-dimethylphenyl moiety probed into the p2 pocket of MCL-1. Following the docking study, a library of 5-aminosalicylic acids was designed and synthesized according to Scheme 3 to validate the in silico model. In addition, the function of the hydroxyl group was re-assessed in this new scaffold by deleting it to generate 3-sulfonamidobenzoic acid derivatives.
Fig. 5 Predicted binding mode of 5-sulfonamidosalicylic acid 6e-OH-5 with MCL-1 (PDB ID: 4HW2) using GOLD: the MCL-1 protein is shown with a solvent-accessible surface (hydrophobic: grey; basic: blue; acidic: red; sulfur: yellow), and the locations of the pockets p1–p4 are indicated; 6e-OH-5 is shown in stick structure (carbon: green; nitrogen: blue; oxygen: red; sulfur: yellow). |
The target molecules were synthesized similarly to in Schemes 1 and 2, with the R1 group constrained as isobutyl. Briefly, esterifications of 3-aminobenzoic acid (18a) and 5-aminosalicylic acid (18b) afforded compounds 19a and 19b, respectively. Reductive aminations with isobutyraldehyde furnished compounds 20a and 20b, which were further functionalized with different sulfonyl chlorides (RSO2Cl) to yield sulfonamides 21aa–21ad, and 21ba–21bc. These compounds were then subjected to saponifications to yield the final molecules 22aa–22ad, and 22ba–22bc. Also, compound 21ad of the 3-sulfonamidobenzoic acid family was further elaborated by an SNAr reaction with 4-chloro-3,5-dimethylphenol, followed by saponification to afford compound 22af. In the case of the 5-sulfonamidosalicylic acid family, compound 21bd was first O-benzylated to give compound 21bd′ before undergoing an SNAr with 4-chloro-3,5-dimethylphenol to afford 21be, and was then subjected to saponification followed by TFA-mediated debenzylation to yield final compound 6e-OH-5.30
The SAR data for both 3-sulfonamidobenzoic acid and 5-sulfonamidosalicylic acid derivatives are reported in Table 3. The pronounced differences of Kis between the benzoic acids (left) and the salicylic acids (right) demonstrate the compelling benefits generated by the hydroxyl group as seen earlier with the 4-aminobenzoic/salicylic acids by comparing Tables 1 and 2. For example, 4-biphenyl 22ac was inactive, yet conversion of the scaffold into the corresponding salicylic acid 22bc resulted in a moderate inhibitor of MCL-1 (Ki = 8.70 μM). Intriguingly, in the benzoic acid class, the 2-naphthyl analogue 22ab (Ki = 61.7 μM) exhibited a stronger binding affinity than the 4-biphenyl counterpart 22ac (Ki > 500 μM), which is in stark contrast to the relationships between the corresponding 4-aminobenzoic acid derivatives 6b and 6c, respectively (Table 1). This suggests that the shifted nitrogen may have changed the binding orientation of the scaffold, potentially now directing the 4-biphenyl moiety at an unfavorable angle towards MCL-1, while guiding the 2-naphthyl group in a favorable direction towards the BH3-binding groove, possibly the p2 pocket. The most potent compounds 22af and 6e-OH-5 carry the 4-chloro-3,5-dimethyl-4-chlorophenoxy group, with Ki = 1.45 μM and 586 nM, respectively, as with other series of compounds. Once again, in light of the data for biphenyls 22ac, in particular, and 22bc, it appears that the ether oxygens of 22af and 6e-OH-5 are pivotal in realizing more effective inhibition of MCL-1, probably through the inherent kink of the ether functional group that promotes better interactions with the terminal (and now more substituted) aryl ring. Although the library of 5-sulfonamidosalicylic acids is limited, the data suggest that the 4-sulfonamidosalicylic acids and the 5-sulfonamidosalicylic acids isomers are comparable in MCL-1 inhibitory activity.
Fig. 7 Plot of chemical shift of MCL-1 protein upon 6e-OH binding, red columns are the amino acid residues that were significantly perturbed (Δ chemical shift >0.3 ppm). |
Fig. 8 NMR chemical shift perturbations of the MCL-1/6e-OH complex mapped onto MCL-1 crystal coordinates (PDB ID: 4HW3) in PyMOL. Residues experiencing chemical shifts of at least 0.3 ppm are shaded red. |
Footnote |
† Electronic supplementary information (ESI) available: All experimental details, and 1H and 13C NMR spectra. See DOI: https://doi.org/10.1039/d2md00277a |
This journal is © The Royal Society of Chemistry 2023 |