Molecular View of Phase Coexistence in Lipid Monolayers

Review of the Interfacial Structure and Properties of Surfactants in Petroleum Production and Geological Storage Systems from a Molecular Scale Perspective

Surfactants play a crucial role in tertiary oil recovery by reducing the interfacial tension between immiscible phases, altering surface wettability, and improving foam film stability. Oil reservoirs have high temperatures and high pressures, making it difficult and hazardous to conduct lab experiments. In this context, molecular dynamics (MD) simulation is a valuable tool for complementing experiments. It can effectively study the microscopic behaviors (such as diffusion, adsorption, and aggregation) of the surfactant molecules in the pore fluids and predict the thermodynamics and kinetics of these systems with a high degree of accuracy. MD simulation also overcomes the limitations of traditional experiments, which often lack the necessary temporal-spatial resolution. Comparing simulated results with experimental data can provide a comprehensive explanation from a microscopic standpoint. This article reviews the state-of-the-art MD simulations of surfactant adsorption and resulting interfacial properties at gas/oil-water interfaces. Initially, the article discusses interfacial properties and methods for evaluating surfactant-formed monolayers, considering variations in interfacial concentration, molecular structure of the surfactants, and synergistic effect of surfactant mixtures. Then, it covers methods for characterizing microstructure at various interfaces and the evolution process of the monolayers' packing state as a function of interfacial concentration and the surfactants' molecular structure. Next, it examines the interactions between surfactants and the aqueous phase, focusing on headgroup solvation and counterion condensation. Finally, it analyzes the influence of hydrophobic phase molecular composition on interactions between surfactants and the hydrophobic phase. This review deepened our understanding of the micro-level mechanisms of oil displacement by surfactants and is beneficial for screening and designing surfactants for oil field applications.

A Review on Investigating the Interactions between Nanoparticles and the Pulmonary Surfactant Monolayer with Coarse-Grained Molecular Dynamics Method

Pulmonary drug delivery has garnered significant attention due to its targeted local lung action, minimal toxic side effects, and high drug utilization. However, the physicochemical properties of inhaled nanoparticles (NPs) used as drug carriers can influence their interactions with the pulmonary surfactant (PS) monolayer, potentially altering the fate of the NPs and impairing the biophysical function of the PS monolayer. Thus, the objective of this review is to summarize how the physicochemical properties of NPs affect their interactions with the PS monolayer. Initially, the definition and properties of NPs, as well as the composition and characteristics of the PS monolayer, are introduced. Subsequently, the coarse-grained molecular dynamics (CGMD) simulation method for studying the interactions between NPs and the PS monolayer is presented. Finally, the implications of the hydrophobicity, size, shape, surface charge, surface modification, and aggregation of NPs on their interactions with the PS monolayer and on the composition of biomolecular corona are discussed. In conclusion, gaining a deeper understanding of the effects of the physicochemical properties of NPs on their interactions with the PS monolayer will contribute to the development of safer and more effective nanomedicines for pulmonary drug delivery.

Numerical investigations of translocation characteristics of nano-silica lunar dust across pulmonary surfactant monolayer

The interactions between nano-silica lunar dust (NSLD) on the moon surface and pulmonary surfactant (PS) monolayer will pose risks to astronaut health in future manned lunar exploration missions, but the specifics of these interactions are unknown. This study investigates them using the coarse-grained molecular dynamics method considering different sizes (5, 10, and 15 nm) and shapes (sphere, ellipsoid, and cube), with special focus on the unique morphology of NSLDs with bugles. The key findings are as follows: (1) The 10 nm and 15 nm NSLDs embed in the PS monolayer through the major sphere of spherical-type, major ellipsoid of ellipsoidal-type, or one edge of cubic-type NSLDs upon contact the PS monolayer. (2) Adsorbed NSLDs cause a higher Sz value (ASz > 0.84), while embedded NSLDs cause a lower Sz value (0.47 < ASz < 0.83) that decreases with an increase in the number of bulges. (3) The embedding process absorbs 50-342 dipalmitoylphosphatidylcholine (DPPC) molecules, reducing the PS monolayer area by 0.21%-6.05%. NSLDs with bulges absorb approximately 9-126 additional DPPC molecules and cause a 0.05%-3.22% reduction in the PS monolayer area compared to NSLDs without bulges. (4) NSLDs move obliquely or vertically within the PS monolayer, displaying two distinct stages with varying velocities. Their movement direction and speed are influenced by the increasing complexity of NSLD with more bulges on them. In general, larger NSLDs with sharper shapes and increasing complex morphology of more bulges cause more significant damages to the PS monolayer. These findings have implications for safeguarding astronaut health in future manned lunar exploration missions.

Anomalous lateral diffusion of lipids during the fluid/gel phase transition of a lipid membrane

A lipid membrane undergoes a phase transition from fluid to gel phase upon changing external thermodynamic conditions, such as decreasing temperature and increasing pressure. Extremophilic organisms face the challenge of preventing this deleterious phase transition. The main focus of their adaptive strategy is to facilitate effective temperature sensing through sensor proteins, relying on the drastic changes in packing density and membrane fluidity during the phase transition. Although the changes in packing density parameters due to the fluid/gel phase transition are studied in detail, the impact on membrane fluidity is less explored in the literature. Understanding the lateral diffusive dynamics of lipids in response to temperature, particularly during the fluid/gel phase transition, is albeit crucial. Here we have simulated the phase transition of a single component lipid membrane composed of dipalmitoylphosphatidylcholine (DPPC) lipids using a coarse-grained (CG) model and studied the changes of the structural and dynamical properties. It is observed that near the phase transition point, both fluid and gel phase domains coexist together. The dynamics remains highly non-Gaussian for a long time even when the mean square displacement reaches the Fickian regime at a much earlier time. This Fickian yet non-Gaussian diffusion (FnGD) is a characteristic of a highly heterogeneous system, previously observed for the lateral diffusion of lipids in raft mimetic membranes having liquid-ordered and liquid-disordered phases co-existing together. We have analyzed the molecular trajectories and calculated the jump-diffusion of the lipids, stemming from sudden jump translations, using a translational jump-diffusion (TJD) approach. An overwhelming contribution of the jump-diffusion of the lipids is observed suggesting anomalous diffusion of lipids during fluid/gel phase transition of the membrane. These results are important in unravelling the intricate nature of lipid diffusion during the phase transition of the membrane and open up a new possibility of investigating the most significant change of membrane properties during phase transition, which can be effectively sensed by proteins.

Detecting dynamic domains and local fluctuations in complex molecular systems via timelapse neighbors shuffling

