Effects of lipid composition on membrane permeation

Emerging Designs and Applications for Biomembrane Biosensors

Nature has inspired the development of biomimetic membrane sensors in which the functionalities of biological molecules, such as proteins and lipids, are harnessed for sensing applications. This review provides an overview of the recent developments for biomembrane sensors compatible with either bulk or planar sensing applications, namely using lipid vesicles or supported lipid bilayers, respectively. We first describe the individual components required for these sensing platforms and the design principles that are considered when constructing them, and we segue into recent applications being implemented across multiple fields. Our goal for this review is to illustrate the versatility of nature's biomembrane toolbox and simultaneously highlight how biosensor platforms can be enhanced by harnessing it.

Dehydration of biological membranes in a non-condensing environment

The study of the dehydration process in a cell membrane allows a better understanding of how water is bound to it. While in prior studies, cell dehydration was commonly analyzed under osmotic stress conditions, in the present work, we focus on the dehydration driven by evaporation in a restricted condensing environment. Using a thermogravimetry method, we studied the dehydration of Escherichia coli through isothermal evaporation in the presence of a gas flux. To figure out the loss of mass in this situation, we first evaluated the dynamics of water evaporation of a suspension of multilamellar liposomes. We found that the evaporation of liposomal suspensions composed of individual lipids is constant, although slightly restricted by the presence of liposomes, while the evaporation of liposomal suspensions composed of a mixture of different lipids follows an exponential decay. This is explained considering that the internal pressure at the air-water interface is proportional to the amount of bound water. The evaporation of water from a biomass sample follows this latter behaviour.

Capric Acid and Myristic Acid Permeability Enhancers in Curved Liposome Membranes

Liposomes are considered as advanced drug delivery systems for cancer treatment. A generation of pH-sensitive liposomes is being developed that use fatty acids (FAs) as a trigger for drug release in tumor tissues. However, FAs are also known to enhance permeability, and it is unclear whether FAs in liposomes may cause drug leakage or premature drug release. The passive permeability of the drug through the membrane of the liposome is thus a crucial factor for timely drug delivery. To investigate how the curvature and lipid composition of liposomes affect their passive permeability, coarse-grained molecular dynamics were performed. The permeability was determined with a counting method. Flat bilayers and three liposomes with varying diameters were studied, which had varying lipid compositions of dipalmitoylphosphatidylcholine, cholesterol, and deprotonated or neutral saturated FAs. The investigated permeants were water and two other small permeants, which have different free energy profiles (solubility) across the membrane. First, for the curvature effect, our results showed that curvature increases the water permeability by reducing the membrane thickness. The permeability increase for water is about a factor of 1.7 for the most curved membranes. However, a high curvature decreases permeability for permeants with free energy profiles that are a mix of wells and barriers in the headgroup region of the membrane. Importantly, the type of experimental setup is expected to play a dominant role in the permeability value, i.e., whether permeants are escaping or entering the liposomes. Second, for the composition effect, FAs decrease both the area per lipid (APL) and the membrane thickness, resulting in permeability increases of up to 55%. Cholesterol has a similar effect on the APL but has the opposite impact on membrane thickness and permeability. Therefore, FAs and cholesterol have opposing effects on permeability, with cholesterol's effect being slightly stronger in our simulated bilayers. As all permeability values were well within a factor of 2, and with liposomes usually being larger and less curved in experimental applications, it can be concluded that the passive drug release from a pH-sensitive liposome does not seem to be significantly affected by the presence of FAs.

Advances in Computational Approaches for Estimating Passive Permeability in Drug Discovery

Passive permeation of cellular membranes is a key feature of many therapeutics. The relevance of passive permeability spans all biological systems as they all employ biomembranes for compartmentalization. A variety of computational techniques are currently utilized and under active development to facilitate the characterization of passive permeability. These methods include lipophilicity relations, molecular dynamics simulations, and machine learning, which vary in accuracy, complexity, and computational cost. This review briefly introduces the underlying theories, such as the prominent inhomogeneous solubility diffusion model, and covers a number of recent applications. Various machine-learning applications, which have demonstrated good potential for high-volume, data-driven permeability predictions, are also discussed. Due to the confluence of novel computational methods and next-generation exascale computers, we anticipate an exciting future for computationally driven permeability predictions.

