The propagator of two nucleons in infinite nuclear matter is evaluated by a diagonalization of the $pphh$ RPA Hamiltonian. This effective Hamiltonian is non-Hermitian and, for specific density domains and partial waves, yields pairs of complex conjugated eigenvalues representing in-medium bound states of two nucleons. The occurrence of these complex poles in the two-particle Greens function is tightly related to the well known BCS pairing approach. It is demonstrated that these complex eigenvalues and the corresponding bound state wavefunctions contain all information about the BCS gap function. This is illustrated by calculations for $^1S_0$ and $^3PF_2$ pairing gaps in neutron matter which essentially coincide with the corresponding gap functions extracted from conventional solutions of the gap equation. Differences between the bound states in the conventional BCS approach and the $pphh$ RPA are arising in the case of $^3SD_1$ channel in symmetric nuclear matter at low densities. These differences are discussed in the context of transition from BEC for quasi-deuterons to the formation of BCS pairing.