It is known that the behavior of many complex systems is controlled by local dynamic rearrangements or fluctuations occurring within them. Complex molecular systems, composed of many molecules interacting with each other in a Brownian storm, make no exception. Despite the rise of machine learning and of sophisticated structural descriptors, detecting local fluctuations and collective transitions in complex dynamic ensembles remains often difficult. Here, we show a machine learning framework based on a descriptor which we name Local Environments and Neighbors Shuffling (LENS), that allows identifying dynamic domains and detecting local fluctuations in a variety of systems in an abstract and efficient way. By tracking how much the microscopic surrounding of each molecular unit changes over time in terms of neighbor individuals, LENS allows characterizing the global (macroscopic) dynamics of molecular systems in phase transition, phases-coexistence, as well as intrinsically characterized by local fluctuations (e.g., defects). Statistical analysis of the LENS time series data extracted from molecular dynamics trajectories of, for example, liquid-like, solid-like, or dynamically diverse complex molecular systems allows tracking in an efficient way the presence of different dynamic domains and of local fluctuations emerging within them. The approach is found robust, versatile, and applicable independently of the features of the system and simply provided that a trajectory containing information on the relative motion of the interacting units is available. We envisage that "such a LENS" will constitute a precious basis for exploring the dynamic complexity of a variety of systems and, given its abstract definition, not necessarily of molecular ones.

Lipid-Coated Nanobubbles in Plants

One of the more surprising occurrences of bulk nanobubbles is in the sap inside the vascular transport system of flowering plants, the xylem. In plants, nanobubbles are subjected to negative pressure in the water and to large pressure fluctuations, sometimes encompassing pressure changes of several MPa over the course of a single day, as well as wide temperature fluctuations. Here, we review the evidence for nanobubbles in plants and for polar lipids that coat them, allowing nanobubbles to persist in this dynamic environment. The review addresses how the dynamic surface tension of polar lipid monolayers allows nanobubbles to avoid dissolution or unstable expansion under negative liquid pressure. In addition, we discuss theoretical considerations about the formation of lipid-coated nanobubbles in plants from gas-filled spaces in the xylem and the role of mesoporous fibrous pit membranes between xylem conduits in creating the bubbles, driven by the pressure gradient between the gas and liquid phase. We discuss the role of surface charges in preventing nanobubble coalescence, and conclude by addressing a number of open questions about nanobubbles in plants.

An in silico osmotic pressure approach allows characterization of pressure-area isotherms of lipid monolayers at low molecular areas

Surface pressure-area isotherms of lipid monolayers at the air-water interface provide essential information about the structure and mechanical behaviour of lipid membranes. These curves can be readily obtained through Langmuir trough measurements and, as such, have been collected for decades in the field of membrane biochemistry. However, it is still challenging to directly observe and understand nanoscopic features of monolayers through such experiments, and molecular dynamics (MD) simulations are generally used to provide a molecular view of such interfaces. In MD simulations, the surface pressure-area (Π-A) isotherms are generally computed using the Kirkwood-Irving formula, that relies on the evaluation of the pressure tensor. This approach, however, has intrinsic limitations when the molecular area in the monolayer is low (typically < 60 Ų per lipid). Recently, an alternative method to compute Π-A isotherms of surfactants, based on the calculation of the three-dimensional osmotic pressure via the implementation of semipermeable barriers was proposed. In this work, we investigate the feasibility of this approach for long-chain surfactants such as phospholipids. We identify some discrepancies between the computed values and experimental results, and we propose a semi-empirical correction based on the molecular structure of the surfactants at the monolayer interface. To validate the potential of this new approach, we simulate several phosphatidylcholine and phosphatidylethanolamine lipids at various temperatures using all-atom and coarse-grained force fields, and we compute the corresponding Π-A isotherms. Our results show that the Π-A isotherms obtained using the new method are in very good agreement with experiments and far superior to the canonical pressure tensor-based method at low molecular areas. This corrected osmotic pressure method allows for accurate characterization of the molecular packing in monolayers in various physical phases.

Effect of naproxen on the model lipid membrane formed on the water-chitosan subphase

Non steroidal anti-inflammatory drugs (NSAIDs) are those of the most common over the counter (OTC) medications widely used by millions of people every day. Unfortunately, despite their popularity those drugs can cause serious side effects in the digestive system (ulcers, bleeding, and pain). These inconveniences are caused by the changes in the structures of the outer phospholipid layers of gastric mucus and mucosa. As a result the H⁺ ions from the stomach acid can pass easily through these natural protective barriers and damage the epithelial cells which causes ulcers and bleeding. Chitosan as a polysaccharide known for its unique biocompatibility, drug delivery possibilities and wound healing effect has been chosen to examine if it can induce the reduction of undesirable effects of naproxen. This paper focuses on the interactions of the naproxen with a model biological membrane with and without the presence of chitosan. Applying the Langmuir technique coupled with the surface potential measurements and the Brewster angle microscope imaging allowed to characterize successfully examined systems in terms of the monolayer compressibility, thickness, stability, electric properties and morphology. The results proved that the presence of naproxen alters the mechanical and electrical properties of the model membrane depending on its surface pressure. Moreover, the addition of chitosan to the lipid-drug system causes significant changes in the properties of the layer, i.e. a reduction of its compressibility, thickness and morphology modification. Nevertheless, chitosan suppresses some changes induced by naproxen such as alteration of the apparent dipole moment and film stability.

The concentration-dependent effect of hydrocortisone on the structure of model lung surfactant monolayer by using an in silico approach

Understanding the adsorption mechanism of corticosteroids in the lung surfactant requires the knowledge of corticosteroid molecular interactions with lung surfactant monolayer (LSM). We employed coarse-grained molecular dynamics simulation to explore the action of hydrocortisone on an LSM comprised of a phospholipid, cholesterol and surfactant protein. The structural and dynamical morphology of the lung surfactant monolayer at different surface tensions were investigated to assess the monolayer compressibility. The simulations were also conducted at the two extreme ends of breathing cycles: exhalation (0 mN m^-1 surface tension) and inhalation (20 mN m^-1 surface tension). The impact of surface tension and hydrocortisone concentration on the monolayer compressibility and stability are significant, resulting the monolayer expansion at higher surface tension. However, at low surface tension, the highly compressed monolayer induces monolayer instability in the presence of the drug due to the accumulation of surfactant protein and drug. The constant area per lipid simulation results demonstrate that the surface pressure-area isotherms show a decrease in area-per-lipid with increased drug concentration. The drug-induced expansion causes considerable instability in the monolayer after a specific drug concentration is attained at inhalation breathing condition, whereas, for exhalation breathing, the monolayer gets more compressed, causing the LSM to collapse. The monolayer collapse occurs for inhalation due to the higher drug concentration, whereas for exhalation due to the accumulation of surfactant proteins and drugs. The findings from this study will aid in enhancing the knowledge of molecular interactions of corticosteroid drugs with lung surfactants to treat respiratory diseases.

