Kinetically-arrested single-polymer nanostructures from amphiphilic mikto-grafted bottlebrushes in solution: a simulation study

Self-assembly prediction of architecture-controlled bottlebrush copolymers in solution using graph convolutional networks

The investigation of bottlebrush copolymer self-assembly in solution involves a comprehensive approach integrating simulation and experimental research, due to their unique physical characteristics. However, the intricate architecture of bottlebrush copolymers and the diverse solvent conditions introduce a wide range of parameter spaces. In this study, we investigated the solution self-assembly behavior of bottlebrush copolymers by combining dissipative particle dynamics (DPD) simulation results and machine learning (ML) including graph convolutional networks (GCNs). The architecture of bottlebrush copolymers is encoded by graphs including connectivity, side chain length, bead types, and interaction parameters of DPD simulation. Using GCN, we accurately predicted the single chain properties of bottlebrush copolymers with over 95% accuracy. Furthermore, phase behavior was precisely predicted using these single chain properties. Shapley additive explanations (SHAP) values of single chain properties to the various self-assembly morphologies were calculated to investigate the correlation between single chain properties and morphologies. In addition, we analyzed single chain properties and phase behavior as a function of DPD interaction parameters, extracting relevant physical properties for vesicle morphology formation. This work paves the way for tailored design in solution of self-assembled nanostructures of bottlebrush copolymers, offering a GCN framework for precise prediction of self-assembly morphologies under various chain architectures and solvent conditions.

Self-assembly of rod-coil-rod block copolymers in a coil-selective solvent: coarse-grained simulation results

The solution self-assembly of amphiphilic polymers provides a versatile approach to design novel nanostructured materials. Multiblock polymers, particularly those composed of liquid crystalline and coil blocks, are of significant interest due to the potential display of nematic ordering in liquid crystalline domains, offering intriguing optical and mechanical properties. In this study, dissipative particle dynamics is used to investigate the solution self-assembly of rod-coil-rod copolymers in a coil-selective solvent. Extensive molecular simulations were conducted to elucidate the impact of polymer composition, concentration and flexibility on the self-assembly behavior. A quantitative analysis was performed to investigate how polymer conformations varied with changes in composition, concentration, and rigidity. Simulation results show that, at small rod compositions, rod-coil-rod polymers self-assemble into micelles at low concentrations, transitioning to network formation as concentration increases. An increase in rod composition leads to the formation of larger aggregates, resulting in cylindrical micelles and membranes. The results reported here also offer insights into the role of flexibility in shaping the self-assembly behavior of rod-coil-rod triblocks in selective solvents, thus, contributing to a comprehensive understanding of the factors governing the formation of diverse structures in the solution self-assembly of triblock copolymers.

Phase Behavior and Conformational Asymmetry near the Comb-to-Bottlebrush Transition in Linear-Brush Block Copolymers

This study explores how conformational asymmetry influences the bulk phase behavior of linear-brush block copolymers. We synthesized 60 diblock copolymers composed of poly(trifluoroethyl methacrylate) as the linear block and poly[oligo(ethylene glycol) methyl ether methacrylate] as the brush block, varying the molecular weight, composition, and side-chain length to introduce different degrees of conformational asymmetry. Using small-angle X-ray scattering, we determined the morphology and phase diagrams for three different side-chain length systems, mainly observing lamellar and cylindrical phases. Increasing the side-chain length of the brush block from three to nine ethylene oxide units introduces sufficient asymmetry between the blocks to alter the phase behavior, shifting the lamellar-to-cylindrical transitions toward lower brush block compositions and transitioning the brush block from the dense comb-like regime to the bottlebrush regime. Coarse-grained simulations support our experimental observations and provide a mapping between the composition and conformational asymmetry. A comparison of our findings to strong stretching theory across multiple phase boundary predictions confirms the transition between the dense comb-like regime and the bottlebrush regime.

Worm-globule transition of amphiphilic pH-responsive heterografted bottlebrushes at air-water interface

