T. J. Cogswell‡
*a,
R. J. Lewisb,
C. Sköldc,
A. Nordqvista,
M. Ahlqvistd and
L. Knerr*a
aMedicinal Chemistry, Early Cardiovascular, Renal and Metabolism, BioPharmaceuticals R&D, AstraZeneca, Gothenburg, Sweden. E-mail: thomas.cogswell@acerta-pharma.com; laurent.knerr@astrazeneca.com
bMedicinal Chemistry, Early Respiratory and Immunology, BioPharmaceuticals R&D, AstraZeneca, Gothenburg, Sweden
cDrug Design and Discovery, Department of Medicinal Chemistry, BMC, Uppsala University, P.O. Box 574, SE751 23 Uppsala, Sweden
dDMPK, Early Cardiovascular, Renal and Metabolism, BioPharmaceuticals R&D, AstraZeneca, Gothenburg, Sweden
First published on 6th November 2024
Conformational control of drug candidates to engineer improved potency and ADME properties is an ongoing area of research. Macrocyclic rings tend to offer a greater degree of rigidity than non-cyclised small molecules, and, as a result they are perfect platforms to instil conformational controls. In this study, the difluoroalkoxyphenyl moiety is examined as a tool to alter the conformation of macrocycles. A fluorinated and non-fluorinated macrocyclic matched pair is compared in terms of conformation preferences and related ADME properties. The synthesised macrocycles are found to give similar major conformations exhibiting a trans amide in the macrocyclic backbone. However, for the fluorinated macrocycle, the major trans amide conformation is in equilibrium with a cis amide minor conformation, seen by 1H NMR in a 4:1 ratio of trans/cis. The conformational fits for the minor fluorinated isomer demonstrate the out of plane preference of the difluoroalkoxy system encouraging the amide within the macrocycle backbone to adopt a cis conformation. The fluorinated macrocycle was less metabolically stable compared to the non-fluorinated, postulated to be a result of the interconversion of trans amide to the cis amide, which potentially could be more readily metabolised.
To capitalize on the preorganisation of macrocyclic structures, there is a significant interest into ways to bias and alter the conformation of macrocycles to optimize their target binding interactions through a rational design approach.8 Common approaches include the incorporation of stereocenters,9 the addition of alkyl groups to backbone atoms,10 and the installation of both saturated and unsaturated heterocycles into the macrocyclic backbone.11 One interesting approach takes advantage of the stereoelectronic interactions that polarized C–F bonds have on neighbouring groups, where orbital alignments favour stabilized conformations.12 Incorporation of fluorine atoms has previously been identified as a means to change conformational preference in many functional systems such as organocatalysis,13 pharmaceuticals,14 and liquid crystals.15
In a seminal publication in 2014, Hunter et al. reported the use of the fluorine gauche effect to demonstrate, for the first time, that stereoselective fluorination could alter the conformation of a peptidic macrocycle.16 Each synthesised isomer of the vicinal fluorinated Unguisin A occupied a unique low energy conformation which differed from the parent non-fluorinated compound. In 2016, O'Hagan and co-workers found that when installing a CF2-group into different positions of the aliphatic backbone of a family of musk lactones, the macrocyclic ring was constrained into a rectangular conformation with the CF2-group always situated at a corner site.17
In a non-macrocyclic system, Durley and co-workers reported an intriguing fluorination effect whilst developing a library of cholesteryl ester transfer protein (CETP) inhibitors.18 When compared to the alkoxyphenyl analogue 1, the incorporation of a fluorinated side chain in 2 led to an 8-fold increase in potency which was attributed to a more optimal configuration adopted within the enzyme active site (Fig. 1A). After studying the conformation through ab initio calculations, they confirmed that the non-fluorinated alkyl chain of 1 aligned with the plane of the aromatic ring whereas the fluorinated chain occupied an orthogonal conformation.19 This rotation of the OCF2 bond allowed the lone pairs of the oxygen to donate into the low lying antibonding orbitals leading to a stabilization of the orthogonal conformation (Fig. 1B). Since this first example, O'Hagan reported that the same conformational effect also occurs when a sulfur linker replaces the oxygen.20
Fig. 1 (A) The structures of CETP inhibitors reported by Durley and co-workers.19 The non-fluorinated analogue 1 has a IC50 of 1.6 μM while the fluorinated analogue 2 has an IC50 of 0.2 μM. The values were determined in a buffer transfer assay. (B) Hyperconjugative π–σ* C–F interaction rationale for the reduced in-plane conformational preference of a fluoroalkoxyaryl group.19 |
A survey of available crystallographic data for difluoromethoxyaryl-containing structures supported by density functional theory (DFT) calculations confirmed that the orthogonal arrangement of two endo fluorines was one of the two major conformations adopted by this module.21 An endo–exo arrangement was also predicted to have a similar energy in the gas phase leading to postulations that, due to this low barrier to rotation, this module has the potential to access different conformations when subjected to different micro-environments.22
Intrigued by the reported impact of the gem-difluorination of alkoxyphenyl group in acyclic systems, we decided to expand these studies by incorporating this fluorinated module into a macrocyclic core. While there is some exemplification of this chemistry in recent patent literature,23 there is still very little data reported on the impact of that design on the macrocycle properties. Therefore, we aimed to expand the toolbox of small, fluorinated units available to increase the conformational diversity and adaptability of macrocyclic scaffolds. Conformational changes through fluorination has been shown to alter physicochemical properties such as polarity, permeability, and solubility, therefore, an in-depth study into these relationships was also conducted.24
Fig. 2 The structure of a BACE-1 inhibitor 3 developed by researchers at Merck for the treatment of Alzheimer's disease (left).25 The target model macrocyclic structures 4a and 4b to be studied herein (right). |
The fluorinated and non-fluorinated macrocycles 4a and 4b were prepared using the synthetic route reported in Scheme 1. The 1,1-difluorinated allylated phenol starting material 5a was synthesised in a previous publication from our laboratories utilizing methodology developed by Ichikawa et al.26,27 The starting materials (5a and 6) were combined through a Heck reaction followed by hydrogenation providing 8a in good yield. Ester hydrolysis and Boc deprotection gave the precyclisation precursor, that upon macrolactamisation yielded the desired macrocyclic scaffold 4a. Using the same synthetic route, the non-fluorinated match pair macrocycle 4b was also prepared.
For macrocycle 4a, two unexpected features of the proton NMR spectrum were immediately apparent (Fig. 3). Firstly, two conformations were visible in an approximate ratio of 4:1 with slow exchange on the NMR timescale. Secondly, the minor conformation exhibited a highly shielded aromatic proton at 4.59 ppm. This signal was 2.2 ppm more shielded than the equivalent proton in the major conformation. Most other protons exhibited differences in chemical shift of 0.2 to 0.9 ppm between the conformations. This conformational signature suggests a significant change in conformation between major and minor forms.
Fig. 3 The stacked 1H NMR spectra in DMSO-d6 of the non-fluorinated macrocycle 4b (top) and fluorinated macrocycle 4a (bottom). |
In contrast, the proton spectrum for non-fluorinated macrocycle 4b did not clearly exhibit slowly exchanging conformations, though a ROESY NMR experiment did reveal the presence of a second conformation (visible as exchange cross peaks from certain protons – see ESI†). The ratio between the major and minor forms was 200:1, a significant difference to that observed in 4a.
Conformational searches for both 4a and 4b were performed that contained examples of both cis and trans amides. After QM optimization, data was fit in the MSpin program using proton–proton NOE distances, coupling constants and calculated proton chemical shifts. The use of chemical shifts was especially important for the minor conformation as only a few proton–proton NOE distances could be accurately determined. Proton-fluorine NOEs were used qualitatively to distinguish between the different conformational fits found for 4a.
For the major conformation of 4a (4amajor), two conformations were found, approximately equally distributed, accounting for about 90% of the population (Fig. 4). The two conformations differ from each other primarily in the relative orientations of the two aromatic rings.
Fig. 4 Conformations 4amajor_#1 and 4amajor_#2 from Mspin making up 90% of the population of 4amajor. |
For the minor conformation of 4a (4aminor), the same procedure was followed to yield a fit composed of two conformers 4aminor_#1 at 65% and 4amajor_#2 at 30% (Fig. 5). Interestingly, the conformers both bear the cis amide in place of the trans, demonstrating a significantly different conformation compared to 4amajor. The highly shielded proton seen in the NMR spectrum (Fig. 3), was found to make an edge/face interaction with the second aromatic ring accounting for the unusual shielding.29
Fig. 5 The conformational fits from Mspin for 4aminor. Two cis amides 4aminor_#1 (left) and 4aminor_#2 at 30% (right). |
A similar approach was followed to find two clusters of similar conformations for the non-fluorinated macrocycle 4b. The two clusters are shown in Fig. 6. On the left are two conformations differing only in the OH bond orientation accounting for 55% of the population and on the right, three conformations accounting for 25% of the population which differ in OH orientation and the amide bond portion of the macrocycle.
The results from Mspin were compared with conformations selected by comparing the sum of the absolute difference between the experimentally derived 1H–1H NOE distances for 4a and 4b respectively (See ESI†) and the corresponding atom distance of a calculated conformational ensemble after QM optimization.
A combination of two conformations was proposed for 4amajor since both fluorine substituents displayed a strong 1H–19F NOE to H11. These two conformations agreed with the solutions found by Mspin (Fig. 4). For 4b an overall best fit was found for a conformation comprising a trans amide bond and the ether substituent in the plane of the phenyl ring, this also best reflected one of the suggested solutions from Mspin (Fig. 6).
The fluorinated analogue 4a is expected to prefer that the O–CF2 bond is out of plane of the bonded aromatic ring. However, this bond is constrained by being part of a macrocyclic ring, thus moving this bond out of plane likely introduces higher energy conformational elements elsewhere in the macrocycle. These conformational compromises can be seen born out in conformations 4amajor_#1, 4amajor_#2, 4aminor_#1 and 4aminor_#2. Conformation 4amajor_#2 is similar to the 4b conformations in that the O–CF2 bond is almost co-planar with the aromatic ring and the ideal chain geometry is only slightly distorted (maximum deviation 18°). In conformation 4amajor_#1 however, the O–CF2 bond is 32° out of plane from the aromatic ring but the ideal chain geometry is more distorted (maximum deviation 37°). When the amide is in a cis geometry (conformations 4aminor_#1 and 4aminor_#2), the O–CF2 bond can move further out of plane (60° and 46°) and the ideal chain geometry is less distorted (maximum deviation 23°). This comes at the expense of some greater deviation from planarity of the amide bond (36° and 27°). 4b adopts a conventional conformation with a trans amide bond and has torsion angles of bonded sp3 atoms (syn, anti, gauche) deviating by no more than 12° from an ideal chain (+60°, 180°, −60°). The macrocyclic O–CH2 bond is within 6° of co-planarity of the aromatic ring except in conformation 4b_#5 when it is 15° out of plane. The amide bond is distorted from planarity by 12–17°. Full details of the torsions observed can be seen in Table S3.†
Using the ratios of trans to cis conformer for 4a and 4b, it is possible to estimate an energy penalty of at least 10 kJ mol−1 for positioning the fluorine atoms approximately in plane with the aromatic ring compared to the proton atoms. This assumes only that the remote fluorine proton exchange has no direct effect on the energetics of the amide bond.
Fluorinated macrocycle 4a | Non-fluorinated macrocycle 4b | |
---|---|---|
logD7.4 | 2.4 | 2.4 |
ePSA (Å2) | 70 | 69 |
Human hepatocytes metabolic stability t1/2 (min) | <2.3 | 10 |
Human liver microsomes metabolic stability t1/2 (min) | <2.3 | 31 |
Chemical stability | 11 | >72 |
Aq. pH1 t1/2 (h) | ||
Chemical stability | >72 | >72 |
Aq. pH7.4 t1/2 (h) | ||
Chemical stability | >72 | >72 |
Aq. pH10 t1/2 (h) | ||
Solubility (μM) | 21 | 243 |
Intrinsic permeability | 96 | 74 |
Caco AB pH 6.5 (inhibition) Papp (1 × 10−6 cm s−1) |
A difference in the metabolic stability was observed between the two systems. In both human hepatocytes and human liver microsomes, the half-life (t1/2) of fluorinated 4a was below a measurable time point whereas the non-fluorinated 4b had a t1/2 of 10 and 31 minutes respectively.
At first glance, this result seems to be counter intuitive as fluorination has been widely reported to improve the metabolic stability of drug candidates.31 To look in more detail, a metabolic stability study was undertaken to identify the metabolites in both the fluorinated and non-fluorinated examples upon incubation in human hepatocytes (Schemes 2 and 3). The amide of 4a was found to be rapidly hydrolysed to the open compound as a major metabolite. Whereas 4b showed little amide hydrolysis and instead was predominantly dealkylated at the phenolic site leading to the phenol derivative as the main metabolite. As expected, the fluorinated macrocycle was completely protected from this dealkylation metabolic pathway. We postulate that the presence of a significant cis amide population due to the conformational impact of the gem-difluorinated system is a major driver in the promotion of metabolism through amide cleavage. This represented a dramatic switch of metabolism pathway upon fluorination and whilst the magnitude is surprising in this case similar effect has been previously reported in the literature.32
To understand the chemical stability profiles of both macrocycles, we complemented the metabolization study with hydrolytic stability studies at pH 1, 7.4 and 10. We could show that, except for the lesser stability of fluorinated macrocycle in acidic conditions (half-life 11 hours at pH 1), both scaffolds were similarly stable with half-lives above 3 days in all conditions tested. Based on pKa measurements for two benzoic acids possessing the fluorinated and non-fluorinated allyloxy groups (see ESI†) we could also verify that the inductive electron withdrawing effect of the fluorinated analogue was modest thus strengthening the evidence for impact on chemical stability of this fluorinated module being limited.
Additionally, we could show that the fluorinated analogue peculiar metabolic stability profile was specific to the macrocyclic scaffold by studying a matched pair of related linear analogues (see ESI†). Interestingly this demonstrated that while the linear non-fluorinated analogue showed comparable stability to its macrocyclic counterpart, the stability of the linear fluorinated analogue was considerably higher than the fluorinated macrocycle based on their half-lives upon incubation with human microsomes or human hepatocytes. The improved stability for the linear fluorinated analogue could not be attributed to a change in lipophilicity as both fluorinated analogues had a logD7.4 of 2.4.
The physicochemical and ADME properties for 4a and 4b were examined, revealing a difference in metabolic stability in both human hepatocytes and human liver microsomes between the two macrocycles. Metabolite identification studies revealed that for 4a, amide cleavage was the predominant metabolic pathway, whereas this was a minor pathway for 4b.
This is the first time this gem-difluorinated alkoxyphenyl conformational effect has been demonstrated experimentally in a macrocyclic system and, interestingly, it was found to have a profound effect on the properties of the compound. The fluorination pattern studied here represents another potential conformational design tool to modulate the conformer population of macrocyclic scaffolds and optimize their interaction with biological targets. Future work to install this group into a diverse array of macrocyclic rings could shed further light on how this effect could be harnessed to tune the conformations of macrocycles. The inclusion of this motif into an active medicinal chemistry regime where target-binding and target-selectivity data are generated represents an important future goal for this area of research.
Footnotes |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d4sc05424e |
‡ Current address: Acerta Pharma B.V., A Member of the AstraZeneca Group, Kloosterstraat 9, 5349 AB, Oss, The Netherlands. |
This journal is © The Royal Society of Chemistry 2024 |