Spontaneous insertion of carbon nanotube bundles inside biomembranes: A hybrid particle-field coarse-grained molecular dynamics study

Bacterial lipids drive compartmentalization on the nanoscale

The design of cellular functions in synthetic systems, inspired by the internal partitioning of living cells, is a constantly growing research field that is paving the way to a large number of new remarkable applications. Several hierarchies of internal compartments like polymersomes, liposomes, and membranes are used to control the transport, release, and chemistry of encapsulated species. However, the experimental characterization and the comprehension of glycolipid mesostructures are far from being fully addressed. Lipid A is indeed a glycolipid and the endotoxic part of Gram-negative bacterial lipopolysaccharide; it is the moiety that is recognized by the eukaryotic receptors giving rise to the modulation of innate immunity. Herein we propose, for the first time, a combined approach based on hybrid Particle-Field (hPF) Molecular Dynamics (MD) simulations and Small Angle X-Ray Scattering (SAXS) experiments to gain a molecular picture of the complex supramolecular structures of lipopolysaccharide (LPS) and lipid A at low hydration levels. The mutual support of data from simulations and experiments allowed the unprecedented discovery of the presence of a nano-compartmentalized phase composed of liposomes of variable size and shape which can be used in synthetic biological applications.

The role of size and nature in nanoparticle binding to a model lung membrane: an atomistic study

Understanding the uptake of nanoparticles (NPs) by different types of cellular membranes plays a pivotal role in the design of NPs for medical applications and in avoiding adverse effects that result in nanotoxicity. Yet, the role of key design parameters, such as the bare NP material, NP size and surface reactivity, and the nature of NP coatings, in membrane remodelling and uptake mechanisms is still very poorly understood, particularly towards the lower range of NP dimensions that are beyond the experimental imaging resolution. The same can be said about the role of a particular membrane composition. Here, we systematically employ biased and unbiased molecular dynamics simulations to calculate the binding energy for three bare materials (Ag/SiO₂/TiO₂) and three NP sizes (1/3/5 nm diameter) with a representative lung surfactant membrane, and to study their binding kinetics. The calculated binding energies show that irrespective of size, Ag nanoparticles bind very strongly to the bilayer, while the NPs made of SiO₂ or TiO₂ experience very low to no binding. The unbiased simulations provide insight into how the NPs and membrane affect each other in terms of the solvent-accessible surface area (SASA) of the NPs and the defect types and fluidity of the membrane. Using these systematic fine-grained results in coarsening procedures will pave the way for simulations considering NP sizes that are well beyond the membrane thickness, i.e. closer to experimental dimensions, for which different binding characteristics and more significant membrane remodelling are expected.

Computer modeling of polymer stars in variable solvent conditions: a comparison of MD simulations, self-consistent field (SCF) modeling and novel hybrid Monte Carlo SCF approach

Computer-aided modeling is a systematic approach to grasp the physics of macromolecules, but it remains essential to know when to trust the results and when not. For a polymer star, we consider three approaches: (i) Molecular Dynamics (MD) simulations and implementing a coarse-grained model, (ii) the self-consistent field approach based on a mean-field approximation and implementing the lattice model due to Scheutjens and Fleer (SF-SCF) and (iii) novel hybrid Monte Carlo self-consistent field (MC-SCF) method, which combines a coarse-grained model driven by a Monte Carlo method and a mean-field representation driven by SF-SCF. We compare the performance of these approaches under a wide range of solvent qualities. The MD approach is formally the most exact but suffers from reasonable convergence. The mean-field approach works similarly in all solvent qualities but is quantitatively least accurate. The MC-SCF hybrid allows us to combine the benefits of the simulation route and the effective performance of SCF. We consider the center-to-end distance R_ce, the radius of gyration R_g² of the star and the polymer density profiles φ(r) of polymer-segments in it. All three methods show a good qualitative agreement one to another. The MC-SCF method is in good agreement with the scaling predictions in the whole range of solvent quality values showing that it grasps the essential physics while remaining computationally in bounds.

