1
Gas-phase structure and conformation of the glycosidase and ceramide glucosyltransferase inhibitor N-benzyl deoxynojirimycin

2
The deoxynojirimycin (DNJ) family of imino sugars are glucose analogues with an NH group replacing the oxygen atom in the pyranose ring.

3
They are powerful inhibitors of glycosidases and ceramide glucosyltransferases.

4
The conformation of N-benzyl-DNJ, isolated in the gas phase, is studied using a combination of resonant two-photon ionization (R2PI), UV/UV hole-burn, and IR ion-dip spectroscopies in conjunction with electronic structure theory calculations.

5
Three distinct conformers, one major and two minor, are present and all three are assigned to structures in which the exocyclic hydroxymethyl group is axial to the piperidine ring (gauche- to the ring nitrogen).

6
This contrasts with the preferred equatorial (gauche+) orientation observed for simple glucosides and may well contribute to the stronger binding of some enzymes to DNJ-based inhibitors compared to their natural glucoside substrates.

Introduction

7
N-Benzyl deoxynojirimycin (N-Bn-DNJ) 1, the N-benzylated derivative of deoxynojirimycin (DNJ) 2, is a mimic of glucose 3 (Scheme 1).

8
DNJ, a member of an important family of imino sugars that function as inhibitors of glycosidases, is an inhibitor of α- and β-glucosidases.1,2

9
Enzymes other than glycosidases are also inhibited by imino sugars; for example, derivatives of DNJ are powerful inhibitors of ceramide-specific glucosyltransferase3,4 whereas DNJ itself is not an effective inhibitor of this enzyme.5

10
X-ray diffraction and NMR studies of DNJ and its N-alkyl derivatives6,7 have revealed differing conformational preferences.

11
In DNJ the hydroxymethyl O(6)H group adopts an equatorial orientation with respect to the piperidine ring, possibly favoured by a weak, intramolecular NH → O6 hydrogen-bond [r(N⋯O) = 2.86 Å, ∠NH–O = 104°]6 but in the N-alkyl derivatives, where this H-bonding is no longer possible, its orientation switches to axial.7

12
It has been suggested that this axial orientation is an important structural requirement for potent α-glucosidase inhibition and the consequent interference in the biosynthesis of N-linked glycoproteins (and anti-viral activity).8

13
The enhanced potency of some α-glucosidases by castanospermine 4, where the O(6)H is forced into an axial conformation by the fused five-membered ring, has lent support to this view.6

14
Despite the potential therapeutic uses of imino sugars, only two crystal structures of biological complexes containing DNJ as a tightly-bound inhibitor have been reported, one to an α- and one to a β-glucosidase.9,10

15
In both, the O(6)H interacts with an acidic residue and is oriented axially to the piperidine ring, which is in a (distorted) 4C1 chair conformation.

16
Interestingly, both also show ordered water molecules within the active site interacting with the ring nitrogen.

17
N-alkylation could serve to displace the bound water molecules, representing an entropic advantage for binding these inhibitors (an advantage also suggested for the increased potency of isofagomine115 over DNJ by Zechel and coworkers10); no complexes containing bound N-alkylated DNJ are available in the protein data bank at this time.

18
The only available structure containing bound castanospermine shows the piperidine ring distorted into a twisted boat conformation which is suggested to be a precursor to the transition state for the cleavage of the glucoside.12

19
Molecular mechanics (MM3) and low-level ab initio (HF/6-31G) calculations for DNJ and castanospermine have supported the assignment of the two structures (established by diffraction and NMR measurements in the condensed phase) to the global minimum conformations predicted for the two molecules isolated in the gas phase.6

20
To date however, there have been no experimental determinations of their conformational landscapes under isolated conditions free of environmental perturbations to compare with (higher-level) theoretical results, or of the ease of populating alternative conformers.

21
Our recent investigations of the conformational structures of the model sugars O-phenyl β-d-gluco- or d-galacto-pyranoside, vaporised by heating or laser ablation and cooled in a free-jet argon expansion, probed via UV mass-selected resonant two-photon ionization (R2PI), resolved through ‘hole-burning’ techniques, and assigned through analysis of their resonant, mid-IR ion-dip spectra,13,14 have prompted their extension to enable the first investigation of an imino sugar in the gas phase.

