Kinetics of lipid bilayer permeation of a series of ionisable drugs and their correlation with human transporter-independent intestinal permeability

Interaction of the Antibiotic Rifampicin with Lipid Membranes

Rifampicin is a broad-spectrum antibiotic, active against several bacterial infections such as tuberculosis. It is a relatively large and structurally complex molecule, including numerous polar groups. Although violating several of Lipinski's rules for efficient intestinal absorption, rifampicin shows good oral bioavailability, permeating through cell membranes in the absorption pathway and those of the target organisms. Some hypotheses have been proposed for its efficient membrane permeation, but the details are mostly unknown. In this work, the interaction of rifampicin with POPC lipid bilayers is studied using experimental biophysics methodologies and atomistic molecular dynamics simulations considering the two most prevalent ionic species at physiological pH, the anionic and the zwitterionic forms. The results show that both ionization forms of rifampicin establish favorable interactions with the membrane lipids, in agreement with the relatively high partition coefficient obtained experimentally. The results from MD simulations and isothermal titration calorimetry using different pH buffers show that the piperazine group inserts deeply in the membrane and is accompanied by a stabilization of its neutral form. The bulky nature of rifampicin and its deep insertion in the membrane lead to a strong perturbation in the lipids local order, decreasing the membrane barrier properties as evaluated from the rate of carboxyfluorescein leaching. Altogether, the comparison between the experimental and MD simulations results provides important insight regarding the rifampicin molecular features responsible for its relatively fast membrane permeation. The lipid POPC used in this study was selected as a simple membrane with relevance for different organisms across all kingdoms. Further studies using more complex lipid compositions will provide details on eventual specificities for rifampicin interaction with the membranes of distinct organisms.

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.

Physical Chemistry of Chloroquine Permeation through the Cell Membrane with Atomistic Detail

We provide a molecular-level description of the thermodynamics and mechanistic aspects of drug permeation through the cell membrane. As a case study, we considered the antimalaria FDA approved drug chloroquine. Molecular dynamics simulations of the molecule (in its neutral and protonated form) were performed in the presence of different lipid bilayers, with the aim of uncovering key aspects of the permeation process, a fundamental step for the drug's action. Free energy values obtained by well-tempered metadynamics simulations suggest that the neutral form is the only permeating protomer, consistent with experimental data. H-bond interactions of the drug with water molecules and membrane headgroups play a crucial role for permeation. The presence of the transmembrane potential, investigated here for the first time in a drug permeation study, does not qualitatively affect these conclusions.

Novel Fluoroquinolones with Possible Antibacterial Activity in Gram-Negative Resistant Pathogens: In Silico Drug Discovery

Antibiotic resistance is a global threat to public health, and the search for new antibacterial therapies is a current research priority. The aim of this in silico study was to test nine new fluoroquinolones previously designed with potential leishmanicidal activity against Campylobacter jejuni, Escherichia coli, Neisseria gonorrhoeae, Pseudomonas aeruginosa, and Salmonella typhi, all of which are considered by the World Health Organization to resistant pathogens of global concern, through molecular docking and molecular dynamics (MD) simulations using wild-type (WT) and mutant-type (MT) DNA gyrases as biological targets. Our results showed that compound 9FQ had the best binding energy with the active site of E. coli in both molecular docking and molecular dynamics simulations. Compound 9FQ interacted with residues of quinolone resistance-determining region (QRDR) in GyrA and GyrB chains, which are important to enzyme activity and through which it could block DNA replication. In addition to compound 9FQ, compound 1FQ also showed a good affinity for DNA gyrase. Thus, these newly designed molecules could have antibacterial activity against Gram-negative microorganisms. These findings represent a promising starting point for further investigation through in vitro assays, which can validate the hypothesis and potentially facilitate the development of novel antibiotic drugs.

Physiologically Based Pharmacokinetic Modeling for Multiple Oral Administration Labetalol in Pregnant Women

BACKGROUND

Labetalol has an irreplaceable role in treating Hypertensive disorders of pregnancy (HDP), a common disease during pregnancy with a prevalence of 5.2-8.2%. However, there were big differences in dosage regimens between various guidelines.

PURPOSE

A physiologically-based pharmacokinetics (PBPK) model was established and validated to evaluate the existing oral dosage regimens, and to compare the difference in plasma concentration between pregnant and non-pregnant women.

METHODS

First, non-pregnant woman models with specific plasma clearance or enzymatic metabolism (UGT1A1, UGT2B7, CYP2C19) were established and validated. For CYP2C19, slow, intermediate, and rapid metabolic phenotypes were considered. Then, a pregnant model with proper structure and parameters adjustment was established and validated against the multiple oral administration data.

