Permeability of membranes in the liquid ordered and liquid disordered phases

PyRETIS 3: Conquering rare and slow events without boundaries

We present and discuss the advancements made in PyRETIS 3, the third instalment of our Python library for an efficient and user-friendly rare event simulation, focused to execute molecular simulations with replica exchange transition interface sampling (RETIS) and its variations. Apart from a general rewiring of the internal code towards a more modular structure, several recently developed sampling strategies have been implemented. These include recently developed Monte Carlo moves to increase path decorrelation and convergence rate, and new ensemble definitions to handle the challenges of long-lived metastable states and transitions with unbounded reactant and product states. Additionally, the post-analysis software PyVisa is now embedded in the main code, allowing fast use of machine-learning algorithms for clustering and visualising collective variables in the simulation data.

Structure and dynamics of cholesterol-mediated aquaporin-0 arrays and implications for lipid rafts

Aquaporin-0 (AQP0) tetramers form square arrays in lens membranes through a yet unknown mechanism, but lens membranes are enriched in sphingomyelin and cholesterol. Here, we determined electron crystallographic structures of AQP0 in sphingomyelin/cholesterol membranes and performed molecular dynamics (MD) simulations to establish that the observed cholesterol positions represent those seen around an isolated AQP0 tetramer and that the AQP0 tetramer largely defines the location and orientation of most of its associated cholesterol molecules. At a high concentration, cholesterol increases the hydrophobic thickness of the annular lipid shell around AQP0 tetramers, which may thus cluster to mitigate the resulting hydrophobic mismatch. Moreover, neighboring AQP0 tetramers sandwich a cholesterol deep in the center of the membrane. MD simulations show that the association of two AQP0 tetramers is necessary to maintain the deep cholesterol in its position and that the deep cholesterol increases the force required to laterally detach two AQP0 tetramers, not only due to protein-protein contacts but also due to increased lipid-protein complementarity. Since each tetramer interacts with four such 'glue' cholesterols, avidity effects may stabilize larger arrays. The principles proposed to drive AQP0 array formation could also underlie protein clustering in lipid rafts.

The antimicrobial fibupeptide lugdunin forms water-filled channel structures in lipid membranes

Recently, a novel cyclo-heptapeptide composed of alternating D,L-amino acids and a unique thiazolidine heterocycle, called lugdunin, was discovered, which is produced by the nasal and skin commensal Staphylococcus lugdunensis. Lugdunin displays potent antimicrobial activity against a broad spectrum of Gram-positive bacteria, including challenging-to-treat methicillin-resistant Staphylococcus aureus (MRSA). Lugdunin specifically inhibits target bacteria by dissipating their membrane potential. However, the precise mode of action of this new class of fibupeptides remains largely elusive. Here, we disclose the mechanism by which lugdunin rapidly destabilizes the bacterial membrane potential using an in vitro approach. The peptide strongly partitions into lipid compositions resembling Gram-positive bacterial membranes but less in those harboring the eukaryotic membrane component cholesterol. Upon insertion, lugdunin forms hydrogen-bonded antiparallel β-sheets by the formation of peptide nanotubes, as demonstrated by ATR-FTIR spectroscopy and molecular dynamics simulations. These hydrophilic nanotubes filled with a water wire facilitate not only the translocation of protons but also of monovalent cations as demonstrated by voltage-clamp experiments on black lipid membranes. Collectively, our results provide evidence that the natural fibupeptide lugdunin acts as a peptidic channel that is spontaneously formed by an intricate stacking mechanism, leading to the dissipation of a bacterial cell's membrane potential.

Spatial Arrangement of the Drug Ibuprofen in a Model Membrane in the Presence of Lipid Rafts

Many pharmaceutical drugs are known to interact with lipid membranes through nonspecific molecular interactions, which affect their therapeutic effect. Ibuprofen is a nonsteroidal anti-inflammatory drug (NSAID) and one of the most commonly prescribed. In the presence of cholesterol, lipid bilayers can separate into nanoscale liquid-disordered and liquid-ordered structures, the latter known as lipid rafts. Here, we study spin-labeled ibuprofen (ibuprofen-SL) in the model membrane consisting of 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC), 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC), and cholesterol in the molar ratio of (0.5-0.5xchol)/(0.5-0.5xchol)/xchol. Electron paramagnetic resonance (EPR) spectroscopy is employed, along with its pulsed version of double electron-electron resonance (DEER, also known as PELDOR). The data obtained indicate lateral lipid-mediated clustering of ibuprofen-SL molecules with a local surface density noticeably larger than that expected for random lateral distribution. In the absence of cholesterol, the data can be interpreted as indicating alternating clustering in two opposing leaflets of the bilayer. In the presence of cholesterol, for xchol ≥ 20 mol %, the results show that ibuprofen-SL molecules have a quasi-regular lateral distribution, with a "superlattice" parameter of ∼3.0 nm. This regularity can be explained by the entrapment of ibuprofen-SL molecules by lipid rafts known to exist in this system with the additional assumption that lipid rafts have a nanoscale substructure.

The Chemical Reactivity of Membrane Lipids

It is well-known that aqueous dispersions of phospholipids spontaneously assemble into bilayer structures. These structures have numerous applications across chemistry and materials science and form the fundamental structural unit of the biological membrane. The particular environment of the lipid bilayer, with a water-poor low dielectric core surrounded by a more polar and better hydrated interfacial region, gives the membrane particular biophysical and physicochemical properties and presents a unique environment for chemical reactions to occur. Many different types of molecule spanning a range of sizes, from dissolved gases through small organics to proteins, are able to interact with membranes and promote chemical changes to lipids that subsequently affect the physicochemical properties of the bilayer. This Review describes the chemical reactivity exhibited by lipids in their membrane form, with an emphasis on conditions where the lipids are well hydrated in the form of bilayers. Key topics include the following: lytic reactions of glyceryl esters, including hydrolysis, aminolysis, and transesterification; oxidation reactions of alkenes in unsaturated fatty acids and sterols, including autoxidation and oxidation by singlet oxygen; reactivity of headgroups, particularly with reactive carbonyl species; and E/Z isomerization of alkenes. The consequences of reactivity for biological activity and biophysical properties are also discussed.

Impact of the Unstirred Water Layer on the Permeation of Small-Molecule Drugs

Over the last two decades, numerous molecular dynamics (MD) simulation-based investigations have attempted to predict the membrane permeability to small-molecule drugs as indicators of their bioavailability, a majority of which utilize the inhomogeneous solubility diffusion (ISD) model. However, MD-based membrane permeability is routinely 3-4 orders of magnitude larger than the values measured with the intestinal perfusion technique. There have been contentious discussions on the sources of the large discrepancies, and the two indisputable, potentially dominant ones are the fixed protonation state of the permeant and the neglect of the unstirred water layer (UWL). Employing six small-molecule drugs of different biopharmaceutical classification system classes, the current MD study relies on the ISD model but introduces the (de)protonation of the permeant by characterizing the permeation free energy of both neutral and charged states. In addition, the role of the UWL as a potential resistance against permeation is explored. The new MD protocol closely mimics the nature of small-molecule permeation and yields estimates that agree well with in vivo intestinal permeability.

How is Membrane Permeation of Small Ionizable Molecules Affected by Protonation Kinetics?

