Synthetic Polymers and Biomembranes. How Do They Interact?: Atomistic Molecular Dynamics Simulation Study of PEO in Contact with a DMPC Lipid Bilayer

Absorption of Polypropylene in Dipalmitoylphosphatidylcholine Membranes: The Role of Molecular Weight and Initial Configuration of Polymer Chains

We study by molecular dynamics simulations the absorption of polypropylene (PP) chains within a dipalmitoylphosphatidylcholine (DPPC) lipid membrane in aqueous solvent. DPPC represents the most abundant phospholipid in biological membranes, while PP is one of the most common synthetic polymers diffused in the anthropic environment. By following in detail the absorption process, and the corresponding structural modification undergone by the membrane, we show how the initial configuration and the PP molecular weight determine the overall behavior of the system. Specifically, if PP chains initially lie on the DPPC surface, they are fully absorbed; likewise, polymers initially included within the membrane cannot escape from. On the other hand, if polymers are placed sufficiently apart from the membrane, they have time to join together and coalesce into a few nanoparticles. At contact, such nanoparticles may completely dissolve (for low molecular weight) and then be absorbed. For high molecular weight, not all of them dissolve, and therefore the system attains a condition in which some of the chains are absorbed, while others form a residual nanoparticle staying outside (but in contact with) the membrane. Such a state─albeit energetically unfavorable with respect to a condition in which all PP chains are absorbed─remains stable, at the least over a substantial simulation time, extending in our study up to 1.6 μs. The tendency for polymers to spontaneously form aggregates, which then prefer to stay in contact with the membrane, is further corroborated by calculation of the DPPC-nanoparticle potential of mean force.

Location of Polyelectrolytes in Swollen Lipid Oligobilayers

Osteoarthritis is caused by degeneration of the cartilage, which covers the bone ends of the joints and is decorated with an oligolamellar phospholipid (PL) bilayer. The gap between the bone ends is filled with synovial fluid mainly containing hyaluronic acid (HA). HA and PLs are supposed to reduce friction and protect the cartilage from wear in joint movement. However, a detailed understanding of the molecular mechanisms of joint lubrication is still missing. Previously, we found that aqueous solutions of HA and poly(allylamine hydrochloride) (PAH), the latter serving as a polymeric analogue to HA, adsorb onto the headgroups of surface-bound 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) oligobilayers and significantly enhance their stability with respect to shear forces, typically occurring in joint movement. We now investigated the precise location of PAH chains across the lipid films in neutron reflectivity measurements, as bridging of the oligobilayers by polyelectrolytes (PEs) might be the cause for their improved mechanical stability. In a first set of experiments, we used hydrogenated PAH and chain-deuterated DMPC (DMPC-d54) to improve the contrast between the lipids and potentially intruding PAH. However, due to difficulties in distinguishing between incorporation of water and PAH, penetration into the lipid chain region could hardly be proven quantitatively. Therefore, we designed a more elaborate experiment based on mixed films of DMPC-d54 and hydrogenated DMPC, which is insensitive to water penetration into the films. Beside facilitating a detailed structural characterization of the oligolamellar system, this elaborate approach showed that PAH adsorbs to the DMPC heads and penetrates the lipid tail strata. No PAH was found in the lipid head strata, which excludes bridging of several lipid bilayers by the PE chains. The data are consistent with the assumption that PAH bridges are formed between the headgroups of two adjacent bilayers and contribute to the enhanced mechanical stability.

A review on microfluidic-assisted nanoparticle synthesis, and their applications using multiscale simulation methods