RESULTS

The predicted labetalol exposure captured the experimental data well. The following simulations with criteria lowering 15 mmHg blood pressure (corresponding to around 108 ng/ml plasma labetalol) found that the maximum daily dosage in the Chinese guideline may be insufficient for some severe HDP patients. Moreover, similar predicted steady-state trough plasma concentration was found between the maximum daily dosage in the American College of Obstetricians and Gynecologists (ACOG) guideline, 800 mg Q8h and a regimen of 200 mg Q6h. Simulations comparing non-pregnant and pregnant women found that the difference in labetalol exposure highly depended on the CYP2C19 metabolic phenotype.

CONCLUSIONS

In summary, this work initially established a PBPK model for multiple oral administration of labetalol for pregnant women. This PBPK model may lead to personalized labetalol medication in the future.

Ligand's Partition to the Lipid Bilayer Should Be Accounted for When Estimating Their Affinity to Proteins

Ligand-protein interactions are usually studied in complex media that also contain lipids. This is particularly relevant for membrane proteins that are always associated with lipid bilayers, but also for water-soluble proteins studied in in vivo conditions. This work addresses the following two questions: (i) How does the neglect of the lipid bilayer influence the apparent ligand-protein affinity? (ii) How can the intrinsic ligand-protein affinity be obtained? Here we present a framework to quantitatively characterize ligand-protein interactions in complex media for proteins with a single binding site. The apparent affinity obtained when following some often-used approximations is also explored, to establish these approximations' validity limits and to allow the estimation of the true affinities from data reported in literature. It is found that an increase in the ligand lipophilicity or in the volume of the lipid bilayer always leads to a decrease in the apparent ligand-protein affinity, both for water-soluble and for membrane proteins. The only exceptions are very polar ligands (excluded from the lipid bilayer) and ligands whose binding affinity to the protein increases supralinearly with ligand lipophilicity. Finally, this work discusses which are the most relevant parameters to consider when exploring the specificity of membrane proteins.

Approaches to Address PK-PD Challenges of Conventional Liposome Formulation with Special Reference to Cancer, Alzheimer's, Diabetes, and Glaucoma: An Update on Modified Liposomal Drug Delivery System

Liposomes nowadays have become a preferential drug delivery system since they provide facilitating properties to drugs, such as improved therapeutic index of encapsulated drug, target and controlled drug delivery, and less toxicity. However, conventional liposomes have shown some disadvantages, such as less drug loading capacity, poor retention, clearance by kidney or reticuloendothelial system, and less release of hydrophilic drugs. Thus, to overcome these disadvantages recently, scientists have explored new approaches and methods, viz., ligand conjugation, polymer coating, and liposome hybrids, including surface-modified liposomes, biopolymer-incorporated liposomes, guest-in-cyclodextrin-in-liposome, liposome-in-hydrogel, liposome-in-film, liposome-in-nanofiber, etc. These approaches have been shown to improve the physiochemical and pharmacokinetic properties of encapsulated drugs. Lately, pharmacokinetic-pharmacodynamic (PK-PD) computational modeling has emerged as a beneficial tool for analyzing the impact of formulation and system-specific factors on the target disposition and therapeutic efficacy of liposomal drugs. There has been an increasing number of liposome-based therapeutic drugs, both FDA approved and those undergoing clinical trials, having application in cancer, Alzheimer's, diabetes, and glaucoma. To meet the continuous demand of health sectors and to produce the desired product, it is important to perform pharmacokinetic studies. This review focuses on the physical, physicochemical, and chemical factors of drugs that influence the target delivery of drugs. It also explains various physiological barriers, such as systemic clearance and extravasation. A novel approach, liposomal-hybrid complex, an innovative approach as a vesicular drug delivery system to overcome limited membrane permeability and bioavailability, has been discussed in the review. Moreover, this review highlights the pharmacokinetic considerations and challenges of poorly absorbed drugs along with the applications of a liposomal delivery system in improving PKPD in various diseases, such as cancer, Alzheimer's, diabetes, and glaucoma.

Calculation of Permeability Coefficients from Solute Equilibration Dynamics: An Assessment of Various Methods

