Interfacial rheology for the assessment of potential health effects of inhaled carbon nanomaterials at variable breathing conditions

Lung surface is the first line of contact between inhaled carbon nanomaterials, CNMs, and the organism, so this is the place where pulmonary health effects begin. The paper analyzes the influence of several CNMs (single- and multi-walled nanotubes with various surface area: 90–1,280 m2/g and aspect ratio: 8–3,750) on the surface-active properties of the lung surfactant, LS, model (Survanta). Effects of CNM concentration (0.1–1 mg/ml) and surface oscillation rate were determined using the oscillating drop method at simulated breathing conditions (2–10 s per cycle, 37 °C). Based on the values of apparent elasticity and viscosity of the interfacial region, new parameters: Sε and Sμ were proposed to evaluate potential effect of particles on the LS at various breathing rates. Some of tested CNMs (e.g., COOH- functionalized short nanotubes) significantly influenced the surfactant dynamics, while the other had weaker effects even at high particle concentration. Analysis of changes in Sε and Sμ provides a new way to evaluate of a possible disturbance of the basic functions of LS. The results show that the expected pulmonary effects caused by inhaled CNMs at variable breathing rate depend not only on particle concentration (inhaled dose) but also on their size, structure and surface properties.

Scientific RepoRtS | (2020) 10:14044 | https://doi.org/10.1038/s41598-020-70909-y www.nature.com/scientificreports/ technique was applied to analyze interactions of nanocarbon particles with the model LS at 0.25 Hz as the typical rate of respiration. However, the question arises if the effect of nanoparticle inhalation will remain similar also at variable breathing pattern, which is a quite common situation in the real-life (people working or doing sport exercises, etc. ). To answer that, in this paper we propose to investigate CNMs interactions with the model LS at dynamic conditions that correspond to different breathing rates (2-10 s per breath). We demonstrate the usefulness of the concepts of interfacial rheology in the determination of such effects, and we propose new parameters to assess them quantitatively.
Quantitative analysis and the physiological role of LS dynamics. Several experimental systems have been used to investigate the dynamic surface-active properties of LS in vitro, including Langmuir trough, oscillating bubble tensiometer, constrained bubble/drop tensiometers 22,27,28 . However, not always appropriate conditions have been assured for such studies (e.g., physiological temperature and the surface area variations relevant to breathing cycle). Due to specific experimental conditions needed to determine dynamic properties of LS at the air/liquid interface, the quantitative analysis of the results also requires special measures. Up to now, the following numerical parameters have been most commonly used: (a) the minimum value of the surface tension, σ min (mN/m), recorded during periodical expansion-contraction cycles of simulated breathing.
(b) the amplitude of surface tension variations measured at such conditions (σ max − σ min ). Clements et al. 29 introduced the stability index, SI, that can be derived from the above-mentioned surface tension values: Periodic variations of the surface tension during oscillations of the air/liquid interfacial area show a time-shift which result in the surface tension hysteresis. The normalized area of the hysteresis loop, HA n , was proposed by Notter et al. 30 as the parameter to quantify this feature: The hysteresis arises due to the time-dependent phenomena (relaxation) that may be attributed both to the intrinsic mechanical properties of the air/liquid interface and to the mass exchange of surface-active molecules between the surface layer and the underlying liquid. These dynamics may be described using 2D rheological formalism, i.e. by determining the apparent surface elasticity and the apparent surface viscosity 31,32 . Several authors studied the surface rheology of air/liquid interfaces with different LS models [33][34][35] , however they did not focus on exact relations to the real breathing conditions (in terms of temperature and surface deformation frequency). In the rheological formalism, the departure from the initial value of the surface tension: Δσ = σ − σ 0 is connected with extensional deformation of the interface, γ: where γ denotes the surface deformation rate. Equation (3) is the Kelvin-Voigt model of a visco-elastic air/liquid interface. The temporary surface dilatational deformation (extension) is defined as: where A denotes the area of the interface at the given time instant t, and A 0 -the initial interfacial area (at t = 0). Rheological parameters of the interface are denoted as ε (N/m-dilatational surface elasticity) and μ (s N/m-dilatational surface viscosity), while Δσ E and Δσ L represent the elastic and viscous contribution of the surface tension deviation, respectively. Small-amplitude harmonic deformation of the interface with the angular frequency ω (rad/s) is described as: and the corresponding visco-elastic response takes the form: where φ (rad) denotes the phase-shift (loss angle) that appears due to the viscous properties of the interfacial region. If φ > 0, the system shows the hysteresis. γ m and Δσ m in Eqs. (5) and (6) denote the amplitudes of surface deformation and surface tension, respectively. For purely elastic response of the interface, the loss angle φ is zero, and, if the viscosity predominates, φ approaches π/2. In general, the loss angle is equal: which indicates that surface tension hysteresis depends both on the viscosity-to-elasticity ratio and on the frequency of surface deformation. The hysteresis is larger when viscosity dominates over elasticity, and when oscillations are slower (i.e., ω is larger).
For surfactants that undergo mass exchange (dynamic adsorption and desorption) with the air/liquid interface (as in LS system), the situation becomes more complicated since both ε and μ may be not constant. In such a case, www.nature.com/scientificreports/ these parameters depend on the current surface deformation γ(t) because the composition the interface is changing during surface dilation or contraction due to the mass exchange, i.e. surfactant diffusion and adsorption/ desorption 31 . This process act towards the relaxation of Δσ. Obviously, surface deformation rate γ(t) also affects the amount of the surfactant transferred between the liquid and the interface at a given time period. Therefore, in the system considered in this work, both surface elasticity and surface viscosity should be considered as not the real (intrinsic) rheological characteristics of the interface but rather as the apparent parameters of the interfacial region. A more comprehensive discussion of these issues may be found elsewhere 31,[33][34][35][36][37] . The phenomenon of surface tension hysteresis has been analyzed and discussed in the relation to the physiological functions of LS 13,26,30,[38][39][40] . It has been postulated that the hysteresis is related both to mechanical aspects of breathing (for instance, to the pressure-volume hysteresis in the lungs) and to the mass transfer phenomena on the pulmonary surface, including the hydrodynamic clearance of inhaled deposits from alveoli 26,38 . The rheological analysis of the interfacial dynamics discussed in the current work facilitates the quantitative assessment of these important processes.

