ing the Hydrophobic Effect with Hydrogen Bonding - ACS Publications

ing the Hydrophobic Effect with Hydrogen Bonding. Huan Yao,†,‡,¶ Hua Ke,†,¶ Xiaobin Zhang,ǁ,¶ San-Jiang Pan,† Ming-Shuang Li,† Liu-Pan Y...
0 downloads 0 Views 3MB Size
Article Cite This: J. Am. Chem. Soc. 2018, 140, 13466−13477

pubs.acs.org/JACS

Molecular Recognition of Hydrophilic Molecules in Water by Combining the Hydrophobic Effect with Hydrogen Bonding Huan Yao,†,‡,¶ Hua Ke,†,¶ Xiaobin Zhang,∥,¶ San-Jiang Pan,† Ming-Shuang Li,† Liu-Pan Yang,† Georg Schreckenbach,∥ and Wei Jiang*,† †

Downloaded via UNIV OF GOTHENBURG on October 17, 2018 at 06:26:56 (UTC). See https://pubs.acs.org/sharingguidelines for options on how to legitimately share published articles.

Department of Chemistry and Shenzhen Grubbs Institute, Southern University of Science and Technology, Shenzhen 518055, China ‡ School of Chemistry and Chemical Engineering, Harbin Institute of Technology, Harbin, 150001, China ∥ Department of Chemistry, University of Manitoba, Winnipeg, Manitoba R3T 2N2, Canada S Supporting Information *

ABSTRACT: During the last half a century, great achievements have been made in molecular recognition in parallel with the invention of numerous synthetic receptors. However, the selective recognition of hydrophilic molecules in water remains a generally accepted challenge in supramolecular chemistry but is commonplace in nature. In an earlier Communication [Huang et al. J. Am. Chem. Soc. 2016, 138, 14550], we reported a pair of endo-functionalized molecular tubes that surprisingly prefer highly hydrophilic molecules over hydrophobic molecules of a similar size and shape. The hydrophobic effect and hydrogen bonding were proposed to be responsible, but their exact roles were not fully elucidated. In this Article, we present a thorough study on the binding behavior of these molecular tubes toward 44 hydrophilic molecules in water. Principal component analysis reveals that the binding strength is weakly correlated to the hydrophobicity, volume, surface area, and dipole moment of guests. Furthermore, molecular dynamics simulations show the hydrophobic effect through releasing the poorly hydrogen-bonded cavity water contributes to the binding of all the hydrophilic molecules, while hydrogen bonding differentiates these molecules and is thus the key to achieve a high selectivity toward certain hydrophilic molecules over other molecules with a similar size and shape. Therefore, a good guest for these molecular tubes should meet the following criteria: the hydrogen-bonding sites should be complementary, and the molecular volume should be large enough to expel all the cavity water but not too large to cause steric hindrance. This rule of thumb may also be used to design a selective receptor for certain hydrophilic molecules. Following these guidelines, a “best-fit” guest was found for the syn-configured molecular tube with a binding constant as high as 106 M−1.

1. INTRODUCTION

interactions with water molecules. It would cost high desolvation energy to recognize them. Since 1967, when Pedersen reported crown ethers,13 macrocyclic receptors have been the major tools in supramolecular chemistry. During the last five decades, numerous synthetic macrocyclic receptors have been reported and intensively investigated,14 including crown ethers,15 cyclodextrins,4 calix[n]arenes,16 cyclotriveratrylenes,17 calixpyrroles,18 cucurbit[n]urils,5 “blue box”,19 heterocalixaromatics,20 pillar[n]arenes,21 “Texas-sized box”,22 bambus[n]uril,23 cyanostar,24 corona[n]arene,25 helicarene,26 biphen[n]arene,27 oxatub[n]arene,28 and many others. Great achievements have been made, which lays the basis for their applications in selfassembly,29 interlocked structures,30 molecular machines,31 supramolecular catalysis,32 molecular sensors,33 and supra-

1

Molecular recognition is the basis of life and the central topic of supramolecular chemistry. In particular, molecular recognition in aqueous solutions2 is attractive, since it may provide simplified models to understand complex biomolecular recognition and the promise for applications in biomedical and environmental sciences. The major tools for molecular recognition in water are water-soluble macrocyclic receptors,3−5 but water-soluble self-assembled hosts,6 foldamers,7 and tweezers8 are used as well. The recognition of cations, anions, and hydrophobic molecules/groups in water has been intensively studied. In contrast, the selective recognition of hydrophilic molecules in water remains a generally accepted challenge in the field of supramolecular chemistry9 and is thus substantially less explored.10−12 This class of molecules, including persistent environmental contaminants, drug molecules, and biomarkers for diseases, is very important. Hydrophilic molecules, usually containing several heteroatoms, are heavily solvated in water through the formation of intensive © 2018 American Chemical Society

Received: August 24, 2018 Published: September 24, 2018 13466

DOI: 10.1021/jacs.8b09157 J. Am. Chem. Soc. 2018, 140, 13466−13477

Article

Journal of the American Chemical Society molecular materials.34 Nevertheless, water-soluble macrocycles are relatively rare.3−5 Most of the known water-soluble hosts are inefficient in the selective recognition of hydrophilic molecules in water. The known macrocyclic and self-assembled hosts can be roughly classified into two categories: the ones with hydrophobic cavities, such as cyclodextrins, and the ones with polar binding sites, such as crown ethers. Surprisingly, some of the hosts with only a hydrophobic cavity have been reported to be able to recognize hydrophilic molecules in water.11,12 Their binding selectivity heavily depends on the size and shape of guests, and rich information in the functional groups is largely ignored. Most importantly, the binding affinity to hydrophilic molecules is usually significantly weaker than that of the hydrophobic guests of similar size and shape. The hosts with polar binding sites can often bind polar molecules in nonpolar solvents, for example, through hydrogen bonding.35 However, the binding is significantly attenuated or even fully eliminated when transferred into water.36 The problem here is that the polar interactions, such as hydrogen bonding, are challenging to be exploited for molecular recognition in water.9e,37 The reasons are as follows: a) water competes with polar binding sites of synthetic receptors, and b) water is very polar (1.85 D, ε = 78), which significantly weakens polar interactions in host−guest complexes. In contrast, hydrophilic molecules/groups can be well recognized by bioreceptors, even with high selectivity. Their binding pockets feature a deep hydrophobic cavity containing polar binding sites, such as hydrogen-bonding sites. In this way, hydrogen-bonding sites are shielded from bulk water, and the surroundings are relatively nonpolar. Together with the hydrophobic effect, hydrogen bonding can significantly enhance the binding affinity and selectivity. This is expressed in the binding pair of avidin and biotin (Figure 1a).38 The polar groups of biotin are complemented by the hydrogenbonding sites in the cavity of avidin, and the residual hydrophobic groups are also satisfied in the deep hydrophobic cavity. Thus, a high binding affinity (1013−1015 M−1) is achieved, making this binding pair very popular in many

