Trypanosoma cruzi glutamate cysteine ligase (TcGCL) is considered a potential drug target to develop novel antichagasic drugs. We have used a variety of computational methods to investigate the interactions between TcGCL with Glutathione (GSH). The three-dimensional structure of TcGCL was constructed by comparative modeling methods using the Saccharomyces cerevisiae glutamate cysteine ligase as template. Molecular dynamics simulations were used to validate the TcGCL model and to analyze the molecular interactions with GSH. Using RMSD clustering, the most prevalent GSH binding modes were identified paying attention to the residues involved in the molecular interactions. The GSH binding modes were used to propose pharmacophore models that can be exploited in further studies to identify novel antichagasic compounds.