PubMed1988
To investigate the biologically active conformation of enkephalin, molecular-dynamics simulations were applied to [Met5]- and [D-Ala2,Met5]-enkephalins. The dynamic trajectory of monomeric extended [Met5]-enkephalin was analysed in terms of relative mobility between respective torsions of backbone chain. After 10 ps of the dynamics simulation, the conformational transition was converged into a stationary state among the beta-bend folded forms, where they are stabilized by several intramolecular hydrogen-bond formations. Similar conformational transition was also observed in the dynamics simulation of [D-Ala2,Met5]enkephalin, which is a more mu-receptor-specific peptide than [Met5]enkephalin. The geometrical correspondence between the monomeric enkephalin conformation in the stationary state and morphine molecule (a mu-specific rigid opiate) was surveyed by virtue of the triangular substructures generated by choosing three functional atoms in each molecule, and good resemblances were observed. On the other hand, the dynamics simulation of the antiparallel extended [Met5]enkephalin dimer showed a trajectory different from that of the monomeric one. Two intermolecular hydrogen bonds at Tyr1 (NH3+)...Met5(CO2-) end residues were held throughout the 100 ps simulation, the dimeric structure being consequently kept. The conformational transition of the backbone chains from the antiparallel extended form to the twisted one took place via an intermediate state. Many conformations revealed during the dynamics simulation showed that the relative orientations of each two Tyr1, Gly3, Phe4 and Met5 residues in the dimer are nearly related by a pseudo-C2-symmetry respectively, and both halves of the dimer structure could be further fitted to the monomeric folded enkephalin conformation. The monomeric and dimeric conformations of enkephalin at their stationary states are discussed in relation to the substrate-specificity for mu- and delta-opioid receptors.