According to the pH-partition hypothesis, the aqueous solution adjacent to a membrane is a mixture of the ionization states of the permeating molecule at fixed Henderson-Hasselbalch concentrations, such that each state passes through the membrane in parallel with its own specific permeability. An alternative view, based on the assumption that the rate of switching ionization states is instantaneous, represents the permeation of ionizable molecules via an effective Boltzmann-weighted average potential (BWAP). Such an assumption is used in constant-pH molecular dynamics simulations. The inhomogeneous solubility-diffusion framework can be used to compute the pH-dependent membrane permeability for each of these two limiting treatments. With biased WTM-eABF molecular dynamics simulations, we computed the potential of mean force and diffusivity of each ionization state of two weakly basic small molecules: nicotine, an addictive drug, and varenicline, a therapeutic for treating nicotine addiction. At pH = 7, the BWAP effective permeability is greater than that determined by pH-partitioning by a factor of 2.5 for nicotine and 5 for varenicline. To assess the importance of ionization kinetics, we present a Smoluchowski master equation that includes explicitly the protonation and deprotonation processes coupled with the diffusive motion across the membrane. At pH = 7, the increase in permeability due to the explicit ionization kinetics is negligible for both nicotine and varenicline. This finding is reaffirmed by combined Brownian dynamics and Markov state model simulations for estimating the permeability of nicotine while allowing changes in its ionization state. We conclude that for these molecules the pH-partition hypothesis correctly captures the physics of the permeation process. The small free energy barriers for the permeation of nicotine and varenicline in their deprotonated neutral forms play a crucial role in establishing the validity of the pH-partitioning mechanism. Essentially, BWAP fails because ionization kinetics are too slow on the time scale of membrane crossing to affect the permeation of small ionizable molecules such as nicotine and varenicline. For the singly protonated state of nicotine, the computational results agree well with experimental measurements (P1 = 1.29 × 10-7 cm/s), but the agreement for neutral (P0 = 6.12 cm/s) and doubly protonated nicotine (P2 = 3.70 × 10-13 cm/s) is slightly worse, likely due to factors associated with the aqueous boundary layer (neutral form) or leaks through paracellular pathways (doubly protonated form).

The asymmetric plasma membrane-A composite material combining different functionalities?: Balancing Barrier Function and Fluidity for Effective Signaling

One persistent puzzle in the life sciences is the asymmetric lipid composition of the cellular plasma membrane: while the exoplasmic leaflet is enriched in lipids carrying predominantly saturated fatty acids, the cytoplasmic leaflet hosts preferentially lipids with (poly-)unsaturated fatty acids. Given the high energy requirements necessary for cells to maintain this asymmetry, the question naturally arises regarding its inherent benefits. In this paper, we propose asymmetry to represent a potential solution for harmonizing two conflicting requirements for the plasma membrane: first, the need to build a barrier for the uncontrolled influx or efflux of substances; and second, the need to form a fluid and dynamic two-dimensional substrate for signaling processes. We hence view here the plasma membrane as a composite material, where the exoplasmic leaflet is mainly responsible for the functional integrity of the barrier and the cytoplasmic leaflet for fluidity. We reinforce the validity of the proposed mechanism by presenting quantitative data from the literature, along with multiple examples that bolster our model.

Capric Acid and Myristic Acid Permeability Enhancers in Curved Liposome Membranes

Liposomes are considered as advanced drug delivery systems for cancer treatment. A generation of pH-sensitive liposomes is being developed that use fatty acids (FAs) as a trigger for drug release in tumor tissues. However, FAs are also known to enhance permeability, and it is unclear whether FAs in liposomes may cause drug leakage or premature drug release. The passive permeability of the drug through the membrane of the liposome is thus a crucial factor for timely drug delivery. To investigate how the curvature and lipid composition of liposomes affect their passive permeability, coarse-grained molecular dynamics were performed. The permeability was determined with a counting method. Flat bilayers and three liposomes with varying diameters were studied, which had varying lipid compositions of dipalmitoylphosphatidylcholine, cholesterol, and deprotonated or neutral saturated FAs. The investigated permeants were water and two other small permeants, which have different free energy profiles (solubility) across the membrane. First, for the curvature effect, our results showed that curvature increases the water permeability by reducing the membrane thickness. The permeability increase for water is about a factor of 1.7 for the most curved membranes. However, a high curvature decreases permeability for permeants with free energy profiles that are a mix of wells and barriers in the headgroup region of the membrane. Importantly, the type of experimental setup is expected to play a dominant role in the permeability value, i.e., whether permeants are escaping or entering the liposomes. Second, for the composition effect, FAs decrease both the area per lipid (APL) and the membrane thickness, resulting in permeability increases of up to 55%. Cholesterol has a similar effect on the APL but has the opposite impact on membrane thickness and permeability. Therefore, FAs and cholesterol have opposing effects on permeability, with cholesterol's effect being slightly stronger in our simulated bilayers. As all permeability values were well within a factor of 2, and with liposomes usually being larger and less curved in experimental applications, it can be concluded that the passive drug release from a pH-sensitive liposome does not seem to be significantly affected by the presence of FAs.

Concentration-dependent mechanism of the binding behavior of ibuprofen to the cell membrane: A molecular dynamic simulation study

Ibuprofen is a commonly used drug for treating headaches, pain, and fever. The lipid bilayer is the primary and most important interface for drugs to interact with biological systems. However, the molecular interactions between ibuprofen and the cell membrane are not well understood. Our findings suggest that the interactions between ibuprofen and the bilayer involve multiple steps and depend on the concentration of the drug. At low concentrations of ibuprofen, it can bind to the surface of the lipid bilayer. The electrostatic and vdW energies of IBU-lipid at 0 ns of the simulation were -22.5 ± 3.2 and -5.9 ± 1.2 kj.mol-1 Fig. 2. In the following, the vdW energy of the IBU-lipid was increased by around -134.6 ± 3.7 kj.mol-1 whereas the electrostatic energy of the IBU-lipid was significantly decreased. This binding is facilitated by electrostatic and vdW interactions between ibuprofen and the head group of lipids. In the second step, ibuprofen is inserted into the lipid bilayer and positioned at the interface between the bilayer and the aqueous phase. In high concentrations of ibuprofen, it moved to the central region of the lipid bilayer. At this concentration, the physical and structural properties of the cell membrane change significantly. Results from the radial distribution function analysis indicate that at low concentrations, ibuprofen molecules are situated close to the head groups of phosphate groups. However, at high concentrations of ibuprofen, these molecules move to the inner side of the lipid bilayer. In addition, our findings indicate that at low concentrations of ibuprofen, these molecules did not significantly alter the physical properties of the cell membrane. In contrast, at high concentrations of ibuprofen, the physical parameters of the hydrocarbon tails, such as thickness, fluidity, and order, changed dramatically. APL parameter for POPC membrane increased slightly to 0.60 and 0.63 nm2 in the presence of low and high concentrations of ibuprofen molecules. The three-step interaction between ibuprofen and the lipid bilayer involves several events, such as the movement of ibuprofen molecules towards the central region of the lipid bilayer and the deformation and alteration of the structural and stability properties of the cell membrane. These effects are observed only at high concentrations of ibuprofen. It appears that the side effects of ibuprofen overdose are related to changes in the properties of the cell membrane and, subsequently, the function of membrane-anchored target proteins.

Cryo-EM elucidates the uroplakin complex structure within liquid-crystalline lipids in the porcine urothelial membrane

