Permeation of nanocrystals across lipid membranes

Numerical investigation of force and deflection of nanoneedle penetration into cell using finite element approach: Parameter study and experimental validation of results

This paper aims to develop a numerical methodology to investigate the penetration process of nanoneedles into cells and the corresponding force and indentation length. The finite element approach via the explicit dynamic method handles convergence difficulties in the nonlinear phenomenon. The cell is modeled as an isotropic elastic hemiellipsoidal shell with a thickness of 200 nm, which represents the lipid membrane and actin cortex, encapsulating cytoplasm that is regarded as an Eulerian body because of its fluid-type behavior. Nanoneedles with diameters 400, 200, and 50 nm are considered for model development based on available experimental data. The Von Mises strain failure criterion is used for rupture detection. A parameter study using 1, 2.5, 5, 7.5, and 10 kPa shows that Young's modulus of the HeLa cell membrane is about 5 kPa. Moreover, a failure strain of 1.2 chosen among 0.2, 0.4, 0.6, 0.8, 1, and 1.2 matches best the experimental data. In addition, a diameter study shows that the relations between force-diameter and indentation length-diameter are linear and polynomial, respectively. Furthermore, regarding the experimental data and by using contour of minimum principal stress around needle and an analytical equation for calculation of buckling force of a woven fabric, we proposed that for a given cell, membrane structural stability-a function of the coupled effect of Young's modulus and actin meshwork size-contributes directly to needle insertion success rate for that type of cell.

Shape matters-the interaction of gold nanoparticles with model lung surfactant monolayers

The lung surfactant monolayer (LSM) forms the main biological barrier for any inhaled particles to enter our bloodstream, including gold nanoparticles (AuNPs) present as air pollutants and under investigation for use in biomedical applications. Understanding the interaction of AuNPs with lung surfactant can assist in understanding how AuNPs enter our lungs. In this study, we use coarse-grained molecular dynamics simulations to investigate the effect of four different shape D AuNPs (spherical, box, icosahedron and rod) on the structure and dynamics of a model LSM, with a particular focus on differences resulting from the shape of the AuNP. Monolayer-AuNP systems were simulated in two different states: the compressed state and the expanded state, representing inhalation and exhalation conditions, respectively. Our results indicate that the compressed state is more affected by the presence of the AuNPs than the expanded state. Our results show that in the compressed state, the AuNPs prevent the monolayer from reaching the close to zero surface tension required for normal exhalation. In the compressed state, all four nanoparticles (NPs) reduce the lipid order parameters and cause a thinning of the monolayer where the particles drag surfactant molecules into the water phase. Comparing the different properties shows no trend concerning which shape has the biggest effect on the monolayer, as shape-dependent effects vary among the different properties. Insights from this study might assist future work of how AuNP shapes affect the LSM during inhalation or exhalation conditions.

Computational Studies of Lipid-Wrapped Gold Nanoparticle Transport Through Model Lung Surfactant Monolayers

Colloidal nanoparticles, such as gold nanoparticles (AuNPs), are promising materials for the delivery of hydrophilic drugs via the pulmonary route. The inhaled nanoparticle drug carriers primarily deposit in lung alveoli and interact with the alveolar surface known as lung surfactants. Therefore, it is vital to understand the interactions of nanocarriers with the surfactant layer. To understand the interactions at the molecular level, here we simulated model lung surfactant monolayers with phospholipid (PL)-wrapped AuNPs at the vacuum-water interface using coarse-grained molecular dynamics simulations. The PL-wrapped AuNPs quickly adsorbed into the surfactant layer, altered the structural properties of the monolayer, and at high concentrations initiated the compressed monolayer to collapse/buckle. Among the surfactant monolayer lipid components, cholesterol adsorbed to the AuNPs preferentially over PL species. The position of the adsorbed PL-AuNPs within the monolayer, and subsequent monolayer perturbation, vary depending on the monolayer phase, monolayer composition, and species of PL used as a ligand. Information provided by these molecular dynamic simulations helps to rationalize why some colloidal nanoparticles work better as nanocarriers than others and aid the design of new ones, to avoid biological toxicity and improve efficacy for pulmonary drug delivery.

Mechanistic Understanding From Molecular Dynamics Simulation in Pharmaceutical Research 1: Drug Delivery

