A boosted unbiased molecular dynamics method for predicting ligands binding mechanisms: probing the binding pathway of dasatinib to Src-kinase

Small molecules such as metabolites and drugs play essential roles in biological processes and pharmaceutical industry. Knowing their interactions with biomacromolecular targets demands a deep understanding of binding mechanisms. Dozens of papers have suggested that discovering of the binding event by means of conventional unbiased molecular dynamics (MD) simulation urges considerable amount of computational resources, therefore, only one who holds a cluster or a supercomputer can afford such extensive simulations. Thus, many researchers who do not own such resources are reluctant to take the benefits of running unbiased molecular dynamics simulation, in full atomistic details, when studying a ligand binding pathway. Many researchers are impelled to be content with biased molecular dynamics simulations which seek its validation due to its intrinsic preconceived framework. In this work, we have presented a workable stratagem to encourage everyone to perform unbiased (unguided) molecular dynamics simulations, in this case a protein-ligand binding process, by typical desktop computers and so achieve valuable results in nanosecond time scale. Here, we have described a dynamical binding’s process of an anticancer drug, the dasatinib, to the c-Src kinase in full atomistic details for the first time, without applying any biasing force or potential which may lead the drug to artificial interactions with the protein. We have attained multiple independent binding events which occurred in the nano-second timescales, surprisingly as little as ∼30 ns. Both the protonated and deprotonated forms of the dasatinib reached the crystallographic binding mode without having any major intermediate state during induction.