Home / Science / Health and medicine / Pharmacology / Drug combination synergy (Bliss, Loewe, CI)Tool #472

Drug synergy calculator: Bliss independence, Loewe additivity, HSA and combination index heatmap from a dose-response matrix

Paste a checkerboard dose matrix of % inhibition for two drugs. Get the excess over Bliss independence, the highest single agent and Loewe additivity as heatmaps, single-drug Hill fits, the Chou-Talalay combination index for every well and summary synergy scores.

Keep it one tap away for everyday use.

Bookmark Drug combination synergy (Bliss, Loewe, CI)?

Press Ctrl + D now. Your browser opens its own bookmark box, where you choose the name and folder.

  1. Matrix
  2. Hill fits
  3. Heatmap
  4. Scores
Nothing is uploadedWorks offline once open

For study and teaching, not medical advice. Check with a qualified professional.

Dose matrix (% inhibition)

Responses are % inhibition (0 = no effect, 100 = full effect); convert viability first (100 − % viability). Fixing the top at 100 % is the median-effect model of Chou and Talalay; fit it when a drug clearly levels off lower. Bliss and HSA use your measured single-drug wells; Loewe and the combination index use the fitted single-drug curves.

Download

What it does

When two drugs are given together, the combined effect can be more than expected (synergy), about as expected (additivity) or less (antagonism). “Expected” needs a reference model, and the field uses several. This page takes a checkerboard experiment (a dose matrix of two drugs, each well’s response as % inhibition) and compares every combination well with three reference models, draws the differences as a heatmap and gives summary scores.

It is for study and teaching, not medical advice. It analyses data from an experiment; it does not suggest doses.

How to use it

  1. Paste the matrix. The first row holds drug B’s doses after a corner cell; each following row starts with a drug A dose, then the responses at each drug B dose. Both drugs need a dose of 0, so the matrix includes each drug alone. Up to 20 × 20 wells. Or use Open a file.
  2. Responses must be % inhibition: 0 means no effect and 100 a full effect. If you measured viability, use 100 − % viability.
  3. Name the drugs and their units for the labels.
  4. Choose the top of the single-drug curves. Fixed at 100 % is the median-effect model of Chou and Talalay. Fitted lets each drug level off lower; if a curve does not level off within your doses, the top is set back to 100 % and the page says so.
  5. Pick the heatmap: excess over Bliss, HSA or Loewe, the combination index (CI), or the observed matrix.

The models

  • Bliss independence: the drugs act independently, so the expected inhibition is EA + EB − EA·EB, using the measured single-drug wells.
  • Highest single agent (HSA): the combination should beat the better of the two drugs alone at those doses.
  • Loewe additivity: each drug behaves like a dilution of the other. The expected effect E solves a / DxA(E) + b / DxB(E) = 1, where DxA and DxB come from Hill curves fitted to each drug alone: E = top × dʰ / (EC50ʰ + dʰ).
  • Combination index: CI = a / DxA + b / DxB at the observed effect. Below 1 is synergy, 1 additivity, above 1 antagonism. It is not defined when the observed effect is beyond a curve’s top.

The excess of a well is observed − expected in % points. Scores are averaged over the combination wells (both doses above 0); the CI summary is a geometric mean.

Worked example

The example is a made-up 6 × 6 matrix. With the top fixed at 100 %, drug A’s curve has an EC50 of 0.8164 µM (Hill slope 1.20) and drug B’s 13.34 nM (slope 0.884). The page shows a mean excess over Bliss of +2.31 (highest +6.23 at 0.3 µM and 1 nM), a mean HSA excess of +9.20, a mean Loewe excess of +5.28 and a combination index of 0.643, with 22 of the 24 wells below 0.9. In the well at 0.3 µM and 3 nM the observed inhibition is 43.8 %, Bliss expects 39.5 % and Loewe 36.6 %, and the CI is 0.75. So the drugs look mildly synergistic by Loewe and CI, and close to additive by Bliss.

Tips

  • The models answer different questions and often disagree, as here. Report which one you used.
  • Use replicate plates and average them before pasting; a single noisy well can make a strong-looking spot.
  • Doses should bracket each drug’s EC50 so the single-drug curves are well defined.

