O. H. Samuli Ollila

Toward Atomistic Resolution Structure of Phosphatidylcholine Headgroup and Glycerol Backbone at Different Ambient Conditions

Phospholipids are essential building blocks of biological membranes. Despite a vast amount of very accurate experimental data, the atomistic resolution structures sampled by the glycerol backbone and choline headgroup in phoshatidylcholine bilayers are not known. Atomistic resolution molecular dynamics simulations have the potential to resolve the structures, and to give an arrestingly intuitive interpretation of the experimental data, but only if the simulations reproduce the data within experimental accuracy. In the present work, we simulated phosphatidylcholine (PC) lipid bilayers with 13 different atomistic models, and compared simulations with NMR experiments in terms of the highly structurally sensitive C–H bond vector order parameters. Focusing on the glycerol backbone and choline headgroups, we showed that the order parameter comparison can be used to judge the atomistic resolution structural accuracy of the models. Accurate models, in turn, allow molecular dynamics simulations to be used as an interpretation tool that translates these NMR data into a dynamic three-dimensional representation of biomolecules in biologically relevant conditions. In addition to lipid bilayers in fully hydrated conditions, we reviewed previous experimental data for dehydrated bilayers and cholesterol-containing bilayers, and interpreted them with simulations. Although none of the existing models reached experimental accuracy, by critically comparing them we were able to distill relevant chemical information: (1) increase of choline order parameters indicates the P–N vector tilting more parallel to the membrane, and (2) cholesterol induces only minor changes to the PC (glycerol backbone) structure. This work has been done as a fully open collaboration, using nmrlipids.blogspot.fi as a communication platform; all the scientific contributions were made publicly on this blog. During the open research process, the repository holding our simulation trajectories and files (https://zenodo.org/collection/user-nmrlipids) has become the most extensive publicly available collection of molecular dynamics simulation trajectories of lipid bilayers.

Protein Shape Change Has a Major Effect on the Gating Energy of a Mechanosensitive Channel

Increasing experimental evidence has shown that membrane protein functionality depends on molecular composition of cell membranes. However, the origin of this dependence is not fully understood. It is reasonable to assume that specific lipid-protein interactions are important, yet more generic effects due to mechanical properties of lipid bilayers likely play a significant role too. Previously it has been demonstrated using models for elastic properties of membranes and lateral pressure profiles of lipid bilayers that the mechanical properties of a lipid bilayer can contribute as much as ∼10 k(B)T to the free energy difference associated with a change in protein conformational state. Here, we extend those previous approaches to a more realistic model for a large mechanosensitive channel (MscL). We use molecular dynamics together with the MARTINI model to simulate the open and closed states of MscL embedded in a DOPC bilayer. We introduce a procedure to calculate the mechanical energy change in the channel gating using a three-dimensional pressure distribution inside a membrane, computed from the molecular dynamics simulations. We decompose the mechanical energy to terms associated with area dilation and shape contribution. Our results highlight that the lateral pressure profile of a lipid bilayer together with the shape change in gating can induce a contribution of ∼30 k(B)T on the gating energy of MscL. This contribution arises largely from the interfacial tension between hydrophobic and hydrophilic regions in a lipid bilayer.

Interfacial tension and surface pressure of high density lipoprotein, low density lipoprotein, and related lipid droplets.

Lipid droplets play a central role in energy storage and metabolism on a cellular scale. Their core is comprised of hydrophobic lipids covered by a surface region consisting of amphiphilic lipids and proteins. For example, high and low density lipoproteins (HDL and LDL, respectively) are essentially lipid droplets surrounded by specific proteins, their main function being to transport cholesterol. Interfacial tension and surface pressure of these particles are of great interest because they are related to the shape and the stability of the droplets and to protein adsorption at the interface. Here we use coarse-grained molecular-dynamics simulations to consider a number of related issues by calculating the interfacial tension in protein-free lipid droplets, and in HDL and LDL particles mimicking physiological conditions. First, our results suggest that the curvature dependence of interfacial tension becomes significant for particles with a radius of ∼5 nm, when the area per molecule in the surface region is <1.4 nm(2). Further, interfacial tensions in the used HDL and LDL models are essentially unaffected by single apo-proteins at the surface. Finally, interfacial tensions of lipoproteins are higher than in thermodynamically stable droplets, suggesting that HDL and LDL are kinetically trapped into a metastable state.

Membrane simulations mimicking acidic pH reveal increased thickness and negative curvature in a bilayer consisting of lysophosphatidylcholines and free fatty acids.

Phospholipids are key components of biological membranes and their lipolysis with phospholipase A(2) (PLA(2)) enzymes occurs in different cellular pH environments. Since no studies are available on the effect of pH on PLA(2)-modified phospholipid membranes, we performed 50-ns atomistic molecular dynamics simulations at three different pH conditions (pH 9.0, 7.5, and 5.5) using a fully PLA(2)-hydrolyzed phosphatidylcholine (PC) bilayer which consists solely of lysophosphatidylcholine and free fatty acid molecules. We found that a decrease in pH results in lateral squeezing of the membrane, i.e. in decreased surface area per headgroup. Thus, at the decreased pH, the lipid hydrocarbon chains had larger S(CD) order parameter values, and also enhanced membrane thickness, as seen in the electron density profiles across the membrane. From the lateral pressure profiles, we found that the values of spontaneous curvature of the two opposing monolayers became negative when the pH was decreased. At low pH, protonation of the free fatty acid headgroups reduces their mutual repulsion and accounts for the pH dependence of all the above-mentioned properties. The altered structural characteristics may significantly affect the overall surface properties of biomembranes in cellular vesicles, lipid droplets, and plasma lipoproteins, play an important role in membrane fission and fusion, and modify interactions between membrane lipids and the proteins embedded within them.

Asymmetric nature of lateral pressure profiles in supported lipid membranes and its implications for membrane protein functions

The study of lipid membrane structure and phase behavior at the nano-scale is of utmost importance due to implications in understanding the interplay of lipids and membrane proteins in biochemical membrane processes. Supported lipid bilayers play a key role in providing the means to understand real biological systems. Yet, with regard to membrane protein activation and their functions, the limitations of supported membranes and the possible artifacts they may induce are weakly understood. Here, we study DPPC (dipalmitoylphosphatidylcholine) bilayers on a weakly hydrophilic substrate using coarse grained simulations and show that the lateral pressure profile of a supported bilayer is distinctly different from the corresponding pressure profile for a free-standing DPPC bilayer. The results indicate that due to the substrate, the lateral pressure profile becomes asymmetric, expressing major peaks at the proximal leaflet of the membrane, implying the membrane to be under strong tension. The results provide a new mechanism to explain malfunctions of transmembrane proteins used in supported bilayers.

Model-free estimation of the effective correlation time for C-H bond reorientation in amphiphilic bilayers: (1)H-(13)C solid-state NMR and MD simulations.

Molecular dynamics (MD) simulations give atomically detailed information on structure and dynamics in amphiphilic bilayer systems on timescales up to about 1 μs. The reorientational dynamics of the C-H bonds is conventionally verified by measurements of (13)C or (2)H nuclear magnetic resonance (NMR) longitudinal relaxation rates R1, which are more sensitive to motional processes with correlation times close to the inverse Larmor frequency, typically around 1-10 ns on standard NMR instrumentation, and are thus less sensitive to the 10-1000 ns timescale motion that can be observed in the MD simulations. We propose an experimental procedure for atomically resolved model-free estimation of the C-H bond effective reorientational correlation time τe, which includes contributions from the entire range of all-atom MD timescales and that can be calculated directly from the MD trajectories. The approach is based on measurements of (13)C R1 and R1ρ relaxation rates, as well as (1)H-(13)C dipolar couplings, and is applicable to anisotropic liquid crystalline lipid or surfactant systems using a conventional solid-state NMR spectrometer and samples with natural isotopic composition. The procedure is demonstrated on a fully hydrated lamellar phase of 1-palmitoyl-2-oleoyl-phosphatidylcholine, yielding values of τe from 0.1 ns for the methyl groups in the choline moiety and at the end of the acyl chains to 3 ns for the g1 methylene group of the glycerol backbone. MD simulations performed with a widely used united-atom force-field reproduce the τe-profile of the major part of the acyl chains but underestimate the dynamics of the glycerol backbone and adjacent molecular segments. The measurement of experimental τe-profiles can be used to study subtle effects on C-H bond reorientational motions in anisotropic liquid crystals, as well as to validate the C-H bond reorientation dynamics predicted in MD simulations of amphiphilic bilayers such as lipid membranes.

Molecular conformation and bilayer pores in a nonionic surfactant lamellar phase studied with (1)h-(13)c solid-state NMR and molecular dynamics simulations.

The structure of the lamellar phase of aqueous pentaethylene glycol mono-n-dodecyl ether (C12E5) surfactant at various temperatures and molar fractions is studied by using united atom molecular dynamics simulations and nuclear magnetic resonance measurements. Namely, the simulation model is used to interpret the magnitude and temperature dependence of experimental C-H order parameter profiles in terms of the molecular conformation and orientation. Our simulations suggest that the low order parameters that are generally measured in poly(ethylene oxide) surfactant bilayers are due to the presence of bilayer pores throughout the entire lamellar phase region.

Acyl Chain Disorder and Azelaoyl Orientation in Lipid Membranes Containing Oxidized Lipids

Oxidized phospholipids occur naturally in conditions of oxidative stress and have been suggested to play an important role in a number of pathological conditions due to their effects on a lipid membrane acyl chain orientation, ordering, and permeability. Here we investigate the effect of the oxidized phospholipid 1-palmitoyl-2-azelaoyl-sn-glycero-3-phosphocholine (PazePC) on a model membrane of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) using a combination of (13)C-(1)H dipolar-recoupling nuclear magnetic resonance (NMR) experiments and united-atom molecular dynamics (MD) simulations. The obtained experimental order parameter SCH profiles show that the presence of 30 mol % PazePC in the bilayer significantly increases the gauche content of the POPC acyl chains, therefore decreasing the thickness of the bilayer, although with no stable bilayer pore formation. The MD simulations reproduce the disordering effect and indicate that the orientation of the azelaoyl chain is highly dependent on its protonation state with acyl chain reversal for fully deprotonated states and a parallel orientation along the interfacial plane for fully protonated states, deprotonated and protonated azelaoyl chains having negative and positive SCH profiles, respectively. Only fully or nearly fully protonated azelaoyl chain are observed in the (13)C-(1)H dipolar-recoupling NMR experiments. The experiments show positive SCH values for the azelaoyl segments confirming for the first time that oxidized chains with polar termini adopt a parallel orientation to the bilayer plane as predicted in MD simulations.