Combining SILCS and Artificial Intelligence for High-Throughput Prediction of the Passive Permeability of Drug Molecules

Membrane permeability of drug molecules plays a significant role in the development of new therapeutic agents. Accordingly, methods to predict the passive permeability of drug candidates during a medicinal chemistry campaign offer the potential to accelerate the drug design process. In this work, we combine the physics-based site identification by ligand competitive saturation (SILCS) method and data-driven artificial intelligence (AI) to create a high-throughput predictive model for the passive permeability of druglike molecules. In this study, we present a comparative analysis of four regression models to predict membrane permeabilities of small druglike molecules; of the tested models, Random Forest was the most predictive yielding an R2 of 0.81 for the independent data set. The input feature vector used to train the developed prediction model includes absolute free energy profiles of ligands through a POPC-cholesterol bilayer based on ligand grid free energy (LGFE) profiles obtained from the SILCS approach. The use of the membrane free energy profiles from SILCS offers information on the physical forces contributing to ligand permeability, while the use of AI yields a more predictive model trained on experimental PAMPA permeability data for a collection of 229 molecules. This combination allows for rapid estimations of ligand permeability at a level of accuracy beyond currently available predictive models while offering insights into the contributions of the functional groups in the ligands to the permeability barrier, thereby offering quantitative information to facilitate rational ligand design.

Analysis of Time-Dependent Pharmacokinetics Using In Vitro-In Vivo Extrapolation and Physiologically Based Pharmacokinetic Modeling

SCR430, a sorafenib derivative, is an investigational drug exhibiting anti-tumor action. This study aimed to have a mechanistic understanding of SCR430's time-dependent pharmacokinetics (TDPK) through an ex vivo study combined with an in vitro-in vivo extrapolation (IVIVE) and physiologically based pharmacokinetic (PBPK) modeling. A non-compartmental pharmacokinetic analysis was performed after intravenous SCR430 administration in female Sprague-Dawley rats for a control group (no treatment), a vehicle group (vehicle only, 14 days, PO), and a repeated-dosing group (SCR430, 30 mg/kg/day, 14 days, PO). In addition, hepatic uptake and metabolism modulation were investigated using isolated hepatocytes from each group of rats. The minimal PBPK model based on IVIVE was constructed to explain SCR430's TDPK. Repeated SCR430 administration decreased the systemic exposure by 4.4-fold, which was explained by increased hepatic clearance (4.7-fold). The ex vivo study using isolated hepatocytes from each group suggested that the increased hepatic uptake (9.4-fold), not the metabolic activity, contributes to the increased hepatic clearance. The minimal PBPK modeling based on an ex vivo study could explain the decreased plasma levels after the repeated doses. The current study demonstrates the TDPK after repeated dosing by hepatic uptake induction, not hepatic metabolism, as well as the effectiveness of an ex vivo approach combined with IVIVE and PBPK modeling to investigate the TDPK.

UV-DIB: label-free permeability determination using droplet interface bilayers

Simple diffusion of molecular entities through a phospholipid bilayer, is a phenomenon of great importance to the pharmaceutical and agricultural industries. Current model lipid systems to probe this typically only employ fluorescence as a readout, thus limiting the range of assessable chemical matter that can be studied. We report a new technology platform, the UV-DIB, which facilitates label free measurement of small molecule translocation rates. This is based upon the coupling of droplet interface bilayer technology with implemented fiber optics to facilitate analysis via ultraviolet spectroscopy, in custom designed PMMA wells. To improve on current DIB technology, the platform was designed to be reusable, with a high sampling rate and a limit of UV detection in the low μM regime. We demonstrate the use of our system to quantify passive diffusion in a reproducible and rapid manner where the system was validated by investigating multiple permeants of varying physicochemical properties across a range of lipid interfaces, each demonstrating differing kinetics. Our system permits the interrogation of structural dependence on the permeation rate of a given compound. We present this ability from two structural perspectives, that of the membrane, and the permeant. We observed a reduction in permeability between pure DOPC and DPhPC interfaces, concurring with literature and demonstrating our ability to study the effects of lipid composition on permeability. In relation to the effects of permeant structure, our device facilitated the rank ordering of various compounds from the xanthine class of compounds, where the structure of each permeant differed by a single group alteration. We found that DIBs were stable up to 5% DMSO, a molecule often used to aid solubilisation of pharmaceutical and agrochemical compounds. The ability of our device to rank-order compounds with such minor structural differences provides a level of precision that is rarely seen in current, industrially applied technologies.

