The Journal of Physical Chemistry B · 2019 · 81 citations · 80 references
Natural SciencesChemistryPhysicsComputational ThermodynamicsPolymer ChemistryThermodynamicsPolymer Phase BehaviorPolymer Physical ChemistryPolymer PhysicsPolymer SimulationPolymer StructurePolymer ThermodynamicsCoiled-coil StructureCritical Solution TemperatureCritical TemperatureLower Critical Solution TemperaturePolymer Coil ConformationPolymer Coil-globule TransitionPolymer Solution ThermodynamicsThermodynamic Analysis
-isopropylacrylamide) (PNIPAM) undergo a phase transition in aqueous solution from a random-coil structural ensemble to a globule structural ensemble at the lower critical solution temperature (LCST). Above this temperature, PNIPAM agglomerates and becomes insoluble, whereas it is soluble below the temperature. Thus, thermosensitive polymers represent essential targets for several applications, e.g., in drug delivery. Although their ability to change structure in response to a temperature alteration is highly relevant for industrial processes, their thermodynamic properties are mostly qualitatively understood, and the quantitative thermodynamic picture is still elusive. In this study, we used a combined atomistic molecular dynamics and well-tempered metadynamics simulation approach to estimate coil-globule transition thermodynamics. An isotactic 30-mer of PNIPAM was investigated over a broad temperature range between 200 and 360 K. The transition from the globule to the random-coil structure was observed with well-tempered metadynamics. For the first time, the free energy surface of PNIPAM was estimated and it is shown that the simulation results are in line with the experimentally observed thermosensitive behavior. Below the LCST, the random-coil ensemble represents the global energy minimum and is thermodynamically favored by 21 ± 9 kJ/mol compared to the globule ensemble; both are separated by a barrier of 49 ± 14 kJ/mol. In contrast, above the LCST, the globule ensemble is thermodynamically favored by 21 ± 8 kJ/mol over the random-coil ensemble. The barrier from random-coil to globule is 17 ± 10 kJ/mol.
80
A smooth particle mesh Ewald method
Ulrich Essmann, L. Perera, Max L. Berkowitz et al. · The Journal of Chemical Physics · 1995 · 22.9K citations
Canonical sampling through velocity rescaling
Giovanni Bussi, Davide Donadio, Michele Parrinello · The Journal of Chemical Physics · 2007 · 18.4K citations · Full text
Natural Sciences, Chemistry, Physics +17
LINCS: A linear constraint solver for molecular simulations
Berk Hess, Henk Bekker, Herman J. C. Berendsen et al. · Journal of Computational Chemistry · 1997 · 17.1K citations
Natural Sciences, Chemistry, Physics +14