Finite element (FE) modeling is commonly used as a method to investigate the influence of medical devices, such as implants and screws and their effects on the biomechanical behavior of the spine. Another simulation method is a multi-body simulation (MBS), where the model is composed of several non-deformable bodies. MBS solvers generally require a very short computing time for dynamic tasks compared to an FE analysis. Considering this computational advantage, in this study, we examine whether parameters whose values are not known a priory can be determined with sufficient accuracy using MBS model. Therefore, we propose a Many-at-a-time sensitivity analysis method that allows approximating these a priory unknown parameters without requiring long simulation times. This method enables a high degree of MBS model optimization to be achieved in an iterative process. The sensitivity analysis method is applied to a simplified screw-vertebra model, consisting of an anterior anchor implant screw and vertebral body of C4. An experiment described in the literature is used as a basis for developing and assessing the potential of the method for sensitivity analysis and to validate the models action. The optimal model parameters for the MBS model were determined to be c=823224N/m for stiffness and d=488Ns/m for damping. The presented method of parameter identification can be used in studies including more complex MBS spine models or to set initial parameter values that are not available as initial values for FE models.