Unveiling host

Scientific Reports volume 12, Article number: 8169 (2022) Cite this article

949 Accesses

8 Citations

Metrics details

We propose a novel scheme of examining the host–guest–solvent interactions in solution from their gas phase structures. By adopting the permethylated β-cyclodextrin (perm β-CD)–protonated L-Lysine non-covalent complex as a prototypical system, we present the infrared multiple photon dissociation (IRMPD) spectrum of the gas phase complex produced by electrospray ionization technique. In order to elucidate the structure of perm β-CD)/LysH+ complex in the gas phase, we carry out quantum chemical calculations to assign the two strong peaks at 3,340 and 3,560 cm−1 in the IRMPD spectrum, finding that the carboxyl forms loose hydrogen bonding with the perm β-CD, whereas the ammonium group of L-Lysine is away from the perm β-CD unit. By simulating the structures of perm β-CD/H+/L-Lysine complex in solution using the supramolecule/continuum model, we find that the extremely unstable gas phase structure corresponds to the most stable conformer in solution.

Host–guest interactions are of fundamental importance as the principle of molecular recognition and self-assembly1,2,3. “Bottom-up” construction of supramolecular functional materials4,5,6,7 based on host–guest interactions has found wide applications for catalysts (artificial enzymes)8,9, sensing10, drug delivery11,12, gene delivery13,14, bioimaging15,16, stimuli-responsive materials17 and photodynamic therapy18,19. Examining the host–guest configurations has been the focus of many intensive studies20,21,22, however, scrutinizing them in solution phase at molecular level is usually a difficult task to achieve. The two components of host–guest complexes are held together with highly specific binding by non-covalent forces, most commonly by hydrogen bonding. The host components are usually large molecules such as cyclodextrins23,24, calix[n]arenes25,26, or cucurbit[n]urils27. The hydrophobic cavity of a macrocyclic host molecule usually encapsulates some part of the guest, whereas the outer rims of the host tend to be hydrophilic. Details for the interactions and bindings between the functional groups in the host and guest are of specific interest in relation to molecular recognition. One of the most intriguing questions concerning the host–guest complexes would be whether their structures observed in gas-phase may reflect the host–guest–solvent interactions in solution or not. If there exists a close relationship between the structures of host–guest complexes in the two phases, then determining their configurations in the gas phase, which might be more unambiguous than the solution phase structures (at least due to the solvent-free environment) may give invaluable information for the detailed interactions between the host and guest (and also solvent molecules). By addressing this highly critical question, we describe here a novel protocol for examining the host–guest–solvent interactions in solution by elucidating the structures of gas phase host–guest complexes. We employ the perm β-CD/LysH+ complex as a prototypical host–guest system to advance the scheme.

First, we obtain the IRMPD spectra of the gas phase perm β-CD/H+/L-Lys complex produced by the ESI/MS procedure. We determine the gas phase structure of perm β-CD/H+/L-Lys by comparing the calculated IR spectra with the experimentally observed IRMPD28,29,30,31,32,33,34,35,36,37,38,39 spectrum of the perm-CD/ L-LysH+ complex. We carry out quantum chemical calculations for the two intense bands at 3340 and 3560 cm−1, assigning them to stretch modes of the side chain ammonium and the carboxyl group, respectively. Second, we show that Gibbs free energy of the experimentally observed gas phase complex (Complex I) obtained by ESI/mass spectrometry (ESI/MS) from solution phase is much higher (by ~ 30 kcal/mol) than that of the lower Gibbs free energy structure (Complex II), demonstrating that the gas phase complex is highly unstable in thermodynamic sense (as was also observed for gas phase perm β-CD/H+/G, G = alanine40, isoleucine41, H2O42). Third, by adopting the supramolecule/continuum approach, we simulate the structures of perm β-CD/H+/L-Lys complex in solution, thereby finding that the relative thermodynamic stability of perm β-CD/H+/L-Lys complex in solution is the complete reversal of that in the gas phase: The extremely higher Gibbs free energy gas phase structure (Complex I) corresponds to the most stable conformer (Complex A) in solution, and the most stable gas phase conformer (Complex II) becomes a less stable one (Complex B) in solution. This clearly indicates that the origin of the high Gibbs free energy gas phase structure (Complex I) is the configuration (Complex A) of perm β-CD/H+/L-Lys in solution, suggesting that the structure of gas phase host–guest complexes produced by ESI/MS/IRMPD technique43 may be used as a very useful and accurate guide for probing host–guest–solvent interaction in solution.

