Martini-3 Coarse-Grained Models for the Bacterial Lipopolysaccharide Outer Membrane of Escherichia coli

The outer lipopolysaccharide (LPS) membrane of Gram-negative bacteria forms the main barrier for transport of antimicrobial molecules into the bacterial cell. In this study we develop coarse-grained models for the outer membrane of Escherichia coli in the Martini-3 framework. The coarse-grained model force field was parametrized and validated using all-atom simulations of symmetric membranes of lipid A and rough LPS as well as a complete asymmetric membrane of LPS with the O-antigen. The bonded parameters were obtained using an iterative refinement procedure with target bonded distributions obtained from all-atom simulations. The membrane thickness, area of the LPS, and density distributions for the different regions as well as the water and ion densities in Martini-3 simulations show excellent agreement with the all-atom data. Additionally the solvent accessible surface area for individual molecules in water was found to be in good agreement. The binding of calcium ions with phosphate and carboxylate moieties of LPS is accurately captured in the Martini-3 model, indicative of the integrity of the highly negatively charged LPS molecules in the outer membranes of Gram-negative bacteria. The melting transition of the coarse-grained lipid A membrane model was found to occur between 300 and 310 K, and the model captured variations in area per LPS, order parameter, and membrane thickness across the melting transition. Our study reveals that the proposed Martini-3 models for LPS are able to capture the physicochemical balance of the complex sugar architecture of the outer membrane of Escherichia coli. The coarse-grained models developed in this study would be useful for determining membrane protein interactions and permeation of potential antimicrobials through bacterial membranes at mesoscopic spatial and temporal scales.

Structure and Elasticity of Mitochondrial Membranes: A Molecular Dynamics Simulation Study

Mitochondria are known as the powerhouse of the cell because they produce energy in the form of adenosine triphosphate. They also have other crucial functions such as regulating apoptosis, calcium homeostasis, and reactive oxygen species production. To perform these diverse functions, mitochondria adopt specific structures and frequently undergo dynamic shape changes, indicating that their mechanical properties play an essential role in their functions. To gain a detailed understanding at the molecular level of the structure and mechanical properties of mitochondria, we carry out atomistic molecular dynamics simulations for three inner mitochondrial membranes and three outer mitochondrial membrane models. These models take into account variations in cardiolipin and cholesterol concentrations as well as the symmetry/asymmetry between the two leaflets. Our simulations allow us to calculate various structural quantities and the bending, twisting, and tilting elastic moduli of the membrane models. Our results indicate that the structures of the inner and outer mitochondrial membranes are quite similar and do not depend much on the variation in lipid compositions. However, the bending modulus of the membranes increases with increasing concentrations of cardiolipin or cholesterol but decreases with a membrane asymmetry. Notably, we found that the dipole potential of the membrane increases with an increasing cardiolipin concentration. Finally, possible roles of cardiolipin in regulating ion and proton currents and maintaining the cristate are discussed in some details.

RE, Mocskos EE, Facelli JC, Pickholz M, Ferraro MB. Molecular dynamics study of the mechanical properties of drug loaded model systems: A comparison of a polymersome with a bilayer

In this work we implement a new methodology to study structural and mechanical properties of systems having spherical and planar symmetries throughout Molecular Dynamics simulations. This methodology is applied here to a drug delivery system based in polymersomes, as an example. The chosen model drug was the local anesthetic prilocaine due to previous parameterization within the used coarse grain scheme. In our approach, mass density profiles (MDPs) are used to obtain key structural parameters of the systems, and pressure profiles are used to estimate the curvature elastic parameters. The calculation of pressure profiles and radial MPDs required the development of specific methods, which were implemented in an in-house built version of the GROMACS 2018 code. The methodology presented in this work is applied to characterize poly(ethylene oxide)-poly(butadiene) polymersomes and bilayers loaded with the model drug prilocaine. Our results show that structural properties of the polymersome membrane could be obtained from bilayer simulations, with significantly lower computational cost compared to whole polymersome simulations, but the bilayer simulations are insufficient to get insights on their mechanical aspects, since the elastic parameters are canceled out for the complete bilayer (as consequence of the symmetry). The simulations of entire polymersomes, although more complex, offer a complementary approach to get insights on the mechanical behavior of the systems.

Coarse-Graining Waters: Unveiling The Effective Hydrophilicity/Hydrophobicity of Individual Protein Atoms and The Roles of Waters' Hydrogens

There have been many coarse-graining methods developed that aim to reduce the sizes of simulated systems and their computational costs. In this work, we develop a new coarse-graining method, called coarse-graining-delta (or δ-CG in short), that reduces the degrees of freedom of the potential energy surface by coarse-graining relative locations of atoms from their unit centers. Our method extends and generalizes the methods used in the coarse-grained normal mode analysis and enables us to study the roles of the individual removed atoms in a system, which have been difficult to study in molecular dynamics simulations. By applying δ-CG to coarse-grain three-point water molecules into single-point solvent particles, we successfully identify the effective hydrophilicity and hydrophobicity of all the individual protein atom types, which collectively correlate well with the known hydrophilic, hydrophobic, and amphipathic characteristics of amino acids. Moreover, our investigation shows that water's hydrogens have two roles in interacting with protein atoms. First, water molecules adjust their poses around different amino acids and their atoms, and the statistical preferences of the hydrogen poses near the atoms determine the effective hydrophilicity and hydrophobicity of the atoms, which have not been successfully addressed before. Second, the collective dynamics of the hydrogens assist the water molecules in escaping from the potential energy wells of the hydrophilic atoms. Our method also shows that coarse-graining a system mathematically leads to breaking antisymmetry of the nonbonded interactions; as a result, two interacting coarse-grained units exert different forces on each other. Our study indicates that the accuracy of coarse-grained force fields, such as the MARTINI force field and the UNRES force field, can be improved in two ways: (i) refining their potential energy functions and coefficients by analyzing the coarse-grained potential energy surface using δ-CG, and (ii) introducing non-antisymmetric interactions.

Recent Advances of Studies on Cell-Penetrating Peptides Based on Molecular Dynamics Simulations

