Abstract
Methylmalonic aciduria and homocystinuria type C protein (MMACHC) is required by the body to metabolize cobalamin (Cbl). Due to its complex structure and cofactor forms, Cbl passes through an extensive series of absorptive and processing steps before being delivered to mitochondrial methyl malonyl-CoA mutase and cytosolic methionine synthase. Depending on the cofactor attached, MMACHC performs either flavin-dependent reductive decyanation or glutathione (GSH)-dependent dealkylation. The alkyl groups of Cbl have to be removed in the presence of GSH to produce intermediates that can later be converted into active cofactor forms. Pathogenic mutations in the GSH binding site, such as R161Q, R161G, R206P, R206W, and R206Q, have been reported to cause Cbl diseases. The impact of these variations on MMACHC's structure and how it affects GSH and Cbl binding at the molecular level is poorly understood. To better understand the molecular basis of this interaction, mutant structures involving the MMACHC-MeCbl-GSH complex were generated using in silico site-directed point mutations and explored using molecular dynamics (MD) simulations. The results revealed that mutations in the key arginine residues disrupt GSH binding by breaking the interactions and reducing the free energy of binding of GSH. Specifically, variations at position 206 appeared to produce weaker GSH binding. The lowered binding affinity for GSH in the variant structures could impact metabolic pathways involving Cbl and its trafficking.