title: Exploring the Origin of Differential Binding Affinities of Human Tubulin Isotypes αβII, αβIII and αβIV for DAMA-Colchicine Using Homology Modelling, Molecular Docking and Molecular Dynamics Simulations.
abstract: Tubulin isotypes are found to play an important role in regulating microtubule dynamics. The isotype composition is also thought to contribute in the development of drug resistance as tubulin isotypes show differential binding affinities for various anti-cancer agents. Tubulin isotypes αβII, αβIII and αβIV show differential binding affinity for colchicine. However, the origin of differential binding affinity is not well understood at the molecular level. Here, we investigate the origin of differential binding affinity of a colchicine analogue N-deacetyl-N-(2-mercaptoacetyl)-colchicine (DAMA-colchicine) for human αβII, αβIII and αβIV isotypes, employing sequence analysis, homology modeling, molecular docking, molecular dynamics simulation and MM-GBSA binding free energy calculations. The sequence analysis study shows that the residue compositions are different in the colchicine binding pocket of αβII and αβIII, whereas no such difference is present in αβIV tubulin isotypes. Further, the molecular docking and molecular dynamics simulations results show that residue differences present at the colchicine binding pocket weaken the bonding interactions and the correct binding of DAMA-colchicine at the interface of αβII and αβIII tubulin isotypes. Post molecular dynamics simulation analysis suggests that these residue variations affect the structure and dynamics of αβII and αβIII tubulin isotypes, which in turn affect the binding of DAMA-colchicine. Further, the binding free-energy calculation shows that αβIV tubulin isotype has the highest binding free-energy and αβIII has the lowest binding free-energy for DAMA-colchicine. The order of binding free-energy for DAMA-colchicine is αβIV ≃ αβII >> αβIII. Thus, our computational approaches provide an insight into the effect of residue variations on differential binding of αβII, αβIII and αβIV tubulin isotypes with DAMA-colchicine and may help to design new analogues with higher binding affinities for tubulin isotypes.
authors.affiliation: Department of Biosciences and Bioengineering, Indian Institute of Technology Bombay, Powai, Mumbai-400076, Maharashtra, India.
authors.forename: Bajarang Vasant
authors.initials: BV
authors.lastname: Kumbhar
authors.affiliation: Department of Biosciences and Bioengineering, Indian Institute of Technology Bombay, Powai, Mumbai-400076, Maharashtra, India.
authors.forename: Anubhaw
authors.initials: A
authors.lastname: Borogaon
authors.affiliation: Department of Biosciences and Bioengineering, Indian Institute of Technology Bombay, Powai, Mumbai-400076, Maharashtra, India.
authors.forename: Dulal
authors.initials: D
authors.lastname: Panda
authors.affiliation: Department of Biosciences and Bioengineering, Indian Institute of Technology Bombay, Powai, Mumbai-400076, Maharashtra, India.
authors.forename: Ambarish
authors.initials: A
authors.lastname: Kunwar
mesh_terms: D000595:Amino Acid Sequence; D001665:Binding Sites; D003078:Colchicine / Q000737:chemistry* / Q000378:metabolism*; D006801:Humans; D008958:Models, Molecular; D062105:Molecular Docking Simulation; D056004:Molecular Dynamics Simulation; D011487:Protein Conformation; D020033:Protein Isoforms; D017386:Sequence Homology, Amino Acid; D014404:Tubulin / Q000737:chemistry* / Q000378:metabolism*
chemical_list: D020033:Protein Isoforms; D014404:Tubulin; D003078:Colchicine