The urothelium, a distinct epithelial tissue lining the urinary tract, serves as an essential component in preserving urinary tract integrity and thwarting infections. The asymmetric unit membrane (AUM), primarily composed of the uroplakin complex, constitutes a critical permeability barrier in fulfilling this role. However, the molecular architectures of both the AUM and the uroplakin complex have remained enigmatic due to the paucity of high-resolution structural data. In this study, we utilized cryo-electron microscopy to elucidate the three-dimensional structure of the uroplakin complex within the porcine AUM. While the global resolution achieved was 3.5 Å, we acknowledge that due to orientation bias, the resolution in the vertical direction was determined to be 6.3 Å. Our findings unveiled that the uroplakin complexes are situated within hexagonally arranged crystalline lipid membrane domains, rich in hexosylceramides. Moreover, our research rectifies a misconception in a previous model by confirming the existence of a domain initially believed to be absent, and pinpointing the accurate location of a crucial Escherichia coli binding site implicated in urinary tract infections. These discoveries offer valuable insights into the molecular underpinnings governing the permeability barrier function of the urothelium and the orchestrated lipid phase formation within the plasma membrane.

Permeability of antioxidants through a lipid bilayer model with coarse-grained simulations

Oxidative stress caused by pollution and lifestyle changes causes an excess of free radicals that react chemically with cell constituents leading to irreversible damage. There are molecules known as antioxidants that reduce the levels of free radicals. Some pigments of fruits and vegetables known as anthocyanins have antioxidant properties. Their interaction with the cell membrane becomes a crucial step in studying these substances. In this research, molecular dynamics simulations, particularly, coarse-grained molecular dynamics (CGMD) were used. This technique aims to replace functional groups with corresponding beads that represent their level of polarity and affinities to other chemical groups. Also, umbrella sampling was carried out to obtain free energy profiles that describe well the orientation and location of antioxidants in a membrane considering Trolox, Cyanidin, Gallic Acid, and Resveratrol molecules to study the structural effects they cause on it. It was concluded in this study that an antioxidant when crossing the membrane does not cause either damage to the structural properties or the loss of packing and stratification of phospholipids. it was also observed that the most reactive part of the molecules could easily approach area A prone to lipid oxidation, which can describe the antioxidant capacity of these molecules.Communicated by Ramaswamy H. Sarma.

Combining SILCS and Artificial Intelligence for High-Throughput Prediction of the Passive Permeability of Drug Molecules

Membrane permeability of drug molecules plays a significant role in the development of new therapeutic agents. Accordingly, methods to predict the passive permeability of drug candidates during a medicinal chemistry campaign offer the potential to accelerate the drug design process. In this work, we combine the physics-based site identification by ligand competitive saturation (SILCS) method and data-driven artificial intelligence (AI) to create a high-throughput predictive model for the passive permeability of druglike molecules. In this study, we present a comparative analysis of four regression models to predict membrane permeabilities of small druglike molecules; of the tested models, Random Forest was the most predictive yielding an R2 of 0.81 for the independent data set. The input feature vector used to train the developed prediction model includes absolute free energy profiles of ligands through a POPC-cholesterol bilayer based on ligand grid free energy (LGFE) profiles obtained from the SILCS approach. The use of the membrane free energy profiles from SILCS offers information on the physical forces contributing to ligand permeability, while the use of AI yields a more predictive model trained on experimental PAMPA permeability data for a collection of 229 molecules. This combination allows for rapid estimations of ligand permeability at a level of accuracy beyond currently available predictive models while offering insights into the contributions of the functional groups in the ligands to the permeability barrier, thereby offering quantitative information to facilitate rational ligand design.

Peptide-lipid hybrid vesicles with stimuli-responsive phase separation for controlled membrane functions

A disulfide-tethered peptide-lipid conjugate self-assembled into a homogeneously distributed peptide-lipid hybrid vesicle. Upon dithiothreitol treatment, the homogeneous peptide-lipid membrane spontaneously divided into lipid-rich and peptide-rich domains, while the vesicle retained its size and shape. Membrane phase separation enhanced temperature-dependent cargo release.

Predictions from First-Principles of Membrane Permeability to Small Molecules: How Useful Are They in Practice?

Predicting from first-principles the rate of passive permeation of small molecules across the biological membrane represents a promising strategy for screening lead compounds upstream in the drug-discovery and development pipeline. One popular avenue for the estimation of permeation rates rests on computer simulations in conjunction with the inhomogeneous solubility-diffusion model, which requires the determination of the free-energy change and position-dependent diffusivity of the substrate along the translocation pathway through the lipid bilayer. In this Perspective, we will clarify the physical meaning of the membrane permeability inferred from such computer simulations, and how theoretical predictions actually relate to what is commonly measured experimentally. We will also examine why these calculations remain both technically challenging and overly computationally expensive, which has hitherto precluded their routine use in nonacademic settings. We finally synopsize possible research directions to meet these challenges, increase the predictive power of physics-based rates of passive permeation, and, by ricochet, improve their practical usefulness.

Unsaturated Lipids Facilitate Partitioning of Transmembrane Peptides into the Liquid Ordered Phase

The affinity of single-pass transmembrane (TM) proteins for ordered membrane phases has been reported to depend on their lipidation, TM length, and lipid accessible surface area. In this work, the raft affinities of the TM domain of the linker for activation of T cells and its depalmitoylated variant are assessed using free energy simulations in a binary bilayer system composed of two laterally patched bilayers of ternary liquid ordered (Lo) and liquid disordered (Ld) phases. These phases are modeled by distinct compositions of distearoylphosphatidylcholine, palmitoyloleoylphosphatidylcholine (POPC), and cholesterol, and the simulations were carried out for 4.5 μs/window. Both peptides are shown to preferentially partition into the Ld phase in agreement with model membrane experiments and previous simulations on ternary lipid mixtures but not with measurements on giant plasma membrane vesicles where the Lo is slightly preferred. However, the 500 ns average relaxation time of lipid rearrangement around the peptide precluded a quantitative analysis of free energy differences arising from peptide palmitoylation and two different lipid compositions. When in the Lo phase, peptides reside in regions rich in POPC and interact preferentially with its unsaturated tail. Hence, the detailed substructure of the Lo phase is an important modulator of peptide partitioning, in addition to the inherent properties of the peptide.

Unveiling Liquid-Crystalline Lipids in the Urothelial Membrane through Cryo-EM

The urothelium, a distinct epithelial tissue lining the urinary tract, serves as an essential component in preserving urinary tract integrity and thwarting infections. The asymmetric unit membrane (AUM), primarily composed of the uroplakin complex, constitutes a critical permeability barrier in fulfilling this role. However, the molecular architectures of both the AUM and the uroplakin complex have remained enigmatic due to the paucity of high-resolution structural data. In this study, we utilized cryo-electron microscopy to elucidate the three-dimensional structure of the uroplakin complex within the porcine AUM. While the global resolution achieved was 3.5 Å, we acknowledge that due to orientation bias, the resolution in the vertical direction was determined to be 6.3 Å. Our findings unveiled that the uroplakin complexes are situated within hexagonally arranged crystalline lipid membrane domains, rich in hexosylceramides. Moreover, our research rectifies a misconception in a previous model by confirming the existence of a domain initially believed to be absent, and pinpointing the accurate location of a crucial Escherichia coli binding site implicated in urinary tract infections. These discoveries offer valuable insights into the molecular underpinnings governing the permeability barrier function of the urothelium and the orchestrated lipid phase formation within the plasma membrane.

Unveiling Liquid-Crystalline Lipids in the Urothelial Membrane through Cryo-EM

