Application of the linear interaction energy method (LIE) to estimate the binding free energy values of Escherichia coli wild-type and mutant arginine repressor C-terminal domain (ArgRc)-l-arginine and ArgRc-l-citrulline protein-ligand complexes

Protein-ligand binding free energy values of wild-type and mutant C-terminal domain of Escherichia coli arginine repressor (ArgRc) protein systems bound to L-arginine or L-citrulline molecules were calculated using the linear interaction energy (LIE) method by molecular dynamics (MD) simulation. The binding behaviour predicted by the dissociation constant (K(d)) calculations from the binding free energy values showed preferences for binding of L-arginine to the wild-type ArgRc but not to the mutant ArgRc(D128N). On the other hand, L-citrulline do not favour binding to wild-type ArgRc but prefer binding to mutant ArgRc(D128N). The dissociation constant for the wild-type ArgRc-L-arginine complex obtained in this study is in agreement with reported experimental results. Our results also support the experimental data for the binding of L-citrulline to the mutant ArgRc(D128N). These showed that LIE method for protein-ligand binding free energy calculation could be applied to the wild-type and the mutant E. coli ArgRc-L-arginine and ArgRc-L-citrulline protein-ligand complexes and possibly to other transcriptional repressor-co-repressor systems as well.