Assessing the Intestinal Permeability of Small Molecule Drugs via Diffusion Motion on a Multidimensional Free Energy Surface

A protocol that accurately assesses the intestinal permeability of small molecule compounds plays an essential role in decreasing the cost and time in inventing a new drug. This manuscript presents a novel computational method to study the passive permeation of small molecule drugs based on the inhomogeneous solubility-diffusion model. The multidimensional free energy surface of the drug transiting through a lipid bilayer is computed with transition-tempered metadynamics that accurately captures the mechanisms of passive permeation. The permeability is computed by following the diffusion motion of the drug molecules along the minimal free energy path found on the multidimensional free energy surface. This computational method is assessed by studying the permeability of five small molecule drugs (ketoprofen, naproxen, metoprolol, propranolol, and salicylic acid). The results demonstrate a remarkable agreement between the computed permeabilities and those measured with the intestinal assay. The in silico method reported in this manuscript also reproduces the permeability measured from the intestinal assay (in vivo) better than the cell-based assays (e.g., PAMPA and Caco-2) do. In addition, the multidimensional free energy surface reveals the interplay between the structure of the small molecule and its permeability, shedding light on strategies of drug optimization.

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.

Molecular dynamics simulations and experimental studies reveal differential permeability of withaferin-A and withanone across the model cell membrane

Poor bioavailability due to the inability to cross the cell membrane is one of the major reasons for the failure of a drug in clinical trials. We have used molecular dynamics simulations to predict the membrane permeability of natural drugs-withanolides (withaferin-A and withanone) that have similar structures but remarkably differ in their cytotoxicity. We found that whereas withaferin-A, could proficiently transverse through the model membrane, withanone showed weak permeability. The free energy profiles for the interaction of withanolides with the model bilayer membrane revealed that whereas the polar head group of the membrane caused high resistance for the passage of withanone, the interior of the membrane behaves similarly for both withanolides. The solvation analysis further revealed that the high solvation of terminal O5 oxygen of withaferin-A was the major driving force for its high permeability; it interacted with the phosphate group of the membrane that led to its smooth passage across the bilayer. The computational predictions were tested by raising and recruiting unique antibodies that react to withaferin-A and withanone. The time-lapsed analyses of control and treated cells demonstrated higher permeation of withaferin-A as compared to withanone. The concurrence between the computation and experimental results thus re-emphasised the use of computational methods for predicting permeability and hence bioavailability of natural drug compounds in the drug development process.

Effects of lipid heterogeneity on model human brain lipid membranes

Cell membranes naturally contain a heterogeneous lipid distribution. However, homogeneous bilayers are commonly preferred and utilised in computer simulations due to their relative simplicity, and the availability of lipid force field parameters. Recently, experimental lipidomics data for the human brain cell membranes under healthy and Alzheimer's disease (AD) conditions were investigated, since disruption to the lipid composition has been implicated in neurodegenerative disorders, including AD [R. B. Chan et al., J. Biol. Chem., 2012, 287, 2678-2688]. In order to observe the effects of lipid complexity on the various bilayer properties, molecular dynamics simulations were used to study four membranes with increasing heterogeneity: a pure POPC membrane, a POPC and cholesterol membrane in a 1 : 1 ratio (POPC-CHOL), and to our knowledge, the first realistic models of a healthy brain membrane and an Alzheimer's diseased brain membrane. Numerous structural, interfacial, and dynamical properties, including the area per lipid, interdigitation, dipole potential, and lateral diffusion of the two simple models, POPC and POPC-CHOL, were analysed and compared to those of the complex brain models consisting of 27 lipid components. As the membranes gain heterogeneity, a number of alterations were found in the structural and dynamical properties, and more significant differences were observed in the lateral diffusion. Additionally, we observed snorkeling behaviour of the lipid tails that may play a role in the permeation of small molecules across biological membranes. In this work, atomistic description of realistic brain membrane models is provided, which can add insight towards the permeability and transport pathways of small molecules across these membrane barriers.

