Behaviour of small solutes and large drugs in a lipid bilayer from computer simulations

Examine stability polyvinyl alcohol-stabilized nanosuspensions to overcome the challenge of poor drug solubility utilizing molecular dynamic simulation

The pharmaceutical industry faces a significant challenge from the low water solubility of nearly 90% of newly developed Active Pharmaceutical Ingredients (APIs). Despite extensive efforts to improve solubility, approximately 40% of these APIs encounter commercialization hurdles, impacting drug efficacy. In this context, a promising strategy will be introduced in which nanosuspensions, particularly polyvinyl alcohol (PVA) as a stabilizer, are applied to increase drug solubility. In this work using molecular dynamics simulations, the nanosuspension of four poorly water-soluble drugs (flurbiprofen, bezafibrate, miconazole, and phenytoin) stabilized with PVA is investigated. The simulation data showed van der Waals energies between polyvinyl alcohol with flurbiprofen and bezafibrate are - 101.12 and - 58.42 kJ/mol, respectively. The results indicate that PVA is an effective stabilizer for these drugs, and superior interactions are obtained with flurbiprofen and bezafibrate. The study also explores the impact of PVA on water molecule diffusion, providing insights into the stability of nanosuspensions. Obtained results also provide valuable insights into hydrogen bond formation, diffusion coefficients, and nanosuspension stability, contributing to the rational design and optimization of pharmaceutical formulations.

Potential of liposomes and lipid membranes for the separation of β-blockers by capillary electromigration and liquid chromatographic techniques

β-Blockers belong to a frequently used class of drugs primarily used to treat heart and circulatory conditions. Here we describe the use of lipid vesicles and liposomes as cell membrane biomimicking models in capillary electromigration (CE) and liquid chromatography (LC) techniques for the investigation of interactions between lipid membranes and β-blockers. In addition to liposomes, the use of commercial intravenous lipid emulsions, and their interactions with β-blockers are also discussed. Different CE and LC instrumental techniques designed for these purposes are introduced. Other methodologies for studying interactions between β-blockers and lipid membranes are also briefly discussed, and the different methodologies are compared. The aim is to give the reader a good overview on the status of the use of liposomes and lipids in CE and LC for studying β-blocker interactions.

Effect of liposome composition on β-blocker interactions studied by capillary electrokinetic chromatography

Liposome capillary electrokinetic chromatography was used to investigate the interactions between three β-blockers of different hydrophobicity and various liposome solutions. The studied β-blockers comprised alprenolol, propranolol, and carvedilol. The composition of the liposome solutions, containing 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine, 1,2-dioleoyl-sn-glycero-3-phos-phoethanolamine, 1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-l-serine, and cholesterol in various molar ratios, was designed by a response surface methodology-central composite design approach. Subsequently, after conducting the liposome capillary electrokinetic chromatography experiments and determining the retention factors from the electrophoretic mobilities of the compounds, and further calculating the distribution coefficients, an analysis of variance was performed. After extracting the statistical models, optimal operational conditions were obtained based on the developed models. To further investigate the interactions between the β-blockers and the liposomes, nanoplasmonic sensing experiments were carried out on two different liposome systems. The overall results demonstrate the strong influence of cholesterol and 1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-l-serine on the distribution coefficients.

Deciphering Ethanol-Driven Swelling, Rupturing, Aggregation, and Fusion of Lipid Vesicles Using Coarse-Grained Molecular Dynamics Simulations

Traditionally, liquid ethanol is known to enhance the permeability of lipid membranes and causes vesicle aggregation and fusion. However, how the amphiphilic ethanol molecules perturb the lipid vesicles to facilitate their aggregation or fusion has not been addressed at any level of molecular simulations. Herein, not only have we developed a coarse-grained (CG) model for liquid ethanol, its aqueous mixture, and hydrated lipid membranes for molecular dynamics (MD) simulations, but also utilized it to delineate the aggregation and fusion of lipid vesicles using CG-MD simulations with multimillion particles. We have systematically parametrized the force-field for pure ethanol and its interactions with hydrated POPC and POPE model lipid membranes. In this process, we have successfully reproduced the bulk ethanol structure and concentration-dependent density of aqueous ethanol. To quantify the interaction of ethanol with lipid membranes, we have reproduced the transfer free energy of the ethanol molecule across the hydrated bilayers, and the concentration-dependent distribution of ethanol molecules across the lipid bilayers. After having acceptable force-field parameters for ethanol-membrane interactions, we have checked the effect of ethanol toward the vesicles comprising POPC lipids. We observe a rapid increase in the size of the POPC lipid vesicles with increasing amounts of ethanol up to 30 mol %. We unambiguously observe swelling and decrease in the thickness of the POPC vesicles with increasing amounts of ethanol up to 30 mol %, beyond which the vesicles begin to lose their integrity and rupture at higher mol % of ethanol. The fusion study of two vesicles demonstrates that fused vesicles can be obtained from 20 to 30 mol % of ethanol provided that they are brought closer than a critical distance at a particular mol %. The multivesicle simulations show that along with the increase in the sizes of vesicles the propensity of vesicle aggregation increases as the mol % of ethanol increases.

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.

