Lipid-Bilayer Permeation of Drug-Like Compounds

Improving the Accuracy of Permeability Data to Gain Predictive Power: Assessing Sources of Variability in Assays Using Cell Monolayers

The ability to predict the rate of permeation of new compounds across biological membranes is of high importance for their success as drugs, as it determines their efficacy, pharmacokinetics, and safety profile. In vitro permeability assays using Caco-2 monolayers are commonly employed to assess permeability across the intestinal epithelium, with an extensive number of apparent permeability coefficient (Papp) values available in the literature and a significant fraction collected in databases. The compilation of these Papp values for large datasets allows for the application of artificial intelligence tools for establishing quantitative structure-permeability relationships (QSPRs) to predict the permeability of new compounds from their structural properties. One of the main challenges that hinders the development of accurate predictions is the existence of multiple Papp values for the same compound, mostly caused by differences in the experimental protocols employed. This review addresses the magnitude of the variability within and between laboratories to interpret its impact on QSPR modelling, systematically and quantitatively assessing the most common sources of variability. This review emphasizes the importance of compiling consistent Papp data and suggests strategies that may be used to obtain such data, contributing to the establishment of robust QSPRs with enhanced predictive power.

Design Principles for Balancing Lipophilicity and Permeability in beyond Rule of 5 Space

An ab initio conformational analysis of oral beyond Rule of 5 (bRo5) drugs was complemented with measured permeability and logP(octanol) to derive design principles conferring oral bioavailability. 3D polar surface area (PSA) thresholds for oral bRo5 drugs coincided with those reported for Ro5 space. The majority of oral bRo5 drugs exceeded the Ro5 logP threshold of 5, reflecting a bias for permeability. Above 500 Da molecular weight (MW), oral drugs and highly permeable Novartis compounds occupy a narrow polarity range (topological or TPSA/MW) of 0.1-0.3 Å2 /Da, whose upper half coincides with the lower 90 percentiles of the Novartis logP set. This TPSA/MW range and 3D PSA below 100 Å2 define the "Rule of ~1 /₅" for balancing lipophilicity and permeability. Neutral TPSA, defined as TPSA minus 3D PSA occurs independent of conformation, intramolecular hydrogen bonds (IMHB) and MW, suggesting it is an intrinsic molecular property. Neutral TPSA increased in the lead optimization (LO) campaigns of three first in class de novo designed bRo5 drugs and may be a useful design parameter in bRo5 space.

Lipophilicity Modulations by Fluorination Correlate with Membrane Partitioning

Bioactive compounds generally need to cross membranes to arrive at their site of action. The octanol-water partition coefficient (lipophilicity, logPOW ) has proven to be an excellent proxy for membrane permeability. In modern drug discovery, logPOW and bioactivity are optimized simultaneously, for which fluorination is one of the relevant strategies. The question arises as to which extent the often subtle logP modifications resulting from different aliphatic fluorine-motif introductions also lead to concomitant membrane permeability changes, given the difference in molecular environment between octanol and (anisotropic) membranes. It was found that for a given compound class, there is excellent correlation between logPOW values with the corresponding membrane molar partitioning coefficients (logKp ); a study enabled by novel solid-state 19 F NMR MAS methodology using lipid vesicles. Our results show that the factors that cause modulation of octanol-water partition coefficients similarly affect membrane permeability.

Effects of Essential Oils and Fragrant Compounds on Appetite: A Systematic Review

Appetite dysregulation is one of the factors contributing to anorexia, bulimia nervosa, obesity, and diabetes. Essential oils or fragrant compounds have been proven to regulate food intake and energy expenditure; hence, this study aimed to summarize their effects on appetite and the underlying mechanisms. The PubMed and Web of Science databases were searched until July 2022. Only two of the 41 studies were performed clinically, and the remaining 39 used animal models. Oral administration was the most common route, and a dosage range of 100-2000 mg/kg for mice or 2-32 mg/kg for rats was applied, with a duration of 12 days to 4 weeks, followed by inhalation (10^-6-10^-3 mg/cage or 10^-9-10^-2 mg/cm³ within 1 h). Approximately 11 essential oil samples and 22 fragrant compounds were found to increase appetite, while 12 essential oils and seven compounds decreased appetite. These fragrant components can exert appetite-regulating effects via leptin resistance, the activity of sympathetic/parasympathetic nerves, or the mRNA expression of neuropeptide Y (NPY)/agouti-related protein (AgRP), cocaine- and amphetamine-regulated transcript (CART)/proopiomelanocortin (POMC) in the hypothalamus. Fragrance memory and cognitive processes may also play roles in appetite regulation. The findings of this study accentuate the potential of essential oils and fragrant compounds to regulate appetite and eating disorders.

