Structures and solvent effects on the 1H and 13C NMR chemical shifts of the family of three pro-lithospermate molecules: DFT study

The solvent effects on the structures and nuclear magnetic resonance spectroscopy of a series of three pro-lithospermate molecules are studied. These molecules include methyl pro-lithospermate, ethyl pro-lithospermate, and n-butyl pro-lithospermate, which exhibit enantiomers and diastereoisomers. The trans-oriented enantiomers were observed by Yu et al. [Nat. Prod. Res. 31(21), 2505-2512 (2017)] and are the subject of the present study. The density functional theory was validated using the gauge-including atomic orbital/B3LYP-D3BJ/6-31++G(d,p) and continuous set of gauge transformation (CSGT)/the long-range corrected Tao, Perdew, Staroverov, and Scuseria's tau-dependent gradient-corrected correlation-exchange functional (LC-TPSS)TPSS/cc-pVTZ methods. The chemical shifts of H-1 and C-13 were calculated in dimethylsulfoxide (DMSO), methanol, and water. It turns out that solvent effects are weak on the geometrical parameters (bond lengths and bond angles). The harmonic oscillator model of the aromaticity index was used to calculate the aromaticity of the subunit rings. With regard to the available experimental results, the chemical shifts of H-1 and C-13 are reproduced better by B3LYP-D3BJ/6-31++G(d,p) than by LC-TPSSTPSS/cc-pVTZ in the DMSO medium. The solvent effects are significant on the C-13 chemical shift, mostly for protic solvents, and insignificant on H-1. The vicinal spin-spin coupling constants of protons were calculated. The coupling constant of trans-olefinic and ortho-aromatic protons fits well with the experiment. (c) 2023 Author(s) All article content, except where otherwise noted, is licensed under a Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

Automated summary

This study investigates the effects of solvents on the structures and nuclear magnetic resonance spectroscopy of three pro-lithospermate molecules. The solvent effects on the geometrical parameters are found to be weak, and better results are obtained by using the B3LYP-D3BJ/6-31++G(d,p) method in DMSO for the calculation of chemical shifts.