applications.39 Analogously, natural lectins use a similar structural arrangement (many hydrogen-bonding sites shielded inside a hydrophobic cavity) to achieve selective recognition of highly hydrophilic carbohydrates.40 However, water-soluble synthetic hosts having polar binding sites in a hydrophobic cavity are very rare41,42 and difficult to synthesize.43 By mimicking the binding pockets of bioreceptors, we proposed endo-functionalized44 molecular tubes to meet this criterion (Figure 1a): positioning hydrogen-bonding donors inside a deep tubular cavity.45 In an earlier attempt, a pair of water-soluble endo-functionalized molecular tubes (1a and 1b, Figure 1b), which was first reported by Glass and coworkers for the recognition of lipid molecules in water through the hydrophobic effect,46 has been used to achieve the selective recognition of a series of highly hydrophilic molecules in water,47−49 including a Group 2B carcinogen and persistent environmental contaminant  1,4-dioxane. The binding behavior of these endo-functionalized molecular tubes is drastically different from that of other macrocyclic and selfassembled hosts. For example, the binding affinity to 1,4dioxane was significantly increased when transferring the compounds from chloroform (10−102 M−1) to water (103− 104 M−1); a very high selectivity was achieved for the highly hydrophilic molecule 1,4-dioxane over other more hydrophobic molecules with a similar size and shape, such as 1,3dioxane and tetrahydropyran.47 This earlier research, reported as a Communication, leaves several unanswered questions. We proposed that the hydrophobic effect and hydrogen bonding are responsible for the high binding affinity and high selectivity toward the highly hydrophilic molecule 1,4-dioxane. However, what are their exact roles? In the literature, quite weak C−H···π interactions have been suggested to be responsible for the binding of hydrophilic molecules in water.11b,d Is this reasonable? Even with molecular tubes 1a and 1b, is hydrogen bonding the major contributor for the binding of all the hydrophilic molecules? What is the role of hydrogen bonding, especially in the high binding selectivity among molecules of similar size and shape? How can the hydrophobic effect contribute to the binding of hydrophilic molecules? The release of the cavity water has been proposed to be responsible for this. However, evidence for the state of the cavity water is needed. Can we achieve a stronger binding (Ka > 105 M−1) with these molecular tubes? Moreover, the binding pockets of these endofunctionalized molecular tubes closely resemble the protein binding pockets. Further in-depth research will not only provide general guidelines for the design of efficient receptors for biologically and environmentally relevant hydrophilic molecules but also help in understanding molecular recognition in complex biological systems. In this research, we would like to answer the abovementioned questions by studying in detail the binding behaviors of the molecular tubes toward 44 small hydrophilic molecules in water. Principal component analysis (PCA) was employed to reveal the factors that affect the binding strength of hydrophilic molecules. The roles of the hydrophobic effect and hydrogen bonding were analyzed, and general guidelines for selectively recognizing hydrophilic molecules in water were obtained. These were further used to find a “best-fit” guest which can tightly bind these molecular tubes.

Figure 1. a) Bioinspired design of endo-functionalized molecular tubes. b) Chemical structures of the molecular tubes. 13467

DOI: 10.1021/jacs.8b09157 J. Am. Chem. Soc. 2018, 140, 13466−13477

Article

Journal of the American Chemical Society

Figure 2. Chemical structures of the small molecules involved in this research as guests. Guest names or abbreviations used in other places can be found here.

2. RESULTS AND DISCUSSION Determination of Binding Constants and Thermodynamic Parameters. The binding targets of this research are focused on neutral hydrophilic molecules. To meet this criterion, 44 guests containing minimal hydrophobic groups were selected (Figure 2). Their solubility data in water can be found in Table 1. Most of them are highly hydrophilic: 14 guests are miscible with water, and 23 guests have a solubility in water of over 0.1 M. These guests can be divided into four classes: a) linear/acyclic molecules, b) nonaromatic heterocycles, c) aromatic heterocycles, and d) aromatic hydrocarbons. The heteroatoms were focused on N, O, and S to emphasize the role of hydrogen bonding, but 1,2-dichloroethane, benzene, and toluene were also used as guests. All these guests undergo upfield chemical shifts in the 1H NMR spectra of their 1:1 mixtures with 1a/1b (Figures S1− S72), suggesting they are inside the cavities of the molecular tubes and experience the shielding effect of the four naphthalenes. Their binding constants were further determined. To be consistent with the earlier research,47 the binding constants and thermodynamic parameters for most of the guests were determined in deionized H2O or D2O. The pH values of the solutions containing 1a and 1b (0.5 mM) are 7.83 and 7.88, respectively. This suggests that most of the guests exist in the neutral form in these solutions (for the pKas of some basic guests, see Table S1). However, the pKas of guests DEA, morpholine, and piperazine are 10.84, 8.50, and 9.73, respectively. To avoid protonating these guests, their binding constants were determined in a phosphate buffer (pH = 12.0). The binding data are listed in Tables 1 and S2−S4. Most of the guests induce a fluorescence enhancement of the hosts. Therefore, fluorescence titrations (FL, Figures S73−S122) were used to determine their binding constants. For the guests that do not induce a fluorescence response, NMR titrations (Figures S130−S185) were performed to afford the binding constants. The binding constants of guests DCE and tetrazoles were obtained by isothermal titration calorimetry (ITC). Of all these guests, only 18 guests’ thermodynamic parameters (Ka, ΔG°, ΔH°, and ΔS°) can be determined by ITC titrations. For the case that the binding constant is smaller than 103 M−1 and

the guests are highly soluble in water, the ITC titration method for low affinity binding introduced by Turnbull and Daranas was adopted.52 A 1:1 binding stoichiometry was assumed for all the guests, which is supported by great nonlinear fittings and Job’s plots of the molecular tubes with the very small guests: acetonitrile, ethyl ether, and ethyl acetate (Figures S123−129). All the titration experiments were repeated three times. The binding constants from different methods are consistent with each other. The binding data for guests 1,4dioxane, 1,3-dioxane, DMF, DMSO, acetone, THF, THP, oxetane, and oxepane, as determined by fluorescence, NMR, or ITC titrations, have been reported earlier47 but were included here for comparison. Principal Component Analysis (PCA). Most of the binding constants are within the range between 10 and 104 M−1. This is very decent when one considers most of these molecules are highly hydrophilic and not usually chosen as guests for molecular recognition. Only three guests (MDO, 1,4-dioxane, and toluene) show binding constants above or close to 104 M−1. Eleven guests (DEA, dimethoxymethane, EA, DMC, 1,4-dioxene, pyrrole, thiophene, thiazole, tetrazole, pyrazine, and benzene) have binding constants in the range of 103−104 M−1. All the other binding constants are below 102 M−1. This suggests that a high binding selectivity exists. To reveal the guest parameters that affect the binding strength, PCA was performed to elucidate the correlations of the logarithm of the binding constants (log Ka) with six molecular descriptors: the logarithm of the octanol:water partition coefficient (log P), dipole moment (μ), molecular volume (V), surface area (S), polarizability (α), and asphericity (Ω A). PCA provides a graphic representation of the correlations that exist in the data set.53 This method has been widely applied in analyzing chemical sensor array data.53 When the correlation coefficient (r) between two parameters is close to 1 or −1, a strong correlation exists, but the two parameters make a positive or negative contribution to each other. PCA has recently been used by Nitschke et al.12 to elucidate the factors governing guest binding strength in the hydrophobic cavity of a self-assembled host. They revealed that log Ka is highly correlated to μ and log P. Lower μ and higher 13468

DOI: 10.1021/jacs.8b09157 J. Am. Chem. Soc. 2018, 140, 13466−13477

Article

Journal of the American Chemical Society

Table 1. Summary of Molecular Descriptors and Binding Constants (Ka) of the Guests with 1a and 1b in Water at 25 °C guest ethanol isopropyl alcohol acetonee diethy ether DEAf acetonitrile DMSOe DMFe N-ethylacetamide nitromethane DCE dimethoxymethane 2,2-dimethoxypropane EA DMC oxetanee THFe THPe oxepanee cyclohexanone thiane 1,3-dioxolane MDO propylene carbonate 1,3-dioxanee 1,4-dioxene 1,4-dioxanee morpholinef piperazinef s-trioxane furan pyrrole thiophene oxazole thiazole imidazole pyrazole 1,2,4-triazole tetrazole pyridine pyrimidine pyrazine benzene toluene

V (Å3)a 76.54 100.13 91.26 123.58 129.91 68.73 102.16 105.78 130.92 74.53 107.87 113.25 152.34 125.99 114.94 88.20 107.45 127.01 148.63 140.93 140.33 96.85 117.76 119.00 115.90 106.94 116.07 124.32 130.42 104.50 93.86 97.36 107.02 85.54 101.76 91.92 91.76 85.74 83.11 111.97 105.06 104.15 116.30 141.22

ΩA a 0.11373 0.03907 0.05329 0.17298 0.17562 0.20998 0.03143 0.08156 0.13133 0.05675 0.21466 0.08744 0.02584 0.12646 0.13347 0.04207 0.04022 0.04092 0.03939 0.06941 0.04774 0.04306 0.06975 0.06329 0.04233 0.04980 0.04253 0.04235 0.04224 0.04396 0.06254 0.06252 0.06994 0.06260 0.07153 0.06258 0.06252 0.06259 0.06255 0.06257 0.06255 0.06280 0.06250 0.08711

S (Å2)a 93.94 113.95 107.99 139.41 142.89 85.62 116.48 121.33 144.55 91.54 121.44 445.61 157.50 140.88 128.20 102.57 119.52 135.38 151.49 147.72 146.85 109.98 131.51 133.46 126.15 120.68 126.70 130.78 134.82 116.30 105.94 110.46 117.95 100.66 113.00 105.26 105.34 100.03 95.43 122.09 117.14 117.32 126.89 147.37

