Protein conformational flexibility modulates kinetics and thermodynamics of drug binding

Structure-kinetic relationships that control the residence time of drug-target complexes: insights from molecular structure and dynamics

Time-dependent target occupancy is a function of both the thermodynamics and kinetics of drug-target interactions. However, while the optimization of thermodynamic affinity through approaches such as structure-based drug design is now relatively straight forward, less is understood about the molecular interactions that control the kinetics of drug complex formation and breakdown since this depends on both the ground and transition state energies on the binding reaction coordinate. In this opinion we highlight several recent examples that shed light on current approaches that are elucidating the factors that control the life-time of the drug-target complex.

Efficiently Photocontrollable or Not? Biological Activity of Photoisomerizable Diarylethenes

Diarylethene derivatives, the biological activity of which can be reversibly changed by irradiation with light of different wavelengths, have shown promise as scientific tools and as candidates for photocontrollable drugs. However, examples demonstrating efficient photocontrol of their biological activity are still relatively rare. This concept article discusses the possible reasons for this situation and presents a critical analysis of existing data and hypotheses in this field, in order to extract the design principles enabling the construction of efficient photocontrollable diarylethene-based molecules. Papers addressing biologically relevant interactions between diarylethenes and biomolecules are analyzed; however, in most published cases, the efficiency of photocontrol in living systems remains to be demonstrated. We hope that this article will encourage further discussion of design principles, primarily among pharmacologists, synthetic and medicinal chemists.

Ligand-Induced Conformational Changes in HSP90 Monitored Time Resolved and Label Free-Towards a Conformational Activity Screening for Drug Discovery

Investigation of protein–ligand interactions is crucial during early drug‐discovery processes. ATR‐FTIR spectroscopy can detect label‐free protein–ligand interactions with high spatiotemporal resolution. Here we immobilized, as an example, the heat shock protein HSP90 on an ATR crystal. This protein is an important molecular target for drugs against several diseases including cancer. With our novel approach we investigated a ligand‐induced secondary structural change. Two specific binding modes of 19 drug‐like compounds were analyzed. Different binding modes can lead to different efficacy and specificity of different drugs. In addition, the k_obs values of ligand dissociation were obtained. The results were validated by X‐ray crystallography for the structural change and by SPR experiments for the dissociation kinetics, but our method yields all data in a single and simple experiment.

Estimation of Drug-Target Residence Times by τ-Random Acceleration Molecular Dynamics Simulations

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.

Kinetics of Huperzine A Dissociation from Acetylcholinesterase via Multiple Unbinding Pathways

The dissociation of huperzine A (hupA) from Torpedo californica acetylcholinesterase ( TcAChE) was investigated by 4 μs unbiased and biased all-atom molecular dynamics (MD) simulations in explicit solvent. We performed our study using memetic sampling (MS) for the determination of reaction pathways (RPs), metadynamics to calculate free energy, and maximum-likelihood estimation (MLE) to recover kinetic rates from unbiased MD simulations. Our simulations suggest that the dissociation of hupA occurs mainly via two RPs: a front door along the axis of the active-site gorge (pwf) and through a new transient side door (pws), i.e., formed by the Ω-loop (residues 67-94 of TcAChE). An analysis of the inhibitor unbinding along the RPs suggests that pws is opened transiently after hupA and the Ω-loop reach a low free-energy transition state characterized by the orientation of the pyridone group of the inhibitor directed toward the Ω-loop plane. Unlike pws, pwf does not require large structural changes in TcAChE to be accessible. The estimated free energies and rates agree well with available experimental data. The dissociation rates along the unbinding pathways are similar, suggesting that the dissociation of hupA along pws is likely to be relevant. This indicates that perturbations to hupA- TcAChE interactions could potentially induce pathway hopping. In summary, our results characterize the slow-onset inhibition of TcAChE by hupA, which may provide the structural and energetic bases for the rational design of the next-generation slow-onset inhibitors with optimized pharmacokinetic properties for the treatment of Alzheimer's disease.

Dynamics revelation of conformational changes and binding modes of heat shock protein 90 induced by inhibitor associations

Heat shock protein 90 (Hsp90) has been an attractive target of potential drug design for antitumor treatment. The current work integrates molecular dynamics (MD) simulations, calculations of binding free energy, and principal component (PC) analysis with scanning of inhibitor–residue interaction to probe the binding modes of inhibitors YK9, YKJ and YKI to Hsp90 and identify the hot spot of the inhibitor–Hsp90 binding. The results suggest that the introductions of two groups G1 and G2 into YKJ and YKI strengthen the binding ability of YKJ and YKI to Hsp90 compared to YK9. PC analysis based MD trajectories prove that inhibitor bindings exert significant effects on the conformational changes, internal dynamics and motion modes of Hsp90, especially for the helix α2 and the loops L1 and L2. The calculations of residue-based free energy decomposition and scanning of the inhibitor–Hsp90 interaction suggest that six residues L107, G108, F138, Y139, W162 and F170 construct the common hot spot of the inhibitor–residue interactions. Moreover the substitutions of the groups G1 and G2 in YKJ and YKI lead to two additional hydrogen bonding interactions and multiple hydrophobic interactions for bindings of YKJ and YKI to Hsp90. This work is also expected to contribute theoretical hints for the design of potent inhibitors toward Hsp90.

Heat shock protein 90 (Hsp90) has been an attractive target of potential drug design for antitumor treatment.