Abstract
The Light-Dependent Protochlorophyllide Oxidoreductase (LPOR) catalyzes a crucial step in chlorophyll biosynthesis: the rare biological photocatalytic reduction of the double C[double bond, length as m-dash]C bond in the precursor, protochlorophyllide (Pchlide). Despite its fundamental significance, limited structural insights into the active complex have hindered understanding of its reaction mechanism. Recently, a high-resolution cryo-EM structure of LPOR in its active conformation challenged our view of pigment binding, residue interactions, and the catalytic process. Surprisingly, this structure contrasts markedly with previous assumptions, particularly regarding the orientation of the bound Pchlide. To gain insights into the substrate binding puzzle, we conducted molecular dynamics simulations, quantum-mechanics/molecular-mechanics (QM/MM) calculations, and site-directed mutagenesis. Two Pchlide binding modes were considered, one aligning with historical proposals (mode A) and another consistent with the recent experimental data (mode B). Binding energy calculations revealed that in contrast to the non-specific interactions found for mode A, mode B exhibits distinct stabilizing interactions that support more thermodynamically favorable binding. A comprehensive analysis incorporating QM/MM-based local energy decomposition unraveled a complex interaction network involving Y177, H319, and the C13(1) carboxy group, influencing the pigment's excited state energy and potentially contributing to substrate specificity. Importantly, our results uniformly favor mode B, challenging established interpretations and emphasizing the need for a comprehensive re-evaluation of the LPOR reaction mechanism in a way that incorporates accurate structural information on pigment interactions and substrate-cofactor positioning in the binding pocket. The results shed light on the intricacies of LPOR's catalytic mechanism and provide a solid foundation for further elucidating the secrets of chlorophyll biosynthesis.