Diblock copolymer bilayers as model for polymersomes: A coarse grain approach

RE, Mocskos EE, Facelli JC, Pickholz M, Ferraro MB. Molecular dynamics study of the mechanical properties of drug loaded model systems: A comparison of a polymersome with a bilayer

In this work we implement a new methodology to study structural and mechanical properties of systems having spherical and planar symmetries throughout Molecular Dynamics simulations. This methodology is applied here to a drug delivery system based in polymersomes, as an example. The chosen model drug was the local anesthetic prilocaine due to previous parameterization within the used coarse grain scheme. In our approach, mass density profiles (MDPs) are used to obtain key structural parameters of the systems, and pressure profiles are used to estimate the curvature elastic parameters. The calculation of pressure profiles and radial MPDs required the development of specific methods, which were implemented in an in-house built version of the GROMACS 2018 code. The methodology presented in this work is applied to characterize poly(ethylene oxide)-poly(butadiene) polymersomes and bilayers loaded with the model drug prilocaine. Our results show that structural properties of the polymersome membrane could be obtained from bilayer simulations, with significantly lower computational cost compared to whole polymersome simulations, but the bilayer simulations are insufficient to get insights on their mechanical aspects, since the elastic parameters are canceled out for the complete bilayer (as consequence of the symmetry). The simulations of entire polymersomes, although more complex, offer a complementary approach to get insights on the mechanical behavior of the systems.

Unraveling the Phase Behavior, Mechanical Stability, and Protein Reconstitution Properties of Polymer-Lipid Hybrid Vesicles

Hybrid vesicles consisting of natural phospholipids and synthetic amphiphilic copolymers have shown remarkable material properties and potential for biotechnology, combining the robustness of polymers with the biocompatibility of phospholipid membranes. To predict and optimize the mixing behavior of lipids and copolymers, as well as understand the interaction between the hybrid membrane and macromolecules like membrane proteins, a comprehensive understanding at the molecular level is essential. This can be achieved by a combination of molecular dynamics simulations and experiments. Here, simulations of POPC and PBD22-b-PEO14 hybrid membranes are shown, uncovering different copolymer configurations depending on the polymer-to-lipid ratio. High polymer concentrations created thicker membranes with an extended polymer conformation, while high lipid content led to the collapse of the polymer chain. High concentrations of polymer were further correlated with a decreased area compression modulus and altered lateral pressure profiles, hypothesized to result in the experimentally observed improvement in membrane protein reconstitution and resistance toward destabilization by detergents. Finally, simulations of a WALP peptide embedded in the bilayer showed that only membranes with up to 50% polymer content favored a transmembrane configuration. These simulations correlate with previous and new experimental results and provide a deeper understanding of the properties of lipid-copolymer hybrid membranes.

Controlling self-assembling co-polymer coatings of hydrophilic polysaccharide substrates via co-polymer block length ratio

HYPOTHESIS

The degree of polymerization of amphiphilic di-block co-polymers, which can be varied with ease in computer simulations, provides a means to control self-assembling di-block co-polymer coatings on hydrophilic substrates.

SIMULATIONS

We examine self-assembly of linear amphiphilic di-block co-polymers on hydrophilic surface via dissipative particle dynamics simulations. The system models a glucose based polysaccharide surface on which random co-polymers of styrene and n-butyl acrylate, as the hydrophobic block, and starch, as the hydrophilic block, forms a film. Such setups are common in e.g. hygiene, pharmaceutical, and paper product applications.

FINDINGS

Variation of the block length ratio (35 monomers in total) reveals that all examined compositions readily coat the substrate. However, strongly asymmetric block co-polymers with short hydrophobic segments are best in wetting the surface, whereas approximately symmetric composition leads to most stable films with highest internal order and well-defined internal stratification. At intermediate asymmetries, isolated hydrophobic domains form. We map the sensitivity and stability of the assembly response for a large variety of interaction parameters. The reported response persists for a wide polymer mixing interactions range, providing general means to tune surface coating films and their internal structure, including compartmentalization.

