Enkephalin: conformational analysis by means of empirical energy calculations

Insulin analogues with modifications in the beta-turn of the B-chain

The beta-turn formed by the amino acid residues 20-23 of the B-chain of insulin has been implicated as an important structural feature of the molecule. In other biologically active peptides, stabilization of beta-turns has resulted in increases in activity. We have synthesized three insulin analogues containing modifications which would be expected to increase the stability of the beta-turn. In two analogues, we have substituted alpha-aminoisobutyric acid (Aib) for the Glu residue normally present in position B21 or for the Arg residue normally present in position B22; in a third compound, we have replaced the Glu residue with its D-isomer. Biological evaluation of these compounds showed that [B21 Aib]insulin displays a potency ca. one-fourth that of natural insulin, while [B22 Aib]insulin is less than 10% as potent. In contrast, [B21 D-Glu]insulin is equipotent with natural insulin. We conclude that the beta-turn region of the insulin molecule normally possesses considerable flexibility, which may be necessary for it to assume a conformation commensurate with high biological activity. If this is the case, [B21 D-Glu]insulin may exhibit a stabilized geometry similar to that of natural insulin when bound to the insulin receptor.

A supermolecule study of the effect of hydration on the conformational behaviour of leucine-enkephalin

A theoretical conformational study was performed on leu-enkephalin in its zwitterionic form, both in vacuo and in the presence of a number, n, of up to 13 water molecules saturating its first hydration shell. The intramolecular energy of enkephalin as well as the intermolecular enkephalin-water and water-water interaction energies were computed with the SIBFA procedure (Sum of Interactions Between Fragments Ab initio computed), which uses additive ab initio multipole systematics and analytical formulas grounded on ab initio SCF computations. Energy minimizations were performed with a polyvalent minimizer, Merlin, with which three distinct derivative and three distinct nonderivative minimizers can be activated in a sequential fashion. Eight different candidate conformations of enkephalin were used as starting points. These conformations are either those found in distinct X-ray structures, or those proposed on the basis of theoretical computations by other authors. In the absence of hydration, they converged towards distinct folded energy-minima, the best four ones being separated by an energy gap of 8.7 kcal/mol. In marked contrast, with up to n = 13, the energetical separation between the six best conformers narrowed down to congruent to 4 kcal/mol. They can be characterized by: (a) either a direct or a water-mediated ammonium-carboxylate interaction; b) either a close proximity (as in morphine) or a large separation between the aromatic rings of Tyr and Phe (intercenter separations of congruent to 4.5 A and congruent to 10.5 A, respectively), with each of the four mutual combinations of (a) and (b) being represented.

Determining minimum energy conformations of polypeptides by dynamic programming

A combinatorial optimization approach is used for solving the multiple-minima problem when determining the low-energy conformations of short polypeptides. Each residue is represented by a finite number of discrete states corresponding to single residue local minima of the energy function. These precomputed values constitute a search table and define the conformational space for discrete minimization by a generalized dynamic programming algorithm that significantly limits the number of intermediate conformations to be generated during the search. Since dynamic programming involves stagewise decisions, it results in buildup-type procedures implemented in two different forms. The first procedure predicts a number of conformations by a completely discrete search and these are subsequently refined by local minimization. The second involves limited continuous local minimization within the combinatorial algorithm, generally restricted to two dihedral angles in a buildup step. Both procedures are tested on 17 short peptides previously studied by other global minimization methods but involving the same potential energy function. The discrete method is extremely fast, but proves to be successful only in 14 of the 17 test problems. The version with limited local minimization finds, however, conformations in all the 17 examples that are close to the ones previously presented in the literature or have lower energies. In addition, results are almost independent of the cutoff energy, the most important parameter governing the search. Although the limited local minimization increases the number of energy evaluations, the method still offers substantial advantages in speed.

