Brominated lipid probes expose structural asymmetries in constricted membranes

Lipids in biological membranes are thought to be functionally organized, but few experimental tools can probe nanoscale membrane structure. Using brominated lipids as contrast probes for cryo-EM and a model ESCRT-III membrane-remodeling system composed of human CHMP1B and IST1, we observed leaflet-level and protein-localized structural lipid patterns within highly constricted and thinned membrane nanotubes. These nanotubes differed markedly from protein-free, flat bilayers in leaflet thickness, lipid diffusion rates and lipid compositional and conformational asymmetries. Simulations and cryo-EM imaging of brominated stearoyl-docosahexanenoyl-phosphocholine showed how a pair of phenylalanine residues scored the outer leaflet with a helical hydrophobic defect where polyunsaturated docosahexaenoyl tails accumulated at the bilayer surface. Combining cryo-EM of halogenated lipids with molecular dynamics thus enables new characterizations of the composition and structure of membranes on molecular length scales.

Lipids in biological membranes are thought to be functionally organized, but few experimental tools can probe nanoscale membrane structure. Using brominated lipids as contrast probes for cryo-EM and a model ESCRT-III membrane-remodeling system composed of human CHMP1B and IST1, we observed leaflet-level and protein-localized structural lipid patterns within highly constricted and thinned membrane nanotubes. These nanotubes differed markedly from protein-free, flat bilayers in leaflet thickness, lipid diffusion rates and lipid compositional and conformational asymmetries. Simulations and cryo-EM imaging of brominated stearoyl-docosahexanenoyl-phosphocholine showed how a pair of phenylalanine residues scored the outer leaflet with a helical hydrophobic defect where polyunsaturated docosahexaenoyl tails accumulated at the bilayer surface. Combining cryo-EM of halogenated lipids with molecular dynamics thus enables new characterizations of the composition and structure of membranes on molecular length scales.
Cells use molecular machines to form and remodel their membrane-defined compartments' compositions, shapes and connections. The regulated activity of these membrane-remodeling machines drives processes such as vesicular traffic and organelle homeostasis. However, the precise mechanisms by which proteins generate mechanical force to catalyze membrane fission, fusion and shape changes remain elusive. Further, the structural evolution of lipids and lipid bilayers during these processes are challenging to study, as is quantifying the energetic contributions of individual lipid species to membrane remodeling. A better understanding of membrane and lipid structure during remodeling is crucial to understanding both the mechanisms of membrane remodeling and, more broadly, the interactions between proteins and membranes.
Molecular-scale insights into the lipid-leaflet, lipid-lipid and lipid-protein dynamics that generate extreme membrane curvature could clarify the mechanisms of membrane remodeling. A thorough understanding of membrane mechanics will account for lipid asymmetry and flip/flop 1,2 , spontaneous lipid curvature 1,3-6 , bending rigidity 2,5,7 , line tensions 8 and protein-generated forces including amino acid insertions 5,9,10 , lateral pressure from protein crowding 11 , shearing forces 7 and lipid-protein 'friction' 12 . The relative contributions of each of these processes to membrane remodeling generally, and to membrane constriction in the present case, are challenging to measure experimentally or explain in structural terms. We need new measures of how membrane structure, composition and protein-generated forces influence membrane properties in vitro and in vivo.
So far, cryo-EM reconstructions of membrane-bound proteins have generated insights into the machines that shape membranes, including dynamin family proteins [13][14][15][16][17] , OPA1 (ref. 18), LPOR 19 , SNAREs 20 , BAR domain containing proteins 21 and ESCRTs 22 . However, these studies have generally failed to reveal molecular-level information about the membrane itself because most lipids display nearly indistinguishable electron scattering 23,24 . Fluid bilayers, moreover, are generally thought to lack a structured pattern at the nanoscale that is recoverable by Article https://doi.org/10.1038/s41594-022-00898-1 40 Å-nearly to the point of fission 22,40 . This system provides an opportunity for probing how peripheral membrane proteins engage distinct lipid species to pattern membrane properties more generally. Additionally, this highly constricted state may be a metastable or intermediate state on pathway to membrane fission and therefore may illuminate the mechanical and molecular properties underlying membrane fission reactions. Previously, we observed that helix α1 of CHMP1B appeared to 'dimple' the outer surface of highly constricted membrane nanotubes, but the properties of this deformation were poorly defined (Fig. 1a,b,e) 22 . To explore the membrane's structure in more detail, we performed coarse-grained-molecular dynamics (CG-MD) simulations of the membrane tubule in the presence of the CHMP1B-IST1 copolymer using the Martini 2.2 CG force field (Fig. 1c,d and Extended Data Fig. 1). We used the previously optimized mixture of stearoyl-docosahexanenoyl-phosphatidylcholine (SDPC), palmitoyl-oleoyl-phosphatidylserine (POPS), di-oleoylphosphatidylinositol-3,5-bisphosphate] (PIP 2 ) and cholesterol (CHOL) for the simulation and subsequent cryo-EM experiments (see Supplementary Fig. 1 for chemical structures) 22 . Simulations were started with a spontaneous tubule assembly MD phase with randomized initial lipid positions between replicates, followed by an extended leaflet equilibration during which lipids flip-flop at artificially induced membrane pores to arrive at consistent inner and outer leaflet compositions across replicates before production data collection (see Methods and image averaging. Therefore, it is typically impossible to distinguish individual lipid species, even approximately, except in high-resolution structures where isolated, static lipids are resolved bound to transmembrane proteins [25][26][27][28] . To overcome this shortcoming, here we describe halogen-labeled lipid probes as contrast-enhancing agents for cryo-EM. We use a model system consisting of human ESCRT-III proteins CHMP1B and IST1 to remodel vesicles containing lipid probes into lipid nanotubes with extremely high curvature 22,29,30 . Using labeled lipids, ESCRT-III proteins, cryo-EM and molecular dynamics simulations, we characterize how ensembles of lipid-protein interactions, lipid conformational changes and the resulting stabilization of strong membrane asymmetries drive membrane constriction and thinning. Leaflet and lipid shape change together with specific interactions between lipids and the ESCRT-III proteins, stabilize membrane curvature in this snapshot of a membrane on the verge of scission.

