Pegylated Phosphine Ligands in Iridium(I) Catalyzed Hydrogen Isotope Exchange Reactions in Aqueous Buffers.

The synthesis of a water-soluble, phosphine-pegylated iridium(I) catalyst and its application in hydrogen isotope exchange (HIE) reactions in buffer is reported. The longer polyethylene glycol side chains on the phosphine increased the water solubility independently from the pH. HIE reactions of polar substrates in protic solvents were studied. DFT calculations gave further insides into the catalytic processes. The scope and limitation of the pegylated catalyst was studied in HIE reactions of several complex compounds in borax buffer at pH 9 and the best conditions were applied in a tritium experiment with the drug telmisartan.

Hydrogen Isotope Exchange by Homogeneous Iridium Catalysis in Aqueous Buffers with Deuterium or Tritium Gas.

We have studied the highly selective homogeneous iridium-catalyzed hydrogen isotope exchange (HIE) with deuterium or tritium gas as an isotope source in water and buffers. With an improved water-soluble Kerr-type catalyst, we have achieved the first insight into applying HIE reactions in aqueous media with varying pH. Density functional theory (DFT) calculations gave consistent insights in the calculated energies of transition states and coordination complexes, further explaining the observed reactivity and guidance on the scope and limitations for HIE reactions in water. Finally, we successfully adapted these findings to tritium chemistry.

SAMPL7 physical property prediction from EC-RISM theory.

Inspired by the successful application of the embedded cluster reference interaction site model (EC-RISM), a combination of quantum-mechanical calculations with three-dimensional RISM theory to predict Gibbs energies of species in solution within the SAMPL6.1 (acidity constants, pKa) and SAMPL6.2 (octanol-water partition coefficients, log P) the methodology was applied to the recent SAMPL7 physical property challenge on aqueous pKa and octanol-water log P values. Not part of the challenge but provided by the organizers, we also computed distribution coefficients log D7.4 from predicted pKa and log P data. While macroscopic pKa predictions compared very favorably with experimental data (root mean square error, RMSE 0.72 pK units), the performance of the log P model (RMSE 1.84) fell behind expectations from the SAMPL6.2 challenge, leading to reasonable log D7.4 predictions (RMSE 1.69) from combining the independent calculations. In the post-submission phase, conformations generated by different methodology yielded results that did not significantly improve the original predictions. While overall satisfactory compared to previous log D challenges, the predicted data suggest that further effort is needed for optimizing the robustness of the partition coefficient model within EC-RISM calculations and for shaping the agreement between experimental conditions and the corresponding model description.

Quantum-mechanical property prediction of solvated drug molecules: what have we learned from a decade of SAMPL blind prediction challenges?

Joint academic-industrial projects supporting drug discovery are frequently pursued to deploy and benchmark cutting-edge methodical developments from academia in a real-world industrial environment at different scales. The dimensionality of tasks ranges from small molecule physicochemical property assessment over protein-ligand interaction up to statistical analyses of biological data. This way, method development and usability both benefit from insights gained at both ends, when predictiveness and readiness of novel approaches are confirmed, but the pharmaceutical drug makers get early access to novel tools for the quality of drug products and benefit of patients. Quantum-mechanical and simulation methods particularly fall into this group of methods, as they require skills and expense in their development but also significant resources in their application, thus are comparatively slowly dripping into the realm of industrial use. Nevertheless, these physics-based methods are becoming more and more useful. Starting with a general overview of these and in particular quantum-mechanical methods for drug discovery we review a decade-long and ongoing collaboration between Sanofi and the Kast group focused on the application of the embedded cluster reference interaction site model (EC-RISM), a solvation model for quantum chemistry, to study small molecule chemistry in the context of joint participation in several SAMPL (Statistical Assessment of Modeling of Proteins and Ligands) blind prediction challenges. Starting with early application to tautomer equilibria in water (SAMPL2) the methodology was further developed to allow for challenge contributions related to predictions of distribution coefficients (SAMPL5) and acidity constants (SAMPL6) over the years. Particular emphasis is put on a frequently overlooked aspect of measuring the quality of models, namely the retrospective analysis of earlier datasets and predictions in light of more recent and advanced developments. We therefore demonstrate the performance of the current methodical state of the art as developed and optimized for the SAMPL6 pKa and octanol-water log P challenges when re-applied to the earlier SAMPL5 cyclohexane-water log D and SAMPL2 tautomer equilibria datasets. Systematic improvement is not consistently found throughout despite the similarity of the problem class, i.e. protonation reactions and phase distribution. Hence, it is possible to learn about hidden bias in model assessment, as results derived from more elaborate methods do not necessarily improve quantitative agreement. This indicates the role of chance or coincidence for model development on the one hand which allows for the identification of systematic error and opportunities toward improvement and reveals possible sources of experimental uncertainty on the other. These insights are particularly useful for further academia-industry collaborations, as both partners are then enabled to optimize both the computational and experimental settings for data generation.

