Bridge columns are subjected to combined actions of axial force, shear force and bending moment during earthquakes, caused by spatially-complex earthquake motions, features of structural configurations and the interaction between input and response characteristics. Combined actions can have significant effects on the force and deformation capacity of RC columns, resulting in unexpected large deformations and extensive damage that in turn influences the performance of bridges as vital components of transportation systems. This paper evaluates the seismic response of three prototype reinforced concrete bridges using comprehensive numerical models that are capable of simulating the complex soil-structural interaction effects and nonlinear behavior of columns. An analytical approach that can capture the shear-flexural interacting behavior is developed to model the realistic nonlinear behavior of RC columns, including the pinching behavior, strength deterioration and stiffness softening due to combined actions of shear force, axial force and bending moment. Seismic response analyses were conducted on the prototype bridges under suites of ground motions. Response quantities of bridges (e.g., drift, acceleration, section force and section moment etc.) are compared and evaluated to identify the effects of vertical motion, structural characteristics and the shear-flexural interaction on seismic demand of bridges.