ESCRT-III protein CHMP1B α1 induces a furrow in the bilayer
When exposed to model membranes in vitro, ESCRT-III proteins can shape membranes into structures with zero 31,32 , positive 22,29,[33][34][35][36] or negative curvature 31,32,[37][38][39] . We recently demonstrated how a pair of human ESCRT-III proteins, CHMP1B and IST1, assemble into helical filaments that act in sequence to remodel membranes into high-curvature nanotubes with an inner diameter of only roughly Supplementary Information for full details). The simulations revealed CHMP1B side chains sterically displacing phospholipid headgroups and exposing the hydrophobic core of the membrane in an apparent 'scoring' of the surface-rather than an elastic dimpling of a continuous headgroup surface as we previously thought 22 ( Fig. 1f and Supplementary Fig. 2). We note that unlike the cryo-EM reconstructions, there is no symmetry applied to the bilayer during the simulations and the tubule fluctuates. Deviations from an ideal cylindrical shape are most obvious at the upper and lower boundary, where the imperfect periodicity of the protein coat results in an exposed membrane. The morphology of the bilayer tube, especially in this region, is sensitive to the number of initial lipids (Extended Data Fig. 1).
To aid visualization of the headgroup exclusion zone, we transformed the protein and lipid positions into cylindrical coordinates and plotted the mean density of headgroups in the outer leaflet (Fig. 1h). The headgroup-excluding furrow aligned with the bulky hydrophobic residues F9 and F13, which protrude from helix α1 on CHMP1B (green dots in Fig. 1h). The bilayer appeared unperturbed away from F9 and F13. Notably, these two residues are the only hydrophobic residues that face the bilayer, with the rest of the bilayer-interacting residues containing cationic side chains. Additionally, given the orientation of CHMP1B helix α1 relative to the curved tubule surface, the pair of phenylalanines are the only residues that protrude into the bilayer (Supplementary Fig. 3a). Despite the 23 Å distance between F9 and F13 pairs on adjacent CHMP1B subunits, the furrow forms a continuous stripe of excluded headgroups in the membrane outer leaflet in simulation ( Fig. 1f,h). This aspect of the simulated structures is in agreement with the experimental structure determined by cryo-EM (Fig. 1e,g), with the EM-derived density map revealing a furrow of low Coulombic potential aligned with the stripe of F9 and F13 positions along the CHMP1B helical filament. Together, the EM maps and CG-MD results indicate that the observed deformation constitutes a stable lipid packing defect due to headgroup displacement by these hydrophobic amino acids of CHMP1B and concomitant exposure of the lipid tails underneath (Fig. 1g,h). Line tension associated with the hydrophobic defect may energetically drive the defects to coalesce into a continuous, helical stripe. Because such an asymmetry in the membrane should alter its mechanical properties, influencing the energetic barriers associated with curvature generation or fission, we sought to further characterize its structure.

Specific lipids enrich at contact sites and change shape
We next examined the molecular structure and dynamics of the lipids throughout the simulated lipid bilayer. Specifically, we analyzed the radii of each of the lipid beads and the angles of each of the lipid tails with respect to the bilayer normal. This analysis revealed pronounced asymmetries in lipid shape correlated with leaflet localization and proximity to the CHMP1B F9, F13 contact site (Fig. 2a-d). The lipid SDPC, which has a polyunsaturated (22:6) tail at the sn-2 position and a fully saturated (18:0) tail at sn-1, showed the strongest asymmetries. While the saturated sn-1 tail typically remained aligned with the membrane normal throughout both leaflets (Fig. 2a-d), the polyunsaturated sn-2 tail of outer leaflet lipids bent backward and radially outward toward  the headgroup region at the F9, F13 contact site (Fig. 2a,c). This shape transformation partially filled the space vacated by the displaced headgroups and enabled the tail to interact directly with the hydrophobic side chains of F9 and F13. Away from the exclusion zone, the polyunsaturated tail adopted a more typical conformations, largely oriented along the membrane normal (Fig. 2b,d).
To quantify lipid conformation, we calculated the lipid tilt angles and radial densities of the phosphate and tail beads of all phospholipids adjacent to and in between the Phe contact sites for both leaflets (Fig. 2e-h and Extended Data Fig. 2). We defined the 'Phe contact site' as the 30° wedge centered through the angular position of F9 and F13, while the 'gaps between the contact sites' were defined as the 30° wedge centered 180° across from the contact sites ( Supplementary Fig. 3). Consistent with our anecdotal observations of tail 'backflipping', the sn-2 tail of SDPC was highly deformed in the outer leaflet at the Phe contact sites with its mean location shifted radially outward and a substantial density beyond both the bilayer phosphate peak, typically the membrane's outer surface, and the Phe side chains (Fig. 2e, r = 60 Å). On the other hand, in the inner leaflet, both tails are largely unperturbed throughout (Fig. 2e,f). The other lipid tails showed more minor degrees of deformation at the Phe contact sites and normal bilayer behavior in the gaps between furrows.
The tilt angle analysis revealed corresponding asymmetries (Fig. 2g,h). In both leaflets, SDPC tails adopted typical values of roughly 30° with respect to the bilayer normal-except for the polyunsaturated sn-2 tail at the furrows, which displayed pronounced tilting with a mean value of 60° and some extreme conformations greater than 90° where the tail bead resided at the membrane surface beyond the radial position of F9 and F13. The other two phospholipid species, POPS and PIP 2 , which have both saturated and monounsaturated lipid tails, did not exhibit a similarly strong location-dependent perturbation (Extended Data Fig. 2a,b).
These observations are consistent with the notion that polyunsaturated lipid tails enable membrane shape plasticity under curvature stress, which previous work has suggested may be a general principle of bilayer deformation 3,22,41-45 . Indeed, the polyunsaturated tail of SDPC is known to be highly flexible. This property could facilitate membrane bending if the lipid changes its mean shape in response to curvature stress [46][47][48][49] . Consistent with this hypothesis, replacing SDPC with monounsaturated palmitoyl-oleoyl-phosphatidylcholine (POPC) precluded high-curvature tubule formation but did not preclude membrane binding (Extended Data Fig. 3). Specifically, CHMP1B and IST1 were unable to constrict POPC-containing liposomes (lipid composition 3) to high curvature, as assessed by negative stain transmission electron microscopy and cryo-EM (Extended Data Fig. 3).
We next sought direct experimental validation of SDPC's localization in the bilayer and whether it undergoes such dramatic shape transformations. To detect individual lipid species using cryo-EM, we synthesized analogs of each lipid with bromine atoms along the unsaturated lipid tail lengths. Due to their more massive nuclei, bromine atoms generate contrast through enhanced electron scattering. The brominated and unbrominated lipids displayed similar bilayer phase behavior and Langmuir monolayer properties, indicating comparable lipid packing and fluidity (Extended Data Fig. 4h,i). Additionally, in bilayers, SDPC-Br and POPS-Br are fluid at room temperature and macroscopically homogeneous, consistent with previous studies 50,51 . CHMP1B and IST1 remodeled brominated vesicles into narrow nanotubes yielding protein structures indistinguishable at roughly 3 Å resolution from protein-bound nanotubes prepared with native lipids (Fig. 3, Extended Data Fig. 3 and Supplementary Fig. 4). Together, these results indicate that halogenation does not meaningfully perturb the properties of these lipids or the bilayers formed from them. We collected high-resolution cryo-EM data and achieved high-resolution reconstructions with vesicles in which each lipid was replaced by its halogenated analog (Table 1, lipid compositions [10][11][12][13][14] or different fractions of SDPC were brominated (Table 1, lipid compositions 10,15,16).
To quantitively compare the reconstructions with brominated and unbrominated lipids, the pixel value distributions were normalized to the CHMP1B intensity from radial averages. We first compared horizontal and vertical slices through each filament and observed differences in the intensities of the lipid bilayer leaflets with brominated lipids. In particular, the SDPC-Br reconstruction showed intense spots of enhanced density adjacent to CHMP1B F9 and F13, consistent with SDPC-Br enrichment at the site (Fig. 3c) and corroborating the CG-MD simulation data showing that lipid tails reached the surface of the bilayer. The other brominated lipids, which are monounsaturated, also accumulated at the membrane-protein contact site, but to a smaller extent. Further, we resolved differences in the positions of different lipid species relative to the Phe contact sites, with the brominated SDPC tail in closer proximity to F9 and F13 both radially and axially than brominated PIP 2 and POPS tails (Fig. 3e,f). These last two anionic lipids more closely approached K16 and R20, two cationic residues previously identified as critical for CHMP1B membrane binding to membranes 22 .
The altered SDPC tail conformation where the bilayer contacts CHMP1B residues F9 and F13 suggests at least two hypotheses. First, lipid headgroup displacement by F9 and F13 creates lipid packing defects and curvature stress that is relieved when lipid tails fill in the packing defects. Second, hydrophobic interactions drive the lipid tails to interact with the nonpolar Phe side chain. These ideas, the first focused on steric packing and the second focused on chemical interactions between lipids and amino acids, are not mutually exclusive. We generated a series of CHMP1B double mutants, F9X + F13X (where X = A, E or L), to test the relative contributions of each possibility.
All of these CHMP1B double mutants still bound to and remodeled native and brominated lipid bilayers, forming membrane-bound copolymers with IST1 that were indistinguishable from the wild-type (WT) CHMP1B/IST1 copolymer in helical symmetry parameters (Fig. 4). F9 and F13 substitution by small nonpolar A side chains eliminated SDPC-Br enrichment at CHMP1B contact sites (in lipid composition 10). SDPC also failed to concentrate in the F to E double mutant. In both the F to A and the F to E double mutants, the protein appeared to induce an elastic deformation of the bilayer without opening an appreciable hydrophobic defect. The observed deformation in the outer membrane was measurably deeper into the bilayer than for WT CHMP1B and did not induce the helical stripe of low Coulombic potential observed in the WT protein. Simulations of the F to A double mutant also showed reduced SDPC tail enrichment and backflipping, as well as increased headgroup density shifted radially inward, at CHMP1B contact sites, consistent with an elastic deformation (Extended Data Fig. 5). Additionally, SDPC sn-2 tails were less perturbed than with WT protein in simulation. By contrast, in the cryo-EM reconstruction, SDPC-Br was observed to accumulate at contact sites with large hydrophobic Leu side chains to a similar but lesser degree than the WT protein, indicating the formation of a hydrophobic defect. Finally, removal of the protein coat in simulations resulted in dramatic reduction of SDPC tail backflipping and recovered mostly normal bilayer structure and lipid conformation (Extended Data Fig. 6). Together, these observations are consistent with chemical interactions between greasy side chains and polyunsaturated tails stabilizing the backflipped lipid tail conformations. On the other hand, the lack of SDPC sn-2 tail enrichment at the highly curved elastic deformations in the F to A and F to E mutants suggests that curvature stress at the contact site does not play a notable role in SDPC accumulation at WT CHMP1B-bilayer contact sites. These results all indicate that hydrophobic interactions, rather than spontaneous curvature, drive local accumulation of polyunsaturated lipid tails at CHMP1B F9 and F13.

