In Silico docking, molecular dynamics and binding energy insights into the bolinaquinone-clathrin terminal domain binding site

Clathrin-mediated endocytosis (CME) is a process that regulates selective internalization of important cellular cargo using clathrin-coated vesicles. Perturbation of this process has been linked to many diseases including cancer and neurodegenerative conditions. Chemical proteomics identified the marine metabolite, 2-hydroxy-5-methoxy-3-(((1S,4aS,8aS)-1,4a,5-Trimethyl-1,2,3,4,4a,7, 8,8a-octahydronaphthalen-2-yl)methyl)cyclohexa-2,5-diene-1,4-dione (bolinaquinone) as a clathrin inhibitor. While being an attractive medicinal chemistry target, the lack of data about bolinaquinone's mode of binding to the clathrin enzyme represents a major limitation for its structural optimization. We have used a molecular modeling approach to rationalize the observed activity of bolinaquinone and to predict its mode of binding with the clathrin terminal domain (CTD). The applied protocol started by global rigid-protein docking followed by flexible docking, molecular dynamics and linear interaction energy calculations. The results revealed the potential of bolinaquinone to interact with various pockets within the CTD, including the clathrin-box binding site. The results also highlight the importance of electrostatic contacts over van der Waals interactions for proper binding between bolinaquinone and its possible binding sites. This study provides a novel model that has the potential to allow rapid elaboration of bolinaquinone analogues as a new class of clathrin inhibitors.