Next Article in Journal
Applying Supercritical Fluid Technology to Prepare Ibuprofen Solid Dispersions with Improved Oral Bioavailability
Previous Article in Journal
Methotrexate-Loaded Solid Lipid Nanoparticles: Protein Functionalization to Improve Brain Biodistribution
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Biopharmaceutical Development of a Bifonazole Multiple Emulsion for Enhanced Epidermal Delivery

by
Joaquim Suñer-Carbó
1,2,
Ana Calpena-Campmany
1,2,
Lyda Halbaut-Bellowa
1,2,
Beatriz Clares-Naveros
2,3,
María José Rodriguez-Lagunas
4,
Elena Barbolini
1,
Joanna Zamarbide-Losada
1 and
Antonio Boix-Montañés
1,2,*
1
Department of Pharmacy and Pharmaceutical Technology and Physical Chemistry, Faculty of Pharmacy and Food Sciences, University of Barcelona, Av. Joan XXIII s/n, 08028 Barcelona, Spain
2
Institute of Nanoscience and Nanotechnology (IN2UB), University of Barcelona, Diagonal 645, 08028 Barcelona, Spain
3
Department of Pharmacy and Pharmaceutical Technology, Faculty of Pharmacy, University of Granada, Campus de la Cartuja, s/n. 18071 Granada, Spain
4
Department of Biochemistry and Physiology, Faculty of Pharmacy and Food Sciences, University of Barcelona, Av. Joan XXIII s/n, 08028 Barcelona, Spain
*
Author to whom correspondence should be addressed.
Submission received: 22 December 2018 / Revised: 29 January 2019 / Accepted: 31 January 2019 / Published: 2 February 2019

Abstract

:
Efficient topical delivery of imidazolic antifungals faces the challenge of overcoming its limited water solubility and its required long-lasting duration of treatments. In this paper, a hydrophilic multiple emulsion (ME) of Bifonazole (BFZ) is shown to maximize its skin retention, minimize its skin permeation, and maintain an acceptable level of being harmless in vivo. The formulations were pharmaceutically characterized and application properties were assessed based on viscosity measurements. Non-Newtonian pseudoplastic shear thinning with apparent thixotropy was observed, facilitating the formulation retention over the skin. The in vitro release profile with vertical diffusion cells showed a predominant square-root release kinetic suggesting an infinite dose depletion from the formulation. Ex vivo human skin permeation and penetration was additionally evaluated. Respective skin permeation was lower than values obtained with a commercial O/W formulation. The combination of amphoteric and non-ionic surfactants increased the bifonazole epidermal accumulation by a factor of twenty. This fact makes the possibility of increasing its current 24 h administration frequency more likely. Eventual alterations of skin integrity caused by the formulations were examined with epidermal histological analysis and in vivo preclinical measurements of skin elasticity and water retrograde permeation. Histological analysis demonstrated that the multiple emulsions were harmless. Additionally, modifications of in vivo skin integrity descriptors were considered as negligible.

Graphical Abstract

1. Introduction

Bifonazole [1-[[1,1’-biphenyl)-4-phenylmethyl]-1H-imidazole) (BFZ) is an imidazolic antifungal drug indicated against skin or mucosal mycoses. Its topical administration requires an intensive drug retention in skin because the pathogenicity of dermatomycoses is related to their ability to gain access to other target tissues [1,2]. Formulations for local therapy are desirable in order to achieve the highest drug penetration deep into the epidermis, which is the site of action. Nevertheless, these drugs are markedly lipophilic and are largely retained in the stratum corneum (SCR) [3], making the treatment of deep-seated epidermal infections more difficult. Topical agents, which are conventionally used for the treatment of skin fungal infections, are usually formulated as creams, lotions, or gels. Several formulation strategies, such as the hydrophilic W/O/W multiple emulsions (ME), can expedite this skin accumulation: They consist of two hydrophilic internal and external phases separated by a lipophilic layer that dissolves hydrophobic drugs and incorporates them into a hydrophilic formulation.
Multiple emulsions have attracted considerable attention in recent years due to their application as potentials systems for drug delivery. In this case, they have been investigated as potential skin delivery systems for bifonazole. Many authors [4] describe MEs as “emulsions of emulsions” since they are formed by phase inversion of a simple emulsion. Thus, at least two surfactants are required: One with low HLB (Hydrophilic Lipophilic Balance) forms the primary W/O emulsion and another surfactant with a high HLB achieves the secondary o/w emulsification. In this sense, polymeric emulsifiers provide a strong, thick, and well-defined film around the W/O interface [5]. Additionally, polysorbates (O/W emulsifiers) and sorbitan esters (W/O emulsifiers) are non-ionic surfactants that enhance the formulation stability and flexibility, and furthermore, they facilitate the skin penetration of the drug [6]. They are stable in acids, do not react with ionic ingredients or actives, and have a history of being safe in their use [7]. In summary, polymeric surfactants in combination with conventional small molecular weight emulsifiers, are suitable emulsifiers for multiple emulsions [8] and are additionally useful to enhance the solubilization of poorly soluble drugs [9].
General advantages of W/O/W multiple emulsions are the selective protection of the entrapped drug, their ability to incorporate different actives in the different emulsion phases and their sustained release effects [10]. They have been widely applied in various fields of pharmacy [11,12], such as vehicles for anticancer agents [13] or prolonged drug delivery systems [14]. In Bifonazole formulations, MEs are more suitable than other similar disperse systems, such as liposomes or micelles: The internal phase of liposomes is hydrophilic and is not able to dissolve significant amounts of lipophilic drugs [15,16] such as Bifonazole, whereas the intermediate lipophilic phase of the ME allows to achieve greater solubilization.
The aim of this study was to evaluate the efficiency of nonionic co-emulsifiers in the outer aqueous phase of a multiple emulsion of Bifonazole for topical application to achieve a lower BFZ human skin permeation and larger skin retention in comparison with commercial O/W formulation. In aggregate, the feasibility of an enhanced (that is lengthened) interval of administration, which maintains an acceptable skin tolerance, has been investigated.

2. Materials and Methods

W/O/W emulsions were prepared with a two-step emulsification method and were physically characterized with optical examination, laser and microscopic droplet size distribution, rheological studies, pH and conductimetric analyses, and they were tested for drug content. In vitro release and ex vivo human skin permeation were also investigated. After a histologic evaluation of its harmlessness, the in vivo skin integrity of the different formulations, including a commercial emulsion (BFZ-CF) was confirmed after their acute administration in humans.

2.1. Substances and Reagents

Bifonazole (BFZ, CAS 60628-96-8), sodium chloride, cetyl palmitate (CP), and polysorbate 80 (Tween® 80, HLB 15) were obtained from Acofarma (Barcelona, Spain). Sorbitan stearate (Span® 60, HLB 4.7) was obtained from Fagron (Terrassa, Spain). The lipophilic surfactant cetyl dimethicone copolyol (Abil® EM90, HLB 5, CDC), an amphoteric surfactant cocamidopropyl betaine (Tego® Betain F, CMB), and the thickening acrylic carbomer (Tego® Carbomer 341ER, TGC) were supplied by Evonik (Essen, Germany). Capric/caprylic triglyceride (Labrafac® Lipophile 1349, LP), were supplied by Gattefossé (Saint-Priest, France). Deionized water (Ω−1 < 0.07 mS/cm) used in all experiments was produced by ourselves in-house in a Milli-Q® system (Millipore Iberica S.A.U., Madrid, Spain). Methanol HPLC was obtained from Fisher Scientific (Madrid, Spain). Ammonium acetate, ethanol, sodium dodecyl sulphate, sodium hydroxide, and potassium monohydrogen phosphate K2HPO4 were obtained from Panreac (Barcelona, Spain). Hank’s balanced salt solution was obtained from Thermo Fisher Scientific (Barcelona, Spain). Reagents for histological preparations were purchased from Sigma and Thermo Fisher Scientific (Barcelona, Spain). Commercial formulation of 1% Bifonazole from Bayer AG (Canesmycospor®, BFZ-CF) was used as a comparable formulation.

2.2. Preparation of Multiple Emulsions

Two, previously optimized, W1/O/W2 multiple emulsions with insoluble imidazolic antifungals [17] were prepared with a two-step emulsification method [18]. Compositions are summarized in Table 1.
Primary W1/O emulsion was prepared by slow addition of the aqueous phase containing the electrolyte (NaCl) to the oil phase containing BFZ (1%, w/w) at 80 ± 2 °C under continuous stirring at 250 rpm until approximately 25 °C. The oil phase was prepared dissolving BFZ in a combination of drug solvent and cosolvent (LP and CP) aided with the lipophilic emulsifying agents (CDC and sorbitan stearate) at 80 ± 2 °C. In the second stage, the obtained primary emulsions were slowly added to the corresponding outer aqueous phases (W2) 50:50 under 250 rpm at room temperature. After complete addition of the primary W1/O emulsion over the external gelified aqueous phase, the resulting mix was paddle stirred for a further 10 min until a homogeneous W1/O/W2 multiple emulsion had been completely formed.
External aqueous phases (W2) had been prepared previously by dispersing the cross-linked TGC polymer, in a co-solvent system of deionized water and the hydrophilic emulsifying agents (CMB and/without polysorbate 80) were neutralized with 10% NaOH solution (10%, w/v) to obtain a pH value of 6.5–7.0.

2.3. Physicochemical Characterization and Stability