Recent years have witnessed an increased interest in the development of nanoparticles (NPs) owing to their potential use in a wide variety of biomedical applications, including drug delivery, imaging agents, gene therapy, and vaccines, where recently, lipid nanoparticle mRNA-based vaccines were developed to prevent SARS-CoV-2 causing COVID-19. NPs typically fall into two broad categories: organic and inorganic. Organic NPs mainly include lipid-based and polymer-based nanoparticles, such as liposomes, solid lipid nanoparticles, polymersomes, dendrimers, and polymer micelles. Gold and silver NPs, iron oxide NPs, quantum dots, and carbon and silica-based nanomaterials make up the bulk of the inorganic NPs. These NPs are prepared using a variety of top-down and bottom-up approaches. Microfluidics provide an attractive synthesis alternative and is advantageous compared to the conventional bulk methods. The microfluidic mixing-based production methods offer better control in achieving the desired size, morphology, shape, size distribution, and surface properties of the synthesized NPs. The technology also exhibits excellent process repeatability, fast handling, less sample usage, and yields greater encapsulation efficiencies. In this article, we provide a comprehensive review of the microfluidic-based passive and active mixing techniques for NP synthesis, and their latest developments. Additionally, a summary of microfluidic devices used for NP production is presented. Nonetheless, despite significant advancements in the experimental procedures, complete details of a nanoparticle-based system cannot be deduced from the experiments alone, and thus, multiscale computer simulations are utilized to perform systematic investigations. The work also details the most common multiscale simulation methods and their advancements in unveiling critical mechanisms involved in nanoparticle synthesis and the interaction of nanoparticles with other entities, especially in biomedical and therapeutic systems. Finally, an analysis is provided on the challenges in microfluidics related to nanoparticle synthesis and applications, and the future perspectives, such as large-scale NP synthesis, and hybrid formulations and devices.

Multiscale Dynamics of Lipid Vesicles in Polymeric Microenvironment

Understanding dynamic and complex interaction of biological membranes with extracellular matrices plays a crucial role in controlling a variety of cell behavior and functions, from cell adhesion and growth to signaling and differentiation. Tremendous interest in tissue engineering has made it possible to design polymeric scaffolds mimicking the topology and mechanical properties of the native extracellular microenvironment; however, a fundamental question remains unanswered: that is, how the viscoelastic extracellular environment modifies the hierarchical dynamics of lipid membranes. In this work, we used aqueous solutions of poly(ethylene glycol) (PEG) with different molecular weights to mimic the viscous medium of cells and nearly monodisperse unilamellar DMPC/DMPG liposomes as a membrane model. Using small-angle X-ray scattering (SAXS), dynamic light scattering, temperature-modulated differential scanning calorimetry, bulk rheology, and fluorescence lifetime spectroscopy, we investigated the structural phase map and multiscale dynamics of the liposome–polymer mixtures. The results suggest an unprecedented dynamic coupling between polymer chains and phospholipid bilayers at different length/time scales. The microviscosity of the lipid bilayers is directly influenced by the relaxation of the whole chain, resulting in accelerated dynamics of lipids within the bilayers in the case of short chains compared to the polymer-free liposome case. At the macroscopic level, the gel-to-fluid transition of the bilayers results in a remarkable thermal-stiffening behavior of polymer–liposome solutions that can be modified by the concentration of the liposomes and the polymer chain length.

Molecular Simulations of PEGylated Biomolecules, Liposomes, and Nanoparticles for Drug Delivery Applications

Since the first polyethylene glycol (PEG)ylated protein was approved by the FDA in 1990, PEGylation has been successfully applied to develop drug delivery systems through experiments, but these experimental results are not always easy to interpret at the atomic level because of the limited resolution of experimental techniques. To determine the optimal size, structure, and density of PEG for drug delivery, the structure and dynamics of PEGylated drug carriers need to be understood close to the atomic scale, as can be done using molecular dynamics simulations, assuming that these simulations can be validated by successful comparisons to experiments. Starting with the development of all-atom and coarse-grained PEG models in 1990s, PEGylated drug carriers have been widely simulated. In particular, recent advances in computer performance and simulation methodologies have allowed for molecular simulations of large complexes of PEGylated drug carriers interacting with other molecules such as anticancer drugs, plasma proteins, membranes, and receptors, which makes it possible to interpret experimental observations at a nearly atomistic resolution, as well as help in the rational design of drug delivery systems for applications in nanomedicine. Here, simulation studies on the following PEGylated drug topics will be reviewed: proteins and peptides, liposomes, and nanoparticles such as dendrimers and carbon nanotubes.

A Multiparametric Fluorescence Probe to Understand the Physicochemical Properties of Small Unilamellar Lipid Vesicles in Poly(ethylene glycol)-Water Medium

