The integral equation formalism polarizable continuum model (IEF-PCM) for solvent effects with the default solvent (water) and solvent parameters, together with the density functional theory method at 6-31G(d) level, was used to optimize molecular structures for polychlorinated biphenyl (PCB) congeners. Four molecular descriptors were selected to develop quantitative structure–activity relationship (QSAR) models for the depuration rate constants (kd) of 63 PCB congeners in a juvenile rainbow trout (Oncorhynchus mykiss). The optimal multiple linear regression (MLR) model has the correlation coefficient R of 0.933 and the root mean square (rms) error of 0.0681 for the total set of 63 PCB congeners. The support vector regression model has R of 0.953 and rms error of 0.0576 for the total set. Both the MLR and SVM QSAR models in this paper were accurate and acceptable compared with other QSAR models for the depuration rate of PCB congeners reported in references. Thus, applying IEF-PCM and B3LYP/6-31G(d) calculations for molecular descriptor derivation of PCB congeners is successful.