Organoleptic characteristics were monitored to detect any visible signs of instability such as creaming, phase separation, or color changes. Droplet homogeneity and size distribution were investigated by optical microscopy (400×) with a Leica DM 1000 LED light microscope (Leica Microsystems, Wetzlar, Germany). Measurements were performed directly and at 1:20 water dilution, just after formulation and monthly during a six-month stability storage period. Droplet size was also analyzed by laser diffractometry with a Malvern Mastersizer 2000 (Malvern Instruments, Worcestershire, UK). The relative distribution in size ranges was calculated based on Mie theory and reported as D[v,0.1], D[v, 0.5], and D[v,0.9], i.e., standard percentile readings of the distribution. The volume-weighted mean diameter, D[4,3] and the volume-surface mean diameter, D[3,2], were also obtained [19]. They were measured (n = 3, water-dispersed, room temperature) at 24 h and then at 180 days after preparation. Conductivities (μS/cm) of the 1:20 diluted W1/O/W2 emulsions were measured (n = 3) with a Crison GLP 30 (Crison Instruments, Alella, Spain) at 24 h and 180 days after preparation. pH was measured (Crison Instruments, Alella, Spain) at the same time (n = 3, room temperature). All results were compared for statistical significance with Student’s t-test (p < 0.05).
The occurrence of destabilization phenomena was assessed by multiple light scattering with a Turbiscan® Lab Expert (Formulaction, L’Union, France). Diluted samples (~35 mL) were placed into cylindrical glass measuring cells, and were completely scanned for light scattering by a reading head that consists of a pulsed near-infrared light source (λ = 880 nm) and two synchronous optical detectors. A pattern of the backscattered light as a function of the sample height was obtained giving a macroscopic fingerprint of the sample at each predetermined time. Measurements (n = 3) were performed at 24 h after preparation and after six months storage at room temperature. Thermal stability tests were also performed under different storage conditions for both primary and multiple emulsions. Test tubes were stored vertically and kept at 4 ± 0.1 °C (in refrigerator), 30 ± 0.1 °C, 40 ± 0.1 °C in stability cabins (Heraeus B5042 E and I 42, Madrid, Spain) and room temperature. Observations were made each week during storage.

2.4. Rheological Properties

Skin application properties (apparent thixotropy and elasticity) were characterised with rotational and oscillatory tests. Rheological measurements were performed with a Haake Rheostress 1® rheometer (Thermo Fisher Scientific, Karlsruhe, Germany) connected to a thermostatic circulator Thermo Haake Phoenix II + Haake C25P. Data were analyzed with Haake Rheowin® Data Manager v. 3.3 software (Thermo Electron Corporation, Karlsruhe, Germany).

2.4.1. Rotational Measurements

Steady-state measurements were addressed with a cone and a mobile upper cone (C60/2° Ti:60 mm diameter, 2° angle). The shear stress (τ) was measured as a function of the shear rate (γ). Viscosity curves (η = f(γ)) and flow curves (τ = f(γ)) were recorded at 25 ± 0.1 °C. The shear rate ramp program included a 3 min ramp-up period from 0 to 100 s−1, 1 min constant shear rate period at 100 s−1, and 3 min ramp-down from 100 to 0 s−1. Representative mathematical models were fitted to flow curves (when they were non-Newtonian) searching for the best descriptive model: Bingham, Ostwald de Waele, Herschel-Bulkley, Casson, and Cross [20]. Selection of the best fitting was based on the correlation coefficient (observed vs. predicted) and chi-square value. The apparent thixotropy (Pa/s) was estimated as the area of hysteresis loop. Steady-state viscosity (η, Pa s) was determined from the constant shear section at 100 s−1 and also from the ramp-up period at 10 s−1.

2.4.2. Dynamic Oscillatory Measurements

Oscillatory tests were performed with a parallel plate-and-plate geometry (Haake PP60 Ti, 60 mm diameter, 0.5 mm gap) using the same rheometer described above. Firstly, the oscillatory stress sweep test was performed at a constant frequency of 1 s−1 from 0.01 to 100 Pa to determine the linear viscoelastic region (LVR). Afterwards, frequency sweep test was carried out between a 0.01 and 10 s−1, at a constant shear rate within the LVR. This was in order to determine the related variations of the storage modulus (G′), loss modulus (G″), phase angle (δ), and the complex viscosity (η*). During each sweep stress test, the changes in storage and loss modulus and the phase angle (G’, G’’ and δ) were plotted as a function of shear stress at a frequency of 1 Hz.

2.5. Drug Content

2.5.1. Drug Extraction

Formulation amounts equivalent to about 5 mg BFZ were mixed with 5 mL of methanol and dispersed in an ultrasound bath for 10 min. Afterwards, the sample was centrifuged at 3000 rpm/10 min and the resulting clear solution was transferred to a vial. After appropriate dilutions with methanol/ pH 7.4 buffer phosphates (75:25), the resulting solution was filtered through 0.45 µm nylon disposable filter (Teknokroma, Barcelona) and measured spectrophotometrically for bifonazole content.

2.5.2. UV Analysis

BFZ concentrations in assay formulations and release samples (see Section 2.6) were measured by UV-spectrophotometry [21] at λ = 254 nm (Thermo® Spectronic Helios Beta, Thermo, Cambridge, UK). Calibration curves were prepared with methanol:pH 7.4 phosphates buffer 75:25 (v/v). Validation according to standard guidelines [22] demonstrated an adequate intraassay and interassay precision (maximum RSD +4.36%) and accuracy (from −3.12 to +2.58%). Linearity ranged between 2.61 and 10.00 μg/mL. The statistically significant determination coefficients (R2) were obtained on each run as a linearity indicator on each run. The recovery in formulation assays was tested in triplicate using a 1% BFZ methanolic solution (MethSol), giving a result of 93.4% (RSD 7.6%) and this was considered acceptable for our purposes.

2.5.3. HPLC Analysis

BFZ concentration in skin permeation samples (see Section 2.7) was quantified with a Waters Alliance 2695 HPLC (Waters, Milford, MA, USA) with a C18 column Brisa LC2, 5μm, 4.6 mm × 250 mm (Teknokroma, Barcelona, Spain). The mobile phase was methanol:ammonium acetate 65 mM (65:35, pH 3.6) [23] at a flow rate of 1 mL/min at 25 °C. Detection was monitored at 220 nm from 1.56 to 100.00 μg/mL. Accuracy was lower than +3.7%. As for the precision, maximum RSD was 0.8% with an acceptable linearity falling in the calibration range. Statistically significant determination coefficients (R2) were obtained on each run as a linearity indicator on each run.

2.5.4. Statistical Analyses

Descriptive statistics, confidence intervals and regressions were calculated with the Graphpad® Prism software v.3 (Graphpad Software, San Diego, CA, USA).

2.6. In Vitro Drug Release

In vitro BFZ drug release was tested from 300 mg donor amounts of each of the following, ME, BFZ-CF and MethSol with 2.54 cm2 vertical diffusion cells [24] using 0.45µm mesh nylon membrane (Teknokroma, Barcelona, Spain). The receptor solution was a well-stirred solution of 12 mL Methanol:K2HPO4 Buffer 0.05M 75:25 (v/v) at 32 ± 0.5 °C. Samples from the receptor compartment (300 µL) were taken during a minimum of 5 h with immediate spectrophotometric analysis and blank reposition. Comparative release profiles were obtained (n = 3) for each formulation.

2.6.1. Drug Release Equations

Zero order, first order, Hyperbola, Weibull, and Higuchi equations [25] were fitted to the cumulative released BFZ amounts using a non-linear regression routine with Graphpad® Prism 3 software (Graphpad Software, San Diego, CA, USA).
The best model was selected based on the minimum value of the Akaike information criterion (AIC) [26] calculated as follows from the residual sum of squares obtained with each regression equation:
AIC = n·ln SSQ + 2P,
where SSQ is the residual sum of squares and P the number of parameters. Best fitting was considered to be as having the minimum AIC value. Descriptive release parameters were: Release rate K (of the best function), released percentages at five hours (%Q5), and efficiency, described as:
Efficiency = AUC05/(Q5·T),
where AUC05 is the individual area under the curve of released amounts and T is the time for the first asymptotic value [27].

2.6.2. Statistical Analyses

Comparisons of parameters were performed by one-way ANOVA with GraphPad® Prism 3 software (Graphpad Software, Sant Diego, CA, USA).

2.7. Drug Permeation—Penetration

BFZ human skin permeation from each formulation was investigated with 0.64 cm2 vertical diffusion cells [24]. Permeation flux and skin retention were calculated for each replicate (n = 3).
Skin was obtained from leftovers of the abdominal region of a 38-year old healthy woman during plastic surgery (Hospital de Barcelona, SCIAS, Barcelona, Spain) following a donation procedure of biologic residues approved by Bioethics Committee of the Barcelona-SCIAS Hospital (reference number: BEC/001/16). The patient provided written informed consent. After chirurgical excision, skin was immediately debrided, immersed in Hanks’ solution, and stored at −20 °C. Prior to the experiment, skin laminar samples (400 μm) were obtained with an electrical dermatome (Aesculap GA 630, Tuttlingen, Germany). Skin specimens were mounted on diffusion cells and a receptor compartment was filled up with Ethanol:Transcutol® P:Water (50:20:30 v/v/v). After temperature equilibration at 32 °C, the skin integrity was verified measuring its transepidermal water loss (TEWL) (Tewameter® TM300, Köln, Germany) just before formulation application (300 mg) and discarding results falling out of the range 4 to 8 g/h/m2 [28].

2.7.1. Drug Permeation

