A molecular dynamics study on opioid activities of biphalin molecule

Molecular dynamics simulations of the biphalin molecule, (Tyr-D-Ala-Gly-Phe-NH)₂, and the active tetrapeptide hydrazide, Tyr-D-Ala-Gly-Phe-NH-NH₂ were performed to investigate the cause of the increased μ and δ receptor binding affinities of the former over the latter. The simulation results demonstrate that the acylation of the two equal tetrapeptide fragments of biphalin produces the constrained hydrazide bridges C₄^a - C₄￠- N₉ - N₁₀ {\hbox{C}}_4^{\alpha } - {{\hbox{C}}_4}\prime - {{\hbox{N}}_9} - {{\hbox{N}}_{{10}}} and N₉ - N₁₀ - C₅￠- C₅^a {{\hbox{N}}_9} - {{\hbox{N}}_{{10}}} - {{\hbox{C}}_5}\prime - {\hbox{C}}_5^{\alpha } , which in turn increase the opportunity of conformations for binding to μ or δ receptors. Meanwhile, the connection of the two active tetrapeptide fragments of biphalin also results in the constrained side chain torsion angle χ² at one of the two residues Phe. This constrained side chain torsion angle not only significantly increases the δ receptor binding affinity but also makes most of the δ receptor binding conformations of biphalin bind to the δ receptor through the fragment containing the mentioned residue Phe.