On the Nature of Guest Complexation in Water: Triggered Wetting-Water-Mediated Binding

Potential distribution theory of alchemical transfer

We present an analytical description of the Alchemical Transfer Method (ATM) for molecular binding using the Potential Distribution Theory (PDT) formalism. ATM models the binding free energy by mapping the bound and unbound states of the complex by translating the ligand coordinates. PDT relates the free energy and the probability densities of the perturbation energy along the alchemical path to the probability density at the initial state, which is the unbound state of the complex in the case of a binding process. Hence, the ATM probability density of the transfer energy at the unbound state is first related by a convolution operation of the probability densities for coupling the ligand to the solvent and coupling it to the solvated receptor-for which analytical descriptions are available-with parameters obtained from maximum likelihood analysis of data from double-decoupling alchemical calculations. PDT is then used to extend this analytical description along the alchemical transfer pathway. We tested the theory on the alchemical binding of five guests to the tetramethyl octa-acid host from the SAMPL8 benchmark set. In each case, the probability densities of the perturbation energy for transfer along the alchemical transfer pathway obtained from numerical calculations match those predicted from the theory and double-decoupling simulations. The work provides a solid theoretical foundation for alchemical transfer, offers physical insights on the form of the probability densities observed in alchemical transfer calculations, and confirms the conceptual and numerical equivalence between the alchemical transfer and double-decoupling processes.

Under the Influence of Water: Molecular Recognition of Organic Hydrophilic Molecules in Water with a Prismarene Host Driven by Hydration Effects

A water-soluble prism[5]arene host can form endo-cavity complexes with hydrophilic organic substances in water by displacing frustrated water molecules from its deep cavity. Water molecules structured at both rims of the prismarene host can mediate hydrogen bonding interactions with the guest. Water-mediated hydrogen bonding interactions were invoked here to elucidate the elevated binding affinities and selectivity of the prismarene host toward hydrophilic organic guests. We show that water at the interface of a host-guest complex can act as an extension of the host structure, facilitating the accommodation of neutral guests within the binding site. This study highlights the crucial role of water in facilitating supramolecular interactions between a deep-cavity prismarene host and organic hydrophilic guests in aqueous medium.

An Affordable Topography-Based Protocol for Assigning a Residue's Character on a Hydropathy (PARCH) Scale

The hydropathy of proteins or quantitative assessment of protein-water interactions has been a topic of interest for decades. Most hydropathy scales use a residue-based or atom-based approach to assign fixed numerical values to the 20 amino acids and categorize them as hydrophilic, hydroneutral, or hydrophobic. These scales overlook the protein's nanoscale topography, such as bumps, crevices, cavities, clefts, pockets, and channels, in calculating the hydropathy of the residues. Some recent studies have included protein topography in determining hydrophobic patches on protein surfaces, but these methods do not provide a hydropathy scale. To overcome the limitations in the existing methods, we have developed a Protocol for Assigning a Residue's Character on the Hydropathy (PARCH) scale that adopts a holistic approach to assigning the hydropathy of a residue. The parch scale evaluates the collective response of the water molecules in the protein's first hydration shell to increasing temperatures. We performed the parch analysis of a set of well-studied proteins that include the following─enzymes, immune proteins, and integral membrane proteins, as well as fungal and virus capsid proteins. Since the parch scale evaluates every residue based on its location, a residue may have very different parch values inside a crevice versus a surface bump. Thus, a residue can have a range of parch values (or hydropathies) dictated by the local geometry. The parch scale calculations are computationally inexpensive and can compare hydropathies of different proteins. The parch analysis can affordably and reliably aid in designing nanostructured surfaces, identifying hydrophilic and hydrophobic patches, and drug discovery.

Accurate Host-Guest Binding Free Energies Using the AMOEBA Polarizable Force Field

A grand challenge of computational biophysics is accurate prediction of interactions between molecules. Molecular dynamics (MD) simulations have recently gained much interest as a tool to directly compute rigorous intermolecular binding affinities. The choice of a fixed point-charge or polarizable multipole force field used in MD is a topic of ongoing discussion. To compare alternative methods, we participated in the SAMPL7 and SAMPL8 Gibb octaacid host-guest challenges to assess the Atomic Multipole Optimized Energetics for Biomolecular Applications (AMOEBA) polarizable multipole force field. Advantages of AMOEBA over fixed charge models include improved representation of molecular electrostatic potentials and better description of water occupying the unligated host cavity. Prospective predictions for 26 host-guest systems exhibit a mean unsigned error vs experiment of 0.848 kcal/mol across all absolute binding free energies, demonstrating excellent agreement between computational and experimental results. In addition, we explore two topics related to the inclusion of ions in MD simulations: use of a neutral co-alchemical protocol and the effect of salt concentration on binding affinity. Use of the co-alchemical method minimally affects computed energies, but salt concentration significantly perturbs our binding results. Higher salt concentration strengthens binding through classical charge screening. In particular, added Na+ ions screen negatively charged carboxylate groups near the binding cavity, thereby diminishing repulsive coulomb interactions with negatively charged guests. Overall, the AMOEBA results demonstrate the accuracy available through a force field providing a detailed energetic description of the four octaacid hosts and 13 charged organic guests. Use of the AMOEBA polarizable atomic multipole force field in conjunction with an alchemical free energy protocol can achieve chemical accuracy in application to realistic molecular systems.