Molecular dynamics simulations of opioid peptide analogs containing multiple conformational restrictions

Molecular dynamics simulations were performed on the potent and slightly mu-receptor selective cyclic dermorphin analog H-Tyr-D-Orn-Phe-Glu-NH2 as well as on analogs containing a conformationally restricted phenylalanine derivative in place of Phe in the 3 position of the peptide sequence. Peptides studied included the potent and highly mu-selective analogs H-Tyr-D-Orn-Aic-Glu-NH2 (Aic = 2-aminoindan-2-carboxylic acid), H-Tyr-D-Orn-Atc-Glu-NH2 (Atc = 2-aminotetralin-2-carboxylic acid) and H-Tyr-D-Orn-D-Atc-Glu-NH2, and the weakly active analog H-Tyr-D-Orn-Tic-Glu-NH2 (Tic = tetrahydroisoquinoline-3-carboxylic acid). Four different starting conformations were chosen for each peptide, and after equilibration each simulation was allowed to proceed for 100 picoseconds at 600 degrees K. The 14-membered ring structures in the Phe-, Aic-, L- and D-Atc-containing analogs showed moderate structural flexibility, while the peptide ring in the Tic-containing analog was more rigid. As theoretically predicted, the phi 3 and psi 3 angles of the Aic-, L- and D-Atc-containing analogs were limited to values of either about +50 degrees or -50 degrees during almost the entire period of the simulations. In the Tic-containing analog the phi 3 and psi 3 angles were 0 degrees and 90 degrees, respectively, and did not change for the entire duration of the simulation. The side chains of the constrained amino acids showed limited movement, but transitions between the allowed conformations did occur on the time scale of the simulations.(ABSTRACT TRUNCATED AT 250 WORDS)