BIOCHEMISTRYVerified Reference Parity

4-Parameter Logistic (4PL)

The gold-standard symmetric sigmoidal bioassay model for quantitative ELISA and ligand binding calibration.

Primary Disciplines:Bioassay AnalysisPharmacologyELISA ImmunoassaysToxicology

Mathematical Formulation

$$y = D + \frac{A - D}{1 + \left(\frac{x}{C}\right)^B}$$
Plaintext:y = D + (A - D) / (1 + (x / C)^B)

Parameters & Physical Meanings

Detailed parameter breakdown, standard units, valid mathematical constraints, and initial guess heuristic algorithms.

SymbolNameUnitConstraintsPhysical InterpretationInitial Guess Heuristic
$A$Zero-Dose Asymptote (Bottom Plateau)OD / LuminescenceunconstrainedTheoretical response signal at zero analyte concentrationMean of low concentration calibrator standards
$B$Hill Slope FactorDimensionlessB > 0Determines steepness of the sigmoidal transition region1.0
$C$Inflection Point (IC50 / EC50)Analyte concentrationC > 0Analyte concentration at the symmetric midpoint of the calibration curveMidpoint of standard calibrator concentration range
$D$Infinite-Dose Asymptote (Top Plateau)OD / LuminescenceD != AMaximum asymptotic signal achieved at saturating analyte concentrationsMean of highest concentration calibrator standards

When to Choose This Model

  • Quantitative ELISA calibration curve generation and back-calculation of unknown sample concentrations
  • Pharmacological dose-inhibition IC50 screening assays
  • Surface Plasmon Resonance (SPR) steady-state affinity quantification

Typical Scientific Applications

  • Clinical ELISA immunoassay calibration
  • Pharmacological IC50 screening
  • Biological potency bioassays
  • Vaccine titer antibody quantification

Expected Fit Profile & Curve Morphology

Shape: sigmoid
Fitted ModelRaw Data

Example Dataset & Expected Fit Output

Synthetic 12-point ELISA optical density (OD450) calibration curve with IC50 = 8.52 ng/mL

R²:0.9998RMSE:0.0180

Expected Converged Parameters

  • A0.075
  • B1.120
  • C8.520
  • D2.890

Sample Experimental Vectors (12 points)

#x (Independent)y (Observed)
10.0500.082
20.1000.095
30.2500.128
40.5000.194
51.0000.324
62.5000.652
75.0001.054
810.0001.582
+ 4 more points available in AltaiPlot preset

Python / SciPy Reference Implementation

Reference library: scipy.optimize (Trust-Region-Reflective).Exact parameterization parity verified.

import numpy as np
from scipy.optimize import curve_fit

def logistic_4p_model(x, A, B, C, D):
    x_safe = np.maximum(x, 1e-15)
    return D + (A - D) / (1.0 + (x_safe / C)**B)

A_init = np.min(y_data)
D_init = np.max(y_data)
C_init = np.median(x_data)
B_init = 1.0

p0 = [A_init, B_init, C_init, D_init]
bounds = ([-np.inf, 0.1, 1e-15, -np.inf], [np.inf, 20.0, np.inf, np.inf])

popt, pcov = curve_fit(logistic_4p_model, x_data, y_data, p0=p0, bounds=bounds)
Numerical Verification Tolerance: Max Relative Error < 3.2e-6Verified: 2026-08-16

Comparative & Alternative Models

Key decision trade-offs between this model and related functional alternatives:

Alternative Model

Hill Equation (Dose-Response)

The definitive non-linear model for sigmoidal cooperative binding and pharmacological dose-response relationships.

$y = y_{\min} + \frac{y_{\max} - y_{\min}}{1 + \left(\frac{K}{x}\right)^n} = y_{\min} + (y_{\max} - y_{\min})\frac{x^n}{K^n + x^n}$

Related Mathematical Models

Models sharing mathematical formulations or physical domains in AltaiPlot:

biochemistry

Hill Equation (Dose-Response)

The definitive non-linear model for sigmoidal cooperative binding and pharmacological dose-response relationships.

$y = y_{\min} + \frac{y_{\max} - y_{\min}}{1 + \left(\frac{K}{x}\right)^n} = y_{\min} + (y_{\max} - y_{\min})\frac{x^n}{K^n + x^n}$
biochemistry

5-Parameter Logistic (5PL)

Asymmetric sigmoidal bioassay model incorporating an asymmetry factor (S) for non-symmetric dose-response curves.

$y = D + \frac{A - D}{\left(1 + \left(\frac{x}{C}\right)^B\right)^S}$
biochemistry

Michaelis-Menten Kinetics

The cornerstone hyperbolic kinetic model describing single-substrate enzyme reaction rates.

$v = \frac{V_{\max} [S]}{K_m + [S]}$

Frequently Asked Questions (4-Parameter Logistic (4PL))

What is the difference between 4PL and 5PL models?

The 4PL model assumes strict point symmetry around the inflection point C. The 5-Parameter Logistic (5PL) model adds an asymmetry parameter (S) to accurately fit bioassays where the approach to the upper asymptote is noticeably slower or faster than the approach to the lower asymptote.

Scientific References & Citations

  • Rodbard, D. (1974). Statistical quality control and routine data processing for radioimmunoassays and immunoradiometric assays. Clinical Chemistry, 20(10), 1255-1270.[Source / DOI]