Computer Aided Modeling of Polymer Reaction Engineering-The Status of Predici, I-Simulation

Sequence Design of Random Heteropolymers as Protein Mimics

Random heteropolymers (RHPs) have been computationally designed and experimentally shown to recapitulate protein-like phase behavior and function. However, unlike proteins, RHP sequences are only statistically defined and cannot be sequenced. Recent developments in reversible-deactivation radical polymerization allowed simulated polymer sequences based on the well-established Mayo-Lewis equation to more accurately reflect ground-truth sequences that are experimentally synthesized. This led to opportunities to perform bioinformatics-inspired analysis on simulated sequences to guide the design, synthesis, and interpretation of RHPs. We compared batches on the order of 10000 simulated RHP sequences that vary by synthetically controllable and measurable RHP characteristics such as chemical heterogeneity and average degree of polymerization. Our analysis spans across 3 levels: segments along a single chain, sequences within a batch, and batch-averaged statistics. We discuss simulator fidelity and highlight the importance of robust segment definition. Examples are presented that demonstrate the use of simulated sequence analysis for in-silico iterative design to mimic protein hydrophobic/hydrophilic segment distributions in RHPs and compare RHP and protein sequence segments to explain experimental results of RHPs that mimic protein function. To facilitate the community use of this workflow, the simulator and analysis modules have been made available through an open source toolkit, the RHPapp.

Modelling altered signalling of G-protein coupled receptors in inflamed environment to advance drug design

We previously reported the successful design, synthesis and testing of the prototype opioid painkiller NFEPP that does not elicit adverse side effects. The design process of NFEPP was based on mathematical modelling of extracellular interactions between G-protein coupled receptors (GPCRs) and ligands, recognizing that GPCRs function differently under pathological versus healthy conditions. We now present an additional and novel stochastic model of GPCR function that includes intracellular dissociation of G-protein subunits and modulation of plasma membrane calcium channels and their dependence on parameters of inflamed and healthy tissue (pH, radicals). The model is validated against in vitro experimental data for the ligands NFEPP and fentanyl at different pH values and radical concentrations. We observe markedly reduced binding affinity and calcium channel inhibition for NFEPP at normal pH compared to lower pH, in contrast to the effect of fentanyl. For increasing radical concentrations, we find enhanced constitutive G-protein activation but reduced ligand binding affinity. Assessing the different effects, the results suggest that, compared to radicals, low pH is a more important determinant of overall GPCR function in an inflamed environment. Future drug design efforts should take this into account.

Toward Green Atom Transfer Radical Polymerization: Current Status and Future Challenges

Reversible-deactivation radical polymerizations (RDRPs) have revolutionized synthetic polymer chemistry. Nowadays, RDRPs facilitate design and preparation of materials with controlled architecture, composition, and functionality. Atom transfer radical polymerization (ATRP) has evolved beyond traditional polymer field, enabling synthesis of organic-inorganic hybrids, bioconjugates, advanced polymers for electronics, energy, and environmentally relevant polymeric materials for broad applications in various fields. This review focuses on the relation between ATRP technology and the 12 principles of green chemistry, which are paramount guidelines in sustainable research and implementation. The green features of ATRP are presented, discussing the environmental and/or health issues and the challenges that remain to be overcome. Key discoveries and recent developments in green ATRP are highlighted, while providing a perspective for future opportunities in this area.

Simulation of the RAFT polymerization in 3D: steric restrictions and incompatibility between species

In this work we developed a RAFT polymerization model taking into account the main reactions of the experimental RAFT process and implemented that model in dissipative particle dynamics (DPD). With...

Chain-Length Dependence of Propagation Rate Coefficient for Methyl Acrylate Polymerization at 25 °C Investigated by PLP-SEC Method

When applying the pulsed laser polymerization–size exclusion chromatography (PLP-SEC) method to determine the propagation rate coefficient (kp) for bulk methyl acrylate at 25 °C, the characteristic ratio L1/L2 of chain...

Pulsed Laser Polymerization – Size Exclusion Chromatography Investigations into Backbiting in Ethylhexyl Acrylate Polymerization

Free radical acrylate polymerization is associated with a considerable number of side reactions occurring during chain growth. Most notably, the secondary chain propagating radicals can undergo intramolecular transfer reactions, forming...

Mechanistic investigation of cyclic ketene acetal radical ring-opening homo- and co-polymerization and preparation of PEO graft copolymers with tunable composition

Radical ring-opening polymerization of cyclic ketene acetals (CKAs) led to homogeneous ester insertion during semi-batch reactions with N-vinyl pyrrolidone (NVP). Investigation of CKA radical transfer allowed for CKA–NVP graft copolymer preparation.

A machine-readable online database for rate coefficients in radical polymerization

An online database created and curated by an IUPAC subcommittee is introduced.

A Thermodynamic Roadmap for the Grafting-through Polymerization of PDMS11MA

Grafting-through atom transfer radical polymerization (ATRP) was used to polymerize a sterically hindered poly(dimethylsiloxane) methacrylate (PDMS₁₁MA, M_n = 1000) macromonomer to high conversion as a function of temperature, solvent, initial monomer concentration, and pressure. Higher polymerization yields were obtained when polymerizations were conducted at (i) lower temperature (T), (ii) in a poor solvent for the side chain, (iii) higher initial monomer concentration ([M]₀), and (iv) higher pressure by mitigating the contribution of the equilibrium monomer concentration ([M]_eq). The enthalpy of polymerization (ΔH_p) and entropy of polymerization (ΔS_p) were more negative in poor solvents. Polymerizations at ambient pressure required higher [M]₀, use of a poor solvent, and lower temperatures to reach higher conversion with good control, whereas high pressure ATRP (HP-ATRP) displayed better control under dilute conditions. Grafting-through polymerization at high P and higher [M]₀ was less controlled, plausibly due to limited solubility and mobility of the copper catalyst in the highly viscous medium.

