Production of liposomes by microfluidics: The impact of post-manufacturing dilution on drug encapsulation and lipid loss

Microfluidic mixing is recognized as a convenient method to produce liposomes for its scalability and reproducibility. Numerous studies have described the effect of process parameters such as flow rate ratios and total flow rate on size and size distribution of vesicles. In this work, we focused our attention on the effect of flow rate ratios on the encapsulation efficiency of liposomes, as we hypothesized that different amount of residual organic solvent could affect the retention of lipophilic drug molecules within the bilayer. In a further step, we investigated how the liposomes integrity and loading were impacted by different methods of solvent removal: direct dialysis and dilution & dialysis. Liposomes were prepared by rapidly mixing an ethanolic solution of lipids and a model drug with buffer in a herringbone micromixer, employing four different flow rate ratios (FRR, 4:1, 7:3, 3:2, 1:1). Quercetin, resveratrol and ascorbyl palmitate were used as model antioxidant drugs with different lipophilicity. Data showed that liposomes produced using lower flow rate ratios (i.e., with more residual ethanol) had lower encapsulation efficiencies as well as a more prominent loss of lipids from the bilayer following purification with direct dialysis. If the amount of residual ethanol was reduced to 5% (dilution & dialysis method), the lipids and drug leakage was prevented. Such effect was correlated with the drug aggregation propensity in different ethanol/water mixtures measured by molecular dynamics simulations. Overall, these results highlight the need to tailor the purification method basing on the molecular properties of the loaded drug to ensure high encapsulation and limit the waste of material

a coarse-grained molecular dynamics description of docetaxel-conjugate release from plga matrices

Despite the extensive use of poly-lactic-glycolic-acid (PLGA) in biomedical applications, computational research on the mesoscopic characterization of PLGA-based delivery systems is limited. In this study, a computational model for PLGA is proposed, developed, and validated for the reproducibility of transport properties that can influence drug release, the rate of which remains difficult to control. For computational efficiency, coarse-grained models of the molecular components under consideration are built using the MARTINI force field version 2.2. The translocation free energy barrier 〖ΔG〗_t^* across the PLGA matrix in the aqueous phase of docetaxel and derivatives of varying size and solubility is predicted via Molecular Dynamics simulations and compared with the experimental release data. The thermodynamic quantity 〖ΔG〗_t^* anticipate, if predicted via simulations, and can help to explain the release kinetics of hydrophobic compounds from the PLGA matrix, albeit within the limit of a drug concentration below a critical aggregation concentration. In addition, this study provides a viable approach for optimizing the particle coating with ligands having desired orientation and stability

Peptide-induced membrane curvature in edge-stabilized open bilayers: A theoretical and molecular dynamics study

Peptide- or protein-induced curvatures of lipid membranes may be studied in molecular dynamics (MD) simulations. In these, membranes are usually modeled as infinitely extended bilayers by using periodic boundary conditions. However, the enforced periodicity results in an underestimation of the bending power of peptides, unless the patch size is much larger than the induced curvature radii. In this letter, we propose a novel approach to evaluate the bending power of a given distribution and/or density of peptides based on the use of flat open-edged lipid patches. To ensure long-lived metastable structures, the patch rim is stabilized in MD simulations by a local enrichment with short-chain lipids. By combining the theory of continuum elastic media with MD simulations, we prove that open-edged patches evolve from a planar state to a closed vesicle, with a transition rate that strongly depends on the concentration of lipid soluble peptides. For close-to-critical values for the patch size and edge energy, the response to even small changes in peptide concentration adopts a transition-like behavior (buckling instability). The usage of open-edged membrane patches amplifies the bending power of peptides, thereby enabling the analysis of the structural properties of membrane-peptide systems. We applied the presented method to investigate the curvature induced by aggregating β -amyloid peptides, unraveling a strong sensitivity of membrane deformation to the peptide concentration

The thermodynamics of simple biomembrane mimetic systems

