Energetics of the Transmembrane Peptide Sorting by Hydrophobic Mismatch.

Hydrophobic mismatch between a lipid membrane and embedded transmembrane peptides or proteins plays a role in their lateral localization and function. Earlier studies have resolved numerous mechanisms through which the peptides and membrane proteins adapt to mismatch, yet the energetics of lateral sorting due to hydrophobic mismatch have remained elusive due to the lack of suitable computational or experimental protocols. Here, we pioneer a molecular dynamics simulation approach to study the sorting of peptides along a membrane thickness gradient. Peptides of different lengths tilt and diffuse along the membrane to eliminate mismatch with a rate directly proportional to the magnitude of mismatch. We extract the 2-dimensional free energy profiles as a function of local thickness and peptide orientation, revealing the relative contributions of sorting and tilting, and suggesting their thermally accessible regimes. Our approach can readily be applied to study other membrane systems of biological interest where hydrophobic mismatch, or membrane thickness in general, plays a role.

Diffusion Analyses along Mean and Gaussian-Curved Membranes with CurD.

Curved cellular membranes are both abundant and functionally relevant. While novel tomography approaches reveal the structural details of curved membranes, their dynamics pose an experimental challenge. Curvature especially affects the diffusion of lipids and macromolecules, yet neither experiments nor continuum models distinguish geometric effects from those caused by curvature-induced changes in membrane properties. Molecular simulations could excel here, yet despite community interest toward curved membranes, tools for their analysis are still lacking. Here, we satisfy this demand by introducing CurD, our novel and openly available implementation of the Vertex-oriented Triangle Propagation algorithm to the study of lipid diffusion along membranes with mean and/or Gaussian curvature. This approach, aided by our highly optimized implementation, computes geodetic distances significantly faster than conventional implementations of path-finding algorithms. Our tool, applied to coarse-grained simulations, allows for the first time the analysis of curvature effects on diffusion at size scales relevant to physiological processes such as endocytosis. Our analyses with different membrane geometries reveal that Gaussian curvature plays a surprisingly small role on lipid motion, whereas mean curvature; i.e., the packing of lipid headgroups largely dictates their mobility.

Exposure to Aldehyde Cherry e-Liquid Flavoring and Its Vaping Byproduct Disrupt Pulmonary Surfactant Biophysical Function.

Over the past decade, there has been a significant rise in the use of vaping devices, particularly among adolescents, raising concerns for effects on respiratory health. Pressingly, many recent vaping-related lung injuries are unexplained by current knowledge, and the overall implications of vaping for respiratory health are poorly understood. This study investigates the effect of hydrophobic vaping liquid chemicals on the pulmonary surfactant biophysical function. We focus on the commonly used flavoring benzaldehyde and its vaping byproduct, benzaldehyde propylene glycol acetal. The study involves rigorous testing of the surfactant biophysical function in Langmuir trough and constrained sessile drop surfactometer experiments with both protein-free synthetic surfactant and hydrophobic protein-containing clinical surfactant models. The study reveals that exposure to these vaping chemicals significantly interferes with the synthetic and clinical surfactant biophysical function. Further atomistic simulations reveal preferential interactions with SP-B and SP-C surfactant proteins. Additionally, data show surfactant lipid-vaping chemical interactions and suggest significant transfer of vaping chemicals to the experimental subphase, indicating a toxicological mechanism for the alveolar epithelium. Our study, therefore, reveals novel mechanisms for the inhalational toxicity of vaping. This highlights the need to reassess the safety of vaping liquids for respiratory health, particularly the use of aldehyde chemicals as vaping flavorings.

Palmitoylation modifies transmembrane adaptor protein PAG for ordered lipid environment: A molecular dynamics simulation study.

We employed all-atom MD simulations to investigate the impact of palmitoylation on the PAG transmembrane peptide within various lipid environments, including the less explored boundary region separating lipid-ordered (Lo) and lipid-disordered (Ld) membrane phases. We found that palmitoylation of the peptide reduces its impact on membrane thickness, particularly within the Lo and boundary environments. Despite their hydrophobic nature, the palmitoyl chains on the peptide did not significantly affect the hydration of the surrounding membrane. Interestingly, the boundary membrane environment was found to be especially compatible with the palmitoylated peptide, suggesting its potential for accumulation in phase boundaries. Our findings highlight the importance of understanding how palmitoylation-modified peptides behave within membranes, with crucial implications for cell signaling and membrane organization. This knowledge may also inform the optimization of lipid membrane-based drug delivery systems, by improving our understanding of how drugs and excipients can be most effectively arranged within these carriers.

