Next Article in Journal
Acute Stress Regulates Sex-Related Molecular Responses in the Human Jejunal Mucosa: Implications for Irritable Bowel Syndrome
Next Article in Special Issue
Endothelial Specific Deletion of Autotaxin Improves Stroke Outcomes
Previous Article in Journal
The Breakthroughs and Caveats of Using Human Pluripotent Stem Cells in Modeling Alzheimer’s Disease
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Computational Approach to Drug Penetration across the Blood-Brain and Blood-Milk Barrier Using Chromatographic Descriptors

Department of Analytical Chemistry, Faculty of Pharmacy, Medical University of Lodz, 90-419 Lodz, Poland
*
Author to whom correspondence should be addressed.
Cells 2023, 12(3), 421; https://doi.org/10.3390/cells12030421
Submission received: 30 November 2022 / Revised: 16 January 2023 / Accepted: 25 January 2023 / Published: 27 January 2023

Abstract

:
Drug penetration through biological barriers is an important aspect of pharmacokinetics. Although the structure of the blood-brain and blood-milk barriers is different, a connection can be found in the literature between drugs entering the central nervous system (CNS) and breast milk. This study was created to reveal such a relationship with the use of statistical modelling. The basic physicochemical properties of 37 active pharmaceutical compounds (APIs) and their chromatographic retention data (TLC and HPLC) were incorporated into calculations as molecular descriptors (MDs). Chromatography was performed in a thin layer format (TLC), where the plates were impregnated with bovine serum albumin to mimic plasma protein binding. Two columns were used in high performance liquid chromatography (HPLC): one with immobilized human serum albumin (HSA), and the other containing an immobilized artificial membrane (IAM). Statistical methods including multiple linear regression (MLR), cluster analysis (CA) and random forest regression (RF) were performed with satisfactory results: the MLR model explains 83% of the independent variable variability related to CNS bioavailability; while the RF model explains up to 87%. In both cases, the parameter related to breast milk penetration was included in the created models. A significant share of reversed-phase TLC retention values was also noticed in the RF model.

1. Introduction

The bioavailability of the blood-brain barrier (BBB) is denoted by the BBB+/− index [1,2]. The classification of BBB+ and BBB− types can be based on calculated log BB values and their comparison with in vivo or in vitro biological test data. Log BB is the most commonly used measure of drug penetration into the central nervous system (CNS), and it is defined as the logarithm of the ratio of drug concentration in the brain to drug concentration in blood [3,4]. Biological experiments determining the level of blood-brain distribution are the best source of information; however, these are extremely time-consuming, costly, and difficult-to-access experiments that require extensive screening.
Another classification which can be found in the literature is CNS+/− [5]. The drug types CNS+ and CNS− do not correspond to the BBB+/− index. CNS inactive drugs include both those that cross the BBB (without appropriate biological targets) and those that do not penetrate the barriers [6,7]. It is more certain to identify BBB+ drugs among CNS+ compounds than to identify BBB- among CNS− [6,8]. There are also studies [9] for which the basis for BBB+ classification was a higher concentration of the drug in the brain than in the blood.
The excretion into breast milk can be described by the M/P parameter, which is the ratio of the drug concentration in milk and plasma. Determining this parameter in vivo is even more difficult than log BB, while its knowledge may increase the safety of pharmacotherapy in breastfeeding women. Difficulties in obtaining the M/P ratio in vivo are mainly due to ethical reasons and the danger of affecting the breastfed infant, therefore most approaches to the subject take the form of case studies. Only a few results are obtained for individual drugs and these may not reflect what is happening in the scale of the entire population. The problem is even greater with drugs newly introduced to the market [10]. The computational methods could be helpful in the estimation of M/P ratios based on API properties and general principles of permeation through biological membranes.
Blood-brain and blood milk-barrier have many differences; the blood-milk barrier is simpler, the main obstacle is the lactocyte membrane. Blood-brain barrier exhibits an additional enzymatic activity and a reinforced structure (Table 1) [11]. CNS-active drugs have been the subject of case studies for safety in the breastfed child; most of them willingly pass into breast milk, and some even have M/P ratios above 1 (e.g., citalopram, escitalopram or sertraline), which indicates a higher concentration in milk than in plasma [12,13].
A small group of active pharmaceutical ingredients APIs (37 compounds) with known biological effects were analysed. The main aim of the study is to investigate the use of chromatographic systems describing the protein binding capacity as analytical models for two phenomena—drug distribution to the CNS and breast milk, and thus to find a relationship between the bioavailability of APIs in the CNS, and the dangers of their use in breastfeeding women. The chemometric analyses were based on the result of an earlier published work [14] in which it was proposed to investigate the aforementioned relationship; however, compounds 137 were never subject to a separate computational experiment. The structures of 37 APIs, along with their applied physicochemical and biological properties, are included in the Supplementary Material (Tables S1–S5).
Plasma protein binding (PPB) is one of the most important pharmacokinetic aspects limiting permeation across membranes and biological barriers. Human serum albumin is the most abundant plasma protein and is responsible for the greatest extent for PPB. Therefore, chromatography-based laboratory models have been developed to estimate the level of HSA binding and its effect on CNS entry. An observation was made of the effect of stationary phase modification with bovine serum albumin (BSA) in thin layer chromatography (TLC) and human serum albumin (HSA) in high performance liquid chromatography (HPLC). BSA is considered to be a cheaper alternative to HSA, with similar binding and physicochemical properties [15].
Due to the different routes of entry of drugs into the brain and milk, an immobilized artificial membrane (IAM) column was also used. IAM chromatography is based on the permeability of a compound through a biological membrane; the phospholipid monolayer is bound to a silica base [16].
Two sets of thin layer affinity chromatography were applied: in normal (NP) and reverse phase (RP), and two mentioned columns in the HPLC technique.