With the ability to transport cargo molecules across cell membranes with low toxicity, cell-penetrating peptides (CPPs) have become promising candidates for next generation peptide-based drug delivery vectors. Over the past three decades since the first CPP was discovered, a great deal of work has been done on the cellular uptake mechanisms and the applications for the delivery of therapeutic molecules, and significant advances have been made. But so far, we still do not have a precise and unified understanding of the structure-activity relationship of the CPPs. Molecular dynamics (MD) simulations provide a method to reveal peptide-membrane interactions at the atomistic level and have become an effective complement to experiments. In this paper, we review the progress of the MD simulations on CPP-membrane interactions, including the computational methods and technical improvements in the MD simulations, the research achievements in the CPP internalization mechanism, CPP decoration and coupling, and the peptide-induced membrane reactions during the penetration process, as well as the comparison of simulated and experimental results.

Membrane Lipid Reshaping Underlies Oxidative Stress Sensing by the Mitochondrial Proteins UCP1 and ANT1

Oxidative stress and ROS are important players in the pathogenesis of numerous diseases. In addition to directly altering proteins, ROS also affects lipids with negative intrinsic curvature such as phosphatidylethanolamine (PE), producing PE adducts and lysolipids. The formation of PE adducts potentiates the protonophoric activity of mitochondrial uncoupling proteins, but the molecular mechanism remains unclear. Here, we linked the ROS-mediated change in lipid shape to the mechanical properties of the membrane and the function of uncoupling protein 1 (UCP1) and adenine nucleotide translocase 1 (ANT1). We show that the increase in the protonophoric activity of both proteins occurs due to the decrease in bending modulus in lipid bilayers in the presence of lysophosphatidylcholines (OPC and MPC) and PE adducts. Moreover, MD simulations showed that modified PEs and lysolipids change the lateral pressure profile of the membrane in the same direction and by the similar amplitude, indicating that modified PEs act as lipids with positive intrinsic curvature. Both results indicate that oxidative stress decreases stored curvature elastic stress (SCES) in the lipid bilayer membrane. We demonstrated that UCP1 and ANT1 sense SCES and proposed a novel regulatory mechanism for the function of these proteins. The new findings should draw the attention of the scientific community to this important and unexplored area of redox biochemistry.

In Cellulo and In Vivo Comparison of Cholesterol, Beta-Sitosterol and Dioleylphosphatidylethanolamine for Lipid Nanoparticle Formulation of mRNA

Lipid Nanoparticles (LNPs) are a leading class of mRNA delivery systems. LNPs are made of an ionizable lipid, a polyethyleneglycol (PEG)-lipid conjugate and helper lipids. The success of LNPs is due to proprietary ionizable lipids and appropriate helper lipids. Using a benchmark lipid (D-Lin-MC3) we compared the ability of three helper lipids to transfect dendritic cells in cellulo and in vivo. Studies revealed that the choice of helper lipid does not influence the transfection efficiency of immortalized cells but, LNPs prepared with DOPE (dioleylphosphatidylethanolamine) and β-sitosterol were more efficient for mRNA transfection in murine dendritic cells than LNPs containing DSPC (distearoylphosphatidylcholine). This higher potency of DOPE and β-sitosterol LNPs for mRNA expression was also evident in vivo but only at low mRNA doses. Overall, these data provide valuable insight for the design of novel mRNA LNP vaccines.

Design and simulation of the liposomal model by using a coarse-grained molecular dynamics approach towards drug delivery goals

The simulated liposome models provide events in molecular biological science and cellular biology. These models may help to understand the cell membrane mechanisms, biological cell interactions, and drug delivery systems. In addition, the liposomes model may resolve specific issues such as membrane transports, ion channels, drug penetration in the membrane, vesicle formation, membrane fusion, and membrane protein function mechanism. One of the approaches to investigate the lipid membranes and the mechanism of their formation is by molecular dynamics (MD) simulations. In this study, we used the coarse-grained MD simulation approach and designed a liposome model system. To simulate the liposome model, we used phospholipids that are present in the structure of natural cell membranes (1,2-Dioleoyl-sn-glycero-3-phosphocholine (DOPC) and 1,2-Dioleoyl-sn-glycero-3-phosphoethanolamine (DOPE)). Simulation conditions such as temperature, ions, water, lipid concentration were performed based on experimental conditions. Our results showed a liposome model (ellipse vesicle structure) during the 2100 ns was formed. Moreover, the analysis confirmed that the stretched and ellipse structure is the best structure that could be formed. The eukaryotic and even the bacterial cells have elliptical and flexible structures. Usually, an elliptical structure is more stable than other assembled structures. The results indicated the assembly of the lipids is directed through short-range interactions (electrostatic interactions and, van der Waals interactions). Total energy (Van der Waals and electrostatic interaction energy) confirmed the designed elliptical liposome structure has suitable stability at the end of the simulation process. Our findings confirmed that phospholipids DOPC and DOPE have a good tendency to form bilayer membranes (liposomal structure) based on their geometric shapes and chemical-physical properties. Finally, we expected the simulated liposomal structure as a simple model to be useful in understanding the function and structure of biological cell membranes. Furthermore, it is useful to design optimal, suitable, and biocompatible liposomes as potential drug carriers.

Chloroplast membrane lipid remodeling protects against dehydration by limiting membrane fusion and distortion

Dehydration damages the structural integrity of the chloroplast membrane and, consequently, the normal photosynthetic function of this organelle. Remodeling of galactolipids by converting monogalactosyl-diacylglycerol (MGDG) to digalactosyl-diacylglycerol (DGDG) and oligo-galactolipids is an effective adaptation strategy for protecting against dehydration damage to the chloroplast membrane. However, detailed molecular mechanisms are missing. In this study, by performing molecular-level simulations of bi-lamellar membranes under various dehydration conditions, we find that MGDG-to-DGDG remodeling protects the chloroplast membrane in a unique manner by simultaneously dictating both the extent and the pattern of fusion stalks formed with the apposed membrane. Specifically, MGDG-rich membranes form elongated stalks at a moderate dehydration level, whereas DGDG-rich membranes form smaller, rounded stalks. Simulations of wild-type and mutant Arabidopsis (Arabidopsis thaliana) outer chloroplast membranes further confirm that the mutant membrane without galactolipid remodeling is more susceptible to membrane fusion due to its higher MGDG content. Our work reveals the underlying physical mechanisms that govern the pattern and extent of membrane fusion structures, paving the way for rational genetic engineering of crops with improved dehydration tolerance.

