We investigate the robustness of some recent results obtained for homogeneous and isotropic cosmological models with conformally coupled scalar fields. For this purpose, we investigate anisotropic homogeneous solutions of the models described by the action $$ S=\int d^4x \sqrt{-g}\left\{F(\phi)R - \partial_a\phi\partial^a\phi -2V(\phi) \right\}, $$ with general \(F(\phi)\) and \(V(\phi)\). We show that such a class of models leads generically to geometrical singularities if for some value of \(\phi\), \(F(\phi)=0\), rendering previous cosmological results obtained for the conformal coupling case highly unstable. We show that stable models can be obtained for suitable choices of \(F(\phi)\) and \(V(\phi)\). Implications for other recent results are also discussed.