The effective fragment potential (EFP) method, a model potential for treating solvent effects and other intermolecular interactions, is interfaced with an electronic structure method, the fragment molecular orbital (FMO) method, that is able to retain high accuracy for ab initio calculations on large molecular systems. The accuracy of the total energies in this novel combined FMO/EFP method is assessed by comparisons with the conventional quantum mechanics (QM)/EFP method. The test cases are water clusters, a peptide, and a dianionic protein (treated with full QM and FMO) combined with water clusters (treated with EFP) at the RHF, B3LYP, and MP2 levels of theory. The basis sets employed range from minimal to augmented double zeta plus polarization. The energy differences between FMO/EFP and the conventional QM/EFP methods are within “chemical accuracy” (1 kcal/mol≈4 kJ/mol).