Comment on the Definition and Labeling of pK50.

We propose a more rigorous definition for the recently introduced concept of pK50. The value of pK50 should be associated not with a "functional group", as originally postulated, but instead with an atom. The proposed clarification is meant to improve the interpretation and labeling of pK50.

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

In contrast to the computational generation of conventional tautomers, the analogous operation that would produce ring-chain tautomers is rarely available in cheminformatics codes. This is partly due to the perceived unimportance of ring-chain tautomerism and partly because specialized algorithms are required to realize the non-local proton transfers that occur during ring-chain rearrangement. Nevertheless, for some types of organic compounds, including sugars, warfarin analogs, fluorescein dyes and some drug-like compounds, ring-chain tautomerism cannot be ignored. In this work, a novel ring-chain tautomer generation algorithm is presented. It differs from previously proposed solutions in that it does not rely on hard-coded patterns of proton migrations and bond rearrangements, and should therefore be more general and maintainable. We deploy this algorithm as part of a workflow which provides an automated solution for tautomer generation and scoring. The workflow identifies protonatable and deprotonatable sites in the molecule using a previously described approach based on rapid micro-pKa prediction. These data are used to distribute the active protons among the protonatable sites exhaustively, at which point alternate resonance structures are considered to obtain pairs of atoms with opposite formal charge. These pairs are connected with a single bond and a 3D undistorted geometry is generated. The scoring of the generated tautomers is performed with a subsequent density functional theory calculation employing an implicit solvent model. We demonstrate the performance of our workflow on several types of organic molecules known to exist in ring-chain tautomeric equilibria in solution. In particular, we show that some ring-chain tautomers not found using previously published algorithms are successfully located by ours.

Generation of Tautomers Using Micro-p Ka's.

Solutions of organic molecules containing one or more heterocycles with conjugated bonds may exist as a mixture of tautomers, but typically only a few of them are significantly populated even though the potential number grows combinatorially with the number of protonation and deprotonation sites. Generating the most stable tautomers from a given input structure is an important and challenging task, and numerous algorithms to tackle it have been proposed in the literature. This work describes a novel approach for tautomer prediction that involves the combined use of molecular mechanics, semiempirical quantum chemistry, and density functional theory. The key idea in our method is to identify the protonation and deprotonation sites using estimated micro-p Ka's for every atom in the molecule as well as in its nearest protonated and deprotonated forms. To generate tautomers in a systematic way with minimal bias, we then consider the full set of tautomers that arise from the combinatorial distribution of all such mobile protons among all protonatable sites, with efficient postprocessing to screen away high-energy species. To estimate the micro-p Ka's, we present a new method designed for the current task, but we emphasize that any alternative method can be used in conjunction with our basic algorithm. Our approach is therefore grounded in the computational prediction of physical properties in aqueous solution, in contrast to other approaches that may rely on the use of hard-coded rules of proton distribution, previously observed tautomerization patterns from a known chemical space, or human input. We present examples of the application of our algorithm to organic and drug-like molecules, with a focus on novel structures where traditional methods are expected to perform worse.

Weighted Averaging Scheme and Local Atomic Descriptor for pKa Prediction Based on Density Functional Theory.

As a continuation of our work on developing a density functional theory-based pKa predictor, we present conceptual improvements to our previously published shell model, which is a hierarchical organization of pKa training sets and which, in principle, covers all chemical space. The improvements concern the way the studied chemical compound is associated with the data points from the training sets. By introducing a new descriptor of the local atomic environment which foregoes dependence on chemical bonding and connectivity, we are able to automatically locate molecules from the training set that are most relevant to the proton dissociation equilibrium under study. This new scheme leads to the prediction of a single pKa value weighted across multiple training sets and thus patches a defect disclosed in the formulation of our previous model. Using the new parametrization approach, the pKa prediction gets rid of outliers reported in previous applications of our approach, eliminates ambiguity in interpreting the results, and improves the overall accuracy. Our new treatment accounts for multiple conformations both on the level of energetics and parametrization. Illustrative results are shown for several types of chemical structures containing guanidine, amidine, amine, and phenol functional groups, and which are representative of practically important large and flexible drug-like molecules. Our method's performance is compared to the performance of other previously published pKa prediction methods. Further possible improvements to the organization of the training sets and the potential application of our new local atomic descriptor to other kinds of parametrizations are discussed.

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.

Discovery of a Novel Class of d-Amino Acid Oxidase Inhibitors Using the Schrödinger Computational Platform.

d-Serine is a coagonist of the N-methyl d-aspartate (NMDA) receptor, a key excitatory neurotransmitter receptor. In the brain, d-serine is synthesized from its l-isomer by serine racemase and is metabolized by the D-amino acid oxidase (DAO, DAAO). Many studies have linked decreased d-serine concentration and/or increased DAO expression and enzyme activity to NMDA dysfunction and schizophrenia. Thus, it is feasible to employ DAO inhibitors for the treatment of schizophrenia and other indications. Powered by the Schrödinger computational modeling platform, we initiated a research program to identify novel DAO inhibitors with the best-in-class properties. The program execution leveraged an hDAO FEP+ model to prospectively predict compound potency. A new class of DAO inhibitors with desirable properties has been discovered from this endeavor. Our modeling technology on this program has not only enhanced the efficiency of structure-activity relationship development but also helped to identify a previously unexplored subpocket for further optimization.

Accurate Quantum Chemical Calculation of Ionization Potentials: Validation of the DFT-LOC Approach via a Large Data Set Obtained from Experiments and Benchmark Quantum Chemical Calculations.