Homeoviscous Adaptation of the Lipid Membrane of a Soil Bacterium Surviving under Diurnal Temperature Variation: A Molecular Simulation Perspective

A recent experiment has reported the lipidome remodeling of a soil-based plant-associated bacterium Methylobacterium extorquens due to diurnal temperature variations. The key adaptation strategy is the headgroup-specific remodeling of the acyl chain. To understand the idiosyncratic adaptation at the molecular level, we simulate the model membrane of the same bacterium using the reported lipidome compositions at four different experimental temperatures. We investigate the temperature-dependent packing density and fluidity of the membrane, the constancy of which is key to the homeoviscous adaptation. The results show that complex lipidome remodeling approximately preserves membrane properties under heat and cold stress. The headgroup-specific remodeling of the acyl chain serves to fine-tune the packing density and fluidity of the membrane at different temperatures. While lipids with strongly interacting headgroups are more abundant at higher temperatures, the lipidome is more dominated by lipids with weaker interacting headgroups at lower temperatures. This adaptation alleviates lipid membrane disruption caused by heat and cold stress. This study provides a molecular picture of the homeoviscous adaptation of the realistic lipid membrane of a soil-based bacterium.

Fluid-gel coexistence in lipid membranes under differential stress

A widely conserved property of many biological lipid bilayers is their asymmetry. In addition to having distinct compositions on its two sides, a membrane can also exhibit different tensions in its two leaflets, a state known as differential stress. Here, we examine how this stress can influence the phase behavior of the constituent lipid monolayers of a single-component membrane. For temperatures sufficiently close to, but still above, the main transition, molecular dynamics simulations show the emergence of finite gel domains within the compressed leaflet. We describe the thermodynamics of this phenomenon by adding two empirical single-leaflet free energies for the fluid-gel transition, each evaluated at its respective asymmetry-dependent lipid density. Finite size effects arising in simulation are included in the theory through a geometry-dependent interfacial term. Our model reproduces the phase coexistence observed in simulation. It could therefore be used to connect the "hidden variable" of differential stress to experimentally observable properties of the main phase transition. These ideas could be generalized to any first-order bilayer phase transition in the presence of asymmetry, including liquid-ordered/liquid-disordered phase separation.

The lung surfactant activity probed with molecular dynamics simulations

The surface of pulmonary alveolar subphase is covered with a mixture of lipids and proteins. This lung surfactant plays a crucial role in lung functioning. It shows a complex phase behavior which can be altered by the interaction with third molecules such as drugs or pollutants. For studying multicomponent biological systems, it is of interest to couple experimental approach with computational modelling yielding atomic-scale information. Simple two, three, or four-component model systems showed to be useful for getting more insight in the interaction between lipids, lipids and proteins or lipids and proteins with drugs and impurities. These systems were studied theoretically using molecular dynamic simulations and experimentally by means of the Langmuir technique. A better understanding of the structure and behavior of lung surfactants obtained from this research is relevant for developing new synthetic surfactants for efficient therapies, and may contribute to public health protection.

Drug Meets Monolayer: Understanding the Interactions of Sterol Drugs with Models of the Lung Surfactant Monolayer Using Molecular Dynamics Simulations

The lung surfactant monolayer (LSM) is a thin layer of lipids and proteins that forms the air/water interface of the alveoli. The primary function of the LSM is to reduce the surface tension at the air/water interface during breathing. The LSM also forms the main biological barrier for any inhaled particles, including drugs, to treat lung diseases. Elucidating the mechanism by which these drugs bind to and absorb into the LSM requires a molecular-level understanding of any drug-induced changes to the morphology, structure, and phase changes of the LSM.Molecular dynamics simulations have been used extensively to study the structure and dynamics of the LSM. The monolayer is usually simulated in at least two states: the compressed state, mimicking exhalation, and the expanded state, mimicking inhalation. In this chapter, we provide detailed instructions on how to set up, run, and analyze coarse-grained MD simulations to study the concentration-dependent effect of a sterol drug on the LSM, both in the expanded and compressed state.

Competitive and/or cooperative interactions of graphene-family materials and benzo[a]pyrene with pulmonary surfactant: a computational and experimental study

BACKGROUND

Airborne nanoparticles can be inhaled and deposit in human alveoli, where pulmonary surfactant (PS) molecules lining at the alveolar air-water interface act as the first barrier against inhaled nanoparticles entering the body. Although considerable efforts have been devoted to elucidate the mechanisms underlying nanoparticle-PS interactions, our understanding on this important issue is limited due to the high complexity of the atmosphere, in which nanoparticles are believed to experience transformations that remarkably change the nanoparticles' surface properties and states. By contrast with bare nanoparticles that have been extensively studied, relatively little is known about the interactions between PS and inhaled nanoparticles which already adsorb contaminants. In this combined experimental and computational effort, we investigate the joint interactions between PS and graphene-family materials (GFMs) with coexisting benzo[a]pyrene (BaP).

RESULTS

Depending on the BaP concentration, molecular agglomeration, and graphene oxidation, different nanocomposite structures are formed via BaPs adsorption on GFMs. Upon deposition of GFMs carrying BaPs at the pulmonary surfactant (PS) layer, competition and cooperation of interactions between different components determines the interfacial processes including BaP solubilization, GFM translocation and PS perturbation. Importantly, BaPs adsorbed on GFMs are solubilized to increase BaP's bioavailability. By contrast with graphene adhering on the PS layer to release part of adsorbed BaPs, more BaPs are released from graphene oxide, which induces a hydrophilic pore in the PS layer and shows adverse effect on the PS biophysical function. Translocation of graphene across the PS layer is facilitated by BaP adsorption through segregating it from contact with PS, while translocation of graphene oxide is suppressed by BaP adsorption due to the increase of surface hydrophobicity. Graphene extracts PS molecules from the layer, and the resultant PS depletion declines with graphene oxidation and BaP adsorption.

CONCLUSION

GFMs showed high adsorption capacity towards BaPs to form nanocomposites. Upon deposition of GFMs carrying BaPs at the alveolar air-water interface covered by a thin PS layer, the interactions of GFM-PS, GFM-BaP and BaP-PS determined the interfacial processes of BaP solubilization, GFM translocation and PS perturbation.