Figure 1a presents the experimental IRMPD spectrum of gas phase perm β-CD–protonated L-Lys complex in the region > 3100 cm−1. Two intense bands at 3340 and 3560 cm−1 are noticeable. The band at 3560 cm−1 is very similar to carboxyl stretch band (~ 3570 cm−1) of the “free” gas phase amino acids44, indicating that it may be assigned a –OH stretch mode of nearly isolated carbonyl group. The theoretical spectrum (calculated at wB97X-D/6-31G* level of theory, using a scale factor of 0.940868 to correct for the DFT errors and anharmonic effects) of Complex I, which fully accounts for the IRMPD spectrum, is also depicted in Fig. 1a along with the normal modes of the absorption bands in Fig. 1b, and the corresponding gas phase perm β-CD/H+/L-Lys structure. The calculated infrared (IR) spectrum for Complex I exhibits two strong bands at 3348 and 3552 cm−1, which is the characteristic feature of the experimental spectrum in Fig. 1a. The two intense bands at 3226 and 3319 cm−1 in the IR spectrum are assigned as the asymmetric stretch of the ammonium and the carboxyl –OH stretch mode, respectively (Fig. 2a). Very weak absorptions at 3393 and at 3323 cm−1 are asymmetric and symmetric –NH2 stretch modes, respectively, but they seem to be so weak that they are almost indistinguishable from background noise in the experimental spectrum.

(a) Theoretical IR spectrum, (b) normal modes, and (c) two views of the structure of Complex I of gas phase perm β-CD–protonated L-Lys. Red sticks represent the positions of experimental IRMPD bands. Hydrogen atoms in permetylated β-cyclodextrin are omitted for clarity. Red: oxygen, blue: nitrogen, gray: carbon, white: hydrogen.

(a) Theoretical IR spectra and (b) structures of lower Gibbs free energy perm β-CD–protonated L-Lys complexes II and III in gas phase.

The presence or absence of hydrogen bonding between the perm β-CD unit and the functional groups (carbonyl, amino, and ammonium), and their locations with respect to the host perm β-CD are the critical features in the gas phase structures of the complexes. In Complex I (see Fig. 1c) the amino and the side chain ammonium groups are far off the perm β-CD unit, whereas the carboxyl loosely interact with the O atoms of the perm β-CD. The weak interactions (ROH = 2.313 Å) between the carboxyl and the perm β-CD unit was also seen in the structure of gas phase perm β-CD/H+/Isoleucine41 system (with RO..HO = 2.053 Å). The calculated length of hydrogen bonding (RO…HO = 2.313 Å), which is much larger than that observed for ordinary hydrogen bonding (~ 1.70 Å), indicates that interactions between the carboxyl and the perm β-CD unit are quite weak in Complex I, in agreement with the high frequency –OH stretch band at 3552 cm−1 in IRMPD spectrum. The origin of this intriguing loose hydrogen bonding of the carboxyl to the perm β-CD can be accounted for by envisaging an intervening water molecule in solution phase, as discussed below.

Of the thermodynamically more stable (lower Gibbs free energy) conformers of gas phase perm β-CD/H+/L-Lys complexes, Fig. 2 presents the two structures (Complex II, III) deserving close inspection. In Complex II and Complex III the ammonium is encapsulated inside the perm β-CD unit, with the carboxyl and the amino group interacting with the host at its rim, respectively. Figure 2 also shows that the calculated IR frequency of carbonyl –OH stretch in Complex II is enormously red-shifted to 3226 cm−1 (the corresponding normal mode is depicted in Supplementary Information (Figure S1)), as the result of strong hydrogen bonding (RO…HO = 1.723 Å, considerably shorter than RO…HO = 2.313 Å in Complex I) between –CO2H and the perm β-CD unit. On the other hand, the corresponding frequency in Complex III is calculated to be blue-shifted (Fig. 2b) from the experimental band to 3607 cm−1. This indicates that in Complex III the carboxyl is completely isolated, free from the influence of the perm β-CD ring, as can be seen from its structure in Fig. 2 (–CO2H is bent outward away from the perm β-CD unit). In this latter conformer, the α–amino group forms a weak hydrogen bond with the CD O atoms (RO…HN = 2.228 Å), giving the high frequency symmetric and asymmetric –NH2 stretch modes at 3335 and 3424 cm−1, respectively. It should be noted that the calculated IR spectra for Complex II and Complex III are in complete disagreement with the experimental IRMPD spectrum in Fig. 1a.

