Molecular Modeling and Molecular Dynamics Simulation Studies on the Selective Binding Mechanism of MTHFD2 Inhibitors

Methylenetetrahydrofolate dehydrogenase 2 (MTHFD2) is a mitochondrial enzyme that plays an important role in purinecarbon metabolism and thymidine biosynthesis. It has attracted broad interest as a novel therapeutic target for cancer. However, a major problem of current MTHFD2 inhibitors is their lack of selectivity and reactivity with its closest isoform, MTHFD1. Recently, the first selective MTHFD2 inhibitor, DS44960156, has been reported and it exhibits a more than 18-fold selectivity for MTHFD2 over MTHFD1. However, mechanism of DS44960156 selective binding to MTHFD2 over MTHFD1 is unknown. In this study, molecular docking, molecular dynamics (MD) simulations, molecular mechanics generalized born/surface area (MM_GBSA) binding free energy calculations, and analysis of the decomposition of binding free energies were used to investigate the selective binding mechanism of DS44960156 to the folate-binding site of MTHFD2 over MTHFD1. The results revealed that contributions from residues Gln100/Gln132, Val55/Asn87, and Gly237/Gly310 in the binding pocket of MTHFD1/MTHFD2 are the key factors responsible for the binding selectivity. These findings explain the selectivity of DS44960156 to MTHFD2 over MTHFD1, and may provide guidance for the future study and design of novel MTHFD2 inhibitors.