We present a linear-scaling symmetry-adapted perturbation theory (SAPT) method that is based on an atomic orbital (AO) formulation of zeroth-order SAPT (SAPT0). The non-dispersive terms are realized with linear-scaling cost using both the continuous fast multipole method (CFMM) and the linear exchange (LinK) approach for integral contractions as well as our efficient Laplace-based coupled-perturbed self-consistent field method (DL-CPSCF) for evaluating response densities. The reformulation of the dispersion term is based on our linear-scaling AO Møller-Plesset second-order perturbation theory (AO-MP2) method, that uses our recently introduced QQR-type screening [S. A. Maurer, D. S. Lambrecht, J. Kussmann, and C. Ochsenfeld, J. Chem. Phys.138, 014101 (2013)] for preselecting numerically significant energy contributions. Similar to scaled opposite-spin MP2, we neglect the exchange-dispersion term in SAPT and introduce a scaling factor for the dispersion term, which compensates for the error and at the same time accounts for basis set incompleteness effects and intramonomer correlation. We show in extensive benchmark calculations that the new scaled-dispersion (sd-)SAPT0 approach provides reliable results for small and large interacting systems where the results with a small 6-31G** basis are roughly comparable to supermolecular MP2 calculations in a triple-zeta basis. The performance of our method is demonstrated with timings on cellulose fragments, DNA systems, and cutouts of a protein-ligand complex with up to 1100 atoms on a single computer core.