Can pyrene probes be used to measure lateral pressure profiles of lipid membranes? Perspective through atomistic simulations.

The lateral pressure profile of lipid bilayers has gained a lot of attention, since changes in the pressure profile have been suggested to shift the membrane protein conformational equilibrium. This relation has been mostly studied with theoretical methods, especially with molecular dynamics simulations, since established methods to measure the lateral pressure profile experimentally have not been available. The only experiments that have attempted to gauge the lateral pressure profile have been done by using di-pyrenyl-phosphatidylcholine (di-pyr-PC) probes. In these experiments, the excimer/monomer fluorescence ratio has been assumed to represent the lateral pressure in the location of the pyrene moieties. Here, we consider the validity of this assumption through atomistic molecular dynamics simulations in a DOPC (dioleoylphosphatidylcholine) membrane, which hosts di-pyr-PC probes with different acyl chain lengths. Based on the simulations, we calculate the pyrene dimerization rate and the lateral pressure at the location of the pyrenes. The dimerization rates are compared with the results of di-pyr-PC probes simulated in vacuum. The comparison indicates that the lateral pressure is not the dominant determinant of the excimer/monomer fluorescence ratio. Thus, the results do not support the usage of di-pyr-PC molecules to measure the shape of the lateral pressure profile. We yet discuss how the probes could potentially be exploited to gain qualitative insight of the changes in pressure profile when lipid composition is altered.

Chapter 2:Lateral Pressure Profiles in Lipid Membranes: Dependence on Molecular Composition

The lateral pressure profile is one of the central quantities of lipid membranes, as it governs membrane elasticity and consequently the activation of certain membrane protein types such as mechanosensitive channels. Nonetheless, the understanding of the dependence of the pressure profile on molecular membrane composition is limited. In this review, we discuss the present understanding of pressure profiles in model membranes, discussing data from both computational and experimental approaches that have dealt with this complex physical observable. We consider pressure profiles in both single- and many-component lipid membranes, paying attention to effects due to polyunsaturated lipids, sterols, and anesthetics. The role of the pressure profile in membrane elasticity and membrane protein activation is given particular attention.