SVM approach for predicting LogP

Collection of Partition Coefficients in Hexadecyltrimethylammonium Bromide, Sodium Cholate, and Lithium Perfluorooctanesulfonate Micellar Solutions: Experimental Determination and Computational Predictions

This study focuses on determining the partition coefficients (logP) of a diverse set of 63 molecules in three distinct micellar systems: hexadecyltrimethylammonium bromide (HTAB), sodium cholate (SC), and lithium perfluorooctanesulfonate (LPFOS). The experimental log p values were obtained through micellar electrokinetic chromatography (MEKC) experiments, conducted under controlled pH conditions. Then, Quantum Mechanics (QM) and machine learning approaches are proposed for the prediction of the partition coefficients in these three micellar systems. In the applied QM approach, the experimentally obtained partition coefficients were correlated with the calculated values for the case of the 15 solvent mixtures. Using Density Function Theory (DFT) with the B3LYP functional, we calculated the solvation free energies of 63 molecules in these 16 solvents. The combined data from the experimental partition coefficients in the three micellar formulations showed that the 1-propanol/water combination demonstrated the best agreement with the experimental partition coefficients for the SC and HTAB micelles. Moreover, we employed the SVM approach and k-means clustering based on the generation of the chemical descriptor space. The analysis revealed distinct partitioning patterns associated with specific characteristic features within each identified class. These results indicate the utility of the combined techniques when we want an efficient and quicker model for predicting partition coefficients in diverse micelles.

Development and test of highly accurate endpoint free energy methods. 2: Prediction of logarithm of n-octanol-water partition coefficient (logP) for druglike molecules using MM-PBSA method

The logarithm of n-octanol-water partition coefficient (logP) is frequently used as an indicator of lipophilicity in drug discovery, which has substantial impacts on the absorption, distribution, metabolism, excretion, and toxicity of a drug candidate. Considering that the experimental measurement of the property is costly and time-consuming, it is of great importance to develop reliable prediction models for logP. In this study, we developed a transfer free energy-based logP prediction model-FElogP. FElogP is based on the simple principle that logP is determined by the free energy change of transferring a molecule from water to n-octanol. The underlying physical method to calculate transfer free energy is the molecular mechanics-Poisson Boltzmann surface area (MM-PBSA), thus this method is named as free energy-based logP (FElogP). The superiority of FElogP model was validated by a large set of 707 structurally diverse molecules in the ZINC database for which the measurement was of high quality. Encouragingly, FElogP outperformed several commonly-used QSPR or machine learning-based logP models, as well as some continuum solvation model-based methods. The root-mean-square error (RMSE) and Pearson correlation coefficient (R) between the predicted and measured values are 0.91 log units and 0.71, respectively, while the runner-up, the logP model implemented in OpenBabel had an RMSE of 1.13 log units and R of 0.67. Given the fact that FElogP was not parameterized against experimental logP directly, its excellent performance is likely to be expanded to arbitrary organic molecules covered by the general AMBER force fields.

A New Formulation for the Concerted Alchemical Calculation of van der Waals and Coulomb Components of Solvation Free Energies

Alchemical free energy calculations via molecular dynamics have been widely used to obtain thermodynamic properties related to protein-ligand binding and solute-solvent interactions. Although soft-core modeling is the most common approach, the linear basis function (LBF) methodology [Naden, L. N.; et al. J. Chem. Theory Comput.2014, 10 (3), 1128; 2015, 11 (6), 2536] has emerged as a suitable alternative. It overcomes the end-point singularity of the scaling method while maintaining essential advantages such as ease of implementation and high flexibility for postprocessing analysis. In the present work, we propose a simple LBF variant and formulate an efficient protocol for evaluating van der Waals and Coulomb components of an alchemical transformation in tandem, in contrast to the prevalent sequential evaluation mode. To validate our proposal, which results from a careful optimization study, we performed solvation free energy calculations and obtained octanol-water partition coefficients of small organic molecules. Comparisons with results obtained via the sequential mode using either another LBF approach or the soft-core model attest to the effectiveness and correctness of our method. In addition, we show that a reaction field model with an infinite dielectric constant can provide very accurate hydration free energies when used instead of a lattice-sum method to model solute-solvent electrostatics.

Grid inhomogeneous solvation theory for cross-solvation in rigid solvents

