The bifunctional of the nonadditive kinetic energy in the reference system of noninteracting electrons ([Formula: see text] [ρA, ρB] = Ts[ρA + ρB] - Ts[ρA] - Ts[ρB]) is the key quantity in orbital-free embedding calculations because they hinge on approximations to[Formula: see text] [ρA,ρB]. Since[Formula: see text] [ρA,ρB] is not linear in ρA, the associated potential (functional derivative)[Formula: see text] [ρ,ρB]/δρ|ρ=ρA(r⃗) changes if ρA varies. In this work, for two approximations to[Formula: see text] [ρA,ρB], which are nonlinear in ρA (gradient-free and gradient-dependent), their linearized versions are constructed, and the resulting changes (linearization errors) in various properties of embedded systems (orbital energies, dipole moments, interaction energies, and electron densities) are analyzed. The considered model embedded systems represent typical nonbonding interactions: van der Waals contacts, hydrogen bonds, complexes involving charged species, and intermolecular complexes of the charge-transfer character. For van der Waals and hydrogen bonded complexes, the linearization of[Formula: see text] [ρA,ρB] affects negligibly the calculated properties. Even for complexes, for which large complexation induced changes of the electron density can be expected, such as the water molecule in the field of a cation, the linearization errors are about 2 orders of magnitude smaller than the interaction induced shifts of the corresponding properties. Linearization of[Formula: see text] [ρA,ρB] is shown to be inadequate for the complexes of a strong charge-transfer character. Compared to gradient-free approximation to[Formula: see text] [ρA,ρB], introduction of gradients increases the linearization error.