Ligand binding and protein dynamics in lactate dehydrogenase

Recent experimental studies suggest that lactate dehydrogenase (LDH) binds its substrate via the formation of a LDH/NADH·substrate encounter complex through a select-fit mechanism, whereby only a minority population of LDH/NADH is binding-competent. In this study, we perform molecular dynamics calculations to explore the variations in structure accessible to the binary complex with a focus on identifying structures that seem likely to be binding-competent and which are in accord with the known experimental characterization of forming binding-competent species. We find that LDH/NADH samples quite a range of protein conformations within our 2.148 ns calculations, some of which yield quite facile access of solvent to the active site. The results suggest that the mobile loop of LDH is perhaps just partially open in these conformations and that multiple open conformations, yielding multiple binding pathways, are likely. These open conformations do not require large-scale unfolding/melting of the binary complex. Rather, open versus closed conformations are due to subtle protein and water rearrangements. Nevertheless, the large heat capacity change observed between binding-competent and binding-incompetent can be explained by changes in solvation and an internal rearrangement of hydrogen bonds. We speculate that such a strategy for binding may be necessary to get a ligand efficiently to a binding pocket that is located fairly deep within the protein's interior.