Quantitative structure-permeability relationships at various pH values for acidic and basic drugs and drug-like compounds

Trends in Molecular Properties, Bioavailability, and Permeability across the Bayer Compound Collection

For oral drugs, medicinal chemists aim to design compounds with high oral bioavailability, of which permeability is a key determinant. Taking advantage of >2000 compounds tested in rat bioavailability studies and >20,000 compounds tested in Caco2 assays at Bayer, we have examined the molecular properties governing bioavailability and permeability. In addition to classical parameters such as logD and molecular weight, we also investigated the relationship between calculated pK_a and permeability. We find that neutral compounds retain permeability up to a molecular weight limit of 700, while stronger acids and bases are restricted to weights of 400-500. We also investigate trends for common properties such as hydrogen bond donors and acceptors, polar surface area, aromatic ring count, and rotatable bonds, including compounds which exceed Lipinski's rule of five (Ro5). These property-structure relationships are combined to provide design guidelines for bioavailable drugs in both traditional and "beyond rule of 5" (bRo5) chemical space.

Development of artificial neural network models to predict the PAMPA effective permeability of new, orally administered drugs active against the coronavirus SARS-CoV-2

Responding to the pandemic caused by SARS-CoV-2, the scientific community intensified efforts to provide drugs effective against the virus. To strengthen these efforts, the "COVID Moonshot" project has been accepting public suggestions for computationally triaged, synthesized, and tested molecules. The project aimed to identify molecules of low molecular weight with activity against the virus, for oral treatment. The ability of a drug to cross the intestinal cell membranes and enter circulation decisively influences its bioavailability, and hence the need to optimize permeability in the early stages of drug discovery. In our present work, as a contribution to the ongoing scientific efforts, we employed artificial neural network algorithms to develop QSAR tools for modelling the PAMPA effective permeability (passive diffusion) of orally administered drugs. We identified a set of 61 features most relevant in explaining drug cell permeability and used them to develop a stacked regression ensemble model, subsequently used to predict the permeability of molecules included in datasets made available through the COVID Moonshot project. Our model was shown to be robust and may provide a promising framework for predicting the potential permeability of molecules not yet synthesized, thus guiding the process of drug design.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13721-023-00410-9.

Physicochemical Profile of Antiandrogen Drug Bicalutamide: Solubility, Distribution, Permeability

The pharmacologically relevant physicochemical properties of the antiandrogen drug bicalutamide (BCL) have been determined for the first time. Solubility in aqueous solution, 1-octanol, n-hexane, and ethanol was measured by the shake flask method in the temperature range of 293.15−313.15 K. The compound was shown to be poorly soluble in aqueous medium and n-hexane; at the same time, an essentially higher solubility in the alcohols was revealed. The following order of molar solubility was determined: ethanol > 1-octanol > water ≈ n-hexane. The solubility was correlated with the Van’t Hoff and Apelblat equations. Evaluation of the Hansen solubility parameters and the atomic group contribution approach of Hoftyzer and Van Krevelen demonstrated consistency with the experimental data and good potential adsorption of bicalutamide. The temperature dependences of the distribution coefficients in the 1-octanol/water and n-hexane/water two-phase systems were measured and discussed in the framework of the thermodynamic approach. The ∆logD parameter determined from the distribution experiment clearly demonstrated the preference of the lipophilic delivery pathways for the compound in the biological media. The overall thermodynamic analysis based on the solubility and distribution results of the present study and the sublimation characteristics published previously has been performed. To this end, the thermodynamic parameters of the dissolution, solvation, and transfer processes were calculated and discussed in view of the solute-solvent interactions. The permeation rate of BCL through the PermeaPad barrier was measured and compared with PAMPA permeability.

Using in vitro ADME data for lead compound selection: An emphasis on PAMPA pH 5 permeability and oral bioavailability

