Next Article in Journal
The Effect of Microwave Pretreatment on Micronutrient Contents, Oxidative Stability and Flavor Quality of Peanut Oil
Next Article in Special Issue
A New Class of Synthetic Flavonolignan-Like Dimers: Still Few Molecules, but with Attractive Properties
Previous Article in Journal
Optimization of the Steam Explosion Pretreatment Effect on Total Flavonoids Content and Antioxidative Activity of Seabuckthom Pomace by Response Surface Methodology
Previous Article in Special Issue
O-Aminoalkyl-O-Trimethyl-2,3-Dehydrosilybins: Synthesis and In Vitro Effects Towards Prostate Cancer Cells
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dermal Delivery of Selected Polyphenols from Silybum marianum. Theoretical and Experimental Study

by
Pavel Kosina
1,†,
Markéta Paloncýová
2,3,†,
Alena Rajnochová Svobodová
1,
Bohumil Zálešák
4,
David Biedermann
5,
Jitka Ulrichová
1 and
Jitka Vostálová
1,*
1
Department of Medical Chemistry and Biochemistry, Faculty of Medicine and Dentistry, Palacký University, Hněvotínská 3, 775 15 Olomouc, Czech Republic
2
Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, tř. 17 listopadu 12, 771 46 Olomouc, Czech Republic
3
KTH Royal Institute of Technology, School of Engineering Sciences in Chemistry, Biotechnology and Health, Department of Theoretical Chemistry & Biology, S-106 91 Stockholm, Sweden
4
Department of Plastic and Aesthetic Surgery, University Hospital Olomouc, I. P. Pavlova 6, 779 00 Olomouc, Czech Republic
5
Institute of Microbiology, Laboratory of Biotransformation, Academy of Sciences of the Czech Republic, Vídeňská 1083, 142 20 Prague, Czech Republic
*
Author to whom correspondence should be addressed.
The contribution as the first author was equal.
Molecules 2019, 24(1), 61; https://doi.org/10.3390/molecules24010061
Submission received: 6 December 2018 / Revised: 20 December 2018 / Accepted: 21 December 2018 / Published: 24 December 2018
(This article belongs to the Special Issue Silymarin and Derivatives: From Biosynthesis to Health Benefits)

Abstract

:
Silymarin is a well-known standardized extract from the seeds of milk thistle (Silybum marianum L., Asteraceae) with a pleiotropic effect on human health, including skin anticancer potential. Detailed characterization of flavonolignans properties affecting interactions with human skin was of interest. The partition coefficients log Pow of main constitutive flavonolignans, taxifolin and their respective dehydro derivatives were determined by a High Performance Liquid Chromatography (HPLC) method and by mathematical (in silico) approaches in n-octanol/water and model lipid membranes. These parameters were compared with human skin intake ex vivo. The experimental log Pow values for individual diastereomers were estimated for the first time. The replacement of n-octanol with model lipid membranes in the theoretical lipophilicity estimation improved the prediction strength. During transdermal transport, all the studied compounds permeated the human skin ex vivo; none of them reached the acceptor liquid. Both experimental/theoretical tools allowed the studied polyphenols to be divided into two groups: low (taxifolin, silychristin, silydianin) vs. high (silybin, dehydrosilybin, isosilybin) lipophilicity and skin intake. In silico predictions can be usefully applied for estimating general lipophilicity trends, such as skin penetration or accumulation predictions. However, the theoretical models cannot yet provide the dermal delivery differences of compounds with very similar physico-chemical properties; e.g., between diastereomers.

Graphical Abstract

1. Introduction

The primary function of the skin is to form a barrier between the body and the environment that is essential for maintaining body homeostasis. The epidermal barrier prevents excessive trans-epidermal loss of water, electrolytes or heat, but also serves as a barrier to the entry of harmful environmental toxic substances and infectious microorganisms or protects the body against damage caused by physical stimuli such as heat or radiation. The outermost layer of the skin, the stratum corneum, although very thin provides the main barrier to the transport of most topically applied compounds. The stratum corneum is formed from denucleated, keratin-filed, non-living cells called corneocytes that are anchored in a lipid matrix (major components include ceramides, fatty acids, and cholesterol). With its specific composition and morphology, the stratum corneum forms a unique biomembrane that is about a thousand times less permeable to water and compounds than other membranes in the body [1,2]. The cutaneous delivery of substances represents an attractive option for the treatment or prevention of various dermatological conditions with minimal side-effects and avoids metabolic biotransformation that may occur in the liver. However, the effective barrier function of the skin limits the use of many compounds as transdermal agents [2]. Therefore, a promising candidate for dermatological application has to cross the stratum corneum barrier and reach viable cells in the deeper layers of the epidermis and dermis.
Silymarin (SM) is a multicomponent extract from the seeds of the milk thistle (Silybum marianum L., Asteraceae). S. marianum is one of the oldest known herbal plants, and has been widely used in traditional European medicine for over two thousand years, especially for treating liver disorders. The main polyphenolic component of SM is the flavonolignan silybin (SB) [3]. Other components in considerable amounts include isosilybin (ISB), silychristin (SC), silydianin (SD) and their flavonoid precursor taxifolin (TA). Dehydrosilybin (DSB) is also present in small amounts [4]. The studied compounds’ structures, including available diastereomers, are presented in Figure 1.
Besides the above-mentioned congeners of silybin, a lot of minor flavonolignans have been identified to date [5,6,7]. Multiple biological actions including antioxidant, anti-inflammatory, immunomodulatory or cytoprotective, make SM and SB attractive for dermatological applications. Therefore, the beneficial effects of SM and SB have been studied on skin carcinogenesis induced chemically [8] or by chronic exposure to UVB light [9], on skin disorders such as contact dermatitis [10], skin hyperpigmentation (melasma) [11], and wound healing [12]. Although many in vivo or clinical studies with SM/SB topical application were conducted, there is minimal information on the skin absorption and permeability (bioavailability) of SB and other SM components.
Efforts to characterize the interactions of compounds with biological membranes have led to an investigation of the hydrophobicity of the utilized molecules [13,14]. Apart from experimental methods, a significant amount of work has been done on in silico predictions in recent years. Hydrophobicity in terms of the widely used n-octanol/water partition coefficient log Pow can be predicted by structure-based tools, such as the freely available AlogP [15], Molinspiration [16] etc. Most of these methods use topological descriptors, or the fragmentation of molecules to functional groups or atoms. From fitting a large number of available experimental data, the effects of the molecular fragments present in molecules on their log Pow can be estimated. These methods are strictly focused on log Pow and work well, especially among molecules similar to the training set. Quantum-calculation based methods or molecular dynamics simulations are another possibility for in silico predictions. These methods are more general and usually not focused on the prediction of a single property. Their performance in terms of log Pow can be worse, but they are not limited to a single case or a certain molecular type. Currently available tools for log Pow prediction have been summarized in a review by Tetko and Poda [17].
The aim of this study was (i) to determine the partition coefficients log Pow of selected polyphenols of SM (SB, ISB, SD, DSB, SC, TA) as well as dehydroderivatives dehydrosilychristin (DSC), dehydrosilydianin (DSD), dehydroisosilybin (DISB), and quercetin (QU) were included in the study (see Figure 1) by High Performance Liquid Chromatography (HPLC) and by mathematical (in silico) approaches in an n-octanol/water system and to compare the obtained values; (ii) to calculate the permeability coefficients of selected SM components for model membranes (bilayer of dioleoylphosphatidylcholine, ceramide or mixture ceramide, cholesterol and lignoceric acid) and to compare them with experimental data obtained using human skin fixed in Franz diffusion cells.

2. Results and Discussions

2.1. Experimental Partition Coefficients log Pow,exp

The partition coefficients (log Pow,exp) of studied polyphenols including available diastereomers were obtained using reverse phase HPLC in the pH range of 4.5–8.5. For evaluating the studied compounds’ hydrophobicity, log Pow,exp of all reference substances were used except for aniline, whose log Pow,exp value (0.9) was conspicuously lower than that of the other reference substances (their individual log Pow and applied amounts are presented in Appendix Table A1). Aniline was only used to evaluate the log Pow,exp of TA at pH 8.5 due to a big negative shift in the TA log Pow,exp value.
The relation between the capacity factor k and log Pow of reference compounds at various pH of mobile phase, employing quadratic equations, was evaluated. The recommended application of the linear equation [14] based on reference compound capacity factors (retention times) and log Pow values for evaluating the test compounds’ log Pow gave only an r2 value of 0.89, see Table A2. The application of the quadratic equation resulted in obtaining more valuable data for the studied compounds’ log Pow,exp evaluation (r2 ~ 0.99), see Table A2. Therefore, the quadratic equations were employed to estimate the log Pow,exp values of the studied compounds.
A minimal influence of temperature (25 °C ~ recommended by OECD107 [13]) vs. 32 °C ~ corresponding to the skin surface temperature) on log Pow,exp values was found (data not shown). Similarly, a minor effect of temperature was found for theoretical (calculation) models (data not shown). The log Pow,exp values of the studied polyphenols, including data for diastereomers, were evaluated in 5 mM phosphate (Table 1), additionally in 5 mM formate and 5 mM acetate buffer as well, in a gradient with methanol at 25 °C with the pH range of 4.5–8.5. A minimal effect of mobile phase composition (phosphate vs. formate vs. acetate) was found (Table 1, Table A3 and Table A4). Application of acetate or formate buffers instead of phosphate allows mass spectrometry identification of studied compounds.
Of all the studied compounds, the lowest log Pow,exp values were found for TA and the highest ones for DISB at all pH values. Among flavonolignans and their dehydro derivatives, SC A exhibited the lowest log Pow,exp value. The hydrophobicity of flavonolignans in increasing order was as follows: SC A < SD < SC B < DSD < SB A < SB B < DSC < ISB A < ISB B < DSB < DISB. Flavonolignan dehydro derivatives were more hydrophobic than the corresponding parent substances. The same effect was observed with the pair of flavonoids TA and QU. The increasing pH value of the mobile phase decreased the log Pow,exp value for all polyphenols (Table 1).