The urothelium, a distinct epithelial tissue lining the urinary tract, serves as an essential component in preserving urinary tract integrity and thwarting infections. The asymmetric unit membrane (AUM), primarily composed of the uroplakin complex, constitutes a critical permeability barrier in fulfilling this role. However, the molecular architectures of both the AUM and the uroplakin complex have remained enigmatic due to the paucity of high-resolution structural data. In this investigation, we employed cryo-electron microscopy to elucidate the three-dimensional structure of the uroplakin complex embedded within the porcine AUM at a resolution of 3.5 Å. Our findings unveiled that the uroplakin complexes are situated within hexagonally arranged crystalline lipid membrane domains, rich in hexosylceramides. Moreover, our research rectifies a misconception in a previous model by confirming the existence of a domain initially believed to be absent, and pinpointing the accurate location of a crucial Escherichia coli binding site implicated in urinary tract infections. These discoveries offer valuable insights into the molecular underpinnings governing the permeability barrier function of the urothelium and the orchestrated lipid phase formation within the plasma membrane.

CHARMM-GUI Membrane Builder: Past, Current, and Future Developments and Applications

Molecular dynamics simulations of membranes and membrane proteins serve as computational microscopes, revealing coordinated events at the membrane interface. As G protein-coupled receptors, ion channels, transporters, and membrane-bound enzymes are important drug targets, understanding their drug binding and action mechanisms in a realistic membrane becomes critical. Advances in materials science and physical chemistry further demand an atomistic understanding of lipid domains and interactions between materials and membranes. Despite a wide range of membrane simulation studies, generating a complex membrane assembly remains challenging. Here, we review the capability of CHARMM-GUI Membrane Builder in the context of emerging research demands, as well as the application examples from the CHARMM-GUI user community, including membrane biophysics, membrane protein drug-binding and dynamics, protein-lipid interactions, and nano-bio interface. We also provide our perspective on future Membrane Builder development.

Molecular simulations and NMR reveal how lipid fluctuations affect membrane mechanics

Lipid bilayers form the main matrix of functional cell membranes, and their dynamics underlie a host of physical and biological processes. Here we show that elastic membrane properties and collective molecular dynamics (MD) are related by the mean-square amplitudes (order parameters) and relaxation rates (correlation times) of lipid acyl chain motions. We performed all-atom MD simulations of liquid-crystalline bilayers that allow direct comparison with carbon-hydrogen (CH) bond relaxations measured with NMR spectroscopy. Previous computational and theoretical approaches have assumed isotropic relaxation, which yields inaccurate description of lipid chain dynamics and incorrect data interpretation. Instead, the new framework includes a fixed bilayer normal (director axis) and restricted anisotropic motion of the CH bonds in accord with their segmental order parameters, enabling robust validation of lipid force fields. Simulated spectral densities of thermally excited CH bond fluctuations exhibited well-defined spin-lattice (Zeeman) relaxations analogous to those in NMR measurements. Their frequency signature could be fit to a simple power-law function, indicative of nematic-like collective dynamics. Moreover, calculated relaxation rates scaled as the squared order parameters yielding an apparent κC modulus for bilayer bending. Our results show a strong correlation with κC values obtained from solid-state NMR studies of bilayers without and with cholesterol as validated by neutron spin-echo measurements of membrane elasticity. The simulations uncover a critical role of interleaflet coupling in membrane mechanics and thus provide important insights into molecular sites of emerging elastic properties within lipid bilayers.

Enzymatic Janus Liposome Micromotors

A liposome-based micromotor system that utilizes regional enzymatic conversion and gas generation to achieve directional motion in water is presented. Constituted mainly of a low-melting lipid and a high-melting lipid together with cholesterol, these liposomes maintain stable Janus configuration at room temperature as a result of lipid liquid-liquid phase separation. Local placement of enzymes such as horseradish peroxidase is realized via affinity binding between avidin and biotin, the latter as a lipid conjugate sorted specifically into one domain of these Janus liposomes as a minor component. In the presence of the substrate, hydrogen peroxide, these enzyme-decorated Janus liposomes undergo directional motion, yielding velocities exceeding thermal diffusion by three folds in some cases. Experimental details on liposome size control, motor assembly, and substrate distribution are presented; effects of key experimental factors on liposome motion, such as substrate concentration and liposome Janus ratio, are also examined. This work thus provides a viable approach to building asymmetrical lipid-assembled, enzyme-attached colloids and, in addition, stresses the importance of asymmetry in achieving particle directional motion.

Path sampling with memory reduction and replica exchange to reach long permeation timescales

Assessing kinetics in biological processes with molecular dynamics simulations remains a computational and conceptual challenge, given the large time and length scales involved. For kinetic transport of biochemical compounds or drug molecules, the permeability through the phospholipid membranes is a key kinetic property, but long timescales are hindering the accurate computation. Technological advances in high-performance computing therefore need to be accompanied by theoretical and methodological developments. In this contribution, the replica exchange transition interface sampling (RETIS) methodology is shown to give perspective toward observing longer permeation pathways. It is first reviewed how RETIS, a path-sampling methodology that gives in principle exact kinetics, can be used to compute membrane permeability. Next, recent and current developments in three RETIS aspects are discussed: several new Monte Carlo moves in the path-sampling algorithm, memory reduction by reducing pathlengths, and exploitation of parallel computing with CPU-imbalanced replicas. Finally, the memory reduction presenting a new replica exchange implementation, coined REPPTIS, is showcased with a permeant needing to pass a membrane with two permeation channels, either representing an entropic or energetic barrier. The REPPTIS results showed clearly that inclusion of some memory and enhancing ergodic sampling via replica exchange moves are both necessary to obtain correct permeability estimates. In an additional example, ibuprofen permeation through a dipalmitoylphosphatidylcholine membrane was modeled. REPPTIS succeeded in estimating the permeability of this amphiphilic drug molecule with metastable states along the permeation pathway. In conclusion, the presented methodological advances allow for deeper insight into membrane biophysics even if the pathways are slow, as RETIS and REPPTIS push the permeability calculations to longer timescales.

Structural Analyses of the Glycolipids in Lipid Rafts

Lipid rafts are usually isolated from cells or tissues using sucrose gradient ultracentrifugation in the presence of detergents such as Triton X-100 at 4 °C. Although detergents should be removed for further structural characterization following fractionation, these compounds are often difficult to completely remove, especially from the glycolipids. In this chapter, we describe a novel method for the fast and convenient removal of detergents from lipid raft glycolipids following fraction and describe the application of this method.

Opportunities for Ivory Nut Residue Valorization as a Source of Nanocellulose Colloidal Suspensions

Ivory nut seeds have been traditionally exploited in Central and South America for obtaining vegetable ivory. The residues from this industry are susceptible to valorization as a source of fatty acids (by organic extraction) and mannans (by alkaline dissolution and regeneration). Nonetheless, cellulose may also be recovered at the end of this fractionation process by acid hydrolysis and functionalization, with associated advantages over other lignocellulosic sources due to the absence of lignin in the endospermic tissue. In this work, various experimental parameters (sulfuric acid concentration, temperature, and hydrolysis time) were investigated to optimize the processing conditions for preparing stable nanocellulose suspensions after ultrasonication. The most stable nanocellulose gel (1 wt% solid content) was obtained after 4-h hydrolysis at 60 °C with 8 M H2SO4 and was characterized by using complementary tech-niques, including dynamic light scattering (DLS), transmission electron microscopy (TEM), X-ray powder diffraction (XRD), nano-fibril sulfation measurements, vibrational and solid-state nuclear magnetic resonance (CP/MAS 13C-NMR) spectroscopies, and thermal analysis. This nanocellulose hydrogel is susceptible to further utilization in various applications and fields, e.g., in agricul-ture for controlling the release of agrochemicals, in pharmaceutics for developing new dosage forms, and in the treatment of wastewater from the textile and paper industries.

