This paper is devoted to dealing with the issue of the estimation of the domain of attraction (DOA) for highly nonlinear neutral complex networks (HNNCNs) with time delays. Firstly, by the Razumikhin approach, we establish several novel lemmas on the estimation of DOA for highly nonlinear neutral differential systems. The cases of bounded non-differentiable delays and unbounded proportional delays are discussed, respectively. Subsequently, by utilizing the proposed lemmas, combining the Lyapunov stability theory and inequality technique, the estimation of DOA on HNNCNs with bounded delays or proportional delays is derived when the chosen control gain is sufficiently large. If initial values start from DOA, then the states of systems will exponentially or polynomially converge to the equilibrium point, which means that the local exponential or polynomial synchronization of HNNCNs is realized. Additionally, the weighted outer-coupling matrix of complex networks is not required to be symmetric, which means that the derived results can be applied to both the undirected networks and directed networks. Finally, several numerical examples are provided to illustrate the feasibility of theoretical findings.