SMAX (M)b d

misc misc misc 8.15 × 10−1 misc misc misc misc 1.32g 1.82 8.61 × 10−2 3.02 × 10−1 7.22 × 10−2 9.08 × 10−1 9.44 × 10−1g misc misc 9.31 × 10−1 2.68 × 10−2 2.34 × 10−1 1.27 × 10−2 3.74 7.94 × 10−1g 1.71 9.94 × 10−1 miscg misc misc misc 1.94 1.47 × 10−1 7.01 × 10−1 3.59 × 10−2 1.55 6.32 × 10−1 2.34 1.20 6.14 3.08 × 10−1 misc 3.58 2.72 2.29 × 10−2 5.71 × 10−3

log Pb −0.30 0.05 −0.24 0.89 0.58 −0.34 −1.35 −1.01 −0.19g −0.33 1.48 0.00 1.38 0.73 0.07g −0.14 0.46 0.82 1.92 0.81 2.28 −0.37 0.08g −0.41 0.18 −0.39g −0.27 −0.86 −1.50 −0.43 1.34 0.75 1.81 0.12 0.44 −0.08 0.26 −0.58 −0.60 0.65 −0.40 −0.26 2.13 2.73

μ (D)c 1.88 1.89 3.09 1.38 0.48 3.88 4.92 4.44 4.10 3.85 0.00 2.38 2.30 1.94 0.11 2.50 2.26 1.83 1.70 3.42 2.37 1.48 1.35 6.16 1.44 0.95 0.00 2.01 0.00 2.65 1.08 1.90 1.02 1.53 1.32 3.97 2.47 2.99 5.77 2.48 2.61 0.00 0.00 0.36

α (Bohr3)c

Ka (1a, M−1)

25.95 36.36 33.22 48.28 50.89 22.47 44.68 40.52 50.95 24.39 40.64 40.67 60.81 47.95 40.37 32.72 43.16 53.97 64.13 61.31 65.96 36.21 46.61 47.59 46.19 43.75 47.23 50.72 54.09 38.87 36.05 39.68 49.12 32.51 45.91 36.16 36.10 32.49 29.72 49.41 45.33 45.97 53.33 65.49

± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ±

9 9 100 290 1500 8 31 310 65 17 270 6700 49 3400 2300 300 230 68 380 35 65 430 22000 51 210 1800 14000 190 210 47 490 1500 4000 130 1400 330 180 210 1800 590 62 2200 1300 9500

h

2 4h 10 39h 120i 3h 4 31 10h 8h 23j 140i 13h 24i 110i 30 25 7 40 5h 6h 39i 730i 5i 21 73i 1500 8i 17i 3h 23i 59i 230i 22i 220i 33h 9i 21i 650j 48i 5i 38i 190i 1700i

Ka (1b, M−1) 43 20 110 330 460 160 130 120 150 110 400 5200 140 2300 1400 240 90 60 150 6 120 260 8100 62 210 1000 3200 130 40 17 150 390 2100 130 420 69 130 61 2100 220 18 1900 750 9900

± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ±

13h 5h 10 88h 11i 43h 10 12 25h 2h 21j 590i 3h 140i 52i 25 10 6 20 1h 9h 5i 79i 5i 21 76i 300 6i 3i 8h 20i 2i 200i 10i 14i 13h 14i 14i 350j 29i 4i 100i 69i 1400i

a Details of how the volume (V), asphericity (ΩA), and surface area (S) were calculated are given in the Supporting Information. bMaximum solubilities in water (SMAX) and values of the logarithm of the octanol:water partition coefficient (log P) were taken from ref 50. cValues were calculated using Gaussian 03 at the RB3LYP/6-31G level of theory. dMiscible with water. eThe binding constants of these guests as determined by NMR titrations have been reported earlier.47 fThe binding constants of these guests were determined in the phosphate buffer (pH = 12). gValues were taken from ref 51. hThe binding constants of these guests were determined by NMR titrations and are listed in Table S3. iThe binding constants of these guests were determined by fluorescence titrations and are listed in Table S2. jThe binding constants of these guests were determined by ITC titrations and are listed in Table S4. All the titration experiments were repeated thrice, and the averaged values with standard deviations are reported here.

log P values lead to higher binding affinities. That is, hydrophobic guests show stronger binding than do hydrophilic molecules. Volume also affects the binding affinities but to a less extent. As shown in Figure 3, the result of the PCA performed on the data (eight variables: log Ka(1a), log Ka(1b), V, S, α, ΩA, μ, and log P) in Table 1 was drawn as a biplot diagram in which an eight-dimensional space was reduced to a two-dimensional

space with two principal components. The values of both the variables (loadings, blue lines) and the components (scores, red solid circles) for each guest are reported in the diagram, and the correlation coefficients are listed in Table S5. An overview regarding how these guest parameters affect the binding affinities may be obtained by closely examining this biplot diagram. 13469

DOI: 10.1021/jacs.8b09157 J. Am. Chem. Soc. 2018, 140, 13466−13477

Article

Journal of the American Chemical Society

Figure 3. PCA biplot diagram obtained by analyzing eight variables (log Ka(1a), log Ka(1b), V, S, α, ΩA, μ, and log P) for each guest. Both loadings (blue lines) and scores (red solid circles) are represented.