Systems biology approaches to study lipidomes in health and disease

Lipids have many important biological roles, such as energy storage sources, structural components of plasma membranes and as intermediates in metabolic and signaling pathways. Lipid metabolism is under tight homeostatic control, exhibiting spatial and dynamic complexity at multiple levels. Consequently, lipid-related disturbances play important roles in the pathogenesis of most of the common diseases. Lipidomics, defined as the study of lipidomes in biological systems, has emerged as a rapidly-growing field. Due to the chemical and functional diversity of lipids, the application of a systems biology approach is essential if one is to address lipid functionality at different physiological levels. In parallel with analytical advances to measure lipids in biological matrices, the field of computational lipidomics has been rapidly advancing, enabling modeling of lipidomes in their pathway, spatial and dynamic contexts. This review focuses on recent progress in systems biology approaches to study lipids in health and disease, with specific emphasis on methodological advances and biomedical applications.

Partitioning of Benzoic Acid into 1,2-Dimyristoyl-sn-glycero-3-phosphocholine and Blood-Brain Barrier Mimetic Bilayers

Using an all-atom explicit water model and replica exchange umbrella sampling simulations, we investigated the molecular mechanisms of benzoic acid partitioning into two model lipid bilayers. The first was formed of 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) lipids, whereas the second was composed of an equimolar mixture of DMPC, 1,2-dimyristoyl-sn-glycero-3-phosphoethanolamine, palmitoylsphingomyelin, and cholesterol to constitute a blood-brain barrier (BBB) mimetic bilayer. Comparative analysis of benzoic acid partitioning into the two bilayers has revealed qualitative similarities. Partitioning into the DMPC and BBB bilayers is thermodynamically favorable although insertion into the former lowers the free energy of benzoic acid by approximately an additional 1 kcal mol^-1. The partitioning energetics for the two bilayers is also largely similar based on the balance of benzoic acid interactions with apolar fatty acid tails, polar lipid headgroups, and water. In both bilayers, benzoic acid retains a considerable number of residual water molecules until reaching the bilayer midplane where it experiences nearly complete dehydration. Upon insertion into the bilayers, benzoic acid undergoes several rotations primarily determined by the interactions with the lipid headgroups. Nonetheless, in addition to the depth of the free energy minimum, the BBB bilayer differs from the DMPC counterpart by a much deeper location of the free energy minimum and the appearance of a high free energy barrier and positioning of benzoic acid near the midplane. Furthermore, DMPC and BBB bilayers exhibit different structural responses to benzoic acid insertion. Taken together, the BBB mimetic bilayer is preferable for an accurate description of benzoic acid partitioning.

Enhanced water permeability across a physiological droplet interface bilayer doped with fullerenes

We measure the water permeability across a physiological lipid bilayer produced by the droplet interface bilayer (DiB) technique. This lipid bilayer can be considered as physiologically relevant because it presents a lipidic composition close to human cell membranes. The measured water permeability coefficients across this lipid bilayer are reported as a function of the cholesterol concentration. It is found that the water permeability coefficients decreased with increasing cholesterol concentration, in agreement with the existing literature. And, consistently, the extracted corresponding activation energies increase with increasing cholesterol concentration in the lipid bilayer. Hence having demonstrated the robustness of the experimental system, we extend this study by exploring the influence of fullerenes on the water permeability of a physiological lipid bilayer. Interestingly, we observe a significant increase of the measured water permeability coefficients across this lipid bilayer for large fullerenes concentration. This enhanced permeability might be related to the conductive properties of fullerenes.

We measure the water permeability across a physiological lipid bilayer produced by the droplet interface bilayer technique.

Molecular dynamics study of solubilization of cyclohexane, benzene, and phenol into mixed micelles composed of sodium dodecyl sulfate and octaethylene glycol monododecyl ether

