Molecular dynamics and (2)H-NMR study of the influence of an amphiphilic peptide on membrane order and dynamics

Best of Two Worlds? How MD Simulations of Amphiphilic Helical Peptides in Membranes Can Complement Data from Oriented Solid-State NMR

The membrane alignment of helical amphiphilic peptides in oriented phospholipid bilayers can be obtained as ensemble and time averages from solid state ²H NMR by fitting the quadrupolar splittings to ideal α-helices. At the same time, molecular dynamics (MD) simulations can provide atomistic insight into peptide-membrane systems. Here, we evaluate the potential of MD simulations to complement the experimental NMR data that is available on three exemplary systems: the natural antimicrobial peptide PGLa and the two designer-made peptides MSI-103 and KIA14, whose sequences were derived from PGLa. Each peptide was simulated for 1 μs in a DMPC lipid bilayer. We calculated from the MD simulations the local angles which define the side chain geometry with respect to the peptide helix. The peptide orientation was then calculated (i) directly from the simulation, (ii) from back-calculated MD-derived NMR splittings, and (iii) from experimental ²H NMR splittings. Our findings are that (1) the membrane orientation and secondary structure of the peptides found in the NMR analysis are generally well reproduced by the simulations; (2) the geometry of the side chains with respect to the helix backbone can deviate significantly from the ideal structure depending on the specific residue, but on average all side chains have the same orientation; and (3) for all of our peptides, the azimuthal rotation angle found from the MD-derived splittings is about 15° smaller than the experimental value.

Interaction of cisplatin and two potential antitumoral platinum(ii) complexes with a model lipid membrane: a combined NMR and MD study

Multinuclear NMR and MD calculations highlighted the different interactions of cisplatin and two potential antitumoral Pt(ii) complexes with a model membrane.

Variation of the lateral mobility of transmembrane peptides with hydrophobic mismatch

A hydrophobic mismatch between protein length and membrane thickness can lead to a modification of protein conformation, function, and oligomerization. To study the role of hydrophobic mismatch, we have measured the change in mobility of transmembrane peptides possessing a hydrophobic helix of various length d(pi) in lipid membranes of giant vesicles. We also used a model system where the hydrophobic thickness of the bilayers, h, can be tuned at will. We precisely measured the diffusion coefficient of the embedded peptides and gained access to the apparent size of diffusing objects. For bilayers thinner than d(pi), the diffusion coefficient decreases, and the derived characteristic sizes are larger than the peptide radii. Previous studies suggest that peptides accommodate by tilting. This scenario was confirmed by ATR-FTIR spectroscopy. As the membrane thickness increases, the value of the diffusion coefficient increases to reach a maximum at h approximately = d(pi). We show that this variation in diffusion coefficient is consistent with a decrease in peptide tilt. To do so, we have derived a relation between the diffusion coefficient and the tilt angle, and we used this relation to derive the peptide tilt from our diffusion measurements. As the membrane thickness increases, the peptides raise (i.e., their tilt is reduced) and reach an upright position and a maximal mobility for h approximately = d(pi). Using accessibility measurements, we show that when the membrane becomes too thick, the peptide polar heads sink into the interfacial region. Surprisingly, this "pinching" behavior does not hinder the lateral diffusion of the transmembrane peptides. Ultimately, a break in the peptide transmembrane anchorage is observed and is revealed by a "jump" in the D values.

Structural and dynamic studies of the γ-M4 trans-membrane domain of the nicotinic acetylcholine receptor

A structural characterization of a synthetic peptide corresponding to the fourth transmembrane domain (M4-TMD) of the gamma-subunit of the nicotinic acetylcholine receptor from Torpedo californica has been undertaken. Solid-state NMR and CD spectroscopy studies indicate that upon reconstitution into lipid vesicles or magnetically aligned lipid bilayers, the synthetic M4-TMD adopts a linear alpha-helical conformation with the helix aligned within 15 degrees of the membrane normal. Furthermore, analysis of the motional averaging of anisotropic interactions present in the solid-state NMR spectra of the reconstituted peptide, indicate that the dynamics of the peptide within the bilayer are highly sensitive to the phase adopted by the lipid bilayer, providing an insight into how the interaction of lipids with this domain may play a important role in the modulation of this receptor by its lipid environment.

Protein–lipid interactions studied with designed transmembrane peptides: role of hydrophobic matching and interfacial anchoring (Review)