2. Materials and Methods

The results of chromatographic experiments, as well as molecular descriptors (MDs) connected to the physicochemical properties of APIs were introduced into the research (Table 2). Physicochemical properties were calculated (Hyperchem; HyperCube Inc., Waterloo, ON, Canada, 2002 and ACD Labs; Advanced Chemistry Development Inc., Toronto, ON, Canada) or collected from online databases: Drugbank [17] and Chembl [18]. Chromatographic data: Rf and log k parameters were obtained experimentally; Rf (retardation factor) is the ratio of the mobile phase distance travelled to the API distance travelled on the TLC plate. HPLC parameter: log k (logarithm of the retention factor), where k is the ratio between the amount of API in the stationary phase to its amount in the mobile phase, and it was obtained from the equation k = (tR − tM)/tM (tR—retention time; tM—void time). TLC and HPLC experiments are detailed in Appendix B.

3. Results

3.1. Correlation Analysis and Multiple Linear Regression

The limit values of the log BB index are different for many of the proposed BBB penetration prediction models [8,9,24]. The log BB = −0.52 is the logical division into BBB+ and BBB–. The logarithmic value of −0.52 corresponds to a 30% brain–to–plasma ratio of the compound [25]. The optimal classification threshold is usually 0 to −1 [25,26]. Studies conducted earlier defined the BBB+ limit as log BB ≥ −0.9 [27]. The calculations in this study use all the rates of BBB penetration through the 37 observed drugs. The parameter used (B2) is based on the log BB and describes the division into BBB +/− [19], and the following limit values were also used: B2 > −0.90 and B2 > −0.52. Another applied descriptor: CNS +/− is described in the literature [17].
In order to combine two pharmacokinetic phenomena—the drug penetration of blood–milk and the blood–brain barrier, the observation of molecular descriptors of 37 compounds was introduced. MDs which, according to the literature, usually affect BBB +/− and M/P are: binding to plasma proteins—PB, the acid–base nature of the compounds—acid/base and PhCharge—physiological charge. The correlation matrix that describes the relationships between these indicators in the study group is presented in Table 3.
The connection of B2, B2 > −0.90 descriptors and PB is clearly visible. Additionally, it was assumed that the best correlation of the binding strength with the plasma proteins is connected with the B2 descriptor. On its basis, the availability of drugs to the brain can be determined. The choice between B2 > −0.90 and B2 > −0.52 is subjective, so the range of drug properties related to the B2 measure was investigated.
In milk penetration correlation analyses (Table 4), APIs indicated as CNS− were removed from a group of 37 cases. Only potentially CNS+ remain, in line with the theory that centrally acting drugs are the most dangerous. There are 25 such CNS+ compounds, but the introduced LLL H, LactMed, and PB indicators further reduced this number (n = 21).
The significant correlation coefficient (inverse relationship) between the bioavailability of the drug in breast milk M/P with the value of PBcode and PB (R = 0.52 and R = −0.42 respectively) indicates that the high level of protein binding reduces the risk of drug penetration into breast milk. The danger of treatment during the lactation period is associated with the nature of biologically active compounds: acids (a; code: −1) are the safest. They strongly bind to human plasma albumin. The bases (b; code: 1) are the least secure. Here, the ratio of milk content (M) to plasma content (P) will be the greatest. The correlation matrix also shows that the relationship between the acid/base, B2 and M/P is very similar (R = 0.54). Such a result confirms the relationship of drug binding power with proteins and their bioavailability in compartments isolated by the blood–brain and blood–milk barriers.
The main task was to study the interrelationships and relationships between the bioavailability of drugs in the CNS and the dangers of their use in women during breastfeeding. Therefore, the MLR analysis was performed for the dependent variable B2. A summary of the regression obtained is presented below. The test result is very good, not only because of the high correlation coefficient R = 0.91, which explains 83% of the total variability of B2 (n = 33), but also due to the introduction of the M/P parameter to the model. The relationship between B2 and M/P is, as expected, directly proportional (Figure 1, Equation (1)).
B2 = −0.33(±0.19) −0.25(±0.02) HA + 0.15(±0.04) Sa + 0.17(±0.08)M/P
  • R = 0.91, R2 = 0.83, F(3,29) = 47.05, p < 0.00000, s = 0.2811; n = 33
  • Q2LOO = 0.78, SDEP = 0.2839, PRESS = 2.893, SPRESS = 0.2875, Q2LMO = 0.80
