Development and evaluation of two-parameter linear free energy models for the prediction of human skin permeability coefficient of neutral organic chemicals

Machine learning for skin permeability prediction: random forest and XG boost regression

Background: Machine learning algorithms that can quickly and easily estimate skin permeability (Kp) are increasingly being used in drug delivery research. The linear free energy relationship (LFER) developed by Abraham is a practical technique for predicting Kp. The permeability coefficients and Abraham solute descriptor values for 175 organic compounds have been documented in the scientific literature.Purpose: The purpose of this project was to use a publicly available dataset to make skin permeability predictions using the random forest and XBoost regression techniques.Methods: We employed Pandas-based methods in JupyterLab to predict permeability coefficient (Kp) from solute descriptors (excess molar refraction [E], combined dipolarity/polarizability [S], overall solute hydrogen bond acidity and basicity [A and B], and the McGowan's characteristic molecular volume [V]).Results: The random forest and XG Boost regression models established statistically significant association between the descriptors and the skin permeability coefficient.

Delivery of small hydrophilic molecules across the stratum corneum: Identification of model systems and parameters to study topical delivery of free amino acids

Free amino acids (FAAs) constitute the largest component (∼40 %) of the so-called natural moisturizing factors of the skin. Their level declines in dry skin conditions and one strategy to overcome this problem may involve the topical delivery of FAAs through appropriate strategy. The objective of the present study was therefore to identify alternative skin models and study the corneocyte-water partition coefficients (KCOR/W) and permeation coefficient (KP) of 18 FAAs. The KCOR/W was studied using standard protocols and the permeation studies were conducted using Franz diffusion cell. The results indicate that the FAAs have high partitioning behavior to the corneocytes. The KCOR/W values of the human COR and that of pig ear skin were better correlated with each other than that of keratin isolated from chicken feathers. The presence of lipid in the stratum corneum (SC), initial concentration of the FAAs, and permeation enhancers affect the KCOR/W. The FAAs have low permeation into the SC which suggests the need for permeation enhancers in designing dosage form containing these compounds. Even though the investigated mathematical models show good prediction of the Kp values, better prediction could be obtained by considering factors such as the possible entrapment of the FAAs by the CORs.

Comparison of in vitro membrane permeabilities of diverse environmental chemicals with in silico predictions

The parallel artificial membrane permeability assay (PAMPA) is widely used for estimating biomembrane permeabilities of experimental drugs in pharmaceutical research. However, there are few reports of studies using PAMPA to measure membrane permeabilities of chemicals of environmental concern (CECs) outside the pharmaceutical domain, many of which differ substantially from drugs in their physicochemical properties. We applied PAMPA methods simulating gastrointestinal (PAMPA-GIT) and blood-brain barrier (PAMPA-BBB) membranes under consistent conditions to 51 CECs, including some pharmaceuticals. A backward stepwise multivariate linear regression was implemented to explore the correlation between the differences of measured permeabilities from PAMPA-GIT and PAMPA-BBB and Abraham solute descriptors. In addition, a previously reported in silico model was evaluated by comparing predicted and measured permeability results. PAMPA-GIT and PAMPA-BBB experimental permeability results agreed relatively well. The backward stepwise multivariate linear regression identified excess molar refraction and polarizability to be significant at the 0.10 level in predicting the differences between PAMPA-GIT and PAMPA-BBB. The in silico model performed well - with predicted permeability of most compounds within two-fold of experimentally measured values. We found that CECs pose experimental challenges to the PAMPA method in terms of having lower solubility and lower stability compared to most drugs.

Chemical Space Covered by Applicability Domains of Quantitative Structure-Property Relationships and Semiempirical Relationships in Chemical Assessments

A significant number of chemicals registered in national and regional chemical inventories require assessments of their potential "hazard" concerns posed to humans and ecological receptors. This warrants knowledge of their partitioning and reactivity properties, which are often predicted by quantitative structure-property relationships (QSPRs) and other semiempirical relationships. It is imperative to evaluate the applicability domain (AD) of these tools to ensure their suitability for assessment purpose. Here, we investigate the extent to which the ADs of commonly used QSPRs and semiempirical relationships cover seven partitioning and reactivity properties of a chemical "space" comprising 81,000+ organic chemicals registered in regulatory and academic chemical inventories. Our findings show that around or more than half of the chemicals studied are covered by at least one of the commonly used QSPRs. The investigated QSPRs demonstrate adequate AD coverage for organochlorides and organobromines but limited AD coverage for chemicals containing fluorine and phosphorus. These QSPRs exhibit limited AD coverage for atmospheric reactivity, biodegradation, and octanol-air partitioning, particularly for ionizable organic chemicals compared to nonionizable ones, challenging assessments of environmental persistence, bioaccumulation capability, and long-range transport potential. We also find that a predictive tool's AD coverage of chemicals depends on how the AD is defined, for example, by the distance of a predicted chemical from the centroid of the training chemicals or by the presence or absence of structural features.

Novel Models for Accurate Estimation of Air-Blood Partitioning: Applications to Individual Compounds and Complex Mixtures of Neutral Organic Compounds

The air-blood partition coefficient (Kab) is extensively employed in human health risk assessment for chemical exposure. However, current Kab estimation approaches either require an extensive number of parameters or lack precision. In this study, we present two novel and parsimonious models to accurately estimate Kab values for individual neutral organic compounds, as well as their complex mixtures. The first model, termed the GC×GC model, was developed based on the retention times of nonpolar chemical analytes on comprehensive two-dimensional gas chromatography (GC×GC). This model is unique in its ability to estimate the Kab values for complex mixtures of nonpolar organic chemicals. The GC×GC model successfully accounted for the Kab variance (R2 = 0.97) and demonstrated strong prediction power (RMSE = 0.31 log unit) for an independent set of nonpolar chemical analytes. Overall, the GC×GC model can be used to estimate Kab values for complex mixtures of neutral organic compounds. The second model, termed the partition model (PM), is based on two types of partition coefficients: octanol to water (Kow) and air to water (Kaw). The PM was able to effectively account for the variability in Kab data (n = 344), yielding an R2 value of 0.93 and root-mean-square error (RMSE) of 0.34 log unit. The predictive power and explanatory performance of the PM were found to be comparable to those of the parameter-intensive Abraham solvation models (ASMs). Additionally, the PM can be integrated into the software EPI Suite, which is widely used in chemical risk assessment for initial screening. The PM provides quick and reliable estimation of Kab compared to ASMs, while the GC×GC model is uniquely suited for estimating Kab values for complex mixtures of neutral organic compounds. In summary, our study introduces two novel and parsimonious models for the accurate estimation of Kab values for both individual compounds and complex mixtures.

Addressing new chemicals of emerging concern (CECs) in an indoor office

Indoor pollutants change over time and place. Exposure to hazardous organics is associated with adverse health effects. This work sampled gaseous organics by Tenax TA tubes in two indoor rooms, i.e., an office set as samples, and the room of chassis dynamometer (RCD) set as backgrounds. Compounds are analyzed by a thermal desorption comprehensive two-dimensional gas chromatography-quadrupole mass spectrometer (TD-GC × GC-qMS). Four new chemicals of emerging concern (CECs) are screened in 469 organics quantified. We proposed a three-step pipeline for CECs screening utilizing GC × GC including 1) non-target scanning of organics with convincing molecular structures and quantification results, 2) statistical analysis between samples and backgrounds to extract useful information, and 3) pixel-based property estimation to evaluate the contamination potential of addressed chemicals. New CECs spotted in this work are all intermediate volatility organic compounds (IVOCs), containing mintketone, isolongifolene, β-funebrene, and (5α)-androstane. Mintketone and sesquiterpenes may be derived from the use of volatile chemical products (VCPs), while (5α)-androstane is probably human-emitted. The occurrence and contamination potential of the addressed new CECs are reported for the first time. Non-target scanning and the measurement of IVOCs are of vital importance to get a full glimpse of indoor organics.

