Synthesis and Structural Characterization of Trivalent Amino Acid Derived Chiral Phosphorus Compounds

Improving the accuracy of 31P NMR chemical shift calculations by use of scaling methods

Calculation of ³¹P NMR chemical shifts for a series of tri- and tetracoordinate phosphorus compounds using several basis sets and density functional theory (DFT) functionals gave a modest fit to experimental chemical shifts, but an excellent linear fit when plotted against the experimental values. The resultant scaling methods were then applied to a variety of "large" compounds previously selected by Latypov et al. and a set of stereoisomeric and unusual compounds selected here. No one method was best for all structural types. For compounds that contain P-P bonds and P-C multiple bonds, the Latypov et al. method using the PBE0 functional was best (mean absolute deviation/root mean square deviation (MAD/RMSD) = 6.9/8.5 ppm and 6.6/8.2 ppm, respectively), but for the full set of compounds gave higher deviations (MAD/RMSD = 8.2/12.3 ppm), and failed by over 60 ppm for a three-membered phosphorus heterocycle. Use of the M06-2X functional for both the structural optimization and NMR chemical shift calculation was best overall for the compounds without P-C multiple bonds (MAD/RMSD = 5.4/7.1 ppm), but failed by 30-49 ppm for compounds having any P-C multiple-bond character. Failures of these magnitudes have not been reported previously for these widely used functionals. These failures were then used to screen a variety of recommended functionals, leading to better overall methods for calculation of these chemical shifts: optimization with the M06-2X functional and NMR calculation with the PBE0 or ωB97x-D functionals gave values for MAD/RMSD = 6.9/8.5 ppm and 6.8/9.1 ppm, respectively, over an experimental chemical shift range of -181 to 356 ppm. Due to the unexplained failures observed, we recommend use of more than one method when looking at novel structures.

A non-Karplus effect: evidence from phosphorus heterocycles and DFT calculations of the dependence of vicinal phosphorus-hydrogen NMR coupling constants on lone-pair conformation

In contrast to literature reports of a Karplus-type curve that correlates (3)J(PH) with phosphorus-hydrogen dihedral angle, a recently reported glycine-derived 1,3,2-oxazaphospholidine (7c) has two hydrogen atoms on the ring with identical PNCH dihedral angles but measured coupling constants of ∼6 and 1.5 Hz. DFT calculations were in accord with these values and suggested that the smaller coupling constant is negative. Experimental evidence of the opposite signs of these coupling constants was obtained by analysis of the ABX NMR spectrum of the new glycine-derived N-p-toluenesulfonyl phosphorus heterocycle 6c. DFT calculations on 6c and on Me(2)NPCl(2) and t-BuPCl(2) were also in accord with NMR data and allowed confirmation of unusual features including a lone pair effect on (3)J(PH), the negative coupling constant, temperature-dependent chemical shifts due to rotation about the sulfonamide S-N bond, and vicinal phosphorus-hydrogen coupling constants over 40 Hz. Calculation of phosphorus-hydrogen coupling constants both as a function of PYCH dihedral angle θ (Y = O, N, C) and lone pair-PYC dihedral angle ω shows similar θ,ω surfaces for (3)J(PH) with a range of (3)J(PH) from -4.4 to +51 Hz and demonstrates the large non-Karplus effect of lone-pair conformation on vicinal phosphorus-hydrogen coupling constants.