Constriction leads to anisotropic inner leaflet thinning
Previous work has predicted that a highly curved bilayer's inner, concave leaflet would be thinner than the outer, convex leaflet based on atomistic simulations 52 . We experimentally measured leaflet thicknesses to test this hypothesis by examining the Coulombic potential (electron scattering) profile from the unbrominated sample and fitting the profile with three Gaussian functions using least-squares 53 .
As shown in Fig. 5a,c, the inner concave leaflet is 3.6 Å thinner than the   outer convex leaflet, in agreement with Yesylevskyy et al. 52 and our own CG-MD simulations of the same system (Fig. 5b,c).
We believe that differences in per-leaflet lipid tail angles seen in simulation may explain the different leaflet thicknesses. In the WT CHMP1B-IST1 with lipid composition 1 simulations, inner leaflet lipid tails sample higher tilt angles than those in the outer leaflet except at the F9, F13 contact site (Fig. 2g,h and Extended Data Fig. 2b), where there is extreme outer leaflet backflipping. Throughout the rest of the bilayer, outer leaflet lipids behave like those in the gap zone (Fig. 2h) and sample lower tilt angles. In contrast, to fill the greater volume nearer the bilayer center versus at the inner leaflet headgroups, inner leaflet lipid tails sample higher tilt angles, effectively shortening those lipids in the radial coordinate and leading to the observed inner leaflet thinning. This effect is also present in simulations with the F9A, F13A mutant (Extended Data Fig. 5d,e) and to a lesser extent in simulation of a protein-free tubule (Extended Data Fig. 6e,f), although that tubule equilibrates to markedly different leaflet compositions than the protein-bound tubules (Extended Data Fig. 6a,c). The increased disorder in the inner leaflet that accompanies higher tilt angles and thinning may facilitate tubule scission by enabling formation of a hemifusion stalk across the tubule lumen.