The steroid mometasone alters protein containing lung surfactant monolayers in a concentration-dependent manner

Mometasone is an investigational anti-inflammatory steroidal drug to treat inflammation via pulmonary administration. For steroid drugs to be effective they need to be adsorbed by lung surfactants, a thin monolayer at the air-water interface in alveoli that reduces surface tension. Information on the molecular-level interactions of the drug with lung surfactants is useful to understand the mechanism of adsorption. In this study, we use coarse-grained molecular dynamics simulation to understand the concentration-dependent effect of mometasone on a lung surfactant monolayer (LSM) composed of lipids and surfactant proteins, under two different breathing conditions (exhalation, at surface tension 0 mNm^-1 and inhalation, surface tension 20-25 mNm^-1). A series of fixed-APL and fixed-surface tension simulations were used to demonstrate that in the absence of drugs, the model LSM reproduces the surface tensions for the compressed and expanded states, as well as compressibility at different surface tensions. In-depth analysis of simulations of a LSM in the presence of five different drug concentrations shows that mometasone alters the structure and dynamics of the LSM in a concentration-dependent manner. Mometasone induces a collapse in the monolayer that is affected by the surfactant protein and surface tension. Overall, these findings suggest that the surfactant proteins, surface tension and drug concentration are all critical components affecting monolayer stability and drug adsorption. The outcomes of this study may be beneficial for a more in-depth understanding of how mometasone is adsorbed by lung surfactants.

Interfacial behavior of phospholipid monolayers revealed by mesoscopic simulation

A mesoscopic model with molecular resolution is presented for dipalmitoyl phosphatidylcholine (DPPC) and palmitoyl oleoyl phosphatidylcholine (POPC) monolayer simulations at the air-water interface using many-body dissipative particle dynamics (MDPD). The parameterization scheme is rigorously based on reproducing the physical properties of water and alkane and the interfacial property of the phospholipid monolayer by comparison with experimental results. Using much less computing cost, these MDPD simulations yield a similar surface pressure-area isotherm as well as similar pressure-related morphologies as all-atom simulations and experiments. Moreover, the compressibility modulus, order parameter of lipid tails, and thickness of the phospholipid monolayer are quantitatively in line with the all-atom simulations and experiments. This model also captures the sensitive changes in the pressure-area isotherms of mixed DPPC/POPC monolayers with altered mixing ratios, indicating that the model is promising for applications with complex natural phospholipid monolayers. These results demonstrate a significant improvement of quantitative phospholipid monolayer simulations over previous coarse-grained models.

Compositional, structural and functional properties of discrete coexisting complexes within bronchoalveolar pulmonary surfactant

Lung surfactant (LS) stabilizes the respiratory surface by forming a film at the alveolar air-liquid interface that reduces surface tension and minimizes the work of breathing. Typically, this surface-active agent has been isolated from animal lungs both for research and biomedical applications. However, these materials are constituted by complex membranous architectures including surface-active and inactive lipid/protein assemblies. In this work, we describe the composition, structure and surface activity of discrete membranous entities that are part of a LS preparation isolated from bronchoalveolar lavages of porcine lungs. Seven different fractions could be resolved from whole surfactant subjected to sucrose density gradient centrifugation. Detailed compositional characterization revealed differences in protein and cholesterol content but no distinct saturated:unsaturated phosphatidylcholine ratios. Moreover, no significant differences were detected regarding apparent hydration at the headgroup region of membranes, as reported by the probe Laurdan, and lipid chain mobility analysed by electron spin resonance (ESR) in spite of the variety of membranous assemblies observed by transmission electron microscopy. In addition, six of the seven separated LS subfractions formed similar, essentially disordered-like, interfacial films and performed efficient surface activity, under physiologically relevant conditions. Altogether, our work show that a LS isolated from porcine lungs is comprised by a heterogenous population of membranous assemblies lacking freshly secreted unused LS complexes sustaining highly dehydrated and ordered membranous assemblies as previously reported. We propose that surfactant subfractions may illustrate intermediates in sequential structural steps within the structural transformations occurring along the respiratory compression-expansion cycles.

Shape matters-the interaction of gold nanoparticles with model lung surfactant monolayers

The lung surfactant monolayer (LSM) forms the main biological barrier for any inhaled particles to enter our bloodstream, including gold nanoparticles (AuNPs) present as air pollutants and under investigation for use in biomedical applications. Understanding the interaction of AuNPs with lung surfactant can assist in understanding how AuNPs enter our lungs. In this study, we use coarse-grained molecular dynamics simulations to investigate the effect of four different shape D AuNPs (spherical, box, icosahedron and rod) on the structure and dynamics of a model LSM, with a particular focus on differences resulting from the shape of the AuNP. Monolayer-AuNP systems were simulated in two different states: the compressed state and the expanded state, representing inhalation and exhalation conditions, respectively. Our results indicate that the compressed state is more affected by the presence of the AuNPs than the expanded state. Our results show that in the compressed state, the AuNPs prevent the monolayer from reaching the close to zero surface tension required for normal exhalation. In the compressed state, all four nanoparticles (NPs) reduce the lipid order parameters and cause a thinning of the monolayer where the particles drag surfactant molecules into the water phase. Comparing the different properties shows no trend concerning which shape has the biggest effect on the monolayer, as shape-dependent effects vary among the different properties. Insights from this study might assist future work of how AuNP shapes affect the LSM during inhalation or exhalation conditions.

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.

A Data-Driven Dimensionality Reduction Approach to Compare and Classify Lipid Force Fields

Molecular dynamics simulations of all-atom and coarse-grained lipid bilayer models are increasingly used to obtain useful insights for understanding the structural dynamics of these assemblies. In this context, one crucial point concerns the comparison of the performance and accuracy of classical force fields (FFs), which sometimes remains elusive. To date, the assessments performed on different classical potentials are mostly based on the comparison with experimental observables, which typically regard average properties. However, local differences of the structure and dynamics, which are poorly captured by average measurements, can make a difference, but these are nontrivial to catch. Here, we propose an agnostic way to compare different FFs at different resolutions (atomistic, united-atom, and coarse-grained), by means of a high-dimensional similarity metrics built on the framework of Smooth Overlap of Atomic Position (SOAP). We compare and classify a set of 13 FFs, modeling 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayers. Our SOAP kernel-based metrics allows us to compare, discriminate, and correlate different FFs at different model resolutions in an unbiased, high-dimensional way. This also captures differences between FFs in modeling nonaverage events (originating from local transitions), for example, the liquid-to-gel phase transition in dipalmitoylphosphatidylcholine (DPPC) bilayers, for which our metrics allows us to identify nucleation centers for the phase transition, highlighting some intrinsic resolution limitations in implicit versus explicit solvent FFs.

