Calorimetric Analysis of the 1:1 Complexes Formed between a Water-soluble Deep-cavity Cavitand, and Cyclic and Acyclic Carboxylic Acids

Guest/Host Complexes of Octa Acid and Amphiphilic Benzylidene-3-methylimidazolidinones Exchange Hosts within the NMR Time Scale

Cavitand octa acid (OA) is established to form a stable capsular assembly with one or two hydrophobic guest molecules (1:2 or 2:2 guest/host complex). Examples are known in which the guest molecule tumbles within the capsule without disrupting the structure of the capsuleplex. This process makes the two OA molecules that form the capsule magnetically equivalent. In this study, we have examined the dynamics of capsules that host amphiphilic benzylidene-3-methylimidazolidinone molecules as guests. In these capsuleplexes, although the guest does not tumble, the two OA molecules become magnetically equivalent because the two OA molecules that form the capsule exchange their positions in the NMR time scale. This is equivalent to the content of the capsule remaining stationary while the capsule swirls around it. Benzylidene-3-methylimidazolidinones form both 1:1 and 1:2 supramolecular complexes with cavitand OA. Two-dimensional NMR, ROESY, and NOESY data suggest that in a 300 ms time scale, the two halves of the capsule exchange between themselves and with free OA. The conclusion drawn here provides valuable information concerning the stability of the OA capsuleplex and cavitandplex that is used as the well-defined space to control the excited-state chemistry and dynamics of confined guest molecules.

Binding free energies in the SAMPL6 octa-acid host-guest challenge calculated with MM and QM methods

We have estimated free energies for the binding of eight carboxylate ligands to two variants of the octa-acid deep-cavity host in the SAMPL6 blind-test challenge (with or without endo methyl groups on the four upper-rim benzoate groups, OAM and OAH, respectively). We employed free-energy perturbation (FEP) for relative binding energies at the molecular mechanics (MM) and the combined quantum mechanical (QM) and MM (QM/MM) levels, the latter obtained with the reference-potential approach with QM/MM sampling for the MM → QM/MM FEP. The semiempirical QM method PM6-DH+ was employed for the ligand in the latter calculations. Moreover, binding free energies were also estimated from QM/MM optimised structures, combined with COSMO-RS estimates of the solvation energy and thermostatistical corrections from MM frequencies. They were performed at the PM6-DH+ level of theory with the full host and guest molecule in the QM system (and also four water molecules in the geometry optimisations) for 10-20 snapshots from molecular dynamics simulations of the complex. Finally, the structure with the lowest free energy was recalculated using the dispersion-corrected density-functional theory method TPSS-D3, for both the structure and the energy. The two FEP approaches gave similar results (PM6-DH+/MM slightly better for OAM), which were among the five submissions with the best performance in the challenge and gave the best results without any fit to data from the SAMPL5 challenge, with mean absolute deviations (MAD) of 2.4-5.2 kJ/mol and a correlation coefficient (R²) of 0.77-0.93. This is the first time QM/MM approaches give binding free energies that are competitive to those obtained with MM for the octa-acid host. The QM/MM-optimised structures gave somewhat worse performance (MAD = 3-8 kJ/mol and R² = 0.1-0.9), but the results were improved compared to previous studies of this system with similar methods.

Bind3P: Optimization of a Water Model Based on Host-Guest Binding Data

We report a water model, Bind3P (Version 0.1), which was obtained by using sensitivity analysis to readjust the Lennard-Jones parameters of the TIP3P model against experimental binding free energies for six host-guest systems, along with pure liquid properties. Tests of Bind3P against >100 experimental binding free energies and enthalpies for host-guest systems distinct from the training set show a consistent drop in the mean signed error, relative to matched calculations with TIP3P. Importantly, Bind3P also yields some improvement in the hydration free energies of small organic molecules and preserves the accuracy of bulk water properties, such as density and the heat of vaporization. The same approach can be applied to more sophisticated water models that can better represent pure water properties. These results lend further support to the concept of integrating host-guest binding data into force field parametrization.

Reversible photoswitching of encapsulated azobenzenes in water

