The vertebral column is the most frequent site of metastatic involvement of the skeleton. Due to the proximity to the spinal cord, from 5% to 10% of all cancer patients develop neurologic manifestations. As a consequence, fracture risk prediction has significant clinical importance. In this study, we model the metastatically involved vertebra so as to parametrically investigate the effects of tumor size, material properties and compressive loading rate on vertebral strength. A two-dimensional axisymmetric finite element model of a spinal motion segment consisting of the first lumbar vertebral body (no posterior elements) and adjacent intervertebral disc was developed to allow the inclusion of a centrally located tumor in the vertebral body. After evaluating elastic, mixed, and poroelastic formulations, we concluded that the poroelastic representation was most suitable for modeling the metastatically involved vertebra's response to compressive load. Maximum principal strains were used to localize regions of potential vertebral trabecular bone failure. Radial and axial vertebral body displacements were used as relative indicators of spinal canal encroachment and endplate failure. Increased tumor size and loading rate, and reduced trabecular bone density all elevated axial and radial displacements and maximum tensile strains. The results of this parametric study suggest that vertebral tumor size and bone density contribute significantly to a patients risk for vertebral fracture and should be incorporated in clinical assessment paradigms.