Monolayers of Cholesterol and Cholesteryl Stearate at the Water/Vapor Interface: A Physico-Chemical Study of Components of the Meibum Layer

Langmuir monolayers containing different amounts of cholesterol and cholesteryl stearate were studied at two different temperatures (24 °C and 35 °C). The main goal was to contribute towards the understanding of how the variations in the chemical composition may affect the physico-chemical properties of these specific lipid monolayers. The model mixture was chosen considering that cholesteryl esters are present in cell membranes and some other biological systems, including human tear lipids. Therefore, an investigation into the effect of the lipid monolayer composition on their interfacial properties may elucidate some of the fundamental reasons for the deficiencies in cell membranes and tear film functioning in vivo. The experimental results have shown that the molar ratio of the mixture plays a crucial role in the modulation of the Langmuir film properties. The condensing effects of the cholesterol and the interactions between the lipids in the monolayer were the main factors altering the monolayer response to dilatational deformation. The modification of the mixture compositions leads to significant changes in the Langmuir films and the mechanical performance, altering the ability of the monolayer to reduce the surface tension and the viscoelastic properties of the monolayers. This suggests that subtle modifications of the biomembrane composition may significantly alter its physiological function.

Structural hallmarks of lung surfactant: Lipid-protein interactions, membrane structure and future challenges

Lung surfactant (LS) is an outstanding example of how a highly regulated and dynamic membrane-based system has evolved to sustain a wealth of structural reorganizations in order to accomplish its biophysical function, as it coats and stabilizes the respiratory air-liquid interface in the mammalian lung. The present review dissects the complexity of the structure-function relationships in LS through an updated description of the lipid-protein interactions and the membrane structures that sustain its synthesis, secretion, interfacial performance and recycling. We also revise the current models and the biophysical techniques employed to study the membranous architecture of LS. It is important to consider that the structure and functional properties of LS are often studied in bulk or under static conditions, in spite that surfactant function is strongly connected with a highly dynamic behaviour, sustained by very polymorphic structures and lipid-lipid, lipid-protein and protein-protein interactions that reorganize in precise spatio-temporal coordinates. We have tried to underline the evidences available of the existence of such structural dynamism in LS. A last important aspect is that the synthesis and assembly of LS is a strongly regulated intracellular process to ensure the establishment of the proper interactions driving LS surface activity, while protecting the integrity of other cell membranes. The use of simplified lipid models or partial natural materials purified from animal tissues could be too simplistic to understand the true molecular mechanisms defining surfactant function in vivo. In this line, we will bring into the attention of the reader the methodological challenges and the questions still open to understand the structure-function relationships of LS at its full biological relevance.

Computational Studies of Lipid-Wrapped Gold Nanoparticle Transport Through Model Lung Surfactant Monolayers

Colloidal nanoparticles, such as gold nanoparticles (AuNPs), are promising materials for the delivery of hydrophilic drugs via the pulmonary route. The inhaled nanoparticle drug carriers primarily deposit in lung alveoli and interact with the alveolar surface known as lung surfactants. Therefore, it is vital to understand the interactions of nanocarriers with the surfactant layer. To understand the interactions at the molecular level, here we simulated model lung surfactant monolayers with phospholipid (PL)-wrapped AuNPs at the vacuum-water interface using coarse-grained molecular dynamics simulations. The PL-wrapped AuNPs quickly adsorbed into the surfactant layer, altered the structural properties of the monolayer, and at high concentrations initiated the compressed monolayer to collapse/buckle. Among the surfactant monolayer lipid components, cholesterol adsorbed to the AuNPs preferentially over PL species. The position of the adsorbed PL-AuNPs within the monolayer, and subsequent monolayer perturbation, vary depending on the monolayer phase, monolayer composition, and species of PL used as a ligand. Information provided by these molecular dynamic simulations helps to rationalize why some colloidal nanoparticles work better as nanocarriers than others and aid the design of new ones, to avoid biological toxicity and improve efficacy for pulmonary drug delivery.

Understanding the Functional Properties of Lipid Heterogeneity in Pulmonary Surfactant Monolayers at the Atomistic Level

Pulmonary surfactant is a complex mixture of lipids and proteins lining the interior of the alveoli, and constitutes the first barrier to both oxygen and pathogens as they progress toward blood circulation. Despite decades of study, the behavior of the pulmonary surfactant at the molecular scale is poorly understood, which hinders the development of effective surfactant replacement therapies, useful in the treatment of several lung-related diseases. In this work, we combined all-atom molecular dynamics simulations, Langmuir trough measurements, and AFM imaging to study synthetic four-component lipid monolayers designed to model protein-free pulmonary surfactant. We characterized the structural and dynamic properties of the monolayers with a special focus on lateral heterogeneity. Remarkably, simulations reproduce almost quantitatively the experimental data on pressure-area isotherms and the presence of lateral heterogeneities highlighted by AFM. Quite surprisingly, the pressure-area isotherms do not show a plateau region, despite the presence of liquid-condensed nanometer-sized domains at surface pressures larger than 20 mN/m. In the simulations, the liquid-condensed domains were small and transient, but they did not coalesce to yield a separate phase. They were only slightly enriched in DPPC and cholesterol, and their chemical composition remained very similar to the overall composition of the monolayer membrane. Instead, they differed from liquid-expanded regions in terms of membrane thickness (in agreement with AFM data), diffusion rates, as well as acyl chain packing and orientation. We hypothesize that such lateral heterogeneities are crucial for lung surfactant function, as they allow both efficient packing, to achieve low surface tension, and sufficient fluidity, critical for rapid adsorption to the air–liquid interface during the breathing cycle.

Modeling Gas-Liquid Interfaces by Dissipative Particle Dynamics: Adsorption and Surface Tension of Cetyl Trimethyl Ammonium Bromide at the Air-Water Interface

