We apply the variational method in Hylleraas coordinates to solve the energy eigenvalue problem for antiprotonic helium molecular systems including (p) over bar He-3(+) and (p) over bar He-4(+). The numerical accuracy on the non-relativistic energies is shown to reach 10 (17), thus the precision of our results is only limited by the width of the metastable states. Expectation values of the Dirac delta operators for these states are also calculated. (C) 2016 Elsevier B.V. All rights reserved.