Neurofilaments (NFs) have been proposed to interact with one another through mutual steric exclusion of their unstructured C-terminal "sidearm" domains, producing order in axonal NF distributions and conferring mechanical strength to the axon. Here we apply theory developed for polymer brushes to examine the relationship between the brush properties of the sidearms and NF organization in axons. We first measure NF-NF radial distribution functions and occupancy probability distributions for adult mice. Interpreting the probability distributions using information theory, we show that the NF distributions may be represented by a single pair potential of mean force. Then, to explore the relationship between model parameters and NF architecture, we conduct two-dimensional Monte Carlo simulations of NF cross-sectional distributions. We impose purely repulsive interaction potentials in which the sidearms are represented as neutral and polyelectrolyte chains. By treating the NFs as telechelic polymer brushes, we also incorporate cross-bridging interactions. Both repulsive potentials are capable of reproducing NF cross-sectional densities and their pair correlations. We find that NF structure is sensitive to changes in brush thickness mediated by chain charge, consistent with the experimental observation that sidearm phosphorylation regulates interfilament spacing. The presence of attractive cross-bridging interactions contributes only modestly to structure for moderate degrees of cross-bridging and leads to NF aggregation for extensive cross-bridging.