At predefined times, 300 μL sample were taken with replacement during a minimum of 32 h. Permeated BFZ amounts per unit area (Qt) were obtained by HPLC-UV as described previously. Permeation extent was described with the last time amount (Q32). Permeation flux (J) was calculated by linear regression between Qt and time in the linear phase.

2.7.2. Drug Penetration

At the end of each experiment, the residual amount of donor formulation was withdrawn and the skin specimens were rinsed with aqueous 0.5% sodium laurylsulphate and later deionized water. Resulting skin specimens were blotted dry, weighed, and minced. Drug extraction was run with fresh receptor solution under ultrasounds (×3). Collected solutions were centrifuged (1500 rpm, 10 min) and immediately analyzed. Drug retention levels were reported as the ratio between the amount of drug in skin and the weight of skin sample (µg/g), and as the amount of drug per surface unit (µg/cm2).

2.7.3. Statistical Analyses

Statistical comparisons (α = 0.05) were performed by one-way ANOVA using the GraphPad® Prism 3 software (Graphpad Software, San Diego, CA, USA).

2.8. Epidermal Histology

Histology of epidermal specimens after the skin permeation-penetration experiments (see Section 2.7) was investigated with optical microscopy. The formulations evaluated were: BFZ-CF, JMLP01BT and MEs placebos (without BFZ). After the permeation experiments, residual contents of donor formulations were mechanically withdrawn and skin samples were set overnight in 4% buffered formaldehyde at room temperature and embedded in paraffin, cut into 5 µm thick sections, stained with hematoxylin-eosin and finally they were observed at ×400 with an Olympus BX41 microscope and Olympus XC50 camera (Olympus, Tokyo, Japan).

2.9. Skin Integrity Assessment

Eventual formulation-induced variations of skin integrity were investigated with simple non-invasive in vivo measurements of skin elasticity, stratum corneum hydration (SCH), and transepidermal water loss (TEWL) [28] in an investigator-blinded parallel clinical study.
The study was based on previous knowledge [29] and conducted in accordance with the Declaration of Helsinki [30] and the Ethics Committee on human experimentation of the University of Barcelona (Reference number: iRB00003099).

2.9.1. Subjects

Ten caucasian male or female volunteers, aging from 20 to 52 years old were randomly entered in the study. After being fully informed, each volunteer signed a document stating they had been informed of their agreement, declaring they had understood the study design, the potential risks, and that they were taking part in the study of their own free will. The volunteers agreed to refrain from using body care cosmetics or moisturizers on the flexor side of the left forearm for one week prior to the measurements. Known previous irritation episodes or allergies to any type of topical formulations, soaps, or surfactants were considered a reason for exclusion. The exclusion criteria were: (a) Three days before the study, having used cosmetics for skin care, humectants, etc., on the arms, (b) known previous skin irritation episodes or allergies to any type of topical formulation (pharmaceutical or cosmetical), and (c) virgin TEWL values outside 4 to 8 g/h/m2 were discarded.

2.9.2. Test Procedure

The individuals stayed in the test room at least 20 min prior to the measurements (room conditions 22 ± 3 °C and 55% relative humidity). Their skin temperature was measured with a Skin Thermometer® ST500 (Courage and Khazaka, Electronic GmbH, Köln, Germany). Circles of 7 cm in diameter were drawn on the left volar forearm of the volunteer, which is the skin area for the application of each formulation. Then, a measurement of each parameter was taken (t0). Afterwards, the formulations were applied in the middle of each circle by means of a soft circular movement with the thumb. One hour after the application of the formulations (t1), a similar set of measurements was obtained at the same skin area. The effects were expressed as mean values of individual differences (t1t0) of each parameter and formulation.

2.9.3. Skin Parameters

The following parameters were tested:
-
Skin elasticity. The effect of the formulation on the elasticity of the upper skin layers was tested with a Cutometer® MPA 580 (Courage and Khazaka, Electronic GmbH, Köln, Germany). This measurement generates a negative pressure, drawing the skin into a probe that leads to a vertical deformation. When the negative pressure is switched off, the skin recovery is characterized [31,32] in terms of skin biomechanical properties.
-
Corneum stratum hydration (SCH). It was performed with a Corneometer® 825 (Courage and Khazaka, Electronic GmbH, Köln, Germany). This measures the capacitance variation of the dielectric properties of epidermic stratum corneum due to changes in skin hydration.
-
Transepidermal Water Loss (TEWL). The retrograde water permeation through skin was measured with a Tewameter® TM 300 (Courage and Khazaka, Electronic GmbH, Köln, Germany). This measures the vapor density gradient across the skin combining temperature and relative humidity sensors located in a hollow cylinder applied on the skin surface.

2.9.4. Statistical Analyses

Individual differences were analyzed with a one-way Anova and a Bonferroni post-test to impute the differences, using graphpad Prism® V. 5.00 (GraphPad Software Inc., San Diego, CA, USA).

3. Results

3.1. Physicochemical Properties and Stability

The prepared simple W1/O emulsions appeared as white, odorless, greasy, homogeneous and with a high consistency. The W1/O/W2 multiple emulsions appeared as white, odorless and homogenous without any signs of precipitation of drug, phase separation, creaming, or any visual change of appearance in any of the samples kept at 5 and 30 °C. At 40 °C, MEs became unstable after 10 weeks storage.
The multi-compartmental structure of both MEs was confirmed by optical microscopy. Multiglobules containing many small droplets in the internal phase are shown in Figure 1. Internal droplets (primary W1/O emulsion) ranged 1 to 5 μm. The diameters of multi-globules ranged between 10–50 μm and were later confirmed by laser diffractometry. An optical microphotograph of BFZ-CF is available in Figure 1. Droplets of O/W BFZ-CF ranged about 5 μm, similarly to the primary W1/O emulsion of the MEs.
Concerning pH and conductivity, MEs remained stable during six months until 30 °C. Conductivities at 24 h and after six months are summarized in Table 2. Statistically significant differences (p < 0.05) were found between the results of both MEs, but not for their differences between zero and six months. No significant pH variations were detected along the stability study.
Laser diffractometry confirmed the results of optical microscopy. Droplet size (D[4,3]) resulted to be 65.2 ± 0.8, 68.5 ± 1.0 and 29.8 ± 9.8µm for JMLP01B0, JMLP01BT, and BFZ-CF, respectively, as also detailed in Table 2.
A combined plot of particle size distributions of JMLP01B0, JMLP01BT (after 24 h of preparation) and BFZ-CF is shown in Figure 2.
The long-term physical stability of the MEs, predicted with the Turbiscan Lab® Expert, used only backscattered light (BS); this was because the formulations were not transparent. No relevant variations of the droplet volume fraction (migration) or mean size (coalescence) were observed (Figure 3).

3.2. Rheological Properties

3.2.1. Rotational Test

The results of steady-state rheology as a function of shear rate are summarized in Table 2. The formulations exhibited non-Newtonian pseudoplastic flow and shear thinning behavior with a consistent decrease in viscosity with increasing shear rate from 0 to 100 s−1. The Cross model showed the best fitting to both ascending and descending stretches. The flow curves (Figure 4) indicated certain apparent thixotropy as the rheograms displayed a moderate hysteresis loop with the downward curve below the upward curve.
Viscosity of BFZ-CF (at 10 and 100 s−1) was markedly higher than MEs values (Table 3 and Figure 5).

3.2.2. Oscillatory Test

Measurements were performed to find the critical stress. As shown in Figure 6, the critical stress was found at stresses below approximately 2 Pa. Following on from these results, a constant shear stress of 1 Pa (50% of the critical value) was selected to perform the frequency sweep tests. Results for both ME formulations revealed, in this frequency range, a prevalence of the elastic over the viscous behavior (G’>G’’), indicated in what is known as the “elastic plateau” and is shown in Table 3 and Figure 6. Additionally, a higher viscosity was observed for JMLP01BT, confirming previous rotational results.
Throughout the stability study, viscosity tended to be lower at end-times (see Table 3). The differences were not statistically significant (p > 0.05). All the cases were at room temperature.

3.3. Drug Content

Formulations JMLP01B0 and JMLP01BT (after 30 days of preparation) contained 90.6 ± 1.8% and 95.8 ± 2.2%, respectively, of the nominal drug content in the oil phase. BFZ-CF contained 94.1 ± 1.7% of the declared content of BFZ (mean, SD).

3.4. Release Test

Figure 7 outlines the BFZ release profiles (mean and SD) for each formulation and a methanolic solution reference. After 5 h, mean percentages (%Q5) were 41.0% and 48.8% for JMLP01B0 and JMLP01BT (asymptotic), respectively.
Based on AIC values, individual profiles were best described with the Higuchi equation in all cases except for the methanolic solution (Weibull). Higuchi release parameters are summarized in Table 4. In the case of the methanolic reference, the sigmoidicity of Weibull function is β = 9.6 ± 5.5 with a td 0.161 ± 0.009 h, as expected from a liquid solution. Based on ANOVA test, Bonferroni’s multiple comparison tests showed differences between JMLP01B0 and the BFZ-CF (p < 0.05). In other cases, the highest efficiency [27] was achieved with the ME containing polysorbate, although differences were not statistically significant between the three formulations (p > 0.05).

3.5. Skin Permeation

Permeation parameters and penetration values are summarized in Table 3. Statistically relevant differences (p < 0.05) were found between JMLP01BT and the other two formulations. The mean skin permeation profiles of BFZ of the three formulations are plotted in Figure 8.

3.6. Histological Analysis