Biological membranes are characterized by a heterogeneous composition, which is not only manifested in the wide variety of their components, but also in aspects like the lateral organization, topology, and conformation of proteins and lipids. In bringing about the correct membrane structure, protein-lipid interactions can be expected to play a prominent role. The extent of hydrophobic matching between transmembrane protein segments and lipids potentially constitutes a versatile director of membrane organization, because a tendency to avoid hydrophobic mismatch could result in compensating adaptations such as tilt of the transmembrane segment or segregation into distinct domains. Also, interfacial interactions between lipid headgroups and the aromatic and charged residues that typically flank transmembrane domains may act as an organizing element. In this review, we discuss the numerous model studies that have systematically explored the influence of hydrophobic matching and interfacial anchoring on membrane structure. Designed peptides consisting of a polyleucine or polyleucine/alanine hydrophobic stretch, which is flanked on both sides by tryptophan or lysine residues, reflect the general layout of transmembrane protein segments. It is shown for phosphatidylcholine bilayers and for other model membranes that these peptides adapt a transmembrane topology without extensive peptide or lipid adaptations under conditions of hydrophobic matching, but that significant rearrangements can result from hydrophobic mismatch. Moreover, these effects depend on the nature of the flanking residues, implying a modulation of the mismatch response by interfacial interactions of the flanking residues. The implications of these model studies for the organization of biomembranes are discussed in the context of recent experiments with more complex systems.

Transmembrane Helix Packing of ErbB/Neu Receptor in Membrane Environment: A Molecular Dynamics Study

Dimerization or oligomerization of the ErbB/Neu receptors are necessary but not sufficient for initiation of receptor signaling. The two intracellular domains must be properly oriented for the juxtaposition of the kinase domains allowing trans-phosphorylation. This suggests that the transmembrane (TM) domain acts as a guide for defining the proper orientation of the intracellular domains. Two structural models, with the two helices either in left-handed or in right-handed coiling have been proposed as the TM domain structure of the active receptor. Because experimental data do not distinguish clearly helix-helix packing, molecular dynamics (MD) simulations are used to investigate the energetic factors that drive Neu TM-TM interactions of the wild and the oncogenic receptor (Val664/Glu mutation) in DMPC or in POPC environments. MD results indicate that helix-lipid interactions in the bilayer core are extremely similar in the two environments and raise the role of the juxtamembrane residues in helix insertion and helix-helix packing. The TM domain shows a greater propensity to adopt a left-handed structure in DMPC, with helices in optimal position for strong inter-helical Hbonds induced by the Glu mutation. In POPC, the right-handed structure is preferentially formed with the participation of water in inter-helical Hbonds. The two structural arrangements of the Neu(TM) helices both with GG4 residue motif in close contact at the interface are permissible in the membrane environment. According to the hypothesis of a monomer-dimer equilibrium of the proteins it is likely that the bilayer imposes structural constraints that favor dimerization-competent structure responsible of the proper topology necessary for receptor activation.

Biophysical studies of the interactions between 14-mer and 21-mer model amphipathic peptides and membranes: Insights on their modes of action

We have investigated the interactions between synthetic amphipathic peptides and zwitterionic model membranes. Peptides with 14 and 21 amino acids composed of leucines and phenylalanines modified by the addition of crown ethers have been synthesized. The 14-mer and 21-mer peptides both possess a helical amphipathic structure as revealed by circular dichroism. To shed light on their mechanism of membrane interaction, different complementary biophysical techniques have been used such as circular dichroism, fluorescence, membrane conductivity measurement and NMR spectroscopy. Results obtained by these different techniques show that the 14-mer peptide is a membrane perturbator that facilitate the leakage of species such as calcein and Na ions, while the 21-mer peptide acts as an ion channel. (31)P solid-state NMR experiments on multilamellar vesicles reveal that the dynamics and/or orientation of the polar headgroups are greatly affected by the presence of the peptides. Similar results have also been obtained in mechanically oriented DLPC and DMPC bilayers where different acyl chain lengths seem to play a role in the interaction. On the other hand, (2)H NMR experiments on multilamellar vesicles demonstrate that the acyl chain order is affected differently by the two peptides. Based on these studies, mechanisms of action are proposed for the 14-mer and 21-mer peptides with zwitterionic membranes.

Insights on the interactions of synthetic amphipathic peptides with model membranes as revealed by 31P and 2H solid-state NMR and infrared spectroscopies

