Determining 3D structure from molecular formula and isotopologue rotational spectra in natural abundance with reflection-equivariant diffusion

J Chem Phys. 2024 Mar 28;160(12):124115. doi: 10.1063/5.0196620.

Abstract

Structure determination is necessary to identify unknown organic molecules, such as those in natural products, forensic samples, the interstellar medium, and laboratory syntheses. Rotational spectroscopy enables structure determination by providing accurate 3D information about small organic molecules via their moments of inertia. Using these moments, Kraitchman analysis determines isotopic substitution coordinates, which are the unsigned |x|, |y|, |z| coordinates of all atoms with natural isotopic abundance, including carbon, nitrogen, and oxygen. While unsigned substitution coordinates can verify guesses of structures, the missing +/- signs make it challenging to determine the actual structure from the substitution coordinates alone. To tackle this inverse problem, we develop Kreed (Kraitchman REflection-Equivariant Diffusion), a generative diffusion model that infers a molecule's complete 3D structure from only its molecular formula, moments of inertia, and unsigned substitution coordinates of heavy atoms. Kreed's top-1 predictions identify the correct 3D structure with near-perfect accuracy on large simulated datasets when provided with substitution coordinates of all heavy atoms with natural isotopic abundance. Accuracy decreases as fewer substitution coordinates are provided, but is retained for smaller molecules. On a test set of experimentally measured substitution coordinates gathered from the literature, Kreed predicts the correct all-atom 3D structure in 25 of 33 cases, demonstrating experimental potential for de novo 3D structure determination with rotational spectroscopy.