
==== Front
J Org Chem
J Org Chem
jo
joceah
The Journal of Organic Chemistry
0022-3263
1520-6904
American Chemical Society

39178339
10.1021/acs.joc.4c00879
Featured Article
Influence of Selective Deoxyfluorination on the Molecular Structure of Type-2 N-Acetyllactosamine
https://orcid.org/0000-0002-3040-8467
Kurfiřt Martin *†⊥
https://orcid.org/0000-0002-9968-1082
Št’astná Lucie Červenková †
https://orcid.org/0000-0002-4495-0070
Dračínský Martin ‡
Pohl Radek ‡
Císařová Ivana §
https://orcid.org/0000-0002-2930-8880
Sýkora Jan #
https://orcid.org/0000-0002-5030-8789
Balouch Martin ∇
https://orcid.org/0009-0002-6038-4264
Baka Michal ∥•
https://orcid.org/0000-0003-1262-2344
Hamala Vojtěch †⊥
https://orcid.org/0000-0003-4462-1469
Cañada F. Javier °□
https://orcid.org/0000-0003-3027-7417
Ardá Ana ¶△
https://orcid.org/0000-0001-5421-8513
Jiménez-Barbero Jesús ¶△▲□
https://orcid.org/0000-0001-5360-1035
Karban Jindřich †
† Institute of Chemical Process Fundamentals, Czech Academy of Sciences, Rozvojová 1/135, CZ-165 00 Praha 6, Czech Republic
‡ Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo náměstí 542/2, CZ-160 00 Praha 6, Czech Republic
§ Department of Inorganic Chemistry, Faculty of Science, Charles University in Prague, Hlavova 8, CZ-128 43 Praha 2, Czech Republic
∥ Institute of Entomology, Biology Centre of the Czech Academy of Sciences, Branišovská 31, 370 05 České Budějovice, Czech Republic
⊥Department of Organic Chemistry, #Department of Analytical Chemistry, ∇Department of Chemical Engineering, and •Department of Food Analysis and Nutrition, University of Chemistry and Technology, Prague, Technická 5, 166 28 Prague 6, Czech Republic
° Centro de Investigaciones Biológicas Margarita Salas, Ramiro de Maeztu 9, 28040 Madrid, Spain
¶ CICbioGUNE, Basque Research & Technology Alliance (BRTA), Bizkaia Technology Park, Building 800, 48162 Derio Bizkaia, Spain
△ Ikerbasque, Basque Foundation for Science, Plaza Euskadi 2, 48013 Bilbao Bizkaia, Spain
▲ Department of Organic and Inorganic Chemistry, Faculty of Science and Technology, University of the Basque Country, EHU-UPV, 48940 Leioa, Spain
□ CIBER de Enfermedades Respiratorias (CIBERES), Avda Monforte de Lemos 3-5, 28029 Madrid, Spain
* Email: kurfirt@icpf.cas.cz
23 08 2024
06 09 2024
89 17 1187511890
12 04 2024
19 08 2024
27 06 2024
© 2024 The Authors. Published by American Chemical Society
2024
The Authors
https://creativecommons.org/licenses/by/4.0/ Permits the broadest form of re-use including for commercial purposes, provided that author attribution and integrity are maintained (https://creativecommons.org/licenses/by/4.0/).

N-Acetyllactosamine is a common saccharide motif found in various biologically active glycans. This motif usually works as a backbone for additional modifications and thus significantly influences glycan conformational behavior and biological activity. In this work, we have investigated the type-2 N-acetyllactosamine scaffold using the complete series of its monodeoxyfluorinated analogs. These glycomimetics have been studied by molecular mechanics, quantum mechanics, X-ray crystallography, and various NMR techniques, which have provided a comprehensive and complete insight into the role of individual hydroxyl groups in the conformational behavior and lipophilicity of N-acetyllactosamine.

Ministerstvo Å kolstvÃ­, MlÃ¡deÅ¾e a TelovÃ½chovy 10.13039/501100001823 90254 IOCB Tech Foundation NA NA Ministerio de Ciencia e InnovaciÃ³n 10.13039/501100004837 PID2021-123781OB-C22 Ministerio de Ciencia e InnovaciÃ³n 10.13039/501100004837 PID2021-123781OB-C21 GrantovÃ¡ Agentura CeskÃ© Republiky 10.13039/501100001824 23-05146S Ministerstvo Å kolstvÃ­, MlÃ¡deÅ¾e a TelovÃ½chovy 10.13039/501100001823 LTC20052 Ministerstvo Å kolstvÃ­, MlÃ¡deÅ¾e a TelovÃ½chovy 10.13039/501100001823 90255 document-id-old-9jo4c00879
document-id-new-14jo4c00879
ccc-price
==== Body
pmcIntroduction

Oligosaccharides and polysaccharides, generally referred to as glycans1 in glycobiology, form a dense three-dimensional layer at the cell surface—glycocalyx—and they also occur abundantly in the extracellular matrix. Thanks to their enormous structural complexity and variability, glycans encode information that is read and translated by lectins and other glycan-binding proteins.2,3 To understand the molecular basis of protein−glycan recognition and exploit this understanding in the design of lectin-targeting therapeutics, glycochemists routinely employ synthetically modified glycans or glycomimetics.4−6 Conformational preferences of both natural and glycomimetic lectin ligands have a great impact on the binding energetics of glycan−protein interactions, because the nature of the glycosidic linkage provides most glycans with a significant degree of flexibility, resulting in highly dynamic conformational equilibria in solution, whereas only a limited subset of conformers is found in glycan−protein complexes.7,8 Indeed, the improvement of the conformational preorganization of a ligand to increase its inhibitory potency is an important factor in the design of glycomimetics.9 It is much more difficult to study the conformational behavior of glycans than of other biomolecules. X-ray crystallography is of limited informative value because the high flexibility of glycans complicates the formation of crystals suitable for a single-crystal X-ray diffraction analysis.10 NMR spectroscopy is very useful, but its routine application is hampered by severe degeneration of chemical shifts in glycans, which in turn requires sophisticated NMR techniques, the use of isotopically labeled glycans or an insertion of a paramagnetic tag.11,12

Over the last decades, the incorporation of fluorine has consistently belonged to the most common strategies in the development of drugs, agrochemicals and other bioactive compounds for research or industry. The replacement of an OH group with fluorine (deoxyfluorination) has served many purposes in the design of glycomimetics.13 A deoxyfluorinated carbohydrate shows a negligible steric deviation from the parent compound because the fluorine atom and the hydroxyl group have a comparable size and electronegativity, and the C−F and C−OH bonds are similar in terms of their length and polarity.13 However, deoxyfluorination removes the hydrogen-bond-donating capacity of the deoxyfluorinated position, reduces its hydrogen-bond-accepting ability, influences the electron density of adjacent atoms and reduces the polar surface area. All of these effects can, in principle, influence the conformational behavior of a deoxyfluorinated glycan. In general, deoxyfluorination has a minimal impact on the ring conformation of hexopyranoses or their simple glycosides as they retain their usual 4C1 conformation.14 The exocyclic fluoromethyl group tends to adopt the gg conformation in gluco-configured fluoroanalogs,15,16 while the gt conformer has been proposed as the most populated for galacto-configured molecules.16−19 In contrast to deoxyfluorinated monosaccharides, the conformational behavior of deoxyfluorinated oligosaccharides has been rarely investigated, which in part results from challenges associated with the synthesis of a set of deoxyfluorinated glycans.

LacNAc 1, also termed as type-2 N-acetyllactosamine (Galβ1−4GlcNAc, Figure 1A), is a biologically highly relevant disaccharide and a ubiquitous building block of various mammalian N- and O-linked glycans.20,21 In N-glycans, the LacNAc motif usually elongates the branching mannoses of the pentasaccharide N-glycan core (Figure 1B), whereas in O-glycans, it is frequently attached to GalNAc or Gal, forming a part of the extended O-glycan cores (Figure 1C). In both of these glycan types, the LacNAc motif can be repeated, forming an oligo- or poly-N-acetyllactosamine chain, which functions as a backbone for additional glycan modifications. Therefore, the type-2 LacNAc is a part and precursor of several glycan-terminating structures, including Lewis blood-group determinants Lewis X, sialyl-Lewis X, and Lewis Y, and type-2 human blood groups H, A and B,22 which are involved in processes such as pathogen adhesion,23 the spread of metastases, or fertilization.24 The LacNAc motif itself acts as a binding ligand for various carbohydrate-binding proteins such as the plant lectins isolated from Solanum tuberosum and Ricinus communis,(25) or, importantly, for galectins,26 which participate in tumor-associated or immune-system-related processes in humans. Moreover, LacNAc and its related structures can also be found in human milk oligosaccharides (HMOs), where they show protective and prebiotic effects.27

Because of the biological relevancy of LacNAc, detailed knowledge of its structure is of high importance. Diverse studies have addressed the conformational behavior of LacNAc 1(29,30) and related β-1,4-linked disaccharides such as lactose (Galβ1−4Glc),17,31−33 cellobiose (Glcβ1−4Glc)34−36 and chitobiose (GlcNAcβ1−4GlcNAc).37−40 These studies are in general agreement that β-1,4-linked disaccharides predominantly adopt the Φ/Ψ g−/g+ conformation of the glycosidic linkage accompanied by minor g−/g− and g+/g+ conformers (Figure 2).41 These conformational minima are mostly governed by stereoelectronic interactions42 and hydrogen bonds.43,44

In this work, we aspire to increase the knowledge of the behavior of LacNAc-derived glycomimetics by studying the influence of monodeoxyfluorination on its conformational states and lipophilicity. To fulfill such aspiration, we have employed a complete series of the systematically monodeoxyfluorinated LacNAcβ1-OMe (LN, 2) analogs 3−8 (Figure 3), which have recently been prepared in our laboratory.45 The methyl β-glycosides (LNs) have been selected for their synthetic accessibility and the presence of β-glycosidic linkage in the repeating poly-N-acetyllactosamine chains of natural glycans. The LNs 2−8 have been studied by a combination of computational and NMR approaches utilizing 13C chemical shifts, J-coupling constants, NOE/ROE correlations and 1H NMR temperature dependences of exchangeable protons. The series has also been studied by single-crystal X-ray diffraction analysis, as we succeeded in obtaining suitable crystals for all deoxyfluorinated LacNAc analogs. The combination of these techniques made it possible to decipher the influence of the systematic monodeoxyfluorination on the conformational behavior of the glycosidic linkage, exocyclic groups at C6/6′ and the acetamido group at C2. The stablest conformers of the compounds 3−8 were subsequently utilized for the calculation of their 1-octanol/water partition coefficients (log P) using the COSMO-RS methodology,46,47 which showed excellent agreement with experimental data.

Results and Discussion

Molecular Mechanics

Initially, the conformational behavior of the LNs 2−8 has been investigated by systematic relaxed potential energy scans of two glycosidic-bond dihedral angles Φ (O5′−C1′−O4−C4) and Ψ (C1′−O4−C4−C3) using molecular mechanics (Figures S1−S7).17,31,48 The LNs 2 and 4−8 showed highly similar Φ/Ψ adiabatic population diagrams with a high occupancy of the Φ/Ψ g−/g+ region, covering 2.05−2.70% of Φ/Ψ conformational space at the 90% probability level (Table 1). In contrast, 3F-LN 3 covers approximately 4.4% of conformational space and exhibits a high population of conformers outside the defined regions (Table 1), predicting its higher conformational flexibility around the glycosidic linkage. Such different computational results probably reflect the disruption of the inter-residue O5′···H−O3 hydrogen bond (Figure 4), commonly found in β-1,4-linked disaccharides.29,38,43,44,49

