An integrated computational approach to the phenomenon of potent and selective inhibition of aurora kinases B and C by a series of 7-substituted indirubins.

A variation of the bromine substitution from 6- to 7-position converts the glycogen synthase kinase-3alpha/beta-(GSK-3-alpha/beta) selective inhibitor 6-bromoindirubin-3'-oxime (6BIO) to a potent inhibitor of Aurora B and C kinases. The novel indirubin analogue 7-bromoindirubin-3'-oxime (7BIO) demonstrated unexpected selectivity against these two kinases since the homologous kinase Aurora A was poorly inhibited. A hypothesis regarding the 7BIO selectivity profile was stated and validated by docking, molecular dynamics, and free energy perturbation calculations. The residue (Thr217AurA, Glu161AurB, Glu127AurC) located in the active site was identified as a major contributor to the enhanced affinity of 7BIO for Aurora B and C versus Aurora A. Furthermore, the docking events of 7BIO and several of its analogues were approached by quantitative models based on semiempirical scoring functions. In the course of model construction and optimization, a number of important factors influencing the quality of each model like the application of force constraints or the sampling method were determined. Among these factors, the presence and treatment of structurally important water molecules had a pronounced impact on the quality of each model. The final model was validated by use of free energy perturbation calculations.