The Tumor Necrosis Factor (TNF) and the TNF receptor (TNFR) superfamilies are composed of 19 ligands and 30 receptors, respectively. The oligomeric properties of ligands, both membrane bound and soluble, has been studied most. However, less is known about the oligomeric properties of TNFRs. Earlier reports identified the extracellular, membrane-distal, cysteine-rich domain as a pre-ligand assembly domain which stabilizes receptor dimers and/or trimers in the absence of ligand. Nevertheless, recent reports based on structural nuclear magnetic resonance (NMR) highlight the intrinsic role of the transmembrane domains to form dimers (p75NTR), trimers (Fas), or dimers of trimers (DR5). Thus, understanding the structural basis of transmembrane oligomerization may shed light on the mechanism for signal transduction and the impact of disease-associated mutations in this region. To this end, here we used an in silico coarse grained molecular dynamics approach with Martini force field to study TNFR transmembrane homotypic interactions. We have first validated this approach studying the three TNFR described by NMR (p75NTR, Fas, and DR5). We have simulated membrane patches composed of 36 helices of the same receptor equidistantly distributed in order to get unbiassed information on spontaneous proteins assemblies. Good agreement was found in the specific residues involved in homotypic interactions and we were able to observe dimers, trimers, and higher-order oligomers corresponding to those reported in NMR experiments. We have, applied this approach to study the assembly of disease-related mutations being able to assess their impact on oligomerization stability. In conclusion, our results showed the usefulness of coarse grained simulations with Martini force field to study in an unbiased manner higher order transmembrane oligomerization.