The constricted bilayer is compositionally heterogenous
Both MD simulations and cryo-EM reconstructions suggested the bilayer composition is not homogeneous, with differences at the CHMP1B F9, F13 contact site versus the gap between the contact sites. By integrating the areas under the radial curve for each leaflet from bromine-labeled reconstructions normalized to the unlabeled reconstruction, we estimated the composition of each leaflet in these two regions (Fig. 6a,b and Extended Data Fig. 4a-d). As shown in Fig. 6c, glycerol inner and outer leaflet planes (pink) for calculating leaflet thicknesses. c, Summary of calculated leaflet and bilayer thicknesses from cryo-EM and CG-MD. Uncertainty is reported as the standard deviation between independent measurements from each half map for cryo-EM, and as the standard deviation between ten simulation replicates for CG-MD. Data for graphs in a and b are available as source data.    Fig. 3a). We hypothesize that such large lipid enrichments can occur because most filaments contain an unremodeled vesicle protruding from the end that can act as a lipid reservoir (for example, Extended Data Fig. 3a). There may also be compositional variation in the initial vesicle population and consequently bias in which vesicles are remodeled and therefore contribute to our cryo-EM reconstructions. We repeated the CG-MD simulations with the composition estimated from cryo-EM (Table 1, lipid composition 17), and this updated simulation displayed analogous asymmetric leaflet thinning and outer leaflet headgroup exclusion displayed by the CG-MD simulations with the initial lipid composition (Extended Data Fig. 7, Supplementary Information Table 8).
Both simulation and experiment indicated that cholesterol was enriched asymmetrically in the inner leaflet versus the outer leaflet at bilayer-protein contact sites. We confirmed the behavior of cholesterol by solving cryo-EM structures with different concentrations of unbrominated cholesterol in the initial lipid vesicles ( Table 1, lipid compositions 4-9). We took advantage of the fact that cholesterol has a lower Coulombic potential than phospholipids. As the concentration of cholesterol in the lipid mixtures increased, the inner leaflet decreased in intensity relative to the outer leaflet. This observation is also consistent with the cholesterol enrichment in the inner leaflet seen in simulation (Extended Data Fig. 8d,e).
We also estimated that PIP 2 was highly enriched not only in the outer leaflet compared with the inner leaflet, but also in the membrane nanotube compared with the initial bulk composition of the vesicle. In light of this, we next investigated the importance of PIP 2 for tubule formation by carrying out membrane-remodeling assays with vesicles containing no PIP 2 ( Table 1, lipid composition 2). In this case, CHMP1B is unable to stably bind to and remodel vesicles, suggesting that electrostatic interactions are indispensable for its functions (Extended Data Fig. 3i,m). Similarly, K16E or R20E mutations on CHMP1B helix α1 ablate membrane binding in the presence of PIP 2 , presumably due to a lack of favorable electrostatic interactions between cationic residues on the protein and anion-rich lipid bilayers ( Supplementary Fig. 5).
Finally, due to the pinning of PIP 2 by K and R residues and of SDPC by F9 and F13 residues, the diffusion coefficients of lipids in the outer leaflet of the bilayer are dramatically slowed in CG-MD simulations (Extended Data Fig. 9), consistent with recent experimental results that used fluorescence recovery after photobleaching to measure lipid lateral diffusion in CHMP1B/IST1 filaments formed using optical tweezers 30 . This reduction in diffusion coefficients has functional implications for membrane fission as membrane shearing becomes possible as the diffusion rate of lipids decreases 7 , and recent work has suggested that CHMP1B/IST1 catalyze membrane fission through a shearing mechanism 30 . The necessity of SDPC and PIP 2 for membrane binding and constriction highlights the importance of local lipid composition as a regulator of membrane remodeling.

Discussion
We combined cryo-EM, heavy-atom-labeled lipid probes and CG-MD simulations to investigate membrane structure and mechanics in ESCRT-III constricted lipid tubules. These efforts revealed how phenylalanine residues from CHMP1B helix 1 displace headgroups to 'score' the surface of the bilayer and create a hydrophobic defect. Flexible, polyunsaturated lipid tails enrich within the furrow created by these Phe residues, where they contort to interact with the hydrophobic side chains. These surface defects coalesce into a helical furrow that may, along with lipid headgroup displacement by F9 and F13, stabilize the bilayer at this high degree of curvature by decreasing lipid density in the convex outer leaflet. We also showed that this high membrane curvature system generates compositionally and dynamically asymmetric leaflets. While lipid tail backflipping is possible for any of the lipid tails, polyunsaturated tails are exceptionally flexible 3,45,48 and therefore have low energy barriers to 'backflipping' 45,49,54 . This role for polyunsaturated lipids in stabilizing curvature and filling in packing defects by undergoing curvature-driven shape changes 43,49 is consistent with previous reports that SDPC facilitates membrane fission 3 and that polyunsaturated tails prefer to reside near the lipid-water interface, especially in curved bilayers 3,47,49 . Additionally, polyunsaturated lipids preferentially interact with some membrane proteins via van der Waals interactions 46,55 .
Together, our findings on compositional variation within the leaflets of ESCRT-III-formed lipid tubules, and the membrane-remodeling assays with different initial lipid compositions, are consistent with lipid compositions having a prominent effect on the energy required to bend a membrane. Membrane-remodeling machines have adapted accordingly to possess organelle-specific properties, and cells likely govern membrane remodeling by regulating local lipid composition through clustering, transport and lipid metabolism 56 .
Structural studies of higher order assemblies of biological macromolecules, which are not defined by the linear sequences of their components, are challenging due to the heterogeneity, disorder and dynamics of these assemblies. Membranes and biomolecular condensates are examples of such systems that are difficult to study with molecular resolution 57-60 . We anticipate that the ability to validate or falsify simulations with heavy-atom cryo-EM will enable new questions to be answered about these critical fluid assemblies, their interactions and the dynamics governing their functions.

Online content
Any methods, additional references, Nature Portfolio reporting summaries, source data, extended data, supplementary information, acknowledgements, peer review information; details of author contributions and competing interests; and statements of data and code availability are available at https://doi.org/10.1038/s41594-022-00898-1.

Synthesis and characterization of bromolipids
Bromolipids were synthesized as previously described 61 . Here, 1-100 mg of each lipid was dissolved in CHCl 3 to 1-10 mg ml −1 . The lipid solutions were stirred on ice, and bromine (stoichiometric with the number of double bonds in the lipid) was added dropwise. The solution was stirred on ice in the dark for 1 h. Solvent and excess bromine were removed by application of vacuum overnight in the dark. Brominated lipids were then aliquoted and stored at −80 °C until needed. To maintain accurate lipid stock concentrations, no further purification was performed. Consequently, the brominated lipids also contain some unbrominated lipid, whose quantities were determined by mass spectrometry and nuclear magnetic resonance (NMR) spectroscopy. Proton NMR ( 1 H NMR) and carbon NMR ( 13 C NMR) spectra were recorded on a Bruker Avance III HD 400 spectrometer in the UCSF NMR Core Facility using Bruker TopSpin 4. NMR spectra were recorded in CDCl 3 , except for PIP 2 -Br that was recorded in 20:10:1 CDCl 3 :MeOD:D 2 O. NMR spectra were analyzed with MestReNova v.14.2.1. Electrospray ionization mass spectra were recorded with a Waters SQD2 mass spectrometer in the Stanford University Mass Spectrometry Facility. See Supplementary Fig. 1 for NMR and mass spectra. The average numbers of bromine atoms per lipid molecule, as measured by NMR and mass spectrometry, are shown in Supplementary Table 1. We measured pressure-area isotherms for pure unbrominated and brominated lipids as described in detail previously 62 . Briefly, we used a KSV NIMA KN 2002 (Biolin Scientific) Langmuir-Blodgett system with a 273 cm 2 Teflon trough and symmetric Delrin barriers controlled by KSV Nima Attension v.2.3. A piece of Whatman no. 1 filter paper was used as a Wilhelmy plate to monitor the surface pressure in the trough. We added water to the clean trough and used a vacuum line attached to a Teflon tip to remove any contaminants from the liquid surface. Lipids (1 mM in chloroform) were spread onto the surface of the water with a glass microsyringe and solvent was allowed to evaporate for 10 min. The barriers were then compressed at a rate of 10 mm min −1 and pressure-area isotherms were collected. All isotherms were recorded at room temperature on a subphase of Milli-Q 18 MΩ water.

