Modeling the self-assembly of lipids and nanotubes in solution: forming vesicles and bicelles with transmembrane nanotube channels

Nanoporous Membranes of Densely Packed Carbon Nanotubes Formed by Lipid-Mediated Self-Assembly

Nanofiltration technology faces the competing challenges of achieving high fluid flux through uniformly narrow pores of a mechanically and chemically stable filter. Supported dense-packed 2D-crystals of single-walled carbon nanotube (CNT) porins with ∼1 nm wide pores could, in principle, meet these challenges. However, such CNT membranes cannot currently be synthesized at high pore density. Here, we use computer simulations to explore lipid-mediated self-assembly as a route toward densely packed CNT membranes, motivated by the analogy to membrane-protein 2D crystallization. In large-scale coarse-grained molecular dynamics (MD) simulations, we find that CNTs in lipid membranes readily self-assemble into large clusters. Lipids trapped between the CNTs lubricate CNT repacking upon collisions of diffusing clusters, thereby facilitating the formation of large ordered structures. Cluster diffusion follows the Saffman-Delbrück law and its generalization by Hughes, Pailthorpe, and White. On longer time scales, we expect the formation of close-packed CNT structures by depletion of the intervening shared annular lipid shell, depending on the relative strength of CNT-CNT and CNT-lipid interactions. Our simulations identify CNT length, diameter, and end functionalization as major factors for the self-assembly of CNT membranes.

Mesoscale Modeling of Agglomeration of Molecular Bottlebrushes: Focus on Conformations and Clustering Criteria

Using dissipative particle dynamics, we characterize dynamics of aggregation of molecular bottlebrushes in solvents of various qualities by tracking the number of clusters, the size of the largest cluster, and an average aggregation number. We focus on a low volume fraction of bottlebrushes in a range of solvents and probe three different cutoff criteria to identify bottlebrushes belonging to the same cluster. We demonstrate that the cutoff criteria which depend on both the coordination number and the length of the side chain allows one to correlate the agglomeration status with the structural characteristics of bottlebrushes in solvents of various qualities. We characterize conformational changes of the bottlebrush within the agglomerates with respect to those of an isolated bottlebrush in the same solvents. The characterization of bottlebrush conformations within the agglomerates is an important step in understanding the relationship between the bottlebrush architecture and material properties. An analysis of three distinct cutoff criteria to identify bottlebrushes belonging to the same cluster introduces a framework to identify both short-lived transient and long-lived agglomerates; the same approach could be further extended to characterize agglomerates of various macromolecules with complex architectures beyond the specific bottlebrush architecture considered herein.

Controlling Degradation and Erosion of Polymer Networks: Insights from Mesoscale Modeling

Understanding and controlling degradation of polymer networks on the mesoscale is critical for a range of applications. We utilize dissipative particle dynamics to capture photocontrolled degradation and erosion processes in hydrogels formed by end-linking of four-arm polyethylene glycol precursors. We demonstrate that the polydispersity and the fraction of broken-off fragments scale with the relative extent of reaction. The reverse gel point measured is close to the value predicted by the bond percolation theory on a diamond lattice. We characterize the erosion process via tracking the mass loss that accounts for the fragments remaining in contact with the percolated network. We quantify the dependence of the mass loss on the extent of reaction and on the properties of the film prior to degradation. These results elucidate the main features of degradation and erosion on the mesoscale and could provide guidelines for future design of degrading materials with dynamically controlled properties.

A hybrid coarse-grained model for structure, solvation and assembly of lipid-like peptides

Reconstituted photosynthetic proteins which are activated upon exposure to solar energy hold enormous potential for powering future solid state devices and solar cells. The functionality and integration of these proteins into such devices has been successfully enabled by lipid-like peptides. Yet, a fundamental understanding of the organization of these peptides with respect to the photosynthetic proteins and themselves remains unknown and is critical for guiding the design of such light-activated devices. This study investigates the relative organization of one such peptide sequence V₆K₂ (V: valine and K: lysine) within assemblies. Given the expansive spatiotemporal scales associated with this study, a hybrid coarse-grained (CG) model which captures the structure, conformation and aggregation of the peptide is adopted. The CG model uses a combination of iterative Boltzmann inversion and force matching to provide insight into the relative organization of V₆K₂ in assemblies. The CG model reproduces the structure of a V₆K₂ peptide sequence along with its all atom (AA) solvation structure. The relative organization of multiple peptides in an assembly, as captured by CG simulations, is in agreement with corresponding results from AA simulations. Also, a backmapping procedure reintroduces the AA details of the peptides within the aggregates captured by the CG model to demonstrate the relative organization of the peptides. Furthermore, a large number of peptides self-assemble into an elongated micelle in the CG simulation, which is consistent with experimental findings. The coarse-graining procedure is tested for transferability to longer peptide sequences, and hence can be extended to other amphiphilic peptide sequences.

