The relationship between the parameters of shear transformation zone (STZ) model, mechanical behavior and serrated plastic flow of Zr65Cu15Ni10Al10 metallic glasses (MGs) were investigated. In this study, mesoscale modeling using the kinetic Monte Carlo methods (kMC), and the bi-dimensional 2-D finite element method (FEM) is implemented in ABAQUS software using the subroutine UMAT to investigate the effects of variation parameters model on mechanical behavior of Zr65Cu15Ni10Al10 MGs ribbon. Using an atomic force microscope (AFM) and a scanning electronic microscope (SEM), we characterize the propagation and mutual interaction of shear bands (SBs) at the surface of notched ribbons. At the interaction sites, the height of SB can be decrease by the offset of the two bands. Emphasis is placed on the effect of strain rate, shear stress and STZ activation volume on the macroscopic mechanical behavior predicted by FE simulation. The relationship between the predicted yield strength and the physical parameter values is analyzed via a stress-field estimation. Our simulation results demonstrate that the variation of parameters that can dominate STZ occurrences has an influence on the stress-strain response of MGs.