Mechanistic Insight into How PEGylation Reduces the Efficacy of pH-Sensitive Liposomes from Molecular Dynamics Simulations

Liposome-based drug delivery systems composed of DOPE stabilized with cholesteryl hemisuccinate (CHMS) have been proposed as a drug delivery mechanism with pH-triggered release as the anionic form (CHSa) is protonated (CHS) at reduced pH; PEGylation is known to decrease this pH sensitivity. In this manuscript, we set out to use molecular dynamics (MD) simulations with a model with all-atom resolution to provide insight into why incorporation of poly(ethyleneglycol) (PEG) into DOPE–CHMS liposomes reduces their pH sensitivity; we also address two additional questions: (1) How CHSa stabilizes DOPE bilayers into a lamellar conformation at a physiological pH of 7.4? and (2) how the change from CHSa to CHS at acidic pH triggers the destabilization of DOPE bilayers? We found that (A) CHSa stabilizes the DOPE lipid membrane by increasing the hydrophilicity of the bilayer surface, (B) when CHSa changes to CHS by pH reduction, DOPE bilayers are destabilized due to a reduction in bilayer hydrophilicity and a reduction in the area per lipid, and (C) PEG stabilizes DOPE bilayers into the lamellar phase, thus reducing the pH sensitivity of the liposomes by increasing the area per lipid through penetration into the bilayer, which is our main focus.

Membrane shape as determinant of protein properties

Membrane lipids play a role in the modulation of a variety of biological processes. This is often achieved through fine-tuned changes in membrane physical and chemical properties. While some membrane physical properties (e.g., curvature, lipid domains, fluidity) have received increased scientific attention over the years, only recently has membrane shape emerged as an active modulator of protein properties. Biological membranes are mostly found organized into a lipid bilayer arrangement, in which the spontaneous shape is an intrinsically flat, planar morphology (in relation to the size of proteins). However, it is known that many cells and organelles have non-planar morphologies. In addition, perturbations in membrane morphology occur in a variety of biological processes. Recent studies have shown that membrane shape can modulate a variety of biological processes by determining protein properties. While membrane shape generation modulates proteins via changes in membrane mechanical properties, membrane shape recognition regulates proteins by providing the optimal surface for interaction. Hence, membranes have evolved an elegant mechanism to couple mesoscopic perturbations to molecular properties and vice-versa. In this review, the regulation of the enzymatic properties of two isoforms of mammalian diacylglycerol kinase, which play important roles in cellular signal transductions, will be used to exemplify the recent advancements in the field of membrane shape recognition, as well as future challenges and perspectives.

Comparison of umbrella sampling and steered molecular dynamics methods for computing free energy profiles of aromatic substrates through phospholipid bilayers

Understanding the permeation of molecules through lipid membranes is fundamental for predicting the cellular uptake of solutes and drug delivery mechanisms. In molecular simulations, the usual approach is to compute the free energy (FE) profile of a molecule across a model lipid bilayer, which can then be used to estimate the permeability of the molecule. Umbrella Sampling (US), which involves carrying out a series of biased simulations along a defined reaction coordinate (usually the bilayer normal direction), is a popular method for the computation of such FE profiles. However, US can be challenging to implement because the results are dependent on the strength of the biasing potential and the spacing of windows along the reaction coordinate, which, in practice, are usually optimized by an inefficient trial and error approach. The Steered Molecular Dynamics implementation of the Jarzynski Equality (JE-SMD) has been identified as an alternative to equilibrium sampling methods for measuring the FE change across a reaction coordinate. In the JE-SMD approach, equilibrium FE values are evaluated from the average of rapid non-equilibrium trajectories, thus avoiding the practical issues that come with US. Here, we use three different corrections of the JE-SMD method to calculate the FE change for the translocation of two aromatic substrates, phenylalanine and toluene, across a lipid bilayer and compare the accuracy and computational efficiency of these approaches to the results obtained using US. We show evidence that when computing the free energy profile, the JE-SMD approach suffers from insufficient sampling convergence of the bilayer environment and is dependent on the characteristic of the aromatic substrate itself. We deduce that, despite its drawbacks, US remains the more viable approach of the two for computing the FE profile.

Lipid-Modified Peptide Nucleic Acids: Synthesis and Application to Programmable Liposome Fusion

Peptide nucleic acids (PNAs) can be modified with aliphatic lipid chains and designed to be water soluble and able to spontaneously insert into phospholipid bilayers. Liposomes with 1.5% negatively charged POPG can be driven to fuse and mix their inner content volumes via functionalization with such lipidated peptide nucleic acids (LiPNAs). During fusion, only low amounts of leakage occur (<5%). We describe here the synthesis and purification of such LiPNAs using an automated peptide synthesizer and the preparation of LiPNA functionalized liposomes. Further, we describe the measurement of LiPNA-induced fusion using a fluorescence-based assay for the content mixing between a liposome population with an encapsulated self-quenching fluorescent dye (SRB) and a buffer-filled liposome population.

Formation and Properties of a Self-Assembled Nanoparticle-Supported Lipid Bilayer Probed through Molecular Dynamics Simulations

We have carried out coarse-grained molecular dynamics (MD) simulations to study the self-assembly procedure of a system of randomly placed lipid molecules, water beads, and a nanoparticle (NP). The self-assembly results in the formation of the nanoparticle-supported lipid bilayer (NPSLBL), with the self-assembly mechanism being driven by events such as the formation of small lipid clusters, merging of the lipid clusters in the vicinity of the NP to form NP-embedded vesicle with a pore, and collapsing of that pore to eventually form the equilibrated NPSLBL system overcoming a large free-energy barrier. Subsequently, we quantify the properties and the configurations of this NPSLBL system. We reveal that unlike our proposition of an equal number of lipid molecules occupying the inner and outer leaflets in a recent report studying the properties of a preassembled lipid bilayer, the equilibrated self-assembled NPSLBL system demonstrates a much larger number of lipid molecules occupying the outer leaflet as compared to the inner leaflet. Second, the thickness of the water layer entrapped between the NP and the inner leaflet shows similar values as predicted by experiments and our previous study. Finally, we reveal that, similar to our previous study, the diffusivity of the lipid molecules in the outer leaflet is larger than that in the inner leaflet but, due to higher temperature employed during our simulations, are even larger than that predicted by our previous study.

