Binding Mechanism of Inhibitors to Heat Shock Protein 90 Investigated by Multiple Independent Molecular Dynamics Simulations and Prediction of Binding Free Energy

Abstract

The heat shock protein (HSP90) has been an import target of drug design in the treatment of human disease. An exploration of the conformational changes in HSP90 can provide useful information for the development of efficient inhibitors targeting HSP90. In this work, multiple independent all-atom molecular dynamics (AAMD) simulations followed by calculations of the molecular mechanics generalized Born surface area (MM-GBSA) were performed to explore the binding mechanism of three inhibitors (W8Y, W8V, and W8S) to HSP90. The dynamics analyses verified that the presence of inhibitors impacts the structural flexibility, correlated movements, and dynamics behavior of HSP90. The results of the MM-GBSA calculations suggest that the selection of GB models and empirical parameters has important influences on the predicted results and verify that van der Waals interactions are the main forces that determine inhibitor–HSP90 binding. The contributions of separate residues to the inhibitor–HSP90 binding process indicate that hydrogen-bonding interactions (HBIs) and hydrophobic interactions play important roles in HSP90–inhibitor identifications. Moreover, residues L34, N37, D40, A41, D79, I82, G83, M84, F124, and T171 are recognized as hot spots of inhibitor–HSP90 binding and provide significant target sites of for the design of drugs related to HSP90. This study aims to contribute to the development of efficient inhibitors that target HSP90 by providing an energy-based and theoretical foundation.

Keywords

• HSP90
• multiple all-atom molecular dynamics simulations
• hot interaction spots
• MM-GBSA
• principal component analysis