In Vitro Membrane Platform for the Visualization of Water Impermeability across the Liquid-Ordered Phase under Hypertonic Conditions

Passive water penetration across the cell membrane by osmotic diffusion is essential for the homeostasis of cell volume, in addition to the protein-assisted active transportation of water. Since membrane components can regulate water permeability, controlling compositional variation during the volume regulatory process is a prerequisite for investigating the underlying mechanisms of water permeation and related membrane dynamics. However, the lack of a viable in vitro membrane platform in hypertonic solutions impedes advanced knowledge of cell volume regulation processes, especially cholesterol-enriched lipid domains called lipid rafts. By reconstituting the liquid-ordered (L_o) domain as a likeness of lipid rafts, we verified suppressed water permeation across the L_o domains, which had yet to be confirmed with experimental demonstrations despite a simulation approach. With the help of direct transfer of the L_o domains from vesicles to supported lipid membranes, the biological roles of lipid composition in suppressed water translocation were experimentally confirmed. Additionally, the improvement in membrane stability under hypertonic conditions was demonstrated based on molecular dynamics simulations.

Planar aggregation of the influenza viral fusion peptide alters membrane structure and hydration, promoting poration

To infect, enveloped viruses employ spike protein, spearheaded by its amphipathic fusion peptide (FP), that upon activation extends out from the viral surface to embed into the target cellular membrane. Here we report that synthesized influenza virus FPs are membrane active, generating pores in giant unilamellar vesicles (GUV), and thus potentially explain both influenza virus' hemolytic activity and the liposome poration seen in cryo-electron tomography. Experimentally, FPs are heterogeneously distributed on the GUV at the time of poration. Consistent with this heterogeneous distribution, molecular dynamics (MD) simulations of asymmetric bilayers with different numbers of FPs in one leaflet show FP aggregation. At the center of FP aggregates, a profound change in the membrane structure results in thinning, higher water permeability, and curvature. Ultimately, a hybrid bilayer nanodomain forms with one lipidic leaflet and one peptidic leaflet. Membrane elastic theory predicts a reduced barrier to water pore formation when even a dimer of FPs thins the membrane as above, and the FPs of that dimer tilt, to continue the leaflet bending initiated by the hydrophobic mismatch between the FP dimer and the surrounding lipid.

Lipid Rafts and Plant Gravisensitivity

The necessity to include plants as a component of a Bioregenerative Life Support System leads to investigations to optimize plant growth facilities as well as a better understanding of the plant cell membrane and its numerous activities in the signaling, transport, and sensing of gravity, drought, and other stressors. The cell membrane participates in numerous processes, including endo- and exocytosis and cell division, and is involved in the response to external stimuli. Variable but stabilized microdomains form in membranes that include specific lipids and proteins that became known as (detergent-resistant) membrane microdomains, or lipid rafts with various subclassifications. The composition, especially the sterol-dependent recruitment of specific proteins affects endo- and exo-membrane domains as well as plasmodesmata. The enhanced saturated fatty acid content in lipid rafts after clinorotation suggests increased rigidity and reduced membrane permeability as a primary response to abiotic and mechanical stress. These results can also be obtained with lipid-sensitive stains. The linkage of the CM to the cytoskeleton via rafts is part of the complex interactions between lipid microdomains, mechanosensitive ion channels, and the organization of the cytoskeleton. These intricately linked structures and functions provide multiple future research directions to elucidate the role of lipid rafts in physiological processes.

The Role of Cryoprotective Agents in Liposome Stabilization and Preservation

To improve liposomes’ usage as drug delivery vehicles, cryoprotectants can be utilized to prevent constituent leakage and liposome instability. Cryoprotective agents (CPAs) or cryoprotectants can protect liposomes from the mechanical stress of ice by vitrifying at a specific temperature, which forms a glassy matrix. The majority of studies on cryoprotectants demonstrate that as the concentration of the cryoprotectant is increased, the liposomal stability improves, resulting in decreased aggregation. The effectiveness of CPAs in maintaining liposome stability in the aqueous state essentially depends on a complex interaction between protectants and bilayer composition. Furthermore, different types of CPAs have distinct effective mechanisms of action; therefore, the combination of several cryoprotectants may be beneficial and novel attributed to the synergistic actions of the CPAs. In this review, we discuss the use of liposomes as drug delivery vehicles, phospholipid–CPA interactions, their thermotropic behavior during freezing, types of CPA and their mechanism for preventing leakage of drugs from liposomes.

Deciphering the origin of the melting profile of unilamellar phosphatidylcholine liposomes by measuring the turbidity of its suspensions

The elucidation of the thermal properties of phosphatidylcholine liposomes is often based on the analysis of the thermal capacity profiles of multilamellar liposomes (MLV), which may qualitatively disagree with those of unilamellar liposomes (LUV). Experiments and interpretation of LUV liposomes is further complicated by aggregation and lamellarization of lipid bilayers in a short time period, which makes it almost impossible to distinguish the signatures of the two types of bilayers. To characterize independently MLV and LUV of 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC), the latter were prepared with the addition of small amounts of 1,2-dipalmitoyl-sn-glycero-3-phosphatidylglycerol (DPPG) which, due to the sterical hindrance and negative charge at a given pH value, cause LUV repellence and contribute to their stability. Differential scanning calorimetry curves and temperature-dependent UV/Vis spectra of the prepared MLV and LUV were measured. Multivariate analysis of spectrophotometric data determined the phase transition temperatures (pretransition at T_p and the main phase transition at T_m), and based on the changes in turbidities, the thickness of the lipid bilayer in LUV was determined. The obtained data suggested that the curvature change is a key distinguishing factor in MLV and LUV heat capacity profiles. By combining the experimental results and those obtained by MD simulations, the interfacial water layer was characterized and its contribution to the thermal properties of LUV was discussed.

Structural and mechanical properties of the red blood cell’s cytoplasmic membrane seen through the lens of biophysics

Red blood cells (RBCs) are the most abundant cell type in the human body and critical suppliers of oxygen. The cells are characterized by a simple structure with no internal organelles. Their two-layered outer shell is composed of a cytoplasmic membrane (RBC_cm) tethered to a spectrin cytoskeleton allowing the cell to be both flexible yet resistant against shear stress. These mechanical properties are intrinsically linked to the molecular composition and organization of their shell. The cytoplasmic membrane is expected to dominate the elastic behavior on small, nanometer length scales, which are most relevant for cellular processes that take place between the fibrils of the cytoskeleton. Several pathologies have been linked to structural and compositional changes within the RBC_cm and the cell’s mechanical properties. We review current findings in terms of RBC lipidomics, lipid organization and elastic properties with a focus on biophysical techniques, such as X-ray and neutron scattering, and Molecular Dynamics simulations, and their biological relevance. In our current understanding, the RBC_cm’s structure is patchy, with nanometer sized liquid ordered and disordered lipid, and peptide domains. At the same time, it is surprisingly soft, with bending rigidities κ of 2–4 k_BT. This is in strong contrast to the current belief that a high concentration of cholesterol results in stiff membranes. This extreme softness is likely the result of an interaction between polyunsaturated lipids and cholesterol, which may also occur in other biological membranes. There is strong evidence in the literature that there is no length scale dependence of κ of whole RBCs.

The bending rigidity of the red blood cell cytoplasmic membrane