Sample preparation
Membrane-remodeling assays were carried out as previously described to create CHMP1B-IST1 copolymer-bound membrane nanotubes 22 . Small unilamellar vesicles (SUVs) were formed by extrusion. To protect the polyunsaturated lipids from oxidation during the formation of vesicles and subsequent membrane-remodeling reaction, we added 300 nmol of the desired lipid mixture to a glass scintillation vial under a nitrogen atmosphere (see Table 1 and Supplementary Table 2 for lipid mixtures). The solvent was evaporated under a stream of nitrogen while rotating the vial. The resulting film was redissolved in 100 µl chloroform and evaporated again. The vial with the dried lipid film was placed under house vacuum in the dark for 2 h. Next, 250 µl of GF buffer (25 mM Tris pH 8.0 and 125 mM NaCl) was added to the scintillation vial and incubated for 10 min in the dark. The lipid film was next resuspended to form multilamellar vesicles by vortexing. SUVs were formed by extrusion 31 times through a polycarbonate membrane with 50 nm pore size using the Avanti Polar Lipids Mini Extruder. SUVs were either used immediately or aliquoted, snap frozen in liquid nitrogen and stored at −80 °C. To form copolymer-bound nanotubes, SUVs (0.5 mg ml −1 in GF buffer) were incubated with 5 mg ml −1 CHMP1B (in GF buffer) for approximately 6 h at room temperature in the dark. IST1 (in GF buffer + 5% glycerol) was added to 10 mg ml −1 and the solution was incubated overnight at room temperature in the dark.
Samples for negative stain EM and cryo-EM were prepared similarly to previous reports 22 . Briefly, for negative stain EM, 4 µl of the membrane-remodeling mixture were applied to glow-discharged, 200 Cu mesh carbon-coated grids (Electron Microscopy Supplies). Grids were stained with 0.75% (w/v) uranyl formate (Structure Probe, Inc.).
All samples were imaged on an FEI Tecnai T12 120 kV electron microscope equipped with a Gatan UltraScan 895 4k CCD camera. For cryo-EM, 4 µl of membrane-remodeling mixture were pipetted onto glow-discharged R1.2/1.3 Quantifoil 200 Cu mesh grids (Quantifoil) in a Mark IV Vitrobot (FEI). After a 10 s wait time at 19 °C and 100% humidity, grids were blotted for 4 s with Whatman no. 1 filter paper and plunged into liquid ethane. Grids were stored under liquid nitrogen until imaged.

Cryo-EM imaging
Cryo-EM data were collected on 300 kV FEI Titan Krios or 200 kV FEI Talos Arctica microscopes. The Krios was equipped with a Gatan Bio-Quantum energy filter and Gatan K3 direct electron detector. The Krios was operated at a nominal magnification of ×105,000 with a total dose of 67 e − /A 2 . The Arctica microscope was operated at a nominal magnification of ×28,000 or ×36,000 with a total dose of 61 e − /A 2 . The Arctica was equipped with a Gatan K3 direct electron detector. The cameras on both microscopes were operated in correlated double sampling and super resolution modes. Samples with lipid compositions 2-8 and 18 were imaged with the Arctica, while samples with lipid compositions 1 and 10-17 were imaged with the Krios. Images were collected with a 50 µm C2 aperture and 100 µm objective aperture. The defocus was varied between −0.6 and −2 µm. Semiautomated data collection was carried out with SerialEM.

Data processing
The software programs used to process and visualize data were compiled and configured by SBGrid 63 .
We processed the cryo-EM data in a manner similar to that reported previously 22 , except that filaments were automatically picked from the micrographs with SPHIRE-crYOLO 64 . The rest of the data processing was performed in RELION v.3.0.8 (ref. 65). All of the cryo-EM datasets were processed with the pipeline shown in Extended Data Fig. 10 unless otherwise noted. Sample motion was corrected with MotionCor2 with dose weighting, and defocus values were estimated with CTFfind4 (refs. 66,67). One round of 2D classification was performed to create a more homogenous data set with similar particle diameters and eliminate particles with poor resolution. A hollow, smooth cylinder was used as an initial model for three-dimensional (3D) auto-refinement without a mask, and 'ignore CTFs until first peak' option for contrast transfer function (CTF) estimation was used. Previously determined helical parameters were used as the initial values in the refinement. Subsequently, these particles underwent a single round of 3D classification to separate 17 and 18 asymmetrical subunits (asu)-per turn structures. Specifically, the number of asymmetrical units, initial twist and initial rise parameters were changed to select for either variant of the filament. During 3D classification, we used a mask that included both the protein and lipid bilayer, whereas further auto-refinements used masks that solely included the protein. High-resolution 3D classes with the desired symmetry were chosen, and each group of classes with the same symmetry went through another round of 3D auto-refinement. Afterward, the refined particles for both the 17 and 18 asu structures went through two rounds of per-particle CTF refinement, with 3D auto-refinement between each round, and the final reconstructions were postprocessed with automatic B factor estimation. The relevant parameters for each step are listed in Supplementary Tables 3-5. We used the CHMP1B-IST1 atomic model previously reported, PDB 6TZ5, which fit all of the 17 subunit per turn reconstructions reported here.
We estimated the local resolution of the lipid bilayer to be roughly 6 Å by gold standard Fourier shell correlation analysis. We created a mask that only encompassed the bilayer and performed an auto-refinement with Fourier shell correlation resolution estimate in RELION. Consequently, we used the same mask to apply a custom 6 Å low-pass filter to the bilayer in the reconstructions while maintaining the high-resolution, roughly 3 Å, features of the protein coat. All the cryo-EM densities shown in the main text have been processed in this https://doi.org/10.1038/s41594-022-00898-1 way. We also confirmed that the application of helical symmetry, normalization and low-pass filter did not alter the results of the analyses (Supplementary Fig. 6).

