The complex inheritance involved in multiple sclerosis (MS) risk has been extensively investigated, but our understanding of MS genetics remains rudimentary. In this study, we explore 51 single nucleotide polymorphisms (SNPs) in 36 candidate genes from the inflammatory pathway and test for gene-gene interactions using complementary case-control, discordant sibling pair, and trio family study designs. We used a sample of 421 carefully diagnosed MS cases and 96 unrelated, healthy controls; discordant sibling pairs from 146 multiplex families; and 275 trio families. We used multifactor dimensionality reduction to explore gene-gene interactions. Based on our analyses, we have identified several statistically significant models including both main effect models and two-locus, three-locus, and four-locus epistasis models that predict MS disease risk with between approximately 61% and 85% accuracy. These results suggest that significant epistasis, or gene-gene interactions, may exist even in the absence of statistically significant individual main effects.