In this review, we outline the growing role that molecular dynamics simulation is able to play as a design tool in drug delivery. We cover both the pharmaceutical and computational backgrounds, in a pedagogical fashion, as this review is designed to be equally accessible to pharmaceutical researchers interested in what this new computational tool is capable of and experts in molecular modeling who wish to pursue pharmaceutical applications as a context for their research. The field has become too broad for us to concisely describe all work that has been carried out; many comprehensive reviews on subtopics of this area are cited. We discuss the insight molecular dynamics modeling has provided in dissolution and solubility, however, the majority of the discussion is focused on nanomedicine: the development of nanoscale drug delivery vehicles. Here we focus on three areas where molecular dynamics modeling has had a particularly strong impact: (1) behavior in the bloodstream and protective polymer corona, (2) Drug loading and controlled release, and (3) Nanoparticle interaction with both model and biological membranes. We conclude with some thoughts on the role that molecular dynamics simulation can grow to play in the development of new drug delivery systems.

The role of SP-B1–25 peptides in lung surfactant monolayers exposed to gold nanoparticles

Lung surfactant monolayer’s (acts as the first line barrier for inhaled nanoparticles) components (lipids and peptides) rearrange themselves by the influence of exposed gold nanoparticles at various stages of the breathing cycle.

Molecular-Level "Observations" of the Behavior of Gold Nanoparticles in Aqueous Solution and Interacting with a Lipid Bilayer Membrane

We use coarse-grained molecular dynamics simulations to "observe" details of interactions between ligand-covered gold nanoparticles and a lipid bilayer model membrane. In molecular dynamics simulations, one puts the individual atoms and groups of atoms of the physical system to be "observed" into a simulation box, specifies the forms of the potential energies of interactions between them (ultimately quantum based), and lets them individually move classically according to Newton's equations of motion, based on the forces arising from the assumed potential energy forms. The atoms that are chemically bonded to each other stay chemically bonded, following known potentials (force fields) that permit internal degrees of freedom (internal rotation, torsion, vibrations), and the interactions between nonbonded atoms are simplified to Lennard-Jones forms (in our case) and coulombic (where electrical charges are present) in which the parameters are previously optimized to reproduce thermodynamic properties or are based on quantum electronic calculations. The system is started out at a reasonable set of coordinates for all atoms or groups of atoms, and then permitted to develop according to the equations of motion, one small step (usually 10 fs time step) at a time, for millions of steps until the system is at a quasi-equilibrium (usually reached after hundreds of nanoseconds). We then let the system play out its motions further for many nanoseconds to observe the behavior, periodically taking snapshots (saving all positions and energies), and post-processing the snapshots to obtain various average descriptions of the system. Alkanethiols of various lengths serve as examples of hydrophobic ligands and methyl-terminated PEG with various numbers of monomer units serve as examples of hydrophilic ligands. Spherical gold particles of various diameters as well as gold nanorods form the core to which ligands are attached. The nanoparticles are characterized at the molecular level, especially the distributions of ligand configurations and their dependence on ligand length, and surface coverage. Self-assembly of the bilayer from an isotropic solution and observation of membrane properties that correspond well to experimental values validate the simulations. The mechanism of permeation of a gold NP coated with either a hydrophobic or a hydrophilic ligand, and its dependence on surface coverage, ligand length, core diameter, and core shape, is investigated. Lipid response such as lipid flip-flops, lipid extraction, and changes in order parameter of the lipid tails are examined in detail. The mechanism of permeation of a PEGylated nanorod is shown to occur by tilting, lying down, rotating, and straightening up. The nature of the information provided by molecular dynamics simulations permits understanding of the detailed behavior of gold nanoparticles interacting with lipid membranes which in turn helps to understand why some known systems work better than others and aids the design of new particles and improvement of methods for preparing existing ones.

Simulated Permeation and Characterization of PEGylated Gold Nanoparticles in a Lipid Bilayer System

PEGylated gold nanoparticles are considered suitable nanocarriers for use in biomedical applications and targeted drug delivery systems. In our previous investigation with the alkanethiol-functionalized gold nanoparticle, we found that permeation across a protein-free phospholipid membrane resulted in damaging effects of lipid displacement and water and ion leakage. In the present study, we carry out a series of coarse-grained molecular simulations to explore permeation of lipid bilayer systems by a PEGylated gold nanoparticle, especially at the bulk-liquid-lipid interface as well as the interface between the two lipid leaflets. Initially, we examine molecular-level details of a PEGylated gold nanoparticle (constructed from cycled annealing) in water and find a distribution of ligand configurations (from mushroom to brush states) present in nanoparticles with medium to high surface coverage. We also find that the characteristic properties of the PEGylated gold nanoparticle do not change when it is placed in a salt solution. In our permeation studies, we investigate events of water and ion penetration as well as lipid translocation while varying the ligand length, nanoparticle surface coverage, and ion concentration gradient of our system. Results from our studies show the following: (1) The number of water molecules in the interior of the membrane during ligand-coated nanoparticle permeation increases with PEGn-SH surface coverage, ligand length, and permeation velocity but is not sensitive to the ion concentration gradient. (2) Lipid molecules do not leave the membrane; instead they complete trans-bilayer lipid flip-flop with longer ligands and higher surface coverages. (3) The lack of formation of stable water pores prevents ion translocation. (4) The PEGylated nanoparticle causes less damage to the membrane overall due to favorable interactions with the lipid headgroups which may explain why experimentalists observe endocytosis of PEGylated nanocarriers in vivo.

Shape Effect on Particle-Lipid Bilayer Membrane Association, Cellular Uptake, and Cytotoxicity

Although computer simulation and cell culture experiments have shown that elongated spherical particles can be taken up into cells more efficiently than spherical particles, experimental investigation on effects of these different shapes over the particle-membrane association has never been reported. Therefore, whether the higher cellular uptake of an elongated spherical particles is a result of a better particle-membrane association as suggested by some calculation works or a consequence of its influence on other cellular trans-membrane components involved in particle translocation process, cannot be concluded. Here, we study the effect of particle shape on the particle-membrane interaction by monitoring the association between particles of various shapes and lipid bilayer membrane of artificial cell-sized liposomes. Among the three shaped lanthanide-doped NaYF4 particles, all with high shape purity and uniformity, similar crystal phase, and surface chemistry, the elongated spherical particle shows the highest level of membrane association, followed by the spherical particle with a similar radius, and the hexagonal prism-shaped particle, respectively. The free energy of membrane curvature calculated based on a membrane indentation induced by a particle association indicates that among the three particle shapes, the elongated spherical particle give the most stable membrane curvature. The elongated spherical particles show the highest cellular uptake into cytosol of human melanoma (A-375) and human liver carcinoma (HepG2) cells when observed through a confocal laser scanning fluorescence microscope. Quantitative study using flow cytometry also gives the same result. The elongated spherical particles also possess the highest cytotoxicity in A-375 and normal skin (WI-38) cell lines, comparing to the other two shaped particles.

Theoretical and computational investigations of nanoparticle-biomembrane interactions in cellular delivery

With the rapid development of nanotechnology, nanoparticles have been widely used in many applications such as phototherapy, cell imaging, and drug/gene delivery. A better understanding of how nanoparticles interact with bio-system (especially cells) is of great importance for their potential biomedical applications. In this review, the current status and perspective of theoretical and computational investigations is presented on the nanoparticle-biomembrane interactions in cellular delivery. In particular, the determining parameters (including the properties of nanoparticles, cell membranes and environments) that govern the cellular uptake of nanoparticles (direct penetration and endocytosis) are discussed. Further, some special attention is paid to their interactions beyond the translocation of nanoparticles across membranes (e.g., nanoparticles escaping from endosome and entering into nucleus). Finally, a summary is given, and the challenging problems of this field in the future are identified.

Surface-functionalized nanoparticle permeation triggers lipid displacement and water and ion leakage

