Search

The binding kinetic properties of potential drugs may significantly influence their subsequent clinical efficacy. Predictions of these properties based on computer simulations provide a useful alternative to their expensive and time-demanding experimental counterparts, even at an early drug discovery stage.Herein, we perform Scaled Molecular Dynamics (ScaledMD) simulations on a set of 27 ligands of HSP90 belonging to more than 7 chemical series in order to estimate their relative residence time. We introduce two new techniques for the analysis and the classification of the simulated unbinding trajectories. The first technique, which helps in estimating the limits of the free energy well around the bound state and the second one, based on a new contact map fingerprint, allows the description and the comparison of the paths that lead to unbinding.Using these analyses, we find that ScaledMD’s relative residence time generally enables the identification of the slowest unbinders. We propose an explanation for the underestimation of the residence times of a subset of compounds and we investigate how the biasing in ScaledMD can affect the mechanistic insights that can be gained from the simulations.

There is a growing appreciation of the importance of drug-target binding kinetics for lead optimization. For G protein-coupled receptors (GPCRs), which mediate signaling over a wide range of time scales, the drug dissociation rate is often a better predictor of in vivo efficacy than binding affinity, although it is more challenging to compute. Here, we assess the ability of the τ-Random Acceleration Molecular Dynamics (τRAMD) approach to reproduce relative residence times and reveal dissociation mechanisms and the effects of allosteric modulation for two important membrane-embedded drug targets: the β2-adrenergic receptor and the muscarinic acetylcholine receptor M2. The dissociation mechanisms observed in the relatively short RAMD simulations (in which molecular dynamics (MD) simulations are performed using an additional force with an adaptively assigned random orientation applied to the ligand) are in general agreement with much more computationally intensive conventional MD and metadynamics simulations. Remarkably, although decreasing the magnitude of the random force generally reduces the number of egress routes observed, the ranking of ligands by dissociation rate is hardly affected and agrees well with experiment. The simulations also reproduce changes in residence time due to allosteric modulation and reveal associated changes in ligand dissociation pathways.

The recent and growing evidence that the efficacy of a drug can be correlated to target binding kinetics has seeded the development of a multitude of novel methods aimed at computing rate constants for receptor-ligand binding processes, as well as gaining an understanding of the binding and unbinding pathways and the determinants of structure-kinetic relationships. These new approaches include various types of enhanced sampling molecular dynamics simulations and the combination of energy-based models with chemometric analysis. We assess these approaches in the light of the varying levels of complexity of protein-ligand binding processes.

The TRAnsient Pockets in Proteins (TRAPP) webserver provides an automated workflow that allows users to explore the dynamics of a protein binding site and to detect pockets or sub-pockets that may transiently open due to protein internal motion. These transient or cryptic sub-pockets may be of interest in the design and optimization of small molecular inhibitors for a protein target of interest. The TRAPP workflow consists of the following three modules: (i) TRAPP structure— generation of an ensemble of structures using one or more of four possible molecular simulation methods; (ii) TRAPP analysis—superposition and clustering of the binding site conformations either in an ensemble of structures generated in step (i) or in PDB structures or trajectories uploaded by the user; and (iii) TRAPP pocket—detection, analysis, and visualization of the binding pocket dynamics and characteristics, such as volume, solvent-exposed area or properties of surrounding residues. A standard sequence conservation score per residue or a differential score per residue, for comparing on- and off-targets, can be calculated and displayed on the binding pocket for an uploaded multiple sequence alignment file, and known protein sequence annotations can be displayed simultaneously. The TRAPP webserver is freely available at http://trapp.h-its.org.