hydrophobic molecules with the appropriate size and shape. However, one should note that this rule is only valid among similar hydrophilic guests. The hydrogen-bonding ability was not included in these molecular descriptors. This rule cannot fully explain the binding behaviors of the molecular tubes which in many cases prefer hydrophilic molecules instead of the hydrophobic molecules with similar size and shape, for example, 1,4-dioxane (14000 M−1 for 1a) vs tetrahydropyran (68 M−1 for 1a), dimethoxymethane (6700 M−1 for 1a) vs diethyl ether (290 M−1 for 1a), and pyrazine (2200 M−1 for 1a) vs pyridine (590 M−1 for 1a). The dipole moment alone could not account for these drastic changes in binding strength. Therefore, further in-depth analyses and more control experiments are needed to paint a complete picture. The Role of the Cavity Water. The driving forces for the binding of hydrophilic molecules in several hosts with only a hydrophobic cavity11,12 were partly attributed to weak noncovalent interactions, such as C−H···π interactions. Is this reasonable? If yes, then why do our molecular tubes47 show an increased binding affinity to 1,4-dioxane when transferred from chloroform to water? Water is more polar than chloroform, and weaker interactions should be observed in water. If no, what is then the driving force? Further research is needed to answer this question. To confirm and test whether enhanced binding is limited to 1,4-dioxane, the binding constants of molecular tubes 2a and 2b with six representative guests (EA, DMC, benzene, MDO, 1,4-dioxane, and pyrazine) in CDCl3 were determined through NMR titrations (Figures S203−S220). All these guests have binding constants above 103 M−1 with these molecular tubes in water. However, their binding constants are substantially weaker in CDCl3. EA, DMC, and benzene show no detectable binding in CDCl3. While the binding constants of MDO, 1,4dioxane, and pyrazine decreased by ca. 2 orders of magnitude from water to chloroform (for 2a: MDO, 41 ± 1 M−1; 1,4-

The binding constants of 1a and 1b are highly correlated with each other (r = 0.862), suggesting the two cavities are rather similar in binding properties. Therefore, the binding data of 1a are mainly used in the following analysis. Generally, there is no strong correlation between the binding constants and these six molecular descriptors (V, S, α, ΩA, μ, and log P). This may be due to the high complexity of this binding system, and a single parameter is not sufficient to significantly affect the binding. However, weak correlations do exist: hydrophobicity (log P, r = 0.250 for 1a), volume (V, r = 0.221 for 1a), surface area (S, r = 0.297 for 1a), and polarizability (α, r = 0.225 for 1a) make positive contributions to the strength of guest binding; the dipole moment (r = −0.476 for 1a) makes a negative contribution, while asphericity (ΩA, r = 0.023 for 1a) is more or less irrelevant. The general trend is similar to that for the self-assembled host reported by Nitschke et al.12 This indicates these endo-functionalized molecular tubes prefer guests that have a large volume, large surface area, high hydrophobicity, and low dipole moment. Examples can be found by closely looking at the score plot in Figure 3 and the data in Table 1. Thirteen of the above-mentioned 14 guests with binding constants above 103 M−1 for 1a have smaller dipole moments than that of the other guests. The exception is tetrazole, for which the reason is still unknown. Additional hydrophobic groups usually lead to an increase in binding affinity, for example, 1,3-dioxolane (430 M−1 for 1a) vs MDO (22000 M−1 for 1a), benzene (1300 M−1 for 1a) vs toluene (9500 M−1 for 1a), and methanol (negligible for 1a, Figure S1) vs ethanol (9 M−1 for 1a). However, the prerequisite is that the guest cannot be too large; otherwise, the binding affinity will be significantly attenuated due to an incongruent fit with the cavity space, for example, dimethoxymethane (6700 M−1 for 1a) vs 2,2-dimethoxypropane (49 M−1 for 1a). This analysis provides a general picture of the binding ability of these molecular tubes: the molecular tubes prefer to bind 13470

DOI: 10.1021/jacs.8b09157 J. Am. Chem. Soc. 2018, 140, 13466−13477

Article

Journal of the American Chemical Society Table 2. Parameters of the Cavity Water of 1a and 1b Obtained from Molecular Dynamics Simulationa 1a hydrophobic region hydrogen-bonding region whole cavity

Nb

mc

1.5 3.0 4.5

1.67 3.20 (2.80)f 2.76

1b PCd

Nb

mc

PCd

48%

2.3e 1.5e 3.8

2.30 2.83 (1.77)f 2.67

35%

a

Only the results with the tip4pEW water model were shown here. For the results with the tip3p water model and the computational details, please see the Supporting Information. The average number of hydrogen bonds per H2O in bulk water is 3.62.54a bAverage number of the water molecules in this region or cavity. cAverage number of hydrogen bonds per H2O in this region or cavity (hydrogen bonds were considered if D−H···A ≤ 130° and D−A ≤ 3.5 Å). dPacking coefficient. eThis was obtained by summing up or averaging the numbers from two hydrophobic regions in 1b. fThe average number of hydrogen bonds excluding N−H···O hydrogen bonds is given in brackets.

dioxane, 62 ± 1 M−1; pyrazine, 11 ± 1 M−1; for 2b: MDO, 48 ± 3 M−1; 1,4-dioxane, 98 ± 5 M−1; and pyrazine, 13 ± 1 M−1). This result is surprising because synthetic hosts with hydrogenbonding sites and an ill-defined cavity show decreased binding affinities to hydrophilic molecules when transferred from nonpolar solvents to water.36 This suggests that water and thus the hydrophobic effect are very important even for the binding of hydrophilic molecules in water by these molecular tubes. How can the hydrophobic effect be important for recognizing hydrophilic molecules in water? We proposed earlier that the cavity water molecules are less hydrogenbonded than the bulk water.47 The release of these water molecules upon complexation is energetically favorable and thus responsible for the binding of hydrophilic molecules in a hydrophobic cavity.54 However, no evidence was provided in this earlier research. Herein, molecular dynamics (MD) simulations were performed to reveal the state of water molecules inside the cavities of 1a and 1b.54a The computed parameters are listed in Tables 2 and S9−S12, and representative snapshots are shown in Figure 4a. The average number of hydrogen bonds per water molecule in the whole cavity (2.76 for 1a and 2.67 for 1b) is markedly less than that of the bulk water (3.62).54a Therefore, the cavity water molecules are indeed poorly hydrogen-bonded. 1a contains ca. 4.5 water molecules in its cavity, while 1b has

ca. 3.8 water molecules. Once these cavity water molecules release to the bulk, more hydrogen bonds should be formed. Therefore, upon complexation, the release of these poorly hydrogen-bonded cavity water should be one major driving force for the binding of hydrophilic molecules. This explains why higher binding constants were achieved in water than in chloroform and why hydrophilic molecules can be bound inside a hydrophobic cavity.10,11 The cavities of 1a and 1b can be further divided into hydrogen-bonding regions (where the amide protons are located) and hydrophobic regions (Figure 4b and Figure S221). On average, each encapsulated water molecule in the hydrogen-bonding regions forms more hydrogen bonds than that in the hydrophobic regions. This is due to the contribution of the inward-directed amide protons. Glass et al. also reported another pair of molecular tubes by substituting the amide groups in 1a and 1b with allyl groups.55 These new molecular tubes have substantially stronger binding affinities (by 2 orders of magnitude) toward the hydrophobic tails of lipids than do the amide counterparts 1a and 1b. Replacing the inward-directed amide groups makes the encapsulated water even poorer in hydrogen bonding. Thus, the release of these cavity water is even more energetically favorable. This may explain the increased binding affinities. In other words, introducing polar binding sites into a hydrophobic cavity decreases its hydrophobicity and thus its binding strength to hydrophobic molecules/groups. However, it is favorable for the binding of hydrophilic molecules with complementary hydrogen-bonding sites. This may contribute to the reversed binding preference of hydrophilic molecules over hydrophobic molecules (see below). Water molecules in the cavity of 1b have a packing coefficient56 of 38% in contrast to the 48% for 1a, which are however less than the value of the bulk water (63%). The cavity of 1a is more well-defined and contains more poorly hydrogen-bonded cavity water molecules than 1b does. These may explain why 1a is generally a better binder than 1b and why the best binding until now was achieved with 1a rather than 1b (MDO, 22000 M−1 (1a) vs 8100 M−1 (1b) and 1,4dioxane, 14000 M−1 (1a) vs 3200 M−1 (1b)). This suggests that the size and shape of the hydrophobic cavity also play an important role in the binding affinity. The Role of Hydrogen Bonding. It was shown in Nitschke’s self-assembled cage with a hydrophobic cavity that hydrophobic guests are preferred over hydrophilic guests.12 For example, tetrahydropyran (740 M−1) is a better guest than 1,3-dioxane (330 M−1) and 1,4-dioxane (150 M−1). However, we reported earlier that 1,4-dioxane (14000 M−1) is a significantly better guest for 1a than are 1,3-dioxane (210 M−1) and tetrahydropyran (68 M−1).47 At first glance, the

Figure 4. a) Representative MD snapshots of 1a and 1b in water. Green dotted lines indicate the hydrogen bonds involving the encapsulated water molecules. b) Hydrophobic and hydrogenbonding regions of 1a. 13471

DOI: 10.1021/jacs.8b09157 J. Am. Chem. Soc. 2018, 140, 13466−13477

Article

Journal of the American Chemical Society

Figure 5. Single-crystal X-ray structures of the complexes of (a) 2a or (b) 2b with the guests. CCDC numbers for these crystals: 1821196− 1821200, 1821203−1821212.

N−H···π interactions between DMF’s π system and the host’s amide proton exist instead of a N−H···O hydrogen bond. This may be because forming a hydrogen bond with the carbonyl oxygen atom of DMF would render the large methyl groups too close to the other amide proton of the host and cause steric repulsion. While for N-ethylacetamide, the host’s NH proton is hydrogen bonded to the guest’s carbonyl oxygen atom, but the steric repulsion between the two amide protons from the host and the guest makes the guest’s amide bond twisted in a nonplanar form. This is quite unusual for an amide bond. Two short and assumably strong hydrogen bonds were detected for guests 1,4-dioxane,47 pyrazine, MDO, DMC, and dimethoxymethane. The two hydrogen-bonding donors in the cavities of the molecular tubes are both satisfied. All these guests have binding constants above 103 M−1. These perfectly matched hydrogen bonds should be responsible for the high binding affinities as well as the high selectivity over the molecules with a similar size and shape but incongruent hydrogen-bonding acceptors. In the absence of guests, the cavity water molecules also form hydrogen bonds with the amide protons in the cavities of 1a and 1b. This is supported by MD simulations (Figure 4a). Binding to a guest would replace these water molecules and destroy these hydrogen bonds. However, new hydrogen bonds with guest molecules will form. Either stronger or weaker hydrogen bonds are formed when replacing the cavity water molecules with a guest. Therefore, this situation must be considered when one evaluates the contribution of hydrogen bonding to the guest binding. NMR titrations of molecular