Adsorption of surfactants at gas-liquid interfaces that causes reduction in the surface tension is a classical problem in colloid and interface science with multiple practical applications in oil and gas recovery, separations, cosmetics, personal care, and biomedicine. Here, we develop an original coarse-grained model of the liquid-gas interface within the conventional dissipative particle dynamics (DPD) framework with the goal of quantitatively predicting the surface tension in the presence of surfactants. As a practical case-study example, we explore the adsorption of the cationic surfactant cetyl trimethyl ammonium bromide (CTAB) on the air-water interface. The gas phase is modeled as a DPD fluid composed of fictitious hard-core "gas" beads with exponentially decaying repulsive potentials to prevent penetration of the liquid phase components. A rigorous parametrization scheme is proposed based on matching the bulk and interfacial properties of water and octane taken as the reference compounds. Quantitative agreement between the simulated and experimental surface tension of CTAB solutions is found for a wide range of bulk surfactant concentrations (∼10^-3 to ∼1 mmol/L) with the reduction of the surface tension from ∼72 mN/m (pure water) to the limiting value of ∼37.5 mN/m at the critical micelle concentration. The gas phase DPD model with the proposed parametrization scheme can be extended and applied to modeling various gas-liquid interfaces with surfactant and lipid monolayers, such as bubble suspensions, foams, froths, etc.

Preferential Protein Partitioning in Biological Membrane with Coexisting Liquid Ordered and Liquid Disordered Phase Behavior: Underlying Design Principles

Several studies now show that certain proteins exhibit selective preference toward liquid ordered (L[Formula: see text]) or toward liquid disordered (L[Formula: see text]) regions of the heterogeneous membrane and some of them have preference for the L[Formula: see text]-L[Formula: see text] interface. Spatially heterogenous organization of lipids, enriched in specific protein molecules, function as platforms for signaling and are involved in several other physiologically critical functions. In this review, we collate together some of the experimental observations of cases where proteins preferentially segregate into different phases and highlight the importance of these preferential localization in terms of underlying functions. We also try to understand the structural features and chemical makeup of the membrane-interacting motifs of these proteins. Finally, we put forth some preliminary analysis on class I viral fusion proteins, some of which are known to partition at the L[Formula: see text]-L[Formula: see text] interface, and through them we try to understand the evolutionary design principles of phase segregating proteins. Put together, this review summarizes the existing studies on preferential partitioning of proteins into different membrane phases while emphasizing the need to understand the molecular design-level features that can help us "engineer" functionally rich peptides and proteins with a programmed membrane partitioning.

Molecular mechanism for bidirectional regulation of CD44 for lipid raft affiliation by palmitoylations and PIP2

The co-localization of Cluster-of-Differentiation-44 protein (CD44) and cytoplasmic adaptors in specific membrane environments is crucial for cell adhesion and migration. The process is controlled by two different pathways: On the one hand palmitoylation keeps CD44 in lipid raft domains and disables the linking to the cytoplasmic adaptor, whereas on the other hand, the presence of phosphatidylinositol-4,5-biphosphate (PIP2) lipids accelerates the formation of the CD44-adaptor complex. The molecular mechanism explaining how CD44 is migrating into and out of the lipid raft domains and its dependence on both palmitoylations and the presence of PIP2 remains, however, elusive. In this study, we performed extensive molecular dynamics simulations to study the raft affinity and translocation of CD44 in phase separated model membranes as well as more realistic plasma membrane environments. We observe a delicate balance between the influence of the palmitoylations and the presence of PIP2 lipids: whereas the palmitoylations of CD44 increases the affinity for raft domains, PIP2 lipids have the opposite effect. Additionally, we studied the association between CD44 and the membrane adaptor FERM in dependence of these factors. We find that the presence of PIP2 lipids allows CD44 and FERM to associate in an experimentally observed binding mode whereas the highly palmitoylated species shows no binding affinity. Together, our results shed light on the sophisticated mechanism on how membrane translocation and peripheral protein association can be controlled by both protein modifications and membrane composition.

Cytoskeleton protein complex involving with association of CD44 and ERMs is critical for cancer-related cellular adhesion and migration. The protein interactions are found to be modulated by chemical modification and membrane microenvironments, but the inherent mechanism is unclear. We obtained molecular dynamic details of CD44 localization switching between raft/non-raft subdomains regulated by palmitoylations and PIP2 molecules. Binding of PIP2 on the palmitoylated CD44 enables it to release from lipid raft, revealing an exceptional role of PIP2 in mediating protein translocation. PIP2 is beneficial for CD44 to associate with the active domain of ERM, in a nearly crystal structure mode. The molecular information will enhance our understanding for PIP2 regulation to protein translocation and membrane association.

The aggregation of striped nanoparticles in mixed phospholipid bilayers

The unique and adjustable properties of nanoparticles present enormous opportunities for their use as targeted drug delivery vectors. For example, nanoparticles functionalized with key surface ligands have been shown to pass through phospholipid bilayers without causing localised disruption. However, the further effects nanoparticles have on multi-component phospholipid bilayers remain unclear. We use coarse-grained computational models to investigate the structural properties of mixed phospholipid bilayers in the presence of ligand-functionalized nanoparticles. Model bilayers are composed of DPPC, DUPC, DFPC and cholesterol, and the nanoparticles are striped with a hydrophobic-ligand band and charged-ligand spherical caps. Our results show that nanoparticles aggregate near unsaturated phospholipid regions, phospholipid bilayer phase-separation is promoted in the presence of nanoparticles, and the heterogeneous components of a phospholipid bilayer play a significant role in the lateral organization of nanoparticles. This study highlights the need for considering the complexity of realistic phospholipid bilayers when optimising ligand functionalized nanoparticles for efficient drug delivery vectors.

Lung surfactant monolayer - A good natural barrier against dibenzo-p-dioxins

The present study deals with interaction of two air pollutants: dibenzodioxin, DD, and its' monochlorinated derivative, 2-chlorodibenzodioxin, 2CLDD, with models of the lung surfactant (LS) system. A monolayer composed of DPPC and POPC in 1:1 molar ratio was used as a model of LS. One component monolayers of DPPC and POPC were also examined, to model the interiors of LC and LE domains in LS, respectively. Molecular dynamics simulations and measurements of surface pressure isotherms, as well as polarization modulation-infrared reflection-absorption spectra were employed to study the influence of dioxins on the monolayers. We demonstrate, that both dioxins adsorb and accumulate in the hydrophobic parts of all three monolayers. DD molecules prefer flat orientation on the surface at large areas. Upon compression, they lift and orient perpendicularly to the monolayer. Flat orientation of DD molecules leads to their large surface area. In consequence they preferentially locate in vicinity of unsaturated chains of POPC - they are small enough to fill void spaces created by kinks in unsaturated chains. 2CLDD orient along monolayer normal already at the largest areas and preference for POPC was not observed for them. In laterally relaxed states, a condensing effect, connected with reduction of surface area available to the lipids was observed for both dioxins. In the case of 2CLDD, additional locally ordering influence of dioxin molecules was detected. In compressed states, the presence of dioxin molecules hinders alignment and uniform ordering of lipid chains.

