Inclusion of explicit electron-proton correlation in the nuclear-electronic orbital approach using Gaussian-type geminal functions.