Self-Assembly of Phosphocholine Derivatives Using the ELBA Coarse-Grained Model: Micelles, Bicelles, and Reverse Micelles

The ELBA coarse-grained force field was originally developed for lipids, and its water model is described as a single-site Lennard-Jones particle with electrostatics modeled by an embedded point-dipole, while other molecules in this force field have a three (or four)-to-one mapping scheme. Here, ELBA was applied to investigate the self-assembly processes of dodecyl-phosphocholine (DPC) micelle, 1,2-dipalmitoyl-sn-glycero-3-phosphocholine/1,2-dihexaoyl-sn-glycero-3-phosphocholine (DPPC/DHPC) bicelles, and DPPC/cyclohexane/water reverse micelles through coarse-grained molecular dynamics (MD) simulations. New parameters were obtained using a simplex algorithm-based calibration procedure to determine the Lennard-Jones parameters for cyclohexane, dodecane, and cyclohexane-dodecane cross-interactions. Density, self-diffusion coefficient, surface tension, and mixture excess volume were found to be in fair agreement with experimental data. These new parameters were used in the simulations, and the obtained structures were analyzed for shape, size, volume, and surface area. Except for the shape of DPC micelles, all other properties match well with available experimental data and all-atom simulations. Remarkably, in agreement with experiments the rodlike shape of the DPPC reverse micelle is well described by ELBA, while all-atom data in the literature predicts a disclike shape. To further check the consistency of the force field in reproducing the correct shapes of reverse micelles, additional simulations were performed doubling the system size. Two distinct reverse micelles were obtained both presenting the rodlike shape and correct aggregation number.

A proteoliposome-based system reveals how lipids control photosynthetic light harvesting

Integral membrane proteins are exposed to a complex and dynamic lipid environment modulated by nonbilayer lipids that can influence protein functions by lipid-protein interactions. The nonbilayer lipid monogalactosyldiacylglycerol (MGDG) is the most abundant lipid in plant photosynthetic thylakoid membranes, but its impact on the functionality of energy-converting membrane protein complexes is unknown. Here, we optimized a detergent-based reconstitution protocol to develop a proteoliposome technique that incorporates the major light-harvesting complex II (LHCII) into compositionally well-defined large unilamellar lipid bilayer vesicles to study the impact of MGDG on light harvesting by LHCII. Using steady-state fluorescence spectroscopy, CD spectroscopy, and time-correlated single-photon counting, we found that both chlorophyll fluorescence quantum yields and fluorescence lifetimes clearly indicate that the presence of MGDG in lipid bilayers switches LHCII from a light-harvesting to a more energy-quenching mode that dissipates harvested light into heat. It is hypothesized that in the in vitro system developed here, MGDG controls light harvesting of LHCII by modulating the hydrostatic lateral membrane pressure profile in the lipid bilayer sensed by LHCII-bound peripheral pigments.

DNA-Mediated Liposome Fusion Observed by Fluorescence Spectrometry

DNA-programmed and controlled fusion of lipid membranes have recently been optimized to reliably mix the contents between two populations of liposomes, each functionalized with complementary lipidated DNA (LiNA) oligomer. In this chapter we describe a procedure for DNA-controlled fusion of liposomes mediated by LiNAs that are designed to force bilayers into close proximity. Using a self-quenching fluorescent dye (Sulforhodamine B) to monitor both the mixing of the internal volumes and leakage of the dye into the outer volume we measure the efficiency of content mixing in the bulk population, allowing for direct comparison between different LiNA designs. By generating samples for calibration corresponding to different amounts of content mixing, the average number of fusion events per labeled liposome can be estimated.

Reducing uncertainty in simulation estimates of the surface tension through a two-scale finite-size analysis: thicker is better

Recent simulation studies of the surface tension γ, and other properties of thin free-standing films, have revealed unexpected finite size effects in which the variance of the properties vary monotonically with the in-plane width of the films, complicating the extrapolation of estimates of film properties to the thermodynamic limit. We carried out molecular dynamics simulations to determine the origin of this phenomenon, and to address the practical problem of developing a more reliable methodology for estimating γ in the thermodynamic limit. We find that there are two distinct finite size effects that must be addressed in a finite size analysis of γ in thin films. The first finite size scale is the in-plane width of the films and the second scale is the simulation cell size in the transverse direction. Increasing the first scale enhances fluctuations in γ, measured by the standard deviation of their distribution, while increasing the second reduces γ fluctuations due to a corresponding increased 'freedom' of the film to fluctuate out of plane. We find that using progressively large simulation cells in the transverse direction, while keeping the film width fixed to an extent in which the full bulk liquid zone is developed, allows us to obtain a smooth extrapolation to the thermodynamic limit, enabling a reduction of the γ uncertainty to a magnitude on the order of 1% for systems having a reasonably large size, i.e., O (1 μm).

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.

Implicit solvent systematic coarse-graining of dioleoylphosphatidylethanolamine lipids: From the inverted hexagonal to the bilayer structure