We studied the interaction between synthetic amphipathic peptides and model membranes by solid-state NMR and infrared spectroscopies. Peptides with 14 and 21 amino acids composed of leucines and phenylalanines modified by the addition of crown ethers were synthesized. The 14-mer and 21-mer peptides both possess a helical amphipathic structure. To shed light on their membrane interaction, (31)P and (2)H solid-state NMR experiments were performed on both peptides in interaction with dimyristoylphosphatidylcholine vesicles in the absence and presence of cholesterol, dimyristoylphosphatidylglycerol vesicles, and oriented bicelles. (31)P NMR experiments on multilamellar vesicles reveal that the dynamics and/or orientation of the polar headgroups are weakly yet markedly affected by the presence of the peptides, whereas (31)P NMR experiments on bicelles indicate no significant changes in the morphology and orientation of the bicelles. On the other hand, (2)H NMR experiments on vesicles reveal that the acyl chain order is affected differently depending on the membrane lipidic composition and on the peptide hydrophobic length. Finally, infrared spectroscopy was used to study the interfacial region of the bilayer. Based on these studies, mechanisms of membrane perturbation are proposed for the 14-mer and 21-mer peptides in interaction with model membranes depending on the bilayer composition and peptide length.

Transmembrane Peptides from Tyrosine Kinase Receptor. Mutation-related Behavior in a Lipid Bilayer Investigated by Molecular Dynamics Simulations

Polar mutations in transmembrane alpha helices may alter the structural details of the hydrophobic sequences and control intermolecular contacts. We have performed molecular dynamics simulations on the transmembrane domain of the proto-oncogenic and the oncogenic forms of the Neu receptor in a fluid DMPC bilayer to test whether the Glu mutation which replaces the Val residue at position 664 may alter the helical structure and its insertion in the membrane. The simulations show that the wild and the mutant forms of the transmembrane domain have a different behavior in the bilayer. The native transmembrane sequence is found to be more flexible than in the presence of the Glu mutation, characterized by a tendency to pi deformation to accommodate the helix length to the membrane thickness. The mutant form of this domain does not evidence helical deformation in the present simulation. Hydrophobic matching is achieved both by a larger helix tilt and a vertical shift of the helix towards the membrane interface, favoring the accessibility of the Glu side chain to the membrane environment. A rapid exchange of hydrogen bond interactions with the surrounding water molecules and the lipid headgroups is observed. The difference in the behavior between the two peptides in a membrane environment was also observed experimentally. Both simulation and experimental results agree with the hypothesis that water may act as an intermediate for the formation of cross links between the facing Glu side chains stabilizing the dimer.

Orientation and Dynamics of Benzyl Alcohol and Benzyl Alkyl Ethers Dissolved in Nematic Lyotropic Liquid Crystals. 2H NMR and Molecular Dynamics Simulations

Most drugs have to cross cell membranes to reach their final target. A better understanding of the distribution, interactions, and dynamics of biologically active molecules in model bilayers is of fundamental importance in understanding drug functioning and design. 2H NMR quadrupole splittings (delta nu(Q)) and longitudinal relaxation times (T1) from the aromatic ring of benzyl alcohol-d5 (C0), a commonly used anesthetic, and a series of linear alkyl benzyl-d5 ethers with chain lengths from 1 to 12 carbon atoms (C1-C12), were measured. The molecules were dissolved in a nematic discotic lyotropic liquid crystal solution made of tetradecyltrimethylammonium chloride (TTAC)/decanol (DeOH)/NaCl/H2O. Values of delta nu(Q) and T1 from 1,1-dideuteriodecanol (15% enriched) and DHO (H2O with 0.2% D2O) were also measured. Delta nu(Q) of DeOH and DHO remained constant throughout the series. The value of delta nu(Q) of the para position of the ring (delta nu(p)) in C1 is 30% smaller than the delta nu(p) of C0. This is attributed to the existence of an H-bond between the alcohol hydroxyl proton and the solvent, which influences the average orientation of the ring. The relaxation data show that T1o,m is always longer than T1p and both decrease with the increase in alkyl chain length. Molecular dynamics simulations of the experimentally studied systems were performed. The aggregate was represented as a bilayer. The distribution, average orientation, and order parameters of the aromatic ring of the guest molecules in the bilayer were examined. Rotational correlation functions of all the C-D bonds and the OH bond from H2O were evaluated, allowing an estimate of the correlation times and T1. According to these results all spins relax in extreme narrowing conditions, except DeOH. Experimental and calculated T1 values differ at most by a factor of 3. However, the order of magnitude and the observed trends are well reproduced by the calculations. The aromatic ring of C0 possesses a unique average orientation in the bilayer. For the ether series, the orientation is modified and the C2 symmetry axis of the aromatic ring is exchanging between two orientations averaging the quadrupole splittings from the ortho and meta positions. The simulation supports the existence of an H-bond between C0 and the solvent not found in the ethers, which should be responsible for the observed differences.