Molecular dynamics calculations of a mixed micelle composed of sodium dodecyl sulfate (SDS) and octaethylene glycol monododecyl ether (C₁₂ E₈ ) were performed for six compositions (SDS/C₁₂ E₈ = 100/0, 80/20, 60/40, 40/60, 20/80, and 0/100) to investigate the composition dependence of the mixed micelle structure and solubilization of cyclohexane, benzene, and phenol molecules by the micelle. The radial density distribution of the hydrophilic polyoxyethylene (POE) group of C₁₂ E₈ as a function of distance from the micelle center is very sharp for micelles with high SDS content because the POE group captures a Na⁺ ion in solution and wraps around it to form a compact crown-ether-like complex. The hydrophobic dodecyl groups of SDS and C₁₂ E₈ were separately distributed in the mixed micelle core. ΔG(r) evaluated for each solute showed that despite the structural changes of the micelle the binding strength of the solute molecules to the micelle did not change significantly. © 2019 Wiley Periodicals, Inc.

Interaction of Antimicrobial Lipopeptides with Bacterial Lipid Bilayers

The resistance of pathogens to traditional antibiotics is currently a global issue of enormous concern. As the discovery and development of new antibiotics become increasingly challenging, synthetic antimicrobial lipopeptides (AMLPs) are now receiving renewed attention as a new class of antimicrobial agents. In contrast to traditional antibiotics, AMLPs act by physically disrupting the cell membrane (rather than targeting specific proteins), thus reducing the risk of inducing bacterial resistance. In this study, we use microsecond-timescale atomistic molecular dynamics simulations to quantify the interaction of a short AMLP (C16-KKK) with model bacterial lipid bilayers. In particular, we investigate how fundamental transmembrane properties change in relation to a range of lipopeptide concentrations. A number of structural, mechanical, and dynamical features are found to be significantly altered in a non-linear fashion. At 10 mol% concentration, lipopeptides have a condensing effect on bacterial bilayers, characterized by a decrease in the area per lipid and an increase in the bilayer order. Higher AMLP concentrations of 25 and 40 mol% destabilize the membrane by disrupting the bilayer core structure, inducing membrane thinning and water leakage. Important transmembrane properties such as the lateral pressure and dipole potential profiles are also affected. Potential implications on membrane function and associated proteins are discussed.

Role of lipids in the thermal plasticity of basidial fungus Favolaschia manipularis

In this study, we examined the lipid composition of two strains of the tropical basidiomycete Favolaschia manipularis (Berk.) Teng, which differ in their adaptive potential to high (35 °C) and low (5 °C) temperatures. The results suggest that adaptation to extreme temperatures involves a change in the molecular composition of sterols, in addition to other well-known mechanisms of regulating membrane thickness and fluidity, such as changes in the lipid unsaturation and in the proportion of bilayer- and non-bilayer-forming lipids. It was demonstrated for the first time that adaptation to high temperature stress in fungi is accompanied by the accumulation of 9(11)-dehydroergosterol and ergosterol peroxide. Furthermore, increased thermal plasticity correlates with high storage lipid (triglycerides) content, accumulation of phosphatidic acid in the membrane, and an equal proportion of bilayer and non-bilayer lipids in the membrane.

Effects of Cholesterol on Water Permittivity of Biomimetic Ion Pair Amphiphile Bilayers: Interplay between Membrane Bending and Molecular Packing

Ion pair amphiphile (IPA), a molecular complex composed of a pair of cationic and anionic amphiphiles, is an inexpensive phospholipid substitute to fabricate vesicles with various pharmaceutical applications. Modulating the physicochemical and permeation properties of IPA vesicles are important for carrier designs. Here, we applied molecular dynamics simulations to examine the cholesterol effects on the structures, mechanics, and water permittivity of hexadecyltrimethylammonium-dodecylsulfate (HTMA-DS) and dodecyltrimethylammonium- hexadecylsulfate (DTMA-HS) IPA bilayers. Structural and mechanical analyses indicate that both IPA systems are in gel phase at 298 K. Adding cholesterol induces alkyl chain ordering around the rigid sterol ring and increases the cavity density within the hydrophilic region of both IPA bilayers. Furthermore, the enhanced alkyl chain ordering and the membrane deformation energy induced by cholesterol increase the permeation free energy penalty. In contrast, cholesterol has minor effects on the water local diffusivities within IPA membranes. Overall, the cholesterol reduces the water permittivity of rigid IPA membranes due to the synergistic effects of increased alkyl chain ordering and enhanced membrane mechanical modulus. The results provide molecular insights into the effects of molecular packing and mechanical deformations on the water permittivity of biomimetic IPA membranes, which is critical for designing IPA vesicular carriers.

