Pattern-free generation and quantum mechanical scoring of ring-chain tautomers

Quantum chemical package Jaguar: A survey of recent developments and unique features

This paper is dedicated to the quantum chemical package Jaguar, which is commercial software developed and distributed by Schrödinger, Inc. We discuss Jaguar's scientific features that are relevant to chemical research as well as describe those aspects of the program that are pertinent to the user interface, the organization of the computer code, and its maintenance and testing. Among the scientific topics that feature prominently in this paper are the quantum chemical methods grounded in the pseudospectral approach. A number of multistep workflows dependent on Jaguar are covered: prediction of protonation equilibria in aqueous solutions (particularly calculations of tautomeric stability and pKa), reactivity predictions based on automated transition state search, assembly of Boltzmann-averaged spectra such as vibrational and electronic circular dichroism, as well as nuclear magnetic resonance. Discussed also are quantum chemical calculations that are oriented toward materials science applications, in particular, prediction of properties of optoelectronic materials and organic semiconductors, and molecular catalyst design. The topic of treatment of conformations inevitably comes up in real world research projects and is considered as part of all the workflows mentioned above. In addition, we examine the role of machine learning methods in quantum chemical calculations performed by Jaguar, from auxiliary functions that return the approximate calculation runtime in a user interface, to prediction of actual molecular properties. The current work is second in a series of reviews of Jaguar, the first having been published more than ten years ago. Thus, this paper serves as a rare milestone on the path that is being traversed by Jaguar's development in more than thirty years of its existence.

Large Computational Survey of Intrinsic Reactivity of Aromatic Carbon Atoms with Respect to a Model Aldehyde Oxidase

Aldehyde oxidase (AOX) and other related molybdenum-containing enzymes are known to oxidize the C-H bonds of aromatic rings. This process contributes to the metabolism of pharmaceutical compounds and, therefore, is of vital importance to drug pharmacokinetics. The present work describes an automated computational workflow and its use for the prediction of intrinsic reactivity of small aromatic molecules toward a minimal model of the active site of AOX. The workflow is based on quantum chemical transition state searches for the underlying single-step oxidation reaction, where the automated protocol includes identification of unique aromatic C-H bonds, creation of three-dimensional reactant and product complex geometries via a templating approach, search for a transition state, and validation of reaction end points. Conformational search on the reactants, products, and the transition states is performed. The automated procedure has been validated on previously reported transition state barriers and was used to evaluate the intrinsic reactivity of nearly three hundred heterocycles commonly found in approved drug molecules. The intrinsic reactivity of more than 1000 individual aromatic carbon sites is reported. Stereochemical and conformational aspects of the oxidation reaction, which have not been discussed in previous studies, are shown to play important roles in accurate modeling of the oxidation reaction. Observations on structural trends that determine the reactivity are provided and rationalized.