Principles for Tuning Hydrophobic Ligand–Receptor Binding Kinetics

Databases of ligand-binding pockets and protein-ligand interactions

Many research groups and institutions have created a variety of databases curating experimental and predicted data related to protein-ligand binding. The landscape of available databases is dynamic, with new databases emerging and established databases becoming defunct. Here, we review the current state of databases that contain binding pockets and protein-ligand binding interactions. We have compiled a list of such databases, fifty-three of which are currently available for use. We discuss variation in how binding pockets are defined and summarize pocket-finding methods. We organize the fifty-three databases into subgroups based on goals and contents, and describe standard use cases. We also illustrate that pockets within the same protein are characterized differently across different databases. Finally, we assess critical issues of sustainability, accessibility and redundancy.

Thermodynamically Optimized Machine-Learned Reaction Coordinates for Hydrophobic Ligand Dissociation

Ligand unbinding is mediated by its free energy change, which has intertwined contributions from both energy and entropy. It is important, but not easy, to quantify their individual contributions to the free energy profile. We model hydrophobic ligand unbinding for two systems, a methane particle and a C60 fullerene, both unbinding from hydrophobic pockets in all-atom water. Using a modified deep learning framework, we learn a thermodynamically optimized reaction coordinate to describe the hydrophobic ligand dissociation for both systems. Interpretation of these reaction coordinates reveals the roles of entropic and enthalpic forces as the ligand and pocket sizes change. In both cases, we observe that the free-energy barrier to unbinding is dominated by entropy considerations. Furthermore, the process of methane unbinding is driven by methane solvation, while fullerene unbinding is driven first by pocket wetting and then fullerene wetting. For both solutes, the direct importance of the distance from the binding pocket to the learned reaction coordinate is present, but low. Our framework and subsequent feature important analysis thus give useful thermodynamic insight into hydrophobic ligand dissociation problems that are otherwise difficult to glean.

Tuning Local Hydration Enables a Deeper Understanding of Protein-Ligand Binding: The PP1-Src Kinase Case

Water plays a key role in biomolecular recognition and binding. Despite the development of several computational and experimental approaches, it is still challenging to comprehensively characterize water-mediated effects on the binding process. Here, we investigate how water affects the binding of Src kinase to one of its inhibitors, PP1. Src kinase is a target for treating several diseases, including cancer. We use biased molecular dynamics simulations, where the hydration of predetermined regions is tuned at will. This computational technique efficiently accelerates the SRC-PP1 binding simulation and allows us to identify several key and yet unexplored aspects of the solvent's role. This study provides a further perspective on the binding phenomenon, which may advance the current drug design approaches for the development of new kinase inhibitors.

Thermodynamics and Kinetics of Drug-Target Binding by Molecular Simulation

Computational studies play an increasingly important role in chemistry and biophysics, mainly thanks to improvements in hardware and algorithms. In drug discovery and development, computational studies can reduce the costs and risks of bringing a new medicine to market. Computational simulations are mainly used to optimize promising new compounds by estimating their binding affinity to proteins. This is challenging due to the complexity of the simulated system. To assess the present and future value of simulation for drug discovery, we review key applications of advanced methods for sampling complex free-energy landscapes at near nonergodicity conditions and for estimating the rate coefficients of very slow processes of pharmacological interest. We outline the statistical mechanics and computational background behind this research, including methods such as steered molecular dynamics and metadynamics. We review recent applications to pharmacology and drug discovery and discuss possible guidelines for the practitioner. Recent trends in machine learning are also briefly discussed. Thanks to the rapid development of methods for characterizing and quantifying rare events, simulation's role in drug discovery is likely to expand, making it a valuable complement to experimental and clinical approaches.

Transient States and Barriers from Molecular Simulations and the Milestoning Theory: Kinetics in Ligand-Protein Recognition and Compound Design