Epidermis treated with the commercial O/W cream (BFZ-CF) presented a slightly detached stratum corneum, followed by a normal stratum granulosum, stratum spinosum, stratum basale, and the dermis (Figure 9a). Skin treated with the MEs placebos, with or without polysorbate (Figure 9b,d), showed an intact stratum corneum with unaltered adjacent structures. The skin in presence of JMLP01BT (Figure 9c) showed a pattern similar to that of its respective ME placebo, with an intact stratum corneum, without abnormal histologic layers detachment. No sign of cellular alteration in stratum spinosum was observed.

3.7. Skin Integrity Assessment

The skin elasticity modified by JMLP01B0 was significantly lower than the values obtained with the other formulations (0.0521 ± 0.003 vs. 0.2232 ± 0.03 or 0.1605 ± 0.04 AU). Otherwise, no variations in Hydration or TEWL before (t0) and after (t1) forearm applications were demonstrated (see Figure 10 and Table 5). In fact, respective variations of TEWL (at the same skin temperature) for both MEs were similar.

4. Discussion

4.1. Physicochemical Properties and Stability

BFZ is practically insoluble in water (0.13 µg/mL, pH 7.4, 32 °C) and sparingly soluble in ethanol [33,34]. Its total dose in the MEs (10 mg/g) was fully dissolved in the lipophilic intermediate phase (17.5% of the total weight) [17,35,36]. The main component of this oil phase, LP, is a synthetic medium chain (C8-10) triacylglyceride which has better interfacial properties in comparison with classical oils. It is a good solvent for lipophilic active pharmaceutical ingredients and is associated with enhanced drug penetration. As described in several papers [37,38], LP stabilizes MEs either alone or in the presence of polysorbate 80. Such a mixture allows greater solubilization of hydrophobic model drugs, and since this oil phase limits the interaction between internal and external aqueous phases, it may influence the drug release pattern depending on the phase-distribution of the drug. As for the pH, the non-ionic surfactant increases slightly the pH of the MEs [39], does not react with ionic substances and has a good history of being safe in its use [4]. Considering the pKa of BFZ (pKa = 6.29) [40] and the pH of Tween® 80 in aqueous solution (pH = 7) [39], a neutral pH in the aqueous phase is convenient to assure about a 50% ionization of drug favoring its transference from the oil phase to the external phase in contact with the skin. Sorbitan monooleate (grade 80) was selected a priori to achieve the best interaction with LP triglyceride in the oil phase and to intensify the subsequent BFZ penetration in stratum corneum.
As for the conductimetric results, the transference of NaCl from the inner layer W1 to W2 is not relevant due to the low conductivity values obtained at the end of the stability study [41]. It suggests the preservation of the phase’s distribution during ageing. An increase in conductivity of W2 would suggest the migration of electrolytes [42] as a result of diffusion or droplet breaking. In our case, the phase’s distribution of this ME is demonstrated as stable.
Droplet size analysis was carried out to detect possible changes in the distribution of the overall size and alterations in the volume mean diameter [43]. However, laser diffractometry only informs about the size distribution of the external globules of MEs [44], considering the dispersed volume fraction as being a simple O/W emulsion. This method, if it were combined with other ones, could give more information and an overview of the behavior of a multiple emulsion [45]. Assuming that BFZ is dissolved in the oil phase, this measurement of stability can be accepted as release indicating. No significant variations (α = 0.05) of droplet size have been detected during storage. Based on D[4,3] and D[3,2], the maintenance of space-fill, interfacial area, and thus, drug transport can be assumed [46]. Within experimental uncertainty both diameters in optical microscopy remained invariant over the whole storage period for each emulsion, indicating the absence of external and internal coalescence for samples kept at room temperature. The analyses of physical destabilization were measured using an optical analyzer Turbiscan® Lab Expert to predict the instability phenomena within backscattering profiles by multiple light scattering. If the %∆BS profiles have a deviation of ≤ 2 %, it can be considered that there are no significant variations in droplet size. Variations ≥10% indicate instability risk [46]. The formulations in this study were milky and opaque, therefore only BS profiles were used to evaluate the physicochemical stability. As shown in Figure 3, no relevant variations of the droplet volume fraction (particle migration) or mean size (coalescence, flocculation) were observed on the MEs. Stability of the MEs can also be assumed to be larger than six months.

4.2. Rheological Properties

Rheological measurements have been used to assess the resulting skin application properties and also its stability [5,47,48]. As is shown in Figure 4a, the steady shear viscosity sharply decreases and there is an increase in shear rate from 0 to 100 s−1. The determination of the disturbance of the W1/O/W2 emulsion microstructure during the shearing process of rheological behaviour (Figure 4b) was evaluated by measuring the area of hysteresis loop (Pa/s) called “apparent thixotropy” [17,49,50]. The thixotropic area provides qualitative information about flow time dependence [51]. In fact, the shearing process could induce a phase inversion of the multiple emulsion to a W/O emulsion type in a semi-solid state and negative thixotropic patterns when they were observed. In the shearing conditions of the present study, the developed multiple emulsions showed a non-negative thixotropic pattern, indicating that the experimental conditions did not lead to phase-inverted emulsion. Nonetheless, this behavior cannot be described as real thixotropy since in multiple emulsions examined by a cone-and-plate viscometer, it is difficult to recover initial structure after being subjected to strong shearing”. In this sense, these flow curves (Figure 4b) indicate certain apparent thixotropy as the rheograms displayed a hysteresis loop.
The best kinetic model (Cross-kinetic equation) provides a general model for pseudoplastic materials and suggests non-Newtonian shear-thinning flow under shear. It demonstrates a good extensibility and spreadability on the skin surface for the MEs. El-Hadidy et al. [40] obtained similar pseudoplasticity of microemulsions containing polysorbate 80. In our case, the presence of polysorbate 80 in the ME has significantly increased viscosity and apparent thixotropy, which are profitable effects in that they increase the physical retention of the formulation [52]. This fact, in line with the effect of carbomer, modulates the consistency of the formulations working towards an optimal application over large skin areas, as other authors have shown [53]. Conversely, the commercial formulation exhibits a greater viscosity (Table 2, Figure 5) requiring a greater apparent thixotropic effect for a proper extension on skin.
The oscillatory sweep stress and frequency sweep tests have been used, basically, to find the storage modulus (G’) that measures the ability of the system to store recoverable elastic energy, and the loss modulus (G’’) that reveals the dynamic viscosity associated with unrecoverable viscosity loss. As shown in Figure 6 and Table 2, all formulations revealed a prevalence of the elastic over the viscous behavior (G’ > G’’). Jiao et al. [54] described these same properties in their formulations containing Tween® 80, remarking its fluidification, as observed with JMLP01B0 and JMLP01BT. It is presumed that this elasticity favors the structural stability and resistance to external forces for longer periods of time [55].

4.3. Release and Skin Permeation

The experimental setup satisfied the required experimental standards as follows: spectrophotometric and chromatographic quantifications of BFZ were previously validated with acceptable linearity, precision, accuracy and recovery. Sink conditions in diffusion cells were guaranteed in all cases [56,57]. Mesh size of nylon membranes was clearly lower than droplet sizes, so one can be assured that the observed profiles are only due to the released drug and skin integrity before experiments [58] was assessed with TEWL measurements.
All release profiles have been described with the Higuchi function (square root), suggesting a depletion-limited release process and reaching asymptotic values rapidly after 5 h (Figure 7). It confirms that the oil intermediate phase acts as a drug reservoir-like with prolonged release towards the external aqueous phase and achieving a pseudo infinite-dose release model.
As initially enunciated, the objective of formulating a ME was the achievement of a lower skin permeation level and higher skin retention rather than a commercial O/W formulation. Skin permeation levels were ca. 10−3 times lower than the corresponding release levels. This difference is explained by the high lipophilicity of imidazolic drugs and their consequent low skin permeability. In fact, the formulation of drugs with limited permeability and solubility (BCS, class IV) offers a specific opportunity for a well-developed topical administration [59] due to their poor systemic absorption in animals and humans [47,48]. In our case, resultant values of BFZ permeation flux (Table 3) are slightly lower than those obtained by Hashiguchi et al. [60] with a BFZ lipophilic suspension through mouse skin (0.3386 µg/cm2/h)).
Currently, the treatment of Tinea versicolor and dermatophytosis [61] with Bifonazole requires its application once every 24 h for at least three weeks [1,2,62,63,64]. Given that these fungi commonly affect the epidermal cornified layers and the dermis [1], the intradermal retention of bifonazole appears as an indicator of its efficient delivery [62]. It is known that imidazolic drugs are prone to reversible bind to queratin in stratum corneum in direct relationship with their lipophilicity [60] prolonging its efficacy. Additionally, the presence of high concentrations of drug with specific surfactants in the formulation contributes to enhance the antimycotic effect [65]. In our case, drug accumulation levels in epidermis achieved with both MEs (Table 3) are higher than those with the commercial cream (0.4 µg/cm2), suggesting the possibility of lengthening the currently used 24 h administration interval up to an advantageous 48 h interval.
Even in the lower layers of the epidermis, the BFZ accumulated is several times higher than the in vitro minimum inhibitory concentrations for dermatophytes, especially along the hair follicles as the preferred route of penetration [66]. This marked intracutaneous distribution has been described, for example, with tioconazole in guinea pig skin, detecting (after 2 h) around 300 µg/g in stratum corneum, but only 25 µg/g in epidermis [67]. These results confirm that if BFZ is relicted in the lipophilic skin regions then its apparent permeation through the dermis can be reduced dramatically [68,69] favoring its antifungal activity in the epidermic tissue.
The efficiency of the drug release profiles from the MEs (cfr.: Equation 2) is highest for the polysorbate-containing formulation (JMLP01BT) and minimum for the commercial formulation (BFZ-CF), which is also directly proportional to the respective skin accumulation levels.
The drug penetration levels from the ME without polysorbate 80 (2.9 µg/cm2) are similar to those obtained by other authors [60] with a BFZ suspension (2.4 µg/cm2) in mouse skin. As for the polysorbate ME formulation, the skin accumulation levels are statistically higher (16.60 µg/cm2) than all the others (Table 3). This formulation dramatically increased the accumulation of BFZ inside the epidermic layers (400 µm thickness). Other authors [70], who directly related the amounts of drug in stratum corneum with the respective uptake of the excipients, have described a similar effect. In fact, the similarity among aliphatic molecules is sometimes associated with a potentiation of the enhancer activities [71]. For example, the effect of polysorbate 80 over the permeation of lipophilic drugs (logP above 4.15) through rat mucosal membranes [72] is related to an increase in diffusion coefficient, and thus, the drug transport, similarly with econazole [16]. A significant increase in skin deposition means improved cutaneous drug availability, and thus, the feasibility of an enlarged interval of administration.

