Isotopic substitution with C-13 on the amide C=O has become an important means of determining localized structural information about peptide conformations with vibrational spectroscopy. Various approaches to the modeling of the interactions between labeled amide sites, specifically for antiparallel two-stranded, P-forming peptides, were investigated, including different force fields [dipole-dipole interaction vs density functional theory (DFT) treatments], basis sets, and sizes of model peptides used for ab initio calculations, as well as employing models of solvation. For these beta-sheet systems the effect of the relative positions of the 13C isotopic labels in each strand on their infrared spectra was investigated. The results suggest that the interaction between labeled amide groups in different strands can be used as an indicator of local P-structure formation, because coupling between close-lying C=O groups on opposing chains leads to the largest frequency shifts, yet some alternate placements can lead to intensity enhancements. The basic character of the coupling interaction between labeled modes on opposing strands is independent of changes in peptide length, water solvent environment, twisting of the sheet structure, and basis set used in the calculations, although the absolute frequencies and detailed coupling magnitudes change under each of these perturbations. In particular, two strands of three amides each contain the basic interactions needed to simulate larger sheets, with the only exception that the C=O groups forming H-bonded rings at the termini can yield different coupling values than central ones of the same structure. Spectral frequencies and intensities were modeled ab initio by DFT primarily at the BPW91/ 6-31 G** level for pairs of three, four, and six amide strands. Comparison to predictions of a classical coupled oscillator model show qualitative but not quantitative agreement with these DFT results.