Modeling the Hydrophobicity of Nanoparticles and Their Interaction with Lipids and Proteins

Development of a bottom-up coarse-grained model for interactions of lipids with TiO 2 nanoparticles

Understanding interactions of inorganic nanoparticles with biomolecules is important in many biotechnology, nanomedicine, and toxicological research, however, the size of typical nanoparticles makes their direct modeling by atomistic simulations unfeasible. Here, we present a bottom-up coarse-graining approach for modeling titanium dioxide (TiO  2 ) nanomaterials in contact with phospholipids that uses the inverse Monte Carlo method to optimize the effective interactions from the structural data obtained in small-scale all-atom simulations of TiO  2 surfaces with lipids in aqueous solution. The resulting coarse-grained models are able to accurately reproduce the structural details of lipid adsorption on different titania surfaces without the use of an explicit solvent, enabling significant computational resource savings and favorable scaling. Our coarse-grained simulations show that small spherical TiO  2 nanoparticles ( r = 2 nm) can only be partially wrapped by a lipid bilayer with phosphoethanolamine headgroups, however, the lipid adsorption increases with the radius of the nanoparticle. The current approach can be used to study the effect of the size and shape of TiO  2 nanoparticles on their interactions with cell membrane lipids, which can be a determining factor in membrane wrapping as well as the recently discovered phenomenon of nanoquarantining, which involves the formation of layered nanomaterial-lipid structures.

How Does pH Affect the Adsorption of Human Serum Protein in the Presence of Hydrophobic and Hydrophilic Nanoparticles at Air-Water and Lipid-Water Interfaces?

This study investigates interaction between hydrophilic (11-mercaptoundecanoic acid (MUA)) and hydrophobic (1-undecanethiol (UDT)) gold nanoparticles (GNPs) with human serum albumin (HSA) protein on air-water and lipid-water interfaces at pH 3 and 7. Vibrational sum frequency generation (VSFG) spectroscopy is used to analyze changes in the intensity of interfacial water molecules and the C-H group of the protein. At the air-water interface, the hydrophobic interaction between the HSA protein and hydrophobic GNPs at pH 3 leads to their accumulation at the interface, resulting in an increased C-H intensity of the protein with a slight decrease in water intensity. Whereas, at pH 7, where the negative charge of the protein results in the reduced surface activity of the HSA compared to pH 3, the interaction between alkyl chain of the hydrophobic GNPs and alkyl group of the protein results in the adsorption of the protein-capped GNPs at the interface. This leads to an increased intensity of the C-H group of protein and water molecules. However, negatively charged hydrophilic GNPs do not induce significant changes in the interfacial water structure or the C-H group of the protein due to the electrostatic force of repulsion with the negatively charged HSA at pH 7. In contrast, at the lipid-water interface, both hydrophobic and hydrophilic GNPs interact with HSA protein, causing disordering of interfacial water molecules at pH 3 and ordering at pH 7. Interestingly, similar behavior of the protein with both types of GNPs results in comparable ordering/disordering at the interface depending on the pH of solution. Furthermore, the VSFG results obtained with the deuterated lipid suggest that changes in ordering and disorder occur due to increased protein adsorption in the presence of GNPs, causing alterations in the membrane structure. These findings give a better understanding of the mechanisms that govern protein-nanoparticle interaction and their consequential effects on the structure, function, and behavior of molecules at the biological membrane interface, which is crucial for developing safe and effective nanoparticle-based therapeutics.

Coarse-grained molecular models of the surface of hair

We present a coarse-grained molecular model of the surface of human hair, which consists of a supported lipid monolayer, in the MARTINI framework. Using coarse-grained molecular dynamics (MD) simulations, we identify a lipid grafting distance that yields a monolayer thickness consistent with both atomistic MD simulations and experimental measurements of the hair surface. Coarse-grained models for fully-functionalised, partially damaged, and fully damaged hair surfaces are created by randomly replacing neutral thioesters with anionic sulfonate groups. This mimics the progressive removal of fatty acids from the hair surface by bleaching and leads to chemically heterogeneous surfaces. Using molecular dynamics (MD) simulations, we study the island structures formed by the lipid monolayers at different degrees of damage in vacuum and in the presence of polar (water) and non-polar (n-hexadecane) solvents. We also use MD simulations to compare the wetting behaviour of water and n-hexadecane droplets on the model surfaces through contact angle measurements, which are compared to experiments using virgin and bleached hair. The model surfaces capture the experimentally-observed transition of the hair surface from hydrophobic (and oleophilic) to hydrophilic (and oleophobic) as the level of bleaching damage increases. By selecting surfaces with specific damage ratios, we obtain contact angles from the MD simulations that are in good agreement with experiments for both solvents on virgin and bleached human hairs. To negate the possible effects of microscale curvature and roughness of real hairs on wetting, we also conduct additional experiments using biomimetic surfaces that are co-functionalised with fatty acids and sulfonate groups. In both the MD simulations and experiments, the cosine of the water contact angle increases linearly with the sulfonate group surface coverage with a similar slope. We expect that the proposed systems will be useful for future molecular dynamics simulations of the adsorption and tribological behaviour of hair, as well as other chemically heterogeneous surfaces.