Lamellar and hexagonal lipid structures are of particular importance in the biological processes such as membrane fusion and budding. Atomistic simulations of formation of these phases and transitions between them are computationally prohibitive, hence development of coarse-grained models is an important part of the methodological development in this area. Here we apply systematic bottom-up coarse-graining to model different phase structures formed by 1,2-dioleoylphosphatidylethanolamine (DOPE) lipid molecules. We started from atomistic simulations of DOPE lipids in water carried out at two different water/lipid molar ratio corresponding to the lamellar L_α and inverted hexagonal H_II structures at low and high lipid concentrations respectively. The atomistic trajectories were mapped to coarse-grained trajectories, in which each lipid was represented by 14 coarse-grained sites. Then the inverse Monte Carlo method was used to compute the effective coarse-grained potentials which for the coarse-grain model reproduce the same structural properties as the atomistic simulations. The potentials derived from the low concentration atomistic simulation were only able to form a bilayer structure, while both L_α and H_II lipid phases were formed in simulations with potentials obtained at high concentration. The typical atomistic configurations of lipids at high concentration combine fragments of both lamellar and non-lamellar structures, that is reflected in the extracted coarse-grained potentials which become transferable and can form a wide range of structures including the inverted hexagonal, bilayer, tubule, vesicle and micellar structures.

Physical properties of model biological lipid bilayers: insights from all-atom molecular dynamics simulations

The physical properties of lipid bilayers are sensitive to the specific type and composition of the lipids that make up the many different types of cell membranes. Studying model bilayers of representative heterogeneous compositions can provide key insights into membrane functionality. In this work, we use atomistic molecular dynamics simulations to characterize key properties in a number of bilayer membranes of varying composition. We first examine basic properties, such as lipid area, volume, and bilayer thickness, of simple, homogeneous bilayers comprising several lipid types, which are prevalent in biological membranes. Such lipids are then used in simulations of heterogeneous systems representative of bacterial, mammalian, and cancer membranes. Our analysis is especially focused on depth-dependent, transmembrane profiles; in particular, we calculate lateral pressure and dipole potential profiles, two fundamental properties which play key roles in a large number of biological functions.

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.

Mechanical properties of drug loaded diblock copolymer bilayers: A molecular dynamics study

In this work, we present results of coarse-grained simulations to study the encapsulation of prilocaine (PLC), both neutral and protonated, on copolymer bilayers through molecular dynamics simulations. Using a previously validated membrane model, we have simulated loaded bilayers at different drug concentrations and at low (protonated PLC) and high (neutral PLC) pH levels. We have characterized key structural parameters of the loaded bilayers in order to understand the effects of encapsulation of PLC on the bilayer structure and mechanical properties. Neutral PLC was encapsulated in the hydrophobic region leading to a thickness increase, while the protonated species partitioned between the water phase and the poly(ethylene oxide)-poly(butadiene) (PBD) interface, relaxing the PBD region and leading to a decrease in the thickness. The tangential pressures of the studied systems were calculated, and their components were decomposed in order to gain insights on their compensation. In all cases, it is observed that the loading of the membrane does not significantly decrease the stability of the bilayer, indicating that the system could be used for drug delivery.

Cell penetrating peptide modulation of membrane biomechanics by Molecular dynamics

The efficacy of a pharmaceutical treatment is often countered by the inadequate membrane permeability, that prevents drugs from reaching their specific intracellular targets. Cell penetrating peptides (CPPs) are able to route across cells' membrane various types of cargo, including drugs and nanoparticles. However, CPPs internalization mechanisms are not yet fully understood and depend on a wide variety of aspects. In this contest, the entry of a CPP into the lipid bilayer might induce molecular conformational changes, including marked variations on membrane's mechanical properties. Understanding how the CPP does influence the mechanical properties of cells membrane is crucial to design, engineer and improve new and existing penetrating peptides. Here, all atom Molecular Dynamics (MD) simulations were used to investigate the interaction between different types of CPPs embedded in a lipid bilayer of dioleoyl phosphatidylcholine (DOPC). In a greater detail, we systematically highlighted how CPP properties are responsible for modulating the membrane bending modulus. Our findings highlighted the CPP hydropathy strongly correlated with penetration of water molecules in the lipid bilayer, thus supporting the hypothesis that the amount of water each CPP can route inside the membrane is modulated by the hydrophobic and hydrophilic character of the peptide. Water penetration promoted by CPPs leads to a local decrease of the lipid order, which emerges macroscopically as a reduction of the membrane bending modulus.

Interfacial properties of the ionic liquid [bmim][triflate] over a wide range of temperatures

We carried out molecular dynamics simulations of the liquid/vacuum equilibrium of the ionic liquid [bmim][triflate] in a wide range of temperatures (323.15 to 573.15 K). The results showed liquid phases with high densities even at temperatures close to the decomposition temperature of the liquid. The density and surface tension behaviors are linear across this wide range of temperatures, which is an extension of the behaviors of these systems at low temperatures, where these properties have been experimentally measured. The interfacial region shows peaks of adsorption of the ions; they are ordered, with the alkyl chains of the [bmim] cations pointing out of the liquid, and the tailing angle of the chains becomes 90° at higher temperatures. The alkyl chains are part of the outermost interfacial region, where intra- and intermolecular tangential forces are in equilibrium; thus, they do not contribute to the total surface tension. Unlike simpler organic liquids, the surface tension is composed of positive normal contributions of intermolecular interactions; these are almost in equilibrium with the negative normal contributions of intramolecular interactions, which are mainly vibrations of the distance and the angle of valence. The pressure profiles show that the molecules are in 'crushed' conformations internally in the bulk liquid and even more so in the normal direction at the interface. The total pressure profiles show values with physical meaning, where the tangential peaks show higher values than normal pressures and give rise to the surface tension. Short cutoff radii for the calculation of intermolecular forces (less than 16.5 Å) produce a system that is not mechanically stable in the region of the bulk liquid (confirmed by radial distribution function calculations); this produces a difference between the normal pressure and the average of the tangential pressures, which affects the calculation of the surface tension due to overestimation by up to 20% when using the global expression, which is extensively used for the calculation of surface tension. The use of a sufficiently long cutoff radius avoids these mechanical balance problems.

The Mechanism by Which Luteolin Disrupts the Cytoplasmic Membrane of Methicillin-Resistant Staphylococcus aureus

