Orientation and Conformational Preference of Leucine-Enkephalin at the Surface of a Hydrated Dimyristoylphosphatidylcholine Bilayer: NMR and MD Simulation

On the deformability of additivated phosphatidylcholine liposomes: Molecular dynamic regimes and membrane elasticity

Liposomes with enhanced elasticity have been proven to increase the efficiency of drug transport across the skin. The understanding of the background physicochemical processes driving the liposome viscoelastic properties is an essential feature for the design of effective formulations involving different lipids and additive molecules. In this work we use field-cycled nuclear magnetic resonance relaxometry to analyze both the mechanical properties of liposome membranes, and their relationship with the involved molecular dynamics. Different liposomal formulations were considered. We show a correlation between the molecular dynamical regime and mesoscopic physical parameters that define the expected deformability of the vesicles. Results strongly suggest that the purity of the used lipids may influence the elastic properties of the membranes in an appreciable way. Common features in the behaviour of the involved dynamic variables were identified by comparing formulations with surfactants of similar molecular weight.

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.

Tailoring the nuclear Overhauser effect for the study of small and medium-sized molecules by solvent viscosity manipulation

The nuclear Overhauser effect (NOE) is a consequence of cross-relaxation between nuclear spins mediated by dipolar coupling. Its sensitivity to internuclear distances has made it an increasingly important tool for the determination of through-space atom proximity relationships within molecules of sizes ranging from the smallest systems to large biopolymers. With the support of sophisticated FT-NMR techniques, the NOE plays an essential role in structure elucidation, conformational and dynamic investigations in liquid-state NMR. The efficiency of magnetization transfer by the NOE depends on the molecular rotational correlation time, whose value depends on solution viscosity. The magnitude of the NOE between ¹H nuclei varies from +50% when molecular tumbling is fast to -100% when it is slow, the latter case corresponding to the spin diffusion limit. In an intermediate tumbling regime, the NOE may be vanishingly small. Increasing the viscosity of the solution increases the motional correlation time, and as a result, otherwise unobservable NOEs may be revealed and brought close to the spin diffusion limit. The goal of this review is to report the resolution of structural problems that benefited from the manipulation of the negative NOE by means of viscous solvents, including examples of molecular structure determination, conformation elucidation and mixture analysis (the ViscY method).

Biophysical approaches for exploring lipopeptide-lipid interactions

In recent years, lipopeptides (LPs) have attracted a lot of attention in the pharmaceutical industry due to their broad-spectrum of antimicrobial activity against a variety of pathogens and their unique mode of action. This class of compounds has enormous potential for application as an alternative to conventional antibiotics and for pest control. Understanding how LPs work from a structural and biophysical standpoint through investigating their interaction with cell membranes is crucial for the rational design of these biomolecules. Various analytical techniques have been developed for studying intramolecular interactions with high resolution. However, these tools have been barely exploited in lipopeptide-lipid interactions studies. These biophysical approaches would give precise insight on these interactions. Here, we reviewed these state-of-the-art analytical techniques. Knowledge at this level is indispensable for understanding LPs activity and particularly their potential specificity, which is relevant information for safe application. Additionally, the principle of each analytical technique is presented and the information acquired is discussed. The key challenges, such as the selection of the membrane model are also been briefly reviewed.

•

A brief overview of topics to understand the generalities of lipopeptide (LP) science.

•

Main analytical techniques used to reveal the interaction and the distorting effect of LP on artificial membranes.

•

Guidelines for selecting of the most adequate membrane models for the given analytical technique.

Selective effect of cell membrane on synaptic neurotransmission

Atomistic molecular dynamics simulations were performed with 13 non-peptidic neurotransmitters (NTs) in three different membrane environments. The results provide compelling evidence that NTs are divided into membrane-binding and membrane-nonbinding molecules. NTs adhere to the postsynaptic membrane surface whenever the ligand-binding sites of their synaptic receptors are buried in the lipid bilayer. In contrast, NTs that have extracellular ligand-binding sites do not have a similar tendency to adhere to the membrane surface. This finding is a seemingly simple yet important addition to the paradigm of neurotransmission, essentially dividing it into membrane-independent and membrane-dependent mechanisms. Moreover, the simulations also indicate that the lipid composition especially in terms of charged lipids can affect the membrane partitioning of NTs. The revised paradigm, highlighting the importance of cell membrane and specific lipids for neurotransmission, should to be of interest to neuroscientists, drug industry and the general public alike.

Supramolecular structure of membrane-associated polypeptides by combining solid-state NMR and molecular dynamics simulations

Elemental biological functions such as molecular signal transduction are determined by the dynamic interplay between polypeptides and the membrane environment. Determining such supramolecular arrangements poses a significant challenge for classical structural biology methods. We introduce an iterative approach that combines magic-angle spinning solid-state NMR spectroscopy and atomistic molecular dynamics simulations for the determination of the structure and topology of membrane-bound systems with a resolution and level of accuracy difficult to obtain by either method alone. Our study focuses on the Shaker B ball peptide that is representative for rapid N-type inactivating domains of voltage-gated K(+) channels, associated with negatively charged lipid bilayers.

Determining the orientation and localization of membrane-bound peptides