4.4. Histological Analyses

While polysorbate 80 can cause anaphylactic reactions (non-immunological) in individuals prone to them, sensitization in selected patients with contact dermatitis is quite rare [73,74]. As a preliminary test for in vivo evaluation, residual skin samples of the in vitro tests were stained and observed in order to for any eventual histological alteration due to this substance. The representative microscopical images (Figure 9) illustrate that neither the presence of polysorbate alters the histological skin structure (placebo with vs. placebo without) and neither does the drug (Figure 9b,c). A for the other the O/W emulsion, it seems that the stratum corneum treated with the commercial emulsion clearly becomes detached from the viable epidermis.

4.5. Skin Integrity

In vivo results of stratum corneum hydration and transepidermal water permeability suggest a tendency to increase hydration and to decrease TEWL Figure 10, Table 4), probably due to an occlusive effect of the ME. This humectant effect is used [75] as an indicator of skin healing in antimycotic treatments [76]. In addition, TEWL is the best indicator of skin integrity [77] and similar effects are observed in both MEs. Thus, the adverse effects of polysorbate 80 over the skin integrity can be considered as negligible.
Regarding elasticity variations (t1t0), significancy (P < 0.05) was observed between the effect of JMLPB0 and the polysorbate-containing ME and the commercial formulation (Table 4). This increase can be considered as non-irritancy relevant because it is not univocally attributable to the presence of polysorbate 80.
In aggregate, the combination of amphoteric and non-ionic surfactants in the ME has not resulted in a substantial alteration of skin integrity. Polysorbates are particularly useful to solubilize hydrophobic organic compounds enhancing its transit through porous substrates [78]. Its use in the formulation of a multiple emulsion has optimized the delivery of BFZ, a highly lipophilic drug, to the epidermic stratum corneum.

5. Conclusions

Bifonazole multiple emulsion with polysorbate achieves a clearly lower skin permeation and a higher skin retention than the conventional cream formulation. The formulation displayed an increase in pH, viscosity and viscoelasticity and a decrease of spreadibility constituting an interesting alternative for dermal application. The polysorbate-containing multiple emulsion has efficiently enhanced the ex vivo penetration of BFZ without altering the in vivo skin integrity and, thus, it optimizes the delivery of bifonazole, a highly lipophilic drug, towards the epidermic stratum corneum.
It has been shown that this type of formulation is likely to be tolerable in terms of biophysical skin-integrity and in terms of histopathological signals. In accordance with this, the proposed antifungal multiple emulsion is the preferred choice for the epidermal delivery of bifonazole with a prolonged interval administration and it is a promising alternative formulation for this imidazolic drug.

Author Contributions

Conceptualization, J.S.-C., A.C.-C., B.C.-N., A.B.-M.; investigation, J.S.-C., J.Z.-L., M.J.R.-L., E.B.; writing-original draft preparation, J.Z.-L., E.B.; writing review and editing, A.C.-C., L.H.-B., J.S.-C., A.B.-M.; funding acquisition L.H.-B.; J.S.-C.

Acknowledgments

