Is inhibition process better described with MD(QM/MM) simulations? The case of urokinase type plasminogen activator inhibitors.

Urokinase plasminogen activator (uPA) is an enzyme involved in cancer growth and metastasis. Therefore, the design of inhibitors of uPA is of high therapeutic value, and several chemical families have been explored, even if none has still emerged, emphasizing the need of a rationalized approach. This work represents a complete computational study of uPA complexed with five inhibitors, which present weak similarities. Molecular dynamics simulations in explicit solvent were conducted, and structural analyses, along with molecular mechanics (MM)/Poisson-Boltzmann surface area free energies estimations, yield precious structure-activity relationships of these inhibitors. Besides, we realized supplemental QM/MM computations that improved drastically the quality of our models providing original information on the hydrogen bonds and charge transfer effects, which are, most often, neglected in other studies. We suggest that these simulations and analyses could be reproduced for other systems involving protein/ligand molecular recognitions.