The B2 index is most strongly represented by the number of hydrogen bond acceptors (HA) and the molecular surface area of the drug (Sa). Interesting results were also obtained from the MLR model of the dependent variable PB (Equation (2)).
PB = 0.91(±0.41) + 0.014(±0.031)logD − 0.15(±0.066)M/P + 0.14(±0.050)Sa − 0.073(±0.027) HA + 0.058(±0.040)eH-eL
R = 0.76, R2 = 0.57, F(5,27) = 7.244, p < 0.00020, s = 0.2417
The repetition of all MDs associated with the B2 value (HA, M/P, Sa) in the resulting PB model suggests the inclusion of the chromatographic data connected with protein binding in explaining the penetration of APIs into the CNS and breast milk. Due to this connection, the chromatographic systems describing the protein binding capacity became analytical models for two pharmacokinetic phenomena—the distribution to the CNS and milk.
Further analyses were based on the chromatographic data: Rf values from TLC (NP and RP mode) and HPLC retention times were collected using two columns: HSA and IAM. Initially, a study of the results of using chromatographic data for a set of CNS +/− cases (n = 37) was performed and presented in the following correlation matrices (Table 4), then all types of chromatographic systems used were checked for the group of bioavailable CNS+ (n = 25) (Table 5).
Both groups of cases (CNS+/− and CNS+ only) clearly associate bioavailability in breast milk with the behavior of compounds tested in a thin layer chromatography environment. They apply both to the stationary phase, obtained by modifying silica gel plates (NP), and to silica gel plates, which are silanized (RP). Slight differences between the groups of cases—CNS+/−, CNS+ result from the size of the groups. The results of experiments using HPLC yielded worse results.
An additional modifier was the Rf value of drugs on non–protein coated plates (control plate—C). The introduction of the Rf/C parameter—NP/C and RP/C—allows the observation of chromatographic effects related only to the presence of the protein in the stationary phase. The modified data are presented separately in the following correlation matrices from various chromatographic experiments (Table 6).
The observation of the correlation of chromatographic data by type of chromatography indicates RP TLC as the system that best describes the characteristics of drug excretion into breast milk in the group of CNS-active APIs. The RP/C descriptor also has the highest correlation with M/P (R = −0.57).

3.2. Random Forest Regression

RF regression analyses were performed using seven MDs. Four models were created, each time using different chromatographic data (Rf values from NP TLC and RP TLC or log k obtained from HPLC-HSA, HPLC-IAM experiments). The dependent variable was the B2 parameter, similar to the MLR analysis. The molecular descriptors that gave the best results and therefore were used in the models are listed in the Table 7:
The predictor importance plot revealed that TLC data (MD no. 7 on the plots), with an indication of the reversed phase system, had the greatest contribution to the creation of the RF model. This impact is comparable to the contribution of log P or HD. HPLC-HSA log k values show a significantly lesser effect on regression, and HPLC-IAM data indicates a negligible influence. (Figure 2).
The most important MD in RF models is the number of hydrogen bond acceptors (HA). Only PB reveals an inversely proportional relationship, which is the expected conclusion. The M/P parameter shows a similar contribution as physicochemical MDs such as PhCharge or log P. This information creates another link between crossing the BBB and the blood-milk barrier.
RF regression models were ultimately built with only TLC data, as reported by the predictor importance plots. R2 and Q2 values oscillated around 0.85–0.87 and 0.78–0.80, respectively (Figure 3 and Figure 4), which is a satisfactory result. The best fit of the model can be seen in the B2 value range from −1.5 to 0; extreme values are less correlated.
With each model calculation in the RF method, different decision trees are involved. Therefore, the RF regression was repeated 20 times to obtain averaged results. Table 8 shows the mean values of R2 and Q2 from the models. The data from NP TLC gives slightly higher results of the model and cross-validation determination coefficients, but it is not a significant difference.

3.3. Cluster Analysis

