Area per lipid and acyl length distributions in fluid phosphatidylcholines determined by (2)H NMR spectroscopy

Single-Molecule Optical Tweezers Study of Protein-Membrane Interactions

Numerous proteins directly or indirectly bind membranes to exert their roles in a wide variety of biological processes. Such membrane binding often occurs in the presence of an external mechanical force. It remains challenging to quantify these interactions using traditional experimental approaches based on a large number of molecules, due to ensemble averaging or the lack of mechanical force. Here we described a new single-molecule approach based on high-resolution optical tweezers to characterize protein-membrane interactions. A single membrane binding protein is attached to the lipid bilayer coated on a silica bead via a flexible polypeptide linker, tethered to another bead via a long DNA handle, and pulled away from the bilayer using optical tweezers. Dynamic protein binding and unbinding is detected by the corresponding changes in the extension of the protein-DNA tether with high spatiotemporal resolution, which reveals the membrane binding affinity, kinetics, and intermediates. We demonstrated the method using C2 domains of extended synaptotagmin 2 (E-Syt2) with a detailed protocol. The method can be widely applied to investigate complex protein-membrane interactions under well-controlled experimental conditions.

Dynamic response on a nanometer scale of binary phospholipid-cholesterol vesicles: Low-frequency Raman scattering insight

Low-frequency Raman spectroscopy was used to study the dynamic response on a nanometer scale of aqueous suspensions of two-component lipid vesicles. Binary mixtures of saturated phospholipid (1,2-dipalmitoyl-sn-glycero-3-phosphocholine, DPPC) and cholesterol are interesting for possible coexistence of solidlike and liquid-ordered phases, while the phase coexistence was not reported for unsaturated phospholipid (1,2-dioleoyl-sn-glycero-3-phosphocholine, DOPC) and cholesterol mixtures. The DOPC-DPPC mixtures represent the well-documented case of coexisting domains of solidlike and liquid-disordered phases. These three series of lipid mixtures are studied here. A broad peak with the maximum in the range of 30-50cm^{-1} and a narrow peak near 10cm^{-1} are observed in the Raman susceptibility of the binary mixtures and attributed to the acousticlike vibrational density of states and layer modes, respectively. Parameters of the broad and narrow peaks are sensitive to lateral and conformational hydrocarbon chain ordering. It was also demonstrated that the low-frequency Raman susceptibility of multicomponent lipid bilayers allows one to determine the phase state of lipid bilayers and distinguish the homogeneous distribution of molecular complexes from coexisting domains with sizes above several nanometers. Thus, the low-frequency Raman spectroscopy provides unique information in studying phase coexistence in lipid bilayers.

Deuteration for High-Resolution Detection of Protons in Protein Magic Angle Spinning (MAS) Solid-State NMR

Proton detection developed in the last 20 years as the method of choice to study biomolecules in the solid state. In perdeuterated proteins, proton dipolar interactions are strongly attenuated, which allows yielding of high-resolution proton spectra. Perdeuteration and backsubstitution of exchangeable protons is essential if samples are rotated with MAS rotation frequencies below 60 kHz. Protonated samples can be investigated directly without spin dilution using proton detection methods in case the MAS frequency exceeds 110 kHz. This review summarizes labeling strategies and the spectroscopic methods to perform experiments that yield assignments, quantitative information on structure, and dynamics using perdeuterated samples. Techniques for solvent suppression, H/D exchange, and deuterium spectroscopy are discussed. Finally, experimental and theoretical results that allow estimation of the sensitivity of proton detected experiments as a function of the MAS frequency and the external B₀ field in a perdeuterated environment are compiled.

Molecular insights into the uptake of SiO2 nanoparticles on phospholipid membrane: Effect of surface properties and particle size

Silica (SiO₂) nanoparticles (NPs) have attracted much attention due to the potential for a wide range of applications and they have been confirmed to be hazardous to humans. Partitioning to phospholipid bilayers is an important way for their bioaccumulation. However, the detailed mechanism of SiO₂ NPs uptake by membrane phospholipids remains uncertain. In this work, molecular dynamics (MD) simulations were employed to study the uptake process of SiO₂ NPs into DPPC bilayers. Results show that these SiO₂ NPs uptake onto DPPC bilayer surface is favorable from the viewpoint of thermodynamics. During the uptake process, the SiO₂ NP needed to adjust the angle of interaction with the DPPC surface until the most stable adsorption configuration was reached. After incorporating into DPPC bilayers, the interaction between PO₄^- group and SiO₂ particle is stronger than -N⁺(CH₃)₃ group and SiO₂. Small SiO₂ NP was found to adsorb to the surface of DPPC bilayer without disturbing the morphology or membrane. In contract, bioaccumulation of large SiO₂ NP to DPPC induced a strong local membrane deformation. In addition, the effect of SiO₂ NP surface functionalization on its interaction with DPPC was also investigated. This molecular-level study reports a complete description of the interaction between SiO₂ NPs and DPPC bilayer, aiming to provide some insights for the further work on the bioaccumulation and hemolytic activity of SiO₂ NPs.

Room temperature XFEL crystallography reveals asymmetry in the vicinity of the two phylloquinones in photosystem I

Photosystem I (PS I) has a symmetric structure with two highly similar branches of pigments at the center that are involved in electron transfer, but shows very different efficiency along the two branches. We have determined the structure of cyanobacterial PS I at room temperature (RT) using femtosecond X-ray pulses from an X-ray free electron laser (XFEL) that shows a clear expansion of the entire protein complex in the direction of the membrane plane, when compared to previous cryogenic structures. This trend was observed by complementary datasets taken at multiple XFEL beamlines. In the RT structure of PS I, we also observe conformational differences between the two branches in the reaction center around the secondary electron acceptors A_1A and A_1B. The π-stacked Phe residues are rotated with a more parallel orientation in the A-branch and an almost perpendicular confirmation in the B-branch, and the symmetry breaking PsaB-Trp673 is tilted and further away from A_1A. These changes increase the asymmetry between the branches and may provide insights into the preferential directionality of electron transfer.

Structure of benzothiadiazine at zwitterionic phospholipid cell membranes

The use of drugs derived from benzothiadiazine, which is a bicyclic heterocyclic benzene derivative, has become a widespread treatment for diseases such as hypertension (treated with diuretics such as bendroflumethiazide or chlorothiazide), low blood sugar (treated with non-diuretic diazoxide), or the human immunodeficiency virus, among others. In this work, we have investigated the interactions of benzothiadiazine with the basic components of cell membranes and solvents, such as phospholipids, cholesterol, ions, and water. The analysis of the mutual microscopic interactions is of central importance to elucidate the local structure of benzothiadiazine as well as the mechanisms responsible for the access of benzothiadiazine to the interior of the cell. We have performed molecular dynamics simulations of benzothiadiazine embedded in three different model zwitterionic bilayer membranes made by dimyristoylphosphatidylcholine, dioleoylphosphatidylcholine, 1,2-dioleoyl-sn-glycero-3-phosphoserine, and cholesterol inside aqueous sodium-chloride solution in order to systematically examine microscopic interactions of benzothiadiazine with the cell membrane at liquid-crystalline phase conditions. From data obtained through radial distribution functions, hydrogen-bonding lengths, and potentials of mean force based on reversible work calculations, we have observed that benzothiadiazine has a strong affinity to stay at the cell membrane interface although it can be fully solvated by water in short periods of time. Furthermore, benzothiadiazine is able to bind lipids and cholesterol chains by means of single and double hydrogen-bonds of different characteristic lengths.

