ABSTRACT:We report an analysis of three schemes for fragment reassociation using multicentered multipolar expansions derived from ab initio quantum wave functions at the Hartree-Fock/6-31G * LCAO level, two of them involving single-bond partitioning in the peptide bond region, and the third one using a partially overlapping procedure based on a methodology proposed by Vigné-Maeder 21 (OME-overlap of multipolar expansions-reassociation method). The effects of different peptide junction treatments in the derivation of molecular electrostatic potentials and molecular electric fields of three peptide sequences are discussed. The results show that the OME reassociation method gives a better and a more homogeneous description of both the potential and the electric field than the other two treatments. We conclude that the OME method is the most indicated for studies involving electrostatic properties of proteins. Our results also indicate that the use of multicentered multipolar expansions coupled to the OME treatment is the best choice in protein studies including solvent effects using, for example, a continuum boundary method to solve the linearized Poisson-Boltzmann equation.