Hybrid Path Integral Monte Carlo/Molecular Dynamics Approach for the Simulation of Rigid Rotating and Translating Molecules