Identification of individuals at risk of bone fractures remains challenging despite recent advances in bone strength assessment. In particular, the future degradation of the microstructure and load adaptation has been disregarded. Bone remodeling simulations have so far been restricted to small-volume samples. Here, we present a large-scale framework for predicting microstructural adaptation in whole human vertebrae. The load-adaptive bone remodeling simulations include estimations of appropriate bone loading of three load cases as boundary conditions with microfinite element analysis. Homeostatic adaptation of whole human vertebrae over a simulated period of 10 years is achieved with changes in bone volume fraction (BV/TV) of less than 5 %. Evaluation on subvolumes shows that simplifying boundary conditions reduces the ability of the system to maintain trabecular structures when keeping remodeling parameters unchanged. By rotating the loading direction, adaptation toward new loading conditions could be induced. This framework shows the possibility of using large-scale bone remodeling simulations toward a more accurate prediction of microstructural changes in whole human bones.