Molecular view of ER membrane remodeling by the Sec61/TRAP translocon.

Protein translocation across the endoplasmic reticulum (ER) membrane is an essential step during protein entry into the secretory pathway. The conserved Sec61 protein-conducting channel facilitates polypeptide translocation and coordinates cotranslational polypeptide-processing events. In cells, the majority of Sec61 is stably associated with a heterotetrameric membrane protein complex, the translocon-associated protein complex (TRAP), yet the mechanism by which TRAP assists in polypeptide translocation remains unknown. Here, we present the structure of the core Sec61/TRAP complex bound to a mammalian ribosome by cryogenic electron microscopy (cryo-EM). Ribosome interactions anchor the Sec61/TRAP complex in a conformation that renders the ER membrane locally thinner by significantly curving its lumenal leaflet. We propose that TRAP stabilizes the ribosome exit tunnel to assist nascent polypeptide insertion through Sec61 and provides a ratcheting mechanism into the ER lumen mediated by direct polypeptide interactions.

Quantitative Comparison against Experiments Reveals Imperfections in Force Fields' Descriptions of POPC-Cholesterol Interactions.

Cholesterol is a central building block in biomembranes, where it induces orientational order, slows diffusion, renders the membrane stiffer, and drives domain formation. Molecular dynamics (MD) simulations have played a crucial role in resolving these effects at the molecular level; yet, it has recently become evident that different MD force fields predict quantitatively different behavior. Although easily neglected, identifying such limitations is increasingly important as the field rapidly progresses toward simulations of complex membranes mimicking the in vivo conditions: pertinent multicomponent simulations must capture accurately the interactions between their fundamental building blocks, such as phospholipids and cholesterol. Here, we define quantitative quality measures for simulations of binary lipid mixtures in membranes against the C-H bond order parameters and lateral diffusion coefficients from NMR spectroscopy as well as the form factors from X-ray scattering. Based on these measures, we perform a systematic evaluation of the ability of commonly used force fields to describe the structure and dynamics of binary mixtures of palmitoyloleoylphosphatidylcholine (POPC) and cholesterol. None of the tested force fields clearly outperforms the others across the tested properties and conditions. Still, the Slipids parameters provide the best overall performance in our tests, especially when dynamic properties are included in the evaluation. The quality evaluation metrics introduced in this work will, particularly, foster future force field development and refinement for multicomponent membranes using automated approaches.

Efficient Simulations of Solvent Asymmetry Across Lipid Membranes Using Flat-Bottom Restraints.

The routinely employed periodic boundary conditions complicate molecular simulations of physiologically relevant asymmetric lipid membranes together with their distinct solvent environments. Therefore, separating the extracellular fluid from its cytosolic counterpart has often been performed using a costly double-bilayer setup. Here, we demonstrate that the lipid membrane and solvent asymmetry can be efficiently modeled with a single lipid bilayer by applying an inverted flat-bottom potential to ions and other solute molecules, thereby restraining them to only interact with the relevant leaflet. We carefully optimized the parameters of the suggested method so that the results obtained using the flat-bottom and double-bilayer approaches become mutually indistinguishable. Then, we apply the flat-bottom approach to lipid bilayers with various compositions and solvent environments, covering ions and cationic peptides to validate the approach in a realistic use case. We also discuss the possible limitations of the method as well as its computational efficiency and provide a step-by-step guide on how to set up such simulations in a straightforward manner.

Signal peptide mimicry primes Sec61 for client-selective inhibition.

Preventing the biogenesis of disease-relevant proteins is an attractive therapeutic strategy, but attempts to target essential protein biogenesis factors have been hampered by excessive toxicity. Here we describe KZR-8445, a cyclic depsipeptide that targets the Sec61 translocon and selectively disrupts secretory and membrane protein biogenesis in a signal peptide-dependent manner. KZR-8445 potently inhibits the secretion of pro-inflammatory cytokines in primary immune cells and is highly efficacious in a mouse model of rheumatoid arthritis. A cryogenic electron microscopy structure reveals that KZR-8445 occupies the fully opened Se61 lateral gate and blocks access to the lumenal plug domain. KZR-8445 binding stabilizes the lateral gate helices in a manner that traps select signal peptides in the Sec61 channel and prevents their movement into the lipid bilayer. Our results establish a framework for the structure-guided discovery of novel therapeutics that selectively modulate Sec61-mediated protein biogenesis.

Reduced binding of apoE4 to complement factor H promotes amyloid-ß oligomerization and neuroinflammation.

