Estimation of skin permeation by liquid chromatography

Dermal absorption is a key process in the drug delivery studies of the pharmaceutical and cosmetic industries, as well as in the fields of dermal toxicology, risk assessment, and the exposure of environmental pollutants. This process is typically described by the skin-water permeability coefficient. However, in vivo determination is laborious and expensive. Thus, in the last few years, the development of prediction models from structure descriptors or subrogation through physico-chemical measurements has gained interest. In the present work, a previous subrogation model based on the chromatographic retention on a common C18 column has been tested for a wide set of drugs with very different chemical nature and having a wide range of permeability values. A total of 65 compounds have been used to establish the correlation between skin permeation and the HPLC retention, corrected by the McGowan volume of the drug. Afterwards it was successfully validated in terms of robustness and prediction ability. Finally, the permeability coefficient was estimated for a set of 29 new drugs, and results compared to the ones obtained by other estimation methods, as well as the available in vitro measured values, with very good agreement

Evaluation of the Ability of PAMPA Membranes to Emulate Biological Processes through the Abraham Solvation Parameter Model

Two parallel artificial membrane permeability assay (PAMPA) systems intended for emulating skin permeability have been characterized through the solvation parameter model of Abraham using multilinear regression analysis. The coefficients of the obtained equations have been compared to the ones already established for other PAMPA membranes using statistical tools. The results indicate that both skin membranes are similar to each other in their physicochemical properties. However, they are different from other PAMPA membranes (e.g., intestinal absorption and blood– brain PAMPAs), mainly in terms of hydrophobicity and hydrogen bonding properties. Next, all PAMPA membranes have been compared to relevant biological processes also characterized through the solvation parameter model. The results highlight that skin-PAMPA membranes are a very good choice to emulate skin permeability

Determination of the retention factor of ionizable compounds in microemulsion electrokinetic chromatography

Determination of the retention factor of ionized compounds in microemulsion electrokinetic chromatography requires two mobility measurements at the same pH: one in the presence of the microemulsion and another in plain buffer. However, it has been observed that in some cases subtracting one mobility from another determined in a different medium leads to negative retention factors, which makes no sense from a chemical point of view. This indicates that there is some error in the process which has a direct impact when retention factors are used for further applications. Here, we evaluate how the components of the microemulsion confer different properties to the buffer medium, particularly varying the viscosity parameter (which is inversely related to mobility). Whereas sodium dodecyl sulfate, the surfactant used in the microemulsion, has little effect on the medium viscosity (only an increase of 5%-6%), the presence of 1-butanol, used as a stabilizer, increases it by around 30%. Meanwhile, heptane, which is used as an oil, provokes a slight decrease. Consequently, the mobilities obtained in the microemulsion system are shifted to higher values (less negative mobilities) compared to mobilities obtained in the aqueous buffer, and so one cannot be directly subtracted from the other. Since the microemulsion-buffer medium cannot be directly reproduced, we propose a correction that takes into account the variation of viscosities. This is determined from the electrophoretic mobility of the benzoate ion. As this ion does not interact with the microemulsion, the ratio of its mobilities (measured in plain buffer and microemulsion) is equivalent to the ratio of viscosities, and can be used as the correction factor for other measurements. Thus, mobilities in buffer and microemulsion media are placed on the same scale, overcoming the errors in retention factor determination

Suitability of skin-PAMPA and chromatographic systems to emulate skin permeation. Influence of pH on skin-PAMPA permeability

The skin permeation, Kp, of a chemical compound is a relevant parameter in fields such as toxicology, exposure to pollutants, or dermal studies of pharmaceutical and cosmetic interest. Nonetheless, its experimental determination is not a trivial task, and for this reason alternative methods to estimate Kp have been developed. This work evaluates the suitability of different methodologies to estimate skin permeation of neutral compounds. Three different approaches have been examined: estimation through the skin-PAMPA (Parallel Artificial Membrane Permeability Assay) permeability, Pe, estimation through the chromatographic retention factor combined with molecular volume, and finally estimation through a quantitative structure–property relationship (QSPR) based on the octanol–water partition coefficient, log Po/w, and molecular volume as descriptors. The three approaches have been tested with the same set of compounds and it has been observed that all of them can be used to estimate Kp with similar results, although the chromatographic method presents slightly improved statistics in addition to the facility of measurement. As many drugs are partially ionised at the pH of skin, the influence of pH in skin-PAMPA permeability has been also studied. To this end, the log Pe vs. pH profiles of a set of 25 compounds of different nature have been determined. As expected, the permeation of neutral forms is higher than the one of ionic forms, and permeation of neutral and ionic species are not governed by the same mechanisms

