Machine Learning for Fast, Quantum Mechanics-Based Approximation of Drug Lipophilicity

Accurately Predicting Barrier Heights for Radical Reactions in Solution Using Deep Graph Networks

Quantitative estimates of reaction barriers and solvent effects are essential for developing kinetic mechanisms and predicting reaction outcomes. Here, we create a new data set of 5,600 unique elementary radical reactions calculated using the M06-2X/def2-QZVP//B3LYP-D3(BJ)/def2-TZVP level of theory. A conformer search is done for each species using TPSS/def2-TZVP. Gibbs free energies of activation and of reaction for these radical reactions in 40 common solvents are obtained using COSMO-RS for solvation effects. These balanced reactions involve the elements H, C, N, O, and S, contain up to 19 heavy atoms, and have atom-mapped SMILES. All transition states are verified by an intrinsic reaction coordinate calculation. We next train a deep graph network to directly estimate the Gibbs free energy of activation and of reaction in both gas and solution phases using only the atom-mapped SMILES of the reactant and product and the SMILES of the solvent. This simple input representation avoids computationally expensive optimizations for the reactant, transition state, and product structures during inference, making our model well-suited for high-throughput predictive chemistry and quickly providing information for (retro-)synthesis planning tools. To properly measure model performance, we report results on both interpolative and extrapolative data splits and also compare to several baseline models. During training and testing, the data set is augmented by including the reverse direction of each reaction and variants with different resonance structures. After data augmentation, we have around 2 million entries to train the model, which achieves a testing set mean absolute error of 1.16 kcal mol-1 for the Gibbs free energy of activation in solution. We anticipate this model will accelerate predictions for high-throughput screening to quickly identify relevant reactions in solution, and our data set will serve as a benchmark for future studies.

Geometric deep learning for molecular property predictions with chemical accuracy across chemical space

Chemical engineers heavily rely on precise knowledge of physicochemical properties to model chemical processes. Despite the growing popularity of deep learning, it is only rarely applied for property prediction due to data scarcity and limited accuracy for compounds in industrially-relevant areas of the chemical space. Herein, we present a geometric deep learning framework for predicting gas- and liquid-phase properties based on novel quantum chemical datasets comprising 124,000 molecules. Our findings reveal that the necessity for quantum-chemical information in deep learning models varies significantly depending on the modeled physicochemical property. Specifically, our top-performing geometric model meets the most stringent criteria for "chemically accurate" thermochemistry predictions. We also show that by carefully selecting the appropriate model featurization and evaluating prediction uncertainties, the reliability of the predictions can be strongly enhanced. These insights represent a crucial step towards establishing deep learning as the standard property prediction workflow in both industry and academia.Scientific contributionWe propose a flexible property prediction tool that can handle two-dimensional and three-dimensional molecular information. A thermochemistry prediction methodology that achieves high-level quantum chemistry accuracy for a broad application range is presented. Trained deep learning models and large novel molecular databases of real-world molecules are provided to offer a directly usable and fast property prediction solution to practitioners.

Towards development of new antimalarial compounds through in silico and in vitro assays

Malaria is one of most widespread infectious disease in world. The antimalarial therapy presents a series of limitations, such as toxicity and the emergence of resistance, which makes the search for new drugs urgent. Thus, it becomes necessary to explore essential and exclusive therapeutic targets of the parasite to achieve selective inhibition. Enoyl-ACP reductase is an enzyme of the type II fatty acid biosynthetic pathway and is responsible for the rate-limiting step in the fatty acid elongation cycle. In this work, we use hierarchical virtual screening and drug repositioning strategies to prioritize compounds for phenotypic assays and molecular dynamics studies. The molecules were tested against chloroquine-resistant W2 strain of Plasmodium falciparum (EC50 between 330.05 and 13.92 µM). Nitrofurantoin was the best antimalarial activity at low micromolar range (EC50 = 13.92 µM). However, a hit compound against malaria must have a biological activity value below 1 µM. A large number of molecules present problems with permeability in biological membranes and reaching an effective concentration in their target's microenvironment. Nitrofurantoin derivatives with inclusions of groups which confer increased lipid solubility (methyl groups, halogens and substituted and unsubstituted aromatic rings) have been proposed. These derivatives were pulled through the lipid bilayer in molecular dynamics simulations. Molecules 14, 18 and 21 presented lower free energy values than nitrofurantoin when crossing the lipid bilayer.