Table 1 lists the relative Gibbs free energies of the gas phase complexes I–III. We find that Complex II and III is much more stable (lower Gibbs free energy) than other complexes by 27–30 kcal/mol. Two structural differences are to be noted between Complex I and Complex II/III): First, in Complex I interactions between the amino and carboxyl with the perm β-CD unit is much weaker than in Complex II/III. Second, in the latter two complexes the side chain –NH3+ is encapsulated inside the cavity of the host perm β-CD, whereas it is essentially isolated in Complex I. Thus, it seems that these stronger interactions between the LysH+ functional groups and the perm β-CD provide Complex II and III far greater thermodynamic stability than Complex I.

From these findings, it is clear that thermal equilibrium has not been achieved for the gas phase perm β-CD/L-LysH+ complexes. The reason for the striking observation of thermodynamically highly unstable conformer (Complex I), not the thermally much more stable conformers (Complex II, III), is extremely intriguing. One plausible proposition is that when a mixture of perm β-CD and L-LysH+ in solution was being transferred into the gas phase, the complex may not fully relax to the global minimum, kinetically trapped in the local minimum (Complex I). This situation may correspond to the case in which the thermodynamically much less stable (higher Gibbs free energy) gas phase conformer (Complex I) may actually be the most stable configuration of perm β-CD/L-LysH+ complex in solution phase under the influence of solvating molecules, and thus the solution phase perm β-CD/L-LysH+ complex that is structurally close to Complex I is observed in the gas phase before thermal equilibrium is attained.

To unravel the origin of the extremely high Gibbs free energy perm β-CD/L-LysH+ non-covalent complex in gas phase, we adopt a model of solution phase host–guest system by employing the supramolecule/continuum approach, by which we treat the solvent in the first shell around the complex as explicit molecules, and other numerous solvent molecules in the second shell and beyond as continuum. By allowing the gas phase complexes (Complex I, II and III) to interact with increasing number of explicit H2O molecules, we find that the model including 16 explicit H2O molecules fully surrounding the β-CD/L-LysH+ complex in the first shell plus water continuum is the best one with regard to the amounts of computational efforts. Figure 3 and Table 1 present the lowest Gibbs free energy solution phase structures (Complex A, B and C) corresponding to the three gas phase conformers (Complex I, II and III, respectively) obtained by extensive search over the potential landscape of the perm β-CD/L-LysH+ system in solution. (Some other low-lying (local minimum Gibbs free energy) structures of are given in Figure S2) We find that Complex A, conforming to the highly unstable gas phase Complex I, is now calculated to be more stable (Gibbs free energy lower by 2.73 at wB97X-D/6-31G*) than Complex B corresponding to the most stable gas phase Complex II, exhibiting a complete reversal of their relative thermodynamic stability in the gas phase by overcoming the huge deficit of ~ 30 kcal/mol. The Gibbs free energy of the solution phase structure C stemming from the gas phase structure III is calculated to be much higher (> 9.9 kcal/mol). This seems to result from the differences in the structural features of Complex A vs. Complex B/C: In all three complexes the carboxyl interacts both with a water molecule and the perm β-CD unit, but the strength of interactions decreases from A to B to C, the H2O–perm β-CD distance increasing from 1.601 to 2.075 and to 2.432 Å. Moreover, two water molecules directly interact with the side chain –NH3+, and are in close contact with other H2O molecules in a chain-like form in Complex A. In contrast, only a more or less isolated H2O molecule binds to the latter functional group inside the cavity of the perm β-CD unit in Complex B/C. Therefore, we conclude that the most stable Complex A in solution is effused by the ESI method into the gas phase, and that the resulting Complex I is observed in gas phase without thermal relaxation to the lower free energy gas phase structures. From the structure of Complex A in Fig. 3, we may now visualize the host–guest–solvent interactions between the perm β-CD, L-LysH+ and solvent molecules in solution: The L-LysH+ guest situates at the rim of the perm-CD host, not encapsulated inside it. The ammonium forms weak hydrogen bonding with the perm β-CD unit, being surrounded by two water molecules. The amino group is enveloped by 3 H2O molecules off the host perm β-CD. The carboxyl weakly interacts with the perm β-CD bridged by a water molecule, also binding to 3 H2O molecules away from the host. These solvating water molecules are linked by a number of other interconnecting water molecules, together forming the first shell around the perm β-CD/L-LysH+ complex. (Details of interactions between the perm β-CD host, the L-LysH+ guest and solvent molecules in solution are described in Figure S3–S5). It would also be worth noticing that the structural features (locations of the guest functional groups relative to the host) of the host–guest complex in solution phase (Complex A) are more or less retained in gas phase (Complex I), although some bond distances seem to relax closer to the host when the solvent molecules are removed from Complex A to form the gas phase Complex I by the ESI procedure: For example, the distance between an ammonium H atom and a perm β-CD O atom in Complex A becomes shorter (from 2.537 to 1.813 Å), and that between the carboxyl H atom and a perm β-CD O atom relaxes from 3.041 to 2.313 Å (Fig. 4).