Tadpole toxicity prediction using chromatographic systems

Toxicity has been emulated in tadpole species through chromatographic systems. The parameter studied to evaluate the non-specific toxicity of a compound is the narcosis concentration (Cnar), which is defined as the concentration needed for the immobilization of the organism. Because experimental investigation with animals is lengthy, costly, technically difficult, and ethically questionable, there is a great interest in developing surrogate physicochemical systems able to emulate biological systems to obtain the same information in a faster, more economic, and easier manner. In order to see which chromatographic systems would be able to emulate tadpole narcosis, both, tadpole narcosis data and data in several chromatographic and electrophoretic systems, were fitted to a linear solvation energy relationship (LSER) model. Thus, by comparison of the models it was possible to see which of the chromatographic systems were more similar to the biological one. The physicochemical systems that best emulate tadpole narcosis were an HPLC system based on an immobilized artificial membrane (IAM) column, and two micellar electrokinetic chromatography (MEKC) systems based on sodium taurocholate (STC) and a mixture of sodium dodecylsulphate (SDS) and Brij 35 as surfactants. A system based on a RP18 HPLC column also was selected for comparison because it is a common column in most analytical laboratories. To establish the models, a set of compounds with known Cnar values were analyzed in the chromatographic, and electrophoretic selected systems and, then, the retention factor (k) was correlated to the concentration of narcosis. Statistics showed that the system based on STC micelles was the best to emulate toxicity in tadpoles. The robustness and predictive ability of the developed models were validated

Determination of the aqueous pKa of very insoluble drugs by capillary electrophoresis: Internal standards for methanol-water extrapolation

A fast determination of acidity constants (pKa) of very insoluble drugs has become a necessity in drug discovery process because it often produces molecules that are highly lipophilic and sparingly soluble in water. In this work the high throughput internal standard capillary electrophoresis (IS-CE) method has been adapted to the determination of pKa of water insoluble compounds by measurement in methanol/aqueous buffer mixtures. For this purpose, the reference pKa values for a set of 46 acid-base compounds of varied structure (internal standards) have been established in methanol-water mixtures at several solvent composition levels (with a maximum of 40% methanol). The IS-CE method has been successfully applied to seven test drugs of different chemical nature with intrinsic solubilities lower than 10-6 M. pKa values have been determined at different methanol/aqueous buffer compositions and afterwards Yasuda-Shedlovsky extrapolation method has been applied to obtain the aqueous pKa. The obtained results have successfully been compared to literature ones obtained by other methods. It is concluded that the IS-CE method allows the determination of aqueous pKa values using low proportions of methanol, becoming then more accurate in the extrapolation procedure than other reference methods

Influence of the acid-base ionization of drugs in their retention in reversed-phase liquid chromatography

The effect of the ionization in the RP-HPLC retention of 66 acid-base compounds, most of them drugs of pharmaceutical interest, is studied. The retention time of the compounds can be related to the pH measured in the mobile phase (pwsH) through the sigmoidal equations derived from distribution of the neutral and ionic forms of the drug into the stationary and mobile phases. Fitting of the obtained retention vs. pH profiles provides the retention times of the ionic and neutral forms and the pKa values of the drugs in the mobile phase (pwsKa). The obtained pwsKa values are linearly correlated to the pKa values in water (pwwKa) with two different correlations, one for neutral acids and another for neutral bases that reflect the different influence of the dielectric constant of the medium in ionization of acids and bases. The retention of the neutral species is well correlated to the octanol-water partition coefficient of the drugs as measure of the lipophilicity of the drug, which affects chromatographic retention. Also, the retention time of the ionized forms is related to the retention time of the neutral forms by two different linear correlations, one for anions and the other for cations. These last correlations point out the different retention behaviour of anions and cations: anions are less retained than cations of the same lipophilicity, as measured by the octanol-water partition coefficient of the neutral form. The different retention behaviour of anionic, cationic and neutral forms is confirmed by the hold-up times obtained from different approaches: pycnometry and retention times of anionic (KBr and KI) and neutral (DMSO) markers. Hold-up times obtained by pycnometric measurements agree with those obtained by retention of neutral markers (0.83-0.85 min), whereas hold-up time for anions is mobile phase pH dependent. At acidic pH it is similar to the hold-up time for neutral markers (0.83 min), but then it decreases with the increase of mobile phase pH to 0.65 min at pH 11. The decrease can be explained by the ionization of the silanols of the column and exclusion of anions by charge repulsion. Although not directly measured, the obtained retention data and correlations indicate hold-up time for cations are similar or slightly lower than hold-up time for neutral compounds (0.77-0.83 min). The model proposed and the correlations obtained can be very useful for its implementation in retention prediction algorithms for optimization of separation purpose