Investigating the interactions of corona-free SWCNTs and cell membrane models using sum-frequency generation

The understanding of the interactions between biomolecules and nanomaterials is of great importance in many areas of nanomedicine and bioapplications. Numerous studies in this area have been performed. However, toxicological aspects involving the interaction between phospholipids and carbon nanotubes (CNTs) remain undefined, especially for those cases in which a protein corona is not formed around the nanomaterial (corona-free nanomaterials). This study focuses on the interaction of Langmuir films of dipalmitoylphosphatidylglycerol (DPPG) and dipalmitoylphosphatidylcholine (DPPC) with corona-free, single-walled CNTs. Surface pressure-area isotherms and sum-frequency generation (SFG) vibrational spectroscopy, a non-linear optical technique used to study surfaces and interfaces, were used to investigate the lipid tail orientation and conformation, aiming to understand the interactions between phospholipids and single walled carbon nanotubes functionalized by carboxylic acid (SWCNTs-COOH) at the air-water interface under low ionic strength conditions. Data from isotherms and SFG spectra revealed that the SWCNT adsorption at the air-water interface is induced by the presence of both lipids, although at a lesser extent for DPPG due to its anionic head group, which could result in repulsion of SWCNTs-COOH that also bear a negative charge. Furthermore, lipid monolayers remained conformationally ordered, indicating insertion of SWCNTs into the lipid monolayer. Our results corroborate previous works and simulations in the literature, but made it possible to perform an in-depth investigation of the interaction of these nanomaterials with components of phospholipid membranes.

Aggregation of Lipid A Variants: A Hybrid Particle-Field Model

Lipid A is one of the three components of bacterial lipopolysaccharides constituting the outer membrane of Gram-negative bacteria, and is recognized to have an important biological role in the inflammatory response of mammalians. Its biological activity is modulated by the number of acyl-chains that are present in the lipid and by the dielectric medium, i.e., the type of counter-ions, through electrostatic interactions. In this paper, we report on a coarse-grained model of chemical variants of Lipid A based on the hybrid particle-field/molecular dynamics approach (hPF-MD). In particular, we investigate the stability of Lipid A bilayers for two different hexa- and tetra-acylated structures. Comparing particle density profiles along bilayer cross-sections, we find good agreement between the hPF-MD model and reference all-atom simulation for both chemical variants of Lipid A. hPF-MD models of constituted bilayers composed by hexa-acylated Lipid A in water are stable within the simulation time. We further validate our model by verifying that the phase behavior of Lipid A/counterion/water mixtures is correctly reproduced. In particular, hPF-MD simulations predict the correct self-assembly of different lamellar and micellar phases from an initially random distribution of Lipid A molecules with counterions in water. Finally, it is possible to observe the spontaneous formation and stability of Lipid A vesicles by fusion of micellar aggregates.

Nanoparticles at Biomimetic Interfaces: Combined Experimental and Simulation Study on Charged Gold Nanoparticles/Lipid Bilayer Interfaces

The poor understanding of the interaction of nanomaterials with biologically relevant interfaces is recognized as one of the major issues currently limiting the development of nanomedicine. The central purpose of this study is to compare experimental (confocal microscopy, fluorescence correlation spectroscopy, X-ray reflectivity) and computational (molecular dynamics simulations) results to thoroughly describe the interaction of cationic gold nanoparticles (AuNPs) with mixed zwitterionic/anionic lipid membranes. The adhesion of AuNPs to the lipid membrane is investigated on different length scales from a structural and dynamical point of view; with this approach, a series of complex phenomena, spanning from lipid extraction, localized membrane disruption, lateral phase separation, and slaved diffusion, are characterized and interpreted from a molecular level to macroscopic observations.

Molecular dynamics simulations of carbon nanotube porins in lipid bilayers

Carbon nanotube porins embedded in lipid membranes are studied by molecular dynamics simulations.

Molecular dynamics simulation strategies for designing carbon-nanotube-based targeted drug delivery