Scission energy and topology of micelles controlled by the molecular structure of additives

We employ coarse-grained (CG) molecular dynamics simulations (MD) to investigate the effects of the molecular structure of additives on the scission energy and morphology of charged micelles. Considering sodium dodecyl sulfate (SDS) as a representative charged surfactant and taking trimethylphenylammonium chloride (TMPAC) and octyltrimethylammonium bromide (OTAB) as oppositely charged additives, we show that the scission energy and topology of micelles vary significantly depending on the molecular structure of the hydrophobic part of the additives. The cyclic aromatic tail of the TMPAC disrupts the core structure of the SDS micelle and hence decreases the micelle scission energy, whereas the linear alkyl tail of the OTAB packs very well with the micelle core and increases the scission energy. Although both the additives have similar head structures, they lead to very different micelle morphologies because of the difference in the shape of their tail structures; ring-like or toroidal shaped micelles are formed in SDS/TMPAC solution whereas bicelle-like structures are formed in SDS/OTAB solution when the additive to surfactant ratio is higher than a certain value.

Gecko adhesion: a molecular-simulation perspective on the effect of humidity

Gecko adhesion is investigated by molecular dynamics simulations. It is known, that the gecko adhesion system shows increased pull-off forces in humid environments. A coarse-grained model of gecko beta keratin, previously developed in our group, is extended and used to elucidate the molecular mechanisms involved in this humidity effect on adhesion. We show that neither the change of the elastic properties of gecko keratin, nor capillary forces, can solely explain the increased pull-off forces of wet gecko keratin. Instead, we establish a molecular picture of gecko adhesion where the interplay between capillary bridges and a mediator effect of water, enhances pull-off forces, consistent with observations in AFM experiments at high humidities. We find that water density is raised locally, in molecular scale asperities of the gecko keratin and that this increase in local water density smoothes the surface-spatula interface. Water, which is absorbed into the keratin, acts as a mediator, and leads during pull-off to the dominant contribution in the van der Waals energy, because the dispersion interactions between water and surface are primarily opposing the pull-off.

Recent advances in fluorescent probes for lipid droplets

Lipid droplets (LDs) have been known as a non-negligible cellular organelle for lipid storage and metabolism. Fluorescent probes for imaging LDs would be paramount for depicting their functions in cells. Although commercially available Nile Red and BODIPY^tm 493/503 have been widely used for labelling LDs, they exhibit unsatisfactory specificity and spectroscopic properties. This feature article reviews the recent advances in organic fluorescent probes for imaging LDs. We first introduce the key points for probe design, including regulating hydrophobicity and enhancing fluorescence quantum yield in LDs. Then, we summarize the structural features and biological applications of some representative LD probes classified by their frameworks. In addition, the current challenges and future research trends for the fluorescent probes of LDs are discussed as well.

Are Plant-Based Carbohydrate Nanoparticles Safe for Inhalation? Investigating Their Interactions with the Pulmonary Surfactant Using Langmuir Monolayers

Nanoparticle carriers show promise for drug delivery, including by inhalation, where the first barrier for uptake in the lungs is the monolayer pulmonary surfactant membrane that coats the air/alveoli interface and is critical to breathing. It is imperative to establish the fate of potential nanocarriers and their effects on the biophysical properties of the pulmonary surfactant. To this end, the impact of the nanoparticle surface charge on the lateral organization, thickness, and recompressibility of Langmuir monolayers of model phospholipid-only and phospholipid-protein mixtures was investigated using native and modified forms of nanophytoglycogen, a carbohydrate-based dendritic polymer extracted from corn as monodisperse nanoparticles. We show that the native (quasi-neutral) and anionic nanophytoglycogens have little impact on the phase behavior and film properties. By contrast, cationic nanophytoglycogen alters the film morphology and increases the hysteresis associated with the work of breathing due to its electrostatic interaction with the anionic phospholipids in the model systems. These findings specifically highlight the importance of surface charge as a selection criterion for inhaled nanoformulations.

Computational investigation of a switchable emulsion stabilized by the mixture of a surfactant and tertiary amine

