Using the numerical model of Scheutjens and Fleer we investigated, on a self-consistent field level, the equilibrium structure of the neurofilament brush formed by the projection domains of NF-H, NF-M and NF-L proteins. Although the actual AA sequences in the projection domains are coarse grained, the different (realistic) solubilities of AA residues and the specific distribution of its intrinsic charges inside the arms of the NF proteins, are taken explicitly into account. We collect strong evidence that the electrostatic interactions are a dominant force that controls the NF brush structure. There exists a remarkable spatial separation of the H, M and L tails. In a dephosphorylated NF we found confined and flower-like conformations for the H and M projection domains, respectively. We demonstrate that the ionization of KSP repeats in NF proteins triggers a conformational transition in the H tail that leads to the expulsion of its terminal (KEP) domain to the periphery of the NF brush. We argue that the phosphorylation of the NF proteins in axons can both increase the interfilament distance and stabilizes cross-bridges between neurofilaments.