<?xml-model href='http://www.tei-c.org/release/xml/tei/custom/schema/relaxng/tei_all.rng' schematypens='http://relaxng.org/ns/structure/1.0'?><TEI xmlns="http://www.tei-c.org/ns/1.0">
	<teiHeader>
		<fileDesc>
			<titleStmt><title level='a'>Secondary Ion Mass Spectrometry of Single Giant Unilamellar Vesicles Reveals Compositional Variability</title></titleStmt>
			<publicationStmt>
				<publisher>Journal of the American Chemical Society</publisher>
				<date>12/20/2023</date>
			</publicationStmt>
			<sourceDesc>
				<bibl> 
					<idno type="par_id">10555419</idno>
					<idno type="doi">10.1021/jacs.3c09039</idno>
					<title level='j'>Journal of the American Chemical Society</title>
<idno>0002-7863</idno>
<biblScope unit="volume">145</biblScope>
<biblScope unit="issue">50</biblScope>					

					<author>Dashiel S Grusky</author><author>Ahanjit Bhattacharya</author><author>Steven G Boxer</author>
				</bibl>
			</sourceDesc>
		</fileDesc>
		<profileDesc>
			<abstract><ab><![CDATA[See Abstract]]></ab></abstract>
		</profileDesc>
	</teiHeader>
	<text><body xmlns="http://www.tei-c.org/ns/1.0" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xmlns:xlink="http://www.w3.org/1999/xlink">