The role of the envelope protein in the stability of a coronavirus model membrane against an ethanolic disinfectant

Ethanol is highly effective against various enveloped viruses and can disable the virus by disintegrating the protective envelope surrounding it. The interactions between the coronavirus envelope (E) protein and its membrane environment play key roles in the stability and function of the viral envelope. By using molecular dynamics simulation, we explore the underlying mechanism of ethanol-induced disruption of a model coronavirus membrane and, in detail, interactions of the E-protein and lipids. We model the membrane bilayer as N-palmitoyl-sphingomyelin and 1-palmitoyl-2-oleoylphosphatidylcholine lipids and the coronavirus E-protein. The study reveals that ethanol causes an increase in the lateral area of the bilayer along with thinning of the bilayer membrane and orientational disordering of lipid tails. Ethanol resides at the head-tail region of the membrane and enhances bilayer permeability. We found an envelope-protein-mediated increase in the ordering of lipid tails. Our simulations also provide important insights into the orientation of the envelope protein in a model membrane environment. At ∼25 mol. % of ethanol in the surrounding ethanol-water phase, we observe disintegration of the lipid bilayer and dislocation of the E-protein from the membrane environment.

Functional Group Distributions, Partition Coefficients, and Resistance Factors in Lipid Bilayers Using Site Identification by Ligand Competitive Saturation

Small molecules such as metabolites and drugs must pass through the membrane of the cell, a barrier primarily comprising phospholipid bilayers and embedded proteins. To better understand the process of passive diffusion, knowledge of the ability of various functional groups to partition across bilayers and the associated energetics would be of utility. In the present study, the site identification by ligand competitive saturation (SILCS) methodology has been applied to sample the distributions of a diverse set of chemical solutes representing the functional groups of small molecules across phospholipid bilayers composed of 0.9:0.1 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine/cholesterol and a mixture of 0.52:0.18:0.3 1,2-dioleoyl-sn-glycero-3-phospho-l-serine/1,2-dioleoyl-sn-glycero-3-phosphocholine/cholesterol used in parallel artificial membrane permeability assay experiments. A combination of oscillating chemical potential grand canonical Monte Carlo and molecular dynamics in the SILCS simulations was applied to achieve solute sampling through the bilayers and surrounding aqueous environment from which the distribution of solutes and the functional groups they represent were obtained. Results show differential distribution of aliphatic versus aromatic groups with the former having increased sampling in the center of the bilayers versus in the region of the glycerol linker for the latter. Variations in the distribution of different polar groups are evident, with large differences between negative acetate and positive methylammonium with accumulation of the polar-neutral and acetate solutes above the bilayer head groups. Conversion of the distributions to absolute free energies allows for a detailed understanding of energetics of functional groups in different regions of the bilayers and for calculation of absolute free-energy profiles of multifunctional drug-like molecules across the bilayers from which partition coefficients and resistance factors suitable for insertion into the homogenous solubility-diffusion equation for calculation of permeability were obtained. Comparisons of the calculated bilayer/solution partition coefficients with 1-octanol/water experimental data for both drug-like molecules and the solutes show overall good agreement, validating the calculated distributions and associated absolute free-energy profiles.

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.

Assessing Barriers for Antimicrobial Penetration in Complex Asymmetric Bacterial Membranes: A Case Study with Thymol

The bacterial cell envelope is a complex multilayered structure evolved to protect bacteria in hostile environments. An understanding of the molecular basis for the interaction and transport of antibacterial therapeutics with the bacterial cell envelope will enable the development of drug molecules to combat bacterial infections and suppress the emergence of drug-resistant strains. Here we report the successful creation of an in vitro supported lipid bilayer (SLB) platform of the outer membrane (OM) of E. coli, an archetypical Gram-negative bacterium, containing the full smooth lipopolysaccharide (S-LPS) architecture of the membrane. Using this platform, we performed fluorescence correlation spectroscopy (FCS) in combination with molecular dynamics (MD) simulations to measure lipid diffusivities and provide molecular insights into the transport of natural antimicrobial agent thymol. Lipid diffusivities measured on symmetric supported lipid bilayers made up of inner membrane lipids show a distinct increase in the presence of thymol as also corroborated by MD simulations. However, lipid diffusivities in the asymmetric OM consisting of only S-LPS are invariant upon exposure to thymol. Increasing the phospholipid content in the LPS-containing outer leaflet improved the penetration toward thymol as reflected in slightly higher relative diffusivity changes in the inner leaflet when compared with the outer leaflet. Free-energy computations reveal the presence of a barrier (∼6 kT) only in the core-saccharide region of the OM for the translocation of thymol while the external O-antigen part is easily traversed. In contrast, thymol spontaneously inserts into the inner membrane. In addition to providing leaflet-resolved penetration barriers in bacterial membranes, we also assess the ability of small molecules to penetrate various membrane components. With rising bacterial resistance, our study opens up the possibility of screening potential antimicrobial drug candidates using these realistic model platforms for Gram-negative bacteria.

