Homology Modeling and Molecular Dynamics Simulations of Trypanosoma cruzi Phosphodiesterase b1

Abstract

Cyclic nucleotide phosphodiesterases have been implicated in the proliferation, differentiation and osmotic regulation of trypanosomatids; in some trypanosomatid species, they have been validated as molecular targets for the development of new therapeutic agents. Because the experimental structure of Trypanosoma cruzi PDEb1 (TcrPDEb1) has not been solved so far, an homology model of the target was created using the structure of Trypanosoma brucei PDEb1 (TbrPDEb1) as a template. The model was refined by extensive enhanced sampling molecular dynamics simulations, and representative snapshots were extracted from the trajectory by combined clustering analysis. This structural ensemble was used to develop a structure-based docking model of the target. The docking accuracy of the model was validated by redocking and cross-docking experiments using all available crystal structures of TbrPDEb1, whereas the scoring accuracy was validated through a retrospective screen, using a carefully curated dataset of compounds assayed against TbrPDEb1 and/or TcrPDEb1. Considering the results from in silico validations, the model may be applied in prospective virtual screening campaigns to identify novel hits, as well as to guide the rational design of potent and selective inhibitors targeting this enzyme.

Keywords: Chagas disease; Trypanosoma cruzi; homology modeling; molecular docking; molecular dynamics; phosphodiesterase.