Homology modeling of T. cruzi and L. major NADH-dependent fumarate reductases: Ligand docking, molecular dynamics validation, and insights on their binding modes

Abstract

Leishmania major and Trypanosoma cruzi are the main causes of leishmaniasis and Chagas disease, two endemic parasitosis identified as neglected diseases by the World Health Organization. Fumarate reductase (FR) is a central enzyme in the conversion of fumarate to succinate, an energy releasing path essential for the survival of these protozoans which is also absent in their mammalian hosts. FR can thus be considered as a good candidate for targeting specific inhibition by new drugs designed against L. major and T. cruzi. The lack of tertiary structures available for LmFR and TcFR has limited until now the possibility of performing structure-based drug design. Here we used homology modeling combined with enzyme-cofactor docking to propose tertiary structures for NADH-dependent LmFR and TcFR using an homologous X-ray crystallographic structure of flavine-adenine dinucleotide (FAD) dependent FR from Shewanella frigidimarina (PDB ID: 1QO8) as template. These models were refined and stabilized with/without substrate in the active site using classical molecular dynamics simulations under quasi-physiological conditions. Structural features relevant for understanding the mechanism of action of the enzyme were also analyzed, with special attention to the hydrogen bond network involving the cofactor and water molecules present at the binding sites. A small set of compounds previously synthesized and assayed for their inhibitory capacity against TcFR ([M(mpo)2] metal complexes with M=PtII, PdII and VIVO and mpo=2-mercaptopyridine N-oxide) and LmFR (licochalcone A) were screened by protein–ligand docking using the NADH-LmFR and NADH-TcFR models here proposed and validated, gaining insight into their binding modes in each enzyme.