Immobilized Keratin HPLC Stationary Phase-A Forgotten Model of Transdermal Absorption: To What Molecular and Biological Properties Is It Relevant?

Chromatographic retention data collected on immobilized keratin (KER) or immobilized artificial membrane (IAM) stationary phases were used to predict skin permeability coefficient (log K_p) and bioconcentration factor (log BCF) of structurally unrelated compounds. Models of both properties contained, apart from chromatographic descriptors, calculated physico-chemical parameters. The log K_p model, containing keratin-based retention factor, has slightly better statistical parameters and is in a better agreement with experimental log K_p data than the model derived from IAM chromatography; both models are applicable primarily to non-ionized compounds.Based on the multiple linear regression (MLR) analyses conducted in this study, it was concluded that immobilized keratin chromatographic support is a moderately useful tool for skin permeability assessment.However, chromatography on immobilized keratin may also be of use for a different purpose-in studies of compounds' bioconcentration in aquatic organisms.

Reliable prediction of sensory irritation threshold values of organic compounds using new models based on linear free energy relationships and GC×GC retention parameters

The human sensory irritation threshold (SIT) is an important biochemical parameter for the exposure assessment of organic air pollutants. First, we recalibrated the Abraham solvation models (ASMs) for 9 SIT endpoints by curating 720 individual experimental SIT values to find an accurate and parsimonious ASM variant, which exhibited root mean square error (RMSE) = 0.174-0.473 log unit. Second, we report linear free energy relationships - henceforth called partition models (PMs) - which exploit the correlations of 9 SIT endpoints with the linear combinations of partition coefficients for octanol-water and air-water systems showing RMSE = 0.221-0.591 log unit. These PMs can easily be integrated into widely used EPI-Suite™ screening tool. The explanatory and predictive performance of PMs were like parameter-intensive ASMs. Third, we present GC × GC models that are based on the retention times of the nonpolar analytes on the comprehensive two-dimensional gas chromatography (GC × GC), which successfully described the SIT variance (R²=0.959-0.996) and depicted a strong predictive power (RMSE = 0.359-0.660 log unit) for an independent set of nonpolar analytes. Taken together, PMs allow easy SIT screening of organic chemicals compared to ASMs. Unlike ASMs, our GC × GC models can be applied to estimate SIT of complex nonpolar mixtures.

