Thomas
Toupy‡
a,
Christopher
Kune‡
b,
Kristof
Van Hecke
c,
Loïc
Quinton
*b and
Jean-Christophe M.
Monbaliu
*a
aCenter for Integrated Technology and Organic Synthesis, MolSys Research Unit, University of Liège, B-4000 Liège (Sart Tilman), Belgium. E-mail: jc.monbaliu@uliege.be; Web: http://www.citos.uliege.be
bMass Spectrometry Laboratory, MolSys Research Unit, University of Liège, B-4000 Liège (Sart Tilman), Belgium
cXStruct, Department of Chemistry, Ghent University, Krijgslaan 281-S3, B-9000 Ghent, Belgium
First published on 22nd November 2021
The peculiar isomer-selective reduction of (α)-hydroxyiminophosphonates (oxime isomers) into (α)-hydroxyaminophosphonate (hydroxylamine) derivatives is presented. A library of 16 (α)-hydroxyiminophosphonates is prepared and studied via a unique multifaceted approach involving the interplay of NMR, XRD, MS, IM-MS and computational chemistry techniques. The combination of NMR, XRD and HPLC enables the seamless separation, identification and quantification of the oxime isomers (E/Z). Tandem MS (MS/MS) enables the determination of the fragmentation patterns for both isomers. Collision energy breakdown curves highlight the order of apparition of the fragments as well as their related energy of fragmentation, demonstrating that the strength of the C–P bond in the Z isomers is much weaker than in the E isomers. Computational chemistry demonstrates that favorable protonation site is isomer-dependent with the phosphonate moiety being the favorable protonation site for the E isomers, while protonation occurs preferentially on the amino moiety for Z isomers regardless of the phosphite source. The combination of these various methods led an unprecedented level of characterization of oxime isomers, providing a better uderstanding of the isomer-dependent behavior of (α)-hydroxyiminophosphonates.
Fig. 1 Comparison between the native amide peptide bond and the principal moieties involved in the backbone of phosphorus analogs of peptides, commonly called phosphonopeptides. |
Since their introduction in 1975 by Martell, they have gained significant attention as bioisosters of native peptides.1 Indeed, introduction of an aminophosphonic acid into a peptide sequence enables several structural and physico-chemical modifications. For instance, the greater lipophilicity of phosphonate compounds compared to their carboxylic analogs is a distinct advantage resulting in an increased overall biodisponibility.2 Hence, numerous examples of phosphonopeptide derivatives have been prepared for their potency as biological agents, including antibacterial agent, antithrombotic agent, inhibition of proteases (metallo-, serine, aspartic, cysteine and HIV) as well as inhibitor of Human Collagenase.3
Due to their renowned interest, numerous synthetic routes have been reported for the preparation of (α)-aminophosphonic acid, such as the three component Kabachnik–Fields reaction,4–10 the Pudovik reaction11–15 or the reduction of (α)-hydroxyiminophosphonate (aka (α)-oxime phosphonate) or (α)-hydroxyaminophosphonate16–19 (Fig. 2a). The main advantage of the latter is that it enables the preparation of the desired product as a free amine compared to the former ones. Two atom economic methods for accessing (α)-aminophosphonic acids and derivatives relies either on (a) a direct nucleophilic amination of (α)-halogenated alkylphosphonates or on (b) the reduction of an intermediate (α)-iminophosphonate derivative. It is worth mentioning that both methods proceed through the formation of an intermediate (α)-hydroxyaminophosphonate (Fig. 2b).
Hydroxyaminophosphonate and hydroxyamino derivatives have found a wide range of synthetic utilities besides being potential intermediates en route toward (α)-aminophosphonic acids. In fact, hydroxyamino derivatives act as versatile amino source for the preparation of various nitrogen-containing molecules.20–23 Besides their common use as N-nucleophiles, they can also serve as electrophilic aminating reagents upon conversion of the hydroxyl moiety by a more potent leaving group.24–26 They have also been used as nitrogen-radical precursors in visible-light photochemistry27 as well as a more suited, non-explosive surrogate for 1-hydroxybenzotriazole (HOBT), which has been the benchmark of coupling agents in peptide chemistry.28 More recently, as antibiotic resistance has become a worldwide issue, hydroxyamino derivatives have shown promising results in the search of new antibacterial agents.29 In light of these numerous utilities, various methods have been developed for their preparation. Among those, the most used one relies on the reduction of the desired oxime (so-called hydroxyimine) derivative. This synthetic route has also been used for the preparation of their phosphono analogs, namely (α)-hydroxyaminophosphonate through reduction of the required (α)-hydroxyiminophosphonate.30 However, recent studies have reported that the yield of this pathway can be strongly influenced by the geometrical isomerism of the oxime moieties (E/Z isomerism, Fig. 3).
Fig. 3 Representation of the geometrical isomerism (E/Z) of the oxime moiety for diethyl (1-(hydroxyimino)-2-phenylethyl)phosphonate 2b. |
For instance, Cramer et al. recently showed that the metal-assisted hydrogenation of the Z isomer of an oxime derivative was much faster than for its E counterpart.31 While studying the reduction of the oxime isomers of (S)-strepetenol A, Blechert et al. also observed a clear difference in both reactivity and selectivity.32 However, in this case, they highlighted that the E-oxime was reduced faster and under milder conditions that the Z-oxime. Nonetheless, both studies showed that under appropriate conditions and reaction time, both isomers reacted to give the desired hydroxylaminated derivative. Very few studies were conducted on the (α)-hydroxyiminophosphonates (i.e., phosphono analogs of (α)-hydroxyimino derivatives) reduction towards (α)-hydroxyamino derivatives. In 2012, Beletskaya et al. showed that under regular metal-assisted hydrogenation, only one of the two oxime isomers underwent reduction, regardless of the nature of the ligand, temperature, hydrogen pressure or reaction time. The nature of the two isomers was determined by 31P NMR due to their difference in chemical shifts. The isomer which remained untouched was attributed to the E-isomer while the Z-isomer completely reacted towards the desired product.30
From these observations, it seems of paramount importance to have a comprehensive understanding of the oxime reduction mechanism involving isomerism for improving the synthesis of hydroxyamino derivatives. This is only reachable with a complete characterization of the oxime isomers. In the 1990s and early 2000s, Eli Breuer and co-workers reported in a series of research papers, as well as in a review, on the preparation, characterization and study of (α)-hydroxyiminophosphonate derivatives.33–40 In their study, they showed that NMR was an invaluable tool for the assessment of the two isomers. In fact, in 31P NMR, the E-isomer resonates at lower field than the Z-isomer, which is backed by the expected shielding effect of the phosphorus by the lone pairs of the oxime oxygen.35 Also consistent with this shielding effect is the upfield shift observed in the 13C NMR spectra of the signal corresponding to the oxime carbon of the E-isomer relative to the Z-isomer. Moreover, a noticeable difference is observed in the coupling constant 1JPC of the oxime carbon with the phosphorus that is higher than 200 Hz for the E-isomer, while for the Z-isomer it usually ranges from 150 to 160 Hz.41 In the case of crystalline compounds, X-ray diffraction (XRD) was proven to be a valuable method for the assessment of the two isomers. During their study, they also discovered a peculiar behavior difference in the degradation of those species. Under thermal stress, the Z-isomer undergoes fragmentation into the corresponding nitrile and dimethyl hydrogen phosphonate through a 4-membered ring intermediate, while the E-isomer did not react. Though this approach proved to be useful in the study of the (α)-hydroxyiminophosphonate isomers and their properties, the main drawback is that it requires isolation of large amounts of the pure isomers, which can be extremely difficult or impossible in some instances.
Besides NMR and XRD, mass spectrometry (MS) is a characterization method widely used to determine the elemental composition and/or partial structural insight using tandem mass spectrometry (MS/MS) such as the atomic connectivity.42 MS is very useful for the analysis of low abundance analytes in complex mixtures, which can be challenging in NMR and XRD. MS can also be hyphenated with other separation methods, such as liquid chromatography (LC-MS), gas chromatography (GC-MS), capillary electrophoresis (CE-MS) and ion mobility (IM-MS). Based on the separation technique involved with MS analyser, isomeric compound can be isolated before their analysis by MS and MS/MS. Moreover, these separation techniques give access to other compound descriptors (e.g., retention time and migration time), increasing the characterization level. For instance, IM-MS has emerged as a well-suited technique to probe the shape of ions in gas-phase based on their mobility or collision cross sections (CCS).43 Separation and characterization of geometrical isomers from IM-MS analysis have been recently reported.44–47
Within the context of synthesizing (α)-aminophosphonic acid analogs of some of the 20 proteogenic amino acids through the reduction of (α)-hydroxyiminophosphonates, a systematic and thorough study of their peculiar behavior was envisioned. Herein is reported the development of a multifaceted approach to characterize (α)-hydroxyiminophosphonates involving the interplay of NMR, XRD, MS, IM-MS and computational chemistry techniques. This study encompasses 16 (α)-hydroxyiminophosphonates synthetized according to a revisited version of Viveros-Ceballos’ classical route for the preparation of (α)-aminophosphonates.19 The reported multifaceted approach leads to a high level of characterization, pushing forward the understanding of the oxime isomers behavior on the reduction mechanism towards hydroxyamino analogs.
The first step involved the Michaelis–Arbuzov reaction between the starting phenylacetyl chloride and triethylphosphite to yield the product as a mixture of tautomers (keto–enol) as reported by Zantour49 and Tishler,50 respectively. In this case, the product was recovered in 91% yield as an 83:17 (keto/enol) mixture as determined by 31P NMR. The second key step relied on the formation of the oxime moiety using hydroxylamine hydrochloride and pyridine as the organic base. The desired (α)-hydroxyiminophosphonate 2b was thus obtained in 66% isolated yield as a mixture of two geometrical isomers (E:71, Z:29). The isomer attribution was conducted using Eli Breuer previous (31P and 13C) NMR results.35 Because 2b is crystalline as a mixture of isomers, XRD analysis was also carried out on the mixture, which reinforced the result from the NMR attribution with a 67:33 E/Z ratio (see ESI section 2.6†). Next, the preparation of a small library of oxime derivatives was envisioned by varying the starting acyl chloride to mimic the side chains of the natural amino acids as well as the phosphite source. The results are displayed in Fig. 4 below. In the case of 15b and 16b, the starting acyl chloride were not commercially available and were prepared from their commercially available carboxylic acid derivatives (see ESI section 2.1.1†).
Fig. 4 Scope of (α)-hydroxyiminophosphonates. The isolated yield and ratio of isomers E/Z (in parentheses) are indicated. All compounds were obtained according to the procedure described in Scheme 1. |
From these results, the effect of the phosphite source, and more specifically of the steric hindrance of the alkyl phosphite, can be highlighted on the ratio of isomers. In fact, comparing the results obtained for 2b and its phosphite derivatives, the more steric hindering triisopropryl phosphite 3b leads to a higher content of E isomer (73.5%) in the mixture compared to triethyl 2b (52.6%) and trimethyl phosphite 1b (45.7%), respectively. The same trend applies for the phosphite derivatives of 4b, 5b, 6b and 7b, 8b, 9b. Looking at the remaining compounds in the scope, which were prepared using triethyl phosphite, it appeared that the nature of the side chain also has an effect on the ratio of isomers. In most cases, the E isomer is predominant, showing that it is more favored compared its Z counterpart. In the case of 10b we determined a 98:2 (E/Z) ratio, which can be attributed to the lack of steric hindrance on the side chain, resulting in almost complete formation of the sole E isomer while for 14b the 94:6 (E/Z) ratio was determined after column purification, which resulted in loss of some of the Z isomer. During this study, HPLC/DAD also emerged as a potent method for the separation and quantification of the (α)-hydroxyiminophosphonate isomers. Indeed, the ratio of isomers obtained by HPLC/DAD matches the one obtained via31P NMR and XRD with a 1–4% margin of error. Moreover, the elution order of the oxime isomers was kept identical, regardless of the structure: for the library of 16 compounds, the Z isomer has a 0.5 min lower retention time than its E counterpart (see ESI Fig. S1†).
To push the characterization of the (α)-hydroxyiminophosphonate isomers one step further, we envisioned using complementary methods such as mass spectrometry. However, regular mass spectrometry gives us the mass-to-charge ratio (m/z) which, in this case, is not sufficient to be able to characterize geometrical isomers. Therefore, the use of tandem MS (MS/MS) was required in order to access partial structural insights such as the atomic connectivity through fragmentation. Although such insights were within reach of MS/MS analysis, isolation of the isomers was still mandatory in order to characterize them. Hence, a Synapt G2 spectrometer coupled with ion-mobility (IM) was used as a mean of separation prior to MS/MS analysis. Even though the separation of the isomers was poorly resolved, this enabled the acquisition of isomeric interference-free MS/MS spectra for the isomers of 2b. The results showed that the CID-mediated fragmentation exhibits different patterns for the two isomers studied (Fig. 5). In fact, the Z isomer was concerned with a first fragmentation occurring at the C–P bond, releasing the phosphorus copula while for the E isomer, the first fragmentation occurred through a McLafferty rearrangement resulting in the loss of one of the ethyl esters. For the remaining compounds of the scope, the IMS resolution of the instrument used was not sufficient to acquire isomeric interference-free spectra. Therefore, HPLC was used for the isomer separation prior to MS/MS analysis. Similar fragmentation pathways were observed on most of the 16 compounds studied. Predominant fragments for the isomers of each compound are reported in ESI Table S2.†
Fig. 5 Isomeric interference-free MS/MS spectra for the isomers of 2b highlighting the difference in their CID-mediated fragmentation pattern. |
From these observations we can highlight that the two oxime isomers have a different behavior on a multitude of facets, including HPLC, NMR, MS/MS and towards their reduction into hydroxylamine derivatives. Therefore, there is a need to better comprehend their behavior, which could extend to a better understanding of their overall reactivity.
Upon analysis of the main fragments and their order of appearance for the scope of isomers, we observed a few trends. The first one relies on the rupture of the C–P bond (between the carbon of the oxime and the phosphonate moiety), while the second and third ones rely on the impact of the isopropyl esters and the ethyl esters on the fragmentation pattern of the isomers, respectively. Regarding the first trend, for the studied compounds, the fragmentation of the C–P bond only occurs for the Z isomers, while it remains untouched for the E isomers except for the methyl ester derivatives 1b, 4b and 7b, where the fragmentation of the C–P bond occurs for both isomers. However, looking at the curve of appearance of the ions resulting from the C–P fragmentation for the isomers of 1b (Fig. 7) (see ESI section 3.1† for the detailed breakdown curves), clearly emphasizes that fragmentation is predominant for the Z isomer. These observations highlight that the strength of the C–P bond in the Z isomer is significantly weaker than in its E counterpart. The second trend concerns the isopropyl ester derivatives 3b, 6b and 9b, which exhibit the same first two fragmentations for both isomers with the successive loss of the two isopropyl esters. This suggests that for these compounds, the P–O isopropyl esters bonds are weaker than the C–P bond. The last trend relies on the ethyl ester derivatives. In fact, we observed that all the Z isomers of the studied ethyl ester derivatives had the same fragmentation pattern regardless of the side chain. Firstly, fragmentation of the C–P bond occurs, releasing diethyl hydrogen phosphate ((HO)P(O)(OEt)2) which then undergoes successive loss of the ethyl esters through McLafferty rearrangement to yield phosphoric acid. Whereas, for the fragmentation of the E isomers, only the successive loss of the two ethyl esters through rearrangement is observed. These observations expand to the whole scope of ethyl ester derivatives except for 14b, 15b and 16b, which have a different behavior compared to the rest of the scope. In fact, these three compounds exhibit a similar fragmentation pattern for both isomers, which can be explained by the presence of heteroatoms along the side chain, potentially influencing the nature and strength of the weakest bond.
Moreover, thanks to the breakdown curve analyses, the monitoring of the parent precursor ion with respect to the collision energy (V) for the two isomers and the mixture of isomers can be used for measuring a valuable information on the strength of the weakest bond, namely, V50. In fact, the V50 corresponds to the collision energy (V), for which we observe 50% fragmentation of the parent precursor. These values highlight the difference in the strength of the weakest bond between isomers. Since CID is an ergodic process, it is important to note that the V50 values estimated from CID are directly related to the number of vibration modes. Therefore, only geometrical isomers (i.e. having the same number of vibration modes) can be compared in terms of V50. We therefore calculated ΔV50 to compare compound behavior upon fragmentation.
The V50 values for each isomers of the scope as well as the mixture of isomers are gathered in Table 1.
Compound | V 50 (E)a (V) | V 50 (Z)a (V) | V 50 (mix)a (V) | ΔV50 (V) |
---|---|---|---|---|
a The V50 values have been extracted from the knee point of the sigmoid equation obtained by fitting the decreasing signal of precursor ion. Random sample consensus (RANSAC) algorithm has been applied to remove potential outliers. b Isolation of the Z isomer via HPLC failed due to its too low content in the mixture of isomers. | ||||
1b | 17.68 | 10.77 | 16.78 | 6.91 |
2b | 15.40 | 12.41 | 15.08 | 2.99 |
3b | 9.56 | 12.20 | 9.59 | −2.64 |
4b | 18.90 | 10.09 | 18.67 | 8.81 |
5b | 15.07 | 11.33 | 15.07 | 3.74 |
6b | 9.36 | 10.88 | 9.41 | −1.52 |
7b | 20.42 | 10.36 | 20.27 | 10.06 |
8b | 15.53 | 12.25 | 15.55 | 3.28 |
9b | 9.48 | 11.13 | 9.48 | −1.65 |
10b | 12.76 | — | — | — |
11b | 16.48 | 10.79 | 16.08 | 5.69 |
12b | 14.67 | 10.67 | 13.81 | 4.00 |
13b | 16.56 | 12.06 | 16.41 | 4.50 |
14b | 8.72 | 8.99 | 8.74 | −0.27 |
15b | 11.56 | 10.13 | 11.42 | 1.43 |
16b | 13.57 | 11.70 | 13.17 | 1.87 |
Upon analysis of the results obtained from the V50 values for the various compounds of the scope, we observe that the nature of the parent alkyl phosphite has an impact on the ΔV50. In fact, if we look at the trimethyl phosphite (P(OMe)3) derivative 1b, a higher ΔV50 is obtained by comparison to the triethyl phosphite (P(OEt)3) derivative 2b, translating in a higher difference in the strength of the weakest bond between the two isomers. Whereas for the triisopropyl phosphite (P(OiPr)3) derivative 3b, a negative ΔV50 is observed. This trend also applies for the derivatives 4b, 5b, 6b and 7b, 8b, 9b. In the cases of 14b, 15b and 16b, which showed similar fragmentation patterns for the two isomers unlike the other triethyl phosphite derivatives, a smaller ΔV50 than for the other P(OEt)3 derivatives is observed, probably due to the presence of heteroatoms along the side chain which drives the fragmentation independently to the isomer form (i.e. same breakdown curves for both isomers). Moreover, the survival yield curves for the mixture of isomers highlight a peculiar phenomenon. In fact, for each compound of the scope, the curve for the mixture of isomers is almost a perfect match of the curve for the E isomer whereas it should consider the weighted contribution of each isomers (i.e. sum of the two sigmoid curves) (see ESI section 3†). This observation suggests that the E/Z ratio is influenced by the MS analysis.
Isomer specific in-source fragmentation and/or difference in protonic affinity likely weighs a lot in the appearance of such phenomenon. The former hypothesis is easily demonstrated by looking at the full MS spectra without collision activation for both isomers. In fact, in the case of 2b, a higher content of fragmentation for the full MS spectrum of the Z isomer compared to the E isomer (Fig. 8) is observed, emphasizing that the Z isomer suffers more of in-source fragmentation than the E isomer.
Fig. 8 Isomeric interference-free full MS spectra without collision activation for the isomers of 2b highlighting a higher content of in-source fragmentation for the Z isomer. |
In order to verify the latter hypothesis, favorable protonation site calculations were conducted on 9 couple of oxime isomers (1b–9b) by DFT at the B3LYP/6-31G(d) level of theory. Calculations showed that each isomer exhibits two stable protonation sites: one in the amine moiety (–N–H+) and the other on the phosphonate moiety (–PO–H+). Standard enthalpy, entropy and free enthalpy variations at 298.15 K (, , , respectively) of the proton transfer reaction (–N–H+ ⇌ –PO–H+) are reported in Table 2. A negative value of indicates that –PO–H+ is more stable than –N–H+. Inversely, a positive value of indicates that –N–H+ is more stable than –PO–H+. The calculated values support that the phosphonate moiety is the more favorable protonation site for E isomers while protonation occurs preferentially on the amino moiety for Z isomers regardless of the phosphite source. Since the favorable protonation site are isomer-dependent, one can consider that the protonic affinity is also dependent on the isomeric form. Based on the DFT calculations, this difference in the protonation sites may explain the difference observed for the breakdown curves of the two isomers as well as their difference in V50, except for compounds 14b, 15b and 16b which display similar behavior and V50 values for both isomers. This might be explained by the presence of heteroatoms along the side chain, potentially changing the nature of the protonation site/protonation behavior.
Compound | Isomer | –N–H+ ⇌ –PO–H+ | ||
---|---|---|---|---|
1b | E | −10.21 | −7.78 | −7.89 |
Z | 70.58 | 11.78 | 67.07 | |
2b | E | −14.27 | −5.07 | −12.75 |
Z | 53.71 | −11.94 | 57.27 | |
3b | E | −18.58 | 6.11 | −20.40 |
Z | 8.02 | −5.79 | 9.74 | |
4b | E | −20.86 | −14.69 | −16.48 |
Z | 6.95 | 4.84 | 5.51 | |
5b | E | −14.44 | −7.35 | −12.25 |
Z | 71.98 | 22.29 | 65.34 | |
6b | E | −24.96 | −7.70 | −22.67 |
Z | 69.26 | 47.58 | 55.07 | |
7b | E | −17.37 | −2.71 | −16.56 |
Z | 81.18 | 23.05 | 74.30 | |
8b | E | −22.61 | 4.63 | −23.99 |
Z | 74.88 | 21.52 | 68.47 | |
9b | E | −36.08 | −8.05 | −33.68 |
Z | 64.52 | 13.59 | 60.46 |
Although MS methods were unable to successfully quantify the isomeric ratio of the investigated oximes (due to overestimation of the E/Z ratio), numerous valuable structural information can be obtained. Indeed, tandem MS (MS/MS) gives access to structural insights such as the atomic connectivity through fragmentation. However, in some cases where the fragmentation pattern is the same for the two isomers, atomic connectivity is not sufficient on its own to characterize these species. Therefore, collision energy breakdown curves can be used as a powerful addition, highlighting the order of apparition of the fragments as well as their related energy of fragmentation. This allows the determination of isomer specific bond strength.
Footnotes |
† Electronic supplementary information (ESI) available: experimental procedures, off-line analysis with structural assignments (NMR, X-ray diffraction and Mass Spectrometry, breakdown curves) and computations. CCDC 2115990–2115991. For ESI and crystallographic data in CIF or other electronic format see DOI: 10.1039/d1qo01564h |
‡ Equal contributions. |
This journal is © the Partner Organisations 2022 |