p-Substituted Tris(2-pyridylmethyl)amines as Ligands for Highly Active ATRP Catalysts: Facile Synthesis and Characterization

A facile and efficient two-step synthesis of p-substituted tris(2-pyridylmethyl)amine (TPMA) ligands to form Cu complexes with the highest activity to date in atom transfer radical polymerization (ATRP) is presented. In the divergent synthesis, p-Cl substituents in tris(4-chloro-2-pyridylmethyl)amine (TPMA^3Cl ) were replaced in one step and high yield by electron-donating cyclic amines (pyrrolidine (TPMA^PYR ), piperidine (TPMA^PIP ), and morpholine (TPMA^MOR )) by nucleophilic aromatic substitution. The [Cu^II (TPMA^NR2 )Br]⁺ complexes exhibited larger energy gaps between frontier molecular orbitals and >0.2 V more negative reduction potentials than [Cu^II (TPMA)Br]⁺ , indicating >3 orders of magnitude higher ATRP activity. [Cu^I (TPMA^PYR )]⁺ exhibited the highest reported activity for Br-capped acrylate chain ends in DMF, and moderate activity toward C-F bonds at room temperature. ATRP of n-butyl acrylate using only 10-25 part per million loadings of [Cu^II (TPMA^NR2 )Br]⁺ exhibited excellent control.

Consequences of End-Group Fidelity for the Quantitative Analysis of Surface Grafting of Polymers

End-functional polymers are employed in a wide range of functionalization reactions, including, among others, polymer surface grafting. Herein, kinetic (PREDICI) simulations are used to investigate to what extent the molar mass distribution (MMD) of α-end-functional polymers derived from RAFT polymerization may differ from the MMD of nonfunctional chains and the overall polymeric material. The results indicate that the MMD of the overall polymeric material (as commonly accessed in experiments) may not provide a good estimate for the MMD of α-end-functional chains, if polymerization conditions deviate strongly from pseudoliving conditions. Careful consideration of this behavior is required when using α-end-functional RAFT polymers in quantitative studies.

New Aspects on the Modeling of Dithiolactone-Mediated Radical Polymerization of Vinyl Monomers

A kinetic model for the dithiolactone-mediated radical polymerization of vinyl monomers based on the persistent radical effect and reversible addition (negligible fragmentation) was used to calculate the polymerization rate and describe molar mass development in the polymerization of methyl methacrylate at 60 °C, using 2,2-azobisisobutyronitrile (AIBN) as an initiator, as well as dihydro-5-phenyl-2(3H)-thiophenethione (DTL1) and dihydro-2(3H)-thiophenethione (DTL2) as controllers. The model was implemented in the PREDICI commercial software. A good agreement between experimental data and model predictions was obtained.

On the Recovery of PLP-Molar Mass Distribution at High Laser Frequencies: A Simulation Study

Due to the inherent difficulties in determination of the degree of branching for polymers produced in pulsed laser polymerization (PLP) experiments, the behavior of the degree of branching and backbiting reaction in high laser frequency and relatively high reaction temperatures have not been well-established. Herein, through a simulation study, the validity of different explanations on the recovery of PLP-molar mass distribution at high laser frequencies is discussed. It is shown that the reduction of the backbiting reaction rate at high laser frequency, and consequent decrease in the degree of branching, is not a necessary condition for recovering the PLP-molar mass distribution. The findings of this work provide simulation support to a previous explanation about the possibility of using high laser frequency for reliable determination of the propagation rate coefficient for acrylic monomers.

Modeling of the Free Radical Copolymerization Kinetics of n-Butyl Acrylate, Methyl Methacrylate and 2-Ethylhexyl Acrylate Using PREDICI®

Kinetic modeling of the bulk free radical copolymerizations of n-butyl acrylate (BA) and 2-ethylhexyl acrylate (EHA); methyl methacrylate (MMA) and EHA; as well as BA, MMA and EHA was performed using the software PREDICI®. Predicted results of conversion versus time, composition versus conversion, and molecular weight development are compared against experimental data at different feed compositions. Diffusion-controlled effects and backbiting for BA were incorporated into the model as they proved to be significant in these polymerizations. The set of estimated global parameters allows one to assess the performance of these copolymerization systems over a wide range of monomer compositions.

Practical Prediction of Heteropolymer Composition and Drift

Composition drift in batch polymerizations is a well-known phenomenon and can lead to composition gradients in polymers synthesized using controlled polymerization methodologies. With known reactivity ratios of monomers, the drift, and thus resultant gradient copolymer, can be designed by adjusting reagent ratios and targeted conversions. Although such prediction is straightforward, it is seldom done, likely due to the perceived difficulty and unfamiliarity for nonspecialists. We seek to remedy this by providing the communities using copolymers with an easy-to-use program called Compositional Drift which is based on the Mayo-Lewis model and the penultimate model of monomer addition, using Monte Carlo methodology. This tool can also be applied to predict composition in nondrifting polymerizations. Herein we supply this tool to the community, showcasing two recent examples of use to guide experimental design and understanding of heteropolymers (RHP).

Modelling and control of the microstructure of comb-like poly(MAA-co-PEGMA) water-soluble copolymers

Modelling and controlling the cumulative copolymer composition of the copolymerisation of water soluble monomers MAA and PEGMAS.