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.

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.

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.

Automated Design of Macrocycles for Therapeutic Applications: From Small Molecules to Peptides and Proteins.

Macrocycles and cyclic peptides are increasingly attractive therapeutic modalities as they often have improved affinity, are able to bind to extended protein surfaces, and otherwise have favorable properties. Macrocyclization of a known binder may stabilize its bioactive conformation and improve its metabolic stability, cell permeability, and in certain cases oral bioavailability. Herein, we present implementation and application of an approach that automatically generates, evaluates, and proposes cyclizations utilizing a library of well-established chemical reactions and reagents. Using the three-dimensional (3D) conformation of the linear molecule in complex with a target protein as the starting point, this approach identifies attachment points, generates linkers, evaluates their geometric compatibility, and ranks the resulting molecules with respect to their predicted conformational stability and interactions with the target protein. As we show here with prospective and retrospective case studies, this procedure can be applied for the macrocyclization of small molecules and peptides and even PROteolysis TArgeting Chimeras (PROTACs) and proteins.

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.

pKa calculations for tautomerizable and conformationally flexible molecules: partition function vs. state transition approach.

Calculations of acidities of molecules with multiple tautomeric and/or conformational states require adequate treatment of the relative energetics of accessible states accompanied by a statistical-mechanical formulation of their contribution to the macroscopic pKa value. Here, we demonstrate rigorously the formal equivalence of two such approaches: a partition function treatment and statistics over transitions between molecular tautomeric and conformational states in the limit of a theory that does not require adjustment by empirical parameters correcting energetic values. However, for a frequently employed correction scheme, linear scaling of (free) energies and regression with respect to reference data taking an additive constant into account, this equivalence breaks down if more than one acid or base state is involved. The consequences of the resulting inconsistency are discussed on our datasets developed for aqueous pKa predictions during the recent SAMPL6 challenge, where molecular state energetics were computed based on the "embedded cluster reference interaction site model" (EC-RISM). This method couples integral equation theory as a solvation model to quantum-chemical calculations and yielded a test set root mean square error of 1.1 pK units from a partition function ansatz. For all practical purposes, the present results indicate that a state transition approach yields comparable accuracy despite the formal theoretical inconsistency, and that an additive regression intercept, which is strictly constant in the limit of large compound mass only, is a valid approximation. Graphical abstract Embedded cluster reference interaction site model-derived vs. experimental pKa for the test set calculated with either the partition function (blue) or the state transition approach (red), using m as a free parameter.

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.

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.

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.

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.

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 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.

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.

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.