Analysis of leaflet compositions
We analyzed the leaflet compositions by comparing the increase in intensity of each leaflet for brominated samples versus the unbrominated reference according to the procedure shown schematically in Extended Data Fig. 4a-d. First, we investigated the structure and composition of the membrane-protein contact sites and the region of the bilayer farthest from the protein (gaps between contact sites). To facilitate this comparison, we rotated each z slice of the cryo-EM reconstructions so that membrane-protein contact sites were aligned axially. We then projected the central 30% of each reconstruction along the z axis, improving the signal to noise ratio. We then calculated radial averages and the area of each leaflet curve above 0, which we defined as the solvent background intensity. The areas for each leaflet of the samples without brominated lipids were subtracted from the areas of the respective leaflets of each of the reconstructions with brominated lipids. This yielded the extra scattering in each leaflet due to each brominated lipid. These values were then normalized to the numbers of bromine atoms per molecule (Supplementary Table 1). The values for each leaflet were summed, and the fraction of the total for each brominated lipid was calculated, such that the composition of each leaflet adds to 100 mol%. Second, the bulk, radially averaged compositions of each leaflet were estimated by applying the same procedure but omitting the rotation alignment of the membrane-protein contact sites. Both analyses were performed separately on each of the half maps, and the reported uncertainties are the standard deviations between the analyses from the two half maps. The estimated bulk leaflet compositions are shown in Supplementary Table 6. We also used a commercially available 19-iodo analog of cholesterol to show that the partitioning of this probe between leaflets is similar to that of brominated cholesterol ( Supplementary Fig. 7).