2.2. Calculated Partition Coefficients log Pow,calc

The computational tools for lipophilicity evaluation (log Pow,calc), AlogP [15], Molinspiration [16] and COSMOtherm 15 [18] showed a slightly different order of the molecular hydrophobicity of the studied compounds compared to the experimental evaluation, but the general trends remained the same (Table 2).
TA was the least lipophilic molecule and DISB was the most hydrophobic. Also, the dehydro derivatives (DISB, DSB, DSC and QU) were more lipophilic than their respective parent substances, which agrees with experimental data. The AlogP tool predicted log Pow,calc values for the studied compounds within the range of 1.1–2.8 and the Molinspiration tool computed log Pow,calc values within the range of 0.7–2.4. Using COSMOtherm 15, a wider range of log Pow,calc values were obtained (1.0–4.2) than with the above-mentioned tools (Table 2). The correlation of the experimental (log Pow,exp) and calculated (log Pow,calc) data in terms of linear fit, correlation coefficient (r2) and mean differences (MD) was evaluated (for data see Appendix Figure A1 and Table A5).
The lowest r2 values were found with COSMOtherm 15 and the highest with Molinspiration (pH 4.5–7.5) and AlogP (pH 8.5). The lowest MD were obtained with AlogP over the whole pH range (Table A5).

2.3. Lipid/Water Partition Coefficients log Plw,calc

Apart from n-octanol/water partition (log Pow,calc), lipid/water partition coefficients (log Plw,calc) in membrane models (dioleoylphosphatidylcholine (DOPC), ceramide (CER) and ceramide, lignoceric acid and cholesterol mixture (CLC)) were calculated using COSMOmic (Table 2). The best correlation between log Plw,calc and log Pow,exp was obtained using the DOPC membrane model, where r2 reached up to 0.84 with respect to the data at pH 4.5 and 5.5. The correlation of log Plw,calc values for CER and CLC membrane models with experimental data was lower than for the DOPC membrane model, as the r2 values varied from 0.75 to 0.82. In terms of MD, the CER membrane model performed best (MD in the range of 0.24 to 0.37 over the whole pH range). In the DOPC membrane model, the MD increased with pH, while in the CLC membrane model the MD decreased with pH. Although working on a very different principle, the CER membrane model appeared to work with a similar level of accuracy to AlogP in terms of MD and together with the DOPC membrane model scored very well in terms of r2 (for data see Appendix Figure A1 and Table A5).

2.4. Skin Permeability of Polyphenols

The skin permeability of polyphenols was studied using Franz diffusion cells with human skin. None of the studied compounds were found in the acceptor liquid, however all were found in the tissue. The skin intake of the studied compounds during transdermal transport was influenced by the ethanol content and pH of the donor liquid. The results are presented in Table 3 and expressed as the amount of individual polyphenols in pmol per g of skin tissue per 24 h.
Apart from DSB, the skin content of all polyphenols was lower when applied in donor liquid with pH 8.5 than pH 6.5. When donor liquid with pH 6.5 was used, the skin intake of polyphenols was in the order: TA < SD < SC < DSB < SB < QU < ISB and at pH 8.5 it was TA < SD < SC < QU < SB < ISB < DSB. For the main silymarin constituent (SB), both diastereomers SB A and SB B were determined. The amount of SB B in skin was found to be slightly higher than that of SB A (data not shown). The skin content of QU dramatically decreased (~8.5 times), when QU was applied in donor liquid with pH 8.5.
Increasing ethanol content had a mostly negative effect on the polyphenol penetration into the skin (Table 3). The skin bioavailability of a substance is usually predicted by its log Pow value [2]. Therefore the log Pow,exp values were compared with the skin intake of polyphenols. At pH 6.5, the correlation coefficient r2 of log Pow,exp and logarithm of skin intake varied between 0.54 to 0.60 (see Appendix BTable A6). However, inspecting the data (Figure 2), DSB was identified as an outlier; after its exclusion the r2 between log Pow,exp, at pH 6.5 and the logarithm of the measured skin intake reached 0.94 (Table A6). The log Pow,exp values at pH 8.5 were not as suitable for estimating the skin intake as the data obtained at pH 6.5. Generally, the highest correlation was observed for ethanol-free systems.

2.5. Calculated Permeability of Polyphenols log Perm

From the free energy profiles on all membrane models (CER, CLC and DOPC, see Figure 3) obtained by COSMOmic, the permeability (log Perm) of polyphenols were calculated (Table 4). The order of free energy profiles of studied polyphenols (Figure 3) as well as a penetration barrier (ΔGpen) and Perm (Table 4) corresponds to lipophilicity (log Plw,calc), calculated by theoretical tool (COSMOmic), see Table 2. The lowest free energy profiles were found for DISB and DSB; and the highest one for TA.

2.6. Calculated Partition Coefficient log Plw,calc and Permeability log Perm Correlation with Polyphenols’ Skin Intake

The experimental skin intake was compared to log Plw,calc and log Perm (for data see Appendix BTable A7). For the skin intake at pH 6.5 (no ethanol) after excluding QU, a better correlation was found with log Perm in CER and CLC membrane models (r2 = 0.95 and 0.94). The DOPC membrane model did not distinguish well between more hydrophobic compounds, the correlation was lower. For the skin intake at pH 8.5, a correlation was found with log Plw,calc for all the membrane models used (Figure 4).
Of the computative log Pow tools used, the best prediction of log Plw,calc, was obtained with the DOPC membrane model (r2 = 0.97). The ethanol content and/or pH increase in the donor liquid influenced the correlation between the skin intake (experimental data) and log Perm or log Plw,calc. For both pH values and ethanol content (5 and 10% (v/v)) a decrease in r2 was observed. At pH 6.5 and an ethanol concentration of 15% (v/v) an increase in r2 was found for all membrane models (for data see Appendix B Table A7). A clear linear dependence was found between log Perm and the skin intake for the less hydrophobic compounds (TA, SD, SC). The calculation was unable to clearly distinguish between the highly hydrophobic compounds, as Perm depended exponentially on the highest free energy. The hydrophobic compounds with high log Plw,calc (affected primarily by the lowest free energy well—see Material and Methods section: Design of computational of hydrophobicity of polyphenols)—did not show any significant dependence of the skin intake on log Perm at the DOPC membrane. For other membrane models, where the penetration free energy barrier is higher than with DOPC (Figure 3, Table 4), a dependence of skin intake on log Perm over the whole range was found, however with two outliers, DSB and QU.

2.7. Discussions

