In this protocol, we describe how to make Loewe additivity-based drug interaction measurements for pairwise and three-way drug combinations.
A synergistic drug combination has a higher efficacy compared to the effects of individual drugs. Checkerboard assays, where drugs are combined in many doses, allow sensitive measurement of drug interactions. However, these assays are costly and do not scale well for measuring interaction among many drugs. Several recent studies have reported drug interaction measurements using a diagonal sampling of the traditional checkerboard assay. This alternative methodology greatly decreases the cost of drug interaction experiments and allows interaction measurement for combinations with many drugs. Here, we describe a protocol to measure the three pairwise interactions and one three-way interaction among three antibiotics in duplicate, in five days, using only three 96-well microplates and standard laboratory equipment. We present representative results showing that the three-antibiotic combination of Levofloxacin + Nalidixic Acid + Penicillin G is synergistic. Our protocol scales up to measure interactions among many drugs and in other biological contexts, allowing for efficient screens for multi-drug synergies against pathogens and tumors.
Drug combinations may exhibit surprisingly high or low effect on a phenotype given the effects of constituent drugs, corresponding to synergistic or antagonistic drug interactions, respectively1,2,3. The use of synergistic combinations may allow dose-escalation for efficacy increase and dose-reduction for side-effect relief. Combination treatments may also apply multiple setbacks to cellular machinery, thereby blocking potential evolutionary escape mechanisms to resistance4. Therefore, combinations of three or more drugs are routinely used in pathogen or cancer treatment5.
Synergy and antagonism are defined by a comparison between the observed effect of a combination versus an expected effect given individual drug effects. Among the models for drug interactions, Loewe additivity is the most stringent and has a well-defined null model (Figure 1)6, and the inferred synergy/antagonism interaction is independent of the drug concentration used6,7. However, the Loewe model is experimentally costly even for a pairwise interaction test. Drug interaction assays traditionally comprise of a 2D matrix of drug concentration combinations (a checkerboard assay) (Figure 2). If 5 doses are used for each drug, then 25 combinations are required, corresponding to one half of a microplate if experiments are conducted in replicate. The cost of this approach prohibits synergy measurement by the Loewe additivity model for multi-drug combinations (Figure 3). For example, to test a 10-way interaction, traditional methods would require more than 100 thousand microplates, barring experimental measurement of high-order synergies by the stringent, well-theorized and concentration-independent Loewe additivity model8.
Current clinical treatments utilize only a fraction of possible drug combinations. For example, the standard treatment of active tuberculosis is a combination of three antibiotics. There are approximately 20 antibiotics used in Mycobacterium tuberculosis (Mtb) treatment. There are 1140 possible 3-way combinations among 20 drugs, each with the potential to have a strong synergy against Mtb. As there has been no cost-effective method to measure drug interactions among many drugs, potentially life-saving synergistic combinations remain untested.
Here, we describe a simple protocol to measure pairwise and three-way drug interactions by sampling only the diagonal of a checkerboard assay (Figure 4 and Figure 5). The underlying concept of sampling the diagonal of a checkerboard experiment was theorized by Berenbaum in his seminal work in 19789. Yet, this approach has only recently been applied to drug synergy screens10,11,12. We present our protocol with Escherichia coli (E. coli) and the growth phenotype. However, we note that the protocol is independent of the biological species and phenotype of interest, and hence may be applied to the measurement of high-order drug synergy in other biological contexts.
NOTE: Any small molecule that inhibits E. coli bacteria growth can be used for the diagonal method. In this protocol, levofloxacin (LEV), nalidixic acid (NAL), and Penicillin G (PNG) will be used as an example, since these drugs show a salient three-way synergy. The workflow of this protocol is shown as Figure 6. Carry out all steps at room temperature. Use fresh aliquots of bacteria and drugs each day. Perform the experiment under biosafety levels appropriate for E. coli.
1. Preparation Steps
2. Serial Dilution Dose-Response Experiment
3. Linear Dilution Dose-Response Experiment
4. Diagonal Drug Interaction Experiment
5. Diagonal Drug Interaction Scores
Previously, we have reported the pairwise interactions among three drugs: LEV, NAL and PNG based on testing in miniaturized checkerboard assays, where two drugs were combined in a 4 x 4 matrix13,14. While NAL and LEV were synergistic, PNG was reported to be antagonistic with both LEV and NAL13,14. Here, we verified these pairwise interactions and measured the three-way interaction among these three drugs using a diagonal assay. Our results demonstrate that LEV+NAL+PNG is a synergistic 3-way antibiotic combination. Schematic representations for the results of the individual experimental procedure sections were given on the right side of Figure 7 and Figure 8. Here, we present and interpret representative raw results from three plate readings, which are given in Figure 9. The top plate reading corresponds to serial and linear dilution experiments conducted in steps 2 and 3. The bottom two plate readings are duplicate interaction plates conducted in step 4.
The raw data in Figure 9 shows that top growth is around 0.55, but there is a 0.05 optical density of the media itself, as observed in the OD600 of the high drug concentrations where there is no growth. Therefore, we define IC50 as (0.55-0.05)/2 = 0.25. For each dose-response, the wells located closest to this value are shown with orange.
The upper half of Figure 9A shows the results from step 2, serial dose-response experiment. The IC50 wells for LEV is at column 10 in two replicates, which correspond to 4 ng/mL. The IC50 for NAL and PNG are at 3 µg/mL and 25 µg/mL, respectively. These concentrations correspond to the 1x concentration shown in Figure 7B. The lower half of Figure 9B shows the results from step 3, linear dose-response experiments. LEV, NAL and PNG's IC50 are found at 0.4x, 0.8x and 1.2x, respectively. These concentrations are assigned as the 1X IC50 for step 4.
Two plates corresponding to two replicate experiments are shown in Figure 9B, where the IC50 wells are shown with orange. In plate 1, all single drugs have their IC50 at 1x concentration. The expected IC50 for the pairwise or three-way combination is calculated by the arithmetic mean of constituent drugs, making expected IC50 for all combinations also 1x concentration. In plate 2, LEV and PNG have their IC50 at the 1x concentration, but NAL IC50 is at 1.2x. The expected IC50 for each combination is defined using the arithmetic means of these IC50 values. For example, the expected IC50 for LEV+NAL and NAL+PNG is 1.1x. The drug interaction score (FIC) for each combination is calculated by dividing observed IC50 with the expected IC50, as shown on the right side of the plates. Inspection of the FIC scores of the two plates demonstrates that LEV+NAL and LEV+NAL+PNG are synergistic, while LEV+PNG and NAL+PNG are antagonistic. FIC scores obtained in the two plates are in agreement, supporting the reliability of the protocol.
Figure 1: Null model definition for Loewe additivity drug interaction model. Two 5 x 5 matrices on a microplate, where drug A is linearly increased in one axis as shown at left, and a top concentration of drug inhibits a quantifiable phenotype (top right). The addition of these matrices is shown in the middle, where the lines connecting equipotent drug A concentrations in each single drug have the same concentration as drug A. When cells are added on this "checkerboard" of drug concentration combinations, it is expected that the recorded phenotype on this line will be equivalent for the wells it connects. In this self-self drug interaction experiment, the isobole depicting a phenotype (shown with a dashed green line) is expected to be linear, defining the additivity null model. Please click here to view a larger version of this figure.
Figure 2: Pairwise drug interactions according to Loewe additivity drug interaction model. When two drugs are combined in a checkerboard assay as in Figure 1, the observed isophenotypic contour may be straight, convex or concave. At right, possible isophenotypic contours (dashed green lines) are superimposed on checkerboard assays. Drug combinations with straight isophenotypic contours are Loewe-additive, as the contours are not different from a self-self drug interaction, which is Loewe-additive by definition. When the isophenotypic contour is significantly concave or convex, the combination is synergistic or antagonistic, respectively. Please click here to view a larger version of this figure.
Figure 3: Three-way drug interactions according to Loewe additivity drug interaction model. Similar to the checkerboard assay for pairwise interactions, three-drugs are combined in a 3D grid (a "checkercube"), where each drug is linearly increased in one axis. If the three drugs were identical, the isophenotypic surface is expected to be flat, defining the additivity for three drug combinations. If the surface is more concave or convex than this Loewe-additive null model, drug combinations are synergistic or antagonistic, respectively. Please click here to view a larger version of this figure.
Figure 4: Diagonal method to measure pairwise drug interactions. For each checkerboard assay, only the regions shown in magenta rectangles are measured. i and ii are linearly increasing single drug concentrations. iii is assessed by making a 1:1 mixture of two drugs and linearly titrating this mixture as if it was a single drug. The FIC is equal to the observed IC50 in the combination divided by the expected IC50 of two single drugs. For the Loewe additivity model, the expected IC50 is approximated by the average IC50 of the two single drugs. A FIC value is 1 for Loewe-additive pairs and is lower or higher than 1 for synergistic or antagonistic pairs, respectively2,12. Please click here to view a larger version of this figure.
Figure 5: Diagonal method to measure three-way drug interactions. For each checkercube assay, only the regions shown are measured. i, ii and iii are linearly increasing single drug concentrations. iv is measured by making a 1:1:1 mixture of three drugs and linearly titrating this mixture as if it were a single drug. The Fractional Inhibitory Concentration is equal to the observed IC50 in the combination divided by the expected IC50 given three single drugs. Please click here to view a larger version of this figure.
Figure 6: Workflow for the diagonal method protocol described herein and the setup details for each microplate. The plate shown in Day 4 is conducted in duplicate. Please click here to view a larger version of this figure.
Figure 7: Serial and linear dilution dose-response experiments. (A) Preparation of serial dilution dose-response for one drug and corresponding final drug concentrations. E. coli cells are added to the plate; growth is recorded after 16 h. Serial IC50, shown in orange, is selected for each drug for use in the following day's linear dilution dose-response experiments. (B) Preparation of linear dilution dose-response for one drug and corresponding final drug concentrations. E. coli cells are added to the plate; growth is recorded after 16 h. For each drug, the selected IC50s that will be used in the following day's drug interaction experiments are shown in orange. Please click here to view a larger version of this figure.
Figure 8: Drug interaction experiments. Preparation of interaction experiments is similar to single drug linear dose-responses, except that a 1:1 or 1:1:1 mixture of drugs is used for two or three drug dose-responses, respectively. Observed IC50s for each dose-response, depicted in orange, are used to calculate FIC scores. Please click here to view a larger version of this figure.
Figure 9: Representative experiment results. Representative results obtained using the described protocol are shown, with details provided in the main text. Please click here to view a larger version of this figure.
The use of drug combinations against pathogens or tumors is an attractive prospect, especially under the circumstances of the drying antibiotic pipeline. However, this potential is hampered by at least two difficulties. The first difficulty is the astronomical number of possible combinations. There are, for example, 4950 possible pairwise combinations among 100 antibiotics. All possible combinations among 100 antibiotics (2100) is on the same order of magnitude with the number of bacteria on Earth (~1030). How to predict strongly synergistic combinations among these possibilities has been the subject of numerous computational studies. The second difficulty is the measurement of high-order drug interactions. Consider that a computational platform may suggest that a certain 10-drug combination is strongly synergistic against a certain pathogen. Traditional methods to test drug interactions is too costly to verify or refute this hypothesis, therefore the study of synergy among many drugs has been outside the boundaries of scientific inquiry. The diagonal method, which was first proposed almost 30 years ago and was used in a few recent synergy screens provide a strong foundation for the first problem, by allowing testing of the interaction among many pairs. It solves the second problem by an informative sampling of the traditional assays and allows the study of high-order drug interactions.
We importantly note that our protocol uses a linear dosing for drug interaction measurements, to provide the sensitivity for detecting even weak interactions. Establishing the right concentration range for linear dosing is a challenging task. By first performing a serial dilution, we make an informed decision about the search space for linear dosing. However, the protocol can be modified to use 2-fold or higher serial dilutions for drug interaction testing. Such a modification would shorten the experiment time and allow the testing of more interactions; however, it would have sensitivity to detect only strongly synergistic or antagonistic interactions.
The protocol we described shows the measurement of pairwise or three-way interactions. A critical aspect of the protocol is that single agents are on the same plate as the combination, to minimize bias due to plate variations. Therefore, the protocol can be adjusted to measure interactions up to 7-way combinations by trivial modifications. Combinations of more than 7 drugs will require more than one 96-well microplate and additional considerations must be taken to ensure correct data integration, such as inter-plate replicates.
A notable limitation of the diagonal method is the restriction that each drug in the assay must inhibit the phenotype of interest. Therefore, the diagonal method is not useful for understanding the interactions among active agents and inert adjuvants. Such 'potentiating' interactions may be studied under alternative models such as Bliss or Highest Single Agent models.
An important consideration for the analysis of high-order drug interactions is the null model choice for the "expected IC50." When two drugs are combined, the combination's effect can only be compared to the single drug effects. When three drugs are combined, the combination's effect can be compared to the single effects or pairwise effects. For example, if all pairwise combinations of three drugs are synergistic, then it may be expected that these drugs will show a three-way synergy. A three-way interaction's deviation from what is expected from pairwise interactions has been recently dubbed "emergent interaction"16,17. For simplicity, our protocol describes the measurement of the three-way combination's "net interaction," which defines the null model as single drug effects. However, the data that is obtained from the protocol may also be used to compute the emergent interaction of the three-way combination. In our analysis, we defined the expected IC50 of the three-way combination as the average of the single drug IC50s. Alternatively, the expected IC50 can be defined as the average of the IC50s of pairwise combinations (~1.1-1.2). When the observed IC50 is divided by this alternative expected IC50, the obtained FIC provides the emergent FIC for the three-way combination, as previously described12. This consideration reveals that LEV+NAL+PNG is more synergistic than what would be expected from the pairwise interactions among three drugs, demonstrating that LEV+NAL+PNG has emergent synergy.
The authors have nothing to disclose.
This work was funded by NIGMS Grant P50GM107618. The authors thank Zohar B. Weinstein for the insightful comments and suggestions on the manuscript.
1.5 mL Semi Micro Cuvette | VWR | 97000-586 | |
1.5 mL Eppendorf Microcentrifuge Tubes | USA Scientific | 4036-3204 | |
1000 µL Tips | Geneseesci | 24830 | |
14 mL Breathable Cell Culture Tube | VWR | 60819-761 | |
20 µL Tips | Geneseesci | 24804 | |
200 µL Tips | Geneseesci | 24815 | |
37 °C Incubator | Panasonic | MIR-262-PA | |
37 °C Shaker Incubator | Thermo Scientific | SHKE8000 | |
5 mL Cell Culture Serological Pipette | VWR | 53300-421 | |
96-well Microplates | VWR | 15705-066 | |
Breathable Sealing Film | USA Scientific | 2920-0010 | |
DMSO | Sigma | 41647 | |
Escherichia coli | ATCC | 700926 | |
Glycerol | Sigma | G9012 | |
LB Broth Powder | RPI | L24065 | |
Levofloxacin | Sigma | 28266 | |
Micropipette | GILSON | PIPETMAN Classic | |
Microplate reader | BioTek | Synergy H1 | |
Multichannel micropipette | VistaLab | 1060 | |
Nalidixic acid | Sigma | N8878 | |
Penicillin G | Sigma | P3032 | |
Pipette Pump | Drummond | 4-000-501 | |
Reagent Reservoir | VWR | 89094-658 | |
Spectrophotometer | BIO-RAD | 1702525 | |
Vortex Mixer | Fisher Scientific | 10-320-807 |