PerMM: A Web Tool and Database for Analysis of Passive Membrane Permeability and Translocation Pathways of Bioactive Molecules

The PerMM web server and database were developed for quantitative analysis and visualization of passive translocation of bioactive molecules across lipid membranes. The server is the first physics-based web tool that calculates membrane binding energies and permeability coefficients of diverse molecules through artificial and natural membranes (phospholipid bilayers, PAMPA-DS, blood-brain barrier, and Caco-2/MDCK cell membranes). It also visualizes the transmembrane translocation pathway as a sequence of translational and rotational positions of a permeant as it moves across the lipid bilayer, along with the corresponding changes in solvation energy. The server can be applied for prediction of permeability coefficients of compounds with diverse chemical scaffolds to facilitate selection and optimization of potential drug leads. The complementary PerMM database allows comparison of computationally and experimentally determined permeability coefficients for more than 500 compounds in different membrane systems. The website and database are freely accessible at https://permm.phar.umich.edu/ .

Physics-Based Method for Modeling Passive Membrane Permeability and Translocation Pathways of Bioactive Molecules

Assessment of permeability is a critical step in the drug development process for selection of drug candidates with favorable ADME properties. We have developed a novel physics-based method for fast computational modeling of passive permeation of diverse classes of molecules across lipid membranes. The method is based on heterogeneous solubility-diffusion theory and operates with all-atom 3D structures of solutes and the anisotropic solvent model of the lipid bilayer characterized by transbilayer profiles of dielectric and hydrogen bonding capacity parameters. The optimal translocation pathway of a solute is determined by moving an ensemble of representative conformations of the molecule through the dioleoyl-phosphatidylcholine (DOPC) bilayer and optimizing their rotational orientations in every point of the transmembrane trajectory. The method calculates (1) the membrane-bound state of the solute molecule; (2) free energy profile of the solute along the permeation pathway; and (3) the permeability coefficient obtained by integration over the transbilayer energy profile and assuming a constant size-dependent diffusivity along the membrane normal. The accuracy of the predictions was evaluated against experimental permeability coefficients measured in pure lipid membranes (for 78 compounds, R² was 0.88 and rmse was 1.15 log units), PAMPA-DS (for 280 compounds, R² was 0.75 and rmse was 1.59 log units), BBB (for 182 compounds, R² was 0.69 and rmse was 0.87 log units), and Caco-2/MDCK assays (for 165 compounds, R² was 0.52 and rmse was 0.89 log units).

Glycyrrhizin-Assisted Transport of Praziquantel Anthelmintic Drug through the Lipid Membrane: An Experiment and MD Simulation

Praziquantel (PZQ) is one of the most widespread anthelmintic drugs. However, the frequent insufficient application of PZQ after oral administration is associated with its low solubility, penetration rate, and bioavailability. In the present study, the permeation of PZQ through a 1,2-dioleoyl- sn-glycero-3-phosphocholine (DOPC) membrane was investigated to probe glycyrrhizin-assisted transport. Glycyrrhizin (or glycyrrhizic acid, GA), a natural saponin, shows the ability to enhance the therapeutic activity of various drugs when it is used as a drug delivery system. However, the molecular mechanism of this effect is still under debate. In the present study, the transport rate was measured experimentally by a parallel artificial membrane permeation assay (PAMPA) and molecular dynamics (MD) simulation with DOPC lipid bilayers. The formation of the noncovalent supramolecular complex of PZQ with disodium salt of GA (Na2GA) in an aqueous solution was proved by the NMR relaxation technique. PAMPA experiments show a strong increase in the amount of the penetrating praziquantel molecules in comparison with a saturated aqueous solution of pure drug used as a control. MD simulation of PZQ penetration through the bilayer demonstrates an increase in permeability into the membrane in the presence of a glycyrrhizin molecule. A decrease in the free energy barrier in the middle of the lipid bilayer was obtained, associated with the hydrogen bond between PZQ and GA. Also, GA reduces the local bilayer surface resistance to penetration of PZQ by rearranging the surface lipid headgroups. This study clarifies the mechanism of increasing the drug's bioavailability in the presence of glycyrrhizin.