<div xmlns="http://www.tei-c.org/ns/1.0"><head>&#9632; INTRODUCTION</head><p>Giant unilamellar vesicles (GUVs) are a commonly used model system to probe lipid phase behavior, <ref type="bibr">[1]</ref><ref type="bibr">[2]</ref><ref type="bibr">[3]</ref><ref type="bibr">[4]</ref> membraneprotein interactions, <ref type="bibr">[5]</ref><ref type="bibr">[6]</ref><ref type="bibr">[7]</ref><ref type="bibr">[8]</ref><ref type="bibr">[9]</ref><ref type="bibr">[10]</ref> and to encapsulate cellular machinery. <ref type="bibr">[11]</ref><ref type="bibr">[12]</ref><ref type="bibr">[13]</ref><ref type="bibr">[14]</ref> GUVs are attractive model systems not only due to their large size (typically 10-20 &#956;m in diameter), which allows for convenient imaging via optical microscopy <ref type="bibr">[15]</ref><ref type="bibr">[16]</ref><ref type="bibr">[17]</ref><ref type="bibr">[18]</ref><ref type="bibr">[19]</ref><ref type="bibr">[20]</ref><ref type="bibr">[21]</ref> but also due to their ease of production.</p><p>GUVs are commonly produced by either gentle hydration or electroformation. Although both methods start with a lipid mixture dried as a film, for gentle hydration the film is dried onto glass, while for electroformation the film is dried onto either platinum electrodes or indium tin oxide slides. <ref type="bibr">22</ref> Films used for gentle hydration are heated above the melting point of the lipid mixture in the presence of either heated aqueous or heated sucrose solution. <ref type="bibr">[23]</ref><ref type="bibr">[24]</ref><ref type="bibr">[25]</ref> The temperature is maintained above the melting point of the lipid mixture for the duration of GUV formation. Films used for electroformation are also rehydrated, typically at low ionic strength, before alternating current is applied to the film, which assists in the formation of GUVs. <ref type="bibr">[26]</ref><ref type="bibr">[27]</ref><ref type="bibr">[28]</ref> This process is similarly conducted at temperatures above the melting point of the lipid mixture.</p><p>Despite their widespread use, there is some evidence pointing to potential issues with GUVs as model systems. In particular, prior work has shown that electroformed GUVs composed of a ternary mixture exhibit significant variation in areas occupied by an L d partitioning fluorescent dye. <ref type="bibr">3</ref> Other work has shown that individual GUVs produced from the same lipid film show significant variations in their phase behavior. <ref type="bibr">29,</ref><ref type="bibr">30</ref> Although work has been done to compare GUV formation methods in terms of their resulting unilamellarity, capacitance, shear viscosity, and a number of defects visible by fluorescence microscopy, <ref type="bibr">31,</ref><ref type="bibr">32</ref> compositional variability has yet to be directly examined. Quantifying this variability is critical, as lipid composition is the primary variable in all GUV-based measurements. Prior indirect measurements have found compositional standard deviations between 1 and 4 mol %. <ref type="bibr">3,</ref><ref type="bibr">33,</ref><ref type="bibr">34</ref> Here, we directly quantify the compositional variability between individual GUVs using stable isotope labeling and secondary ion mass spectrometry (SIMS). Stable isotope labeling is not only non-perturbative (although perdeuterated lipids exhibit lower transition temperatures) <ref type="bibr">35,</ref><ref type="bibr">36</ref> but also allows for the concentration of labeled lipids to be directly quantified.</p><p>SIMS can be used to measure isotopic ratios present in a sample with high precision. <ref type="bibr">37,</ref><ref type="bibr">38</ref> The Cameca NanoSIMS 50L functions by collisional sputtering of a freeze-dried bilayer containing isotopically labeled lipids with primary cesium ions (Cs + ). This process ejects secondary ions, which are then separated by a mass analyzer, allowing up to seven species to be detected simultaneously. The high sensitivity and mass precision (being able to resolve species such as <ref type="bibr">12</ref> C 1 H -and 13 C -which have similar masses) <ref type="bibr">38,</ref><ref type="bibr">39</ref> of the NanoSIMS 50L allows for compositional information on individual bilayers to be obtained. <ref type="bibr">38,</ref><ref type="bibr">40</ref> For the experiments reported here, GUVs formed via either gentle hydration or electroformation were exposed to NanoSIMS substrates (10 nm SiO 2 -coated Si) where they rupture to form supported bilayer patches whose area is proportional to the surface area of the parent GUV. <ref type="bibr">41,</ref><ref type="bibr">42</ref> This process is depicted in Figure <ref type="figure">1</ref>. In parallel, continuous supported bilayers (SLBs) were formed by conventional small vesicle fusion. SIMS can then be used to examine the compositional variability of the resulting SLBs formed by different methods. The resulting compositional variability is assessed either via external calibration curves, which determine absolute concentrations of isotopically labeled lipids in an SLB, or by examining the relative concentrations of two isotopically labeled lipids. This second method avoids relying on the accuracy of external standards and is not subject to the potential contamination of the NanoSIMS substrate. Both methods demonstrate that patches formed from individual GUVs show composition variation with standard deviations on the order of 1-4 mol % and that GUVs prepared via electroformation are less variable (standard deviations around 1-2 mol %) than GUVs formed via gentle hydration (standard deviations around 2-4 mol %). These measurements are consistent with prior indirect estimates of GUV composition. <ref type="bibr">3,</ref><ref type="bibr">33,</ref><ref type="bibr">34</ref> Additionally, quantification of the average concentration of cholesterol in GUVs formed via gentle hydration and electroformation suggests that GUVs formed by electroformation have significantly less (around 5 mol %) cholesterol on average than GUVs formed via gentle hydration.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>&#9632; MATERIALS AND METHODS</head><p>All natural abundance lipids, cholesterol, and 2 H 31 -POPC (1palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) were purchased from Avanti Polar Lipids. 2 H 7 -cholesterol was purchased from Cayman Chemical. Texas Red 1,2-dihexadecanoyl-sn-glycero-3phosphoethanolamine (TR-DHPE) was purchased from Thermo Fisher Scientific. Four inch &lt;100&gt; p-type silicon wafers (9.5 nm SiO 2 ) were purchased from Silicon Quest International and were diced to 5 &#215; 5 mm to fit in the NanoSIMS sample holder. NanoSIMS substrates were patterned with a chrome grid (5 nm height, 5 &#956;m width) with 25, 50, or 100 &#956;m 2 dimensions via photolithography to facilitate correlative imaging by fluorescence microscopy. All solvents were purchased from Fisher. <ref type="bibr">13</ref> C 18 -POPC, <ref type="bibr">13</ref> C 18 -DSPC (1,2-distearoyl-snglycero-3-phosphocholine), and <ref type="bibr">15</ref> N-POPC were synthesized as previously described. <ref type="bibr">38</ref>, <ref type="bibr">43</ref>  <ref type="bibr">13</ref> C 27 -cholesterol was isolated as previously described. <ref type="bibr">44</ref> Structures for the labeled lipids used in this study can be seen in Figure <ref type="figure">2</ref>.</p><p>Master Stock Solutions. Master stock solutions were prepared with the lipid mixture of interest dissolved in chloroform in sufficient quantities such that multiple batches of GUVs could be produced from the same vial. Master stocks were made by first adding several hundred microliters of chloroform to a 2 mL glass vial. Lipids were then added to the vial by withdrawing the appropriate volume from a pure lipid stock solution in chloroform and then injecting the volume beneath the chloroform into the master stock vial. This ensured that each lipid remained fully dissolved within the master stock. All master stocks also contained 0.1 mol % TR-DHPE so GUV patches could be examined via fluorescence microscopy once ruptured on patterned  substrates. Critically, any comparison between different methods of GUV formation was done by using lipid films produced from one of these master stocks. Also note that all lipid films were dried down at room temperature (23 &#177; 1 &#176;C). Other work has suggested that films dried down at higher temperatures may yield more homogeneous GUVs <ref type="bibr">33</ref> . This effect has not been explored further in the current work.</p><p>Gentle Hydration. Films were dried in 2 mL glass vials from Fisher. First, 200 &#956;L of chloroform was added to a glass vial. 50 nmol of lipid in chloroform was taken from a master stock and added beneath the 200 &#956;L of chloroform in the vial such that the lipid mixture remained dissolved. The glass vial was then vortexed, bath sonicated, and vortexed again for 30 s each. The film was then dried under a stream of argon. The vial was then placed in a desiccator overnight to remove any residual solvent. Films were then rehydrated in 0.5 mL of submicrometer-filtered 500 mM sucrose and heated to 65 &#177; 0.5 &#176;C, above the melting temperature of DSPC (54.4 &#176;C), for 15 h.</p><p>Note that while the heating temperature is above the melting point of the mixture, lower compositional variability may be attained by heating the mixture a further 10 &#176;C above the melting point to approximately 75 &#176;C, as has been noted by prior work. <ref type="bibr">33,</ref><ref type="bibr">45</ref> Electroformation. The electroformation chamber and platinum electrodes were thoroughly cleaned before the lipid films were dried on the electrodes. The chamber was first bath sonicated at 56 &#176;C in a mixture of 7&#215; detergent, ethanol, and deionized water in a 1:3:3 ratio. The setup was then rinsed in deionized water for 20 min before being further bath sonicated in deionized water and rinsed again with deionized water for another 20 min. The chamber was then sonicated in ethanol at room temperature. After the chamber was removed from the ethanol, the setup was dried immediately and kept in a desiccator until use. Films were formed by directly spreading 66 nmol of lipid taken from a master stock onto the platinum electrodes. After lipids were spread on the electrodes, the chamber was kept in a desiccator overnight to remove residual chloroform. The chamber was then sealed with clean glass slides and vacuum grease. Once sealed, the chamber was filled with 1.5 mL of submicrometer-filtered and degassed 500 mM sucrose before being heated to 65 &#177; 0.5 &#176;C. GUVs were electroformed at 10 Hz, 3 V (peak to peak) for 2 h and then at 1 Hz, 3 V for another 30 min.</p><p>Vesicle Extrusion. Glass test tubes were filled with 200 &#956;L of chloroform. 28 nmol of lipid from a master stock was then injected beneath chloroform in the test tube. Lipid films were then dried under argon onto the test tube sides before being desiccated overnight. Films were then resuspended in 1&#215; phosphate-buffered saline (137 mM NaCl, 2.7 mM KCl, 8 mM NaPO 4 , pH 7.2, submicrometerfiltered) and vortexed for 1 min. The buffer with the resuspended lipids was then passed through a membrane with 100 nm pore size 61 times while being heated to 65 &#176;C to form SUVs (small unilamellar vesicles).</p><p>NanoSIMS Sample Preparation. GUVs were ruptured onto silicon substrates to form SLB patches after being allowed to briefly cool for 30 &#177; 5 min to room temperature (23 &#177; 1 &#176;C). Silicon substrates were plasma cleaned for 10 min after which they were submerged in phosphate buffer (240 mM NaCl, 10 mM NaH 2 PO 4 at pH 7.4). Then GUVs were deposited over the submerged substrates and allowed to incubate until approximately 10-15% of the surface was covered in SLB patches. GUV deposition was observed with a Nikon Eclipse 80i epifluorescence microscope equipped with an Andor Clara camera. The substrate and bilayers were then extensively washed with Milli-Q water. Although GUVs, particularly those formed via gentle hydration, can be multivesicular, upon GUV rupture to form an SLB patch, the internal vesicles are liberated and rinsed away. Cleaned substrates were flash frozen in liquid nitrogen and then subjected to low pressure for at least 12 h to sublimate any vitreous ice.</p><p>To form continuous bilayers on NanoSIMS substrates, 100 nm vesicles were incubated over plasma-cleaned substrates for 1 min before being washed extensively with Milli-Q water. Continuous bilayers on substrates were found to be more susceptible to dewetting during flash freezing so substrates with continuous bilayers were removed from Milli-Q water with the bilayer facing upside-down. This ensured that a drop of water remained in contact with the bilayer at all times before the substrate was flash frozen in liquid nitrogen. Once frozen, substrates were subjected to low pressure using the same method as for GUV patch samples. Lipid monolayers were formed using a KSV NIMA KN 2002 (Biolin Scientific, Stockholm, Sweden) Langmuir trough (273 cm 2 ) at room temperature (23 &#177; 1 &#176;C). Whatman filter paper was used as a Wilhelmy plate to monitor surface pressure. Lipid mixtures dissolved in chloroform were spread on water (&gt;18 M&#937; from Milli-Q system) within the clean trough using a glass syringe. The chloroform was left to evaporate for 10 min and the barriers were compressed at 10 mm/min until the surface pressure reached 32 mN/m. Plasma-cleaned NanoSIMS substrates were glued to a glass slide and pulled through the air-water interface at a rate of 1 mm/min, while the surface pressure was maintained at 32 mN/m. Lipid monolayers were not subjected to freeze-drying, as unlike SLBs, lipid monolayers are stable in air. Both lipid monolayers and freezedried lipid bilayers were stored in a desiccator when not being analyzed via NanoSIMS.</p><p>NanoSIMS Analysis. Analysis was performed on a Cameca NanoSIMS 50L instrument at Stanford University. Correlative fluorescence imaging helped facilitate the selection of bilayer patches so that debris on the substrate surface were avoided. Bilayer patches were imaged with a 2pA 133 Cs + primary beam. Ten 25 &#215; 25 &#956;m scans (256 &#215; 256 pixels, 1 ms dwell time) were collected, which is enough to remove all of the deposited material on the surface. Secondary ion detectors were set to 2 H -, 12 C -, 13 C -, <ref type="bibr">12</ref> C 2 H -, 12 C 2 2 H -, <ref type="bibr">13</ref> C 2 H -, and 13 C 2 2 H -for samples containing 2 H-or <ref type="bibr">13</ref> C-labeled lipids or cholesterol. If 15 N-labeled POPC was contained within SLB patches, secondary ion detectors were set to 2 H -, 12 C -, 13 C -, <ref type="bibr">12</ref> C 2 H -, 12 C 14 N -, 13 C 15 N -, and <ref type="bibr">13</ref> C 2 2 H -. Standard samples (the calibration curves described below) were regularly analyzed to ensure that isotope ratios were reproducible from session to session.</p><p>Data Analysis. Images were analyzed using ImageJ (National Institutes of Health, USA) with the OpenMIMS plugin (National Resource for Mass Spectrometry, Harvard University USA). Planes were summed, and regions of interest were manually selected in order to exclude any debris on the sample. Total counts within each region of interest were determined via the "Tomography" tab. These counts were then used to determine the ratios of interest (typically 2 H -/ 13 C -, 13 C -/( 13 C -+ 12 C -), and 2 H -/( 13 C -+ 12 C -)). Calculating these ratios allows for the size of the analyzed patches to be taken into account and allows for further quantification via external calibration standards.</p><p>Calibration Standards. Concentration calibration standards were made from lipid mixtures dissolved in chloroform containing a known mol % of labeled lipid. Calibration standards used to quantify labeled lipid concentrations in ternary SLBs contained 20 mol % cholesterol, as this is present in the ternary mixture. Calibration curves without cholesterol were also prepared so that labeled lipid concentrations in mixtures without cholesterol could quantified. The prepared lipid mixtures dissolved in chloroform were spread on plasma-cleaned NanoSIMS substrates. The chloroform was then allowed to evaporate to form a lipid film. Concentration calibration standards were kept in a desiccator until use.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>&#9632; RESULTS AND DISCUSSION</head><p>Comparing GUV Preparation Methods via Double Labeling. First, we assessed the variability in relative concentrations of labeled lipids for different GUV preparation methods. Although direct concentration quantification via external calibration curves can give estimates of the concentration of a labeled species, this method is subject to inaccuracies. These experiments can be skewed by surface contamination and are heavily reliant on the accuracy of the calibration standards. In order to avoid these issues, and inspired by prior work, <ref type="bibr">46</ref> initial experiments examined the relative change in ion counts from patch to patch resulting Journal of the American Chemical Society from two lipids with different isotopic labels. This was done with a POPC:DSPC:CHOL 40:40:20 mixture. This mixture was chosen as GUV patches with this composition do not display macroscopic phase separation between lipid components <ref type="bibr">1,</ref><ref type="bibr">4,</ref><ref type="bibr">16</ref> within the 50 nm lateral resolution of the NanoSIMS primary ion beam (Figure <ref type="figure">S1</ref>; there is nanoscale separation present in these SLB patches which can be detected by ion recombination <ref type="bibr">43</ref> ). This mixture is also well-studied in GUVs and GUV patches. <ref type="bibr">4,</ref><ref type="bibr">43,</ref><ref type="bibr">47</ref> In order to determine which pair of labeled lipids is most sensitive to relative changes in concentration, two ternary samples, one with the composition 13 C 18 -DSPC:POPC: 2 H 31 -POPC:CHOL 20:20:20:20:20 and the other with the composition DSPC: <ref type="bibr">1 3</ref> C 1 8 -POPC: 2 H 3 1 -POPC:CHOL 40:20:20:20 were prepared. Additionally, a sample with the composition 2 H 31 -POPC: <ref type="bibr">13</ref> C 18 -POPC:POPC 20:20:60 was also prepared (samples with this composition are denoted as pure POPC). <ref type="bibr">13</ref> C -/( 13 C -+ 12 C -) and 2 H -/( 13 C -+ 12 C -) ratios were then measured for 30 GUV patches formed via gentle hydration for all three samples.</p><p>Figure <ref type="figure">3</ref> shows the correlation between the 2 H -/( 13 C -+ 12 C -) and <ref type="bibr">13</ref> C -/( 13 C -+ 12 C -) ratios for each labeling scheme. The 2 H -/( 13 C -+ 12 C -) ratio tracks the amount of 2 H-labeled lipid in the bilayer, while the 13 C -/( 13 C -+ 12 C -) ratio tracks the amount of <ref type="bibr">13</ref> C-labeled lipid in the bilayer. Figure <ref type="figure">3</ref> demonstrates that the correlation between the ratios is a function of both the location of the isotopic label and the overall composition of the bilayer. The ternary mixture containing both 2 H 31 -POPC and 13 C 18 -POPC displays a positive correlation between the lipid concentrations. Conversely, the ternary mixture containing 2 H 31 -POPC and <ref type="bibr">13</ref> C 18 -DSPC displays a clear negative correlation between the Correlations between the 2 H -/( 13 C -+ 12 C -) and <ref type="bibr">13</ref> C -/( 13 C -+ 12 C -) ratios for each lipid mixture. The 2 H -/( 13 C -+ 12 C -) ratio tracks the amount of 2 H-labeled lipid in the bilayer, while the 13 C -/( 13 C -+ 12 C -) ratio tracks the amount of <ref type="bibr">13</ref> C-labeled lipid in the bilayer. (A) Correlation between 2 H 31-POPC and 13 C 18 -DSPC ratios in the ternary mixture DSPC: <ref type="bibr">13</ref>   Measured 2 H -/ 13 C -ratios for different methods. GUV patches formed by gentle hydration and electroformation were compared to continuous bilayers formed from SUVs via their 2 H -/ 13 C -ratios. All samples were produced from a master stock with nominal composition DSPC: 13 C 18 -DSPC:POPC: 2 H 31 -POPC:CHOL 20:20:20:20:20. (A) Displays the 2 H -/ 13 C -ratio distributions for each preparation method. Significance was determined via an F-test. For this and subsequent plots, each point represents a measurement made on a single GUV patch or corral containing an SLB. Thirty bilayers were examined for each sample. For this and all subsequent plots, * p &#8804; 0.05, ** p &#8804; 0.01, *** p &#8804; 0.001, and **** p &#8804; 0.0001. The displayed curves are normal distributions calculated by using the standard deviation and average from the GUV patch measurements. Note that these have been overlaid to guide the reader and that the area under the curve has not been normalized. A dot plot representation of the distributions can be seen in Figure <ref type="figure">S21</ref>. (B) Calculated standard deviations of the 2 H -/ 13 C -ratio for each preparation method.</p><p>concentrations of the two components. These trends can be attributed to preferential interactions between <ref type="bibr">13</ref> C 18 -POPC and 2 H 31 -POPC and unfavorable interactions between <ref type="bibr">13</ref> C 18 -DSPC and 2 H 31 -POPC within ternary bilayers.</p><p>As shown in Figure <ref type="figure">3C</ref>, pure POPC SLBs containing 2 H 31 -POPC and 13 C 18 -POPC did not show any clear correlation in concentration between the differently labeled lipids. 2 H 31 -POPC, <ref type="bibr">13</ref> C 18 -POPC, and natural abundance POPC within the pure POPC mixture only differ from each other in terms of isotopic labeling. As a result, when the concentration of one labeled component is higher in a GUV, the extent to which it replaces either the other labeled or unlabeled component is essentially random.</p><p>Based on the results from Figure <ref type="figure">3</ref>, a ternary mixture containing <ref type="bibr">13</ref> C 18 -DSPC and 2 H 31 -POPC is expected to be the most sensitive to relative concentration changes between the two labeled lipids, as the concentrations of these two components are anticorrelated. Therefore, this mixture was used to study the compositional variability of GUVs formed by electroformation or gentle hydration.</p><p>GUVs with this composition were generated from the same master stock by either gentle hydration or electroformation. The same master stock was also used to generate 100 nm SUVs which were then ruptured onto NanoSIMS substrates to form continuous bilayers within the corrals of the patterned NanoSIMS substrate. Since these continuous bilayers are formed from hundreds of SUVs, the relative concentrations of 13 C 18 -DSPC and 2 H 31 -POPC should be more consistent than in SLB patches formed from the rupture of a single GUV. The measured 2 H -/ 13 C -ratio for 30 individual bilayer patches or 30 corrals containing continuous bilayers is shown in Figure <ref type="figure">4</ref>. Based on the measured ratios, gentle hydration has the highest standard deviation, &#177;1.6, in relative concentration. The standard deviation for GUVs formed via electroformation is significantly (determined via F-test) lower at &#177;1.1. The normality of each distribution assessed by F-test as is further discussed in Section 2 of the Supporting Information. Continuous bilayers formed from SUVs have the lowest standard deviation in relative concentration, at &#177;0.7. Replicate samples produced from films dried from the same master stock also suggest that electroformation is less compositionally variable (lower measured standard deviations) than gentle hydration (Figure <ref type="figure">S4</ref>).</p><p>The raw data for all of the distributions in the main text and the Supporting Information can be seen in Section 16 of the Supporting Information. Additionally, each distribution presented is replotted as a dot plot in Section 17 of the Supporting Information to provide a second way to visualize the data.</p><p>Absolute Concentration Variability. Although measuring the 2 H -/ 13 C -ratio is useful for comparing the variability in different methods while avoiding concerns regarding surface contamination and calibration accuracy (further discussed in Section 4 of the Supporting Information), quantification of absolute concentration is useful for determining how much the mol % of a particular lipid varies from GUV to GUV. Therefore, external calibration curves, such as those shown in Figure <ref type="figure">5A</ref>, were used to relate quantitative ion ratios to labeled lipid concentrations. This allowed the concentration of each labeled species within one GUV patch to be determined. This analysis was performed for the GUV samples discussed in Figure <ref type="figure">4</ref>.</p><p>As shown in Figure <ref type="figure">5B</ref>,C, GUV patches produced via electroformation showed consistently lower standard deviations in both <ref type="bibr">13</ref> C 18 -DSPC and 2 H 31 -POPC concentrations than those formed via gentle hydration as demonstrated by the calculated standard deviations and supported by a corresponding F-test. To determine the lower bound on variability (i.e., how much of the variation in concentration measurements is due to instrument noise and surface contamination), a monolayer with overall composition 2 H 31 -POPC: 13 C 18 -POPC:-POPC 20:20:60 (pure POPC composition) was formed on a NanoSIMS substrate and compared to both GUV samples. Since a monolayer should be compositionally homogeneous across the substrate surface, this measurement can be used to approximate how much of the variability (the measured standard deviation) is due to the experimental method. As seen in Figure <ref type="figure">5C</ref>, measurements on pure POPC monolayers resulted in a lower standard deviation in 2 H 31 -POPC concentrations than that measured in GUVs formed via electroformation or gentle hydration (&#177;1.4 mol % for monolayers, &#177;2.5 mol % for electroformed GUVs, and &#177;3.7 mol % for GUVs formed by gentle hydration). However, these monolayers cannot be compared via the 2 H -/ 13 C -ratio as was done in Figure <ref type="figure">4</ref> because the correlation between 2 H 31 -POPC and 13 C 18 -POPC concentrations in pure POPC is dramatically different from the correlation between 2 H 31 -POPC and 13 C 18 -DSPC in ternary mixtures (Figure <ref type="figure">3</ref>). Additionally, ternary monolayers formed from the ternary master stock containing 2 H 31 -POPC and 13 C 18 -DSPC display macroscale separation (Figure <ref type="figure">S6</ref>) and are therefore not comparable to the GUV samples, as is further discussed in Section 5 of the Supporting Information. Therefore, monolayers can be compared to GUV samples only via absolute concentrations. Additional analysis suggests that the size of the region selected for quantification within the 25 &#215; 25 &#956;m analysis region does not significantly impact the distributions shown in Figures <ref type="figure">4</ref> and <ref type="figure">5</ref> (further discussed in Section 6 of the Supporting Information). Therefore, the intrinsic signal-to-noise ratio of a single patch is unlikely to substantially impact the GUV variability observed for different formation methods. It should be noted that the absolute values measured for isotope-labeled lipid concentrations demonstrate deviation from the nominal concentrations in the master stock solution. In some cases, the labeled lipid concentration deviates by 5-6 mol % from the expected concentration. This is likely the result of multiple factors including substrate surface contamination by <ref type="bibr">12</ref> C, the accuracy of the external calibration curves, and error involved in preparing labeled lipid stock solutions (which are then used to make master stocks).</p><p>Sources of Variability. In order to further explore sources of the observed GUV variability, GUVs with pure POPC composition ( 2 H 31 -POPC: <ref type="bibr">13</ref> C 18 -POPC:POPC 20:20:60) were formed by either gentle hydration or electroformation. Thirty GUV patches with this composition were analyzed for each method, and their absolute concentrations were determined via external calibration curves. In Figure <ref type="figure">6</ref>, the calculated 2 H 31 -POPC concentrations in pure POPC samples are compared to the 2 H 31 -POPC concentrations in the previously discussed ternary mixture. Ternary GUV patches formed either by electroformation or gentle hydration have higher measured standard deviations in 2 H 31 -POPC concentrations relative to pure POPC patches formed by the same method (Figure <ref type="figure">6C</ref>). This lower variability for the pure POPC GUVs can also be seen for the <ref type="bibr">13</ref> C 18 -labeled lipids (Figure <ref type="figure">S9</ref>). These results suggest that more complex lipid compositions lead to considerably more compositional variability. It also suggests that the observed compositional variability is not solely due to residual contamination of either the platinum electrodes used for electroformation or the glass vials used for gentle hydration. If these surfaces had significant contamination, pure POPC patches would have compositional variability comparable to that of ternary patches. Comparison of 2 H 31 -POPC concentration variability between pure POPC patches formed via either electroformation or gentle hydration did not show a significant difference in compositional variability (Figure <ref type="figure">S10</ref>). This suggests that the higher standard deviations observed in GUV patches formed via gentle hydration (Figures <ref type="figure">4</ref> and <ref type="figure">5</ref>) are not due to higher residual contamination on glass vials relative to the platinum electrodes but rather result from the complexity of the ternary mixture and the method of GUV formation. Additionally, the average 2 H -, 12 C -, and <ref type="bibr">13</ref> C - counts were compared between regions of interest with and without a bilayer to gauge the level of contamination on the substrate surface (Figure <ref type="figure">S11</ref>). Regions of interest containing bilayer showed considerably higher signal on all detectors than regions of interest in an exposed substrate. Further analysis was also performed to correct the observed variabilities in ternary GUV patches for noise due to sample preparation and analysis (discussed in Section 10 of the Supporting Information). Both pure POPC and ternary GUVs were formed by either electroformation (A) or gentle hydration (B). For both methods, the patches composed of pure POPC display significantly less 2 H 31 -POPC concentration variability than that of ternary patches formed using the same method. Dot plot representations of these distributions can be seen in Figure <ref type="figure">S23</ref>. (C) Displays the calculated standard deviations for each sample.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>Journal of the American Chemical Society</head><p>To further examine potential sources of variability, lipid films from the master stock containing 13 C 18 -DSPC and 2 H 31 -POPC were dried directly on NanoSIMS substrates and imaged. No significant separation within the 50 nm lateral resolution of the NanoSIMS was observed between 13 C 18 -DSPC and 2 H 31 -POPC within the film (Figure <ref type="figure">S13</ref>). However, atomic recombination experiments demonstrated that there is nanoscale separation between POPC and DSPC within the film (Figure <ref type="figure">S14</ref>). These preferential interactions within the film may contribute to lipid sorting while GUVs are being formed.</p><p>Additional experiments also suggest that the size of the GUV patch analyzed is not correlated with either 13 C 18 -DSPC or 2 H 31 -POPC concentrations (Section 12 of the Supporting Information). Furthermore, additional data suggest that two different GUVs typically do not rupture to form one patch (Section 13 of the Supporting Information). This suggests that most patches are the product of a single GUV and that minimal compositional averaging between GUVs occurs.</p><p>The measured standard deviations as well as the corrected standard deviations (Section 10 of the Supporting Information) agree relatively well with indirect measurements of GUVto-GUV compositional variation. Prior work using fluorescence microscopy to quantify the area fraction occupied by optically resolvable domains approximated a &#177;2.1 mol % standard deviation in electroformed DOPC/eSM/CHOL GUVs. Similarly, estimates from GUV transition temperatures measure a standard deviation of approximately &#177;2 45 or &#177;4 mol % for electroformed GUVs. These results are in relatively good agreement with the direct measurements presented here, where we measure standard deviations around 1-2 mol % for electroformed GUVs and 2-4 mol % for GUVs formed via gentle hydration. However, it should be noted that we used DSPC for this study, which is a higher melting point lipid and may lend itself to slightly greater variability due to the high temperatures (T m = 54.4 &#176;C) needed to keep this lipid fluid.</p><p>Quantifying Cholesterol Concentrations. Cholesterol concentration was also examined in the GUV patches. Unlike phospholipids, Cholesterol localized significantly to the edges of GUV patches formed by either electroformation or gentle hydration as seen in Figures <ref type="figure">7A</ref> and <ref type="figure">S18</ref>. This observation is consistent with work from other groups, who have noted that L o domains can localize to the edge of GUV patches <ref type="bibr">42,</ref><ref type="bibr">48</ref> although other work has suggested that L d domains can also partition to patch edges, <ref type="bibr">49</ref> suggesting that this effect may depend on the composition examined. However, for the composition examined here, we observe cholesterol (presumably in L o domains) partitioning to the edges of bilayer patches. This partitioning makes quantification of cholesterol variability in GUVs significantly more challenging as the relative ratio of edge to center within the analyzed region needs to be considered as well as how the overall concentration of cholesterol in a patch may affect its partitioning between the edge and center. As a result, the variability in cholesterol concentration from GUV to GUV was not assessed. Instead, the average cholesterol concentration was determined for different GUV formation methods. Ternary GUV patches were formed via gentle hydration and electroformation using a master stock with nominal composition DSPC:POPC: 2 H 3 1 -POPC: 1 3 C 2 7 -CHOL 40:20:20:20. The average concentration was then calculated by either including or excluding the cholesterol-rich edges. Regardless of the analysis method used, the average concentration of cholesterol was lower in electroformed GUV patches, as can be seen in Figure <ref type="figure">7C</ref>. The average 2 H 31 -POPC concentration was the same between the two methods regardless of the analysis method (Figure <ref type="figure">S19</ref>). Average cholesterol concentration was also examined in GUV patches where all three components are isotopically labeled (Figure <ref type="figure">S18</ref>). These experiments also demonstrate lower average cholesterol concentrations in the electroformed GUVs.</p><p>Although there appears to be cholesterol partitioning to the edge of bilayer patches, the increase in cholesterol concentration when the edge is included in the analysis relative to when it is excluded appears to be moderate. There is only a statistically significant difference between the analysis methods (p = 0.03), for GUVs formed via gentle hydration (Figure <ref type="figure">7C</ref>). Replicate samples (Figure <ref type="figure">S20</ref>) do not reproduce this </p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>Journal of the American Chemical Society</head><p>difference between the analysis methods. These replicates are of critical importance, as they suggest that the changes in cholesterol concentration are likely not an artifact resulting from the challenges associated with measuring absolute concentrations, as discussed at the end of the previous section.</p><p>The source of the difference in cholesterol concentrations between GUVs formed by electroformation and those formed by gentle hydration may be the result of the alternating current applied during electroformation. While the phospholipids present in the ternary mixture are zwitterionic and potentially more affected by the alternating current applied during electroformation, cholesterol is neutral and therefore may be less responsive to the applied current. This would result in a lower incorporation into electroformed GUVs. However, the experiments performed here cannot provide a definitive explanation or mechanism for the difference in cholesterol concentrations between the two methods. Additionally, it should be noted that prior measurements suggest that hydrated films <ref type="bibr">50,</ref><ref type="bibr">51</ref> have comparable cholesterol solubility limits to electroformed films. <ref type="bibr">52</ref> These results imply that both methods lead to equivalent cholesterol concentrations and disagree with the results presented here.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>&#9632; CONCLUSIONS</head><p>GUVs are a widely used model system for probing lipidprotein and lipid-lipid <ref type="bibr">53,</ref><ref type="bibr">54</ref> interactions. Despite the widespread use of GUVs, little work has been done to probe GUVto-GUV compositional variation. This is likely due to the lack of methods to accurately assess the concentration of a given lipid within a single GUV. The Cameca NanoSIMS 50L allows for high-precision determination of the concentration of individual lipid species via non-perturbative stable isotope labeling.</p><p>Here, we demonstrate that variability on the order of 1-4 mol % is present in GUVs composed of a ternary DSPC:POPC:CHOL mixture. It is shown that GUVs formed via electroformation have considerably less compositional variability than those formed via gentle hydration (GUVs formed by electroformation have concentration standard deviations of 1-2 mol % compared to the 2-4 mol % of GUVs formed via gentle hydration). This is true regardless of whether the relative change in concentration between two labeled species is calculated or if the variability in the absolute concentration is determined via external calibration curves. Although the mechanism behind the lower variability seen in electroformed GUVs is unclear, it is clear that ternary mixtures are far more variable than pure mixtures and that preferential interactions between certain lipids are present in the films used to form ternary GUVs.</p><p>Although the variability in cholesterol concentration is not examined here, the average concentration of labeled cholesterol was compared between different methods. This demonstrated that electroformed GUVs have a lower average cholesterol concentration. Additionally, it is worth noting that other potential disadvantages have been reported with electroformation that are not discussed here. <ref type="bibr">55</ref> Therefore, while electroformation may yield more compositionally uniform GUVs, there are downsides to the method that must be considered.</p><p>&#9632; ASSOCIATED CONTENT * s&#305; Supporting Information The Supporting Information is available free of charge at <ref type="url">https://pubs.acs.org/doi/10.1021/jacs.3c09039</ref>. Choice of lipid mixture; normality verification for Ftests; replicate electroformation and gentle hydration samples; verification of absolute calibration accuracy; ternary monolayers; effect of analysis area on concentration distributions; 13 C-labeled lipid comparisons in pure and ternary mixtures; comparison of pure POPC mixtures with different methods; detector counts on bilayer and exposed substrates; correcting ternary GUV compositional variability; nanoscale heterogeneity in lipid films; relation between GUV patch size and 13 C 18 -DSPC and 2 H 31 -POPC concentrations; assessing the extent of GUV mixing during patch formation; tracking cholesterol concentration in triply labeled bilayers; assessing average 2 H 31 -POPC concentration in 13 C 27 -CHOL containing patches; tabulated raw GUV composition data; dot plot representations of GUV composition distributions (PDF) &#9632; AUTHOR INFORMATION Corresponding Author Steven G. Boxer -Department of Chemistry, Stanford University, Stanford, California 94305-5012, United States; orcid.org/0000-0001-9167-4286; Email: sboxer@ stanford.edu</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head n="1.">Choice of Lipid Mixture</head><p>For this study, a ternary mixture with an overall composition of 40:40:20 DSPC:POPC:CHOL was used. The ternary phase diagram for this mixture is shown in Figure <ref type="figure">S1</ref>. Measurements for this diagram were made at 23&#176;C. This mixture was chosen because it is well-studied in GUVs 1 and in bilayer patches formed from GUVs on SiO 2 2,3 . Additionally, since the mixture does not display macroscopic phase separation within the 50nm lateral resolution of the NanoSIMS, bilayer patches formed from GUVs appear uniform, simplifying quantitative analysis. It should be noted that although macroscopic domain formation is not apparent, nanoscale domain formation has been documented in this mixture using methods that surpass both the optical diffraction limit and the lateral resolution of the NanoSIMS <ref type="bibr">[1]</ref><ref type="bibr">[2]</ref><ref type="bibr">[3]</ref> . Note that such phase diagrams are based on the assumption that the composition of GUVs match the nominal composition of the components dissolved in organic solvents prior to GUV formation. Note that while the nominal composition is placed within the L d + L o coexistence region of the diagram (denoted by the red star in figure <ref type="figure">S1</ref>). As a result, some GUVs may contain some degree of solid L&#946;, which although not visible via NanoSIMS, may still be present and have implications for GUV deposition, and corresponding compositional variability.  Data compared via F-test were assessed for normality, as a fundamental assumption of F-tests is that the data examined should be normally distributed. The normality of a given distribution can be examined via quantile-quantile plots (QQ-plots), which plot experimental values against predicted values assuming a normal distribution. Several sample QQ-plots can be seen in Figure <ref type="figure">S2</ref>, and they display a linear relationship between actual and predicted values, indicating that the data are normally distributed. Plots were made for all samples, with the majority being highly linear, similar to the sample plots shown in Figure <ref type="figure">S2</ref>. The only distribution with visible deviation from linearity in the QQ-plot is the 13 C 18 -DSPC concentration distribution (See figure <ref type="figure">5B</ref>) for ternary DSPC/POPC/CHOL GUVs formed via gentle hydration. The plot for this distribution can be seen in Figure <ref type="figure">S3</ref>, which suggests that the distribution is S5 somewhat right skewed. It should therefore be noted that the application of an F-test to this distribution may be less valid than for the remaining data sets. </p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head n="3.">Replicate Electroformation and Gentle Hydration Samples</head><p>Replicate electroformation and gentle hydration samples were examined to see if the difference in variability between the two methods was reproducible. The results from these replicate samples can be seen in Figure <ref type="figure">S4</ref>. These samples were produced using the same master stock as the samples shown in Figure <ref type="figure">4</ref> of the main text. Relative concentration changes were assessed via the 2 H -/ 13 C -ratio, which again demonstrated significantly less variability in electroformed GUV patches. Absolute concentration quantification via external calibration curves showed that the 2 H 31 -POPC concentration was also more homogenous in electroformed GUV patches, although the 13 C 18 -DSPC concentration did not display a significant difference in variability between the two methods. </p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>S6</head></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head n="4.">Verification of Absolute Calibration Accuracy</head><p>The accuracy of the external calibration curves for <ref type="bibr">13</ref> C 18 -POPC and 13 C 18 -DSPC was estimated by comparing measured <ref type="bibr">13</ref> C -/( 13 C -+ 12 C -) ion ratios to expected ratios. Expected ratios were calculated by determining the number of <ref type="bibr">13</ref> C and <ref type="bibr">12</ref> C carbons present in a sample. For example, in a sample containing 10 mol % 13 C 18 -POPC (remainder being natural abundance POPC), out of 100 lipid molecules 10 of them are <ref type="bibr">13</ref> C 18 -POPC, leading to a total of 180 intentionally <ref type="bibr">13</ref> C-labeled carbons. Additionally, 1.1% of the remaining carbons will be <ref type="bibr">13</ref> C at natural abundance, adding another 44.2 <ref type="bibr">13</ref> C carbons. Dividing 224.2 by the total number of carbons (4200), yields the predicted 13 C -/( 13 C -+ 12 C -) ion ratio of 0.0533. These calculations were carried out for calibration curves containing varying concentrations of 13 C 18 -POPC and <ref type="bibr">13</ref> C 18 -DSPC with the resulting values tabulated next to experimental values in Figure <ref type="figure">S5</ref>. Figure <ref type="figure">S5</ref> demonstrates that there is reasonable agreement between the predicted values and the experimental data. </p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>S7</head></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head n="5.">Ternary Monolayers</head><p>Ternary monolayers made using the 13 C 18 -DSPC and 2 H 31 -POPC ternary master stock were formed on a LB trough as described in the materials and methods. The advantage of examining a ternary monolayer with this composition is that it could be more homogenous than a GUV sample and it would have the same anti-correlation between <ref type="bibr">13</ref> C-and 2 H-labeled lipids as the GUV samples analyzed in Figure <ref type="figure">3</ref> and <ref type="figure">4</ref> of the main text. This would allow for the ternary monolayer to be compared to the GUV samples via the 2 H -/ 13 C - ratio. However, this comparison was not made due to macroscale separation present in ternary monolayers with this composition, as can be seen in Figure <ref type="figure">S6</ref>. This separation between saturated and unsaturated lipids is likely due to the fact that these monolayers are not heated while being produced and therefore do not allow for significant mixing between lipids. Due to this macroscale separation, the measured 2 H -/ 13 C -ratio will depend heavily on the relative area occupied by the POPC and DSPC rich regions within a given S8 analysis region. Due to this confounding variable, a direct comparison between ternary monolayers and ternary GUVs was not made. </p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head n="6.">Effect of Analysis Area on Calculated Distributions</head><p>A potential concern with the methodology presented here is that smaller GUV patches may have intrinsically lower signal to noise. To assess if the size of a given GUV patch, and therefore the size of the region analyzed, influences the resulting concentration distributions, the gentle hydration and electroformation samples discussed in Figures <ref type="figure">4</ref> and <ref type="figure">5</ref> of the main text were further investigated. GUV patches were re-analyzed using a fixed region of interest so that any effect due to the size of analysis region on the resulting measurements is controlled for. A 30x30 pixel region of interest (2.9x2.9 &#956;m) was defined for each GUV patch and the corresponding 2 H -/ 13 C -ratios as well as 2 H 31 -POPC and 13 C 18 -DSPC concentrations were calculated (the latter two used the same set of calibration curves shown in Figure <ref type="figure">5</ref> of the main text to determine labeled lipid concentrations) within this region. This was done for each of the GUV patches analyzed in Figures <ref type="figure">4</ref> and <ref type="figure">5</ref> of the main text. The resulting fixed area distributions were then compared to the distributions calculated using the whole GUV patch (the same distributions shown in Figures <ref type="figure">4</ref> and <ref type="figure">5</ref> of the main text), with the results shown in Figure <ref type="figure">S7</ref>. For GUV patches formed by either gentle hydration or electroformation, there is not a significant difference (determined via F-test) between distributions calculated using the full GUV patch and those calculated using the fixed 2.9x2.9 &#956;m region of interest. This suggests that the total area selected for 2 H -/ 13 C -and labeled lipid concentration analysis does not have a significant impact on the measured distributions and their corresponding variability. However, there is a consistent increase in the standard deviations of distributions calculated with the fixed area S9 analysis relative to those where the full patch was used (Figure <ref type="figure">S7 G</ref> and <ref type="figure">H</ref>). Although this change is not statistically significant, it does suggest that using smaller analysis regions typically leads to slight increases in variability. This is presumably due to smaller analysis regions having less signal. GUV patches formed by either gentle hydration or electroformation were further compared using the fixed area method as shown in Figure <ref type="figure">S8</ref>. The fixed area method recapitulates the result that GUVs formed via electroformation have less compositional variability than those formed via gentle hydration. This is evident using both relative changes in labeled lipid concentration ( 2 H -/ 13 C -ratio) and absolute labeled lipid concentrations ( 13 C 18 -DSPC and 2 H 31 -POPC concentrations).</p><p>Figure <ref type="figure">S7</ref>. Effect of analysis area on measured distributions. 2 H -/ 13 C -ratios, <ref type="bibr">13</ref> C 18 -DSPC concentrations and 2 H 31 -POPC concentrations were determined within a fixed 30x30 pixel region (2.9x2.9 &#956;m). This analysis was performed on the same GUV patches analyzed in Figures <ref type="figure">4</ref> and <ref type="figure">5</ref> of the main text. The distributions calculated using a fixed area are plotted next to distributions calculated for the same GUV patches but using their full area. The resulting distributions for GUVs made via gentle hydration (A-C) and electroformation (D-F), with corresponding standard deviations seen in (G) and (H) are displayed. The variability in the distributions remained the same (determined via F-test) regardless of which analysis method was used. Dot plot representations of these distributions can be seen in Figure <ref type="figure">S25</ref>.  <ref type="figure">4</ref> and <ref type="figure">5</ref> of the main text) had their 2 H -/ 13 C -ratios, <ref type="bibr">13</ref> C 18 -DSPC concentrations and 2 H 31 -POPC concentrations calculated within a fixed 30x30 pixel area (2.9x2.9 &#956;m). The calculated 2 H -/ 13 C -ratio and concentration distributions are compared for each method (A-C) with the corresponding standard deviations shown in (D). Dot plot representations of these distributions can be seen in Figure <ref type="figure">S26</ref>.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head n="7.">13 C-labeled Lipid Comparisons in Pure and Ternary Mixtures</head><p>A pure POPC mixture ( 13 C 18 -POPC: 2 H 31 -POPC:POPC 20:20:60) was used to form GUVs either via electroformation or gentle hydration. In addition to comparing 2 H 31 -POPC concentration variability in the pure POPC sample to that of the ternary mixture (Figure <ref type="figure">6</ref> of the main text) the concentration variability of the 13 C 18 -labeled lipid ( 13 C 18 -POPC for the pure POPC sample and 13 C 18 -DSPC in the ternary mixture) was also compared. This comparison can be seen in Figure <ref type="figure">S9</ref> for GUVs formed via gentle hydration and electroformation. In both cases, the pure POPC sample displayed significantly less variability, recapitulating the results in Figure <ref type="figure">6</ref> of the main text.  <ref type="bibr">13</ref> C-labeled lipid concentrations in pure POPC and ternary mixtures. <ref type="bibr">13</ref> C-lipid concentration comparisons between pure and ternary mixtures formed by gentle hydration (A) or electroformation (B). For the pure POPC sample, the <ref type="bibr">13</ref> C-labeled lipid is 13 C 18 -POPC, whereas for the ternary mixture, <ref type="bibr">13</ref> C 18 -DSPC is the <ref type="bibr">13</ref> C-labeled lipid. Dot plot representations of these distributions can be seen in Figure <ref type="figure">S27</ref>. (C) Standard deviations for the distributions shown in (A) and (B).</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head n="8.">Comparison of Pure POPC Mixtures with different Methods</head><p>Figure <ref type="figure">S10</ref>. 2 H 31 -POPC concentrations in pure POPC samples. 2 H 31 -POPC concentrations were compared between samples formed either by electroformation or gentle hydration. The resulting 2 H 31 -POPC concentrations distributions can be seen in (A) and the corresponding standard deviations can be seen in (B). Dot plot representations of these distributions can be seen in Figure <ref type="figure">S28</ref>.</p><p>The variability in pure POPC samples ( 13 C 18 -POPC: 2 H 31 -POPC:POPC 20:20:60) was examined for patches formed via electroformation and gentle hydration. As seen in Figure <ref type="figure">S10</ref>, the 2 H 31 -POPC concentration variability is not significantly different between the two methods as assessed via F-test. This suggest that the methods themselves do not have significant differences in terms of surface contamination, and that any difference in variability between gentle hydration and electroformation in ternary mixtures is not the result of one method being more prone to contamination of GUVs or systematic error in the NanoSIMS analysis. Instead, the differences in variability are due to the method of GUV preparation itself.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head n="9.">Detector Counts on Bilayers and Exposed Substrate</head><p>To further explore if significant contamination is present on the surfaces on which bilayers were deposited, the raw ion counts on several detectors were examined. 30x30 pixel (2.9x2.9 &#956;m) regions of interest (ROIs) were defined either on exposed substrate or on bilayer patches. For at least 23 ROIs on either bilayer or exposed substrate, the total counts on the detectors set to 2 H -, <ref type="bibr">13</ref> C -and 12 C -were determined and then averaged. This analysis was done for GUV patches containing 13 C 18 -DSPC and 2 H 31 -POPC. GUV patches formed via gentle hydration or electroformation were analyzed (data taken from the same patches discussed in Figures <ref type="figure">4</ref> and <ref type="figure">5</ref> of the main text), with the results displayed in Figure <ref type="figure">S11</ref>. For all detectors, average counts were significantly higher in regions containing bilayer than those without bilayer (assessed via t-test). This suggests that there is minimal contamination of the surface, particularly for 2 H -and 13 C -, which have 32 and 9 times more counts on bilayer than on exposed substrate respectively. 12 C -counts on a bilayer are only three times higher those on exposed substrate, suggesting some level of <ref type="bibr">12</ref> C is deposited on the surface even where there is not bilayer present, which is the major limitation of the absolute quantification method detailed in the main text. However, the 12 C -counts appear to be consistent from sample to sample, as the raw <ref type="bibr">12</ref> C -counts on exposed substrate are similar between the exposed gentle hydration and electroformation substrates. Additionally, as suggested by Figure <ref type="figure">6</ref>, Figure <ref type="figure">S9</ref>, and Figure <ref type="figure">S10</ref>, sample contamination alone cannot explain the large variability seen in labeled lipid concentrations of ternary samples. This suggests that while there is 12 C -background, this background is relatively constant from sample to sample and not a major factor driving the observed variability in ternary mixtures.</p><p>Additionally, it should be noted that although the natural abundance of 13 C -is 1.1%, the average 13 C -/( 13 C - + 12 C -) ratio on regions without bilayer is 0.039 and 0.024 for gentle hydration and electroformation respectively. This corresponds to 3.9% and 2.4% of the carbon background on the substrate being labeled. This suggests that some debris is deposited on the substrate, presumably during the rupturing and freezedrying process, leading to <ref type="bibr">13</ref> C enrichment on the substrate from <ref type="bibr">13</ref> C-labeled lipids.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>S13</head></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>Figure S11</head><p>. Comparison of detector counts on bilayer and exposed substrate. Raw counts on detectors set to 2 H -, <ref type="bibr">13</ref> C - and 12 C -were compared between regions of interest on bilayer and on exposed substrate. The top row displays counts on each detector for bilayer patches formed from GUVs via gentle hydration and regions on the same sample where the substrate is exposed (no bilayer present). The bottom row displays the equivalent data collected from an electroformation sample. Both samples had the same nominal lipid composition (DSPC: <ref type="bibr">13</ref>  </p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head n="10.">Correcting Ternary GUV Compositional Variability</head><p>In order to approximate the absolute compositional variability (excluding variability resulting from sample preparation and NanoSIMS analysis), the standard deviations of ternary GUVs formed via electroformation or gentle hydration were corrected using several different methods. The expression used to correct the ternary samples can be seen in figure <ref type="figure">S12A</ref>, where the observed ternary standard deviations are adjusted using the standard deviations from different control samples. The standard deviations from four different control samples were used as proxies for the variability due to sample preparation and instrumentation: continuous bilayers formed via rupture of 100nm diameter extruded ternary SUVs, monolayers composed of pure POPC with 20 mol % 2 H 31 -POPC and GUVs with the sample pure POPC composition formed via either gentle hydration or electroformation. The first sample can be used as a proxy for the variability in sample preparation and analysis for 2 H -/ 13 C -ratio, as well as <ref type="bibr">13</ref> C 18 -DSPC and 2 H 31 -POPC concentrations. Continuous bilayers formed from SUVs should be more compositionally homogenous than individual GUV patches, as each bilayer within a corral of a NanoSIMS substrate is S14 formed from the rupture of hundreds of SUVs, and therefore should represent the average composition more accurately. This can be seen in Figure <ref type="figure">S12B</ref>, where the variability in ternary SLBs formed from extruded SUVs is consistently lower than variability in ternary GUV patches formed via either gentle hydration or electroformation. Correcting the standard deviations of the ternary GUV samples using the ternary SUV bilayers yields the table shown in Figure <ref type="figure">S12C</ref>. However, since the preparation of continuous bilayers is not identical to that of bilayer patches, additional corrections were performed using pure POPC samples ( 2 H 31 -POPC: <ref type="bibr">13</ref> C 18 -POPC:POPC 20:20:60). Since these samples have only one lipid component, they are expected to be far more uniform. Additionally, pure POPC samples formed via electroformation and gentle hydration have the same sample preparation protocol as their corresponding ternary mixtures, and therefore are better at approximating the variability inherent to sample preparation. However, these samples can only be compared via their 2 H 31 -POPC concentrations, as they do not contain 13 C 18 -DSPC, and the corresponding 2 H -/ 13 C -ratio is not comparable between pure and ternary mixtures (Figure <ref type="figure">3</ref> of the main text). Corrections were also performed using pure POPC monolayers, which may be the best approximation for instrument variability, given that they should be essentially homogenous. However, since monolayers have distinctly different sample preparation from GUV patches, correcting ternary GUV patch variability with pure POPC monolayer variability does not account for variation introduced from sample preparation.</p><p>The corrections to 2 H 31 -POPC concentration variability can be seen in figure <ref type="figure">S12D</ref> and are similar regardless of which pure POPC sample is used. The results in Figures <ref type="figure">S12C</ref> and <ref type="figure">D</ref> can be taken as approximations of GUV variability due to composition and method, rather than variability from sample preparation and NanoSIMS analysis. It should also be noted that since the measured variability in 2 H 31 -POPC concentrations in continuous bilayers produced from ternary SUVs is nearly identical to that of the other pure POPC samples, the covariance between sample preparation, instrument noise and the observed variability is likely minimal (ternary mixtures do not have inherently higher variation due to sample preparation or instrumentation). </p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head n="11.">Nanoscale Heterogeneity in Lipid Films</head><p>Since GUVs are produced from a dried down lipid film, one possible source of GUV heterogeneity is the lipid film itself. Lipids films made from the ternary master stock containing <ref type="bibr">13</ref> C 18 -DSPC and 2 H 31 -POPC were examined for macroscale separation. As can be seen in Figure <ref type="figure">S13</ref>, there is no resolvable macroscopic separation between <ref type="bibr">13</ref> C 18 -DSPC and 2 H 31 -POPC within the ternary film. The NanoSIMS 50L has a 50nm lateral resolution limit, so it is possible that preferential interactions between lipids within the film take place over a smaller length scale. To address this possibility, recombination between <ref type="bibr">13</ref> C -and 2 H - to form triply labeled acetylide ( 13 C 2 2 H -) was examined. When a sample containing <ref type="bibr">13</ref> C-and 2 H-labels on different lipids is rastered by the primary Cs + , secondary <ref type="bibr">13</ref> C -and 2 H -ions are ejected and then recombine to form <ref type="bibr">13</ref> C 2 2 H -(figure <ref type="figure">S14A</ref> and <ref type="figure">B</ref>). This method, which has be examined in detail previously with <ref type="bibr">13</ref> C 15 N - or <ref type="bibr">13</ref> C 2 2 H -recombinant ions <ref type="bibr">[2]</ref><ref type="bibr">[3]</ref><ref type="bibr">[4]</ref> , is sensitive to the average distance between the two different isotopic labels on the order of 1-3nm. The recombination ratio, &#8476;( 13 C 2 2 H -) = <ref type="bibr">13</ref> C 2 2 H -/( 13 C 2 2 H -+ <ref type="bibr">13</ref> C 2 H -+ 12 C 2 2 H -) for a ternary film containing <ref type="bibr">13</ref> C 18 -DSPC and 2 H 31 -POPC can be compared to the &#8476;( 13 C 2 2 H -) for a ternary film containing <ref type="bibr">13</ref> C 18 -POPC and 2 H 31 -POPC. This comparison provides information on which pair of labeled lipids is closer together on average. A higher &#8476;( 13 C 2 2 H -) corresponds to a smaller average distance between the differently labeled lipids, while a lower &#8476;( 13 C 2 2 H -) corresponds to a larger average distance between the differently labeled lipids. As can be seen in Figure <ref type="figure">S14C</ref>, the measured &#8476;( 13 C 2 2 H -) in films containing <ref type="bibr">13</ref> C 18 -and 2 H 31 -POPC is higher than that of films containing 13 C 18 -DSPC and 2 H 31 -POPC (assessed via ttest). This suggests that, on the nanoscale, the two differently labeled POPC lipids are closer together on average than labeled DSPC and POPC. Although this may not be the source of the observed GUV heterogeneity, it provides a possible explanation for the correlations observed in Figure <ref type="figure">3</ref> of the main text.</p><p>POPC lipids interact strongly with each other relative to POPC and DSPC, which is presumably why <ref type="bibr">13</ref>    </p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head n="12.">Relation between GUV Patch Size and 13 C 18 -DSPC and 2 H 31 -POPC Concentrations</head><p>In order to assess other possible sources of GUV variability, the relationship between patch area and the measured concentrations of 2 H 31 -POPC and 13 C 18 -DSPC was assessed. Ternary GUV patches that could fit entirely into the 25x25 &#956;m 2 raster were analyzed and their 2 H 31 -POPC and 13 C 18 -DSPC concentration determined via external calibration curves. This method was applied to ten GUV patches formed via gentle hydration (Figure <ref type="figure">S15 A-B</ref>) or electroformation (Figure <ref type="figure">S15 C-D</ref>) with no correlation seen between patch area and labeled lipid concentrations. This suggests that there is minimal effect of patch area on the resulting labeled lipid concentrations. However, it should be noted that this analysis is restricted to GUVs that can fit within the analysis region and may not apply to larger GUVs. Additionally, since the S18 history of each GUV patch is not known, the area of a given patch may not be a direct reflection of the size of GUV that produced the observed patch. </p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head n="13.">Assessing the Extent of GUV Mixing during Patch Formation</head><p>Although it is generally assumed that each GUV patch results from the rupture and spreading of a single GUV, it is possible that two or more GUVs may rupture near each other and combine, which would average their two concentrations together. This would result in a decrease in the observed variability in labeled lipid concentrations. To assess the extent to which GUV patches can combine, ternary GUVs with the compositions DSPC:POPC: 2 H 31 -POPC:CHOL 40:20:20:20 ( 2 H 31 -POPC GUVs) and DSPC: <ref type="bibr">13</ref>  The 2 H -/ 13 C -ratio was then measured for 15 ( 2 H 31 -POPC GUVs and 13 C 18 -DSPC GUVs) or 30 ( 2 H 31 -POPC GUVs + 13 C 18 -DSPC GUVs) patches. The resulting ratios can be seen in Figure <ref type="figure">S16</ref>. Both the 2 H 31 -POPC and 13 C 18 -DSPC GUV patches have high and low ratios respectively and are monomodal, as only one labeled lipid is present in each of these GUV samples. For the sample prepared using a mixture of 2 H 31 -POPC GUVs and 13 C 18 -DSPC GUVs ( 2 H 31 -POPC GUVs + 13 C 18 -DSPC GUVs), the resulting ratios have two major clusters, or populations. Each population has a ratio similar to that of the samples prepared using only 2 H 31 -POPC or <ref type="bibr">13</ref> C 18 -DSPC GUVs, suggesting that the majority of the patches are not substantially cross-contaminated with another labeled species. This suggests that the majority of individual GUV patches do not result from the averaging of multiple GUVs. However, a few patches were observed with some level of cross-contamination and have 2 H -/ 13 C -ratios between the two major populations. These patches contain approximately 3-5 mol % of the non-dominant labeled species, suggesting that some averaging between the 13 C 18 -DSPC and 2 H 31 -POPC GUVs occurred in these patches. Additionally, the two dominant populations in the mixed sample appear to be slightly broadened, suggesting that these patches may have some minimal level of contamination from the other label. This may result from debris formed during the freeze-drying process (see section 9 of the supporting information) which may lead to minor contamination (less than 1 mol %) of patches. This merging of GUVs to form single patches has been observed in prior work <ref type="bibr">5</ref> .</p><p>It should also be noted that one GUV may lead to the formation of several patches 6 . It cannot be established in this work if different patches resulting from the same parent GUV have different compositions, as these patches cannot be distinguished. As a result, it cannot be determined how this affects the observed variability in GUV patch composition. Given that one GUV may lead to multiple patches, this makes the use of a mixture with nanoscale separation, such as the mixture used here, a better choice for compositional analysis, as GUVs with nanoscale separation are likely less susceptible to large-scale compositional changes during GUV rupture. S20 Figure S16. Extent of mixing during GUV patch formation. 2 H -/ 13 C -Ratios for GUV patches formed with GUVs containing either 2 H 31 -POPC or 13 C 18 -DSPC.When GUVs containing only 2 H 31 -POPC or only <ref type="bibr">13</ref> C 18 -DSPC were deposited on NanoSIMS substrates, high and low 2 H -/ 13 C -ratios were measured respectively. When a mixture of the differently labeled GUVs were deposited on a NanoSIMS substrate, the resulting distribution of ratios is predominantly bimodal, with ratios reflecting that of the GUV patches containing only one labeled lipid.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head n="14.">Tracking Cholesterol Concentrations in Triply Labeled Bilayers</head><p>As further evidence for the changes in cholesterol concentration seen in Figure <ref type="figure">7</ref>, triply labeled bilayer patches were formed from GUVs using a master stock with nominal composition 20:20:20:20:20 POPC: <ref type="bibr">15</ref> N-POPC:DSPC: 13 C 18 -DSPC: 2 H 7 -CHOL. This mixture allows for every lipid component to be tracked. Calibration curves for 15 N-POPC and 2 H 7 -CHOL are shown in Figure <ref type="figure">S17</ref>. As can be seen in Figure <ref type="figure">S18A</ref>, each isotopically labeled species can be visualized and quantified. The average concentration across 8-9 bilayers was measured for patches formed via electroformation or gentle hydration as well as continuous bilayers formed via rupture of SUVs. Figure <ref type="figure">S18B</ref> shows that neither 13 C 18 -DSPC nor <ref type="bibr">15</ref> N-POPC are significantly different in terms of concentration between the three methods, but the average 2 H 7 -cholesterol concentration is significantly lower in electroformed GUVs. The 15 N-POPC concentration is also significantly higher in GUVs formed by electroformation relative to those formed by gentle hydration (p = 0.037). Since electroformed GUVs contain less cholesterol, the concentration of other lipids must be higher, explaining why there is a slight increase in the average 15 N-POPC concentration. However, it should be noted that a total of 9 average concentrations were compared (three different labeled species compared between the three different methods), therefore if the required significance level is adjusted to match the number of comparisons using the Bonferroni correction, (&#945;/n, where n = the number of tests conducted), S21 this difference would no longer be statistically significant. The changes in cholesterol concentration remain statistically significant regardless of whether or not this adjustment is performed.  S22 15. Assessing Average 2 H 31 -POPC Concentration in 13 C 27 -CHOL Containing Patches</p><p>The average concentration of 2 H 31 -POPC was also measured for the GUV patches analyzed in Figure <ref type="figure">7</ref> of the main text. Analysis was conducted both including and excluding edges (same regions of interest as were used in Figure <ref type="figure">7</ref> of the main text). Figure <ref type="figure">S19</ref> shows that regardless of whether edges were included or excluded, the average 2 H 31 -POPC remained constant between GUV patches formed via electroformation or gentle hydration. It may be expected that the 2 H 31 -POPC concentration in electroformed GUVs will be higher, given that the lower average cholesterol concentration in these GUVs must lead to higher concentrations of other species; however, there is only a 3-5 mol % decrease in cholesterol concentration in GUVs formed via electroformation. Labeled 2 H 31 -POPC only makes up 25% of the remaining mixture so only a 1 mol % increase could reasonably be expected in 2 H 31 -POPC concentration in electroformed GUVs relative to GUVs formed via gentle hydration. Since this is a relatively small change, it is not detectable within the error of the experiment. Additionally, there is a statistically significant increase (p = 0.0489) in 2 H 31 -POPC between analyses including and excluding the edges of bilayer patches formed by gentle hydration (Figure <ref type="figure">S19</ref>). Since the analyses with the edge excluded show a higher average 2 H 31 -POPC concentration, this can be rationalized by noting that cholesterol appears to be enriched at the edges of bilayers patches, therefore another lipid must have a lower concentration. However, this increase in 2 H 31 -POPC concentration does not appear to be consistent for both GUVs formed via gentle hydration or electroformation and is not reproduced in replicates (Figure <ref type="figure">S20B</ref>). This suggests that the change in POPC concentrations between the edge and center is minimal within the sensitivity of the experiment.</p><p>Given that other GUV patch samples display disagreements between their measured concentrations and the nominal concentrations of their respective stock solutions, these replicates suggest that these differences in cholesterol concentration are reproducible and likely not the result of error in the measurements. However, it is worth noting that there is error in absolute concentration for cholesterol as well, likely due to aforementioned surface contamination, calibration curve accuracy and error in the concentration of lipid stock solutions.  </p><note type="other">S23</note></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head n="16.">Tabulated Raw GUV Composition Data</head><p>The tables below record 2 H -/ 13 C -ratios, as well as labeled lipid compositions for all figures in the main text and supporting information. Ratios and compositions are reported to two decimal places, although it should be noted that the second decimal place is presumably less accurate, however for the purposes of reporting the data, it has been left in. Tables have been organized according to the figures they correspond to (as denoted in the table title) and are listed in order of appearance.</p></div><note xmlns="http://www.tei-c.org/ns/1.0" place="foot" xml:id="foot_0"><p>https://doi.org/10.1021/jacs.3c09039 J. Am. Chem. Soc. 2023, 145, 27521-27530</p></note>
			<note xmlns="http://www.tei-c.org/ns/1.0" place="foot" xml:id="foot_1"><p>Journal of the American Chemical Society</p></note>
		</body>
		</text>
</TEI>