Study of lipid heterogeneity on bilayer membranes using molecular dynamics simulations

Human cell membranes consist of various lipids that are essential for their structure and function. It typically comprises phosphatidylcholine (POPC), phosphatidylethanolamine (POPE), phosphatidylserine (POPS), sphingomyelin (PSM), and cholesterol (CHL). Several experimental and computational techniques have been employed to characterize the composition of human cell membranes, however, CHL enriched membrane is still not clearly understood through these techniques. Molecular dynamics simulation results illustrated the biophysical properties of heterogeneous membranes based on the lipid composition as well as the concentration of lipids, exclusively for CHL and PSM. Herein, we have investigated the structure-function relationships of lipids comparatively to delineate the effect of heterogeneity on the biophysical properties of different membranes. It has been observed that the significant fraction of CHL (i.e., ~33% in ternary, ~25% in quaternary, and ~16% in senary type bilayers) in combination with other lipids introduced compactness, and increased the thickness of the membrane. The analysis of lipid mass density stated that the density of lipid head group, phosphate_, and glycerol-ester in presence of CHL with or without PSM is an underlying reason for membrane ordering. Results also revealed that the presence of POPI and POPS are the reasons for an adequate drop in the ordering of lipid chain, particularly on POPE chain. The self-interaction of CHL, PSM, POPE and the interaction of CHL and POPC with POPE seem to determine the structure and function of the heterogeneous membrane. Our findings provide a qualitative understanding of the effect of membrane heterogeneity on the physiological properties of membranes. The structures inspected in this study would help to select the heterogeneous bilayer model to mimic the human cell membranes to analyse or characterize the membrane-associated phenomena.

Impact of Selected Small-Molecule Kinase Inhibitors on Lipid Membranes

Small-molecule protein kinase inhibitors are used for the treatment of various diseases. Although their effect(s) on the respective kinase are generally quite well understood, surprisingly, their interaction with membranes is only barely investigated; even though these drugs necessarily come into contact with the plasma and intracellular membranes. Using biophysical methods such as NMR, ESR, and fluorescence spectroscopy in combination with lipid vesicles, we studied the membrane interaction of the kinase inhibitors sunitinib, erlotinib, idelalisib, and lenvatinib; these drugs are characterized by medium log p values, a parameter reflecting the overall hydrophobicity of the molecules, which is one important parameter to predict the interaction with lipid membranes. While all four molecules tend to embed in a similar region of the lipid membrane, their presence has different impacts on membrane structure and dynamics. Most notably, sunitinib, exhibiting the lowest log p value of the four inhibitors, effectively influences membrane integrity, while the others do not. This shows that the estimation of the effect of drug molecules on lipid membranes can be rather complex. In this context, experimental studies on lipid membranes are necessary to (i) identify drugs that may disturb membranes and (ii) characterize drug-membrane interactions on a molecular level. Such knowledge is important for understanding the efficacy and potential side effects of respective drugs.

The dynamic face of lipid membranes

Cell membranes - primarily composed of lipids, sterols, and proteins - form a dynamic interface between living cells and their environment. They act as a mechanical barrier around the cell while selectively facilitating material transport, signal transduction, and various other functions necessary for the cell viability. The complex functionality of cell membranes and the hierarchical motions and responses they exhibit demand a thorough understanding of the origin of different membrane dynamics and how they are influenced by molecular additives and environmental cues. These dynamic modes include single-molecule diffusion, thermal fluctuations, and large-scale membrane deformations, to name a few. This review highlights advances in investigating structure-driven dynamics associated with model cell membranes, with a particular focus on insights gained from neutron scattering and spectroscopy experiments. We discuss the uniqueness of neutron contrast variation and its remarkable potential in probing selective membrane structure and dynamics on spatial and temporal scales over which key biological functions occur. We also present a summary of current and future opportunities in synergistic combinations of neutron scattering with molecular dynamics (MD) simulations to gain further understanding of the molecular mechanisms underlying complex membrane functions.

Endoplasmic reticulum phospholipid scramblase activity revealed after protein reconstitution into giant unilamellar vesicles containing a photostable lipid reporter

Transbilayer movement of phospholipids in biological membranes is mediated by a diverse set of lipid transporters. Among them are scramblases that facilitate a rapid bi-directional movement of lipids without metabolic energy input. Here, we established a new fluorescence microscopy-based assay for detecting phospholipid scramblase activity of membrane proteins upon their reconstitution into giant unilamellar vesicles formed from proteoliposomes by electroformation. The assay is based on chemical bleaching of fluorescence of a photostable ATTO-dye labeled phospholipid with the membrane-impermeant reductant sodium dithionite. We demonstrate that this new methodology is suitable for the study of the scramblase activity of the yeast endoplasmic reticulum at single vesicle level.

Improved Protocol to Tackle the pH Effects on Membrane-Inserting Peptides

Many important biological pathways rely on membrane-interacting peptides or proteins, which can alter the biophysical properties of the cell membrane by simply adsorbing to its surface to undergo a full insertion process. To study these phenomena with atomistic detail, model peptides have been used to refine the current computational methodologies. Improvements have been made with force-field parameters, enhanced sampling techniques to obtain faster sampling, and the addition of chemical-physical properties, such as pH, whose influence dramatically increases at the water/membrane interface. The pH (low) insertion peptide (pHLIP) is a peptide that inserts across a membrane bilayer depending on the pH due to the presence of a key residue (Asp14) whose acidity-induced protonation triggers the whole process. The complex nature of these peptide/membrane interactions resulted in sampling limitations of the protonation and configurational space albeit using state-of-the-art methods such as the constant-pH molecular dynamics. To address this issue and circumvent those limitations, new simulations were performed with our newly developed pH-replica exchange method using wild-type (wt)-pHLIP in different 2-oleoyl-1-palmitoyl-sn-glycero-3-phosphocholine membrane sizes. This technique provided enhanced sampling and allowed for the calculation of more complete Asp14 pK_a profiles. The conformational heterogeneity derived from strong electrostatic interactions between Asp14 and the lipid phosphate groups was identified as the source of most pK_a variability. In spite of these persistent and harder-to-equilibrate phosphate interactions, the pK_a values at deeper regions (6.0-6.2) still predicted the experimental pK of insertion (6.0) since the electrostatic perturbation decays as the residue inserts further into the membrane. We also observed that reducing the system size leads to membrane deformations where it increasingly loses the ability to accommodate the pHLIP-induced perturbations. This indicates that large membrane patches, such as 256 or even 352 lipids, are needed to obtain stable and more realistic pHLIP/membrane systems. These results strengthen our method pK_a predictive and analytical capabilities to study the intricate play of electrostatic effects of the peptide/membrane interface, granting confidence for future applications in similar systems.

MLLPA: A Machine Learning-assisted Python module to study phase-specific events in lipid membranes

