The quark-meson coupling model, first developed by Guichon and Thomas, is a nuclear mean field model, in which the quark substructure of nucleons is taken into account. This automatically provides for a mechanism of nuclear saturation.