Functionalized nanoparticles (NPs) are considered suitable carriers for targeted drug delivery systems. However, the ion and water leakage induced by permeation of these nanoparticles is a challenge in these drug delivery methods because of cytotoxic effects of some ions. In this study, we have carried out a series of coarse-grained molecular dynamics simulations to investigate the effect of length of ligands on permeation of a nanoparticle across a protein-free phospholipid bilayer membrane. Water and ion penetration as well as incidence of lipid flip-flop events and loss of lipid molecules from the membrane are explored in this study while varying the nanoparticle size, length of ligand, ion concentration gradient, pressure differential across the membrane, and nanoparticle permeation velocity. Some results from our studies include (1) the number of water molecules in the interior of the membrane during ligand-coated nanoparticle permeation increases with nanoparticle size, ligand length, pressure differential, and permeation velocity but is not sensitive to the ion concentration gradient; (2) some lipid molecules leave the membrane by being entangled with ligands of the NP instead of completing the flip-flop that permits them to rejoin the membrane, thereby leading to fewer flip-flop events; and (3) the formation of water columns or water "fingers" provides a mechanism of ion transport across lipid bilayer membranes, but such ion penetration events are less likely for sodium ions than chloride ions and less likely for nanoparticles with longer-ligands.

Electrical method to quantify nanoparticle interaction with lipid bilayers

Understanding as well as rapidly screening the interaction of nanoparticles with cell membranes is of central importance for biological applications such as drug and gene delivery. Recently, we have shown that "striped" mixed-monolayer-coated gold nanoparticles spontaneously penetrate a variety of cell membranes through a passive pathway. Here, we report an electrical approach to screen and readily quantify the interaction between nanoparticles and bilayer lipid membranes. Membrane adsorption is monitored through the capacitive increase of suspended planar lipid membranes upon fusion with nanoparticles. We adopt a Langmuir isotherm model to characterize the adsorption of nanoparticles by bilayer lipid membranes and extract the partition coefficient, K, and the standard free energy gain by this spontaneous process, for a variety of sizes of cell-membrane-penetrating nanoparticles. We believe that the method presented here will be a useful qualitative and quantitative tool to determine nanoparticle interaction with lipid bilayers and consequently with cell membranes.

Nanoparticle permeation induces water penetration, ion transport, and lipid flip-flop

Nanoparticles are generally considered excellent candidates for targeted drug delivery. However, ion leakage and cytotoxicity induced by nanoparticle permeation is a potential problem in such drug delivery schemes because of the toxic effect of many ions. In this study, we have carried out a series of coarse-grained molecular dynamics simulations to investigate the water penetration, ion transport, and lipid molecule flip-flop in a protein-free phospholipid bilayer membrane during nanoparticle permeation. The effect of ion concentration gradient, pressure differential across the membrane, nanoparticle size, and permeation velocity have been examined in this work. Some conclusions from our studies include (1) The number of water molecules in the interior of the membrane during the nanoparticle permeation increases with the nanoparticle size and the pressure differential across the membrane but is unaffected by the nanoparticle permeation velocity or the ion concentration gradient. (2) Ion transport is sensitive to the size of nanoparticle as well as the ion concentration gradient between two sides of the membrane; no anion/cation selectivity is observed for small nanoparticle permeation, while anions are preferentially translocated through the membrane when the size of nanoparticle is large enough. (3) Incidences of lipid molecule flip-flop increases with the size of nanoparticle and ion concentration gradient and decreases with the pressure differential and the nanoparticle permeation velocity.

Molecular dynamics simulations of creatine kinase and adenine nucleotide translocase in mitochondrial membrane patch

Interaction between mitochondrial creatine kinase (MtCK) and adenine nucleotide translocase (ANT) can play an important role in determining energy transfer pathways in the cell. Although the functional coupling between MtCK and ANT has been demonstrated, the precise mechanism of the coupling is not clear. To study the details of the coupling, we turned to molecular dynamics simulations. We introduce a new coarse-grained molecular dynamics model of a patch of the mitochondrial inner membrane containing a transmembrane ANT and an MtCK above the membrane. The membrane model consists of three major types of lipids (phosphatidylcholine, phosphatidylethanolamine, and cardiolipin) in a roughly 2:1:1 molar ratio. A thermodynamics-based coarse-grained force field, termed MARTINI, has been used together with the GROMACS molecular dynamics package for all simulated systems in this work. Several physical properties of the system are reproduced by the model and are in agreement with known data. This includes membrane thickness, dimension of the proteins, and diffusion constants. We have studied the binding of MtCK to the membrane and demonstrated the effect of cardiolipin on the stabilization of the binding. In addition, our simulations predict which part of the MtCK protein sequence interacts with the membrane. Taken together, the model has been verified by dynamical and structural data and can be used as the basis for further studies.