Molecular Ordering in Lipid Monolayers: An Atomistic Simulation

We report on atomistic simulations of DPPC lipid monolayers using the CHARMM36 lipid force field (and also the Slipid force field as a control case), combined with a four-point OPC water model. The entire two-phase region where domains of the "liquid-condensed" (LC) phase coexist with domains of the "liquid-expanded" (LE) phase has been explored. The simulations are long enough that the complete phase-transition stage, with two domains coexisting in the monolayer, is reached in all cases. Also, system sizes used are larger than those in previous works. As expected, domains of the minority phase are elongated, emphasizing the importance of anisotropic van der Waals and/or electrostatic dipolar interactions in the monolayer plane. The molecular structure is quantified in terms of distribution functions for the hydrocarbon chains and the PN dipoles. In contrast to previous work, where average distributions are calculated, distributions are here extracted for each of the coexisting phases by first identifying lipid molecules that belong to either LC or LE regions. In the case of the CHARMM36 force field, the three-dimensional distributions show that the average tilt angle of the chains with respect to the normal outward direction is (39.0 ± 0.1)° in the LC phase and (48.1 ± 0.5)° in the LC phase. In the case of the PN dipoles, the distributions indicate a tilt angle of (110.8 ± 0.5)° in the LC phase and (112.5 ± 0.5)° in the LE phase. These results are quantitatively different from those in previous works, which indicated a smaller normal component of the PN dipole. Also, the distributions of the monolayer-projected chains and PN dipoles have been calculated. Chain distributions peak along a particular direction in the LC domains, while they are uniform in the LE phase. Long-range ordering associated with the projected PN dipoles is absent in both phases. These results strongly suggest that LC domains do not exhibit dipolar ordering in the plane of the monolayer, the effect of these components being averaged out at short distances. Therefore, the only relevant component of the molecular dipoles, with regard to both intra- and long-range interdomain interactions, is normal to the monolayer. Also, the local orientation of chain projections is almost constant in LC domains and points in the direction along which domains are elongated, suggesting that the line tension driving the phase transition might be anisotropic with respect to the interfacial domain boundary.

Surfactant Charge Modulates Structure and Stability of Lipase-Embedded Monolayers at Marine-Relevant Aerosol Surfaces

Lipases, as well as other enzymes, are present and active within the sea surface microlayer (SSML). Upon bubble bursting, lipases partition into sea spray aerosol (SSA) along with surface-active molecules such as lipids. Lipases are likely to be embedded in the lipid monolayer at the SSA surface and thus have the potential to influence SSA interfacial structure and chemistry. Elucidating the structure of the lipid monolayer at SSA interfaces and how this structure is altered upon interaction with a protein system like lipase is of interest, given the importance of how aerosols interact with sunlight, influence cloud formation, and provide surfaces for chemical reactions. Herein, we report an integrated experimental and computational study of Burkholderia cepacia lipase (BCL) embedded in a lipid monolayer and highlight the important role of electrostatic, rather than hydrophobic, interactions as a driver for monolayer stability. Specifically, we combine Langmuir film experiments and molecular dynamics (MD) simulations to examine the detailed interactions between the zwitterionic dipalmitoylphosphatidylcholine (DPPC) monolayer and BCL. Upon insertion of BCL from the underlying subphase into the lipid monolayer, it is shown that BCL permeates and largely disorders the monolayer while strongly interacting with zwitterionic DPPC molecules, as experimentally observed by Langmuir adsorption curves and infrared reflectance absorbance spectroscopy. Explicitly solvated, all-atom MD is then used to provide insights into inter- and intramolecular interactions that drive these observations, with specific attention to the formation of salt bridges or ionic-bonding interactions. We show that after insertion into the DPPC monolayer, lipase is maintained at high surface pressures and in large BCL concentrations by forming a salt-bridge-stabilized lipase-DPPC complex. In comparison, when embedded in an anionic monolayer at low surface pressures, BCL preferentially forms intramolecular salt bridges, reducing its total favorable interactions with the surfactant and partitioning out of the monolayer shortly after injection. Overall, this study shows that the structure and dynamics of lipase-embedded SSA surfaces vary based on surface charge and pressure and that these variations have the potential to differentially modulate the properties of marine aerosols.

A Simulation Study on the Interaction Between Pollutant Nanoparticles and the Pulmonary Surfactant Monolayer

A good understanding of the mechanism of interaction between inhaled pollutant nanoparticles (NPs) and the pulmonary surfactant monolayer is useful to study the impact of fine particulate matter on human health. In this work, we established coarse-grained models of four representative NPs with different hydrophilicity properties in the air (i.e., CaSO₄, C, SiO₂, and C₆H₁₄O₂ NPs) and the pulmonary surfactant monolayer. Molecular dynamic simulations of the interaction during exhalation and inhalation breathing states were performed. The effects of NP hydrophilicity levels, NP structural properties, and cholesterol content in the monolayer on the behaviors of NP embedment or the transmembrane were analyzed by calculating the changes in potential energy, NP displacement, monolayer orderliness, and surface tension. Results showed that NPs can inhibit the ability of the monolayer to adjust surface tension. For all breathing states, the hydrophobic C NP cannot translocate across the monolayer and had the greatest influence on the structural properties of the monolayer, whereas the strongly hydrophilic SiO₂ and C₆H₁₄O₂ NPs can cross the monolayer with little impact. The semi-hydrophilic CaSO₄ NP can penetrate the monolayer only during the inhalation breathing state. The hydrophilic flaky NP shows the best penetration ability, followed by the rod-shaped NP and spherical NP in turn. An increase in cholesterol content of the monolayer led to improved orderliness and decreased fluidity of the membrane system due to enhanced intermolecular forces. Consequently, difficulty in crossing the monolayer increased for the NPs.

Interfacial ordering of tristearin induced by glycerol monooleate and PGPR: A coarse-grained molecular dynamics study