Methicillin-resistant Staphylococcus aureus (MRSA) is one of the most versatile human pathogens. Luteolin (LUT) has anti-MRSA activity by disrupting the MRSA cytoplasmic membrane. However, the mechanism by which luteolin disrupts the membrane remains unclear. Here, we performed differential scanning calorimetry (DSC) and all-atomic molecular dynamics (AA-MD) simulations to investigate the interactions and effects of LUT on model membranes composed of phosphatidylcholine (PC) and phosphatidylglycerol (PG). We detected the transition thermodynamic parameters of dipalmitoylphosphatidylcholine (DPPC) liposomes, dipalmitoylphosphatidylglycerol (DPPG) liposomes, and liposomes composed of both DPPC and DPPG at different LUT concentrations and showed that LUT molecules were located between polar heads and the hydrophobic region of DPPC/DPPG. In the MD trajectories, LUT molecules ranging from 5 to 50 had different effects on the membranes thickness, fluidity and ordered property of lipids, and lateral pressure of lipid bilayers composed of dioleoylphosphatidylcholine (DOPC) and dioleoylphosphatidylglycerol (DOPG). Also, most LUT molecules were distributed in the region between the phosphorus atoms and C9 atoms of DOPC and DOPG. On the basis of the combination of these results, we conclude that the distinct effects of LUT on lipid bilayers composed of PCs and PGs may elucidate the mechanism by which LUT disrupts the cytoplasmic membrane of MRSA.

Effects of High Pressure on Phospholipid Bilayers

The response of lipid membranes to changes in external pressure is important for many biological processes, and it can also be exploited for technological applications. In this work, we employ all-atom molecular dynamics simulations to characterize the changes in the physical properties of phospholipid bilayers brought about by high pressure (1000 bar). In particular, we study how the response differs, in relation to different chain unsaturation levels, by comparing monounsaturated 1-palmitoyl-2-oleoyl-phosphatidylcholine (POPC) and biunsaturated dioleoyl-phosphatidylcholine (DOPC) bilayers. Various structural, mechanical, and dynamical features are found to be altered by the pressure increase in both bilayers. Notably, for most properties, including bilayer area and thickness, lipid order parameters, lateral pressure profile, and curvature frustration energy, we observe significantly more pronounced effects for monounsaturated POPC than biunsaturated DOPC. Possible biological implications of the results obtained are discussed, especially in relation to how different lipids can control the structure and function of membrane proteins.

Membrane Curvature and Lipid Composition Synergize To Regulate N-Ras Anchor Recruitment

Proteins anchored to membranes through covalently linked fatty acids and/or isoprenoid groups play crucial roles in all forms of life. Sorting and trafficking of lipidated proteins has traditionally been discussed in the context of partitioning to membrane domains of different lipid composition. We recently showed that membrane shape/curvature can in itself mediate the recruitment of lipidated proteins. However, exactly how membrane curvature and composition synergize remains largely unexplored. Here we investigated how three critical structural parameters of lipids, namely acyl chain saturation, headgroup size, and acyl chain length, modulate the capacity of membrane curvature to recruit lipidated proteins. As a model system we used the lipidated minimal membrane anchor of the GTPase, N-Ras (tN-Ras). Our data revealed complex synergistic effects, whereby tN-Ras binding was higher on planar DOPC than POPC membranes, but inversely higher on curved POPC than DOPC membranes. This variation in the binding to both planar and curved membranes leads to a net increase in the recruitment by membrane curvature of tN-Ras when reducing the acyl chain saturation state. Additionally, we found increased recruitment by membrane curvature of tN-Ras when substituting PC for PE, and when decreasing acyl chain length from 14 to 12 carbons (DMPC versus DLPC). However, these variations in recruitment ability had different origins, with the headgroup size primarily influencing tN-Ras binding to planar membranes whereas the change in acyl chain length primarily affected binding to curved membranes. Molecular field theory calculations recapitulated these findings and revealed lateral pressure as an underlying biophysical mechanism dictating how curvature and composition synergize to modulate recruitment of lipidated proteins. Our findings suggest that the different compositions of cellular compartments could modulate the potency of membrane curvature to recruit lipidated proteins and thereby synergistically regulate the trafficking and sorting of lipidated proteins.

On the dynamics of molecular self-assembly and the structural analysis of bilayer membranes using coarse-grained molecular dynamics simulations

We present a molecular dynamics simulation study of the self-assembly of coarse-grained lipid molecules from unbiased random initial configurations. Our lipid model is based on a well-tried CG polymer model with an additional potential that mimics the hydrophobic properties of lipid tails. We find that several stages of self-organization of lipid clusters are involved in the dynamics of bilayer formation and that the resulting equilibrium structures sensitively depend on the strength of hydrophobic interactions hc of the lipid tails and on temperature T. The obtained stable lipid membranes are quantitatively analyzed with respect to their local structure and their degree of order. At equilibrium, we obtain self-stabilizing bilayer membrane structures that exhibit a bending stiffness κB and compression modulus KC comparable to experimental measurements under physiological conditions. We present a phase diagram of our lipid model which covers a sol-gel transition, a liquid (or gel-like) phase including stable bilayer structures and vesicle formation, as well as a quasi-crystalline phase. We also determine the exact conditions for temperature T and degree of hydrophobicity hc for stable bilayer formation including closed vesicles.

Model parameters for simulation of physiological lipids

Coarse grain simulation of proteins in their physiological membrane environment can offer insight across timescales, but requires a comprehensive force field. Parameters are explored for multicomponent bilayers composed of unsaturated lipids DOPC and DOPE, mixed-chain saturation POPC and POPE, and anionic lipids found in bacteria: POPG and cardiolipin. A nonbond representation obtained from multiscale force matching is adapted for these lipids and combined with an improved bonding description of cholesterol. Equilibrating the area per lipid yields robust bilayer simulations and properties for common lipid mixtures with the exception of pure DOPE, which has a known tendency to form nonlamellar phase. The models maintain consistency with an existing lipid-protein interaction model, making the force field of general utility for studying membrane proteins in physiologically representative bilayers.

An overview of molecular dynamics simulations of oxidized lipid systems, with a comparison of ELBA and MARTINI force fields for coarse grained lipid simulations

