Aqueous Drug Solubility: What Do We Measure, Calculate and QSPR Predict?

Development of an active-site titrant for SARS-CoV-2 main protease as an indispensable tool for evaluating enzyme kinetics

A titrant for the SARS-CoV-2 main protease (Mpro) was developed that enables, for the first time, the exact determination of the concentration of the enzymatically active Mpro by active-site titration. The covalent binding mode of the tetrapeptidic titrant was elucidated by the determination of the crystal structure of the enzyme-titrant complex. Four fluorogenic substrates of Mpro, including a prototypical, internally quenched Dabcyl-EDANS peptide, were compared in terms of solubility under typical assay conditions. By exploiting the new titrant, key kinetic parameters for the Mpro-catalyzed cleavage of these substrates were determined.

Targeted Protein Degradation: Advances, Challenges, and Prospects for Computational Methods

The therapeutic approach of targeted protein degradation (TPD) is gaining momentum due to its potentially superior effects compared with protein inhibition. Recent advancements in the biotech and pharmaceutical sectors have led to the development of compounds that are currently in human trials, with some showing promising clinical results. However, the use of computational tools in TPD is still limited, as it has distinct characteristics compared with traditional computational drug design methods. TPD involves creating a ternary structure (protein-degrader-ligase) responsible for the biological function, such as ubiquitination and subsequent proteasomal degradation, which depends on the spatial orientation of the protein of interest (POI) relative to E2-loaded ubiquitin. Modeling this structure necessitates a unique blend of tools initially developed for small molecules (e.g., docking) and biologics (e.g., protein-protein interaction modeling). Additionally, degrader molecules, particularly heterobifunctional degraders, are generally larger than conventional small molecule drugs, leading to challenges in determining drug-like properties like solubility and permeability. Furthermore, the catalytic nature of TPD makes occupancy-based modeling insufficient. TPD consists of multiple interconnected yet distinct steps, such as POI binding, E3 ligase binding, ternary structure interactions, ubiquitination, and degradation, along with traditional small molecule properties. A comprehensive set of tools is needed to address the dynamic nature of the induced proximity ternary complex and its implications for ubiquitination. In this Perspective, we discuss the current state of computational tools for TPD. We start by describing the series of steps involved in the degradation process and the experimental methods used to characterize them. Then, we delve into a detailed analysis of the computational tools employed in TPD. We also present an integrative approach that has proven successful for degrader design and its impact on project decisions. Finally, we examine the future prospects of computational methods in TPD and the areas with the greatest potential for impact.

Bioinformatics Tools for the Analysis of Active Compounds Identified in Ranunculaceae Species

The chemical compounds from extracts of three Ranunculaceae species, Aconitum toxicum Rchb., Anemone nemorosa L. and Helleborus odorus Waldst. & Kit. ex Willd., respectively, were isolated using the HPLC purification technique and analyzed from a bioinformatics point of view. The classes of compounds identified based on the proportion in the rhizomes/leaves/flowers used for microwave-assisted extraction and ultrasound-assisted extraction were alkaloids and phenols. Here, the quantifying of pharmacokinetics, pharmacogenomics and pharmacodynamics helps us to identify the actual biologically active compounds. Our results showed that (i) pharmacokinetically, the compounds show good absorption at the intestinal level and high permeability at the level of the central nervous system for alkaloids; (ii) regarding pharmacogenomics, alkaloids can influence tumor sensitivity and the effectiveness of some treatments; (iii) and pharmacodynamically, the compounds of these Ranunculaceae species bind to carbonic anhydrase and aldose reductase. The results obtained showed a high affinity of the compounds in the binding solution at the level of carbonic anhydrases. Carbonic anhydrase inhibitors extracted from natural sources can represent the path to new drugs useful both in the treatment of glaucoma, but also of some renal, neurological and even neoplastic diseases. The identification of natural compounds with the role of inhibitors can have a role in different types of pathologies, both associated with studied and known receptors such as carbonic anhydrase and aldose reductase, as well as new pathologies not yet addressed.

Blinded Predictions and Post Hoc Analysis of the Second Solubility Challenge Data: Exploring Training Data and Feature Set Selection for Machine and Deep Learning Models

Accurate methods to predict solubility from molecular structure are highly sought after in the chemical sciences. To assess the state of the art, the American Chemical Society organized a "Second Solubility Challenge" in 2019, in which competitors were invited to submit blinded predictions of the solubilities of 132 drug-like molecules. In the first part of this article, we describe the development of two models that were submitted to the Blind Challenge in 2019 but which have not previously been reported. These models were based on computationally inexpensive molecular descriptors and traditional machine learning algorithms and were trained on a relatively small data set of 300 molecules. In the second part of the article, to test the hypothesis that predictions would improve with more advanced algorithms and higher volumes of training data, we compare these original predictions with those made after the deadline using deep learning models trained on larger solubility data sets consisting of 2999 and 5697 molecules. The results show that there are several algorithms that are able to obtain near state-of-the-art performance on the solubility challenge data sets, with the best model, a graph convolutional neural network, resulting in an RMSE of 0.86 log units. Critical analysis of the models reveals systematic differences between the performance of models using certain feature sets and training data sets. The results suggest that careful selection of high quality training data from relevant regions of chemical space is critical for prediction accuracy but that other methodological issues remain problematic for machine learning solubility models, such as the difficulty in modeling complex chemical spaces from sparse training data sets.

Intrinsic Aqueous Solubility: Mechanistically Transparent Data-Driven Modeling of Drug Substances

