Investigation of the effect of homocysteinylation of substance P on its binding to the NK1 receptor using molecular dynamics simulation

Abstract

Substance P is a neurotransmitter or modulator in both the central and peripheral nervous systems. In this work, modifications of the lysine in SP by homocysteine and an acetyl group as well as the conformational dynamics of the native and modified SP peptides and their complexes with the NK1 receptor were studied via MD simulation. It was found that modifying SP stabilizes the peptide structure, but the modified SP peptides are less likely to bind to the NK1 receptor, so the resulting complexes are less stable. The RMSD of native SP (~0.33 nm) is about twice as large as that of the modified SP peptides (~0.18 nm), while the RMSD for the receptor complexed with native SP is ~0.3 nm, and that for the receptor complexed with either of the modified peptides is ~0.35 nm, which demonstrates the high stability of the modified SP peptides as well as the receptor complexed with native SP. Such behavior was also observed in other structural analyses. The binding free energies of the native and modified SP peptides with the NK1 receptor were also compared. The ΔG_bind values for the binding of homocysteinylated SP to the NK1 receptor and the binding of the acetylated SP and native SP to the NK1 receptor were -38.89, -64.46, and - 264.52 kJ mol^-1, respectively. Modification of the lysine of SP decreases the binding affinity of the peptide to the NK1 receptor. In other words, homocysteinylation or acetylation of SP leads to weaker interactions of the peptide with the NK1 receptor compared to those between native SP and NK1. We propose that this phenomenon leads to increased levels of homocysteinylated SP in plasma in many diseases such as breast cancer. Graphical abstract Substance P (SP) is a neuropeptide which binds to the NK1 receptor. SP is of great pharmacological interest, as agonists and antagonists of SP can potentially be used to treat many chronic diseases. Therefore, in this work, the lysine (LYS) in SP was theoretically modified with a homocysteine or acetyl group to explore the effects of such a modification on the binding affinity of this peptide with the NK1 receptor and the structural dynamics of the resulting complex.

Keywords: Acetyl; Homocysteine; Homocysteinylation; MD simulation; NK1 receptor; Substance P.