In general, the interaction conditions between the discrete stones are not taken into account by structural engineers during the modeling and analyzing of historical stone bridges. However, many structural damages in the stone bridges occur due to ignoring the interaction conditions between discrete stones. In this study, it is aimed to examine the seismic behavior of a historical stone bridge by considering the interaction stiffness parameters between stone elements. For this purpose, Tokatli historical stone arch bridge was built in 1179 in Karab??k-Turkey, is chosen for three-dimensional (3D) seismic analyses. Firstly, the 3D finite-difference model of the Tokatli stone bridge is created using the FLAC3D software. During the modeling processes, the Burger-Creep material model which was not used to examine the seismic behavior of historical stone bridges in the past is utilized. Furthermore, the free-field and quiet non-reflecting boundary conditions are defined to the lateral and bottom boundaries of the bridge. Thanks to these boundary conditions, earthquake waves do not reflect in the 3D model. After each stone element is modeled separately, stiffness elements are defined between the stone elements. Three situations of the stiffness elements are considered in the seismic analyses; a) for only normal direction b) for only shear direction c) for both normal and shear directions. The earthquake analyses of the bridge are performed for these three different situations of the bridge. The far-fault and near-fault conditions of 1989 Loma Prieta earthquake are taken into account during the earthquake analyses. According to the seismic analysis results, the directions of the stiffness parameters seriously changed the earthquake behavior of the Tokatli bridge. Moreover, the most critical stiffness parameter is determined for seismic analyses of historical stone arch bridges.