Selective and reversible disruption of mitochondrial inner membrane protein complexes by lipophilic cations

Triphenylphosphonium (TPP) derivatives are commonly used to target chemical into mitochondria. We show that alkyl-TPP cause reversible, dose- and hydrophobicity-dependent alterations of mitochondrial morphology and function and a selective decrease of mitochondrial inner membrane proteins including subunits of the respiratory chain complexes, as well as components of the mitochondrial calcium uniporter complex. The treatment with alkyl-TPP resulted in the cleavage of the pro-fusion and cristae organisation regulator Optic atrophy-1. The structural and functional effects of alkyl-TPP were found to be reversible and not merely due to loss of membrane potential. A similar effect was observed with the mitochondria-targeted antioxidant MitoQ.

A story of peptides, lipophilicity and chromatography - back and forth in time

Peptides, as part of the beyond the rule of 5 (bRo5) chemical space, represent a unique class of pharmaceutical compounds. Because of their exceptional position in the chemical space between traditional small molecules (molecular weight (MW) < 500 Da) and large therapeutic proteins (MW > 5000 Da), peptides became promising candidates for targeting challenging binding sites, including even targets traditionally considered as undruggable - e.g. intracellular protein-protein interactions. However, basic knowledge about physicochemical properties that are important for a drug to be membrane permeable is missing but would enhance the drug discovery process of bRo5 molecules. Consequently, there is a demand for quick and simple lipophilicity determination methods for peptides. In comparison to the traditional lipophilicity determination methods via shake flask and in silico prediction, chromatography-based methods could have multiple benefits such as the requirement of low analyte amount, insensitivity to impurities and high throughput. Herein we elucidate the role of peptide lipophilicity and different lipophilicity values. Further, we summarize peptide analysis via common chromatographic techniques, in specific reversed phase liquid chromatography, hydrophilic interaction liquid chromatography and supercritical fluid chromatography and their role in drug discovery and development process.

Membrane Permeability of Modified Butyltriphenylphosphonium Cations

The alkyltriphenylphosphonium (TPP) group is the most widely used vector targeted to mitochondria. Previously, the length of the alkyl linker was varied as well as structural modifications in the TPP phenyl rings to obtain the optimal therapeutic effect of a pharmacophore conjugated with a lipophilic cation. In the present work, we synthesized butyltriphenylphosphonium cations halogenated and methylated in phenyl rings (C₄TPP-X) and measured electrical current through a planar lipid bilayer in the presence of C₄TPP-X. The permeability of C₄TPP-X varied in the range of 6 orders of magnitude and correlates well with the previously measured translocation rate constant for dodecyltriphenylphosphonium analogues. The partition coefficient of the butyltriphenylphosphonium analogues obtained by calculating the difference in the free energy of cation solvation in water and octane using quantum chemical methods correlates well with the permeability values. Using an ion-selective electrode, a lower degree of accumulation of analogues with halogenated phenyl groups was found on isolated mitochondria of rat liver, which is in agreement with their permeability decrease. Our results indicate the translocation of the butyltriphenylphosphonium cations across the hydrophobic membrane core as rate-limiting stage in the permeability process rather than their binding/release to/from the membrane.

Recent Experimental Developments in Studying Passive Membrane Transport of Drug Molecules

