FT-IR spectroscopic study on conformational equilibria of [Leu]5-enkephalin in DMSO and 2H2O solutions

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.

Ultrafast vibrational dynamics of the tyrosine ring mode and its application to enkephalin insertion into phospholipid membranes as probed by two-dimensional infrared spectroscopy

Enkephalins are small opioid peptides whose binding conformations are catalyzed by phospholipid membranes. Binding to opioid receptors is determined by the orientation of tyrosine and phenylalanine side chains. In this work, we investigate the effects of different charged phospholipid headgroups on the insertion of the tyrosine side chain into a lipid bilayer using a combination of 2D IR spectroscopy, anharmonic DFT calculations, and third order response function modeling. The insertion is probed by using the ∼1515 cm^-1 tyrosine ring breathing mode, which we found exhibits rich vibrational dynamics on the picosecond timescale. These dynamics include rapid intramolecular vibrational energy redistribution (IVR), where some of the energy ends up in a dark state that shows up as an anharmonically shifted combination band. The waiting-time dependent 2D IR spectra also show an unusual line shape distortion that affects the extraction of the frequency-frequency correlation function (FFCF), which is the dynamic observable of interest that reflects the tyrosine side chain's insertion into the lipid bilayer. We proposed three models to account for this distortion: a hot-state exchange model, a local environment dependent IVR model, and a coherence transfer model. A qualitative analysis of these models suggests that the local environment dependent IVR rate best explains the line shape distortion, while the coherence transfer model best reproduced the effects on the FFCF. Even with these complex dynamics, we found that the tyrosine ring mode's FFCF is qualitatively correlated with the degree of insertion expected from the different phospholipid headgroups.

Interactions of tyrosine in Leu-enkephalin at a membrane-water interface: an ultrafast two-dimensional infrared study combined with density functional calculations and molecular dynamics simulations

The interactions of neuropeptides and membranes play an important role in peptide hormone function. Our current understanding of peptide-membrane interactions remains limited due to the paucity of experimental techniques capable of probing such interactions. In this work, we study the nature of opioid peptide-membrane interactions using ultrafast two-dimensional infrared (2D IR) spectroscopy. The high temporal resolution of 2D IR is particularly suited for studying highly flexible opioid peptides. We investigate the location of the tyrosine (Tyr) side chain of leucine-enkephalin (Lenk) in lipid bilayer membranes by measuring spectral diffusion of the phenolic ring vibrational mode in three different systems: Lenk in lipid bilayer membranes (bicelles), Lenk in deuterated water, and p-cresol in deuterated water. Frequency-frequency correlation functions obtained from waiting-time-dependent 2D IR spectra reveal an ultrafast decaying component with an approximately 1 ps time constant that is common for all three systems. On the basis of density functional theory calculations and molecular dynamics simulations, this spectral diffusion component is attributed to hydrogen-bond dynamics of the phenolic hydroxyl group interacting with bulk water. Unlike p-cresol in water, both Lenk systems exhibit static spectral inhomogeneity, which can be attributed to conformational distributions of Lenk that do not interconvert within 4 ps. Our results suggest that the Tyr side chain of Lenk in bicelles is located at the water-abundant region at the membrane-water interface and not embedded into the hydrophobic core.

Reaction of naphthalene-2,3-dicarboxaldehyde with enkephalins for LC-fluorescence and LC-MS analysis: conformational studies by molecular modeling and H/D exchange mass spectrometry

A new labeling method compatible with both laser-induced fluorescence (LIF) and MS detection for enkephalins, which uses naphthalene-2,3-dicarboxaldehyde (NDA) and a new nucleophilic agent (N,N-dimethylaminoethanethiol) is described. When the derivative is separated via reverse phase HPLC and detected via MS, two different peaks with similar exact mass but different fluorescence and fragmentation properties are obtained. To interpret these results, molecular modeling and H/D exchange mass spectrometry studies were investigated to test the hypothesis that the peak obtained by LC/LIF/MS analysis depends on the site of protonation of the labeled enkephalins. The peptides labeled with NDA and N,N-dimethylaminoethanethiol were separated on a reverse phase C18 column with a gradient of aqueous 0.1% formic acid and acetonitrile. In mass spectrometry, two peaks are observed with the same exact mass for each molecule while only one peak is detected using fluorescence. Tandem mass spectrometry experiments of ion m/z 809.5 were performed on each chromatographic peak; the first peak (which is not observed by LIF detection) gives a fragment corresponding to the loss of the aminothiol side chain while no fragmentation is observed on the second peak, which was detected by fluorescence. The hypothesis is that each peak represents the labeled enkephalin with different sites of protonation. According to this hypothesis, three fundamental conformations that were closed to the unlabeled leucine-enkephalin were obtained by molecular modeling: a beta-turn like conformation with two hydrogen bonds, a 3(10)-helix with an H bond, and finally, the extended form without any intramolecular interactions. H/D exchange mass spectrometry experiments with D(2)O and d(2-)formic acid as eluent was used to determine which conformation is involved in each peak.