The crystal structures of [Met5]enkephalin and a third form of [Leu5]enkephalin: observations of a novel pleated beta-sheet

Exploring Conformational Preferences of Leu-enkephalin Using the Conformational Search and Double-Hybrid DFT Energy Calculations

The conformational preferences of Leu-enkephalin (Leu-Enk) were explored by the conformational search and density functional theory (DFT) calculations. By a combination of low-energy conformers of each residue, the initial structures of the neutral Leu-Enk were generated and optimized using the ECEPP3 force field in the gas phase. These structures were reoptimized at the HF/3-21G(d) and M06-2X levels of theory with 6-31G(d) and 6-31+G(d) basis functions. We finally located the 139 structures with the relative energy <10 kcal mol^-1 in the gas phase, from which the structures of the corresponding zwitterionic Leu-Enk were generated and reoptimized at the M06-2X/6-31+G(d) level of theory using the implicit solvation model based on density (SMD) in water. The conformational preferences of Leu-Enk were analyzed using Gibbs free energies corrected by single-point energies calculated at the double-hybrid DSD-PBEP86-D3BJ/def2-TZVP level of theory in the gas phase and in water. The neutral Leu-Enk dominantly adopted a folded structure in the gas phase stabilized by three H-bonds with a βII'-bend-like motif at the Gly3-Phe4 sequence and a close contact between the side chains of Phe4 and Leu5. The zwitterionic Leu-Enk exhibited a folded structure in water stabilized by three H-bonds with double β-bends such as a βII' bend at the Gly2-Gly3 sequence and a βI bend at the Gly3-Phe4 sequence. The calculated ensemble-averaged distance between C_Gly2 ^α and C_Leu5 ^α of the zwitterionic Leu-Enk in water is consistent with the value estimated from the simulated annealing using the distance constraints derived from nuclear Overhauser effect spectroscopy (NOESY) spectra in water. Interestingly, the preferred conformations of the neutral and zwitterionic Leu-Enk are new folded structures not predicted by earlier computational studies. According to the refined model of the zwitterionic Leu-Enk bound to δ-opioid receptor (δOR), there were favorable interactions of the terminal charged groups of Leu-Enk with the side chains of charged residues of δOR as well as a favorable C_Aryl···H interaction of the Phe4 residue of Leu-Enk with Trp284 of δOR. Hence, these favorable interactions would induce the folded structure of the zwitterionic Leu-Enk with double β-bends isolated in water into the "bioactive conformation" like an extended structure when binding to δOR.

Co-Evolution of Opioid and Adrenergic Ligands and Receptors: Shared, Complementary Modules Explain Evolution of Functional Interactions and Suggest Novel Engineering Possibilities

Cross-talk between opioid and adrenergic receptors is well-characterized and involves second messenger systems, the formation of receptor heterodimers, and the presence of extracellular allosteric binding regions for the complementary ligand; however, the evolutionary origins of these interactions have not been investigated. We propose that opioid and adrenergic ligands and receptors co-evolved from a common set of modular precursors so that they share binding functions. We demonstrate the plausibility of this hypothesis through a review of experimental evidence for molecularly complementary modules and report unexpected homologies between the two receptor types. Briefly, opioids form homodimers also bind adrenergic compounds; opioids bind to conserved extracellular regions of adrenergic receptors while adrenergic compounds bind to conserved extracellular regions of opioid receptors; opioid-like modules appear in both sets of receptors within key ligand-binding regions. Transmembrane regions associated with homodimerization of each class of receptors are also highly conserved across receptor types and implicated in heterodimerization. This conservation of multiple functional modules suggests opioid–adrenergic ligand and receptor co-evolution and provides mechanisms for explaining the evolution of their crosstalk. These modules also suggest the structure of a primordial receptor, providing clues for engineering receptor functions.

Multiple, consecutive, fully-extended 2.0₅-helix peptide conformation

The peptide 2.0(5)-helix does exist. It has been experimentally authenticated both in the crystalline state (by X-ray diffraction) and in solution (by several spectroscopic techniques). It is the most common conformation for C(α)-tetrasubstituted α-amino acids with at least two atoms in each side chain, provided that cyclization on the C(α)-atom is absent. X-Ray diffraction has allowed a detailed description of its geometrical and three-dimensional (3D)-structural features. The infrared absorption and the nuclear magnetic resonance parameters characteristics of this multiple, consecutive, fully-extended structure have been described. Conformational energy calculations are in agreement with the experimental findings. As the contribution per amino acid residue to the length of this helix is the longest possible, its exploitation as a molecular spacer is quite promising. However, it is a rather fragile 3D-structure and particularly sensitive to solvent polarity. Interestingly, in such a case, it may reversibly convert to the much shorter 3(10)-helix, thus generating an attractive molecular spring.

Fine refinement of solid-state molecular structures of Leu- and Met-enkephalins by NMR crystallography