The addition of three N-terminal histidines to β-lactamase inhibitor protein was shown experimentally to increase its binding potency to an Au(111) surface substantially but the binding mechanism was not resolved. Here, we propose a complete adsorption mechanism for this fusion protein by means of a multi-scale simulation approach and free energy calculations. We find that adsorption is a three-step process: (i) recognition of the surface predominantly by the histidine fusion peptide and formation of an encounter complex facilitated by a reduced dielectric screening of water in the interfacial region, (ii) adsorption of the protein on the surface and adoption of a specific binding orientation, and (iii) adaptation of the protein structure on the metal surface accompanied by induced fit. We anticipate that the mechanistic features of protein adsorption to an Au(111) surface revealed here can be extended to other inorganic surfaces and proteins and will therefore aid the design of specific protein-surface interactions.

Drug-target residence time (τ), one of the main determinants of drug efficacy, remains highly challenging to predict computationally and, therefore, is usually not considered in the early stages of drug design. Here, we present an efficient computational method, τ-random acceleration molecular dynamics (τRAMD), for the ranking of drug candidates by their residence time and obtaining insights into ligand-target dissociation mechanisms. We assessed τRAMD on a data set of 70 diverse drug-like ligands of the N-terminal domain of HSP90α, a pharmaceutically important target with a highly flexible binding site, obtaining computed relative residence times with an accuracy of about 2.3τ for 78% of the compounds and less than 2.0τ within congeneric series. Analysis of dissociation trajectories reveals features that affect ligand unbinding rates, including transient polar interactions and steric hindrance. These results suggest that τRAMD will be widely applicable as a computationally efficient aid to improving drug residence times during lead optimization.

Structure-based drug design has often been restricted by the rather static picture of protein–ligand complexes presented by crystal structures, despite the widely accepted importance of protein flexibility in biomolecular recognition. Here we report a detailed experimental and computational study of the drug target, human heat shock protein 90, to explore the contribution of protein dynamics to the binding thermodynamics and kinetics of drug-like compounds. We observe that their binding properties depend on whether the protein has a loop or a helical conformation in the binding site of the ligand-bound state. Compounds bound to the helical conformation display slow association and dissociation rates, high-affinity and high cellular efficacy, and predominantly entropically driven binding. An important entropic contribution comes from the greater flexibility of the helical relative to the loop conformation in the ligand-bound state. This unusual mechanism suggests increasing target flexibility in the bound state by ligand design as a new strategy for drug discovery., An understanding of the dynamics of drug binding and unbinding processes is important for drug discovery. Here, the authors give insights into the binding mechanism of small drug-like molecules to human Hsp90 by combining thermodynamics and kinetics studies as well as molecular dynamics simulations.

Simulations of the long-time scale motions of a ligand binding pocket in a protein may open up new perspectives for the design of compounds with steric or chemical properties differing from those of known binders. However, slow motions of proteins are difficult to access using standard molecular dynamics (MD) simulations and are thus usually neglected in computational drug design. Here, we introduce two nonequilibrium MD approaches to identify conformational changes of a binding site and detect transient pockets associated with these motions. The methods proposed are based on the rotamerically induced perturbation (RIP) MD approach, which employs perturbation of side-chain torsional motion for initiating large-scale protein movement. The first approach, Langevin-RIP (L-RIP), entails a series of short Langevin MD simulations, each starting with perturbation of one of the side-chains lining the binding site of interest. L-RIP provides extensive sampling of conformational changes of the binding site. In less than 1 ns of MD simulation with L-RIP, we observed distortions of the α-helix in the ATP binding site of HSP90 and flipping of the DFG loop in Src kinase. In the second approach, RIPlig, a perturbation is applied to a pseudoligand placed in different parts of a binding pocket, which enables flexible regions of the binding site to be identified in a small number of 10 ps MD simulations. The methods were evaluated for four test proteins displaying different types and degrees of binding site flexibility. Both methods reveal all transient pocket regions in less than a total of 10 ns of simulations, even though many of these regions remained closed in 100 ns conventional MD. The proposed methods provide computationally efficient tools to explore binding site flexibility and can aid in the functional characterization of protein pockets, and the identification of transient pockets for ligand design.