An important mechanical property of cells is the membrane bending modulus, κ. In the case of red blood cells (RBCs) there is a composite membrane consisting of a cytoplasmic membrane and an underlying spectrin network. Literature values of κ are puzzling, as they are reported over a wide range, from 5 kBT to 230 kBT. To disentangle the contribution of the cytoplasmic membrane from the spectrin network, we investigated the bending of red blood cell cytoplasmic membranes (RBCcm) in the absence of spectrin and adenosine triphosphate (ATP). We used a combination of X-ray diffuse scattering (XDS), neutron spin-echo (NSE) spectrometry and Molecular Dynamics (MD) simulations. Our results indicate values of κ of order 4 kBT to 6 kBT, relatively small compared to literature values for most single component lipid bilayers. We suggest two ways this relative softness might confer biological advantage.

Dynamics of the Acinetobacter baumannii inner membrane under exogenous polyunsaturated fatty acid stress

Exogenous polyunsaturated fatty acids (PUFAs) are readily incorporated into the synthesis pathways of A. baumannii membrane phospholipids, where they contribute to reduced bacterial fitness and increased antimicrobial susceptibility. Here we examine the impact of PUFA membrane modification on membrane organisation and biophysical properties using coarse grained MARTINI simulations of chemically representative membrane models developed from mass-spectrometry datasets of an untreated, arachidonic acid (AA) treated and docosahexaenoic acid (DHA) treated A. baumannii membranes. Enzymatic integration of AA or DHA into phospholipids of the A. baumannii membrane resulted in modulation of membrane biophysical properties. Membrane thickness decreased slightly following PUFA treatment, concomitant with changes in the lateral area per lipid of each lipid headgroup class. PUFA treatment resulted in a decrease in membrane ordering and an increase in lipid lateral diffusion. Changes in lateral membrane organisation were observed in the PUFA treated membranes, with a concurrent increase in ordered cardiolipin domains and disordered PUFA-containing domains. Notably, separation between ordered and disordered domains was enhanced and was more pronounced for DHA relative to AA, providing a possible mechanism for greater antimicrobial action of DHA relative to AA observed experimentally. Furthermore, the membrane active antimicrobial, pentamidine, preferentially adsorbs to cardiolipin domains of the A. baumannii model membranes. This interaction, and membrane penetration of pentamidine, was enhanced following PUFA treatment. Cumulatively, this work explores the wide-ranging effects of PUFA incorporation on the A. baumannii membrane and provides a molecular basis for bacterial inner membrane disruption by PUFAs.

Phospholipid-porphyrin conjugates: deciphering the driving forces behind their supramolecular assemblies

Phospholipid-porphyrin conjugates (PL-Por) are nowadays considered as a unique class of building blocks that can self-assemble into supramolecular structures that possess multifunctional properties and enhanced optoelectronics characteristics compared to their disassembled counterparts. However, despite their versatile properties, little is known about the impact of the packing parameter of PL-Por conjugates on their assembling mechanism and their molecular organization inside these assemblies. To gain a better understanding on their assembling properties, we synthesized two new series of PL-Por conjugates with different alkyl sn2-chain lengths linked via an amide bond to either pheophorbide-a (Ph_xLPC) or pyropheophorbide-a (Pyr_xLPC). By combining a variety of experimental techniques with molecular dynamics (MD) simulations, we investigated both the assembling and optical properties of the PL-Por either self-assembled or when incorporated into lipid bilayers. We demonstrated that independently of the linker length, Ph_xLPC assembled into closed ovoid structures, whereas Pyr_xLPC formed rigid open sheets. Interestingly, Pyr_xLPC assemblies displayed a significant red shift and narrowing of the Q-band indicating the formation of ordered J-aggregates. The MD simulations highlighted the central role of the interaction between porphyrin cores rather than the length difference between the two phospholipid chains in controlling the structure of the lipid bilayer membranes and thus their optical properties. Indeed, while Ph_xLPC have the tendency to form inter-leaflet π-stacked dimers, Pyr_xLPC conjugates formed dimers within the same leaflet. Altogether, this work could be used as guidelines for the design of new PL-Por conjugates that self-assemble into bilayer-like supramolecular structures with tunable morphology and optical properties.

Computational development of a phase-sensitive membrane raft probe

Derivatives of the widely used 1,6-diphenyl-1,3,5-hexatriene molecular probe have been considered using a multiscale approach involving spin-flip time-dependent density functional theory, classical molecular dynamics and hybrid quantum mechanics/molecular mechanics. We identify a potential probe of membrane phase (i.e. to preferentially detect liquid-ordered regions of lipid bilayers), which exhibits restricted access to a conical intersection in the liquid-ordered phase but is freely accessible in less ordered molecular environments. The characteristics of this probe also mark it as a candidate for an aggregation induced emission fluorophore.

Membrane thickness, lipid phase and sterol type are determining factors in the permeability of membranes to small solutes

Cell membranes provide a selective semi-permeable barrier to the passive transport of molecules. This property differs greatly between organisms. While the cytoplasmic membrane of bacterial cells is highly permeable for weak acids and glycerol, yeasts can maintain large concentration gradients. Here we show that such differences can arise from the physical state of the plasma membrane. By combining stopped-flow kinetic measurements with molecular dynamics simulations, we performed a systematic analysis of the permeability of a variety of small molecules through synthetic membranes of different lipid composition to obtain detailed molecular insight into the permeation mechanisms. While membrane thickness is an important parameter for the permeability through fluid membranes, the largest differences occur when the membranes transit from the liquid-disordered to liquid-ordered and/or to gel state, which is in agreement with previous work on passive diffusion of water. By comparing our results with in vivo measurements from yeast, we conclude that the yeast membrane exists in a highly ordered and rigid state, which is comparable to synthetic saturated DPPC-sterol membranes.

Membrane permeability of small molecules depends on the composition of the lipid bilayer. Here, authors compare permeability measured on membranes in different physical states and conclude that the yeast membrane exists in a highly ordered phase.

Biosafety-inspired structural optimization of triazolium ionic liquids based on structure-toxicity relationships

Ionic liquids (ILs), owing to their low vapor pressure and excellent solvating ability, are being increasingly applied in various industries to replace highly toxic organic solvents. They mainly pollute aquatic environment and soils, directly endangering eco-environment and human health. Therefore, it is critical to understand and optimize structural motifs of ILs with reduced toxicity. Considering human oral exposure is the major route, our investigations employed a human cell panel (modeling oral exposures) including human stomach (GES-1), intestinal (FHC), liver (HepG2) and kidney (HEK293) cells using a series of experimental and computational approaches to explore the cytotoxicity and molecular mechanism of ILs. We discovered that the cytotoxicity of triazolium and imidazolium ILs was human cell line-dependent with cytotoxicity in an order of FHC > GES-1 > HepG2 > HEK293. For this reason, a toxicity assay using a single cell line was highly inappropriate. Compared to anions (Br^-, OTs^-, OTMBS^-) we tested, the cation of ILs played a major role in causing cytotoxicity. Ionic liquids with cations having longer hydrophobic sidechains (IL09 vs. IL01) readily insert into cell membranes with enhanced membrane and lipidomic perturbations, induce cytotoxicity by triggering cell cycle arrest and apoptosis. Reducing sidechain length and incorporating three nitrogen atoms (triazolium) instead of two (imidazolium) in the cation core alleviated cytotoxicity by reducing cell membrane perturbations and cell function interference. These findings provide important guiding principles for the design of the next-generation of "green" and safe ILs.