Molecular dynamics simulations were performed to investigate the CO2-responsiveness of an oil-in-water (O/W) emulsion stabilized by sodium oleate (NaOA) with a tertiary amine additive, named pentamethyl diethylenetriamine (PMA). The simulated results were in accordance with the experimental observations. That is, the surfactant NaOA itself can stabilize dodecane/water emulsions in aqueous solution, while the CO2-reponsiveness was strongly related to the added PMA. The electroneutral PMA molecules preferred to be located in the core region of the droplets. Thus, under the same conditions, the size of the droplet containing PMA is predictably larger than that without PMA. The increased extent of the charged surfactant headgroups distribution can increase the electrostatic repulsion between the droplets in the emulsion solution, which is the important reason why a much more stable emulsion is obtained by adding PMA. When PMA molecules were protonated to PMA2+ by bubbling CO2, they migrated from the interior to the surface of the droplets under electrostatic attraction, forming ion pairs with OA-. The binding between PMA2+ and OA- made the distribution of the surfactants very concentrated on the droplet surface, leading to large hydrophobic areas exposed to water. Besides, the hydration interactions of OA- headgroups decreased because they were covered by PMA2+. The calculated potential of mean force (PMF) confirmed that the electrostatic repulsion between droplets was crucial for the emulsion stabilization.

Effect of concentration of PEG coated gold nanoparticle on lung surfactant studied with coarse-grained molecular dynamics simulations

The surface modification of nanoparticles can not only change the physical and chemical properties of particles, such as the hydrophilic and hydrophobic properties and surface charges of nanoparticles to a certain extent, but also bring new functions to nanoparticles, such as membrane permeability and targeting. Inhaled nanoparticles (NPs) are experienced by the first biological barrier inside the alveolus known as lung surfactant (LS), consisting of phospholipids and proteins in the form of the monolayer at the air-water interface. Inhaled NPs can reach deep into the lungs and interfere with the biophysical properties of the lung components. The interaction mechanisms of bare gold nanoparticles (AuNPs) with the LS monolayer are not well understood. Coarse-grained molecular dynamics simulations were carried out to have a study on the interactions of PEG coated AuNPs with LS monolayers. It was observed that the interactions of AuNPs and LS components make the monolayer structure deform and change the biophysical properties of LS monolayer. The results also indicate that AuNPs with high concentrations hinder the lowering of the LS surface tension and reduce lateral mobility of lipids. Overall, the simulation results can provide guidance for the design of ligand protected NPs as drug carriers and can identify the nanoparticles potential side effect on lung surfactant.

Size-dependent interaction of hydrophilic/hydrophobic ligand functionalized cationic and anionic nanoparticles with lipid bilayers

We study the nature of nanoparticle (NPs)-membrane interaction as a function of nanoparticle size for different functionalization using molecular dynamics simulation. Zinc sulphide quantum dots of size, 2 nm and 4 nm are used as model NPs, and DLPC and DPPC lipid bilayers are used as model membranes. We use coarse-grained polarizable MARTINI model (MPW) to simulate the NPs and lipid bilayers. Our simulation results show that uncharged bare NPs penetrate the lipid bilayers and embed themselves within the hydrophobic core of the bilayer both in the gel and fluid phases. NPs of size 4 nm are shown to disrupt the bilayer. The bilayer recovers from the damages caused by smaller NPs of size 2 nm. In case of either purely hydrophilic or hybrid (with hydrophilic/hydrophobic ratio of 2:1) ligand-functionalized NPs of smaller size (shell size 2 nm), only cationic NPs bind to the bilayer. However, for larger NPs with a shell size of 4 nm, both anionic and cationic hybrid functionalized NPs bind to the bilayer. The performance of standard Martini (SM) force field for the charged NP/bilayer systems has also been tested and compared with the results obtained using MPW model. Although the overall trend that the cationic NPs interact strongly with the bilayers than their anionic counterparts has been captured correctly using SM, the adsorption behaviour of the functionalized NPs differ significantly in the SM force field. The interaction of anionic NPs with both fluid and gel bilayers has been observed to be least accurately represented in the SM force field.

Molecular Dynamics Simulation of Pure n-Alkanes and Their Mixtures at Elevated Temperatures Using Atomistic and Coarse-Grained Force Fields

The properties of higher n-alkanes and their mixtures is a topic of significant interest for the oil and chemical industry. However, the experimental data at high temperatures are scarce. The present study focuses on simulating n-dodecane, n-octacosane, their binary mixture at a n-dodecane mole fraction of 0.3, and a model mixture of the commercially available hydrocarbon wax SX-70 to evaluate the performance of several force fields on the reproduction of properties such as liquid densities, surface tension, and viscosities. Molecular dynamics simulations over a broad temperature range from 323.15 to 573.15 K were employed in examining a broad set of atomistic molecular models assessed for the reproduction of experimental data. The well-established united atom TraPPE (TraPPE-UA) was compared against the all atom optimized potentials for liquid simulations (OPLS) reparametrization for long n-alkanes, L-OPLS, as well as Lipid14 and MARTINI force fields. All models qualitatively reproduce the temperature dependence of the aforementioned properties, but TraPPE-UA was found to reproduce liquid densities most accurately and consistently over the entire temperature range. TraPPE-UA and MARTINI were very successful in reproducing surface tensions, and L-OPLS was found to be the most accurate in reproducing the measured viscosities as compared to the other models. Our simulations show that these widely used force fields originating from the world of biomolecular simulations are suitable candidates in the study of n-alkane properties, both in the pure and mixture states.