Two views of the structures (Complex A, B and C) of perm β-CD–protonated L-Lys complexes in solution, based on the gas phase (a) Complex I, (b) Complex II, and (c) Complex III, respectively. Pale blue background represents solvent continuum.

Comparison of perm-CD–LysH+ configurations in solution (Complex A, wB97X-D/6-31G*) and in gas phase (Complex I, wB97X-D/6-311G*).

In summary, our present study gives an unambiguous clue to the origin of the high Gibbs free energy perm β-CD/H+/Ala, Ile, L-Lys and H2O non-covalent complexes experimentally observed in gas phase: Their source is the perm β-CD/guest complexes in solution from which the gas phase complexes are produced by the ESI/MS before reaching thermal relaxation to equilibrium. The IRMPD and quantum chemical methods were employed here to determine the structures of the gas phase host–guest complexes, but any other techniques may also be useful. These findings offer a profound implication and a novel scheme with a view to solution phase configurations of host–guest complexes: The structures of the gas phase host–guest complexes produced and identified by the ESI/IRMPD spectroscopy carry important and accurate information for their configurations in solution phase. That is, by “freezing” the host–guest complexes and determining their structures in the gas phase, it is possible to ‘fingerprint’ the host–guest–solvent configurations in solution.

Infrared multiple photon dissociation spectroscopy experiments28,29,30,31,32,33,34,35,36,37,38,39 were performed on a 4.7 T ESI-Fourier transform ion cyclotron resonance (FT-ICR) mass spectrometer (Ionspec, Lake Forest, CA, USA). The detailed IRMPD spectroscopy setup was described in detail elsewhere42. Non-covalent complexes of permethylated β-CD , protonated L-Lys and water, i.e., (β-CD + L-Lys + H)+, were generated by electrospraying a solution of 20 mM permethylated β-CD and L-Lys in 49:49:2 (v/v/v) H2O:MeOH:AcOH. A permethylated β-CD was chosen as a main target of study instead of native β-CD in order to simplify the IR spectrum, and thus to facilitate the interpretation of the IR spectrum. Non-covalent complexes were isolated in the ICR cell using stored waveform inverse Fourier transform (SWIFT) waveforms and stored in the ICR cell for 25 s. The stored non-covalent complexes were irradiated using IR laser light pulses (10 Hz). A photodissociation spectrum was obtained by recording the abundance of the (β-CD + L-Lys + H+) peak with (A) and without (A0) infrared laser irradiation as a function of the irradiated laser wavenumber in cm−1; ion-depletion measurement. The abundance of the precursor non-covalent complex ion was monitored in 5 nm steps. The abundance (A and A0) of the precursor complex ion was read directly from the data acquisition software without any data processing such as apodization or zero-filling.

An IRMPD spectrum was obtained in the wavenumber region from 3050 to 3750 cm−1, using a commercial OPO laser (single stage optical parametric oscillator with an unseeded cavity design and a bulk potassium titanyl arsenate (KTA) crystal, IROPO, OPOTek, Carlsbad, CA, USA). For the IRMPD spectrum, the recorded depletion of precursor complex ions was normalized using the measured (and smoothed) laser power assuming a linear dependence. The IR laser irradiation time was adjusted using a mechanical shutter (UNIBLITZ, Vincent Associates, Rochester, NY, USA) that was synchronized with an FT-ICR TTL signal with some delay time. The IR laser light was collimated using a series of irises with adjustable apertures and was carefully aligned to maximize the overlap between ions in the ICR cell and the laser light. The IR laser was focused using a CaF2 lens with a 1 m focal length.