This study presents a novel computational approach to study molecular recognition and binding kinetics for drug-like compounds dissociating from a flexible protein system. The intermediates and their free energy profile during ligand association and dissociation processes control ligand-protein binding kinetics and bring a more complete picture of ligand-protein binding. The method applied the milestoning theory to extract kinetics and thermodynamics information from running short classical molecular dynamics (MD) simulations for frames from a given dissociation path. High-dimensional ligand-protein motions (3N-6 degrees of freedom) during ligand dissociation were reduced by use of principal component modes for assigning more than 100 milestones, and classical MD runs were allowed to travel multiple milestones to efficiently obtain ensemble distribution of initial structures for MD simulations and estimate the transition time and rate during ligand traveling between milestones. We used five pyrazolourea ligands and cyclin-dependent kinase 8 with cyclin C (CDK8/CycC) as our model system as well as metadynamics and a pathway search method to sample dissociation pathways. With our strategy, we constructed the free energy profile for highly mobile biomolecular systems. The computed binding free energy and residence time correctly ranked the pyrazolourea ligand series, in agreement with experimental data. Guided by a barrier of a ligand passing an αC helix and activation loop, we introduced one hydroxyl group to parent compounds to design our ligands with increased residence time and validated our prediction by experiments. This work provides a novel and robust approach to investigate dissociation kinetics of large and flexible systems for understanding unbinding mechanisms and designing new small-molecule drugs with desired binding kinetics.

Variational implicit-solvent predictions of the dry-wet transition pathways for ligand-receptor binding and unbinding kinetics

Ligand-receptor binding and unbinding are fundamental biomolecular processes and particularly essential to drug efficacy. Environmental water fluctuations, however, impact the corresponding thermodynamics and kinetics and thereby challenge theoretical descriptions. Here, we devise a holistic, implicit-solvent, multimethod approach to predict the (un)binding kinetics for a generic ligand-pocket model. We use the variational implicit-solvent model (VISM) to calculate the solute-solvent interfacial structures and the corresponding free energies, and combine the VISM with the string method to obtain the minimum energy paths and transition states between the various metastable ("dry" and "wet") hydration states. The resulting dry-wet transition rates are then used in a spatially dependent multistate continuous-time Markov chain Brownian dynamics simulation and the related Fokker-Planck equation calculations of the ligand stochastic motion, providing the mean first-passage times for binding and unbinding. We find the hydration transitions to significantly slow down the binding process, in semiquantitative agreement with existing explicit-water simulations, but significantly accelerate the unbinding process. Moreover, our methods allow the characterization of nonequilibrium hydration states of pocket and ligand during the ligand movement, for which we find substantial memory and hysteresis effects for binding vs. unbinding. Our study thus provides a significant step forward toward efficient, physics-based interpretation and predictions of the complex kinetics in realistic ligand-receptor systems.

Drug-Target Association Kinetics in Drug Discovery

The important role of ligand-receptor binding kinetics in drug design and discovery is increasingly recognized by the drug research community. Over the past decade, accumulating evidence has shown that optimizing the ligand's dissociation rate constant can lead to desirable duration of in vivo target occupancy and, hence, improved pharmacodynamic properties. However, the association rate constant as a pharmacological principle remains less investigated, whereas it can play an equally important role in the selection of drug candidates. This review provides a compilation and discussion of otherwise scarce and dispersed information on this topic, bringing to light the importance of drug-target association in kinetics-directed drug design and discovery.

Transfer Free Energies and Partitioning of Small Molecules in Collapsed PNIPAM Polymers

