Search

Chagas' disease is a parasitic infection widely distributed throughout Latin America, with devastating consequences in terms of human morbidity and mortality. Cruzain, the major cysteine protease from Trypanosoma cruzi, is an attractive target for antitrypanosomal chemotherapy. In the present work, classical two-dimensional quantitative structure-activity relationships (2D QSAR) and hologram QSAR (HQSAR) studies were performed on a training set of 45 thiosemicarbazone and semicarbazone derivatives as inhibitors of T. cruzi cruzain. Significant statistical models (HQSAR, q(2) = 0.75 and r(2) = 0.96; classical QSAR, q(2) = 0.72 and r(2) = 0.83) were obtained, indicating their consistency for untested compounds. The models were then used to evaluate an external test set containing 10 compounds which were not included in the training set, and the predicted values were in good agreement with the experimental results (HQSAR, r(2)(pred) = 0.95; classical QSAR, r(2)(pred) = 0.91), indicating the existence of complementary between the two ligand-based drug design techniques.

We report here the first crystal structure of a stable isosteric analogue of 1,3-bisphospho-d-glyceric acid (1,3-BPGA) bound to the catalytic domain of Trypanosoma cruzi glycosomal glyceraldehyde-3-phosphate dehydrogenase (gGAPDH) in which the two phosphoryl moieties interact with Arg249. This complex possibly illustrates a step of the catalytic process by which Arg249 may induce compression of the product formed, allowing its expulsion from the active site. Structural modifications were introduced into this isosteric analogue and the respective inhibitory effects of the resulting diphosphorylated compounds on T. cruzi and Trypanosoma brucei gGAPDHs were investigated by enzymatic inhibition studies, fluorescence spectroscopy, site-directed mutagenesis, and molecular modelling. Despite the high homology between the two trypanomastid gGAPDHs (> 95%), we have identified specific interactions that could be used to design selective irreversible inhibitors against T. cruzi gGAPDH.

Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) catalyzes the reversible oxidative phosphorylation of d-glyceraldehyde 3-phosphate (GAP) into d-glycerate 1,3-bisphosphate (1,3-diPG) in the presence of NAD(+) and inorganic phosphate (P(i)). Within the active site, two anion-binding sites were ascribed to the binding of the C3 phosphate of GAP (P(s)) and to the binding of the attacking phosphate ion (P(i)). The role played by these two sites in the catalytic mechanism in connection with the functional role of coenzyme exchange (NADH-NAD(+) shuttle) has been investigated by several studies leading to the C3 phosphate flipping model proposed by Skarzynski et al. [Skarzynski, T., Moody, P. C., and Wonacott, A. J. (1987) J. Mol. Biol. 193, 171-187]. This model has not yet received direct confirmation. To gain further insight into the role of both sites, we synthesized irreversible inhibitors which form with the essential cysteine residue a thioacyl enzyme analogue of the catalytic intermediate. Here we report the refined glycosomal Trypanosoma cruzi GAPDH in complex with a covalently bound GAP analogue at an improved resolution of 2.0-2.5 A. For this holo-thioacyl enzyme complex, a flip-flop movement is clearly characterized, the change from the P(i) to the P(s) binding site being correlated with the coenzyme exchange step: the weaker interaction of the intermediate when bound at the P(s) site with the cofactor allows its release and also the binding of the inorganic phosphate for the next catalytic step. This result gives strong experimental support for the generally accepted flip-flop model of the catalytic mechanism in GAPDH.

Drug design strategies based on Comparative Molecular Field Analysis (CoMFA) have been used to predict the activity of new compounds. The major advantage of this approach is that it permits the analysis of a large number of quantitative descriptors and uses chemometric methods such as partial least squares (PLS) to correlate changes in bioactivity with changes in chemical structure. Because it is often difficult to rationalize all variables affecting the binding affinity of compounds using CoMFA solely, the program GRID was used to describe ligands in terms of their molecular interaction fields, MIFs. The program VolSurf that is able to compress the relevant information present in 3D maps into a few descriptors can treat these GRID fields. The binding affinities of a new set of compounds consisting of 13 coumarins, for one of which the three-dimensional ligand-enzyme bound structure is known, were studied. A final model based on the mentioned programs was independently validated by synthesizing and testing new coumarin derivatives. By relying on our knowledge of the real physical data (i.e., combining crystallographic and binding affinity results), it is also shown that ligand-based design agrees with structure-based design. The compound with the highest binding affinity was the coumarin chalepin, isolated from Rutaceae species, with an IC50 value of 55.5 microM towards the enzyme glyceraldehyde-3-phosphate dehydrogenase (gGAPDH) from glycosomes of the parasite Trypanosoma cruzi, the causative agent of Chagas' disease. The proposed models from GRID MIFs have revealed the importance of lipophilic interactions in modulating the inhibition, but without excluding the dependence on stereo-electronic properties as found from CoMFA fields.