This analysis was performed to compare the variability of the blood–brain barrier permeation parameters and the chromatographic data describing them.
The parameter CNS+/− was used as an indicator of the bioavailability of drugs in the brain. As an indicator of pharmacotherapy safety in the feeding period, the M/P code parameter was used, which was created by scaling the M/P parameter. The M/Pcode parameter ranges from 1–4. The values 1, 2, 3 and 4 have cases where the determined M/P level does not exceed 0.3, 0.8, 1.1 and 4, respectively. The logical conclusion also remains that drugs with a lower concentration in milk than in blood can be considered safer (M/P−). If the M/P value exceeds 1, the use of pharmacotherapy in a nursing woman is inadvisable or prohibited (M/P+). The study group was observed with the possibility of dividing it into two or four clusters: drugs CNS+ and CNS− and M/Pcode 1–4, according to the level of safety of use.
The application of two clusters (Appendix A Figure A1) is not effective. The CNS+/− takes the value of 0.5 and 0.8 for clusters with values 0 and 1, and all M/P values, in both clusters that are below 1 (instead of <1 and >1); this means that both clusters gather mixed cases. Such analysis is not effective. This is due to the considerable complexity of both phenomena, which are not obviously defined, even though biomedical models.
The splitting of the cases into four clusters improved the result (Figure 5). In the case of M/P values, the four clusters split into two groups—clusters 2 and 4 have an M/P− < 1, and clusters 1 and 3 are M/P+ > 1. At the same time, clusters 1 and 3 showed M/P+ (3.4–3.8) values, and clusters 2 and 4 corresponded to M/P− (1.1–2.2). The definition of the CNS+/− value strictly follows this division. An interesting observation is the distribution of the values of the chromatographic data. In both TLC and HPLC formats, all extreme data (the smallest and largest Rf and log k values) correspond to clusters 1 and 3. The intermediate clusters are 2 and 4.
Assuming that the inaccuracy of the division of cases into M/P+ and M/P− may be caused by a particularly difficult and ambiguous differentiation of CNS+/− APIs, a grouping analysis of k–means excluding CNS− drugs was performed.
An M/P value of less than 1.0 shows that only minimal amounts of drugs are transferred into the milk; these types of drugs are classified as low risk (LR). Drugs with an M/P value of 1.0 or more may be present in breast milk at a higher concentration than in the mother’s plasma and are classified as high risk (HR). Thus, this leads to more drug being transported into the infant’s body and, consequently, to side effects [12].
As the results of these experiments depend on various uncontrolled variables such as laboratory conditions, geographic region and lactation time, the reliability of these experimental data is questionable. In addition, there are many drugs for which M/P ratios have not been determined, so it is necessary and useful to develop some theoretical methods such as the Quantitative Structure Property Relationship (QSPR) for predicting M/P values. In QSPR approaches, the chemical structures of compounds are quantitatively correlated with their biological activity [28].
In the current analysis, an additional parameter, HR/LR, was introduced, where HR—high risk (code value 1) of passing the drug into milk. The second value is LR–low risk (code with a value of 0). The result of the analysis, with the application of the 4 CNS+ APIs clusters, are shown in the graph below (Figure 6).

4. Discussion

The inversely proportional relationship between the strength of API and protein binding and its potential presence in the CNS has been demonstrated in the presented models. There is also a strong similarity between the acid/base and PB descriptors, and M/P milk penetration (correlation 0.54 and −0.42) or B2 brain penetration (correlation 0.54 and −0.31). This confirms the relationship between API-protein binding and their bioavailability in compartments isolated by blood-brain and blood-milk barriers.
The study reveals that the chromatographic systems can be analytical models for two different pharmacokinetic phenomena—distribution into the CNS and breast milk. The idea of analytical models for imaging pharmacokinetic properties can be based in both cases on the observation of the affinity of APIs to plasma proteins. An analysis of the chromatographic data with the division into the types of chromatography indicates TLC as the system that better describes the characteristics of drug penetration into breast milk when compared with liquid chromatography. The RP TLC Rf value was the most useful in an RF regression model with B2 as the BBB penetration parameter; RP/C has the highest correlation (−0.57) with the M/P descriptor.
M/P was involved in both the MLR and RF regression models of BBB penetration, and the relationship was directly proportional. This may lead to the assumption that the use of CNS+ drugs in breastfeeding women is more dangerous. However, due to difficulties in experimentally determining the M/P values, there is uncertainty in the reported M/P values which may influence the results of QSPR modelling. Therefore, classifying drugs into the high-risk (HR) and low-risk (LR) classes is more appropriate than accurately predicting their M/P values.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/cells12030421/s1, Table S1. Biological properties of 37 APIs. Table S2. Physicochemical properties of 37 APIs. Table S3. Chromatographic data and their derivatives obtained from NP TLC and RP TLC experiments. Table S4. Chromatographic data and their derivatives obtained from the HPLCHSA i HPLCIAM experiments. Table S5. Structures of 37 APIs.

Author Contributions

Conceptualization, B.E. and W.K.; methodology, B.E.; software, B.E. and W.K.; validation, B.E. and W.K.; formal analysis, B.E. and W.K.; investigation, R.A.; resources, R.A.; data curation, B.E. and W.K.; writing—original draft preparation, W.K.; writing—review and editing, B.E.; visualization, W.K.; supervision, B.E.; project administration, B.E. and W.K.; funding acquisition, B.E. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by an internal grant from the Medical University of Lodz number 503/3-016-03/503-31-001.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Appendix A