Despite the use of silymarin in dermatological preparations, a limited amount of information about the hydrophobicity of silymarin’s polyphenols is available. Our study focused on characterizing the hydrophobicity of selected flavonolignans (including diastereomers and dehydro derivatives) and the flavonoids taxifolin and quercetin and their interaction with model membranes and human skin using theoretical and experimental tools.
For compounds with a possible dermal application, it is essential to be able to cross the lipids layer barrier of the stratum corneum and be soluble in both hydrophobic and hydrophilic environments. These ambivalent properties can be evaluated by the n-octanol/water partition coefficient [2]. The official OECD method, based on the partition of a studied compound between the water and n-octanol phases (n-octanol/water partition coefficient; log Pow) measured by the shake flask method, imitates the compound’s interaction with the biomembrane very well [13]. In this way the log Pow also characterizes the partitioning of a compound between the lipophilic stratum corneum and the underlying hydrophilic living cells of the epidermis [2]. Today, instead of the shake flask method, another OECD Guideline, Test No. 117 using isocratic HPLC elution was applied [14], but in gradient elution mode [19]. The methods for log Pow evaluation are based on a comparison of the log Pow values of a set of reference compounds and their retention times (capacity factors) in a C18 chromatographic column with retention times and capacity factors of studied compounds [14]. The application of the HPLC method for estimating the log Pow values of polyphenols has several advantages: a stable temperature during analysis, the application of small amounts of pure studied compounds as well as their mixtures with similar structure and/or similar physico-chemical properties. By the HPLC gradient elution method, the mixture of main constitutive S. marianum flavonolignans (SB, SC, SD, ISB), including diastereomers for SB, SC and ISB, and their respective dehydro derivatives (eleven individual compounds) as well as TA and its dehydro derivative QU could be separated in one run and the obtained data were suitable for the evaluation of log Pow,exp values (Table 1).
Data for individual diastereomers of SB, SC and ISB as well as for DISB are published here for the first time. The log Pow,calc values that we found for QU and TA (Table 2) are in agreement with reports by Maroziene et al. [20]. The log Pow,exp value of QU corresponds to the data by Rothwell et al. [21]. Similarly, the log Pow,exp values for the pair of flavonolignans DSB and SB obtained by the shake flask method published by Gažák et al. [22] agree with our data obtained by the HPLC method. Here we further found that the hydrophobicity of SM’s compounds decreased with increasing pH in all applied solutions (mobile phases), which corresponds to the results for TA, SB, ISB, SC, SD obtained by the shake flask method [23]. A small effect of pH on log Pow,exp was observed at pH from 4.5 to 6.5. An obvious decrease in hydrophobicity (log Pow,exp) was found from pH 6.5 to 8.5 (Table 1). This phenomenon is connected with the ionization of ionizable groups in the studied compounds. The pKa1 value of SB, SC, SD [24], TA [25] and QU [26] is close to neutral pH (~7.0) and pKa2 is ~8.5. The strongest decrease in log Pow in the case of TA at pH 8.5 relates to its pKa values (pKa1 7.1; pKa2 8.6; pKa3 8.59; pKa4 11.82) [25] as well. Generally, the neutral form of the molecule is more lipophilic than the deprotonated one. The order of the log Pow,exp of the studied polyphenols (TA, SB, SC, SD and ISB) was not the same as that found by Zeng et al. [23]. It could be caused by the experimental limitations of the shake flask method. Studied compounds can occur in both phases and also in an interphase; this phenomenon is eliminated in the reverse phase chromatography method used in our study.
The experimental evaluation of compound hydrophobicity is time-consuming, experimentally (the shake flask method) and financially (reverse phase HPLC based approaches) demanding and also unfavorable to the environment (does not follow green chemistry approaches) [19]. Currently a large number of sophisticated computational tools, based on various parameters, are available [17,27,28]. These in silico methods are helpful for screening a large set of compounds and for basic behavior classification in terms of compound hydrophobicity and bioavailability, reducing the amount of in vivo experiments or clinical trials [29]. On the other hand, the theoretical methods have some limitations, including the inability to properly simulate all situations in a real biological system (i.e., skin). In our study the tools COSMOtherm 15, AlogP and Molinspiration were used. COSMOtherm 15 is generally focused on the partitioning between fluid phases and not directly on log Pow, and does not consider the pH effect; therefore, it scored worst of the calculation methods used. At pH 4.5–7.5, the highest correlation was obtained with Molinspiration and at pH 8.5 with AlogP. This may be caused by AlogP’s use of electrotopological state descriptors [30]. This is an atomic approach that also considers the valence state of the atoms, which may better reflect the experimental molecule ionization at pH 8.5 (corresponding to the pKa of studied polyphenols). AlogP and Molinspiration tools are based on a large set of experimental log Pow values and both are applicable around neutral pH. None of the theoretical tools used were able to reflect the behavior of the set of polyphenols over the whole pH range (4.5–8.5). According to the correlations obtained, the choice of optimal log Pow calculation tool should consider the desired real experimental conditions ~ the real applicative conditions of the compound. When the studied polyphenols are in ionized form (pH over pKa) only the AlogP tool is able to estimate log Pow,calc values. At pH below pKa, Molinspiration described the system best. The replacement of n-octanol with lipid membrane models (DOPC, CER or CLC) in COSMOmic improved the prediction strength. COSMOtherm and COSMOmic inputs are a 3D structure of the compound (conformational variability is considered) [18] and thus differentiation between diastereomers A and B is possible. The differences in the hydrophobicity of the diastereomers were too small to distinguish them, and this ability remains a benefit of experimental evaluation. ALogP and Molinspiration tools do not work with diastereomerism.
Although n-octanol/water partitioning log Pow is taken as a gold standard for lipophilicity estimation [19], considering the lipophilicity directly in terms of the lipid/water partition coefficient log Plw predicts it better. However, using theoretical approaches to characterize the transdermal delivery in vivo is still complicated due to the complexity of the systems in terms of composition or the real interaction between human skin tissue and a topically applied compound. Therefore, both approaches were applied in our study. Many papers have demonstrated that increasing lipophilicity increases the skin permeation of compounds, and a log Pow of 2–3 seems to be optimal. It is likely that these molecules with intermediate lipophilicity can permeate via both the lipid and polar microenvironments in the intercellular route [31]. Accordingly, most of the studied polyphenols are good candidates for dermal application, especially at acidic pH (4.5–6.5), see Table 1. Our transdermal transport data demonstrated that none of the studied compounds were able to penetrate through the skin, as none of them were detected in the acceptor liquid. However, all the polyphenols were found in the human skin (Table 3). In this way the human skin works as a trap and only the skin intake (accumulation) of polyphenols can be measured. QU skin intake should be taken with care, as the presented results on human skin were affected by QU instability in aqueous solutions that was previously reported [32]. Thus, with QU it is not possible to ensure a constant QU concentration in the donor liquid during experiments (24 h) as with the other studied polyphenols. Due to QU decomposition, the amount in the skin is most likely significantly reduced. A higher instability of QU occurs in alkaline pH [32], which corresponds with a lower amount of QU in the skin at the higher pH of 8.5 in our experiments (Table 3). Therefore, correlations of the polyphenol skin intake with log Perm and log Plw,calc were performed without QU.
Although the long-term practical use of silymarin, flavonolignan skin delivery is poorly documented. The previous work [33] that studied the in vitro delivery of SB, SD and SC using native, chemically, and physically modified mouse skin does not agree with our results (Table 3). In the mouse models, the flavonolignans were able to penetrate through the skin into the acceptor liquid (SB > SD > SC). However, the human skin has a higher barrier function (more robust stratum corneum and thicker epidermis) than mouse skin and so mouse skin is not a fully accepted model to mimic the penetration of compounds through the human skin [34]. Nevertheless, SB deposition in mouse skin was notably higher than those of SC and SD [33], which corresponds to our data (Table 3). The order of polyphenol permeation through CER and CLC membrane models (Table 4) also partially agrees with the published data on mouse skin [33]. According to our data on the experimental skin uptake, the studied polyphenols can be divided into two groups: low intake (TA, SC, SD) and high intake (SB, DSB, ISB). This division also corresponded to their hydrophobicity based on theoretical data (log Pow,calc, log Plw,calc) and experimental log Pow values as well. Our results further showed a significant effect of pH on the skin intake of flavonolignans and flavonoids. All these properties may be important for designing dermal preparations containing these polyphenols.

3. Materials and Methods

3.1. Studied Compounds

Pure flavonolignans (SB, SC, SD and ISB) and their respective oxidized products (dehydrosilychristin (DSC), dehydrosilydianin (DSD), dehydroisosilybin (DISB) and dehydrosilybin (DSB)) were prepared at the Institute of Microbiology, Academy of Sciences of the Czech Republic, Prague, Czech Republic according to previously published studies. For SB and ISB, individual diastereomers were available [35]. SC (natural mixture of diastereomers) and SD were isolated from silymarin by Sephadex LH-20 column chromatography (Sigma Aldrich, Prague, Czech Republic) as described previously [36]. DSB was prepared as described by Gažák et al. [37]. The preparation of other dehydro derivatives was published elsewhere [38]. Flavonoids TA and quercetin (QU) were purchased from Sigma Aldrich (Prague, Czech Republic).

Studied Compound Solutions for Hydrophobicity Evaluation

The stock solutions of the studied compounds (1 g L−1) were prepared in methanol. To minimize the coelution of analytes during HPLC separation, two mixtures of the studied compounds, at the final concentration of 10 mg L−1, were prepared. Mixture 1 contained SD, diastereomer A of silybin (SB A), diastereomer A of isosilybin (ISB A), DSC, DSD and DISB. Mixture 2 contained TA, mixture of diastereomers A and B of silychristin (SC A and B), QU, diastereomer B of silybin (SB B), diastereomer B of isosilybin (ISB B) and DSB.

3.2. Chemicals

Aniline, phenol, 4-chloroaniline, nitrobenzene, benzene, trichloroethylene, toluene, chlorobenzene, and naphthalene, all p.a. grade, were obtained from Lach-Ner (Neratovice, Czech Republic). Buffer components, Dulbecco’s Modified Eagle medium without phenol red, and other chemicals were obtained from Sigma-Aldrich (Prague, Czech Republic). Methanol and acetonitrile, both HPLC gradient grades, were obtained from Merck (Darmstadt, Germany). All solutions were prepared using reverse-osmosis deionized water (Ultrapur, Watrex, Prague, Czech Republic). Nitrogen (99.999%) was obtained from SIAD Czech (Prague, Czech Republic) and helium (99.999%) from Linde Gas (Prague, Czech Republic).

Reference Compounds for Hydrophobicity Evaluation

Reference compounds (aniline, phenol, 4-chloroaniline, nitrobenzene, benzene, trichloroethylene, toluene, chlorobenzene, and naphthalene) for partition coefficient (log Pow) evaluation by chromatographic method were dissolved in methanol and diluted to the appropriate concentration (Table A1). Finally, the mixture of reference compounds was injected into the chromatographic column.

3.3. HPLC System for Hydrophobicity Evaluation