reversed binding selectivity was surprising and puzzling. Compared to the hosts with only a hydrophobic cavity, these molecular tubes have additional hydrogen-bonding sites in their hydrophobic cavity. We therefore proposed earlier47 that hydrogen bonding should be responsible for this surprising behavior and the high selectivity to the highly hydrophilic molecule 1,4-dioxane. However, is hydrogen bonding significantly involved in the binding of all the other hydrophilic molecules? These hydrogen bonds are shielded in a hydrophobic cavity and do not interfere with the bulk water molecules. Therefore, crystal structures may help reveal the role of hydrogen bonding in the binding of all these hydrophilic molecules. The slow evaporation of the CH2Cl2 solutions of 2a or 2b in the presence of guests affords single crystals of the complexes, which are suitable for X-ray crystallography. The single-crystal structures are shown in Figure 5. None of these crystals were grown from water, but the binding mode of guests in the host cavity should be similar to that in water since most of the solvent molecules in the cavity are anyway expelled upon guest encapsulation. Single hydrogen bond (for guests EA, DEA, and Nethylacetamide) or N−H···π interactions (for guests DMF, pyrrole, thiazole, imidazole, and benzene) were detected in most of these crystal structures. These guests usually show decent binding affinities (102−103 M−1 for most of them) with 1a or 1b in water. It is noted that some of these guests contain hydrogen-bonding acceptors, but no hydrogen bonds are formed. This may be due to steric reasons. For example, weak 13472

DOI: 10.1021/jacs.8b09157 J. Am. Chem. Soc. 2018, 140, 13466−13477

Article

Journal of the American Chemical Society

hydrophobic effect, but the hydrophobic group should not cause any steric repulsion with the cavity walls. Thermodynamic Consequence. How would this unique combination of hydrogen-bonding sites with a hydrophobic cavity affect the binding thermodynamics? The thermodynamic parameters (ΔG°, ΔH°, and ΔS°) for 18 guests were available from ITC titrations (Figure 7, and Table S4). For

tube 1a in H2O:D2O (9:1) by six representative guests (EA, DMC, benzene, MDO, 1,4-dioxane, and pyrazine) were performed (Figures S223−S228) to reveal this effect. Since H2O is dominant in this solution, the amide protons are not fully exchanged with deuterium and thus visible in NMR spectra. The chemical shifts of the amide protons can be used as an estimate of the changes in the strength of these hydrogen bonds.57 Stronger hydrogen bonds would cause a downfield shift of these proton signals, while an upfield shift suggests the formation of weaker hydrogen bonds. As shown in Figure 6, the amide protons shift downfield with the addition of MDO, 1,4-dioxane, or pyrazine, suggesting

Figure 7. Thermodynamic data (ΔH° and TΔS°) for the binding of neutral guests with 1a (red dots) or 1b (blue triangles).

most of the guests, the binding is generally dominated by the enthalpic contribution with a slightly favorable or unfavorable entropic contribution. However, for benzene and toluene, both enthalpy and entropy are favorable. The cavity water molecules are less hydrogen-bonded than the bulk water molecules. Upon release into the bulk, more hydrogen bonds should be formed. This was suggested to be the cause of the dominated enthalpy in the guest binding of a well-defined hydrophobic cavity.54a Nevertheless, a more complex mechanism is also possible.43 Weak enthalpy−entropy compensation (Figure 7) is observed for the binding data of both 1a and 1b. These data are more scattered than those of other synthetic systems.60 With respect to the cause, not only the tightness of the binding should be considered but also the entropy of solvation may change. The release of cavity water from the hosts and the desolvation of hydrophilic guests could further complicate the situation. For different hydrophilic guest molecules, their solvation and desolvation are different, which is rather distinct from hydrophobic molecules. Together, this may explain why the enthalpy−entropy compensation is weak here. In addition, the strong enthalpy−entropy compensation in other supramolecular systems was suggested to be possibly due to large uncertainties in the enthalpy measurement.61 That is, a weak enthalpy−entropy compensation should be more common. A clear solvent isotope effect was also observed. ITC titrations of 1a and 1b by the six representative guests (EA, DMC, benzene, MDO, 1,4-dioxane, and pyrazine) in both H2O and D2O were performed (Tables 3 and S13, and Figures S229−S252). Generally, solvent isotope effects (KD/KH = 1.1− 1.5, ΔΔGD/H = −0.2 to −1.0 kJ/mol) were observed for all the cases. However, for the guests (EA, DMC, and benzene) that form weaker hydrogen bonds than those of the encapsulated water molecules, the binding is less exothermic in D2O than in H2O (ΔΔHD/H = 0−4.7 kJ/mol), and a large enthalpy− entropy compensation was seen; while for the guests forming

Figure 6. Changes of the chemical shift of the amide protons of 1a (0.5 mM, H2O/D2O = 9:1) with the addition of guests EA, DMC, benzene, MDO, 1,4-dioxane, or pyrazine.

the formation of stronger hydrogen bonds than those with the encapsulated water molecules. However, with guests EA, DMC, or benzene,58 the amide protons shift upfield, indicating the formation of weaker hydrogen bonds with these guests. That is, from the point of view of hydrogen bonding, the binding is favorable for guests MDO, 1,4-dioxane, and pyrazine but unfavorable for guests EA, DMC, and benzene. This is in line with the binding studies in CDCl3. Considering the high desolvation penalty of EA and DMC, the binding is even more unfavorable. Consequently, the binding to EA and DMC should mainly come from the hydrophobic effect, that is, the release of the cavity water. While for guests MDO, 1,4-dioxane, and pyrazine, hydrogen bonding and the hydrophobic effect through releasing the cavity water both contribute significantly. This explains why a high binding selectivity was achieved for these guests over that of the other guests with similar size and shape. This is consistent with the statement by Biedermann et al. on how to design a high-affinity receptor with high selectivity.59 To summarize the findings discussed in the above three sections, we conclude a good guest for these endo-functionalized molecular tubes should meet the following criteria: a) the guest should contain two hydrogen-bonding acceptor atoms that are located in the appropriate orientation and distance to complement the hydrogen-bonding donors in the cavities. An additional heteroatom is better eliminated to minimize the dipole moment; b) when the hydrogen-bonding sites are satisfied, a large hydrophobic group is needed to expel all the cavity water and to maximize the contribution of 13473

DOI: 10.1021/jacs.8b09157 J. Am. Chem. Soc. 2018, 140, 13466−13477

Article

Journal of the American Chemical Society

Table 3. ΔG°, ΔH°, and −TΔS° (kJ mol−1) of the Complexes of 1a (Left) and 1b (Right) with the Six Representative Guestsa 1a EA DMC benzene MDO 1,4-dioxane pyrazine

H2O D2O H2O D2O H2O D2O H2O D2O H2O D2O H2O D2O

1b

ΔG°

ΔH°

−TΔS°

ΔG°

ΔH°

−TΔS°

−20.2 −21.2 −18.7 −19.5 −17.7 −18.6 −24.3 −25.0 −23.5 −24.1 −19.8 −20.2

−21.6 −20.5 −24.0 −21.7 −12.2 −7.5 −25.0 −25.5 −25.6 −26.5 −19.9 −20.4

1.4 −0.7 5.3 2.2 −5.5 −11.1 0.7 0.5 2.1 2.4 0.1 0.2

−18.9 −19.9 −17.8 −18.6 −17.1 −17.3 −21.9 −22.8 −20.6 −21.7 −18.4 −18.7

−26.1 −23.3 −23.9 −20.0 −10.1 −10.2 −23.8 −24.1 −25.5 −26.2 −18.1 −21.1

7.2 3.4 6.1 1.4 −7.0 −7.1 1.9 1.3 4.9 4.5 −0.3 2.4

a

The data were determined by ITC titrations in D2O or H2O.

Table 4. Association Constants Ka (M−1) and Thermodynamic Parameters (ΔG°, ΔH°, and −TΔS°; kJ mol−1) of Molecular Tubes 1a and 1b with 3 and 4 in H2O at 25 °C as Determined by Fluorescence (FL) and ITC Titrations