Figure A1. Mean values graph of the chromatographic data for two clusters.
Figure A1. Mean values graph of the chromatographic data for two clusters.
Cells 12 00421 g0a1

Appendix B

Chromatographic Experiments

Drugs were subjected to an affinity chromatography analysis for human and bovine plasma proteins in column (HPLC) and thin-layer (TLC) formats, respectively. A total of 37 active substances were isolated with methanol extraction from pharmaceutical preparations. The purity of the isolated substances was checked by TLC chromatography and densitometric scanning. The obtained APIs were then dissolved in 99.8% methanol to give 1 mg/mL solutions which were then used in TLC and HPLC chromatography.
For TLC chromatography, 20 × 20 cm glass plates from Merck that were covered with silica gel and a fluorescent indicator were used. Normal phase (NP) plates were used with regular Merck TLC Silica gel 60 F254 plates, while in the reverse phase (RP), silanized plates RP-2: Merck TLC Silica gel 60 RP-2 F254, silanized. The mobile phase was buffered with ammonium acetate to pH 7.4 (LACH-NER, Ammonium acetate). Composition of the mobile phase: - water, methanol and acetonitrile, with HPLC gradient class (60:20:20 (v/v/v)).
The stationary phase of the plates, both NP and RP-2, was modified with an aqueous solution of bovine serum albumin BSA (Sigma Aldrich bovine serum albumin, lyophilized powder). The unfolded protein-impregnated plates and control plates were scanned with a Desaga CD 60 densitometer. The Rf values (described in the article as NP and RP, respectively) were determined. Analytical wavelengths were selected individually for each API by means of multi-wavelength scanning (200 to 300 nm).
In the HPLC format, a chromatography column immobilized on the resin with human serum albumin HSA (Daicel: CHIRALPAK®HSA, 5 µm; 4 × 10 mm) and IAM artificial biological membrane columns (Regis Technologies Inc., Morton Grove, IL, USA: IAM.PC.DD.2, 10 µm; 4.6 × 10 mm) was used.
Organic solvents (acetonitrile and methanol) and water (J.T. Baker HPLC) were used in HPLC. Ammonium acetate (LACH-NER) was used to prepare the acetate buffer (HPLCHSA), while the phosphate buffer (HPLCIAM) was prepared using (Sigma, Phosphate buffered saline, tablets). The mobile phase for HPLCHSA was a mixture of 10 mM acetate buffer pH 7.4, acetonitrile and methanol in the ratio 85:10:5 (v/v/v). The mobile phase in the case of HPLCIAM was a mixture of 10 mM phosphate buffer pH 7.4 and acetonitrile in the ratio 80:20 (v/v). The parameter that was collected from the HPLC was the logarithm of the retention factor (log k), described as log k(HSA) or log k(IAM).