The SAMPL6 challenge on predicting octanol-water partition coefficients from EC-RISM theory.

Results are reported for octanol-water partition coefficients (log P) of the neutral states of drug-like molecules provided during the SAMPL6 (Statistical Assessment of Modeling of Proteins and Ligands) blind prediction challenge from applying the "embedded cluster reference interaction site model" (EC-RISM) as a solvation model for quantum-chemical calculations. Following the strategy outlined during earlier SAMPL challenges we first train 1- and 2-parameter water-free ("dry") and water-saturated ("wet") models for n-octanol solvation Gibbs energies with respect to experimental values from the "Minnesota Solvation Database" (MNSOL), yielding a root mean square error (RMSE) of 1.5 kcal mol-1 for the best-performing 2-parameter wet model, while the optimal water model developed for the pKa part of the SAMPL6 challenge is kept unchanged (RMSE 1.6 kcal mol-1 for neutral compounds from a model trained on both neutral and ionic species). Applying these models to the blind prediction set yields a log P RMSE of less than 0.5 for our best model (2-parameters, wet). Further analysis of our results reveals that a single compound is responsible for most of the error, SM15, without which the RMSE drops to 0.2. Since this is the only compound in the challenge dataset with a hydroxyl group we investigate other alcohols for which Gibbs energy of solvation data for both water and n-octanol are available in the MNSOL database to demonstrate a systematic cause of error and to discuss strategies for improvement.

C-H Functionalization-Prediction of Selectivity in Iridium(I)-Catalyzed Hydrogen Isotope Exchange Competition Reactions.

An assessment of the C-H activation catalyst [(COD)Ir(IMes)(PPh3 )]PF6 (COD=1,5-cyclooctadiene, IMes=1,3-bis(2,4,6-trimethylphenyl)imidazol-2-ylidene) in the deuteration of phenyl rings containing different functional directing groups is divulged. Competition experiments have revealed a clear order of the directing groups in the hydrogen isotope exchange (HIE) with an iridium (I) catalyst. Through DFT calculations the iridium-substrate coordination complex has been identified to be the main trigger for reactivity and selectivity in the competition situation with two or more directing groups. We postulate that the competition concept found in this HIE reaction can be used to explain regioselectivities in other transition-metal-catalyzed functionalization reactions of complex drug-type molecules as long as a C-H activation mechanism is involved.

Computational Investigation of Drug Phototoxicity: Photosafety Assessment, Photo-Toxophore Identification, and Machine Learning.

One of the most appreciated capabilities of computational toxicology is to support the design of pharmaceuticals with reduced toxicological hazard. To this end, we have strengthened our drug photosafety assessments by applying novel computer models for the anticipation of in vitro phototoxicity and human photosensitization. These models are typically used in pharmaceutical discovery projects as part of the compound toxicity assessments and compound optimization methods. To ensure good data quality and aiming at models with global applicability we separately compiled and curated highly chemically diverse data sets from 3T3 NRU phototoxicity reports (450 compounds) and clinical photosensitization alerts (1419 compounds) which are provided as supplements. The latter data gives rise to a comprehensive list of explanatory fragments for visual guidance, termed phototoxophores, by application of a Bayesian statistics approach. To extend beyond the domain of well sampled fragments we applied machine learning techniques based on explanatory descriptors such as pharmacophoric fingerprints or, more important, accurate electronic energy descriptors. Electronic descriptors were extracted from quantum chemical computations at the density functional theory (DFT) level. Accurate UV/vis spectral absorption descriptors and pharmacophoric fingerprints turned out to be necessary for predictive computer models, which were both derived from Deep Neural Networks but also the simpler Random Decision Forests approach. Model accuracies of 83-85% could typically be reached for diverse test data sets and other company in-house data, while model sensitivity (the capability of correctly detecting toxicants) was even better, reaching 86%-90%. Importantly, a computer model-triggered response-map allowed for graphical/chemical interpretability also in the case of previously unknown phototoxophores. The photosafety models described here are currently applied in a prospective manner for the hazard identification, prioritization, and optimization of newly designed molecules.