X-ray Crystallography

We have succeeded in obtaining crystals of the disaccharides 3−8 suitable for X-ray diffraction analysis (Figures S8−S13). As expected, the LNs 3−8 crystallized in the Φ/Ψ g−/g+ conformation in accordance with the global minimum found in molecular mechanics, corresponding to the crystal structure of the parental nonfluorinated LacNAc 1.29 To gain a clear insight into the effects of individual fluorine incorporations, we compared the lengths of inter-residue O5′···H−O3 hydrogen bonds and the dihedral angles Φ, Ψ, Φ2, Ω1, Ω2, Ω3 and Ω4 (Figures 2, 4 and Table 2). The dihedral angle Φ2 [O5−C1−O−CH3] describes the conformation of the OMe glycosidic linkage, the torsions Ω1 [O5−C5−C6−O6] and Ω2 [O5′−C5′−C6′−O6′] describe the orientation of exocyclic hydroxyls, and the remaining torsions Ω3 [C1−C2−N−C] and Ω4 [C2−N−C−C] define the conformation of the acetamido group.

Interestingly, the compounds 4−7 showed analogous trends in the length of the O5′···H−O3 hydrogen bond and the glycoside torsions Φ, Φ2 and, with minor exceptions, also Ψ (Table 2). Compared to α-d-LacNAc29 (O5′···O3 = 2.79 Å), 2′F-LN 5 and 6F-LN 4 showed longer O5′···O3 distances, namely 3.04 and 3.02 Å, respectively. This trend seems to indicate that the deoxyfluorination of the positions 2′ and 6 somehow weakens the inter-residue hydrogen bond. It is tempting to speculate that the longer O5′···O3 distances are related to an increase in flexibility around the glycosidic linkage. In contrast, the fluorination of the 3′-position (6) and 6′-position (8) has significantly decreased the O5′··· O3 distance to 2.68 and 2.71 Å, respectively. Again, it is tempting to speculate that this shortening of the O5′···O3 distance indicates some rigidification of the glycosidic linkage. 4′F-LN (7) has shown approximately the same O5′··· O3 distance as 1 (2.82 Å), which implies that 4′-fluorination has not influenced the glycoside conformational behavior. The significant variations between the three-dimensional structures of the limiting compounds 2′F-LN 5 and 3′F-LN 6 are illustrated by their overlay, as depicted in Figure 5A.

In the series 5 → 4 → 7 → 8 → 6, the absolute values of Φ torsions gradually decreased, which can be illustrated as a dihedral closing (Figure 6A). This closing significantly correlated with the decreasing O5′···O3 distance (R2 = 0.91, Figure 6B), indicating the interdependence of these two molecular descriptors. Similarly, Ψ torsion exhibited a gradual decrease with the decreasing O5′···O3 distance (Figure 6A), also with a high coefficient of determination (R2 = 0.84, Figure 6C). Surprisingly, the absolute value of Φ2 torsion also showed a gradual increase in the series 5 → 4 → 7 → 8 → 6, represented by its opening (Figure 6A), strongly correlating with the O5′···O3 distance (R2 = 0.87, Figure 6D). Therefore, we may hypothesize that the conformation around the glycosidic linkage in β-1,4-linked oligosaccharides is somehow influenced by the modulation of the O5′···H−O3 hydrogen bond strength of a neighboring disaccharide unit. In 3F-LN 3, the value of Ψ was the lowest in the series 3−8. This exceptional behavior of 3F-LN 3 in comparison with the remaining compounds illustrates that the precise geometry within the Φ/Ψ g−/g+ conformational region of LacNAc is somehow modulated by the disruption of the O5′···H−O3 hydrogen bond. The effect of 3-fluorination can be illustrated by the overlay of the crystal structures of 3F-LN 3 and LacNAc 1 (Figure 5B), which shows a change of the relative position of the gluco- and galacto-configured rings upon fluorination.

The torsion Ω1, describing the conformation of the exocyclic 6-hydroxymethyl/fluoromethyl groups in the GlcNAc rings, was independent of the deoxyfluorination pattern, with values being close to −60° (Table 2 and Figure 7A, gg conformation). Similarly, the torsion Ω2, describing the conformation of the exocyclic 6′-hydroxymethyl/fluoromethyl groups in the Gal moieties, showed similar values, close to 60° (Table 2 and Figure 7B, gt conformation) in compounds 4−8. Intriguingly, 3F-LN 3 had Ω2 = 173°, which is dramatically different from those of the remaining disaccharides and close to the staggered tg orientation (Table 2, Figures 5B and 7C). Although the crystallization of 3 as the 6′-tg conformer could be accidental because both gt and tg rotamers are present in Gal moieties in solution, the increased preference for the 6′-tg conformation was independently determined by NMR in DMSO-d6 solution, discussed in the next section. The values of the remaining torsions, Ω3 and Ω4, were consistent with the (Z)-anti conformation of the acetamido group (Figure 7D), which is regarded as stabler than the second possible (Z)-syn conformation (Figure 7E).50

NMR Analysis

To understand the behavior of the LNs 2−8 in solution, we performed a quantum-mechanics-assisted NMR analysis. The analysis was based on the evaluation of chemical shifts, indirect J-coupling constants, qualitative NOEs (ROEs), and the temperature dependence of the chemical shifts of exchangeable protons.52 This investigation was carried out in DMSO-d6 solution as it enables facile observation of exchangeable protons.

Chemical Shifts of Skeletal Protons and Carbons

Deoxyfluorination shifted the geminal protons by approximately 1 ppm downfield and the directly attached 13C carbons by more than 20 ppm downfield in all cases (Table S2, green cells). Similarly, vicinal protons separated from the fluorine atom by three single bonds were shifted downfield by 0.18−0.41 ppm. (Table S2A, blue cells). Conversely, 13C carbons separated from the fluorine by two or three bonds were generally shifted upfield by 0.3−5.9 ppm, although there were a few exceptions of a minor downfield shift by 0.1−0.7 ppm (Table S2B, blue cells). These observations are also in agreement with the data reported for other fluorinated sugars.53,54

Interestingly, both 2′F-LN 5 and 3F-LN 3 showed selective 1H and 13C long-range fluorine-induced perturbation of chemical shifts (Δδ Table S2, yellow cells). In 5, the protons H6proR and H6proS, separated from 2′-fluorine by seven bonds, were shifted upfield by 0.07 and 0.09 ppm, respectively. No other proton in 5, separated from fluorine by six (H6′proR, H6′proS, H3, H5) or seven bonds (H2), showed Δδ of such magnitude. Additionally, the carbon C6 exhibited an upfield shift by 0.7 ppm, the highest value among the carbons separated from 2′-fluorine by five or six bonds. These long-range Δδs indicate that the fluorine at the 2′-position influences the carbon and protons at the 6-position via a through-space effect. This conclusion is supported by the relatively short distance between the fluorine atom and H6proR (2.64 Å) or C6 (3.35 Å) in the X-ray structure of 2′F-LN 5 (Figure 8A). Furthermore, we have detected a through-space 2′-fluorine−carbon C6 coupling in the disaccharide 2′F-LN 5. This coupling has been described in detail elsewhere.55 Similarly, 3F-LN 3 has shown a significant long-range Δδs of H5′ and C5′, also associated with a short distance between the 3-fluorine nucleus and H5′ (3.37 Å) or C5′ (3.58 Å) in the crystal structure of 3F-LN 3 (Figure 8B). Therefore, these observed long-range Δδs in 5 and 3 are in accordance with the preference for the Φ/Ψ g−/g+ conformation in DMSO-d6 solution.

13C NMR Chemical Shift-Based Conformational Analysis

To support the interpretation of the NMR data, the geometries of the major conformational minima of the LNs 2−8 were optimized by means of density functional theory (DFT) to generate the conformers A−L (Figure 9A). These geometries had been selected to probe the conformations around the glycosidic linkages, the exocyclic hydroxyls, and the acetamido moieties systematically. The comparison of the zero-point-corrected free energies (Figures S14 and S15) had been identified as inappropriate for conformational analysis because of the tendency of DFT methods to overestimate the strength of hydrogen bonds.56,57 Therefore, the conformational analysis was based on the correlation between the DFT-calculated and the experimentally obtained NMR parameters (Tables S14−S20). Initially, we focused on 13C chemical shifts given their sensitivity to the conformational changes of the glycosidic linkage and the exocyclic CH2OH/F groups, and their relative insensitivity to subtle alterations in the orientation of exchangeable protons.58−60 The deviations between a set of experimental and calculated 13C chemical shifts were expressed as a mean-square deviation (MSD, Figure 9B).

Generally speaking, the calculated 13C NMR chemical shifts of the conformer Φ/Ψ g−/g+ (A) showed significantly better agreement (lower MSDs) with the experimentally obtained 13C NMR shifts than those predicted for the conformers Φ/Ψ g−/g− (H) and Φ/Ψ g+/g+ (I). This trend in MSDs confirms the preference for the Φ/Ψ g−/g+ glycoside conformation in all compounds in accordance with X-ray analysis and the reported data.17 Moreover, as expected, the predicted 13C NMR shifts of the 6-tg conformers J, K and L exhibited significantly worse agreement (higher MSDs) with the experimental 13C NMR shifts than those of the 6-gg conformers A, D and E, and the 6-gt conformers B, C and F. This observation further confirms the general preference for the 6-gg and 6-gt orientations over the 6-tg orientation of the GlcNAc exocyclic CH2OH group.17,18 Similarly, the 6′-gg conformers (E and F) exhibited higher MSDs than the 6′-tg (A and B) and 6′-gt (C and D) rotamers, also in accordance with the expectations for Gal moieties.17,18 The acetamido group showed a preference for the (Z)-anti conformation (A) over the (Z)-syn conformation (G) in accordance with the X-ray analysis and reported data for acetamide sugar moieties.50

Apart from these general trends, the methodology enabled us to identify subtle differences in the conformational behavior of the LNs 2−8, caused by selective deoxyfluorination. The fluorinated disaccharides 6F-LN 4, 2′F-LN 5 and 4′F-LN 7 provided comparable MSDs in all the conformers A−L as LN 2, indicating a minimal impact of the fluorination of these positions on the conformational behavior. In comparison with the LN 2, 3F-LN 3 showed higher MSD values of the conformers Φ/Ψ g−/g+ (A−F), and lower MSD values of the conformers Φ/Ψ g−/g− (H) and Φ/Ψ g+/g+ (I), which is consistent with the 3F-LN 3 being more conformationally flexible than the LN 2, probably because of the disruption of the O5′···H−O3 hydrogen bond. The 3′F-LN 6 exhibited higher MSDs of the Φ/Ψ g−/g+ conformers A−D than the LN 2, also implying the presence of additional conformational dynamics in solution. For the 6′F-LN 8, the best fit between the experimental and calculated shifts was deduced for the conformer D, with the 6′-gt conformation of the galactopyranoside exocyclic fluoromethyl group. This fact is in agreement with previous observations,18 which reported that the fluorine atom increases the gauche effect, thus stabilizing the gt conformer over the tg conformer in the galacto-configuration (Figure 7 and 9).

In summary, the correlation between the DFT-calculated and the experimentally determined 13C chemical shifts has provided a detailed insight into the conformational behavior of the LNs 2−8 in solution concerning the conformation of the glycosidic linkage, exocyclic hydroxyls, and acetamido group. Fittingly, the results obtained are also in agreement with the corresponding crystal structures.

J-Coupling-Based Conformational Analysis