We use coarse-grained molecular dynamics simulations to study the effect of surfactant structure on the ordering of bulk tristearin at an oil-water interface. In the absence of surfactant, tristearin acyl chains are marginally aligned normal to the interface. The surfactant glycerol monooleate (GMO), a common small-molecule monoacylglycerol (MW: 357 g/mol), preferentially adsorbs to the oil-water interface, displacing more of the tristearin as its concentration increases. The tristearin that remains at the interface is closely aligned normal to the interface. Adjacent to the interface, bulk tristearin increasingly aligns with its acyl chains entwined with the GMO acyl chain, which also preferentially aligns normal to the interface. In contrast, polyglycerol polyricinoleate (PGPR), a bulkier, polymeric surfactant (MW: 1398 g/mol for a molecule with five monomers), both displaces tristearin from the interface and reduces the alignment of the molecules that remain. We suggest that the similar fatty acid moieties of GMO (oleic acid) and tristearin (stearic acid) lead to liquid-state association and alignment, the latter of which can then serve as a template onto which tristearin crystals can nucleate. Conversely, by both displacing tristearin from the interface and reducing alignment below that of the surfactant-free system, PGPR eliminates the possibility of tristearin interfacial crystallisation.

Protein Partitioning into Ordered Membrane Domains: Insights from Simulations

Cellular membranes are laterally organized into domains of distinct structures and compositions by the differential interaction affinities between various membrane lipids and proteins. A prominent example of such structures are lipid rafts, which are ordered, tightly packed domains that have been widely implicated in cellular processes. The functionality of raft domains is driven by their selective recruitment of specific membrane proteins to regulate their interactions and functions; however, there have been few general insights into the factors that determine the partitioning of membrane proteins between coexisting liquid domains. In this work, we used extensive coarse-grained and atomistic molecular dynamics simulations, potential of mean force calculations, and conceptual models to describe the partitioning dynamics and energetics of a model transmembrane domain from the linker of activation of T cells. We find that partitioning between domains is determined by an interplay between protein-lipid interactions and differential lipid packing between raft and nonraft domains. Specifically, we show that partitioning into ordered domains is promoted by preferential interactions between peptides and ordered lipids, mediated in large part by modification of the peptides by saturated fatty acids (i.e., palmitoylation). Ordered phase affinity is also promoted by elastic effects, specifically hydrophobic matching between the membrane and the peptide. Conversely, ordered domain partitioning is disfavored by the tight molecular packing of the lipids therein. The balance of these dominant drivers determines partitioning. In the case of the wild-type linker of activation of T cells transmembrane domain, these factors combine to yield enrichment of the peptide at L_o/L_d interfaces. These results define some of the general principles governing protein partitioning between coexisting membrane domains and potentially explain previous disparities among experiments and simulations across model systems.

How to Determine Lipid Interactions in Membranes from Experiment Through the Ising Model

The determination and the meaning of interactions in lipid bilayers are discussed and interpreted through the Ising model. Originally developed to understand phase transitions in ferromagnetic systems, the Ising model applies equally well to lipid bilayers. In the case of a membrane, the essence of the Ising model is that each lipid is represented by a site on a lattice and that the interaction of each site with its nearest neighbors is represented by an energy parameter ω. To calculate the thermodynamic properties of the system, such as the enthalpy, the Gibbs energy, and the heat capacity, the partition function is derived. The calculation of the configurational entropy factor in the partition function, however, requires approximations or the use of Monte Carlo (MC) simulations. Those approximations are described. Ultimately, MC simulations are used in combination with experiment to determine the interaction parameters ω in lipid bilayers. Several experimental approaches are described, which can be used to obtain interaction parameters. They include nearest-neighbor recognition, differential scanning calorimetry, and Förster resonance energy transfer. Those approaches are most powerful when used in combination of MC simulations of Ising models. Lipid membranes of different compositions are discussed, which have been studied with these approaches. They include mixtures of cholesterol, saturated (ordered) phospholipids, and unsaturated (disordered) phospholipids. The interactions between those lipid species are examined as a function of molecular properties such as the degree of unsaturation and the acyl chain length. The general rule that emerges is that interactions between different lipids are usually unfavorable. The exception is that cholesterol interacts favorably with saturated (ordered) phospholipids. However, the interaction of cholesterol with unsaturated phospholipids becomes extremely unfavorable as the degree of unsaturation increases.

Atomistic Model for Nearly Quantitative Simulations of Langmuir Monolayers

Lung surfactant and a tear film lipid layer are examples of biologically relevant macromolecular structures found at the air-water interface. Because of their complexity, they are often studied in terms of simplified lipid layers, the simplest example being a Langmuir monolayer. Given the profound biological significance of these lipid assemblies, there is a need to understand their structure and dynamics on the nanoscale, yet there are not many techniques able to provide this information. Atomistic molecular dynamics simulations would be a tool fit for this purpose; however, the simulation models suggested until now have been qualitative instead of quantitative. This limitation has mainly stemmed from the challenge to correctly describe the surface tension of water with simulation parameters compatible with other biomolecules. In this work, we show that this limitation can be overcome by using the recently introduced four-point OPC water model, whose surface tension for water is demonstrated to be quantitatively consistent with experimental data and which is also shown to be compatible with the commonly employed lipid models. We further establish that the approach of combining the OPC four-point water model with the CHARMM36 lipid force field provides nearly quantitative agreement with experiments for the surface pressure-area isotherm for POPC and DPPC monolayers, also including the experimentally observed phase coexistence in a DPPC monolayer. The simulation models reported in this work pave the way for nearly quantitative atomistic studies of lipid-rich biological structures at air-water interfaces.

Composition Fluctuations in Lipid Bilayers

Cell membranes contain multiple lipid and protein components having heterogeneous in-plane (lateral) distribution. Nanoscale rafts are believed to play an important functional role, but their phase state—domains of coexisting phases or composition fluctuations—is unknown. As a step toward understanding lateral organization of cell membranes, we investigate the difference between nanoscale domains of coexisting phases and composition fluctuations in lipid bilayers. We simulate model lipid bilayers with the MARTINI coarse-grained force field on length scales of tens of nanometers and timescales of tens of microseconds. We use a binary and a ternary mixture: a saturated and an unsaturated lipid, or a saturated lipid, an unsaturated lipid, and cholesterol, respectively. In these mixtures, the phase behavior can be tuned from a mixed state to a coexistence of a liquid-crystalline and a gel, or a liquid-ordered and a liquid-disordered phase. Transition from a two-phase to a one-phase state is achieved by raising the temperature and adding a hybrid lipid (with a saturated and an unsaturated chain). We analyze the evolution of bilayer properties along this transition: domains of two phases transform to fluctuations with local ordering and compositional demixing. Nanoscale domains and fluctuations differ in several properties, including interleaflet overlap and boundary length. Hybrid lipids show no enrichment at the boundary, but decrease the difference between the coexisting phases by ordering the disordered phase, which could explain their role in cell membranes.