A quantitative infrared determination of acyl chain conformation in gramicidin/dipalmitoylphosphatidylcholine mixtures

Coupling of Membrane Nanodomain Formation and Enhanced Electroporation near Phase Transition

Biological cells are enveloped by a heterogeneous lipid bilayer that prevents the uncontrolled exchange of substances between the cell interior and its environment. In particular, membranes act as a continuous barrier for salt and macromolecules to ensure proper physiological functions within the cell. However, it has been shown that membrane permeability strongly depends on temperature and, for phospholipid bilayers, displays a maximum at the transition between the gel and fluid phase. Here, extensive molecular dynamics simulations of dipalmitoylphosphatidylcholine bilayers were employed to characterize the membrane structure and dynamics close to phase transition, as well as its stability with respect to an external electric field. Atomistic simulations revealed the dynamic appearance and disappearance of spatially related nanometer-sized thick ordered and thin interdigitating domains in a fluid-like bilayer close to the phase transition temperature (Tm). These structures likely represent metastable precursors of the ripple phase that vanished at increased temperatures. Similarly, a two-phase bilayer with coexisting gel and fluid domains featured a thickness minimum at the interface because of splaying and interdigitating lipids. For all systems, application of an external electric field revealed a reduced bilayer stability with respect to pore formation for temperatures close to Tm. Pore formation occurred exclusively in thin interdigitating membrane nanodomains. These findings provide a link between the increased membrane permeability and the structural heterogeneity close to phase transition.

Membrane phase transition during heating and cooling: molecular insight into reversible melting

With increasing temperature, lipid bilayers undergo a gel-fluid phase transition, which plays an essential role in many physiological phenomena. In the present work, this first-order phase transition was investigated for variable heating and cooling rates for a dipalmitoylphosphatidylcholine (DPPC) lipid bilayer by means of atomistic molecular dynamics simulations. Alternative methods to track the melting temperature [Formula: see text] are compared. The resulting [Formula: see text] is shown to be independent of the scan rate for small heating rates (0.05-0.3 K/ns) implying reversible melting, and increases for larger heating (0.3-4 K/ns) or cooling rates (2-0.1 K/ns). The reported dependency of the melting temperature on the heating rate is in perfect agreement with a two-state kinetic rate model as suggested previously. Expansion and shrinkage, as well as the dynamics of melting seeds is described. The simulations further exhibit a relative shift between melting seeds in opposing membrane leaflets as predicted from continuum elastic theory.

Electric field driven changes of a gramicidin containing lipid bilayer supported on a Au(111) surface

Langmuir-Blodgett and Langmuir-Schaeffer methods were employed to deposit a mixed bilayer consisting of 90% of 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) and 10% of gramicidin (GD), a short 15 residue ion channel forming peptide, onto a Au(111) electrode surface. This architecture allowed us to investigate the effect of the electrostatic potential applied to the electrode on the orientation and conformation of DMPC molecules in the bilayer containing the ion channel. The charge density data were determined from chronocoulometry experiments. The electric field and the potential across the membrane were determined through the use of charge density curves. The magnitudes of potentials across the gold-supported biomimetic membrane were comparable to the transmembrane potential acting on a natural membrane. The information regarding the orientation and conformation of DMPC and GD molecules in the bilayer was obtained from photon polarization modulation infrared reflection absorption spectroscopy (PMIRRAS) measurements. The results show that the bilayer is adsorbed, in direct contact with the metal surface, when the potential across the interface is more positive than -0.4 V and is lifted from the gold surface when the potential across the interface is more negative than -0.4 V. This change in the state of the bilayer has a significant impact on the orientation and conformation of the phospholipid and gramicidin molecules. The potential induced changes in the membrane containing peptide were compared to the changes in the structure of the pure DMPC bilayer determined in earlier studies.

Infrared studies of protein-induced perturbation of lipids in lipoproteins and membranes

The paper reviews the main recent publications concerning infrared (IR) spectroscopy as applied to the study of lipid-protein interactions in model and cell membranes, lipoproteins, and related systems (e.g. lung surfactant). The review focuses mainly on transmission IR. Based on the available data, a number of general conclusions are presented on the perturbations caused by proteins on either the hydrocarbon chains, the polar headgroups or the interface region. Lipid-protein interactions in native cell membranes do not reveal significant differences from what is observed in semisynthetic model systems.

Structure, energetics, and dynamics of lipid-protein interactions: A molecular dynamics study of the gramicidin A channel in a DMPC bilayer