Insight into the forces governing a system is essential for understanding its behavior and function. Thermodynamic investigations provide a wealth of information that is not, or is hardly, available from other methods. This article reviews thermodynamic approaches and assays to measure collective properties such as heat adsorption / emission and volume variations. These methods can be successfully applied to the study of lipid vesicles (liposomes) and biological membranes. With respect to instrumentation, differential scanning calorimetry, pressure perturbation calorimetry, isothermal titration calorimetry, dilatometry, and acoustic techniques aimed at measuring the isothermal and adiabatic processes, two- and three-dimensional compressibilities are considered. Applications of these techniques to lipid systems include the measurement of different thermodynamic parameters and a detailed characterization of thermotropic, barotropic, and lyotropic phase behavior. The membrane binding and / or partitioning of solutes (proteins, peptides, drugs, surfactants, ions, etc.) can also be quantified and modeled. Many thermodynamic assays are available for studying the effect of proteins and other additives on membranes, characterizing non-ideal mixing, domain formation, bilayer stability, curvature strain, permeability, solubilization, and fusion. Studies of membrane proteins in lipid environments elucidate lipid–protein interactions in membranes. Finally, a plethora of relaxation phenomena toward equilibrium thermodynamic structures can be also investigated. The systems are described in terms of enthalpic and entropic forces, equilibrium constants, heat capacities, partial volume changes, volume and area compressibility, and so on, also shedding light on the stability of the structures and the molecular origin and mechanism of the structural changes

Anomalous viscosity effect in the early stages of the ion-assisted adhesion/fusion event between lipid bilayers: A theoretical and computational study

The effect of viscosity on the encounter rate of two interacting membranes was investigated by combining a non-equilibrium Fokker-Planck model together with extensive Molecular Dynamics (MD) calculations. The encounter probability and stabilization of transient contact points represent the preliminary steps toward short-range adhesion and fusion of lipid leaflets. To strengthen our analytical model, we used a Coarse Grained MD method to follow the behavior of two charged palmitoyl oleoyl phosphatidylglycerol membranes embedded in a electrolyte-containing box at different viscosity regimes. Solvent friction was modulated by varying the concentration of a neutral, water-soluble polymer, polyethylene glycol, while contact points were stabilized by divalent ions that form bridges among juxtaposed membranes. While a naïve picture foresees a monotonous decrease of the membranes encounter rate with solvent viscosity, both the analytical model and MD simulations show a complex behavior. Under particular conditions, the encounter rate could exhibit a maximum at a critical viscosity value or for a critical concentration of bridging ions. These results seem to be confirmed by experimental observations taken from the literature

Anomalous viscosity effect in the early stages of the ion-assisted adhesion/fusion event between lipid bilayers:A theoretical and computational study

<p>The effect of viscosity on the encounter rate of two interacting membranes was investigated by combining a non-equilibrium Fokker-Planck model together with extensive Molecular Dynamics (MD) calculations. The encounter probability and stabilization of transient contact points represent the preliminary steps toward short-range adhesion and fusion of lipid leaflets. To strengthen our analytical model, we used a Coarse Grained MD method to follow the behavior of two charged palmitoyl oleoyl phosphatidylglycerol membranes embedded in a electrolyte-containing box at different viscosity regimes. Solvent friction was modulated by varying the concentration of a neutral, water-soluble polymer, polyethylene glycol, while contact points were stabilized by divalent ions that form bridges among juxtaposed membranes. While a naive picture foresees a monotonous decrease of the membranes encounter rate with solvent viscosity, both the analytical model and MD simulations show a complex behavior. Under particular conditions, the encounter rate could exhibit a maximum at a critical viscosity value or for a critical concentration of bridging ions. These results seem to be confirmed by experimental observations taken from the literature. (C) 2013 AIP Publishing LLC.</p>

Zoledronate Derivatives as Potential Inhibitors of Uridine Diphosphate-Galactose Ceramide Galactosyltransferase 8:A Combined Molecular Docking and Dynamic Study