To complement the 13C-based methodology discussed above, the conformations of the exocyclic hydroxymethyl/fluoromethyl groups were investigated using vicinal J-couplings between the diastereotopic protons at the 6- and 6′-positions and the corresponding vicinal protons at the 5- and 5′-positions.18 The conformation of these groups can be defined as a fast dynamic equilibrium between three staggered conformers denoted as gg, gt and tg (Figure 7). Each of these conformers is characterized by the specific values of the limiting 3JH5-H6proR and 3JH5-H6proS coupling constants, which makes it possible to calculate the gg/gt/tg ratio from the corresponding experimental 3J coupling constants (Table 3). We have calculated the limiting couplings of gg/gt/tg conformers by DFT (Table S21) and used them for the estimation of the populations of individual conformations (fgg, fgt, ftg values). The assignment of proR and proS protons is described in the Supporting Information

For the GlcNAc units, the LNs 2 and 4−8 exhibited the couplings 3JH5-H6proR (4.6−5.1 Hz) and 3JH5-H6proS (1.8−2.4 Hz), corresponding to the strong preference of 6-gg (63−66%) over 6-gt (30−35%) conformers, and a negligible population of the 6-tg conformer (≤5%). The 3F-LN 3 showed decreased 3JH5-H6proR coupling (4.0 Hz), consistent with a slightly enhanced population of the 6-gg conformer (71%). The negligible population of the tg conformer in the GlcNAc unit is explained by an unfavorable steric 1,3-interaction (syn-pentane-type) between the O6H and O4H hydroxyl groups and the absence of a stabilizing stereoelectronic gauche effect.18

For the Gal units, the compounds 4−7 exhibited comparable values of the diagnostic couplings 3JH5′-H6′proR (6.7−7.3 Hz) and 3JH5′-H6′proS (5.2−5.9 Hz), corresponding to the gg/gt/tg conformer ratio of ca. 15/55/30 in accordance with the previously published conformational analysis of lactose analogs.17 The 6′F-LN 8 showed the largest difference between the coupling constants 3JH5′-H6′proR (7.9 Hz) and 3JH5′-H6′proS (3.2 Hz), corresponding to the gg/gt/tg ratio of 17/76/7. The strong preference for the 6′-gt is a result of the enhancement of the gauche effect via fluorine electronegativity, which is in agreement with the 13C-based conformational analysis described above and reported data.18 Interestingly, the 3F-LN 3 was the only analog showing a larger 3JH5′-H6′proS coupling constant (7.1 Hz) than 3JH5′-H6′proR (5.9 Hz), corresponding to the highest population of the 6′-tg conformer (50%). Such an observation is in accordance with the above-discussed X-ray analysis because 3F-LN 3 crystallized in the 6′-tg conformation (Figure 5B). Furthermore, the fact that the conformation around C6′H2OH is affected by the fluorination of the 3-position indicates an interaction between the hydroxyls at the 3- and 6′-positions, consistent with the temperature-dependent 1H NMR analysis discussed below.

The conformational behavior of an acetamido group in carbohydrates is usually described as a dynamic equilibrium between (Z)-anti and (Z)-syn rotamers (Figure 7), where (Z)-anti is considered the stabler geometry.50 However, a qualitative analysis of the JH2-NH couplings could indicate the prevalence of the (Z)-syn conformation in solution, because the experimental JH2-NH coupling values were 8.5−9.1 Hz (Table 3), closer to the JH2-NH ≈ 8.3 Hz values estimated by DFT calculations for the (Z)-syn conformer, and relatively far from the limiting coupling value JH2-NH ≈ 11.6 Hz, computed for the (Z)-anti geometry (Table S22). To explain such a discrepancy, the conformational behavior of the acetamido group was investigated by DFT in the presence of an explicit DMSO molecule (Figure S18). Interestingly, the explicit solvation significantly influenced the acetamide geometry and the calculated JH2-NH coupling. Therefore, the implicit solvation model is not an appropriate approximation of solvation effects concerning acetamido protons, which are involved in hydrogen bonding with the solvent molecules.61 In fact, the analysis of NOE contacts between the NH amide proton and H1, H2 and H3 ring protons (see below) was crucial for the proper definition of the preferential conformation of the acetamide.

Analysis of ROESY NMR Spectra

Complementarily to the 13C NMR-based and J-coupling-based results, the compounds 2−8 were studied using 1H-1H ROESY NMR to estimate the spatial proximity of protons, which can be related to the presence of specific conformations (Figure 10).

The LNs 3−8 exhibited strong H1′/H4 contacts, confirming the expected dominance of the Φ/Ψ g−/g+ conformation (Table 4, Figure 10). The H1′-H3/H5 contacts were significantly weaker, yet detectable, indicating the presence of a minor population of the Φ/Ψ g−/g− conformation in all fluoro-analogs. The presence of the H2′-H4 contact, characteristic of the Φ/Ψ g+/g+ conformation, could only be assessed for the 2′F-LN 5 and 3′F-LN 6 because of the signal overlap in the other compounds. Nevertheless, none of those compounds showed a clear H2′-H4 NOE correlation, indicating that the population of the Φ/Ψ g+/g+ conformer is marginal, if present at all.

The LNs 4−8 showed a stronger H4/H6proR NOE than H4/H6proS, indicating a preference for the 6-gt conformer over the 6-tg in the GlcNAc unit, in agreement with the DFT analysis of both 13C chemical shifts and J-couplings described above. In LN 2 and 3F-LN 3, these contacts could not be analyzed separately because of the overlap of the protons H6proR and H6proS. A similar issue with the overlapping protons H6′proR and H6′proS was found in the compounds 2 and 4−6. In 3F-LN 3, 4′F-LN 7 and 6′F-LN 8, however the introduction of the fluorine atom caused the separation of the signals of the H6′proR and H6′proS protons, enabling the identification of a stronger NOE for H4′/H6′proS than for H4′/H6′proR (Figure S17). This confirmed the preference of the 6′-gt over 6′-gg conformation, which was also in agreement with the DFT analysis of both 13C chemical shifts and J-couplings. Furthermore, the LNs 3−8 showed a stronger NOE correlation between the protons (N)H and H1/H3 than between the protons (N)H and H2, indicating the expected preference for the (Z)-anti over (Z)-syn acetamide conformation. The observation of weak (N)H-(O3)H contacts indicated that the acetamide Ω3 torsion angle slightly deviated from the ideal arrangement (the antiperiplanar geometry of H2 and NH protons), in accordance with DFT analysis with explicit DMSO solvation (Figure S18).

Analysis of Hydroxyl-Solvent Interactions by 1H NMR and HSQC-TOCSY

The 1H chemical shifts of exchangeable protons exhibit a linear dependence on temperature, which can be employed for the identification of intra- and intermolecular hydrogen bonds in oligosaccharides.43,44 Strongly solvated hydroxyl protons show higher values of −Δδ/ΔT temperature coefficients than those engaged in intramolecular hydrogen bonds.62,63

The nonfluorinated LN 2 shows marked differences in the temperature coefficients −Δδ/ΔT of individual hydroxyls as determined by variable temperature NMR in DMSO-d6 (Table 5). The lowest coefficient (2.0 ppb/K) has been observed for O3−H, which corresponds to the presence of the O5′···H−O3 hydrogen bond in agreement with X-ray analysis. On the other hand, the protons residing at the hydroxyls O2′−H and O3′−H have higher temperature coefficients (6.8 and 7.2 ppb/K respectively), which confirms their significant solvation. Such results are consistent with previous reports describing that the hydroxyl groups in vicinal equatorial diols cannot form direct intramolecular hydrogen bonds with each other, although they can form strong bidentate hydrogen bonds with solvent molecules.64,65 The temperature coefficients of the primary hydroxyls O6−H and O6′−H (5.9 and 4.9 ppb/K, respectively) are quite similar, although the observed values indicate a higher engagement of the hydroxyl O6′−H in intramolecular hydrogen bonding.

Our series 3−8 provides a unique opportunity to probe intramolecular hydrogen bonding patterns in the LN system (Table 5). First of all, the NH protons on acetamido groups show only negligible differences in −Δδ/ΔT coefficients, suggesting that NHAc solvation is not affected by monofluorination. The presence of fluorine generally results in a decrease in the −Δδ/ΔT temperature-dependence coefficients of the vicinal hydroxyls. Taking into consideration that DMSO-d6 could only act as a hydrogen-bond acceptor,66 deoxyfluorination possibly decreases the hydrogen-bond-donating capacity of the corresponding equatorial vicinal hydroxyl, decreasing the stability of the O−H···O=SMe2 hydrogen bond, in accordance with previous reports.67−69 When compared to the nonfluorinated LN 2, both 2′F-LN 5 and 4′F-LN 7 show basically the same decrease of the O3′−H coefficients −Δδ/ΔT by 1.2 ppb/K, indicating that the fluorine configuration—axial or equatorial—does not play a role in reducing the hydrogen-bond-donating capacity of the equatorial vicinal hydroxyls.

Apart from the vicinal effect, the fluorine-induced changes of −Δδ/ΔT coefficients also exhibited small, yet detectable, long-range effects, reflecting the subtle alteration of the LN molecular structure. Primarily, the coefficients of O3−H showed a dependence on the position of fluorine. The highest coefficient was observed for 2′F-LN 5 (2.6 ppb/K) whereas the lowest one for 3′F-LN 6 (1.7 ppb/K). The remaining fluorinated disaccharides 4, 7 and 8 were in between these values, with their coefficients comparable to LN 2 (2.0 ppb/K). Fittingly, the −Δδ/ΔT coefficients of O3−H correlated with the X-ray-determined O5′···O3 distance (Figure 11, R2 = 0.81), strongly indicating that the trends observed in the crystal structures are also relevant for DMSO-d6 solution. Accordingly, 2′F-LN 5, with the longest O5′···O3 distance according to the X-ray, showed the highest temperature coefficient for O3−H, suggesting a weaker inter-residual O5′···H−O3 hydrogen bond than LN 2 also in solution. Conversely, both 3′F-LN 6 and 6′F-LN 8, with the shortest O5′···O3 distance in solid state, seem to have stronger O5′···H−O3 hydrogen bonds in solution than LN 2 according to their lower −Δδ/ΔT. The linear relationship between the −Δδ/ΔT coefficients of O3−H and the O5′···O3 distance (Figure 11A) may be employed for an estimation of the intramolecular hydrogen-bond length in other β-1,4-linked oligosaccharides, in cases where an X-ray analysis is not possible.

Valuable information has also been extracted from the temperature coefficients of O2′−H, with the lowest −Δδ/ΔT coefficient (4.9 ppb/K) observed for 3′F-LN 6, which corresponds to a decrease by 1.9 ppb/K in comparison with LN 2. This decrease is higher than one would expect from solely the fluorine vicinal effect (decrease by 1.2 ppb/K for 3′−OH in 5 or 7). Similarly, 6′F-LN 8 features the second lowest −Δδ/ΔT coefficient of O2′−H (5.5 ppb/K), which is decreased by 1.3 ppb/K relative to LN 2, despite the separation of the proton concerned by seven bonds from fluorine. Such a decrease strongly suggests the presence of a transient O2′−H···O6 hydrogen bond in both 3′F-LN 6 and 6′F-LN 8 in DMSO-d6 solution, which can also be correlated with stronger O5′···H−O3 hydrogen bonds implying a slight reduction of glycoside linkage conformational flexibility in these two analogs. Fittingly, the presence of the O2′−H···O6 hydrogen bond in 3′F-LN 6 (Figure 12A) has been independently confirmed by the observation of HSQC-TOCSY cross correlation between O2′−H and C6 (Figure 12B). In this experiment, no correlation between O2′−H and C4 or C5 has been observed, confirming the presence of a proton/proton coupling transferred through the hydrogen bond.