The microscopic details of lipid-protein interactions are examined using molecular dynamics simulations of the gramicidin A channel embedded in a fully hydrated dimyristoyl phosphatidylcholine (DMPC) bilayer. A novel construction protocol was used to assemble the initial configurations of the membrane protein complex for the simulations. Three hundred systems were constructed with different initial lipid placement and conformations. Seven systems were simulated with molecular dynamics. One system was simulated for a total of 600 psec, four were simulated for 300 psec, and two for 100 psec. Analysis of the resulting trajectories shows that the bulk solvent-membrane interface region is much broader than traditionally pictured in simplified continuum theories: its width is almost 15 angstroms. In addition, lipid-protein interactions are far more varied, both structurally and energetically, than is usually assumed: the total interaction energy between the gramicidin A and the individual lipids varies from 0 to -50 kcal/mol. The deuterium quadrupolar splittings of the lipid acyl chains calculated from the trajectories are in good agreement with experimental data. The lipid chains in direct contact with the GA are ordered but the effect is not uniform due to the irregular surface of the protein. Energy decompositions shows that the most energetically favorable interactions between lipid and protein involve nearly equal contributions from van der Waals and electrostatic interactions. The tryptophans, located near the bulk-membrane interface, appear to be particularly important in mediating both hydrogen bonding interactions with the lipid glycerol backbone and water and also in forming favorable van der Waals contacts with the hydrocarbon chains. In contrast, the interactions of the leucine residues with the lipids, also located near the interface, are dominated by van der Waals interactions with the hydrocarbon lipid chains.

Molecular dynamics simulation of the gramicidin channel in a phospholipid bilayer

A molecular dynamics simulation of the gramicidin A channel in an explicit dimyristoyl phosphatidylcholine bilayer was generated to study the details of lipid-protein interactions at the microscopic level. Solid-state NMR properties of the channel averaged over the 500-psec trajectory are in excellent agreement with available experimental data. In contrast with the assumptions of macroscopic models, the membrane/solution interface region is found to be at least 12 A thick. The tryptophan side chains, located within the interface, are found to form hydrogen bonds with the ester carbonyl groups of the lipids and with water, suggesting their important contribution to the stability of membrane proteins. Individual lipid-protein interactions are seen to vary from near 0 to -50 kcal/mol. The most strongly interacting conformations are short-lived and have a nearly equal contribution from both van der Waals and electrostatic energies. This approach for performing molecular dynamics simulations of membrane proteins in explicit phospholipid bilayers should help in studying the structure, dynamics, and energetics of lipid-protein interactions.

Solvent history dependence of gramicidin-lipid interactions: a Raman and infrared spectroscopic study

We have investigated the interactions between gramicidin and a model membrane composed of one phospholipid, dimyristoylphosphatidylcholine, as a function of the cosolubilization solvent and incubation time used in the sample preparation. Three organic solvents have been used; trifluoroethanol, a mixture of methanol/chloroform (1:1 v/v), and ethanol. Using Fourier transform infrared spectroscopy, we have demonstrated that the conformation adopted by gramicidin in the membrane is dependent upon the cosolubilization solvent used, and, only with trifluoroethanol, it is possible to incorporate gramicidin entirely as a beta 6.3-helix. Moreover, Raman spectroscopy results indicate that the orientation of the tryptophan side chains in gramicidin and their interaction with the hydrocarbon chains and the carbonyl groups of the lipids are also dependent on the cosolubilization solvent. On the other hand, the effect of the incorporation of gramicidin on the thermotropism of the lipid bilayer was found to be dependent upon the conformation of gramicidin in the lipid bilayers.

Acyl chain conformational ordering in 1,2 dipalmitoylphosphatidylethanolamine. Integration of FT-IR and 2H NMR results

The extent of trans-gauche isomerization at the 4 and 4' positions of the acyl chains of fully hydrated 4,4,4',4'-d4 1,2-dipalmitoylphosphatidylethanolamine (4-d4 DPPE) bilayers was quantitatively evaluated from the infrared (IR) intensity of the CD2 rocking modes. About 20% gauche conformers were observed at 72 degrees C (above Tm), while at 23 degrees C, well below Tm, about 4% were noted. The order parameter SC-D was determined from 2H nuclear magnetic resonance (NMR) quadrupolar splittings. SC-D is the product of a segmental order parameter (S gamma), which depends on conformational order, and a chain order parameter (S alpha) which depends on slower motions such as chain wobble. The IR-determined percentage of gauche forms was converted into a segmental order parameter and factored out of the measured value for SC-D to yield an estimate of S alpha = 0.59 for L alpha phase DPPE. A comparison with S alpha for 1,2-dipalmitoylphosphatidylcholine (DPPC) suggests that increased wobble is responsible for enhanced motional averaging of the quadrupolar splittings in the latter at a similar reduced temperature. The extent of conformational disordering [at the 4(4') position] is essentially unchanged between the two molecules. The current study demonstrates the advantage of integrating quantitative IR with 2H NMR data, for elucidation of the contributions of the individual motions that average the NMR quadrupolar splittings.