Membrane permeability plays an important role in oral drug absorption. Caco-2 and Madin-Darby Canine Kidney (MDCK) cell culture systems have been widely used for assessing intestinal permeability. Since most drugs are absorbed passively, Parallel Artificial Membrane Permeability Assay (PAMPA) has gained popularity as a low-cost and high-throughput method in early drug discovery when compared to high-cost, labor intensive cell-based assays. At the National Center for Advancing Translational Sciences (NCATS), PAMPA pH 5 is employed as one of the Tier I absorption, distribution, metabolism, and elimination (ADME) assays. In this study, we have developed a quantitative structure activity relationship (QSAR) model using our ∼6500 compound PAMPA pH 5 permeability dataset. Along with ensemble decision tree-based methods such as Random Forest and eXtreme Gradient Boosting, we employed deep neural network and a graph convolutional neural network to model PAMPA pH 5 permeability. The classification models trained on a balanced training set provided accuracies ranging from 71% to 78% on the external set. Of the four classifiers, the graph convolutional neural network that directly operates on molecular graphs offered the best classification performance. Additionally, an ∼85% correlation was obtained between PAMPA pH 5 permeability and in vivo oral bioavailability in mice and rats. These results suggest that data from this assay (experimental or predicted) can be used to rank-order compounds for preclinical in vivo testing with a high degree of confidence, reducing cost and attrition as well as accelerating the drug discovery process. Additionally, experimental data for 486 compounds (PubChem AID: 1645871) and the best models have been made publicly available (https://opendata.ncats.nih.gov/adme/).

Validating ADME QSAR Models Using Marketed Drugs

Problems with drug ADME are responsible for many clinical failures. By understanding the ADME properties of marketed drugs and modeling how chemical structure contributes to these inherent properties, we can help new projects reduce their risk profiles. Kinetic aqueous solubility, the parallel artificial membrane permeability assay (PAMPA), and rat liver microsomal stability constitute the Tier I ADME assays at the National Center for Advancing Translational Sciences (NCATS). Using recent data generated from in-house lead optimization Tier I studies, we update quantitative structure-activity relationship (QSAR) models for these three endpoints and validate in silico performance against a set of marketed drugs (balanced accuracies range between 71% and 85%). Improved models and experimental datasets are of direct relevance to drug discovery projects and, together with the prediction services that have been made publicly available at the ADME@NCATS web portal (https://opendata.ncats.nih.gov/adme/), provide important tools for the drug discovery community. The results are discussed in light of our previously reported ADME models and state-of-the-art models from scientific literature.Graphical Abstract[Figure: see text].

A role of flavonoids in cytochrome c-cardiolipin interactions

The processes preceding the detachment of cytochrome c (cyt c) from the inner mitochondrial membrane in intrinsic apoptosis involve peroxidation of cardiolipin (CL) catalyzed by cyt c-CL complex. In the present work, we studied the effect of 17 dietary flavonoids on the peroxidase activity of cyt c bound to liposomes. Specifically, we explored the relationship between peroxidase activity and flavonoids' (1) potential to modulate cyt c unfolding, (2) effect on the oxidation state of heme iron, (3) membrane permeability, (4) membrane binding energy, and (5) structure. The measurements revealed that flavones, flavonols, and flavanols were the strongest, while isoflavones were the weakest inhibitors of the oxidation. Flavonoids' peroxidase inhibition activity correlated positively with their potential to suppress Trp-59 fluorescence in cyt c as well as the number of OH groups. Hydrophilic flavonoids, such as catechin, having the lowest membrane permeability and the strongest binding with phosphocholine (PC) based on the quantum chemical calculations exhibited the strongest inhibition of Amplex Red (AR) peroxidation, suggesting a membrane-protective function of flavonoids at the surface. The results of the present research specify basic principles for the design of molecules that will control the catalytic oxidation of lipids in mitochondrial membranes. These principles take into account the number of hydroxyl groups and hydrophilicity of flavonoids.

Real-Time Parallel Artificial Membrane Permeability Assay Based on Supramolecular Fluorescent Artificial Receptors

Parallel artificial membrane permeability assay (PAMPA) is a screening tool for the evaluation of drug permeability across various biological membrane systems in a microplate format. In PAMPA, a drug candidate is allowed to pass through the lipid layer of a particular well during an incubation period of, typically, 10–16 h. In a second step, the samples of each well are transferred to a UV-Vis–compatible microplate and optically measured (applicable only to analytes with sufficient absorbance) or sampled by mass-spectrometric analysis. The required incubation period, sample transfer, and detection methods jointly limit the scalability of PAMPA to high-throughput screening format. We introduce a modification of the PAMPA method that allows direct fluorescence detection, without sample transfer, in real time (RT-PAMPA). The method employs the use of a fluorescent artificial receptor (FAR), composed of a macrocycle in combination with an encapsulated fluorescent dye, administered in the acceptor chamber of conventional PAMPA microplates. Because the detection principle relies on the molecular recognition of an analyte by the receptor and the associated fluorescence response, concentration changes of any analyte that binds to the receptor can be monitored (molecules with aromatic residues in the present example), regardless of the spectroscopic properties of the analyte itself. Moreover, because the fluorescence of the (upper) acceptor well can be read out directly by fluorescence in a microplate reader, the permeation of the drug through the planar lipid layer can be monitored in real time. Compared with the traditional assay, RT-PAMPA allows not only quantification of the permeability characteristics but also rapid differentiation between fast and slow diffusion events.

Modelling of antiproliferative activity measured in HeLa cervical cancer cells in a series of xanthene derivatives

Cancer remains one of the leading causes of death in humans, and new drug substances are therefore being developed. Thus, the anti-cancer activity of xanthene derivatives has become an important topic in the development of new and potent anti-cancer drug substances. Previously published novel series of xanthen-3-one and xanthen-1,8-dione derivatives have been synthesized in one of our laboratories and showed anti-proliferative activity in HeLa cancer cell lines. This series serves as a good basis to develop quantitative structure-activity relationship (QSAR), to study the relations between anti-proliferative activity and chemical structures. A QSAR model has been derived that relies only on two-dimensional molecular descriptors, providing mechanistic insight into the anti-proliferative activity of xanthene derivatives. The model is validated internally and externally and additionally with the set of inactive compounds of the original data, confirming model applicability for the design and discovery of novel xanthene derivatives. The QSAR model is available at the QsarDB repository (http://dx.doi.10.15152/QDB.237).

In Silico Prediction of PAMPA Effective Permeability Using a Two-QSAR Approach

Oral administration is the preferred and predominant route of choice for medication. As such, drug absorption is one of critical drug metabolism and pharmacokinetics (DM/PK) parameters that should be taken into consideration in the process of drug discovery and development. The cell-free in vitro parallel artificial membrane permeability assay (PAMPA) has been adopted as the primary screening to assess the passive diffusion of compounds in the practical applications. A classical quantitative structure–activity relationship (QSAR) model and a machine learning (ML)-based QSAR model were derived using the partial least square (PLS) scheme and hierarchical support vector regression (HSVR) scheme to elucidate the underlying passive diffusion mechanism and to predict the PAMPA effective permeability, respectively, in this study. It was observed that HSVR executed better than PLS as manifested by the predictions of the samples in the training set, test set, and outlier set as well as various statistical assessments. When applied to the mock test, which was designated to mimic real challenges, HSVR also showed better predictive performance. PLS, conversely, cannot cover some mechanistically interpretable relationships between descriptors and permeability. Accordingly, the synergy of predictive HSVR and interpretable PLS models can be greatly useful in facilitating drug discovery and development by predicting passive diffusion.

Prediction of Membrane Permeation of Drug Molecules by Combining an Implicit Membrane Model with Machine Learning

Lipid membrane permeation of drug molecules was investigated with Heterogeneous Dielectric Generalized Born (HDGB)-based models using solubility-diffusion theory and machine learning. Free energy profiles were obtained for neutral molecules by the standard HDGB and Dynamic HDGB (DHDGB) to account for the membrane deformation upon insertion of drugs. We also obtained hybrid free energy profiles where the neutralization of charged molecules was taken into account upon membrane insertion. The evaluation of the predictions was done against experimental permeability coefficients from Parallel Artificial Membrane Permeability Assays (PAMPA), and effects of partial charge sets, CGenFF, AM1-BCC, and OPLS, on the performance of the predictions were discussed. (D)HDGB-based models improved the predictions over the two-state implicit membrane models, and partial charge sets seemed to have a strong impact on the predictions. Machine learning increased the accuracy of the predictions, although it could not outperform the physics-based approach in terms of correlations.

pH-permeability profiles for drug substances: Experimental detection, comparison with human intestinal absorption and modelling

The influence of pH on human intestinal absorption is frequently not considered in early drug discovery studies in the modelling and subsequent prediction of intestinal absorption for drug candidates. To bridge this gap, in this study, experimental membrane permeability data were measured for current and former drug substances with a parallel artificial membrane permeability assay (PAMPA) at different pH values (3, 5, 7.4 and 9). The presented data are in good agreement with human intestinal absorption, showing a clear influence of pH on the efficiency of intestinal absorption. For the measured data, simple and general quantitative structure-activity relationships (QSARs) were developed for each pH that makes it possible to predict the pH profiles for passive membrane permeability (i.e., a pH-permeability profile), and these predictions coincide well with the experimental data. QSARs are also proposed for the data series of highest and intrinsic membrane permeability. The molecular descriptors in the models were analysed and mechanistically related to the interaction pattern of permeability in membranes. In addition to the regression models, classification models are also proposed. All models were successfully validated and blind tested with external data. The models are available in the QsarDB repository (http://dx.doi.org/10.15152/QDB.203).

Learning how to use IAM chromatography for predicting permeability

The interest for IAM (Immobilized Artificial Membranes) chromatography in the prediction of drug permeability is increasing. Here we firstly set-up a dataset of 253 molecules including neutral and ionized drugs and few organic compounds for which we either measured or retrieved from the literature IAM.PC.DD2 log K_w^IAM data. Then we applied block relevance (BR) analysis to extract from PLS models the relative contribution of intermolecular forces governing log K_w^IAM and Δlog K_w^IAM (a combined descriptor calculated from log K_w^IAM). Finally, the relationship between log K_w^IAM, Δlog K_w^IAM and passive permeability determined in both PAMPA and MDCK-LE systems was looked for. Models provided the basis for a rational application of IAM chromatography in permeability prediction.

Chemical structure and correlation analysis of HIV-1 NNRT and NRT inhibitors and database-curated, published inhibition constants with chemical structure in diverse datasets

Human immunodeficiency virus (HIV-1) reverse transcriptase is a major target for designing anti-HIV drugs. Developed inhibitors are divided into non-nucleoside analog reverse-transcriptase inhibitors (NNRTIs) and nucleoside analog reverse-transcriptase inhibitors (NRTIs) depending on their mechanism. Given that many inhibitors have been studied and for many of them binding affinity constants have been calculated, it is beneficial to analyze the chemical landscape of these families of inhibitors and correlate these inhibition constants with molecular structure descriptors. For this, the HIV-1 RT data was retrieved from the ChEMBL database, carefully curated, and original literature verified, grouped into NRTIs and NNRTIs, analyzed using a hierarchical scaffold classification method and modelled with best multi-linear regression approach. Analysis of the HIV-1 NNRTIs subset results in ten different common structural parent types of oxazepanone, piperazinone, pyrazine, oxazinanone, diazinanone, pyridine, pyrrole, diazepanone, thiazole, and triazine. The same analysis for HIV-1 NRTIs groups structures into four different parent types of uracil, pyrimide, pyrimidione, and imidazole. Each scaffold tree corresponding to the parent types has been carefully analyzed and examined, and changes in chemical structure favorable to potency and stability are highlighted. For both subsets, descriptive and predictive QSAR models are derived, discussed and externally validated, revealing general trends in relationships between molecular structure and binding affinity constants in structurally diverse datasets. Data and QSAR models are available at the QsarDB repository (http://dx.doi.org/10.15152/QDB.202).

Quantitative structure-permeability relationships at various pH values for neutral and amphoteric drugs and drug-like compounds

Human intestinal absorption is a key property for orally administered drugs and is dependent on pH. This study focuses on neutral and amphoteric compounds and their membrane permeabilities across the range of pH values found in the human intestine. The membrane permeability values for 15 neutral and 60 amphoteric compounds at pH 3, 5, 7.4 and 9 were measured using the parallel artificial membrane permeability assay (PAMPA). For each data series the quantitative structure-permeability relationships were developed and analysed. The results show that the membrane permeability of neutral compounds is attributed to a single structural characteristic, the hydrogen bond donor ability. Amphoteric compounds are more complex because of their chemical constitution, and therefore require three-parameter models to describe and predict membrane permeability. Analysis of the models for amphoteric compounds reveals that membrane permeability depends on multiple structural characteristics: the partition coefficient, hydrogen bond properties and the shape of the molecules. In addition to conventional validation strategies, two external compounds (isradipine and omeprazole) were tested and revealed very good agreement of pH profiles between experimental and predicted membrane permeability for all of the developed models. Selected QSAR models are available at the QsarDB repository (http://dx.doi.org/10.15152/QDB.184).