Machine Learning-assisted Lipid Phase Analysis (MLLPA) is a new Python 3 module developed to analyze phase domains in a lipid membrane based on lipid molecular states. Reading standard simulation coordinate and trajectory files, the software first analyze the phase composition of the lipid membrane by using machine learning tools to label each individual molecules with respect to their state, and then decompose the simulation box using Voronoi tessellations to analyze the local environment of all the molecules of interest. MLLPA is versatile as it can read from multiple format (e.g., GROMACS, LAMMPS) and from either all-atom (e.g., CHARMM36) or coarse-grain models (e.g., Martini). It can also analyze multiple geometries of membranes (e.g., bilayers, vesicles). Finally, the software allows for training with more than two phases, allowing for multiple phase coexistence analysis.

Altered Membrane Mechanics Provides a Receptor-Independent Pathway for Serotonin Action

Serotonin, an important signaling molecule in humans, has an unexpectedly high lipid membrane affinity. The significance of this finding has evoked considerable speculation. Here we show that membrane binding by serotonin can directly modulate membrane properties and cellular function, providing an activity pathway completely independent of serotonin receptors. Atomic force microscopy shows that serotonin makes artificial lipid bilayers softer, and induces nucleation of liquid disordered domains inside the raft-like liquid-ordered domains. Solid-state NMR spectroscopy corroborates this data at the atomic level, revealing a homogeneous decrease in the order parameter of the lipid chains in the presence of serotonin. In the RN46A immortalized serotonergic neuronal cell line, extracellular serotonin enhances transferrin receptor endocytosis, even in the presence of broad-spectrum serotonin receptor and transporter inhibitors. Similarly, it increases the membrane binding and internalization of oligomeric peptides. Our results uncover a mode of serotonin-membrane interaction that can potentiate key cellular processes in a receptor-independent fashion.

Understanding the influence of lipid bilayers and ligand molecules in determining the conformational dynamics of somatostatin receptor 2

Somatostatin receptor 2 (SSTR2) is a G-protein coupled receptor (GPCR) that controls numerous cellular processes including cell-to-cell signaling. In this study, we report how the lipid and ligand molecules influence the conformational dynamics of the membrane-bound SSTR2. Molecular simulations of different holo and apoenzyme complexes of SSTR2 in the presence and absence of a lipid bilayer were performed, observed, and correlated with previously reported studies. We identified the important SSTR2 residues that take part in the formation of the SSTR2-ligand complex. On analyzing the molecular simulation trajectories, we identified that the residue D3.32 is crucial in determining the bioactive conformation of SSTR2 ligands in the binding site. Based on the results, we suggest that designing a novel SSTR2 ligand with an H-bond donor group at the R1 position, and hydrophobic groups at R2 and R3 might have higher activity and SSTR2-selectivity. We analyzed the simulated systems to identify other important structural features involved in SSTR2-ligand binding and to observe the different conformational changes that occur in the protein after the ligand binding. Additionally, we studied the conformational dynamics of N- and C-terminal regions of SSTR2 in the presence and absence of the lipid bilayer. Both the systems were compared to understand the influence of lipid molecules in the formation of secondary structural domains by these extracellular regions. The comparative study revealed that the secondary structural elements formed by C-terminal residues in presence of lipid molecules is crucial for the functioning of SSTR2. Our study results highlight the structural complexities involved in the functioning of SSTR upon binding with the ligands in the presence and absence of lipid bilayer, which is essential for designing novel drug targets.

Understanding the influence of lipid bilayers and ligand molecules in determining the conformational dynamics of somatostatin receptor 2

Somatostatin receptor 2 (SSTR2) is a G-protein coupled receptor (GPCR) that controls numerous cellular processes including cell-to-cell signaling. In this study, we report how the lipid and ligand molecules influence the conformational dynamics of the membrane-bound SSTR2. Molecular simulations of different holo and apoenzyme complexes of SSTR2 in the presence and absence of a lipid bilayer were performed, observed, and correlated with previously reported studies. We identified the important SSTR2 residues that take part in the formation of the SSTR2-ligand complex. On analyzing the molecular simulation trajectories, we identified that the residue D3.32 is crucial in determining the bioactive conformation of SSTR2 ligands in the binding site. Based on the results, we suggest that designing a novel SSTR2 ligand with an H-bond donor group at the R1 position, and hydrophobic groups at R2 and R3 might have higher activity and SSTR2-selectivity. We analyzed the simulated systems to identify other important structural features involved in SSTR2-ligand binding and to observe the different conformational changes that occur in the protein after the ligand binding. Additionally, we studied the conformational dynamics of N- and C-terminal regions of SSTR2 in the presence and absence of the lipid bilayer. Both the systems were compared to understand the influence of lipid molecules in the formation of secondary structural domains by these extracellular regions. The comparative study revealed that the secondary structural elements formed by C-terminal residues in presence of lipid molecules is crucial for the functioning of SSTR2. Our study results highlight the structural complexities involved in the functioning of SSTR upon binding with the ligands in the presence and absence of lipid bilayer, which is essential for designing novel drug targets.

The Bacteriostatic Activity of 2-Phenylethanol Derivatives Correlates with Membrane Binding Affinity

The hydrophobic tails of aliphatic primary alcohols do insert into the hydrophobic core of a lipid bilayer. Thereby, they disrupt hydrophobic interactions between the lipid molecules, resulting in a decreased lipid order, i.e., an increased membrane fluidity. While aromatic alcohols, such as 2-phenylethanol, also insert into lipid bilayers and disturb the membrane organization, the impact of aromatic alcohols on the structure of biological membranes, as well as the potential physiological implication of membrane incorporation has only been studied to a limited extent. Although diverse targets are discussed to be causing the bacteriostatic and bactericidal activity of 2-phenylethanol, it is clear that 2-phenylethanol severely affects the structure of biomembranes, which has been linked to its bacteriostatic activity. Yet, in fungi some 2-phenylethanol derivatives are also produced, some of which appear to also have bacteriostatic activities. We showed that the 2-phenylethanol derivatives phenylacetic acid, phenyllactic acid, and methyl phenylacetate, but not Tyrosol, were fully incorporated into model membranes and affected the membrane organization. Furthermore, we observed that the propensity of the herein-analyzed molecules to partition into biomembranes positively correlated with their respective bacteriostatic activity, which clearly linked the bacteriotoxic activity of the substances to biomembranes.

Microscopic Interactions of Melatonin, Serotonin and Tryptophan with Zwitterionic Phospholipid Membranes

The interactions at the atomic level between small molecules and the main components of cellular plasma membranes are crucial for elucidating the mechanisms allowing for the entrance of such small species inside the cell. We have performed molecular dynamics and metadynamics simulations of tryptophan, serotonin, and melatonin at the interface of zwitterionic phospholipid bilayers. In this work, we will review recent computer simulation developments and report microscopic properties, such as the area per lipid and thickness of the membranes, atomic radial distribution functions, angular orientations, and free energy landscapes of small molecule binding to the membrane. Cholesterol affects the behaviour of the small molecules, which are mainly buried in the interfacial regions. We have observed a competition between the binding of small molecules to phospholipids and cholesterol through lipidic hydrogen-bonds. Free energy barriers that are associated to translational and orientational changes of melatonin have been found to be between 10-20 kJ/mol for distances of 1 nm between melatonin and the center of the membrane. Corresponding barriers for tryptophan and serotonin that are obtained from reversible work methods are of the order of 10 kJ/mol and reveal strong hydrogen bonding between such species and specific phospholipid sites. The diffusion of tryptophan and melatonin is of the order of 10-7 cm2/s for the cholesterol-free and cholesterol-rich setups.