Molecular Dynamics Simulation of Small Molecules Interacting with Biological Membranes

Cell membranes protect and compartmentalise cells and their organelles. The semi-permeable nature of these membranes controls the exchange of solutes across their structure. Characterising the interaction of small molecules with biological membranes is critical to understanding of physiological processes, drug action and permeation, and many biotechnological applications. This review provides an overview of how molecular simulations are used to study the interaction of small molecules with biological membranes, with a particular focus on the interactions of water, organic compounds, drugs and short peptides with models of plasma cell membrane and stratum corneum lipid bilayers. This review will not delve on other types of membranes which might have different composition and arrangement, such as thylakoid or mitochondrial membranes. The application of unbiased molecular dynamics simulations and enhanced sampling methods such as umbrella sampling, metadynamics and replica exchange are described using key examples. This review demonstrates how state-of-the-art molecular simulations have been used successfully to describe the mechanism of binding and permeation of small molecules with biological membranes, as well as associated changes to the structure and dynamics of these membranes. The review concludes with an outlook on future directions in this field.

Model-Informed Drug Discovery and Development Strategy for the Rapid Development of Anti-Tuberculosis Drug Combinations

The increasing emergence of drug-resistant tuberculosis requires new effective and safe drug regimens. However, drug discovery and development are challenging, lengthy and costly. The framework of model-informed drug discovery and development (MID3) is proposed to be applied throughout the preclinical to clinical phases to provide an informative prediction of drug exposure and efficacy in humans in order to select novel anti-tuberculosis drug combinations. The MID3 includes pharmacokinetic-pharmacodynamic and quantitative systems pharmacology models, machine learning and artificial intelligence, which integrates all the available knowledge related to disease and the compounds. A translational in vitro-in vivo link throughout modeling and simulation is crucial to optimize the selection of regimens with the highest probability of receiving approval from regulatory authorities. In vitro-in vivo correlation (IVIVC) and physiologically-based pharmacokinetic modeling provide powerful tools to predict pharmacokinetic drug-drug interactions based on preclinical information. Mechanistic or semi-mechanistic pharmacokinetic-pharmacodynamic models have been successfully applied to predict the clinical exposure-response profile for anti-tuberculosis drugs using preclinical data. Potential pharmacodynamic drug-drug interactions can be predicted from in vitro data through IVIVC and pharmacokinetic-pharmacodynamic modeling accounting for translational factors. It is essential for academic and industrial drug developers to collaborate across disciplines to realize the huge potential of MID3.

Intestinal Permeability and Drug Absorption: Predictive Experimental, Computational and In Vivo Approaches

The main objective of this review is to discuss recent advancements in the overall investigation and in vivo prediction of drug absorption. The intestinal permeability of an orally administered drug (given the value P_eff) has been widely used to determine the rate and extent of the drug’s intestinal absorption (F_abs) in humans. Preclinical gastrointestinal (GI) absorption models are currently in demand for the pharmaceutical development of novel dosage forms and new drug products. However, there is a strong need to improve our understanding of the interplay between pharmaceutical, biopharmaceutical, biochemical, and physiological factors when predicting F_abs and bioavailability. Currently, our knowledge of GI secretion, GI motility, and regional intestinal permeability, in both healthy subjects and patients with GI diseases, is limited by the relative inaccessibility of some intestinal segments of the human GI tract. In particular, our understanding of the complex and highly dynamic physiology of the region from the mid-jejunum to the sigmoid colon could be significantly improved. One approach to the assessment of intestinal permeability is to use animal models that allow these intestinal regions to be investigated in detail and then to compare the results with those from simple human permeability models such as cell cultures. Investigation of intestinal drug permeation processes is a crucial biopharmaceutical step in the development of oral pharmaceutical products. The determination of the intestinal P_eff for a specific drug is dependent on the technique, model, and conditions applied, and is influenced by multiple interactions between the drug molecule and the biological membranes.

Effects of absorption-modifying excipients on jejunal drug absorption in simulated fasted and fed luminal conditions

Oral administration of drug products is the preferred administration route. In recent decades there has been an increase in drug candidates with low solubility and/or low permeability. To increase the possibility of oral administration for the poorly permeating drugs, the use of absorption modifying excipients (AMEs) has been proposed. These types of AMEs may also affect the regulatory assessment of a novel drug delivery system if they affect the absorption of a drug from any of the four BCS classes. The effects of AMEs have previously been investigated in various animal models, including the single-pass intestinal perfusion (SPIP) in rats. To further improve the biorelevance and the in vivo predictiveness of the SPIP model, four compounds (atenolol, enalaprilat, ketoprofen, metoprolol) were perfused in fasted or fed state simulated intestinal fluid (FaSSIF or FeSSIF) together with the AMEs N-acetyl-cysteine, caprate, or sodium dodecyl sulfate. For the highly soluble and poorly permeating compounds enalaprilat and atenolol (BCS class III), the flux was increased the most by the addition of SDS in both FaSSIF and FeSSIF. For ketoprofen (BCS class II), the flux decreased in the presence of all AMEs in at least one of the perfusion media. The flux of metoprolol (BCS class I) was not affected by any of the excipients in none of simulated prandial states. The changes in magnitude in the absorption of the compounds were in general smaller in FeSSIF than in FaSSIF. This may be explained by a reduced free concentration AMEs in FeSSIF. Further, the results in FeSSIF were similar to those from intrajejunal bolus administration in rat in a previous study. This suggests that the biorelevance of the SPIP method may be increased when investigating the effects of AMEs, by the addition of intraluminal constituents representative to fasted and/or fed state to the inlet perfusate.