This paper presents a methodology that allows the fine refinement of the crystal and molecular structure for compounds for which the data deposited in the crystallographic bases are of poor quality. Such species belong to the group of samples with molecular disorder. In the Cambridge Crystallographic Data Center (CCDC), there are approximately 22,000 deposited structures with an R-factor over 10. The powerful methodology we present employs crystal data for Leu-enkephalin (two crystallographic forms) with R-factor values of 14.0 and 8.9 and for Met-enkephalin (one form) with an R-factor of 10.5. NMR crystallography was employed in testing the X-ray data and the quality of the structure refinement. The GIPAW (gauge invariant projector augmented wave) method was used to optimize the coordinates of the enkephalins and to compute NMR parameters. As we reveal, this complementary approach makes it possible to generate a reasonable set of new coordinates that better correlate to real samples. This methodology is general and can be employed in the study of each compound possessing magnetically active nuclei.

Preparation and evaluation at the delta opioid receptor of a series of linear leu-enkephalin analogues obtained by systematic replacement of the amides

Leu-enkephalin analogues, in which the amide bonds were sequentially and systematically replaced either by ester or N-methyl amide bonds, were prepared using classical organic chemistry as well as solid phase peptide synthesis (SPPS). The peptidomimetics were characterized using competition binding, ERK1/2 phosphorylation, receptor internalization, and contractility assays to evaluate their pharmacological profile over the delta opioid receptor (DOPr). The lipophilicity (LogD7.4) and plasma stability of the active analogues were also measured. Our results revealed that the last amide bond can be successfully replaced by either an ester or an N-methyl amide bond without significantly decreasing the biological activity of the corresponding analogues when compared to Leu-enkephalin. The peptidomimetics with an N-methyl amide function between residues Phe and Leu were found to be more lipophilic and more stable than Leu-enkephalin. Findings from the present study further revealed that the hydrogen-bond donor properties of the fourth amide of Leu-enkephalin are not important for its biological activity on DOPr. Our results show that the systematic replacement of amide bonds by isosteric functions represents an efficient way to design and synthesize novel peptide analogues with enhanced stability. Our findings further suggest that such a strategy can also be useful to study the biological roles of amide bonds.

MOLS sampling and its applications in structural biophysics

This review describes the MOLS method and its applications. This computational method has been developed in our laboratory primarily to explore the conformational space of small peptides and identify features of interest, particularly the minima, i.e., the low energy conformations. A systematic "brute-force" search through the vast conformational space for such features faces the insurmountable problem of combinatorial explosion, whilst other techniques, e.g., Monte Carlo searches, are somewhat limited in their region of exploration and may be considered inexhaustive. The MOLS method, on the other hand, uses a sampling technique commonly employed in experimental design theory to identify a small sample of the conformational space that nevertheless retains information about the entire space. The information is extracted using a technique that is a variant of the self-consistent mean field technique, which has been used to identify, for example, the optimal set of side-chain conformations in a protein. Applications of the MOLS method to understand peptide structure, predict the structures of loops in proteins, predict three-dimensional structures of small proteins, and arrive at the best conformation, orientation, and positions of a small molecule ligand in a protein receptor site have all yielded satisfactory results.

Exploring the Backbone of Enkephalins To Adjust Their Pharmacological Profile for the δ-Opioid Receptor

The role of each of the four amide bonds in Leu(5)-enkephalin was investigated by systematically and sequentially replacing each with its corresponding trans-alkene. Six Leu(5)-enkephalin analogs based on six dipeptide surrogates and two Met(5)-enkephalin analogs were synthesized and thoroughly tested using a δ-opioid receptor internalization assay, an ERK1/2 activation assay, and a competition binding assay to evaluate their biological properties. We observed that an E-alkene can efficiently replace the first amide bond of Leu(5)- and Met(5)-enkephalin without significantly affecting biological activity. By contrast, the second amide bond was found to be highly sensitive to the same modification, suggesting that it is involved in biologically essential intra- or intermolecular interactions. Finally, we observed that the affinity and activity of analogs containing an E-alkene at either the third or fourth position were partially reduced, indicating that these amide bonds are less important for these intra- or intermolecular interactions. Overall, our study demonstrates that the systematic and sequential replacement of amide bonds by E-alkene represents an efficient way to explore peptide backbones.

NCAD, a database integrating the intrinsic conformational preferences of non-coded amino acids

Peptides and proteins find an ever-increasing number of applications in the biomedical and materials engineering fields. The use of non-proteinogenic amino acids endowed with diverse physicochemical and structural features opens the possibility to design proteins and peptides with novel properties and functions. Moreover, non-proteinogenic residues are particularly useful to control the three-dimensional arrangement of peptidic chains, which is a crucial issue for most applications. However, information regarding such amino acids--also called non-coded, non-canonical, or non-standard--is usually scattered among publications specialized in quite diverse fields as well as in patents. Making all these data useful to the scientific community requires new tools and a framework for their assembly and coherent organization. We have successfully compiled, organized, and built a database (NCAD, Non-Coded Amino acids Database) containing information about the intrinsic conformational preferences of non-proteinogenic residues determined by quantum mechanical calculations, as well as bibliographic information about their synthesis, physical and spectroscopic characterization, conformational propensities established experimentally, and applications. The architecture of the database is presented in this work together with the first family of non-coded residues included, namely, alpha-tetrasubstituted alpha-amino acids. Furthermore, the NCAD usefulness is demonstrated through a test-case application example.

