Recently, experimental results have demonstrated that the load carrying capacity of the human spine substantially increases under the follower load condition. Thus, it is essential to prove that a follower load can be generated in vivo by activating the appropriate muscles in order to demonstrate the possibility that the stability of the spinal column could be maintained through a follower load mechanism. The aim of this study was to analyze the coordination of the trunk muscles in order to understand the role of the muscles in generating the follower load. A three-dimensional finite element model of the lumbar spine was developed from T12 to S1 and 117 pairs of trunk muscles (58 pairs of superficial muscles and 59 pairs of deep muscles) were considered. The follower load concept was mathematically represented as an optimization problem. The muscle forces required to generate the follower load were predicted by solving the optimization problem. The corresponding displacements and rotations at all nodes were estimated along with the follower forces, shear forces, and joint moments acting on those nodes. In addition, the muscle forces and the corresponding responses were investigated when the activations of the deep muscles or the superficial muscles were restricted to 75% of the maximum activation, respectively. Significantly larger numbers of deep muscles were involved in the generation of the follower load than the number of superficial muscles, regardless of the restriction on muscle activation. The shear force and the resultant joint moment are more influenced by the change in muscle activation in the superficial muscles. A larger number of deep trunk muscles were activated in order to maintain the spinal posture in the lumbar spine. In addition, the deep muscles have a larger capability to reduce the shear force and the resultant joint moment with respect to the perturbation of the external load or muscle fatigue compared to the superficial muscles.