Enhancing Multi-species Liver Microsomal Stability Prediction through Artificial Intelligence

Liver microsomal stability, a crucial aspect of metabolic stability, significantly impacts practical drug discovery. However, current models for predicting liver microsomal stability are based on limited molecular information from a single species. To address this limitation, we constructed the largest public database of compounds from three common species: human, rat, and mouse. Subsequently, we developed a series of classification models using both traditional descriptor-based and classic graph-based machine learning (ML) algorithms. Remarkably, the best-performing models for the three species achieved Matthews correlation coefficients (MCCs) of 0.616, 0.603, and 0.574, respectively, on the test set. Furthermore, through the construction of consensus models based on these individual models, we have demonstrated their superior predictive performance in comparison with the existing models of the same type. To explore the similarities and differences in the properties of liver microsomal stability among multispecies molecules, we conducted preliminary interpretative explorations using the Shapley additive explanations (SHAP) and atom heatmap approaches for the models and misclassified molecules. Additionally, we further investigated representative structural modifications and substructures that decrease the liver microsomal stability in different species using the matched molecule pair analysis (MMPA) method and substructure extraction techniques. The established prediction models, along with insightful interpretation information regarding liver microsomal stability, will significantly contribute to enhancing the efficiency of exploring practical drugs for development.

CMMS-GCL: cross-modality metabolic stability prediction with graph contrastive learning

MOTIVATION

Metabolic stability plays a crucial role in the early stages of drug discovery and development. Accurately modeling and predicting molecular metabolic stability has great potential for the efficient screening of drug candidates as well as the optimization of lead compounds. Considering wet-lab experiment is time-consuming, laborious, and expensive, in silico prediction of metabolic stability is an alternative choice. However, few computational methods have been developed to address this task. In addition, it remains a significant challenge to explain key functional groups determining metabolic stability.

RESULTS

To address these issues, we develop a novel cross-modality graph contrastive learning model named CMMS-GCL for predicting the metabolic stability of drug candidates. In our framework, we design deep learning methods to extract features for molecules from two modality data, i.e. SMILES sequence and molecule graph. In particular, for the sequence data, we design a multihead attention BiGRU-based encoder to preserve the context of symbols to learn sequence representations of molecules. For the graph data, we propose a graph contrastive learning-based encoder to learn structure representations by effectively capturing the consistencies between local and global structures. We further exploit fully connected neural networks to combine the sequence and structure representations for model training. Extensive experimental results on two datasets demonstrate that our CMMS-GCL consistently outperforms seven state-of-the-art methods. Furthermore, a collection of case studies on sequence data and statistical analyses of the graph structure module strengthens the validation of the interpretability of crucial functional groups recognized by CMMS-GCL. Overall, CMMS-GCL can serve as an effective and interpretable tool for predicting metabolic stability, identifying critical functional groups, and thus facilitating the drug discovery process and lead compound optimization.

AVAILABILITY AND IMPLEMENTATION

The code and data underlying this article are freely available at https://github.com/dubingxue/CMMS-GCL.

Ensemble modeling with machine learning and deep learning to provide interpretable generalized rules for classifying CNS drugs with high prediction power

The trade-off between a machine learning (ML) and deep learning (DL) model's predictability and its interpretability has been a rising concern in central nervous system-related quantitative structure-activity relationship (CNS-QSAR) analysis. Many state-of-the-art predictive modeling failed to provide structural insights due to their black box-like nature. Lack of interpretability and further to provide easy simple rules would be challenging for CNS-QSAR models. To address these issues, we develop a protocol to combine the power of ML and DL to generate a set of simple rules that are easy to interpret with high prediction power. A data set of 940 market drugs (315 CNS-active, 625 CNS-inactive) with support vector machine and graph convolutional network algorithms were used. Individual ML/DL modeling methods were also constructed for comparison. The performance of these models was evaluated using an additional external dataset of 117 market drugs (42 CNS-active, 75 CNS-inactive). Fingerprint-split validation was adopted to ensure model stringency and generalizability. The resulting novel hybrid ensemble model outperformed other constituent traditional QSAR models with an accuracy of 0.96 and an F1 score of 0.95. With the power of the interpretability provided with this protocol, our model laid down a set of simple physicochemical rules to determine whether a compound can be a CNS drug using six sub-structural features. These rules displayed higher classification ability than classical guidelines, with higher specificity and more mechanistic insights than just for blood-brain barrier permeability. This hybrid protocol can potentially be used for other drug property predictions.

Characterization of phase I and phase II metabolites of hop (Humulus lupulus L.) bitter acids: In vitro and in vivo metabolic profiling by UHPLC-Q-Orbitrap

Bitter acids are a class of prenylated phloroglucinol derivatives present in Humulus lupulus L., known for their multiple healthy properties, nevertheless, research regarding their metabolism and stability is lacking. This study was aimed to elucidate the metabolic stability of hop α- and β-acids and characterize I and II phase metabolites in vitro and in vivo. For this purpose, an ultra high performance liquid chromatography-high resolution mass spectrometry (UHPLC-HRMS) method was developed and validated. Mice liver microsomes were used to assess metabolic stability; in vitro t_1/2 and clearance values were calculated, showing a moderate metabolism for α-acids (_avgt_1/2: 120.01 min, _avgCL_int 11.96 μL/min/mg), while β-acids were metabolized faster (_avgt_1/2: 103.01 min, _avgCL_int: 13.83 μL/min/mg). I and II phase metabolites were characterized both in in vitro, and in vivo, in mouse plasma and urine after oral administration. A combined full scan/data dependent/precursor ion list-triggered neutral loss (FS/dd-MS²/PIL-tNL) strategy was used to detect unknown and expected metabolites. As a result, 33 compounds were detected, including novel metabolites, such as 9 potential oxidized metabolites of humulones (M6-M14), and 10 glucuronide conjugates of α-acids, comprising 7 glucuronide derivatives of oxidized phase I metabolites (M26-M32). The proposed method extends the current knowledge regarding metabolization of hop α- and β-acids and could be applied for the comprehension of the metabolic fate of this class of compounds in different species, as well as for in vivo pharmacokinetic studies.