Entropy Rules: Molecular Dynamics Simulations of Model Oligomers for Thermoresponsive Polymers

We attempted to attain atomic-scale insights into the mechanism of the heat-induced phase transition of two thermoresponsive polymers containing amide groups, poly(N-isopropylacrylamide) (PNIPAM) and poly(2-isopropyl-2-oxazoline) (PIPOZ), and we succeeded in reproducing the existence of lower critical solution temperature (LCST). The simulation data are in accord with experimental findings. We found out that the entropy has an important contribution to the thermodynamics of the phase separation transition. Moreover, after decomposing further the entropy change to contributions from the solutes and from the solvent, it appeared out that the entropy of the solvent has the decisive share for the lowering of the free energy of the system when increasing the temperature above the LCST. Our conclusion is that the thermoresponsive behavior is driven by the entropy of the solvent. The water molecules structured around the functional groups of the polymer that are exposed to contact with the solvent in the extended conformation lower the enthalpy of the system, but at certain temperature the extended conformation of the polymer collapses as a result of dominating entropy gain from “released” water molecules. We stress also on the importance of using more than one reference molecule in the simulation box at the setup of the simulation.

[1]  C. Atilgan,et al.  Multiscale modeling of poly(2-isopropyl-2-oxazoline) chains in aqueous solution , 2017 .

[2]  Berk Hess,et al.  GROMACS: High performance molecular simulations through multi-level parallelism from laptops to supercomputers , 2015 .

[3]  H. G. Schild.  Poly(N-isopropylacrylamide): experiment, theory and application , 1992 .

[4]  Toshiyuki Shikata,et al.  Contrary hydration behavior of N-isopropylacrylamide to its polymer, P(NIPAm), with a lower critical solution temperature. , 2007, The journal of physical chemistry. B.

[5]  Vincent D Ustach,et al.  Direct Phase Equilibrium Simulations of NIPAM Oligomers in Water. , 2016, The journal of physical chemistry. B.

[6]  Patrick K. Quoika,et al.  Coil-Globule Transition Thermodynamics of Poly(N-isopropylacrylamide). , 2019, The journal of physical chemistry. B.

[7]  P. Kollman,et al.  Automatic atom type and bond type perception in molecular mechanical calculations. , 2006, Journal of molecular graphics & modelling.

[8]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[9]  J. Cerdà,et al.  Particle-particle particle-mesh method for dipolar interactions: on error estimates and efficiency of schemes with analytical differentiation and mesh interlacing. , 2011, The Journal of chemical physics.

[10]  Berk Hess,et al.  LINCS: A linear constraint solver for molecular simulations , 1997 .

[11]  M. Stevens,et al.  Study of the Polymer Length Dependence of the Single Chain Transition Temperature in Syndiotactic Poly(N-isopropylacrylamide) Oligomers in Water , 2012 .

[12]  Hoover,et al.  Canonical dynamics: Equilibrium phase-space distributions. , 1985, Physical review. A, General physics.

[13]  Afsaneh Lavasanifar,et al.  Amphiphilic block copolymers for drug delivery. , 2003, Journal of pharmaceutical sciences.

[14]  J. Schlitter.  Estimation of absolute and relative entropies of macromolecules using the covariance matrix , 1993 .

[15]  G. Odriozola,et al.  Conformation change of an isotactic poly (N-isopropylacrylamide) membrane: Molecular dynamics. , 2017, The Journal of chemical physics.

[16]  Sanket A. Deshmukh,et al.  Role of solvation dynamics and local ordering of water in inducing conformational transitions in poly(N-isopropylacrylamide) oligomers through the LCST. , 2012, The journal of physical chemistry. B.

[17]  E. Spohr,et al.  Molecular Dynamics Investigation of the Thermo‐Responsive Polymer Poly(N‐isopropylacrylamide) , 2012 .

[18]  M. Karplus,et al.  Method for estimating the configurational entropy of macromolecules , 1981 .

[19]  X. Qiu,et al.  Dissecting the Mechanism of the Heat-Induced Phase Separation and Crystallization of Poly(2-isopropyl-2-oxazoline) in Water through Vibrational Spectroscopy and Molecular Orbital Calculations , 2012 .

[20]  D. van der Spoel,et al.  GROMACS: A message-passing parallel molecular dynamics implementation , 1995 .

