Molecular docking and transcriptomic analysis reveal the mechanism of myosin-derived peptides activating bitter receptor of hT2R1

To better understand the bitterness effect and molecule mechanism of myosin-derived peptides activating bitter receptors, the interaction between myosin-derived peptides of dry-cured ham and bitter receptors was investigated by molecular docking and molecular dynamics simulation; the signal transduction mechanism of myosin-derived peptides was explored by HEK-293T cells using calcium imaging and transcriptomics analysis. Lower CDOCKER energy was observed during the interaction between myosin-derived peptides and hT2R1 by molecular docking compared with hT2R4, hT2R5, hT2R8, hT2R14 and hT2R16. Hydrogen bonds and hydrophobic interaction were the most important interaction forces which stabilized the interaction of hT2R1 and myosin-derived peptides. Compared with LEKEKSELK and TEELEEAKK, the RMSF values and EC₅₀ values of HVLATLGEK were lower, indicating that hT2R1 was more sensitive to HVLATLGEK stimulation. Transcriptomics and KEGG analyses showed that 767 differentially expressed genes were found and mainly involved in cAMP signaling pathway, neuroactive ligand-receptor interaction, calcium signaling pathway and MAPK signaling pathway after stimulating of HVLATLGEK. Protein-protein interaction network further demonstrated that DDIT3, FOS, FOSB, MYC, EGR1 and CCN2 were the key genes to connect the six functional clusters including ligand-receptor interaction and signal transduction.