The formation of a strong O6′−H···O4′ or O4′−H···O6′ hydrogen bond, which is in principle possible in galacto-configured carbohydrates, has been excluded because the fluorination of the 4′-position (7) or the 6′-position (8) did not cause any significant increase in the solvation of O6′-H or O4′-H, respectively. This is in accordance with a very minor population of the 6′-gg conformation, as determined by both 13C- and J-coupling-based NMR analyses (Table 3), which would be stabilized by these hydrogen bonds. Conversely, the DFT calculations suggested the existence of a weak O6′−H···O3 hydrogen bond (Figure 11B), which could explain the somehow lower solvation of O6′−H (−Δδ/ΔT = 4.9 ppb/K) in comparison with that of O6−H (−Δδ/ΔT = 5.9 ppb/K) in LN 2. The presence of such a hydrogen bond is further supported by the increased O6′−H solvation in 3F-LN 3 (−Δδ/ΔT = 5.7, Table 5), where this hydrogen bond is disabled by 3-deoxyfluorination. Such an observation is also consistent with the prevalence of the 6′-gt conformer in 4−8 and the preference for the 6′-tg in 3F-LN 3 (Tables 2, 3 and 4).

Relationships between Molecular Structure and Lipophilicity

Lipophilicity is a key parameter for the design of new pharmacologically active compounds.70 The lipophilicity of the fluorinated LNs 3−8, expressed as an octanol/water partition coefficient (log P), was previously determined by the shake-flask procedure using the measurement of 19F NMR resonances in the octanol and water phases (Table 6).45 In this work, the log P values of the compounds 3−8 have been estimated via the COSMO-RS methodology,46 using the stablest conformers of the LNs 3−8. The calculation has furnished data in excellent agreement with the shake-flask experiment (Table 6, Figure 13A).

Both the experimental and computational COSMO-RS analysis has identified 2′F-LN 5 as the most lipophilic and 6F-LN 4 as the second most lipophilic molecule in the series 3−8. Conversely, 3F-LN 3 has been identified as the least lipophilic in the series. The remaining compounds 6, 7 and 8 have shown similar intermediate lipophilicity (Figures 13A and 13B). Although DMSO and water solvation should be different, it is tempting to analyze the lipophilicity data by means of the determined −Δδ/ΔT temperature coefficient. Nevertheless, this analysis should be taken with caution. The relatively low lipophilicity of 3F-LN 3 can be rationalized by the previously discussed temperature dependence of exchangeable-proton chemical shifts. Such analysis has confirmed the involvement of the hydroxyl O3−H in the inter-residual O5′···H−O3 hydrogen bond, which should decrease its hydration. Therefore, 3-deoxyfluorination to provide 3F-LN 3 does not increase lipophilicity to the same extent as deoxyfluorination in other positions. This observation also indirectly indicates that the inter-residual O5′···H−O3 hydrogen bond persists even in water solution.71

Conversely, the relatively high lipophilicity of the compounds 2′F-LN 5 and 6F-LN 4 is consistent with the high −Δδ/ΔT temperature coefficients of the protons O2′−H and O6−H in the LN 2 (6.8 and 5.9 ppb/K, respectively), indicating their strong solvation and explaining the significant increase in lipophilicity upon their deoxyfluorination. Interestingly, the relatively high increase in lipophilicity found for 2′F-LN 5 and 6F-LN 4 is also reflected in their crystal structures because only those two LN analogs have crystallized without bound water molecules (Figures S8−S13). Moreover, the crystal structure of 6′F-LN 8 contained crystal-bound water tethering the hydroxyls O2′−H and O6−H (Figure 13C). Such bidentate O2′−H/O6−H solvation is obviously weakened in both 2′F-LN 5 and 6F-LN 4, consistent with their relatively high lipophilicity.

Conclusions

We have herein analyzed the complete series of monofluorinated LacNAcβ1-OMe analogs to investigate the role of individual hydroxyls in the conformational behavior and lipophilicity of N-acetyllactosamine. These analogs have been studied by molecular and quantum mechanics, X-ray crystallography and NMR-based techniques, which have revealed that some intramolecular inter-residue hydrogen bonds (O5′···H−O3, O2′−H···O6 and O6′−H···O3) fine-tune N-acetyllactosamine conformational behavior. The deoxyfluorination of LacNAcβ1-OMe at the 6- and 4′-positions has no significant effect on the conformation. The deoxyfluorination of the 2′-position has resulted in the weakening of the O5′···H−O3 inter-residual hydrogen bond, likely leading to a subtle flexibilization of the glycosidic linkage. The most significant glycosidic-linkage flexibilization was observed upon the deoxyfluorination of the 3-position, which was attributed to the disruption of the O5′···H−O3 hydrogen bond. The 3-fluoro-N-acetyllactosamine analog also exhibited an anomaly in the conformational behavior of the exocyclic hydroxyl at position 6′, likely caused by a disruption of the O6′−H···O3 hydrogen bond. In contrast, 3′- and 6′- deoxyfluorination resulted in a subtle glycosidic-linkage conformational rigidification, probably associated with the strengthening of the O5′···H−O3 and O2′−H···O6 hydrogen bonds in solution. Such conformational effects probably influence biological activity. For example, we previously reported that 2′F-LacNAcβ1-OMe binds human galectin-1 with significantly higher affinity than nonfluorinated LacNAcβ1-OMe.45 As 2′F-LacNAcβ1-OMe showed the weakening of the O5′···H−O3 hydrogen bond, it is tempting to speculate that glycoside bond conformational flexibilization accounts for the affinity increase. However, as protein−ligand recognition is a complex process, further investigation would be required to confirm such speculation.

COSMO-RS in-silico calculations have confirmed the previously published lipophilicity of the investigated deoxyfluorinated LacNAcβ1-OMe analogs.45 NMR and X-ray analysis have enabled a partial rationalization of the observed lipophilicity trends, indicating that hydroxyls at 2′- and 6-positions are significantly more solvated in water solution then the hydroxyl at 3-position. Our detailed conformational analysis and lipophilicity rationalization could serve as a guideline for the synthesis of novel tailored LacNAc-derived glycomimetics.

Experimental Section

Molecular Mechanics

The systematic relaxed scan of glycoside torsions Φ (H1′−C1′−O4−C4) and Ψ (C1′−O4−C4−H4) was performed using the MacroModel72,73 software implemented in Maestro modeling environment (Maestro, Schrödinger, LLC, New York, version 13.8.135, MMshare Version 6.4.135, Release 2023−4, Platform Linux-x86_64). For the discussion, the glycosidic angles were described following the definition based on heavy atoms, Φ (O5′−C1′−O4−C4) and Ψ (C1′−O4−C4−C3) for comparison with experimental data derived from crystallographic structures. The topology of investigated systems 2−8 was constructed by automatic approach using MM3 force field74 with the use of implicit GBSA75 water solvation. This solvation model was selected to allow comparison with the literature data.48 Each scan generated a grid of torsions Φ and Ψ in an interval −300° to 60° with a step of 10°, therefore 1369 structures. Four combinations of exocyclic chain conformations 6/6′: gt/gg, tg/gg, gt/gt, tg/gt were considered, generating the 5476 structures in total for each investigated system 2−8. Every generated structure was minimized using the default PRCG (Polak-Ribiere conjugate Gradient) minimization76 with maximally 2500 iterations and 0.05 convergence threshold. The final adiabatic Φ/Ψ energy diagrams were built using the lowest energy structure for each obtained Φ/Ψ point. The potential energy maps were transformed into the Φ/Ψ probability distribution by Boltzmann distribution. To compare individual compounds 2−8, the obtained data were statistically analyzed in Microsoft excel Professional Plus 2019 and plotted through GNU Octave 8.3.0.77

X-ray Crystallography

Crystals of compounds 3−8, suitable for X-ray analysis, were prepared by the vapor diffusion technique78 using the demineralized water as the solvent and tert-butanol as the antisolvent. The starting concentration of the compound was approximately 20 mM. Demineralized water with conductivity ≤0.1 μS was obtained using ROTFM-5SV reverse osmosis. All crystals showed melting points of approximately 350 °C accompanied by decomposition, preventing their precise determination. The diffraction data was collected on a Bruker D8 VENTURE Kappa Duo PHOTON 100 CMOS with monochromated Mo/Cu−Kα radiation. Structures were solved by direct methods (SHELXT79) and refined by full-matrix least-squares on F2 values (SHELXL79 or CRYSTALS80). All heavy atoms were refined anisotropically. Hydrogen atoms were localized from the expected geometry and difference electron density maps and usually were refined isotropically. ORTEP-381 was used for the presentation of the structure. The crystallographic data for the structures reported in this paper have been deposited in the Cambridge Crystallographic Data Center as a supplementary publication. Copies of the data can be obtained free of charge on application to CCDC, e-mail: deposit@ccdc.cam.ac.uk. BIOVIA Discovery Studio Visualizer v21.1.0.20298 was used to generate the associated graphics.

NMR Analysis

NMR spectra for the characterization of compounds 2−8 in DMSO-d6 were acquired using Bruker Avance 400 (1H at 400.1 MHz, 19F at 376.4 MHz, 13C at 100.6 MHz) or Bruker Avance 500 (1H at 500.1 MHz, 19F at 470.5 MHz, 13C at 125.8 MHz) at 25 °C. The 1H and 13C NMR spectra were referenced to the DMSO-d6, 2.50/39.52 ppm. The 19F NMR spectra were referenced to the line of the external standard hexafluorobenzene (δ/ppm −163.86 in DMSO-d6). Structural assignments were made with additional information from gCOSY, gHSQC, gHMBC and H,C-gHSQC-TOSCY). Diagnostic JH6′proR-H5 and JH6′proS-H5 couplings were extracted by 1D selective proton gradient ROESY, 1D selective gradient TOCSY or selective homonuclear decoupling using Bruker Avance 400 MHz, Bruker Avance III HD 600 MHz or Bruker Avance III HD 850 MHz. Bruker Avance III HD NMR spectrometers (600 and 850 MHz) were equipped with the inverse triple resonance cryo-probe with ATM module (5 mm CPTCI 1H/13C/15N/D Z-GRD). The analysis of the spatial proton contacts and the exchange contacts with residual water in DMSO-d6 was performed on Bruker Avance 400 via H,H-ROESY (roesyphpp.2) at 25 °C with 2 s of relaxation delay D1 and 200 ms pulse for ROESY spinlock (P15). The 1H temperature dependence experiments of exchangeable protons were performed on Bruker Avance 400 in the temperature range of 25−50 °C using a common zg30 pulse sequence for proton acquisition. All experiments were processed using Topspin 4.1.3 or MestreNova 12.0.3 software.

DFT Calculations

DFT calculations were performed in Gaussian16-B.0182 at the B3LYP-D3/6-311++g(d,p)83−86 level of theory with diffuse and polarization functions on all atoms and D3 empirical dispersion correction according to Grimme.87,88 All calculations were performed with implicit DMSO solvation utilizing a polarizable conductor calculation model (C-PCM).89,90 All additional calculation variables such as convergence criteria were set to default. In all optimized structures, we calculated vibration frequencies to evaluate the character of the stationary points (confirmation of the minimum). The optimized geometries were used for the calculation of 13C NMR shifts and J-couplings using the Gauge Independent Atomic Orbital (GIAO) method91−95 including the Fermi contact contribution96 to the nuclear spin−spin coupling. The isotropic parts of the 13C magnetic shielding tensors (σISO) were converted into chemical shifts (δtheor) using the linear regression method.97,98 The method is based on the linear dependence between experimentally measured chemical shifts (δexp) and theoretically obtained σISO. First, the linear regression coefficients (the slope and the intercept) between σISO and δexp have been obtained using the slope and intercept functions in Microsoft Excel Professional Plus 2019. The values of δexp were used as x_array and σISO as y_array. The slope and intercept values were later used for the calculation of δtheor according to the equation. 1.1

