The voltage-gated sodium channel (VGSC) is a complex, which is composed of one pore-forming α subunit and at least one β subunit. Up to now, five β subunits are known: β1/β1A, β1B, β2, β3, and β4, encoded by four genes (SCN1B∼SCN4B). It is critical to have a deep understanding of the interaction between β1 and β3 subunits, two subunits which frequently appear in many diseases concurrently. In this study, we had screened out the new template of β1 subunit for homology modelling, which shares higher similarity to β3. Docking studies of the β1 and β3 homology model were conducted, and likely β1 and β3 binding loci were investigated. The results revealed that β1–β3 is more likely to form a di-polymer than β1–β1 based on molecular interaction analysis, including potential energy analysis, Van der Waals (VDW) energy analysis and electrostatic energy analysis, and in addition, consideration of the hydrogen bonds and hydrophobic contacts that are involved. Based on these analyses, the residues His122 and Lys140 of β1 and Glu 66, Asn 131, Asp 118, Glu 120, Glu133, Asn135, Ser 137 of β3 were predicted to play a functional role.