Monoacylglycerol lipase (MAGL) can regulate the endocannabinoid system and thus becomes a target of antidepressant drugs. In this paper, molecular docking and molecular dynamics simulations, combined with binding free energy calculation, were employed to investigate the inhibitory mechanism and binding modes of four aryl formyl piperidine derivative inhibitors with different 1-substituents to MAGL. The results showed that in the four systems, the main four regions where the enzyme bound to the inhibitor included around the head aromatic ring, the head carbonyl oxygen, the tail amide bond, and the tail benzene ring. The significant conformational changes in the more flexible lid domain of the enzyme were caused by 1-substituted group differences of inhibitors and resulted in different degrees of flipping in the tail of the inhibitor. The flipping led to a different direction of the tail amide bond and made a greater variation in its interaction with some of the charged residues in the enzyme, which further contributed to a different swing of the tail benzene ring. If the swing is large enough, it can weaken the binding strength of the head carbonyl oxygen to its nearby residues, and even the whole inhibitor with the enzyme so that the inhibition decreases.
Keywords: MM/PBSA; depression; inhibitor; molecular dynamics simulations; monoacylglycerol lipase.