Does the Lipid Bilayer Orchestrate Access and Binding of Ligands to Transmembrane Orthosteric/Allosteric Sites of G Protein-Coupled Receptors?

The ligand-binding sites of many G protein-coupled receptors (GPCRs) are situated around and deeply embedded within the central pocket formed by their seven transmembrane-spanning α-helical domains. Generally, these binding sites are assumed accessible to endogenous ligands from the aqueous phase. Recent advances in the structural biology of GPCRs, along with biophysical and computational studies, suggest that amphiphilic and lipophilic molecules may gain access to these receptors by first partitioning into the membrane and then reaching the binding site via lateral diffusion through the lipid bilayer. In addition, several crystal structures of class A and class B GPCRs bound to their ligands offer unprecedented details on the existence of lipid-facing allosteric binding sites outside the transmembrane helices that can only be reached via lipid pathways. The highly organized structure of the lipid bilayer may direct lipophilic or amphiphilic drugs to a specific depth within the bilayer, changing local concentration of the drug near the binding site and affecting its binding kinetics. Additionally, the constraints of the lipid bilayer, including its composition and biophysical properties, may play a critical role in "pre-organizing" ligand molecules in an optimal orientation and conformation to facilitate receptor binding. Despite its clear involvement in molecular recognition processes, the critical role of the membrane in binding ligands to lipid-exposed transmembrane binding sites remains poorly understood and warrants comprehensive investigation. Understanding the mechanistic basis of the structure-membrane interaction relationship of drugs will not only provide useful insights about receptor binding kinetics but will also enhance our ability to take advantage of the apparent membrane contributions when designing drugs that target transmembrane proteins with improved efficacy and safety. In this minireview, we summarize recent structural and computational studies on membrane contributions to binding processes, elucidating both lipid pathways of ligand access and binding mechanisms for several orthosteric and allosteric ligands of class A and class B GPCRs.

Effects of lipid composition on membrane permeation

Passive permeation through lipid membranes is an essential process in biology. In vivo membranes typically consist of mixtures of lamellar and nonlamellar lipids. Lamellar lipids are characterized by their tendency to form lamellar sheet-like structures, which are predominant in nature. Nonlamellar lipids, when isolated, instead form more geometrically complex nonlamellar phases. While mixed lamellar/nonlamellar lipid membranes tend to adopt the ubiquitous lamellar bilayer structure, the presence of nonlamellar lipids is known to have profound effects on key membrane properties, such as internal distributions of stress and elastic properties, which in turn may alter related biological processes. This work focuses on one such process, i.e., permeation, by utilising atomistic molecular dynamics simulations in order to obtain transfer free energy profiles, diffusion profiles and permeation coefficients for a series of thirteen small molecules and drugs. Each permeant is tested on two bilayer membranes of different lipid composition, i.e., purely lamellar and mixed lamellar/nonlamellar. Our results indicate that the presence of nonlamellar lipids reduces permeation for smaller molecules (molecular weight < 100) but facilitates it for the largest ones (molecular weight > 100). This work represents an advancement towards the development of more realistic in silico permeability assays, which may have a substantial future impact in the area of rational drug design.

Molecular transport through membranes: Accurate permeability coefficients from multidimensional potentials of mean force and local diffusion constants

Estimating the permeability coefficient of small molecules through lipid bilayer membranes plays an important role in the development of effective drug candidates. In silico simulations can produce acceptable relative permeability coefficients for a series of small molecules; however, the absolute permeability coefficients from simulations are usually off by orders of magnitude. In addition to differences between the lipid bilayers used in vitro and in silico, the poor convergence of permeation free energy profiles and over-simplified diffusion models have contributed to these discrepancies. In this paper, we present a multidimensional inhomogeneous solubility-diffusion model to study the permeability of a small molecule drug (trimethoprim) passing through a POPC (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) lipid bilayer. Our approach improves the permeation model in three ways: First, the free energy profile (potential of mean force, PMF) is two-dimensional in two key coordinates rather than simply one-dimensional along the direction normal to the bilayer. Second, the 2-D PMF calculation has improved convergence due to application of the recently developed transition-tempered metadynamics with randomly initialized replicas, while third, the local diffusivity coefficient was calculated along the direction of the minimum free energy path on the two-dimensional PMF. The permeability is then calculated as a line integral along the minimum free energy path of the PMF. With this approach, we report a considerably more accurate permeability coefficient (only 2-5 times larger than the experimental value). We also compare our approach with the common practice of computing permeability coefficients based only on the translation of the center of mass of the drug molecule. Our paper concludes with a discussion of approaches for minimizing the computational cost for the purpose of more rapidly screening a large number of drug candidate molecules.