Results and discussion
The first results obtained during drop oscillations of LS sample without particles at various frequencies revealed that both surface elasticity, ε, and viscosity, μ, depend linearly on the oscillation period (T = 2π/ω) in the studied surface deformation range, Fig. 1. These relationships can be expressed as: where S ε (N/ms) and S µ (N/m) denote slopes of linear functions. They inform about the sensitivity of ε and μ to changes in the deformation rate corresponding to variable breathing frequency. Both parameters of interception, I ε (N/m) and I µ (s N/m), indicate the extrapolated values of surface elasticity or viscosity at T approaching 0, assuming that Eqs. (8) and (9) remain valid up to the infinite deformation frequency. S ε and S µ can be considered the sensitive numerical indicators of the dynamic surface tension-lowering activity of LS. As seen from Fig. 1, for pure LS (i.e. the control sample -no particles in the system) S µ and S ε are equal 2.95 mN/m and − 2.02 mN/m s, respectively. These numbers show that slower oscillations (i.e., higher T values) correspond to low surface elasticity and increased surface viscosity. Such findings fully agree with the theory since, in general, the viscosity/ dissipation/relaxation are always prevailing at longer time-scale of deformation process 31,34 .
The morphology of studied CNMs is presented on SEM pictures (Fig. 2). Differences in the external diameter of nanotubes MW-1, MW-2, MW-f and SW are clearly seen ( Fig. 2a-d, respectively) and agree with the data listed in Table 2. Both types of nanohorn particles NH and NH-f have similar morphology (note different scale in Fig. 2e,f), however they show various tendency to agglomerate. NH are agglomerated to a greater extent than NH-f, what results in the difference in SSA of these particles (Table 2).
CNMs added to LS cause a shift of both linear relationships ε(T) and μ(T)- Fig. 3. The relationships of surface elasticity and surface viscosity vs deformation period T are altered to a different extent, which clearly indicates a wide spectrum of effects induced by various nanoparticles as a function of concentration. Graph in Fig. 3 also confirms that it is possible to quantify the influence of CNMs on LS interfacial dynamics under variable breathing conditions by comparing the slopes S ε and S µ .  Fig. 4. The window of quasi-constant S ε and S µ values assumed as ± 25% of these values for pure LS (Survanta), is also plotted in this graph (Fig. 3) to indicate clearly the concentration of a given CNM that corresponds to the most substantial alterations in LS dynamic surface activity. This relatively wide window embeds the scatter of data that is typically found in the multi-component system during the measurements of the dynamic surface tension during surface oscillations 32 . It is therefore assumed that results beyond this range may suggest that the dynamic surfactant activity has been altered to the degree that may have consequences for health. After analyzing of all presented data, the potential impact of tested carbon nanoparticles on the LS system is summarized in Table 1.
It is seen that not all tested CNMs induce direct changes in the dynamic activity of the surfactant, even at high nanoparticle concentration (~ 1 mg/ml). For instance, not unequivocal effects were found for multi-walled carbon nanotubes MW-1 and MW-f. Although these nanoparticles are characterized by the similar length (10-30 μm), they differ both in the aspect ratio, AR (1,250-3,750 and 333-1,500, respectively) and the specific surface area, SSA (445 and 190 m 2 /g, respectively). This may suggest that these two parameters are not critical for the effects of carbon nanotubes on LS dynamics at the air/liquid interface. In contrast, data obtained for other CNMs suggest Using the slopes S ε and S μ as indicators of the sensitivity of the apparent rheological parameters of the interface to variations in the oscillation (breathing) rate suggests that strong effects are caused by two types of nanoparticles of the relatively high SSA: NH-f (SSA ≈ 1,300 m 2 /g, AR = 8-25) and SW (SSA ≈ 950 m 2 /g, AR = 1,500) when their concentration exceeds 0.25-0.5 mg/ml. The analogous trend was found for NH particles at similar concentration, although these particles have much lower values of SSA and AR (≈ 400 m 2 /g, and ≤ 25, respectively). The observed effects may be caused by adsorption of LS molecules to the nanoparticles resulting in surfactant depletion, both on the interface and in the aqueous phase 22,41,42 . It should be noted though that the reported SSA values (Table 2) may be not equal to the area accessible for the adsorption of LS molecules. This issue will be discussed later.
Partly hydrophilic COOH-functionalized short nanotubes NH-f (nanohorns) with high SSA and low AR (~ 1,280 m 2 /g and 8-25, respectively), strongly modify dynamic interfacial properties of LS. In contrast to pure LS or the surfactant in the presence of other nanomaterials (where S μ remains either almost constant or increases at high CNM concentration), a sharp decrease in S μ is observed for NH-f particles at their concentrations above 0.5 mg/ml. It suggests a strong disturbance of LS dynamics at variable surface oscillation frequency. Reduction of S μ to almost zero at high concentrations of NH-f particles means that changes in the frequency have no effect on the viscous response of LS interface. Accordingly, the surface tension hysteresis becomes independent on the breathing rate.
It is a completely opposite dynamic response of the air/liquid interface than in the normal ("healthy") situation depicted in Fig. 1, where slower oscillations evidently increase the viscous character of the interface. Therefore, this effect of NH-f particles will decrease the σ-A hysteresis. As stated earlier, large hysteresis of the dynamic surface tension is essential for normal breathing mechanics and mass exchange in the lungs 13,38 , so the above findings should be considered negative for LS functions.
Another interesting observation is that MW-2 nanotubes induce the maximum change in LS dynamics not at the highest (i.e. 1 mg/ml) but at a lower particle concentration (0.5-0.75 mg/ml). A similar effect has been already observed in the earlier studies 22 and it can be attributed to the reduction of effective surface area of CNMs available for LS adsorption by particle aggregation favored by high particle concentration. This behavior of rod-like structures on the air/water interface and in the aqueous phase has been documented by other researchers 43,44 . It may be postulated that the assembling of nanoparticles in the liquid phase and on the air/liquid interface is also enhanced by low AR of particles thanks to reduction of steric restrictions and easier (faster) particle diffusion, Fig. 5. Such an assumption is partly supported by the data of Allegri et al. 45 who were able to explain the suppression of the measured toxic effects by enhanced aggregation of short carbon nanotubes. This mechanism can also explain why the impact of MW-2 (at 0.25-0.75 mg/ml concentration range) is stronger than the effect of MW-1 particles, despite MW-1 have significantly higher SSA. Table 1. Nanoparticle concentration range c (mg/ml) that corresponds to the strongest effects of CNMs on both parameters of the rheological response of LS surface: S ε and S μ .  www.nature.com/scientificreports/ Another explanation of these effects may be proposed after noting that SSA which are reported in Table 2 (in the Methods section) and considered in the ongoing discussion, were determined by nitrogen adsorption method. Such SSA values obviously overestimate the effective surface area of CNMs that is available for adsorption of large LS molecules. For instance, the inner surface of hydrophobic nanotubes with small internal diameter (such as MW-1) is poorly wetted by water, and this may considerably limit penetration of the liquid to the inner volume of MW-1 nanotubes. Moreover, the access of relatively large surfactant molecules or their aggregates to these internal spaces is geometrically restricted. Therefore, the surfaces inside the thinnest nanotubes do not contribute to the mass-exchange and sorption processes that are decisive for the final effects of CNMs on LS dynamics at the air/liquid interface. It is also possible that different CNMs might be dispersed to a different degree when sonicated (see the Methods), which will also cause that their effective surface area is lower than declared. Inter-particle adhesion forces responsible for this effect depend on the CNM surface properties, dimensions and concentration, i.e. on the same factors which are important for aggregation.