Summary of simulations
Five sets of simulations were conducted. The first set of membrane tubule simulations produced most of the data discussed in the main text (all main text simulation data except Fig. 6). This set consisted of ten replicate simulations using a CG Martini 2.2 model of two turns of the IST1-CHMP1B protein coat (PDB 6TZ5) using a lipid molar composition corresponding to starting lipid composition 1 (Table 1). Full results are presented in Extended Data Fig. 2. The second set was identical to the first, except in silico mutagenesis was used to make the F9A + F13A double mutant of CHMP1B and only three replicates were run (full results in Extended Data Fig. 5). The third set were run using the experimentally derived membrane tubule composition, again with WT CHMP1B (Table 1, composition 17). This more highly anionic lipid mixture was substantially less stable than composition 1, and hence the equilibration procedures and run parameters were adjusted as described in the Supplementary Information (full results in Extended Data Fig. 7). The fourth simulation was run to compare features of the tubules formed in the presence of protein, from simulation set 1, with a tubule of equal curvature but formed without protein. A tubule was formed from spontaneous assembly using a lipid mixture corresponding to composition 1 and run through equilibration and production in the absence of protein. Full results are in Extended Data Fig. 6. The fifth set consisted of one replicate of a simple flat-bilayer simulation corresponding to lipid composition 1. Lipids were allowed to equilibrate between leaflets for several microseconds by opening a pore in the bilayer. See Supplementary Information for a detailed description of this procedure. 1 and 2 (CHMP1B/IST1 WT (1) and F9A + F13A CHMP1B mutant  (2) with composition 1). The structure of two turns of the WT CHMP1B-IST1 protein coat (34 subunits of each protein, PDB 6TZ5) was centered in a 30 × 30 × 18 nm simulation box and converted to a CG Martini 2.2 representation. For the F9A + F13A simulations, the protein structure was first edited to contain the mutations and then coarse grained. The protein was fully position restrained throughout simulation to maintain the cryo-EM resolved protein structure in relation to the membrane. Martini model lipids were chosen to correspond to those used in cryo-EM: SDPC (Martini name PUPC), POPS (Martini name POPS) and CHOL (Martini name CHOL) were directly available from the Martini website. A di-oleoyl PIP 2 model was constructed from the Martini model 1-palmitoyl-2-oleoyl PIP 2 (Martini name POP2) by changing parameters for the second bead of the CG palmitoyl tail from the saturated C type to the unsaturated D type. The total number of lipids was tuned using a set of initial tubule self-assembly simulations as described in Supplementary Information with results in Extended Data Fig. 1, arriving at a total of 1,300 with 754 SDPC, 234 POPS, 234 CHOL and 78 PIP 2 . Details on lipid placement are described in Supplementary  Information. Additionally, the scripts used, as well as a set of example input and output files, for the initial random placement of lipids within the lumen of the protein coat have been deposited on Zenodo (https:// doi.org/10.5281/zenodo.7232344). Instructions for their use are in Supplementary Information. The simulation box was then solvated using standard Martini water and ions to neutralize overall system charge and reach 150 mM NaCl. Simulations were run using Gromacs 2018.8, using Martini-recommended run parameters. Reaction-field electrostatics were used with a 1.1 nm Coulomb cutoff. The velocity rescale thermostat was used to maintain a temperature of 320 K, with separate temperature coupling groups for protein, lipids and solvent. The timestep and pressure control method differed between equilibration and production phases and are described in detail in Supplementary Information. One round of minimization was performed before dynamics, which proceeded in three general steps-spontaneous tubule assembly, leaflet equilibration and production-described in Supplementary Information. The results described in the main text are from analysis on the last 2.4 µs of the production phase of each of the ten replicates for the WT simulations. Three replicates were run for each mutant simulation, with analysis also on the last 2.4 µs of the production phase of each replicate. See Supplementary Table 7 for the leaflet compositions of each replicate. Pymol v.2.3.0 was used to visualize simulation frames.

Set 3 (WT CHMP1B/IST1 with lipid composition 17).
Simulation of WT CHMP1B/IST1 with lipid composition 17 required modifications to the protocol due to instabilities in the membrane that resulted in persistent bubbling across a range of run parameters. These simulations were ultimately run entirely in NVT, as only turning off pressure control reliably avoided bubbling. Full details on the run parameters, other tested combinations and anticipated effects on the results are in Supplementary Information. The local composition results described in main text (Fig. 6) are from analysis on the last 2.4 µs of the production phase of each of the three replicates. See Supplementary Table 7 for the leaflet compositions of each replicate. Set 4 (protein-free tubule simulation). One simulation of a lipid composition 1 tubule of the same curvature without the protein was run. All parameters used and overall setup were identical to simulation set 1, with the exception that isotropic pressure coupling was used, and the leaflet equilibration phase was extended to 30 µs. Production lasted 3 µs, with the last 2.4 µs used for analysis. With semi-isotropic pressure coupling on, and without the protein coat maintaining the tubule shape, the system quickly relaxed to an expanded tubule radius to minimize the curvature energy of the membrane. Therefore, the simulation was run using isotropic pressure coupling, which allowed for fluctuations in the tubule length while maintaining the desired curvature.

Simulation analysis
Leaflet assignment was performed with single linkage agglomerative clustering on the Cartesian coordinates of the phospholipid phosphate beads at every frame, unambiguously producing two clusters corresponding to the inner and outer leaflets ( Supplementary Fig. 2). For each frame, cholesterols were assigned to a leaflet by determining the nearest phospholipid phosphate to the cholesterol hydroxyl group and assigning the sterol to the same leaflet. See Supplementary Table 8 for leaflet and bilayer thicknesses for each simulated condition.
Cylindrical coordinates (r, θ) were determined relative to the long axis of the membrane tubule, which was aligned with the z axis of the box. For each frame we determined the optimal position of the tubule center (x 0 , y 0 ) in the xy plane that optimally placed all CHMP1B F9 + F13 (or F9A + F13A) residues along the protein filament equidistant from the cylindrical center line using a least-squares fitting procedure. All atomic coordinates for that snapshot were then calculated in that reference frame. For the protein-free tubule, the mid-line of the cylinder was taken as the z axis of the box.
Headgroup density calculations were performed using the NC3 bead of SDPC, CNO bead of POPS and the C1 bead of PIP 2 . A 2D Gaussian kernel density estimate on the set of (θ, z) coordinates from each trajectory was then calculated to compute the heatmaps. The grid spacing for the kernel density estimate was set of 0.835 Å in z and 1° in θ, to correspond with the EM density analysis of the tubule surface. To avoid edge effects, positions within the upper and lower 15% of the borders in z and θ were mirrored to effectively expand the box, and the kernel density for that expanded box was calculated and then trimmed back down to the unmirrored box dimensions for plotting. Densities for individual trajectories were subsequently averaged and normalized to span the range from 0 to 1.
Radial density profiles for each lipid bead type were calculated for each leaflet within different zones of the membrane surface relative to the protein coat ( Supplementary Information and Supplementary Fig.  3) using a Gaussian kernel density estimate on the r coordinates from each trajectory. The bandwidth was set to 0.5 Å, and densities were calculated from 0 to 80 Å with a 0.2 Å grid spacing. Densities were then normalized by 2πrΔr to account for the increasing effective volume at greater r, and then multiplied by the number of beads of the specified type within that zone. Thus, densities reflect enrichment or depletion of different lipid species and components. Replicate profiles were averaged to arrive at the final density profile.
Lipid tilt vectors for each tail were calculated by first computing lipid orientation vectors for each tail defined as the vector from the lipid phosphate bead to the last bead of each tail. The bilayer normal was assumed to be the radial component of the cylindrical coordinate system, using +r for the inner leaflet and −r for the outer leaflet. Hence, tilt values of 0° correspond to lipids with configurations normal to the local membrane plane regardless of leaflet assignment. Lipids were assigned to zones based on the location of the phosphate bead. Per zone angle distributions for each lipid type and tail were calculated across each replicate using a Gaussian kernel density estimate spanning 0 to 180° with a bandwidth of 0.05° and a grid spacing of 0.9°. Replicate profiles were averaged to arrive at the final density profile.
Local compositions per zone were calculated by gathering all nonheadgroup beads within each zone and counting the number of beads present from each lipid species. The total number of beads per species was then normalized by the number of beads that make up a whole lipid of each individual lipid type in the model, excluding the headgroups, to get a fractional molar number of each type of lipid present (that is, 12 beads for SDPC, 11 for POPS, eight for CHOL and 11 for PIP 2 ). Headgroups were excluded since they have an outsized effect for PIP 2 lipids (five beads per headgroup versus one for SDPC and POPS). The normalized molar values for each lipid type were then summed and used to calculate local molar compositions per zone and leaflet. Compositions from independent trajectories were averaged to arrive at final reported values.

Reporting summary
Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.

Data availability
The cryo-EM maps have been deposited into the Electron Microscopy Data Bank (accession numbers EMD-27991, EMD-28694-28719 and EMD-28722). The motion-corrected cryo-EM micrographs have been deposited in the EMPIAR database (https://www.ebi.ac.uk/empiar/) with accession number EMPIAR-11277. PDB 6TZ5 was used for simulations and data analysis. Source data are provided with this paper.

Code availability
Code for setting up and analyzing molecular dynamics simulations has been deposited on Zenodo (https://doi.org/10.5281/zenodo.7232344).  , c); whereas an over-packed tubule bulges outward away from the protein (e) or fills the tubule lumen with lipids (f). A lipid count of 1300 (d) was used for production simulations since this resulted in a consistent tubule radius along the z axis with a clearly defined and solvated tubule lumen. These simulations employed a protein coat (PDB ID 3JC1) with opposite handedness to the coat used for all simulations reported in the main text (PDB ID 6TZ5), but nearly identical dimensions. These simulations used an isotropic Berendsen barostat for pressure coupling and were run with a 0.03 ps timestep for 6 µs, while all other parameters were identical to those described elsewhere in Methods). (g, h) Leaflet compositions over time from one replicate of WT CHMP1B-IST1 composition 1 simulation during leaflet equilibration (procedure described in SI Section 13). Leaflet assignment was determined using a simple r cutoff (r < 42 Å inner, else outer) applied to the lipid center of mass, since the clustering algorithm performed poorly in the presence of the equilibration pore due to the connection between the two leaflets. PIP 2 was found to typically require the longest time to equilibrate, likely due to its low total number in the lipid mix and strong electrostatic attraction to the protein coat. Equilibration was run for 9 µs for composition 1 simulations and 6 µs for composition 17 simulations, as composition 17 replicates more quickly converged to similar leaflet compositions. Data for graphs in g, h are available as source data. Fig. 4 | Compositional analysis of lipid bilayers in cryo-EM reconstructions and biophysical characterization of brominated lipids. (a) Analysis of the structure and composition of the membrane-protein contact sites and the gaps between contact sites. To improve the signal to noise ratio, each Z slice of the cryo-EM reconstruction was rotated to align it with the previous frame, resulting in an 'uncoiled' nanotube that was projected in Z. (b) Radial averages for 30 ° at the contact sites (left) and between the contact sites (right), where the protein is farthest from the bilayer. Radial averages for cryo-EM reconstructions at the contact sites (c) and at the gaps between contact sites (d). Cryo-EM reconstructions from samples in which different fractions of the SDPC was brominated. (e) Quantification of the excess intensity in each leaflet, with normalization per leaflet. Reducing the SDPC-Br concentration results in attenuated Br scattering in each leaflet. The reduction is not linear with SDPC-Br concentration, suggesting that SDPC-Br is somewhat enriched in the lipid bilayer nanotubes relative to SDPC. Data are plotted as the mean from independent calculations from each of the two half maps, and error bars represent the standard deviation from the independent half maps. (f) Visualization of the SDPC-Br titration along a Z line scan of the outer leaflet from vertical slices of the cryo-EM reconstructions (g). (h, i) Langmuir pressure-area isotherms show that brominated lipids behave similarly to their unbrominated analogs in lipid monolayers. Although brominated lipids are less compressible than the unsaturated lipids, they are much more compressible than saturated lipids like DSPC. Data for graphs in c-e, g-i are available as source data.

Extended Data
Article https://doi.org/10.1038/s41594-022-00898-1 Extended Data Fig. 5 | Lipid bilayer structure from CG-MD simulations with F9A + F13A CHMP1B and lipid composition 1. (a) Average outer leaflet headgroup density represented as a two-dimensional (θ, z) density from three independent replicates. F9A + F13A positions denoted by green dots. (b) Snapshot of the outer tubule surface, with protein and solvent hidden (color coding same as Extended Data Fig. 1a-f). The smaller alanine side chains appear less able to maintain a persistent hydrophobic defect than the WT (main text Fig. 1h). (c, d) Full radial density (c) and tail angle data (d) at the mutated alanine contact site (left panels) and the gap between the contact site (right panels). Formatting corresponds to Extended Data Fig. 2a, b. At the mutated alanine contact site in the outer leaflet, phospholipid headgroups are pushed inward radially, consistent with a putative elastic deformation of the bilayer, while the polyunsaturated tail of SDPC displays reduced backflipping compared to WT (Extended Data Fig. 2a). For tail angle distributions (d), outer leaflet lipid tails are more tilted at the alanine contact site than elsewhere, though less than at WT F9 + F13 tubule site (Extended Data Fig. 2b). (e) Bilayer and leaflet thicknesses determined from mean radial lipid positions, formatting corresponds to Extended Data Fig. 2c. Bilayer and leaflet thicknesses in the center table; mean locations of midplane and glycerol planes in the right table. Leaflet thickness are nearly identical to the WT composition 1 simulation values (Extended Data Fig. 2c), though each plane is shifted to slightly lower values of r. (f, g) Mean local compositions at the alanine contact site (f) and at the gap between the contact site (g). The ± value is the standard deviation across the 3 independent replicates. (h) Local compositions per leaflet per zone, plotted as the axial z position of the zone to correspond with Fig. 3f (main text). Color coding and formatting identical to Extended Data Fig. 2f. Error bars represent ± 1 standard deviation between replicates. The variations in local composition appear to follow similar trends as in the WT composition 1 simulations (Extended Data Fig. 2f), though at smaller magnitude. Data for graphs in a, c-e, h are available as source data. Snapshot of the outer tubule surface from the last frame from the production stage with solvent hidden (see Extended Data Fig. 1 for color coding). (c) Overall leaflet compositions calculated from production data. Relative to the WT CHMP1B-IST1 lipid composition 1 simulations (Table S7), there is stronger SDPC partitioning to the outer leaflet, stronger CHOL partitioning to the inner leaflet, and an inversion of the PIP 2 partitioning now to the inner rather than the outer leaflet. d-e) Full radial density (d) and tail angle data (e). Radial density profiles correspond to different segments of the lipid, color coded as shown on the CG representations at the left. See Extended Data Fig. 2 for color coding.
Outer leaflet densities are drawn as solid lines; inner leaflet densities are drawn as dashed lines. Note y-axis scale varies by lipid type due to abundance. In the absence of the protein, lipids are relatively normally oriented and there is no extreme backflipping even by the polyunsaturated tail of SDPC. That tail remains slightly shifted away from the bilayer center due to its flexibility but does not show substantial density past the phosphate peak (compare to Extended Data Fig. 2a). Tail angle distributions (e) reflect this as well. Color scheme and line styles follow the same scheme as (d), and all distributions are normalized to 1. (f) Characterization of bilayer and leaflet thicknesses. Bilayer and leaflet thicknesses are shown in the center table, and mean locations of the midplane and glycerol planes are at the right. Data for graphs in a, d-f are available as source data. Fig. 7 | Lipid bilayer structure from CG-MD simulations with wild type CHMP1B-IST1 and lipid composition 17 simulations. (a) Average outer leaflet headgroup patterning represented as a two-dimensional (θ, z) density from three independent replicates. (b) Representative snapshot of the outer tubule surface, with protein and solvent hidden (see Extended Data Fig. 2 for color coding). The scale and persistence of the hydrophobic defect appear similar to the WT simulations using Composition 1 (Fig. 1h). (c, d) Full radial density (c) and tail angle data (d) at the phenylalanine contact site (left panels) and the gap between the contact site (right panels). Radial density profiles correspond to different segments of the lipid, color coded as shown on the CG representations at the left. See Extended Data Fig. 2 for color coding. Outer leaflet densities are drawn as solid lines; inner leaflet densities are drawn as dashed lines. Note y-axis scale varies by lipid type due to abundance. Main structural features seen in the WT protein + lipid composition 1 simulations (Extended Data Fig. 2a, b) are preserved despite the varied composition. Tail angle distributions (d) behave similarly as well. Color scheme and line styles follow the same scheme as (c), and all distributions are normalized to 1. (e) Bilayer and leaflet thicknesses determined from mean radial lipid positions. Bilayer and leaflet thicknesses are summarized in the center table, and mean locations of midplane and glycerol planes are in the right table. The leaflets similarly show differential thinning, and the bilayer is only marginally thicker than that of the WT composition 1 simulations (Extended Data Fig. 2c), perhaps indicating that the reduction in overall SDPC content for composition 17 produces a less flexible, thicker membrane. (f) Local compositions per leaflet per zone, plotted as the axial z position of the zone to correspond with Fig. 3f (main text). See Fig. S3 for zone definitions. Error bars represent ± 1 standard deviation between replicates. Data for graphs in a, c-e, f are available as source data. Inner leaflet values are shown with hollow markers and outer leaflet markers are filled. Horizontal error bars at left are ± 1 standard deviation of the mean time-averaged radial coordinate across ten trajectories, and at right are ± 1 standard deviation of radial coordinates across time in the protein-free trajectory. (d) Diffusion rates in the presence of protein are similar between the WT and mutant CHMP1B simulations, with outer leaflet lipids-especially PIP 2showing dramatically slowed diffusion rates relative to inner leaflet lipids. This reduction in diffusion is likely due to electrostatic interactions with positively charged CHMP1B residues, which transiently pin charged lipids like PIP 2 , and the headgroup exclusion region at the membrane-protein contact site (Fig. 1h), which largely prevents lipids from crossing the helical furrow formed by these residues. The slight increase in diffusion for outer leaflet SDPC and POPS in the F9A + F13A simulations may reflect the lessened headgroup exclusion (Extended Data Fig. 5a). With the protein removed, each lipid shows faster diffusion in the outer leaflet, as expected from the tubule geometry 68 . We expect that the reduction in inner leaflet diffusion rates for the no-protein case relative to the with-protein simulations is attributable to higher amounts of saturated tails in the inner leaflet (Extended Data Fig. 6c vs. Table S7). Faster outer leaflet diffusion in the absence of protein suggests a friction-like effect from the protein. Finally, increased diffusion rates in the outer leaflet of the protein-free tubule relative to the flat bilayer suggests that the increase in area-per-lipid increases diffusion rates. Data for graphs in a-c are available as source data.