Solid-state NMR spectroscopy

Solid-state nuclear magnetic resonance (NMR) spectroscopy is an atomic-level method used to determine the chemical structure, three-dimensional structure, and dynamics of solids and semi-solids. This Primer summarizes the basic principles of NMR as applied to the wide range of solid systems. The fundamental nuclear spin interactions and the effects of magnetic fields and radiofrequency pulses on nuclear spins are the same as in liquid-state NMR. However, because of the anisotropy of the interactions in the solid state, the majority of high-resolution solid-state NMR spectra is measured under magic-angle spinning (MAS), which has profound effects on the types of radiofrequency pulse sequences required to extract structural and dynamical information. We describe the most common MAS NMR experiments and data analysis approaches for investigating biological macromolecules, organic materials, and inorganic solids. Continuing development of sensitivity-enhancement approaches, including 1H-detected fast MAS experiments, dynamic nuclear polarization, and experiments tailored to ultrahigh magnetic fields, is described. We highlight recent applications of solid-state NMR to biological and materials chemistry. The Primer ends with a discussion of current limitations of NMR to study solids, and points to future avenues of development to further enhance the capabilities of this sophisticated spectroscopy for new applications.

Phase Separation in Lipid Lamellae Result from Ceramide Conformations and Lateral Packing Structure

Intercellular lipids in the stratum corneum protect the living body from invasion by allergens and pathogens, and also suppresses water evaporation within the body. It is important to understand how differences in the microstructure of intercellular lipids arise. This microstructure is affected by lipid composition. Studies using intercellular lipid models have reported the formation of two phases with different short lamellar periodicities. However, the details of the packing structure characteristics of the two phases observed in these intercellular lipid models are unclear. Our previous report revealed that different short periodicity phases coexist in the N-(α-hydroxyoctadecanoyl)-dihydrosphingosine (CER[ADS]), cholesterol (CHOL), and palmitic acid (PA) complex model. In this study, the characteristics of the packing structure of two phases with different short lamellar periodicities, which were observed in the intercellular lipid model (CER[ADS]/CHOL/PA) that we used previously, were adjusted for models with different lipid compositions. The characteristics of the packed and lamellar structures have been determined by temperature-scanning small-angle X-ray scattering and wide-angle X-ray diffraction measurements simultaneously. These differences in lamellar structure were thought to be caused by differences in ceramides (CER) conformation between the hairpin and the V-shape type. The lamellar structure of the V-shaped CER conformation has a low orthorhombic ratio. The above results suggest that an increase in the ratio of CER with the V-shaped structure causes the lamellar structure to have low orthorhombic ratio, thereby contributing to a decrease in the bilayer's barrier function.

Implementation of a self-consistent slab model of bilayer structure in the SasView suite

Slab models are simple and useful structural descriptions which have long been used to describe lyotropic lamellar phases, such as lipid bilayers. Typically, slab models assume a midline symmetry and break a bilayer structure into three pieces, a central solvent-free core and two symmetric outer layers composed of the soluble portion of the amphiphile and associated solvent. This breakdown matches reasonably well to the distribution of neutron scattering length density and therefore is a convenient and common approach for the treatment of small-angle scattering data. Here, an implementation of this model within the SasView software suite is reported. The implementation is intended to provide physical consistency through the area per amphiphile molecule and number of solvent molecules included within the solvent-exposed outer layer. The proper use of this model requires knowledge of (or good estimates for) the amphiphile and solvent molecule volume and atomic composition, ultimately providing a self-consistent data treatment with only two free parameters: the lateral area per amphiphile molecule and the number of solvent molecules included in the outer region per amphiphile molecule. The use of this code is demonstrated in the fitting of standard lipid bilayer data sets, obtaining structural parameters consistent with prior literature and illustrating the typical and ideal cases of fitting for neutron scattering data obtained using single or multiple contrast conditions. While demonstrated here for lipid bilayers, this model is intended for general application to block copolymers, surfactants, and other lyotropic lamellar phase structures for which a slab model is able to reasonably estimate the neutron scattering length density/electron-density profile of inner and outer layers of the lamellae.

Interplay of cholesterol, membrane bilayers and the AT1R: A cholesterol consensus motif on AT1R is revealed

Hypertension, mediated by the Angiotensin II receptor type 1 (AT1R), is still the major cause of premature death despite the discovery of novel therapeutics, highlighting the importance of an in depth understanding of the drug-AT1R recognition mechanisms coupled with the impact of the membrane environment on the interaction of drugs with AT1R. Herein, we examine the interplay of cholesterol-lipid-candesartan and the AT1R using Molecular Dynamics simulations of a model membrane consisting of 60:40 mol%. DPPC:cholesterol, candesartan and the AT1R, mimicking the physiological cholesterol concentration in sarcolemma membranes. The simulations of the model membrane of 60:40 mol%. DPPC:cholesterol were further validated using DOSY NMR experiments. Interestingly, our results suggest a significant role of cholesterol in the AT1R function imposed through a Cholesterol Consensus Motif (CCM) in the receptor, which could be crucial in the drug binding process. Candesartan diffusion towards AT1R through incorporation into lipid bilayers, appears to be retarded by the presence of cholesterol. However, its direct approach towards AT1R may be facilitated through the mobility induced on the N-terminus by the cholesterol binding on the CCM these novel insights could pave the way towards the development of more potent pharmaceutical agents to combat hypertension more effectively.

Binding of the small-molecule kinase inhibitor ruxolitinib to membranes does not disturb membrane integrity

Ruxolitinib is a small-molecule protein kinase inhibitor, which is used as a therapeutic agent against several diseases. Due to its anti-inflammatory impact, ruxolitinib has also been considered recently for usage in the treatment of Covid-19. While the specific effects of ruxolitinib on Janus kinases (JAK) is comparatively well investigated, its (unspecific) impact on membranes has not been studied in detail so far. Therefore, we characterized the interaction of this drug with lipid membranes employing different biophysical approaches. Ruxolitinib incorporates into the glycerol region of lipid membranes causing an increase in disorder of the lipid chains. This binding, however, has only marginal influence on the structure and integrity of membranes as found by leakage and permeation assays.

•

The drug ruxolitinib binds to lipid membranes.

•

Ruxolitinib intercalates into the glycerol region of the bilayer.

•

The membrane binding of ruxolitinib causes an increased disorder of the lipid chains.

•

The membrane binding of ruxolitinib has no influence on the membrane integrity.

Serotonin Alters the Phase Equilibrium of a Ternary Mixture of Phospholipids and Cholesterol

