Molecular basis for benzimidazole resistance from a novel β-tubulin binding site model

Abstract

Benzimidazole-2-carbamate derivatives (BzCs) are the most commonly used antiparasitic drugs for the treatment of protozoan and helminthic infections. BzCs inhibit the microtubule polymerization mechanism through binding selectively to the β-tubulin subunit in which mutations have been identified that lead to drug resistance. Currently, the lack of crystallographic structures of β-tubulin in parasites has limited the study of the binding site and the analysis of the resistance to BzCs. Recently, our research group has proposed a model to explain the interaction between the BzCs and a binding site in the β-tubulin. Herein, we report the homology models of two susceptible (Haemonchus contortus and Giardia intestinalis) parasites and one unsusceptible (Entamoeba histolytica) generated using the structure of the mammal Ovis aries, considered as a low susceptible organism, as a template. Additionally, the mechanism by which the principal single point mutations Phe167Tyr, Glu198Ala and Phe200Tyr could lead to resistance to BzCs is analyzed. Molecular docking and molecular dynamics studies were carried out in order to evaluate the models and the ligand–protein complexes’ behaviors. This study represents a first attempt towards understanding, at the molecular level, the structural composition of β-tubulin in all organisms, also suggesting possible resistance mechanisms. Furthermore, these results support the importance of benzimidazole derivative optimization in order to design more potent and selective (less toxic) molecules for the treatment of parasitic diseases.