The APOE4 variant of apolipoprotein E (apoE) is the most prevalent genetic risk allele associated with late-onset Alzheimer's disease (AD). ApoE interacts with complement regulator factor H (FH), but the role of this interaction in AD pathogenesis is unknown. Here we elucidate the mechanism by which isoform-specific binding of apoE to FH alters Aß1-42-mediated neurotoxicity and clearance. Flow cytometry and transcriptomic analysis reveal that apoE and FH reduce binding of Aß1-42 to complement receptor 3 (CR3) and subsequent phagocytosis by microglia which alters expression of genes involved in AD. Moreover, FH forms complement-resistant oligomers with apoE/Aß1-42 complexes and the formation of these complexes is isoform specific with apoE2 and apoE3 showing higher affinity to FH than apoE4. These FH/apoE complexes reduce Aß1-42 oligomerization and toxicity, and colocalize with complement activator C1q deposited on Aß plaques in the brain. These findings provide an important mechanistic insight into AD pathogenesis and explain how the strongest genetic risk factor for AD predisposes for neuroinflammation in the early stages of the disease pathology.

Protein Crowding and Cholesterol Increase Cell Membrane Viscosity in a Temperature Dependent Manner.

Shear viscosity of lipid membranes dictates how fast lipids, proteins, and other membrane constituents travel along the membrane and rotate around their principal axis, thus governing the rates of diffusion-limited reactions taking place at membranes. In this framework, the heterogeneity of biomembranes indicates that cells could regulate these rates via varying local viscosities. Unfortunately, experiments to probe membrane viscosity under various conditions are tedious and error prone. Molecular dynamics simulations provide an attractive alternative, especially given that recent theoretical developments enable the elimination of finite-size effects in simulations. Here, we use a variety of different equilibrium methods to extract the shear viscosities of lipid membranes from both coarse-grained and all-atom molecular dynamics simulations. We systematically probe the variables relevant for cellular membranes, namely, membrane protein crowding, cholesterol concentration, and the length and saturation level of lipid acyl chains, as well as temperature. Our results highlight that in their physiologically relevant ranges, protein concentration, cholesterol concentration, and temperature have significantly larger effects on membrane viscosity than lipid acyl chain length and unsaturation level. In particular, the crowding with proteins has a significant effect on the shear viscosity of lipid membranes and thus on the diffusion occurring in the membranes. Our work also provides the largest collection of membrane viscosity values from simulation to date, which can be used by the community to predict the diffusion coefficients or their trends via the Saffman-Delbrück description. Additionally, it is worth emphasizing that diffusion coefficients extracted from simulations exploiting periodic boundary conditions must be corrected for the finite-size effects prior to comparison with experiment, for which the present collection of viscosity values can readily be used. Finally, our thorough comparison to experiments suggests that there is room for improvement in the description of bilayer dynamics provided by the present force fields.

Surfactant Proteins SP-B and SP-C in Pulmonary Surfactant Monolayers: Physical Properties Controlled by Specific Protein-Lipid Interactions.

The lining of the alveoli is covered by pulmonary surfactant, a complex mixture of surface-active lipids and proteins that enables efficient gas exchange between inhaled air and the circulation. Despite decades of advancements in the study of the pulmonary surfactant, the molecular scale behavior of the surfactant and the inherent role of the number of different lipids and proteins in surfactant behavior are not fully understood. The most important proteins in this complex system are the surfactant proteins SP-B and SP-C. Given this, in this work we performed nonequilibrium all-atom molecular dynamics simulations to study the interplay of SP-B and SP-C with multicomponent lipid monolayers mimicking the pulmonary surfactant in composition. The simulations were complemented by z-scan fluorescence correlation spectroscopy and atomic force microscopy measurements. Our state-of-the-art simulation model reproduces experimental pressure-area isotherms and lateral diffusion coefficients. In agreement with previous research, the inclusion of either SP-B and SP-C increases surface pressure, and our simulations provide a molecular scale explanation for this effect: The proteins display preferential lipid interactions with phosphatidylglycerol, they reside predominantly in the lipid acyl chain region, and they partition into the liquid expanded phase or even induce it in an otherwise packed monolayer. The latter effect is also visible in our atomic force microscopy images. The research done contributes to a better understanding of the roles of specific lipids and proteins in surfactant function, thus helping to develop better synthetic products for surfactant replacement therapy used in the treatment of many fatal lung-related injuries and diseases.

Nanoscale membrane curvature sorts lipid phases and alters lipid diffusion.