There are no specific funding sources or universitary fellowship for this experimental work. Article processing charges were partially supported by University of Barcelona.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Quatresooz, P.; Piérard-Franchimont, C.; Arrese, J.E.; Piérard, G.E. Clinicopathologic presentations of dermatomycoses in cancer patients. J. Eur. Acad. Derm. Venereol. 2008, 22, 907–917. [Google Scholar] [CrossRef] [PubMed]
  2. SmPC Canesmycospor®. Agencia Española de Medicamentos y Productos Sanitarios. Bayer Hispania SL. Revisited on 02/2015. Available online: https://cima.aemps.es/cima/publico/lista.html (accessed on 29 December 2018).
  3. Schaefer, H.; Stüttgen, G. Absolute concentrations of an antimycotic agent, Econazole, in the human skin after local application. Arzneim. Forsch. Drug Res. 1976, 26, 432–435. [Google Scholar]
  4. Seifriz, W. Studies in emulsion III. Double reversal of oil emulsions occasioned by the same electrolyte. J. Phys. Chem. 1925, 29, 738–749. [Google Scholar] [CrossRef]
  5. Grossiord, J.L.; Seiller, M. Applications. In Multiple Emulsion: Structure, Properties and Applications; Grossiord, J.L., Seiller, M., Eds.; Editions de Santé: Paris, France, 1998; p. 169. [Google Scholar]
  6. Nokhodchi, A.; Shokri, J.; Dashbolaghi, A.; Hassan-Zadeh, D.; Ghafourian, T.; Barzegar-Jalali, M. The enhancement effect of surfactants on the penetration of lorazepam through rat skin. Int. J. Pharm. 2003, 250, 359–369. [Google Scholar] [CrossRef]
  7. Firooz, A.; Nafisi, S.; Maibach, H.I. Novel drug delivery strategies for improving econazole antifungal action. Int. J. Pharm. 2015, 495, 599–607. [Google Scholar] [CrossRef] [PubMed]
  8. Garti, N. Double emulsions-scope, limitations and new achievements. Colloids Surf. 1997, 123–124, 233–246. [Google Scholar] [CrossRef]
  9. Kaur, G.; Mehta, S.K. Developments of Polysorbate (Tween) based microemulsions: Preclinical drug delivery, toxicity and antimicrobial applications. Int. J. Pharm. 2017, 529, 134–160. [Google Scholar] [CrossRef]
  10. Silva, A.; Grossiord, J.L.; Puiseux, F.; Seiller, M. Insulin in W/O/W multiple emulsions: Preparation, characterization and determination of stability towards proteases in vitro. J. Microencap. 1997, 14, 311–319. [Google Scholar]
  11. McClements, D.J. Advances in fabrication of emulsions with enhanced functionality using structural design principles. Curr. Opin. Colloid Interface Sci. 2012, 17, 235–245. [Google Scholar] [CrossRef]
  12. Tang, S.Y.; Sivakumar, M.; Ng, A.M.H.; Shridharan, P. Anti-inflammatory and analgesic activity of novel oral aspirin-loaded nanoemulsion and nano multiple emulsion formulations generated using ultrasound cavitation. Int. J. Pharm. 2012, 430, 299–306. [Google Scholar] [CrossRef]
  13. Omotosho, J.A.; Whateley, A.T.; Florence, A.T.; Bell, G. Release of cytotoxic agents from multiple W/O/W Emulsions. J. Pharm. Pharmacol. 1987, 39, 38P. [Google Scholar]
  14. Nakhare, S.; Vyas, S.P. Prolonged release multiple emulsion based system bearing rifampicin: In vitro characterisation. Drug Dev. Ind. Pharm. 1995, 21, 869–878. [Google Scholar] [CrossRef]
  15. Li, M.; Du, C.; Guo, N.; Teng, Y.; Meng, X.; Sun, H.; Li, S.; Yu, P.; Galons, H. Composition design and medical application of liposomes. Eur. J. Med. Chem. 2019, 164, 640–653. [Google Scholar] [CrossRef] [PubMed]
  16. Bachhav, Y.G.; Mondon, K.; Kalia, Y.N.; Gurny, R.; Moller, M. Novel micelle formulation to increase cutaneous bioavailability of azole antifungals. J. Control. Release 2011, 153, 126–132. [Google Scholar] [CrossRef] [PubMed]
  17. Suñer, J.; Calpena, A.C.; Clares, B.; Cañadas, C.; Halbaut, L. Development of Clotrimazole Multiple W/O/W Emulsions as Vehicles for Drug Delivery: Effects of Additives on Emulsion Stability. AAPS PharmSciTech 2017, 18, 539–550. [Google Scholar] [CrossRef] [PubMed]
  18. Matsumoto, S.; Kita, Y.; Yonezawa, D. An attempt at preparing water-in-oil-in water multiple-phase emulsions. J. Colloid Interface. Sci. 1976, 57, 353–361. [Google Scholar] [CrossRef]
  19. Malvern Intruments Ltd. Mastersizer 2000; User Manual. Worcestershire, UK, 2007. Available online: https://www.malvernpanalytical.com/en (accessed on 29 December 2018).
  20. Park, E.K.; Song, K.W. Rheological evaluation of petroleum jelly as a base material in ointment and cream formulations, steady shear flow behaviour. Arch. Pharm. Res. 2010, 33, 141–150. [Google Scholar] [CrossRef]
  21. Popović, G.; Čakar, M.; Agbaba, D. Determination of bifonazole in creams containing methyl- and propyl p-hydroxybenzoate by derivative spectrophotometric method. J. Pharm. Biomed. Anal. 2003, 33, 131–136. [Google Scholar] [CrossRef]
  22. European Medicines Agency (EMEA). Guideline on Validation of Bioanalytical Methods, Document Reference EMEA/CHMP/EWP/192217/2009; EMEA: London, UK, 2011. [Google Scholar]
  23. Čudina, O.A.; Čomor, M.I.; Janković, I.A. Simultaneous Determination of Bifonazole and Benzyl Alcohol in Pharmaceutical Formulations by Reverse-Phase HPLC. Chromatographia 2005, 61, 415–418. [Google Scholar] [CrossRef]
  24. Keshary, P.R.; Huang, Y.C.; Chien, Y.W. Mechanism of transdermal controlled nitroglycerin administration (III), Control of skin permeation rate and optimization. Drug Dev. Ind. Pharm. 1985, 11, 1213–1253. [Google Scholar] [CrossRef]
  25. Costa, P.; Sousa, J.M. Modeling and comparison of dissolution profiles. Eur. J. Pharm. Sci. 2001, 13, 123–133. [Google Scholar] [CrossRef]
  26. Yamaoka, K.; Nakagawa, T.; Uno, T. Application of Akaike’s information criterion (AIC) in the evaluation of linear pharmacokinetic equations. J. Pharmacokinet. Biopharm. 1978, 6, 165–175. [Google Scholar] [CrossRef] [PubMed]
  27. Khan, K.A. The concept of dissolution efficiency. J. Pharm. Pharmacol. 1975, 27, 48–49. [Google Scholar] [CrossRef]
  28. Oestmann, E.; Lavrijsen, A.P.; Hermans, J.; Ponec, M. Skin barrier function in healthy volunteers as assessed by transepidermal water loss and vascular response to hexyl nicotinate, intra- and inter-individual variability. Br. J. Dermatol. 1993, 128, 130–136. [Google Scholar] [CrossRef] [PubMed]
  29. Torres, E.; Suner-Carbó, J.; Aróztegui, M.; Halbaut, L.; Barbé, C. Propuesta de protocolo de analisis sensorial en productos semisólidos: Cremas i geles. Not. Cosmet. Perfum. 2001, 260, 5–11. [Google Scholar]
  30. World Medical Association. World medical association declaration of Helsinki: Ethical principles for medical research involving human subjects. JAMA 2013, 310, 2191–2194. [Google Scholar] [CrossRef]
  31. Bonaparte, J.P.; Ellis, D.; Chung, J. The effect of probe to skin contact force on Cutometer MPA 580 measurements. J. Med. Eng. Technol. 2013, 37, 208–212. [Google Scholar] [CrossRef]
  32. Neto, P.; Ferreira, M.; Bahia, F.; Costa, P. Improvement of the methods for skin mechanical properties evaluation through correlation between different techniques and factor analysis. Skin Res. Technol. 2013, 19, 405–416. [Google Scholar] [CrossRef]
  33. Council of Europe. The European Pharmacopoeia Commission. In European Pharmacopoeia, 8th ed.; Bifonazole monograph 04/2012:1395; Council of Europe: Strasbourg Cedex, France, 2016; p. 1670. [Google Scholar]
  34. Elsayed, M.A. Development of topical therapeutics for management of onychomycosis and other nail disorders, A pharmaceutical perspective. J. Control. Release 2015, 199, 132–144. [Google Scholar] [CrossRef]
  35. Suñer, J.; Boix, A.; Halbaut, L.; Velázquez, N.; Zamarbide, J.; Bozal-de-Febrer, N.; Calpena, A.C. Skin permeation of Econazole nitrate formulated in an enhanced hydrophilic multiple emulsion. Mycoses 2017, 60, 166–177. [Google Scholar] [CrossRef]
  36. Hameyer, P.; Jenni, K.R. Emulsifiers for multiple emulsions. Cosm. Toil. 1996, 111, 39–48. [Google Scholar]
  37. Devani, M.; Ashford, M.; Craig, D.Q.M. The emulsification and solubilisation properties of polyglycolysed oils in self-emulsifying formulations. J. Pharm. Pharmacol. 2004, 56, 307–316. [Google Scholar] [CrossRef] [PubMed]
  38. Yang, J.H.; Kim, Y.I.; Kim, K.M. Preparation and evaluation of aceclofenac microemulsion for transdermal delivery system. Arch. Pharm. Res. 2002, 25, 534–540. [Google Scholar] [CrossRef] [PubMed]
  39. Shalaeva, M.; Kenseth, J.; Lombard, F.; Batin, A. Measurement of dissociation constants (pKa Values) of organic compounds by multiplexed capillary electrophoresis using aqueous and cosolvent buffers. J. Pharm. Sci. 2008, 97, 2581–2606. [Google Scholar] [CrossRef] [PubMed]
  40. El-Hadidy, G.N.; Ibrahim, H.K.; Mohamed, M.I.; El-Milligi, M.F. Microemulsions as vehicles for topical administration of voriconazole, formulation and in vitro evaluation. Drug Dev. Ind. Pharm. 2012, 38, 64–72. [Google Scholar] [CrossRef] [PubMed]
  41. Muguet, V.; Seiller, M.; Barratt, G.; Ozer, O.; Marty, J.P.; Grossiord, J.L. Formulation of shear rate sensitive multiple emulsions. J. Control. Release 2001, 70, 37–49. [Google Scholar] [CrossRef]
  42. Tedajo, G.M.; Seiller, M.; Prognon, P.; Grossiord, J.L. pH compartmented w/o/w multiple emulsion, a diffusion study. J. Control. Release 2001, 75, 45–53. [Google Scholar] [CrossRef]
  43. Marku, D.; Wahlgren, M.; Rayner, M.; Sjöö, M.; Timgren, A. Characterization of starch Pickering emulsions for potential applications in topical formulations. Int. J. Pharm. 2012, 428, 1–7. [Google Scholar] [CrossRef]
  44. Dickinson, E.; Evison, J.; Owusu, R.K. Preparation of fine protein-stabilized water-in-oil-in water emulsions. Food Hydrocolloids 1991, 5, 481–485. [Google Scholar] [CrossRef]
  45. Robbe-Tomine, L.; Le Hen-Ferrenbach, C.; Pouget, T.; Tranchant, J.F. Multiple emulsions visualization methods and particle size analysis. In Multiple Emulsion: Structure, Properties and Applications; Grossiord, J.L., Seiller, M., Eds.; Editions de Sante: Paris, France, 1998. [Google Scholar]
  46. Celia, C.; Trapasso, E.; Cosco, D.; Paolino, D.; Fresta, M. Turbiscan® Lab Expert analysis of the stability of ethosomes and ultradeformable liposomes containing a bilayer fluidizing agent. Colloids Surf. B Biointerfaces 2009, 72, 155–160. [Google Scholar] [CrossRef]
  47. Vasiljevic, D.; Vuleta, G.; Primorac, M. The characterization of the semi-solid W/O/W emulsions with low concentrations of the primary polymeric emulsifier. Int. J. Cosm. Sci. 2005, 27, 81–87. [Google Scholar] [CrossRef] [PubMed]
  48. Krishnaiah, Y.S.; Xu, X.; Rahman, Z.; Yang, Y.; Katragadda, U.; Lionberger, R.; Peters, J.R.; Uhl, K.; Khan, M.A. Development of performance matrix for generic product equivalence of acyclovir topical creams. Int. J. Pharm. 2014, 475, 110–122. [Google Scholar] [CrossRef] [PubMed]
  49. Soriano-Ruiz, J.L.; Suñer-Carbó, J.; Calpena-Campmany, A.C.; Bozal-de Febrer, N.; Halbaut-Bellowa, L.; Boix-Montañés, A.; Souto, E.B.; Clares-Naveros, B. Clotrimazole multiple W/O/W emulsion as anticandidal agent: Characterization and evaluation on skin and mucosae. Colloids Surf. B Biointerfaces 2019, 175, 166–174. [Google Scholar] [CrossRef] [PubMed]
  50. Kovács, A.; Erős, I.; Csóka, I. Optimization and development of stable w/o/w cosmetic multiple emulsions by means of the Quality by Design approach. Int. J. Cosmet. Sci. 2016, 38, 128–138. [Google Scholar] [CrossRef] [PubMed]
  51. Kawashima, Y.; Hino, T.; Takeuchi, H.; Niwa, T.; Horibe, K. Rheological study of w/o/w emulsion by a cone-and-plate viscometer: Negative thixotropy and shear-induced phase inversion. Int. J. Pharm. 1991, 72, 65–77. [Google Scholar] [CrossRef]
  52. Silva, A.C.; Amaral, M.H.; Gonzalez-Mira, E.; Santos, D.; Ferreira, D. Solid lipid nanoparticles (SLN) based hydrogels as potential carriers for oral transmucosal delivery of Risperidone, Preparation and characterization studies. Colloids Surf. B Biointerfaces 2012, 93, 241–248. [Google Scholar] [CrossRef] [PubMed]
  53. Valenta, C.; Schultz, K. Influence of carrageenan on the rheology and skin permeation of microemulsion formulations. J. Control. Release 2004, 95, 257–265. [Google Scholar] [CrossRef] [PubMed]
  54. Jiao, J.; Rhodes, D.G.; Burgess, D.J. Multiple emulsion stability, Pressure balance and Interfacial Film strength. J. Colloid Interface Sci. 2002, 250, 444–450. [Google Scholar] [CrossRef] [PubMed]
  55. Korhonen, M.; Niskanen, H.; Kiesvaara, J.; Yliruusi, J. Determination of optimal combination of surfactants in creams using rheology measurements. Int. J. Pharm. 2000, 197, 143–151. [Google Scholar] [CrossRef]
  56. Cañadas-Enrich, C.; Abrego, G.; Alvarado, H.L.; Calpena-Campmany, A.C.; Boix-Montañes, A. Pranoprofen quantification in ex vivo corneal and scleral permeation samples: Analytical validation. J. Pharm. Biomed. Anal. 2018, 160, 109–118. [Google Scholar] [CrossRef]
  57. European Medicines Agency. Draft Guideline on Quality and Equivalence of Topical Products; CHMP/QWP/708282/2018; European Medicines Agency: London, UK, 2018. [Google Scholar]
  58. Klang, V.; Schwarz, J.C.; Haberfeld, S.; Xiao, P.; Wirth, M.; Valenta, C. Skin integrity testing and monitoring of in vitro tape stripping by capacitance-based sensor imaging. Skin Res. Technol. 2013, 19, e259–e272. [Google Scholar] [CrossRef] [PubMed]
  59. Fromtling, R.A. Overview of medically important antifungal azole derivatives. Clin. Microbiol. Rev. 1988, 1, 187–217. [Google Scholar] [CrossRef] [PubMed]
  60. Hashiguchi, T.; Kodama, A.; Ryu, A.; Otagiri, M. Retention capacity of topical imidazole antifungal agents in the skin. Int. J. Pharm. 1998, 161, 195–204. [Google Scholar] [CrossRef]
  61. Fromtling, R.A. Bifonazole (MycosporRp) an update. Drugs Today 1985, 21, 401–404. [Google Scholar]
  62. Plempel, M.; Regel, E.; Büchel, K.H. Antimycotic efficacy of Bifonazole in vitro and in vivo. Arzneimittel-Forschung 1983, 33, 517–524. [Google Scholar]
  63. Polak, A. Antifungal activity of four antifungal drugs in the cutaneous retention time test. Sabouraudia J. Med. Vet. Mycol. 1984, 22, 501–503. [Google Scholar] [CrossRef]
  64. Ritter, W.; Siefert, H.M. Biological disposition and percutaneous absorption of bifonazole in animals and man. In Recent Trends in the Discovery, Development and Evaluation of Antifungal Agents; Fromtling, R.A., Ed.; JR Prous Publishers: Barcelona, Spain, 1987; pp. 383–405. [Google Scholar]
  65. Beggs, W.H.; Hughes, C.E. Exploitation of the direct cell damaging action of antifungal azoles. Diagn. Microb. Infect. Dis. 1987, 6, 1–3. [Google Scholar] [CrossRef]
  66. Patzschke, K.; Ritter, W.; Siefert, H.M.; Weber, H.; Wegner, L.A. Pharmacokietic studies following systemic and topical administration of [14C]-Bifonazole in Man. Arzneimittel-Forschung 1983, 33, 745–750. [Google Scholar]
  67. Sobue, S.; Sekiguchi, K. Difference in Percutaneous Absorption and Intracutaneous Distribution in Guinea Pigs among Topical Antifungal Drugs (Tioconazole Solution, Tioconazole Cream, Miconazole Nitrate Solution and Bifonazole Solution). Biol. Pharm. Bull. 2004, 27, 1428–1432. [Google Scholar] [CrossRef]
  68. Klang, V.; Haberfeld, S.; Hartl, A.; Valenta, C. Effect of γ-cyclodextrin on the in vitro skin permeation of a steroidal drug from nanoemulsions, Impact of experimental setup. Int. J. Pharm. 2012, 423, 535–542. [Google Scholar] [CrossRef]
  69. Baker, E.J.; Hadgraft, J. In vitro percutaneous absorption of arildone, a highly lipophilic drug, and the apparent no-effect of the penetration enhancer Azone in excised human skin. Pharm. Res. 1995, 12, 993–997. [Google Scholar] [CrossRef] [PubMed]
  70. Mohammed, D.; Matts, P.J.; Hadgraft, J.; Lane, M.E. In Vitro–In Vivo Correlation in Skin Permeation. Pharm. Res. 2014, 31, 394–400. [Google Scholar] [CrossRef] [PubMed]
  71. Lopez-Castellano, A.; Cortell-Ivars, C.; Lopez-Carballo, G.; Herraez-Dominguez, M. The influence of Span® 20 on stratum corneum lipids in Langmuir monolayers, comparison with Azone®. Int. J. Pharm. 2000, 203, 245–253. [Google Scholar] [CrossRef]
  72. Salimi, A.; Hedayatipour, N.; Moghimipour, E. The Effect of Various Vehicles on the Naproxen Permeability through Rat Skin, A Mechanistic Study by DSC and FT-IR Techniques. Adv. Pharm. Bull. E 2016, 6, 9–16. [Google Scholar] [CrossRef]
  73. Tosti, A.; Guerra, L.; Morelli, R.; Barda, F. Prevalence and sources of sensitization to emulsifiers, a clinical study. Contact Dermat. 1990, 23, 68–72. [Google Scholar] [CrossRef]
  74. Coors, E.A.; Seybold, H.; Merk, H.F.; Mahler, V. Polysorbate 80 in medical products and nonimmunologic anaphylactoid reactions. Ann. Allergy Asthma Immunol. 2005, 95, 593–599. [Google Scholar] [CrossRef]
  75. Smits, J.; Weibel, M.; Herbst, N. Hydro-Gain®. Un Sistema humectante de origen vegetal que estimula la hidratación de la piel y fortalece la barrera lipídica. Not. Cosmet. Perfum. 2015, 342, 16–24. [Google Scholar]
  76. Nielsen, P.G. The Importance of the Vehicle in the Treatment of Dermatophytosis in Hereditary Palmo-Plantar Keratoderma. Mycoses 1984, 27, 227–230. [Google Scholar] [CrossRef]
  77. Nguyen, H.X.; Puri, A.; Banga, A.K. Methods to simulate rubbing of topical formulation for in vitro skin permeation studies. Int. J. Pharm. 2017, 519, 22–33. [Google Scholar] [CrossRef]
  78. Cheng, M.; Zeng, G.; Huang, D.; Yang, C.; Lai, C.; Zhang, C.; Liu, Y. Advantages and challenges of Tween 80 surfactant-enhanced technologies for the remediation of soils contaminated with hydrophobic organic compounds. Chem. Eng. J. 2017, 314, 98–113. [Google Scholar] [CrossRef]