References

  1. Meng, F.; Xi, Y.; Huang, J.; Ayers, P.W. A curated diverse molecular database of blood-brain barrier permeability with chemical descriptors. Sci. Data 2021, 8, 289. [Google Scholar] [CrossRef] [PubMed]
  2. Dichiara, M.; Amata, B.; Turnaturi, R.; Marrazzo, A.; Amata, E. Tuning Properties for Blood-Brain Barrier Permeation: A Statistics-Based Analysis. ACS Chem. Neurosci. 2020, 11, 34–44. [Google Scholar] [CrossRef] [PubMed]
  3. Subramanian, G.; Kitchen, D.B. Computational models to predict blood-brain barrier permeation and CNS activity. J. Comput. Aided. Mol. Des. 2003, 17, 643–664. [Google Scholar] [CrossRef] [PubMed]
  4. Gleeson, M.P. Generation of a set of simple, interpretable ADMET rules of thumb. J. Med. Chem. 2008, 51, 817–834. [Google Scholar] [CrossRef]
  5. Yoon, C.H.; Kim, S.J.; Shin, B.S.; Lee, K.C.; Yoo, S.D. Rapid screening of blood-brain barrier penetration of drugs using the immobilized artificial membrane phosphatidylcholine column chromatography. J. Biomol. Screen. 2006, 11, 13–20. [Google Scholar] [CrossRef] [Green Version]
  6. Adenot, M.; Lahana, R. Blood-Brain Barrier Permeation Models: Discriminating between Potential CNS and Non-CNS Drugs Including P-Glycoprotein Substrates. J. Chem. Inf. Comput. Sci. 2004, 44, 239–248. [Google Scholar] [CrossRef] [PubMed]
  7. Van De Waterbeemd, H.; Camenisch, G.; Folkers, G.; Chretien, J.R.; Raevsky, O.A. Estimation of blood-brain barrier crossing of drugs using molecular size and shape, and H-bonding descriptors. J. Drug Target. 1998, 6, 151–165. [Google Scholar] [CrossRef]
  8. Zhao, Y.H.; Abraham, M.H.; Ibrahim, A.; Fish, P.V.; Cole, S.; Lewis, M.L.; De Groot, M.J.; Reynolds, D.P. Predicting penetration across the blood-brain barrier from simple descriptors and fragmentation schemes. J. Chem. Inf. Model. 2007, 47, 170–175. [Google Scholar] [CrossRef]
  9. Cruciani, G.; Pastor, M.; Guba, W. VolSurf: A new tool for the pharmacokinetic optimization of lead compounds. In European Journal of Pharmaceutical Sciences; Elsevier: Amsterdam, The Netherlands, 2000; Volume 11. [Google Scholar]
  10. Nauwelaerts, N.; Deferm, N.; Smits, A.; Bernardini, C.; Lammens, B.; Gandia, P.; Panchaud, A.; Nordeng, H.; Bacci, M.L.; Forni, M.; et al. A comprehensive review on non-clinical methods to study transfer of medication into breast milk—A contribution from the ConcePTION project. Biomed. Pharmacother. 2021, 136, 111038. [Google Scholar] [CrossRef]
  11. Kadry, H.; Noorani, B.; Cucullo, L. A blood–brain barrier overview on structure, function, impairment, and biomarkers of integrity. Fluids Barriers CNS 2020, 17, 69. [Google Scholar] [CrossRef]
  12. Gentile, S.; Rossi, A.; Bellantuono, C. SSRIs during breastfeeding: Spotlight on milk-to-plasma ratio. Arch. Womens. Ment. Health 2007, 10, 39–51. [Google Scholar] [CrossRef] [PubMed]
  13. Davanzo, R.; Dal Bo, S.; Bua, J.; Copertino, M.; Zanelli, E.; Matarazzo, L. Antiepileptic drugs and breastfeeding. Ital. J. Pediatr. 2013, 39, 50. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Wanat, K.; Khakimov, B.; Brzezińska, E. Comparison of statistical methods for predicting penetration capacity of drugs into human breast milk using physicochemical, pharmacokinetic and chromatographic descriptors. SAR QSAR Environ. Res. 2020, 31, 457–475. [Google Scholar] [CrossRef]
  15. Maier, R.; Fries, M.R.; Buchholz, C.; Zhang, F.; Schreiber, F. Human versus Bovine Serum Albumin: A Subtle Difference in Hydrophobicity Leads to Large Differences in Bulk and Interface Behavior. Cryst. Growth Des. 2021, 21, 5451–5459. [Google Scholar] [CrossRef]
  16. Luco, J.M.; Salinas, A.P.; Torriero, A.A.J.; Vázquez, R.N.; Raba, J.; Marchevsky, E. Immobilized Artificial Membrane Chromatography: Quantitative Structure-Retention Relationships of Structurally Diverse Drugs. J. Chem. Inf. Comput. Sci. 2003, 43, 2129–2136. [Google Scholar] [CrossRef] [PubMed]
  17. Drugbank. Available online: https://www.drugbank.ca/drugs/DB01174 (accessed on 9 September 2022).
  18. Chembl Database. Available online: https://www.ebi.ac.uk/chembl/ (accessed on 20 September 2022).
  19. Clark, D.E. Rapid calculation of polar molecular surface area and its application to the prediction of transport phenomena. 2. Prediction of blood-brain barrier penetration. J. Pharm. Sci. 1999, 88, 815–821. [Google Scholar] [CrossRef]
  20. Drugs and Lactation Database (LactMed). Available online: https://toxnet.nlm.nih.gov/pda/lactmed.htm (accessed on 11 November 2019).
  21. Hale, T.W. Medications and Mother’s Milk, 15th ed.; Pharmasoft Medical Publishing: Amarillo, TX, USA, 2012. [Google Scholar]
  22. Agatonovic-Kustrin, S.; Tucker, I.G.; Zecevic, M.; Zivanovic, L.J. Prediction of drug transfer into human milk from theoretically derived descriptors. Anal. Chim. Acta 2000, 418, 181–195. [Google Scholar] [CrossRef]
  23. Katritzky, A.R.; Dobchev, D.A.; Hür, E.; Fara, D.C.; Karelson, M. QSAR treatment of drugs transfer into human breast milk. Bioorganic Med. Chem. 2005, 13, 1623–1632. [Google Scholar] [CrossRef]
  24. Li, H.; Yap, C.W.; Ung, C.Y.; Xue, Y.; Cao, Z.W.; Chen, Y.Z. Effect of selection of molecular descriptors on the prediction of blood-brain barrier penetrating and nonpenetrating agents by statistical learning methods. J. Chem. Inf. Model. 2005, 45, 1376–1384. [Google Scholar] [CrossRef]
  25. Fan, Y.; Unwalla, R.; Denny, R.A.; Di, L.; Kerns, E.H.; Diller, D.J.; Humblet, C. Insights for predicting blood-brain barrier penetration of CNS targeted molecules using QSPR approaches. J. Chem. Inf. Model. 2010, 50, 1123–1133. [Google Scholar] [CrossRef]
  26. Lobell, M.; Molnár, L.; Keserü, G.M. Recent advances in the prediction of blood-brain partitioning from molecular structure. J. Pharm. Sci. 2003, 92, 360–370. [Google Scholar] [CrossRef] [PubMed]
  27. Sobańska, A.W.; Wanat, K.; Brzezińska, E. Prediction of the Blood-Brain Barrier Permeability Using RP-18 Thin Layer Chromatography. Open Chem. 2019, 17, 43–56. [Google Scholar] [CrossRef]
  28. Peterson, Y.K.; Wang, X.S.; Casey, P.J.; Tropsha, A. Discovery of geranylgeranyltransferase-I inhibitors with novel scaffolds by the means of quantitative structure-activity relationship modeling, virtual screening, and experimental validation. J. Med. Chem. 2009, 52, 4210–4220. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Scatter plot of the relationships between the predicted and observed variables B2 in the MLR model.