Intrinsic aqueous solubility is a foundational property for understanding the chemical, technological, pharmaceutical, and environmental behavior of drug substances. Despite years of solubility research, molecular structure-based prediction of the intrinsic aqueous solubility of drug substances is still under active investigation. This paper describes the authors’ systematic data-driven modelling in which two fit-for-purpose training data sets for intrinsic aqueous solubility were collected and curated, and three quantitative structure–property relationships were derived to make predictions for the most recent solubility challenge. All three models perform well individually, while being mechanistically transparent and easy to understand. Molecular descriptors involved in the models are related to the following key steps in the solubility process: dissociation of the molecule from the crystal, formation of a cavity in the solvent, and insertion of the molecule into the solvent. A consensus modeling approach with these models remarkably improved prediction capability and reduced the number of strong outliers by more than two times. The performance and outliers of the second solubility challenge predictions were analyzed retrospectively. All developed models have been published in the QsarDB.org repository according to FAIR principles and can be used without restrictions for exploring, downloading, and making predictions.

Multi-channel GCN ensembled machine learning model for molecular aqueous solubility prediction on a clean dataset

This study constructed a new aqueous solubility dataset and a solubility regression model which was ensembled by GCN and machine learning models. Aqueous solubility is a key physiochemical property of small molecules in drug discovery. In the past few decades, there have been many studies about solubility prediction. However, many of these studies have high root mean squared error (RMSE). Meanwhile, their dataset always contains salt compounds and solubility data obtained from different experimental conditions. In this paper, we constructed a clean dataset with 2609 compounds, which was small but contains only solubility records without salts at the same temperatures (25 °C). Here, we applied graph convolutional neural network (GCN) to construct an aqueous solubility prediction model. To enhance the performance of the model, the molecular MACCS key fingerprints and physiochemical descriptors were also combined with the GCN model to build a multi-channel model. Additionally, the authors also built two machine learning models (support vector regression and gradient boost decision tree) and assembled them to the GCN model to improve the root mean squared error (RMSE = 0.665). Finally, comparative experiments have shown that our framework achieved the best performance on ESOL dataset (RMSE_val = 0.56, RMSE_test = 0.44) and surpassed four established software on aqueous solubility prediction of new compounds.

Contemporary Medicinal Chemistry Strategies for the Discovery and Development of Novel HIV-1 Non-nucleoside Reverse Transcriptase Inhibitors

Currently, HIV-1 non-nucleoside reverse transcriptase inhibitors (NNRTIs) are a major component of the highly active anti-retroviral therapy (HAART) regimen. However, the occurrence of drug-resistant strains and adverse reactions after long-term usage have inevitably compromised the clinical application of NNRTIs. Therefore, the development of novel inhibitors with distinct anti-resistance profiles and better pharmacological properties is still an enormous challenge. Herein, we summarize state-of-the-art medicinal chemistry strategies for the discovery of potent NNRTIs, such as structure-based design strategies, contemporary computer-aided drug design, covalent-binding strategies, and the application of multi-target-directed ligands. The strategies described here will facilitate the identification of promising HIV-1 NNRTIs.

Pushing the limits of solubility prediction via quality-oriented data selection

Accurate prediction of the solubility of chemical substances in solvents remains a challenge. The sparsity of high-quality solubility data is recognized as the biggest hurdle in the development of robust data-driven methods for practical use. Nonetheless, the effects of the quality and quantity of data on aqueous solubility predictions have not yet been scrutinized. In this study, the roles of the size and the quality of data sets on the performances of the solubility prediction models are unraveled, and the concepts of actual and observed performances are introduced. In an effort to curtail the gap between actual and observed performances, a quality-oriented data selection method, which evaluates the quality of data and extracts the most accurate part of it through statistical validation, is designed. Applying this method on the largest publicly available solubility database and using a consensus machine learning approach, a top-performing solubility prediction model is achieved.

ADME prediction with KNIME: In silico aqueous solubility consensus model based on supervised recursive random forest approaches

In-silico prediction of aqueous solubility plays an important role during the drug discovery and development processes. For many years, the limited performance of in-silico solubility models has been attributed to the lack of high-quality solubility data for pharmaceutical molecules. However, some studies suggest that the poor accuracy of solubility prediction is not related to the quality of the experimental data and that more precise methodologies (algorithms and/or set of descriptors) are required for predicting aqueous solubility for pharmaceutical molecules. In this study a large and diverse database was generated with aqueous solubility values collected from two public sources; two new recursive machine-learning approaches were developed for data cleaning and variable selection, and a consensus model based on regression and classification algorithms was created. The modeling protocol, which includes the curation of chemical and experimental data, was implemented in KNIME, with the aim of obtaining an automated workflow for the prediction of new databases. Finally, we compared several methods or models available in the literature with our consensus model, showing results comparable or even outperforming previous published models.

Fundamental aspects of DMPK optimization of targeted protein degraders

Targeted protein degraders are an emerging modality. Their properties fall outside the traditional small-molecule property space and are in the 'beyond rule of 5' space. Consequently, drug discovery programs focused on developing orally bioavailable degraders are expected to face complex drug metabolism and pharmacokinetics (DMPK) challenges compared with traditional small molecules. Nevertheless, little information is available on the DMPK optimization of oral degraders. Therefore, in this review, we discuss our experience of these DMPK optimization challenges and present methodologies and strategies to overcome the hurdles dealing with this new small-molecule modality, specifically in developing oral degraders to treat cancer.