stronger hydrogen bonds, the opposite was observed (ΔΔHD/H = −3.0 to −0.3 kJ/mol), and the magnitude of enthalpy and entropy changes is substantially smaller. The solvent isotope effect for these molecular tubes is larger when compared to that of cyclodextrins (KD/KH ≤ 1.25)62 and smaller when compared to that of cucurbiturils (KD/KH ≥ 1.5).59 Although the same cavities are involved, and similar cavity water molecules are released, the hydrogen-bonding strength among the cavity water is different in D2O and H2O. It is wellknown that the zero-point vibrational energy of the O−D bonds of D2O is lower than for the O−H bonds of H2O. Therefore, D2O forms stronger hydrogen bonds than does H2O.63 In addition, the amide protons in the host cavities are labile and would be fully exchanged to deuterium in D2O. The hydrogen bonds between host and guest would also be influenced through the H/D exchange of the amide protons. All these factors together contribute to the observed solvent isotope effect and explain the different thermodynamic origins of the two groups of guests. The “Best-Fit” Guest. Following the criteria for the binding with these molecular tubes, a “best-fit” guest was searched for. From the crystal structures (Figure 4), these hydrophilic guests only occupy the hydrogen-bonding region and part of the hydrophobic region of the cavity. Additional water may coexist with these guest molecules in the cavity.49 Therefore, a higher binding affinity may be achieved with a guest to fully occupy the cavity space and thus expel all the cavity water and to simultaneously satisfy the hydrogenbonding sites.48 1,4-Dioxane and MDO with complementary hydrogen-bonding sites are the only two guests whose binding constants are larger than 104 M−1 (with 1a). Therefore, these two guests were used as the structural basis to search for the “best-fit” guest. Since the unfilled region is hydrophobic, the fusion of a hydrophobic group to these guests was attempted to improve the binding affinities. Thus, guests 3 and 4 (Table 4) were selected, and their binding constants were determined (Figures S253−S262 and Table 4). However, guest 3 does not give a significantly higher binding affinity, presumably due to steric hindrance with the cavity walls. In contrast, guest 4 achieved a binding constant of 106 M−1 with 1a, which is better than that of the epoxide 2-methyl-3-phenyloxirane (105 M−1) reported earlier.48 This binding constant is rather high when one considers that the guest is small and contains a hydrophilic acetal group. The binding is again mainly contributed by

FL

ITC

Ka × 10

4

3 4

1a 1b 1a 1b

1.6 2.9 220 51

± ± ± ±

0.3 0.1 18 0.1

4

Ka × 10

ΔG°

ΔH°

−TΔS°

± ± ± ±

0.1 0.1 0.1 2.7

−23.1 −23.5 −35.2 −32.4

−38.8 −38.7 −45.1 −52.7

15.7 15.2 9.9 20.3

1.1 1.3 140 46

enthalpy with an unfavorable entropy. Guest 4 not only satisfies the hydrogen-bonding sites in the cavity of 1a but also structurally complements the hydrophobic region and expels all the cavity water molecules. From another point of view, guest 4 can be considered to be obtained by attaching an additional polar acetal group to benzene. Competition experiments (Figure S263) between 4 and benzene (Ka = 1300 M−1, for 1a) in the presence of 1a indicate that 1a overwhelmingly prefers guest 4 over benzene. This suggests that the additional polar acetal group is crucial for the high binding affinity of 4. In addition, this indicates that these endo-functionalized molecular tubes are also good receptors for molecules with both hydrophilic and hydrophobic groups, which is the common structural feature for the majority of organic molecules. This also suggests strong binding for synthetic receptors is not limited to highly hydrophobic guests,3b and a well-solubilized molecule may also provide nanomolar binding provided a synthetic host with complementary cavity size and inward-directed binding sites is available.

3. CONCLUSIONS In this Article, we reported a systematic study on the binding behavior of the endo-functionalized molecular tubes to 44 hydrophilic molecules in water by using NMR, fluorescence, and ITC titrations, X-ray crystallography, and molecular dynamics simulations. Principal component analysis reveals that the binding strength is weakly correlated to the hydrophobicity, volume, surface area, and dipole moment of 13474

DOI: 10.1021/jacs.8b09157 J. Am. Chem. Soc. 2018, 140, 13466−13477

Article

Journal of the American Chemical Society guests. Furthermore, the hydrophobic effect through the release of the poorly hydrogen-bonded water molecules in a well-defined cavity is the major driving force. Hydrogen bonding is the key factor in determining the high binding selectivity to certain hydrophilic molecules over other molecules with similar size and shape: a) introducing hydrogen-bonding sites make the cavity less hydrophobic and thus decreases the binding affinity toward hydrophobic guests with similar size and shape; b) among hydrophilic molecules, only the ones with complementary hydrogenbonding sites form stronger hydrogen bonds than does the encapsulated water. Therefore, weak noncovalent interactions, such as C−H···π interactions, should not contribute much to the binding of hydrophilic molecules in water. These findings clarify the role of the hydrophobic effect and noncovalent interactions in the binding of hydrophilic molecules and echo the computational research on protein active sites.64 The unique feature of combining hydrogen-bonding sites into a hydrophobic cavity also leads to thermodynamic consequences: a dominated enthalpy, a weak enthalpy− entropy compensation, and a positive solvent isotope effect. In addition, a good guest for these molecular tubes should meet the following criteria: a) the heteroatoms in the guest should be complementary to the hydrogen-bonding sites in the cavity in terms of orientation, distance, and number; b) on this basis, the hydrophobic group of the guest should be large enough to expel all the cavity water molecules and maximize the hydrophobic effect but not too large to cause any steric repulsion with the cavity walls; the dipole moment should be low. Following these guidelines, a binding constant as high as 106 M−1 was achieved with 2-phenyl-1,3-dioxolane (4). This suggests extremely strong binding with synthetic receptors is not limited to highly hydrophobic guest molecules3b but may also be achieved with hydrophilic molecules. Vice versa, the experience gained here can be used to design an efficient and selective receptor for a certain biologically or environmentally relevant hydrophilic molecule. First, a deep hydrophobic cavity with complementary size and shape is needed. Second, the heteroatoms or polar groups in the hydrophilic molecule should be satisfied by hydrogen-bonding sites or other polar binding sites. More complementary hydrogen bonds or polar interactions would lead to a higher selectivity. Last, these polar interactions should be located inside the deep hydrophobic cavity to experience a nonpolar environment and avoid water competition. The endo-functionalized molecular tubes presented here are unprecedented in terms of binding ability and binding mechanism and are complementary to conventional macrocyclic receptors. Therefore, we believe the present research lays a basis for their further applications, for example, in sensing,47,48 catalysis, self-assembly, molecular machinery, and materials science. Moreover, the cavities of these molecular tubes closely resemble protein binding pockets. The conclusions obtained here should also be helpful in understanding molecular recognition in complex biological systems.





Experimental details, 1H NMR spectra of the complexes, computational details, NMR titration data, fluorescence titration data, ITC titration data, and single-crystal X-ray data (PDF) Crystallographic data for pyrazine@2a (CIF) Crystallographic data for DMF@2a (CIF) Crystallographic data for pyrrole@2a (CIF) Crystallographic data for thiazole@2a (CIF) Crystallographic data for imidazole@2a (CIF) Crystallographic data for benzene@2a (CIF) Crystallographic data for EA@2b (CIF) Crystallographic data for N-ethylacetamide@2b (CIF) Crystallographic data for DMC@2b (CIF) Crystallographic data for dimethoxymethane@2b (CIF) Crystallographic data for DEA@2b (CIF) Crystallographic data for MDO@2b (CIF) Crystallographic data for 1,4-dioxene@2b (CIF) Crystallographic data for pyrazine@2b (CIF) Crystallographic data for pyrrole@2b (CIF)

AUTHOR INFORMATION

Corresponding Author

*[email protected]. ORCID

Georg Schreckenbach: 0000-0002-4614-0901 Wei Jiang: 0000-0001-7683-5811 Author Contributions ¶

H.Y., H.K., and X.-B.Z. contributed equally.

Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS This research was financially supported by the National Natural Science Foundation of China (Nos. 21572097, 21772083, 21822104; W.J.), the Shenzhen Special Funds (KQJSCX20170728162528382, JCYJ20170307105848463; W.J.), the Shenzhen Nobel Prize Scientists Laboratory Project (C17213101), and the Natural Sciences and Engineering Research Council of Canada (NSERC, Discovery Grant; G.S.). We thank SUSTC-MCPC for the support with instruments.



REFERENCES

