Based on the spatial conformations of human interleukin-6 (hIL-6) derived from nuclear magnetic resonance analysis and human interleukin-6 receptor (hIL-6R) modeled with homology modeling method using human growth hormone receptor as template, the interaction between hIL-6 and its receptor (hIL-6R) is studied with docking program according to the surface electrostatic potential analysis and spatial conformation complement. The stable region structure composed of hIL-6 and hIL-6R is obtained on the basis of molecular mechanism optimization and molecular dynamics simulation. The binding domain between hIL-6 and hIL-6R is predicted theoretically. Furthermore, the especial binding sites that influence the interaction between hIL-6 and hIL-6R are confirmed. The results lay a theoretical foundation for confirming the active regions of hIL-6 and designing novel antagonist with computer-guided techniques.