Characterizing hydration sites in protein-ligand complexes towards the design of novel ligands.

Water is an essential part of protein binding sites and mediates interactions to ligands. Its displacement by ligand parts affects the free binding energy of resulting protein-ligand complexes. Therefore the characterization of solvation properties is important for design. Of particular interest is the propensity of localized water to be favorably displaced by a ligand. This review discusses two popular computational approaches addressing these questions, namely WaterMap based on statistical mechanics analysis of MD simulations and 3D RISM based on integral equation theory of liquids. The theoretical background and recent applications in structure-based design will be presented.

The SAMPL6 challenge on predicting aqueous pKa values from EC-RISM theory.

The "embedded cluster reference interaction site model" (EC-RISM) integral equation theory is applied to the problem of predicting aqueous pKa values for drug-like molecules based on an ensemble of tautomers. EC-RISM is based on self-consistent calculations of a solute's electronic structure and the distribution function of surrounding water. Following-up on the workflow developed after the SAMPL5 challenge on cyclohexane-water distribution coefficients we extended and improved the methodology by taking into account exact electrostatic solute-solvent interactions taken from the wave function in solution. As before, the model is calibrated against Gibbs energies of hydration from the "Minnesota Solvation Database" and a public dataset of acidity constants of organic acids and bases by adjusting in total 4 parameters, among which only 3 are relevant for predicting pKa values. While the best-performing training model yields a root-mean-square error (RMSE) of 1 pK unit, the corresponding test set prediction on the full SAMPL6 dataset of macroscopic pKa values using the same level of theory exhibits slightly larger error (1.7 pK units) than the best test set model submitted (1.7 pK units for corresponding training set vs. test set performance of 1.6). Post-submission analysis revealed a number of physical optimization options regarding the numerical treatment of electrostatic interactions and conformational sampling. While the experimental test set data revealed after submission was not used for reparametrizing the methodology, the best physically optimized models consequentially result in RMSEs of 1.5 if only improved electrostatic interactions are considered and of 1.1 if, in addition, conformational sampling accounts for quantum-chemically derived rankings. We conclude that these numbers are probably near the ultimate accuracy achievable with the simple 3-parameter model using a single or the two best-ranking conformations per tautomer or microstate. Finally, relations of the present macrostate approach to microstate pKa results are discussed and some illustrative results for microstate populations are presented.

Highly Selective Directed Iridium-Catalyzed Hydrogen Isotope Exchange Reactions of Aliphatic Amides.

For the first time, we describe highly selective homogeneous iridium-catalyzed hydrogen isotope exchange (HIE) of unactivated C(sp3 ) centers in aliphatic amides. When using the commercially available Kerr catalyst, the HIE with a series of common antibody-drug conjugate (ADC) linker side chains proceeds with high yields, high regioselectivity, and with deuterium incorporation up to 99 %. The method is fully translatable to the specific requirements of tritium chemistry and its effectiveness was demonstrated by direct tritium labelling of a maytansinoid. The scope of the method can be extended to simple amino acids, with high HIE activity observed for glycine and alanine. In di- and tripeptides, a very interesting protecting-group-dependent tunable selectivity was observed. DFT calculations gave insight into the energies of the transition states, thereby explaining the observed selectivity and the influence of the amino acid protecting groups.

Thermodynamic Characterization of Hydration Sites from Integral Equation-Derived Free Energy Densities: Application to Protein Binding Sites and Ligand Series.