A central quantity in the design of functional hydrogels used as nanocarrier systems, for instance, for drug delivery or adaptive nanocatalysis, is the partition ratio, which quantifies the uptake of a molecular substance by the polymer matrix. By employing all-atom molecular dynamics simulations, we study the solvation and partitioning (with respect to bulk water) of small subnanometer-sized solutes in a dense matrix of collapsed poly( N-isopropylacrylamide) (PNIPAM) polymers above the lower critical solution temperature in aqueous solution. We examine the roles of the solute's polarity and its size on the solubility properties in the thermoresponsive polymer. We show that the transfer free energies of nonpolar solutes from bulk water into the polymer are favorable and scale in a good approximation with the solute's surface area. Even for small solute size variation, partitioning can vary over orders of magnitude. A polar nature of the solute, on the other hand, generally opposes the transfer, at least for alkyl solutes. Finally, we find a strong correlation between the transfer free energies in the gel and the adsorption free energies on a single extended polymer chain, which enables us to relate the partition ratios in the swollen and collapsed state of a PNIPAM gel.

Role of Molecular Interactions and Protein Rearrangement in the Dissociation Kinetics of p38α MAP Kinase Type-I/II/III Inhibitors

Understanding the governing factors of fast or slow inhibitor binding/unbinding assists in developing drugs with preferred kinetic properties. For inhibitors with the same binding affinity targeting different binding sites of the same protein, the kinetic behavior can profoundly differ. In this study, we investigated unbinding kinetics and mechanisms of fast (type-I) and slow (type-II/III) binders of p38α mitogen-activated protein kinase, where the crystal structures showed that type-I and type-II/III inhibitors bind to pockets with different conformations of the Asp-Phe-Gly (DFG) motif. The work used methods that combine conventional molecular dynamics (MD), accelerated molecular dynamics (AMD) simulations, and the newly developed pathway search guided by internal motions (PSIM) method to find dissociation pathways. The study focuses on revealing key interactions and molecular rearrangements that hinder ligand dissociation by using umbrella sampling and post-MD processing to examine changes in free energy during ligand unbinding. As anticipated, the initial dissociation steps all require breaking interactions that appeared in crystal structures of the bound complexes. Interestingly, for type-I inhibitors such as SB2, p38α keeps barrier-free conformational fluctuation in the ligand-bound complex and during ligand dissociation. In contrast, with a type-II/III inhibitor such as BIRB796, with the rearrangements of p38α in its bound state, ligand unbinding features energetically unfavorable protein-ligand concerted movement. Our results also show that the type-II/III inhibitors preferred dissociation pathways through the allosteric channel, which is consistent with an existing publication. The study suggests that the level of required protein rearrangement is one major determining factor of drug binding kinetics in p38α systems, providing useful information for development of inhibitors.

Kinetics of human cannabinoid 1 (CB1) receptor antagonists: Structure-kinetics relationships (SKR) and implications for insurmountable antagonism

While equilibrium binding affinities and in vitro functional antagonism of CB1 receptor antagonists have been studied in detail, little is known on the kinetics of their receptor interaction. In this study, we therefore conducted kinetic assays for nine 1-(4,5-diarylthiophene-2-carbonyl)-4-phenylpiperidine-4-carboxamide derivatives and included the CB1 antagonist rimonabant as a comparison. For this we newly developed a dual-point competition association assay with [³H]CP55940 as the radioligand. This assay yielded Kinetic Rate Index (KRI) values from which structure-kinetics relationships (SKR) of hCB1 receptor antagonists could be established. The fast dissociating antagonist 6 had a similar receptor residence time (RT) as rimonabant, i.e. 19 and 14 min, respectively, while the slowest dissociating antagonist (9) had a very long RT of 2222 min, i.e. pseudo-irreversible dissociation kinetics. In functional assays, 9 displayed insurmountable antagonism, while the effects of the shortest RT antagonist 6 and rimonabant were surmountable. Taken together, this study shows that hCB1 receptor antagonists can have very divergent RTs, which are not correlated to their equilibrium affinities. Furthermore, their RTs appear to define their mode of functional antagonism, i.e. surmountable vs. insurmountable. Finally, based on the recently resolved hCB1 receptor crystal structure, we propose that the differences in RT can be explained by a different binding mode of antagonist 9 from short RT antagonists that is able to displace unfavorable water molecules. Taken together, these findings are of importance for future design and evaluation of potent and safe hCB1 receptor antagonists.