The carbon nanotube (CNT)-based target-specific delivery of drugs, or other molecular cargo, has emerged as one of the most promising biomedical applications of nanotechnology. To achieve efficient CNT-based drug delivery, the interactions between the drug, CNT and biomolecular target need to be properly optimized. Recent advances in the computer-aided molecular design tools, in particular molecular dynamics (MD) simulation studies, offer an appropriate low-cost approach for such optimization. This review highlights the various potential MD approaches for the simulation of CNT interactions with cell membranes while emphasizing various methods of cellular internalization and toxicities of CNTs to build new strategies for designing rational CNT-based targeted drug delivery to circumvent the limitations associated with the various clinically available nonspecific therapeutic agents.

Carbon Nanotubes Mediate Fusion of Lipid Vesicles

The fusion of lipid membranes is opposed by high energetic barriers. In living organisms, complex protein machineries carry out this biologically essential process. Here we show that membrane-spanning carbon nanotubes (CNTs) can trigger spontaneous fusion of small lipid vesicles. In coarse-grained molecular dynamics simulations, we find that a CNT bridging between two vesicles locally perturbs their lipid structure. Their outer leaflets merge as the CNT pulls lipids out of the membranes, creating an hourglass-shaped fusion intermediate with still intact inner leaflets. As the CNT moves away from the symmetry axis connecting the vesicle centers, the inner leaflets merge, forming a pore that completes fusion. The distinct mechanism of CNT-mediated membrane fusion may be transferable, providing guidance in the development of fusion agents, e.g., for the targeted delivery of drugs or nucleic acids.

Surface-Enhanced Raman Scattering to Evaluate Nanomaterial Cytotoxicity on Living Cells

The increasing number of reports about false positive or negative results from conventional cytotoxicity assays of nanomaterials (NMs) suggests that more reliable NM toxicity assessment methods should be developed. Here, we report a novel approach for nanotoxicity evaluation based on surface-enhanced Raman spectroscopy (SERS). Three model NMs were tested on two model cell lines and the results were validated by WST-1 cytotoxicity assay and annexin V-FITC/propidium iodide (PI) staining as apoptosis-necrosis assay. The localization of nanoparticles (NPs) in the cells and the cellular conditions upon NP incubation were visualized by transmission electron microscopy (TEM) and enhanced dark-field (EDF) microscopy. SERS revealed a broader view on the consequences of cell-NM interactions compared to the conventional cytotoxicity assays where only one aspect of toxicity can be measured by one assay type. The results suggest that SERS can significantly contribute to the cytotoxicity evaluation bypassing NM or assay component-related complications with less effort.

Structure of Carbon Nanotube Porins in Lipid Bilayers: An in Situ Small-Angle X-ray Scattering (SAXS) Study

Carbon nanotube porins (CNTPs), small segments of carbon nanotubes capable of forming defined pores in lipid membranes, are important future components for bionanoelectronic devices as they could provide a robust analog of biological membrane channels. In order to control the incorporation of these CNT channels into lipid bilayers, it is important to understand the structure of the CNTPs before and after insertion into the lipid bilayer as well as the impact of such insertion on the bilayer structure. Here we employed a noninvasive in situ probe, small-angle X-ray scattering, to study the integration of CNT porins into dioleoylphosphatidylcholine bilayers. Our results show that CNTPs in solution are stabilized by a monolayer of lipid molecules wrapped around their outer surface. We also demonstrate that insertion of CNTPs into the lipid bilayer results in decreased bilayer thickness with the magnitude of this effect increasing with the concentration of CNTPs.

Erratum. Computational insight in the role of fusogenic lipopeptides at the onset of liposome fusion