FAQ

Is ZIP included? Not yet; Bliss, HSA, Loewe and CI are.

Why are the single-drug wells 0 on the Bliss map? Bliss uses them as its reference, so they match by definition.

Is my data uploaded? No. Everything runs in this page.

Where the numbers come from

  • Bliss independence: C. I. Bliss, Ann. Appl. Biol. 26 (1939) 585. Loewe additivity: S. Loewe and H. Muischnek, Arch. Exp. Pathol. Pharmakol. 114 (1926) 313. Combination index: T.-C. Chou and P. Talalay, Adv. Enzyme Regul. 22 (1984) 27; M. C. Berenbaum, Pharmacol. Rev. 41 (1989) 93. Single-drug wells for Bliss and HSA as in SynergyFinder: A. Ianevski et al., Bioinformatics 33 (2017) 2413. Formulas only; no data or tables copied.
  • Calculations written by fibub. Checked in the tests against an independent Python implementation (SciPy curve_fit Hill fits with several starts, brentq for the Loewe surface) on the example and 20 random matrices: Hill parameters agree to 3 × 10⁻⁵, every well to 3 × 10⁻⁴ % points; also against sham combinations (a drug with itself gives CI = 1).
  • The example matrix is made up.

Standard scientific values, such as constants, atomic weights and melting points, are facts found in every textbook; we name where ours were checked. If you are an author or publisher and want a value, a credit or a table changed or removed, tell us and we will act quickly.

Last updated 2026-11-09.

Related tools

Enzyme kinetics fitter #112 Paste substrate concentrations and initial rates, with or without an inhibitor, and fit Michaelis–Menten, substrate inhibition, Hill, or competitive, uncompetitive, non-competitive and mixed inhibition by non-linear least squares. Get Km, Vmax, Ki, kcat and kcat/Km with errors, a model comparison by AICc, and the Lineweaver–Burk, Eadie–Hofstee and Hanes–Woolf values side by side. LC50 probit analysis #151 Find the LC50 or EC50 from a toxicity test where a number of organisms die or respond at each concentration. Fit a probit or logit line by maximum likelihood, correct for deaths in the control, and get LC10, LC50, LC90 or any other level with 95 % fiducial limits, the slope, a goodness-of-fit test and plots. Binding curve fitter #114 Paste ligand concentrations with total and non-specific binding, or a titration signal, and fit one site, one site with non-specific binding, two sites, Hill or tight binding with ligand depletion. Get Kd, Bmax and the Hill coefficient with errors, a model comparison by AICc, the Scatchard plot and an IC50 to Ki converter. Dose-response and ELISA curve fitter (4PL) #29 Fit a four-parameter logistic curve to dose-response or ELISA data. Get the EC50 or IC50 with its 95 % interval and the Hill slope, and read unknown samples off the curve. Pharmacokinetics simulator (teaching) #473 See what a regimen does to an invented drug: one- or two-compartment model, IV bolus, IV infusion or oral absorption, repeated at an interval. Draws the concentration-time curve and gives Cmax and Cmin at steady state, the average level, the accumulation ratio and the half-lives. A teaching tool: it never works out an amount to give. Non-compartmental PK analysis (NCA) #474 Paste a measured concentration-time profile and get Cmax, Tmax, AUC0-t by linear or linear-up/log-down trapezoids, the terminal slope lambda_z by best fit or chosen points, half-life, AUC0-inf, % extrapolated, AUMC, MRT, and CL/F and Vz/F from the amount given in your study. Benchmark dose (BMD) calculator #475 Fit five dose-response models to quantal (yes/no) toxicity data by maximum likelihood, get the benchmark dose at a chosen extra risk and its lower limit (BMDL) by profile likelihood, compare models by AIC and goodness of fit, and plot them. Publication plot maker #27 Paste columns or open a CSV or Excel file and get a clean scientific figure with error bars, sized for a journal column. Download vector SVG or PDF, or PNG at 300 or 600 dpi. Your data stays in your browser.