Permethylated β-CD was purchased from Sejin Chemical Industry Co. (Seoul, Korea) and was used without purification. Water and methanol were of HPLC grade (Fischer, Seoul, Korea), acetic acid and L-Lys were purchased from Sigma (Seoul, Korea).

All calculations were conducted using the wB97X-D functional, which may treat weak interactions (hydrogen bonding included) very well, with 6-311G** and 6-31G** basis set implemented in Gaussian09 suite of programs45. For the structures in solution we also carried out single point calculations using the wB97X-D/6-311G** method at the wB97X-D/6-31G** optimized geometry. We obtained the conformers of gas phase perm β-CD/H+/L-Lys based on the structures of the perm β-CD unit described in our previous work42, wherein the NVT ensemble molecular dynamics simulations and PM6 calculations were performed for the gas phase complexes, using the energy window of 35.4 kcal/mol for the lowest 30 conformers. A scale factor of 0.940868 was used to best fit the calculated IR frequencies to experimental spectrum by the least square procedure. All structures were obtained by verifying all real vibrational frequencies. For the corresponding structures of perm β-CD/LysH+ complex in solution, we used the supramolecule/continuum approach, in which the solvent molecules directly interacting with –CO2H, -NH2 and side chain ammonium groups in the first shell around the complex are treated as explicit molecules, and other numerous H2O molecules in the second shell and beyond as water continuum by the SMD method46. Although we used as solvent the mixture of water, methanol, and acetic acid (49:49:2 volume ratio), we modelled our present system as the perm β-CD/LysH+ complex in water for simplicity and by invoking two presumptions: First, methanol would preferentially evaporate as the ESI droplets becomes smaller due to higher volatility. Second, H2O molecules would mostly constitute the first shell around the perm β-CD/LysH+ complex because the two -OH groups per H2O molecule would form much stronger hydrogen bonding with the functional groups (ammonium, amino and carboxyl) than CH3OH would. We proceeded by adding water molecules to the functional groups in perm β-CD/H+/L-Lys complex, beginning with 4 H2O molecules, searching for the lowest Gibbs free energy configurations at each step. The Gibbs free energies were calculated by the standard options (bulk electrostatic contribution from a self-consistent reaction field treatment by superposition of nuclear-centered spheres, and short-range interactions between the solute and solvent molecules in the first shell obtained by geometry-dependent atomic surface tensions, etc.) for the keyword SCRF = (SMD, solvent = water). Comparisons are made with those optimized from initial configurations with by randomly added water molecules using the VMD (visual molecular dynamics) program.