The ability to measure the passive membrane permeation of drug-like molecules is of fundamental biological and pharmaceutical importance. Of significance, passive diffusion across the cellular membranes plays an effective role in the delivery of many pharmaceutical agents to intracellular targets. Hence, approaches for quantitative measurement of membrane permeability have been the topics of research for decades, resulting in sophisticated biomimetic systems coupled with advanced techniques. In this review, recent developments in experimental approaches along with theoretical models for quantitative and real-time analysis of membrane transport of drug-like molecules through mimetic and living cell membranes are discussed. The focus is on time-resolved fluorescence-based, surface plasmon resonance, and second-harmonic light scattering approaches. The current understanding of how properties of the membrane and permeant affect the permeation process is discussed.

Permeation Studies across Symmetric and Asymmetric Membranes in Microdroplet Arrays

We investigated the permeation of molecules across lipid membranes on an open microfluidic platform. An array of droplet pairs was created by spotting aqueous droplets, dispersed in a lipid oil solution, onto a plate with cavities surrounded by a hydrophobic substrate. Droplets in two adjacent cavities come in contact and form an artificial lipid bilayer, called a droplet interface bilayer (DIB). The method allows for monitoring permeation of fluorescently tagged compounds from a donor droplet to an acceptor droplet. A mathematical model was applied to describe the kinetics and determine the permeation coefficient. We also demonstrate that permeation kinetics can be followed over a series of droplets, all connected via DIBs. Moreover, by changing the lipid oil composition after spotting donor droplets, we were able to create asymmetric membranes that we used to mimic the asymmetry of the cellular plasma membrane. Finally, we developed a protocol to separate and extract the droplets for label-free analysis of permeating compounds by liquid chromatography-mass spectrometry. Our versatile platform has the potential to become a new tool for the screening of drug membrane permeability in the future.

Volatile molecules for COVID-19: A possible pharmacological strategy?

COVID-19 is a novel coronavirus disease with a higher incidence of bilateral pneumonia and pleural effusion. The high pulmonary tropism and contagiousness of the virus, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), have stimulated new approaches to combat its widespread diffusion. In developing new pharmacological strategies, the chemical characteristic of volatility can add therapeutic value to the hypothetical drug candidate. Volatile molecules are characterized by a high vapor pressure and are consequently easily exhaled by the lungs after ingestion. This feature could be exploited from a pharmacological point of view, reaching the site of action in an uncommon way but allowing for drug delivery. In this way, a hypothetical molecule for COVID-19 should have a balance between its lung exhalation characteristics and both antiviral and anti-inflammatory pharmacological action. Here, the feasibility, advantages, and disadvantages of a therapy based on oral administration of possible volatile drugs for COVID-19 will be discussed. Both aerosolized antiviral therapy and oral intake of volatile molecules are briefly reviewed, and an evaluation of 1,8-cineole is provided in view of a possible clinical use and also for asymptomatic COVID-19.

Updating the portfolio of physicochemical descriptors related to permeability in the beyond the rule of 5 chemical space