Krabbe's disease is a neurodegenerative disorder caused by deficiency of galactocerebrosidase activity that affects the myelin sheath of the nervous system, involving dysfunctional metabolism of sphingolipids. It has no cure. Because substrate inhibition therapy has been shown to be effective in some human lysosomal storage diseases, we hypothesize that a substrate inhibition therapeutic approach might be appropriate to allow correction of the imbalance between formation and breakdown of glycosphingolipids and to prevent pathological storage of psychosine. The enzyme responsible for the biosynthesis of galactosylceramide and psychosine is uridine diphosphate-galactose ceramide galactosyltransferase (2-hydroxyacylsphingosine 1-b-galactosyltransferase; UGT8; EC 2.4.1.45), which catalyzes the transferring of galactose from uridine diphosphate-galactose to ceramide or sphingosine, an important step of the biosynthesis of galactosphingolipids. Because some bisphosphonates have been identified as selective galactosyltransferase inhibitors, we verify the binding affinity to a generated model of the enzyme UGT8 and investigate the molecular mechanisms of UGT8-ligand interactions of the bisphosphonate zoledronate by a multistep framework combining homology modeling, molecular docking, and molecular dynamics simulations. From structural information on UGTs' active site stereochemistry, charge density, and access through the hydrophobic environment, the molecular docking procedure allowed us to identify zoledronate as a potential inhibitor of human ceramide galactosyltransferase. More importantly, zoledronate derivates were designed through computational modeling as putative new inhibitors. Experiments in vivo and in vitro have been planned to verify the possibility of using zoledronate and/or the newly identified inhibitors of UGT8 for a substrate inhibition therapy useful for treatment of Krabbe's disease and/or other lysosomal disorders

A unifying framework for amyloid-mediated membrane damage: The lipid-chaperon hypothesis

Over the past thirty years, researchers have highlighted the role played by a class of proteins or polypeptides that forms pathogenic amyloid aggregates in vivo, including i) the amyloid Abeta peptide, which is known to form senile plaques in Alzheimer's disease; ii) alpha-synuclein, responsible for Lewy body formation in Parkinson's disease and iii) IAPP, which is the protein component of type 2 diabetes-associated islet amyloids. These proteins, known as intrinsically disordered proteins (IDPs), are present as highly dynamic conformational ensembles. IDPs can partially (mis) fold into (dys) functional conformations and accumulate as amyloid aggregates upon interaction with other cytosolic partners such as proteins or lipid membranes. In addition, an increasing number of reports link the toxicity of amyloid proteins to their harmful effects on membrane integrity. Still, the molecular mechanism underlying the amyloidogenic proteins transfer from the aqueous environment to the hydrocarbon core of the membrane is poorly understood. This review starts with a historical overview of the toxicity models of amyloidogenic proteins to contextualize the more recent lipid-chaperone hypothesis. Then, we report the early molecular-level events in the aggregation and ion-channel pore formation of Abeta, IAPP, and alpha-synuclein interacting with model membranes, emphasizing the complexity of these processes due to their different spatial-temporal resolutions. Next, we underline the need for a combined experimental and computational approach, focusing on the strengths and weaknesses of the most commonly used techniques. Finally, the last two chapters highlight the crucial role of lipid-protein complexes as molecular switches among ion-channel-like formation, detergent-like, and fibril formation mechanisms and their implication in fighting amyloidogenic diseases.Comment: 45 pages, 3 figure

Out of Equilibrium Divergence of Dissipation in an Oscillating Bubble Coated by Surfactants

We report measurements of the relaxation and resonance frequency of forced oscillating bubbles covered by a layer of surface-active molecules, the anionic surfactant sodium dodecyl sulfate (SDS). Less systematic investigations have been also carried out on neutral and cationic surfactants. A divergence of the viscous damping is observed at a very low bulk concentration. Subtle variations in the resonance peak are also measured. Bubble oscillations are driven by an electric field and measured with a sensitive interferometric technique. Results are interpreted with a model which takes care of the coupling between the dynamics of fluid surface oscillations and the properties of a surfactant monolayer in the vicinity of the phase transition from a gas-like distribution to a liquid-like assembly (the so-called gas–LE transition). Important charge effects are also considered. The basic assumptions of the model (cooperative adsorption of the surfactant at the air–water interface and coupling between the shape of the deformed surface and the local surfactant concentration) have been fully confirmed by extensive coarse-grained molecular dynamics simulations on model systems