Grid Inhomogeneous Solvation Theory (GIST) has proven useful to calculate localized thermodynamic properties of water around a solute. Numerous studies have leveraged this information to enhance structure-based binding predictions. We have recently extended GIST toward chloroform as a solvent to allow the prediction of passive membrane permeability. Here, we further generalize the GIST algorithm toward all solvents that can be modeled as rigid molecules. This restriction is inherent to the method and is already present in the inhomogeneous solvation theory. Here, we show that our approach can be applied to various solvent molecules by comparing the results of GIST simulations with thermodynamic integration (TI) calculations and experimental results. Additionally, we analyze and compare a matrix consisting of 100 entries of ten different solvent molecules solvated within each other. We find that the GIST results are highly correlated with TI calculations as well as experiments. For some solvents, we find Pearson correlations of up to 0.99 to the true entropy, while others are affected by the first-order approximation more strongly. The enthalpy-entropy splitting provided by GIST allows us to extend a recently published approach, which estimates higher order entropies by a linear scaling of the first-order entropy, to solvents other than water. Furthermore, we investigate the convergence of GIST in different solvents. We conclude that our extension to GIST reliably calculates localized thermodynamic properties for different solvents and thereby significantly extends the applicability of this widely used method.

Energy-entropy prediction of octanol-water logP of SAMPL7 N-acyl sulfonamide bioisosters

Partition coefficients quantify a molecule's distribution between two immiscible liquid phases. While there are many methods to compute them, there is not yet a method based on the free energy of each system in terms of energy and entropy, where entropy depends on the probability distribution of all quantum states of the system. Here we test a method in this class called Energy Entropy Multiscale Cell Correlation (EE-MCC) for the calculation of octanol-water logP values for 22 N-acyl sulfonamides in the SAMPL7 Physical Properties Challenge (Statistical Assessment of the Modelling of Proteins and Ligands). EE-MCC logP values have a mean error of 1.8 logP units versus experiment and a standard error of the mean of 1.0 logP units for three separate calculations. These errors are primarily due to getting sufficiently converged energies to give accurate differences of large numbers, particularly for the large-molecule solvent octanol. However, this is also an issue for entropy, and approximations in the force field and MCC theory also contribute to the error. Unique to MCC is that it explains the entropy contributions over all the degrees of freedom of all molecules in the system. A gain in orientational entropy of water is the main favourable entropic contribution, supported by small gains in solute vibrational and orientational entropy but offset by unfavourable changes in the orientational entropy of octanol, the vibrational entropy of both solvents, and the positional and conformational entropy of the solute.

ClassicalGSG: Prediction of log P using classical molecular force fields and geometric scattering for graphs

This work examines methods for predicting the partition coefficient (log P) for a dataset of small molecules. Here, we use atomic attributes such as radius and partial charge, which are typically used as force field parameters in classical molecular dynamics simulations. These atomic attributes are transformed into index-invariant molecular features using a recently developed method called geometric scattering for graphs (GSG). We call this approach "ClassicalGSG" and examine its performance under a broad range of conditions and hyperparameters. We train ClassicalGSG log P predictors with neural networks using 10,722 molecules from the OpenChem dataset and apply them to predict the log P values from four independent test sets. The ClassicalGSG method's performance is compared to a baseline model that employs graph convolutional networks. Our results show that the best prediction accuracies are obtained using atomic attributes generated with the CHARMM generalized force field and 2D molecular structures.

Artificial intelligence and machine learning-aided drug discovery in central nervous system diseases: State-of-the-arts and future directions

Neurological disorders significantly outnumber diseases in other therapeutic areas. However, developing drugs for central nervous system (CNS) disorders remains the most challenging area in drug discovery, accompanied with the long timelines and high attrition rates. With the rapid growth of biomedical data enabled by advanced experimental technologies, artificial intelligence (AI) and machine learning (ML) have emerged as an indispensable tool to draw meaningful insights and improve decision making in drug discovery. Thanks to the advancements in AI and ML algorithms, now the AI/ML-driven solutions have an unprecedented potential to accelerate the process of CNS drug discovery with better success rate. In this review, we comprehensively summarize AI/ML-powered pharmaceutical discovery efforts and their implementations in the CNS area. After introducing the AI/ML models as well as the conceptualization and data preparation, we outline the applications of AI/ML technologies to several key procedures in drug discovery, including target identification, compound screening, hit/lead generation and optimization, drug response and synergy prediction, de novo drug design, and drug repurposing. We review the current state-of-the-art of AI/ML-guided CNS drug discovery, focusing on blood-brain barrier permeability prediction and implementation into therapeutic discovery for neurological diseases. Finally, we discuss the major challenges and limitations of current approaches and possible future directions that may provide resolutions to these difficulties.

Bioconcentration potential and microbial toxicity of onium cations in photoacid generators