Efficient molecular switching in confined spaces is critical for the successful development of artificial molecular machines. However, molecular switching events often entail large structural changes and therefore require conformational freedom, which is typically limited under confinement conditions. Here, we investigated the behavior of azobenzene-the key building block of light-controlled molecular machines-in a confined environment that is flexible and can adapt its shape to that of the bound guest. To this end, we encapsulated several structurally diverse azobenzenes within the cavity of a flexible, water-soluble coordination cage, and investigated their light-responsive behavior. Using UV/Vis absorption spectroscopy and a combination of NMR methods, we showed that each of the encapsulated azobenzenes exhibited distinct switching properties. An azobenzene forming a 1:1 host-guest inclusion complex could be efficiently photoisomerized in a reversible fashion. In contrast, successful switching in inclusion complexes incorporating two azobenzene guests was dependent on the availability of free cages in the system, and it involved reversible trafficking of azobenzene between the cages. In the absence of extra cages, photoswitching was either suppressed or it involved expulsion of azobenzene from the cage and consequently its precipitation from the solution. This finding was utilized to develop an information storage medium in which messages could be written and erased in a reversible fashion using light.

The Thermodynamics of Anion Complexation to Nonpolar Pockets

The interactions between nonpolar surfaces and polarizable anions lie in a gray area between the hydrophobic and Hofmeister effects. To assess the affinity of these interactions, NMR and ITC were used to probe the thermodynamics of eight anions binding to four different hosts whose pockets each consist primarily of hydrocarbon. Two classes of host were examined: cavitands and cyclodextrins. For all hosts, anion affinity was found to follow the Hofmeister series, with associations ranging from 1.6-5.7 kcal mol-1. Despite the fact that cavitand hosts 1 and 2 possess intrinsic negative electrostatic fields, it was determined that these more enveloping hosts generally bound anions more strongly. The observation that the four hosts each possess specific anion affinities that cannot be readily explained by their structures, points to the importance of counter cations and the solvation of the "empty" hosts, free guests, and host-guest complexes, in defining the affinity.

Tuning the Binding Dynamics of a Guest-Octaacid Capsule through Noncovalent Anchoring

Hydrophobic or hydrophilic substituents have different effects on the binding dynamics of pyrene derivatives with a 2:1 capsule formed from two octaacid cavitands, showing a subtle interplay of different kinetic factors. Anchoring of the methyl group of 1-methylpyrene within one cavitand slowed the association and dissociation dynamics of the 1:1 complex by at least 1000 times when compared to the 1:1 complex for pyrene. This slow down for the transient formation of the 1:1 complex is responsible for the overall increase in stability of the 2:1 complex without affecting the overall capsule dissociation. For 1-pyrenemethanol, its residence time in the 2:1 capsule is shorter compared to that of pyrene despite both guests having similar equilibrium constants for the binding of the second cavitand, suggesting that the hydroxymethyl substituent close to the equatorial region of the capsule can interact with water during the partial opening of the capsule.

Predicting Binding Free Energies: Frontiers and Benchmarks

Binding free energy calculations based on molecular simulations provide predicted affinities for biomolecular complexes. These calculations begin with a detailed description of a system, including its chemical composition and the interactions among its components. Simulations of the system are then used to compute thermodynamic information, such as binding affinities. Because of their promise for guiding molecular design, these calculations have recently begun to see widespread applications in early-stage drug discovery. However, many hurdles remain in making them a robust and reliable tool. In this review, we highlight key challenges of these calculations, discuss some examples of these challenges, and call for the designation of standard community benchmark test systems that will help the research community generate and evaluate progress. In our view, progress will require careful assessment and evaluation of new methods, force fields, and modeling innovations on well-characterized benchmark systems, and we lay out our vision for how this can be achieved.

Binding free energies in the SAMPL5 octa-acid host-guest challenge calculated with DFT-D3 and CCSD(T)