The HPLC chromatographic system Dionex UltiMate 3000 (Dionex Corp., Sunnyvale, CA, USA) consisting of a degasser (SRD-3400, 4 DEGASSER CH), a binary pump (HPG-3400SD), an autosampler (WPS-3000 TSL ANALYTICAL), a column compartment (TCC-3000RS) and a diode array detector (DAD-3000, 190–400 nm) was equipped with a Purospher Star RP-18e, 55 × 2 mm, 3 µm chromatographic column (Merck, Darmstadt, Germany). The log Pow values of reference compounds, flavonolignans, TA and QU were determined in gradient elution with mobile phase A: 5 mM phosphate (or 5 mM formate or 5 mM acetate buffer, respectively) with 5% MeOH (v/v); in gradient with the mobile phase B: 100% (v/v) MeOH (linear gradient elution (% B, v/v): 0 min (0% B), 22 min (65% B), 25 min (65% (v/v) B), 25.1 min (0% (v/v) B), 30 min (0% (v/v) B). The pH of mobile phase A was set to 4.5, 5.5, 6.5, 7.5, and 8.5 for each utilized buffer. The flow rate was 0.3 mL min−1, the injection volume was 10 µL. The temperature of the autosampler as well as the column oven was set to 25 or 32 °C. Detection was carried out at 254 nm. The dead and retention time of the reference compounds as well as studied compounds were measured. The analysis was carried out at least twice.

3.4. Hydrophobicity Evaluation

The hydrophobicity (lipophilicity) of compounds was evaluated as the partition coefficient (Pow), defined as the ratio of the equilibrium concentrations of the dissolved substance in n-octanol and water [14].
log P o w = log   C   n o c t a n o l C   w a t e r
The gradient reverse phase HPLC method [19] using a C18 column was used to estimate the log Pow of each studied compound by comparing its retention time (tr) with the tr of reference substances with known log Pow. The capacity factor (k) for the reference compound is calculated using the tr of reference compound and the dead time of analysis (t0):
k = t R t 0 t 0
A relation between k and log Pow is usually shown in the following equation, which is used for evaluating the log Pow of a studied compound [19]:
log P o w = a + b × log   k
However, for a more accurate relation between k and log Pow, a polynomial equation of degree 2 (quadratic equation) was applied.
k = a × ( log P o w ) 2 + b × log   P o w + c
For evaluating the effect of pH on Pow value, mobile phases A with various pH levels (4.5; 5.5; 6.5; 7.5, and 8.5) were utilized. Two temperatures (25 and 32 °C) were used during the HPLC analysis.

3.5. In Vitro Skin Penetration (Transdermal Transport)

Breast tissue specimens were obtained from healthy women undergoing plastic surgery at the Department of Plastic and Aesthetic Surgery (University Hospital in Olomouc). The use of skin tissue complied with the Ethics Committee of the University Hospital in Olomouc and Faculty of Medicine and Dentistry, Palacký University, Olomouc (date: 6.4.2009, ref. number: 41/09). All patients had given their written informed consent. The skin fragments were transported in phosphate buffered saline (PBS) containing antibiotics (penicillin (500 mg mL−1), streptomycin (500 U mL−1) and amphotericin B (1.25 mg mL−1)). The skin was then washed with PBS three times and used as a membrane in Franz cells. The stratum corneum was oriented into the donor compartment. Liquid Dulbecco’s Modified Eagle medium without phenol red and antibiotics was used as the acceptor. The donor compartment was filled with the studied compounds (2 mL; 50 μM) in phosphate buffer (50 mM, pH = 6.5 or 8.5) with various concentrations of ethanol (0, 5, 10 or 15% (v/v)). The available diffusion area between the donor and acceptor compartment was 1.77 cm2. The donor compartment and sampling port was closed with parafilm to protect those compartments against evaporation during the experiment. The Franz cells were incubated in a water bath to guarantee a temperature of the donor liquid of 32 °C (temperature of on the human skin surface) and the stirring rate of acceptor liquid was 200 rpm. To minimize photodecomposition of test compounds, the incubation was carried out in the dark.
After 24 h, the aliquots of donor and acceptor liquids were collected. The skin was removed from the Franz cell and washed with PBS and cleaned and dried with a cotton mull. The skin tissue was weighed, cut with scissors into small pieces, mixed with an extraction mixture (ACN/MeOH/H2O, 50:40:10, v/v/v) in the ratio 1:3, and homogenized for 3 min at 21,500 rpm using a ULTRA-TURRAX® T 25 basic IKA® with a S 25 N-10 G dispersing element (IKA®-Werke Staufen, Germany). The homogenate was centrifuged (12,000× g, 10 min, 4 °C). The supernatant was then diluted with HPLC mobile phase A in the ratio 1:1 (an 8-fold final dilution). The donor liquid was finally diluted with HPLC mobile phase A in the ratio 1:9. The acceptor liquid was diluted with HPLC mobile phase A in the ratio 1:1. After centrifugation (12,000× g, 10 min, 4 °C), samples were analyzed by HPLC (see below).

3.6. HPLC/HESI-MS of Selected Flavonolignans, Taxifolin and Quercetin in Biological Matrices

The HPLC chromatographic system was the same as for the hydrophobicity study (Section 3.3) with the following modifications: the monolithic chromatographic column Chromolith Performance RP-18e, 100 × 2 mm and monolithic guard column RP-18e, 5 mm × 2 mm (Merck, Darmstadt, Germany) were utilized for the analysis of individual selected flavonolignans/flavonoids in gradient elution with mobile phase A: MeOH:H2O:CH3COOH (37:63:0.5); in a gradient with mobile phase B: 100% MeOH (linear gradient elution (%, v): 0–1 min (0% B), 1–3 min (100% B), 3–5 min (100% B), 5–5.1 min (0% B), 5.1–8.5 min (0% B)). The flow rate was 0.3 mL min−1, the injection volume was 10 µL. Before HPLC/HESI-MS analysis, all samples were kept in dark. The temperature of the autosampler was set to 10 °C, and the column oven was set to 30 °C. The HPLC system was on-line connected to the quadrupole ion-trap MS instrument LCQ Fleet (Thermo Scientific, Waltham, MA, USA) operating in a negative as well as a positive heated electron spray ionization (HESI) mode. The HESI-MS parameters selected were: heater temperature 100 °C, capillary temperature 295 °C; in negative/positive mode: spray voltage 2.9/4.0 kV, and capillary voltage −20/+31 V. Nitrogen was used as the sheath, auxiliary and sweep gas, and helium was used as the collision gas. The sheath, auxiliary and sweep gas flow rates were 35, 5, and 1 (in arbitrary units). The MS spectrum was monitored in the range 150–1000 m/z for both ionization modes. The MS2 fragments were isolated from the following parent ions for the negative mode: 301, 303, 479, and 481 m/z, respectively; for the positive mode: 303, 305, 481, and 483 m/z, respectively. The normalized collision energy was 25. The quantification of selected flavonolignans, TA, and QU was carried out according to a calibration curve in the range of 0.1–100 ng per injection.

3.7. Design of Computational Evaluation of Hydrophobicity of Polyphenols

In order to obtain a computational prediction of the hydrophobicity of flavonolignans, TA, and QU, the n-octanol-water partition coefficients were calculated (log Pow,calc) with the commercial software COSMOtherm 15 [18] based on quantum-chemical calculations and by the freely available online tools AlogP [15] and Molinspiration [16]. AlogP and Molinspiration are based on the analysis of available experimental data; specific conditions such as temperature, pH, etc. cannot be chosen. Data can be interpreted as an estimation of log Pow,calc for mean conditions in which the original thousands of data were obtained. Details about the computational methods used can be found; e.g., in the review by Tetko and Poda [17]. In contrast, in COSMOtherm 15, the parameters for the calculation can be set up. Log Pow,calc as a partition coefficient between two phases—water (unionized form) and wet n-octanol (n-octanol saturated with water; 27.4% of mass is formed by water, rest is n-octanol)—were calculated at 25 and 32 °C. More details about the calculations can be found in Appendix C. The COSMOmic [39] tool was used for lipid/water partition coefficient calculation (log Plw,calc). Free energy profiles along the axis perpendicular to the membrane plane from water to the middle of the membrane (Figure 5) were calculated. This can be used further for calculations of log Plw,calc or permeability (Perm):
P l w , c a l c = 0 n ( e Δ G ( z ) R T ρ ( z ) w a t e r ρ ( n ) w a t e r ) d z × A P L M l i p i d s m u
P e r m = 1 / o u t s i d e z e Δ G ( z ) R T d z
where ΔG(z) is the free energy at depth z, R is the molar gas constant, T is the thermodynamic temperature, ρ(z)water is the water density at depth z and ρ(n)water is the density of bulk water—the integration runs from the middle of the membrane (z = 0) through the membrane separated into layers parallel to the membrane plane up to the nth layer situated in bulk water. Thanks to the multiplying factor (where APL is the area per lipid, Mlipids is the molecular weight of lipids, mu is the atomic mass constant) the partition coefficient log Plw,calc is in the units used in the experimental work; i.e., kg(lipid)/L(water).
A fluid dioleoylphosphatidylcholine (DOPC) membrane model was used for lipids, as well as a skin-related solid membrane model. The DOPC membrane model is well established and the calculated partitioning into the DOPC membrane reproduces the experiment well [40]. In order to obtain structural inputs for the COSMOmic tool, the lipid bilayers were composed not only with DOPC, but also with ceramide NS24 (CER), or an equimolar mixture of CER, lignoceric acid and cholesterol (CLC). From the free energy profiles on the membranes obtained by COSMOmic, a penetration barrier (ΔGpen) and Perm were also extracted. Detail information about calculation models are mentioned in Appendix C.

3.8. Statistical Analysis

Experimental (log Pow,exp) and theoretical (log Pow,calc, log Plw,calc) hydrophobicity data were compared in order to observe correlations and differences. For the correlations, MS Excel was used to create a correlation matrix with the Analysis Toolpak add-on providing a correlation coefficient between all the calculated and measured data. Further, a linear regression was done between the data measured in the phosphate buffer at different pH and the calculated data (also in MS Excel). Apart from correlation coefficients and linear regression equations, the differences between the experimental and calculated values were evaluated as their mean differences MD.
1 / N i N | log P o w ( l w ) , c a l c , i log P o w   , e x p , i |

4. Conclusions

In conclusion, log Pow values for selected silymarin’s flavonolignans, including their respective diastereomers and dehydro derivatives as well as flavonoids taxifolin and quercetin were experimentally determined. For the first time the data for individual diastereomers were determined. The ex vivo permeation of flavonolignans, taxifolin and quercetin through the human skin showed that none of the studied compounds reach the acceptor liquid, and all of them accumulate in the tissue. The computational prediction of log Pow values and penetration through model lipid membranes for flavonolignans, and their respective dehydro derivatives, taxifolin and quercetin were carried out and compared with experimental values. Although we observed a reliable prediction of log Pow values (r2 > 0.8), none of the models worked for all the experimental conditions (at various pH). Theoretical tools were unable to clearly distinguish between stereoisomers. The use of lipid membranes for the theoretical estimation of lipophilicity provides a more reliable lipophilicity prediction than log Pow. Although in silico models cannot yet predict the dermal delivery differences of compounds with very similar structural motifs (for example optical isomers) with the desired level of certainty, they can be usefully applied for estimating general lipophilicity trends and skin accumulation predictions. Therefore, the theoretically received data should be still confirmed experimentally.
Both experimental and theoretical tools enabled the studied polyphenols to be divided into two groups based on their skin intake and lipophilicity (taxifolin, silychristin, silydianin vs. silybin, dehydrosilybin, isosilybin). This behavior should be considered during the development of dermatological preparations containing silymarin or its pure flavonolignans as biologically active components.

Author Contributions

P.K., M.P., A.R.S., B.Z. and J.V.; methodology, P.K., M.P., A.R.S. and J.V.; writing—original draft preparation, M.P.; formal analysis, B.Z. and D.B.; resources, J.V. and J.U.; funding acquisition, J.V.; project administration, D.B. and J.U.; writing—review and editing.

Funding

This work was funded by the Czech Grant Agency grants 15-10897S and P208/12/G016 and by the Institutional Support of Palacký University in Olomouc RVO 61989592.

Acknowledgments

We wish to thank Eva Směšná and Iveta Hatalová (University Hospital in Olomouc) for assistance with skin tissue donor recruitment and Benjamin John Watson-Jones for language correction.

Conflicts of Interest

The authors report no conflicts of interest. The authors alone are responsible for the content and writing of the paper.

Appendix A

Table A1. Reference compounds’ log Pow values and composition of reference compound mixture.
Table A1. Reference compounds’ log Pow values and composition of reference compound mixture.
Reference Compoundlog PowAmount in Mixture
Aniline0.92 µl L−1
Phenol1.510 ng L−1
4-Chloroaniline1.810 ng L−1
Nitrobenzene1.910 ng L−1
Benzene2.130 µl L−1
Trichloroethylene2.450 µl L−1
Toluene2.730 µl L−1
Chlorobenzene2.850 µl L−1
Naphthalene3.65 ng L−1
Table A2. Relation between capacity factor (k) and log Pow of reference compounds in phosphate buffer.
Table A2. Relation between capacity factor (k) and log Pow of reference compounds in phosphate buffer.
pHLinear Equitationr2Quadratic Equitationr2
4.5y = 15.564x − 9.21410.8970y = −7.6607x2 + 54.46x − 55.2690.9953
5.5y = 14.935x − 6.56720.8942y = −7.5071x2 + 53.051x − 51.6980.9964
6.5y = 14.855x − 6.42290.8965y = −7.3626x2 + 52.238x − 50.6850.9962
7.5y = 14.895x − 6.22640.8924y = −7.5201x2 + 53.077x − 51.4350.9954
8.5y = 14.99x − 6.73670.8925y = −7.5849x2 + 53.501x − 52.3350.9959
The separation of reference compounds was carried out in 5 mM phosphate buffer with MeOH (5%, v/v), pH was set in the range of 4.5–8.5 (HPLC mobile phase A) in a linear gradient with MeOH (100%; HPLC mobile phase B) at 25 °C. For further conditions, see Materials and Methods section.
Table A3. Hydrophobicity (log Pow,exp) of studied compounds in formate buffers.
Table A3. Hydrophobicity (log Pow,exp) of studied compounds in formate buffers.
Compoundlog Pow,exp
pH
4.55.56.57.58.5
Taxifolin1.776 ± 0.0001.763 ± 0.0001.704 ± 0.0001.572 ± 0.0020.714 ± 0.007
Silychristin A2.059 ± 0.0072.040 ± 0.0002.003 ± 0.0011.925 ± 0.0021.714 ± 0.001
Silydianin2.088 ± 0.0032.071 ± 0.0002.025 ± 0.0001.935 ± 0.0011.734 ± 0.000
Silychristin B2.114 ± 0.0072.096 ± 0.0002.060 ± 0.0001.986 ± 0.0021.778 ± 0.006
Dehydrosilydianin2.187 ± 0.0012.175 ± 0.0002.154 ± 0.0002.104 ± 0.0001.872 ± 0.001
Quercetin2.279 ± 0.0012.266 ± 0.0002.248 ± 0.0012.203 ± 0.0011.895 ± 0.000
Silybin A2.318 ± 0.0012.302 ± 0.0002.272 ± 0.0002.200 ± 0.0011.920 ± 0.000
Silybin B2.350 ± 0.0012.335 ± 0.0002.304 ± 0.0012.229 ± 0.0031.959 ± 0.000
Dehydrosilychristin2.433 ± 0.0022.421 ± 0.0002.404 ± 0.0052.383 ± 0.0002.183 ± 0.001
Isosilybin A2.437 ± 0.0012.422 ± 0.0002.397 ± 0.0002.328 ± 0.0002.059 ± 0.000
Isosilybin B2.459 ± 0.0012.444 ± 0.0002.415 ± 0.0002.348 ± 0.0012.081 ± 0.001
Dehydrosilybin2.947 ± 0.0032.936 ± 0.0012.928 ± 0.0012.898 ± 0.0072.661 ± 0.010
Dehydroisosilybin3.161 ± 0.0093.143 ± 0.0023.134 ± 0.0043.096 ± 0.0042.815 ± 0.004
The letters A or B in the name of a compound mean that diastereomers of the compound exist and they were separated by a chromatographic system under the applied experimental conditions. Data are presented as mean ± S.D., n = 4.
Table A4. Hydrophobicity (log Pow,exp) of studied compounds in acetate buffers.
Table A4. Hydrophobicity (log Pow,exp) of studied compounds in acetate buffers.
Compoundlog Pow,exp
pH
4.55.56.57.58.5
Taxifolin1.774 ± 0.0001.765 ± 0.0011.710 ± 0.0011.577 ± 0.0060.783 ± 0.013
Silychristin A2.051 ± 0.0012.042 ± 0.0001.999 ± 0.0011.920 ± 0.0071.733 ± 0.000
Silydianin2.083 ± 0.0012.075 ± 0.0052.017 ± 0.0001.928 ± 0.0081.757 ± 0.000
Silychristin B2.106 ± 0.0012.097 ± 0.0012.054 ± 0.0001.980 ± 0.0071.793 ± 0.000
Dehydrosilydianin2.183 ± 0.0012.178 ± 0.0022.151 ± 0.0012.089 ± 0.0081.894 ± 0.001
Quercetin2.273 ± 0.0012.268 ± 0.0002.245 ± 0.0012.202 ± 0.0091.919 ± 0.003
Silybin A2.315 ± 0.0002.307 ± 0.0062.257 ± 0.0012.185 ± 0.0081.942 ± 0.001
Silybin B2.347 ± 0.0012.337 ± 0.0012.289 ± 0.0012.217 ± 0.0061.979 ± 0.000
Dehydrosilychristin2.429 ± 0.0022.424 ± 0.0012.405 ± 0.0042.390 ± 0.0042.203 ± 0.003
Isosilybin A2.434 ± 0.0012.428 ± 0.0052.377 ± 0.0012.310 ± 0.0102.079 ± 0.000
Isosilybin B2.455 ± 0.0012.446 ± 0.0002.396 ± 0.0012.330 ± 0.0102.102 ± 0.000
Dehydrosilybin2.943 ± 0.0012.940 ± 0.0042.909 ± 0.0032.888 ± 0.0152.677 ± 0.000
Dehydroisosilybin3.153 ± 0.0023.148 ± 0.0073.110 ± 0.0023.074 ± 0.0212.840 ± 0.005
The letters A or B in the name of a compound mean that diastereomers of the compound exist and they were separated by a chromatographic system under the applied experimental conditions. Data are presented as mean ± S.D., n = 4.

Appendix B

Figure A1. Relation between (A) calculated log Pow,calc and experimental log Pow,exp measured in phosphate buffer at pH 4.5–8.5; (B) calculated log Plw,calc values and experimental log Pow,exp measured in phosphate buffer at pH 4.5–8.5.
Figure A1. Relation between (A) calculated log Pow,calc and experimental log Pow,exp measured in phosphate buffer at pH 4.5–8.5; (B) calculated log Plw,calc values and experimental log Pow,exp measured in phosphate buffer at pH 4.5–8.5.
Molecules 24 00061 g0a1
Table A5. Correlation coefficient (r2) and mean differences (MD) of n-octanol/water partition coefficient (log Pow,calc) to log Pow,exp measured at different pH and calculated lipid/water partition coefficient (log Plw,calc) in model membranes.
Table A5. Correlation coefficient (r2) and mean differences (MD) of n-octanol/water partition coefficient (log Pow,calc) to log Pow,exp measured at different pH and calculated lipid/water partition coefficient (log Plw,calc) in model membranes.
ParameterpHlog Pow,calc log Plw,calc
COSMOthermMolinspirationAlogP CERDOPCCLC
r24.50.680.830.74 0.820.840.79
5.50.680.830.74 0.820.840.79
6.50.670.870.74 0.790.810.75
7.50.660.870.77 0.810.820.76
8.50.610.710.87 0.780.770.75
MD4.50.830.780.21 0.260.490.69
5.50.840.780.21 0.260.500.69
6.50.900.660.19 0.240.580.62
7.51.020.470.27 0.320.750.51
8.51.080.350.32 0.370.900.44
Table A6. Correlation coefficients (r2) between experimental n-octanol/water partition coefficients (log Pow,exp) and skin intake at pH 6.5 and 8.5.
Table A6. Correlation coefficients (r2) between experimental n-octanol/water partition coefficients (log Pow,exp) and skin intake at pH 6.5 and 8.5.
Log of Skin Intake, pH 6.5 (without DSB)Log of Skin Intake, pH 8.5
% of Ethanol in Donor Liquid% of Ethanol in Donor Liquid
pH051015051015
Log Pow,exp6.50.54 (0.94)0.55 (0.88)0.57 (0.91)0.60 (0.85)----
8.5----0.770.630.830.67
- Only results obtained at the same pH were included in the correlation.
Table A7. Correlation coefficients (r2) between skin intake and calculated permeability (log Perm), lipid/water partition coefficient (log Plw,calc) and penetration barrier (ΔGpen) in membrane models.
Table A7. Correlation coefficients (r2) between skin intake and calculated permeability (log Perm), lipid/water partition coefficient (log Plw,calc) and penetration barrier (ΔGpen) in membrane models.
Skin Intake, log Units, pH 6.5 (without QU) Skin Intake, log Units, pH 8.5
% of Ethanol in Donor Liquid% of Ethanol in Donor Liquid
051015051015
log PermCER0.950.930.820.940.940.940.760.63
DOPC0.870.690.610.750.680.790.570.40
CLC0.940.920.830.940.950.930.790.66
log Plw,calcCER0.610.660.870.580.800.520.740.75
DOPC0.840.870.850.860.970.830.810.76
CLC0.760.810.900.760.920.710.840.79
ΔGpen (kcal/mol)CER0.840.800.560.890.820.940.650.53
DOPC0.820.890.640.880.820.870.500.44
CLC0.830.790.560.880.820.920.610.52
CER—ceramide, DOPC—dioleoylphosphatidylcholine, CLC—ceramide, lignoceric acid and cholesterol mixture. The cells are coloured according to the r2 values based on the displayed colour scale. The linear regression equations for skin intake at pH 8.5 on model membranes were: CER: log (Intake to skin) = 0.76 * log Plw,calc,CER + 2.43. DOPC: log (Intake to skin) = 0.63 * log Plw,calc,DOPC + 2.36. CLC: log (Intake to skin) = 0.54 * log Plw,calc,CLC + 3.20.

Appendix C

Calculation Details:

In order to obtain a computational prediction of the hydrophobicity of flavonolignans, TA, and QU, the n-octanol-water partition coefficients were calculated (log Pow,calc) with the commercial software COSMOtherm 15 software [18] based on quantum-chemical calculations and by the freely available online tools AlogP [15] and Molinspiration [16].
First, the molecular structures were downloaded as pdf files from Pubchem [41]. For COSMOtherm calculations, the same workflow was used as recently in Navrátilová et al. [42]. Individual conformers were generated using the Schrodinger software package [43,44] and for each conformer its vacuum energy and sigma-surfaces were calculated. These were calculated at the BP/TZVP-fine level of theory both for vacuum and cosmo solvent calculation with Turbomole 6.3 [45]. For the Log Pow calculation, COSMOtherm provides a pre-defined setup for partitioning between two fluid phases, where octanol is modelled as a wet octanol, meaning 27.4% of its mass is water, reflecting experimentally more relevant conditions than using dry 100% octanol. The partitioning calculation was performed at 25 and 32 °C. Generally, COSMOtherm calculates the thermodynamic properties of fluid phases, including the partition of substances between variously composed phases. However, a full description of COSMOtherm’s principles is outside the scope of this paper and was described; e.g., by Klamt [46].
For lipid/water partition coefficient calculations, we used COSMOmic [39]. The semi-continuous tool of COSMOtherm uses a lipid bilayer structure/s as an input, separates them into layers and calculates the partition coefficient of the drug into each location (into each layer taking into account the possibility of interactions of the drug with neighboring layers by its rotation). The partition coefficients into individual locations can be easily converted into a free energy profile or into the general lipid/water partition coefficient log Plw. Fluid dioleoylphosphatidylcholine (DOPC) membrane models as well as skin-related solid membrane models were used. The DOPC membrane model is well established and the calculated overall partitioning into phosphatidylcholine membranes reproduces the experiment well. The solid membrane models, ceramide-based ones, do not reflect the basic COSMOtherm requirement for work with fluid phases, but can be also used for estimation drug/lipid interactions, taking into account the need for free energy barrier re-evaluation.
Our previous MD simulations of lipid bilayers consisted of DOPC, ceramide NS24 (CER) and an equimolar mixture of ceramide NS24, lignoceric acid and cholesterol (referred to below as the CLC-mixture) were used. The distribution of individual lipid and water atoms during the last 100 ns of simulations and the averaged atomic distributions of each of the membranes were monitored. Further, cosmo surfaces (the inducted charge distribution of surface segments) on each of the lipid types (BP/TZVP-fine level of theory) were calculated and COSMOmic was run. The permeability coefficient log Perm was also extracted from the free energy profiles on the membranes obtained with COSMOmic.

References

  1. Abd, E.; Yousef, S.A.; Pastore, M.N.; Telaprolu, K.; Mohammed, Y.H.; Namjoshi, S.; Grice, J.E.; Roberts, M.S. Skin models for the testing of transdermal drugs. Clin. Pharmacol. 2016, 8, 163–176. [Google Scholar] [CrossRef] [PubMed]
  2. N’Da, D.D. Prodrug strategies for enhancing the percutaneous absorption of drugs. Molecules 2014, 19, 20780–20807. [Google Scholar] [CrossRef] [PubMed]
  3. Bijak, M. Silybin, a Major Bioactive Component of Milk Thistle (Silybum marianum L. Gaernt.)—Chemistry, Bioavailability, and Metabolism. Molecules 2017, 22, 1942. [Google Scholar] [CrossRef] [PubMed]
  4. Chambers, C.S.; Holečková, V.; Petrásková, L.; Biedermann, D.; Valentová, K.; Buchta, M.; Křen, V. The silymarin composition and why does it matter??? Food Res. Int. 2017, 100, 339–353. [Google Scholar] [CrossRef] [PubMed]
  5. Biedermann, D.; Buchta, M.; Holečková, V.; Sedlák, D.; Valentová, K.; Cvačka, J.; Bednárová, L.; Křenková, A.; Kuzma, M.; Škuta, C.; et al. Silychristin: Skeletal Alterations and Biological Activities. J. Nat. Prod. 2016, 79, 3086–3092. [Google Scholar] [CrossRef]
  6. Chambers, C.S.; Valentova, K.; Kren, V. “Non-Taxifolin” Derived Flavonolignans: Phytochemistry and Biology. Curr. Pharm. Des. 2015, 21, 5489–5500. [Google Scholar] [CrossRef] [PubMed]
  7. Csupor, D.; Csorba, A.; Hohmann, J. Recent advances in the analysis of flavonolignans of Silybum marianum. J. Pharm. Biomed. Anal. 2016, 130, 301–317. [Google Scholar] [CrossRef] [PubMed]
  8. Singh, R.P.; Agarwal, R. Mechanisms and preclinical efficacy of silibinin in preventing skin cancer. Eur. J. Cancer 2005, 41, 1969–1979. [Google Scholar] [CrossRef] [PubMed]
  9. Vaid, M.; Katiyar, S.K. Molecular mechanisms of inhibition of photocarcinogenesis by silymarin, a phytochemical from milk thistle (Silybum marianum L. Gaertn.) (Review). Int. J. Oncol. 2010, 36, 1053–1060. [Google Scholar] [CrossRef] [PubMed]
  10. Han, M.H.; Yoon, W.K.; Lee, H.; Han, S.B.; Lee, K.; Park, S.K.; Yang, K.H.; Kim, H.M.; Kang, J.S. Topical application of silymarin reduces chemical-induced irritant contact dermatitis in BALB/c mice. Int. Immunopharmacol. 2007, 7, 1651–1658. [Google Scholar] [CrossRef] [PubMed]
  11. Altaei, T. The treatment of melasma by silymarin cream. BMC Dermatol. 2012, 12, 1–6. [Google Scholar] [CrossRef]
  12. Sharifi, R.; Rastegar, H.; Kamalinejad, M.; Dehpour, A.R.; Tavangar, S.M.; Paknejad, M.; Mehrabani Natanzi, M.; Ghannadian, N.; Akbari, M.; Pasalar, P. Effect of topical application of silymarin (Silybum marianum) on excision wound healing in albino rats. Acta Med. Iran. 2012, 50, 583–588. [Google Scholar] [PubMed]
  13. OECD Guideline for the Testing of Chemicals. Test No. 107: Partition Coefficient (n octanol/water), Shake Flask Method. 27 July 1995. Available online: https://www.oecd-ilibrary.org/environment/test-no-107-partition-coefficient-n-octanol-water-shake-flask-method_9789264069626-en (accessed on 20 May 2018).
  14. OECD Guidelines for the Testing of Chemicals, Section 1, Physical-Chemical properties. Test No. 117: Partition Coefficient (n-octanol/water), HPLC Method. 23 November 2004. Available online: https://www.oecd-ilibrary.org/environment/test-no-117-partition-coefficient-n-octanol-water-hplc-method_9789264069824-en (accessed on 20 May 2018).
  15. Tetko, I.V.; Gasteiger, J.; Todeschini, R.; Mauri, A.; Livingstone, D.; Ertl, P.; Palyulin, V.A.; Radchenko, E.V.; Zefirov, N.S.; Makarenko, A.S.; et al. Virtual computational chemistry laboratory-design and description. J. Comput. Aided Mol. Des. 2005, 19, 453–463. [Google Scholar] [CrossRef]
  16. Molinspiration Cheminformatics. Available online: http://www.molinspiration.com/ (accessed on 20 May 2018).
  17. Tetko, I.V.; Poda, G.I. Prediction of Log P with Property-Based Methods. In Molecular Drug Properties: Measurement and Prediction; Wiley-VCH: Weinheim, Germany, 2007; pp. 381–406. ISBN 978-3-527-31755-4. [Google Scholar]
  18. Eckert, F.; Klamt, A. Cosmologic GmbH & Co. KG: Leverkusen, Germany, 2013.
  19. Valkó, K. Application of high-performance liquid chromatography based measurements of lipophilicity to model biological distribution. J. Chromatogr. A 2004, 1037, 299–310. [Google Scholar] [CrossRef] [PubMed]
  20. Maroziene, A.; Kliukiene, R.; Sarlauskas, J.; Cenas, N. Inhibition of phthalocyanine-sensitized photohemolysis of human erythrocytes by polyphenolic antioxidants: Description of quantitative structure-activity relationships. Cancer Lett. 2000, 157, 39–44. [Google Scholar] [CrossRef]
  21. Rothwell, J.A.; Day, A.J.; Morgan, M.R. Experimental determination of octanol-water partition coefficients of quercetin and related flavonoids. J. Agric. Food Chem. 2005, 53, 4355–4360. [Google Scholar] [CrossRef]
  22. Gazák, R.; Svobodová, A.; Psotová, J.; Sedmera, P.; Prikrylová, V.; Walterová, D.; Kren, V. Oxidised derivatives of silybin and their antiradical and antioxidant activity. Bioorg. Med. Chem. 2004, 12, 5677–5687. [Google Scholar] [CrossRef]
  23. Zeng, Q.P.; Liu, Z.H.; Huang, A.W.; Zhang, J.; Song, H.T. Preparation and characterization of silymarin synchronized-release microporous osmotic pump tablets. Drug Des. Devel. Ther. 2016, 10, 519–531. [Google Scholar] [CrossRef]
  24. Meloun, M.; Burkoňová, D.; Syrový, T.; Vrána, A. Thermodynamic dissociation constants of silychristin, silybin, silydianin and mycophenolate by the regression analysis of spectrophotometric data. Anal. Chim Acta. 2003, 486, 125–141. [Google Scholar] [CrossRef]
  25. Shubina, V.S.; Shatalina, Y.V. Absorption Spectroscopy Study of Acid-Base and Metal-Binding Properties of Flavanones. J. Appl. Spectrosc. 2013, 80, 761–766. [Google Scholar] [CrossRef]
  26. Herrero-Martínez, J.M.; Sanmartin, M.; Rosés, M.; Bosch, E.; Ràfols, C. Determination of dissociation constants of flavonoids by capillary electrophoresis. Electrophoresis 2005, 26, 1886–1895. [Google Scholar] [CrossRef] [PubMed]
  27. Garrido, N.M.; Queimada, A.J.; Jorge, M.; Macedo, E.A.; Economou, I.G. 1-Octanol/Water Partition Coefficients of N-Alkanes from Molecular Simulations of Absolute Solvation Free Energies. J. Chem. Theory Comput. 2009, 5, 2436–2446. [Google Scholar] [CrossRef] [PubMed]
  28. Sun, H. (Ed.) Quantitative Structure-Property Relationships Models for Lipophilicity and Aqueous Solubility. In A Practical Guide to Rational Drug Design; Elsevier: Cambridge, UK, 2016; pp. 193–223. [Google Scholar]
  29. Ekins, S.; Mestres, J.; Testa, B. In Silico Pharmacology for Drug Discovery: Methods for Virtual Ligand Screening and Profiling. Br. J. Pharmacol. 2007, 152, 9–20. [Google Scholar] [CrossRef] [PubMed]
  30. Mannhold, R.; Poda, G.I.; Ostermann, C.; Tetko, I.V. Calculation of molecular lipophilicity: State-of-the-art and comparison of log P methods on more than 96,000 compounds. J. Pharm. Sci. 2009, 98, 861–893. [Google Scholar] [CrossRef] [PubMed]
  31. Benson, H.A.E. Skin Structure, Function, and Permeation. In Topical and Transdermal Drug Delivery: Principles and Practice; Benson, H.A.E., Watkinson, A.C., Eds.; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2012. [Google Scholar]
  32. Wang, W.; Sun, C.; Mao, L.; Ma, P.; Liu, F.; Yang, L.; Gao, Y. The biological activities, chemical stability, metabolism and delivery systems of quercetin: A review. Trends Food Sci. Technol. 2016, 56, 21–38. [Google Scholar] [CrossRef]
  33. Hung, C.F.; Lin, Y.K.; Zhang, L.W.; Chang, C.H.; Fang, J.Y. Topical delivery of silymarin constituents via the skin route. Acta Pharmacol. Sin. 2010, 31, 118–126. [Google Scholar] [CrossRef] [PubMed]
  34. Jung, E.C.; Maibach, H.I. Animal models for percutaneous absorption. J. Appl. Toxicol. 2015, 35, 1–10. [Google Scholar] [CrossRef]
  35. Monti, D.; Gažák, R.; Marhol, P.; Biedermann, D.; Purchartová, K.; Fedrigo, M.; Riva, S.; Křen, V. Enzymatic Kinetic Resolution of Silybin Diastereoisomers. J. Nat. Prod. 2010, 73, 613–619. [Google Scholar] [CrossRef]
  36. Křenek, K.; Marhol, P.; Peikerová, Ž.; Křen, V.; Biedermann, D. Preparatory separation of the silymarin flavonolignans by Sephadex LH-20 gel. Food Res. Int. 2014, 65, 115–120. [Google Scholar] [CrossRef]
  37. Gažák, R.; Trouillas, P.; Biedermann, D.; Fuksová, K.; Marhol, P.; Kuzma, M.; Křen, V. Base-catalyzed oxidation of silybin and isosilybin into 2,3-dehydro derivatives. Tetrahedron Lett. 2013, 54, 315–317. [Google Scholar] [CrossRef]
  38. Pyszková, M.; Biler, M.; Biedermann, D.; Valentová, K.; Kuzma, M.; Vrba, J.; Ulrichová, J.; Sokolová, R.; Mojovic, M.; Popovic-Bijelic, A.; et al. Flavonolignan 2,3-dehydroderivatives: Preparation, antiradical and cytoprotective activity. Free Rad. Biol. Med. 2016, 90, 114–125. [Google Scholar] [CrossRef]
  39. Klamt, A.; Huniar, U.; Spycher, S.; Keldenich, J. COSMOmic: A mechanistic approach to the calculation of membrane-water partition coefficients and internal distributions within membranes and micelles. J. Phys. Chem. B 2008, 112, 12148–12157. [Google Scholar] [CrossRef] [PubMed]
  40. Ingram, T.; Storm, S.; Kloss, L.; Mehling, T.; Jakobtorweihen, S.; Smirnova, I.V. Prediction of micelle/water and liposome/water partition coefficients based on molecular dynamics simulations, COSMO-RS, and COSMOmic. Langmuir. 2013, 29, 3527–3537. [Google Scholar] [CrossRef] [PubMed]
  41. Pubchem. Available online: https://pubchem.ncbi.nlm.nih.gov/ (accessed on 20 May 2018).
  42. Navrátilová, V.; Paloncýová, M.; Berka, K.; Mise, S.; Haga, Y.; Matsumura, C.; Sakaki, T.; Inui, H.; Otyepka, M. Molecular Insights into the Role of a Distal F240A Mutation That Alters CYP1A1 Activity towards Persistent Organic Pollutants. BBA Gen. Subj. 2017, 1861, 2852–2860. [Google Scholar] [CrossRef]
  43. Schrödinger Release 2015-1, Maestro. Schrödinger; LLC: New York, NY, USA, 2015.
  44. Small-Molecule Drug Discovery Suite 2015-4, Schrödinger; LLC: New York, NY, USA, 2016.
  45. TURBOMOLE: Program Package for ab initio Electronic Structure Calculations. Available online: http://www.turbomole.com (accessed on 20 May 2018).
  46. Klamt, A. The COSMO and COSMO-RS Solvation Models. Wiley Interdiscip. Rev. Comput. Mol. Sci. 2011, 1, 699–709. [Google Scholar] [CrossRef]
Sample Availability: Studied compounds are available from the authors.
Figure 1. Structure of studied polyphenols.
Figure 1. Structure of studied polyphenols.
Molecules 24 00061 g001
Figure 2. Relation of compounds’ skin intake to their lipophilicity. The relation of skin intake to experimental n-octanol/water partition coefficient (log Pow,exp) at pH 6.5 is presented on the left; the relation of skin intake to log Pow,exp at pH 8.5 is shown on the right.
Figure 2. Relation of compounds’ skin intake to their lipophilicity. The relation of skin intake to experimental n-octanol/water partition coefficient (log Pow,exp) at pH 6.5 is presented on the left; the relation of skin intake to log Pow,exp at pH 8.5 is shown on the right.
Molecules 24 00061 g002
Figure 3. Free energy profiles of studied compound in ceramide (CER), dioleoylphosphatidylcholine (DOPC) and ceramide, lignoceric acid and cholesterol mixture (CLC) membrane models.
Figure 3. Free energy profiles of studied compound in ceramide (CER), dioleoylphosphatidylcholine (DOPC) and ceramide, lignoceric acid and cholesterol mixture (CLC) membrane models.
Molecules 24 00061 g003
Figure 4. Relation between compounds’ skin intake and calculated permeability or partition coefficients in lipid membrane models. The relation between the logarithms of the compounds’ skin intake at pH 6.5 and calculated permeability (log Perm) in lipid membrane models is presented on the left; the logarithms of the compounds’ skin intake at pH 8.5 and calculated lipid membrane models partition coefficients (log Plw,calc) is shown on the right. CER—ceramide, DOPC—dioleoylphosphatidylcholine, CLC—ceramide, lignoceric acid and cholesterol mixture.
Figure 4. Relation between compounds’ skin intake and calculated permeability or partition coefficients in lipid membrane models. The relation between the logarithms of the compounds’ skin intake at pH 6.5 and calculated permeability (log Perm) in lipid membrane models is presented on the left; the logarithms of the compounds’ skin intake at pH 8.5 and calculated lipid membrane models partition coefficients (log Plw,calc) is shown on the right. CER—ceramide, DOPC—dioleoylphosphatidylcholine, CLC—ceramide, lignoceric acid and cholesterol mixture.
Molecules 24 00061 g004
Figure 5. Structure of model membrane and free energy profile of penetration barrier. Lipid membrane constituents: the hydrocarbon chains in the membrane core are represented as cyan and white sticks; oxygen, phosphorus and nitrogen atoms in the DOPC membrane are shown as red, olive and blue balls; water molecules are represented as grey balls. The graph shows the free energy profile of the penetration barrier (ΔGpen) and the free energy at depth z ΔG(z).
Figure 5. Structure of model membrane and free energy profile of penetration barrier. Lipid membrane constituents: the hydrocarbon chains in the membrane core are represented as cyan and white sticks; oxygen, phosphorus and nitrogen atoms in the DOPC membrane are shown as red, olive and blue balls; water molecules are represented as grey balls. The graph shows the free energy profile of the penetration barrier (ΔGpen) and the free energy at depth z ΔG(z).
Molecules 24 00061 g005
Table 1. Hydrophobicity (log Pow,exp) of studied compounds in phosphate buffers.
Table 1. Hydrophobicity (log Pow,exp) of studied compounds in phosphate buffers.
Compoundlog Pow,exp
pH
4.55.56.57.58.5
Taxifolin1.761 ± 0.0171.763 ± 0.0021.674 ± 0.0171.486 ± 0.0260.802 ± 0.046
Silychristin A2.045 ± 0.0102.042 ± 0.0031.952 ± 0.0171.786 ± 0.0211.705 ± 0.003
Silydianin2.076 ± 0.0112.077 ± 0.0101.964 ± 0.0191.805 ± 0.0171.742 ± 0.004
Silychristin B2.100 ± 0.0082.097 ± 0.0032.003 ± 0.0181.837 ± 0.0201.758 ± 0.002
Dehydrosilydianin2.178 ± 0.0092.181 ± 0.0082.108 ± 0.0141.948 ± 0.0241.843 ± 0.004
Quercetin2.269 ± 0.0102.266 ± 0.0062.197 ± 0.0151.979 ± 0.0301.778 ± 0.005
Silybin A2.314 ± 0.0082.310 ± 0.0102.180 ± 0.0241.972 ± 0.0251.872 ± 0.002
Silybin B2.346 ± 0.0062.336 ± 0.0092.203 ± 0.0191.999 ± 0.0231.904 ± 0.002
Dehydrosilychristin2.425 ± 0.0032.420 ± 0.0072.352 ± 0.0162.162 ± 0.0292.025 ± 0.010
Isosilybin A2.432 ± 0.0022.422 ± 0.0082.289 ± 0.0242.081 ± 0.0231.992 ± 0.002
Isosilybin B2.458 ± 0.0052.444 ± 0.0102.306 ± 0.0262.100 ± 0.0222.014 ± 0.002
Dehydrosilybin2.948 ± 0.0032.932 ± 0.0192.788 ± 0.0362.508 ± 0.0382.353 ± 0.019
Dehydroisosilybin3.154 ± 0.0023.131 ± 0.0032.931 ± 0.0472.619 ± 0.0392.462 ± 0.026
The letters A or B in the name of compounds mean that diastereomers of the compound exist and they were separated by a chromatographic system under the applied experimental conditions. Data are presented as mean ± S.D., n = 4.
Table 2. Hydrophobicity of studied compounds predicted by computational tools as octanol/water partition coefficient (log Pow,calc) and lipid/water partition coefficient (log Plw,calc) in membrane models.
Table 2. Hydrophobicity of studied compounds predicted by computational tools as octanol/water partition coefficient (log Pow,calc) and lipid/water partition coefficient (log Plw,calc) in membrane models.
Compoundlog Pow,calc log Plw,calc
COSMOthermMolinspirationAlogP CERDOPCCLC
Taxifolin1.010.711.07 1.141.490.22
Silychristin A2.551.262.13 1.642.351.10
Silydianin1.291.211.84 2.282.251.51
Silychristin B2.511.262.13 1.762.351.06
Dehydrosilydianin1.071.211.84 2.132.081.34
Quercetin2.171.681.81 1.802.170.79
Silybin A3.301.472.35 2.213.061.93
Silybin B3.201.472.35 2.192.901.83
Dehydrosilychristin3.612.242.57 2.263.091.50
Isosilybin A3.561.472.35 2.823.442.47
Isosilybin B3.251.472.35 2.463.172.16
Dehydrosilybin4.172.442.80 2.783.682.57
Dehydroisosilybin4.182.442.81 3.423.973.09
CER—ceramide, DOPC—dioleoylphosphatidylcholine, CLC—ceramide, lignoceric acid and cholesterol mixture. The letters A or B in the name of a compound mean that diastereomers of the compound exists and they were calculated separately in the computational tools.
Table 3. Skin intake of studied compounds evaluated by transdermal transport.
Table 3. Skin intake of studied compounds evaluated by transdermal transport.
Compound Concentration (pmol g−1 of Skin)
pH 6.5 pH 8.5
EtOH (%, v/v) EtOH (%, v/v)
051015 051015
Taxifolin3101423720682228 2099151812022022
Silychristin16717865554315690 5358516629123275
Silydianin6997660678093457 4966220462576129
Silybin53093397441598431531 1767626824132587311
Quercetin56230786222862849162 65866185109352699
Isosilybin77643716384737529113 319852013185917969
Dehydrosilybin45105451452239935810 50248303522847726721
The studied compounds were applied to the skin membrane at a concentration of 50 µM in the donor liquid and incubated at 32 °C. After 24 h, the skin was collected for evaluating the content of the studied compounds.
Table 4. Permeability (log Perm) and penetration barrier (ΔGpen) calculated from COSMOmic free energy profiles in membrane models.
Table 4. Permeability (log Perm) and penetration barrier (ΔGpen) calculated from COSMOmic free energy profiles in membrane models.
Compoundlog Perm ΔGpen (kcal/mol)
CERDOPCCLC CERDOPCCLC
Taxifolin−4.56−1.44−4.67 9.395.307.78
Silychristin−2.41−0.35−2.53 6.964.695.96
Silydianin−3.40−0.95−3.31 9.296.027.98
Silybin A−1.01−0.28−1.07 5.493.595.08
Silybin B−1.05−0.30−1.13 5.773.455.07
Quercetin−3.11−0.25−3.23 8.244.376.41
Isosilybin A−0.51−0.29−0.56 5.972.885.10
Dehydrosilybin−0.28−0.26−0.24 4.552.884.07
CER—ceramide, DOPC—dioleoylphosphatidylcholine, CLC—ceramide, lignoceric acid and cholesterol mixture. The letters A or B for silybin mean that calculations were done for its respective diastereomers.

Share and Cite

MDPI and ACS Style

Kosina, P.; Paloncýová, M.; Rajnochová Svobodová, A.; Zálešák, B.; Biedermann, D.; Ulrichová, J.; Vostálová, J. Dermal Delivery of Selected Polyphenols from Silybum marianum. Theoretical and Experimental Study. Molecules 2019, 24, 61. https://doi.org/10.3390/molecules24010061

AMA Style

Kosina P, Paloncýová M, Rajnochová Svobodová A, Zálešák B, Biedermann D, Ulrichová J, Vostálová J. Dermal Delivery of Selected Polyphenols from Silybum marianum. Theoretical and Experimental Study. Molecules. 2019; 24(1):61. https://doi.org/10.3390/molecules24010061

Chicago/Turabian Style

Kosina, Pavel, Markéta Paloncýová, Alena Rajnochová Svobodová, Bohumil Zálešák, David Biedermann, Jitka Ulrichová, and Jitka Vostálová. 2019. "Dermal Delivery of Selected Polyphenols from Silybum marianum. Theoretical and Experimental Study" Molecules 24, no. 1: 61. https://doi.org/10.3390/molecules24010061

APA Style

Kosina, P., Paloncýová, M., Rajnochová Svobodová, A., Zálešák, B., Biedermann, D., Ulrichová, J., & Vostálová, J. (2019). Dermal Delivery of Selected Polyphenols from Silybum marianum. Theoretical and Experimental Study. Molecules, 24(1), 61. https://doi.org/10.3390/molecules24010061

Article Metrics

Back to TopTop