Abstract
Vanillin holds significant importance as a flavoring agent in various industries, including food, pharmaceuticals, and cosmetics. The CoA-dependent pathway for the biosynthesis of vanillin from ferulic acid involved feruloyl-CoA synthase (Fcs) and enoyl-CoA hydratase/lyase (Ech). In this research, the Fcs and Ech were derived from Streptomyces sp. strain V-1. The sequence conservation and structural features of Ech were analyzed by computational techniques including sequence alignment and molecular dynamics simulation. After detailed study for the major binding modes and key amino acid residues between Ech and substrates, a series of mutations (F74W, A130G, A130G/T132S, R147Q, Q255R, ΔT90, ΔTGPEIL, ΔN1-11, ΔC260-287) were obtained by rational design. Finally, the yield of vanillin produced by these mutants was verified by whole-cell catalysis. The results indicated that three mutants, F74W, Q147R, and ΔN1-11, showed higher yields than wild-type Ech. Molecular dynamics simulations and residue energy decomposition identified the basic residues K37, R38, K561, and R564 as the key residues affecting the free energy of binding between Ech and feruloyl-coenzyme A (FCA). The large changes in electrostatic interacting and polar solvating energies caused by the mutations may lead to decreased enzyme activity. This study provides important theoretical guidance as well as experimental data for the biosynthetic pathway of vanillin.