An official website of the United States government
We apply the Alchemical Transfer Method (ATM) and a bespoke fixed partial charge force field to the SAMPL9 bCD host-guest binding free energy prediction challenge that comprises a combination of complexes formed between five phenothiazine guests and two cyclodextrin hosts. Multiple chemical forms, competing binding poses, and computational modeling challenges pose significant obstacles to obtaining reliable computational predictions for these systems. The phenothiazine guests exist in solution as racemic mixtures of enantiomers related by nitrogen inversions that bind the hosts in various binding poses, each requiring an individual free energy analysis. Due to the large size of the guests and the conformational reorganization of the hosts, which prevent a direct absolute binding free energy route, binding free energies are obtained by a series of absolute and relative binding alchemical steps for each chemical species in each binding pose. Metadynamics-accelerated conformational sampling was found to be necessary to address the poor convergence of some numerical estimates affected by conformational trapping. Despite these challenges, our blinded predictions quantitatively reproduced the experimental affinities for the beta-cyclodextrin host and, to a lesser extent, those with a methylated derivative. The work illustrates the challenges of obtaining reliable free energy data in in-silico drug design for even seemingly simple systems and introduces some of the technologies available to tackle them.
more »
« less
Cucurbiturils mimicked by low polarizability solvents with pre-formed cavities: an empirical model to predict hydrocarbon selectivity
Relative binding affinities of a series of nine rigid hydrocarbons towards the cavity formed by a portion of the inner wall of cucurbit[8]uril (CB[8]) and a positive auxiliary guest were determined by competitive 19 F NMR titrations in deuterium oxide. The corresponding free binding energies were corrected by the hydrocarbon computed solvation energies to obtain their free energies of transfer from the gas phase to the CB[8]/auxiliary guest cavity. These energies correlate linearly with the hydrocarbon static polarizabilities, thereby suggesting that the selectivity is driven, perhaps exclusively, by dispersive interactions between the hydrocarbons and the tailor-made cavity, regardless of the degree of unsaturation of the guests. The free energies of transfer also correlate linearly with the energy released upon introduction of the hydrocarbon into a pre-formed cavity extruded from a solvent (benzene) selected to mimic the polarity and polarizability of the CB[8]/auxiliary probe cavity – and this, with a unity slope. Among other features, this empirical model also accurately predicts the relative binding affinities of various rigid hydrocarbons to CB[6] and CB[7], as well as noble gases to CB[5], when the macrocycles are mimicked with pre-formed cavities in perfluorohexane or perfluorohexane/benzene mixtures, both being notoriously non-polar and non-polarizable environments.
more »
« less
- PAR ID:
- 10355232
- Date Published:
- Journal Name:
- Chemical Science
- Volume:
- 13
- Issue:
- 15
- ISSN:
- 2041-6520
- Page Range / eLocation ID:
- 4388 to 4396
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Improper reaction coordinates can pose significant problems for path-based binding free energy calculations. Particularly, omission of long timescale motions can lead to over-estimation of the energetic barriers between the bound and unbound states. Many methods exist to construct the optimal reaction coordinate using a pre-defined basis set of features. Although simulations are typically conducted in explicit solvent, the solvent atoms are often excluded by these feature sets—resulting in little being known about their role in reaction coordinates, and ultimately, their role in determining (un)binding rates and free energies. In this work, analysis is done on an extensive set of host-guest unbinding trajectories, working to characterize differences between high and low probability unbinding trajectories with a focus on solvent-based features, including host-ion interactions, guest-ion interactions and location-dependent ion densities. We find that differences in ion densities as well as guest-ion interactions strongly correlate with differences in the probabilities of reactive paths that are used to determine free energies of (un)binding and play a significant role in the unbinding process.more » « less
-
Abstract Cucurbit[n]urils (CB[n]s) are cyclic macrocycles with rich host‐guest chemistry. In many cases, guest binding in CB[n]s results in host structural deformations. Unfortunately, measuring such deformations remains a major challenge, with only a handful of manual estimations reported in the literature. To address this challenge, we have developed the public program ElliptiCB[n], which is available on GitHub, that provides a robust and automated method for measuring the elliptical deformations in CB[n] hosts. We outline the development and validation of this approach, apply ElliptiCB[n] to measure the ellipticity of the 1113 available CB[n] structures from the Cambridge Structural Database (CSD), and directly investigate the structural deformations of CB[5], CB[6], CB[7], CB[8], and CB[10] hosts. We also report the general landscape of accessible CB[n] elliptical deformations and compare ellipticity distributions across CB[n] hosts and host‐guest complexes. We found that in almost all cases guest binding significantly impacts the distribution of host ellipticity distributions and that these distributions are dissimilar across host‐guest complexes of differently sized CB[n]s. We anticipate that this work will provide a useful approach for understanding of the flexibility of CB[n] hosts and will also enable future measurement and standardization of ellipticity measurements of CB[n]s.more » « less
-
Self-assembled Fe4 L 6 cage complexes with variable internal functions can be synthesized from a 2,7-dibromocarbazole ligand scaffold, which orients six functional groups to the cage interior. Both ethylthiomethylether and ethyldimethylamino groups can be incorporated. The cages show strong ligand-centered fluorescence emission and a broad range of guest binding properties. Coencapsulation of neutral organic guests is favored in the larger, unfunctionalized cage cavity, whereas the thioether cage has a more sterically hindered cavity that favors 1 : 1 guest binding. Binding affinities up to 10 6 M −1 in CH3 CN are seen. The dimethylamino cage is more complex, as the internal amines display partial protonation and can be deprotonated by amine bases. This amine cage displays affinity for a broad range of neutral organic substrates, with affinities and stoichiometries comparable to that of the similarly sized thioether cage. These species show that simple variations in ligand backbone allow variations in the number and type of functions that can be displayed towards the cavity of self-assembled hosts, which will have applications in biomimetic sensing, catalysis and molecular recognition.more » « less
-
Abstract Cooperativity plays a critical role in self‐assembly and molecular recognition. A rigid aromatic oligoamide macrocycle with a cyclodirectional backbone binds with DABCO‐based cationic guests in a 2 : 1 ratio in high affinities (Ktotal≈1013 M−2) in the highly polar DMF. The host–guest binding also exhibits exceptionally strong positive cooperativity quantified by interaction factors α that are among the largest for synthetic host–guest systems. The unusually strong positive cooperativity, revealed by isothermal titration calorimetry (ITC) and fully corroborated by mass spectrometry, NMR and computational studies, is driven by guest‐induced stacking of the macrocycles and stabilization from the alkyl end chains of the guests, interactions that appear upon binding the second macrocycle. With its tight binding driven by extraordinary positive cooperativity, this host–guest system provides a tunable platform for studying molecular interactions and for constructing stable supramolecular assemblies.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1039/D1SC06728A
