Computational blind docking approach was used for mapping of possible binding sites in L-type pyruvate kinase subunit for peptides, RRASVA and the phosphorylated derivative RRAS(Pi)VA, which model the phosphorylatable N-terminal regulatory domain of the enzyme. In parallel, the same docking analysis was done for both substrates of this enzyme, phosphoenolpyruvate (PEP) and adenosine diphosphate (ADP), and for docking of fructose 1,6-bisphosphate (FBP), which is the allosteric activator of the enzyme. The binding properties of the entire surface of the protein were scanned and several possible binding sites were identified in domains A and C of the protein, while domain B revealed no docking sites for peptides or for substrates or the allosteric regulator. It was found that the docking sites of different ligands were partially overlapping, pointing to the possibility that some regulatory effects, observed in the case of L-type pyruvate kinase, may be caused by the competition of different ligands for the same binding sites. 
Computational simulation of ligand docking to L-type pyruvate kinase subunit