Unsaturated and saturated phospholipids tend to laterally segregate, especially in the presence of cholesterol. Small molecules such as neurotransmitters, toxins, drugs etc. possibly modulate this lateral segregation. The small aromatic neurotransmitter serotonin (5-HT) has been found to bind to membranes. We studied the lipid structure and packing of a ternary membrane mixture consisting of palmitoyl-oleoyl-phosphatidylcholine, palmitoyl-sphingomyelin, and cholesterol at a molar ratio of 4/4/2 in the absence and in the presence of 5-HT, using a combination of solid-state 2H NMR, atomic force microscopy, and atomistic molecular dynamics (MD) simulations. Both NMR and MD report formation of a liquid ordered (L o ) and a liquid disordered (L d ) phase coexistence with small domains. Lipid exchange between the domains was fast such that single component 2H NMR spectra are detected over a wide temperature range. A drastic restructuring of the domains was induced when 5-HT is added to the membranes at a 9 mol% concentration relative to the lipids. 2H NMR spectra of all components of the mixture showed two prominent contributions indicative of molecules of the same kind residing both in the disordered and the ordered phase. Compared to the data in the absence of 5-HT, the lipid chain order in the disordered phase was further decreased in the presence of 5-HT. Likewise, addition of serotonin increased lipid chain order within the ordered phase. These characteristic lipid chain order changes were confirmed by MD simulations. The 5-HT-induced larger difference in lipid chain order results in more pronounced differences in the hydrophobic thickness of the individual membrane domains. The correspondingly enlarged hydrophobic mismatch between ordered and disordered phases is assumed to increase the line tension at the domain boundary, which drives the system into formation of larger size domains. These results not only demonstrate that small membrane binding molecules such as neurotransmitters have a profound impact on essential membrane properties. It also suggests a mechanism by which the interaction of small molecules with membranes can influence the function of membrane proteins and non-cognate receptors. Altered membrane properties may modify lateral sorting of membrane protein, membrane protein conformation, and thus influence their function as suspected for neurotransmitters, local anesthetics, and other small drug molecules.

Modulation of Membrane Fluidity to Control Interfacial Water Structure and Dynamics in Saturated and Unsaturated Phospholipid Vesicles

The structure and dynamics of interfacial water in biological systems regulate the biochemical reactions. But, it is still enigmatic how the behavior of the interfacial water molecule is controlled. Here, we have investigated the effect of membrane fluidity on the structure and dynamics of interfacial water molecules in biologically relevant phopholipid vesicles. This study delineates that modulation of membrane fluidity through interlipid separation and unsaturation not only mitigate membrane rigidity but also disrupt the strong hydrogen bond (H-bond) network around the lipid bilayer interface. As a result, a disorder in H-bonding between water molecules arises several layers beyond the first hydration shell of the polar headgroup, which essentially modifies the interfacial water structure and dynamics. Furthermore, we have also provided evidence of increasing transportation through these modulated membranes, which enhance the membrane mediated isomerization reaction rate.

Conformation and Orientation of Branched Acyl Chains Responsible for the Physical Stability of Diphytanoylphosphatidylcholine

Diphytanoylphosphatidylcholine (DPhPC) is a synthetic phospholipid in which two methyl-branched acyl chains are introduced into the glycerol moiety, mimicking phospholipids of eukaryotic and eubacterial origins. The lipid bilayers of DPhPC reproduce the outstanding physical properties of methyl-branched lipids that occur in archaeal membranes. DPhPC is commonly used as the base lipid in biophysical experiments, particularly for recording ion-channel currents. However, the dynamics of lipid molecules that induces their useful physical properties is still unclear. In this study, we examined the conformation and orientation of the methyl-branched acyl chain of DPhPC in a membrane using ²H nuclear magnetic resonance (NMR) measurements of the synthetic lipid with a high stereochemical purity and molecular dynamics (MD) simulations. Deuterium-labeled 3',3'-CD₃,D-DPhPC (2) and 7',7'-CD₃,D-DPhPC (3) showed the characteristic quadrupole splitting width in the ²H NMR spectra, which corresponded to the bent orientation reported for the archaeal lipid PGP-Me [Yamagami, M., et al. (2019) Biochemistry 58, 3869-3879]. However, MD simulations, which reproduced the ²H NMR results well, unveiled the unknown features of DPhPC in the membrane; DPhPC has a chain-specific average orientation, where two bent orientations with upward and downward methyl groups occur at positions C3 and C7 of the sn-1 and sn-2 chains of DPhPC, respectively. These MD and NMR results reveal that these two bent orientations define the average orientation of DPhPC for the shallow part of the acyl chains, which is considered to be an important factor in the stability of DPhPC membranes.

Microsecond molecular dynamics simulation of the adsorption and penetration of oil droplets on cellular membrane

The hazardous effects of petroleum contaminants in the soil and water environment are highly associated with their interactions with cellular membranes, but our understanding on the molecular-level mechanisms for the adsorption and penetration of heavy oil mixture on cellular membrane is very limited. In this study, microsecond molecular dynamics simulations were performed to gain insights into the morphological evolution and penetration dynamics of the multi-component and single-component oil droplets on the dipalmitoylphosphatidylcholine lipid membrane. Results highlighted the inhibition effect of the resins on the penetration of alkanes and aromatics, because they would form net structure making it difficult to release the latter two components from the oil droplet to the membrane. It also demonstrated the obviously different patterns of penetration between alkanes and aromatics. The overall steps for the toluene penetration included detachment from oil droplet, dispersion in water, adsorption on membrane surface, structure adjustment and penetration into membrane. By contrast, the step of dispersion in water was not necessary for the alkanes' penetration. Instead, it relied on the adsorption of the whole oil droplet on the membrane surface which resulted in the formation of pores on the membrane surface by local structure deformation in the lipid head group regions.

Lateral heterogeneity and domain formation in cellular membranes

As early as the development of the fluid mosaic model for cellular membranes, researchers began observing the telltale signs of lateral heterogeneity. Over the decades this has led to the development of the lipid raft hypothesis and the ensuing controversy that has unfolded, as a result. Here, we review the physical concepts behind domain formation in lipid membranes, both of their structural and dynamic origins. This, then leads into a discussion of coarse-grained, phenomenological approaches that describe the wide range of phases associated with lipid lateral heterogeneity. We use these physical concepts to describe the interaction between raft-lipid species, such as long-chain saturated lipids, sphingomyelin, and cholesterol, and non-raft forming lipids, such as those with short acyl chains or unsaturated fatty acids. While debate has persisted on the biological relevance of lipid domains, recent research, described here, continues to identify biological roles for rafts and new experimental approaches have revealed the existence of lipid domains in living systems. Given the recent progress on both the biological and structural aspects of raft formation, the research area of membrane lateral heterogeneity will not only expand, but will continue to produce exciting results.

The Sphingosine and Acyl Chains of Ceramide [NS] Show Very Different Structure and Dynamics That Challenge Our Understanding of the Skin Barrier

The lipid phase of the uppermost human skin layer is thought to comprise highly rigid lipids in an orthorhombic phase state to protect the body against the environment. By synthesizing sphingosine-d28 deuterated N-lignoceroyl-d-erythro-sphingosine (ceramide [NS]), we compare the structure and dynamics of both chains of that lipid in biologically relevant mixtures using X-ray diffraction, 2 H NMR analysis, and infrared spectroscopy. Our results reveal a substantial fraction of sphingosine chains in a fluid and dynamic phase state at physiological temperature. These findings prompt revision of our current understanding of the skin lipid barrier, where an extended ceramide [NS] conformation is preferred and a possible domain structure is proposed. Mobile lipid chains may be crucial for skin elasticity and the translocation of physiologically important molecules.

Membrane Interaction of Ibuprofen with Cholesterol-Containing Lipid Membranes