Methods

Spectroscopy

22
Gas-phase N-Bn-DNJ (synthetic details provided as electronic supplementary information (ESI)) was produced by placing the sample in an oven attached to the expansion side of a pulsed nozzle (General Valve, 0.8 mm diameter, 4 bar argon backing pressure, 10 Hz pulse rate).

23
The molecules, vaporized at 175 °C, were thus entrained and cooled in the pulsed expanding argon jet.

24
The beam was skimmed and then intersected by infrared (IR) and/or ultraviolet (UV) laser beam(s) tuned to the desired wavelength(s) for excitation and ionization in the acceleration region of a time-of-flight mass spectrometer (RM Jordan) [using a resonant two-photon ionization (R2PI) scheme] enabling mass selective detection of the ions thus produced.

25
Tuneable UV laser radiation for the R2PI experiments was generated by frequency doubling the output of a Nd:YAG-pumped pulsed dye laser [GCR-170/Laser Analytical Systems (LAS)] operated with a mixture of Coumarin 153 and Coumarin 500 dye (∼0.4–1 mJ pulse−1).

26
This laser was employed as the scanned, ‘burn’ laser in UV/UV hole-burn experiments: the fixed ‘probe’ pulse, fired ∼200 ns after the burn laser, was generated by a second, excimer-pumped pulsed dye laser (Lambda Physik FL3002), tuned to successive features of interest in the R2PI spectrum.

27
Resonant IR ion dip scheme spectra were recorded in an analogous manner with the burn laser generating tuneable pulsed IR radiation (∼3000–4000 cm−1, 0.4 cm−1 bandwidth, 4.5 mJ pulse−1), through difference frequency mixing the output of a Nd:YAG laser (Powerlite Precision 8000) and that of a pumped dye laser (ND6000, Continuum) in a LiNbO3 crystal: the UV probe pulse was provided by the LAS laser system, fired ∼150 ns after the IR pulse.

Computation

28
Starting structures for geometry optimization were generated using molecular mechanics (MMFF94c and MM3 forcefields) and conformation searching with the Monte Carlo multiple minimization method implemented in the MacroModel software (MacroModel v.8.5, Schrödinger, LLC);15 1000 structures were generated.

29
A total of 40 conformers were selected for subsequent electronic structure theory calculations; these included all structures within 15 kJ mol−1 of the most stable in the conformational search, as well as other assorted structures, such as ones with the benzyl group axial to the (4C1) piperidine ring, 1C4 conformers and boat conformers.

30
The selected conformers were further optimized using B3LYP hybrid-method density functional theory and the 6-31+G(d) basis set.

31
Analytical frequency and thermochemistry calculations were performed on the B3LYP geometry-optimized conformations, as were single-point energy calculations, using perturbation theory (MP2) and the 6-311+G(d,p) basis set.

32
All electronic structure theory calculations were performed using the GAUSSIAN software16.

Spectroscopy and structure

33
The lowest lying conformers of N-Bn-DNJ computed ab initio are shown in Fig. 1 The global minimum structure is calculated to lie at an energy 4.5 kJ mol−1 below its nearest neighbour, see Table 1.

34
All the conformers are associated with a gauche (G−) [O6–C6–C5–N] conformation, creating an axially oriented O(6)H group (regardless of the orientation of the N-benzyl substituent): the lowest lying equatorially oriented structures, gauche (G+) or trans (T), all lie at relative energies ≥9 kJ mol−1 above the global minimum.

35
The six most stable conformers differ primarily by the orientation of the benzyl ring (I vs. II vs. III) and the orientation of the co-operative ‘daisy chain’ sequence of hydrogen bonds formed by the hydroxyls on the ring carbons [(a) vs. (b): in conformers (a) they adopt a gauche orientation with respect to the Cn+1 carbon, while conformers of type (b) present a trans orientation].

36
Conformers IIa,b and IIIa,b are stabilized by a hydrogen bond between OH6 and the ring nitrogen [r(OH6⋯N) ∼ 2.2–2.4 Å, ∠OH–N ∼115°], an interaction which is considerably weaker in the conformers Ia,b [r(OH6⋯N) ∼ 2.8–3.0 Å, ∠OH–N ∼95°] because of the position of the benzyl ring; on the other hand conformer Ia, the global minimum structure, displays the longest ‘daisy chain’ sequence of hydroxyl H-bonds.

37
Fig. 2 shows the one-colour R2PI spectrum of jet-cooled N-Bn-DNJ (top trace) and below it three distinct UV hole-burn spectra, labeled α, β and β′; they account for all of the intense peaks observed in the R2PI spectrum and correspond to three different conformers of the molecule.

38
The UV contour and vibronic spacings displayed by the principal conformer, α, are quite distinct from the (nearly identical) spectra displayed by the two minor conformers, β and β′.

39
In addition to their S1 ← S0 origin bands, centred around 37 280 cm−1 (α) and 37 500 cm−1 (β/β′), all three hole-burn spectra display a series of more intense vibronic bands located ∼530 cm−1 higher in energy (associated with the excitation of a benzyl ring mode, a feature commonly seen in one colour R2PI spectra of compounds containing this chromophore17).

40
Each band accommodates a series of much lower frequency vibronic progressions with spacings of ∼37.5 and 41 cm−1 (conformer α) and ∼21 and 24 cm−1 (conformers β and β′).

41
These spectra suggest that conformers β and β′ are extremely similar to each other while conformer α is somewhat different.

42
The mass of N-Bn-DNJ is 253 u but in all three conformers more than 90% of the observed ion signal is at mass 222, a loss of 31 u.

43
Spectra measured in the fragment and parent ion mass channels show the same features, though conformers β and β′ have a greater ratio of fragment to parent ions than conformer α.

44
The fragmentation, in all probability, corresponds to the cleavage of the C5–C6 bond and loss of the hydroxymethyl radical.

45
Although the benzyl ring is ionized in the R2PI scheme used, the ionization energy of tertiary amines is ∼1 eV lower than benzyl compounds18 and migration of the charge from the benzyl ring chromophore to the amine will endow the molecular ion with significant internal energy and promote its fragmentation, further promoted by excitation of the ring deformation mode at ∼530 cm−1.

46
The difference in fragmentation ratios among the conformers may simply be due to the larger excess energy in the vibrationally excited ions, but it may also be due to differences in Franck–Condon overlap between the excited and ionized states in the different conformers α and β/β′.

47
Ion-dip spectra of the three conformers, measured in the O–H stretching region of the infrared are shown in Fig. 3.

48
In each case all four O–H modes are visible although the central pair lie very close in energy and therefore appear as a partially resolved doublet.

49
The close similarity between the minor conformers, β and β′, is again apparent; each one presents a very similar vibrational spectrum with the O–H band at lowest wavenumber more strongly shifted than for conformer α.

50
The O–H stretching modes are very sensitive to the local environment and provide good indicators of conformation.

51
In the monosaccharides13,14O-phenyl β-d-gluco- and d-galacto-pyranoside, vibrational bands appearing in the region 3630–3650 cm−1 are associated with a sequence of weakly bonded hydroxyl groups arranged in an equatorial fashion around the pyranoside sugar ring.

52
Stronger H-bonded interactions shift the O–H bands to lower wavenumber and broaden them; those associated with relatively ‘free’ (not H-bond donating) OH groups appear at slightly higher wave number.

53
All three conformers of N-Bn-DNJ exhibit one relatively free O–H mode, two bands consistent with a weakly hydrogen bonded and one band associated with a more strongly bonded OH group.

54
In conformers β and β′, particularly, this latter band is considerably shifted (∼80 cm−1 below the central doublet) indicating a strong perturbation of one of the OH stretches in both these conformers.

55
As the hydroxyl groups attached directly to the ring (in its 4C1 chair conformation) are geometrically constrained and cannot form strong interactions with each other, this implies a strong interaction between the hydroxymethyl group and most likely, the ring nitrogen in the two conformers, β and β′.