Figure 1. Optical photomicrographs at 400× magnification of JMLP01B0 (a,c), JMLP01BT (b,d) and BFZ-CF (e) during the storage period (a,b: At 24 h of preparation; and c,d: At six months of preparation).
Figure 1. Optical photomicrographs at 400× magnification of JMLP01B0 (a,c), JMLP01BT (b,d) and BFZ-CF (e) during the storage period (a,b: At 24 h of preparation; and c,d: At six months of preparation).
Pharmaceutics 11 00066 g001
Figure 2. Droplet size distribution for JMLP01B0, JMLP01BT, and BFZ-CF as volume percentages.
Figure 2. Droplet size distribution for JMLP01B0, JMLP01BT, and BFZ-CF as volume percentages.
Pharmaceutics 11 00066 g002
Figure 3. Backscattering profiles (%∆BS, left ordinates) recorded at room temperature of JMLP01B0 (a) and JMLP01BT (b) at different times (0 to 24 h, right ordinates). Sample cell length (0 to 50 mm) is shown on abscissa.
Figure 3. Backscattering profiles (%∆BS, left ordinates) recorded at room temperature of JMLP01B0 (a) and JMLP01BT (b) at different times (0 to 24 h, right ordinates). Sample cell length (0 to 50 mm) is shown on abscissa.
Pharmaceutics 11 00066 g003
Figure 4. Viscosity (a) and flow (b) curves of MEs as functions of shear rate after 24 h of preparation.
Figure 4. Viscosity (a) and flow (b) curves of MEs as functions of shear rate after 24 h of preparation.
Pharmaceutics 11 00066 g004
Figure 5. Viscosity (blue line) and Flow curve (red line) of BFZ-CF as a function of shear rate.
Figure 5. Viscosity (blue line) and Flow curve (red line) of BFZ-CF as a function of shear rate.
Pharmaceutics 11 00066 g005
Figure 6. Storage modulus (G’), loss modulus (G’’) and phase angle (δ) of JMLP01B0 and JMLP01BT during a sweep stress test (top). Storage modulus (G’), loss modulus (G’’) and complex viscosity (η*) versus frequency during a frequency sweep test for JMLP01B0 and JMLP01BT (bottom) at 24 h of preparation.
Figure 6. Storage modulus (G’), loss modulus (G’’) and phase angle (δ) of JMLP01B0 and JMLP01BT during a sweep stress test (top). Storage modulus (G’), loss modulus (G’’) and complex viscosity (η*) versus frequency during a frequency sweep test for JMLP01B0 and JMLP01BT (bottom) at 24 h of preparation.
Pharmaceutics 11 00066 g006
Figure 7. Release profiles (mean values and standard deviations) of BFZ (squares JMP01BT, triangles up JMP01B0, triangles down BFZ-CF, rhombus 1% methanolic solution).
Figure 7. Release profiles (mean values and standard deviations) of BFZ (squares JMP01BT, triangles up JMP01B0, triangles down BFZ-CF, rhombus 1% methanolic solution).
Pharmaceutics 11 00066 g007
Figure 8. BFZ human skin permeated amounts (mean values and standard deviations) for the BFZ-CF (squares), JMLP01B0 (triangles down) and JMLP01BT (triangles up).
Figure 8. BFZ human skin permeated amounts (mean values and standard deviations) for the BFZ-CF (squares), JMLP01B0 (triangles down) and JMLP01BT (triangles up).
Pharmaceutics 11 00066 g008
Figure 9. Histological analysis of skin structure. Hematoxylin and eosin staining of skin treated with: (a) BFZ-CF, (b) ME reference free of polysorbate, (d) ME reference with polysorbate, and (c) JMLP01BT (400× magnification). Sc = stratum corneum, sg = stratum granulosum, ss = stratum spinosum, sb = stratum basale, d = dermis, and * = loss of stratum corneum.
Figure 9. Histological analysis of skin structure. Hematoxylin and eosin staining of skin treated with: (a) BFZ-CF, (b) ME reference free of polysorbate, (d) ME reference with polysorbate, and (c) JMLP01BT (400× magnification). Sc = stratum corneum, sg = stratum granulosum, ss = stratum spinosum, sb = stratum basale, d = dermis, and * = loss of stratum corneum.
Pharmaceutics 11 00066 g009
Figure 10. Stratum corneum Hydration (SCH) at t0 (a) and t1 (b), Transepidermal Water Loss (TEWL) at t0 (c) and at t1 (d). Variations are expressed as standard deviations.
Figure 10. Stratum corneum Hydration (SCH) at t0 (a) and t1 (b), Transepidermal Water Loss (TEWL) at t0 (c) and at t1 (d). Variations are expressed as standard deviations.
Pharmaceutics 11 00066 g010
Table 1. Formulation and percentage composition (w/w) of BFZ MEs without polysorbate (“–B0”) and with polysorbate, suffixed with “–BT”.
Table 1. Formulation and percentage composition (w/w) of BFZ MEs without polysorbate (“–B0”) and with polysorbate, suffixed with “–BT”.
ComponentsPercentage Composition (w/w)
JMLP01B0JMLP01BT
Oil phase (O)
Bifonazole1.001.00
Capric/caprylic triglyceride (Labrafac® Lipophile 1349)11.0011.00
Cetyl palmitate2.002.00
Cetyl dimethicone copolyol (Abil® EM 90)1.501.50
Sorbitan stearate (Span® 60)2.002.00
Internal aqueous phase (W1)
Sodium chloride0.250.25
Purified water at pH 6.632.2532.25
External aqueous phase (W2)
Carbomer (Tego® Carbomer 341)0.200.20
Cocamidopropyl betaine (Tego® Betaine F)0.700.70
Polysorbate 80 (Tween® 80)-1.00
Purified water at pH 6.649.1048.10
Table 2. pH, conductivity and particle size parameters by laser diffraction, for JMLP01B0, JMLP01BT and BFZ-CF at 24 h (first lines) and 6 months (second lines) (mean ± SD, n = 3).
Table 2. pH, conductivity and particle size parameters by laser diffraction, for JMLP01B0, JMLP01BT and BFZ-CF at 24 h (first lines) and 6 months (second lines) (mean ± SD, n = 3).
SamplepHConductivity (μm/s−1)Droplet Size Analysis (μm)
D[4,3]D[3,2]D[v,0.1]D[v,0.5]D[v,0.9]
JMLP01B05.80 ± 0.04
6.01 ± 0.01
309.5 ± 6.5
336.3 ± 7.1
65.2 ± 0.8
75.2 ± 0.7
9.7 ± 0.2
10.5 ± 0.7
3.4 ± 0.3
4.1 ± 0.7
35.8 ± 0.3
33.8 ± 0.5
172.9 ± 6.3
189.2 ± 7.1
JMLP01BT6.11 ± 0.02
6.30 ± 0.04
257.3 ± 3.3
289.0 ± 3.5
68.5 ± 1.0
60.8 ± 0.9
12.8 ± 0.1
12.1 ± 0.8
6.0 ± 0.2
6.3 ± 0.8
40.4 ± 5.6
49.5 ± 1.8
168.1 ± 3.4
165.02 ± 4.1
BFZ-CF6.09 ± 0.04
N.A.
199.1 ± 6.3
N.A.
29.8 ± 9.8
N.A.
5.8 ± 0.5
N.A.
1.8 ± 0.2
N.A.
18.0 ± 2.1
N.A.
59.0 ± 13.9
N.A.
(N.A. Not applies).
Table 3. Mean viscosity values (mPa.s) at 10 s−1 and 100 s−1 (mean ± SD, n = 3) and oscillatory test of storage modulus (G’), loss modulus (G’’), and complex viscosity (η*) at 0.01 s−1 and 10 s−1.
Table 3. Mean viscosity values (mPa.s) at 10 s−1 and 100 s−1 (mean ± SD, n = 3) and oscillatory test of storage modulus (G’), loss modulus (G’’), and complex viscosity (η*) at 0.01 s−1 and 10 s−1.
FormulationFrequency (s1)Viscosity (mPa.s)Oscillatory Measurements (at 24 h)
at 24 h180 daysG’(Pa)G’’(Pa)η*(Pa.s)
JMLP01B00.01N.A.N.A.62.0819.471035.00
101745.0 ± 16.71705.0 ± 15.785.0927.451.42
100237.5 ± 4.7205.8 ± 2.6N.A.N.A.N.A.
JMLP01BT0.01N.A.N.A.81.9930.291391
102176.0 ± 24.62093.5 ± 19.9167.6060.972.84
100295.8 ± 6.8261.2 ± 4.1N.A.N.A.N.A.
BFZ-CF0.01N.A.N.A.3736187566,530
1019,301.7 ± 99.719,240.5 ± 116.915,9404160262.20
1001901.3 ± 30.91963.0 ± 23.7N.A.N.A.N.A.
N.A.: Not Apply.
Table 4. Drug release parameters (K and efficiency), skin permeation flux and skin penetration of BFZ from both MEs and the commercial formulation (mean± SD). Statistical differences are indicated with superscript letters.
Table 4. Drug release parameters (K and efficiency), skin permeation flux and skin penetration of BFZ from both MEs and the commercial formulation (mean± SD). Statistical differences are indicated with superscript letters.
FormulationK HiguchiAUC05/(Q5·T)Permeation FluxSkin Penetration
(µg·h−1/2)(h)(µg/cm2h)(µg/g)(µg/cm2)
JMLP01B0539.3 ± 71.50.644 ± 0.0720.137 ± 0.008 b13.12 ± 1.24 b2.87 ± 0.36
JMLP01BT690.4 ± 97.8 a0.716 ± 0.0260.064 ± 0.012106.20 ± 8.7316.60 ± 7.09
BFZ-CF780.0 ± 87.50.692 ± 0.0980.132 ± 0.009 b4.91 ± 0.72 b0.38 ± 0.07 b
a Statistical differences with JMLP01B0. b Statistical differences with JMLP01BT.
Table 5. Differences (t1t0) of in vivo skin integrity tests with JMLP01B0, JMLP01BT, and BFZ-CF (mean ± SD, n = 10). Statistical differences between formulations are indicated with superscripts.
Table 5. Differences (t1t0) of in vivo skin integrity tests with JMLP01B0, JMLP01BT, and BFZ-CF (mean ± SD, n = 10). Statistical differences between formulations are indicated with superscripts.
FormulationsSkin Elasticity
Mean Differences (AU)
Skin Hydration
Mean Differences (AU)
TEWL
Mean Differences (g/h·m2)
JMLP01B00.0521 ± 0.003 *8.1 ± 7.83.6 ± 1.4
JMLP01BT0.2232 ± 0.033.6 ± 1.43.6 ± 1.9
BFZ-CF0.1605 ± 0.046.0 ± 6.51.5 ± 2.3
AU: Arbitrary units. * Statistical differences with JMLP01BT and BFZ-CF.