Conformational features responsible for the binding of cyclic analogues of enkephalin to opioid receptors. II. Models of mu- and delta-receptor-bound structures for analogues containing Phe4

Models of mu- and delta-receptor-bound backbone conformations of enkephalin cyclic analogues containing Phe4 were determined by comparing geometrical similarity among the previously found low-energy backbone structures of [D-Cys2,Cys5]-enkephalinamide, [D-Cys2,D-Cys5]-enkephalinamide, [D-Pen2,L-Pen5]-enkephalin and [D-Pen2,D-Pen5]-enkephalin. The present mu-receptor-bound conformation resembles a beta-I bend in the peptide backbone centred on the Gly3-Phe4 region. Two slightly different models were found for the delta-receptor-bound conformation; both of them are more extended than the mu-receptor-bound conformation and include a gamma-turn (or a gamma-like turn) on the Gly3 residue. Energetically favourable rotamers of Tyr and Phe side chains were also determined for the mu- and delta-conformations. The present models of mu- and delta-conformations share geometrical similarity with the low-energy structures of Leu-enkephalin and the Tyr-D-Lys-Gly-Phe-analogue.

Transferred nuclear Overhauser effect analyses of membrane-bound enkephalin analogues by 1H nuclear magnetic resonance: correlation between activities and membrane-bound conformations

Leu-enkephalin, [D-Ala2]Leu-enkephalin, and [D-Ala2]Leu-enkephalinamide (agonists) and [L-Ala2]Leu-enkephalin (inactive analogue) bind to lipid bilayer consisting of phosphatidylcholine and phosphatidylserine. The conformations that these compounds assume, once bound to perdeuterated phospholipid bilayer, have been shown to be unique, as shown by the transferred nuclear Overhauser effect (TRNOE) of 1H NMR spectroscopy. In addition, their location in the bilayer was analyzed by TRNOE in the presence of spin-labeled phospholipids. These analyses showed a clear relationship between the activity and the peptide-membrane interaction. The three active peptides, when bound to membranes, adopt the same conformation, characterized by a type II' beta-turn around Gly3-Phe4 and a gamma-turn around Gly2 (or D-Ala2). The inactive analogue, [L-Ala2]Leu-enkephalin, displayed a completely different TRNOE pattern corresponding to a different conformation in the membrane-bound state. The tyrosine residue of the active compounds is not inserted into the interior of membrane, but it is inserted into the bilayer for the L-Ala2 analogue. According to these results, [L-Ala2]Leu-enkephalin may be explained to be inactive because the mode of binding to the membranes is different from that of active compounds.

Applications of simulated annealing to peptides

We report the application of a new conformation searching algorithm called simulated annealing to the location of the global minimum energy conformation of peptides. Simulated annealing is a Metropolis Monte Carlo approach to conformation generation in which both the energy and temperature dependence of the Boltzmann distribution guides the search for the global minimum. Both uphill and downhill moves are possible, which allows the molecule to escape from local minima. Applications to the 20 natural amino acid "dipeptide models" as well as to polyalanines up to Ala80 are very successful in finding the lowest energy conformation. A history file of the simulated annealing process allows reconstruction and examination of the random walk around conformation space. A separate program, Conf-Gen, reads the history file and extracts all low-energy conformations visited during the run.

The multiple-minima problem in the conformational analysis of polypeptides. III. An electrostatically driven Monte Carlo method: tests on enkephalin

The three-dimensional conformation of Met-enkephalin, corresponding to the lowest minimum of the empirical potential energy function ECEPP/2 (empirical conformational energy program for peptides), has been determined using a new algorithm, viz. the Electrostatically Driven Monte Carlo Method. This methodology assumes that a polypeptide or protein molecule is driven toward the native structure by the combined action of electrostatic interactions and stochastic conformational changes associated with thermal movements. These features are included in the algorithm that produces a Monte Carlo search in the conformational hyperspace of the polypeptide, using electrostatic predictions and a random sampling technique to locate low-energy conformations. In addition, we have incorporated an alternative mechanism that allows the structure to escape from some conformational regions representing metastable local energy minima and even from regions of the conformational space with great stability. In 33 test calculations on Met-enkephalin, starting from arbitrary or completely random conformations, the structure corresponding to the global energy minimum was found in all the cases analyzed, with a relatively small search of the conformational space. Some of these starting conformations were right or left-handed alpha-helices, characterized by good electrostatic interactions involving their backbone peptide dipoles; nevertheless, the procedure was able to convert such locally stable structures to the global-minimum conformation.

Synthesis of 17O isotope labeled Leu-enkephalin and 17O n.m.r

Oxygen-17 isotope was introduced into the alpha-carboxyl group of glycine, 1-phenylalanine, 1-leucine and 1-tyrosine by acid catalyzed exchange of 17O from H2O(17) or by acid hydrolysis of respective amino acid methyl esters in H2O(17). Quantitative enrichment of glycine was achieved by acid hydrolysis of amino acetonitrile in H2O(17). For alpha-amino protection in amino acids t-butoxycarbonyl (Boc) group was employed for 17O labeled enkephalin synthesis. Five analogues of Leu-enkephalins (I-V) labeled with 17O at different amino acid residues were synthesized by solid phase method. 17O n.m.r. spectra were measured at 24.4 and 67.8 MHz for Leu-enkephalins 17O labeled at Gly2 and Phe4 positions. A downfield shift was observed for 17O labeled Gly2 Leu-enkephalin upon heating. This shift is indicative of the rupture of intramolecular hydrogen bonds. The preliminary results confirm the hypothesis that an intramolecular hydrogen bond exists between the carbonyl group of Gly2 and NH group of Leu5.

On the biologically active structures of cholecystokinin, little gastrin, and enkephalin in the gastrointestinal system

The biologically active conformations of a series of four peptides [four cholecystokinin (CCK)-related peptides and enkephalin] in their interactions with gastrointestinal receptors have been deduced using conformational computational analysis. The two peptides that interact exclusively with peripheral-type CCK receptors are the heptapeptide COOH-terminal fragment from CCK (CCK-7) and the analogous sequence from cerulein (CER-7) in which threonine replaces the methionine proximal to the NH2 terminus. The two peptides that interact exclusively with the gastrin receptor in the stomach are the active COOH-terminal fragment of little gastrin and the COOH-terminal tetrapeptide sequence common to all of these peptides, CCK-4. We find that preferred conformations for the peripherally active peptides CCK-7 and CER-7 are principally beta-bends, whereas little gastrin and CCK-4 are fundamentally helical. In the class of lowest energy structures for both CCK-7 and CER-7, the aromatic rings of the tyrosine and phenylalanine lie close to one another whereas the tryptophan indole ring points in the opposite direction. This structure is superimposable on the structures of a set of rigid indolyl benzodiazepine derivatives that interact with complete specificity and high affinity with peripheral CCK receptors further suggesting that the computed beta-bends are the biologically active conformation. The biologically active conformation for CCK-4 and the little gastrin hexapeptide has also been deduced. By excluding conformations common to CCK-7 and CCK-4, which do not bond to each other's receptors, and then by selecting conformations in common to CCK-4 and the gastrin-related hexapeptide, which do bind to each other's receptors, we deduce that the biologically active conformation at the gastrin receptor is partly helical and one in which the indole of tryptophan and the aromatic ring of phenylalanine are close to one another while the methionine and aspartic acid side chains point in the opposite direction. These major differences in preferred structures between the common CCK-7/CER-7 peptides and the common CCK-4/little gastrin peptides explain the mutually exclusive activities of these two sets of peptides. We have observed that [Met]enkephalin strongly antagonizes the action of the naturally occurring peripherally active CCK-8 (CCK-7 with an NH2-terminal aspartic acid residue added). The computed lowest energy structures for this opiate peptide closely resemble key features of the computed CCK-7/CER-7 structure, further supporting the proposed structure.

Solvent dependent conformational statistics of enkephalin and angiotensin II

A screened electrostatic potential model of hydration, recently proposed, is used to simulate the statistical behaviour of enkephalin and angiotensin II in solution. Curves of the end-to-end distances as well as calculated values of n.m.r. coupling constants are given and compared with experiments and results obtained from other theoretical studies. The satisfactory agreement between calculated values and the corresponding experimental data indicates that the present theoretical approach is a simple and efficient way to take into account solvent effects on the conformation of biological molecules.

NOE data at 500 MHz reveal the proximity of phenyl and tyrosine rings in enkephalin

Met5-enkephalin-a pentapeptide (Tyr-Gly-Gly-Phe-Met)-can exist in two possible folded arrangements with a rigid two-hydrogen-bonded network. In one arrangement, a Gly 2-Gly 3 beta-bend is formed and in the other a Gly 3-Phe 4 beta-bend. The two conformations are distinguished by the spatial relation of Tyr 1 and Phe 4: in the Gly 2-Gly 3 beta-bend, Tyr 1 and Phe 4 can be brought close to each other while in the Gly 3-Phe 4 beta-bend they are far apart (greater than 5 A). We have utilized one-dimensional (1D) nuclear Overhauser effect (NOE) measurements between the ring protons of Tyr 1 and Phe 4 to determine their proximity. The NOE data clearly show that a pair protons, one each from Tyr 1 and Phe 4, are as close as 3.3 A while other inter-proton distances are beyond 4.5 A. Therefore, we propose the presence of a Gly 2-Gly 3 beta-bend (in which Tyr 1 and Phe 4 are spatially close) for Met5-enkephalin in solution. The structure of Met5-enkephalin in solution is very similar to the single crystal structure of Leu5-enkephalin and tends to explain the biological activity data of several modified enkephalins.

