Open Access Article
Shuo Qian
*a,
Honghai Zhanga,
Wellington Leitea,
Andrew E. Whitten
b,
Piotr Zolnierczuk
a,
Yue Yuana and
Qiu Zhanga
aOak Ridge National Laboratory, Oak Ridge, TN, USA. E-mail: qians@ornl.gov; Tel: +1 865 241 1934
bAustralian Nuclear Science and Technology Organisation, Lucas Heights, NSW 2234, Australia
First published on 12th March 2026
Plastic waste is now pervasive in the environment, breaking down into microplastics and nanoplastics under many environmental conditions. These particles have been found in various ecosystems and even in human tissues, raising significant environmental and health concerns. In this study, we investigated the interaction of polystyrene nanoplastics, with and without surface modifications, on biomembrane structures using contrast-matching small-angle neutron scattering and neutron spin echo spectroscopy. The neutron contrast matching enabled the selective study of biomembranes in the presence of nanoplastics. Two model membranes were employed: a simple zwitterionic bilayer (i.e., dimyristoylphosphatidylcholine [DMPC]) and an Escherichia coli lipid extract as a bacterial membrane model. The results show profound membrane disruptions, including possible thinning, vesicle fragmentation, lipid monolayer formation, and inter-vesicle aggregation, with the more severe effects observed in DMPC membranes. Notably, E. coli membranes exhibited greater resilience, suggesting that natural membranes with diverse lipid compositions may reduce susceptibility to perturbation by extracellular nanoplastics. These findings highlight potential risks posed by environmental nanoplastic particles to biological membranes, with insights into molecular-level interactions and the environmental toxicity of nanoplastics. This work provides a foundation for future studies into nanoplastic–biomembrane interactions and their broader implications for health and environment using neutrons.
Environmental significanceNanoplastics represent a pervasive and emerging environmental threat with the potential to disrupt critical biological processes. This study provides a molecular-level understanding of how polystyrene nanoplastics interact with biomembranes, revealing significant structural disruptions, including possible membrane thinning, vesicle fragmentation, aggregate formation, and lipid monolayer formation. By employing advanced neutron scattering techniques, this research bridges and narrows a critical knowledge gap in the nanoscale behavior of nanoplastics with biomembranes. These findings highlight the potential for nanoplastics to impair cellular integrity, amplifying their environmental toxicity, while also suggesting that natural membranes with diverse lipids may be less susceptible to such perturbations. These insights are vital for assessing the ecological and health risks of nanoplastics and guiding mitigation strategies to address their impact. |
Recent studies highlight the pervasive nature of microplastics and nanoplastics: they have been found in various ecosystems and even in human tissues. For instance, researchers have detected nanoplastics in human placentas, raising concerns about their potential influence on fetal development.7 Additionally, nanoplastics have been observed crossing the blood–brain barrier in animal models, suggesting possible neurological effects.8 Nanoplastics promote in vitro aggregation of α-synuclein, which is associated with dementias such as Parkinson's disease.9 In vitro studies have shown that micro- and nanosized plastic particles can accumulate in the tissues of various living organisms, including mussels, shrimp, and fish.10–12 These results all point to the presence of plastic particles in the biosphere as a potentially significant environmental and health issue, motivating the urgent need to understand the interaction of plastic particles at different biological levels.13
In this study, we focus on how nanoplastic particles (NPs) with different surface modifications interact and modify biomembrane molecular structures under various conditions using contrast-matching small-angle neutron scattering (SANS), along with additional insights into membrane dynamics from neutron spin echo (NSE) spectroscopy. NPs have been found to have greater potential for accumulation and translocation within host tissues11,14,15 than their micro-sized counterpart. Many studies to date have shown that particle size is a primary factor for plastic particles in entering and distributing within tissues and cells. The significantly increased surface area of smaller plastic particles can also adsorb and transport organic pollutants,16 transforming them into vehicles for toxic substances within host organisms. They can act as Trojan horses, facilitating the uptake of other environmental contaminants such as heavy metals and persistent organic pollutants.17–20 Biomembranes are important semipermeable barriers that separate cells or organelles from their surrounding environment and are crucial interfaces for many biological processes. Recent evidence suggests that plastic particles can enter cells via various pathways, including endocytosis and direct membrane penetration.21–25 Some studies have also shown the destabilizing effect on lipid membranes;26–31 however, direct experimental insights into how NPs alter membrane structure during their interaction are currently lacking.
SANS is a bulk technique capable of providing hierarchical structure information across different length scales from a few nanometers to a few hundred nanometers simultaneously. It is a key method for probing from trans-bilayer membrane structure to overall vesicle morphology in solutions.32–37 Compared to other structural techniques such as X-ray scattering, electron and light microscopy, the large difference in neutron scattering cross-section between hydrogen isotopes (protium and deuterium) allows adjustment of contrast between solute and solvent by sample deuteration and/or changing solvent D2O/H2O ratios.38,39 This allows contrast matching, by which scattering from the NP structure is eliminated, enabling selective observation of membrane structure. Neutron contrast matching has proven highly useful in elucidating complex, multicomponent systems.32,40 At the same time, biomembranes have shown extraordinary viscoelastic properties that are required for essential biological functions. NSE spectroscopy measures the neutron energy change in the orders of a few nano-electronvolts (neVs) as it scatters through a sample. The energy change reveals sub-microsecond dynamic motions on length scales from a few to hundreds of nanometers, allowing the derivation of various elastic-mechanical properties, such as the membrane bending modulus and bilayer area compressibility.41–44 These are important membrane characteristics that occur under stress or strain during interactions with other molecules. Changes in membrane rigidity can influence key biological functions and may even compromise the integrity of the membrane as a whole. These changes also corroborate the membrane structure changes revealed by SANS.
Two model membrane systems were selected for this study: 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) as a simple zwitterionic model bilayer to represent charge neutral membranes, and Escherichia coli (E. coli) total lipid extract as a model for bacterial membranes with anionic lipids that are ubiquitous in the environment.45 Notably, fully deuterated E. coli cells were grown, and their total lipid was extracted and used to enhance neutron contrast. Both lipids were used to form large unilamellar vesicles (LUVs) with relatively narrow size distribution, widely used as model systems in biophysical research.41,46,47 We caution readers that the model systems used here represent much simplified biomembranes compared to those of live cells. For instance, lipopolysaccharide and membrane proteins are absent in these models. Their potential influence on the extent and outcome of the interaction was not considered in this study.
We recognize that micro- and nanoplastics in the environment represent highly heterogeneous and dynamic systems. Their chemical composition, surface functionalization, degradation pathways, and environmental conditions—including pH, ionic strength, and interactions with natural organic matter—span a vast parameter space that is difficult to capture in a single laboratory study.48,49 Environmental exposure scenarios are also highly variable; for example, nanoparticles that are inhaled may experience very different conditions than those that are ingested, and factors such as salt composition and natural polymers can strongly influence particle–membrane interactions. In this context, commercially available spherical polystyrene (PS) NPs of uniform size were selected as model particles for this study. Although PS accounts for only a small fraction of global polymer production,50 spherical NPs of uniform size are produced commercially with a variety of surface functionalization, such as unmodified, or aminated (AM) (positively charged) and carboxylated (CA) (negatively charged). Several studies have therefore used them to provide valuable insights into the environmental impact of NPs even though they may not fully capture the wider heterogeneity of environmental NPs.9,14,15,27 This work demonstrates how neutron scattering techniques can provide mechanistic insight into nanoplastic–membrane interactions under solution conditions, highlighting their potential to complement environmental and toxicological studies.
To produce deuterated E. coli total lipid, E. coli BL21 (DE3) strain was inoculated into 5 mL of Luria–Bertani media from glycerol stock and incubated at 37 °C overnight. The cell cultures were then adapted to grow in D2O Enfors-based minimal media with D-glycerol as the carbon source. The cells underwent three exchanges of media to ensure complete adaptation from H2O to D2O-based media. Then the culture was grown at 37 °C overnight in the D2O Enfors minimal media51 supplemented with D-glycerol with initial inoculation at an optical density at 600 nm of approximately 0.1. The cells were harvested by centrifugation at 6000×g for 30 min and then washed with 1% NaCl three times followed by freeze-drying for 48 h before lipid extraction. The E. coli total lipid was extracted using the modified Bligh and Dyer method.52 In brief, 250 mL chloroform, 500 mL methanol, and 200 mL water were sequentially added to the freeze-dried cells (1.1 g) to achieve a final chloroform
:
methanol
:
water ratio of 1
:
2
:
0.8 (v/v/v). The cells were vortexed for 15 s immediately following the addition of each solvent and allowed to stand for approximately 18 h at room temperature with stirring. This process was followed by the addition of chloroform (250 mL) and water (250 mL) to obtain a final chloroform
:
methanol
:
water ratio of 1
:
1
:
0.9 (v/v/v) to separate the aqueous and organic fractions of the mixture. A separatory funnel was used to collect the lower chloroform phase, which contains the lipid extract. The chloroform–methanol mixture was then evaporated under reduced pressure using a rotovap (Buchi, R-205). The lipid extract was dried under vacuum overnight to remove traces of organic solvents.
At Bio-SANS, the neutron wavelength was set to 0.6 nm with a wavelength spread of Δλ/λ = 15%. The collimation distance was 17.42 m. Two detector configurations were used, with sample-to-detector distances of 1.1 m (wing detector) and 15.5 m (main detector). A neutron beam with a diameter of 14 mm was used to illuminate the samples.
At Bilby, the time-of-flight mode was employed with a wavelength range of 0.2–1.0 nm and a wavelength resolution of Δλ/λ = 10.8%. The collimation distance was 10.75 m. The sample-to-detector distances were 11.0 m (main detector), 2.50 m (upper curtain), 2.99 m (lower curtain), 1.50 m (left curtain), and 1.50 m (right curtain). A neutron beam with a diameter of 12.5 mm was used to illuminate the samples.
The lipid solutions were transferred into quartz circular cells with a 1 mm beam path (Hellma, Germany) for measurements. The samples were kept at 30 °C during the experiment. The acquired data were processed using the facility-provided data reduction software (drt-SANS version 1.10.2 for Bio-SANS; Mantid version 6.8.0 for Bilby). The processed data were corrected for transmission and background scattering from solvent and quartz cells, and normalized to an absolute intensity scale relative to the primary beam, resulting in 1D scattering curves (I vs. Q), where I is the scattering intensity arising from the contrast between the lipids and the solvent, and Q is the scattering vector from Q = 4π
sin
θ/λ, where θ is half of the scattering angle, and λ is the neutron wavelength. The instrument resolution was provided with the data and was used in the SANS data analysis to account for the smearing effect. This ensures that the impact of the differences in instrument configurations from different instruments was mitigated.
To determine the contrast matching point (CMP) of NPs, the scattering intensity I at Qmin = 0.03 nm−1 was measured at different D2O ratios (v/v). From a linear regression of
vs. D2O ratio (v/v), the CMP, in terms of D2O ratio, was determined from where the line intercepts zero intensity.38 Additional details are included in the SI.
The trans-bilayer membrane structure, dominated by the hydrocarbon chain high in deuterium, was modeled as a single slab of neutron scattering length density (SLD), defined as the total coherent neutron scattering length of a molecule divided by the molecular volume. It was fitted with an isotropic lamellar model. The scattering intensity I(Q) for the randomly oriented bilayer in a vesicle is I(Q) = scale × (2πP(Q)/Q2d) + background, where the form factor is P(Q) = (2Δρ2/Q2) × (1 − cos(Qd)). The total layer thickness is given by d, and Δρ is the SLD difference (the contrast) between the lipid and the solvent.55 The polydispersity of the thickness d is assumed to be a Gaussian distribution implemented in SASView. In the case of the inter-bilayer peak presented in the data, a Caillè structure factor is used as implemented in SASView to determine the thickness. Note the uncertainties reported in the tables reflect fitting uncertainty rather than an assertion of experimental resolution in the structure.
A generalized Guinier–Porod function is used to fit low-Q data when applicable. This empirical model provides the size and dimensionality of the sample, including shape intermediates between typical well-defined 1D, 2D, or 3D shapes (e.g., long rod, thin/extended sheet, or globular particle, respectively).56 In the case of a possible globular 3D shape within the Q range, the Guinier function is used to obtain the radius of gyration Rg by fitting I(Q) = scale × exp{(−Q2Rg2/3)} + background.57 The fitting parameters and details are shown in the SI.
The data analysis follows Zilman and Granek,60 using a freestanding single-membrane fluctuation model to fit the intermediate scattering function
, where t is the Fourier time. Then the Q-dependent relaxation rate was obtained from fitting Γ(Q). The effective bending modulus
can be obtained using the equation from Watson and Brown:61
, where γ is assumed to be 1 when
≫ kBT, which is typical for lipid bilayers near room temperature,44,60,62 and ηD2O is the viscosity of D2O at 25 °C.
| Diameter (nm) | PDI | ζ-potential (mV) | |
|---|---|---|---|
| PS | 35.3 | 0.09 | −36.37 |
| CA | 42.1 | 0.04 | −31.75 |
| AM | 52.6 | 0.03 | 27.83 |
From the contrast variation series measured at different D2O ratios, we determined that the CMP for NPs is 26 ± 2% D2O (Fig. 1), consistent with polystyrene SLD.64 Since polystyrene is the major constituent of all PS NPs, only AM PS was examined in the CMP measurements. Limited SANS beam time led us to select the largest particles to maximize scattering intensity. SANS and scanning electron microscopy confirmed the spherical morphology of the particle (see SI for details). All the SANS samples with LUVs were prepared in 26% D2O to eliminate the scattering from NPs and highlight deuterated lipid D54-DMPC and E. coli lipid with CMP of approximately 95% D2O (dominated by the D54 chain) and higher than 100% D2O,65 respectively.
![]() | ||
| Fig. 1 Contrast match point of the PS nanoparticle. The intersection of the linear fitting of the data points to the zero-intensity line indicates that the PS CMP is 26 ± 2% D2O. | ||
The high-Q region shows the more detailed structure across the lipid bilayer. The DMPC-only scattering profile displays a sharp minimum near Q = 2 nm−1, indicating a well-defined bilayer–water interface and a slab-like deuterated hydrocarbon chain within the bilayer. Upon addition of NPs, all minima are reduced to blunt kinks. The much shallower slope approaching the kinks further shows that the bilayer–water interface is blurrier because the Porod exponent deviates to much less than 4.66 Based on previous experience with the SANS data analysis, the blunting of the sharp minimum resulting from a bilayer indicates either aggravating asymmetry in the lipid composition across the bilayer or increasing thickness heterogeneity.35,39,47 In this case, only DMPC is presented so we attributed the change to thickness heterogeneity, for example, the polydispersity of thickness. Previous studies have shown that plastic particles can induce mechanical stretching of membranes, resulting in thinning of the bilayer.26
The effects of the different NP types and concentrations on overall vesicle morphology also vary. For AM, at a very high concentration (2%), the scattering curve starts to flatten at low Q, suggesting the formation of particles smaller than the original LUVs, possibly smaller vesicles or heterogeneous aggregates (Fig. 2A). Fitting these data with the Guinier–Porod function yields an Rg of approximately 14.7 nm, much smaller than the extruded LUV radius of about 50.0 nm. By contrast, the Porod exponent is close to 2, indicating that the large thin sheet remains present. So, the features of smaller particles and thin sheet bilayers coexist in the sample. The possible smaller vesicles or aggregates are not multilamellar as there is no visible multilamellar diffraction peak as in some of E. coli lipid (Fig. 5A).
The model-free assessment of the data points to a wide range of resulting structures from NP–membrane interaction. Given the variety of features observed in the scattering profiles, particularly in the low-Q region, no single model can adequately describe the entire Q range in the SANS data for all the samples. We focused on the high-Q data (>0.4 nm−1), which encoded the thickness heterogeneity, to identify a common trend among the different NPs at different concentrations. Typically, a core-multi-shell model can be used to describe a detailed lipid bilayer structure accounting for the overall size of the vesicle, thicknesses of headgroups, and hydrocarbon chains.35,36,67 However, the overparameterization and inter-dependency of these parameters cause overfitting and ambiguity in interpretation. Hence, the more general and less ambiguous lamellar model was chosen to fit the high-Q data.55 The model describes a single thickness in a lamellar structure distributed in a random orientation with a polydispersity factor, which is used to gauge heterogeneity in thickness.
The DMPC-only sample shows a thickness of 2.8 nm with negligible thickness polydispersity (Table 2). Adding AM significantly reduces the bilayer thickness and significantly increases the thickness polydispersity (Fig. 3A). For example, 0.5% of AM thins the bilayer to 2.5 nm with about 19% polydispersity. Then 1% AM further thins the bilayer and increases polydispersity further. The trend partially reverses at 2% AM with slightly increased thickness and decreased polydispersity and then at 1% AM, coincident with the possible formation of smaller particles such as smaller vesicles or aggregates. If the thinning effect is caused by the stretching of membrane in the presence of AM NPs,26,68 then the reduction may indicate that the tension caused by AM particles is partially relieved as those smaller particles start to form (Fig. 3B). The dimension variable obtained from the Guinier–Porod analysis is 0.38, corresponding to a dimensionality parameter 3 − s = 2.62.56 This intermediate value between a two-dimensional membrane sheet and a three-dimensional particle suggests that the vesicle membranes have deviated from smooth bilayers toward partially space-filling, disordered mass-fractal structures. Such behavior is consistent with membrane corrugation, folding or nanoparticle-induced aggregation or cluster formation on or within the bilayer. The reduced Rg further supports the formation of smaller vesicles or aggregates. For the disrupted LUVs, this structural evolution may produce heterogeneous, partially space-filling assemblies that exist as nanoparticle–lipid aggregate complexes or smaller vesicles, or other yet-unidentified structures. Meanwhile, the Porod exponent remaining close to 2 over a broad Q range indicates that a substantial fraction of lipids are still in the form of vesicular bilayers.
| Thickness (nm) | Thickness polydispersity (%) | Note | |
|---|---|---|---|
| DMPC lipid only | 2.8 ± 0.1 | Negligible | |
| AM 0.5% | 2.5 ± 0.2 | 20 ± 2 | |
| AM 1% | 2.2 ± 0.1 | 30 ± 5 | |
| AM 2% | 2.5 ± 0.1 | 15 ± 4 | Low Q is fitted with the Guinier–Porod function |
| CA 0.5% | 2.4 ± 0.1 | 20 ± 1 | |
| CA 1% | 2.2 ± 0.1 | 29 ± 3 | |
| CA 2% | 1.8 ± 0.1 | 52 ± 8 | |
| PS 0.5% | 2.9 ± 0.1 | 32 ± 1 | Low Q is fitted with the Guinier function |
| PS 1% | 2.2 ± 0.1 | 34 ± 5 | |
| PS 2% | 1.5 ± 0.2 | 76 ± 14 | |
CA also induces thickness changes, although less drastically than AM. As shown in Fig. 2B, the power law exponents of 2, still indicative of a thin sheet, are the same as the DMPC-only sample at all concentrations. Smaller particles do not form as no Guinier region is shown within the Q range observed. The lamellar model fitting results show a generally decreasing bilayer thickness and increasing polydispersity as the CA concentration increases (Fig. 3A). The thickness changes to 1.8 nm, and polydispersity is greater than 50%, indicating a significant disturbance of CA on the membrane.
The unmodified PS NPs exhibited the most significant change (Fig. 2C). At only 0.5% concentration, much smaller particles are formed, as indicated by the leveling scattering profile at low Q. This can be fitted with the Guinier function. The result indicates that the Rg of the smaller particles is about 10 nm. Then, at higher concentrations of 1% and 2%, the scattering profiles revert to a power law decay with an exponent of about 2. However, the identical undulation near 0.1 nm−1 (highlighted by arrows in Fig. 2C) from both 1% and 2% samples indicates a mixture of smaller particles like the ones presented in the 0.5% PS sample. Fitting at the high-Q region reveals a significant decrease in lamellar thickness at 1% and 2% PS. The thickness at 2% PS is almost halved compared with the DMPC-only sample, accompanied by greatly increased polydispersity to more than 75%. This result suggests great disruption of the bilayer by the PS NPs to form a monolayer, likely on the surface of hydrophobic PS NPs (Fig. 3C).30 Given the large polydispersity, LUVs with bilayers, heterogenous aggregates and smaller PS-monolayer particles likely coexist in the sample. However, fitting with the Guinier function for accurate particle sizes combined with the power law cannot yield a satisfactory result, indicating a more complex mixture of different structures in the sample. Overall, large thin sheets remain, as the power law exponent of 2 suggests, but some of the sheets are monolayer instead of bilayer. We reiterate here that the drastic change and the complexity as shown in this case prevented us from applying a single model for fitting the data in the whole Q range. In some cases, as mentioned above, the samples become mixtures of vesicular and non-vesicular particles, so a typical vesicle model is inappropriate for extracting the thickness information. The more general lamellar thickness with polydispersity avoids these issues and provides consistent interpretation, at least for the high-Q region.
This significant effect can also be seen from NSE results (Fig. 4). Adding 1% AM or CA to the DMPC LUVs changes the effective bending modulus of the DMPC bilayer about 20% to 30%, either rigidifying or softening the bilayer, respectively. However, a large decrease in membrane rigidity occurs with 1% PS NPs added in DMPC. The effective bending changes from 203 ± 8 kBT to 24 ± 2 kBT, almost an order of magnitude smaller. This result suggests some fundamental changes to the membrane, corroborating the likely formation of monolayers indicated by the drastic thickness reduction. The intrinsic bending modulus and the bilayer area compressibility modulus can be derived from the effective bending modulus. However, these analyses are based on lipid bilayer models following Zilman–Granek theory, rather than on monolayers, for which no corresponding model has been derived to our knowledge.60 The detailed analysis is presented in the SI.
The greater change is likely driven by the hydrophobic nature of PS, which favors interaction with the hydrophobic core of the bilayer. This occurs even in the presence of initiator anions surrounding the PS stock solution, as indicated by the negative ζ-potential (Table 1). We hypothesize that these anions from water-soluble initiators are more easily dispersed during interactions with the bilayer headgroups, compared with the covalently bound functional groups present on AM and CA NPs. Another possible contributing factor is the actual NP-to-vesicle ratio. For the 0.5% w/v samples, we estimate that the ratio of AM to DMPC vesicles is approximately 1
:
2. Based on the smaller sizes of CA and PS measured by DLS, the corresponding ratios are roughly 1
:
1 and 3
:
2, respectively. The higher relative concentration of PS NPs may exacerbate their interaction with the membranes and contribute to the observed differences here. The pronounced effects observed with PS NPs highlight their potential to influence biological membranes. The larger, charged and more hydrophilic nature of AM and CA appears to mitigate some of these interactions compared with the unmodified PS NPs. But the resulting perturbations—such as membrane thinning (Fig. 3A) and the formation of smaller particles (Fig. 3B)—suggest that all three types of nanoplastics can cause significant changes to the membrane structure.
The LUVs formed from the E. coli lipid also exhibit a well-defined unilamellar bilayer, as evidenced by the low-Q data indicating an extended bilayer thin sheet structure (Fig. 5). In contrast to the DMPC bilayer, the first minimum between 1 and 2 nm−1 is less prominent, indicating that some degree of thickness polydispersity originated from its complex composition with lipids of different chain length. Fitting this region, the lamellar model shows a much thicker bilayer of 4.1 nm with about 18% thickness polydispersity (Table 3). This result is consistent with its higher degree of deuteration—including the headgroups—and the presence of both long- and short-chain lipids, which together contribute to a thicker lamellar structure.65 Interestingly, the addition of NPs led to sharper minima around the region, suggesting a likely decrease in thickness polydispersity; therefore, the interaction yields a more rigid and uniform membrane as NPs round up different lipids. At the same time, the addition of NPs reduces the bilayer thickness for all the samples, as evidenced by the shifts in the positions of the first minima to higher Q values. This effect was observed even at the lowest NP concentrations tested, demonstrating the significant influence of both unmodified and surface-modified NPs.
| Thickness (nm) | Thickness polydispersity (%) | Note | |
|---|---|---|---|
| E. coli lipid only | 4.1 ± 0.1 | 18 ± 1 | |
| AM 0.5% | 3.4 ± 0.1 | Negligible | Caillè structure factor is used to combine with the lamellar form factor |
| AM 1% | 3.3 ± 0.1 | Negligible | |
| AM 2% | 3.7 ± 0.1 | Negligible | |
| CA 0.5% | 4.1 ± 0.1 | 15.9 ± 1.1 | |
| CA 1% | 3.8 ± 0.1 | Negligible | |
| CA 2% | 3.6 ± 0.1 | Negligible | |
| PS 0.5% | 3.6 ± 0.1 | Negligible | |
| PS 1% | 3.5 ± 0.1 | Negligible | |
| PS 2% | 3.4 ± 0.1 | Negligible | |
The lamellar model fitting results also substantiate this thinning effect. For example, adding AM reduces the bilayer thickness from 4.1 nm to between about 3.2 and 3.7 nm, with negligible polydispersity factors to fit the data. Adding CA and PS NPs reduces the bilayer thickness gradually in a concentration dependent manner. No halved reduction in thickness occurs for any of the samples, therefore no monolayer formation, unlike in some of DMPC samples. This result indicates that the bilayer structure is preserved.
Although most NP additions preserve the unilamellar structure, the addition of 0.5% AM induces inter-vesicle attraction, evidenced by a prominent peak at Q = 1 nm−1 (corresponding to a repeating spacing of 6 nm as used by the Caillè structure factor). This is possible for positively charged AM NPs to bring the negatively charged E. coli bilayer in close proximity (Fig. 3D). Even in the 1% AM sample, a small shoulder peak can be seen near Q = 0.8 nm−1 on top of the lamellar form factor (indicated by the arrow in Fig. 5A). However, increasing AM concentration may balance out the overall charge distribution around the E. coli lipid LUVs, so the inter-vesicle aggregation is greatly reduced in the 1% AM sample and eventually disappears in the 2% AM sample. There might be other possibilities of aggregates, but regardless, the integrity of the bilayer is intact judging from the thickness change and higher concentration data.
At low and medium Q, the changes in the E. coli lipid bilayer are also less drastic than in DMPC. The extended thin sheet with the power law exponent of 2 occurs in all the samples. No smaller particles form, even with PS NPs, unlike DMPC. The unilamellar structure of the vesicles is well preserved. The complex lipid composition of the E. coli membrane, including charged lipids such as PG and lipids with high intrinsic curvature such as PE, as well as the diversity in chain lengths and saturation, appears to partially mitigate the influence of NPs compared with the simpler DMPC system. However, the membrane thinning in all samples even at low NP concentrations, and possible inter-vesicle aggregation in some cases, still significantly affect the membranes.
The membrane systems employed include a basic phospholipid system with a zwitterionic DMPC lipid and a more compositionally complex membrane with an anionic E. coli lipid, enabling comparison of PS NP interactions under different conditions. The data reveal significant disturbances in both membrane types at the molecular level, regardless of NP concentration, charge, or surface modification. A few different possible modes of membrane disruption were identified at the molecular level.
The most prevalent effect is membrane thinning (Fig. 3A). Although overall vesicle shape remained unchanged, as observed in several DMPC and E. coli lipid samples, membrane surface adsorption generated tension that thinned the bilayer considerably.68 This thinning is likely to affect membrane rigidity, potentially influencing various membrane-associated processes (e.g., membrane raft and membrane proteins). The second possible mode of disruption is vesicle fragmentation into smaller vesicles and/or other heterogeneous particles. This mode was pronounced at high NP concentrations, particularly in AM-DMPC samples (Fig. 3B). The third possible mode involves lipid monolayer formation induced by interactions with hydrophobic PS NPs (Fig. 3C). Fragmentation and monolayer formation are both destructive to vesicle and bilayer integrity.
Interestingly, the more biologically relevant E. coli lipid membrane exhibited some degree of resilience to these destructive effects, showing fewer signs of damage. This result suggests that membranes composed of a diverse range of lipids—especially with charged lipids and different chain lengths as seen in real cells—may be more resistant to NP-induced damage. Nonetheless, another possible mode of structural alteration—inter-vesicle aggregation (Fig. 3D)—was noted at some low NP concentrations and poses potential risks around membranes. For example, organelles and transport vesicles within cells could be similarly affected by such aggregation. It should be noted that additional components in live E. coli cells, such as lipopolysaccharides and membrane proteins, will certainly influence these interactions. Due to the absence of these molecules in our model systems, their further impact will need to be examined using more sophisticated model membrane systems.
Nevertheless, these data underscore the potential for NPs to disrupt the delicate balance of membrane structures and functions. The observed structural and dynamic changes may impair membrane rigidity, protein functions, and vesicular transport. Significant biological consequences could occur, raising concerns about the environmental and health risks associated with nanoplastics.
The spherical morphology of the NPs used in this study may underestimate the effects that sharper-edged or irregularly shaped particles might exert on membranes. On the other hand, the NP concentrations employed here are several orders of magnitude higher than typical environmental nanoplastic levels. While these elevated concentrations may greatly amplify NP–membrane interactions, they primarily represent localized events that provide molecular-level insights into the interaction mechanisms. Additional surface properties could also exacerbate NP–membrane interactions, and these variables merit further exploration. Although a standardized approach to nanoplastic studies in this context has not yet been established, the present work demonstrates the utility of SANS contrast matching and offers key insights into NP–membrane interactions that are not accessible by other structural techniques such as cryo-EM.27 These findings not only advance the understanding of nanoplastics' potential to cause molecular damage but also lay the groundwork for future investigations into the biological impacts of environmentally relevant NPs at the molecular level, especially by neutron scattering techniques.
Supplementary information (SI): additional SANS and NSE analysis, addtional SEM and cryo-EM viewgraphs. See DOI: https://doi.org/10.1039/d5en00747j.
| This journal is © The Royal Society of Chemistry 2026 |