Water molecules play an essential role for mediating interactions between ligands and protein binding sites. Displacement of specific water molecules can favorably modulate the free energy of binding of protein-ligand complexes. Here, the nature of water interactions in protein binding sites is investigated by 3D RISM (three-dimensional reference interaction site model) integral equation theory to understand and exploit local thermodynamic features of water molecules by ranking their possible displacement in structure-based design. Unlike molecular dynamics-based approaches, 3D RISM theory allows for fast and noise-free calculations using the same detailed level of solute-solvent interaction description. Here we correlate molecular water entities instead of mere site density maxima with local contributions to the solvation free energy using novel algorithms. Distinct water molecules and hydration sites are investigated in multiple protein-ligand X-ray structures, namely streptavidin, factor Xa, and factor VIIa, based on 3D RISM-derived free energy density fields. Our approach allows the semiquantitative assessment of whether a given structural water molecule can potentially be targeted for replacement in structure-based design. Finally, PLS-based regression models from free energy density fields used within a 3D-QSAR approach (CARMa - comparative analysis of 3D RISM Maps) are shown to be able to extract relevant information for the interpretation of structure-activity relationship (SAR) trends, as demonstrated for a series of serine protease inhibitors.

The hpCADD NDDO Hamiltonian: Parametrization.

A neglect of diatomic differential overlap (NDDO) Hamiltonian has been parametrized as an electronic component of a polarizable force field. Coulomb and exchange potentials derived directly from the NDDO Hamiltonian in principle can be used with classical potentials, thus forming the basis for a new generation of efficiently applicable multipolar polarizable force fields. The new hpCADD Hamiltonian uses force-field-like atom types and reproduces the electrostatic properties (dipole moment, molecular electrostatic potential) and Koopmans' theorem ionization potentials closely, as demonstrated for a large training set and an independent test set of small molecules. The Hamiltonian is not intended to reproduce geometries or total energies well, as these will be controlled by the classical force-field potentials. In order to establish the hpCADD Hamiltonian as an electronic component in force-field-based calculations, we tested its performance in combination with the 3D reference interaction site model (3D RISM) for aqueous solutions. Comparison of the resulting solvation free energies for the training and test sets to atomic charges derived from standard procedures, exact solute-solvent electrostatics based on high-level quantum-chemical reference data, and established semiempirical Hamiltonians demonstrates the advantages of the hpCADD parametrization.

Characterization of RA839, a Noncovalent Small Molecule Binder to Keap1 and Selective Activator of Nrf2 Signaling.

The activation of the transcription factor NF-E2-related factor 2 (Nrf2) maintains cellular homeostasis in response to oxidative stress by the regulation of multiple cytoprotective genes. Without stressors, the activity of Nrf2 is inhibited by its interaction with the Keap1 (kelch-like ECH-associated protein 1). Here, we describe (3S)-1-[4-[(2,3,5,6-tetramethylphenyl) sulfonylamino]-1-naphthyl]pyrrolidine-3-carboxylic acid (RA839), a small molecule that binds noncovalently to the Nrf2-interacting kelch domain of Keap1 with a Kd of ∼6 µM, as demonstrated by x-ray co-crystallization and isothermal titration calorimetry. Whole genome DNA arrays showed that at 10 µM RA839 significantly regulated 105 probe sets in bone marrow-derived macrophages. Canonical pathway mapping of these probe sets revealed an activation of pathways linked with Nrf2 signaling. These pathways were also activated after the activation of Nrf2 by the silencing of Keap1 expression. RA839 regulated only two genes in Nrf2 knock-out macrophages. Similar to the activation of Nrf2 by either silencing of Keap1 expression or by the reactive compound 2-cyano-3,12-dioxooleana-1,9-dien-28-oic acid methyl ester (CDDO-Me), RA839 prevented the induction of both inducible nitric-oxide synthase expression and nitric oxide release in response to lipopolysaccharides in macrophages. In mice, RA839 acutely induced Nrf2 target gene expression in liver. RA839 is a selective inhibitor of the Keap1/Nrf2 interaction and a useful tool compound to study the biology of Nrf2.

The SAMPL5 challenge for embedded-cluster integral equation theory: solvation free energies, aqueous pK a, and cyclohexane-water log D.