The datasets generated and/or analysed during the current study are available in the Computational Chemistry Results Repository (https://www.iochem-bd.org/), https://doi.org/10.19061/iochem-bd-6-124.

Robert, F. How far can we push chemical self-assembly?. Science 309, 95 (2005).

Article Google Scholar

Percec, V. et al. Self-assembly of amphiphilic dendritic dipeptides into helical pores. Nature 430, 764–768 (2004).

Article ADS CAS PubMed Google Scholar

Ma, X. & Zhao, Y. Biomedical applications of supramolecular systems based on host–guest interactions. Chem. Rev. 115, 7794–7839 (2015).

Article CAS PubMed Google Scholar

Yang, H., Yuan, B., Zhang, X. & Scherman, O. A. Supramolecular chemistry at interfaces: Host–guest interactions for fabricating multifunctional biointerfaces. Acc. Chem. Res. 47, 2106–2115 (2014).

Article CAS PubMed Google Scholar

Zhan, W., Wei, T., Yu, Q. & Chen, H. Fabrication of supramolecular bioactive surfaces via β-cyclodextrin-based host–guest interactions. ACS Appl. Mater. Interfaces 10, 36585–36601 (2018).

Article CAS PubMed Google Scholar

Kim, C., Agasti, S. S., Zhu, Z., Isaacs, L. & Rotello, V. M. Recognition-mediated activation of therapeutic gold nanoparticles inside living cells. Nat. Chem. 2, 962–966 (2010).

Article CAS PubMed PubMed Central Google Scholar

Schalley, C. A. Supramolecular chemistry goes gas phase: the mass spectrometric examination of noncovalent interactions in host–guest chemistry and molecular recognition. Int. J. Mass Spectrom. 194, 11–39 (2000).

Article CAS Google Scholar

Ortega-Caballero, F., Rousseau, C., Christensen, B., Petersen, T. E. & Bols, M. Remarkable supramolecular catalysis of glycoside hydrolysis by a cyclodextrin cyanohydrin. J. Am. Chem. Soc. 127, 3238–3239 (2005).

Article CAS PubMed Google Scholar

Ortega-Caballero, F., Bjerre, J., Laustsen, L. S. & Bols, M. Four orders of magnitude rate increase in artificial enzyme-catalyzed aryl glycoside hydrolysis. J. Org. Chem. 70, 7217–7226 (2005).

Article CAS PubMed Google Scholar

Shahgaldian, P. & Pieles, U. Cyclodextrin derivatives as chiral supramolecular receptors for enantioselective sensing. Sensors 6, 593–615 (2006).

Article ADS CAS PubMed Central Google Scholar

Hu, Q.-D., Tang, G.-P. & Chu, P. K. Cyclodextrin-based host–guest supramolecular nanoparticles for delivery: from design to applications. Acc. Chem. Res. 47, 2017–2025 (2014).

Article CAS PubMed Google Scholar

Uekama, K., Hirayama, F. & Irie, T. Cyclodextrin drug carrier systems. Chem. Rev. 98, 2045–2076 (1998).

Article CAS PubMed Google Scholar

Miyata, K., Nishiyama, N. & Kataoka, K. Rational design of smart supramolecular assemblies for gene delivery: chemical challenges in the creation of artificial viruses. Chem. Soc. Rev. 41, 2562–2574 (2012).

Article CAS PubMed Google Scholar

Davis, M. E. et al. Evidence of RNAi in humans from systemically administered siRNA via targeted nanoparticles. Nature 464, 1067–1070 (2010).

Article ADS CAS PubMed PubMed Central Google Scholar

Carvalho, C. P., Uzunova, V. D., Da Silva, J. P., Nau, W. M. & Pischel, U. A photoinduced pH jump applied to drug release from cucurbit [7] uril. Chem. Commun. 47, 8793–8795 (2011).

Article CAS Google Scholar

Jiao, D. et al. Supramolecular peptide amphiphile vesicles through host–guest complexation. Angew. Chemie Int. Ed. 51, 9633–9637 (2012).

Article CAS Google Scholar

Zhang, J. et al. Dual stimuli-responsive supramolecular self-assemblies based on the host–guest interaction between β-cyclodextrin and azobenzene for cellular drug release. Mol. Pharm. 17, 1100–1113 (2020).

Article CAS PubMed Google Scholar

Qu, D.-H., Wang, Q.-C., Zhang, Q.-W., Ma, X. & Tian, H. Photoresponsive host–guest functional systems. Chem. Rev. 115, 7543–7588 (2015).

Article CAS PubMed Google Scholar

Kitagishi, H. et al. Cellular uptake of octaarginine-conjugated tetraarylporphyrin included by per-O-methylated β-cyclodextrin. Org. Biomol. Chem. 11, 3203–3211 (2013).

Article CAS PubMed Google Scholar

Cescutti, P., Garozzo, D. & Rizzo, R. Effect of methylation of β-cyclodextrin on the formation of inclusion complexes with aromatic compounds. An ionspray mass spectrometry investigation. Carbohydr. Res. 302, 1–6 (1997).

Article CAS Google Scholar

Sun, W., Cui, M., Liu, S., Song, F. & Elkin, Y. N. Electrospray ionization mass spectrometry of cyclodextrin complexes with amino acids in incubated solutions and in eluates of gel permeation chromatography. Rapid Commun. mass Spectrom. 12, 2016–2022 (1998).

3.0.CO;2-1" data-track-action="article reference" href="https://doi.org/10.1002%2F%28SICI%291097-0231%2819981230%2912%3A24%3C2016%3A%3AAID-RCM430%3E3.0.CO%3B2-1" aria-label="Article reference 21" data-doi="10.1002/(SICI)1097-0231(19981230)12:243.0.CO;2-1">Article ADS CAS Google Scholar

Mele, A. & Malpezzi, L. Noncovalent association phenomena of 2, 5-dihydroxybenzoic acid with cyclic and linear oligosaccharides. A matrix-assisted laser desorption/ionization time-of-flight mass spectrometric and X-ray crystallographic study. J. Am. Soc. Mass Spectrom. 11, 228–236 (2000).

Article CAS PubMed Google Scholar

Armstrong, D. W., Ward, T. J., Armstrong, R. D. & Beesley, T. E. Separation of drug stereoisomers by the formation of beta-cyclodextrin inclusion complexes. Science (80-.) 232, 1132–1135 (1986).

Article ADS CAS Google Scholar

Lee, J., Lee, S.-S., Lee, S. & Oh, H. B. Noncovalent complexes of cyclodextrin with small organic molecules: Applications and insights into host–guest interactions in the gas phase and condensed phase. Molecules 25, 4048 (2020).

Article CAS PubMed Central Google Scholar

Nimse, S. B. & Kim, T. Biological applications of functionalized calixarenes. Chem. Soc. Rev. 42, 366–386 (2013).

Article CAS PubMed Google Scholar

Guo, D.-S. & Liu, Y. Supramolecular chemistry of p-sulfonatocalix [n] arenes and its biological applications. Acc. Chem. Res. 47, 1925–1934 (2014).

Article CAS PubMed Google Scholar

El-Sheshtawy, H. S. et al. A supramolecular approach for enhanced antibacterial activity and extended shelf-life of fluoroquinolone drugs with cucurbit [7] uril. Sci. Rep. 8, 1–10 (2018).

Article CAS Google Scholar

Oomens, J., Sartakov, B. G., Meijer, G. & Von Helden, G. Gas-phase infrared multiple photon dissociation spectroscopy of mass-selected molecular ions. Int. J. Mass Spectrom. 254, 1–19 (2006).

Article CAS Google Scholar

Bush, M. F., O’Brien, J. T., Prell, J. S., Saykally, R. J. & Williams, E. R. Infrared spectroscopy of cationized arginine in the gas phase: direct evidence for the transition from nonzwitterionic to zwitterionic structure. J. Am. Chem. Soc. 129, 1612–1622 (2007).

Article CAS PubMed PubMed Central Google Scholar

Polfer, N. C. Infrared multiple photon dissociation spectroscopy of trapped ions. Chem. Soc. Rev. 40, 2211–2221 (2011).

Article CAS PubMed Google Scholar

Azargun, M. & Fridgen, T. D. Guanine tetrads: an IRMPD spectroscopy, energy resolved SORI-CID, and computational study of M (9-ethylguanine) 4+(M = Li, Na, K, Rb, Cs) in the gas phase. Phys. Chem. Chem. Phys. 17, 25778–25785 (2015).

Article CAS PubMed Google Scholar

Oh, H. et al. Secondary and tertiary structures of gaseous protein ions characterized by electron capture dissociation mass spectrometry and photofragment spectroscopy. Proc. Natl. Acad. Sci. 99, 15863–15868 (2002).

Article ADS CAS PubMed PubMed Central Google Scholar

Yao, H., Steill, J. D., Oomens, J. & Jockusch, R. A. Infrared multiple photon dissociation action spectroscopy and computational studies of mass-selected gas-phase fluorescein and 2′,7′-dichlorofluorescein ions. J. Phys. Chem. A 115, 9739–9747 (2011).

Article CAS PubMed Google Scholar

Cheng, R., Loire, E. & Fridgen, T. D. Hydrogen bonding in alkali metal cation-bound i-motif-like dimers of 1-methyl cytosine: An IRMPD spectroscopic and computational study. Phys. Chem. Chem. Phys. 21, 11103–11110 (2019).

Article CAS PubMed Google Scholar

Seo, J. et al. Retention of native protein structures in the absence of solvent: a coupled ion mobility and spectroscopic study. Angew. Chemie Int. Ed. 55, 14173–14176 (2016).

Article CAS Google Scholar

Bush, M. F. et al. Infrared spectroscopy of cationized lysine and ε-N-methyllysine in the gas phase: effects of alkali-metal ion size and proton affinity on zwitterion stability. J. Phys. Chem. A 111, 7753–7760 (2007).

Article CAS PubMed Google Scholar

Scuderi, D., Lepere, V., Piani, G., Bouchet, A. & Zehnacker-Rentien, A. Structural characterization of the UV-induced fragmentation products in an ion trap by infrared multiple photon dissociation spectroscopy. J. Phys. Chem. Lett. 5, 56–61 (2014).

Article CAS PubMed Google Scholar

Sen, A., Le Barbu-Debus, K., Scuderi, D. & Zehnacker-Rentien, A. Mass spectrometry study and infrared spectroscopy of the complex between camphor and the two enantiomers of protonated alanine: the role of higher-energy conformers in the enantioselectivity of the dissociation rate constants. Chirality 25, 436–443 (2013).

Article CAS PubMed Google Scholar

Burt, M. et al. Assessing the impact of anion–π effects on phenylalanine ion structures using IRMPD spectroscopy. Phys. Chem. Chem. Phys. 16, 24223–24234 (2014).

Article CAS PubMed Google Scholar

Lee, S.-S. et al. Chiral differentiation of d-and l-alanine by permethylated β-cyclodextrin: IRMPD spectroscopy and DFT methods. Phys. Chem. Chem. Phys. 19, 14729–14737 (2017).

Article CAS PubMed Google Scholar

Lee, S.-S. et al. Chiral differentiation of D-and L-isoleucine using permethylated β-cyclodextrin: Infrared multiple photon dissociation spectroscopy, ion-mobility mass spectrometry, and DFT calculations. Phys. Chem. Chem. Phys. 20, 30428–30436 (2018).

Article CAS PubMed Google Scholar

Lee, S.-S. et al. Infrared multiple photon dissociation spectroscopy and density functional theory (DFT) studies of protonated permethylated β-cyclodextrin–water non-covalent complexes. Phys. Chem. Chem. Phys. 16, 8376–8383 (2014).

Article CAS PubMed Google Scholar

Bokatzian, S. S., Stover, M. L., Plummer, C. E., Dixon, D. A. & Cassady, C. J. An experimental and computational investigation into the gas-phase acidities of tyrosine and phenylalanine: Three structures for deprotonated tyrosine. J. Phys. Chem. B 118, 12630–12643 (2014).

Article CAS PubMed Google Scholar

Linder, R., Nispel, M., Häber, T. & Kleinermanns, K. Gas-phase FT-IR-spectra of natural amino acids. Chem. Phys. Lett. 409, 260–264 (2005).

Article ADS CAS Google Scholar

Frisch, M., et al. Gaussian 09, revision D. 01 (2009).

Marenich, A. V., Cramer, C. J. & Truhlar, D. G. Universal solvation model based on solute electron density and on a continuum model of the solvent defined by the bulk dielectric constant and atomic surface tensions. J. Phys. Chem. B 113, 6378–6396 (2009).

Article CAS PubMed Google Scholar

Download references

This work was supported by grants from the National Research Foundation of Korea funded by the Ministry of Education (NRF-2019R1F1A1057609, 2018R1A6A1A03024940, 2021R1A2C2007397), and the KISTI Supercomputing Center.

These authors contributed equally: Hyoju Choi and Young-Ho Oh.

Department of Applied Chemistry, Kyung Hee University, Gyeonggi, 17104, Republic of Korea

Hyoju Choi, Young-Ho Oh, Sung-Sik Lee & Sungyul Lee

Department of Chemistry, Sogang University, Seoul, 121-742, Republic of Korea

Soojin Park & Han Bin Oh

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

H.C. and Y.-H.O. carried out calculations for the structures of the complexes. S.P. produced the IRMPD spectra of the complex in gas phase. S.-S.L. helped to obtain the structures of the gas phase complexes. H.B.O. designed and performed the experimental procedures (ESI/MS/IRMPD). S.L. devised the overall concept of the study and the modeling of the system.

Correspondence to Han Bin Oh or Sungyul Lee.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Choi, H., Oh, YH., Park, S. et al. Unveiling host–guest–solvent interactions in solution by identifying highly unstable host–guest configurations in thermal non-equilibrium gas phase. Sci Rep 12, 8169 (2022). https://doi.org/10.1038/s41598-022-12226-0

Download citation

Received: 07 February 2022

Accepted: 05 May 2022

Published: 17 May 2022

DOI: https://doi.org/10.1038/s41598-022-12226-0

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.