The Martini Model in Materials Science

The Martini model, a coarse-grained force field initially developed with biomolecular simulations in mind, has found an increasing number of applications in the field of soft materials science. The model's underlying building block principle does not pose restrictions on its application beyond biomolecular systems. Here, the main applications to date of the Martini model in materials science are highlighted, and a perspective for the future developments in this field is given, particularly in light of recent developments such as the new version of the model, Martini 3.

Stable polymer bilayers for protein channel recordings at high guanidinium chloride concentrations

The use of chaotropic reagents is common in biophysical characterization of biomolecules. When the study involves transmembrane protein channels, the stability of the protein channel and supporting bilayer membrane must be considered. In this letter, we show that planar bilayers composed of poly(1,2-butadiene)-b-poly(ethylene oxide) diblock copolymer are stable and leak-free at high guanidinium chloride concentrations, in contrast to diphytanoyl phosphatidylcholine bilayers, which exhibit deleterious leakage under similar conditions. Furthermore, insertion and functional analysis of channels such as α-hemolysin and MspA are straightforward in these polymer membranes. Finally, we demonstrate that α-hemolysin channels maintain their structural integrity at 2 M guanidinium chloride concentrations using blunt DNA hairpins as molecular reporters.

Molecular simulation of the shape deformation of a polymersome

Vesicles composed of diblock copolymers, or polymersomes, have proven to possess numerous applications ranging from drug delivery to catalytically driven nano-motors. The shape of a polymersome can be responsive to external stimuli, such as light or solvent. Molecular dynamics simulations reveal that the shape change upon the contraction of the inner volume of a polymersome vesicle occurs in two separate regimes-a stretching regime and a bending regime. The barrier is shown to be dependent on the solvent environment. These results suggest that tailoring the bending modulus of polymer membranes can be used as a design methodology to engineer new stimuli-responsive vesicles.

Study of the Lamellar and Micellar Phases of Pluronic F127: A Molecular Dynamics Approach

In this work, we analyzed the behavior of Pluronic F127 through molecular dynamics simulations at the coarse-grain level, focusing on the micellar and lamellar phases. To this aim, two initial polymer conformations were considered, S-shape and U-shape, for both simulated phases. Through the simulations, we were able to examine the structural and mechanical properties that are difficult to access through experiments. Since no transition between S and U shapes was observed in our simulations, we inferred that all single co-polymers had memory of their initial configuration. Nevertheless, most copolymers had a more complex amorphous structure, where hydrophilic beads were part of the lamellar-like core. Finally, an overall comparison of the micellar a lamellar phases showed that the lamellar thickness was in the same order of magnitude as the micelle diameter (approx. 30 nm). Therefore, high micelle concentration could lead to lamellar formation. With this new information, we could understand lamellae as orderly packed micelles.

Mechanical properties of drug loaded diblock copolymer bilayers: A molecular dynamics study

In this work, we present results of coarse-grained simulations to study the encapsulation of prilocaine (PLC), both neutral and protonated, on copolymer bilayers through molecular dynamics simulations. Using a previously validated membrane model, we have simulated loaded bilayers at different drug concentrations and at low (protonated PLC) and high (neutral PLC) pH levels. We have characterized key structural parameters of the loaded bilayers in order to understand the effects of encapsulation of PLC on the bilayer structure and mechanical properties. Neutral PLC was encapsulated in the hydrophobic region leading to a thickness increase, while the protonated species partitioned between the water phase and the poly(ethylene oxide)-poly(butadiene) (PBD) interface, relaxing the PBD region and leading to a decrease in the thickness. The tangential pressures of the studied systems were calculated, and their components were decomposed in order to gain insights on their compensation. In all cases, it is observed that the loading of the membrane does not significantly decrease the stability of the bilayer, indicating that the system could be used for drug delivery.