Atomistic Picture of Fluorescent Probes with Hydrocarbon Tails in Lipid Bilayer Membranes: An Investigation of Selective Affinities and Fluorescent Anisotropies in Different Environmental Phases

By reverting to spectroscopy, changes in the biological environment of a fluorescent probe can be monitored and the presence of various phases of the surrounding lipid bilayer membranes can be detected. However, it is currently not always clear in which phase the probe resides. The well-known orange 1,1'-dioctadecyl-3,3,3',3'-tetramethylindodicarbo-cyanine perchlorate (DiI-C18(5)) fluorophore, for instance, and the new, blue BODIPY (4,4-difluoro-4-bora-3 a,4 a-diaza- s-indacene) derivative were experimentally seen to target and highlight identical parts of giant unilamellar vesicles of various compositions, comprising mixtures of dipalmitoylphosphatidylcholine (DPPC), dioleoylphosphatidylcholine (DOPC), sphingomyelin (SM), and cholesterol (Chol). However, it was not clear which of the coexisting membrane phases were visualized (Bacalum et al., Langmuir. 2016, 32, 3495). The present study addresses this issue by utilizing large-scale molecular dynamics simulations and the z-constraint method, which allows evaluating Gibbs free-energy profiles. The current calculations give an indication why, at room temperature, both BODIPY and DiI-C18(5) probes prefer the gel (S_o) phase in DOPC/DPPC (2:3 molar ratio) and the liquid-ordered (L_o) phase in DOPC/SM/Chol (1:2:1 molar ratio) mixtures. This study highlights the important differences in orientation and location and therefore in efficiency between the probes when they are used in fluorescence microscopy to screen various lipid bilayer membrane phases. Dependent on the lipid composition, the angle between the transition-state dipole moments of both probes and the normal to the membrane is found to deviate clearly from 90°. It is seen that the DiI-C18(5) probe is located in the headgroup region of the SM/Chol mixture, in close contact with water molecules. A fluorescence anisotropy study also indicates that DiI-C18(5) gives rise to a distinctive behavior in the SM/Chol membrane compared to the other considered membranes. The latter behavior has not been seen for the studied BODIPY probe, which is located deeper in the membrane.

Formulation and Molecular Dynamics Simulations of a Fusidic Acid Nanosuspension for Simultaneously Enhancing Solubility and Antibacterial Activity

The aim of the present study was to formulate a nanosuspension (FA-NS) of fusidic acid (FA) to enhance its aqueous solubility and antibacterial activity. The nanosuspension was characterized using various in vitro, in silico, and in vivo techniques. The size, polydispersity index, and zeta potential of the optimized FA-NS were 265 ± 2.25 nm, 0.158 ± 0.026, and -16.9 ± 0.794 mV, respectively. The molecular dynamics simulation of FA and Poloxamer-188 showed an interaction and binding energy of -74.42 kJ/mol and -49.764 ± 1.298 kJ/mol, respectively, with van der Waals interactions playing a major role in the spontaneous binding. There was an 8-fold increase in the solubility of FA in a nanosuspension compared to the bare drug. The MTT assays showed a cell viability of 75-100% confirming the nontoxic nature of FA-NS. In vitro antibacterial activity revealed a 16- and 18-fold enhanced activity against Staphylococcus aureus (SA) and methicillin-resistant SA (MRSA), respectively, when compared to bare FA. Flowcytometry showed that MRSA cells treated with FA-NS had almost twice the percentage of dead bacteria in the population, despite having an 8-fold lower MIC in comparison to the bare drug. The in vivo skin-infected mice showed a 76-fold reduction in the MRSA load for the FA-NS treated group compared to that of the bare FA. These results show that the nanosuspension of antibiotics can enhance their solubility and antibacterial activity simultaneously.

Composition Dependence of Water Permeation Across Multicomponent Gel-Phase Bilayers

The permeability of multicomponent phospholipid bilayers in the gel phase is investigated via molecular dynamics simulation. The physical role of the different molecules is probed by comparing multiple mixed-component bilayers containing distearylphosphatidylcholine (DSPC) with varying amounts of either the emollient isostearyl isostearate or long-chain alcohol (dodecanol, octadecanol, or tetracosanol) molecules. Permeability is found to depend on both the tail packing density and hydrogen bonding between lipid headgroups and water. Whereas the addition of emollient or alcohol molecules to a gel-phase DSPC bilayer can increase the tail packing density, it also disturbed the hydrogen-bonding network, which in turn can increase interfacial water dynamics. These phenomena have opposing effects on bilayer permeability, which is found to depend on the balance between enhanced tail packing and decreased hydrogen bonding.

Sequence-Dependent Interfacial Adsorption and Permeation of Dipeptides across Phospholipid Membranes