56
Infrared spectra calculated for each of the six lowest-lying conformers (see Fig. 1) are also shown in Fig. 3.

57
Conformers Ia and Ib, differing in the orientation of the ring hydroxyl groups (note the effect on the relative positions of the alternately ‘free’ or hydrogen-bonded O–H stretch modes, σ4 and σ2, in each one of the three ‘a,b’ conformer pairs) and conformer IIa all present very similar IR spectra and each could provide a reasonable match for the recorded spectrum of the principal conformer α.

58
On balance however, its assignment to the structure Ia is strongly favoured for the following reasons:

59
(i) Calculations at 0 K place conformer Ia at the global minimum, lying 4.5 kJ mol−1 below the energy of conformer IIa and 7.5 kJ mol−1 below Ib.

60
(ii) There is a very close correspondence between the spacings in the vibronic progressions in its R2PI spectrum, ∼37.5 and 41.0 cm−1 (in the S1 state) and the in-plane and out-of-plane benzyl ring wag modes of conformer Ia calculated (in the S0 state) at 37.6 and 42.3 cm−1.

61
The equivalent lowest wavenumber modes in conformer IIa are calculated at 20.6 and 24.1 cm−1.

62
The assignment of the two minor conformers, β and β′, is less straightforward; the strong shift of the lowest wavenumber OH band provides an important clue but each one of the conformer pairs IIa/b and IIIa/b exhibits a hydrogen bond between the hydroxymethyl group and the ring nitrogen which significantly lowers the calculated wave number of the O–H6 band, σ6.

63
Nonetheless, one of the remaining conformers, β or β′, is most likely to be associated with the conformation IIb, since this is the second most stable conformer at 0 K and its calculated IR spectrum provides an excellent match to the measured IR spectra of both β and β′.

64
Assignment of the partner conformer to the structure IIIb (which differs from IIb only in the orientation of the benzyl group) is tempting but at 0 K it is calculated to lie ∼4 kJ mol−1 higher in energy than conformer IIa [and at the temperature of the oven IIa becomes the lowest energy conformer while IIIb is the highest of the six (Table 1)].

65
Assignment of conformers β/β′ to the IIa/IIb pair is favoured by their similar low frequency vibration/internal rotation modes (calculated to be a twist of the benzyl group) lying at 24.1 and 20.6 cm−1, respectively (in the ground state, S0); these match very well the measured (S1) vibronic progressions of 24 and 21 cm−1.

66
In contrast, normal mode analysis of the conformers IIIa/b show the lowest energy mode (also a benzyl twist though with some benzyl wag coupled in) at 34.5 and 35.3 cm−1.

67
Additional evidence favouring assignment to the IIa/IIb pair is provided by the extreme similarity of their R2PI spectra.

68
The environment of the benzyl ring chromophore of conformer IIb is closer to that of conformer IIa (which differs in the orientation of the chain of hydroxyl groups) than conformer IIIb (which differs in the orientation of the benzyl ring).

69
On balance, therefore the assignment of β and β′ to conformers IIa/IIb is the most favoured.

70
Most importantly, both the major (α) and the two minor (β, β′) conformers are associated with structures in which the OH6 group is axial (G−) to the ring nitrogen.

71
The conserved G− orientation of the hydroxymethyl group of N-Bn-DNJ differs from the preferred orientation in simple glucosides.

72
In 1-O-phenyl β-d-gluco-pyranoside, only ∼25% of the observed population adopts a G− orientation; the most populated conformation (∼70%) is G+.13

73
In the corresponding galacto-pyranoside, the preference for G+ is even stronger (∼90%).14

74
This innate preference for the G− orientation in the imino sugars may well contribute to their strong binding compared to the natural glucoside substrate.

75
Although the present investigation has focused on a neutral N-alkylated imino sugar, which is not thought to be protonated in its anti-ceramide glucosyltransferase activity, the ring nitrogen atom of DNJ-based inhibitors is probably protonated in the active sites of glycosidases.5

76
A study of the conformation of protonated DNJ in the gas phase may provide further insight into the mechanism of such inhibition.