We have tried to calculate the free energy for the binding of six small ligands to two variants of the octa-acid deep cavitand host in the SAMPL5 blind challenge. We employed structures minimised with dispersion-corrected density-functional theory with small basis sets and energies were calculated using large basis sets. Solvation energies were calculated with continuum methods and thermostatistical corrections were obtained from frequencies calculated at the HF-3c level. Care was taken to minimise the effects of the flexibility of the host by keeping the complexes as symmetric and similar as possible. In some calculations, the large net charge of the host was reduced by removing the propionate and benzoate groups. In addition, the effect of a restricted molecular dynamics sampling of structures was tested. Finally, we tried to improve the energies by using the DLPNO-CCSD(T) approach. Unfortunately, results of quite poor quality were obtained, with no correlation to the experimental data, systematically too positive affinities (by ~50 kJ/mol) and a mean absolute error (after removal of the systematic error) of 11-16 kJ/mol. DLPNO-CCSD(T) did not improve the results, so the accuracy is not limited by the energy function. Instead, four likely sources of errors were identified: first, the minimised structures were often incorrect, owing to the omission of explicit solvent. They could be partly improved by performing the minimisations in a continuum solvent with four water molecules around the charged groups of the ligands. Second, some ligands could bind in several different conformations, requiring sampling of reasonable structures. Third, there is an indication the continuum-solvation model has problems to accurately describe the binding of both the negatively and positively charged guest molecules. Fourth, different methods to calculate the thermostatistical corrections gave results that differed by up to 30 kJ/mol and there is an indication that HF-3c overestimates the entropy term. In conclusion, it is a challenge to calculate binding affinities for this octa-acid system with quantum-mechanical methods.

Binding of carboxylate and trimethylammonium salts to octa-acid and TEMOA deep-cavity cavitands

In participation of the fifth statistical assessment of modeling of proteins and ligands (SAMPL5), the strength of association of six guests (3-8) to two hosts (1 and 2) were measured by ¹H NMR and ITC. Each host possessed a unique and well-defined binding pocket, whilst the wide array of amphiphilic guests possessed binding moieties that included: a terminal alkyne, nitro-arene, alkyl halide and cyano-arene groups. Solubilizing head groups for the guests included both positively charged trimethylammonium and negatively charged carboxylate functionality. Measured association constants (K _a ) covered five orders of magnitude, ranging from 56 M^-1 for guest 6 binding with host 2 up to 7.43 × 10⁶ M^-1 for guest 6 binding to host 1.

Preparative scale and convenient synthesis of a water-soluble, deep cavitand

Cavitands are established tools of supramolecular chemistry and molecular recognition, and they are finding increasing application in sensing and sequestration of physiologically relevant molecules in aqueous solution. The synthesis of a water-soluble, deep cavitand is described. The route comprises six (linear) steps from commercially available precursors, and it relies on the fourfold oligomeric cyclization reaction of resorcinol with 2,3-dihydrofuran that leads to the formation of a shallow resorcinarene framework; condensation with aromatic panels, which deepens the hydrophobic binding cavity; construction of rigid urea functionalities on the upper rim; and the introduction of the water-solubilizing methylimidazolium groups on the lower rim. Late intermediates of the synthesis can be used in the preparation of congener cavitands with different properties and applications, and a sample of such a synthetic procedure is included in this protocol. Emphasis is placed on scaled-up reactions and on purification procedures that afford materials in high yield and avoid chromatographic purification. This protocol provides improvements over previously described procedures, and it enables the preparation of sizable amounts of deep cavitands: 7 g of a water-soluble cavitand can be prepared from resorcinol in 13 working days.

From steroids to aqueous supramolecular chemistry: an autobiographical career review

The focus of my group's research is aqueous supramolecular chemistry; we try to understand how chemical entities interact with water and consequently how they interact with each other. This personal history recounts my career experiences that led to his involvement with this fascinating area of science.

ITC and NMR Analysis of the Encapsulation of Fatty Acids within a Water-Soluble Cavitand and its Dimeric Capsule