Heterografted molecular bottlebrushes (MBBs) with side chains composed of poly(n-butyl acrylate) (PnBA) and pH-responsive poly(2-(N,N-diethylamino)ethyl methacrylate) (PDEAEMA, pKa = 7.4) have been shown to be efficient, robust, and responsive emulsifiers. However, it remains unknown how they respond to external stimuli at interfaces. In this work, the shape-changing behavior of six hetero- and homografted MBBs at air-water interfaces in response to pH changes and lateral compression was investigated using a Langmuir-Blodgett trough and atomic force microscopy. At a surface pressure of 0.5 mN m-1, PDEAEMA-containing MBBs showed no worm-globule transitions when the pH was increased from 4.0 to 10.0, at which PDEAEMA becomes insoluble in water. Upon lateral compression at pH 4.0, MBBs with a mole fraction of PDEAEMA side chains (xPDEAEMA) < 0.50 underwent pronounced worm-globule shape transitions; there was an increasing tendency for bottlebrushes to become connected with increasing xPDEAEMA. At xPDEAEMA = 0.76, the molecules remained wormlike even at high compression. These observations were presumably caused by the increased electrostatic repulsion between protonated PDEAEMA side chains in the subphase with increasing xPDEAEMA, hindering the shape change. At pH 10.0, MBBs with xPDEAEMA < 0.50 showed a lower tendency to change their wormlike morphologies upon compression than at pH 4.0. No shape transition was observed when xPDEAEMA > 0.50, attributed to the relatively high affinity toward water and the rigidity of PDEAEMA. This study revealed the shape-changing behavior of amphiphilic pH-responsive MBBs at air-water interfaces, which could be useful for future design of multicomponent MBBs for potential applications.

Amphiphilic Heterograft Copolymers Bearing Biocompatible/Biodegradable Grafts

The amphiphilic heterograft copolymers bearing biocompatible/biodegradable grafts [poly(2-methyl-2-oxazoline-co-2-pentyl-2-oxazoline)-g-poly(d-l-lactic acid)/poly(2-ethyl-2-oxazoline)] were synthesized successfully by the combination of cationic ring-opening polymerization and click chemistry via the ⟨"grafting to"⟩ approach. The challenge of this synthesis was to graft together hydrophobic and hydrophilic chains on a hydrophilic platform based on PMeOx. The efficiency of grafting depends on the chemical nature of the grafts and of the length of the macromolecular chains. The self-assembly of these polymers in aqueous media was investigated by DLS, cryo-TEM, and SANS. The results demonstrated that different morphologies were obtained from nanospheres and vesicles to filaments depending on the hydrophilic weight ratio in the heterograft copolymer varying from 0.38 until 0.84. As poly(2-ethyl-2-oxazoline) is known to be thermoresponsive, the influence of temperature rise on the nanoassembly stability was studied in water and in a physiological medium. SANS and DLS measurements during a temperature ramp allowed to show that nanoassemblies start to self-assemble in "raspberry like" primary structures at 50 °C, and these structures grow and get denser as the temperature is increased further. These amphiphilic heterograft copolymers may include hydrophobic drugs and should find important applications for biomedical applications which require stealth properties.

Exploring conformations of comb-like polymers with varying grafting density in dilute solutions

Comb-like polymers have shown potential as advanced materials for a diverse palette of applications due to the tunability of their polymer architecture. To date, however, it still remains a challenge to understand how the conformational properties of these polymers arise from the interplay of their architectural parameters. In this work, extensive simulations were performed using dissipative particle dynamics to investigate the effect of grafting density, backbone length, and sidechain length on the conformations of comb-like polymers immersed in a good solvent. To quantify the effect of these architectural parameters on polymer conformations, we computed the asphericity, radius of gyration, and backbone and sidechain end-to-end distances. Bond-bond correlation functions and effective Kuhn lengths were computed to quantify the topological stiffness induced by sidechain-sidechain interactions. Simulation results reveal that the effective Kuhn length increases as grafting density and sidechain length increase, in agreement with previous experimental and theoretical studies. This increase in stiffness results in comb-like polymers adopting extended conformations as grafting density and sidechain length increase. Simulation results regarding the radius of gyration of comb-like polymers as a function of grafting density are compared with scaling theory predictions based on a free energy proposed by Morozova and Lodge [ACS Macro Lett. 6, 1274-1279 (2017)] and scaling arguments by Tang et al. [Macromolecules 55, 8668-8675 (2022)].

Effect of Block Sequence on the Solution Self-Assembly of Symmetric ABCBA Pentablock Polymers in a Selective Solvent

Solution self-assembly of multiblock polymers offers a platform to create complex functional self-assembled nanostructures. However, a complete understanding of the effect of the different single-molecule-level parameters and solution conditions on the self-assembled morphology is still lacking. In this work, we have used dissipative particle dynamics to investigate the solution self-assembly of symmetric ABCBA linear pentablock polymers in a selective solvent and examined the effect of the block sequence, composition, and polymer concentration on the final morphology and polymer conformations. We confirmed that block sequence has an effect on the self-assembled morphologies, and it has a strong influence on polymer conformations that give place to physical gels for the sequence where the solvophilic block is located in the middle of the macromolecule. Our results are summarized in terms of morphology diagrams in the composition-concentration parameter space.

Morphological Transitions in Micelles of Amphiphilic Bottlebrushes upon the Adsorption and Compression at the Liquid Interface

Densely grafted comb-like macromolecules (bottlebrushes) with alternating solvophobic and solvophilic side chains were studied in a selective solvent and at the liquid interface using mesoscopic computer simulations. The effects of backbone length and copolymer composition were considered. While self-assembly in solution revealed only spherical aggregates for all ar-chitectures studied, adsorption onto the liquid interface in particular cases resulted in morpho-logical changes, with worm-like aggregates or a continuous monolayer observed. In turn, the compression of macromolecules at the interface also leads to morphological transitions, includ-ing the formation of a mesh-like percolated structure. The obtained results may be useful for the preparation of solid nanoparticles of anisotropic shape or nanostructured ultra-thin copolymer films.

Atypical vesicles and membranes with monolayer and multilayer structures formed by graft copolymers with diblock side-chains: nonlamellar structures and curvature-enhanced permeability

Graft copolymers with diblock side-chains A_m(-graft-B₃A_y)_n in a selective solvent have been reported to self-assemble into vesicles, but the structure is expected to differ distinctly from those of lipid bilayers. Surprisingly, the number of alternating hydrophobic A-block and hydrophilic B-block layers in the vesicle can vary from a monolayer to multilayers such as the hepta-layer, subject to the same copolymer concentration. The area density of the copolymer layer is not uniform across the membrane. This structural difference among different layers is attributed to the neighboring environment and the curvature of the layer. Because of the unusual polymer conformations, nonlamellar structures of polymersomes are formed, and they are much more intricate than those of liposomes. In fact, a copolymer can contribute to a single or two hydrophilic layers, and it can provide up to three hydrophobic layers. The influence of the backbone length (m) and side-chain length (y) and the permeation dynamics are also studied. The thickness of hydrophobic layers is found to increase with increasing side-chain length but is not sensitive to the backbone length. Although the permeation time increases with the layer number for planar membranes, the opposite behavior is observed for spherical vesicles owing to the curvature-enhanced permeability associated with Laplace pressure.

Influences of the Periodicity in Molecular Architecture on the Phase Diagrams and Microphase Transitions of the Janus Double-Brush Copolymer with a Loose Graft

The backbone of the Janus double-brush copolymer may break during long-term service, but whether this breakage affects the self-assembled phase state and microphase transitions of the material is still unknown. For the Janus double-brush copolymers with a periodicity in molecular architecture ranging from 1 to 10, the influences of the architectural periodicity on their phase diagrams and order–disorder transitions (ODT) were investigated by the self-consistent mean field theory (SCFT). In total, nine microphases with long-range order were found. By comparing the phase diagrams between copolymers of different periodicity, a decrease in periodicity or breakage along the copolymer backbone had nearly no influence on the phase diagrams unless the periodicity was too short to be smaller than 3. For copolymers with neutral backbones, a decrease in periodicity or breakage along the copolymer backbone reduced the critical segregation strengths of the whole copolymer at ODT. The equations for the critical segregation strengths at ODT, the architectural periodicity, and the volume fraction of the backbone were established for the Janus double-brush copolymers. The theoretical calculations were consistent with the previous theoretical, experimental, and simulation results.

DPD Modelling of the Self- and Co-Assembly of Polymers and Polyelectrolytes in Aqueous Media: Impact on Polymer Science

This review article is addressed to a broad community of polymer scientists. We outline and analyse the fundamentals of the dissipative particle dynamics (DPD) simulation method from the point of view of polymer physics and review the articles on polymer systems published in approximately the last two decades, focusing on their impact on macromolecular science. Special attention is devoted to polymer and polyelectrolyte self- and co-assembly and self-organisation and to the problems connected with the implementation of explicit electrostatics in DPD numerical machinery. Critical analysis of the results of a number of successful DPD studies of complex polymer systems published recently documents the importance and suitability of this coarse-grained method for studying polymer systems.

Effects of temperature on chaotropic anion-induced shape transitions of star molecular bottlebrushes with heterografted poly(ethylene oxide) and poly(N,N-dialkylaminoethyl methacrylate) side chains in acidic water

This article reports a study of the effects of temperature on chaotropic anion (CA)-induced star-globule shape transitions in acidic water of three-arm star bottlebrushes composed of heterografted poly(ethylene oxide) (PEO) and either poly(2-(N,N-dimethylamino)ethyl methacrylate) (PDMAEMA) or poly(2-(N,N-diethylamino)ethyl methacrylate) (PDEAEMA) (the brushes denoted as SMB-11 and -22, respectively). The brush polymers were synthesized by grafting alkyne-end-functionalized PEO and PDMAEMA or PDEAEMA onto an azide-bearing three-arm star backbone polymer using the copper(i)-catalyzed alkyne-azide cycloaddition reaction. Six anions were studied for their effects on the conformations of SMB-11 and -22 in acidic water: super CAs [Fe(CN)6]3- and [Fe(CN)6]4-, moderate CAs PF6- and ClO4-, weak CA I-, and for comparison, kosmotropic anion SO42-. At 25 °C, the addition of super and moderate CAs induced shape transitions of SMB-11 and -22 in pH 4.50 water from a starlike to a collapsed globular state stabilized by PEO side chains, which was driven by the ion pairing of protonated tertiary amine groups with CAs and the chaotropic effect. The shape changes occurred at much lower salt concentrations for super CAs than moderate CAs. Upon heating from near room temperature to 70 °C, the super CA-collapsed brushes remained in the globular state, whereas the moderate CA-collapsed brushes underwent reversible globule-to-star shape transitions. The transition temperature increased with increasing salt concentration and was found to be higher for SMB-22 at the same salt concentration, presumably caused by the chaotropic effect. In contrast, I- and SO42- had small effects on the conformations of SMB-11 and -22 at 25 °C in the studied salt concentration range, and only small and gradual size variations were observed upon heating to 70 °C. The results reported here may have potential uses in the design of stimuli-responsive systems for substance encapsulation and release.

Shape-Changing Bottlebrush Polymers

Bottlebrush polymers (BBPs), composed of relatively short polymeric side chains densely grafted on a polymer backbone, exhibit many unique characteristics and hold promise for a variety of applications. This Perspective focuses on environmentally induced shape-changing behavior of BBPs at interface and in solution, particularly worm/star-globule shape transitions. While BBPs with a single type of homopolymer or random copolymer side chains have been shown to undergo pronounced worm-to-globule shape changes in response to external stimuli, the collapsed brushes are unstable and prone to aggregation. By introducing a second, solvophilic polymer into the side chains, either as a distinct type of side chain or as the outer block of block copolymer side chains, the collapsed brushes not only are stabilized but also create unimolecular micellar nanostructures, which can be used for, e.g., encapsulation and delivery of substances. The current challenges in the design, synthesis, and characterization of stimuli-responsive shape-changing BBPs are discussed.

Using reactive dissipative particle dynamics to understand local shape manipulation of polymer vesicles

Biological cells have long been of interest to researchers due to their capacity to actively control their shape. Accordingly, there is significant interest in generating simplified synthetic protocells that can alter their shape based on an externally or internally generated stimulus. To date, most progress has been made towards controlling the global shape of a protocell, whereas less is known about generating a local shape change. Here, we seek to better understand the possible mechanisms for producing local morphological changes in a popular protocell system, the block copolymer vesicle. Accordingly, we have combined Dissipative Particle Dynamics (DPD) and the Split Reactive Brownian Dynamics algorithm (SRBD) to produce a simulation tool that is capable of modeling the dynamics of self-assembled polymer structures as they undergo chemical reactions. Using this Reactive DPD or RDPD method, we investigate local morphological change driven by either the microinjection of a stimulus or an enzymatically-produced stimulus. We find that sub-vesicle-scale morphological change can be induced by either a solvent stimulus that swells the vesicle membrane, or by a reactant stimulus that alters the chemistry of the block polymer in the membrane corona. Notably, the latter method results in a more persistent local deformation than the former, which we attribute to the slower diffusion of polymer chains relative to the solvent. We quantify this deformation and show that it can be modulated by altering the interaction parameter of the parts of the polymer chain that are affected by the stimulus.

Effect of side-chain length on solute encapsulation by amphiphilic heterografted brush copolymers

Anisotropic nanomaterials are non-spherical structures that possess unique shape-dependent physicochemical properties and functionalities. Inspired by the abundance of filamentous entities in nature, cylindrical nanostructures have gained significant attention due to their unique performance. Herein, we discuss the effect of side-chain length on the encapsulation properties of amphiphilic heterografted bottlebrushes. We observed that by grafting a long hydrophilic block to the double-brush, we were able to restrict solute-induced conformational changes, thus producing drug-loaded anisotropic carriers. Unimolecular encapsulation in brushes was solute-dependent as shown here for probucol and rose bengal lactone. Stabilization with an amphiphilic diblock copolymer-consisting of the same type of blocks as those comprising the heterografted brush-served to explain the solute-dependent behavior observed for brushes, suggesting that solutes with a higher propensity to nucleation could be more effectively stabilized by the anisotropic carrier in a unimolecular worm-like construct.