Despite the widespread utilization of onium salts as photoacid generators (PAGs) in semiconductor photolithography, their environmental, health, and safety (EHS) properties remain poorly understood. The present work reports the bioconcentration potential of five representative onium species (four sulfonium and one iodonium compound) by determining the octanol-water partition coefficient (P_OW) and lipid membrane affinity coefficient (K_MA); microbial toxicity was evaluated using the bioluminescent bacterium Aliivibrio fischeri (Microtox bioassay). Four of the oniums exhibited varying degrees of hydrophobic (lipophilic) partitioning (log P_OW: 0.08-4.12; K_MA: 1.70-5.62). A strong positive linear correlation was observed between log P_OW and K_MA (K_MA = log P_OW + 1.76, R² = 0.99). The bioconcentration factors (log BCF) estimated from P_OW and K_MA for the four oniums ranged from 0.13 to 3.67 L kg^-1. Bis-(4-tert-butyl phenyl)-iodonium and triphenylsulfonium had 50% inhibitory concentrations (IC₅₀) of 4.8 and 84.6 μM, whereas the IC₅₀ values of the other three oniums were not determined because these values were higher than their aqueous solubility. Given the increased regulatory scrutiny regarding the use and potential health impacts from onium PAGs, this study fulfills critical knowledge gaps concerning the EHS properties of PAG oniums, enabling more comprehensive evaluation of their environmental impacts and potential risk management strategies.

Solvation Thermodynamics in Different Solvents: Water-Chloroform Partition Coefficients from Grid Inhomogeneous Solvation Theory

Reliable information on partition coefficients plays a key role in drug development, as solubility decisively affects bioavailability. In a physicochemical context, the partition coefficient of a solute between two different solvents can be described as a function of solvation free energies. Hence, substantial scientific efforts have been made toward accurate predictions of solvation free energies in various solvents. The grid inhomogeneous solvation theory (GIST) facilitates the calculation of solvation free energies. In this study, we introduce an extended version of the GIST algorithm, which enables the calculation for chloroform in addition to water. Furthermore, GIST allows localization of enthalpic and entropic contributions. We test our approach by calculating partition coefficients between water and chloroform for a set of eight small molecules. We report a Pearson correlation coefficient of 0.96 between experimentally determined and calculated partition coefficients. The capability to reliably predict partition coefficients between water and chloroform and the possibility to localize their contributions allow the optimization of a compound's partition coefficient. Therefore, we presume that this methodology will be of great benefit for the efficient development of pharmaceuticals.

Concepts of Artificial Intelligence for Computer-Assisted Drug Discovery

Artificial intelligence (AI), and, in particular, deep learning as a subcategory of AI, provides opportunities for the discovery and development of innovative drugs. Various machine learning approaches have recently (re)emerged, some of which may be considered instances of domain-specific AI which have been successfully employed for drug discovery and design. This review provides a comprehensive portrayal of these machine learning techniques and of their applications in medicinal chemistry. After introducing the basic principles, alongside some application notes, of the various machine learning algorithms, the current state-of-the art of AI-assisted pharmaceutical discovery is discussed, including applications in structure- and ligand-based virtual screening, de novo drug design, physicochemical and pharmacokinetic property prediction, drug repurposing, and related aspects. Finally, several challenges and limitations of the current methods are summarized, with a view to potential future directions for AI-assisted drug discovery and design.

Large-scale ligand-based predictive modelling using support vector machines

The increasing size of datasets in drug discovery makes it challenging to build robust and accurate predictive models within a reasonable amount of time. In order to investigate the effect of dataset sizes on predictive performance and modelling time, ligand-based regression models were trained on open datasets of varying sizes of up to 1.2 million chemical structures. For modelling, two implementations of support vector machines (SVM) were used. Chemical structures were described by the signatures molecular descriptor. Results showed that for the larger datasets, the LIBLINEAR SVM implementation performed on par with the well-established libsvm with a radial basis function kernel, but with dramatically less time for model building even on modest computer resources. Using a non-linear kernel proved to be infeasible for large data sizes, even with substantial computational resources on a computer cluster. To deploy the resulting models, we extended the Bioclipse decision support framework to support models from LIBLINEAR and made our models of logD and solubility available from within Bioclipse.

Calculating Partition Coefficients of Small Molecules in Octanol/Water and Cyclohexane/Water

Partition coefficients describe how a solute is distributed between two immiscible solvents. They are used in drug design as a measure of a solute's hydrophobicity and a proxy for its membrane permeability. We calculate partition coefficients from transfer free energies using molecular dynamics simulations in explicit solvent. Setup is done by our new Solvation Toolkit which automates the process of creating input files for any combination of solutes and solvents for many popular molecular dynamics software packages. We calculate partition coefficients between octanol/water and cyclohexane/water with the Generalized AMBER Force Field (GAFF) and the Dielectric Corrected GAFF (GAFF-DC). With similar methods in the past we found a root-mean-squared error (RMSE) of 6.3 kJ/mol in hydration free energies which would correspond to an error of around 1.6 log units in partition coefficients if solvation free energies in both solvents were estimated with comparable accuracy. Here we find an overall RMSE of about 1.2 log units with both force fields. Results from GAFF and GAFF-DC seem to exhibit systematic biases in opposite directions for calculated cyclohexane/water partition coefficients.