Many naturally occurring bioactive peptides bind to biological membranes. Studying and elucidating the mode of interaction is often an essential step to understand their molecular and biological functions. To obtain the complete orientation and immersion depth of such compounds in the membrane or a membrane-mimetic system, a number of methods are available, which are separated in this review into four main classes: solution NMR, solid-state NMR, EPR and other methods. Solution NMR methods include the Nuclear Overhauser Effect (NOE) between peptide and membrane signals, residual dipolar couplings and the use of paramagnetic probes, either within the membrane-mimetic or in the solvent. The vast array of solid state NMR methods to study membrane-bound peptide orientation and localization includes the anisotropic chemical shift, PISA wheels, dipolar waves, the GALA, MAOS and REDOR methods and again the use of paramagnetic additives on relaxation rates. Paramagnetic additives, with their effect on spectral linewidths, have also been used in EPR spectroscopy. Additionally, the orientation of a peptide within a membrane can be obtained by the anisotropic hyperfine tensor of a rigidly attached nitroxide label. Besides these magnetic resonance techniques a series of other methods to probe the orientation of peptides in membranes has been developed, consisting of fluorescence-, infrared- and oriented circular dichroism spectroscopy, colorimetry, interface-sensitive X-ray and neutron scattering and Quartz crystal microbalance.

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.

Molecular dynamics simulations of a lipovitellin-derived amphiphilic beta-sheet homologous to apoB-100 beta-sheets at a hydrophobic decane-water interface

Lipovitellin, an egg-yolk lipoprotein, transports lipids in a pocket surrounded by amphiphilic beta-sheets. Its X-ray structure provides possibilities to study interactions between lipophilic beta-sheets and lipids at the atomic level. Here, we studied a 67-residue-long amphiphilic beta-sheet of lipovitellin previously suggested a suitable working model for studies of the lipid-binding behaviour of amphiphilic beta-sheet regions in apolipoprotein B-100 (apoB-100). We performed four molecular dynamics simulations with different starting configurations to define characteristics of the amphiphilic beta-sheet model at a decane-water interface. In each simulation the model beta-sheet bound keenly to the decane layer via its hydrophobic surface. The structural profiles showed unchanged secondary structure of the beta-sheet during the attachment. Also, aromatic side chains, especially tryptophans and tyrosines, mediated the attachment to the hydrophobic layer and influenced the orientation of the decane molecules that are in contact with the beta-sheet. In conclusion, the present simulations reveal high affinity of a lipovitellin-derived amphiphilic beta-sheet to a hydrophobic decane layer. They lay thereby the basis for further studies of the interaction between amphiphilic beta-sheets and lipids in complex molecular systems, like LDL particles, in which the large apoB-100 is the main protein component.

Study of the interaction of GFG tripeptide with cesium perfluorooctanoate micelles by means of NMR spectroscopy and MD simulations

The interaction of glycyl-phenylalanyl-glycine (GFG) with bilayers formed by cesium perfluorooctanoate (CsPFO) in water was investigated in the isotropic phase by means of 1H NMR and molecular dynamics (MD) simulations. Details on the preferential location of the different residues of GFG were obtained from selective variations of chemical shift with peptide concentration and of line width in the presence of the paramagnetic ion Mn2+. The analysis of 1H NMR spectra recorded at different concentrations and temperatures allowed the association constant and the enthalpy change upon binding to be evaluated. MD simulations highlighted the hydrogen bonds formed between the different GFG functional groups and the micelle. Both NMR and MD gave indications of high affinity of GFG with the micelle, with the N-terminal residue anchoring on the surface via hydrogen bonds with the micelle COO(-) groups.

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.

Structure, topology, and tilt of cell-signaling peptides containing nuclear localization sequences in membrane bilayers determined by solid-state NMR and molecular dynamics simulation studies

Cell-signaling peptides have been extensively used to transport functional molecules across the plasma membrane into living cells. These peptides consist of a hydrophobic sequence and a cationic nuclear localization sequence (NLS). It has been assumed that the hydrophobic region penetrates the hydrophobic lipid bilayer and delivers the NLS inside the cell. To better understand the transport mechanism of these peptides, in this study, we investigated the structure, orientation, tilt of the peptide relative to the bilayer normal, and the membrane interaction of two cell-signaling peptides, SA and SKP. Results from CD and solid-state NMR experiments combined with molecular dynamics simulations suggest that the hydrophobic region is helical and has a transmembrane orientation with the helical axis tilted away from the bilayer normal. The influence of the hydrophobic mismatch, between the hydrophobic length of the peptide and the hydrophobic thickness of the bilayer, on the tilt angle of the peptides was investigated using thicker POPC and thinner DMPC bilayers. NMR experiments showed that the hydrophobic domain of each peptide has a tilt angle of 15 +/- 3 degrees in POPC, whereas in DMPC, 25 +/- 3 degree and 30 +/- 3 degree tilts were observed for SA and SKP peptides, respectively. These results are in good agreement with molecular dynamics simulations, which predict a tilt angle of 13.3 degrees (SA in POPC), 16.4 degrees (SKP in POPC), 22.3 degrees (SA in DMPC), and 31.7 degrees (SKP in DMPC). These results and simulations on the hydrophobic fragment of SA or SKP suggest that the tilt of helices increases with a decrease in bilayer thickness without changing the phase, order, and structure of the lipid bilayers.