High Adsorption of α-Glucosidase on Polymer Brush-Modified Anisotropic Particles Acquired by Electrospraying-A Combined Experimental and Simulation Study

In this particular contribution, we aim to immobilize a model enzyme such as α-glucosidase onto poly(DMAEMA) [poly(2-dimethyl amino ethyl methacrylate)] brush-modified anisotropic (cup- and disc-shaped) biocompatible polymeric particles. The anisotropic particles comprising a blend of PLA [poly(lactide)] and poly(MMA-co-BEMA) [poly((methyl methacrylate)-co-(2-(2-bromopropionyloxy) ethyl methacrylate)] were acquired by electrospraying, a scalable and convenient technique. We have also demonstrated the role of a swollen polymer brush grafted on the surface of cup-/disc-shaped particles via surface-initiated atom transfer radical polymerization in immobilizing an unprecedentedly high loading of enzyme [441 mg/g (cup)-589 mg/g (disc) of particles, 15-20 times higher than that of the literature-reported system] as compared to non-brush-modified particles. Circular dichroism spectroscopy was used to predict the structural changes of the enzyme upon immobilization onto the carrier particles. An enormously high amount of enzymes with preserved activity (∼85 ± 13% for cups and ∼78 ± 15% for discs) was found to adhere onto brush-modified particles at pH 7 via electrostatic adsorption. These findings were further explored at the atomistic level using a coarse-grained dissipative particle dynamics simulation approach, which exhibited excellent correlation with experimental results. In addition, accelerated particle separation was also achieved via magnetic force-induced aggregation within 20 min (without a centrifuge) by incorporating magnetic nanoparticles into disc-shaped particles while electrojetting. This further strengthens the technical feasibility of the process, which holds immense potential to be applied for various enzymes intended for several applications.

Native-Based Dissipative Particle Dynamics Approach for α-Helical Folding

We developed a dissipative particle dynamics (DPD) approach that captures polyalanine folding into a stable helical conformation. Within the proposed native-based approach, the DPD parameters are derived based on the contact map constructed from the molecular dynamics (MD) simulations. We show that the proposed approach reproduces the folding of polypeptides of various lengths, including bundle formation for sufficiently long polypeptides. The proposed approach also allows one to capture the folding of the helical segments of the lysozyme. With further development of computationally efficient native-based DPD approaches for folding, modeling of a range of biomaterials incorporating α-helical segments could be extended to time and length scales far beyond those accessible in molecular dynamics simulations.

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.

Penetration of nanoparticles across a lipid bilayer: effects of particle stiffness and surface hydrophobicity

The cellular uptake of nanoparticles (NPs) has drawn significant attention due to their great importance and potential in drug delivery, bioimaging, and specific targeting. Here, we conduct a computational study on the translocation process of soft nanoparticles with different elasticities and surface hydrophobicities through a lipid bilayer membrane. It is shown that the translocation abilities of hydrophilic NPs can be enhanced by increasing their stiffness, while the penetrability of hydrophobic NPs is weakened by increasing the particle stiffness. The free energy analysis indicates that rigid hydrophilic NPs and soft hydrophobic NPs encounter lower energy barriers during penetration. In direct translocation, different deformation modes are observed for NPs with different surface hydrophobicities during cellular internalization. Further, deformation analysis demonstrates that hydrophilic NPs are flattened in the membrane plane, while hydrophobic NPs are elongated along the membrane norm during penetration. We conclude that the elasticity of NPs has an obvious impact on their ability to penetrate across the lipid bilayer membrane through different morphological responses of hydrophilic and hydrophobic NPs. These results shed light on the coupled effects of particle elasticity and surface hydrophobicity on the cellular uptake of elastic NPs, which may provide useful guidelines for designing effective nanocarrier systems for drug delivery.

Indentation of Graphene-Covered Atomic Force Microscopy Probe Across a Lipid Bilayer Membrane: Effect of Tip Shape, Size, and Surface Hydrophobicity

Understanding the interaction of graphene with cell membranes is crucial to the development of graphene-based biological applications and the management of graphene safety issues. To help reveal the key factors controlling the interaction between graphene and cell membranes, here we adopt the dissipative particle dynamics method to analyze the evolution of interaction force and free energy as the graphene-covered atomic force microscopy (AFM) probe indents across a lipid bilayer. The simulation results show that the graphene-covered AFM probe can cause severe deformation of the cell membrane which drives the lipid molecule to adsorb and diffuse at the surface of graphene. The breakthrough force and free energy are calculated to study the effects of the tip shape, size, and surface hydrophobicity on the piercing behaviors of graphene-covered AFM. In addition, the deformation of cell membrane can decrease the dependency of the breakthrough force on the tip shape. The analysis of surface functionalization suggests that the horizontal patterns on graphene can change the preferred orientation in the penetration process, but the vertical patterns on graphene may disrupt the cell membrane. What's more, the bending stiffness of graphene has little influence on the penetration process as graphene pierces into the cell membrane. These results provide useful guidelines for the molecular design of graphene materials with controllable cell penetrability.

Molecular dynamics simulations of carbon nanotube porins in lipid bilayers

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

Perturbation of the pulmonary surfactant monolayer by single-walled carbon nanotubes: a molecular dynamics study

Single-walled carbon nanotubes (SWCNTs) are at present synthesized on a large scale with a variety of applications. The increasing likelihood of exposure to SWCNTs, however, puts human health at a high risk. As the front line of the innate host defense system, the pulmonary surfactant monolayer (PSM) at the air-water interface of the lungs interacts with the inhaled SWCNTs, which in turn inevitably perturb the ultrastructure of the PSM and affect its biophysical functions. Here, using molecular dynamics simulations, we demonstrate how the diameter and length of SWCNTs critically regulate their interactions with the PSM. Compared to their diameters, the inhalation toxicity of SWCNTs was found to be largely affected by their lengths. Short SWCNTs with lengths comparable to the monolayer thickness are found to vertically insert into the PSM with no indication of translocation, possibly leading to accumulation of SWCNTs in the PSM with prolonged retention and increased inflammation potentials. The perturbation also comes from the forming water pores across the PSM. Longer SWCNTs are found to horizontally insert into the PSM during inspiration, and they can be wrapped by the PSM during deep expiration via a tube diameter-dependent self-rotation. The potential toxicity of longer SWCNTs comes from severe lipid depletion and the PSM-rigidifying effect. Our findings could help reveal the inhalation toxicity of SWCNTs, and pave the way for the safe use of SWCNTs as vehicles for pulmonary drug delivery.

Flow-Induced Shape Reconfiguration, Phase Separation, and Rupture of Bio-Inspired Vesicles

The structural integrity of red blood cells and drug delivery carriers through blood vessels is dependent upon their ability to adapt their shape during their transportation. Our goal is to examine the role of the composition of bio-inspired multicomponent and hairy vesicles on their shape during their transport through in a channel. Through the dissipative particle dynamics simulation technique, we apply Poiseuille flow in a cylindrical channel. We investigate the effect of flow conditions and concentration of key molecular components on the shape, phase separation, and structural integrity of the bio-inspired multicomponent and hairy vesicles. Our results show the Reynolds number and molecular composition of the vesicles impact their flow-induced deformation, phase separation on the outer monolayer due to the Marangoni effect, and rupture. The findings from this study could be used to enhance the design of drug delivery and tissue engineering systems.

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.

Modeling Interactions between Multicomponent Vesicles and Antimicrobial Peptide-Inspired Nanoparticles

We examine the interaction between peptide-inspired nanoparticles, or nanopins, and multicomponent vesicles using the dissipative particle dynamics simulation technique. We study the role of nanopin architecture and cholesterol concentration on the binding of the nanopins to the lipid bilayer, their insertion, and postembedding self-organization. We find the insertion to be triggered by enthalpically unfavorable interactions between the hydrophilic solvent and the lipophilic components of the nanopins. The nanopins are observed to form aggregates in solution, insert into the bilayer, and disassemble into the individual nanopins following the insertion process. We examine factors that influence the orientation of the nanopins in the host vesicle. We report the length of the hydrophilic segment of the nanopins to regulate their orientation within the clusters before the embedding process and in the bilayer, after the postinsertion disassembly of the aggregates. The orientation angle distribution for a given nanopin architecture is found to be driven by energy minimization. In addition, higher concentration of cholesterol is observed to constrain the orientation of the nanopins. We also report thermal fluctuations to induce transverse diffusion of nanopins with specific architectures. The incidence of transverse diffusion is observed to decrease with the concentration of cholesterol. Our results can provide guidelines for designing peptide-inspired nanoparticles or macromolecules that can interface with living cells to serve as sensors for applications in medicine, sustainability, and energy.

Self-Assembly and Critical Aggregation Concentration Measurements of ABA Triblock Copolymers with Varying B Block Types: Model Development, Prediction, and Validation

The dissipative particle dynamics (DPD) simulation technique is a coarse-grained (CG) molecular dynamics-based approach that can effectively capture the hydrodynamics of complex systems while retaining essential information about the structural properties of the molecular species. An advantageous feature of DPD is that it utilizes soft repulsive interactions between the beads, which are CG representation of groups of atoms or molecules. In this study, we used the DPD simulation technique to study the aggregation characteristics of ABA triblock copolymers in aqueous medium. Pluronic polymers (PEG-PPO-PEG) were modeled as two segments of hydrophilic beads and one segment of hydrophobic beads. Tyrosine-derived PEG5K-b-oligo(desaminotyrosyl tyrosine octyl ester-suberate)-b-PEG5K (PEG5K-oligo(DTO-SA)-PEG5K) block copolymers possess alternate rigid and flexible components along the hydrophobic oligo(DTO-SA) chain, and were modeled as two segments of hydrophilic beads and one segment of hydrophobic, alternate soft and hard beads. The formation, structure, and morphology of the initial aggregation of the polymer molecules in aqueous medium were investigated by following the aggregation dynamics. The dimensions of the aggregates predicted by the computational approach were in good agreement with corresponding results from experiments, for the Pluronic and PEG5K-oligo(DTO-SA)-PEG5K block copolymers. In addition, DPD simulations were utilized to determine the critical aggregation concentration (CAC), which was compared with corresponding results from an experimental approach. For Pluronic polymers F68, F88, F108, and F127, the computational results agreed well with experimental measurements of the CAC measurements. For PEG5K-b-oligo(DTO-SA)-b-PEG5K block polymers, the complexity in polymer structure made it difficult to directly determine their CAC values via the CG scheme. Therefore, we determined CAC values of a series of triblock copolymers with 3-8 DTO-SA units using DPD simulations, and used these results to predict the CAC values of triblock copolymers with higher molecular weights by extrapolation. In parallel, a PEG5K-b-oligo(DTO-SA)-b-PEG5K block copolymer was synthesized, and the CAC value was determined experimentally using the pyrene method. The experimental CAC value agreed well with the CAC value predicted by simulation. These results validate our CG models, and demonstrate an avenue to simulate and predict aggregation characteristics of ABA amphiphilic triblock copolymers with complex structures.

Prediction of solvent-induced morphological changes of polyelectrolyte diblock copolymer micelles

Self-assembly processes of polyelectrolyte block copolymers are ubiquitous in industrial and biological processes; understanding their physical properties can also provide insights into the design of polyelectrolyte materials with novel and tailored properties. Here, we report systematic analysis on how the ionic strength of the solvent and the length of the polyelectrolyte block affect the self-assembly and morphology of the polyelectrolyte block copolymer materials by constructing a salt-dependent morphological phase diagram using an implicit solvent ionic strength (ISIS) method for dissipative particle dynamics (DPD) simulations. This diagram permits the determination of the conditions for the morphological transition into a specific shape, namely vesicles or lamellar aggregates, wormlike/cylindrical micelles, and spherical micelles. The scaling behavior for the size of spherical micelles is predicted, in terms of radius of gyration (R(g,m)) and thickness of corona (Hcorona), as a function of solvent ionic strength (c(s)) and polyelectrolyte length (NA), which are R(g,m) ∼ c(s)(-0.06)N(A)(0.54) and Hcorona ∼ c(s)(-0.11)N(A)(0.75). The simulation results were corroborated through AFM and static light scattering measurements on the example of the self-assembly of monodisperse, single-stranded DNA block-copolynucleotides (polyT50-b-F-dUTP). Overall, we were able to predict the salt-responsive morphology of polyelectrolyte materials in aqueous solution and show that a spherical-cylindrical-lamellar change in morphology can be obtained through an increase in solvent ionic strength or a decrease of polyelectrolyte length.

Molecular modeling of membrane responses to the adsorption of rotating nanoparticles: promoted cell uptake and mechanical membrane rupture

Recently, a unique dynamic magnetic field was developed to induce the rotational movement of superparamagnetic iron oxide nanoparticles. This technique has been applied to remotely control both cellular internalization and apoptosis. Therefore, a thorough understanding of how a lipid membrane responds to the introduction of rotating NPs is quite important to promote the applications of this technique in a variety of biomedical area. Here, we performed Dissipative Particle Dynamics (DPD) simulations to systematically investigate the interaction mechanism between lipid membranes and rotating NPs. Two kinds of membrane responses are observed. One is the promoted cell uptake and the other is the mechanical membrane rupture. The promoting effect of NP rotation on the cell uptake is ascribed to the enhanced membrane monolayer protrusion, which can wrap the NP from the top side. Meanwhile, the rotating NP exerts a shearing force on the membrane. Accordingly, the membrane undergoes a local distortion around the NP. If the shearing force exceeds a critical value, the local membrane distortion develops into a mechanical rupture. A number of factors, like NP size, NP shape, ligand density and rotation speed, are critical in both of the above membrane responses.

Molecular modeling of interaction between lipid monolayer and graphene nanosheets: implications for pulmonary nanotoxicity and pulmonary drug delivery

Understanding how nanoparticles interact with the pulmonary surfactant monolayer (PSM) is of great importance for safe applications in biomedicine and for evaluation of both health and environment impacts.

Hierarchical multi-scale simulations of adhesion at polymer–metal interfaces: dry and wet conditions

Multi-scale simulations are performed to study the adhesion properties of different polymer–metal interfaces in the absence and presence of water.

Synthetic biology outside the cell: linking computational tools to cell-free systems

As mathematical models become more commonly integrated into the study of biology, a common language for describing biological processes is manifesting. Many tools have emerged for the simulation of in vivo synthetic biological systems, with only a few examples of prominent work done on predicting the dynamics of cell-free synthetic systems. At the same time, experimental biologists have begun to study dynamics of in vitro systems encapsulated by amphiphilic molecules, opening the door for the development of a new generation of biomimetic systems. In this review, we explore both in vivo and in vitro models of biochemical networks with a special focus on tools that could be applied to the construction of cell-free expression systems. We believe that quantitative studies of complex cellular mechanisms and pathways in synthetic systems can yield important insights into what makes cells different from conventional chemical systems.

Phase segregation in bio-inspired multi-component vesicles encompassing double tail phospholipid species

Our aim is to investigate the phase segregation and the structure of multi-component bio-inspired phospholipid vesicles via dissipative particle dynamics. The chemical distinction in the phospholipid species arises due to different head and tail group moieties, and molecular stiffness of the hydrocarbon tails. The individual amphiphilic phospholipid molecular species are represented by a hydrophilic head group and two hydrophobic tails. The distinct chemical nature of the moieties is modeled effectively via soft repulsive interaction parameters, and the molecular rigidity is tuned via suitable three-body potential constants. We demonstrate the formation of a stable hybrid vesicle through the self-assembly of the amphiphilic phospholipid molecules in the presence of a hydrophilic solvent. We investigate and characterize the phase segregation and the structure of the binary vesicles for different phospholipid mixtures. Our results demonstrate macroscopic phase separation for phospholipid mixtures composed of species with different hydrocarbon tail groups. We also investigate the relationship between the phase segregation and thermodynamic variables such as interfacial line tension and surface tension, and obtain correspondence between existing theory and experiments, and our simulation results. We report variations in the molecular chain stiffness to have negligible contributions to the phase segregation in the mixed bilayer, and to demonstrate shape transformations of the hybrid vesicle. Our results can be used to design novel bio-inspired hybrid vehicles for potential applications in biomedicine, sensing, imaging and sustainability.

Self-healing vesicles deposit lipid-coated Janus particles into nanoscopic trenches

Using dissipative particle dynamics (DPD) simulations, we model the interaction between nanoscopic lipid vesicles and Janus nanoparticles localized on an adhesive substrate in the presence of an imposed flow. The system is immersed in a hydrophilic solution, and the hydrophilic substrate contains nanoscopic trenches, which are either step- or wedge-shaped. The fluid-driven vesicle successfully picks up Janus particles on the substrate, transports these particles as cargo along the surface, and then drops off the particles into the trenches. For Janus particles with a relatively large hydrophobic region, lipids from the bilayer membrane become detached from the vesicle and bound to the hydrophobic domain of the deposited particle. While the detachment of these lipids rips the vesicle, it provides a coating that effectively shields the hydrophobic portion of the nanoparticle from the outer solution. After the particle has been dropped off, the torn vesicle undergoes structural rearrangement, reforming into a closed structure that resembles its original shape. In effect, the vesicle displays pronounced adaptive behavior, shedding lipids to form a protective coating around the particle and undergoing a self-healing process after the particle has been deposited. This responsive, adaptive behavior is observed in cases involving both the step- and wedge-shaped trenches, but the step trench is more effective at inducing particle drop off. The results reveal that the introduction of grooves or trenches into a hydrophilic surface can facilitate the targeted delivery of amphiphilic particles by self-healing vesicles, which could be used for successive delivery events.

Nano-pipette directed transport of nanotube transmembrane channels and hybrid vesicles

Using computational modeling, we simulate the interactions between a nanopipette and transmembrane, end-functionalized nanotubes that are localized within flat bilayers or nanoscopic vesicles. The functional groups (hairs) provide a "handle" for the moving pipette to controllably pick up and move the nanotubes to specific locations in the flat membrane, or the hybrid vesicle to specified regions on a surface. The ability to localize these hybrid vesicles on surfaces paves the way for creating nanoreactor arrays in fluidic devices.

"Zero-dimensional" single-walled carbon nanotubes

The shorter, the more dispersible: An iterative, emulsion-based shortening technique has been used to reduce the length of single-walled carbon nanotubes (SWNTs) to the same order of magnitude as their diameter (ca. 1 nm), thus achieving an effectively "zero-dimensional" structure with improved dispersibility and, after hydroxylation, long-term water solubility. Finally, zero-dimensional SWNTs were positively identified using mass spectrometry for the first time.

Novel Ultrathin Membranes Composed of Organic Ions

Until recently, construction of bilayers was an exclusive mission of nature. It requires careful choice of compounds, whose delicate interplay between head group attraction and chain repulsion engenders a truly unique balance over a narrow temperature range. We report the investigation of artificial bilayers composed of long-chained organic ions, such as dodecyltrimethylammonium (DMA(+)) and perfluorooctaonate (PFO(-)). Various ratios of DMA/PFO surfactants result in bilayers of different stability, thickness, area per molecule, and density profiles. In our quest for water filtration, we incorporated aquaporin protein into the DMA/PFO bilayer but did not observe sufficient stability of the system. We discuss further steps to utilize these surfactant bilayers as highly selective, salt-impermeable membranes.

Harnessing fluid-driven vesicles to pick up and drop off Janus particles

Using dissipative particle dynamics (DPD) simulations, we model the interaction between nanoscopic lipid vesicles and Janus nanoparticles in the presence of an imposed flow. Both the vesicle and Janus nanoparticles are localized on a hydrophilic substrate and immersed in a hydrophilic solution. The fluid-driven vesicle successfully picks up Janus particles on the substrate and transports these particles as cargo along the surface. The vesicle can carry up to four particles as its payload. Hence, the vesicles can act as nanoscopic "vacuum cleaners", collecting nanoscopic debris localized on the floors of the fluidic devices. Importantly, these studies reveal how an imposed flow can facilitate the incorporation of nanoparticles into nanoscale vesicles. With the introduction of a "sticky" domain on the substrate, the vesicles can also robustly drop off and deposit the particles on the surface. The controlled pickup and delivery of nanoparticles via lipid vesicles can play an important step in the bottom-up assembly of these nanoparticles within small-scale fluidic devices.

DPD Simulation of Protein Conformations: From α-Helices to β-Structures

We suggest a coarse-grained model for DPD simulations of polypeptides in solutions. The model mimics hydrogen bonding that stabilizes α-helical and β-structures using dissociable Morse bonds between quasiparticles representing the peptide groups amenable to hydrogen bonding. We demonstrate the capabilities of the model by simulating transitions between coil-like, globular, α-helical, and β-hairpin configurations of model peptides, varying Morse potential parameters, the hydrophobicities of residue side chains, and pH, which determines the charges of residue side chains. We construct a model triblock polypeptide mimicking the sequence of residues α-synuclein at two different pHs. The conformations of this model polypeptide depend on pH similarly to the behavior observed experimentally. The suggested approach to accounting for hydrogen bond formation within the general DPD framework may make the DPD method a competitive alternative to CGMD for modeling equilibrium and dynamic properties of proteins and polypeptides, especially during their transport in confined environments.

Surface-structure-regulated penetration of nanoparticles across a cell membrane

The cell uptake rate of nanoparticles (NPs) coated with mixed hydrophilic/hydrophobic ligands is known to be strongly influenced by the ligand pattern on the nanoparticle surface. To help reveal the physical mechanism behind this intriguing phenomenon, here we perform dissipative particle dynamics simulations to analyze the evolution of free energy as the ligand-coated NPs pierce through a lipid bilayer. Four characteristic ligand patterns are considered: striated NPs with alternating hydrophilic and hydrophobic groups compared to NPs with randomly mixed ligands at the same hydrophilic to hydrophobic ratio, as well as NPs coated with homogeneous hydrophilic or hydrophobic ligands. The free energy analysis indicates that among the four ligand patterns under study, the striated NP encounters the lowest energy barrier during translocation across the membrane. Further analysis reveals that the translocation of the striated NP is facilitated by the constraint of its rotational degree of freedom by the anisotropic ligand pattern, which prevented the free energy of the system from sinking to a deeper valley as the NP passes through the hydrophobic core of the bilayer. Finally, the critical forces required for almost instant penetration of these patterned NPs across the bilayer are calculated and shown to be consistent with the free energy analysis. These findings provide useful guidelines for the molecular design of patterned NPs for controllable cell penetrability.

Nanochannel with uniform and Janus surfaces: shear thinning and thickening in surfactant solution

On basis of molecular simulation of confined surfactant solutions, we show that by adding chemical patterns on the inner surface of nanochannels dynamical properties of the confined surfactant solutions could be modified from shear thinning to shear thickening. To this end, we select uniformly hydrophobic and hydrophilic surfaces as well as a stripe-patterned Janus surface as three prototype confining surfaces of nanochannels. In all three nanochannels, when the surfactant solution is under relatively low shear rates, it shears thin. Under moderate shear rates, a sharp decrease in the shear viscosity could occur due to surfactant morphology transition. Under relatively high shear rates, a shear-thinning-to-thickening transition can emerge due to the tendency of stratification normal to the confining surface. Our simulation study offers a guide to steering dynamic properties of surfactant fluids in nanofluidic devices through engineering surfaces of nanochannels by design.

Controlled release of nanoparticles and macromolecules from responsive microgel capsules

Using a mesoscale computational model, we probe the release of nanoparticles and linear macromolecules from hollow microgel capsules that swell and deswell in response to external stimuli. Our simulations reveal that responsive microcapsules can be effectively utilized for steady and pulsatile release of encapsulated solutes. Swollen gel capsules allow steady, diffusive release of nanoparticles and polymer chains, whereas gel deswelling causes burst-like discharge of solutes driven by an outward flow of the solvent enclosed within a shrinking capsule. We demonstrate that this hydrodynamic release can be regulated by introducing rigid microscopic rods in the capsule interior. Thus, our findings disclose an efficient approach for controlled release from stimuli-responsive microcarriers that could be useful for designing advanced drug delivery systems.