Biological membranes and model lipid systems containing high amounts of unsaturated lipids and sterols are subject to chemical and/or photo-induced lipid oxidation, which leads to the creation of exotic oxidized lipid products (OxPLs). OxPLs are known to have significant physiological impact in cellular systems and also affect physical properties of both biological and model lipid bilayers. In this paper we (i) provide a perspective on the existing literature on simulations of lipid bilayer systems containing oxidized lipid species as well as the main related experimental results, (ii) describe our new data of all-atom and coarse-grained simulations of hydroperoxidized lipid monolayer and bilayer systems and (iii) provide a comparison of the MARTINI and ELBA coarse grained force fields for lipid bilayer systems. We show that the better electrostatic treatment of interactions in ELBA is able to resolve previous conflicts between experiments and simulations. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Dependence of norfloxacin diffusion across bilayers on lipid composition

Antibiotic resistance is a growing concern in medicine and raises the need to develop and design new drug molecules that can efficiently inhibit bacterial replication. Spurring the passive uptake of the drug molecules is an obvious solution. However our limited understanding of drug-membrane interactions due to the presence of an overwhelming variety of lipids constituting cellular membranes and the lack of facile tools to probe the bio-physical interactions between drugs and lipids imposes a major challenge towards developing new drug molecules that can enter the cell via passive diffusion. Here, we used a label-free micro-fluidic platform combined with giant unilamellar lipid vesicles to investigate the permeability of membranes containing mixtures of DOPE and DOPG in DOPC, leading to a label-free measurement of passive membrane-permeability of autofluorescent antibiotics. A fluoroquinolone drug, norfloxacin was used as a case study. Our results indicate that the diffusion of norfloxacin is strongly dependent on the lipid composition which is not expected from the traditional octanol-lipid partition co-efficient assay. The anionic lipid, DOPG, slows the diffusion process whereas the diffusion across liposomes containing DOPE increases with higher DOPE concentration. Our findings emphasise the need to investigate drug-membrane interactions with focus on the specificity of drugs to lipids for efficient drug delivery, drug encapsulation and targeted drug-delivery.

Predicting Partition Coefficients with a Simple All-Atom/Coarse-Grained Hybrid Model

The solvation free energy is an essential quantity in force field development and in numerous applications. Here, we present the estimation of solvation free energies in polar (water, hexanol, octanol, and nonanol) and in apolar (hexane, octane, and nonane) media. The estimates are produced using molecular dynamics simulations employing a simple all-atom/coarse-grained hybrid model (AA/ELBA) and are therefore very efficient. More than 150 solutes were taken from the Minnesota solvation database and represent small, organic molecules. The mean absolute deviation for the different solvents ranges between 2.0 and 4.1 kJ/mol, and the correlation coefficient ranges between 0.78 and 0.99, indicating that the predictions are accurate. Outliers are identified, and potential avenues for improvements are discussed. Furthermore, partition coefficients between water and the organic solvents were estimated, and the percentage of the predictions that has the correct sign ranges between 74% (for octane) and 92% (for octanol and hexanol). Finally, membrane/water partition coefficients are replaced with hexane/water and octanol/water partition coefficients, and the latter is found to be as accurate as the expensive membrane calculations, indicating a wider application area.

Effects of Lipid Composition on Bilayer Membranes Quantified by All-Atom Molecular Dynamics

Biological bilayer membranes typically contain varying amounts of lamellar and nonlamellar lipids. Lamellar lipids, such as dioleoylphosphatidylcholine (DOPC), are defined by their tendency to form the lamellar phase, ubiquitous in biology. Nonlamellar lipids, such as dioleoylphosphatidylethanolamine (DOPE), prefer instead to form nonlamellar phases, which are mostly nonbiological. However, nonlamellar lipids mix with lamellar lipids in biomembrane structures that remain overall lamellar. Importantly, changes in the lamellar vs nonlamellar lipid composition are believed to affect membrane function and modulate membrane proteins. In this work, we employ atomistic molecular dynamics simulations to quantify how a range of bilayer properties are altered by variations in the lamellar vs nonlamellar lipid composition. Specifically, we simulate five DOPC/DOPE bilayers at mixing ratios of 1/0, 3/1, 1/1, 1/3, and 0/1. We examine properties including lipid area and bilayer thickness, as well as the transmembrane profiles of electron density, lateral pressure, electric field, and dipole potential. While the bilayer structure is only marginally altered by lipid composition changes, dramatic effects are observed for the lateral pressure, electric field, and dipole potential profiles. Possible implications for membrane function are discussed.

A simple and transferable all-atom/coarse-grained hybrid model to study membrane processes

We present an efficient all-atom/coarse-grained hybrid model and apply it to membrane processes. This model is an extension of the all-atom/ELBA model applied previously to processes in water. Here, we improve the efficiency of the model by implementing a multiple-time step integrator that allows the atoms and the coarse-grained beads to be propagated at different timesteps. Furthermore, we fine-tune the interaction between the atoms and the coarse-grained beads by computing the potential of mean force of amino acid side chain analogs along the membrane normal and comparing to atomistic simulations. The model was independently validated on the calculation of small-molecule partition coefficients. Finally, we apply the model to membrane peptides. We studied the tilt angle of the Walp23 and Kalp23 helices in two different model membranes and the stability of the glycophorin A dimer. The model is efficient, accurate, and straightforward to use, as it does not require any extra interaction particles, layers of atomistic solvent molecules or tabulated potentials, thus offering a novel, simple approach to study membrane processes.

A computer model of a polyunsaturated monogalactolipid bilayer

1,2-di-O-acyl-3-O-β-D-galactopyranosyl-sn-glycerol (MGDG) is the main lipid component of thylakoid membranes of higher plants and algae. This monogalactolipid is thought of as a non-bilayer lipid but actually it can form both lamellar and nonlamellar phases. In this study, molecular dynamics (MD) simulations of the fully hydrated di-18:3 MGDG bilayer in the lamellar phase were carried out at 310 and 295 K for 200 and 450 ns, respectively, using the GROMACS 4 software package and OPLS-AA force field. At both temperatures, the lamellar phase of the systems was stable. The pure di-18:3 MGDG bilayer is the first step towards creating a computer model of the lipid matrix of the thylakoid membrane and the main aim of this study was to validate the computer model of di-18:3 MGDG in the bilayer and also to assess the properties of the bilayer. However, only a few of the predicted properties could be compared with those derived experimentally and in other MD simulations because of insufficient amount of such data. Thus, direct validation of the MGDG bilayer proved difficult. Therefore, in the validation process also an indirect approach was used, in which a computer model of the 1,2-dioleoyl-sn-glycero-3-phosphatidylcholine (DOPC) bilayer simulated at the same temperatures using the same force field as the MGDG bilayer was assessed. Successful validation of the DOPC bilayer parameterized in the OPLS-AA force field and similar properties of the MGDG molecules in the pure 18:3 MGDG and in binary 18:3 MGDG-PC bilayers indicate that the computer model of the MDGD molecule is faithful and the MGDG bilayer is representative on the time scales covered in these MD simulations.