Methods
The CNMs studied in this work are characterized in Table 2. They have been carefully selected to cover a possibly broad range of dimensions, structure, specific surface area and surface composition.
The morphology of the particles was studied using SEM (Ultra Plus, Zeiss, Germany) after coating with Cr in Q150T ES sputter coater (Quorum Technologies, UK). The dynamic surface tension of a model lung surfactant (Survanta-Abbott Labs., France) was measured during sinusoidal oscillations of a pendant drop (PAT-1M drop shape tensiometer, Sinterface Technologies, Germany). Survanta is an animal-derived (bovine) mixture that contains 25 mg/ml phospholipids (including 11.0-15.5 mg/ml disaturated phosphatidylcholine), 0.5-1.75 mg/ ml triglycerides, 1.4-3.5 mg/ml free fatty acids, and less than 1.0 mg/ml SP-B and SP-C protein) in sterile 0.9% sodium chloride solution. Mixtures of LS with each type of CNM were prepared by mixing Survanta with presonicated suspensions of nanoparticles in ultra-pure water (ultrasonic energy: 4.8 kJ/ml). This allowed to obtain the mixtures with five CNM concentrations: c = 0.1, 0.25, 0.5, 0.75 and 1 mg/ml. The lowest concentrations correspond to the predicted values after human exposure to highly contaminated dusty environment 40 . The highest values (0.5-1 mg/ml) of nanoparticle concentrations are less realistic, however they allow to shed light on the mechanism of the physicochemical interactions which are less easily detected when CNMs concentration is low. The amount of Survanta and aqueous suspensions of CNMs were adjusted to obtain the same LS concentration in each final sample (2.5 mg/ml), identical to the LS concentration in the control sample without CNMs. The surface tension at the beginning of drop pulsation, i.e. after system equilibration (600 s) was in the range of 35-50 mN/m, depending on CNMs type and concentration. Droplets with the initial interfacial area of 16 mm 2 were pulsated at physiological temperature 36.8 ± 0.2 °C at 20% oscillation amplitude (+ 10% and -10% of

conclusions
The basic physicochemical functions of the lung surfactant are related to the dynamic surface phenomena in the pulmonary alveoli during the breathing cycle. Presented in vitro results confirm the possibility of disturbance of these phenomena by carbon nanomaterials (short or long, single-walled or multi-walled nanotubes) that may be present in the LS system after inhalation at high exposures. The paper introduces new quantitative parameters which help to quantify such effects at conditions that correspond to variable breathing frequency which reflect different levels of physical activity of the body.
It was shown that the influence of CNMs on LS dynamics can be assessed by tracing the changes of the apparent surface rheological parameters of air/liquid interface at different rates of surface deformation. The CNMs studied here alter the surface activity of LS to a different degree under conditions that correspond to variable respiration frequencies. These changes are very specific, i.e. they depend not only on CNM concentration but also on other particle characteristics, and, in general, they are difficult to be predicted a priori. However, the measured phenomena are of high importance in the assessment of a possible impact of inhaled carbon nanomaterials on the respiratory functions in the real-life conditions where changes in the level of physical activity (i.e., in the breathing rate) are common.

Slower diffusion
Faster diffusion The results show that various CNMs differently affect the LS dynamics, and that the effects depend in a distinctive way on the concentration of a given nanomaterial in the system. It has been confirmed that the specific surface area of carbon nanoparticles is one of the key factors since the interaction with LS are partly governed by adsorption of surfactant molecules to the nanomaterials. However, it was also shown that the nominal SSA values cannot be unequivocally correlated with the effect on LS. It is explained by differences in CNMs dispersion in the aqueous phase and by noting that SSA values measured by sorption of small molecules (nitrogen) overstate the CNMs surface area accessible for large surfactant molecules and their aggregates. The discussed effect was observed, among others, for multi-walled carbon nanotubes SW. It was also found that in some cases (i.e., MW-2) nanoparticles probably aggregate stronger above a certain concentration, and this phenomenon further reduces the particle surface area available for LS adsorption. The discussed effect is more important for nanoparticles with low aspect ratio, since they diffuse and aggregate more easily, and form more packed aggregated structures.
The most marked effect of CNMs on LS surface dynamics was observed in the case of functionalized (i.e. hydrophilic) carbon nanohorns with high SSA (MW-f). The presence of such nanoparticles in LS has a unique effect on the parameters S ε and S μ that have been introduced to trace the relationship between ε (surface elasticity) or μ (surface viscosity), and the time-scale of surface deformations, i.e. T (oscillation period). The viscous response of the model alveolar surface is completely suppressed by MW-f nanoparticles at high concentrations, so they reduce the surface tension hysteresis even at a slow deformation (breathing) rates. This suggests that LS will behave abnormally during the respiratory cycle, which may impair the essential pulmonary functions, including ventilation mechanics and the mass transfer.
The results not only indicate that inhaled carbon nanotubes and nanohorns with different surface properties may be responsible for health hazard due to specific physicochemical phenomena on the lung surface, but they also show that these effects may be dependent on respiratory dynamics that is related to different physical activity of the body. Such a dependency of interfacial effects in the LS system on the breathing rate provides a new information required for a more complete analysis of potential risks from inhaled carbon nanoparticles.

Data availability
Raw experimental data available on request. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creat iveco mmons .org/licen ses/by/4.0/.