We investigate permeation of three blocked dipeptides with different side chain polarity across a phospholipid membrane and their behavior at the water-membrane interface by way of molecular dynamics simulations. Hydrophilic serine-serine dipeptide is found to desorb from the interface to aqueous phase, whereas hydrophobic phenylalanine-leucine and amphiphilic serine-leucine tend to accumulate at the interface with a free energy minimum of -3 kcal/mol. All three dipeptides exhibit free energy barriers to permeation across the membrane located at the center of the bilayer. The height of the barrier is strongly sequence dependent and increases with the dipeptide polarity. It is equal to 3.5, 6.4, and 10.0 kcal/mol for phenylalanine-leucine, serine-leucine, and serine-serine, respectively. The corresponding permeability coefficients are equal to 4.6 × 10^-3, 4.5 × 10^-5, and 8.7 × 10^-8 cm/s. The apparent insensitivity of membrane permeability to hydrophobicity of dipeptides, found in some experiments, is attributed to neglecting corrections for unstirred water layers near membrane surface, which are significant for hydrophobic species. Different hydrophobicity of the dipeptides also influences their conformations and orientations, both at the interface and inside the membrane. In particular, penetration of hydrophilic serine-serine dipeptide causes the formation of water-filled defects in the bilayer. These results are relevant to the delivery of peptide-based therapeutic agents.

Physical insights into the blood-brain barrier translocation mechanisms

The number of individuals suffering from diseases of the central nervous system (CNS) is growing with an aging population. While candidate drugs for many of these diseases are available, most of these pharmaceutical agents cannot reach the brain rendering most of the drug therapies that target the CNS inefficient. The reason is the blood-brain barrier (BBB), a complex and dynamic interface that controls the influx and efflux of substances through a number of different translocation mechanisms. Here, we present these mechanisms providing, also, the necessary background related to the morphology and various characteristics of the BBB. Moreover, we discuss various numerical and simulation approaches used to study the BBB, and possible future directions based on multi-scale methods. We anticipate that this review will motivate multi-disciplinary research on the BBB aiming at the design of effective drug therapies.

Effect of Size and Surface Charge of Gold Nanoparticles on their Skin Permeability: A Molecular Dynamics Study

Molecular level understanding of permeation of nanoparticles through human skin establishes the basis for development of novel transdermal drug delivery systems and design and formulation of cosmetics. Recent experiments suggest that surface coated nano-sized gold nanoparticles (AuNPs) can penetrate the rat and human skin. However, the mechanisms by which these AuNPs penetrate are not well understood. In this study, we have carried out coarse grained molecular dynamics simulations to explore the permeation of dodecanethiol coated neutral hydrophobic AuNPs of different sizes (2–5 nm) and surface charges (cationic and anionic) through the model skin lipid membrane. The results indicate that the neutral hydrophobic AuNPs disrupted the bilayer and entered in it with in ~200 ns, while charged AuNPs were adsorbed on the bilayer headgroup. The permeation free energy calculation revealed that at the head group of the bilayer, a very small barrier existed for neutral hydrophobic AuNP while a free energy minimum was observed for charged AuNPs. The permeability was maximum for neutral 2 nm gold nanoparticle (AuNP) and minimum for 3 nm cationic AuNP. The obtained results are aligned with recent experimental findings. This study would be helpful in designing customized nanoparticles for cosmetic and transdermal drug delivery application.

Structural effects and translocation of spontaneous membrane-translocating peptides with POPC bilayer

Martini coarse-grained force field simulations have been carried out to estimate the free energy profiles of the spontaneous membrane-translocating peptide TP2 and one negative control peptide ONEG with POPC as the model bilayer. The results show that the free energy minimum of TP2 is [Formula: see text]20[Formula: see text]kJ/mol lower than that of ONEG. In addition, the minimum of TP2 shifts slightly to the bilayer center compared with ONEG. The translocation barrier height for TP2 and ONEG are 119.0[Formula: see text]kJ/mol and 155.7[Formula: see text]kJ/mol, respectively. The lower central energy barrier of TP2 facilitates the transition between two leaflets of POPC. Both translocating peptides induce the formation of funnel-shaped structures at the bilayer center, but TP2 has a more compact structure and brings less perturbation compared with ONEG. Subsequently all atom molecular simulations testify the findings. It is indicated that compared with its negative control ONEG, TP2 binds better with lipid and penetrates deeper into bilayer with less perturbation to the bilayer structure. Our findings may shed light on the design and virtual screening of spontaneous membrane-translocating peptides.

Algorithm for Designing Nanoscale Supramolecular Therapeutics with Increased Anticancer Efficacy