Buffer-induced swelling and vesicle budding in binary lipid mixtures of dioleoylphosphatidylcholine:dioleoylphosphatidylethanolamine and dioleoylphosphatidylcholine:lysophosphatidylcholine using small-angle X-ray scattering and ³¹P static NMR

A large variety of data exists on lipid phase behavior; however, it is mostly in nonbuffered systems over nonbiological temperature ranges. We present biophysical data on lipid mixtures of dioleoylphosphatidylcholine (DOPC), dioleoylphosphatidylethanolamine (DOPE), and lysophosphatidylcholine (LysoPC) examining their behaviors in excess water and buffer systems over the temperature range 4-34 °C. These mixtures are commonly used to investigate the effects of spontaneous curvature on integral membrane proteins. Using small-angle X-ray scattering (SAXS) and (31)P NMR, we observed lamellar and vesicle phases, with the buffer causing an increase in the layer spacing. Increasing amounts of DOPE in a DOPC bilayer decreased the layer spacing of the mesophase, while the opposite trend was observed for increasing amounts of LysoPC. (31)P static NMR was used to analyze the DOPC:LysoPC samples to investigate the vesicle sizes present, with evidence of vesicle budding observed at LysoPC concentrations above 30 mol %. NMR line shapes were fitted using an adapted program accounting for the distortion of the lipids within the magnetic field. The distortion of the vesicle, because of magnetic susceptibility, varied with LysoPC content, and a discontinuity was found in both the water and buffer samples. Generally, the distortion increased with LysoPC content; however, at a ratio of DOPC:LysoPC 60:40, the sample showed a level of distortion of the vesicle similar to that of pure DOPC. This implies an increased flexibility in the membrane at this point. Commonly, the assumption is that for increasing LysoPC concentration there is a reduction in membrane tension, implying that estimations of membrane tension based on spontaneous curvature assumptions may not be accurate.

Calculation and visualization of atomistic mechanical stresses in nanomaterials and biomolecules

Many biomolecules have machine-like functions, and accordingly are discussed in terms of mechanical properties like force and motion. However, the concept of stress, a mechanical property that is of fundamental importance in the study of macroscopic mechanics, is not commonly applied in the biomolecular context. We anticipate that microscopical stress analyses of biomolecules and nanomaterials will provide useful mechanistic insights and help guide molecular design. To enable such applications, we have developed Calculator of Atomistic Mechanical Stress (CAMS), an open-source software package for computing atomic resolution stresses from molecular dynamics (MD) simulations. The software also enables decomposition of stress into contributions from bonded, nonbonded and Generalized Born potential terms. CAMS reads GROMACS topology and trajectory files, which are easily generated from AMBER files as well; and time-varying stresses may be animated and visualized in the VMD viewer. Here, we review relevant theory and present illustrative applications.

Direct Mixing of Atomistic Solutes and Coarse-Grained Water

We present a new dual-resolution approach for coupling atomistic and coarse-grained models in molecular dynamics simulations of hydrated systems. In particular, a coarse-grained point dipolar water model is used to solvate molecules represented with standard all-atom force fields. A unique characteristic of our methodology is that the mixing of resolutions is direct, meaning that no additional or ad hoc scaling factors, intermediate regions, or extra sites are required. To validate the methodology, we compute the hydration free energy of 14 atomistic small molecules (analogs of amino acid side chains) solvated by the coarse-grained water. Remarkably, our predictions reproduce the experimental data as accurately as the predictions from state-of-the-art fully atomistic simulations. We also show that the hydration free energy of the coarse-grained water itself is in comparable or better agreement with the experimental value than the predictions from all but one of the most common multisite atomistic models. The coarse-grained water is then applied to solvate a typical atomistic protein containing both α-helix and β-strand elements. Moreover, parallel tempering simulations are performed to investigate the folding free energy landscape of a representative α helical and a β hairpin structure. For the simulations considered in this work, our dual-resolution method is found to be 3 to 6 times more computationally efficient than corresponding fully atomistic approaches.

Preclinical development and ocular biodistribution of gemini-DNA nanoparticles after intravitreal and topical administration: towards non-invasive glaucoma gene therapy

Gene therapy could offer improvement in the treatment of glaucoma compared to the current standard of lowering intraocular pressure. We have developed and characterized non-viral gemini surfactant-phospholipid nanoparticles (GL-NPs) for intravitreal and topical administration. Optimized GL-NPs (size range 150-180 nm) were biocompatible with rat retinal ganglion (RGC-5) cells with >95% viability by PrestoBlue™ assay. GL-NPs carrying Cy5-labeled plasmid DNA demonstrated distinct trafficking behavior and biodisposition within the eye in vivo after intravitreal or topical application with respect to pathways of movement and physicochemical stability. After intravitreal injection in mice, GL-NPs localized within the nerve fiber layer of the retina, whereas after topical application, GL-NPs were located in several anterior chamber tissues, including the limbus, iris and conjunctiva. GL-NPs were thermodynamically stable in the vitreous and tear fluid and were trafficked as single, non-aggregated particles after both types of administration.

FROM THE CLINICAL EDITOR

In this paper, the development and characterization of non-viral gemini surfactant-phospholipid nanoparticles is reported with the goal of establishing a gene delivery system that addresses glaucoma in a non-invasive fashion. The authors found that after topical application, the concentration of these nanoparticles was higher in anterior chamber-related components of the eye, whereas intra-vitreal administration resulted in accumulation in the retinal nerve fibre layer.