Homology modeling and in silico site directed mutagenesis of pyruvate ferredoxin oxidoreductase from Clostridium thermocellum

Pyruvate ferredoxin oxidoreductase is the crucial enzyme that involves in bioethanol synthesis pathway of Clostridium thermocellum. It is an ethanologenic organism but has been investigated less on its enzyme structure. The amino acid sequence of Pyruvate ferredoxin oxidoreductase was derived from UNIPROT and the screened crystal structure was taken as the template for homology modeling using MODELLER 9V11. The model was loop refined and was validated using RMSD, ProSA and PROCHECK. The docking and per residue interaction studies were carried out to elucidate the interaction energies of amino acid residues with pyruvate. To enhance the binding of pyruvate with the enzyme, mutation studies were carried out by replacing Thr31 as it had a less interaction energy. Out of 10 mutants, T31N, T31Q and T31G were selected using potential energy and the residual energy calculations. Five nanoseconds explicit MD simulations were run for apo, wild type and mutants T31N, T31Q and T31G using Desmond. RMSD, RMSF, distance plots and H-bonds analysis proved T31G to be a favorable mutant for binding of pyruvate. Thus, modeling PFOR would help in profound understanding of its structural clefts and mutation studies would aid in improving the enzyme efficiency.