In the chemical world, evolution is mirrored in the origin of nanoscale supramolecular structures from molecular subunits. The complexity of function acquired in a supramolecular system over a molecular subunit can be harnessed in the treatment of cancer. However, the design of supramolecular nanostructures is hindered by a limited atomistic level understanding of interactions between building blocks. Here, we report the development of a computational algorithm, which we term Volvox after the first multicellular organism, that sequentially integrates quantum mechanical energy-state- and force-field-based models with large-scale all-atomistic explicit water molecular dynamics simulations to design stable nanoscale lipidic supramolecular structures. In one example, we demonstrate that Volvox enables the design of a nanoscale taxane supramolecular therapeutic. In another example, we demonstrate that Volvox can be extended to optimizing the ratio of excipients to form a stable nanoscale supramolecular therapeutic. The nanoscale taxane supramolecular therapeutic exerts greater antitumor efficacy than a clinically used taxane in vivo. Volvox can emerge as a powerful tool in the design of nanoscale supramolecular therapeutics for effective treatment of cancer.

The Free Energy of Small Solute Permeation through the Escherichia coli Outer Membrane Has a Distinctly Asymmetric Profile

Permeation of small molecules across cell membranes is a ubiquitous process in biology and is dependent on the principles of physical chemistry at the molecular level. Here we use atomistic molecular dynamics simulations to calculate the free energy of permeation of a range of small molecules through a model of the outer membrane of Escherichia coli, an archetypical Gram-negative bacterium. The model membrane contains lipopolysaccharide (LPS) molecules in the outer leaflet and phospholipids in the inner leaflet. Our results show that the energetic barriers to permeation through the two leaflets of the membrane are distinctly asymmetric; the LPS headgroups provide a less energetically favorable environment for organic compounds than do phospholipids. In summary, we provide the first reported estimates of the relative free energies associated with the different chemical environments experienced by solutes as they attempt to cross the outer membrane of a Gram-negative bacterium. These results provide key insights for the development of novel antibiotics that target these bacteria.

Two Relations to Estimate Membrane Permeability Using Milestoning

Prediction of passive permeation rates of solutes across lipid bilayers is important to drug design, toxicology, and other biological processes such as signaling. The inhomogeneous solubility-diffusion (ISD) equation is traditionally used to relate the position-dependent potential of mean force and diffusivity to the permeability coefficient. The ISD equation is derived via the Smoluchowski equation and assumes overdamped system dynamics. It has been suggested that the complex membrane environment may exhibit more complicated damping conditions. Here we derive a variant of the inhomogeneous solubility diffusion equation as a function of the mean first passage time (MFPT) and show how milestoning, a method that can estimate kinetic quantities of interest, can be used to estimate the MFPT of membrane crossing and, by extension, the permeability coefficient. We further describe a second scheme, agnostic to the damping condition, to estimate the permeability coefficient from milestoning results or other methods that compute a probability of membrane crossing. The derived relationships are tested using a one-dimensional Langevin dynamics toy system confirming that the presented theoretical methods can be used to estimate permeabilities given simulation and milestoning results.

Study of procaine and tetracaine in the lipid bilayer using molecular dynamics simulation

Despite available experimental results, the molecular mechanism of action of local anesthetics upon the nervous system and contribution of the cell membrane to the process are still controversial. In this work, molecular dynamics simulations were performed to investigate the effect of two clinically used local anesthetics, procaine and tetracaine, on the structure and dynamics of a fully hydrated dimyristoylphosphatidylcholine lipid bilayer. We focused on comparing the main effects of uncharged and charged drugs on various properties of the lipid membrane: mass density distribution, diffusion coefficient, order parameter, radial distribution function, hydrogen bonding, electrostatic potential, headgroup angle, and water dipole orientation. To compare the diffusive nature of anesthetic through the lipid membrane quantitatively, we investigated the hexadecane/water partition coefficient using expanded ensemble simulation. We predicted the permeability coefficient of anesthetics in the following order: uncharged tetracaine > uncharged procaine > charged tetracaine > charged procaine. We also shown that the charged forms of drugs are more potent in hydrogen bonding, disturbing the lipid headgroups, changing the orientation of water dipoles, and increasing the headgroup electrostatic potential more than uncharged drugs, while the uncharged drugs make the lipid diffusion faster and increase the tail order parameter. The results of these simulation studies suggest that the different forms of anesthetics induce different structural modifications in the lipid bilayer, which provides new insights into their molecular mechanism.

Naratriptan aggregation in lipid bilayers: perspectives from molecular dynamics simulations

In order to understand the interaction between naratriptan and a fully hydrated bilayer of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidyl-choline (POPC), we carried out molecular dynamics simulations. The simulations were performed considering neutral and protonated ionization states, starting from different initial conditions. At physiological pH, the protonated state of naratriptan is predominant. It is expected that neutral compounds could have larger membrane partition than charged compounds. However, for the specific case of triptans, it is difficult to study neutral species in membranes experimentally, making computer simulations an interesting tool. When the naratriptan molecules were originally placed in water, they partitioned between the bilayer/water interface and water phase, as has been described for similar compounds. From this condition, the drugs displayed low access to the hydrophobic environment, with no significant effects on bilayer organization. The molecules anchored in the interface, due mainly to the barrier function of the polar and oriented lipid heads. On the other hand, when placed inside the bilayer, both neutral and protonated naratriptan showed self-aggregation in the lipid tail environment. In particular, the protonated species exhibited a pore-like structure, dragging water through this environment. Graphical Abstract Different behaviour of Naratriptan and Sumatriptan, when the drugs were originally placed in the lipid core.

