Open Access Article
Mohammad A. Haquea,
Riku Imamurab,
George A. Browna,
Venkata R. Krishnamurthia,
Isabelle I. Niyonshutib,
Tiffany Marcellea,
Leanne E. Mathurinb,
Jingyi Chen
*bcd and
Yong Wang
*acd
aDepartment of Physics, University of Arkansas, Fayetteville, Arkansas 72701, USA. E-mail: yongwang@uark.edu
bDepartment of Chemistry and Biochemistry, University of Arkansas, Fayetteville, Arkansas 72701, USA. E-mail: chenj@uark.edu
cMicroelectronics and Photonics Graduate Program, University of Arkansas, Fayetteville, Arkansas 72701, USA
dCell and Molecular Biology Program, University of Arkansas, Fayetteville, Arkansas 72701, USA
First published on 12th December 2017
Silver nanoparticles (AgNPs) are well known to exhibit antimicrobial effects through plausible interactions with proteins and/or deoxyribonucleic acid inside the bacteria. Yet a quantitative understanding on the antimicrobial activities of AgNPs remains obscure. Here we conducted in-depth kinetic growth assays and colony-forming unit (CFU) assays on Escherichia coli (E. coli) cultured in AgNPs or Ag ion-containing growth media. Compared to the Ag-absent culture medium, it was found that the growth rate of the bacteria remained unaffected but the lag time of the bacterial growth was extended due to the presence of AgNPs or Ag ions. From the CFU-based time-kill curves, we observed that fractions of E. coli were killed exponentially in the presence of AgNPs or Ag ions. Based on the experimental data, a quantitative model was established to describe the antimicrobial activity of AgNPs. The predictions from this model agree well with the experimental results. We also showed that the parameters in our model as well as their dependence on the concentrations of Ag and bacteria could, in turn, be determined experimentally. It is expected that our quantitative model and the associated parameters provide an alternative means to minimum inhibitory concentration values for characterizing the antimicrobial activities of Ag ions and AgNPs.
Much effort has been made and various methods have been exploited to investigate how the metal nanoparticles suppress the growth of, or kill, bacteria.3,5,7,20,23,24 It has been found that the shape, size, and surface modifications of metal nanoparticles play important roles in their antimicrobial activities.2,12,15,23,24 Descriptive models have been proposed to explain how metal nanoparticles suppress the growth of bacteria.7,11,22,24 Despite the exciting progresses, little quantitative understanding on the antimicrobial mechanism of metal nanoparticles exists, leaving many basic questions unclear and thus hindering further antimicrobial applications of these nanoparticles. In addition, difficulties currently exist for quantifying the antibiotic activities of metal ions and nanoparticles.
In this work, we focused on developing a quantitative understanding on the antimicrobial activities of silver nanoparticles (AgNPs) and silver ions (Ag ions) using kinetic growth assays and colony-forming-unit (CFU) assays in combination with quantitative modelling. We investigated the kinetic growth curves of Escherichia coli (E. coli) with AgNPs or Ag ions in the growth media, which showed that the growth rate of the bacteria was not affected but the lag time of the bacterial growth was elongated in the presence of AgNPs or Ag ions. In addition, CFU assays and time-kill curves were used to monitor the killing kinetics of bacteria in the presence of AgNPs or Ag ions. Based on these experimental observations, we developed a quantitative model to describe the antibacterial activity of AgNPs. Our model explains the experimental results very well. Furthermore, we quantified the dependence of the parameters in our model on the concentrations of Ag and bacteria experimentally. As an alternative to the minimum inhibitory concentration (MIC) values, our model and its parameters allow us to characterize the antimicrobial activities of Ag ions and AgNPs quantitatively.
000, Sigma-Aldrich), and 4 mL of 282 mM silver trifluoroacetate (AgTFA, Alfa Aesar) were sequentially added to the reaction flask. Once the LSPR peak reached ∼430 nm (∼35 min after addition of AgTFA), the reaction was quenched in an ice bath. Upon cooling, the product was collected by adding acetone to the reaction solution at a ratio of 5
:
1 and centrifuging at 6000 rcf for 10 min. The resulting pellet was purified twice with H2O, collected by centrifugation at 20
000 rcf for 10 min, and resuspended in 10 mL of H2O for future use.For the growth-curve measurements with varying Ag concentrations, the initial OD600 of E. coli was set to 0.05, while the concentrations of AgNPs were 0 (negative control), 1, 5, 10, 20, 30, 40 μg mL−1, or the concentrations of AgNO3 were 0 (negative control), 1, 2, 5, 7, 10 μg mL−1 (the corresponding concentrations of Ag ions were then 0.63 × [AgNO3]). For the experiments with varying initial OD600 values of bacteria, the concentrations of AgNPs and AgNO3 were fixed at 20 μg mL−1 and 5 μg mL−1, respectively, while the initial OD600 values of E. coli were 0.001, 0.005, 0.01, 0.05, 0.1, and 0.2.
000 rcf for another 15 minutes. From the second centrifugation, the supernatants were collected and used for measurement of the concentration of Ag ions by inductively coupled plasma mass spectrometry (ICP-MS, Thermal Scientific). Calibration curves for ICP-MS were obtained by AgNO3 solutions (0–100 μg L−1) in 2% (v/v) nitric acid.
, can be used to describe the growth of bacteria, where k is the growth rate.38 With unlimited resources, k is a constant and the growth of the bacteria is exponential.38 However, under our experimental conditions where growth medium is limited, the growth rate k depends the cell number or cell density,
, where k0 is the initial growth rate, and A is the asymptotic cell density.38 The equation of the bacterial growth becomes
, resulting in a logistic growth curve,
![]() | (1) |
At t = 0, the cell number starts with n0, while as t → ∞, the cell number approaches the asymptotic value A.38 Note that the cell number is equivalent to the cell density or approximately the OD600 readings,39–41 and thus the OD600 values of a bacterial culture and the cell density follow the same equation (eqn (1)).
To account for the antimicrobial activity of AgNPs and Ag ions, we modelled that E. coli were knocked into a suppressed state, S, upon the addition of Ag (in the form of AgNPs or Ag ions). The treated bacteria went through two possible pathways: some bacteria are killed (into the dead state, D), while some bacteria adapt to the new environment, wake up and go back to the active state (A). After waking-up, the active bacteria then undergo the normal cycles of growth and reproduction. We refer to this model as the Suppressed-Active-Dead (SAD) model. The general idea of our model is that bacteria might develop adaptation to the new Ag-containing environment. Similar ideas of acclimation have been developed to model the lag phase of bacterial growth.42–44 In addition, a previously reported active but nonculturable (ABNC) state45 provides a possible biological account for the S state in our model.
Next, conversions between the three states (S → A and S → D) were introduced into the model. If na, ns, and nd are the numbers of cells in the active, suppressed, and dead states, respectively, then the SAD model can be quantitatively described by the following differential equations
![]() | (2) |
It was observed that, after 12 hours of suppression, E. coli with AgNPs at 40 μg mL−1 grew again, which indicates that the commonly defined MIC for the E. coli under the described experimental conditions is above 40 μg mL−1, consistent with some of the previously reported MIC values.13,53,54 The OD600 value of this sample reached its asymptotic value before 30 hours (Fig. 2B), while the asymptotic value was similar to that of the wild type and the samples with AgNPs at lower concentrations (Fig. 2B). Simultaneous control experiments ruled out the possibility of any contamination. Therefore, it is likely that the suppression on the growth of E. coli by AgNPs is only transient under these experimental conditions. In addition, we explored whether the transiency of the antibacterial effect of AgNPs is unique, by measuring the kinetic growth of the same bacteria with Ag ions in the growth medium (provided from AgNO3). Silver ions have been reported to cause DNA condensation and inactivate certain proteins in bacteria.45,55 We performed the same growth-curve measurements on the same bacteria with AgNO3 at various concentrations (0, 1, 2, 5, 7, 10 μg mL−1). The antibacterial effects of Ag ions were detected in our experiments (Fig. 2C). More importantly, we found that, similar to AgNPs, the suppression of bacterial growth by the Ag ions is also transient with AgNO3 in the media at concentrations up to 10 μg mL−1 (Fig. 2D): the bacteria that were completely suppressed in the first 12 hours, but they eventually grew at later time points (Fig. 2D). On the other hand, at higher concentrations of Ag ions (≥30 μg mL−1 AgNO3), no growth of E. coli was observed for up to 10 days.
Interestingly, we found that the suppressing effect of AgNPs and Ag ions on bacteria is not through reducing the maximum specific growth rate of the bacteria. Instead, they elongated the lag time of the bacterial growth significantly. This can be seen directly from the growth curves of bacteria shown in Fig. 2. For example, in the presence of 20 μg mL−1 AgNPs, the suppression on the bacterial growth is obvious; however, the growth rate of the bacteria is very similar to that of the wild type (0 μg mL−1) or with lower concentrations of AgNPs (1, 5 and 10 μg mL−1) as shown in Fig. 2A and 2B. Similar results were observed for Ag ions: although 5 μg mL−1 Ag ions delayed the growth of the bacteria for at least 4 hours, the growth rate of the bacteria was not affected significantly by the Ag ions (Fig. 2C and 2D).
To quantify these observations, we extracted the maximum specific growth rate (μm) and the lag time (λ) from the growth curves by fitting them with the Gompertz model,40,41
. As shown in Fig. 3A and B, the maximum specific growth rates μm of the bacteria were the same (within error) with AgNPs or Ag ions at various concentrations in the growth media. In contrast, the lag time of the bacterial growth λ increased significantly as the concentration of AgNPs or Ag ions in the growth media went up (Fig. 3C and D). In addition, it was found that the elongation of the lag time increased quickly as the concentration of AgNPs or Ag ions goes up: the dose dependence of the lag time can be fitted quite well with quadratic functions (Fig. S3†). Furthermore, it is worthwhile to point out that a previous study on the effect of some antibiotic drugs (ciprofloxacin, tetracycline, etc.) on E. coli showed that the maximum specific growth rate of E. coli decreases as the concentration of the antibiotics increased while the lag time did not change significantly.40 The difference between our results with AgNPs and the previous experiments with antibiotics indicates that the antimicrobial mechanism of AgNPs and Ag ions is possibly not the same as that of those antibiotic drugs.40
It has been controversial regarding the roles of AgNPs vs. Ag ions released from AgNPs in the antimicrobial activity of AgNPs.22,36,50,56 To address this question, we determined the concentration of Ag ions released from dissolution of AgNPs in both water and LB growth medium. We measured that the dissolution percentage of our PVP-coated AgNPs ranges from 0.3% to 1.4% in LB growth medium, and from 0.4% to 2.7% in water (Table S1†). The values are consistent with the range of typical values (<10%) reported previously.36,37,50,57 The low concentration of dissolved Ag ions in the AgNPs solutions, <0.3 μg mL−1 Ag ions (equivalently, <0.5 μg mL−1 AgNO3), is likely to be responsible for why AgNPs showed lower toxicity than Ag ions in our experiments (Fig. 2). On the other hand, based on the growth curves in the presence of AgNO3 (Fig. 2), the concentration of dissolved Ag ions is too low to result in any significant suppression on the bacterial growth by the dissolution of AgNPs alone. Therefore, our data indicates that particle-specific effects are likely to play a role in the toxicity of our PVP-coated AgNPs. This observation supports the conclusions from ref. 36 and 56, but is in contradiction to ref. 22 and 50.
, or n(t) ∼ e−αt, where α is a parameter quantifying the killing rate of the antimicrobial agents.
The killing rate α can be estimated experimentally from the initial portion of the time-kill curves, which was presented in Fig. 4B. Under our experimental conditions where the concentration of AgNPs ranged from 10 to 80 μg mL−1, the killing rate was roughly independent of the concentration of AgNPs (Fig. 4B). Therefore, we combined the data for AgNPs at different concentrations and fitted them with a single curve, ln
n(t)=−αt + c, which yielded α = 0.77 ± 0.10 h−1. The R2 value for the fitting over the combined data is 0.89, indicating that the use of a single killing rate is reasonable (Fig. 4B).
It is also noted that the exponential growth phases (i.e., the increasing portions of the curves) were shifted to the right in the presence of AgNPs (Fig. 4A and S1C†), which indicated the time for preparing the bacteria to grow exponentially was longer in the presence of AgNPs than without AgNPs. In other words, the apparent lag phase was elongated due to the addition of AgNPs, consistent with our observations from the kinetic growth assays (Fig. 2 and 3).
We note that the counts of colonies in CFU assays dropped down significantly in the time-kill curves when the bacteria were treated with Ag (Fig. 4), but the OD600 values in the kinetic growth curves did not decrease significantly upon the addition of AgNPs or Ag ions (Fig. 5). The discrepancy is due to the difference on what exactly the two assays measure. In the CFU assays, only the live bacteria appear as colonies that are counted. In contract, OD600 readings do not distinguish whether the bacteria are alive or dead. As long as the bacterial cells are not lysed, they contribute to the measured OD600 values. Therefore, the observed discrepancy indirectly indicates that most of the bacteria killed by AgNPs or Ag ions were not lysed, which is consistent with a model proposed by Zhou et al.24 and other studies using TEM imaging and X-ray microanalysis.3,55
Our developed model, termed as the Suppressed-Active-Dead (SAD) model, was described in the Section of Materials and Methods, and shown in Fig. 6A. We solved the SAD model was solved numerically, and it turned out that the SAD model predicted all the key features in our observations of the aforementioned several types of experiments. For example, several examples of growth curves of E. coli with different wake-up rates β, but with a constant killing rate α, were shown in Fig. 6B. Note that all the bacterial cells (nOD = na + ns + nd) were included for plotting the growth curves, as all of them contributed to the OD600 readings, assuming the dead cells were not lysed. The shape of the growth curves in our model was sigmoid, the same as the logistic function of the “wild type” (eqn (1), blue curve in Fig. 6B). Second, it was observed that the exponential growth phases were basically parallel to that of the “wild type” for all the different wake-up rates, indicating that the maximum specific growth rate μm was neither affected by the introduction of the suppressed state, nor by the wake-up rate β. In contrast, the lag time was extended in the SAD model, and the lag time extension was longer for smaller wake-up rates, which presumably correspond to higher concentrations of Ag.
In addition, our SAD model reproduced the time-kill curves based on CFU assays, as shown in Fig. 6C, where the numbers of live bacterial cells (active or suppressed, nCFU = na + ns) were plotted as a function of time with various wake-up rates. The CFU counts from our SAD model decreased in the first few hours and then went up, exactly the same as our experimental observations (Fig. 4A). Furthermore, it was shown from the model that the turning point moved lower vertically and the effective killing time became longer for smaller “wake-up” rates (corresponding to higher concentrations of silver), both of which are consistent with our experimental time-kill curves (Fig. 4A). Interestingly, although a constant killing rate (α) was used to predict the time-kill curves, the killing rates appeared to be slightly different (orange curve vs. purple one in Fig. 6C). Therefore, although the measured decreasing slope for AgNPs at 10 μg mL−1 (green squares in Fig. 4A) appeared to be different from that for the AgNPs at 40 μg mL−1 (blue triangles in Fig. 4A), it is not necessary to assume different killing rates under our experimental conditions. On the other hand, we point out that, in general, the killing rate is expected to depend on concentration of Ag, as well as the physiochemical properties of AgNPs such as the size, shape and surface modifications.
Moreover, the experimental results of halted growth of E. coli in the exponential phase due to the addition of AgNPs or Ag ions could be predicted from our SAD model. For the prediction, the bacterial population was first calculated according to our SAD model with an initial condition that all the bacteria were in the active state. When the OD600 reached 0.5, Ag was added to the culture virtually by forcing all the bacteria to go into the suppressed state. Then the cell density of the culture (nOD = na + ns + nd) was calculated again as a function of time according to the SAD model. As shown in Fig. 6D, our SAD model accurately predicted the halted growth of the bacteria: with smaller wake-up rates (or equivalently higher concentrations of Ag) the growth of bacteria became slower or halted. Again, the prediction of the SAD model is consistent with our experimental observations (Fig. 5).
Our model also predicts that the bacteria halted in the exponential phase would resume to grow with the same growth rate after they wake up (Fig. S4A†). To verify this prediction, we repeated the experiments in Fig. 5 for a longer time: the growth of bacteria with Ag ions added in the exponential phase (10 μg mL−1 AgNO3) was monitored for 36 hours. It was confirmed that, after a long lag time, the bacterial culture resumed growth with the same specific growth rate μm as the control (Fig. S4B†), consistent with our prediction (Fig. 6D and S4A†).
It is worthwhile to point out that our SAD model encompasses both bacteriostatic and bactericidal effects of antimicrobial agents. Note that although the definitions of “bacteriostatic” and “bactericidal” seem straightforward (i.e., inhibiting vs. killing), the two pure categories of antimicrobial agents do not exist in reality:52 bactericidal drugs usually fail to kill every single microorganism within 24 hours, and most so-called bacteriostatic drugs do kill some bacteria within 24 hours.52 The practical difference between the two categories lies mainly on the quantitative side (i.e., below or above 99.9% killing). In principle, our SAD model can predict both bactericidal and bacteriostatic effects by tuning the two parameters, α and β. If α were high enough, the antimicrobial effect would become bactericidal; if β were high enough, the effect would be bacteriostatic. In addition, it is expected that the SAD model can be extended to more complicated models to account for the existence of persisters, which are a small fraction of quiescent bacteria (i.e., the suppressed state in our SAD model) that survive and regrow.59
As discussed above, although the killing rate α is expected to depend generally on various factors, it can be treated as a constant under our experimental conditions. Therefore, in this study to avoid over-fitting, we assumed a constant killing rate when estimating the wake-up rate β from measured kinetic growth curves (Fig. 2). Another rationale behind this assumption is that changing the killing rate alone predicted that the turning points move towards lower left (Fig. S5†), conflicting with experimental data in which the turning points shifted towards lower right (Fig. 4A). Therefore, we conclude that the wake-up rate is the major effecting parameter under the current experimental conditions.
The idea for estimating the wake-up rate β is to optimize it in the model (eqn (2), with α fixed at 0.77 in the presence of AgNPs) such that the predicted growth curves from the model matches the measured ones. The optimization process is similar to curve fitting through least-square regression, i.e., by minimizing the squared residuals
where yi are measured values and
are fitted values. However, instead of estimating the fitted values
via an explicit function f(xi,β) in least-square curve-fitting, the “fitted” y-values in the optimization here were estimated by numerically solving our model (eqn (2)).
It is possible to determine the dependence of the wake-up rate on the concentration of AgNPs in the growth media by estimating the wake-up rates β for each and every measured growth curve for E. coli in the presence of AgNPs at different concentrations. Note that each growth curve was “fitted” independently, and no prior information about [AgNP] was needed for the “fitting”. However, the independently “fitted” wake-up rates β were clearly dependent on [AgNP] as shown in Fig. 7A. It was found that the logarithm of the wake-up rate was linear to the concentration of AgNPs, ln
β = κs × [AgNP] + θs where κs is the slope and θs is the intercept of the ln
β ∼ [AgNP] plot. Fitting the β ∼ [AgNP] curve resulted in κs = −0.40 ± 0.03 and θs = 4.5±0.7, with a R2 value of 0.98 (Fig. 7C). The negative slope, κs, indicates that AgNPs at higher concentrations would result in longer lag phases, consistent with our experimental measurements (Fig. 3). We note that, in a previous model, based on energy flux, for the toxicity of metal-based nanoparticles, a parameter termed as the target acclimation energy density was responsible to characterize the dynamics of bacterial acclimation to new environments,42,43 and thus relevant to our wake-up rate β. Our results that β depends on the concentration of AgNPs indicate that the target acclimation energy density in the models by Klanjscek et al.42,43 is likely to correlate with exposure to AgNPs.
Similarly, we determined the dependence of the wake-up rate β on the initial concentration of E. coli, OD0600. For this purpose, we performed kinetic growth assays for E. coli with a constant concentration of AgNPs (20 μg mL−1) but different initial concentrations of the bacteria. We observed that the growth curves were shifted to the right at lower initial bacterial concentrations (Fig. S6†). Fitting of the growth curves showed that the maximum specific growth rate didn't change much for different OD0600 (Fig. S6 and S7†), but the lag times increased when the initial OD600 of the bacterial cultures decreased (Fig. S6 and S7†). These experimental results could be reproduced qualitatively using our SAD model (Fig. S8†). More importantly, we obtained the wake-up rates at varying OD0600 values from this set of growth curves as described above. Interestingly, the logarithm of the wake-up rate is again roughly linear, ln
β = κd × OD0600 + θd, where κd = 43 ± 8 (the slope) and θd = −5.3 ± 0.9 (the intercept) were fitted parameters from Fig. 7B (R2 = 0.90). The positive slope, κd, indicates that lower initial bacterial concentrations would result in longer lag phases, consistent with our experimental measurements (Fig. S7†).
Taking together, we obtained that the wake-up rate was approximately exponential to the concentrations,
![]() | (3) |
![]() | (4) |
The dependences of the wake-up rate on the concentration of AgNPs and initial concentration of bacteria in the culture provided us an additional way to verify our SAD model. This is because both dependences are from eqn (4) but they are from two independent sets of kinetic growth curves (varying [AgNP] vs. varying OD0600). By comparing eqn (3) and (4), we have
![]() | (5) |
![]() | (6) |
As the parameters θs and θd were obtained with OD0600 = 0.05 and [AgNP] = 20 μg mL−1, respectively, if our SAD model is correct, we would expect that
| θs − θd = 0.05κd − 20κs | (7) |
Plugging in the measured values, we found that the left hand side of eqn (7) is θs − θd = 9.8 ± 1.1 and the right hand side is 0.05κd − 20κs = 10.1 ± 0.7. Therefore, both sides of eqn (7) are indeed the same (within error), which strongly justifies the SAD model.
We also note that the exponential dependence of β on the concentration of AgNPs in this study was supported by experimental results reported by others previously.11,22 Using the measured relation between β and [AgNP], the SAD model produced a sigmoid relation between the killing percentage of bacteria and the concentration of AgNPs (Fig. S9†), qualitatively consistent with the experimental results by Xiu et al.22
Our model provides a phenomenological description and quantification of the antimicrobial activity of AgNPs and Ag ions. The model by itself does not explain how the bacteria were knocked out into the suppressed state by the AgNPs or Ag ions, or what is the suppressed state on molecular basis, or how the bacteria adapt to AgNPs and “wake up” to multiply again. On one hand, it is expected that our model can be extended by integrating mechanistic modeling of the toxicity of AgNPs and Ag ions. For example, the models by Klanjscek et al.42,43 based on energy flux provide a plausible direction. On the other hand, to answer these mechanistic questions, further experimental, mechanistic investigations at single cell and single molecule levels are required. Several studies in the literature using TEM and X-ray microanalysis showed that bacteria were not lysed by AgNPs or Ag ions;3,24,55 instead, the AgNPs or Ag ions caused condensation of DNA, lowered expression level of proteins, aggregation of proteins, or other physical changes.3,24,55 These observations indirectly support our hypothesis of the suppressed state, and indicate the possibilities of bacterial adaptation to AgNPs in the growth medium and the reactivation of the bacteria. To elucidate the molecular mechanism of AgNPs and other metal nanoparticles on bacteria, recent emerging super-resolution fluorescence imaging on live single cells at the molecular level60–65 might be of great help in the future.
Additionally, we emphasize that the killing and wake-up rates in our model generally depend on the physicochemical properties of AgNPs such as size, shape, and surface modifications, which have been shown to play a role in their antimicrobial activities.2,12,15 We expect that our model is ready to take into account various physicochemical properties, and thus applicable to a broad range of metal nanoparticles and ions. Systematic future investigations will address these aspects.
Footnote |
| † Electronic supplementary information (ESI) available. See DOI: 10.1039/c7ra10495b |
| This journal is © The Royal Society of Chemistry 2017 |