The precise spatiotemporal control of nanoscale membrane shape and composition is the result of a complex interplay of individual and collective molecular behaviors. Here, we employed single-molecule localization microscopy and computational simulations to observe single-lipid diffusion and sorting in model membranes with varying compositions, phases, temperatures, and curvatures. Supported lipid bilayers were created over 50-nm-radius nanoparticles to mimic the size of naturally occurring membrane buds, such as endocytic pits and the formation of viral envelopes. The curved membranes recruited liquid-disordered lipid phases while altering the diffusion and sorting of tracer lipids. Disorder-preferring fluorescent lipids sorted to and experienced faster diffusion on the nanoscale curvature only when embedded in a membrane capable of sustaining lipid phase separation at low temperatures. The curvature-induced sorting and faster diffusion even occurred when the sample temperature was above the miscibility temperature of the planar membrane, implying that the nanoscale curvature could induce phase separation in otherwise homogeneous membranes. Further confirmation and understanding of these results are provided by continuum and coarse-grained molecular dynamics simulations with explicit and spontaneous curvature-phase coupling, respectively. The curvature-induced membrane compositional heterogeneity and altered dynamics were achieved only with a coupling of the curvature with a lipid phase separation. These cross-validating results demonstrate the complex interplay of lipid phases, molecular diffusion, and nanoscale membrane curvature that are critical for membrane functionality.

Ionic Strength and Solution Composition Dictate the Adsorption of Cell-Penetrating Peptides onto Phosphatidylcholine Membranes.

Adsorption of arginine-rich positively charged peptides onto neutral zwitterionic phosphocholine (PC) bilayers is a key step in the translocation of those potent cell-penetrating peptides into the cell interior. In the past, we have shown both theoretically and experimentally that polyarginines adsorb to the neutral PC-supported lipid bilayers in contrast to polylysines. However, comparing our results with previous studies showed that the results often do not match even at the qualitative level. The adsorption of arginine-rich peptides onto 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) may qualitatively depend on the actual experimental conditions where binding experiments have been performed. In this work, we systematically studied the adsorption of R9 and K9 peptides onto the POPC bilayer, aided by molecular dynamics (MD) simulations and fluorescence cross-correlation spectroscopy (FCCS) experiments. Using MD simulations, we tested a series of increasing peptide concentrations, in parallel with increasing Na+ and Ca2+ salt concentrations, showing that the apparent strength of adsorption of R9 decreases upon the increase of peptide or salt concentration in the system. The key result from the simulations is that the salt concentrations used experimentally can alter the picture of peptide adsorption qualitatively. Using FCCS experiments with fluorescently labeled R9 and K9, we first demonstrated that the binding of R9 to POPC is tighter by almost 2 orders of magnitude compared to that of K9. Finally, upon the addition of an excess of either Na+ or Ca2+ ions with R9, the total fluorescence correlation signal is lost, which implies the unbinding of R9 from the PC bilayer, in agreement with our predictions from MD simulations.

Accurate Simulations of Lipid Monolayers Require a Water Model with Correct Surface Tension.

Lipid monolayers provide our lungs and eyes their functionality and serve as proxy systems in biomembrane research. Therefore, lipid monolayers have been studied intensively including using molecular dynamics simulations, which are able to probe their lateral structure and interactions with, e.g., pharmaceuticals or nanoparticles. However, such simulations have struggled in describing the forces at the air-water interface. Particularly, the surface tension of water and long-range van der Waals interactions have been considered critical, but their importance in monolayer simulations has been evaluated only separately. Here, we combine the recent C36/LJ-PME lipid force field that includes long-range van der Waals forces with water models that reproduce experimental surface tensions to elucidate the importance of these contributions in monolayer simulations. Our results suggest that a water model with correct surface tension is necessary to reproduce experimental surface pressure-area isotherms and monolayer phase behavior. The latter includes the liquid expanded and liquid condensed phases, their coexistence, and the opening of pores at the correct area per lipid upon expansion. Despite these improvements of the C36/LJ-PME with certain water models, the standard cutoff-based CHARMM36 lipid model with the 4-point OPC water model still provides the best agreement with experiments. Our results emphasize the importance of using high-quality water models in applications and parameter development in molecular dynamics simulations of biomolecules.

The Two Faces of the Liquid Ordered Phase.