We performed an extensive computational study to obtain insight in the molecular mechanisms that take place prior to membrane fusion. We focused on membrane-anchored hybrid macromolecules (lipid–polymer–oligopeptide) that mimic biological SNARE proteins in terms of liposome fusion characteristics [H. Robson Marsden et al., 2009]; efficient micro-second simulation was enabled by combining validated MARTINI force fields for the molecular building blocks in coarse-grained molecular dynamics (CGMD). We find that individual peptide domains in the hybrid macromolecules bind and partially integrate parallel to the membrane surface, in agreement with experimental findings. By varying several experimental design parameters, we observe that peptide domains remain in the solvent phase only in two cases: (1) for solitary lipopeptides (low concentration), below a threshold area per lipid in the membrane, and (2) when the lipopeptide concentration is high enough for the peptide domains to self-assemble into tetrameric homo-complexes. The peptide-membrane binding is not affected by solvent-induced peptide unfolding, which we mimicked by relaxing the usual MARTINI helix constraints. Remarkably, in this case, a reverse transition to a helical secondary structure is observed after binding, highlighting the role of the membrane as a template (partitioning-folding coupling). Our findings undermine the current view of the initial stages towards fusion, in which membranes are thought to be kept in close apposition via dimerization of individual complementary peptides in the solvent phase. Although we did not study actual fusion, our simulations show that the formation of homomers, which is suppressed in experimental peptide pair design and therefore believed to be insignificant for fusion, by peptides anchored to the same membrane does play a key role in this locking mechanism and potentially also in membrane destabilization that precedes fusion.

Computational insight in the role of fusogenic lipopeptides at the onset of liposome fusion

We performed an extensive computational study to obtain insight in the molecular mechanisms that take place prior to membrane fusion. We focused on membrane-anchored hybrid macromolecules (lipid-polymer-oligopeptide) that mimic biological SNARE proteins in terms of liposome fusion characteristics [H. Robson Marsden et al., 2009]; efficient micro-second simulation was enabled by combining validated MARTINI force fields for the molecular building blocks in coarse-grained molecular dynamics (CGMD). We find that individual peptide domains in the hybrid macromolecules bind and partially integrate parallel to the membrane surface, in agreement with experimental findings. By varying several experimental design parameters, we observe that peptide domains remain in the solvent phase only in two cases: (1) for solitary lipopeptides (low concentration), below a threshold area per lipid in the membrane, and (2) when the lipopeptide concentration is high enough for the peptide domains to self-assemble into tetrameric homo-complexes. The peptide-membrane binding is not affected by solvent-induced peptide unfolding, which we mimicked by relaxing the usual MARTINI helix constraints. Remarkably, in this case, a reverse transition to a helical secondary structure is observed after binding, highlighting the role of the membrane as a template (partitioning-folding coupling). Our findings undermine the current view of the initial stages towards fusion, in which membranes are thought to be kept in close apposition via dimerization of individual complementary peptides in the solvent phase. Although we did not study actual fusion, our simulations show that the formation of homomers, which is suppressed in experimental peptide-pair design and therefore believed to be insignificant for fusion, by peptides anchored to the same membrane does play a key role in this locking mechanism and potentially also in membrane destabilization that precede fusion.

Generation of Well-Relaxed All-Atom Models of Large Molecular Weight Polymer Melts: A Hybrid Particle-Continuum Approach Based on Particle-Field Molecular Dynamics Simulations

A procedure based on Molecular Dynamics (MD) simulations employing soft potentials derived from self-consistent field (SCF) theory (named MD-SCF) able to generate well-relaxed all-atom structures of polymer melts is proposed. All-atom structures having structural correlations indistinguishable from ones obtained by long MD relaxations have been obtained for poly(methyl methacrylate) (PMMA) and poly(ethylene oxide) (PEO) melts. The proposed procedure leads to computational costs mainly related on system size rather than to the chain length. Several advantages of the proposed procedure over current coarse-graining/reverse mapping strategies are apparent. No parametrization is needed to generate relaxed structures of different polymers at different scales or resolutions. There is no need for special algorithms or back-mapping schemes to change the resolution of the models. This characteristic makes the procedure general and its extension to other polymer architectures straightforward. A similar procedure can be easily extended to the generation of all-atom structures of block copolymer melts and polymer nanocomposites.