Oxygen Storage in Stacked Phospholipid Membranes Under an Oxygen Gradient as a Model for Myelin Sheaths

Axons in the brain and peripheral nervous system are enveloped by myelin sheaths, which are composed of stacked membrane bilayers containing large fractions of cholesterol, phospholipids, and glycolipids. The oxygen availability to the nearby oxygen consuming cytochrome c oxidase in the mitochondria is essential for the well-functioning of a cell. By constructing a rate network model based on molecular dynamics simulations, and solving it for steady-state conditions, this work calculates the oxygen storage in stacked membranes under an oxygen gradient. It is found that stacking membranes increases the oxygen storage capacity, indicating that myelin can function as an oxygen reservoir. However, it is found that the storage enhancement levels out for stacks with a large number of bilayers, suggesting why myelin sheaths consist of only 10-300 membranes rather than thousands. The presence of additional water between the stacked bilayers, as seen in cancer cells, is shown to diminish myelin oxygen storage enhancement.

Designing amphiphilic Janus nanoparticles with tunable lipid raft affinity via molecular dynamics simulation

Due to the differential interactions among lipids and proteins, the plasma membrane can segregate into a series of functional nanoscale membrane domains ("lipid rafts"), which are essential in multiple biological processes such as signaling transduction, protein trafficking and endocytosis. On the other hand, Janus nanoparticles (NPs) have shown great promise in various biomedical applications due to their asymmetric characteristics and can integrate different surface properties and thus synergetic functions. Hence, in this work, we aim to design an amphiphilic Janus NP to target and regulate lipid rafts via tuning its surface ligand amphiphilicity using coarse-grained molecular dynamics (MD) simulations. Our μs-scale free coarse-grained MD simulations as well as umbrella sampling free energy calculations indicated that the hydrophobicity of the hydrophobic surface ligands not only determined the lateral membrane partitioning thermodynamics of Janus NPs in phase-separated lipid membranes, but also the difficulty in their insertion into different membrane domains of the lipid membrane. These two factors jointly regulated the lipid raft affinity of Janus NPs. Meanwhile, the hydrophilicity of the hydrophilic surface ligands could affect the insertion ability of Janus NPs. Besides, the ultra-small size could ensure the membrane-bound behavior of Janus NPs without disrupting the overall structure and phase separation kinetics of the lipid membrane. These results may provide valuable insights into the design of functional NPs targeting and controllably regulating lipid rafts.

Blood bank storage of red blood cells increases RBC cytoplasmic membrane order and bending rigidity

Blood banks around the world store blood components for several weeks ensuring its availability for transfusion medicine. Red blood cells (RBCs) are known to undergo compositional changes during storage, which may impact the cells' function and eventually the recipients' health. We extracted the RBC's cytoplasmic membrane (RBCcm) to study the effect of storage on the membranes' molecular structure and bending rigidity by a combination of X-ray diffraction (XRD), X-ray diffuse scattering (XDS) and coarse grained Molecular Dynamics (MD) simulations. Blood was stored in commercial blood bags for 2 and 5 weeks, respectively and compared to freshly drawn blood. Using mass spectrometry, we measured an increase of fatty acids together with a slight shift towards shorter tail lengths. We observe an increased fraction (6%) of liquid ordered (lo) domains in the RBCcms with storage time, and an increased lipid packing in these domains, leading to an increased membrane thickness and membrane order. The size of both, lo and liquid disordered (ld) lipid domains was found to decrease with increased storage time by up to 25%. XDS experiments reveal a storage dependent increase in the RBCcm's bending modulus κ by a factor of 2.8, from 1.9 kBT to 5.3 kBT. MD simulations were conducted in the absence of proteins. The results show that the membrane composition has a small contribution to the increased bending rigidity and suggests additional protein-driven mechanisms.

Site of Cholesterol Oxidation Impacts Its Localization and Domain Formation in the Neuronal Plasma Membrane

Cholesterol is integral to the structure of mammalian cell membranes. Oxidation of cholesterol alters how it behaves in the membrane and influences the membrane biophysical properties. Elevated levels of oxidized cholesterol are associated with neurodegenerative diseases such as Alzheimer's disease, Parkinson's disease, multiple sclerosis, and Huntington's disease. Previous work has investigated the impact of oxidized cholesterol in the context of simple model membrane systems. However, there is a growing body of literature that shows that complex membranes possessing physiological phospholipid distributions have different properties from those of binary or trinary model membranes. In the current work, the impact of oxidized cholesterol on the biophysical properties of a complex neuronal plasma membrane is investigated using coarse-grained Martini molecular dynamics simulations. Comparison of the native neuronal membrane to neuronal membranes containing 10% tail-oxidized or 10% head-oxidized cholesterol shows that the site of oxidization changes the behavior of the oxidized cholesterol in the membrane. Furthermore, species-specific domain formation is observed between each oxidized cholesterol and minor lipid classes. Although both tail-oxidized and head-oxidized cholesterols modulate the biophysical properties of the membrane, smaller changes are observed in the complex neuronal membrane than seen in the previous work on simple binary or trinary model membranes. This work highlights the presence of compensatory effects of lipid diversity in the complex neuronal membrane. Overall, this study improves our molecular-level understanding of the effects of oxidized cholesterol on the properties of neuronal tissue and emphasizes the importance of studying membranes with realistic lipid compositions.

Factors determining barrier properties to oxygen transport across model and cell plasma membranes based on EPR spin-label oximetry

This review is motivated by the exciting new area of radiation therapy using a phenomenon termed FLASH in which oxygen is thought to have a central role. Well-established principles of radiation biology and physics suggest that if oxygen has a strong role, it should be the level at the DNA. The key aspect discussed is the rate of oxygen diffusion. If oxygen freely diffuses into cells and rapidly equilibrates, then measurements in the extracellular compartment would enable FLASH to be investigated using existing methodologies that can readily measure oxygen in the extracellular compartment. EPR spin-label oximetry allows evaluation of the oxygen permeability coefficient across lipid bilayer membranes. It is established that simple fluid phase lipid bilayers are not barriers to oxygen transport. However, further investigations indicate that many physical and chemical (compositional) factor can significantly decrease this permeation. In biological cell plasma membranes, the lipid bilayer forms the matrix in which integral membrane proteins are immersed, changing organization and properties of the lipid matrix. To evaluate oxygen permeability coefficients across these complex membranes, oxygen permeation across all membrane domains and components must be considered. In this review, we consider many of the factors that affect (decrease) oxygen permeation across cell plasma membranes. Finally, we address the question, can the plasma membrane of the cell form a barrier to the free diffusion of oxygen into the cell interior? If there is a barrier then this must be considered in the investigations of the role of oxygen in FLASH.

Atomistic simulations modify interpretation of spin-label oximetry data. Part 1: intensified water-lipid interfacial resistances

