Methodology based on the path integral formulation of quantum mechanics offers a computationally tractable way to include nuclear quantum effects in atomistic simulations of extended condensed-phase systems. Here, we will review some of the required theoretical and algorithmic ingredients, with focus on ab initio molecular dynamics simulations. While all static quantities can be calculated numerically exactly, dynamical information can be obtained from these simulations approximately within the framework of centroid or ring polymer molecular dynamics. Equilibrium isotope effects are accessible through the thermal expectation values of the quantum kinetic energy. Ring polymer contraction is a controlled approximation which allows for an efficient evaluation of physical interactions in a path integral system. It is especially useful with our recent extension to ab initio interactions, where it accelerates practical simulations by up to two orders of magnitude. Importantly, it can be combined with centroid and ring polymer molecular dynamics to obtain approximate dynamics from path integral simulations in a computationally efficient way. We will illustrate the use of this methodology with results of ab initio path integral molecular dynamics simulations of aqueous solutions.