Conformational features of linear and cyclic enkephalins. A computational study

A theoretical conformational study using the CICADA program package (J. Mol. Struct. (Theochem), 337 (1995) 17) was performed for two linear enkephalins, Leu-enkephalin and Met-enkephalin, and two cyclic analogues, DLFE and DPDPE. The conformational flexibilities of whole molecules and selected torsions were calculated.

The low energy conformers obtained were compared with structures obtained by spectroscopic methods. The mutual space positions of key elements for receptor recognition were analyzed. Conformations were clustered using RMS deviation computed for selected atoms. The different conformational behavior of aromatic rings in cyclic analogues of enkephalins was observed. While aromatic rings of cyclic analogues exhibit different conformational behavior, the linear enkephalins show similar behavior in these key parts.

Hydrogen bonds predicted by spectroscopic measurements were confirmed by our calculations. Also very specific conformational features, like concerted conformational movement, were analyzed.     

Conformational free energy surface of the linear DPDPE peptide: cost of pre-organization for disulfide bond formation

We have calculated the free energy differences between four conformers of the linear form of the opioid pentapeptide DPDPE in aqueous solution. The conformers are Cyc, representing the structure adopted by the linear peptide prior to disulfide bond formation, β _C and β _E , two slightly different β-turns previously identified in unconstrained molecular dynamics simulations, and Ext, an extended structure. Our simulations indicate that β _E is the most stable of the studied conformers of linear DPDPE in aqueous solution, with β _C , Cyc and Ext having free energies higher by 2.3, 6.3, and 28.2 kcal/mol, respectively. The free energy differences of 4.0 kcal/mol between β _C and Cyc, and 6.3 kcal/mol between β _E and Cyc, reflect the cost of pre-organizing the linear peptide into a conformation conducive for disulfide bond formation. Such a conformational change is a pre-requisite for the chemical reaction of S–S bond formation to proceed. The relatively low population of the cyclic-like structure agrees qualitatively with observed lower potency and different receptor specificity of the linear form relative to the cyclic peptide, and with previous unconstrained simulation results. Free energy component analysis indicates that the moderate stability difference of 4.0–6.3 kcal/mol between the β-turns and the cyclic-like structure results from cancellation of two large opposing effects. In accord with intuition, the relaxed β-turns have conformational strain 43–45 kcal/mol lower than the Cyc structure. However, the cyclic-like conformer interacts with water about 39 kcal/mol strongly than the open β-turns. Our simulations are the first application of the recently developed multidimensional conformational free energy thermodynamic integration (CFTI) protocol to a solvated system, with fast convergence of the free energy obtained by fixing all flexible dihedrals. Additionally, the availability of the CFTI multidimensional free energy gradient leads to a new decomposition scheme, giving the contribution of each fixed dihedral to the overall free energy change and providing additional insight into the microscopic mechanisms of the studied processes. Received: 20 April 1998 / Accepted: 9 September 1998 / Published online: 7 December 1998