FDAPT (2-formyl-5-(4'-N,N-dimethylaminophenyl)thiophene) efficiently senses the minimum alteration of lipid bilayer microenvironment with all six different fluorescence parameters namely emission wavelength, fluorescence intensity, steady-state anisotropy, and their corresponding time-dependent parameters (Sahu et al., J. Phys. Chem. B 2018, 122, 7308-7318). In the present work, the effect of poly(ethylene glycol) on the small unilamellar vesicle is demonstrated with the emission behavior of the FDAPT probe. A medium and a high molecular weight PEG were chosen to perturb the lipid vesicles. The alteration of the bilayer polarity, water content inside bilayer, lipid packing density in the perturbed vesicles reflect significant changes in different fluorescence parameters of FDAPT probe. The effect of PEG on the unilamellar vesicle was rationalized with the alteration of the emission behavior, fluorescence lifetime, steady-state anisotropy and anisotropy decay of the probe. The simple and convenient fluorescence measurements provide new insights into the effect of PEG on the packing density, water volume, micro polarity, and microviscosity of the small unilamellar vesicle. The physiological understanding was extended to rationalize the cryoprotecting behavior of PEG.

Electrospun poly(3-hydroxybutyrate-co-3-hydroxyvalerate)/polyethylene oxide (PEO) microfibers reinforced with ZnO nanocrystals for antibacterial and antibiofilm wound dressing applications

Biocompatible and biodegradable polymers for designing wound dressing materials.

Structural Modifications of DPPC Bilayers upon Inclusion of an Antibacterial Cationic Bolaamphiphile

The emergence of antibiotic-resistant bacterial strains has fostered fundamental research to develop alternative antimicrobial strategies. Among the several systems proposed so far, the association complexes (nanoplexes) formed by transcription factor decoys (TFDs), i.e., short oligonucleotides targeting a crucial bacterial transcription factor, and a bolaform cationic amphiphile, 10,10'-(dodecane-1,12-diyl)-bis-(9-amino-1,2,3,4-tetrahydroacridinium) chloride (12-bis-THA), have demonstrated their potential in vitro and in vivo. The application of these nanoplexes is hampered by a scarce colloidal stability, which can be addressed by including the bolaamphiphile in a liposomal carrier, which is then associated to the TFD. The present study reports an investigation on the effects of 12-bis-THA on the structure of synthetic lipid bilayers to assess the morphology of the mixed assemblies, gain insight into the location of the host within the bilayer, and determine the loading capacity of the carrier. Our results demonstrate that 12-bis-THA promptly inserts within 1,2-dipalmitoyl- sn-glycero-3-phosphocholine (DPPC) bilayers, bending its C-12 spacer chain to adopt a conelike shape and shifting the gel-liquid crystalline transition of the chains to lower temperatures. The host liposomal structure is retained for a bolaamphiphile concentration of up to 3.2% mol to DPPC, whereas higher concentrations lead to the destabilization by means of a detergency-like mechanism, with the simultaneous existence of different lamellar-based structures, such as liposomes, bicelles, and rafts, in which DPPC and 12-bis-THA could be present in different molar ratios. Overall, these results shed light on the interaction of the bolaamphiphile with a lipid bilayer and provide valuable insight to better formulate the antimicrobial amphiphile in liposomal carriers to circumvent the colloidal instability of nanoplexes.

Molecular Dynamics Simulations of a Poly(ethylene glycol)-Grafted Polyamide Membrane and Its Interaction with a Calcium Alginate Gel

Molecular dynamics simulations are carried out to investigate the antifouling property of a polyethylene glycol (PEG)-grafted polyamide (PA) membrane. Our specific interest is the computational study of the interaction between a grafted PEG coating and an alginate gel foulant by a steered molecular dynamics approach. Simulation results show that the PEG coating can hold a tightly bound hydration water layer. When the alginate gel is dragged to approach the PEG coating surface, a strong repulsive hydration force is observed due to the compression of this hydration layer. Detailed calculations of the potential of mean force (PMF) show that the repulsive interaction between the alginate gel and the hydration water layer around the PEG coating has a dominant contribution to the total repulsive PMF. We have also studied the effect of the PEG coverage on the membrane-foulant interactions. We find that the alginate gel has a strong tendency to drift to the uncovered PA membrane surfaces (namely, the PEG hollows). However, direct attachment of the gel to the PA membrane surface can be avoided if the gel size is slightly larger than the PEG hollow site.

Modeling the effect of nano-sized polymer particles on the properties of lipid membranes

The interaction between polymers and biological membranes has recently gained significant interest in several research areas. On the biomedical side, dendrimers, linear polyelectrolytes, and neutral copolymers find application as drug and gene delivery agents, as biocidal agents, and as platforms for biological sensors. On the environmental side, plastic debris is often disposed of in the oceans and gets degraded into small particles; therefore concern is raising about the interaction of small plastic particles with living organisms. From both perspectives, it is crucial to understand the processes driving the interaction between polymers and cell membranes. In recent times progress in computer technology and simulation methods has allowed computational predictions on the molecular mechanism of interaction between polymeric materials and lipid membranes. Here we review the computational studies on the interaction between lipid membranes and different classes of polymers: dendrimers, linear charged polymers, polyethylene glycol (PEG) and its derivatives, polystyrene, and some generic models of polymer chains. We conclude by discussing some of the technical challenges in this area and future developments.

Coarse-grained models for aqueous polyethylene glycol solutions

A new coarse-grained force field is developed for polyethylene glycol (PEG) in water. The force field is based on the MARTINI model but with the big multipole water (BMW) model for the solvent. The polymer force field is reparameterized using the MARTINI protocol. The new force field removes the ring-like conformations seen in simulations of short chains with the MARTINI force field; these conformations are not observed in atomistic simulations. We also investigate the effect of using parameters for the end-group that are different from those for the repeat units, with the MARTINI and BMW/MARTINI models. We find that the new BMW/MARTINI force field removes the ring-like conformations seen in the MARTINI models and has more accurate predictions for the density of neat PEG. However, solvent-separated-pairs between chain ends and slow dynamics of the PEG reflect its own artifacts. We also carry out fine-grained simulations of PEG with bundled water clusters and show that the water bundling can lead to ring-like conformations of the polymer molecules. The simulations emphasize the pitfalls of coarse-graining several molecules into one site and suggest that polymer-solvent systems might be a stringent test for coarse-grained force fields.

Cell surface engineering with polyelectrolyte multilayer thin films

Layer-by-layer assembly of polyelectrolyte multilayer (PEM) films represents a bottom-up approach for re-engineering the molecular landscape of cell surfaces with spatially continuous and molecularly uniform ultrathin films. However, fabricating PEMs on viable cells has proven challenging owing to the high cytotoxicity of polycations. Here, we report the rational engineering of a new class of PEMs with modular biological functionality and tunable physicochemical properties which have been engineered to abrogate cytotoxicity. Specifically, we have discovered a subset of cationic copolymers that undergoes a conformational change, which mitigates membrane disruption and facilitates the deposition of PEMs on cell surfaces that are tailorable in composition, reactivity, thickness, and mechanical properties. Furthermore, we demonstrate the first successful in vivo application of PEM-engineered cells, which maintained viability and function upon transplantation and were used as carriers for in vivo delivery of PEMs containing biomolecular payloads. This new class of polymeric film and the design strategies developed herein establish an enabling technology for cell transplantation and other therapies based on engineered cells.

Molecular dynamics simulation of interactions between a sodium dodecyl sulfate micelle and a poly(ethylene oxide) polymer

We have performed atomistic molecular dynamics simulations of an anionic sodium dodecyl sulfate (SDS) micelle and a nonionic poly(ethylene oxide) (PEO) polymer in aqueous solution. The micelle consisted of 60 surfactant molecules, and the polymer chain lengths varied from 20 to 40 monomers. The force field parameters for PEO were adjusted by using 1,2-dimethoxymethane (DME) as a model compound and matching its hydration enthalpy and conformational behavior to experiment. Excellent agreement with previous experimental and simulation work was obtained through these modifications. The simulated scaling behavior of the PEO radius of gyration was also in close agreement with experimental results. The SDS-PEO simulations show that the polymer resides on the micelle surface and at the hydrocarbon-water interface, leading to a selective reduction in the hydrophobic contribution to the solvent-accessible surface area of the micelle. The association is mainly driven by hydrophobic interactions between the polymer and surfactant tails, while the interaction between the polymer and sulfate headgroups on the micelle surface is weak. The 40-monomer chain is mostly wrapped around the micelle, and nearly 90% of the monomers are adsorbed at low PEO concentration. Simulations were also performed with multiple 20-monomer chains, and gradual addition of polymer indicates that about 120 monomers are required to saturate the micelle surface. The stoichiometry of the resulting complex is in close agreement with experimental results, and the commonly accepted "beaded necklace" structure of the SDS-PEO complex is recovered by our simulations.