Conformational States of E7010 Is Complemented by Microclusters of Water Inside the α,β-Tubulin Core

The alpha,beta-tubulin is the building block of microtubules, which is associated with and dissociated from the microtubular architecture complying with the dynamic instability of the microtubules. This dynamic instability has a direct relation with the spindle formation by the microtubules and cell division kinetics. E7010 is one of the promising ligands of an alpha,beta-tubulin protein that binds at the core of this protein and can diminish the protein's ability to fit to a growing microtubule, thus frustrating cell division. Although X-ray crystallography has reported a specific binding conformation of E7010 in PDB, molecular dynamics (MD) simulations have revealed two other conformational states of the ligand capable of binding to tubulin with stabilities close to that state reported in PDB. To rationalize this quasidegeneracy of ligand binding modes, MD simulations have further revealed that the understanding of the mechanism of E7010-tubulin binding remains incomplete unless the role of water molecules to bridge this interaction is taken into consideration, a very critical insight that was not visible from the PDB structure. Further, these water molecules differ from the standard examples of bridging waters which generally exist as isolated water molecules between the receptor and the ligand. In the present case, the water molecules sandwiched between ligand and protein, sequestered from the bulk solvent, integrate with each other by an H-bonds network forming a group, which appear as microclusters of water. The structural packing with the ligand binding pocket and the bridging interactions between protein and ligand take place through such clusters. The presence of this microcluster of water is not just cosmetic, instead they have a crucial impact on the ligand binding thermodynamics. Only with the explicit consideration of these water clusters in the binding energy calculations (MMGBSA) is the stability of the native mode of ligand binding reported in PDB rationalized. At the same time, two other binding modes are elucidated to be quasi-degenerate with the native state and that indicates the further possibility in gaining more entropic stabilization of the complex. The role of such bridging water clusters to enhance the protein-ligand interaction will be insightful for designing the next generation prospective compounds in the field of cancer therapeutics.