The role of membrane cholesterol in cellular function and dysfunction has been the subject of much inquiry. A few studies have suggested that cholesterol may slow oxygen diffusive transport, altering membrane physical properties and reducing oxygen permeability. The primary experimental technique used in recent years to study membrane oxygen transport is saturation-recovery electron paramagnetic resonance (EPR) oximetry, using spin-label probes targeted to specific regions of a lipid bilayer. The technique has been used, in particular, to assess the influence of cholesterol on oxygen transport and membrane permeability. The reliability of such EPR recordings at the water-lipid interface near the phospholipid headgroups has been challenged by all-atom molecular dynamics (MD) simulation data that show substantive agreement with spin-label probe measurements throughout much of the bilayer. This work uses further MD simulations, with an updated oxygen model, to determine the location of the maximum resistance to permeation and the rate-limiting barrier to oxygen permeation in 1-palmitoyl,2-oleoylphosphatidylcholine (POPC) and POPC/cholesterol bilayers at 25 and 35°C. The current simulations show a spike of resistance to permeation in the headgroup region that was not detected by EPR but was predicted in early theoretical work by Diamond and Katz. Published experimental nuclear magnetic resonance (NMR) oxygen measurements provide key validation of the MD models and indicate that the positions and relative magnitudes of the phosphatidylcholine resistance peaks are accurate. Consideration of the headgroup-region resistances predicts bilayer permeability coefficients lower than estimated in EPR studies, giving permeabilities lower than the permeability of unstirred water layers of the same thickness. Here, the permeability of POPC at 35°C is estimated to be 13 cm/s, compared with 10 cm/s for POPC/cholesterol and 118 cm/s for simulation water layers of similar thickness. The value for POPC is 12 times lower than estimated from EPR measurements, while the value for POPC/cholesterol is 5 times lower. These findings underscore the value of atomic resolution models for guiding the interpretation of experimental probe-based measurements.

Data-driven equation for drug-membrane permeability across drugs and membranes

Drug efficacy depends on its capacity to permeate across the cell membrane. We consider the prediction of passive drug-membrane permeability coefficients. Beyond the widely recognized correlation with hydrophobicity, we additionally consider the functional relationship between passive permeation and acidity. To discover easily interpretable equations that explain the data well, we use the recently proposed sure-independence screening and sparsifying operator (SISSO), an artificial-intelligence technique that combines symbolic regression with compressed sensing. Our study is based on a large in silico dataset of 0.4 × 10⁶ small molecules extracted from coarse-grained simulations. We rationalize the equation suggested by SISSO via an analysis of the inhomogeneous solubility-diffusion model in several asymptotic acidity regimes. We further extend our analysis to the dependence on lipid-membrane composition. Lipid-tail unsaturation plays a key role but surprisingly contributes stepwise rather than proportionally. Our results are in line with previously observed changes in permeability, suggesting the distinction between liquid-disordered and liquid-ordered permeation. Together, compressed sensing with analytically derived asymptotes establish and validate an accurate, broadly applicable, and interpretable equation for passive permeability across both drug and lipid-tail chemistry.

Improving the accuracy and convergence of drug permeation simulations via machine-learned collective variables

Understanding the permeation of biomolecules through cellular membranes is critical for many biotechnological applications, including targeted drug delivery, pathogen detection, and the development of new antibiotics. To this end, computer simulations are routinely used to probe the underlying mechanisms of membrane permeation. Despite great progress and continued development, permeation simulations of realistic systems (e.g., more complex drug molecules or biologics through heterogeneous membranes) remain extremely challenging if not intractable. In this work, we combine molecular dynamics simulations with transition-tempered metadynamics and techniques from the variational approach to conformational dynamics to study the permeation mechanism of a drug molecule, trimethoprim, through a multicomponent membrane. We show that collective variables (CVs) obtained from an unsupervised machine learning algorithm called time-structure based Independent Component Analysis (tICA) improve performance and substantially accelerate convergence of permeation potential of mean force (PMF) calculations. The addition of cholesterol to the lipid bilayer is shown to increase both the width and height of the free energy barrier due to a condensing effect (lower area per lipid) and increase bilayer thickness. Additionally, the tICA CVs reveal a subtle effect of cholesterol increasing the resistance to permeation in the lipid head group region, which is not observed when canonical CVs are used. We conclude that the use of tICA CVs can enable more efficient PMF calculations with additional insight into the permeation mechanism.

The role of plasmalogens, Forssman lipids, and sphingolipid hydroxylation in modulating the biophysical properties of the epithelial plasma membrane

A coarse-grain model of the epithelial plasma membrane was developed from high-resolution lipidomic data and simulated using the MARTINI force field to characterize its biophysical properties. Plasmalogen lipids, Forssman glycosphingolipids, and hydroxylated Forssman glycosphingolipids and sphingomyelin were systematically added to determine their structural effects. Plasmalogen lipids have a minimal effect on the overall biophysical properties of the epithelial plasma membrane. In line with the hypothesized role of Forssman lipids in the epithelial apical membrane, the introduction of Forssman lipids initiates the formation of glycosphingolipid-rich nanoscale lipid domains, which also include phosphatidylethanolamine (PE), sphingomyelin (SM), and cholesterol (CHOL). This decreases the lateral diffusion in the extracellular leaflet, as well as the area per lipid of domain forming lipids, most notably PE. Finally, hydroxylation of the Forssman glycosphingolipids and sphingomyelin further modulates the lateral organization of the membrane. Through comparison to the previously studied average and neuronal plasma membranes, the impact of membrane lipid composition on membrane properties was characterized. Overall, this study furthers our understanding of the biophysical properties of complex membranes and the impact of lipid diversity in modulating membrane properties.

Sampling efficiency of the counting method for permeability calculations estimated with the inhomogeneous solubility-diffusion model

Permeability is a key property in various fields such as membrane technology for chemical separation and transport of substances through cell membranes. At the molecular scale, the counting method uses the number of membrane crossings in a conventional unbiased molecular dynamics simulation to predict the permeability. This contribution investigates under which conditions the counting method has insufficient statistics. An equation is derived for a compartmental model based on the inhomogeneous solubility-diffusion (Smoluchowski) model, giving insight into how the flux correlates with the solubility of permeants. This equation shows that a membrane crossing is a rare event not only when the membrane forms a large free energy barrier but also when the membrane forms a deep free energy well that traps permeants. Such a permeant trap has a high permeability; yet, the counting method suffers from poor statistics. To illustrate this, coarse-grained MD was run for 16 systems of dipalmitoylphosphatidylcholine bilayer membranes with different permeant types. The composition rule for permeability is shown to also hold for fluxes, and it is highlighted that the considered thickness of the membrane causes uncertainty in the permeability calculation of highly permeable membranes. In conclusion, a high permeability in itself is not an effective indicator of the sampling efficiency of the counting method, and caution should be taken for permeants whose solubility varies greatly over the simulation box. A practical consequence relevant in, e.g., drug design is that a drug with high membrane permeability might get trapped by membranes thus reducing its efficacy.

Lipid Oxidation: Role of Membrane Phase-Separated Domains

Lipid oxidation is associated with several inflammatory and neurodegenerative diseases, but many questions to unravel its effects on biomembranes are still open due to the complexity of the topic. For instance, recent studies indicated that phase-separated domains can have a significant effect on membrane function. It is reported that domain interfaces are "hot spots" for pore formation, but the underlying mechanisms and the effect of oxidation-induced phase separation on membranes remain elusive. Thus, to evaluate the permeability of the membrane coexisting of liquid-ordered (Lo) and liquid-disordered (Ld) domains, we performed atomistic molecular dynamics simulations. Specifically, we studied the membrane permeability of nonoxidized or oxidized homogeneous membranes (single-phase) and at the Lo/Ld domain interfaces of heterogeneous membranes, where the Ld domain is composed of either oxidized or nonoxidized lipids. Our simulation results reveal that the addition of only 1.5% of lipid aldehyde molecules at the Lo/Ld domain interfaces of heterogeneous membranes increases the membrane permeability, whereas their addition at homogeneous membranes does not have any effect. This study is of interest for a better understanding of cancer treatment methods based on oxidative stress (causing among others lipid oxidation), such as plasma medicine and photodynamic therapy.