The overall and detailed elucidation (including the stereochemical aspects) of enzymatic mechanisms requires the access to all reliable information related to the natural isotopic fractionation of both precursors and products. Natural abundance deuterium (NAD) 2D-NMR experiments in polypeptide liquid-crystalline solutions are a new, suitable tool for analyzing site-specific deuterium isotopic distribution profiles. Here this method is utilized for analyzing saturated C14 to C18 fatty acid methyl esters (FAMEs), which are challenging because of the crowding of signals in a narrow spectral region. Experiments in achiral and chiral oriented solutions were performed. The spectral analysis is supplemented by the theoretical prediction of quadrupolar splittings as a function of the geometry and flexibility of FAMEs, based on a novel computational methodology. This allows us to confirm the spectral assignments, while providing insights into the mechanism of solute ordering in liquid-crystalline polypeptide solutions. This is found to be dominated by steric repulsions between FAMEs and polypeptides.