X-ray diffraction studies of enkephalins. Crystal structure of [(4'-bromo) Phe4,Leu5]enkephalin

In order to investigate the structure-activity relationship of [Leu5]- and [Met5]enkephalins, [(4'-bromo)Phe4, Leu5]-, [(4'-bromo)Phe4, Met5]- and [Met5] enkephalins were synthesized and crystallized. The crystal structure of [(4'-bromo) Phe4, Leu5]- enkephalin was determined by X-ray diffraction method using the heavy atom method and refined to R = 0.092 by the least-squares method. The molecule in this crystal took essentially the same type I' beta-turn conformation found in [Leu5]enkephalin [Smith & Griffin (1978) Science 199, 1214-1216). On the other hand, the preliminary three-dimensional Patterson analyses showed that the most probable conformations of [(4'-bromo)Phe4,Met5]- and [Met5]enkephalins are both the dimeric extended forms. Based on these insights, the biologically active conformation of enkephalin was discussed in relation to the mu- and delta-receptors.

Crystal structure of leucine-enkephalin

The crystal structure of leucine-enkephalin has been determined in a crystal form that has four independent enkephalin molecules and much water and dimethylformamide solvent in the asymmetric unit. All four enkephalins have extended peptide backbones with the Tyr, Phe and Leu side chains above and below the plane of the backbone. There is evidence that this extended conformation may provide an acceptable model for enkephalin binding to opiate mu-receptors.

Stabilization of beta-turn conformations in enkephalins. alpha-Aminoisobutyric acid analogs

Stereochemical constraints have been introduced into the enkephalin backbone by substituting alpha-aminoisobutyryl (Aib) residues at positions 2 and 3, instead of Gly. 1H n.m.r. studies of Tyr-Aib-Gly-Phe-Met-NH2, Tyr-Aib-Aib-Phe-Met-NH2 and Tyr-Gly-Aib-Phe-Met-NH2 demonstrate the occurrence of folded, intramolecularly hydrogen bonded structures in organic solvents. Similar conformations are also favoured in the corresponding t-butyloxycarbonyl protected tetrapeptides, which lack the Tyr residue. A beta-turn centred at positions 2 and 3 is proposed for the Aib2-Gly3 analog. In the Gly2-Aib3 analog, the beta-turn has Aib3-Phe4 as the corner residues. The Aib2-Aib3 analog adopts a consecutive beta-turn or 3(10) helical conformation. High in vivo biological activity is observed for the Aib2 and Aib2-Aib3 analogs, while the Aib3 peptide is significantly less active.

The design of biologically active polypeptides

On the one hand considerable advances are being made in the computer-aided calculation of polypeptide structure and behaviour, while on the other even very simple reasoning has lead to some apparent successes in the design and synthesis of biologically active peptides. While this looks extremely promising, there is a case for arguing that progress made in both these areas is to some degree illusory. Such criticisms as can be made, however, turn out to be highly constructive, suggesting a profitable and still workable approach to the problem of rational polypeptide design.

270-MHz 1H nuclear-magnetic-resonance study of met-enkephalin in solvent mixtures. Conformational transition from dimethylsulphoxide to water

1H spectra at 270 MHz of zwitterionic Met-enkephalin pentapeptide (Tyr-Gly-Gly-Phe-Met) in (C2H3)2SO/water mixtures are reported and discussed in terms of solvent-induced conformational transitions. The analysis of the chemical shifts, line widths, coupling constants and rotamer populations around chi 1 and chi 2 suggests that the conformational properties of Met-enkephalin in the two solvents are quite different. In aqueous solution, the preferred structure, characterized by the absence of intramolecularly hydrogen-bonded NH groups and head-to-tail interactions, very likely is an equilibrium of unfolded conformations with approximately equal energy. In (C2H3)2SO, the preferred structure is folded, with the Met-5 NH intramolecularly bonded and the Gly-3 NH protected from the solvent, while the Gly-2 and Phe-4 amide protons are solvent exposed. A conformational transition of Met-enkephalin from the intramolecularly bonded to the unbonded one takes place at about 40 mol-% water in (C2H3)2SO, involving the Met-5 NH proton and the Tyr-Gly-Gly fragment. The Phe-4 and Met-5 phi angles do not change appreciably, which suggest that an inversion at the Gly-3 residue of the folded form, responsible for the conformational transition, does not affect the C-terminal moiety. At about 70 mol-% water in (C2H3)2SO a change in the solvent mixture properties affects the chi 1 rotamer populations and the ring dynamics of the aromatic side chains. The line broadening of the Tyr-1 delta and epsilon proton resonances indicates a specific interaction of the N-terminal ring with the solvent.

Theoretical conformational analysis of enkephalin analogues related to fluorescence and n.m.r. measurements in aqueous solution

The properties related to non-radiative energy transfer of a number of enkephalin analogues with tryptophan substituted for phenylalanine in position 4 and n.m.r. 3JNH-C alpha H coupling constants of corresponding [Phe4]-enkephalin analogues are being derived from semi-empirical conformational energy. The molecules considered contain a glycyl, a D-alanyl or an L-alanyl as second residue; two of the compounds are N-methylated at position 4 or 5. The [Trp4]-enkephalin analogues and the corresponding [Phe4]-enkephalin analogues display nearly parallel affinities in the opiate receptor binding assay (Schiller et al. (1). The comparison of computed and experimental properties shows that an ensemble of conformers is a satisfactory representation of the state of these molecules in water.

Conformational energy calculations on enkephalins and enkephalin analogs. Classification of conformations to different configurational types

Conformational energy calculations were carried out on the peptide enkaphalins (ENK) and selected analogs to find those conformers of low energy. The analogs studied include [D-Ala2]Enk-NH2, [D-Ala2]Enk, [D-Met2, Pro5]Enk-NH2, [D-Ala2, D-Phe5]Enk, [D-Ala2, D-Leu5]Enk, [D-Ala2, (N-Me)Phe4, Met5] Enk-NH2 and [D-Ala2, (N-Me)Met5]Enk-NH2. When the low-energy conformers for all the analogs are compared, different allowed backbone conformations are found which orient the functional side-chains such that three classifications of structures appear. Each classification shows a unique configuration of side-chain positions in space even though different backbone conformations are found within each classification.

Reverse turns in peptides and proteins

The evidence that reverse turns frequently occur as structural components of proteins, as well as of linear and cyclic peptides, is overwhelming. This review summarizes and examines critically the experimental evidence derived from physical methods such as 1H and 13C nuclear magnetic resonance spectroscopy, spin-lattice relaxation time, circular dichroism, IR spectroscopy, and X-ray crystallography. Secondly, theoretical evidence obtained from energy calculations, which rely on empirical energy functions, and correlative methods, which rely on algorithms based on the frequency of occurrence of amino acids, is evaluated. In particular, those theoretical studies for which complementary physical studies have been completed are emphasized. Finally, examples of structure-function relationships involving reverse turns and their biological recognition are demonstrated.

Comparative theoretical conformational analysis of the insulin A6-A11 peptides

The insulins of pig, ox, horse and goat differ only by the cyclic A6-A11 peptide which is thought to be an antigenic determinant of the molecule. The structure of the four peptides is investigated by conformational energy calculations in order to verify whether a common backbone conformation can be found for all four species, the antigenic difference being consequently due only to the difference in the side chains exposed to the solvent, or whether each sequence gives rise to a preferential backbone conformation, which would lead to the conclusion that the antigenic difference is conveyed by a more pronounced difference in the molecular structure. From the results of the study on the isolated peptides, it appears that an energetically favourable backbone conformation common to all four species can be found. This conformation is compared with the structure deduced from the X-ray diffraction data available for pig insulin.

Intramolecularly hydrogen-bonded peptide conformations

Over the past few years the possible occurrence of intramolecularly hydrogen-bonded structures in linear and cyclic peptides has attracted increasing attention. In this review emphasis is given to solid-state studies, particularly by X-ray diffraction and infrared absorption techniques. Conformational energy calculations are also considered. The discussion is focused both on model peptides and biological activity polypeptide molecules. The tetrapeptide system (Formula: see text), examined allows one to discuss the extended C5 structure and the various folded conformations, namely the C7 (gamma-turn), C8, C10 (beta-turn), C11, and C13 conformations. The four latter forms may include cis peptide configurations. The oxy-analogs to the C7, C10, and C13 conformations and structures containing bifurcated hydrogen bonds are also discussed. The last sections describe intramolecularly hydrogen-bonded peptide structures involving: (1) a side-chain group, (2) the N-protecting group (in synthetic model compounds), and (3) a beta-amino acid.

Synthetic enkephalins. Addicting properties and conformational studies in solution

The addicting properties of [Leu5]enkephalin in mice are conserved in the LSer3 analogue and lost both in the L-Ser2 analogue and in all the L-Cha4 derivatives of the above peptides. Fluorescence measurements in water show the presence of hydrogen-bonded tyrosul OH groups in [Leu5]-enkephalin and in its LSer2 analogue. The Phe4/Cha replacements do not influence these equilibria, but they affect the near u.v. dichroism of the hydrogen bonded tyrosyl residues. In the peptide absorption region in water solution, only [Leu5]-enkephalin and its cyclohexylalnyl derivative show a positive dichroism towards high frequencies, which is maintained in 8 M ura. No clear relation is found between conformation(s) in solution and biological activity. A II' beta-turn, with residues in positions 2 and 3 at the corners, is suggested for the conformation of enkephalin bound to the receptors involved in the bioassay here used.