We report here NMR and ITC studies of the binding of ionizable guests (carboxylate acids) to a deep-cavity cavitand. These studies reveal that the shortest guests favored 1:1 complex formation, but the longer the alkyl chain the more the 2:1 host-guest capsule is favored. For intermediate-sized guests, the equilibrium between these two states is controlled by pH; at low values the capsule containing the carboxylic acid guest is favored, whereas as the pH is raised deprotonation of the guest favors the 1:1 complex. Interestingly, for one host-guest pair the energy required to de-cap the 2:1 capsular complex and form the 1:1 complex is sufficient to shift the pK_a of the guest by ~ 3-4 orders of magnitude (4.1-5.4 kcal mol^-1). The two largest guests examined form stable 2:1 capsules, with in both cases the guest adopting a relatively high energy J-shaped motif. Furthermore, these 2:1 complexes are sufficiently stable that at high pH guest deprotonation occurs without de-capping of the capsule.

Molecular containers assembled through the hydrophobic effect

This review focuses on molecular containers formed by assembly processes driven by the hydrophobic effect, and summarizes the progress made in the field over the last ten years. This small but growing facet of supramolecular chemistry discusses three classes of molecules used by researchers to investigate how self-assembly can be applied to form discrete, mono-dispersed, and structurally well-defined supramolecular entities. The approaches demonstrate the importance of preorganization of arrays of rigid moieties to define a specific form predisposed to bind, fold, or assemble. As the examples demonstrate, studying these systems and their properties is teaching us how to control supramolecular chemistry in water, shedding light on aspects of aqueous solutions chemistry, and illustrating novel applications that harness the unique properties of the hydrophobic effect.

BEDAM binding free energy predictions for the SAMPL4 octa-acid host challenge

The binding energy distribution analysis method (BEDAM) protocol has been employed as part of the SAMPL4 blind challenge to predict the binding free energies of a set of octa-acid host-guest complexes. The resulting predictions were consistently judged as some of the most accurate predictions in this category of the SAMPL4 challenge in terms of quantitative accuracy and statistical correlation relative to the experimental values, which were not known at the time the predictions were made. The work has been conducted as part of a hands-on graduate class laboratory session. Collectively the students, aided by automated setup and analysis tools, performed the bulk of the calculations and the numerical and structural analysis. The success of the experiment confirms the reliability of the BEDAM methodology and it shows that physics-based atomistic binding free energy estimation models, when properly streamlined and automated, can be successfully employed by non-specialists.

Differentiation of small alkane and alkyl halide constitutional isomers via encapsulation

Previously we have demonstrated that host 1 is capable of hydrocarbon gas separation by selective sequestration of butane from a mixture with propane in the headspace above a solution of the host (C. L. D. Gibb, B. C. Gibb, J. Am. Chem. Soc., 2006, 128, 16498-16499). Expanding on the idea of using this host as a means to affect guest discrimination, we report here on NMR studies of the binding of constitutional isomers of two classes of small molecules: hexanes and chloropentanes. Our results indicate that even with these relatively straightforward classes of molecules, guest binding is complicated. Overall, host 1 displays a preference to bind guests with a X-C(R2)-C(R2)-Me (X = Cl or Me) structure, and more generally, a preference for branched guests rather than highly flexible, unbranched derivatives. The complexity of binding of these isomers is magnified when considering molecular differentiation between pairs of guests. In such cases, different guests demonstrated different propensities to self-sort; some self-sort exclusively, while others show very little propensity to do so. However, whereas the pool of guests reveals some general correlations between binding strength and structure, no obvious relationship between structure and degree of self-sorting was observed.

Simulation optimization of spherical non-polar guest recognition by deep-cavity cavitands

Biomimetic deep-cavity cavitand hosts possess unique recognition and encapsulation properties that make them capable of selectively binding a range of non-polar guests within their hydrophobic pocket. Adamantane based derivatives which snuggly fit within the pocket of octa-acid deep cavity cavitands exhibit some of the strongest host binding. Here we explore the roles of guest size and attractiveness on optimizing guest binding to form 1:1 complexes with octa-acid cavitands in water. Specifically we simulate the water-mediated interactions of the cavitand with adamantane and a range of simple Lennard-Jones guests of varying diameter and attractive well-depth. Initial simulations performed with methane indicate hydrated methanes preferentially reside within the host pocket, although these guests frequently trade places with water and other methanes in bulk solution. The interaction strength of hydrophobic guests increases with increasing size from sizes slightly smaller than methane to Lennard-Jones guests comparable in size to adamantane. Over this guest size range the preferential guest binding location migrates from the bottom of the host pocket upwards. For guests larger than adamantane, however, binding becomes less favorable as the minimum in the potential-of-mean force shifts to the cavitand face around the portal. For a fixed guest diameter, the Lennard-Jones well-depth is found to systematically shift the guest-host potential-of-mean force to lower free energies, however, the optimal guest size is found to be insensitive to increasing well-depth. Ultimately our simulations show that adamantane lies within the optimal range of guest sizes with significant attractive interactions to match the most tightly bound Lennard-Jones guests studied.