We present TRAPP (TRAnsient Pockets in Proteins), a new automated software platform for tracking, analysis, and visualization of binding pocket variations along a protein motion trajectory or within an ensemble of protein structures that may encompass conformational changes ranging from local side chain fluctuations to global backbone motions. TRAPP performs accurate grid-based calculations of the shape and physicochemical characteristics of a binding pocket for each structure and detects the conserved and transient regions of the pocket in an ensemble of protein conformations. It also provides tools for tracing the opening of a particular subpocket and residues that contribute to the binding site. TRAPP thus enables an assessment of the druggability of a disease-related target protein taking its flexibility into account.

Protein-nanoparticle associations have important applications in nanoscience and nanotechnology such as targeted drug delivery and theranostics. However, the mechanisms by which proteins recognize nanoparticles and the determinants of specificity are still poorly understood at the microscopic level. Gold is a promising material in nanoparticles for nanobiotechnology applications because of the ease of its functionalization and its tunable optical properties. Ubiquitin is a small, cysteine-free protein (ubiquitous in eukaryotes) whose binding to gold nanoparticles has been characterized recently by nuclear magnetic resonance (NMR). To reveal the molecular basis of these protein-nanoparticle interactions, we performed simulations at multiple levels (ab initio quantum mechanics, classical molecular dynamics and Brownian dynamics) and compared the results with experimental data (circular dichroism and NMR). The results provide a model of the ensemble of structures constituting the ubiquitin-gold surface complex, and insights into the driving forces for the binding of ubiquitin to gold nanoparticles, the role of nanoparticle surfactants (citrate) in the association process, and the origin of the perturbations in the NMR chemical shifts.

Due to the contribution of drug-target binding kinetics to drug efficacy, there is a high level of interest in developing methods to predict drug-target binding kinetic parameters. During the review period, a wide range of enhanced sampling molecular dynamics simulation-based methods has been developed for computing drug-target binding kinetics and studying binding and unbinding mechanisms. Here, we assess the performance of these methods considering two benchmark systems in detail: mutant T4 lysozyme-ligand complexes and a large set of N-HSP90-inhibitor complexes. The results indicate that some of the simulation methods can already be usefully applied in drug discovery or lead optimization programs but that further studies on more high-quality experimental benchmark datasets are necessary to improve and validate computational methods., Figure 3 was improved. A definition of PIB was included. Reference to WE was added (ref. 20), reference to RAMD was corrected (ref. 43)

The dissociation of ligands from proteins and other biomacromolecules occurs over a wide range of timescales. For most pharmaceutically relevant inhibitors, these timescales are far beyond those that are accessible by conventional molecular dynamics (MD) simulation. Consequently, to explore ligand egress mechanisms and compute dissociation rates, it is necessary to enhance the sampling of ligand unbinding. Random Acceleration MD (RAMD) is a simple method to enhance ligand egress from a macromolecular binding site, which enables the exploration of ligand egress routes without prior knowledge of the reaction coordinates. Furthermore, the tauRAMD procedure can be used to compute the relative residence times of ligands. When combined with a machine-learning analysis of protein-ligand interaction fingerprints (IFPs), molecular features that affect ligand unbinding kinetics can be identified. Here, we describe the implementation of RAMD in GROMACS 2020, which provides significantly improved computational performance, with scaling to large molecular systems. For the automated analysis of RAMD results, we developed MD-IFP, a set of tools for the generation of IFPs along unbinding trajectories and for their use in the exploration of ligand dynamics. We demonstrate that the analysis of ligand dissociation trajectories by mapping them onto the IFP space enables the characterization of ligand dissociation routes and metastable states. The combined implementation of RAMD and MD-IFP provides a computationally efficient and freely available workflow that can be applied to hundreds of compounds in a reasonable computational time and will facilitate the use of tauRAMD in drug design., Comment: 35 pages, 5 figures in the main text and 9 figures in the supplementary information. The following article has been submitted to The Journal of Chemical Physics. After it is published, it will be found at https://publishing.aip.org/resources/librarians/products/journals/