Model ion channels: gramicidin and alamethicin

We have discussed in some detail a variety of experimental studies which were designed to elucidate the conformational and dynamic properties of gramicidin and alamethicin. Although the behavior of these peptides is by no means fully characterized, these studies have already permitted aspects of ion channel activity to be understood in molecular terms. Studies with gramicidin in a variety of organic solutions have revealed conformational heterogeneity of this peptide; at least five major isomers exist, several of which have been characterized in detail using NMR spectroscopy and X-ray crystallography. When added to lipid membranes gramicidin undergoes a further conformational conversion. Although the conformation of gramicidin in membranes is not as well characterized as the solution conformation(s) and an X-ray structure is not yet available, detailed data, particularly from solid-state NMR studies, continue to become available and a right-handed beta 6.3 helical conformation of the peptide backbone is now generally accepted. Two of these beta 6.3 helices joined at their N-termini are believed to form the conducting channel. The conformational behavior of the side-chains of gramicidin in the membrane-bound form is not well established and several NMR, CD, fluorescence and theoretical studies are now focussed on this. Although the side-chains do not directly contact the permeating ions, they can have distinct effects on conductance and selectivity by altering the electrostatic environment sensed by the ion. The dynamics of both side-chain and backbone conformations of gramicidin appear critical to a detailed understanding of the ion transport process in this channel. As the description of the membrane-bound conformation of gramicidin becomes more detailed, simulations of ion transport using computational methods are likely to improve and will further our understanding of the processes of ion transport. As well as internal motion of the backbone and side-chains, gramicidin undergoes rotational and translational motion in the plane of the membrane. These motions do not appear to be essential for the process of ion transport but can affect channel lifetime since lifetime is determined by the rate of association and dissociation of gramicidin monomers. Gramicidin-membrane interactions are also likely to be involved in the frequency of occurrence of channel subconductance states, the frequency of channel flickering and fundamentally in the stability of the membrane-bound gramicidin conformation. Alamethicin forms channels in membranes which are strongly voltage-dependent. The molecular origin of voltage-dependent conductances has been a fundamental problem in biophysics for many years.(ABSTRACT TRUNCATED AT 400 WORDS)

CD2 rocking modes as quantitative infrared probes of one-, two-, and three-bond conformational disorder in dipalmitoylphosphatidylcholine and dipalmitoylphosphatidylcholine/cholesterol mixtures

The use of CD2 rocking modes in the IR spectrum as quantitative probes of phospholipid conformational disorder has recently been described for aqueous dispersions of 1,2-dipalmitoylphosphatidylcholine (DPPC) and DPPC/cholesterol mixtures [Mendelsohn et al. (1989) Biochemistry 28, 8934-8939; Davies et al. (1990) Biochemistry 29, 4368-4373]. Initial studies focused at the 4, 6, and 10 acyl chain positions of DPPC. In the current work, the method is extended to the 2, 3, 12, and 13 positions. Conformational disorder in the L alpha phase is approximately the same (about 20% gauche) at positions 4, 10, and 13, but an unexpected higher value is observed (about 30%) at the 6 position. Cholesterol (33 mol%) restricts gauche rotamer formation by factors ranging from 6 to 9 at positions 4 and 6, respectively, to 1.5-2 at positions 10, 12, and 13. Quantitative analysis for the DPPC/cholesterol "liquid-ordered" phase indicates the occurrence of 1.2 gauche bonds/chain, a marked reduction from the 3.6-4.2 gauche bonds/chain for DPPC alone. Proximity to the ester moiety at acyl chain position 3 perturbs the vibrational coupling patterns of the CD2 rocking modes and eliminates their sensitivity to conformational change. In addition, the feasibility of a method based on the conformation-dependent coupling between CD2 rocking frequencies of two successive CD2 groups for the quantitative detection of specific, position-dependent king (gtg') and isolated gauche (gtt) conformers is demonstrated. Finally, comparisons between IR measurements and explicit theoretical predictions of acyl chain conformational order are presented.