Subsequently, the mean square deviations (MSD) between the experimental 13C chemical shifts δexp and the theoretically obtained 13C chemical shifts δtheor conformer were calculated. BIOVIA Discovery Studio Visualizer v21.1.0.20298 was used to generate the associated graphics.

COSMO-RS Calculations

Selected DFT-minimized geometries representing the most stable conformers of compounds 3−8 (See ESI II) were utilized for the calculation of Log P using the conductor-like screening model for real solvents (COSMO-RS) method,46 which is a semiempirical approach that shows high accuracy for partition prediction.47 Each conformer selected by the previous approach was used for the subsequent DFT/B−P/TZVP vacuum99 and COSMO water optimization with fine grid option using Turbomole 6.3.100 The .cosmo files for each conformer were obtained from this procedure and subsequently used for the partitioning calculation. Using the COSMOtherm X18 software and .cosmo files of all calculated conformers, the chemical potential of the molecules in water and water-saturated octanol (27.4 molar % of water in octanol) was calculated, furnishing the partition coefficients.

Figure 1 A) The structure of N-acetyllactosamine 1 with atom numbering and its representation according to the symbol nomenclature for glycans (SNFG).28 B) An example of a complex N-glycan on mature glycoproteins containing two LacNAc units.20 C) An extended core 2 O-GalNAc glycan from human respiratory mucins containing one LacNAc unit.21

Figure 2 A) The typical glycoside linkage conformational minima of the LacNAc system Φ/Ψ g−/g+, g−/g− and g+/g+. B) The definition of the torsions Φ and Ψ, based on heavy atoms, and the corresponding terminology of the ideal staggered Φ and Ψ rotamers.

Figure 3 Structures of the LacNAc analogs (LNs) 2−8 investigated in this study.

Figure 4 Spatial structure of the LN 2, showing the inter-residue O5′···H−O3 hydrogen bond and the torsions Φ2 and Ω1−Ω4.

Figure 5 Overlay of the crystal structures A) 2′F-LN 5 and 3′F-LN 6; B) α-d-LacNAc 1 and 3F-LN 3.

Figure 6 A) A summary of the changes of the torsions Φ, Ψ and Φ2 in the series 4−8, with the decreasing O5′···O3 distance. B)−D) Correlations between Φ, Ψ and Φ2 torsions and the O5′···O3 distance.

Figure 7 A)−C) Three possible staggered conformers, gg, gt and tg, of the exocyclic hydroxymethyl group in hexopyranoses according to the IUPAC rules.51 D)−E) The definition of two possible conformations, (Z)-anti and (Z)-syn, on the acetamido group.

Figure 8 A) The crystal structure of 2′F-LN 5 with a depiction of the F···H6proR/C6 distances. B) The crystal structure of 3F-LN 3 with a depiction of the F···H5′/C6′ distances.

Figure 9 A) The structures of the conformers A−L, represented by the LN 2. B) The mean-square deviations (MSDs) between calculated and experimentally obtained 13C NMR shifts for the conformers A−L of 2−8. The conformers H and I of 6′F-LN 8 represent the only exception from the series, as their 6′-hydroxymethyl conformation was set to gt. t.

Figure 10 Possible conformations around the glycosidic linkage, the exocyclic hydroxymethyl groups and the acetamido group, together with the corresponding expected exclusive NOE/ROE spatial contacts for every geometry.

Figure 11 A) The correlation between the O3−H temperature coefficients −Δδ/ΔT and the X-ray-determined O5′···O3 distance. B) The structure of LN 2 with a depiction of O5′···H−O3 and O6′−H···O3 hydrogen bonds.

Figure 12 A) The DFT-optimized geometry of 3′F-LN 6 showing the O2′−H···O6 hydrogen bond. B) The HSQC-TOCSY spectrum of 3′F-LN 6, depicting the O2′−H/C6 correlation.

Figure 13 A) The correlation between the experimentally determined and the calculated values of log P. B) The structure of the LN 2 with the colored labeling of individual hydroxyls illustrating the lipophilicity trends. C) The crystal structure of 6′F-LN 8 with a bridging water molecule.

Table 1 MM3 Molecular-Mechanics Conformational Analysis of the LacNAc Analogs 2−8

 	 	Distributions [%]b	
Compound	D90a	Φ/Ψg−/g+c	Φ/Ψg−/g−d	Φ/Ψg+/g+e	
LN 2	2.56%	98.99	0.07	0.22	
3F-LN 3	4.38%	92.45	1.16	1.23	
6F-LN 4	2.48%	98.84	0.09	0.32	
2′F-LN 5	2.41%	97.31	0.15	0.34	
3′F-LN 6	2.19%	88.87	0.56	9.38	
4′F-LN 7	2.70%	98.52	0.18	0.15	
6′F-LN 8	2.05%	97.15	0.19	0.62	
a D90 = relative Φ/Ψ conformational space corresponding to the 90% Boltzmann distribution of conformers.

b The distribution of defined Φ/Ψ: g−/g+, g−/g−and g+/g+ conformations.

c Φ/Ψ = (−110°− −30°)/(60°−160°).

d Φ/Ψ = (−110°− −30°)/(−90°− −30°).

e Φ/Ψ = (30°−90°)/(60°−160°).

Table 2 Comparison of O5′···O3 Distances and Torsions Φ, Ψ, Φ2, and Ω1−Ω4 in the Series of 3−8 with LacNAc 1h

a [O5′−C1′−O4−C4].

b [C1′−O4−C4−C3].

c [O5−C1−O−C].

d [O5−C5−C6−O6/F6].

e [O5′−C5′−C6′−O6′/′F6′].

f [C1−C2−N−C].

g [C2−N−C−C].

h The color code indicates the relative value of a given descriptor: blue = the value is higher in comparison with α-d-LacNAc 1,29 orange = the value is lower in comparison with α-d-LacNAc 1,29 green = a depiction of the trend in Φ2 cannot be compared to α-d-LacNAc 1(29) (α-configured hydroxyl group at C1).

Table 3 Experimental 3JH5−H6proR/H6proS, 3JH5′-H6′proR/H6′proS, and 3JH2-NH Couplings [Hz] and the Calculated Distribution f [%] of Side Chain Conformers gg/gt/tg18

 	GlcNAc	Gal	NHAc	
Compound	3JH5-H6proR	3JH5-H6proS	fgg	fgt	ftg	3JH5′-H6′proR	3JH5′-H6′proS	fgg	fgt	ftg	3JH2-NH	
LN 2	4.7a	2.2a	65	30	5	-b	-b	-b	-b	-b	8.6	
3F-LN 3	4.0c	2.4c	71	29	0	5.9d	7.1d	21	29	50	9.1	
6F-LN 4	4.6	1.8	64	35	1	6.7e	5.2e	23	54	23	8.5	
2′F-LN 5	5.1	1.9	65	35	0	6.8e	5.4e	20	53	27	8.6	
3′F-LN 6	5.0a	2.3a	63	34	4	7.3f	5.7f	11	56	33	8.6	
4′F-LN 7	4.8a	2.4a	65	30	5	7.3	5.9	12	51	37	8.6	
6′F-LN 8	4.8	2.2	66	32	1	7.9g	3.2g	17	76	7	8.6	
a Obtained using 1D selective homonuclear decoupling of O6-H.

b The couplings could not be determined due to a spectral overlap.

c Obtained using 1D selective NH−H5 TOCSY transfer.

d Obtained using 1D selective homonuclear decoupling of O6′-H.

e Estimated from deconvolution of H5′ using 1D selective H1′−H5′ ROESY transfer.

f Estimated from deconvolution of H5′ using 1D selective O4′H−H5′ TOCSY transfer.

g Confirmed by 1H measurement at 31 °C.

Table 4 Diagnostic ROESY Contactsa

contact	H1′-H4	H1′-H3/H5	H2′-H4	H4-H6proR	H4-H6proS	H4-H6′proR	H4-H6′proS	NH-H1/H3	NH-H2	NH-O3H	
conformation	Φ/Ψg−/g+	Φ/Ψ g−/g−	Φ/Ψ g+/g+	6-gt	6-tg	6′-gg	6′-gt	(Z)-anti	(Z)-syn	(Z)-antic	
LN 2	-b	-b	-b	-b	-b	-b	-b	S	W	W	
3F-LN 3	S	W	-b	-b	-b	VW	W	S	W	-d	
6F-LN 4	S	W	-	W	No	-b	-b	S	W	W	
2′F-LN 5	S	W	No	W	No	-b	-b	S	W	W	
3′F-LN 6	S	W	No	W	VW	-b	-b	S	W	W	
4′F-LN 7	S	W	-b	W	VW	VW	S	S	W	W	
6′F-LN 8	S	W	-b	M	W	No	S	S	W	W	
a Contact strength: S = strong, M = medium, W = weak, VW = very weak, No = contact not observed, - = analysis could not be performed.

b The detection of the contact was not possible due to the signal overlap.

c The contact indicates that the acetamide (Z)-anti geometry is slightly deviated from the ideal antiperiplanar arrangement.

d O3H is substituted by fluorine.

Table 5 Temperature Coefficients −Δδ/ΔT ppb/K of Exchangeable Protons in 2−8a

Compound	NHAc	OH-C3	OH-C2′	OH-C3′	OH-C4′	OH-C6′	OH-C6	
LN 2	4.2	2.0	6.8	7.2	-b	4.9	5.9	
3F-LN 3	4.0	-	6.0	6.8	5.8	5.7	5.3	
6F-LN 4	4.2	2.2	6.2	7.5	5.9	4.7	-	
2′F-LN 5	4.2	2.6	-	6.0	6.1	4.4	5.9	
3′F-LN 6	4.1	1.7	4.9	-	5.2	4.9	5.6	
4′F-LN 7	4.1	2.1	6.0	6.0	-	4.4	5.5	
6′F-LN 8	4.1	1.9	5.5	7.2	5.7	-	5.4	
a The temperature coefficients were determined in DMSO-d6, ppb = ppm·10−3.

b The analysis could not be performed due to the signal overlap.

Table 6 Lipophilicities of the LN Analogs 3−8a

Compound	3F-LN 3	6F-LN 4	2′F-LN 5	3′F-LN 6	4′F-LN 7	6′F-LN 8	
log Pexp.	−3.45	−2.95	−2.92	−3.09	−3.10	−3.28	
log Pcalc	−3.62	−2.84	−2.53	−3.16	−3.08	−3.07	
a The lipophilicity of 3−8 has been calculated from DFT-optimized global-minimum geometries using the COSMO-RS methodology.

Data Availability Statement

The data underlying this study are available in the published article and its Supporting Information

Supporting Information Available

The Supporting Information is available free of charge at https://pubs.acs.org/doi/10.1021/acs.joc.4c00879.The Supporting Information I contains outputs of molecular-mechanics calculations (energy and population diagrams); X-ray data of compounds 3−8; NMR assignment of compounds 3−8 in DMSO-d6; and DFT-calculated NMR parameters of compounds 3−8 (PDF).

The Supporting Information II contains copies of 1H, 13C, 19F, 1H−1H gCOSY, 1H−13C gHSQC, 1H−13C gHMBC, 1H−13C gHSQC-TOCSY, 1H−1H gROESY and selective 1D proton NMR-spectra for all investigated compounds, cartesian coordinates, computed total energy values, computed sum of electronic and thermal free energies, and number of imaginary frequencies of all optimized geometries (DOCX).

Supplementary Material

jo4c00879_si_001.pdf

jo4c00879_si_002.docx

