1
|
Yagasaki T, Matubayasi N. Molecular Dynamics Simulations of Concentrated and Dilute Aqueous Solutions of Poly( N-Isopropylacrylamide) Using a Modified OPLS-AA Model. J Phys Chem B 2025. [PMID: 40372464 DOI: 10.1021/acs.jpcb.5c00789] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/16/2025]
Abstract
We perform molecular dynamics (MD) simulations of aqueous poly(N-isopropylacrylamide) (PNIPAM) solutions using the OPLS-AA model and its modified version in combination with the TIP4P/2005 water model. The original OPLS-AA model underestimates the demixing temperature of concentrated solutions and the coil-globule transition temperature in dilute solutions, whereas the modified model reproduces both. It was reported that the transition from the globule state to the coil state never occurred in MD simulations with the original OPLS-AA model. This problem is avoided by using the modified model because of the fast dynamics near the coil-globule transition temperature, which is higher for the modified model than for the original model by ∼60 K. The distribution functions of two structural order parameters calculated from well equilibrated trajectories of the modified model clearly show that the coil-globule transition is not a discontinuous first-order transition, as suggested in previous MD studies, but a continuous transition as predicted from mean field theory. We also find that the effect of polymer-polymer hydrogen bonding on the dynamics of aggregated polymer chains is less significant than assumed in experimental studies.
Collapse
Affiliation(s)
- Takuma Yagasaki
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka 560-8531, Japan
| | - Nobuyuki Matubayasi
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka 560-8531, Japan
| |
Collapse
|
2
|
Mandalaparthy V, van der Vegt NFA. A generic model for pH-sensitive collapse of hydrophobic polymers. Phys Chem Chem Phys 2025; 27:6984-6993. [PMID: 40104906 DOI: 10.1039/d4cp04756g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/20/2025]
Abstract
The hydrophobic effect is an important contributor to the stability of proteins and may be influenced by many factors including the pH of the solution. To simplify the study of pH effects on proteins, we parameterize biologically motivated titratable monomers which we insert into the sequence of a hydrophobic polymer and study via constant pH molecular dynamics (MD) simulations. We calculate the potential of mean force of the polymer as a function of its radius of gyration at different pH values and observe that the collapsed state of the polymer is destabilized when the titratable monomer is more charged (high pH for an acid and low pH for a base). Further, the extent of the destabilization is influenced by the position of the titratable monomer along the polymer sequence. The pKa value of the titratable monomer is also observed to be sensitive to polymer conformation, in agreement with protein studies. We further study a zwitterionic polymer with an acidic and a basic monomer in the same sequence which presents a pH-dependent hairpin formation. Our model provides a simplified yet powerful framework to study pH effects on the hydrophobic effect, providing insights into mechanisms governing the behavior of intrinsically disordered proteins (IDPs) and pH-sensitive drug delivery, among other applications.
Collapse
Affiliation(s)
- Varun Mandalaparthy
- Department of Chemistry, Technical University of Darmstadt, 64287 Darmstadt, Germany.
| | - Nico F A van der Vegt
- Department of Chemistry, Technical University of Darmstadt, 64287 Darmstadt, Germany
| |
Collapse
|
3
|
Mortara L, Mukhina T, Chaimovich H, Brezesinski G, van der Vegt NFA, Schneck E. Anion Competition at Positively Charged Surfactant Monolayers. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2024; 40:6949-6961. [PMID: 38502024 DOI: 10.1021/acs.langmuir.3c04003] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 03/20/2024]
Abstract
Interactions of anions with hydrophobic surfaces of proteins and water-soluble polymers depend on the ability of the ions to shed their hydration shells. At positively charged surfactant monolayers, the interactions of anions are less well understood. Due to the interplay of electrostatic surface forces, hydration effects, and ion-ion interactions in the electrostatic double layer, a comprehensive microscopic picture remains elusive. Herein, we study the interactions of chloride, bromide, and a mixture of these two anions at the aqueous interface of dihexadecyldimethylammonium (DHDA+) and dioctadecyldimethylammonium (DODA+) cationic monolayers. Using molecular dynamics simulations and three surface-sensitive X-ray scattering techniques, we demonstrate that bromide interacts preferentially over chloride with both monolayers. The structure of the two monolayers and their interfacial electron density profiles obtained from the simulations quantitatively reproduce the experimental data. We observe that chloride and bromide form contact ion pairs with the quaternary ammonium groups on both monolayers. However, ion pairing with bromide leads to a greater reduction in the number of water molecules hydrating the anion, resulting in more energetically stable ion pairs. This leads to long-range (>3 nm) lateral correlations between bromide ions on the structured DODA+ monolayer. These observations indicate that ion hydration is the dominant factor determining the interfacial electrolyte structure.
Collapse
Affiliation(s)
- Laura Mortara
- Chemistry Institute, University of São Paulo, São Paulo, SP 05508-000, Brazil
- Physics Department, Technical University of Darmstadt, Darmstadt 64289, Germany
| | - Tetiana Mukhina
- Physics Department, Technical University of Darmstadt, Darmstadt 64289, Germany
| | - Hernan Chaimovich
- Chemistry Institute, University of São Paulo, São Paulo, SP 05508-000, Brazil
| | - Gerald Brezesinski
- Physics Department, Technical University of Darmstadt, Darmstadt 64289, Germany
| | | | - Emanuel Schneck
- Physics Department, Technical University of Darmstadt, Darmstadt 64289, Germany
| |
Collapse
|
4
|
Higuchi H, Ikeda-Fukazawa T. Interactions between Water and a Hydrophobic Polymer. J Phys Chem B 2024; 128:1927-1935. [PMID: 38369787 DOI: 10.1021/acs.jpcb.3c07440] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/20/2024]
Abstract
To investigate the mechanisms of interactions between a hydrophobic polymer and water, molecular dynamics calculations and Raman spectroscopic measurements of cis-1,4-polyisoprene (PI)-water systems were performed. The results show that PI in water undergoes a coil-globule transition at around 248 K. The transition is attributed to changes in the density and diffusivity of water. The volume expansion of the supercooled liquid water induces the coil structure of PI. The phase separation of PI from water with an increase in the self-diffusion coefficient of water molecules results in the globule structure of PI. The self-diffusion coefficient of free water with PI is larger than that of pure water because PI has an effect to decrease the hydrogen-bonding strength of water. The result suggests that the effects of the coexisting water are important factors governing the physical and chemical properties of hydrophobic polymers.
Collapse
Affiliation(s)
- Hikaru Higuchi
- Department of Applied Chemistry, Meiji University, Kawasaki 214-8571, Japan
| | | |
Collapse
|
5
|
Zhao Y, Bharadwaj S, Myers RL, Okur HI, Bui PT, Cao M, Welsh LK, Yang T, Cremer PS, van der Vegt NFA. Solvation Behavior of Elastin-like Polypeptides in Divalent Metal Salt Solutions. J Phys Chem Lett 2023; 14:10113-10118. [PMID: 37921693 DOI: 10.1021/acs.jpclett.3c02476] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/04/2023]
Abstract
The effects of CaCl2 and MgCl2 on the cloud point temperature of two different elastin-like polypeptides (ELPs) were studied using a combination of cloud point measurements, molecular dynamics simulations, and infrared spectroscopy. Changes in the cloud point for the ELPs in aqueous divalent metal cation solutions were primarily governed by two competing interactions: the cation-amide oxygen electrostatic interaction and the hydration of the cation. In particular, Ca2+ cations can more readily shed their hydration shells and directly contact two amide oxygens by the formation of ion bridges. By contrast, Mg2+ cations were more strongly hydrated and preferred to partition toward the amide oxygens along with their hydration shells. In fact, although hydrophilic ELP V5A2G3 was salted-out at low concentrations of MgCl2, it was salted-in at higher salt concentrations. By contrast, CaCl2 salted the ELP sharply out of solution at higher salt concentrations because of the bridging effect.
Collapse
Affiliation(s)
- Yani Zhao
- Department of Chemistry, Technical University of Darmstadt, 64287 Darmstadt, Germany
| | - Swaminath Bharadwaj
- Department of Chemistry, Technical University of Darmstadt, 64287 Darmstadt, Germany
- Department of Chemical Engineering, Shiv Nadar Institution of Eminence, Gautam Buddha Nagar, Uttar Pradesh 201314, India
| | - Ryan L Myers
- Department of Chemistry, Penn State University, University Park, Pennsylvania 16802, United States
| | - Halil I Okur
- Department of Chemistry, Penn State University, University Park, Pennsylvania 16802, United States
| | - Pho T Bui
- Department of Chemistry, Penn State University, University Park, Pennsylvania 16802, United States
| | - Mengrui Cao
- Department of Chemistry, Penn State University, University Park, Pennsylvania 16802, United States
| | - Lauren K Welsh
- Department of Chemistry, Penn State University, University Park, Pennsylvania 16802, United States
| | - Tinglu Yang
- Department of Chemistry, Penn State University, University Park, Pennsylvania 16802, United States
| | - Paul S Cremer
- Department of Chemistry, Penn State University, University Park, Pennsylvania 16802, United States
- Department of Biochemistry and Molecular Biology, Penn State University, University Park, Pennsylvania 16802, United States
| | - Nico F A van der Vegt
- Department of Chemistry, Technical University of Darmstadt, 64287 Darmstadt, Germany
| |
Collapse
|
6
|
Seki T, Yu CC, Chiang KY, Greco A, Yu X, Matsumura F, Bonn M, Nagata Y. Ions Speciation at the Water-Air Interface. J Am Chem Soc 2023; 145:10622-10630. [PMID: 37139910 DOI: 10.1021/jacs.3c00517] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/05/2023]
Abstract
In typical aqueous systems, including naturally occurring sweet and salt water and tap water, multiple ion species are co-solvated. At the water-air interface, these ions are known to affect the chemical reactivity, aerosol formation, climate, and water odor. Yet, the composition of ions at the water interface has remained enigmatic. Here, using surface-specific heterodyne-detected sum-frequency generation spectroscopy, we quantify the relative surface activity of two co-solvated ions in solution. We find that more hydrophobic ions are speciated to the interface due to the hydrophilic ions. Quantitative analysis shows that the interfacial hydrophobic ion population increases with decreasing interfacial hydrophilic ion population at the interface. Simulations show that the solvation energy difference between the ions and the intrinsic surface propensity of ions determine the extent of an ion's speciation by other ions. This mechanism provides a unified view of the speciation of monatomic and polyatomic ions at electrolyte solution interfaces.
Collapse
Affiliation(s)
- Takakazu Seki
- Max Planck Institute for Polymer Research, Ackermannweg 10, Mainz 55128, Germany
- Graduate School of Science and Technology, Hirosaki University, Hirosaki 036-8561, Aomori, Japan
| | - Chun-Chieh Yu
- Max Planck Institute for Polymer Research, Ackermannweg 10, Mainz 55128, Germany
| | - Kuo-Yang Chiang
- Max Planck Institute for Polymer Research, Ackermannweg 10, Mainz 55128, Germany
| | - Alessandro Greco
- Max Planck Institute for Polymer Research, Ackermannweg 10, Mainz 55128, Germany
| | - Xiaoqing Yu
- Max Planck Institute for Polymer Research, Ackermannweg 10, Mainz 55128, Germany
| | - Fumiki Matsumura
- Max Planck Institute for Polymer Research, Ackermannweg 10, Mainz 55128, Germany
| | - Mischa Bonn
- Max Planck Institute for Polymer Research, Ackermannweg 10, Mainz 55128, Germany
| | - Yuki Nagata
- Max Planck Institute for Polymer Research, Ackermannweg 10, Mainz 55128, Germany
| |
Collapse
|
7
|
Lu H, Wang D, Huang D, Feng L, Zhang H, Zhu P. Product from sessile droplet evaporation of PNIPAM/water system above LCST: A block or micro/nano-particles? J Colloid Interface Sci 2023; 634:769-781. [PMID: 36565619 DOI: 10.1016/j.jcis.2022.12.097] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/07/2022] [Revised: 12/16/2022] [Accepted: 12/18/2022] [Indexed: 12/24/2022]
Abstract
PNIPAM as a stimuli-responsive polymer has generated extreme interests due to its versatile applications. However, there is no research report on whether PNIPAM micro/nano-particles can be extracted from its suspension after phase separation. In the present work, LCST-type phase separation in self-synthesized PNIPAM/water system was investigated in depth by dividing the DLS testing process into four stages. In addition to quenching duration, temperature rise process, quenching temperature and PNIPAM concentration all have a great influence on particle size of the suspension. Meanwhile, the steady-state rheology and dynamic viscoelasticity results show that PNIPAM micro/nano-particles in the suspension are "soft" that can deform. Finally, FE-SEM was used to observe the morphology of dehydrated PNIPAM extracted by sessile droplet evaporation under different conditions. The results indicate that these "soft" particles are easier to fuse together, do not have sufficient mechanical strength to maintain their spherical morphology after dehydration. But the above fusion can be suppressed by adjusting evaporation conditions to acquire smaller PNIPAM particles which have sufficient mechanical properties to keep their basic particle morphology. Further, by changing evaporation pressure to positive or negative pressure, dehydrated PNIPAM micro/nano-particles with excellent uniformity and separation can be obtained. This work will provide guidance for extracting micro/nano-particles from polymer/diluent systems with LCST.
Collapse
Affiliation(s)
- Hongwei Lu
- College of Textile Science and Engineering (International Institute of Silk), Zhejiang Sci-Tech University, Hangzhou, Zhejiang 310018, China
| | - Danling Wang
- Zhongce Rubber Group Co., Ltd, Hangzhou, Zhejiang 310018, China
| | - Daye Huang
- Zhongce Rubber Group Co., Ltd, Hangzhou, Zhejiang 310018, China
| | - Luyao Feng
- College of Textile Science and Engineering (International Institute of Silk), Zhejiang Sci-Tech University, Hangzhou, Zhejiang 310018, China
| | - Huapeng Zhang
- College of Textile Science and Engineering (International Institute of Silk), Zhejiang Sci-Tech University, Hangzhou, Zhejiang 310018, China
| | - Peng Zhu
- College of Textile Science and Engineering (International Institute of Silk), Zhejiang Sci-Tech University, Hangzhou, Zhejiang 310018, China.
| |
Collapse
|
8
|
Mandalaparthy V, Noid WG. A simple theory for interfacial properties of dilute solutions. J Chem Phys 2022; 157:034703. [DOI: 10.1063/5.0098579] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Recent studies suggest that cosolute mixtures may exert significant non-additive effects upon protein stability. The corresponding liquid–vapor interfaces may provide useful insight into these non-additive effects. Accordingly, in this work, we relate the interfacial properties of dilute multicomponent solutions to the interactions between solutes. We first derive a simple model for the surface excess of solutes in terms of thermodynamic observables. We then develop a lattice-based statistical mechanical perturbation theory to derive these observables from microscopic interactions. Rather than adopting a random mixing approximation, this dilute solution theory (DST) exactly treats solute–solute interactions to lowest order in perturbation theory. Although it cannot treat concentrated solutions, Monte Carlo (MC) simulations demonstrate that DST describes the interactions in dilute solutions with much greater accuracy than regular solution theory. Importantly, DST emphasizes a fundamental distinction between the “intrinsic” and “effective” preferences of solutes for interfaces. DST predicts that three classes of solutes can be distinguished by their intrinsic preference for interfaces. While the surface preference of strong depletants is relatively insensitive to interactions, the surface preference of strong surfactants can be modulated by interactions at the interface. Moreover, DST predicts that the surface preference of weak depletants and weak surfactants can be qualitatively inverted by interactions in the bulk. We also demonstrate that DST can be extended to treat surface polarization effects and to model experimental data. MC simulations validate the accuracy of DST predictions for lattice systems that correspond to molar concentrations.
Collapse
Affiliation(s)
- Varun Mandalaparthy
- Department of Chemistry, Penn State University, University Park, State College, Pennsylvania 16802, USA
| | - W. G. Noid
- Department of Chemistry, Penn State University, University Park, State College, Pennsylvania 16802, USA
| |
Collapse
|