Coexisting liquid ordered (Lo) and liquid disordered (Ld) lipid phases in synthetic and plasma membrane-derived vesicles are commonly used to model the heterogeneity of biological membranes, including their putative ordered rafts. However, raft-associated proteins exclusively partition to the Ld and not the Lo phase in these model systems. We believe that the difference stems from the different microscopic structures of the lipid rafts at physiological temperature and the Lo phase studied at room temperature. To probe this structural diversity across temperatures, we performed atomistic molecular dynamics simulations, differential scanning calorimetry, and fluorescence spectroscopy on Lo phase membranes. Our results suggest that raft-associated proteins are excluded from the Lo phase at room temperature due to the presence of a stiff, hexagonally packed lipid structure. This structure melts upon heating, which could lead to the preferential solvation of proteins by order-preferring lipids. This structural transition is manifested as a subtle crossover in membrane properties; yet, both temperature regimes still fulfill the definition of the Lo phase. We postulate that in the compositionally complex plasma membrane and in vesicles derived therefrom, both molecular structures can be present depending on the local lipid composition. These structural differences must be taken into account when using synthetic or plasma membrane-derived vesicles as a model for cellular membrane heterogeneity below the physiological temperature.

Inverse Conformational Selection in Lipid-Protein Binding.

Interest in lipid interactions with proteins and other biomolecules is emerging not only in fundamental biochemistry but also in the field of nanobiotechnology where lipids are commonly used, for example, in carriers of mRNA vaccines. The outward-facing components of cellular membranes and lipid nanoparticles, the lipid headgroups, regulate membrane interactions with approaching substances, such as proteins, drugs, RNA, or viruses. Because lipid headgroup conformational ensembles have not been experimentally determined in physiologically relevant conditions, an essential question about their interactions with other biomolecules remains unanswered: Do headgroups exchange between a few rigid structures, or fluctuate freely across a practically continuous spectrum of conformations? Here, we combine solid-state NMR experiments and molecular dynamics simulations from the NMRlipids Project to resolve the conformational ensembles of headgroups of four key lipid types in various biologically relevant conditions. We find that lipid headgroups sample a wide range of overlapping conformations in both neutral and charged cellular membranes, and that differences in the headgroup chemistry manifest only in probability distributions of conformations. Furthermore, the analysis of 894 protein-bound lipid structures from the Protein Data Bank suggests that lipids can bind to proteins in a wide range of conformations, which are not limited by the headgroup chemistry. We propose that lipids can select a suitable headgroup conformation from the wide range available to them to fit the various binding sites in proteins. The proposed inverse conformational selection model will extend also to lipid binding to targets other than proteins, such as drugs, RNA, and viruses.

Nonconverged Constraints Cause Artificial Temperature Gradients in Lipid Bilayer Simulations.

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

Anisotropic diffusion of membrane proteins at experimental timescales.

Single-particle tracking (SPT) experiments of lipids and membrane proteins provide a wealth of information about the properties of biomembranes. Careful analysis of SPT trajectories can reveal deviations from ideal Brownian behavior. Among others, this includes confinement effects and anomalous diffusion, which are manifestations of both the nanoscale structure of the underlying membrane and the structure of the diffuser. With the rapid increase in temporal and spatial resolution of experimental methods, a new aspect of the motion of the particle, namely, anisotropic diffusion, might become relevant. This aspect that so far received only little attention is the anisotropy of the diffusive motion and may soon provide an additional proxy to the structure and topology of biomembranes. Unfortunately, the theoretical framework for detecting and interpreting anisotropy effects is currently scattered and incomplete. Here, we provide a computational method to evaluate the degree of anisotropy directly from molecular dynamics simulations and also point out a way to compare the obtained results with those available from SPT experiments. In order to probe the effects of anisotropic diffusion, we performed coarse-grained molecular dynamics simulations of peripheral and integral membrane proteins in flat and curved bilayers. In agreement with the theoretical basis, our computational results indicate that anisotropy can persist up to the rotational relaxation time [τ=(2Dr)-1], after which isotropic diffusion is observed. Moreover, the underlying topology of the membrane bilayer can couple with the geometry of the particle, thus extending the spatiotemporal domain over which this type of motion can be detected.

Martini 3: a general purpose force field for coarse-grained molecular dynamics.

The coarse-grained Martini force field is widely used in biomolecular simulations. Here we present the refined model, Martini 3 ( http://cgmartini.nl ), with an improved interaction balance, new bead types and expanded ability to include specific interactions representing, for example, hydrogen bonding and electronic polarizability. The updated model allows more accurate predictions of molecular packing and interactions in general, which is exemplified with a vast and diverse set of applications, ranging from oil/water partitioning and miscibility data to complex molecular systems, involving protein-protein and protein-lipid interactions and material science applications as ionic liquids and aedamers.

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.