The current article presents the study of neutral stochastic functional differential equations driven by G-Brownian motion in the phase space C_{q}((-infty ,0];mathbb{R}^{n}). The mean-square boundedness of solutions has been derived. The convergence of solutions with different initial data has been investigated. The boundedness and convergence of solution maps have been obtained. In addition, the L^{2}_{G} and exponential estimates of solutions have been determined.