Multi-objective goal-directed optimization of de novo stable organic radicals for aqueous redox flow batteries

Sowndarya, S. S. V.; Law, J.; Tripp, C.; Duplyakin, D.; Skordilis, E.; Biagioni, D.; Paton, R. S.; St. John, P. C. Nat. Mach. Intell. 2022, 7, 720–730

Mechanistic Studies Yield Improved Protocols for Base-Catalyzed anti-Markovnikov Alcohol Addition Reactions

Luo, C.; Alegre-Requena, J. V.; Sujansky, S. J.; Pajk, S.; Gallegos, L. C.; Paton, R. S.; Bandar, J. S. J. Am. Chem. Soc.  2022, 144, 22, 9586–9596

Homologation of Electron-Rich Benzyl Bromide Derivatives via Diazo C–C Bond Insertion

Modak, A.; Alegre-Requena, J. V.; Lescure, L.; Rynders, K. J.; Paton, R. S.; Race, N. J. Am. Chem. Soc. 2022, 144, 86–92


CASCADE stands for ChemicAShift CAlculation with DEep learning. It is a stereochemically-aware graph network for the prediction of NMR chemical shifts. Model training was performed against 8,000 DFT structures followed by transfer learning with experimental  spectra. A web-server has been created to access CASCADE predictions from SMILES or by drawing structures in the graphical interface. An automated workflow executes 3D structure embedding and MMFF conformer searching. The full ensemble of optimized conformations are passed to a trained graph neural network to predict the NMR chemical shifts (in ppm) for C and H atoms. The underlying datasets used for training and the Python code to run CASCADE from the command line have been made available.



DBSTEP is a python package for obtaining DFT-Based Steric Parameters from 3-dimensional chemical structures. It can parse the outputs from most computational chemistry programs and other common molecular structure file formats. Steric properties can either be obtained exactly or by using a Cartesian grid, the latter approach being amenable to the featurization of a molecular isodensity surface (DBSTEP can process wavefunction files) rather than using classical atomic radii. Currently,  traditional Sterimol parameters (L, Bmin, Bmax) and percent buried volume parameters are implemented, as well as  our novel steric parameter vectors Sterimol2vec and vol2vec. This package is designed for use on the command line or alternatively implemented in a Python script for use in a computational workflow to collect steric parameters.

Importance of Engineered and Learned Molecular Representations in Predicting Organic Reactivity, Selectivity, and Chemical Properties

Gallegos, L. C.; Luchini, G.; St. John, P. C.; Kim, S.; Paton, R. S. Acc. Chem. Res. 2021, 54, 827–836



A Python program to compute quasi-harmonic thermochemical data and potential energy surface diagrams from frequency calculations at a given temperature/concentration, corrected for the effects of vibrational scaling-factors. All (electronic, translational, rotational and vibrational) partition functions are recomputed and can be correct to any temperature or concentration. The first public version of GoodVibes was released in 2016 and it has undergone several revisions since, during which time it has been used by many groups around the world. The program is described in the publication: GoodVibes: automated thermochemistry for heterogeneous computational chemistry data



A program to generate Boltzmann-weighted Sterimol Steric Parameters for conformationally-flexible substituents that integrates with PyMol. The program contains an automated computational workflow which computes multidimensional Sterimol parameters. For flexible molecules or substituents, the program will generate & optimize a conformational ensemble, and produce Boltzmann-weighted Sterimol parameters. It has been developed as a PyMol plugin and can be run from within the graphical user interface. The wSterimol code is described in more detail in Conformational Effects on Physical-Organic Descriptors – the Case of Sterimol Steric Parameters