Predicting the rate at which substances permeate membrane barriers in vivo is crucial for drug development. Permeability coefficients obtained from in vitro studies are valuable for this goal. These are normally determined by following the dynamics of solute equilibration between two membrane-separated compartments. However, the correct calculation of permeability coefficients from such data is not always straightforward. To address these problems, here we develop a kinetic model for solute permeation through lipid membrane barriers that includes the two membrane leaflets as compartments in a four-compartment model. Accounting for solute association with the membrane allows assessing various methods in a wide variety of conditions. The results showed that the often-used expression Papp= β × r/3 is inapplicable to very large or very small vesicles, to moderately or highly lipophilic solutes, or when the development of a significant pH gradient opposes the solute’s flux. We establish useful relationships that overcome these limitations and allow predicting permeability in compartmentalised in vitro or in vivo systems with specific properties. Finally, from the parameters for the interaction of the solute with the membrane barrier, we defined an intrinsic permeability coefficient that facilitates quantitative comparisons between solutes.

Mechanistic Understanding from Molecular Dynamics in Pharmaceutical Research 2: Lipid Membrane in Drug Design

We review the use of molecular dynamics (MD) simulation as a drug design tool in the context of the role that the lipid membrane can play in drug action, i.e., the interaction between candidate drug molecules and lipid membranes. In the standard "lock and key" paradigm, only the interaction between the drug and a specific active site of a specific protein is considered; the environment in which the drug acts is, from a biophysical perspective, far more complex than this. The possible mechanisms though which a drug can be designed to tinker with physiological processes are significantly broader than merely fitting to a single active site of a single protein. In this paper, we focus on the role of the lipid membrane, arguably the most important element outside the proteins themselves, as a case study. We discuss work that has been carried out, using MD simulation, concerning the transfection of drugs through membranes that act as biological barriers in the path of the drugs, the behavior of drug molecules within membranes, how their collective behavior can affect the structure and properties of the membrane and, finally, the role lipid membranes, to which the vast majority of drug target proteins are associated, can play in mediating the interaction between drug and target protein. This review paper is the second in a two-part series covering MD simulation as a tool in pharmaceutical research; both are designed as pedagogical review papers aimed at both pharmaceutical scientists interested in exploring how the tool of MD simulation can be applied to their research and computational scientists interested in exploring the possibility of a pharmaceutical context for their research.

Measuring the oral bioavailability of protein hydrolysates derived from food sources: A critical review of current bioassays

BACKGROUND

Food proteins are a source of hydrolysates with potentially useful biological attributes. Bioactive peptides from food-derived proteins are released from hydrolysates using exogenous industrial processes or endogenous intestinal enzymes. Current in vitro permeability assays have limitations in predicting the oral bioavailability (BA) of bioactive peptides in humans. There are also difficulties in relating the low blood levels of food-derived bioactive peptides detected in preclinical in vivo models to pharmacodynamic read-outs relevant for humans.

SCOPE AND APPROACH

In this review, we describe in vitro assays of digestion, permeation, and metabolism as indirect predictors of the potential oral BA of hydrolysates and their constituent bioactive peptides. We discuss the relationship between industrial hydrolysis processes and the oral BA of hydrolysates and their peptide by-products.

KEY FINDINGS

Hydrolysates are challenging for analytical detection methods due to capacity for enzymatic generation of peptides with novel sequences and also new modifications of these peptides during digestion. Mass spectrometry and peptidomics can improve the capacity to detect individual peptides released from complex hydrolysates in biological milieu.

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.

Synthesis and Structure-Affinity Relationship of Small Molecules for Imaging Human CD80 by Positron Emission Tomography

The costimulatory molecule CD80 is an early marker for immune activation. It is upregulated on activated antigen-presenting cells. We aimed at developing a tracer for imaging CD80 by positron emission tomography (PET). Novel CD80 ligands were synthesized and tested by SPR for affinity to human CD80 (hCD80) and displacement of endogenous ligands. Several compounds bound with one-digit nanomolar affinity to hCD80 and displaced CTLA-4 and CD28 at nanomolar concentrations. A structure-affinity relationship study revealed relevant moieties for strong affinity to hCD80 and positions for further modifications. Lead compound MT107 (7f) was radiolabeled with carbon-11. In vitro, [¹¹C]MT107 showed specific binding to hCD80-positive tissue and high plasma protein binding. In vivo, [¹¹C]MT107 accumulated in liver, gall bladder, and intestines but only scarcely in hCD80-positive xenografts. The unfavorable in vivo performance may result from high plasma protein binding and extensive biliary excretion.

Dynamic Protonation Dramatically Affects the Membrane Permeability of Drug-like Molecules