An estimate of the numbers and density of low-energy structures (or decoys) in the conformational landscape of proteins

BACKGROUND

The conformational energy landscape of a protein, as calculated by known potential energy functions, has several minima, and one of these corresponds to its native structure. It is however difficult to comprehensively estimate the actual numbers of low energy structures (or decoys), the relationships between them, and how the numbers scale with the size of the protein.

METHODOLOGY

We have developed an algorithm to rapidly and efficiently identify the low energy conformers of oligo peptides by using mutually orthogonal Latin squares to sample the potential energy hyper surface. Using this algorithm, and the ECEPP/3 potential function, we have made an exhaustive enumeration of the low-energy structures of peptides of different lengths, and have extrapolated these results to larger polypeptides.

CONCLUSIONS AND SIGNIFICANCE

We show that the number of native-like structures for a polypeptide is, in general, an exponential function of its sequence length. The density of these structures in conformational space remains more or less constant and all the increase appears to come from an expansion in the volume of the space. These results are consistent with earlier reports that were based on other models and techniques.

Exploring conformational space using a mean field technique with MOLS sampling

The computational identification of all the low energy structures of a peptide given only its sequence is not an easy task even for small peptides,due to the multiple-minima problem and combinatorial explosion. We have developed an algorithm, called the MOLS technique,that addresses this problem, and have applied it to a number of different aspects of the study of peptide and protein structure. Conformational studies of oligopeptides, including loop sequences in proteins have been carried out using this technique. In general the calculations identified all the folds determined by previous studies,and in addition picked up other energetically favorable structures. The method was also used to map the energy surface of the peptides. In another application, we have combined the MOLS technique, using it to generate a library of low energy structures of an oligopeptide, with a genetic algorithm to predict protein structures. The method has also been applied to explore the conformational space of loops in protein structures.Further, it has been applied to the problem of docking a ligand in its receptor site, with encouraging results.

Solid-state IR-LD spectroscopic and theoretical analysis of glycine-containing peptides and their hydrochlorides

As part of an investigation on the coordination ability of peptides, structural analyses of the solid di-, tri, and tetrapeptides glycyl-glycine (GG), glycyl-glycyl-glycine (GGG), glycyl-glycyl-glycyl-glycine (GGGG), and their protonated hydrochlorides glycyl-glycine.HCl (GGH), glycyl-glycyl-glycine.HCl (GGGH), and glycyl-glycyl-glycyl-glycine.HCl (GGGGH) have been carried out. The quantum chemical calculations (Hartree-Fock/6-31++G**) and linear-dichroic infrared (IR-LD) spectroscopy predict a near to linear structure of the pure ligands, but the experimental IR-LD data are in accordance with a cross-linked disposition of amide fragments in the protonated forms.

Conformational stability and three-dimensional model of the δ-opioid pharmacophore for the extended antiparallel dimer structure of Met-enkephalin in water

The conformational stability of the extended antiparallel dimer structure of Met-enkephalin in water was analyzed by examining the hydration structure of enkephalin using molecular dynamics simulations. The result shows that, despite of the hydrophicility of the terminal atoms in the pentapeptide, the main contributor for the stability of the dimer in water is the four intermolecular hydrogen bonds between the Gly(2) and Phe(4) groups. The three-dimensional model of the delta-opioid pharmacophore for this dimer structure was also established. Such a model was demonstrated to match the delta-opioid pharmacophore query derived from the non-peptides SIOM, TAN-67, and OMI perfectly. This result thus strongly supports the assumption that the dimer structure of Met-enkephalin is a possible delta-receptor binding conformation.

Solid state linear-dichroic infrared spectral and theoretical analysis of methionine-containing tripeptides

Stereo structural characterization of methionine-containing tripeptides glycine-glycine-methionine (Gly-Gly-Met) and methionine-glycine-glycine (Met-Gly-Gly) is carried out by means of linear-dichroic infrared spectral analysis of oriented as nematic liquid crystal suspension solid samples. Thus obtained data are compared with theoretical ones received by ab initio calculations (Hartree-Fock level of theory and 6-31G** basis set). An IR-spectral characterization of both tripeptides is included as well.

Investigation of interaction of Leu-enkephalin with lipid membranes

Enkephalins are peptides with morphine-like activity. To achieve their biological function, they must be transported from an aqueous phase to the lipid-rich environment of their membrane bound receptor proteins. In our study, zeta potential (ZP) method was used to detect the association of Leu-enkephalin and Leu-enkephalinamide with phospholipid liposomes constituted from egg-phosphatidylcholine (EPC), dioleoyl-phosphatidylethanolamine (DOPE), cholesterol (Chol), sphingomyelin (SM) as well as soybean phospholipid (SBPL). Transfer of the peptides over lipid membranes was examined by electrophysiology technique (ET) and fluorescence spectroscopy (FS), and further confirmed using 4-fluoro-7-nitrobenzofurazan (NBD-F) labeled Leu-enkephalin (NBD-F-enkephalin) with confocal laser scanning microscopy method (CLSM). Results of zeta potential showed that enkephalinamide associated with lipid membranes and gradually saturated on the membranes either hydrophobically or electrostatically or both. Data from electrophysiology technique indicated that Leu-enkephalin could cause transmembrane currents, suggesting the transfer of peptides across lipid membranes. Transfer examined by fluorescence spectroscopy implied that it could be separated into three steps, adsorption, transportation and desorption, which was afterward reaffirmed by confocal laser scanning microscopy. Transfer efficiencies of enkephalin across SBPL, EPC/DOPE, EPC/DOPE/SM, EPC/SM and EPC/Chol lipid bilayer membranes were evaluated with ET and CLSM experiments. Results showed that the addition of either sphingomyelin or cholesterol, or negatively charged lipid in lipid membrane composition could lower the transfer efficiency.

Structural and spectroscopic analysis of hydrogensquarates of glycine-containing tripeptides

The hydrogensquarates of glycine-containing tripeptides glycylglycylglycine (H-Gly-Gly-Gly-OH), glycylglycylmethionine (H-Gly-Gly-Met-OH), and methionylglycylglycine (H-Met-Gly-Gly-OH) have been characterized structurally. Quantum chemical ab initio calculations, solid-state linear-dichroic infrared (IR-LD) spectroscopy, 1H and 13C NMR data, ESI-MS, HPLC-MS/MS, TGV, and DSC methods were employed. The structures consist in a positively charged peptide moiety and a negative hydrogensquarate anion (HSq), stabilized by strong intermolecular hydrogen bonds.

Nearest-neighbor nonparametric method for estimating the configurational entropy of complex molecules

A method for estimating the configurational (i.e., non-kinetic) part of the entropy of internal motion in complex molecules is introduced that does not assume any particular parametric form for the underlying probability density function. It is based on the nearest-neighbor (NN) distances of the points of a sample of internal molecular coordinates obtained by a computer simulation of a given molecule. As the method does not make any assumptions about the underlying potential energy function, it accounts fully for any anharmonicity of internal molecular motion. It provides an asymptotically unbiased and consistent estimate of the configurational part of the entropy of the internal degrees of freedom of the molecule. The NN method is illustrated by estimating the configurational entropy of internal rotation of capsaicin and two stereoisomers of tartaric acid, and by providing a much closer upper bound on the configurational entropy of internal rotation of a pentapeptide molecule than that obtained by the standard quasi-harmonic method. As a measure of dependence between any two internal molecular coordinates, a general coefficient of association based on the information-theoretic quantity of mutual information is proposed. Using NN estimates of this measure, statistical clustering procedures can be employed to group the coordinates into clusters of manageable dimensions and characterized by minimal dependence between coordinates belonging to different clusters.

The role of crystallography in drug design

Structure and function are intimately related. Nowhere is this more important than the area of bioactive molecules. It has been shown that the enantioselectivity of an enzyme is directly related to its chirality. X-ray crystallography is the only method for determining the "absolute" configuration of a molecule and is the most comprehensive technique available to determine the structure of any molecule at atomic resolution. Results from crystallographic studies provide unambiguous, accurate, and reliable 3-dimensional structural parameters, which are prerequisites for rational drug design and structure-based functional studies.

Conformational studies on enkephalins using the MOLS technique

Conformational studies of two linear enkephalin molecules, Met-enkephalin and Leu-enkephalin, have been carried using the mutually orthogonal Latin squares (MOLS) technique with the ECEPP/3 force field. This technique was developed recently in our laboratory to perform an unbiased search of the conformational space of peptides and to locate low energy conformations. The present study identified all the folds predicted by other studies, and in addition picked up other energetically favorable structures. The results suggest that the peptide backbone exists as a mixture of folded and unfolded forms (approximately 50% each). The study also provides information on the distribution of the low energy conformations that we have classified on the basis of structural motifs, backbone hydrogen-bonding patterns, and root mean square deviations in atomic positions.

A multidimensional 1H NMR investigation of the conformation of methionine-enkephalin in fast-tumbling bicelles

Enkephalins are pentapeptides found in the central nervous system. It is believed that these neuropeptides interact with the nerve cell membrane to adopt a conformation suitable for their binding to an opiate receptor. In this work, we have determined the three-dimensional structure of methionine-enkephalin (Menk) in fast-tumbling bicelles using multidimensional (1)H NMR. Bicelles were selected as model membranes because both their bilayer organization and composition resemble those of natural biomembranes. The effect of the membrane composition on the peptide conformation was explored using both zwitterionic (PC bicelles) and negatively charged bicelles (Bic/PG). Pulsed field gradient experiments allowed the determination of the proportion of Menk bound to the model membranes. Approximately 60% of the water-soluble enkephalin was found to associate to the bicellar systems. Structure calculations from torsion angle and NOE-based distance constraints suggest the presence of both micro - and delta-selective conformers of Menk in each system and slightly different conformers in PC bicelles and Bic/PG. As opposed to previous studies of enkephalins in membrane mimetic systems, our results show that these opiate peptides could adopt several conformations in a membrane environment, which is consistent with the flexibility and poor selectivity of enkephalins.

Structural alterations of enkephalins in the presence of GM1 ganglioside micelles

The enkephalins are endogenous neurotransmitters and bind with high affinity at the delta-receptor. Gangliosides, the major glycans of nerve cells, known to interact both with receptors and ligands on the cell surface, have been implicated to modulate the actions of opioid receptors by allosteric regulation (Wu, G.; Lu, Z. H.; Wei, T. J.; Howells, R. D.; Christoffers, K.; Leeden R. W. Ann NY Acad Sci 1998, 845, 126-138). We have studied the interactions between enkephalins and monosialylated ganglioside GM1 using NMR spectroscopy and fluorescence. The structural models of enkephalins in the presence of GM1 micelles were generated using two-dimensional (1)H-ROESY experiments along with restrained molecular dynamics simulations. We report a conformational alteration of enkephalins in the presence of GM1 micelles.

The probability distribution of side-chain conformations in [Leu] and [Met]enkephalin determines the potency and selectivity to μ and δ opiate receptors

The structure of enkephalin, a small neuropeptide with five amino acids, has been simulated on computers using molecular dynamics. Such simulations exhibit a few stable conformations, which also have been identified experimentally. The simulations provide the possibility to perform cluster analysis in the space defined by potentially pharmacophoric measures such as dihedral angles, side-chain orientation, etc. By analyzing the statistics of the resulting clusters, the probability distribution of the side-chain conformations may be determined. These probabilities allow us to predict the selectivity of [Leu]enkephalin and [Met]enkephalin to the known mu- and delta-type opiate receptors to which they bind as agonists. Other plausible consequences of these probability distributions are discussed in relation to the way in which they may influence the dynamics of the synapse.

Enhanced sampling of the molecular potential energy surface using mutually orthogonal latin squares: application to peptide structures

The computational identification of the optimal three-dimensional fold of even a small peptide chain from its sequence, without reference to other known structures, is a complex problem. There have been several attempts at solving this by sampling the potential energy surface of the molecule in a systematic manner. Here we present a new method to carry out the sampling, and to identify low energy conformers of the molecule. The method uses mutually orthogonal Latin squares to select (of the order of) n(2) points from the multidimensional conformation space of size m(n), where n is the number of dimensions (i.e., the number of conformational variables), and m specifies the fineness of the search grid. The sampling is accomplished by first calculating the value of the potential energy function at each one of the selected points. This is followed by analysis of these values of the potential energy to obtain the optimal value for each of the n-variables separately. We show that the set of the n-optimal values obtained in this manner specifies a low energy conformation of the molecule. Repeated application of the method identifies other low energy structures. The computational complexity of this algorithm scales as the fourth power of the size of the molecule. We applied this method to several small peptides, such as the neuropeptide enkephalin, and could identify a set of low energy conformations for each. Many of the structures identified by this method have also been previously identified and characterized by experiment and theory. We also compared the best structures obtained for the tripeptide (Ala)(3) by the present method, with those obtained by an exhaustive grid search, and showed that the algorithm is successful in identifying all the low energy conformers of this molecule.

Assemblies of Alzheimer's peptides A beta 25-35 and A beta 31-35: reverse-turn conformation and side-chain interactions revealed by X-ray diffraction

Alzheimer's beta amyloid protein (A beta) is a 39 to 43 amino acid peptide that is a major component in the neuritic plaques of Alzheimer's disease (AD). The assemblies constituted from residues 25-35 (A beta 25-35), which is a sequence homologous to the tachykinin or neurokinin class of neuropeptides, are neurotoxic. We used X-ray diffraction and electron microscopy to investigate the structure of the assemblies formed by A beta 25-35 peptides and of various length sequences therein, and of tachykinin-like analogues. Most solubilized peptides after subsequent drying produced diffraction patterns characteristic of beta-sheet structure. Moreover, the peptides A beta 31-35 (Ile-Ile-Gly-Leu-Met) and tachykinin analogue A beta(Phe(31))31-35 (Phe-Ile-Gly-Leu-Met) gave powder diffraction patterns to 2.8A Bragg spacing. The observed reflections were indexed by an orthogonal unit cell having dimensions of a=9.36 A, b=15.83 A, and c=20.10 A for the native A beta 31-35 peptide, and a=9.46 A, b=16.22 A, and c=11.06 A for the peptide having the Ile31Phe substitution. The initial model was a beta strand where the hydrogen bonding, chain, and intersheet directions were placed along the a, b, and c axes. An atomic model was fit to the electron density distribution, and subsequent refinement resulted in R factors of 0.27 and 0.26, respectively. Both peptides showed a reverse turn at Gly33 which results in intramolecular hydrogen bonding between the antiparallel chains. Based on previous reports that antagonists for the tachykinin substance P require a reverse turn, and that A beta is cytotoxic when it is oligomeric or fibrillar, we propose that the tachykinin-like A beta 31-35 domain is a turn exposed at the A beta oligomer surface where it could interact with the ligand-binding site of the tachykinin G-protein-coupled receptor.

The structure of an endomorphin analogue incorporating 1-aminocyclohexane-1-carboxlylic acid for proline is similar to the beta-turn of Leu-enkephalin

Endomorphin (EM2, Tyr-Pro-Phe-Phe-NH(2)) can assume various conformations related to cis/trans-rotamers of the amide linkage of Tyr-Pro. To control isomerization, restricted or flexible components have been introduced at the Pro position. We focused on [Chx(2)]EM2, an EM2 analogue substituting 1-aminocyclohexane-1-carboxlylic acid (Chx) for Pro. X-ray diffraction analysis revealed that [Chx(2)]EM2 is folded into the trans-form of Tyr-Chx. The manner of folding resembled that seen in D-TIPP, an EM analogue incorporating tetrahydroisoquinoline carboxylic acid, as well as the beta-turn of Leu-enkephalin. Selectivity for the opioid mu-receptor was fairly well conserved by [Chx(2)]EM, suggesting that the folded form is important for mu-selectivity.

A conformational analysis of leucine enkephalin as a function of pH

The conformations of Leu enkephalin in aqueous solution have been investigated as a function of pH using molecular dynamics simulations. The simulations suggest the peptide backbone exists as a mixture of folded and unfolded forms (approximately 50% each) at neutral pH, but is always unfolded at low or high pH. The folded form at neutral pH possesses a 2 --> 5 hydrogen bond and a close head to tail separation. No significant intramolecular hydrogen bonding of the carbonyl oxygens was observed in either the folded or unfolded forms of the peptide. Analysis of the Gly carbonyl oxygens and terminal groups indicated that, while the conformational population distribution of Leu enkephalin did vary noticeably as a function of pH, their hydration was essentially independent of pH and in agreement with the available NMR data. Further study indicated that the unfolded state of the peptide was not random in nature and consisted of one major unfolded backbone arrangement stabilized by a persistent hydrophobic interaction between the side chains of Tyr and Leu.

Cross-checking of nanoelectrospray ionization mass spectrometry and computer simulation for the evaluation of the interaction strength of non-covalently bound enkephalins in solution

Nanoelectrospray ionization mass spectrometry (nanoESI-MS) and computer simulation were applied to the characterization of non-covalent interactions of [Leu5]-enkephalin (LE) and its optical isomers, [D-Tyr1, Leu5]-enkephalin (Y-LE), [D-Phe4, Leu5]-enkephalin (F-LE) and [D-Tyr1, D-Phe4, Leu5]-enkephalin (YF-LE). The dimer formation tendencies of the optical isomers of LE were evaluated by nanoESI-MS using quadruply deuterated LE (H2N-Tyr-(2,2-d2)Gly-(2,2-d2)Gly-Phe-Leu-COOH, d4-LE) as an internal standard. The relative interaction strengths of the optical isomers of LE were estimated to be Y-LE < F-LE < LE < YF-LE. Geometry optimization calculations were performed for interactions in vacuo and in water using a semi-empirical SCF method (PM3). The initial coordinate of the dimer structure of LE was taken from that obtained from single-crystalline x-ray diffraction analysis. Estimates of the interaction strengths of the dimer complexes were based on the heats of formation of a dimer complex (Hd) and the corresponding monomers (Hm) using the equation DeltaH = Hd - 2Hm. The values of DeltaH obtained from the calculations for interactions in water decreased in the order Y-LE > F-LE > LE > YF-LE. Since the smaller values of DeltaH correspond to stronger interactions between peptides, the results from computer simulations were qualitatively consistent with those obtained from the nanoESI experiments. The possibility of cross-checking these independent techniques was demonstrated using medium-sized molecules of biological importance. The agreement of the results from the two techniques suggested that nanoESI experiments, at least qualitatively, reflected the relative interaction strengths of non-covalently bound enkephalins in aqueous solution.

Conformational analysis of [Met5]-enkephalin: solvation and ionization considerations

[Met5]-Enkephalin has the sequence Tyr-Gly-Gly-Phe-Met. Only the extended conformation of the peptide has been observed by X-ray crystallography. Nuclear magnetic resonance spectroscopy supports the presence of a turn at Gly 3 and Phe 4 in dimethyl sulfoxide. In this study, the peptide conformational states and thermodynamic properties are understood in terms of ionization state and solvent environment. In the calculation, final conformations obtained from multiple independent Monte Carlo simulated annealing conformational searches are starting points for molecular dynamics simulations. In an aqueous environment given by the use of solvation free energy and the zwitterionic state, dominant structural motifs computed are G-P Type II' bend, G-G Type II' bend, and G-G Type I' bend motifs, in order of increasing free energy. In the calculation of the peptide with neutral N- and C-termini and solvation free energy, the extended conformer dominates (by at least a factor of 2.5), and the conformation of another low free energy conformer superimposes well on the pharmacophoric groups of morphine. Neutralization of charge and solvation induce and stabilize the extended conformation, respectively. A mechanism of inter-conversion between the extended conformer and three bent conformers is supported by phi/psi-scatter plots, and by the conformer relative free energies. An estimate of the entropy change of receptor unbinding is 8.3 cal K-1 mol-1, which gives rise to a -2.5 kcal/mol entropy contribution to the free energy of unbinding at 25 degrees C. The conformational analysis methodology described here should be useful in studies on short peptides and flexible protein surface loops that have important biological implications.

Structural studies of opioid peptides: a review of recent progress in x-ray diffraction studies

The solid state structures of many opioid peptide agonists have been elucidated by x-ray diffraction analysis. Recently, the first structure of an opioid peptide antagonist has been determined. Theoretically, linear peptides can have many different backbone conformations, yet early x-ray studies (1983-1987) on enkephalin and its analogues showed only two different backbone conformations: extended and single beta-bend. In 1989 enkephalin was observed in a third conformation, a double beta-bend. Since that time diffraction studies have been completed on the rationally designed linear opioid peptide agonists DTLET (Tyr-D-Thr-Gly-Phe-Leu-Thr) and DADLE (D-Ala2,D-Leu5-enkephalin) as well as on several cyclic enkephalin analogues including DPDPE (Tyr-[D-Pen-Gly-Phe-D-Pen]) and JOM-13 (Tyr-[D-Cys-Phe-D-Pen]). The most recent review of the x-ray studies on this class of compounds was written in 1988. This paper will update that review to include the results of studies completed since that time.

Characteristic molecular packing in the crystal structure of tert-butoxycarbonyl-L-phenylalanyl-L-methionine methyl ester

The molecular conformation and association of the peptide Boc-L-Phe-L-Met-OMe have been studied in the solid state by X-ray diffraction. The peptide crystallizes in the orthorhombic system, space group P2(1)2(1)2(1), with cell parameters of a = 9.821(2), b = 25.394(6), c = 28.714(8) A, V = 7161(3) A3. The structure has been solved by direct methods and refined to a final R of 0.079 for 5464 independent reflections with Fo > or = sigma(Fo). The crystal consists of three independent molecular conformations per asymmetric unit. Respective peptide backbones adopt an extended conformation with the side-chains of Phe and Met residues being arranged below and above the backbone chains. Contrary to the sheet structure most frequently observed in the crystal packing of the extended peptide conformations, three independent molecules lie spirally along the c-axis and form a pin-wheel-like crystal packing. The sheet structures formed by two of three independent molecules are almost at right angles to the backbone of the remaining molecule. This molecular packing mode would provide a possible interaction model between the intersecting beta-sheet structure and single-strand structure of polypeptide.

Observation of a sterically unfavorable side-chain conformation in a leucyl residue: crystal and molecular structure of L-leucyl-L-leucine.DMSO solvate

The crystal structure of a dipeptide L-leucyl-L-leucine (C12H24N2O3) has been determined. The crystals are monoclinic, space group P2(1), with a = 5.434(4) A, b = 15.712(7) A, c = 11.275(2) A, beta = 100.41(1) degrees, and Z = 2. The crystals contain one molecule of dimethyl sulfoxide (DMSO) as solvent of crystallization for each dipeptide molecule. The structure has been solved by direct methods and refined to a final R index of 0.059 for 920 reflections (sin theta/lambda < or = 0.60 A-1) with I > or = 2 sigma (I). The trans peptide unit shows substantial degree of non-planarity (delta omega = 14 degrees). The peptide backbone adopts an extended conformation with torsion angles of psi 1 = 138(1) degrees, omega 1 = 166(1) degrees, phi 2 = -149.3(7) degrees, psi 21 = 164.2(7) degrees, and psi 22 = -15(1) degrees. For the first leucyl residue, the side-chain conformation is specified by the torsion angles 1 chi 1 = 176.7(7) degrees, 1 chi 21 = 62(1) degrees, 1 chi 22 = -177.4(8) degrees; the second leucyl residue adopts a sterically unfavorable conformation with 2 chi 1 = 61(1) degrees, 2 chi 21 = 97(1) degrees, and 2 chi 22 = -151(1) degrees. The packing involves head-to-tail interaction of peptide molecules and segregation of polar and nonpolar regions. The DMSO molecule is strongly hydrogen bonded to the terminal NH3+ group.

Crystal structure of deltakephalin: a delta-selective opioid peptide with a novel beta-bend-like conformation

The solid-state structure of deltakephalin (Tyr-DThr-Gly-Phe-Leu-Thr) has been determined by single-crystal X-ray diffraction. Deltakephalin (DTLET) is a synthetic opioid peptide which differs from enkephalin in that a D-Thr has been substituted for Gly2 and a sixth residue, L-Thr, has been added. Clear colorless plates obtained using vapor diffusion and macro-seeding crystallization techniques were monoclinic; space group C2 with a = 27.389(5), b = 9.205(2), c = 16.788(2) A, beta = 98.87(2) degrees and V = 4181.4(14) A3. The asymmetric unit contained one molecule of DTLET and six molecules of water, giving a calculated density of 1.28 g cm-3. The crystal structure revealed that DTLET has a pseudo type I' beta-bend which is stabilized by an intramolecular side-chain to backbone hydrogen bond. This is the first reported observation of a pseudo beta-bend conformation in a solid-state structure of an enkephalin analog.

Conserved and novel structural characteristics of enantiomorphic Leu-enkephalin. X-ray crystal analysis of Leu-enkephalin enantiomer, L-Tyr-Gly-Gly-L-Phe-L-Leu and D-Tyr-Gly-Gly-D-Phe-D-Leu

The crystal of the Leu-enkephalin racemate (L-Tyr-Gly-Gly-L-Phe-L-Leu and D-Tyr-Gly-Gly-D-Phe-D-Leu) was obtained as a centrosymmetric space group. Crystal data: C28H37N5O7 x 1.5H2O, Mw = 582.6, triclinic, space group P1, a = 11.176(3), b = 16.115(3), c = 10.204(4) A, alpha = 92.41(3), beta = 104.86(2), gamma = 85.35(2)degrees, V = 1770(1)A3, Z = 2; F(000) = 640, mu(CuK alpha) = 6.50 cm-1, D chi = 1.081 g cm-3. The structure was determined by X-ray diffraction. The conformation of the Leu-enkephalin racemate was classified into the extended form which has been often observed in natural enkephalin. The symmetry-related molecules were connected by hydrogen bonds and arranged in an antiparallel fashion. The molecular packing showed a sheet structure similar to that of natural enkephalin.

Investigation of conformational equilibrium of polypeptides by internal coordinate stochastic dynamics. Met5-enkephalin

The equilibrium population of different conformational states of a polypeptide can in principle be obtained by a very long molecular dynamics simulation. The method of internal coordinate molecular dynamics earlier developed in this laboratory (A.K. Mazur and R.A. Abagyan J. Biomol. Struct. Dyn. 6,833 (1989)) allows one to use time steps much larger than usual for computing molecular trajectories. It is shown here that the sampling of the conformational space can be additionally enhanced by adding a random component to the set of forces applied to atoms. We describe the algorithms by which the random force is introduced and also a special method which excludes the fast rotation of polar hydrogens from equations of motion but keeps them movable. As a result the task stated in the title becomes realistic. Internal coordinate stochastic dynamics is applied for scanning the conformational space of the pentapeptide Met5-enkephalin which is a common test example widely used in theoretical studies. A large number of conformational transitions is observed during the 20 ns simulation starting from the global energy minimum thus allowing us to arrive at a nearly Boltzmann distribution of populations of conformational states. A few states are found which are distinguished by high apparent configurational entropy which turn out to correspond well to experimentally observed conformations of enkephalins.

Necessary conditions for avoiding incorrect polypeptide folds in conformational search by energy minimization

Low energy conformations have been generated for melittin, pancreatic polypeptide, and ribonuclease S-peptide, both in the vicinity of x-ray structures by energy refinement and by an unconstrained search over the entire conformational space. Since the structural polymorphism of these medium-sized peptides in crystal and solution is moderate, comparing the calculated conformation to x-ray and nmr data provides information on local and global behavior of potential functions. Local analysis includes standardization calculations, which show that models with standard geometry can approximate good resolution x-ray data with less than 0.5 A rms deviation (RMSD). However, the atomic coordinates are shifted up to 2 A RMSD by local energy minimization, and thus 2 A is generally the smallest RMSD value one can target in a conformational search using the same energy evaluation models. The unconstrained search was performed by a buildup-type method based on dynamic programming. To accelerate the generation of structures in the conformational search, we used the ECEPP potential, defined in terms of standard polypeptide geometry. A number of low energy conformations were further refined by relaxing the assumption of standard bond lengths and bond angles through the use of the CHARMM potential, and the hydrophobic folding energies of Eisenberg and McLachlan were calculated. Each conformation is described in terms of the RMSD from the native, hydrogen-bonding structure, solvent-accessible surface area, and the ratio of surfaces corresponding to nonpolar and polar residues. The unconstrained search finds conformations that are different from the native, sometimes substantially, and in addition, have lower conformational energies than the native. The origin of deviations is different for each of the three peptides, but in all examples the refined x-ray structures have lower energies than the calculated incorrect folds when (1) the assumption of standard bond lengths and bond angles is relaxed; (2) a small and constant effective dielectric permittivity (epsilon < 10) is used; and (3) the hydrophobic folding energy is incorporated into the potential.