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

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.

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.

Vibrational Spectroscopic Map, Vibrational Spectroscopy, and Intermolecular Interaction

Vibrational spectroscopy is an essential tool in chemical analyses, biological assays, and studies of functional materials. Over the past decade, various coherent nonlinear vibrational spectroscopic techniques have been developed and enabled researchers to study time-correlations of the fluctuating frequencies that are directly related to solute-solvent dynamics, dynamical changes in molecular conformations and local electrostatic environments, chemical and biochemical reactions, protein structural dynamics and functions, characteristic processes of functional materials, and so on. In order to gain incisive and quantitative information on the local electrostatic environment, molecular conformation, protein structure and interprotein contacts, ligand binding kinetics, and electric and optical properties of functional materials, a variety of vibrational probes have been developed and site-specifically incorporated into molecular, biological, and material systems for time-resolved vibrational spectroscopic investigation. However, still, an all-encompassing theory that describes the vibrational solvatochromism, electrochromism, and dynamic fluctuation of vibrational frequencies has not been completely established mainly due to the intrinsic complexity of intermolecular interactions in condensed phases. In particular, the amount of data obtained from the linear and nonlinear vibrational spectroscopic experiments has been rapidly increasing, but the lack of a quantitative method to interpret these measurements has been one major obstacle in broadening the applications of these methods. Among various theoretical models, one of the most successful approaches is a semiempirical model generally referred to as the vibrational spectroscopic map that is based on a rigorous theory of intermolecular interactions. Recently, genetic algorithm, neural network, and machine learning approaches have been applied to the development of vibrational solvatochromism theory. In this review, we provide comprehensive descriptions of the theoretical foundation and various examples showing its extraordinary successes in the interpretations of experimental observations. In addition, a brief introduction to a newly created repository Web site (http://frequencymap.org) for vibrational spectroscopic maps is presented. We anticipate that a combination of the vibrational frequency map approach and state-of-the-art multidimensional vibrational spectroscopy will be one of the most fruitful ways to study the structure and dynamics of chemical, biological, and functional molecular systems in the future.

Unperturbed Detection of the Dynamic Structure in the Hydrophobic Core of Trp-Cage via Two-Dimensional Infrared Spectroscopy

The tyrosine ring mode is an intrinsic non-perturbing site-specific infrared reporter for conformational dynamics within protein systems. This transition is influenced by direct and indirect interactions associated with the electron-donating ability and the hydrophobicity of the surrounding molecules. Utilizing an intrinsic tyrosine moiety, two-dimensional infrared spectra of Trp-cage, often called the "hydrogen atom" of protein folding, were measured in the folded and denatured states to uncover the dynamics of the hydrophobic core. The vibrational lifetimes and the correlation decays of the tyrosine ring mode showed significant changes upon both temperature and chemical denaturation of the Trp-cage miniprotein, indicating important structural features of the hydrophobic core and its dynamics. The observed Trp6-Tyr3 interactions are in good agreement with the prior studies of the folded state, but they reach beyond the static structure. These stacking interactions and orientations fluctuate on the picosecond time scale as measured through the spectral dephasing within a dehydrated environment.

Spectroscopic and computational study of melittin, cecropin A, and the hybrid peptide CM15

Antimicrobial peptides (AMPs), such as cecropin A from silk moth, are key components of the innate immune system. They are effective defensive weapons against invading pathogens, yet they do not target host eukaryotic cells. In contrast, peptide toxins, such as honeybee melittin, are nondiscriminating and target both eukaryotic and prokaryotic cells. An AMP-toxin hybrid peptide that is composed of cecropin A and melittin (CM15) improves upon the antimicrobial activity of cecropin A without displaying the nonspecific, hemolytic properties of melittin. Here we report fluorescence and UV resonance Raman spectra of melittin, cecropin A, and CM15 with the goal of elucidating peptide-membrane interactions that help guide specificity. We have probed the potency for membrane disruption, local environment and structure of the single tryptophan residue, backbone conformation near the peptide hinge, and amide backbone structure of the peptides in lipid environments that mimic eukaryotic and prokaryotic membranes. These experimental results suggest that melittin inserts deeply into the bilayer, whereas cecropin A remains localized to the lipid headgroup region. A surprising finding is that CM15 is a potent membrane-disruptor despite its largely unfolded conformation. A molecular dynamics analysis complements these data and demonstrates the ability of CM15 to associate favorably with membranes as an unfolded peptide. This combined experimental–computational study suggests that new models for peptide–membrane interactions should be considered.