Unraveling the conformational dynamics of glycerol 3-phosphate dehydrogenase, a nicotinamide adenine dinucleotide-dependent enzyme of Leishmania mexicana.

Allosteric changes modulate the enzymatic activity, leading to activation or inhibition of the molecular target. Understanding the induced fit accommodation mechanism of a ligand in its lowest-free energy state and the subsequent conformational changes induced in the protein are important questions for drug design. In the present study, molecular dynamics (MD) simulations, binding free energy calculations, and principal component analysis (PCA) were applied to analyze the glycerol-3-phosphate dehydrogenase of Leishmania mexicana (LmGPDH) conformational changes induced by its cofactor and substrate binding. GPDH is a nicotinamide adenine dinucleotide (NAD)-dependent enzyme, which has been reported as an interesting target for drug discovery and development against leishmaniasis. Despite its relevance for glycolysis and pentose phosphate pathways, the structural flexibility and conformational motions of LmGPDH in complex with NADH and dihydroxyacetone phosphate (DHAP) remain unexplored. Here, we analyzed the conformational dynamics of the enzyme-NADH complex (cofactor), and the enzyme-NADH-DHAP complex (adduct), mapped the hydrogen-bond interactions for the complexes and pointed some structural determinants of the enzyme that emerge from these contacts to NADH and DHAP. Finally, we proposed a consistent mechanism for the conformational changes on the first step of the reversible redox conversion of dihydroxyacetone phosphate to glycerol 3-phosphate, indicating key residues and interactions that could be further explored in drug discovery.

Computational study of conformational changes in human 3-hydroxy-3-methylglutaryl coenzyme reductase induced by substrate binding.

The enzyme 3-hydroxy-3-methylglutaryl coenzyme A reductase (HMGR) is mainly involved in the regulation of cholesterol biosynthesis. HMGR catalyses the reduction of 3-hydroxy-3-methylglutaryl coenzyme A (HMG-CoA) to mevalonate at the expense of two NADPH molecules in a two-step reversible reaction. In the present study, we constructed a model of human HMGR (hHMGR) to explore the conformational changes of HMGR in complex with HMG-CoA and NADPH. In addition, we analysed the complete sequence of the Flap domain using molecular dynamics (MD) simulations and principal component analysis (PCA). The simulations revealed that the Flap domain plays an important role in catalytic site activation and substrate binding. The apo form of hHMGR remained in an open state, while a substrate-induced closure of the Flap domain was observed for holo hHMGR. Our study also demonstrated that the phosphorylation of Ser872 induces significant conformational changes in the Flap domain that lead to a complete closure of the active site, suggesting three principal conformations for the first stage of hHMGR catalysis. Our results were consistent with previous proposed models for the catalytic mechanism of hHMGR. Communicated by Ramaswamy H. Sarma.