Permeability (P_m) across biological membranes is of fundamental importance and a key factor in drug absorption, distribution, and development. Although the majority of drugs will be charged at some point during oral delivery, our understanding of membrane permeation by charged species is limited. The canonical model assumes that only neutral molecules partition into and passively permeate across membranes, but there is mounting evidence that these processes are also facile for certain charged species. However, it is unknown whether such ionizable permeants dynamically neutralize at the membrane surface or permeate in their charged form. To probe protonation-coupled permeation in atomic detail, we herein apply continuous constant-pH molecular dynamics along with free energy sampling to study the permeation of a weak base propranolol (PPL), and evaluate the impact of including dynamic protonation on P_m. The simulations reveal that PPL dynamically neutralizes at the lipid-tail interface, which dramatically influences the permeation free energy landscape and explains why the conventional model overestimates the assigned intrinsic permeability. We demonstrate how fixed-charge-state simulations can account for this effect, and propose a revised model that better describes pH-coupled partitioning and permeation. Our results demonstrate how dynamic changes in protonation state may play a critical role in the permeation of ionizable molecules, including pharmaceuticals and drug-like molecules, thus requiring a revision of the standard picture.

Combined Ion-Trapping and Mass Balance Models To Describe the pH-Dependent Uptake and Toxicity of Acidic and Basic Pharmaceuticals in Zebrafish Embryos ( Danio rerio)

The aim of the current study was to understand and develop models to predict the pH-dependent toxicity of ionizable pharmaceuticals in embryos of the zebrafish Danio rerio. We found a higher uptake and toxicity with increasing neutral fraction of acids (diclofenac, genistein, naproxen, torasemide, and warfarin) and bases (metoprolol and propranolol). Simple mass balance models accounting for the partitioning to lipids and proteins in the zebrafish embryo were found to be suitable to predict the bioconcentration after 96 h of exposure if pH values did not differ much from the internal pH of 7.55. For other pH values, a kinetic ion-trap model for the zebrafish embryo explained the pH dependence of biouptake and toxicity. The total internal lethal concentrations killing 50% of the zebrafish embryos (ILC₅₀) were calculated from the measured BCF and LC₅₀. The resulting ILC₅₀ were independent of external pH. Critical membrane concentrations were deduced by an internal mass balance model, and apart from diclofenac, whose specific toxicity in fish had already been established, all pharmaceuticals were confirmed to act as baseline toxicants in zebrafish.

Permeation characteristics of tetracyclines in parallel artificial membrane permeation assay

The purpose of the present study was to characterize the passive permeation of tetracyclines in the parallel artificial membrane permeation assay (PAMPA). Tetracyclines exist as zwitterion at physiological pH. The PAMPA membrane was prepared by impregnating a phospholipid/decane solution to a filter support. The permeation coefficient (P_e) of tetracycline (TC) was markedly affected by the lipid composition of the PAMPA membrane. No permeation was observed when phospholipid was not added (pure decane membrane, P_e < 0.05 × 10^-6 cm/sec). With the addition of 2 % PC, little or no increase in P_e was observed. The addition of 1 % PE increased the P_e value more than tenfold. The addition of 2 % soybean lecithin containing phosphatidylinositol (PI) and phosphatidic acid (PA) increased the P_e value to above 4 × 10^-6 cm/sec. The P_e value was further increased to 15 × 10^-6 cm/sec by increasing the concentration of soybean lecithin from 2 to 10 %. The P_e value showed pH and temperature dependence, whereas it was not affected by the ionic strength, TC concentration, and ion-pair transport inhibitors. A weak correlation was observed between the P_e values and octanol-buffer distribution coefficients of tetracyclines. These results suggest that inter-molecular interactions between TC and PE, PI and/or PA facilitate the passive diffusion of TC across the PAMPA membrane.

Partition of Amphiphilic Molecules to Lipid Bilayers by ITC: Low-Affinity Solutes

A protocol is developed to allow the accurate characterization of partition to lipid bilayers for solutes with low affinity, using isothermal titration calorimetry. The methodology proposed is suitable for studies using complex membranes, such as intact biomembranes or whole cells. In the method developed, the association is characterized at increasing solute concentrations. This allows the characterization of solute partition into unperturbed membranes, as well as effects induced by high solute concentrations. Most druglike molecules are expected to interact with low-to-moderate affinity with relevant cell membranes. This is due to both the need for a relatively high aqueous solubility of the drug and the poor binding properties of the cell membranes. The methodology is applied to characterize the interaction of antibiotic Rifampicin with 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine and with lipid bilayers representative of bacterial membranes.