We predict cyclohexane-water distribution coefficients (log D 7.4) for drug-like molecules taken from the SAMPL5 blind prediction challenge by the "embedded cluster reference interaction site model" (EC-RISM) integral equation theory. This task involves the coupled problem of predicting both partition coefficients (log P) of neutral species between the solvents and aqueous acidity constants (pK a) in order to account for a change of protonation states. The first issue is addressed by calibrating an EC-RISM-based model for solvation free energies derived from the "Minnesota Solvation Database" (MNSOL) for both water and cyclohexane utilizing a correction based on the partial molar volume, yielding a root mean square error (RMSE) of 2.4 kcal mol-1 for water and 0.8-0.9 kcal mol-1 for cyclohexane depending on the parametrization. The second one is treated by employing on one hand an empirical pK a model (MoKa) and, on the other hand, an EC-RISM-derived regression of published acidity constants (RMSE of 1.5 for a single model covering acids and bases). In total, at most 8 adjustable parameters are necessary (2-3 for each solvent and two for the pK a) for training solvation and acidity models. Applying the final models to the log D 7.4 dataset corresponds to evaluating an independent test set comprising other, composite observables, yielding, for different cyclohexane parametrizations, 2.0-2.1 for the RMSE with the first and 2.2-2.8 with the combined first and second SAMPL5 data set batches. Notably, a pure log P model (assuming neutral species only) performs statistically similarly for these particular compounds. The nature of the approximations and possible perspectives for future developments are discussed.

1,1-Dioxo-5,6-dihydro-[4,1,2]oxathiazines, a novel class of 11ß-HSD1 inhibitors for the treatment of diabetes.

Racemic cis-1,1-dioxo-5,6-dihydro-[4,1,2]oxathiazine derivative 4a was isolated as an impurity in a sample of a hit from a HTS campaign on 11ß-hydroxysteroid dehydrogenase type 1 (11ß-HSD1). After separation by chiral chromatography the 4a-S, 8a-R enantiomer of compound 4a was identified as the true, potent enzyme inhibitor. The cocrystal structure of 4a with human and murine 11ß-HSD1 revealed the unique binding mode of the oxathiazine series. SAR elucidation and optimization in regard to metabolic stability led to monocyclic tetramethyloxathiazines as exemplified by compound 21g.

Discovery of SAR184841, a potent and long-lasting inhibitor of 11ß-hydroxysteroid dehydrogenase type 1, active in a physiopathological animal model of T2D.

Starting from 11ß-HSD1 inhibitors that were active ex vivo but with Cyp 3A4 liability, we obtained a new series of adamantane ureas displaying potent inhibition of both human and rodent 11ß-HSD1 enzymes, devoid of Cyp 3A4 interactions, and rationally designed to provide long-lasting inhibition in target tissues. Final optimizations lead to SAR184841 with good oral pharmacokinetic properties showing in vivo activity and improvement of metabolic parameters in a physiopathological model of type 2 diabetes.

Quantum mechanics-based properties for 3D-QSAR.

We have used a set of four local properties based on semiempirical molecular orbital calculations (electron density (ρ), hydrogen bond donor field (HDF), hydrogen bond acceptor field (HAF), and molecular lipophilicity potential (MLP)) for 3D-QSAR studies to overcome the limitations of the current force field-based molecular interaction fields (MIFs). These properties can be calculated rapidly and are thus amenable to high-throughput industrial applications. Their statistical performance was compared with that of conventional 3D-QSAR approaches using nine data sets (angiotensin converting enzyme inhibitors (ACE), acetylcholinesterase inhibitors (AchE), benzodiazepine receptor ligands (BZR), cyclooxygenase-2 inhibitors (COX2), dihydrofolate reductase inhibitors (DHFR), glycogen phosphorylase b inhibitors (GPB), thermolysin inhibitors (THER), thrombin inhibitors (THR), and serine protease factor Xa inhibitors (fXa)). The 3D-QSAR models generated were tested thoroughly for robustness and predictive ability. The average performance of the quantum mechanical molecular interaction field (QM-MIF) models for the nine data sets is better than that of the conventional force field-based MIFs. In the individual data sets, the QM-MIF models always perform better than, or as well as, the conventional approaches. It is particularly encouraging that the relative performance of the QM-MIF models improves in the external validation. In addition, the models generated showed statistical stability with respect to model building procedure variations such as grid spacing size and grid orientation. QM-MIF contour maps reproduce the features important for ligand binding for the example data set (factor Xa inhibitors), demonstrating the intuitive chemical interpretability of QM-MIFs.

Impact of Applicability Domains to Generative Artificial Intelligence.

