Thermodynamically Optimized Machine-learned Reaction Coordinates for Hydrophobic Ligand Dissociation
arXiv:2310.03819 · doi:10.1021/acs.jpcb.3c08304
Abstract
Ligand unbinding is mediated by the free energy change, which has intertwined contributions from both energy and entropy. It is important but not easy to quantify their individual contributions. We model hydrophobic ligand unbinding for two systems, a methane particle and a C60 fullerene, both unbinding from hydrophobic pockets in all-atom water. By using a modified deep learning framework, we learn a thermodynamically optimized reaction coordinate to describe hydrophobic ligand dissociation for both systems. Interpretation of these reaction coordinates reveals the roles of entropic and enthalpic forces as ligand and pocket sizes change. Irrespective of the contrasting roles of energy and entropy, we also find that for both the systems the transition from the bound to unbound states is driven primarily by solvation of the pocket and ligand, independent of ligand size. Our framework thus gives useful thermodynamic insight into hydrophobic ligand dissociation problems that are otherwise difficult to glean.
arXiv comment27 pages; 5 figures
References in corpus (18)
- PLUMED 2: New feathers for an old bird
- Well-Tempered Metadynamics: A Smoothly Converging and Tunable Free-Energy Method
- Identification of slow molecular order parameters for Markov model construction
- Sampling rare switching events in biochemical networks
- Time-lagged autoencoders: Deep learning of slow collective variables for molecular kinetics
- Caliber based spectral gap optimization of order parameters (SGOOP) for sampling complex molecular systems
- Data-driven model reduction and transfer operator approximation
- Molecular enhanced sampling with autoencoders: On-the-fly collective variable discovery and accelerated free energy landscape exploration
- Variational Encoding of Complex Dynamics
- A variational conformational dynamics approach to the selection of collective variables in metadynamics
- State Predictive Information Bottleneck
- The role of water and steric constraints in the kinetics of cavity-ligand unbinding
- Driving and characterizing nucleation of urea and glycine polymorphs in water
- How wet should be the reaction coordinate for ligand unbinding?
- Assessing models of force-dependent unbinding rates via infrequent metadynamics
- Is the Local Ion Density Sufficient to Drive NaCl Nucleation from the Melt and Aqueous Solution?
- On the Role of Solvent in Hydrophobic Cavity-ligand Recognition Kinetics
- Thermodynamics-inspired Explanations of Artificial Intelligence