Converging free energies of binding in cucurbit[7]uril and octa-acid host-guest systems from SAMPL4 using expanded ensemble simulations

Molecular containers such as cucurbit[7]uril (CB7) and the octa-acid (OA) host are ideal simplified model test systems for optimizing and analyzing methods for computing free energies of binding intended for use with biologically relevant protein-ligand complexes. To this end, we have performed initially blind free energy calculations to determine the free energies of binding for ligands of both the CB7 and OA hosts. A subset of the selected guest molecules were those included in the SAMPL4 prediction challenge. Using expanded ensemble simulations in the dimension of coupling host-guest intermolecular interactions, we are able to show that our estimates in most cases can be demonstrated to fully converge and that the errors in our estimates are due almost entirely to the assigned force field parameters and the choice of environmental conditions used to model experiment. We confirm the convergence through the use of alternative simulation methodologies and thermodynamic pathways, analyzing sampled conformations, and directly observing changes of the free energy with respect to simulation time. Our results demonstrate the benefits of enhanced sampling of multiple local free energy minima made possible by the use of expanded ensemble molecular dynamics and may indicate the presence of significant problems with current transferable force fields for organic molecules when used for calculating binding affinities, especially in non-protein chemistries.

Exhaustive docking and solvated interaction energy scoring: lessons learned from the SAMPL4 challenge

We continued prospective assessments of the Wilma-solvated interaction energy (SIE) platform for pose prediction, binding affinity prediction, and virtual screening on the challenging SAMPL4 data sets including the HIV-integrase inhibitor and two host-guest systems. New features of the docking algorithm and scoring function are tested here prospectively for the first time. Wilma-SIE provides good correlations with actual binding affinities over a wide range of binding affinities that includes strong binders as in the case of SAMPL4 host-guest systems. Absolute binding affinities are also reproduced with appropriate training of the scoring function on available data sets or from comparative estimation of the change in target's vibrational entropy. Even when binding modes are known, SIE predictions lack correlation with experimental affinities within dynamic ranges below 2 kcal/mol as in the case of HIV-integrase ligands, but they correctly signaled the narrowness of the dynamic range. Using a common protein structure for all ligands can reduce the noise, while incorporating a more sophisticated solvation treatment improves absolute predictions. The HIV-integrase virtual screening data set consists of promiscuous weak binders with relatively high flexibility and thus it falls outside of the applicability domain of the Wilma-SIE docking platform. Despite these difficulties, unbiased docking around three known binding sites of the enzyme resulted in over a third of ligands being docked within 2 Å from their actual poses and over half of the ligands docked in the correct site, leading to better-than-random virtual screening results.

Binding of cyclic carboxylates to octa-acid deep-cavity cavitand

As part of the fourth statistical assessment of modeling of proteins and ligands (sampl.eyesopen.com) prediction challenge, the strength of association of nine guests (1-9) binding to octa-acid host was determined by a combination of (1)H NMR and isothermal titration calorimetry. Association constants in sodium tetraborate buffered (pH 9.2) aqueous solution ranged from 5.39 × 10(2) M(-1) in the case of benzoate 1, up to 3.82 × 10(5) M(-1) for trans-4-methylcyclohexanoate 7. Overall, the free energy difference between the free energies of complexation of these weakest and strongest binding guests was ΔΔG° = 3.88 kcal mol(-1). Based on a multitude of previous studies, the anticipated order of strength of binding was close to that which was actually obtained. However, the binding of guest 3 (4-ethylbenzoate) was considerably stronger than initially estimated.