Share and Cite

MDPI and ACS Style

Suñer-Carbó, J.; Calpena-Campmany, A.; Halbaut-Bellowa, L.; Clares-Naveros, B.; Rodriguez-Lagunas, M.J.; Barbolini, E.; Zamarbide-Losada, J.; Boix-Montañés, A. Biopharmaceutical Development of a Bifonazole Multiple Emulsion for Enhanced Epidermal Delivery. Pharmaceutics 2019, 11, 66. https://0-doi-org.brum.beds.ac.uk/10.3390/pharmaceutics11020066

AMA Style

Suñer-Carbó J, Calpena-Campmany A, Halbaut-Bellowa L, Clares-Naveros B, Rodriguez-Lagunas MJ, Barbolini E, Zamarbide-Losada J, Boix-Montañés A. Biopharmaceutical Development of a Bifonazole Multiple Emulsion for Enhanced Epidermal Delivery. Pharmaceutics. 2019; 11(2):66. https://0-doi-org.brum.beds.ac.uk/10.3390/pharmaceutics11020066

Chicago/Turabian Style

Suñer-Carbó, Joaquim, Ana Calpena-Campmany, Lyda Halbaut-Bellowa, Beatriz Clares-Naveros, María José Rodriguez-Lagunas, Elena Barbolini, Joanna Zamarbide-Losada, and Antonio Boix-Montañés. 2019. "Biopharmaceutical Development of a Bifonazole Multiple Emulsion for Enhanced Epidermal Delivery" Pharmaceutics 11, no. 2: 66. https://0-doi-org.brum.beds.ac.uk/10.3390/pharmaceutics11020066

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