ABSTRACT
Antimicrobial peptides (AMPs) are ancient and conserved across the tree of life. Their efficacy over evolutionary time has been largely attributed to their mechanisms of killing. Yet, the understanding of their pharmacodynamics both in vivo and in vitro is very limited. This is, however, crucial for applications of AMPs as drugs and also informs the understanding of the action of AMPs in natural immune systems. Here, we selected six different AMPs from different organisms to test their individual and combined effects in vitro. We analyzed their pharmacodynamics based on the Hill function and evaluated the interaction of combinations of two and three AMPs. Interactions of AMPs in our study were mostly synergistic, and three-AMP combinations displayed stronger synergism than two-AMP combinations. This suggests synergism to be a common phenomenon in AMP interaction. Additionally, AMPs displayed a sharp increase in killing within a narrow dose range, contrasting with those of antibiotics. We suggest that our results could lead a way toward better evaluation of AMP application in practice and shed some light on the evolutionary consequences of antimicrobial peptide interactions within the immune system of organisms.
INTRODUCTION
Combinations of drugs can result in three different forms of interactions: synergism, additivity, and antagonism (1–4); i.e., the effect of two drugs combined is stronger, equal, and weaker than that of the individual drug in the equivalent dose, respectively. Combination treatment is supposed to potentially eliminate resistant strains, delay the evolution of drug resistance, reduce the dosage of individual drugs, and hence, diminish side effects (3, 5, 6). A few recent studies, however, report that the success of combination therapy is context dependent, particularly when targeting both sensitive and resistant strains with a combination of drugs of unknown interaction (7–9). These results demonstrated that synergistic drug pairs can efficiently eradicate bacteria but exacerbate selection of resistance, while antagonistic drug pairs showed the reverse trends.
Various methods have been developed to address the efficacy of mostly two-way drug combinations (1, 2). One of the most commonly used approaches in both theoretical and applied research is Loewe additivity (2, 9–11). Here, the effect of two drugs in combination is determined by the sum of ratios of concentrations of drugs in combination divided by concentration of drugs used individually. Note that both the individual drug concentrations and the combined concentrations have the same effect on bacterial growth; we call these concentrations isoeffective concentrations. Theoretically, if the isoeffective concentrations of equivalent effect level achieved in a matrix of gradients of concentrations were connected by line, a concave line represents synergism, while a convex line represents antagonism (2, 12, 13). Recently, a mechanism-free approach was used successfully to predict the outcome of three antibiotics on the interaction between all three possible two-way combinations (14, 15), but the results do not particularly address the question about the nature of interaction (synergism, additivity, or antagonism). How these approaches can be used for a new class of antimicrobials, antimicrobial peptides (AMPs), is basically unknown. Studies on combinatorial effects of antimicrobial peptides, especially within a pharmacodynamics framework, are scarce (16, 17).
Antimicrobial peptides (AMPs), which form an important component of immune defenses in multicellular organisms (18, 19), have been proposed and are being used as new antibiotic drugs. Some AMPs are already commercially available and ready to be applied in clinical practice to replace or accompany conventional antibiotics (20). Additionally, they are supposed to be less likely to induce resistance and mutagenesis in the natural environment, although resistant strains can be obtained under intensive selection in the laboratory (21–23). When AMPs are used in medical applications, they necessarily interact with the patient's own AMPs. Some experimental studies have addressed the effect of individual pairs of AMPs within the context of innate immunity. Coexpressed AMPs on frog skin, PGLa and magainin-2, are synergistic when applied to both Escherichia coli and tumor cells (16). Moreover, AMPs from mammals (17) and insects (24, 25) were shown to synergize. Hence, understanding general principles of AMP interaction will also contribute to our understanding of interactions of AMPs as immune effectors.
Here, we take a pharmacodynamic approach to study the combination effects of AMPs and with combinations of two and three AMPs. Pharmacodynamics capture the functional relationship between drug dosage and bacterial growth or death. We use a modeling approach based on the Hill function (26–28). This model estimates four parameters: MIC, κ, ψmin, and ψmax (Fig. 1A). The minimal concentration at which antibiotic substances can inhibit growth of bacteria is MIC; κ depicts the steepness of the curve relating bacterial growth to drug concentration (Fig. 1B); ψmin and ψmax represent the minimum and maximum growth rates of bacteria, respectively. We studied the pairwise and three-way interactions of AMPs using this pharmacodynamic approach embedded in a Loewe additivity framework. According to the work of Loewe, this was achieved by using either one-half (pairwise) or one-third (three-way) of the concentration of each individual drug (see Materials and Methods) (10). We examined the nature of the interactions, synergism or antagonism, and the concentration dependency of the killing. Using a derivative of the human AMP LL-37 enabled us to study interactions between AMPs expressed by patients' innate immunity and AMPs employed as drugs.
Schematic illustration of four parameters, MIC, ψmax, ψmin, and κ, predicted by the Hill function. The MIC is estimated by the lowest concentration that inhibits the growth of the whole treated bacterium population. ψmax and ψmin represent the maximal and minimal growth rates of bacteria under gradients of drug treatment, respectively. κ predicts the shape and slope of the pharmacodynamic curve; the higher the κ value, the steeper the pharmacodynamic curve.
MATERIALS AND METHODS
Bacteria and media.Escherichia coli MG1655 was grown in Luria-Bertani (LB) broth at 37°C with aeration at 220 rpm in 50-ml tubes. Two hundred microliters of overnight culture was resuspended into 15 ml fresh LB broth, cultured under the same conditions for an additional 2 h, and then used for subsequent assays. Mueller-Hinton (MH) broth was used for the assay of MICs and time-killing curves.
AMPs and antibiotics.We used six different AMPs from different classes of organisms that are commercially available (AnaSpec): cecropin A (Cec) (insect), LL 19-27 (LL) (mammal), melittin (Mel) (insect), pexiganan (Pex) (synthesized AMP, an analog of magainin II; a kind gift of Michael Zasloff), indolicidin (Ind) (mammal), and apidaecin (Api) (insect) (see Table S1 in the supplemental material). These AMPs are effective on either Gram-positive or Gram-negative bacteria (see reviews in references 29 and 30). However, some of these AMPs, e.g., melittin and pexiganan, have anticancer effects (16, 31, 32), which means that they are potentially toxic to human cells, such as erythrocytes. Recently, some low-cell-toxic and serum-stable AMPs have also been under development (33, 34). All these AMPs were dissolved in distilled water with an initial concentration of 1 mg/ml, 5 mg/ml, 10 mg/ml, 1 mg/ml, 1 mg/ml, and 25 mg/ml, respectively, as stock solutions. All antibiotics—ampicillin, ciprofloxacin, gentamicin, kanamycin, neomycin, rifabutin, spectinomycin, and tetracycline—were also dissolved in distilled water and made into 10-mg/ml stock solutions. All the solutions of AMPs and antibiotics were stored at −20°C in a dark environment.
MIC determination.According to a standard protocol (35), stock solutions of AMPs were diluted in MH broth and then diluted in 96-well plates with a 2-fold gradient, that is, from 0.25 μg/ml to 128 μg/ml. All the gradients of antibiotics were from 0.02 μg/ml to 50 μg/ml, except that the gradient of ciprofloxacin was from 0.002 μg/ml to 1 μg/ml. Approximately 5 × 105 log-phase bacteria were added to each well. A positive control containing MH broth and bacteria and a negative control containing only MH broth were included in each plate, and plates were incubated at 37°C overnight.
Measuring killing curves.To estimate killing curves of each AMP and all possible combinations of AMPs, 100× MICs of AMPs were combined as a volume ratio of 1:1 and 1:1:1 in two-AMP combinations and three-AMP combinations; hence, the concentrations of individual drugs are halved or reduced by two-thirds. Thus, Loewe additivity would result in a MIC of any one combination equal to the MIC of the individual drugs (see equations 5 and 6). Thus, 21 two-AMP combinations and 20 three-AMP combinations were generated. An AMP(s) was diluted, starting with 100× MIC, in a 96-well plate to form a 2-fold gradient of concentrations, and 2 × 106 log-phase bacteria were added to a total volume of 100 μl. The plates were incubated at 37°C. Killing was assessed within 1 h, as killing by AMPs is very fast (36, 37). Ten microliters of a mixture of AMPs and bacteria was taken out every 20 min and then immediately diluted in saline solution and plated on the solid agar plates. These solid agar plates were transferred into a 37°C incubator and cultured overnight for CFU determination. The limit of detection in our system is 100 CFUs.
Modeling killing curves.To model the killing curve, the relationship between the concentration of AMP(s) and the killing and/or growth rate of exposed bacteria, we used a Hill function (26):
Determining the effect of combination.Based on the Hill function and isobologram analysis, we obtained the isoeffective concentrations of single drugs and of combinations which achieved a given percentage of their maximal effects or fraction level. The Loewe additivity model defines the additive effect of isoeffective combinations of drugs that result in a certain effect. For example, the combination of drug A and drug B in the isoeffective concentrations, which are CisoA and CisoB, can achieve a level of effect which can also be achieved individually by drug A or drug B with a concentration of CA or CB, respectively. Mathematically, the combination effect of drug A and drug B is defined as follows:
RESULTS
Killing and pharmacodynamic curves.We tested the in vitro effects of single AMPs, two-AMP combinations, and three-AMP combinations on E. coli. All killing curves were obtained by counting viable CFU after treatment (see Fig. S1 in the supplemental material). In most cases, the number of surviving bacteria drastically decreased as a function of time at higher concentrations while slightly increasing at lower concentrations. Killing occurred very quickly at higher concentrations in our system (i.e., bacterial densities below the limit of detection).
The four pharmacodynamic parameters, MIC, κ, ψmax, and ψmin, were estimated by the MCMC method using the generalized linear regression fitted killing rate as a function of concentrations of AMP(s) (Fig. 2 and 3; also see Table S2 in the supplemental material). Notably, all the single AMPs and two- and three-AMP combinations had almost the same ψmin (analysis of variance [ANOVA], ψmin, F1,39 = 1.855, P = 0.181) (Fig. 3; see also Table S2 in the supplemental material). ψmax values were also identical in different treatments as the growth rate of bacteria in low concentrations of AMP(s) was presumably close to the natural growth rate. Two pharmacodynamic parameters, MIC and κ, varied among different treatments, with three-AMP combinations having the lowest MICs and the highest κ values (ANOVA, MIC, F1,39 = 6.647, P = 0.0138; κ, F1,39 = 7.447, P = 0.00935) (Fig. 3; see also Table S2 in the supplemental material). All the treatments showed nearly the same pharmacodynamic trend: a sharp decrease of net bacterial growth with an increasing concentration of AMP(s) as depicted by κ.
Pharmacodynamic curves of AMPs. The pharmacodynamic curves of AMPs were obtained by fitting killing curves to the Hill function (see equation 4). Combinations of two or three AMPs were differentiated. The curves illustrate the effects (reflected as net bacterial growth rate) of increasing the concentrations of AMP(s). The ribbon represents the 95% confidence interval.
Variations of MIC, κ, and ψmin in the Hill function predicted by the MCMC method. Results showed that these parameters vary among combinations with different numbers of AMPs. MICs declined with increasing numbers of AMPs in combination (ANOVA, F1,39 = 6.647, P = 0.0138); combinations with three AMPs had the highest κ values (ANOVA, F1,39 = 7.447, P = 0.00935). ψmin (Psimin) did not show significant differences among single AMPs and two- and three-AMP combinations (ANOVA, F1,39 = 1.855, P = 0.181).
Most AMP combinations are synergistic, but synergy is stronger in three-AMP combinations.To determine the interaction of AMPs, we used Loewe additivity (see equations 5 and 6). The combination index was calculated for concentrations between 5% and 95% of the maximal effect (equation 3). For two-AMP combinations, we found that most of the two-AMP combinations (67%) were completely synergistic (combination indexes were lower than 1) within the effect range, except for the combination of apidaecin and LL 19-27 (ApiLL), which was antagonistic across the whole range; the combinations PexApi and IndApi were antagonistic in low-concentration combinations but synergistic in high-concentration combinations (Fig. 4; also see Fig. S2 in the supplemental material). However, the combinations CecApi and MelApi had a reverse trend, as they were synergistic in lower-concentration combinations and antagonistic in higher-concentration combinations (Fig. 4; see also Fig. S2 in the supplemental material). Eighty-five percent of three-AMP combinations were completely synergistic within the effect range while the combination LLPexApi was completely antagonistic; LLIndApi showed synergistic effects in lower-concentration combinations and antagonistic effects in higher-concentration combinations, but MelIndApi had the reverse trend (Fig. 4; see also Fig. S2 in the supplemental material).
Combination index of AMPs applied at concentrations which can achieve 50% of their maximal effect (E50). At E50, all the combinations with Api (except the combination of Api and Mel) showed antagonistic effects in two-AMP combinations, but only two combinations, ApiIndMel and ApiLLPex, showed antagonistic effects in three-AMP combinations. The gradient of colors represents different levels of each interaction.
Another interesting finding is that three-AMP combinations generally have stronger effects than do two-AMP combinations at a given fraction level within the effect range. The average combination indexes of three-AMP combinations were 30% lower than those of two-AMP combinations (Student's t test, t = 8.2016, df = 606.57, P = 1.42e−15) (Fig. 5). We observed no differences between effects of different fractions for the three-way interactions (ANOVA, F1,661 = 1.332, P = 0.2488) (Fig. 5).
Combination index of fraction level within the effect range. Three-AMP combinations are more significantly synergistic than two-AMP combinations (Student's t test, t = 8.2016, df = 606.57, P = 1.42e−15). The combination index did not vary within different effect ranges in both two-AMP and three-AMP combinations (ANOVA, F1,661 = 1.332, P = 0.2488). Black dots denote outliers.
Relationship between κ values and selection.We compared the κ values of different combinations of AMPs and between AMPs and antibiotics. κ values are higher the more AMPs that are combined. We also found that κ values of AMPs are significantly higher (ANOVA, F1,77 = 150.5, P < 0.001) (Fig. 6) than those of antibiotics, for data obtained in both our laboratory and other laboratories (ANOVA, F1,36 = 1.591, P = 0.215) (Fig. 6).
Comparison of κ values of AMPs and antibiotics in our experiment and other, similar experiments. With similar experimental methods, conditions of measurements, and mathematical models, κ values of antibiotics in our experiment are similar to those of antibiotics in other experiments (ANOVA, F1,36 = 1.591, P = 0.215). However, κ values of AMPs are significantly higher than those of antibiotics both in our experiment and in other experiments (ANOVA, F1,77 = 150.5, P < 0.001). Data in boxes “Antibiotics (1),” “Antibiotics (2),” and “Antibiotics (3)” are from references 26, 40, and 11, respectively.
DISCUSSION
Pharmacodynamic approaches have been frequently applied to conventional antibiotics (2, 11, 26, 40). A good understanding of how antimicrobial peptides eradicate bacteria in complex systems not only relies on the molecular mechanisms of killing but, importantly, necessitates investigation of pharmacodynamics in vitro, as done here and in vivo (C. Zanchi, P. R. Johnston, and J. Rolff, unpublished data). Generally, the maximal killing values were almost identical in treatments with all the AMPs and their two- and three-way combinations, which means that high concentrations of an AMP(s) may eradicate bacteria with similar efficiencies. Due to fast killing of AMPs and the limit of detection in our system, the real maximal killing rate might be masked at higher concentrations, e.g., concentrations in which the limit of detection is reached within 20 min. However, the MIC and κ significantly varied among single AMPs and two- and three-AMP combinations. As numbers of AMPs increased in combination, the MIC of that combination decreased, with the lowest value seen in three-AMP combinations, and κ was much higher in three-AMP combinations (Fig. 3). More AMPs combined with lower MICs demonstrate that the absolute quantity of AMP needs to be decreased to achieve the same killing. Higher κ values in combined AMPs result in a drastic decrease in bacterial killing rate within a narrow range of concentrations of the AMP(s). The combination of AMPs might improve the efficiency of bacterial killing. Taken together, the decreasing MIC and increasing κ values in combinations with increasing numbers of AMPs suggest that synergism is common in AMP combinations (17, 24, 41).
We observed broad synergistic effects in almost all the two- and three-way combinations. Although some AMPs, like apidaecin, had a relatively weak effect with a high MIC, the killing could still be enhanced by adding one or more AMPs with stronger individual effects. Synergism, albeit not within a pharmacodynamics framework, has been reported for 2-way combinations of antibiotics (8, 42), AMPs (16, 17, 24), antimicrobial peptoids (43), antibiotics and AMPs (44, 45), and AMPs and antimicrobial peptoids (43). The AMPs, originating from different species in our experiment, showed robust synergism, which suggests a general effect.
The molecular mechanisms of interaction, especially antagonisms, of AMPs are largely unknown. As most AMPs target the membrane of pathogens, their interactions are unlikely to directly disrupt the metabolic network in the cell like certain antibiotics. A recent study suggested that synergism was caused by the conjugation of coapplied AMPs, which form a supermolecule and better-stabilized pores (41). This is also confirmed by chemically conjoined synthesized peptides (46). Furthermore, pore-forming peptides can also assist other coapplied transmembrane AMPs to quickly invade bacterial cells and substantially interrupt the metabolism (47).
In our pharmacodynamic model, the important parameter κ depicts the steepness of the pharmacodynamic curve and is a measure of the sensitivity of the response of the bacteria to changes in the concentrations of the antimicrobial substances. A steeper pharmacodynamic curve with higher κ values illustrates that bacteria are very sensitive to the change of concentrations of AMPs and antibiotics, which means that the given antibiotic substance (e.g., combinations of AMPs) has a narrower range of concentrations exerting selection on bacteria.
Additionally, κ value could be an important indicator of resistance selection of given antibiotic agents. Traditionally, the presence of antimicrobial substances above the MIC is thought to favor resistant strains. The mutant selection window (MSW) is defined as the difference in the MICs of a resistant and a susceptible strain (48, 49). Thus, the MSW can be specifically defined as the range between the concentration killing all the sensitive strains and the concentration killing all the resistant strains. Additionally, MSW also can be defined as a range of concentrations which can de novo select mutant strains from a completely sensitive population (50–52). Higher κ values in combinations of AMPs denote a steeper pharmacodynamic curve, which means that the range of concentrations selecting resistance—the MSW—can be narrowed. Especially, the sub-MIC part of the MSW is predicted to be very small for high κ values. A previous theoretical study also demonstrated that the synergistic contribution of the immune system can potentially narrow the mutant selection window of antibiotics (53). We observed a synergistic interaction in combinations of AMPs that mirrors, in the case of LL 19-27, interactions between the immune system and drugs. Higher κ values of AMPs than of antibiotics might partially explain the fact that bacteria are unlikely to develop resistance to AMPs in nature, although resistant strains can emerge under intensive selection in the laboratory (21, 54).
Conclusion.Our study suggests that the synergistic effect between AMPs may be a common phenomenon, as we observed strong synergistic interactions in two-AMP and three-AMP combinations. Interestingly, these three-AMP combinations are even more synergistic than two-AMP combinations. If synergistic interactions of AMPs are ubiquitous, than two practical implications arise: (i) AMPs that strongly synergize with host AMPs should be utilized and (ii) combinations provide the opportunity to reduce side effects, as they lead to an overall reduction in dosage. In the context of innate immunity, selection should favor organisms producing AMP cocktails. This can be considered a cost-efficient way of reducing bacterial loads in a host (19, 55).
Long-lasting coexpression of combinations of AMPs has been recorded in Xenopus laevis (56) and Tenebrio molitor (19), where it is correlated with metabolic suppression. Thus, evolving a more efficient killing system based on a relatively energy-constrained system, which expresses only a limited number of AMPs, is necessary and practical. A function of synergism among AMPs is one of the possible ways to mitigate the costs.
Our results have some implications for the applied use of AMPs as drugs. The production of AMPs is currently expensive (20). The broad synergism observed in our experiment means that combined applications of AMPs could also reduce the consumption of total AMPs just as in the immune system, which could eventually save costs of treatment and reduce toxicity. As humans express AMPs such as LL-37 in their innate immune system, synergisms between these AMPs and AMPs applied as drugs should be taken into account. In our study, the human AMP derivative LL 17-29 synergized with almost all combinations of AMPs. Though resistance to single AMPs evolves readily in vitro, it is might be less likely under combinations (54). It is possible that in some situations combinations delay the development of resistance in medical practice, as pathogens could pay a higher cost to evolve resistance to multidrug treatment (57–60).
ACKNOWLEDGMENTS
We are grateful to Anto Raja Dominic for assistance in statistics and to Alexandro Rodríguez-Rojas for experimental suggestions. We thank Vera Vollenweider and Baydaa El Shazely for comments on an early draft.
G.Y., D.Y.B., R.R.R., and J.R. conceived and designed the experiments. G.Y. performed the experiments. G.Y. and D.Y.B. analyzed the data. D.Y.B., R.R.R., and J.R. contributed reagents/materials/analysis tools. G.Y. and J.R. wrote the paper.
FOOTNOTES
- Received 7 October 2015.
- Returned for modification 22 November 2015.
- Accepted 20 December 2015.
- Accepted manuscript posted online 4 January 2016.
Supplemental material for this article may be found at http://dx.doi.org/10.1128/AAC.02434-15.
- Copyright © 2016 Yu et al.
This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.