Feasibility of the estimation of octanol-water distribution coefficients of acidic drugs by microemulsion electrokinetic chromatography

Previous studies have shown that a microemulsion electrokinetic chromatography (MEEKC) system can estimate the logarithm of the octanol-water partition coefficient (log Po/w) of neutral solutes. In the present work, the applicability of the method to partially and fully ionized acids has been evaluated. Naproxen, a monoprotic acid, has been used as test solute. The retention factor (k) of this compound has been measured in MEEKC at several values of pH and the retention factor-pH profile has been established. As log Po/w correlates with log kMEEKC for neutral compounds, this correlation has been used to estimate the logarithm of the octanol-water partition coefficient of the neutral (log Po/w(HA)), and the fully ionized (log P o/w(A-)) forms of naproxen. Then, the logarithm of the octanol-water distribution coefficient (log Do/w) of the partially ionized form of the acid has been estimated. The comparison of the estimated values with the ones obtained experimentally using the classical procedures, such as the shake-flask method, shows differences under 0.4 log Do/w units either if the acid is partially ionized or in its neutral form in the most part of the pH range. However, the method overestimates the log Do/w of the highly (>99.5 %) or fully ionized form of naproxe

Estimation of the octanol-water distribution coefficient of acidic compounds by microemulsion electrokinetic chromatography

The feasibility of extending the determination of the lipophilicity of partially ionized acids (log Do/w) by microemulsion electrokinetic chromatography (MEEKC) is tested. Theoretical considerations predict that a linear log Do/w vs. log k correlation can be obtained only when the neutral and ionic forms of an acid follow the same correlation equation and the slope of the correlation is unity. In practice, since the lipophilicity of the neutral acid is much higher than that of the ionic form and the correlation slope is not very different from 1, the general linear correlation for neutral compounds can be applied across most of the ionization range of the acid. The linear correlation between log Po/w and log k of 20 neutral solutes (calibration curve) has been established and extended to 6 acids used as models, tested across their full ionization range. log Do/w-pH, and log k-pH profiles have been obtained for these 6 acids, and plotted log Do/w against log k for any acid at any degree of ionization. Furthermore, the log Do/w of the acids has been estimated from the calibration curve and log k-pH profile, and compared to values in the literature determined using reference methods such as the shake-flask one. Accurate values have been obtained using the MEEKC method when the acids are in their neutral form or partially ionized (ionization degree, α < 0.995). However, this parameter is overestimated when the acids are highly or fully ionized (α ≈ 1). Finally, in order to test the applicability of this method, we have applied the same procedure to estimate log Do/w at pH = 7.4 (blood physiological pH) of a set of 30 additional compounds (including partially and fully ionized acids). The results at this pH foll

Optimization of experimental conditions for skin-PAMPA measurements

In recent years, the parallel artificial membrane permeability assay (PAMPA) has been extended for prediction of skin permeation by developing an artificial membrane which mimics the stratum corneum structure, skin-PAMPA. In the present work, the different parameters affecting skin-PAMPA permeability, such as incubation time and stirring, have been studied to establish ideal assay conditions to generate quality data for a screening of active pharmaceutical ingredients (API) in early stage drug discovery. Another important parameter is membrane retention, which shows dependence on lipophilicity when compounds are in their neutral form. Furthermore, the stability of the membrane has been investigated at different pH values, especially at basic pHs. Finally, a good correlation between human skin permeability and skin-PAMPA permeability, with a large dataset (n = 46), has been established. The optimized assay conditions were an incubation time of 4 hours with stirring in a pH below 8. With all these considerations the thickness of the aqueous boundary layer is decreased as much as possible and the membrane stability is guaranteed