This paper presents a consistent strategy for parametrization of coarse-grained models of chain molecules in dissipative particle dynamics (DPD), where the soft-core DPD interaction parameters are fitted to the activities in solutions of reference compounds that represent different fragments of target molecules. The intercomponent parameters are matched either to the infinite dilution activity coefficients in binary solutions or to the solvent activity in polymer solutions. The respective calibration relationships between activity and intercomponent interaction parameter are constructed from the results of Monte Carlo simulation of the coarse-grained solutions of reference compounds. The chain conformation is controlled by the near neighbor and second neighbor bond potentials, which are parametrized by fitting the intramolecular radial distribution functions of the coarse-grained chains to the respective atomistic molecular dynamics simulations. The consistency, accuracy, and transferability of the proposed parametrization strategy is demonstrated drawing on the example of nonionic surfactants of the poly(ethylene oxide) alkyl ether (CnEm) family. The lengths of tail and head sequences are varied (n = 8-12 and m = 3-9), so that the critical micelle concentration ranges from 10 to 0.1 m.M. The surfactants are modeled at different coarse-graining levels using DPD beads of different diameters. We found consistent agreement with experimental data for the critical micelle concentration and aggregation number, especially for surfactants with relatively long hydrophilic segments. Depending on the system, we observed surfactant aggregation into spheroidal, elongated, or core shell micelles, as well as into irregular agglomerates. Using the models at different coarse-graining levels for the same molecules, we found that the smaller the bead size the better is agreement with experimental data.