Figure 1. Scatter plot of the relationships between the predicted and observed variables B2 in the MLR model.
Cells 12 00421 g001
Figure 2. Model predictor importance plots. MDs are listed in accordance with Table 5. MD no. 7 is the value obtained in chromatographic experiments: Rf or log k.
Figure 2. Model predictor importance plots. MDs are listed in accordance with Table 5. MD no. 7 is the value obtained in chromatographic experiments: Rf or log k.
Cells 12 00421 g002
Figure 3. Scatter plot of predicted vs. observed B2 values in RF model, using RP TLC data. R2 = coefficient of determination; Q2 = coefficient of determination for the cross-validated models; RMSECV = root-mean-square error of cross-validation.
Figure 3. Scatter plot of predicted vs. observed B2 values in RF model, using RP TLC data. R2 = coefficient of determination; Q2 = coefficient of determination for the cross-validated models; RMSECV = root-mean-square error of cross-validation.
Cells 12 00421 g003
Figure 4. Scatter plot of predicted vs. observed B2 values in an RF model using NP TLC data. R2 = coefficient of determination; Q2 = coefficient of determination for the cross-validated models; RMSECV = root-mean-square error of cross-validation.
Figure 4. Scatter plot of predicted vs. observed B2 values in an RF model using NP TLC data. R2 = coefficient of determination; Q2 = coefficient of determination for the cross-validated models; RMSECV = root-mean-square error of cross-validation.
Cells 12 00421 g004
Figure 5. Mean value graph of the chromatographic data for four clusters: CNS+, CNS−, M/P+ and M/P−.
Figure 5. Mean value graph of the chromatographic data for four clusters: CNS+, CNS−, M/P+ and M/P−.
Cells 12 00421 g005
Figure 6. Graph of the mean values of the chromatographic data for four clusters of CNS+ drugs.
Figure 6. Graph of the mean values of the chromatographic data for four clusters of CNS+ drugs.
Cells 12 00421 g006
Table 1. Differences and similarities between blood-brain barrier and blood milk barrier.
Table 1. Differences and similarities between blood-brain barrier and blood milk barrier.
Blood-Brain BarrierBlood-Milk Barrier
structuredouble phospholipid membrane with additional tight (TJs) and adherence junctions (AJs) in the basement membrane
  • double phospholipid membrane
  • the first days of lactation: large cellular spaces are formed between lactocytes, which can be easily penetrated by compounds from blood vessels and lymph
functionprotection of the CNS against the penetration of unnecessary and toxic substances and selective supply of necessary ones (glucose, nutrients, oxygen)passing nutrients into milk, and retaining unnecessary substances (mainly macromolecular ones)
additional activity
  • P-glycoprotein activity
  • specific enzymatic activity—breakdown of chemical compounds
  • system of active removal of substances (brain-to-blood efflux), carried out by transporters
not observed
mechanisms of penetration
  • simple diffusion
  • membrane transporters
  • endocytosis (often via receptors)
  • mainly simple diffusion
  • secretory route (beginning of the lactation, depends on the structure of the mammary gland)