Deciphering the membrane interaction of drug molecules is important for improving drug delivery, cellular uptake, and the understanding of side effects of a given drug molecule. For the anti-inflammatory drug ibuprofen, several studies reported contradictory results regarding the impact of ibuprofen on cholesterol-containing lipid membranes. Here, we investigated membrane localization and orientation as well as the influence of ibuprofen on membrane properties in POPC/cholesterol bilayers using solid-state NMR spectroscopy and other biophysical assays. The presence of ibuprofen disturbs the molecular order of phospholipids as shown by alterations of the ²H and ³¹P-NMR spectra of the lipids, but does not lead to an increased membrane permeability or changes of the phase state of the bilayer. ¹H MAS NOESY NMR results demonstrate that ibuprofen adopts a mean position in the upper chain/glycerol region of the POPC membrane, oriented with its polar carbonyl group towards the aqueous phase. This membrane position is only marginally altered in the presence of cholesterol. A previously reported result that ibuprofen is expelled from the membrane interface in cholesterol-containing DMPC bilayers could not be confirmed.

Optimization of Slipids Force Field Parameters Describing Headgroups of Phospholipids

The molecular mechanics force field Slipids developed in a series of works by Jämbeck and Lyubartsev (J. Phys. Chem. B 2012, 116, 3164-3179; J. Chem. Theory Comput. 2012, 8, 2938-2948) generally provides a good description of various lipid bilayer systems. However, it was also found that order parameters of C-H bonds in the glycerol moiety of the phosphatidylcholine headgroup deviate significantly from NMR results. In this work, the dihedral force field parameters have been reparameterized in order to improve the agreement with experiment. For this purpose, we have computed energies for a large amount of lipid headgroup conformations using density functional theory on the B3P86/cc-pvqz level and optimized dihedral angle parameters simultaneously to provide the best fit to the quantum chemical energies. The new parameter set was validated for three lipid bilayer systems against a number of experimental properties including order parameters, area per lipid, scattering form factors, bilayer thickness, area compressibility and lateral diffusion coefficients. In addition, the order parameter dependence on cholesterol content in the POPC bilayer was investigated. It is shown that the new force field significantly improves agreement with the experimental order parameters for the lipid headgroup while keeping good agreement with other experimentally measured properties.

Scaling relationships for the elastic moduli and viscosity of mixed lipid membranes

The elastic and viscous properties of biological membranes play a vital role in controlling cell functions that require local reorganization of the membrane components as well as dramatic shape changes such as endocytosis, vesicular trafficking, and cell division. These properties are widely acknowledged to depend on the unique composition of lipids within the membrane, yet the effects of lipid mixing on the membrane biophysical properties remain poorly understood. Here, we present a comprehensive characterization of the structural, elastic, and viscous properties of fluid membranes composed of binary mixtures of lipids with different tail lengths. We show that the mixed lipid membrane properties are not simply additive quantities of the single-component analogs. Instead, the mixed membranes are more dynamic than either of their constituents, quantified as a decrease in their bending modulus, area compressibility modulus, and viscosity. While the enhanced dynamics are seemingly unexpected, we show that the measured moduli and viscosity for both the mixed and single-component bilayers all scale with the area per lipid and collapse onto respective master curves. This scaling links the increase in dynamics to mixing-induced changes in the lipid packing and membrane structure. More importantly, the results show that the membrane properties can be manipulated through lipid composition the same way bimodal blends of surfactants, liquid crystals, and polymers are used to engineer the mechanical properties of soft materials, with broad implications for understanding how lipid diversity relates to biomembrane function.

Protein crowding in the inner mitochondrial membrane

The inner membrane of mitochondria is known for its low lipid-to-protein ratio. Calculations based on the size and the concentration of the principal membrane components, suggest about half of the hydrophobic volume of the membrane is occupied by proteins. Such high degree of crowding is expected to strain the hydrophobic coupling between proteins and lipids unless stabilizing mechanisms are in place. Both protein supercomplexes and cardiolipin are likely to be critical for the integrity of the inner mitochondrial membrane because they reduce the energy penalty of crowding.

Molecular-Level Understanding of the Somatostatin Receptor 1 (SSTR1)-Ligand Binding: A Structural Biology Study Based on Computational Methods

Somatostatin receptor 1 (SSTR1), a subtype of somatostatin receptors, is involved in various signaling mechanisms in different parts of the human body. Like most of the G-protein-coupled receptors (GPCRs), the available information on the structural features of SSTR1 responsible for the biological activity is scarce. In this study, we report a molecular-level understanding of SSTR1-ligand binding, which could be helpful in solving the structural complexities involved in SSTR1 functioning. Based on a three-dimensional quantitative structure-activity relationship (3D-QSAR) study using comparative molecular field analysis (CoMFA) and comparative molecular similarity index analysis (CoMSIA), we have identified that an electronegative, less-bulkier, and hydrophobic atom substitution can substantially increase the biological activity of SSTR1 ligands. A density functional theory (DFT) study has been followed to study the electron-related properties of the SSTR1 ligands and to validate the results obtained via the 3D-QSAR study. 3D models of SSTR1-ligand systems have been embedded in lipid-lipid bilayer membranes to perform molecular dynamics (MD) simulations. Analysis of the MD trajectories reveals important information about the crucial residues involved in SSTR1-ligand binding and various conformational changes in the protein that occur after ligand binding. Additionally, we have identified the probable ligand-binding site of SSTR1 and validated it using MD. We have also studied the favorable conditions that are essential for forming the most stable and lowest-energy bioactive conformation of the ligands inside the binding site. The results of the study could be useful in constructing more potent and novel SSTR1 antagonists and agonists.

How cholesterol stiffens unsaturated lipid membranes

Cholesterol is an integral component of eukaryotic cell membranes and a key molecule in controlling membrane fluidity, organization, and other physicochemical parameters. It also plays a regulatory function in antibiotic drug resistance and the immune response of cells against viruses, by stabilizing the membrane against structural damage. While it is well understood that, structurally, cholesterol exhibits a densification effect on fluid lipid membranes, its effects on membrane bending rigidity are assumed to be nonuniversal; i.e., cholesterol stiffens saturated lipid membranes, but has no stiffening effect on membranes populated by unsaturated lipids, such as 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC). This observation presents a clear challenge to structure-property relationships and to our understanding of cholesterol-mediated biological functions. Here, using a comprehensive approach-combining neutron spin-echo (NSE) spectroscopy, solid-state deuterium NMR (²H NMR) spectroscopy, and molecular dynamics (MD) simulations-we report that cholesterol locally increases the bending rigidity of DOPC membranes, similar to saturated membranes, by increasing the bilayer's packing density. All three techniques, inherently sensitive to mesoscale bending fluctuations, show up to a threefold increase in effective bending rigidity with increasing cholesterol content approaching a mole fraction of 50%. Our observations are in good agreement with the known effects of cholesterol on the area-compressibility modulus and membrane structure, reaffirming membrane structure-property relationships. The current findings point to a scale-dependent manifestation of membrane properties, highlighting the need to reassess cholesterol's role in controlling membrane bending rigidity over mesoscopic length and time scales of important biological functions, such as viral budding and lipid-protein interactions.

Role of the lipid bilayer in outer membrane protein folding in Gram-negative bacteria

β-Barrel outer membrane proteins (OMPs) represent the major proteinaceous component of the outer membrane (OM) of Gram-negative bacteria. These proteins perform key roles in cell structure and morphology, nutrient acquisition, colonization and invasion, and protection against external toxic threats such as antibiotics. To become functional, OMPs must fold and insert into a crowded and asymmetric OM that lacks much freely accessible lipid. This feat is accomplished in the absence of an external energy source and is thought to be driven by the high thermodynamic stability of folded OMPs in the OM. With such a stable fold, the challenge that bacteria face in assembling OMPs into the OM is how to overcome the initial energy barrier of membrane insertion. In this review, we highlight the roles of the lipid environment and the OM in modulating the OMP-folding landscape and discuss the factors that guide folding in vitro and in vivo We particularly focus on the composition, architecture, and physical properties of the OM and how an understanding of the folding properties of OMPs in vitro can help explain the challenges they encounter during folding in vivo Current models of OMP biogenesis in the cellular environment are still in flux, but the stakes for improving the accuracy of these models are high. OMP folding is an essential process in all Gram-negative bacteria, and considering the looming crisis of widespread microbial drug resistance it is an attractive target. To bring down this vital OMP-supported barrier to antibiotics, we must first understand how bacterial cells build it.

Interaction of the small-molecule kinase inhibitors tofacitinib and lapatinib with membranes

Lapatinib and tofacitinib are small-molecule kinase inhibitors approved for the treatment of advanced or metastatic breast cancer and rheumatoid arthritis, respectively. So far, the mechanisms which are responsible for their activities are not entirely understood. Here, we focus on the interaction of these drug molecules with phospholipid membranes, which has not yet been investigated before in molecular detail. Owing to their lipophilic characteristics, quantitatively reflected by large differences of the partition equilibrium between water and octanol phases (expressed by logP values), rather drastic differences in the membrane interaction of both molecules have to be expected. Applying experimental (nuclear magnetic resonance, fluorescence and ESR spectroscopy) and theoretical (molecular dynamics simulations) approaches, we found that lapatinib and tofacitinib bind to lipid membranes and insert into the lipid-water interface of the bilayer. For lapatinib, a deeper embedding into the membrane bilayer was observed than for tofacitinib implying different impacts of the molecules on the bilayer structure. While for tofacitinib, no influence to the membrane structure was found, lapatinib causes a membrane disturbance, as concluded from an increased permeability of the membrane for polar molecules. These data may contribute to a better understanding of the cellular uptake mechanism(s) and the side effects of the drugs.

Atomistic Molecular Dynamics Simulations of Propofol and Fentanyl in Phosphatidylcholine Lipid Bilayers

Atomistic molecular dynamics (MD) and steered MD simulations in combination with umbrella sampling methodology were utilized to study the general anesthetic propofol and the opioid analgesic fentanyl and their interaction with lipid bilayers, which is not yet fully understood. These molecules were inserted into two different fully hydrated phospholipid bilayers, namely, dioleoylphosphatidylcholine (DOPC) and dipalmitoylphosphatidylcholine (DPPC), to investigate the effects that these drugs have on the bilayer. We determined the role of the lipid chain length and saturation on the behavior of the two drugs. Pure, fully hydrated DOPC and DPPC bilayers were also simulated, and the results were in excellent agreement with the experimental values. Various structural and mechanical properties of each system, such as the area per lipid, area compressibility modulus, order parameter, lateral lipid diffusion, hydrogen bonds, and radial distribution functions, have been calculated to assess how the drug molecules affect the different bilayers. From the calculated results, we show that fentanyl and propofol generally follow similar trends in each bilayer but adopt different favorable positions close to the headgroup/chain interface at the carbonyl groups. Propofol was shown to selectively form hydrogen bonds at the carbonyl carbon in each bilayer, whereas fentanyl interacts with water molecules at the headgroup interface. From the calculated free-energy profiles, we determined that both molecules show a preference for the low-density, low-order acyl chain region of the bilayers and both significantly preferred the DOPC bilayer with propofol and fentanyl having energy minima at -6.66 and -43.07 kcal mol^-1, respectively. This study suggests that different chain lengths and levels of saturation directly affect the properties of these two important molecules, which are seen to work together to control anesthesia in surgical applications.

Structural Properties of Inverted Hexagonal Phase: A Hybrid Computational and Experimental Approach

Inverted/reverse hexagonal (H_II) phases are of special interest in several fields of research, including nanomedicine. We used molecular dynamics (MD) simulation to study H_II systems composed of dioleoylphosphatidylethanolamine (DOPE) and palmitoyloleoylphosphatidylethanolamine (POPE) at several hydration levels and temperatures. The effect of the hydration level on several H_II structural parameters, including deuterium order parameters, was investigated. We further used MD simulations to estimate the maximum hydrations of DOPE and POPE H_II lattices at several given temperatures. Finally, the effect of acyl chain unsaturation degree on the H_II structure was studied via comparing the DOPE with POPE H_II systems. In addition to MD simulations, we used deuterium nuclear magnetic resonance (²H NMR) and small-angle X-ray scattering (SAXS) experiments to measure the DOPE acyl chain order parameters, lattice plane distances, and the water core radius in H_II phase DOPE samples at several temperatures in the presence of excess water. Structural parameters calculated from MD simulations are in excellent agreement with the experimental data. Dehydration decreases the radius of the water core. An increase in hydration level slightly increased the deuterium order parameter of lipids acyl chains, whereas an increase in temperature decreased it. Lipid cylinders undulated along the cylinder axis as a function of hydration level. The maximum hydration levels of PE H_II phases at different temperatures were successfully predicted by MD simulations based on a single experimental measurement for the lattice plane distance in the presence of excess water. An increase in temperature decreases the maximum hydration and consequently the radius of the water core and lattice plane distances. Finally, DOPE formed H_II structures with a higher curvature compared to POPE, as expected. We propose a general protocol for constructing computational H_II systems that correspond to the experimental systems. This protocol could be used to study H_II systems composed of molecules other than the PE systems used here and to improve and validate force field parameters by using the target data in the H_II phase.

Ultrafast Dynamics at the Lipid-Water Interface: DMSO Modulates H-Bond Lifetimes

Dimethyl sulfoxide (DMSO) is a common cosolvent and cryopreservation agent used to freeze cells and tissues. DMSO alters the H-bond structure of water, but its interactions with biomolecules and, specifically, with biological interfaces remain poorly understood. Here we investigate the effects of DMSO on the H-bond dynamics at the lipid-water interface using a combination of ultrafast two-dimensional infrared (2D IR) spectroscopy and molecular dynamics simulations. Ester carbonyl absorption spectra show that DMSO dehydrates the interface, and simulations show that the area per lipid is decreased. Ultrafast 2D IR spectra measure the time scales of frequency fluctuations at the ester carbonyl positions located precisely between the hydrophobic and hydrophilic regions of the membrane. 2D IR measurements show that low DMSO concentrations (<10 mol %) induce ∼40% faster H-bond dynamics compared with pure water, whereas increased concentrations (>10-20 mol %) once again slow down the dynamics. This slow-fast-slow trend is described in terms of two different solvation regimes. Below 10 mol %, DMSO weakens the interfacial H bond, leading to faster "bulk-like" dynamics, whereas above 10 mol %, water molecules become "relatively immobilized" as the H-bond networks becoming disrupted by the H-bond donor/acceptor imbalance at the interface. These studies are an important step toward characterizing the environments around lipid membranes, which are essential to numerous biological processes.

Gold nanoparticles interacting with synthetic lipid rafts: an AFM investigation

Inorganic nanoparticles (NPs) represent promising examples of engineered nanomaterials, providing interesting biomedical solutions in several fields, like therapeutics and diagnostics. Despite the extensive number of investigations motivated by their remarkable potential for nanomedicinal applications, the interactions of NPs with biological interfaces are still poorly understood. The effect of NPs on living organisms is mediated by biological barriers, such as the cell plasma membrane, whose lateral heterogeneity is thought to play a prominent role in NPs adsorption and uptake pathways. In particular, biological membranes feature the presence of rafts, that is segregated lipid micro and/or nanodomains in the so-called liquid ordered phase (L_o ), immiscible with the surrounding liquid disordered phase (L_d ). Rafts are involved in various biological functions and act as sites for the selective adsorption of materials on the membrane. Indeed, the thickness mismatch present along their boundaries generates energetically favourable conditions for the adsorption of NPs. Despite its clear implications in NPs internalisation processes and cytotoxicity, a direct proof of the selective adsorption of NPs along the rafts' boundaries is still missing to date. Here we use multicomponent supported lipid bilayers (SLBs) as reliable synthetic models, reproducing the nanometric lateral heterogeneity of cell membranes. After being characterised by atomic force microscopy (AFM) and neutron reflectivity (NR), multidomain SLBs are challenged by prototypical inorganic nanoparticles, that is citrated gold nanoparticles (AuNPs), under simplified and highly controlled conditions. By exploiting AFM, we demonstrate that AuNPs preferentially target lipid phase boundaries as adsorption sites. The herein reported study consolidates and extends the fundamental knowledge on NPs-membrane interactions, which constitute a key aspect to consider when designing NPs-related biomedical applications. LAY DESCRIPTION: Inorganic nanoparticles (NPs) represent promising examples of engineered nanomaterials, offering interesting biomedical solutions in multiple fields like therapeutics and diagnostics. Despite being extensively investigated due to their remarkable potential for nanomedicinal applications, the interaction of NPs with biological systems is in several cases still poorly understood. The interaction of NPs with living organisms is mediated by biological barriers, such as the cell plasma membrane. Supported lipid bilayers (SLBs) represent suitable synthetic membrane models for studying the physicochemical properties of natural interfaces and their interaction with inorganic nanomaterials under simplified and controlled conditions. Recently, multicomponent SLBs were developed in order to mimic the lateral heterogeneity of most biological membranes. In particular, biological membranes feature the presence of rafts, that is segregated lipid micro and/or nanodomains, enriched in cholesterol, sphingomyelin, saturated glycerophospholipids and glycosphingolipids: these lipids segregate in the so-called liquid-ordered phase (L_o ), characterised by a high molecular packing degree, which promotes the phase separation from the surrounding liquid-crystalline (disordered, L_d ) phase, where the intermolecular mobility is increased. Rafts are thought to participate in the formation and targeting of nano-sized biogenic lipid vesicles and are also actively involved in multiple membrane processes. Indeed, L_o -L_d phase boundaries represent high energy areas, providing active sites for the preferential adsorption of external material. The selective adsorption of NPs along the phase boundaries of rafted membranes has been theorised and indirectly probed by different research groups; however, a direct proof of this phenomenon is still missing to date. We herein exploit atomic force microscopy (AFM) to directly visualise the preferential adsorption of gold nanoparticles (AuNPs) along the phase boundaries of multicomponent SLBs (previously characterised by neutron reflectivity), obtained from synthetic vesicles containing both an L_d and an L_o phase. The quantitative localisation and morphometry of AuNPs adsorbed on the SLB reveal important information on their interaction with the lipid matrix and directly prove the already theorised differential NPs-lipid interaction at the phase boundaries. The presented results could help the development of future NP-based applications, involving NPs adsorption on membranes with nanoscale phase segregations.

A New Lipid Force Field (FUJI)

To explore inhomogeneous and anisotropic systems such as lipid bilayers, the Lennard-Jones particle mesh Ewald (LJ-PME) method has been applied without a conventional isotropic dispersion correction. As the popular AMBER and CHARMM lipid force fields were developed using a cutoff scheme, their lipid bilayers unacceptably shrink when using the LJ-PME method. In this study, a new all-atom lipid force field (FUJI) was developed on the basis of the AMBER force-field scheme including the Lipid14 van der Waals parameters. Point charges were calculated using the restrained electrostatic potentials of many lipid conformers. Further, torsion energy profiles were calculated using high-level ab initio molecular orbitals (LCCSD(T)/Aug-cc-pVTZ//LMP2/Aug-cc-pVTZ), following which the molecular mechanical dihedral parameters were derived through a fast Fourier transform. By incorporation of these parameters into a new lipid force field without fitting experimental data, the desired lipid characteristics such as the area per lipid and lateral diffusion coefficients were obtained through GROMACS molecular dynamics simulations using the LJ-PME method and virtual hydrogen sites. The calculated area per lipid and lateral diffusion coefficients showed satisfactory agreement with experimental data. Furthermore, the electron-density profiles along the membrane normal were calculated for pure lipid bilayers, and the resulting membrane thicknesses agreed well with the experimental values. As the new lipid force field is compatible with FUJI for protein and small molecules, the new FUJI force field will offer accurate modeling for complex systems consisting of various membrane proteins and lipids.

Carotenoids promote lateral packing and condensation of lipid membranes

Carotenoids are pigment molecules that protect biomembranes against degradation and may be involved in the formation of functional bacterial membrane microdomains. Little is known on whether different types of carotenoids have different effects on the membrane or if there is any concentration dependence of these effects. In this work, we present results from molecular dynamics simulations of phospholipid bilayers containing different amounts of either β-carotene or zeaxanthin. Both β-carotene and zeaxanthin show the ability to laterally condense the membrane lipids and reduce their inter-leaflet interactions. With increasing concentrations, both carotenoids increase the bilayer thickness and rigidity. The results reveal that carotenoids have similar effects to cholesterol on regulating the behavior of fluid-phase membranes, suggesting that they could function as sterol substitutes and confirming their potential role in the formation of functional membrane domains.

Hybrid particle-field molecular dynamics under constant pressure

Hybrid particle-field methods are computationally efficient approaches for modeling soft matter systems. So far, applications of these methodologies have been limited to constant volume conditions. Here, we reformulate particle-field interactions to represent systems coupled to constant external pressure. First, we show that the commonly used particle-field energy functional can be modified to model and parameterize the isotropic contributions to the pressure tensor without interfering with the microscopic forces on the particles. Second, we employ a square gradient particle-field interaction term to model non-isotropic contributions to the pressure tensor, such as in surface tension phenomena. This formulation is implemented within the hybrid particle-field molecular dynamics approach and is tested on a series of model systems. Simulations of a homogeneous water box demonstrate that it is possible to parameterize the equation of state to reproduce any target density for a given external pressure. Moreover, the same parameterization is transferable to systems of similar coarse-grained mapping resolution. Finally, we evaluate the feasibility of the proposed approach on coarse-grained models of phospholipids, finding that the term between water and the lipid hydrocarbon tails is alone sufficient to reproduce the experimental area per lipid in constant-pressure simulations and to produce a qualitatively correct lateral pressure profile.