Effects of lipid composition on membrane distribution and permeability of natural quinones

Natural quinones are amphiphilic molecules that function as mobile charge carriers in biological energy transduction. Their distribution and permeation across membranes are important for binding to enzymatic complexes and for proton translocation. Here, we employ molecular dynamics simulations and free energy calculations with a carefully calibrated classical force-field to probe quinone distribution and permeation in a multi-component bilayer trying to mimic the composition of membranes involved in bioenergetic processes. Ubiquinone, ubiquinol, plastoquinone and menaquinone molecules with short and long isoprenoid tails are simulated. We find that penetration of water molecules bound to the polar quinone head increases considerably in the less ordered and porous bilayer formed by di-linoleoyl (18:2) phospholipids, resulting in a lower free energy barrier for quinone permeation and faster transversal diffusion. In equilibrium, quinone and quinol heads localize preferentially near lipid glycerol groups, but do not perform specific contacts with lipid polar heads. Quinone distribution is not altered significantly by the quinone head, tail and lipid composition in comparison to a single-component bilayer. This study highlights the role of lipid acyl chain unsaturation for permeation and transversal diffusion of polar molecules across biological membranes.

Molecular Dynamics Simulations of Membrane Permeability

This Review illustrates the evaluation of permeability of lipid membranes from molecular dynamics (MD) simulation primarily using water and oxygen as examples. Membrane entrance, translocation, and exit of these simple permeants (one hydrophilic and one hydrophobic) can be simulated by conventional MD, and permeabilities can be evaluated directly by Fick's First Law, transition rates, and a global Bayesian analysis of the inhomogeneous solubility-diffusion model. The assorted results, many of which are applicable to simulations of nonbiological membranes, highlight the limitations of the homogeneous solubility diffusion model; support the utility of inhomogeneous solubility diffusion and compartmental models; underscore the need for comparison with experiment for both simple solvent systems (such as water/hexadecane) and well-characterized membranes; and demonstrate the need for microsecond simulations for even simple permeants like water and oxygen. Undulations, subdiffusion, fractional viscosity dependence, periodic boundary conditions, and recent developments in the field are also discussed. Last, while enhanced sampling methods and increasingly sophisticated treatments of diffusion add substantially to the repertoire of simulation-based approaches, they do not address directly the critical need for force fields with polarizability and multipoles, and constant pH methods.

Structural and barrier properties of the skin ceramide lipid bilayer: a molecular dynamics simulation study

Skin provides excellent protection against the harsh external environment and foreign substances. The lipid matrix of the stratum corneum, which contains various kinds of ceramides, plays a major role in the barrier function of the skin. Here we report a study of the effects of ceramide type on the structural and transport properties of ceramide bilayers using molecular dynamics (MD) simulations. Specifically, the effects of headgroup chemistry (number and positions of hydroxyl groups) and tail structure (unsaturation of the sphingoid moiety) on the structural and transport properties of various ceramide bilayers at 310 K were analyzed. Theoretical results for structural properties such as area per lipid, bilayer thickness, lateral arrangement, order parameter, and hydrogen bonding are reported here and compared with corresponding experimental data. Our study revealed that the presence of a double bond disrupts the bilayer packing, which leads to a low area compressibility modulus, a large area per lipid, and low bilayer thickness. Furthermore, the effect of structural changes on water permeation was studied using steered MD simulations. Water permeation was found to be influenced by headgroup polarity, chain packing, and the ability of the water to hydrogen bond with the ceramides. The molecular-level information obtained from the current study should aid the design of mixed bilayer systems with desired properties and provide the basis for the development of higher order coarse-grained models.