Table 2. List of molecular descriptors (MDs) used in statistical analyses.
Table 2. List of molecular descriptors (MDs) used in statistical analyses.
MDDescriptionSource
acid/baseacidic, basic or neutral nature of an APIChembl
B2computational parameter, determines penetration through the blood-brain barrier: log bb = 0.547–0.016 PSA[19]
CNS+/−describes the bioavailability in the CNSDrugbank
eL-eHionisation capacityHyperchem
HAhydrogen bonds acceptorsHyperchem
HDhydrogen bonds donorsHyperchem
LactMeddescribes the bioavailability in the breastmilk (code: 0–2)[20]
LLL HHale’s toxicity scale (code: 1–5)[21]
log Ddistribution coefficientACD Labs
log kHSAlogarithm of retention factor from HPLCHSAHPLC
log kIAMlogarithm of retention factor from HPLCIAMHPLC
log Ppartition coefficientHyperchem
log U/Ddescribes the extent of ionisation, calculated from pKa, according to the equations: pKa—pH (acids) or
pH—pKa (bases)
ACD Labs
M/Pmilk-to-plasma ratio; ratio between API concentration in milk and plasma[21,22,23]
M/PcodeM/P converted to categorical values (code: 1–4)
NP; NP/Cthe Rf from NP BSA-modified plate; Rf from modified plate/Rf from clear one (control)TLC
PBprotein bindingDrugbank
PBcodePB converted to categorical values (code: 1–5)
PhChargecharge of a compound in physiological environmentDrugbank
PSApolar surface area of a moleculeHyperchem
RP; RP/Cthe Rf from RP-2 BSA-modified plate; Rf from modified plate/Rf from clear one (control)TLC
Sasurface area of a moleculeHyperchem
Table 3. Correlation matrix of the most important molecular descriptors for crossing the BBB.
Table 3. Correlation matrix of the most important molecular descriptors for crossing the BBB.
n = 37CNS+/−B2 > −0.9B2PhChargeAcid/BaseB2 > −0.52PB
CNS+/−1.00
B2 > −0.90.361.00
B20.570.771.00
PhCharge0.01−0.15−0.181.00
acid/base0.280.100.310.711.00
B2 > −0.520.530.780.730.000.231.00
PB−0.18−0.41−0.45−0.07−0.15−0.261.00
Table 4. Correlation matrix in the CNS+ group.
Table 4. Correlation matrix in the CNS+ group.
n = 21B2LactMedLLL Hlog U/DM/PPBPBcodeAcid/Base
B21.00
LactMed0.0851.00
LLL H0.140.0991.00
log U/D−0.14−0.25−0.281.00
M/P0.250.160.060−0.181.00
PB−0.31−0.067−0.10−0.061−0.421.00
PBcode−0.170.330.28−0.0450.52−0.821.00
acid/base0.540.300.21−0.0900.540.0320.301.00
Table 5. Internal correlations between chromatographic parameters for CNS+/− cases and CNS+ only.
Table 5. Internal correlations between chromatographic parameters for CNS+/− cases and CNS+ only.
CNS+/−, n = 37
NPRPlogk(HSA)logk(IAM)
M/P−0.30−0.27−0.041−0.031
CNS+, n = 25
NPRPlogk(HSA)logk(IAM)
M/P−0.33−0.42−0.0590.24
Table 6. Internal correlations between chromatographic parameters and control for CNS + APIs.
Table 6. Internal correlations between chromatographic parameters and control for CNS + APIs.
CNS+, n=25
NPNP/CRPRP/C
M/P–0.33–0.35–0.42–0.57
Table 7. Molecular descriptors used in random forest regression.
Table 7. Molecular descriptors used in random forest regression.
No.Molecular DescriptorNo.Molecular Descriptor
1.M/P5.HD
2.PB6.log P
3.PhCharge7.chromatographic data: Rf or log k
4.HA
Table 8. RF modelling results, models were created 20 times, and each time different decision trees were involved.
Table 8. RF modelling results, models were created 20 times, and each time different decision trees were involved.
NP TLCRP TLC
R2 = 0.8615; Q2 = 0.7915R2 = 0.8484; Q2 = 0.7850
CV = 55.725CV = 56.105
RMSECV = 0.3279RMSECV = 0.3302
p = 5.087 × 10−17p = 3.454 × 10−16
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Karolina, W.; Agata, R.; Elżbieta, B. Computational Approach to Drug Penetration across the Blood-Brain and Blood-Milk Barrier Using Chromatographic Descriptors. Cells 2023, 12, 421. https://doi.org/10.3390/cells12030421

AMA Style

Karolina W, Agata R, Elżbieta B. Computational Approach to Drug Penetration across the Blood-Brain and Blood-Milk Barrier Using Chromatographic Descriptors. Cells. 2023; 12(3):421. https://doi.org/10.3390/cells12030421

Chicago/Turabian Style

Karolina, Wanat, Rojek Agata, and Brzezińska Elżbieta. 2023. "Computational Approach to Drug Penetration across the Blood-Brain and Blood-Milk Barrier Using Chromatographic Descriptors" Cells 12, no. 3: 421. https://doi.org/10.3390/cells12030421

APA Style

Karolina, W., Agata, R., & Elżbieta, B. (2023). Computational Approach to Drug Penetration across the Blood-Brain and Blood-Milk Barrier Using Chromatographic Descriptors. Cells, 12(3), 421. https://doi.org/10.3390/cells12030421

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop