Structural alterations of enkephalins in the presence of GM1 ganglioside micelles

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 CGly2 α and CLeu5 α 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 CAryl···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.

Conformational heterogeneity of a leucine enkephalin analogue in aqueous solution and sodium dodecyl sulfate micelles: comparison of time-resolved FRET and molecular dynamics simulations

We have undertaken time-resolved Förster resonance energy transfer (FRET) and molecular dynamics simulations to analyze conformations and conformational heterogeneity of an analogue of leucine enkephalin in solution and in the presence of sodium dodecyl sulfate (SDS) micelles. Enkephalins are opioid pentapeptides that interact with opioid receptors in the central nervous system. We used time-correlated single-photon counting to detect energy transfer between the N-terminal tyrosine and a tryptophan residue substituted for phenylalanine at the 4 position. FRET from Tyr to Trp was measured over a temperature range from 5 to 55 degrees C in aqueous solution. By taking into account Tyr rotamer interconversion rates measured previously, we determined average distances between Tyr and Trp for the two populated rotameric conformations of Tyr. Molecular dynamics simulations (100 ns) support this analysis and indicate extensive conformational heterogeneity. The simulations also predict that the FRET orientational factor is correlated with the Tyr-Trp separation. Failure to account for the correlation between orientation and distance results in errors that appear to be largely offset in the leucine enkephalin analogue (YGGWL) by a weighting bias inherent in the R(-6) dependence of the energy-transfer rate. The Tyr lifetimes decrease upon titration of the peptides with SDS, indicating formation of compact conformations of the peptide in the micelle environment. This result is consistent with the conjecture that the lipid environment may induce formation of bioactive conformations of the peptide.

Evidence for effect of GM1 on opioid peptide conformation: NMR study on leucine enkephalin in ganglioside-containing isotropic phospholipid bicelles

Enkephalins are endogenous neuropeptides that have opioid-like activities and compete with morphines for the receptor binding. The binding of these neuropeptides to membrane appears crucial since enkephalins interact with the nerve cell membranes to achieve bioactive conformations that fit onto multiple receptor sites (micro, delta, and kappa). Using NMR spectroscopy, we have determined the solution structure of the small opiate pentapeptide leucine enkephalin in the presence of isotropic phospholipid bicelles: phosphocholine bicelles (DMPC:CHAPS 1:4) and phosphocholine bicelles doped with ganglioside GM1 (DMPC:CHAPS:GM1 1:4:0.3). Bicelles containing GM1 were found to interact strongly with leucine enkephalin, whereas a somewhat weaker interaction was observed in the case of bicelles without GM1. Structure calculation from torsion angles, chemical shifts, and NOE-based distance constraints explored that the peptide could flexibly switch between several mu- and delta-selective conformations in both the bicelles though micro-selective conformations turned out to be geometrically preferred in each bicellar system. A detailed analysis of the structures presented supports the variance over the singly associated conformation of enkephalin in nerve cell membranes.

Structural transition of a 15 amino acid residue peptide induced by GM1

The ganglioside GM1-binding peptide, p3, with a sequence of VWRLLAPPFSNRLLP, displayed a clear structural alteration depending on the presence or absence of GM1 micelles. The three-dimensional structures of the p3 peptide in the free and GM1 bound states were analyzed using two-dimensional NMR spectroscopic experiments with distance-restrained simulated annealing calculations. The NMR experiments for the p3 peptide alone indicated that the peptide has two conformers derived from the exchange of cis and trans forms at Pro(7)-Pro(8). Further study with theoretical modeling revealed that the p3 peptide has a curb conformation without regular secondary structure. On the other hand, the NMR studies for the p3 peptide with the GM1 micelles elucidated a trans conformer and gave a structure stabilized by hydrophobic interactions of beta- and helical turns. Based on these structural investigations, tryptophan, a core residue of the hydrophobic cluster, might be an essential residue for the recognition of the GM1 saccharides. The dynamic transition of the p3 peptide may play an important role in the function of GM1 as a multiple receptor as in the traditional pathway of the infection by cholera toxin.

A pulsed field gradient NMR diffusion investigation of enkephalin peptide-sodium dodecyl sulfate micelle association

Pulsed field gradient NMR (PFG-NMR) diffusion experiments were used to investigate the binding of leucine and methionine enkephalin peptides to anionic sodium dodecyl sulfate (SDS) micelles. The study was undertaken to investigate the mechanism of interaction between enkephalin peptides and SDS micelles and to determine if NMR-derived association constants, K(eq), can predict the elution order in electrokinetic chromatography (EKC). In EKC, peptides are separated on the basis of their interactions with micelles. The Leu-enkephalin peptide-micelle association constant increased from 130 +/- 8 to 1459 +/- 57 and 1744 +/- 64 M(-1), respectively, when an Arg or Lys was added to the C-terminus. The association constant of Leu-enkephalinamide was approximately equal to that of Leu-enkephalin-Arg. Substitution of Phe4 with a Trp or Gly2 with an Ala in the Leu-enkephalin peptides also increased the micelle binding affinity. These results confirm that the interaction of Leu-enkephalin peptides with SDS micelles is largely electrostatic and that the non-polar amino acid side chains interact with the hydrophobic micelle core. The peptide-micelle association constants for the cationic Met-enkephalin peptides were also greater than their zwitterionic counterparts. For example, the Met-enkephalin K(eq) value was 162 +/- 9 M(-1), while the association constants for Met-enkephalin-Arg, Met-enkephalin-Lys, and Met-enkephalinamide were, respectively, 674 +/- 31, 426 +/- 23, and 453 +/- 27 M(-1). In both Met-enkephalin and Met-enkephalinamide, replacing Gly2 with an Ala did not significantly increase the association constant. These results confirm that with the Met-enkephalin peptides, there was little or no interaction of the amino acid side chains with the micelle core. For both the Leu-enkephalin and Met-enkephalin peptides, the association constants were consistent with EKC results, in that the peptides with smaller K(eq) values were found to elute before those with larger association constants.

Interaction and structural study of kinin peptide bradykinin and ganglioside monosialylated 1 micelle

Partitioning of small proteins and peptides from the aqueous to membrane phase is often coupled with folding. In this work we examine the binding and folding of the kinin peptide, bradykinin (BK), in the presence of the ganglioside monosialylated 1 (GM1) micelle. Using two-dimensional NMR techniques, we have shown that at low concentration, GM1 micelle is able to induce a turn conformation to BK. A pulsed-field gradient diffusion NMR study indicated that the peptide partitions into the GM1 micelle with a DeltaG(part) of -3.14 +/- 0.03 kcal/mol. A saturation transfer difference (STD) NMR study indicated that the binding is mostly through hydrophobic residues.

Anomalous behavior of water around sodium dodecyl sulphate micelles

The dynamics, structural properties, and energetics of hydration water around a sodium dodecyl sulphate micelle have been investigated using molecular dynamics simulation. A clear revelation of the slow dynamics of the hydration water has been made by separate measurements of the rotational and translational properties. Calculated diffusion coefficients fall within the range of experimentally observed quantities. The water-micelle head group (MHG) hydrogen bond is more stable (by an amount approximately 7.0 kcal/mol) compared to the water-water hydrogen bond. The difference in stability of the water monomers forming different numbers of hydrogen bonds (n=0,1,2) with the MHG has clearly been shown from the analyses of their rotational relaxation, residence times, as well as the energy of interaction with different components of the system. The singly hydrogen-bonded water species is the most abundant and stable. The entropy plays the key role in controlling the relative abundance of the different species.

Conformational alteration of bradykinin in presence of GM1 micelle

We report here the interaction of bradykinin with ganglioside GM1 by circular dichroism, steady-state fluorescence, and one-dimensional 1H NMR spectroscopy. Circular dichroism spectroscopy is indicative of a turn formation of bradykinin backbone in the presence of GM1 micelle. The fluorescence quenching efficiencies of iodide and acrylamide are substantially reduced, indicating a shielding of phenylalanine residue of bradykinin from aqueous environment. Significant line broadening of NMR resonances of bradykinin, suggestive of motional restriction, is observed.