Acidic proteins found in mineralized tissues act as nature's crystal engineers, where they play a key role in promoting or inhibiting the growth of minerals such as hydroxyapatite (HAP), Ca-10(PO4)(6)(OH)(2), the main mineral component of bone and teeth. Key to understanding the structural basis of protein-crystal recognition and protein control of hard tissue growth is the nature of interactions between the protein side chains and the crystal surface. In an earlier work we have measured the proximity of the lysine (K6) side chain in an SN-15 peptide fragment of the salivary protein statherin adsorbed to the Phosphorus-rich surface of HAP using solid-state NMR recoupling experiments. N-15 {P-31} rotational echo double resonance (REDOR) NMR data on the side-chain nitrogen in K6 gave rise to three different models of protein-surface interaction to explain the experimental data acquired. In this work we extend the analysis of the REDOR data by examining the contribution of interactions between surface phosphorus atoms to the observed N-15 REDOR decay. We performed P-31-P-31 recoupling experiments in HAP and (NH4)(2)HPO4 (DHP) to explore the nature of dipolar coupled P-31 spin networks. These studies indicate that extensive networks of dipolar coupled P-31 spins can be represented as stronger effective dipolar couplings, the existence of which must be included in the analysis of REDOR data. We carried out N-15 {P-31} REDOR in the case of DHP to determine how the size of the dephasing spin network influences the interpretation of the REDOR data. Although use of an extended P-31 coupled spin network simulates the REDOR data well, a simplified P-31 dephasing system composed of two spins with a larger dipolar coupling also simulates the REDOR data and only perturbs the heteronuclear couplings very slightly. The P-31-P-31 dipolar couplings between phosphorus nuclei in HAP can be replaced by an effective dipolar interaction of 600 Hz between two P-31 spins. We incorporated this coupling and applied the above approach to reanalyze the N-15 {P-31} REDOR of the lysine side chain approaching the HAP surface and have refined the binding models proposed earlier. We obtain N-15-P-31 distances between 3.3 and 5 angstrom from these models that are indicative of the possibility of a lysine-phosphate hydrogen bond.