Simulation studies of protein-induced bilayer deformations, and lipid-induced protein tilting, on a mesoscopic model for lipid bilayers with embedded proteins

Biological membranes are complex and highly cooperative structures. To relate biomembrane structure to their biological function it is often necessary to consider simpler systems. Lipid bilayers composed of one or two lipid species, and with embedded proteins, provide a model system for biological membranes. Here we present a mesoscopic model for lipid bilayers with embedded proteins, which we have studied with the help of the dissipative particle dynamics simulation technique. Because hydrophobic matching is believed to be one of the main physical mechanisms regulating lipid-protein interactions in membranes, we considered proteins of different hydrophobic length (as well as different sizes). We studied the cooperative behavior of the lipid-protein system at mesoscopic time- and lengthscales. In particular, we correlated in a systematic way the protein-induced bilayer perturbation, and the lipid-induced protein tilt, with the hydrophobic mismatch (positive and negative) between the protein hydrophobic length and the pure lipid bilayer hydrophobic thickness. The protein-induced bilayer perturbation was quantified in terms of a coherence length, xi(P), of the lipid bilayer hydrophobic thickness profile around the protein. The dependence on temperature of xi(P), and the protein tilt-angle, were studied above the main-transition temperature of the pure system, i.e., in the fluid phase. We found that xi(P) depends on mismatch, i.e., the higher the mismatch is, the longer xi(P) becomes, at least for positive values of mismatch; a dependence on the protein size appears as well. In the case of large model proteins experiencing extreme mismatch conditions, in the region next to the so-called lipid annulus, there appears an undershooting (or overshooting) region where the bilayer hydrophobic thickness is locally lower (or higher) than in the unperturbed bilayer, depending on whether the protein hydrophobic length is longer (or shorter) than the pure lipid bilayer hydrophobic thickness. Proteins may tilt when embedded in a too-thin bilayer. Our simulation data suggest that, when the embedded protein has a small size, the main mechanism to compensate for a large hydrophobic mismatch is the tilt, whereas large proteins react to negative mismatch by causing an increase of the hydrophobic thickness of the nearby bilayer. Furthermore, for the case of small, peptidelike proteins, we found the same type of functional dependence of the protein tilt-angle on mismatch, as was recently detected by fluorescence spectroscopy measurements.

Bilayer conformation of fusion peptide of influenza virus hemagglutinin: a molecular dynamics simulation study

Unraveling the conformation of membrane-bound viral fusion peptides is essential for understanding how those peptides destabilize the bilayer topology of lipids that is important for virus-cell membrane fusion. Here, molecular dynamics (MD) simulations were performed to investigate the conformation of the 20 amino acids long fusion peptide of influenza hemagglutinin of strain X31 bound to a dimyristoyl phosphatidylcholine (DMPC) bilayer. The simulations revealed that the peptide adopts a kinked conformation, in agreement with the NMR structures of a related peptide in detergent micelles. The peptide is located at the amphipathic interface between the headgroups and hydrocarbon chains of the lipid by an energetically favorable arrangement: The hydrophobic side chains of the peptides are embedded into the hydrophobic region and the hydrophilic side chains are in the headgroup region. The N-terminus of the peptide is localized close to the amphipathic interface. The molecular dynamics simulations also revealed that the peptide affects the surrounding bilayer structure. The average hydrophobic thickness of the lipid phase close to the N-terminus is reduced in comparison with the average hydrophobic thickness of a pure dimyristoyl phosphatidylcholine bilayer.

Molecular dynamics (MD) investigations of preformed structures of the transmembrane domain of the oncogenic Neu receptor dimer in a DMPC bilayer

The critical Val/Glu mutation in the membrane spanning domain of the rat Neu receptor confers the ability for ligand-independent signaling and leads to increased dimerization and transforming ability. There is evidence that the two transmembrane interacting helices play a role in receptor activation by imposing orientation constraints to the intracellular tyrosine kinase domains. By using MD simulations we have attempted to discriminate between correct and improper helix-helix packing by examining the structural and energetic properties of preformed left-handed and right-handed structures in a fully hydrated DMPC bilayer. The best energetic balance between the residues at the helix-helix interface and the residues exposed to the lipids is obtained for helices in symmetrical left-handed interactions packed together via Glu side chain/Ala backbone interhelical hydrogen bonds. Analyses demonstrate the importance of the ATVEG motif in helix-helix packing and point to additional contacting residues necessary for association. Our findings, all consistent with experimental data, suggest that a symmetrical left-handed structure of the helices could be the transmembrane domain configuration that promotes receptor activation and transformation. The present study may provide further insight into signal transduction mechanisms of the ErbB/Neu receptors.