Predicting skin permeability using HuskinDB

A freely accessible database has recently been released that provides measurements available in the literature on human skin permeation data, known as the ‘Human Skin Database – HuskinDB’. Although this database is extremely useful for sourcing permeation data to help with toxicity and efficacy determination, it cannot be beneficial when wishing to consider unlisted, or novel compounds. This study undertakes analysis of the data from within HuskinDB to create a model that predicts permeation for any compound (within the range of properties used to create the model). Using permeability coefficient (K_p) data from within this resource, several models were established for K_p values for compounds of interest by varying the experimental parameters chosen and using standard physicochemical data. Multiple regression analysis facilitated creation of one particularly successful model to predict K_p through human skin based only on three chemical properties. The model transforms the dataset from simply a resource of information to a more beneficial model that can be used to replace permeation testing for a wide range of compounds.

Exploring the role of octanol-water partition coefficient and Henry's law constant in predicting the lipid-water partition coefficients of organic chemicals

Partition coefficients for storage lipid-water (logKlw) and phospholipid-water (logKpw) phases are key parameters to understand the bioaccumulation and toxicity of organic contaminants. However, the published experimental databases of these properties are dwarfs and current estimation approaches are cumbersome. Here, we present partition models that exploit the correlations of logKlw, and of logKpw with the linear combinations of the octanol-water partition coefficient (logKow) and the dimensionless Henry's law constant (air-water partition coefficient, logKaw). The calibrated partition models successfully describe the variations in logKlw data (n = 305, R2 = 0.971, root-mean-square-error (rmse) = 0.375), and in logKpw data (n = 131, R2 = 0.953, rmse = 0.413). With the inputs of logKow and logKaw estimated from the U.S. EPA's EPI Suite, our models of logKlw and logKpw have exhibited rmse = 0.52 with respect to experimental values indicating suitability of these models for inclusion in the EPI Suite. Our models perform similar to or better than the previously reported models such as one parameter partition models, Abraham solvation models, and models based on quantum-chemical calculations. Taken together, our models are robust, easy-to-use, and provide insight into variations of logKlw and logKpw in terms of hydrophobicity and volatility trait of chemicals.

IAM Chromatographic Models of Skin Permeation

Chromatographic retention factor log k_IAM obtained from IAM HPLC chromatography with buffered aqueous mobile phases and calculated molecular descriptors (surface area—S_a; molar volume—V_M; polar surface area—PSA; count of freely rotable bonds—FRB; H-bond acceptor count—HA; energy of the highest occupied molecular orbital—E_HOMO; energy of the lowest unoccupied orbital—E_LUMO; and polarizability—α) obtained for a group of 160 structurally unrelated compounds were tested in order to generate useful models of solutes’ skin permeability coefficient log K_p. It was established that log k_IAM obtained in the conditions described in this study is not sufficient as a sole predictor of the skin permeability coefficient. Simple put, potentially useful models based on log k_IAM and readily available calculated descriptors, accounting for 85 to 91% of the total variability, were generated using Multiple Linear Regression (MLR).The models proposed in the study were tested on a group of 20 compounds with known experimental log K_p values.

Investigation of partition coefficients and fingerprints of atmospheric gas- and particle-phase intermediate volatility and semi-volatile organic compounds using pixel-based approaches

Ambient gas- and particle-phase intermediate volatility and semi-volatile organic compounds (I/SVOCs) of Beijing were analyzed by a thermal desorption comprehensive two-dimensional gas chromatography quadrupole mass spectrometry (TD-GC × GC-qMS). A pixel-based scheme combing the integration-based approach was applied for partition coefficients estimation and fingerprints identification. Blob-by-blob recognition was firstly utilized to characterize I/SVOCs from the molecular level. 412 blobs in gas-phase and 460 blobs in particle-phase were resolved, covering a total response of 47.5% and 43.5%. A large pool of I/SVOCs was found with a large diversity of chemical classes in both gas- and particle-phase. Acids (8.5%), b-alkanes (5.8%), n-alkanes (C8-C25, 5.3%), and aromatics (4.4%) were dominant in gas-phase while esters (7.0%, including volatile chemical product compounds, VCPs), n-alkanes (C9-C34, 5.7%), acids (4.6%), and siloxanes (3.6%) were abundant in particle-phase. Air pollutants were then evaluated by a two-parameter linear free energy relationship (LFER) model, which could be further implemented in the two-dimensional volatility basis set (2D-VBS) model. Multiway principal component analysis (MPCA) and partial least squares-discriminant analysis (PLS-DA) implied that naphthalenes, phenol, propyl-benzene isomers, and oxygenated volatile organic compounds (OVOCs) were key components in the gas-phase under different pollution levels. This work gives more insight into property estimation and fingerprints identification for complex ambient samples.