Density functional theory (DFT) is known to often fail when calculating thermodynamic values, such as ionization potentials (IPs), due to nondynamical error (i.e., the self-interaction term). Localized orbital corrections (LOCs), derived from assigning corresponding corrections for the atomic orbitals, bonds, and paired and unpaired electrons, are utilized to correct the IPs calculated from DFT. Some of the assigned parameters, which are physically due to the contraction of and change of the environment around a bond, depend on identifying the location in the molecule from which the electron is removed using differences in the charge density between neutral and oxidized species. In our training set, various small organic and inorganic molecules from the literature with the reported experimental IP were collected using the NIST database. For certain molecules with uncertain or no experimental measurements, we obtain the IP using coupled cluster theory and auxiliary field quantum Monte Carlo. After applying these corrections, as generated by least-squares regression, LOC reduces the mean absolute deviation (MAD) of the training set from 0.143 to 0.046 eV (R2 = 0.895), and LOC reduces the MAD of the test set from 0.192 to 0.097 eV (R2 = 0.833).

Empirical Conversion of pK a Values between Different Solvents and Interpretation of the Parameters: Application to Water, Acetonitrile, Dimethyl Sulfoxide, and Methanol.

An empirical conversion method (ECM) that transforms pK a values of arbitrary organic compounds from one solvent to the other is introduced. We demonstrate the method's usefulness and performance on pK a conversions involving water and organic solvents acetonitrile (MeCN), dimethyl sulfoxide (Me2SO), and methanol (MeOH). We focus on the pK a conversion from the known reference value in water to the other three organic solvents, although such a conversion can also be performed between any pair of the considered solvents. The ECM works with an additive parameter that is specific to a solvent and a molecular family (essentially characterized by a functional group that is titrated). We formally show that the method can be formulated with a single additive parameter, and that the extra multiplicative parameter used in other works is not required. The values of the additive parameter are determined from known pK a data, and their interpretation is provided on the basis of physicochemical concepts. The data set of known pK a values is augmented with pK a values computed with the recently introduced electrostatic transform method, whose validity is demonstrated. For a validation of our method, we consider pK a conversions for two data sets of titratable compounds. The first data set involves 81 relatively small molecules belonging to 19 different molecular families, with the pK a data available in all four considered solvents. The second data set involves 76 titratable molecules from 5 additional molecular families. These molecules are typically larger, and their experimental pK a values are available only in Me2SO and water. The validation tests show that the agreement between the experimental pK a data and the ECM predictions is generally good, with absolute errors often on the order of 0.5 pH units. The presence of a few outliers is rationalized, and observed trends with respect to molecular families are discussed.

Automated Transition State Search and Its Application to Diverse Types of Organic Reactions.

Transition state search is at the center of multiple types of computational chemical predictions related to mechanistic investigations, reactivity and regioselectivity predictions, and catalyst design. The process of finding transition states in practice is, however, a laborious multistep operation that requires significant user involvement. Here, we report a highly automated workflow designed to locate transition states for a given elementary reaction with minimal setup overhead. The only essential inputs required from the user are the structures of the separated reactants and products. The seamless workflow combining computational technologies from the fields of cheminformatics, molecular mechanics, and quantum chemistry automatically finds the most probable correspondence between the atoms in the reactants and the products, generates a transition state guess, launches a transition state search through a combined approach involving the relaxing string method and the quadratic synchronous transit, and finally validates the transition state via the analysis of the reactive chemical bonds and imaginary vibrational frequencies as well as by the intrinsic reaction coordinate method. Our approach does not target any specific reaction type, nor does it depend on training data; instead, it is meant to be of general applicability for a wide variety of reaction types. The workflow is highly flexible, permitting modifications such as a choice of accuracy, level of theory, basis set, or solvation treatment. Successfully located transition states can be used for setting up transition state guesses in related reactions, saving computational time and increasing the probability of success. The utility and performance of the method are demonstrated in applications to transition state searches in reactions typical for organic chemistry, medicinal chemistry, and homogeneous catalysis research. In particular, applications of our code to Michael additions, hydrogen abstractions, Diels-Alder cycloadditions, carbene insertions, and an enzyme reaction model involving a molybdenum complex are shown and discussed.

Multiconformation, Density Functional Theory-Based pKa Prediction in Application to Large, Flexible Organic Molecules with Diverse Functional Groups.

We consider the conformational flexibility of molecules and its implications for micro- and macro-pKa. The corresponding formulas are derived and discussed against the background of a comprehensive scientific and algorithmic description of the latest version of our computer program Jaguar pKa, a density functional theory-based pKa predictor, which is now capable of acting on multiple conformations explicitly. Jaguar pKa is essentially a complex computational workflow incorporating research and technologies from the fields of cheminformatics, molecular mechanics, quantum mechanics, and implicit solvation models. The workflow also makes use of automatically applied empirical corrections which account for the systematic errors resulting from the neglect of explicit solvent interactions in the algorithm's implicit solvent model. Applications of our program to large, flexible organic molecules representing several classes of functional groups are shown, with a particular emphasis in illustrations laid on drug-like molecules. It is demonstrated that a combination of aggressive conformational search and an explicit consideration of multiple conformations nearly eliminates the dependence of results on the initially chosen conformation. In certain cases this leads to unprecedented accuracy, which is sufficient for distinguishing stereoisomers that have slightly different pKa values. An application of Jaguar pKa to proton sponges, the pKa of which are strongly influenced by steric effects, showcases the advantages that pKa predictors based on quantum mechanical calculations have over similar empirical programs.