Lactoferricins impair the cytosolic membrane of Escherichia coli within a few seconds and accumulate inside the cell
Abstract
We report the real-time response of Escherichia coli to lactoferricin-derived antimicrobial peptides (AMPs) on length scales bridging microscopic cell sizes to nanoscopic lipid packing using millisecond time-resolved synchrotron small-angle X-ray scattering. Coupling a multiscale scattering data analysis to biophysical assays for peptide partitioning revealed that the AMPs rapidly permeabilize the cytosolic membrane within less than 3 s—much faster than previously considered. Final intracellular AMP concentrations of ∼80–100 mM suggest an efficient obstruction of physiologically important processes as the primary cause of bacterial killing. On the other hand, damage of the cell envelope and leakage occurred also at sublethal peptide concentrations, thus emerging as a collateral effect of AMP activity that does not kill the bacteria. This implies that the impairment of the membrane barrier is a necessary but not sufficient condition for microbial killing by lactoferricins. The most efficient AMP studied exceeds others in both speed of permeabilizing membranes and lowest intracellular peptide concentration needed to inhibit bacterial growth.
Editor's evaluation
This article presents groundbreaking data on the effects of a family of antimicrobial peptides on bacterial cells, obtained by time-resolved small-angle X-ray and neutron scattering experiments coupled to stopped-flow mixing. Application of this approach to cells is highly innovative. The main result is that the peptides reach the cytosol in a few seconds, where their accumulation at high concentrations finally kills the bacteria.
https://doi.org/10.7554/eLife.72850.sa0Introduction
Progress in designing antibiotics with novel key-lock mechanisms is not keeping pace with the worldwide growing number of (multi) resistant bacterial strains, encouraging significant research efforts in promising alternatives such as antimicrobial peptides (AMPs) (Lohner, 2001). AMPs are part of the natural innate immune system and provide a first line of defense against pathogens. Their advantage compared to conventional antibiotics relies on a rapid impairment of the barrier function of the bacterial envelope by unspecific physical interactions, often coupled to an ensuing targeting of bacterial DNA or ribosomes (for review, see, e.g., Wimley and Hristova, 2011; Lohner, 2017; Malanovic et al., 2020; Stella et al., 2021).
Membrane-active AMPs contain specific sequences of cationic and apolar amino acids, granting high affinity to the hydrophobic core of lipid membranes and selectivity toward the negatively charged surfaces of bacterial envelopes. However, despite intense research for several decades, a comprehensive understanding of the specific series of events that pertain to the bactericidal or bacteriostatic activity of AMPs is still elusive. To a large extent, this is due to the persisting challenge of merging results from in vitro studies with those obtained from lipid membrane mimics, often leading to significant controversies (Wimley and Hristova, 2011). This is nurtured, on the one hand, by difficulties in engineering lipid model systems of sufficiently high complexity to mimic the diverse physicochemical properties of bacterial membranes. On the other hand, the complexity of live bacteria challenges experimental and computational techniques to obtain quantitative results on the molecular level. For example, cryogenic transmission electron microscopy (TEM) provides high subcellular spatial resolution, but might give misleading information due to artifacts that potentially originate from staining or invasive sample preparation. Moreover, structural kinetics occurring in the seconds time scale are yet not accessible to cryo-TEM on cells, but would be needed to unravel the sequence of events induced by AMP activity. Kinetic experiments, using high-speed atomic force microscopy, for example, showed an AMP-induced corrugation of the bacterial outer surface within the first 13 s after addition of the AMP, but were limited by the intrinsic time resolution of the technique (Fantner et al., 2010). Additionally, such topological experiments do not provide insight into concurring intracellular changes. Video microscopy, combined with fluorescence labeling schemes for peptides or cellular content, in turn provided the appropriate spatiotemporal resolution to differentiate AMP activity in different cells within several tens of seconds (see Choi et al., 2016 for review). Importantly, such experiments reported that AMPs often reach the cytoplasm within a few minutes, suggesting that the final target for arresting bacterial growth or killing is not the cytoplasmic membrane (Sochacki et al., 2011). However, fluorescence labeling may easily tweak the delicate balance of macromolecular interactions and thus affect experimental observations.
Elastic X-ray and neutron scattering experiments are well-established noninvasive techniques to interrogate the structural properties of biological or biobased matter without the need to resort to any bulky labels. However, quantitative analysis of complex biologically relevant systems is challenging and has hampered progress in applying these techniques to live cells (Semeraro et al., 2021a). To this end, a full analysis of cells necessitates to account for structural features at diverse hierarchical levels. We recently refined an analytical multiscale scattering model of live Escherichia coli, making extensive use of the different sensitivities of X-rays and neutrons to matter, including H/D contrast variation (Semeraro et al., 2021b). This allowed us to detail the bacterial hierarchical structure on length scales bridging four orders of magnitude, that is, spanning from bacterial size to the molecular packing of lipopolysaccharides (LPS) in the outer leaflet of the outer membrane. Here, we use this model, taking advantage of the fact that the full breadth of structural information is encoded in a single scattering pattern, and exploit millisecond time-resolved synchrotron (ultra) small-angle X-ray scattering (USAXS/SAXS) to study the response of E. coli to three well-characterized lactoferricin-derived AMPs: LF11-215 (FWRIRIRR-NH2), LF11-324 (PFFWRIRIRR-NH2), and O-LF11-215 (octanoyl-FWRIRIRR-NH2) (Zweytick et al., 2011; Zweytick et al., 2014; Sánchez-Gómez et al., 2015; Marx et al., 2021b).
Joining these elastic scattering experiments with TEM and assays for determining peptide partitioning as a function of peptide activity enabled us to gain unprecedented insight into the peptide-induced sequence of events. We observed that the studied peptides are able to cause a loss of cytoplasmic content just within a few seconds. Coupling this finding to the derived final (after 1 hr) intracellular peptide concentrations at full bacterial growth inhibition (~80-100 mM) suggests a peptide uptake on similar time scales, and hence much faster than previously reported for other AMPs, such as rhodamine-labeled LL-37 (Sochacki et al., 2011). The most effective AMP presently studied, LF11-324, excelled others by most swiftly permeabilizing the cytoplamic membrane and killing the bacteria at the lowest cytosolic/periplasmic concentration. Consistent with previous studies on these AMPs (Zweytick et al., 2011), we also observed severe damage of the bacterial cell envelope—here defined in terms of loss of LPS packing, loss of positional correlations between outer and inner membranes and vesiculation/tubulation. Interestingly, this structural impairment of the bacterial envelope also occurred at peptide concentrations far below the minimum inhibitory concentration (MIC), but on different (slower) time scales than leakage. Leakage kinetics increased with peptide concentration up to the MIC. Yet, the overall loss of cytoplasmic molecules was not affected by AMP concentration, that is, it was equal at MIC and below (even at 1% growth inhibition!). The primary cause of bactericidal or bacteriostatic activity of the presently studied peptides is thus not a damage of the structural integrity of the cell wall, but appears to be due to a fast and efficient impairment of a series of here not further detailed physiological processes occurring within the intracellular space.
Results
Defining structural reference states of AMP activity in E. coli
Unraveling the timeline of structural events occurring in E. coli due to lactoferricin activity by USAXS/SAXS necessitates a detailed prior characterization of two reference states: (i) neat bacteria before AMP administration (‘initial state’) and (ii) AMP-affected/killed bacteria (‘end state’). Here, end state refers to 1 hr of incubation of bacteria at a given AMP concentration. Both states can be treated as quasi equilibrium structures, enabling a detailed characterization through a joint application of X-ray and neutron scattering. This allowed us to constrain parameters in the analysis of the time-resolved USAXS/SAXS data described further below. We have recently reported initial-state structures of E. coli ATCC 25922 (and other strains) at different hierarchical length scales—including size of bacteria, distance between inner and outer membranes, and LPS packing density—in terms of a multiscale analytical model using joint USAXS/SAXS and (very) small-angle neutron scattering (VSANS)/SANS experiments (Semeraro et al., 2021b). Importantly, while this analysis revealed a sensitivity of X-rays/neutrons to structural features of the cellular envelope or bacteria size, we also found that neither flagella or fimbriae, nor macromolecules located in the cytosol (DNA, ribosomes, proteins, etc.) contribute discernible scattering patterns and hence can be detected.
The same methodological concept was applied here to reveal the end-state structure of E. coli. Further, and analogously to Semeraro et al., 2021b, we coupled the USAXS/SAXS and SANS analysis to TEM in order to remove ambiguities in adjustable parameters originating from ensemble averaging (Figure 1). In order to couple scattering experiments to AMP susceptibility assays, it is important to scale peptide concentrations appropriately. Firstly, because SAXS/SANS experiments require up to 104 times higher bacterial concentrations than standard growth-inhibition experiments. Secondly, because the total number of AMPs partitioning into bacteria (and hence also their antimicrobial activity) depends on cell concentration in a nontrivial manner (Marx et al., 2021b). We consequently determined the MICs of all here-studied AMPs as a function of cell number density prior to scattering experiments; the corresponding data are reported in Figure 1—figure supplement 1.
Small-angle scattering (SAS) patterns of initial and end states showed distinct differences, which are highlighted in Figure 1—figure supplement 2E and F. Figure 1—figure supplement 2E, in particular, demonstrates a decrease of intensity at very low q-values (as approximation of forward scattering), and a faster intensity decay at high for the end state, coupled to a loss of two intensity wiggles at and (see also Figure 1A). The comparison of SANS data at 90% D2O (Figure 1—figure supplement 2F) instead shows a smoothening of an intensity shoulder at coupled to a slight peak-shift toward lower q-values for the end state. The application of our previously developed analytical model for E. coli (Semeraro et al., 2021b) allowed us to gain detailed insight into the ultrastructural changes pertaining to these differences in scattered intensities. To achieve this goal, we jointly analyzed USAXS/SAXS and SANS data at five different contrasts (10, 30, 40, 50, and 90% D2O; see Figure 1—figure supplement 2). Taking AMP-induced changes of bacterial ultrastructure as observed by TEM into account then leads to the following observations.
Membrane ruffling, mainly originating from increased fluctuations of cytoplasmic membranes clearly seen by TEM, can be modeled by changing the average distance between inner and outer membranes and its fluctuation (Figure 1B), which readily accounts for the above-described change of the SANS feature in Figure 1—figure supplement 2F. Consistent with previous reports (Sochacki et al., 2011), we also observed an overall shrinking of bacterial size. This correlates with a loss of periplasmic or cytosolic material. Both effects explain the changes in USAXS intensity at very low q-values (Figure 1—figure supplement 2E).
The scattering shoulder in USAXS/SAXS data at originates from positional correlations between LPS oligosaccharide cores, as revealed by our previous in-depth analysis (Semeraro et al., 2021b). Its smearing out in the presence of peptide consequently indicates a loss of LPS packing correlations, either because of a decrease of the number of LPS molecules on the outer surface or due to increased membrane roughness or waviness. This effect occurs on sub-nanometer length scales and is thus hidden from TEM. Notably, the end-state scattering patterns in the presence of either LF11-324 or LF11-215 were identical, despite their distinct MIC values (data not shown). Modeling their scattering data revealed an additional contribution, which we could associate with the help of TEM to the formation of outer membrane vesicles (OMVs), having an average diameter of ~30 nm (Figure 1A and B, Figure 1—figure supplement 2E). Instead, the separation of the nucleotide region from the nucleotide-free cytosol revealed by TEM (Figure 1A), and in agreement with a previous report on same AMPs (Zweytick et al., 2011), is not observed in our scattering data. As noted above, SAS is not sensitive to macromolecules and aggregates thereof inside the cell (Semeraro et al., 2021b). Finally, O-LF11-215 led to additional effects. First of all, O-LF11-215, because of its increased hydrophobicity, forms aggregates in buffer as clearly shown by SAXS (Appendix 1). In particular, incubating O-LF11-215 with bacteria at concentrations needed for SAS experiments led to the formation of macroscopic aggregates, which impeded the measurements of end states by SAS. TEM experiments, however, showed the formation of extramembranous tubes (Figure 1C).
All end-state SAS data were fitted with our analytical model using a Monte Carlo genetic-selection algorithm; see Figure 1—figure supplement 4 for a graphical scheme of all adjustable parameters. Compared to standard least-square algorithms, this allowed us to obtain robust parameter values from probability density functions and to judge possible parameter correlations (for details, see ‘Materials and methods’). Based on our previous modeling of neat E. coli (Semeraro et al., 2021b), the individual SANS contrasts (Figure 1—figure supplement 2A) are motivated as follows: (i) the distance between the inner and outer membranes and its fluctuations, , is accentuated at 90, 50, and 40% D2O; (ii) 30–50% D2O highlight the scattering length density (SLD) balance between the cytosol , periplasm , and the peptidoglycan layer ; and (iii) 10% D2O yields a similar contrast than X-ray and thus provides a direct and independent control for USAXS/SAXS experiments. Table 1 summarizes the results of this analysis for LF11-324 by listing the changes between initial and end states. As discussed above, our analysis also reports overall bacterial size changes. E. coli are reasonably well represented in reciprocal space using an ellipsoid with multiple shells to account for the different compartments (Semeraro et al., 2021b). Since our smallest probed scattering vector, (corresponding to distances of about 2 µm), does not allow to probe the full length of the bacteria, we report size changes in terms of changes of the minor radius of the ellipsoid and, in particular, that of the cytoplasmic space . The major radius was coupled to by fixing the major to minor radius ratio, , to a value of 2 (see Table 2). Finally, the USAXS/SAXS scattering shoulder at reports on positional correlations between LPS oligosaccharide cores (Semeraro et al., 2021b). In our analytical model, this is captured by the number of oligosaccharide cores (OS), . To characterize the peptide-induced changes, we define the LPS packing parameter, , as the ratio between the number of OS (i.e., LPS) in the current and initial state (). Other parameters of the analytical model, such as the thicknesses of the hydrophobic and hydrophilic layers of the cytoplasmic and outer membranes, as well as their average SLDs, the thickness of the peptidoglycan layer, and the radius of gyration of the oligosaccharides, were fixed at the same values as reported in Semeraro et al., 2021b (see also Tables 2 and 3). We additionally checked for contributions of peptides to the SLDs of each bacterial compartment, including cytoplasmic (CM) and outer membranes (OM). However, the associated changes were found to be insignificant.
The observed LF11-324-induced changes to E. coli show, upon the application of our model, a decrease of , along with the increase of , signifying leakage of mainly low-weight molecules and ions from the cytoplasm (Figure 1—figure supplement 2B and C), whose scattering dominate both and (Semeraro et al., 2021b). Surprisingly, final and values did not depend on peptide concentration, that is, were equal within experimental uncertainty at MIC and sub-MICs (see Figure 2B). The overall increase of the distance between inner and outer membranes, , corresponds to about 18%; its fluctuations increase by almost a factor of two. The associated disorder in the periplasm then also accounts for the drop of the contrast of the peptidoglycan layer by about 2.6%. The LPS packing parameter, , instead drops dramatically by 30–40%. Finally, the observed cell shrinkage amounts to ∼5%, which leads to a decrease of cell surface of approximately 2 ×106 nm2. Apparently this is at least in part compensated by OMV formation; indeed the overall cell-surface decrease is on the same order of magnitude as the estimate of the total surface of all OMVs of nm2 as obtained from analyzing SAS data (see Appendix 2).
Kinetics: Time-resolved USAXS/SAXS
The structural transitions from initial to end state were followed by USAXS/SAXS at millisecond time resolution. Stopped-flow mixing ensured thorough and rapid re-dispersion (mixing time of 50 ms) of peptides and bacteria (Figure 2—figure supplement 1) and led to immediate changes of scattering patterns. Using the initial and end-state values obtained for the adjustable parameters as constraints allowed us to fit kinetic scattering data as detailed in the ‘Materials and methods’. A close examination of the results obtained for the three LF11s (Figure 2, Figure 2—figure supplement 2, and Figure 2—figure supplement 3) allowed to discern AMP concentration-independent parameters (, , , and ) from AMP concentration-dependent parameters (, , and ).
First, we focus on LF11-324 and AMP concentration-independent parameters. The packing of LPS started to decrease at after mixing (Figure 2A). Changes of , , and in turn are largely decoupled from this remodeling of the outer membrane, with a common onset of 2–10 min after peptide addition (Figure 2D–F). AMP concentration-dependent parameters instead showed an increasing delay of changes upon decreasing the amount of administered peptide (Figure 2B and C). In particular, exhibited a pronounced increase already at at highest peptide concentration shifting to min at the presently lowest studied peptide concentration. The onset times of decrease of do not appear to be directly correlated with these changes, but also shifted progressively to later times with decreasing AMPs. Because of the dominant contribution of low-molecular-weight molecules and ions to both and (Semeraro et al., 2021b), we surmise that these changes are due to a leakage of inner and outer membranes. Briefly, the cytoplasmic content diffuses first in the periplasmic space and, simultaneously, material from the periplasm leaks out of the cell. This process leads to the rather simultaneous decrease of and increase of . Because of differences in the individual onsets of these trends at a given AMP concentration, we can only associate a time range for the beginning of permeation of the cytoplasmic membrane: 3–10 s at , 10–20 s at , and 50–120 s at . Despite these differences, and as noted above, (and ) reached the same final values for all three AMP concentrations. Finally, the drop of started at 20–50 s for and , and >10 min for (Figure 2C).
Interestingly, LF11-215 led to almost identical kinetics for the concentration-independent parameters (, , and in Figure 2—figure supplement 2A–C and in Figure 2—figure supplement 3A). This suggests that LF11-215 and LF11-324 remodel the membranes and affect the macroscopic stability of the cell envelope in a very similar fashion, although we found no decrease of for LF11-215 at (Figure 2—figure supplement 3A). Importantly, and started to change much later for this peptide ( at , Figure 2—figure supplement 3B). The overall changes of and indicate, however, that the AMP concentration-induced loss of cytoplasmic content is a common feature of LF11-324 and LF11-215.
The kinetics of the cell-envelope damage caused by O-LF11-215 were found to proceed analogously to both LF11-324 and LF11-215 (Figure 2—figure supplement 2D–F). However, while changes of and occurred at about similar times than for LF11-324 (Figure 2B and H), (Δt ~ 2 min, Figure 2G) and (Δt ~ 10 min, Figure 2I) changed much later. The initial increase of is likely to be an artifact, probably due to the convolution of scattering originating from cells and peptide aggregates (see Figure 2—figure supplement 1D).
LF11-215/LF11-324-induced OMV formation was discernible for Δt > 1–2 min (Figure 3A) and appears to be peptide- and concentration-independent within the first 10 min. Finally, the peptide cluster term, introduced for the analysis of O-LF11-215 (see ‘Materials and methods’), enabled us to estimate that a large increase of peptide uptake starts after about 2 min, which, however, did not depend on peptide concentration (Figure 3B).
Peptide partitioning and cooperativity
We applied a previously detailed assay for AMP partitioning in E. coli based on growth inhibition (Marx et al., 2021b). A statistical analysis of the corresponding data in terms of cumulative distribution functions (CDFs, see Appendix 3) allowed us to derive probability density functions (PDFs), describing bacterial growth inhibition as a function of cell concentration , including the minimum AMP concentration needed for inhibiting a given percentage of E. coli, ICx (Figure 4); note that . We obtained Poisson-like PDFs of ‘killing’ events as a function of peptide concentration from the sigmoidal shapes of the measured CDFs (Appendix 3). Besides the specific mode of action of a peptide, the shape, width, and position of the PDFs are a result of intrinsic statistical variations and fluctuations of the systems’ properties, for example, bacterial size, stage of mitosis, local concentration of partitioned peptides, biological variability within the populations, etc.
In agreement with our previous report for E. coli K12 (Marx et al., 2021b), the MIC of LF11-324 is lower than that of LF11-215 at all cell concentrations. For O-LF11-215, MIC values matched those of LF11-324 at low cell concentrations, but increased strongly with , finally superseding that of LF11-215 and becoming immeasurably high due to the above-described peptide aggregation (Figure 1—figure supplement 1). The growth inhibition probabilities for LF11-324 and LF11-215 peak close to the IC50s (Figure 4). The distributions are much sharper, that is, have a smaller full-width-at-half-maximum (FWHM) of the probability distributions, , for LF11-324 than for LF11-215. This suggests an increased ‘cooperativity’ of killing for LF11-324, in the sense of a two-state transition model between alive and dead bacteria (see Appendix 3). increased with cell concentration, for example, from at CFU/ml to at CFU/ml for LF11-324. Significant noise levels in growth inhibition data for O-LF11-215 impeded a determination of killing probabilities at inhibitory concentrations . However, data retrieved at higher inhibitory concentrations suggest that the probability distributions roughly match those of LF11-324 at low cell concentrations, but become broader than that of LF11-215 at high cell content (Figure 4—figure supplement 1). This is a signature of loss of killing efficacy at high , most likely due to peptide self-aggregation as discussed above.
Next, we derived for each IC (or growth-inhibited fraction, ) the average number of cell-partitioned peptides, , and the effective partitioning coefficient, , applying a previously reported thermodynamic formalism (Marx et al., 2021b; see also Equation 1). increased for all three peptides with , although the changes of were smallest for LF11-324, followed by LF11-215 and O-LF11-215 (Figure 5A). These results were confirmed independently also by Trp-fluorescence spectroscopy (see Appendix 4).
This partitioning behavior was also mirrored in the dependence of . was nearly constant for LF11-324, increased only slightly for LF11-215, and showed the largest variation for O-LF11-215, reaching about 2.5 times higher levels than the other two peptides (Figure 5B). The approximate equal values of LF11-324 and LF11-215 for demonstrate that both peptides partition about equally well into E. coli, not only at the MIC, but in a wide range of values.
ζ-potential measurements were performed ∼5 min and 1 hr (end state) after adding AMPs. Results did not show significant differences for the two time points. This enabled us to further discriminate the activities of LF11-324 and LF11-215. We observed an initial increase of ( refers to the reference system, i.e., neat bacteria) at low peptide concentrations, reaching plateau values of for LF11-215 and for LF11-324 for (Figure 5—figure supplement 1). As detailed in Marx et al., 2021b, these data allowed us to calculate the maximum number of peptides associated to the LPS leaflet, . We found and for LF11-215 and LF11-324, respectively, in the -plateau regions. Normalizing these results by the number of cell-partitioned peptides () reveals an overall decrease of with (see Figure 5C). Note that the -potential analysis provides (with the applied mathematical and physical constrains) upper-boundary values for (see ‘Materials and methods’). That is, . Thus, our conservative ratio strikingly demonstrates that most of the peptides are located within the intracellular compartments of E. coli at the MIC. The fraction of outer-leaflet-partitioned peptides increased toward lower , and somewhat stronger for LF11-215, but does not exceed 10%. An analogous analysis for O-LF11-215 was impeded by the peptide aggregates, whose sizes were in the same order or even larger than that of bacteria (Figure 5—figure supplement 2).
Discussion
The entire set of the above-detailed results reveals a complex scenario that needs to be evaluated with utmost care. To this end, we stress that all techniques used in this work report ensemble averages. That is, they account for the entire bacterial population within the samples. For example, an LPS packing value of —observed for growth-inhibited fractions of 1%, as well as 99.9% after 1 hr (Figure 2A)—signifies that the ensemble average shows the same loss of LPS packing under both conditions, irrespective of individual differences in some isolated cells.
A particularly striking result for the LF11-215/LF11-324 end states is that the peptide-induced effects are similar and independent of peptide concentration (Figures 1 and 2). That is, even at growth-inhibited fractions of just 1%, we observed much the same cellular permeabilization and structural changes of the bacterial ultrastructure as at quasi fully growth-inhibited E. coli (see also Figure 4 and Figure 4—figure supplement 1). Here, it is important to bear in mind that, unlike our susceptibility assays, we did not add a growth medium after incubation with peptides in SAS experiments. Our combined SAS and growth inhibition data provide unambiguous evidence that bacteria are able to recover at sub-MIC concentrations from a severe collateral damage of their cell envelope. Consequently, this damage cannot be the primary killing cause of bacteria. Our structural studies only revealed peptide-specific and concentration-dependent effects for the kinetics of , , and (Figure 2). It is unlikely though that antimicrobial activity is specific to the rate of bacteria size changes. In turn, the different onsets of changes of and may provide some clues to events that finally inhibit the growth of E. coli. In what follows, we will thus focus on the sequence of events that can be deduced from our study for LF11-324 and LF11-215. Before doing so, we note that also O-LF11-215 caused comparable variations of the abovementioned parameters on similar time scales. In this case, however, an equally detailed analysis was impeded by the propensity of O-LF11-215 to aggregate in buffer solution.
Figure 6 provides a scheme of the timeline of events revealed by this study, taking also into account our complementary experiments on peptide partitioning. Within a few seconds, LF11-324 and LF11-215 cause changes of the LPS packing density (Δt ~ 10 s), as well as periplasmic and cytoplasmic SLDs. Since ions and small molecules (<1 kDa) are the major contributors to its X-ray SLD (Semeraro et al., 2021b), the decrease of will be dominated by a leakage of these entities (e.g., ATP). Molecules leaking from the cytosol first diffuse into the periplasm and then further into extracellular space, leading to a net increase of . Leaking into extracellular space follows from the observation that final levels do not reach those of , despite the much larger reservoir for cytosolic molecules. Further, the initial SLDs of buffer and periplasm are comparable, also explaining why our technique is not directly detecting outer membrane leakage. Hence, either observed change of or is due to a permeabilization of both cytoplasmic and outer membranes.
For LF11-324, the permeabilization of the cytoplasmic membrane occurred as fast as 3–10 s after mixing at [P] = 1.2× MIC. Dropping peptide concentration led to a slowing down of this effect (10–20 s for [P] = 0.7× MIC, and 50–120 s for [P] = 0.3× MIC), but did not affect the final cytoplasmic density, that is, the overall loss of material (Figure 2). AMPs need to translocate all the way through the cell wall in order to induce such effects, implying that peptide translocation of the outer membrane proceeds on time scales faster than or . We stress that the usage of the term ‘translocate’ does not imply that AMPs are able to pass through inner and outer membranes without any noticeable effect on membrane structure, such as the transient formation of pores, nor does it exclude their presence. It merely refers to a generic, unspecified uptake of the peptides beyond the resolution of the present experiments. Further, since we are not able to discriminate membrane peptide adsorption kinetics from the onset of leakage of cellular content, we cannot comment on a potential ‘carpet mechanism,’ proposed from analogous AMP partitioning studies in bacteria (Roversi et al., 2014; Savini et al., 2020). We note, however, that LF11-215 was found to partition into artificial lipid membranes, without noticeable effects on lipid flip-flop (Marx et al., 2021a). We detail this discussion further below.
The drop of is a natural consequence of the loss of cellular content, but occurs at somewhat later times. This delay is most likely due to the stored elastic energies of the peptidoglycan layer, which will initially resist rapid deformations (Yang et al., 2018). Note, however, that the peptidoglycan properties are likely affected by direct interactions with peptides (Zhu et al., 2019). Remarkably, cytoplasmic membrane leakage in the MIC range occurs ~1 min later for LF11-215 (Figure 2—figure supplement 3).
Pronounced differences between LF11-324 and LF11-215 were also observed from their efficacies as a function of cell concentration (Figure 4). At equal , growth-inhibition probability distributions are much narrower for LF11-324. Apparently, this increased ‘cooperativity’ correlates with the peptide’s ability to permeabilize the cytoplasmic membrane faster. It is further illuminating to discuss the total amount of peptide penetrating into the intracellular space. Both peptides saturate the outer LPS leaflet already at concentrations lower than 0.3× MIC (Figure 5—figure supplement 1), corresponding to about one AMP per 4–5 LPS molecules as upper-bound value (see in ‘Materials and methods’). Thus, consistent with Figure 2A and B, peptides penetrate the outer membrane already when the effects on bacterial growth are still very small. Assuming that the inner leaflet of the outer membrane, as well as the two leaflets of the inner membrane, host to first-order similar amounts of LF11-324 and LF11-215, we calculate that ≥92% of all peptides are located in the intracellular space at the respective MICs. Using our USAXS/SAXS/DLS data for the average cytoplasmic and periplasmic volumes, we further arrive at huge intracellular concentrations of ~80 mM for LF11-324 and ~100 mM for LF11-215. The relative difference in cytosolic/periplasmic concentrations between the two peptides corresponds to about the relative difference in MICs (Figure 1—figure supplement 1). Our finding that only a minor fraction of the peptides is located within the membranes remains valid even at sub-MIC values, offering an attractive explanation for the similar and concentration-independent average cell-wall damage-induced LF11-324 and LF11-215. Specifically, we propose that this damage (LPS packing, OMV formation, membrane ruffling and detachment) is a consequence of the small fraction of peptides that saturate the membranes and, possibly, the peptidoglycan layer. Following the thermodynamic principles of the partitioning framework and ζ-potential measurements, this small fraction is supposed to be rather constant above AMP concentrations of about 0.2× MIC. Along these lines, recent in-cell FRET experiments showed that only a small fraction of the peptides partitioned in a cell is actually interacting with the bacterial surface (Kaji et al., 2021).
About 10 years ago, Wimley, 2010 already discussed the consequences of having 107 to 109 AMPs per cell, suggesting a ‘reservoir’ of non-membrane-bound peptides that would outnumber proteins, ATP, and other metabolites. Cytosolic targets were also confirmed by our TEM data (and those of others; Hammer et al., 2010; Scheenstra et al., 2019) showing a collapsed nucleoid region (Figure 1). For instance, this could be due to a competition mechanism with polyamines as putative stabilizers of the functional architecture of the DNA ring (Hou et al., 2001). In support of the hypothesis of cytosolic targets, recent solid-state 31P-NMR measurements of E. coli in the presence of AMPs revealed increased dynamics of DNA and RNA phosphate groups correlated to TEM observations of collapsed nucleoid volume (Overall et al., 2019). In addition, however, AMP interactions with other negatively charged metabolites and macromolecules in the cytosol, such as ribosomes and proteins, should be considered as potentially detrimental to bacteria (Zhu et al., 2019). The about 1.4 times lower intracellular concentration of LF11-324 then is a signature of a higher potency compared to LF11-215, possibly hampering a number of physiological functions (Scocchi et al., 2016), or inducing an apoptosis-like mechanism (Dwyer et al., 2012) due to interactions with any of the above-named cytosolic molecules. The exact killing mechanism of the here-studied AMPs, however, evades current experimental design.
It follows from the considerations above that bacteria have an increased probability to recover from the peptides’ attack if the intracellular concentrations of LF11-324 and LF11-215 fall below those reported above (i.e., and mM, respectively). Thus, both the ability of LF11-324 to swiftly permeabilize the cytoplasmic membrane by translocating through the cell envelope and its higher propensity to interfere with physiological processes contribute to its higher ‘cooperativity’ in killing E. coli (Figure 4). Our experimental setup is not sensitive to directly observe either transient membrane pores or other membrane defects. Note, however, that previous studies in lipid-only mimics of bacterial membranes showed only weak membrane remodeling effects of lactoferricin derivatives compared to other peptides (Zweytick et al., 2011; Marx et al., 2021b). We propose that the here-studied LF11s follow the interfacial activity hypothesis, valid for both antimicrobial (Wimley, 2010) and cell-penetrating peptides (Kauffman et al., 2015). That is, they are able to locally and transiently induce disorder in the lipids’ hydrocarbon chain regime, which sufficiently affects the permeability barrier of the bilayer against ions, small polar molecules, and, eventually, the peptides themselves. This might be also a consequence of the rather short amino acid sequence of the lactoferricin derivatives, implying a too small length (~1 nm) to span the whole membrane thickness (3–4 nm). In turn, the highly flexible secondary structure and an amphipathic moment being aligned along the peptide’s backbone (Zorko et al., 2009; Zweytick et al., 2011; Zweytick et al., 2014) should allow LF11-324 and LF11-215 to translocate membranes at higher rates than observed for linear peptides (Ulmschneider, 2017; Kabelka and Vácha, 2018). Note that peptide translocation can also lead to transient membrane leakage events (Ulmschneider, 2017), even with negligible AMP-induced lipid flip-flop (Marx et al., 2021a). The higher hydrophobicity of O-LF11-215 should increase the likelihood of remaining membrane bound, which might build up differential membrane curvature stress and lead to the observed formation of membrane tubules. We also note that the different leakage kinetics for the LF11 peptides suggest a coupling to translocation kinetics, which in turn depends on membrane partitioning of the AMPs. Indeed, recently reported data for cytoplasmic membrane mimics of cardiolipin, phosphatidylethanolamine, and phosphatidylglycerol show a somewhat faster membrane partitioning for LF11-324 than for LF11-215 (Marx et al., 2021b).
In conclusion, the superior time resolution and sensitivity to ensemble averaged structural changes from cellular size to molecular packing of synchrotron USAXS/SAXS allowed us to demonstrate, upon combination with advanced data modeling and complementary techniques for peptide partitioning, that LF11-based AMPs are able to reach the cytosolic membrane of bacteria on the seconds time scale and subsequently accumulate at high concentrations (80–100 mM) in all bacterial compartments, including the cytosol. Based on differences observed for LF11-324 and LF11-215, two factors emerge as key components for their antimicrobial efficacy: (i) a fast translocation through inner and outer membranes, and (ii) an efficient shutdown of here not further detailed vital chemical/biochemical processes, that is the lower the number of ‘needed’ AMPs the better (here achieved by LF11-324). The latter conclusion is corroborated by the most striking fact that the final cell-wall damage does not depend on peptide concentration. That is, AMPs inhibiting only 1% of bacteria induce on average—besides some single isolated variations—as much membrane ruffling, detachment of inner and outer membranes, leakage, etc., as AMPs at full levels of growth inhibition. Hence, the collateral damage of the cell-wall structure, including leakage of cellular content, is a ‘by-product’ of AMP activity. Instead, we suggest that the primary cause of growth inhibition is driven by interactions of the AMPs with yet-to-be-determined cytosolic/periplasmic molecules. Most likely candidates are the polyanionic DNA, RNA, ribosomes, and proteins (Zhu et al., 2019), or charged metabolites. We strongly emphasize that the ability of the here-studied peptides to impair the barrier function of bacterial membranes remains an unquestioned highly important property. However, at the same time our study shows that this effect alone is not sufficient to explain all of our data.
It is currently not clear whether the present findings can also be extended to other AMPs. Yet, the importance of the cell-penetrating ability was similarly discussed for, for example, peptides Sub3 and LL-37 (Torcato et al., 2013; Zhu et al., 2019). Further, a number of AMPs have been reported to show partitioning behavior in bacteria and accumulation in the intracellular volume analogously to the here-studied lactoferricin derivatives (see, e.g., Sochacki et al., 2011; Zhu et al., 2019; Savini et al., 2020; Loffredo et al., 2021; Kaji et al., 2021, and references therein). Le et al., 2017 gave a more complete account of potential intracellular targets of AMPs—either of fundamental or secondary importance for bactericidal activity. We thus propose that the combination of membrane translocation speed and efficient impairment of physiological processes are generic factors that should be considered in designing future AMPs to combat infectious diseases. This also implies a widening of the pure focus on membrane activity of AMPs currently applied in many studies.
Materials and methods
Samples
E. coli ATCC 25922 were provided by the American Type Culture Collection (Manassas, VA). Freeze-dried peptide powders of LF11-215 (H-FWRIRIRR-NH2), LF11-324 (H-PFFWRIRIRR-NH2), and O-LF11-215 (octanoyl-FWRIRIRR-NH2), purity >95%, were purchased from the Polypeptide Laboratories (San Diego, CA). Lysogeny broth (LB)-agar and LB medium were obtained from Carl Roth (Karlsruhe, Germany). All the other chemicals were purchased from Sigma-Aldrich (Vienna, Austria).
Bacterial cultures
Request a detailed protocolBacterial colonies of E. coli ATCC 25922 were grown in LB-agar plates at 37°C. Overnight cultures (ONCs) were prepared by inoculating a single colony in 3 ml LB-medium in sterile polypropylene conical tubes (15 ml), allowing for growth under aerobic conditions for 12–16 hr in a shaking incubator at 37°C. Main cultures (MCs) were then prepared by diluting ONCs in 10 ml LB-medium in 50 ml sterile polypropylene conical tubes. Bacteria in the MCs grew under the same conditions applied to ONCs up to the middle of the exponential growth phase. Cells were then immediately washed twice and resuspended in nutrient-free and isotonic phosphate-buffered saline (PBS) solution (phosphate buffer 20 mM, NaCl 130 mM) at pH 7.4.
AMP samples
Request a detailed protocolLF11-324 and O-LF11-215 peptides displayed a weak solubility in PBS. AMP stocks (including LF11-215) were then prepared by adding acetic acid and DMSO, up to 0.3 and 3% vol/vol, respectively. Peptide stock solutions were diluted for measurements yielding a final concentration of 0.01% acetic acid and 0.1% vol/vol dimethyl sulfoxide (DMSO) (final pH 7.2). Hence, possible effects of DMSO on the cell envelope, as observed for model membrane structures (Gironi et al., 2020), can be neglected. Control USAXS/SAXS experiments adding a similar amount DMSO and acetic acid to E. coli also showed no discernible effect of the organic solvent (data not shown). Similarly, control experiments were performed to exclude effects on bacterial growth. Stock concentrations were determined by measuring the absorption band of the Trp residues with the spectrophotometer NanoDrop ND-1000 (Thermo Fisher Scientific, Waltham, MA). Peptide stock solutions were stored in silanized glass tubes until use.
Antimicrobial activity and partitioning modeling
The antimicrobial activity of the AMPs on E. coli was tested in the bacterial concentration range of 5 × 105 to 109 CFU/ml using a modified susceptibility microdilution assay (Jorgensen and Ferraro, 2009). Cell suspensions were incubated with AMPs in buffer at a given peptide and cell concentration for 2 hr at 37°C (control samples were incubated in buffer only). Cell growth was monitored upon addition of double concentrated LB-medium for about 20 hr using a Bioscreen C MBR (Oy Growth Curves Ab, Helsinki, Finland).
Partitioning modeling
Request a detailed protocolThe analysis of the inhibited fraction of cells, , as a function of the total concentration of peptides, , enabled the extraction of ICx values as a function of , as detailed in Appendix 3. Following previously reported methods to determine peptide partition in liposomes (White et al., 1998) or bacterial systems (Savini et al., 2017), we derived the number of cell partitioned AMPs for a given ICx value using
where and are the concentrations of AMPs dispersed in the aqueous phase and segregated into the cells, respectively; is the Avogadro’s constant; is the concentration of water molecules in bulk (55.3 M at 37°C); is the effective mole-fraction partitioning coefficient; and is the average number of peptide monomers that are partitioned within a single cell. Details of this analysis are reported in Marx et al., 2021b.
A similar approach was exploited in the case of O-LF11-215 peptide clusters in solution. In this case, the total peptide concentration is given by , where is the molar concentration of aggregates, each of them consisting of an average number of peptides . We also define the aggregate fraction and assume the equilibrium state . The definition of the molar partitioning coefficient (Marx et al., 2021b) refers to the balance of concentration of free peptides in bulk and partitioned peptide into the cells. Hence, its bare definition is unaffected by the presence of clusters. Finally, it is trivial to show that in this case Equation 1 becomes
where the fitting parameter is an upper boundary value for the actual number of partitioned peptides per cell .
ζ-potential, cell size, and outer leaflet distribution of peptides
Request a detailed protocolζ-potential and dynamic light scattering (DLS) measurements were carried out with the Zetasizer Nano ZSP (Malvern Panalytical, Malvern, UK). E. coli suspensions were incubated with different concentrations of AMPs in buffer for 1 hr at 37°C prior to each measurements. Control samples (no AMPs) were suspended and incubated in buffer. Measurements were also repeated just after mixing with peptides (about 5 min waiting time for sample loading, data not shown), exhibiting negligible variation from samples incubated for 1 hr. A concentration of 107 CFU/ml provides the optimal compromise between high signal-to-noise ratio and low multiple-scattering bias. The AMP concentrations were centered on the MIC values previously determined with the susceptibility microdilution assays, spanning from about 0.2× to 2.5× MIC. For ζ-potential measurements, the voltage for the electrodes was set to 4 V, such that currents did not exceed 1 mA, because of the high conductivity of the PBS buffer. Further, measurements were paused between repetitions for 180 s. This prevented heat productions, leading to sample denaturation and accumulation on the electrodes. The experiments were repeated three times, and, due to the low sensitivity of such a setup, each of them consisted of a minimum of six measurements (see also Marx et al., 2021b). For each system, ζ-potential values and associated errors were given by the medians and the median absolute deviations, respectively, averaging over at least 18 repetitions. The same number of scans was also used to obtain mean and standard deviation values for the hydrodynamic diameter, , of the cells.
From the measured ζ-potential values, we calculated the maximum number of peptides that partition into the outer LPS leaflet, , as reported recently (Marx et al., 2021b)
where (Wiese et al., 1998) and (Zweytick et al., 2011) are the nominal charges of LPS and LF11-215 or LF11-324 AMPs, respectively; ζ and are the ζ-potential and total surface values of the system upon addition of peptides; and and S0 are the respective reference values (no AMPs). is the estimated number of LPS molecules, where (Kim et al., 2016; Micciulla et al., 2019) is the lateral area per LPS molecule. The prefactor originates from considering a maximum surface coverage of 90% by LPS molecules (Seltmann and Holst, 2002). S0 was derived from DLS measurements, approximating the bacterial shape by a cylinder , considering that the hydrodynamic radius is approximately equivalent to the radius of gyration for micron-sized objects. Then fixing the radius at ~400 nm (Semeraro et al., 2021b) and retrieving the length from the above relation for one obtains .
Fluorescence spectroscopy
Fluorescence spectroscopy experiments were done with the Cary Eclipse Fluorescence Spectrophotometer (Varian/Agilent Technologies, Palo Alto, CA). The excitation wavelength was set to λ = 280 nm (which corresponds to the maximum intensity of the absorption/excitation band of Trp), and emission spectra were acquired in the λ-range between 290 and 500 nm, with the Trp emission band peak being expected to lie around 330–350 nm. Samples were loaded in quartz cuvettes of 1 cm path length. The background was subtracted from every Trp spectrum prior to further analysis.
Peptide solubility
Request a detailed protocolTrp emission allowed determining whether LF11 peptides form aggregates or not in the MIC range. Spectra of LF11-only samples at [P] = 100 µg/ml were fitted with the log-normal-like function (Burstein and Emelyanenko, 1996; Ladokhin et al., 2000):
where and I0 are, respectively, wavelength and intensity of the emission peak; is the FWHM of the band; is a skewness parameter (fixed at an optimum value of 1.36 after testing; and ).
Both LF11-215 and LF11-324 (see Appendix 1) showed a peak at about nm and nm. This is consistent with a location of the Trp residues in polar chemical environments having full mobility and thus suggests that these AMPs are monomeric (Burstein and Emelyanenko, 1996). In contrast, the acylated O-LF11-215 showed a significant blue shift related to a location of Trp within apolar surroundings (Ladokhin et al., 2000), indicating the formation of peptide aggregates.
Peptide partitioning
Request a detailed protocolAnalogously to the partitioning analysis performed for lipid-only membranes (Marx et al., 2021b), the Trp emission of bacteria/AMPs mixtures were treated as a two-component signal, one coming from the peptides in the aqueous phase, and the second one from AMPs interacting with the cells.
Bacterial suspensions were incubated with different concentrations of AMPs in buffer for 1 hr at 37°C (incubator Thermomixer C, Eppendorf, Germany). Reference samples (no AMPs) were suspended and also incubated in PBS. Experiments were performed at cell concentrations of 5 × 107, 108, and 5 × 108 CFU/ml, and AMPs amounts equal to 0.5×, 1×, and 2× MIC (each experiment was repeated three times). Fluorescence intensities were background-subtracted using the bacteria-only reference spectra at the corresponding concentrations. This enabled us to subtract the average signal from the aromatic residues in the cells, and the scattering arising from the high concentration of the cell suspensions. Emission spectra of bare bacterial systems surprisingly showed a single and clean Trp-emission band, peaking at nm ( nm). The peak intensity increased linearly from 4.9 (arb. units) at 107 cells/ml to 122 (arb. units) at 5 × 108 cells/ml; then turbidity effects had to be accounted for. LF11-215 and LF11-324 showed peaks of comparable intensities even in the high cell-concentration regime, enabling us to retrieve background-subtracted spectra of optimal signal-to-noise ratio. Spectra were analyzed with a linear combination of two independent bands (see Equation 4) and , referring to AMPs in bulk (W) and partitioned into the lipid bilayer (B). and were fixed to the reference values obtained by analyzing spectra from pure AMPs, and , , , and were freely adjusted. LF11-only solutions were measured to calibrate their intensity dependence in buffer. Then, the retrieved values were converted to the concentration of dissociated peptides . This allowed us to obtain the so-called peptide bound fraction as (see Appendix 4). The aggregation of O-LF11-215 led to low values (see Appendix 1) and precluded a similar analysis for these peptides.
Transmission electron microscopy
Request a detailed protocolE. coli suspensions at 5 × 108 CFU/ml were mixed with peptides at the corresponding MICs and 0.5× MICs in PBS, and incubated for 1 hr at 37°C (incubator Thermomixer C, Eppendorf, Germany). Control samples (no AMPs) were suspended and incubated in PBS. Cell suspensions were centrifuged at 1300 × for 4 min in 1.5 ml Eppendorf tubes and resuspended (fixed) in 3% vol/vol glutaraldehyde and 0.1 M cacodylate buffer to bring the rapid cessation of biological activity and to preserve the structure of the cell. After fixation, the samples were washed and postfixed in 1% vol/vol OsO4 in 0.1 M cacodylate buffer. Dehydration was carried out in an ascending ethanol series followed by two steps with propylene oxide (Hayat, 1989) and embedded in Agar Low Viscosity Resin (Agar Scientific, Stansted, UK). Ultrathin sections (70–80 nm) were prepared on an Ultramicrotome UC6 (Leica Microsystems, Vienna, Austria) equipped with a 35° Ultra Diamondknife (Diatome, Nidau, Switzerland). The grids were poststained with Uranyless (Science Services, Munich, Germany) and lead citrate according to Reynolds (Hayat, 1989). TEM images were acquired with Tecnai T12 at 120 kV (TFS, Warmond, Netherlands).
Small-angle scattering
(Ultra-) small-angle X-ray scattering
Request a detailed protocolUSAXS/SAXS measurements were performed on the TRUSAXS beamline (ID02) at the European Synchrotron Research Facility (ESRF), Grenoble, France. The instrument uses a monochromatic beam (λ = 0.0995 nm) that is collimated in a pinhole configuration. Measurements were performed with sample-to-detector distances of 30.8 and 3.0 m, covering a q-range of 0.001–2.5 nm–1 (Narayanan et al., 2018). Two-dimensional scattering patterns were acquired on a Rayonix MX170 detector, normalized to absolute scale, and azimuthally averaged to obtain the corresponding one-dimensional USAXS/SAXS profiles. The normalized cumulative background from the buffer, sample cell, and instrument were subtracted to obtain the final . Bacterial samples ( CFU/ml) were incubated with peptides in buffer for 1 hr at 37°C and directly measured in a quartz capillaries of 2 mm diameter (37°C), mounted on a flow-through setup in order to maximize the precision of the background subtraction. Time-resolved experiments were instead performed with a stopped-flow rapid mixing device (SFM-3/4 Biologic, Seyssinet-Pariset, France), with 50 ms mixing of bacterial and peptides stock suspensions (37°C), and enabling data acquisition after a kinetic time of about 2.5 ms (Narayanan et al., 2014). A total of 50 frames was recorded for each experiment with an exposure time of 0.05 s and a logarithmic time-spacing ranging from 17.5 ms to about 10 min. Radiation damage tests were performed on reference systems prior to setting this X-ray exposure times. The scattering intensities were further corrected for sedimentation and background scattering from the stopped-flow cell.
Contrast-variation small-angle neutron scattering
Request a detailed protocolSANS experiments were performed on the D11 instrument at the Institut Laue-Langevin (ILL), Grenoble, France, with a multiwire 3He detector of 256 × 256 pixels (3.75 × 3.75 mm2). Four different setups (sample-to-detector distances of 2, 8, 20.5, and 39 m with corresponding collimations of 5.5, 8, 20.5, and 40.5 m), at a wavelength λ = 0.56 nm (Δλ/λ = 9%), covered a q-range of 0.014–3 nm–1. Two distinct E. coli suspensions were incubated with peptides LF11-215 and LF11-324 in buffer for 1 hr at 37°C. The bacterial concentration during the incubation was 109 CFU/ml, and both peptide concentrations were in the range of the measured MICs. Samples were then washed and resuspended in five different PBS solutions, containing 10, 30, 40, 50, or 90 wt% D2O. Samples ( CFU/ml) were contained in quartz Hellma 120-QS banjo-shaped cuvettes of 2 mm pathway and measured at 37°C. Cuvettes were mounted on a rotating sample holder, which prevented the bacteria from sedimenting. Data were reduced with the Lamp program from ILL, performing flat field, solid angle, dead time, and transmission correction. Further data were normalized by the incident neutron flux (via a monitor) and corrected by the contribution from an empty cell. Experimental setup information and data are available at https://doi.ill.fr/10.5291/ILL-DATA.8-03-910 .
Note that the present experimental time (~2 hr) is much shorter than the onset of cell lysis (Zweytick et al., 2011). As control, SANS measurements were repeated at extended times after mixing with the AMP in order to test for sample stability (data not shown), in terms of shape, cell-wall structure, and densities. The scattering intensities after 2, 4, 6, and 8 hr were comparable (with the exception of a weak decrease of between 2 and 4 hr). Further, a comparison with TEM and SAXS data suggests that the peptide-induced cell damage does not progress any further after 1 hr.
Data analysis: Peptide clusters
Request a detailed protocolReference O-LF11-215 samples were measured in the MIC range to investigate the microstructure of the peptide clusters. The SAXS pattern of O-LF11-215 was fitted with the equation
where the term in brackets is related to the structure of the aggregates, and the exponential Guinier function accounts for the form factor of their subunits of radius of gyration (Zemb and Lindner, 2002), and is forward scattering intensity. The function is the Porod law that describes the high-q asymptotic trend of scattering signal from the aggregates (Glatter et al., 1982), where is a scaling factor that depends on the surface properties of the aggregates, and is related to their fractal dimension (Sorensen, 2001; see Appendix 1).
Data analysis: Bacterial modeling
Request a detailed protocolX-ray and neutron scattering data were jointly analyzed with a recently reported analytical scattering model (Semeraro et al., 2021b). USAXS/SAXS patterns of end states displayed an excess scattering contribution between q ∼ 0.1–0.2 nm–1 in the case of LF11-215 and LF11-324, not visible in the corresponding SANS patterns. Note that SANS experiments were conducted on samples that were washed and resuspended in different D2O-containing buffer, while SAXS data were acquired immediately after 1 hr incubation with peptides. Together with the observation of OMVs by TEM (Figure 1), this suggests that the additional scattering contribution in SAXS data could be due to freely diffusing OMVs in the suspension medium.
All scattering data were fitted with the analytical functions
where is the scattering form factor for E. coli, as reported in Semeraro et al., 2021b, and is the form factor of a sphere of radius , with j1 being the normalized spherical Bessel function of order 1. The prefactor is the OMV forward scattering, where , , and are, respectively, their number, volume, and SLD difference to the buffer. The choice of a simple spherical form factor was driven by its simplicity for checking whether bacteria and OMVs have to be considered as noninteracting scatterers or not. Tests using an interaction cross-term approximating budding OMVs by spheres decorating a larger surface (Larson-Smith et al., 2010) did not result in significant contributions, confirming the dominance of freely diffusing OMVs. Note also that the values were independent of the shape of the normalized form factor as they include our estimation of and (Appendix 2). Finally, in the case of SANS data, instrumental smearing was taken into account. Data were fitted with a convolution of and a Gaussian function with q-dependent width values, as provided by the reduction tools at D11. In USAXS/SAXS data, the smearing effect was negligible.
After thorough testing, the analysis of SAS data (end states and kinetics) was conducted using only seven adjustable parameters describing . These were the number of LPS molecules, ; the cytoplasm radius, ; the SLDs of the cytoplasmic, , and periplamic space, ; the periplasmic average thickness, , and its deviation, ; and the SLD of the peptidoglycan layer, . Additionally, and were fitted for scattering intensities in the presence of LF11-215 and LF11-324, while was used and adjusted in the case of O-LF11-215. Other parameters of were fixed according to the O-LF11-215 alone systems (see Table 1). This allowed us to fully describe the scattering-pattern variations upon peptide activity. Other parameters, including those accounting for the structure of inner and outer membranes, were fixed to the references values (see complete description in Semeraro et al., 2021b; all fixed parameters are listed in Tables 2 and 3).
The scattering intensities of O-LF11-215-aggregates were comparable to that of bacteria in the high q-range (Figure 2—figure supplement 1D). While this affected the quality of the ultrastructural parameters, it enabled at the same time the investigation of the kinetics of the AMP uptake. Indeed, by assuming that O-LF11-215 is primarily forming aggregates in solution, at ms can be converted to the total known peptide concentration . Hence, the further assumption that peptides leaving the clusters are directly partitioning into the cell allows to convert the difference to . It follows that
where is the number of O-LF11-215 partitioned within the volume of a single cell that can be obtained time-resolved USAXS/SASX data.
Finally, time-resolved USAXS/SAXS data were fitted using the parameters of the initial (see Semeraro et al., 2021b) and end states as boundaries and guide to refine the minimization (see Figure 2—figure supplement 4). This was accomplished by means of a genetic selection algorithm exploiting >300 repetitions of converging fittings (see details in Semeraro et al., 2021b). Fit results were filtered using an ad hoc threshold of , which enabled the exclusion of local minima. Filtered results were used to build up parameter distributions that in turn provided mean values, associated errors (standard deviation), and correlation coefficients. Exemplary distributions are shown in Figure 2—figure supplement 5 and Figure 2—figure supplement 6. Variations in and at ms are due to lower signal-to-noise ratio available in time-resolved measurements.
Appendix 1
Clusters of acylated peptide O-LF11-215
Peptide clusters formed by O-LF11-215 were investigated by Trp fluorescence and USAXS/SAXS. Trp spectra displayed a nm and an FWHM of about 67 nm, which can be related to a heterogeneous distribution of Trp with apolar surroundings (Ladokhin et al., 2000). In addition, O-LF11-215 exhibited Trp phosphorescence emission at about 450 nm, which usually is not measurable due to its dynamic quenching by oxygen and impurities in aqueous suspensions (Cioni and Strambini, 2002). Its presence suggests that a significant portion of Trp residues are buried in hydrophobic cores, with no access to the solvent and with a local high viscosity.
USAXS/SAXS data in the low q-range (qmin ~ 0.005 nm–1) exhibited a featureless decay of intensity with a slope of . This slope value is typical for mass fractals, that is, highly branched objects with high surface-to-volume ratio, while suggests a minimum aggregate size of . Furthermore, a Guinier term is needed to fit the shoulder at about nm–1 corresponding to an average radius of gyration nm. Note that this feature also does not vanish for different choices of scaling constants for background subtraction. Interestingly, this value is way too high to describe O-LF11-215 monomers, whose expected radius of gyration would be lt1 nm. Possibly, peptide monomers create smaller aggregates of mean size nm, which in turn form a heterogeneous and branched supramolecular structure with the characteristics of a mass fractal.
Appendix 2
Estimating the scattering contribution of OMVs
The prefactor of the scattering contribution from extracellular, independent objects used in Equation 6 is , where is the number of OMVs, is the volume of an OMV, and is average SLD difference to the buffer. This definition of forward scattering (per single cell) is valid for every noninteracting object. Hence, to validate the assumption that this scattering contribution is related to OMVs, it is interesting to calculate possible , , and values. Note that even if modeling OMVs as homogeneous spheres appears as a crude first-order approximation, can be associated to its radius (within ∼10% confidence). Assuming the same lipid asymmetry as in the outer membrane, the inner leaflet of OMVs can be mimicked by a 3:1 mole mixture of palmitoyl-oleoyl-phosphatidylethanolamine (POPE) and palmitoyl-oleoyl-phosphatidylglygerol (POPG), respectively (De Siervo, 1969; Lohner et al., 2008; Leber et al., 2018). The SLD membrane profiles of these lipids have been thoroughly investigated (Kučerka et al., 2012, Kučerka et al., 2015). The outer leaflet might instead be dominated by LPS, whose lipid A possesses about six short C14:0 chains (Kim et al., 2016), and the polar region can be approximated as two PG units, in terms of molecular volume and SLD. In addition, LPS inner and outer core volumes and SLDs can be calculated from Heinrichs et al., 1998 and Müller-Loennies et al., 2003, neglecting O-antigen chains for simplicity. Gathering all this information, similarly to the membrane structure estimation in Semeraro et al., 2021b, the vesicles would have a membrane thickness of 4.1 nm and an average SLD of 1.1 × 10–3 nm–2 (volume-weighted averages). The lumen of OMVs can be quite diversely composed (Beveridge, 1999). We tentatively assigned the SLD of the periplasmic space of the end-state system, that is, 9.68 × 10–4 nm–2.
Together with a buffer SLD of 9.47 × 10–4 nm–2, a measured 1500 ± 200 nm and a 15.4 ± 0.6 nm, for instance, then leads to the estimate = 1200 ± 400 and a total lipid surface (inner and outer leaflets of all OMVs) of (4 ± 2 )× 106 nm2.
Appendix 3
Statistics of bacterial inhibition
Assuming that the AMP-induced delayed bacterial growth is entirely due to a lower number density of survived cells (Marx et al., 2021b), the inhibited fraction of cells, , as a function of peptide and cell concentrations was fitted with a heuristic approach. Specifically, we used the sigmoidal Gompertz function
where [P] is the total peptide molar concentration, and b and c and are related, respectively, to the position and width of the sigmoidal. Physically this can be interpreted analogously to a transition between two states; alive and dead bacteria. The width of this transition is associated with its cooperativity, that is, the sharper the transition the higher the cooperativity. This does, however, not imply any cooperative peptide–peptide or peptide–bacteria interaction. It is a mere measure of efficiency to inhibit bacterial growth that will be influenced by several stochastic processes, such as cell-to-cell variations of the average number of partitioned peptides.
Mathematically, can be associated to a CDF, that is, it measures the probability of finding a certain number (or fraction) of inhibited cells at a given peptide concentration. This allows to derive the PDF by calculating the derivative, , as well as the inverse CDF, , that is, the -quantiles. This maps values to the inhibitory peptide concentrations ICx, where is the corresponding inhibited bacterial percentage; by definition, . In addition, the set of and values as a function of can be interpolated to obtain, for example, a continuous trend of ICx as a function of .
Note that, by definition, the derived PDFs express the probability of having a ‘killing’ event at a given [P].
Appendix 4
Tryptophan fluorescence
The native fluorescence of the single Trp residue present in LF11 peptides was exploited to validate the partitioning investigation through the equi-activity analysis. For every system, emission signals from partitioned peptides exhibited a weak blue shift, with values in the range 346–354 nm for LF11-215, and 340–350 nm for LF11-324. values showed no significant variations, instead. Interestingly, these values are consistent with a scenario in which a significant amount of partitioned peptides are heterogeneously dispersed in a polar environment and in a configuration allowing full dynamics of the Trp residues (Burstein and Emelyanenko, 1996; Ladokhin et al., 2000). This is consistent with the values calculated via ζ-potential measurements, suggesting that a relevant portion of partitioned AMPs are still in monomeric state in the intracellular space.
Given that , values extracted via Trp fluorescence were compared with those obtained from the antimicrobial activity assays at [P] = MIC and MIC × 0.5 (LF11-215; data not shown). These two independent methods gave comparable , confirming the validity of values.
Data availability
The current manuscript is a biophysical study, reporting data analysis of scattering data and peptide partitioning assays in vitro. All relevant data are included and plotted in the manuscript. In addition SANS raw data are accessible (http://dx.doi.org/10.5291/ILL-DATA.8-03-910). The modelling code for data analysis consists of a standard chi-squared minimization algorithm. The implemented analytical functions are described in the Methods and Materials section of the manuscript.
-
Institut Laue-LangevinAntimicrobial peptide induced phase separation in E. coli membrane mimetic systems.https://doi.org/10.5291/ILL-DATA.8-03-910
References
-
Structures of gram-negative cell walls and their derived membrane vesiclesJournal of Bacteriology 181:4725–4733.https://doi.org/10.1128/JB.181.16.4725-4733.1999
-
Log-normal description of fluorescence spectra of organic fluorophoresPhotochemistry and Photobiology 64:316–320.https://doi.org/10.1111/j.1751-1097.1996.tb02464.x
-
Lights, camera, action! Antimicrobial peptide mechanisms imaged in space and timeTrends in Microbiology 24:111–122.https://doi.org/10.1016/j.tim.2015.11.004
-
Tryptophan phosphorescence and pressure effects on protein structureBiochimica et Biophysica Acta 1595:116–130.https://doi.org/10.1016/s0167-4838(01)00339-9
-
BookPrinciples and Techniques of Electron Microscopy - Biological ApplicationsCambridge University Press.https://doi.org/10.1007/978-1-349-09857-6
-
Antimicrobial susceptibility testing: a review of general principles and contemporary practicesClinical Infectious Diseases 49:1749–1755.https://doi.org/10.1086/647952
-
Mechanism matters: A taxonomy of cell penetrating peptidesTrends in Biochemical Sciences 40:749–764.https://doi.org/10.1016/j.tibs.2015.10.004
-
Bilayer properties of lipid a from various gram-negative bacteriaBiophysical Journal 111:1750–1760.https://doi.org/10.1016/j.bpj.2016.09.001
-
Molecular structures of fluid phosphatidylethanolamine bilayers obtained from simulation-to-experiment comparisons and experimental scattering density profilesThe Journal of Physical Chemistry. B 119:1947–1956.https://doi.org/10.1021/jp511159q
-
How to measure and analyze tryptophan fluorescence in membranes properly, and why bother?Analytical Biochemistry 285:235–245.https://doi.org/10.1006/abio.2000.4773
-
Small angle scattering model for Pickering emulsions and raspberry particlesJournal of Colloid and Interface Science 343:36–41.https://doi.org/10.1016/j.jcis.2009.11.033
-
Intracellular targeting mechanisms by antimicrobial peptidesAntimicrobial Agents and Chemotherapy 61:1–16.https://doi.org/10.1128/AAC.02340-16
-
Synergism of antimicrobial frog peptides couples to membrane intrinsic curvature strainBiophysical Journal 114:1945–1954.https://doi.org/10.1016/j.bpj.2018.03.006
-
BookDevelopment of Novel Antimicrobial Agents: Emerging Strategies (1st edition)Norfolk, UK: Horizon Scientific Press.
-
BookLiposome-Based Biomembrane Mimetic Systems: Implications for Lipid–Peptide InteractionsIn: Leitmannova Liu A, editors. Advances in Planar Lipid Bilayers and Liposomes. Academic Press. pp. 103–137.https://doi.org/10.1016/S1554-4516(07)06005-X
-
Membrane-active antimicrobial peptides as template structures for novel antibiotic agentsCurrent Topics in Medicinal Chemistry 17:508–519.https://doi.org/10.2174/1568026616666160713122404
-
Experimental concepts for linking the biological activities of antimicrobial peptides to their molecular modes of actionBiochimica et Biophysica Acta. Biomembranes 1862:183275.https://doi.org/10.1016/j.bbamem.2020.183275
-
Antimicrobial peptide activity in asymmetric bacterial membrane mimicsFaraday Discussions 232:435–447.https://doi.org/10.1039/D1FD00039J
-
Bridging the antimicrobial activity of two lactoferricin derivatives in E. coli and lipid-only membranesFrontiers in Medical Technology 3:625975.https://doi.org/10.3389/fmedt.2021.625975
-
BookProbing the self-assembly of unilamellar vesicles using time-resolved saxsIn: Narayanan T, editors. In Advances in Planar Lipid Bilayers and Liposomes. Elsevier. pp. 171–196.https://doi.org/10.1016/B978-0-12-418698-9.00007-1
-
A multipurpose instrument for time-resolved ultra-small-angle and coherent X-ray scatteringJournal of Applied Crystallography 51:1511–1524.https://doi.org/10.1107/S1600576718012748
-
In situ monitoring of bacteria under antimicrobial stress using 31P solid-state NMRInternational Journal of Molecular Sciences 20:E181.https://doi.org/10.3390/ijms20010181
-
How many antimicrobial peptide molecules kill a bacterium? The case of PMAP-23ACS Chemical Biology 9:2003–2007.https://doi.org/10.1021/cb500426r
-
Binding of an antimicrobial peptide to bacterial cells: Interaction with different species, strains and cellular componentsBiochimica et Biophysica Acta. Biomembranes 1862:183291.https://doi.org/10.1016/j.bbamem.2020.183291
-
Non-membrane permeabilizing modes of action of antimicrobial peptides on bacteriaCurrent Topics in Medicinal Chemistry 16:76–88.https://doi.org/10.2174/1568026615666150703121009
-
Evolution of the analytical scattering model of live Escherichia coliJournal of Applied Crystallography 54:473–485.https://doi.org/10.1107/S1600576721000169
-
Light scattering by fractal aggregates: A reviewAerosol Science and Technology 35:648–687.https://doi.org/10.1080/02786820117868
-
Antimicrobial peptides: successes, challenges and unanswered questionsThe Journal of Membrane Biology 239:27–34.https://doi.org/10.1007/s00232-011-9343-0
-
BookNeutrons, X-Rays and Light: Scattering Methods Applied to Soft Condensed MatterNorth-Holland delta series. Elsevier.
-
Expression, purification and structural studies of a short antimicrobial peptideBiochimica et Biophysica Acta 1788:314–323.https://doi.org/10.1016/j.bbamem.2008.10.015
-
Studies on lactoferricin-derived Escherichia coli membrane-active peptides reveal differences in the mechanism of N-acylated versus nonacylated peptidesThe Journal of Biological Chemistry 286:21266–21276.https://doi.org/10.1074/jbc.M110.195412
Article and author information
Author details
Funding
Austrian Science Fund (P 30921)
- Karl Lohner
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
ESRF – The European Synchrotron and the Institut Laue–Langevin (ILL) are acknowledged for provision of SAXS (proposals LS-2513 and LS-2869) and SANS (exp. 8-03-910) beamtimes. The authors are grateful to T Narayanan for his invaluable support, and thank JM Devos, D Marquardt, and M Pachler for their support during the proof-of-concept experiments (LS-2513), and the biological support laboratory at EMBL Grenoble for providing access to the laboratory equipment for bacterial sample preparation. The authors also acknowledge N Malanovic for sharing her expertise about bacterial cultures and S Keller for the fruitful discussions. Finally, the authors thank the whole staff of ID02 and the D11 for support and availability.
Copyright
© 2022, Semeraro et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 951
- views
-
- 266
- downloads
-
- 14
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Computational and Systems Biology
- Physics of Living Systems
Explaining biodiversity is a fundamental issue in ecology. A long-standing puzzle lies in the paradox of the plankton: many species of plankton feeding on a limited variety of resources coexist, apparently flouting the competitive exclusion principle (CEP), which holds that the number of predator (consumer) species cannot exceed that of the resources at a steady state. Here, we present a mechanistic model and demonstrate that intraspecific interference among the consumers enables a plethora of consumer species to coexist at constant population densities with only one or a handful of resource species. This facilitated biodiversity is resistant to stochasticity, either with the stochastic simulation algorithm or individual-based modeling. Our model naturally explains the classical experiments that invalidate the CEP, quantitatively illustrates the universal S-shaped pattern of the rank-abundance curves across a wide range of ecological communities, and can be broadly used to resolve the mystery of biodiversity in many natural ecosystems.
-
- Computational and Systems Biology
- Physics of Living Systems
Planar cell polarity (PCP) – tissue-scale alignment of the direction of asymmetric localization of proteins at the cell-cell interface – is essential for embryonic development and physiological functions. Abnormalities in PCP can result in developmental imperfections, including neural tube closure defects and misaligned hair follicles. Decoding the mechanisms responsible for PCP establishment and maintenance remains a fundamental open question. While the roles of various molecules – broadly classified into “global” and “local” modules – have been well-studied, their necessity and sufficiency in explaining PCP and connecting their perturbations to experimentally observed patterns have not been examined. Here, we develop a minimal model that captures the proposed features of PCP establishment – a global tissue-level gradient and local asymmetric distribution of protein complexes. The proposed model suggests that while polarity can emerge without a gradient, the gradient not only acts as a global cue but also increases the robustness of PCP against stochastic perturbations. We also recapitulated and quantified the experimentally observed features of swirling patterns and domineering non-autonomy, using only three free model parameters - the rate of protein binding to membrane, the concentration of PCP proteins, and the gradient steepness. We explain how self-stabilizing asymmetric protein localizations in the presence of tissue-level gradient can lead to robust PCP patterns and reveal minimal design principles for a polarized system.