Molecular Dynamics Simulation Study of Permeation of Molecules through Skin Lipid Bilayer

Stratum Corneum (SC), the outermost layer of skin, is mainly responsible for skin's barrier function. The complex lipid matrix of SC determines these barrier properties. In this study, the lipid matrix is modeled as an equimolar mixture of ceramide (CER), cholesterol (CHOL), and free fatty acid (FFA). The permeation of water, oxygen, ethanol, acetic acid, urea, butanol, benzene, dimethyl sulfoxide (DMSO), toluene, phenol, styrene, and ethylbenzene across this layer is studied using a constrained MD simulations technique. Several long constrained simulations are performed at a skin temperature of 310 K under NPT conditions. The free energy profiles and diffusion coefficients along the bilayer normal have been calculated for each molecule. Permeability coefficients are also calculated and compared with experimental data. The main resistance for the permeation of hydrophilic and hydrophobic permeants has been found to be in the interior of the lipid bilayer and near the lipid-water interface, respectively. The obtained permeability is found to be a few orders of magnitude higher than experimental values for hydrophilic molecules while for hydrophobic molecules more discrepancy was observed. Overall, the qualitative ranking is consistent with the experiments.

In silico pharmacology: Drug membrane partitioning and crossing

Over the past decade, molecular dynamics (MD) simulations have become particularly powerful to rationalize drug insertion and partitioning in lipid bilayers. MD simulations efficiently support experimental evidences, with a comprehensive understanding of molecular interactions driving insertion and crossing. Prediction of drug partitioning is discussed with respect to drug families (anesthetics; β-blockers; non-steroidal anti-inflammatory drugs; antioxidants; antiviral drugs; antimicrobial peptides). To accurately evaluate passive permeation coefficients turned out to be a complex theoretical challenge; however the recent methodological developments based on biased MD simulations are particularly promising. Particular attention is paid to membrane composition (e.g., presence of cholesterol), which influences drug partitioning and permeation. Recent studies concerning in silico models of membrane proteins involved in drug transport (influx and efflux) are also reported here. These studies have allowed gaining insight in drug efflux by, e.g., ABC transporters at an atomic resolution, explicitly accounting for the mandatory forces induced by the surrounded lipid bilayer. Large-scale conformational changes were thoroughly analyzed.

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.

Simulation-Based Approaches for Determining Membrane Permeability of Small Compounds

Predicting the rate of nonfacilitated permeation of solutes across lipid bilayers is important to drug design, toxicology, and signaling. These rates can be estimated using molecular dynamics simulations combined with the inhomogeneous solubility-diffusion model, which requires calculation of the potential of mean force and position-dependent diffusivity of the solute along the transmembrane axis. In this paper, we assess the efficiency and accuracy of several methods for the calculation of the permeability of a model DMPC bilayer to urea, benzoic acid, and codeine. We compare umbrella sampling, replica exchange umbrella sampling, adaptive biasing force, and multiple-walker adaptive biasing force for the calculation of the transmembrane PMF. No definitive advantage for any of these methods in their ability to predict the membrane permeability coefficient Pm was found, provided that a sufficiently long equilibration is performed. For diffusivities, a Bayesian inference method was compared to a generalized Langevin method, both being sensitive to chosen parameters and the slow relaxation of membrane defects. Agreement within 1.5 log units of the computed Pm with experiment is found for all permeants and methods. Remaining discrepancies can likely be attributed to limitations of the force field as well as slowly relaxing collective movements within the lipid environment. Numerical calculations based on model profiles show that Pm can be reliably estimated from only a few data points, leading to recommendations for calculating Pm from simulations.

Comparative exploration of hydrogen sulfide and water transmembrane free energy surfaces via orthogonal space tempering free energy sampling

Hydrogen sulfide (H2 S), a commonly known toxic gas compound, possesses unique chemical features that allow this small solute molecule to quickly diffuse through cell membranes. Taking advantage of the recent orthogonal space tempering (OST) method, we comparatively mapped the transmembrane free energy landscapes of H2 S and its structural analogue, water (H2 O), seeking to decipher the molecular determinants that govern their drastically different permeabilities. As revealed by our OST sampling results, in contrast to the highly polar water solute, hydrogen sulfide is evidently amphipathic, and thus inside membrane is favorably localized at the interfacial region, that is, the interface between the polar head-group and nonpolar acyl chain regions. Because the membrane binding affinity of H2 S is mainly governed by its small hydrophobic moiety and the barrier height inbetween the interfacial region and the membrane center is largely determined by its moderate polarity, the transmembrane free energy barriers to encounter by this toxic molecule are very small. Moreover when H2 S diffuses from the bulk solution to the membrane center, the above two effects nearly cancel each other, so as to lead to a negligible free energy difference. This study not only explains why H2 S can quickly pass through cell membranes but also provides a practical illustration on how to use the OST free energy sampling method to conveniently analyze complex molecular processes. © 2015 Wiley Periodicals, Inc.