(1) (a) Molecular Recognition. In Supramolecular Chemistry: From Molecules to Nanomaterials; Gale, P. A., Steed, J. W., Eds.; Wiley: Chichester, U.K., 2012; Vol. 3. (b) Chatterji, D. Basics of Molecular Recognition; CRC Press: Boca Raton, FL, 2016. (2) (a) Oshovsky, G. V.; Reinhoudt, D. N.; Verboom, W. Angew. Chem., Int. Ed. 2007, 46, 2366. (b) Kataev, E. A.; Muller, C. Tetrahedron 2014, 70, 137. (3) (a) Murray, J.; Kim, K.; Ogoshi, T.; Yao, W.; Gibb, B. C. Chem. Soc. Rev. 2017, 46, 2479. (b) Liu, W.; Samanta, S. K.; Smith, B. D.; Isaacs, L. Chem. Soc. Rev. 2017, 46, 2391. (c) Chen, Y.; Huang, F.; Li, Z.-T.; Liu, Y. Sci. China: Chem. 2018, 61, 979. (4) (a) Rekharsky, M. V.; Inoue, Y. Chem. Rev. 1998, 98, 1875. (b) Liu, Y.; Chen, Y. Acc. Chem. Res. 2006, 39, 681. (c) Crini, G. Chem. Rev. 2014, 114, 10940. (5) (a) Lee, J. W.; Samal, S.; Selvapalam, N.; Kim, H.-J.; Kim, K. Acc. Chem. Res. 2003, 36, 621. (b) Lagona, J.; Mukhopadhyay, P.; Chakrabarti, S.; Isaacs, L. Angew. Chem., Int. Ed. 2005, 44, 4844. (c) Urbach, A.; Ramalingam, V. Isr. J. Chem. 2011, 51, 664. (d) Assaf, K. I.; Nau, W. M. Chem. Soc. Rev. 2015, 44, 394. (e) Barrow, S. J.; Kasera, S.; Rowland, M. J.; del Barrio, J.; Scherman, O. A. Chem. Rev. 2015, 115, 12320.

ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/jacs.8b09157. 13475

DOI: 10.1021/jacs.8b09157 J. Am. Chem. Soc. 2018, 140, 13466−13477

Article

Journal of the American Chemical Society (6) (a) Jordan, J. H.; Gibb, B. C. Chem. Soc. Rev. 2015, 44, 547. (b) Yoshizawa, M.; Klosterman, J. K.; Fujita, M. Angew. Chem., Int. Ed. 2009, 48, 3418. (c) Brown, C. J.; Toste, F. D.; Bergman, R. G.; Raymond, K. N. Chem. Rev. 2015, 115, 3012. (d) McConnell, A. J.; Wood, C. S.; Neelakandan, P. P.; Nitschke, J. R. Chem. Rev. 2015, 115, 7729. (7) Hill, D. J.; Mio, M. J.; Prince, R. B.; Hughes, T. S.; Moore, J. S. Chem. Rev. 2001, 101, 3893. (8) Klärner, F.-G.; Schrader, T. Acc. Chem. Res. 2013, 46, 967. (9) (a) Schneider, H.-J.; Yatsimirsky, A. K. Chem. Soc. Rev. 2008, 37, 263. (b) Synthetic Receptors for Biomolecules: Design Principles and Applications; Smith, B. D., Ed.; The Royal Society of Chemistry: Cambridge, U.K., 2015. (c) Lee, H. H. L.; Lee, J. W.; Jang, Y.; Ko, Y. H.; Kim, K.; Kim, H. I. Angew. Chem., Int. Ed. 2016, 55, 8249. (d) Zhao, Y. ChemPhysChem 2013, 14, 3878. (e) Davis, A. P.; Kubik, S.; Cort, A. D. Org. Biomol. Chem. 2015, 13, 2499. (10) (a) Kato, Y.; Conn, M. M.; Rebek, J., Jr. J. Am. Chem. Soc. 1994, 116, 3279. (b) Torneiro, M.; Still, W. C. J. Am. Chem. Soc. 1995, 117, 5887. (c) Schmuck, C. Synlett 2011, 2011, 1798. (d) Ferrand, Y.; Crump, M. P.; Davis, A. P. Science 2007, 318, 619. (e) Mooibroek, T. J.; Casas-Solvas, J. M.; Harniman, R. L.; Renney, C. M.; Carter, T. S.; Crump, M. P.; Davis, A. P. Nat. Chem. 2016, 8, 69. (11) (a) Biedermann, F.; Uzunova, V. D.; Scherman, O. A.; Nau, W. M.; De Simone, A. J. Am. Chem. Soc. 2012, 134, 15318. (b) Yamashina, M.; Akita, M.; Hasegawa, T.; Hayashi, S.; Yoshizawa, M. Sci. Adv. 2017, 3, e1701126. (c) Kusaba, S.; Yamashina, M.; Akita, M.; Kikuchi, T.; Yoshizawa, M. Angew. Chem., Int. Ed. 2018, 57, 3706. (d) Kobayashi, K.; Asakawa, Y.; Kato, Y.; Aoyama, Y. J. Am. Chem. Soc. 1992, 114, 10307. (12) Smulders, M. M. J.; Zarra, S.; Nitschke, J. R. J. Am. Chem. Soc. 2013, 135, 7039. (13) Pedersen, C. J. J. Am. Chem. Soc. 1967, 89, 7017. (14) Liu, Z.; Nalluri, S. K. M.; Stoddart, F. J. Chem. Soc. Rev. 2017, 46, 2459. (15) Gokel, G. W.; Leevy, M.; Weber, M. E. Chem. Rev. 2004, 104, 2723. (16) Calixarenes and Beyond; Neri, P., Sessler, J. L., Wang, M.-X., Eds.; Springer: Heidelberg, 2016; DOI: 10.1007/978-3-319-31867-7. (17) Hardie, M. J. Chem. Soc. Rev. 2010, 39, 516. (18) Gale, P. A.; Anzenbacher, P., Jr.; Sessler, J. L. Coord. Chem. Rev. 2001, 222, 57. (19) (a) Odell, B.; Reddington, M. V.; Slawin, A. M. Z.; Spencer, N.; Stoddart, J. F.; Williams, D. J. Angew. Chem., Int. Ed. Engl. 1988, 27, 1547. (b) Dale, E. J.; Vermeulen, N. A.; Juríček, M.; Barnes, J. C.; Young, R. M.; Wasielewski, M. R.; Stoddart, J. F. Acc. Chem. Res. 2016, 49, 262. (20) Wang, M.-X. Acc. Chem. Res. 2012, 45, 182. (21) Ogoshi, T.; Yamagishi, T.-a.; Nakamoto, Y. Chem. Rev. 2016, 116, 7937. (22) Rambo, B. M.; Gong, H.-Y.; Oh, M.; Sessler, J. L. Acc. Chem. Res. 2012, 45, 1390. (23) Svec, J.; Necas, M.; Sindelar, V. Angew. Chem., Int. Ed. 2010, 49, 2378. (24) Lee, S.; Chen, C.-H.; Flood, A. H. Nat. Chem. 2013, 5, 704. (25) (a) Guo, Q.-H.; Fu, Z.-D.; Zhao, L.; Wang, M.-X. Angew. Chem., Int. Ed. 2014, 53, 13548. (b) Wang, M.-X. Sci. China: Chem. 2018, 61, 993. (26) Zhang, G.-W.; Li, P.-F.; Meng, Z.; Wang, H.-X.; Han, Y.; Chen, C.-F. Angew. Chem., Int. Ed. 2016, 55, 5304. (27) Chen, H.; Fan, J.; Hu, X.; Ma, J.; Wang, S.; Li, J.; Yu, Y.; Jia, X.; Li, C. Chem. Sci. 2015, 6, 197. (28) Jia, F.; He, Z.; Yang, L.-P.; Pan, Z.-S.; Yi, M.; Jiang, R.-W.; Jiang, W. Chem. Sci. 2015, 6, 6731. (29) Yu, G.; Jie, K.; Huang, F. Chem. Rev. 2015, 115, 7240. (30) (a) Wenz, G.; Han, B.-H.; Müller, A. Chem. Rev. 2006, 106, 782. (b) Xue, M.; Yang, Y.; Chi, X.; Yan, X.; Huang, F. Chem. Rev. 2015, 115, 7398. (31) (a) Balzani, V.; Credi, A.; Raymo, F. M.; Stoddart, J. F. Angew. Chem., Int. Ed. 2000, 39, 3348. (b) Kay, E. R.; Leigh, D. A.; Zerbetto,

