We combine conditioning techniques with sparse grid quadrature rules to develop a computationally efficient method to approximate marginal, but not necessarily univariate, posterior quantities, yielding approximate Bayesian inference via Sparse grid Quadrature Evaluation (BISQuE) for hierarchical models. BISQuE reformulates posterior quantities as weighted integrals of conditional quantities, such as densities and expectations. Sparse grid quadrature rules allow computationally efficient approximation of high dimensional integrals, which appear in hierarchical models with many hyperparameters. BISQuE reduces computational effort relative to standard, Markov chain Monte Carlo methods by at least two orders of magnitude on several applied and illustrative models. We also briefly discuss using BISQuE to apply Integrated Nested Laplace Approximations (INLA) to models with more hyperparameters than is currently practical.