The subgrid-scale closure method developed by Germano et al. [Phys. Fluids A 3, 1760 (1991)] is modified by use of a least squares technique to minimize the difference between the closure assumption and the resolved stresses. This modification removes a source of singularity and is believed to improve the method's applicability.