An official website of the United States government
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
Local Ion Densities can Influence Transition Paths of Molecular Binding
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
- Award ID(s):
- 1761320
- PAR ID:
- 10469090
- Publisher / Repository:
- Frontiers
- Date Published:
- Journal Name:
- Frontiers in Molecular Biosciences
- Volume:
- 9
- ISSN:
- 2296-889X
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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
-
Dibenzo-18-crown-6 (DB18C6) is a single-crown ether that can act as a host for a guest ion. In an effort to illuminate the relationships among structure, dynamics, and thermodynamics of ligand binding in a simple model for understanding the affinity and specificity of ligand interactions, nuclear magnetic resonance (NMR) experiments and density functional theory (DFT) were used to study the interaction of DB18C6 with ammonium ion. 1H-NMR was used to follow the titration of DB18C6 with ammonium chloride in deuterated methanol, a solvent chosen for its amphipathic character. Ammonium ion binds strongly to DB18C6 with a dissociation equilibrium constant at least as low as ~ 10 - 6 M. DFT calculations were used to identify optimized conformations of bound and free DB18C6 and to estimate its binding energy with ammonium ion in implicit solvent. An approach is described that accounts for geometry relaxation in addition to solvation correction and basis set superposition error; to our knowledge, this is the first such report that includes the energy difference from optimizing species geometry. The lowest-energy conformer of free DB18C6 in implicit methanol acquires an open, W-shaped structure that is also the lowest-energy conformer found for the DB18C6-ammonium ion complex. These results form a foundation for further studies of this system by molecular dynamics simulations.more » « less
-
We propose a free energy calculation method for receptor–ligand binding, which have multiple binding poses that avoids exhaustive enumeration of the poses. For systems with multiple binding poses, the standard procedure is to enumerate orientations of the binding poses, restrain the ligand to each orientation, and then, calculate the binding free energies for each binding pose. In this study, we modify a part of the thermodynamic cycle in order to sample a broader conformational space of the ligand in the binding site. This modification leads to more accurate free energy calculation without performing separate free energy simulations for each binding pose. We applied our modification to simple model host–guest systems as a test, which have only two binding poses, by using a single decoupling method (SDM) in implicit solvent. The results showed that the binding free energies obtained from our method without knowing the two binding poses were in good agreement with the benchmark results obtained by explicit enumeration of the binding poses. Our method is applicable to other alchemical binding free energy calculation methods such as the double decoupling method (DDM) in explicit solvent. We performed a calculation for a protein–ligand system with explicit solvent using our modified thermodynamic path. The results of the free energy simulation along our modified path were in good agreement with the results of conventional DDM, which requires a separate binding free energy calculation for each of the binding poses of the example of phenol binding to T4 lysozyme in explicit solvent. © 2019 Wiley Periodicals, Inc.more » « less
-
Molecular tweezers using aromatic interactions for binding normally work best in polar instead of nonpolar solvents due to the strong solvophobic effect in the binding. Inspired by biological receptors that utilize “delocalized binding interactions” remote from the binding interface to strengthen guest-binding, we constructed molecular tweezers that have a reversed solvent effect. As the direct aromatic binding interactions were weakened by nonpolar solvent, guest-triggered intrahost interactions between two strategically placed carboxylic acids became stronger and contributed to the binding.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.3389/fmolb.2022.858316