Quantum mechanical-based strategies in drug discovery: Finding the pace to new challenges in drug design

The expansion of the chemical space to tangible libraries containing billions of synthesizable molecules opens exciting opportunities for drug discovery, but also challenges the power of computer-aided drug design to prioritize the best candidates. This directly hits quantum mechanics (QM) methods, which provide chemically accurate properties, but subject to small-sized systems. Preserving accuracy while optimizing the computational cost is at the heart of many efforts to develop high-quality, efficient QM-based strategies, reflected in refined algorithms and computational approaches. The design of QM-tailored physics-based force fields and the coupling of QM with machine learning, in conjunction with the computing performance of supercomputing resources, will enhance the ability to use these methods in drug discovery. The challenge is formidable, but we will undoubtedly see impressive advances that will define a new era.

PICKAPEP: An application for parameter calculation and visualization of cyclized and modified peptidomimetics

The interest in peptides and especially in peptidomimetic structures has risen enormously in the past few years. Novel modification strategies including nonnatural amino acids, sophisticated cyclization strategies, and side chain modifications to improve the pharmacokinetic properties of peptides are continuously arising. However, a calculator tool accompanying the current development in peptide sciences towards modified peptides is missing. Herein, we present the application PICKAPEP, enabling the virtual construction and visualization of peptidomimetics ranging from well-known cyclized and modified peptides such as ciclosporin A up to fully self-designed peptide-based structures with custom amino acids. Calculated parameters include the molecular weight, the water-octanol partition coefficient, the topological polar surface area, the number of rotatable bonds, and the peptide SMILES code. To our knowledge, PICKAPEP is the first tool allowing users to add custom amino acids as building blocks and also the only tool giving the possibility to process large peptide libraries and calculate parameters for multiple peptides at once. We believe that PICKAPEP will support peptide researchers in their work and will find wide application in current as well as future peptide drug development processes. PICKAPEP is available open source for Windows and Mac operating systems (https://urldefense.com/v3/__https://www.research-collection.ethz.ch/handle/20.500.11850/681174__;!!N11eV2iwtfs!qt5f_2lNd6IZUDH1HVSVwg0zYzS8-nFazQ8c61jS5GaD5vkVS5C3igyfh3haJRnaX8ugW7o9VWUiCihPqcptmaWoqwYf9LvZTQ$).

In silico ADME/tox comes of age: twenty years later

In the early 2000s pharmaceutical drug discovery was beginning to use computational approaches for absorption, distribution, metabolism, excretion and toxicity (ADME/Tox, also known as ADMET) prediction. This emphasis on prediction was an effort to reduce the risk of later stage failures from ADME/Tox.Much has been written in the intervening twenty plus years and significant expenditure has occurred in companies developing these in silico capabilities which can be gleaned from publications. It is therefore an appropriate time to briefly reflect on what was proposed then and what the reality is today.20 years ago, we tended to optimise bioactivity and perhaps one ADME/Tox property at a time. Previously pharmaceutical companies needed a whole infrastructure for models - in silico and in vitro experts, IT, champions on a project team, educators and management support. Now we are in the age of generative de novo design where bioactivity and many ADME/Tox properties can be optimised and large language model technologies are available.There are also some challenges such as the focus on very large molecules which may be outside of current ADME/Tox models.We provide an opportunity to look forward with the increasing public data for ADME/Tox as well as expanded types of algorithms available.

Topological regression as an interpretable and efficient tool for quantitative structure-activity relationship modeling

Quantitative structure-activity relationship (QSAR) modeling is a powerful tool for drug discovery, yet the lack of interpretability of commonly used QSAR models hinders their application in molecular design. We propose a similarity-based regression framework, topological regression (TR), that offers a statistically grounded, computationally fast, and interpretable technique to predict drug responses. We compare the predictive performance of TR on 530 ChEMBL human target activity datasets against the predictive performance of deep-learning-based QSAR models. Our results suggest that our sparse TR model can achieve equal, if not better, performance than the deep learning-based QSAR models and provide better intuitive interpretation by extracting an approximate isometry between the chemical space of the drugs and their activity space.

Development of a Novel In Silico Classification Model to Assess Reactive Metabolite Formation in the Cysteine Trapping Assay and Investigation of Important Substructures

Predicting whether a compound can cause drug-induced liver injury (DILI) is difficult due to the complexity of drug mechanism. The cysteine trapping assay is a method for detecting reactive metabolites that bind to microsomes covalently. However, it is cumbersome to use 35S isotope-labeled cysteine for this assay. Therefore, we constructed an in silico classification model for predicting a positive/negative outcome in the cysteine trapping assay. We collected 475 compounds (436 in-house compounds and 39 publicly available drugs) based on experimental data performed in this study, and the composition of the results showed 248 positives and 227 negatives. Using a Message Passing Neural Network (MPNN) and Random Forest (RF) with extended connectivity fingerprint (ECFP) 4, we built machine learning models to predict the covalent binding risk of compounds. In the time-split dataset, AUC-ROC of MPNN and RF were 0.625 and 0.559 in the hold-out test, restrictively. This result suggests that the MPNN model has a higher predictivity than RF in the time-split dataset. Hence, we conclude that the in silico MPNN classification model for the cysteine trapping assay has a better predictive power. Furthermore, most of the substructures that contributed positively to the cysteine trapping assay were consistent with previous results.

Prospective de novo drug design with deep interactome learning

De novo drug design aims to generate molecules from scratch that possess specific chemical and pharmacological properties. We present a computational approach utilizing interactome-based deep learning for ligand- and structure-based generation of drug-like molecules. This method capitalizes on the unique strengths of both graph neural networks and chemical language models, offering an alternative to the need for application-specific reinforcement, transfer, or few-shot learning. It enables the "zero-shot" construction of compound libraries tailored to possess specific bioactivity, synthesizability, and structural novelty. In order to proactively evaluate the deep interactome learning framework for protein structure-based drug design, potential new ligands targeting the binding site of the human peroxisome proliferator-activated receptor (PPAR) subtype gamma are generated. The top-ranking designs are chemically synthesized and computationally, biophysically, and biochemically characterized. Potent PPAR partial agonists are identified, demonstrating favorable activity and the desired selectivity profiles for both nuclear receptors and off-target interactions. Crystal structure determination of the ligand-receptor complex confirms the anticipated binding mode. This successful outcome positively advocates interactome-based de novo design for application in bioorganic and medicinal chemistry, enabling the creation of innovative bioactive molecules.

Design Principles for Balancing Lipophilicity and Permeability in beyond Rule of 5 Space

An ab initio conformational analysis of oral beyond Rule of 5 (bRo5) drugs was complemented with measured permeability and logP(octanol) to derive design principles conferring oral bioavailability. 3D polar surface area (PSA) thresholds for oral bRo5 drugs coincided with those reported for Ro5 space. The majority of oral bRo5 drugs exceeded the Ro5 logP threshold of 5, reflecting a bias for permeability. Above 500 Da molecular weight (MW), oral drugs and highly permeable Novartis compounds occupy a narrow polarity range (topological or TPSA/MW) of 0.1-0.3 Å2 /Da, whose upper half coincides with the lower 90 percentiles of the Novartis logP set. This TPSA/MW range and 3D PSA below 100 Å2 define the "Rule of ~1 /₅" for balancing lipophilicity and permeability. Neutral TPSA, defined as TPSA minus 3D PSA occurs independent of conformation, intramolecular hydrogen bonds (IMHB) and MW, suggesting it is an intrinsic molecular property. Neutral TPSA increased in the lead optimization (LO) campaigns of three first in class de novo designed bRo5 drugs and may be a useful design parameter in bRo5 space.

Enabling late-stage drug diversification by high-throughput experimentation with geometric deep learning

Late-stage functionalization is an economical approach to optimize the properties of drug candidates. However, the chemical complexity of drug molecules often makes late-stage diversification challenging. To address this problem, a late-stage functionalization platform based on geometric deep learning and high-throughput reaction screening was developed. Considering borylation as a critical step in late-stage functionalization, the computational model predicted reaction yields for diverse reaction conditions with a mean absolute error margin of 4-5%, while the reactivity of novel reactions with known and unknown substrates was classified with a balanced accuracy of 92% and 67%, respectively. The regioselectivity of the major products was accurately captured with a classifier F-score of 67%. When applied to 23 diverse commercial drug molecules, the platform successfully identified numerous opportunities for structural diversification. The influence of steric and electronic information on model performance was quantified, and a comprehensive simple user-friendly reaction format was introduced that proved to be a key enabler for seamlessly integrating deep learning and high-throughput experimentation for late-stage functionalization.

Exploring protein-ligand binding affinity prediction with electron density-based geometric deep learning

Rational structure-based drug design relies on accurate predictions of protein-ligand binding affinity from structural molecular information. Although deep learning-based methods for predicting binding affinity have shown promise in computational drug design, certain approaches have faced criticism for their potential to inadequately capture the fundamental physical interactions between ligands and their macromolecular targets or for being susceptible to dataset biases. Herein, we propose to include bond-critical points based on the electron density of a protein-ligand complex as a fundamental physical representation of protein-ligand interactions. Employing a geometric deep learning model, we explore the usefulness of these bond-critical points to predict absolute binding affinities of protein-ligand complexes, benchmark model performance against existing methods, and provide a critical analysis of this new approach. The models achieved root-mean-squared errors of 1.4-1.8 log units on the PDBbind dataset, and 1.0-1.7 log units on the PDE10A dataset, not indicating significant advantages over benchmark methods, and thus rendering the utility of electron density for deep learning models context-dependent. The relationship between intermolecular electron density and corresponding binding affinity was analyzed, and Pearson correlation coefficients r > 0.7 were obtained for several macromolecular targets.

Chemprop: A Machine Learning Package for Chemical Property Prediction

Deep learning has become a powerful and frequently employed tool for the prediction of molecular properties, thus creating a need for open-source and versatile software solutions that can be operated by nonexperts. Among the current approaches, directed message-passing neural networks (D-MPNNs) have proven to perform well on a variety of property prediction tasks. The software package Chemprop implements the D-MPNN architecture and offers simple, easy, and fast access to machine-learned molecular properties. Compared to its initial version, we present a multitude of new Chemprop functionalities such as the support of multimolecule properties, reactions, atom/bond-level properties, and spectra. Further, we incorporate various uncertainty quantification and calibration methods along with related metrics as well as pretraining and transfer learning workflows, improved hyperparameter optimization, and other customization options concerning loss functions or atom/bond features. We benchmark D-MPNN models trained using Chemprop with the new reaction, atom-level, and spectra functionality on a variety of property prediction data sets, including MoleculeNet and SAMPL, and observe state-of-the-art performance on the prediction of water-octanol partition coefficients, reaction barrier heights, atomic partial charges, and absorption spectra. Chemprop enables out-of-the-box training of D-MPNN models for a variety of problem settings in fast, user-friendly, and open-source software.

Identifying opportunities for late-stage C-H alkylation with high-throughput experimentation and in silico reaction screening

Enhancing the properties of advanced drug candidates is aided by the direct incorporation of specific chemical groups, avoiding the need to construct the entire compound from the ground up. Nevertheless, their chemical intricacy often poses challenges in predicting reactivity for C-H activation reactions and planning their synthesis. We adopted a reaction screening approach that combines high-throughput experimentation (HTE) at a nanomolar scale with computational graph neural networks (GNNs). This approach aims to identify suitable substrates for late-stage C-H alkylation using Minisci-type chemistry. GNNs were trained using experimentally generated reactions derived from in-house HTE and literature data. These trained models were then used to predict, in a forward-looking manner, the coupling of 3180 advanced heterocyclic building blocks with a diverse set of sp3-rich carboxylic acids. This predictive approach aimed to explore the substrate landscape for Minisci-type alkylations. Promising candidates were chosen, their production was scaled up, and they were subsequently isolated and characterized. This process led to the creation of 30 novel, functionally modified molecules that hold potential for further refinement. These results positively advocate the application of HTE-based machine learning to virtual reaction screening.

The rise of automated curiosity-driven discoveries in chemistry

The quest for generating novel chemistry knowledge is critical in scientific advancement, and machine learning (ML) has emerged as an asset in this pursuit. Through interpolation among learned patterns, ML can tackle tasks that were previously deemed demanding to machines. This distinctive capacity of ML provides invaluable aid to bench chemists in their daily work. However, current ML tools are typically designed to prioritize experiments with the highest likelihood of success, i.e., higher predictive confidence. In this perspective, we build on current trends that suggest a future in which ML could be just as beneficial in exploring uncharted search spaces through simulated curiosity. We discuss how low and 'negative' data can catalyse one-/few-shot learning, and how the broader use of curious ML and novelty detection algorithms can propel the next wave of chemical discoveries. We anticipate that ML for curiosity-driven research will help the community overcome potentially biased assumptions and uncover unexpected findings in the chemical sciences at an accelerated pace.

Making the Case for Quantum Mechanics in Predictive Toxicology─Nearly 100 Years Too Late?

The use of quantum mechanics (QM) has long been the norm to study covalent-binding phenomena in chemistry and biochemistry. The pharmaceutical industry leverages QM models explicitly in covalent drug discovery and implicitly to characterize short-range interactions in noncovalent binding. Predictive toxicology has resisted widespread adoption of QM, including in the pharmaceutical industry, despite its obvious relevance to the metabolic processes in the upstream of adverse outcome pathways and advances in both QM methods and computational resources, which support fit-for-purpose applications in reasonable timeframes. Here, we make the case for embracing QM as an indispensable part of a toxicologist's toolkit. We argue that QM provides the necessary orthogonality to alert-based expert systems and traditional QSARs, consistent with calls for animal-free integrated testing strategies for safety assessments of commercial chemicals. We outline existing roadblocks to this transition, including the need to train model developers in QM and the shift toward service-based toxicity models that utilize high-performance computing clusters. Lastly, we describe recent examples of successful implementations of QM in hazard assessments and propose how in silico toxicology can be further advanced by integrating QM with artificial intelligence.

Open-Source Machine Learning in Computational Chemistry

The field of computational chemistry has seen a significant increase in the integration of machine learning concepts and algorithms. In this Perspective, we surveyed 179 open-source software projects, with corresponding peer-reviewed papers published within the last 5 years, to better understand the topics within the field being investigated by machine learning approaches. For each project, we provide a short description, the link to the code, the accompanying license type, and whether the training data and resulting models are made publicly available. Based on those deposited in GitHub repositories, the most popular employed Python libraries are identified. We hope that this survey will serve as a resource to learn about machine learning or specific architectures thereof by identifying accessible codes with accompanying papers on a topic basis. To this end, we also include computational chemistry open-source software for generating training data and fundamental Python libraries for machine learning. Based on our observations and considering the three pillars of collaborative machine learning work, open data, open source (code), and open models, we provide some suggestions to the community.

Structure-based drug design with geometric deep learning

Structure-based drug design uses three-dimensional geometric information of macromolecules, such as proteins or nucleic acids, to identify suitable ligands. Geometric deep learning, an emerging concept of neural-network-based machine learning, has been applied to macromolecular structures. This review provides an overview of the recent applications of geometric deep learning in bioorganic and medicinal chemistry, highlighting its potential for structure-based drug discovery and design. Emphasis is placed on molecular property prediction, ligand binding site and pose prediction, and structure-based de novo molecular design. The current challenges and opportunities are highlighted, and a forecast of the future of geometric deep learning for drug discovery is presented.