For stability analysis of the lumbar spine, the hypothesis presented is that the disc has stress sensors driving feedback mechanism, which could react to the imposed loads by adjusting the contraction of the muscles. Fusion in the motion segment of the lumbar spinal column is believed to alter the stability of the spinal column. To identify this effect finite element (FE) models combined with optimization technique was applied and quantify the role of each muscle and reaction forces in the spinal column with respect to the fusion level. The musculoskeletal FE model was consisted with detailed whole lumbar spine, pelvis, sacrum, coccyx and simplified trunk model. Vertebral body and pelvis were modeled as a rigid body and the rib cage was constructed with rigid truss element for the computational efficiency. Spinal fusion model was applied to L3-L4, L4-L5, L5-S1 (single level) and L3-L5 (two levels) segments. Muscle architecture with 46 local muscles was used as acting directions. Minimization of the nucleus pressure deviation and annulus fiber average axial stress deviation was selected for cost function. As a result, spinal fusion produced reaction changes at each motion segment as well as contribution of each muscle. Longissimus thoracis and psoas major muscle showed dramatic changes for the cases of L5-S1 and L3-L5 level fusion. Muscle force change at each muscle also generated relatively high nucleus pressure not only at the adjacent level but at another level, which can explain disc degeneration pattern observed in clinical study.