The authors declare no competing financial interest.

Acknowledgments

The financial support from the Czech Science Foundation (grant No. 23-05146S) and from the Ministry of Education, Youth and Sports (grant No. LTC20052) is gratefully acknowledged. Martin Balouch would like to thank Karel Berka (UP Olomouc) for the access to the COSMOTherm software. Vojtech Hamala is grateful for the financial support of the Martina Roeselová Memorial Fellowship granted by the IOCB Tech Foundation. Computational resources have been provided by the e-INFRA CZ project (ID: 90254), supported by the Ministry of Education, Youth and Sports of the Czech Republic and by the ELIXIR-CZ project (ID: 90255), part of the international ELIXIR infrastructure. Jesús Jiménez-Barbero, Ana Ardá and F. Javier Cañada acknowledge the financial support from Spanish Ministry of Science, Innovation and Universities under grants PID2021-123781OB-C21 and C22.
==== Refs
References

Glycans (G02645) in IUPAC Compendium of Chemical Terminology, 3rd ed. International Union of Pure and Applied Chemistry, 2006. Online version 3.0.1, 201910.1351/goldbook.G02645.
Buffone A. Jr. ; Weaver V. M. Don’t Sugarcoat It: How Glycocalyx Composition Influences Cancer Progression. J. Cell Biol. 2020, 219 (1 ), e201910070 10.1083/jcb.201910070.31874115
Ambrosi M. ; Cameron N. R. ; Davis B. G. Lectins: Tools for the Molecular Understanding of the Glycocode. Org. Biomol. Chem. 2005, 3 (9 ), 1593–1608. 10.1039/b414350g.15858635
Leusmann S. ; Ménová P. ; Shanin E. ; Titz A. ; Rademacher C. Glycomimetics for the Inhibition and Modulation of Lectins. Chem. Soc. Rev. 2023, 52 (11 ), 3663–3740. 10.1039/D2CS00954D.37232696
Tamburrini A. ; Colombo C. ; Bernardi A. Design and Synthesis of Glycomimetics: Recent Advances. Med. Res. Rev. 2020, 40 (2 ), 495–531. 10.1002/med.21625.31361344
Hevey R. Bioisosteres of Carbohydrate Functional Groups in Glycomimetic Design. Biomimetics 2019, 4 (3 ), 53 10.3390/biomimetics4030053.31357673
Angulo J. ; Zimmer J. ; Imberty A. ; Prestegard J. H. Structural Biology of Glycan Recognition. In Essentials of Glycobiology, 4th ed.; Varki A. , Cummings R. D. , Esko J. D. , Stanley P. , Hart G. W. , Aebi M. , Mohnen D. , Kinoshita T. , Packer N. H. , Prestegard J. H. , Schnaar R. L. , Seeberger P. H. , Eds.; Cold Spring Harbor Laboratory Press: Cold Spring Harbor (NY), 2022; pp 403−418 10.1101/glycobiology.4e.30.
Verteramo M. L. ; Stenström O. ; Ignjatović M. M. ; Caldararu O. ; Olsson M. A. ; Manzoni F. ; Leffler H. ; Oksanen E. ; Logan D. T. ; Nilsson U. J. ; Ryde U. ; Akke M. Interplay between Conformational Entropy and Solvation Entropy in Protein−Ligand Binding. J. Am. Chem. Soc. 2019, 141 (5 ), 2012–2026. 10.1021/jacs.8b11099.30618244
Hevey R. Strategies for the Development of Glycomimetic Drug Candidates. Pharmaceuticals 2019, 12 (2 ), 55 10.3390/ph12020055.30978966
Yu Y. ; Delbianco M. Conformational Studies of Oligosaccharides. Chem.—Eur. J. 2020, 26 (44 ), 9814–9825. 10.1002/chem.202001370.32329095
Widmalm G. Glycan Shape, Motions, and Interactions Explored by NMR Spectroscopy. JACS Au 2024, 4 (1 ), 20–39. 10.1021/jacsau.3c00639.38274261
Dal Colle M. C. S. ; Fittolani G. ; Delbianco M. Synthetic Approaches to Break the Chemical Shift Degeneracy of Glycans. ChemBioChem. 2022, 23 (24 ), e202200416 10.1002/cbic.202200416.36005282
Hevey R. The Role of Fluorine in Glycomimetic Drug Design. Chem.—Eur. J. 2021, 27 , 2240–2253. 10.1002/chem.202003135.32901973
Linclau B. ; Ardá A. ; Reichardt N.-C. ; Sollogoub M. ; Unione L. ; Vincent S. P. ; Jiménez-Barbero J. Fluorinated Carbohydrates as Chemical Probes for Molecular Recognition Studies. Current Status and Perspectives. Chem. Soc. Rev. 2020, 49 (12 ), 3863–3888. 10.1039/C9CS00099B.32520059
St-Gelais J. ; Côté É. ; Lainé D. ; Johnson P. A. ; Giguère D. Addressing the Structural Complexity of Fluorinated Glucose Analogues: Insight into Lipophilicities and Solvation Effects. Chem.—Eur. J. 2020, 26 (59 ), 13499–13506. 10.1002/chem.202002825.32652740
Hamala V. ; Červenková Št’astná L. ; Kurfiřt M. ; Cuřínová P. ; Dračínský M. ; Karban J. Synthesis of Multiply Fluorinated N-Acetyl-D-Glucosamine and D-Galactosamine Analogs via the Corresponding Deoxyfluorinated Glucosazide and Galactosazide Phenyl Thioglycosides. Beilstein J. Org. Chem. 2021, 17 (1 ), 1086–1095. 10.3762/bjoc.17.85.34093878
Fernández P. ; Jiménez-Barbero J. The Conformation of Some Halodeoxy Analogues of Methyl β-Lactoside in D2O and DMSO-D6 Solutions. J. Carbohydr. Chem. 1994, 13 (2 ), 207–233. 10.1080/07328309408009189.
Bock K. ; Duus J. Ø. A Conformational Study of Hydroxymethyl Groups in Carbohydrates Investigated by 1H NMR Spectroscopy. J. Carbohydr. Chem. 1994, 13 (4 ), 513–543. 10.1080/07328309408011662.
Denavit V. ; Lainé D. ; Bouzriba C. ; Shanina E. ; Gillon É. ; Fortin S. ; Rademacher C. ; Imberty A. ; Giguère D. Stereoselective Synthesis of Fluorinated Galactopyranosides as Potential Molecular Probes for Galactophilic Proteins: Assessment of Monofluorogalactoside−LecA Interactions. Chem.—Eur. J. 2019, 25 , 4478–4490. 10.1002/chem.201806197.30690814
Stanley P. ; Moremen K. W. ; Lewis N. E. ; Taniguchi N. ; Aebi M. N-Glycans. In Essentials of Glycobiology; Varki A. , Cummings R. D. , Esko J. D. , Stanley P. , Hart G. W. , Aebi M. , Mohnen D. , Kinoshita T. , Packer N. H. , Prestegard J. H. , Schnaar R. L. , Seeberger P. H. , Eds.; Cold Spring Harbor Laboratory Press: Cold Spring Harbor (NY), 2022 10.1101/glycobiology.4e.9.
Brockhausen I. ; Wandall H. H. ; Hagen K. G. T. ; Stanley P. O-GalNAc Glycans. In Essentials of Glycobiology; Varki A. , Cummings R. D. , Esko J. D. , Stanley P. , Hart G. W. , Aebi M. , Mohnen D. , Kinoshita T. , Packer N. H. , Prestegard J. H. , Schnaar R. L. , Seeberger P. H. , Eds.; Cold Spring Harbor Laboratory Press: Cold Spring Harbor (NY), 2022 10.1101/glycobiology.4e.10.
Magalhães A. ; Ismail M. N. ; Reis C. A. Sweet Receptors Mediate the Adhesion of the Gastric Pathogen Helicobacter Pylori: Glycoproteomic Strategies. Expert Rev. Proteomics 2010, 7 (3 ), 307–310. 10.1586/epr.10.18.20536299
Mitchell E. ; Houles C. ; Sudakevitz D. ; Wimmerova M. ; Gautier C. ; Pérez S. ; Wu A. M. ; Gilboa-Garber N. ; Imberty A. Structural Basis for Oligosaccharide-Mediated Adhesion of Pseudomonas Aeruginosa in the Lungs of Cystic Fibrosis Patients. Nat. Struct. Mol. Biol. 2002, 9 (12 ), 918–921. 10.1038/nsb865.
Kerr C. L. ; Hanna W. F. ; Shaper J. H. ; Wright W. W. Lewis X-Containing Glycans Are Specific and Potent Competitive Inhibitors of the Binding of ZP3 to Complementary Sites on Capacitated, Acrosome-Intact Mouse Sperm1. Biol. Reprod. 2004, 71 (3 ), 770–777. 10.1095/biolreprod.103.023812.15128590
Bojar D. ; Meche L. ; Meng G. ; Eng W. ; Smith D. F. ; Cummings R. D. ; Mahal L. K. A Useful Guide to Lectin Binding: Machine-Learning Directed Annotation of 57 Unique Lectin Specificities. ACS Chem. Biol. 2022, 17 (11 ), 2993–3012. 10.1021/acschembio.1c00689.35084820
Cummings R. D. ; Liu F.-T. ; Rabinovich G. A. ; Stowell S. R. ; Vasta G. R. Galectins. In Essentials of Glycobiology; Varki A. , Cummings R. D. , Esko J. D. , Stanley P. , Hart G. W. , Aebi M. , Mohnen D. , Kinoshita T. , Packer N. H. , Prestegard J. H. , Schnaar R. L. , Seeberger P. H. , Eds.; Cold Spring Harbor Laboratory Press: Cold Spring Harbor (NY), 2022.10.1101/glycobiology.4e.36.
Arnold J. W. ; Whittington H. D. ; Dagher S. F. ; Roach J. ; Azcarate-Peril M. A. ; Bruno-Barcena J. M. Safety and Modulatory Effects of Humanized Galacto-Oligosaccharides on the Gut Microbiome. Front. Nutr. 2021, 8 , 640100 10.3389/fnut.2021.640100.33898497
Varki A. ; Cummings R. D. ; Aebi M. ; Packer N. H. ; Seeberger P. H. ; Esko J. D. ; Stanley P. ; Hart G. ; Darvill A. ; Kinoshita T. ; Prestegard J. J. ; Schnaar R. L. ; Freeze H. H. ; Marth J. D. ; Bertozzi C. R. ; Etzler M. E. ; Frank M. ; Vliegenthart J. F. ; Lütteke T. ; Perez S. ; Bolton E. ; Rudd P. ; Paulson J. ; Kanehisa M. ; Toukach P. ; Aoki-Kinoshita K. F. ; Dell A. ; Narimatsu H. ; York W. ; Taniguchi N. ; Kornfeld S. Symbol Nomenclature for Graphical Representations of Glycans. Glycobiology 2015, 25 (12 ), 1323–1324. 10.1093/glycob/cwv091.26543186
Longchambon F. ; Ohanessian J. ; Gillier-Pandraud H. ; Duchet D. ; Jacquinet J.-C. ; Sinaÿ P. Structure de la N-acétyllactosamine (acétamido-2 désoxy-2 O-β-d-galactopyrannosyl-4a-d-Glucopyrannose). Acta Cryst. Section B 1981, 37 (3 ), 601–607. 10.1107/S0567740881003658.
García-Aparicio V. ; Sollogoub M. ; Blériot Y. ; Colliou V. ; André S. ; Asensio J. L. ; Cañada F. J. ; Gabius H.-J. ; Sinaÿ P. ; Jiménez-Barbero J. The Conformation of the C-Glycosyl Analogue of N-Acetyl-Lactosamine in the Free State and Bound to a Toxic Plant Agglutinin and Human Adhesion/Growth-Regulatory Galectin-1. Carbohydr. Res. 2007, 342 (12 ), 1918–1928. 10.1016/j.carres.2007.02.034.17408600
Cocinero E. J. ; Gamblin D. P. ; Davis B. G. ; Simons J. P. The Building Blocks of Cellulose: The Intrinsic Conformational Structures of Cellobiose, Its Epimer, Lactose, and Their Singly Hydrated Complexes. J. Am. Chem. Soc. 2009, 131 (31 ), 11117–11123. 10.1021/ja903322w.19722675
Rivera-Sagredo A. ; Jiménez-Barbero J. ; Martín-Lomas M. N.m.r. Studies of the Conformation of Analogues of Methyl β-Lactoside in Methyl Sulfoxide-D6. Carbohydr. Res. 1991, 221 (1 ), 37–47. 10.1016/0008-6215(91)80047-Q.1816924
Erdélyi M. ; d’Auvergne E. ; Navarro-Vázquez A. ; Leonov A. ; Griesinger C. Dynamics of the Glycosidic Bond: Conformational Space of Lactose. Chem.—Eur. J. 2011, 17 , 9368–9376. 10.1002/chem.201100854.21755545
Schnupf U. ; Momany F. A. Rapidly Calculated DFT Relaxed Iso-Potential ϕ/ψ Maps: β-Cellobiose. Cellulose 2011, 18 (4 ), 859–887. 10.1007/s10570-011-9537-8.
Larsson E. A. ; Staaf M. ; Söderman P. ; Höög C. ; Widmalm G. Determination of the Conformational Flexibility of Methyl α-Cellobioside in Solution by NMR Spectroscopy and Molecular Simulations. J. Phys. Chem. A 2004, 108 (18 ), 3932–3937. 10.1021/jp049714c.
Tvaroška I. Theoretical Studies on the Conformation of Saccharides. VIII. Solvent Effect on the Stability of β-Cellobiose Conformers. Biopolymers 1984, 23 (10 ), 1951–1960. 10.1002/bip.360231011.
Yui T. ; Kobayashi H. ; Kitamura S. ; Imada K. Conformational Analysis of Chitobiose and Chitosan. Biopolymers 1994, 34 (2 ), 203–208. 10.1002/bip.360340206.
Mo F. ; Jensen L. H. The crystal structure of a β-(1→ 4) linked disaccharide, α-N,N’-diacetylchitobiose monohydrate. Acta Cryst. B 1978, 34 (5 ), 1562–1569. 10.1107/S0567740878006081.
Fettke A. ; Peikow D. ; Peter M. G. ; Kleinpeter E. Synthesis and Conformational Analysis of Glycomimetic Analogs of Thiochitobiose. Tetrahedron 2009, 65 (22 ), 4356–4366. 10.1016/j.tet.2009.03.067.
Mo F. On the Conformational Variability of the N-Acetylglucosamine β-(1→4) Linked Dimer. Crystal and Molecular Structure of Beta-N,N’-Diacetylchitobiose Trihydrate. Acta Chem. Scand. A 1979, 33 , 207–218. 10.3891/acta.chem.scand.33a-0207.
Martín-Pastor M. ; Canales A. ; Corzana F. ; Asensio J. L. ; Jiménez-Barbero J. Limited Flexibility of Lactose Detected from Residual Dipolar Couplings Using Molecular Dynamics Simulations and Steric Alignment Methods. J. Am. Chem. Soc. 2005, 127 (10 ), 3589–3595. 10.1021/ja043445m.15755180
Zhang W. ; Turney T. ; Meredith R. ; Pan Q. ; Sernau L. ; Wang X. ; Hu X. ; Woods R. J. ; Carmichael I. ; Serianni A. S. Conformational Populations of β-(1→4) O-Glycosidic Linkages Using Redundant NMR J-Couplings and Circular Statistics. J. Phys. Chem. B 2017, 121 (14 ), 3042–3058. 10.1021/acs.jpcb.7b02252.28296420
Bernet B. ; Vasella A. Intra- and Intermolecular H-Bonds of Alcohols in DMSO, 1H-NMR Analysis of Inter-Residue H-Bonds in Selected Oligosaccharides: Cellobiose, Lactose, N,N′-Diacetylchitobiose, Maltose, Sucrose, Agarose, and Hyaluronates. Helv. Chim. Acta 2000, 83 (9 ), 2055–2071. 10.1002/1522-2675(20000906)83:9<2055::AID-HLCA2055>3.0.CO;2-C.
Rönnols J. ; Engström O. ; Schnupf U. ; Säwén E. ; Brady J. W. ; Widmalm G. Inter-Residual Hydrogen Bonding in Carbohydrates Unraveled by NMR Spectroscopy and Molecular Dynamics Simulations. ChemBioChem. 2019, 20 (19 ), 2519–2528. 10.1002/cbic.201900301.31066963
Kurfiřt M. ; Dračínský M. ; Červenková Št’astná L. ; Cuřínová P. ; Hamala V. ; Hovorková M. ; Bojarová P. ; Karban J. Selectively Deoxyfluorinated N-Acetyllactosamine Analogues as 19F NMR Probes to Study Carbohydrate-Galectin Interactions. Chem.—Eur. J. 2021, 27 (51 ), 13040–13051. 10.1002/chem.202101752.34216419
Klamt A. ; Eckert F. COSMO-RS: A Novel and Efficient Method for the a Priori Prediction of Thermophysical Data of Liquids. Fluid Phase Equilib. 2000, 172 (1 ), 43–72. 10.1016/S0378-3812(00)00357-5.
Bannan C. C. ; Burley K. H. ; Chiu M. ; Shirts M. R. ; Gilson M. K. ; Mobley D. L. Blind Prediction of Cyclohexane-Water Distribution Coefficients from the SAMPL5 Challenge. J. Comput. Aided Mol. Des. 2016, 30 (11 ), 927–944. 10.1007/s10822-016-9954-8.27677750
Vidal P. ; Jiménez-Barbero J. ; Espinosa J. F. Conformational Flexibility around the Gal-β-(1 → 3)-Glc Linkage: Experimental Evidence for the Existence of the Anti-ψ Conformation in Aqueous Solution. Carbohydr. Res. 2016, 433 , 36–40. 10.1016/j.carres.2016.06.009.27434833
O’Dell W. B. ; Baker D. C. ; McLain S. E. Structural Evidence for Inter-Residue Hydrogen Bonding Observed for Cellobiose in Aqueous Solution. PLoS One 2012, 7 (10 ), e45311 10.1371/journal.pone.0045311.23056199
Fowler P. ; Bernet B. ; Vasella A. A 1H-NMR Spectroscopic Investigation of the Conformation of the Acetamido Group in Some Derivatives of N-Acetyl-d-Allosamine and -d-Glucosamine. Helv. Chim. Acta 1996, 79 (1 ), 269–287. 10.1002/hlca.19960790127.
Grindley T. B. Structure and Conformation of Carbohydrates. In Glycoscience: Chemistry and Chemical Biology; Fraser-Reid B. O. , Tatsuta K. , Thiem J. , Eds.; Springer: Berlin, Heidelberg, 2008; pp 3−55 10.1007/978-3-540-30429-6_1.
Battistel M. D. ; Azurmendi H. F. ; Yu B. ; Freedberg D. I. NMR of Glycans: Shedding New Light on Old Problems. Prog. Nucl. Magn. Reson. Spectrosc. 2014, 79 , 48–68. 10.1016/j.pnmrs.2014.01.001.24815364
Csuk R. ; Glänzer B. I. N.M.R. Spectroscopy of Fluorinated Monosaccharides. In Adv. Carbohydr. Chem. Biochem.; Tipson R. S. , Horton D. , Eds.; Academic Press, 1988; Vol. 46 , pp 73−177 10.1016/S0065-2318(08)60165-1.
Michalik M. ; Hein M. ; Frank M. NMR Spectra of Fluorinated Carbohydrates. Carbohydr. Res. 2000, 327 (1 ), 185–218. 10.1016/S0008-6215(99)00323-7.10968682
Chaloupecká E. ; Kurfiřt M. ; Červenková Št’astná L. ; Karban J. ; Dračínský M. Exploring Long-Range Fluorine−Carbon J-Coupling for Conformational Analysis of Deoxyfluorinated Disaccharides: A Combined Computational and NMR Study. Bioorg. Chem. 2024, 147 , 107388 10.1016/j.bioorg.2024.107388.38678775
Bork N. ; Du L. ; Reiman H. ; Kurtén T. ; Kjaergaard H. G. Benchmarking Ab Initio Binding Energies of Hydrogen-Bonded Molecular Clusters Based on FTIR Spectroscopy. J. Phys. Chem. A 2014, 118 (28 ), 5316–5322. 10.1021/jp5037537.24988143
Oswald S. ; Suhm M. A. Soft Experimental Constraints for Soft Interactions: A Spectroscopic Benchmark Data Set for Weak and Strong Hydrogen Bonds. Phys. Chem. Chem. Phys. 2019, 21 (35 ), 18799–18810. 10.1039/C9CP03651B.31453998
Sergeyev I. ; Moyna G. Determination of the Three-Dimensional Structure of Oligosaccharides in the Solid State from Experimental 13C NMR Data and Ab Initio Chemical Shift Surfaces. Carbohydr. Res. 2005, 340 (6 ), 1165–1174. 10.1016/j.carres.2005.02.022.15797132
Swalina C. W. ; Zauhar R. J. ; DeGrazia M. J. ; Moyna G. Derivation of 13C Chemical Shift Surfaces for the Anomeric Carbons of Oligosaccharides and Glycopeptides Using Ab Initio Methodology. J. Biomol. NMR 2001, 21 (1 ), 49–61. 10.1023/A:1011928919734.11693568
Krivdin L. B. Computational NMR of Carbohydrates: Theoretical Background, Applications, and Perspectives. Molecules 2021, 26 (9 ), 2450 10.3390/molecules26092450.33922318
Dračínský M. ; Bouř P. Computational Analysis of Solvent Effects in NMR Spectroscopy. J. Chem. Theory Comput. 2010, 6 (1 ), 288–299. 10.1021/ct900498b.26614339
Poppe L. ; Stuike-Prill R. ; Meyer B. ; van Halbeek H. The Solution Conformation of Sialyl-α(2→6)-Lactose Studied by Modern NMR Techniques and Monte Carlo Simulations. J. Biomol. NMR 1992, 2 (2 ), 109–136. 10.1007/BF01875524.1422148
Charisiadis P. ; Kontogianni V. G. ; Tsiafoulis C. G. ; Tzakos A. G. ; Siskos M. ; Gerothanassis I. P. 1H-NMR as a Structural and Analytical Tool of Intra- and Intermolecular Hydrogen Bonds of Phenol-Containing Natural Products and Model Compounds. Molecules 2014, 19 (9 ), 13643–13682. 10.3390/molecules190913643.25185070
Klein R. A. Electron Density Topological Analysis of Hydrogen Bonding in Glucopyranose and Hydrated Glucopyranose. J. Am. Chem. Soc. 2002, 124 (46 ), 13931–13937. 10.1021/ja0206947.12431125
Klein R. A. Hydrogen bonding in diols and binary diol−water systems investigated using DFT methods. II. Calculated infrared OH-stretch frequencies, force constants, and NMR chemical shifts correlate with hydrogen bond geometry and electron density topology. A reevaluation of geometrical criteria for hydrogen bonding. J. Comput. Chem. 2003, 24 (9 ), 1120–1131. 10.1002/jcc.10256.12759911
Oh K.-I. ; Rajesh K. ; Stanton J. F. ; Baiz C. R. Quantifying Hydrogen-Bond Populations in Dimethyl Sulfoxide/Water Mixtures. Angew. Chem., Int. Ed. 2017, 56 (38 ), 11375–11379. 10.1002/anie.201704162.
Linclau B. ; Graton J. ; Le Questel J.-Y. 8 - Influence of Fluorination on Alcohol Hydrogen-Bond Donating Properties. In Progress in Fluorine Science, Fluorine in Life Sciences: Pharmaceuticals, Medicinal Diagnostics, and Agrochemicals; Haufe G. , Leroux F. R. , Eds.; Progress in Fluorine Science; Academic Press, 2019; pp 301−324 10.1016/B978-0-12-812733-9.00008-8.
Graton J. ; Compain G. ; Besseau F. ; Bogdan E. ; Watts J. M. ; Mtashobya L. ; Wang Z. ; Weymouth-Wilson A. ; Galland N. ; Le Questel J.-Y. ; Linclau B. Influence of Alcohol β-Fluorination on Hydrogen-Bond Acidity of Conformationally Flexible Substrates. Chem.—Eur. J. 2017, 23 (12 ), 2811–2819. 10.1002/chem.201604940.27906491
Pietruś W. ; Kurczab R. ; Kafel R. ; Machalska E. ; Kalinowska-Tłuścik J. ; Hogendorf A. ; Żylewski M. ; Baranska M. ; Bojarski A. J. How Can Fluorine Directly and Indirectly Affect the Hydrogen Bonding in Molecular Systems? − A Case Study for Monofluoroanilines. Spectrochim. Acta, Part A 2021, 252 , 119536 10.1016/j.saa.2021.119536.
Linclau B. ; Wang Z. ; Compain G. ; Paumelle V. ; Fontenelle C. Q. ; Wells N. ; Weymouth-Wilson A. Investigating the Influence of (Deoxy)Fluorination on the Lipophilicity of Non-UV-Active Fluorinated Alkanols and Carbohydrates by a New Log P Determination Method. Angew. Chem., Int. Ed. 2016, 55 (2 ), 674–678. 10.1002/anie.201509460.
Wang D. ; Ámundadóttir M. L. ; van Gunsteren W. F. ; Hünenberger P. H. Intramolecular Hydrogen-Bonding in Aqueous Carbohydrates as a Cause or Consequence of Conformational Preferences: A Molecular Dynamics Study of Cellobiose Stereoisomers. Eur. Biophys. J. 2013, 42 (7 ), 521–537. 10.1007/s00249-013-0901-5.23660988
Mohamadi F. ; Richards N. G. J. ; Guida W. C. ; Liskamp R. ; Lipton M. ; Caufield C. ; Chang G. ; Hendrickson T. ; Still W. C. Macromodel—an Integrated Software System for Modeling Organic and Bioorganic Molecules Using Molecular Mechanics. J. Comput. Chem. 1990, 11 (4 ), 440–467. 10.1002/jcc.540110405.
Watts K. S. ; Dalal P. ; Tebben A. J. ; Cheney D. L. ; Shelley J. C. Macrocycle Conformational Sampling with MacroModel. J. Chem. Inf. Model. 2014, 54 (10 ), 2680–2696. 10.1021/ci5001696.25233464
Allinger N. L. ; Yuh Y. H. ; Lii J. H. Molecular mechanics. The MM3 force field for hydrocarbons. 1. J. Am. Chem. Soc. 1989, 111 (23 ), 8551–8566. 10.1021/ja00205a001.
Onufriev A. V. ; Case D. A. Generalized Born Implicit Solvent Models for Biomolecules. Annu. Rev. Biophys. 2019, 48 , 275–296. 10.1146/annurev-biophys-052118-115325.30857399
Polak E. ; Ribiere G. Note sur la convergence de méthodes de directions conjuguées. ESAIM: Math. Modell. Numer. Anal. 1969, 3 (R1 ), 35–43. 10.1051/m2an/196903R100351.
Eaton J. W. ; Bateman D. ; Hauberg S. ; Wehbring R. GNU Octave Version 8.3.0 Manual: A High-Level Interactive Language for Numerical Computations; 2023. https://www.gnu.org/software/octave/doc/v8.3.0/ (accessed 2024-04-10).
Mittal A. ; Malhotra D. ; Jain P. ; Kalia A. ; Shunmugaperumal T. Studies on Aspirin Crystals Generated by a Modified Vapor Diffusion Method. AAPS PharmSciTech 2016, 17 (4 ), 988–994. 10.1208/s12249-015-0432-8.26729528
Sheldrick G. M. SHELXL-2017/1, Program for the Solution of Crystal Structures; University of Göttingen, Germany.
Cooper R. I. ; Thompson A. L. ; Watkin D. J. CRYSTALS Enhancements: Dealing with Hydrogen Atoms in Refinement. J. Appl. Crystallogr. 2010, 43 (5 ), 1100–1107. 10.1107/S0021889810025598.
Farrugia L. J. ORTEP-3 for Windows - a Version of ORTEP-III with a Graphical User Interface (GUI). J. Appl. Crystallogr. 1997, 30 (5 ), 565–565. 10.1107/S0021889897003117.
Frisch M. J. ; Trucks G. W. ; Schlegel H. B. ; Scuseria G. E. ; Robb M. A. ; Cheeseman J. R. ; Scalmani G. ; Barone V. ; Petersson G. A. ; Nakatsuji H. ; Li X. ; Caricato M. ; Marenich A. V. ; Bloino J. ; Janesko B. G. ; Gomperts R. ; Mennucci B. ; Hratchian H. P. ; Ortiz J. V. ; Izmaylov A. F. ; Sonnenberg J. L. ; Williams-Young D. , Ding F. ; Lipparini F. ; Egidi F. ; Goings J. ; Peng B. ; Petrone A. ; Henderson T. ; Ranasinghe D. ; Zakrzewski V. G. ; Gao J. ; Rega N. ; Zheng G. ; Liang W. ; Hada M. ; Ehara M. ; Toyota K. ; Fukuda R. ; Hasegawa J. ; Ishida M. ; Nakajima T. ; Honda Y. ; Kitao O. ; Nakai H. ; Vreven T. ; Throssell K. ; Montgomery J. A. Jr. ; Peralta J. E. ; Ogliaro F. ; Bearpark M. J. ; Heyd J. J. ; Brothers E. N. ; Kudin K. N. ; Staroverov V. N. ; Keith T. A. ; Kobayashi R. ; Normand J. ; Raghavachari K. ; Rendell A. P. ; Burant J. C. ; Iyengar S. S. ; Tomasi J. ; Cossi M. ; Millam J. M. ; Klene M. ; Adamo C. ; Cammi R. ; Ochterski J. W. ; Martin R. L. ; Morokuma K. ; Farkas O. ; Foresman J. B. ; Fox D. J. Gaussian 16, Revision B.01; Gaussian, Inc., Wallingford, CT, 2016.
Lee C. ; Yang W. ; Parr R. G. Development of the Colle-Salvetti Correlation-Energy Formula into a Functional of the Electron Density. Phys. Rev. B 1988, 37 (2 ), 785–789. 10.1103/PhysRevB.37.785.
Becke A. D. Density-functional Thermochemistry. III. The Role of Exact Exchange. J. Chem. Phys. 1993, 98 (7 ), 5648–5652. 10.1063/1.464913.
Vosko S. H. ; Wilk L. ; Nusair M. Accurate Spin-Dependent Electron Liquid Correlation Energies for Local Spin Density Calculations: A Critical Analysis. Can. J. Phys. 1980, 58 (8 ), 1200–1211. 10.1139/p80-159.
Stephens P. J. ; Devlin F. J. ; Chabalowski C. F. ; Frisch M. J. Ab Initio Calculation of Vibrational Absorption and Circular Dichroism Spectra Using Density Functional Force Fields. J. Phys. Chem. 1994, 98 (45 ), 11623–11627. 10.1021/j100096a001.
Grimme S. ; Antony J. ; Ehrlich S. ; Krieg H. A Consistent and Accurate Ab Initio Parametrization of Density Functional Dispersion Correction (DFT-D) for the 94 Elements H-Pu. J. Chem. Phys. 2010, 132 (15 ), 154104 10.1063/1.3382344.20423165
Grimme S. ; Hansen A. ; Brandenburg J. G. ; Bannwarth C. Dispersion-Corrected Mean-Field Electronic Structure Methods. Chem. Rev. 2016, 116 (9 ), 5105–5154. 10.1021/acs.chemrev.5b00533.27077966
Barone V. ; Cossi M. Quantum Calculation of Molecular Energies and Energy Gradients in Solution by a Conductor Solvent Model. J. Phys. Chem. A 1998, 102 (11 ), 1995–2001. 10.1021/jp9716997.
Cossi M. ; Rega N. ; Scalmani G. ; Barone V. Energies, Structures, and Electronic Properties of Molecules in Solution with the C-PCM Solvation Model. J. Comput. Chem. 2003, 24 (6 ), 669–681. 10.1002/jcc.10189.12666158
London F. Théorie quantique des courants interatomiques dans les combinaisons aromatiques. J. Phys. Radium 1937, 8 (10 ), 397–409. 10.1051/jphysrad:01937008010039700.
McWeeny R. Perturbation Theory for the Fock-Dirac Density Matrix. Phys. Rev. 1962, 126 (3 ), 1028–1034. 10.1103/PhysRev.126.1028.
Ditchfield R. Self-Consistent Perturbation Theory of Diamagnetism. Mol. Phys. 1974, 27 (4 ), 789–807. 10.1080/00268977400100711.
Wolinski K. ; Hinton J. F. ; Pulay P. Efficient Implementation of the Gauge-Independent Atomic Orbital Method for NMR Chemical Shift Calculations. J. Am. Chem. Soc. 1990, 112 (23 ), 8251–8260. 10.1021/ja00179a005.
Cheeseman J. R. ; Trucks G. W. ; Keith T. A. ; Frisch M. J. A Comparison of Models for Calculating Nuclear Magnetic Resonance Shielding Tensors. J. Chem. Phys. 1996, 104 (14 ), 5497–5509. 10.1063/1.471789.
Deng W. ; Cheeseman J. R. ; Frisch M. J. Calculation of Nuclear Spin−Spin Coupling Constants of Molecules with First and Second Row Atoms in Study of Basis Set Dependence. J. Chem. Theory Comput. 2006, 2 (4 ), 1028–1037. 10.1021/ct600110u.26633062
Lodewyk M. W. ; Siebert M. R. ; Tantillo D. J. Computational Prediction of 1H and 13C Chemical Shifts: A Useful Tool for Natural Product, Mechanistic, and Synthetic Organic Chemistry. Chem. Rev. 2012, 112 (3 ), 1839–1862. 10.1021/cr200106v.22091891
Bally T. ; Rablen P. R. Quantum-Chemical Simulation of 1H NMR Spectra. 2. Comparison of DFT-Based Procedures for Computing Proton−Proton Coupling Constants in Organic Molecules. J. Org. Chem. 2011, 76 (12 ), 4818–4830. 10.1021/jo200513q.21574622
Weigend F. ; Ahlrichs R. Balanced Basis Sets of Split Valence, Triple Zeta Valence and Quadruple Zeta Valence Quality for H to Rn: Design and Assessment of Accuracy. Phys. Chem. Chem. Phys. 2005, 7 (18 ), 3297–3305. 10.1039/b508541a.16240044
TURBOMOLE 6.3, a Development of University of Karlsruhe and Forschungszentrum Karlsruhe GmbH, 1989−2007; TURBOMOLE GmbH, 2011.
