2.1. Variable Temperature and Concentration 1H NMR Chemical Shifts of Carboxylic Protons and 1D 1H NMR Transient NOE in CDCl3
The chemical shifts of the carboxylic protons, δ(COOH), and phenol OH group, δ(OH), are very informative criteria for the investigation of various types of hydrogen bond interactions [
28,
29,
30,
31]. δ(COOH) and δ(OH) are deshielded in the presence of hydrogen bond interactions and linear correlations between
1H NMR chemical shifts and hydrogen bond distances have been reported [
30,
31]. Temperature has also a significant effect, thus, by increasing the temperature, the
1H NMR chemical shifts are shielded due to breaking of hydrogen bond interactions (negative temperature coefficients, Δδ/ΔΤ). The
1H NMR resonances of the COOH groups display broad signals at room temperature in CDCl
3. The broadening is mainly due to intermolecular proton exchange of the COOH group with the residual H
2O in CDCl
3 solution. The use of low concentrations (c < 100 mM) has a profound effect on proton exchange rate, which results in excessive line broadening and variable chemical shifts. The use of activated molecular shifts in the bottom of the NMR tube, but outside the active volume of the NMR coil, resulted in a significant reduction in the line widths which allowed the accurate determination of the chemical shifts and Δδ/ΔΤ values.
δ(COOH) chemical shifts at 298 K, Δδ/ΔΤ (ppb K
-1), and statistical analysis (coefficient of linear regression R
2 and intercept) of the data of
Figure 1 are shown in
Table 1. The temperature-dependent changes of the chemical shifts of the free fatty acids investigated are linear and the derived Δδ/ΔT values, with R
2 > 0.992, cover a range of -42.74 to -29.52 ppb K
-1. These values are significantly larger, in absolute terms, than those obtained in the liquid state for caproleic acid, oleic acid, α-linolenic acid, EPA and DHA (-16.43 to -10.32 ppb K
-1) [
28] (
Table 1) and semi-fluorinated oleic, elaidic and stearic acids [
32]. This shows that, by increasing the temperature, the intermolecular hydrogen bonds are more readily broken in CDCl
3 solution than those in the liquid state.
Numerous investigations of various carboxylic acids in CCl
4 and CHCl
3 were interpreted in terms of mixtures of cyclic and linear dimers, cyclic and linear trimers and monomers [
33,
34,
35,
36,
37,
38,
39]. For long chain carboxylic acids, such as in FFAs, the formation of centro-symmetric hydrogen bond species through carboxylic groups appears to be the major structural mode. Thus, the single crystal X-ray structural analysis of linoleic acid, α-linolenic acid and arachidonic acid [
40] showed the formation of centro-symmetric cyclic hydrogen bonds, which deviate from planarity by 26.7
o, with short O
…O distances of 2.67 Å.
Figure 1 and the data of
Table 1 demonstrate that caproleic acid and oleic acid and the ω-3 ALA and EPA form intermolecular hydrogen bond interactions, since the chemical shifts of the carboxylic protons are strongly deshielded (11.17 to 10.39 ppm, at 298 K) (
Table 1). In caproleic acid, oleic acid, ALA and EPA the hydrogen bond species through carboxylic groups, therefore, are the major components in CDCl
3 solution. This is in agreement with literature data [
41] of the minor presence (1% to 3%) of the monomeric species in the liquid state for octanoic, nonanoic, decanoic and undecanoic acids in the temperature range of 280 K to 360 K.
The chemical shifts of the carboxylic groups of CA, OA, ALA, and EPA in CDCl
3 (
Table 1) are slightly more shielded by 1.17 to 0.14 ppm, relative to those in the liquid state [
28]. This can be attributed to the major role of the centro-symmetric cyclic dimers relative to contributions of other components of the equilibrium mixtures in both liquid state and CDCl
3 solution. Detailed dilution studies of caproleic acid in the range of 400 mM to 1 mM showed a very significant shielding in the concentration range below 15 mM due to increased contribution of the monomeric species. Thus, at 10 mM, the chemical shift of caproleic acid is ~ 8.6 ppm, while that of oleic acid, at 2 mM, is ~ 9.3 ppm. Further research in needed to determine the precise values of dimer-to-monomer dissociation constants, which apparently depend on the length of the side chain and the presence of multiple cis double bonds, as in the case of ω-3 fatty acids, which result in a significant ‘kink’ into the chain (see discussion below).
DHA is a particular case since the chemical shift of the carboxylic group is strongly shielded (δ = 9.07 ppm at 298 K). The chemical shift at 298 K is very similar to that in the liquid state (8.60 ppm) [
28] and the Δδ/ΔΤ value (-29.52 ppb K
-1) is larger to that in the liquid state. It can, therefore, be concluded that for the dimeric DHA in CDCl
3, a structural mode of intermolecular hydrogen bonds through carboxylic groups and an intermolecular hydrogen bond between the carboxylic group of one molecule and the terminal double bond of the second molecule of DHA, plays a significant role, as in the case of the liquid state [
28]. The OH
…π hydrogen bond has been suggested to have significant structural roles in bioorganic chemistry [
42,
43] and biochemistry [
44,
45].
1D transient NOE experiments were performed for caproleic acid (CA), oleic acid (OA), α-linolenic acid (ALA), EPA, and DHA using various concentrations (100 mM, 50 mM, and 20 mM) in CDCl
3 and various mixing times, τ
m. The NOE grows during the period τ
m, starting from zero [
33].
Figure 2 shows 1D NOE NMR spectra of oleic acid (OA) and α-linolenic acid (ALA) (concentration = 20 mM), using various τ
m values with selective excitation of the CH
3– group. Even for a short τ
m = 100 ms, there are weak NOE connectivities with the H2, H3 protons which are antiphase with respect to the irradiated CH
3– group. This is due to the formation of low molecular weight hydrogen-bonded species with τ
c values within the extreme narrowing condition (ω
οτ
c << 1) in the concentration range of 100 to 20 mM. By increasing τ
m, an approximately linear increase in the amplitude of the NOE signal intensities is observed which shows that the NOE is due to, through space, proximity of the CH
3– group and the C
H2–C
H2–COOH protons in the hydrogen bond species, rather than due to spin diffusion through the chain of the CH bonds.
Similar results were obtained with EPA (
Figure 3A). The magnitude of all the NOE signal intensities of DHA, however, is significantly reduced relative to those of OA, ALA and EPA. This can be attributed to the formation of low molecular weight hydrogen-bonded aggregates in the range of minimum NOE signal intensities, i.e., ω
οτ
c ~ 1.
1D transient NOE NMR spectra of the caproleic acid (CA), with selective excitation of the α-CH
2 protons, is shown in the
Supplementary Figure S1(A). As in the case of OA, ALA, EPA, and DHA the NOE connectivities are anti-phase with respect to the α-CH
2 group, due to the formation of low molecular weight hydrogen bond aggregates with τ
c values within the extreme narrowing condition (ω
οτ
c << 1). The magnitude of NOEs, however, with the terminal CH(9) = CH
2(10) protons was significantly less than those observed between α-CH
2 and the terminal CH
3- group of OA, ALA and EPA. This can be attributed to the minor formation of hydrogen bond interdigitated aggregates.
2.2. Variable Temperature 1H NMR Chemical Shifts of Carboxylic Protons and 1D 1H NMR Transient NOE in DMSO-d6
Exchange broadening due to intermolecular proton exchange between COOH groups and residual H
2O, can be significantly eliminated in DMSO-d
6 due to its strong hydrogen bond and solvation ability. δ(COOH) and Δδ/ΔT values can, therefore, be determined accurately. The chemical shifts of the carboxylic protons, δ(COOH), in DMSO-d
6 solution (c = 20 mM) are very similar and appear in a very narrow chemical shift range for all the FFAs (11.94-12.08 ppm) and are more deshielded relative to those in CDCl
3 (
Table 1 and
Figure 1). This shows that the centro-symmetric cyclic dimers do not exist in DMSO-d
6 due to the strong hydrogen bond and solvation ability of the DMSO molecules. In DHA, the flip-flop process between the classical intermolecular centro-symmetric bonds through the carboxylic groups and an intermolecular hydrogen bond between the carboxylic group of one molecule and the terminal double bond of the second molecule of DHA is also eliminated in DMSO solution. Further confirmation was also obtained from the Δδ/ΔT values in DMSO-d
6 (-6.62 to -7.72 ppb K
-1) which are significantly smaller, in absolute terms, than those in CDCl
3. This demonstrates that the effect of increasing the temperature results in significantly less pronounced breaking of hydrogen bond interactions in DMSO-d
6, relative to those in CDCl
3 solution.
The great hydrogen bond and solvation ability of DMSO is clearly demonstrated from variable temperature experiments of an equimolar mixture of caproleic acid and DMSO-d6. The chemical shift of the carboxylic proton at 298 K (δ = 11.90 ppm) and its temperature coefficient (Δδ/ΔΤ = -6.77 ppb K-1) clearly show the elimination of the centro-symmetric cyclic dimers through the carboxylic groups.
1D transient NOE experiments were performed for the FFAs in DMSO-d
6 with concentration c=20mM.
Figure 4 shows NOE NMR spectra of OA and ALA using various τ
m values with selective excitation of the terminal CH
3- group. Even for the relatively short τ
m = 100 ms, there are NOEs with the H2 and H3 protons which are antiphase with respect to the CH
3- group. This is due to the formation of low molecular weight hydrogen-bonded aggregates with τ
c values within the extreme narrowing condition (ω
oτ
m << 1). By increasing τ
m an increase in the amplitude of the NOE connectivities is observed which can be attributed to, through space, proximity of the CH
3–group and the C
H2–C
H2–COOH protons in the hydrogen bond species, rather than due to spin diffusion through the chain of the CH bonds.
Similar results were obtained with EPA and DHA (
Figure 5). Selective excitation of the terminal CH
3- group results in anti-phase NOE connectivities with H2, H3, even for the relatively short mixing time τ
m = 100 ms. This demonstrates the proximity, through space, of the CH
3- group and the C
H2–C
H2–COOH protons in the low molecular weight hydrogen bond interdigitated aggregates, within the extreme narrowing condition (ω
oτ
m << 1).
1D transient NOE NMR spectra of caproleic acid (CA), using various τ
m values with selective excitation of α-CH
2 protons, are shown in the
Supplementary Figure S1(B). The magnitude of the anti-phase NOEs, with the terminal CH(9)=CH
2(10) protons, was found to be significantly less than those observed between α-CH
2 and the terminal CH
3- groups of OA, ALA, EPA and DHA. This can be attributed to the minor formation of hydrogen bond interdigitated species.
2.3. DFT Calculations in CHCl3 – Comparison with the Liquid State
Computational approaches have been proved very successful in elucidating structural and spectroscopic experimental data of free fatty acids in the liquid state [
27,
28]. Moreover, this approach can be used as a predictive tool in biotechnology for predesigned properties of functional free fatty acid aggregates by tuning their interatomic interactions in organic materials [
46]. Based on the state of the FFA carboxylic proton, it can be determined if the FFA in the proper solution can be used as a transport or catalytic medium [
47]. The present computations were designed to investigate possible inter- and intramolecular interactions that justify the experimental δ(COOH) and 1D NOE NMR results, presented above. Caproleic acid was investigated in the dimeric structure forming O–H
…O=C centro-symmetric hydrogen bonds (
Figure 6a), in the cyclic trimeric (
Figure 6b) and linear trimeric (
Figure 6c) structures in implicit solvation (IEFPCM-chloroform). In the centro-symmetric dimeric structure (
Figure 6a), the dihedral angle defined by the four oxygen atoms of the carboxylic groups is only 0.8
o, the (O)H
…O(C) and O
…O hydrogen bond distances are 1.66 and 2.65 Å, respectively, and the O–H
…O bond angle is indicative of a nearly linear (178.0
o) hydrogen bond interaction. These values can be compared with the O
…O distance of 2.67 Å and deviation from planarity of 26.7
o in the single crystal X-ray structure of linolenic acid, α-linolenic acid and arachidonic acid [
40]. The experimental chemical shifts of caproleic acid (δ = 11.08 ppm at 298 K,
Table 1) are rather indistinguishable on the basis of the structures of
Figure 6a,b (13.6 ppm and 12.9/11.2/10/7 ppm, respectively,
Table 2). In the linear aggregate structure
Figure 6c, the presence of a carboxylic group which does not participate in hydrogen bond interactions (12.2/11.2/6.8 ppm), results in an average chemical shift of 10.4 ppm. A minor contribution of the structural model 6c, therefore, could account for the deviation of the experimental data from the computational data of the structures 6a and 6b. Moreover, the hydrophobic effect generated by the carbon chains in 6c, seems to play an antagonistic role with respect to the cyclic structure 6b.
Computations were also performed with the tetrameric caproleic acid, in a parallel orientation similar to the single crystal X-ray structures of free fatty acids [
40] and in an antiparallel orientation, in agreement with the experimental weak NOE data of the through-space proximity of the α-CH
2 and the terminal CH(9)=CH
2(10) olefinic protons. Similar methodology was used for the interpretation of the NOEs observed in the liquid state for CA, OA, ALA, and EPA [
28]. The calculated chemical shifts of the carboxylic proton for the tetrameric CA, in the parallel configuration vary between 14.3 and 13.0 ppm while in the antiparallel configuration between 13.8 and 13.2 ppm. The chemical shift difference of 1.3 ppm observed for the parallel arrangement can be attributed to the two interacting cyclic hydrogen bonds.
2.4. DFT Calculations in DMSO
The DFT calculated
1H NMR chemical shifts of the carboxylic protons with a discrete solvation molecule of DMSO were investigated in the case of a single molecule of CA, CA dimer with parallel and antiparallel arrangements (
Figure 7 and
Table 3). The representative molecular system is a caproleic acid molecule interacting with a DMSO molecule, explicitly present in the design, while the DMSO solvent is present implicitly (
Figure 7a). To this interacting pair, another one, identical to the first, was added and oriented parallel and antiparallel to it (
Figure 7b,c). These configurations were chosen to explore possible interactions between DMSO and the proton of the carboxylic group or the double bond of the caproleic acid and the proton of the carboxylic group. The results presented in
Table 3 indicate that the orientations of
Figure 7 produce practically indistinguishable δ(COOH) chemical shifts with values ranging from 13.4 to 14.2 ppm. In all cases very strong hydrogen bond interactions of the carboxylic protons with the oxygen of the DMSO molecule were observed with OH
…O distances of 1.59 to 1.62 Å and bond angles of 168.5
o to 171.2
o. These hydrogen bond distances are significantly shorter than those observed in the centro-symmetric hydrogen bond interactions through the carboxylic groups with OH
...O distances of 1.66 Å.
The results of comparing the complexation energy of the caproleic dimer in
Figure 6a and the caproleic acid-DMSO complex in
Figure 7a are very informative. For the structure 6a, the complexation energy is -21.2 kcal/mole, while it is -18.0 kcal/mole for the 7a. Given that the centro-symmetric hydrogen bond is double while in the caproleic-DMSO complex, only one hydrogen bond is formed, DMSO seems to be the most potent antagonist for this interaction.
Similar results were obtained with the α-linolenic acid. The OH
…O hydrogen bond distance (1.64 Å), the O–H
…O bond angle (168.5
o) and the COOH chemical shift (δ = 13.43 ppm) are indicative of a very strong intermolecular hydrogen bond with a single solvation molecule of DMSO (
Table 3 and
Supplementary Figure S2).