Computational prediction and design of membrane protein–protein interactions facilitate biomedical engineering and biotechnological applications. Due to their antimicrobial activity, human defensins play an important role in the innate immune system. Human defensins are attractive pharmaceutical targets due to their small size, broad activity spectrum, reduced immunogenicity, and resistance to proteolysis. Protein engineering based modification of defensins can improve their pharmaceutical properties. Here we present an approach to computationally probe defensins’ oligomerization states in the membrane. First, we develop a novel docking and rescoring algorithm. Then, on the basis of the 3D structure of Sapecin, an insect defensin, and a model of its antimicrobial ion-channel, we optimize the parameters of our empirical scoring function. Finally, we apply our docking program and scoring function to the hBD-2 (human β-defensin-2) molecule and obtain structures of four possible oligomers. These results can be used in higher level simulations.