skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


Title: Calculating Binding Free Energies in Model Host–Guest Systems with Unrestrained Advanced Sampling
Host–guest interactions are important to the design of pharmaceuticals and, more broadly, to soft materials as they can enable targeted, strong, and specific interactions between molecules. The binding process between the host and guest may be classified as a “rare event” when viewing the system at atomic scales, such as those explored in molecular dynamics simulations. To obtain equilibrium binding conformations and dissociation constants from these simulations, it is essential to resolve these rare events. Advanced sampling methods such as the adaptive biasing force (ABF) promote the occurrence of less probable configurations in a system, therefore facilitating the sampling of essential collective variables that characterize the host–guest interactions. Here, we present the application of ABF to a rod–cavitand coarse-grained model of host–guest systems to acquire the potential of mean force. We show that the employment of ABF enables the computation of the configurational and thermodynamic properties of bound and unbound states, including the free energy landscape. Moreover, we identify important dynamic bottlenecks that limit sampling and discuss how these may be addressed in more general systems.  more » « less
Award ID(s):
1751988
PAR ID:
10508419
Author(s) / Creator(s):
; ;
Publisher / Repository:
American Chemical Society
Date Published:
Journal Name:
Journal of Chemical Theory and Computation
Volume:
20
Issue:
9
ISSN:
1549-9618
Page Range / eLocation ID:
3927 to 3934
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; (, Journal of Computational Chemistry)
    Abstract The prediction of (un)binding rates and free energies is of great significance to the drug design process. Although many enhanced sampling algorithms and approaches have been developed, there is not yet a reliable workflow to predict these quantities. Previously we have shown that free energies and transition rates can be calculated by directly simulating the binding and unbinding processes with our variant of the WE algorithm “Resampling of Ensembles by Variation Optimization”, or “REVO”. Here, we calculate binding free energies retrospectively for three SAMPL6 host‐guest systems and prospectively for a SAMPL9 system to test a modification of REVO that restricts its cloning behavior in quasi‐unbound states. Specifically, trajectories cannot clone if they meet a physical requirement that represents a high likelihood of unbinding, which in the case of this work is a center‐of‐mass to center‐of‐mass distance. The overall effect of this change was difficult to predict, as it results in fewer unbinding events each of which with a much higher statistical weight. For all four systems tested, this new strategy produced either more accurate unbinding free energies or more consistent results between simulations than the standard REVO algorithm. This approach is highly flexible, and any feature of interest for a system can be used to determine cloning eligibility. These findings thus constitute an important improvement in the calculation of transition rates and binding free energies with the weighted ensemble method. 
    more » « less
  2. ; ; ; ; ; ; ; (, Chemical Science)
    Science still does not have the ability to accurately predict the affinity that ligands have for proteins. In an attempt to address this, the Statistical Assessment of Modeling of Proteins and Ligands (SAMPL) series of blind predictive challenges is a community-wide exercise aimed at advancing computational techniques as standard predictive tools in rational drug design. In each cycle, a range of biologically relevant systems of different levels of complexity are selected to test the latest modeling methods. As part of this on-going exercise, and as a step towards understanding the important factors in context dependent guest binding, we challenged the computational community to determine the affinity of a series of negatively and positively charged guests to two constitutionally isomeric cavitand hosts: octa-acid 1 , and exo -octa acid 2 . Our affinity determinations, combined with molecular dynamics simulations, reveal asymmetries in affinities between host–guest pairs that cannot alone be explained by simple coulombic interactions, but also point to the importance of host–water interactions. Our work reveals the key facets of molecular recognition in water, emphasizes where improvements need to be made in modelling, and shed light on the complex problem of ligand-protein binding in the aqueous realm. 
    more » « less
  3. ; (, Frontiers in Molecular Biosciences)
    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
  4. ; ; ; ; ; ; (, Physical Chemistry Chemical Physics)
    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
  5. ; ; ; ; ; ; ; (, Angewandte Chemie International Edition)
    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
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email