Boris
Gouilleux
*a,
François-Marie
Moussallieh
a and
Philippe
Lesot
*ab
aUniversité Paris-Saclay, UFR d'Orsay, RMN en Milieu Orienté, ICMMO, UMR CNRS 8182, Bât. HM1, 17-19, Avenue des Sciences, 91400, Orsay, France
bCentre National de la Recherche Scientifique (CNRS), 3, Rue Michel Ange, F-75016 Paris, France. E-mail: boris.gouilleux@universite-paris-saclay.fr; philippe.lesot@universite-paris-saclay.fr
First published on 17th April 2024
Controlling the enantiomeric purity of chiral drugs is of paramount importance in pharmaceutical chemistry. Isotropic 1H NMR spectroscopy involving chiral agents is a widely used method for discriminating enantiomers and quantifying their relative proportions. However, the relatively weak spectral separation of enantiomers (1H Δδiso(R, S)) in frequency units at low and moderate magnetic fields, as well as the lack of versatility of a majority of those agents with respect to different chemical functions, may limit the general use of this approach. In this article, we investigate the analytical potential of 19F NMR in anisotropic chiral media for the enantiomeric analysis of fluorinated active pharmaceutical ingredients (API) via two residual anisotropic NMR interactions: the chemical shift anisotropy (19F-RCSA) and dipolar coupling ((19F–19F)-RDC). Lyotropic chiral liquid crystals (CLC) based on poly-γ-benzyl-L-glutamate (PBLG) show an interesting versatility and adaptability to enantiodiscrimination as illustrated for two chiral drugs, Flurbiprofen® (FLU) and Efavirenz® (EFA), which have very different chemical functions. The approach has been tested on a routine 300 MHz NMR spectrometer equipped with a standard probe (5 mm BBFO probe) in a high-throughput context (i.e., ≈10 s of NMR experiments) while the performance for enantiomeric excess (ee) measurement is evaluated in terms of trueness and precision. The limits of detection (LOD) determined were 0.17 and 0.16 μmol ml−1 for FLU and EFA, respectively, allow working in dilute conditions even with such a short experimental duration. The enantiodiscrimination capabilities are also discussed with respect to experimental features such as CLC composition and temperature.
Several approaches have been developed for enantiomeric discrimination and the determination of enantiomeric excess (ee), so far, such as chiroptical tools (VCD, ROA, …), gas or liquid chiral chromatography (GC, HPLC), electrophoresis, mass spectrometry and NMR spectroscopy.5 Although chromatographic tools are widely used for this purpose, NMR spectroscopy can also be successfully used for the enantiomeric analysis of chiral API.6–13 These NMR methods rely on the interaction between the analyte to be studied and a chiral auxiliary (CA) to form: (i) either stable diastereoisomers by covalent bonding or (ii) diastereomeric adducts via intermolecular interaction.14–16 In isotropic solutions, these CAs can be either chiral derivatizing agents (CDAs) as Mosher's acid,17 chiral metallic complexes,18,19 chiral lanthanide shift reagents,20,21 or chiral solvating agents (CSAs) such as cyclodextrin-based agents.22 The conversion of enantiomers to diastereomeric entities leads to spectral discriminations through variations of isotropic chemical shift (Δδiso(R, S)). Basically, weak spectral enantio-separations (in frequency units) generally measured on routine 1H NMR (<9.4 T) can be increased by recording spectra on very high-field NMR instruments (if accessible). Other approaches using 19F NMR on fluorinated solutes23 or even fluorinated CSA/CDA as recently reported,24–27 can provide elegant tools. However, the need of designing a tailored CSA/CDA (not always commercially available) for a specific chemical group remains the main challenge and may limit the general use of these approaches.
The concept of enantiospecific non-covalent interactions between the R/S-isomers and a chiral (enantiopure) selector in isotropic solvents can be extended to anisotropic media, such as chiral weakly aligning systems (e.g., helical polymer-based lyotropic liquid crystals). In this case, enantiomers exhibit a difference of orientational order in average, these molecular orientations being generally expressed by the Saupe tensor Sαβ.28 From NMR point of view, the spectral enantiomeric discriminations are observed through a difference of residual anisotropic NMR interactions, such as residual chemical shift anisotropy (RCSA), residual dipolar coupling (RDC) or residual quadrupolar coupling (RQC), whose amplitudes depend on the orientational order of the analyte.29–34
Compared to NMR methods in isotropic solvents supported by CAs, NMR in chiral weakly aligning media offers at least three analytical advantages that can expand the role of NMR spectroscopy in chiral analysis of chiral drugs. The first and the most interesting one is the diversity and the magnitude of these anisotropic observables (RDC, RQC, RCSA) associated with any magnetically active nucleus (whether abundant or not), and that generally allows for effective spectral enantiodiscriminations. In addition, the magnitude of RQCs and RDCs (unlike a chemical shift difference) is independent of the magnetic field strength (B0), making easier to transfer methods to spectrometers operating at lower fields. The second one results from the shape recognition mechanisms that play a central role in enantiodiscrimination mechanisms in CLCs.35 As a result, no specific chemical functions or group on the analyte (unlike CSAs/CDAs) is required, so that all CLCs are potentially able to differentiate enantiomers whatever the chemical characteristics of the chiral analyte. The third advantage is based on the absence of possible differences in affinity between the analyte and the chiral agent which can lead to significant discrepancies between the real value and the measured value of the ee%.24,27De facto, no correction coefficient must be introduced to determine the ee% from the measurement of peak surfaces for each enantiomer. The last but not the least advantage is that several chiral polymers (as polypeptides) are commercially available, and no preliminary synthetical step is required to explore their potential.
Although these properties have benefited organic chemistry for the analysis of small molecules, NMR in CLCs has been little explored for the study of pharmaceutical compounds, so far. The main reason for this lack of interest is the inherent complexity of anisotropy 1H 1D NMR spectra of medium-sized molecules due to the numerous long-range (1H–1H)-RDCs that broaden spectral patterns. Other exploitable heteronuclear experiments (e.g., 13C and 2H) are possible but may suffer from insufficient sensitivity to deal with samples of 5–20 mg of API at high enantiomeric excess in a high-throughput framework with routine NMR equipment (e.g., 300–400 MHz spectrometers and conventional probes). In this context, 19F–{1H} NMR appears to be a promising alternative since in the case of monofluorinated or trifluoromethylated chiral compounds (F-R* or CF3-R*) encountered in pharmacopeia,36 a dramatic simplification in spectral analysis is expected, with a suitable/reasonable sensitivity (≈86% of 1H NMR) and a significant chemical shift dispersion of 19F resonances. This advantage has recently been exploited for the chiral analysis in isotropic solvents involving either fluorinated analytes (or even a mixture of chiral analytes),7,23 or newly developed chiral agents including a fluorine tag.24–26
From anisotropic NMR point of view, 19F nuclei are rather sensitive nuclear spies to observe molecular orientational order differences induced by weakly chiral oriented media. These differences can be observed through two residual anisotropic interactions for R and S isomers: the 19F-RCSA, also noted, 19F σaniso(R, S), and the (19F–19F)-RDC, noted |DFF(R, S)|. Based on 19F-RCSA, the spectral enantiomeric discrimination is measured through a difference of frequency (in Hz or ppm) between each enantiomer:
|Δνaniso(R, S)| = |νaniso(R) − νaniso(S)| | (1) |
![]() | (2) |
Expressed in a molecular frame (a, b, c axes), the terms σiso and σaniso can be derived from the shielding tensor as:
![]() | (3) |
In the above equations ν and σ are related to the 19F nuclei while γF is the fluorine gyromagnetic ratio.
Based on (19F–19F)-RDC, the spectral discrimination between enantiomers is measured as:
|ΔDFF(R, S)| = |DFF(R) − DFF(S)| | (4) |
![]() | (5) |
The order parameter SFF of the internuclear vector rFF between the coupled 19F nuclei can be derived as:
![]() | (6) |
Previous analysis of around twenty small fluorinated chiral compounds (monofluorinated or including a –CF3 moiety) using poly-γ-benzyl-L-glutamate (PBLG) as a chiral mesophase has shown that variations in the 19F anisotropic chemical shift and 19F–19F total homonuclear coupling lie in a range from 0.01 to 0.350 ppm and from 1.4 to 118 Hz, respectively.37–39 Based on this valuable enantiodiscrimination ability, we investigate for the first time the analytical potential of this approach to determine enantiomeric excesses (ee) for fluorinated chiral pharmaceutical ingredients using a routine 300 MHz (7.05 T) spectrometer equipped with a conventional 5 mm NMR probe, a commonly encountered setup in both academic and industrial laboratories. The key role of the mesophase composition and temperature on the enantiomeric resolution is investigated and presented in detail. The analytical performance of the herein method is then evaluated in terms of precision (here short-term repeatability), trueness (i.e., deviation from the expected ee% values) and limits of detection. Two different APIs are targeted in this study: (i) Flurbiprofen® (FLU), a chiral non-steroidal anti-inflammatory drug where the (S)-form possesses most of the beneficial anti-inflammatory activity,40 and (ii) Efavirenz® (EFA), a chiral drug whose (S)-isomer is a non-nucleoside inhibitor of human immunodeficiency virus (HIV) reverse transcriptase (see Scheme 1).41 These two bioactive chiral molecules involve various chemical functions: FLU is a biphenyl propionic carboxylic acid while EFA is mainly characterized by a carbamate function into a benzoxazine ring, as well as a difference of fluorine moiety: a single fluorine in FLUversus a trifluoromethyl group (–CF3) in EFA. In addition, the fluorinated group used for the enantiomeric discrimination is directly bound to the stereogenic center in EFA whereas the latter is located further away in the case of FLU. Interestingly, while EFA undergoes restricted degrees of freedom, FLU is a more flexible molecule which can experience conformational changes upon interaction with PBLG side chains.42 This may lead to different chiral recognition mechanisms between these two APIs with respect to the CLC. From an analytical point of view, all these chemical characteristic differences make FLU and EFA very well-suited analytes to test the versatility of PBLG mesophase when applied to the analysis of pharmaceutical compounds.
The oriented samples were prepared directly in the NMR tube by mixing the chiral analyte (<30 mg), the PBLG (nearby 100 mg) and the organic co-solvent (CHCl3). The amount of chloroform (about 600 mg) was here adjusted so that the total weight% of PBLG (i.e., mass of PBLG/total mass) was maintained at around 14%, and the sample length reached 4 cm. The exact sample compositions are given in Table SI-1 in ESI.† The 5 mm NMR tube was then fire-sealed to avoid CHCl3 evaporation over time and several low-speed centrifugation cycles of the sample (e.g., 500 rpm during 20 s) were carried out to remove matter gradients. The tube was inverted between each centrifugation.
![]() | (7) |
The former situation is observed on a racemic sample of FLU (see Fig. 2a) where the signals of enantiomers are separated by a difference of chemical shift anisotropy: 19F |Δδaniso(R, S)| = 0.051 ppm, while the (R)- and (S)-isomers of a racemic sample of EFA are discriminated by variations of 19F-RCSA: 19F |Δδaniso(R, S)| = 0.085 ppm and (19F–19F)-RDC: |ΔTFF(R, S)| = 12.1 Hz as shown in Fig. 2b. The linewidth is of the order of 2–3 Hz at half-maximum in these spectra, which underlines the satisfactory resolution obtained in such lyotropic CLCs. Note the weak lines marked by black stars in Fig. 2a do not arise from fluorinated impurities, but rather to satellite lines owing to short and long-range (13C–19F)-RDCs. Assignment of these signals to the multiple 13C-isotopomers is detailed in ESI-III.†
![]() | ||
Fig. 2 282.4 MHz 19F–{1H} NMR spectra on 20 mg of racemic FLU (a) and EFA (b) dissolved in a PBLG/CHCl3 based lyotropic liquid crystal. The assignment of the stereodescriptors R and S relies on further experiments performed on scalemic samples. The black stars at the bottom of NMR lines correspond to 13C–19F satellite lines (see details in ESI-III†). |
The potential role of the analyte concentration on the resulting values of 19F-RCSA and (19F–19F)-RDC has been little investigated so far, since it is generally considered as a guest molecule with minor impacts on the orienting properties of the mesophase. To explore this effect, anisotropic 19F–{1H} NMR experiments are carried out on a series of samples with a growing amount of (S)-EFA from 5.8 mg (0.046 mol L−1) to 30.3 mg (0.24 mol L−1), while the mass percentage of PBLG is kept constant nearby 14% (see Table SI-1†). The spectra of Fig. 3 show a large increase of both 19F-RSCA and (19F–19F)-RDC as the concentration in (S)-EFA decreases. Besides, the offset from the 19F δiso value varies from 0.080 to 0.560 ppm (i.e., +600%), while |TFF| increases from 13.9 to 131.7 Hz (i.e., +847%). Such a large variation in RDC quantities is not calculated for the co-solvent CHCl3, whose the associated total coupling |TCH| only varies from 309 to 315 Hz (i.e., 1.9% variation) for samples from 0.16 to 0.045 mol L−1 in EFA, which rules out a destabilization of the mesophase with the increasing amount of solute. In addition, analysis of 1H self-diffusion experiments (1H DOSY NMR)44,45 under isotropic conditions has also precluded any possibility of concentration-dependent self-aggregation leading to a conformational change in EFA. Details are given in ESI-VII.†
Two arguments can be put forward to explain the variation in (19F–19F)-RDC in methyl group by simply considering the order parameter associated with the C–CF3 vector:
![]() | (8) |
(i) A change of SC–CF3 induced by a reorientation of the C–CF3 axis of the solute with a subsequent modification of the angle between the C–CF3 vector and the magnetic field B0. This can lead to a significant increase of RDC, especially if
is initially closed to the magic angle (i.e., 54.7°), (ii) a modification of the ensemble average 〈…〉 as a function of the solute concentration. This ensemble average can be seen, in a first approximation, as a fast exchanging system where:
![]() | (9) |
In a second step, it is interesting to investigate the impact of the solute concentration when the two enantiomers are present in the chiral liquid crystal. At first, 19F–{1H} NMR experiments are performed on a series of racemic samples prepared with various amounts of FLU, from 4.8 to 30.0 mg, corresponding to a 0.025–0.153 mol L−1 concentration range for each enantiomer. Table 1 shows clearly the effect of solute concentration on the enantiomeric discrimination with an increase of 19F |Δδaniso(R, S)| from 0.040 to 0.072 ppm for the highest and lowest concentration, respectively. In this case, diluted samples lead to a better spectral separation, which is an interesting analytical feature as the sensitivity of 19F NMR enables experiments at relatively low concentrations. Furthermore, 19F–{1H} spectra are recorded on three different scalemic samples (ee(R) = 9.6, 58.2 and 88.7%) prepared with a fixed total amount of FLU (i.e., 30 mg). As seen in Table 1, the variation of 19F-RCSA progressively increases as a function of the enantiomeric excess: 19F |Δδaniso(R, S)| = 0.054 ppm for ee(R) = 9.6% and 0.082 ppm for ee(R) = 89.5%. This effect of the enantiomeric excess on the resulting 19F–{1H} spectra is also visible for EFA as pointed out in Fig. SI-3† where in one hand 19F |Δδaniso(R, S)| = 0.095 ppm and |ΔTFF(R, S)| = 7.8 Hz for ee(S) = 71.5%, while 19F |Δδaniso(R, S)| = 0.109 ppm and |ΔTFF(R, S)| = 4.8 Hz for ee(S) = 90.5% in another hand.
m(FLU)a (mg) | ee(R) (%) | 19F |Δδaniso(R, S)| (ppm) |
---|---|---|
a Total mass of FLU: m(R-FLU) + m(S-FLU). b All experiments were carried out with 14% (wt) of PBLG and at 297 K. | ||
30.0 | 0.0 | 0.0403 |
20.2 | 0.0 | 0.0507 |
14.8 | 0.0 | 0.0549 |
10.2 | 0.0 | 0.0706 |
4.8 | 0.0 | 0.0720 |
30.0 | 0.0 | 0.0403 |
29.9 | 9.6 | 0.0543 |
30.8 | 58.2 | 0.0670 |
30.4 | 88.7 | 0.0820 |
![]() | ||
Fig. 4 Effect of the temperature on the enantiomeric discrimination by variation of 19F-RCSA evaluated on a 30 mg racemic sample of FLU dissolved in PBLG/CHCl3. |
A more intriguing effect is observed in the case of scalemic samples of EFA (e.g. ee(S) = 90.5%) where 19F-RCSA and (19F–19F)-RDC evolve differently from one enantiomer to another with respect to the temperature. While |TFF(S)| (major enantiomer) increases linearly with this parameter, as observed in the racemic sample, the evolution of |TFF(R)| is more complex with a drop close to zero Hz at 310 K before increasing again at higher temperatures (see Fig. 6). Fundamentally, this non-linear variation suggests a modification of the average orientation of the (R)-isomer within the liquid-crystalline phase, leading to a change in enantiomeric discrimination. In addition to this enantiospecific variation of (19F–19F)-RDC, the 19F δaniso(S) and 19F δaniso(R) evolve linearly with the temperature, but at different rates resulting in 19F |Δδaniso(R, S)| = 0.116 ppm at 300 K versus 0.099 ppm at 320 K.
![]() | ||
Fig. 6 Effect of the temperature monitored on a 30 mg scalemic sample of EFA (ee(S) = 90.5%) dissolved in PBLG/CHCl3. In contrast to the racemic series, different evolutions of 19F-RCSA and (19F–19F)-RDC are observed from one enantiomer to another leading to variations of signal patterns. The variation of both anisotropic parameters is given in Table SI-3.† |
Overall, these results highlight the significant impact of solute concentration, enantiomeric ratio and temperature on the enantiomeric resolution of these 19F NMR experiments in lyotropic chiral liquid crystals.
Another possible bias may occur when 1H-decoupling in the heteronuclear experience is applied and leads to a differential nuclear Overhauser effect (nOe) between analytes. To limit this effect, it is recommended to apply the 1H decoupling only during the 19F signal detection period (“inverse gated” experiment). However, these two usual requirements may be alleviated in the case of relative quantification of enantiomers. Indeed, relaxation properties are similar between the two enantiomers and no significant differences of relaxation times have been so far reported in lyotropic liquid crystals, at the best of our knowledge. Taniso1(19F) have been determined for the R and S-isomers of both FLU (Taniso1(19F) = 1.2 s) and EFA (Taniso1(19F) = 0.5 s) by the inversion–recuperation method and no differences (greater than the fitting error) have ever been observed between enantiomers. This allows to work with a smaller repetition time (TR ≈ 1.25 × Taniso1(19F)) which provides the optimal SNR per time unit without a systematic quantitative bias. Interestingly, this condition therefore improves the method's detection limit for a given experimental duration.
19F–{1H} NMR experiments (NS = 8, i.e. Texp ≈ 10 s) were carried out on both racemic and scalemic series to evaluate the accuracy of enantiomeric excess measurements in different spectral situations. The experimental ee% are derived as in eqn (7) and the signal integration is performed by a line-fitting procedure (see Fig. 7) as described in the Experimental section. The deviation from the expected values of the measurements (i.e., the trueness) is assessed on racemic samples since the expected ee = 0% does not suffer from potential weighting uncertainties during the sample preparation. The precision (here defined as short term repeatability) is evaluated by the standard deviation (SD%) computed on series of at least seven successive NMR experiments. Scalemic samples are also considered to evaluate the performance of the method on more challenging cases, such as ee ≈ 90%, where large differences of signal area between lines imperfectly resolved are detrimental for an accurate ee% determination.
In the case of FLU (see Table 2), the trueness evaluated on the racemic sample is lower than 1% (0.8%), which validates the absence of significant biases from partial saturation (TR ≈ 1.25 × T1(19F)) whenever a trueness of 1% is targeted. The SNR achieved under these experimental conditions (Texp ≈ 10 s) is of the order of 103 and remains higher than 300 even for the minor enantiomer at ee = 88.7%. The resulting precision (i.e., SD%) is better than 0.5% for all the FLU samples (see Table 2). It could also be noted that precision is not affected by the ee% to be measured, which suggests that the uncertainty mainly relies on the line-fitting procedure, rather than a lack of SNR for the minor enantiomer. For scalemic samples, the trueness remains lower than 1% until ee(R) = 58.2%, but increasing significantly at ee(R) = 88.7% where a deviation of −4.2% is observed. This deviation can arise from an imperfect Lorentzian lineshape of experimental signals causing weak line-fitting errors, which become detrimental for the trueness at very high ee%, whenever enantiomeric signals are not perfectly baseline-separated. In the case of EFA (see Table 3), the trueness evaluated on the racemic sample is better than 0.5%, which demonstrates again the absence of biases of the proposed method. This good performance is also extended in the case of scalemic samples since the measured values of ee% lie in the acceptable interval with deviation lower than 0.5% from the expected values. Table 3 also gives the SNR values associated with the (R)- and (S)-isomers for a racemic and two scalemic samples of EFA. Such fast 19F NMR experiments enable a minimal SNR value of about 200 for the minor R-isomer in a sample with ee(S) = 90.5%. This leads to standard deviations lower than 0.5%, even at the highest enantiomeric excess (i.e., ee(S) = 90.5%) highlighting the high precision of this fast 19F NMR method.
ee Theo.(R)b (%) | ee Exp.(R)c (%) | Deviationd (%) | SDe (%) | SNR (S) | SNR (R) |
---|---|---|---|---|---|
a The total mass of API is nearby 30 mg in all samples. Measurements are achieved by 19F–{1H} NMR experiment of ≈10 s (NS = 8) at 7.05 T (ν0 = 282.4 MHz). b Values given with 95% confidence intervals due to weighing uncertainties for scalemic samples. c Mean value of seven successive experiments. d Computed as eeExp − eeTheo. e Standard deviation computed on seven successive experiments. | |||||
0.0 | 0.8 | 0.8 | 0.32 | 1854 | 1839 |
9.6 ± 0.7 | 9.5 | −0.1 | 0.40 | 1691 | 2130 |
58.2 ± 0.8 | 57.3 | −0.9 | 0.22 | 891 | 3019 |
88.7 ± 1.1 | 84.5 | −4.2 | 0.23 | 337 | 3465 |
ee Theo.(S)b (%) | ee Exp.(S)c (%) | Deviationd (%) | SDe (%) | SNR (S) | SNR (R) |
---|---|---|---|---|---|
The footnotes (a–e) are identical to footnotes of Table 2. | |||||
0.0 | −0.4 | −0.4 | 0.37 | 1762 | 1743 |
71.5 ± 0.9 | 71.7 | 0.2 | 0.50 | 2816 | 517 |
90.5 ± 1.1 | 90.9 | 0.4 | 0.35 | 3686 | 191 |
![]() | (10) |
![]() | ||
Fig. 8 Estimated LOD according to the magnetic field B0 of the spectrometer with fixed experimental durations Texp (a) and as a function of Texp for given B0 values (b). The LOD curves are extrapolated viaeqn (10) from the experimental LOD of EFA (black cross). |
Results in terms of spectral enantiodiscriminations are very promising with potential applications in pharmaceutical analytical laboratories. First, the determination of the ee% is obtained with a precision better than 0.5% and a trueness better than 1% (at the exception of one sample) in a high-throughput context (experimental durations ≤ 10 s) and routine NMR setup operating at moderate magnetic fields (e.g., 7.05 T, conventional NMR probe) with a limit of detection of 42 μg mL−1 (i.e., 0.17 μmol mL−1) for FLU, and 50 μg mL−1 (i.e., 0.16 μmol mL−1) for EFA, in such experimental conditions. Second, the use of anisotropic NMR interactions such as the 19F–19F residual dipolar couplings of methyl groups, whose amplitude is not related to the magnetic field strength, eases the method transfer at lower magnetic fields. This aspect could make enantiomeric analysis with benchtop NMR spectrometers possible in the near future.
This first study on API has also explored the significant role of the liquid crystal composition, especially the chiral analyte concentration on the resulting 19F-RCSA and (19F–19F)-RDC values, and showing higher values of 19F anisotropic observables on diluted samples as well as variations of the enantiomeric discrimination according to the ee%. The analysis of the effect of temperature indicates a slight improvement of the spectral discrimination of enantiomers of FLU at 297 K compared to higher temperatures whereas this parameter should be adjusted to avoid detrimental peak-overlaps in the case of EFA. The whole understanding of the complex interplay between mesophase composition and sample temperature on the average orientation of analyte and discrimination mechanisms, and thus on the measured residual anisotropic observables, is still an on-going challenge under investigation.
Footnote |
† Electronic supplementary information (ESI) available: (i) NMR sample composition, ee(%) and uncertainty; (ii) assignment of 13C–19F satellites on 19F–{1H} spectra of FLU in PBLG/CHCl3; (iii) further 19F–{1H} spectra of EFA at different ee% or temperature; (v) limits of detection (LOD). See DOI: https://doi.org/10.1039/d4an00237g |
This journal is © The Royal Society of Chemistry 2024 |