Guest packing motifs within a supramolecular nanocapsule and a covalent analogue

Two hosts that utilize the hydrophobic effect to assemble and/or encapsulate guest molecules were studied. The hosts, octa-acid (OA) and hexalene diamine-linked octa-acid (HOA), were shown to complex a broad range of n-alkanes up to n-hexacosane (C26H54). A combination of (1)H NMR, NMR diffusion, COSY, and NOESY experiments revealed four different guest packing motifs, depending on the size of the guest and the nature of the host. As a function of guest size, smooth transitions from one motif to the next were observed and allowed qualification of their relative stabilities. Furthermore, although the two hosts engender ostensibly identical encapsulation environments, their different assembly properties lead to quite distinct packing-motif profiles, i.e., how the motifs change as a function of guest size.

Anion binding to hydrophobic concavity is central to the salting-in effects of Hofmeister chaotropes

For over 120 years it has been appreciated that certain salts (kosmotropes) cause the precipitation of proteins, while others (chaotropes) increase their solubility. The cause of this "Hofmeister effect" is still unclear, especially with the original concept that kosmotropic anions "make" water structure and chaotropes "break" it being countered by recent studies suggesting otherwise. Here, we present the first direct evidence that chaotropic anions have an affinity for hydrophobic concavity and that it is competition between a convex hydrophobe and the anion for a binding site that leads to the apparent weakening of the hydrophobic effect by chaotropes. In combination, these results suggest that chaotropes primarily induce protein solubilization by direct binding to concavity in the molten globule state of a protein.

Nature of supramolecular complexes controlled by the structure of the guest molecules: formation of octa acid based capsuleplex and cavitandplex

Factors that govern inclusion of organic molecules within octa acid (OA), a synthetic deep cavity cavitand, have been delineated by examining the complexation behavior of a number of organic molecules with varying dimensions and functionalities with OA. The formation of two types of complexes has been noted: the one which we call cavitandplex is a partially open complex in which a part of the guest molecule remains exposed to water, and the other termed capsuleplex is formed through assembly of two OA molecules. In capsuleplex, the guest is protected from water. Generally, guest molecules that possess ionic head groups form cavitandplex, and all others form capsuleplex. Capsuleplex may contain one or two guest molecules within the capsule. Small organic molecules (<10 A in length) may form both 2:1 and 2:2 capsuleplex, while longer ones (>12 A) preferentially form 2:1 capsuleplex. Extensive 1H NMR experiments have been carried out to characterize host-guest complexes. In the absence of the guest, OA tends to aggregate in water. The extent of aggregation depends on the concentration of OA and the presence of salts in solution. We expect the information obtained from this study to be of great value in predicting the nature of complexes with a given guest and facilitating appropriate guest chosen by researchers.

Dendronized supramolecular nanocapsules: pH independent, water-soluble, deep-cavity cavitands assemble via the hydrophobic effect

At neutral pH, dendronized deep-cavity cavitands were shown to form supramolecular nanocapsules via assembly around a range of guest molecules.

Water inside a hydrophobic cavitand molecule

The structure and dynamics of water inside a water-soluble, bowl-shaped cavitand molecule with a hydrophobic interior are studied using molecular dynamics computer simulations. The simulations find that the number of inside water molecules is about 4.5, but it fluctuates from being completely empty to full on a time scale of tens of nanoseconds. The transition from empty to full is energetically favorable and entropically unfavorable. The water molecules inside have fewer hydrogen bonds than the bulk and in general weaker interactions; the lower energy results from the nearest-neighbor interactions with the cavitand atoms and the water molecules at the entrance of the cavitand, interactions that are lost upon dewetting. An analysis of translational and rotational motion suggests that the lower entropy of the inside water molecules is due to decreased translational entropy, which outweighs an increased orientational entropy. The cavitand molecule acts as a host binding hydrophobic guests, and dewetting can be induced by the presence of a hydrophobic guest molecule about 3 A above the entrance. At this position, the guest displaces the water molecules which stabilize the inside water molecules and the empty cavitand becomes more stable than the full.