Molecular generative artificial intelligence is drawing significant attention in the drug design community, with several experimentally validated proof of concepts already published. Nevertheless, generative models are known for sometimes generating unrealistic, unstable, unsynthesizable, or uninteresting structures. This calls for methods to constrain those algorithms to generate structures in drug-like portions of the chemical space. While the concept of applicability domains for predictive models is well studied, its counterpart for generative models is not yet well-defined. In this work, we empirically examine various possibilities and propose applicability domains suited for generative models. Using both public and internal data sets, we use generative methods to generate novel structures that are predicted to be actives by a corresponding quantitative structure-activity relationships model while constraining the generative model to stay within a given applicability domain. Our work looks at several applicability domain definitions, combining various criteria, such as structural similarity to the training set, similarity of physicochemical properties, unwanted substructures, and quantitative estimate of drug-likeness. We assess the structures generated from both qualitative and quantitative points of view and find that the applicability domain definitions have a strong influence on the drug-likeness of generated molecules. An extensive analysis of our results allows us to identify applicability domain definitions that are best suited for generating drug-like molecules with generative models. We anticipate that this work will help foster the adoption of generative models in an industrial context.

3D-QSAR based on quantum-chemical molecular fields: toward an improved description of halogen interactions.

Current 3D-QSAR methods such as CoMFA or CoMSIA make use of classical force-field approaches for calculating molecular fields. Thus, they can not adequately account for noncovalent interactions involving halogen atoms like halogen bonds or halogen-π interactions. These deficiencies in the underlying force fields result from the lack of treatment of the anisotropy of the electron density distribution of those atoms, known as the "σ-hole", although recent developments have begun to take specific interactions such as halogen bonding into account. We have now replaced classical force field derived molecular fields by local properties such as the local ionization energy, local electron affinity, or local polarizability, calculated using quantum-mechanical (QM) techniques that do not suffer from the above limitation for 3D-QSAR. We first investigate the characteristics of QM-based local property fields to show that they are suitable for statistical analyses after suitable pretreatment. We then analyze these property fields with partial least-squares (PLS) regression to predict biological affinities of two data sets comprising factor Xa and GABA-A/benzodiazepine receptor ligands. While the resulting models perform equally well or even slightly better in terms of consistency and predictivity than the classical CoMFA fields, the most important aspect of these augmented field-types is that the chemical interpretation of resulting QM-based property field models reveals unique SAR trends driven by electrostatic and polarizability effects, which cannot be extracted directly from CoMFA electrostatic maps. Within the factor Xa set, the interaction of chlorine and bromine atoms with a tyrosine side chain in the protease S1 pocket are correctly predicted. Within the GABA-A/benzodiazepine ligand data set, PLS models of high predictivity resulted for our QM-based property fields, providing novel insights into key features of the SAR for two receptor subtypes and cross-receptor selectivity of the ligands. The detailed interpretation of regression models derived using improved QM-derived property fields thus provides a significant advantage by revealing chemically meaningful correlations with biological activity and helps in understanding novel structure-activity relationship features. This will allow such knowledge to be used to design novel molecules on the basis of interactions additional to steric and hydrogen-bonding features.

Development of in silico filters to predict activation of the pregnane X receptor (PXR) by structurally diverse drug-like molecules.

The pregnane X receptor (PXR), a member of the nuclear hormone superfamily, regulates the expression of several enzymes and transporters involved in metabolically relevant processes. The significant induction of CYP450 enzymes by PXR, in particular CYP3A4, might significantly alter the metabolism of prescribed drugs. In order to early identify molecules in drug discovery with a potential to activate PXR as antitarget, we developed fast and reliable in silico filters by ligand-based QSAR techniques. Two classification models were established on a diverse dataset of 434 drug-like molecules. A second augmented set allowed focusing on interesting regions in chemical space. These classifiers are based on decision trees combined with a genetic algorithm based variable selection to arrive at predictive models. The classifier for the first dataset on 29 descriptors showed good performance on a test set with a correct classification of both 100% for PXR activators and non-activators plus 87% for activators and 83% for non-activators in an external dataset. The second classifier then correctly predicts 97% activators and 91% non-activators in a test set and 94% for activators and 64% non-activators in an external set of 50 molecules, which still qualifies for application as a filter focusing on PXR activators. Finally a quantitative model for PXR activation for a subset of these molecules was derived using a regression-tree approach combined with GA variable selection. This final model shows a predictive r(2) of 0.774 for the test set and 0.452 for an external set of 33 molecules. Thus, the combination of these filters consistently provide guidelines for lowering PXR activation in novel candidate molecules.