Beyond rule of 5 (bRo5) molecules are attracting significant interest in modern drug discovery mostly because many novel targets require large and more flexible structures. The main aim of this paper is the identification of ad hoc bRo5 physicochemical descriptors of ionization, lipophilicity, polarity and chameleonicity and their measurement. We used different methods to collect ionization (pK_a measures and log k'80 PLRP-S trends), lipophilicity (in octanol/water, in apolar systems and in biomimetic environments), polarity (Δlog P_oct-tol, EPSA and Δlog K_W^IAM) and chameleonicity (ChameLogD) descriptors for 26 bRo5 drugs. A second aim was to check the relationship between physicochemical descriptors and permeability for a subset of compounds for which solid permeability values are reported in the literature. Results showed that the physicochemical profile in the bRo5 chemical space is often experimentally accessible, albeit more tools are required to overcome limitations of individual methods. For the investigated compounds, permeability is governed by Δlog P_oct-tol and preliminary data support that chameleonicity could also have an impact.

5-Fluorouracil-Complete Insight into Its Neutral and Ionised Forms

5-Fluorouracil (5FU), a common anti-cancer drug, occurs in four tautomeric forms and possesses two potential sites of both protonation and deprotonation. Tautomeric and resonance structures of the ionized forms of 5FU create the systems of connected equilibriums. Since there are contradictory reports on the ionized forms of 5FU in the literature, complex theoretical studies on neutral, protonated and deprotonated forms of 5FU, based on the broad spectrum of DFT methods, are presented. These indicate that the O4 oxygen is more willingly protonated than the O2 oxygen and the N1 nitrogen is more willingly deprotonated than the N3 nitrogen in a gas phase. Such preferences are due to advantageous charge delocalization of the respective ions, which is demonstrated by the NBO and ESP analyses. In an aqueous phase, stability differences between respective protonated and deprotonated forms of 5FU are significantly diminished due to the competition between the mesomeric effect and solvation. The calculated pK_a values of the protonated, neutral and singly deprotonated 5FU indicate that 5FU does not exist in the protonated and double-deprotonated forms in the pH range of 0–14. The neutral form dominates below pH 8 and the N1 deprotonated form dominates above pH 8.

Impact of Molecular Partitioning and Partial Equilibration on the Estimation of Diffusion Coefficients from Release Experiments

The present work addresses the effect of partial equilibration and molecular partitioning on the interpretation of release experiments. In this regard, it is shown how release profiles and the values of extracted transport parameters are affected by the time protocol chosen for sample collection by considering a series of experiments where the latter is systematically varied. Caffeine is investigated as a main model drug because of its similar affinity for water and lipids, while monolinolein-based lipid cubic phases are chosen as host matrices because of their wide employment in release studies. Our findings point to a progressive decline in diffusion rate upon increasing the time step, that is, the gap in time between two consecutive pickups, which is a signature of increasing equilibration of caffeine concentration between the lipidic mesophase and the water phase. Furthermore, the amount of released molecules at the first pickup displays negligible changes for large time steps, indicating complete equilibration in such cases. A model is introduced based on Fick's diffusion which goes beyond the assumption of perfect-sink conditions, a common feature of the typical theoretical approaches hitherto developed. The model is shown to account quantitatively for the experimental data and is subsequently employed to clarify the interplay of the adopted release protocol with the various transport parameters in determining the final outcome of the release process. Particularly, two additional molecular drugs are considered, namely glucose and proflavine, which are, respectively, more hydrophilic and hydrophobic than caffeine, thus allowing elucidating the role of molecular partitioning.

A look at the smelly side of physiology: transport of short chain fatty acids

Fermentative organs such as the caecum, the colon, and the rumen have evolved to produce and absorb energy rich short chain fatty acids (SCFA) from otherwise indigestible substrates. Classical models postulate diffusional uptake of the undissociated acid (HSCFA). However, in net terms, a major part of SCFA absorption occurs with uptake of Na⁺ and resembles classical, coupled electroneutral NaCl transport. Considerable evidence suggests that the anion transporting proteins expressed by epithelia of fermentative organs are poorly selective and that their main function may be to transport acetate^-, propionate^-, butyrate^- and HCO₃^- as the physiologically relevant anions. Apical uptake of SCFA thus involves non-saturable diffusion of the undissociated acid (HSCFA), SCFA^-/HCO₃^- exchange via DRA (SLC26A3) and/or SCFA^--H⁺ symport (MCT1, SLC16A1). All mechanisms lead to cytosolic acidification with stimulation of Na⁺/H⁺ exchange via NHE (SLC9A2/3). Basolaterally, Na⁺ leaves via the Na⁺/K⁺-ATPase with recirculation of K⁺. Na⁺ efflux drives the transport of SCFA^- anions through volume-regulated anion channels, such as maxi-anion channels (possibly SLCO2A1), LRRC8, anoctamins, or uncoupled exchangers. When luminal buffering is inadequate, basolateral efflux will increasingly involve SCFA^-/ HCO₃^- exchange (AE1/2, SCL4A1/2), or efflux of SCFA^- with H⁺ (MCT1/4, SLC16A1/3). Furthermore, protons can be basolaterally removed by NHE1 (SCL9A1) or NBCe1 (SLC4A4). The purpose of these transport proteins is to maximize the amount of SCFA transported from the tightly buffered ingesta while minimizing acid transport through the epithelium. As known from the rumen for many decades, a disturbance of these processes is likely to cause severe colonic disease.

The Flip-Flop Diffusion Mechanism across Lipids in a Hybrid Bilayer Membrane

In this study, we examine the mechanism of flip-flop diffusion of proton carriers across the lipid layer of a hybrid bilayer membrane (HBM). The HBM consists of a lipid monolayer appended on top of a self-assembled monolayer containing a Cu-based O2 reduction catalyst on a Au electrode. The flip-flop diffusion rates of the proton carriers dictate the kinetics of O2 reduction by the electrocatalyst. By varying both the tail lengths of the proton carriers and the lipids, we find the combinations of lengths that maximize the flip-flop diffusion rate. These experimental results combined with biophysical modeling studies allow us to propose a detailed mechanism for transmembrane flip-flop diffusion in HBM systems, which involves the bending of the alkyl tail of the proton carrier as the rate-determining step. Additional studies with an unbendable proton carrier further validate these mechanistic findings.

When barriers ignore the "rule-of-five"

Why are a few drugs with properties beyond the rule of 5 (bRo5) absorbed across the intestinal mucosa while most other bRo5 compounds are not? Are such exceptional bRo5 compounds exclusively taken up by carrier-mediated transport or are they able to permeate the lipid bilayer (passive lipoidal diffusion)? Our experimental data with liposomes indicate that tetracycline, which violates one rule of the Ro5, and rifampicin, violating three of the rules, significantly permeate a phospholipid bilayer with kinetics similar to labetalol and metoprolol, respectively. Published data from experimental work and molecular dynamics simulations suggest that the formation of intramolecular H-bonds and the possibility to adopt an elongated shape besides the presence of a significant fraction of net neutral species facilitate lipid bilayer permeation. As an alternative to lipid bilayer permeation, carrier proteins can be targeted to improve absorption, with the potential drawbacks of drug-drug interactions and non-linear pharmacokinetics.

Structural perturbation of a dipalmitoylphosphatidylcholine (DPPC) bilayer by warfarin and its bolaamphiphilic analogue: A molecular dynamics study

Compounds with nominally similar biological activity may exhibit differential toxicity due to differences in their interactions with cell membranes. Many pharmaceutical compounds are amphiphilic and can be taken up by phospholipid bilayers, interacting strongly with the lipid-aqueous interface whether or not subsequent permeation through the bilayer is possible. Bolaamphiphilic compounds, which possess two hydrophilic ends and a hydrophobic linker, can likewise undergo spontaneous uptake by bilayers. While membrane-spanning bolaamphiphiles can stabilize membranes, small molecules with this characteristic have the potential to create membrane defects via disruption of bilayer structure and dynamics. When compared to the amphiphilic therapeutic anticoagulant, warfarin, the bolaamphiphilic analogue, brodifacoum, exhibits heightened toxicity that goes beyond superior inhibition of the pharmacological target enzyme. We explore, herein, the consequences of anticoagulant accumulation in a dipalmitoylphosphatidylcholine (DPPC) bilayer. Coarse-grained molecular dynamics simulations reveal that permeation of phospholipid bilayers by brodifacoum causes a disruption of membrane barrier function that is driven by the bolaamphiphilic nature and size of this molecule. We find that brodifacoum partitioning into bilayers causes membrane thinning and permeabilization and promotes lipid flip-flop - phenomena that are suspected to play a role in triggering cell death. These phenomena are either absent or less pronounced in the case of the less toxic, amphiphilic compound, warfarin.

Polar interactions drug/phospholipids estimated by IAM-HPLC vs cultured cell line passage data: Their relationships and comparison of their effectiveness in predicting drug human intestinal absorption

The relationships between data of passage through Caco-2 cultured cell lines (log Papp), taken from the literature, for 38 structurally unrelated compounds and both n-octanol lipophilicity parameters (log P(N) and log D(7.4)) and phospholipid affinity indexes were investigated. Phospholipid affinity(log k W(IAM)) was experimentally determined by HPLC on two different phospholipid stationary phases and the polar/electrostatic interaction component drug/phospholipids (Δ log k W(IAM)) was calculated according to a method we previously proposed. Log Papp moderately related to lipophilicity values measured at pH 7.4 (log D(7.4)), according to a parabolic pattern, but poorly related with log k W(IAM). Furthermore, a significant inverse linear relationship with Δ l og k W(IAM) values was only observed for the analytes with m.w. >300 Da, for which paracellular diffusion can be considered a minor transport route in vivo. Indeed, it has been reported that Caco-2 passage data also encode secondary passage mechanisms, which participate in a different extent to the jejunal absorption in vivo and cannot be directly equated to the corresponding human in situ log Peff values, unless a normalization is performed. In an attempt to elucidate this issue, 47 structurally unrelated compounds whose cultured cell line passage data were corrected for the effects of the aqueous boundary layer and paracellular permeability, so as to express transcellular intrinsic permeability, log P 0(Caco-2/MDCK), were also considered. Highly significant inverse linear relationships were observed between log P 0(Caco-2/MDCK) and Δlog k W(IAM) values from both IAM.PC.MG (r(2)=0.765) and IAM.PC.DD2 (r(2)=0.806) stationary phases whereas the relationships with either lipophilicity in n-octanol or log k W(IAM) values were very poor. The results of the present study, in complete agreement with those of our recent study on the relationships between jejunal absorption data measured in situ and Δ log k W(IAM) values, confirm the soundness of Δ log k W(IAM) parameters in the prediction of the intestinal absorption of drugs. From a mechanistic point of view, they suggest that the polar/electrostatic forces between drugs and phospholipids play a major role in the passage through biomembranes.

Quantitative aspects of drug permeation across in vitro and in vivo barriers

The kinetics of permeation across epithelial and endothelial cell sheets and across cell membranes is determinant for the pharmacokinetics of a drug. In vitro transport experiments with cultured cells or artificial barriers have tremendously improved the predictability of the in vivo behaviour of tested compounds. This article focuses on the parameters and calculation methods that are used to describe permeation quantitatively, with a focus on in vitro experiments and the prediction of intestinal absorption and blood-brain barrier passage. It shows under which in vitro experimental conditions standard calculations are adequate and under which conditions equations should be adapted to the experimental details. The impact of volume differences between donor and receiver compartments, pH gradients, addition of albumin, accumulation in the barrier and unidirectional transport by an efflux transporter on the results is shown in simulations. The article should make researchers aware of experimental factors that affect the outcome of a permeation experiment and how to account for this during data analysis. Finally, strategies to predict the in vivo behaviour of a compound based on the in vitro data are discussed. The goal of the article is to support researchers in choosing experimental conditions and calculation methods that deliver appropriate and reproducible results in permeation studies in vitro.

Influence of pH-dependent aquatic toxicity of ionizable pharmaceuticals on risk assessments over environmental pH ranges

Due to variation in pH, ionizable pharmaceuticals entering aquatic environments experience different degrees of ionization, which may affect toxicity. Using data from toxicity testing at only neutral pH may potentially under- or overestimate actual toxicity at pH ranges found in natural aquatic environments. Here we show relative pH-dependent acute toxicity to Daphnia magna for the pharmaceutical weak bases fluoxetine and sertraline, as well as the weak acids naproxen, diclofenac, ibuprofen and ketoprofen. A probabilistic modelling approach using the pH-dependent toxicity data for D. magna and an environmental pH distribution based on over 4000 European running waters from 21 countries predicted that environmental toxicity for the weak bases may be underestimated if pH 7 was assumed. The model predicted median underestimation by a factor of 3, with 90% of the model results ranging from 1 to 6. Consequently, due to the slightly basic nature of the European running waters, weak acid toxicity may be overestimated by a factor of 2. Predicted median toxicity was 0.5 of that assuming pH 7 with 90% of the results ranging from 0.03 to 5. Because aquatic pH exhibits large variation both within and between countries, we advise the use of site-specific risk assessments for ionizable pharmaceuticals in making informed water management decisions.

Relating Caco-2 permeability to molecular properties using block relevance analysis

New insights into physicochemical drivers for Caco-2 permeability have been gained using a modelling technique known as block relevance analysis.

Complex Interplay between the P-Glycoprotein Multidrug Efflux Pump and the Membrane: Its Role in Modulating Protein Function

Multidrug resistance in cancer is linked to expression of the P-glycoprotein multidrug transporter (Pgp, ABCB1), which exports many structurally diverse compounds from cells. Substrates first partition into the bilayer and then interact with a large flexible binding pocket within the transporter's transmembrane regions. Pgp has been described as a hydrophobic vacuum cleaner or an outwardly directed drug/lipid flippase. Recent X-ray crystal structures have shed some light on the nature of the drug-binding pocket and suggested routes by which substrates can enter it from the membrane. Detergents have profound effects on Pgp function, and several appear to be substrates. Biochemical and biophysical studies in vitro, some using purified reconstituted protein, have explored the effects of the membrane environment. They have demonstrated that Pgp is involved in a complex relationship with its lipid environment, which modulates the behavior of its substrates, as well as various functions of the protein, including ATP hydrolysis, drug binding, and drug transport. Membrane lipid composition and fluidity, phospholipid headgroup and acyl chain length all influence Pgp function. Recent studies focusing on thermodynamics and kinetics have revealed some important principles governing Pgp-lipid and substrate-lipid interactions, and how these affect drug-binding and transport. In some cells, Pgp is associated with cholesterol-rich microdomains, which may modulate its functions. The relationship between Pgp and cholesterol remains an open question; however, it clearly affects several aspects of its function in addition to substrate-membrane partitioning. The action of Pgp modulators appears to depend on their membrane permeability, and membrane fluidizers and surfactants reverse drug resistance, likely via an indirect mechanism. A detailed understanding of how the membrane affects Pgp substrates and Pgp's catalytic cycle may lead to new strategies to combat clinical drug resistance.

A liposomal fluorescence assay to study permeation kinetics of drug-like weak bases across the lipid bilayer

Lipid bilayer permeation is considered the major route for in vivo barrier passage of drugs. Despite this fact, no technique is currently available to measure the kinetics of permeation across a single lipid bilayer of structurally unrelated drug-like solutes. We developed a liposomal fluorescence assay capable to determine permeation kinetics of basic drug-like solutes across lipid bilayers. The assay is based on the hypothesis that permeation of a weak base along a concentration gradient results in net proton release at the cis-side and net proton capture at the trans-side of the bilayer. The resulting pH changes were monitored with pH-sensitive fluorophores: Test compounds were incubated with liposomes containing a pH-sensitive fluorophore at the bilayer surfaces or in the aqueous lumen and fluorescence changes were monitored with a stopped-flow apparatus in solution or by total internal reflection fluorescence microscopy with surface-captured liposomes on a microfluidic platform. Incubation with lipophilic basic drugs resulted in the expected fluorescence changes while incubation with compounds without basic functionality or high polarity did not affect fluorescence. Kinetics of fluorescence changes followed bi-exponential functions. Logarithmic permeation coefficients (logPermapp) determined in solution and by microfluidics technology showed a good correlation (r(2)=0.94, n=7) and logPermapp increased with increasing lipophilicity. Neither diffusion in the aqueous phase nor partitioning into the bilayer was rate-limiting. PEGylation of 2% of the liposomal lipids reduced Permapp by a factor ~300. In conclusion, the presented liposomal fluorescence assay is capable to determine permeation kinetics of weak basic drug-like solutes across lipid bilayers. The method is adaptable to microfluidics technology for high-throughput measurements and can potentially be modified to work for weak acid solutes.

Modeling bioavailability to organs protected by biological barriers

Computational pharmacokinetic (PK) modeling gives access to drug concentration vs. time profiles in target organs and allows better interpretation of clinical observations of therapeutic or toxic effects. Physiologically-based PK (PBPK) models in particular, based on mechanistic descriptions of the body anatomy and physiology, may also help to extrapolate in vitro or animal data to human.

Once in the systemic circulation, a chemical has access to the microvasculature of every organ or tissue. However, its penetration in the brain, retina, thymus, spinal cord, testis, placenta,… may be limited or even fully prevented by dynamic physiological blood-tissue barriers. Those barriers are both physical (involving tight junctions between adjacent cells) and biochemical (involving metabolizing enzymes and transporters).

On those cases, correct mechanistic characterization of the passage (or not) of molecules through the barrier can be crucial for improved PBPK modeling and prediction.

In parallel, attempts to understand and quantitatively characterize the processes involved in drug penetration of physiological barriers have led to the development of several in vitro experimental models. Data from such assays are very useful to calibrate PBPK models.

We review here those in vitro and computational models, highlighting the challenges and perspectives for in vitro and computational models to better assess drug availability to target tissues.

Physicochemical characterization of a PEGylated liposomal drug formulation using capillary electrophoresis

In this work, the applicability of using CE to perform a physicochemical characterization of a PEGylated liposomal drug formulation of the anti-cancer agent oxaliplatin was examined. Characterization of the liposomal drug formulation using CE instrumentation encompassed: determination of the electrophoretic mobilities, size determination by Taylor dispersion analysis and interaction studies. Electrophoretic mobilities determined by CE were compared with the results obtained by laser Doppler electrophoresis, which were found to be subject to larger variation. Average hydrodynamic diameters of the liposome preparations, as determined by Taylor dispersion analysis, were in the range of 61-84 nm and were compared with the results obtained by dynamic light scattering. Interactions between oxaliplatin (and paracetamol) and the PEGylated liposome were non-detectable by CE frontal analysis as well as by liposome electrokinetic chromatography. In contrast, for the more lipophilic compound propranolol, apparent liposome-aqueous phase distribution coefficients (D(lip) ) were successfully determined by both electrokinetic chromatography (log D(lip) =2.10) and by CE frontal analysis (log D(lip) =2.14). It is envisioned that CE and capillary-based techniques, including Taylor dispersion analysis, will be useful tools for the characterization of nanoparticulate (e.g. liposomal) drug formulations.

The pH-dependent toxicity of basic pharmaceuticals in the green algae Scenedesmus vacuolatus can be explained with a toxicokinetic ion-trapping model

Several previous studies revealed that pharmaceuticals with aliphatic amine function exhibit a higher toxicity toward algae than toward other aquatic organisms. Here we investigated the pH-dependent toxicity of the five basic pharmaceuticals fluoxetine, its metabolite norfluoxetine, propranolol, lidocaine, and trimipramine. For all of them, the toxicity increased with increasing pH when aqueous effect concentrations were considered. Since these pharmaceuticals contain a basic amine group that is protonated and thus positively charged at physiological pH and because algae are capable of biological homeostasis, i.e., pH inside the algal cell remains virtually independent of variable external pH, the speciation of aliphatic amines can be different inside the algal cell compared to the external medium. Therefore, we hypothesized that the high toxicity of aliphatic amines in algae is a toxicokinetic effect caused by speciation and not a toxicodynamic effect caused by a specific mode of toxic action. This hypothesis also implies that internal effect concentrations are independent on external pH. On this basis we developed a simple toxicokinetic model, which assumes that only the neutral molecule is bioavailable and can pass the plasma membrane. This assumption is likely to be valid at pH values down to two units below the acidity constant (pK(a)). For lower pH values a more complex model would have to be evoked that includes, an, albeit smaller, permeability of the charged species. For pH>pK(a)-2, we can safely assume that the outer membrane serves as insulator and that the charged species is formed inside the cell according to the pH in the cytoplasm. Thus this toxicokinetic model is an ion-trapping model. The input parameters of this model are the measured aqueous effect concentrations determined as a function of pH and the membrane-water partitioning, which was modelled by the liposome-water partition coefficients of the neutral and cationic species. They were deduced from experimentally determined liposome-water distribution ratios at various pH values measured with an equilibrium dialysis method. The modelled internal effect concentrations were independent of the external pH and effective membrane burdens were in the same range as for other baseline toxicants found in the literature for algae, daphnids and fish. These results confirm that the higher algal toxicity of pharmaceuticals with an aliphatic amine group can be explained by a toxicokinetic effect and that these pharmaceuticals do not exhibit a specific mode of action in algae but act as baseline toxicants.