Molecular dynamics simulations of the transmembrane domain of the oncogenic ErbB2 receptor dimer in a DMPC bilayer

Molecular dynamics simulations of an atomic model of the transmembrane domain of the oncogenic ErbB2 receptor dimer embedded in an explicit dimyristoylphosphatidylcholine (DMPC) bilayer were performed for more than 4 ns. The oncogenic Glu mutation in the membrane spanning segment plays a major role in tyrosine kinase activity and receptor dimerization, and is thought to be partly responsible for the structure of the transmembrane domain of the active receptor. MD results show that the interactions between the two transmembrane helices are characteristic of a left-handed packing as previously demonstrated from in vacuo simulations. Moreover, MD results reveal the absence of persistent hydrogen bonds between the Glu side chains in a membrane environment, which raise the question of the ability for Glu alone to stabilize the TM domain of the ErbB2 receptor. Interestingly the formation of the alpha-pi motif in the two ErbB2 transmembrane helices confirms the concept of intrinsic sequence-induced conformational flexibility. From a careful analysis of our MD results, we suggest that the left-handed helix-helix packing could be the key to correctly orient the intracellular domain of the activated receptor dimer. The prediction of such interactions from computer simulations represents a new step towards the understanding of signaling mechanisms.

Lipid-protein interactions in biological membranes: a structural perspective

Lipid molecules bound to membrane proteins are resolved in some high-resolution structures of membrane proteins. An analysis of these structures provides a framework within which to analyse the nature of lipid-protein interactions within membranes. Membrane proteins are surrounded by a shell or annulus of lipid molecules, equivalent to the solvent layer surrounding a water-soluble protein. The lipid bilayer extends right up to the membrane protein, with a uniform thickness around the protein. The surface of a membrane protein contains many shallow grooves and protrusions to which the fatty acyl chains of the surrounding lipids conform to provide tight packing into the membrane. An individual lipid molecule will remain in the annular shell around a protein for only a short period of time. Binding to the annular shell shows relatively little structural specificity. As well as the annular lipid, there is evidence for other lipid molecules bound between the transmembrane alpha-helices of the protein; these lipids are referred to as non-annular lipids. The average thickness of the hydrophobic domain of a membrane protein is about 29 A, with a few proteins having significantly smaller or greater thicknesses than the average. Hydrophobic mismatch between a membrane protein and the surrounding lipid bilayer generally leads to only small changes in membrane thickness. Possible adaptations in the protein to minimise mismatch include tilting of the helices and rotation of side chains at the ends of the helices. Packing of transmembrane alpha-helices is dependent on the chain length of the surrounding phospholipids. The function of membrane proteins is dependent on the thickness of the surrounding lipid bilayer, sometimes on the presence of specific, usually anionic, phospholipids, and sometimes on the phase of the phospholipid.

1H and (13)C NMR of multilamellar dispersions of polyunsaturated (22:6) phospholipids

The polyunsaturated fatty acid docosahexaenoic acid (DHA) makes up approximately 50% of the lipid chains in the retinal rod outer segment disk membranes and a large fraction of the lipid chains in the membranes of neuronal tissues. There is an extensive literature concerned with the dietary requirements for essential fatty acids and the importance of DHA to human health, but relatively little research has been done on the physical properties of this important molecule. Using (1)H and (13)C MAS NMR measurements of dispersions of 1-palmitoyl-2-docosahexaenoyl-phosphatidylcholine in excess phosphate buffer, we have unambiguously assigned most of the resonances in both the (1)H and (13)C NMR spectra. We were able to use cross-polarization spectroscopy to follow the transfer of polarization from specific (1)H nuclei not only to their directly bonded (13)C but also to those (13)C that are in close proximity, even though they are not directly bonded. Cross-peaks in two-dimensional cross-polarization spectra revealed a close association between the choline headgroup and at least part of the DHA chain but not with the palmitate chain. Finally, we examined the dynamics of the different parts of this lipid molecule, using rotating frame spin-lattice relaxation measurements, and found that methylene groups of both chains experience important motions with correlation times in the 10-micros range, with those for the palmitate chain being approximately 50% longer than those of the DHA chain. The choline headgroup and the chain terminal groups have significantly shorter correlation times, and that part of the dipolar interaction that is fluctuating at these correlation times is significantly smaller for these groups than it is for the palmitate and DHA chain methylenes.