Molecular mechanism of HIV-1 TAT peptide and its conjugated gold nanoparticles translocating across lipid membranes

The trans-acting activator of transcription (TAT) peptide, which is derived from human immunodeficiency virus-1 (HIV-1), has been widely used as an effective nanocarrier to transport extracellular substances into cells. However, the underlying translocation mechanism of TAT peptide across cell membranes still remains controversial. Besides, the molecular process of TAT peptide facilitating the transport of extracellular substances into cells is largely unknown. In this study, we explore the interactions of TAT peptides and their conjugated gold nanoparticles with lipid membranes by coarse-grained molecular dynamics simulations. It is found that the TAT peptides can hardly penetrate through the membrane at low peptide concentrations; after the concentration increases to a threshold value, they can cross the membrane through an induced nanopore due to the transmembrane electrostatic potential difference. The translocation of TAT peptides is mainly caused by the overall structural changes of membranes. Furthermore, we demonstrate that the translocation of gold nanoparticles (AuNPs) across the membrane is significantly affected by the number of grafted TAT peptides on the particle surface. The transmembrane efficiency of AuNPs may even be reduced when a small number of peptides modify them; whereas, when the number of grafted peptides increases to a certain value, the TAT-AuNP complex can translocate across the membrane in a pore-mediated way. Based on our findings, an effective strategy has been proposed to enhance the delivery efficiency of AuNPs. The present study can improve our understanding of the interactions between TAT peptides and cell membranes; it may also give some insightful suggestions on the design and development of nanocarriers with high efficiency for the delivery of nanoparticles and drugs.

Stretch and Breakage of Wormlike Micelles under Uniaxial Strain: A Simulation Study and Comparison with Experimental Results

We use coarse-grained (CG) molecular dynamics simulations to determine the effect of uniaxial strain on the stress, scission stress, and scission energy of solutions of wormlike micelles of cetyltrimethylammonium chloride/sodium salicylate (NaSal). We find that the breaking stress, stretch modulus, and scission energy of the charged micelles are nonmonotonic functions of oppositely charged hydrotrope (NaSal) concentration. While the stretch modulus shows a peak at a value of surfactant-to-hydrotrope concentration ratio ( R) close to unity as expected due to neutralization of head-group charge at R = 1, the breaking stress and scission energy produce a peak at R < 1.0 because of thinning of the micelle diameter with increased R. The breaking stress from the simulations depends on the rate of deformation and roughly agrees with the experimental values of Rothstein ( J. Rheol. 2003 , 47 , 1227 ) after extrapolation to the much lower experimental rates. The method and results can be used to predict the effects of flow and mechanical stress on rates of micellar breakage, which is important in the rheology of wormlike micellar solutions.

Computational Investigations of the Interaction between the Cell Membrane and Nanoparticles Coated with a Pulmonary Surfactant

When inhaled nanoparticles (NPs) come into the deep lung, they develop a biomolecular corona by interacting with the pulmonary surfactant. The adsorption of the phospholipids and proteins gives a new biological identity to the NPs, which may alter their subsequent interactions with cells and other biological entities. Investigations of the interaction between the cell membrane and NPs coated with such a biomolecular corona are important in understanding the role of the biofluids on cellular uptake and estimating the dosing capacity and the nanotoxicology of NPs. In this paper, using dissipative particle dynamics, we investigate how the physicochemical properties of the coating pulmonary surfactant lipids and proteins affect the membrane response for inhaled NPs. We pinpoint several key factors in the endocytosis of lipid NPs, including the deformation of the coating lipids, coating lipid density, and ligand-receptor binding strength. Further studies reveal that the deformation of the coating lipids consumes energy but on the other hand promotes the coating ligands to bind with receptors more tightly. The coating lipid density controls the amount of the ligands as well as the hydrophobicity of the lipid NPs, thus affecting the endocytosis kinetics through the specific and nonspecific interactions. It is also found that the hydrophobic surfactant proteins associated with lipids can accelerate the endocytosis process of the NPs, but the endocytosis efficiency mainly depends on the density of the coating surfactant lipids. These findings can help understand how the pulmonary surfactant alters the biocompatibility of the inhaled NPs and provide some guidelines in designing an NP complex for efficient pulmonary drug delivery.