Research on controllable degradation of sulfonylurea herbicides

Through studying structure, bioassay and soil degradation tri-factor relationship, potential controllable degradation of SU was firstly explored and summarized.

Synthesis and evaluation of novel prenylated chalcone derivatives as anti-leishmanial and anti-trypanosomal compounds

Chalcones form a class of compounds that belong to the flavonoid family and are widely distributed in plants. Their simple structure and the ease of preparation make chalcones attractive scaffolds for the synthesis of a large number of derivatives enabling the evaluation of the effects of different functional groups on biological activities. In this Letter, we report the successful synthesis of a series of novel prenylated chalcones via Claisen-Schmidt condensation and the evaluation of their effect on the viability of the Trypanosomatidae parasites Leishmania amazonensis, Leishmania infantum and Trypanosoma cruzi.

Prediction and interpretation of the lipophilicity of small peptides

Peptide-based drug discovery has considerably expanded and solid in silico tools for the prediction of physico-chemical properties of peptides are urgently needed. In this work we tested some combinations of descriptors/algorithms to find the best model to predict [Formula: see text] of a series of peptides. To do that we evaluate the models statistical performances but also their skills in providing a reliable deconvolution of the balance of intermolecular forces governing the partitioning phenomenon. Results prove that a PLS model based on VolSurf+ descriptors is the best tool to predict [Formula: see text] of neutral and ionised peptides. The mechanistic interpretation also reveals that the inclusion in the chemical structure of a HBD group is more efficient in decreasing lipophilicity than the inclusion of a HBA group.

Design, synthesis, and biological evaluation of novel trifluoromethyl indoles as potent HIV-1 NNRTIs with an improved drug resistance profile

A novel series of trifluoromethyl indole derivatives have been designed, synthesized and evaluated for anti-HIV-1 activities in MT-2 cells. The hydrophobic constant, acute toxicity, carcinogenicity and mutagenicity were predicted. Trifluoromethyl indoles 10i and 10k showed extremely promising activities against WT HIV-1 with IC50 values at the low nanomolar level, similar to efavirenz, better than nevirapine, and also possessed higher potency towards the drug-resistant mutant strain Y181C than nevirapine. Preliminary SAR and docking studies of detailed binding mode provided some insights for discovery of more potent NNRTIs.

Binary classification of aqueous solubility using support vector machines with reduction and recombination feature selection

Aqueous solubility is recognized as a critical parameter in both the early- and late-stage drug discovery. Therefore, in silico modeling of solubility has attracted extensive interests in recent years. Most previous studies have been limited in using relatively small data sets with limited diversity, which in turn limits the predictability of derived models. In this work, we present a support vector machines model for the binary classification of solubility by taking advantage of the largest known public data set that contains over 46 000 compounds with experimental solubility. Our model was optimized in combination with a reduction and recombination feature selection strategy. The best model demonstrated robust performance in both cross-validation and prediction of two independent test sets, indicating it could be a practical tool to select soluble compounds for screening, purchasing, and synthesizing. Moreover, our work may be used for comparative evaluation of solubility classification studies ascribe to the use of completely public resources.

In silico log P prediction for a large data set with support vector machines, radial basis neural networks and multiple linear regression

Oil/water partition coefficient (log P) is one of the key points for lead compound to be drug. In silico log P models based solely on chemical structures have become an important part of modern drug discovery. Here, we report support vector machines, radial basis function neural networks, and multiple linear regression methods to investigate the correlation between partition coefficient and physico-chemical descriptors for a large data set of compounds. The correlation coefficient r(2) between experimental and predicted log P for training and test sets by support vector machines, radial basis function neural networks, and multiple linear regression is 0.92, 0.90, and 0.88, respectively. The results show that non-linear support vector machines derives statistical models that have better prediction ability than those of radial basis function neural networks and multiple linear regression methods. This indicates that support vector machines can be used as an alternative modeling tool for quantitative structure-property/activity relationships studies.

Prediction of mutagenic toxicity by combination of Recursive Partitioning and Support Vector Machines

The study of prediction of toxicity is very important and necessary because measurement of toxicity is typically time-consuming and expensive. In this paper, Recursive Partitioning (RP) method was used to select descriptors. RP and Support Vector Machines (SVM) were used to construct structure-toxicity relationship models, RP model and SVM model, respectively. The performances of the two models are different. The prediction accuracies of the RP model are 80.2% for mutagenic compounds in MDL's toxicity database, 83.4% for compounds in CMC and 84.9% for agrochemicals in in-house database respectively. Those of SVM model are 81.4%, 87.0% and 87.3% respectively.