Dataset Information

Interpretation of fluctuation spectra in lipid bilayer simulations.

ABSTRACT: Atomic resolution and coarse-grained simulations of dimyristoylphosphatidylcholine lipid bilayers were analyzed for fluctuations perpendicular to the bilayer using a completely Fourier-based method. We find that the fluctuation spectrum of motions perpendicular to the bilayer can be decomposed into just two parts: 1), a pure undulation spectrum proportional to q(-4) that dominates in the small-q regime; and 2), a molecular density structure factor contribution that dominates in the large-q regime. There is no need for a term proportional to q(-2) that has been postulated for protrusion fluctuations and that appeared to have been necessary to fit the spectrum for intermediate q. We suggest that earlier reports of such a term were due to the artifact of binning and smoothing in real space before obtaining the Fourier spectrum. The observability of an intermediate protrusion regime from the fluctuation spectrum is discussed based on measured and calculated material constants.

SUBMITTER: Brandt EG

PROVIDER: S-EPMC3149257 | biostudies-literature | 2011 May

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Interpretation of fluctuation spectra in lipid bilayer simulations.

Brandt Erik G EG Braun Anthony R AR Sachs Jonathan N JN Nagle John F JF Edholm Olle O

Biophysical journal 20110501 9

Atomic resolution and coarse-grained simulations of dimyristoylphosphatidylcholine lipid bilayers were analyzed for fluctuations perpendicular to the bilayer using a completely Fourier-based method. We find that the fluctuation spectrum of motions perpendicular to the bilayer can be decomposed into just two parts: 1), a pure undulation spectrum proportional to q(-4) that dominates in the small-q regime; and 2), a molecular density structure factor contribution that dominates in the large-q regim ...[more]

PMID: 21539777

Similar Datasets

Project description:Molecular dynamics (MD) simulations have become an indispensable tool to investigate phase separation in model membrane systems. In particular, simulations based on coarse-grained (CG) models have found widespread use due to their increased computational efficiency, allowing for simulations of multicomponent lipid bilayers undergoing phase separation into liquid-ordered and liquid-disordered domains. Here, we show that a significant temperature difference between molecule types can artificially arise in CG MD membrane simulations with the standard Martini simulation parameters in GROMACS. In particular, the linear constraint solver (LINCS) algorithm does not converge with its default settings, resulting in serious temperature differences between molecules in a time step-dependent manner. We demonstrate that the underlying reason for this behavior is the presence of highly constrained moieties, such as cholesterol. Their presence can critically impact numerous structural and dynamic membrane properties obtained from such simulations. Furthermore, any preference of these molecules toward a certain membrane phase can lead to spatial temperature gradients, which can amplify the degree of phase separation or even induce it in compositions that would otherwise mix well. We systematically investigated the effect of the integration time step and LINCS settings on membrane properties. Our data show that for cholesterol-containing membranes, a time step of 20 fs should be combined with at least lincs_iter = 2 and lincs_order = 12, while using a time step of 30 fs requires at least lincs_iter = 3 and lincs_order = 12 to bring the temperature differences to a level where they do not perturb central membrane properties. Moreover, we show that in cases where stricter LINCS settings are computationally too demanding, coupling the lipids in multiple groups to the temperature bath offers a practical workaround to the problem, although the validity of this approach should be further verified. Finally, we show that similar temperature gradients can also emerge in atomistic simulations using the CHARMM force field in combination with settings that allow for a 5 fs integration step.

Project description:Spin-labeled lipids are commonly used as fluorescence quenchers in studies of membrane penetration of dye-labeled proteins and peptides using depth-dependent quenching. Accurate calculations of depth of the fluorophore rely on the use of several spin labels placed in the membrane at various positions. The depth of the quenchers (spin probes) has to be determined independently; however, experimental determination of transverse distributions of spin probe depths is difficult. In this Article, we use molecular dynamics (MD) simulations to study the membrane behavior and depth distributions of spin-labeled phospholipids in a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayer. To probe different depths within the bilayer, a series containing five Doxyl-labeled lipids (n-Doxyl PC) has been studied, in which a spin moiety was covalently attached to nth carbon atoms (where n = 5, 7, 10, 12, and 14) of the sn-2 stearoyl chain of the host phospholipid. Our results demonstrate that the chain-attached spin labels are broadly distributed across the model membrane and their environment is characterized by a high degree of mobility and structural heterogeneity. Despite the high thermal disorder, the depth distributions of the Doxyl labels were found to correlate well with their attachment positions, indicating that the distribution of the spin label within the model membrane is dictated by the depth of the nth lipid carbon atom and not by intrinsic properties of the label. In contrast, a much broader and heterogeneous distribution was observed for a headgroup-attached Tempo spin label of Tempo-PC lipids. MD simulations reveal that, due to the hydrophobic nature, a Tempo moiety favors partitioning from the headgroup region deeper into the membrane. Depending on the concentration of Tempo-PC lipids, the probable depth of the Tempo moiety could span a range from 14.4 to 18.2 Å from the membrane center. Comparison of the MD-estimated immersion depths of Tempo and n-Doxyl labels with their suggested experimental depth positions allows us to review critically the possible sources of error in depth-dependent fluorescence quenching studies.

Project description:Solid-state NMR spectroscopy has emerged as an excellent tool to study the structure and dynamics of membrane proteins under native-like conditions in lipid bilayers. One of the key considerations in experimental design is the uniaxial rotational diffusion of the protein that can affect the NMR spectral observables. In this regard, temperature plays a fundamental role in modulating the phase properties of the lipids, which directly influences the rotational diffusion rate of the protein in the bilayer. In fact, it is well established that below the main phase transition temperature of the lipid bilayer the protein's motion is significantly slowed while above this critical temperature the rate is increased. In this article, we carried out a systematic comparison of the signal intensity and spectral resolution as a function of temperature using magic-angle-spinning (MAS) solid-state NMR spectroscopy. These observables were directly correlated with the relative fluidity of the lipid bilayer as inferred from differential scanning calorimetry (DSC). We applied our hybrid biophysical approach to two polytopic membrane proteins from the small multidrug resistance family (EmrE and SugE) reconstituted into model membrane lipid bilayers (DMPC-14:0 and DPPC-16:0). From these experiments, we conclude that the rotational diffusion giving optimal spectral resolution occurs at a bilayer fluidity of ~5%, which corresponds to the percentage of lipids in the fluid or liquid-crystalline fraction. At the temperature corresponding to this critical value of fluidity, there is sufficient mobility to reduce inhomogeneous line broadening that occurs at lower temperatures. A greater extent of fluidity leads to faster uniaxial rotational diffusion and a sigmoidal-type reduction in the NMR signal intensity, which stems from intermediate-exchange dynamics where the motion has a similar frequency as the NMR observables (i.e., dipolar couplings and chemical shift anisotropy). These experiments provide insight into the optimal temperature range and corresponding bilayer fluidity to study membrane proteins by solid-state NMR spectroscopy.

Dataset Information

Interpretation of fluctuation spectra in lipid bilayer simulations.

Publications

Interpretation of fluctuation spectra in lipid bilayer simulations.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets