New Insight into Cluster Aggregation Mechanism during Polymerization-Induced Self-Assembly by Molecular Dynamics Simulation

Modern Trends in Polymerization-Induced Self-Assembly

Polymerization-induced self-assembly (PISA) is a powerful and versatile technique for producing colloidal dispersions of block copolymer particles with desired morphologies. Currently, PISA can be carried out in various media, over a wide range of temperatures, and using different mechanisms. This method enables the production of biodegradable objects and particles with various functionalities and stimuli sensitivity. Consequently, PISA offers a broad spectrum of potential commercial applications. The aim of this review is to provide an overview of the current state of rational synthesis of block copolymer particles with diverse morphologies using various PISA techniques and mechanisms. The discussion begins with an examination of the main thermodynamic, kinetic, and structural aspects of block copolymer micellization, followed by an exploration of the key principles of PISA in the formation of gradient and block copolymers. The review also delves into the main mechanisms of PISA implementation and the principles governing particle morphology. Finally, the potential future developments in PISA are considered.

Stimuli-responsive viscosity modifiers

Stimuli responsive viscosity modifiers entail an important class of materials which allow for smart material formation utilizing various stimuli for switching such as pH, temperature, light and salinity. They have seen applications in the biomedical space including tissue engineering and drug delivery, wherein stimuli responsive hydrogels and polymeric vessels have been extensively applied. Applications have also been seen in other domains like the energy sector and automobile industry, in technologies such as enhanced oil recovery. The chemistry and microstructural arrangements of the aqueous morphologies of dissolved materials are usually sensitive to the aforementioned stimuli which subsequently results in rheological sensitivity as well. Herein, we overview different structures capable of viscosity modification as well as go over the rheological theory associated with classical systems studied in literature. A detailed analysis allows us to explore correlations between commonly discussed models such as molecular packing parameter, tube reptation and stress relaxation with structural and rheological changes. We then present five primary mechanisms corresponding to stimuli responsive viscosity modification: (i) packing parameter modification via functional group conditioning and (ii) via dynamic bond formation, (iii) mesh formation by interlinking of network nodes, (iv) viscosity modification by chain conformation changes and (v) viscosity modification by particle jamming. We also overview several recent examples from literature that employ the concepts discussed to create novel classes of intriguing stimuli responsive structures and their corresponding rheological properties. Furthermore, we also explore systems that are responsive to multiple stimuli which can provide enhanced functionality and versatility by providing multi-level and precise actuation. Such systems have been used for programmed site-specific drug delivery.

Does freezing induce self-assembly of polymers? A molecular dynamics study

This work investigates the freezing-induced self-assembly (FISA) of polyvinyl alcohol (PVA) and PVA-like polymers using molecular dynamics simulations. In particular, the effect of the degree of supercooling, degree of polymerization, polymer type, and initial local concentration on the FISA was studied. It was found that the preeminent factor responsible for FISA is not the diffusion of the polymers away from the nucleating ice front, but the increase in the polymer's local concentration upon freezing of the solvent (water). At a higher degree of supercooling, the polymers are engulfed by the growing ice front, impeding their diffusion into the supercooled solution and finally inhibiting their self-assembly. Conversely, at a relatively lower degree of supercooling, the rate of diffusion of the polymers into the supercooled solution is higher, which increases their local concentration and results in FISA. FISA was also observed to depend on the polymer-solvent interactions. Strongly favorable solute-solvent interactions hinder the self-assembly, whereas unfavorable solute-solvent interactions promote the self-assembly. The polymer and aggregate morphology were investigated using the radius of gyration, end-to-end distance, and asphericity analysis. This study brings molecular insights into the quintessential factors governing self-assembly via freezing of the solvent, which is a novel self-assembly technique especially suitable for biomedical applications.

Computational mesoscale framework for biological clustering and fractal aggregation

Hierarchical clustering due to diffusion and reaction is a widespread occurrence in natural phenomena, displaying fractal behavior with non-integer size scaling. The study of this phenomenon has garnered interest in both biological systems such as morphogenesis and blood clotting, and synthetic systems such as colloids and polymers. The modeling of biological clustering can be difficult, as it can occur on a variety of scales and involve multiple mechanisms, necessitating the use of various methods to capture its behavior. Here, we propose a novel framework, the generalized-mesoscale-clustering (GMC), for the study of complex hierarchical clustering phenomena in biological systems. The GMC framework incorporates the effects of hydrodynamic interactions, bonding, and surface tension, and allows for the analysis of both static and dynamic states of cluster development. The framework is applied to a range of biological clustering mechanisms, with a focus on blood-related clustering from fibrin network formation to platelet aggregation. Our study highlights the importance of a comprehensive characterization of the structural properties of the cluster, including fractal dimension, pore-scale diffusion, initiation time, and consolidation time, in fully understanding the behavior of biological clustering systems. The GMC framework also provides the potential to investigate the temporal evolution and mechanical properties of the clusters by tracking bond density and including hydrodynamic interactions.

Morphology Control via RAFT Emulsion Polymerization-Induced Self-Assembly: Systematic Investigation of Core-Forming Blocks

Polymerization-induced self-assembly (PISA) is a useful formulation for readily obtaining nanoparticles from block copolymers in situ. Reversible addition-fragmentation chain-transfer (RAFT) emulsion polymerization is utilized as one of the PISA formulations. Various factors have so far been investigated for obtaining nonspherical particles via RAFT emulsion polymerization, such as the steric structure of the shell, the glass-transition temperature (T _g) of the core-forming block, and the water solubility of the core-forming monomer. This study focuses on core-forming blocks without changing the structure of the shell-forming block. In particular, we elucidate the balance between T _g for the core-forming block and the water solubility of the core monomer. A series of alkyl methacrylates, such as methyl methacrylate (MMA), ethyl methacrylate (EMA), and n-propyl methacrylate (PrMA), are emulsion-polymerized in the presence of a poly[poly(ethylene glycol) methyl ether methacrylate] (PPEGMA) macromolecular chain-transfer agent via the RAFT process. The resulting in situ morphology changes to form shapes such as spheres, worms (toroids), and vesicles are systematically investigated. The properties of the core that determine whether a morphological change occurs from spheres are (i) the solubility of the core-forming monomer in water, (ii) the relationship between T _g for the core-forming block and the polymerization temperature, and (iii) the hydrophobic core volume, which changes the packing parameter. These factors allow prediction of the block copolymer morphology produced during RAFT emulsion polymerization of other methacrylates such as n-butyl methacrylate (BuMA), tetrahydrofurfuryl methacrylate (THFMA) with physical properties of the homopolymer (poly(tetrahydrofurfuryl methacrylate) (PTHFMA)) between those for poly(MMA) (PMMA) and PBuMA, and 1-adamantyl methacrylate (ADMA) with low monomer solubility in water and high T _g of the homopolymer (PADMA).

Coarse-Grained Molecular Dynamics Simulations of Paclitaxel-Loaded Polymeric Micelles

A continuous manufacturing technology based on coaxial turbulent jet in coflow was previously developed to produce paclitaxel-loaded polymeric micelles. Herein, coarse-grained molecular dynamics (CG-MD) simulations were implemented to better understand the effect of the material attributes (i.e., the drug-polymer ratio and the ethanol concentration) and process parameters (i.e., temperature) on the self-assembly process of polymeric micelles as well as to provide molecular details on micelle instability. An all-atom (AA) poly (ethylene glycol)-poly (lactic acid) (PEG-PLA) polymer model was developed as the reference for parameterizing a coarse-grained (CG) model, and the AA polymer model was further validated with experimental glass transition temperature (T_g). The model transferability was verified by comparing structural properties between the AA and CG models. The CG model was further validated with experimental data, including micelle particle size measurements and drug encapsulation efficiency. Furthermore, the encapsulation of paclitaxel into the polymeric micelles was included in the simulations, taking into consideration the interactions between the paclitaxel and the polymers. The results from various points of view demonstrated a strong dependence of the shape of the micelles on the drug encapsulation, with micelles transitioning from spherical to ellipsoidal structures with an increasing paclitaxel amount. Simulation data were also used to identify the critical aggregation number (i.e., the number of polymer and drug molecules required for transition from one shape to another). Improved micellar structural stability was found with a larger micellar size and less solvent accessibility. Lastly, an evaluation was performed on the micellar dissociation free energy using a steered molecular dynamics simulation over a range of temperatures and ethanol concentrations. These simulations revealed that at higher ethanol and temperature conditions, micelles become destabilized, resulting in greater paclitaxel release. The increased drug release was determined to originate from the solvation of the hydrophobic core, which promoted micellar swelling and an associated reduction in hydrophobic interactions, leading to a loosely packed micellar structure.

Selective solvent conditions influence sequence development and supramolecular assembly in step-growth copolymerization

Sequence control in synthetic copolymers remains a tantalizing objective in polymer science due to the influence of sequence on material properties and self-organization. A greater understanding of sequence development throughout the polymerization process will aid the design of simple, generalizable methods to control sequence and tune supramolecular assembly. In previous simulations of solution-based step-growth copolymerizations, we have shown that weak, non-bonding attractions between monomers of the same type can produce a microphase separation among the lengthening nascent oligomers and thereby alter sequence. This work explores the phenomenon further, examining how effective attractive interactions, mediated by a solvent selective for one of the reacting species, impact the development of sequence and the supramolecular assembly in a simple A-B copolymerization. We find that as the effective attractions between monomers increase, an emergent self-organization of the reactants causes a shift in reaction kinetics and sequence development. When the solvent-mediated interactions are selective enough, the simple mixture of A and B monomers oligomerize and self-assemble into structures characteristic of amphiphilic copolymers. The composition and morphology of these structures and the sequences of their chains are sensitive to the relative balance of affinities between the comonomer species. Our results demonstrate the impact of differing A-B monomer-solvent affinities on sequence development in solution-based copolymerizations and are of consequence to the informed design of synthetic methods for sequence controlled amphiphilic copolymers and their aggregates.

Synthesis of low glass transition temperature worms comprising a poly(styrene-stat-n-butyl acrylate) core segment via polymerization-induced self-assembly in RAFT aqueous emulsion polymerization

Synthesis of nanodimensional polymeric worms of low glass transition temperature using aqueous polymerization-induced self-assembly.

Improving the self-assembly of bioresponsive nanocarriers by engineering doped nanocarbons: a computational atomistic insight

Here, molecular dynamics (MD) simulations were employed to explore the self-assembly of polymers and docetaxel (DTX) as an anticancer drug in the presence of nitrogen, phosphorous, and boron-nitrogen incorporated graphene and fullerene. The electrostatic potential and the Gibbs free energy of the self-assembled materials were used to optimize the atomic doping percentage of the N- and P-doped formulations at 10% and 50%, respectively. Poly lactic-glycolic acid (PLGA)- polyethylene glycol (PEG)-based polymeric nanoparticles were assembled in the presence of nanocarbons in the common (corresponding to the bulk environment) and interface of organic/aqueous solutions (corresponding to the microfluidic environment). Assessment of the modeling results (e.g., size, hydrophobicity, and energy) indicated that among the nanocarbons, the N-doped graphene nanosheet in the interface method created more stable polymeric nanoparticles (PNPs). Energy analysis demonstrated that doping with nanocarbons increased the electrostatic interaction energy in the self-assembly process. On the other hand, the fullerene-based nanocarbons promoted van der Waals intramolecular interactions in the PNPs. Next, the selected N-doped graphene nanosheet was utilized to prepare nanoparticles and explore the physicochemical properties of the nanosheets in the permeation of the resultant nanoparticles through cell-based lipid bilayer membranes. In agreement with the previous results, the N-graphene assisted PNP in the interface method and was translocated into and through the cell membrane with more stable interactions. In summary, the present MD simulation results demonstrated the success of 2D graphene dopants in the nucleation and growth of PLGA-based nanoparticles for improving anticancer drug delivery to cells, establishing new promising materials and a way to assess their performance that should be further studied.

Simulation of polymerization induced phase separation in model thermosets

Polymerization induced phase separation (PIPS) in a three component thermoset is studied using molecular dynamics simulations of a new coarse-grained thermoset model. The system includes two crosslinker molecules, which differ in their glass transition temperatures (T_g) and chain length and thus have the potential for phase separation. One crosslinker has a high T_g corresponding to a rubbery behavior, and simulations were performed for a short length (4 beads) and a long length (33 beads). The resin and other crosslinker have low T_g. A coarse-grained model is developed with these features and with interaction parameters determined so that for either rubbery crosslinker length, the system is in the liquid state at the cure temperature. For sufficiently slow reaction rates, the long rubbery molecule exhibits PIPS into a bicontinuous array of nanoscale domains, but the short one does not, reproducing recent experimental results. The simulations demonstrate that the reaction rates must be slow enough to allow diffusion to yield phase separation. Particularly, the reaction rate corresponding to the secondary amine must be very slow, else the structure of crosslinked clusters and the substantially increased diffusion time will prevent PIPS.

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.

Elucidating preparation-structure relationships for the morphology evolution during the RAFT dispersion polymerization of N-acryloyl thiomorpholine

Depending on the conditions, nearly monodisperse spherical micelles or complex morphologies are formed during a polymerization induced self-assembly (PISA) process based on the water-soluble monomer N-acryloylthiomorpholine.

In situ SAXS studies of a prototypical RAFT aqueous dispersion polymerization formulation: monitoring the evolution in copolymer morphology during polymerization-induced self-assembly

Small-angle X-ray scattering (SAXS) is used to characterize the in situ formation of diblock copolymer spheres, worms and vesicles during reversible addition-fragmentation chain transfer (RAFT) aqueous dispersion polymerization of 2-hydroxypropyl methacrylate at 70 °C using a poly(glycerol monomethacrylate) steric stabilizer. 1H NMR spectroscopy indicates more than 99% HPMA conversion within 80 min, while transmission electron microscopy and dynamic light scattering studies are consistent with the final morphology being pure vesicles. Analysis of time-resolved SAXS patterns for this prototypical polymerization-induced self-assembly (PISA) formulation enables the evolution in copolymer morphology, particle diameter, mean aggregation number, solvent volume fraction, surface density of copolymer chains and their mean inter-chain separation distance at the nanoparticle surface to be monitored. Furthermore, the change in vesicle diameter and membrane thickness during the final stages of polymerization supports an 'inward growth' mechanism.

Poly(ethylene glycol)-b-poly(vinyl acetate) block copolymer particles with various morphologies via RAFT/MADIX aqueous emulsion PISA

The polymerization-induced self-assembly (PISA) of amphiphilic diblock copolymers of poly(ethylene glycol)-b-poly(vinyl acetate) in water was achieved through macromolecular design via interchange of xanthate (MADIX) polymerization in emulsion.

Unravelling the formation of BAB block copolymer assemblies during PISA in water

BAB triblock copolymers prepared by PISA in water self-assemble into a transient network of bridged micelles. The slowdown of the exchange of B blocks between micelles during PISA is highlighted as well as the parameters affecting the polymerization.

Synthesis and direct assembly of linear–dendritic copolymers via CuAAC click polymerization-induced self-assembly (CPISA)

A one-pot method was developed for in situ preparation of linear–dendritic copolymer assemblies via click polymerization-induced self-assembly (CPISA).