F. Angew. Chem., Int. Ed. 2007, 46, 72. (c) Erbas-Cakmak, S.; Leigh, D. A.; McTernan, C. T.; Nussbaumer, A. L. Chem. Rev. 2015, 115, 10081. (32) Breslow, R. Acc. Chem. Res. 1995, 28, 146. (33) de Silva, A. P.; Gunaratne, H. Q. N.; Gunnlaugsson, T.; Huxley, A. J. M.; McCoy, C. P.; Rademacher, J. T.; Rice, T. E. Chem. Rev. 1997, 97, 1515. (34) (a) Zheng, B.; Wang, F.; Dong, S.; Huang, F. Chem. Soc. Rev. 2012, 41, 1621. (b) Guo, D.-S.; Liu, Y. Chem. Soc. Rev. 2012, 41, 5907. (c) Avestro, A.-J.; Belowich, M. E.; Stoddart, J. F. Chem. Soc. Rev. 2012, 41, 5881. (d) Harada, A.; Takashima, Y.; Nakahata, M. Acc. Chem. Res. 2014, 47, 2128. (e) Qu, D.-H.; Wang, Q.-C.; Zhang, Q.W.; Ma, X.; Tian, H. Chem. Rev. 2015, 115, 7543. (f) Ma, X.; Zhao, Y. Chem. Rev. 2015, 115, 7794. (35) For one representative example, see: Chandramouli, N.; Ferrand, Y.; Lautrette, G.; Kauffmann, B.; Mackereth, C. D.; Laguerre, M.; Dubreuil, D.; Huc, I. Nat. Chem. 2015, 7, 334. (36) For two examples in which the binding affinities decrease by 4 orders of magnitude when transferred from chloroform to water, see: (a) Allott, C.; Adams, H.; Hunter, C. A.; Thomas, J. A.; Bernad, P. L., Jr.; Rotger, C. Chem. Commun. 1998, 2449. (b) Klein, E.; Ferrand, Y.; Barwell, N. P.; Davis, A. P. Angew. Chem., Int. Ed. 2008, 47, 2693. (37) Anslyn, E. V.; Dougherty, D. A. Modern Physical Organic Chemistry; University Science Books: Sausalito, 2006; pp 230−232. (38) Weber, P. C.; Ohlendorf, D. H.; Wendoloski, J. J.; Salemme, F. Science 1989, 243, 85. (39) Avidin-Biotin Interactions: Methods and Applications; McMahon, R. J., Ed.; Humana Press: Totowa, NJ, 2008; DOI: 10.1007/978-159745-579-4. (40) Vyas, N. K.; Vyas, M. N.; Quiocho, F. A. Science 1988, 242, 1290. (41) (a) Butterfield, S. M.; Rebek, J., Jr. J. Am. Chem. Soc. 2006, 128, 15366. (b) Verdejo, B.; Gil-Ramírez, G.; Ballester, P. J. Am. Chem. Soc. 2009, 131, 3178. (c) Escobar, L.; Díaz-Moscoso, A.; Ballester, P. Chem. Sci. 2018, 9, 7186. (42) For an example on hydrogen-bonded assembly shielded in a hydrophobic cavity, see: Sawada, T.; Yoshizawa, M.; Sato, S.; Fujita, M. Nat. Chem. 2009, 1, 53. (43) Cremer, P. S.; Flood, A. H.; Gibb, B. C.; Mobley, D. L. Nat. Chem. 2018, 10, 8. (44) For an early water-soluble cyclophane with an inwardly directed ketone group, see: (a) Carcanague, D. R.; Knobler, C. B.; Diederich, F. J. Am. Chem. Soc. 1992, 114, 1515. For nonwater-soluble capsules with functionalized interiors, see: (b) Adriaenssens, L.; Ballester, P. Chem. Soc. Rev. 2013, 42, 3261. (45) (a) Huang, G.-B.; He, Z.; Cai, C.-X.; Pan, F.; Yang, D.; Rissanen, K.; Jiang, W. Chem. Commun. 2015, 51, 15490. (b) Huang, G.-B.; Valkonen, A.; Rissanen, K.; Jiang, W. Chem. Commun. 2016, 52, 9078. (46) Shorthill, B. J.; Avetta, C. T.; Glass, T. E. J. Am. Chem. Soc. 2004, 126, 12732. (47) Huang, G.-B.; Wang, S.-H.; Ke, H.; Yang, L.-P.; Jiang, W. J. Am. Chem. Soc. 2016, 138, 14550. (48) Wang, L.-L.; Chen, Z.; Liu, W.-E.; Ke, H.; Wang, S.-H.; Jiang, W. J. Am. Chem. Soc. 2017, 139, 8436. (49) Paul, R.; Paul, S. Phys. Chem. Chem. Phys. 2018, 20, 16540. (50) (a) CRC Handbook of Chemistry and Physics, 93rd ed.; Haynes, W. M., Ed.; CRC Press: Boca Raton, FL, 2012. (b) CRC Handbook of Chemistry and Physics, Internet Version 2005; Lide, D. R., Ed.; CRC Press: Boca Raton, FL, 2005; http://hbcponline.com (accessed Sept 26, 2018). (c) CRC Handbook of physical properties of organic chemicals; Howard, P. H., Meylan, W. M., Eds.; CRC Press: Boca Raton, FL, 1997. (51) https://scifinder.cas.org/scifinder/view/scifinder/ scifinderExplore.jsf (accessed Sept 26, 2018). (52) Turnbull, W. B.; Daranas, A. H. J. Am. Chem. Soc. 2003, 125, 14859. (53) Jurs, P. C.; Bakken, G. A.; McClelland, H. E. Chem. Rev. 2000, 100, 2649. 13476

DOI: 10.1021/jacs.8b09157 J. Am. Chem. Soc. 2018, 140, 13466−13477

Article

Journal of the American Chemical Society (54) (a) Biedermann, F.; Nau, W. M.; Schneider, H.-J. Angew. Chem., Int. Ed. 2014, 53, 11158. (b) Persch, E.; Dumele, O.; Diederich, F. Angew. Chem., Int. Ed. 2015, 54, 3290. (c) Hillyer, M. B.; Gibb, B. C. Annu. Rev. Phys. Chem. 2016, 67, 307. (d) Metherell, A. J.; Cullen, W.; Williams, N. H.; Ward, M. D. Chem. - Eur. J. 2018, 24, 1554. (55) Avetta, C. T.; Shorthill, B. J.; Ren, C.; Glass, T. E. J. Org. Chem. 2012, 77, 851. (56) Mecozzi, S.; Rebek, J., Jr. Chem. - Eur. J. 1998, 4, 1016. (57) Ajami, D.; Tolstoy, P. M.; Dube, H.; Odermatt, S.; Koeppe, B.; Guo, J.; Limbach, H.-H.; Rebek, J., Jr. Angew. Chem., Int. Ed. 2011, 50, 528. (58) For benzene, the significant upfield shift of the NH protons may be partly due to the shielding effect of benzene, because they are located on the top of benzene and form N-H···π interactions (see the crystal structure in Figure 5). (59) Biedermann, F.; Vendruscolo, M.; Scherman, O. A.; De Simone, A.; Nau, W. M. J. Am. Chem. Soc. 2013, 135, 14879. (60) (a) Houk, K. N.; Leach, A. G.; Kim, S. P.; Zhang, X. Angew. Chem., Int. Ed. 2003, 42, 4872. (b) Leung, D. H.; Bergman, R. G.; Raymond, K. N. J. Am. Chem. Soc. 2008, 130, 2798. (61) Chodera, J. D.; Mobley, D. L. Annu. Rev. Biophys. 2013, 42, 121. (62) Rekharsky, M. V.; Inoue, Y. J. Am. Chem. Soc. 2002, 124, 12361. (63) (a) Chervenak, M. C.; Toone, E. J. J. Am. Chem. Soc. 1994, 116, 10533. (b) Lopez, M. M.; Makhatadze, G. I. Biophys. Chem. 1998, 74, 117. (64) Young, T.; Abel, R.; Kim, B.; Berne, B. J.; Friesner, R. A. Proc. Natl. Acad. Sci. U. S. A. 2007, 104, 808.

13477

DOI: 10.1021/jacs.8b09157 J. Am. Chem. Soc. 2018, 140, 13466−13477