[21]  G. Assche,et al.  Demixing and remixing kinetics of poly(2-isopropyl-2-oxazoline) (PIPOZ) aqueous solutions studied by modulated temperature differential scanning calorimetry , 2010 .

[22]  Alexander D. MacKerell,et al.  All-atom empirical potential for molecular modeling and dynamics studies of proteins. , 1998, The journal of physical chemistry. B.

[23]  C. Marques,et al.  Molecular dynamics study of the LCST transition in aqueous poly(N-n-propylacrylamide). , 2018, Physical chemistry chemical physics : PCCP.

[24]  Martin Kröger,et al.  Poly(N-isopropylacrylamide) Phase Diagrams: Fifty Years of Research. , 2015, Angewandte Chemie.

[25]  Jun Soo Kim,et al.  Collapse-Swelling Transitions of a Thermoresponsive, Single Poly(N-isopropylacrylamide) Chain in Water. , 2016, The journal of physical chemistry. B.

[26]  G. Longhi,et al.  Molecular dynamics simulation of a model oligomer for poly(N-isopropylamide) in water , 2004 .

[27]  S. Nosé.  A unified formulation of the constant temperature molecular dynamics methods , 1984 .

[28]  Shengtong Sun,et al.  Conformational changes in the heat-induced crystallization of poly(2-isopropyl-2-oxazoline) in the solid state. , 2015, Physical chemistry chemical physics : PCCP.

[29]  M. Stevens,et al.  Single chain structure of a poly(N-isopropylacrylamide) surfactant in water. , 2015, The journal of physical chemistry. B.

[30]  P. Kollman,et al.  A well-behaved electrostatic potential-based method using charge restraints for deriving atomic char , 1993 .

[31]  M. Textor,et al.  Polyoxazolines for nonfouling surface coatings--a direct comparison to the gold standard PEG. , 2012, Macromolecular rapid communications.

[32]  Berk Hess,et al.  GROMACS 3.0: a package for molecular simulation and trajectory analysis , 2001 .

[33]  X. Daura,et al.  Entropy calculations on a reversibly folding peptide: Changes in solute free energy cannot explain folding behavior , 2001, Proteins.

[34]  Axel H. E. Müller,et al.  Thermosensitive water-soluble copolymers with doubly responsive reversibly interacting entities , 2007 .

[35]  Wim F Vranken,et al.  ACPYPE - AnteChamber PYthon Parser interfacE , 2012, BMC Research Notes.

[36]  A. Lonardi,et al.  Toward modeling thermoresponsive polymer networks: a molecular dynamics simulation study of N-isopropyl acrylamide co-oligomers. , 2010, Journal of Physical Chemistry B.

[37]  Carsten Kutzner,et al.  GROMACS 4:  Algorithms for Highly Efficient, Load-Balanced, and Scalable Molecular Simulation. , 2008, Journal of chemical theory and computation.

[38]  N. V. D. van der Vegt,et al.  Improved Temperature Behavior of PNIPAM in Water with a Modified OPLS Model. , 2019, The journal of physical chemistry. B.

[39]  J. Andrew McCammon,et al.  Absolute Single-Molecule Entropies from Quasi-Harmonic Analysis of Microsecond Molecular Dynamics: Correction Terms and Convergence Properties , 2009, Journal of chemical theory and computation.

[40]  Gerrit Groenhof,et al.  GROMACS: Fast, flexible, and free , 2005, J. Comput. Chem..

[41]  K Schulten,et al.  VMD: visual molecular dynamics. , 1996, Journal of molecular graphics.

[42]  Young-Jin Kim,et al.  Thermo-responsive polymers and their application as smart biomaterials. , 2017, Journal of materials chemistry. B.

[43]  D. Moscatelli,et al.  Thermo-responsive polymers: Applications of smart materials in drug delivery and tissue engineering. , 2019, Materials science & engineering. C, Materials for biological applications.

[44]  E. Chiessi,et al.  Influence of Tacticity on Hydrophobicity of Poly(N-isopropylacrylamide): A Single Chain Molecular Dynamics Simulation Study. , 2016, The journal of physical chemistry. B.

[45]  R. Cordone,et al.  Molecular dynamics simulation of aqueous solutions of 26-unit segments of p(NIPAAm) and of p(NIPAAm) "doped" with amino acid based comonomers. , 2008, The journal of physical chemistry. B.